Molecular dynamics simulations using temperature-enhanced essential dynamics replica exchange

Today's standard molecular dynamics simulations of moderately sized biomolecular systems at full atomic resolution are typically limited to the nanosecond timescale and therefore suffer from limited conformational sampling. Efficient ensemble-preserving algorithms like replica exchange (REX) may alleviate this problem somewhat but are still computationally prohibitive due to the large number of degrees of freedom involved. Aiming at increased sampling efficiency, we present a novel simulation method combining the ideas of essential dynamics and REX. Unlike standard REX, in each replica only a selection of essential collective modes of a subsystem of interest (essential subspace) is coupled to a higher temperature, with the remainder of the system staying at a reference temperature, T(0). This selective excitation along with the replica framework permits efficient approximate ensemble-preserving conformational sampling and allows much larger temperature differences between replicas, thereby considerably enhancing sampling efficiency. Ensemble properties and sampling performance of the method are discussed using dialanine and guanylin test systems, with multi-microsecond molecular dynamics simulations of these test systems serving as references.     

Molecular dynamics simulations

Molecular dynamics simulations have become a standard tool for the investigation of biomolecules. Simulations are performed of ever bigger systems using more realistic boundary conditions and better sampling due to longer sampling times. Recently, realistic simulations of systems as complex as transmembrane channels have become feasible. Simulations aid our understanding of biochemical processes and give a dynamic dimension to structural data; for example, the transformation of harmless prion protein into the disease-causing agent has been modeled.     

Molecular dynamics simulations of a DMSO/water mixture using the AMBER force field

Due to its protective properties of biological samples at low temperatures and under desiccation, dimethyl sulfoxide (DMSO) in aqueous solutions has been studied widely by many experimental approaches and molecular dynamics (MD) simulations. In the case of the latter, AMBER is among the most commonly used force fields for simulations of biomolecular systems; however, the parameters for DMSO published by Fox and Kollman in 1998 have only been tested for pure liquid DMSO. We have conducted an MD simulation study of DMSO in a water mixture and computed several structural and dynamical properties such as of the mean density, self-diffusion coefficient, hydrogen bonding and DMSO and water ordering. The AMBER force field of DMSO is seen to reproduce well most of the experimental properties of DMSO in water, with the mixture displaying strong and specific water ordering, as observed in experiments and multiple other MD simulations with other non-polarizable force fields.

Open image in new window

Graphical abstract Hydration structure within hydrogen-bonding distance around a DMSOmolecule

     

Molecular dynamics simulations of ion separation in nano-channel water flows using an electric field

ABSTRACT

Removal of undesired substances from water is a field of investigation recently focused at the nanoscale. Towards this direction, molecular dynamics simulations are conducted in this paper to investigate unwanted ion removal in nanochannel flows. The simulation method incorporates a Poiseuille-like water/ion flow system at the nanoscale where an electric field, of various magnitudes in the range of E?=?0.25–1.5?V/Å, is applied perpendicular to the flow, leading anions and cations close to the wall regions, similar to the Capacitive De-Ionization method. The time needed for ions to reach equilibrium, i.e. to flow in the region near the walls while pure water flows in the channel interior, is t?=?1.3?ns when E?=?1.5?V/Å and t?=?4.0?ns when E?=?0.25?V/Å, showing a dependency on the value of the electric field. Calculations on density, velocity, and temperature values report on fluid properties to be used in the proposed desalination configuration and could act as a basis to guide novel technological applications and extend to higher scales.     

Influence of point defects on the electronic properties of boron nitride nanosheets

Density functional theory was utilized to study the electronic properties of boron nitride (BN) sheets, taking into account the presence of defects. The structure considered consisted of a central hexagon surrounded by alternating pentagons (three) and heptagons (three). The isocoronene cluster model with an armchair edge was used with three different chemical compositions. In the first structure, three B–B bonds were formed where one B in the dimer was part of the central hexagon. In the second structure, three N–N–N bonds were formed at the periphery of the cluster, around the central hexagon. In the third structure, three N–N bonds were formed in a similar fashion to the first model. Our results indicated that the third structure was the most stable configuration; this exhibited planar geometry, semiconductor behavior, and ionic character. To explore the effects of doping, we replaced B and N atoms with C atoms, considering different atomic positions in the central hexagon. When an N atom was replaced with a C atom, the new structure was a semiconductor, but when a B atom was replaced with a C atom, the new structure was a semimetal. At the same time, the polarity increased, inducing covalent behavior. Replacing two N atoms with two C atoms also resulted in a semiconductor, while replacing two B atoms with two C atoms yielded a semimetal; in both cases the bonding was covalent. When three B (three N) atoms of the central hexagon were replaced with three C atoms, the new structure exhibited a transition to a conductor (remained a semiconductor) with low polarity. When monovacancies (N) and divacancies (B and N) were inserted into the lattice, the system was transformed into a covalent semiconductor. Finally, the electrostatic potential surface was calculated in order to explore intermolecular properties such as the charge distribution, which showed how the reactivity of the boron nitride sheets was affected by doping and orbital hybridization.     

Molecular dynamics simulations of interaction between sub-bituminous coal and water

The high moisture content of sub-bituminous coal is associated with the interactions between coal and water. Because of complex composition and structure, the graphite surface modified by hydroxyl, carboxyl and carbonyl groups was used to represent the surface model of sub-bituminous coal according to XPS results. Density profiles for oxygen atoms and hydrogen atoms indicate that the coal surface properties affect the structural and dynamic characteristics of the interfacial water molecules. The interfacial water exhibits much more ordering than bulk water. The results of radial distribution functions, mean square displacement and local self-diffusion coefficient for water molecule related to three oxygen moieties confirmed that the water molecules prefer to absorb with carboxylic groups, and adsorption of water molecules at the hydroxy and carbonyl is similar.     

Molecular dynamics simulations of Leu-enkephalin in water and DMSO.

The structure of Leu-enkephalin (L-Enk) and Met-enkephalin (M-Enk) have frequently been studied, in particular by nuclear magnetic resonance spectroscopy. After more than 20 years of research, it was concluded that enkephalins have no preferred structure in aqueous solution, but that they may have in other solvents. We have performed molecular dynamics simulations of zwitterionic L-Enk in water, and zwitterionic as well as neutral L-Enk dimethyl sulfoxide (DMSO). In water the peptide is very flexible, although there seems to be a preference for compact conformations. In DMSO, the peptide forms a clear salt bridge in the zwitterionic form, but has no preferred conformation in the neutral form. This difference in conformation may provide an explanation for measurements in DMSO in which multiple conformations were found to exist. In this paper we introduce a new formulation for a dihedral angle autocorrelation function, and apply it to study side-chain dynamics in L-Enk. We find that the side-chain dynamics of the large Tyr and Phe residues cannot be adequately sampled in 2.0-ns simulations, while this does seem to be possible for the smaller Leu side chain.     

Molecular dynamics simulations of water within models of ion channels.

The transbilayer pores formed by ion channel proteins contain extended columns of water molecules. The dynamic properties of such waters have been suggested to differ from those of water in its bulk state. Molecular dynamics simulations of ion channel models solvated within and at the mouths of their pores are used to investigate the dynamics and structure of intra-pore water. Three classes of channel model are investigated: a) parallel bundles of hydrophobic (Ala20) alpha-helices; b) eight-stranded hydrophobic (Ala10) antiparallel beta-barrels; and c) parallel bundles of amphipathic alpha-helices (namely, delta-toxin, alamethicin, and nicotinic acetylcholine receptor M2 helix). The self-diffusion coefficients of water molecules within the pores are reduced significantly relative to bulk water in all of the models. Water rotational reorientation rates are also reduced within the pores, particularly in those pores formed by alpha-helix bundles. In the narrowest pore (that of the Ala20 pentameric helix bundle) self-diffusion coefficients and reorientation rates of intra-pore waters are reduced by approximately an order of magnitude relative to bulk solvent. In Ala20 helix bundles the water dipoles orient antiparallel to the helix dipoles. Such dipole/dipole interaction between water and pore may explain how water-filled ion channels may be formed by hydrophobic helices. In the bundles of amphipathic helices the orientation of water dipoles is modulated by the presence of charged side chains. No preferential orientation of water dipoles relative to the pore axis is observed in the hydrophobic beta-barrel models.     

Molecular dynamics simulations of xDNA

xDNA is a modified DNA, which contains natural as well as expanded bases. Expanded bases are generated by the addition of a benzene spacer to the natural bases. A set of AMBER force‐field parameters were derived for the expanded bases and the structural dynamics of the xDNA decamer ( xT5 ′ G xT A xC xG C xA xG T3′ ) · ( xA5′ C T xG C G xT A xC A3′) was explored using a 22 ns molecular dynamics simulation in explicit solvent. During the simulation, the duplex retained its Watson‐Crick base‐pairing and double helical structure, with deviations from the starting B‐form geometry towards A‐form; the deviations are mainly in the backbone torsion angles and in the helical parameters. The sugar pucker of the residues were distributed among a variety of modes; C2′ endo, C1′ exo, O4′ endo, C4′ exo, C2′ exo, and C3′ endo. The enhanced stacking interactions on account of the modification in the bases could help to retain the duplex nature of the helix with minor deviations from the ideal geometry. In our simulation, the xDNA showed a reduced minor groove width and an enlarged major groove width in comparison with the NMR structure. Both the grooves are larger than that of standard B‐DNA, but major groove width is larger than that of A‐DNA with almost equal minor groove width. The enlarged groove widths and the possibility of additional hydration in the grooves makes xDNA a potential molecule for various applications. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 351–360, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Molecular dynamics simulations of carrabiose

Molecular mechanics calculations have been performed for the disaccharide carrabiose, one of the repeat units of β-carrageenan, as a general model for the (1→4)-linkage in the carrageenans. An adiabatic conformational energy map for this unsulfated molecule was prepared by constrained energy minimization and compared to a previously reported rigid-residue energy map for the sulfated molecule and to a similar adiabatic map for neocarrabiose, the related (1→3)-linked dimer repeat unit of β-carrageenan. Molecular dynamics simulations of this molecule in vacuo and in an aqueous (TIP3P) solution were calculated, and the observed motions were found to be generally consistent with the vacuum adiabatic energy map. Unlike the case observed in previous simulations of neocarrabiose, little salvation shift in the molecular conformation was observed for carrabiose. From the dynamics, the linkage was observed to be relatively flexible, as has been inferred from experiment on sulfated carrageenan polymers. © 1997 John Wiley & Sons, Inc.     

Molecular dynamics simulations of metalloproteins

Molecular dynamics simulations are now commonly applied to metalloproteins, despite the challenges introduced by the presence of metal ions. Force field parameters are nowadays available also for these 'exotic' atoms and several biological systems have been successfully studied. Some of the most relevant results and methodological advancements are reviewed.     

Adsorption and immobilisation of human insulin on graphene monoxide,silicon carbide and boron nitride nanosheets investigated by molecular dynamics simulation

The adsorption and immobilisation of human insulin onto the bio-compatible nanosheets including graphene monoxide, silicon carbide and boron nitride nanosheets were studied by molecular dynamics simulation at the temperature of 310 K. After equilibration, heating and 100 ns production molecular dynamic runs, it was found that the insulin was adsorbed and immobilised onto the considered surfaces in a native folded state. The structural parameters, including root-mean-square deviation and fluctuation, surface accessible solvent area, radius of gyration (R_g) and the distance between the centre of the mass of immobilised protein and the surface of the considered nanosheets, were measured, analysed and discussed. The energetics of the studied systems such as the interaction energy between protein and nanosheet was also measured and addressed. The discussions were centred on the structural and energetic parameters of the protein and nanosheets, including charge density, hydrophobicity, hydrophilicity and residue polarity. The results also showed that the active site of C-termini of chain B played an important role in the adsorption process and this could be helpful in the protection of insulin in its smart delivery and release applications.     

Molecular dynamics simulations in photosynthesis

Photosynthesis Research - Photosynthesis is regulated by a dynamic interplay between proteins, enzymes, pigments, lipids, and cofactors that takes place on a large spatio-temporal scale. Molecular...     

Molecular dynamics insights into protein-glycosaminoglycan systems from microsecond-scale simulations

Heparin is a key player in cell signaling via its physical interactions with protein targets in the extracellular matrix. However, basic molecular level understanding of these highly biologically relevant intermolecular interactions is still incomplete. In this study, for the first time, microsecond-scale MD simulations are reported for a complex between fibroblast growth factor 1 and heparin. We rigorously analyze this molecular system in terms of the conformational space, structural, energetic, and dynamic characteristics. We reveal that the conformational selection mechanism of binding denotes a recognition specificity determinant. We conclude that the length of the simulation could be crucial for evaluation of some of the analyzed parameters. Our data provide novel significant insights into the interactions in the fibroblast growth factor 1 complex with heparin, in particular, and into the physical-chemical nature of protein-glycosaminoglycan systems in general, which have potential applicability for biomaterials development in the area of regenerative medicine.