Coupled semiempirical quantum mechanics and molecular mechanics (QM/MM) calculations on the aqueous solvation free energies of ionized molecules

The aqueous solvation free energies of ionized molecules were computed using a coupled quantum mechanical and molecular mechanical (QM/MM) model based on the AM1, MNDO, and PM3 semiempirical molecular orbital methods for the solute molecule and the TIP3P molecular mechanics model for liquid water. The present work is an extension of our model for neutral solutes where we assumed that the total free energy is the sum of components derived from the electrostatic/polarization terms in the Hamiltonian plus an empirical “nonpolar” term. The electrostatic/polarization contributions to the solvation free energies were computed using molecular dynamics (MD) simulation and thermodynamic integration techniques, while the nonpolar contributions were taken from the literature. The contribution to the electrostatic/polarization component of the free energy due to nonbonded interactions outside the cutoff radii used in the MD simulations was approximated by a Born solvation term. The experimental free energies were reproduced satisfactorily using variational parameters from the vdW terms as in the original model, in addition to a parameter from the one-electron integral terms. The new one-electron parameter was required to account for the short-range effects of overlapping atomic charge densities. The radial distribution functions obtained from the MD simulations showed the expected H-bonded structures between the ionized solute molecule and solvent molecules. We also obtained satisfactory results by neglecting both the empirical nonpolar term and the electronic polarization of the solute, i.e., by implementing a nonpolarization model. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 1028–1038, 1999     

Comparative quantum mechanics/molecular mechanics (QM/MM) and density functional theory calculations on the oxo-iron species of taurine/alpha-ketoglutarate dioxygenase

We present here the first quantum mechanical/molecular mechanics (QM/MM) studies of taurine/alpha-ketoglutarate dioxygenase (TauD) enzymes. Our studies are focused on the chemical properties of the oxo-iron species and the effect of the protein environment on its structural and electronic behavior. Although the active site region of TauD is very polar with many key hydrogen bonding interactions and salt bridges, the actual effect of the protein environment on the ordering and relative energies of the possible spin state structures is found to be quite small. Optimized geometries are very close to ones observed with density functional theory models that did not take the protein environment into consideration. The calculations show that protonation of the histidine ligands of iron is essential to reproduce the correct electronic representations of the enzyme. Hydroxylation studies of taurine by the oxo-iron active species predict that it is a very efficient catalyst that reacts with substrates via low reaction barriers.     

Importance of accurate charges in molecular docking: quantum mechanical/molecular mechanical (QM/MM) approach

The extent to which accuracy of electric charges plays a role in protein-ligand docking is investigated through development of a docking algorithm, which incorporates quantum mechanical/molecular mechanical (QM/MM) calculations. In this algorithm, fixed charges of ligands obtained from force field parameterization are replaced by QM/MM calculations in the protein environment, treating only the ligands as the quantum region. The algorithm is tested on a set of 40 cocrystallized structures taken from the Protein Data Bank (PDB) and provides strong evidence that use of nonfixed charges is important. An algorithm, dubbed "Survival of the Fittest" (SOF) algorithm, is implemented to incorporate QM/MM charge calculations without any prior knowledge of native structures of the complexes. Using an iterative protocol, this algorithm is able in many cases to converge to a nativelike structure in systems where redocking of the ligand using a standard fixed charge force field exhibits nontrivial errors. The results demonstrate that polarization effects can play a significant role in determining the structures of protein-ligand complexes, and provide a promising start towards the development of more accurate docking methods for lead optimization applications.     

Simulations of liquid ammonia based on the combined quantum mechanical/molecular mechanical (QM/MM) approach

Two combined quantum mechanics/molecular mechanics (QM/MM) molecular dynamics simulations, namely, HF/MM and B3LYP/MM, have been performed to investigate the local structure and dynamics of liquid ammonia. The most interesting region, a sphere containing a central reference molecule and all its nearest surrounding molecules (first coordination shell), was treated by the Hartree-Fock (HF) and hybrid density functional B3LYP methods, whereas the rest of the system was described by the classical pair potentials. On the basis of both HF and B3LYP methods, it is observed that the hydrogen bonding in this peculiar liquid is weak. The structure and dynamics of this liquid are suggested to be determined by the steric packing effects, rather than by the directional hydrogen bonding interactions. Compared to previous empirical as well as Car-Parrinello (CP) molecular dynamics studies, our QM/MM simulations provide detailed information that is in better agreement with experimental data.     

A QM/MM investigation of the activation and catalytic mechanism of Fe-only hydrogenases

Fe-only hydrogenases are enzymes that catalyze dihydrogen production or oxidation, due to the presence of an unusual Fe(6)S(6) cluster (the so-called H-cluster) in their active site, which is composed of a Fe(2)S(2) subsite, directly involved in catalysis, and a classical Fe(4)S(4) cubane cluster. Here, we present a hybrid quantum mechanical and molecular mechanical (QM/MM) investigation of the Fe-only hydrogenase from Desulfovibrio desulfuricans, in order to unravel key issues regarding the activation of the enzyme from its completely oxidized inactive state (Hoxinact) and the influence of the protein environment on the structural and catalytic properties of the H-cluster. Our results show that the Fe(2)S(2) subcluster in the Fe(II)Fe(II) redox state - which is experimentally observed for the completely oxidized form of the enzyme - binds a water molecule to one of its metal centers. The computed QM/MM energy values for water binding to the diferrous subsite are in fact over 70 kJ mol(-1); however, the affinity toward water decreases by 1 order of magnitude after a one-electron reduction of H(ox)(inact), thus leading to the release of coordinated water from the H-cluster. The investigation of a catalytic cycle of the Fe-only hydrogenase that implies formation of a terminal hydride ion and a di(thiomethyl)amine (DTMA) molecule acting as an acid/base catalyst indicates that all steps have reasonable reaction energies and that the influence of the protein on the thermodynamic profile of H(2) production catalysis is not negligible. QM/MM results show that the interactions between the Fe(2)S(2) subsite and the protein environment could give place to structural rearrangements of the H-cluster functional for catalysis, provided that the bidentate ligand that bridges the iron atoms in the binuclear subsite is actually a DTMA residue.     

Ab initio quantum mechanical charge field (QMCF) molecular dynamics: a QM/MM – MD procedure for accurate simulations of ions and complexes

A new formalism for quantum mechanical / molecular mechanical (QM/MM) dynamics of chemical species in solution has been developed, which does not require the construction of any other potential functions except those for solvent–solvent interactions, maintains all the advantages of large simulation boxes and ensures the accuracy of ab initio quantum mechanics for all forces acting in the chemically most relevant region. Interactions between solute and more distant solvent molecules are incorporated by a dynamically adjusted force field corresponding to the actual molecular configuration of the simulated system and charges derived from the electron distribution in the solvate. The new formalism has been tested with some examples of hydrated ions, for which accurate conventional ab initio QM/MM simulations have been previously performed, and the comparison shows equivalence and in some aspects superiority of the new method. As this simulation procedure does not require any tedious construction of two-and three-body interaction potentials inherent to conventional QM/MM approaches, it opens the straightforward access to ab initio molecular dynamics simulations of any kind of solutes, such as metal complexes and other composite species in solution.     

The elusive oxidant species of cytochrome P450 enzymes: characterization by combined quantum mechanical/molecular mechanical (QM/MM) calculations

The primary oxidant of cytochrome P450 enzymes, Compound I, is hard to detect experimentally; in the case of cytochrome P450(cam), this intermediate does not accumulate in solution during the catalytic cycle even at temperatures as low as 200 K (ref 4). Theory can play an important role in characterizing such elusive species. We present here combined quantum mechanical/molecular mechanical (QM/MM) calculations of Compound I of cytochrome P450(cam) in the full enzyme environment as well as density functional studies of the isolated QM region. The calculations assign the ground state of the species, quantify the effect of polarization and hydrogen bonding on its properties, and show that the protein environment and its specific hydrogen bonding to the cysteinate ligand are crucial for sustaining the Fe-S bond and for preventing the full oxidation of the sulfur.     

A QM/MM study on the catalytic mechanism of pyruvate decarboxylase

Pyruvate decarboxylase (PDC) is a typical thiamin diphosphate (ThDP)-dependent enzyme with widespread applications in industry. Though studies regarding the reaction mechanism of PDC have been reported, they are mainly focused on the formation of ThDP ylide and some elementary steps in the catalytic cycle, studies about the whole catalytic cycle of PDC are still not completed. In these previous studies, a major controversy is whether the key active residues (Glu473, Glu50′, Asp27′, His113′, His114′) are protonated or ionized during the reaction. To explore the catalytic mechanism and the role of key residues in the active site, three whole-enzyme models were considered, and the combined QM/MM calculations on the nonoxidative decarboxylation of pyruvate to acetaldehyde catalyzed by PDC were performed. According to our computational results, the fundamental reaction pathways, the complete energy profiles of the whole catalytic cycle, and the specific role of key residues in the common steps were obtained. It is also found that the same residue with different protonation states will lead to different reaction pathways and energy profiles. The mechanism derived from the model in which the residues (Glu473, Glu50′, Asp27′, His113′, His114′) are in their protonated states is most consistent with experimental observations. Therefore, extreme care must be taken when assigning the protonation states in the mechanism study. Because the experimental determination of protonation state is currently difficult, the combined QM/MM method provides an indirect means for determining the active-site protonation state.     

Solvents effects on the mechanism of cellulose hydrolysis: A QM/MM study

This article reports a combined quantum mechanics/molecular mechanics (QM/MM) investigation on the acid hydrolysis of cellulose in water using two different models, cellobiose and a 40‐unit cellulose chain. The explicitly treated solvent molecules strongly influence the conformations, intramolecular hydrogen bonds, and exoanomeric effects in these models. As these features are largely responsible for the barrier to cellulose hydrolysis, the present QM/MM results for the pathways and reaction intermediates in water are expected to be more realistic than those from a former density functional theory (DFT) study with implicit solvent (CPCM). However, in a qualitative sense, there is reasonable agreement between the DFT/CPCM and QM/MM predictions for the reaction mechanism. Differences arise mainly from specific solute–solvent hydrogen bonds that are only captured by QM/MM and not by DFT/CPCM. © 2015 Wiley Periodicals, Inc.     

A molecular mechanics (MM3) study of fluorinated hydrocarbons

A molecular mechanics study of small saturated hydrocarbons (up to C-6), substituted by up to six fluorines was carried out with the MM3 force field. Perfluorobutane and Teflon were also studied. A parameter set was developed for use in the calculation of bond lengths, bond angles, torsion angles, conformational energies, barriers to rotation, dipole moments, moments of inertia and vibrational frequencies for these compounds. The results are in good agreement with experiment when only one or two fluorines are present, but some rather large discrepancies were noted when the F/H ratio becomes high. These can be taken into account only by using a force field more complicated than MM3. Some of the requirements of such a force field are delineated. Some pertinent ab initio results are also reported in this article.     

Gold(I) in liquid ammonia: ab initio QM/MM molecular dynamics simulation

Structure and dynamics investigations of Au(I) ion in liquid ammonia have been performed by means of a molecular dynamics simulation based on ab initio quantum mechanical/molecular mechanical forces, where the first solvation shell was treated by quantum mechanics at Hartree-Fock level. The outer region of the system was described using a newly constructed classical three-body corrected potential. A rigid structure of the first solvation shell was observed with an average Au-N distance of 2.15 A and a coordination number of 2.0.     

螯合吸附材料PAO/SiO2对重金属离子的螯合吸附行为

将丙烯腈接枝聚合在微米级硅胶微粒表面,经偕胺肟化转变,制得了接枝有聚偕胺肟(PAO)的复合型螯合吸附材料PAO/SiO2。本文重点考察了螯合吸附材料PAO/SiO2对几种重金属离子的螯合吸附行为,深入地研究了吸附机理。研究结果表明,偕胺肟基团与重金属离子之间的静电作用与配位螯合作用的协同,导致PAO/SiO2对重金属离子产生强的螯合吸附作用。在可抑制金属离子水解的pH范围内,介质的pH值越高,PAO/SiO2的螯合吸附能力越强;PAO/SiO2对性质不同的金属离子的吸附性能是有差别的,吸附容量的顺序为Cu2+Ni2+Pb2+Cd2+。     

Characterization of structure and dynamics of an aqueous scandium(iii) ion by an extended ab initio QM/MM molecular dynamics simulation

Hydration structure and dynamics of an aqueous Sc(iii) solution were characterized by means of an extended ab initio quantum mechanical/molecular dynamical (QM/MM) molecular dynamics simulation at Hartree-Fock level. A monocapped trigonal prismatic structure composed of seven water molecules surrounding scandium(iii) ion was proposed by the QM/MM simulation including the quantum mechanical effects for the first and second hydration shells. The mean Sc(iii)-O bond length of 2.14 ? was identified for six prism water molecules with one capping water located at around 2.26 ?, reproducing well the X-ray diffraction data. The Sc(iii)-O stretching frequency of 432 cm(-1) corresponding to a force constant of 130 N m(-1), evaluated from the enlarged QM/MM simulation, is in good agreement with the experimentally determined value of 430 cm(-1) (128 N m(-1)). Various water exchange processes in the second hydration shell of the hydrated Sc(iii) ion predict a mean ligand residence time of 7.3 ps.     

Dynamic adsorption and desorption of heavy metal ions on poly(acrylaminophosphonic-carboxyl-hydrazide) chelating fiber

The dynamic adsorption and desorption properties, including the effect of pH value and flow rate on the adsorption, eluent acidity and volume, eluting velocity and re-use, of Cu(II), Pb(II), Zn(II), Cd(II), Mn(II), Ni(II), Co(II) and Hg(II) ions on the column loaded with poly(acrylaminophosphonic-carboxyl-hydrazide) chelating fiber were investigated. The recovery of Mn(II), Co(II), Cd(II), Ni(II) and Zn(II) ions in the presence of Na, K, Ca and Mg ions was examined. The preconcentration of trace amounts of Mn(II), Co(II), Cd(II), Ni(II) and Zn(II) ions from model solution samples was carried out with satisfactory results. The amount of the metal ions detected after preconcentration and recovery by this technique was basically in agreement with the added amount. The method is rapid, precise and simple. Received: 15 October 1997 / Revised: 17 March 1998 / Accepted: 20 March 1998     

Dynamic adsorption and desorption of heavy metal ions on poly(acrylaminophosphonic-carboxyl-hydrazide) chelating fiber

The dynamic adsorption and desorption properties, including the effect of pH value and flow rate on the adsorption, eluent acidity and volume, eluting velocity and re-use, of Cu(II), Pb(II), Zn(II), Cd(II), Mn(II), Ni(II), Co(II) and Hg(II) ions on the column loaded with poly(acrylaminophosphonic-carboxyl-hydrazide) chelating fiber were investigated. The recovery of Mn(II), Co(II), Cd(II), Ni(II) and Zn(II) ions in the presence of Na, K, Ca and Mg ions was examined. The preconcentration of trace amounts of Mn(II), Co(II), Cd(II), Ni(II) and Zn(II) ions from model solution samples was carried out with satisfactory results. The amount of the metal ions detected after preconcentration and recovery by this technique was basically in agreement with the added amount. The method is rapid, precise and simple. Received: 15 October 1997 / Revised: 17 March 1998 / Accepted: 20 March 1998     

Catalytic mechanism and metal specificity of bacterial peptide deformylase: a density functional theory QM/MM study

Bacterial peptide deformylase (PDF) represents a novel class of mononuclear iron peptidase, and has an intriguing metal preference different from most other metalloproteases. Using a hybrid density functional theory (B3LYP) QM/MM method, we have theoretically investigated its catalytic mechanism and metal specificity by studying both Fe2+-PDF and Zn2+-PDF. In both forms of PDF, the conserved Glu133 residue is protonated in the reactant complex, and acts as a general acid during the reaction. The initial reaction step is the nucleophilic attack of the metal-bound hydroxide on the carbonyl carbon of the substrate. Our calculations indicate that the metal ion in Fe2+-PDF is always pentacoordinated during the reaction process, while that in Zn2+-PDF is only tetrahedrally coordinated and not bound to the substrate in the reactant complex. This difference in their metal coordination is suggested to account for the lower activity of Zn2+-PDF in comparison with Fe2+-PDF.