- Award ID(s):
- 1807740
- NSF-PAR ID:
- 10110480
- Date Published:
- Journal Name:
- Molecular simulation
- ISSN:
- 0892-7022
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
Schneidman-Duhovny, Dina (Ed.)We present OpenAWSEM and Open3SPN2, new cross-compatible implementations of coarse-grained models for protein (AWSEM) and DNA (3SPN2) molecular dynamics simulations within the OpenMM framework. These new implementations retain the chemical accuracy and intrinsic efficiency of the original models while adding GPU acceleration and the ease of forcefield modification provided by OpenMM’s Custom Forces software framework. By utilizing GPUs, we achieve around a 30-fold speedup in protein and protein-DNA simulations over the existing LAMMPS-based implementations running on a single CPU core. We showcase the benefits of OpenMM’s Custom Forces framework by devising and implementing two new potentials that allow us to address important aspects of protein folding and structure prediction and by testing the ability of the combined OpenAWSEM and Open3SPN2 to model protein-DNA binding. The first potential is used to describe the changes in effective interactions that occur as a protein becomes partially buried in a membrane. We also introduced an interaction to describe proteins with multiple disulfide bonds. Using simple pairwise disulfide bonding terms results in unphysical clustering of cysteine residues, posing a problem when simulating the folding of proteins with many cysteines. We now can computationally reproduce Anfinsen’s early Nobel prize winning experiments by using OpenMM’s Custom Forces framework to introduce a multi-body disulfide bonding term that prevents unphysical clustering. Our protein-DNA simulations show that the binding landscape is funneled towards structures that are quite similar to those found using experiments. In summary, this paper provides a simulation tool for the molecular biophysics community that is both easy to use and sufficiently efficient to simulate large proteins and large protein-DNA systems that are central to many cellular processes. These codes should facilitate the interplay between molecular simulations and cellular studies, which have been hampered by the large mismatch between the time and length scales accessible to molecular simulations and those relevant to cell biology.more » « less
-
Simulating the dynamics of ions near polarizable nanoparticles (NPs) using coarse-grained models is extremely challenging due to the need to solve the Poisson equation at every simulation timestep. Recently, a molecular dynamics (MD) method based on a dynamical optimization framework bypassed this obstacle by representing the polarization charge density as virtual dynamic variables and evolving them in parallel with the physical dynamics of ions. We highlight the computational gains accessible with the integration of machine learning (ML) methods for parameter prediction in MD simulations by demonstrating how they were realized in MD simulations of ions near polarizable NPs. An artificial neural network–based regression model was integrated with MD simulation and predicted the optimal simulation timestep and optimization parameters characterizing the virtual system with 94.3% success. The ML-enabled auto-tuning of parameters generated accurate dynamics of ions for ≈ 10 million steps while improving the stability of the simulation by over an order of magnitude. The integration of ML-enhanced framework with hybrid Open Multi-Processing / Message Passing Interface (OpenMP/MPI) parallelization techniques reduced the computational time of simulating systems with thousands of ions and induced charges from thousands of hours to tens of hours, yielding a maximum speedup of ≈ 3 from ML-only acceleration and a maximum speedup of ≈ 600 from the combination of ML and parallel computing methods. Extraction of ionic structure in concentrated electrolytes near oil–water emulsions demonstrates the success of the method. The approach can be generalized to select optimal parameters in other MD applications and energy minimization problems.more » « less
-
For decades, explorations with ground state, thermal reactions combined with pseudophase kinetic models and methods for interpreting the results have provided insights into the properties of the different regions of homogeneous association colloids. More recent successful determination of antioxidant (AO) distributions by this approach is providing new insights into AO efficiency in opaque, well‐mixed two‐phase intact emulsions and eliminating the need to separate the phases. The chemical probe reacts with AOs exclusively in the interfacial region of the emulsion, permitting simplification of the kinetic treatment, and determining its distribution between the oil, interfacial, and aqueous regions. AO distributions are obtained from the two partition constants,
and , of the AOs between the oil‐interfacial and aqueous‐interfacial regions, respectively. and values are obtained by fitting the observed rate constant, k obs, versus surfactant concentration profiles with an overall kinetic approach or model we call the “pseudophase chemical kinetic method.” However, because emulsions break up and reform, and reactants and other components diffuse at various time scales within and between the oil, interfacial, and water regions,k obscould also depend on reactant diffusion coefficients. Here we demonstrate that reactant diffusion is generally orders of magnitude faster than most thermal reactions and reactant distributions between the multiple oil, aqueous, and interfacial droplets and regions of emulsions are in dynamic equilibrium throughout the multiphase systems during the time course of the reaction. Thus, kinetic probes are powerful tools for determining structure‐reactivity, for example, the HLB, relationships governing AO distributions and efficiencies in emulsions.Practical applications: The analysis presented here demonstrates that one of the basic assumptions of the pseudophase chemical kinetic model that we have developed has a solid foundation in the properties of emulsions. That is, we can determine the distributions of reactants between oil (O), interfacial (I), and aqueous (W) regions of the emulsions because the diffusivity coefficients of reactants within emulsions are orders of magnitude greater than the rate of the reaction between the antioxidant and the 4‐hexadecylbenzenediazonium probe. Consequently, we can use the same kinetic model in emulsions as we have used in homogeneous microemulsions. The method permits determination of the partition constants of many antioxidants between the O‐I and W‐I regions of the emulsions and from them their distributions. The method provides new insights into the relationships between antioxidant hydrophobic‐lipophilic balance (HLB) and its efficiency in emulsions and a natural explanation for the cut‐off effect observed with increasing antioxidant HLB.Interpreting chemical reactivity in emulsions requires that reactive components
A andB be in dynamic equilibrium, that separate second order rate constants,k defined for the oil, interfacial, and aqueous regions, subscripts O, I, and W, and for simplicity, the volume of each region depends on the added volumes of oil, surfactant, and water. -
Abstract Predicting microbial metabolic rates and emergent biogeochemical fluxes remains challenging due to the many unknown population dynamical, physiological and reaction‐kinetic parameters and uncertainties in species composition. Here, we show that the need for these parameters can be eliminated when population dynamics and reaction kinetics operate at much shorter time scales than physical mixing processes. Such scenarios are widespread in poorly mixed water columns and sediments. In this ‘fast‐reaction‐transport’ (FRT) limit, all that is required for predictions are chemical boundary conditions, the physical mixing processes and reaction stoichiometries, while no knowledge of species composition, physiology or population/reaction kinetic parameters is needed. Using time‐series data spanning years 2001–2014 and depths 180–900 m across the permanently anoxic Cariaco Basin, we demonstrate that the FRT approach can accurately predict the dynamics of major electron donors and acceptors (Pearson
r ≥ 0.9 in all cases). Hence, many microbial processes in this system are largely transport limited and thus predictable regardless of species composition, population dynamics and kinetics. Our approach enables predictions for many systems in which microbial community dynamics and kinetics are unknown. Our findings also reveal a mechanism for the frequently observed decoupling between function and taxonomy in microbial systems. -
As the demand for PET plastic products continues to grow, developing effective processes to reduce their pollution is of critical importance. Pyrolysis, a promising technology to produce lighter and recyclable components from wasted plastic products, has therefore received considerable attention. In this work, the rapid pyrolysis of PET was studied by using reactive molecular dynamics (MD) simulations. Mechanisms for yielding gas species were unraveled, which involve the generation of ethylene and TPA radicals from ester oxygen−alkyl carbon bond dissociation and condensation reactions to consume TPA radicals with the products of long chains containing a phenyl benzoate structure and CO2. As atomistic simulations are typically conducted at the time scale of a few nanoseconds, a high temperature (i.e. >1000 K) is adopted for accelerated reaction events. To apply the results from MD simulations to practical pyrolysis processes, a kinetic model based on a set of ordinary differential equations was established, which is capable of describing the key products of PET pyrolysis as a function of time and temperature. It was further exploited to determine the optimal reaction conditions for low environmental impact. Overall, this study conducted a detailed mechanism study of PET pyrolysis and established an effective kinetic model for the main species. The approach presented herein to extract kinetic information such as detailed kinetic constants and activation energies from atomistic MD simulations can also be applied to related systems such as the pyrolysis of other polymers.more » « less