Adaptive resolution simulation of an atomistic DNA molecule in MARTINI salt solution

We present a dual-resolution model of a deoxyribonucleic acid (DNA) molecule in a bathing solution, where we concurrently couple atomistic bundled water and ions with the coarse-grained MAR- TINI model of the solvent. We use our fine-grained salt solution model as a solvent in the inner shell surrounding the DNA molecule, whereas the solvent in the outer shell is modeled by the coarse-grained model. The solvent entities can exchange between the two domains and adapt their resolution accordingly. We critically asses the performance of our multiscale model in adaptive resolution simulations of an infinitely long DNA molecule, focusing on the structural characteristics of the solvent around DNA. Our analysis shows that the adaptive resolution scheme does not produce any noticeable artifacts in comparison to a refer- ence system simulated in full detail. The effect of using a bundled-SPC model, required for multiscaling, compared to the standard free SPC model is also evaluated. Our multiscale approach opens the way for large scale applications of DNA and other biomolecules which require a large solvent reservoir to avoid boundary effects

Polyhedral vesicles

Polyhedral vesicles with a large bending modulus of the membrane such as the gel phase lipid membrane were studied using a Brownian dynamics simulation. The vesicles exhibit various polyhedral morphologies such as tetrahedron and cube shapes. We clarified two types of line defects on the edges of the polyhedrons: cracks of both monolayers at the spontaneous curvature of monolayer $C_{\text {0}}<0$, and a crack of the inner monolayer at $C_{\text {0}}\ge0$. Around the latter defect, the inner monolayer curves positively. Our results suggested that the polyhedral morphology is controlled by $C_{\text {0}}$.Comment: 4 pages, 5 figure

Hydrophobic and ionic-interactions in bulk and confined water with implications for collapse and folding of proteins

Water and water-mediated interactions determine thermodynamic and kinetics of protein folding, protein aggregation and self-assembly in confined spaces. To obtain insights into the role of water in the context of folding problems, we describe computer simulations of a few related model systems. The dynamics of collapse of eicosane shows that upon expulsion of water the linear hydrocarbon chain adopts an ordered helical hairpin structure with 1.5 turns. The structure of dimer of eicosane molecules has two well ordered helical hairpins that are stacked perpendicular to each other. As a prelude to studying folding in confined spaces we used simulations to understand changes in hydrophobic and ionic interactions in nano droplets. Solvation of hydrophobic and charged species change drastically in nano water droplets. Hydrophobic species are localized at the boundary. The tendency of ions to be at the boundary where water density is low increases as the charge density decreases. Interaction between hydrophobic, polar, and charged residue are also profoundly altered in confined spaces. Using the results of computer simulations and accounting for loss of chain entropy upon confinement we argue and then demonstrate, using simulations in explicit water, that ordered states of generic amphiphilic peptide sequences should be stabilized in cylindrical nanopores

The Martini Coarse-Grained Force Field

The Martini force field is a coarse-grained force field suited for molecular dynamics simulations of biomolecular systems. The force field has been parameterized in a systematic way, based on the reproduction of partitioning free energies between polar and apolar phases of a large number of chemical compounds. In this chapter the methodology underlying the force field is presented together with details of its parameterization and limitations. Then currently available topologies are described with a short overview of the key elements of their parameterization. These include the new polarizable Martini water model. A set of three selected ongoing studies using the Martini force field is presented. Finally the latest lines of development are discussed

Effect of Undulations on Surface Tension in Simulated Bilayers

To understand the effect of the finite size of simulation cells on the equilibrium properties of bilayers, an extensive series of glycerolmonoolein bilayer molecular dynamics simulations in which the surface area and system size were systematically changed have been conducted. Systems ranging from 200 to 1800 lipids were simulated, covering length scales up to 20 nm. The dependence of the surface tension on the area per lipid is shown, although long simulation times were needed (up to 40 ns) to obtain reliable estimates. As the size of simulated patches increases, long wavelength undulatory modes appear with a concomitant increase in the area compressibility due to coupling of undulation modes to area fluctuations. Both the undulatory intensities and the peristaltic intensities of the bilayer fluctuations can be fitted in the long wavelength limit to continuum model predictions. The effect of system size on surface tension appears to depend on the stress conditions.

Molecular Dynamics Simulations of Mixed Micelles Modeling Human Bile

Extensive molecular dynamics (MD) simulations of binary systems of phospholipids and bile salts, a model for human bile, have been performed. Recent progress in hardware and software development allows simulation of the spontaneous aggregation of the constituents into small mixed micelles, in agreement with experimental observations. The MD simulations reveal the structure of these micelles at atomic detail. The phospholipids are packed radially with their headgroups at the surface and the hydrophobic tails pointing toward the micellar center. The bile salts act as wedges between the phospholipid headgroups, with their hydrophilic sides exposed to the aqueous environment. The structure of the micelles strongly resembles the previously proposed radial shell model. Simulations including small fractions of cholesterol reveal how cholesterol is solubilized inside these mixed micelles without changing their overall structure.