Peptide interfacial adsorption is kinetically limited by the thermodynamic stability of self association

We present a study of the adsorption of two peptides at the octane–water interface. The first peptide, Lac21, exists in mixed monomer–tetramer equilibrium in bulk solution with an appreciable monomer concentration. The second peptide, Lac28, exists as a tetramer in solution, with minimal exposed hydrophobic surface. A kinetic limitation to interfacial adsorption exists for Lac28 at moderate to high surface coverage that is not observed for Lac21. We estimate the potential energy barrier for Lac28 adsorption to be 42 kJ/mol and show that this is comparable to the expected free energy barrier for tetramer dissociation. This finding suggests that, at moderate to high surface coverage, adsorption is kinetically limited by the availability of interfacially active monomeric “domains” in the subinterfacial region. We also show how the commonly used empirical equation for protein adsorption dynamics can be used to estimate the potential energy barrier for adsorption. Such an approach is shown to be consistent with a formal description of diffusion–adsorption, provided a large potential energy barrier exists. This work demonstrates that the dynamics of interfacial adsorption depend on protein thermodynamic stability, and hence structure, in a quantifiable way.

The 3′ untranslated region of a rice α-amylase gene functions as a sugar-dependent mRNA stability determinant

In plants, sugar feedback regulation provides a mechanism for control of carbohydrate allocation and utilization among tissues and organs. The sugar repression of α-amylase gene expression in rice provides an ideal model for studying the mechanism of sugar feedback regulation. We have shown previously that sugar repression of α-amylase gene expression in rice suspension cells involves control of both transcription rate and mRNA stability. The α-amylase mRNA is significantly more stable in sucrose-starved cells than in sucrose-provided cells. To elucidate the mechanism of sugar-dependent mRNA turnover, we have examined the effect of αAmy3 3′ untranslated region (UTR) on mRNA stability by functional analyses in transformed rice suspension cells. We found that the entire αAmy3 3′ UTR and two of its subdomains can independently mediate sugar-dependent repression of reporter mRNA accumulation. Analysis of reporter mRNA half-lives demonstrated that the entire αAmy3 3′ UTR and the two subdomains each functioned as a sugar-dependent destabilizing determinant in the turnover of mRNA. Nuclear run-on transcription analysis further confirmed that the αAmy3 3′ UTR and the two subdomains did not affect the transcription rate of promoter. The identification of sequence elements in the α-amylase mRNA that dictate the differential stability has very important implications for the study of sugar-dependent mRNA decay mechanisms.

The thermodynamic origin of the stability of a thermophilic ribozyme

Understanding the mechanism of thermodynamic stability of an RNA structure has significant implications for the function and design of RNA. We investigated the equilibrium folding of a thermophilic ribozyme and its mesophilic homologue by using hydroxyl radical protection, small-angle x-ray scattering, and circular dichroism. Both RNAs require Mg2+ to fold to their native structures that are very similar. The stability is measured as a function of Mg2+ and urea concentrations at different temperatures. The enhanced stability of the thermophilic ribozyme primarily is derived from a tremendous increase in the amount of structure formed in the ultimate folding transition. This increase in structure formation and cooperativity arises because the penultimate and the ultimate folding transitions in the mesophilic ribozyme become linked into a single transition in the folding of the thermophilic ribozyme. Therefore, the starting point, or reference state, for the transition to the native, functional thermophilic ribozyme is significantly less structured. The shift in the reference state, and the resulting increase in folding cooperativity, is likely due to the stabilization of selected native interactions that only form in the ultimate transition. This mechanism of using a less structured intermediate and increased cooperativity to achieve higher functional stability for tertiary RNAs is fundamentally different from that commonly proposed to explain the increased stability of thermophilic proteins.

The challenge of paleoecological stasis: reassessing sources of evolutionary stability.

The paleontological record of the lower and middle Paleozoic Appalachian foreland basin demonstrates an unprecedented level of ecological and morphological stability on geological time scales. Some 70-80% of fossil morphospecies within assemblages persist in similar relative abundances in coordinated packages lasting as long as 7 million years despite evidence for environmental change and biotic disturbances. These intervals of stability are separated by much shorter periods of ecological and evolutionary change. This pattern appears widespread in the fossil record. Existing concepts of the evolutionary process are unable to explain this uniquely paleontological observation of faunawide coordinated stasis. A principle of evolutionary stability that arises from the ecosystem is explored here. We propose that hierarchical ecosystem theory, when extended to geological time scales, can explain long-term paleoecological stability as the result of ecosystem organization in response to high-frequency disturbance. The accompanying stability of fossil morphologies results from "ecological locking," in which selection is seen as a high-rate response of populations that is hierarchically constrained by lower-rate ecological processes. When disturbance exceeds the capacity of the system, ecological crashes remove these higher-level constraints, and evolution is free to proceed at high rates of directional selection during the organization of a new stable ecological hierarchy.

Types of evolutionary stability and the problem of cooperation.

The evolutionary stability of cooperation is a problem of fundamental importance for the biological and social sciences. Different claims have been made about this issue: whereas Axelrod and Hamilton's [Axelrod, R. & Hamilton, W. (1981) Science 211, 1390-1398] widely recognized conclusion is that cooperative rules such as "tit for tat" are evolutionarily stable strategies in the iterated prisoner's dilemma (IPD), Boyd and Lorberbaum [Boyd, R. & Lorberbaum, J. (1987) Nature (London) 327, 58-59] have claimed that no pure strategy is evolutionarily stable in this game. Here we explain why these claims are not contradictory by showing in what sense strategies in the IPD can and cannot be stable and by creating a conceptual framework that yields the type of evolutionary stability attainable in the IPD and in repeated games in general. Having established the relevant concept of stability, we report theorems on some basic properties of strategies that are stable in this sense. We first show that the IPD has "too many" such strategies, so that being stable does not discriminate among behavioral rules. Stable strategies differ, however, on a property that is crucial for their evolutionary survival--the size of the invasion they can resist. This property can be interpreted as a strategy's evolutionary robustness. Conditionally cooperative strategies such as tit for tat are the most robust. Cooperative behavior supported by these strategies is the most robust evolutionary equilibrium: the easiest to attain, and the hardest to disrupt.

An autonomous folding unit mediates the assembly of two-stranded coiled coils

Subunit oligomerization of many proteins is mediated by coiled-coil domains. Although the basic features contributing to the thermodynamic stability of coiled coils are well understood, the mechanistic details of their assembly have not yet been dissected. Here we report a 13-residue sequence pattern that occurs with limited sequence variations in many two-stranded coiled coils and that is absolutely required for the assembly of the Dictyostelium discoideum actin-bundling protein cortexillin I and the yeast transcriptional activator GCN4. The functional relationship between coiled-coil “trigger” sequences was manifested by replacing the intrinsic trigger motif of GCN4 with the related sequence from cortexillin I. We demonstrate that these trigger sequences represent autonomous helical folding units that, in contrast to arbitrarily chosen heptad repeats, can mediate coiled-coil formation. Aside from being of general interest for protein folding, trigger motifs should be of particular importance in the protein de novo design.

A retro-inverso peptide corresponding to the GH loop of foot-and-mouth disease virus elicits high levels of long-lasting protective neutralizing antibodies

Peptides corresponding to the immunodominant loop located at residues 135–158 on capsid protein VP1 of foot-and-mouth disease virus (FMDV) generally elicit high levels of anti-peptide and virus-neutralizing antibodies. In some instances, however, the level of neutralizing antibodies is low or even negligible, even though the level of anti-peptide antibodies is high. We have shown previously that the antigenic activity of peptide 141–159 of VP1 of a variant of serotype A can be mimicked by a retro-inverso (all-d retro or retroenantio) peptide analogue. This retro-inverso analogue induced greater and longer-lasting antibody titers than did the corresponding l-peptide. We now show that a single inoculation of the retro-inverso analogue elicits high levels of neutralizing antibodies that persist longer than those induced against the corresponding l-peptide and confer substantial protection in guinea pigs challenged with the cognate virus. In view of the high stability to proteases of retro-inverso peptide analogues and their enhanced immunogenicity, these results have practical relevance in designing potential peptide vaccines.

Identification of a Human VPF/VEGF 3′ Untranslated Region Mediating Hypoxia-induced mRNA Stability

Hypoxia is a prominent feature of malignant tumors that are characterized by angiogenesis and vascular hyperpermeability. Vascular permeability factor/vascular endothelial growth factor (VPF/VEGF) has been shown to be up-regulated in the vicinity of necrotic tumor areas, and hypoxia potently induces VPF/VEGF expression in several tumor cell lines in vitro. Here we report that hypoxia-induced VPF/VEGF expression is mediated by increased transcription and mRNA stability in human M21 melanoma cells. RNA-binding/electrophoretic mobility shift assays identified a single 125-bp AU-rich element in the 3′ untranslated region that formed hypoxia-inducible RNA-protein complexes. Hypoxia-induced expression of chimeric luciferase reporter constructs containing this 125-bp AU-rich hypoxia stability region were significantly higher than constructs containing an adjacent 3′ untranslated region element without RNA-binding activity. Using UV-cross-linking studies, we have identified a series of hypoxia-induced proteins of 90/88 kDa, 72 kDa, 60 kDa, 56 kDa, and 46 kDa that bound to the hypoxia stability region element. The 90/88-kDa and 60-kDa species were specifically competed by excess hypoxia stability region RNA. Thus, increased VPF/VEGF mRNA stability induced by hypoxia is mediated, at least in part, by specific interactions between a defined mRNA stability sequence in the 3′ untranslated region and distinct mRNA-binding proteins in human tumor cells.

Staphylococcal protein A: Unfolding pathways, unfolded states, and differences between the B and E domains

We present multiple native and denaturation simulations of the B and E domains of the three-helix bundle protein A, totaling 60 ns. The C-terminal helix (H3) consistently denatures later than either of the other two helices and contains residual helical structure in the denatured state. These results are consistent with experiments suggesting that the isolated H3 fragment is more stable than H1 and H2 and that H3 forms early in folding. Interestingly, the denatured state of the B domain is much more compact than that of the E domain. This sequence-dependent effect on the dimensions of the denatured state and the lack of correlation with structure suggest that the radius of gyration can be a misleading reaction coordinate for unfolding/folding. Various unfolding and refolding events are observed in the denaturation simulations. In some cases, the transitions are facilitated through interactions with other portions of the protein—contact-assisted helix formation. In the native simulations, the E domain is very stable: after 6 ns, the Cα root-mean-square deviation from the starting structure is less than 1.4 Å. In contrast, the native state of the B domain deviates more and its inter-helical angles fluctuate. In apparent contrast, we note that the B domain is thermodynamically more stable than the E domain. The simulations suggest that the increased stability of the B domain may be due to heightened mobility, and therefore entropy, in the native state and decreased mobility/entropy in the more compact denatured state.

A quantitative method for evaluating the stabilities of nucleic acids

We report a general method for screening, in solution, the impact of deviations from canonical Watson-Crick composition on the thermodynamic stability of nucleic acid duplexes. We demonstrate how fluorescence resonance energy transfer (FRET) can be used to detect directly free energy differences between an initially formed “reference” duplex (usually a Watson-Crick duplex) and a related “test” duplex containing a lesion/alteration of interest (e.g., a mismatch, a modified, a deleted, or a bulged base, etc.). In one application, one titrates into a solution containing a fluorescently labeled, FRET-active, reference duplex, an unlabeled, single-stranded nucleic acid (test strand), which may or may not compete successfully to form a new duplex. When a new duplex forms by strand displacement, it will not exhibit FRET. The resultant titration curve (normalized fluorescence intensity vs. logarithm of test strand concentration) yields a value for the difference in stability (free energy) between the newly formed, test strand-containing duplex and the initial reference duplex. The use of competitive equilibria in this assay allows the measurement of equilibrium association constants that far exceed the magnitudes accessible by conventional titrimetric techniques. Additionally, because of the sensitivity of fluorescence, the method requires several orders of magnitude less material than most other solution methods. We discuss the advantages of this method for detecting and characterizing any modification that alters duplex stability, including, but not limited to, mutagenic lesions. We underscore the wide range of accessible free energy values that can be defined by this method, the applicability of the method in probing for a myriad of nucleic acid variations, such as single nucleotide polymorphisms, and the potential of the method for high throughput screening.

Absence of stable intermediates on the folding pathway of barnase

Barnase is one of the few protein models that has been studied extensively for protein folding. Previous studies led to the conclusion that barnase folds through a very stable submillisecond intermediate (≈3 kcal/mol). The structure of this intermediate was characterized intensively by using a protein engineering approach. This intermediate has now been reexamined with three direct and independent methods. (i) Hydrogen exchange experiments show very small protection factors (≈2) for the putative intermediate, indicating a stability of ≈0.0 kcal/mol. (ii) Denaturant-dependent unfolding of the putative intermediate is noncooperative and indicates a stability less than 0.0 kcal/mol. (iii) The logarithm of the unfolding rate constant of native barnase vs. denaturant concentrations is not linear. Together with the measured rate (“I” to N), this nonlinear behavior accounts for almost all of the protein stability, leaving only about 0.3 kcal/mol that could be attributed to the rapidly formed intermediate. Other observations previously interpreted to support the presence of an intermediate are now known to have alternative explanations. These results cast doubts on the previous conclusions on the nature of the early folding state in barnase and therefore should have important implications in understanding the early folding events of barnase and other proteins in general.

A meanfield approach to the thermodynamics of a protein-solvent system with application to the oligomerization of the tumor suppressor p53

The thermodynamic stability and oligomerization status of the tumor suppressor p53 tetramerization domain have been studied experimentally and theoretically. A series of hydrophilic mutations at Met-340 and Leu-344 of human p53 were designed to disrupt the hydrophobic dimer–dimer interface of the tetrameric oligomerization domain of p53 (residues 325–355). Meanfield calculations of the free energy of the solvated mutants as a function of interdimer distance were compared with experimental data on the thermal stability and oligomeric state (tetramer, dimer, or equilibrium mixture of both) of each mutant. The calculations predicted a decreasing stability and oligomeric state for the following amino acids at residue 340: Met (tetramer) > Ser Asp, His, Gln, > Glu, Lys (dimer), whereas the experimental results showed the following order: Met (tetramer) > Ser > Gln > His, Lys > Asp, Glu (dimers). For residue 344, the calculated trend was Leu (tetramer) > Ala > Arg, Gln, Lys (dimer), and the experimental trend was Leu (tetramer) > Ala, Arg, Gln, Lys (dimer). The discrepancy for the lysine side chain at residue 340 is attributed to the dual nature of lysine, both hydrophobic and charged. The incorrect prediction of stability of the mutant with Asp at residue 340 is attributed to the fact that within the meanfield approach, we use the wild-type backbone configuration for all mutants, but low melting temperatures suggest a softening of the α-helices at the dimer–dimer interface. Overall, this initial application of meanfield theory toward a protein-solvent system is encouraging for the application of the theoretical model to more complex systems.

Stimulation of bacteriophage T4 middle transcription by the T4 proteins MotA and AsiA occurs at two distinct steps in the transcription cycle

The bacteriophage T4 encodes proteins that are responsible for tightly regulating mRNA synthesis throughout phage development in Escherichia coli. The three classes of T4 promoters (early, middle, and late) are utilized sequentially by the host RNA polymerase as a result of phage-induced modifications. One such modification is the tight binding of the T4 AsiA protein to the σ70 subunit of the RNA polymerase. This interaction is pivotal for the transition between T4 early and middle transcription, since it both inhibits recognition of host and T4 early promoters and stimulates T4 middle mode synthesis. The activation of T4 middle transcription also requires the T4 MotA protein, bound specifically to its recognition sequence, the “Mot box,” which is centered at position −30 of these promoters. Accordingly, the two T4 proteins working in concert are sufficient to effectively switch the transcription specificity of the RNA polymerase holoenzyme. Herein, we investigate the mechanism of transcription activation and report that, while the presence of MotA and AsiA increases the initial recruitment of RNA polymerase to a T4 middle promoter, it does not alter the intrinsic stability of the discrete complexes formed. In addition, we have characterized the RNA polymerase-promoter species by UV laser footprinting and followed their evolution from open into initiating complexes. These data, combined with in vitro transcription assays, indicate that AsiA and MotA facilitate promoter escape, thereby stimulating the production of full-length transcripts.

Site-specifically located 8-amino-2′-deoxyguanosine: thermodynamic stability and mutagenic properties in Escherichia coli

2-Nitropropane (2-NP), an important industrial solvent and a component of cigarette smoke, is mutagenic in bacteria and carcinogenic in rats. 8-Amino-2′-deoxyguanosine (8-amino-dG) is one of the types of DNA damage found in liver, the target organ in 2-NP-treated rats. To investigate the thermodynamic properties of 8-amino-dG opposite each of the four DNA bases, we have synthesized an 11mer, d(CCATCG*CTACC), in which G* represents the modified base. By annealing a complementary DNA strand to this modified 11mer, four sets of duplexes were generated each containing one of the four DNA bases opposite the lesion. Circular dichroism studies indicated that 8-amino-dG did not alter the global helical properties of natural right-handed B-DNA. The thermal stability of each duplex was examined by UV melting measurements and compared with its unmodified counterpart. For the unmodified 11mer, the relative stability of the complementary DNA bases opposite G was in the order C > T > G > A, as determined from their –ΔG° values. The free energy change of each modified duplex was lower than its unmodified counterpart, except for the G*:G pair that exhibited a higher melting transition and a larger –ΔG° than the G:G duplex. Nevertheless, the stability of the modified 11mer duplex also followed the order C > T > G > A when placed opposite 8-amino-dG. To explore if 8-amino-dG opposite another 8-amino-dG has any advantage in base pairing, a G*:G* duplex was evaluated, which showed that the stability of this duplex was similar to the G*:G duplex. Mutagenesis of 8-amino-dG in this sequence context was studied in Escherichia coli, which showed that the lesion is weakly mutagenic (mutation frequency ∼10–3) but still can induce a variety of targeted and semi-targeted mutations.

The Interaction between Cold and Light Controls the Expression of the Cold-Regulated Barley Gene cor14b and the Accumulation of the Corresponding Protein1

We report the expression of the barley (Hordeum vulgare L.) COR (cold-regulated) gene cor14b (formerly pt59) and the accumulation of its chloroplast-localized protein product. A polyclonal antibody raised against the cor14b-encoded protein detected two chloroplast COR proteins: COR14a and COR14b. N-terminal sequencing of COR14a and expression of cor14b in Arabidopsis plants showed that COR14a is not encoded by the cor14b sequence, but it shared homology with the wheat (Triticum aestivum L.) WCS19 COR protein. The expression of cor14b was strongly impaired in the barley albino mutant an, suggesting the involvement of a plastidial factor in the control of gene expression. Low-level accumulation of COR14b was induced by cold treatment in etiolated plants, although cor14b expression and protein accumulation were enhanced after a short light pulse. Light quality was a determining factor in regulating gene expression: red or blue but not far-red or green light pulses were able to promote COR14b accumulation in etiolated plants, suggesting that phytochrome and blue light photoreceptors may be involved in the control of cor14b gene expression. Maximum accumulation of COR14b was reached only when plants were grown and/or hardened under the standard photoperiod. The effect of light on the COR14b stability was demonstrated by using transgenic Arabidopsis. These plants constitutively expressed cor14b mRNAs regardless of temperature and light conditions; nevertheless, green plants accumulated about twice as much COR14b protein as etiolated plants.