Quantum Mechanical/molecular Mechanical Molecular Dynamics Research Articles

Tryptophan to Tryptophan Hole Hopping in an Azurin Construct.

Electron transfer (ET) between neutral and cationic tryptophan residues in the azurin construct [ReI(H126)(CO)3(dmp)](W124)(W122)CuI (dmp = 4,7-Me2-1,10-phenanthroline) was investigated by Born-Oppenheimer quantum-mechanics/molecular mechanics/molecular dynamics (QM/MM/MD) simulations. We focused on W124•+ ← W122 ET, which is the middle step of the photochemical hole-hopping process *ReII(CO)3(dmp•-) ← W124 ← W122 ← CuI, where sequential hopping amounts to nearly 10,000-fold acceleration over single-step tunneling (ACS Cent. Sci. 2019, 5, 192-200). In accordance with experiments, UKS-DFT QM/MM/MD simulations identified forward and reverse steps of W124•+ ↔ W122 ET equilibrium, as well as back ET ReI(CO)3(dmp•-) → W124•+ that restores *ReII(CO)3(dmp•-). Strong electronic coupling between the two indoles (≥40 meV in the crossing region) makes the productive W124•+ ← W122 ET adiabatic. Energies of the two redox states are driven to degeneracy by fluctuations of the electrostatic potential at the two indoles, mainly caused by water solvation, with contributions from the protein dynamics in the W122 vicinity. ET probability depends on the orientation of Re(CO)3(dmp) relative to W124 and its rotation diminishes the hopping yield. Comparison with hole hopping in natural systems reveals structural and dynamics factors that are important for designing efficient hole-hopping processes.

The solvation dynamics of CO2 by quantum mechanical molecular dynamics

The solvation dynamics of CO2 in an aqueous solution were investigated using quantum mechanical molecular mechanical molecular dynamics (QM/MM-MD) simulations. It is demonstrated that the formation of H2CO3 occurs through direct reactions between CO2 and nH2O, with extremely high activation barriers in the gas phase. However, in a solution, the activation energy decreases as the number of H2O molecules increases. Specifically, for the CO2 − H2O system, the activation energy is about 32 kcal/mol, while for the CO2 − 2H2O and CO2 − 3H2O systems, it decreases to 28 kcal/mol and 15 kcal/mol, respectively. These findings suggest that the solvation of CO2 in a solution favors a step-wise mechanism.

Photoinduced hole hopping through tryptophans in proteins

Hole hopping through tryptophan/tyrosine chains enables rapid unidirectional charge transport over long distances. We have elucidated structural and dynamical factors controlling hopping speed and efficiency in two modified azurin constructs that include a rhenium(I) sensitizer, Re(His)(CO)3(dmp)+, and one or two tryptophans (W1, W2). Experimental kinetics investigations showed that the two closely spaced (3 to 4 Å) intervening tryptophans dramatically accelerated long-range electron transfer (ET) from CuI to the photoexcited sensitizer. In our theoretical work, we found that time-dependent density-functional theory (TDDFT) quantum mechanics/molecular mechanics/molecular dynamics (QM/MM/MD) trajectories of low-lying triplet excited states of ReI(His)(CO)3(dmp)+-W1(-W2) exhibited crossings between sensitizer-localized (*Re) and charge-separated [ReI(His)(CO)3(dmp•-)/(W1•+ or W2•+)] (CS1 or CS2) states. Our analysis revealed that the distances, angles, and mutual orientations of ET-active cofactors fluctuate in a relatively narrow range in which the cofactors are strongly coupled, enabling adiabatic ET. Water-dominated electrostatic field fluctuations bring *Re and CS1 states to a crossing where *Re(CO)3(dmp)+←W1 ET occurs, and CS1 becomes the lowest triplet state. ET is promoted by solvation dynamics around *Re(CO)3(dmp)+(W1); and CS1 is stabilized by Re(dmp•-)/W1•+ electron/hole interaction and enhanced W1•+ solvation. The second hop, W1•+←W2, is facilitated by water fluctuations near the W1/W2 unit, taking place when the electrostatic potential at W2 drops well below that at W1•+ Insufficient solvation and reorganization around W2 make W1•+←W2 ET endergonic, shifting the equilibrium toward W1•+ and decreasing the charge-separation yield. We suggest that multiscale TDDFT/MM/MD is a suitable technique to model the simultaneous evolution of photogenerated excited-state manifolds.

Mechanism of Nucleotidyltransfer Reaction and Role of Mg2+ Ion in Sugar Nucleotidyltransferases

Sugar nucleotidyl transferases (SNTs) catalyze nucleotidyltransfer reactions to form sugar-nucleotides and pyrophosphate in the presence of two Mg2+ ions (Mg2+A and Mg2+B). We unveil the mechanism and free energetics of nucleotidyl transfer reaction in an SNT called GlmU through hybrid quantum mechanics-molecular mechanics molecular dynamics simulations and free energy calculations. The study identifies the roles of the active site residues and the Mg2+ ions in catalyzing the reaction. Of great significance, we are able to compare the free energy barrier for the reaction with that for the Mg2+-assisted release of the product (i.e., pyrophosphate) into the solution, shedding light on the general mechanistic and kinetic aspects of catalysis by SNTs.

Reaction Coordinate, Free Energy, and Rate of Intramolecular Proton Transfer in Human Carbonic Anhydrase II.

The role of structure and dynamics of an enzyme has been investigated at three different stages of its function including the chemical event it catalyzes. A one-pot computational method has been designed for each of these stages on the basis of classical and/or quantum mechanical-molecular mechanical molecular dynamics and transition path sampling simulations. For a pair of initial and final states A and B separated by a high free-energy barrier, using a two-stage selection process, several collective variables (CVs) are identified that can delineate A and B. However, these CVs are found to exhibit strong cross-coupling over the transition paths. A set of mutually orthogonal order parameters is then derived from these CVs and an optimal reaction coordinate, r, determined applying half-trajectory likelihood maximization along with a Bayesian information criterion. The transition paths are also used to project the multidimensional free energy surface and barrier crossing dynamics along r. The proposed scheme has been applied to the rate-determining intramolecular proton transfer reaction of the well-known enzyme human carbonic anhydrase II. The potential of mean force, F( r), in the absence of the chemical step is found to reproduce earlier results on the equilibrium population of two side-chain orientations of key residue His-64. Estimation of rate constants, k, from mean first passage times for the three different stages of catalysis shows that the rate-determining step of intramolecular proton transfer occurs with k ≃ 1.0 × 106 s-1, in close agreement with known experimental results.

Gauche effects of glucopyranose by QM/MM-MD simulations

Both ab initio cluster calculations and quantum mechanics/molecular mechanics molecular dynamics (QM/MM-MD) simulations were utilized in the conformational study of glucopyranose. The relative populations of gauche–gauche (GG), gauche–trans (GT), and trans–gauche (TG) by QM/MM-MD closely reproduced the experimental observation where the population of TG is much less than that of the others. Radial distribution functions revealed the less solvated TG as compared to GG and GT. The local structural analysis further showed that a bridging water between O6 and O5 of glucopyranose is missing especially in TG. In short, our newly developed QM/MM-MD method yielded a consistent result with previous studies, in which the gauche effect in glucopyranose comes from the efficient solvation of GG and GT by bridging water. In addition, detailed structural information for the gauche effect has been deduced from our simulations.

PUPIL: A Software Integration System for Multi-Scale QM/MM-MD Simulations and Its Application to Biomolecular Systems.

PUPIL (Program for User Package Interfacing and Linking) implements a distinctive multi-scale approach to hybrid quantum mechanical/molecular mechanical molecular dynamics (QM/MM-MD) simulations. Originally developed to interface different external programs for multi-scale simulation with applications in the materials sciences, PUPIL is finding increasing use in the study of complex biological systems. Advanced MD techniques from the external packages can be applied readily to a hybrid QM/MM treatment in which the forces and energy for the QM region can be computed by any of the QM methods available in any of the other external packages. Here, we give a survey of PUPIL design philosophy, main features, and key implementation decisions, with an orientation to biomolecular simulation. We discuss recently implemented features which enable highly realistic simulations of complex biological systems which have more than one active site that must be treated concurrently. Examples are given.

Polyaniline emeraldine salt in the amorphous solid state: polaron versus bipolaron.

The polaronic and bipolaronic forms of polyaniline emeraldine salt (PAni-ES) in the amorphous solid state have been simulated using classical molecular dynamics (MD) and hybrid quantum mechanical/molecular mechanical-molecular dynamics (QM/MM-MD) approaches. It should be remarked that the electronic state of PAni-ES has been theoretically investigated in the gas phase, solution phase, and crystalline state, but this is the first study in the amorphous solid state, which is the most typical for this conducting polymer. MD simulations were carried out using force-field parametrizations explicitly developed for polaronic and bipolaronic models. QM/MM-MD calculations were performed using a quantum mechanical zone defined by four repeat units. In addition of the structural and electronic characteristics of the two forms of PAni-ES, MD and QM/MM-MD simulations indicate that the bipolaronic is the most stable state of amorphous PAni-ES. Complementary studies have been carried out using different experimental techniques. Although the morphology and topography of doped and undoped PAni are very similar, comparison of their UV-vis spectra supports the preference toward the bipolaronic form of PAni-ES.

Solvent effects on the ultrafast nonradiative deactivation mechanisms of thymine in aqueous solution: Excited-state QM/MM molecular dynamics simulations

On-the-fly excited-state quantum mechanics/molecular mechanics molecular dynamics (QM/MM-MD) simulations of thymine in aqueous solution are performed to investigate the role of solvent water molecules on the nonradiative deactivation process. The complete active space second-order perturbation theory (CASPT2) method is employed for a thymine molecule as the QM part in order to provide a reliable description of the excited-state potential energies. It is found that, in addition to the previously reported deactivation pathway involving the twisting of the C-C double bond in the pyrimidine ring, another efficient deactivation pathway leading to conical intersections that accompanies the out-of-plane displacement of the carbonyl group is observed in aqueous solution. Decay through this pathway is not observed in the gas phase simulations, and our analysis indicates that the hydrogen bonds with solvent water molecules play a key role in stabilizing the potential energies of thymine in this additional decay pathway.

Quantum mechanical/molecular mechanical molecular dynamic simulation of zinc(II) ion in water

The hydration structure of the Zn(II) ion has been studied by classical molecular dynamics (CMD) and quantum mechanical/molecular mechanical molecular dynamics (QM/MM-MD) simulations, using newly constructed ab initio two-body potentials and a three-body correction function. The average first-shell hydration structure is a 6-coordinate octahedron, and the structure parameters obtained by the QM/MM-MD simulation are in reasonable agreement with the experimental values. The first-shell angular distribution functions show a satisfactory consistency for both simulations, suggesting that the classical potential functions well describe the O–Zn–O angular distortion in the CMD simulation, while the quantum mechanical treatment of the Zn(II) ion and its first hydration shell in the QM/MM-MD simulation increases the tilting probability of the coordinating H 2O molecules. Dynamic properties such as vibrational and rotational motions of water molecules were evaluated using velocity autocorrelation functions (VACFs).

Lithium(I) ion hydration: a QM/MM-MD study

The hydration structure of the lithium(I) ion has been studied by a combined quantum mechanical/molecular mechanical molecular dynamics (QM/MM-MD) method at the HF and B3LYP level of theory. Water molecules of both the first and the second hydration shell have been included in the QM treatment. The effect of differently sized QM regions is discussed. The structure obtained by B3LYP is quite insensitive to changes in QM size whereas the HF simulation predicts a weak second hydration shell structure with a large QM region.

Molecular dynamics simulations of Hg2+ in aqueous solution including N-body effects

Pair-potential, three-body corrected and combined ab initio quantum mechanics/molecular mechanics molecular dynamics (QM/MM-MD) simulations have been performed for Hg2+ in water. The ion’s hydration structure was evaluated in terms of radial distribution functions, coordination numbers, and angular distributions, together with changes in ligand structure upon binding of water molecules. An octahedral solvate complex with coordination number of 6/22 for first/second shell results from QM/MM-MD simulations, with average Hg2+–O distance 2.42 Å for the first shell. The pair-potential simulation gives strongly wrong results, the inclusion of three-body effects corrects many, but not all of the errors. The hydration enthalpy of −553±3 kcal/mol obtained from the QM/MM-MD simulation seems a reasonable value for dilute solution, when compared with experimental estimations.

Dynamical properties of the water molecules in the hydration shells of Fe(II) and Fe(III) ions: ab initio QM/MM molecular dynamics simulations

Dynamical properties, librational and vibrational motions of water molecules in the first and second hydration shells of the Fe(II) and Fe(III) ion were evaluated by means of velocity autocorrelation functions obtained by combined quantum mechanical/molecular mechanical molecular dynamics (QM/MM-MD) simulations. The frequencies of rotation around three principal axes and the frequencies of intramolecular vibration of the water molecules in the first hydration shells obtained from the simulations are blue-shifted for both ions compared to those observed experimentally for liquid water. The intramolecular geometry of water molecules in the quantum mechanically treated region (ion plus first hydration shell) shows shorter O–H bonds and wider H–O–H angles than the bulk solvent.