Thermal stability of apolipoprotein A-I in high-density lipoproteins by molecular dynamics

Deciphering structural aspects of reverse cholesterol transport: mapping the knowns and unknowns

Atherosclerosis is a major contributor to the onset and progression of cardiovascular disease (CVD). Cholesterol-loaded foam cells play a pivotal role in forming atherosclerotic plaques. Induction of cholesterol efflux from these cells may be a promising approach in treating CVD. The reverse cholesterol transport (RCT) pathway delivers cholesteryl ester (CE) packaged in high-density lipoproteins (HDL) from non-hepatic cells to the liver, thereby minimising cholesterol load of peripheral cells. RCT takes place via a well-organised interplay amongst apolipoprotein A1 (ApoA1), lecithin cholesterol acyltransferase (LCAT), ATP binding cassette transporter A1 (ABCA1), scavenger receptor-B1 (SR-B1), and the amount of free cholesterol. Unfortunately, modulation of RCT for treating atherosclerosis has failed in clinical trials owing to our lack of understanding of the relationship between HDL function and RCT. The fate of non-hepatic CEs in HDL is dependent on their access to proteins involved in remodelling and can be regulated at the structural level. An inadequate understanding of this inhibits the design of rational strategies for therapeutic interventions. Herein we extensively review the structure-function relationships that are essential for RCT. We also focus on genetic mutations that disturb the structural stability of proteins involved in RCT, rendering them partially or completely non-functional. Further studies are necessary for understanding the structural aspects of RCT pathway completely, and this review highlights alternative theories and unanswered questions.

Molecular dynamics simulations of lipid nanodiscs

A lipid nanodisc is a discoidal lipid bilayer stabilized by proteins, peptides, or polymers on its edge. Nanodiscs have two important connections to structural biology. The first is associated with high-density lipoprotein (HDL), a particle with a variety of functionalities including lipid transport. Nascent HDL (nHDL) is a nanodisc stabilized by Apolipoprotein A-I (APOA1). Determining the structure of APOA1 and its mimetic peptides in nanodiscs is crucial to understanding pathologies related to HDL maturation and designing effective therapies. Secondly, nanodiscs offer non-detergent membrane-mimicking environments and greatly facilitate structural studies of membrane proteins. Although seemingly similar, natural and synthetic nanodiscs are different in that nHDL is heterogeneous in size, due to APOA1 elasticity, and gradually matures to become spherical. Synthetic nanodiscs, in contrast, should be homogenous, stable, and size-tunable. This report reviews previous molecular dynamics (MD) simulation studies of nanodiscs and illustrates convergence and accuracy issues using results from new multi-microsecond atomistic MD simulations. These new simulations reveal that APOA1 helices take 10-20 μs to rearrange on the nanodisc, while peptides take 2 μs to migrate from the disc surfaces to the edge. These systems can also become kinetically trapped depending on the initial conditions. For example, APOA1 was trapped in a biologically irrelevant conformation for the duration of a 10 μs trajectory; the peptides were similarly trapped for 5 μs. It therefore remains essential to validate MD simulations of these systems with experiments due to convergence and accuracy issues. This article is part of a Special Issue entitled: Emergence of Complex Behavior in Biomembranes edited by Marjorie Longo.

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.

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.

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.

A robust all-atom model for LCAT generated by homology modeling

LCAT is activated by apoA-I to form cholesteryl ester. We combined two structures, phospholipase A2 (PLA2) that hydrolyzes the ester bond at the sn-2 position of oxidized (short) acyl chains of phospholipid, and bacteriophage tubulin PhuZ, as C- and N-terminal templates, respectively, to create a novel homology model for human LCAT. The juxtaposition of multiple structural motifs matching experimental data is compelling evidence for the general correctness of many features of the model: i) The N-terminal 10 residues of the model, required for LCAT activity, extend the hydrophobic binding trough for the sn-2 chain 15-20 Å relative to PLA2. ii) The topography of the trough places the ester bond of the sn-2 chain less than 5 Å from the hydroxyl of the catalytic nucleophile, S181. iii) A β-hairpin resembling a lipase lid separates S181 from solvent. iv) S181 interacts with three functionally critical residues: E149, that regulates sn-2 chain specificity, and K128 and R147, whose mutations cause LCAT deficiency. Because the model provides a novel explanation for the complicated thermodynamic problem of the transfer of hydrophobic substrates from HDL to the catalytic triad of LCAT, it is an important step toward understanding the antiatherogenic role of HDL in reverse cholesterol transport.

An experimentally robust model of monomeric apolipoprotein A-I created from a chimera of two X-ray structures and molecular dynamics simulations

High-density lipoprotein (HDL) retards atherosclerosis by accepting cholesterol from the artery wall. However, the structure of the proposed acceptor, monomeric apolipoprotein A-I (apoA-I), the major protein of HDL, is poorly understood. Two published models for monomeric apoA-I used cross-linking distance constraints to derive best fit conformations. This approach has limitations. (i) Cross-linked peptides provide no information about secondary structure. (ii) A protein chain can be folded in multiple ways to create a best fit. (iii) Ad hoc folding of a secondary structure is unlikely to produce a stable orientation of hydrophobic and hydrophilic residues. To address these limitations, we used a different approach. We first noted that the dimeric apoA-I crystal structure, (Δ185-243)apoA-I, is topologically identical to a monomer in which helix 5 forms a helical hairpin, a monomer with a hydrophobic cleft running the length of the molecule. We then realized that a second crystal structure, (Δ1-43)apoA-I, contains a C-terminal structure that fits snuggly via aromatic and hydrophobic interactions into the hydrophobic cleft. Consequently, we combined these crystal structures into an initial model that was subjected to molecular dynamics simulations. We tested the initial and simulated models and the two previously published models in three ways: against two published data sets (domains predicted to be helical by H/D exchange and six spin-coupled residues) and against our own experimentally determined cross-linking distance constraints. We note that the best fit simulation model, superior by all tests to previously published models, has dynamic features of a molten globule with interesting implications for the functions of apoA-I.

The structure of human apolipoprotein A-IV as revealed by stable isotope-assisted cross-linking, molecular dynamics, and small angle x-ray scattering

Apolipoprotein (apo)A-IV plays important roles in dietary lipid and glucose metabolism, and knowledge of its structure is required to fully understand the molecular basis of these functions. However, typical of the entire class of exchangeable apolipoproteins, its dynamic nature and affinity for lipid has posed challenges to traditional high resolution structural approaches. We previously reported an x-ray crystal structure of a dimeric truncation mutant of apoA-IV, which showed a unique helix-swapping molecular interface. Unfortunately, the structures of the N and C termini that are important for lipid binding were not visualized. To build a more complete model, we used chemical cross-linking to derive distance constraints across the full-length protein. The approach was enhanced with stable isotope labeling to overcome ambiguities in determining molecular span of the cross-links given the remarkable similarities in the monomeric and dimeric apoA-IV structures. Using 51 distance constraints, we created a starting model for full-length monomeric apoA-IV and then subjected it to two modeling approaches: (i) molecular dynamics simulations and (ii) fitting to small angle x-ray scattering data. This resulted in the most detailed models yet for lipid-free monomeric or dimeric apoA-IV. Importantly, these models were of sufficient detail to direct the experimental identification of new functional residues that participate in a "clasp" mechanism to modulate apoA-IV lipid affinity. The isotope-assisted cross-linking approach should prove useful for further study of this family of apolipoproteins in both the lipid-free and -bound states.

Mechanistic picture for conformational transition of a membrane transporter at atomic resolution

During their transport cycle, ATP-binding cassette (ABC) transporters undergo large-scale conformational changes between inward- and outward-facing states. Using an approach based on designing system-specific reaction coordinates and using nonequilibrium work relations, we have performed extensive all-atom molecular dynamics simulations in the presence of explicit membrane/solvent to sample a large number of mechanistically distinct pathways for the conformational transition of MsbA, a bacterial ABC exporter whose structure has been solved in multiple functional states. The computational approach developed here is based on (i) extensive exploration of system-specific biasing protocols (e.g., using collective variables designed based on available low-resolution crystal structures) and (ii) using nonequilibrium work relations for comparing the relevance of the transition pathways. The most relevant transition pathway identified using this approach involves several distinct stages reflecting the complex nature of the structural changes associated with the function of the protein. The opening of the cytoplasmic gate during the outward- to inward-facing transition of apo MsbA is found to be disfavored when the periplasmic gate is open and facilitated by a twisting motion of the nucleotide-binding domains that involves a dramatic change in their relative orientation. These results highlight the cooperativity between the transmembrane and the nucleotide-binding domains in the conformational transition of ABC exporters. The approach introduced here provides a framework to study large-scale conformational changes of other membrane transporters whose computational investigation at an atomic resolution may not be currently feasible using conventional methods.

Conservation of apolipoprotein A-I's central domain structural elements upon lipid association on different high-density lipoprotein subclasses

The antiatherogenic properties of apolipoprotein A-I (apoA-I) are derived, in part, from lipidation-state-dependent structural elements that manifest at different stages of apoA-I's progression from lipid-free protein to spherical high-density lipoprotein (HDL). Previously, we reported the structure of apoA-I's N-terminus on reconstituted HDLs (rHDLs) of different sizes. We have now investigated at the single-residue level the conformational adaptations of three regions in the central domain of apoA-I (residues 119-124, 139-144, and 164-170) upon apoA-I lipid binding and HDL formation. An important function associated with these residues of apoA-I is the activation of lecithin:cholesterol acyltransferase (LCAT), the enzyme responsible for catalyzing HDL maturation. Structural examination was performed by site-directed tryptophan fluorescence and spin-label electron paramagnetic resonance spectroscopies for both the lipid-free protein and rHDL particles 7.8, 8.4, and 9.6 nm in diameter. The two methods provide complementary information about residue side chain mobility and molecular accessibility, as well as the polarity of the local environment at the targeted positions. The modulation of these biophysical parameters yielded new insight into the importance of structural elements in the central domain of apoA-I. In particular, we determined that the loosely lipid-associated structure of residues 134-145 is conserved in all rHDL particles. Truncation of this region completely abolished LCAT activation but did not significantly affect rHDL size, reaffirming the important role of this structural element in HDL function.

Structural basis for distinct functions of the naturally occurring Cys mutants of human apolipoprotein A-I

HDL removes cell cholesterol and protects against atherosclerosis. ApoA-I provides a flexible structural scaffold and an important functional ligand on the HDL surface. We propose structural models for apoA-I(Milano) (R173C) and apoA-I(Paris) (R151C) mutants that show high cardioprotection despite low HDL levels. Previous studies established that two apoA-I molecules encircle HDL in an antiparallel, helical double-belt conformation. Recently, we solved the atomic structure of lipid-free Δ(185-243)apoA-I and proposed a conformational ensemble for apoA-I(WT) on HDL. Here we modify this ensemble to understand how intermolecular disulfides involving C173 or C151 influence protein conformation. The double-belt conformations are modified by belt rotation, main-chain unhinging around Gly, and Pro-induced helical bending, and they are verified by comparison with previous experimental studies and by molecular dynamics simulations of apoA-I(Milano) homodimer. In our models, the molecular termini repack on various-sized HDL, while packing around helix-5 in apoA-I(WT), helix-6 in apoA-I(Paris), or helix-7 in apoA-I(Milano) homodimer is largely conserved. We propose how the disulfide-induced constraints alter the protein conformation and facilitate dissociation of the C-terminal segment from HDL to recruit additional lipid. Our models unify previous studies of apoA-I(Milano) and demonstrate how the mutational effects propagate to the molecular termini, altering their conformations, dynamics, and function.

MD simulations suggest important surface differences between reconstituted and circulating spherical HDL

Since spheroidal HDL particles (sHDL) are highly dynamic, molecular dynamics (MD) simulations are useful for obtaining structural models. Here we use MD to simulate sHDL with stoichiometries of reconstituted and circulating particles. The hydrophobic effect during simulations rapidly remodels discoidal HDL containing mixed lipids to sHDL containing a cholesteryl ester/triglyceride (CE/TG) core. We compare the results of simulations of previously characterized reconstituted sHDL particles containing two or three apoA-I created in the absence of phospholipid transfer protein (PLTP) with simulations of circulating human HDL containing two or three apoA-I without apoA-II. We find that circulating sHDL compared with reconstituted sHDL with the same number of apoA-I per particle contain approximately equal volumes of core lipid but significantly less surface lipid monolayers. We conclude that in vitro reconstituted sHDL particles contain kinetically trapped excess phospholipid and are less than ideal models for circulating sHDL particles. In the circulation, phospholipid transfer via PLTP decreases the ratio of phospholipid to apolipoprotein for all sHDL particles. Further, sHDL containing two or three apoA-I adapt to changes in surface area by condensation of common conformational motifs. These results represent an important step toward resolving the complicated issue of the protein and lipid stoichiometry of circulating HDL.

Validation of previous computer models and MD simulations of discoidal HDL by a recent crystal structure of apoA-I

HDL is a population of apoA-I-containing particles inversely correlated with heart disease. Because HDL is a soft form of matter deformable by thermal fluctuations, structure determination has been difficult. Here, we compare the recently published crystal structure of lipid-free (Δ185-243)apoA-I with apoA-I structure from models and molecular dynamics (MD) simulations of discoidal HDL. These analyses validate four of our previous structural findings for apoA-I: i) a baseline double belt diameter of 105 Å ii) central α helixes with an 11/3 pitch; iii) a "presentation tunnel" gap between pairwise helix 5 repeats hypothesized to move acyl chains and unesterified cholesterol from the lipid bilayer to the active sites of LCAT; and iv) interchain salt bridges hypothesized to stabilize the LL5/5 chain registry. These analyses are also consistent with our finding that multiple salt bridge-forming residues in the N-terminus of apoA-I render that conserved domain "sticky." Additionally, our crystal MD comparisons led to two new hypotheses: i) the interchain leucine-zippers previously reported between the pair-wise helix 5 repeats drive lipid-free apoA-I registration; ii) lipidation induces rotations of helix 5 to allow formation of interchain salt bridges, creating the LCAT presentation tunnel and "zip-locking" apoA-I into its full LL5/5 registration.

Conformational coupling of the nucleotide-binding and the transmembrane domains in ABC transporters

Basic architecture of ABC transporters includes two transmembrane domains (TMDs) and two nucleotide-binding domains (NBDs). Although the transport process takes place in the TMDs, which provide the substrate translocation pathway across the cell membrane and control its accessibility between the two sides of the membrane, the energy required for the process is provided by conformational changes induced in the NBDs by binding and hydrolysis of ATP. Nucleotide-dependent conformational changes in the NBDs, therefore, need to be coupled to structural changes in the TMDs. Using molecular dynamics simulations, we have investigated the structural elements involved in the conformational coupling between the NBDs and the TMDs in the Escherichia coli maltose transporter, an ABC importer for which an intact structure is available both in inward-facing and outward-facing conformations. The prevailing model of coupling is primarily based on a single structural motif, known as the coupling helices, as the main structural element for the NBD-TMD coupling. Surprisingly, we find that in the absence of the NBDs the coupling helices can be conformationally decoupled from the rest of the TMDs, despite their covalent connection. That is, the structural integrity of the coupling helices and their tight coupling to the core of the TMDs rely on the contacts provided by the NBDs. Based on the conformational and dynamical analysis of the simulation trajectories, we propose that the core coupling elements in the maltose transporter involve contributions from several structural motifs located at the NBD-TMD interface, namely, the EAA loops from the TMDs, and the Q-loop and the ENI motifs from the NBDs. These three structural motifs in small ABC importers show a high degree of correlation in motion and mediate the necessary conformational coupling between the core of TMDs and the helical subdomains of NBDs. A comprehensive analysis of the structurally known ABC transporters shows a high degree of conservation of the identified 3-motif coupling elements only in the subfamily of small ABC importers, suggesting a distinct mode of NBD-TMD coupling from the other two major ABC transporter folds, namely large ABC importers and ABC exporters.

Crystal structure of C-terminal truncated apolipoprotein A-I reveals the assembly of high density lipoprotein (HDL) by dimerization

Apolipoprotein A-I (apoA-I) plays important structural and functional roles in plasma high density lipoprotein (HDL) that is responsible for reverse cholesterol transport. However, a molecular understanding of HDL assembly and function remains enigmatic. The 2.2-Å crystal structure of Δ(185-243)apoA-I reported here shows that it forms a half-circle dimer. The backbone of the dimer consists of two elongated antiparallel proline-kinked helices (five AB tandem repeats). The N-terminal domain of each molecule forms a four-helix bundle with the helical C-terminal region of the symmetry-related partner. The central region forms a flexible domain with two antiparallel helices connecting the bundles at each end. The two-domain dimer structure based on helical repeats suggests the role of apoA-I in the formation of discoidal HDL particles. Furthermore, the structure suggests the possible interaction with lecithin-cholesterol acyltransferase and may shed light on the molecular details of the effect of the Milano, Paris, and Fin mutations.

Local/bulk determinants of conformational stability of exchangeable apolipoproteins

GuHCl-induced denaturation of human plasma apoA-I, apoA-II, apoA-IV, apoE3 and three recombinant apoE isoforms in solution and discoidal complexes with phosphatidylcholine (only plasma proteins) was studied. The protein conformational stability (ΔG(H(2)O)) and a slope of linear dependence of free energy of unfolding on GuHCl concentration (m-value) were estimated with the three equilibrium schemes. The data for all proteins, except apoA-II, fit with the three-state model, thus evidencing two-domain structure. The predicted folding rate of the four apoE in solution correlated with conformational stability. The dependence disappeared at the inclusion of apoA-I and apoA-IV into analysis and the m-values, adjusted for residue number in helices (m(rh)), differed between those for apoE and apoA-I/apoA-IV. However, the m(rh)-values for six proteins correlated positively with the fractional change in accessible surface area at unfolding for Phe, Lys and Asn, while negatively for Arg, Ala and Gly residues. The difference between the adjusted ΔG(rh)(H(2)O) values for apolipoproteins in complexes and in solution decreased at the increase of reduced temperature (T(obs)-T(t))/T(t). The induction of intrinsic disorder by arginine residues may be of primary importance in metabolism and function of exchangeable apolipoproteins, while their stability in nascent discoidal HDL is controlled by the physical state of phosphatidylcholine.

Alteration of negatively charged residues in the 89 to 99 domain of apoA-I affects lipid homeostasis and maturation of HDL

In this study, we investigated the role of positively and negatively charged amino acids within the 89-99 region of apolipoprotein A-I (apoA-I), which are highly conserved in mammals, on plasma lipid homeostasis and the biogenesis of HDL. We previously showed that deletion of the 89-99 region of apoA-I increased plasma cholesterol and phospholipids, but it did not affect plasma triglycerides. Functional studies using adenovirus-mediated gene transfer of two apoA-I mutants in apoA-I-deficient mice showed that apoA-I[D89A/E91A/E92A] increased plasma cholesterol and caused severe hypertriglyceridemia. HDL levels were reduced, and approximately 40% of the apoA-I was distributed in VLDL/IDL. The HDL consisted of mostly spherical and a few discoidal particles and contained preβ1 and α4-HDL subpopulations. The lipid, lipoprotein, and HDL profiles generated by the apoA-I[K94A/K96A] mutant were similar to those of wild-type (WT) apoA-I. Coexpression of apoA-I[D89A/E91A/E92A] and human lipoprotein lipase abolished hypertriglyceridemia, restored in part the α1,2,3,4 HDL subpopulations, and redistributed apoA-I in the HDL2/HDL3 regions, but it did not prevent the formation of discoidal HDL particles. Physicochemical studies showed that the apoA-I[D89A/E91A/E92A] mutant had reduced α-helical content and effective enthalpy of thermal denaturation, increased exposure of hydrophobic surfaces, and increased affinity for triglyceride-rich emulsions. We conclude that residues D89, E91, and E92 of apoA-I are important for plasma cholesterol and triglyceride homeostasis as well as for the maturation of HDL.

"Sticky" and "promiscuous", the yin and yang of apolipoprotein A-I termini in discoidal high-density lipoproteins: a combined computational-experimental approach

Apolipoprotein (apo) A-I-containing lipoproteins in the form of high-density lipoproteins (HDL) are inversely correlated with atherosclerosis. Because HDL is a soft form of condensed matter easily deformable by thermal fluctuations, the molecular mechanisms for HDL remodeling are not well understood. A promising approach to understanding HDL structure and dynamics is molecular dynamics (MD). In the present study, two computational strategies, MD simulated annealing (MDSA) and MD temperature jump, were combined with experimental particle reconstitution to explore molecular mechanisms for phospholipid- (PL-) rich HDL particle remodeling. The N-terminal domains of full-length apoA-I were shown to be "sticky", acting as a molecular latch largely driven by salt bridges, until, at a critical threshold of particle size, the associated domains released to expose extensive hydrocarbon regions of the PL to solvent. The "sticky" N-termini also associate with other apoA-I domains, perhaps being involved in N-terminal loops suggested by other laboratories. Alternatively, the overlapping helix 10 C-terminal domains of apoA-I were observed to be extremely mobile or "promiscuous", transiently exposing limited hydrocarbon regions of PL. Based upon these models and reconstitution studies, we propose that separation of the N-terminal domains, as particles exceed a critical size, triggers fusion between particles or between particles and membranes, while the C-terminal domains of apoA-I drive the exchange of polar lipids down concentration gradients between particles. This hypothesis has significant biological relevance since lipid exchange and particle remodeling are critically important processes during metabolism of HDL particles at every step in the antiatherogenic process of reverse cholesterol transport.

Sequence conservation of apolipoprotein A-I affords novel insights into HDL structure-function

We performed alignment of apolipoprotein A-I (apoA-I) sequences from 31 species of animals. We found there is specific conservation of salt bridge-forming residues in the first 30 residues of apoA-I and general conservation of a variety of residue types in the central domain, helix 2/3 to helix 7/8. In the lipid-associating domain, helix 7 and helix 10 are the most and least conserved helixes, respectively. Furthermore, eight residues are completely conserved: P66, R83, P121, E191, and P220, and three of seven Tyr residues in human apoA-I, Y18, Y115, and Y192, are conserved. Residue Y18 appears to be important for assembly of HDL. E191-Y192 represents the only completely conserved pair of adjacent residues in apoA-I; Y192 is a preferred target for site-specific oxidative modification within atheroma, and molecular dynamic simulations suggest that the conserved pair E191-Y192 is in a solvent-exposed loop-helix-loop. Molecular dynamics testing of human apoA-I showed that M112 and M148 interact with Y115, a microenvironment unique to human apoA-I. Finally, conservation of Arg residues in the α11/3 helical wheel position 7 supports several possibilities: interactions with adjacent phospholipid molecules and/or oxidized lipids and/or binding of antioxidant enzymes through cation-π orbital interactions. We conclude that sequence alignment of apoA-I provides unique insights into apoA-I structure-function relationship.

Capturing Functional Motions of Membrane Channels and Transporters with Molecular Dynamics Simulation

Conformational changes of proteins are involved in all aspects of protein function in biology. Almost all classes of proteins respond to changes in their environment, ligand binding, and interaction with other proteins and regulatory agents through undergoing conformational changes of various degrees and magnitudes. Membrane channels and transporters are the major classes of proteins that are responsible for mediating efficient and selective transport of materials across the cellular membrane. Similar to other proteins, they take advantage of conformational changes to make transitions between various functional states. In channels, large-scale conformational changes are mostly involved in the process of "gating", i.e., opening and closing of the pore of the channel protein in response to various signals. In transporters, conformational changes constitute various steps of the conduction process, and, thus, are more closely integrated in the transport process. Owing to significant progress in developing highly efficient parallel algorithms in molecular dynamics simulations and increased computational resources, and combined with the availability of high-resolution, atomic structures of membrane proteins, we are in an unprecedented position to use computer simulation and modeling methodologies to investigate the mechanism of function of membrane channels and transporters. While the entire transport cycle is still out of reach of current methodologies, many steps involved in the function of transport proteins have been characterized with molecular dynamics simulations. Here, we present several examples of such studies from our laboratory, in which functionally relevant conformational changes of membrane channels and transporters have been characterized using extended simulations.

Assessment of the validity of the double superhelix model for reconstituted high density lipoproteins: a combined computational-experimental approach

For several decades, the standard model for high density lipoprotein (HDL) particles reconstituted from apolipoprotein A-I (apoA-I) and phospholipid (apoA-I/HDL) has been a discoidal particle ∼100 Å in diameter and the thickness of a phospholipid bilayer. Recently, Wu et al. (Wu, Z., Gogonea, V., Lee, X., Wagner, M. A., Li, X. M., Huang, Y., Undurti, A., May, R. P., Haertlein, M., Moulin, M., Gutsche, I., Zaccai, G., Didonato, J. A., and Hazen, S. L. (2009) J. Biol. Chem. 284, 36605-36619) used small angle neutron scattering to develop a new model they termed double superhelix (DSH) apoA-I that is dramatically different from the standard model. Their model possesses an open helical shape that wraps around a prolate ellipsoidal type I hexagonal lyotropic liquid crystalline phase. Here, we used three independent approaches, molecular dynamics, EM tomography, and fluorescence resonance energy transfer spectroscopy (FRET) to assess the validity of the DSH model. (i) By using molecular dynamics, two different approaches, all-atom simulated annealing and coarse-grained simulation, show that initial ellipsoidal DSH particles rapidly collapse to discoidal bilayer structures. These results suggest that, compatible with current knowledge of lipid phase diagrams, apoA-I cannot stabilize hexagonal I phase particles of phospholipid. (ii) By using EM, two different approaches, negative stain and cryo-EM tomography, show that reconstituted apoA-I/HDL particles are discoidal in shape. (iii) By using FRET, reconstituted apoA-I/HDL particles show a 28-34-Å intermolecular separation between terminal domain residues 40 and 240, a distance that is incompatible with the dimensions of the DSH model. Therefore, we suggest that, although novel, the DSH model is energetically unfavorable and not likely to be correct. Rather, we conclude that all evidence supports the likelihood that reconstituted apoA-I/HDL particles, in general, are discoidal in shape.

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.

Isolation, characterization, and stability of discretely-sized nanolipoprotein particles assembled with apolipophorin-III

BACKGROUND

Nanolipoprotein particles (NLPs) are discoidal, nanometer-sized particles comprised of self-assembled phospholipid membranes and apolipoproteins. NLPs assembled with human apolipoproteins have been used for myriad biotechnology applications, including membrane protein solubilization, drug delivery, and diagnostic imaging. To expand the repertoire of lipoproteins for these applications, insect apolipophorin-III (apoLp-III) was evaluated for the ability to form discretely-sized, homogeneous, and stable NLPs.

METHODOLOGY

Four NLP populations distinct with regards to particle diameters (ranging in size from 10 nm to >25 nm) and lipid-to-apoLp-III ratios were readily isolated to high purity by size exclusion chromatography. Remodeling of the purified NLP species over time at 4 degrees C was monitored by native gel electrophoresis, size exclusion chromatography, and atomic force microscopy. Purified 20 nm NLPs displayed no remodeling and remained stable for over 1 year. Purified NLPs with 10 nm and 15 nm diameters ultimately remodeled into 20 nm NLPs over a period of months. Intra-particle chemical cross-linking of apoLp-III stabilized NLPs of all sizes.

CONCLUSIONS

ApoLp-III-based NLPs can be readily prepared, purified, characterized, and stabilized, suggesting their utility for biotechnological applications.

Myeloperoxidase: an oxidative pathway for generating dysfunctional high-density lipoprotein

Accumulation of low-density lipoprotein (LDL)-derived cholesterol by artery wall macrophages triggers atherosclerosis, the leading cause of cardiovascular disease. Conversely, high-density lipoprotein (HDL) retards atherosclerosis by promoting cholesterol efflux from macrophages by the membrane-associated ATP-binding cassette transporter A1 (ABCA1) pathway. HDL has been proposed to lose its cardioprotective effects in subjects with atherosclerosis, but the underlying mechanisms are poorly understood. One potential pathway involves oxidative damage by myeloperoxidase (MPO), a heme enzyme secreted by human artery wall macrophages. We used mass spectrometry to demonstrate that HDL isolated from patients with established cardiovascular disease contains elevated levels of 3-chlorotyrosine and 3-nitrotyrosine, two characteristic products of MPO. When apolipoprotein A-I (apoA-I), the major HDL protein, was oxidized by MPO, its ability to promote cellular cholesterol efflux by ABCA1 was impaired. Moreover, oxidized apoA-I was unable to activate lecithin:cholesterol acyltransferase (LCAT), which rapidly converts free cholesterol to cholesteryl ester, a critical step in HDL maturation. Biochemical studies implicated tyrosine chlorination and methionine oxygenation in the loss of ABCA1 and LCAT activity by oxidized apoA-I. Oxidation of specific residues in apoA-I inhibited two key steps in cholesterol efflux from macrophages, raising the possibility that MPO initiates a pathway for generating dysfunctional HDL in humans.

Pressure perturbation calorimetry of apolipoproteins in solution and in model lipoproteins

High-density lipoproteins (HDLs) are complexes of lipids and proteins (termed apolipoproteins) that remove cell cholesterol and protect from atherosclerosis. Apolipoproteins contain amphipathic alpha-helices that have high content (> or = 1/3) and distinct distribution of charged and apolar residues, adopt molten globule-like conformations in solution, and bind to lipid surfaces. We report the first pressure perturbation calorimetry (PPC) study of apolipoproteins. In solution, the main HDL protein, apoA-I, shows relatively large volume contraction, DeltaV(unf) = -0.33%, and an apparent reduction in thermal expansivity upon unfolding, Deltaalpha(unf) < or = 0, which has not been observed in other proteins. We propose that these values are dominated by increased charged residue hydration upon alpha-helical unfolding, which may result from disruption of multiple salt bridges. At 5 degrees C, apoA-I shows large thermal expansion coefficient, alpha(5 degrees) = 15.10(-4) K(-1), that rapidly declines upon heating from 5 to 40 degrees C, alpha(40 degrees) - alpha(5 degrees) = -4.10(-4) K(-1); apolipoprotein C-I shows similar values of alpha(5 degrees) and alpha(40 degrees). These values are larger than in globular proteins. They indicate dominant effect of charged residue hydration, which may modulate functional apolipoprotein interactions with a broad range of their protein and lipid ligands. The first PPC analysis of a protein-lipid complex is reported, which focuses on the chain melting transition in model HDL containing apoA-I or apoC-I, dimyristoyl phosphatidylcholine, and 0-20% cholesterol. The results may provide new insights into volumetric properties of HDL that modulate metabolic lipoprotein remodeling during cholesterol transport.

Dynamics of activation of lecithin:cholesterol acyltransferase by apolipoprotein A-I

The product of transesterification of phospholipid acyl chains and unesterified cholesterol (UC) by the enzyme lecithin:cholesterol acyltransferase (LCAT) is cholesteryl ester (CE). Activation of LCAT by apolipoprotein (apo) A-I on nascent (discoidal) high-density lipoproteins (HDL) is essential for formation of mature (spheroidal) HDL during the antiatherogenic process of reverse cholesterol transport. Here we report all-atom and coarse-grained (CG) molecular dynamics (MD) simulations of HDL particles that have major implications for mechanisms of LCAT activation. Both the all-atom and CG simulations provide support for a model in which the helix 5/5 domains of apoA-I create an amphipathic "presentation tunnel" that exposes methyl ends of acyl chains at the bilayer center to solvent. Further, CG simulations show that UC also becomes inserted with high efficiency into the amphipathic presentation tunnel with its hydroxyl moiety (UC-OH) exposed to solvent; these results are consistent with trajectory analyses of the all-atom simulations showing that UC is being concentrated in the vicinity of the presentation tunnel. Finally, consistent with known product inhibition of CE-rich HDL by CE, CG simulations of CE-rich spheroidal HDL indicate partial blockage of the amphipathic presentation tunnel by CE. These results lead us to propose the following working hypothesis. After attachment of LCAT to discoidal HDL, the helix 5/5 domains in apoA-I form amphipathic presentation tunnels for migration of hydrophobic acyl chains and amphipathic UC from the bilayer to the phospholipase A2-like and esterification active sites of LCAT, respectively. This hypothesis is currently being tested by site-directed mutagenesis.

Structures of discoidal high density lipoproteins: a combined computational-experimental approach

Conversion of discoidal phospholipid (PL)-rich high density lipoprotein (HDL) to spheroidal cholesteryl ester-rich HDL is a central step in reverse cholesterol transport. A detailed understanding of this process and the atheroprotective role of apolipoprotein A-I (apoA-I) requires knowledge of the structure and dynamics of these various particles. This study, combining computation with experimentation, illuminates structural features of apoA-I allowing it to incorporate varying amounts of PL. Molecular dynamics simulated annealing of PL-rich HDL models containing unesterified cholesterol results in double belt structures with the same general saddle-shaped conformation of both our previous molecular dynamics simulations at 310 K and the x-ray structure of lipid-free apoA-I. Conversion from a discoidal to a saddle-shaped particle involves loss of helicity and formation of loops in opposing antiparallel parts of the double belt. During surface expansion caused by the temperature-jump step, the curved palmitoyloleoylphosphatidylcholine bilayer surfaces approach planarity. Relaxation back into saddle-shaped structures after cool down and equilibration further supports the saddle-shaped particle model. Our kinetic analyses of reconstituted particles demonstrate that PL-rich particles exist in discrete sizes corresponding to local energetic minima. Agreement of experimental and computational determinations of particle size/shape and apoA-I helicity provide additional support for the saddle-shaped particle model. Truncation experiments combined with simulations suggest that the N-terminal proline-rich domain of apoA-I influences the stability of PL-rich HDL particles. We propose that apoA-I incorporates increasing PL in the form of minimal surface bilayers through the incremental unwinding of an initially twisted saddle-shaped apoA-I double belt structure.

The helix bundle: a reversible lipid binding motif

Apolipoproteins are the protein components of lipoproteins that have the innate ability to inter convert between a lipid-free and a lipid-bound form in a facile manner, a remarkable property conferred by the helix bundle motif. Composed of a series of four or five amphipathic alpha-helices that fold to form a helix bundle, this motif allows the en face orientation of the hydrophobic faces of the alpha-helices in the protein interior in the lipid-free state. A conformational switch then permits helix-helix interactions to be substituted by helix-lipid interactions upon lipid binding interaction. This review compares the apolipoprotein high-resolution structures and the factors that trigger this switch in insect apolipophorin III and the mammalian apolipoproteins, apolipoprotein E and apolipoprotein A-I, pointing out the commonalities and key differences in the mode of lipid interaction. Further insights into the lipid-bound conformation of apolipoproteins are required to fully understand their functional role under physiological conditions.

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.