The dual-basin landscape in GFP folding
Edited by Alan R. Fersht, University of Cambridge, Cambridge, United Kingdom, and approved July 2, 2008
Abstract
Recent experimental studies suggest that the mature GFP has an unconventional landscape composed of an early folding event with a typical funneled landscape, followed by a very slow search and rearrangement step into the locked, active chromophore-containing structure. As we have shown previously, the substantial difference in time scales is what generates the observed hysteresis in thermodynamic folding. The interconversion between locked and the soft folding structures at intermediate denaturant concentrations is so slow that it is not observed under the typical experimental observation time. Simulations of a coarse-grained model were used to describe the fast folding event as well as identify native-like intermediates on energy landscapes enroute to the fluorescent native fold. Interestingly, these simulations reveal structural features of the slow dynamic transition to chromophore activation. Experimental evidence presented here shows that the trapped, native-like intermediate has structural heterogeneity in residues previously linked to chromophore formation. We propose that the final step of GFP folding is a “locking” mechanism leading to chromophore formation and high stability. The combination of previous experimental work and current simulation work is explained in the context of a dual-basin folding mechanism described above.
Sign up for PNAS alerts.
Get alerts for new articles, or get an alert when an article is cited.
GFP has become a common reagent in biochemistry laboratories because of its ability to fold and form a visually fluorescent chromophore through autocatalytic cyclization and dehydration/oxidation reactions (1). Recent studies on the folding of GFP have shown multiple kinetic phases (2, 3) and high barriers to folding that lead to a slow “drift” of the unfolding equilibrium transition curve (4). In fact, a “superfolder” variant of GFP (sfGFP) (5) has an apparent hysteresis in equilibrium folding that is observed in the experiment (6, 7) shown in Fig. 1A. Considering the apparent nonthermodynamic behavior observed in experiment, can we use simulation to probe the free-energy landscape and explain these results?
Fig. 1.
Proteins whose folding has been optimized by evolution have a sufficiently smooth, funneled energy landscape, with an ensemble of folding routes directed toward the native structural basin (8–10). Coarse-grained models with potentials designed to bias the protein toward the native basin, structure-based (Gō) models (11), have been able to capture the main structural features of the folding mechanism observed in experiments. Because these models are energetically unfrustrated, the agreement to experiments demonstrate that typical well folding proteins have such small amounts of residual energetic frustration that the observed structural heterogeneity observed in folding is mostly determined by geometrical properties of the native structure.
Although it is a good folder, GFP has a folding landscape that differs from most proteins because of a structural rearrangement needed to produce an active chromophore. GFP with an active chromophore exhibits an apparent hysteresis, where folding and unfolding occur on the order of minutes to hours, whereas the hysteresis zone collapses on the order of months (6). If the active chromophore is mutated out, the protein exhibits typical fully reversible thermodynamics. Here, we propose that folding and unfolding of GFP occur in a typical folding funnel, with a proposed additional basin. Interconversion between basins via a local high barrier transition associated with the chromophore occurs, leading to the hysteresis observed in experimental folding.
The current work uses a combination of experimental and simulation techniques to probe the free-energy landscape of GFP. The simulation utilizes a topology-based model, using a recent method to improve sampling on a large, high contact-order protein (12). Simulation results focus on the U-to-N transition and are consistent with previous experimental work on the folding of GFP, picking up the frustration, intermediates, and high barriers previously observed (2–4, 13–16). Previous experimental results are presented with the simulations to explain some consequences of a dual basin in folding, and how it can lead to the nonthermodynamic behavior observed in experiment. The dual basins in the funnel lead to a barrier controlled, fast-folding step and a subsequent prefactor-controlled slow search step to lock the barrel for subsequent chromophore isomerization.
Results
Structure of GFP.
The structure of the GFP family of proteins exists as an 11-stranded β-barrel, with an unusual kinked α-helix running through the center, which contains a unique backbone rearrangement to make the fluorescent chromophore. The helix makes contacts with all of the strands of the barrel and makes contacts over the entire sequence, causing extremely high contact order. The helix contacts are highlighted to guide the eye on a contact map in supporting information (SI) Fig. S1 Bottom Right. Protein folding has been linked to the formation of the chromophore, because the protein is in its folded form (17) before kinking of the helix and chromophore formation (18). Our coarse-grained model captures some structural aspects of the chromophore that are important in describing folding through a noncanonical kink in the central α-helix. Recent work on GFP-family chromophore isomerization using crystallography and molecular dynamics has revealed chromophore isomerization in fluorescent proteins (19–21) and major differences in the interior hydrogen-bond networks between chromophore isomerizations (22–24).
Kinetic Simulations and Intermediates.
Kinetic analyses (see Materials and Methods) showed two-state folding observed in 70%, on-pathway intermediates in 26%, and several misfolded (off-pathway) species (4%) in our analysis of 250 folding trajectories. A graph of the probability mass function (PMF) of observing the protein at a given number of native contacts (Q) during all trajectories versus Q is shown in Fig. 2. Local minima at particular Q values indicate unfolded, native, and intermediate species. The unfolded state corresponds to low values of Q, whereas the folded state corresponds to a value Q≈672 (total native contacts) as indicated. We consider an intermediate to be populated if the protein samples a specific Q for 50 time steps. For simplicity of presentation, we focus on trajectories that sample two-state or intermediate states during folding. Trajectories that misfold are discussed in the next section.
Fig. 2.
We observe two multiple routes in folding. Folding always progresses toward the native state in the routes containing intermediates. Structural information on the intermediates observed during each route is given in Fig. 3. In the first route, there are two major intermediates, IK500 and IK550, which are seen in 15.0% and 13.2% of trajectories, respectively, as well as a minor intermediate, IK450 in 5.2% of trajectories. IK450, IK500, and IK550 are unfolded at their N termini, missing three, two, and one β-strand, respectively, similar to results observed in single-molecule pulling experiments (13, 14) and simulation (16). In addition, IK450 is also destabilized in the contacts between the helix and the entire barrel. Intermediate IK600 is observed on the minor route and is missing a C-terminal strand from the barrel, similar to observations in self-organized polymer simulations (25).
Fig. 3.
Experimental kinetic folding results have multiple intermediates and routes (2, 3, 6). Our simulation work supports previous experimental work as well as single-molecule pulling simulation with evidence of bifurcations in the free-energy landscape (16). Furthermore, these pathways are identified by unfolding from the C versus the N terminus, with unfolding from the C terminus the minor route, occurring in 28% of pulling simulation. Structural analysis of these intermediates shows structured areas of the protein surrounding islands of unformed contacts, denoted in figures by red circles (Fig. 3). These contacts map onto the structure near the cis-prolyl bond M88-P89 (see Discussion).
Misfolding Is Observed in Gō-Model Kinetic Simulation.
As opposed to the common perception, native structure-based models can also capture misfolding events. In this study, misfolded structures are observed 4% of the time in kinetic simulations and can be grouped into two categories (Fig. S2.). The first category (MF71 in Fig. S2.) has the central α-helix unable to pack into the formed barrel, and the second category (MF199 in Fig. S2.) is missing the first two β-strands, similar to IK450, except the N terminus is “knotted” between strands, creating entropic traps. Thus, energetically smooth Gō models can capture topological frustration that may lead to off-pathway misfolding, as observed in this study.
Thermodynamic Simulations and Intermediates.
To improve thermodynamic sampling and probe equilibrium intermediates for a large, high contact-order protein, we use a method that weights the potential energy to lower the free-energy barrier during simulation as described (12). After simulation, the results are reweighted, and multiple trajectories are analyzed independently to obtain the standard deviation of calculated parameters (n = 10). The free energy as a function of Q, as determined by the weighted-histogram analysis method (WHAM) is presented in Fig. 4C. Consistent with experimental studies, our results predict a large and broad thermodynamic barrier. One clear thermodynamic intermediate (I475) is observed.
Fig. 4.
Intermediate I475 has missing β-strands on the N-terminal end (Fig. 5), and is similar to the kinetic intermediate IK500. Intermediate I475, like IK500, also shows contacts around residue 89 with a low probability of being made (Fig. 5, circles). The simulation intermediate is missing contacts in strands 1–3, and contacts near the prolines are unformed, whereas the residues not making contacts sample an extended strand conformation. Even in the final transition state (Q = 0.75), one strand is missing, and the proline residues have not made all contacts. In the native basin (Q = 0.80), the last strand has folded, and the proline residues indicated are making contacts (Fig. S3). These simulation results support the presence of an intermediate in equilibrium folding, with unfolded areas in the structure that correspond to areas implicated in hysteresis. Furthermore, the ability to observe an intermediate with frustration in the prolines (see below) is remarkable, even with a Cα model that does not explicitly model the chromophore and proline isomerization.
Fig. 5.
Nonthermodynamic Behavior in Experiment.
The hysteresis observed in the folding of GFP (Fig. 1A) is evidence of a kinetic process controlling unfolding during a “thermodynamic” experiment. Although the zones of reversible folding and unfolding ([Gdn-HCl] < 1.5 M),([Gdn-HCl] > 5.0 M) fold and unfold on the order of minutes to hours, the hysteretic unfolding process unfolds on the order of months (5). A chevron analysis shows transition rates in the hysteresis zone are magnitudes slower than those predicted by the folding and unfolding arms. In other words, hysteretic processes are no longer dominated by folding but by local fluctuations involving the chromophore (6). Unfolding kinetics in Zone III (Fig. 1) as a function of isomerization time show a hidden phase, consistent with isomerization (6). The dependence of hysteresis on the presence of the chromophore is confirmed with mutagenesis, either by mutating a catalytic residue for chromophore formation (R96A) or removing structural requirements for chromophore formation (Y74M/M88Y). Removal of the chromophore shows thermodynamic results similar to the refolding transition in sfGFP (6).
A slow conversion between two states with differing stabilities would lead to this result. From a folding funnel perspective, this is evidence of two native-like basins with differing stability. In the hysteresis zone, unfolding is dominated by local fluctuations controlling a transition between native {Nnat} and native-like {Niso} states over months. Refolding is dominated by the {Niso} to unfolded {U} transition and entails a global, concerted process. We propose that the nonthermodynamic behavior observed in the thermodynamic folding of GFP can be explained with a dual-basin folding funnel, with native basin {Nnat}, and native-like basin {Niso} (see Discussion).
Previous work linked two prolines to the chromophore as a structural requirement for proper chromophore formation because mutation abolished both chromophore formation and hysteresis (6). Because we can trap the {Niso} form, we perform experiments on a hysteretically trapped species of GFP using a residue-level probe, NMR spectroscopy. Fig. 6A shows a 1H-15N HSQC for sfGFP (red) overlaid with the trapped species of sfGFP (blue) discussed previously (6). The chemical shift dispersion observed in the two spectra is nearly identical, evidence of a native-like structure (Fig. S5); however, three peaks show subtle differences, as shown in Fig. 6 B and C. The third peak, which is not assigned, lies at chemical shift 1H 7.99 ppm and 15N 125.5 ppm. The assigned, shifted residues lie near the prolines, which are linked to the chromophore. Once again, although proline isomerization is not included in a simple Gō model, our simulations seem to pick structural frustration in these areas (circles in Figs. 3 and 5). This area may be prone to geometrical frustration because of the complex topology of GFP, which leads to these two prolines folding in disfavored conformation (26).
Fig. 6.
Discussion
A Dual-Basin Landscape Explains Hysteresis.
Unlike usual proteins, GFP with an active chromophore folds in a dual-basin energy landscape. Folding occurs by a “typical” first step in a funnel-like landscape (10), followed by a slow search through isomerization states and structural fluctuations toward the locked active basin {Nnat}. The first transition toward the soft basin {Niso} is investigated by using a coarse-grained model, which is sufficient to extract the full ensemble of states involved in this process (11). This initial folding structure is flexible and fluctuates easily among several local minima. The flexibility of {Niso} is apparent from the different native-like intermediates observed during simulations. These intermediate-like states proceed through the slow transition into the locked, native basin of GFP {Nnat}.
GFP must form a β-barrel then close and lock the barrel around the central chromophore containing α-helix to achieve the proper native fold. Although the chromophore is able to isomerize unhindered in solution without the structural restraints of the surrounding barrel, isomerization is limited in the folded state unless there is a hula twist (22, 23, 27). This conformational transition is responsible for the slow search toward the second basin {Niso}. Interestingly, destabilization in nearby barrel strands is observed in backbone dynamics studies (4, 20) and may be related to chromophore isomerization, consistent with our simulation results and recent experimental results (28). Finally, although both cis and trans forms appear to be compatible within the folded barrel, the hydrogen-bonded network is altered between these states (24), leading to two folded forms of GFP with different stabilities.
Under strong refolding conditions (Fig. 1, Zone I), the bias toward the native state is sufficient to drive the full folding event. Under these highly native conditions, there is enough stabilizing energy to lock the protein into the native, fluorescent state, similar to results incorporating the C terminus during ClpXP degradation (29). Similarly, in strongly unfolding conditions (Fig. 1, Zone III), the high stability of the locked, native form is overcome, and the direct {Nnat} → {U} route is observed.
The hysteresis zone (Fig. 1, Zone II), must take into account both the stability of the different basins and the difference in time scales of the rates in the first and second steps. Hysteresis is a result of the first transition ({Nnat} ⇔ {Niso}) being much slower than the second one ({Niso} ⇔ {U}) during unfolding and folding. Under denaturing conditions in Zone II (Fig. 1), {Nnat} and {Niso} are less stable than {U}, but {Nnat} is kinetically stable (Fig. 1C, Zone II). If the protein is in the {Nnat} state, it will stay in that state until unlikely, random fluctuations crack the surface and move it toward {Niso}. Once this unlikely event occurs, unfolding is downhill toward {U}. These random fluctuations happen rarely because the structure of GFP limits conformational changes, leading to the slow drift of the unfolding transition as long as experimental observation times are shorter than the rare transitions.
The hysteresis observed at “equilibrium” is a consequence of the disparate time scales. Rare fluctuations occur between {Nnat} and {Niso} on an extremely slow time scale (slower than experiments and folding), but once {Niso} is populated, the transition to {U} is on a much faster time scale. However, during refolding, {U} will never populate {Niso} because of the change in stability. The reason hysteresis is observed in experiment, is because the 96-hour “equilibration” time lies between the folding time scale of minutes and the rare fluctuation time scale of months. This leads to the result that the refolding transition is in true equilibrium, but the unfolding transition is limited by a very slow kinetic process (See Fig. 1).
Experiment and Theory Support the Dual-Basin Energy Landscape.
Simulation results show that the {Niso} basin is flexible and composed of multiple conformations, many of which have a partially folded β-barrel around the central α-helix. Kinetic simulations show native-like intermediates, many of which have difficulty getting the last strand into the barrel (Fig. 3). Also, residues in the lid of the barrel (circles on Figs. 3 and 5) fail to form their native contacts until GFP is almost fully folded. Similarly, in thermodynamic simulations, the last step is the inclusion of the N-terminal β-strand into the fold to complete the β-barrel. In fact, the final transition state is still missing a strand; only in the final step does the barrel fully form, consistent with other simulations using a self-organized polymer model (25). These native-like intermediates observed in the simulation of GFP folding suggest that they are possible states involved en route to the active form. Jumps between these native-like intermediates combined with isomerization are responsible for the slow search in the approach to the final locking and activating transition.
During refolding, the barrel needs to cross a threshold of stability needed to fold, lock, and hold the isomerization of the chromophore and prolines in their native {Nnat} states. Refolding thus occurs in a concerted fashion under highly stabilizing conditions (Fig. 1., Zone I). Presumably, this locking during refolding entails simultaneously incorporating the final strand (simulation data) and capturing and holding the proper chromophore and proline isomerizations under highly refolding conditions (experimental data). Interestingly, the final strand appears to be able to form and hold an “improperly” isomerized chromophore under transition conditions, trapping the {Niso} form, as observed in 1D 1H NMR (6).
We observed subtle structural differences between {Niso} and {Nnat} in previous NMR experiments (6). In the present study, we expand upon these results with 2D 1H-15N HSQC analysis to pinpoint areas of structural heterogeneity between {Nnat} and {Niso} (Fig. 6). The changes we see in the 1H-15N HSQC spectrum map to residues adjacent to the prolines previously linked to chromophore formation (6). Although our model does not include proline isomerization, on a structural level, simulation and experiment may be compared (11, 30–33). The results from our simulations show late folding at the residues that correspond to these prolines and their neighbors (circles, Figs. 3 and 5 and Fig. S3).
Conclusions
We explore the hysteresis of GFP with the aid of a coarse-grained structure-based model. Our simulations suggest a mechanism that explains the hysteresis in the context of a dual-basin landscape consisting of an active, native form {Nnat} and a form with improper isomerizations {Niso}. GFP exhibits a fast folding step in a usual funnel landscape, followed by a slow search step toward the active chromophore basin. It is throughout this slow search where structural locking (isomerization) and chromophore activation takes place. During this second step, folding is no longer limited by a conventional folding mechanism but the required slow search. Because hysteresis is abolished in the absence of the formed chromophore (6), we propose that chromophore isomerization is key to the process that limits the unfolding kinetics in the {Nnat} → {Niso} transition and that hysteresis is a consequence of the barrel needing to form and lock down and properly position the chromophore.
Materials and Methods
The Model.
GFP was modeled by using only the Cα atoms as determined from the crystal structure 1EMA (34). The chromophore residue was broken into three Cα residues and was not modeled in a cyclized state for simulation, because the focus of this study was global folding and not chromophore formation. The energy function, a Gō-type model, is based only on the native structure and is considered energetically unfrustrated (35). The energy potential has been described in detail (31, 36).
Simulations.
Molecular dynamics (MD) simulations are used to obtain folding-unfolding transitions. For this study, all simulations were run using an unmodified sander_classic in the AMBER6 package (37), except during thermodynamic sampling runs.
Kinetic Runs.
Kinetics were initiated by unfolding the native structure with random initial velocities at 2.4 TF for 0.5–1.5 ns, and then cooling at 1.2 TF for 1 ns. Folding was initiated by running simulations for 100 ns at 0.83 TF with the final velocities from the 1.2-TF step. Statistics were attained by simulating 250 trajectories and analyzing them in five parts. Folding temperature TF was 418 K and determined by using enhanced sampling.
Thermodynamic Sampling.
Thermodynamic simulations used a modified potential to increase sampling by reweighting the force using the density of states, as described in detail by Gosavi et al. (12). The parameters determined for the well are Edepth = 13.8, Emid = 195, and σ = 92. Trajectories were simulated for 1.4 μs at temperatures between 0.998 TF and 1.005 TF. The density of states was reweighted from the modified potential, and the WHAM algorithm (38) was used to determine folding temperature and free energy from all simulations. Ten trajectories were analyzed by using WHAM to determine statistics of calculated parameters.
Experimental Equilibrium Folding Transition.
Equilibrium transitions were prepared as described (6). Chromophore fluorescence spectra were obtained by using the Fluoromax-2 spectrofluorimeter (Spex Edison) with an excitation wavelength of 450 nm and the emission measured at 508 nm unless otherwise noted. Equilibrium unfolding experiments were run with protein samples (0.1 mg/ml in 12.5 mM Tris·HCl at pH 6.8 and varying Gdn-HCl concentrations ranging from 0 to 7.5 M, at a temperature of 22°C. Experiments were equilibrated for 96 h. Refolding curves were run by equilibrating sfGFP in 7.6 M Gdn-HCl and subsequently refolded by dilution to varying final Gdn-HCl concentrations and equilibrated for 96 h. The fraction folded by chromophore fluorescence was monitored at 508 nm, after excitation at 450 nm (3-nm slit width).
NMR of Trapped Species.
The trapped sfGFP species was prepared as described (6). All NMR experiments were performed at 315 K on a Bruker Avance800 equipped with a triple-resonance triple-axis gradient probe. Final conditions for each sample were 30 mM phosphate buffer (pH 6.8), and 100 mM NaCl, 95% H2O/5% 2H2O. Data were processed by using NMRPipe (39). Two-dimensional 1H-15N HSQC spectra were acquired for sfGFP and trapped sfGFP. Spectra were acquired with a 17-ppm spectral width using a 2D HSQC pulse program (40). Assignments for sfGFP were determined by using previous assignments for GFPuv (41) and several 3D experiments including HNCA (42), HNCACB (43), TROSY HNCA (44, 45), TROSY HNCACB (44, 45), and CBCACONH (43).
Acknowledgments.
We thank the members of the laboratory for discussions of the work, and the Keck Institute for computational support. This work was funded by Molecular Biophysics Training Grant GM 08326 (to B.T.A.), Los Alamos National Laboratory (P.A.J.), National Institutes of Health Grants GM 54038 (to P.A.J. and J.N.O.) and DK 54441 (to P.A.J.), the Center for Theoretical Biophysics (Grant PHY-0822283 with additional support from MCB-0543906), and the Camille and Henry Dreyfus Start Up Award (to J.M.F.).
Supporting Information
Supporting Information (PDF)
Supporting Information
- Download
- 1.70 MB
Recent work by Orte et al. identifies parallel routes of unfolding for a GFP variant consistent with our work.
References
1
RY Tsien, The green fluorescent protein. Annu Rev Biochem 67, 509–544 (1998).
2
S Enoki, K Saeki, K Maki, K Kuwajima, Acid denaturation and refolding of green fluorescent protein. Biochemistry 43, 14238–14248 (2004).
3
H Fukuda, M Arai, K Kuwajima, Folding of green fluorescent protein and the cycle3 mutant. Biochemistry 39, 12025–12032 (2000).
4
JR Huang, TD Craggs, J Christodoulou, SE Jackson, Stable intermediate states and high energy barriers in the unfolding of GFP. J Mol Biol 370, 356–371 (2007).
5
JD Pedelacq, S Cabantous, T Tran, TC Terwilliger, GS Waldo, Engineering and characterization of a superfolder green fluorescent protein. Nat Biotechnol 24, 79–88 (2006).
6
BT Andrews, AR Schoenfish, M Roy, G Waldo, PA Jennings, The rough energy landscape of superfolder GFP is linked to the chromophore. J Mol Biol 373, 476–490 (2007).
7
B Wielgus-Kutrowska, M Narczyk, A Buszko, A Bzowska, PL Clark, Folding and unfolding of a non-fluorescent mutant of green fluorescent protein. J Phys Condens Matter 19, 285223 (2007).
8
JD Bryngelson, JN Onuchic, ND Socci, PG Wolynes, Funnels, pathways, and the energy landscape of protein folding: A synthesis. Proteins 21, 167–195 (1995).
9
PE Leopold, M Montal, JN Onuchic, Protein folding funnels: A kinetic approach to the sequence-structure relationship. Proc Natl Acad Sci USA 89, 8721–8725 (1992).
10
JN Onuchic, Z Luthey-Schulten, PG Wolynes, Theory of protein folding: The energy landscape perspective. Annu Rev Phys Chem 48, 545–600 (1997).
11
C Clementi, PA Jennings, JN Onuchic, How native-state topology affects the folding of dihydrofolate reductase and interleukin-1β. Proc Natl Acad Sci USA 97, 5871–5876 (2000).
12
S Gosavi, LL Chavez, PA Jennings, JN Onuchic, Topological frustration and the folding of interleukin-1β. J Mol Biol 357, 986–996 (2006).
13
H Dietz, F Berkemeier, M Bertz, M Rief, Anisotropic deformation response of single protein molecules. Proc Natl Acad Sci USA 103, 12724–12728 (2006).
14
H Dietz, M Rief, Exploring the energy landscape of GFP by single-molecule mechanical experiments. Proc Natl Acad Sci USA 101, 16192–16197 (2004).
15
S Enoki, et al., The equilibrium unfolding intermediate observed at pH 4 and its relationship with the kinetic folding intermediates in green fluorescent protein. J Mol Biol 361, 969–982 (2006).
16
M Mickler, et al., Revealing the bifurcation in the unfolding pathways of GFP by using single-molecule experiments and simulations. Proc Natl Acad Sci USA 104, 20268–20273 (2007).
17
BG Reid, GC Flynn, Chromophore formation in green fluorescent protein. Biochemistry 36, 6786–6791 (1997).
18
DP Barondeau, CD Putnam, CJ Kassmann, JA Tainer, ED Getzoff, Mechanism and energetics of green fluorescent protein chromophore synthesis revealed by trapped intermediate structures. Proc Natl Acad Sci USA 100, 12111–12116 (2003).
19
S Habuchi, et al., Evidence for the isomerization and decarboxylation in the photoconversion of the red fluorescent protein DsRed. J Am Chem Soc 127, 8977–8984 (2005).
20
MH Seifert, et al., Slow exchange in the chromophore of a green fluorescent protein variant. J Am Chem Soc 124, 7932–7942 (2002).
21
W Weber, V Helms, JA McCammon, PW Langhoff, Shedding light on the dark and weakly fluorescent states of green fluorescent proteins. Proc Natl Acad Sci USA 96, 6177–6182 (1999).
22
MC Chen, CR Lambert, JD Urgitis, M Zimmer, Photoisomerization of green fluorescent protein and the dimensions of the chromophore cavity. Chem Phys 270, 157–164 (2001).
23
SL Maddalo, M Zimmer, The role of the protein matrix in green fluorescent protein fluorescence. Photochem Photobiol 82, 367–372 (2006).
24
R Nifosi, V Tozzini, Cis-trans photolsomerization of the chromophore in the green fluorescent protein variant E(2)GFP: A molecular dynamics study. Chem Phys 323, 358–368 (2006).
25
C Hyeon, RI Dima, D Thirumalai, Pathways and kinetic barriers in mechanical unfolding and refolding of RNA and proteins. Structure (London) 14, 1633–1645 (2006).
26
U Reimer, et al., Side-chain effects on peptidyl-prolyl cis/trans isomerisation. J Mol Biol 279, 449–460 (1998).
27
SS Patnaik, S Trohalaki, R Pachter, Molecular modeling of green fluorescent protein: Structural effects of chromophore deprotonation. Biopolymers 75, 441–452 (2004).
28
H Mizuno, et al., Light-dependent regulation of structural flexibility in a photochromic fluorescent protein. Proc Natl Acad Sci USA 105, 9927–9932 (2008).
29
A Martin, TA Baker, RT Sauer, Protein unfolding by a AAA+ protease is dependent on ATP-hydrolysis rates and substrate energy landscapes. Nat Struct Mol Biol 15, 139–145 (2008).
30
RD Schaeffer, A Fersht, V Daggett, Combining experiment and simulation in protein folding: closing the gap for small model systems. Curr Opin Struct Biol 18, 4–9 (2008).
31
JM Finke, JN Onuchic, Equilibrium and kinetic folding pathways of a TIM barrel with a funneled energy landscape. Biophys J 89, 488–505 (2005).
32
J Kubelka, WA Eaton, J Hofrichter, Experimental tests of villin subdomain folding simulations. J Mol Biol 329, 625–630 (2003).
33
Z Gu, MK Rao, WR Forsyth, JM Finke, CR Matthews, Structural analysis of kinetic folding intermediates for a TIM barrel protein, indole-3-glycerol phosphate synthase, by hydrogen exchange mass spectrometry and Go model simulation. J Mol Biol 374, 528–546 (2007).
34
M Ormo, et al., Crystal structure of the Aequorea victoria green fluorescent protein. Science 273, 1392–1395 (1996).
35
N Go, Theoretical studies of protein folding. Annu Rev Biophys Bioeng 12, 183–210 (1983).
36
C Clementi, H Nymeyer, JN Onuchic, Topological and energetic factors: What determines the structural details of the transition state ensemble and “en-route” intermediates for protein folding? An investigation for small globular proteins. J Mol Biol 298, 937–953 (2000).
37
DA Pearlman, et al., Amber, a package of computer-programs for applying molecular mechanics, normal-mode analysis, molecular-dynamics and free-energy calculations to simulate the structural and energetic properties of molecules. Comput Phys Commun 91, 1–41 (1995).
38
S Kumar, D Bouzida, RH Swendsen, PA Kollman, JM Rosenberg, The weighted histogram analysis method for free-energy calculations on biomolecules. 1. The method. J Comput Chem 13, 1011–1021 (1992).
39
F Delaglio, et al., NMRPipe: a multidimensional spectral processing system based on UNIX pipes. J Biomol NMR 6, 277–293 (1995).
40
V Sklenar, M Piotto, R Leppik, V Saudek, Gradient-tailored water suppression for H-1-N-15 Hsqc experiments optimized to retain full sensitivity. J Magn Reson Ser A 102, 241–245 (1993).
41
J Georgescu, T Rehm, J Wiehler, B Steipe, TA Holak, Backbone H(N), N, C(alpha) and C(beta) assignment of the GFPuv mutant. J Biomol NMR 25, 161–162 (2003).
42
LE Kay, GY Xu, T Yamazaki, enhanced-sensitivity triple-resonance spectroscopy with minimal H2O saturation. J Magn Reson Ser A 109, 129–133 (1994).
43
DR Muhandiram, LE Kay, gradient-enhanced triple-resonance 3-dimensional NMR experiments with improved sensitivity. J Magn Reson Ser B 103, 203–216 (1994).
44
A Eletsky, A Kienhofer, K Pervushin, TROSY NMR with partially deuterated proteins. J Biomol NMR 20, 177–180 (2001).
45
M Salzmann, K Pervushin, G Wider, H Senn, K Wuthrich, TROSY in triple-resonance experiments: New perspectives for sequential NMR assignment of large proteins. Proc Natl Acad Sci USA 95, 13585–13590 (1998).
46
A Orte, TD Craggs, SS White, SE Jackson, D Klenerman, Evidence of an intermediate and parallel pathways in protein unfolding from single-molecule fluorescence. J Am Chem Soc 130, 7898–7907 (2008).
Information & Authors
Information
Published in
Classifications
Copyright
© 2008 by The National Academy of Sciences of the USA.
Submission history
Received: April 28, 2008
Published online: August 26, 2008
Published in issue: August 26, 2008
Keywords
Acknowledgments
We thank the members of the laboratory for discussions of the work, and the Keck Institute for computational support. This work was funded by Molecular Biophysics Training Grant GM 08326 (to B.T.A.), Los Alamos National Laboratory (P.A.J.), National Institutes of Health Grants GM 54038 (to P.A.J. and J.N.O.) and DK 54441 (to P.A.J.), the Center for Theoretical Biophysics (Grant PHY-0822283 with additional support from MCB-0543906), and the Camille and Henry Dreyfus Start Up Award (to J.M.F.).
Notes
This article is a PNAS Direct Submission.
This article contains supporting information online at www.pnas.org/cgi/content/full/0804039105/DCSupplemental.
Authors
Competing Interests
The authors declare no conflict of interest.
Metrics & Citations
Metrics
Citation statements
Altmetrics
Citations
If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.
Cited by
Loading...
View Options
View options
PDF format
Download this article as a PDF file
DOWNLOAD PDFGet Access
Login options
Check if you have access through your login credentials or your institution to get full access on this article.
Personal login Institutional LoginRecommend to a librarian
Recommend PNAS to a LibrarianPurchase options
Purchase this article to get full access to it.