Dataset Information

Chemical Reactions in Classical Molecular Dynamics.

ABSTRACT: An algorithm capable of incorporating multi-step reaction mechanisms into atomistic molecular dynamics (MD) simulations using traditional fixed valence force fields is proposed and implemented within the framework of LAMMPS (Large-scale Atomic Molecular Massively Parallel Simulator). This extension, referred to as fix bond/react, enables bonding topology modifications during a running MD simulation using pre- and post-reaction bonding templates to carry out a pre-specified reaction. Candidate reactants are first identified by interatomic separation, followed by the application of a generalized topology matching algorithm to confirm they match the pre-reaction template. This is followed by a topology conversion to match the post-reaction template and a dynamic relaxation to minimize high energy configurations. Two case studies, the condensation polymerization of nylon 6,6 and the formation of a highly-crosslinked epoxy, are simulated to demonstrate the robustness, stability, and speed of the algorithm. Improvements which could increase its utility are discussed.

SUBMITTER: Gissinger JR

PROVIDER: S-EPMC7608055 | biostudies-literature | 2017 Oct

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Chemical Reactions in Classical Molecular Dynamics.

Gissinger Jacob R JR Jensen Benjamin D BD Wise Kristopher E KE

Polymer 20170918

An algorithm capable of incorporating multi-step reaction mechanisms into atomistic molecular dynamics (MD) simulations using traditional fixed valence force fields is proposed and implemented within the framework of LAMMPS (Large-scale Atomic Molecular Massively Parallel Simulator). This extension, referred to as <i>fix bond/react</i>, enables bonding topology modifications during a running MD simulation using pre- and post-reaction bonding templates to carry out a pre-specified reaction. Candi ...[more]

PMID: 33149370

Similar Datasets

Project description:Hybrid quantum/classical molecular dynamics simulations of the two proton transfer reactions catalyzed by ketosteroid isomerase are presented. The potential energy surfaces for the proton transfer reactions are described with the empirical valence bond method. Nuclear quantum effects of the transferring hydrogen increase the rates by a factor of approximately 8, and dynamical barrier recrossings decrease the rates by a factor of 3-4. For both proton transfer reactions, the donor-acceptor distance decreases substantially at the transition state. The carboxylate group of the Asp38 side chain, which serves as the proton acceptor and donor in the first and second steps, respectively, rotates significantly between the two proton transfer reactions. The hydrogen-bonding interactions within the active site are consistent with the hydrogen bonding of both Asp99 and Tyr14 to the substrate. The simulations suggest that a hydrogen bond between Asp99 and the substrate is present from the beginning of the first proton transfer step, whereas the hydrogen bond between Tyr14 and the substrate is virtually absent in the first part of this step but forms nearly concurrently with the formation of the transition state. Both hydrogen bonds are present throughout the second proton transfer step until partial dissociation of the product. The hydrogen bond between Tyr14 and Tyr55 is present throughout both proton transfer steps. The active site residues are more mobile during the first step than during the second step. The van der Waals interaction energy between the substrate and the enzyme remains virtually constant along the reaction pathway, but the electrostatic interaction energy is significantly stronger for the dienolate intermediate than for the reactant and product. Mobile loop regions distal to the active site exhibit significant structural rearrangements and, in some cases, qualitative changes in the electrostatic potential during the catalytic reaction. These results suggest that relatively small conformational changes of the enzyme active site and substrate strengthen the hydrogen bonds that stabilize the intermediate, thereby facilitating the proton transfer reactions. Moreover, the conformational and electrostatic changes associated with these reactions are not limited to the active site but rather extend throughout the entire enzyme.

Project description:For the successful development and application of lubricants, a full understanding of their complex nanoscale behavior under a wide range of external conditions is required, but this is difficult to obtain experimentally. Nonequilibrium molecular dynamics (NEMD) simulations can be used to yield unique insights into the atomic-scale structure and friction of lubricants and additives; however, the accuracy of the results depend on the chosen force-field. In this study, we demonstrate that the use of an accurate, all-atom force-field is critical in order to; (i) accurately predict important properties of long-chain, linear molecules; and (ii) reproduce experimental friction behavior of multi-component tribological systems. In particular, we focus on n-hexadecane, an important model lubricant with a wide range of industrial applications. Moreover, simulating conditions common in tribological systems, i.e., high temperatures and pressures (HTHP), allows the limits of the selected force-fields to be tested. In the first section, a large number of united-atom and all-atom force-fields are benchmarked in terms of their density and viscosity prediction accuracy of n-hexadecane using equilibrium molecular dynamics (EMD) simulations at ambient and HTHP conditions. Whilst united-atom force-fields accurately reproduce experimental density, the viscosity is significantly under-predicted compared to all-atom force-fields and experiments. Moreover, some all-atom force-fields yield elevated melting points, leading to significant overestimation of both the density and viscosity. In the second section, the most accurate united-atom and all-atom force-field are compared in confined NEMD simulations which probe the structure and friction of stearic acid adsorbed on iron oxide and separated by a thin layer of n-hexadecane. The united-atom force-field provides an accurate representation of the structure of the confined stearic acid film; however, friction coefficients are consistently under-predicted and the friction-coverage and friction-velocity behavior deviates from that observed using all-atom force-fields and experimentally. This has important implications regarding force-field selection for NEMD simulations of systems containing long-chain, linear molecules; specifically, it is recommended that accurate all-atom potentials, such as L-OPLS-AA, are employed.

Dataset Information

Chemical Reactions in Classical Molecular Dynamics.

Publications

Chemical Reactions in Classical Molecular Dynamics.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets