| 41. |
- Cao, Lili, et al.
(författare)
-
Quantum refinement with multiple conformations : Application to the P-cluster in nitrogenase
- 2020
-
Ingår i: Acta Crystallographica Section D: Structural Biology. - 2059-7983. ; 76, s. 1145-1156
-
Tidskriftsartikel (refereegranskat)abstract
- X-ray crystallography is the main source of atomistic information on the structure of proteins. Normal crystal structures are obtained as a compromise between the X-ray scattering data and a set of empirical restraints that ensure chemically reasonable bond lengths and angles. However, such restraints are not always available or accurate for nonstandard parts of the structure, for example substrates, inhibitors and metal sites. The method of quantum refinement, in which these empirical restraints are replaced by quantum-mechanical (QM) calculations, has previously been suggested for small but interesting parts of the protein. Here, this approach is extended to allow for multiple conformations in the QM region by performing separate QM calculations for each conformation. This approach is shown to work properly and leads to improved structures in terms of electron-density maps and real-space difference density Z-scores. It is also shown that the quality of the structures can be gauged using QM strain energies. The approach, called ComQumX-2QM, is applied to the P-cluster in two different crystal structures of the enzyme nitrogenase, i.e. an Fe 8 S 7 Cys 6 cluster, used for electron transfer. One structure is at a very high resolution (1.0 Å) and shows a mixture of two different oxidation states, the fully reduced P N state (Fe 8 2+, 20%) and the doubly oxidized P 2+ state (80%). In the original crystal structure the coordinates differed for only two iron ions, but here it is shown that the two states also show differences in other atoms of up to 0.7 Å. The second structure is at a more modest resolution, 2.1 Å, and was originally suggested to show only the one-electron oxidized state, P 1+. Here, it is shown that it is rather a 50/50% mixture of the P 1+ and P 2+ states and that many of the Fe - Fe and Fe - S distances in the original structure were quite inaccurate (by up to 0.8 Å). This shows that the new ComQumX-2QM approach can be used to sort out what is actually seen in crystal structures with dual conformations and to give locally improved coordinates.
|
|
| 42. |
- Cao, Lili, et al.
(författare)
-
What Is the Structure of the E4 Intermediate in Nitrogenase?
- 2020
-
Ingår i: Journal of Chemical Theory and Computation. - : American Chemical Society (ACS). - 1549-9618 .- 1549-9626. ; 16, s. 1936-1952
-
Tidskriftsartikel (refereegranskat)abstract
- Nitrogenase is the only enzyme that can cleave the strong triple bond inN2. The active site contains a complicated MoFe7S9C cluster. It is believed that itneeds to accept four protons and electrons, forming the E4 state, before it can bind N2.However, there is no consensus on the atomic structure of the E4 state. Experimentalstudies indicate that it should contain two hydride ions bridging two pairs of Fe ions,and it has been suggested that both hydride ions as well as the two protons bind on thesame face of the cluster. On the other hand, density functional theory (DFT) studieshave indicated that it is energetically more favorable with either three hydride ions or with a triply protonated carbide ion, depending on the DFT functional. We have performed a systematic combined quantum mechanical and molecular mechanical (QM/MM) study of possible E4 states with two bridging hydride ions. Our calculations suggest that the most favorable structure has hydride ions bridging the Fe2/6 and Fe3/7 ion pairs. In fact, such structures are 14 kJ/mol more stable than structures with three hydride ions, showing that pure DFT functionals give energetically most favorable structures in agreement with experiments. An important reason for this finding is that we have identified a new type of broken-symmetry state that involves only two Fe ions with minority spin, in contrast to the previously studied states with three Fe ions with minority spin. The energetically best structures have the two hydride ions on different faces of the FeMo cluster, whereas better agreement with ENDOR data is obtained if they are on the same face; such structures are only 6−22 kJ/mol less stable.
|
|
| 43. |
- Chalupsky, Jakub, et al.
(författare)
-
Multireference ab initio calculations on reaction intermediates of the multicopper oxidases
- 2006
-
Ingår i: Inorganic Chemistry. - : American Chemical Society (ACS). - 1520-510X .- 0020-1669. ; 45:26, s. 11051-11059
-
Tidskriftsartikel (refereegranskat)abstract
- The multicopper oxidases (MCOs) couple the four-electron reduction of dioxygen to water with four one-electron oxidations of various substrates. Extensive spectroscopic studies have identified several intermediates in the MCO catalytic cycle, but they have not been able to settle the structures of three of the intermediates, viz. the native intermediate (NI), the peroxy intermediate (PI), and the peroxy adduct (PA). The suggested structures have been further refined and characterized by quantum mechanical/molecular mechanical (QM/MM) calculations. In this paper, we try to establish a direct link between theory and experiment, by calculating spectroscopic parameters for these intermediates using multireference wave functions from the multistate CASPT2 and MRDDCI2 methods. Thereby, we have been able to reproduce low-spin ground states (S = 0 or S = 1/2) for all the MCO intermediates, as well as a low-lying (similar to 150 cm(-1)) doublet state and a doublet-quartet energy gap of similar to 780 cm(-1) for the NI. Moreover, we reproduce the zero-field splitting (similar to 70 cm(-1)) of the ground E-2 state in a D-3 symmetric hydroxy-bridged trinuclear Cu(II) model of the NI and obtain a quantitatively correct quartet-doublet splitting (164 cm(-1)) for a mu 3-oxo-bridged trinuclear Cu( II) cluster. All results support the suggestion that the NI has an O-2-atom in the center of the trinuclear cluster, whereas both the PI and PA have an O-2(2-) ion in the center of the cluster, in agreement with the QM/MM results and spectroscopic measurements.
|
|
| 44. |
- Ciancetta, Antonella, et al.
(författare)
-
A QM/MM study of the binding of RAPTA ligands to cathepsin B
- 2011
-
Ingår i: Journal of Computer-Aided Molecular Design. - : Springer Science and Business Media LLC. - 1573-4951 .- 0920-654X. ; 25:8, s. 729-742
-
Tidskriftsartikel (refereegranskat)abstract
- We have carried out quantum mechanical (QM) and QM/MM (combined QM and molecular mechanics) calculations, as well as molecular dynamics (MD) simulations to study the binding of a series of six RAPTA (Ru(II)-arene-1,3,5-triaza-7-phosphatricyclo-[3.3.1.1] decane) complexes with different arene substituents to cathepsin B. The recently developed QM/MM-PBSA approach (QM/MM combined with Poisson-Boltzmann solvent-accessible surface area solvation) has been used to estimate binding affinities. The QM calculations reproduce the antitumour activities of the complexes with a correlation coefficient (r(2)) of 0.35-0.86 after a conformational search. The QM/MM-PBSA method gave a better correlation (r(2) = 0.59) when the protein was fixed to the crystal structure, but more reasonable ligand structures and absolute binding energies were obtained if the protein was allowed to relax, indicating that the ligands are strained when the protein is kept fixed. In addition, the best correlation (r(2) = 0.80) was obtained when only the QM energies were used, which suggests that the MM and continuum solvation energies are not accurate enough to predict the binding of a charged metal complex to a charged protein. Taking into account the protein flexibility by means of MD simulations slightly improves the correlation (r(2) = 0.91), but the absolute energies are still too large and the results are sensitive to the details in the calculations, illustrating that it is hard to obtain stable predictions when full flexible protein is included in the calculations.
|
|
| 45. |
- Cirri, Damiano, et al.
(författare)
-
Computationally enhanced X-ray diffraction analysis of a gold(III) complex interacting with the human telomeric DNA G-quadruplex. Unravelling non-unique ligand positioning
- 2022
-
Ingår i: International Journal of Biological Macromolecules. - : Elsevier BV. - 0141-8130. ; 211, s. 506-513
-
Tidskriftsartikel (refereegranskat)abstract
- The crystal structure of the human telomeric DNA Tel24 G-quadruplex (Tel24 = TAG3(T2AG3)3T) in complex with the novel [AuL] species (with L = 2,4,6-tris(2-pyrimidyl)-1,3,5-triazine - TPymT-α) was solved by a novel joint molecular mechanical (MM)/quantum mechanical (QM) innovative approach. The quantum-refinement crystallographic method (crystallographic refinement enhanced with quantum mechanical calculation) was adapted to treat the [AuL]/G-quadruplex structure, where each gold complex in the binding site was found spread over four equally occupied positions. The four positions were first determined by docking restrained to the crystallographically determined metal ions' coordinates. Then, the quantum refinement method was used to resolve the poorly defined density around the ligands and improve the crystallographic determination, revealing that the binding preferences of this metallodrug toward Tel24 G-quadruplex arise from a combined effect of pyrimidine stacking, metal–guanine interactions and charge–charge neutralizing action of the π-acid triazine. The occurrence of interaction in solution with the Tel24 G-quadruplex DNA was further proved through DNA melting experiments, which showed a slight destabilisation of the quadruplex upon adduct formation.
|
|
| 46. |
- De Kerpel, Jan O A, et al.
(författare)
-
Geometric and Electronic Structure of Co(II)-Substituted Azurin
- 1999
-
Ingår i: The Journal of Physical Chemistry Part B. - : American Chemical Society (ACS). - 1520-5207 .- 1520-6106. ; 103:39, s. 8375-8382
-
Tidskriftsartikel (refereegranskat)abstract
- The molecular and electronic structures of Co(II)-substituted azurin have been investigated using several realistic models of the metal coordination sphere. The geometry of the models was optimized using the hybrid density functional B3LYP method and compared to the structures obtained for similar Cu(II) models. It is found that Co(II) prefers a distorted tetrahedral structure with four strong bonds to two histidine nitrogens, the cysteine sulphur, and the backbone carbonyl group. This is in contrast to Cu(II), where two weak axial bonds to methionine and the backbone oxygen are found, combined with three strong bonds to the histidines and cysteine in the equatorial plane of a trigonal bipyramidal structure. The optimal structure of the models conforms with experimental crystal data, indicating that the active-site structure in these proteins is determined by the preferences of the metal ion and its ligand and not by protein strain. The electronic structure and spectrum of the Co(imidazole) 2(SH)(SH) 2(HCONH 2) + model have been investigated in detail using multiconfigurational second-order perturbation theory based on a complete active-space wavefunction (CASPT2). Nine ligand-field transitions and six S cys → Co charge-transfer transitions have been calculated, and all experimentally observed absorption bands in the absorption spectrum of Co(II) azurin have been assigned. It is shown that the Co-S cys bond is more ionic than the Cu-S cys bond and that this causes the blue shift and weakening of the charge-transfer states in the spectrum of Co(II)-substituted azurin compared to native copper protein.
|
|
| 47. |
- De Kerpel, Jan O A, et al.
(författare)
-
Protein strain in blue copper proteins studied by free energy perturbations
- 1999
-
Ingår i: Proteins: Structure, Function and Genetics. - 0887-3585. ; 36:2, s. 157-174
-
Tidskriftsartikel (refereegranskat)abstract
- Free energy perturbations have been performed on two blue copper proteins, plastocyanin and nitrite reductase. By changing the copper coordination geometry, force constants, and charges, we have estimated the maximum energy with which the proteins may distort the copper coordination sphere. By comparing this energy with the quantum chemical energy cost for the same perturbation on the isolated copper complex, various hypotheses about protein strain have been tested. The calculations show that the protein can only modify the copper-methionine bond length by a modest amount of energy -- <5 kJ/mol--and they lend no support to the suggestion that the quite appreciable difference in the copper coordination geometry encountered in the two proteins is a result of the proteins enforcing different Cu- methionine bond lengths. On the contrary, this bond is very flexible, and neither the geometry nor the electronic structure change appreciably when the bond length is changed. Moreover, the proteins are rather indifferent to the length of this bond. Instead, the Cu(II) coordination geometries in the two proteins represent two distinct minima on the potential surface of the copper ligand sphere, characterized by different electronic structures, a tetragonal, mainly or-bonded, structure in nitrite reductase and a trigonal, π-bonded, structure in plastocyanin. In vacuum, the structures have almost the same energy, and they are stabilized in the proteins by a combination of geometric and electrostatic interactions. Plastocyanin favors the bond lengths and electrostatics of the trigonal structure, whereas in nitrite reductase, the angles are the main discriminating factor.
|
|
| 48. |
- De Kerpel, Jan O A, et al.
(författare)
-
Theoretical study of the structural and spectroscopic properties of stellacyanin
- 1998
-
Ingår i: The Journal of Physical Chemistry Part B. - : American Chemical Society (ACS). - 1520-5207 .- 1520-6106. ; 102:23, s. 4638-4647
-
Tidskriftsartikel (refereegranskat)abstract
- The electronic spectrum of the azurin Met121Gln mutant, a model of the blue copper protein stellacyanin, has been studied by ab initio multiconfigurational second-order perturbation theory (the CASPT2 method), including the effect of the protein and solvent by point charges. The six lowest electronic transitions have been calculated and assigned with an error of less than 2400 cm-1. The ground-state singly occupied orbital is found to be a predominantly π antibonding orbital involving Cu3d and Scys3pπ. However, it also contains a significant amount (18%) of Cu-Scys σ antibonding character. This σ interaction is responsible for the appearance in the absorption spectrum of a band at 460 nm, with a significantly higher intensity than observed for other, axial, type 1 copper proteins (i.e., plastocyanin or azurin). The π-σ mixing is caused by the axial glutamine ligand binding at a much shorter distance to copper than the corresponding methionine ligand in the normal blue copper proteins. This explains why, based on its spectral properties, stellacyanin is classified among the rhombic type 1 copper proteins, although its structure is clearly trigonal, as it is for the axial proteins. Calculations have also been performed on structures where the glutamine model coordinates to the copper ion via the deprotonated N∈ atom instead of the O∈ atom. However, the resulting transition energies do not resemble the experimental spectrum obtained at normal or elevated pH. Thus, the results do not confirm the suggestion that the coordinating atom of glutamine changes at high pH.
|
|
| 49. |
- de Visser, Sam P., et al.
(författare)
-
Computational modelling of oxygenation processes in enzymes and biomimetic model complexes
- 2014
-
Ingår i: Chemical Communications. - 1364-548X. ; 50, s. 262-282
-
Forskningsöversikt (refereegranskat)abstract
- With computational resources becoming more efficient and more powerful and at the same time cheaper, computational methods have become more and more popular for studies on biochemical and biomimetic systems. Although large efforts from the scientific community have gone into exploring the possibilities of computational methods on large biochemical systems, such studies are not without pitfalls and often cannot be routinely done but require expert execution. In this review we summarize and highlight advances in computational methodology and its application to enzymatic and biomimetic model systems. In particular, we emphasize on topical and state-of-the-art methodologies that are able to either reproduce experimental findings, e.g., spectroscopic parameters and rate constants, accurately or give predictions on short-lived intermediates and fast reaction processes in nature. Moreover, we give examples of processes where certain computational methods dramatically fail.
|
|
| 50. |
- Delcey, Mickaël G., 1988-, et al.
(författare)
-
Accurate calculations of geometries and singlet-triplet energy differences for active-site models of [NiFe] hydrogenase
- 2014
-
Ingår i: Physical Chemistry, Chemical Physics - PCCP. - : Royal Society of Chemistry (RSC). - 1463-9076 .- 1463-9084. ; 16:17, s. 7927-7938
-
Tidskriftsartikel (refereegranskat)abstract
- We have studied the geometry and singlet-triplet energy difference of two mono-nuclear Ni2+ models related to the active site in [NiFe] hydrogenase. Multiconfigurational second-order perturbation theory based on a complete active-space wavefunction with an active space of 12 electrons in 12 orbitals, CASPT2(12,12), reproduces experimental bond lengths to within 1 pm. Calculated singlet-triplet energy differences agree with those obtained from coupled-cluster calculations with single, double and (perturbatively treated) triple excitations (CCSD(T)) to within 12 kJ mol(-1). For a bimetallic model of the active site of [NiFe] hydrogenase, the CASPT2(12,12) results were compared with the results obtained with an extended active space of 22 electrons in 22 orbitals. This is so large that we need to use restricted active-space theory (RASPT2). The calculations predict that the singlet state is 48-57 kJ mol(-1) more stable than the triplet state for this model of the Ni-Sl(a) state. However, in the [NiFe] hydrogenase protein, the structure around the Ni ion is far from the square-planar structure preferred by the singlet state. This destabilises the singlet state so that it is only similar to 24 kJ mol(-1) more stable than the triplet state. Finally, we have studied how various density functional theory methods compare to the experimental, CCSD(T), CASPT2, and RASPT2 results. Semi-local functionals predict the best singlet-triplet energy differences, with BP86, TPSS, and PBE giving mean unsigned errors of 12-13 kJ mol(-1) (maximum errors of 25-31 kJ mol(-1)) compared to CCSD(T). For bond lengths, several methods give good results, e. g. TPSS, BP86, and M06, with mean unsigned errors of 2 pm for the bond lengths if relativistic effects are considered.
|
|
