Search

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

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

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

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

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.

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.

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.

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.

Dimerization of the transmembrane domain of glycophorin A is mediated by a seven residue motif LIxxGVxxGVxxT through a combination of van der Waals and hydrogen bonding interactions. One of the unusual features of the motif is the large number of beta-branched amino acids that may limit the entropic cost of dimerization by restricting side-chain motion in the monomeric transmembrane helix. Deuterium NMR spectroscopy is used to characterize the dynamics of fully deuterated Val80 and Val84, two essential amino acids of the dimerization motif. Deuterium spectra of the glycophorin A transmembrane dimer were obtained using synthetic peptides corresponding to the transmembrane sequence containing either perdeuterated Val80 or Val84. These data were compared with spectra of monomeric glycophorin A peptides deuterated at Val84. In all cases, the deuterium line shapes are characterized by fast methyl group rotation with virtually no motion about the C(alpha)-C(beta) bond. This is consistent with restriction of the side chain in both the monomer and dimer due to intrahelical packing interactions involving the beta-methyl groups, and indicates that there is no energy cost associated with dimerization due to loss of conformational entropy. In contrast, deuterium NMR spectra of Met81 and Val82, in the lipid interface, reflected greater motional averaging and fast exchange between different side-chain conformers.

Molecular interactions between transmembrane alpha-helices can be explored using global searching molecular dynamics simulations (GSMDS), a method that produces a group of probable low energy structures. We have shown previously that the correct model in various homooligomers is always located at the bottom of one of various possible energy basins. Unfortunately, the correct model is not necessarily the one with the lowest energy according to the computational protocol, which has resulted in overlooking of this parameter in favor of experimental data. In an attempt to use energetic considerations in the aforementioned analysis, we used global searching molecular dynamics simulations on three homooligomers of different sizes, the structures of which are known. As expected, our results show that even when the conformational space searched includes the correct structure, taking together simulations using both left and right handedness, the correct model does not necessarily have the lowest energy. However, for the models derived from the simulation that uses the correct handedness, the lowest energy model is always at, or very close to, the correct orientation. We hypothesize that this should also be true when simulations are performed using homologous sequences, and consequently lowest energy models with the right handedness should produce a cluster around a certain orientation. In contrast, using the wrong handedness the lowest energy structures for each sequence should appear at many different orientations. The rationale behind this is that, although more than one energy basin may exist, basins that do not contain the correct model will shift or disappear because they will be destabilized by at least one conservative (i.e. silent) mutation, whereas the basin containing the correct model will remain. This not only allows one to point to the possible handedness of the bundle, but can be used to overcome ambiguities arising from the use of homologous sequences in the analysis of global searching molecular dynamics simulations. In addition, because clustering of lowest energy models arising from homologous sequences only happens when the estimation of the helix tilt is correct, it may provide a validation for the helix tilt estimate.

The transmembrane helix of glycophorin A contains a seven-residue motif, LIxxGVxxGVxxT, that mediates protein dimerization. Threonine is the only polar amino acid in this motif with the potential to stabilize the dimer through hydrogen-bonding interactions. Polarized Fourier transform infrared spectroscopy is used to establish a robust protocol for incorporating glycophorin A transmembrane peptides into membrane bilayers. Analysis of the dichroic ratio of the 1655-cm(-1) amide I vibration indicates that peptides reconstituted by detergent dialysis have a transmembrane orientation with a helix crossing angle of <35 degrees. Solid-state nuclear magnetic resonance spectroscopy is used to establish high resolution structural restraints on the conformation and packing of Thr-87 in the dimer interface. Rotational resonance measurement of a 2.9-A distance between the gamma-methyl and backbone carbonyl carbons of Thr-87 is consistent with a gauche- conformation for the chi1 torsion angle. Rotational-echo double-resonance measurements demonstrate close packing (4.0 +/- 0.2 A) of the Thr-87 gamma-methyl group with the backbone nitrogen of Ile-88 across the dimer interface. The short interhelical distance places the beta-hydroxyl of Thr-87 within hydrogen-bonding range of the backbone carbonyl of Val-84 on the opposing helix. These results refine the structure of the glycophorin A dimer in membrane bilayers and highlight the complementary role of small and polar residues in the tight association of transmembrane helices in membrane proteins.

To assess local elasticity in the red cell's spectrin-actin network, nano-particles were tethered to actin nodes and their constrained thermal motions were tracked. Cells were both immobilized and controllably deformed by aspiration into a micropipette. Since the network is well-appreciated as soft, thermal fluctuations even in an unstressed portion of network were expected to be many tens of nanometers based on simple equipartition ideas. Real-time particle tracking indeed reveals such root-mean-squared motions for 40-nm fluorescent beads either tethered to actin directly within a cell ghost or connected to actin from outside a cell via glycophorin. Moreover, the elastically constrained displacements are significant on the scale of the network's internodal distance of approximately 60-80 nm. Surprisingly, along the aspirated projection-where the network is axially extended by as much as twofold or more-fluctuations in the axial direction are increased by almost twofold relative to motions in the unstressed network. The molecular basis for such strain softening is discussed broadly in terms of force-driven transitions. Specific considerations are given to 1) protein dissociations that reduce network connectivity, and 2) unfolding kinetics of a localized few of the red cell's approximately 10(7) spectrin repeats.

Transmembrane helices are no longer believed to be just hydrophobic segments that exist solely to anchor proteins to a lipid bilayer, but rather they appear to have the capacity to specify function and structure. Specific interactions take place between hydrophobic segments within the lipid bilayer whereby subtle mutations that normally would be considered innocuous can result in dramatic structural differences. That such specificity takes place within the lipid bilayer implies that it may be possible to identify the most favorable interaction surface of transmembrane alpha-helices based on computational methods alone, as shown in this study. Herein, an attempt is made to map the energy surface of several transmembrane helix-helix interactions for several homo-oligomerizing proteins, where experimental data regarding their structure exist (glycophorin A, phospholamban, Influenza virus A M2, Influenza virus C CM2, and HIV vpu). It is shown that due to symmetry constraints in homo-oligomers the computational problem can be simplified. The results obtained are mostly consistent with known structural data and may additionally provide a view of possible alternate and intermediate configurations.

Solution small angle x-ray scattering can be used to study the association of transmembrane proteins solubilized in detergent micelles. We have used the alpha-helical transmembrane domain of the human erythrocyte glycophorin A (GpA) fused to the carboxyl terminus of monomeric staphylococcal nuclease (SN/GpA) as a model system for study. By matching the average electron density of the detergent micelles to that of the buffer solution, the micelle contribution to the small angle scattering vanishes, and the molecular weight and the radius of gyration of the proteins can be determined. SN/GpA has been found to dimerize in a zwitterionic detergent micelle, N-dodecyl-N,N-(dimethylammonio)butyrate (DDMAB), whose average electron density naturally matches the electron density of an aqueous buffer. The dimerization occurs through the transmembrane domains of GpA. With the aid of the nuclease domain scattering, the orientation of the helices within a dimer can be determined to be parallel by radius of gyration analysis. The association constant of a mutant (G83I) that weakens the GpA dimerization has been determined to be 24 microM in the DDMAB environment. The experimental methods established here could be used to apply solution small angle x-ray scattering to studying the association and interactions of other membrane proteins.

Cholesterol has a concentration-dependent effect on membrane organization. It is able to control the membrane permeability by inducing conformational ordering of the lipid chains. A systematic investigation of lipid bilayer permeability is described in the present work. It takes advantage of the transmembrane potential difference modulation induced in vesicles when an external electric field is applied. The magnitude of this modulation is under the control of the membrane electrical permeability. When brought to a critical value by the external field, the membrane potential difference induces a new membrane organization. The membrane is then permeable and prone to solubilized membrane protein back-insertion. This is obtained for an external field strength, which depends on membrane native permeability. This approach was used to study the cholesterol effect on phosphatidylcholine bilayers. Studies have been performed with lipids in gel and in fluid states. When cholesterol is present, it does not affect electropermeabilization and electroinsertion in lipids in the fluid state. When lipids are in the gel state, cholesterol has a dose-dependent effect. When present at 6% (mol/mol), cholesterol prevents electropermeabilization and electroinsertion. When cholesterol is present at more than 12%, electropermeabilization and electroinsertion are obtained under milder field conditions. This is tentatively explained by a cholesterol-induced alteration of the hydrophobic barrier of the bilayer core. Our results indicate that lipid membrane permeability is affected by the cholesterol content.

Oriented multilamellar systems containing phospholipids and peptides have been formed on a germanium internal reflection element. Attenuated total reflection infrared spectra have been recorded and the linear dichroism of peptide amide I and amide II bands measured. Using peptides for which the orientation had been previously studied under similar experimental conditions by 15N solid-state nuclear magnetic resonance spectroscopy, important conclusions were drawn on the approach to be used to derive secondary structure orientation in a membrane from dichroic ratios. In particular, it is shown that the influence of the film thickness and refractive index on the orientation determination can be evaluated from the value of RATRiso, i.e., the dichroic ratio of a dipole oriented at the magic angle or with isotropic mobility. A series of peptides was used to test the validity of our suggestions on various helix orientations in the membrane. These include magainin 2 and hydrophobic (hPhi20) model peptides, the transmembrane segment of glycophorin (GLY), and LAH4, a designed peptide antibiotic that changes between a transmembrane and an in-plane orientation in a pH-dependent manner.