Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2015 Sep 8;112(36):11241-6.
doi: 10.1073/pnas.1503909112. Epub 2015 Aug 21.

Protein folding guides disulfide bond formation

Affiliations

Protein folding guides disulfide bond formation

Meng Qin et al. Proc Natl Acad Sci U S A. .

Abstract

The Anfinsen principle that the protein sequence uniquely determines its structure is based on experiments on oxidative refolding of a protein with disulfide bonds. The problem of how protein folding drives disulfide bond formation is poorly understood. Here, we have solved this long-standing problem by creating a general method for implementing the chemistry of disulfide bond formation and rupture in coarse-grained molecular simulations. As a case study, we investigate the oxidative folding of bovine pancreatic trypsin inhibitor (BPTI). After confirming the experimental findings that the multiple routes to the folded state contain a network of states dominated by native disulfides, we show that the entropically unfavorable native single disulfide [14-38] between Cys14 and Cys38 forms only after polypeptide chain collapse and complete structuring of the central core of the protein containing an antiparallel β-sheet. Subsequent assembly, resulting in native two-disulfide bonds and the folded state, involves substantial unfolding of the protein and transient population of nonnative structures. The rate of [14-38] formation increases as the β-sheet stability increases. The flux to the native state, through a network of kinetically connected native-like intermediates, changes dramatically by altering the redox conditions. Disulfide bond formation between Cys residues not present in the native state are relevant only on the time scale of collapse of BPTI. The finding that formation of specific collapsed native-like structures guides efficient folding is applicable to a broad class of single-domain proteins, including enzyme-catalyzed disulfide proteins.

Keywords: disulfide proteins; early collapse; enzyme-catalyzed folding; native-like interactions; nonnative interactions.

PubMed Disclaimer

Conflict of interest statement

The authors declare no conflict of interest.

Figures

Fig. 1.
Fig. 1.
(A) Ribbon diagram of the native structure of the 58 residue BPTI containing three disulfide bonds (marked in yellow) between residues Cys5 and Cys55 [5–55], Cys14 and Cys38 [14–38], and Cys30 and Cys51 [30–51], respectively. The antiparallel β-sheet is in red. (B) Simplified representation of the secondary tructure of BPTI and the three native disulfide bonds in BPTI. (C) Variables βO and βR mimicking the redox conditions. Small βO (βR) represent strongly oxidizing (reducing) condition. The star with βO = 1.0 and βR = 1.5 is used in most of the simulations. These values are a mixture of mildly oxidizing and reducing condition. (D) Distribution of fraction of native contacts obtained from high-temperature simulations.
Fig. S1.
Fig. S1.
The mean values of the three structural parameters dα, θαk, and nα obtained from equilibrium simulations at kBT = 0.9ϵ for BPTI containing with different combination of S–S bonds. Here, dα is the distance between the two cysteines in native state, θαk (k=1,2) is the orientation angle defined in Fig. S2, and nα is the number of residues within a spherical shell with radius R (Fig. S3). The label α = 1, 2, and 3 refer to the three disulfide bonds [5–55], [14–38], and [30–51], respectively. For example, for [5–55] (Left), these values are obtained from the simulations of BPTI mutants containing [5–55], N′, NSHSH, and N.
Fig. S2.
Fig. S2.
Definition of the orientation angle θ1 and θ2. The indices i − 1, i, i + 1 and j − 1, j, i + 1 are the two different peptide segments such that i and j are the two cysteines that can form a disulfide bond in the native state; θ1 is the angle between the vector connecting the cysteine residues i and j (i < j) and the covalent bond connecting Cysi and the neighboring Cαi1. A similar definition holds for θ2.
Fig. S3.
Fig. S3.
The definition of the number of residues (nα) within a spherical shell with radius R, drawn from O, the center of the αth S–S bond.
Fig. 2.
Fig. 2.
The folding pathway, represented as a network of native-like states, connecting the fully unfolded state (R) to the folded (N) state with three disulfide bonds. For βO = 1.0 and βR = 1.5, the early event produces predominantly [14–38]. Subsequently, there is a bifurcation in the pathway with this intermediate rearranging to [5–55] and [30–51]. The percentages indicate the dominant route to N from R. Representative structures of all of the relevant states are shown. The numbers in parentheses were obtained from simulations with βO = 2.0 and βR = 4.0, which mimics the redox conditions used by Weissman and Kim (24). In red are the fluxes through the various native intermediates to the native state obtained from simulations that consider native and nonnative disulfide bond formation. The qualitative agreement between the two simulations is striking.
Fig. 3.
Fig. 3.
(A) Dependence of the mean first passage time (black line) for forming [14–38] as the oxidizing condition is changed from being strong (small βO) to weak (large βO). The red line gives the yield of [14–38]. (B) Time-dependent decay, Puα(t), of the three native single disulfide species. (Inset) lnPuα(t); the lines are linear fits.
Fig. S4.
Fig. S4.
(A) A trajectory showing the time-dependent changes during the [14–38] → [30–51] transition. The distance d1438 increases (red) abruptly around t 9× 103τ and is accompanied by a decrease in d3051 (blue). During the transition, Q remains low (scale is on the right), which is also illustrated by the snapshot below. (B) Same as A except the trajectories and conformations are for the [14–38] → [5–55] transition. Note that the scales for the distance as well as Q are different in A and B. (C) Distribution of the radius of gyration for the ensemble of R′ structures sampled during the [14–38] → [30–51] transition, an example of which is shown in A. (D) Same as C except the distribution is for the [14–38] → [5–55] transition.
Fig. S5.
Fig. S5.
Distribution of Rg using an ensemble of unfolded conformations generated at high temperature. The set of initial conformations with Rg> 20 Å and fraction of native contacts Q< 0.1 (Fig. 1D) is used to initiate oxidative refolding. The initial conformations are devoid of any disulfide bonds.
Fig. 4.
Fig. 4.
(A) An example of one of 2,000 oxidative folding trajectories showing the route to the folded state from a fully reduced starting conformation with low Q. The blue curve shows the decrease in the radius of gyration, Rg (scale is on the right). The gray lines show formation of various disulfide species labeled on the left. Snapshots (a–h) show some of the conformations sampled in the trajectory. (B) Plot of the distance between the distance between residues 14 and 38, d1438, as a function of the fraction, Qβ, of contacts in the β-hairpin shows that the hairpin (formed between β-strands from Ile18 to Asn24, and from Leu29 to Tyr35) is fully structured before [14–38] formation. (C) Distribution of Qβ when d1438 = 5.7 Å (the distance at which [S–S] bond forms) for the first time in a folding trajectory.
Fig. S6.
Fig. S6.
A second folding trajectory is shown to illustrate that except for changes in the value of the first passage time, the behavior is similar to the one shown in Fig. 4A. Some of the sampled conformations (a–h) are explicitly shown. This figure, together with the one in Fig. 4A, shows that structures similar to the native state of BPTI are populated in the folding process relatively early.
Fig. 5.
Fig. 5.
(A) Illustration of the dynamics of rearrangement from N′ to NSHSH using distance between the different Cys residues. The colors are illustrated in the figure. In this transition, BPTI samples both the nonnative species [30–51, 5–14] (structure on the left) and [30–51, 5–38] (conformation on the right). (B) Same as A except this trajectory describes the N* NSHSH transition.
Fig. S7.
Fig. S7.
Distributions of the radius of gyration (Rg) for the three native two-disulfide intermediates calculated from simulations emphasizing only native-like interactions. The distributions are narrow and are roughly peaked at the value of Rg corresponding to the native state, thus confirming the near native structures in these intermediates.
Fig. 6.
Fig. 6.
Time for forming [14–38] relative to the wild-type as function of λ (defined in Eqs. S1 and S2). The stability of the β-hairpin increases (decreases) as λ increases (decreases). The relative times are given for two variants. One of them (black line) is for the WT (λ = 1.0) and the other blue is for a pseudo mutant, in which disulfide bonds other than [14–38] cannot form. The red curves show the yield of [14–38] for the WT.
Fig. S8.
Fig. S8.
Kinetics of formation of three single-disulfide bonds for type I model (A), type II model (B), and type III model (C). The three models are described in detail in the main text. The mean first-passage time values, quoted in the main text, are obtained using single exponential fits shown as solid lines in A–C.
Fig. S9.
Fig. S9.
Folding pathways for an in silico mutant created by setting ESS=0 (Eq. S1) for only the [14–38] disulfide bond. Although the mutant does fold clearly, the resulting pathways bear no resemblance to the ones found for the wild type (Fig. 2).
Fig. S10.
Fig. S10.
Kinetics of formation of [14–38] (black) and the entropically favored nonnative [5–14] (red). We also show the decay of Rg in blue (scale on the right). The results show that both the native [14–38] and the nonnative [5–14] only form after considerable compaction of the polypeptide chain. In addition, even with nonnative interactions [14–38] forms slightly ahead of [5–14] with substantially greater yield. The simulations were performed using βO=1.0 and βR=1.5.

References

    1. Anfinsen CB. Principles that govern the folding of protein chains. Science. 1973;181(4096):223–230. - PubMed
    1. Anfinsen CB, Scheraga HA. Experimental and theoretical aspects of protein folding. Adv Protein Chem. 1975;29:205–300. - PubMed
    1. Schröder M, Kaufman RJ. The mammalian unfolded protein response. Annu Rev Biochem. 2005;74:739–789. - PubMed
    1. Welker E, Wedemeyer WJ, Narayan M, Scheraga HA. Coupling of conformational folding and disulfide-bond reactions in oxidative folding of proteins. Biochemistry. 2001;40(31):9059–9064. - PubMed
    1. Camacho CJ, Thirumalai D. Theoretical predictions of folding pathways by using the proximity rule, with applications to bovine pancreatic trypsin inhibitor. Proc Natl Acad Sci USA. 1995;92(5):1277–1281. - PMC - PubMed

Publication types

LinkOut - more resources