On the correlation of S-S and C-S raman bands frequencies and intensities with conformations of disulfide bridges in proteins

We have studied dependence of the S-S and C-S raman bands frequencies and intensities on the conformation of disulfide bridges. These correlations have been analyzed on the basis of the normal coordinates treatment in the valence force field and from comparison of the raman and X-ray data for the model compounds and proteins. The found regularity of the intensity variations with conformation make it possible to distinguish these spectral perturbations from those associated with the change in the number of disulfide bridges in proteins. Particularly, this was demonstrated on the gamma-irradiated insulin solutions spectra. We also present some data appropriate for estimation of the content of disulfide bridges in proteins using internal intensity standards.     

Disulfide connectivity prediction using recursive neural networks and evolutionary information

MOTIVATION: We focus on the prediction of disulfide bridges in proteins starting from their amino acid sequence and from the knowledge of the disulfide bonding state of each cysteine. The location of disulfide bridges is a structural feature that conveys important information about the protein main chain conformation and can therefore help towards the solution of the folding problem. Existing approaches based on weighted graph matching algorithms do not take advantage of evolutionary information. Recursive neural networks (RNN), on the other hand, can handle in a natural way complex data structures such as graphs whose vertices are labeled by real vectors, allowing us to incorporate multiple alignment profiles in the graphical representation of disulfide connectivity patterns. RESULTS: The core of the method is the use of machine learning tools to rank alternative disulfide connectivity patterns. We develop an ad-hoc RNN architecture for scoring labeled undirected graphs that represent connectivity patterns. In order to compare our algorithm with previous methods, we report experimental results on the SWISS-PROT 39 dataset. We find that using multiple alignment profiles allows us to obtain significant prediction accuracy improvements, clearly demonstrating the important role played by evolutionary information. AVAILABILITY: The Web interface of the predictor is available at http://neural.dsi.unifi.it/cysteines     

Learning about protein folding via potential functions

Over the last few years we have developed an empirical potential function that solves the protein structure recognition problem: given the sequence for an n-residue globular protein and a collection of plausible protein conformations, including the native conformation for that sequence, identify the correct, native conformation. Having determined this potential on the basis of only some 6500 native/nonnative pairs of structures for 58 proteins, we find it recognizes the native conformation for essentially all compact, soluble, globular proteins having known native conformations in comparisons with 10⁴ to 10⁶ reasonable alternative conformations apiece. In this sense, the potential encodes nearly all the essential features of globular protein conformational preference. In addition it “knows” about many additional factors in protein folding, such as the stabilization of multimeric proteins, quaternary structure, the role of disulfide bridges and ligands, proproteins vs. processed proteins, and minimal strand lengths in globular proteins. Comparisons are made with other sorts of protein folding problems, and applications in protein conformational determination and prediction are discussed. © 1994 Wiley-Liss, Inc.     

Secretion of Human Serum Albumin by Kluyveromyces lactis Overexpressing KlPDI1 and KlERO1

The control of protein conformation during translocation through the endoplasmic reticulum is often a bottleneck for heterologous protein production. The core pathway of the oxidative folding machinery includes two conserved proteins: Pdi1p and Ero1p. We increased the dosage of the genes encoding these proteins in the yeast Kluyveromyces lactis and evaluated the secretion of heterologous proteins. KlERO1, an orthologue of Saccharomyces cerevisiae ERO1, was cloned by functional complementation of the ts phenotype of an Scero1 mutant. The expression of KlERO1 was induced by treatment of the cells with dithiothreitol and by overexpression of human serum albumin (HSA), a disulfide bond-rich protein. Duplication of either PDI1 or ERO1 led to a similar increase in HSA yield. Duplication of both genes accelerated the secretion of HSA and improved cell growth rate and yield. Increasing the dosage of KlERO1 did not affect the production of human interleukin 1β, a protein that has no disulfide bridges. The results confirm that the ERO1 genes of S. cerevisiae and K. lactis are functionally similar even though portions of their coding sequence are quite different and the phenotypes of mutants overexpressing the genes differ. The marked effects of KlERO1 copy number on the expression of heterologous proteins with a high number of disulfide bridges suggests that control of KlERO1 and KlPDI1 is important for the production of high levels of heterologous proteins of this type.     

A Molecular Dynamics Study of Human Defensins HBD-1 and HNP-3 in Water

Abstract

Mammalian defensins are crucial components of the innate immune system. They are characterized by three disulfide bridges and exhibit broad spectrum antibacterial activity. The spacing between the cysteines and disulfide connectivities in the two classes of defensins, the α- and β-forms, are different. The structural motif of 3 β-strands appears to be conserved in α and β-defensins despite differences in disulfide connectivities and spacing between cysteines. In this study, Molecular Dynamics Simulations (MDS) have been carried out to study the conformational behavior of α- and β-defensins with and without disulfide bridges. Our results indicate that β-strands in the C-terminal region of HBD-1 and HNP-3 do not unfold during the course of MDS. The segment adopting α-helix in HBD-1 unfolds early during the simulations. The backbone hydrogen bonds in HBD-1 and HNP-3 are broken during MDS. When the disulfide bonds are absent, the N-terminal β-strand unfolds by 20 ns but β-strands are observed in the C-terminal region of HNP-3. HBD-1, without disulfide bridges, unfolds to a greater extent during the course of the MDS. Examination of distances between sulfur atoms of cysteines without disulfide bridges during the simulations indicate that there is no specific preference for native disulfide bridges, which could be the reason for the experimental observation of non-native disulfide bridge formation during chemical synthesis of human α- and β-defensins. Since defensins with non-native disulfide bridges are biologically active, the exact three dimensional structures observed for native HBD-1 and HNP-3 does not appear to be essential for exhibiting antibacterial activity.     

Contribution of allosteric disulfide in the structural regulation of membrane-bound tissue factor–factor VIIa binary complex

Abstract

Two distinct populations, active and cryptic forms of tissue factor (TF), reside on the cell surface. Apart from phospholipid contribution, various models have been introduced to explain decryption/encryption of TF. The proposed model, the switching of Cys186–Cys209 bond of TF, has become the matter of controversy. However, it is well accepted that this disulfide has an immense influence upon ligand factor VIIa (FVIIa) for its binding. However, molecular level understanding for this remains unveiled due to lack of detailed structural information. In this regard, we have performed the molecular dynamic study of membrane-bound TF/TF–FVIIa in both the forms (±Cys186–Cys209 allosteric disulfide bond), individually. Dynamic study depicts that disulfide bond provides structural rigidity of TF in both free and ligand-bound forms. This disulfide bond also governs the conformation of FVIIa structure as well as the binding affinity of FVIIa toward TF. Significant differences in lipid–protein interaction profiles of both the forms of TF in the complex were observed. Two forms of TF, oxidized and reduced, have different structural conformation and behave differentially toward its ligand FVIIa. This disulfide bond not only alters the conformation of GLA domain of FVIIa in the vicinity but allosterically regulates the conformation of the distantly located FVIIa protease domain. We suggest that the redox status of the disulfide bond also governs the lipid-mediated interactions with both TF and FVIIa.

Communicated by Ramaswamy H. Sarma     

The crystallographically determined structures of atypical strained disulfides engineered into subtilisin

The geometries of two disulfide bridges genetically engineered into subtilisin have been characterized by x-ray crystallography to determine the structural and energetic constraints involved in introducing disulfide bonds into proteins. Both disulfide bridges (Cys-24-Cys-87 and Cys-22-Cys-87) exhibit atypical sets of dihedral angles compared to those for other reported disulfide structures in proteins. The geometric trends for naturally occurring disulfides in protein crystal structures are examined. Comparison of the disulfide-containing mutant protein structures with the wild-type structure shows that, in both cases, disulfide incorporation is accommodated by relatively minor changes in local main-chain conformation. The Cys-22-Cys-87 disulfide has two high energy dihedral angles (X2 = 121 degrees, X2' = 143 degrees). Both disulfides produce short non-bonded contacts with the main-chain.     

Derivation of rules for comparative protein modeling from a database of protein structure alignments

We describe a database of protein structure alignments as well as methods and tools that use this database to improve comparative protein modeling. The current version of the database contains 105 alignments of similar proteins or protein segments. The database comprises 416 entries, 78,495 residues, 1,233 equivalent entry pairs, and 230,396 pairs of equivalent alignment positions. At present, the main application of the database is to improve comparative modeling by satisfaction of spatial restraints implemented in the program MODELLER (?ali A, Blundell TL, 1993, J Mol Biol 234:779–815). To illustrate the usefulness of the database, the restraints on the conformation of a disulfide bridge provided by an equivalent disulfide bridge in a related structure are derived from the alignments; the prediction success of the disulfide dihedral angle classes is increased to approximately 80%, compared to approximately 55% for modeling that relies on the stereochemistry of disulfide bridges alone. The second example of the use of the database is the derivation of the probability density function for comparative modeling of the cis/trans isomerism of the proline residues; the prediction success is increased from 0% to 82.9% for cis-proline and from 93.3% to 96.2% for trans-proline. The database is available via electronic mail.     

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°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 β-sheets. Much as in the structures of other short-chain thaumatin-like proteins, no α-helix is present. The circular dichroism spectrum of TLXI confirms the absence of α-helices and the presence of antiparallel β-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°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.     

X-ray diffraction studies on muscle regulation

Using the intensity of the outer part of the second actin layer line as an indicator of thin filament conformation in vertebrate muscle we were able to identify the four different states of rest, and the three states induced by the presence of Ca²⁺ ions, rigor bridge attachment and actively cycling bridges, respectively. These findings are in qualitative agreement with a number of biochemical studies by Eisenberg and Greene and others, indicating that activation of the thin filament depends both on Ca²⁺ ions and crossbridge binding. Yet quantitatively, the biochemical data and our structural data are contradictory. Whereas the biochemical studies suggest a strong coupling between structural changes of the thin filament and the ATPase activity, the structural studies indicate that this is not necessarily the case.Troponin molecules also change their conformation upon activation depending on both Ca²⁺ ions and crossbridge binding as demonstrated by the early part of the time course of the thin filament meridional reflections in contracting frog muscle.Low ionic strength which has been shown by Brenner and collaborators to increase weakly binding crossbridges in relaxed rabbit psoas muscle does not influence the intensity of the second actin layer line in this muscle. Yet in contracting frog muscle the increase of the second actin layer line increases very rapidly in one step, suggesting that weak binding bridges which are attached to actin prior to force production may indeed influence the thin filament conformation. It therefore appears that weakly bound bridges in the low ionic strength state do not have the same effect on the thin filament conformation as weakly bound bridges in an actively contracting muscle.Arthropod muscles like the thin filament regulated lobster muscles differ from vertebrate muscle in not showing an increase of the second layer line during contraction, which may have to do with differences in crossbridge attachment. The myosin-regulated molluscan muscle ABRM shows a large increase on the second actin layer line upon phasic contraction and a much smaller increase in catch or rigor, indicating that actively cycling bridges influence the thin filament conformation differently than catch or rigor bridges.Several pieces of evidence which we have briefly outlined in this paper suggest that the thin filament conformational changes we have observed do not arise solely from tropomyosin movements and that conformational changes of actin domains should be considered.     

Modification of protein stability by introduction of disulfide bridges and prolines: geometric criteria for mutation sites

We define geometrical parameters to characterize disulfide bridges using x-ray crystal structure data on small molecules and use them to suggest replacements of amino acids by cysteines in order to introduce disulfide bridges to increase thermal stability in proteins. We also define geometric parameters to identify target amino acids for replacements by prolines in order to conserve desired structural attributes in the vicinity of disulfide mutations leading to further structural and thermal stability of proteins. The geometric criteria are applied to the serine protease, subtilisin, to model stereochemically favorable disulfide mutants without altering the active site geometry, implying conservation of native biological activity.     

Synthesis, solution structure, binding activity, and cGMP activation of human guanylin and its disulfide isomer

Guanylin is a recently isolated peptide consisting of 15 amino acid residues with four cysteines, which may form two intramolecular disulfide bridges, and stimulates intestinal membrane guanylate cyclase. The position of the disulfide linkages of guanylin was predicted from its structural similarity to a heat stable enterotoxin which is thought to be responsible for secretory diarrhoea. Both guanylin, with disulfide positions 4–12 and 7–15, and its disulfide isomer, with disulfides positions 4–15 and 7–12, were chemically synthesized by the solid-phase method and purified. Two specific disulfides were selectively formed and confirmed by sequencing, mass spectrometry and high-performance liquid chromatography in combination with enzymatic cleavage. The structure of both isomers has been investigated in solution by ¹H nuclear magnetic resonance spectroscopy. Guanylin exists as a mixture of two stable conformations which have compact spiral structures, from comparison with literature data. In contrast, the disulfide isomer of guanylin shows only a single conformation with an elongated curved plate-like structure. Binding assays were performed using labelled guanylin with membranes obtained from rat jejunum. Both disulfide isomers were investigated by the cGMP assay. Both binding and cGMP assays indicated that the relevant form of disulfide bridges in the intact guanylin was as predicted.     

Kiwi protein inhibitor of pectin methylesterase amino-acid sequence and structural importance of two disulfide bridges.

A protein acting as a powerful inhibitor of plant pectin methylesterase was isolated from kiwi (Actinidia chinensis) fruit. The complete amino-acid sequence of the pectin methylesterase inhibitor (PMEI) was determined by direct protein analysis. The sequence comprises 152 amino-acid residues, accounting for a molecular mass of 16 277 Da. The far-UV CD spectrum indicated a predominant alpha-helix conformation in the secondary structure. The protein has five cysteine residues but neither tryptophan nor methionine. Analysis of fragments obtained after digestion of the protein alkylated without previous reduction identified two disulfide bridges connecting Cys9 with Cys18, and Cys74 with Cys114; Cys140 bears a free thiol group. A database search pointed out a similarity between PMEI and plant invertase inhibitors. In particular, the four Cys residues, which in PMEI are involved in the disulfide bridges, are conserved. This allows us to infer that also in the homologous proteins, whose primary structure was deduced only by cDNA sequencing, those cysteine residues are engaged in two disulfide bridges, and constitute a common structural motif. The comparison of the sequence of these inhibitors confirms the existence of a novel class of proteins with moderate but significant sequence conservation, comprising plant proteins acting as inhibitors of sugar metabolism enzymes, and probably involved in various steps of plant development.     

An empirical approach to protein conformation stability and flexibility

Experimental measurements of disulfide bond stability at various stages of protein folding are considered in terms of the effective concentrations of the thiol groups relative to each other; values of up to 10⁷M are observed, so that intramolecular interactions within the interior of a protein are much more stable, and provide greater stability to the folded conformation, than those on the surface or in a flexible segment. Intramolecular interactions can have substantially lower free energies than intermolecular, for solely entropic reasons; this implies that polar interactions, such as hydrogen bonds and salt bridges, can provide net stabilization to a folded conformation, in spite of the unfolded protein having intermolecular interactions with the solvent. These considerations can account for the lower free energy and enthalpy of the folded state and are useful for considering protein flexibility.     

Structure-activity relationship studies of gomesin: importance of the disulfide bridges for conformation, bioactivities, and serum stability

Gomesin is an antimicrobial peptide isolated from hemocytes of the Brazilian spider Acanthoscurria gomesiana that contains two disulfide bridges Cys(2-15)/Cys(6-11) and presents a beta-hairpin structure. To investigate the role of the disulfide bridges on gomesin conformation, bioactivities, and serum stability, structure-activity relationship (SAR) studies were conducted. Initially, gomesin and variants lacking one or both disulfide bridges were synthesized. CD studies showed that the gomesin structure is very rigid independently of the solvent environment. On the other hand, the linearized analogues adopted secondary structures according to the environment, while the monocyclic disulfide-bridged peptides had a tendency to adopt a turn structure. The absence of one or both bridges resulted in a decrease in the antimicrobial and hemolytic activities. In addition, serum stability studies revealed that, contrasting to gomesin that was stable even after 48 h of incubation, the linearized analogues were rapidly degraded. The replacement of the disulfide bounds by lactam bridges led to monocyclic and bicyclic compounds. SAR studies indicated that the monocyclic lactam-bridged analogues tend to assume a alpha-helical structure being less potent, hemolytic, and serum stable than the wild-type gomesin. On the other hand, the bicyclic lactam/disulfide-bridged analogues displayed a similar conformation and degradation kinetics identical to gomesin. However, the antimicrobial activity appeared to be dependent on the lactam bridge position and size. These findings indicated that (i) the secondary structure plays a pivotal role for the full activity of gomesin; (ii) the antimicrobial and hemolytic activities of gomesin are correlated events; (iii) while at least one of the disulfide bridges is needed for the maintenance of a significant antimicrobial activity of gomesin, both bridges are required for high serum stability and optimal conformation; and finally (iv) the best analogue obtained was the bicyclo (2-15,6-11)[Glu2, Cys(6,11), Lys15]-Gm since it is as stable and potent as gomesin.     

S-glutathionylation of human glyceraldehyde-3-phosphate dehydrogenase and possible role of Cys152-Cys156 disulfide bridge in the active site of the protein

BackgroundWe previously showed that glyceraldehyde-3-phosphate dehydrogenase (GAPDH) is S-glutathionylated in the presence of H₂O₂ and GSH. S-glutathionylation was shown to result in the formation of a disulfide bridge in the active site of the protein. In the present work, the possible biological significance of the disulfide bridge was investigated.MethodsHuman recombinant GAPDH with the mutation C156S (hGAPDH_C156S) was obtained to prevent the formation of the disulfide bridge. Properties of S-glutathionylated hGAPDH_C156S were studied in comparison with those of the wild-type protein hGAPDH.ResultsS-glutathionylation of hGAPDH and hGAPDH_C156S results in the reversible inactivation of the proteins. In both cases, the modification results in corresponding mixed disulfides between the catalytic Cys152 and GSH. In the case of hGAPDH, the mixed disulfide breaks down yielding Cys152-Cys156 disulfide bridge in the active site. In hGAPDH_C156S, the mixed disulfide is stable. Differential scanning calorimetry method showed that S-glutathionylation leads to destabilization of hGAPDH molecule, but does not affect significantly hGAPDH_C156S. Reactivation of S-glutathionylated hGAPDH in the presence of GSH and glutaredoxin 1 is approximately two-fold more efficient compared to that of hGAPDH_C156S.ConclusionsS-glutathionylation induces the formation of Cys152-Cys156 disulfide bond in the active site of hGAPDH, which results in structural changes of the protein molecule. Cys156 is important for reactivation of S-glutathionylated GAPDH by glutaredoxin 1.General significanceThe described mechanism may be important for interaction between GAPDH and other proteins and ligands, involved in cell signaling.     

Low thermodynamic but high kinetic stability of an antifreeze protein from Rhagium mordax

The equilibrium heat stability and the kinetic heat tolerance of a recombinant antifreeze protein (AFP) from the beetle Rhagium mordax (RmAFP1) are studied through differential scanning calorimetry and circular dichroism spectroscopy. In contrast to other insect AFPs studied with this respect, the RmAFP1 has only one disulfide bridge. The melting temperature, T_m, of the protein is determined to be 28.5°C (pH 7.4), which is much lower than most of those reported for AFPs or globular proteins in general. Despite its low melting temperature, both biophysical and activity measurements show that the protein almost completely refolds into the native state after repeated exposure of 70°C. RmAFP1 thus appears to be kinetically stable even far above its melting temperature. Thermodynamically, the insect AFPs seem to be dividable in three groups, relating to their content of disulfide bridges and widths of the ice binding motifs; high melting temperature AFPs (high disulfide content, TxT motifs), low melting temperature but high refolding capability AFPs (one disulfide bridge, TxTxTxT motifs) and irreversibly unfolded AFPs at low temperatures (no disulfide bridges, TxTxTxTxT motifs). The property of being able to cope with high temperature exposures may appear peculiar for proteins which strictly have their effect at subzero temperatures. Different aspects of this are discussed.     

Cysteine Residues and the Structure of the Rat Renal Proximal Tubular Type II Sodium Phosphate Cotransporter (Rat NaPi IIa)

The rat renal Na/P _i cotransporter type IIa (rat NaP_i IIa) is a 637 amino acid protein containing 12 cysteine residues. We examined the effect of different cysteine modifying methanethiosulfonate (MTS)-reagents and the disulfide bond reducing agent tris(2-carboxyethyl)phosphine (TCEP) on the transport activity of wild-type and 12 single cysteine substitution mutants of rat NaPi IIa expressed in Xenopus laevis oocytes. The transport activity of the wild-type protein was resistant to three membrane impermeant MTS-reagents (MTSEA, MTSET and MTSES). In contrast, membrane permeant methyl methanethiosulfonate (MMTS) and TCEP inhibited the transport activity of both the wild-type, as well as all the single mutant proteins. This indicated the existence of more than one functionally important cysteine residue, not accessible extracellularly, and at least 2 disulfide bridges. To identify the disulfide bridges, three double mutants lacking 2 of the 3 cysteine residues predicted to be extracellular in different combinations were examined. This led to the identification of one disulfide bridge between C306 and C334; reconsideration of the topological model predictions suggested a second disulfide bridge between C225 and C520. Evaluation of a fourth double mutant indicated that at least one of two disulfide bridges (C306 and C334; C225 and C520) has to be formed to allow the surface expression of a functional cotransporter. A revised secondary structure is proposed which includes two partially repeated motifs that are connected by disulfide bridges formed between cysteine pairs C306-C334 and C225-C520. Received: 13 December 1999/Revised: 31 March 2000     

Conformational changes of α-lactalbumin induced by the stepwise reduction of its disulfide bridges: The effect of the disulfide bridges on the structural stability of the protein in sodium dodecyl sulfate solution

Four disulfide bridges of bovineα-lactalbumin (α-lact) were selectively reduced to obtain its derivatives with three, two, and zero disulfide bridges (designated as 3SS, 2SS, and OSSα-lact, respectively). The original helicity was almost maintained in 3SSα-lact missing only the Cys6-Cysl20 bridge. Upon the reduction of both Cys28-Cys111 and Cys6-Cys120 bridges, various changes occurred in the protein. In particular, the maximum fluorescence of 1-anilinonaphthalene-8-sulfonic acid was observed in this stage. Upon the reduction of all disulfide bridges, the hydrophobic box of the protein, formed by Trp60, Ile95, Tyr103, and Trp104, was disrupted and an internal helical structure was destroyed. The conformation of each derivative was examined mainly in a solution of sodium dodecyl sulfate. In the surfactant solution, the helicity increased from 33% to 37% in 3SSα-lact, from 26% to 31% in 2SSα-lact, and from 18% to 37% in OSSα-lact, as against from 34% to 44% in intactα-lact. On the other hand, the tryptophan fluorescence of each derivative was affected in very low surfactant concentrations, suggesting that the tertiary structure considerably changed prior to the secondary structural change in the surfactant solution.     

Conformational study of a native monodisulfide bridge analogue of EETI II

Trypsin inhibitor EETI II, possessing six cysteinesengaged in three disulfide bridges, shares a commonstructural motif with other proteins of differentorigins and functions. To understand the principlesthat govern folding of this largely distributed basicscaffold, mainly composed of a small triple-stranded-sheet, we have studied different stages in thefolding of EETI II. The conformational properties ofa synthetic analogue of EETI II possessing only onenative (15-27) disulfide bridge were investigated withthe combined use of ¹H NMR and molecularmodelling. Although two native-like reverse turns wereobserved, formation of -sheet could not beevidenced in the one disulfide analogue, while themotif has been shown to be present in a foldingintermediate with two native disulfide bridges (9-21and 15-27). These results suggest that the structuralmotif requires stabilisation by two disulfide bridges