Two-state Manner Research Articles

Barrierless Transition Identified during Folding of Barstar by using Time-Resolved FRET from 5-Flurotryptophan

Despite a large volume of theoretical support for barrierless (downhill) folding of proteins, the extent of experimental observation is rather limited. This limitation is even more striking when we address the kinetic process of folding rather than equilibrium titrations of folding-unfolding transition. In this work we use the fluorescence of a single 5-Fluorotryptophan (F-Trp) in barstar as the energy donor and study the evolution of intramolecular distances and distance distributions during the kinetics of folding. Enhanced homogeneity of fluorescence decay kinetics of F-Trp when compared to that of Trp enables unequivocal interpretation of distribution of fluorescence lifetimes obtained from the Maximum Entropy Method (MEM) as arising from folded and unfolded conformations of the protein. Of the two intramolecular distances monitored during the folding process, W53-C82 distance showed continuous decrease while W53-C40 distance behaved in a two-state manner. Implication of these observations to the folding process is discussed.

The complex folding behavior of HIV-1-protease monomer revealed by optical-tweezer single-molecule experiments and molecular dynamics simulations

We have used optical tweezers and molecular dynamics simulations to investigate the unfolding and refolding process of a stable monomeric form of HIV-1-protease (PR). We have characterized the behavior under tension of the native state (N), and that of the ensemble of partially folded (PF) conformations the protein visits en route to N, which collectively act as a long-lived state controlling the slow kinetic phase of the folding process. Our results reveal a rich network of unfolding events, where the native state unfolds either in a two-state manner or by populating an intermediate state I, while the PF state unravels through a multitude of pathways, underscoring its structural heterogeneity. Refolding of mechanically denatured HIV-1-PR monomers is also a multiple-pathway process. Molecular dynamics simulations allowed us to gain insight into possible conformations the protein adopts along the unfolding pathways, and provide information regarding possible structural features of the PF state.

Ca2+ binding enhanced mechanical stability of an archaeal crystallin.

Structural topology plays an important role in protein mechanical stability. Proteins with β-sandwich topology consisting of Greek key structural motifs, for example, I27 of muscle titin and 10FNIII of fibronectin, are mechanically resistant as shown by single-molecule force spectroscopy (SMFS). In proteins with β-sandwich topology, if the terminal strands are directly connected by backbone H-bonding then this geometry can serve as a “mechanical clamp”. Proteins with this geometry are shown to have very high unfolding forces. Here, we set out to explore the mechanical properties of a protein, M-crystallin, which belongs to β-sandwich topology consisting of Greek key motifs but its overall structure lacks the “mechanical clamp” geometry at the termini. M-crystallin is a Ca2+ binding protein from Methanosarcina acetivorans that is evolutionarily related to the vertebrate eye lens β and γ-crystallins. We constructed an octamer of crystallin, (M-crystallin)8, and using SMFS, we show that M-crystallin unfolds in a two-state manner with an unfolding force ∼90 pN (at a pulling speed of 1000 nm/sec), which is much lower than that of I27. Our study highlights that the β-sandwich topology proteins with a different strand-connectivity than that of I27 and 10FNIII, as well as lacking “mechanical clamp” geometry, can be mechanically resistant. Furthermore, Ca2+ binding not only stabilizes M-crystallin by 11.4 kcal/mol but also increases its unfolding force by ∼35 pN at the same pulling speed. The differences in the mechanical properties of apo and holo M-crystallins are further characterized using pulling speed dependent measurements and they show that Ca2+ binding reduces the unfolding potential width from 0.55 nm to 0.38 nm. These results are explained using a simple two-state unfolding energy landscape.

Probing the origins of two-state folding

Many protein systems fold in a two-state manner. Random models, however, rarely display two-state kinetics and thus such behavior should not be accepted as a default. While theories for the prevalence of two-state kinetics have been presented, none sufficiently explain the breadth of experimental observations. A model, making minimal assumptions, is introduced that suggests two-state behavior is likely for any system with an overwhelmingly populated native state. We show two-state folding is a natural consequence of such two-state thermodynamics, and is strengthened by increasing the population of the native state. Further, the model exhibits hub-like behavior, with slow interconversions between unfolded states. Despite this, the unfolded state equilibrates quickly relative to the folding time. This apparent paradox is readily understood through this model. Finally, our results compare favorable with measurements of folding rates as a function of chain length and Keq, providing new insight into these relations.

High Thermal Stability and Unique Trimer Formation of Cytochromec′ from ThermophilicHydrogenophilus thermoluteolus

Sequence analysis indicated that thermophilic Hydrogenophilus thermoluteolus cytochrome c' (PHCP) and its mesophilic homolog, Allochromatium vinosum cytochrome c' (AVCP), closely resemble each other in a phylogenetic tree of the cytochrome c' family, with 55% sequence identity. The denaturation temperature of PHCP was 87 °C, 35 °C higher than that of AVCP. Furthermore, PHCP exhibited a larger enthalpy change value during its thermal denaturation than AVCP. While AVCP was dimeric, as observed previously, PHCP was trimeric, and this was the first observation as a cytochrome c'. Dissociation of trimeric PHCP and its protein denaturation reversibly occurred at the same time in a two-state transition manner. Therefore, PHCP is enthalpically more stable than AVCP, perhaps due to its unique trimeric form, in addition to the lower number of Gly residues in its putative α-helical regions.

Thermal unfolding of the N-terminal region of p53 monitored by circular dichroism spectroscopy

It has been estimated that 30% of eukaryotic protein and 70% of transcription factors are intrinsically disordered (ID). The biochemical significance of proteins that lack stable tertiary structure, however, is not clearly understood, largely owing to an inability to assign well-defined structures to specific biological tasks. In an attempt to investigate the structural character of ID protein, we have measured the circular dichroism spectrum of the N-terminal region of p53 over a range of temperatures and solution conditions. p53 is a well-studied transcription factor that has a proline-rich N-terminal ID region containing two activation domains. High proline content is a property commonly associated with ID, and thus p53 may be a good model system for investigating the biochemical importance of ID. The spectra presented here suggest that the N-terminal region of p53 may adopt an ordered structure under physiological conditions and that this structure can be thermally unfolded in an apparent two-state manner. The midpoint temperature for this thermal unfolding of the N-terminal region of p53 was at the near-physiological temperature of 39°C, suggesting the possibility of a physiological role for the observed structural equilibrium.

Multiple States of Dimeric Aggregates Formed by (Amido–ethynyl)helicene Bidomain Compound and (Amido–ethynyl–amido)helicene Tridomain Compound

An (amido-ethynyl)helicene bidomain compound and an (amido-ethynyl-amido)helicene tridomain compound were synthesized. The multidomain compounds were designed on the basis of previous findings that amido and ethynyl oligomers form dimeric aggregates with properties orthogonal to each other. Four aggregate states of multidomain compounds, namely, all-dimer, amido-dimer, ethynyl-dimer, and random-coil states, were obtained in different solvents, which were analyzed by circular dichroism (CD), UV/Vis, (1)H NMR, and IR spectroscopy; vapor pressure osmometry (VPO); dynamic light scattering (DLS); and atomic force microscopy (AFM). The amido and ethynyl domains independently aggregated and disaggregated in a two-state manner. Reversible structural changes occurred for a tridomain compound between the ethynyl-dimer/random-coil state and the all-dimer/amido-dimer state with heating and cooling. Two structural change processes with different properties were obtained using a single compound.

Denaturant-dependent folding of GFP

We use molecular simulations using a coarse-grained model to map the folding landscape of Green Fluorescent Protein (GFP), which is extensively used as a marker in cell biology and biotechnology. Thermal and Guanidinium chloride (GdmCl) induced unfolding of a variant of GFP, without the chromophore, occurs in an apparent two-state manner. The calculated midpoint of the equilibrium folding in GdmCl, taken into account using the Molecular Transfer Model (MTM), is in excellent agreement with the experiments. The melting temperatures decrease linearly as the concentrations of GdmCl and urea are increased. The structural features of rarely populated equilibrium intermediates, visible only in free energy profiles projected along a few order parameters, are remarkably similar to those identified in a number of ensemble experiments in GFP with the chromophore. The excellent agreement between simulations and experiments show that the equilibrium intermediates are stabilized by the chromophore. Folding kinetics, upon temperature quench, show that GFP first collapses and populates an ensemble of compact structures. Despite the seeming simplicity of the equilibrium folding, flux to the native state flows through multiple channels and can be described by the kinetic partitioning mechanism. Detailed analysis of the folding trajectories show that both equilibrium and several kinetic intermediates, including misfolded structures, are sampled during folding. Interestingly, the intermediates characterized in the simulations coincide with those identified in single molecule pulling experiments. Our predictions, amenable to experimental tests, show that MTM is a practical way to simulate the effect of denaturants on the folding of large proteins.

Thermodynamic Characterization of RNA 2 × 3 Nucleotide Internal Loops

To better elucidate RNA structure-function relationships and to improve the design of pharmaceutical agents that target specific RNA motifs, an understanding of RNA primary, secondary, and tertiary structure is necessary. The prediction of RNA secondary structure from sequence is an intermediate step in predicting RNA three-dimensional structure. RNA secondary structure is typically predicted using a nearest neighbor model based on free energy parameters. The current free energy parameters for 2 × 3 nucleotide loops are based on a 23-member data set of 2 × 3 loops and internal loops of other sizes. A database of representative RNA secondary structures was searched to identify 2 × 3 nucleotide loops that occur in nature. Seventeen of the most frequent 2 × 3 nucleotide loops in this database were studied by optical melting experiments. Fifteen of these loops melted in a two-state manner, and the associated experimental ΔG°(37,2×3) values are, on average, 0.6 and 0.7 kcal/mol different from the values predicted for these internal loops using the predictive models proposed by Lu, Turner, and Mathews [Lu, Z. J., Turner, D. H., and Mathews, D. H. (2006) Nucleic Acids Res. 34, 4912-4924] and Chen and Turner [Chen, G., and Turner, D. H. (2006) Biochemistry 45, 4025-4043], respectively. These new ΔG°(37,2×3) values can be used to update the current algorithms that predict secondary structure from sequence. To improve free energy calculations for duplexes containing 2 × 3 nucleotide loops that still do not have experimentally determined free energy contributions, an updated predictive model was derived. This new model resulted from a linear regression analysis of the data reported here combined with 31 previously studied 2 × 3 nucleotide internal loops. Most of the values for the parameters in this new predictive model are within experimental error of those of the previous models, suggesting that approximations and assumptions associated with the derivation of the previous nearest neighbor parameters were valid. The updated predictive model predicts free energies of 2 × 3 nucleotide internal loops within 0.4 kcal/mol, on average, of the experimental free energy values. Both the experimental values and the updated predictive model can be used to improve secondary structure prediction from sequence.

An expanded view of the protein folding landscape of PDZ domains

Most protein domains fold in an apparently co-operative and two-state manner with only the native and denatured states significantly populated at any experimental condition. However, the protein folding energy landscape is often rugged and different transition states may be rate limiting for the folding reaction under different conditions, as seen for the PDZ protein domain family. We have here analyzed the folding kinetics of two PDZ domains and found that a previously undetected third transition state is rate limiting under conditions that stabilize the native state relative to the denatured state. In light of these results, we have re-analyzed previous folding data on PDZ domains and present a unified folding mechanism with three distinct transition states separated by two high-energy intermediates. Our data show that sequence composition tunes the relative stabilities of folding transition states within the PDZ family, while the overall mechanism is determined by topology. This model captures the kinetic folding mechanism of all PDZ domains studied to date.

GB1 Is Not a Two-State Folder: Identification and Characterization of an On-Pathway Intermediate

The folding pathway of the small α/β protein GB1 has been extensively studied during the past two decades using both theoretical and experimental approaches. These studies provided a consensus view that the protein folds in a two-state manner. Here, we reassessed the folding of GB1, both by experiments and simulations, and detected the presence of an on-pathway intermediate. This intermediate has eluded earlier experimental characterization and is distinct from the collapsed state previously identified using ultrarapid mixing. Failure to identify the presence of an intermediate affects some of the conclusions that have been drawn for GB1, a popular model for protein folding studies.

Unzip Single Protein Zippers using Optical Tweezers

Basic leucine zippers (bZIPs) are dimeric alpha-helical domains found in many transcription factors. They comprise coiled coil regions (zippers) that mediate specific protein dimerization, and basic regions that undergo disorder-to-order transitions to grip their cognate DNA binding sites. It remains unclear how folding of the two regions and binding of DNA are coupled. We investigated the detailed folding kinetics of a model bZIP domain derived from GCN4 transcriptional factor, using high-resolution optical tweezers. When pulled from the same amino- or carboxyl-terminals, the coiled coil alone folds and unfolds rapidly in a two-state manner, with little energy barrier. The folding kinetics of the full GCN4 bZIP domain is similar to the coiled coil alone. However, the presence of its binding site (AP1) induces a third new state corresponding to the bZIP-DNA complex. The probability of this new state critically depends upon the AP1 concentration, which reaches a maximum at an intermediate AP1 concentration. Results from our single-molecule manipulation experiments reveal a complex folding energy landscape even for a small protein domain that can be further tuned by their binding ligands.

Role of Protein Stabilizers on the Conformation of the Unfolded State of Cytochrome c and Its Early Folding Kinetics: INVESTIGATION AT SINGLE MOLECULAR RESOLUTION

An insight into the conformation and dynamics of unfolded and early intermediate states of a protein is essential to understand the mechanism of its aggregation and to design potent inhibitor molecules. Fluorescence correlation spectroscopy has been used to study the effects of several model protein stabilizers on the conformation of the unfolded state and early folding dynamics of tetramethyl rhodamine-labeled cytochrome c from Saccharomyces cerevisiae at single molecular resolution. Special attention has been given to arginine, which is a widely used stabilizer for improving refolding yield of different proteins. The value of the hydrodynamic radius (r(H)) obtained by analyzing the intensity fluctuations of the diffusing molecules has been found to increase in a two-state manner as the protein is unfolded by urea. The results further show that the presence of arginine and other protein stabilizers favors a relatively structured conformation of the unfolded states (r(H) of 29 A) over an extended one (r(H) of 40 A), which forms in their absence. Also, the time constant of a kinetic component (tau(R)) of about 30 micros has been observed by analyzing the correlation functions, which represents formation of a collapsed state. This time constant varies with urea concentration representing an inverted Chevron plot that shows a roll-over and behavior in the absence of arginine. To the best of our knowledge, this is one of the first applications of fluorescence correlation spectroscopy to study direct folding kinetics of a protein.

Engineering a two-helix bundle protein for folding studies

The SAP domain from the Saccharomyces cerevisiae THO1 protein contains a hydrophobic core and just two α-helices. It could provide a system for studying protein folding that bridges the gap between studies on isolated helices and those on larger protein domains. We have engineered the SAP domain for protein folding studies by inserting a tryptophan residue into the hydrophobic core (L31W) and solved its structure. The helical regions had a backbone root mean-squared deviation of 0.9 Å from those of wild type. The mutation L31W destabilised wild type by 0.8 ± 0.1 kcal mol−1. The mutant folded in a reversible, apparent two-state manner with a microscopic folding rate constant of around 3700 s−1 and is suitable for extended studies of folding.

Dynamic force spectroscopy of DNA hairpins: II. Irreversibility and dissipation

We investigate irreversibility and dissipation in single molecules that cooperativelyfold/unfold in a two-state manner under the action of mechanical force. We apply paththermodynamics to derive analytical expressions for the average dissipated work and theaverage hopping number in two-state systems. It is shown how these quantities only dependon two parameters that characterize the folding/unfolding kinetics of the molecule: thefragility and the coexistence hopping rate. The latter has to be rescaled to take intoaccount the appropriate experimental set-up. Finally we carry out pulling experimentswith optical tweezers in a specifically designed DNA hairpin that shows two-statecooperative folding. We then use these experimental results to validate our theoreticalpredictions.

Coupling of Domain Swapping to Kinetic Stability in a Thioredoxin Mutant

The thioredoxin (Trx) fold is a small monomeric domain that is ubiquitous in redox-active enzymes. Trxs are characterized by a typical WCGPC active-site sequence motif. A single active-site mutation of the tryptophan to an alanine in Staphylococcus aureus Trx converts the oxidized protein into a biologically inactive domain-swapped dimer. While the monomeric protein unfolds reversibly in a two-state manner, the oxidized dimeric form is kinetically stable and converts to the monomeric form upon refolding. After reduction, the half-life of the dimer decreases many orders of magnitude to ∼4.3 h, indicating that the active-site disulfide between Cys29 and Cys32 is an important determinant for the kinetics of unfolding. We propose kinetic stability as a possible evolutionary strategy in the evolution of multimeric proteins from their monomeric ancestors by domain swapping, which, for this biologically inactive Trx mutant, turned out to be an evolutionary dead end.

Structural Insight into RNA Hairpin Folding Intermediates

Hairpins are a ubiquitous secondary structure motif in RNA molecules. Despite their simple structure, there is some debate over whether they fold in a two-state or multi-state manner. We have studied the folding of a small tetraloop hairpin using a serial version of replica exchange molecular dynamics on a distributed computing environment. On the basis of these simulations, we have identified a number of intermediates that are consistent with experimental results. We also find that folding is not simply the reverse of high-temperature unfolding and suggest that this may be a general feature of biomolecular folding.

Up and down states

Operationally, up and down states refer to the observation that neurons have two preferred subthreshold membrane potentials, both subthreshold for action potential generation. The most common method of detecting these states is the so-called all-points histogram (Cowan and Wilson, 1994; Wilson and Kawaguchi, 1996; Pare et al., 1998). This method, borrowed from the study of single channel currents, plots the frequency of occurrence of various values of membrane potential for every point in a digitized record. Figure 1 shows a typical example for a striatal spiny neuron, and for a cortical pyramidal cell in layer V, recorded simultaneously. The histogram to the left of each one shows the amount of time the cell spends at each value of membrane potential. Both cells toggle between two preferred membrane potentials, one very hyperpolarized (Down state), and one more depolarized (Up state). In both cells, the Up state is only a few millivolts from the action potential threshold. Usually, membrane potential fluctuations around the Up state are of higher amplitude, whereas the Down state is relatively free of noise. Figure 1 Up and down states. Up and Down states have been most often studied in animals anesthetized with urethane or other anesthetics that induce slow coherent oscillations in the cortex similar to those seen during sleep or anesthesia (e.g. Steriade et al., 1993; Steriade et al., 2001; Mahon et al., 2006; Destexhe et al. 2007). Because of this, Up and Down states are sometimes used as a synonym for slow oscillations. That usage is avoided here, and, Up and Down states will refer only to the set of cellular and network properties that causes neurons to respond to synaptic input in a two-state manner. These cellular properties arise from ionic conductances that are always present in the cell, and can continue to influence cellular activity in other circumstances not associated with slow oscillations. For example, much of the work on Up and Down states has been obtained in tissue slices, which are neither anethetized nor asleep, and often not showing any slow oscillations. Neurons may exhibit two-state behavior because of their intrinsic properties, or because they are in a network that imposes it on them, or both, and may be expressed as a part of a variety of activity states.

Folding thermodynamics and kinetics of the leucine‐rich repeat domain of the virulence factor Internalin B

Although the folding of alpha-helical repeat proteins has been well characterized, much less is known about the folding of repeat proteins containing beta-sheets. Here we investigate the folding thermodynamics and kinetics of the leucine-rich repeat (LRR) domain of Internalin B (InlB), an extracellular virulence factor from the bacterium Lysteria monocytogenes. This domain contains seven tandem leucine-rich repeats, of which each contribute a single beta-strand that forms a continuous beta-sheet with neighboring repeats, and an N-terminal alpha-helical capping motif. Despite its modular structure, InlB folds in an equilibrium two-state manner, as reflected by the identical thermodynamic parameters obtained by monitoring its sigmoidal urea-induced unfolding transition by different spectroscopic probes. Although equilibrium two-state folding is common in alpha-helical repeat proteins, to date, InlB is the only beta-sheet-containing repeat protein for which this behavior is observed. Surprisingly, unlike other repeat proteins exhibiting equilibrium two-state folding, InlB also folds by a simple two-state kinetic mechanism lacking intermediates, aside from the effects of prolyl isomerization on the denatured state. However, like other repeat proteins, InlB also folds significantly more slowly than expected from contact order. When plotted against urea, the rate constants for the fast refolding and single unfolding phases constitute a linear chevron that, when fitted with a kinetic two-state model, yields thermodynamic parameters matching those observed for equilibrium folding. Based on these kinetic parameters, the transition state is estimated to comprise 40% of the total surface area buried upon folding, indicating that a large fraction of the native contacts are formed in the rate-limiting step to folding.

Effect of bilayer morphology on the subgel phase formation

We examined how morphology of bilayer assemblies affects the kinetics of the subgel phase formation in dimyristoylphosphatidylglycerol (DMPG) bilayers, which change their morphology depending on NaCl concentration. Quantitative analysis of the kinetics revealed that in flat sheet-like structures (bilayer sheets) the subgel phase forms in a simple two-state manner with the relaxation time of about 3 min at −10 °C while in vesicles it forms much slower under a multi-step process. Freeze-etch electron microscopic observations suggested that the kinetics of the subgel phase formation is directly correlated with the morphology of bilayer assemblies. It is likely that the bilayer sheet structure is more favorable to the subgel phase formation in DMPG bilayers than the vesicular structure.