Structure-based prediction of the stability of transmembrane helix-helix interactions: the sequence dependence of glycophorin A dimerization

A Rigorous Framework for Calculating Protein-Protein Binding Affinities in Membranes

Calculating the binding free energy of integral transmembrane (TM) proteins is crucial for understanding the mechanisms by which they recognize one another and reversibly associate. The glycophorin A (GpA) homodimer, composed of two α-helical segments, has long served as a model system for studying TM protein reversible association. The present work establishes a methodological framework for calculating the binding affinity of the GpA homodimer in the heterogeneous environment of a membrane. Our investigation carefully considered a variety of protocols, including the appropriate choice of the force field, rigorous standardization reflecting the experimental conditions, sampling algorithm, anisotropic environment, and collective variables, to accurately describe GpA dimerization via molecular dynamics-based approaches. Specifically, two strategies were explored: (i) an unrestrained potential mean force (PMF) calculation, which merely enhances sampling along the separation of the two binding partners without any restraint, and (ii) a so-called "geometrical route", whereby the α-helices are progressively separated with imposed restraints on their orientational, positional, and conformational degrees of freedom to accelerate convergence. Our simulations reveal that the simplified, unrestrained PMF approach is inadequate for the description of GpA dimerization. Instead, the geometrical route, tailored specifically to GpA in a membrane environment, yields excellent agreement with experimental data within a reasonable computational time. A dimerization free energy of -10.7 kcal/mol is obtained, in fairly good agreement with available experimental data. The geometrical route further helps elucidate how environmental forces drive association before helical interactions stabilize it. Our simulations also brought to light a distinct, long-lived spatial arrangement that potentially serves as an intermediate state during dimer formation. The methodological advances in the generalized geometrical route provide a powerful tool for accurate and efficient binding-affinity calculations of intricate TM protein complexes in inhomogeneous environments.

In Search of a Molecular View of Peptide-Lipid Interactions in Membranes

Lipid bilayer membranes are often represented as a continuous nonpolar slab with a certain thickness bounded by two more polar interfaces. Phenomena such as peptide binding to the membrane surface, folding, insertion, translocation, and diffusion are typically interpreted on the basis of this view. In this Perspective, I argue that this membrane representation as a hydrophobic continuum solvent is not adequate to understand peptide-lipid interactions. Lipids are not small compared to membrane-active peptides: their sizes are similar. Therefore, peptide diffusion needs to be understood in terms of free volume, not classical continuum mechanics; peptide solubility or partitioning in membranes cannot be interpreted in terms of hydrophobic mismatch between membrane thickness and peptide length; peptide folding and translocation, often involving cationic peptides, can only be understood if realizing that lipids adapt to the presence of peptides and the membrane may undergo considerable lipid redistribution in the process. In all of those instances, the detailed molecular interactions between the peptide residues and the lipid components are essential to understand the mechanisms involved.

Thermodynamic analysis of the GASright transmembrane motif supports energetic model of dimerization

The GASright motif, best known as the fold of the glycophorin A transmembrane dimer, is one of the most common dimerization motifs in membrane proteins, characterized by its hallmark GxxxG-like sequence motifs (GxxxG, AxxxG, GxxxS, and similar). Structurally, GASright displays a right-handed crossing angle and short interhelical distance. Contact between the helical backbones favors the formation of networks of weak hydrogen bonds between Cα-H carbon donors and carbonyl acceptors on opposing helices (Cα-H···O=C). To understand the factors that modulate the stability of GASright, we previously presented a computational and experimental structure-based analysis of 26 predicted dimers. We found that the contributions of van der Waals packing and Cα-H hydrogen bonding to stability, as inferred from the structural models, correlated well with relative dimerization propensities estimated experimentally with the in vivo assay TOXCAT. Here we test this model with a quantitative thermodynamic analysis. We used Förster resonance energy transfer (FRET) to determine the free energy of dimerization of a representative subset of seven of the 26 original TOXCAT dimers using FRET. To overcome the technical issue arising from limited sampling of the dimerization isotherm, we introduced a globally fitting strategy across a set of constructs comprising a wide range of stabilities. This strategy yielded precise thermodynamic data that show strikingly good agreement between the original propensities and ΔG° of association in detergent, suggesting that TOXCAT is a thermodynamically driven process. From the correlation between TOXCAT and thermodynamic stability, the predicted free energy for all the 26 GASright dimers was calculated. These energies correlate with the in silico ΔE scores of dimerization that were computed on the basis of their predicted structure. These findings corroborate our original model with quantitative thermodynamic evidence, strengthening the hypothesis that van der Waals and Cα-H hydrogen bond interactions are the key modulators of GASright stability.

How physical forces drive the process of helical membrane protein folding

Protein folding is a fundamental process of life with important implications throughout biology. Indeed, tens of thousands of mutations have been associated with diseases, and most of these mutations are believed to affect protein folding rather than function. Correct folding is also a key element of design. These factors have motivated decades of research on protein folding. Unfortunately, knowledge of membrane protein folding lags that of soluble proteins. This gap is partly caused by the greater technical challenges associated with membrane protein studies, but also because of additional complexities. While soluble proteins fold in a homogenous water environment, membrane proteins fold in a setting that ranges from bulk water to highly charged to apolar. Thus, the forces that drive folding vary in different regions of the protein, and this complexity needs to be incorporated into our understanding of the folding process. Here, we review our understanding of membrane protein folding biophysics. Despite the greater challenge, better model systems and new experimental techniques are starting to unravel the forces and pathways in membrane protein folding.

Complete genome of the thermophilic purple sulfur Bacterium Thermochromatium tepidum compared to Allochromatium vinosum and other Chromatiaceae

The complete genome sequence of the thermophilic purple sulfur bacterium Thermochromatium tepidum strain MC^T (DSM 3771^T) is described and contrasted with that of its mesophilic relative Allochromatium vinosum strain D (DSM 180^T) and other Chromatiaceae. The Tch. tepidum genome is a single circular chromosome of 2,958,290 base pairs with no plasmids and is substantially smaller than the genome of Alc. vinosum. The Tch. tepidum genome encodes two forms of RuBisCO and contains nifHDK and several other genes encoding a molybdenum nitrogenase but lacks a gene encoding a protein that assembles the Fe-S cluster required to form a functional nitrogenase molybdenum-iron cofactor, leaving the phototroph phenotypically Nif^-. Tch. tepidum contains genes necessary for oxidizing sulfide to sulfate as photosynthetic electron donor but is genetically unequipped to either oxidize thiosulfate as an electron donor or carry out assimilative sulfate reduction, both of which are physiological hallmarks of Alc. vinosum. Also unlike Alc. vinosum, Tch. tepidum is obligately phototrophic and unable to grow chemotrophically in darkness by respiration. Several genes present in the Alc. vinosum genome that are absent from the genome of Tch. tepidum likely contribute to the major physiological differences observed between these related purple sulfur bacteria that inhabit distinct ecological niches.

Small Residues Inhibit Homo-Dimerization of the Human Carbonic Anhydrase XII Transmembrane Domain

Amino acids with small side chains and motifs of small residues in a distance of four are rather abundant in human single-span transmembrane helices. While interaction of such helices appears to be common, the role of the small residues in mediating and/or stabilizing transmembrane helix oligomers remains mostly elusive. Yet, the mere existence of (small)xxx(small) motifs in transmembrane helices is frequently used to model dimeric TM helix structures. The single transmembrane helix of the human carbonic anhydrases XII contains a large number of amino acids with small side chains, and critical involvement of these small amino acids in dimerization of the transmembrane domain has been suggested. Using the GALLEX assay, we show here that the transmembrane domain indeed forms a strong transmembrane helix oligomer within a biological membrane. However, single or multiple mutations of small residue(s) to isoleucine almost always increased, rather than decreased, the interaction propensities. Reduction of helix flexibility and of protein–lipid contacts caused by a reduced lipid accessible surface area likely results in stabilization of helix–helix interactions within the membrane.

Structure-activity relationship of presenilin in γ-secretase-mediated intramembrane cleavage

Genetic research on familial cases of Alzheimer disease have identified presenilin (PS) as an important membrane protein in the pathomechanism of this disease. PS is the catalytic subunit of γ-secretase, which is responsible for the generation of amyloid-β peptide deposited in the brains of Alzheimer disease patients. γ-Secretase is an atypical protease composed of four membrane proteins (i.e., presenilin, nicastrin, anterior pharynx defective-1 (Aph-1), and presenilin enhancer-2 (Pen-2)) and mediates intramembrane proteolysis. Numerous investigations have been conducted toward understanding the structural features of γ-secretase components as well as the cleavage mechanism of γ-secretase. In this review, we summarize our current understanding of the structure and activity relationship of the γ-secretase complex.

Disrupting the transmembrane domain-mediated oligomerization of protein tyrosine phosphatase receptor J inhibits EGFR-driven cancer cell phenotypes

Receptor protein tyrosine phosphatases (RPTPs) play critical regulatory roles in mammalian signal transduction. However, the structural basis for the regulation of their catalytic activity is not fully understood, and RPTPs are generally not therapeutically targetable. This knowledge gap is partially due to the lack of known natural ligands or selective agonists of RPTPs. Contrary to what is known from structure-function studies of receptor tyrosine kinases (RTKs), RPTP activities have been reported to be suppressed by dimerization, which may prevent RPTPs from accessing their RTK substrates. We report here that homodimerization of protein tyrosine phosphatase receptor J (PTPRJ, also known as DEP-1) is regulated by specific transmembrane (TM) residues. We found that disrupting these interactions destabilizes homodimerization of full-length PTPRJ in cells, reduces the phosphorylation of the known PTPRJ substrate epidermal growth factor receptor (EGFR) and of other downstream signaling effectors, antagonizes EGFR-driven cell phenotypes, and promotes substrate access. We demonstrate these observations in human cancer cells using mutational studies and identified a peptide that binds to the PTPRJ TM domain and represents the first example of an allosteric agonist of RPTPs. The results of our study provide fundamental structural and functional insights into how PTPRJ activity is tuned by TM interactions in cells. Our findings also open up opportunities for developing peptide-based agents that could be used as tools to probe RPTPs' signaling mechanisms or to manage cancers driven by RTK signaling.

Biologically Active Ultra-Simple Proteins Reveal Principles of Transmembrane Domain Interactions

Specific interactions between the helical membrane-spanning domains of transmembrane proteins play central roles in the proper folding and oligomerization of these proteins. However, the relationship between the hydrophobic amino acid sequences of transmembrane domains and their functional interactions is in most cases unknown. Here, we use ultra-simple artificial proteins to systematically study the sequence basis for transmembrane domain interactions. We show that most short homopolymeric polyleucine transmembrane proteins containing single amino acid substitutions can activate the platelet-derived growth factor β receptor or the erythropoietin receptor in cultured mouse cells, resulting in cell transformation or proliferation. These proteins displayed complex patterns of activity that were markedly affected by seemingly minor sequence differences in the ultra-simple protein itself or in the transmembrane domain of the target receptor, and the effects of these sequence differences are not additive. In addition, specific leucine residues along the length of these proteins are required for activity, and the positions of these required leucines differ based on the identity and position of the central substituted amino acid. Our results suggest that these ultra-simple proteins use a variety of molecular mechanisms to activate the same target and that diversification of transmembrane domain sequences over the course of evolution minimized off-target interactions.

Predicting the helix-helix interactions from correlated residue mutations

Helix-helix interactions are crucial in the structure assembly, stability and function of helix-rich proteins including many membrane proteins. In spite of remarkable progresses over the past decades, the accuracy of predicting protein structures from their amino acid sequences is still far from satisfaction. In this work, we focused on a simpler problem, the prediction of helix-helix interactions, the results of which could facilitate practical protein structure prediction by constraining the sampling space. Specifically, we started from the noisy 2D residue contact maps derived from correlated residue mutations, and utilized ridge detection to identify the characteristic residue contact patterns for helix-helix interactions. The ridge information as well as a few additional features were then fed into a machine learning model HHConPred to predict interactions between helix pairs. In an independent test, our method achieved an F-measure of ∼60% for predicting helix-helix interactions. Moreover, although the model was trained mainly using soluble proteins, it could be extended to membrane proteins with at least comparable performance relatively to previous approaches that were generated purely using membrane proteins. All data and source codes are available at http://166.111.152.91/Downloads.html or https://github.com/dpxiong/HHConPred.

NMR Investigation of Structures of G-protein Coupled Receptor Folding Intermediates

Folding of G-protein coupled receptors (GPCRs) according to the two-stage model (Popot, J. L., and Engelman, D. M. (1990) Biochemistry 29, 4031-4037) is postulated to proceed in 2 steps: partitioning of the polypeptide into the membrane followed by diffusion until native contacts are formed. Herein we investigate conformational preferences of fragments of the yeast Ste2p receptor using NMR. Constructs comprising the first, the first two, and the first three transmembrane (TM) segments, as well as a construct comprising TM1-TM2 covalently linked to TM7 were examined. We observed that the isolated TM1 does not form a stable helix nor does it integrate well into the micelle. TM1 is significantly stabilized upon interaction with TM2, forming a helical hairpin reported previously (Neumoin, A., Cohen, L. S., Arshava, B., Tantry, S., Becker, J. M., Zerbe, O., and Naider, F. (2009) Biophys. J. 96, 3187-3196), and in this case the protein integrates into the hydrophobic interior of the micelle. TM123 displays a strong tendency to oligomerize, but hydrogen exchange data reveal that the center of TM3 is solvent exposed. In all GPCRs so-far structurally characterized TM7 forms many contacts with TM1 and TM2. In our study TM127 integrates well into the hydrophobic environment, but TM7 does not stably pack against the remaining helices. Topology mapping in microsomal membranes also indicates that TM1 does not integrate in a membrane-spanning fashion, but that TM12, TM123, and TM127 adopt predominantly native-like topologies. The data from our study would be consistent with the retention of individual helices of incompletely synthesized GPCRs in the vicinity of the translocon until the complete receptor is released into the membrane interior.

Conserved transmembrane glycine residues in the Shigella flexneri polysaccharide co-polymerase protein WzzB influence protein-protein interactions

The O antigen (Oag) component of lipopolysaccharides (LPS) is crucial for virulence and Oag chain-length regulation is controlled by the polysaccharide co-polymerase class 1 (PCP1) proteins. Crystal structure analyses indicate that structural conservation among PCP1 proteins is highly maintained, however the mechanism of Oag modal-chain-length control remains to be fully elucidated. Shigella flexneri PCP1 protein WzzBSF confers a modal-chain length of 10-17 Oag repeat units (RUs), whereas the Salmonella enterica Typhimurium PCP1 protein WzzBST confers a modal-chain length of ~16-28 Oag RUs. Both proteins share >70 % overall sequence identity and contain two transmembrane (TM1 and TM2) regions, whereby a conserved proline-glycine-rich motif overlapping the TM2 region is identical in both proteins. Conserved glycine residues within TM2 are functionally important, as glycine to alanine substitutions at positions 305 and 311 confer very short Oag modal-chain length (~2-6 Oag RUs). In this study, WzzBSF was co-expressed with WzzBST in S. flexneri and a single intermediate modal-chain length of ~11-21 Oag RUs was observed, suggesting the presence of Wzz:Wzz interactions. Interestingly, co-expression of WzzBSF with WzzBG305A/G311A conferred a bimodal LPS Oag chain length (despite over 99 % protein sequence identity), and we hypothesized that the proteins fail to interact. Co-purification assays detected His6-WzzBSF co-purifying with FLAG-tagged WzzBST but not with FLAG-tagged WzzBG305A/G311A, supporting our hypothesis. These data indicate that the conserved glycine residues in TM2 are involved in Wzz:Wzz interactions, and provide insight into key interactions that drive Oag modal length control.

His499 Regulates Dimerization and Prevents Oncogenic Activation by Asparagine Mutations of the Human Thrombopoietin Receptor

Ligand binding to the extracellular domain of the thrombopoietin receptor (TpoR) imparts a specific orientation on the transmembrane (TM) and intracellular domains of the receptors that is required for physiologic activation via receptor dimerization. To map the inactive and active dimeric orientations of the TM helices, we performed asparagine (Asn)-scanning mutagenesis of the TM domains of the murine and human TpoR. Substitution of Asn at only one position (S505N) activated the human receptor, whereas Asn substitutions at several positions activated the murine receptor. Second site mutational studies indicate that His(499) near the N terminus of the TM domain is responsible for protecting the human receptor from activation by Asn mutations. Structural studies reveal that the sequence preceding His(499) is helical in the murine receptor but non-helical in peptides corresponding to the TM domain of the inactive human receptor. The activating S505N mutation and the small molecule agonist eltrombopag both induce helix in this region of the TM domain and are associated with dimerization and activation of the human receptor. Thus, His(499) regulates the activation of human TpoR and provides additional protection against activating mutations, such as oncogenic Asn mutations in the TM domain.

Transmembrane recognition of the semaphorin co-receptors neuropilin 1 and plexin A1: coarse-grained simulations

The cancer associated class 3 semaphorins require direct binding to neuropilins and association to plexins to trigger cell signaling. Here, we address the role of the transmembrane domains of neuropilin 1 and plexin A1 for the dimerization of the two receptors by characterizing the assembly in lipid bilayers using coarse-grained molecular dynamics simulations. From experimental evidence using a two-hybrid system showing the biochemical association of the two receptors transmembrane domains, we performed molecular simulations in DOPC and POPC demonstrating spontaneously assembly to form homodimers and heterodimers with a very high propensity for right-handed packing of the helices. Inversely, left-handed packing was observed with a very low propensity. This mode of packing was observed uniquely when the plexin A1 transmembrane domain was involved in association. Potential of mean force calculations were used to predict a hierarchy of self-association for the monomers: the two neuropilin 1 transmembrane domains strongly associated, neuropilin 1 and plexin A1 transmembrane domains associated less and the two plexin A1 transmembrane domains weakly but significantly associated. We demonstrated that homodimerization and heterodimerization are driven by GxxxG motifs, and that the sequence context modulates the packing mode of the plexin A1 transmembrane domains. This work presents major advances towards our understanding of membrane signaling platforms assembly through membrane domains and provides exquisite information for the design of antagonist drugs defining a novel class of therapeutic agents.

Compensatory mutants of the bovine papillomavirus E5 protein and the platelet-derived growth factor β receptor reveal a complex direct transmembrane interaction

The 44-amino-acid E5 protein of bovine papillomavirus is a dimeric transmembrane protein that exists in a stable complex with the platelet-derived growth factor (PDGF) β receptor, causing receptor activation and cell transformation. The transmembrane domain of the PDGF β receptor is required for complex formation, but it is not known if the two proteins contact one another directly. Here, we studied a PDGF β receptor mutant containing a leucine-to-isoleucine substitution in its transmembrane domain, which prevents complex formation with the wild-type E5 protein in mouse BaF3 cells and inhibits receptor activation by the E5 protein. We selected E5 mutants containing either a small deletion or multiple substitution mutations that restored binding to the mutant PDGF β receptor, resulting in receptor activation and growth factor independence. These E5 mutants displayed lower activity with PDGF β receptor mutants containing other transmembrane substitutions in the vicinity of the original mutation, and one of them cooperated with a receptor mutant containing a distal mutation in the juxtamembrane domain. These results provide strong genetic evidence that the transmembrane domains of the E5 protein and the PDGF β receptor contact one another directly. They also demonstrate that different mutations in the E5 protein allow it to tolerate the same mutation in the PDGF β receptor transmembrane domain and that a mutation in the E5 protein can allow it to tolerate different mutations in the PDGF β receptor. Thus, the rules governing direct interactions between transmembrane helices are complex and not restricted to local interactions.

Homozygosity for a null allele of SMIM1 defines the Vel-negative blood group phenotype

The Vel antigen is present on red blood cells (RBCs) from all humans except rare Vel-negative individuals who can form antibodies to Vel in response to transfusion or pregnancy. These antibodies may cause severe hemolytic reactions in blood recipients. We combined SNP profiling and transcriptional network modeling to link the Vel-negative phenotype to SMIM1, located in a 97-kb haplotype block on chromosome 1p36. This gene encodes a previously undiscovered, evolutionarily conserved transmembrane protein expressed on RBCs. Notably, 35 of 35 Vel-negative individuals were homozygous for a frameshift deletion of 17 bp in exon 3. Functional studies using antibodies raised against SMIM1 peptides confirmed a null phenotype in RBC membranes, and SMIM1 overexpression induced Vel expression. Genotype screening estimated that ~1 of 17 Swedish blood donors is a heterozygous deletion carrier and ~1 of 1,200 is a homozygous deletion knockout and enabled identification of Vel-negative donors. Our results establish SMIM1 as a new erythroid gene and Vel as a new blood group system.

Structural differences between thermophilic and mesophilic membrane proteins

The evolutionary adaptations of thermophilic water-soluble proteins required for maintaining stability at high temperature have been extensively investigated. Little is known about the adaptations in membrane proteins, however. Here, we compare many properties of mesophilic and thermophilic membrane protein structures, including side-chain burial, packing, hydrogen bonding, transmembrane kinks, loop lengths, hydrophobicity, and other sequence features. Most of these properties are quite similar between mesophiles and thermophiles although we observe a slight increase in side-chain burial and possibly a slight decrease in the frequency of transmembrane kinks in thermophilic membrane protein structures. The most striking difference is the increased hydrophobicity of thermophilic transmembrane helices, possibly reflecting more stringent hydrophobicity requirements for membrane partitioning at high temperature. In agreement with prior work examining transmembrane sequences, we find that thermophiles have an increase in small residues (Gly, Ala, Ser, and Val) and a strong suppression of Cys. We also find a relative dearth of most strongly polar residues (Asp, Asn, Glu, Gln, and Arg). These results suggest that in thermophiles, there is significant evolutionary pressure to offload destabilizing polar amino acids, to decrease the entropy cost of side chain burial, and to eliminate thermally sensitive amino acids.

Computational design of membrane proteins

This article reviews the recent successes of computational protein design techniques applied to integral membrane proteins. This emerging area is still handicapped by significant difficulties in the experimental characterization of the stability and structure of the designed proteins. Nevertheless, by focusing on oligomeric complexes of single-span transmembrane (TM) peptides with detectable activity, the computational design of membrane proteins has already produced very exciting results. The 'take-home message' is that optimization of van der Waals packing and hydrogen bonding (both 'canonical' and weak Cα-H⋯O bonds) can produce functional structures of remarkable stability and specificity in the membrane.

Toward understanding protocell mechanosensation

Mechanosensitive (MS) channels can prevent bacterial bursting during hypo-osmotic shocks by responding to increases in lateral tension at the membrane level through an integrated and coordinated opening mechanism. Mechanical regulation in protocells could have been one of the first mechanisms to evolve in order to preserve their integrity against changing environmental conditions. How has the rich functional diversity found in present cells been created throughout evolution, and what did the primordial MS channels look like? This review has been written with the aim of identifying which factors may have been important for the appearance of the first osmotic valve in a prebiotic context, and what this valve may have been like. It highlights the mechanical properties of lipid bilayers, the association of peptides as aggregates in membranes, and the conservation of sequence motifs as central aspects to understand the evolution of proteins that gate below the tension required for spontaneous pore formation and membrane rupture. The arguments developed here apply to both MscL and MscS homologs, but could be valid to mechano-susceptible proteins in general.

Role of transmembrane segment 5 and extracellular loop 3 in the homodimerization of human ABCC1

Resistance to multiple anticancer agents is a major obstacle in the successful treatment of cancers. Overexpression of some ATP-binding cassette (ABC) membrane transporters such as ABCC1 has been shown to be a major contributor of multidrug resistance (MDR) in both laboratory cell line models and the clinical setting. ABCC1 has been thought to function as a homodimer with a putative dimerization domain located in the first 281 amino acid residues, including MSD0 and L0 domains. In this study, we further mapped in detail the dimerization site and placed it in TM5 and ECL3 in MSD0 using co-expression and co-immunoprecipitation of a series of deletion constructs. TM5 and ECL3 in one subunit appear to interact with TM5 and ECL3 in the opposing subunit in a sequence-independent manner, but their physical location together with the hydrophobicity of TM5 and the length of ECL3 appears to be important contributors to the dimerization ability of ABCC1.

Lipid-modulated sequence-specific association of glycophorin A in membranes

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.

The role of individual amino acids in the dimerization of CR4 and ACR4 transmembrane domains

CRINKLY4 is a growth factor-like plant receptor kinase designated as CR4 in Zea mays and ACR4 in Arabidopsis. Using the TOXCAT system, a genetic assay that measures helix interactions in a natural membrane environment, we have previously demonstrated that the dimerization potential of the ACR4 transmembrane (TM) domain is significantly weaker than that of CR4 TM domain, even though 13 of the 24 residues are identical. Neither of the TM domains contain the GxxxG motif that has been shown to be important for the dimerization of the TM segments of several receptors. To further investigate the relationship between protein sequence and dimerization potential, we (a) mutated each of the 11 differing residues in the CR4 TM domain to the corresponding residue of ACR4 (b) made reciprocal mutations in ACR4 and (c) made hybrids consisting of half CR4 and half ACR4 TM domains. Our results suggest that most mutations in ACR4 or CR4 TM domains have low to moderate effects on the dimerization potential and that residues in the N-terminal half of the CR4 TM domain are important for dimerization.

Protein folding in membranes

Separation of cells and organelles by bilayer membranes is a fundamental principle of life. Cellular membranes contain a baffling variety of proteins, which fulfil vital functions as receptors and signal transducers, channels and transporters, motors and anchors. The vast majority of membrane-bound proteins contain bundles of alpha-helical transmembrane domains. Understanding how these proteins adopt their native, biologically active structures in the complex milieu of a membrane is therefore a major challenge in today's life sciences. Here, we review recent progress in the folding, unfolding and refolding of alpha-helical membrane proteins and compare the molecular interactions that stabilise proteins in lipid bilayers. We also provide a critical discussion of a detergent denaturation assay that is increasingly used to determine membrane-protein stability but is not devoid of conceptual difficulties.

Single-spanning transmembrane domains in cell growth and cell-cell interactions: More than meets the eye?

As a whole, integral membrane proteins represent about one third of sequenced genomes, and more than 50% of currently available drugs target membrane proteins, often cell surface receptors. Some membrane protein classes, with a defined number of transmembrane (TM) helices, are receiving much attention because of their great functional and pharmacological importance, such as G protein-coupled receptors possessing 7 TM segments. Although they represent roughly half of all membrane proteins, bitopic proteins (with only 1 TM helix) have so far been less well characterized. Though they include many essential families of receptors, such as adhesion molecules and receptor tyrosine kinases, many of which are excellent targets for biopharmaceuticals (peptides, antibodies, et al.). A growing body of evidence suggests a major role for interactions between TM domains of these receptors in signaling, through homo and heteromeric associations, conformational changes, assembly of signaling platforms, etc. Significantly, mutations within single domains are frequent in human disease, such as cancer or developmental disorders. This review attempts to give an overview of current knowledge about these interactions, from structural data to therapeutic perspectives, focusing on bitopic proteins involved in cell signaling.

Helix-helix interaction patterns in membrane proteins

Membrane-spanning α-helices represent major sites of protein-protein interaction in membrane protein oligomerization and folding. As such, these interactions may be of exquisite specificity. Specificity often rests on a complex interplay of different types of residues forming the helix-helix interfaces via dense packing and different non-covalent forces, including van der Waal’s forces, hydrogen bonding, charge-charge interactions, and aromatic interactions. These interfaces often contain complex residue motifs where the contribution of constituent amino acids depends on the context of the surrounding sequence. Moreover, transmembrane helix-helix interactions are increasingly recognized as being dynamic and dependent on the functional state of a given protein.

Intermonomer hydrogen bonds enhance GxxxG-driven dimerization of the BNIP3 transmembrane domain: roles for sequence context in helix-helix association in membranes

We determined the sequence dependence of human BNIP3 transmembrane domain dimerization using the biological assay TOXCAT. Mutants in which intermonomer hydrogen bonds between Ser172 and His173 are abolished show moderate interaction, indicating that side-chain hydrogen bonds contribute to dimer stability but are not essential to dimerization. Mutants in which a GxxxG motif composed of Gly180 and Gly184 has been abolished show little or no interaction, demonstrating the critical nature of the GxxxG motif to BNIP3 dimerization. These findings show that side-chain hydrogen bonds can enhance the intrinsic dimerization of a GxxxG motif and that sequence context can control how hydrogen bonds influence helix-helix interactions in membranes. The dimer interface mapped by TOXCAT mutagenesis agrees closely with the interfaces observed in the NMR structure and inferred from mutational analysis of dimerization on SDS-PAGE, showing that the native dimer structure is retained in detergents. We show that TOXCAT and SDS-PAGE give complementary and consistent information about BNIP3 transmembrane domain dimerization: TOXCAT is insensitive to mutations that have modest effects on self-association in detergents but readily discriminates among mutations that completely disrupt detergent-resistant dimerization. The close agreement between conclusions reached from TOXCAT and SDS-PAGE data for BNIP3 suggests that accurate estimates of the relative effects of mutations on native-state protein-protein interactions can be obtained even when the detergent environment is strongly disruptive.

Similar energetic contributions of packing in the core of membrane and water-soluble proteins

A major driving force for water-soluble protein folding is the hydrophobic effect, but membrane proteins cannot make use of this stabilizing contribution in the apolar core of the bilayer. It has been proposed that membrane proteins compensate by packing more efficiently. We therefore investigated packing contributions experimentally by observing the energetic and structural consequences of cavity creating mutations in the core of a membrane protein. We observed little difference in the packing energetics of water and membrane soluble proteins. Our results imply that other mechanisms are employed to stabilize the structure of membrane proteins.

Structural imperatives impose diverse evolutionary constraints on helical membrane proteins

The amino acid sequences of transmembrane regions of helical membrane proteins are highly constrained, diverging at slower rates than their extramembrane regions and than water-soluble proteins. Moreover, helical membrane proteins seem to fall into fewer families than water-soluble proteins. The reason for the differential restrictions on sequence remains unexplained. Here, we show that the evolution of transmembrane regions is slowed by a previously unrecognized structural constraint: Transmembrane regions bury more residues than extramembrane regions and soluble proteins. This fundamental feature of membrane protein structure is an important contributor to the differences in evolutionary rate and to an increased susceptibility of the transmembrane regions to disease-causing single-nucleotide polymorphisms.

Mixing and matching detergents for membrane protein NMR structure determination

One major obstacle to membrane protein structure determination is the selection of a detergent micelle that mimics the native lipid bilayer. Currently, detergents are selected by exhaustive screening because the effects of protein-detergent interactions on protein structure are poorly understood. In this study, the structure and dynamics of an integral membrane protein in different detergents is investigated by nuclear magnetic resonance (NMR) and electron paramagnetic resonance (EPR) spectroscopy and small-angle X-ray scattering (SAXS). The results suggest that matching of the micelle dimensions to the protein's hydrophobic surface avoids exchange processes that reduce the completeness of the NMR observations. Based on these dimensions, several mixed micelles were designed that improved the completeness of NMR observations. These findings provide a basis for the rational design of mixed micelles that may advance membrane protein structure determination by NMR.

Lactose permease as a paradigm for membrane transport proteins (Review)

Our structural knowledge of the major facilitator superfamily (MFS) has dramatically increased in the past year with three structures of proteins from the MFS (oxalate/formate antiporter; lactose/proton symporter and the P(i)/glycerol-3-phosphate antiporter). All three structures revealed 12 transmembrane helices forming two distinct domains and could imply that members of the MFS have preserved both secondary as well as tertiary structural elements during evolution. Lactose permease, a particularly well-studied member of the MFS, has been extensively explored by a number of molecular biological, biochemical and biophysical approaches. In this review, we take a closer look at the structure of LacY and incorporate a wealth of biochemical and biophysical data in order to propose a possible mechanism for lactose/proton symport. In addition, we make some brief comparisons between the structures of LacY and GlpT.

Interaction and conformational dynamics of membrane-spanning protein helices

Within 1 or 2 decades, the reputation of membrane-spanning alpha-helices has changed dramatically. Once mostly regarded as dull membrane anchors, transmembrane domains are now recognized as major instigators of protein-protein interaction. These interactions may be of exquisite specificity in mediating assembly of stable membrane protein complexes from cognate subunits. Further, they can be reversible and regulatable by external factors to allow for dynamic changes of protein conformation in biological function. Finally, these helices are increasingly regarded as dynamic domains. These domains can move relative to each other in different functional protein conformations. In addition, small-scale backbone fluctuations may affect their function and their impact on surrounding lipid shells. Elucidating the ways by which these intricate structural features are encoded by the amino acid sequences will be a fascinating subject of research for years to come.

RHYTHM--a server to predict the orientation of transmembrane helices in channels and membrane-coils

RHYTHM is a web server that predicts buried versus exposed residues of helical membrane proteins. Starting from a given protein sequence, secondary and tertiary structure information is calculated by RHYTHM within only a few seconds. The prediction applies structural information from a growing data base of precalculated packing files and evolutionary information from sequence patterns conserved in a representative dataset of membrane proteins ('Pfam-domains'). The program uses two types of position specific matrices to account for the different geometries of packing in channels and transporters ('channels') or other membrane proteins ('membrane-coils'). The output provides information on the secondary structure and topology of the protein and specifically on the contact type of each residue and its conservation. This information can be downloaded as a graphical file for illustration, a text file for analysis and statistics and a PyMOL file for modeling purposes. The server can be freely accessed at: URL: http://proteinformatics.de/rhythm.

Fast predictions of thermodynamics and kinetics of protein-protein recognition from structures: from molecular design to systems biology

The increasing call for an overall picture of the interactions between the components of a biological system that give rise to the observed function is often summarized by the expression systems biology. Both the interpretative and predictive capabilities of holistic models of biochemical systems, however, depend to a large extent on the level of physico-chemical knowledge of the individual molecular interactions making up the network. This review is focused on the structure-based quantitative characterization of protein-protein interactions, ubiquitous in any biochemical pathway. Recently developed, fast and effective computational methods are reviewed, which allow the assessment of kinetic and thermodynamic features of the association-dissociation processes of protein complexes, both in water soluble and membrane environments. The performance and the accuracy of fast and semi-empirical structure-based methods have reached comparable levels with respect to the classical and more elegant molecular simulations. Nevertheless, the broad accessibility and lower computational cost provide the former methods with the advantageous possibility to perform systems-level analyses including extensive in silico mutagenesis screenings and large-scale structural predictions of multiprotein complexes.

Amyloid-beta peptide (Abeta) neurotoxicity is modulated by the rate of peptide aggregation: Abeta dimers and trimers correlate with neurotoxicity

Alzheimer's disease is an age-related neurodegenerative disorder with its toxicity linked to the generation of amyloid-beta peptide (Abeta). Within the Abeta sequence, there is a systemic repeat of a GxxxG motif, which theoretical studies have suggested may be involved in both peptide aggregation and membrane perturbation, processes that have been implicated in Abeta toxicity. We synthesized modified Abeta peptides, substituting glycine for leucine residues within the GxxxG repeat motif (GSL peptides). These GSL peptides undergo beta-sheet and fibril formation at an increased rate compared with wild-type Abeta. The accelerated rate of amyloid fibril formation resulted in a decrease in the presence of small soluble oligomers such as dimeric and trimeric forms of Abeta in solution, as detected by mass spectrometry. This reduction in the presence of small soluble oligomers resulted in reduced binding to lipid membranes and attenuated toxicity for the GSL peptides. The potential role that dimer and trimer species binding to lipid plays in Abeta toxicity was further highlighted when it was observed that annexin V, a protein that inhibits Abeta toxicity, specifically inhibited Abeta dimers from binding to lipid membranes.

Structural and functional characterization of nonstructural protein 2 for its role in hepatitis C virus assembly

The hepatitis C virus (HCV) is a flavivirus replicating in the cytoplasm of infected cells. The HCV genome is a single-stranded RNA encoding a polyprotein that is cleaved by cellular and viral proteases into 10 different products. While the structural proteins core protein, envelope protein 1 (E1) and E2 build up the virus particle, most nonstructural (NS) proteins are required for RNA replication. One of the least studied proteins is NS2, which is composed of a C-terminal cytosolic protease domain and a highly hydrophobic N-terminal domain. It is assumed that the latter is composed of three trans-membrane segments (TMS) that tightly attach NS2 to intracellular membranes. Taking advantage of a system to study HCV assembly in a hepatoma cell line, in this study we performed a detailed characterization of NS2 with respect to its role for virus particle assembly. In agreement with an earlier report ( Jones, C. T., Murray, C. L., Eastman, D. K., Tassello, J., and Rice, C. M. (2007) J. Virol. 81, 8374-8383 ), we demonstrate that the protease domain, but not its enzymatic activity, is required for infectious virus production. We also show that serine residue 168 in NS2, implicated in the phosphorylation and stability of this protein, is dispensable for virion formation. In addition, we determined the NMR structure of the first TMS of NS2 and show that the N-terminal segment (amino acids 3-11) forms a putative flexible helical element connected to a stable alpha-helix (amino acids 12-21) that includes an absolutely conserved helix side in genotype 1b. By using this structure as well as the amino acid conservation as a guide for a functional study, we determined the contribution of individual amino acid residues in TMS1 for HCV assembly. We identified several residues that are critical for virion formation, most notably a central glycine residue at position 10 of TMS1. Finally, we demonstrate that mutations in NS2 blocking HCV assembly can be rescued by trans-complementation.

Free energy calculations applied to membrane proteins

Selected applications of free energy calculations to the realm of membrane proteins are reviewed. The theoretical underpinnings of these calculations are described, focusing on free energy perturbation and the use of thermodynamic integration to determine free energy changes along well-delineated order parameters. Current strategies for improving the reliability of free energy calculations, while making them somewhat more affordable are outlined. Application of the free energy methodology to understand the structure and function of membrane proteins is illustrated in three concrete examples: The binding of an agonist ligand to a G protein-coupled receptor, the assisted transport of a small permeant through a membrane channel, and the recognition and association of transmembrane alpha-helical domains.

Subunit b-dimer of the Escherichia coli ATP synthase can form left-handed coiled-coils

One remaining challenge to our understanding of the ATP synthase concerns the dimeric coiled-coil stator subunit b of bacterial synthases. The subunit b-dimer has been implicated in important protein interactions that appear necessary for energy conservation and that may be instrumental in energy conservation during rotary catalysis by the synthase. Understanding the stator structure and its interactions with the rest of the enzyme is crucial to the understanding of the overall catalytic mechanism. Controversy exists on whether subunit b adopts a classic left-handed or a presumed right-handed dimeric coiled-coil and whether or not staggered pairing between nonhomologous residues in the homodimer is required for intersubunit packing. In this study we generated molecular models of the Escherichia coli subunit b-dimer that were based on the well-established heptad-repeat packing exhibited by left-handed, dimeric coiled-coils by employing simulated annealing protocols with structural restraints collected from known structures. In addition, we attempted to create hypothetical right-handed coiled-coil models and left- and right-handed models with staggered packing in the coiled-coil domains. Our analyses suggest that the available structural and biochemical evidence for subunit b can be accommodated by classic left-handed, dimeric coiled-coil quaternary structures.

Transmembrane domain interactions control biological functions of neuropilin-1

Neuropilin-1 (NRP1) is a transmembrane receptor playing a pivotal role in the control of semaphorins and VEGF signaling pathways. The exact mechanism controlling semaphorin receptor complex formation is unknown. A structural analysis and modeling of NRP1 revealed a putative dimerization GxxxG motif potentially important for NRP1 dimerization and oligomerization. Our data show that this motif mediates the dimerization of the transmembrane domain of NRP1 as demonstrated by a dimerization assay (ToxLuc assay) performed in natural membrane and FRET analysis. A synthetic peptide derived from the transmembrane segment of NRP1 abolished the inhibitory effect of Sema3A. This effect depends on the capacity of the peptide to interfere with NRP1 dimerization and the formation of oligomeric complexes. Mutation of the GxxxG dimerization motif in the transmembrane domain of NRP1 confirmed its biological importance for Sema3A signaling. Overall, our results shed light on an essential step required for semaphorin signaling and provide novel evidence for the crucial role of transmembrane domain of bitopic protein containing GxxxG motif in the formation of receptor complexes that are a prerequisite for cell signaling.

A mutational study of transmembrane helix–helix interactions

Diverse methods have been developed and applied in the recent years to study interaction of transmembrane alpha-helices and often interaction of single transmembrane helices is followed on SDS-gels. Here we compare two measurements of the stability of a transmembrane helix-helix interaction, and the stability of the PsbF transmembrane helix dimer was determined in a biological membrane as well as in SDS. The observations described in this study demonstrate that the environment, in which a transmembrane helix interaction is studied, can be very critical and detergent properties can significantly influence transmembrane helix interactions, especially, when the transmembrane domain contains strongly polar residues.

Comparison of stability predictions and simulated unfolding of rhodopsin structures

Developing a better mechanistic understanding of membrane protein folding is urgently needed because of the discovery of an increasing number of human diseases, where membrane protein instability and misfolding is involved. Towards this goal, we investigated folding and stability of 7-transmembrane (TM) helical bundles by computational methods. We compared the results of three different algorithms for predicting changes in stability of proteins against an experimental mutation dataset obtained for bacteriorhodopsin (BR) and mammalian rhodopsin and find that 61.6% and 70.6% of the mutation results can potentially be explained by known local contributors to the stability of the folded state of BR and mammalian rhodopsin, respectively. To obtain further information on the predicted folding pathway of 7-TM proteins, we conducted simulated thermal unfolding experiments of all available rhodopsin structures with resolution better than 3 angstroms using the Floppy Inclusions and Rigid Substructure Topography (FIRST) method (Jacobs, D. J., A. J. Rader, L. A. Kuhn and M. F. Thorpe [2001] Proteins 44, 150) described previously for a single mammalian rhodopsin structure (Rader et al. [2004] PNAS 101, 7246). In statistical comparison we found that structures of mammalian rhodopsin have a stability core that is characterized by long-range interactions involving amino acids close in space but distant in sequence comprising positions from both extracellular loop and TM regions. In contrast, BR-simulated unfolding does not reveal such a core but is dominated by interactions within individual and groups of TM helices, consistent with the two-stage hypothesis of membrane protein folding. Similar results were obtained for halo- and sensory rhodopsins as for BRs. However, the average folding core energies of sensory rhodopsins were in between those observed for mammalian rhodopsins and BRs hinting at a possible evolution of these structures toward a rhodopsin-like behavior. These results support the conclusion that although the two-stage model can explain the mechanisms of folding and stability of BR, it fails to account for the folding and stability of mammalian rhodopsin, even though the two proteins are structurally related.

Aggregation of polyalanine in a hydrophobic environment

The dimerization of polyalanine peptides in a hydrophobic environment was explored using replica exchange molecular dynamics simulations. A nonpolar solvent (cyclohexane) was used to mimic, among other hydrophobic environments, the hydrophobic interior of a membrane in which the peptides are fully embedded. Our simulations reveal that while the polyalanine monomer preferentially adopts a beta-hairpin conformation, dimeric phases exist in an equilibrium between random coil, alpha-helical, beta-sheet, and beta-hairpin states. A thermodynamic characterization of the dimeric phases reveals that electric dipole-dipole interactions and optimal side-chain packing stabilize alpha-helical conformations, while hydrogen bond interactions favor beta-sheet conformations. Possible pathways leading to the formation of alpha-helical and beta-sheet dimers are discussed.

Molecular dynamics simulation approach for the prediction of transmembrane helix–helix heterodimers assembly

Computational methods are useful to identify favorable structures of transmembrane (TM) helix oligomers when experimental data are not available or when they cannot help to interpret helix-helix association. We report here a global search method using molecular dynamics (MD) simulations to predict the structures of transmembrane homo and heterodimers. The present approach is based only on sequence information without any experimental data and is first applied to glycophorin A to validate the protocol and to the HER2-HER3 heterodimer receptor. The method successfully reproduces the experimental structures of the TM domain of glycophorin A (GpA(TM)) with a root mean square deviation of 1.5 A. The search protocol identifies three energetically stable models of the TM domain of HER2-HER3 receptor with favorable helix-helix arrangement, including right-handed and left-handed coiled-coils. The predicted TM structures exhibit the GxxxG-like motif at the dimer interface which is presumed to drive receptor oligomerization. We demonstrate that native structures of TM domain can be predicted without quantitative experimental data. This search protocol could help to predict structures of the TM domain of HER heterodimer family.

In Silico Screening of Mutational Effects on Transmembrane Helix Dimerization: Insights from Rigid-Body Docking and Molecular Dynamics Simulations

In this study, a docking-based protocol has been probed for its ability to predict the effects of 32 single and double mutations on glycophorin A (GpA) homodimerization. Rigid-body docking simulations have been paralleled by molecular dynamics (MD) simulations in implicit membrane. The rigid-body docking-based approach proved effective in reconstituting the native architecture of the GpA dimer for the wild type and the wild-type-like mutants. The good correlative models between the intermolecular interaction descriptors derived both by rigid-body docking simulations and by MD simulations and experimental relative free energies support the assumption that the mutation-induced changes in the association free energy of GpA helices are essentially ascribed to differences in the packing interactions, whereas almost all the variations in the entropic contributions to the association free energy are constant and/or negligible. The MD-based models achieved provide insights into the structural determinants for disruptive and restoring mutational effects. The computational approaches presented in this study are fast and effective, and constitute simple and promising tools for in silico screening of mutational effects on the association properties of integral membrane proteins.