Stabilization of Prebiotic Vesicles by Peptides Depends on Sequence and Chirality: A Mechanism for Selection of Protocell-Associated Peptides.

Cells require oligonucleotides and polypeptides with specific, homochiral sequences to perform essential functions, but it is unclear how such oligomers were selected from random sequences at the origin of life. Cells were probably preceded by simple compartments such as fatty acid vesicles, and oligomers that increased the stability, growth, or division of vesicles could have thereby increased in frequency. We therefore tested whether prebiotic peptides alter the stability or growth of vesicles composed of a prebiotic fatty acid. We find that three of 15 dipeptides tested reduce salt-induced flocculation of vesicles. All three contain leucine, and increasing their length increases the efficacy. Also, leucine-leucine but not alanine-alanine increases the size of vesicles grown by multiple additions of micelles. In a molecular simulation, leucine-leucine docks to the membrane, with the side chains inserted into the hydrophobic core of the bilayer, while alanine-alanine fails to dock. Finally, the heterochiral forms of leucine-leucine, at a high concentration, rapidly shrink the vesicles and make them leakier and less stable to high pH than the homochiral forms do. Thus, prebiotic peptide-membrane interactions influence the flocculation, growth, size, leakiness, and pH stability of prebiotic vesicles, with differential effects due to sequence, length, and chirality. These differences could lead to a population of vesicles enriched for peptides with beneficial sequence and chirality, beginning selection for the functional oligomers that underpin life.

Rapid Formation of Peptide/Lipid Coaggregates by the Amyloidogenic Seminal Peptide PAP248-286.

Protein/lipid coassembly is an understudied phenomenon that is important to the function of antimicrobial peptides as well as the pathological effects of amyloid. Here, we study the coassembly process of PAP248-286, a seminal peptide that displays both amyloid-forming and antimicrobial activity. PAP248-286 is a peptide fragment of prostatic acid phosphatase and has been reported to form amyloid fibrils, known as semen-derived enhancer of viral infection (SEVI), that enhance the viral infectivity of human immunodeficiency virus. We find that in addition to forming amyloid, PAP248-286 much more readily assembles with lipid vesicles into peptide/lipid coaggregates that resemble amyloid fibrils in some important ways but are a distinct species. The formation of these PAP248-286/lipid coaggregates, which we term "messicles," is controlled by the peptide:lipid (P:L) ratio and by the lipid composition. The optimal P:L ratio is around 1:10, and at least 70% anionic lipid is required for coaggregate formation. Once formed, messicles are not disrupted by subsequent changes in P:L ratio. We propose that messicles form through a polyvalent assembly mechanism, in which a critical surface density of PAP248-286 on liposomes enables peptide-mediated particle bridging into larger species. Even at ∼50-fold lower PAP248-286 concentrations, messicles form at least 10-fold faster than amyloid fibrils. It is therefore possible that some or all of the biological activities assigned to SEVI, the amyloid form of PAP248-286, could instead be attributed to a PAP248-286/lipid coaggregate. More broadly speaking, this work could provide a potential framework for the discovery and characterization of nonamyloid peptide/lipid coaggregates by other amyloid-forming proteins and antimicrobial peptides.

Modulated and spiral surface patterns on deformable lipid vesicles.

We investigate the behavior of two-dimensional systems that exhibit a transition between homogeneous and spatially inhomogeneous phases, which have spherical topology, and whose mechanical properties depend on the local value of the order parameter. One example of such a system is multicomponent lipid bilayer vesicles, which serve as a model to study cellular membranes. Under certain conditions, such bilayers separate into coexisting liquid-ordered and liquid-disordered regions. When arranged into the shape of small vesicles, this phase coexistence can result in spatial patterns that are more complex than the basic two-domain configuration encountered in typical bulk systems. The difference in bending rigidity between the liquid-ordered and liquid-disordered regions couples the shape of the vesicle to the local composition. We show that this interplay gives rise to a rich phase diagram that includes homogeneous, separated, and axisymmetric modulated phases that are divided by regions of spiral patterns in the surface morphology.

Mechanisms of DNA hybridization: Transition path analysis of a simulation-informed Markov model.

Complementary deoxyribonucleic acid (DNA) strands in solution reliably hybridize to form stable duplexes. We study the kinetics of the hybridization process and the mechanisms by which two initially isolated strands come together to form a stable double helix. We adopt a multi-step computational approach. First, we perform a large number of Brownian dynamics simulations of the hybridization process using the coarse-grained oxDNA2 model. Second, we use these simulations to construct a Markov state model of DNA dynamics that uses a state decomposition based on the inter-strand hydrogen bonding pattern. Third, we take advantage of transition path theory to obtain quantitative information about the thermodynamic and dynamic properties of the hybridization process. We find that while there is a large ensemble of possible hybridization pathways, there is a single dominant mechanism in which an initial base pair forms close to either end of the nascent double helix, and the remaining bases pair sequentially in a zipper-like fashion. We also show that the number of formed base pairs by itself is insufficient to describe the transition state of the hybridization process.

Spatially Resolving the Condensing Effect of Cholesterol in Lipid Bilayers.

We study the effect of cholesterol on the structure of dipalmitoylphosphatidylcholine phospholipid bilayers. Using extensive molecular dynamics computer simulations at atomistic resolution, we observe and quantify several structural changes upon increasing cholesterol content that are collectively known as the condensing effect: a thickening of the bilayer, an increase in lipid tail order, and a decrease in lateral area. We also observe a change in leaflet interdigitation and a lack thereof in the distributions of dipalmitoylphosphatidylcholine headgroup orientations. These results, obtained over a wide range of cholesterol mole fractions, are then used to calibrate the analysis of phospholipid properties in bilayers containing a single cholesterol molecule per leaflet, which we perform in a spatially resolved way. We find that a single cholesterol molecule affects phospholipids in its first and second solvation shells, which puts the range of this effective interaction to be on the order of 1-2 nm. We also observe a tendency of phospholipids to orient their polar headgroups toward the cholesterol, which provides additional support for the umbrella model of bilayer organization.

Relating the structure factors of two-dimensional materials in planar and spherical geometries.

Scattering structure factors provide essential insight into material properties and are routinely obtained in experiments, computer simulations, and theoretical analyses. Different approaches favor different geometries of the material. In case of lipid bilayers, scattering experiments can be performed on spherical vesicles, while simulations and theory often consider planar membrane patches. We derive an approximate relationship between the structure functions of such a material in planar and spherical geometries. We illustrate the usefulness of this relationship in a case study of a Gaussian material that supports both homogeneous and microemulsion phases. Within its range of applicability, this relationship enables a model-free comparison of structure factors of the same material in different geometries.

Phase diagrams of multicomponent lipid vesicles: Effects of finite size and spherical geometry.

We study the phase behavior of multicomponent lipid bilayer vesicles that can exhibit intriguing morphological patterns and lateral phase separation. We use a modified Landau-Ginzburg model capable of describing spatially uniform phases, microemulsions, and modulated phases on a spherical surface. We calculate its phase diagram for multiple vesicle sizes using analytical and numerical techniques as well as Monte Carlo simulations. Consistent with previous studies on planar systems, we find that thermal fluctuations move phase boundaries, stabilizing phases of higher disorder. We also show that the phase diagram is sensitive to the size of the system at small vesicle radii. Such finite size effects are likely relevant in experiments on small, unilamellar vesicles and should be considered in their comparison to theoretical and simulation results.

Macroscopic phase separation, modulated phases, and microemulsions: a unified picture of rafts.

We simulate a simple phenomenological model describing phase behavior in a multicomponent membrane, a model capable of producing macroscopic phase separation, modulated phases, and microemulsions, all of which have been discussed in terms of raft phenomena. We show that one effect of thermal fluctuations on the mean-field phase diagram is that it permits a direct transition between either one of the coexisting liquid phases to a microemulsion. This implies that one system exhibiting phase separation can be related to a similar system exhibiting the heterogeneities characteristic of a microemulsion. The two systems could differ in their average membrane composition or in the relative compositions of their exoplasmic and cytoplasmic leaves. The model provides a unified description of these raft-associated phenomena.

Building Markov state models with solvent dynamics.

BACKGROUND: Markov state models have been widely used to study conformational changes of biological macromolecules. These models are built from short timescale simulations and then propagated to extract long timescale dynamics. However, the solvent information in molecular simulations are often ignored in current methods, because of the large number of solvent molecules in a system and the indistinguishability of solvent molecules upon their exchange. METHODS: We present a solvent signature that compactly summarizes the solvent distribution in the high-dimensional data, and then define a distance metric between different configurations using this signature. We next incorporate the solvent information into the construction of Markov state models and present a fast geometric clustering algorithm which combines both the solute-based and solvent-based distances. RESULTS: We have tested our method on several different molecular dynamical systems, including alanine dipeptide, carbon nanotube, and benzene rings. With the new solvent-based signatures, we are able to identify different solvent distributions near the solute. Furthermore, when the solute has a concave shape, we can also capture the water number inside the solute structure. Finally we have compared the performances of different Markov state models. The experiment results show that our approach improves the existing methods both in the computational running time and the metastability. CONCLUSIONS: In this paper we have initiated an study to build Markov state models for molecular dynamical systems with solvent degrees of freedom. The methods we described should also be broadly applicable to a wide range of biomolecular simulation analyses.

The effect of spin polarization on double electron-electron resonance (DEER) spectroscopy.

Double electron-electron resonance (DEER) spectroscopy measures the distribution of distances between two electron spins in the nanometer range, often on doubly spin-labeled proteins, via the modulation of a refocused spin echo by the dipolar interaction between the spins. DEER is commonly conducted under conditions where the polarization of the spins is small. Here, we examine the DEER signal under conditions of high spin polarization, thermally obtainable at low temperatures and high magnetic fields, and show that the signal acquires a polarization-dependent out-of-phase component both for the intramolecular and intermolecular contributions. For the latter, this corresponds to a phase shift of the spin echo that is linear in the pump pulse position. We derive a compact analytical form of this phase shift and show experimental measurements using monoradical and biradical nitroxides at several fields and temperatures. The effect highlights a novel aspect of the fundamental spin physics underlying DEER spectroscopy.

Effects of Salts on the Solvation of Hydrophobic Objects in Water.

The solvation of large, hydrophobic objects in water is facilitated by the formation of a low-density region surrounding the solute that is separated from the bulk liquid by an interface, which has a structure that resembles that between a liquid and its vapor. We study the effect of dissolved sodium chloride on the thermodynamics of solvation and on the solvent structure surrounding hydrophobic solutes in the size regime where this interface is not yet fully formed. Using biased Molecular Dynamics computer simulations, we calculate solvation free energies and orientational distributions of water molecules at different salt concentrations and solute sizes. We find that while the effects of sodium chloride on thermodynamic properties are small, the ions' response to the presence of a hydrophobic solute differs significantly from that of the water. Our findings provide mechanistic insight into how our understanding of hydrophobic solvation in water can be extended to electrolyte solutions.

Large-scale simulations of fluctuating biological membranes.

We present a simple, and physically motivated, coarse-grained model of a lipid bilayer, suited for micron scale computer simulations. Each approximately 25 nm(2) patch of bilayer is represented by a spherical particle. Mimicking forces of hydrophobic association, multiparticle interactions suppress the exposure of each sphere's equator to its implicit solvent surroundings. The requirement of high equatorial density stabilizes two-dimensional structures without necessitating crystalline order, allowing us to match both the elasticity and fluidity of natural lipid membranes. We illustrate the model's versatility and realism by characterizing a membrane's response to a prodding nanorod.

Effect of alcohol on the phase separation in model membranes.

The discovery of coexisting liquid-ordered and liquid-disordered phases in multicomponent lipid bilayers has received widespread attention due to its potential relevance for biological systems. One of the many open questions is how the presence of additional components affects the nature of the coexisting phases. Of particular interest is the addition of alcohols because their anesthetic properties may arise from modulating bilayer behavior. We use coarse-grained Molecular Dynamics simulations to gain insight into the partitioning preferences of linear n-alcohols into ordered and disordered bilayers alongside their effects on local membrane structure. We find that alcohols cause only small changes to membrane composition alongside a lack of significant effects on membrane thickness and lipid tail order. Cholesterol and n-alcohol trans-bilayer motion is measured and found to be near or within the range of previous atomistic results. The cholesterol flip-flop rates increase with both n-alcohol length and concentration for octanol, dodecanol, and hexadecanol, indicating a decrease in lipid order. Umbrella sampling simulations of removing cholesterol from tertiary membranes find no significant difference with or without n-alcohols at various concentrations. Simulations of a phase-separated bilayer show that octanol preferentially partitions into the liquid-disordered phase in a ratio of approximately 3:1 over the liquid-ordered phase. Furthermore, partition coefficients of alcohol in single-phase membranes show a preference for longer alcohols (dodecanol and hexadecanol) to partition preferentially into the liquid-ordered phase, while decreasing the length of the alcohol reverses this trend. Our work tests experimental results while also investigating the ability for coarse-grained MARTINI simulations to capture minute differences in model membrane spatial arrangements on the nanoscale level.

Negatively Charged Lipids Exhibit Negligible Effects on the Water Repellency of Montmorillonite Films.

Amphiphilic molecules can alter the wettability of soil minerals. To determine how the headgroup chemistry of amphiphiles determines these effects, we investigate a system of the clay montmorillonite with long-chain phospholipids. We use phosphatidylglycerol (PG) phospholipids to contrast with our previous work using phosphatidylethanolamine (PE) lipids. Zwitterionic PE lipids can sorb to the negatively charged montmorillonite surface, whereas negatively charged PG lipids cannot. Employing a suite of techniques from molecular dynamics, atomic force microscopy, fluorescence microscopy, and contact angle measurements, we define sample characteristics from molecular-scale structure to the macroscopic wettability. We find that PG lipids do not significantly alter montmorillonite wetting characteristics, such as the contact angle, flow viscosity, and the characteristic time scale for droplet imbibition. On comparing PE and PG lipid/clay films, we find that, among the phospholipids compared, they must have three characteristics to change clay/lipid film wettability: they must bind to the mineral surface, be solid at room temperature, and have a relatively continuous distribution throughout the film.

Laterally Resolved Small-Angle Scattering Intensity from Lipid Bilayer Simulations: An Exact and a Limited-Range Treatment.

When combined, molecular simulations and small-angle scattering experiments are able to provide molecular-scale resolution of structure. Separately, scattering experiments provide only intermingled pair correlations between atoms, while molecular simulations are limited by model quality and the relatively short time scales that they can access. Their combined strength relies on agreement between the experimental spectra and those computed by simulation. To date, computing the neutron spectra from a molecular simulation of a lipid bilayer is straightforward only if the structure is approximated by laterally averaging the in-plane bilayer structure. However, this neglects all information about lateral heterogeneity, e.g., clustering of components in a lipid mixture. This paper presents two methods for computing the scattering intensity of simulated bilayers with in-plane heterogeneity, enabling a full treatment of both the transverse and lateral bilayer structure for the first time. The first method, termed the Dirac Brush, computes the exact spectra including spurious artifacts resulting from using information from neighboring periodic cells to account for the long-range structure of the bilayer. The second method, termed PFFT, applies a mean-field treatment in the field far from a scattering element, resulting in a correlation range that can be tuned (eliminating correlations with neighboring periodic images), but with computational cost that prohibits obtaining the exact (Dirac Brush) spectra. Following their derivation, the two methods are applied to a coarse-grained molecular simulation of a bilayer inhomogeneity, demonstrating the contributions of lateral correlations to the resulting spectra.

Connecting wettability, topography, and chemistry in a simple lipid-montmorillonite system.

HYPOTHESIS: While soil water repellency causes a variety of undesirable environmental effects, the underlying mechanism is unknown. We investigate the coupled effects of chemical characteristics and surface topology in a simple model system of two lipids, DSPE (1,2-distearoyl-sn-glycero-3-phosphoethanolamine) and DOPE (1,2-dioleoyl-sn-glycero-3-phosphoethanolamine), and a clay substrate. These closely-related lipids allowed the study of how a small change in chemical structure influences the surface hydrophobicity. EXPERIMENTS: Techniques ranging from molecular (simulations) to nanoscopic (atomic force microscopy) to microscopic (fluorescence microscopy) to macroscopic (contact angle measurements) were used to explore interactions at all length scales. The wettability was assessed from initial contact angle and time-dependent changes in droplet shape. FINDINGS: The lipid distribution depended on the lipid's melting temperature: solid lipids did not spread evenly through the film, while liquid ones did. However, the initial contact angle did not change appreciably with the addition of DSPE or DOPE. Only DSPE heated above its melting temperature induced significant changes. In addition to the initial contact angle, quantitative variables extracted from the change in droplet shape over time correlated with the film topography or lipid distribution. These results define a new quantitative approach to investigating partially-wettable soils and provide a potential rationale for why clays can remediate water-repellent soils.

DNA Base Pair Mismatches Induce Structural Changes and Alter the Free-Energy Landscape of Base Flip.

Double-stranded DNA may contain mismatched base pairs beyond the Watson-Crick pairs, guanine-cytosine and adenine-thymine. Such mismatches bear adverse consequences for human health. We utilize molecular dynamics and metadynamics computer simulations to study the structure and dynamics of both matched and mismatched base pairs. We discover significant differences between the matched and mismatched pairs in structure and base flip work profiles. Mismatched pairs shift more in the plane normal to the DNA strand and exhibit more noncanonical structures, including the e-motif. We discuss the potential implications on the mismatch repair enzymes' detection of DNA mismatches.

Identifying the Onset of Phase Separation in Quaternary Lipid Bilayer Systems from Coarse-Grained Simulations.

Understanding the (de)mixing behavior of multicomponent lipid bilayers is an important step toward unraveling the nature of spatial composition heterogeneities in cellular membranes and their role in biological function. We use coarse-grained molecular dynamics simulations to study the composition phase diagram of a quaternary mixture of phospholipids and cholesterol. This mixture is known to exhibit both uniform and coexisting phases. We compare and combine different statistical measures of membrane structure to identify the onset of phase coexistence in composition space. An important element in our approach is the dependence of composition heterogeneities on the size of the system. While homogeneous phases can be structured and display long correlation lengths, the hallmark behavior of phase coexistence is the scaling of the apparent correlation length with system size. Because the latter cannot be easily varied in simulations, our method instead uses information obtained from observation windows of different sizes to accurately distinguish phase coexistence from structured homogeneous phases. This approach is built on very general physical principles, and will be beneficial to future studies of the phase behavior of multicomponent lipid bilayers.

Free Energy Calculations of Membrane Permeation: Challenges Due to Strong Headgroup-Solute Interactions.

Understanding how different classes of molecules move across biological membranes is a prerequisite to predicting a solute's permeation rate, which is a critical factor in the fields of drug design and pharmacology. We use biased molecular dynamics computer simulations to calculate and compare the free energy profiles of translocation of several small molecules across 1,2-dioleoyl- sn-glycero-3-phosphocholine (DOPC) lipid bilayers as a first step toward determining the most efficient method for free energy calculations. We study the translocation of arginine, a sodium ion, alanine, and a single water molecule using the metadynamics, umbrella sampling, and replica exchange umbrella sampling techniques. Within the fixed lengths of our simulations, we find that all methods produce similar results for charge-neutral permeants, but not for polar or positively charged molecules. We identify the long relaxation time scale of electrostatic interactions between lipid headgroups and the solute to be the principal cause of this difference and show that this slow process can lead to an erroneous dependence of computed free energy profiles on the initial system configuration. We demonstrate the use of committor analysis to validate the proper sampling of the presumed transition state, which in our simulations is achieved only in replica exchange calculations. On the basis of these results we provide some useful guidance to perform and evaluate free energy calculations of membrane permeation.

Segue between favorable and unfavorable solvation.

Solvation of small and large clusters are studied by simulation, considering a range of solvent-solute attractive energy strengths. Over a wide range of conditions, both for solvation in the Lennard-Jones liquid and in the SPC model of water, it is shown that the mean solvent density varies linearly with changes in solvent-solute adhesion or attractive energy strength. This behavior is understood from the perspective of Weeks' theory of solvation (Annu. Rev. Phys. Chem. 2002, 53, 533) and supports theories based upon that perspective.