Effect of cut-off distance used in molecular dynamics simulations on fluid properties

The effect of cut-off distance used in molecular dynamics (MD) simulations on fluid properties was studied systematically in both canonical (NVT) and isothermal–isobaric (NPT) ensembles. Results show that the cut-off distance in the NVT ensemble plays little role in determining the equilibrium structure of fluid if the ensemble has a high density. However, pressures calculated in the same NVT ensembles strongly depend on the cut-off distance used. In the NPT ensemble, cut-off distance plays a key role in determining fluid equilibrium structure, density and self-diffusion coefficient. The characteristic of the radial distribution function of fluid in NPT ensembles depending on the cut-off distance used in MD simulations means that the WCA theory (a perturbation theory developed by Weeks, Chandler and Andersen) is not suitable for NPT ensembles because the assumption (the effect of the attractive force in determining the liquid structure is negligible) used in the WCA theory is not valid. The dependence of fluid properties on the cut-off distance also indicates that using the WCA potential (the repulsive part of the intermolecular potential proposed in the WCA theory) to calculate fluid transport in heterogeneous systems could lead to significant errors or incorrect results.     

Concentration effect of cimetidine with POPC bilayer: a molecular dynamics simulation study

All-atom molecular dynamics simulations have been performed on cimetidine in the presence of a palmitoyloleoylphosphatidylcholine (POPC) bilayer. The free energy profile of a single cimetidine molecule passing across POPC bilayer displays a minimum at the interface of bilayer and water. Ten cimetidine molecules were inserted into POPC bilayer to obtain an 8 mol % drug model, and molecular dynamics results showed that cimetidine molecules reside at the polar region of POPC bilayer with sulphur atoms directing to the hydrophobic region. By comparing the one drug model with 8 mol % drug model, one can see that the central barrier to cross the membrane increases while the free energy in bulk water decreases, indicating that the ability of cimetidine passing across the POPC bilayer weakens at increased concentration. In addition, the free energy minimum shifts closer to the hydrophobic core. Our results indicate that with the increased drug concentration, it is more difficult for cimetidine to enter and pass across POPC bilayer.     

Unphysical nucleation of diamond in the extended cutoff Tersoff potential

In simulations of carbon materials it is common practice to view the coefficients of the cutoff function as free parameters which can be optimised according to the system of interest. This work examines a common modification to the widely used Tersoff potential in which the coefficient of the upper cutoff is increased to improve the properties of amorphous carbon. Using molecular dynamics simulations, we show that this so-called extended cutoff Tersoff model leads to nucleation of diamond nanocrystals during annealing of amorphous carbon. By varying the density of the system, and examining the radial distribution function in conjunction with the modified cutoff function, we demonstrate that this behaviour is unphysical, and does not represent a new pathway for synthesising diamond. Viewed from a broader perspective, this observation provides a cautionary tale against altering the parameters of empirical potentials without fully considering the wider implications.     

Dynamic simulations of stimuli-responsive switching of azobenzene derivatives in self-assembled monolayers: reactive rotation potential and switching functions

Although the force field (FF)-based molecular dynamics (MD) simulation has been widely applied to rationalise the experimental observations and measurements in chemistry, physics, materials and life science for years, traditional FF suffers from the incapability for describing chemical reactions, which are crucial in many important transformation processes. In order to simulate the collective switching process in azobenzene-based self-assemble monolayers on Au(111) surface, reactive MD simulations with alternative FF were implemented. The classic torsion function has been modified to depict the diabatic potential energy curves for cis and trans isomers, respectively. A switching function is further introduced to connect two N = N rotation functions, and the surface hopping between the cis and trans curves is allowed. By using the reactive rotation potential and switching function, the collective effect of numerous reaction centres and the influence of environment on the quantum yield in the complex system were explored at mesoscopic dimension and timescales. The reactive FF may be also applicable for other complicated systems containing stilbene derivatives. Limitation and perspective for further developments for the other complicated reactions are also addressed.     

Langevin dynamics simulations reveal biologically relevant folds arising from the incorporation of a torsional potential

The conformational behaviour of polymer chains has been examined using Langevin dynamics simulation techniques. Polymer chains were modelled as “beads” undergoing Brownian motion in a defined potential that accounted for stretching, bending and solvation energies. As expected, the competition between chain stiffness and solvent interactions was found to yield standard swollen or collapsed configurations in good or poor solvents, respectively. However, when a torsional term was introduced into the model, additional biologically relevant conformations such as helices, sheets, turns and hairpins naturally arose.     

Stabilization of Microtubules by Taxane Diterpenoids: Insight from Docking and MD simulations

Microtubules are formed from the molecules of tubulin, whose dynamics is important for many functions in a cell, the most dramatic of which is mitosis. Taxol is known to interact within a specific site on tubulin and also believed to block cell-cycle progression during mitosis by binding to and stabilizing microtubules. Along with the tremendous potential that taxol has shown as an anticancer drug, clinical problems exist with solubility, toxicity, and development of drug resistance. The crystal structure of taxane diterpenoids, namely, 10, 13-deacetyl-abeo-baccatin-IV (I), 5-acetyl-2-deacetoxydecinnamoyl-taxinine-0.29hydrate (II), 7, 9-dideacetyltaxayuntin (III), and Taxawallin-K (IV), are very similar to the taxol molecule. Considerable attention has been given to such molecules whose archetype is taxol but do not posses long aliphatic chains, to be developed as a substitute for taxol with fewer side effects. In the present work, the molecular docking of these taxane diterpenoids has been carried out with the tubulin alpha-beta dimer (1TUB) and refined microtubule structure (1JFF) using Glide-XP, in order to assess the potential of tubulin binding of these cytotoxic agents. Results show that all the ligands dock into the classical taxol binding site of tubulin. Taxol shows the best binding capabilities. On the basis of docking energy and interactions, apart from taxol, molecule II has a better tendency of binding with 1TUB while molecule I shows better binding capability with bovine tubulin 1JFF. To validate the binding capabilities, molecular dynamics (MD) simulations of the best docked complexes of ligands with 1JFF have been carried out for 15.0 ns using DESMOND. Average RMSD variations and time line study of interactions and contacts indicate that these complexes remain stable during the course of the dynamics. However, taxol and molecule II prevail over other taxoids.

Electronic supplementary material

The online version of this article (doi:10.1007/s10867-014-9369-5) contains supplementary material, which is available to authorized users.     

Assembly of lipoprotein particles revealed by coarse-grained molecular dynamics simulations

High-density lipoproteins (HDL) function as cholesterol transporters, facilitating the removal of excess cholesterol from the body. Due to the heterogeneity of native HDL particles (both in size and shape), the details on how these protein-lipid particles form and the structure they assume in their lipid-associated states are not well characterized. We report here a study of the self-assembly of discoidal HDL particles using coarse-grained (CG) molecular dynamics. The microsecond simulations reveal the self-assembly of HDL particles from disordered protein-lipid complexes to form structures containing many of the features of the generally accepted double-belt model for discoidal HDL particles. HDL assembly is found to proceed in two broad steps, aggregation of proteins and lipids driven by the hydrophobic effect which occurs on a approximately 1 micros time scale, followed by the optimization of the protein structure driven by increasingly specific protein-protein interactions.     

Ligand unbinding pathways from the vitamin D receptor studied by molecular dynamics simulations

Molecular dynamics simulation techniques have been used to study the unbinding pathways of 1α,25-dihydroxyvitamin D₃ from the ligand-binding pocket of the vitamin D receptor (VDR). The pathways observed in a large number of relatively short (<200 ps) random acceleration molecular dynamics (RAMD) trajectories were found to be in fair agreement, both in terms of pathway locations and deduced relative preferences, compared to targeted molecular dynamics (TMD) and streered molecular dynamics simulations (SMD). However, the high-velocity ligand expulsions of RAMD tend to favor straight expulsion trajectories and the observed relative frequencies of different pathways were biased towards the probability of entering a particular exit channel. Simulations indicated that for VDR the unbinding pathway between the H1–H2 loop and the β-sheet between H5 and H6 is more favorable than the pathway located between the H1–H2 loop and H3. The latter pathway has been suggested to be the most likely unbinding path for thyroid hormone receptors (TRs) and a likely path for retinoic acid receptor. Ligand entry/exit through these two pathways would not require displacement of H12 from its agonistic position. Differences in the packing of the H1, H2, H3 and β-sheet region explain the changed relative preference of the two unbinding pathways in VDR and TRs. Based on the crystal structures of the ligand binding domains of class 2 nuclear receptors, whose members are VDR and TRs, this receptor class can be divided in two groups according to the packing of the H1, H2, H3 and β-sheet region. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Studies on the structural stability of rabbit prion probed by molecular dynamics simulations of its wild-type and mutants

Prion diseases are invariably fatal and highly infectious neurodegenerative diseases that affect humans and animals. Rabbits are the only mammalian species reported to be resistant to infection from prion diseases isolated from other species (Vorberg et al., 2003). Fortunately, the NMR structure of rabbit prion (124-228) (PDB entry 2FJ3), the NMR structure of rabbit prion protein mutation S173N (PDB entry 2JOH) and the NMR structure of rabbit prion protein mutation I214V (PDB entry 2JOM) were released recently. This paper studies these NMR structures by molecular dynamics simulations. Simulation results confirm the structural stability of wild-type rabbit prion, and show that the salt bridge between D177 and R163 greatly contributes to the structural stability of rabbit prion protein.     

Potential of mean force of the hepatitis C virus core protein–monoclonal 19D9D6 antibody interaction

Antigen–antibody interactions are critical for understanding antigen–antibody associations in immunology. To shed further light on this question, we studied a dissociation of the 19D9D6-HCV core protein antibody complex structure. However, forced separations in single molecule experiments are difficult, and therefore molecular simulation techniques were applied in our study. The stretching, that is, the distance between the center of mass of the HCV core protein and the 19D9D6 antibody, has been studied using the potential of mean force calculations based on molecular dynamics and the explicit water model. Our simulations indicate that the 7 residues Gly70, Gly72, Gly134, Gly158, Glu219, Gln221 and Tyr314, the interaction region (antibody), and the 14 interprotein molecular hydrogen bonds might play important roles in the antigen–antibody interaction, and this finding may be useful for protein engineering of this antigen–antibody structure. In addition, the 3 residues Gly134, Gly158 and Tyr314 might be more important in the development of bioactive antibody analogs.     

Molecular dynamics (MD) simulations for the prediction of chiral discrimination of N-acetylphenylalanine enantiomers by cyclomaltoheptaose (beta-cyclodextrin, beta-CD) based on the MM-PBSA (molecular mechanics-Poisson-Boltzmann surface area) approach

Molecular dynamics (MD) simulations were performed for the prediction of chiral discrimination of N-acetylphenylalanine enantiomers by cyclomaltoheptaose (beta-cyclodextrin, beta-CD). Binding free energies and various conformational properties were obtained using by the MM-PBSA (molecular mechanics Poisson-Boltzmann/surface area) approach. The calculated relative difference (DeltaDeltabinding) of binding free energy was in fine agreement with the experimentally determined value. The difference of rotameric distributions of guest N-acetylphenylalanine enantiomers complexed with the host, beta-CD, was observed after the conformational analyses, suggesting that the conformational changes of guest captured within host cavity would be a decisive factor for enantiodifferentiation at a molecular level.     

Effects of deformability and thermal motion of lipid membrane on electroporation: by molecular dynamics simulations

Effects of mechanical properties and thermal motion of POPE lipid membrane on electroporation were studied by molecular dynamics simulations. Among simulations in which specific atoms of lipids were artificially constrained at their equilibrium positions using a spring with force constant of 2.0 kcal/(mol Ų) in the external electric field of 1.4 kcal/(mol Å e), only constraint on lateral motions of lipid tails prohibited electroporation while non-tail parts had little effects. When force constant decreased to 0.2 kcal/(mol Ų) in the position constraints on lipid tails in the external electric field of 2.0 kcal/(mol Å e), water molecules began to enter the membrane. Position constraints of lipid tails allow water to penetrate from both sides of membrane. Thermal motion of lipids can induce initial defects in the hydrophobic core of membrane, which are favorable nucleation sites for electroporation. Simulations at different temperatures revealed that as the temperature increases, the time taken to the initial pore formation will decrease.     

Prediction of the stability of coiled coils using molecular dynamics simulations

We have performed 40–80 ns-long molecular dynamics (MD) simulations of the GCN4 leucine zipper and synthetic coiled coils using the GROMOS96 (43a2) and OPLS-AA force fields, with the aim of predicting coiled coil stability. Starting with an initial configuration of two peptides placed in an ideal coiled coil configuration, we find that changing the amino acid sequence modestly or decreasing peptide length can lead to a decrease in the final α-helicity of coiled coils, although for peptides as long or longer than 16 residues, the values of helicity do not decrease to the low values seen in the experimental results of Lumb et al. (Biochemistry. 1994, 33, 7361–7367) or of Su et al. (Biochemistry. 1994, 33, 15501–15510), presumably because the simulations are not long enough. We find, however, that helicity correlates positively with the number of close hydrophobic interactions between the two peptides, showing that stable coiled coils in the simulations are tightly packed. The minimum interhelical distances are 0.50–0.66 nm for charged groups, indicating that favorable charge interactions are also important for the stability of the coiled coil.     

Towards a quantitative rationalization of multicomponent glass properties by means of molecular dynamics simulations

This review summarizes the achievements obtained by making use of molecular dynamics (MD) simulations in the elucidation of the structure of multicomponent glasses exerting bioactive properties. Emphasis on critical aspects of MD simulations for oxide glasses treatment is given. The potentiality of the quantitative structure–property relationships (QSPR) analysis as a tool for interpretative and predictive purposes is highlighted.     

Molecular dynamics simulations of palmitate entry into the hydrophobic pocket of the fatty acid binding protein

The entry of substrate into the active site is the first event in any enzymatic reaction. However, due to the short time interval between the encounter and the formation of the stable complex, the detailed steps are experimentally unobserved. In the present study, we report a molecular dynamics simulation of the encounter between palmitate molecule and the Toad Liver fatty acid binding protein, ending with the formation of a stable complex resemblance in structure of other proteins of this family. The forces operating on the system leading to the formation of the tight complex are discussed.     

Probing the interactions of the solvated electron with DNA by molecular dynamics simulations: bromodeoxyuridine substituted DNA

Solvated electrons () are produced during water radiolysis and can interact with biological substrates, including DNA. To augment DNA damage, radiosensitizers such as bromo-deoxyuridine (BUdR), often referred to as an “electron affinic radiosensitizer”, are incorporated in place of isosteric thymidine. However, little is known about the primary interactions of with DNA. In the present study we addressed this problem by applying molecular modeling and molecular dynamics (MD) simulations to a system of normal (BUdR·A)-DNA and a hydrated electron, where the excess electron was modeled as a localized (H₂O)₆ anionic cluster. Our goals were to evaluate the suitability of the MD simulations for this application; to characterize the motion of around DNA (e.g., diffusion coefficients); to identify and describe configurational states of close localization to DNA; and to evaluate the structural dynamics of DNA in the presence of . The results indicate that has distinct space-preferences for forming close contacts with DNA and is more likely to interact directly with nucleotides other than BUdR. Several classes of DNA - contact sites, all within the major groove, were distinguished depending on the structure of the intermediate water layer H-bonding pattern (or its absence, i.e., a direct H-bonding of with DNA bases). Large-scale structural perturbations were identified during and after the approached the DNA from the major groove side, coupled with deeper penetration of sodium counterions in the minor groove. Figure A rare configuration showing direct interaction between the solvated electron and DNA, where (yellow) and N7(A16) are H-bonded. The close approach from the major groove side invokes deep Na⁺ (magenta) penetration into the minor DNA groove (Fig. 7a).     

The molecular weight cut-off of microcapsules is determined by the reaction between alginate and polylysine

Mammalian cells encapsulated in alginate-polylysine microcapsules are used as artificial organs in cancer research and in biotechnology. These applications require microcapsules with a reproducible mol. wt. cut-off. The high cost of the polycation, polylysine, requires an efficient preparation procedure. This article shows that the overall reported contact time of 5 minutes at ambient conditions should be increased several times in order to reach a maximal binding between the calcium alginate beads and 0.1% (w/v) polylysine solutions. An increase of the polylysine concentration from 0.0125% to 0.8% (w/v) resulted in a faster maximal binding, but the amount of polylysine bound increased also. Immersion of calcium alginate beads with a diameter of 750 mum, prepared from 1 mL alginate, in 30 mL of a 0.8% (w/v) polylysine solution, resulted in a polylysine spill of more than 89%. The time required to reach a maximal binding was related to the reaction temperature. The interaction zone between calcium alginate beads and fluorescein isothiocyanate-labeled polylysine solutions was visualized with a confocal laser scanning microscope as a function of time. Microcapsules, prepared at 40 degrees C with 0.1% (w/v) polylysine solutions with mol. wts. between 12 and 249.2 kD, were permeable for fluorescein isothiocyanate-labeled dextran, mol. wt. 4.7, but not for 40.5 kD. Higher polylysine concentrations resulted in a membrane with a mol. wt. cut-off lower than 4.7 kD. (c) 1993 John Wiley & Sons, Inc.