Assembly of lipids and proteins into lipoprotein particles

The structure, self-assembly and dynamics of lipid nanodiscs revealed by computational approaches

Nanodisc technology is increasingly being used in structural, biochemical and biophysical studies of membrane proteins. The computational approaches have revealed many important features of nanodisc assembly, structures and dynamics. Therefore, we reviewed the application of computational approaches, especially molecular modeling and molecular dyncamics (MD) simulations, to characterize nanodiscs, including the structural models, assembly and disassembly, protocols for modeling, structural properties and dynamics, and protein-lipid interactions in nanodiscs. More amazing computational studies about nanodiscs are looked forward to in the future.

Styrene-Maleic Acid Copolymer Nanodiscs to Determine the Shape of Membrane Proteins

Lipid nanodiscs can be used to solubilize functional membrane proteins (MPs) in nativelike environments. Thus, they are promising reagents that have been proven useful to characterize MPs. Both protein and non-protein molecular belts have shown promise to maintain the structural integrity of MPs in lipid nanodiscs. Small-angle neutron scattering (SANS) can be used to determine low-resolution structures of proteins in solution, which can be enhanced through the use of contrast variation methods. We present theoretical contrast variation SANS results for protein and styrene-maleic acid copolymer (SMA) belt 1,2-dimyristoyl-sn-glycero-3-phosphorylcholine (DMPC) nanodiscs with and without additional bound or transmembrane proteins. The predicted scattering properties are derived from atomistic molecular dynamics simulations to account for conformational fluctuations, and we determine deuterium-labeling conditions such that SANS intensity profiles only include contributions from the scattering of the MP of interest. We propose strategies to tune the neutron scattering length densities (SLDs) of the SMA and DMPC using selective deuterium labeling such that the SLD of the nanodisc becomes homogeneous and its scattering can essentially be eliminated in solvents containing an appropriate amount of D₂O. These finely tuned labeled polymer-based nanodiscs are expected to be useful to extract the size and molecular shape information of MPs using SANS-based contrast variation experiments, and they can be used with MPs of any molecular weight.

Nanodiscs: A toolkit for membrane protein science

Membrane proteins are involved in numerous vital biological processes, including transport, signal transduction and the enzymes in a variety of metabolic pathways. Integral membrane proteins account for up to 30% of the human proteome and they make up more than half of all currently marketed therapeutic targets. Unfortunately, membrane proteins are inherently recalcitrant to study using the normal toolkit available to scientists, and one is most often left with the challenge of finding inhibitors, activators and specific antibodies using a denatured or detergent solubilized aggregate. The Nanodisc platform circumvents these challenges by providing a self-assembled system that renders typically insoluble, yet biologically and pharmacologically significant, targets such as receptors, transporters, enzymes, and viral antigens soluble in aqueous media in a native-like bilayer environment that maintain a target's functional activity. By providing a bilayer surface of defined composition and structure, Nanodiscs have found great utility in the study of cellular signaling complexes that assemble on a membrane surface. Nanodiscs provide a nanometer scale vehicle for the in vivo delivery of amphipathic drugs, therapeutic lipids, tethered nucleic acids, imaging agents and active protein complexes. This means for generating nanoscale lipid bilayers has spawned the successful use of numerous other polymer and peptide amphipathic systems. This review, in celebration of the Anfinsen Award, summarizes some recent results and provides an inroad into the current and historical literature.

Dark, Ultra-Dark and Ultra-Bright Nanodiscs for membrane protein investigations

We describe the construction, expression and purification of three new membrane scaffold proteins (MSP) for use in assembling Nanodiscs. These new MSPs have a variety of luminescent properties for use in combination with several analytical methods. "Dark" MSP has no tryptophan residues, "Ultra-Dark" replaces both tryptophan and tyrosine with non-fluorescent side chains, and "Ultra-Bright" adds additional tryptophans to the parent membrane scaffold protein to provide a dramatic increase in native tryptophan fluorescence. All MSPs were used to successfully assemble Nanodiscs nominally 10 nm in diameter, and the resultant bilayer structure was characterized. An example of the usefulness of these new scaffold proteins is provided.

Structure and dynamics of a nanodisc by integrating NMR, SAXS and SANS experiments with molecular dynamics simulations

Nanodiscs are membrane mimetics that consist of a protein belt surrounding a lipid bilayer, and are broadly used for characterization of membrane proteins. Here, we investigate the structure, dynamics and biophysical properties of two small nanodiscs, MSP1D1ΔH5 and ΔH4H5. We combine our SAXS and SANS experiments with molecular dynamics simulations and previously obtained NMR and EPR data to derive and validate a conformational ensemble that represents the structure and dynamics of the nanodisc. We find that it displays conformational heterogeneity with various elliptical shapes, and with substantial differences in lipid ordering in the centre and rim of the discs. Together, our results reconcile previous apparently conflicting observations about the shape of nanodiscs, and pave the way for future integrative studies of larger complex systems such as membrane proteins embedded in nanodiscs.

Functional Reconstitution of HlyB, a Type I Secretion ABC Transporter, in Saposin-A Nanoparticles

Type I secretion systems (T1SS) are ubiquitous transport machineries in Gram-negative bacteria. They comprise a relatively simple assembly of three membrane-localised proteins: an inner-membrane complex composed of an ABC transporter and a membrane fusion protein, and a TolC-like outer membrane component. T1SS transport a wide variety of substrates with broad functional diversity. The ABC transporter hemolysin B (HlyB), for example, is part of the hemolysin A-T1SS in Escherichia coli. In contrast to canonical ABC transporters, an accessory domain, a C39 peptidase-like domain (CLD), is located at the N-terminus of HlyB and is essential for secretion. In this study, we have established an optimised purification protocol for HlyB and the subsequent reconstitution employing the saposin-nanoparticle system. We point out the negative influence of free detergent on the basal ATPase activity of HlyB, studied the influence of a lysolipid or lipid matrix on activity and present functional studies with the full-length substrate proHlyA in its folded and unfolded states, which both have a stimulatory effect on the ATPase activity.

CHARMM-GUI Nanodisc Builder for modeling and simulation of various nanodisc systems

Nanodiscs are discoidal protein-lipid complexes that have wide applications in membrane protein studies. Modeling and simulation of nanodiscs are challenging due to the absence of structures of many membrane scaffold proteins (MSPs) that wrap around the membrane bilayer. We have developed CHARMM-GUI Nanodisc Builder (http://www.charmm-gui.org/input/nanodisc) to facilitate the setup of nanodisc simulation systems by modeling the MSPs with defined size and known structural features. A total of 11 different nanodiscs with a diameter from 80 to 180 Å are made available in both the all-atom CHARMM and two coarse-grained (PACE and Martini) force fields. The usage of the Nanodisc Builder is demonstrated with various simulation systems. The structures and dynamics of proteins and lipids in these systems were analyzed, showing similar behaviors to those from previous all-atom and coarse-grained nanodisc simulations. We expect the Nanodisc Builder to be a convenient and reliable tool for modeling and simulation of nanodisc systems. © 2019 Wiley Periodicals, Inc.

Molecular simulations of self-assembly processes in metal-organic frameworks: Model dependence

Molecular simulation is a powerful tool for investigating microscopic behavior in various chemical systems, where the use of suitable models is critical to successfully reproduce the structural and dynamic properties of the real systems of interest. In this context, molecular dynamics simulation studies of self-assembly processes in metal-organic frameworks (MOFs), a well-known class of porous materials with interesting chemical and physical properties, are relatively challenging, where a reasonably accurate representation of metal-ligand interactions is anticipated to play an important role. In the current study, we both investigate the performance of some existing models and introduce and test new models to help explore the self-assembly in an archetypal Zn-carboxylate MOF system. To this end, the behavior of six different Zn-ion models, three solvent models, and two ligand models was examined and validated against key experimental structural parameters. To explore longer time scale ordering events during MOF self-assembly via explicit solvent simulations, it is necessary to identify a suitable combination of simplified model components representing metal ions, organic ligands, and solvent molecules. It was observed that an extended cationic dummy atom (ECDA) Zn-ion model combined with an all-atom carboxylate ligand model and a simple dipolar solvent model can reproduce characteristic experimental structures for the archetypal MOF system. The successful use of these models in extensive sets of molecular simulations, which provide key insights into the self-assembly mechanism of this archetypal MOF system occurring during the early stages of this process, has been very recently reported.

Nanodiscs in Membrane Biochemistry and Biophysics

Membrane proteins play a most important part in metabolism, signaling, cell motility, transport, development, and many other biochemical and biophysical processes which constitute fundamentals of life on the molecular level. Detailed understanding of these processes is necessary for the progress of life sciences and biomedical applications. Nanodiscs provide a new and powerful tool for a broad spectrum of biochemical and biophysical studies of membrane proteins and are commonly acknowledged as an optimal membrane mimetic system that provides control over size, composition, and specific functional modifications on the nanometer scale. In this review we attempted to combine a comprehensive list of various applications of nanodisc technology with systematic analysis of the most attractive features of this system and advantages provided by nanodiscs for structural and mechanistic studies of membrane proteins.

The power, pitfalls and potential of the nanodisc system for NMR-based studies

The choice of a suitable membrane mimicking environment is of fundamental importance for the characterization of structure and function of membrane proteins. In this respect, usage of the lipid bilayer nanodisc technology provides a unique potential for nuclear magnetic resonance (NMR)-based studies. This review summarizes the recent advances in this field, focusing on (i) the strengths of the system, (ii) the bottlenecks that may be faced, and (iii) promising capabilities that may be explored in future studies.

Computational studies of plasma lipoprotein lipids

Plasma lipoproteins are macromolecular assemblies of proteins and lipids found in the blood. The lipid components of lipoproteins are amphipathic lipids such as phospholipids (PLs), and unesterified cholesterols (UCs) and hydrophobic lipids such as cholesteryl esters (CEs) and triglycerides (TGs). Since lipoproteins are soft matter supramolecular assemblies easily deformable by thermal fluctuations and they also exist in varying densities and protein/lipid components, a detailed understanding of their structure/function is experimentally difficult. Molecular dynamics (MD) simulation has emerged as a particularly promising way to explore the structure and dynamics of lipoproteins. The purpose of this review is to survey the current status of computational studies of the lipid components of the lipoproteins. Computational studies aim to explore three levels of complexity for the 3-dimensional structural dynamics of lipoproteins at various metabolic stages: (i) lipoprotein particles consist of protein with minimal lipid; (ii) lipoprotein particles consist of PL-rich discoidal bilayer-like lipid particles; (iii) mature circulating lipoprotein particles consist of CE-rich or TG-rich spheroidal lipid-droplet-like particles. Due to energy barriers involved in conversion between these species, other biomolecules also participate in lipoprotein biological assembly. For example: (i) lipid-poor apolipoprotein A-I (apoA-I) interacts with ATP-binding cassette transporter A1 (ABCA1) to produce nascent discoidal high density lipoprotein (dHDL) particles; (ii) lecithin-cholesterol acyltransferase (LCAT) mediates the conversion of UC to CE in dHDL, driving spheroidal HDL (sHDL) formation; (iii) transfer proteins, cholesterol ester transfer protein (CETP) and phospholipid transfer protein (PLTP), transfer both CE and TG and PL, respectively, between lipoprotein particles. Computational studies have the potential to explore different lipoprotein particles at each metabolic stage in atomistic detail. This review discusses the current status of computational methods including all-atom MD (AAMD), coarse-grain MD (CGMD), and MD-simulated annealing (MDSA) and their applications in lipoprotein structural dynamics and biological assemblies. Results from MD simulations are discussed and compared across studies in order to identify key findings, controversies, issues and future directions. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Encapsulated membrane proteins: A simplified system for molecular simulation

Over the past 50years there has been considerable progress in our understanding of biomolecular interactions at an atomic level. This in turn has allowed molecular simulation methods employing full atomistic modelling at ever larger scales to develop. However, some challenging areas still remain where there is either a lack of atomic resolution structures or where the simulation system is inherently complex. An area where both challenges are present is that of membranes containing membrane proteins. In this review we analyse a new practical approach to membrane protein study that offers a potential new route to high resolution structures and the possibility to simplify simulations. These new approaches collectively recognise that preservation of the interaction between the membrane protein and the lipid bilayer is often essential to maintain structure and function. The new methods preserve these interactions by producing nano-scale disc shaped particles that include bilayer and the chosen protein. Currently two approaches lead in this area: the MSP system that relies on peptides to stabilise the discs, and SMALPs where an amphipathic styrene maleic acid copolymer is used. Both methods greatly enable protein production and hence have the potential to accelerate atomic resolution structure determination as well as providing a simplified format for simulations of membrane protein dynamics. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Structural Insights into High Density Lipoprotein: Old Models and New Facts

The physiological link between circulating high density lipoprotein (HDL) levels and cardiovascular disease is well-documented, albeit its intricacies are not well-understood. An improved appreciation of HDL function and overall role in vascular health and disease requires at its foundation a better understanding of the lipoprotein's molecular structure, its formation, and its process of maturation through interactions with various plasma enzymes and cell receptors that intervene along the pathway of reverse cholesterol transport. This review focuses on summarizing recent developments in the field of lipid free apoA-I and HDL structure, with emphasis on new insights revealed by newly published nascent and spherical HDL models constructed by combining low resolution structures obtained from small angle neutron scattering (SANS) with contrast variation and geometrical constraints derived from hydrogen-deuterium exchange (HDX), crosslinking mass spectrometry, electron microscopy, Förster resonance energy transfer, and electron spin resonance. Recently published low resolution structures of nascent and spherical HDL obtained from SANS with contrast variation and isotopic labeling of apolipoprotein A-I (apoA-I) will be critically reviewed and discussed in terms of how they accommodate existing biophysical structural data from alternative approaches. The new low resolution structures revealed and also provided some answers to long standing questions concerning lipid organization and particle maturation of lipoproteins. The review will discuss the merits of newly proposed SANS based all atom models for nascent and spherical HDL, and compare them with accepted models. Finally, naturally occurring and bioengineered mutations in apoA-I, and their impact on HDL phenotype, are reviewed and discuss together with new therapeutics employed for restoring HDL function.

Force Field Development for Lipid Membrane Simulations

With the rapid development of computer power and wide availability of modelling software computer simulations of realistic models of lipid membranes, including their interactions with various molecular species, polypeptides and membrane proteins have become feasible for many research groups. The crucial issue of the reliability of such simulations is the quality of the force field, and many efforts, especially in the latest several years, have been devoted to parametrization and optimization of the force fields for biomembrane modelling. In this review, we give account of the recent development in this area, covering different classes of force fields, principles of the force field parametrization, comparison of the force fields, and their experimental validation. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Bicelles and Other Membrane Mimics: Comparison of Structure, Properties, and Dynamics from MD Simulations

The increased interest in studying membrane proteins has led to the development of new membrane mimics such as bicelles and nanodiscs. However, only limited knowledge is available of how these membrane mimics are affected by embedded proteins and how well they mimic a lipid bilayer. Herein, we present molecular dynamics simulations to elucidate structural and dynamic properties of small bicelles and compare them to a large alignable bicelle, a small nanodisc, and a lipid bilayer. Properties such as lipid packing and properties related to embedding both an α-helical peptide and a transmembrane protein are investigated. The small bicelles are found to be very dynamic and mainly assume a prolate shape substantiating that small bicelles cannot be regarded as well-defined disclike structures. However, addition of a peptide results in an increased tendency to form disc-shaped bicelles. The small bicelles and the nanodiscs show increased peptide solvation and difference in peptide orientation compared to embedding in a bilayer. The large bicelle imitated a bilayer well with respect to both curvature and peptide solvation, although peripheral binding of short tailed lipids to the embedded proteins is observed, which could hinder ligand binding or multimer formation.

Toward a Coarse-Grained Protein Model Coupled with a Coarse-Grained Solvent Model: Solvation Free Energies of Amino Acid Side Chains

Recently, we reported that molecular dynamics (MD) simulations using a coarse-grained (CG) peptide model coupled with a CG water model are able to reproduce many of the structural and thermodynamic features of short peptides with nonpolar side chains at 10(3) times the normal speed (JCTC, 2007, 3, 2146-2161). To further develop a CG protein model for MD simulations, we systematically parametrized the side chains of all 20 naturally occurring amino acids. We developed the parameters by fitting the dihedral potentials of 13 small molecules, the densities and self-solvation free energies of liquids of eight organic molecules, and the hydration free energies of 35 small organic molecules. In a set of 11 classes of compounds (105 in total) including alkanes, alcohols, ethers, ketones/aldehydes, amines, amides, aromatics, carboxylic acids, sulfides/thiols, alkyl ammoniums, and carboxylate ions, the average error in the calculated hydration free energies compared with experimental results is about 1.4 kJ/mol. The average error in the calculated transfer free energies of the 19 side-chain analogues of amino acids from cyclohexane to water is about 2.2 kJ/mol. These results are comparable to the results of all-atom models.

Structure and Dynamics of Phospholipid Nanodiscs from All-Atom and Coarse-Grained Simulations

We investigated structural and dynamical properties of nanodiscs comprising dimyristoylphosphatidylcholine (DMPC) lipids and major scaffold protein MSP1Δ(1-22) from human apolipoprotein A-1 using combined all-atom and coarse-grained (CG) molecular dynamics (MD) simulations. The computational efficiency of the Martini-CG force field enables the spontaneous self-assembly of lipids and scaffold proteins into stable nanodisc structures on time scales up to tens of microseconds. Subsequent all-atom and CG-MD simulations reveal that the lipids in the nanodisc have lower configurational entropy and higher acyl tail order than in a lamellar bilayer phase. These altered average properties arise from rather differential behavior of lipids, depending on their location in the nanodisc. Since the scaffold proteins exert constrictive forces from the outer rim of the disc toward its center, lipids at the center of the nanodisc are highly ordered, whereas annular lipids that are in contact with the MSP proteins are remarkably disordered due to perturbed packing. Although specific differences between all-atom and CG simulations are also evident, the results obtained at both levels of resolution are in overall good agreement with each other and provide atomic level interpretations of recent experiments. Thus, the present study highlights the applicability of multiscale simulation approaches for nanodisc systems and opens the way for future applications, including the study of nanodisc-embedded membrane proteins.

Reconstitution of membrane proteins: a GPCR as an example

Membrane proteins are the gatekeepers to the cell and are essential to the function of all cells, controlling the flow of molecules and information across the cell membrane. Much effort has been put into the development of systems for studying membrane proteins in simplified environments that nevertheless mimic their native lipid environment. After isolation and production of purified membrane proteins in detergent, it is often necessary to reconstitute them into a lipid structure such as liposome, nanodisc, or lipodisq. Each of these has the advantage of returning the protein to a defined lipid environment, and the choice of system depends on the application. Regardless of the system to be used, the fundamental process involves the removal of detergent and incorporation of the protein into a stable lipid system. This chapter details methodologies we have developed, mainly focussed on the model G protein-coupled receptor (GPCR) neurotensin receptor 1, and the GPCR-homologue and model, bacteriorhopdopsin.

Tunable nucleation time of functional sphingomyelinase--lipid features studied by membrane array statistic tool

Aggregation or assembly of lipids and proteins could significantly change the proteins' function. A peripheral membrane enzyme, sphingomyelinase (SMase), has been reported to be able to assemble to a functional feature with its lipid substrate, sphingomyelin (SM), and its lipid product, ceramide (Cer). SMase seems to processes its substrate more effectively in this feature. Here, we report that the functional feature has a tunable formation time. The peculiar behavior is that the feature formation has a time lag depending on the membrane composition. We hypothesized that the time lag is due to the significant nucleation energy barrier when the feature phase forms in its metastable parent phase in the 2-D lipid membrane. To study the stochastic nucleation of the feature, we built a corralled lipid membrane platform with numerous isolated membrane systems in parallel to capture the nucleation statistics. Using the high-throughput approach and the appropriate experimental design to circumvent the interplay of the complicated phase segregation in membranes induced by SMase, we found that the nucleation rate of the feature can be tuned by the supersaturation of the enzyme, the lipid substrate, and the lipid product, in the fluid phase of the membrane. The correlation between the supersaturation and the nucleation rate can be well described by the classical nucleation theory equation, suggesting that the feature formation follows the nucleation process with a certain component ratio specified in the equation. The certain relative component ratio suggests that the feature may have certain organization instead of being random aggregation. In addition, our finding suggests that nucleation could serve as a time lag control mechanism in this enzymatic system, and ways to reduce nucleation energy barrier could be used to shorten the aggregation time lag and vice versa.

Multiscale molecular dynamics simulations of membrane proteins

The time and length scales accessible by biomolecular simulations continue to increase. This is in part due to improvements in algorithms and computing performance, but is also the result of the emergence of coarse-grained (CG) potentials, which complement and extend the information obtainable from fully detailed models. CG methods have already proven successful for a range of applications that benefit from the ability to rapidly simulate spontaneous self-assembly within a lipid membrane environment, including the insertion and/or oligomerization of a range of "toy models," transmembrane peptides, and single- and multi-domain proteins. While these simplified approaches sacrifice atomistic level detail, it is now straightforward to "reverse map" from CG to atomistic descriptions, providing a strategy to assemble membrane proteins within a lipid environment, prior to all-atom simulation. Moreover, recent developments have been made in "dual resolution" techniques, allowing different molecules in the system to be modeled with atomistic or CG resolution simultaneously.

The future of molecular dynamics simulations in drug discovery

Molecular dynamics simulations can now track rapid processes—those occurring in less than about a millisecond—at atomic resolution for many biologically relevant systems. These simulations appear poised to exert a significant impact on how new drugs are found, perhaps even transforming the very process of drug discovery. We predict here future results we can expect from, and enhancements we need to make in, molecular dynamics simulations over the coming 25 years, and in so doing set out several Grand Challenges for the field. In the context of the problems now facing the pharmaceutical industry, we ask how we can best address drug discovery needs of the next quarter century using molecular dynamics simulations, and we suggest some possible approaches.

Conformational transitions in the membrane scaffold protein of phospholipid bilayer nanodiscs

Phospholipid bilayer nanodiscs are model membrane systems that provide an environment where membrane proteins are highly stable and monodisperse without the use of detergents or liposomes. Nanodiscs consist of a discoidal phospholipid bilayer encircled by two copies of an amphipathic alpha helical membrane scaffold protein, which is modeled from apolipoprotein A-1. Hydrogen exchange mass spectrometry was used to probe the structure and dynamics of the scaffold protein in the presence and absence of lipid. On nanodisc self-assembly, the entire scaffold protein gained significant protection from exchange, consistent with a large, protein-wide, structural rearrangement. This protection was short-lived and the scaffold protein was highly deuterated within 2 h. Several regions of the scaffold protein, in both the lipid-free and lipid-associated states, displayed EX1 unfolding kinetics. The rapid deuteration of the scaffold protein and the presence of correlated unfolding events both indicate that nanodiscs are dynamic rather than rigid bodies in solution. This work provides a catalog of the expected scaffold protein peptic peptides in a nanodisc-hydrogen exchange mass spectrometry experiment and their deuterium uptake signatures, data that can be used as a benchmark to verify correct assembly and nanodisc structure. Such reference data will be useful control data for all hydrogen exchange mass spectrometry experiments involving nanodiscs in which transmembrane or lipid-associated proteins are the primary molecule(s) of interest.

Formation of salt bridges mediates internal dimerization of myosin VI medial tail domain

The unconventional motor protein, myosin VI, is known to dimerize upon cargo binding to its C-terminal end. It has been shown that one of its tail domains, called the medial tail domain, is a dimerization region. The domain contains an unusual pattern of alternating charged residues and a few hydrophobic residues. To reveal the unknown dimerization mechanism of the medial tail domain, we employed molecular dynamics and single-molecule experimental techniques. Both techniques suggest that the formation of electrostatic-based interhelical salt bridges between oppositely charged residues is a key dimerization factor. For the dimerization to occur, the two identical helices within the dimer do not bind in a symmetric fashion, but rather with an offset of about one helical repeat. Calculations of the dimer-dissociation energy find the contribution of hydrophobic residues to the dimerization process to be minor; they also find that the asymmetric homodimer state is energetically favorable over a state of separate helices.

Role of lipids in spheroidal high density lipoproteins

We study the structure and dynamics of spherical high density lipoprotein (HDL) particles through coarse-grained multi-microsecond molecular dynamics simulations. We simulate both a lipid droplet without the apolipoprotein A-I (apoA-I) and the full HDL particle including two apoA-I molecules surrounding the lipid compartment. The present models are the first ones among computational studies where the size and lipid composition of HDL are realistic, corresponding to human serum HDL. We focus on the role of lipids in HDL structure and dynamics. Particular attention is paid to the assembly of lipids and the influence of lipid-protein interactions on HDL properties. We find that the properties of lipids depend significantly on their location in the particle (core, intermediate region, surface). Unlike the hydrophobic core, the intermediate and surface regions are characterized by prominent conformational lipid order. Yet, not only the conformations but also the dynamics of lipids are found to be distinctly different in the different regions of HDL, highlighting the importance of dynamics in considering the functionalization of HDL. The structure of the lipid droplet close to the HDL-water interface is altered by the presence of apoA-Is, with most prominent changes being observed for cholesterol and polar lipids. For cholesterol, slow trafficking between the surface layer and the regimes underneath is observed. The lipid-protein interactions are strongest for cholesterol, in particular its interaction with hydrophobic residues of apoA-I. Our results reveal that not only hydrophobicity but also conformational entropy of the molecules are the driving forces in the formation of HDL structure. The results provide the first detailed structural model for HDL and its dynamics with and without apoA-I, and indicate how the interplay and competition between entropy and detailed interactions may be used in nanoparticle and drug design through self-assembly.

Conformational analysis of membrane proteins in phospholipid bilayer nanodiscs by hydrogen exchange mass spectrometry

The study of membrane protein structure and enzymology has traditionally been hampered by the inherent insolubility of membrane proteins in aqueous environments and experimental challenges in emulating an in vivo lipid environment. Phospholipid bilayer nanodiscs have recently been shown to be of great use for the study of membrane proteins since they offer a controllable, stable, and monodisperse model membrane with a nativelike lipid bilayer. Here we report the integration of nanodiscs with hydrogen exchange (HX) mass spectrometry (MS) experiments, thereby allowing for analysis of the native conformation of membrane proteins. gamma-Glutamyl carboxylase (GGCX), an approximately 94 kDa transmembrane protein, was inserted into nanodiscs and labeled with deuterium oxide under native conditions. Analytical parameters including sample-handling and chromatographic separation were optimized to measure the incorporation of deuterium into GGCX. Coupling nanodisc technology with HX MS offers an effective approach for investigating the conformation and dynamics of membrane proteins in their native environment and is therefore capable of providing much needed insight into the function of membrane proteins.

Congruency between biophysical data from multiple platforms and molecular dynamics simulation of the double-super helix model of nascent high-density lipoprotein

The predicted structure and molecular trajectories from >80 ns molecular dynamics simulation of the solvated Double-Super Helix (DSH) model of nascent high-density lipoprotein (HDL) were determined and compared with experimental data on reconstituted nascent HDL obtained from multiple biophysical platforms, including small angle neutron scattering (SANS) with contrast variation, hydrogen-deuterium exchange tandem mass spectrometry (H/D-MS/MS), nuclear magnetic resonance spectroscopy (NMR), cross-linking tandem mass spectrometry (MS/MS), fluorescence resonance energy transfer (FRET), electron spin resonance spectroscopy (ESR), and electron microscopy. In general, biophysical constraints experimentally derived from the multiple platforms agree with the same quantities evaluated using the simulation trajectory. Notably, key structural features postulated for the recent DSH model of nascent HDL are retained during the simulation, including (1) the superhelical conformation of the antiparallel apolipoprotein A1 (apoA1) chains, (2) the lipid micellar-pseudolamellar organization, and (3) the solvent-exposed Solar Flare loops, proposed sites of interaction with LCAT (lecithin cholesteryl acyltransferase). Analysis of salt bridge persistence during simulation provides insights into structural features of apoA1 that forms the backbone of the lipoprotein. The combination of molecular dynamics simulation and experimental data from a broad range of biophysical platforms serves as a powerful approach to studying large macromolecular assemblies such as lipoproteins. This application to nascent HDL validates the DSH model proposed earlier and suggests new structural details of nascent HDL.

Engineering extended membrane scaffold proteins for self-assembly of soluble nanoscale lipid bilayers

High-density lipoproteins (HDLs) play an important role in human health through the metabolism and trafficking of cholesterol as well as providing the feedstocks for steroid hormone biosynthesis. These particles contain proteins, primarily Apo-AI and phospholipid and progress through various structural forms including 'lipid-poor', 'discoidal' and 'spherical' entities as cholesterol esters and lipid are incorporated. The discoidal form of HDL is stabilized in solution by two encircling belts of Apo-AI. Previous protein engineering of the Apo-AI sequence has led to a series of amphipathic helical proteins, termed membrane scaffold proteins (MSPs), which have shown great value in assembling nanoscale soluble membrane bilayers, termed Nanodiscs, of homogeneous size and composition and in the assembly of numerous integral membrane proteins for biophysical and biochemical investigations. In this communication we document a protein engineering approach to generate and optimize an extended polypeptide MSP, which will self-assemble phospholipids into larger Nanodiscs with diameters of 16-17 nm. We extensively characterize these structures by size exclusion chromatography and solution X-ray scattering.

Conformation of dimeric apolipoprotein A-I milano on recombinant lipoprotein particles

Apolipoprotein A-I Milano (apoA-I(Milano)) is a naturally occurring human mutation of wild-type apolipoprotein A-I (apoA-I(WT)) having cystine substituted for arginine(173). Two molecules of apo-I(WT) form disks with phospholipid having a defined relationship between the apoA-I(WT) molecules. ApoA-I(Milano) forms cystine homodimers that would not allow the protein to adopt the conformation reported for apoA-I(WT). The conformational constraints for dimeric apoA-I(Milano) recombinant high-density lipoprotein (rHDL) disks made with phospholipid were deduced from a combination of chemical cross-linking and mass spectrometry. Lysine-selective homobifunctional cross-linkers were reacted with homogeneous rHDL having diameters of 78 and 125 A. After reduction, cross-linked apoA-I(Milano) was separated from monomeric apoprotein by gel electrophoresis and then subjected to in-gel trypsin digest. Cross-linked peptides were confirmed by MS/MS sequencing. The cross-links provided distance constraints that were used to refine models of lipid-bound dimeric apoA-I(Milano). These studies suggest that a single dimeric apoA-I(Milano) on 78 A diameter rHDL girdles the edge of a phospholipid disk assuming a "belt" conformation similar to the "belt" region of apoA-I(WT) on rHDL. However, the C-terminal end of dimeric apoA-I(Milano) wraps around the periphery of the particle to shield the fatty acid chains from water rather than folding back onto the "belt" as does apoA-I(WT). The two apoA-I(Milano) dimers on a 125 A diameter rHDL do not encircle the periphery of a phospholipid disk but appear to reside on the surface of a laminar micelle.

Coarse-grained ions without charges: reproducing the solvation structure of NaCl in water using short-ranged potentials

A coarse-grained model of NaCl in water is presented where the ions are modeled without charge to avoid computationally challenging electrostatics. A monatomic model of water [V. Molinero and E. B. Moore, J. Phys. Chem. B 113, 4008 (2009)] is used as the basis for this coarse-grain approach. The ability of Na(+) to disrupt the native tetrahedral arrangement of water molecules, and of Cl(-) to integrate within this organization, is preserved in this mW-ion model through parametrization focused on water's solvation of these ions. This model successfully reproduces the structural effect of ions on water, referenced to observations from experiments and atomistic molecular dynamics simulations, while using extremely short-ranged potentials. Without Coulomb interactions the model replicates details of the ion-water structure such as distinguishing contact and solvent-separated ion pairs and the free energy barriers between them. The approach of mimicking ionic effects with short-ranged interactions results in performance gains of two orders of magnitude compared to Ewald methods. Explored over a broad range of salt concentration, the model reproduces the solvation structure and trends of diffusion relative to atomistic simulations and experimental results. The functional form of the mW-ion model can be parametrized to represent other electrolytes. With increased computational efficiency and reliable structural fidelity, this model promises to be an asset for accessing significantly longer simulation time scales with an explicit solvent in a coarse-grained system involving, for example, polyelectrolytes such as proteins, nucleic acids, and fuel-cell membranes.

Membrane protein assembly into Nanodiscs

Nanodiscs are soluble nanoscale phospholipid bilayers which can self-assemble integral membrane proteins for biophysical, enzymatic or structural investigations. This means for rendering membrane proteins soluble at the single molecule level offers advantages over liposomes or detergent micelles in terms of size, stability, ability to add genetically modifiable features to the Nanodisc structure and ready access to both sides of the phospholipid bilayer domain. Thus the Nanodisc system provides a novel platform for understanding membrane protein function. We provide an overview of the Nanodisc approach and document through several examples many of the applications to the study of the structure and function of integral membrane proteins.

Double superhelix model of high density lipoprotein

High density lipoprotein (HDL), the carrier of so-called "good" cholesterol, serves as the major athero-protective lipoprotein and has emerged as a key therapeutic target for cardiovascular disease. We applied small angle neutron scattering (SANS) with contrast variation and selective isotopic deuteration to the study of nascent HDL to obtain the low resolution structure in solution of the overall time-averaged conformation of apolipoprotein AI (apoA-I) versus the lipid (acyl chain) core of the particle. Remarkably, apoA-I is observed to possess an open helical shape that wraps around a central ellipsoidal lipid phase. Using the low resolution SANS shapes of the protein and lipid core as scaffolding, an all-atom computational model for the protein and lipid components of nascent HDL was developed by integrating complementary structural data from hydrogen/deuterium exchange mass spectrometry and previously published constraints from multiple biophysical techniques. Both SANS data and the new computational model, the double superhelix model, suggest an unexpected structural arrangement of protein and lipids of nascent HDL, an anti-parallel double superhelix wrapped around an ellipsoidal lipid phase. The protein and lipid organization in nascent HDL envisages a potential generalized mechanism for lipoprotein biogenesis and remodeling, biological processes critical to sterol and lipid transport, organismal energy metabolism, and innate immunity.

Simulations of membrane tubulation by lattices of amphiphysin N-BAR domains

Membrane compartments of manifold shapes are found in cells, often sculpted by cellular proteins. In particular, proteins of the BAR domain superfamily participate in membrane-sculpting processes in vivo and reshape also in vitro low-curvature membrane liposomes into high-curvature tubes and vesicles. Here we show by means of computer simulations totaling over 1 millisecond, how lattices involving parallel rows of amphiphysin N-BAR domains sculpt flat membranes into tubes. A highly detailed, dynamic picture of the 100-microsecond formation of membrane tubes by lattices of N-BAR domains is obtained. Lattice types inducing a wide range of membrane curvatures, with radii approximately 15-100 nm, are explored. The results suggest that multiple lattice types are viable for efficient membrane bending.

Maturation of high-density lipoproteins

Human high-density lipoproteins (HDLs) are involved in the transport of cholesterol. The mechanism by which HDL assembles and functions is not well understood owing to a lack of structural information on circulating spherical HDL. Here, we report a series of molecular dynamics simulations that describe the maturation of discoidal HDL into spherical HDL upon incorporation of cholesterol ester as well as the resulting atomic level structure of a mature circulating spherical HDL particle. Sixty cholesterol ester molecules were added in a stepwise fashion to a discoidal HDL particle containing two apolipoproteins wrapped around a 160 dipalmitoylphosphatidylcholine lipid bilayer. The resulting matured particle, captured in a coarse-grained description, was then described in a consistent all-atom representation and analysed in chemical detail. The simulations show that maturation results from the formation of a highly dynamic hydrophobic core comprised of cholesterol ester surrounded by phospholipid and protein; the two apolipoprotein strands remain in a belt-like conformation as seen in the discoidal HDL particle, but with flexible N- and C-terminal helices and a central region stabilized by salt bridges. In the otherwise flexible lipoproteins, a less mobile central region provides an ideal location to bind lecithin cholesterol acyltransferase, the key enzyme that converts cholesterol to cholesterol ester during HDL maturation.

Transport-related structures and processes of the nuclear pore complex studied through molecular dynamics

Nuclear pore complexes (NPCs) are selectively gated pathways between nucleoplasm and cytoplasm. Whereas small molecules can diffuse freely through NPCs, large molecules (>40 kD) can pass only when bound to transport receptors. The NPC central channel is filled with disordered proteins, rich in phenylalanine-glycine (FG) repeats, referred to as FG-nups. Our simulations, carried out at coarse-grained and all-atom levels, show that arrays of FG-nups tethered to a planar surface, at an FG-repeat density found in the NPC, form dynamic brush-like structures of multiprotein bundles, whereas individual FG-nups form dynamic globular structures. More than half of the FG-repeats are found on the surface of the bundles, offering a favorable environment for transport receptors. Binding to FG-repeats and a sliding motion of NTF2 induced by binding and unbinding to phenylalanines were observed when adding this transport receptor into one of the brush-like structures.

Four-scale description of membrane sculpting by BAR domains

BAR domains are proteins that sense and sculpt curved membranes in cells, furnishing a relatively well-studied example of mechanisms employed in cellular morphogenesis. We report a computational study of membrane bending by BAR domains at four levels of resolution, described by 1), all-atom molecular dynamics; 2), residue-based coarse-graining (resolving single amino acids and lipid molecules); 3), shape-based coarse-graining (resolving overall protein and membrane shapes); and 4), a continuum elastic membrane model. Membrane sculpting performed by BAR domains collectively is observed in agreement with experiments. Different arrangements of BAR domains on the membrane surface are found to lead to distinct membrane curvatures and bending dynamics.

Three-dimensional models of HDL apoA-I: implications for its assembly and function

The purpose of this review is to highlight recent advances toward the refinement of a three-dimensional structure for lipid-bound apolipoprotein A-I (apoA-I) on recombinant HDL. Recently, X-ray crystallography has yielded a new structure for full-length, lipid-free apoA-I. Although this approach has not yet been successful in solving the three-dimensional structure of lipid-bound apoA-I, analysis of the X-ray structures has been of immense help in the interpretation of structural data obtained from other methods that yield structural information. Recent studies emphasize the use of mass spectrometry to unambiguously identify cross-linked peptides or to quantify solvent accessibility using hydrogen-deuterium exchange. The combination of mass spectrometry, molecular modeling, molecular dynamic analysis, and small-angle X-ray diffraction has provided additional structural information on apoA-I folding that complements previous approaches.

Quantifying size distributions of nanolipoprotein particles with single-particle analysis and molecular dynamic simulations

Self-assembly of purified apolipoproteins and phospholipids results in the formation of nanometer-sized lipoprotein complexes, referred to as nanolipoprotein particles (NLPs). These bilayer constructs are fully soluble in aqueous environments and hold great promise as a model system to aid in solubilizing membrane proteins. Size variability in the self-assembly process has been recognized for some time, yet limited studies have been conducted to examine this phenomenon. Understanding the source of this heterogeneity may lead to methods to mitigate heterogeneity or to control NLP size, which may be important for tailoring NLPs for specific membrane proteins. Here, we have used atomic force microscopy, ion mobility spectrometry, and transmission electron microscopy to quantify NLP size distributions on the single-particle scale, specifically focusing on assemblies with 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and a recombinant apolipoprotein E variant containing the N-terminal 22 kDa fragment (E422k). Four discrete sizes of E422k/DMPC NLPs were identified by all three techniques, with diameters centered at approximately 14.5, 19, 23.5, and 28 nm. Computer simulations suggest that these sizes are related to the structure and number of E422k lipoproteins surrounding the NLPs and particles with an odd number of lipoproteins are consistent with the double-belt model, in which at least one lipoprotein adopts a hairpin structure.

Molecular models need to be tested: the case of a solar flares discoidal HDL model

In the absence of atomic structures of high-density lipoproteins in their lipid-bound states, many molecular models have been produced based on experimental data. Using molecular dynamics, we show that a recently proposed "solar-flares" model of discoidal high-density lipoprotein is implausible. Our simulations show a collapse of the protruding solar-flare loops and a notable protein rearrangement due to an energetically unfavorable orientation of the hydrophobic protein surface toward the aqueous solvent.