Biomedical subjects

Daniel Herschlag

Publications and source records attributed to Daniel Herschlag.

At least 19 recordsLinked to original sources

PubMed2006

Direct measurement of the full, sequence-dependent folding landscape of a nucleic acid.

Nucleic acid hairpins provide a powerful model system for understanding macromolecular folding, with free-energy landscapes that can be readily manipulated by changing the hairpin sequence. The full shapes of energy landscapes for the reversible folding of DNA hairpins under controlled loads exerted by an optical force clamp were obtained by deconvolution from high-resolution, single-molecule trajectories. The locations and heights of the energy barriers for hairpin folding could be tuned by adjusting the number and location of G:C base pairs, and the presence and position of folding intermediates were controlled by introducing single-nucleotide mismatches.

Base Pair Mismatch↗

PubMed2006

Structural transitions and thermodynamics of a glycine-dependent riboswitch from Vibrio cholerae.

Riboswitches are complex folded RNA domains found in noncoding regions of mRNA that regulate gene expression upon small molecule binding. Recently, Breaker and coworkers reported a tandem aptamer riboswitch (VCI-II) that binds glycine cooperatively. Here, we use hydroxyl radical footprinting and small-angle X-ray scattering (SAXS) to study the conformations of this tandem aptamer as a function of Mg(2+) and glycine concentration. We fit a simple three-state thermodynamic model that describes the energetic coupling between magnesium-induced folding and glycine binding. Furthermore, we characterize the structural conformations of each of the three states: In low salt with no magnesium present, the VCI-II construct has an extended overall conformation, presumably representing unfolded structures. Addition of millimolar concentrations of Mg(2+) in the absence of glycine leads to a significant compaction and partial folding as judged by hydroxyl radical protections. In the presence of millimolar Mg(2+) concentrations, the tandem aptamer binds glycine cooperatively. The glycine binding transition involves a further compaction, additional tertiary packing interactions and further uptake of magnesium ions relative to the state in high Mg(2+) but no glycine. Employing density reconstruction algorithms, we obtain low resolution 3-D structures for all three states from the SAXS measurements. These data provide a first glimpse into the structural conformations of the VCI-II aptamer, establish rigorous constraints for further modeling, and provide a framework for future mechanistic studies.

Aptamers, Nucleotide↗

PubMed2006

The paradoxical behavior of a highly structured misfolded intermediate in RNA folding.

Like many structured RNAs, the Tetrahymena group I ribozyme is prone to misfolding. Here we probe a long-lived misfolded species, referred to as M, and uncover paradoxical aspects of its structure and folding. Previous work indicated that a non-native local secondary structure, termed alt P3, led to formation of M during folding in vitro. Surprisingly, hydroxyl radical footprinting, fluorescence measurements with site-specifically incorporated 2-aminopurine, and functional assays indicate that the native P3, not alt P3, is present in the M state. The paradoxical behavior of alt P3 presumably arises because alt P3 biases folding toward M, but, after commitment to this folding pathway and before formation of M, alt P3 is replaced by P3. Further, structural and functional probes demonstrate that the misfolded ribozyme contains extensive native structure, with only local differences between the two states, and the misfolded structure even possesses partial catalytic activity. Despite the similarity of these structures, re-folding of M to the native state is very slow and is strongly accelerated by urea, Na+, and increased temperature and strongly impeded by Mg2+ and the presence of native peripheral contacts. The paradoxical observations of extensive native structure within the misfolded species but slow conversion of this species to the native state are readily reconciled by a model in which the misfolded state is a topological isomer of the native state, and computational results support the feasibility of this model. We speculate that the complex topology of RNA secondary structures and the inherent rigidity of RNA helices render kinetic traps due to topological isomers considerably more common for RNA than for proteins.

Animals↗

PubMed2006

Structural and functional comparisons of nucleotide pyrophosphatase/phosphodiesterase and alkaline phosphatase: implications for mechanism and evolution.

The rapid expansion of the amount of genomic and structural data has provided many examples of enzymes with evolutionarily related active sites that catalyze different reactions. Functional comparisons of these active sites can provide insight into the origins of the enormous catalytic proficiency of enzymes and the evolutionary changes that can lead to different enzyme activities. The alkaline phosphatase (AP) superfamily is an ideal system to use in making such comparisons given the extensive data available on both nonenzymatic and enzymatic phosphoryl transfer reactions. Some superfamily members, such as AP itself, preferentially hydrolyze phosphate monoesters, whereas others, such as nucleotide pyrophosphatase/phosphodiesterase (NPP), preferentially hydrolyze phosphate diesters. We have measured rate constants for NPP-catalyzed hydrolysis of phosphate diesters and monoesters. NPP preferentially catalyzes diester hydrolysis by factors of 10(2)-10(6), depending on the identity of the diester substrate. To identify features of the NPP active site that could lead to preferential phosphate diester hydrolysis, we have determined the structure of NPP in the absence of ligands and in complexes with vanadate and AMP. Comparisons to existing structures of AP reveal bimetallo cores that are structurally indistinguishable, but there are several distinct structural features outside of the conserved bimetallo site. The structural and functional data together suggest that some of these distinct functional groups provide specific substrate binding interactions, whereas others tune the properties of the bimetallo active site itself to discriminate between phosphate diester and monoester substrates.

Adenosine Monophosphate↗

PubMed2006

Interdomain communication in DNA topoisomerase II. DNA binding and enzyme activation.

Topoisomerase II catalyzes the ATP-dependent transport of a DNA segment (T-DNA) through a transient double strand break in another DNA segment (G-DNA). A fundamental mechanistic question is how the individual steps in this process are coordinated. We probed communication between the DNA binding sites and the individual enzymatic activities, ATP hydrolysis, and DNA cleavage. We employed short DNA duplexes to control occupancy at the two binding sites of wild-type enzyme and a variant with a G-DNA site mutation. The DNA concentration dependence of ATP hydrolysis and a fluorescence anisotropy assay provided thermodynamic information about DNA binding. The results suggest that G-DNA binds with higher affinity than T-DNA. Enzyme with only G-DNA bound is competent to cleave DNA, indicating that T-DNA is dispensable for DNA cleavage. The ATPase activity of enzyme bound solely to G-DNA is partially stimulated. Full stimulation requires binding of T-DNA. Both DNA binding sites therefore signal to the ATPase domains. The results support and extend current mechanistic models for topoisomerase II-catalyzed DNA transport and provide a framework for future mechanistic dissection.

Adenosine Triphosphatases↗

PubMed2006

Nanomechanical measurements of the sequence-dependent folding landscapes of single nucleic acid hairpins.

Nucleic acid hairpins provide a powerful model system for probing the formation of secondary structure. We report a systematic study of the kinetics and thermodynamics of the folding transition for individual DNA hairpins of varying stem length, loop length, and stem GC content. Folding was induced mechanically in a high-resolution optical trap using a unique force clamp arrangement with fast response times. We measured 20 different hairpin sequences with quasi-random stem sequences that were 6-30 bp long, polythymidine loops that were 3-30 nt long, and stem GC content that ranged from 0% to 100%. For all hairpins studied, folding and unfolding were characterized by a single transition. From the force dependence of these rates, we determined the position and height of the energy barrier, finding that the transition state for duplex formation involves the formation of 1-2 bp next to the loop. By measuring unfolding energies spanning one order of magnitude, transition rates covering six orders of magnitude, and hairpin opening distances with subnanometer precision, our results define the essential features of the energy landscape for folding. We find quantitative agreement over the entire range of measurements with a hybrid landscape model that combines thermodynamic nearest-neighbor free energies and nanomechanical DNA stretching energies.

Base Composition↗

PubMed2006

Testing electrostatic complementarity in enzyme catalysis: hydrogen bonding in the ketosteroid isomerase oxyanion hole.

A longstanding proposal in enzymology is that enzymes are electrostatically and geometrically complementary to the transition states of the reactions they catalyze and that this complementarity contributes to catalysis. Experimental evaluation of this contribution, however, has been difficult. We have systematically dissected the potential contribution to catalysis from electrostatic complementarity in ketosteroid isomerase. Phenolates, analogs of the transition state and reaction intermediate, bind and accept two hydrogen bonds in an active site oxyanion hole. The binding of substituted phenolates of constant molecular shape but increasing pK(a) models the charge accumulation in the oxyanion hole during the enzymatic reaction. As charge localization increases, the NMR chemical shifts of protons involved in oxyanion hole hydrogen bonds increase by 0.50-0.76 ppm/pK(a) unit, suggesting a bond shortening of 0.02 A/pK(a) unit. Nevertheless, there is little change in binding affinity across a series of substituted phenolates (DeltaDeltaG = -0.2 kcal/mol/pK(a) unit). The small effect of increased charge localization on affinity occurs despite the shortening of the hydrogen bonds and a large favorable change in binding enthalpy (DeltaDeltaH = -2.0 kcal/mol/pK(a) unit). This shallow dependence of binding affinity suggests that electrostatic complementarity in the oxyanion hole makes at most a modest contribution to catalysis of 300-fold. We propose that geometrical complementarity between the oxyanion hole hydrogen-bond donors and the transition state oxyanion provides a significant catalytic contribution, and suggest that KSI, like other enzymes, achieves its catalytic prowess through a combination of modest contributions from several mechanisms rather than from a single dominant contribution.

Anions↗

PubMed2006

Genome-wide identification of mRNAs associated with the translational regulator PUMILIO in Drosophila melanogaster.

Genome-wide identification of RNAs associated with RNA-binding proteins is crucial for deciphering posttranscriptional regulatory systems. PUMILIO is a member of the evolutionary conserved Puf-family of RNA-binding proteins that repress gene expression posttranscriptionally. We generated transgenic flies expressing affinity-tagged PUMILIO under the control of an ovary-specific promoter, and we purified PUMILIO from whole adult flies and embryos and analyzed associated mRNAs by using DNA microarrays. Distinct sets comprising hundreds of mRNAs were associated with PUMILIO at the two developmental stages. Many of these mRNAs encode functionally related proteins, supporting a model for coordinated regulation of posttranscriptional modules by specific RNA-binding proteins. We identified a characteristic sequence motif in the 3'-untranslated regions of mRNAs associated with PUMILIO, and the sufficiency of this motif for interaction with PUMILIO was confirmed by RNA pull-down experiments with biotinylated synthetic RNAs. The RNA motif strikingly resembles the one previously identified for Puf3p, one of five Saccharomyces cerevisiae Puf proteins; however, proteins encoded by the associated mRNAs in yeast and Drosophila do not appear to be related. The results suggest extensive posttranscriptional regulation by PUMILIO and uncover evolutionary features of this conserved family of RNA-binding proteins.

3' Untranslated Regions↗

PubMed2006

Alkaline phosphatase mono- and diesterase reactions: comparative transition state analysis.

Enzyme-catalyzed phosphoryl transfer reactions have frequently been suggested to proceed through transition states that are altered from their solution counterparts. Previous work with Escherichia coli alkaline phosphatase (AP), however, suggests that this enzyme catalyzes the hydrolysis of phosphate monoesters through a loose, dissociative transition state, similar to that in solution. AP also exhibits catalytic promiscuity, with a low level of phosphodiesterase activity, despite the tighter, more associative transition state for phosphate diester hydrolysis in solution. Because AP is evolutionarily optimized for phosphate monoester hydrolysis, it is possible that the active site environment alters the transition state for diester hydrolysis to be looser in its bonding to the incoming and outgoing groups. To test this possibility, we have measured the nonenzymatic and AP-catalyzed rate of reaction for a series of substituted methyl phenyl phosphate diesters. The values of beta(lg) and additional data suggest that the transition state for AP-catalyzed phosphate diester hydrolysis is indistinguishable from that in solution. Instead of altering transition state structure, AP catalyzes phosphoryl transfer reactions by recognizing and stabilizing transition states similar to those in aqueous solution. The AP active site therefore has the ability to recognize different transition states, a property that could assist in the evolutionary optimization of promiscuous activities.

Alkaline Phosphatase↗

PubMed2005

Functional identification of catalytic metal ion binding sites within RNA.

The viability of living systems depends inextricably on enzymes that catalyze phosphoryl transfer reactions. For many enzymes in this class, including several ribozymes, divalent metal ions serve as obligate cofactors. Understanding how metal ions mediate catalysis requires elucidation of metal ion interactions with both the enzyme and the substrate(s). In the Tetrahymena group I intron, previous work using atomic mutagenesis and quantitative analysis of metal ion rescue behavior identified three metal ions (MA, MB, and MC) that make five interactions with the ribozyme substrates in the reaction's transition state. Here, we combine substrate atomic mutagenesis with site-specific phosphorothioate substitutions in the ribozyme backbone to develop a powerful, general strategy for defining the ligands of catalytic metal ions within RNA. In applying this strategy to the Tetrahymena group I intron, we have identified the pro-SP phosphoryl oxygen at nucleotide C262 as a ribozyme ligand for MC. Our findings establish a direct connection between the ribozyme core and the functionally defined model of the chemical transition state, thereby extending the known set of transition-state interactions and providing information critical for the application of the recent group I intron crystallographic structures to the understanding of catalysis.

Animals↗

PubMed2005

Structural specificity conferred by a group I RNA peripheral element.

Like proteins, structured RNAs must specify a native conformation that is more stable than all other possible conformations. Local structure is much more stable for RNA than for protein, so it is likely that the principal challenge for RNA is to stabilize the native structure relative to misfolded and partially folded intermediates rather than unfolded structures. Many structured RNAs contain peripheral structural elements, which surround the core elements. Although it is clear that peripheral elements stabilize structure within RNAs that contain them, it has not yet been explored whether they specifically stabilize the native states relative to alternative folds. A two-piece version of the group I intron RNA from Tetrahymena is used here to show that the peripheral element P5abc binds to the native conformation of the rest of the RNA 50,000 times more tightly than it binds to a long-lived misfolded conformation. Thus, P5abc stabilizes the native conformation by approximately 6 kcal/mol relative to this misfolded conformation. Further, activity measurements show that for the RNA lacking P5abc, the native conformation is only marginally preferred over the misfolded conformation (<0.5 kcal/mol), indicating that the peripheral structure of this RNA is required to achieve a significant thermodynamic preference for the native state. Such "structural specificity" may be a general function of RNA peripheral domains.

Animals↗

PubMed2005

Alkaline phosphatase catalysis is ultrasensitive to charge sequestered between the active site zinc ions.

Escherichia coli alkaline phosphatase (AP) is a prototypical bimetalloenzyme, facilitating catalysis of phosphate monoester hydrolysis with two Zn2+ metal ions that are only 4 A apart. In the reaction's transition state, one of the nonbridging oxygen atoms of the transferred group appears to interact directly with the Zn2+ ion metallocluster. To determine the importance and the energetic properties of this interaction, we systematically varied the charge on this oxygen atom, exploiting the ability of AP to catalyze reactions of different classes of substrates. We observed that the AP catalytic proficiency correlates very well (R2 = 0.98) with the charge on this oxygen atom, over 8 orders of magnitude of catalytic proficiency. The slope of this linear correlation (31 +/- 2 kcal/mol per unit charge) is extraordinarily steep, indicating that AP greatly discriminates between differentially charged substrates. We suggest that this discrimination arises via an electrostatic interaction with the bimetallocluster. The dependence of the AP catalytic proficiency on the nonbridging oxygen charge is much larger than charge perturbation effects observed previously for other proteins. We propose that AP uses folding energy to position the two Zn2+ metal ions in close proximity, thereby creating an active site with a high electrostatic potential that is extraordinarily sensitive to the charge that "solvates" the metallocluster. The sensitivity of enzyme energetics to systematic variation in electrostatic properties provides a powerful measure of the active site environment. Future work comparing the sensitivity of related enzymes that have been optimized to catalyze different reactions will help reveal how natural selection has tuned related active sites to favor different reactions.

Alkaline Phosphatase↗

PubMed2005

Determining the Mg2+ stoichiometry for folding an RNA metal ion core.

The folding and catalytic function of RNA molecules depend on their interactions with divalent metal ions, such as magnesium. As with every molecular process, the most basic knowledge required for understanding the close relationship of an RNA with its metal ions is the stoichiometry of the interaction. Unfortunately, inventories of the numbers of divalent ions associated with unfolded and folded RNA states have been unattainable. A common approach has been to interpret Hill coefficients fit to folding equilibria as the number of metal ions bound upon folding. However, this approach is vitiated by the presence of diffusely associated divalent ions in a dynamic ion atmosphere and by the likelihood of multiple transitions along a folding pathway. We demonstrate that the use of molar concentrations of background monovalent salt can alleviate these complications. These simplifying solution conditions allow a precise determination of the stoichiometry of the magnesium ions involved in folding the metal ion core of the P4-P6 domain of the Tetrahymena group I ribozyme. Hill analysis of hydroxyl radical footprinting data suggests that the P4-P6 RNA core folds cooperatively upon the association of two metal ions. This unexpectedly small stoichiometry is strongly supported by counting magnesium ions associated with the P4-P6 RNA via fluorescence titration and atomic emission spectroscopy. By pinpointing the metal ion stoichiometry, these measurements provide a critical but previously missing step in the thermodynamic dissection of the coupling between metal ion binding and RNA folding.

Catalysis↗

PubMed2005

Promiscuous catalysis by the tetrahymena group I ribozyme.

Catalytic promiscuity, the ability of an enzyme to catalyze alternative reactions, has been suggested to have played an important role in the evolution of new catalytic activities in protein enzymes. Similarly, promiscuous activities may have been advantageous in an earlier RNA world. The Tetrahymena Group I ribozyme naturally catalyzes the site-specific guanosine attack on an anionic phosphate diester and has been shown to also catalyze aminoacyl transfer to water, albeit with a small rate acceleration (<10-fold). This inefficient catalysis could be due to the differences in charge and/or geometry requirements for the two reactions. Herein, we describe a new promiscuous activity of this ribozyme, the site-specific guanosine attack on a neutral phosphonate diester. This alternative substrate lacks the negative charge at the reaction center but, in contrast to the aminoacyl substrate, can undergo nucleophilic attack with the same geometry as the natural substrate. Our results show that the neutral phosphonate reaction is catalyzed about 1 x 106-fold, substantially better than the acyl transfer but far below the normal anionic substrate. We conclude that both charge and geometry are important factors for catalysis of the normal reaction and that promiscuous catalytic activities of ribozymes could have been created or enhanced by reorienting and swapping RNA domains.

Animals↗

PubMed2005

Dissecting eukaryotic translation and its control by ribosome density mapping.

Translation of an mRNA is generally divided into three stages: initiation, elongation and termination. The relative rates of these steps determine both the number and position of ribosomes along the mRNA, but traditional velocity sedimentation assays for the translational status of mRNA determine only the number of bound ribosomes. We developed a procedure, termed Ribosome Density Mapping (RDM), that uses site-specific cleavage of polysomal mRNA followed by separation on a sucrose gradient and northern analysis, to determine the number of ribosomes associated with specified portions of a particular mRNA. This procedure allows us to test models for translation and its control, and to examine properties of individual steps of translation in vivo. We tested specific predictions from the current model for translational control of GCN4 expression in yeast and found that ribosomes were differentially associated with the uORFs elements and coding region under different growth conditions, consistent with this model. We also mapped ribosome density along the ORF of several mRNAs, to probe basic kinetic properties of translational steps in yeast. We found no detectable decline in ribosome density between the 5' and 3' ends of the ORFs, suggesting that the average processivity of elongation is very high. Conversely, there was no queue of ribosomes at the termination site, suggesting that termination is not very slow relative to elongation and initiation. Finally, the RDM results suggest that less frequent initiation of translation on mRNAs with longer ORFs is responsible for the inverse correlation between ORF length and ribosomal density that we observed in a global analysis of translation. These results provide new insights into eukaryotic translation in vivo.

DNA-Binding Proteins↗

PubMed2005

Probing counterion modulated repulsion and attraction between nucleic acid duplexes in solution.

Understanding biological and physical processes involving nucleic acids, such as the binding of proteins to DNA and RNA, DNA condensation, and RNA folding, requires an understanding of the ion atmosphere that surrounds nucleic acids. We have used a simple model DNA system to determine how the ion atmosphere modulates interactions between duplexes in the absence of specific metal ion-binding sites and other complicated interactions. In particular, we have tested whether the Coulomb repulsion between nucleic acids can be reversed by counterions to give a net attraction, as has been proposed recently for the rapid collapse observed early in RNA folding. The conformation of two DNA duplexes tethered by a flexible neutral linker was determined in the presence of a series of cations by small angle x-ray scattering. The small angle x-ray scattering profiles of two control molecules with distinct shapes (a continuous duplex and a mimic of the compact DNA) were in good agreement with predictions, establishing the applicability of this approach. Under low-salt conditions (20 mM Na+), the tethered duplexes are extended because of a Coulombic repulsion estimated to be 2-5 kT/bp. Addition of high concentrations of Na+ (1.2 M), Mg2+ (0.6 M), and spermidine3+ (75 mM) resulted in electrostatic relaxation to a random state. These results indicate that a counterion-induced attractive force between nucleic acid duplexes is not significant under physiological conditions. An upper limit on the magnitude of the attractive potential under all tested ionic conditions is estimated.

DNA↗

PubMed2005

SAFA: semi-automated footprinting analysis software for high-throughput quantification of nucleic acid footprinting experiments.

Footprinting is a powerful and widely used tool for characterizing the structure, thermodynamics, and kinetics of nucleic acid folding and ligand binding reactions. However, quantitative analysis of the gel images produced by footprinting experiments is tedious and time-consuming, due to the absence of informatics tools specifically designed for footprinting analysis. We have developed SAFA, a semi-automated footprinting analysis software package that achieves accurate gel quantification while reducing the time to analyze a gel from several hours to 15 min or less. The increase in analysis speed is achieved through a graphical user interface that implements a novel methodology for lane and band assignment, called "gel rectification," and an optimized band deconvolution algorithm. The SAFA software yields results that are consistent with published methodologies and reduces the investigator-dependent variability compared to less automated methods. These software developments simplify the analysis procedure for a footprinting gel and can therefore facilitate the use of quantitative footprinting techniques in nucleic acid laboratories that otherwise might not have considered their use. Further, the increased throughput provided by SAFA may allow a more comprehensive understanding of molecular interactions. The software and documentation are freely available for download at http://safa.stanford.edu.

Algorithms↗

PubMed2004

Principles of RNA compaction: insights from the equilibrium folding pathway of the P4-P6 RNA domain in monovalent cations.

Counterions are required for RNA folding, and divalent metal ions such as Mg(2+) are often critical. To dissect the role of counterions, we have compared global and local folding of wild-type and mutant variants of P4-P6 RNA derived from the Tetrahymena group I ribozyme in monovalent and in divalent metal ions. A remarkably simple picture of the folding thermodynamics emerges. The equilibrium folding pathway in monovalent ions displays two phases. In the first phase, RNA molecules that are initially in an extended conformation enforced by charge-charge repulsion are relaxed by electrostatic screening to a state with increased flexibility but without formation of long-range tertiary contacts. At higher concentrations of monovalent ions, a state that is nearly identical to the native folded state in the presence of Mg(2+) is formed, with tertiary contacts that involve base and backbone interactions but without the subset of interactions that involve specific divalent metal ion-binding sites. The folding model derived from these and previous results provides a robust framework for understanding the equilibrium and kinetic folding of RNA.

Animals↗