A general method for determining helix packing in membrane proteins in situ: Helices I and II are close to helix VII in the lactose permease of Escherichia coli

It was previously shown that coexpression of the lactose permease of Escherichia coli in two contiguous fragments leads to functional complementation. We demonstrate here that site-directed thiol crosslinking of coexpressed permease fragments can be used to determine helix proximity in situ without the necessity of purifying the permease. After coexpression of the six N-terminal (N6) and six C-terminal (C6) transmembrane helices, each with a single Cys residue, crosslinking was carried out in native membranes and assessed by the mobility of anti-C-terminal-reactive polypeptides on immunoblots. A Cys residue at position 242 or 245 (helix VII) forms a disulfide with a Cys residue at either position 28 or 29 (helix I), but not with a Cys residue at position 27, which is on the opposite face of helix I, thereby indicating that the face of helix I containing Pro-28 and Phe-29 is close to helix VII. Similarly, a Cys residue at position 242 or 245 (helix VII) forms a disulfide with a Cys residue at either position 52 or 53 (helix II), but not with a Cys residue at position 54. Furthermore, low-efficiency crosslinking is observed between a Cys residue at position 52 or 53 and a Cys residue at position 361 (helix XI). The results indicate that helix VII lies in close proximity to both helices I and II and that helix II is also close to helix XI. The method should be applicable to a number of different polytopic membrane proteins.

HUNT: launch of a full-length cDNA database from the Helix Research Institute

The Helix Research Institute (HRI) in Japan is releasing 4356 HUman Novel Transcripts and related information in the newly established HUNT database. The institute is a joint research project principally funded by the Japanese Ministry of International Trade and Industry, and the clones were sequenced in the governmental New Energy and Industrial Technology Development Organization (NEDO) Human cDNA Sequencing Project. The HUNT database contains an extensive amount of annotation from advanced analysis and represents an essential bioinformatics contribution towards understanding of the gene function. The HRI human cDNA clones were obtained from full-length enriched cDNA libraries constructed with the oligo-capping method and have resulted in novel full-length cDNA sequences. A large fraction has little similarity to any proteins of known function and to obtain clues about possible function we have developed original analysis procedures. Any putative function deduced here can be validated or refuted by complementary analysis results. The user can also extract information from specific categories like PROSITE patterns, PFAM domains, PSORT localization, transmembrane helices and clones with GENIUS structure assignments. The HUNT database can be accessed at http://www.hri.co.jp/HUNT.

BAliBASE (Benchmark Alignment dataBASE): enhancements for repeats, transmembrane sequences and circular permutations

BAliBASE is specifically designed to serve as an evaluation resource to address all the problems encountered when aligning complete sequences. The database contains high quality, manually constructed multiple sequence alignments together with detailed annotations. The alignments are all based on three-dimensional structural superpositions, with the exception of the transmembrane sequences. The first release provided sets of reference alignments dealing with the problems of high variability, unequal repartition and large N/C-terminal extensions and internal insertions. Here we describe version 2.0 of the database, which incorporates three new reference sets of alignments containing structural repeats, trans­membrane sequences and circular permutations to evaluate the accuracy of detection/prediction and alignment of these complex sequences. BAliBASE can be viewed at the web site http://www-igbmc.u-strasbg.fr/BioInfo/BAliBASE2/index.html or can be downloaded from ftp://ftp-igbmc.u-strasbg.fr/pub/BAliBASE2/.

Detecting the conformational change of transmembrane signaling in a bacterial chemoreceptor by measuring effects on disulfide cross-linking in vivo.

Transmembrane signaling by bacterial chemoreceptors is thought to involve relative movement among the four transmembrane helices of the homodimer. We assayed that movement by measuring effects of ligand occupancy on rates of oxidative cross-linking between cysteines introduced into neighboring helices of the transmembrane domain of chemoreceptor Trg from Escherichia coli. Measurements were done on chemoreceptors in their native environment, intact cells that were motile and chemotactically responsive. Receptor occupancy did not appear to cause drastic rearrangement of the four-helix structure since, among 67 cysteine pairs tested, the same 19 exhibited oxidative cross-linking in the presence or absence of saturating chemoattractant. However, occupancy did cause subtle changes that were detected as effects on rates of cross-linking. Among the seven disulfides appropriate for measurements of initial rates of formation, ligand occupancy had significant and different effects on all three cross-links that connected the two helices within a subunit but had minimal effects on the four that spanned the packing interface between subunits. This constitutes direct evidence that the conformational change of transmembrane signaling involves significant movement within a subunit and minimal movement between subunits, a pattern deduced from several previous studies and now documented directly. Among possible modes of movement between the two helices of a subunit, axial sliding of one helix relative to the other was the conformational change that best accounted for the observed effects on cross-linking.

Prediction of the stability of DNA triplexes.

We present rules that allow one to predict the stability of DNA pyrimidine.purine.pyrimidine (Y.R.Y) triple helices on the basis of the sequence. The rules were derived from van't Hoff analysis of 23 oligonucleotide triplexes tested at a variety of pH values. To predict the enthalpy of triplex formation (delta H degrees), a simple nearest-neighbor model was found to be sufficient. However, to accurately predict the free energy of the triplex (delta G degrees), a combination model consisting of five parameters was needed. These parameters were (i) the delta G degrees for helix initiation, (ii) the delta G degrees for adding a T-A.T triple, (iii) the delta G degrees for adding a C(+)-G.C triple, (iv) the penalty for adjacent C bases, and (v) the pH dependence of the C(+)-G.C triple's stability. The fitted parameters are highly consistent with thermodynamic data from the basis set, generally predicting both delta H degrees and delta G degrees to within the experimental error. Examination of the parameters points out several interesting features. The combination model predicts that C(+) -G.C. triples are much more stabilizing than T-A.T triples below pH 7.0 and that the stability of the former increases approximately equal to 1 kcal/mol per pH unit as the pH is decreased. Surprisingly though, the most stable sequence is predicted to be a CT repeat, as adjacent C bases partially cancel the stability of one another. The parameters successfully predict tm values from other laboratories, with some interesting exceptions.

Molecular mechanism of transmembrane signaling by the aspartate receptor: a model.

The aspartate receptor of bacterial chemotaxis is representative of a large class of membrane-spanning receptors found in prokaryotic and eukaryotic organisms. These receptors, which regulate histidine kinase pathways and possess two putative transmembrane helices per subunit, appear to control a wide variety of cellular processes. The best characterized subgroup of the two-helix receptor class is the homologous family of chemosensory receptors from Escherichia coli and Salmonella typhimurium, including the aspartate receptor. This receptor binds aspartate, an attractant, in the periplasmic compartment and undergoes an intramolecular, transmembrane conformational change, thereby modulating the autophosphorylation rate of a bound histidine kinase in the cytoplasm. Here, we analyze recent results from x-ray crystallographic, solution 19F NMR, and engineered disulfide studies probing the aspartate-induced structural change within the periplasmic and transmembrane regions of the receptor. Together, these approaches provide evidence that aspartate binding triggers a "swinging-piston" displacement of the second membrane-spanning helix, which is proposed to communicate the signal across the bilayer.

Three distinct structural environments of a transmembrane domain in the inwardly rectifying potassium channel ROMK1 defined by perturbation.

To probe the protein environment of an ion channel, we have perturbed the structure of a transmembrane domain by substituting side chains with those of two different sizes by using site-specific mutagenesis. We have used Trp and Ala as a high- and a low-impact perturbation probe, respectively, to replace each of 18 consecutive residues within the putative second transmembrane segment, M2, of an inwardly rectifying potassium channel, ROMK1. Our rationale is that a change in the channel function as a consequence of these mutations at a particular position will reflect the structural environment of the altered side chain. Each position can then be assigned to one of three classes of environments, as grated by different levels of perturbation: very tolerant (channel functions with both Trp and Ala substitutions), tolerant (function preserved with Ala but not with Trp substitution), and intolerant (either Ala or Trp substitution destroys function). We identify the very tolerant environment as being lipid-facing, tolerant as protein-interior-facing, and intolerant as pore-facing. We observe a strikingly ordered pattern of perturbation of all three environmental classes. This result indicates that M2 is a straight alpha-helix.

Transmembrane signaling characterized in bacterial chemoreceptors by using sulfhydryl cross-linking in vivo.

Transmembrane signaling by bacterial chemoreceptors is thought to involve conformational changes within a stable homodimer. We investigated the functional consequences of constraining movement between pairs of helices in the four-helix structure of the transmembrane domain of chemoreceptor Trg. Using a family of cysteine-containing receptors, we identified oxidation treatments for intact cells that catalyzed essentially complete sulfhydryl cross-linking at selected positions and yet left flagellar and sensory functions largely unperturbed. Constraining movement by cross-links between subunits had little effect on tactic response, but constraining movement between transmembrane segments of the monomer drastically reduced function. We deduce that transmembrane signaling requires substantial movement between transmembrane helices of a monomer but not between interacting helices across the interface between subunits.

Expression and biochemical analysis of the entire HIV-2 gp41 ectodomain: determinants of stability map to N- and C-terminal sequences outside the 6-helix bundle core

The folding of HIV gp41 into a 6-helix bundle drives virus-cell membrane fusion. To examine the structural relationship between the 6-helix bundle core domain and other regions of gp41, we expressed in Escherichia coli, the entire ectodomain of HIV-2(ST) gp41 as a soluble, trimeric maltose-binding protein (MBP)/gp41 chimera. Limiting proteolysis indicated that the Cys-591-Cys-597 disulfide-bonded region is outside a core domain comprising two peptides, Thr-529-Trp-589 and Val-604-Ser-666. A biochemical examination of MBP/gp41 chimeras encompassing these core peptides; indicated that the N-terminal polar segment, 521-528, and C-terminal membrane-proximal segment, 658-666, cooperate in stabilizing the ectodomain. A functional interaction between sequences outside the gp41 core may contribute energy to membrane fusion. (C) 2004 Published by Elsevier B.V. on behalf of the Federation of European Biochemical Societies.

A key role for transmembrane prolines in calcitonin receptor-like receptor agonist binding and signalling:implications for family B G-protein-coupled receptors

Calcitonin receptor like-receptor is a family B G-protein coupled receptor (GPCR). It requires receptor activity modifying protein (RAMP) 1 to give a calcitonin gene-related peptide (CGRP) receptor. Little is known of how members of this receptor family function. Proline residues often form important kinks in alpha-helices. Therefore, all proline residues within the transmembrane helices of the receptor (Pro241, Pro244 in helix 4, Pro275 in helix 5, Pro321 and Pro331 in helix 6) were mutated to alanine. Pro241 Pro275, and Pro321 are highly conserved throughout all family B GPCRs. The binding of CGRP and its ability to stimulate cAMP production were investigated in mutant and wild-type receptors after transient transfection into COS-7 cells with RAMP1. The P321A mutation significantly decreased the pEC(50) for CGRP and reduced its affinity but did not change cell-surface expression. Antagonist binding [CGRP(8-37) and 1-piperidinecarboxamide N-[2-[[5amino-1-[[4-(4-pyridinyl)-1-piperazinyl]carbonyl]pentyl]amino]-1-[(3 5-dibromo-4-hydroxyphenyl)methyl]-2-oxoethyl]-4-(1,4-dihydro-2-oxo-3(2H)-quina zolinyl) (BIBN4096BS)] was little altered by the mutation. Adrenomedullin-mediated signaling was disrupted when P321A was coexpressed with RAMP1, RAMP2, or RAMP3. The P331A mutant produced a moderate reduction in CGRP binding and receptor activation. Mutation of the other residues had no effect on receptor function. Thus, Pro321 and Pro331 are required for agonist binding and receptor activation. Modeling suggested that Pro321 induces a bend in helix 6, bringing its C terminus near that of helix 3, as seen in many family A GPCRs. This is abolished in P321A. P321A-I325P predicted to restore this conformation, showed wild-type activation. Modeling can also rationalize the effects of transmembrane proline mutants previously reported for another family B GPCR, the VPAC(1) receptor.

Beta barrel trans-membrane proteins:enhanced prediction using a Bayesian approach

Membrane proteins, which constitute approximately 20% of most genomes, form two main classes: alpha helical and beta barrel transmembrane proteins. Using methods based on Bayesian Networks, a powerful approach for statistical inference, we have sought to address beta-barrel topology prediction. The beta-barrel topology predictor reports individual strand accuracies of 88.6%. The method outlined here represents a potentially important advance in the computational determination of membrane protein topology.

Prediction of peptide and protein propensity for amyloid formation

This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Machine learning tools for protein annotation: the cases of transmembrane β-barrel and myristoylated proteins

Biology is now a “Big Data Science” thanks to technological advancements allowing the characterization of the whole macromolecular content of a cell or a collection of cells. This opens interesting perspectives, but only a small portion of this data may be experimentally characterized. From this derives the demand of accurate and efficient computational tools for automatic annotation of biological molecules. This is even more true when dealing with membrane proteins, on which my research project is focused leading to the development of two machine learning-based methods: BetAware-Deep and SVMyr. BetAware-Deep is a tool for the detection and topology prediction of transmembrane beta-barrel proteins found in Gram-negative bacteria. These proteins are involved in many biological processes and primary candidates as drug targets. BetAware-Deep exploits the combination of a deep learning framework (bidirectional long short-term memory) and a probabilistic graphical model (grammatical-restrained hidden conditional random field). Moreover, it introduced a modified formulation of the hydrophobic moment, designed to include the evolutionary information. BetAware-Deep outperformed all the available methods in topology prediction and reported high scores in the detection task. Glycine myristoylation in Eukaryotes is the binding of a myristic acid on an N-terminal glycine. SVMyr is a fast method based on support vector machines designed to predict this modification in dataset of proteomic scale. It uses as input octapeptides and exploits computational scores derived from experimental examples and mean physicochemical features. SVMyr outperformed all the available methods for co-translational myristoylation prediction. In addition, it allows (as a unique feature) the prediction of post-translational myristoylation. Both the tools here described are designed having in mind best practices for the development of machine learning-based tools outlined by the bioinformatics community. Moreover, they are made available via user-friendly web servers. All this make them valuable tools for filling the gap between sequential and annotated data.

Exploring Iris Colour Prediction And Ancestry Inference In Admixed Populations Of South America.

New DNA-based predictive tests for physical characteristics and inference of ancestry are highly informative tools that are being increasingly used in forensic genetic analysis. Two eye colour prediction models: a Bayesian classifier - Snipper and a multinomial logistic regression (MLR) system for the Irisplex assay, have been described for the analysis of unadmixed European populations. Since multiple SNPs in combination contribute in varying degrees to eye colour predictability in Europeans, it is likely that these predictive tests will perform in different ways amongst admixed populations that have European co-ancestry, compared to unadmixed Europeans. In this study we examined 99 individuals from two admixed South American populations comparing eye colour versus ancestry in order to reveal a direct correlation of light eye colour phenotypes with European co-ancestry in admixed individuals. Additionally, eye colour prediction following six prediction models, using varying numbers of SNPs and based on Snipper and MLR, were applied to the study populations. Furthermore, patterns of eye colour prediction have been inferred for a set of publicly available admixed and globally distributed populations from the HGDP-CEPH panel and 1000 Genomes databases with a special emphasis on admixed American populations similar to those of the study samples.

Petroleomics By Electrospray Ionization Ft-icr Mass Spectrometry Coupled To Partial Least Squares With Variable Selection Methods: Prediction Of The Total Acid Number Of Crude Oils.

Negative-ion mode electrospray ionization, ESI(-), with Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR MS) was coupled to a Partial Least Squares (PLS) regression and variable selection methods to estimate the total acid number (TAN) of Brazilian crude oil samples. Generally, ESI(-)-FT-ICR mass spectra present a power of resolution of ca. 500,000 and a mass accuracy less than 1 ppm, producing a data matrix containing over 5700 variables per sample. These variables correspond to heteroatom-containing species detected as deprotonated molecules, [M - H](-) ions, which are identified primarily as naphthenic acids, phenols and carbazole analog species. The TAN values for all samples ranged from 0.06 to 3.61 mg of KOH g(-1). To facilitate the spectral interpretation, three methods of variable selection were studied: variable importance in the projection (VIP), interval partial least squares (iPLS) and elimination of uninformative variables (UVE). The UVE method seems to be more appropriate for selecting important variables, reducing the dimension of the variables to 183 and producing a root mean square error of prediction of 0.32 mg of KOH g(-1). By reducing the size of the data, it was possible to relate the selected variables with their corresponding molecular formulas, thus identifying the main chemical species responsible for the TAN values.