Prediction of the translocon-mediated membrane insertion free energies of protein sequences

In silico evaluation of the influence of the translocon on partitioning of membrane segments

Background

The locations of the TM segments inside the membrane proteins are the consequence of a cascade of several events: the localizing of the nascent chain to the membrane, its insertion through the translocon, and the conformation adopted to reach its stable state inside the lipid bilayer. Even though the hydrophobic h-region of signal peptides and a typical TM segment are both composed of mostly hydrophobic side chains, the translocon has the ability to determine whether a given segment is to be inserted into the membrane. Our goal is to acquire robust biological insights into the influence of the translocon on membrane insertion of helices, obtained from the in silico discrimination between signal peptides and transmembrane segments of bitopic proteins. Therefore, by exploiting this subtle difference, we produce an optimized scale that evaluates the tendency of each amino acid to form sequences destined for membrane insertion by the translocon.

Results

The learning phase of our approach is conducted on carefully chosen data and easily converges on an optimal solution called the PMIscale (Potential Membrane Insertion scale). Our study leads to two striking results. Firstly, with a very simple sliding-window prediction method, PMIscale enables an efficient discrimination between signal peptides and signal anchors. Secondly, PMIscale is also able to identify TM segments and to localize them within protein sequences.

Conclusions

Despite its simplicity, the localization method based on PMIscale nearly attains the highest level of TM topography prediction accuracy as the current state-of-the-art prediction methods. These observations confirm the prominent role of the translocon in the localization of TM segments and suggest several biological hypotheses about the physical properties of the translocon.

Protein translocation across the ER membrane

Protein translocation into the endoplasmic reticulum (ER) is the first and decisive step in the biogenesis of most extracellular and many soluble organelle proteins in eukaryotic cells. It is mechanistically related to protein export from eubacteria and archaea and to the integration of newly synthesized membrane proteins into the ER membrane and the plasma membranes of eubacteria and archaea (with the exception of tail anchored membrane proteins). Typically, protein translocation into the ER involves cleavable amino terminal signal peptides in precursor proteins and sophisticated transport machinery components in the cytosol, the ER membrane, and the ER lumen. Depending on the hydrophobicity and/or overall amino acid content of the precursor protein, transport can occur co- or posttranslationally. The respective mechanism determines the requirements for certain cytosolic transport components. The two mechanisms merge at the level of the ER membrane, specifically, at the heterotrimeric Sec61 complex present in the membrane. The Sec61 complex provides a signal peptide recognition site and forms a polypeptide conducting channel. Apparently, the Sec61 complex is gated by various ligands, such as signal peptides of the transport substrates, ribosomes (in cotranslational transport), and the ER lumenal molecular chaperone, BiP. Binding of BiP to the incoming polypeptide contributes to efficiency and unidirectionality of transport. Recent insights into the structure of the Sec61 complex and the comparison of the transport mechanisms and machineries in the yeast Saccharomyces cerevisiae, the human parasite Trypanosoma brucei, and mammals have various important mechanistic as well as potential medical implications. This article is part of a Special Issue entitled Protein translocation across or insertion into membranes.

Contribution of charged and polar residues for the formation of the E1-E2 heterodimer from Hepatitis C Virus

The transmembrane domains of the envelope glycoprotein E1 and E2 have crucial multifunctional roles in the biogenesis of hepatitis C virus. We have performed molecular dynamics simulations to investigate a structural model of the transmembrane segments of the E1-E2 heterodimer. The simulations support the key role of the Lys370-Asp728 ion pair for mediating the E1-E2 heterodimerization. In comparison to these two residues, the simulation results also reveal the differential effect of the conserved Arg730 residue that has been observed in experimental studies. Furthermore, we discovered the formation of inter-helical hydrogen bonds via Asn367 that stabilize dimer formation. Simulations of single and double mutants further demonstrate the importance of the ion-pair and polar interactions between the interacting helix monomers. The conformation of the E1 fragment in the simulation of the E1-E2 heterodimer is in close agreement with an NMR structure of the E1 transmembrane segment. The proposed model of the E1-E2 heterodimer supports the postulated cooperative insertion of both helices by the translocon complex into the bilayer.

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

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

MINS2: revisiting the molecular code for transmembrane-helix recognition by the Sec61 translocon

To be fully functional, membrane proteins should not only fold, but also get inserted into the membrane, which is mediated by the Sec61 translocon. Recent experimental studies have attempted to elucidate how the Sec61 translocon accomplishes this delicate task by measuring the translocon-mediated membrane insertion free energies of 357 systematically designed peptides. On the basis of this data set, we have developed MINS2, a novel sequence-based computational method for predicting the membrane insertion free energies of protein sequences. A benchmark analysis of MINS2 shows that MINS2 signi.cantly outperforms previously proposed methods. Importantly, the application of MINS2 to known membrane protein structures shows that a better prediction of membrane insertion free energies does not lead to a better prediction of transmembrane segments of polytopic membrane proteins.

AVAILABILITY

A web server for MINS2 is publicly available at http://service.bioinformatik.uni-saarland.de/mins.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.