Structural imperatives impose diverse evolutionary constraints on helical membrane proteins

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.

C-terminal amino acids in the type I transmembrane domain of L-type lectin VIP36 affect γ-secretase susceptibility

Regulated intramembrane proteolysis (RIP) is a two-step processing mechanism for transmembrane proteins consisting of ectodomain shedding (shedding), which removes the extracellular domain through juxtamembrane processing and intramembrane proteolysis, which processes membrane-anchored shedding products within the transmembrane domain. RIP irreversibly converts one transmembrane protein into multiple soluble proteins that perform various physiological functions. The only requirement for the substrate of γ-secretase, the major enzyme responsible for intramembrane proteolysis of type I transmembrane proteins, is the absence of a large extracellular domain, and it is thought that γ-secretase can process any type I membrane protein as long as it is shed. In the present study, we showed that the shedding susceptible type I membrane protein VIP36 (36 kDa vesicular integral membrane protein) and its homolog, VIPL, have different γ-secretase susceptibilities in their transmembrane domains. Analysis of the substitution mutants suggested that γ-secretase susceptibility is regulated by C-terminal amino acids in the transmembrane domain. We also compared the transmembrane domains of several shedding susceptible membrane proteins and found that each had a different γ-secretase susceptibility. These results suggest that the transmembrane domain is not simply a stretch of hydrophobic amino acids but is an important element that regulates membrane protein function by controlling the lifetime of the membrane-anchored shedding product.

Lipid Bilayer Strengthens the Cooperative Network of a Membrane-Integral Enzyme

Lipid bilayer provides a two-dimensional hydrophobic solvent milieu for membrane proteins in cells. Although the native bilayer is widely recognized as an optimal environment for folding and function of membrane proteins, the underlying physical basis remains elusive. Here, employing the intramembrane protease GlpG of Escherichia coli as a model, we elucidate how the bilayer stabilizes a membrane protein and engages the protein's residue interaction network compared to the nonnative hydrophobic medium, micelles. We find that the bilayer enhances GlpG stability by promoting residue burial in the protein interior compared to micelles. Strikingly, while the cooperative residue interactions cluster into multiple distinct regions in micelles, the whole packed regions of the protein act as a single cooperative unit in the bilayer. Molecular dynamics (MD) simulation indicates that lipids less efficiently solvate GlpG than detergents. Thus, the bilayerinduced enhancement of stability and cooperativity likely stems from the dominant intraprotein interactions outcompeting the weak lipid solvation. Our findings reveal a foundational mechanism in the folding, function, and quality control of membrane proteins. The enhanced cooperativity benefits function facilitating propagation of local structural perturbation across the membrane. However, the same phenomenon can render the proteins' conformational integrity vulnerable to missense mutations causing conformational diseases1,2.

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.

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.

Transmembrane Helices Are an Over-Presented and Evolutionarily Conserved Source of Major Histocompatibility Complex Class I and II Epitopes

Cytolytic T cell responses are predicted to be biased towards membrane proteins. The peptide-binding grooves of most alleles of histocompatibility complex class I (MHC-I) are relatively hydrophobic, therefore peptide fragments derived from human transmembrane helices (TMHs) are predicted to be presented more often as would be expected based on their abundance in the proteome. However, the physiological reason of why membrane proteins might be over-presented is unclear. In this study, we show that the predicted over-presentation of TMH-derived peptides is general, as it is predicted for bacteria and viruses and for both MHC-I and MHC-II, and confirmed by re-analysis of epitope databases. Moreover, we show that TMHs are evolutionarily more conserved, because single nucleotide polymorphisms (SNPs) are present relatively less frequently in TMH-coding chromosomal regions compared to regions coding for extracellular and cytoplasmic protein regions. Thus, our findings suggest that both cytolytic and helper T cells are more tuned to respond to membrane proteins, because these are evolutionary more conserved. We speculate that TMHs are less prone to mutations that enable pathogens to evade T cell responses.

The conserved serine transporter SdaC moonlights to enable self recognition

Cells can use self recognition to achieve cooperative behaviors. Self-recognition genes are thought to principally evolve in tandem with partner self-recognition alleles. However, other constraints on protein evolution could exist. Here, we have identified an interaction outside of self-recognition loci that could constrain the sequence variation of a self-recognition protein. We show that during collective swarm expansion in Proteus mirabilis, self-recognition signaling co-opts SdaC, a serine transporter. Serine uptake is crucial for bacterial survival and colonization. Single-residue variants of SdaC reveal that self recognition requires an open conformation of the protein; serine transport is dispensable. A distant ortholog from Escherichia coli is sufficient for self recognition; however, a paralogous serine transporter, YhaO, is not. Thus, SdaC couples self recognition and serine transport, likely through a shared molecular interface. Self recognition proteins may follow the framework of a complex interaction network rather than an isolated two-protein system. Understanding molecular and ecological constraints on self-recognition proteins lays the groundwork for insights into the evolution of self recognition and emergent collective behaviors. Importance Bacteria can receive secret messages from kin during migration. For Proteus mirabilis, these messages are necessary for virulence in multi-species infections. We show that a serine transporter-conserved among gamma-enterobacteria- enables self recognition. Molecular co-option of nutrient uptake could limit the sequence variation of these message proteins. SdaC is the primary transporter for L-serine, a vital metabolite for colonization during disease. Unlike many self-recognition receptors, SdaC is sufficiently conserved between species to achieve recognition. The predicted open conformation is shared by transport and recognition. SdaC reveals the interdependence of communication and nutrient acquisition. As the broader interactions of self-recognition proteins are studied, features shared among microbial self-recognition systems, such as Dictyostelium spp. and Neurospora spp., could emerge.

Principal component analysis of alpha-helix deformations in transmembrane proteins

α-helices are deformable secondary structural components regularly observed in protein folds. The overall flexibility of an α-helix can be resolved into constituent physical deformations such as bending in two orthogonal planes and twisting along the principal axis. We used Principal Component Analysis to identify and quantify the contribution of each of these dominant deformation modes in transmembrane α-helices, extramembrane α-helices, and α-helices in soluble proteins. Using three α-helical samples from Protein Data Bank entries spanning these three cellular contexts, we determined that the relative contributions of these modes towards total deformation are independent of the α-helix's surroundings. This conclusion is supported by the observation that the identities of the top three deformation modes, the scaling behaviours of mode eigenvalues as a function of α-helix length, and the percentage contribution of individual modes on total variance were comparable across all three α-helical samples. These findings highlight that α-helical deformations are independent of cellular location and will prove to be valuable in furthering the development of flexible templates in de novo protein design.

Altering CLC stoichiometry by reducing non-polar side-chains at the dimerization interface

CLC-ec1 is a Cl^-/H⁺ antiporter that forms stable homodimers in lipid bilayers, with a free energy of -10.9 kcal/mol in 2:1 POPE/POPG lipid bilayers. The dimerization interface is formed by four transmembrane helices: H, I, P and Q, that are lined by non-polar side-chains that come in close contact, yet it is unclear as to whether their interactions drive dimerization. To investigate whether non-polar side-chains are required for dimer assembly, we designed a series of constructs where side-chain packing in the dimer state is significantly reduced by making 4-5 alanine substitutions along each helix (H-ala, I-ala, P-ala, Q-ala). All constructs are functional and three purify as stable dimers in detergent micelles despite the removal of significant side-chain interactions. On the other hand, H-ala shows the unique behavior of purifying as a mixture of monomers and dimers, followed by a rapid and complete conversion to monomers. In lipid bilayers, all four constructs are monomeric as examined by single-molecule photobleaching analysis. Further study of the H-helix shows that the single mutation L194A is sufficient to yield monomeric CLC-ec1 in detergent micelles and lipid bilayers. X-ray crystal structures of L194A reveal the protein re-assembles to form dimers, with a structure that is identical to wild-type. Altogether, these results demonstrate that non-polar membrane embedded side-chains play an important role in defining dimer stability, but the stoichiometry is highly contextual to the solvent environment. Furthermore, we discovered that L194 is a molecular hot-spot for defining dimerization of CLC-ec1.

Molecular evolution of a collage of cholesterol interaction motifs in transmembrane helix V of the serotonin1A receptor

The human serotonin_1A receptor is a representative member of the superfamily of G protein-coupled receptors (GPCRs) and an important drug target for neurological disorders. Using a combination of biochemical, biophysical and molecular dynamics simulation approaches, we and others have shown that membrane cholesterol modulates the organization, dynamics and function of vertebrate serotonin_1A receptors. Previous studies have shown that the cytoplasmic portion of transmembrane helix V (TM V) and the extramembraneous intracellular loop 3 are critical for G-protein coupling, phosphorylation and desensitization of the receptor. We have recently resolved a collage of putative cholesterol interaction motifs from the amino acid sequence overlapping this region. In this paper, we explore the sequence plasticity of this fragment that may have adapted to altered membrane lipidome, after vertebrates evolved from primordial invertebrates. Since invertebrates have lower levels of membrane cholesterol relative to vertebrates, we compared TM V sequence fragments from invertebrate serotonin₁ receptors with vertebrate orthologs to infer the sequence plasticity in TM V. We report that the average number of cholesterol interaction motifs in TM V for diverse phyla represents an increasing trend that could mirror vertebrate evolution from primordial invertebrates. By statistical modeling, we propose that the collage of cholesterol interaction motifs in TM V of the human serotonin_1A receptor may have evolved from rudimentary collages, reminiscent of primordial invertebrate orthologs. Taken together, we propose that a repertoire of cholesterol-philic nonsynonymous substitutions may have enhanced collage complexity in TM V during vertebrate evolution.

Structural cavities are critical to balancing stability and activity of a membrane-integral enzyme

Packing interaction is a critical driving force in the folding of helical membrane proteins. Despite the importance, packing defects (i.e., cavities including voids, pockets, and pores) are prevalent in membrane-integral enzymes, channels, transporters, and receptors, playing essential roles in function. Then, a question arises regarding how the two competing requirements, packing for stability vs. cavities for function, are reconciled in membrane protein structures. Here, using the intramembrane protease GlpG of Escherichia coli as a model and cavity-filling mutation as a probe, we tested the impacts of native cavities on the thermodynamic stability and function of a membrane protein. We find several stabilizing mutations which induce substantial activity reduction without distorting the active site. Notably, these mutations are all mapped onto the regions of conformational flexibility and functional importance, indicating that the cavities facilitate functional movement of GlpG while compromising the stability. Experiment and molecular dynamics simulation suggest that the stabilization is induced by the coupling between enhanced protein packing and weakly unfavorable lipid desolvation, or solely by favorable lipid solvation on the cavities. Our result suggests that, stabilized by the relatively weak interactions with lipids, cavities are accommodated in membrane proteins without severe energetic cost, which, in turn, serve as a platform to fine-tune the balance between stability and flexibility for optimal activity.

De novo protein design, a retrospective

Proteins are molecular machines whose function depends on their ability to achieve complex folds with precisely defined structural and dynamic properties. The rational design of proteins from first-principles, or de novo, was once considered to be impossible, but today proteins with a variety of folds and functions have been realized. We review the evolution of the field from its earliest days, placing particular emphasis on how this endeavor has illuminated our understanding of the principles underlying the folding and function of natural proteins, and is informing the design of macromolecules with unprecedented structures and properties. An initial set of milestones in de novo protein design focused on the construction of sequences that folded in water and membranes to adopt folded conformations. The first proteins were designed from first-principles using very simple physical models. As computers became more powerful, the use of the rotamer approximation allowed one to discover amino acid sequences that stabilize the desired fold. As the crystallographic database of protein structures expanded in subsequent years, it became possible to construct proteins by assembling short backbone fragments that frequently recur in Nature. The second set of milestones in de novo design involves the discovery of complex functions. Proteins have been designed to bind a variety of metals, porphyrins, and other cofactors. The design of proteins that catalyze hydrolysis and oxygen-dependent reactions has progressed significantly. However, de novo design of catalysts for energetically demanding reactions, or even proteins that bind with high affinity and specificity to highly functionalized complex polar molecules remains an importnant challenge that is now being achieved. Finally, the protein design contributed significantly to our understanding of membrane protein folding and transport of ions across membranes. The area of membrane protein design, or more generally of biomimetic polymers that function in mixed or non-aqueous environments, is now becoming increasingly possible.

Genetic intolerance analysis as a tool for protein science

Recent advances in whole genome and exome sequencing have dramatically increased the database of human gene variations. There are now enough sequenced human exomes and genomes to begin to identify gene variations that are notable because they are NOT observed in sequenced human genomes, apparently because they are subject to "purifying selection", exemplifying genetic intolerance. Such "dysprocreative" gene variations are embryonic lethal or prevent reproduction through any one of a number of possible mechanisms. Here we review an emerging quantitative approach, "Missense Tolerance Ratio" (MTR) analysis, that is used to assess protein-encoding gene (cDNA) sequence intolerance to missense mutations based on analysis of the >100 K and growing number of currently available human genome and exome sequences. This approach is already useful for analyzing intolerance to mutations in cDNA segments with a resolution on the order of 90 bases. Moreover, as the number of sequenced genomes/exomes increases by orders of magnitude it may eventually be possible to assess mutational tolerance in a statistically robust manner at or near single site resolution. Here we focus on how cDNA intolerance analysis complements other bioinformatic methods to illuminate structure-folding-function relationships for the encoded proteins. A set of disease-linked membrane proteins is employed to provide examples.

Packing of apolar side chains enables accurate design of highly stable membrane proteins

The features that stabilize the structures of membrane proteins remain poorly understood. Polar interactions contribute modestly, and the hydrophobic effect contributes little to the energetics of apolar side-chain packing in membranes. Disruption of steric packing can destabilize the native folds of membrane proteins, but is packing alone sufficient to drive folding in lipids? If so, then membrane proteins stabilized by this feature should be readily designed and structurally characterized-yet this has not been achieved. Through simulation of the natural protein phospholamban and redesign of variants, we define a steric packing code underlying its assembly. Synthetic membrane proteins designed using this code and stabilized entirely by apolar side chains conform to the intended fold. Although highly stable, the steric complementarity required for their folding is surprisingly stringent. Structural informatics shows that the designed packing motif recurs across the proteome, emphasizing a prominent role for precise apolar packing in membrane protein folding, stabilization, and evolution.

Prediction of disease-associated mutations in the transmembrane regions of proteins with known 3D structure

Being able to assess the phenotypic effects of mutations is a much required capability in precision medicine. However, most of the currently available structure-based methods actually predict stability changes caused by mutations rather than their pathogenic potential. There are also no dedicated methods to predict damaging mutations specifically in transmembrane proteins. In this study we developed and applied a machine-learning approach to discriminate between disease-associated and benign point mutations in the transmembrane regions of proteins with known 3D structure. The method, called BorodaTM (BOosted RegressiOn trees for Disease-Associated mutations in TransMembrane proteins), was trained on sequence-, structure-, and energy-derived descriptors. When compared with the state-of-the-art methods, BorodaTM is superior in classifying point mutations in transmembrane regions. Using BorodaTM we have conducted a large-scale mutation analysis in the transmembrane regions of human proteins with known 3D structures. For each protein we generated structural models for all point mutations by replacing each residue to 19 possible residue types. We classified ~1.8 millions point mutations as benign or diseased-associated and made all predictions available as a Web-server at https://www.iitm.ac.in/bioinfo/MutHTP/boroda.php.

Accurate computational design of multipass transmembrane proteins

The computational design of transmembrane proteins with more than one membrane-spanning region remains a major challenge. We report the design of transmembrane monomers, homodimers, trimers, and tetramers with 76 to 215 residue subunits containing two to four membrane-spanning regions and up to 860 total residues that adopt the target oligomerization state in detergent solution. The designed proteins localize to the plasma membrane in bacteria and in mammalian cells, and magnetic tweezer unfolding experiments in the membrane indicate that they are very stable. Crystal structures of the designed dimer and tetramer-a rocket-shaped structure with a wide cytoplasmic base that funnels into eight transmembrane helices-are very close to the design models. Our results pave the way for the design of multispan membrane proteins with new functions.

Design of self-assembling transmembrane helical bundles to elucidate principles required for membrane protein folding and ion transport

Ion transporters and channels are able to identify and act on specific substrates among myriads of ions and molecules critical to cellular processes, such as homeostasis, cell signalling, nutrient influx and drug efflux. Recently, we designed Rocker, a minimalist model for Zn²⁺/H⁺ co-transport. The success of this effort suggests that de novo membrane protein design has now come of age so as to serve a key approach towards probing the determinants of membrane protein folding, assembly and function. Here, we review general principles that can be used to design membrane proteins, with particular reference to helical assemblies with transport function. We also provide new functional and NMR data that probe the dynamic mechanism of conduction through Rocker.This article is part of the themed issue 'Membrane pores: from structure and assembly, to medicine and technology'.

Systematic Analyses and Prediction of Human Drug Side Effect Associated Proteins from the Perspective of Protein Evolution

Identification of various factors involved in adverse drug reactions in target proteins to develop therapeutic drugs with minimal/no side effect is very important. In this context, we have performed a comparative evolutionary rate analyses between the genes exhibiting drug side-effect(s) (SET) and genes showing no side effect (NSET) with an aim to increase the prediction accuracy of SET/NSET proteins using evolutionary rate determinants. We found that SET proteins are more conserved than the NSET proteins. The rates of evolution between SET and NSET protein primarily depend upon their noncomplex (protein complex association number = 0) forming nature, phylogenetic age, multifunctionality, membrane localization, and transmembrane helix content irrespective of their essentiality, total druggability (total number of drugs/target), m-RNA expression level, and tissue expression breadth. We also introduced two novel terms—killer druggability (number of drugs with killing side effect(s)/target), essential druggability (number of drugs targeting essential proteins/target) to explain the evolutionary rate variation between SET and NSET proteins. Interestingly, we noticed that SET proteins are younger than NSET proteins and multifunctional younger SET proteins are candidates of acquiring killing side effects. We provide evidence that higher killer druggability, multifunctionality, and transmembrane helices support the conservation of SET proteins over NSET proteins in spite of their recent origin. By employing all these entities, our Support Vector Machine model predicts human SET/NSET proteins to a high degree of accuracy (∼86%).

Reversible Unfolding of Rhomboid Intramembrane Proteases

Denaturant-induced unfolding of helical membrane proteins provides insights into their mechanism of folding and domain organization, which take place in the chemically heterogeneous, anisotropic environment of a lipid membrane. Rhomboid proteases are intramembrane proteases that play key roles in various diseases. Crystal structures have revealed a compact helical bundle with a buried active site, which requires conformational changes for the cleavage of transmembrane substrates. A dimeric form of the rhomboid protease has been shown to be important for activity. In this study, we examine the mechanism of refolding for two distinct rhomboids to gain insight into their secondary structure-activity relationships. Although helicity is largely abolished in the unfolded states of both proteins, unfolding is completely reversible for HiGlpG but only partially reversible for PsAarA. Refolding of both proteins results in reassociation of the dimer, with a 90% regain of catalytic activity for HiGlpG but only a 70% regain for PsAarA. For both proteins, a broad, gradual transition from the native, folded state to the denatured, partly unfolded state was revealed with the aid of circular dichroism spectroscopy as a function of denaturant concentration, thus arguing against a classical two-state model as found for many globular soluble proteins. Thermal denaturation has irreversible destabilizing effects on both proteins, yet reveals important functional details regarding substrate accessibility to the buried active site. This concerted biophysical and functional analysis demonstrates that HiGlpG, with a simple six-transmembrane-segment organization, is more robust than PsAarA, which has seven predicted transmembrane segments, thus rendering HiGlpG amenable to in vitro studies of membrane-protein folding.

Membrane Proteins Are Dramatically Less Conserved than Water-Soluble Proteins across the Tree of Life

Membrane proteins are crucial in transport, signaling, bioenergetics, catalysis, and as drug targets. Here, we show that membrane proteins have dramatically fewer detectable orthologs than water-soluble proteins, less than half in most species analyzed. This sparse distribution could reflect rapid divergence or gene loss. We find that both mechanisms operate. First, membrane proteins evolve faster than water-soluble proteins, particularly in their exterior-facing portions. Second, we demonstrate that predicted ancestral membrane proteins are preferentially lost compared with water-soluble proteins in closely related species of archaea and bacteria. These patterns are consistent across the whole tree of life, and in each of the three domains of archaea, bacteria, and eukaryotes. Our findings point to a fundamental evolutionary principle: membrane proteins evolve faster due to stronger adaptive selection in changing environments, whereas cytosolic proteins are under more stringent purifying selection in the homeostatic interior of the cell. This effect should be strongest in prokaryotes, weaker in unicellular eukaryotes (with intracellular membranes), and weakest in multicellular eukaryotes (with extracellular homeostasis). We demonstrate that this is indeed the case. Similarly, we show that extracellular water-soluble proteins exhibit an even stronger pattern of low homology than membrane proteins. These striking differences in conservation of membrane proteins versus water-soluble proteins have important implications for evolution and medicine.

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 membrane- and soluble-protein helix-helix interactome: similar geometry via different interactions

α Helices are a basic unit of protein secondary structure and therefore the interaction between helices is crucial to understanding tertiary and higher-order folds. Comparing subtle variations in the structural and sequence motifs between membrane and soluble proteins sheds light on the different constraints faced by each environment and elucidates the complex puzzle of membrane protein folding. Here, we demonstrate that membrane and water-soluble helix pairs share a small number of similar folds with various interhelical distances. The composition of the residues that pack at the interface between corresponding motifs shows that hydrophobic residues tend to be more enriched in the water-soluble class of structures and small residues in the transmembrane class. The latter group facilitates packing via sidechain- and backbone-mediated hydrogen bonds within the low-dielectric membrane milieu. The helix-helix interactome space, with its associated sequence preferences and accompanying hydrogen-bonding patterns, should be useful for engineering, prediction, and design of protein structure.

Comprehensive, structurally-informed alignment and phylogeny of vertebrate biogenic amine receptors

Biogenic amine receptors play critical roles in regulating behavior and physiology in both vertebrates and invertebrates, particularly within the central nervous system. Members of the G-protein coupled receptor (GPCR) family, these receptors interact with endogenous bioamine ligands such as dopamine, serotonin, and epinephrine, and are targeted by a wide array of pharmaceuticals. Despite the clear clinical and biological importance of these receptors, their evolutionary history remains poorly characterized. In particular, the relationships among biogenic amine receptors and any specific evolutionary constraints acting within distinct receptor subtypes are largely unknown. To advance and facilitate studies in this receptor family, we have constructed a comprehensive, high-quality sequence alignment of vertebrate biogenic amine receptors. In particular, we have integrated a traditional multiple sequence approach with robust structural domain predictions to ensure that alignment columns accurately capture the highly-conserved GPCR structural domains, and we demonstrate how ignoring structural information produces spurious inferences of homology. Using this alignment, we have constructed a structurally-partitioned maximum-likelihood phylogeny from which we deduce novel biogenic amine receptor relationships and uncover previously unrecognized lineage-specific receptor clades. Moreover, we find that roughly 1% of the 3039 sequences in our final alignment are either misannotated or unclassified, and we propose updated classifications for these receptors. We release our comprehensive alignment and its corresponding phylogeny as a resource for future research into the evolution and diversification of biogenic amine receptors.

Elucidating the genotype-phenotype relationships and network perturbations of human shared and specific disease genes from an evolutionary perspective

To date, numerous studies have been attempted to determine the extent of variation in evolutionary rates between human disease and nondisease (ND) genes. In our present study, we have considered human autosomal monogenic (Mendelian) disease genes, which were classified into two groups according to the number of phenotypic defects, that is, specific disease (SPD) gene (one gene: one defect) and shared disease (SHD) gene (one gene: multiple defects). Here, we have compared the evolutionary rates of these two groups of genes, that is, SPD genes and SHD genes with respect to ND genes. We observed that the average evolutionary rates are slow in SHD group, intermediate in SPD group, and fast in ND group. Group-to-group evolutionary rate differences remain statistically significant regardless of their gene expression levels and number of defects. We demonstrated that disease genes are under strong selective constraint if they emerge through edgetic perturbation or drug-induced perturbation of the interactome network, show tissue-restricted expression, and are involved in transmembrane transport. Among all the factors, our regression analyses interestingly suggest the independent effects of 1) drug-induced perturbation and 2) the interaction term of expression breadth and transmembrane transport on protein evolutionary rates. We reasoned that the drug-induced network disruption is a combination of several edgetic perturbations and, thus, has more severe effect on gene phenotypes.

Membrane proteins can have high kinetic stability

Approximately 10% of water-soluble proteins are considered kinetically stable with unfolding half-lives in the range of weeks to millenia. These proteins only rarely sample the unfolded state and may never unfold on their respective biological time scales. It is still not known whether membrane proteins can be kinetically stable, however. Here we examine the subunit dissociation rate of the trimeric membrane enzyme, diacylglycerol kinase, from Escherichia coli as a proxy for complete unfolding. We find that dissociation occurs with a half-life of at least several weeks, demonstrating that membrane proteins can remain locked in a folded state for long periods of time. These results reveal that evolution can use kinetic stability to regulate the biological function of membrane proteins, as it can for soluble proteins. Moreover, it appears that the generation of kinetic stability could be a viable target for membrane protein engineering efforts.

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.

Ca2+ signalling early in evolution--all but primitive

Early in evolution, Ca(2+) emerged as the most important second messenger for regulating widely different cellular functions. In eukaryotic cells Ca(2+) signals originate from several sources, i.e. influx from the outside medium, release from internal stores or from both. In mammalian cells, Ca(2+)-release channels represented by inositol 1,4,5-trisphosphate receptors and ryanodine receptors (InsP3R and RyR, respectively) are the most important. In unicellular organisms and plants, these channels are characterised with much less precision. In the ciliated protozoan, Paramecium tetraurelia, 34 molecularly distinct Ca(2+)-release channels that can be grouped in six subfamilies, based on criteria such as domain structure, pore, selectivity filter and activation mechanism have been identified. Some of these channels are genuine InsP3Rs and some are related to RyRs. Others show some--but not all--features that are characteristic for one or the other type of release channel. Localisation and gene silencing experiments revealed widely different--yet distinct--localisation, activation and functional engagement of the different Ca(2+)-release channels. Here, we shall discuss early evolutionary routes of Ca(2+)-release machinery in protozoa and demonstrate that detailed domain analyses and scrutinised functional analyses are instrumental for in-depth evolutionary mapping of Ca(2+)-release channels in unicellular organisms.

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.

Membrane environment imposes unique selection pressures on transmembrane domains of G protein-coupled receptors

We have investigated the influence of the plasma membrane environment on the molecular evolution of G protein-coupled receptors (GPCRs), the largest receptor family in Metazoa. In particular, we have analyzed the site-specific rate variation across the two primary structural partitions, transmembrane (TM) and extramembrane (EM), of these membrane proteins. We find that TM domains evolve more slowly than do EM domains, though TM domains display increased rate heterogeneity relative to their EM counterparts. Although the majority of residues across GPCRs experience strong to weak purifying selection, many GPCRs experience positive selection at both TM and EM residues, albeit with a slight bias towards the EM. Further, a subset of GPCRs, chemosensory receptors (including olfactory and taste receptors), exhibit increased rates of evolution relative to other GPCRs, an effect which is more pronounced in their TM spans. Although it has been previously suggested that the TM's low evolutionary rate is caused by their high percentage of buried residues, we show that their attenuated rate seems to stem from the strong biophysical constraints of the membrane itself, or by functional requirements. In spite of the strong evolutionary constraints acting on the TM spans of GPCRs, positive selection and high levels of evolutionary rate variability are common. Thus, biophysical constraints should not be presumed to preclude a protein's ability to evolve.

Sequence alignment of viral channel proteins with cellular ion channels

Sequence alignment is an important tool for identifying regions of similarities among proteins and for, thus, establishing functional and structural relationships between different proteins. Here, alignments of transmembrane domains (TMDs) of viral channel forming proteins with host ion channels and toxins are evaluated. The following representatives of polytopic viral channel proteins are chosen: (i) p7 of HCV and 2B of Polio virus (two TMDs) and (ii) 3a of SARS-CoV (three TMDs). Using ClustalW2, each of the TMDs of the viral channels is aligned, and the overlap is mapped onto structural models of the host channels and toxins focusing on the pore-lining TMDs. The analysis reveals that p7 and 2B TMDs align with the pore-facing TMD of MscL, and 3a-TMDs align with those of ligand-gated ion channels. Possible implications concerning the mechanism of function of the viral proteins are discussed.

Architectural and thermodynamic principles underlying intramembrane protease function

Intramembrane proteases hydrolyze peptide bonds within the membrane as a signaling paradigm universal to all life forms and with implications in disease. Deciphering the architectural strategies supporting intramembrane proteolysis is an essential but unattained goal. We integrated a new, quantitative and high-throughput thermal light-scattering technology, reversible equilibrium un/refolding, and quantitative protease assays to interrogate rhomboid architecture with 151 purified variants. Rhomboid proteases maintain low intrinsic thermodynamic stability (ΔG=2.1-4.5kcal/mol) resulting from a multitude of generally-weak transmembrane packing interactions, making them highly-responsive to their environment. Stability is consolidated by two buried glycines and several packing leucines, with a few multifaceted hydrogen bonds strategically-deployed to two peripheral regions. Opposite these regions lie transmembrane segment 5 and connected loops that are notably exempt of structural responsibility, suggesting intramembrane proteolysis involves considerable but localized protein dynamics. Our analyses provide a comprehensive ‘heat map’ of the physio-chemical anatomy underlying membrane-immersed enzyme function at unprecedented resolution.

EvoluCode: Evolutionary Barcodes as a Unifying Framework for Multilevel Evolutionary Data

Evolutionary systems biology aims to uncover the general trends and principles governing the evolution of biological networks. An essential part of this process is the reconstruction and analysis of the evolutionary histories of these complex, dynamic networks. Unfortunately, the methodologies for representing and exploiting such complex evolutionary histories in large scale studies are currently limited. Here, we propose a new formalism, called EvoluCode (Evolutionary barCode), which allows the integration of different evolutionary parameters (eg, sequence conservation, orthology, synteny …) in a unifying format and facilitates the multilevel analysis and visualization of complex evolutionary histories at the genome scale. The advantages of the approach are demonstrated by constructing barcodes representing the evolution of the complete human proteome. Two large-scale studies are then described: (i) the mapping and visualization of the barcodes on the human chromosomes and (ii) automatic clustering of the barcodes to highlight protein subsets sharing similar evolutionary histories and their functional analysis. The methodologies developed here open the way to the efficient application of other data mining and knowledge extraction techniques in evolutionary systems biology studies. A database containing all EvoluCode data is available at: http://lbgi.igbmc.fr/barcodes.

Calcium-release channels in paramecium. Genomic expansion, differential positioning and partial transcriptional elimination

The release of Ca²⁺ from internal stores is a major source of signal Ca²⁺ in almost all cell types. The internal Ca²⁺ pools are activated via two main families of intracellular Ca²⁺-release channels, the ryanodine and the inositol 1,4,5-trisphosphate (InsP₃) receptors. Among multicellular organisms these channel types are ubiquitous, whereas in most unicellular eukaryotes the identification of orthologs is impaired probably due to evolutionary sequence divergence. However, the ciliated protozoan Paramecium allowed us to prognosticate six groups, with a total of 34 genes, encoding proteins with characteristics typical of InsP₃ and ryanodine receptors by BLAST search of the Paramecium database. We here report that these Ca²⁺-release channels may display all or only some of the characteristics of canonical InsP₃ and ryanodine receptors. In all cases, prediction methods indicate the presence of six trans-membrane regions in the C-terminal domains, thus corresponding to canonical InsP₃ receptors, while a sequence homologous to the InsP₃-binding domain is present only in some types. Only two types have been analyzed in detail previously. We now show, by using antibodies and eventually by green fluorescent protein labeling, that the members of all six groups localize to distinct organelles known to participate in vesicle trafficking and, thus, may provide Ca²⁺ for local membrane-membrane interactions. Whole genome duplication can explain radiation within the six groups. Comparative and evolutionary evaluation suggests derivation from a common ancestor of canonical InsP₃ and ryanodine receptors. With one group we could ascertain, to our knowledge for the first time, aberrant splicing in one thoroughly analyzed Paramecium gene. This yields truncated forms and, thus, may indicate a way to pseudogene formation. No comparable analysis is available for any other, free-living or parasitic/pathogenic protozoan.

Predicting residue-residue contacts and helix-helix interactions in transmembrane proteins using an integrative feature-based random forest approach

Integral membrane proteins constitute 25-30% of genomes and play crucial roles in many biological processes. However, less than 1% of membrane protein structures are in the Protein Data Bank. In this context, it is important to develop reliable computational methods for predicting the structures of membrane proteins. Here, we present the first application of random forest (RF) for residue-residue contact prediction in transmembrane proteins, which we term as TMhhcp. Rigorous cross-validation tests indicate that the built RF models provide a more favorable prediction performance compared with two state-of-the-art methods, i.e., TMHcon and MEMPACK. Using a strict leave-one-protein-out jackknifing procedure, they were capable of reaching the top L/5 prediction accuracies of 49.5% and 48.8% for two different residue contact definitions, respectively. The predicted residue contacts were further employed to predict interacting helical pairs and achieved the Matthew's correlation coefficients of 0.430 and 0.424, according to two different residue contact definitions, respectively. To facilitate the academic community, the TMhhcp server has been made freely accessible at http://protein.cau.edu.cn/tmhhcp.

The rhomboid protease family: a decade of progress on function and mechanism

Summary

Rhomboid proteases are the largest family of enzymes that hydrolyze peptide bonds within the cell membrane. Although discovered to be serine proteases only a decade ago, rhomboid proteases are already considered to be the best understood intramembrane proteases. The presence of rhomboid proteins in all domains of life emphasizes their importance but makes their evolutionary history difficult to chart with confidence. Phylogenetics nevertheless offers three guiding principles for interpreting rhomboid function. The near ubiquity of rhomboid proteases across evolution suggests broad, organizational roles that are not directly essential for cell survival. Functions have been deciphered in only about a dozen organisms and fall into four general categories: initiating cell signaling in animals, facilitating bacterial quorum sensing, regulating mitochondrial homeostasis, and dismantling adhesion complexes of parasitic protozoa. Although in no organism has the full complement of rhomboid function yet been elucidated, links to devastating human disease are emerging rapidly, including to Parkinson's disease, type II diabetes, cancer, and bacterial and malaria infection. Rhomboid proteases are unlike most proteolytic enzymes, because they are membrane-immersed; understanding how the membrane immersion affects their function remains a key challenge.

Tolerance to changes in membrane lipid composition as a selected trait of membrane proteins

Membrane lipid composition can vary dramatically across the three domains of life and even within single organisms. Here we review evidence that the lipid-exposed surfaces of membrane proteins have generally evolved to maintain correct structure and function in the face of major changes in lipid composition. Such tolerance has allowed evolution to extensively remodel membrane lipid compositions during the emergence of new species without having to extensively remodel the associated membrane proteins. The tolerance of membrane proteins also permits single-cell organisms to vary their membrane lipid composition in response to their changing environments and allows dynamic and organelle-specific variations in the lipid compositions of eukaryotic cells. Membrane protein structural biology has greatly benefited from this seemingly intrinsic property of membrane proteins: the majority of structures determined to date have been characterized under model membrane conditions that little resemble those of native membranes. Nevertheless, with a few notable exceptions, most experimentally determined membrane protein structures appear, to a good approximation, to faithfully report on native structure.

Stable folding core in the folding transition state of an alpha-helical integral membrane protein

Defining the structural features of a transition state is important in understanding a folding reaction. Here, we use Φ-value and double mutant analyses to probe the folding transition state of the membrane protein bacteriorhodopsin. We focus on the final C-terminal helix, helix G, of this seven transmembrane helical protein. Φ-values could be derived for 12 amino acid residues in helix G, most of which have low or intermediate values, suggesting that native structure is disrupted at these amino acid positions in the transition state. Notably, a cluster of residues between E204 and M209 all have Φ-values close to zero. Disruption of helix G is further confirmed by a low Φ-value of 0.2 between residues T170 on helix F and S226 on helix G, suggesting the absence of a native hydrogen bond between helices F and G. Φ-values for paired mutations involved in four interhelical hydrogen bonds revealed that all but one of these bonds is absent in the transition state. The unstructured helix G contrasts with Φ-values along helix B that are generally high, implying native structure in helix B in the transition state. Thus helix B seems to constitute part of a stable folding nucleus while the consolidation of helix G is a relatively late folding event. Polarization of secondary structure correlates with sequence position, with a structured helix B near the N terminus contrasting with an unstructured C-terminal helix G.

Viral channel forming proteins - modeling the target

The cellular and subcellular membranes encounter an important playground for the activity of membrane proteins encoded by viruses. Viral membrane proteins, similar to their host companions, can be integral or attached to the membrane. They are involved in directing the cellular and viral reproduction, the fusion and budding processes. This review focuses especially on those integral viral membrane proteins which form channels or pores, the classification to be so, modeling by in silico methods and potential drug candidates. The sequence of an isolate of Vpu from HIV-1 is aligned with host ion channels and a toxin. The focus is on the alignment of the transmembrane domains. The results of the alignment are mapped onto the 3D structures of the respective channels and toxin. The results of the mapping support the idea of a 'channel-pore dualism' for Vpu.