Spatial geometric arrangements of disulfide-crosslinked loops in nonplanar proteins

The method presented earlier [T. Kikuchi, G. Némethy, and H.A. Scheraga, (1986) J. Comput. Chem. 7 , 67] for the classification of patterns of the three-dimensional folding of a covalently crosslinked polypeptide chain has been extended to nonplanar proteins. The procedure described earlier was applicable only to proteins termed planar, i.e., with a connexity of the crosslinks (e.g., disulfide bonds) that can be represented in a planar diagram. The procedure described in the present work is applicable to any (planar or nonplanar) pattern of crosslinking. The classification is based on a systematic and objective method of enumeration of spatial geometric arrangements of loops (SGAL) using no information other than the location of the disulfide bonds in the amino acid sequence. Various SGAL classes correspond to the presence of different ways of mutual penetration of loops, called thrustings and entanglements. Information on SGAL classes can be of use in structural predictions of folding patterns of proteins.     

The role of conformational flexibility on protein supercharging in native electrospray ionization

Effects of covalent intramolecular bonds, either native disulfide bridges or chemical crosslinks, on ESI supercharging of proteins from aqueous solutions were investigated. Chemically modifying cytochrome c with up to seven crosslinks or ubiquitin with up to two crosslinks did not affect the average or maximum charge states of these proteins in the absence of m-nitrobenzyl alcohol (m-NBA), but the extent of supercharging induced by m-NBA increased with decreasing numbers of crosslinks. For the model random coil polypeptide reduced/alkylated RNase A, a decrease in charging with increasing m-NBA concentration attributable to reduced surface tension of the ESI droplet was observed, whereas native RNase A electrosprayed from these same solutions exhibited enhanced charging. The inverse relationship between the extent of supercharging and the number of intramolecular crosslinks for folded proteins, as well as the absence of supercharging for proteins that are random coils in aqueous solution, indicate that conformational restrictions induced by the crosslinks reduce the extent of supercharging. These results provide additional evidence that protein and protein complex supercharging from aqueous solution is primarily due to partial or significant unfolding that occurs as a result of chemical and/or thermal denaturation induced by the supercharging reagent late in the ESI droplet lifetime.     

Performance assessment of different constraining potentials in computational structure prediction for disulfide-bridged proteins

The presence of disulfide bonds in proteins has very important implications on the three-dimensional structure and folding of proteins. An adequate treatment of disulfide bonds in de-novo protein simulations is therefore very important. Here we present a computational study of a set of small disulfide-bridged proteins using an all-atom stochastic search approach and including various constraining potentials to describe the disulfide bonds. The proposed potentials can easily be implemented in any code based on all-atom force fields and employed in simulations to achieve an improved prediction of protein structure. Exploring different potential parameters and comparing the structures to those from unconstrained simulations and to experimental structures by means of a scoring function we demonstrate that the inclusion of constraining potentials improves the quality of final structures significantly. For some proteins (1KVG and 1PG1) the native conformation is visited only in simulations in presence of constraints. Overall, we found that the Morse potential has optimal performance, in particular for the β-sheet proteins.     

Dynamic redox environment-intensified disulfide bond shuffling for protein refolding in vitro: molecular simulation and experimental validation

One challenge in protein refolding is to dissociate the non-native disulfide bonds and promote the formation of native ones. In this study, we present a coarse-grained off-lattice model protein containing disulfide bonds and simulate disulfide bond shuffling during the folding of this model protein. Introduction of disulfide bonds in the model protein led to enhanced conformational stability but reduced foldability in comparison to counterpart protein without disulfide bonds. The folding trajectory suggested that the model protein retained the two-step folding mechanism in terms of hydrophobic collapse and structural rearrangement. The disulfide bonds located in the hydrophobic core were formed before the collapsing step, while the bonds located on the protein surface were formed during the rearrangement step. While a reductive environment at the initial stage of folding favored the formation of native disulfide bonds in the hydrophobic core, an oxidative environment at a later stage of folding was required for the formation of disulfide bonds at protein surface. Appling a dynamic redox environment, that is, one that changes from reductive to oxidative, intensified disulfide bond shuffling and thus resulted in improved recovery of the native conformation. The above-mentioned simulation was experimentally validated by refolding hen-egg lysozyme at different urea concentrations and oxidized glutathione/reduced glutathione (GSSG/GSH) ratios, and an optimal redox environment, in terms of the GSSG to GSH ratio, was identified. The implementation of a dynamic redox environment by tuning the GSSG/GSH ratio further improved the refolding yield of lysozyme, as predicted by molecular simulation.     

A new method for rapid characterization of the folding pathways of multidisulfide-containing proteins

Oxidative folding is a composite process that consists of both the conformational folding to the native three-dimensional structure and the regeneration of the native disulfide bonds of a protein, frequently involving over 100 disulfide intermediate species. Understanding the oxidative folding pathways of a multiple-disulfide-containing protein is a very difficult task that often requires years of devoted research due to the high complexity of the process and the very similar features of the large number of intermediates. Here we developed a method for rapidly delineating the major features of the oxidative folding pathways of a protein. The method examines the temperature dependence of the oxidative folding rate of the protein in combination with reduction pulses. Reduction pulses expose the presence of structured intermediates along the pathways. The correlation between the regeneration rate at different temperatures and the stability of the structured intermediates reveals the role that the intermediates play in determining the pathway. The method was first tested with bovine pancreatic ribonuclease A whose folding pathways were defined earlier. Then, it was explored to discern some of the major features of the folding pathways of its homologue, frog Onconase. The results suggest that the stability of the three-dimensional structure of the native protein is a major determinant of the folding rate in oxidative folding.     

Improving the accuracy of predicting disulfide connectivity by feature selection

Disulfide bonds are primary covalent cross‐links formed between two cysteine residues in the same or different protein polypeptide chains, which play important roles in the folding and stability of proteins. However, computational prediction of disulfide connectivity directly from protein primary sequences is challenging due to the nonlocal nature of disulfide bonds in the context of sequences, and the number of possible disulfide patterns grows exponentially when the number of cysteine residues increases. In the previous studies, disulfide connectivity prediction was usually performed in high‐dimensional feature space, which can cause a variety of problems in statistical learning, such as the dimension disaster, overfitting, and feature redundancy. In this study, we propose an efficient feature selection technique for analyzing the importance of each feature component. On the basis of this approach, we selected the most important features for predicting the connectivity pattern of intra‐chain disulfide bonds. Our results have shown that the high‐dimensional features contain redundant information, and the prediction performance can be further improved when these high‐dimensional features are reduced to a lower but more compact dimensional space. Our results also indicate that the global protein features contribute little to the formation and prediction of disulfide bonds, while the local sequential and structural information play important roles. All these findings provide important insights for structural studies of disulfide‐rich proteins. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Development of a novel method to populate native disulfide-bonded intermediates for structural characterization of proteins: implications for the mechanism of oxidative folding of RNase A

RNase A, a model protein for oxidative folding studies, has four native disulfide bonds. The roles of des [40-95] and des [65-72], the two native-like structured three-disulfide-bonded intermediates populated between 8 and 25 degrees C during the oxidative folding of RNase A, are well characterized. Recent work focuses on both the formation of these structured disulfide intermediates from their unstructured precursors and on the subsequent oxidation of the structured species to form the native protein. The major obstacles in this work are the very low concentration of the precursor species and the difficulty of isolating some of the structured intermediates. Here, we demonstrate a novel method that enables the native disulfide-bonded intermediates to be populated and studied regardless of whether they have stable structure and/or are present at low concentrations during the oxidative folding or reductive unfolding process. The application of this method enabled us to populate and, in turn, study the key intermediates with two native disulfide bonds on the oxidative folding pathway of RNase A; it also facilitated the isolation of des [58-110] and des [26-84], the other two native-like structured des species whose isolation had thus far not been possible.     

蛋白质中二硫键的定位及其质谱分析*

二硫键是一种常见的蛋白质翻译后修饰,对稳定蛋白质的空间结构,保持及调节其生物活性等都有着非常重要的作用。因此,确定二硫键在蛋白质中的位置是全面了解含二硫键蛋白化学结构的重要方面。在众多实验方法中,现代质谱技术因其操作简单、快速、灵敏等优点而成为分析二硫键的重要手段。本文介绍了目前主要的定位二硫键的方法,以及质谱在二硫键定位分析中的应用与进展。     

Native conformational tendencies in unfolded polypeptides: development of a novel method to assess native conformational tendencies in the reduced forms of multiple disulfide-bonded proteins

Oxidative folding is the concomitant formation of the native disulfide bonds and the native tertiary structure from the reduced and unfolded polypeptide. Of interest is the inherent conformational tendency (bias) present in the reduced polypeptide to dictate the formation of the full set of native disulfide bonds. Here, by application of a novel tool, we have been able to assess this "native conformational tendency" present in reduced and unfolded bovine pancreatic ribonuclease A (RNase A). The essence of this method lies in the ability of the oxidant [Pt(en)(2)Cl(2)](2+) (where "en" is ethylenediamine) to oxidize disulfide bonds under conditions in which both reduction and disulfide reshuffling, which are essential for rearranging non-native disulfide bonds, are extremely slow. When applied to RNase A, the method revealed little or no bias toward formation of the full native set of disulfide bonds in the fully reduced protein.     

Noncovalent cross-links in context with other structural and functional elements of proteins

Proteins are heteropolymers with evolutionary selected native sequences of residues. These native sequences code for unique and stable 3D structures indispensable for biochemical activity and for proteolysis resistance, the latter which guarantees an appropriate lifetime for the protein in the protease rich cellular environment. Cross-links between residues close in space but far in the primary structure are required to maintain the folded structure of proteins. Some of these cross-links are covalent, most frequently disulfide bonds, but the majority of the cross-links are sets of cooperative noncovalent long-range interactions. In this paper we focus on special clusters of noncovalent long-range interactions: the Stabilization Centers (SCs). The relation between the SCs and secondary structural elements as well as the relation between SCs and functionally important regions of proteins are presented to show a detailed picture of these clusters, which are believed to be primarily responsible for major aspects of protein stability.     

diSBPred: A machine learning based approach for disulfide bond prediction

The protein disulfide bond is a covalent bond that forms during post-translational modification by the oxidation of a pair of cysteines. In protein, the disulfide bond is the most frequent covalent link between amino acids after the peptide bond. It plays a significant role in three-dimensional (3D) ab initio protein structure prediction (aiPSP), stabilizing protein conformation, post-translational modification, and protein folding. In aiPSP, the location of disulfide bonds can strongly reduce the conformational space searching by imposing geometrical constraints. Existing experimental techniques for the determination of disulfide bonds are time-consuming and expensive. Thus, developing sequence-based computational methods for disulfide bond prediction becomes indispensable. This study proposed a stacking-based machine learning approach for disulfide bond prediction (diSBPred). Various useful sequence and structure-based features are extracted for effective training, including conservation profile, residue solvent accessibility, torsion angle flexibility, disorder probability, a sequential distance between cysteines, and more. The prediction of disulfide bonds is carried out in two stages: first, individual cysteines are predicted as either bonding or non-bonding; second, the cysteine-pairs are predicted as either bonding or non-bonding by including the results from cysteine bonding prediction as a feature.The examination of the relevance of the features employed in this study and the features utilized in the existing nearest neighbor algorithm (NNA) method shows that the features used in this study improve about 7.39 % in jackknife validation balanced accuracy. Moreover, for individual cysteine bonding prediction and cysteine-pair bonding prediction, diSBPred provides a 10-fold cross-validation balanced accuracy of 82.29 % and 94.20 %, respectively. Altogether, our predictor achieves an improvement of 43.25 % based on balanced accuracy compared to the existing NNA based approach. Thus, diSBPred can be utilized to annotate the cysteine bonding residues of protein sequences whose structures are unknown as well as improve the accuracy of the aiPSP method, which can further aid in experimental studies of the disulfide bond and structure determination.     

Modulation of structure and dynamics by disulfide bond formation in unfolded states

During oxidative folding, the formation of disulfide bonds has profound effects on guiding the protein folding pathway. Until now, comparatively little is known about the changes in the conformational dynamics in folding intermediates of proteins that contain only a subset of their native disulfide bonds. In this comprehensive study, we probe the conformational landscape of non-native states of lysozyme containing a single native disulfide bond utilizing nuclear magnetic resonance (NMR) spectroscopy, small-angle X-ray scattering (SAXS), circular dichroism (CD) data, and modeling approaches. The impact on conformational dynamics varies widely depending on the loop size of the single disulfide variants and deviates significantly from random coil predictions for both NMR and SAXS data. From these experiments, we conclude that the introduction of single disulfides spanning a large portion of the polypeptide chain shifts the structure and dynamics of hydrophobic core residues of the protein so that these regions exhibit levels of order comparable to the native state on the nanosecond time scale.     

Flexible Folding: Disulfide-Containing Peptides and Proteins Choose the Pathway Depending on the Environments

In the last few decades, development of novel experimental techniques, such as new types of disulfide (SS)-forming reagents and genetic and chemical technologies for synthesizing designed artificial proteins, is opening a new realm of the oxidative folding study where peptides and proteins can be folded under physiologically more relevant conditions. In this review, after a brief overview of the historical and physicochemical background of oxidative protein folding study, recently revealed folding pathways of several representative peptides and proteins are summarized, including those having two, three, or four SS bonds in the native state, as well as those with odd Cys residues or consisting of two peptide chains. Comparison of the updated pathways with those reported in the early years has revealed the flexible nature of the protein folding pathways. The significantly different pathways characterized for hen-egg white lysozyme and bovine milk α-lactalbumin, which belong to the same protein superfamily, suggest that the information of protein folding pathways, not only the native folded structure, is encoded in the amino acid sequence. The application of the flexible pathways of peptides and proteins to the engineering of folded three-dimensional structures is an interesting and important issue in the new realm of the current oxidative protein folding study.     

Topology-dependent protein folding rates analyzed by a stereochemical model

It is an experimental fact that gross topological parameters of the native structure of small proteins presenting two-state kinetics, as relative contact order chi, correlate with the logarithm of their respective folding rate constant kappa(f). However, reported results show specific cases for which the (chi,log kappa(f)) dependence does not follow the overall trend of the entire collection of experimental data. Therefore, an interesting point to be clarified is to what extent the native topology alone can explain these exceptional data. In this work, the structural determinants of the folding kinetics are investigated by means of a 27-mer lattice model, in that each native is represented by a compact self-avoiding (CSA) configuration. The hydrophobic effect and steric constraints are taken as basic ingredients of the folding mechanism, and each CSA configuration is characterized according to its composition of specific patterns (resembling basic structural elements such as loops, sheets, and helices). Our results suggest that (i) folding rate constants are largely influenced by topological details of the native structure, as configurational pattern types and their combinations, and (ii) global parameters, as the relative contact order, may not be effective to detect them. Distinct pattern types and their combinations are determinants of what we call here the "content of secondary-type" structure (sigma) of the native: high sigma implies a large kappa(f). The largest part of all CSA configurations presents a mix of distinct structural patterns, which determine the chixlog kappa(f) linear dependence: Those structures not presenting a proper chi-dependent balance of patterns have their folding kinetics affected with respect to the pretense linear correlation between chi and log kappa(f). The basic physical mechanism relating sigma and kappa(f) involves the concept of cooperativity: If the native is composed of patterns producing a spatial order rich in effective short-range contacts, a properly designed sequence undertakes a fast folding process. On the other hand, the presence of some structural patterns, such as long loops, may reduce substantially the folding performance. This fact is illustrated through natives having a very similar topology but presenting a distinct folding rate kappa(f), and by analyzing structures having the same chi but different sigma.     

A protein folding degree measure and its dependence on crystal packing, protein size, secondary structure, and domain structural class

Comparing two or more protein structures with respect to their degree of folding is common practice in structural biology despite the fact that there is no scale for a folding degree. Here we introduce a formal definition of a folding degree, capable of quantitative characterization. This enables ordering among protein chains based on their degree of folding. The folding degree of a data set of 152 representative nonhomologous proteins is then studied. We demonstrate that the variation in the folding degree seen for this data set is not due to crystallization artifacts or experimental conditions, such as resolution, refinement protocol, pH, or temperature. A good linear relationship is observed between the folding degree and the percentages of secondary structures in the protein. The folding degree is able to account for the small changes produced in the structure due to crystal packing and temperature. Automating the classification of proteins into their respective structural domain classes, namely mainly-alpha, mainly-beta, and alpha-beta, is also possible.     

Energy minimization of rigid-geometry polypeptides with exactly closed disulfide loops

A method has been developed for minimizing the energy of a polypeptide with rigid geometry while keeping all disulfide loops closed exactly. Exact closure of disulfide loops implies that some dihedral angles become implicit functions of the remaining dihedral angles in the polypeptide; the efficacy of the method is related to the manner in which the implicitly defined dihedral angles are chosen. The method has been used to find minimum-energy conformations of bovine pancreatic trypsin inhibitor, ribonuclease A, crambin, the defensin HNP3 dimer, and ω-conotoxin. For the first two proteins, the starting conformations for energy minimization had been derived previously from crystal structures using pseudopotentials to keep the disulfide loops almost closed. Starting conformations for the remaining three proteins were derived from their crystal or NMR structures by similar procedures. In all cases, the energy-minimized structures had a significantly and, in some cases, substantially, lower energy than the starting structures. The RMS deviations between the exactly closed energy- minimized structures and the crystal or NMR structures from which they were derived ranged from 0.9 Å to 1.9 Å, suggesting that the computed structures can serve as “regularized” native structures for these proteins. The energy of a ribonuclease derivative lacking the 65–72 disulfide bridge was minimized using the procedure; the result showed that this derivative has a low-energy structure with a conformation very close to that of native ribonuclease, and is consistent with its postulated role in the folding of ribonuclease. These results offer strong support for the validity of the rigid-geometry model in the studies of the conformational energy of proteins. © 1997 by John Wiley & Sons, Inc.     

A statistical model for predicting protein folding rates from amino acid sequence with structural class information

Prediction of protein folding rates from amino acid sequences is one of the most important challenges in molecular biology. In this work, I have related the protein folding rates with physical-chemical, energetic and conformational properties of amino acid residues. I found that the classification of proteins into different structural classes shows an excellent correlation between amino acid properties and folding rates of two- and three-state proteins, indicating the importance of native state topology in determining the protein folding rates. I have formulated a simple linear regression model for predicting the protein folding rates from amino acid sequences along with structural class information and obtained an excellent agreement between predicted and experimentally observed folding rates of proteins; the correlation coefficients are 0.99, 0.96 and 0.95, respectively, for all-alpha, all-beta and mixed class proteins. This is the first available method, which is capable of predicting the protein folding rates just from the amino acid sequence with the aid of generic amino acid properties and structural class information.     

Molecular chaperones, folding catalysts, and the recovery of active recombinant proteins fromE. coli

The high-level expression of recombinant gene products in the gramnegative bacteriumEscherichia coli often results in the misfolding of the protein of interest and its subsequent degradation by cellular proteases or its deposition into biologically inactive aggregates known as inclusion bodies. It has recently become clear that in vivo protein folding is an energy-dependent process mediated by two classes of folding modulators. Molecular chaperones, such as the DnaK-DnaJ-GrpE and GroEL-GroES systems, suppress off-pathway aggregation reactions and facilitate proper folding through ATP-coordinated cycles of binding and release of folding intermediates. On the other hand, folding catalysts (foldases) accelerate rate-limiting steps along the protein folding pathway such as thecis/trans isomerization of peptidyl-prolyl bonds and the formation and reshuffling of disulfide bridges. Manipulating the cytoplasmic folding environment by increasing the intracellular concentration of all or specific folding modulators, or by inactivating genes encoding these proteins, holds great promise in facilitating the production and purification of heterologous proteins. Purified folding modulators and artificial systems that mimic their mode of action have also proven useful in improving the in vitro refolding yields of chemically denatured polypeptides. This review examines the usefulness and limitations of molecular chaperones and folding catalysts in both in vivo and in vitro folding processes.     

A small-molecule catalyst of protein folding in vitro and in vivo

BACKGROUND: The formation of native disulfide bonds between cysteine residues often limits the rate and yield of protein folding. The enzyme protein disulfide isomerase (PDI) catalyzes the interchange of disulfide bonds in substrate proteins. The two -Cys-Gly-His-Cys- active sites of PDI provide a thiol that has a low pKa value and a disulfide bond of high reduction potential (Eo'). RESULTS: A synthetic small-molecule dithiol, (+/-)-trans-1,2-bis(2-mercaptoacetamido)cyclohexane (BMC), has a pKa value of 8.3 and an Eo' value of -0.24 V. These values are similar to those of the PDI active sites. BMC catalyzes the activation of scrambled ribonuclease A, an inactive enzyme with non-native disulfide bonds, and doubles the yield of active enzyme. A monothiol analog of BMC, N-methylmercaptoacetamide, is a less efficient catalyst than BMC. BMC in the growth medium of Saccharomyces cerevisiae cells increases by > threefold the heterologous secretion of Schizosaccharomyces pombe acid phosphatase, which has eight disulfide bonds. This effect is similar to that from the overproduction of PDI in the S. cerevisiae cells, indicating that BMC, like PDI, can catalyze protein folding in vivo. CONCLUSIONS: A small-molecule dithiol with a low thiol pKa value and high disulfide Eo' value can mimic PDI by catalyzing the formation of native disulfide bonds in proteins, both in vitro and in vivo.