| 151. |
- Kaukonen, Markus, et al.
(författare)
-
Proton transfer at metal sites in proteins studied by quantum mechanical free-energy perturbations
- 2008
-
Ingår i: Journal of Chemical Theory and Computation. - : American Chemical Society (ACS). - 1549-9618 .- 1549-9626. ; 4:6, s. 985-1001
-
Tidskriftsartikel (refereegranskat)abstract
- Catalytic metal sites in enzymes frequently have second-sphere carboxylate groups that neutralize the charge of the site and share protons with first-sphere ligands. This gives rise to an ambiguity concerning the position of this proton, which has turned out to be hard to settle with experimental, as well as theoretical, methods. We study three such proton-transfer reactions in two proteins and show that, in [Ni,Fe] hydrogenase, the bridging Cys-546 ligand is deprotonated by His-79, whereas in oxidized copper nitrite reductase, the His-100 ligand is neutral and the copper-bound water molecule is deprotonated by Asp-98. We show that these reactions strongly depend on the electrostatic interactions with the surrounding protein and solvent, because there is a large change in the dipole moment of the active site (2-6 D). Neither vacuum quantum mechanical (QM) calculations with large models, a continuum solvent, or a Poisson-Boltzmann treatment of the surroundings, nor combined QM and molecular mechanics (QM/MM) optimizations give reliable estimates of the proton-transfer energies (mean absolute deviations of over 20 kJ/mol). Instead, QM/MM free-energy perturbations are needed to obtain reliable estimates of the reaction energies. These calculations also indicate what interactions and residues are important for the energy, showing how the quantum system may be systematically enlarged. With such a procedure, results with an uncertainty of similar to 10 kJ/mol can be obtained, provided that a proper QM method is used.
|
|
| 152. |
- Kaukonen, Markus, et al.
(författare)
-
QM/MM-PBSA method to estimate free energies for reactions in proteins
- 2008
-
Ingår i: The Journal of Physical Chemistry Part B. - : American Chemical Society (ACS). - 1520-5207 .- 1520-6106. ; 112:39, s. 12537-12548
-
Tidskriftsartikel (refereegranskat)abstract
- We have developed a method to estimate free energies of reactions in proteins, called QM/MM-PBSA. It estimates the internal energy of the reactive site by quantum mechanical (QM) calculations, whereas bonded, electrostatic, and van der Waals interactions with the surrounding protein are calculated at the molecular mechanics (MM) level. The electrostatic part of the solvation energy of the reactant and the product is estimated by solving the Poisson-Boltzmann (PB) equation, and the nonpolar part of the solvation energy is estimated from the change in solvent-accessible surface area (SA). Finally, the change in entropy is estimated from the vibrational frequencies. We test this method for five proton-transfer reactions in the active sites of [Ni,Fe] hydrogenase and copper nitrite. reductase. We show that QM/MM-PBSA reproduces the results of a strict QM/MM free-energy perturbation method with a mean absolute deviation (MAD) of 8-10 kJ/mol if snapshots from molecular dynamics simulations are used and 4-14 kJ/mol if a single QM/MM structure is used. This is appreciably better than the original QM/MM results or if the QM energies are supplemented with a point-charge model, a self-consistent reaction field, or a PB model of the protein and the solvent, which give MADs of 22-36 kJ/mol for the same test set.
|
|
| 153. |
- Kelpšas, Vinardas, et al.
(författare)
-
Neutron structures of Leishmania mexicana triosephosphate isomerase in complex with reaction-intermediate mimics shed light on the proton-shuttling steps
- 2021
-
Ingår i: IUCrJ. - 2052-2525. ; 8:Pt 4, s. 633-643
-
Tidskriftsartikel (refereegranskat)abstract
- Triosephosphate isomerase (TIM) is a key enzyme in glycolysis that catalyses the interconversion of glyceraldehyde 3-phosphate and dihydroxy-acetone phosphate. This simple reaction involves the shuttling of protons mediated by protolysable side chains. The catalytic power of TIM is thought to stem from its ability to facilitate the deprotonation of a carbon next to a carbonyl group to generate an enediolate intermediate. The enediolate intermediate is believed to be mimicked by the inhibitor 2-phosphoglycolate (PGA) and the subsequent enediol intermediate by phosphoglycolohydroxamate (PGH). Here, neutron structures of Leishmania mexicana TIM have been determined with both inhibitors, and joint neutron/X-ray refinement followed by quantum refinement has been performed. The structures show that in the PGA complex the postulated general base Glu167 is protonated, while in the PGH complex it remains deprotonated. The deuteron is clearly localized on Glu167 in the PGA-TIM structure, suggesting an asymmetric hydrogen bond instead of a low-barrier hydrogen bond. The full picture of the active-site protonation states allowed an investigation of the reaction mechanism using density-functional theory calculations.
|
|
| 154. |
- Kongsted, Jacob, et al.
(författare)
-
An improved method to predict the entropy term with the MM/PBSA approach.
- 2009
-
Ingår i: Journal of Computer-Aided Molecular Design. - : Springer Science and Business Media LLC. - 1573-4951 .- 0920-654X. ; 23, s. 63-71
-
Tidskriftsartikel (refereegranskat)abstract
- A method is suggested to calculate improved entropies within the MM/PBSA approach (molecular mechanics combined with Poisson-Boltzmann and surface area calculations) to estimate protein-ligand binding affinities. In the conventional approach, the protein is truncated outside ~8 A from the ligand. This system is freely minimised using a distance-dependent dielectric constant (to simulate the removed protein and solvent). However, this can lead to extensive changes in the molecular geometry, giving rise to a large standard deviation in this term. In our new approach, we introduce a buffer region ~4 A outside the truncated protein (including solvent molecules) and keep it fixed during the minimisation. Thereby, we reduce the standard deviation by a factor of 2-4, ensuring that the entropy term no longer limits the precision of the MM/PBSA predictions. The new method is tested for the binding of seven biotin analogues to avidin, eight amidinobenzyl-indole-carboxamide inhibitors to factor Xa, and two substrates to cytochrome P450 3A4 and 2C9. It is shown that it gives more stable results and often improved predictions of the relative binding affinities.
|
|
| 155. |
- Kongsted, Jacob, et al.
(författare)
-
How accurate are continuum solvation models for drug-like molecules?
- 2009
-
Ingår i: Journal of Computer-Aided Molecular Design. - : Springer Science and Business Media LLC. - 1573-4951 .- 0920-654X. ; 23:7, s. 395-409
-
Tidskriftsartikel (refereegranskat)abstract
- We have estimated the hydration free energy for 20 neutral drug-like molecules, as well as for three series of 6-11 inhibitors to avidin, factor Xa, and galectin-3 with four different continuum solvent approaches (the polarised continuum method the Langevin dipole method, the finite-difference solution of the Poisson equation, and the generalised Born method), and several variants of each, giving in total 24 different methods. All four types of methods have been thoroughly calibrated for a number of experimentally known small organic molecules with a mean absolute deviation (MAD) of 1-6 kJ/mol for neutral molecules and 4-30 kJ/mol for ions. However, for the drug-like molecules, the accuracy seems to be appreciably worse. The reason for this is that drug-like molecules are more polar than small organic molecules and that the uncertainty of the methods is proportional to the size of the solvation energy. Therefore, the accuracy of continuum solvation methods should be discussed in relative, rather than absolute, terms. In fact, the mean unsigned relative deviations of the best solvation methods, 0.09 for neutral and 0.05 for ionic molecules, correspond to 2-20 kJ/mol absolute error for the drug-like molecules in this investigation, or 2-3,000 in terms of binding constants. Fortunately, the accuracy of all methods can be improved if only relative energies within a series of inhibitors are considered, especially if all of them have the same net charge. Then, all except two methods give MADs of 2-5 kJ/mol (corresponding to an uncertainty of a factor of 2-7 in the binding constant). Interestingly, the generalised Born methods typically give better results than the Poison-Boltzmann methods.
|
|
| 156. |
- Kongsted, Jacob, et al.
(författare)
-
Prediction and rationalization of the pH dependence of the activity and stability of family 11 xylanases
- 2007
-
Ingår i: Biochemistry. - : American Chemical Society (ACS). - 0006-2960 .- 1520-4995. ; 46:47, s. 13581-13592
-
Tidskriftsartikel (refereegranskat)abstract
- This paper presents a study of the pH dependence of the activity and stability of a set of family 11 xylanases for which X-ray structures are available, using the PROPKA approach. The xylanases are traditionally divided into basic and acidic xylanases, depending on whether the catalytic acid is hydrogen bonded to an Asn or Asp residue. Using X-ray structures, the predicted pH values of optimal activity of the basic xylanases are in the range of 5.2-6.9, which is in reasonable agreement with the available experimental values of 5-6.5. In the case of acidic xylanases, there are only four X-ray structures available, and using these structures, the predicted pHs of optimal activity are in the range of 4.2-5.0, compared to an observed range of 2-4.6. The influence of dynamical fluctuations of the protein structure is investigated for Bacillus agaradhaerens and Aspergillus kawachii xylanase using molecular dynamics (MD) simulations to provide snapshots from which average values can be computed. This decreases the respective predicted pH optima from 6.2-6.7 and 4.8 to 5.3 +/- 0.3 and 4.0 +/- 0.2, respectively, which are in better agreement with the observed values of 5.6 and 2, respectively. The change is primarily due to structural fluctuations of an Arg residue near the catalytic nucleophile, which lowers its pK(a) value compared to using the X-ray structure. The MD simulations and some X-ray structures indicate that this Arg residue can form a hydrogen bond to the catalytic base, and it is hypothesized that this hydrogen bond is stabilized by an additional hydrogen bond to another Glu residue present only in acidic xylanases. Formation of such a hydrogen bond is predicted to lower the pH optimum of A. kawachii xylanase to 2.9 +/- 0.3, which is in reasonable agreement, with the observed value of 2. The predicted pH of optimal stability is in excellent agreement with the pH value at which the melting temperature (T-m) is greatest. Some correlation is observed between the pH-dependent free energy of unfolding and T-m, suggesting that the thermostability of the xylanases is partly due to a difference in residues with shifted pK(a) values. Thus, the thermostability of xylanases (and proteins in general) can perhaps be increased by mutations that introduce ionizable residues with pK(a) values significantly lower than standard values.
|
|
| 157. |
- Kvassman, Jan, et al.
(författare)
-
Mechanism of glyceraldehyde‐3‐phosphate transfer from aldolase to glyceraldehyde‐3‐phosphate dehydrogenase
- 1988
-
Ingår i: European Journal of Biochemistry. - : Wiley. - 0014-2956 .- 1432-1033. ; 172:2, s. 427-431
-
Tidskriftsartikel (refereegranskat)abstract
- The catalytic interaction of glyceraldehyde‐3‐phosphate dehydrogenase with glyceraldehydes‐3‐phosphate has been examined by transient‐state kinetic methods. The results confirm previous reports that the apparent Km for oxidative phosphorylation of glyceraldehydes‐3‐phosphate decreases at least 50‐fold when the substrate is generated in a coupled reaction system through the action of aldolase on fructose 1,6‐bisphosphate, but lend no support to the proposal that glyceraldehydes 3‐phosphate is directly transferred between the two enzymes without prior release to the reaction medium. A theoretical analysis is presented which shows that the kinetic behaviour of the coupled two‐enzyme system is compatible in all respects tested with a free‐diffusion mechanism for the transfer of glyceraldehydes 3‐phosphate from the producing enzyme to the consuming one.
|
|
| 158. |
- Källrot, Niklas, et al.
(författare)
-
Theoretical study of structure of catalytic copper site in nitrite reductase
- 2005
-
Ingår i: International Journal of Quantum Chemistry. - : Wiley. - 0020-7608. ; 102:5, s. 520-541
-
Tidskriftsartikel (refereegranskat)abstract
- The catalytic copper site in nitrite reductase contains a Cu2+ ion bound to three histidine (His) ligands and a solvent molecule. Sites from various sources show a conspicuous variation in the structure. In some proteins, it is close to tetrahedral (even more so than are blue copper proteins), whereas in other proteins, it has a structure more similar to that expected for a type 2 copper site. We have studied this site with a number of theoretical methods, ranging from vacuum optimizations, combined quantum and molecular mechanics (QM/MM) optimization, quantum refinement (X-ray crystallography supplemented by quantum chemical calculations), and accurate energy calculations. We show that the difference in the structure arises from a movement of the solvent molecule and that this movement is determined by a compromise between its hydrogen bond interactions and the intrinsic preferences of the copper site. If the solvent molecule is deprotonated, the two structures have a similar energy, whereas if it is protonated, the more tetrahedral structure is energetically favorable. Neither of the structures involves a pi interaction as in the blue copper proteins; instead, both are strongly distorted tetragonal structures with sigma bonds to all four ligands. We have also examined the position of hydrogen atoms shared between second-sphere carboxylate groups and the first-sphere solvent molecule and one of the His ligands. In the oxidized state, the structure with the solvent deprotonate(d) but the His residue protonated seems to be most stable. (c) 2004 Wiley Periodicals, Inc.
|
|
| 159. |
- Larsson, Ernst D., et al.
(författare)
-
Is density functional theory accurate for lytic polysaccharide monooxygenase enzymes
- 2020
-
Ingår i: Dalton Transactions. - : Royal Society of Chemistry (RSC). - 1477-9226 .- 1477-9234. ; 49:5, s. 1501-1512
-
Tidskriftsartikel (refereegranskat)abstract
- The lytic polysaccharide monooxygenase (LPMO) enzymes boost polysaccharide depolymerization through oxidative chemistry, which has fueled the hope for more energy-efficient production of biofuel. We have recently proposed a mechanism for the oxidation of the polysaccharide substrate (E. D. Hedegård and U. Ryde, Chem. Sci., 2018, 9, 3866-3880). In this mechanism, intermediates with superoxide, oxyl, as well as hydroxyl (i.e. [CuO2]+, [CuO]+ and [CuOH]2+) cores were involved. These complexes can have both singlet and triplet spin states, and both spin-states may be important for how LPMOs function during catalytic turnover. Previous calculations on LPMOs have exclusively been based on density functional theory (DFT). However, different DFT functionals are known to display large differences for spin-state splittings in transition-metal complexes, and this has also been an issue for LPMOs. In this paper, we study the accuracy of DFT for spin-state splittings in superoxide, oxyl, and hydroxyl intermediates involved in LPMO turnover. As reference we employ multiconfigurational perturbation theory (CASPT2).
|
|
| 160. |
- Lewin, Marcus, et al.
(författare)
-
Hybrids between Rubus caesius and Rubus sect. Corylifolii (Rosaceae) and their relation to R. cyclomorphus, R. tiliaster, R. glauciformis, R. slesvicensis and R. firmus
- 2022
-
Ingår i: Nordic Journal of Botany. - : Wiley. - 0107-055X .- 1756-1051. ; 2022:12
-
Tidskriftsartikel (refereegranskat)abstract
- Using flow cytometry and microsatellite DNA analysis (MSDA), we studied hybrids between Rubus caesius and various other Rubus species, with an emphasis on Sweden. We show that hybrids between Rubus caesius and Rubus sect. Corylifolii arise easily. They show a large variation in morphology, but can normally be recognised by a number of characters. They are typically hexaploids, but ~13% of the hybrids are pentaploids and ~8% tetraploids. With MSDA, they are harder to identify than hybrids with R. idaeus, partly because all Corylifolii species have themselves arisen from R. caesius hybrids and therefore share many alleles with R. caesius, partly because hybridisation with R. caesius seems to give rise to some variation in microsatellite regions. Hybrids with species of R. sect. Rubus are much rarer in Sweden and we have only identified one such case. We show that R. cyclomorphus and R. tiliaster are not proper apomictic species, at least not in Sweden, but rather a collection of genetically unrelated hybrids between R. caesius and R. raduloides or R. camptostachys. MSDA analysis of R. slesvicensis and R. firmus shows that these species are problematic. We identify two distinct taxa, one from Skåne in Sweden, which we describe as a new species, R. ruboculus and another from the neighbourhood of Viborg and Schleswig, although the latter is probably not identical with R. slesvicensis s.s. All the other samples form a diverse group of putative R. caesius hybrids at three different ploidy levels. Rubus glauciformis is uniform in Småland, Öland and Blekinge, but becomes harder to distinguish from various hybrids in Skåne. We conclude that once R. caesius and R. idaeus hybrids are properly identified, along with a small number of new apomictic species (often with a local distribution), the genus Rubus does not pose any serious taxonomical problems in Sweden.
|
|
