Simulations Of Coarse-grained Model Research Articles

Constraint Release in Entangled Binary Blends of Linear Polymers: A Molecular Dynamics Study

We present extensive molecular dynamics simulations of the dynamics of diluted long probe chains entangled with a matrix of shorter chains. The chain lengths of both components are above the entanglement strand length, and the ratio of their lengths is varied over a wide range to cover the crossover from the chain reptation regime to tube Rouse motion regime of the long probe chains. Reducing the matrix chain length results in a faster decay of the dynamic structure factor of the probe chains, in good agreement with recent neutron spin echo experiments. The diffusion of the long chains, measured by the mean square displacements of the monomers and the centers of mass of the chains, demonstrates a systematic speed-up relative to the pure reptation behavior expected for monodisperse melts of sufficiently long polymers. On the other hand, the diffusion of the matrix chains is only weakly perturbed by the diluted long probe chains. The simulation results are qualitatively consistent with the theoretical predictions based on constraint release Rouse model, but a detailed comparison reveals the existence of a broad distribution of the disentanglement rates, which is partly confirmed by an analysis of the packing and diffusion of the matrix chains in the tube region of the probe chains. A coarse-grained simulation model based on the tube Rouse motion model with incorporation of the probability distribution of the tube segment jump rates is developed and shows results qualitatively consistent with the fine scale molecular dynamics simulations. However, we observe a breakdown in the tube Rouse model when the short chain length is decreased to around N-S = 80, which is roughly 3.5 times the entanglement spacing N-e(P) = 23. The location of this transition may be sensitive to the chain bending potential used in our simulations.

Coarse‐grained force field for the nucleosome from self‐consistent multiscaling

A coarse-grained simulation model for the nucleosome is developed, using a methodology modified from previous work on the ribosome. Protein residues and DNA nucleotides are represented as beads, interacting through harmonic (for neighboring) or Morse (for nonbonded) potentials. Force-field parameters were estimated by Boltzmann inversion of the corresponding radial distribution functions obtained from a 5-ns all-atom molecular dynamics (MD) simulation, and were refined to produce agreement with the all-atom MD simulation. This self-consistent multiscale approach yields a coarse-grained model that is capable of reproducing equilibrium structural properties calculated from a 50-ns all-atom MD simulation. This coarse-grained model speeds up nucleosome simulations by a factor of 10(3) and is expected to be useful in examining biologically relevant dynamical nucleosome phenomena on the microsecond timescale and beyond.

Coarse-grained lattice model simulations of sequence-structure fitness of a ribosome-inactivating protein

Many realistic protein-engineering design problems extend beyond the computational limits of what is considered practical when applying all-atom molecular-dynamics simulation methods. Lattice models provide computationally robust alternatives, yet most are regarded as too simplistic to accurately capture the details of complex designs. We revisit a coarse-grained lattice simulation model and demonstrate that a multiresolution modeling approach of reconstructing all-atom structures from lattice chains is of sufficient accuracy to resolve the comparability of sequence-structure modifications of the ricin A-chain (RTA) protein fold. For a modeled structure, the unfolding-folding transition temperature was calculated from the heat capacity using either the potential energy from the lattice model or the all-atom CHARMM19 force-field plus a generalized Born solvent approximation. We found, that despite the low-resolution modeling of conformational states, the potential energy functions were capable of detecting the relative change in the thermodynamic transition temperature that distinguishes between a protein design and the native RTA fold in excellent accord with reported experimental studies of thermal denaturation. A discussion is provided of different sequences fitted to the RTA fold and a possible unfolding model.

Aggregation of POSS Monomers in Liquid Hexane: A Molecular-Simulation Study

Polyhedral oligomeric silsesquioxanes (POSS) are multifunctional molecules that can be employed as building blocks to develop nanocomposite materials whose mechanical properties often improve upon those of traditional polymeric materials. We report here molecular simulation results for the effective potential of mean force between octamethyl POSS monomers and between POSS monomers in which one methyl group has been substituted by a linear alkane chain of nine carbon atoms in liquid normal hexane at 300 and 400 K. The results are discussed and compared to available data for the effective interactions between octamethyl POSS monomers in normal hexadecane. Our results show that the effective short-ranged POSS-POSS attraction is significantly weaker in hexane than it is in hexadecane, perhaps explaining why normal hexane is often the solvent of choice for the preparation of POSS-containing materials. Additionally, we provide results for the radial distribution functions between selected sites in the POSS monomers that can be used both to understand the association between POSS monomers in solution and to parametrize coarse-grained simulation models. Such models will be used to study the formation of POSS-containing supramolecular structures such as lamellae or micelles that are currently not accessible by atomistic simulation and can be compared to experimental observations.

Modeling real dynamics in the coarse-grained representation of condensed phase systems

This work presents a systematic multiscale methodology to provide a more faithful representation of real dynamics in coarse-grained molecular simulation models. The theoretical formalism is based on the recently developed multiscale coarse-graining (MS-CG) method [S. Izvekov and G. A. Voth, J. Phys. Chem. B. 109, 2469 (2005); J. Chem. Phys. 123, 134105 (2005)] and relies on the generalized Langevin equation approach and its simpler Langevin equation limit. The friction coefficients are determined in multiscale fashion from the underlying all-atom molecular dynamics simulations using force-velocity and velocity-velocity correlation functions for the coarse-grained sites. The diffusion properties in the resulting CG Brownian dynamics simulations are shown to be quite accurate. The time dependence of the velocity autocorrelation function is also well-reproduced relative to the all-atom model if sufficient resolution of the CG sites is implemented.

A comparison of united atom, explicit atom, and coarse-grained simulation models for poly(ethylene oxide)

We compare static and dynamic properties obtained from three levels of modeling for molecular dynamics simulation of poly(ethylene oxide) (PEO). Neutron scattering data are used as a test of each model's accuracy. The three simulation models are an explicit atom (EA) model (all the hydrogens are taken into account explicitly), a united atom (UA) model (CH(2) and CH(3) groups are considered as a single unit), and a coarse-grained (CG) model (six united atoms are taken as one bead). All three models accurately describe the PEO static structure factor as measured by neutron diffraction. Dynamics are assessed by comparison to neutron time of flight data, which follow self-motion of protons. Hydrogen atom motion from the EA model and carbon/oxygen atom motion from the UA model closely follow the experimental hydrogen motion, while hydrogen atoms reinserted in the UA model are too fast. The EA and UA models provide a good description of the orientation properties of C-H vectors measured by nuclear magnetic resonance experiments. Although dynamic observables in the CG model are in excellent agreement with their united atom counterparts, they cannot be compared to neutron data because the time after which the CG model is valid is greater than the neutron decay times.

Molecular View of Hexagonal Phase Formation in Phospholipid Membranes

Important biological processes, such as vesicle fusion or budding, require the cell matrix to undergo a transition from a lamellar to a nonlamellar state. Although equilibrium properties of membranes are amenable to detailed theoretical studies, collective rearrangements involved in phase transitions have thus far only been modeled on a qualitative level. Here, for the first time, the complete transition pathway from a multilamellar to an inverted hexagonal phase is elucidated at near-atomic detail using a recently developed coarse-grained molecular dynamics simulation model. Insight is provided into experimentally inaccessible data such as the molecular structure of the intermediates and the kinetics involved. Starting from multilamellar configurations, the spontaneous formation of stalks between the bilayers is observed on a nanosecond timescale at elevated temperatures or reduced hydration levels. The stalks subsequently elongate in a cooperative manner leading to the formation of an inverted hexagonal phase. The rate of stalk elongation is ∼0.1 nm ns −1. Within a narrow hydration/temperature/composition range the stalks appear stable and rearrange into the rhombohedral phase.