PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of jlrJournal of Lipid Research
 
J Lipid Res. 2012 September; 53(9): 1851–1863.
PMCID: PMC3413226

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

Abstract

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.

Keywords: apolipoprotein A-I, all-atom simulation, Molecular Dynamics, X-ray crystallography, α11/3 helix

HDL is a supramolecular assembly of lipid and protein that represents a nanoscale “soft” form of condensed matter easily deformable by thermal fluctuations (1). Because the “soft” lipid environment in HDL modulates apoA-I conformation in dynamic ways, a detailed understanding of HDL structure-function requires new approaches.

The mature sequence of human apoA-I has a length of 243 amino acid residues and is encoded by exons 3 and 4 of the apolipoprotein gene in chromosome 11. The common lipid-associating motif in apoA-I is the amphipathic α helix (2, 3). Exon 3 encodes for residues 1–43 in human apoA-I. This domain, commonly referred to as the globular domain (4), contains an N-terminal segment of 10 residues, followed by three 11-mer amino acid tandem repeats, designated G1, G2, and G3, respectively.

Barker and Dayhoff (5), Fitch (6), and McLachlan (7) independently observed that apoA-I contains multiple repeats of 22 amino acids, each of which is a tandem array of two 11-mers. DNA sequencing work confirmed the existence of a 22-mer periodicity in apoA-I (8). The portion of apoA-I that follows the globular domain, residues 44–241 in human apoA-I, referred to as the lipid-associating domain (4), contains eight 22-mer and two 11-mer tandem repeats (called helixes 1–10), each of which has the helical wheel signature of an amphipathic α helix. The majority of the tandem helical repeats are punctuated by prolines. The fundamental unit of repeat in exon 4 is not an 11-mer but a 22-mer that is made up of two 11-mers (6, 8, 9). It has been suggested that the different apolipoprotein genes arose from a common ancestral gene by gene duplications of 11 residue codons (10).

Both negative-stain (11) and cryo electron (12) microscopy of apoA-I reconstituted with phospholipid reveal discoidal HDL particles the thickness of a phospholipid bilayer that are similar to nascent HDL particles (13). The first tangible experimental evidence for the conformation of apoA-I on the disc edge was determination of a 4 Å resolution solution-phase X-ray structure for residues 44–243 of lipid-free apoA-I reported by Borhani et al. in 1997 (14) that suggested an antiparallel double-belt model. In the first experimental test of the belt model, Koppaka et al. (15) used polarized attenuated total internal reflection Fourier-transform infrared spectroscopy to conclude that the result unambiguously supported a belt model. Theoretical considerations of the geometric and physical chemical nature of the apoA-I double-belt model resulted in the publication by our lab of an atomic resolution antiparallel double-belt amphipathic helical model for discoidal HDL (16) with a helix–helix registry termed LL5/5 (implying juxtaposition of pairwise helix 5 antiparallel repeats) surrounding the edge of a bilayer disc. This model provided sound theoretical rationale for use of the lipid-free X-ray crystal structure as a valid model for lipid-associated apoA-I. Five laboratories subsequently have studied reconstituted discoidal HDL using a variety of physical chemical methods, and the results have been consistent with our double-belt model (1721).

One approach to coping with the dynamic nature of HDL is the use of molecular dynamics (MD) that produces a computer simulation of physical movements of atoms and molecules. The atoms and molecules are allowed to interact for a period of time, giving a view of the motion of the atoms. MD uses high-performance computing to predict successive atomic positions and velocities of a molecule or molecular assembly by numerically solving Newton's laws of motion (F = ma). The hydrophobic effect results naturally from the statistical thermodyamics of the probability of ordering of water molecules at nonpolar interaces versus the much higher probability of water forming hydrogen bonds in every direction. A number of computer alogrithms have been developed [e.g., NAMD (22) and GROMACS (23)] that generate a trajectory of atomic motion over time that acts upon atoms and bonds. Forces between the particles and potential energy are defined by molecular mechanics force fields. Temperature is simulated by random assignment of a Boltzmann velocity distribution, equivalent to the desired temperature, to all of the atoms of the simulation. Because electrostatic interactions are accounted for by using fixed partial charges, chemical reactions involving chemical bonds being created or broken cannot be simulated (for a review, see Ref. 24).

Although conventional wisdom has been that time constraints and potential energy barriers are major obstacles to the use of MD simulations for determining protein conformations, we hypothesized that four profound constraints would be imposed on the conformations of lipid-associated apolipoproteins by lipid: the low dielectric environment, the hydrophobic effect, the amphipathic helix, and the planar bilayer environment (16, 25), such that MD simulations of HDL would produce useful structural information. Accordingly, by application of MD simulations to our detailed antiparallel double-belt model for discoidal HDL (16), we developed atomistic models of the dynamic interactions of apoA-I with itself and the “soft” lipid components of discoidal HDL (2531).

To increase exploration of conformational space of these HDL particles through all-atom MD, we developed a protocol that involves the use of temperature jumps to 500 K and unrestrained simulated annealing, approaches generally frowned upon by the MD community. Why did we think that MD combined with temperature jumps and/or unrestrained simulated annealing, generally considered to be inadequate for determination of globular protein structure, was capable of determining HDL structure? The answer, we hypothesized, was in the constraints noted above imposed by lipid on conformations of lipid-associated proteins (16, 25), constraints that are lacking in globular proteins.

At the time that our various models were published (2531), no high-resolution structural information on apoA-I was available for comparison. This situation changed with the recent publication of a high-resolution crystal structure of C-terminally truncated lipid-free apoA-I by Mei and Atkinson (32). Examination of this structure confirms many of the conformational features of apoA-I discoidal HDL previously proposed by us using modeling and MD simulations. These conformational features include structural details of the geometry of the double-belt model (16, 26, 28, 31), the existence of the α11/3 helix structural motif (16), and three of our previously reported key MD findings with biological relevance for the double-belt model for discoidal HDL, the LCAT presentation tunnel (26), the importance of interchain salt bridges in double-helix registration (16, 25, 28, 33), and the presence of a ”sticky” N-terminus (26, 30).

METHODS

The all-atom MD and molecular dynamics-simulated annealing (MDSA) protocol were performed using NAMD (34) as described (25, 27). Each system was solvated with the solvate plug-in of Visual Molecular Dynamics (VMD) (35).The TIP3P water model was used (36). Each system was also ionized and charge-neutralized with NaCl to 0.15 M with the autoionize plug-in of VMD. The CHARMM 22 (37, 38) and 27 (39, 40) force fields were used for protein and lipid molecules, respectively. For simulations at a fixed temperature, velocity reassignments occurred every 1 ns to prevent the “flying ice cube” effect (41). The MDSA protocol consisted of 10 ns equilibration at 500 K, 10 ns nonlinearly cooling to 310 K, and 10 ns equilibration at 310 K (2630). The initial MD at 310 K followed by the MDSA protocol was replicated to give an ensemble of 16 trajectories and final structures of discoidal HDL particles with POPC and unesterified cholesterol (UC) with stoichiometries of POPC:UC:apoA-I of 160:24:2 (26, 28, 30).

Other simulations previously published and used in this paper are: i) an ensemble of six discoidal particles with POPC:(Δ1-40)apoA-I stoichiometry of 100:2 subjected to temperature jump MD (500 K), which expanded the disc diameter to a size comparable to 160:2 (25, 27); ii) an ensemble of two discoidal particles with POPC:UC:apoA-I stoichiometry of 240:32:2 subjected to temperature jump MD (500 K) (26, 28, 30); and iii) an ensemble of four discoidal particles with the same stoichiometry as the ensemble of sixteen but simulated prior to that larger set (26, 28, 30). These two sets are also referred to as an ensemble of twenty 160:24:2 particles (26, 28, 30).

RESULTS

Comparison of apoA-I in the crystal structure of the C-terminally truncated lipid-free apoA-I to apoA-I in our discoidal HDL models and simulations demonstrates certain similar or identical conformational features of both crystal and simulated antiparallel apoA-I dimers, several that were observed originally in our MD simulations and proposed to have important biological functions (16, 25, 26, 28).

Double-belt model for apoA-I in discoidal HDL

Although multiple laboratories (1721) have used a variety of physical chemical methods that support our original detailed double-belt model for discoidal HDL, no direct structural data, other than the low-resolution X-ray crystal structure for (Δ1-43)apoA-I published by Borhani et al. (14), have been available to provide atomic resolution information about specific features of the double belt. The Mei and Atkinson structure (32) provides high-resolution atomic resolution information applicable for half of the double-belt model discoidal HDL model (16). As shown in Fig. 1, residues 79–178 in the unit cell dimer AB create a hemi circle of continuous antiparallel double helical domains. This portion of the crystal structure represents the complete pairwise antiparallel helix 4/6, 5/5, and 6/4 repeats and portions of the antiparallel helix 2-3/7 and 7/3-2 pairs. The diameter of this half circle is 105–106 Å, a value very close to diameters of the major reconstituted discoidal HDL particle created by disc formation from dimyristoylphospha-tidylcholine (DMPC) and POPC.

Fig. 1.
Molecular graphics cartoon image of the antiparallel double-belt continuous helical structure formed by residues 78–178 in the crystal structure of (Δ185–243)apoA-I (32). Helix 5, green; prolines, yellow; residues E78, space-filling. ...

We showed by gradient gel electrophoresis (GGE) that the most-common reconstituted DMPC:apoA-I particle, which we called R2-2 (termed the 96 Å particle by other groups), reaches a constant Stokes diameter of 105.5 (± 0.7) Å (n = 6) at a DMPC:apoA-I ratio of 180:2 and that a comparable reconstituted POPC:apoA-I particle, also called R2-2, reaches a slightly smaller constant Stokes diameter of 102 Å at a POPC:apoA-I ratio of 160:2 (28, 33).

Simple geometric calculations using the formula for the area of a circle, A = πr2, support these two particle sizes. For DMPC, assuming a molecular surface area of 63 Å2 (42), a bilayer disc containing 180 DMPC has a diameter of 85.0 Å. Adding the thickness of an α helix (10 Å) on each side, the full DMPC discoidal HDL particle has a diameter of 105.0 Å., close to the diameter measured by GGE of 105.5 Å (28, 33). For POPC, assuming a molecular surface area of 65 Å2 (28), a bilayer disc containing 160 POPC has a diameter of 81.4 Å. Adding the thickness of an α helix (10 Å) on each side, the full POPC discoidal HDL particle has a diameter of 101.4 Å, also close to the diameter measured by GGE of 102 Å (28).

Modified ideal α helix with a pitch of three turns per 11 residues (α11/3 helix)

In 1999, we introduced the concept of the α11/3 amphipathic helix. We proposed this slight conformational variation on the ideal α18/5 helical motif as a mechanism to create an amphipathic helical ring with a continuous hydrophobic edge on the inside of the ring, ideal for encirclement of the hydrophobic edge of a bilayer disc (16). We further suggested that only the central residues of the lipid-associating domain, helixes 2–9, residues 66–220, formed an α11/3 amphipathic helical conformation (16).

This conformation was used to create all our starting discoidal HDL models for MD simulations. However, this slight variation on the α helix structure is almost indistinguishable from an idealized α helix and has been difficult to confirm, even in the final MD-simulated discoidal structures, because of the curvature of the helical segments around the disc.

We examined the high-resolution crystal structure of (Δ185–243)apoA-I (32) for evidence of the α11/3 helix. In Fig. 2A, we present two molecular graphic views of a straight 12 residue α helical segment (residues 144–155; Fig. 1) of helix 6 for comparison with plots of the sequence in the form of α11/3 and α18/5 helical wheel diagrams (Fig. 2B, C, respectively). We avoided inclusion of P143 because of the obvious helix bends induced in the double belt by prolines (Fig. 1). Although there are slight irregularities in the α helicity of residues 144–155, this 12-residue segment clearly has an α11/3-helical pitch; note that residues 1 and 12 lie in the same plane along the long axis of the helix.

Fig. 2.
Comparison of the helical pitch of residues 144–155 in helix 6 of the crystal structure of (Δ185–243)apoA-I (32) with wheel diagrams of α11/3 and α18/5 helixes. A: Upper panel, backbone (solid gray lines) and transparent ...

Presentation tunnel hypothesis

Our MD simulations of discoidal HDL consistently result in a gap forming between the pairwise helix 5 repeats (26, 31) of the antiparallel molecular double-belt structure (16), a gap not present in the idealized starting structure with extended side chains (Fig. 3). The unique and consistent conformation of the residues forming the gap suggested to us a molecular basis for the activation by lipid-associated apoA-I of the plasma enzyme LCAT: after attachment of LCAT to discoidal HDL, the pairwise helix 5 repeats in apoA-I create an amphipathic “presentation” tunnel for migration of hydrophobic acyl chains and amphipathic UC from the bilayer to the phospholipase A2-like and esterification active sites of LCAT, respectively (26, 31). Examination of our ensemble of 16 MDSA simulations (26, 30, 43) shows that, in 14 of the 16, the methyl ends of POPC acyl chains (POPC-Me) from the bilayer center are inserted into the presentation tunnel and exposed to solvent (Fig. 4A). In the other two examples of the ensemble, the hydroxyl moiety of UC (UC-OH) is inserted into the presentation tunnel and exposed to solvent (Fig. 4B).

Fig. 3.
Comparison of the pairwise helix 5 repeats of the starting perfectly circular double-belt model used for the 16 MDSA simulations of discoidal HDL (30). A: Cross-eyed stereo of the helix 5 pair in the apoA-I double belt used for all 16 all-atom MDSA simulations ...
Fig. 4.
Representative examples of MD simulations showing lipids in the amphipathic presentation tunnel of the pairwise helix 5 domain. A: Cross-eyed stereo of the exposure to solvent in one of the 16 members of the 160:24:2 MDSA ensemble of a POPC-Me at the ...

Our earlier MD simulations of discoidal HDL composed of POPC without UC show that the presentation tunnel forms even in the absence of UC (27). The representative example shown in Fig. 4C of a presentation tunnel formed in discoidal HDL particles containing POPC without UC illustrates the key conformational features of the presentation tunnel (26, 31): i) two clusters of salt bridges form between residues E125, K133, E136, and K140 that pull the two K133 acyl chains away from the center of the pairwise helix 5 repeats to create a gap; ii) the low-bulk residues G129 and A130, located between the two clusters of salt bridges enlarge the gap or tunnel in the pairwise helix 5 repeats; larger residues would obstruct the gap; iii) the resulting tunnel is lined on its long edges by K133 and E125 and on its narrow edges by G129 and A130; G129 faces the solvent, A130 the lipid; iv) in essentially every simulation of discoidal HDL, either one (Fig. 4A) or two (Fig. 4C) terminal methyl groups or UC (Fig. 4B) are inserted into the tunnel gap (26).

Comparison of the crystal structure of C-terminally truncated lipid-free apoA-I (32) in Fig. 5A, D with a representative MD structure of discoidal HDL in Fig. 5B, C and 5E, F shows that both structures form salt bridges between residues K133, E136, and K140, and in both structures, the low-bulk residues G129 and A130 (26) are arranged between the salt bridges to create a gap or tunnel between the pairwise helix 5 repeats; this similarity was noted by Mei and Atkinson (32).

Fig. 5.
Comparison of the pairwise helix 5 domain of the recent crystal structure of C-terminally truncated apoA-I (32) with a representative example of 16 MD simulations of discoidal HDL (30). A: Face-on cross-eyed stereo view of the C-terminally truncated (Δ185–243)apoA-I ...

There is one pronounced difference between the protein crystal and the discoidal HDL simulation structures regarding pairwise helix 5 interactions: E125 in the crystal structure forms no interchain salt bridges with K133, whereas E125–K133 interchain salt bridges are generally formed in the discoidal HDL simulation structures (compare Fig. 5A, B). The reason is simple: in the antiparallel double helical protein-alone crystal structure, the hydrophobic faces of interchain pairwise helix 5 repeats remove themselves from solvent by associating with one another (Fig. 6A), but they rotate ~50° around the long axis of each helix to associate with the edges of the discoidal lipid bilayer (Fig. 6B). This rotation has several effects: i) A130, located in wheel position 10 in the center of the hydrophobic face (16, 44), rotates out of the tunnel gap to point directly at the lipid, and G129 rotates into the gap (compare Fig. 6A, B); this rotation results in an elongation of the presentation tunnel (compare Fig. 5D, E); ii) E125, K133, E136, and K140 on antiparallel helixes rotate toward one another (Fig. 6), allowing E125 to form interchain salt bridges with K133; iii) L126 (starred residues in Fig. 5A) rotates out of the gap and faces the lipid; this consensus conformation, two salt bridge clusters separated by low-bulk residues, that creates and stabilizes the tunnel is convincingly illustrated by molecular alignment of the pairwise helix 5 domains of the 16 simulations of the MD ensemble shown in Fig. 7.

Fig. 6.
Stick figure cross-eyed stereo comparisons of the relative rotations of the pairwise helix 5 repeats in the C-terminally truncated (Δ185–243)apoA-I dimeric crystal structure (32) and in a representative example of discoidal HDL subjected ...
Fig. 7.
Alignment of pairwise helix 5 domains (residues 121–143) of the apoA-I double belt from the 16 simulations of the MDSA ensemble of discoidal HDL (30). A: View of polar (solvent exposed) face of aligned pairwise helix 5 domains. B: View of nonpolar ...

Stabilization of LL5/5 chain registry by salt bridges

Another finding from our MD simulations of discoidal HDL confirmed by the crystal structure is formation of wp2↔wp2 (wp refers to α11/3 wheel positions within each 11-mer tandem repeat) interchain “solvent-inaccessible” salt bridges (16, 25, 28), bonds strengthened by the low dielectric environment (45) of the lipid matrix of the bilayer disc that we suggested represented stabilizing elements for the LL5/5 double-belt registry. Comparison of side views of the crystal and MD structures (Fig. 8A, B, respectively) shows that two of the three pairs of “solvent-inaccessible” salt bridges (E111↔H155, D89↔R177) seen in MD simulations are also seen in the crystal structure, noted by Mei and Atkinson (32). The third pair, E78↔R188, is not present in the crystal structure because of C-terminal truncation at residue N184 (32). The high frequency of formation of the first two pairs of “solvent-inaccessible” salt bridges is indicated by alignment of the 16 simulations of the MD ensemble between residues 78–120 and 143–188 (Fig. 8C). The third salt bridge pair, E78↔R188, in a more-mobile region of the double belt, forms less frequently (for example, it has not formed in Fig. 8B). As noted elsewhere (27, 28), this terminal pair of “solvent-inaccessible” salt bridges, E78↔R188, unlike the two more-central ones, are not only often broken during MD simulations but are poorly conserved evolutionarily (44).

Fig. 8.
Validation of proposed solvent salt bridge alignment of apoA-I double belt in the LL5/5 registry. A: Cross-eyed stereo side view of residues 79–178 of the C-terminally truncated (Δ185–243)apoA-I dimeric crystal structure (32). ...

We have also proposed that wp5↔wp9 interchain solvent-accessible salt bridges play a role in interchain stability (16, 25, 28). Residues putatively involved in formation of solvent-accessible salt bridges are also shown in Fig. 8A, B. Moving circumferentially from the central pairwise helix 5 repeats, no wp5↔wp9 interchain salt bridges are found in the crystal structure until E169↔K96 and R173↔D92 at the beginning of the pairwise helix 7/3 junctions. We have found that solvent-accessible interchain salt bridges form less readily (16, 25, 28) than “solvent-inaccessible” salt bridges during MD simulations (for example, E162↔K103 has not formed in the example shown in Fig. 8B).

It is clear from examination of the salt bridge cluster, E169↔K96, R173↔D92, and E78↔H177, present in both the crystal and the MD simulations shown in Fig. 8A, B, that E78↔H177 is behind the cluster E169↔K96 and R173↔D92 and thus less accessible to solvent. A reasonable question, however, is whether “solvent-inaccessible” salt bridges should be considered “solvent-inaccessible” or simply “less-solvent-accessible” than those we term “solvent accessible.”

Because we have shown that the pairwise helix 5 repeats in the MD-simulated discoidal HDL relative to the crystal structure are rotated ~50° toward the discoidal bilayer, we looked at the rotational positioning of the “solvent-inaccessible” salt bridges in the crystal structure versus those in MD-simulated discoidal HDL (Fig. 9). Surprisingly, in the helix pairwise helix 4/6 repeats, there is very little difference between the rotational alignments of the helical pairs in the crystal structure versus the MD-simulated discoidal HDL. In Fig. 9, the Cα atoms of wp10 are indicated by magenta spheres. The rotation of helix 6 is slightly different between the crystal and the simulation; in the simulation, helix 6 is rotated 15–20° toward the discoidal bilayer compared with helix 6 in the crystal. In any case, the relative rotations of the “solvent-inaccessible” salt bridge (E111↔H155) in the crystal structure versus that in the MD-simulated discoidal HDL are similar (Fig. 9).

Fig. 9.
Backbone cross-eyed stereo comparisons of the relative rotations of the pairwise helix 4/6 segments in the C-terminally truncated (Δ185–243)apoA-I dimeric crystal structure (32) and a representative example of discoidal HDL subjected to ...

The tertiary structure of the crystal dimer can explain the difference in relative rotational alignments of the pairwise helix 5 and helix 4/6 repeats. Only the helix 5 repeats are completely exposed to solvent; the helix 4/6 repeats form a “four-helix” bundle with the N-terminal domains (32), indicated by two circles in Fig. 9A, that partially mimic the effects of a lipid bilayer.

Figure 9 also provides evidence of the existence of the α11/3 helical pitch in the simulation as well as the crystal structures. In Fig. 9, the wp2 residue pairs present in helix 4 (Y100 and E111) and helix 6 (L144 and H155) are indicated by arrows. In both helixes, in both the crystal and simulation structures, the wp2 Cα atoms lie on a plane through the long axis of the helix, unequivocally supporting the existence of the α11/3 in these two tandem helical segments in both crystal and simulation structures.

Molecular stickiness of the N-terminal domain

Finally we showed from MD simulations (27, 28, 30) that the N-terminal domains (residues 1–22) act as a form of molecular “Velcro” to avidly stick to their N-terminal companions in the double belt (Fig. 5E) as well as other portions of apoA-I (30), a property that we hypothesize stabilizes multiple conformations of apoA-I on HDL (30). The stickiness of N-terminal pairs during MD simulations appears to be largely driven by strategically placed moieties primed to form salt bridges (27, 28, 30), particularly the acidic D1 and E2 and the basic R10 and N-terminal −NH3+ (Fig. 10). The three residues, D1, E2, and R10 are highly conserved, as are many of the first 22 residues of the N-terminal domain (44).

Fig.10.
Cross-eyed stereo image of salt bridge formation between acidic and basic moieties of the N-terminal 10 residue pairs of antiparallel double-belt apoA-I in a discoidal particle after 6 ns of MD simulation (30). This image was the result of a high-temperature ...

Although the two N-terminal residues D1 and E2, because of excess motion, are not resolved in the crystal structure (32), the Velcro hypothesis is consistent with analysis of the X-ray data (Fig. 11). The N-terminal domains of each apoA-I dimer, e.g., chain A in unit cell AB, contact adjacent apoA-I dimers, e.g., chain D in unit cell CD, at the juncture of helixes 6 and 7 (Fig. 11A). Chain D contains two highly conserved basic residues, R160 and R171 (30), and one highly conserved acidic residue, D168 (30), that are extended toward and into close proximity with the last resolved N-terminal residue in chain A, P3 (Fig. 11A, B). Residues R160, R171, and D168 in chain D are 9.8 Å, 10.4 Å, and 6.8 Å from P3 in chain A, respectively. Based upon the geometry of Fig. 11, this conformational proximity suggests that interaction of chain A with chain D in Fig. 11B involves interaction of the N-terminal region of chain A with residues R160, R171, and perhaps D168 of chain D. This interaction is almost certainly through salt bridge formation between the terminal two residues of D1 and E2 of chain A, and perhaps its N-terminal −NH3+ with residues R160, R171, and perhaps D168 of chain D. The lack of resolution of D1 and E2 in the crystal structure and the distances of R160 and R171 from P3 are consistent with alternate salt bridge pairing of D1/E2 with R160/R171 that would introduce what would appear to be motion in the terminal residues of chain A. Modeling (data not shown) shows that for D1 to salt bridge alternately with R160 and R171, residue Q5 must make a substantial rotation about its [var phi], ψ bonds. The b-factor analysis of Cα atoms in chain A illustrated in Fig. 11C indicates that consistent with alternate salt bridge pairing and the sticky N-terminal hypothesis, Q5 is more mobile (yellow) than adjacent residues P4 and S6 (both green), whereas P3 is more mobile than all three (red).

Fig.11.
Cross-eyed stereo images of the interaction of the N-terminus of 1 apoA-I in 1 U cell dimer (chain A) with the helix 6-helix 7 junction of another apoA-I in a second unit cell dimer (chain D) from the C-terminally truncated (Δ185–243)apoA-I ...

DISCUSSION

At present, no direct experimental methods exist for determining atomic resolution structures of biologically important nanoscale particles of soft condensed matter, such as intact HDL and other lipoproteins. The similarity reported here between key conformational elements from the crystal structure of lipid-free apoA-I and apoA-I in an ensemble of discoidal HDL subjected to MD simulation (25, 28) supports the proposition that virtual experimentation by MD simulation provides a reasonably robust method for determining details of the structure of soft protein-lipid assemblies such as discoidal HDL particles. Our simulations of discoidal HDL were successful because we developed a detailed model for the starting protein registry of the antiparallel apoA-I chains, the detailed antiparallel double-belt model (16). The studies reported here confirm two elements of that starting model: i) a lowest-free-energy particle diameter of 105–106 Å and ii) an α11/3 helical pitch. If we had started with a random registry and an ideal α18/5 helical pitch, convergence to the LL5/5 registry would have been beyond current computational capabilities. It is likely that biosynthesis of nascent discoidal HDL particles will depend critically upon self-registry of the lipid-free apoA-I (46). If this is true, the C-terminally truncated crystal structure of Mei and Atkinson (32) will prove useful in elucidating the structure of circulating lipid-free/lipid-poor (pre-β-HDL) apoA-I.

In the crystal structure, residues 79–178 in the independent dimer AB create essentially a perfect half circle of continuous antiparallel double helical domains with a diameter of 105–106 Å (Fig. 1). This is the portion of the apoA-I double belt in discoidal HDL that remains relatively immobile during MD simulations (30) and is largely responsible for determining the diameter of the discoidal particle. This feature: i) supports the relevance of the crystal structure for discoidal HDL; ii) supports our MD simulation results; and iii) suggests that the conformation of residues 78–178 in the Mei-Atkinson crystal reflects the lowest-free-energy state for apoA-I in a reconstituted discoidal HDL structure. Supporting the lowest-free-energy concept, the diameter of this most-prevalent particle varies only 3 Å between reconstituted discoidal HDL made from two very different phospholipids, DMPC and POPC (28).

Many other laboratories call this most-prevalent of reconstituted discoidal HDL particles the 96 Å particle (19, 47, 48). Because this particle varies in diameter with both the lipid:protein molar ratio and the type of phospholipid, we prefer the less-ambiguous nomenclature R2-2 (30, 33). We also believe that 96 Å is not the correct diameter of this particle, for reasons we have discussed elsewhere (33). The simple geometric calculations in the RESULTS section using the formula for the area of a circle, A = πr2, giving a diameter of 102–111 Å, depending upon the phospholipid, also support this conclusion.

In none of our simulations do we see evidence for the Belt-Buckle model of Bhat et al. (47). Because of the unusual stickiness of the N-terminal domains of the double-belt apoA-I to one another that we see in our MD simulations (30), except when particles are too large to maintain N-terminal/N-terminal salt bridge connections (30), we never see the kind of movement of the N-terminal domain that these authors suggest. It could, of course, be that we do not explore sufficient conformational space in the time frame available to our MD simulations. However, we, are in favor of another explanation: attachment of a long cross-linking reagent to the N-terminal end of apoA-I on the discs creates serious artifacts: i) attachment increases the flexibility of the N-terminus, resulting in cross-linking to lysine residues outside the actual cross-linking distance; ii) some of the cross-links are to the −NH3+ moiety of the N-terminus itself, and because that basic moiety is part of the sticky salt bridge network, as is K12, N-terminal stickiness is destroyed, and that portion flops around to cross-linking to lysine residues outside the actual cross-linking distance; iii) or both.

We would like to comment on the previously published crystal structure for untruncated apoA-I (49). There are substantial reasons to think that this is an invalid structure. For details, see http://main.uab.edu/Sites/reporter/articles/71570/, “UAB's Statement on Protein Data Bank Issues.”

We proposed the α11/3 helical pitch, a slight variation on the ideal α18/5 helical motif, as a mechanism to create an amphipathic helical ring with a continuous hydrophobic edge on the inside of the ring, ideal for encirclement of the hydrophobic edge of a bilayer disc (16). The α11/3 helix was mentioned by Linus Pauling in his original description of the α helix when he observed that the pitch of the α helix could vary from 3.6 (18/5) to 3.67 (11/3) residues per turn without significant energy penalty (50). The result presented in this paper unequivocally support the existence of the α11/3 helix in both crystal and simulation structures of apoA-I. The α11/3 amphipathic helix has been suggested for several lipid-associating proteins (apoA-I, apoA-IV, and α-synuclein), but the results presented here are the first solid structural evidence for this variant of the ideal α18/5 helix.

The results in this paper provide validation of the LCAT presentation tunnel hypothesis. Certainly both the crystal and MD simulation structures of apoA-I show that the pairwise helix 5 repeats are uniquely designed to create a gap that exposes both acyl chains and UC. This result of the N-terminally truncated apoA-I crystal structure fully validates the LCAT presentation tunnel hypothesis, has enormous biological implications for HDL function, and provides a powerful working model for confirmation by future in-solution experimental studies.

We show that the presence or absence of UC has no effect on the formation of the gap between the pairwise helix 5 repeats. The major difference between discoidal HDL with and without UC that we find in our MD simulations is discoidal bilayer flexibility and stability; UC-containing discs are less flexible and more stable than those without UC (data not shown).

The presence in the C-terminally truncated crystal structure of wp2↔wp2 interchain “solvent-inaccessible” salt bridges supports our hypothesis that these represent stabilizing elements for the LL5/5 double-belt registry (16, 25, 28). Some portions of the crystal structure also show the presence of interchain solvent-accessible salt bridges, suggesting that they too may play a role in LL5/5 interchain stabilization. Based upon the analysis shown in Fig. 9, in the future, we will refer to the wp2↔wp2 interchain salt bridges as less-solvent-accessable and the wp5↔wp9 interchain salt bridges as more-solvent-accessable.

Because a reasonable assumption is that helix 5 initiates the double-helix registry, the lack of any interchain salt bridges in the pairwise helix 5 repeats in the crystal structure requires discussion. We looked for other interchain interactions that could initiate the LL5/5 registry in the lipid-free apoA-I structure. We observed one pair of aromatic residues, Y100 and Y166, which might be close enough to forming π-π stacks, and juxtaposition of one pair of methionine residues, M112 and M148. None of these features of the crystal structure seemed necessary or sufficient to induce registry of the pairwise helix 5 repeats in the lipid-free state.

However, as noted by Mei and Atkinson (32), the pairwise helix 5 repeats form impressive leucine zipper motifs between opposing chains at both ends: L141/L122/L137/L126 and perhaps L134 (Fig. 12). These multiple residue leucine zipper motifs at opposite ends of the pairwise helix 5 repeats should be inconsequential once the domain is associated with lipid but seem highly likely to be the principal motif for aligning the pairwise antiparallel helix 5 repeats in the lipid-free state. Upon association with lipid, rotation of each pairwise helix 5 by ~50° would induce interchain salt bridges, create the presentation tunnel, and “zip-lock” the double belt into the LL5/5 rotamer registration.

Fig.12.
Cross-eyed stereo images of the pairwise helix 5 repeats showing the presence of four- to five-residue leucine zippers at both ends in the C-terminally truncated (Δ185–243)apoA-I dimeric crystal structure (32). Helix 5, cartoon green; ...

Finally, although MD simulations convince us of the likely importance of the N-terminal stickiness of apoA-I in stabilizing both discoidal and spheroidal HDL particles, because the most-important residues in induction of stickiness, D1 and E2, are not visible in the crystal structure, our arguments for evidence of stickiness are indirect and less than conclusive.

In conclusion, we have demonstrated that molecular modeling and MD simulations are able to provide accurate detailed structural information about soft protein-lipid assemblies such as HDL, even when high temperature jump MD simulated annealing is used to increase exploration of conformational space.

Acknowledgments

The authors thank UAB Information Technology's Research Computing group and the UAB Department of Mechanical Engineering for use of the high-performance computer cluster Cheaha that they jointly maintain.

Footnotes

Abbreviations:

DMPC
dimyristoylphosphatidylcholine
GGE
gradient gel electrophoresis
MD
molecular dynamics
MDSA
molecular dynamics-simulated annealing
POPC-Me
methyl ends of POPC acyl chains
RMSD
root mean squared deviation
UC
unesterified cholesterol
UC-OH
hydroxyl moiety of UC
VMD
Visual Molecular Dynamics
wp
α11/3 wheel positions within each 11-mer tandem repeat

This work was supported by National Institutes of Health Grants HL-34343 and HL-102515. Its contents are solely the responsibility of the authors and do not necessarily represent the official views of the National Institutes of Health or other granting agencies.

REFERENCES

1. Zlotnick A. 2004. Viruses and the physics of soft condensed matter. Proc. Natl. Acad. Sci. USA. 101: 15549–15550. [PubMed]
2. Segrest J. P., Jackson R. L., Morrisett J. D., Gotto A. M., Jr 1974. A molecular theory of lipid-protein interactions in the plasma lipoproteins. FEBS Lett. 38: 247–258. [PubMed]
3. Segrest J. P., Garber D. W., Brouillette C. G., Harvey S. C., Anantharamaiah G. M. 1994. The amphipathic alpha helix: a multifunctional structural motif in plasma apolipoproteins. Adv. Protein Chem. 45: 303–369. [PubMed]
4. Segrest J. P., Jones M. K., De Loof H., Brouillette C. G., Venkatachalapathi Y. V., Anantharamaiah G. M. 1992. The amphipathic helix in the exchangeable apolipoproteins: a review of secondary structure and function. J. Lipid Res. 33: 141–166. [PubMed]
5. Barker W. C., Dayhoff M. O. 1977. Evolution of lipoproteins deduced from protein sequence data. Comp. Biochem. Physiol. B. 57: 309–315. [PubMed]
6. Fitch W. M. 1977. Phylogenies constrained by the crossover process as illustrated by human hemoglobins and a thirteen-cycle, eleven-amino-acid repeat in human apolipoprotein A-I. Genetics. 86: 623–644. [PubMed]
7. McLachlan A. D. 1977. Repeated helical pattern in apolipoprotein-A-I. Nature. 267: 465–466. [PubMed]
8. Karathanasis S. K., Zannis V. I., Breslow J. L. 1983. Isolation and characterization of the human apolipoprotein A-I gene. Proc. Natl. Acad. Sci. USA. 80: 6147–6151. [PubMed]
9. Luo C. C., Li W. H., Moore M. N., Chan L. 1986. Structure and evolution of the apolipoprotein multigene family. J. Mol. Biol. 187: 325–340. [PubMed]
10. Li W. H., Tanimura M., Luo C. C., Datta S., Chan L. 1988. The apolipoprotein multigene family: biosynthesis, structure, structure-function relationships, and evolution. J. Lipid Res. 29: 245–271. [PubMed]
11. Tall A. R., Lange Y. 1978. Interaction of cholesterol, phospholipid and apoprotein in high density lipoprotein recombinants. Biochim. Biophys. Acta. 513: 185–197. [PubMed]
12. van Antwerpen R., Chen G. C., Pullinger C. R., Kane J. P., LaBelle M., Krauss R. M., Luna-Chavez C., Forte T. M., Gilkey J. C. 1997. Cryo-electron microscopy of low density lipoprotein and reconstituted discoidal high density lipoprotein: imaging of the apolipoprotein moiety. J. Lipid Res. 38: 659–669. [PubMed]
13. Dory L., Boquet L. M., Hamilton R. L., Sloop C. H., Roheim P. S. 1985. Heterogeneity of dog interstitial fluid (peripheral lymph) high density lipoproteins: implications for a role in reverse cholesterol transport. J. Lipid Res. 26: 519–527. [PubMed]
14. Borhani D. W., Rogers D. P., Engler J. A., Brouillette C. G. 1997. Crystal structure of truncated human apolipoprotein A-I suggests a lipid-bound conformation. Proc. Natl. Acad. Sci. USA. 94: 12291–12296. [PubMed]
15. Koppaka V., Silvestro L., Engler J. A., Brouillette C. G., Axelsen P. H. 1999. The structure of human lipoprotein A-I. Evidence for the “belt” model. J. Biol. Chem. 274: 14541–14544. [PubMed]
16. Segrest J. P., Jones M. K., Klon A. E., Sheldahl C. J., Hellinger M., De Loof H., Harvey S. C. 1999. A detailed molecular belt model for apolipoprotein A-I in discoidal high density lipoprotein. J. Biol. Chem. 274: 31755–31758. [PubMed]
17. Tricerri M. A., Behling Agree A. K., Sanchez S. A., Jonas A. 2000. Characterization of apolipoprotein A-I structure using a cysteine-specific fluorescence probe. Biochemistry. 39: 14682–14691. [PubMed]
18. Davidson W. S., Hilliard G. M. 2003. The spatial organization of apolipoprotein A-I on the edge of discoidal high density lipoprotein particles: a mass spectrometry study. J. Biol. Chem. 278: 27199–27207. [PubMed]
19. Martin D. D., Budamagunta M. S., Ryan R. O., Voss J. C., Oda M. N. 2006. Apolipoprotein A-I assumes a “looped belt” conformation on reconstituted high density lipoprotein. J. Biol. Chem. 281: 20418–20426. [PubMed]
20. Bhat S., Sorci-Thomas M. G., Tuladhar R., Samuel M. P., Thomas M. J. 2007. Conformational adaptation of apolipoprotein A-I to discretely sized phospholipid complexes. Biochemistry. 46: 7811–7821. [PMC free article] [PubMed]
21. Wu Z., Wagner M. A., Zheng L., Parks J. S., Shy J. M., III, Smith J. D., Gogonea V., Hazen S. L. 2007. The refined structure of nascent HDL reveals a key functional domain for particle maturation and dysfunction. Nat. Struct. Mol. Biol. 14: 861–868. [PubMed]
22. Phillips J. C., Braun R., Wang W., Gumbart J., Tajkhorshid E., Villa E., Chipot C., Skeel R. D., Kale L., Schulten K. 2005. Scalable molecular dynamics with NAMD. J. Comput. Chem. 26: 1781–1802. [PMC free article] [PubMed]
23. Hess B., Kutzner C., van der Spoel D., Lindahl E. 2008. GROMACS 4: algorithms for highly efficient, load-balanced, and scalable molecular simulation. J. Chem. Theory Comput. 4: 435–447. [PubMed]
24. Rueda M., Ferrer-Costa C., Meyer T., Perez A., Camps J., Hospital A., Gelpi J. L., Orozco M. 2007. A consensus view of protein dynamics. Proc. Natl. Acad. Sci. USA. 104: 796–801. [PubMed]
25. Catte A., Patterson J. C., Jones M. K., Jerome W. G., Bashtovyy D., Su Z., Gu F., Chen J., Aliste M. P., Harvey S. C., et al. 2006. Novel changes in discoidal high density lipoprotein morphology: a molecular dynamics study. Biophys. J. 90: 4345–4360. [PubMed]
26. Jones M. K., Catte A., Li L., Segrest J. P. 2009. Dynamics of activation of lecithin:cholesterol acyltransferase by apolipoprotein A-I. Biochemistry. 48: 11196–11210. [PMC free article] [PubMed]
27. Jones M. K., Catte A., Patterson J. C., Gu F., Chen J., Li L., Segrest J. P. 2009. Thermal stability of apolipoprotein A-I in high-density lipoproteins by molecular dynamics. Biophys. J. 96: 354–371. [PubMed]
28. Gu F., Jones M. K., Chen J., Patterson J. C., Catte A., Jerome W. G., Li L., Segrest J. P. 2010. Structures of discoidal high density lipoproteins: a combined computational-experimental approach. J. Biol. Chem. 285: 4652–4665. [PMC free article] [PubMed]
29. Jones M. K., Zhang L., Catte A., Li L., Oda M. N., Ren G., Segrest J. P. 2010. Assessment of the validity of the double superhelix model for reconstituted high density lipoproteins: a combined computational-experimental approach. J. Biol. Chem. 285: 41161–41171. [PMC free article] [PubMed]
30. Jones M. K., Gu F., Catte A., Li L., Segrest J. P. 2011. “Sticky” and “promiscuous,” the yin and yang of apolipoprotein A-I termini in discoidal high-density lipoproteins: a combined computational-experimental approach. Biochemistry. 50: 2249–2263. [PMC free article] [PubMed]
31. Li L., Li S., Jones M. K., Segrest J. P. 2012. Rotational and hinge dynamics of discoidal high density lipoproteins probed by interchain disulfide bond formation. Biochim. Biophys. Acta. 1821: 481–489. [PMC free article] [PubMed]
32. Mei X., Atkinson D. 2011. Crystal structure of C-terminal truncated apolipoprotein A-I reveals the assembly of high density lipoprotein (HDL) by dimerization. J. Biol. Chem. 286: 38570–38582. [PMC free article] [PubMed]
33. Li L., Chen J., Mishra V. K., Kurtz J. A., Cao D., Klon A. E., Harvey S. C., Anantharamaiah G. M., Segrest J. P. 2004. Double belt structure of discoidal high density lipoproteins: molecular basis for size heterogeneity. J. Mol. Biol. 343: 1293–1311. [PubMed]
34. Kale´ L., Skeel R., Bhandarkar M., Brunner R., Gursoy A., Krawetz N., Phillips J., Shinozaki A., Varadarajan K., Schulten K. 1999. NAMD2: greater scalability for parallel molecular dynamics. J. Comput. Phys. 151: 283–312.
35. Humphrey W., Dalke A., Schulten K. 1996. VMD: visual molecular dynamics. J. Mol. Graph. 14: 33–38, 27–28. [PubMed]
36. Jorgensen W. L., Chandrasekhar J., Madura J. D. 1983. Comparison of simple potential functions for simulating liquid water. J. Chem. Phys. 79: 926–935.
37. Brooks B. R., Bruccoleri R. E., Olafson B. D., States D. J., Swaminathan S., Karplus M. 1983. CHARMM: a program for macromolecular energy, minimization, and dynamics calculations. J. Comput. Chem. 4: 187–217.
38. MacKerell A. D., Jr, Bashford D., Bellot M., Dunbrack R. L., Jr, Evanseck J., Field M. J., Fischer S., Gao J., Guo H., Ha S., et al. 1998. All-atom empirical potential for molecular modeling and dynamics studies of proteins. J. Phys. Chem. B. 102: 3586–3616. [PubMed]
39. Schlenkrich M., Brickmann J., MacKerell A., Jr, Karplus M. 1996. An empirical potential energy function for phospholipids: criteria for parameter optimization and applications. In Biological Membranes: A Molecular Perspective from Computation and Experiment. K. M. Merz and B. Roux, editors. Birkhauser, Boston. 31–81.
40. Feller S. E., Pastor R. W. 1997. Length scales of lipid dynamics and molecular dynamics. Pac. Symp. Biocomput. 142–150.
41. Harvey S. C., Tan R. K. Z., Cheatham T. E. 1998. The flying ice cube: velocity rescaling in molecular dynamics leads to violation of energy equipartition. J. Comput. Chem. 19: 726–740.
42. Hogberg C. J., Lyubartsev A. P. 2006. A molecular dynamics investigation of the influence of hydration and temperature on structural and dynamical properties of a dimyristoylphosphatidylcholine bilayer. J. Phys. Chem. B. 110: 14326–14336. [PubMed]
43. Catte A., Patterson J. C., Bashtovyy D., Jones M. K., Gu F., Li L., Rampioni A., Sengupta D., Vuorela T., Niemela P., et al. 2008. Structure of spheroidal HDL particles revealed by combined atomistic and coarse-grained simulations. Biophys. J. 94: 2306–2319. [PubMed]
44. Bashtovyy D., Jones M. K., Anantharamaiah G. M., Segrest J. P. 2011. Sequence conservation of apolipoprotein A-I affords novel insights into HDL structure-function. J. Lipid Res. 52: 435–450. [PMC free article] [PubMed]
45. Blundell T., Barlow D., Borkakoti N., Thornton J. 1983. Solvent-induced distortions and the curvature of alpha-helices. Nature. 306: 281–283. [PubMed]
46. Jayaraman S., Cavigiolio G., Gursky O. 2012. Folded functional lipid-poor apolipoprotein A-I obtained by heating of high-density lipoproteins: relevance to high-density lipoprotein biogenesis. Biochem. J. 442: 703–712. [PMC free article] [PubMed]
47. Bhat S., Sorci-Thomas M. G., Alexander E. T., Samuel M. P., Thomas M. J. 2005. Intermolecular contact between globular N-terminal fold and C-terminal domain of apoA-I stabilizes its lipid-bound conformation: studies employing chemical cross-linking and mass spectrometry. J. Biol. Chem. 280: 33015–33025. [PubMed]
48. Davidson W. S., Thompson T. B. 2007. The structure of apolipoprotein A-I in high density lipoproteins. J. Biol. Chem. 282: 22249–22253. [PubMed]
49. Ajees A. A., Anantharamaiah G. M., Mishra V. K., Hussain M. M., Murthy H. M. 2006. Crystal structure of human apolipoprotein A-I: insights into its protective effect against cardiovascular diseases. Proc. Natl. Acad. Sci. USA. 103: 2126–2131. [PubMed]
50. Pauling L., Corey R. B. 1951. Atomic coordinates and structure factors for two helical configurations of polypeptide chains. Proc. Natl. Acad. Sci. USA. 37: 235–240. [PubMed]
51. Segrest J. P., Jones M. K., Mishra V. K., Anantharamaiah G. M. 2002. Experimental and computational studies of the interactions of amphipathic peptides with lipid surfaces. Current Topics in Membranes 52: 391–429.
52. Bernstein H. J. 2000. Recent changes to RasMol, recombining the variants. Trends Biochem. Sci. 25: 453–455. [PubMed]

Articles from Journal of Lipid Research are provided here courtesy of American Society for Biochemistry and Molecular Biology