Exploring polymorphism: the case of benzene.
Explore the source record for details and available documents.
Benzene↗
Biomedical subjects
Michele Parrinello
Publications and source records attributed to Michele Parrinello.
Charge localization in DNA fibers.
We study by first-principles molecular dynamics the mechanism of electron hole (positive charge) localization in a laboratory realizable radical cation Z DNA crystal. We find that at room temperature structural deformation does not provide an efficient localization mechanism. Instead, we find evidence for the importance of changes in the protonation state for stabilizing the radical defect.
DNA, Z-Form↗
Assessing the accuracy of metadynamics.
Metadynamics is a powerful technique that has been successfully exploited to explore the multidimensional free energy surface of complex polyatomic systems and predict transition mechanisms in very different fields, ranging from chemistry and solid-state physics to biophysics. We here derive an explicit expression for the accuracy of the methodology and provide a way to choose the parameters of the method in order to optimize its performance.
Journal Article↗
A recipe for the computation of the free energy barrier and the lowest free energy path of concerted reactions.
The recently introduced hills method (Proc. Natl. Acad. Sci. U.S.A. 2002, 99, 12562) is a powerful tool to compute the multidimensional free energy surface of intrinsically concerted reactions. We have extended this method by focusing our attention on localizing the lowest free energy path that connects the stable reactant and product states. This path represents the most probable reaction mechanism, similar to the zero temperature intrinsic reaction coordinate, but also includes finite temperature effects. The transformation of the multidimensional problem to a one-dimensional reaction coordinate allows for accurate convergence of the free energy profile along the lowest free energy path using standard free energy methods. Here we apply the hills method, our lowest free energy path search algorithm, and umbrella sampling to the prototype S(N)2 reaction. The hills method replaces the in many cases difficult problem of finding a good reaction coordinate with choosing relatively simple collective variables, such as the bond lengths of the broken and formed chemical bonds. The second part of the paper presents a guide to using the hills method, in which we test and fine-tune the method for optimal accuracy and efficiency using the umbrella sampling results as a reference.
Journal Article↗
Topological defects and bulk melting of hexagonal ice.
We use classical molecular dynamics combined with the recently developed metadynamics method [Laio, A.; Parrinello, M. Proc. Natl. Acad. Sci. U.S.A. 2002, 99, 20] to study the process of bulk melting in hexagonal ice. Our simulations show that bulk melting is mediated by the formation of topological defects which preserve the coordination of the tetrahedral network. Such defects cluster to form a defective region involving about 50 molecules with a surprisingly long lifetime. The subsequent formation of coordination defects triggers the transition to the liquid state.
Ice↗
Flexible docking in solution using metadynamics.
We apply our recently developed metadynamics method to the docking of ligands on flexible receptors in water solution. This method mimics the real dynamics of a ligand exiting or entering an enzyme and in so doing reconstructs the free energy surface. We apply it to four docking cases: beta-trypsin/benzamidine, beta-trypsin/chlorobenzamidine, immunoglobulin McPC-603/phosphocholine, and cyclin-dependent kinase 2/staurosporine. In every case studied, the method is able to predict the docked geometry and the free energy of docking. Its added value with respect to many other available methods is that it reconstructs the complete free energy surface, including all the relevant minima and the barriers between them.
Benzamidines↗
The influence of temperature and density functional models in ab initio molecular dynamics simulation of liquid water.
The performance of density functional theory methods for the modeling of condensed aqueous systems is hard to predict and validation by ab initio molecular simulation of liquid water is absolutely necessary. In order to assess the reliability of these tests, the effect of temperature on the structure and dynamics of liquid water has been characterized with 16 simulations of 20 ps in the temperature range of 280-380 K. We find a pronounced influence of temperature on the pair correlation functions and on the diffusion constant including nonergodic behavior on the time scale of the simulation in the lower temperature range (which includes ambient temperature). These observations were taken into account in a consistent comparison of a series of density functionals (BLYP, PBE, TPSS, OLYP, HCTH120, HCTH407). All simulations were carried out using an ab initio molecular dynamics approach in which wave functions are represented using Gaussians and the density is expanded in an auxiliary basis of plane waves. Whereas the first three functionals show similar behavior, it is found that the latter three functionals yield more diffusive dynamics and less structure.
Journal Article↗
Azulene-to-naphthalene rearrangement: the Car-Parrinello metadynamics method explores various reaction mechanisms.
We studied the thermal intramolecular and radical rearrangement of azulene to naphthalene by employing a novel metadynamics method based on Car-Parrinello molecular dynamics. We demonstrate that relatively short simulations can provide us with several possible reaction mechanisms for the rearrangement. We show that different choices of the collective coordinates can steer the reaction along different pathways, thus offering the possibility of choosing the most probable mechanism. We consider herein three intramolecular mechanisms and two radical pathways. We found the norcaradiene pathway to be the preferable intramolecular mechanism, whereas the spiran mechanism is the favored radical route. We obtained high activation energies for all the intramolecular pathways (81.5-98.6 kcal mol(-1)), whereas the radical routes have activation energies of 24-39 kcal mol(-1). The calculations have also resulted in elementary steps and intermediates not yet considered. A few attractive features of the metadynamics method in studying chemical reactions are pointed out.
Journal Article↗
Influence of DNA structure on the reactivity of the guanine radical cation.
Oxidative damage of DNA via radical cation formation is a common cause of mutagenesis, cancer and of the physiological changes associated with aging. By using state-of-the-art ab initio molecular dynamics simulations, we study the mechanism that guides the first steps of this process. In the mechanism proposed here, guanine, which among the bases has the lowest oxidation potential, and the phosphate backbone play a crucial role. We found that the rate limiting step is the water protolysis. We illuminate the role of the local environment in considerably lowering the barrier. Of particular relevance in this respect is the role of the phosphate backbone.
Cytosine↗
Correlations among hydrogen bonds in liquid water.
By performing computer simulations of water with the TIP5P potential we show that structures formed by two or more hydrogen bonds affect the dynamical and static properties of water, especially in the vicinity of freezing temperature. In particular, the short time correlation between two coupled hydrogen bonds cannot be predicted assuming the statistical independence of the single hydrogen bonds. This introduces an additional relaxation time of approximately 9 ps close to the freezing point. We also find that the time persistence of structures formed by several hydrogen bonds (the first solvation shell) correlates with the local density, which is smaller around water molecules with a long-living environment.
Computer Simulation↗
A minimum free energy reaction path for the E2 reaction between fluoro ethane and a fluoride ion.
The prototype binuclear elimination (E2) reaction illustrates the mechanism of a large number of biochemical and industrial applied processes but has received surprisingly little attention in theoretical studies compared to, for example, the substitution (SN2) reaction. This is due to its concerted mechanism, which requires an independent description of the three bonds that are being formed or broken. In this work, we have taken the advantage of a new and promising methodology to efficiently sample intrinsically multidimensional free-energy surfaces. We locate the lowest free-energy reaction path in the 3D configurational space and use this finite-temperature intrinsic reaction coordinate in an umbrella sampling scheme to access the temperature contributions to high accuracy. The small increase of the barrier and the decrease of the overall endothermicity for the E2 reaction due to entropic contributions is non-trivial. Moreover, our strategy to efficiently handle multiple reaction coordinates could be a great benefit to many chemistry-related fields, such as enzyme catalysis, reactions in solution, and nucleation processes.
Journal Article↗
Effective force fields for condensed phase systems from ab initio molecular dynamics simulation: a new method for force-matching.
A novel least-squares fitting approach is presented to obtain classical force fields from trajectory and force databases produced by ab initio (e.g., Car-Parrinello) molecular dynamics (MD) simulations. The method was applied to derive effective nonpolarizable three-site force fields for liquid water at ambient conditions from Car-Parrinello MD simulations in the Becke-Lee-Yang-Parr approximation to the electronic density functional theory. The force-matching procedure includes a fit of short-ranged nonbonded forces, bonded forces, and atomic partial charges. The various parameterizations of the water force field differ by an enforced smooth cut-off applied to the short-ranged interaction term. These were obtained by fitting to the trajectory and force data produced by Car-Parrinello MD simulations of systems of 32 and 64 H(2)O molecules. The new water force fields were developed assuming both flexible or rigid molecular geometry. The simulated structural and self-diffusion properties of liquid water using the fitted force fields are in close agreement with those observed in the underlying Car-Parrinello MD simulations. The resulting empirical models compare to experiment much better than many conventional simple point charge (SPC) models. The fitted potential is also shown to combine well with more sophisticated intramolecular potentials. Importantly, the computational cost of the new models is comparable to that for SPC-like potentials.
Journal Article↗
Hydrogen bond driven chemical reactions: Beckmann rearrangement of cyclohexanone oxime into epsilon-caprolactam in supercritical water.
Recent experiments have shown that supercritical water (SCW) has the ability to accelerate and make selective synthetic organic reactions, thus replacing the common but environmentally harmful acid and basic catalysts. In an attempt to understand the intimate mechanism behind this observation, we analyze, via first-principles molecular dynamics, the Beckmann rearrangement of cyclohexanone oxime into epsilon-caprolactam in supercritical water, for which accurate experimental evidence has been reported. Differences in the wetting of the hydrophilic parts of the solute, enhanced by SCW, and the disrupted hydrogen bond network are shown to be crucial in triggering the reaction and in making it selective. Furthermore, the enhanced concentrations of H(+) in SCW play an important role in starting the reaction.
Journal Article↗
Reconstructing the density of states by history-dependent metadynamics.
We present a novel method for the calculation of the energy density of states D(E) for systems described by classical statistical mechanics. The method builds on an extension of a recently proposed strategy that allows the free-energy profile of a canonical system to be recovered within a preassigned accuracy [Proc. Natl. Acad. Sci. U.S.A. 99, 12562 (2002)]]. The method allows a good control over the error on the recovered system entropy. This fact is exploited to obtain D(E) more efficiently by combining measurements at different temperatures. The accuracy and efficiency of the method are tested for the two-dimensional Ising model (up to size 50 x 50) by comparison with both exact results and previous studies. This method is a general one and should be applicable to more realistic model systems.
Journal Article↗
Dispersion corrections to density functionals for water aromatic interactions.
We investigate recently published methods for extending density functional theory to the description of long-range dispersive interactions. In all schemes an empirical correction consisting of a C6r(-6) term is introduced that is damped at short range. The coefficient C6 is calculated either from average molecular or atomic polarizabilities. We calculate geometry-dependent interaction energy profiles for the water benzene cluster and compare the results with second-order Møller-Plesset calculations. Our results indicate that the use of the B3LYP functional in combination with an appropriate mixing rule and damping function is recommended for the interaction of water with aromatics.
Journal Article↗
Microscopic Mechanism of Antibiotics Translocation through a Porin.
OmpF from the outer membrane of Escherichia coli is a general porin considered to be the main pathway for beta-lactam antibiotics. The availability of a high-resolution crystal structure of OmpF and new experimental techniques at the single-molecule level have opened the way to the investigation of the microscopic mechanisms that allow the passage of antibiotics through bacterial pores. We applied molecular dynamics simulations to investigate the translocation process of ampicillin (Amp) through OmpF. Using a recent algorithm capable of accelerating molecular dynamics simulations we have been able to obtain a reaction path for the translocation of Amp through OmpF. The mechanism of passage depends both on the internal degrees of freedom of Amp and on interactions of Amp with OmpF. Understanding this mechanism would help us design more efficient antibiotics and shed light on nature's way of devising channels able to enhance the transport of molecules through membranes.
Algorithms↗
Hydrogen bonding in water.
Computer simulations using density functional theory based ab initio path integral molecular dynamics have been carried out to investigate hydrogen bonding in water under ambient conditions. Structural predictions for both H2O and D2O, which include the effects of zero-point energy, thermal motion, and many body polarization effects, are contrasted with classical simulations that ignore nuclear quantum effects. The calculated effect of H/D isotope substitution on the water structure is much smaller than the difference between the classical and quantum path integral results, and is in excellent agreement with the measured H/D difference data from both neutron and x-ray scattering.
Journal Article↗
Efficient exploration of reactive potential energy surfaces using Car-Parrinello molecular dynamics.
The possibility of observing chemical reactions in ab initio molecular dynamics runs is severely hindered by the short simulation time accessible. We propose a new method for accelerating the reaction process, based on the ideas of the extended Lagrangian and coarse-grained non-Markovian metadynamics. We demonstrate that by this method it is possible to simulate reactions involving complex atomic rearrangements and very large energy barriers in runs of a few picoseconds.
Journal Article↗