Efficiency of DNA replication in the polymerase chain reaction

A detailed quantitative kinetic model for the polymerase chain reaction (PCR) is developed, which allows us to predict the probability of replication of a DNA molecule in terms of the physical parameters involved in the system. The important issue of the determination of the number of PCR cycles during which this probability can be considered to be a constant is solved within the framework of the model. New phenomena of multimodality and scaling behavior in the distribution of the number of molecules after a given number of PCR cycles are presented. The relevance of the model for quantitative PCR is discussed, and a novel quantitative PCR technique is proposed.

Limited internal friction in the rate-limiting step of a two-state protein folding reaction

Small, single-domain proteins typically fold via a compact transition-state ensemble in a process well fitted by a simple, two-state model. To characterize the rate-limiting conformational changes that underlie two-state folding, we have investigated experimentally the effects of changing solvent viscosity on the refolding of the IgG binding domain of protein L. In conjunction with numerical simulations, our results indicate that the rate-limiting conformational changes of the folding of this domain are strongly coupled to solvent viscosity and lack any significant “internal friction” arising from intrachain collisions. When compared with the previously determined solvent viscosity dependencies of other, more restricted conformational changes, our results suggest that the rate-limiting folding transition involves conformational fluctuations that displace considerable amounts of solvent. Reconciling evidence that the folding transition state ensemble is comprised of highly collapsed species with these and similar, previously reported results should provide a significant constraint for theoretical models of the folding process.

Cyclin A activates the DNA polymerase δ-dependent elongation machinery in vitro: A parvovirus DNA replication model

Replication of the single-stranded linear DNA genome of parvovirus minute virus of mice (MVM) starts with complementary strand synthesis from the 3′-terminal snap-back telomere, which serves as a primer for the formation of double-stranded replicative form (RF) DNA. This DNA elongation reaction, designated conversion, is exclusively dependent on cellular factors. In cell extracts, we found that complementary strand synthesis was inhibited by the cyclin-dependent kinase inhibitor p21WAF1/CIP1 and rescued by the addition of proliferating cell nuclear antigen, arguing for the involvement of DNA polymerase (Pol) δ in the conversion reaction. In vivo time course analyses using synchronized MVM-infected A9 cells allowed initial detection of MVM RF DNA at the G1/S phase transition, coinciding with the onset of cyclin A expression and cyclin A-associated kinase activity. Under in vitro conditions, formation of RF DNA was efficiently supported by A9 S cell extracts, but only marginally by G1 cell extracts. Addition of recombinant cyclin A stimulated DNA conversion in G1 cell extracts, and correlated with a concomitant increase in cyclin A-associated kinase activity. Conversely, a specific antibody neutralizing cyclin A-dependent kinase activity, abolished the capacity of S cell extracts for DNA conversion. We found no evidence for the involvement of cyclin E in the regulation of the conversion reaction. We conclude that cyclin A is necessary for activation of complementary strand synthesis, which we propose as a model reaction to study the cell cycle regulation of the Pol δ-dependent elongation machinery.

Sensitive detection of DNA polymorphisms by the serial invasive signal amplification reaction

The invasive signal amplification reaction has been previously developed for quantitative detection of nucleic acids and discrimination of single-nucleotide polymorphisms. Here we describe a method that couples two invasive reactions into a serial isothermal homogeneous assay using fluorescence resonance energy transfer detection. The serial version of the assay generates more than 107 reporter molecules for each molecule of target DNA in a 4-h reaction; this sensitivity, coupled with the exquisite specificity of the reaction, is sufficient for direct detection of less than 1,000 target molecules with no prior target amplification. Here we present a kinetic analysis of the parameters affecting signal and background generation in the serial invasive signal amplification reaction and describe a simple kinetic model of the assay. We demonstrate the ability of the assay to detect as few as 600 copies of the methylene tetrahydrofolate reductase gene in samples of human genomic DNA. We also demonstrate the ability of the assay to discriminate single base differences in this gene by using 20 ng of human genomic DNA.

How does a β-hairpin fold/unfold? Competition between topology and heterogeneity in a solvable model

We study the competition between topological effects and sequence inhomogeneities in determining the thermodynamics and the un/folding kinetics of a β-hairpin. Our work utilizes a new exactly solvable model that allows for arbitrary configurations of native contacts. In general, the competition between heterogeneity and topology results in a crossover of the dominant transition state. Interestingly, near this crossover, the single reaction coordinate picture can be seriously misleading. Our results also suggest that inferring the folding pathway from unfolding simulations is not always justified.

“Coarse” stability and bifurcation analysis using time-steppers: A reaction-diffusion example

Evolutionary, pattern forming partial differential equations (PDEs) are often derived as limiting descriptions of microscopic, kinetic theory-based models of molecular processes (e.g., reaction and diffusion). The PDE dynamic behavior can be probed through direct simulation (time integration) or, more systematically, through stability/bifurcation calculations; time-stepper-based approaches, like the Recursive Projection Method [Shroff, G. M. & Keller, H. B. (1993) SIAM J. Numer. Anal. 30, 1099–1120] provide an attractive framework for the latter. We demonstrate an adaptation of this approach that allows for a direct, effective (“coarse”) bifurcation analysis of microscopic, kinetic-based models; this is illustrated through a comparative study of the FitzHugh-Nagumo PDE and of a corresponding Lattice–Boltzmann model.

Relationship between the oxidation potential and electron spin density of the primary electron donor in reaction centers from Rhodobacter sphaeroides

The primary electron donor in bacterial reaction centers is a dimer of bacteriochlorophyll a molecules, labeled L or M based on their proximity to the symmetry-related protein subunits. The electronic structure of the bacteriochlorophyll dimer was probed by introducing small systematic variations in the bacteriochlorophyll–protein interactions by a series of site-directed mutations that replaced residue Leu M160 with histidine, tyrosine, glutamic acid, glutamine, aspartic acid, asparagine, lysine, and serine. The midpoint potentials for oxidation of the dimer in the mutants showed an almost continuous increase up to ≈60 mV compared with wild type. The spin density distribution of the unpaired electron in the cation radical state of the dimer was determined by electron–nuclear–nuclear triple resonance spectroscopy in solution. The ratio of the spin density on the L side of the dimer to the M side varied from ≈2:1 to ≈5:1 in the mutants compared with ≈2:1 for wild type. The correlation between the midpoint potential and spin density distribution was described using a simple molecular orbital model, in which the major effect of the mutations is assumed to be a change in the energy of the M half of the dimer, providing estimates for the coupling and energy levels of the orbitals in the dimer. These results demonstrate that the midpoint potential can be fine-tuned by electrostatic interactions with amino acids near the dimer and show that the properties of the electronic structure of a donor or acceptor in a protein complex can be directly related to functional properties such as the oxidation–reduction midpoint potential.

Structural model of cytochrome b559 in photosystem II based on a mutant with genetically fused subunits

Photosystem II is a reaction center protein complex located in photosynthetic membranes of plants, algae, and cyanobacteria. Using light energy, photosystem II catalyzes the oxidation of water and the reduction of plastoquinone, resulting in the release of molecular oxygen. A key component of photosystem II is cytochrome b559, a membrane-embedded heme protein with an unknown function. The cytochrome is unusual in that a heme links two separate polypeptide subunits, α and β, either as a heterodimer (αβ) or as two homodimers (α2 and β2). To determine the structural organization of cytochrome b559 in the membrane, we used site-directed mutagenesis to fuse the coding regions of the two respective genes in the cyanobacterium Synechocystis sp. PCC 6803. In this construction, the C terminus of the α subunit (9 kDa) is attached to the N terminus of the β subunit (5 kDa) to form a 14-kDa αβ fusion protein that is predicted to have two membrane-spanning α-helices with antiparallel orientations. Cells containing the αβ fusion protein grow photoautotrophically and assemble functional photosystem II complexes. Optical spectroscopy shows that the αβ fusion protein binds heme and is incorporated into photosystem II. These data support a structural model of cytochrome b559 in which one heme is coordinated to an α2 homodimer and a second heme is coordinated to a β2 homodimer. In this model, each photosystem II complex contains two cytochrome b559 hemes, with the α2 heme located near the stromal side of the membrane and the β2 heme located near the lumenal side.

Potency of Michael reaction acceptors as inducers of enzymes that protect against carcinogenesis depends on their reactivity with sulfhydryl groups

Induction of phase 2 enzymes and elevations of glutathione are major and sufficient strategies for protecting mammals and their cells against the toxic and carcinogenic effects of electrophiles and reactive forms of oxygen. Inducers belong to nine chemical classes and have few common properties except for their ability to modify sulfhydryl groups by oxidation, reduction, or alkylation. Much evidence suggests that the cellular “sensor” molecule that recognizes the inducers and signals the enhanced transcription of phase 2 genes does so by virtue of unique and highly reactive sulfhydryl functions that recognize and covalently react with the inducers. Benzylidene-alkanones and -cycloalkanones are Michael reaction acceptors whose inducer potency is profoundly increased by the presence of ortho- (but not other) hydroxyl substituent(s) on the aromatic ring(s). This enhancement correlates with more rapid reactivity of the ortho-hydroxylated derivatives with model sulfhydryl compounds. Proton NMR spectroscopy provides no evidence for increased electrophilicity of the β-vinyl carbons (the presumed site of nucleophilic attack) on the hydroxylated inducers. Surprisingly, these ortho-hydroxyl groups display a propensity for extensive intermolecular hydrogen bond formation, which may raise the reactivity and facilitate addition of mercaptans, thereby raising inducer potencies.

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.

Ras-catalyzed hydrolysis of GTP: a new perspective from model studies.

Despite the biological and medical importance of signal transduction via Ras proteins and despite considerable kinetic and structural studies of wild-type and mutant Ras proteins, the mechanism of Ras-catalyzed GTP hydrolysis remains controversial. We take a different approach to this problem: the uncatalyzed hydrolysis of GTP is analyzed, and the understanding derived is applied to the Ras-catalyzed reaction. Evaluation of previous mechanistic proposals from this chemical perspective suggests that proton abstraction from the attacking water by a general base and stabilization of charge development on the gamma-phosphoryl oxygen atoms would not be catalytic. Rather, this analysis focuses attention on the GDP leaving group, including the beta-gamma bridge oxygen of GTP, the atom that undergoes the largest change in charge in going from the ground state to the transition state. This leads to a new catalytic proposal in which a hydrogen bond from the backbone amide of Gly-13 to this bridge oxygen is strengthened in the transition state relative to the ground state, within an active site that provides a template complementary to the transition state. Strengthened transition state interactions of the active site lysine, Lys-16, with the beta-nonbridging phosphoryl oxygens and a network of interactions that positions the nucleophilic water molecule and gamma-phosphoryl group with respect to one another may also contribute to catalysis. It is speculated that a significant fraction of the GAP-activated GTPase activity of Ras arises from an additional interaction of the beta-gamma bridge oxygen with an Arg side chain that is provided in trans by GAP. The conclusions for Ras and related G proteins are expected to apply more widely to other enzymes that catalyze phosphoryl (-PO(3)2-) transfer, including kinases and phosphatases.

Simple model of protein folding kinetics.

A simple model of the kinetics of protein folding is presented. The reaction coordinate is the "correctness" of a configuration compared with the native state. The model has a gap in the energy spectrum, a large configurational entropy, a free energy barrier between folded and partially folded states, and a good thermodynamic folding transition. Folding kinetics is described by a master equation. The folding time is estimated by means of a local thermodynamic equilibrium assumption and then is calculated both numerically and analytically by solving the master equation. The folding time has a maximum near the folding transition temperature and can have a minimum at a lower temperature.

The folding of GroEL-bound barnase as a model for chaperonin-mediated protein folding.

We have analyzed the pathway of folding of barnase bound to GroEL to resolve the controversy of whether proteins can fold while bound to chaperonins (GroEL or Cpn60) or fold only after their release into solution. Four phases in the folding were detected by rapid-reaction kinetic measurements of the intrinsic fluorescence of both wild type and barnase mutants. The phases were assigned from their rate laws, sensitivity to mutations, and correspondence to regain of catalytic activity. At high ratios of denatured barnase to GroEL, 4 mol of barnase rapidly bind per 14-mer of GroEL. At high ratios of GroEL to barnase, 1 mol of barnase binds with a rate constant of 3.5 x 10(7) s-1.M-1. This molecule then refolds with a low rate constant that changes on mutation in parallel with the rate constant for the folding in solution. This rate constant corresponds to the regain of the overall catalytic activity of barnase and increases 15-fold on the addition of ATP to a physiologically relevant value of approximately 0.4 s-1. The multiply bound molecules of barnase that are present at high ratios of GroEL to barnase fold with a rate constant that is also sensitive to mutation but is 10 times higher. If the 110-residue barnase can fold when bound to GroEL and many moles can bind simultaneously, then smaller parts of large proteins should be able to fold while bound.

A multimer model for P680, the primary electron donor of photosystem II.

We consider a model of the photosystem II (PS II) reaction center in which its spectral properties result from weak (approximately 100 cm-1) excitonic interactions between the majority of reaction center chlorins. Such a model is consistent with a structure similar to that of the reaction center of purple bacteria but with a reduced coupling of the chlorophyll special pair. We find that this model is consistent with many experimental studies of PS II. The similarity in magnitude of the exciton coupling and energetic disorder in PS II results in the exciton states being structurally highly heterogeneous. This model suggests that P680, the primary electron donor of PS II, should not be considered a dimer but a multimer of several weakly coupled pigments, including the pheophytin electron acceptor. We thus conclude that even if the reaction center of PS II is structurally similar to that of purple bacteria, its spectroscopy and primary photochemistry may be very different.

A secondary-structure model for the self-cleaving region of Neurospora VS RNA.

Neurospora VS RNA performs an RNA-mediated self-cleavage reaction whose products contain 2',3'-cyclic phosphate and 5'-hydroxyl termini. This reaction is similar to those of hammerhead, hairpin, and hepatitis delta virus ribozymes; however, VS RNA is not similar in sequence to these other self-cleaving motifs. Here we propose a model for the secondary structure of the self-cleaving region of VS RNA, supported by site-directed mutagenesis and chemical modification structure probing data. The secondary structure of VS RNA is distinct from those of the other naturally occurring RNA self-cleaving domains. In addition to a unique secondary structure, several Mg-dependent interactions occur during the folding of VS RNA into its active tertiary conformation.