Molecular structures of fluid phosphatidylethanolamine bilayers obtained from simulation-to-experiment comparisons and experimental scattering density profiles

Cholesterol Promotes Protein Binding by Affecting Membrane Electrostatics and Solvation Properties

Binding of the retroviral structural protein Gag to the cellular plasma membrane is mediated by the protein's matrix (MA) domain. Prominent among MA-PM interactions is electrostatic attraction between the positively charged MA domain and the negatively charged plasma membrane inner leaflet. Previously, we reported that membrane association of HIV-1 Gag, as well as purified Rous sarcoma virus (RSV) MA and Gag, depends strongly on the presence of acidic lipids and is enhanced by cholesterol (Chol). The mechanism underlying this enhancement was unclear. Here, using a broad set of in vitro and in silico techniques we addressed molecular mechanisms of association between RSV MA and model membranes, and investigated how Chol enhances this association. In neutron scattering experiments with liposomes in the presence or absence of Chol, MA preferentially interacted with preexisting POPS-rich clusters formed by nonideal lipid mixing, binding peripherally to the lipid headgroups with minimal perturbation to the bilayer structure. Molecular dynamics simulations showed a stronger MA-bilayer interaction in the presence of Chol, and a large Chol-driven increase in lipid packing and membrane surface charge density. Although in vitro MA-liposome association is influenced by disparate variables, including ionic strength and concentrations of Chol and charged lipids, continuum electrostatic theory revealed an underlying dependence on membrane surface potential. Together, these results conclusively show that Chol affects RSV MA-membrane association by making the electrostatic potential at the membrane surface more negative, while decreasing the penalty for lipid headgroup desolvation. The presented approach can be applied to other viral and nonviral proteins.

Cholesterol Promotes Protein Binding by Affecting Membrane Electrostatics and Solvation Properties

Binding of the retroviral structural protein Gag to the cellular plasma membrane is mediated by the protein's matrix (MA) domain. Prominent among MA-PM interactions is electrostatic attraction between the positively charged MA domain and the negatively charged plasma membrane inner leaflet. Previously, we reported that membrane association of HIV-1 Gag, as well as purified Rous sarcoma virus (RSV) MA and Gag, depends strongly on the presence of acidic lipids and is enhanced by cholesterol (Chol). The mechanism underlying this enhancement was unclear. Here, using a broad set of in vitro and in silico techniques we addressed molecular mechanisms of association between RSV MA and model membranes, and investigated how Chol enhances this association. In neutron scattering experiments with liposomes in the presence or absence of Chol, MA preferentially interacted with preexisting POPS-rich clusters formed by nonideal lipid mixing, binding peripherally to the lipid headgroups with minimal perturbation to the bilayer structure. Molecular dynamics simulations showed a stronger MA-bilayer interaction in the presence of Chol, and a large Chol-driven increase in lipid packing and membrane surface charge density. Although in vitro MA-liposome association is influenced by disparate variables, including ionic strength and concentrations of Chol and charged lipids, continuum electrostatic theory revealed an underlying dependence on membrane surface potential. Together, these results conclusively show that Chol affects RSV MA-membrane association by making the electrostatic potential at the membrane surface more negative, while decreasing the penalty for lipid headgroup desolvation. The presented approach can be applied to other viral and nonviral proteins.

Characterization of Ciprofloxacin Permeation Pathways across the Porin OmpC Using Metadynamics and a String Method

The rapid spreading of antimicrobial resistance in Gram-negative bacteria has become a major threat for humans as well as animals. As one of the main factors involved, the permeability of the outer membrane has attracted a great deal of attention recently. However, the knowledge regarding the translocation mechanisms for most available antibiotics is so far rather limited. Here, a theoretical study concerning the diffusion route of ciprofloxacin across the outer membrane porin OmpC from E. coli is presented. To this end, we establish a protocol to characterize meaningful permeation pathways by combining metadynamics with the zero-temperature string method. It was found that the lowest-energy pathway requires a reorientation of ciprofloxacin in the extracellular side of the porin before reaching the constriction region with its carboxyl group ahead. Several affinity sites have been identified, and their metastability has been evaluated using unbiased simulations. Such a detailed understanding is potentially very helpful in guiding the development of next generation antibiotics.

Drude Polarizable Force Field for Molecular Dynamics Simulations of Saturated and Unsaturated Zwitterionic Lipids

Additive force fields are designed to account for induced electronic polarization in a mean-field average way, using effective empirical fixed charges. The limitation of this approximation is cause for serious concerns, particularly in the case of lipid membranes, where the molecular environment undergoes dramatic variations over microscopic length scales. A polarizable force field based on the classical Drude oscillator offers a practical and computationally efficient framework for an improved representation of electrostatic interactions in molecular simulations. Building on the first-generation Drude polarizable force field for the dipalmitoylphosphatidylcholine 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) molecule, the present effort was undertaken to improve this initial model and expand the force field to a wider range of phospholipid molecules. New lipids parametrized include dimyristoylphosphatidylcholine (DMPC), dilauroylphosphatidylcholine (DLPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), dipalmitoylphosphatidylethanolamine (DPPE), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE), and 1,2-dioleoyl-sn-glycero-3-phosphoethanolamine (DOPE). The iterative optimization protocol employed in this effort led to lipid models that achieve a good balance between reproducing quantum mechanical data on model compound representative of phospholipids and reproducing a range of experimental condensed phase properties of bilayers. A parametrization strategy based on a restrained ensemble-maximum entropy methodology was used to help accurately match the experimental NMR order parameters in the polar headgroup region. All the parameters were developed to be compatible with the remainder of the Drude polarizable force field, which includes water, ions, proteins, DNA, and selected carbohydrates.

Complex biomembrane mimetics on the sub-nanometer scale

Biomimetic lipid vesicles are indispensable tools for gaining insight into the biophysics of cell physiology on the molecular level. The level of complexity of these model systems has steadily increased, and now spans from domain-forming lipid mixtures to asymmetric lipid bilayers. Here, we review recent progress in the development and application of elastic neutron and X-ray scattering techniques for studying these systems in situ and under physiologically relevant conditions on the nanometer to sub-nanometer length scales. In particular, we focus on: (1) structural details of coexisting liquid-ordered and liquid-disordered domains, including their thickness and lipid packing mismatch as a function of a size transition from nanoscopic to macroscopic domains; (2) membrane-mediated protein partitioning into lipid domains; (3) the role of the aqueous medium in tuning interactions between membranes and domains; and (4) leaflet-specific structure in asymmetric bilayers and passive lipid flip-flop.

Membrane Disruption Mechanism of a Prion Peptide (106-126) Investigated by Atomic Force Microscopy, Raman and Electron Paramagnetic Resonance Spectroscopy

A fragment of the human prion protein spanning residues 106-126 (PrP106-126) recapitulates many essential properties of the disease-causing protein such as amyloidogenicity and cytotoxicity. PrP106-126 has an amphipathic characteristic that resembles many antimicrobial peptides (AMPs). Therefore, the toxic effect of PrP106-126 could arise from a direct association of monomeric peptides with the membrane matrix. Several experimental approaches are employed to scrutinize the impacts of monomeric PrP106-126 on model lipid membranes. Porous defects in planar bilayers are observed by using solution atomic force microscopy. Adding cholesterol does not impede defect formation. A force spectroscopy experiment shows that PrP106-126 reduces Young's modulus of planar lipid bilayers. We use Raman microspectroscopy to study the effect of PrP106-126 on lipid atomic vibrational dynamics. For phosphatidylcholine lipids, PrP106-126 disorders the intrachain conformation, while the interchain interaction is not altered; for phosphatidylethanolamine lipids, PrP106-126 increases the interchain interaction, while the intrachain conformational order remains similar. We explain the observed differences by considering different modes of peptide insertion. Finally, electron paramagnetic resonance spectroscopy shows that PrP106-126 progressively decreases the orientational order of lipid acyl chains in magnetically aligned bicelles. Together, our experimental data support the proposition that monomeric PrP106-126 can disrupt lipid membranes by using similar mechanisms found in AMPs.

A novel pathway for amyloids self-assembly in aggregates at nanomolar concentration mediated by the interaction with surfaces

A limitation of the amyloid hypothesis in explaining the development of neurodegenerative diseases is that the level of amyloidogenic polypeptide in vivo is below the critical concentration required to form the aggregates observed in post-mortem brains. We discovered a novel, on-surface aggregation pathway of amyloidogenic polypeptide that eliminates this long-standing controversy. We applied atomic force microscope (AFM) to demonstrate directly that on-surface aggregation takes place at a concentration at which no aggregation in solution is observed. The experiments were performed with the full-size Aβ protein (Aβ42), a decapeptide Aβ(14-23) and α-synuclein; all three systems demonstrate a dramatic preference of the on-surface aggregation pathway compared to the aggregation in the bulk solution. Time-lapse AFM imaging, in solution, show that over time, oligomers increase in size and number and release in solution, suggesting that assembled aggregates can serve as nuclei for aggregation in bulk solution. Computational modeling performed with the all-atom MD simulations for Aβ(14-23) peptide shows that surface interactions induce conformational transitions of the monomer, which facilitate interactions with another monomer that undergoes conformational changes stabilizing the dimer assembly. Our findings suggest that interactions of amyloidogenic polypeptides with cellular surfaces play a major role in determining disease onset.

Lipid bilayer thickness determines cholesterol's location in model membranes

Cholesterol is an essential biomolecule of animal cell membranes, and an important precursor for the biosynthesis of certain hormones and vitamins. It is also thought to play a key role in cell signaling processes associated with functional plasma membrane microdomains (domains enriched in cholesterol), commonly referred to as rafts. In all of these diverse biological phenomena, the transverse location of cholesterol in the membrane is almost certainly an important structural feature. Using a combination of neutron scattering and solid-state ²H NMR, we have determined the location and orientation of cholesterol in phosphatidylcholine (PC) model membranes having fatty acids of different lengths and degrees of unsaturation. The data establish that cholesterol reorients rapidly about the bilayer normal in all the membranes studied, but is tilted and forced to span the bilayer midplane in the very thin bilayers. The possibility that cholesterol lies flat in the middle of bilayers, including those made from PC lipids containing polyunsaturated fatty acids (PUFAs), is ruled out. These results support the notion that hydrophobic thickness is the primary determinant of cholesterol's location in membranes.

An extensive simulation study of lipid bilayer properties with different head groups, acyl chain lengths, and chain saturations

Previous MD simulations of six phosphocholine (PC) lipid bilayers demonstrated the accuracy of the CHARMM36 force field (C36FF) for PC bilayer simulation at varied temperatures (BBA-Biomembranes, 1838 (2014): 2520-2529). In this work, we further examine the accuracy of C36FF over a wide temperature range for a broader range of lipid types such as various head groups (phosphatidic acid (PA), PC, phosphoethanolamine (PE), phosphoglycerol (PG), and phosphoserine (PS)), and tails (saturated, mono-, mixed- and poly-unsaturated acyl chains with varied chain lengths). The structural properties (surface area per lipid (SA/lip), overall bilayer thickness, hydrophobic thickness, headgroup-to-headgroup thickness, deuterium order parameter (S_CD), and spin-lattice relaxation time (T₁)) obtained from simulations agree well with nearly all available experimental data. Our analyses indicate that PS lipids have the most inter-lipid hydrogen bonds, while PG lipids have the most intra-lipid hydrogen bonds, which play the main role in their low SA/lip in PS lipids and low thicknesses in PG lipids, respectively. PS, PE, and PA lipids have the largest contact clusters with on average 5-8 lipids per cluster, while PC and PG have clusters of 4 lipids based on a cutoff distance of 6.5Å. PS lipids have much slower lipid wobble (i.e., higher correlation time) than other head groups at a given temperature as the hydrogen bonded network significantly reduces a lipid's mobility, and the rate of lipid wobble increases dramatically as temperature increases. These in-depth analyses facilitate further understanding of lipid bilayers at the atomic level.

Simulation of lipid bilayer self-assembly using all-atom lipid force fields

In this manuscript we expand significantly on our earlier communication by investigating the bilayer self-assembly of eight different types of phospholipids in unbiased molecular dynamics (MD) simulations using three widely used all-atom lipid force fields. Irrespective of the underlying force field, the lipids are shown to spontaneously form stable lamellar bilayer structures within 1 microsecond, the majority of which display properties in satisfactory agreement with the experimental data. The lipids self-assemble via the same general mechanism, though at formation rates that differ both between lipid types, force fields and even repeats on the same lipid/force field combination. In addition to zwitterionic phosphatidylcholine (PC) and phosphatidylethanolamine (PE) lipids, anionic phosphatidylserine (PS) and phosphatidylglycerol (PG) lipids are represented. To our knowledge this is the first time bilayer self-assembly of phospholipids with negatively charged head groups is demonstrated in all-atom MD simulations.

On scattered waves and lipid domains: detecting membrane rafts with X-rays and neutrons

In order to understand the biological role of lipids in cell membranes, it is necessary to determine the mesoscopic structure of well-defined model membrane systems. Neutron and X-ray scattering are non-invasive, probe-free techniques that have been used extensively in such systems to probe length scales ranging from angstroms to microns, and dynamics occurring over picosecond to millisecond time scales. Recent developments in the area of phase separated lipid systems mimicking membrane rafts will be presented, and the underlying concepts of the different scattering techniques used to study them will be discussed in detail.

Structural Significance of Lipid Diversity as Studied by Small Angle Neutron and X-ray Scattering

We review recent developments in the rapidly growing field of membrane biophysics, with a focus on the structural properties of single lipid bilayers determined by different scattering techniques, namely neutron and X-ray scattering. The need for accurate lipid structural properties is emphasized by the sometimes conflicting results found in the literature, even in the case of the most studied lipid bilayers. Increasingly, accurate and detailed structural models require more experimental data, such as those from contrast varied neutron scattering and X-ray scattering experiments that are jointly refined with molecular dynamics simulations. This experimental and computational approach produces robust bilayer structural parameters that enable insights, for example, into the interplay between collective membrane properties and its components (e.g., hydrocarbon chain length and unsaturation, and lipid headgroup composition). From model studies such as these, one is better able to appreciate how a real biological membrane can be tuned by balancing the contributions from the lipid's different moieties (e.g., acyl chains, headgroups, backbones, etc.).

Biomembranes research using thermal and cold neutrons

In 1932 James Chadwick discovered the neutron using a polonium source and a beryllium target (Chadwick, 1932). In a letter to Niels Bohr dated February 24, 1932, Chadwick wrote: "whatever the radiation from Be may be, it has most remarkable properties." Where it concerns hydrogen-rich biological materials, the "most remarkable" property is the neutron's differential sensitivity for hydrogen and its isotope deuterium. Such differential sensitivity is unique to neutron scattering, which unlike X-ray scattering, arises from nuclear forces. Consequently, the coherent neutron scattering length can experience a dramatic change in magnitude and phase as a result of resonance scattering, imparting sensitivity to both light and heavy atoms, and in favorable cases to their isotopic variants. This article describes recent biomembranes research using a variety of neutron scattering techniques.