Andrea Catte

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.

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.

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.

Thermal Stability of Apolipoprotein A-I in High-Density Lipoproteins by Molecular Dynamics

Apolipoprotein (apo) A-I is an unusually flexible protein whose lipid-associated structure is poorly understood. Thermal denaturation, which is used to measure the global helix stability of high-density lipoprotein (HDL)-associated apoA-I, provides no information about local helix stability. Here we report the use of temperature jump molecular dynamics (MD) simulations to scan the per-residue helix stability of apoA-I in phospholipid-rich HDL. When three 20 ns MD simulations were performed at 500 K on each of two particles created by MD simulations at 310 K, bilayers remained intact but expanded by 40%, and total apoA-I helicity decreased from 95% to 72%. Of significance, the conformations of the overlapping N- and C-terminal domains of apoA-I in the particles were unusually mobile, exposing hydrocarbon regions of the phospholipid to solvent; a lack of buried interhelical salt bridges in the terminal domains correlated with increased mobility. Nondenaturing gradient gels show that 40% expansion of the phospholipid surface of 100:2 particles by addition of palmitoyloleoylphosphatidylcholine exceeds the threshold of particle stability. As a unifying hypothesis, we propose that the terminal domains of apoA-I are phospholipid concentration-sensitive molecular triggers for fusion/remodeling of HDL particles. Since HDL remodeling is necessary for cholesterol transport, our model for remodeling has substantial biomedical implications.

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.

Surface Density-Induced Pleating of a Lipid Monolayer Drives Nascent High-Density Lipoprotein Assembly

Biogenesis of high-density lipoproteins (HDL) is coupled to the transmembrane protein, ATP-binding cassette transporter A1 (ABCA1), which transports phospholipid (PL) from the inner to the outer membrane monolayer. Using a combination of computational and experimental approaches, we show that increased outer lipid monolayer surface density, driven by excess PL or membrane insertion of amphipathic helices, results in pleating of the outer monolayer to form membrane-attached discoidal bilayers. Apolipoprotein (apo)A-I accelerates and stabilizes the pleats. In the absence of apoA-I, pleats collapse to form vesicles. These results mimic cells overexpressing ABCA1 that, in the absence of apoA-I, form and release vesicles. We conclude that the basic driving force for nascent discoidal HDL assembly is a PL pump-induced surface density increase that produces lipid monolayer pleating. We then argue that ABCA1 forms an extracellular reservoir containing an isolated pressurized lipid monolayer decoupled from the transbilayer density buffering of cholesterol.

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.

Multinuclear NMR Investigation of the NaDNA/Ethidium Bromide Anisotropic System

Abstract A combined use of 31P, 23Na, 2H and 17O NMR spectroscopies and polarized light microscopy has been employed to investigate the effect of the ethidium bromide (EB) binding on the liquid crystalline phase of concentrated double stranded DNA solutions. The optical textures and the 31P and 23Na NMR spectra of the DNA anisotropic solutions show that the intercalation of EB induces significant modifications either in the arrangements of the DNA rods and the surrounding ionic atmosphere. On the contrary, no indication of significant changes of the orientational order of the water molecules around DNA emerges from the water 2H and 17O NMR spectra.

Direct prediction of EPR spectra from lipid bilayers : understanding structure and dynamics in biological membranes.

Abstract Of the many biophysical techniques now being brought to bear on studies of membranes, electron paramagnetic resonance (EPR) of nitroxide spin probes was the first to provide information about both mobility and ordering in lipid membranes. Here, we report the first prediction of variable temperature EPR spectra of model lipid bilayers in the presence and absence of cholesterol from the results of large scale fully atomistic molecular dynamics (MD) simulations. Three types of structurally different spin probes were employed in order to study different parts of the bilayer. Our results demonstrate very good agreement with experiment and thus confirm the accuracy of the latest lipid force fields. The atomic resolution of the simulations allows the interpretation of the molecular motions and interactions in terms of their impact on the sensitive EPR line shapes. Direct versus indirect effects of cholesterol on the dynamics of spin probes are analysed. Given the complexity of structural organisation in lipid bilayers, the advantage of using a combined MD‐EPR simulation approach is two‐fold. Firstly, prediction of EPR line shapes directly from MD trajectories of actual phospholipid structures allows unambiguous interpretation of EPR spectra of biological membranes in terms of complex motions. Secondly, such an approach provides an ultimate test bed for the up‐to‐date MD simulation models employed in the studies of biological membranes, an area that currently attracts great attention.

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.