Kinetic stability as a mechanism for protease longevity

The folding of the extracellular serine protease, α-lytic protease (αLP; EC 3.4.21.12) reveals a novel mechanism for stability that appears to lead to a longer functional lifetime for the protease. For αLP, stability is based not on thermodynamics, but on kinetics. Whereas this has required the coevolution of a pro region to facilitate folding, the result has been the optimization of native-state properties independent of their consequences on thermodynamic stability. Structural and mutational data lead to a model for catalysis of folding in which the pro region binds to a conserved β-hairpin in the αLP C-terminal domain, stabilizing the folding transition state and the native state. The pro region is then proteolytically degraded, leaving the active αLP trapped in a metastable conformation. This metastability appears to be a consequence of pressure to evolve properties of the native state, including a large, highly cooperative barrier to unfolding, and extreme rigidity, that reduce susceptibility to proteolytic degradation. In a test of survival under highly proteolytic conditions, homologous mammalian proteases that have not evolved kinetic stability are much more rapidly degraded than αLP. Kinetic stability as a means to longevity is likely to be a mechanism conserved among the majority of extracellular bacterial pro-proteases and may emerge as a general strategy for intracellular eukaryotic proteases subject to harsh conditions as well.

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.

Optimizing the stability of single-chain proteins by linker length and composition mutagenesis

Linker length and composition were varied in libraries of single-chain Arc repressor, resulting in proteins with effective concentrations ranging over six orders of magnitude (10 μM–10 M). Linkers of 11 residues or more were required for biological activity. Equilibrium stability varied substantially with linker length, reaching a maximum for glycine-rich linkers containing 19 residues. The effects of linker length on equilibrium stability arise from significant and sometimes opposing changes in folding and unfolding kinetics. By fixing the linker length at 19 residues and varying the ratio of Ala/Gly or Ser/Gly in a 16-residue-randomized region, the effects of linker flexibility were examined. In these libraries, composition rather than sequence appears to determine stability. Maximum stability in the Ala/Gly library was observed for a protein containing 11 alanines and five glycines in the randomized region of the linker. In the Ser/Gly library, the most stable protein had seven serines and nine glycines in this region. Analysis of folding and unfolding rates suggests that alanine acts largely by accelerating folding, whereas serine acts predominantly to slow unfolding. These results demonstrate an important role for linker design in determining the stability and folding kinetics of single-chain proteins and suggest strategies for optimizing these parameters.

Conditional Synergism between Cryptochrome 1 and Phytochrome B Is Shown by the Analysis of phyA, phyB, and hy4 Simple, Double, and Triple Mutants in Arabidopsis

Wild-type or phyA, phyB, or hy4 mutant Arabidopsis seedlings lacking phytochrome A (phyA), phytochrome B (phyB), or cryptochrome 1 (cry1), respectively, and the double and triple mutants were used in combination with blue-light treatments given simultaneously with red or far-red light. We investigated the interaction between phytochromes and cry1 in the control of hypocotyl growth and cotyledon unfolding. Under conditions deficient for cry1 (short exposures to blue light) or phyB (far-red background), these photoreceptors acted synergistically: Under short exposures to blue light (3 h/d) added to a red-light background, cry1 activity required phyB (e.g. the hy4 mutant was taller than the wild type but the phyBhy4 mutant was not taller than the phyB mutant). Under prolonged exposures to blue light (24 h/d) added to a far-red light background, phyB activity required cry1 (e.g. the phyAphyB mutant was taller than the phyA mutant but the phyAphyBhy4 mutant was not taller than the phyAhy4 mutant). Under more favorable light inputs, i.e. prolonged exposures to blue light added to a red-light background, the effects of cry1 and phyB were independent. Thus, the synergism between phyB and cry1 is conditional. The effect of cry1 was not reduced by the phyA mutation under any tested light condition. Under continuous blue light the triple mutant phyAphyBhy4 showed reduced hypocotyl growth inhibition and cotyledon unfolding compared with the phyAphyB mutant. The action of cry1 in the phyAphyB double mutant was higher under the red-light than the far-red-light background, indicating a synergistic interaction between cry1 and phytochromes C, D, or E; however, a residual action of cry1 independent of any phytochrome is likely to occur.

Genetic Interactions between Phytochrome A, Phytochrome B, and Cryptochrome 1 during Arabidopsis Development1

Single, double, and triple null combinations of Arabidopsis mutants lacking the photoreceptors phytochrome (phy) A (phyA-201), phyB (phyB-5), and cryptochrome (cry) 1 (hy4-2.23n) were examined for de-etiolation responses in high-fluence red, far-red, blue, and broad-spectrum white light. Cotyledon unhooking, unfolding, and expansion, hypocotyl growth, and the accumulation of chlorophylls and anthocyanin in 5-d-old seedlings were measured under each light condition and in the dark. phyA was the major photoreceptor/effector for most far-red-light responses, although phyB and cry1 modulated anthocyanin accumulation in a phyA-dependent manner. phyB was the major photoreceptor in red light, although cry1 acted as a phyA/phyB-dependent modulator of chlorophyll accumulation under these conditions. All three photoreceptors contributed to most blue light deetiolation responses, either redundantly or additively; however, phyB acted as a modulator of cotyledon expansion dependent on the presence of cry1. As reported previously, flowering time in long days was promoted by phyA and inhibited by phyB, with each suppressing the other's effect. In addition to the effector/modulator relationships described above, measurements of hypocotyls from blue-light-grown seedlings demonstrated phytochrome activity in blue light and cry1 activity in a phyAphyB mutant background.

Different Phototransduction Kinetics of Phytochrome A and Phytochrome B in Arabidopsis thaliana1

The kinetics of phototransduction of phytochrome A (phyA) and phytochrome B (phyB) were compared in etiolated Arabidopsis thaliana seedlings. The responses of hypocotyl growth, cotyledon unfolding, and expression of a light-harvesting chlorophyll a/b-binding protein of the photosystem II gene promoter fused to the coding region of β-glucuronidase (used as a reporter enzyme) were mediated by phyA under continuous far-red light (FR) and by phyB under continuous red light (R). The seedlings were exposed hourly either to n min of FR followed by 60 minus n min in darkness or to n min of R, 3 min of FR (to back-convert phyB to its inactive form), and 57 minus n min of darkness. For the three processes investigated here, the kinetics of phototransduction of phyB were faster than that of phyA. For instance, 15 min R h−1 (terminated with a FR pulse) were almost as effective as continuous R, whereas 15 min of FR h−1 caused less than 30% of the effect of continuous FR. This difference is interpreted in terms of divergence of signal transduction pathways downstream from phyA and phyB.

Cytochrome c′ folding triggered by electron transfer: Fast and slow formation of four-helix bundles

Reduced (FeII) Rhodopseudomonas palustris cytochrome c′ (Cyt c′) is more stable toward unfolding ([GuHCl]1/2 = 2.9(1) M) than the oxidized (FeIII) protein ([GuHCl]1/2 = 1.9(1) M). The difference in folding free energies (ΔΔGf° = 70 meV) is less than half of the difference in reduction potentials of the folded protein (100 mV vs. NHE) and a free heme in aqueous solution (≈−150 mV). The spectroscopic features of unfolded FeII–Cyt c′ indicate a low-spin heme that is axially coordinated to methionine sulfur (Met-15 or Met-25). Time-resolved absorption measurements after CO photodissociation from unfolded FeII(CO)–Cyt c′ confirm that methionine can bind to the ferroheme on the microsecond time scale [kobs = 5(2) × 104 s−1]. Protein folding was initiated by photoreduction (two-photon laser excitation of NADH) of unfolded FeIII–Cyt c′ ([GuHCl] = 2.02–2.54 M). Folding kinetics monitored by heme absorption span a wide time range and are highly heterogeneous; there are fast-folding (≈103 s−1), intermediate-folding (102–101 s−1), and slow-folding (10−1 s−1) populations, with the last two likely containing methionine-ligated (Met-15 or Met-25) ferrohemes. Kinetics after photoreduction of unfolded FeIII–Cyt c′ in the presence of CO are attributable to CO binding [1.4(6) × 103 s−1] and FeII(CO)–Cyt c′ folding [2.8(9) s−1] processes; stopped-flow triggered folding of FeIII–Cyt c′ (which does not contain a protein-derived sixth ligand) is adequately described by a single kinetics phase with an estimated folding time constant of ≈4 ms [ΔGf° = −33(3) kJ mol−1] at zero denaturant.

beta-Lactamase binds to GroEL in a conformation highly protected against hydrogen/deuterium exchange.

Escherichia coli RTEM beta-lactamase reversibly forms a stable complex with GroEL, devoid of any enzymatic activity, at 48 degrees C. When beta-lactamase is diluted from this complex into denaturant solution, its unfolding rate is identical to that from the native state, while the unfolding rate from the molten globule state is too fast to be measured. Electrospray mass spectrometry shows that the rate of proton exchange in beta-lactamase in the complex at 48 degrees C is slower than in the absence of GroEL at the same temperature, and resembles the exchange of the native state at 25 degrees C. Similarly, the final number of protected deuterons is higher in the presence of GroEL than in its absence. We conclude that, for beta-lactamase, a state with significant native structure is bound to GroEL. Thus, different proteins are recognized by GroEL in very different states, ranging from totally unfolded to native-like, and this recognition may depend on which state can provide sufficient accessible hydrophobic amino acids in a suitably clustered arrangement. Reversible binding of native-like states with hydrophobic patches may be an important property of GroEL to protect the cell from aggregating protein after heat-shock.

What determines the folding of the chromatin fiber?

In this review, we attempt to summarize, in a critical manner, what is currently known about the processes of condensation and decondensation of chromatin fibers. We begin with a critical analysis of the possible mechanisms for condensation, considering both old and new evidence as to whether the linker DNA between nucleosomes bends or remains straight in the condensed structure. Concluding that the preponderance of evidence is for straight linkers, we ask what other fundamental process might allow condensation, and argue that there is evidence for linker histone-induced contraction of the internucleosome angle, as salt concentration is raised toward physiological levels. We also ask how certain specific regions of chromatin can become decondensed, even at physiological salt concentration, to allow transcription. We consider linker histone depletion and acetylation of the core histone tails, as possible mechanisms. On the basis of recent evidence, we suggest a unified model linking targeted acetylation of specific genomic regions to linker histone depletion, with unfolding of the condensed fiber as a consequence.

Evidence for aminoacylation-induced conformational changes in human mitochondrial tRNAs.

Analysis by acid polyacrylamide/urea gel electrophoresis of 14 individual mitochondrial tRNAs (mt-tRNAs) from human cells has revealed a variable decrease in mobility of the aminoacylated relative to the nonacylated form, with the degree of separation of the two forms not being correlated with the mass, polar character, or charge of the amino acid. Separation of the charged and uncharged species has been found to be independent of tRNA denaturation, being observed also in the absence of urea. In another approach, electrophoresis through a perpendicular denaturing gradient gel of several individual mt-tRNAs has shown a progressive unfolding of the tRNA with increasing denaturant concentration, which is consistent with an initial disruption of tertiary interactions, followed by the sequential melting of the four stems of the cloverleaf structure. A detailed analysis of the unfolding process of charged and uncharged tRNALys and tRNALeu(UUR) has revealed that the separation of the two forms of these tRNAs persisted throughout the almost entire range of denaturant concentrations used and was lost upon denaturation of the last helical domain(s), which most likely included the amino acid acceptor stem. These observations strongly suggest that the electrophoretic retardation of the charged species reflects an aminoacylation-induced conformational change of the 3'-end of these mt-tRNAs, with possible significant implications in connection with the known role of the acceptor end in tRNA interactions with the ribosomal peptidyl transferase center and the elongation factor Tu.

The N-terminal domain of Escherichia coli enzyme I of the phosphoenolpyruvate/glycose phosphotransferase system: molecular cloning and characterization.

The bacterial phosphoenolpyruvate/glycose phosphotransferase system (PTS) comprises a group of proteins that catalyze the transfer of the phosphoryl group from phosphoenolpyruvate (PEP) to sugars concomitant with their translocation. The first two steps of the phosphotransfer sequence are PEP <--> Enzyme I (EI) <--> HPr (the histidine-containing phosphocarrier protein). We have proposed that many functions of the PTS are regulated by EI, which undergoes a monomer/dimer transition. EI monomer (63.5 kDa) comprises two major domains: a flexible C-terminal domain (EI-C) and a protease-resistant, structurally stable N-terminal domain (EI-N) containing the active site His. Trypsin treatment of Salmonella typhimurium EI yielded EI-N, designated EI-N(t). Homogeneous recombinant Escherichia coli EI-N [i.e., EI-N(r)], has now been prepared in quantity, shows the expected thermodynamic unfolding properties and, similarly to EI-N(t), is phosphorylated by phospho-HPr, but not by PEP. In addition, binding of EI-N(r) to HPr was studied by isothermal titration calorimetry: K/a = 1.4 x 10(5) M(-1) and delta H = +8.8 kcal x mol(-1). Both values are comparable to those for HPr binding to intact EI. Fluorescence anisotropy [dansyl-EI-N(r)] and gel filtration of EI-N(r) show that it does not dimerize. These results emphasize the role of EI-C in dimerization and the regulation of intact EI.

A new isoleucine substitution of Val-20 in transthyretin tetramers selectively impairs dimer-dimer contacts and causes systemic amyloidosis.

The most frequent form of inherited amyloidoses is associated with mutations in the transthyretin (TTR) gene coding for 127-amino acid residues of four identical, noncovalently linked subunits that form a pair of dimers in the plasma protein complex. Amyloid fibrils containing the variant and to a lesser extent the wild-type form of the TTR molecule are deposited in various organs, including peripheral nerves and the myocardium, with polyneuropathy and cardiomyopathy as major clinical manifestations. So far, more than 40 distinct amino acid substitutions distributed throughout the TTR sequence over 30 positions have been found to be correlated with an increased amyloidogenicity of TTR. Most of these amyloidogenic amino acid substitutions are suspected to alter the conformation and stability of the monomer. Here we identify and characterize by protein and DNA analysis a novel amyloidogenic Val-20 to Ile mutation in a German three-generation family. The index patient suffered from severe amyloid cardiomyopathy at the age of 60. Conformational stability and unfolding behavior of the Ile-20 monomer in urea gradients was found to be almost indistinguishable from that of wild-type TTR. In contrast, tetramer stability was significantly reduced in agreement with the expected change in the interactions between the two opposing dimers via the side chain of Ile-20. Our observations provide strong evidence for the view that amyloidogenic amino acid substitutions in TTR facilitate the conversion of tetrameric TTR complexes into those conformational intermediates of the TTR folding pathway that have an intrinsic amyloidogenic potential.

Toward a mechanism for GroEL.GroES chaperone activity: an ATPase-gated and -pulsed folding and annealing cage.

Free GroEL binds denatured proteins very tightly: it retards the folding of barnase 400-fold and catalyzes unfolding fluctuations in native barnase and its folding intermediate. GroEL undergoes an allosteric transition from its tight-binding T-state to a weaker binding R-state on the cooperative binding of nucleotides (ATP/ADP) and GroES. The preformed GroEL.GroES.nucleotide complex retards the folding of barnase by only a factor of 4, and the folding rate is much higher than the ATPase activity that releases GroES from the complex. Binding of GroES and nucleotides to a preformed GroEL.denatured-barnase complex forms an intermediately fast-folding complex. We propose the following mechanism for the molecular chaperone. Denatured proteins bind to the resting GroEL.GroES.nucleotide complex. Fast-folding proteins are ejected as native structures before ATP hydrolysis. Slow-folding proteins enter chaperoning cycles of annealing and folding after the initial ATP hydrolysis. This step causes transient release of GroES and formation of the GroEL.denatured-protein complexes with higher annealing potential. The intermediately fast-folding complex is formed on subsequent rebinding of GroES. The ATPase activity of GroEL.GroES is thus the gatekeeper that selects for initial entry of slow-folding proteins to the chaperone action and then pumps successive transitions from the faster-folding R-states to the tighter-binding/stronger annealing T-states. The molecular chaperone acts as a combination of folding cage and an annealing machine.

Least activation path for protein folding: investigation of staphylococcal nuclease folding by stopped-flow circular dichroism.

Is the pathway of protein folding determined by the relative stability of folding intermediates, or by the relative height of the activation barriers leading to these intermediates? This is a fundamental question for resolving the Levinthal paradox, which stated that protein folding by a random search mechanism would require a time too long to be plausible. To answer this question, we have studied the guanidinium chloride (GdmCl)-induced folding/unfolding of staphylococcal nuclease [(SNase, formerly EC 3.1.4.7; now called microbial nuclease or endonuclease, EC 3.1.31.1] by stopped-flow circular dichroism (CD) and differential scanning microcalorimetry (DSC). The data show that while the equilibrium transition is a quasi-two-state process, kinetics in the 2-ms to 500-s time range are triphasic. Data support the sequential mechanism for SNase folding: U3 <--> U2 <--> U1 <--> N0, where U1, U2, and U3 are substates of the unfolded protein and N0 is the native state. Analysis of the relative population of the U1, U2, and U3 species in 2.0 M GdmCl gives delta-G values for the U3 --> U2 reaction of +0.1 kcal/mol and for the U2 --> U1 reaction of -0.49 kcal/mol. The delta-G value for the U1 --> N0 reaction is calculated to be -4.5 kcal/mol from DSC data. The activation energy, enthalpy, and entropy for each kinetic step are also determined. These results allow us to make the following four conclusions. (i) Although the U1, U2, and U3 states are nearly isoenergetic, no random walk occurs among them during the folding. The pathway of folding is unique and sequential. In other words, the relative stability of the folding intermediates does not dictate the folding pathway. Instead, the folding is a descent toward the global free-energy minimum of the native state via the least activation path in the vast energy landscape. Barrier avoidance leads the way, and barrier height limits the rate. Thus, the Levinthal paradox is not applicable to the protein-folding problem. (ii) The main folding reaction (U1 --> N0), in which the peptide chain acquires most of its free energy (via van der Waals' contacts, hydrogen bonding, and electrostatic interactions), is a highly concerted process. These energy-acquiring events take place in a single kinetic phase. (iii) U1 appears to be a compact unfolded species; the rate of conversion of U2 to U1 depends on the viscosity of solution. (iv) All four relaxation times reported here depend on GdmCl concentrations: it is likely that none involve the cis/trans isomerization of prolines. Finally, a mechanism is presented in which formation of sheet-like chain conformations and a hydrophobic condensation event precede the main-chain folding reaction.

The Drosophila melanogaster dodo (dod) gene, conserved in humans, is functionally interchangeable with the ESS1 cell division gene of Saccharomyces cerevisiae.

We have sequenced the region of DNA adjacent to and including the flightless (fli) gene of Drosophila melanogaster and molecularly characterized four transcription units within it, which we have named tweety (twe), flightless (fli), dodo (dod), and penguin (pen). We have performed deletion and transgenic analysis to determine the consequences of the quadruple gene removal. Only the flightless gene is vital to the organism; the simultaneous absence of the other three allows the overriding majority of individuals to develop to adulthood and to fly normally. These gene deletion results are evaluated in the context of the redundancy and degeneracy inherent in many genetic networks. Our cDNA analyses and data-base searches reveal that the predicted dodo protein has homologs in other eukaryotes and that it is made up of two different domains. The first, designated WW, is involved in protein-protein interactions and is found in functionally diverse proteins including human dystrophin. The second is involved in accelerating protein folding and unfolding and is found in Escherichia coli in a new family of peptidylprolyl cis-trans isomerases (PPIases; EC 5.2.1.8). In eukaryotes, PPIases occur in the nucleus and the cytoplasm and can form stable associations with transcription factors, receptors, and kinases. Given this particular combination of domains, the dodo protein may well participate in a multisubunit complex involved in the folding and activation of signaling molecules. When we expressed the dodo gene product in Saccharomyces cerevisiae, it rescued the lethal phenotype of the ESS1 cell division gene.