Observation of strange kinetics in protein folding

Highly nonexponential folding kinetics in aqueous solution have been observed during temperature jump-induced refolding of two proteins, yeast phosphoglycerate kinase and a ubiquitin mutant. The observations are most easily interpreted in terms of downhill folding, which posits a heterogeneous ensemble of structures en route to the folded state. The data are also reconciled with exponential kinetics measured under different experimental conditions and with titration experiments indicating cooperative folding.

Magnesium-dependent folding of self-splicing RNA: Exploring the link between cooperativity, thermodynamics, and kinetics

Folding of the Tetrahymena self-splicing RNA into its active conformation involves a set of discrete intermediate states. The Mg2+-dependent equilibrium transition from the intermediates to the native structure is more cooperative than the formation of the intermediates from the unfolded states. We show that the degree of cooperativity is linked to the free energy of each transition and that the rate of the slow transition from the intermediates to the native state decreases exponentially with increasing Mg2+ concentration. Monovalent salts, which stabilize the folded RNA nonspecifically, induce states that fold in less than 30 s after Mg2+ is added to the RNA. A simple model is proposed that predicts the folding kinetics from the Mg2+-dependent change in the relative stabilities of the intermediate and native states.

Stretching single-domain proteins: Phase diagram and kinetics of force-induced unfolding

Single-molecule force spectroscopy reveals unfolding of domains in titin on stretching. We provide a theoretical framework for these experiments by computing the phase diagrams for force-induced unfolding of single-domain proteins using lattice models. The results show that two-state folders (at zero force) unravel cooperatively, whereas stretching of non-two-state folders occurs through intermediates. The stretching rates of individual molecules show great variations reflecting the heterogeneity of force-induced unfolding pathways. The approach to the stretched state occurs in a stepwise “quantized” manner. Unfolding dynamics and forces required to stretch proteins depend sensitively on topology. The unfolding rates increase exponentially with force f till an optimum value, which is determined by the barrier to unfolding when f = 0. A mapping of these results to proteins shows qualitative agreement with force-induced unfolding of Ig-like domains in titin. We show that single-molecule force spectroscopy can be used to map the folding free energy landscape of proteins in the absence of denaturants.

Calcium release from the endoplasmic reticulum of higher plants elicited by the NADP metabolite nicotinic acid adenine dinucleotide phosphate

Higher plants share with animals a responsiveness to the Ca2+ mobilizing agents inositol 1,4,5-trisphosphate (InsP3) and cyclic ADP-ribose (cADPR). In this study, by using a vesicular 45Ca2+ flux assay, we demonstrate that microsomal vesicles from red beet and cauliflower also respond to nicotinic acid adenine dinucleotide phosphate (NAADP), a Ca2+-releasing molecule recently described in marine invertebrates. NAADP potently mobilizes Ca2+ with a K1/2 = 96 nM from microsomes of nonvacuolar origin in red beet. Analysis of sucrose gradient-separated cauliflower microsomes revealed that the NAADP-sensitive Ca2+ pool was derived from the endoplasmic reticulum. This exclusively nonvacuolar location of the NAADP-sensitive Ca2+ pathway distinguishes it from the InsP3- and cADPR-gated pathways. Desensitization experiments revealed that homogenates derived from cauliflower tissue contained low levels of NAADP (125 pmol/mg) and were competent in NAADP synthesis when provided with the substrates NADP and nicotinic acid. NAADP-induced Ca2+ release is insensitive to heparin and 8-NH2-cADPR, specific inhibitors of the InsP3- and cADPR-controlled mechanisms, respectively. However, NAADP-induced Ca2+ release could be blocked by pretreatment with a subthreshold dose of NAADP, as previously observed in sea urchin eggs. Furthermore, the NAADP-gated Ca2+ release pathway is independent of cytosolic free Ca2+ and therefore incapable of operating Ca2+-induced Ca2+ release. In contrast to the sea urchin system, the NAADP-gated Ca2+ release pathway in plants is not blocked by L-type channel antagonists. The existence of multiple Ca2+ mobilization pathways and Ca2+ release sites might contribute to the generation of stimulus-specific Ca2+ signals in plant cells.

Structural basis for heterogeneous kinetics: Reengineering the hairpin ribozyme

The RNA cleavage reaction catalyzed by the hairpin ribozyme shows biphasic kinetics, and chase experiments show that the slow phase of the reaction results from reversible substrate binding to an inactive conformational isomer. To investigate the structural basis for the heterogeneous kinetics, we have developed an enzymatic RNA modification method that selectively traps substrate bound to the inactive conformer and allows the two forms of the ribozyme-substrate complex to be separated and analyzed by using both physical and kinetic strategies. The inactive form of the complex was trapped by the addition of T4 RNA ligase to a cleavage reaction, resulting in covalent linkage of the 5′ end of the substrate to the 3′ end of the ribozyme and in selective and quantitative ablation of the slow kinetic phase of the reaction. This result indicates that the inactive form of the ribozyme-substrate complex can adopt a conformation in which helices 2 and 3 are coaxially stacked, whereas the active form does not have access to this conformation, because of a sharp bend at the helical junction that presumably is stabilized by inter-domain tertiary contacts required for catalytic activity. These results were used to improve the activity of the hairpin ribozyme by designing new interfaces between the two domains, one containing a non-nucleotidic orthobenzene linkage and the other replacing the two-way junction with a three-way junction. Each of these modified ribozymes preferentially adopts the active conformation and displays improved catalytic efficiency.

Implications of macromolecular crowding for signal transduction and metabolite channeling

The effect of different total enzyme concentrations on the flux through the bacterial phosphoenolpyruvate:carbohydrate phosphotransferase system (PTS) in vitro was determined by measuring PTS-mediated carbohydrate phosphorylation at different dilutions of cell-free extract of Escherichia coli. The dependence of the flux on the protein concentration was more than linear but less than quadratic. The combined flux–response coefficient of the four enzymes constituting the glucose PTS decreased slightly from values of ≈1.8 with increasing protein concentrations in the assay. Addition of the macromolecular crowding agents polyethylene glycol (PEG) 6000 and PEG 35000 led to a sharper decrease in the combined flux–response coefficient, in one case to values of ≈1. PEG 6000 stimulated the PTS flux at lower protein concentrations and inhibited the flux at higher protein concentrations, with the transition depending on the PEG 6000 concentration. This suggests that macromolecular crowding decreases the dissociation rate constants of enzyme complexes. High concentrations of the microsolute glycerol did not affect the combined flux–response coefficient. The data could be explained with a kinetic model of macromolecular crowding in a two-enzyme group-transfer pathway. Our results suggest that, because of the crowded environment in the cell, the different PTS enzymes form complexes that live long on the time-scale of their turnover. The implications for the metabolic behavior and control properties of the PTS, and for the effect of macromolecular crowding on nonequilibrium processes, are discussed.

Kinetics of formation of hypoxanthine containing base pairs by HIV-RT: RNA template effects on the base substitution frequencies

Hypoxanthine (H), the deamination product of adenine, has been implicated in the high frequency of A to G transitions observed in retroviral and other RNA genomes. Although H·C base pairs are thermodynamically more stable than other H·N pairs, polymerase selection may be determined in part by kinetic factors. Therefore, the hypoxanthine induced substitution pattern resulting from replication by viral polymerases may be more complex than that predicted from thermodynamics. We have examined the steady-state kinetics of formation of base pairs opposite template H in RNA by HIV-RT, and for the incorporation of dITP during first- and second-strand synthesis. Hypoxanthine in an RNA template enhances the k2app for pairing with standard dNTPs by factors of 10–1000 relative to adenine at the same sequence position. The order of base pairing preferences for H in RNA was observed to be H·C >> H·T > H·A > H·G. Steady-state kinetics of insertion for all possible mispairs formed with dITP were examined on RNA and DNA templates of identical sequence. Insertion of dITP opposite all bases occurs 2–20 times more frequently on RNA templates. This bias for higher insertion frequencies on RNA relative to DNA templates is also observed for formation of mispairs at template A. This kinetic advantage afforded by RNA templates for mismatches and pairing involving H suggests a higher induction of mutations at adenines during first-strand synthesis by HIV-RT.

Solvent effects on the energy landscapes and folding kinetics of polyalanine

The effect of a solvation on the thermodynamics and kinetics of polyalanine (Ala12) is explored on the basis of its energy landscapes in vacuum and in an aqueous solution. Both energy landscapes are characterized by two basins, one associated with α-helical structures and the other with coil and β-structures of the peptide. In both environments, the basin that corresponds to the α-helical structure is considerably narrower than the basin corresponding to the β-state, reflecting their different contributions to the entropy of the peptide. In vacuum, the α-helical state of Ala12 constitutes the native state, in agreement with common helical propensity scales, whereas in the aqueous medium, the α-helical state is destabilized, and the β-state becomes the native state. Thus solvation has a dramatic effect on the energy landscape of this peptide, resulting in an inverted stability of the two states. Different folding and unfolding time scales for Ala12 in hydrophilic and hydrophobic chemical environments are caused by the higher entropy of the native state in water relative to vacuum. The concept of a helical propensity has to be extended to incorporate environmental solvent effects.

Kinetics of protein import into isolated Xenopus oocyte nuclei

An in vitro assay for nucleocytoplasmic transport was established in which signal-dependent protein import is reproduced faithfully by isolated purified nuclei. The assay permits the precise quantification of import kinetics and the discrimination between translocation through the nuclear envelope and intranuclear transport. Nuclei were manually isolated from Xenopus oocytes and after manual purification incubated with a medium containing a green fluorescent transport substrate, karyopherins α2 and β1, a red fluorescent control substrate, an energy mix and, for keeping an osmotic balance, 20% (wt/vol) BSA. Import of transport substrates into the nucleus and exclusion of the control substrate were monitored simultaneously by two-color confocal microscopy. Two widely differing import substrates were used: the recombinant protein P4K [480 kDa, four nuclear localization sequences (NLSs) per P4K tetramer], and NLS-BSA (90 kDa, 15 NLSs). The measurements suggested that import, at the specific conditions used in this study, consisted of two consecutive processes: (i) the rapid equilibration of the concentration difference across the nuclear envelope, a process involving binding and translocation of substrate by the nuclear pore complex , and (ii) the dissipation of the intranuclear concentration difference by diffusion.

Genetically altered AMPA-type glutamate receptor kinetics in interneurons disrupt long-range synchrony of gamma oscillation

Gamma oscillations synchronized between distant neuronal populations may be critical for binding together brain regions devoted to common processing tasks. Network modeling predicts that such synchrony depends in part on the fast time course of excitatory postsynaptic potentials (EPSPs) in interneurons, and that even moderate slowing of this time course will disrupt synchrony. We generated mice with slowed interneuron EPSPs by gene targeting, in which the gene encoding the 67-kDa form of glutamic acid decarboxylase (GAD67) was altered to drive expression of the α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) glutamate receptor subunit GluR-B. GluR-B is a determinant of the relatively slow EPSPs in excitatory neurons and is normally expressed at low levels in γ-aminobutyric acid (GABA)ergic interneurons, but at high levels in the GAD-GluR-B mice. In both wild-type and GAD-GluR-B mice, tetanic stimuli evoked gamma oscillations that were indistinguishable in local field potential recordings. Remarkably, however, oscillation synchrony between spatially separated sites was severely disrupted in the mutant, in association with changes in interneuron firing patterns. The congruence between mouse and model suggests that the rapid time course of AMPA receptor-mediated EPSPs in interneurons might serve to allow gamma oscillations to synchronize over distance.

The Xanthophyll Cycle Modulates the Kinetics of Nonphotochemical Energy Dissipation in Isolated Light-Harvesting Complexes, Intact Chloroplasts, and Leaves of Spinach1

We analyzed the kinetics of nonphotochemical quenching of chlorophyll fluorescence (qN) in spinach (Spinacia oleracea) leaves, chloroplasts, and purified light-harvesting complexes. The characteristic biphasic pattern of fluorescence quenching in dark-adapted leaves, which was removed by preillumination, was evidence of light activation of qN, a process correlated with the de-epoxidation state of the xanthophyll cycle carotenoids. Chloroplasts isolated from dark-adapted and light-activated leaves confirmed the nature of light activation: faster and greater quenching at a subsaturating transthylakoid pH gradient. The light-harvesting chlorophyll a/b-binding complexes of photosystem II were isolated from dark-adapted and light-activated leaves. When isolated from light-activated leaves, these complexes showed an increase in the rate of quenching in vitro compared with samples prepared from dark-adapted leaves. In all cases, the quenching kinetics were fitted to a single component hyperbolic function. For leaves, chloroplasts, and light-harvesting complexes, the presence of zeaxanthin was associated with an increased rate constant for the induction of quenching. We discuss the significance of these observations in terms of the mechanism and control of qN.

Non-Arrhenius kinetics for the loop closure of a DNA hairpin

Intramolecular chain diffusion is an elementary process in the conformational fluctuations of the DNA hairpin-loop. We have studied the temperature and viscosity dependence of a model DNA hairpin-loop by FRET (fluorescence resonance energy transfer) fluctuation spectroscopy (FRETfs). Apparent thermodynamic parameters were obtained by analyzing the correlation amplitude through a two-state model and are consistent with steady-state fluorescence measurements. The kinetics of closing the loop show non-Arrhenius behavior, in agreement with theoretical prediction and other experimental measurements on peptide folding. The fluctuation rates show a fractional power dependence (β = 0.83) on the solution viscosity. A much slower intrachain diffusion coefficient in comparison to that of polypeptides was derived based on the first passage time theory of SSS [Szabo, A., Schulten, K. & Schulten, Z. (1980) J. Chem. Phys. 72, 4350–4357], suggesting that intrachain interactions, especially stacking interaction in the loop, might increase the roughness of the free energy surface of the DNA hairpin-loop.

Kinetics and thermodynamics of T cell receptor– autoantigen interactions in murine experimental autoimmune encephalomyelitis

In the current study, cellular and molecular approaches have been used to analyze the biophysical nature of T cell receptor (TCR)–peptide MHC (pMHC) interactions for two autoreactive TCRs. These two TCRs recognize the N-terminal epitope of myelin basic protein (MBP1–11) bound to the MHC class II protein, I-Au, and are associated with murine experimental autoimmune encephalomyelitis. Mice transgenic for the TCRs have been generated and characterized in other laboratories. These analyses indicate that the mice either develop encephalomyelitis spontaneously (172.10 TCR) or only if immunized with autoantigen in adjuvant (1934.4 TCR). Here, we show that the 172.10 TCR binds MBP1–11:I-Au with a 4–5-fold higher affinity than the 1934.4 TCR. Consistent with the higher affinity, 172.10 T hybridoma cells are significantly more responsive to autoantigen than 1934.4 cells. The interaction of the 172.10 TCR with cognate ligand is more entropically unfavorable than that of the 1934.4 TCR, indicating that the 172.10 TCR undergoes greater conformational rearrangements upon ligand binding. The studies therefore suggest a correlation between the strength and plasticity of a TCR–pMHC interaction and the frequency of spontaneous disease in the corresponding TCR transgenic mice. The comparative analysis of these two TCRs has implications for understanding autoreactive T cell recognition and activation.

A statistical mechanical model for β-hairpin kinetics

Understanding the mechanism of protein secondary structure formation is an essential part of the protein-folding puzzle. Here, we describe a simple statistical mechanical model for the formation of a β-hairpin, the minimal structural element of the antiparallel β-pleated sheet. The model accurately describes the thermodynamic and kinetic behavior of a 16-residue, β-hairpin-forming peptide, successfully explaining its two-state behavior and apparent negative activation energy for folding. The model classifies structures according to their backbone conformation, defined by 15 pairs of dihedral angles, and is further simplified by considering only the 120 structures with contiguous stretches of native pairs of backbone dihedral angles. This single sequence approximation is tested by comparison with a more complete model that includes the 215 possible conformations and 15 × 215 possible kinetic transitions. Finally, we use the model to predict the equilibrium unfolding curves and kinetics for several variants of the β-hairpin peptide.