Prediction of disulfide connectivity in proteins.

MOTIVATION: A major problem in protein structure prediction is the correct location of disulfide bridges in cysteine-rich proteins. In protein-folding prediction, the location of disulfide bridges can strongly reduce the search in the conformational space. Therefore the correct prediction of the disulfide connectivity starting from the protein residue sequence may also help in predicting its 3D structure. RESULTS: In this paper we equate the problem of predicting the disulfide connectivity in proteins to a problem of finding the graph matching with the maximum weight. The graph vertices are the residues of cysteine-forming disulfide bridges, and the weight edges are contact potentials. In order to solve this problem we develop and test different residue contact potentials. The best performing one, based on the Edmonds-Gabow algorithm and Monte-Carlo simulated annealing reaches an accuracy significantly higher than that obtained with a general mean force contact potential. Significantly, in the case of proteins with four disulfide bonds in the structure, the accuracy is 17 times higher than that of a random predictor. The method presented here can be used to locate putative disulfide bridges in protein-folding. AVAILABILITY: The program is available upon request from the authors. CONTACT: Casadio@alma.unibo.it; Piero@biocomp.unibo.it.     

Stereoselective disulfide formation stabilizes the local peptide conformation in nisin mimics

Nisin is a polymacrocyclic peptide antimicrobial with high activity against Gram-positive bacteria. Lanthionine and methyllanthionine bridges, closing the macrocycles, are stabilized by thioether bonds, formed between cysteines and dehydrated serine or threonine. The role of polypeptide backbone conformation in the formation of macrocycles A and B within cysteine mutants of nisin residues 1?12 is investigated here by molecular dynamics simulations. Enantiomeric combinational space of Cys3 and Cys7 and of Cys8 and Cys11 is examined for the preference of disulfide bond formation over helical turn formation within this region. A clear preference for spontaneous disulfide formation and closure of rings 3,7 and 8,11 is demonstrated for the D-Cys3, D-Cys7, L-Cys8, L-Cys11 nisin homologue, while interlinked rings A and B are obtained through disulfide bridges between L-Cys3 and D-Cys8 and between D-Cys7 and D-Cys11. This study offers a simple designer approach to solid phase synthesis of macrocyclic peptides and lantibiotic analogues.     

A method for the characterization of foldings in protein ribbon models

The ribbon model of chain macromolecules is a useful tool for analyzing some of the targe-scale shape features of these complex systems. Up to now, the ribbon model has been used mostly to produce graphical displays, which are usually analyzed by visual inspection. In this work we suggest a computational method for characterizing automatically, in a concise and algebraic fashion, some of the important shape features of these ribbon models. The procedure is based on a graph-theoretical and knot-theoretical characterization of three well-defined projections of a space curve associated with the ribbon. The labeled graphs can be characterized by the handedness of the crossovers in the ribbon that are the vertices of the graph. The method can be used to provide a fully algebraic representation of the changes occurring when a molecule, such as a protein, undergoes conformational rearrangements (folding), as well as to provide a shape comparison for a pair of related molecular ribbons. This algebraic representation is well suited for easy storage, retrieval, and computer manipulation of the information on the ribbon's shape. Illustrative examples of the method are provided.     

Topology of globular proteins

This paper inquires whether it is reasonable to expect the native structure of proteins to be “knotted”. To this end, some topological properties of polypeptides containing disulfide bridges are discussed using notions from mathematical knot theory and graph theory. The probability of occurrence of knots in random cyclic polymers is calculated as a function of chain length by elementary Monte Carlo methods. The implications of this for protein renaturation and for determining the tertiary structure of proteins are discussed.     

Ribosomal protein L7/L12 cross-links to proteins in separate regions of the 50 S ribosomal subunit of Escherichia coli

The 50 S ribosomal subunits from Escherichia coli were modified by reaction with 2-iminothiolane under conditions in which 65 sulfhydryl groups, about 2/protein, were added per subunit. Earlier work showed that protein L7/L12 was modified more extensively than the average but that nearly all 50 S proteins contained sulfhydryl groups. Mild oxidation led to the formation of disulfide protein-protein cross-links. These were fractionated by urea gel electrophoresis and then analyzed by diagonal gel electrophoresis. Cross-linked complexes containing two, three, and possibly four copies of L7/L12 were evident. Cross-links between L7/L12 and other ribosomal proteins were also formed. These proteins were identified as L5, L6, L10, L11, and, in lower yield, L9, L14, and L17. The yields of cross-links to L5, L6, L10, and L11 were comparable to the most abundant cross-links formed. Similar experiments were performed with 70 S ribosomes. Protein L7/L12 in 70 S ribosomes was cross-linked to proteins L6, L10, and L11. The strong L7/L12-L5 cross-link found in 50 S subunits was absent in 70 S ribosomes. No cross-links between 30 S proteins and L7/L12 were observed.     

Contributions of basic residues to ribosomal protein L11 recognition of RNA

The C-terminal domain of ribosomal protein L11, L11-C76, binds in the distorted minor groove of a helix within a 58 nucleotide domain of 23 S rRNA. To study the electrostatic component of RNA recognition in this protein, arginine and lysine residues have been individually mutated to alanine or methionine residues at the nine sequence positions that are conserved as basic residues among bacterial L11 homologs. In measurements of the salt dependence of RNA-binding, five of these mutants have a reduced value of - partial differentiallog(K(obs))/ partial differentiallog[KCl] as compared to the parent L11-C76 sequence, indicating that these residues interact with the RNA electrostatic field. These five residues are located at the perimeter of the RNA-binding surface of the protein; all five of them form salt bridges with phosphates in the crystal structure of the complex. A sixth residue, Lys47, was found to make an electrostatic contribution to binding when measurements were made at pH 6.0, but not at pH 7.0; its pK in the free protein must be <6.5. The unusual behavior of Lys47 is explained by its burial in the hydrophobic core of the free protein, and unburial in the RNA-bound protein, where it forms a salt bridge with a phosphate. The contributions of these six residues to the electrostatic component of binding are not additive; thus the magnitude of the salt dependence cannot be used to count the number of ionic interactions in this complex. By interacting with irregular features of the RNA backbone, including an S-turn, these basic residues contribute to the specificity of L11 for its target site.     

Ribosomal protein L7/L12 has a helix-turn-helix motif similar to that found in DNA-binding regulatory proteins.

Inspection of the structure of the C-terminal domain of ribosomal protein L7/L12 (1) reveals a helix-turn-helix motif similar to the one found in many DNA-binding regulatory proteins (2-5). The 19 alpha-carbon atoms of the L7/L12 alpha-helices superimpose on the DNA binding helices of CAP and cro with root-mean-square distances between corresponding alpha carbons of 1.45 and 1.55 A, respectively. These helices in L7/L12 are within a patch of highly conserved residues on the surface of L7/L12 whose role is as yet uncertain. We raise the possibility that they may constitute a binding site for nucleic acids, most probably RNA. Consistent with this hypothesis are calculations of the electrostatic charge potential surrounding the protein, which show a region of positive potential centered on the first of these helices.     

Homologous proteins with different folds: the three-dimensional structures of domains 1 and 6 of the multiple Kazal-type inhibitor LEKTI

We have determined the solution structures of recombinant domain 1 and native domain 6 of the multi-domain Kazal-type serine proteinase inhibitor LEKTI using multi-dimensional NMR spectroscopy. While two of the 15 potential inhibitory LEKTI domains contain three disulfide bonds typical of Kazal-type inhibitors, the remaining 13 domains have only two of these disulfide bridges. Therefore, they may represent a novel type of serine proteinase inhibitor. The first and the sixth LEKTI domain, which have been isolated from human blood ultrafiltrate, belong to this group. In spite of sharing the same disulfide pattern and a sequence identity of about 35% from the first to the fourth cysteine, the two proteins show different structures in this region. The three-dimensional structure of domain 6 consists of two helices and a beta-hairpin structure, and closely resembles the three-dimensional fold of classical Kazal-type serine proteinase inhibitors including the inhibitory binding loop. Domain 6 has been shown to be an efficient, but non-permanent serine proteinase inhibitor. The backbone geometry of its canonical loop is not as well defined as the remaining structural elements, providing a possible explanation for its non-permanent inhibitory activity. We conclude that domain 6 belongs to a subfamily of classical Kazal-type inhibitors, as the third disulfide bond and a third beta-strand are missing. The three-dimensional structure of domain 1 shows three helices and a beta-hairpin, but the central part of the structure differs remarkably from that of domain 6. The sequence adopting hairpin structure in domain 6 exhibits helical conformation in domain 1, and none of the residues within the putative P3 to P3' stretch features backbone angles that resemble those of the canonical loop of known proteinase inhibitors. No proteinase has been found to be inhibited by domain 1. We conclude that domain 1 adopts a new protein fold and is no canonical serine proteinase inhibitor.     

1H 2D NMR and distance geometry study of the folding of Ecballium elaterium trypsin inhibitor, a member of the squash inhibitors family

The solution conformation of synthetic Ecballium elaterium trypsin inhibitor II, a 28-residue peptide with 3 disulfide bridges, has been studied by 1H 2D NMR measurements. Secondary structure elements were determined: a miniantiparallel beta-sheet Met 7-Cys 9 and Gly 25-Cys 27, a beta-hairpin 20-28 with beta-turn 22-25, and two tight turns Asp 12-Cys 15 and Leu 16-Cys 19. A set of interproton distance restraints deduced from two-dimensional nuclear Overhauser enhancement spectra and 13 phi backbone torsion angles restraints were used as the basis of three-dimensional structure computations including disulfide bridges arrangement by using distance geometry calculations. Computations for the 15 possible S-S linkage combinations lead to the proposal of the array 2-19, 9-21, 15-27 as the most probable structure for EETI II.     

The unfolding pathway and conformational stability of potato carboxypeptidase inhibitor

The unfolding and denaturation curves of potato carboxypeptidase inhibitor (PCI) were investigated using the technique of disulfide scrambling. In the presence of denaturant and thiol initiator, the native PCI denatures by shuffling its native disulfide bonds and converts to form a mixture of scrambled PCI that consists of 9 out of a possible 14 isomers. The denaturation curve is determined by the fraction of native PCI converted to scrambled isomers under increasing concentrations of denaturant. The concentration of guanidine thiocyanate, guanidine hydrochloride, and urea required to denature 50% of the native PCI was found to be 0.7, 1.45, and 8 m, respectively. The PCI unfolding curve was constructed through the analysis of structures of scrambled isomers that were denatured under increasing concentrations of denaturant. These results reveal the existence of structurally defined unfolding intermediates and a progressive expansion of the polypeptide chain. The yield of the beads-form isomer (Cys(8)-Cys(12), Cys(18)-Cys(24), and Cys(27)-Cys(34)) as a fraction of total denatured PCI was shown to be directly proportional to the strength of the denaturing condition. Furthermore, the PCI sequence was unable to fold quantitatively into a single native structure. Under physiological conditions, the scrambled isomers of PCI that constitute about 4% of the protein were in equilibrium with native PCI.     

Structure of the C-terminal domain of the ribosomal protein L7/L12 from Escherichia coli at 1.7 A

The structure of a C-terminal fragment of the ribosomal protein L7/L12 from Escherichia coli has been refined using crystallographic data to 1.7 A resolution. The R-value is 17.4%. Six residues at the N terminus are too disordered in the structure to be localized. These residues are probably part of a hinge in the complete L7/L12 molecule. The possibility that a 2-fold crystallographic axis is a molecular 2-fold axis is discussed. A patch of invariant residues on the surface of the dimer is probably involved in functional interactions with elongation factors.     

NMR solution structure of alpha-conotoxin ImI and comparison to other conotoxins specific for neuronal nicotinic acetylcholine receptors

Alpha-Conotoxins, peptides produced by predatory species of Conus marine snails, are potent antagonists of nicotinic acetylcholine receptors (nAChRs), ligand-gated ion channels involved in synaptic transmission. We determined the NMR solution structure of the smallest known alpha-conotoxin, ImI, a 12 amino acid peptide that binds specifically to neuronal alpha7-containing nAChRs in mammals. Calculation of the structure was based on a total of 80 upper distance constraints and 31 dihedral angle constraints resulting in 20 representative conformers with an average pairwise rmsd of 0.44 A from the mean structure for the backbone atoms N, Calpha, and C' of residues 2-11. The structure of ImI is characterized by two compact loops, defined by two disulfide bridges, which form distinct subdomains separated by a deep cleft. Two short 310-helical regions in the first loop are followed by a C-terminal beta-turn in the second. The two disulfide bridges and Ala 9 form a rigid hydrophobic core, orienting the other amino acid side chains toward the surface. Comparison of the three-dimensional structure of ImI to those of the larger, 16 amino acid alpha-conotoxins PnIA, PnIB, MII, and EpI-also specific for neuronal nAChRs-reveals remarkable similarity in local backbone conformations and relative solvent-accessible surface areas. The core scaffold is conserved in all five conotoxins, whereas the residues in solvent-exposed positions are highly variable. The second helical region, and the specific amino acids that the helix exposes to solvent, may be particularly important for binding and selectivity. This comparative analysis provides a three-dimensional structural basis for interpretation of mutagenesis data and structure-activity relationships for ImI as well other neuronal alpha-conotoxins.     

Expression of a synthetic gene encoding potato carboxypeptidase inhibitor using a bacterial secretion vector.

A synthetic gene encoding the 39-amino-acid (aa) potato carboxypeptidase inhibitor IIa (PCI-IIa) has been constructed and expressed using the secretion vector, pIN-III-ompA-3, fused in frame to the OmpA signal peptide-encoding sequence. Recombinant Escherichia coli secreted a PCI with 10 additional aa at the N terminus (rePCI + 10). These extra aa were removed by site-directed mutagenesis giving a PCI with no additional aa (rePCI), as shown by fast atom bombardment mass spectrometry (M(r) 4295). The two forms of rePCI were found almost exclusively in the culture medium, not in the periplasmic space, as would be expected from OmpA signal peptide fusions. Both rePCI + 10 and rePCI are biologically active and react strongly with serum raised against PCI from potato. A method for the purification of rePCI to homogeneity has been developed. The purified rePCI shows a Ki for carboxypeptidase A within the range of the natural PCI-IIa (1.5-2.7 nM). These results indicate that both rePCI + 10 and rePCI are properly folded and that their three disulfide bridges are correctly formed. Together with previous reports, our results show that fusion to a secretion signal peptide is an effective way of producing small proteins containing disulfide bridges in a biologically active form.     

Primary structure of Escherichia coli ribosomal protein L31.

Protein L31 from the 50S ribosomal subunit of Escherichia coli was manually sequenced by the dansyl-Edman method. Owing to the availability of only small quantities of purified L31, sequencing methods were scaled down such that the entire primary structure could be determined with 700 microgram of protein. The techniques employed are described in detail. The protein consists of a single chain of 62 amino acids, with a calculated molecular weight of 6967. Four half-cystine residues were identified at positions 16, 18, 37, and 40. Evidence is presented that suggests that these residues form two disulfide bridges in the protein, as isolated.     

The NMR structure of Escherichia coli ribosomal protein L25 shows homology to general stress proteins and glutaminyl-tRNA synthetases.

The structure of the Escherichia coli ribosomal protein L25 has been determined to an r.m.s. displacement of backbone heavy atoms of 0.62 +/- 0.14 A by multi-dimensional heteronuclear NMR spectroscopy on protein samples uniformly labeled with 15N or 15N/13C. L25 shows a new topology for RNA-binding proteins consisting of a six-stranded beta-barrel and two alpha-helices. A putative RNA-binding surface for L25 has been obtained by comparison of backbone 15N chemical shifts for L25 with and without a bound cognate RNA containing the eubacterial E-loop that is the site for binding of L25 to 5S ribosomal RNA. Sequence comparisons with related proteins, including the general stress protein, CTC, show that the residues involved in RNA binding are highly conserved, thereby providing further confirmation of the binding surface. Tertiary structure comparisons indicate that the six-stranded beta-barrels of L25 and of the tRNA anticodon-binding domain of glutaminyl-tRNA synthetase are similar.     

The acidic ribosomal phosphoprotein of eukaryotes and its relationship to ribosomal proteins L7 and L12 of Escherichia coli.

The acidic ribosomal phosphoprotein, Lgamma, of Krebs II ascites cells was further characterized and compared with proteins L7 and L12 of Escherichia coli. Ribosomal protein Lgamma was selectively removed from 60S sibosomal subunits by 50% ethanol and 1M-NH4Cl, and antibodies raised against protein Lgamma cross-reacted with E. coli protein L12 in immunodiffusion experiments. These and other, previously reported, data support the proposal that the uekaryotic counterpart of E. coli proteins L7 and L12 is phosphorylated.     

Disulfide connectivity of human immunoglobulin G2 structural isoforms

In this communication we present the detailed disulfide structure of IgG2 molecules. The consensus structural model of human IgGs represents the hinge region positioned as a flexible linker connecting structurally isolated Fc and Fab domains. IgG2 molecules are organized differently from that model and exhibit multiple structural isoforms composed of (heavy chain-light chain-hinge) covalent complexes. We describe the precise connection of all the disulfide bridges and show that the IgG2 C H1 and C-terminal C L cysteine residues are either linked to each other or to the two upper hinge cysteine residues specific to the IgG2 subclass. A defined arrangement of these disulfide bridges is unique to each isoform. Mutation of a single cysteine residue in the hinge region eliminates these natural complexes. These results show that IgG2 structure is significantly different from the conventionally accepted immunoglobulin structural model and may help to explain some of the unique biological activity attributed only to this subclass.     

Structure of human lysosomal membrane glycoprotein 1. Assignment of disulfide bonds and visualization of its domain arrangement

The amino acid sequence of one of the major lysosomal membrane glycoproteins, lysosome-associated membrane protein 1 (lamp-1), was deduced from its cDNA sequence (Fukuda, M., Viitala, J., Matteson, J., and Carlsson, S. R. (1988) J. Biol. Chem. 263, 18920-18928). This amino acid sequence suggests that lamp-1 contains a hinge-like structure and could form disulfide bridges that are observed in the immunoglobulin superfamily. To test this possibility, we have determined the positions of the disulfide bridges by isolating and sequencing cystine-containing peptides which contain disulfide bridges. The results indicate that disulfide arrangement of lamp-1 is different from that of immunoglobulins. Each molecule contains, in total, four loops formed by disulfide bonds, and each loop contains 36-39 amino acid residues. However, none of the disulfide bonds connects two domains that are separated by a hinge-like structure. The results indicate that the hinge region has no ordered structure, and the relative positions of the two domains can be altered in space. Examination of the ultrastructure of lamp-1 by electron microscopy showed that the hinge-like structure actually functions as a hinge. These results indicate that the lamp-1 molecule represents a novel family of glycoproteins with unique structural properties.     

Biochemical and structural characterization of TLXI, the Triticum aestivum L. thaumatin-like xylanase inhibitor

Thaumatin-like xylanase inhibitors (TLXI) are recently discovered wheat proteins. They belong to the family of the thaumatin-like proteins and inhibit glycoside hydrolase family 11 endoxylanases commonly used in different cereal based (bio)technological processes. We here report on the biochemical characterisation of TLXI. Its inhibition activity is temperature- and pH-dependent and shows a maximum at approximately 40 degrees C and pH 5.0. The TLXI structure model, generated with the crystal structure of thaumatin as template, shows the occurrence of five disulfide bridges and three beta-sheets. Much as in the structures of other short-chain thaumatin-like proteins, no alpha-helix is present. The circular dichroism spectrum of TLXI confirms the absence of alpha-helices and the presence of antiparallel beta-sheets. All ten cysteine residues in TLXI are involved in disulfide bridges. TLXI is stable for at least 120 min between pH 1-12 and for at least 2 hours at 100 degrees C, making it much more stable than the other two xylanase inhibitors from wheat, i.e. Triticum aestivum xylanase inhibitor (TAXI) and xylanase inhibitor protein (XIP). This high stability can probably be ascribed to the high number of disulfide bridges, much as seen for other thaumatin-like proteins.     

Structure of the Nidogen Binding LE Module of the Laminin γ1 Chain in Solution

The structure of the single LE module between residues 791 and 848 of the laminin γ1 chain, which contains the high affinity binding site for nidogen, has been probed using NMR methods. The module folds into an autonomous domain which has a stable and unique three-dimensional (3D) structure in solution. The 3D structure was determined on the basis of 362 interproton distance constraints derived from nuclear Overhauser enhancement measurements and 39 π angles, supplemented by 5 ψ and 22 χ₁angles. The main features of the NMR structures are two-stranded antiparallel β-sheets which are separated by loops and cross-connected by four disulfide bridges. The N-terminal segment which contains the first three disulfide bridges is similar to epidermal growth factor. The C-terminal segment has an S-like backbone profile with a crossover at the last disulfide bridge and comprises two three-residue long β-strands that form an antiparallel β-sheet. The LE module possesses an exposed nidogen binding loop that projects away from the main body of the protein. The side-chains of three amino acids which are crucial for binding (Asp, Asn, Val) are all exposed at the domain surface. An inactivating Asn-Ser mutation in this region showed the same 3D structure indicating that these three residues, and possibly an additional Tyr in an adjacent loop, provide direct contacts in the interaction with nidogen.