Molecular dynamics study of high-density liquid water using a modified central-force potential

QM/MM investigations of organic chemistry oriented questions

About 35 years after its first suggestion, QM/MM became the standard theoretical approach to investigate enzymatic structures and processes. The success is due to the ability of QM/MM to provide an accurate atomistic picture of enzymes and related processes. This picture can even be turned into a movie if nuclei-dynamics is taken into account to describe enzymatic processes. In the field of organic chemistry, QM/MM methods are used to a much lesser extent although almost all relevant processes happen in condensed matter or are influenced by complicated interactions between substrate and catalyst. There is less importance for theoretical organic chemistry since the influence of nonpolar solvents is rather weak and the effect of polar solvents can often be accurately described by continuum approaches. Catalytic processes (homogeneous and heterogeneous) can often be reduced to truncated model systems, which are so small that pure quantum-mechanical approaches can be employed. However, since QM/MM becomes more and more efficient due to the success in software and hardware developments, it is more and more used in theoretical organic chemistry to study effects which result from the molecular nature of the environment. It is shown by many examples discussed in this review that the influence can be tremendous, even for nonpolar reactions. The importance of environmental effects in theoretical spectroscopy was already known. Due to its benefits, QM/MM can be expected to experience ongoing growth for the next decade.In the present chapter we give an overview of QM/MM developments and their importance in theoretical organic chemistry, and review applications which give impressions of the possibilities and the importance of the relevant effects. Since there is already a bunch of excellent reviews dealing with QM/MM, we will discuss fundamental ingredients and developments of QM/MM very briefly with a focus on very recent progress. For the applications we follow a similar strategy.

An extended ab initio QM/MM MD approach to structure and dynamics of Zn(II) in aqueous solution

Structural and dynamical properties of Zn(II) in aqueous solution were investigated, based on an ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulation at double-zeta Hartree-Fock quantum mechanical level including the first and second hydration shells into the QM region. The inclusion of the second shell in the QM region resulted in significant changes in the properties of the hydrate. The first shell coordination number was found to be 6, the second shell consists of approximately 14 water molecules. The structural properties were determined in terms of RDF, ADF, tilt and theta angle distributions, while dynamics were characterized by mean ligand residence times, ion-ligand stretching frequencies and the vibrational and librational motions of water ligands.

Dynamics of ionic and hydrophobic solutes in water-methanol mixtures of varying composition

We have carried out a series of molecular-dynamics simulations of water-methanol mixtures containing either an ionic or a neutral atomic solute to investigate the effects of composition of the mixture on the diffusion of these solutes. Altogether, we have considered 17 different systems of varying composition ranging from pure water to pure methanol. The diffusion coefficients of ionic solutes are found to show nonideal behavior with variation of composition of the solvent mixture. The extent of nonideality of the solute diffusion is found to be similar to the nonideality that is observed for the diffusion and orientational relaxation of water and methanol molecules in these mixtures and is attributed to the enhanced stability of the hydrogen bonds and formation of interspecies complexes in the mixtures. The neutral solute shows characteristics of hydrophobic solvation and its diffusion decreases monotonically with increase of methanol concentration. The present simulation results are compared with those of experiments wherever available.

Structure and ultrafast dynamics of liquid water: A quantum mechanics/molecular mechanics molecular dynamics simulations study

A quantum mechanics/molecular mechanics molecular dynamics simulation was performed for liquid water to investigate structural and dynamical properties of this peculiar liquid. The most important region containing a central reference molecule and all nearest surrounding molecules (first coordination shell) was treated by Hartree-Fock (HF), post-Hartree-Fock [second-order Moller-Plesset perturbation theory (MP2)], and hybrid density functional B3LYP [Becke's three parameter functional (B3) with the correlation functional of Lee, Yang, and Parr (LYP)] methods. In addition, another HF-level simulation (2HF) included the full second coordination shell. Site to site interactions between oxygen-oxygen, oxygen-hydrogen, and hydrogen-hydrogen atoms of all ab initio methods were compared to experimental data. The absence of a second peak and the appearance of a shoulder instead in the gO-O graph obtained from the 2HF simulation is notable, as this feature has been observed so far only for pressurized or heated water. Dynamical data show that the 2HF procedure compensates some of the deficiency of the HF one-shell simulation, reducing the difference between correlated (MP2) and HF results. B3LYP apparently leads to too rigid structures and thus to an artificial slow down of the dynamics.

Water under pressure

Water under pressure is investigated by first principles molecular dynamics, with a focus on the changes in hydrogen bonding and the oxygen network in the nondissociative regime. At a pressure of 10 GPa and a temperature of 600 K, which is close to the freezing point, no appreciable molecular dissociation is observed in the simulations. However, the structure of water is substantially altered from that at ambient conditions. The liquid exhibits a much larger coordination of oxygen atoms, an essential weakening of hydrogen bonding, and sizable changes in the density of electronic states close to the Fermi level. Our results provide new structural data for direct comparison with future experiments.

Molecular mechanisms of pressure induced conformational changes in BPTI

We have performed a 800 ps molecular dynamics simulation of bovine pancreatic trypsin inhibitor (BPTI) in water coupled to a pressure bath at 1, 10,000, 15,000, and 20,000 bar. The simulation reproduces quite well the experimental behavior of the protein under high pressure. The protein keeps its globular form, but adopts a different conformation with a very small reduction in volume. Some residues in the hydrophobic core become exposed to water and a large part of the secondary structure of the protein, (60% of the sheet structure and 40% of the helical structure) is denatured between 10 and 15 kbar. This is in good agreement with experimental data (Goossens, K., et al. Eur. J. Biochem, 236:254-262, 1996) that show denaturation of BPTI between 8 and 14 kbar. A further increase of the pressure results in a freezing of the protein as deduced from the large decrease of the mobility of the residues. During the simulation, the normal structure of water changes from an ice Ih-like to an ice VI-like structure, while keeping the liquid state. The driving force of the high pressure induced conformational transition seems be the higher compressibility of the water compared with the protein. This produces a change in the solvent properties and leads to penetration of the solvent into the hydrophobic core.