TMKink: a method to predict transmembrane helix kinks

SCOT: Rethinking the classification of secondary structure elements

MOTIVATION

Secondary structure classification is one of the most important issues in structure-based analyses due to its impact on secondary structure prediction, structural alignment and protein visualization. There are still open challenges concerning helix and sheet assignments which are currently not addressed by a single multi-purpose software.

RESULTS

We introduce SCOT (Secondary structure Classification On Turns) as a novel secondary structure element assignment software which supports the assignment of turns, right-handed α-, 310- and π-helices, left-handed α- and 310-helices, 2.27- and polyproline II helices, β-sheets and kinks. We demonstrate that the introduction of helix Purity values enables a clear differentiation between helix classes. SCOT's unique strengths are highlighted by comparing it to six state-of-the-art methods (DSSP, STRIDE, ASSP, SEGNO, DISICL and SHAFT). The assignment approaches were compared concerning geometric consistency, protein structure quality and flexibility dependency and their impact on secondary structure element-based structural alignments. We show that only SCOT's combination of hydrogen bonds, geometric criteria and dihedral angles enables robust assignments independent of the structure quality and flexibility. We demonstrate that this combination and the elaborate kink detection lead to SCOT's clear superiority for protein alignments. As the resulting helices and strands are provided in a PDB conform output format, they can immediately be used for structure alignment algorithms. Taken together, the application of our new method and the straight-forward visualization using the accompanying PyMOL scripts enable the comprehensive analysis of regular backbone geometries in proteins.

AVAILABILITY AND IMPLEMENTATION

https://this-group.rocks.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Early adolescent brain markers of late adolescent academic functioning

Academic performance in adolescence strongly influences adult prospects. Intelligence quotient (IQ) has historically been considered a strong predictor of academic performance. Less objectively explored have been morphometric features. We analyzed brain MRI morphometry metrics in early adolescence (age 12-14 years) as quantitative predictors of academic performance over high school using a naïve Bayesian classifier approach with n = 170 subjects. Based on the mean GPA, subjects were divided into high (GPA ≥3.54; n = 87) and low (GPA <3.54; n = 83) academic performers. Covariance analysis was performed to look at the influence of subject demographics. We examined predictive features from the 343 available regions (surface areas, cortical thickness, and subcortical volumes) and applied 4 algorithms for selection and reduction of attributes using Weka. Cortical thickness measures performed better than surface areas or subcortical volumes as predictors of academic performance. We identified 15 cortical thickness regions most predictive of academic performance, three of which have not been described in the literature predictive of academic performance. These were in the left hemisphere fusiform, bilateral insula, and left hemisphere paracentral regions. Prediction had a sensitivity of 0.65 and specificity of 0.73 with independent validation. Follow-up independent t-test analyses between high and low academic achievers on 10 of 15 regions showed between-group significance at the p < 0.05 level. High achievers demonstrated thicker cortices than low achievers. These newly identified regions may help pinpoint new targets for further study in understanding the developing adolescent brain in the classroom setting.

Curvature induction and membrane remodeling by FAM134B reticulon homology domain assist selective ER-phagy

FAM134B/RETREG1 is a selective ER-phagy receptor that regulates the size and shape of the endoplasmic reticulum. The structure of its reticulon-homology domain (RHD), an element shared with other ER-shaping proteins, and the mechanism of membrane shaping remain poorly understood. Using molecular modeling and molecular dynamics (MD) simulations, we assemble a structural model for the RHD of FAM134B. Through MD simulations of FAM134B in flat and curved membranes, we relate the dynamic RHD structure with its two wedge-shaped transmembrane helical hairpins and two amphipathic helices to FAM134B functions in membrane-curvature induction and curvature-mediated protein sorting. FAM134B clustering, as expected to occur in autophagic puncta, amplifies the membrane-shaping effects. Electron microscopy of in vitro liposome remodeling experiments support the membrane remodeling functions of the different RHD structural elements. Disruption of the RHD structure affects selective autophagy flux and leads to disease states.

FAM134B/RETREG1 is a selective ER-phagy receptor that regulates the size and shape of the endoplasmic reticulum. Here authors use molecular modeling and molecular dynamics simulations to assemble a structural model for the reticulon-homology domain of FAM134B.

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

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

Chiral Ramachandran Plots II: General Trends and Protein Chirality Spectra

The degree of chirality of protein backbone residues is used to enrich the Ramachandran plot (RP) and create three-dimensional chiral RPs with much more structural information. Detailed comparative analysis of the four classical RPs (general, glycine, proline, and pre-proline) is provided, including statistical analysis of quantitative chirality distributions in the maps and in the secondary structures. Our results show that points with outlier chirality levels represent special transitional points in the folded protein such as α-helix kinks, twists of β-strands, and transition points between secondary structures. A protein chirality spectrum in which the degree of chirality of each residue is plotted against the sequence number explores these special points. More than 65000 residues extracted from 200 high-quality proteins are used for this study, which shows that quantitative chirality is a general and useful structural parameter for protein conformational studies.

A statistical model for helices with applications

Motivated by a cutting edge problem related to the shape of α -helices in proteins, we formulate a parametric statistical model, which incorporates the cylindrical nature of the helix. Our focus is to detect a "kink," which is a drastic change in the axial direction of the helix. We propose a statistical model for the straight α -helix and derive the maximum likelihood estimation procedure. The cylinder is an accepted geometric model for α -helices, but our statistical formulation, for the first time, quantifies the uncertainty in atom positions around the cylinder. We propose a change point technique "Kink-Detector" to detect a kink location along the helix. Unlike classical change point problems, the change in direction of a helix depends on a simultaneous shift of multiple data points rather than a single data point, and is less straightforward. Our biological building block is crowdsourced data on straight and kinked helices; which has set a gold standard. We use this data to identify salient features to construct Kink-detector, test its performance and gain some insights. We find the performance of Kink-detector comparable to its computational competitor called "Kink-Finder." We highlight that identification of kinks by visual assessment can have limitations and Kink-detector may help in such cases. Further, an analysis of crowdsourced curved α -helices finds that Kink-detector is also effective in detecting moderate changes in axial directions.

Backbone Hydrogen Bond Strengths Can Vary Widely in Transmembrane Helices

Although backbone hydrogen bonds in transmembrane (TM) helices have the potential to be very strong due to the low dielectric and low water environment of the membrane, their strength has never been assessed experimentally. Moreover, variations in hydrogen bond strength might be necessary to facilitate the TM helix breaking and bending that is often needed to satisfy functional imperatives. Here we employed equilibrium hydrogen/deuterium fractionation factors to measure backbone hydrogen bond strengths in the TM helix of the amyloid precursor protein (APP). We find an enormous range of hydrogen bond free energies, with some weaker than water–water hydrogen bonds and some over 6 kcal/mol stronger than water–water hydrogen bonds. We find that weak hydrogen bonds are at or near preferred γ-secretase cleavage sites, suggesting that the sequence of APP and possibly other cleaved TM helices may be designed, in part, to make their backbones accessible for cleavage. The finding that hydrogen bond strengths in a TM helix can vary widely has implications for membrane protein function, dynamics, evolution, and design.

Molecular Properties of Globin Channels and Pores: Role of Cholesterol in Ligand Binding and Movement

Globins contain one or more cavities that control or affect such functions as ligand movement and ligand binding. Here we report that the extended globin family [cytoglobin (Cygb); neuroglobin (Ngb); myoglobin (Mb); hemoglobin (Hb) subunits Hba(α); and Hbb(β)] contain either a transmembrane (TM) helix or pore-lining region as well as internal cavities. Protein motif/domain analyses indicate that Ngb and Hbb each contain 5 cholesterol- binding (CRAC/CARC) domains and 1 caveolin binding motif, whereas the Cygb dimer has 6 cholesterol-binding domains but lacks caveolin-binding motifs. Mb and Hba each exhibit 2 cholesterol-binding domains and also lack caveolin-binding motifs. The Hb αβ-tetramer contains 14 cholesterol-binding domains. Computer algorithms indicate that Cygb and Ngb cavities display multiple partitions and C-terminal pore-lining regions, whereas Mb has three major cavities plus a C-terminal pore-lining region. The Hb tetramer exhibits a large internal cavity but the subunits differ in that they contain a C-terminal TM helix (Hba) and pore-lining region (Hbb). The cavities include 43 of 190 Cygb residues, 38 of 151 of Ngb residues, 55 of 154 Mb residues, and 137 of 688 residues in the Hb tetramer. Each cavity complex includes 6 to 8 residues of the TM helix or pore-lining region and CRAC/CARC domains exist within all cavities. Erythrocyte Hb αβ-tetramers are largely cytosolic but also bind to a membrane anion exchange protein, "band 3," which contains a large internal cavity and 12 TM helices (5 being pore-lining regions). The Hba TM helix may be the erythrocyte membrane "band 3" attachment site. "Band 3" contributes 4 caveolin binding motifs and 10 CRAC/CARC domains. Cholesterol binding may create lipid-disordered phases that alter globin cavities and facilitate ligand movement, permitting ion channel formation and conformational changes that orchestrate anion and ligand (O2, CO2, NO) movement within the large internal cavities and channels of the globins.

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.

Examining the Conservation of Kinks in Alpha Helices

Kinks are a structural feature of alpha-helices and many are known to have functional roles. Kinks have previously tended to be defined in a binary fashion. In this paper we have deliberately moved towards defining them on a continuum, which given the unimodal distribution of kink angles is a better description. From this perspective, we examine the conservation of kinks in proteins. We find that kink angles are not generally a conserved property of homologs, pointing either to their not being functionally critical or to their function being related to conformational flexibility. In the latter case, the different structures of homologs are providing snapshots of different conformations. Sequence identity between homologous helices is informative in terms of kink conservation, but almost equally so is the sequence identity of residues in spatial proximity to the kink. In the specific case of proline, which is known to be prevalent in kinked helices, loss of a proline from a kinked helix often also results in the loss of a kink or reduction in its kink angle. We carried out a study of the seven transmembrane helices in the GPCR family and found that changes in kinks could be related both to subfamilies of GPCRs and also, in a particular subfamily, to the binding of agonists or antagonists. These results suggest conformational change upon receptor activation within the GPCR family. We also found correlation between kink angles in different helices, and the possibility of concerted motion could be investigated further by applying our method to molecular dynamics simulations. These observations reinforce the belief that helix kinks are key, functional, flexible points in structures.

Evolution of the α-Subunit of Na/K-ATPase from Paramecium to Homo sapiens: Invariance of Transmembrane Helix Topology

Na/K-ATPase is a key plasma membrane enzyme involved in cell signaling, volume regulation, and maintenance of electrochemical gradients. The α-subunit, central to these functions, belongs to a large family of P-type ATPases. Differences in transmembrane (TM) helix topology, sequence homology, helix–helix contacts, cell signaling, and protein domains of Na/K-ATPase α-subunit were compared in fungi (Beauveria), unicellular organisms (Paramecia), primitive multicellular organisms (Hydra), and vertebrates (Xenopus, Homo sapiens), and correlated with evolution of physiological functions in the α-subunit. All α-subunits are of similar length, with groupings of four and six helices in the N- and C-terminal regions, respectively. Minimal homology was seen for protein domain patterns in Paramecium and Hydra, with high correlation between Hydra and vertebrates. Paramecium α-subunits display extensive disorder, with minimal helix contacts. Increases in helix contacts in Hydra approached vertebrates. Protein motifs known to be associated with membrane lipid rafts and cell signaling reveal significant positional shifts between Paramecium and Hydra vulgaris, indicating that regional membrane fluidity changes occur during evolution. Putative steroid binding sites overlapping TM-3 occurred in all species. Sites associated with G-protein-receptor stimulation occur both in vertebrates and amphibia but not in Hydra or Paramecia. The C-terminus moiety “KETYY,” necessary for the Na⁺ activation of pump phosphorylation, is not present in unicellular species indicating the absence of classical Na⁺/K⁺-pumps. The basic protein topology evolved earliest, followed by increases in protein domains and ordered helical arrays, correlated with appearance of α-subunit regions known to involve cell signaling, membrane recycling, and ion channel formation.

Transmembrane helices in "classical" nuclear reproductive steroid receptors: a perspective

Abstract Steroid receptors of the nuclear receptor superfamily are proposed to be either: 1) located in the cytosol and moved to the cell nucleus upon activation, 2) tethered to the inside of the plasma membrane, or 3) retained in the nucleus until free steroid hormone enters and activates specific receptors. Using computational methods to analyze peptide receptor topology, we find that the “classical” nuclear receptors for progesterone (PRB/PGR), androgen (ARB/AR) and estrogen (ER1/ESR1) contain two transmembrane helices (TMH) within their ligand-binding domains (LBD).The MEMSAT-SVM algorithm indicates that ARB and ER2 (but not PRB or ER1) contain a pore-lining (channel-forming) region which may merge with other pore-lining regions to form a membrane channel. ER2 lacks a TMH, but contains a single pore-lining region. The MemBrain algorithm predicts that PRB, ARB and ER1 each contain one TMH plus a half TMH separated by 51 amino acids.ER2 contains two half helices. The TM-2 helices of ARB, ER1 and ER2 each contain 9-13 amino acid motifs reported to translocate the receptor to the plasma membrane, as well as cysteine palmitoylation sites. PoreWalker analysis of X-ray crystallographic data identifies a pore or channel within the LBDs of ARB and ER1 and predicts that 70 and 72 residues are pore-lining residues, respectively. The data suggest that (except for ER2), cytosolic receptors become anchored to the plasma membrane following synthesis. Half-helices and pore-lining regions in turn form functional ion channels and/or facilitate passive steroid uptake into the cell. In perspective, steroid-dependent insertion of “classical” receptors containing pore-lining regions into the plasma membrane may regulate permeability to ions such as Ca²⁺, Na⁺ or K⁺, as well as facilitate steroid translocation into the nucleus.

Sequence-structure relationship study in all-α transmembrane proteins using an unsupervised learning approach

Transmembrane proteins (TMPs) are major drug targets, but the knowledge of their precise topology structure remains highly limited compared with globular proteins. In spite of the difficulties in obtaining their structures, an important effort has been made these last years to increase their number from an experimental and computational point of view. In view of this emerging challenge, the development of computational methods to extract knowledge from these data is crucial for the better understanding of their functions and in improving the quality of structural models. Here, we revisit an efficient unsupervised learning procedure, called Hybrid Protein Model (HPM), which is applied to the analysis of transmembrane proteins belonging to the all-α structural class. HPM method is an original classification procedure that efficiently combines sequence and structure learning. The procedure was initially applied to the analysis of globular proteins. In the present case, HPM classifies a set of overlapping protein fragments, extracted from a non-redundant databank of TMP 3D structure. After fine-tuning of the learning parameters, the optimal classification results in 65 clusters. They represent at best similar relationships between sequence and local structure properties of TMPs. Interestingly, HPM distinguishes among the resulting clusters two helical regions with distinct hydrophobic patterns. This underlines the complexity of the topology of these proteins. The HPM classification enlightens unusual relationship between amino acids in TMP fragments, which can be useful to elaborate new amino acids substitution matrices. Finally, two challenging applications are described: the first one aims at annotating protein functions (channel or not), the second one intends to assess the quality of the structures (X-ray or models) via a new scoring function deduced from the HPM classification.

Topological analysis of the Escherichia coli WcaJ protein reveals a new conserved configuration for the polyisoprenyl-phosphate hexose-1-phosphate transferase family

WcaJ is an Escherichia coli membrane enzyme catalysing the biosynthesis of undecaprenyl-diphosphate-glucose, the first step in the assembly of colanic acid exopolysaccharide. WcaJ belongs to a large family of polyisoprenyl-phosphate hexose-1-phosphate transferases (PHPTs) sharing a similar predicted topology consisting of an N-terminal domain containing four transmembrane helices (TMHs), a large central periplasmic loop, and a C-terminal domain containing the fifth TMH (TMH-V) and a cytosolic tail. However, the topology of PHPTs has not been experimentally validated. Here, we investigated the topology of WcaJ using a combination of LacZ/PhoA reporter fusions and sulfhydryl labelling by PEGylation of novel cysteine residues introduced into a cysteine-less WcaJ. The results showed that the large central loop and the C-terminal tail both reside in the cytoplasm and are separated by TMH-V, which does not fully span the membrane, likely forming a "hairpin" structure. Modelling of TMH-V revealed that a highly conserved proline might contribute to a helix-break-helix structure in all PHPT members. Bioinformatic analyses show that all of these features are conserved in PHPT homologues from Gram-negative and Gram-positive bacteria. Our data demonstrate a novel topological configuration for PHPTs, which is proposed as a signature for all members of this enzyme family.

Computational modeling of membrane proteins

The determination of membrane protein (MP) structures has always trailed that of soluble proteins due to difficulties in their overexpression, reconstitution into membrane mimetics, and subsequent structure determination. The percentage of MP structures in the protein databank (PDB) has been at a constant 1-2% for the last decade. In contrast, over half of all drugs target MPs, only highlighting how little we understand about drug-specific effects in the human body. To reduce this gap, researchers have attempted to predict structural features of MPs even before the first structure was experimentally elucidated. In this review, we present current computational methods to predict MP structure, starting with secondary structure prediction, prediction of trans-membrane spans, and topology. Even though these methods generate reliable predictions, challenges such as predicting kinks or precise beginnings and ends of secondary structure elements are still waiting to be addressed. We describe recent developments in the prediction of 3D structures of both α-helical MPs as well as β-barrels using comparative modeling techniques, de novo methods, and molecular dynamics (MD) simulations. The increase of MP structures has (1) facilitated comparative modeling due to availability of more and better templates, and (2) improved the statistics for knowledge-based scoring functions. Moreover, de novo methods have benefited from the use of correlated mutations as restraints. Finally, we outline current advances that will likely shape the field in the forthcoming decade.

The pore-lining regions in cytochrome c oxidases: A computational analysis of caveolin, cholesterol and transmembrane helix contributions to proton movement

Cytochrome c oxidase (CcO) is the terminal enzyme in the electron transfer chain. CcO catalyzes a four electron reduction of O2 to water at a catalytic site formed by high-spin heme (a3) and copper atoms (CuB). While it is recognized that proton movement is coupled to oxygen reduction, the proton channel(s) have not been well defined. Using computational methods developed to study protein topology, membrane channels and 3D packing arrangements within transmembrane (TM) helix arrays, we find that subunit-1 (COX-1), subunit-2 (COX-2) and subunit-3 (COX-3) contribute 139, 46 and 25 residues, respectively, to channel formation between the mitochondrial matrix and intermembrane space. Nine of 12 TM helices in COX-1, both helices in COX-2 and 5 of the 6 TM helices in COX-3 are pore-lining regions (possible channel formers). Heme a3 and the CuB sites (as well as the CuA center of COX-2) are located within the channel that includes TM-6, TM-7, TM-10 and TM-11 of COX-1 and are associated with multiple cholesterol and caveolin-binding (CB) motifs. Sequence analysis identifies five CB motifs within COX-1, two within COX-2 and four within COX-3; each caveolin containing a pore-lining helix C-terminal to a TM helix-turn-helix. Channel formation involves interaction between multiple pore-lining regions within protein subunits and/or dimers. PoreWalker analysis lends support to the D-channel model of proton translocation. Under physiological conditions, caveolins may introduce channel formers juxtaposed to those in COX-1, COX-2 and COX-3, which together with cholesterol may form channel(s) essential for proton translocation through the inner mitochondrial membrane.

SKINK: a web server for string kernel based kink prediction in α-helices

MOTIVATION

The reasons for distortions from optimal α-helical geometry are widely unknown, but their influences on structural changes of proteins are significant. Hence, their prediction is a crucial problem in structural bioinformatics. Here, we present a new web server, called SKINK, for string kernel based kink prediction. Extending our previous study, we also annotate the most probable kink position in a given α-helix sequence.

AVAILABILITY AND IMPLEMENTATION

The SKINK web server is freely accessible at http://biows-inf.zdv.uni-mainz.de/skink. Moreover, SKINK is a module of the BALL software, also freely available at www.ballview.org.

High-resolution modeling of transmembrane helical protein structures from distant homologues

Eukaryotic transmembrane helical (TMH) proteins perform a wide diversity of critical cellular functions, but remain structurally largely uncharacterized and their high-resolution structure prediction is currently hindered by the lack of close structural homologues. To address this problem, we present a novel and generic method for accurately modeling large TMH protein structures from distant homologues exhibiting distinct loop and TMH conformations. Models of the adenosine A2AR and chemokine CXCR4 receptors were first ranked in GPCR-DOCK blind prediction contests in the receptor structure accuracy category. In a benchmark of 50 TMH protein homolog pairs of diverse topology (from 5 to 12 TMHs), size (from 183 to 420 residues) and sequence identity (from 15% to 70%), the method improves most starting templates, and achieves near-atomic accuracy prediction of membrane-embedded regions. Unlike starting templates, the models are of suitable quality for computer-based protein engineering: redesigned models and redesigned X-ray structures exhibit very similar native interactions. The method should prove useful for the atom-level modeling and design of a large fraction of structurally uncharacterized TMH proteins from a wide range of structural homologues.

Helix kinks are equally prevalent in soluble and membrane proteins

Helix kinks are a common feature of α-helical membrane proteins, but are thought to be rare in soluble proteins. In this study we find that kinks are a feature of long α-helices in both soluble and membrane proteins, rather than just transmembrane α-helices. The apparent rarity of kinks in soluble proteins is due to the relative infrequency of long helices (≥20 residues) in these proteins. We compare length-matched sets of soluble and membrane helices, and find that the frequency of kinks, the role of Proline, the patterns of other amino acid around kinks (allowing for the expected differences in amino acid distributions between the two types of protein), and the effects of hydrogen bonds are the same for the two types of helices. In both types of protein, helices that contain Proline in the second and subsequent turns are very frequently kinked. However, there are a sizeable proportion of kinked helices that do not contain a Proline in either their sequence or sequence homolog. Moreover, we observe that in soluble proteins, kinked helices have a structural preference in that they typically point into the solvent.

Computational prediction of kink properties of helices in membrane proteins

We have combined molecular dynamics simulations and fold identification procedures to investigate the structure of 696 kinked and 120 unkinked transmembrane (TM) helices in the PDBTM database. Our main aim of this study is to understand the formation of helical kinks by simulating their quasi-equilibrium heating processes, which might be relevant to the prediction of their structural features. The simulated structural features of these TM helices, including the position and the angle of helical kinks, were analyzed and compared with statistical data from PDBTM. From quasi-equilibrium heating processes of TM helices with four very different relaxation time constants, we found that these processes gave comparable predictions of the structural features of TM helices. Overall, 95 % of our best kink position predictions have an error of no more than two residues and 75 % of our best angle predictions have an error of less than 15°. Various structure assessments have been carried out to assess our predicted models of TM helices in PDBTM. Our results show that, in 696 predicted kinked helices, 70 % have a RMSD less than 2 Å, 71 % have a TM-score greater than 0.5, 69 % have a MaxSub score greater than 0.8, 60 % have a GDT-TS score greater than 85, and 58 % have a GDT-HA score greater than 70. For unkinked helices, our predicted models are also highly consistent with their crystal structure. These results provide strong supports for our assumption that kink formation of TM helices in quasi-equilibrium heating processes is relevant to predicting the structure of TM helices.

Architecture of the hepatitis C virus E1 glycoprotein transmembrane domain studied by NMR

Oligomerization of hepatitis C viral envelope proteins E1 and E2 is essential to virus fusion and assembly. Although interactions within the transmembrane (TM) domains of these glycoproteins have proven contributions to the E1/E2 heterodimerization process and consequent infectivity, there is little structural information on this entry mechanism. Here, as a first step towards our long-term goal of understanding the interaction between E1 and E2 TM-domains, we have expressed, purified and characterized E1-TM using structural biomolecular NMR methods. An MBP-fusion expression system yielded sufficient quantities of pure E1-TM, which was solubilized in two membrane-mimicking environments, SDS- and LPPG-micelles, affording samples amenable to NMR studies. Triple resonance assignment experiments and relaxation measurements provided information on the secondary structure and global fold of E1-TM in these environments. In SDS micelles E1-TM adopts a helical conformation, with helical stretches at residues 354-363 and 371-379 separated by a more flexible segment of residues 364-370. In LPPG micelles a helical conformation was observed for residues 354-377 with greater flexibility in the 366-367 dyad, suggesting LPPG provides a more native environment for the peptide. Replacement of key positively charged residue K370 with an alanine did not affect the secondary structure of E1-TM but did change the relative positioning within the micelle of the two helices. These results lay the foundation for structure determination of E1-TM and a molecular understanding of how E1-TM flexibility enhances its interaction with E2-TM during heterodimerization and membrane fusion.

The glove-like structure of the conserved membrane protein TatC provides insight into signal sequence recognition in twin-arginine translocation

In bacteria, two signal-sequence-dependent secretion pathways translocate proteins across the cytoplasmic membrane. Although the mechanism of the ubiquitous general secretory pathway is becoming well understood, that of the twin-arginine translocation pathway, responsible for translocation of folded proteins across the bilayer, is more mysterious. TatC, the largest and most conserved of three integral membrane components, provides the initial binding site of the signal sequence prior to pore assembly. Here, we present two crystal structures of TatC from the thermophilic bacteria Aquifex aeolicus at 4.0 Å and 6.8 Å resolution. The membrane architecture of TatC includes a glove-shaped structure with a lipid-exposed pocket predicted by molecular dynamics to distort the membrane. Correlating the biochemical literature to these results suggests that the signal sequence binds in this pocket, leading to structural changes that facilitate higher order assemblies.

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.

Assembly of the transmembrane domain of E. coli PhoQ histidine kinase: implications for signal transduction from molecular simulations

The PhoQP two-component system is a signaling complex essential for bacterial virulence and cationic antimicrobial peptide resistance. PhoQ is the histidine kinase chemoreceptor of this tandem machine and assembles in a homodimer conformation spanning the bacterial inner membrane. Currently, a full understanding of the PhoQ signal transduction is hindered by the lack of a complete atomistic structure. In this study, an atomistic model of the key transmembrane (TM) domain is assembled by using molecular simulations, guided by experimental cross-linking data. The formation of a polar pocket involving Asn202 in the lumen of the tetrameric TM bundle is crucial for the assembly and solvation of the domain. Moreover, a concerted displacement of the TM helices at the periplasmic side is found to modulate a rotation at the cytoplasmic end, supporting the transduction of the chemical signal through a combination of scissoring and rotational movement of the TM helices.

Two-component systems (TCSs) are signaling complexes essential for bacterial survival and virulence. PhoQ is the histidine kinase chemoreceptor of the PhoQ-PhoP tandem machine that detects the concentration of cationic species at the inner membrane of Gram-negative bacteria. A full understanding of the PhoQ signal transduction mechanism is currently hindered by the lack of a complete atomistic structure. Here, by using molecular simulations integrated with cross-linking disulfide scanning data, we present the first structural model of the transmembrane (TM) portion of PhoQ from E. coli. Its structural and dynamic features induce a concerted displacement of the TM helices at the periplasmic side, which modulates a rotation at the cytoplasmic end. This supports the idea that signal transduction is promoted through a combination of scissoring and rotational movements of the TM helices. This complex mechanism is the key to understanding how the chemical stimuli sensed by the periplasmic sensor domain trigger, via the relay of the HAMP domain, the histidine auto-phosphorylation and kinase/phosphatase activity at the cytoplasmic end.

G protein-coupled receptors--recent advances

The years 2000 and 2007 witnessed milestones in current understanding of G protein-coupled receptor (GPCR) structural biology. In 2000 the first GPCR, bovine rhodopsin, was crystallized and the structure was solved, while in 2007 the structure of β(2)-adrenergic receptor, the first GPCR with diffusible ligands, was determined owing to advances in microcrystallization and an insertion of the fast-folding lysozyme into the receptor. In parallel with those crystallographic studies, the biological and biochemical characterization of GPCRs has advanced considerably because those receptors are molecular targets for many of currently used drugs. Therefore, the mechanisms of activation and signal transduction to the cell interior deduced from known GPCRs structures are of the highest importance for drug discovery. These proteins are the most diversified membrane receptors encoded by hundreds of genes in our genome. They participate in processes responsible for vision, smell, taste and neuronal transmission in response to photons or binding of ions, hormones, peptides, chemokines and other factors. Although the GPCRs share a common seven-transmembrane α-helical bundle structure their binding sites can accommodate thousands of different ligands. The ligands, including agonists, antagonists or inverse agonists change the structure of the receptor. With bound agonists they can form a complex with a suitable G protein, be phosphorylated by kinases or bind arrestin. The discovered signaling cascades invoked by arrestin independently of G proteins makes the GPCR activating scheme more complex such that a ligand acting as an antagonist for G protein signaling can also act as an agonist in arrestin-dependent signaling. Additionally, the existence of multiple ligand-dependent partial activation states as well as dimerization of GPCRs result in a 'microprocessor-like' action of these receptors rather than an 'on-off' switch as was commonly believed only a decade ago.

A transmembrane domain and GxxxG motifs within L2 are essential for papillomavirus infection

During cellular invasion, human papillomavirus type 16 (HPV16) must transfer its viral genome (vDNA) across the endosomal membrane prior to its accumulation at nuclear PML bodies for the establishment of infection. After cellular uptake, the capsid likely undergoes pH-dependent disassembly within the endo-/lysosomal compartment, thereby exposing hidden domains in L2 that facilitate membrane penetration of L2/vDNA complexes. In an effort to identify regions of L2 that might physically interact with membranes, we have subjected the L2 sequence to multiple transmembrane (TM) domain prediction algorithms. Here, we describe a conserved TM domain within L2 (residues 45 to 67) and investigate its role in HPV16 infection. In vitro, the predicted TM domain adopts an alpha-helical structure in lipid environments and can function as a real TM domain, although not as efficiently as the bona fide TM domain of PDGFR. An L2 double point mutant renders the TM domain nonfunctional and blocks HPV16 infection by preventing endosomal translocation of vDNA. The TM domain contains three highly conserved GxxxG motifs. These motifs can facilitate homotypic and heterotypic interactions between TM helices, activities that may be important for vDNA translocation. Disruption of some of these GxxxG motifs resulted in noninfectious viruses, indicating a critical role in infection. Using a ToxR-based homo-oligomerization assay, we show a propensity for this TM domain to self-associate in a GxxxG-dependent manner. These data suggest an important role for the self-associating L2 TM domain and the conserved GxxxG motifs in the transfer of vDNA across the endo-/lysosomal membrane.

Statistical analyses and computational prediction of helical kinks in membrane proteins

We have carried out statistical analyses and computer simulations of helical kinks for TM helices in the PDBTM database. About 59 % of 1562 TM helices showed a significant kink, and 38 % of these kinks are associated with prolines in a range of ±4 residues. Our analyses show that helical kinks are more populated in the central region of helices, particularly in the range of 1-3 residues away from the helix center. Among 1,053 helical kinks analyzed, 88 % of kinks are bends (change in helix axis without loss of helical character) and 12 % are disruptions (change in helix axis and loss of helical character). It is found that proline residues tend to cause larger kink angles in helical bends, while this effect is not observed in helical disruptions. A further analysis of these kinked helices suggests that a kinked helix usually has 1-2 broken backbone hydrogen bonds with the corresponding N-O distance in the range of 4.2-8.7 Å, whose distribution is sharply peaked at 4.9 Å followed by an exponential decay with increasing distance. Our main aims of this study are to understand the formation of helical kinks and to predict their structural features. Therefore we further performed molecular dynamics (MD) simulations under four simulation scenarios to investigate kink formation in 37 kinked TM helices and 5 unkinked TM helices. The representative models of these kinked helices are predicted by a clustering algorithm, SPICKER, from numerous decoy structures possessing the above generic features of kinked helices. Our results show an accuracy of 95 % in predicting the kink position of kinked TM helices and an error less than 10° in the angle prediction of 71.4 % kinked helices. For unkinked helices, based on various structure similarity tests, our predicted models are highly consistent with their crystal structure. These results provide strong supports for the validity of our method in predicting the structure of TM helices.

Molecular dissection of TatC defines critical regions essential for protein transport and a TatB-TatC contact site

The twin arginine transport (Tat) system transports folded proteins across the prokaryotic cytoplasmic membrane and the plant thylakoid membrane. TatC is the largest and most conserved component of the Tat machinery. It forms a multisubunit complex with TatB and binds the signal peptides of Tat substrates. Here we have taken a random mutagenesis approach to identify substitutions in Escherichia coli TatC that inactivate protein transport. We identify 32 individual amino acid substitutions that abolish or severely compromise TatC activity. The majority of the inactivating substitutions fall within the first two periplasmic loops of TatC. These regions are predicted to have conserved secondary structure and results of extensive amino acid insertion and deletion mutagenesis are consistent with these conserved elements being essential for TatC function. Three inactivating substitutions were identified in the fifth transmembrane helix of TatC. The inactive M205R variant could be suppressed by mutations affecting amino acids in the transmembrane helix of TatB. A physical interaction between TatC helix 5 and the TatB transmembrane helix was confirmed by the formation of a site-specific disulphide bond between TatC M205C and TatB L9C variants. This is the first molecular contact site mapped to single amino acid level between these two proteins.

Shifting hydrogen bonds may produce flexible transmembrane helices

The intricate functions of membrane proteins would not be possible without bends or breaks that are remarkably common in transmembrane helices. The frequent helix distortions are nevertheless surprising because backbone hydrogen bonds should be strong in an apolar membrane, potentially rigidifying helices. It is therefore mysterious how distortions can be generated by the evolutionary currency of random point mutations. Here we show that we can engineer a transition between distinct distorted helix conformations in bacteriorhodopsin with a single-point mutation. Moreover, we estimate the energetic cost of the conformational transitions to be smaller than 1 kcal/mol. We propose that the low energy of distortion is explained in part by the shifting of backbone hydrogen bonding partners. Consistent with this view, extensive backbone hydrogen bond shifts occur during helix conformational changes that accompany functional cycles. Our results explain how evolution has been able to liberally exploit transmembrane helix bending for the optimization of membrane protein structure, function, and dynamics.

String kernels and high-quality data set for improved prediction of kinked helices in α-helical membrane proteins

The reasons for distortions from optimal α-helical geometry are widely unknown, but their influences on structural changes of proteins are significant. Hence, their prediction is a crucial problem in structural bioinformatics. For the particular case of kink prediction, we generated a data set of 132 membrane proteins containing 1014 manually labeled helices and examined the environment of kinks. Our sequence analysis confirms the great relevance of proline and reveals disproportionately high occurrences of glycine and serine at kink positions. The structural analysis shows significantly different solvent accessible surface area mean values for kinked and nonkinked helices. More important, we used this data set to validate string kernels for support vector machines as a new kink prediction method. Applying the new predictor, about 80% of all helices could be correctly predicted as kinked or nonkinked even when focusing on small helical fragments. The results exceed recently reported accuracies of alternative approaches and are a consequence of both the method and the data set.