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.     

Be(II) in aqueous solution--an extended ab initio QM/MM MD study

An ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulation at double-zeta restricted Hartree-Fock (RHF) level was performed at 293.15 K, including first and second hydration shell in the QM region to study the structural and dynamical properties of the Be(II)-hydrate in aqueous solution. The first tetrahedrally arranged hydration shell, with the four water molecules located at a mean Be-O distance of 1.61 A, is highly inert with respect to ligand exchange processes. The second shell, however, consisting in average of approximately 9.2 water ligands at a mean Be-O distance of 3.7 A and the third shell at a mean Be-O distance of 5.4 A with approximately 19 ligands rapidly exchange water molecules between them and with the bulk, respectively. Other structural parameters such as radial and angular distribution functions (RDF and ADF) and tilt- and theta-angle distributions were also evaluated. The dynamics of the hydrate were studied in terms of ligand mean residence times (MRTs) and librational and vibrational frequencies. The mean residence times for second shell and third shell ligands were determined as 4.8 and 3.2 ps, respectively. The Be-O stretching frequency of 658 cm(-1), associated with a force constant of 147 N m(-1) could be overestimated but it certainly reflects the exceptional stability of the ion-ligand bond in the first hydration shell.     

Pseudobond ab initio QM/MM approach and its applications to enzyme reactions

This perspective article mainly focuses on the development and applications of a pseudobond ab initio QM/MM approach to study enzyme reactions. The following aspects of methodology development are discussed: the approaches for the QM/MM covalent boundary problem, an efficient iterative optimization procedure, the methods to determine enzyme reaction paths, and the approaches to calculate free energy change in enzyme reactions. Several applications are described to illustrate the capability of the methods. Finally, future directions are discussed.     

Polarization and charge-transfer effects in aqueous solution via ab initio QM/MM simulations

Combined ab initio quantum mechanical and molecular mechanical (QM/MM) simulations coupled with the block-localized wave function energy decomposition (BLW-ED) method have been conducted to study the solvation of two prototypical ionic systems, acetate and methylammonium ions in aqueous solution. Calculations reveal that the electronic polarization between the targeted solutes and water is the primary many-body effect, whereas the charge-transfer term only makes a small fraction of the total solute-solvent interaction energy. In particular, the polarization effect is dominated by the solvent (water) polarization.     

How does the cAMP-dependent protein kinase catalyze the phosphorylation reaction: an ab initio QM/MM study

We have carried out density functional theory QM/MM calculations on the catalytic subunit of cAMP-dependent protein kinase (PKA). The QM/MM calculations indicate that the phosphorylation reaction catalyzed by PKA is mainly dissociative, and Asp166 serves as the catalytic base to accept the proton delivered by the substrate peptide. Among the key interactions in the active site, the Mg(2+) ions, glycine rich loop, and Lys72 are found to stabilize the transition state through electrostatic interactions. On the other hand, Lys168, Asn171, Asp184, and the conserved waters bound to Mg(2+) ions do not directly contribute to lower the energy barrier of the phosphorylation reaction, and possible roles for these residues are proposed. The QM/MM calculations with different QM/MM partition schemes or different initial structures yield consistent results. In addition, we have carried out 12 ns molecular dynamics simulations on both wild type and K168A mutated PKA, respectively, to demonstrate that the catalytic role of Lys168 is to keep ATP and substrate peptide in the near-attack reactive conformation.     

Role of the catalytic triad and oxyanion hole in acetylcholinesterase catalysis: an ab initio QM/MM study

The initial step of the acylation reaction catalyzed by acetylcholinesterase (AChE) has been studied by a combined ab initio quantum mechanical/molecular mechanical (QM/MM) approach. The reaction proceeds through the nucleophilic addition of the Ser203 O to the carbonyl C of acetylcholine, and the reaction is facilitated by simultaneous proton transfer from Ser203 to His447. The calculated potential energy barrier at the MP2(6-31+G) QM/MM level is 10.5 kcal/mol, consistent with the experimental reaction rate. The third residue of the catalytic triad, Glu334, is found to be essential in stabilizing the transition state through electrostatic interactions. The oxyanion hole, formed by peptidic NH groups from Gly121, Gly122, and Ala204, is also found to play an important role in catalysis. Our calculations indicate that, in the AChE-ACh Michaelis complex, only two hydrogen bonds are formed between the carbonyl oxygen of ACh and the peptidic NH groups of Gly121 and Gly122. As the reaction proceeds, the distance between the carbonyl oxygen of ACh and NH group of Ala204 becomes smaller, and the third hydrogen bond is formed both in the transition state and in the tetrahedral intermediate.     

Structure and dynamics of hydrated NH(4) (+): an ab initio QM/MM molecular dynamics simulation

A combined ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulation has been performed to investigate solvation structure and dynamics of NH(4) (+) in water. The most interesting region, the sphere includes an ammonium ion and its first hydration shell, was treated at the Hartree-Fock level using DZV basis set, while the rest of the system was described by classical pair potentials. On the basis of detailed QM/MM simulation results, the solvation structure of NH(4) (+) is rather flexible, in which many water molecules are cooperatively involved in the solvation shell of the ion. Of particular interest, the QM/MM results show fast translation and rotation of NH(4) (+) in water. This phenomenon has resulted from multiple coordination, which drives the NH(4) (+) to translate and rotate quite freely within its surrounding water molecules. In addition, a "structure-breaking" behavior of the NH(4) (+) is well reflected by the detailed analysis on the water exchange process and the mean residence times of water molecules surrounding the ion.     

A minimal implementation of the AMBER-GAUSSIAN interface for ab initio QM/MM-MD simulation

For applying to a number of theoretical methodologies based on an ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics method connecting AMBER9 with GAUSSIAN03, we have developed an AMBER-GAUSSIAN interface (AG-IF), which can be one of the simplest architectures. In the AG-IF, only a few subroutines addition is necessary to retrieve the QM/MM energy and forces, obtained by GAUSSIAN, for solving a set of Newtonian equations of motion in AMBER. It is, therefore, easy to be modified for individual applications since AG-IF utilizes most of those functions originally equipped not only in AMBER but also in GAUSSIAN. In the present minimal implementation, only AMBER is modified, whereas GAUSSIAN is left unchanged. Moreover, a different method of calculating electrostatic forces of MM atoms interacting with QM region is proposed. Using the AG-IF, we also demonstrate three examples of application: (i) the QM versus MM comparison in the radial distribution function, (ii) the free energy gradient method, and (iii) the charge from interaction energy and forces.     

Mixed ab initio QM/MM modeling using frozen orbitals and tests with alanine dipeptide and tetrapeptide

Methodology is discussed for mixed ab initio quantum mechanics/molecular mechanics modeling of systems where the quantum mechanics (QM) and molecular mechanics (MM) regions are within the same molecule. The ab initio QM calculations are at the restricted Hartree–Fock level using the pseudospectral method of the Jaguar program while the MM part is treated with the OPLS force fields implemented in the IMPACT program. The interface between the QM and MM regions, in particular, is elaborated upon, as it is dealt with by “breaking” bonds at the boundaries and using Boys-localized orbitals found from model molecules in place of the bonds. These orbitals are kept frozen during QM calculations. Results from tests of the method to find relative conformational energies and geometries of alanine dipeptides and alanine tetrapeptides are presented along with comparisons to pure QM and pure MM calculations. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 1468–1494, 1999     

Structure and dynamics of Au+ ion in aqueous solution: ab initio QM/MM MD simulations

Structure and dynamics of hydrated Au(+) have been investigated by means of molecular dynamics simulations based on ab initio quantum mechanical molecular mechanical forces at Hartree-Fock level for the treatment of the first hydration shell. The outer region of the system was described using a newly constructed classical three-body corrected potential. The structure was evaluated in terms of radial and angular distribution functions and coordination number distributions. Water exchange processes between coordination shells and bulk indicate a very labile structure of the first hydration shell whose average coordination number of 4.7 is a mixture of 3-, 4-, 5-, 6-, and 7-coordinated species. Fast water exchange reactions between first and second hydration shell occur, and the second hydration shell is exceptionally large. Therefore, the mean residence time of water molecules in the first hydration shell (5.6 ps/7.5 ps for t*= 0.5 ps/2.0 ps) is shorter than that in the second shell (9.4 ps/21.2 ps for t*= 0.5 ps/2.0 ps), leading to a quite specific picture of a "structure-breaking" effect.     

Reliable treatment of electrostatics in combined QM/MM simulation of macromolecules

A robust approach for dealing with electrostatic interactions for spherical boundary conditions has been implemented in the QM/MM framework. The development was based on the generalized solvent boundary potential (GSBP) method proposed by Im et al. [J. Chem. Phys. 114, 2924 (2001)], and the specific implementation was applied to the self-consistent-charge density-functional tight-binding approach as the quantum mechanics (QM) level, although extension to other QM methods is straightforward. Compared to the popular stochastic boundary-condition scheme, the new protocol offers a balanced treatment between quantum mechanics/molecular mechanics (QM/MM) and MM/MM interactions; it also includes the effect of the bulk solvent and macromolecule atoms outside of the microscopic region at the Poisson-Boltzmann level. The new method was illustrated with application to the enzyme human carbonic anhydrase II and compared to stochastic boundary-condition simulations using different electrostatic treatments. The GSBP-based QM/MM simulations were most consistent with available experimental data, while conventional stochastic boundary simulations yielded various artifacts depending on different electrostatic models. The results highlight the importance of carefully treating electrostatics in QM/MM simulations of biomolecules and suggest that the commonly used truncation schemes should be avoided in QM/MM simulations, especially in simulations that involve extensive conformational samplings. The development of the GSBP-based QM/MM protocol has opened up the exciting possibility of studying chemical events in very complex biomolecular systems in a multiscale framework.     

Ab initio and QM/MM study of electron addition on the disulfide bond in thioredoxin

Thioredoxin controls the intracellular redox potential through a disulfide/dithiol couple. Under conditions of oxidative stress, this protein functions via one-electron exchange, in which formation of the disulfide radical anion occurs. Combined quantum mechanical (QM) and molecular mechanical (MM) calculations using two- and three-level ONIOM schemes were performed on the thioredoxin (Trx) protein of Chlamydomonas reinhardtii in its oxidized-disulfide and one-electron-reduced forms. In both cases, the active site disulfide moiety was described at the MP2(fc)/6-31+G(d) level, and larger regions of varying sizes around the active site were described at the B3LYP/6-31+G(d) level. The remainder of the 112 residues and 33 water molecules of the crystal structure (PDB entry 1EP7) were described by the AMBER force field. Adiabatic electron affinities were calculated for the disulfide bond in all systems. Separate QM or QM/QM calculations were performed on the QM regions to establish the role of the remainder of the protein on the active site properties. The radical anion species becomes more stable as the number of amide groups in the vicinity increases. One-electron reduction potentials were calculated for the small molecule models, and approximated for the protein for which the values are similar to the experimental one (approximately 0 V). This high reduction potential is due to interaction with the charged end of Lys40, as indicated by mutation in silico to norleucine. The inclusion of the protonated Asp30 side chain and a water molecule in the QM region leads to an increase in the electron affinity. Proton transfer from the Asp30 side chain to the Cys39 sulfur in the radical anion species is strongly disfavored. The radical anion is more stable than the protonated form, which is consistent with experimental results.     

The Jahn-Teller effect of the TiIII ion in aqueous solution: extended ab initio QM/MM molecular dynamics simulations.

Combined ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations, including only the first and the first and second hydration shells in the QM region, were performed for TiIII in aqueous solution. The hydration structure of TiIII is discussed in terms of radial distribution functions, coordination-number distributions and several angle distributions. Dynamical properties, such as librational and vibrational motions and TiIII-O vibrations, were evaluated. A fast dynamical Jahn-Teller effect of TiIII(aq) was observed in the QM/MM simulations, in particular when the second hydration shell was included into the QM region. The results justify the computational effort required for the inclusion of the second hydration shell into the QM region and show the importance of this effort for obtaining accurate hydration-shell geometries, dynamical properties, and details of the Jahn-Teller effect.     

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.