Search

Over the last two decades, an increasing number of studies has been devoted to a deeper understanding of the molecular process involved in the binding of various agonists and antagonists to active and inactive conformations of β2-adrenergic receptor (β2 AR). The 3.2 A X-ray crystal structure of human β2 AR active-state in combination with the endogenous low affinity agonist adrenaline offers an ideal starting structure for studying the binding of various catecholamines to adrenergic receptors. We show that molecular docking of L-DOPA and Droxidopa into rigid and flexible β2 AR models leads for both ligands to binding anchor sites comparable to those experimentally reported for adrenaline, namely D113/N312 and S203/S204/S207 side chains. Both ligands have a hydrogen bond network that is extremely similar to those of noradrenaline and dopamine. Interestingly, re-docking neutral and protonated versions of adrenaline to rigid and flexible β2 AR models results in binding poses that are more energetically stable and distinct from the X-ray crystal structure. Similarly, lowest energy conformations of noradrenaline and dopamine generated by docking into flexible β2 AR models had binding free energies lower than those of best poses in rigid receptor models. Furthermore, our findings show that L-DOPA and Droxidopa molecules have binding affinities comparable to those predicted for adrenaline, noradrenaline, and dopamine, which are consistent with previous experimental and computational findings and supported by the MD simulations of β2AR-ligand complexes performed here. SIGNIFICANCE The β2-adrenergic receptor (β2 AR), which belongs to the vast family of Guanine nucleotide-binding protein coupled receptors (GPCRs), is a transmembrane protein that is activated by the catecholamines norepinephrine (noradrenaline) and epinephrine (adrenaline). Using molecular docking methods, we discovered that the exogenous ligands L-DOPA and Droxidopa have the same hydrogen bond binding sites as endogenous agonists such as adrenaline, noradrenaline, and dopamine. Furthermore, all of the catecholamines studied had different hydrophobic binding sites in the receptor. Our findings might help researchers to better understand how existing and new drugs chemically similar to Droxidopa, which is used to treat Parkinson's disease, interact with β2 AR.

Summary 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.

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.

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 A 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-A 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.

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.

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.