Helical packing patterns in membrane and soluble proteins

Scaling up Functional Analyses of the G Protein-Coupled Receptor Rhodopsin

Eukaryotic cells use G protein-coupled receptors (GPCRs) to convert external stimuli into internal signals to elicit cellular responses. However, how mutations in GPCR-coding genes affect GPCR activation and downstream signaling pathways remain poorly understood. Approaches such as deep mutational scanning show promise in investigations of GPCRs, but a high-throughput method to measure rhodopsin activation has yet to be achieved. Here, we scale up a fluorescent reporter assay in budding yeast that we engineered to study rhodopsin's light-activated signal transduction. Using this approach, we measured the mutational effects of over 1200 individual human rhodopsin mutants, generated by low-frequency random mutagenesis of the GPCR rhodopsin (RHO) gene. Analysis of the data in the context of rhodopsin's three-dimensional structure reveals that transmembrane helices are generally less tolerant to mutations compared to flanking helices that face the lipid bilayer, which suggest that mutational tolerance is contingent on both the local environment surrounding specific residues and the specific position of these residues in the protein structure. Comparison of functional scores from our screen to clinically identified rhodopsin disease variants found many pathogenic mutants to be loss of function. Lastly, functional scores from our assay were consistent with a complex counterion mechanism involved in ligand-binding and rhodopsin activation. Our results demonstrate that deep mutational scanning is possible for rhodopsin activation and can be an effective method for revealing properties of mutational tolerance that may be generalizable to other transmembrane proteins.

Infrared Multiphoton Dissociation Enables Top-Down Characterization of Membrane Protein Complexes and G Protein-Coupled Receptors

Membrane proteins are challenging to analyze by native mass spectrometry (MS) as their hydrophobic nature typically requires stabilization in detergent micelles that are removed prior to analysis via collisional activation. There is however a practical limit to the amount of energy which can be applied, which often precludes subsequent characterization by top-down MS. To overcome this barrier, we have applied a modified Orbitrap Eclipse Tribrid mass spectrometer coupled to an infrared laser within a high-pressure linear ion trap. We show how tuning the intensity and time of incident photons enables liberation of membrane proteins from detergent micelles. Specifically, we relate the ease of micelle removal to the infrared absorption of detergents in both condensed and gas phases. Top-down MS via infrared multiphoton dissociation (IRMPD), results in good sequence coverage enabling unambiguous identification of membrane proteins and their complexes. By contrasting and comparing the fragmentation patterns of the ammonia channel with two class A GPCRs, we identify successive cleavage of adjacent amino acids within transmembrane domains. Using gas-phase molecular dynamics simulations, we show that areas prone to fragmentation maintain aspects of protein structure at increasing temperatures. Altogether, we propose a rationale to explain why and where in the protein fragment ions are generated.

Infrared Multiphoton Dissociation Enables Top-Down Characterization of Membrane Protein Complexes and G Protein-Coupled Receptors

Membrane proteins are challenging to analyze by native mass spectrometry (MS) as their hydrophobic nature typically requires stabilization in detergent micelles that are removed prior to analysis via collisional activation. There is however a practical limit to the amount of energy which can be applied, which often precludes subsequent characterization by top-down MS. To overcome this barrier, we have applied a modified Orbitrap Eclipse Tribrid mass spectrometer coupled to an infrared laser within a high-pressure linear ion trap. We show how tuning the intensity and time of incident photons enables liberation of membrane proteins from detergent micelles. Specifically, we relate the ease of micelle removal to the infrared absorption of detergents in both condensed and gas phases. Top-down MS via infrared multiphoton dissociation (IRMPD), results in good sequence coverage enabling unambiguous identification of membrane proteins and their complexes. By contrasting and comparing the fragmentation patterns of the ammonia channel with two class A GPCRs, we identify successive cleavage of adjacent amino acids within transmembrane domains. Using gas-phase molecular dynamics simulations, we show that areas prone to fragmentation maintain aspects of protein structure at increasing temperatures. Altogether, we propose a rationale to explain why and where in the protein fragment ions are generated.

3D interaction homology: The hydrophobic residues alanine, isoleucine, leucine, proline and valine play different structural roles in soluble and membrane proteins

The aliphatic hydrophobic amino acid residues—alanine, isoleucine, leucine, proline and valine—are among the most common found in proteins. Their structural role in proteins is seemingly obvious: engage in hydrophobic interactions to stabilize secondary, and to a lesser extent, tertiary and quaternary structure. However, favorable hydrophobic interactions involving the sidechains of these residue types are generally less significant than the unfavorable set arising from interactions with polar atoms. Importantly, the constellation of interactions between residue sidechains and their environments can be recorded as three-dimensional maps that, in turn, can be clustered. The clustered average map sets compose a library of interaction profiles encoding interaction strengths, interaction types and the optimal 3D position for the interacting partners. This library is backbone angle-dependent and suggests solvent and lipid accessibility for each unique interaction profile. In this work, in addition to analysis of soluble proteins, a large set of membrane proteins that contained optimized artificial lipids were evaluated by parsing the structures into three distinct components: soluble extramembrane domain, lipid facing transmembrane domain, core transmembrane domain. The aliphatic residues were extracted from each of these sets and passed through our calculation protocol. Notable observations include: the roles of aliphatic residues in soluble proteins and in the membrane protein’s soluble domains are nearly identical, although the latter are slightly more solvent accessible; by comparing maps calculated with sidechain-lipid interactions to maps ignoring those interactions, the potential extent of residue-lipid and residue-interactions can be assessed and likely exploited in structure prediction and modeling; amongst these residue types, the levels of lipid engagement show isoleucine as the most engaged, while the other residues are largely interacting with neighboring helical residues.

Resemblance-Ranking Peptide Library to Screen for Binders to Antibodies on a Peptidomic Scale

A novel resemblance-ranking peptide library with 160,000 10-meric peptides was designed to search for selective binders to antibodies. The resemblance-ranking principle enabled the selection of sequences that are most similar to the human peptidome. The library was synthesized with ultra-high-density peptide arrays. As proof of principle, screens for selective binders were performed for the therapeutic anti-CD20 antibody rituximab. Several features in the amino acid composition of antibody-binding peptides were identified. The selective affinity of rituximab increased with an increase in the number of hydrophobic amino acids in a peptide, mainly tryptophan and phenylalanine, while a total charge of the peptide remained relatively small. Peptides with a higher affinity exhibited a lower sum helix propensity. For the 30 strongest peptide binders, a substitutional analysis was performed to determine dissociation constants and the invariant amino acids for binding to rituximab. The strongest selective peptides had a dissociation constant in the hundreds of the nano-molar range. The substitutional analysis revealed a specific hydrophobic epitope for rituximab. To show that conformational binders can, in principle, be detected in array format, cyclic peptide substitutions that are similar to the target of rituximab were investigated. Since the specific binders selected via the resemblance-ranking peptide library were based on the hydrophobic interactions that are widespread in the world of biomolecules, the library can be used to screen for potential linear epitopes that may provide information about the cross-reactivity of antibodies.

A unified evolutionary origin for the ubiquitous protein transporters SecY and YidC

BACKGROUND

Protein transporters translocate hydrophilic segments of polypeptide across hydrophobic cell membranes. Two protein transporters are ubiquitous and date back to the last universal common ancestor: SecY and YidC. SecY consists of two pseudosymmetric halves, which together form a membrane-spanning protein-conducting channel. YidC is an asymmetric molecule with a protein-conducting hydrophilic groove that partially spans the membrane. Although both transporters mediate insertion of membrane proteins with short translocated domains, only SecY transports secretory proteins and membrane proteins with long translocated domains. The evolutionary origins of these ancient and essential transporters are not known.

RESULTS

The features conserved by the two halves of SecY indicate that their common ancestor was an antiparallel homodimeric channel. Structural searches with SecY's halves detect exceptional similarity with YidC homologs. The SecY halves and YidC share a fold comprising a three-helix bundle interrupted by a helical hairpin. In YidC, this hairpin is cytoplasmic and facilitates substrate delivery, whereas in SecY, it is transmembrane and forms the substrate-binding lateral gate helices. In both transporters, the three-helix bundle forms a protein-conducting hydrophilic groove delimited by a conserved hydrophobic residue. Based on these similarities, we propose that SecY originated as a YidC homolog which formed a channel by juxtaposing two hydrophilic grooves in an antiparallel homodimer. We find that archaeal YidC and its eukaryotic descendants use this same dimerisation interface to heterodimerise with a conserved partner. YidC's sufficiency for the function of simple cells is suggested by the results of reductive evolution in mitochondria and plastids, which tend to retain SecY only if they require translocation of large hydrophilic domains.

CONCLUSIONS

SecY and YidC share previously unrecognised similarities in sequence, structure, mechanism, and function. Our delineation of a detailed correspondence between these two essential and ancient transporters enables a deeper mechanistic understanding of how each functions. Furthermore, key differences between them help explain how SecY performs its distinctive function in the recognition and translocation of secretory proteins. The unified theory presented here explains the evolution of these features, and thus reconstructs a key step in the origin of cells.

Saturation Mutagenesis of the Transmembrane Region of HokC in Escherichia coli Reveals Its High Tolerance to Mutations

Cells adapt to different stress conditions, such as the antibiotics presence. This adaptation sometimes is achieved by changing relevant protein positions, of which the mutability is limited by structural constrains. Understanding the basis of these constrains represent an important challenge for both basic science and potential biotechnological applications. To study these constraints, we performed a systematic saturation mutagenesis of the transmembrane region of HokC, a toxin used by Escherichia coli to control its own population, and observed that 92% of single-point mutations are tolerated and that all the non-tolerated mutations have compensatory mutations that reverse their effect. We provide experimental evidence that HokC accumulates multiple compensatory mutations that are found as correlated mutations in the HokC family multiple sequence alignment. In agreement with these observations, transmembrane proteins show higher probability to present correlated mutations and are less densely packed locally than globular proteins; previous mutagenesis results on transmembrane proteins further support our observations on the high tolerability to mutations of transmembrane regions of proteins. Thus, our experimental results reveal the HokC transmembrane region high tolerance to loss-of-function mutations that is associated with low sequence conservation and high rate of correlated mutations in the HokC family sequences alignment, which are features shared with other transmembrane proteins.

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.

Principles and Methods in Computational Membrane Protein Design

After decades of progress in computational protein design, the design of proteins folding and functioning in lipid membranes appears today as the next frontier. Some notable successes in the de novo design of simplified model membrane protein systems have helped articulate fundamental principles of protein folding, architecture and interaction in the hydrophobic lipid environment. These principles are reviewed here, together with the computational methods and approaches that were used to identify them. We provide an overview of the methodological innovations in the generation of new protein structures and functions and in the development of membrane-specific energy functions. We highlight the opportunities offered by new machine learning approaches applied to protein design, and by new experimental characterization techniques applied to membrane proteins. Although membrane protein design is in its infancy, it appears more reachable than previously thought.

Membrane Barrels Are Taller, Fatter, Inside-Out Soluble Barrels

Up-and-down β-barrel topology exists in both the membrane and soluble environment. By comparing features of these structurally similar proteins, we can determine what features are particular to the environment rather than the fold. Here we compare structures of membrane β-barrels to soluble β-barrels and evaluate their relative size, shape, amino acid composition, hydrophobicity, and periodicity. We find that membrane β-barrels are generally larger than soluble β-barrels, with more strands per barrel and more amino acids per strand, making them wider and taller. We also find that membrane β-barrels are inside-out soluble β-barrels. The inward region of membrane β-barrels has similar hydrophobicity to the outward region of soluble β-barrels, and the outward region of membrane β-barrels has similar hydrophobicity to the inward region of the soluble β-barrels. Moreover, even though both types of β-barrel have been assumed to have strands with amino acids that alternate in direction and hydrophobicity, we find that the membrane β-barrels have more regular alternation than soluble β-barrels. These features give insight into how membrane barrels maintain their fold and function in the membrane.

Targeting Nsp9 as an anti-SARS-CoV-2 strategy

Non-structural protein 9 (Nsp9) plays a key role in viral replication of coronavirus and represents a promising target for anti-SARS-CoV-2 strategies.

Mutations in transmembrane proteins: diseases, evolutionary insights, prediction and comparison with globular proteins

Membrane proteins are unique in that they interact with lipid bilayers, making them indispensable for transporting molecules and relaying signals between and across cells. Due to the significance of the protein's functions, mutations often have profound effects on the fitness of the host. This is apparent both from experimental studies, which implicated numerous missense variants in diseases, as well as from evolutionary signals that allow elucidating the physicochemical constraints that intermembrane and aqueous environments bring. In this review, we report on the current state of knowledge acquired on missense variants (referred to as to single amino acid variants) affecting membrane proteins as well as the insights that can be extrapolated from data already available. This includes an overview of the annotations for membrane protein variants that have been collated within databases dedicated to the topic, bioinformatics approaches that leverage evolutionary information in order to shed light on previously uncharacterized membrane protein structures or interaction interfaces, tools for predicting the effects of mutations tailored specifically towards the characteristics of membrane proteins as well as two clinically relevant case studies explaining the implications of mutated membrane proteins in cancer and cardiomyopathy.

Membrane receptor activation mechanisms and transmembrane peptide tools to elucidate them

Single-pass membrane receptors contain extracellular domains that respond to external stimuli and transmit information to intracellular domains through a single transmembrane (TM) α-helix. Because membrane receptors have various roles in homeostasis, signaling malfunctions of these receptors can cause disease. Despite their importance, there is still much to be understood mechanistically about how single-pass receptors are activated. In general, single-pass receptors respond to extracellular stimuli via alterations in their oligomeric state. The details of this process are still the focus of intense study, and several lines of evidence indicate that the TM domain (TMD) of the receptor plays a central role. We discuss three major mechanistic hypotheses for receptor activation: ligand-induced dimerization, ligand-induced rotation, and receptor clustering. Recent observations suggest that receptors can use a combination of these activation mechanisms and that technical limitations can bias interpretation. Short peptides derived from receptor TMDs, which can be identified by screening or rationally developed on the basis of the structure or sequence of their targets, have provided critical insights into receptor function. Here, we explore recent evidence that, depending on the target receptor, TMD peptides cannot only inhibit but also activate target receptors and can accommodate novel, bifunctional designs. Furthermore, we call for more sharing of negative results to inform the TMD peptide field, which is rapidly transforming into a suite of unique tools with the potential for future therapeutics.

Design of a Helical-Stabilized, Cyclic, and Nontoxic Analogue of the Peptide Cm-p5 with Improved Antifungal Activity

Following the information obtained by a rational design study, a cyclic and helical-stabilized analogue of the peptide Cm-p5 was synthetized. The cyclic monomer showed an increased activity in vitro against Candida albicans and Candida parapsilosis, compared to Cm-p5. Initially, 14 mutants of Cm-p5 were synthesized following a rational design to improve the antifungal activity and pharmacological properties. Antimicrobial testing showed that the activity was lost in each of these 14 analogues, suggesting, as a main conclusion, that a Glu-His salt bridge could stabilize Cm-p5 helical conformation during the interaction with the plasma membrane. A derivative, obtained by substitution of Glu and His for Cys, was synthesized and oxidized with the generation of a cyclic monomer with improved antifungal activity. In addition, two dimers were generated during the oxidation procedure, a parallel and antiparallel one. The dimers showed a helical secondary structure in water, whereas the cyclic monomer only showed this conformation in SDS. Molecular dynamic simulations confirmed the helical stabilizations for all of them, therefore indicating the possible essential role of the Glu-His salt bridge. In addition, the antiparallel dimer showed a moderate activity against Pseudomonas aeruginosa and a significant activity against Listeria monocytogenes. Neither the cyclic monomer nor the dimers were toxic against macrophages or THP-1 human cells. Due to its increased capacity for fungal control compared to fluconazole, its low cytotoxicity, together with a stabilized α-helix and disulfide bridges, that may advance its metabolic stability, and in vivo activity, the new cyclic Cm-p5 monomer represents a potential systemic antifungal therapeutic candidate.

Dimerization of Coronavirus nsp9 with Diverse Modes Enhances Its Nucleic Acid Binding Affinity

Coronaviruses pose serious health threats to humans and other animals. Understanding the mechanisms of their replication has important implications for global health and economic stability. Nonstructural protein 9 (nsp9) is an essential RNA binding protein for coronavirus replication. However, the mechanisms of the dimerization and nucleic acid binding of nsp9 remain elusive. Here, we report four crystal structures, including wild-type porcine delta coronavirus (PDCoV) nsp9, PDCoV nsp9-ΔN7 (N-terminal 7 amino acids deleted), wild-type porcine epidemic diarrhea virus (PEDV) nsp9, and PEDV nsp9-C59A mutant. These structures reveal the diverse dimerization forms of coronavirus nsp9. We first found that the N-finger of nsp9 from PDCoV plays a critical role in dimerization. Meanwhile, PEDV nsp9 is distinguished by the presence of a disulfide bond in the dimer interface. Interestingly, size exclusion chromatography and analytical ultracentrifugation analyses indicate that the PDCoV nsp9-ΔN7 and PEDV nsp9-C59A mutants are monomeric in solution. In addition, electrophoretic mobility shift assays and microscale thermophoresis analysis indicate that the monomeric forms of PDCoV nsp9 and PEDV nsp9 still have nucleic acid binding affinity, although it is lower than that of the wild type. Our results show that the diverse dimerization forms of coronavirus nsp9 proteins enhance their nucleic acid binding affinity.IMPORTANCE Coronaviruses cause widespread respiratory, gastrointestinal, and central nervous system diseases in humans and other animals, threatening human health and causing economic loss. Coronavirus nsp9, a member of the replication complex, is an important RNA binding subunit in the RNA-synthesizing machinery of all coronaviruses. However, the mechanisms of the dimerization and nucleic acid binding of nsp9 remain elusive. In this study we determined the nsp9 crystal structures of PDCoV and PEDV. We first found that the N-finger of nsp9 from PDCoV plays a critical role in dimerization. Meanwhile, PEDV nsp9 is distinguished by the presence of a disulfide bond in the dimer interface. This study provides a structural and functional basis for understanding the mechanism of dimerization and shows that the diverse dimerization modes of coronavirus nsp9 proteins enhance their nucleic acid binding affinity. Importantly, these findings may provide a new insight for antiviral drug development.

The ins and outs of vesicular monoamine transporters

Yaffe et al. review structure-guided studies that have provided insight into the mechanism of proton-monoamine antiport by VMATs.

The H⁺-coupled vesicular monoamine transporter (VMAT) is a transporter essential for life. VMAT mediates packaging of the monoamines serotonin, dopamine, norepinephrine, and histamine from the neuronal cytoplasm into presynaptic vesicles, which is a key step in the regulated release of neurotransmitters. However, a detailed understanding of the mechanism of VMAT function has been limited by the lack of availability of high-resolution structural data. In recent years, a series of studies guided by homology models has revealed significant insights into VMAT function, identifying residues that contribute to the binding site and to specific steps in the transport cycle. Moreover, to characterize the conformational transitions that occur upon binding of the substrate and coupling ion, we have taken advantage of the unique and powerful pharmacology of VMAT as well as of mutants that affect the conformational equilibrium of the protein and shift it toward defined conformations. This has allowed us to identify an important role for the proton gradient in driving a shift from lumen-facing to cytoplasm-facing conformations.

Exploring conformational states and helical packings in the P2X receptor transmembrane domain by molecular dynamics simulation

The P2X receptor is a trimeric transmembrane protein that acts as an ATP-gated ion channel. Its transmembrane domain (TMD) contains only six helices and three of them, the M2 helices, line the ion conduction pathway. Here, using molecular dynamics simulation, I identify four conformational states of the TMD that are associated with four types of packing between M2 helices. Packing in the extracellular half of the M2 helix produces closed conformations, while packing in the intracellular half produces both open and closed conformations. State transition is observed and supports a mechanism where iris-like twisting of the M2 helices switches the location of helical packing between the extracellular and the intracellular halves of the helices. In addition, this twisting motion alters the position and orientation of residue side-chains relative to the pore and therefore influences the pore geometry and possibly ion permeation. Helical packing, on the other hand, may restrict the twisting motion and generate discrete conformational states.

The Cytokinin Oxidase/Dehydrogenase CKX1 Is a Membrane-Bound Protein Requiring Homooligomerization in the Endoplasmic Reticulum for Its Cellular Activity

Degradation of the plant hormone cytokinin is controlled by cytokinin oxidase/dehydrogenase (CKX) enzymes. The molecular and cellular behavior of these proteins is still largely unknown. In this study, we show that CKX1 is a type II single-pass membrane protein that localizes predominantly to the endoplasmic reticulum (ER) in Arabidopsis (Arabidopsis thaliana). This indicates that this CKX isoform is a bona fide ER protein directly controlling the cytokinin, which triggers the signaling from the ER. By using various approaches, we demonstrate that CKX1 forms homodimers and homooligomers in vivo. The amino-terminal part of CKX1 was necessary and sufficient for the protein oligomerization as well as for targeting and retention in the ER. Moreover, we show that protein-protein interaction is largely facilitated by transmembrane helices and depends on a functional GxxxG-like interaction motif. Importantly, mutations rendering CKX1 monomeric interfere with its steady-state localization in the ER and cause a loss of the CKX1 biological activity by increasing its ER-associated degradation. Therefore, our study provides evidence that oligomerization is a crucial parameter regulating CKX1 biological activity and the cytokinin concentration in the ER. The work also lends strong support for the cytokinin signaling from the ER and for the functional relevance of the cytokinin pool in this compartment.

Comparison of NMR and crystal structures of membrane proteins and computational refinement to improve model quality

Membrane proteins are challenging to study and restraints for structure determination are typically sparse or of low resolution because the membrane environment that surrounds them leads to a variety of experimental challenges. When membrane protein structures are determined by different techniques in different environments, a natural question is "which structure is most biologically relevant?" Towards answering this question, we compiled a dataset of membrane proteins with known structures determined by both solution NMR and X-ray crystallography. By investigating differences between the structures, we found that RMSDs between crystal and NMR structures are below 5 Å in the membrane region, NMR ensembles have a higher convergence in the membrane region, crystal structures typically have a straighter transmembrane region, have higher stereo-chemical correctness, and are more tightly packed. After quantifying these differences, we used high-resolution refinement of the NMR structures to mitigate them, which paves the way for identifying and improving the structural quality of membrane proteins.

Structural basis for dimerization and RNA binding of avian infectious bronchitis virus nsp9

The potential for infection by coronaviruses (CoVs) has become a serious concern with the recent emergence of Middle East respiratory syndrome and severe acute respiratory syndrome (SARS) in the human population. CoVs encode two large polyproteins, which are then processed into 15–16 nonstructural proteins (nsps) that make significant contributions to viral replication and transcription by assembling the RNA replicase complex. Among them, nsp9 plays an essential role in viral replication by forming a homodimer that binds single‐stranded RNA. Thus, disrupting nsp9 dimerization is a potential anti‐CoV therapy. However, different nsp9 dimer forms have been reported for alpha‐ and beta‐CoVs, and no structural information is available for gamma‐CoVs. Here we determined the crystal structure of nsp9 from the avian infectious bronchitis virus (IBV), a representative gamma‐CoV that affects the economy of the poultry industry because it can infect domestic fowl. IBV nsp9 forms a homodimer via interactions across a hydrophobic interface, which consists of two parallel alpha helices near the carboxy terminus of the protein. The IBV nsp9 dimer resembles that of SARS‐CoV nsp9, indicating that this type of dimerization is conserved among all CoVs. This makes disruption of the dimeric interface an excellent strategy for developing anti‐CoV therapies. To facilitate this effort, we characterized the roles of six conserved residues on this interface using site‐directed mutagenesis and a multitude of biochemical and biophysical methods. We found that three residues are critical for nsp9 dimerization and its abitlity to bind RNA.

PDB Code(s): 5C94

Evidence for Amino Acid Snorkeling from a High-Resolution, In Vivo Analysis of Fis1 Tail-Anchor Insertion at the Mitochondrial Outer Membrane

Proteins localized to mitochondria by a carboxyl-terminal tail anchor (TA) play roles in apoptosis, mitochondrial dynamics, and mitochondrial protein import. To reveal characteristics of TAs that may be important for mitochondrial targeting, we focused our attention upon the TA of the Saccharomyces cerevisiaeFis1 protein. Specifically, we generated a library of Fis1p TA variants fused to the Gal4 transcription factor, then, using next-generation sequencing, revealed which Fis1p TA mutations inhibited membrane insertion and allowed Gal4p activity in the nucleus. Prompted by our global analysis, we subsequently analyzed the ability of individual Fis1p TA mutants to localize to mitochondria. Our findings suggest that the membrane-associated domain of the Fis1p TA may be bipartite in nature, and we encountered evidence that the positively charged patch at the carboxyl terminus of Fis1p is required for both membrane insertion and organelle specificity. Furthermore, lengthening or shortening of the Fis1p TA by up to three amino acids did not inhibit mitochondrial targeting, arguing against a model in which TA length directs insertion of TAs to distinct organelles. Most importantly, positively charged residues were more acceptable at several positions within the membrane-associated domain of the Fis1p TA than negatively charged residues. These findings, emerging from the first high-resolution analysis of an organelle targeting sequence by deep mutational scanning, provide strong, in vivo evidence that lysine and arginine can “snorkel,” or become stably incorporated within a lipid bilayer by placing terminal charges of their side chains at the membrane interface.

Helix perturbations in membrane proteins assist in inter-helical interactions and optimal helix positioning in the bilayer

Transmembrane (TM) helices in integral membrane proteins are primarily α-helical in structure. Here we analyze 1134 TM helices in 90 high resolution membrane proteins and find that apart from the widely prevalent α-helices, TM regions also contain stretches of 3₁₀ (3 to 8 residues) and π-helices (5 to 19 residues) with distinct sequence signatures. The various helix perturbations in TM regions comprise of helices with kinked geometry, as well as those with an interspersed 3₁₀/π-helical fragment and show high occurrence in a few membrane proteins. Proline is frequently present at sites of these perturbations, but it is neither a necessary nor a sufficient requirement. Helix perturbations are also conserved within a family of membrane proteins despite low sequence identity in the perturbed region. Furthermore, a perturbation influences the geometry of the TM helix, mediates inter-helical interactions within and across protein chains and avoids hydrophobic mismatch of the helix termini with the bilayer. An analysis of π-helices in the TM regions of the heme copper oxidase superfamily shows that interspersed π-helices can vary in length from 6 to 19 amino acids or be entirely absent, depending upon the protein function. The results presented here would be helpful for prediction of 3₁₀ and π-helices in TM regions and can assist the computational design of membrane proteins.

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 .

Solution NMR Structure Determination of Polytopic α-Helical Membrane Proteins: A Guide to Spin Label Paramagnetic Relaxation Enhancement Restraints

Solution nuclear magnetic resonance structures of polytopic α-helical membrane proteins require additional restraints beyond the traditional Nuclear Overhauser Effect (NOE) restraints. Several methods have been developed and this review focuses on paramagnetic relaxation enhancement (PRE). Important aspects of spin labeling, PRE measurements, structure calculations, and structural quality are discussed.

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.

Folding and stability of the aquaglyceroporin GlpF: Implications for human aqua(glycero)porin diseases

Aquaporins are highly selective polytopic transmembrane channel proteins that facilitate the permeation of water across cellular membranes in a large diversity of organisms. Defects in aquaporin function are associated with common diseases, such as nephrogenic diabetes insipidus, congenital cataract and certain types of cancer. In general, aquaporins have a highly conserved structure; from prokaryotes to humans. The conserved structure, together with structural dynamics and the structural framework for substrate selectivity is discussed. The folding pathway of aquaporins has been a topic of several studies in recent years. These studies revealed that a conserved protein structure can be reached by following different folding pathways. Based on the available data, we suggest a complex folding pathway for aquaporins, starting from the insertion of individual helices up to the formation of the tetrameric aquaporin structure. The consequences of some known mutations in human aquaporin-encoding genes, which most likely affect the folding and stability of human aquaporins, are discussed.

Life at the border: adaptation of proteins to anisotropic membrane environment

This review discusses main features of transmembrane (TM) proteins which distinguish them from water-soluble proteins and allow their adaptation to the anisotropic membrane environment. We overview the structural limitations on membrane protein architecture, spatial arrangement of proteins in membranes and their intrinsic hydrophobic thickness, co-translational and post-translational folding and insertion into lipid bilayers, topogenesis, high propensity to form oligomers, and large-scale conformational transitions during membrane insertion and transport function. Special attention is paid to the polarity of TM protein surfaces described by profiles of dipolarity/polarizability and hydrogen-bonding capacity parameters that match polarity of the lipid environment. Analysis of distributions of Trp resides on surfaces of TM proteins from different biological membranes indicates that interfacial membrane regions with preferential accumulation of Trp indole rings correspond to the outer part of the lipid acyl chain region-between double bonds and carbonyl groups of lipids. These "midpolar" regions are not always symmetric in proteins from natural membranes. We also examined the hydrophobic effect that drives insertion of proteins into lipid bilayer and different free energy contributions to TM protein stability, including attractive van der Waals forces and hydrogen bonds, side-chain conformational entropy, the hydrophobic mismatch, membrane deformations, and specific protein-lipid binding.

Terminal residue hydrophobicity modulates transmembrane helix-helix interactions

Central to the formation of tertiary structure in membrane protein folding is the presence of amino acid sequence motifs (such as "small-XXX-small" segments) in the TM segments that promote interaction-compatible surfaces through which the TM α-helices interact. Here, we sought to elucidate additional factors that may work in tandem to dictate the ultimate interaction fate of TM-embedded segments. In this context, we used proteolipid protein (PLP), the major protein from central nervous system myelin for which mutant-dependent non-native oligomerization has been implicated in neurological disorders, to explore the specific effects of TM boundary residues (the membrane entry and exit points), keying on the secondary structure and self-association of peptides corresponding to the PLP TM2 α-helix (wild-type sequence ⁶⁶AFQYVIYGTASFFFLYGALLLAEGF⁹⁰). Using gel electrophoresis, circular dichroism, and Förster resonance energy transfer in the membrane-mimetic detergent sodium dodecyl sulfate (SDS), we found that mutation of F90 to residues such as A, I, L, or V maintains the onset of TM2-TM2 dimerization, whereas mutation to E, G, Q, N, S, or T abrogates dimer formation. We attribute this sensitivity to changes in local hydrophobicity, viz., a decrease in hydrophobicity reduces local lipid-peptide interactions, which in turn disrupts peptide α-helicity and hence the effectiveness of an incipient interaction-compatible surface. Our results show that the secondary structure and oligomeric state of PLP TM2 Lys-tagged peptides are significantly modulated by the specific nature of their C-terminal boundary residue, thus providing insight as to how point mutations, particularly where they produce disease states, can compromise the folding process.

Probabilistic grammatical model for helix-helix contact site classification

Background

Hidden Markov Models power many state‐of‐the‐art tools in the field of protein bioinformatics. While excelling in their tasks, these methods of protein analysis do not convey directly information on medium‐ and long‐range residue‐residue interactions. This requires an expressive power of at least context‐free grammars. However, application of more powerful grammar formalisms to protein analysis has been surprisingly limited.

Results

In this work, we present a probabilistic grammatical framework for problem‐specific protein languages and apply it to classification of transmembrane helix‐helix pairs configurations. The core of the model consists of a probabilistic context‐free grammar, automatically inferred by a genetic algorithm from only a generic set of expert‐based rules and positive training samples. The model was applied to produce sequence based descriptors of four classes of transmembrane helix‐helix contact site configurations. The highest performance of the classifiers reached AUCROC of 0.70. The analysis of grammar parse trees revealed the ability of representing structural features of helix‐helix contact sites.

Conclusions

We demonstrated that our probabilistic context‐free framework for analysis of protein sequences outperforms the state of the art in the task of helix‐helix contact site classification. However, this is achieved without necessarily requiring modeling long range dependencies between interacting residues. A significant feature of our approach is that grammar rules and parse trees are human‐readable. Thus they could provide biologically meaningful information for molecular biologists.

MP:PD--a data base of internal packing densities, internal packing defects and internal waters of helical membrane proteins

The membrane protein packing database (MP:PD) (http://proteinformatics.charite.de/mppd) is a database of helical membrane proteins featuring internal atomic packing densities, cavities and waters. Membrane proteins are not tightly packed but contain a considerable number of internal cavities that differ in volume, polarity and solvent accessibility as well as in their filling with internal water. Internal cavities are supposed to be regions of high physical compressibility. By serving as mobile hydrogen bonding donors or acceptors, internal waters likely facilitate transition between different functional states. Despite these distinct functional roles, internal cavities of helical membrane proteins are not well characterized, mainly because most internal waters are not resolved by crystal structure analysis. Here we combined various computational biophysical techniques to characterize internal cavities, reassign positions of internal waters and calculate internal packing densities of all available helical membrane protein structures and stored them in MP:PD. The database can be searched using keywords and entries can be downloaded. Each entry can be visualized in Provi, a Jmol-based protein viewer that provides an integrated display of low energy waters alongside membrane planes, internal packing density, hydrophobic cavities and hydrogen bonds.

Robust driving forces for transmembrane helix packing

The packing structures of transmembrane helices are traditionally attributed to patterns in residues along the contact surface. In this view, besides keeping the helices confined in the membrane, the bilayer has only a minor effect on the helices structure. Here, we use two different approaches to show that the lipid environment has a crucial effect in determining the cross-angle distribution of packed helices. We analyzed structural data of a membrane proteins database. We show that the distribution of cross angles of helix pairs in this database is statistically indistinguishable from the cross-angle distribution of two noninteracting helices imbedded in the membrane. These results suggest that the cross angle is, to a large extent, determined by the tilt angle of the individual helices. We test this hypothesis using molecular simulations of a coarse-grained model that contains no specific residue interactions. These simulations reproduce the same cross-angle distribution as found in the database. As the tilt angle of a helix is dominated by hydrophobic mismatch between the protein and surrounding lipids, our results indicate that hydrophobic mismatch is the dominant factor guiding the transmembrane helix packing. Other short-range forces might then fine-tune the structure to its final configuration.

Protein Structure in Membrane Domains

Of great interest to the academic and pharmaceutical research communities, helical transmembrane proteins are characterized by their ability to dissolve and fold in lipid bilayers—properties conferred by polypeptide spans termed transmembrane domains (TMDs). The apolar nature of TMDs necessitates the use of membrane-mimetic solvents for many structure and folding studies. This review examines the relationship between TMD structure and solvent environment, focusing on principles elucidated largely in membrane-mimetic environments with single-TMD protein and peptide models. Following a brief description of TMD sequence and conformational characteristics gleaned from the structural database, we present an overview of the conceptual models used to study folding in vitro. The impact of sequence and solvent context on the incorporation of TMDs into membranes, and its role in measurements of TMD self-assembly strengths, is then described. We conclude with a discussion of the nonspecific effects of membrane components on TMD stability.

An amino acid packing code for α-helical structure and protein design

This work demonstrates that all packing in α-helices can be simplified to repetitive patterns of a single motif: the knob-socket. Using the precision of Voronoi Polyhedra/Delauney Tessellations to identify contacts, the knob-socket is a four-residue tetrahedral motif: a knob residue on one α-helix packs into the three-residue socket on another α-helix. The principle of the knob-socket model relates the packing between levels of protein structure: the intra-helical packing arrangements within secondary structure that permit inter-helix tertiary packing interactions. Within an α-helix, the three-residue sockets arrange residues into a uniform packing lattice. Inter-helix packing results from a definable pattern of interdigitated knob-socket motifs between two α-helices. Furthermore, the knob-socket model classifies three types of sockets: (1) free, favoring only intra-helical packing; (2) filled, favoring inter-helical interactions; and (3) non, disfavoring α-helical structure. The amino acid propensities in these three socket classes essentially represent an amino acid code for structure in α-helical packing. Using this code, we used a novel yet straightforward approach for the design of α-helical structure to validate the knob-socket model. Unique sequences for three peptides were created to produce a predicted amount of α-helical structure: mostly helical, some helical, and no helix. These three peptides were synthesized, and helical content was assessed using CD spectroscopy. The measured α-helicity of each peptide was consistent with the expected predictions. These results and analysis demonstrate that the knob-socket motif functions as the basic unit of packing and presents an intuitive tool to decipher the rules governing packing in protein structure.

Comparative properties of caveolar and noncaveolar preparations of kidney Na+/K+-ATPase

To evaluate previously proposed functions of renal caveolar Na⁺/K⁺-ATPase, we modified the standard procedures for the preparation of the purified membrane-bound kidney enzyme, separated the caveolar and noncaveolar pools, and compared their properties. While the subunits of Na⁺/K⁺-ATPase (α,β,γ) constituted most of the protein content of the noncaveolar pool, the caveolar pool also contained caveolins and major caveolar proteins annexin-2 tetramer and E-cadherin. Ouabain-sensitive Na⁺/K⁺-ATPase activities of the two pools had similar properties and equal molar activities, indicating that the caveolar enzyme retains its ion transport function and does not contain nonpumping enzyme. As minor constituents, both caveolar and noncaveolar pools also contained Src, EGFR, PI3K, and several other proteins known to be involved in stimulous-induced signaling by Na⁺/K⁺-ATPase, indicating that signaling function is not limited to the caveolar pool. Endogenous Src was active in both pools but was not further activated by ouabain, calling into question direct interaction of Src with native Na⁺/K⁺-ATPase. Chemical cross-linking, co-immunoprecipitation, and immunodetection studies showed that in the caveolar pool, caveolin-1 oligomers, annexin-2 tetramers, and oligomers of the α,β,γ-protomers of Na⁺/K⁺-ATPase form a large multiprotein complex. In conjunction with known roles of E-cadherin and the β-subunit of Na⁺/K⁺-ATPase in cell adhesion and noted intercellular β,β-contacts within the structure of Na⁺/K⁺-ATPase, our findings suggest that interacting caveolar Na⁺/K⁺-ATPases located at renal adherens junctions maintain contact of two adjacent cells, conduct essential ion pumping, and are capable of locus-specific signaling in junctional cells.

Stability of mitochondrial membrane proteins in terrestrial vertebrates predicts aerobic capacity and longevity

The cellular energy produced by mitochondria is a fundamental currency of life. However, the extent to which mitochondrial (mt) performance (power and endurance) is adapted to habitats and life strategies of vertebrates is not well understood. A global analysis of mt genomes revealed that hydrophobicity (HYD) of mt membrane proteins (MMPs) is much lower in terrestrial vertebrates than in fishes and shows a strong negative correlation with serine/threonine composition (STC). Here, we present evidence that this systematic feature of MMPs was crucial for the evolution of large terrestrial vertebrates with high aerobic capacity. An Arrhenius-type equation gave positive correlations between STC and maximum life span (MLS) in terrestrial vertebrates (with a few exceptions relating to the lifestyle of small animals with a high resting metabolic rate [RMR]) and negative correlations in secondary marine vertebrates, such as cetaceans and alligators (which returned from land to water, utilizing buoyancy with increased body size). In particular, marked STC increases in primates (especially hominoids) among placentals were associated with very high MLS values. We connected these STC increases in MMPs with greater stability of respiratory complexes by estimating the degradation of the Arrhenius plot given by accelerating mtRMR up to mt maximum metabolic rate. Both mtRMR and HYD in terrestrial vertebrates decreased with increasing body mass. Decreases in mtRMR raise MMP stability when high mobility is not required, whereas decreased HYD may weaken this stability under the hydrophobic environment of lipid bilayer. High maximal metabolic rates (5–10 RMR), which we postulate require high MMP mobility, presumably render MMPs more unstable. A marked rise in STC may therefore be essential to stabilize MMPs, perhaps as dynamic supercomplexes, via hydrogen bonds associated with serine/threonine motifs.

TMKink: a method to predict transmembrane helix kinks

A hallmark of membrane protein structure is the large number of distorted transmembrane helices. Because of the prevalence of bends, it is important to not only understand how they are generated but also to learn how to predict their occurrence. Here, we find that there are local sequence preferences in kinked helices, most notably a higher abundance of proline, which can be exploited to identify bends from local sequence information. A neural network predictor identifies over two-thirds of all bends (sensitivity 0.70) with high reliability (specificity 0.89). It is likely that more structural data will allow for better helix distortion predictors with increased coverage in the future. The kink predictor, TMKink, is available at http://tmkinkpredictor.mbi.ucla.edu/.

TMPad: an integrated structural database for helix-packing folds in transmembrane proteins

α-Helical transmembrane (TM) proteins play an important role in many critical and diverse biological processes, and specific associations between TM helices are important determinants for membrane protein folding, dynamics and function. In order to gain insights into the above phenomena, it is necessary to investigate different types of helix-packing modes and interactions. However, such information is difficult to obtain because of the experimental impediment and a lack of a well-annotated source of helix-packing folds in TM proteins. We have developed the TMPad (TransMembrane Protein Helix-Packing Database) which addresses the above issues by integrating experimentally observed helix–helix interactions and related structural information of membrane proteins. Specifically, the TMPad offers pre-calculated geometric descriptors at the helix-packing interface including residue backbone/side-chain contacts, interhelical distances and crossing angles, helical translational shifts and rotational angles. The TMPad also includes the corresponding sequence, topology, lipid accessibility, ligand-binding information and supports structural classification, schematic diagrams and visualization of the above structural features of TM helix-packing. Through detailed annotations and visualizations of helix-packing, this online resource can serve as an information gateway for deciphering the relationship between helix–helix interactions and higher levels of organization in TM protein structure and function. The website of the TMPad is freely accessible to the public at http://bio-cluster.iis.sinica.edu.tw/TMPad.

Statistical analysis and exposure status classification of transmembrane beta barrel residues

Several computational methods exist for the identification of transmembrane beta barrel proteins (TMBs) from sequence. Some of these methods also provide the transmembrane (TM) boundaries of the putative TMBs. The aim of this study is to (1) derive the propensities of the TM residues to be exposed to the lipid bilayer and (2) to predict the exposure status (i.e. exposed to the bilayer or hidden in protein structure) of TMB residues. Three novel propensity scales namely, BTMC, BTMI and HTMI were derived for the TMB residues at the hydrophobic core region of the outer membrane (OM), the lipid-water interface regions of the OM, and for the helical membrane proteins (HMPs) residues at the lipid-water interface regions of the inner membrane (IM), respectively. Separate propensity scales were derived for monomeric and functionally oligomeric TMBs. The derived propensities reflect differing physico-chemical properties of the respective membrane bilayer regions and were employed in a computational method for the prediction of the exposure status of TMB residues. Based on the these propensities, the conservation indices and the frequency profile of the residues, the transmembrane residues were classified into buried/exposed with an accuracy of 77.91% and 80.42% for the residues at the membrane core and the interface regions, respectively. The correlation of the derived scales with different physico-chemical properties obtained from the AAIndex database are also discussed. Knowledge about the residue propensities and burial status will be useful in annotating putative TMBs with unknown structure.

Driving forces for transmembrane alpha-helix oligomerization

We present what we believe to be a novel statistical contact potential based on solved structures of transmembrane (TM) alpha-helical bundles, and we use this contact potential to investigate the amino acid likelihood of stabilizing helix-helix interfaces. To increase statistical significance, we have reduced the full contact energy matrix to a four-flavor alphabet of amino acids, automatically determined by our methodology, in which we find that polarity is a more dominant factor of group identity than is size, with charged or polar groups most often occupying the same face, whereas polar/apolar residue pairs tend to occupy opposite faces. We found that the most polar residues strongly influence interhelical contact formation, although they occur rarely in TM helical bundles. Two-body contact energies in the reduced letter code are capable of determining native structure from a large decoy set for a majority of test TM proteins, at the same time illustrating that certain higher-order sequence correlations are necessary for more accurate structure predictions.

Interpretation of 2H-NMR experiments on the orientation of the transmembrane helix WALP23 by computer simulations

Orientation, dynamics, and packing of transmembrane helical peptides are important determinants of membrane protein structure, dynamics, and function. Because it is difficult to investigate these aspects by studying real membrane proteins, model transmembrane helical peptides are widely used. NMR experiments provide information on both orientation and dynamics of peptides, but they require that motional models be interpreted. Different motional models yield different interpretations of quadrupolar splittings (QS) in terms of helix orientation and dynamics. Here, we use coarse-grained (CG) molecular dynamics (MD) simulations to investigate the behavior of a well-known model transmembrane peptide, WALP23, under different hydrophobic matching/mismatching conditions. We compare experimental (2)H-NMR QS (directly measured in experiments), as well as helix tilt angle and azimuthal rotation (not directly measured), with CG MD simulation results. For QS, the agreement is significantly better than previously obtained with atomistic simulations, indicating that equilibrium sampling is more important than atomistic details for reproducing experimental QS. Calculations of helix orientation confirm that the interpretation of QS depends on the motional model used. Our simulations suggest that WALP23 can form dimers, which are more stable in an antiparallel arrangement. The origin of the preference for the antiparallel orientation lies not only in electrostatic interactions but also in better surface complementarity. In most cases, a mixture of monomers and antiparallel dimers provides better agreement with NMR data compared to the monomer and the parallel dimer. CG MD simulations allow predictions of helix orientation and dynamics and interpretation of QS data without requiring any assumption about the motional model.

Predicting transmembrane helix packing arrangements using residue contacts and a force-directed algorithm

Alpha-helical transmembrane proteins constitute roughly 30% of a typical genome and are involved in a wide variety of important biological processes including cell signalling, transport of membrane-impermeable molecules and cell recognition. Despite significant efforts to predict transmembrane protein topology, comparatively little attention has been directed toward developing a method to pack the helices together. Here, we present a novel approach to predict lipid exposure, residue contacts, helix-helix interactions and finally the optimal helical packing arrangement of transmembrane proteins. Using molecular dynamics data, we have trained and cross-validated a support vector machine (SVM) classifier to predict per residue lipid exposure with 69% accuracy. This information is combined with additional features to train a second SVM to predict residue contacts which are then used to determine helix-helix interaction with up to 65% accuracy under stringent cross-validation on a non-redundant test set. Our method is also able to discriminate native from decoy helical packing arrangements with up to 70% accuracy. Finally, we employ a force-directed algorithm to construct the optimal helical packing arrangement which demonstrates success for proteins containing up to 13 transmembrane helices. This software is freely available as source code from http://bioinf.cs.ucl.ac.uk/memsat/mempack/.

Alpha-helical transmembrane proteins constitute a significant proportion of the proteins encoded in a typical genome and are involved in a wide variety of important biological processes. Many common diseases including diabetes, hypertension and epilepsy have been related to transmembrane protein dysfunction, therefore they represent one of the most important classes of protein for pharmaceutical intervention. However, due to the experimental difficulties of structure determination, this class of protein is severely under-represented in structural databases. Here, we present a novel approach that is able to predict lipid exposure, residue contacts, helix-helix interactions and finally the optimal helical packing arrangement of a transmembrane protein. Under stringent cross-validation, our approach demonstrates a significant improvement in prediction over existing software. This method can be used to gain insights into transmembrane protein folding and enhance the quality of ab initio modelling, while providing testable hypotheses for a variety of studies including protein design, mutagenesis and thermostability experiments.

[Prediction of the spatial structure of proteins: emphasis on membrane targets]

Knowledge of the spatial structure of proteins is a prerequisite for both awareness of their functional mechanisms and the framework for rational drug discovery and design. Meanwhile, direct structural determination is often hampered or impractical due to the complexity, expensiveness, and limited capabilities of experimental techniques. These issues are especially pronounced for integral membrane proteins. On numerous occasions, the theoretical prediction of protein structures may facilitate the process by exploiting physical or empirical principles. This paper surveys modern techniques for the prediction of the spatial structure of proteins using computer algorithms, and the main emphasis is placed on the most "complex" targets - membrane proteins (MPs). The first part of the review describes de novo methods based on empirical physical principles; in the second part, a comparative modeling philosophy, which accounts for the structure of related proteins, is described. Special focus is made regarding pharmacologically relevant classes of G-coupled receptors, receptor tyrosine ki-nases, and other MPs. Algorithms for the assessment of the models quality and potential fields of application of computer models are discussed.

Critical review of general guidelines for membrane proteins model building and analysis

Membrane proteins play major roles in many biological processes such as signalling, transport, etc. They have been shown to be involved in the development of many diseases and have become important drug targets per se. The understanding of their functional properties may be facilitated if a 3D structure is available. However, in the case of membrane proteins, only a few 3D structures have been solved to date. Bioinformatics and molecular modelling approaches are thus powerful alternatives to fill the gap between the sequence and the structure. Here, a review of the most recent approaches is proposed together with guidelines on how to use them. In addition, insofar as important biological processes require conformational changes, we discuss some interesting methods aimed at exploring the dynamic behaviour of proteins in their membrane environment. The paper ends with a brief description of useful approaches for determining oligomerisation or ligand binding sites.

An intramembrane aromatic network determines pentameric assembly of Cys-loop receptors

Cys-loop receptors are pentameric ligand-gated ion channels (pLGICs) that mediate fast synaptic transmission. Here functional pentameric assembly of truncated fragments comprising the ligand-binding N-terminal ectodomains and the first three transmembrane helices, M1-M3, of both the inhibitory glycine receptor (GlyR) alpha1 and the 5HT(3)A receptor subunits was found to be rescued by coexpressing the complementary fourth transmembrane helix, M4. Alanine scanning identified multiple aromatic residues in M1, M3 and M4 as key determinants of GlyR assembly. Homology modeling and molecular dynamics simulations revealed that these residues define an interhelical aromatic network, which we propose determines the geometry of M1-M4 tetrahelical packing such that nascent pLGIC subunits must adopt a closed fivefold symmetry. Because pLGIC ectodomains form random nonstoichiometric oligomers, proper pentameric assembly apparently depends on intersubunit interactions between extracellular domains and intrasubunit interactions between transmembrane segments.