A QM/MM study of the complexes formed by aluminum and iron with serum transferrin at neutral and acidic pH

Serum transferrin (sTf) transports iron in serum and internalizes in cells via receptor mediated endocytosis. Additionally, sTf has been identified as the predominant aluminum carrier in serum. Some questions remain unclear about the exact mechanism for the metal release or whether the aluminum and iron show the same binding mode during the entire process. In the present work, simulation techniques at quantum and atomic levels have been employed in order to gain access into a molecular level understanding of the metal-bound sTf complex, and to describe the binding of Al(III) and Fe(III) ions to sTf. First, hybrid quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulations were carried out in order to analyze the dynamics of the aluminum-loaded complex, taking into account the different pH conditions in blood and into the cell. Moreover, the complexes formed by transferrin with Al(III) and Fe(III) were optimized with high level density functional theory (DFT)/MM methods. All these results indicate that the interaction mode of Al(III) and Fe(III) with sTf change upon different pH conditions, and that the coordination of Al(III) and Fe(III) is not equivalent during the metal intake, transport and release processes. Our results emphasize the importance of the pH on the metal binding and release mechanism and suggest that Al(III) can follow the iron pathway to get access into cells, although once there, it may show a different binding mode, leading to a different mechanism for its release.     

QM/MM free energy simulations: recent progress and challenges

Due to the higher computational cost relative to pure molecular mechanical (MM) simulations, hybrid quantum mechanical/molecular mechanical (QM/MM) free energy simulations particularly require a careful consideration of balancing computational cost and accuracy. Here, we review several recent developments in free energy methods most relevant to QM/MM simulations and discuss several topics motivated by these developments using simple but informative examples that involve processes in water. For chemical reactions, we highlight the value of invoking enhanced sampling technique (e.g. replica-exchange) in umbrella sampling calculations and the value of including collective environmental variables (e.g. hydration level) in metadynamics simulations; we also illustrate the sensitivity of string calculations, especially free energy along the path, to various parameters in the computation. Alchemical free energy simulations with a specific thermodynamic cycle are used to probe the effect of including the first solvation shell into the QM region when computing solvation free energies. For cases where high-level QM/MM potential functions are needed, we analyse two different approaches: the QM/MM–MFEP method of Yang and co-workers and perturbative correction to low-level QM/MM free energy results. For the examples analysed here, both approaches seem productive although care needs to be exercised when analysing the perturbative corrections.     

QM/MM investigation of structure and spectroscopic properties of a vanadium-containing peroxidase

We present a comprehensive analysis of the most likely ground state configuration of the resting state of vanadium dependent chloroperoxidase (VCPO) based on quantum mechanics/molecular mechanics (QM/MM) evaluations of ground state properties, UV-vis spectra and NMR chemical shifts. Within the QM/MM framework, density functional theory (DFT) calculations are used to characterize the resting state of VCPO via time-dependent density functional theory (TD-DFT) calculations of electronic excitation energies and NMR chemical shifts. Comparison with available experimental data allows us to determine the most likely protonation state of VCPO, a state which results in a doubly protonated axial oxygen, a site largely stabilized by hydrogen bonds. We found that the bulk of the protein that is beyond the immediate layer surrounding the cofactor, has an important electrostatic effect on the absorption maximum. Through examination of frontier orbitals, we analyze the nature of two bound water molecules and the extent to which relevant residues in the active site influence the spectroscopy calculations.     

QM/MM studies on the catalytic mechanism of Phenylethanolamine N-methyltransferase

Epinephrine is a naturally occurring adrenomedullary hormone that transduces environmental stressors into cardiovascular actions. As the only route in the catecholamine biosynthetic pathway, Phenylethanolamine N-methyltransferase (PNMT) catalyzes the synthesis of epinephrine. To elucidate the detailed mechanism of enzymatic catalysis of PNMT, combined quantum-mechanical/molecular-mechanical (QM/MM) calculations were performed. The calculation results reveal that this catalysis contains three elementary steps: the deprotonation of protonated norepinphrine, the methyl transferring step and deprotonation of the methylated norepinphrine. The methyl transferring step was proved to be the rate-determining step undergoing a SN2 mechanism with an energy barrier of 16.4 kcal/mol. During the whole catalysis, two glutamic acids Glu185 and Glu219 were proved to be loaded with different effects according to the calculations results of the mutants. These calculation results can be used to explain the experimental observations and make a good complementarity for the previous QM study.     

QM/MM study of the mechanism of enzymatic limonene 1,2-epoxide hydrolysis

Limonene 1,2-epoxide hydrolase (LEH) is completely different from those of classic epoxide hydrolases (EHs) which catalyze the hydrolysis of epoxides to vicinal diols. A novel concerted general acid catalysis step involving the Asp101-Arg99-Asp132 triad is proposed to play an important role in the mechanism. Combined quantum-mechanical/molecular-mechanical (QM/MM) calculations gave activation barriers of 16.9 and 25.1 kcal/mol at the B3LYP/6-31G(d,p)//CHARMM level for nucleophilic attack on the more and less substituted epoxide carbons, respectively. Furthermore, the important roles of residues Arg99, Tyr53 and Asn55 on mutated LEH were evaluated by QM/MM-scanned energy mapping. These results may provide an explanation for site-directed mutagenesis.     

Asymmetric biodegradation of the nerve agents Sarin and VX by human dUTPase: chemometrics,molecular docking and hybrid QM/MM calculations

Organophosphorus compounds (OP) nerve agents are among the most toxic chemical substances known. Their toxicity is due to their ability to bind to acetylcholinesterase. Currently, some enzymes, such as phosphotriesterase, human serum paraoxonase 1 and diisopropyl fluorophosphatase, capable of degrading OP, have been characterized. Regarding the importance of bioremediation methods for detoxication of OP, this work aims to study the interaction modes between the human human deoxyuridine triphosphate nucleotidohydrolase (dUTPase) and Sarin and VX, considering their R_p and S_p enantiomers, to evaluate the asymmetric catalysis of those compounds. In previous work, this enzyme has shown good potential to degrade phosphotriesters, and based on this characteristic, we have applied the human dUTPase to the OP degradation. Molecular docking, chemometrics and mixed quantum and molecular mechanics calculations have been employed, showing a good interaction between dUTPase and OP. Two possible reaction mechanisms were tested, and according to our theoretical results, the catalytic degradation of OP by dUTPase can take place via both mechanisms, beyond being stereoselective, that is, dUTPase cleaves one enantiomer preferentially in relation to other. Chemometric techniques provided excellent assistance for performing this theoretical investigation. The dUTPase study shows importance by the fact of it being a human enzyme.

Communicated by Ramaswamy H. Sarma     

Combined QM/MM calculations of active-site vibrations in binding process of P450cam to putidaredoxin

Combined QM/MM calculations of the active-site of cytochrome P450cam have been performed before and after the binding of P450cam to putidaredoxin. The calculations were carried out for both a 5-coordinated and a 6-coordinated active-site of cytochrome P450cam, with either a water molecule or a carbon monoxide molecule as a 6th distal ligand. An experimentally observed increase in the Fe-S stretching frequency that occurs after cytochrome P450cam binds to putidaredoxin, has been reproduced in our study. Experimentally observed changes in the Fe-C and C-O vibration frequencies that occur after binding of both proteins, have also been reproduced in our study. The computed increase of the Fe-S and Fe-C stretching frequencies is correlated with a corresponding decrease of the Fe-S and Fe-C interatomic distances. According to our calculations, for the active-site with carbon monoxide in the triplet electronic state, the binding process increases the spin densities on the iron and sulfur atoms, which changes the Fe-C and C-O stretching frequencies in opposite directions, in agreement with experimental data.     

QM/MM study of D-fructose in aqueous solution

The QM/MM molecular dynamics methodology was applied to the study of the two main D-fructose tautomers present in aqueous solution, beta-D-fructofuranose and beta-D-fructopyranose. The solute was treated at the AM1 semi-empirical level, and for the solvent water molecules we used the TIP3P potential. We analyzed the structure of the water molecules around the hydroxyl groups to explain the differences in sweet taste between the two tautomers.     

Role of the arginine finger in Ras.RasGAP revealed by QM/MM calculations

In the Ras.Ras.GAP complex, hydrolysis of guanosine triphosphate is strongly accelerated GAP as compared to Ras alone. This is largely attributed to the arginine finger R789(GAP) pointing to AlF(x) in the transition state analogue. We performed QM/MM simulations where triphosphate was treated using the quantum mechanical method of density functional theory, while the protein complex and water environment were described classically using MD. Compared to Ras, the crucial electron shift, bond stretching and distortion towards an eclipsed gamma-to-beta orientation are much more pronounced. The arginine finger is shown to act by displacing water out of the binding niche. The resulting enhanced electrostatic field catalyses the cleavage step.     

Binding free energy calculation with QM/MM hybrid methods for Abl-Kinase inhibitor

We report a Quantum mechanics/Molecular Mechanics–Poisson-Boltzmann/ Surface Area (QM/MM-PB/SA) method to calculate the binding free energy of c-Abl human tyrosine kinase by combining the QM and MM principles where the ligand is treated quantum mechanically and the rest of the receptor by classical molecular mechanics. To study the role of entropy and the flexibility of the protein ligand complex in a solvated environment, molecular dynamics calculations are performed using a hybrid QM/MM approach. This work shows that the results of the QM/MM approach are strongly correlated with the binding affinity. The QM/MM interaction energy in our reported study confirms the importance of electronic and polarization contributions, which are often neglected in classical MM-PB/SA calculations. Moreover, a comparison of semi-empirical methods like DFTB-SCC, PM3, MNDO, MNDO-PDDG, and PDDG-PM3 is also performed. The results of the study show that the implementation of a DFTB-SCC semi-empirical Hamiltonian that is derived from DFT gives better results than other methods. We have performed such studies using the AMBER molecular dynamic package for the first time. The calculated binding free energy is also in agreement with the experimentally determined binding affinity for c-Abl tyrosine kinase complex with Imatinib.     

QM/MM and free-energy simulations of deacylation reaction catalysed by sedolisin,a serine-carboxyl peptidase

Quantum mechanical/molecular mechanical free-energy simulations were performed to understand the deacylation reaction catalysed by sedolisin (a serine-carboxyl peptidase) and to elucidate the catalytic mechanism and the role of the active-site residues during the process. The results given here demonstrate that Asp170 may act as a general acid/base catalyst for the deacylation reaction. It is also shown that the electrostatic oxyanion hole interactions involving Asp170 may be less effective in transition state stabilisation for the deacylation step in the sedolisin-catalysed reaction compared to the general acid/base mechanism. The proton transfer processes during the enzyme-catalysed process were examined, and their role in the catalysis was discussed.     

QM/MM study of the insertion of metal ion into protoporphyrin IX by ferrochelatase

Ferrochelatase catalyzes the metallation of protoporphyrin IX in the terminal step of heme biosynthesis. Mutations in the ferrochelatase gene can lead to the disease erythropoietic porphyria. The catalyzing mechanism of ferrochelatase is still not fully understood. In this paper, we have studied the insertion of Fe²⁺ into the protoporphyrin IX ring by Bacillussubtilis ferrochelatase using combined quantum mechanical and molecular mechanics (QM/MM) calculations. Geometries were optimized at the BP86/6-31G∗ level and energies were calculated at the B3LYP/TZVP level. The overall process involves the stepwise displacement of Glu-264, His-183, and a water molecule from Fe²⁺, and the removal of two protons from the porphyrin ring. The rate-determining step is the cleavage of the bond between the oxygen atom of Glu-264 and Fe²⁺, concomitant with the formation of the first Fe-N bond. It has an energy barrier of 57 kJ mol⁻¹. The porphyrin ring is only slightly distorted in the enzyme active site. The residue Tyr-13 plays a key role for the catalytic process extracting two protons from protoporphyrin IX.     

Characterization of a complete cycle of acetylcholinesterase catalysis by <Emphasis Type="Italic">ab initio</Emphasis> QM/MM modeling

The reaction mechanism of acetylcholine hydrolysis by acetylcholinesterase, including both acylation and deacylation stages from the enzyme-substrate (ES) to the enzyme-product (EP) molecular complexes, is examined by using an ab initio type quantum mechanical – molecular mechanical (QM/MM) approach. The density functional theory PBE0/aug-6–31+G* method for a fairly large quantum part trapped inside the native protein environment, and the AMBER force field parameters in the molecular mechanical part are employed in computations. All reaction steps, including the formation of the first tetrahedral intermediate (TI1), the acylenzyme (EA) complex, the second tetrahedral intermediate (TI2), and the EP complex, are modeled at the same theoretical level. In agreement with the experimental rate constants, the estimated activation energy barrier of the deacylation stage is slightly higher than that for the acylation phase. The critical role of the non-triad Glu202 amino acid residue in orienting lytic water molecule and in stabilizing the second tetrahedral intermediate at the deacylation stage of the enzymatic process is demonstrated. Figure The computed energy diagram for the reaction path from the enzyme – substrate complex (ES) to the enzyme-product complex (EP).     

A comparative study of galactose oxidase and active site analogs based on QM/MM Car-Parrinello simulations

A parallel study of the radical copper enzyme galactose oxidase (GOase) and a low molecular weight analog of the active site was performed with dynamical density functional and mixed quantum-classical calculations. This combined approach enables a direct comparison of the properties of the biomimetic and the natural systems throughout the course of the catalytic reaction. In both cases, five essential forms of the catalytic cycle have been investigated: the resting state in its semi-reduced (catalytically inactive) and its oxidized (catalytically active) form, A ^semi and A ^ox, respectively; a protonated intermediate B; the transition state for the rate-determining hydrogen abstraction step C, and its product D. For A and B the electronic properties of the biomimetic compound are qualitatively very similar to the ones of the natural target. However, in agreement with the experimentally observed difference in catalytic activity, the calculated activation energy for the hydrogen abstraction step is distinctly lower for GOase (16 kcal/mol) than for the mimetic compound (21 kcal/mol). The enzymatic transition state is stabilized by a delocalization of the unpaired spin density over the sulfur-modified equatorial tyrosine Tyr272, an effect that for geometric reasons is essentially absent in the biomimetic compound. Further differences between the mimic and its natural target concern the structure of the product of the abstraction step, which is characterized by a weakly coordinated aldehyde complex for the latter and a tightly bound linear complex for the former. Received 14 October 1999 · Accepted: 19 January 2000     

Proton transfer function of carbonic anhydrase: Insights from QM/MM simulations

Recent QM/MM analyses of proton transfer function of human carbonic anhydrase II (CAII) are briefly reviewed. The topics include a preliminary analysis of nuclear quadrupole coupling constant calculations for the zinc ion and more detailed analyses of microscopic pK_a of the zinc-bound water and free energy profile for the proton transfer. From a methodological perspective, our results emphasize that performing sufficient sampling is essential to the calculation of all these quantities, which reflects the well solvated nature of CAII active site. From a mechanistic perspective, our analyses highlight the importance of electrostatics in shaping the energetics and kinetics of proton transfer in CAII for its function. We argue that once the pK_a for the zinc-bound water is modulated to be in the proper range (~ 7.0), proton transfer through a relatively well solvated cavity towards/from the protein surface (His64) does not require any major acceleration. Therefore, although structural details like the length of the water wire between the donor and acceptor groups still may make a non-negligible contribution, our computational results and the framework of analysis suggest that the significance of such “fine-tuning” is likely secondary to the modulation of pK_a of the zinc-bound water. We encourage further experimental analysis with mutation of (charged) residues not in the immediate neighborhood of the zinc ion to quantitatively test this electrostatics based framework; in particular, Φ analysis based on these mutations may shed further light into the relative importance of the classical Grotthus mechanism and the “proton hole” pathway that we have proposed recently for CAII.     

Molecular electrostatic potentials in aromatic substituted 4-hydroxyquino-2-lones: Glycine/NMDA receptor antagonists

Hydroxyquinolone derivatives have proven to be useful for inhibition at the glycine binding site of N-methyl-D-aspartate (NMDA) receptor. In this work the electronic structure, molecular electrostatic potential (MESP) and vibrational characteristics of a set of C₃ substituted 4-hydroxyquino-2-lone (HQ) derivatives, which act as Glycine/NMDA receptor antagonists, have been investigated using the density functional calculations. In the optimized structures a substituent at the C₃ site of HQ tends to adopt a helical structure. MESP investigations reveal that the ligands showing better inhibition activity should possess electron-rich regions extending over the substituent and carbonyl group of HQ. A correlation of inhibitory activity to the molecular electrostatic potential topography at the carbonyl oxygen as well as to the molecular electron density topography turns out to be a significant output of the investigation. Figure Quantam chemical approach has been employed to understand the reactivity of a set of hydroxyquinolone derivatives known for their inhibition activity towards Glycine/NMDA receptor. Molecular electrostatic potential topography has been used as a tool to understand the reactivity pattern