Search

Martini 3 is a widely used coarse-grained simulation method for large-scale biomolecular simulations. It can be combined with a Go̅ model to realistically describe higher-order protein structures while allowing the folding and unfolding events. However, as of today, this method has largely been used only for individual monomers. In this article, we describe how the Go̅ model can be implemented within the framework of Martini 3 for a multimer system, taking into account both intramolecular and intermolecular interactions in an oligomeric protein system. We demonstrate the method by showing how it can be applied to both structural stability maintenance and assembly/disassembly of protein oligomers, using aquaporin tetramer, insulin dimer, and amyloid-β fibril as examples. We find that addition of intermolecular Go̅ potentials stabilizes the quaternary structure of proteins. The strength of the Go̅ potentials can be tuned so that the internal fluctuations of proteins match the behavior of atomistic simulation models, however, the results also show that the use of too strong intermolecular Go̅ potentials weakens the chemical specificity of oligomerization. The Martini-Go̅ model presented here enables the use of Go̅ potentials in oligomeric molecular systems in a computationally efficient and parallelizable manner, especially in the case of homopolymers, where the number of identical protein monomers is high. This paves the way for coarse-grained simulations of large protein complexes, such as viral protein capsids and prion fibrils, in complex biological environments.

In molecular simulations, periodic boundary conditions are typically used to avoid surface effects occurring at the boundaries of the simulation box. A consequence of this is that molecules and assemblies may appear split over the boundaries. Broken molecular assemblies make it difficult to interpret, analyze, and visualize molecular simulation data. We present a general and fast algorithm that repairs molecular assemblies that are broken due to periodic boundary conditions. The open source method presented here, MDVWhole, works for all translation-only crystallographic periodic boundary conditions. The method consumes little memory and can fix the visualization of the assembly of millions of particles in a few seconds. Thus, it is suitable for processing both single simulation frames and long trajectories with millions of points.

Shape defines the structure and function of cellular membranes. In cell division, the cell membrane deforms into a "dumbbell" shape, while organelles such as the autophagosome exhibit "stomatocyte" shapes. Bottom-up in vitro reconstitution of protein machineries that stabilize or resolve the membrane necks in such deformed liposome structures is of considerable interest to characterize their function. Here we develop a DNA-nanotechnology-based approach that we call the synthetic membrane shaper (SMS), where cholesterol-linked DNA structures attach to the liposome membrane to reproducibly generate high yields of stomatocytes and dumbbells. In silico simulations confirm the shape-stabilizing role of the SMS. We show that the SMS is fully compatible with protein reconstitution by assembling bacterial divisome proteins (DynaminA, FtsZ:ZipA) at the catenoidal neck of these membrane structures. The SMS approach provides a general tool for studying protein binding to complex membrane geometries that will greatly benefit synthetic cell research.

Lipid nanoemulsions are promising nanomaterials for drug delivery applications in food, pharmaceutical and cosmetic industries. Despite the noteworthy commercial interest, little is known about their supramolecular organization, especially about how such multicomponent formulations interact with cell membranes. In the present work, coarse-grained molecular dynamics simulations have been employed to study the self-assembly of a 15-component lipid nanoemulsion droplet containing vitamins A and E for skin delivery. Our results display aspects of the unique "onion-like" agglomeration between the chemical constituents in the different layers of the lipid nanodroplet. Vitamin E molecules are more concentrated in the center of the droplet together with other hydrophobic constituents such as the triglycerides with long tails. On the other hand, vitamin A occupies an intermediate layer between the core and the co-emulsifier surface of the nanodroplet, together with lecithin phospholipids. Coarse-grained molecular dynamics simulations were also performed to provide insight into the first steps involved in absorption and penetration of the nanodroplet through skin membrane models, representing an intracellular (hair follicle infundibulum) and intercellular pathway (stratum corneum) through the skin. Our data provide a first view on the complex organization of commercial nanoemulsion and its interaction with skin membranes. We expect our results to open the way towards the rational design of such nanomaterials.

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.

The precise assembly and engineering of molecular machines capable of handling biomolecules play crucial roles in most single-molecule methods. In this work we use components from all three domains of life to fabricate an integrated multiprotein complex that controls the unfolding and threading of individual proteins across a nanopore. This 900 kDa multicomponent device was made in two steps. First, we designed a stable and low-noise β-barrel nanopore sensor by linking the transmembrane region of bacterial protective antigen to a mammalian proteasome activator. An archaeal 20S proteasome was then built into the artificial nanopore to control the unfolding and linearized transport of proteins across the nanopore. This multicomponent molecular machine opens the door to two approaches in single-molecule protein analysis, in which selected substrate proteins are unfolded, fed to into the proteasomal chamber and then addressed either as fragmented peptides or intact polypeptides., (© 2021. The Author(s), under exclusive licence to Springer Nature Limited.)

The coarse-grained Martini force field is widely used in biomolecular simulations. Here we present the refined model, Martini 3 ( http://cgmartini.nl ), with an improved interaction balance, new bead types and expanded ability to include specific interactions representing, for example, hydrogen bonding and electronic polarizability. The updated model allows more accurate predictions of molecular packing and interactions in general, which is exemplified with a vast and diverse set of applications, ranging from oil/water partitioning and miscibility data to complex molecular systems, involving protein-protein and protein-lipid interactions and material science applications as ionic liquids and aedamers.

Martini is a coarse-grained (CG) force field suitable for molecular dynamics (MD) simulations of (bio)molecular systems. It is based on mapping of two to four heavy atoms to one CG particle. The effective interactions between the CG particles are parametrized to reproduce partitioning free energies of small chemical compounds between polar and apolar phases. In this chapter, a summary of the key elements of this CG force field is presented, followed by an example of practical application: a lipoplex-membrane fusion experiment. Formulated as hands-on practice, this chapter contains guidelines to build CG models of important biological systems, such as asymmetric bilayers and double-stranded DNA. Finally, a series of notes containing useful information, limitations, and tips are described in the last section.