Pyrophosphorolysis-activated polymerization (PAP): application to allele-specific amplification

To measure mutation load or to detect minimal residual disease, a robust method for identifying one mutant allele in the range of 10(6)-10(9) wild-type alleles would be advantageous. Herein, we present evidence that pyrophosphorolysis-activated polymerization (PAP) has the potential to provide a highly specific and robust method of allele-specific amplification if DNA polymerases with higher pyrophosphorolysis activity can be found or engineered. In PAP, pyrophosphorolysis and polymerization by DNA polymerase are coupled serially by utilizing a pyrophosphorolysis-activatable oligonucleotide (P*). P*, which is an allele-specific oligonucleotide with a dideoxynucleotide at the 3' terminus, can be activated by pyrophosphorolysis to remove the 3' terminal dideoxynucleotide in the presence of pyrophosphate (PPi) and the complementary strand of the allelic template; then the activated P* can be extended by DNA polymerization. Specificity results from both pyrophosphorolysis and polymerization because significant nonspecific amplification requires the combination of mismatch pyrophosphorolysis and misincorporation by the DNA polymerase, which is an extremely rare event. Proof of principle has been achieved with a polymorphic site within the human D1 dopamine receptor gene. The effects of the dideoxyoligonucleotide sequences, DNA polymerases, PPi concentrations, allele-specific templates, pH and dNTP concentrations were examined.     

Pyrophosphorolysis-activatable oligonucleotides may facilitate detection of rare alleles, mutation scanning and analysis of chromatin structures

Pyrophosphorolysis-activated polymerization (PAP) was initially developed to enhance the specificity of allele-specific PCR for detection of known mutations in the presence of a great excess of wild-type allele. The high specificity of PAP derives from the serial coupling of pyrophosphorolysis-mediated activation of a pyrophosphorolysis-activatable oligonucleotide (P*) followed by extension of the activated oligonucleotide. Herein, we demonstrate that genetically engineered DNA polymerases greatly improve the efficiency of PAP, making it a practical technique for detection of rare mutations. We also show that P* oligonucleotides have the novel and unexpected property of high sensitivity to mismatches throughout at least the 16 3′-terminal nucleotides. Thus, PAP constitutes a technology platform of potential utility whenever high specificity is required along the length of an oligonucleotide.     

P1 nuclease cleavage is dependent on length of the mismatches in DNA

P1 nuclease is one of the most extensively used single-strand DNA specific nucleases in molecular biology. In modern biology, it is used as an enzymatic probe to detect altered DNA conformations. It is well documented that P1 cleaves single-stranded nucleic acids and single-stranded DNA regions. The fact that P1 can act under a wide range of conditions, including physiological pH and temperature make it the most commonly used enzymatic probe in DNA structural studies. Surprisingly, to this date, there is no study to characterize the influence of length of mismatches on P1 sensitivity. Using a series of radioactively labeled oligomeric DNA substrates-containing mismatches, we find that P1 nuclease cleavage is dependent on the length of mismatches. P1 does not cleave DNA when there is a single-base mismatch. P1 cleavage efficiency is optimum when mismatch length is 3 nt or more. Changing the position of the mismatches also does not make any difference in cleavage efficacy. These novel findings on P1 properties have implications for its use in DNA structure and DNA repair studies.     

Detection of extremely rare alleles by bidirectional pyrophosphorolysis-activated polymerization allele-specific amplification (Bi-PAP-A): measurement of mutation load in mammalian tissues

Pyrophosphorolysis-activated polymerization (PAP) was developed to detect extremely rare mutations in complex genomes. In theory, PAP can detect a copy of a single base mutation present in 3 x 10(11) copies of the wild-type allele. In practice, the selectivity of detection is limited by a bypass reaction involving a polymerase extension error from the unblocked oligonucleotide annealed to the opposing strand. Bidirectional PAP allele-specific amplification (Bi-PAP-A) is a novel method that uses two opposing 3'-terminal blocked pyrophosphorolysis-activatable oligonucleotides (P*s) with one nucleotide overlap at their 3' termini. This eliminates the problematic bypass reaction. The selectivity of Bi-PAP-A was examined using lambda phage DNA as a model system. Bi-PAP-A selectively detected two copies of a rare mutated allele in the presence of at least 2 x 10(9) copies of the wild-type lambda phage DNA. Bi-PAP-A was then applied to direct detection of spontaneous somatic mutations in the mouse genome at a frequency as low as 3 x 10(-9). A 370-fold variation in the frequency of a specific somatic mutation among different mouse samples was found, suggesting clonal expansion of mutation occurring during early development and a hyper-Poisson variance. Bi-PAP-A is a rapid, general, and automatable method for the detection of rare mutations.     

Multiplex dosage pyrophosphorolysis-activated polymerization: application to the detection of heterozygous deletions

Large heterozygous chromosomal deletions and gene duplications are important classes of mutations that are generally missed by standard PCR amplification and sequencing. Multiplex dosage pyrophosphorolysis-activated polymerization (MD-PAP), a derivative of PAP, was utilized to detect these types of mutations. PAP is a method for nucleic acid amplification in which 3' blocked oligonucleotides (P*) are activated by pyrophosphorolysis when annealed to the target template and subsequently extended. A key advantage to this technology is that PAP reactions produce little or no primer-dimer or false priming. As a result of this enhanced specificity, MD-PAP is easy to optimize. Herein, we utilize MD-PAP to determine gene dosage of each exon of the human factor IX gene by comparison with one endogenous internal control from the ATM gene. Estimated dosage is proportional to the actual template copy number over a minimum dynamic range from 1 to 16 copies. A blinded analysis detected 100% of 43 heterozygous deletions of exons in the human factor IX gene.     

Conformational determinants of tandem GU mismatches in RNA: insights from molecular dynamics simulations and quantum mechanical calculations

Structure and energetic properties of base pair mismatches in duplex RNA have been the focus of numerous investigations due to their role in many important biological functions. Such efforts have contributed to the development of models for secondary structure prediction of RNA, including the nearest-neighbor model. In RNA duplexes containing GU mismatches, 5'-GU-3' tandem mismatches have a different thermodynamic stability than 5'-UG-3' mismatches. In addition, 5'-GU-3' mismatches in some sequence contexts do not follow the nearest-neighbor model for stability. To characterize the underlying atomic forces that determine the structural and thermodynamic properties of GU tandem mismatches, molecular dynamics (MD) simulations were performed on a series of 5'-GU-3' and 5'-UG-3' duplexes in different sequence contexts. Overall, the MD-derived structural models agree well with experimental data, including local deviations in base step helicoidal parameters in the region of the GU mismatches and the model where duplex stability is associated with the pattern of GU hydrogen bonding. Further analysis of the simulations, validated by data from quantum mechanical calculations, suggests that the experimentally observed differences in thermodynamic stability are dominated by GG interstrand followed by GU intrastrand base stacking interactions that dictate the one versus two hydrogen bonding scenarios for the GU pairs. In addition, the inability of 5'-GU-3' mismatches in different sequence contexts to all fit into the nearest-neighbor model is indicated to be associated with interactions of the central four base pairs with the surrounding base pairs. The results emphasize the role of GG and GU stacking interactions on the structure and thermodynamics of GU mismatches in RNA.     

Characterization of polymerase chain reaction amplification of specific alleles

Under certain conditions, polymerase chain reaction (PCR) can be used to differentially amplify one allele over another. To characterize the phenomenon, we have made a series of PCR primers and determined whether differential amplification could be detected after agarose gel electrophoresis. Two allele pairs were examined; one pair represents a transversion and one pair represents a transition. The following conclusions emerge: (i) when the 3' or the 3' penultimate base of the oligonucleotide mismatched an allele, no amplification product could be detected; (ii) when the mismatches were 3 and 4 bases from the 3' end of the primer, differential amplification was still observed, but only at certain concentrations of magnesium chloride; (iii) the mismatched allele can be detected in the presence of a 40-fold excess of the matched allele; (iv) primers as short as 13 nucleotides were effective; and (v) the specificity of the amplification could be overwhelmed by greatly increasing the concentration of target DNA.     

Generic expansion of the substrate spectrum of a DNA polymerase by directed evolution

DNA polymerases recognize their substrates with exceptionally high specificity, restricting the use of unnatural nucleotides and the applications they enable. We describe a strategy to expand the substrate range of polymerases. By selecting for the extension of distorting 3' mismatches, we obtained mutants of Taq DNA polymerase that not only promiscuously extended mismatches, but had acquired a generic ability to process a diverse range of noncanonical substrates while maintaining high catalytic turnover, processivity and fidelity. Unlike the wild-type enzyme, they bypassed blocking lesions such as an abasic site, a thymidine dimer or the base analog 5-nitroindol and performed PCR amplification with complete substitution of all four nucleotide triphosphates with phosphorothioates or the substitution of one with the equivalent fluorescent dye-labeled nucleotide triphosphate. Such 'unfussy' polymerases have immediate utility, as we demonstrate by the generation of microarray probes with up to 20-fold brighter fluorescence.     

Star-PAP control of BIK expression and apoptosis is regulated by nuclear PIPKIα and PKCδ signaling

BIK protein is an initiator of mitochondrial apoptosis, and BIK expression is induced by proapoptotic signals, including DNA damage. Here, we demonstrate that 3' end processing and expression of BIK mRNA are controlled by the nuclear PI4,5P(2)-regulated poly(A) polymerase Star-PAP downstream of DNA damage. Nuclear PKCδ is a key mediator of apoptosis, and DNA damage stimulates PKCδ association with the Star-PAP complex where PKCδ is required for Star-PAP-dependent BIK expression. PKCδ binds the PI4,5P(2)-generating enzyme PIPKIα, which is essential for PKCδ interaction with the Star-PAP complex, and PKCδ activity is directly stimulated by PI4,5P(2). Features in the BIK 3' UTR uniquely define Star-PAP specificity and may block canonical PAP activity toward BIK mRNA. This reveals a nuclear phosphoinositide signaling nexus where PIPKIα, PI4,5P(2), and PKCδ regulate Star-PAP control of BIK expression and induction of apoptosis. This pathway is distinct from the Star-PAP-mediated oxidative stress pathway indicating signal-specific regulation of mRNA 3' end processing.     

Specificity and enzymatic mechanism of the editing exonuclease of Escherichia coli DNA polymerase III.

Exonucleolytic editing is a major contributor to the fidelity of DNA replication by the multisubunit DNA polymerase (pol) III holoenzyme. To investigate the source of editing specificity, we have studied the isolated exonuclease subunit, epsilon, and the pol III core subassembly, which carries the epsilon, theta, and alpha (polymerase) subunits. Using oligonucleotides with specific terminal mismatches, we have found that both epsilon and pol III core preferentially excise a mispaired 3' terminus and therefore have intrinsic editing specificity. For both epsilon and pol III core, exonuclease activity is much more effective with single-strand DNA; with a double-strand DNA, the exonuclease is strongly temperature-dependent. We conclude that the epsilon subunit of pol III holoenzyme is itself a specific editing exonuclease and that the source of specificity is the greater melting capacity of a mispaired 3' terminus.     

Sequence specificity of the P6 pairing for splicing of the group I td intron of phage T4.

Seventeen non-directed td- (thymidylate synthase-deficient) splicing-defective mutations isolated in phage T4 were localized within the catalytic core of the ribozyme. All of the mutations occur in conserved structural elements that form part of the td intron core secondary structure. Remarkably, seven of the seventeen independently isolated mutations clustered in the dinucleotide 5' element (P6[5']) of the putative two-base-pair P6 stem. An analysis of this region was undertaken by site-directed mutagenesis of the plasmid-borne td gene, leading to the following findings: First, the short P6 pairing in the td secondary structure model was verified with appropriate P6[5'] and P6[3'] compensatory mutations. Second, all P6[5'] and P6[3'] mutants are defective in the first step of splicing, guanosine-dependent 5' splice site cleavage, whereas their activity at the 3' splice site is variable. Third, residual in vitro splicing activity of the mutants altered on only one side of the P6 pairing is correlated with the ability to form an alternative two-base-pair P6 stem. Fourth, the degree to which the compensatory mutants have their splicing activity restored is highly condition-dependent. Restoration of phenotype of the compensatory P6[5']:[3'] constructs is weak under stringent in vitro conditions as well as in vivo. This sequence specificity is consistent with phylogenetic conservation of the P6 pairing elements in group I introns, and suggests either structural constraints on the P6 stem or a dual function of one or both pairing elements.     

Subsite preferences of pepstatin-insensitive carboxyl proteinases from prokaryotes: kumamolysin, a thermostable pepstatin-insensitive carboxyl proteinase

Kumamolysin, a carboxyl proteinase from Bacillus novosp. MN-32, is characterized by its thermostability and insensitivity to aspartic proteinase inhibitors such as pepstatin, diazoacetyl-DL-norleucine methylester, and 1,2-epoxy-3-(p-nitro-phenoxy)propane. Here, its substrate specificity was elucidated using two series of synthetic chromogenic substrates: P(5)-P(4)-P(3)-P(2)-Phe*Nph (p-nitrophenylalanine: *cleavage site)-P(2)'-P(3)', in which the amino acid residues at the P(5)-P(2), P(2)' and P(3)' positions were systematically substituted. Among 74 substrates, kumamolysin was shown to hydrolyze Lys-Pro-Ile-Pro-Phe-Nph-Arg-Leu most effectively. The kinetic parameters of this peptide were K(m) = 41+/-5 microM, k(cat) = 176+/- 10 s(-1), and k(cat)/K(m) = 4.3+/-0.6 mM(-1) x s(-1). These systematic analyses revealed the following features: (i) Kumamolysin had a unique preference for the P(2) position. Kumamolysin preferentially hydrolyzed peptides having an Ala or Pro residue at the P(2) position; this was also observed for the pepstatin-insensitive carboxyl proteinase from Bacillus coagulans J-4 [J-4; Shibata et al. (1998) J. Biochem. 124, 642-647]. Other carboxyl proteinases, including Pseudomonas sp. 101 pepstatin-insensitive carboxyl proteinase (PCP) and Xanthomonas sp. T-22 pepstatin-insensitive carboxyl proteinase (XCP), preferred peptides having hydrophobic and bulky amino acid residue such as Leu at the P(2) position. (ii) Kumamolysin preferred such charged amino acid residues as Glu or Arg at the P(2)' position, suggesting that the S(2)' subsite of kumamolysin is occupied by hydrophilic residues, similar to that of PCP, XCP, and J-4. In general, the S(2)' subsite of pepstatin-sensitive carboxyl proteinases (aspartic proteinases) is hydrophobic in nature. Thus, the hydrophilic nature of the S(2)' subsite was confirmed to be a distinguishing feature of pepstatin-insensitive carboxyl proteinases from prokaryotes.     

Solution structure of the ActD-5'-CCGTT3GTGG-3' complex: drug interaction with tandem G.T mismatches and hairpin loop backbone

Binding of actinomycin D (ActD) to the seemingly single-stranded DNA (ssDNA) oligomer 5'-CCGTT3 GTGG-3' has been studied in solution using high-resolution nuclear magnetic resonance (NMR) techniques. A strong binding constant (8 x 10(6) M(-1)) and high quality NMR spectra have allowed us to determine the initial DNA structure using distance geometry as well as the final ActD-5'-CCGTT3 GTGG-3' complex structure using constrained molecular dynamics calculations. The DNA oligomer 5'-CCGTT3GTGG-3' in the complex forms a hairpin structure with tandem G.T mismatches at the stem region next to a loop of three stacked thymine bases pointing toward the major groove. Bipartite T2O-GH1 and T2O-G2NH2 hydrogen bonds were detected for the G.T mismatches that further stabilize this unusual DNA hairpin. The phenoxazone chromophore of ActD intercalates nicely between the tandem G.T mismatches in essentially one major orientation. Additional hydrophobic interactions between the ActD quinoid amino acid residues with the loop T5-T6-T7 backbone protons were also observed. The hydrophobic G-phenoxazone-G interaction in the ActD-5'-CCGTT3GTGG-3' complex is more robust than that of the classical ActD- 5'-CCGCT3GCGG-3' complex, consistent with the roughly 2-fold stronger binding of ActD to the 5'-CCGTT3GTGG-3' sequence than to its 5'-CCG CT3GCGG-3' counterpart. Stabilization by ActD of a hairpin containing non-canonical stem base pairs further strengthens the notion that ActD or other related compounds may serve as a sequence- specific ssDNA-binding agent that inhibits human immunodeficiency virus (HIV) and other retroviruses replicating through ssDNA intermediates.     

Structure of human prostatic acid phosphatase gene.

Two cDNA clones containing the complete protein-coding sequence of 1,188 nucleotides as well as the 5' and 3' non-coding regions of human prostatic acid phosphatase (PAP) were isolated and sequenced. The size of PAP mRNAs from benign prostate hyperplasia and cancerous prostate was estimated to be 3.2Kb, indicating that the 3' downstream polyadenylation signal was used. Several genomic clones containing parts of the human PAP gene were isolated and the nucleotide sequence of ten exons and their flanking regions was determined. The protein-coding sequence of the human PAP gene was interrupted by nine introns. The positions of all nine introns present in the human PAP gene were homologous to those of the first nine introns in the human lysosomal acid phosphatase (LAP) gene. However, the last (11th) exon of the LAP gene encoding the COOH-terminal domain, which includes a transmembrane segment, was found to be absent in human PAP gene. Southern blot analysis of ten mammalian genomic DNAs gave multiple EcoRI fragments. The data of human genomic DNAs were consistent with the total length of the PAP gene of at least 50 kilobases.     

DNA triplex formation of oligonucleotide analogues consisting of linker groups and octamer segments that have opposite sugar-phosphate backbone polarities

The DNA oligomer analogues 3'd(CTTTCTTT)5'-P4-5'd(TTCTTCTT)3' (IV), 5'd-(TTTCTTTC)3'-P2-3'd(CTTTCTTT)5' (V), and 5'd(TTTCTTTC)3'-P2-3'd(CTTTCTTT)5'-P4-5'd-(TTCTTCTT)3' (VI) (P2 = P*P and P4 = P*P*P*P, where P = phosphate and * = 1,3-propanediol) have been synthesized. These oligomers consist of a linker group or groups and homopyrimidine oligonucleotides which have opposite sugar-phosphate backbone polarities. These oligomer analogues are designed to form triplexes with a duplex, 5'd(AAAGAAAGCCCTTTCTTTAAGAAGAA)3'.5'd(TTCTTCTTAAA- GAAAGGGCTTTCTTT)3' (I), which contains small homopurine clusters alternately located in both strands. The length of the linker groups, P2 and P4, was based upon a computer modeling analysis. Triplex formation by the unlinked octamers 5'd(TTCTTCTT)3' (II) and 5'd(TTTCTTTC)3' (III) and the linked oligomer analogues IV-VI with the target duplex was studied by thermal denaturation at pH 5.2. The order of stabilities of triplex formation by these oligomers was I-V much much greater than I-IV greater than I-(II, III). The mixture of I and VI showed two transitions corresponding to the dissociation of the third strand. The higher transition corresponded to the dissociation of 3'-3'-linked octamer segments, and the lower one corresponded to the dissociation of 5'-5'-linked octamer segments. The Tm of the latter transition was higher than that of the I-IV triplex; thus the triplex formed by the 5'-5'-linked octamer segment was stabilized by the triplex formed by the 3'-3'-linked octamer segments in the I-VI triplex. Triplex formation of this system was also studied in the presence of ethidium bromide.(ABSTRACT TRUNCATED AT 250 WORDS)     

Pyrophosphorolysis by Type II DNA polymerases: implications for pyrophosphorolysis-activated polymerization

We find that Type II DNA polymerases can catalyze pyrophosphorolysis, the reverse reaction of DNA polymerization. This property is applied utilizing pyrophosphorolysis-activated polymerization (PAP), a method of nucleic acid amplification using serial coupling of pyrophosphorolysis and polymerization. PAP can be used for ultrarare allele detection (detection of minimal residual disease and cancer risk assessment through measurement of mutation load) and for microarray-based scanning for unknown mutations. Herein, we show that Type II DNA polymerases efficiently catalyze template-dependent pyrophosphorolysis to activate oligonucleotides blocked at their 3' termini with acyclonucleotides in which a 2-hydroxyethoxymethyl group substitutes for the 2'-deoxyribofuranosyl sugar. Type II archeon DNA polymerases Vent (exo-) and Pfu (exo-) can be utilized for PAP or a bidirectional form of PAP with acyclonucleotide-blocked oligonucleotides, but not with dideoxynucleotide-blocked oligonucleotides. In contrast, a Type I DNA polymerase, TaqFS, can utilize either acyclonucleotide-blocked or dideoxynucleotide-blocked oligonucleotides. These findings expand the potential of nascent PAP technology.     

Modulation of the Pyrococcus abyssi NucS endonuclease activity by replication clamp at functional and structural levels

Pyrococcus abyssi NucS is the founding member of a new family of structure-specific DNA endonucleases that interact with the replication clamp proliferating cell nuclear antigen (PCNA). Using a combination of small angle x-ray scattering and surface plasmon resonance analyses, we demonstrate the formation of a stable complex in solution, in which one molecule of the PabNucS homodimer binds to the outside surface of the PabPCNA homotrimer. Using fluorescent labels, PCNA is shown to increase the binding affinity of NucS toward single-strand/double-strand junctions on 5' and 3' flaps, as well as to modulate the cleavage specificity on the branched DNA structures. Our results indicate that the presence of a single major contact between the PabNucS and PabPCNA proteins, together with the complex-induced DNA bending, facilitate conformational flexibility required for specific cleavage at the single-strand/double-strand DNA junction.     

Evidence for a SAL1-PAP chloroplast retrograde pathway that functions in drought and high light signaling in Arabidopsis

Compartmentation of the eukaryotic cell requires a complex set of subcellular messages, including multiple retrograde signals from the chloroplast and mitochondria to the nucleus, to regulate gene expression. Here, we propose that one such signal is a phosphonucleotide (3'-phosphoadenosine 5'-phosphate [PAP]), which accumulates in Arabidopsis thaliana in response to drought and high light (HL) stress and that the enzyme SAL1 regulates its levels by dephosphorylating PAP to AMP. SAL1 accumulates in chloroplasts and mitochondria but not in the cytosol. sal1 mutants accumulate 20-fold more PAP without a marked change in inositol phosphate levels, demonstrating that PAP is a primary in vivo substrate. Significantly, transgenic targeting of SAL1 to either the nucleus or chloroplast of sal1 mutants lowers the total PAP levels and expression of the HL-inducible ASCORBATE PEROXIDASE2 gene. This indicates that PAP must be able to move between cellular compartments. The mode of action for PAP could be inhibition of 5' to 3' exoribonucleases (XRNs), as SAL1 and the nuclear XRNs modulate the expression of a similar subset of HL and drought-inducible genes, sal1 mutants accumulate XRN substrates, and PAP can inhibit yeast (Saccharomyces cerevisiae) XRNs. We propose a SAL1-PAP retrograde pathway that can alter nuclear gene expression during HL and drought stress.     

Catalytic residues and substrate specificity of scytalidoglutamic peptidase, the first member of the eqolisin in family (G1) of peptidases

Scytalidoglutamic peptidase (SGP) is the first-discovered member of the eqolisin family of peptidases with a unique structure and a presumed novel catalytic dyad (E136 and Q53) [Fujinaga et al., PNAS 101 (2004) 3364-3369]. Mutants of SGP, E136A, Q53A, and Q53E lost both the autoprocessing and enzymatic activities of the wild-type enzyme. Coupled with the results from the structural analysis of SGP, Glu136 and Gln53 were identified as the catalytic residues. The substrate specificity of SGP is unique, particularly, in the preference at the P3 (basic amino acid), P1' (small a.a.), and P3' (basic a.a.) positions. Superior substrates and inhibitors have been synthesized for kinetic studies based on the results reported here. kcat, Km, and kcat/Km of SGP for D-Dap(MeNHBz)-GFKFF*ALRK(Dnp)-D-R-D-R were 34.8 s-1, 0.065 microM, and 535 microM-1 s-1, respectively. Ki of Ac-FKF-(3S,4S)-phenylstatinyl-LR-NH2 for SGP was 1.2x10(-10) M. Taken together, we can conclude that SGP has not only structural and catalytic novelties but also a unique subsite structure.