The A2453-C2499 wobble base pair in Escherichia coli 23S ribosomal RNA is responsible for pH sensitivity of the peptidyltransferase active site conformation

Peptide bond formation, catalyzed by the ribosomal peptidyltransferase, has long been known to be sensitive to monovalent cation concentrations and pH. More recently, we and others have shown that residue A2451 in the peptidyltransferase center of the Escherichia coli 50S ribosomal subunit changes conformation in response to alterations in pH, depending on ionic conditions and temperature. Two wobble pairs, A2453-C2499 and A2450-C2063, have been proposed as potential candidates to convey pH-dependent flexibility to the peptidyltransferase center. Each is presumed to possess a near-neutral pK_a, and both lie in proximity to A2451. We show through mutagenesis and chemical probing that the identity of the A2453-C2499 base pair, but not the A2450-C2063 base pair, is critical for the pH-dependent structural rearrangement of A2451. We conclude that, while the A2453-C2499 base pair may be important for maintaining the structure of the active site in the E.coli peptidyltransferase center, its lack of conservation makes it, and consequently its near-neutral pK_a, unlikely to contribute to function during peptide bond formation.     

Constraints on error rate revealed by computational study of G•U tautomerization in translation

In translation, G•U mismatch in codon-anticodon decoding is an error hotspot likely due to transition of G•U from wobble (wb) to Watson-Crick (WC) geometry, which is governed by keto/enol tautomerization (wb-WC reaction). Yet, effects of the ribosome on the wb-WC reaction and its implications for decoding mechanism remain unclear. Employing quantum-mechanical/molecular-mechanical umbrella sampling simulations using models of the ribosomal decoding site (A site) we determined that the wb-WC reaction is endoergic in the open, but weakly exoergic in the closed A-site state. We extended the classical ‘induced-fit’ model of initial selection by incorporating wb-WC reaction parameters in open and closed states. For predicted parameters, the non-equilibrium exoergic wb-WC reaction is kinetically limited by the decoding rates. The model explains early observations of the WC geometry of G•U from equilibrium structural studies and reveals discrimination capacity for the working ribosome operating at non-equilibrium conditions. The equilibration of the exoergic wb-WC reaction counteracts the equilibration of the open-closed transition of the A site, constraining the decoding accuracy and potentially explaining the persistence of the G•U as an error hotspot. Our results unify structural and mechanistic views of codon-anticodon decoding and generalize the ‘induced-fit’ model for flexible substrates.     

Impact of base pair identity 5′ to the spliceosomal branch site adenosine on branch site conformation

The branch site helix from Saccharomyces cerevisiae with pseudouridine (ψ) incorporated in a phylogenetically conserved position of U2 snRNA features an extrahelical branch site adenosine (A) that forms a base triple interaction with the minor groove edge of a widely conserved purine_{U2 strand}-pyrimidine_{intron strand} (R_U2-Y_intron) base pair two positions upstream. In these studies, NMR spectra of a duplex in which 2-aminopurine (2ap), a fluorescent analog of adenine lacking the proposed hydrogen bond donor, was substituted for the branch site A, indicated that the substitution does not alter the extrahelical position of the branch site residue; thus, it appears that a hydrogen bond between the adenine amino group and the R-Y pair is not obligatory for stabilization of the extrahelical conformation. In contrast, reversal of the orientation of A_U2-U_intron to U_U2-A_intron resulted in an intrahelical position for the branch site A or 2ap. Fluorescence intensity of 2ap substituted for the branch site A with the original R_U2-Y_intron orientation (AU or GC) was high, consistent with an extrahelical position, whereas fluorescence in helices with the reversed R-Y orientation, or with a mismatched pair (A-U → G•A or U•C), was markedly quenched, implying that the residue was stacked in the helix. The A 5′ to the branch site residue was not extrahelical in any of the duplexes. These findings suggest that the R_U2-Y_intron base pair orientation in the ψ-dependent branch site helix plays an important role in positioning the branch site A for recognition and/or function.     

Ion-induced folding of a kink turn that departs from the conventional sequence

Kink turns (k-turns) are important structural motifs that create a sharp axial bend in RNA. Most conform to a consensus in which a three-nucleotide bulge is followed by consecutive G•A and A•G base pairs, and when these G•A pairs are modified in vitro this generally leads to a failure to adopt the k-turn conformation. Kt-23 in the 30S ribosomal subunit of Thermus thermophilus is a rare exception in which the bulge-distal A•G pair is replaced by a non-Watson–Crick A•U pair. In the context of the ribosome, Kt-23 adopts a completely conventional k-turn geometry. We show here that this sequence is induced to fold into a k-turn structure in an isolated RNA duplex by Mg²⁺ or Na⁺ ions. Therefore, the Kt-23 is intrinsically stable despite lacking the key A•G pair; its formation requires neither tertiary interactions nor protein binding. Moreover, the Kt-23 k-turn is stabilized by the same critical hydrogen-bonding interactions within the core of the structure that are found in more conventional sequences such as the near-consensus Kt-7. T. thermophilus Kt-23 has two further non-Watson–Crick base pairs within the non-canonical helix, three and four nucleotides from the bulge, and we find that the nature of these pairs influences the ability of the RNA to adopt k-turn conformation, although the base pair adjacent to the A•U pair is more important than the other.     

Puromycin binding to the donor (P-) site of Escherichia coli ribosomes

Puromycin inhibits the interaction of peptidyl-tRNA analogues AcPhe-tRNA_ox-red^Phe, AcPhe-tRNA^Phe and fMet-tRNA_f^Met with the donor (P-) site of Escherichia coli ribosomes. affects almost equally both the rate of the binding and the equilibrium of the system. This means that the effect is due to direct competition for the P-site, but not due to the indirect influence via the acceptor (A-) site. The inhibition was observed also in 30 S ribosomal subunits, therefore the puromycin binding site is situated far from the peptidyl transferase center. Quantitative measurements show that the affinity of puromycin for its new ribosomal binding site is similar to its affinity for the acceptor site of the peptidyl transferase center.     

Effect of 6-thioguanine on the stability of duplex DNA

The incorporation of 6-thioguanine (S6G) into DNA is a prerequisite for its cytotoxic action, but duplex structure is not significantly perturbed by the presence of the lesion [J. Bohon and C. R. de los Santos (2003) Nucleic Acids Res., 31, 1331–1338]. It is therefore possible that the mechanism of cytotoxicity relies on a loss of stability rather than a pathway involving direct structural recognition. The research described here focuses on the changes in thermodynamic properties of duplex DNA owing to the introduction of S6G as well as the kinetic properties of base pairs involving S6G. Replacement of a guanine in a G•C pair by S6G results in ~1 kcal/mol less favorable Gibbs free energy of duplex formation at 37°C. S6G•T and G•T mismatch-containing duplexes have almost identical Gibbs free energy at 37°C, with values ~3 kcal/mol less favorable than that of the control. Base pair stability is affected by S6G. The lifetime of the normal G•C base pair is ~125 ms, whereas that of the G•T mismatch is below the detection limit. The lifetimes of S6G•C and S6G•T pairs are ~7 and 2 ms, respectively, demonstrating that, although S6G significantly decreases the stability of the pairing with cytosine, it slightly increases that of a mismatch.     

An innate twist between Crick’s wobble and Watson-Crick base pairs

Non-Watson-Crick pairs like the G·U wobble are frequent in RNA duplexes. Their geometric dissimilarity (nonisostericity) with the Watson-Crick base pairs and among themselves imparts structural variations decisive for biological functions. Through a novel circular representation of base pairs, a simple and general metric scheme for quantification of base-pair nonisostericity, in terms of residual twist and radial difference that can also envisage its mechanistic effect, is proposed. The scheme is exemplified by G·U and U·G wobble pairs, and their predicable local effects on helical twist angle are validated by MD simulations. New insights into a possible rationale for contextual occurrence of G·U and other non-WC pairs, as well as the influence of a G·U pair on other non-Watson-Crick pair neighborhood and RNA-protein interactions are obtained from analysis of crystal structure data. A few instances of RNA-protein interactions along the major groove are documented in addition to the well-recognized interaction of the G·U pair along the minor groove. The nonisostericity-mediated influence of wobble pairs for facilitating helical packing through long-range interactions in ribosomal RNAs is also reviewed.     

Biochemical defects in minor spliceosome function in the developmental disorder MOPD I

Biallelic mutations of the human RNU4ATAC gene, which codes for the minor spliceosomal U4atac snRNA, cause the developmental disorder, MOPD I/TALS. To date, nine separate mutations in RNU4ATAC have been identified in MOPD I patients. Evidence suggests that all of these mutations lead to abrogation of U4atac snRNA function and impaired minor intron splicing. However, the molecular basis of these effects is unknown. Here, we use a variety of in vitro and in vivo assays to address this question. We find that only one mutation, 124G>A, leads to significantly reduced expression of U4atac snRNA, whereas four mutations, 30G>A, 50G>A, 50G>C and 51G>A, show impaired binding of essential protein components of the U4atac/U6atac di-snRNP in vitro and in vivo. Analysis of MOPD I patient fibroblasts and iPS cells homozygous for the most common mutation, 51G>A, shows reduced levels of the U4atac/U6atac.U5 tri-snRNP complex as determined by glycerol gradient sedimentation and immunoprecipitation. In this report, we establish a mechanistic basis for MOPD I disease and show that the inefficient splicing of genes containing U12-dependent introns in patient cells is due to defects in minor tri-snRNP formation, and the MOPD I-associated RNU4ATAC mutations can affect multiple facets of minor snRNA function.     

Different Biochemical Properties Explain Why Two Equivalent Gα Subunit Mutants Cause Unrelated Diseases

There is an increasing number of disease-associated Gα mutations identified from genome-wide sequencing campaigns or targeted efforts. Albright''s Hereditary Osteodystrophy (AHO) was the first inherited disease associated with loss-of-function mutations in a G protein (Gαs) and other studies revealed gain-of-function Gα mutations in cancer. Here we attempted to solve the apparent quandary posed by the fact that the same mutation in two different G proteins appeared associated with both AHO and cancer. We first confirmed the presence of an inherited Gαs-R265H mutation from a previously described clinical case report of AHO. This mutation is structurally analogous to Gαo-R243H, an oncogenic mutant with increased activity in vitro and in cells due to rapid nucleotide exchange. We found that, contrary to Gαo-R243H, Gαs-R265H activity is compromised due to greatly impaired nucleotide binding in vitro and in cells. We obtained equivalent results when comparing another AHO mutation in Gαs (D173N) with a counterpart cancer mutation in Gαo (D151N). Gαo-R243H binds nucleotides efficiently under steady-state conditions but releases GDP much faster than the WT protein, suggesting diminished affinity for the nucleotide. These results indicate that the same disease-linked mutation in two different G proteins affects a common biochemical feature (nucleotide affinity) but to a different grade depending on the G protein (mild decrease for Gαo and severe for Gαs). We conclude that Gαs-R265H has dramatically impaired nucleotide affinity leading to the loss-of-function in AHO whereas Gαo-R243H has a mild decrease in nucleotide affinity that causes rapid nucleotide turnover and subsequent hyperactivity in cancer.     

Investigation of the ribosomal peptidyl transferase center using a photoaffinity label

The synthesis of a peptidyl-tRNA photoaffinity analog, 2-nitro-4-azidophenoxy-4′-phenylacetyl-phenylalanyl-tRNA^Phe is described. Covalent attachment of this analog to Escherichia coli 70 S ribosomes requires poly(U)-stimulated binding prior to photolysis. Peptidyl site binding is indicated by the ability of puromycin to release the peptidyl moiety from non-photolyzed samples. Covalently attached 2-nitro-4-azidophenoxy-4-phenylacetyl-Phe-tRNA^Phe can subsequently participate in peptidyl transfer with [³H]Phe-tRNA^Phe bound at the aminoacyl site. This means that the covalent reaction does not produce sufficient distortion of the peptidyl site and its bound tRNA to inactivate the peptidyl transference. If peptidyl transfer with [³H]Phe-tRNA^Phe is allowed to proceed before photolysis, covalent reaction can still occur. In all cases, the main reaction products are two 50 S ribosomal proteins, L11 and L18. The results strongly indicate that these two proteins either form part of the peptidyl transferase center or are located adjacent to it. Presumably, α-halocarbonyl affinity reagents used previously failed to identify these two proteins because they lack easily accessible, reactive nucleophilic groups.     

NMR structure and Mg2+ binding of an RNA segment that underlies the L7/L12 stalk in the E.coli 50S ribosomal subunit

Helix 42 of Domain II of Escherichia coli 23S ribosomal RNA underlies the L7/L12 stalk in the ribosome and may be significant in positioning this feature relative to the rest of the 50S ribosomal subunit. Unlike the Haloarcula marismortui and Deinococcus radiodurans examples, the lower portion of helix 42 in E.coli contains two consecutive G•A oppositions with both adenines on the same side of the stem. Herein, the structure of an analog of positions 1037–1043 and 1112–1118 in the helix 42 region is reported. NMR spectra and structure calculations support a cis Watson–Crick/Watson–Crick (cis W.C.) G•A conformation for the tandem (G•A)₂ in the analog and a minimally perturbed helical duplex stem. Mg²⁺ titration studies imply that the cis W.C. geometry of the tandem (G•A)₂ probably allows O6 of G20 and N1 of A4 to coordinate with a Mg²⁺ ion as indicated by the largest chemical shift changes associated with the imino group of G20 and the H8 of G20 and A4. A cross-strand bridging Mg²⁺ coordination has also been found in a different sequence context in the crystal structure of H.marismortui 23S rRNA, and therefore it may be a rare but general motif in Mg²⁺ coordination.     

Puromycin interacts with the donor (P) site of Escherichia coli ribosomes

Puromycin inhibits the interaction of peptidyl-tRNA analogs AcPhe-tRNA Phe ox-red, AcPhe-tRNA Phe and FMet-tRNA f Met with the donor (P) site of Escherichia coli ribosomes. It affects both template-free and poly(U)-dependent systems. The inhibition is apparently due to direct competition for the P-site. On isolated 30S ribosomal subunits it was shown that the puromycin binding site is situated far from the peptidyl transferase center. Quantitative measurements of the inhibition revealed that the affinity constant of puromycin for the P-site is not less than its affinity for the A-moiety of the peptidyl transferase center [1.1 divided by 3.8) X 10(3) M-1).     

Quadruplex formation by both G-rich and C-rich DNA strands of the C9orf72 (GGGGCC)8?(GGCCCC)8 repeat: effect of CpG methylation

Unusual DNA/RNA structures of the C9orf72 repeat may participate in repeat expansions or pathogenesis of amyotrophic lateral sclerosis and frontotemporal dementia. Expanded repeats are CpG methylated with unknown consequences. Typically, quadruplex structures form by G-rich but not complementary C-rich strands. Using CD, UV and electrophoresis, we characterized the structures formed by (GGGGCC)8 and (GGCCCC)8 strands with and without 5-methylcytosine (5mCpG) or 5-hydroxymethylcytosine (5hmCpG) methylation. All strands formed heterogenous mixtures of structures, with features of quadruplexes (at pH 7.5, in K⁺, Na⁺ or Li⁺), but no feature typical of i-motifs. C-rich strands formed quadruplexes, likely stabilized by G•C•G•C-tetrads and C•C•C•C-tetrads. Unlike G•G•G•G-tetrads, some G•C•G•C-tetrad conformations do not require the N7-Guanine position, hence C9orf72 quadruplexes still formed when N7-deazaGuanine replace all Guanines. 5mCpG and 5hmCpG increased and decreased the thermal stability of these structures. hnRNPK, through band-shift analysis, bound C-rich but not G-rich strands, with a binding preference of unmethylated > 5hmCpG > 5mCpG, where methylated DNA-protein complexes were retained in the wells, distinct from unmethylated complexes. Our findings suggest that for C-rich sequences interspersed with G-residues, one must consider quadruplex formation and that methylation of quadruplexes may affect epigenetic processes.     

Functional correction of CFTR mutations in human airway epithelial cells using adenine base editors

Mutations in the CFTR gene that lead to premature stop codons or splicing defects cause cystic fibrosis (CF) and are not amenable to treatment by small-molecule modulators. Here, we investigate the use of adenine base editor (ABE) ribonucleoproteins (RNPs) that convert A•T to G•C base pairs as a therapeutic strategy for three CF-causing mutations. Using ABE RNPs, we corrected in human airway epithelial cells premature stop codon mutations (R553X and W1282X) and a splice-site mutation (3849 + 10 kb C > T). Following ABE delivery, DNA sequencing revealed correction of these pathogenic mutations at efficiencies that reached 38–82% with minimal bystander edits or indels. This range of editing was sufficient to attain functional correction of CFTR-dependent anion channel activity in primary epithelial cells from CF patients and in a CF patient-derived cell line. These results demonstrate the utility of base editor RNPs to repair CFTR mutations that are not currently treatable with approved therapeutics.     

Direct estimation of the mitochondrial DNA mutation rate in Drosophila melanogaster

Mitochondrial DNA (mtDNA) variants are widely used in evolutionary genetics as markers for population history and to estimate divergence times among taxa. Inferences of species history are generally based on phylogenetic comparisons, which assume that molecular evolution is clock-like. Between-species comparisons have also been used to estimate the mutation rate, using sites that are thought to evolve neutrally. We directly estimated the mtDNA mutation rate by scanning the mitochondrial genome of Drosophila melanogaster lines that had undergone approximately 200 generations of spontaneous mutation accumulation (MA). We detected a total of 28 point mutations and eight insertion-deletion (indel) mutations, yielding an estimate for the single-nucleotide mutation rate of 6.2 × 10⁻⁸ per site per fly generation. Most mutations were heteroplasmic within a line, and their frequency distribution suggests that the effective number of mitochondrial genomes transmitted per female per generation is about 30. We observed repeated occurrences of some indel mutations, suggesting that indel mutational hotspots are common. Among the point mutations, there is a large excess of G→A mutations on the major strand (the sense strand for the majority of mitochondrial genes). These mutations tend to occur at nonsynonymous sites of protein-coding genes, and they are expected to be deleterious, so do not become fixed between species. The overall mtDNA mutation rate per base pair per fly generation in Drosophila is estimated to be about 10× higher than the nuclear mutation rate, but the mitochondrial major strand G→A mutation rate is about 70× higher than the nuclear rate. Silent sites are substantially more strongly biased towards A and T than nonsynonymous sites, consistent with the extreme mutation bias towards A+T. Strand-asymmetric mutation bias, coupled with selection to maintain specific nonsynonymous bases, therefore provides an explanation for the extreme base composition of the mitochondrial genome of Drosophila.     

Essential mechanisms in the catalysis of peptide bond formation on the ribosome

Peptide bond formation is the main catalytic function of the ribosome. The mechanism of catalysis is presumed to be highly conserved in all organisms. We tested the conservation by comparing mechanistic features of the peptidyl transfer reaction on ribosomes from Escherichia coli and the Gram-positive bacterium Mycobacterium smegmatis. In both cases, the major contribution to catalysis was the lowering of the activation entropy. The rate of peptide bond formation was pH independent with the natural substrate, amino-acyl-tRNA, but was slowed down 200-fold with decreasing pH when puromycin was used as a substrate analog. Mutation of the conserved base A2451 of 23 S rRNA to U did not abolish the pH dependence of the reaction with puromycin in M. smegmatis, suggesting that A2451 did not confer the pH dependence. However, the A2451U mutation alters the structure of the peptidyl transferase center and changes the pattern of pH-dependent rearrangements, as probed by chemical modification of 23 S rRNA. A2451 seems to function as a pivot point in ordering the structure of the peptidyl transferase center rather than taking part in chemical catalysis.     

Structural insights of non-canonical U?U pair and Hoogsteen interaction probed with Se atom

Unlike DNA, in addition to the 2′-OH group, uracil nucleobase and its modifications play essential roles in structure and function diversities of non-coding RNAs. Non-canonical U•U base pair is ubiquitous in non-coding RNAs, which are highly diversified. However, it is not completely clear how uracil plays the diversifing roles. To investigate and compare the uracil in U-A and U•U base pairs, we have decided to probe them with a selenium atom by synthesizing the novel 4-Se-uridine (^SeU) phosphoramidite and Se-nucleobase-modified RNAs (^SeU-RNAs), where the exo-4-oxygen of uracil is replaced by selenium. Our crystal structure studies of U-A and U•U pairs reveal that the native and Se-derivatized structures are virtually identical, and both U-A and U•U pairs can accommodate large Se atoms. Our thermostability and crystal structure studies indicate that the weakened H-bonding in U-A pair may be compensated by the base stacking, and that the stacking of the trans-Hoogsteen U•U pairs may stabilize RNA duplex and its junction. Our result confirms that the hydrogen bond (O4^…H-C5) of the Hoogsteen pair is weak. Using the Se atom probe, our Se-functionalization studies reveal more insights into the U•U interaction and U-participation in structure and function diversification of nucleic acids.     

The critical role of the universally conserved A2602 of 23S ribosomal RNA in the release of the nascent peptide during translation termination

The ribosomal peptidyl transferase center is responsible for two fundamental reactions, peptide bond formation and nascent peptide release, during the elongation and termination phases of protein synthesis, respectively. We used in vitro genetics to investigate the functional importance of conserved 23S rRNA nucleotides located in the peptidyl transferase active site for transpeptidation and peptidyl-tRNA hydrolysis. While mutations at A2451, U2585, and C2063 (E. coli numbering) did not significantly affect either of the reactions, substitution of A2602 with C or its deletion abolished the ribosome ability to promote peptide release but had little effect on transpeptidation. This indicates that the mechanism of peptide release is distinct from that of peptide bond formation, with A2602 playing a critical role in peptide release during translation termination.     

Non-enzymic binding of formylmethionyl-transfer RNAf to Artemia salina ribosomes

40 S ribosomal subunits of Artemia salina embryos can bind formylmethionyl-transfer RNA_f non-enzymically, i.e., in the absence of initiation factors. This, like the enzymic reaction, is largely AUG-dependent. Much more fMet-tRNA_f is bound by 80 S ribosomes but, in this case, a large fraction (about two-thirds) of the binding is AUG-independent. Whereas the AUG-dependent binding is very sensitive to edeine, a potent initiation inhibitor, the AUG-independent binding is resistant to this antibiotic. Virtually all of the bound fMet-tRNA_f is in all cases capable of reacting with puromycin to form fMet-puromycin; hence the bound aminoaoyl-tRNA is in the peptidyl (donor) site of the 80 S ribosome. Non-acylated tRNAs also bind to this site with high affinity in a codon-independent reaction and block the 80 S binding of fMet-tRNA_f. The properties of the peptidyl site are consistent with a non-decoding site which harbors the initiator aminoacyl-tRNA, when the 80 S initiation complex is formed, and to which every molecule of tRNA remains temporarily attached following peptide bond synthesis.