Backmapping triangulated surfaces to coarse-grained membrane models

Molecular Dynamics Simulations of Curved Lipid Membranes

Eukaryotic cells contain membranes with various curvatures, from the near-plane plasma membrane to the highly curved membranes of organelles, vesicles, and membrane protrusions. These curvatures are generated and sustained by curvature-inducing proteins, peptides, and lipids, and describing these mechanisms is an important scientific challenge. In addition to that, some molecules can sense membrane curvature and thereby be trafficked to specific locations. The description of curvature sensing is another fundamental challenge. Curved lipid membranes and their interplay with membrane-associated proteins can be investigated with molecular dynamics (MD) simulations. Various methods for simulating curved membranes with MD are discussed here, including tools for setting up simulation of vesicles and methods for sustaining membrane curvature. The latter are divided into methods that exploit scaffolding virtual beads, methods that use curvature-inducing molecules, and methods applying virtual forces. The variety of simulation tools allow researcher to closely match the conditions of experimental studies of membrane curvatures.

Asynchronous Reciprocal Coupling of Martini 2.2 Coarse-Grained and CHARMM36 All-Atom Simulations in an Automated Multiscale Framework

The appeal of multiscale modeling approaches is predicated on the promise of combinatorial synergy. However, this promise can only be realized when distinct scales are combined with reciprocal consistency. Here, we consider multiscale molecular dynamics (MD) simulations that combine the accuracy and macromolecular flexibility accessible to fixed-charge all-atom (AA) representations with the sampling speed accessible to reductive, coarse-grained (CG) representations. AA-to-CG conversions are relatively straightforward because deterministic routines with unique outcomes are achievable. Conversely, CG-to-AA conversions have many solutions due to a surge in the number of degrees of freedom. While automated tools for biomolecular CG-to-AA transformation exist, we find that one popular option, called Backward, is prone to stochastic failure and the AA models that it does generate frequently have compromised protein structure and incorrect stereochemistry. Although these shortcomings can likely be circumvented by human intervention in isolated instances, automated multiscale coupling requires reliable and robust scale conversion. Here, we detail an extension to Multiscale Machine-learned Modeling Infrastructure (MuMMI), including an improved CG-to-AA conversion tool called sinceCG. This tool is reliable (∼98% weakly correlated repeat success rate), automatable (no unrecoverable hangs), and yields AA models that generally preserve protein secondary structure and maintain correct stereochemistry. We describe how the MuMMI framework identifies CG system configurations of interest, converts them to AA representations, and simulates them at the AA scale while on-the-fly analyses provide feedback to update CG parameters. Application to systems containing the peripheral membrane protein RAS and proximal components of RAF kinase on complex eight-component lipid bilayers with ∼1.5 million atoms is discussed in the context of MuMMI.

Modelling structural properties of cyanine dye nanotubes at coarse-grained level

Self-assembly is a ubiquitous process spanning from biomolecular aggregates to nanomaterials. Even though the resulting aggregates can be studied through experimental techniques, the dynamic pathways of the process and the molecular details of the final structures are not necessarily easy to resolve. Consequently, rational design of self-assembling aggregates and their properties remains extremely challenging. At the same time, modelling the self-assembly with computational methods is not trivial, because its spatio-temporal scales are usually beyond the limits of all-atom based simulations. The use of coarse-grained (CG) models can alleviate this limitation, but usually suffers from the lack of optimised parameters for the molecular constituents. In this work, we describe the procedure of parametrizing a CG Martini model for a cyanine dye (C8S3) that self-assembles into hollow double-walled nanotubes. First, we optimised the model based on quantum mechanics calculations and all-atom reference simulations, in combination with available experimental data. Then, we conducted random self-assembly simulations, and the performance of our model was tested on preformed assemblies. Our simulations provide information on the time-dependent local arrangement of this cyanine dye, when aggregates are being formed. Furthermore, we provide guidelines for designing and optimising parameters for similar self-assembling nanomaterials.

Key aspects of the past 30 years of protein design

Proteins are the workhorse of life. They are the building infrastructure of living systems; they are the most efficient molecular machines known, and their enzymatic activity is still unmatched in versatility by any artificial system. Perhaps proteins' most remarkable feature is their modularity. The large amount of information required to specify each protein's function is analogically encoded with an alphabet of just ∼20 letters. The protein folding problem is how to encode all such information in a sequence of 20 letters. In this review, we go through the last 30 years of research to summarize the state of the art and highlight some applications related to fundamental problems of protein evolution.

Integrating cellular and molecular structures and dynamics into whole-cell models

A complete description of the state of the cell requires knowledge of its size, shape, components, intracellular reactions, and interactions with its environment-all of these as a function of time and cell growth. Adding to this list is the need for theoretical models and simulations that integrate and help to interpret this daunting amount of experimental data. It seems like an overwhelming list of requirements, but progress is being made on many fronts. In this review, we discuss the current challenges and problems in obtaining sufficient information about each aspect of a dynamical whole-cell model (DWCM) for simple and well-studied bacterial systems.

Specific interactions of peripheral membrane proteins with lipids: what can molecular simulations show us?

Peripheral membrane proteins (PMPs) can reversibly and specifically bind to biological membranes to carry out functions such as cell signalling, enzymatic activity, or membrane remodelling. Structures of these proteins and of their lipid-binding domains are typically solved in a soluble form, sometimes with a lipid or lipid headgroup at the binding site. To provide a detailed molecular view of PMP interactions with the membrane, computational methods such as molecular dynamics (MD) simulations can be applied. Here, we outline recent attempts to characterise these binding interactions, focusing on both intracellular proteins, such as phosphatidylinositol phosphate (PIP)-binding domains, and extracellular proteins such as glycolipid-binding bacterial exotoxins. We compare methods used to identify and analyse lipid-binding sites from simulation data and highlight recent work characterising the energetics of these interactions using free energy calculations. We describe how improvements in methodologies and computing power will help MD simulations to continue to contribute to this field in the future.

The ugly, bad, and good stories of large-scale biomolecular simulations

Molecular modeling of large biomolecular assemblies exemplifies a disruptive area holding both promises and contentions. Propelled by peta and exascale computing, several simulation methodologies have now matured into user-friendly tools that are successfully employed for modeling viruses, membranous nano-constructs, and key pieces of the genetic machinery. We present three unifying biophysical themes that emanate from some of the most recent multi-million atom simulation endeavors. Despite connecting molecular changes with phenotypic outcomes, the quality measures of these simulations remain questionable. We discuss the existing and upcoming strategies for constructing representative ensembles of large systems, how new computing technologies will boost this area, and make a point that integrative modeling guided by experimental data is the future of biomolecular computations.

Assembly and Analysis of Cell-Scale Membrane Envelopes

The march toward exascale computing will enable routine molecular simulation of larger and more complex systems, for example, simulation of entire viral particles, on the scale of approximately billions of atoms─a simulation size commensurate with a small bacterial cell. Anticipating the future hardware capabilities that will enable this type of research and paralleling advances in experimental structural biology, efforts are currently underway to develop software tools, procedures, and workflows for constructing cell-scale structures. Herein, we describe our efforts in developing and implementing an efficient and robust workflow for construction of cell-scale membrane envelopes and embedding membrane proteins into them. A new approach for construction of massive membrane structures that are stable during the simulations is built on implementing a subtractive assembly technique coupled with the development of a structure concatenation tool (fastmerge), which eliminates overlapping elements based on volumetric criteria rather than adding successive molecules to the simulation system. Using this approach, we have constructed two "protocells" consisting of MARTINI coarse-grained beads to represent cellular membranes, one the size of a cellular organelle and another the size of a small bacterial cell. The membrane envelopes constructed here remain whole during the molecular dynamics simulations performed and exhibit water flux only through specific proteins, demonstrating the success of our methodology in creating tight cell-like membrane compartments. Extended simulations of these cell-scale structures highlight the propensity for nonspecific interactions between adjacent membrane proteins leading to the formation of protein microclusters on the cell surface, an insight uniquely enabled by the scale of the simulations. We anticipate that the experiences and best practices presented here will form the basis for the next generation of cell-scale models, which will begin to address the addition of soluble proteins, nucleic acids, and small molecules essential to the function of a cell.

Polyply; a python suite for facilitating simulations of macromolecules and nanomaterials

Molecular dynamics simulations play an increasingly important role in the rational design of (nano)-materials and in the study of biomacromolecules. However, generating input files and realistic starting coordinates for these simulations is a major bottleneck, especially for high throughput protocols and for complex multi-component systems. To eliminate this bottleneck, we present the polyply software suite that provides 1) a multi-scale graph matching algorithm designed to generate parameters quickly and for arbitrarily complex polymeric topologies, and 2) a generic multi-scale random walk protocol capable of setting up complex systems efficiently and independent of the target force-field or model resolution. We benchmark quality and performance of the approach by creating realistic coordinates for polymer melt simulations, single-stranded as well as circular single-stranded DNA. We further demonstrate the power of our approach by setting up a microphase-separated block copolymer system, and by generating a liquid-liquid phase separated system inside a lipid vesicle.

Coarse-Grain Simulations of Membrane-Adsorbed Helical Peptides

The amphipathic α-helix is a common motif for peptide adsorption to membranes. Many physiologically relevant events involving membrane-adsorbed peptides occur over time and size scales readily accessible to coarse-grain molecular dynamics simulations. This methodological suitability, however, comes with a number of pitfalls. Here, I exemplify a multi-step adsorption equilibration procedure on the antimicrobial peptide Magainin 2. It involves careful control of peptide freedom to promote optimal membrane adsorption before other interactions are allowed. This shortens preparation times prior to production simulations while avoiding divergence into unrealistic or artifactual configurations.

Coupling Bulk Phase Separation of Disordered Proteins to Membrane Domain Formation in Molecular Simulations on a Bespoke Compute Fabric

Phospholipid membranes surround the cell and its internal organelles, and their multicomponent nature allows the formation of domains that are important in cellular signalling, the immune system, and bacterial infection. Cytoplasmic compartments are also created by the phase separation of intrinsically disordered proteins into biomolecular condensates. The ubiquity of lipid membranes and protein condensates raises the question of how three-dimensional droplets might interact with two-dimensional domains, and whether this coupling has physiological or pathological importance. Here, we explore the equilibrium morphologies of a dilute phase of a model disordered protein interacting with an ideal-mixing, two-component lipid membrane using coarse-grained molecular simulations. We find that the proteins can wet the membrane with and without domain formation, and form phase separated droplets bound to membrane domains. Results from much larger simulations performed on a novel non-von-Neumann compute architecture called POETS, which greatly accelerates their execution compared to conventional hardware, confirm the observations. Reducing the wall clock time for such simulations requires new architectures and computational techniques. We demonstrate here an inter-disciplinary approach that uses real-world biophysical questions to drive the development of new computing hardware and simulation algorithms.

Sequential Voxel-Based Leaflet Segmentation of Complex Lipid Morphologies

As molecular dynamics simulations increase in complexity, new analysis tools are necessary to facilitate interpreting the results. Lipids, for instance, are known to form many complicated morphologies, because of their amphipathic nature, becoming more intricate as the particle count increases. A few lipids might form a micelle, where aggregation of tens of thousands could lead to vesicle formation. Millions of lipids comprise a cell and its organelle membranes, and are involved in processes such as neurotransmission and transfection. To study such phenomena, it is useful to have analysis tools that understand what is meant by emerging entities such as micelles and vesicles. Studying such systems at the particle level only becomes extremely tedious, counterintuitive, and computationally expensive. To address this issue, we developed a method to track all the individual lipid leaflets, allowing for easy and quick detection of topological changes at the mesoscale. By using a voxel-based approach and focusing on locality, we forego costly geometrical operations without losing important details and chronologically identify the lipid segments using the Jaccard index. Thus, we achieve a consistent sequential segmentation on a wide variety of (lipid) systems, including monolayers, bilayers, vesicles, inverted hexagonal phases, up to the membranes of a full mitochondrion. It also discriminates between adhesion and fusion of leaflets. We show that our method produces consistent results without the need for prefitting parameters, and segmentation of millions of particles can be achieved on a desktop machine.

Accelerated Simulations of Molecular Systems through Learning of Effective Dynamics

Simulations are vital for understanding and predicting the evolution of complex molecular systems. However, despite advances in algorithms and special purpose hardware, accessing the time scales necessary to capture the structural evolution of biomolecules remains a daunting task. In this work, we present a novel framework to advance simulation time scales by up to 3 orders of magnitude by learning the effective dynamics (LED) of molecular systems. LED augments the equation-free methodology by employing a probabilistic mapping between coarse and fine scales using mixture density network (MDN) autoencoders and evolves the non-Markovian latent dynamics using long short-term memory MDNs. We demonstrate the effectiveness of LED in the Müller-Brown potential, the Trp cage protein, and the alanine dipeptide. LED identifies explainable reduced-order representations, i.e., collective variables, and can generate, at any instant, all-atom molecular trajectories consistent with the collective variables. We believe that the proposed framework provides a dramatic increase to simulation capabilities and opens new horizons for the effective modeling of complex molecular systems.

Structural Insights into the Human Mitochondrial Pyruvate Carrier Complexes

Pyruvate metabolism requires the mitochondrial pyruvate carrier (MPC) proteins to transport pyruvate from the intermembrane space through the inner mitochondrial membrane to the mitochondrial matrix. The lack of the atomic structures of MPC hampers the understanding of the functional states of MPC and molecular interactions with the substrate or inhibitor. Here, we develop the de novo models of human MPC complexes and characterize the conformational dynamics of the MPC heterodimer formed by MPC1 and MPC2 (MPC1/2) by computational simulations. Our results reveal that functional MPC1/2 prefers to adopt an inward-open conformation, with the carrier open to the matrix side, whereas the outward-open states are less populated. The energy barrier for pyruvate transport in MPC1/2 is low enough, and the inhibitor UK5099 blocks the pyruvate transport by stably binding to MPC1/2. Notably, consistent with experimental results, the MPC1 L79H mutation significantly alters the conformations of MPC1/2 and thus fails for substrate transport. However, the MPC1 R97W mutation seems to retain the transport activity. The present de novo models of MPC complexes provide structural insights into the conformational states of MPC complexes and mechanistic understanding of interactions between the substrate/inhibitor and MPC proteins.

CG2AT2: an Enhanced Fragment-Based Approach for Serial Multi-scale Molecular Dynamics Simulations

Coarse-grained molecular dynamics provides a means for simulating the assembly and interactions of macromolecular complexes at a reduced level of representation, thereby allowing both longer timescale and larger sized simulations. Here, we describe an enhanced fragment-based protocol for converting macromolecular complexes from coarse-grained to atomistic resolution, for further refinement and analysis. While the focus is upon systems that comprise an integral membrane protein embedded in a phospholipid bilayer, the technique is also suitable for membrane-anchored and soluble protein/nucleotide complexes. Overall, this provides a method for generating an accurate and well-equilibrated atomic-level description of a macromolecular complex. The approach is evaluated using a diverse test set of 11 system configurations of varying size and complexity. Simulations are assessed in terms of protein stereochemistry, conformational drift, lipid/protein interactions, and lipid dynamics.

Computational Approaches to Explore Bacterial Toxin Entry into the Host Cell

Many bacteria secrete toxic protein complexes that modify and disrupt essential processes in the infected cell that can lead to cell death. To conduct their action, these toxins often need to cross the cell membrane and reach a specific substrate inside the cell. The investigation of these protein complexes is essential not only for understanding their biological functions but also for the rational design of targeted drug delivery vehicles that must navigate across the cell membrane to deliver their therapeutic payload. Despite the immense advances in experimental techniques, the investigations of the toxin entry mechanism have remained challenging. Computer simulations are robust complementary tools that allow for the exploration of biological processes in exceptional detail. In this review, we first highlight the strength of computational methods, with a special focus on all-atom molecular dynamics, coarse-grained, and mesoscopic models, for exploring different stages of the toxin protein entry mechanism. We then summarize recent developments that are significantly advancing our understanding, notably of the glycolipid–lectin (GL-Lect) endocytosis of bacterial Shiga and cholera toxins. The methods discussed here are also applicable to the design of membrane-penetrating nanoparticles and the study of the phenomenon of protein phase separation at the surface of the membrane. Finally, we discuss other likely routes for future development.

Development of accurate coarse-grained force fields for weakly polar groups by an indirect parameterization strategy

Coarse-grained (CG) molecular dynamics simulations are widely used to predict morphological structures and interpret mechanisms of mesoscopic behavior between the scope of traditional experiments and all-atom simulations. However, most current CG force fields (FFs) are not precise enough, especially for polar molecules or functional groups. A main obstacle in developing accurate CG FFs for polar molecules is the freezing problem met at room temperature. In this work, we introduce an indirect parametrization strategy for weakly polar groups by considering their short-chain homologs to avoid freezing. Here, a polar group containing three to four heavy atoms is mapped into one CG bead that is connected to one alkyl bead composed of three or four carbons. The CG beads interact via 4-parameter nonbonded Morse potentials and harmonic bonded potentials. An efficient meta-multilinear interpolation parameterization algorithm, as recently developed by us, is used to rigorously optimize the force parameters. Satisfactory accuracy is witnessed in terms of the density, heat of vaporization, surface tension, and solvation free energy of the homologs of twelve polar molecules, all deviating from the experiment by less than 5%. The transferability of the current FF is indicated by the predicted density, heat of vaporization, and end-to-end distance distributions of fatty acid methyl esters composed of multiple functional groups parameterized in this work.

Molecular Dynamics for Antimicrobial Peptide Discovery

Although antimicrobial resistance is an increasingly significant public health concern, there have only been two new classes of antibiotics approved for human use since the 1960s. Understanding the mechanisms of action of antibiotics is critical for novel antibiotic discovery, but novel approaches are needed that do not exclusively rely on experiments. Molecular dynamics simulation is a computational tool that uses simple models of the atoms in a system to discover nanoscale insights into the dynamic relationship between mechanism and biological function. Such insights can lay the framework for elucidating the mechanism of action and optimizing antibiotic templates. Antimicrobial peptides represent a promising solution to escalating antimicrobial resistance, given their lesser tendency to induce resistance than that of small-molecule antibiotics. Simulations of these agents have already revealed how they interact with bacterial membranes and the underlying physiochemical features directing their structure and function. In this minireview, we discuss how traditional molecular dynamics simulation works and its role and potential for the development of new antibiotic candidates with an emphasis on antimicrobial peptides.

Simulating realistic membrane shapes

Biological membranes exhibit diversity in their shapes and complexity in chemical compositions that are linked to many cellular functions. These two central features of biomembranes have been the subject of numerous simulation studies, using a diverse range of computational techniques. Currently, the field is able to capture this complexity at increasing levels of realism and connect the microscopic view on protein-lipid interactions to cellular morphologies at the level of entire organelles. Here we highlight recent advances in this topic, identify current bottlenecks, and sketch possible ways ahead.

Not only in silico drug discovery: Molecular modeling towards in silico drug delivery formulations

The use of methods at molecular scale for the discovery of new potential active ligands, as well as previously unknown binding sites for target proteins, is now an established reality. Literature offers many successful stories of active compounds developed starting from insights obtained in silico and approved by Food and Drug Administration (FDA). One of the most famous examples is raltegravir, a HIV integrase inhibitor, which was developed after the discovery of a previously unknown transient binding area thanks to molecular dynamics simulations. Molecular simulations have the potential to also improve the design and engineering of drug delivery devices, which are still largely based on fundamental conservation equations. Although they can highlight the dominant release mechanism and quantitatively link the release rate to design parameters (size, drug loading, et cetera), their spatial resolution does not allow to fully capture how phenomena at molecular scale influence system behavior. In this scenario, the "computational microscope" offered by simulations at atomic scale can shed light on the impact of molecular interactions on crucial parameters such as release rate and the response of the drug delivery device to external stimuli, providing insights that are difficult or impossible to obtain experimentally. Moreover, the new paradigm brought by nanomedicine further underlined the importance of such computational microscope to study the interactions between nanoparticles and biological components with an unprecedented level of detail. Such knowledge is a fundamental pillar to perform device engineering and to achieve efficient and safe formulations. After a brief theoretical background, this review aims at discussing the potential of molecular simulations for the rational design of drug delivery systems.

Optimizing Gō-MARTINI Coarse-Grained Model for F-BAR Protein on Lipid Membrane

Coarse-grained (CG) molecular dynamics (MD) simulations allow us to access much larger length and time scales than atomistic MD simulations, providing an attractive alternative to the conventional simulations. Based on the well-known MARTINI CG force field, the recently developed Gō-MARTINI model for proteins describes large-amplitude structural dynamics, which has not been possible with the commonly used elastic network model. Using the Gō-MARTINI model, we conduct MD simulations of the F-BAR Pacsin1 protein on lipid membrane. We observe that structural changes of the non-globular protein are largely dependent on the definition of the native contacts in the Gō model. To address this issue, we introduced a simple cutoff scheme and tuned the cutoff distance of the native contacts and the interaction strength of the Lennard-Jones potentials in the Gō-MARTINI model. With the optimized Gō-MARTINI model, we show that it reproduces structural fluctuations of the Pacsin1 dimer from atomistic simulations. We also show that two Pacsin1 dimers properly assemble through lateral interaction on the lipid membrane. Our work presents a first step towards describing membrane remodeling processes in the Gō-MARTINI CG framework by simulating a crucial step of protein assembly on the membrane.

Power Failure of Mitochondria and Oxidative Stress in Neurodegeneration and Its Computational Models

The brain needs more energy than other organs in the body. Mitochondria are the generator of vital power in the living organism. Not only do mitochondria sense signals from the outside of a cell, but they also orchestrate the cascade of subcellular events by supplying adenosine-5'-triphosphate (ATP), the biochemical energy. It is known that impaired mitochondrial function and oxidative stress contribute or lead to neuronal damage and degeneration of the brain. This mini-review focuses on addressing how mitochondrial dysfunction and oxidative stress are associated with the pathogenesis of neurodegenerative disorders including Alzheimer's disease, amyotrophic lateral sclerosis, Huntington's disease, and Parkinson's disease. In addition, we discuss state-of-the-art computational models of mitochondrial functions in relation to oxidative stress and neurodegeneration. Together, a better understanding of brain disease-specific mitochondrial dysfunction and oxidative stress can pave the way to developing antioxidant therapeutic strategies to ameliorate neuronal activity and prevent neurodegeneration.

Wrapping Up Viruses at Multiscale Resolution: Optimizing PACKMOL and SIRAH Execution for Simulating the Zika Virus

Simulating huge biomolecular complexes of million atoms at relevant biological time scales is becoming accessible to the broad scientific community. That proves to be crucial for urgent responses against emergent diseases in real time. Yet, there are still issues to sort regarding the system setup so that molecular dynamics (MD) simulations can be run in a simple and standard way. Here, we introduce an optimized pipeline for building and simulating enveloped virus-like particles (VLP). First, the membrane packing problem is tackled with new features and optimized options in PACKMOL. This allows preparing accurate membrane models of thousands of lipids in the context of a VLP within a few hours using a single CPU. Then, the assembly of the VLP system is done within the multiscale framework of the coarse-grained SIRAH force field. Finally, the equilibration protocol provides a system ready for production MD simulations within a few days on broadly accessible GPU resources. The pipeline is applied to study the Zika virus as a test case for large biomolecular systems. The VLP stabilizes at approximately 0.5 μs of MD simulation, reproducing correlations greater than 0.90 against experimental density maps from cryo-electron microscopy. Detailed structural analysis of the protein envelope also shows very good agreement in root-mean-square deviations and B-factors with the experimental data. The level of details attained shows for the first time a possible role for anionic phospholipids in stabilizing the envelope. Combining an efficient and reliable setup procedure with an accurate coarse-grained force field provides a valuable pipeline for simulating arbitrary viral systems or subcellular compartments, paving the way toward whole-cell simulations.

Computational and Experimental Approaches to Investigate Lipid Nanoparticles as Drug and Gene Delivery Systems

Lipid nanoparticles (LNPs) have been widely applied in drug and gene delivery. More than twenty years ago, DoxilTM was the first LNPs-based drug approved by the US Food and Drug Administration (FDA). Since then, with decades of research and development, more and more LNP-based therapeutics have been used to treat diverse diseases, which often offer the benefits of reduced toxicity and/or enhanced efficacy compared to the active ingredients alone. Here, we provide a review of recent advances in the development of efficient and robust LNPs for drug/gene delivery. We emphasize the importance of rationally combining experimental and computational approaches, especially those providing multiscale structural and functional information of LNPs, to the design of novel and powerful LNP-based delivery systems.

Double membrane formation in heterogeneous vesicles

Lipids are capable of forming a variety of structures, including multi-lamellar vesicles. Layered lipid membranes are found in cell organelles, such as autophagosomes and mitochondria. Here, we present a mechanism for the formation of a double-walled vesicle (i.e., two lipid bilayers) from a unilamellar vesicle through the partitioning and phase separation of a small molecule. Using molecular dynamics simulations, we show that double membrane formation proceeds via a nucleation and growth process - i.e., after a critical concentration of the small molecules, a patch of double membrane nucleates and grows to cover the entire vesicle. We discuss the implications of this mechanism and theoretical approaches for understanding the evolution and formation of double membranes.

Extending the Martini Coarse-Grained Force Field to N-Glycans

Glycans play a vital role in a large number of cellular processes. Their complex and flexible nature hampers structure-function studies using experimental techniques. Molecular dynamics (MD) simulations can help in understanding dynamic aspects of glycans if the force field parameters used can reproduce key experimentally observed properties. Here, we present optimized coarse-grained (CG) Martini force field parameters for N-glycans, calibrated against experimentally derived binding affinities for lectins. The CG bonded parameters were obtained from atomistic (ATM) simulations for different glycan topologies including high mannose and complex glycans with various branching patterns. In the CG model, additional elastic networks are shown to improve maintenance of the overall conformational distribution. Solvation free energies and octanol-water partition coefficients were also calculated for various N-glycan disaccharide combinations. When using standard Martini nonbonded parameters, we observed that glycans spontaneously aggregated in the solution and required down-scaling of their interactions for reproduction of ATM model radial distribution functions. We also optimized the nonbonded interactions for glycans interacting with seven lectin candidates and show that a relatively modest scaling down of the glycan-protein interactions can reproduce free energies obtained from experimental studies. These parameters should be of use in studying the role of glycans in various glycoproteins and carbohydrate binding proteins as well as their complexes, while benefiting from the efficiency of CG sampling.