Coacervate formation studied by explicit solvent coarse-grain molecular dynamics with the Martini model

Complex coacervates are liquid–liquid phase separated systems, typically containing oppositely charged polyelectrolytes. They are widely studied for their functional properties as well as their potential involvement in cellular compartmentalization as biomolecular condensates. Diffusion and partitioning of solutes into a coacervate phase are important to address because their highly dynamic nature is one of their most important functional characteristics in real-world systems, but are difficult to study experimentally or even theoretically without an explicit representation of every molecule in the system. Here, we present an explicit-solvent, molecular dynamics coarse-grain model of complex coacervates, based on the Martini 3.0 force field. We demonstrate the accuracy of the model by reproducing the salt dependent coacervation of poly-lysine and poly-glutamate systems, and show the potential of the model by simulating the partitioning of ions and small nucleotides between the condensate and surrounding solvent phase. Our model paves the way for simulating coacervates and biomolecular condensates in a wide range of conditions, with near-atomic resolution.

Martini 3 force field can capture the experimental trends of complex coacervates and can be extended to gain physical insight on the mechanisms that drive the formation of LLPS.     

Comparing Dimerization Free Energies and Binding Modes of Small Aromatic Molecules with Different Force Fields

Dimerization free energies are fundamental quantities that describe the strength of interaction of different molecules. Obtaining accurate experimental values for small molecules and disentangling the conformations that contribute most to the binding can be extremely difficult, due to the size of the systems and the small energy differences. In many cases, one has to resort to computational methods to calculate such properties. In this work, we used molecular dynamics simulations in conjunction with metadynamics to calculate the free energy of dimerization of small aromatic rings, and compared three models from popular online servers for atomistic force fields, namely G54a7, CHARMM36 and OPLS. We show that, regardless of the force field, the profiles for the dimerization free energy of these compounds are very similar. However, significant care needs to be taken when studying larger molecules, since the deviations from the trends increase with the size of the molecules, resulting in force field dependent preferred stacking modes; for example, in the cases of pyrene and tetracene. Our results provide a useful background study for using topology builders to model systems which rely on stacking of aromatic moieties, and are relevant in areas ranging from drug design to supramolecular assembly.     

Association behaviour of glucitol amine gemini surfactants

The association behaviour of a number of glucitol amine gemini surfactants has been investigated by means of molecular dynamics and self-consistent-field calculations. We have shown that the titratable head group of the surfactant is responsible for a micelle-to-membrane transition when changing the pH. Furthermore, the association structure of this group of surfactants is shown to be very sensitive to ionic strength. The combination of a charged head group, a spacer, and the hydrophilic glucitol side chains is responsible for the possible structural transitions in the associates as a function of ionic strength and pH. Received: 12 January 2002 and Received in final form 10 March 2002     

Cholesterol shows preference for the interior of polyunsaturated lipid membranes

Recent neutron scattering experiments showed a striking manifestation of the aversion between polyunsaturated fatty acid (PUFA)-containing lipids and cholesterol. Selectively deuterated cholesterol/ 1,2-diarachidonylphosphatidylcholine (DAPC) samples revealed that the hydroxyl of the sterol resides at the center of the bilayer. Here we use a recently parametrized coarse grain simulation model to shed light on these puzzling experimental observations. Using a simulation setup in close correspondence to the experimental conditions, we reproduce the experimental neutron scattering profiles to a large extent. The simulations allow us to analyze the behavior of cholesterol in detail; we show that the interaction of cholesterol with the PUFA chains of DAPC leads to a fast flip-flop rate for the sterol and an increased preference of the sterol for the unusual location embedded between the monolayer leaflets.     

Structural determinants of the supramolecular organization of G protein-coupled receptors in bilayers

The G protein-coupled receptor (GPCR) rhodopsin self-assembles into supramolecular structures in native bilayers, but the structural determinants of receptor oligomerization are not known. We carried out multiple self-assembly coarse-grained molecular dynamics (CGMD) simulations of model membranes containing up to 64 molecules of the visual receptor rhodopsin over time scales reaching 100 μs. The simulations show strong preferential interaction modes between receptors. Two primary modes of receptor-receptor interactions are consistent with umbrella sampling/potential of mean force (PMF) calculations as a function of the distance between a pair of receptors. The preferential interfaces, involving helices (H) 1/8, 4/5 and 5, present no energy barrier to forming a very stable receptor dimer. Most notably, the PMFs show that the preferred rhodopsin dimer exists in a tail-to-tail conformation, with the interface comprising transmembrane H1/H2 and amphipathic H8 at the extracellular and cytoplasmic surfaces, respectively. This dimer orientation is in line with earlier electron microscopy, X-ray, and cross-linking experiments of rhodopsin and other GPCRs. Less stable interfaces, involving H4 and H6, have a free energy barrier for desolvation (delipidation) of the interfaces and appear to be designed to stabilize "lubricated" (i.e., lipid-coated) dimers. The overall CGMD strategy used here is general and can be applied to study the homo- and heterodimerization of GPCRs and other transmembrane proteins. Systematic extension of the work will deepen our understanding of the forces involved in the membrane organization of integral membrane proteins.     

Molecular dynamics simulation of the formation, structure, and dynamics of small phospholipid vesicles

Here, we use coarse grained molecular dynamics (MD) simulations to study the spontaneous aggregation of dipalmitoylphosphatidylcholine (DPPC) lipids into small unilamellar vesicles. We show that the aggregation process occurs on a nanosecond time scale, with bicelles and cuplike vesicles formed at intermediate stages. Formation of hemifused vesicles is also observed at higher lipid concentration. With either 25% dipalmitoylphosphatidylethanolamine (DPPE) or lysoPC mixed into the system, the final stages of the aggregation process occur significantly faster. The structure of the spontaneously formed vesicles is analyzed in detail. Microsecond simulations of isolated vesicles reveal significant differences in the packing of the lipids between the inner and outer monolayers, and between PC, PE, and lysoPC. Due to the small size of the vesicles they remain almost perfectly spherical, undergoing very limited shape fluctuations or bilayer undulations. The lipid lateral diffusion rate is found to be faster in the outer than in the inner monolayer. The water permeability coefficient of the pure DPPC vesicles is of the order of 10(-)(3) cm s(-)(1), in agreement with experimental measurements.     

Simulation of domain formation in DLPC-DSPC mixed bilayers

Binary mixtures of two phosphatidylcholines of different chain lengths are simulated in the bilayer state. We find a phase transition between a liquid state and a gel state at all concentrations. This phase transition is characterized by the area per lipid headgroup, the order parameter, and a change in dynamics. At concentrations with a majority of the longer lipid, we find phase separation into a gel and a liquid state in a small temperature window. This leads to a strong dynamic heterogeneity. Experimental phase transition temperatures are reproduced semiquantitatively. We see a clear shift in the phase transition to higher temperatures with increasing concentration of the longer lipid.     

The mechanism of vesicle fusion as revealed by molecular dynamics simulations

We describe molecular dynamics simulations elucidating the molecular details of the process of fusion for small lipid vesicles. The simulations are based on a coarse grained (CG) lipid model that accurately represents the lamellar state of a variety of phospholipids and enables us to observe intermediate stages during fusion at near atomic detail. Simulations were conducted on a variety of systems containing common phospholipids such as phosphatidylcholine (PC), phosphatidylethanolamine (PE), lysoPC, and mixtures of the above. The fusion intermediates found are in general agreement with the stalk-pore mechanism. Transient pores sometimes form adjacent to the stalk, however, resulting in the mixing of lipids from the outer and inner monolayers. The speed of stalk formation and the opening of the fusion pore can be modulated by altering the lipid composition in qualitative agreement with experimental observations.     

Application of mean field boundary potentials in simulations of lipid vesicles

A method is presented to enhance the efficiency of simulations of lipid vesicles. The method increases computational speed by eliminating water molecules that either surround the vesicle or reside in the interior of the vesicle, without altering the properties of the water at the membrane interface. Specifically, mean field force approximation (MFFA) boundary potentials are used to replace both the internal and external excess bulk solvent. In addition to reducing the cost of simulating preformed vesicles, the molding effect of the boundary potentials also enhances the formation and equilibration of vesicles from random solutions of lipid in water. Vesicles with diameters in the range from 20 to 60 nm were obtained on a nanosecond time scale, without any noticeable effect of the boundary potentials on their structure.     

Gating motions in voltage-gated potassium channels revealed by coarse-grained molecular dynamics simulations

Voltage-gated potassium (Kv) channels are ubiquitous transmembrane proteins involved in electric signaling of excitable tissues. A fundamental property of these channels is the ability to open or close in response to changes in the membrane potential. To date, their structure-based activation mechanism remains unclear, and there is a large controversy on how these gates function at the molecular level, in particular, how movements of the voltage sensor domain are coupled to channel gating. So far, all mechanisms proposed for this coupling are based on the crystal structure of the open voltage-gated Kv1.2 channel and structural models of the closed form based on electrophysiology experiments. Here, we use coarse-grain (CG) molecular dynamics simulations that allow conformational changes from the open to the closed form of the channel (embedded in its membrane environment) to be followed. Despite the low specificity of the CG force field, the obtained closed structure satisfies several experimental constraints. The overall results suggest a gating mechanism in which a lateral displacement the S4-S5 linker leads to a closing of the gate. Only a small up-down movement of the S4 helices is noticed. Additionally, the study suggests a peculiar upward motion of the intracellular tetramerization domain of the channel, hence providing a molecular view on how this domain may further regulate conduction in Kv channels.