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.     

Multiscale modeling of molecular structure and optical properties of complex supramolecular aggregates

Supramolecular aggregates of synthetic dye molecules offer great perspectives to prepare biomimetic functional materials for light-harvesting and energy transport. The design is complicated by the fact that structure–property relationships are hard to establish, because the molecular packing results from a delicate balance of interactions and the excitonic properties that dictate the optics and excited state dynamics, in turn sensitively depend on this packing. Here we show how an iterative multiscale approach combining molecular dynamics and quantum mechanical exciton modeling can be used to obtain accurate insight into the packing of thousands of cyanine dye molecules in a complex double-walled tubular aggregate in close interaction with its solvent environment. Our approach allows us to answer open questions not only on the structure of these prototypical aggregates, but also about their molecular-scale structural and energetic heterogeneity, as well as on the microscopic origin of their photophysical properties. This opens the route to accurate predictions of energy transport and other functional properties.

Multiscale modeling resolves the molecular structure of a synthetic light-harvesting complex, unraveling the microscopic origin of its photophysical properties.

Supramolecular structures may self-assemble from a variety of building blocks, resulting in a wide range of advanced materials with attractive biomimetic, sensing, catalytic, optoelectronic and photonic functionalities.^1–10 The close-packed nanoscale organization of the individual molecules within a supramolecular system, held together via noncovalent interactions, gives rise to the aggregate''s (collective) properties. Assemblies consisting of dye molecules often exhibit unique collective optical properties and are of interest for opto-electronic applications as well as artificial light-harvesting complexes that mimic natural antenna systems of photosynthetic bacteria and plants.^11–13 For example, chlorosomal antenna complexes of photosynthetic green sulfur bacteria are self-assembled into multilayer tubular structures having bacteriochlorophyll pigments as building blocks.^14–16 The structure of these antenna complexes and the underlying molecular arrangement ensures that the process of light-harvesting and excitation energy transport is very efficient, even under extremely low light conditions.^17,18 The quest to recreate such efficiency under laboratory conditions has sparked numerous studies of synthetic self-assembled systems mimicking natural chlorosomes, e.g. using porphyrins,¹⁹ zinc chlorin,²⁰ and cyanine dyes.²¹ Of particular interest are the tubular aggregates of 3,3′-bis(2-sulfopropyl)-5,5′,6,6′-tetrachloro-1,1′-dioctylbenzimidacarbocyanine (C8S3).^22–25 Cryo-TEM reveals a hierarchy of supramolecular architectures, including double-walled nanotubes; under certain conditions, bundles of nanotubes arise.²⁶ Thus, this system allows for the occurrence of electronic excitation energy transport at various levels: within one wall, between walls of one tube, and between different tubes, similar to the situation in natural systems.^27,28To understand how such supramolecular systems work, as well as propose design rules for new materials, it is essential to determine the relationship between molecular structure and optical properties. Current experimental techniques, however, are unable to resolve the structure at the molecular level. This, in combination with the sensitivity of spectral properties to the details of the molecular packing, leads to a crucial role for theoretical modeling.²⁹ For example, molecular dynamics (MD) simulations have been used to predict the molecular packing within a variety of supramolecular assemblies.^30–34 However, synthetic amphiphiles with aromatic groups, such as cyanine dyes—often used to prepare aggregates with optical functionality—tend to fall into kinetic traps during spontaneous self-assembly simulations and the packing of the aromatic chromophores remains highly disordered on the accessible time scale, leading to predicted (optical) spectra that are not consistent with experimental data.³⁵ This problem can be overcome by building assemblies based upon proposed architectures and assessing their stability in relatively short MD simulations.^36–38 The drawback of this approach is the requirement of a thorough understanding of what to use as a starting point and how to validate the structure. In any case, proper validation requires the modeling of the optical spectra of the obtained structure, and finally, comparing it to the experiment. The demanding character of such methods explains why an important role is played by phenomenological modeling, in which a molecular packing is guessed and the optics is obtained from parametrizing an exciton model that describes the collective excited states of the assembly with interactions dictated by the guessed packing. By comparing the calculated spectra to experimental ones, the structure and exciton model may be fine-tuned. While this method has been successful in describing spectra,^23,39 it is limited in its predictive power and also lacks access to essential microscopic parameters, such as tuning of the optical excitation energies imposed by the environment, disorder in these energies and structural heterogeneity.In this work, we use an advanced multiscale approach to determine structure–optical property relationships for the C8S3 double-walled nanotubes, guided by comparison to experiments. The optical spectrum of these aggregates, in which multiple exciton peaks may be discerned, suggests a rather complex underlying molecular packing. This fact, combined with their sheer size going up to many thousands of molecules, makes these systems exceptionally challenging to resolve and leaves important questions concerning structure–function relationships unanswered or under debate, for instance the origin of the splitting between the two lowest-energy spectral bands.^23,38 Here, we answer these questions by iteratively combining MD simulations to capture the details of molecular packing and structural disorder, an exciton Hamiltonian approach to calculate optical signatures, and explicit microelectrostatic calculations to estimate energetic disorder and solvent shifts. Previous attempts to reveal the structure of cyanine-based nanotubes were limited to small-scale system sizes,^37,38 modeling optical features phenomenologically rather than using atomistic information³⁸ or featuring simpler, single-walled systems.³⁷ In addition to answering important questions for the C8S3 double-walled nanotubes, our study opens the way to explain and predict at an unprecedented level of detail the functional properties of other highly complex molecular materials.     

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.     

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.     

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.     

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.     

Reconstruction of atomistic details from coarse‐grained structures

We present an algorithm to reconstruct atomistic structures from their corresponding coarse‐grained (CG) representations and its implementation into the freely available molecular dynamics (MD) program package GROMACS. The central part of the algorithm is a simulated annealing MD simulation in which the CG and atomistic structures are coupled via restraints. A number of examples demonstrate the application of the reconstruction procedure to obtain low‐energy atomistic structural ensembles from their CG counterparts. We reconstructed individual molecules in vacuo (NCQ tripeptide, dipalmitoylphosphatidylcholine, and cholesterol), bulk water, and a WALP transmembrane peptide embedded in a solvated lipid bilayer. The first examples serve to optimize the parameters for the reconstruction procedure, whereas the latter examples illustrate the applicability to condensed‐phase biomolecular systems. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

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.