The effects of disulfide bonds on the denatured state of barnase

Conformational dynamics in peptide toxins: Implications for receptor interactions and molecular design

Peptide toxins are potent and often exquisitely selective probes of the structure and function of ion channels and receptors, and are therefore of significant interest to the pharmaceutical and biotech industries as both pharmacological tools and therapeutic leads. The three-dimensional structures of peptide toxins are essential as a basis for understanding their structure-activity relationships and their binding to target receptors, as well as in guiding the design of analogues with modified potency and/or selectivity for key targets. NMR spectroscopy has played a key role in elucidating the structures of peptide toxins and probing their structure-function relationships. In this article, we highlight the additional important contribution of NMR to characterising the dynamics of peptide toxins. We also compare the information available from NMR measurements with that afforded by molecular dynamics simulations. We describe several examples of the importance of dynamics measurements over a range of timescales for understanding the structure-function relationships of peptide toxins and their receptor engagement. Peptide toxins that inhibit the voltage-gated potassium channel K_V1.3 with pM affinities display different degrees of conformational flexibility, even though they contain multiple disulfide bonds, and this flexibility can affect the relative orientation of residues that have been shown to be critical for channel binding. Information on the dynamic properties of peptide toxins is important in the design of analogues or mimetics where receptor-bound structures are not available.

Modulating charge patterning and ionic strength as a strategy to induce conformational changes in intrinsically disordered proteins

We present an analytical theory to describe conformational changes as a function of salt for polymers with a given sequence of charges. We apply this model to describe Intrinsically Disordered Proteins (IDPs) by explicitly accounting for charged residues and their exact placement in the primary sequence while approximating the effect of non-electrostatic interactions at a mean-field level by effective short-range (two body and three-body) interaction parameters. The effect of ions is introduced by treating electrostatic interactions within Debye-Huckle approximation. Using typical values of the short-range mean-field parameters derived from all-atom Monte Carlo simulations (at zero salt), we predict the conformational changes as a function of salt concentration. We notice that conformational transitions in response to changes in ionic strength strongly depend on sequence specific charge patterning. For example, globule to coil transition can be observed upon increasing salt concentration, in stark contrast to uniformly charged polyelectrolyte theories based on net charge only. In addition, it is possible to observe non-monotonic behavior with salt as well. Drastic differences in salt-induced conformational transitions is also evident between two doubly phosphorylated sequences-derived from the same wild type sequence-that only differ in the site of phosphorylation. Similar effects are also predicted between two sequences derived from the same parent sequence differing by a single site mutation where a negative charge is replaced by a positive charge. These effects are purely a result of charge decoration and can only be understood in terms of metrics based on specific placement of charges, and cannot be explained by models based on charge composition alone. Identifying sequences and hot spots within sequences-for post translational modification or charge mutation-using our high-throughput theory will yield fundamental insights into design and biological regulation mediated by phosphorylation and/or local changes in salt concentration.

Effect of an Imposed Contact on Secondary Structure in the Denatured State of Yeast Iso-1-cytochrome c

There is considerable evidence that long-range interactions stabilize residual protein structure under denaturing conditions. However, evaluation of the effect of a specific contact on structure in the denatured state has been difficult. Iso-1-cytochrome c variants with a Lys54 → His mutation form a particularly stable His-heme loop in the denatured state, suggestive of loop-induced residual structure. We have used multidimensional nuclear magnetic resonance methods to assign ¹H and ¹⁵N backbone amide and ¹³C backbone and side chain chemical shifts in the denatured state of iso-1-cytochrome c carrying the Lys54 → His mutation in 3 and 6 M guanidine hydrochloride and at both pH 6.4, where the His54-heme loop is formed, and pH 3.6, where the His54-heme loop is broken. Using the secondary structure propensity score, with the 6 M guanidine hydrochloride chemical shift data as a random coil reference state for data collected in 3 M guanidine hydrochloride, we found residual helical structure in the denatured state for the 60s helix and the C-terminal helix, but not in the N-terminal helix in the presence or absence of the His54-heme loop. Non-native helical structure is observed in two regions that form Ω-loops in the native state. There is more residual helical structure in the C-terminal helix at pH 6.4 when the loop is formed. Loop formation also appears to stabilize helical structure near His54, consistent with induction of helical structure observed when His-heme bonds form in heme-peptide model systems. The results are discussed in the context of the folding mechanism of cytochrome c.

Different dynamics and pathway of disulfide bonds reduction of two human defensins, a molecular dynamics simulation study

Human defensins are a class of antimicrobial peptides that are crucial components of the innate immune system. Both human α defensin type 5 (HD5) and human β defensin type 3 (hBD-3) have 6 cysteine residues which form 3 pairs of disulfide bonds in oxidizing condition. Disulfide bond linking is important to the protein structure stabilization, and the disulfide bond linking and breaking order have been shown to influence protein function. In this project, microsecond long molecular dynamics simulations were performed to study the structure and dynamics of HD5 and hBD-3 wildtype and analogs which have all 3 disulfide bonds released in reducing condition. The structure of hBD-3 was found to be more dynamic and flexible than HD5, based on RMSD, RMSF, and radius of gyration calculations. The disulfide bridge breaking order of HD5 and hBD-3 in reducing condition was predicted by two kinds of methods, which gave consistent results. It was found that the disulfide bonds breaking pathways for HD5 and hBD-3 are very different. The breaking of disulfide bonds can influence the dimer interface by making the dimer structure less stable for both kinds of defensin. In order to understand the difference in dynamics and disulfide bond breaking pathway, hydrophilic and hydrophobic accessible surface areas (ASA), buried surface area between cysteine pairs, entropy of cysteine pairs, and internal energy were calculated. Comparing to the wildtype, hBD-3 analog is more hydrophobic, while HD5 is more hydrophilic. For hBD-3, the disulfide breaking is mainly entropy driven, while other factors such as the solvation effects may take the major role in controlling HD5 disulfide breaking pathway. Proteins 2017; 85:665-681. © 2016 Wiley Periodicals, Inc.

Insights from molecular dynamics simulations for computational protein design

A grand challenge in the field of structural biology is to design and engineer proteins that exhibit targeted functions. Although much success on this front has been achieved, design success rates remain low, an ever-present reminder of our limited understanding of the relationship between amino acid sequences and the structures they adopt. In addition to experimental techniques and rational design strategies, computational methods have been employed to aid in the design and engineering of proteins. Molecular dynamics (MD) is one such method that simulates the motions of proteins according to classical dynamics. Here, we review how insights into protein dynamics derived from MD simulations have influenced the design of proteins. One of the greatest strengths of MD is its capacity to reveal information beyond what is available in the static structures deposited in the Protein Data Bank. In this regard simulations can be used to directly guide protein design by providing atomistic details of the dynamic molecular interactions contributing to protein stability and function. MD simulations can also be used as a virtual screening tool to rank, select, identify, and assess potential designs. MD is uniquely poised to inform protein design efforts where the application requires realistic models of protein dynamics and atomic level descriptions of the relationship between dynamics and function. Here, we review cases where MD simulations was used to modulate protein stability and protein function by providing information regarding the conformation(s), conformational transitions, interactions, and dynamics that govern stability and function. In addition, we discuss cases where conformations from protein folding/unfolding simulations have been exploited for protein design, yielding novel outcomes that could not be obtained from static structures.

Role of Proteome Physical Chemistry in Cell Behavior

We review how major cell behaviors, such as bacterial growth laws, are derived from the physical chemistry of the cell’s proteins. On one hand, cell actions depend on the individual biological functionalities of their many genes and proteins. On the other hand, the common physics among proteins can be as important as the unique biology that distinguishes them. For example, bacterial growth rates depend strongly on temperature. This dependence can be explained by the folding stabilities across a cell’s proteome. Such modeling explains how thermophilic and mesophilic organisms differ, and how oxidative damage of highly charged proteins can lead to unfolding and aggregation in aging cells. Cells have characteristic time scales. For example, E. coli can duplicate as fast as 2–3 times per hour. These time scales can be explained by protein dynamics (the rates of synthesis and degradation, folding, and diffusional transport). It rationalizes how bacterial growth is slowed down by added salt. In the same way that the behaviors of inanimate materials can be expressed in terms of the statistical distributions of atoms and molecules, some cell behaviors can be expressed in terms of distributions of protein properties, giving insights into the microscopic basis of growth laws in simple cells.

Implication of disulfide bridge induced thermal reversibility, structural and functional stability for luciferase

Firefly luciferase is a relatively unstable protein and commonly loses its activity at room temperature because of structural changes. The structural and functional stability of this protein is critical for its enzymatic applications. Different approaches are applied to increase the stability of this enzyme such as designing of covalent cross-links (disulfide bonds). In this study, luciferase mutants containing one or two disulfide bonds were compared to the native protein for their for their structural, thermodynamic, and functional properties. Mutant forms of P. Pyralis luciferase A²⁹⁶C-A³²⁶C and A²⁹⁶C-A³²⁶C/P⁴⁵¹C-V⁴⁶⁹C were used. Thermodynamic and biophysical studies were carried out using UV-Vis, fluorescence, circular dichroism, luminescence spectroscopy and differential scanning calorimetry (DSC). We observed that both mutant forms of the protein were more stable than the wild-type enzyme. However, the single disulfide bond containing mutant was structurally and functionally more stable than the mutant protein containing two disulfide bonds. Furthermore, the enzymatic activity of the single disulfide bond containing mutant protein was 7-folds greater than the wild type and the double disulfide bond proteins. The A²⁹⁶C-A³²⁶C mutation also increased the reversibility and disaggregation of the protein. The enhanced activity of the single disulfide bond mutant protein was contributed to the expansion of its active site cleft, which was confirmed by bioinformatics tools.

The contribution of two disulfide bonds in the trypsin binding domain of horsegram (Dolichos biflorus) Bowman-Birk inhibitor to thermal stability and functionality

The major Bowman-Birk inhibitor (BBIs) of horsegram (Dolichos biflorus) HGI-III, contains seven interweaving disulfides and is extremely stable to high temperatures. The contributions of two disulfide bonds in the trypsin domain to thermal stability and functionality were evaluated using disulfide deletion variants of wild type protein. Thermal denaturation kinetics, differential scanning calorimetry and urea denaturation studies indicate that the absence of either of the two disulfides destabilizes the protein significantly. C20-C66 contributes substantially to both thermal stability and controls trypsin and chymotrypsin inhibitor activity. These two disulfides act in synergy as deletion of both disulfides leads to a complete loss of thermal stability. The data indicate that the two subdomains are not entirely independent of each other. Long range interactions, between the domains are facilitated by C20-C66. The deletion of the disulfide bonds also increased proteolytic susceptibility in a manner similar to the decreased thermal stability. From this study of rHGI a prototype of legume BBIs in can be concluded that among the array of seven evolutionarily conserved disulfide bonds, the disulfide C20-C66 that connects a residue in the trypsin domain with a residue at the border of the same domain plays a dominant role in maintaining functional and structural stability.

Disulfide bonding in protein biophysics

It has been known for many decades that cell surface, soluble-secreted, and extracellular matrix proteins are generally rich in disulfide bonds, but only more recently has the functional diversity of disulfide bonding in extracellular proteins been appreciated. In addition to the classic mechanisms by which disulfide bonds enhance protein thermodynamic stability, disulfides in certain configurations contribute particular mechanical properties to proteins that sense and respond to tensile forces. Disulfides may help warp protein folds for the evolution of new functions, or they may fasten aggregation-prone flaps of polypeptide to protein surfaces to prevent fibrilization or oligomerization. Disulfides can also be used to package and secure macromolecular cargo for intercellular transport. A series of case studies illustrating diverse biophysical roles of disulfide bonding are reviewed, with a focus on proteins functioning in the extracellular environment.

Multi-state proteins: approach allowing experimental determination of the formation order of structure elements in the green fluorescent protein

The most complex problem in studying multi-state protein folding is the determination of the sequence of formation of protein intermediate states. A far more complex issue is to determine at what stages of protein folding its various parts (secondary structure elements) develop. The structure and properties of different intermediate states depend in particular on these parts. An experimental approach, named μ-analysis, which allows understanding the order of formation of structural elements upon folding of a multi-state protein was used in this study. In this approach the same elements of the protein secondary structure are “tested” by substitutions of single hydrophobic amino acids and by incorporation of cysteine bridges. Single substitutions of hydrophobic amino acids contribute to yielding information on the late stages of protein folding while incorporation of ss-bridges allows obtaining data on the initial stages of folding. As a result of such an μ-analysis, we have determined the order of formation of beta-hairpins upon folding of the green fluorescent protein.

Dual beneficial effect of interloop disulfide bond for single domain antibody fragments

The antigen-binding fragment of functional heavy chain antibodies (HCAbs) in camelids comprises a single domain, named the variable domain of heavy chain of HCAbs (VHH). The VHH harbors remarkable amino acid substitutions in the framework region-2 to generate an antigen-binding domain that functions in the absence of a light chain partner. The substitutions provide a more hydrophilic, hence more soluble, character to the VHH but decrease the intrinsic stability of the domain. Here we investigate the functional role of an additional hallmark of dromedary VHHs, i.e. the extra disulfide bond between the first and third antigen-binding loops. After substituting the cysteines forming this interloop cystine by all 20 amino acids, we selected and characterized several VHHs that retain antigen binding capacity. Although VHH domains can function in the absence of an interloop disulfide bond, we demonstrate that its presence constitutes a net advantage. First, the disulfide bond stabilizes the domain and counteracts the destabilization by the framework region-2 hallmark amino acids. Second, the disulfide bond rigidifies the long third antigen-binding loop, leading to a stronger antigen interaction. This dual beneficial effect explains the in vivo antibody maturation process favoring VHH domains with an interloop disulfide bond.

Physical limits of cells and proteomes

What are the physical limits to cell behavior? Often, the physical limitations can be dominated by the proteome, the cell's complement of proteins. We combine known protein sizes, stabilities, and rates of folding and diffusion, with the known protein-length distributions P(N) of proteomes (Escherichia coli, yeast, and worm), to formulate distributions and scaling relationships in order to address questions of cell physics. Why do mesophilic cells die around 50 °C? How can the maximal growth-rate temperature (around 37 °C) occur so close to the cell-death temperature? The model shows that the cell's death temperature coincides with a denaturation catastrophe of its proteome. The reason cells can function so well just a few degrees below their death temperature is because proteome denaturation is so cooperative. Why are cells so dense-packed with protein molecules (about 20% by volume)? Cells are packed at a density that maximizes biochemical reaction rates. At lower densities, proteins collide too rarely. At higher densities, proteins diffuse too slowly through the crowded cell. What limits cell sizes and growth rates? Cell growth is limited by rates of protein synthesis, by the folding rates of its slowest proteins, and--for large cells--by the rates of its protein diffusion. Useful insights into cell physics may be obtainable from scaling laws that encapsulate information from protein knowledge bases.

How do thermophilic proteins and proteomes withstand high temperature?

We attempt to understand the origin of enhanced stability in thermophilic proteins by analyzing thermodynamic data for 116 proteins, the largest data set achieved to date. We compute changes in entropy and enthalpy at the convergence temperature where different driving forces are maximally decoupled, in contrast to the majority of previous studies that were performed at the melting temperature. We find, on average, that the gain in enthalpy upon folding is lower in thermophiles than in mesophiles, whereas the loss in entropy upon folding is higher in mesophiles than in thermophiles. This implies that entropic stabilization may be responsible for the high melting temperature, and hints at residual structure or compactness of the denatured state in thermophiles. We find a similar trend by analyzing a homologous set of proteins classified based only on the optimum growth temperature of the organisms from which they were extracted. We find that the folding free energy at the temperature of maximal stability is significantly more favorable in thermophiles than in mesophiles, whereas the maximal stability temperature itself is similar between these two classes. Furthermore, we extend the thermodynamic analysis to model the entire proteome. The results explain the high optimal growth temperature in thermophilic organisms and are in excellent quantitative agreement with full thermal growth rate data obtained in a dozen thermophilic and mesophilic organisms.

Analysis on conservation of disulphide bonds and their structural features in homologous protein domain families

Background

Disulphide bridges are well known to play key roles in stability, folding and functions of proteins. Introduction or deletion of disulphides by site-directed mutagenesis have produced varying effects on stability and folding depending upon the protein and location of disulphide in the 3-D structure. Given the lack of complete understanding it is worthwhile to learn from an analysis of extent of conservation of disulphides in homologous proteins. We have also addressed the question of what structural interactions replaces a disulphide in a homologue in another homologue.

Results

Using a dataset involving 34,752 pairwise comparisons of homologous protein domains corresponding to 300 protein domain families of known 3-D structures, we provide a comprehensive analysis of extent of conservation of disulphide bridges and their structural features. We report that only 54% of all the disulphide bonds compared between the homologous pairs are conserved, even if, a small fraction of the non-conserved disulphides do include cytoplasmic proteins. Also, only about one fourth of the distinct disulphides are conserved in all the members in protein families. We note that while conservation of disulphide is common in many families, disulphide bond mutations are quite prevalent. Interestingly, we note that there is no clear relationship between sequence identity between two homologous proteins and disulphide bond conservation. Our analysis on structural features at the sites where cysteines forming disulphide in one homologue are replaced by non-Cys residues show that the elimination of a disulphide in a homologue need not always result in stabilizing interactions between equivalent residues.

Conclusion

We observe that in the homologous proteins, disulphide bonds are conserved only to a modest extent. Very interestingly, we note that extent of conservation of disulphide in homologous proteins is unrelated to the overall sequence identity between homologues. The non-conserved disulphides are often associated with variable structural features that were recruited to be associated with differentiation or specialisation of protein function.

Engineered disulfide bonds increase active-site local stability and reduce catalytic activity of a cold-adapted alkaline phosphatase

Alkaline phosphatase is an extracellular enzyme that is membrane-bound in eukaryotes but resides in the periplasmic space of bacteria. It normally carries four cysteine residues that form two disulfide bonds, for instance in the APs of Escherichia coli and vertebrates. An AP variant from a Vibrio sp. has only one cysteine residue. This cysteine is second next to the nucleophilic serine in the active site. We have individually modified seven residues to cysteine that are on two loops predicted to be within a 5 A radius. Four of them formed a disulfide bond to the endogenous cysteine. Thermal stability was monitored by circular dichroism and activity measurements. Global stability was similar to the wild-type enzyme. However, a significant increase in heat-stability was observed for the disulfide-containing variants using activity as a measure, together with a large reduction in catalytic rates (k(cat)) and a general decrease in Km values. The results suggest that a high degree of mobility near the active site and in the helix carrying the endogenous cysteine is essential for full catalytic efficiency in the cold-adapted AP.

Thermodynamics of protein denatured states

Recent work on the thermodynamics of protein denatured states is providing insight into the stability of residual structure and the conformational constraints that affect the disordered states of proteins. Current data from native state hydrogen exchange and the pH dependence of protein stability indicate that residual structure can modulate the stability of the denatured state by up to 4 kcal mol(-1). NMR structural data have emphasized the role of hydrophobic clusters in stabilizing denatured state residual structures, however recent results indicate that electrostatic interactions, both favorable and unfavorable, are also important modulators of the stability of the denatured state. Thermodynamics methods that take advantage of histidine-heme ligation chemistry have also been developed to probe the conformational constraints that act on denatured states. These methods have provided insights into the role of excluded volume, chain stiffness, and loop persistence in modulating the conformational preferences of highly disordered proteins. New insights into protein folding and novel methods to manipulate protein stability are emerging from this work.

Influence of the Internal Disulfide Bridge on the Folding Pathway of the CL Antibody Domain

Disulfide bridges are one of the most important factors stabilizing the native structure of a protein. Whereas the basis for their stabilizing effect is well understood, their role in a protein folding reaction still seems to require further attention. We used the constant domain of the antibody light chain (C(L)), a representative of the ubiquitous immunoglobulin (Ig)-superfamily, to delineate the kinetic role of its single buried disulfide bridge. Independent of its redox state, the monomeric C(L) domain adopts a typical Ig-fold under native conditions and does not retain significant structural elements when unfolded. Interestingly, its folding pathway is strongly influenced by the disulfide bridge. The more stable oxidized protein folds via a highly structured on-pathway intermediate, whereas the destabilized reduced protein populates a misfolded off-pathway species on its way to the native state. In both cases, the formation of the intermediate species is shown to be independent of the isomerization state of the Tyr(141)-Pro(142) bond. Our results demonstrate that the internal disulfide bridge in an antibody domain restricts the folding pathway by bringing residues of the folding nucleus into proximity thus facilitating the way to the native state.

Mechanism of stabilization of Bacillus circulans xylanase upon the introduction of disulfide bonds

The introduction of disulfide bonds has been used as a strategy to enhance the stability of Bacillus circulans xylanase. The transition temperature of the S100C/N148C (DS1), V98C/A152C (DS2), and A1GC/G187,C188 (cXl) in comparison to the wild type was increased by 5.0, 4.1 and 3.8 degrees C, respectively. Interestingly, a combination of two disulfide bonds of DS1 and cXl (cDS1, circular disulfide 1) led to a 12 degrees C increase in the transition temperature. Importantly, an increase in the melting point and DeltaDeltaG values of the cDS1 mutant was cooperative. These results suggest that the mechanism of stabilization by disulfide bonds under irreversible denaturation condition is achieved through: (1) a change in the rate-limiting step on the denaturation pathway; (2) destabilizing the unfolded state without affecting the relative rate constants on the denaturation pathway (like cXl mutant); and (3) or combination of the two (cDS1 mutant).

Improving the stability of an antibody variable fragment by a combination of knowledge-based approaches: validation and mechanisms

Numerous approaches have been described to obtain variable fragments of antibodies (Fv or scFv) that are sufficiently stable for their applications. Here, we combined several knowledge-based methods to increase the stability of pre-existing scFvs by design. Firstly, the consensus sequence approach was used in a non-stringent way to predict a large basic set of potentially stabilizing mutations. These mutations were then prioritized by other methods of design, mainly the formation of additional hydrogen bonds, an increase in the hydrophilicity of solvent exposed residues, and previously described mutations in other antibodies. We validated this combined method with antibody mAbD1.3, directed against lysozyme. Fourteen potentially stabilizing mutations were designed and introduced into scFvD1.3 by site-directed mutagenesis, either individually or in combinations. We characterized the effects of the mutations on the thermodynamic stability of scFvD1.3 by experiments of unfolding with urea, monitored by spectrofluorometry, and tested the additivity of their effects by double-mutant cycles. We also quantified the individual contributions of the resistance to denaturation ([urea](1/2)) and cooperativity of unfolding (m) to the variations of stability and the energy of coupling between mutations by a novel approach. Most mutations (75%) were stabilizing and none was destabilizing. The progressive recombination of the mutations into the same molecule of scFvD1.3 showed that their effects were mostly additive or synergistic, provided a large overall increase in protein stability (9.1 kcal/mol), and resulted in a highly stable scFvD1.3 derivative. The mechanisms of the mutations and of their combinations involved variations in the resistance to denaturation, cooperativity of unfolding, and likely residual structures of the denatured state, which was constrained by two disulfide bonds. This combined method should be applicable to any recombinant antibody fragment, through a single step of mutagenesis.

Engineering Stabilising β-Sheet Interactions into a Conformationally Flexible Region of the Folding Transition State of Ubiquitin

Protein engineering studies suggest that the transition state for the folding of ubiquitin is highly polarised towards the N-terminal part of the sequence and involves a nucleus of residues within the beta-hairpin (residues 1-17) and main alpha-helix (residues 23-34). In contrast, the observation of small phi-values for residues in the C-terminal portion of the sequence (residues 35-76), coupled with a folding topology that results in a much higher contact order, suggests that fast folding of ubiquitin is dependent upon configurational flexibility in the C-terminal part of the polypeptide chain to ensure passage down a relatively smooth folding funnel to the native state. We show that the introduction of a small mini-hairpin motif as an extension of the native 43-50 hairpin stabilises local interactions in the C-terminal part of the sequence, resulting largely in a deceleration of the unfolding kinetics without perturbing the apparent two-state folding mechanism. However, a single-point Leu-->Phe substitution within the engineered hairpin sequence leads to the premature collapse of the denatured ensemble through the stabilisation of non-native interactions and the population of a compact intermediate. Non-linear effects in the kinetic data at low concentrations of denaturant suggest that the collapsed state, which is further stabilised in the presence of cosmotropic salts, may subsequently fold directly to the native state through a "triangular" reaction scheme involving internal rearrangement rather than unfolding and refolding.

Kinetics of Loop Formation and Breakage in the Denatured State of Iso-1-cytochrome c

The earliest events in protein folding involve the formation of simple loops. Observing the rates of loop closure under denaturing conditions can provide direct insight into the relative probability and sequence determinants for formation of loops of different sizes. The persistence of these initial contacts is equally important for efficient folding, so measurement of rates of loop breakage under denaturing conditions is also essential. We have used stopped-flow and continuous-flow methods to measure the rates of histidine-heme loop formation and breakage in the denatured state of iso-1-cytochrome c (in the presence of 3 M guanidine HCl). The data indicate that the mechanism for forming loops is a two-step process, the first step being the deprotonation of the histidine, and the second step being the binding of the histidine to the heme. This mechanism makes it possible to extract both the rate constants of formation, k(f), and breakage, k(b), of loops from the pH dependence of the observed rate constant, k(obs). To determine the dependence of k(f) and k(b) on loop size, we have carried out kinetic measurements for seven single surface histidine variants of iso-1-cytochrome c. A scaling factor (the dependence of k(f) on log[loop size]) of approximately -1.8 is observed for loop formation, similar to that observed in other systems. The magnitude of k(b) varies from 30 s(-1) to 300 s(-1), indicating that the stability of different loops varies considerably. The implications of the kinetics of loop formation and breakage in the denatured state for the mechanism of protein folding are discussed.

Prediction of disulfide connectivity from protein sequences

The difficulties in predicting disulfide connectivity from protein sequences lie in the nonlocal properties of the disulfide bridges that involve cysteine pairs at large sequence separation. Though some progress has been recently made in the prediction of disulfide connectivity, the current methods predict less than half of the disulfide patterns for the data set sharing less than 30% sequence identity. In this report, we use the support vector machines based on sequence features such as the coupling between the local sequence environments of cysteine pair, the cysteines sequence separations, and the global sequence descriptor, such as amino acid content. Our approach is able to predict 55% of the disulfide patterns of proteins with two to five disulfide bridges, which is 11-26% higher than other methods in the literature.

Quantitative measurement of protein stability from unfolding equilibria monitored with the fluorescence maximum wavelength

The fluorescence of tryptophan is used as a signal to monitor the unfolding of proteins, in particular the intensity of fluorescence and the wavelength of its maximum lambda(max). The law of the signal is linear with respect to the concentrations of the reactants for the intensity but not for lambda(max). Consequently, the stability of a protein and its variation upon mutation cannot be deduced directly from measurements made with lambda(max). Here, we established a rigorous law of the signal for lambda(max). We then compared the stability DeltaG(H(2)O) and coefficient of cooperativity m for a two-state equilibrium of unfolding, monitored with lambda(max), when the rigorous and empirical linear laws of the signal are applied. The corrective terms involve the curvature of the emission spectra at their lambda(max) and can be determined experimentally. The rigorous and empirical values of the cooperativity coefficient m are equal within the experimental error for this parameter. In contrast, the rigorous and empirical values of the stability DeltaG(H(2)O) generally differ. However, they are equal within the experimental error if the curvatures of the spectra for the native and unfolded states are identical. We validated this analysis experimentally using domain 3 of the envelope glycoprotein of the dengue virus and the single-chain variable fragment (scFv) of antibody mAbD1.3, directed against lysozyme.

Cooperativity of the oxidization of cysteines in globular proteins

Based on the 639 non-homologous proteins with 2910 cysteine-containing segments of well-resolved three-dimensional structures, a novel approach has been proposed to predict the disulfide-bonding state of cysteines in proteins by constructing a two-stage classifier combining a first global linear discriminator based on their amino acid composition and a second local support vector machine classifier. The overall prediction accuracy of this hybrid classifier for the disulfide-bonding state of cysteines in proteins has scored 84.1% and 80.1%, when measured on cysteine and protein basis using the rigorous jack-knife procedure, respectively. It shows that whether cysteines should form disulfide bonds depends not only on the global structural features of proteins but also on the local sequence environment of proteins. The result demonstrates the applicability of this novel method and provides comparable prediction performance compared with existing methods for the prediction of the oxidation states of cysteines in proteins.

Biological functions of the disulfides in bovine pancreatic deoxyribonuclease

We characterized the biochemical functions of the small nonessential (C101-C104) and the large essential (C173-C209) disulfides in bovine pancreatic (bp) DNase using alanine mutants [brDNase(C101A)] and [brDNase(C173A) and brDNase(C209A)], respectively. We also characterized the effects of an additional third disulfide [brDNase(F192C/A217C)]. Without the Ca(2+) protection, bpDNase and brDNase(C101A) were readily inactivated by trypsin, whereas brDNase(F192C/A217C) remained active. With Ca(2+), all forms of DNase, except for brDNase(C101A), were protected against trypsin. All forms of DNase, after being dissolved in 6 M guanidine-HCl, were fully reactivated by diluting into a Ca(2+)-containing buffer. However, when diluted into a Ca(2+)-free buffer, bpDNase and brDNase(C101A) remained inactive, but 60% of the bpDNase activity was restored with brDNase(F192C/A217C). When heated, bpDNase was inactivated at a transition temperature of 65 degrees C, brDNase(C101A) at 60 degrees C, and brDNase(F192C/A217C) at 73 degrees C, indicating that the small disulfide, albeit not essential for activity, is important for the structural integrity, and that the introduction of a third disulfide can further stabilize the enzyme. When pellets of brDNase(C173A) and brDNase(C209A) in inclusion bodies were dissolved in 6 M guanidine-HCl and then diluted into a Ca(2+)-containing buffer, 10%-18% of the bpDNase activity was restored, suggesting that the "essential" disulfide is not absolutely crucial for enzymatic catalysis. Owing to the structure-based sequence alignment revealing homology between the "nonessential" disulfide of bpDNase and the active-site motif of thioredoxin, we measured 39% of the thioredoxin-like activity for bpDNase based on the rate of insulin precipitation (DeltaA650nm/min). Thus, the disulfides in bpDNase not only play the role of stabilizing the protein molecule but also may engage in biological functions such as the disulfide/dithiol exchange reaction.

Structure-based design of a fluorimetric redox active peptide probe

Structure-based iterative design was used to prepare a disulfide-containing nonapeptide as a fluorimetric probe for chemical and biochemical disulfide forming and breaking reactions. The peptide is composed entirely of natural amino acids and exhibits a marked (42%) change in fluorescence between its oxidized and its reduced states. The probe is easily synthesized and highly water soluble and exhibits well-behaved kinetics on reduction with the reductant tris-carboxyethylphosphine. The reduced peptide is an excellent substrate of the enzyme quiescin-sulfhydryl oxidase and may find utility in the characterization of other disulfide oxidoreductases.

Prediction of the disulfide-bonding state of cysteines in proteins based on dipeptide composition

In this paper, a novel approach has been introduced to predict the disulfide-bonding state of cysteines in proteins by means of a linear discriminator based on their dipeptide composition. The prediction is performed with a newly enlarged dataset with 8114 cysteine-containing segments extracted from 1856 non-homologous proteins of well-resolved three-dimensional structures. The oxidation of cysteines exhibits obvious cooperativity: almost all cysteines in disulfide-bond-containing proteins are in the oxidized form. This cooperativity can be well described by protein's dipeptide composition, based on which the prediction accuracy of the oxidation form of cysteines scores as high as 89.1% and 85.2%, when measured on cysteine and protein basis using the rigorous jack-knife procedure, respectively. The result demonstrates the applicability of this new relatively simple method and provides superior prediction performance compared with existing methods for the prediction of the oxidation states of cysteines in proteins.

Conformational Properties of the Iso-1-Cytochrome c Denatured State: Dependence on Guanidine Hydrochloride Concentration

Production of seven single surface histidine variants of yeast iso-1-cytochrome c allowed measurement of the apparent pK(a), pK(a)(obs), for histidine-heme loop formation for loops of nine to 83 amino acid residues under varying denaturing conditions (2 M to 6 M guanidine hydrochloride, gdnHCl). A linear correlation between pK(a)(obs) and the log of the loop size is expected for a random coil, pK(a)(obs) proportional to k log(n), where k is a scaling factor and n is the number of monomers in the loop. For small loops of nine, 16, and 22 monomers, no dependence of pK(a)(obs) on loop size was observed at any denaturant concentration indicating effects from chain stiffness. For larger loops of 37, 56, 72, and 83 monomers, the dependence of pK(a)(obs) on log(n) was linear and the slope of that dependence decreased with increasing concentration of denaturant. The scaling factor obtained at 5 M and 6 M gdnHCl for the larger loop sizes was approximately -2.0, close to the value of -2.2 expected for a random coil with excluded volume. However, scaling factors obtained under less harsh denaturing conditions (2 M to 4.5 M gdnHCl) deviated strongly from that expected for a random coil, being in the range -3 to -4. The gdnHCl dependence of pK(a)(obs) at each loop size was also evaluated to obtain denaturant m-values. Short loops where chain stiffness dominates had similar m-values of approximately 0.25 kcal/mol M. For larger loops m-values decrease with increasing loop size indicating that less hydrophobic area is sequestered when larger loops form. It is known that the earliest events in protein folding involve the formation of simple loops. The data from these studies provide direct insight into the relative probability with which loops of different sizes will form, as well as the factors which affect loop formation.

Prediction of the bonding states of cysteines Using the support vector machines based on multiple feature vectors and cysteine state sequences

The support vector machine (SVM) method is used to predict the bonding states of cysteines. Besides using local descriptors such as the local sequences, we include global information, such as amino acid compositions and the patterns of the states of cysteines (bonded or nonbonded), or cysteine state sequences, of the proteins. We found that SVM based on local sequences or global amino acid compositions yielded similar prediction accuracies for the data set comprising 4136 cysteine-containing segments extracted from 969 nonhomologous proteins. However, the SVM method based on multiple feature vectors (combining local sequences and global amino acid compositions) significantly improves the prediction accuracy, from 80% to 86%. If coupled with cysteine state sequences, SVM based on multiple feature vectors yields 90% in overall prediction accuracy and a 0.77 Matthews correlation coefficient, around 10% and 22% higher than the corresponding values obtained by SVM based on local sequence information.

The Kinetic Pathway of Folding of Barnase

To search for folding intermediates, we have examined the folding and unfolding kinetics of wild-type barnase and four representative mutants under a wide range of conditions that span two-state and multi-state kinetics. The choice of mutants and conditions provided in-built controls for artifacts that might distort the interpretation of kinetics, such as the non-linearity of kinetic and equilibrium data with concentration of denaturant. We measured unfolding rate constants over a complete range of denaturant concentration by using by 1H/2H-exchange kinetics under conditions that favour folding, conventional stopped-flow methods at higher denaturant concentrations and continuous flow. Under conditions that favour multi-state kinetics, plots of the rate constants for unfolding against denaturant concentration fitted quantitatively to the equation for three-state kinetics, with a sigmoid component for a change of rate determining step, as did the refolding kinetics. The position of the transition state on the reaction pathway, as measured by solvent exposure (the Tanford beta value) also moved with denaturant concentration, fitting quantitatively to the same equations with a change of rate determining step. The sigmoid behaviour disappeared under conditions that favoured two-state kinetics. Those data combined with direct structural observations and simulation support a minimal reaction pathway for the folding of barnase that involves two detectable folding intermediates. The first intermediate, I(1), is the denatured state under physiological conditions, D(Phys), which has native-like topology, is lower in energy than the random-flight denatured state U and is suggested by molecular dynamics simulation of unfolding to be on-pathway. The second intermediate, I(2), is high energy, and is proven by the change in rate determining step in the unfolding kinetics to be on-pathway. The change in rate determining step in unfolding with structure or environment reflects the change in partitioning of this intermediate to products or starting materials.

Relationship between protein structures and disulfide-bonding patterns

We found that that disulfide-bonding patterns can be used to discriminate structure similarity. Our method, based on the hierarchical clustering scheme, is applicable to proteins with two or more disulfide bonds and is able to detect the structural similarities of proteins of low sequence identities (<25%). Our results show the surprisingly close relationship between disulfide-bonding patterns and proteins structures. Our findings should be useful in protein structure modeling.

Prediction of the disulfide-bonding state of cysteines in proteins at 88% accuracy

The task of predicting the cysteine-bonding state in proteins starting from the residue chain is addressed by implementing a new hybrid system that combines a neural network and a hidden Markov model (hidden neural network). Training is performed using 4136 cysteine-containing segments extracted from 969 nonhomologous proteins of well-resolved three-dimensional structure. After a 20-fold cross-validation procedure, the efficiency of the prediction scores as high as 88% and 84%, when measured on cysteine and protein basis, respectively. These results outperform previously described methods for the same task.

Unspecific hydrophobic stabilization of folding transition states

Here we present a method for determining the inference of non-native conformations in the folding of a small domain, alpha-spectrin Src homology 3 domain. This method relies on the preservation of all native interactions after Tyr/Phe exchanges in solvent-exposed, contact-free positions. Minor changes in solvent exposure and free energy of the denatured ensemble are in agreement with the reverse hydrophobic effect, as the Tyr/Phe mutations slightly change the polypeptide hydrophilic/hydrophobic balance. Interestingly, more important Gibbs energy variations are observed in the transition state ensemble (TSE). Considering the small changes induced by the H/OH replacements, the observed energy variations in the TSE are rather notable, but of a magnitude that would remain undetected under regular mutations that alter the folded structure free energy. Hydrophobic residues outside of the folding nucleus contribute to the stability of the TSE in an unspecific nonlinear manner, producing a significant acceleration of both unfolding and refolding rates, with little effect on stability. These results suggest that sectors of the protein transiently reside in non-native areas of the landscape during folding, with implications in the reading of phi values from protein engineering experiments. Contrary to previous proposals, the principle that emerges is that non-native contacts, or conformations, could be beneficial in evolution and design of some fast folding proteins.

Predicting the disulfide bonding state of cysteines using protein descriptors

Knowledge of the disulfide bonding state of the cysteines of proteins is of major interest in designing numerous molecular biology experiments, or in predicting their three-dimensional structure. Previous methods using the information gained from aligned sets of sequences have reached up to 82% of success in predicting the oxidation state of cysteines. In the present study, we assess the relative efficiency of different descriptors in predicting the cysteine disulfide bonding states. Our results suggest that the information on the residues flanking the cysteines is less informative about the disulfide bonding state than about the amino acid content of the whole protein. Using a combination of logistic functions learned with subsets of proteins homogeneous in terms of their amino acid content, we propose a simple prediction approach, starting from a single sequence, that reaches success rates close to 84%. This score can be improved by avoiding predictions regarding cysteines for which the decision is not well marked. For example, we obtain a score close to 87% correct prediction when we exclude predicting 10% of the cysteines.

Lack of definable nucleation sites in the rate-limiting transition state of barnase under native conditions

It has been shown that the burst-phase (submillisecond) intermediate of barnase, if it exists, can be only marginally more stable than the fully unfolded state at pH 6.3 and 25 degrees C. In the study reported here, no stable burst-phase intermediate could be detected, even in the presence of stabilizing salt (0.4 M Na(2)SO(4)). These results suggest that a burst-phase intermediate with even marginal stability does not exist. The absence of such an intermediate in turn suggests the need for re-examination of the rate-limiting transition state (RLTS) under native conditions, which was previously characterized by using a three-state model with a stable intermediate and protein engineering. Surprisingly, mutations throughout the structure of barnase do not significantly affect the folding rate, suggesting a lack of specific favorable interactions among the side-chains in the RLTS. This RLTS is clearly different from that previously characterized under denaturing conditions, indicating that changes take place in the RLTS under native and denaturing conditions. The occurrence of such changes is further supported by the observation that the unfolding rate constants of barnase and its mutants were divergent or convergent as a function of denaturant concentrations. Consistent with changes in the RLTS, a re-analysis of data from native-state hydrogen exchange studies has shown that the logarithm of the unfolding rate constant inflects down under low concentrations of denaturant. Here, we discuss in detail the question of whether changes in the RLTS involve a kinetically silent intermediate that occurs after the initial RLTS.