Carbon and nitrogen isotopic variability in foxtail millet (Setaria italica) with watering regime.

RATIONALE: Carbonised plant remains are analysed for reconstruction of past climates and agricultural regimes. Several recent studies have used C4 plants to address related questions, and correlations between modern C4 plant δ13 C values and rainfall have been found. The millets were important food crops in prehistoric Eurasia, yet little is known about causes of isotopic variation within millet species. Previous research has shown there to be significant isotopic variation between millet accessions. Here we compare isotope ratios from plants grown under different watering regimes. This allows for a consideration of whether or not Setaria italica is a good proxy for environmental reconstruction. METHODS: We compare stable isotope ratios of Setaria italica plants grown in a controlled environment chamber with different watering regimes. We compare the carbon isotope ratios of leaves and grains, and the nitrogen isotope ratios of grains, from 12 accessions of Setaria italica. RESULTS: We find significant isotopic variability between watering regimes. Carbon isotope ratios are positively correlated with water availability, and on average vary by 1.9‰ and 1.7‰ for leaves and grains, respectively. Grain nitrogen isotope ratios also vary with watering regime; however, the highest isotope ratios are found with the 130-mL watering regime. CONCLUSIONS: The carbon isotope ratios of Setaria italica are strongly correlated with water availability. However, the correlation is the opposite to that seen in studies of C3 plants. The difference in isotopic ratio due to watering regime is comparable with that seen between different accessions; thus distinguishing between changing varieties of Setaria italica and changing climate is problematic. In terms of grain nitrogen isotope ratios, the highest δ15 N values were not associated with the lowest watering regime. Again, δ15 N variation is comparable with that which would be expected from an aridity effect or a manuring effect, and thus distinguishing between these factors is probably problematic.

Cooked starchy food in hearths ca. 120 kya and 65 kya (MIS 5e and MIS 4) from Klasies River Cave, South Africa.

Plant carbohydrates currently constitute 55-80% of the modern human diet (FAO and WHO, 1997) and some of today's key global health issues are associated with excessive carbohydrate consumption. However, starch carbohydrate is still a poorly understood element of modern human diet and our past starch diet may provide insights for future research. Despite an archaeological narrative that links our early hominin ancestors to a diet that is rich in roots and tubers, there is little deep time archaeological evidence of human plant starch consumption. Geneticists hypothesise that the duplication of starch digestion genes in early Homo sapiens (∼300 kya), is an adaptive response to an increased starch diet. Here we offer the earliest evidence of identified fragments of charred starch plant tissue (parenchyma) from cave and rock shelter hearths dated to Marine Isotope Stage (MIS) 5e and MIS 4, from the Middle Stone Age (MSA) site of Klasies River main site, South Africa (34.06°S, 24.24°E).

An Evaluation of the Crystal Structure of C-terminal Truncated Apolipoprotein A-I in Solution Reveals Structural Dynamics Related to Lipid Binding.

Apolipoprotein (apo) A-I mediates many of the anti-atherogenic functions attributed to high density lipoprotein. Unfortunately, efforts toward a high resolution structure of full-length apoA-I have not been fruitful, although there have been successes with deletion mutants. Recently, a C-terminal truncation (apoA-I(Δ185-243)) was crystallized as a dimer. The structure showed two helical bundles connected by a long, curved pair of swapped helical domains. To compare this structure to that existing under solution conditions, we applied small angle x-ray scattering and isotope-assisted chemical cross-linking to apoA-I(Δ185-243) in its dimeric and monomeric forms. For the dimer, we found evidence for the shared domains and aspects of the N-terminal bundles, but not the molecular curvature seen in the crystal. We also found that the N-terminal bundles equilibrate between open and closed states. Interestingly, this movement is one of the transitions proposed during lipid binding. The monomer was consistent with a model in which the long shared helix doubles back onto the helical bundle. Combined with the crystal structure, these data offer an important starting point to understand the molecular details of high density lipoprotein biogenesis.

Intraspecific carbon and nitrogen isotopic variability in foxtail millet (Setaria italica).

RATIONALE: Isotopic palaeodietary studies generally focus on bone collagen from human and/or animal remains. While plant remains are rarely analysed, it is known that plant isotope values can vary as a result of numerous factors, including soil conditions, the environment and type of plant. The millets were important food crops in prehistoric Eurasia, yet little is known about the isotopic differences within millet species. METHODS: Here we compare the stable isotope ratios within and between Setaria italica plants grown in a controlled environment chamber. Using homogenised samples, we compare carbon isotope ratios of leaves and grains, and nitrogen isotope ratios of grains, from 29 accessions of Setaria italica. RESULTS: We find significant isotopic variability within single leaves and panicles, and between leaves and panicles within the same plant, which must be considered when undertaking plant isotope studies. We find that the leaves and grains from the different accessions have a ca 2‰ range in δ(13) C values, while the nitrogen isotope values in the grains have a ca 6‰ range. We also find an average offset of 0.9‰ between leaves and grains in their δ(13) C values. CONCLUSIONS: The variation found is large enough to have archaeological implications and within- and between-plant isotope variability should be considered in isotope studies. The range in δ(15) N values is particularly significant as it is larger than the typical values quoted for a trophic level enrichment, and as such may lead to erroneous interpretations of the amount of animal protein in human or animal diets. It is therefore necessary to account for the variability in plant stable isotope values during palaeodietary reconstructions. Copyright © 2016 John Wiley & Sons, Ltd.

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.

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.

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.

Waxy phenotype evolution in the allotetraploid cereal broomcorn millet: mutations at the GBSSI locus in their functional and phylogenetic context.

Waxy mutants, in which endosperm starch contains ~100% amylopectin rather than the wild-type composition of ~70% amylopectin and ~30% amylose, occur in many domesticated cereals. The cultivation of waxy varieties is concentrated in east Asia, where there is a culinary preference for glutinous-textured foods that may have developed from ancient food processing traditions. The waxy phenotype results from mutations in the GBSSI gene, which catalyzes amylose synthesis. Broomcorn or proso millet (Panicum miliaceum L.) is one of the world's oldest cultivated cereals, which spread across Eurasia early in prehistory. Recent phylogeographic analysis has shown strong genetic structuring that likely reflects ancient expansion patterns. Broomcorn millet is highly unusual in being an allotetraploid cereal with fully waxy varieties. Previous work characterized two homeologous GBSSI loci, with multiple alleles at each, but could not determine whether both loci contributed to GBSSI function. We first tested the relative contribution of the two GBSSI loci to amylose synthesis and second tested the association between GBSSI alleles and phylogeographic structure inferred from simple sequence repeats (SSRs). We evaluated the phenotype of all known GBSSI genotypes in broomcorn millet by assaying starch composition and protein function. The results showed that the GBSSI-S locus is the major locus controlling endosperm amylose content, and the GBSSI-L locus has strongly reduced synthesis capacity. We genotyped 178 individuals from landraces from across Eurasia for the 2 GBSSI and 16 SSR loci and analyzed phylogeographic structuring and the geographic and phylogenetic distribution of GBSSI alleles. We found that GBSSI alleles have distinct spatial distributions and strong associations with particular genetic clusters defined by SSRs. The combination of alleles that results in a partially waxy phenotype does not exist in landrace populations. Our data suggest that broomcorn millet is a system in the process of becoming diploidized for the GBSSI locus responsible for grain amylose. Mutant alleles show some exchange between genetic groups, which was favored by selection for the waxy phenotype in particular regions. Partially waxy phenotypes were probably selected against-this unexpected finding shows that better understanding is needed of the human biology of this phenomenon that distinguishes cereal use in eastern and western cultures.

Reticulate evolution in Panicum (Poaceae): the origin of tetraploid broomcorn millet, P. miliaceum.

Panicum miliaceum (broomcorn millet) is a tetraploid cereal, which was among the first domesticated crops, but is now a minor crop despite its high water use efficiency. The ancestors of this species have not been determined; we aimed to identify likely candidates within the genus, where phylogenies are poorly resolved. Nuclear and chloroplast DNA sequences from P. miliaceum and a range of diploid and tetraploid relatives were used to develop phylogenies of the diploid and tetraploid species. Chromosomal in situ hybridization with genomic DNA as a probe was used to characterize the genomes in the tetraploid P. miliaceum and a tetraploid accession of P. repens. In situ hybridization showed that half the chromosomes of P. miliaceum hybridized more strongly with labelled genomic DNA from P. capillare, and half with labelled DNA from P. repens. Genomic DNA probes differentiated two sets of 18 chromosomes in the tetraploid P. repens. Our phylogenetic data support the allotetraploid origin of P. miliaceum, with the maternal ancestor being P. capillare (or a close relative) and the other genome being shared with P. repens. Our P. repens accession was also an allotetraploid with two dissimilar but closely related genomes, the maternal genome being similar to P. sumatrense. Further collection of Panicum species, particularly from the Old World, is required. It is important to identify why the water-efficient P. miliaceum is now of minimal importance in agriculture, and it may be valuable to exploit the diversity in this species and its ancestors.

Wheat in the Mediterranean revisited--tetraploid wheat landraces assessed with elite bread wheat Single Nucleotide Polymorphism markers.

BACKGROUND: Single Nucleotide Polymorphism (SNP) panels recently developed for the assessment of genetic diversity in wheat are primarily based on elite varieties, mostly those of bread wheat. The usefulness of such SNP panels for studying wheat evolution and domestication has not yet been fully explored and ascertainment bias issues can potentially affect their applicability when studying landraces and tetraploid ancestors of bread wheat. We here evaluate whether population structure and evolutionary history can be assessed in tetraploid landrace wheats using SNP markers previously developed for the analysis of elite cultivars of hexaploid wheat. RESULTS: We genotyped more than 100 tetraploid wheat landraces and wild emmer wheat accessions, some of which had previously been screened with SSR markers, for an existing SNP panel and obtained publically available genotypes for the same SNPs for hexaploid wheat varieties and landraces. Results showed that quantification of genetic diversity can be affected by ascertainment bias but that the effects of ascertainment bias can at least partly be alleviated by merging SNPs to haplotypes. Analyses of population structure and genetic differentiation show strong subdivision between the tetraploid wheat subspecies, except for durum and rivet that are not separable. A more detailed population structure of durum landraces could be obtained than with SSR markers. The results also suggest an emmer, rather than durum, ancestry of bread wheat and with gene flow from wild emmer. CONCLUSIONS: SNP markers developed for elite cultivars show great potential for inferring population structure and can address evolutionary questions in landrace wheat. Issues of marker genome specificity and mapping need, however, to be addressed. Ascertainment bias does not seem to interfere with the ability of a SNP marker system developed for elite bread wheat accessions to detect population structure in other types of wheat.

Volumetric determination of apolipoprotein stoichiometry of circulating HDL subspecies.

Although HDL is inversely correlated with coronary heart disease, elevated HDL-cholesterol is not always protective. Additionally, HDL has biological functions that transcend any antiatherogenic role: shotgun proteomics show that HDL particles contain 84 proteins (latest count), many correlating with antioxidant and anti-inflammatory properties of HDL. ApoA-I has been suggested to serve as a platform for the assembly of these protein components on HDL with specific functions - the HDL proteome. However, the stoichiometry of apoA-I in HDL subspecies is poorly understood. Here we use a combination of immunoaffinity chromatography data and volumetric analysis to evaluate the size and stoichiometry of LpA-I and LpA-I,A-II particles. We conclude that there are three major LpA-I subspecies: two major particles, HDL[4] in the HDL3 size range (d = 85.0 ± 1.2 Å) and HDL[7] in the HDL2 size range (d = 108.5 ± 3.8 Å) with apoA-I stoichiometries of 3 and 4, respectively, and a small minor particle, HDL[1] (d = 73.8 ± 2.1Å) with an apoA-I stoichiometry of 2. Additionally, we conclude that the molar ratio of apolipoprotein to surface lipid is significantly higher in circulating HDL subspecies than in reconstituted spherical HDL particles, presumably reflecting a lack of phospholipid transfer protein in reconstitution protocols.

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.

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.

Rotational and hinge dynamics of discoidal high density lipoproteins probed by interchain disulfide bond formation.

To develop a detailed double belt model for discoidal HDL, we previously scored inter-helical salt bridges between all possible registries of two stacked antiparallel amphipathic helical rings of apolipoprotein (apo) A-I. The top score was the antiparallel apposition of helix 5 with 5 followed closely by appositions of helix 5 with 4 and helix 5 with 6. The rationale for the current study is that, for each of the optimal scores, a pair of identical residues can be identified in juxtaposition directly on the contact edge between the two antiparallel helical belts of apoA-I. Further, these residues are always in the '9th position' in one of the eighteen 11-mer repeats that make up the lipid-associating domain of apoA-I. To illustrate our terminology, 129j (LL5/5) refers to the juxtaposition of the Cα atoms of G129 (in a '9th position') in the pairwise helix 5 domains. We reasoned that if identical residues in the double belt juxtapositions were mutated to a cysteine and kept under reducing conditions during disc formation, we would have a precise method for determining registration in discoidal HDL by formation of a disulfide-linked apoA-I homodimer. Using this approach, we conclude that 129j (LL5/5) is the major rotamer orientation for double belt HDL and propose that the small ubiquitous gap between the pairwise helix 5 portions of the double belt in larger HDL discoidal particles is significantly dynamic to hinge off the disc edge under certain conditions, e.g., in smaller particles or perhaps following binding of the enzyme LCAT. This article is part of a Special Issue entitled Advances in High Density Lipoprotein Formation and Metabolism: A Tribute to John F. Oram (1945-2010).

Genome sequencing reveals the genetic architecture of heterostyly and domestication history of common buckwheat.

Common buckwheat, Fagopyrum esculentum, is an orphan crop domesticated in southwest China that exhibits heterostylous self-incompatibility. Here we present chromosome-scale assemblies of a self-compatible F. esculentum accession and a self-compatible wild relative, Fagopyrum homotropicum, together with the resequencing of 104 wild and cultivated F. esculentum accessions. Using these genomic data, we report the roles of transposable elements and whole-genome duplications in the evolution of Fagopyrum. In addition, we show that (1) the breakdown of heterostyly occurs through the disruption of a hemizygous gene jointly regulating the style length and female compatibility and (2) southeast Tibet was involved in common buckwheat domestication. Moreover, we obtained mutants conferring the waxy phenotype for the first time in buckwheat. These findings demonstrate the utility of our F. esculentum assembly as a reference genome and promise to accelerate buckwheat research and breeding.

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.

The earliest evidence of millet as a staple crop: New light on neolithic foodways in North China.

There is a growing body of archaeobotanical evidence for the harvesting of millet in Eurasia prior to 5,000 cal. BC. Yet direct evidence for the extent of millet consumption in this time period is rare. This contradiction may be due to millet crops making only a minor contribution to the diet before 5,000 BC. In this article, drawing from recent excavations in North China, we present evidence for millet crops making a substantial contribution to human and animal diets in periods, which correspond chronologically with the time depth of the archaeobotanical record. We infer that in eastern Inner Mongolia, human adoption of millets, which may or may be not related to substantial agriculture, happened at the Early Neolithic, with direct dates between 5,800 and 5,300 cal. BC.

ABCA1 is an extracellular phospholipid translocase.

Production of high density lipoprotein (HDL) requires ATP-binding cassette transporter A1 (ABCA1) to drive phospholipid (PL) from the plasma membrane into extracellular apolipoprotein A-I. Here, we use simulations to show that domains of ABCA1 within the plasma membrane remove PL from the membrane's outer leaflet. In our simulations, after the lipid diffuses into the interior of ABCA1's outward-open cavity, PL extracted by the gateway passes through a ring-shaped domain, the annulus orifice, which forms the base of an elongated hydrophobic tunnel in the transporter's extracellular domain. Engineered mutations in the gateway and annulus strongly inhibit lipid export by ABCA1 without affecting cell-surface expression levels. Our finding that ABCA1 extracts lipid from the outer face of the plasma membrane and forces it through its gateway and annulus into an elongated hydrophobic tunnel contrasts with the alternating access model, which proposes that ABCA1 flops PL substrate from the inner leaflet to the outer leaflet of the membrane. Consistent with our model, ABCA1 lacks the charged amino acid residues in the transmembrane domain found in the floppase members of the ABC transporter family.

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.

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