Thermodynamic stability of wild-type and mutant p53 core domain

Some 50% of human cancers are associated with mutations in the core domain of the tumor suppressor p53. Many mutations are thought just to destabilize the protein. To assess this and the possibility of rescue, we have set up a system to analyze the stability of the core domain and its mutants. The use of differential scanning calorimetry or spectroscopy to measure its melting temperature leads to irreversible denaturation and aggregation and so is useful as only a qualitative guide to stability. There are excellent two-state denaturation curves on the addition of urea that may be analyzed quantitatively. One Zn2+ ion remains tightly bound in the holo-form of p53 throughout the denaturation curve. The stability of wild type is 6.0 kcal (1 kcal = 4.18 kJ)/mol at 25°C and 9.8 kcal/mol at 10°C. The oncogenic mutants R175H, C242S, R248Q, R249S, and R273H are destabilized by 3.0, 2.9, 1.9, 1.9, and 0.4 kcal/mol, respectively. Under certain denaturing conditions, the wild-type domain forms an aggregate that is relatively highly fluorescent at 340 nm on excitation at 280 nm. The destabilized mutants give this fluorescence under milder denaturation conditions.

Direct interaction of flagellin termini essential for polymorphic ability of flagellar filament

We report the structures of flagellar filaments reconstituted from various flagellins with small terminal truncations. Flagellins from Salmonella typhimurium strains SJW1103 (wild type), SJW1660, and SJW1655 were used, which form a left-handed supercoil, the L- and R-type straight forms, respectively. Structure analyses were done by electron cryomicroscopy and helical image reconstruction with a help of x-ray fiber diffraction for determining precise helical symmetries. Truncation of either terminal region, irrespective of the original flagellin species, results in a straight filament having a helical symmetry distinct either from the L- or R-type. This filament structure is named Lt-type. Although the local subunit packing is similar in all three types, a close comparison shows that the Lt-type packing is almost identical to the R-type but distinct from the L-type, which demonstrates the strong two-state preference of the subunit interactions. The structure clearly suggests that both termini are located in the inner tube of the concentric double-tubular structure of the filament core, and their proper interaction is responsible for the correct folding of fairly large terminal regions that form the inner tube. The double tubular structure appears to be essential for the polymorphic ability of flagellar filaments, which is required for the swimming–tumbling of bacterial taxis.

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.

Memory landscapes of single-enzyme molecules

Immobilized single horseradish peroxidase enzymes were observed by confocal fluorescence spectroscopy during catalysis of the oxidation reaction of the nonfluorescent dihydrorhodamine 6G substrate into the highly fluorescent product rhodamine 6G. By extracting only the non-Markovian behavior of the spectroscopic two-state process of enzyme-product complex formation and release, memory landscapes were generated for single-enzyme molecules. The memory landscapes can be used to discriminate between different origins of stretched exponential kinetics that are found in the first-order correlation analysis. Memory landscapes of single-enzyme data shows oscillations that are expected in a single-enzyme system that possesses a set of transient states. Alternative origins of the oscillations may not, however, be ruled out. The data and analysis indicate that substrate interaction with the enzyme selects a set of conformational substates for which the enzyme is active.

Geometric constraints explain much of the species richness pattern in African birds

The world contains boundaries (e.g., continental edge for terrestrial taxa) that impose geometric constraints on the distribution of species ranges. Thus, contrary to traditional thinking, the expected species richness pattern in absence of ecological or physiographical factors is unlikely to be uniform. Species richness has been shown to peak in the middle of a bounded one-dimensional domain, even in the absence of ecological or physiographical factors. Because species ranges are not linear, an extension of the approach to two dimensions is necessary. Here we present a two-dimensional null model accounting for effects of geometric constraints. We use the model to examine the effects of continental edge on the distribution of terrestrial animals in Africa and compare the predictions with the observed pattern of species richness in birds endemic to the continent. Latitudinal, longitudinal, and two-dimensional patterns of species richness are predicted well from the modeled null effects alone. As expected, null effects are of high significance for wide ranging species only. Our results highlight the conceptual significance of an until recently neglected constraint from continental shape alone and support a more cautious analysis of species richness patterns at this scale.

Influence of cisplatin intrastrand crosslinking on the conformation, thermal stability, and energetics of a 20-mer DNA duplex.

cis-Diamminedichloroplatinum(II) (cisplatin) is a widely used anticancer drug that binds to and crosslinks DNA. The major DNA adduct of the drug results from coordination of two adjacent guanine bases to platinum to form the intrastrand crosslink cis-[Pt(NH3)2[d(GpG)-N7(1), -N7(2)]] (cis-Pt-GG). In the present study, spectroscopic and calorimetric techniques were employed to characterize the influence of this crosslink on the conformation, thermal stability, and energetics of a site-specifically platinated 20-mer DNA duplex. CD spectroscopic and thermal denaturation data revealed that the crosslink alters the structure of the host duplex, consistent with a shift from a B-like to an A-like conformation; lowers its thermal stability by approximately 9 degrees C; and reduces its thermodynamic stability by 6.3 kcal/mol at 25 degrees C, most of which is enthalpic in origin; but it does not alter the two-state melting behavior exhibited by the parent, unmodified duplex, despite the significant crosslink-induced changes noted above. The energetic consequences of the cis-Pt-GG crosslink are discussed in relation to the structural perturbations it induces in DNA and to how these crosslink-induced perturbations might modulate protein binding.

Slip complexity in dynamic models of earthquake faults.

We summarize recent evidence that models of earthquake faults with dynamically unstable friction laws but no externally imposed heterogeneities can exhibit slip complexity. Two models are described here. The first is a one-dimensional model with velocity-weakening stick-slip friction; the second is a two-dimensional elastodynamic model with slip-weakening friction. Both exhibit small-event complexity and chaotic sequences of large characteristic events. The large events in both models are composed of Heaton pulses. We argue that the key ingredients of these models are reasonably accurate representations of the properties of real faults.

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.

Kinetic traps in lysozyme folding.

Folding of lysozyme from hen egg white was investigated by using interrupted refolding experiments. This method makes use of a high energy barrier between the native state and transient folding intermediates, and, in contrast to conventional optical techniques, it enables one to specifically monitor the amount of native molecules during protein folding. The results show that under strongly native conditions lysozyme can refold on parallel pathways. The major part of the lysozyme molecules (86%) refold on a slow kinetic pathway with well-populated partially folded states. Additionally, 14% of the molecules fold faster. The rate constant of formation of native molecules on the fast pathway corresponds well to the rate constant expected for folding to occur by a two-state process without any detectable intermediates. The results suggest that formation of the native state for the major fraction of lysozyme molecules is retarded compared with the direct folding process. Partially structured intermediates that transiently populate seem to be kinetically trapped in a conformation that can only slowly reach the native structure.

Submillisecond folding of monomeric lambda repressor.

The folding kinetics of a truncated form of the N-terminal domain of phage lambda repressor [lambda 6-85] has been investigated by using the technique of dynamic NMR. lambda 6-85 has been shown previously to fold in a purely two-state fashion. This allows the determination of folding and unfolding rates from simulation of the exchange-broadened aromatic resonances of Tyr-22. The folding kinetics were determined over a range of 1.35 to 3.14 M urea. The urea dependence of both folding and unfolding rate constants is exponential, suggesting that the rate-determining step is invariant at the urea concentrations studied. The folding and unfolding rates extrapolated to 0 M urea at 37 degrees C are 3600 +/- 400 s-1 and 27 +/- 6 s-1, respectively. The observed lambda 6-85 folding rate constant exceeds that of other fast-folding globular proteins by a factor of 14-54. The urea dependence of the folding and unfolding rate constants suggests that the transition state of the rate-determining step is considerably more exposed to solvent than previously studied protein-folding transition states. The surprising rapidity of lambda 6-85 folding and unfolding may be the consequence of its all-helical secondary structure. These kinetic results clearly demonstrate that all of the fundamental events of protein folding can occur on the submillisecond time scale.