Jason D. Perlmutter

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.

Mechanisms of Virus Assembly

Viruses are nanoscale entities containing a nucleic acid genome encased in a protein shell called a capsid and in some cases are surrounded by a lipid bilayer membrane. This review summarizes the physics that govern the processes by which capsids assemble within their host cells and in vitro. We describe the thermodynamics and kinetics for the assembly of protein subunits into icosahedral capsid shells and how these are modified in cases in which the capsid assembles around a nucleic acid or on a lipid bilayer. We present experimental and theoretical techniques used to characterize capsid assembly, and we highlight aspects of virus assembly that are likely to receive significant attention in the near future.

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.

Viral genome structures are optimal for capsid assembly

Understanding how virus capsids assemble around their nucleic acid (NA) genomes could promote efforts to block viral propagation or to reengineer capsids for gene therapy applications. We develop a coarse-grained model of capsid proteins and NAs with which we investigate assembly dynamics and thermodynamics. In contrast to recent theoretical models, we find that capsids spontaneously ‘overcharge’; that is, the negative charge of the NA exceeds the positive charge on capsid. When applied to specific viruses, the optimal NA lengths closely correspond to the natural genome lengths. Calculations based on linear polyelectrolytes rather than base-paired NAs underpredict the optimal length, demonstrating the importance of NA structure to capsid assembly. These results suggest that electrostatics, excluded volume, and NA tertiary structure are sufficient to predict assembly thermodynamics and that the ability of viruses to selectively encapsidate their genomic NAs can be explained, at least in part, on a thermodynamic basis. DOI: http://dx.doi.org/10.7554/eLife.00632.001

All-atom and coarse-grained molecular dynamics simulations of a membrane protein stabilizing polymer.

Amphipathic polymers called amphipols (APols) have been developed as an alternative to detergents for stabilizing membrane proteins (MPs) in aqueous solutions. APols provide MPs with a particularly mild environment and, as a rule, keep them in a native functional state for longer periods than do detergents. Amphipol A8-35, a derivative of polyacrylate, is widely used and has been particularly well studied experimentally. In aqueous solutions, A8-35 molecules self-assemble into well-defined globular particles with a mass of ∼40 kDa and a R(g) of ∼2.4 nm. As a first step towards describing MP/A8-35 complexes by molecular dynamics (MD), we present three sets of simulations of the pure APol particle. First, we performed a series of all-atom MD (AAMD) simulations of the particle in solution, starting from an arbitrary initial configuration. Although AAMD simulations result in stable cohesive particles over a 45 ns simulation, the equilibration of the particle organization is limited. This motivated the use of coarse-grained MD (CGMD), allowing us to investigate processes on the microsecond time scale, including de novo particle assembly. We present a detailed description of the parametrization of the CGMD model from the AAMD simulations and a characterization of the resulting CGMD particles. Our third set of simulations utilizes reverse coarse-graining (rCG), through which we obtain all-atom coordinates from a CGMD simulation. This allows a higher-resolution characterization of a configuration determined by a long-timescale simulation. Excellent agreement is observed between MD models and experimental, small-angle neutron scattering data. The MD data provides new insight into the structure and dynamics of A8-35 particles, which is possibly relevant to the stabilizing effects of APols on MPs, as well as a starting point for modeling MP/A8-35 complexes.

Pathways for Virus Assembly around Nucleic Acids

Understanding the pathways by which viral capsid proteins assemble around their genomes could identify key intermediates as potential drug targets. In this work, we use computer simulations to characterize assembly over a wide range of capsid protein-protein interaction strengths and solution ionic strengths. We find that assembly pathways can be categorized into two classes, in which intermediates are either predominantly ordered or disordered. Our results suggest that estimating the protein-protein and the protein-genome binding affinities may be sufficient to predict which pathway occurs. Furthermore, the calculated phase diagrams suggest that knowledge of the dominant assembly pathway and its relationship to control parameters could identify optimal strategies to thwart or redirect assembly to block infection. Finally, analysis of simulation trajectories suggests that the two classes of assembly pathways can be distinguished in single-molecule fluorescence correlation spectroscopy or bulk time-resolved small-angle X-ray scattering experiments.

Tumor Necrosis Factor-related Apoptosis-inducing Ligand (TRAIL) Induces Death Receptor 5 Networks That Are Highly Organized

Background: Whether ligand-induced clusters of DR5 have a specific structural organization is unknown. Results: Ligand binding results in the formation of death receptor dimers that exist within high molecular weight networks. Conclusion: Ligand-induced DR5 clusters are highly organized networks formed through dimerization of receptor trimers. Significance: The biophysical character of DR5 networks may have implications for future rational design of DR5-targeted therapeutics. Recent evidence suggests that TNF-related apoptosis-inducing ligand (TRAIL), a death-inducing cytokine with anti-tumor potential, initiates apoptosis by re-organizing TRAIL receptors into large clusters, although the structure of these clusters and the mechanism by which they assemble are unknown. Here, we demonstrate that TRAIL receptor 2 (DR5) forms receptor dimers in a ligand-dependent manner at endogenous receptor levels, and these receptor dimers exist within high molecular weight networks. Using mutational analysis, FRET, fluorescence microscopy, synthetic biochemistry, and molecular modeling, we find that receptor dimerization relies upon covalent and noncovalent interactions between membrane-proximal residues. Additionally, by using FRET, we show that the oligomeric structure of two functional isoforms of DR5 is indistinguishable. The resulting model of DR5 activation should revise the accepted architecture of the functioning units of DR5 and the structurally homologous TNF receptor superfamily members.

Molecular Dynamics Simulations of a Membrane Protein/Amphipol Complex

Amphipathic polymers known as “amphipols” provide a highly stabilizing environment for handling membrane proteins in aqueous solutions. A8-35, an amphipol with a polyacrylate backbone and hydrophobic grafts, has been extensively characterized and widely employed for structural and functional studies of membrane proteins using biochemical and biophysical approaches. Given the sensitivity of membrane proteins to their environment, it is important to examine what effects amphipols may have on the structure and dynamics of the proteins they complex. Here we present the first molecular dynamics study of an amphipol-stabilized membrane protein, using Escherichia coli OmpX as a model. We begin by describing the structure of the complexes formed by supplementing OmpX with increasing amounts of A8-35, in order to determine how the amphipol interacts with the transmembrane and extramembrane surfaces of the protein. We then compare the dynamics of the protein in either A8-35, a detergent, or a lipid bilayer. We find that protein dynamics on all accessible length scales is restrained by A8-35, which provides a basis to understanding some of the stabilizing and functional effects of amphipols that have been experimentally observed.