Dataset Information

A protocol for preparing explicitly solvated systems for stable molecular dynamics simulations.

ABSTRACT: Before beginning the production phase of molecular dynamics simulations, i.e., the phase that produces the data to be analyzed, it is often necessary to first perform a series of one or more preparatory minimizations and/or molecular dynamics simulations in order to ensure that subsequent production simulations are stable. This is particularly important for simulations with explicit solvent molecules. Despite the preparatory minimizations and simulations being ubiquitous and essential for stable production simulations, there are currently no general recommended procedures to perform them and very few criteria to decide whether the system is capable of producing a stable simulation trajectory. Here, we propose a simple and well-defined ten step simulation preparation protocol for explicitly solvated biomolecules, which can be applied to a wide variety of system types, as well as a simple test based on the system density for determining whether the simulation is stabilized.

SUBMITTER: Roe DR

PROVIDER: S-EPMC7413747 | biostudies-literature | 2020 Aug

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

A protocol for preparing explicitly solvated systems for stable molecular dynamics simulations.

Roe Daniel R DR Brooks Bernard R BR

The Journal of chemical physics 20200801 5

Before beginning the production phase of molecular dynamics simulations, i.e., the phase that produces the data to be analyzed, it is often necessary to first perform a series of one or more preparatory minimizations and/or molecular dynamics simulations in order to ensure that subsequent production simulations are stable. This is particularly important for simulations with explicit solvent molecules. Despite the preparatory minimizations and simulations being ubiquitous and essential for stable ...[more]

PMID: 32770927

Similar Datasets

Project description:We investigated excess electron solvation dynamics in (NH3)n- ammonia clusters in the n = 8-32 size range by performing finite temperature molecular dynamics simulations. In particular, we focused on three possible scenarios. The first case is designed to model electron attachment to small neutral ammonia clusters (n ≤ ∼10) that form hydrogen-bonded chains. The excess electron is bound to the clusters via dipole bound states, and persists with a VDE of ∼160 meV at 100 K for the n = 8 cluster. The coupled nuclear and electronic relaxation is fast (within ∼100 fs) and takes place predominantly by intermolecular librational motions and the intramolecular umbrella mode. The second scenario illustrates the mechanism of excess electron attachment to cold compact neutral clusters of medium size (8 ≤ n ≤ 32). The neutral clusters show increasing tendency with size to bind the excess electron on the surface of the clusters in weakly bound, diffuse, and highly delocalized states. Anionic relaxation trajectories launched from these initial states provide no indication for excess electron stabilization for sizes n < 24. Excess electrons are likely to autodetach from these clusters. The two largest investigated clusters (n = 24 and 32) can accommodate the excess electron in electronic states that are mainly localized on the surface, but they may be partly embedded in the cluster. In the last 500 fs of the simulated trajectories, the VDE of the solvated electron fluctuates around ∼200 meV for n = 24 and ∼500 meV for n = 32, consistent with the values extrapolated from the experimentally observed linear VDE-n-1/3 trend. In the third case, we prepared neutral ammonia cluster configurations, including an n = 48 cluster, that contain possible electron localization sites within the interior of the cluster. Excess electrons added to these clusters localize in cavities with high VDEs up to 1.9 eV for n = 48. The computed VDE values for larger clusters are considerably higher than the experimentally observed photoelectric threshold energy for the solvated electron in bulk ammonia (∼1.4 eV). Molecular dynamics simulations launched from these initial cavity states strongly indicate, however, that these cavity structures exist only for ∼200 fs. During the relaxation the electron leaves the cavity and becomes delocalized, while the cluster loses its initial compactness.

Project description:Alkali (Li(+), Na(+), K(+), Rb(+), and Cs(+)) and halide (F(-), Cl(-), Br(-), and I(-)) ions play an important role in many biological phenomena, roles that range from stabilization of biomolecular structure, to influence on biomolecular dynamics, to key physiological influence on homeostasis and signaling. To properly model ionic interaction and stability in atomistic simulations of biomolecular structure, dynamics, folding, catalysis, and function, an accurate model or representation of the monovalent ions is critically necessary. A good model needs to simultaneously reproduce many properties of ions, including their structure, dynamics, solvation, and moreover both the interactions of these ions with each other in the crystal and in solution and the interactions of ions with other molecules. At present, the best force fields for biomolecules employ a simple additive, nonpolarizable, and pairwise potential for atomic interaction. In this work, we describe our efforts to build better models of the monovalent ions within the pairwise Coulombic and 6-12 Lennard-Jones framework, where the models are tuned to balance crystal and solution properties in Ewald simulations with specific choices of well-known water models. Although it has been clearly demonstrated that truly accurate treatments of ions will require inclusion of nonadditivity and polarizability (particularly with the anions) and ultimately even a quantum mechanical treatment, our goal was to simply push the limits of the additive treatments to see if a balanced model could be created. The applied methodology is general and can be extended to other ions and to polarizable force-field models. Our starting point centered on observations from long simulations of biomolecules in salt solution with the AMBER force fields where salt crystals formed well below their solubility limit. The likely cause of the artifact in the AMBER parameters relates to the naive mixing of the Smith and Dang chloride parameters with AMBER-adapted Aqvist cation parameters. To provide a more appropriate balance, we reoptimized the parameters of the Lennard-Jones potential for the ions and specific choices of water models. To validate and optimize the parameters, we calculated hydration free energies of the solvated ions and also lattice energies (LE) and lattice constants (LC) of alkali halide salt crystals. This is the first effort that systematically scans across the Lennard-Jones space (well depth and radius) while balancing ion properties like LE and LC across all pair combinations of the alkali ions and halide ions. The optimization across the entire monovalent series avoids systematic deviations. The ion parameters developed, optimized, and characterized were targeted for use with some of the most commonly used rigid and nonpolarizable water models, specifically TIP3P, TIP4P EW, and SPC/E. In addition to well reproducing the solution and crystal properties, the new ion parameters well reproduce binding energies of the ions to water and the radii of the first hydration shells.

Dataset Information

A protocol for preparing explicitly solvated systems for stable molecular dynamics simulations.

Publications

A protocol for preparing explicitly solvated systems for stable molecular dynamics simulations.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets