Accelerating cellulose nanocomposite design through coarse-grained molecular dynamics simulation and rule of mixture analysis

Multi-scale Ensemble Modeling of Modular Proteins with Intrinsically Disordered Linker Regions: Application to p53

Coarse-grained (CG) molecular dynamics (MD) simulations are performed to study the insertion of cyclic peptide nanotubes into cell membranes and to examine whether cyclic peptide nanotubes can function as an ion channel and thereby as an antibacterial agent. To do so, the two coarse-grained (CG) models for lipid molecules and for proteins developed by Marrink et al. (J. Phys. Chem. B 2004, 108, 750) and by Shih et al. (J. Phys. Chem. B 2006, 110, 3674), respectively, were extended and modified. These CG models were verified by performing CG MD and all-atom (AA) MD simulations for a cyclic peptide nanotube, 8 x cyclo[(-d-Ala-l-Glu-d-Ala-l-Gln-)2], in water and by comparing the results from the two simulations. Comparison between static and dynamic (water transport) properties obtained from both simulations shows good agreement. To study nanotube insertion, a CG cyclic peptide nanotube, 8 x cyclo[(-Trp-d-Leu-)4], was prepared above the surface of a CG DPPC lipid bilayer, restrained with constraints, and equilibrated, and then a series of CG MD simulations were carried out by lifting the constraints imposed on the nanotube. The CG MD simulations show that the cyclic peptide nanotube spontaneously inserts into and reorients inside the lipid bilayer. After insertion, the long axis of the cyclic peptide nanotube is aligned approximately perpendicular to the bilayer plane indicating that the nanotube can function as an ion channel and as an antibacterial agent. Tilt structures of the cyclic peptide nanotubes inside the lipid bilayer are found to be in agreement with experiment and earlier AA simulations. Lipid flip-flop, a migration of lipid molecules from one leaflet to the other leaflet of the lipid bilayer, is also observed from the CG MD simulations. Finally, the CG MD simulations reveal that a lipid headgroup can be inserted into the cyclic peptide nanotube. This process is confirmed by an AA MD simulation.

<p>The use of computational simulation to study the dynamics and interactions of macromolecules has become an important tool in the field of biochemistry. A common method to perform these simulations is to use all-atom explicit-solvent molecular dynamics (MD). However, due to the limitations in computational power currently available, this method is not practical for simulating large-scale biomolecular systems on long timescales. An alternative is to perform implicit-solvent Brownian dynamics (BD) simulations using a coarse grained (CG) model that allows for increased computational efficiency. However, if simulations using the CG model are not realistic, then the gain in computational efficiency from using a CG model is not worthwhile.</p> <p>This thesis describes the derivation of a set of bonded and nonbonded CG potential functions for use in implicit-solvent BD simulations of proteins derived from all-atom explicit-solvent MD simulations of amino acids. To determine which force field and water model to use in the MD simulations, Chapter II describes 1 Μs all-atom explicit-solvent MD simulations of glycine, asparagine, phenylalanine, and valine solutions at 50, 100, 200 and 300 mg/ml concentrations performed using eight different force field and water model combinations. To evaluate the accuracy of the force fields at high solute concentrations, the density, viscosity, and dielectric increments of the four amino acids were calculated from the simulations and compared to experimental results. Additionally, the change in the strength of hydrophobic and electrostatic interactions with increasing solute concentration was calculated for each force field and water model combination. As a result of this study, the Amber ff99SB-ILDN force field and TIP4P-Ew explicit-solvent water model were chosen for all subsequent MD simulations. Chapter III describes the derivation of CG bonded potential functions from 1 Μs all-atom explicit-solvent MD simulations of each of the twenty amino acids, including a separate simulation for protonated histidine. The angle and dihedral probability distributions sampled during the MD simulations were used to optimize the bonded potential functions using the iterative Boltzmann inversion (IBI) method. Chapter IV describes the derivation of CG nonbonded potential functions from 1 Μs all-atom explicit-solvent MD simulations of every possible pairing of the amino acids (231 different systems). The radial distribution functions calculated from these MD simulations were used to optimize a set of nonbonded CG potential functions using the IBI method. The optimized set of bonded and nonbonded potential functions, which is termed COFFDROP (COarse-grained Force Field for Dynamic Representation Of Proteins), quantitatively reproduced all of the calculated MD distributions. To determine if COFFDROP would be useful for simulations of bimolecular systems, Chapter V describes the testing of the transferability of the force field. First, COFFDROP was used to simulate concentrated amino acid solutions. The clustering of the solutes in these simulations was directly compared with results from corresponding all-atom explicit-solvent MD simulations and found to be in excellent agreement. Next, BD simulations of 9.2 mM solutions of the small protein villin headpiece were performed. The proteins aggregated during these simulations, which is in agreement with results from MD simulation but in disagreement with experiment. After scaling the strength of COFFDROP's nonbonded potential functions by a factor of 0.8 and rerunning the BD simulations, the amount of aggregation was comparable to experimental observations. Based on these results, COFFDROP is likely to be applicable in CG BD simulations of large, highly concentrated, biomolecular systems.</p>

Structure-based coarse-grained Gō-like models have been used extensively in deciphering protein folding mechanisms because of their simplicity and tractability. Meanwhile, explicit-solvent molecular dynamics (MD) simulations with physics-based all-atom force fields have been applied successfully to simulate folding/unfolding transitions for several small, fast-folding proteins. To explore the degree to which coarse-grained Gō-like models and their extensions to incorporate nonnative interactions are capable of producing folding processes similar to those in all-atom MD simulations, here we systematically compare the computed unfolded states, transition states, and transition paths obtained using coarse-grained models and all-atom explicit-solvent MD simulations. The conformations in the unfolded state in common Gō models are more extended, and are thus more in line with experiment, than those from all-atom MD simulations. Nevertheless, the structural features of transition states obtained by the two types of models are largely similar. In contrast, the folding transition paths are significantly more sensitive to modeling details. In particular, when common Gō-like models are augmented with nonnative interactions, the predicted dimensions of the unfolded conformations become similar to those computed using all-atom MD. With this connection, the large deviations of all-atom MD from simple diffusion theory are likely caused in part by the presence of significant nonnative effects in folding processes modelled by current atomic force fields. The ramifications of our findings to the application of coarse-grained modeling to more complex biomolecular systems are discussed.

Dynamics of Nucleosome Tails Studied by All-Atom and Coarse-Grained MD Simulations

Mechanics of Cadherin Unbinding using Coarse-Grained Models

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 Modeling of the Actin Filament Derived from Atomistic-Scale Simulations

Dynamic Coarse-Graining of DNA Melting Kinetics: Enthalpic and Entropic Effects of Cooperative Base-Pair Dynamics

Modeling the Open-to-closed Transition of Adenylate Kinase: All-atom Molecular Dynamics Simulations and a Double-Well Network Model

Bacterial cytoskeletal protein FtsZ binds and hydrolyzes GTP, and assembles into dynamic filaments that are essential for cell division. Here, we used a multi-scale computational strategy that combined all-atom molecular dynamics (MD) simulations and coarse-grained models to reveal the conformational dynamics of assembled FtsZ. We found that the top end of a filament is highly dynamic and can undergo T-to-R transitions in both GTP- and GDP-bound states. We observed several subcategories of nucleation related dimer species, which leading to a feasible nucleation pathway. In addition, we observed that FtsZ filament exhibits noticeable amounts of twisting, indicating a substantial helicity of the FtsZ filament. These results agree with the previously models and experimental data. Anisotropy network model (ANM) analysis revealed a polymerization enhanced assembly cooperativity, and indicated that the cooperative motions in FtsZ are encoded in the structure. Taken together, our study provides a molecular-level understanding of the diversity of the structural states of FtsZ and the relationships among polymerization, hydrolysis, and cooperative assembly, which should shed new light on the molecular basis of FtsZ’s cooperativity.

Amphipathic polymers called amphipols (APols) have been developed as an alternative to detergents for stabilizing membrane proteins (MPs) in aqueous solutions. APols provide MPs with a particularly mild environment and, as a rule, keep them in a native functional state for longer periods than do detergents. Amphipol A8-35, a derivative of polyacrylate, is widely used and has been particularly well studied experimentally. In aqueous solutions, A8-35 molecules self-assemble into well-defined globular particles with a mass of ∼40 kDa and a R(g) of ∼2.4 nm. As a first step towards describing MP/A8-35 complexes by molecular dynamics (MD), we present three sets of simulations of the pure APol particle. First, we performed a series of all-atom MD (AAMD) simulations of the particle in solution, starting from an arbitrary initial configuration. Although AAMD simulations result in stable cohesive particles over a 45 ns simulation, the equilibration of the particle organization is limited. This motivated the use of coarse-grained MD (CGMD), allowing us to investigate processes on the microsecond time scale, including de novo particle assembly. We present a detailed description of the parametrization of the CGMD model from the AAMD simulations and a characterization of the resulting CGMD particles. Our third set of simulations utilizes reverse coarse-graining (rCG), through which we obtain all-atom coordinates from a CGMD simulation. This allows a higher-resolution characterization of a configuration determined by a long-timescale simulation. Excellent agreement is observed between MD models and experimental, small-angle neutron scattering data. The MD data provides new insight into the structure and dynamics of A8-35 particles, which is possibly relevant to the stabilizing effects of APols on MPs, as well as a starting point for modeling MP/A8-35 complexes.

In this paper, we propose coarse-grained single-site (CGSS), wall-CO(2), and CO(2)-CO(2) interaction potential models to study the structure of carbon dioxide under confinement. The CGSS potentials are used in an empirical potential based quasi-continuum theory, EQT, to compute the center-of-mass density and potential profiles of CO(2) confined inside different size graphite slit pores. Results obtained from EQT are compared with those obtained from all-atom molecular dynamics (AA-MD) simulations, and are found to be in good agreement with each other. Though these CGSS interaction potentials are primarily developed and parameterized for EQT, they are also used to perform coarse-grained molecular dynamics (CG-MD) simulations. The results obtained from CG-MD simulations are also found to be in reasonable agreement with AA-MD simulation results.

Cell surface calreticulin (CRT) binding to thrombospondin-1 (TSP1), regulates cell adhesion, migration, anoikis resistance, and collagen production. Due to the essential role of membrane microdomains in CRT-mediated focal adhesion disassembly, we previously studied the effect of raft-like bilayers on TSP1–CRT interactions with all-atom molecular dynamics (AAMD) simulations. However, the simulated systems of protein on the surface of the bilayer(s) in the explicit solvent are too large for long timescale AAMD simulations due to computational expense. In this study, we adopted a multiscale modeling approach of combining AAMD, coarse-grained molecule dynamics (CGMD), and reversed AAMD (REV AAMD) simulations to investigate the interactions of single CRT or of the TSP1–CRT complex with a membrane microdomain at microsecond timescale. Results showed that CRT conformational stabilization by binding of TSP1 in AAMD simulation was undetectable in CGMD simulation, but it was recovered in REV AAMD simulation. Similarly, interactions of the CRT N-domain and TSP1 with the membrane microdomain were lost in CGMD simulations but they were re-gained in the REV AAMD simulations. There was the higher coordination of the CRT P-domain in the TSP1–CRT complex with the lipid components of membrane microdomain compared to that of single CRT, which could directly affect the conformation of CRT and further mediate CRT recruitment of LDL receptor-related protein for signaling events. This study provides structural and molecular insights into TSP1–CRT interactions in a membrane microdomain environment and demonstrates the feasibility of using multiscale simulations to investigate the interactions between protein and membrane microdomains at a long timescale.

RNA folding prediction is very meaningful and challenging. The molecular dynamics simulation (MDS) of all atoms (AA) is limited to the folding of small RNA molecules. At present, most of the practical models are coarse grained (CG) model, and the coarse-grained force field (CGFF) parameters usually depend on known RNA structures. However, the limitation of the CGFF is obvious that it is difficult to study the modified RNA. Based on the 3 beads model (AIMS_RNA_B3), we proposed the AIMS_RNA_B5 model with three beads representing a base and two beads representing the main chain (sugar group and phosphate group). We first run the all atom molecular dynamic simulation (AAMDS), and fit the CGFF parameter with the AA trajectory. Then perform the coarse-grained molecular dynamic simulation (CGMDS). AAMDS is the foundation of CGMDS. CGMDS is mainly to carry out the conformation sampling based on the current AAMDS state and improve the folding speed. We simulated the folding of three RNAs, which belong to hairpin, pseudoknot and tRNA respectively. Compared to the AIMS_RNA_B3 model, the AIMS_RNA_B5 model is more reasonable and performs better.