Jonathan N. Sachs

Lipid Bilayer Structure Determined by the Simultaneous Analysis of Neutron and X-Ray Scattering Data

Quantitative structures were obtained for the fully hydrated fluid phases of dioleoylphosphatidylcholine (DOPC) and dipalmitoylphosphatidylcholine (DPPC) bilayers by simultaneously analyzing x-ray and neutron scattering data. The neutron data for DOPC included two solvent contrasts, 50% and 100% D(2)O. For DPPC, additional contrast data were obtained with deuterated analogs DPPC_d62, DPPC_d13, and DPPC_d9. For the analysis, we developed a model that is based on volume probability distributions and their spatial conservation. The models design was guided and tested by a DOPC molecular dynamics simulation. The model consistently captures the salient features found in both electron and neutron scattering density profiles. A key result of the analysis is the molecular surface area, A. For DPPC at 50 degrees C A = 63.0 A(2), whereas for DOPC at 30 degrees C A = 67.4 A(2), with estimated uncertainties of 1 A(2). Although A for DPPC agrees with a recently reported value obtained solely from the analysis of x-ray scattering data, A for DOPC is almost 10% smaller. This improved method for determining lipid areas helps to reconcile long-standing differences in the values of lipid areas obtained from stand-alone x-ray and neutron scattering experiments and poses new challenges for molecular dynamics simulations.

Amphipols From A to Z

Amphipols (APols) are short amphipathic polymers that can substitute for detergents to keep integral membrane proteins (MPs) water soluble. In this review, we discuss their structure and solution behavior; the way they associate with MPs; and the structure, dynamics, and solution properties of the resulting complexes. All MPs tested to date form water-soluble complexes with APols, and their biochemical stability is in general greatly improved compared with MPs in detergent solutions. The functionality and ligand-binding properties of APol-trapped MPs are reviewed, and the mechanisms by which APols stabilize MPs are discussed. Applications of APols include MP folding and cell-free synthesis, structural studies by NMR, electron microscopy and X-ray diffraction, APol-mediated immobilization of MPs onto solid supports, proteomics, delivery of MPs to preexisting membranes, and vaccine formulation.

Atomistic simulations of biologically realistic transmembrane potential gradients.

We present all-atom molecular dynamics simulations of biologically realistic transmembrane potential gradients across a DMPC bilayer. These simulations are the first to model this gradient in all-atom detail, with the field generated solely by explicit ion dynamics. Unlike traditional bilayer simulations that have one bilayer per unit cell, we simulate a 170 mV potential gradient by using a unit cell consisting of three salt-water baths separated by two bilayers, with full three-dimensional periodicity. The study shows that current computational resources are powerful enough to generate a truly electrified interface, as we show the predicted effect of the field on the overall charge distribution. Additionally, starting from Poissons equation, we show a new derivation of the double integral equation for calculating the potential profile in systems with this type of periodicity.

The Methionine-aromatic Motif Plays a Unique Role in Stabilizing Protein Structure

Background: The interaction between methionine and aromatic residues in protein complexes is poorly understood. Results: The Met-aromatic motif is prevalent in known protein structures and stabilizes TNF ligand-receptor binding interactions. Conclusion: The Met sulfur-aromatic binding motif provides additional stabilization over purely hydrophobic interactions and at longer distances. Significance: This motif is prevalent and may be associated with a number of mutation- and age-associated diseases. Of the 20 amino acids, the precise function of methionine (Met) remains among the least well understood. To establish a determining characteristic of methionine that fundamentally differentiates it from purely hydrophobic residues, we have used in vitro cellular experiments, molecular simulations, quantum calculations, and a bioinformatics screen of the Protein Data Bank. We show that approximately one-third of all known protein structures contain an energetically stabilizing Met-aromatic motif and, remarkably, that greater than 10,000 structures contain this motif more than 10 times. Critically, we show that as compared with a purely hydrophobic interaction, the Met-aromatic motif yields an additional stabilization of 1–1.5 kcal/mol. To highlight its importance and to dissect the energetic underpinnings of this motif, we have studied two clinically relevant TNF ligand-receptor complexes, namely TRAIL-DR5 and LTα-TNFR1. In both cases, we show that the motif is necessary for high affinity ligand binding as well as function. Additionally, we highlight previously overlooked instances of the motif in several disease-related Met mutations. Our results strongly suggest that the Met-aromatic motif should be exploited in the rational design of therapeutics targeting a range of proteins.

The Effect of Cholesterol on Short- and Long-Chain Monounsaturated Lipid Bilayers as Determined by Molecular Dynamics Simulations and X-Ray Scattering

We investigate the structure of cholesterol-containing membranes composed of either short-chain (diC14:1PC) or long-chain (diC22:1PC) monounsaturated phospholipids. Bilayer structural information is derived from all-atom molecular dynamics simulations, which are validated via direct comparison to x-ray scattering experiments. We show that the addition of 40 mol % cholesterol results in a nearly identical increase in the thickness of the two different bilayers. In both cases, the chain ordering dominates over the hydrophobic matching between the length of the cholesterol molecule and the hydrocarbon thickness of the bilayer, which one would expect to cause a thinning of the diC22:1PC bilayer. For both bilayers there is substantial headgroup rearrangement for lipids directly in contact with cholesterol, supporting the so-called umbrella model. Importantly, in diC14:1PC bilayers, a dynamic network of hydrogen bonds stabilizes long-lived reorientations of some cholesterol molecules, during which they are found to lie perpendicular to the bilayer normal, deep within the bilayers hydrophobic core. Additionally, the simulations show that the diC14:1PC bilayer is significantly more permeable to water. These differences may be correlated with faster cholesterol flip-flop between the leaflets of short-chain lipid bilayers, resulting in an asymmetric distribution of cholesterol molecules. This asymmetry was observed experimentally in a case of unilamellar vesicles (ULVs), and reproduced through a set of novel asymmetric simulations. In contrast to ULVs, experimental data for oriented multilamellar stacks does not show the asymmetry, suggesting that it results from the curvature of the ULV bilayers.

Interleaflet Interaction and Asymmetry in Phase Separated Lipid Bilayers: Molecular Dynamics Simulations

In order to investigate experimentally inaccessible, molecular-level detail regarding interleaflet interaction in membranes, we have run an extensive series of coarse-grained molecular dynamics simulations of phase separated lipid bilayers. The simulations are motivated by differences in lipid and cholesterol composition in the inner and outer leaflets of biological membranes. Over the past several years, this phenomenon has inspired a series of experiments in model membrane systems which have explored the effects of lipid compositional asymmetry in the two leaflets. The simulations are directed at understanding one potential consequence of compositional asymmetry, that being regions of bilayers where liquid-ordered (L(o)) domains in one leaflet are opposite liquid-disordered (L(d)) domains in the other leaflet (phase asymmetry). The simulated bilayers are of two sorts: 1) Compositionally symmetric leaflets where each of the two leaflets contains an identical, phase separated (L(o)/L(d)) mixture of cholesterol, saturated and unsaturated phospholipid; and 2) Compositionally asymmetric leaflets, where one leaflet contains a phase separated (L(o)/L(d)) mixture while the other contains only unsaturated lipid, which on its own would be in the L(d) phase. In addition, we have run simulations where the lengths of the saturated lipid chains as well as the mole ratios of the three lipid components are varied. Collectively, we report on three types of interleaflet coupling within a bilayer. First, we show the effects of compositional asymmetry on acyl chain tilt and order, lipid rotational dynamics, and lateral diffusion in regions of leaflets that are opposite L(o) domains. Second, we show substantial effects of compositional asymmetry on local bilayer curvature, with the conclusion that phase separated leaflets resist curvature, while inducing large degrees of curvature in an opposing L(d) leaflet. Finally, in compositionally symmetric, phase separated bilayers, we find phase asymmetry (domain antiregistration) between the two leaflets occurs as a consequence of mismatched acyl chain-lengths in the saturated and unsaturated lipids.

Interpretation of Fluctuation Spectra in Lipid Bilayer Simulations

Atomic resolution and coarse-grained simulations of dimyristoylphosphatidylcholine lipid bilayers were analyzed for fluctuations perpendicular to the bilayer using a completely Fourier-based method. We find that the fluctuation spectrum of motions perpendicular to the bilayer can be decomposed into just two parts: 1), a pure undulation spectrum proportional to q(-4) that dominates in the small-q regime; and 2), a molecular density structure factor contribution that dominates in the large-q regime. There is no need for a term proportional to q(-2) that has been postulated for protrusion fluctuations and that appeared to have been necessary to fit the spectrum for intermediate q. We suggest that earlier reports of such a term were due to the artifact of binning and smoothing in real space before obtaining the Fourier spectrum. The observability of an intermediate protrusion regime from the fluctuation spectrum is discussed based on measured and calculated material constants.

Curvature dynamics of α-Synuclein familial Parkinson disease mutants. Molecular simulations of the Micelle-and Bilayer-bound forms

α-Synuclein remains a protein of interest due to its propensity to form fibrillar aggregates in neurodegenerative disease and its putative function in synaptic vesicle regulation. Herein, we present a series of atomistic molecular dynamics simulations of wild-type α-synuclein and three Parkinson disease familial mutants (A30P, A53T, and E46K) in two distinct environments. First, in order to match recent NMR experiments, we have simulated each protein bound to an SDS detergent micelle. Second, in order to connect more closely to the true biological environment, we have simulated the proteins bound to a 1,2-dioleoyl-sn-glycero-3-phosphoserine lipid bilayer. In the micelle-bound case, we find that the wild type and all of the variants of α-synuclein flatten the underlying micelle, decreasing its surface area. A30P is known to lessen α-synuclein/membrane affinity and, consistent with experiment, destabilizes the simulated secondary structure. In the case of A53T, our simulations reveal a range of stabilizing hydrogen bonds that form with the threonine. In both environments, the E46K mutation, which is known to increase bilayer affinity, leads to an additional hydrogen bond between the protein and either the detergent or lipid. Simulations indicate that αS and its variants are less dynamic in the bilayer than in the micelle. Furthermore, the simulations of the mutants suggest how changes in the structure and dynamics of α-synuclein may affect its biological role.Alpha-synuclein remains a protein of interest due to its propensity to form fibrillar aggregates in neurodegenerative disease and its putative function in synaptic vesicle regulation. Herein, we present a series of atomistic molecular dynamics simulations of wild-type alpha-synuclein and three Parkinson disease familial mutants (A30P, A53T, and E46K) in two distinct environments. First, in order to match recent NMR experiments, we have simulated each protein bound to an SDS detergent micelle. Second, in order to connect more closely to the true biological environment, we have simulated the proteins bound to a 1,2-dioleoyl-sn-glycero-3-phosphoserine lipid bilayer. In the micelle-bound case, we find that the wild type and all of the variants of alpha-synuclein flatten the underlying micelle, decreasing its surface area. A30P is known to lessen alpha-synuclein/membrane affinity and, consistent with experiment, destabilizes the simulated secondary structure. In the case of A53T, our simulations reveal a range of stabilizing hydrogen bonds that form with the threonine. In both environments, the E46K mutation, which is known to increase bilayer affinity, leads to an additional hydrogen bond between the protein and either the detergent or lipid. Simulations indicate that alphaS and its variants are less dynamic in the bilayer than in the micelle. Furthermore, the simulations of the mutants suggest how changes in the structure and dynamics of alpha-synuclein may affect its biological role.

α-Synuclein induces both positive mean curvature and negative Gaussian curvature in membranes.

Using a combination of X-ray scattering, fluorescence correlation spectroscopy, coarse-grained molecular dynamics (MD) simulations and potential of mean force calculations, we have explored the membrane remodeling effects of monomeric α-synuclein (αS). Our initial findings from multiple approaches are that αS (1) causes a significant thinning of the bilayer and (2) stabilizes positive mean curvature, such that the maximum principle curvature matches that of synaptic vesicles, αS-induced tubules, and the synthetic lipid vesicles to which the protein binds most tightly. This suggests that αS binding to synaptic vesicles likely stabilizes their intrinsic curvature. We then show that αS induces local negative Gaussian curvature, an effect that occurs in regions of αS shown previously via NMR and corroborated by MD simulation to have significant conformational flexibility. The induction of negative Gaussian curvature, which has implications for all curvature-sensing and curvature-generating amphipathic α-helices, supports a hypothesis that connects helix insertion to fusion and fission of vesicles, processes that have recently been linked to αS function. Then, in an effort to explain these biophysical properties of αS, we promote an intrinsic curvature-field model that recasts long-range protein-protein interactions in terms of the interactions between the local curvature fields generated by lipid-protein complexes.

Determination of Electron Density Profiles and Area from Simulations of Undulating Membranes

The traditional method for extracting electron density and other transmembrane profiles from molecular dynamics simulations of lipid bilayers fails for large bilayer systems, because it assumes a flat reference surface that does not take into account long wavelength undulations. We have developed what we believe to be a novel set of methods to characterize these undulations and extract the underlying profiles in the large systems. Our approach first obtains an undulation reference surface for each frame in the simulation and subsequently isolates the long-wavelength undulations by filtering out the intrinsic short wavelength modes. We then describe two methods to obtain the appropriate profiles from the undulating reference surface. Most combinations of methods give similar results for the electron density profiles of our simulations of 1024 DMPC lipids. From simulations of smaller systems, we also characterize the finite size effect related to the boundary conditions of the simulation box. In addition, we have developed a set of methods that use the undulation reference surface to determine the true area per lipid which, due to undulations, is larger than the projected area commonly reported from simulations.