Next Article in Journal
A Preliminary Study of Cu Exposure Effects upon Alzheimer’s Amyloid Pathology
Next Article in Special Issue
The Future of Protein Secondary Structure Prediction Was Invented by Oleg Ptitsyn
Previous Article in Journal
Hydroxylamine Analogue of Agmatine: Magic Bullet for Arginine Decarboxylase
Previous Article in Special Issue
How A Novel Scientific Concept Was Coined the “Molten Globule State”
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Biomolecules
Volume 10
Issue 3
10.3390/biom10030407
biomolecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

The Molten Globule, and Two-State vs. Non-Two-State Folding of Globular Proteins

by
Kunihiro Kuwajima
1,2
1
Department of Physics, School of Science, the University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan
2
School of Computational Sciences, Korea Institute for Advanced Study (KIAS), Seoul 02455, Korea
Biomolecules 2020, 10(3), 407; https://doi.org/10.3390/biom10030407
Submission received: 19 February 2020 / Revised: 3 March 2020 / Accepted: 6 March 2020 / Published: 6 March 2020
(This article belongs to the Special Issue Protein Folding and Quality Control Mechanisms - In Memory of Prof. Oleg B. Ptitsyn (1929-1999))
Browse Figures
Review Reports Versions Notes

Abstract

:
From experimental studies of protein folding, it is now clear that there are two types of folding behavior, i.e., two-state folding and non-two-state folding, and understanding the relationships between these apparently different folding behaviors is essential for fully elucidating the molecular mechanisms of protein folding. This article describes how the presence of the two types of folding behavior has been confirmed experimentally, and discusses the relationships between the two-state and the non-two-state folding reactions, on the basis of available data on the correlations of the folding rate constant with various structure-based properties, which are determined primarily by the backbone topology of proteins. Finally, a two-stage hierarchical model is proposed as a general mechanism of protein folding. In this model, protein folding occurs in a hierarchical manner, reflecting the hierarchy of the native three-dimensional structure, as embodied in the case of non-two-state folding with an accumulation of the molten globule state as a folding intermediate. The two-state folding is thus merely a simplified version of the hierarchical folding caused either by an alteration in the rate-limiting step of folding or by destabilization of the intermediate.

1. Introduction

Elucidation of the molecular mechanisms of protein folding is a fundamental problem in biophysics and molecular biology, although almost 60 years have elapsed since Anfinsen and his coworkers discovered the genetic control of the native tertiary structure of proteins [1]. In the interim, there have been many experimental studies on protein folding (refs. [2,3,4] and references cited therein). From these studies, it is now clear that there are two types of folding behavior, i.e., non-two-state folding, involving at least one folding intermediate, and two-state folding without any detectable intermediate during kinetic refolding from the fully unfolded (U) state to the native (N) state [5,6]. Th thermodynamic states, like N, U, and the folding intermediate, of a protein surrounded by solvent water are large ensembles of many microstates, but the different thermodynamic states are well separated from each other by a free-energy barrier higher than the thermal fluctuation energy. Understanding the relationships between the two-state and non-two-state folding may thus be essential in order to fully elucidate the molecular mechanisms of protein folding. This article will discuss these relationships on the basis of available data on the correlations of the rate constant of folding with various structure-based properties, which are determined primarily by the backbone structure (topology) of proteins.
If we assume that a protein folds into the N state by a random search of all possible conformations, the time required for successful folding would be unrealistically long for even a small protein. This is known as Levinthal’s paradox [7] (see also [2,8]). Traditionally, the presence of a specific pathway of folding was thought to provide a specific and unique way to get from U to N in a reasonable amount of time [7], and hence, detection and characterization of specific intermediates, present along the pathway of folding, were basic approaches for elucidating the mechanisms of folding. The molten globule (MG) state was proposed as such a specific intermediate of folding [9], and experimentally, the MG state was observed as an equilibrium unfolding intermediate in certain proteins, and also as a transient folding intermediate, which accumulates at an early stage of kinetic refolding [10,11]. The transient MG state was observed not only in the proteins that show the equilibrium MG state but also in the proteins that show the two-state unfolding transition at equilibrium. Oleg B. Ptitsyn, our colleagues, and I thus conducted a cooperative research project supported by the Human Frontier Science Program (HFSP) for three years from 1993. The title of the project was “Role of the Molten Globule State in the Folding of Globular Protein.” In the next section (Section 2) of this article, the MG state in protein folding is thus described.
The two-state folding behavior, in which both equilibria and kinetics of folding/unfolding occur only between N and U, was discovered for a number of small globular proteins, usually with less than 100 amino acid residues [12,13], and this discovery changed the focus of many researchers in regard to the two-state proteins. These initial efforts suggested that the two-state folding, which required no stable intermediates for complete and successful folding to the N state, was simpler than the non-two-state folding. Theoretical and computational studies of protein folding have proposed an energy landscape theory (a funnel model) of folding, and the two-state folding behavior has been interpreted reasonably well by this model [14]. In Section 3, the two-state vs. non-two-state folding of globular proteins will be described. This section mainly deals with the significance of a burst-phase intermediate, which is accumulated within the dead time of kinetic folding experiments, because the burst-phase intermediate has often been observed in non-two-state proteins [11]. Sub-millisecond continuous-flow techniques [15] and more recently in-gel folding techniques [16] have revealed that the burst-phase intermediate is a real folding intermediate separated from the U state by a free-energy barrier, so that there are two-types of globular proteins in terms of the folding behavior, i.e., two-state proteins and non-two-state proteins.
In Section 4, the relationships between the two-state and non-two-state folding reactions will be discussed. The logarithmic rate constant of folding, ln(kf), shows significant correlations with structure-based properties, which are determined by the backbone structure. Recently, various structure-based properties that correlate with ln(kf) were reported, and the correlations of ln(kf) with the structure-based properties were found to be quite similar between two-state and non-two-state proteins, indicating that the physical mechanisms behind two-state and non-two-state folding are essentially identical (see Section 4). From these results, it is concluded that the general mechanism of protein folding can be described using a two-stage hierarchical model, with the two-state folding being merely a simplified version of the hierarchical folding.

2. Role of The Molten Globule State in Protein Folding

The molten globule (MG) state is an intermediate ensemble between the N and the U states of globular proteins. It has the following characteristics: (1) the presence of a substantial amount of native-like secondary structure; (2) the absence of most of the specific tertiary structure associated with the tight packing of side chains; (3) the presence of a native-like backbone fold (topology), and the compactness of the overall shape of the molecule, with a radius only 10%–30% larger than that in the N state; and (4) the presence of a loosely organized hydrophobic core, and the heterogeneity of the three-dimensional structure, in which certain subdomains of the molecule are more organized than others [10,11]. For more recent and comprehensive information on the MG state, see [17]. The loosely organized hydrophobic core is often characterized by binding experiments using a hydrophobic fluorescent probe, 1-anilino-8-naphthalene sulfonate (ANS) [18]. The equilibrium MG state has been observed for a number of globular proteins under mildly denaturing conditions (e.g., at acidic pH, at a high pressure or in the presence of an intermediate concentration of a strong denaturant such as guanidinium chloride (GdmCl) or urea) [10,11,19]. The role of the MG state in kinetic folding of globular proteins was an intriguing issue when our project started [20], and Oleg proposed that the MG state is a general intermediate in protein folding [9]. In fact, transient kinetic intermediates were detected and characterized for a number of globular proteins, including α-lactalbumin [21], lysozyme [21,22], cytochrome c [23], ribonuclease A [24] and apo-myoglobin [25], by a kinetic circular dichroism (CD) technique and a pulsed hydrogen/deuterium (H/D)-exchange method combined with two dimensional (2D) NMR spectroscopy. Close similarity between the equilibrium MG state and the kinetic folding intermediate thus characterized was demonstrated for certain proteins by coincidence of the equilibrium unfolding transition curve of the MG state and the pre-equilibrium unfolding transition curve of the kinetic intermediate [26,27], and by close similarity of the H/D-exchange protection profile between the equilibrium and kinetic intermediates [25].
To further characterize the equilibrium and kinetic MG states, Ptitsyn’s group and the Japanese members of our cooperative project carried out joint experiments, and Gennady V. Semisotnov (Institute of Protein Research, Russia) often visited Japan to participate in these efforts. Hiroshi Kihara (Kansai Medical Univ.), Yoshiyuki Amemiya (Univ. Tokyo), Kazumoto Kimura (Dokkyo Univ.), and students of our laboratories at that time were also involved in these experiments. We utilized a synchrotron radiation facility at the High Energy Accelerator Research Organization, Tsukuba, Japan to characterize the MG state of proteins by a small angle X-ray scattering (SAXS) technique, which gives us the information about the size and shape of a protein molecule in solution [29]. Figure 1 shows small-angle X-ray scattering (SAXS) patterns and Kratky plots of native and fully unfolded carbonic anhydrase [28]. The SAXS patterns are represented by the scattering intensity, I(h), as a function of the scattering vector, h, where h = ( 4 π sin θ ) / λ ; 2θ is the scattering angle, and λ is the wavelength of the X-ray. The Kratky plots are given by I(h) × h2 as a function of h [29]. From Figure 1, we can see large differences in the scattering intensity and the shape of scattering curves between the globular N state and the coil-like U state of the protein [28]. Such differences in the scattering properties encouraged us to investigate the time-resolved SAXS during kinetic refolding. If the early kinetic folding intermediate is identical to the MG state, it would be possible to directly observe a compact shape of the protein molecule, an important characteristic of the MG state, at an early stage of kinetic folding. We investigated the kinetic refolding reactions of β-lactoglobulin and α-lactalbumin induced by a denaturant (urea or GdmCl) concentration jump using a stopped-flow SAXS apparatus [30,31]. Both the proteins formed a compact globular structure with a radius of gyration (Rg) identical to that of the MG state, although the intermediate of β-lactoglobulin exhibited a non-native α-helical structure [30,32].
Figure 2 shows time-dependent changes in the Rg value during refolding of α-lactalbumin and the Kratky plots of the protein in the early stages of refolding [31]. Changes in the Rg2 values were fitted to a single-exponential function, because Rg2 (instead of Rg) has a linear dependence upon the fractional populations of the individual states [31,33]. The rate constant thus obtained was 0.49 (±0.07) s−1, which is in agreement with the rate measured by other spectroscopic techniques. The Rg value after complete refolding was 15.5 (±0.1) Å, and the Rg value at zero time of refolding obtained by extrapolation of the fitting curve to zero time was 17.5 (±0.2) Å. Because this protein forms the transient folding intermediate within the dead time (~7 ms) of the stopped-flow measurement, the folding intermediate has an Rg of 17.5 (±0.2) Å, which is in good agreement with the known Rg value of the equilibrium MG state [31]. The Kratky plot obtained by averaging the scattering pattern within 10–20 ms of refolding is coincident with that of the equilibrium MG state, indicating that both the size and shape of the transient folding intermediate were identical to those of the equilibrium MG state.
In the SAXS experiments of Figure 2, we used a 2D charge-coupled device (CCD)-based X-ray detector with a beryllium-windowed X-ray image intensifier, by which the S/N ratio of the scattering data was dramatically improved as compared with the data obtained by a one-dimensional detector [31]. This 2D CCD-based detector was developed by Amemiya and coworkers [34], and detailed procedures for data correction and the analysis software were described by Ito et al. [35]. The 2D CCD-based detector was used for nearly a decade in the early 2000s to characterize compact folding intermediates of other globular proteins [36,37,38,39,40,41,42], but at present, a more sensitive PILATUS pixel detector is widely used for the SAXS experiments [43,44,45].

3. Two-State vs. Non-Two-State Folding of Globular Proteins

In 1991, just two years before the start of the above-mentioned project supported by HFSP, Jackson and Fersht [12] in Cambridge reported that barley chymotrypsin inhibitor 2 (CI2), a small monomeric protein of 83 residues (the first 18 residues are unstructured), appears to be a rare example in which both equilibria and kinetics are described by a two-state model. As far as I know, this was the first report on observation of the two-state folding of globular proteins. The criteria for the two-state folding are (i) observation of a single equilibrium unfolding transition as a function of an external parameter such as denaturant (GdmCl or urea) concentration or temperature; (ii) linear free-energy relationships of the equilibrium unfolding free energy and the activation free energies for kinetic folding and unfolding with respect to the denaturant concentration; (iii) single-exponential kinetics for both folding and unfolding with the observation of the full change of a structural probe, fluorescence, CD ellipticity or any other index used to monitor the conformational transition, between U and N; (iv) consistency between the equilibrium and the kinetic parameters of folding/unfolding; and (v) agreement between the van’t Hoff enthalpy and the calorimetric enthalpy of equilibrium thermal unfolding. Any deviations from the above criteria suggest the presence of an intermediate between N and U [2]. However, when a proline peptide bond undergoes a slow cis-trans isomerization in the U state, the folding kinetics become complex in spite of the absence of any intermediates between N and U [2,46], and such a case may also be classified as a two-state folding.
Although Jackson and Fersht [12] reported that CI2 was a rare example, later many small proteins, usually with less than 100 amino-acid residues, were shown to fold with a similar simple two-state mechanism, and in 1998, Jackson reported more than 20 examples of two-state folders, which show wide variation in folding rates from microseconds to seconds [13]. Interestingly, the kinetic refolding reaction of protein L, a typical two-state protein, measured by the time-resolved SAXS has shown that the chain collapse occurs concomitantly with the formation of the N state, although there is some disagreement between SAXS and single-molecule Förster resonance energy transfer experiments [33,47]. My colleagues in KIAS (the Korea Institute for Advanced Study) and I recently constructed a standardized protein folding database (PFDB) with temperature correction (http://lee.kias.re.kr/~bala/PFDB/), in which the ln(kf) values of all listed proteins were calculated at the standard temperature (25 °C) [48]. Among the total 141 proteins registered in PFDB, 89 are two-state proteins, and hence two-state proteins are no longer rare examples. However, this does not necessarily mean that a majority of naturally occurring proteins are of the two-state type. Because two-state proteins are apparently simpler than non-two-state proteins, researchers in the protein folding field tend to investigate two-state proteins, which would account for the large number of two-state proteins registered in PFDB. Among the 89 two-state proteins registered in PFDB, 53 are composed of an isolated domain of a multi-domain protein. The average chain length of the two-state proteins registered in PFDB is 92.6, while the average chain length of the non-two-state proteins is 139.4 [48].
The folding intermediate appears not to be required for the complete folding of two-state proteins, which has raised questions as to the significance of the folding intermediates previously detected and characterized in non-two-state proteins [49,50]. The formation of an early kinetic intermediate often occurs too quickly to measure directly by a conventional rapid reaction technique such as the stopped-flow method, and the process occurs in a burst phase within the dead time of the measurement [11,51,52]. This situation raises a question in regard to the burst-phase—namely, whether an observed burst-phase signal reflects a real folding event. To answer this question, the burst-phase fluorescence or CD signals at an early stage of kinetic refolding of cytochrome c and ribonuclease A were thus investigated by Sosnick et al. [53,54,55,56]. The refolding reactions were induced by stopped-flow concentration jumps of GdmCl, and the burst-phase signals thus obtained were compared with the signals of non-folding polypeptides: fragments 1–80 and 1–65 for cytochrome c and a disulfide broken form for ribonuclease A. Rather surprisingly, the signals of the intact proteins and the non-folding polypeptides were coincident with each other. It was thus concluded that the sub-millisecond burst-phase signals may not reflect the fast formation of the folding intermediate, but rather may reflect some solvent-dependent contraction/collapse of the U state, since the water solution at a low concentration of GdmCl is a poor solvent for the polypeptides [53,54,55], and the initial barrier mechanism in protein folding was proposed [56]. Hydrophobic amino-acid residues, usually buried inside the molecule in the N state, are exposed to solvent water in the U state, and hence the water solution is thought to be a poor solvent for the unfolded polypeptide [57].
Roder’s group [58,59] and Takahashi’s group [39,60] thus investigated the sub-millisecond processes of folding reactions of cytochrome c and ribonuclease A by continuous-flow techniques, which have a higher time resolution than the conventional stopped-flow method. In 1998, Shastry and Roder [58] found that the kinetics of the cytochrome c folding measured by tryptophan (Trp 59) fluorescence exhibits a major exponential process with a time constant of ~60 μs, independent of initial conditions and heme ligation state, indicating that a common free energy barrier is encountered during the formation of an early compact intermediate. Therefore, there are at least two distinct stages in the cytochrome c folding, with the compact folding intermediate accumulating in the first stage, followed by the rate-limiting formation of the N state. Interestingly, Qui et al. [61] found by the laser temperature-jump method that not only full-length cytochrome c but also non-folding peptides 1–65 and 1–80 displayed similar exponential decays, indicating that the free-energy barrier to collapse was similar in the folding and non-folding sequences. Because the non-folding peptides contain significant portions of the natural cytochrome c sequence, it might be possible that such sequences can form a similar compact structure [62]. For ribonuclease A, Welker et al. [59] reported that the very fast folding U state (Uvf), which contains only native proline isomers, exhibits an exponential decay in fluorescence with a time constant of ~80 μs, whereas the equilibrium U state (Ueq), which contains nonnative proline isomers, did not show this exponential decay, indicating that the isomerization state of proline peptide bonds can have a major impact on the early stages of the kinetic folding. It may be noted that the burst phase previously observed in ribonuclease A by Houry et al. [52] using a stopped-flow double-jump apparatus was for Uvf and not for Ueq. Kimura et al. [39] have shown that the sub-millisecond folding intermediate of intact ribonuclease A strongly binds to ANS and has a compact structure, characteristics of the MG state, but that disulfide-reduced ribonuclease A shows neither any significant binding to ANS nor such compaction of the structure, indicating the difference in the sub-millisecond conformations of ribonuclease A and its reduced form.
Recently, Okabe et al. [16] employed an elegant technique to delineate the solution burst-phase events of folding by encapsulating the proteins in silica gels. Within a porous silica gel with a large water content, large scale protein motions are dramatically slowed, but the in-gel N state can retain its solution properties [63,64,65,66]. In particular, the folding of a protein in such a gel is dramatically slowed, and this makes it possible to directly observe the entire folding process without substantially perturbing the folding pathway [67,68,69,70]. Figure 3 shows (A) representative far-UV CD spectra of cytochrome c during the refolding and (B) a typical kinetic refolding curve measured by the CD ellipticity at 220 nm for the initial resolvable phase of the in-gel refolding; the protein that was initially acid unfolded at pH 1.8 was refolded at pH 4.5 and 25°C in a silica gel containing 80% (w/w) water [16]. The CD spectroscopy in the far-UV region provides a very effective tool to characterize the secondary structure of a protein, and in particular, the CD ellipticity at 220 nm gives a measure of the α-helical content of the protein [71]. The CD spectra in the N and U states of cytochrome c in the gel are nearly identical to the corresponding spectra in solution, indicating that the conformations in the N and U states are not influenced by the gel matrix (Figure 3A). As shown in Figure 3B, the initial phase of the in-gel folding exhibits an exponential decay with a time constant of ~30 s. Because the fastest phase of the protein in solution at pH 4.5 and 22 °C has a time constant of 57 μs [58], the refolding rate was at least 5.3 × 105 times reduced in the in-gel folding compared with the refolding in solution. The CD spectrum of the early kinetic intermediate observed at 3 min, 6 times the time constant (~30 s), in the in-gel folding (Figure 3A) is very similar to the CD spectrum at 400 μs, seven times the time constant (57 μs), in the free solution refolding reported by Akiyama et al. [60], suggesting that the same transient folding intermediate accumulates at the initial stage in the in-gel and the free solution refolding reactions. It may also be noted that the acid compact state of cytochrome c stabilized by KCl, known also as the MG state [72,73], is an equilibrium counterpart of a late folding intermediate separated from U by the rate-limiting free-energy barrier [74], and hence it is different from the intermediate shown in Figure 3 and the MG state in this article.
Figure 4 shows kinetic refolding curves measured by the CD ellipticity at 220 nm for the initial resolvable phase of the in-gel refolding reactions of four proteins, equine β-lactoglobulin, human tear lipocalin, bovine α-lactalbumin, and hen-egg white lysozyme [16]. The proteins were initially unfolded in the silica gel immersed into 6 M GdmCl, and the refolding was initiated in the gel by a GdmCl concentration jump produced by washing the protein-embedded gel with a refolding buffer solution. From Figure 4, all four proteins exhibit an exponential decay and accumulate a transient folding intermediate. The rate constants measured at different wavelengths were coincident with each other (insets of Figure 4), and the CD spectrum of the transient intermediate for each protein was obtained from the kinetic refolding curves measured at different wavelengths. From comparisons of the transient CD spectrum thus obtained with the CD spectrum of the burst-phase intermediate in free refolding in solution and also with the CD spectrum of the equilibrium MG state, if available, for each protein, it is concluded that the folding intermediate formed at the fastest initial stage in the in-gel refolding is equivalent to the intermediate formed in the burst phase in the free refolding in solution. Figure 4 also shows that the CD ellipticity value for each protein obtained by extrapolation of the exponential decay curve to zero time does not coincide with the value in the U state, indicating the presence of a burst phase before the exponential decay. However, the U state values shown in Figure 4 are those at 6 M GdmCl, and when we assume a linear dependence of the U state values on GdmCl concentration, the zero-time ellipticity becomes much closer to the values in the U state under the refolding condition. The burst-phase signals in Figure 4 may represent the rapid collapse of the protein molecules caused by changing the solvent environment from a good solvent (6 M GdmCl) to a poor solvent (~0 M GdmCl).

4. The Relationships between the Two-State and Non-Two-State Folding Reactions

From the above results, it now appears that there are two types of globular proteins in terms of the folding behavior, i.e., two-state proteins and non-two-state proteins. The two-state folding reactions under native conditions are simply represented by Equation (1), where the backward reaction is neglected because it is negligible as compared with the forward reaction:
U → k f N
On the other hand, the non-two-state proteins often accumulate the MG state at an early stage in kinetic refolding from the U state [11]. The simplest reaction scheme of the non-two-state folding is thus given by:
U ⇌ k IU k UI I → k f N
where I is the folding intermediate that has the characteristics of the MG state, kUI is the microscopic rate constant of formation of I from U, and kIU is the microscopic rate constant of formation of U from I. For the non-two-state proteins that accumulate a sufficiently stable I state during folding, kUI is much larger than kIU, and hence, kUI is eventually equivalent to the apparently observed rate constant of formation of I; this observed rate constant is denoted by kI in the following. Understanding the relationship between the two types of proteins may be important for fully elucidating the molecular mechanisms of protein folding. Although Equations (1) and (2) represent the simplest schemes for the two-state and the non-two-state folding with a single defined pathway, the folding reactions along multiple parallel pathways have also been observed in real proteins [2,75], but the arguments and conclusions in this article may be valid also for the multiple-pathway folding.
For two-state proteins, significant correlations between ln(kf) and various structure-based properties have been well documented in the late 1990s and the early 2000s. The structure-based properties include the relative contact order (RCO) [76,77,78], the long-range order (LRO) [79,80], the number of sequence-distant native pairs (Qd) [81,82], the chain topology parameter [83], the total contact distance [84], the cliquishness [85] and the relative logarithmic contact order [86]. All these structure-based properties are determined by the number of contacts between two residues in close proximity, with their Cα–Cα distance less than 6–8 Å, and by the separation lengths between the contacting residues along the primary sequence. These structure-based properties thus depend on the backbone structure (topology) of a protein, with a more complex backbone topology corresponding to a slower the folding rate. Originally, it was proposed that the correlations of ln(kf) with the above structure-based properties (except cliquishness) are present only for two-state proteins, and for non-two-state proteins, whose folding reactions may be complicated by the rate of escape from the folding intermediates trapped on a ragged energy landscape of folding, such simple correlations of ln(kf) with the structure-based properties may be absent [82].
However, later it became clear that non-two-state proteins also show similar correlations with the structure-based properties. The ln(kf) values of non-two-state proteins show a significant correlation with the chain length (L), i.e., the number of amino acid residues, and also with the absolute contact order (ACO), which is given by ACO = RCO × L [87,88,89]. Similar correlations of ln(kf) with ACO were observed for two-state and non-two-state proteins, but the correlation with L was worse for two-state proteins than for non-two-state proteins [87,89]; see also [90]. Among the above-mentioned structure-based properties for two-state proteins, three properties, LRO, Qd and the cliquishness, were also reported to exhibit significant correlations for non-two-state proteins [85,89,91]; LRO and Qd are closely related as LRO ≈ Qd/L. More sophisticated structure-based properties, well correlating with ln(kf) for both two-state and non-two-state proteins, have also been proposed [92,93,94,95,96,97,98,99,100]. These structure-based properties include the n-order contact distance [93], the geometric contact number, which is the number of nonlocal contacts well packed by a Voronoi criterion [94], the inter-residue interaction parameter, which considers the distances between all residue pairs in a protein [97], the average topological information, which is derived from probability densities of all residue pairs based on a self-avoiding random walk model [99], the entanglement of the native backbone structure [98], the logarithmic ACO, computed on shortcut networks of the native structure [100], the prediction by an algorithm for the prediction of folding and unfolding rates (PREFUR) using L and the structural class as only protein-specific input [96], and the effective chain length, evaluated from the number of helical residues and the number of helices [92] or from the amino acid composition of a protein [95]. More strict physical models to evaluate ln(kf) for two-state and non-two-state proteins from the free-energy profiles obtained by simple statistical thermodynamic calculations were also reported [101,102,103,104]. Garbuzynskiy et al. [105] have shown that physical theory, which describes scaling relationships between L and the activation free energy of folding, together with biological constraints outlines a “golden triangle” limiting the possible range of ln(kf) for single-domain globular proteins. The ln(kf) values of most of the globular proteins, including both two-state and non-two-state folders, are distributed within this narrow triangle.
Figure 5A shows the dependence of the logarithmic rate constant of folding on Qd for two-state and non-two-state proteins, taken from Kamagata et al. [89]. Qd in Figure 5A represents the number of sequence-distant native pairs, where the two residues in each pair are separated more than 12 residues along the primary sequence, and the Cα– Cα distance in space of the two residues is within 6 Å in the PDB structure of the protein [82]. Following the topomer search model of folding, proposed by Makarov and Plaxco [82], the logarithm of (kf/Qd) is plotted against Qd in Figure 5A, and the ln(kf/Qd) values are shown for 18 two-state proteins and 22 non-two-state proteins. Among the 22 non-two-state proteins, the kI values were also determined experimentally for 10 proteins, and hence, the values of ln(kI/Qd) are also presented for these proteins. Among the remaining 12 non-two-state proteins, seven exhibited the burst-phase intermediate, and five showed only a roll-over in the chevron plot; see [2] for comprehensible descriptions of the roll-over and the chevron plot. From Figure 5A, both the value of ln(kf/Qd) and its dependence on Qd in the case of the two-state proteins are very similar to those found in the plots of ln(kI/Qd) and ln(kf/Qd) against Qd for the non-two-state proteins. Such similarities between two-state and non-two-state proteins were also found rather generally in the relationships between ln(kf) and structure-based properties (see Figure 5B), which include ACO [87,88,89], LRO [91], the cliquishness [85], the n-order contact distance [93], the geometric contact number [94], the inter-residue interaction parameter [97], the average topological information [99], the entanglement of the native backbone structure [98], and the other structure-based properties [92,95,96,100]. As to the golden triangle for scaling of protein folding rates, proposed by Garbuzynskiy et al. [105], we also find similar distributions of ln(kf) between the two-state and the non-two-state proteins. These similarities thus clearly demonstrate that the physical mechanisms of folding are not essentially different between the two types of proteins.
If the mechanisms behind two-state and non-two-state folding are essentially identical, what accounts for the differences between the two-state and the non-two-state proteins? Figure 5A, which includes the data for both kf and kI for non-two-state proteins, may provide a reasonable answer to this question. From Figure 5A, the slope of the plot of ln(kf/Qd) against Qd is steeper (i.e., more negative) for the two-state proteins than for the non-two-state proteins. This trend can also be found in the plots of ln(kf) against the other structure-based properties, which usually represent the backbone topological complexity of a protein (Figure 5B) [93,94,97]. The steeper dependence of the plots against Qd and the other structure-based properties for the two-state proteins may be interpreted in terms of a shift in the rate-limiting step, caused by a change in the value of each structure-based property. From Figure 5A, the kf values of the two-state proteins are larger than those of the non-two-state proteins at Qd values less than 25, but rather coincide with the kI values of the non-two-state proteins. This apparently suggests that the U → I process of Equation (2) has become rate-limiting in these two-state proteins. At Qd values between 25 and 70, the two-state kf shows variations in either the non-two-state kI or kf, and at Qd values larger than 70, the two-state kf values rather coincide with the non-two-state kf values. It is thus suggested that the I → N process of Equation (2) is now rate-limiting in the two-state proteins with a Qd larger than 70.
From the above results, a two-stage hierarchical model is proposed as a general mechanism of protein folding. In this model, protein folding occurs in a hierarchical manner, reflecting the hierarchy of the native three-dimensional structure, as embodied in the case of non-two-state folding with an accumulation of the MG state as a folding intermediate. In this model, the protein folding process is divided into two stages: the first stage is formation of the MG state (U → I in Equation (2)); and the second stage is formation of the N state from the MG state (I → N). The two-state folding is thus merely a simplified version of the hierarchical folding caused either by an alteration in the rate-limiting step of folding or by destabilization of the intermediate. Therefore, there are two scenarios for interpreting the two-state folding behavior. First, when the backbone topology is sufficiently simple and organized by contacts between residues nearby in the primary sequence, the protein can easily form the specific conformation of side-chain packing, because the number of possible specific interactions is so small that the formation of the backbone structure rather uniquely determine the side-chain packing interactions. In this case, the first stage of the hierarchical folding becomes rate-limiting, leading to the two-state kinetic behavior of the protein. This first scenario is thus consistent with the initial barrier mechanism proposed by Krantz et al. [56]. The intermediate, even if it exists, is located in the native side after the free-energy barrier of folding (see Figure 6B). Such an intermediate is often called a hidden intermediate, and can be detected and characterized by the native-state H/D-exchange technique combined with 2D NMR spectroscopy [106]. In the second scenario, the backbone structure of the protein is organized by long-sequence-distant contacts, and hence there are too many ways to make the side-chain packing interactions quickly, keeping the second stage of the hierarchical folding rate-limiting, but destabilization of the intermediate leads to the two-state kinetic behavior. Although in this case, we apparently observe the two-state kinetics, a rigorous kinetic analysis may still follow Equation (2), and the apparently observed rate constant of folding (i.e., kf in Equation (1)) should be equivalent to ( k UI / ( k UI + k IU ) ) × k f in Equation (2). The existence of an unstable high-energy intermediate was expected experimentally for certain two-state proteins from the unfolding-limb or the refolding-limb curvature of the chevron plot [107], and this is consistent with the two-state behavior in the second scenario. The free-energy profiles of the non-two-state and the two-state folding reactions are shown in Figure 6. Here, panels (A), (B), and (C) represent the hierarchical three-state folding, the two-state folding in the first scenario, and the two-state folding in the second scenario, respectively. From Figure 6, we can clearly see that the non-two-state folding reaction and the two-state folding reactions in the first and the second scenarios are closely related with each other, and a subtle change in the free-energy landscape of folding may result in a change in the apparent folding behavior. In fact, certain proteins that normally fold in the non-two-state manner can fold with two-state kinetics, depending on the solution conditions, such as changes in pH or temperature [108,109]. There are also cases in which mutations lead to the switching between two-state and non-two-state kinetics [110,111]. All these observations are reasonably understood in terms of the two-stage hierarchical model, so that the model may provide a unified picture of the mechanisms of protein folding.

5. Conclusions

1. In classic experimental studies on protein folding, the detection and characterization of kinetic folding intermediates were basic approaches to elucidate the molecular mechanisms of folding. The MG state was found and characterized as such a folding intermediate. The MG state in protein folding was thus briefly described. The characterization of the folding intermediate by time-resolved SAXS techniques, which was originally started by a collaboration with Ptitsyn’s group, led us to conclude that the intermediate thus characterized is equivalent to the MG state.
2. Due to the observation of the two-state folding behavior for a number of small globular proteins, the two-state folding engendered much interest among researchers in the field of protein folding. Apparently, the folding intermediate is not required for complete folding of two-state proteins, and this raised questions as to the significance of the folding intermediates, such as the previously observed and characterized MG state. This was especially the case for the burst-phase intermediate that accumulates within the dead time of conventional rapid reaction techniques such as the stopped-flow method.
3. The use of the continuous-flow techniques and more recently the in-gel folding techniques made it possible to directly observe the previously unobservable solution burst-phase process of folding. The solution burst-phase intermediate was shown to be a real folding intermediate separated from the U state by a free-energy barrier. Therefore, it is now clear that there are two types of globular proteins in terms of the folding behavior, i.e., two-state proteins and non-two-state proteins.
4. There are significant correlations between the logarithmic rate constant of folding, ln(kf), and various structure-based properties such as RCO, LRO, Qd and so on. Although it was originally proposed that the correlation of ln(kf) with the structure-based properties exists only for two-state proteins, later it became clear that non-two-state proteins also show similar correlations with the structure-based properties, and hence the physical mechanisms behind two-state and non-two-state folding are essentially identical.
5. Finally, the two-stage hierarchical model was proposed as a general mechanism of protein folding. In this model, protein folding occurs in a hierarchical manner, reflecting the hierarchy of the native three-dimensional structure, as embodied in the case of non-two-state folding with an accumulation of the MG state as a folding intermediate. The two-state folding is thus merely a simplified version of the hierarchical folding caused either by an alteration in the rate-limiting step of folding or by destabilization of the intermediate. The two-stage hierarchical model also suggests that the non-two-state folding and the two-state folding are closely related with each other, and a subtle change in the free-energy landscape of folding may result in a change in the apparent folding behavior.

Funding

This work was supported by JSPS (Japan Society for the Promotion of Science) KAKENHI Grant Number JP16K07314.

Acknowledgments

The author is grateful to Valentina E. Bychkova, Heinrich Roder, Munehito Arai and Kiyoto Kamagata for their careful reading of this article and valuable comments.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Epstein, C.J.; Goldberger, R.F.; Anfinsen, C.B. The Genetic Control of Tertiary Protein Structure: Studies with Model Systems. Cold Spring Harb. Symp. Quant. Biol. 1963, 28, 439–449. [Google Scholar] [CrossRef]
  2. Morris, E.; Searle, M. Overview of Protein Folding Mechanisms: Experimental and Theoretical Approaches to Probing Energy Landscapes. Curr. Protoc. Protein Sci. 2012, 68, 28.2.1–28.2.22. [Google Scholar] [CrossRef]
  3. Sinha, K.K.; Udgaonkar, J.B. Early events in protein folding. Curr. Sci. 2009, 96, 1053–1070. [Google Scholar]
  4. Bergasa-Caceres, F.; Haas, E.; Rabitz, H. Nature’s Shortcut to Protein Folding. J. Phys. Chem. B 2019, 123, 4463–4476. [Google Scholar] [CrossRef]
  5. Barrick, D. What have we learned from the studies of two-state folders, and what are the unanswered questions about two-state protein folding? Phys. Biol. 2009, 6, 015001. [Google Scholar] [CrossRef] [Green Version]
  6. Tsytlonok, M.; Itzhaki, L.S. The how’s and why’s of protein folding intermediates. Arch. Biochem. Biophys. 2013, 531, 14–23. [Google Scholar] [CrossRef]
  7. Levinthal, C. Are there pathways for protein folding? J. Chim. Phys. 1968, 65, 44–45. [Google Scholar] [CrossRef]
  8. Wetlaufer, D.B. Nucleation, rapid folding, and globular intrachain regions in proteins. Proc. Natl. Acad. Sci. USA 1973, 70, 697–701. [Google Scholar] [CrossRef] [Green Version]
  9. Ptitsyn, O.; Pain, R.; Semisotnov, G.; Zerovnik, E.; Razgulyaev, O. Evidence for a molten globule state as a general intermediate in protein folding. FEBS Lett. 1990, 262, 20–24. [Google Scholar] [CrossRef] [Green Version]
  10. Ptitsyn, O. Molten Globule and Protein Folding. Adv. Protein Chem. 1995, 47, 83–229. [Google Scholar]
  11. Arai, M.; Kuwajima, K. Role of the molten globule state in protein folding. Adv. Protein Chem. 2000, 53, 209–282. [Google Scholar]
  12. Jackson, S.E.; Fersht, A.R. Folding of chymotrypsin inhibitor 2. 1. Evidence for a two-state transition. Biochemistry 1991, 30, 10428–10435. [Google Scholar] [CrossRef] [PubMed]
  13. Jackson, S.E. How do small single-domain proteins fold? Fold. Des. 1998, 3, R81–R91. [Google Scholar] [CrossRef] [Green Version]
  14. Onuchic, J.N.; Nymeyer, H.; García, A.E.; Chahine, J.; Socci, N.D. The energy landscape theory of protein folding: Insights into folding mechanisms and scenarios. Adv. Protein Chem. 2000, 53, 87–152. [Google Scholar]
  15. Roder, H.; Maki, K.; Cheng, H. Early Events in Protein Folding Explored by Rapid Mixing Methods. Chem. Rev. 2006, 106, 1836–1861. [Google Scholar] [CrossRef]
  16. Okabe, T.; Tsukamoto, S.; Fujiwara, K.; Shibayama, N.; Ikeguchi, M. Delineation of Solution Burst-Phase Protein Folding Events by Encapsulating the Proteins in Silica Gels. Biochemistry 2014, 53, 3858–3866. [Google Scholar] [CrossRef]
  17. Bychkova, V.E.; Semisotnov, G.V.; Balobanov, V.A.; Finkelstein, A.V. The Molten Globule Concept: 45 Years Later. Biochemistry 2018, 83, S33–S47. [Google Scholar] [CrossRef]
  18. Semisotnov, G.V.; Rodionova, N.A.; Razgulyaev, O.I.; Uversky, V.N.; Gripas’, A.F.; Gilmanshin, R.I. Study of the “molten globule” intermediate state in protein folding by a hydrophobic fluorescent probe. Biopolymers 1991, 31, 119–128. [Google Scholar] [CrossRef]
  19. Kamatari, Y. High-pressure NMR spectroscopy for characterizing folding intermediates and denatured states of proteins. Methods 2004, 34, 133–143. [Google Scholar] [CrossRef]
  20. Kuwajima, K. The molten globule state as a clue for understanding the folding and cooperativity of globular-protein structure. Proteins: Struct. Funct. Bioinform. 1989, 6, 87–103. [Google Scholar] [CrossRef]
  21. Kuwajima, K.; Hiraoka, Y.; Ikeguchi, M.; Sugai, S. Comparison of the transient folding intermediates in lysozyme and alpha-lactalbumin. Biochemisty 1985, 24, 874–881. [Google Scholar] [CrossRef]
  22. E Radford, S.; Dobson, C.M.; Evans, P.A. The folding of hen lysozyme involves partially structured intermediates and multiple pathways. Nature 1992, 358, 302–307. [Google Scholar] [CrossRef]
  23. Roder, H.; Elöve, G.A.; Englander, S.W. Structural characterization of folding intermediates in cytochrome c by H-exchange labelling and proton NMR. Nature 1988, 335, 700–704. [Google Scholar] [CrossRef] [Green Version]
  24. Udgaonkar, J.B.; Baldwin, R.L. NMR evidence for an early framework intermediate on the folding pathway of ribonuclease A. Nature 1988, 335, 694–699. [Google Scholar] [CrossRef]
  25. Jennings, P.; Wright, P.E. Formation of a molten globule intermediate early in the kinetic folding pathway of apomyoglobin. Science 1993, 262, 892–896. [Google Scholar] [CrossRef]
  26. Ikeguchi, M.; Kuwajima, K.; Mitani, M.; Sugai, S. Evidence for identity between the equilibrium unfolding intermediate and a transient folding intermediate: A comparative-study of the folding reactions of α-lactalbumin and lysozyme. Biochemistry 1996, 25, 6965–6972. [Google Scholar] [CrossRef]
  27. Arai, M.; Kuwajima, K. Rapid formation of a molten globule intermediate in refolding of α-lactalbumin. Fold. Des. 1996, 1, 275–287. [Google Scholar] [CrossRef] [Green Version]
  28. Semisotnov, G.V.; Kihara, H.; Kotova, N.V.; Kimura, K.; Amemiya, Y.; Wakabayashi, K.; Serdyuk, I.N.; Timchenko, A.A.; Chiba, K.; Nikaido, K.; et al. Protein Globularization During Folding. A Study by Synchrotron Small-angle X-ray Scattering. J. Mol. Biol. 1996, 262, 559–574. [Google Scholar] [CrossRef]
  29. Glatter, O.; Kratky, O. Small Angle X-ray Scattering; Academic Press: London, UK, 1982. [Google Scholar]
  30. Arai, M.; Ikura, T.; Semisotnov, G.V.; Kihara, H.; Amemiya, Y.; Kuwajima, K. Kinetic refolding of β-lactoglobulin. Studies by synchrotron X-ray scattering, and circular dichroism, absorption and fluorescence spectroscopy 1 1Edited by P. E. Wright. J. Mol. Biol. 1998, 275, 149–162. [Google Scholar] [CrossRef]
  31. Arai, M.; Ito, K.; Inobe, T.; Nakao, M.; Maki, K.; Kamagata, K.; Kihara, H.; Amemiya, Y.; Kuwajima, K. Fast Compaction of α-Lactalbumin During Folding Studied by Stopped-flow X-ray Scattering. J. Mol. Biol. 2002, 321, 121–132. [Google Scholar] [CrossRef]
  32. Kuwajima, K.; Yamaya, H.; Sugai, S. The Burst-phase Intermediate in the Refolding of β-Lactoglobulin Studied by Stopped-flow Circular Dichroism and Absorption Spectroscopy. J. Mol. Biol. 1996, 264, 806–822. [Google Scholar] [CrossRef]
  33. Plaxco, K.; Millett, I.S.; Segel, D.J.; Doniach, S.; Baker, D. Chain collapse can occur concomitantly with the rate-limiting step in protein folding. Nat. Struct. Biol. 1999, 6, 554–556. [Google Scholar]
  34. Amemiya, Y.; Ito, K.; Yagi, N.; Asano, Y.; Wakabayashi, K.; Ueki, T.; Endo, T. Large-aperture TV detector with a beryllium-windowed image intensifier for x-ray diffraction. Rev. Sci. Instruments 1995, 66, 2290–2294. [Google Scholar] [CrossRef]
  35. Ito, K.; Kamikubo, H.; Yagi, N.; Amemiya, Y. Correction Method and Software for Image Distortion and Nonuniform Response in Charge-Coupled Device-Based X-ray Detectors Utilizing X-ray Image intensifier. Jpn. J. Appl. Phys. 2005, 44, 8684–8691. [Google Scholar] [CrossRef]
  36. Akiyama, S.; Takahashi, S.; Kimura, T.; Ishimori, K.; Morishima, I.; Nishikawa, Y.; Fujisawa, T. Conformational landscape of cytochrome c folding studied by microsecond-resolved small-angle x-ray scattering. Proc. Natl. Acad. Sci. USA 2002, 99, 1329–1334. [Google Scholar] [CrossRef] [Green Version]
  37. Uzawa, T.; Akiyama, S.; Kimura, T.; Takahashi, S.; Ishimori, K.; Morishima, I.; Fujisawa, T. Collapse and search dynamics of apomyoglobin folding revealed by submillisecond observations of α-helical content and compactness. Proc. Natl. Acad. Sci. USA 2004, 101, 1171–1176. [Google Scholar] [CrossRef] [Green Version]
  38. Larios, E.; Li, J.; Schulten, K.; Kihara, H.; Gruebele, M. Multiple Probes Reveal a Native-like Intermediate During Low-temperature Refolding of Ubiquitin. J. Mol. Biol. 2004, 340, 115–125. [Google Scholar] [CrossRef]
  39. Kimura, T.; Akiyama, S.; Uzawa, T.; Ishimori, K.; Morishima, I.; Fujisawa, T.; Takahashi, S. Specifically Collapsed Intermediate in the Early Stage of the Folding of Ribonuclease A. J. Mol. Biol. 2005, 350, 349–362. [Google Scholar] [CrossRef]
  40. Kimura, T.; Uzawa, T.; Ishimori, K.; Morishima, I.; Takahashi, S.; Konno, T.; Akiyama, S.; Fujisawa, T. Specific collapse followed by slow hydrogen-bond formation of β-sheet in the folding of single-chain monellin. Proc. Natl. Acad. Sci. USA 2005, 102, 2748–2753. [Google Scholar] [CrossRef] [Green Version]
  41. Arai, M.; Kondrashkina, E.; Kayatekin, C.; Matthews, C.R.; Iwakura, M.; Bilsel, O. Microsecond Hydrophobic Collapse in the Folding of Escherichia coli Dihydrofolate Reductase, an α/β-Type Protein. J. Mol. Biol. 2007, 368, 219–229. [Google Scholar] [CrossRef]
  42. Konuma, T.; Kimura, T.; Matsumoto, S.; Goto, Y.; Fujisawa, T.; Fersht, A.R.; Takahashi, S. Time-Resolved Small-Angle X-ray Scattering Study of the Folding Dynamics of Barnase. J. Mol. Biol. 2011, 405, 1284–1294. [Google Scholar] [CrossRef] [PubMed]
  43. Brönnimann, C.; Eikenberry, E.F.; Henrich, B.; Horisberger, R.; Huelsen, G.; Pohl, E.; Schmitt, B.; Schulze-Briese, C.; Suzuki, M.; Tomizaki, T.; et al. The PILATUS 1M detector. J. Synchrotron Radiat. 2006, 13, 120–130. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Blanchet, C.; Spilotros, A.; Schwemmer, F.; Graewert, M.A.; Kikhney, A.; Jeffries, C.M.; Franke, D.; Mark, D.; Zengerle, R.; Cipriani, F.; et al. Versatile sample environments and automation for biological solution X-ray scattering experiments at the P12 beamline (PETRA III, DESY). J. Appl. Crystallogr. 2015, 48, 431–443. [Google Scholar] [CrossRef] [PubMed]
  45. Osaka, K.; Matsumoto, T.; Taniguchi, Y.; Inoue, D.; Sato, M.; Sano, N. High-throughput and automated SAXS/USAXS experiment for industrial use at BL19B2 in SPring-8. In AIP Conference Proceedings; AIP Publishing: New York, NY, USA, 2016; Volume 1741, p. 30003. [Google Scholar]
  46. Brandts, J.F.; Halvorson, H.R.; Brennan, M. Consideration of the possibility that the slow step in protein denaturation reactions is due to cis-trans isomerism of proline residues. Biochemistry 1975, 14, 4953–4963. [Google Scholar] [CrossRef]
  47. Yoo, T.Y.; Meisburger, S.; Hinshaw, J.; Pollack, L.; Haran, G.; Sosnick, T.R.; Plaxco, K. Small-Angle X-ray Scattering and Single-Molecule FRET Spectroscopy Produce Highly Divergent Views of the Low-Denaturant Unfolded State. J. Mol. Biol. 2012, 418, 226–236. [Google Scholar] [CrossRef] [Green Version]
  48. Manavalan, B.; Kuwajima, K.; Lee, J. PFDB: A standardized protein folding database with temperature correction. Sci. Rep. 2019, 9, 1588. [Google Scholar] [CrossRef]
  49. Silow, M.; Oliveberg, M. Transient aggregates in protein folding are easily mistaken for folding intermediates. Proc. Natl. Acad. Sci. USA 1997, 94, 6084–6086. [Google Scholar] [CrossRef] [Green Version]
  50. Went, H.M.; Benítez-Cardoza, C.G.; Jackson, S.E. Is an intermediate state populated on the folding pathway of ubiquitin? FEBS Lett. 2004, 567, 333–338. [Google Scholar] [CrossRef] [Green Version]
  51. Colón, W.; Elöve, G.A.; Wakem, L.P.; Sherman, F.; Roder, H. Side Chain Packing of the N- and C-Terminal Helices Plays a Critical Role in the Kinetics of Cytochromec Folding. Biochem. 1996, 35, 5538–5549. [Google Scholar] [CrossRef]
  52. Houry, W.A.; Rothwarf, D.M.; Scheraga, H.A. Circular Dichroism Evidence for the Presence of Burst-Phase Intermediates on the Conformational Folding Pathway of Ribonuclease A. Biochemistry 1996, 35, 10125–10133. [Google Scholar] [CrossRef]
  53. Sosnick, T.R.; Mayne, L.; Englander, S.W. Molecular collapse: The rate-limiting step in two-state cytochrome c folding. Proteins: Struct. Funct. Bioinform. 1996, 24, 413–426. [Google Scholar] [CrossRef]
  54. Sosnick, T.R.; Shtilerman, M.D.; Mayne, L.; Englander, S.W. Ultrafast signals in protein folding and the polypeptide contracted state. Proc. Natl. Acad. Sci. USA 1997, 94, 8545–8550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Qi, P.X.; Sosnick, T.R.; Englander, S.W. The burst phase in ribonuclease A folding and solvent dependence of the unfolded state. Nat. Genet. 1998, 5, 882–884. [Google Scholar] [CrossRef] [PubMed]
  56. Krantz, B.A.; Mayne, L.; Rumbley, J.; Englander, S.W.; Sosnick, T.R. Fast and Slow Intermediate Accumulation and the Initial Barrier Mechanism in Protein Folding. J. Mol. Biol. 2002, 324, 359–371. [Google Scholar] [CrossRef]
  57. Tanford, C. Protein denaturation. C. Theoretical models for the mechanism of denaturation. Adv. Protein Chem. 1970, 24, 1–95. [Google Scholar]
  58. Shastry, M.; Roder, H. Evidence for barrier-limited protein folding kinetics on the microsecond time scale. Nat. Struct. Biol. 1998, 5, 385–392. [Google Scholar] [CrossRef]
  59. Welker, E.; Maki, K.; Shastry, M.C.R.; Juminaga, D.; Bhat, R.; Scheraga, H.A.; Roder, H. Ultrarapid mixing experiments shed new light on the characteristics of the initial conformational ensemble during the folding of ribonuclease A. Proc. Natl. Acad. Sci. USA 2004, 101, 17681–17686. [Google Scholar] [CrossRef] [Green Version]
  60. Akiyama, S.; Takahashi, S.; Ishimori, K.; Morishima, I. Stepwise formation of α-helices during cytochrome c folding. Nat. Struct Biol. 2000, 7, 514–520. [Google Scholar]
  61. Qiu, L.; Zachariah, C.; Hagen, S.J. Fast Chain Contraction during Protein Folding: “Foldability” and Collapse Dynamics. Phys. Rev. Lett. 2003, 90, 168103. [Google Scholar] [CrossRef] [Green Version]
  62. Goldbeck, R.A.; Chen, E.; Kliger, D.S. Early Events, Kinetic Intermediates and the Mechanism of Protein Folding in Cytochrome c. Int. J. Mol. Sci. 2009, 10, 1476–1499. [Google Scholar] [CrossRef] [Green Version]
  63. Ellerby, L.; Nishida, C.; Nishida, F.; Yamanaka, S.; Dunn, B.; Valentine, J.; Zink, J. Encapsulation of proteins in transparent porous silicate glasses prepared by the sol-gel method. Science 1992, 255, 1113–1115. [Google Scholar] [CrossRef] [PubMed]
  64. Shibayama, N.; Saigo, S. Fixation of the Quaternary Structures of Human Adult Haemoglobin by Encapsulation in Transparent Porous Silica Gels. J. Mol. Biol. 1995, 251, 203–209. [Google Scholar] [CrossRef] [PubMed]
  65. Shibayama, N.; Saigo, S. Kinetics of the Allosteric Transition in Hemoglobin within Silicate Sol−Gels. J. Am. Chem. Soc. 1999, 121, 444–445. [Google Scholar] [CrossRef]
  66. Ronda, L.; Bruno, S.; Campanini, B.; Mozzarelli, A.; Abbruzzetti, S.; Viappiani, C.; Cupane, A.; Levantino, M.; Bettati, S. Immobilization of Proteins in Silica Gel: Biochemical and Biophysical Properties. Curr. Org. Chem. 2015, 19, 1653–1668. [Google Scholar] [CrossRef] [Green Version]
  67. Samuni, U.; Navati, M.S.; Juszczak, L.J.; Dantsker, D.; Yang, M.; Friedman, J.M. Unfolding and Refolding of Sol−Gel Encapsulated Carbonmonoxymyoglobin: An Orchestrated Spectroscopic Study of Intermediates and Kinetics. J. Phys. Chem. B 2000, 104, 10802–10813. [Google Scholar] [CrossRef]
  68. Peterson, E.S.; Leonard, E.F.; Foulke, J.A.; Oliff, M.C.; Salisbury, R.D.; Kim, D.Y. Folding Myoglobin within a Sol-Gel Glass: Protein Folding Constrained to a Small Volume. Biophys. J. 2008, 95, 322–332. [Google Scholar] [CrossRef] [Green Version]
  69. Shibayama, N. Slow Motion Analysis of Protein Folding Intermediates within Wet Silica Gels. Biochem. 2008, 47, 5784–5794. [Google Scholar] [CrossRef]
  70. Shibayama, N. Circular dichroism study on the early folding events of β-lactoglobulin entrapped in wet silica gels. FEBS Lett. 2008, 582, 2668–2672. [Google Scholar] [CrossRef] [Green Version]
  71. Woody, R.W. Circular Dichroism of Protein-Folding Intermediates. Methods in Enzymology 2004, 380, 242–285. [Google Scholar]
  72. Ohgushi, M.; Wada, A. ‘Molten-globule state’: A compact form of globular proteins with mobile side-chains. FEBS Lett. 1983, 164, 21–24. [Google Scholar] [CrossRef] [Green Version]
  73. Nakamura, S.; Seki, Y.; Katoh, E.; Kidokoro, S.-I. Thermodynamic and Structural Properties of the Acid Molten Globule State of Horse Cytochromec. Biochemistry 2011, 50, 3116–3126. [Google Scholar] [CrossRef] [PubMed]
  74. Colón, W.; Roder, H.; Col, W. Kinetic intermediates in the formation of the cytochrome c molten globule. Nat. Struct. Biol. 1996, 3, 1019–1025. [Google Scholar] [CrossRef] [PubMed]
  75. Kamagata, K.; Sawano, Y.; Tanokura, M.; Kuwajima, K. Multiple parallel-pathway folding of proline-free Staphylococcal nuclease. J. Mol. Biol. 2003, 332, 1143–1153. [Google Scholar] [CrossRef] [PubMed]
  76. Plaxco, K.; Simons, K.T.; Baker, D. Contact order, transition state placement and the refolding rates of single domain proteins. J. Mol. Biol. 1998, 277, 985–994. [Google Scholar] [CrossRef] [Green Version]
  77. Baker, D. A surprising simplicity to protein folding. Nature 2000, 405, 39–42. [Google Scholar] [CrossRef]
  78. Plaxco, K.; Simons, K.T.; Ruczinski, I.; Baker, D. Topology, stability, sequence, and length: defining the determinants of two-state protein folding kinetics. Biochemistry 2000, 39, 11177–11183. [Google Scholar] [CrossRef]
  79. Gromiha, M.M.; Selvaraj, S. Comparison between long-range interactions and contact order in determining the folding rate of two-state proteins: Application of long-range order to folding rate prediction. J. Mol. Biol. 2001, 310, 27–32. [Google Scholar] [CrossRef]
  80. Harihar, B.; Selvaraj, S. Refinement of the long-range order parameter in predicting folding rates of two-state proteins. Biopolymer 2009, 91, 928–935. [Google Scholar] [CrossRef]
  81. Makarov, D.E.; Keller, C.A.; Plaxco, K.; Metiu, H. How the folding rate constant of simple, single-domain proteins depends on the number of native contacts. Proc. Natl. Acad. Sci. USA 2002, 99, 3535–3539. [Google Scholar] [CrossRef] [Green Version]
  82. Makarov, D.E.; Plaxco, K. The topomer search model: A simple, quantitative theory of two-state protein folding kinetics. Protein Sci. 2003, 12, 17–26. [Google Scholar] [CrossRef] [Green Version]
  83. Nölting, B.; Schälike, W.; Hampel, P.; Grundig, F.; Gantert, S.; Sips, N.; Bandlow, W.; Qi, P.X. Structural determinants of the rate of protein folding. J. Theor. Biol. 2003, 223, 299–307. [Google Scholar] [CrossRef]
  84. Zhou, H.; Zhou, Y. Folding rate prediction using total contact distance. Biophys. J. 2002, 82, 458–463. [Google Scholar] [CrossRef] [Green Version]
  85. Micheletti, C. Prediction of folding rates and transition-state placement from native-state geometry. Proteins: Struct. Funct. Bioinform. 2003, 51, 74–84. [Google Scholar] [CrossRef] [PubMed]
  86. Dixit, P.D.; Weikl, T. A simple measure of native-state topology and chain connectivity predicts the folding rates of two-state proteins with and without crosslinks. Proteins: Struct. Funct. Bioinform. 2006, 64, 193–197. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  87. Galzitskaya, O.V.; Garbuzynskiy, S.O.; Ivankov, D.N.; Finkelstein, A. Chain length is the main determinant of the folding rate for proteins with three-state folding kinetics. Proteins: Struct. Funct. Bioinform. 2003, 51, 162–166. [Google Scholar] [CrossRef]
  88. Ivankov, D.N.; Garbuzynskiy, S.O.; Alm, E.; Plaxco, K.; Baker, D.; Finkelstein, A. Contact order revisited: Influence of protein size on the folding rate. Protein Sci. 2003, 12, 2057–2062. [Google Scholar] [CrossRef] [Green Version]
  89. Kamagata, K.; Arai, M.; Kuwajima, K. Unification of the folding mechanisms of non-two-state and two-state proteins. J. Mol. Biol. 2004, 339, 951–965. [Google Scholar] [CrossRef]
  90. Istomin, A.Y.; Jacobs, N.J.; Livesay, D.R. On the role of structural class of a protein with two-state folding kinetics in determining correlations between its size, topology, and folding rate. Protein Sci. 2007, 16, 2564–2569. [Google Scholar] [CrossRef] [Green Version]
  91. Harihar, B.; Selvaraj, S. Analysis of rate-limiting long-range contacts in the folding rate of three-state and two-state Proteins. Protein Pept. Lett. 2011, 18, 1042–1052. [Google Scholar] [CrossRef]
  92. Ivankov, D.N.; Finkelstein, A. Prediction of protein folding rates from the amino acid sequence-predicted secondary structure. Proc. Natl. Acad. Sci. USA 2004, 101, 8942–8944. [Google Scholar] [CrossRef] [Green Version]
  93. Zhang, L.; Sun, T. Folding rate prediction using n-order contact distance for proteins with two- and three-state folding kinetics. Biophys. Chem. 2005, 113, 9–16. [Google Scholar] [CrossRef] [PubMed]
  94. Ouyang, Z.; Liang, J. Predicting protein folding rates from geometric contact and amino acid sequence. Protein Sci. 2008, 17, 1256–1263. [Google Scholar] [CrossRef] [PubMed]
  95. Chang, L.; Wang, J.; Wang, W. Composition-based effective chain length for prediction of protein folding rates. Phys. Rev. E 2010, 82, 12. [Google Scholar] [CrossRef] [PubMed]
  96. De Sancho, D.; Muñoz, V. Integrated prediction of protein folding and unfolding rates from only size and structural class. Phys. Chem. Chem. Phys. 2011, 13, 17030. [Google Scholar] [CrossRef]
  97. Huang, S.; Huang, J.T. Inter-residue interaction is a determinant of protein folding kinetics. J. Theor. Biol. 2013, 317, 224–228. [Google Scholar] [CrossRef]
  98. Baiesi, M.; Orlandini, E.; Seno, F.; Trovato, A. Exploring the correlation between the folding rates of proteins and the entanglement of their native states. J. Phys. A Math. Theor. 2017, 50, 504001. [Google Scholar] [CrossRef] [Green Version]
  99. Censoni, L.; Martıínez, L. Prediction of kinetics of protein folding with non-redundant contact information. Bioinformatic 2018, 34, 4034–4038. [Google Scholar] [CrossRef]
  100. Khor, S. Folding with a protein’s native shortcut network. Proteins: Struct. Funct. Bioinform. 2018, 86, 924–934. [Google Scholar] [CrossRef] [Green Version]
  101. Muñoz, V.; Eaton, W.A. A simple model for calculating the kinetics of protein folding from three-dimensional structures. Proc. Natl. Acad. Sci. USA 1999, 96, 11311–11316. [Google Scholar] [CrossRef] [Green Version]
  102. Alm, E.; Morozov, A.; Kortemme, T.; Baker, D. Simple physical models connect theory and experiment in protein folding kinetics. J. Mol. Biol. 2002, 322, 463–476. [Google Scholar] [CrossRef]
  103. Garbuzynskiy, S.O.; Finkelstein, A.; Galzitskaya, O.V. Outlining folding nuclei in globular proteins. J. Mol. Biol. 2004, 336, 509–525. [Google Scholar] [CrossRef] [PubMed]
  104. Galzitskaya, O.V.; Glyakina, A.V. Nucleation-based prediction of the protein folding rate and its correlation with the folding nucleus size. Proteins: Struct. Funct. Bioinform. 2012, 80, 2711–2727. [Google Scholar] [CrossRef] [PubMed]
  105. Garbuzynskiy, S.O.; Ivankov, D.N.; Bogatyreva, N.; Finkelstein, A. Golden triangle for folding rates of globular proteins. Proc. Natl. Acad. Sci. USA 2012, 110, 147–150. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  106. Bai, Y.; Feng, H.; Zhou, Z. Population and structure determination of hidden folding intermediates by native-state hydrogen exchange-directed protein engineering and nuclear magnetic resonance. Methods Mol. Biol. 2007, 350, 69–81. [Google Scholar]
  107. Sanchez, I.E.; Kiefhaber, T. Evidence for sequential barriers and obligatory intermediates in apparent two-state protein folding. J. Mol. Biol. 2003, 325, 367–376. [Google Scholar] [CrossRef]
  108. Dalby, P.A.; Clarke, J.; Johnson, C.M.; Fersht, A.R. Folding intermediates of wild-type and mutants of barnase. II. correlation of changes in equilibrium amide exchange kinetics with the population of the folding intermediate. J. Mol. Biol. 1998, 276, 647–656. [Google Scholar] [CrossRef]
  109. McCallister, E.L.; Alm, E.; Baker, D. Critical role of β-hairpin formation in protein G folding. Nat. Struct. Biol. 2000, 7, 669–673. [Google Scholar] [CrossRef]
  110. Spudich, G.M.; Miller, E.J.; Marqusee, S. Destabilization of the Escherichia coli RNase H Kinetic Intermediate: Switching Between a Two-state and Three-state Folding Mechanism. J. Mol. Biol. 2004, 335, 609–618. [Google Scholar] [CrossRef]
  111. Friel, C.T.; Beddard, G.S.; E Radford, S. Switching Two-state to Three-state Kinetics in the Helical Protein Im9 via the Optimisation of Stabilising Non-native Interactions by Design. J. Mol. Biol. 2004, 342, 261–273. [Google Scholar] [CrossRef]
Biomolecules 10 00407 g001 550
Figure 1. Synchrotron small-angle X-ray scattering (SAXS) patterns (A) and Kratky plots (B) for bovine carbonic anhydrase (1) in the N state (0.05 M Tris-HCl (pH 8)) and (2) in the U state (0.05 M Tris-HCl (pH 8), 8.5 M urea). Reproduced with permission from [28].
Figure 1. Synchrotron small-angle X-ray scattering (SAXS) patterns (A) and Kratky plots (B) for bovine carbonic anhydrase (1) in the N state (0.05 M Tris-HCl (pH 8)) and (2) in the U state (0.05 M Tris-HCl (pH 8), 8.5 M urea). Reproduced with permission from [28].
Biomolecules 10 00407 g001
Biomolecules 10 00407 g002 550
Figure 2. Time-dependent changes in the Rg during the refolding of α-lactalbumin (A), and the Kratky plots of the protein during the refolding (B). The refolding was induced by a concentration jump of GdmCl from 4 M to 0.77 M in the presence of 1 mM CaCl2 (pH 7.0 and 4.5 °C). (A) The time-dependent changes in Rg2 were fitted to a single-exponential function, and the fitting curve of Rg thus obtained is shown as a thick red line. The dotted line shows the Rg value after the refolding is completed, and errors are the standard errors in fitting of Guinier plots. (B) The Kratky plots averaged between 10 and 20 ms after the initiation of refolding (red thick line), between 10 and 110 ms (purple line), between 10 and 960 ms (green dotted line), and between 8 and 11 s (black thick line). Blue filled circles show the Kratky plot of the equilibrium MG state at pH 2. The scattering intensities of each plot are normalized with the respective zero-angle intensities. Reproduced with permission from [31].
Figure 2. Time-dependent changes in the Rg during the refolding of α-lactalbumin (A), and the Kratky plots of the protein during the refolding (B). The refolding was induced by a concentration jump of GdmCl from 4 M to 0.77 M in the presence of 1 mM CaCl2 (pH 7.0 and 4.5 °C). (A) The time-dependent changes in Rg2 were fitted to a single-exponential function, and the fitting curve of Rg thus obtained is shown as a thick red line. The dotted line shows the Rg value after the refolding is completed, and errors are the standard errors in fitting of Guinier plots. (B) The Kratky plots averaged between 10 and 20 ms after the initiation of refolding (red thick line), between 10 and 110 ms (purple line), between 10 and 960 ms (green dotted line), and between 8 and 11 s (black thick line). Blue filled circles show the Kratky plot of the equilibrium MG state at pH 2. The scattering intensities of each plot are normalized with the respective zero-angle intensities. Reproduced with permission from [31].
Biomolecules 10 00407 g002
Biomolecules 10 00407 g003 550
Figure 3. Representative far-UV CD spectra recorded during the folding of cytochrome c at 25 °C and pH 4.5 in a silica gel containing 80% (w/w) water. (A) The times indicated are those after initiation of the folding reaction. The red circles correspond to the mean residue ellipticity [θ] obtained by the extrapolation of kinetic curves to time zero. The solid black lines represent the spectra of the N and U states in gel. The open circles and triangles are the [θ] values of the spectra of the N and U states, respectively, in solution. Note that the spectra of the solution N and U states were recorded using the same folding and unfolding conditions as were used for the gel measurements. (B) A representative kinetic curve plotted using the [θ] at 220 nm for the initial resolvable cytochrome c in-gel folding phase at 25 °C. The times are those after initiation of the folding reaction. The thick red line is the fitted curve drawn using a multi-exponential decay function. The dotted red curve is the single exponential curve. The dotted black line in the upper panel is the [θ] at 220 nm for the U state. The residuals are shown in the lower panel. The inset shows the rate constants of the initial phase calculated for different wavelengths. Reproduced with permission from [16].
Figure 3. Representative far-UV CD spectra recorded during the folding of cytochrome c at 25 °C and pH 4.5 in a silica gel containing 80% (w/w) water. (A) The times indicated are those after initiation of the folding reaction. The red circles correspond to the mean residue ellipticity [θ] obtained by the extrapolation of kinetic curves to time zero. The solid black lines represent the spectra of the N and U states in gel. The open circles and triangles are the [θ] values of the spectra of the N and U states, respectively, in solution. Note that the spectra of the solution N and U states were recorded using the same folding and unfolding conditions as were used for the gel measurements. (B) A representative kinetic curve plotted using the [θ] at 220 nm for the initial resolvable cytochrome c in-gel folding phase at 25 °C. The times are those after initiation of the folding reaction. The thick red line is the fitted curve drawn using a multi-exponential decay function. The dotted red curve is the single exponential curve. The dotted black line in the upper panel is the [θ] at 220 nm for the U state. The residuals are shown in the lower panel. The inset shows the rate constants of the initial phase calculated for different wavelengths. Reproduced with permission from [16].
Biomolecules 10 00407 g003
Biomolecules 10 00407 g004 550
Figure 4. Representative plots of the [θ] at 220 nm during the U → I transitions at 25 °C for (A) equine β-lactoglobulin, (B) human tear lipocalin, (C) bovine α-lactalbumin, and (D) hen egg-white lysozyme in the gels. The times are those after initiation of the folding reaction. The thick red lines are each a fit Table 220. nm. The residuals are shown under each kinetic plot. The insets show the rate constants for the U → I phase calculated for each wavelength, with the horizontal lines being the average values. Reproduced with permission from [16]
Figure 4. Representative plots of the [θ] at 220 nm during the U → I transitions at 25 °C for (A) equine β-lactoglobulin, (B) human tear lipocalin, (C) bovine α-lactalbumin, and (D) hen egg-white lysozyme in the gels. The times are those after initiation of the folding reaction. The thick red lines are each a fit Table 220. nm. The residuals are shown under each kinetic plot. The insets show the rate constants for the U → I phase calculated for each wavelength, with the horizontal lines being the average values. Reproduced with permission from [16]
Biomolecules 10 00407 g004
Biomolecules 10 00407 g005 550
Figure 5. A comparison of the folding rates of two-state and non-two-state proteins. (A) The ln(kf/Qd) and ln(kI/Qd) values are plotted against Qd. Filled green squares (■) represent the ln(kf/Qd) values of two-state proteins. Filled red triangles (▲) represent the ln(kI/Qd) values of 10 non-two-state proteins for which the formation of the intermediate was kinetically observed. Filled circles (●) represent the ln(kf/Qd) values for 22 non-two-state proteins, i.e., the above 10 proteins plus the 12 proteins that show a burst-phase intermediate or only a roll-over in the chevron plot. The green, red, and blue continuous lines represent the best linear fit for ln(kf/Qd) for the 18 two-state proteins, ln(kI/Qd) for the 10 non-two-state proteins, and ln(kf/Qd) for the 22 non-two-state proteins, respectively. Reproduced and modified with permission from [89]. (B) Relationship between different structure-based properties and ln(kf) of two-state (open squares) and non-two-state (filled squares) proteins. (i) Relative contact order, RCO (R = −0.15); (ii) absolute contact order, ACO (R = −0.77); (iii) chain length, L (R = −0.72); and (iv) geometric contact number, nα (R = −0.83), where R is the correlation coefficient. Reproduced with permission from [94].
Figure 5. A comparison of the folding rates of two-state and non-two-state proteins. (A) The ln(kf/Qd) and ln(kI/Qd) values are plotted against Qd. Filled green squares (■) represent the ln(kf/Qd) values of two-state proteins. Filled red triangles (▲) represent the ln(kI/Qd) values of 10 non-two-state proteins for which the formation of the intermediate was kinetically observed. Filled circles (●) represent the ln(kf/Qd) values for 22 non-two-state proteins, i.e., the above 10 proteins plus the 12 proteins that show a burst-phase intermediate or only a roll-over in the chevron plot. The green, red, and blue continuous lines represent the best linear fit for ln(kf/Qd) for the 18 two-state proteins, ln(kI/Qd) for the 10 non-two-state proteins, and ln(kf/Qd) for the 22 non-two-state proteins, respectively. Reproduced and modified with permission from [89]. (B) Relationship between different structure-based properties and ln(kf) of two-state (open squares) and non-two-state (filled squares) proteins. (i) Relative contact order, RCO (R = −0.15); (ii) absolute contact order, ACO (R = −0.77); (iii) chain length, L (R = −0.72); and (iv) geometric contact number, nα (R = −0.83), where R is the correlation coefficient. Reproduced with permission from [94].
Biomolecules 10 00407 g005
Biomolecules 10 00407 g006 550
Figure 6. The free-energy profiles of the non-two-state and the two-state folding reactions. (A) Hierarchical three-state folding, (B) two-state folding in the first scenario, and (C) two-state folding in the second scenario. The horizontal broken line in each panel shows the free-energy level of U. Arrows in panel (B) and (C) indicate the places where the free-energy level has been changes as compared with the level in panel (A).
Figure 6. The free-energy profiles of the non-two-state and the two-state folding reactions. (A) Hierarchical three-state folding, (B) two-state folding in the first scenario, and (C) two-state folding in the second scenario. The horizontal broken line in each panel shows the free-energy level of U. Arrows in panel (B) and (C) indicate the places where the free-energy level has been changes as compared with the level in panel (A).
Biomolecules 10 00407 g006

Share and Cite

MDPI and ACS Style

Kuwajima, K. The Molten Globule, and Two-State vs. Non-Two-State Folding of Globular Proteins. Biomolecules 2020, 10, 407. https://doi.org/10.3390/biom10030407

AMA Style

Kuwajima K. The Molten Globule, and Two-State vs. Non-Two-State Folding of Globular Proteins. Biomolecules. 2020; 10(3):407. https://doi.org/10.3390/biom10030407

Chicago/Turabian Style

Kuwajima, Kunihiro. 2020. "The Molten Globule, and Two-State vs. Non-Two-State Folding of Globular Proteins" Biomolecules 10, no. 3: 407. https://doi.org/10.3390/biom10030407

APA Style

Kuwajima, K. (2020). The Molten Globule, and Two-State vs. Non-Two-State Folding of Globular Proteins. Biomolecules, 10(3), 407. https://doi.org/10.3390/biom10030407

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop