Search

Anion exchanger-1 (AE1) is the main erythroid Cl - /HCO 3 - transporter that supports CO 2 transport. Glycophorin A (GPA), a component of the AE1 complexes, facilitates AE1 expression and anion transport, but Glycophorin B (GPB) does not. Here, we dissected the structural components of GPA/GPB involved in glycophorin-AE1 trafficking by comparing them with three GPB variants-GPBhead (lacking the transmembrane domain [TMD]), GPBtail (mainly the TMD), and GP.Mur (glycophorin B-A-B hybrid). GPB-derived GP.Mur bears an O-glycopeptide that encompasses the R18 epitope, which is present in GPA but not GPB. By flow cytometry, AE1 expression in the control erythrocytes increased with the GPA-R18 expression; GYP.Mur + / + erythrocytes bearing both GP.Mur and GPA expressed more R18 epitopes and more AE1 proteins. In contrast, heterologously expressed GPBtail and GPB were predominantly localized in the Golgi apparatus of HEK-293 cells, whereas GBhead was diffuse throughout the cytosol, suggesting that glycophorin transmembrane encoded an ER/Golgi retention signal. AE1 coexpression could reduce the ER/Golgi retention of GPB, but not of GPBtail or GPBhead. Thus, there are forward-trafficking and transmembrane-driven ER/Golgi retention signals encoded in the glycophorin sequences. How the balance between these opposite trafficking signals could affect glycophorin sorting into AE1 complexes and influence erythroid anion transport remains to be explored.

In the erythrocyte membrane, the interactions between glycophorin A (GPA) and Band 3 are associated strongly with the biological function of the membrane and several blood disorders. In this work, using coarse-grained molecular-dynamics simulations, we systematically investigate the effects of cholesterol and phosphatidylinositol-4,5-bisphosphate (PIP2) on the interactions of GPA with Band 3 in the model erythrocyte membranes. We examine the dynamics of the interactions of GPA with Band 3 in different lipid bilayers on the microsecond time scale and calculate the binding free energy between GPA and Band 3. The results indicate that cholesterols thermodynamically favor the binding of GPA to Band 3 by increasing the thickness of the lipid bilayer and by producing an effective attraction between the proteins due to the depletion effect. Cholesterols also slow the kinetics of the binding of GPA to Band 3 by reducing the lateral mobility of the lipids and proteins and may influence the binding sites between the proteins. The anionic PIP2 lipids prefer binding to the surface of the proteins through electrostatic attraction between the PIP2 headgroup and the positively charged residues on the protein surface. Ions in the solvent facilitate PIP2 aggregation, which promotes the binding of GPA to Band 3., (Copyright © 2022 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Low plasma levels of Proprotein Convertase Subtilisin/Kexin 9 (PCSK9) are associated with decreased low-density lipoprotein (LDL) cholesterol and a reduced risk of cardiovascular disease. PCSK9 binds to the epidermal growth factor-like repeat A (EGFA) domain of LDL receptors (LDLR), very low-density lipoprotein receptors (VLDLR), apolipoprotein E receptor 2 (ApoER2), and lipoprotein receptor-related protein 1 (LRP1) and accelerates their degradation, thus acting as a key regulator of lipid metabolism. Antibody and RNAi-based PCSK9 inhibitor treatments lower cholesterol and prevent cardiovascular incidents in patients, but their high-cost hampers market penetration. We sought to develop a safe, long-term and one-time solution to treat hyperlipidemia. We created a cDNA encoding a chimeric protein in which the extracellular N- terminus of red blood cells (RBCs) specific glycophorin A was fused to the LDLR EGFA domain and introduced this gene into mouse bone marrow hematopoietic stem and progenitor cells (HSPCs). Following transplantation into irradiated mice, the animals produced RBCs with the EGFA domain (EGFA-GPA RBCs) displayed on their surface. These animals showed significantly reduced plasma PCSK9 (66.5% decrease) and reduced LDL levels (40% decrease) for as long as 12 months post-transplantation. Furthermore, the EGFA- GPA mice remained lean for life and maintained normal body weight under a high-fat diet. Hematopoietic stem cell gene therapy can generate red blood cells expressing an EGFA-glycophorin A chimeric protein as a practical and long-term strategy for treating chronic hyperlipidemia and obesity., Competing Interests: H.F.L. and H.L.P. served as scientific advisors and have/had equity in Rubius, a biotechnology company that seeks to exploit engineered red blood cells as therapeutics but does not provide financial support for the technology described in this paper. This does not alter our adherence to PLOS ONE policies on sharing data and materials.

Transmembrane helix association is a fundamental step in the folding of helical membrane proteins. The prototypical example of this association is formation of the glycophorin dimer. While its structure and stability have been well-characterized experimentally, the detailed assembly mechanism is harder to obtain. Here, we use all-atom simulations within phospholipid membrane to study glycophorin association. We find that initial association results in the formation of a non-native intermediate, separated by a significant free energy barrier from the dimer with a native binding interface. We have used transition-path sampling to determine the association mechanism. We find that the mechanism of the initial bimolecular association to form the intermediate state can be mediated by many possible contacts, but seems to be particularly favoured by formation of non-native contacts between the C-termini of the two helices. On the other hand, the contacts which are key to determining progression from the intermediate to the native state are those which define the native binding interface, reminiscent of the role played by native contacts in determining folding of globular proteins. As a check on the simulations, we have computed association and dissociation rates from the transition-path sampling. We obtain results in reasonable accord with available experimental data, after correcting for differences in native state stability. Our results yield an atomistic description of the mechanism for a simple prototype of helical membrane protein folding., Competing Interests: The authors have declared that no competing interests exist.

The biconcave disk shape of the mammalian red blood cell (RBC) is unique to the RBC and is vital for its circulatory function. Due to the absence of a transcellular cytoskeleton, RBC shape is determined by the membrane skeleton, a network of actin filaments cross-linked by spectrin and attached to membrane proteins. While the physical properties of a uniformly distributed actin network interacting with the lipid bilayer membrane have been assumed to control RBC shape, recent experiments reveal that RBC biconcave shape also depends on the contractile activity of nonmuscle myosin IIA (NMIIA) motor proteins. Here, we use the classical Helfrich-Canham model for the RBC membrane to test the role of heterogeneous force distributions along the membrane and mimic the contractile activity of sparsely distributed NMIIA filaments. By incorporating this additional contribution to the Helfrich-Canham energy, we find that the RBC biconcave shape depends on the ratio of forces per unit volume in the dimple and rim regions of the RBC. Experimental measurements of NMIIA densities at the dimple and rim validate our prediction that (a) membrane forces must be non-uniform along the RBC membrane and (b) the force density must be larger in the dimple than the rim to produce the observed membrane curvatures. Furthermore, we predict that RBC membrane tension and the orientation of the applied forces play important roles in regulating this force-shape landscape. Our findings of heterogeneous force distributions on the plasma membrane for RBC shape maintenance may also have implications for shape maintenance in different cell types., Competing Interests: The authors have declared that no competing interests exist.

New methods for quantifying extracellular vesicles (EVs) in complex biofluids are critically needed. We report the development of a new technology combining size exclusion chromatography (SEC), a commonly used EV purification technique, with fluorescence detection of specifically labelled EVs. The resulting platform, Flu-SEC, demonstrates a linear response to concentration of specific EVs and could form the basis of a system with phenotyping capability. Flu-SEC was validated using red blood cell derived EVs (REVs), which provide an ideal EV model with monodisperse size distribution and high EV concentration. Microfluidic Resistive Pulse Sensing (MRPS) was used to accurately determine the size distribution and concentration of REVs. Anti-CD235a antibody, specific to glycophorin A, and the more general wheat germ agglutinin (WGA), were selected to label REVs. The results show the quantitative power of Flu-SEC: a highly linear fluorescence response over a wide range of concentrations. Moreover, the Flu-SEC technique reports the ratio of EV-bound and free-antibody molecules, an important metric for determining optimal labelling conditions for other applications. Flu-SEC represents an orthogonal tool to single-particle fluorescent methods such as flow cytometry and fluorescent NTA, for the quantification and phenotyping of EVs.

The Band 3 (AE1, SLC4A1) membrane protein is found in red blood cells and in kidney where it functions as an electro-neutral chloride/bicarbonate exchanger. In this study, we have used molecular dynamics simulations to provide the first realistic model of the dimeric membrane domain of human Band 3 in an asymmetric lipid bilayer containing a full complement of phospholipids, including phosphatidylinositol 4,5-bisphosphate (PIP2) and cholesterol, and its partner membrane protein Glycophorin A (GPA). The simulations show that the annular layer in the inner leaflet surrounding Band 3 was enriched in phosphatidylserine and PIP2 molecules. Cholesterol was also enriched around Band 3 but also at the dimer interface. The interaction of these lipids with specific sites on Band 3 may play a role in the folding and function of this anion transport membrane protein. GPA associates with Band 3 to form the Wright (Wr) blood group antigen, an interaction that involves an ionic bond between Glu658 in Band 3 and Arg61 in GPA. We were able to recreate this complex by performing simulations to allow the dimeric transmembrane portion of GPA to interact with Band 3 in a model membrane. Large-scale simulations showed that the GPA dimer can bridge Band 3 dimers resulting in the dynamic formation of long strands of alternating Band 3 and GPA dimers., Competing Interests: The authors have declared that no competing interests exist.

Lectin-glycan interactions play important roles in many biological systems, but the nature of glycoprotein counter-receptors expressed on cell membranes is often poorly understood. To help overcome this problem, we developed a method based on proximity labeling technology. Using a peroxidase-coupled lectin, addition of H2O2 and tyramide-biotin substrates leads to generation of short-range biotin radicals that biotinylate proteins in the immediate vicinity of the bound lectin, which can subsequently be identified. As a proof-of-principle, sialoadhesin-horseradish peroxidase-human IgG1 Fc recombinant protein constructs were precomplexed with anti-Fc antibodies, bound to human erythrocytes and reacted with H2O2 and tyramide-SS-biotin. The erythrocyte membrane protein with strongest biotinylation was identified as glycophorin A, in agreement with early studies using lectin overlay and reglycosylation approaches. As a further test of the method, the plant lectin MAL II was conjugated with horseradish peroxidase and used in proximity labeling of human erythrocytes. Glycophorin A was again selectively labeled, which is consistent with previous reports that MAL II has high affinity for glycophorin. This method could be applied to other lectins to identify their membrane counter-receptors., (© The Author 2017. Published by Oxford University Press.)

The assembly of integral membrane protein complexes is frequently supported by transmembrane domain (TMD) interactions. Here, we present the BLaTM assay that measures homotypic as well as heterotypic TMD-TMD interactions in a bacterial membrane. The system is based on complementation of β-lactamase fragments genetically fused to interacting TMDs, which confers ampicillin resistance to expressing cells. We validated BLaTM by showing that the assay faithfully reports known sequence-specific interactions of both types. In a practical application, we used BLaTM to screen a focussed combinatorial library for heterotypic interactions driven by electrostatic forces. The results reveal novel patterns of ionizable amino acids within the isolated TMD pairs. Those patterns indicate that formation of heterotypic TMD pairs is most efficiently supported by closely spaced ionizable residues of opposite charge. In addition, TMD heteromerization can apparently be driven by hydrogen bonding between basic or between acidic residues.

Analysis of intact glycopeptides by mass spectrometry is essential to determining the microheterogeneity of protein glycosylation. Higher-energy collisional dissociation (HCD) fragmentation of glycopeptides generates mono- or disaccharide ions called oxonium ions that carry information about the structure of the fragmented glycans. Here, we investigated the link between glycan structures and the intensity of oxonium ions in the spectra of glycopeptides and utilized this information to improve the identification of glycopeptides in biological samples. Tandem spectra of glycopeptides from fetuin, glycophorin A, ovalbumin and gp120 tryptic digests were used to build a spectral database of N- and O-linked glycopeptides. Logistic regression was applied to this database to develop model to distinguish between the spectra of N- and O-linked glycopeptides. Remarkably, the developed model was found to reliably distinguish between the N- and O-linked glycopeptides using the spectral features of the oxonium ions using verification spectral set. Finally, the performance of the developed predictive model was evaluated in HILIC enriched glycopeptides extracted from human serum. The results showed that pre-classification of tandem spectra based on their glycosylation type improved the identification of N-linked glycopeptides. The developed model facilitates interpretation of tandem mass spectrometry data for assignment of glycopeptides.

Insertion of helix-forming segments into the membrane and their association determines the structure, function, and expression levels of all plasma membrane proteins. However, systematic and reliable quantification of membrane-protein energetics has been challenging. We developed a deep mutational scanning method to monitor the effects of hundreds of point mutations on helix insertion and self-association within the bacterial inner membrane. The assay quantifies insertion energetics for all natural amino acids at 27 positions across the membrane, revealing that the hydrophobicity of biological membranes is significantly higher than appreciated. We further quantitate the contributions to membrane-protein insertion from positively charged residues at the cytoplasm-membrane interface and reveal large and unanticipated differences among these residues. Finally, we derive comprehensive mutational landscapes in the membrane domains of Glycophorin A and the ErbB2 oncogene, and find that insertion and self-association are strongly coupled in receptor homodimers.

Membranes are central for cells as borders to the environment or intracellular organelle definition. They are composed of and harbor different molecules like various lipid species and sterols, and they are generally crowded with proteins. The membrane system is very dynamic and components show lateral, rotational and translational diffusion. The consequence of the latter is that phase separation can occur in membranes in vivo and in vitro. It was documented that molecular dynamics simulations of an idealized plasma membrane model result in formation of membrane areas where either saturated lipids and cholesterol (liquid-ordered character, Lo) or unsaturated lipids (liquid-disordered character, Ld) were enriched. Furthermore, current discussions favor the idea that proteins are sorted into the liquid-disordered phase of model membranes, but experimental support for the behavior of isolated proteins in native membranes is sparse. To gain insight into the protein behavior we built a model of the red blood cell membrane with integrated glycophorin A dimer. The sorting and the dynamics of the dimer were subsequently explored by coarse-grained molecular dynamics simulations. In addition, we inspected the impact of lipid head groups and the presence of cholesterol within the membrane on the dynamics of the dimer within the membrane. We observed that cholesterol is important for the formation of membrane areas with Lo and Ld character. Moreover, it is an important factor for the reproduction of the dynamic behavior of the protein found in its native environment. The protein dimer was exclusively sorted into the domain of Ld character in the model red blood cell plasma membrane. Therefore, we present structural information on the glycophorin A dimer distribution in the plasma membrane in the absence of other factors like e.g. lipid anchors in a coarse grain resolution.

Glycophorins C and D are highly glycosylated integral sialoglycoproteins of human red blood cell membranes carrying the Gerbich blood group antigens. The O- and N-glycosidic chains of the major erythrocyte glycoprotein (Lisowska E. 2001, Antigenic properties of human glycophorins - an update. Adv Exp Med Biol, 491:155-169; Tomita M and Marchesi VT. 1975, Amino-acid sequence and oligosaccharide attachment sites of human erythrocyte glycophorin. Proc Natl Acad Sci USA, 72:2964-2968.) are well characterized but the structure of GPC N-glycans has remained unknown. This problem became important since it was reported that GPC N-glycans play an essential role in the interaction with Plasmodium falciparum EBA-140 merozoite ligand. The elucidation of these structures seems essential for full characterization of the GPC binding site for the EBA-140 ligand. We have employed detailed structural analysis using sequential mass spectrometry to show that many GPC N-glycans contain H2 antigen structures and several contain polylactosamine structures capped with fucose. The results obtained indicate structural heterogeneity of the GPC N-glycans and show the existence of structural elements not found in glycophorin A N-glycans. Our results also open a possibility of new interpretation of the data concerning the binding of P. falciparum EBA-140 ligand to GPC. We hypothesize that preferable terminal fucosylation of N-glycosidic chains containing repeating lactosamine units of the GPC Gerbich variant could be an explanation for why the EBA-140 ligand does not react with GPC Gerbich and an indication that the EBA-140 interaction with GPC is distinctly dependent on the GPC N-glycan structure., (© The Author 2014. Published by Oxford University Press. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Plasmodium falciparum invasion of human red blood cells (RBCs) is an intricate process requiring a number of distinct ligand-receptor interactions at the merozoite-erythrocyte interface. Merozoite surface protein 1 (MSP1), a highly abundant ligand coating the merozoite surface in all species of malaria parasites, is essential for RBC invasion and considered a leading candidate for inclusion in a multiple-subunit vaccine against malaria. Our previous studies identified an interaction between the carboxyl-terminus of MSP1 and RBC band 3. Here, by employing phage display technology, we report a novel interaction between the amino-terminus of MSP1 and RBC glycophorin A (GPA). Mapping of the binding domains established a direct interaction between malaria MSP1 and human GPA within a region of MSP1 known to potently inhibit P falciparum invasion of human RBCs. Furthermore, a genetically modified mouse model lacking the GPA- band 3 complex in RBCs is completely resistant to malaria infection in vivo. These findings suggest an essential role of the MSP1-GPA-band 3 complex during the initial adhesion phase of malaria parasite invasion of RBCs., (© 2015 by The American Society of Hematology.)

As traditional detergents might destabilize or even denature membrane proteins, amphiphilic polymers have moved into the focus of membrane-protein research in recent years. Thus far, Amphipols are the best studied amphiphilic copolymers, having a hydrophilic backbone with short hydrophobic chains. However, since stabilizing as well as destabilizing effects of the Amphipol belt on the structure of membrane proteins have been described, we systematically analyze the impact of the most commonly used Amphipol A8-35 on the structure and stability of a well-defined transmembrane protein model, the glycophorin A transmembrane helix dimer. Amphipols are not able to directly extract proteins from their native membranes, and detergents are typically replaced by Amphipols only after protein extraction from membranes. As Amphipols form mixed micelles with detergents, a better understanding of Amphipol-detergent interactions is required. Therefore, we analyze the interaction of A8-35 with the anionic detergent sodium dodecyl sulfate and describe the impact of the mixed-micelle-like system on the stability of a transmembrane helix dimer. As A8-35 may highly stabilize and thereby rigidify a transmembrane protein structure, modest destabilization by controlled addition of detergents and formation of mixed micellar systems might be helpful to preserve the function of a membrane protein in Amphipol environments.

The hydrophobic thickness of membranes, which is manly defined by fatty acids, influences the packing of transmembrane domains of proteins and thus can modulate the activity of these proteins. We analyzed the dynamics of the dimerization of Glycophorin A (GpA) by molecular dynamics simulations to describe the fatty acid dependence of the transmembrane region assembly. GpA represents a well-established model for dimerization of single transmembrane helices containing a GxxxG motif in vitro and in silico. We performed simulations of the dynamics of the NMR-derived dimer as well as self-assembly simulations of monomers in membranes composed of different fatty acid chains and monitored the formed interfaces and their transitions. The observed dimeric interfaces, which also include the one known from NMR, are highly dynamic and converted into each other. The frequency of interface formation and the preferred transitions between interfaces similar to the interface observed by NMR analysis strongly depend on the fatty acid used to build the membrane. Molecular dynamic simulations after adaptation of the helix topology parameters to better represent NMR derived structures of single transmembrane helices yielded an enhanced occurrence of the interface determined by NMR in molecular dynamics simulations. Taken together we give insights into the influence of fatty acids and helix conformation on the dynamics of the transmembrane domain of GpA.

The malaria parasite, Plasmodium falciparum, and related parasites use a variety of proteins with Duffy-Binding Like (DBL) domains to bind glycoproteins on the surface of host cells. Among these proteins, the 175 kDa erythrocyte binding antigen, EBA-175, specifically binds to glycophorin A on the surface of human erythrocytes during the process of merozoite invasion. The domain responsible for glycophorin A binding was identified as region II (RII) which contains two DBL domains, F1 and F2. The crystal structure of this region revealed a dimer that is presumed to represent the glycophorin A binding conformation as sialic acid binding sites and large cavities are observed at the dimer interface. The dimer interface is largely composed of two loops from within each monomer, identified as the F1 and F2 β-fingers that contact depressions in the opposing monomers in a similar manner. Previous studies have identified a panel of five monoclonal antibodies (mAbs) termed R215 to R218 and R256 that bind to RII and inhibit invasion of erythrocytes to varying extents. In this study, we predict the F2 β-finger region as the conformational epitope for mAbs, R215, R217, and R256, and confirm binding for the most effective blocking mAb R217 and R215 to a synthetic peptide mimic of the F2 β-finger. Localization of the epitope to the dimerization and glycan binding sites of EBA-175 RII and site-directed mutagenesis within the predicted epitope are consistent with R215 and R217 blocking erythrocyte invasion by Plasmodium falciparum by preventing formation of the EBA-175- glycophorin A complex.

Disrupting erythrocyte invasion by Plasmodium falciparum is an attractive approach to combat malaria. P. falciparum EBA-175 (PfEBA-175) engages the host receptor Glycophorin A (GpA) during invasion and is a leading vaccine candidate. Antibodies that recognize PfEBA-175 can prevent parasite growth, although not all antibodies are inhibitory. Here, using x-ray crystallography, small-angle x-ray scattering and functional studies, we report the structural basis and mechanism for inhibition by two PfEBA-175 antibodies. Structures of each antibody in complex with the PfEBA-175 receptor binding domain reveal that the most potent inhibitory antibody, R217, engages critical GpA binding residues and the proposed dimer interface of PfEBA-175. A second weakly inhibitory antibody, R218, binds to an asparagine-rich surface loop. We show that the epitopes identified by structural studies are critical for antibody binding. Together, the structural and mapping studies reveal distinct mechanisms of action, with R217 directly preventing receptor binding while R218 allows for receptor binding. Using a direct receptor binding assay we show R217 directly blocks GpA engagement while R218 does not. Our studies elaborate on the complex interaction between PfEBA-175 and GpA and highlight new approaches to targeting the molecular mechanism of P. falciparum invasion of erythrocytes. The results suggest studies aiming to improve the efficacy of blood-stage vaccines, either by selecting single or combining multiple parasite antigens, should assess the antibody response to defined inhibitory epitopes as well as the response to the whole protein antigen. Finally, this work demonstrates the importance of identifying inhibitory-epitopes and avoiding decoy-epitopes in antibody-based therapies, vaccines and diagnostics.

Detergents might affect membrane protein structures by promoting intramolecular interactions that are different from those found in native membrane bilayers, and fine-tuning detergent properties can be crucial for obtaining structural information of intact and functional transmembrane proteins. To systematically investigate the influence of the detergent concentration and acyl-chain length on the stability of a transmembrane protein structure, the stability of the human glycophorin A transmembrane helix dimer has been analyzed in lyso-phosphatidylcholine micelles of different acyl-chain length. While our results indicate that the transmembrane protein is destabilized in detergents with increasing chain-length, the diameter of the hydrophobic micelle core was found to be less crucial. Thus, hydrophobic mismatch appears to be less important in detergent micelles than in lipid bilayers and individual detergent molecules appear to be able to stretch within a micelle to match the hydrophobic thickness of the peptide. However, the stability of the GpA TM helix dimer linearly depends on the aggregation number of the lyso-PC detergents, indicating that not only is the chemistry of the detergent headgroup and acyl-chain region central for classifying a detergent as harsh or mild, but the detergent aggregation number might also be important., (Copyright © 2012 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Interaction of protein 4.1 (4.1R) with the transmembrane protein glycophorin C (GPC) regulates the functions of erythrocyte membrane. Fluorescence correlation spectroscopy (FCS) was used to define the interaction of EGFP-4.1R with DsRed-GPC on transport vesicles (TVs) by measuring their fluctuation in living cells. DsRed-GPC expressed in HeLa cells was delivered to the plasma membrane through slow vesicle transport. EGFP-4.1R, which freely diffused in the cytosol when expressed alone, diffused slowly when co-expressed with DsRed-GPC, indicating association of EGFP-4.1R with TVs. Fluorescence cross-correlation spectroscopy (FCCS) showed direct interaction of EGFP-4.1R with DsRed-GPC on TVs. The present study demonstrates that 4.1R binds to GPC on TVs in living cells., (Copyright © 2012 Federation of European Biochemical Societies. Published by Elsevier B.V. All rights reserved.)

The vast majority of membrane proteins are anchored to biological membranes through hydrophobic α-helices. Sequence analysis of high-resolution membrane protein structures show that ionizable amino acid residues are present in transmembrane (TM) helices, often with a functional and/or structural role. Here, using as scaffold the hydrophobic TM domain of the model membrane protein glycophorin A (GpA), we address the consequences of replacing specific residues by ionizable amino acids on TM helix insertion and packing, both in detergent micelles and in biological membranes. Our findings demonstrate that ionizable residues are stably inserted in hydrophobic environments, and tolerated in the dimerization process when oriented toward the lipid face, emphasizing the complexity of protein-lipid interactions in biological membranes.

Malaria has been a very strong selection pressure in recent human evolution, particularly in Africa. Of the one million deaths per year due to malaria, more than 90% are in sub-Saharan Africa, a region with high levels of genetic variation and population substructure. However, there have been few studies of nucleotide variation at genetic loci that are relevant to malaria susceptibility across geographically and genetically diverse ethnic groups in Africa. Invasion of erythrocytes by Plasmodium falciparum parasites is central to the pathology of malaria. Glycophorin A (GYPA) and B (GYPB), which determine MN and Ss blood types, are two major receptors that are expressed on erythrocyte surfaces and interact with parasite ligands. We analyzed nucleotide diversity of the glycophorin gene family in 15 African populations with different levels of malaria exposure. High levels of nucleotide diversity and gene conversion were found at these genes. We observed divergent patterns of genetic variation between these duplicated genes and between different extracellular domains of GYPA. Specifically, we identified fixed adaptive changes at exons 3-4 of GYPA. By contrast, we observed an allele frequency spectrum skewed toward a significant excess of intermediate-frequency alleles at GYPA exon 2 in many populations; the degree of spectrum distortion is correlated with malaria exposure, possibly because of the joint effects of gene conversion and balancing selection. We also identified a haplotype causing three amino acid changes in the extracellular domain of glycophorin B. This haplotype might have evolved adaptively in five populations with high exposure to malaria., (Copyright © 2011 The American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

In criminal investigations, there are three stages involved when studying bloodstains: search and orientation, confirmation, and individualization. Confirmatory tests have two aims: to show that the stain contains a human biological fluid and to confirm the type of biological fluid. The need to determine the nature of the evidence is reflected in the latest bibliography, where the possibility of employing mRNA and miRNA markers for this purpose is proposed. While these new proposals are being investigated, the kits for determining human hemoglobin currently provide a simple solution for resolving this issue. With these kits, the possibility of obtaining false positives and false negatives is well known. However, recently, a new problem has been detected. This involves the interference caused by new cleaning products that contain sodium percarbonate (or active oxygen) when determining human hemoglobin. With the aim to resolve this problem, this work studied the ability of the human glycophorin A test to determine human blood in samples that have been treated with active oxygen. Our results show that the human glycophorin A test has a greater resistance to the destructive effect of the new detergents containing active oxygen; consequently, it provides an alternative to be taken into consideration in the confirmatory diagnoses of bloodstains.

Deuterium NMR spectroscopy was used to study how the positioning of a dimerization motif within a transbilayer polypeptide influences its orientation and dynamics in bilayers. Three polypeptide variants comprising glycophorin A transmembrane (GpATM) dimerization motifs incorporated into lysine-terminated poly-leucine-alanine helices were mixed into 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine multilamellar vesicles. The variants differed in orientation of the motif segment around the helix axis with respect to the peptide ends. Polypeptides were labeled with methyl-deuterated alanines at positions that were identically situated relative to the peptide ends (Ala-20 and Ala-22) and at two positions within the motif. An analysis of quadrupole splittings revealed similar tilts and orientations of the peptide ends for all three variants, suggesting that average orientations were dominated by interactions at the bilayer surface. For one variant, however, fast orientational fluctuations about the helix axis were significantly smaller. This may indicate some perturbation of peptide dynamics and conformation by interactions that are sensitive to the motif orientation relative to the peptide ends. For the variant that displayed distinct dynamics, one orientation consistent with observed splittings corresponded to the motif being situated such that its two glycines were particularly accessible to adjacent peptides., (Copyright © 2011 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

In this study, we addressed the straightforward reconstitution of red blood cell (RBC) membranes on the surface of cationic organic polymer beads. The RBC membrane-bead complex was obtained by the incubation of white, unsealed rat RBC ghost membranes with a nonporous quaternary ammonium-type anion-exchange polymer bead with a 350-550 microm diameter. Confocal microscopic observations using a fluorescence membrane probe revealed that the RBC membranes were reconstituted on the outer surface of the bead without any remarkable structural gaps in the membrane. The absence of activity of two peripheral enzymes that latently reside on the cytoplasmic face of the RBC membranes demonstrated that the orientation of the RBC membranes immobilized on the beads was asymmetric as well as that in the native state. The RBC membrane-polymer bead complex was incubated with a primary antibody that is directed against the amino-terminal extracellular domain of the integral protein glycophorin A (GPA). The resulting complex was further incubated with a fluorescent secondary antibody and then subjected to confocal microscopic observations. Fluorescence resulting in the binding of the secondary antibody was found on the surface of the complex, which indicates that the amino-terminal extracellular domain of GPA is exposed to the surface of the complex. In addition, the anion uptake function of the most abundant integral protein anion-exchanger 1 (AE1) immobilized on the polymer beads was inhibited by pretreatment with its specific inhibitor 4,4'-diisothiocyano-2,2'-stilbene disulfonate, as is observed for the intact RBCs. Based on all these results, the RBC membranes were thought to be reconstituted on the ionic polymer beads by our one-pot procedure while maintaining the orientation and functions of the membrane proteins to some extent., (Copyright 2010 Elsevier Inc. All rights reserved.)

Protein association in lipid membranes is a complex process with thermodynamics directed by a multitude of different factors. Amino-acid sequence is a molecular parameter that affects dimerization as shown by limited directed mutations along the transmembrane domains. Membrane-mediated interactions are also important although details of such contributions remain largely unclear. In this study, we probe directly the free energy of association of Glycophorin A by means of extensive parallel Monte Carlo simulations with recently developed methods and a model that accounts for sequence-specificity while representing lipid membranes faithfully. We find that lipid-induced interactions are significant both at short and intermediate separations. The ability of molecules to tilt in a specific hydrophobic environment extends their accessible interfaces, leading to intermittent contacts during protein recognition. The dimer with the lowest free energy is largely determined by the favorable lipid-induced attractive interactions at the closest distance. Finally, the coarse-grained model employed herein, together with the extensive sampling performed, provides estimates of the free energy of association that are in excellent agreement with existing data., (Copyright 2010 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

The presence of erythropoietin (Epo) receptors on cells besides red blood cell precursors, such as cancer cells or megakaryocyte precursors, can lead to side effects during Epo therapy including enhanced tumor growth and platelet production. It would be ideal if the action of Epo could be limited to erythroid precursors. To address this issue, we constructed single-chain variable fragment (scFv)-Epo fusion proteins that used the anti-glycophorin 10F7 scFv amino-terminal to Epo analogues that would have minimal activity alone. We introduced the Epo mutations N147A, R150A and R150E, which progressively weakened receptor affinity in the context of Epo alone, as defined by cell proliferation assays using TF-1 or UT-7 cells. Fusion of these mutant proteins to the 10F7 scFv significantly rescued the activity of the mutant proteins, but had a relatively small effect on wild-type Epo. For example, fusion to the 10F7 scFv enhanced the activity of Epo(R150A) by 10- to 27-fold, while a corresponding fusion to wild-type Epo enhanced its activity only up to 2.7-fold. When glycophorin was blocked by antibody competition or reduced by siRNA-mediated inhibition of expression, the activity of 10F7 scFv-Epo(R150A) was correspondingly reduced, while such inhibition had essentially no effect on the activity of 10F7 scFv-Epo(wild-type). In addition, potent stimulation of Epo receptors by 10F7 scFv-Epo(R150A) was observed in long-term proliferation and viability assays. Taken together, these results indicate that a combination of targeting and affinity modulation can be used to engineer forms of Epo with enhanced cell-type specificity.

There is experimental evidence that the transmembrane fragment spanning amino acids 70-86 of glycophorin A, GpA70-86, forms amyloid fibrils and the inhibitor RGTFEGKF prevents GpA70-86 fibril formation at an equimolar ratio. Both the GpA70-86 and inhibitor peptides contain a GxxxG motif as found in many amyloid proteins such as the Alzheimer's amyloid beta-peptide and prion protein. To explore the intrinsic, early interaction and inhibition mechanism, we have determined the structures of GpA70-86 in the absence and presence of the inhibitor by means of extensive molecular dynamics simulations in explicit solvent. Consistent with experiments on the fibrils, our simulations show that the two GxxxG motifs interact significantly at the monomer level. They go, however, one step beyond by indicating that the inhibitor has a significant impact on the global structure of GpA70-86, but a limited influence on the conformations of the GxxxG motif. Implications of our simulations on amyloid propagation of proteins containing GxxxG motifs are discussed.

GYPC encodes two erythrocyte surface sialoglycoproteins in humans, glycophorin C and glycophorin D (GPC and GPD), via initiation of translation at two start codons on a single transcript. The malaria-causing parasite Plasmodium falciparum uses GPC as a means of invasion into the human red blood cell. Here, we examine the molecular evolution of GYPC among the Hominoidea (Greater and Lesser Apes) and also the pattern of polymorphism at the locus in a global human sample. We find an excess of nonsynonymous divergence among species that appears to be caused solely by accelerated evolution of GYPC in the human lineage. Moreover, we find that the ability of GYPC to encode both GPC and GPD is a uniquely human trait, caused by the evolution of the GPC start codon in the human lineage. The pattern of polymorphism among humans is consistent with a hitchhiking event at the locus, suggesting that positive natural selection affected GYPC in the relatively recent past. Because GPC is exploited by P. falciparum for invasion of the red blood cell, we hypothesize that selection for evasion of P. falciparum has caused accelerated evolution of GYPC in humans (relative to other primates) and that this positive selection has continued to act in the recent evolution of our species. These data suggest that malaria has played a powerful role in shaping molecules on the surface of the human red blood cell. In addition, our examination of GYPC reveals a novel mechanism of protein evolution: co-option of untranslated region (UTR) sequence following the formation of a new start codon. In the case of human GYPC, the ancestral protein (GPD) continues to be produced through leaky translation. Because leaky translation is a widespread phenomenon among genes and organisms, we suggest that co-option of UTR sequence may be an important source of protein innovation.

The GxxxG sequence motif mediates the association of transmembrane (TM) helices by providing a site of close contact between them. However, it is not sufficient for strong association. For example, both bacteriophage M13 major coat protein (MCP) and human erythrocyte protein glycophorin A (GpA) contain a GxxxG motif in their TM domains and form a homodimer, but the association affinity of MCP, measured by the ToxCAT in vivo assay, is dramatically weaker than that of GpA. Even when all interfacial residues of MCP were substituted for those of GpA (MCP-GpA), association remained significantly weaker than in GpA. Here we provide an explanation for these experimental observations using molecular dynamics simulations in an implicit membrane (IMM1-GC). The association free energies of GpA29 (GpA with 29 residues all from the wild-type sequence), GpA15p11 (GpA with 15 residues from the wild-type sequence plus 11 flanking residues from the ToxCAT construct), MCP, and MCP-GpA TM helices were calculated and compared. MCP and MCP-GpA have the same flanking residues used in the ToxCAT assay as those in GpA15p11, but the position of the flanking residues relative to the GxxxG motif is different. The calculated association free energies follow experimental observations: the association affinity of MCP-GpA falls between those of GpA15p11 and MCP wild-type. MCP exhibits an equally strong interhelical interaction in the TM domain. A major reason for the weaker association of MCP in the calculations was the noninterfacial residue Lys-40, which in the dimer structure is forced to be buried in the membrane interior. To alleviate the desolvation cost, in MCP and MCP-GpA dimers, Lys-40 gets deprotonated. A second factor that modulates association affinity is the flanking residues. Thanks to them, GpA15p11 exhibits a much stronger association affinity than GpA29. The positioning of the flanking residues is also important, as evidenced by the difference in association affinity between MCP and MCP-GpA on one hand and GpA15p11 on the other. Thus, residues outside the contact interface can exert a significant influence on transmembrane helix association affinity.

The dimerizations of membrane proteins, Outer Membrane Phospholipase A (OMPLA) and glycophorin A (GPA), have been simulated by an adapted Brownian Dynamics program. To mimic the membrane protein environment, we introduced a hybrid electrostatic potential map of membrane and water for electrostatic interaction calculations. We added a van der Waals potential term to the force field of the current version of the BD program to simulate the short-range interactions of the two monomers. We reduced the BD sampling space from three dimensions to two dimensions to improve the efficiency of BD simulations for membrane proteins. The OMPLA and GPA dimers predicted by our 2D-BD simulation and structural refinement is in good agreement with the experimental structures. The adapted 2D-BD method could be used for prediction of dimerization of other membrane proteins, such as G protein-coupled receptors, to help better understanding of the structures and functions of membrane proteins.

Glycophorins are the major sialoglycoproteins in red blood cell membranes, possessing various physiological and pathological roles. We examined membrane glycoproteins in canine red cells and cloned cDNAs for two major glycophorins, glycophorins A (GPA) and C (GPC) from bone marrow cells. Periodic acid-Schiff staining and immunoblotting analyses showed that canine red cell membranes contained several glycoproteins immunoreactive to an anti-bovine GPC antibody, whereas the most abundant sialoglycoproteins, the candidates for GPA, did not react with an anti-human GPA antibody. The amino acid sequences of the extracellular domains of GPA and GPC had no significant homology to those from other mammalian species, including humans, and had O-linked and/or N-linked glycosylation sites. On the other hand, the C-terminal cytoplasmic domain and/or the transmembrane helices of GPA and GPC were conserved among species, indicating some functional significance of those regions in red cell membranes that include dimerization of GPA in the membrane-spanning region, and association of GPC with membrane skeletal proteins through binding with protein 4.1 and p55 in the cytoplasmic domain. These findings provide insights for clinical studies to evaluate the involvement of GPA and GPC in the pathogenesis of red cell diseases.

Objective: To identify the binding site on glycophorin A (GPA) for EBA-175 to provide clue for developing short peptide vaccine and therapeutic agents against Plasmodium falciparum., Methods: With the recombinant protein of EBA-175 as the target molecule, the mimetic peptides of GPA were screened from a 12-mer random peptide library. Three rounds of biopanning were carried out, and enzyme-linked immunosorbent assay (ELISA), competitive ELISA, Dot-ELISA and Western blotting used to evaluate the binding between the phage-borne peptides and EBA-175. The insert DNA sequences of positive clones were determined and their amino acid sequences deduced., Results: Thirty clones from the third round were randomly selected, of which 27 were found positive by sandwich ELISA. Competitive ELISA proved that most of the phage-borne peptides could competitively inhibit the binding of antibody (EBA-175 Ab) with EBA-175. Analysis of DNA and amino acid sequences indicated that 24 positive phage clones contained the conservative sequence of IRR, which was highly homologous with the 114-116 amino acids of GPA., Conclusion: These phage-displayed peptides can bind with EBA-175, and the amino acid sequence IRR might play an important role in the binding between EBA-175 and GPA.

Integral membrane proteins pose a major challenge for protein-structure prediction because only approximately 100 high-resolution structures are available currently, thereby impeding the development of rules or empirical potentials to predict the packing of transmembrane alpha-helices. However, when an intermediate-resolution electron microscopy (EM) map is available, it can be used to provide restraints which, in combination with a suitable computational protocol, make structure prediction feasible. In this work we present such a protocol, which proceeds in three stages: 1), generation of an ensemble of alpha-helices by flexible fitting into each of the density rods in the low-resolution EM map, spanning a range of rotational angles around the main helical axes and translational shifts along the density rods; 2), fast optimization of side chains and scoring of the resulting conformations; and 3), refinement of the lowest-scoring conformations with internal coordinate mechanics, by optimizing the van der Waals, electrostatics, hydrogen bonding, torsional, and solvation energy contributions. In addition, our method implements a penalty term through a so-called tethering map, derived from the EM map, which restrains the positions of the alpha-helices. The protocol was validated on three test cases: GpA, KcsA, and MscL.

The tertiary interactions between amide-I vibrators on the separate helices of transmembrane helix dimers were probed by ultrafast 2D vibrational photon echo spectroscopy. The 2D IR approach proves to be a useful structural method for the study of membrane-bound structures. The 27-residue human erythrocyte protein Glycophorin A transmembrane peptide sequence: KKITLIIFG(79)VMAGVIGTILLISWG(94)IKK was labeled at G(79) and G(94) with (13)C=(16)O or (13)C=(18)O. The isotopomers and their 50:50 mixtures formed helical dimers in SDS micelles whose 2D IR spectra showed components from homodimers when both helices had either (13)C=(16)O or (13)C=(18)O substitution and a heterodimer when one had (13)C=(16)O substitution and the other had (13)C=(18)O substitution. The cross-peaks in the pure heterodimer 2D IR difference spectrum and the splitting of the homodimer peaks in the linear IR spectrum show that the amide-I mode is delocalized across a pair of helices. The excitation exchange coupling in the range 4.3-6.3 cm(-1) arises from through-space interactions between amide units on different helices. The angle between the two Gly(79) amide-I transition dipoles, estimated at 103 degrees from linear IR spectroscopy and 110 degrees from 2D IR spectroscopy, combined with the coupling led to a structural picture of the hydrophobic interface that is remarkably consistent with results from NMR on helix dimers. The helix crossing angle in SDS is estimated at 45 degrees. Two-dimensional IR spectroscopy also sets limits on the range of geometrical parameters for the helix dimers from an analysis of the coupling constant distribution.

Molecular dynamics (MD) simulations of proteins in a lipid bilayer environment are usually undertaken with one or a few starting structures. Here we report a search protocol for systematically exploring the possible interactions in helical bundle transmembrane proteins, a frequently occurring structural motif. The search protocol correctly identifies the experimentally known structure of the dimeric human glycophorin A transmembrane domain as the lowest energy structure among five different models without any prior assumptions, whilst an identical in vacuo search fails to identify the correct structure. The lowest energy structure from the search in a lipid bilayer has a root mean square deviation of 1.1A to the experimental structure. We have applied the same search protocol to the unknown transmembrane structure of the oncogenic mutant ErbB-2 protein, a member of the family of epidermal growth factor receptors. Resulting structures show the role of glutamic acid hydrogen bonding and close helical packing. Water molecules may also play a key role in stabilisation of the transmembrane helix association.

A large number of experimental studies have been devoted to quantifying the interaction between transmembrane (TM) helices in detergent micelles and, more recently, in bilayers. Theoretical calculation of association free energy of TM helices would be useful for predicting the propensity of given sequences to oligomerize and for understanding the difference between association in micelles and in bilayers. In this article, the theoretical foundation for calculating the standard association free energy of TM helices is laid out and is applied to glycophorin A in both micelles and bilayers. The standard association free energy is decomposed into the effective energy, translational, rotational, and conformational entropy terms. The effective energy of association is obtained by molecular dynamics simulations in an implicit membrane model. The translational and rotational entropy of association is calculated from the probability distribution of the translational and rotational degrees of freedom obtained from the molecular dynamics simulations. The side-chain conformational entropy of association is estimated from the probability distribution obtained by rigid rotation of all side-chain dihedral angles. The calculated standard association free energy of glycophorin A in N-dodecylphosphocholine micelles is in good agreement with the experimental value. The translational entropy cost is larger, whereas the rotational entropy cost is smaller in bilayers than in micelles. The standard association free energy in 1,2-dimyristoyl-sn-glycero-3-phosphocholine bilayers is calculated to be approximately 1.3 kcal/mol more favorable than in N-dodecylphosphocholine micelles, consistent with available experimental data.

Plasmodium vivax uses a single member of the Duffy binding-like (DBL) receptor family to invade erythrocytes and is not found in West Africa where its erythrocyte ligand, the Duffy blood group antigen, is missing. In contrast, Plasmodium falciparum expresses four members of the DBL family, and remarkably, single-point mutations of two of these receptors (BAEBL and JESEBL) bind to entirely different erythrocyte ligands, greatly expanding the range of erythrocytes that P. falciparum can invade. In this article, we describe the molecular basis of the binding specificity for one BAEBL variant (VSTK) that binds to glycophorin C. We demonstrate that soluble glycophorin C completely blocks the binding of BAEBL (VSTK) to human erythrocytes, requiring 0.7 microM for 50% inhibition, a concentration similar to that required by glycophorin A to block the binding of erythrocyte-binding antigen 175 to erythrocytes. BAEBL (VSTK) does not bind to Gerbich-negative erythrocytes that express a truncated form of glycophorin C because it lacks exon 3. The N-linked oligosaccharide of Gerbich-negative glycophorin C has a markedly different composition than the wild-type glycophorin C. Moreover, removal of the N-linked oligosaccharide from the wild-type glycophorin C eliminates its ability to inhibit binding of BAEBL (VSTK) to erythrocytes. These findings are consistent with the ligand for BAEBL (VSTK) being, in part, the N-linked oligosaccharide and suggest that single-point mutations in BAEBL allow P. falciparum to recognize oligosaccharides on different erythrocyte surface glycoproteins or glycolipids, greatly increasing its invasion range.

Self-association of two hydrophobic alpha-helices is studied via unrestrained Monte Carlo (MC) simulations in a hydrophobic slab described by an effective potential. The system under study represents two transmembrane (TM) segments of human glycophorin A (GpA), which form homo-dimers in membranes. The influence of TM electrostatic potential, thickness and hydrophobicity degree of lipid bilayer is investigated. It is shown that the membrane environment stabilizes alpha-helical conformation of GpA monomers, induces their TM insertion and facilitates inter-helical contacts. Head-to-head orientation of the helices is promoted by the voltage difference across the membrane. Subsequent "fine-tuned" assembling of the dimer is mediated by van der Waals interactions. Only the models of dimer, calculated in a hydrophobic slab with applied voltage agree with experimental data, while simulations in vacuo or without TM voltage fail to give reasonable results. The moderate structural heterogeneity of GpA dimers (existence of several groups of states with close energies) is proposed to reflect their equilibrium dynamics in membrane-mimics. The calculations performed for GpA mutants G83A and G86L permit rationalization of mutagenesis data for them. The results of Monte Carlo simulations critically depend on the parameters of the membrane model: adequate description of helix association is achieved in the water-cyclohexane-water system with the membrane thickness 30-34 A, while in membranes with different hydrophobicities and thickness unrealistic conformations of the dimer are found. The computational approach permits efficient prediction of TM helical oligomers based solely on the sequences of interacting peptides.

One of the hallmark features of the integrin receptors is the ability to transmit signals bidirectionally through the cell membrane. The transmembrane integrin domains are pivotal to the signaling events. An understanding of the signaling mechanism requires structural information. Here, we report a structural model of the transmembrane and part of the cytosolic domains of the alphaIIbbeta3 integrin in its resting state. The model was obtained computationally by a restrained conformational search of helix-helix interactions. It agrees with one published NMR structure of the cytoplasmic complex and can put many experimental findings on structural grounds. According to our model, integrins form an intricately designed coiled-coil structure in the resting state. The conserved Glycophorin A (GpA)-like sequence motif of the alpha, but not the beta, subunit, is in the interface of this model. Based on our calculations and other data, a signaling mechanism that involves a transient GpA-like structure is proposed.

The principles that govern the folding and packing of membrane proteins are still not completely understood. In the present work, we have revisited the glycophorin A (GpA) dimerisation motif that mediates transmembrane (TM) helix association, one of the best-suited models of membrane protein oligomerisation. By using artificial polyleucine TM segments we have demonstrated in this study that a pattern of only five amino acids (GVxxGVxxT) promotes specific dimerisation. Further, we have used this minimised GpA motif to assess the influence of hydrophobic matching on the TM helix packing process in detergent micelles and found that this factor modulates helix-helix association and/or dissociation between TM fragments.

Projection of transmembrane helices using a Uniform B-spline Algorithm is a tool for the visualization of interactions between helices in membrane proteins. It allows the user to generate projections of 3D helices, no matter what their deviations from a canonical helix might be. When associated with adapted coloring schemes it facilitates the comprehension of helix-helix interactions. Examples of transmembrane proteins were chosen to illustrate the advantages that this method provides. In the glycophorin A dimer we can easily appreciate the structural features behind homodimerisation. Using the structure of the fumarate reductase we analyze the contact surfaces inside a helical bundle and thanks to structures from a molecular dynamics simulation we see how modifications in structure and electrostatics relate to their interaction. We propose the use of this tool as an aid to the visualization and analysis of transmembrane helix surfaces and properties.

Insertion and formation of membrane proteins involves the interaction of protein helices with one another in lipid environments. Researchers have studied glycophorin A (GpA) transmembrane helices embedded in sodium dodecyl sulfate (SDS) micelles to identify contacts significant for helix dimerization. However, a detailed picture of the conformation and dynamics of the GpA-SDS system cannot be obtained solely through experiment. Molecular dynamics simulations of SDS and a GpA dimer can provide an atomic-level picture of SDS aggregation and helix association. We report 2.5-ns simulations of GpA wild-type and mutants in a preformed micelle as well as a 32-ns simulation showing the formation of a complete micelle around wild-type GpA from an initially random placement of SDS molecules in an aqueous environment. In the latter case, an initial instability of GpA helices in water is reversed after the helices become surrounded by SDS. The properties of the spontaneously formed micelle surrounding the GpA are indistinguishable from those of the preformed micelle surrounding the GpA dimer.

Homo- and hetero-oligomeric interactions between the transmembrane (TM) helices of integrin alpha and beta subunits may play an important role in integrin activation and clustering. As a first step to understanding these interactions, we used the TOXCAT assay to measure oligomerization of the wild-type alpha(IIb) TM helix and single-site TM domain mutants. TOXCAT measures the oligomerization of a chimeric protein containing a TM helix in the Escherichia coli inner membrane via the transcriptional activation of the gene for chloramphenicol acetyltransferase. We found the amount of chloramphenicol acetyltransferase induced by the wild-type alpha(IIb) TM helix was approximately half that induced by the strongly dimerizing TM helix of glycophorin A, confirming that the alpha(IIb) TM domain oligomerizes in biological membranes. Mutating each of the alpha(IIb) TM domain residues to either Ala, Leu, Ile, or Val revealed that a GXXXG motif mediates oligomerization. Further, we found that the residue preceding each glycine contributed to the oligomerization interface, as did the residue at position i + 4 after the second Gly of GXXXG. Thus, the sequence XXVGXXGGXXXLXX is critical for oligomerization of alpha(IIb) TM helix. These data were used to generate an atomic model of the alpha(IIb) homodimer, revealing a family of structures with right-handed crossing angles of 40 degrees to 60 degrees, consistent with a 4.0-residue periodicity, and with an interface rotated by 50 degrees relative to glycophorin A. Thus, although the alpha(IIb) TM helix makes use of the GXXXG framework, neighboring residues have evolved to engineer its dimerization interface, enabling it to subserve specific and specialized functions.

The presented work describes a structural model for integrin homooligomerization, focusing on the transmembrane domains. The two noncovalently linked integrin subunits, alpha and beta, were previously shown to homodimerize or homotrimerize, respectively. Our work is based on published mutational work that induced homotrimerization of beta3 integrins. The mutations provided structural restraints for the creation of a structural model of the beta3 homotrimer by a computational search of the conformational space of homomeric interactions of the beta3 integrin. Additionally, we explored possible conformations of the alphaIIb integrin homodimer, for which no unique solution was found. Two possible models of signal transduction, involving two different alphaIIb conformations, are discussed. One of the possible homodimeric alphaIIb conformations is GpA like, which is in line with experimental evidence. Based on our here-presented structural models and on recent experiments, we will argue that most probably the heteromeric alpha/beta transmembrane complex separates in the course of clustering.

Sequence motifs are responsible for ensuring the proper assembly of transmembrane (TM) helices in the lipid bilayer. To understand the mechanism by which the affinity of a common TM-TM interactive motif is controlled at the sequence level, we compared two well characterized GXXXG motif-containing homodimers, those formed by human erythrocyte protein glycophorin A (GpA, high-affinity dimer) and those formed by bacteriophage M13 major coat protein (MCP, low affinity dimer). In both constructs, the GXXXG motif is necessary for TM-TM association. Although the remaining interfacial residues (underlined) in GpA (LIXXGVXXGVXXT) differ from those in MCP (VVXXGAXXGIXXF), molecular modeling performed here indicated that GpA and MCP dimers possess the same overall fold. Thus, we could introduce GpA interfacial residues, alone and in combination, into the MCP sequence to help decrypt the determinants of dimer affinity. Using both in vivo TOXCAT assays and SDS-PAGE gel migration rates of synthetic peptides derived from TM regions of the proteins, we found that the most distal interfacial sites, 12 residues apart (and approximately 18 A in structural space), work in concert to control TM-TM affinity synergistically.

During erythroblast enucleation, nuclei surrounded by plasma membrane separate from erythroblast cytoplasm. A key aspect of this process is sorting of erythroblast plasma membrane components to reticulocytes and expelled nuclei. Although it is known that cytoskeletal elements actin and spectrin partition to reticulocytes, little is understood about molecular mechanisms governing plasma membrane protein sorting. We chose glycophorin A (GPA) as a model integral protein to begin investigating protein-sorting mechanisms. Using immunofluorescence microscopy and Western blotting we found that GPA sorted predominantly to reticulocytes. We hypothesized that the degree of skeletal linkage might control the sorting pattern of transmembrane proteins. To explore this hypothesis, we quantified the extent of GPA association to the cytoskeleton in erythroblasts, young reticulocytes, and mature erythrocytes using fluorescence imaged microdeformation (FIMD) and observed that GPA underwent dramatic reorganization during terminal differentiation. We discovered that GPA was more connected to the membrane cytoskeleton, either directly or indirectly, in erythroblasts and young reticulocytes than in mature cells. We conclude that skeletal protein association can regulate protein sorting during enucleation. Further, we suggest that the enhanced rigidity of reticulocyte membranes observed in earlier investigations results, at least in part, from increased connectivity of GPA with the spectrin-based skeleton.

G protein-coupled receptors (GPCRs) can form dimeric or oligomeric complexes in vivo. However, the functions and mechanisms of oligomerization remain poorly understood for most GPCRs, including the alpha-factor receptor (STE2 gene product) of the yeast Saccharomyces cerevisiae. Here we provide evidence indicating that alpha-factor receptor oligomerization involves a GXXXG motif in the first transmembrane domain (TM1), similar to the transmembrane dimerization domain of glycophorin A. Results of fluorescence resonance energy transfer, fluorescence microscopy, endocytosis assays of receptor oligomerization in living cells, and agonist binding assays indicated that amino acid substitutions affecting the glycine residues of the GXXXG motif impaired alpha-factor receptor oligomerization and biogenesis in vivo but did not significantly impair agonist binding affinity. Mutant receptors exhibited signaling defects that were not due to impaired cell surface expression, indicating that oligomerization promotes alpha-factor receptor signal transduction. Structure-function studies suggested that the GXXXG motif in TM1 of the alpha-factor receptor promotes oligomerization by a mechanism similar to that used by the GXXXG dimerization motif of glycophorin A. In many mammalian GPCRs, motifs related to the GXXXG sequence are present in TM1 or other TM domains, suggesting that similar mechanisms are used by many GPCRs to form dimers or oligomeric arrays.

Exploiting recent developments in generalized Born (GB) electrostatics theory, we have reformulated the calculation of the self-electrostatic solvation energy to account for the influence of biological membranes. Consistent with continuum Poisson-Boltzmann (PB) electrostatics, the membrane is approximated as an solvent-inaccessible infinite planar low-dielectric slab. The present membrane GB model closely reproduces the PB electrostatic solvation energy profile across the membrane. The nonpolar contribution to the solvation energy is taken to be proportional to the solvent-exposed surface area (SA) with a phenomenological surface tension coefficient. The proposed membrane GB/SA model requires minor modifications of the pre-existing GB model and appears to be quite efficient. By combining this implicit model for the solvent/bilayer environment with advanced computational sampling methods, like replica-exchange molecular dynamics, we are able to fold and assemble helical membrane peptides. We examine the reliability of this model and approach by applications to three membrane peptides: melittin from bee venom, the transmembrane domain of the M2 protein from Influenza A (M2-TMP), and the transmembrane domain of glycophorin A (GpA). In the context of these proteins, we explore the role of biological membranes (represented as a low-dielectric medium) in affecting the conformational changes in melittin, the tilt of transmembrane peptides with respect to the membrane normal (M2-TMP), helix-to-helix interactions in membranes (GpA), and the prediction of the configuration of transmembrane helical bundles (GpA). The present method is found to perform well in each of these cases and is anticipated to be useful in the study of folding and assembly of membrane proteins as well as in structure refinement and modeling of membrane proteins where a limited number of experimental observables are available.

We have examined the role of the environment on the interactions between transmembrane helices using, as a model system, the dimerization of the glycophorin A transmembrane helix. In this study we have focused on micellar environments and have examined a series of detergents that include a range of alkyl chain lengths, combined with ionic, zwitterionic, and nonionic headgroups. For each we have measured how the apparent equilibrium constant depends on the detergent concentration. In two detergents we also measured the thermal sensitivity of the equilibrium constant, from which we derive the van't Hoff enthalpy and entropy. We show that several simple models are inadequate for explaining our results; however, models that include the effect of detergent concentration on detergent binding are able to account for our measurements. Our analysis suggests that the effects of detergents on helix association are due to a pair of opposing effects: an enthalpic effect, which drives association as the detergent concentration is increased and which is sensitive to the chemical nature of the detergent headgroup, opposed by an entropic effect, which drives peptide dissociation as the detergent concentration is raised. Our results also indicate that the monomer-monomer interface is relatively hydrophilic and that association within detergent micelles is driven by the enthalpy change. The wide variations in glycophorin a dimmer, stability with the detergent used, together with the realization that this results from the balance between two opposing effects, suggests that detergents might be selected that drive association rather than dissociation of peptide dimers.