Computational differentiation of N-terminal signal peptides and transmembrane helices

Signal peptides and transmembrane helices both contain a stretch of hydrophobic amino acids. This common feature makes it difficult for signal peptide and transmembrane helix predictors to correctly assign identity to stretches of hydrophobic residues near the N-terminal methionine of a protein sequence. The inability to reliably distinguish between N-terminal transmembrane helix and signal peptide is an error with serious consequences for the prediction of protein secretory status or transmembrane topology. In this study, we report a new method for differentiating protein N-terminal signal peptides and transmembrane helices. Based on the sequence features extracted from hydrophobic regions (amino acid frequency, hydrophobicity, and the start position), we set up discriminant functions and examined them on non-redundant datasets with jackknife tests. This method can incorporate other signal peptide prediction methods and achieve higher prediction accuracy. For Gram-negative bacterial proteins, 95.7% of N-terminal signal peptides and transmembrane helices can be correctly predicted (coefficient 0.90). Given a sensitivity of 90%, transmembrane helices can be identified from signal peptides with a precision of 99% (coefficient 0.92). For eukaryotic proteins, 94.2% of N-terminal signal peptides and transmembrane helices can be correctly predicted with coefficient 0.83. Given a sensitivity of 90%, transmembrane helices can be identified from signal peptides with a precision of 87% (coefficient 0.85). The method can be used to complement current transmembrane protein prediction and signal peptide prediction methods to improve their prediction accuracies. (C) 2003 Elsevier Inc. All rights reserved.

Light-induced exposure of the cytoplasmic end of transmembrane helix seven in rhodopsin

A key step in signal transduction in the visual cell is the light-induced conformational change of rhodopsin that triggers the binding and activation of the guanine nucleotide-binding protein. Site-directed mAbs against bovine rhodopsin were produced and used to detect and characterize these conformational changes upon light activation. Among several antibodies that bound exclusively to the light-activated state, an antibody (IgG subclass) with the highest affinity (Ka ≈ 6 × 10−9 M) was further purified and characterized. The epitope of this antibody was mapped to the amino acid sequence 304–311. This epitope extends from the central region to the cytoplasmic end of the seventh transmembrane helix and incorporates a part of a highly conserved NPXXY motif, a critical region for signaling and agonist-induced internalization of several biogenic amine and peptide receptors. In the dark state, no binding of the antibody to rhodopsin was detected. Accessibility of the epitope to the antibody correlated with formation of the metarhodopsin II photointermediate and was reduced significantly at the metarhodopsin III intermediate. Further, incubation of the antigen–antibody complex with 11-cis-retinal failed to regenerate the native rhodopsin chromophore. These results suggest significant and reversible conformational changes in close proximity to the cytoplasmic end of the seventh transmembrane helix of rhodopsin that might be important for folding and signaling.

Protonatable residues at the cytoplasmic end of transmembrane helix-2 in the signal transducer HtrI control photochemistry and function of sensory rhodopsin I.

Neutral residue replacements were made of 21 acidic and basic residues within the N-terminal half of the Halobacterium salinarium signal transducer HtrI [the halobacterial transducer for sensory rhodopsin I (SRI)] by site-specific mutagenesis. The replacements are all within the region of HtrI that we previously concluded from deletion analysis to contain sites of interaction with the phototaxis receptor SRI. Immunoblotting shows plasmid expression of the htrI-sopI operon containing the mutations produces SRI and mutant HtrI in cells at near wild-type levels. Six of the HtrI mutations perturb photochemical kinetics of SRI and one reverses the phototaxis response. Substitution with neutral amino acids of Asp-86, Glu-87, and Glu-108 accelerate, and of Arg-70, Arg-84, and Arg-99 retard, the SRI photocycle. Opposite effects on photocycle rate cancel in double mutants containing one replaced acidic and one replaced basic residue. Laser flash spectroscopy shows the kinetic perturbations are due to alteration of the rate of reprotonation of the retinylidene Schiff base. All of these mutations permit normal attractant and repellent signaling. On the other hand, the substitution of Glu-56 with the isosteric glutamine converts the normally attractant effect of orange light to a repellent signal in vivo at neutral pH (inverted signaling). Low pH corrects the inversion due to Glu-56 -> Gln and the apparent pK of the inversion is increased when arginine is substituted at position 56. The results indicate that the cytoplasmic end of transmembrane helix-2 and the initial part of the cytoplasmic domain contain interaction sites with SRI. To explain these and previous results, we propose a model in which (i) the HtrI region identified here forms part of an electrostatic bonding network that extends through the SRI protein and includes its photoactive site; (ii) alteration of this network by photoisomerization-induced Schiff base deprotonation and reprotonation shifts HtrI between attractant and repellent conformations; and (iii) HtrI mutations and extracellular pH alter the equilibrium ratios of these conformations.

Structure and function in rhodopsin: correct folding and misfolding in point mutants at and in proximity to the site of the retinitis pigmentosa mutation Leu-125-->Arg in the transmembrane helix C.

L125R is a mutation in the transmembrane helix C of rhodopsin that is associated with autosomal dominant retinitis pigmentosa. To probe the orientation of the helix and its packing in the transmembrane domain, we have prepared and studied the mutations E122R, I123R, A124R, S127R, L125F, and L125A at, and in proximity to, the above mutation site. Like L125R, the opsin expressed in COS-1 cells from E122R did not bind 11-cis-retinal, whereas those from I123R and S127R formed the rhodopsin chromophore partially. A124R opsin formed the rhodopsin chromophore (lambda max 495 nm) in the dark, but the metarhodopsin II formed on illumination decayed about 6.5 times faster than that of the wild type and was defective in transducin activation. The mutant opsins from L125F and L125A bound 11-cis-retinal only partially, and in both cases, the mixtures of the proteins produced were separated into retinal-binding and non-retinal-binding (misfolded) fractions. The purified mutant rhodopsin from L125F showed lambda max at 500 nm, whereas that from L125A showed lambda max at 503 nm. The mutant rhodopsin L125F showed abnormal bleaching behavior and both mutants on illumination showed destabilized metarhodopsin II species and reduced transducin activation. Because previous results have indicated that misfolding in rhodopsin is due to the formation of a disulfide bond other than the normal disulfide bond between Cys-110 and Cys-187 in the intradiscal domain, we conclude from the misfolding in mutants L125F and L125A that the folding in vivo in the transmembrane domain is coupled to that in the intradiscal domain.

SVMtm: Support vector machines to predict transmembrane segments

A new method has been developed for prediction of transmembrane helices using support vector machines. Different coding schemes of protein sequences were explored, and their performances were assessed by crossvalidation tests. The best performance method can predict the transmembrane helices with sensitivity of 93.4% and precision of 92.0%. For each predicted transmembrane segment, a score is given to show the strength of transmembrane signal and the prediction reliability. In particular, this method can distinguish transmembrane proteins from soluble proteins with an accuracy of similar to99%. This method can be used to complement current transmembrane helix prediction methods and can be Used for consensus analysis of entire proteomes . The predictor is located at http://genet.imb.uq.edu.au/predictors/ SVMtm. (C) 2004 Wiley Periodicals, Inc.

Molecular basis for reduced estrone sulfate transport and altered modulator sensitivity of transmembrane helix (TM) 6 and TM17 mutants of multidrug resistance protein 1 (ABCC1)

Multidrug resistance protein 1 (MRP1) confers drug resistance and also mediates cellular efflux of many organic anions. MRP1 also transports glutathione (GSH); furthermore, this tripeptide stimulates transport of several substrates, including estrone 3-sulfate. We have previously shown that mutations of Lys(332) in transmembrane helix (TM) 6 and Trp(1246) in TM17 cause different substrate-selective losses in MRP1 transport activity. Here we have extended our characterization of mutants K332L and W1246C to further define the different roles these two residues play in determining the substrate and inhibitor specificity of MRP1. Thus, we have shown that TM17-Trp(1246) is crucial for conferring drug resistance and for binding and transport of methotrexate, estradiol glucuronide, and estrone 3-sulfate, as well as for binding of the tricyclic isoxazole inhibitor N-[3-(9-chloro-3-methyl-4-oxo-4H-isoxazolo-[4,3-c]quinolin-5-yl)-cyclohexylmethyl]-benzamide (LY465803). In contrast, TM6-Lys(332) is important for enabling GSH and GSH-containing compounds to serve as substrates (e.g., leukotriene C(4)) or modulators (e.g., S-decyl-GSH, GSH disulfide) of MRP1 and, further, for enabling GSH (or S-methyl-GSH) to enhance the transport of estrone 3-sulfate and increase the inhibitory potency of LY465803. On the other hand, both mutants are as sensitive as wild-type MRP1 to the non-GSH-containing inhibitors (E)-3-[[[3-[2-(7-chloro-2-quinolinyl)ethenyl]phenyl][[3-(dimethylamino)-3-oxopropyl]thio]methyl]thio]-propanoic acid (MK571), 1-[2-hydroxy-3-propyl-4-[4-(1H-tetrazol-5-yl)butoxy]phenyl]-ethanone (LY171883), and highly potent 6-[4'-carboxyphenylthio]-5[S]-hydroxy-7[E], 11[Z]14[Z]-eicosatetrenoic acid (BAY u9773). Finally, the differing abilities of the cysteinyl leukotriene derivatives leukotriene C(4), D(4), and F(4) to inhibit estradiol glucuronide transport by wild-type and K332L mutant MRP1 provide further evidence that TM6-Lys(332) is involved in the recognition of the gamma-Glu portion of substrates and modulators containing GSH or GSH-like moieties.

Predicting the solvent accessibility of transmembrane residues from protein sequence

In this study, we propose a novel method to predict the solvent accessible surface areas of transmembrane residues. For both transmembrane alpha-helix and beta-barrel residues, the correlation coefficients between the predicted and observed accessible surface areas are around 0.65. On the basis of predicted accessible surface areas, residues exposed to the lipid environment or buried inside a protein can be identified by using certain cutoff thresholds. We have extensively examined our approach based on different definitions of accessible surface areas and a variety of sets of control parameters. Given that experimentally determining the structures of membrane proteins is very difficult and membrane proteins are actually abundant in nature, our approach is useful for theoretically modeling membrane protein tertiary structures, particularly for modeling the assembly of transmembrane domains. This approach can be used to annotate the membrane proteins in proteomes to provide extra structural and functional information.

Catechol-binding serines of beta(2)-adrenergic receptors control the equilibrium between active and inactive receptor states.

The binding free energy for the interaction between serines 204 and 207 of the fifth transmembrane helix of the beta(2)-adrenergic receptor (beta(2)-AR) and catecholic hydroxyl (OH) groups of adrenergic agonists was analyzed using double mutant cycles. Binding affinities for catecholic and noncatecholic agonists were measured in wild-type and mutant receptors, carrying alanine replacement of the two serines (S204A, S207A beta(2)-AR), a constitutive activating mutation, or both. The free energy coupling between the losses of binding energy attributable to OH deletion from the ligand and from the receptor indicates a strong interaction (nonadditivity) as expected for a direct binding between the two sets of groups. However, we also measured a significant interaction between the deletion of OH groups from the receptor and the constitutive activating mutation. This suggests that a fraction of the decrease in agonist affinity caused by serine mutagenesis may involve a shift in the conformational equilibrium of the receptor toward the inactive state. Direct measurements using a transient transfection assay confirm this prediction. The constitutive activity of the (S204A, S207A) beta(2)-AR mutant is 50 to 60% lower than that of the wild-type beta(2)-AR. We conclude that S204 and S207 do not only provide a docking site for the agonist, but also control the equilibrium of the receptor between active (R*) and inactive (R) forms.

Protein structure prediction servers at University College London

A number of state-of-the-art protein structure prediction servers have been developed by researchers working in the Bioinformatics Unit at University College London. The popular PSIPRED server allows users to perform secondary structure prediction, transmembrane topology prediction and protein fold recognition. More recent servers include DISOPRED for the prediction of protein dynamic disorder and DomPred for domain boundary prediction.

The PSIPRED protein structure prediction server

The PSIPRED protein structure prediction server allows users to submit a protein sequence, perform a prediction of their choice and receive the results of the prediction both textually via e-mail and graphically via the web. The user may select one of three prediction methods to apply to their sequence: PSIPRED, a highly accurate secondary structure prediction method; MEMSAT 2, a new version of a widely used transmembrane topology prediction method; or GenTHREADER, a sequence profile based fold recognition method.

A transmembrane form of annexin XII detected by site-directed spin labeling

Previous studies of the annexin family of Ca2+ binding proteins identified a soluble monomer in the absence of Ca2+ and a trimer adsorbed on the membrane surface in the presence of Ca2+. On the basis of site-directed spin-labeling studies of annexin XII at low pH, we now report a membrane-inserted form of the protein with a dramatically different structure. The data suggest that upon insertion a continuous transmembrane α-helix is reversibly formed from a helix–loop–helix motif in the solution structure. Other regions with similar membrane-insertion potential were identified in the amino acid sequence, and we propose that the corresponding helices come together to form an aqueous pore that mediates the ion channel activity reported for several annexins.

Polar residues drive association of polyleucine transmembrane helices

Although many polar residues are directly involved in transmembrane protein functions, the extent to which they contribute to more general structural features is still unclear. Previous studies have demonstrated that asparagine residues can drive transmembrane helix association through interhelical hydrogen bonding [Choma, C., Gratkowski, H., Lear, J. D. & DeGrado, W. F. (2000) Nat. Struct. Biol. 7, 161–166; and Zhou, F. X., Cocco, M. J., Russ, W. P., Brunger, A. T. & Engelman, D. M. (2000) Nat. Struct. Biol. 7, 154–160]. We have studied the ability of other polar residues to promote helix association in detergent micelles and in biological membranes. Our results show that polyleucine sequences with Asn, Asp, Gln, Glu, and His, residues capable of being simultaneously hydrogen bond donors and acceptors, form homo- or heterooligomers. In contrast, polyleucine sequences with Ser, Thr, and Tyr do not associate more than the polyleucine sequence alone. The results therefore provide experimental evidence that interactions between polar residues in the helices of transmembrane proteins may serve to provide structural stability and oligomerization specificity. Furthermore, such interactions can allow structural flexibility required for the function of some membrane proteins.

Photolabeling reveals the proximity of the alpha-neurotoxin binding site to the M2 helix of the ion channel in the nicotinic acetylcholine receptor.

A photoactivatable derivative of neurotoxin II from Naja naja oxiana containing a 125I-labeled p-azidosalicylamidoethyl-1,3'-dithiopropyl label at Lys-25 forms a photo-induced cross-link with the delta subunit of the membrane-bound Torpedo californica nicotinic acetylcholine receptor (AChR). The cross-linked radioactive receptor peptide was isolated by reverse-phase HPLC after tryptic digestion of the labeled delta subunit. The sequence of this peptide, delta-(260-277), and the position of the label at Ala-268 were established by matrix-assisted laser-desorption-ionization mass spectrometry based on the molecular mass and on post-source decay fragment analysis. With the known dimensions of the AChR molecule, of the photolabel, and of alpha-neurotoxin, finding the cross-link at delta Ala-268 (located in the upper part of the channel-forming transmembrane helix M2) means that the center of the alpha-neurotoxin binding site is situated at least approximately 40 A from the extracellular surface of the AChR, proximal to the channel axis.

Long-term treatment of familial male-limited precocious puberty (testotoxicosis) with cyproterone acetate or ketoconazole

Background Familial male-limited precocious puberty (FMPP) or testotoxicosis is a rare gonadotrophin-independent form of sexual precocity caused by constitutively activating mutations of the LH receptor. Several clinical therapeutic approaches have been reported for this disorder, but with a paucity of long-term outcome data. Objective To evaluate the long-term treatment of testotoxicosis with cyproterone acetate or ketoconazole. Design A multicentric retrospective clinical study. Patients Ten boys from eight unrelated Brazilian families who carried known LH-receptor activating mutations were treated with 70 mg/m(2) cyproterone acetate (n = 5) or 10 mg/kg ketoconazole (n = 5) for a mean period of 5 and 8 years, respectively. Measurements Chronological and bone ages, bone age/chronological age ratio, target height (TH) range, adult height, basal and GnRH-stimulated gonadotrophin levels and basal testosterone levels were assessed. Results Growth velocity decreased significantly during treatment with cyproterone acetate or ketoconazole when compared to pretreatment value in each group (P < 0.05). Bone age/chronological age ratio decreased significantly after cyproterone acetate or ketoconazole therapy. Basal testosterone levels were significantly lower in patients undergoing ketoconazole compared to cyproterone acetate treatment [0.6 +/- 0.3 nmol/l (42 +/- 21 ng/dl) vs. 5.6 +/- 4.0 nmol/l (392 +/- 280 ng/dl); P < 0.05], as expected. Secondary gonadotrophin-dependent precocious puberty occurred at a similar frequency (40%) in both groups. Five patients have attained adult height and two patients have already reached 90% of their adult height. Two of them achieved their TH range and one patient, for whom TH was not available, had an adult height of 0.3 SDS. Four boys (two in each group) did not attain their TH range. Conclusion Long-term treatment with cyproterone acetate or ketoconazole resulted in similar outcomes without important side-effects in boys with testotoxicosis. However, both therapies showed limited efficacy in attaining normal adult height.

A gating mutation in the internal pore of ASIC1a.

Using a substituted cysteine accessibility scan, we have investigated the structures that form the internal pore of the acid-sensing ion channel 1a. We have identified the amino acid residues Ala-22, Ile-33, and Phe-34 in the amino terminus and Arg-43 in the first transmembrane helix, which when mutated into cysteine, were modified by intracellular application of MTSET, resulting in channel inhibition. The inhibition of the R43C mutant by internal MTSET requires opening of the channel. In addition, binding of Cd2+ ions to R43C slows the channel inactivation. This indicates that the first transmembrane helix undergoes conformational changes during channel inactivation. The effect of Cd2+ on R43C can be obtained with Cd2+ applied at either the extracellular or the intracellular side, indicating that R43C is located in the channel pore. The block of the A22C, I33C, and F34C mutants by MTSET suggests that these residues in the amino terminus of the channel also participate to the internal pore.