RedMDStream: Parameterization and Simulation Toolbox for Coarse-Grained Molecular Dynamics Models

Coarse-grained (CG) models in molecular dynamics (MD) are powerful tools to simulate the dynamics of large biomolecular systems on micro- to millisecond timescales. However, the CG model, potential energy terms, and parameters are typically not transferable between different molecules and problems. So parameterizing CG force fields, which is both tedious and time-consuming, is often necessary. We present RedMDStream, a software for developing, testing, and simulating biomolecules with CG MD models. Development includes an automatic procedure for the optimization of potential energy parameters based on metaheuristic methods. As an example we describe the parameterization of a simple CG MD model of an RNA hairpin.     

Phase Transitions in Coarse-Grained Lipid Bilayers Containing Cholesterol by Molecular Dynamics Simulations

Coarse-grained simulations of model membranes containing mixtures of phospholipid and cholesterol molecules at different concentrations and temperatures have been performed. A random mixing without tendencies for segregation or formation of domains was observed on spatial scales corresponding to a few thousand lipids and timescales up to several microseconds. The gel-to-liquid crystalline phase transition is successively weakened with increasing amounts of cholesterol without disappearing completely even at a concentration of cholesterol as high as 60%. The phase transition temperature increases slightly depending on the cholesterol concentration. The gel phase system undergoes a transition with increasing amounts of cholesterol from a solid-ordered phase into a liquid-ordered one. In the solid phase, the amplitude of the oscillations in the radial distribution function decays algebraically with a prefactor that goes to zero at the solid-liquid transition.     

Molecular Dynamics Simulation of Structural Phase Transition of AlPO4 Induced by Pressure

Abstract

The mechanism of pressure-induced phase transition of AlPO₄ has been investigated by means of a molecular dynamics method of constant temperature and pressure. A new crystalline phase with space group C2, which has not yet been experimentally found, appears by an instantaneous compression of 60, 70 and 80 GPa at 300 K. At high temperature (2500 K) and pressure (58 GPa), another new phase of AlPo₄ (y-phase), which is composed of PO₆ and AlO₆ octahedra, has been observed.     

The Structure of Liquid Methanol: A Molecular Dynamics Study Using Three-Site Models

Abstract

The structure of liquid methanol at 298.15 K is investigated by performing molecular dynamics (MD) simulations in NVE ensemble using two 3-site force field models. The simulated structural results are compared with the recent neutron diffraction (ND) results obtained at the partial pair distribution function (pdf) level by employing H/D substitution on the hydroxyl hydrogen, Ho. Overall agreement is found between the simulated and experimental total intermolecular radial distribution functions (rdfs). The ability of the 3-site model simulations to satisfactorily reproduce experimental X—X (X = C, O or H- a methyl hydrogen) intermolecular partial distribution function, dominated by contributions from the methyl group. demonstrates that the methyl group does not participate in any bonding in the liquid. However, a comparison between the simulated and experimental Ho—Ho and X—Ho functions reveals that discrepancies still exist at a quantitative level.     

Study on Photoresist Transfer in Rotational Near-Field Photolithography Using Molecular Dynamics Simulation

Rotational near-field photolithography (RNFP) has a great potential for nanostructure fabrication. However, photoresist may be transferred between head and disk resulting in failure of the photolithography process. Two models of photoresist transfer from a rotating disk to a photolithography head were developed using molecular dynamics simulation. The simulation results revealed that photoresist transfer is mainly caused by the strong attraction between the hydrogen atoms in the photoresist and the chromium of the head. When the distance between disk and head is less than the critical flying height, the transfer amount increases with the decrease of rotational speed and flying height. Maintaining the flying height above the critical value is the most effective way to ensure the non-transfer of photoresist.

     

Rotation-Vibration Coupling in Liquid Water Investigated by Molecular Dynamics Simulation

Abstract

The standard rotation-vibration decoupling approximation for liquid water is tested using a method based on the split of the molecular dipole fluctuation correlation function into pure rotational, pure vibrational and crossed terms. Our results corroborate the total decoupling between vibrational and rotational motions in liquid water at ambient conditions and provide a tool able to test the reliability of this approximation for molecular liquids at different states.     

Potential Therapeutic Agents on Alzheimer's Disease through Molecular Docking and Molecular Dynamics Simulation Study of Plant-Based Compounds

Globally Alzheimer's disease (AD) is a highly complex, heterogeneous, and multifactorial neurological disease. AD is categorized clinically through a steady loss in memory and progressive decline of cognitive function. So far, there is no effective cure is available for the treatment of AD. Here, we identified Plant-based compounds (PBCs) from seven therapeutic plants through pharmacophore and pharmacokinetics approaches. Subsequently, we retrieved 65 AD associated proteins by Text Mining approach .We observed the interactions between 39 PBCs with 65 AD-associated targets by using molecular docking. Further, we carried out Molecular dynamics simulation analysis to predict the steady binding of top drug-target complexes. The entire MD simulation results analysis was evidence that seven drug-target complexes consistently interacted during the in silico experiment. The top complexes were the target CHLE interacted with 2 PBCs (Pseudojujubogenin and Anahygrine), target VDAC1 interacted with Withanolide R, target THOP1 interacted with Withaolide R, target AOFB interacted with 2 PBCs (Nardostachysin and Viscosalactone B), and target ACHE interacted with the drug (12-Deoxywithastramonolide). These PBCs have stably and flexibly interacted at the protein‘s active site region. Our results suggest that these PBCs and targets are potential therapeutic candidates for molecular development in AD.     

Solution and Crystal Molecular Dynamics Simulation Study of m4-Cyanovirin-N Mutants Complexed with Di-Mannose

Cyanovirin-N (CVN) is a highly potent anti-HIV carbohydrate-binding agent that establishes its microbicide activity through interaction with mannose-rich glycoprotein gp120 on the virion surface. The m4-CVN and P51G-m4-CVN mutants represent simple models for studying the high-affinity binding site, B^M. A recently determined 1.35 Å high-resolution structure of P51G-m4-CVN provided details on the di-mannose binding mechanism, and suggested that the Arg-76 and Glu-41 residues are critical components of high mannose specificity and affinity. We performed molecular-dynamics simulations in solution and a crystal environment to study the role of Arg-76. Network analysis and clustering were used to characterize the dynamics of Arg-76. The results of our explicit solvent solution and crystal simulations showed a significant correlation with conformations of Arg-76 proposed from x-ray crystallographic studies. However, the crystal simulation showed that the crystal environment strongly biases conformational sampling of the Arg-76 residue. The solution simulations demonstrated no conformational preferences for Arg-76, which would support its critical role as the residue that locks the ligand in the bound state. Instead, a comparative analysis of trajectories from >50 ns of simulation for two mutants revealed the existence of a very stable eight-hydrogen-bond network between the di-mannose ligand and predominantly main-chain atoms. This network may play a key role in the specific recognition and strong binding of mannose oligomers in CVN and its homologs.     

Investigate the Binding of Catechins to Trypsin Using Docking and Molecular Dynamics Simulation

To explore the inhibitory mechanism of catechins for digestive enzymes, we investigated the binding mode of catechins to a typical digestive enzyme-trypsin and analyzed the structure-activity relationship of catechins, using an integration of molecular docking, molecular dynamics simulation and binding free energy calculation. We found that catechins with different structures bound to a conservative pocket S1 of trypsin, which is comprised of residues 189–195, 214–220 and 225–228. In the trypsin-catechin complexes, Asp189 by forming strong hydrogen bonding, and Gln192, Trp215 and Gly216 through hydrophobic interactions, all significantly contribute to the binding of catechins. The number and the position of hydroxyl and aromatic groups, the structure of stereoisomers, and the orientation of catechins in the binding pocket S1 of trypsin all affect the binding affinity. The binding affinity is in the order of Epigallocatechin gallate (EGCG) > Epicatechin gallate (ECG) > Epicatechin (EC) > Epigallocatechin (EGC), and 2R-3R EGCG shows the strongest binding affinity out of other stereoisomers. Meanwhile, the synergic conformational changes of residues and catechins were also analyzed. These findings will be helpful in understanding the knowledge of interactions between catechins and trypsin and referable for the design of novel polyphenol based functional food and nutriceutical formulas.     

Promotion of Crystal Phase Transitions by Mass-of-cell Control in Molecular Dynamics Simulations: Phase Transitions in Benzene Crystals

Abstract

The promotion of crystal phase transitions in molecular dynamics (MD) simulations was realized by controlling the momentum of the MD cell. It was implemented by increasing the mass or velocity of the MD cell instantaneously during simulations within the framework of the constant-pressure method by Parrinello and Rahman. This method induced phase transitions in benzene crystals which have not been obtained in conventional MD simulations. This method is useful for the global search of stable (and metastable) crystal structures.     

Molecular Dynamics Simulation of Alkali Silicates Based on the Quantum Mechanical Potential Surfaces

Abstract

Two potential parameter sets for alkali silicates were derived on the basis of ab-initio MO calculations. One is a model containing completely ionic alkali (model I), and another is that derived from cluster calculations (model II). These sets were tested against the crystal, glass, and liquid of metasilicates. The model II can reproduce these structures well under constant pressure conditions, and is found to be better than model I as a whole.     

Using a Hydrophobic Contact Potential to Evaluate Native and Near-native Folds Generated by Molecular Dynamics Simulations

There are several knowledge-based energy functions that can distinguish the native fold from a pool of grossly misfolded decoys for a given sequence of amino acids. These decoys, which are typically generated by mounting, or “threading”, the sequence onto the backbones of unrelated protein structures, tend to be non-compact and quite different from the native structure: the root-mean-squared (RMS) deviations from the native are commonly in the range of 15 to 20 Å. Effective energy functions should also demonstrate a similar recognition capability when presented with compact decoys that depart only slightly in conformation from the correct structure (i.e. those with RMS deviations of ∼5 Å or less). Recently, we developed a simple yet powerful method for native fold recognition based on the tendency for native folds to form hydrophobic cores. Our energy measure, which we call the hydrophobic fitness score, is challenged to recognize the native fold from 2000 near-native structures generated for each of five small monomeric proteins. First, 1000 conformations for each protein were generated by molecular dynamics simulation at room temperature. The average RMS deviation of this set of 5000 was 1.5 Å. A total of 323 decoys had energies lower than native; however, none of these had RMS deviations greater than 2 Å. Another 1000 structures were generated for each at high temperature, in which a greater range of conformational space was explored (4.3 Å average RMS deviation). Out of this set, only seven decoys were misrecognized. The hydrophobic fitness energy of a conformation is strongly dependent upon the RMS deviation. On average our potential yields energy values which are lowest for the population of structures generated at room temperature, intermediate for those produced at high temperature and highest for those constructed by threading methods. In general, the lowest energy decoy conformations have backbones very close to native structure. The possible utility of our method for screening backbone candidates for the purpose of modelling by side-chain packing optimization is discussed.     

Exploring the pH-Dependent Substrate Transport Mechanism of FocA Using Molecular Dynamics Simulation

FocA belongs to the formate-nitrate transporter family and plays an essential role in the export and uptake of formate in organisms. According to the available crystal structures, the N-terminal residues of FocA are structurally featureless at physiological conditions but at reduced pH form helices to harbor the cytoplasmic entrance of the substrate permeation pathway, which apparently explains the cessation of electrical signal observed in electrophysiological experiments. In this work, we found by structural analysis and molecular dynamics simulations that those N-terminal helices cannot effectively preclude the substrate permeation. Equilibrium simulations and thermodynamic calculations suggest that FocA is permeable to both formate and formic acid, the latter of which is transparent to electrophysiological studies as an electrically neutral species. Hence, the cease of electrical current at acidic pH may be caused by the change of the transported substrate from formate to formic acid. In addition, the mechanism of formate export at physiological pH is discussed.     

A Molecular Dynamics Simulation Study of Coaxial Stacking in RNA

Abstract

We report on unrestrained molecular dynamics simulations of an RNA tetramer binding to a tetra-nucleotide overhang at the 5′-end of an RNA hairpin (nicked structure) and of the corresponding continuous hairpin with Na⁺ as counterions. The simulations lead to stable structures and in this way a structural model for the coaxially stacked RNA hairpin is generated. The stacking interface in the coaxially stacked nicked hairpin structure is characterized by a reduced twist and shift and a slightly increased propeller twist as compared to the continuous system. This leads to an increased overlap between C22 and G23 in the stacking interface of the nicked structure. In the simulations the continuous RNA hairpin has an almost straight helical axis. On the other hand, the corresponding axis for the nicked structure exhibits a marked kink of 39°. The stacking interface exhibits no increased flexibility as compared to the corresponding base pair step in the continuous structure.     

Peptide Aggregation and Pore Formation in a Lipid Bilayer: A Combined Coarse-Grained and All Atom Molecular Dynamics Study

We present a simulation study where different resolutions, namely coarse-grained (CG) and all-atom (AA) molecular dynamics simulations, are used sequentially to combine the long timescale reachable by CG simulations with the high resolution of AA simulations, to describe the complete processes of peptide aggregation and pore formation by alamethicin peptides in a hydrated lipid bilayer. In the 1-μs CG simulations the peptides spontaneously aggregate in the lipid bilayer and exhibit occasional transitions between the membrane-spanning and the surface-bound configurations. One of the CG systems at t = 1 μs is reverted to an AA representation and subjected to AA simulation for 50 ns, during which water molecules penetrate the lipid bilayer through interactions with the peptide aggregates, and the membrane starts leaking water. During the AA simulation significant deviations from the α-helical structure of the peptides are observed, however, the size and arrangement of the clusters are not affected within the studied time frame. Solid-state NMR experiments designed to match closely the setup used in the molecular dynamics simulations provide strong support for our finding that alamethicin peptides adopt a diverse set of configurations in a lipid bilayer, which is in sharp contrast to the prevailing view of alamethicin oligomers formed by perfectly aligned helical alamethicin peptides in a lipid bilayer.     

Using the GELFREE 8100 Fractionation System for Molecular Weight-Based Fractionation with Liquid Phase Recovery

The GELFREE 8100 Fractionation System is a novel protein fractionation system designed to maximize protein recovery during molecular weight based fractionation. The system is comprised of single-use, 8-sample capacity cartridges and a benchtop GELFREE Fractionation Instrument. During separation, a constant voltage is applied between the anode and cathode reservoirs, and each protein mixture is electrophoretically driven from a loading chamber into a specially designed gel column gel. Proteins are concentrated into a tight band in a stacking gel, and separated based on their respective electrophoretic mobilities in a resolving gel. As proteins elute from the column, they are trapped and concentrated in liquid phase in the collection chamber, free of the gel. The instrument is then paused at specific time intervals, and fractions are collected using a pipette. This process is repeated until all desired fractions have been collected. If fewer than 8 samples are run on a cartridge, any unused chambers can be used in subsequent separations.This novel technology facilitates the quick and simple separation of up to 8 complex protein mixtures simultaneously, and offers several advantages when compared to previously available fractionation methods. This system is capable of fractionating up to 1mg of total protein per channel, for a total of 8mg per cartridge. Intact proteins over a broad mass range are separated on the basis of molecular weight, retaining important physiochemical properties of the analyte. The liquid phase entrapment provides for high recovery while eliminating the need for band or spot cutting, making the fractionation process highly reproducible¹.Download video file.^{(89M, mp4)}     

Using Cross-links to Study Ribosomal Dynamics

Abstract After publications of 3-D models of a static ribosome and its large and small subunits, one of the next tasks is to recognize movable ribosomal elements responsible for mechanical shifts during protein synthesis. Statistic analysis of available cross-linking data allowed us to reveal three well separated groups of motions in the ribosome: I, mean magnitude of 10?; II, most abundant, centered at 20? and of wide dispersion, and III, sparsely populated, with large distances up to 95?. The last group, III, comprises elements, like the L7/12-stalk and the L1- protuberance, that adopt different positions in crystallographic or electron micrographic structures, and neighboring hairpins 88 and 89, indicating mobility. We demonstrate that the cross- linking method can be applied to study ribosomal dynamics, including large-scale functional movements and, in particular, to estimate which structures participate in molecular switches.     

From Molecular to Process Simulation: Novel Approaches to the Prediction of Phase Equilibria and PVT Behavior Based on Molecular/Quantum Mechanics and Molecular Dynamics Simulations

Abstract

Actually, in modern process simulators, more than 75% of the code implemented is dedicated to physical properties estimation, calculation and predictions. Data banks storing pure component parameters and binary interaction parameters for phase equilibrium calculations are extensively used and continuously implemented in actual process simulators. This gives an idea of the important role physical properties availability plays in process simulation.

In this paper we propose a new way for coupling molecular and process simulation. The basic machinery is to resort to molecular/quantum mechanics and molecular dynamics simulation techniques for generating the parameters of some equations of state that will subsequently be used for the prediction of phase equilibria and PVT behavior of small and polymeric molecules as well. This information, in turn, will be used as input in the process simulator, thus creating a final and well-defined bridge between molecular and process simulations in chemical engineering.     

A Molecular Dynamics Simulation Study of Rigid and Non-rigid Hard Dumb-bells

Abstract

We present a comparative study, using molecular dynamics, of systems of diatomic, hard dumb-bell, molecules in which the interatomic distance is either constrained to a fixed value or is allowed to vary freely between preset limits. A significant improvement in simulation effciency can be attained by allowing the bond length to vary. We find that thermodynamic properties, and some time correlation functions, are only slightly affected by the removal of the rigid bond-length constraint. The atomic velocity correlation function responds dramatically at short times to changes in the degree of non-rigidity, but at long times these differences are much less important.