Balanced branching in transcription termination

The theory of stochastic transcription termination based on free-energy competition [von Hippel, P. H. & Yager, T. D. (1992) Science 255, 809–812 and von Hippel, P. H. & Yager, T. D. (1991) Proc. Natl. Acad. Sci. USA 88, 2307–2311] requires two or more reaction rates to be delicately balanced over a wide range of physical conditions. A large body of work on glasses and large molecules suggests that this balancing should be impossible in such a large system in the absence of a new organizing principle of matter. We review the experimental literature of termination and find no evidence for such a principle, but do find many troubling inconsistencies, most notably, anomalous memory effects. These effects suggest that termination has a deterministic component and may conceivably not be stochastic at all. We find that a key experiment by Wilson and von Hippel [Wilson, K. S. & von Hippel, P. H. (1994) J. Mol. Biol. 244, 36–51] thought to demonstrate stochastic termination was an incorrectly analyzed regulatory effect of Mg2+ binding.

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.

Quantification of the hydrophobic interaction by simulations of the aggregation of small hydrophobic solutes in water

The hydrophobic interaction, the tendency for nonpolar molecules to aggregate in solution, is a major driving force in biology. In a direct approach to the physical basis of the hydrophobic effect, nanosecond molecular dynamics simulations were performed on increasing numbers of hydrocarbon solute molecules in water-filled boxes of different sizes. The intermittent formation of solute clusters gives a free energy that is proportional to the loss in exposed molecular surface area with a constant of proportionality of 45 ± 6 cal/mol⋅Å2. The molecular surface area is the envelope of the solute cluster that is impenetrable by solvent and is somewhat smaller than the more traditional solvent-accessible surface area, which is the area transcribed by the radius of a solvent molecule rolled over the surface of the cluster. When we apply a factor relating molecular surface area to solvent-accessible surface area, we obtain 24 cal/mol⋅Å2. Ours is the first direct calculation, to our knowledge, of the hydrophobic interaction from molecular dynamics simulations; the excellent qualitative and quantitative agreement with experiment proves that simple van der Waals interactions and atomic point-charge electrostatics account for the most important driving force in biology.

Biochemical evolution III: Polymerization on organophilic silica-rich surfaces, crystal–chemical modeling, formation of first cells, and geological clues

Catalysis at organophilic silica-rich surfaces of zeolites and feldspars might generate replicating biopolymers from simple chemicals supplied by meteorites, volcanic gases, and other geological sources. Crystal–chemical modeling yielded packings for amino acids neatly encapsulated in 10-ring channels of the molecular sieve silicalite-ZSM-5-(mutinaite). Calculation of binding and activation energies for catalytic assembly into polymers is progressing for a chemical composition with one catalytic Al–OH site per 25 neutral Si tetrahedral sites. Internal channel intersections and external terminations provide special stereochemical features suitable for complex organic species. Polymer migration along nano/micrometer channels of ancient weathered feldspars, plus exploitation of phosphorus and various transition metals in entrapped apatite and other microminerals, might have generated complexes of replicating catalytic biomolecules, leading to primitive cellular organisms. The first cell wall might have been an internal mineral surface, from which the cell developed a protective biological cap emerging into a nutrient-rich “soup.” Ultimately, the biological cap might have expanded into a complete cell wall, allowing mobility and colonization of energy-rich challenging environments. Electron microscopy of honeycomb channels inside weathered feldspars of the Shap granite (northwest England) has revealed modern bacteria, perhaps indicative of Archean ones. All known early rocks were metamorphosed too highly during geologic time to permit simple survival of large-pore zeolites, honeycombed feldspar, and encapsulated species. Possible microscopic clues to the proposed mineral adsorbents/catalysts are discussed for planning of systematic study of black cherts from weakly metamorphosed Archaean sediments.

Monitoring the GAP catalyzed H-Ras GTPase reaction at atomic resolution in real time

The molecular reaction mechanism of the GTPase-activating protein (GAP)-catalyzed GTP hydrolysis by Ras was investigated by time resolved Fourier transform infrared (FTIR) difference spectroscopy using caged GTP (P3-1-(2-nitro)phenylethyl guanosine 5′-O-triphosphate) as photolabile trigger. This approach provides the complete GTPase reaction pathway with time resolution of milliseconds at the atomic level. Up to now, one structural model of the GAP⋅Ras⋅GDP⋅AlFx transition state analog is known, which represents a “snap shot” along the reaction-pathway. As now revealed, binding of GAP to Ras⋅GTP shifts negative charge from the γ to β phosphate. Such a shift was already identified by FTIR in GTP because of Ras binding and is now shown to be enhanced by GAP binding. Because the charge distribution of the GAP⋅Ras⋅GTP complex thus resembles a more dissociative-like transition state and is more like that in GDP, the activation free energy is reduced. An intermediate is observed on the reaction pathway that appears when the bond between β and γ phosphate is cleaved. In the intermediate, the released Pi is strongly bound to the protein and surprisingly shows bands typical of those seen for phosphorylated enzyme intermediates. All these results provide a mechanistic picture that is different from the intrinsic GTPase reaction of Ras. FTIR analysis reveals the release of Pi from the protein complex as the rate-limiting step for the GAP-catalyzed reaction. The approach presented allows the study not only of single proteins but of protein–protein interactions without intrinsic chromophores, in the non-crystalline state, in real time at the atomic level.

Configurational diffusion down a folding funnel describes the dynamics of DNA hairpins

Elucidating the mechanism of folding of polynucleotides depends on accurate estimates of free energy surfaces and a quantitative description of the kinetics of structure formation. Here, the kinetics of hairpin formation in single-stranded DNA are measured after a laser temperature jump. The kinetics are modeled as configurational diffusion on a free energy surface obtained from a statistical mechanical description of equilibrium melting profiles. The effective diffusion coefficient is found to be strongly temperature-dependent in the nucleation step as a result of formation of misfolded loops that do not lead to subsequent zipping. This simple system exhibits many of the features predicted from theoretical studies of protein folding, including a funnel-like energy surface with many folding pathways, trapping in misfolded conformations, and non-Arrhenius folding rates.

Binding of the cellulose-binding domain of exoglucanase Cex from Cellulomonas fimi to insoluble microcrystalline cellulose is entropically driven.

Isothermal titration microcalorimetry is combined with solution-depletion isotherm data to analyze the thermodynamics of binding of the cellulose-binding domain (CBD) from the beta-1,4-(exo)glucanase Cex of Cellulomonas fimi to insoluble bacterial microcrystalline cellulose. Analysis of isothermal titration microcalorimetry data against two putative binding models indicates that the bacterial microcrystalline cellulose surface presents two independent classes of binding sites, with the predominant high-affinity site being characterized by a Langmuir-type Ka of 6.3 (+/-1.4) x 10(7) M-1 and the low-affinity site by a Ka of 1.1 (+/-0.6) x 10(6) M-1. CBDCex binding to either site is exothermic, but is mainly driven by a large positive change in entropy. This differs from protein binding to soluble carbohydrates, which is usually driven by a relatively large exothermic standard enthalpy change for binding. Differential heat capacity changes are large and negative, indicating that sorbent and protein dehydration effects make a dominant contribution to the driving force for binding.

Mutations in copper-zinc superoxide dismutase that cause amyotrophic lateral sclerosis alter the zinc binding site and the redox behavior of the protein.

A series of mutant human and yeast copper-zinc superoxide dismutases has been prepared, with mutations corresponding to those found in familial amyotrophic lateral sclerosis (ALS; also known as Lou Gehrig's disease). These proteins have been characterized with respect to their metal-binding characteristics and their redox reactivities. Replacement of Zn2+ ion in the zinc sites of several of these proteins with either Cu2+ or Co2+ gave metal-substituted derivatives with spectroscopic properties different from those of the analogous derivative of the wild-type proteins, indicating that the geometries of binding of these metal ions to the zinc site were affected by the mutations. Several of the ALS-associated mutant copper-zinc superoxide dismutases were also found to be reduced by ascorbate at significantly greater rate than the wild-type proteins. We conclude that similar alterations in the properties of the zinc binding site can be caused by mutations scattered throughout the protein structure. This finding may help to explain what is perhaps the most perplexing question in copper-zinc superoxide dismutase-associated familial ALS-i.e., how such a diverse set of mutations can result in the same gain of function that causes the disease.

On the mechanism of action of ribonuclease A: relevance of enzymatic studies with a p-nitrophenylphosphate ester and a thiophosphate ester.

It has been reported that His-119 of ribonuclease A plays a major role as an imidazolium ion acid catalyst in the cyclization/cleavage of normal dinucleotides but that it is not needed for the cyclization/cleavage of 3'-uridyl p-nitrophenyl phosphate. We see that this is also true for simple buffer catalysis, where imidazole (as in His-12 of the enzyme), but not imidazolium ion, plays a significant catalytic role with the nitrophenyl substrate, but both are catalytic for normal dinucleotides such as uridyluridine. Rate studies show that the enzyme catalyzes the cyclization of the nitrophenylphosphate derivative 47,000,000 times less effectively (kcat/kuncat) than it does uridyladenosine, indicating that approximately 50% of the catalytic free energy change is lost with this substrate. This suggests that the nitrophenyl substrate is not correctly bound to take full advantage of the catalytic groups of the enzyme and is thus not a good guide to the mechanism used by normal nucleotides. The published data on kinetic effects with ribonuclease A of substituting thiophosphate groups for the phosphate groups of normal substrates has been discussed elsewhere, and it was argued that these effects are suggestive of the classical mechanism for ribonuclease action, not the novel mechanism we have recently proposed. The details of these rate effects, including stereochemical preferences in the thiophosphate series, can be invoked as support for our newer mechanism.

The mode of action and the structure of a herbicide in complex with its target: binding of activated hydantocidin to the feedback regulation site of adenylosuccinate synthetase.

(+)-Hydantocidin, a recently discovered natural spironucleoside with potent herbicidal activity, is shown to be a proherbicide that, after phosphorylation at the 5' position, inhibits adenylosuccinate synthetase, an enzyme involved in de novo purine synthesis. The mode of binding of hydantocidin 5'-monophosphate to the target enzyme was analyzed by determining the crystal structure of the enzyme-inhibitor complex at 2.6-A resolution. It was found that adenylosuccinate synthetase binds the phosphorylated compound in the same fashion as it does adenosine 5'-monophosphate, the natural feedback regulator of this enzyme. This work provides the first crystal structure of a herbicide-target complex reported to date.

Persistent confusion of total entropy and chemical system entropy in chemical thermodynamics.

The change in free energy with temperature at constant pressure of a chemical reaction is determined by the sum (dS) of changes in entropy of the system of reagents, dS(i), and the additional entropy change of the surroundings, dS(H), that results from the enthalpy change, W. A faulty identification of the total entropy change on reaction with dS(i) has been responsible for the attribution of general validity to the expressions (d deltaG/dT)p = -deltaS(i) and d(deltaG/T)/d(1/T)= deltaH, which are found in most textbooks and in innumerable papers.

Critical contribution of beta chain residue 57 in peptide binding ability of both HLA-DR and -DQ molecules.

Position 57 in the beta chain of HLA class II molecules maintains an Asp/non-Asp dimorphism that has been conserved through evolution and is implicated in susceptibility to some autoimmune diseases. The latter effect may be due to the influence of this residue on the ability of class II alleles to bind specific pathogenic peptides. We utilized highly homologous pairs of both DR and DQ alleles that varied at residue 57 to investigate the impact of this dimorphism on binding of model peptides. Using a direct binding assay of biotinylated peptides on whole cells expressing the desired alleles, we report several peptides that bind differentially to the allele pairs depending on the presence or absence of Asp at position 57. Peptides with negatively charged residues at anchor position 9 bind well to alleles not containing Asp at position 57 in the beta chain but cannot bind well to homologous Asp-positive alleles. By changing the peptides at the single residue predicted to interact with this position 57, we demonstrate a drastically altered or reversed pattern of binding. Ala analog peptides confirm these interactions and identify a limited set of interaction sites between the bound peptides and the class II molecules. Clarification of the impact of specific class II polymorphisms on generating unique allele-specific peptide binding "repertoires" will aid in our understanding of the development of specific immune responses and HLA-associated diseases.

Chaperonin-facilitated protein folding: optimization of rate and yield by an iterative annealing mechanism.

We develop a heuristic model for chaperonin-facilitated protein folding, the iterative annealing mechanism, based on theoretical descriptions of "rugged" conformational free energy landscapes for protein folding, and on experimental evidence that (i) folding proceeds by a nucleation mechanism whereby correct and incorrect nucleation lead to fast and slow folding kinetics, respectively, and (ii) chaperonins optimize the rate and yield of protein folding by an active ATP-dependent process. The chaperonins GroEL and GroES catalyze the folding of ribulose bisphosphate carboxylase at a rate proportional to the GroEL concentration. Kinetically trapped folding-incompetent conformers of ribulose bisphosphate carboxylase are converted to the native state in a reaction involving multiple rounds of quantized ATP hydrolysis by GroEL. We propose that chaperonins optimize protein folding by an iterative annealing mechanism; they repeatedly bind kinetically trapped conformers, randomly disrupt their structure, and release them in less folded states, allowing substrate proteins multiple opportunities to find pathways leading to the most thermodynamically stable state. By this mechanism, chaperonins greatly expand the range of environmental conditions in which folding to the native state is possible. We suggest that the development of this device for optimizing protein folding was an early and significant evolutionary event.

Dressed polyions, counterion condensation, and adsorption excess in polyelectrolyte solutions.

The phenomenon of Manning-Oosawa counterion condensation is given an explicit statistical mechanical and qualitative basis via a dressed polyelectrolyte formalism in connection with the topology of the electrostatic free-energy surface and is derived explicitly in terms of the adsorption excess of ions about the polyion via the nonlinear Poisson-Boltzmann equation. The approach is closely analogous to the theory of ion binding in micelles. Our results not only elucidate a Poisson-Boltzmann analysis, which shows that a fraction of the counterions lie within a finite volume around the polyion even if the volume of the system tends towards infinity, but also provide a direct link between Manning's theta-the number of condensed counterions for each polyion site-and a statistical thermodynamic quantity, namely, the adsorption excess per monomer.

Progress toward chemical accuracy in the computer simulation of condensed phase reactions.

We describe a procedure for the generation of chemically accurate computer-simulation models to study chemical reactions in the condensed phase. The process involves (i) the use of a coupled semiempirical quantum and classical molecular mechanics method to represent solutes and solvent, respectively; (ii) the optimization of semiempirical quantum mechanics (QM) parameters to produce a computationally efficient and chemically accurate QM model; (iii) the calibration of a quantum/classical microsolvation model using ab initio quantum theory; and (iv) the use of statistical mechanical principles and methods to simulate, on massively parallel computers, the thermodynamic properties of chemical reactions in aqueous solution. The utility of this process is demonstrated by the calculation of the enthalpy of reaction in vacuum and free energy change in aqueous solution for a proton transfer involving methanol, methoxide, imidazole, and imidazolium, which are functional groups involved with proton transfers in many biochemical systems. An optimized semiempirical QM model is produced, which results in the calculation of heats of formation of the above chemical species to within 1.0 kcal/mol (1 kcal = 4.18 kJ) of experimental values. The use of the calibrated QM and microsolvation QM/MM (molecular mechanics) models for the simulation of a proton transfer in aqueous solution gives a calculated free energy that is within 1.0 kcal/mol (12.2 calculated vs. 12.8 experimental) of a value estimated from experimental pKa values of the reacting species.