Rabbit liver tRNA1Val:II. unusual secondary structure of T psi C stem and loop due to a U54:A60 base pair.

In contrast to all other known tRNAs, mammalian tRNA1Val contains two adenosines A59 and A60, opposite to U54 and psi 55 in the U psi CG sequence of the T psi C loop, which could form unusual A:U (or A: psi pairs in addition to the five "normal" G:C pairs. In order to measure the number of G:C and A:U (A: psi) pairs in the T psi C stem, we prepared the 30 nucleotide long 3'-terminal fragment of this tRNA by "m7G-cleavage". From differentiated melting curves and temperature jump experiments it was concluded that the T psi C stem in this fragment is in fact extended by an additional A60:U54 pair. A dimer of this fragment with 14 base pairs was characterized by gel electrophoresis and by the same physical methods. An additional A:U pair in the tRNA1Val fragment does not necessarily mean that this is also true for intact tRNA. However, we showed that U54 is far less available for enzymatic methylation in mammalian tRNA1Val compared to tRNA from T-E. coli. This clear difference in U54 reactivity, together with the identification of an extra A60:U54 pair in the U psi CG containing fragment suggests the presence of a 6 base pair T psi C stem and a 5 nucleotide T psi C loop in this tRNA.     

Specificity and Catalytic Mechanism in Family 5 Uracil DNA Glycosylase

UDGb belongs to family 5 of the uracil DNA glycosylase (UDG) superfamily. Here, we report that family 5 UDGb from Thermus thermophilus HB8 is not only a uracil DNA glycosyase acting on G/U, T/U, C/U, and A/U base pairs, but also a hypoxanthine DNA glycosylase acting on G/I, T/I, and A/I base pairs and a xanthine DNA glycosylase acting on all double-stranded and single-stranded xanthine-containing DNA. Analysis of potentials of mean force indicates that the tendency of hypoxanthine base flipping follows the order of G/I > T/I, A/I > C/I, matching the trend of hypoxanthine DNA glycosylase activity observed in vitro. Genetic analysis indicates that family 5 UDGb can also act as an enzyme to remove uracil incorporated into DNA through the existence of dUTP in the nucleotide pool. Mutational analysis coupled with molecular modeling and molecular dynamics analysis reveals that although hydrogen bonding to O₂ of uracil underlies the UDG activity in a dissociative fashion, Tth UDGb relies on multiple catalytic residues to facilitate its excision of hypoxanthine and xanthine. This study underscores the structural and functional diversity in the UDG superfamily.     

Identification and assignment of base pairs in three helical stems of Torulopsis utilis ribosomal 5S RNA and its RNase T1 and RNase T2 cleaved fragments via 500-MHz proton homonuclear overhauser enhancements

Imino proton resonances in the downfield region (10-14 ppm) of the 500-MHz 1H NMR spectrum of Torulopsis utilis 5S RNA are identified (A X U, G X C, or G X U) and assigned to base pairs in helices I, IV, and V via analysis of homonuclear Overhauser enhancements (NOE) from intact T. utilis 5S RNA, its RNase T1 and RNase T2 digested fragments, and a second yeast (Saccharomyces cerevisiae) 5S RNA whose nucleotide sequence differs at only six residues from that of T. utilis 5S RNA. The near-identical chemical shifts and NOE behavior of most of the common peaks from these four RNAs strongly suggest that helices I, IV, and V retain the same conformation after RNase digestion and that both T. utilis and S. cerevisiae 5S RNAs share a common secondary and tertiary structure. Of the four G X U base pairs identified in the intact 5S RNA, two are assigned to the terminal stem (helix I) and the other two to helices IV and V. Seven of the nine base pairs of the terminal stem have been assigned. Our experimental demonstration of a G X U base pair in helix V supports the 5S RNA secondary structural model of Luehrsen and Fox [Luehrsen, K. R., & Fox, G.E. (1981) Proc. Natl. Acad. Sci. U.S.A. 78, 2150-2154]. Finally, the base-pair proton peak assigned to the terminal G X U in helix V of the RNase T2 cleaved fragment is shifted downfield from that in the intact 5S RNA, suggesting that helices I and V may be coaxial in intact T. utilis 5S RNA.     

Ectopic expression of AID in a non-B cell line triggers A:T and G:C point mutations in non-replicating episomal vectors

Somatic hypermutation (SHM) of immunoglobulin genes is currently viewed as a two step process initiated by the deamination of deoxycytidine (C) to deoxyuridine (U), catalysed by the activation induced deaminase (AID). Phase 1 mutations arise from DNA replication across the uracil residue or the abasic site, generated by the uracil-DNA glycosylase, yielding transitions or transversions at G:C pairs. Phase 2 mutations result from the recognition of the U:G mismatch by the Msh2/Msh6 complex (MutS Homologue), followed by the excision of the mismatched nucleotide and the repair, by the low fidelity DNA polymerase eta, of the gap generated by the exonuclease I. These mutations are mainly focused at A:T pairs. Whereas in activated B cells both G:C and A:T pairs are equally targeted, ectopic expression of AID was shown to trigger only G:C mutations on a stably integrated reporter gene. Here we show that when using non-replicative episomal vectors containing a GFP gene, inactivated by the introduction of stop codons at various positions, a high level of EGFP positive cells was obtained after transient expression in Jurkat cells constitutively expressing AID. We show that mutations at G:C and A:T pairs are produced. EGFP positive cells are obtained in the absence of vector replication demonstrating that the mutations are dependent only on the mismatch repair (MMR) pathway. This implies that the generation of phase 1 mutations is not a prerequisite for the expression of phase 2 mutations.     

Effects of flanking G . C base pairs on internal Watson-Crick, G . U, and nonbonded base pairs within a short ribonucleic acid duplex

A series of pentaribonucleotides, ApGpXpGpU (where X identical to A, G, C, or U), was synthesized to investigate the effects of flanking G . C pairs on internal Watson-Crick, G . U, and nonbonded base pairs. Sequences ApGpApCpU (Tm = 26 degrees C) and ApGpCpCpU (Tm = 25 degrees C) were each found to form a duplex with non-base-paired internal residues that stacked with the rest of the sequence but were not looped out. ApGpGpCpU also forms a duplex (Tm = 30 degrees C) but with dangling terminal nonbonded adenosines rather than internal nonbonded guanosines. ApGpUpCpU prefers a stacked single-strand conformation. In addition, contribution to duplex stability from an internal A . U or G . C base pair is enhanced by 6 degrees C when flanked by G . C base pairs as compared to A . U base pairs. G . C base pairs flanking an internal G . U base pair were found to be more tolerant to the altered conformation of a G . U pair and result in an increase to stability comparable with that found for an internal A . U base pair.     

Spectroscopic evidence for the uneven distribution of adenine and uracil residues in ribosomal ribonucleic acid of Drosophila melanogaster and of Plasmodium knowlesi and its possible evolutionary significance.

RNA was isolated from subribosomal particles of the malaria parasite Plasmodium knowlesi. The nucleotide composition (mole fraction) of the principal species was obtained (S-rRNA, 0.295A, 0.36U, 0.25G, 0.105C: L-rRNA, 0.326A, 0.31U, 0.228G, 0.144C). Ribosomal RNA was also isolated from Drosophila melanogaster. Optical properties of these A + U-rich species were measured. In all four cases analysis of the hypochromic effect revealed that adenine and uracil residues tended to form clusters along the polynucleotide chain. A substantial fraction of residues was located in bihelical regions of approx. 50% G-C base pairs or in regions of approx. 30-35% G-C base pairs. The possible evolutionary significance of these results was considered on the basis of comparison with properties of rRNA from bacteria (Escherichia coli) and a mammal (rabbit reticulocyte).     

The identity determinants required for the discrimination between tRNAGlu and tRNAAsp by glutamyl-tRNA synthetase from Escherichia coli.

We previously elucidated the major determinant set for Escherichia coli tRNAGlu identity (U34, U35, C36, A37, G1*C72, U2*A71, U11*A24, U13*G22**Alpha46, and Delta47) and showed that the set is sufficient to switch the identity of tRNAGln to Glu [Sekine, S., Nureki, O., Sakamoto, K., Niimi, T., Tateno, M., Go, M., Kohno, T., Brisson, A., Lapointe, J. & Yokoyama, S. (1996) J. Mol. Biol. 256, 685-700]. In the present study, we attempted to switch the identity of tRNAAsp, which has a sequence similar to that of tRNAGlu, and consequently possesses many nucleotide residues corresponding to the Glu identity determinants (U35, C36, A37, G1*C72, and U11*A24). A simple transplantation of the rest of the major determinants (U34, U2*A71, U13*G22**Alpha46, and Delta47) to the framework of tRNAAsp did not result in a sufficient switch of the tRNAAsp identity to Glu. To confer an optimal glutamate accepting activity to tRNAAsp, two other elements, C4*G69 in the middle of the acceptor stem and C12*G23**C9 in the augmented D helix, were required. Consistently, the two base pairs, C4*G69 and C12*G23, in tRNAGlu had been shown to exist in the interface with glutamyl-tRNA synthetase (GluRS) by phosphate-group footprinting. We also found the two elements in the framework of tRNAGln, and determined that their contributions successfully changed the identity of tRNAGln to Glu in the previous study. By the identity-determinant set (C4*G69 and C12*G23**C9 in addition to U34, U35, C36, A37, G1*C72, U2*A71, U11*A24, U13*G22**Alpha46, and Delta47) the activity of GluRS was optimized and efficient discrimination from the noncognate tRNAs was achieved.     

Incorporation of dUTP does not mediate mutation of A:T base pairs in Ig genes in vivo

Activation-induced cytidine deaminase (AID) protein initiates Ig gene mutation by deaminating cytosines, converting them into uracils. Excision of AID-induced uracils by uracil-N-glycosylase is responsible for most transversion mutations at G:C base pairs. On the other hand, processing of AID-induced G:U mismatches by mismatch repair factors is responsible for most mutation at Ig A:T base pairs. Why mismatch processing should be error prone is unknown. One theory proposes that long patch excision in G1-phase leads to dUTP-incorporation opposite adenines as a result of the higher G1-phase ratio of nuclear dUTP to dTTP. Subsequent base excision at the A:U base pairs produced could then create non-instructional templates leading to permanent mutations at A:T base pairs (1). This compelling theory has remained untested. We have developed a method to rapidly modify DNA repair pathways in mutating mouse B cells in vivo by transducing Ig knock-in splenic mouse B cells with GFP-tagged retroviruses, then adoptively transferring GFP⁺ cells, along with appropriate antigen, into primed congenic hosts. We have used this method to show that dUTP-incorporation is unlikely to be the cause of AID-induced mutation of A:T base pairs, and instead propose that A:T mutations might arise as an indirect consequence of nucleotide paucity during AID-induced DNA repair.     

Hypermutation at A/T sites during G.U mismatch repair in vitro by human B-cell lysates

Somatic hypermutation in the variable regions of immunoglobulin genes is required to produce high affinity antibody molecules. Somatic hypermutation results by processing G.U mismatches generated when activation-induced cytidine deaminase (AID) deaminates C to U. Mutations at C/G sites are targeted mainly at deamination sites, whereas mutations at A/T sites entail error-prone DNA gap repair. We used B-cell lysates to analyze salient features of somatic hypermutation with in vitro mutational assays. Tonsil and hypermutating Ramos B-cells convert C-->U in accord with AID motif specificities, whereas HeLa cells do not. Using tonsil cell lysates to repair a G.U mismatch, A/T and G/C targeted mutations occur about equally, whereas Ramos cell lysates make fewer mutations at A/T sites (approximately 24%) compared with G/C sites (approximately 76%). In contrast, mutations in HeLa cell lysates occur almost exclusively at G/C sites (> 95%). By recapitulating two basic features of B-cell-specific somatic hypermutation, G/C mutations targeted to AID hot spot motifs and elevated A/T mutations dependent on error-prone processing of G.U mispairs, these cell free assays provide a practical method to reconstitute error-prone mismatch repair using purified B-cell proteins.     

Structure of SL4 RNA from the HIV-1 packaging signal.

The NMR-based structure is described for an RNA model of stem-loop 4 (SL4) from the HIV-1 major packaging domain. The GAGA tetraloop adopts a conformation similar to the classic GNRA form, although there are differences in the details. The type II tandem G.U pairs have a combination of wobble and bifurcated hydrogen bonds where the uracil 2-carbonyl oxygen is hydrogen-bonded to both G,H1 and G,H2. There is the likelihood of a Na(+) ion coordinated to the four carbonyl oxygens in the major groove for these G.U pairs and perhaps to the N7 lone pairs of the G bases as well. A continuous stack of five bases extends over nearly the whole length of the stem to the base of the loop in the RNA 16mer: C15/U14/G13/G5/C6. There is no evidence for a terminal G.A pair; instead, G1 appears quite unrestrained, and A16 stacks on both C15 and G2. Residues G2 through G5 exhibit broadened resonances, especially G3 and U4, suggesting enhanced mobility for the 5'-side of the stem. The structure shows G2/G3/U4 stacking along the same strand, nearly isolated from interaction with the other bases. This is probably an important factor in the signal broadening and apparent mobility of these residues and the low stability of the 16mer hairpin against thermal denaturation.     

Theoretical studies on protein-nucleic acid interactions. II. Hydrogen bonding of amino acid side chains with bases and base pairs of nucleic acids

With a view to understanding the role of hydrogen bonds in the recognition of nucleic acids by proteins, hydrogen bonding between the bases and base pairs of nucleic acids and the amino acids (Asn, Gln, Asp and Glu, and charged residues Arg⁺, Glu^?, and Asp^?) has been studied by a second-order perturbation theory. Binding energies have been calculated for all possible configurations involving a pair of hydrogen bonds between the base (or base pair) and the amino acid residue. Our results show that the hydrogen bonding in these cases has a large contribution from electrostatic interaction. In general, the charged amino acids, compared to the uncharged ones, form more stable complexes with bases or base pairs. The hydrogen-bond energies are an order of magnitude smaller than the Coulombic interaction energies between basic amino acids (Lys⁺, Arg⁺, and His⁺) and the phosphate groups of nucleic acids. The stabilities of the complexes of amino acids Asn, Gln, Asp, and Glu with bases are in the order: G–X > C–X > A–X U–X or T–X, and G · C–X > A · T(U)–X, where X is one of these amino acid residues. It has been shown that Glu^? and Asp^? can recognize guanine in single-stranded nucleic acids; Arg⁺ can recognize G · C base pairs from A · T base pairs in double-stranded structures.     

Structural basis for type I and type II deficiencies of antithrombotic plasma protein C: Patterns revealed by three-dimensional molecular modelling of mutations of the protease domain

Familial deficiency of protein C is associated with inherited thrombophilia. To explore how specific missense mutations might cause observed clinical phenotypes, known protein C missense mutations were mapped onto three-dimensional homology models of the protein C protease domain, and the implications for domain folding and structure were evaluated. Most Type I missense mutations either replaced internal hydrophobic residues (I201T, L223F, A259V, A267T, A346T, A346V, G376D) or nearby interacting residues (I403M, T298M, Q184H), thus disrupting the packing of internal hydrophobic side chains, or changed hydrophilic residues, thus disrupting ion pairs (N256D, R178W). Mutations (P168L, R169W) at the activation site destabilized the region containing the activation peptide structure. Most Type II mutations involved solvent-exposed residues and were clustered either in a positively charged region (R147W, R157Q, R229Q, R352W) or were located in or near the active site region (S252N, D359N, G381S, G391S, H211Q). The cluster of arginines 147, 157, 229, and 352 may identify a functionally important exosite. Identification of the spatial relationships of natural mutations in the protein C model is helpful for understanding manifestations of protein C deficiency and for identification of novel, functionally important molecular features and exosites. © 1994 John Wiley & Sons, Inc.     

Zipper-like Watson-Crick base-pairs

A series of DNA heptadecamers containing the DNA analogues of RNA E-like 5'-d(GXA)/(AYG)-5' motifs (X/Y is complementary T/A, A/T, C/G, or G/C pair) were studied using nuclear magnetic resonance (NMR) methodology and distance geometry (DG)/molecular dynamics (MD) approaches. Such oligomers reveal excellent resolution in NMR spectra and exhibit many unusual nuclear Overhauser effects (NOEs) that allow for good characterization of an unusual zipper-like conformation with zipper-like Watson-Crick base-pairs; the potential canonical X.Y H-bonding is not present, and the central X/Y pairs are transformed instead into inter-strand stacks that are bracketed by sheared G.A base-pairs. Such phenomenal structural change is brought about mainly through two backbone torsional angle adjustments, i.e. delta from C2'-endo to C3'-endo for the sugar puckers of unpaired residues and gamma from gauche(+) to trans for the following 3'-adenosine residues. Such motifs are analogous to the previously studied (GGA)(2) motif presumably present in the human centromeric (TGGAA)(n) tandem repeat sequence. The novel zipper-like motifs are only 4-7 deg. C less stable than the (GGA)(2) motif, suggesting that inter-strand base stacking plays an important role in stabilizing unusual nucleic acid structures. The discovery that canonical Watson-Crick G.C or A.T hydrogen-bonded pairs can be transformed into stacking pairs greatly increases the repertoire for unusual nucleic acid structural motifs.     

Crystal structure of human selenocysteine tRNA

Selenocysteine (Sec) is the 21st amino acid in translation. Sec tRNA (tRNA^Sec) has an anticodon complementary to the UGA codon. We solved the crystal structure of human tRNA^Sec. tRNA^Sec has a 9-bp acceptor stem and a 4-bp T stem, in contrast with the 7-bp acceptor stem and the 5-bp T stem in the canonical tRNAs. The acceptor stem is kinked between the U6:U67 and G7:C66 base pairs, leading to a bent acceptor-T stem helix. tRNA^Sec has a 6-bp D stem and a 4-nt D loop. The long D stem includes unique A14:U21 and G15:C20a pairs. The D-loop:T-loop interactions include the base pairs G18:U55 and U16:U59, and a unique base triple, U20:G19:C56. The extra arm comprises of a 6-bp stem and a 4-nt loop. Remarkably, the D stem and the extra arm do not form tertiary interactions in tRNA^Sec. Instead, tRNA^Sec has an open cavity, in place of the tertiary core of a canonical tRNA. The linker residues, A8 and U9, connecting the acceptor and D stems, are not involved in tertiary base pairing. Instead, U9 is stacked on the first base pair of the extra arm. These features might allow tRNA^Sec to be the target of the Sec synthesis/incorporation machineries.     

Mismatch-induced conformational distortions in polymerase beta support an induced-fit mechanism for fidelity

Molecular dynamics simulations of DNA polymerase (pol) beta complexed with different incorrect incoming nucleotides (G x G, G x T, and T x T template base x incoming nucleotide combinations) at the template-primer terminus are analyzed to delineate structure-function relationships for aberrant base pairs in a polymerase active site. Comparisons, made to pol beta structure and motions in the presence of a correct base pair, are designed to gain atomically detailed insights into the process of nucleotide selection and discrimination. In the presence of an incorrect incoming nucleotide, alpha-helix N of the thumb subdomain believed to be required for pol beta's catalytic cycling moves toward the open conformation rather than the closed conformation as observed for the correct base pair (G x C) before the chemical reaction. Correspondingly, active-site residues in the microenvironment of the incoming base are in intermediate conformations for non-Watson-Crick pairs. The incorrect incoming nucleotide and the corresponding template residue assume distorted conformations and do not form Watson-Crick bonds. Furthermore, the coordination number and the arrangement of ligands observed around the catalytic and nucleotide binding magnesium ions are mismatch specific. Significantly, the crucial nucleotidyl transferase reaction distance (P(alpha)-O3') for the mismatches between the incoming nucleotide and the primer terminus is not ideally compatible with the chemical reaction of primer extension that follows these conformational changes. Moreover, the extent of active-site distortion can be related to experimentally determined rates of nucleotide misincorporation and to the overall energy barrier associated with polymerase activity. Together, our studies provide structure-function insights into the DNA polymerase-induced constraints (i.e., alpha-helix N conformation, DNA base pair bonding, conformation of protein residues in the vicinity of dNTP, and magnesium ions coordination) during nucleotide discrimination and pol beta-nucleotide interactions specific to each mispair and how they may regulate fidelity. They also lend further support to our recent hypothesis that additional conformational energy barriers are involved following nucleotide binding but prior to the chemical reaction.     

Reverse polarization in amino acid and nucleotide substitution patterns between human-mouse orthologs of two compositional extrema.

Genome-wide analysis of sequence divergence patterns in 12,024 human-mouse orthologous pairs reveals, for the first time, that the trends in nucleotide and amino acid substitutions in orthologs of high and low GC composition are highly asymmetric and polarized to opposite directions. The entire dataset has been divided into three groups on the basis of the GC content at third codon sites of human genes: high, medium, and low. High-GC orthologs exhibit significant bias in favor of the replacements, Thr --> Ala, Ser --> Ala, Val --> Ala, Lys --> Arg, Asn --> Ser, Ile --> Val etc., from mouse to human, whereas in low-GC orthologs, the reverse trends prevail. In general, in the high-GC group, residues encoded by A/U-rich codons of mouse proteins tend to be replaced by the residues encoded by relatively G/C-rich codons in their human orthologs, whereas the opposite trend is observed among the low-GC orthologous pairs. The medium-GC group shares some trends with high-GC group and some with low-GC group. The only significant trend common in all groups of orthologs, irrespective of their GC bias, is (Asp)(Mouse) --> (Glu)(Human) replacement. At the nucleotide level, high-GC orthologs have undergone a large excess of (A/T)(Mouse) --> (G/C)(Human) substitutions over (G/C)(Mouse) --> (A/T)(Human) at each codon position, whereas for low-GC orthologs, the reverse is true.     

The longest (A+T) and (G+C) blocks in the human and other genomes

We analysed complete or almost complete nucleotide sequences of the human, chimp, mouse, rat, chicken, dog, and other genomes to find that they contain extremely long (A+T) a (G+C) blocks that do not occur at all in the corresponding randomized sequences. The longest is an (A+T) block containing 1040 consecutive AT pairs that occurs in the 16th human chromosome. The longest human (G+C) block has 261 bp in length. About a half of the longest blocks occur in introns. The (A+T) blocks are discrete units whereas the (G+C) blocks are diffuse. They are imbedded in the genome through connectors longer than 1 kilobase where the (G+C) content gradually decreases to the value of 50%. Remarkably, the (A+T) as well as (G+C) blocks are substantially shorter in the chimp genome. Chicken is characteristic by very long (G+C) blocks that are even longer than in the human genome. Though much shorter, long (G+C) and especially (A+T) blocks occur in lower organisms as well, which means that AT and GC pair clustering is an ancient property that has evolved into large scales in higher eukaryote genomes and the human genome in particular. Very long (A+T) and (G+C) blocks confer specific biophysical properties on DNA that are likely to influence genome folding in cell nuclei and its functional properties.