Why are polar residues within the membrane core evolutionary conserved?

Experimental and computational approaches for membrane protein insertion and topology determination

Membrane proteins play pivotal roles in a wide array of cellular processes and constitute approximately a quarter of the protein-coding genes across all organisms. Despite their ubiquity and biological significance, our understanding of these proteins remains notably less comprehensive compared to their soluble counterparts. This disparity in knowledge can be attributed, in part, to the inherent challenges associated with employing specialized techniques for the investigation of membrane protein insertion and topology. This review will center on a discussion of molecular biology methodologies and computational prediction tools designed to elucidate the insertion and topology of helical membrane proteins.

Diverse genetic causes of amenorrhea in an ethnically homogeneous cohort and an evolving approach to diagnosis

RESEARCH QUESTION

Premature ovarian insufficiency (POI) is characterised by amenorrhea associated with elevated follicle stimulating hormone (FSH) under the age of 40 years and affects 1-3.7% women. Genetic factors explain 20-30% of POI cases, but most causes remain unknown despite genomic advancements.

DESIGN

We used whole exome sequencing (WES) in four Iranian families, validated variants via Sanger sequencing, and conducted the Acyl-cLIP assay to measure HHAT enzyme activity.

RESULTS

Despite ethnic homogeneity, WES revealed diverse genetic causes, including a novel homozygous nonsense variant in SYCP2L, impacting synaptonemal complex (SC) assembly, in the first family. Interestingly, the second family had two independent causes for amenorrhea - the mother had POI due to a novel homozygous loss-of-function variant in FANCM (required for chromosomal stability) and her daughter had primary amenorrhea due to a novel homozygous GNRHR (required for gonadotropic signalling) frameshift variant. WES analysis also provided cytogenetic insights. WES revealed one individual was in fact 46, XY and had a novel homozygous missense variant of uncertain significance in HHAT, potentially responsible for complete sex reversal although functional assays did not support impaired HHAT activity. In the remaining individual, WES indicated likely mosaic Turners with the majority of X chromosome variants having an allelic balance of ∼85% or ∼15%. Microarray validated the individual had 90% 45,XO.

CONCLUSIONS

This study demonstrates the diverse causes of amenorrhea in a small, isolated ethnic cohort highlighting how a genetic cause in one individual may not clarify familial cases. We propose that, in time, genomic sequencing may become a single universal test required for the diagnosis of infertility conditions such as POI.

Effects of stand-alone polar residue on membrane protein stability and structure

Helical membrane proteins generally have a hydrophobic nature, with apolar side chains comprising the majority of the transmembrane (TM) helices. However, whenever polar side chains are present in the TM domain, they often exert a crucial role in structural interactions with other polar residues, such as TM helix associations and oligomerization. Moreover, polar residues in the TM region also often participate in protein functions, such as the Schiff base bonding between Lys residues and retinal in rhodopsin-like membrane proteins. Although many studies have focused on these functional polar residues, our understanding of stand-alone polar residues that are energetically unfavored in TM helixes is limited. Here, we adopted bacteriorhodopsin (bR) as a model system and systematically mutated 17 of its apolar Leu or Phe residues to polar Asn. Stability measurements of the resulting mutants revealed that all of these polar substitutions reduced bR stability to various extents, and the extent of destabilization of each mutant bR is also correlated to different structural factors, such as the relative accessible surface area and membrane depth of the mutation site. Structural analyses of these Asn residues revealed that they form sidechain-to-backbone hydrogen bonds that alleviate the unfavorable energetics in hydrophobic and apolar surroundings. Our results indicate that membrane proteins are able to accommodate certain stand-alone polar residues in the TM region without disrupting overall structures.

An ATP13A1-assisted topogenesis pathway for folding multi-spanning membrane proteins

Many multi-spanning membrane proteins contain poorly hydrophobic transmembrane domains (pTMDs) protected from phospholipid in mature structure. Nascent pTMDs are difficult for translocon to recognize and insert. How pTMDs are discerned and packed into mature, muti-spanning configuration remains unclear. Here, we report that pTMD elicits a post-translational topogenesis pathway for its recognition and integration. Using six-spanning protein adenosine triphosphate-binding cassette transporter G2 (ABCG2) and cultured human cells as models, we show that ABCG2's pTMD2 can pass through translocon into the endoplasmic reticulum (ER) lumen, yielding an intermediate with inserted yet mis-oriented downstream TMDs. After translation, the intermediate recruits P5A-ATPase ATP13A1, which facilitates TMD re-orientation, allowing further folding and the integration of the remaining lumen-exposed pTMD2. Depleting ATP13A1 or disrupting pTMD-characteristic residues arrests intermediates with mis-oriented and exposed TMDs. Our results explain how a "difficult" pTMD is co-translationally skipped for insertion and post-translationally buried into the final correct structure at the late folding stage to avoid excessive lipid exposure.

G Protein-Coupled Receptors Regulated by Membrane Potential

G protein-coupled receptors (GPCRs) are involved in a vast majority of signal transduction processes. Although they span the cell membrane, they have not been considered to be regulated by the membrane potential. Numerous studies over the last two decades have demonstrated that several GPCRs, including muscarinic, adrenergic, dopaminergic, and glutamatergic receptors, are voltage regulated. Following these observations, an effort was made to elucidate the molecular basis for this regulatory effect. In this review, we will describe the advances in understanding the voltage dependence of GPCRs, the suggested molecular mechanisms that underlie this phenomenon, and the possible physiological roles that it may play.

Intramembrane client recognition potentiates the chaperone functions of calnexin

One-third of the human proteome is comprised of membrane proteins, which are particularly vulnerable to misfolding and often require folding assistance by molecular chaperones. Calnexin (CNX), which engages client proteins via its sugar-binding lectin domain, is one of the most abundant ER chaperones, and plays an important role in membrane protein biogenesis. Based on mass spectrometric analyses, we here show that calnexin interacts with a large number of nonglycosylated membrane proteins, indicative of additional nonlectin binding modes. We find that calnexin preferentially bind misfolded membrane proteins and that it uses its single transmembrane domain (TMD) for client recognition. Combining experimental and computational approaches, we systematically dissect signatures for intramembrane client recognition by calnexin, and identify sequence motifs within the calnexin TMD region that mediate client binding. Building on this, we show that intramembrane client binding potentiates the chaperone functions of calnexin. Together, these data reveal a widespread role of calnexin client recognition in the lipid bilayer, which synergizes with its established lectin-based substrate binding. Molecular chaperones thus can combine different interaction modes to support the biogenesis of the diverse eukaryotic membrane proteome.

Protein Design: From the Aspect of Water Solubility and Stability

Water solubility and structural stability are key merits for proteins defined by the primary sequence and 3D-conformation. Their manipulation represents important aspects of the protein design field that relies on the accurate placement of amino acids and molecular interactions, guided by underlying physiochemical principles. Emulated designer proteins with well-defined properties both fuel the knowledge-base for more precise computational design models and are used in various biomedical and nanotechnological applications. The continuous developments in protein science, increasing computing power, new algorithms, and characterization techniques provide sophisticated toolkits for solubility design beyond guess work. In this review, we summarize recent advances in the protein design field with respect to water solubility and structural stability. After introducing fundamental design rules, we discuss the transmembrane protein solubilization and de novo transmembrane protein design. Traditional strategies to enhance protein solubility and structural stability are introduced. The designs of stable protein complexes and high-order assemblies are covered. Computational methodologies behind these endeavors, including structure prediction programs, machine learning algorithms, and specialty software dedicated to the evaluation of protein solubility and aggregation, are discussed. The findings and opportunities for Cryo-EM are presented. This review provides an overview of significant progress and prospects in accurate protein design for solubility and stability.

Experimental determination and data-driven prediction of homotypic transmembrane domain interfaces

•

Homotypic TMD interfaces identified by different techniques share strong similarities.

•

The GxxxG motif is the feature most strongly associated with interfaces.

•

Other features include conservation, polarity, coevolution, and depth in the membrane

•

The role of each of each feature strongly depends on the individual protein.

•

Machine-learning helps predict interfaces from evolutionary sequence data

Interactions between their transmembrane domains (TMDs) frequently support the assembly of single-pass membrane proteins to non-covalent complexes. Yet, the TMD-TMD interactome remains largely uncharted. With a view to predicting homotypic TMD-TMD interfaces from primary structure, we performed a systematic analysis of their physical and evolutionary properties. To this end, we generated a dataset of 50 self-interacting TMDs. This dataset contains interfaces of nine TMDs from bitopic human proteins (Ire1, Armcx6, Tie1, ATP1B1, PTPRO, PTPRU, PTPRG, DDR1, and Siglec7) that were experimentally identified here and combined with literature data. We show that interfacial residues of these homotypic TMD-TMD interfaces tend to be more conserved, coevolved and polar than non-interfacial residues. Further, we suggest for the first time that interface positions are deficient in β-branched residues, and likely to be located deep in the hydrophobic core of the membrane. Overrepresentation of the GxxxG motif at interfaces is strong, but that of (small)xxx(small) motifs is weak. The multiplicity of these features and the individual character of TMD-TMD interfaces, as uncovered here, prompted us to train a machine learning algorithm. The resulting prediction method, THOIPA (www.thoipa.org), excels in the prediction of key interface residues from evolutionary sequence data.

Conserved Glu-47 and Lys-50 residues are critical for UDP-N-acetylglucosamine/UMP antiport activity of the mouse Golgi-associated transporter Slc35a3

Nucleotide sugar transporters (NSTs) regulate the flux of activated sugars from the cytosol into the lumen of the Golgi apparatus where glycosyltransferases use them for the modification of proteins, lipids, and proteoglycans. It has been well-established that NSTs are antiporters that exchange nucleotide sugars with the respective nucleoside monophosphate. Nevertheless, information about the molecular basis of ligand recognition and transport is scarce. Here, using topology predictors, cysteine-scanning mutagenesis, expression of GFP-tagged protein variants, and phenotypic complementation of the yeast strain Kl3, we identified residues involved in the activity of a mouse UDP-GlcNAc transporter, murine solute carrier family 35 member A3 (mSlc35a3). We specifically focused on the putative transmembrane helix 2 (TMH2) and observed that cells expressing E47C or K50C mSlc35a3 variants had lower levels of GlcNAc-containing glycoconjugates than WT cells, indicating impaired UDP-GlcNAc transport activity of these two variants. A conservative substitution analysis revealed that single or double substitutions of Glu-47 and Lys-50 do not restore GlcNAc glycoconjugates. Analysis of mSlc35a3 and its genetic variants reconstituted into proteoliposomes disclosed the following: (i) all variants act as UDP-GlcNAc/UMP antiporters; (ii) conservative substitutions (E47D, E47Q, K50R, or K50H) impair UDP-GlcNAc uptake; and (iii) substitutions of Glu-47 and Lys-50 dramatically alter kinetic parameters, consistent with a critical role of these two residues in mSlc35a3 function. A bioinformatics analysis revealed that an EXXK motif in TMH2 is highly conserved across SLC35 A subfamily members, and a 3D-homology model predicted that Glu-47 and Lys-50 are facing the central cavity of the protein.

Folding and Misfolding of Human Membrane Proteins in Health and Disease: From Single Molecules to Cellular Proteostasis

Advances over the past 25 years have revealed much about how the structural properties of membranes and associated proteins are linked to the thermodynamics and kinetics of membrane protein (MP) folding. At the same time biochemical progress has outlined how cellular proteostasis networks mediate MP folding and manage misfolding in the cell. When combined with results from genomic sequencing, these studies have established paradigms for how MP folding and misfolding are linked to the molecular etiologies of a variety of diseases. This emerging framework has paved the way for the development of a new class of small molecule "pharmacological chaperones" that bind to and stabilize misfolded MP variants, some of which are now in clinical use. In this review, we comprehensively outline current perspectives on the folding and misfolding of integral MPs as well as the mechanisms of cellular MP quality control. Based on these perspectives, we highlight new opportunities for innovations that bridge our molecular understanding of the energetics of MP folding with the nuanced complexity of biological systems. Given the many linkages between MP misfolding and human disease, we also examine some of the exciting opportunities to leverage these advances to address emerging challenges in the development of therapeutics and precision medicine.

A network of chaperones prevents and detects failures in membrane protein lipid bilayer integration

A fundamental step in membrane protein biogenesis is their integration into the lipid bilayer with a defined orientation of each transmembrane segment. Despite this, it remains unclear how cells detect and handle failures in this process. Here we show that single point mutations in the membrane protein connexin 32 (Cx32), which cause Charcot-Marie-Tooth disease, can cause failures in membrane integration. This leads to Cx32 transport defects and rapid degradation. Our data show that multiple chaperones detect and remedy this aberrant behavior: the ER-membrane complex (EMC) aids in membrane integration of low-hydrophobicity transmembrane segments. If they fail to integrate, these are recognized by the ER-lumenal chaperone BiP. Ultimately, the E3 ligase gp78 ubiquitinates Cx32 proteins, targeting them for degradation. Thus, cells use a coordinated system of chaperones for the complex task of membrane protein biogenesis, which can be compromised by single point mutations, causing human disease.

A Ceramide-Regulated Element in the Late Endosomal Protein LAPTM4B Controls Amino Acid Transporter Interaction

Membrane proteins are functionally regulated by the composition of the surrounding lipid bilayer. The late endosomal compartment is a central site for the generation of ceramide, a bioactive sphingolipid, which regulates responses to cell stress. The molecular interactions between ceramide and late endosomal transmembrane proteins are unknown. Here, we uncover in atomistic detail the ceramide interaction of Lysosome Associated Protein Transmembrane 4B (LAPTM4B), implicated in ceramide-dependent cell death and autophagy, and its functional relevance in lysosomal nutrient signaling. The ceramide-mediated regulation of LAPTM4B depends on a sphingolipid interaction motif and an adjacent aspartate residue in the protein's third transmembrane (TM3) helix. The interaction motif provides the preferred contact points for ceramide while the neighboring membrane-embedded acidic residue confers flexibility that is subject to ceramide-induced conformational changes, reducing TM3 bending. This facilitates the interaction between LAPTM4B and the amino acid transporter heavy chain 4F2hc, thereby controlling mTORC signaling. These findings provide mechanistic insights into how transmembrane proteins sense and respond to ceramide.

Helix formation and stability in membranes

In this article we review current understanding of basic principles for the folding of membrane proteins, focusing on the more abundant alpha-helical class. Membrane proteins, vital to many biological functions and implicated in numerous diseases, fold into their active conformations in the complex environment of the cell bilayer membrane. While many membrane proteins rely on the translocon and chaperone proteins to fold correctly, others can achieve their functional form in the absence of any translation apparatus or other aides. Nevertheless, the spontaneous folding process is not well understood at the molecular level. Recent findings suggest that helix fraying and loop formation may be important for overall structure, dynamics and regulation of function. Several types of membrane helices with ionizable amino acids change their topology with pH. Additionally we note that some peptides, including many that are rich in arginine, and a particular analogue of gramicidin, are able passively to translocate across cell membranes. The findings indicate that a final protein structure in a lipid-bilayer membrane is sequence-based, with lipids contributing to stability and regulation. While much progress has been made toward understanding the folding process for alpha-helical membrane proteins, it remains a work in progress. This article is part of a Special Issue entitled: Emergence of Complex Behavior in Biomembranes edited by Marjorie Longo.

Glucagon-Like Peptide-1 and Its Class B G Protein-Coupled Receptors: A Long March to Therapeutic Successes

The glucagon-like peptide (GLP)-1 receptor (GLP-1R) is a class B G protein-coupled receptor (GPCR) that mediates the action of GLP-1, a peptide hormone secreted from three major tissues in humans, enteroendocrine L cells in the distal intestine, α cells in the pancreas, and the central nervous system, which exerts important actions useful in the management of type 2 diabetes mellitus and obesity, including glucose homeostasis and regulation of gastric motility and food intake. Peptidic analogs of GLP-1 have been successfully developed with enhanced bioavailability and pharmacological activity. Physiologic and biochemical studies with truncated, chimeric, and mutated peptides and GLP-1R variants, together with ligand-bound crystal structures of the extracellular domain and the first three-dimensional structures of the 7-helical transmembrane domain of class B GPCRs, have provided the basis for a two-domain-binding mechanism of GLP-1 with its cognate receptor. Although efforts in discovering therapeutically viable nonpeptidic GLP-1R agonists have been hampered, small-molecule modulators offer complementary chemical tools to peptide analogs to investigate ligand-directed biased cellular signaling of GLP-1R. The integrated pharmacological and structural information of different GLP-1 analogs and homologous receptors give new insights into the molecular determinants of GLP-1R ligand selectivity and functional activity, thereby providing novel opportunities in the design and development of more efficacious agents to treat metabolic disorders.

Identifying ionic interactions within a membrane using BLaTM, a genetic tool to measure homo- and heterotypic transmembrane helix-helix interactions

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.

Key interactions by conserved polar amino acids located at the transmembrane helical boundaries in Class B GPCRs modulate activation, effector specificity and biased signalling in the glucagon-like peptide-1 receptor

Class B GPCRs can activate multiple signalling effectors with the potential to exhibit biased agonism in response to ligand stimulation. Previously, we highlighted key TM domain polar amino acids that were crucial for the function of the GLP-1 receptor, a key therapeutic target for diabetes and obesity. Using a combination of mutagenesis, pharmacological characterisation, mathematical and computational molecular modelling, this study identifies additional highly conserved polar residues located towards the TM helical boundaries of Class B GPCRs that are important for GLP-1 receptor stability and/or controlling signalling specificity and biased agonism. This includes (i) three positively charged residues (R3.30227, K4.64288, R5.40310) located at the extracellular boundaries of TMs 3, 4 and 5 that are predicted in molecular models to stabilise extracellular loop 2, a crucial domain for ligand affinity and receptor activation; (ii) a predicted hydrogen bond network between residues located in TMs 2 (R2.46176), 6 (R6.37348) and 7 (N7.61406 and E7.63408) at the cytoplasmic face of the receptor that is important for stabilising the inactive receptor and directing signalling specificity, (iii) residues at the bottom of TM 5 (R5.56326) and TM6 (K6.35346 and K6.40351) that are crucial for receptor activation and downstream signalling; (iv) residues predicted to be involved in stabilisation of TM4 (N2.52182 and Y3.52250) that also influence cell signalling. Collectively, this work expands our understanding of peptide-mediated signalling by the GLP-1 receptor.

Ionization Properties of Histidine Residues in the Lipid Bilayer Membrane Environment

We address the critically important ionization properties of histidine side chains of membrane proteins, when exposed directly to lipid acyl chains within lipid bilayer membranes. The problem is important for addressing general principles that may underlie membrane protein function. To this end, we have employed a favorable host peptide framework provided by GWALP23 (acetyl-GGALW(5)LALALALALALALW(19)LAGA-amide). We inserted His residues into position 12 or 14 of GWALP23 (replacing either Leu(12) or Leu(14)) and incorporated specific [(2)H]Ala labels within the helical core sequence. Solid-state (2)H NMR spectra report the folding and orientation of the core sequence, revealing marked differences in the histidine-containing transmembrane helix behavior between acidic and neutral pH conditions. At neutral pH, the GWALP23-H12 and GWALP23-H14 helices exhibit well defined tilted transmembrane orientations in dioleoylphosphatidylcholine (DOPC)and dilauroylphosphatidylcholine (DLPC) bilayer membranes. Under acidic conditions, when His(12) is protonated and charged, the GWALP23-H12 helix exhibits a major population that moves to the DOPC bilayer surface and a minor population that occupies multiple transmembrane states. The response to protonation of His(14) is an increase in helix tilt, but GWALP23-H14 remains in a transmembrane orientation. The results suggest pKa values of less than 3 for His(12) and about 3-5 for His(14) in DOPC membranes. In the thinner DLPC bilayers, with increased water access, the helices are less responsive to changes in pH. The combined results enable us to compare the ionization properties of lipid-exposed His, Lys, and Arg side chains in lipid bilayer membranes.

Type IV pilins regulate their own expression via direct intramembrane interactions with the sensor kinase PilS

Type IV pili are important virulence factors for many pathogens, including Pseudomonas aeruginosa Transcription of the major pilin gene-pilA-is controlled by the PilS-PilR two-component system in response to unknown signals. The absence of a periplasmic sensing domain suggested that PilS may sense an intramembrane signal, possibly PilA. We suggest that direct interactions between PilA and PilS in the inner membrane reduce pilA transcription when PilA levels are high. Overexpression in trans of PilA proteins with diverse and/or truncated C termini decreased native pilA transcription, suggesting that the highly conserved N terminus of PilA was the regulatory signal. Point mutations in PilA or PilS that disrupted their interaction prevented autoregulation of pilA transcription. A subset of PilA point mutants retained the ability to interact with PilS but could no longer decrease pilA transcription, suggesting that interaction between the pilin and sensor kinase is necessary but not sufficient for pilA autoregulation. Furthermore, PilS's phosphatase motif was required for the autoregulation of pilA transcription, suggesting that under conditions where PilA is abundant, the PilA-PilS interaction promotes PilR dephosphorylation and thus down-regulation of further pilA transcription. These data reveal a clever bacterial inventory control strategy in which the major subunit of an important P. aeruginosa virulence factor controls its own expression.

Allostery and Biased Agonism at Class B G Protein-Coupled Receptors

Class B G protein-coupled receptors (GPCRs) respond to paracrine or endocrine peptide hormones involved in control of bone homeostasis, glucose regulation, satiety, and gastro-intestinal function, as well as pain transmission. These receptors are targets for existing drugs that treat osteoporosis, hypercalcaemia, Paget's disease, type II diabetes, and obesity and are being actively pursued as targets for numerous other diseases. Exploitation of class B receptors has been limited by difficulties with small molecule drug discovery and development and an under appreciation of factors governing optimal therapeutic efficacy. Recently, there has been increasing awareness of novel attributes of GPCR function that offer new opportunity for drug development. These include the presence of allosteric binding sites on the receptor that can be exploited as drug binding pockets and the ability of individual drugs to enrich subpopulations of receptor conformations to selectively control signaling, a phenomenon termed biased agonism. In this review, current knowledge of biased signaling and small molecule allostery within class B GPCRs is discussed, highlighting areas that have progressed significantly over the past decade, in addition to those that remain largely unexplored with respect to these phenomena.

Accurate Prediction of Contact Numbers for Multi-Spanning Helical Membrane Proteins

Prediction of the three-dimensional (3D) structures of proteins by computational methods is acknowledged as an unsolved problem. Accurate prediction of important structural characteristics such as contact number is expected to accelerate the otherwise slow progress being made in the prediction of 3D structure of proteins. Here, we present a dropout neural network-based method, TMH-Expo, for predicting the contact number of transmembrane helix (TMH) residues from sequence. Neuronal dropout is a strategy where certain neurons of the network are excluded from back-propagation to prevent co-adaptation of hidden-layer neurons. By using neuronal dropout, overfitting was significantly reduced and performance was noticeably improved. For multi-spanning helical membrane proteins, TMH-Expo achieved a remarkable Pearson correlation coefficient of 0.69 between predicted and experimental values and a mean absolute error of only 1.68. In addition, among those membrane protein-membrane protein interface residues, 76.8% were correctly predicted. Mapping of predicted contact numbers onto structures indicates that contact numbers predicted by TMH-Expo reflect the exposure patterns of TMHs and reveal membrane protein-membrane protein interfaces, reinforcing the potential of predicted contact numbers to be used as restraints for 3D structure prediction and protein-protein docking. TMH-Expo can be accessed via a Web server at www.meilerlab.org .

A Hydrogen-Bonded Polar Network in the Core of the Glucagon-Like Peptide-1 Receptor Is a Fulcrum for Biased Agonism: Lessons from Class B Crystal Structures

The glucagon-like peptide 1 (GLP-1) receptor is a class B G protein-coupled receptor (GPCR) that is a key target for treatments for type II diabetes and obesity. This receptor, like other class B GPCRs, displays biased agonism, though the physiologic significance of this is yet to be elucidated. Previous work has implicated R2.60(190), N3.43(240), Q7.49(394), and H6.52(363) as key residues involved in peptide-mediated biased agonism, with R2.60(190), N3.43(240), and Q7.49(394) predicted to form a polar interaction network. In this study, we used novel insight gained from recent crystal structures of the transmembrane domains of the glucagon and corticotropin releasing factor 1 (CRF1) receptors to develop improved models of the GLP-1 receptor that predict additional key molecular interactions with these amino acids. We have introduced E6.53(364)A, N3.43(240)Q, Q7.49(394)N, and N3.43(240)Q/Q7.49(394)N mutations to probe the role of predicted H-bonding and charge-charge interactions in driving cAMP, calcium, or extracellular signal-regulated kinase (ERK) signaling. A polar interaction between E6.53(364) and R2.60(190) was predicted to be important for GLP-1- and exendin-4-, but not oxyntomodulin-mediated cAMP formation and also ERK1/2 phosphorylation. In contrast, Q7.49(394), but not R2.60(190)/E6.53(364) was critical for calcium mobilization for all three peptides. Mutation of N3.43(240) and Q7.49(394) had differential effects on individual peptides, providing evidence for molecular differences in activation transition. Collectively, this work expands our understanding of peptide-mediated signaling from the GLP-1 receptor and the key role that the central polar network plays in these events.

SCMMTP: identifying and characterizing membrane transport proteins using propensity scores of dipeptides

BACKGROUND

Identifying putative membrane transport proteins (MTPs) and understanding the transport mechanisms involved remain important challenges for the advancement of structural and functional genomics. However, the transporter characters are mainly acquired from MTP crystal structures which are hard to crystalize. Therefore, it is desirable to develop bioinformatics tools for the effective large-scale analysis of available sequences to identify novel transporters and characterize such transporters.

RESULTS

This work proposes a novel method (SCMMTP) based on the scoring card method (SCM) using dipeptide composition to identify and characterize MTPs from an existing dataset containing 900 MTPs and 660 non-MTPs which are separated into a training dataset consisting 1,380 proteins and an independent dataset consisting 180 proteins. The SCMMTP produced estimating propensity scores for amino acids and dipeptides as MTPs. The SCMMTP training and test accuracy levels respectively reached 83.81% and 76.11%. The test accuracy of support vector machine (SVM) using a complicated classification method with a low possibility for biological interpretation and position-specific substitution matrix (PSSM) as a protein feature is 80.56%, thus SCMMTP is comparable to SVM-PSSM. To identify MTPs, SCMMTP is applied to three datasets including: 1) human transmembrane proteins, 2) a photosynthetic protein dataset, and 3) a human protein database. MTPs showing α-helix rich structure is agreed with previous studies. The MTPs used residues with low hydration energy. It is hypothesized that, after filtering substrates, the hydrated water molecules need to be released from the pore regions.

CONCLUSIONS

SCMMTP yields estimating propensity scores for amino acids and dipeptides as MTPs, which can be used to identify novel MTPs and characterize transport mechanisms for use in further experiments.

AVAILABILITY

http://iclab.life.nctu.edu.tw/iclab_webtools/SCMMTP/.

Marginally hydrophobic transmembrane α-helices shaping membrane protein folding

Cells have developed an incredible machinery to facilitate the insertion of membrane proteins into the membrane. While we have a fairly good understanding of the mechanism and determinants of membrane integration, more data is needed to understand the insertion of membrane proteins with more complex insertion and folding pathways. This review will focus on marginally hydrophobic transmembrane helices and their influence on membrane protein folding. These weakly hydrophobic transmembrane segments are by themselves not recognized by the translocon and therefore rely on local sequence context for membrane integration. How can such segments reside within the membrane? We will discuss this in the light of features found in the protein itself as well as the environment it resides in. Several characteristics in proteins have been described to influence the insertion of marginally hydrophobic helices. Additionally, the influence of biological membranes is significant. To begin with, the actual cost for having polar groups within the membrane may not be as high as expected; the presence of proteins in the membrane as well as characteristics of some amino acids may enable a transmembrane helix to harbor a charged residue. The lipid environment has also been shown to directly influence the topology as well as membrane boundaries of transmembrane helices-implying a dynamic relationship between membrane proteins and their environment.

The safety dance: biophysics of membrane protein folding and misfolding in a cellular context

Most biological processes require the production and degradation of proteins, a task that weighs heavily on the cell. Mutations that compromise the conformational stability of proteins place both specific and general burdens on cellular protein homeostasis (proteostasis) in ways that contribute to numerous diseases. Efforts to elucidate the chain of molecular events responsible for diseases of protein folding address one of the foremost challenges in biomedical science. However, relatively little is known about the processes by which mutations prompt the misfolding of α-helical membrane proteins, which rely on an intricate network of cellular machinery to acquire and maintain their functional structures within cellular membranes. In this review, we summarize the current understanding of the physical principles that guide membrane protein biogenesis and folding in the context of mammalian cells. Additionally, we explore how pathogenic mutations that influence biogenesis may differ from those that disrupt folding and assembly, as well as how this may relate to disease mechanisms and therapeutic intervention. These perspectives indicate an imperative for the use of information from structural, cellular, and biochemical studies of membrane proteins in the design of novel therapeutics and in personalized medicine.

Mechanisms of integral membrane protein insertion and folding

The biogenesis, folding, and structure of α-helical membrane proteins (MPs) are important to understand because they underlie virtually all physiological processes in cells including key metabolic pathways, such as the respiratory chain and the photosystems, as well as the transport of solutes and signals across membranes. Nearly all MPs require translocons--often referred to as protein-conducting channels--for proper insertion into their target membrane. Remarkable progress toward understanding the structure and functioning of translocons has been made during the past decade. Here, we review and assess this progress critically. All available evidence indicates that MPs are equilibrium structures that achieve their final structural states by folding along thermodynamically controlled pathways. The main challenge for cells is the targeting and membrane insertion of highly hydrophobic amino acid sequences. Targeting and insertion are managed in cells principally by interactions between ribosomes and membrane-embedded translocons. Our review examines the biophysical and biological boundaries of MP insertion and the folding of polytopic MPs in vivo. A theme of the review is the under-appreciated role of basic thermodynamic principles in MP folding and assembly. Thermodynamics not only dictates the final folded structure but also is the driving force for the evolution of the ribosome-translocon system of assembly. We conclude the review with a perspective suggesting a new view of translocon-guided MP insertion.

Influence of Pathogenic Mutations on the Energetics of Translocon-Mediated Bilayer Integration of Transmembrane Helices

Aberrant protein folding and assembly contribute to a number of diseases, and efforts to rationalize how pathogenic mutations cause this phenomenon represent an important imperative in biochemical research. However, for α-helical membrane proteins, this task is complicated by the fact that membrane proteins require intricate machinery to achieve structural and functional maturity under cellular conditions. In this work, we utilized the ΔG predictor algorithm ( www.dgpred.cbr.su.se ) to survey 470 known pathogenic mutations occurring in five misfolding-prone α-helical membrane proteins for their predicted effects on the translocon-mediated membrane integration of transmembrane helices, a critical step in biosynthesis and folding of nascent membrane proteins. The results suggest that about 10 % of these mutations are likely to have adverse effects on the topogenesis of nascent membrane proteins. These results suggest that the misfolding of a modest but nonetheless significant subset of pathogenic variants may begin at the translocon. Potential implications for therapeutic design and personalized medicine are discussed.

Folding of Aquaporin 1: multiple evidence that helix 3 can shift out of the membrane core

The folding of most integral membrane proteins follows a two-step process: initially, individual transmembrane helices are inserted into the membrane by the Sec translocon. Thereafter, these helices fold to shape the final conformation of the protein. However, for some proteins, including Aquaporin 1 (AQP1), the folding appears to follow a more complicated path. AQP1 has been reported to first insert as a four-helical intermediate, where helix 2 and 4 are not inserted into the membrane. In a second step, this intermediate is folded into a six-helical topology. During this process, the orientation of the third helix is inverted. Here, we propose a mechanism for how this reorientation could be initiated: first, helix 3 slides out from the membrane core resulting in that the preceding loop enters the membrane. The final conformation could then be formed as helix 2, 3, and 4 are inserted into the membrane and the reentrant regions come together. We find support for the first step in this process by showing that the loop preceding helix 3 can insert into the membrane. Further, hydrophobicity curves, experimentally measured insertion efficiencies and MD-simulations suggest that the barrier between these two hydrophobic regions is relatively low, supporting the idea that helix 3 can slide out of the membrane core, initiating the rearrangement process.

Why is the biological hydrophobicity scale more accurate than earlier experimental hydrophobicity scales?

The recognition of transmembrane helices by the translocon is primarily guided by the average hydrophobicity of the potential transmembrane helix. However, the exact hydrophobicity of each amino acid can be identified in several different ways. The free energy of transfer for amino acid analogues between a hydrophobic media, for example, octanol and water can be measured or obtained from simulations, the hydrophobicity can also be estimated by statistical properties from known transmembrane segments and finally the contribution of each amino acid type for the probability of translocon recognition has recently been measured directly. Although these scales correlate quite well, there are clear differences between them and it is not well understood which scale represents neither the biology best nor what the differences are. Here, we try to provide some answers to this by studying the ability of different scales to recognize transmembrane helices and predict the topology of transmembrane proteins. From this analysis it is clear that the biological hydrophobicity scale as well scales created from statistical analysis of membrane helices perform better than earlier experimental scales that are mainly based on measurements of amino acid analogs and not directly on transmembrane helix recognition. Using these results we identified the properties of the scales that perform better than other scales. We find, for instance, that the better performing scales consider proline more hydrophilic. This shows that transmembrane recognition is not only governed by pure hydrophobicity but also by the helix preferences for amino acids, as proline is a strong helix breaker.

Quality control of integral membrane proteins by assembly-dependent membrane integration

Cell-surface multiprotein complexes are synthesized in the endoplasmic reticulum (ER), where they undergo cotranslational membrane integration and assembly. The quality control mechanisms that oversee these processes remain poorly understood. We show that less hydrophobic transmembrane (TM) regions derived from several single-pass TM proteins can enter the ER lumen completely. Once mislocalized, they are recognized by the Hsp70 chaperone BiP. In a detailed analysis for one of these proteins, the αβT cell receptor (αβTCR), we show that unassembled ER-lumenal subunits are rapidly degraded, whereas specific subunit interactions en route to the native receptor promote membrane integration of the less hydrophobic TM segments, thereby stabilizing the protein. For the TCR α chain, both complete ER import and subunit assembly depend on the same pivotal residue in its TM region. Thus, membrane integration linked to protein assembly allows cellular quality control of membrane proteins and connects the lumenal ER chaperone machinery to membrane protein biogenesis.

GPCR activation: protonation and membrane potential

GPCR proteins represent the largest family of signaling membrane proteins in eukaryotic cells. Their importance to basic cell biology, human diseases, and pharmaceutical interventions is well established. Many crystal structures of GPCR proteins have been reported in both active and inactive conformations. These data indicate that agonist binding alone is not sufficient to trigger the conformational change of GPCRs necessary for binding of downstream G-proteins, yet other essential factors remain elusive. Based on analysis of available GPCR crystal structures, we identified a potential conformational switch around the conserved Asp2.50, which consistently shows distinct conformations between inactive and active states. Combining the structural information with the current literature, we propose an energy-coupling mechanism, in which the interaction between a charge change of the GPCR protein and the membrane potential of the living cell plays a key role for GPCR activation.

Quantitative residue-level structure-evolution relationships in the yeast membrane proteome

Membrane proteins exist in distinctly different environments than do soluble proteins, resulting in differences between their respective biophysical and evolutionary properties. In comparison with soluble proteins, relatively little is known about how the unique biophysical properties of membrane proteins affect their evolutionary properties at the residue level. In particular, transmembrane (TM) regions of membrane proteins tend to be more conserved than regions outside of the membrane (extramembrane [EM] regions), but the mechanisms underlying this phenomenon are not well understood. Here, we combine homology-based high-resolution three-dimensional protein models with rigorous evolutionary rate calculations to quantitatively assess residue-level structure-evolution relationships in the yeast membrane proteome. We find that residue evolutionary rate increases linearly with decreasing residue burial, regardless of the hydrophobic or hydrophilic nature of the solvent environment. This finding supports a direct relationship between a residue's selective constraint and the extent of its packing interactions with neighboring residues, independent of hydrophobic effects. Most importantly, for a fixed degree of burial, residues from TM regions tend to evolve more slowly than residues from EM regions. We attribute this difference to the increased importance of packing constraints and the decreased importance of hydrophobic effects in TM regions. This additional selective constraint on TM residues plays a dominant role in explaining why TM regions evolve more slowly than EM regions. In addition to revealing the universality of the linear relationship between residue burial and selective constraint across solvent environments, our work highlights the distinct residue-level evolutionary consequences imposed by the unique biophysical properties of the membrane environment.

A computationally designed water-soluble variant of a G-protein-coupled receptor: the human mu opioid receptor

G-protein-coupled receptors (GPCRs) play essential roles in various physiological processes, and are widely targeted by pharmaceutical drugs. Despite their importance, studying GPCRs has been problematic due to difficulties in isolating large quantities of these membrane proteins in forms that retain their ligand binding capabilities. Creating water-soluble variants of GPCRs by mutating the exterior, transmembrane residues provides a potential method to overcome these difficulties. Here we present the first study involving the computational design, expression and characterization of water-soluble variant of a human GPCR, the human mu opioid receptor (MUR), which is involved in pain and addiction. An atomistic structure of the transmembrane domain was built using comparative (homology) modeling and known GPCR structures. This structure was highly similar to the subsequently determined structure of the murine receptor and was used to computationally design 53 mutations of exterior residues in the transmembrane region, yielding a variant intended to be soluble in aqueous media. The designed variant expressed in high yield in Escherichia coli and was water soluble. The variant shared structural and functionally related features with the native human MUR, including helical secondary structure and comparable affinity for the antagonist naltrexone (K_d  = 65 nM). The roles of cholesterol and disulfide bonds on the stability of the receptor variant were also investigated. This study exemplifies the potential of the computational approach to produce water-soluble variants of GPCRs amenable for structural and functionally related characterization in aqueous solution.

Polar transmembrane interactions drive formation of ligand-specific and signal pathway-biased family B G protein-coupled receptor conformations

Recently, the concept of ligand-directed signaling--the ability of different ligands of an individual receptor to promote distinct patterns of cellular response--has gained much traction in the field of drug discovery, with the potential to sculpt biological response to favor therapeutically beneficial signaling pathways over those leading to harmful effects. However, there is limited understanding of the mechanistic basis underlying biased signaling. The glucagon-like peptide-1 receptor is a major target for treatment of type-2 diabetes and is subject to ligand-directed signaling. Here, we demonstrate the importance of polar transmembrane residues conserved within family B G protein-coupled receptors, not only for protein folding and expression, but also in controlling activation transition, ligand-biased, and pathway-biased signaling. Distinct clusters of polar residues were important for receptor activation and signal preference, globally changing the profile of receptor response to distinct peptide ligands, including endogenous ligands glucagon-like peptide-1, oxyntomodulin, and the clinically used mimetic exendin-4.

Polar/Ionizable residues in transmembrane segments: effects on helix-helix packing

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.

A model-structure of a periplasm-facing state of the NhaA antiporter suggests the molecular underpinnings of pH-induced conformational changes

The Escherichia coli NhaA antiporter couples the transport of H⁺ and Na⁺ (or Li⁺) ions to maintain the proper pH range and Na⁺ concentration in cells. A crystal structure of NhaA, solved at pH 4, comprises 12 transmembrane helices (TMs), arranged in two domains, with a large cytoplasm-facing funnel and a smaller periplasm-facing funnel. NhaA undergoes conformational changes, e.g. after pH elevation to alkaline ranges, and we used two computational approaches to explore them. On the basis of pseudo-symmetric features of the crystal structure, we predicted the structural architecture of an alternate, periplasm-facing state. In contrast to the crystal structure, the model presents a closed cytoplasmic funnel, and a periplasmic funnel of greater volume. To examine the transporter functional direction of motion, we conducted elastic network analysis of the crystal structure and detected two main normal modes of motion. Notably, both analyses predicted similar trends of conformational changes, consisting of an overall rotational motion of the two domains around a putative symmetry axis at the funnel centers, perpendicular to the membrane plane. This motion, along with conformational changes within specific helices, resulted in closure at the cytoplasmic end and opening at the periplasmic end. Cross-linking experiments, performed between segments on opposite sides of the cytoplasmic funnel, revealed pH-dependent interactions consistent with the proposed conformational changes. We suggest that the model-structure and predicted motion represent alkaline pH-induced conformational changes, mediated by a cluster of evolutionarily conserved, titratable residues, at the cytoplasmic ends of TMs II, V, and IX.

Orientational preferences of neighboring helices can drive ER insertion of a marginally hydrophobic transmembrane helix

α-helical integral membrane proteins critically depend on the correct insertion of their transmembrane α helices into the lipid bilayer for proper folding, yet a surprisingly large fraction of the transmembrane α helices in multispanning integral membrane proteins are not sufficiently hydrophobic to insert into the target membrane by themselves. How can such marginally hydrophobic segments nevertheless form transmembrane helices in the folded structure? Here, we show that a transmembrane helix with a strong orientational preference (N(cyt)-C(lum) or N(lum)-C(cyt)) can both increase and decrease the hydrophobicity threshold for membrane insertion of a neighboring, marginally hydrophobic helix. This effect helps explain the "missing hydrophobicity" in polytopic membrane proteins.

Efficiency of detergents at maintaining membrane protein structures in their biologically relevant forms

High-resolution structural analysis of membrane proteins by X-ray crystallography or solution NMR spectroscopy often requires their solubilization in the membrane-mimetic environments of detergents. Yet the choice of a detergent suitable for a given study remains largely empirical. In the present work, we considered the micelle-crystallized structures of lactose permease (LacY), the sodium/galactose symporter (vSGLT), the vitamin B(12) transporter (BtuCD), and the arginine/agmatine antiporter (AdiC). Representative transmembrane (TM) segments were selected from these proteins based on their relative contact(s) with water, lipid, and/or within the protein, and were synthesized as Lys-tagged peptides. Each peptide was studied by circular dichroism and fluorescence spectroscopy in water, and in the presence of the detergents sodium dodecylsulfate (SDS, anionic); n-dodecyl phosphatidylcholine (DPC, zwitterionic); n-dodecyl-β-d-maltoside (DDM, neutral); and n-octyl-β-d-glucoside (OG, neutral, varying acyl tail length). We found that (i) the secondary structures of the TM segments were statistically indistinguishable in the four detergents studied; and (ii) a strong correlation exists between the extent of helical structure of each individual TM segment in detergents with its helicity level as it exists in the full-length protein, indicating that helix adoption is fundamentally the same in both environments. The denaturing properties of so-called 'harsh' detergents may thus largely be due to their interactions with non-membranous regions of proteins. Given the consistency of structural features observed for each TM segment in a variety of micellar media, the overall results suggest that the structure likely corresponds to its relevant biological form in the intact protein in its native lipid bilayer environment.