Decarboxylation Mechanism of iso-Orotate Decarboxylase Revisited.

iso-Orotate decarboxylase (IDCase), which is involved in the thymidine salvage pathway, has attracted considerable interest owing to its chemical similarity to a hypothetical DNA decarboxylase in mammals. Although valuable insights into the active DNA demethylation of 5-methyl-cytosine can be obtained from the decarboxylation mechanism of 5-carboxyl-uracil (5caU) catalyzed by IDCase, this mechanism remains under debate. In this study, the catalytic mechanism of 5caU decarboxylation by IDCase was studied using hybrid quantum mechanics/molecular mechanics (QM/MM) methodologies and density functional theory (DFT) calculations with a truncated model. The calculations supported a mechanism involving three sequential stages: activation of the 5caU substrate via proton transfer from an arginine (R262') to the carboxyl group of 5caU, formation of a tetrahedral intermediate, and decarboxylation of the tetrahedral intermediate to generate uracil as the product. The reaction pathways and structures obtained using the QM/MM and DFT methods coincided with each other. These simulations provided detailed insights into the unique mechanism of IDCase, clarifying various unresolved issues, such as the critical role of R262'. In addition, aspartate D323 was found to act as a general base in the tetrahedral intermediate formation step and a general acid in the later C-C bond cleavage step.

Structural and functional properties of uridine 5'-monophosphate synthase from Coffea arabica.

In higher eukaryotes and plants, the last two sequential steps in the de novo biosynthesis of uridine 5'-monophosphate (UMP) are catalyzed by a bifunctional natural chimeric protein called UMP synthase (UMPS). In higher plants, UMPS consists of two naturally fused enzymes: orotate phosphoribosyltransferase (OPRTase) at N-terminal and orotidine-5'-monophosphate decarboxylase (ODCase) at C-terminal. In this work, we obtained the full functional recombinant protein UMPS from Coffea arabica (CaUMPS) and studied its structure-function relationships. A biochemical and structural characterization of a plant UMPS with its two functional domains is described together with the presentation of the first crystal structure of a plant ODCase at 1.4 Å resolution. The kinetic parameters measured of CaOPRTase and CaODCase domains were comparable to those reported. The crystallographic structure revealed that CaODCase is a dimer that conserves the typical fold observed in other ODCases from prokaryote and eukaryote with a 1-deoxy-ribofuranose-5'-phosphate molecule bound in the active site of one subunit induced a closed conformation. Our results add to the knowledge of one of the key enzymes of the de novo biosynthesis of pyrimidines in plant metabolism and open the door to future applications.

Human uridine 5'-monophosphate synthase stores metabolic potential in inactive biomolecular condensates.

Human uridine 5'-monophosphate synthase (HsUMPS) is a bifunctional enzyme that catalyzes the final two steps in de novo pyrimidine biosynthesis. The individual orotate phosphoribosyl transferase and orotidine monophosphate domains have been well characterized, but little is known about the overall structure of the protein and how the organization of domains impacts function. Using a combination of chromatography, electron microscopy, and complementary biophysical methods, we report herein that HsUMPS can be observed in two structurally distinct states, an enzymatically active dimeric form and a nonactive multimeric form. These two states readily interconvert to reach an equilibrium that is sensitive to perturbations of the active site and the presence of substrate. We determined that the smaller molecular weight form of HsUMPS is an S-shaped dimer that can self-assemble into relatively well-ordered globular condensates. Our analysis suggests that the transition between dimer and multimer is driven primarily by oligomerization of the orotate phosphoribosyl transferase domain. While the cellular distribution of HsUMPS is unaffected, quantification by mass spectrometry revealed that de novo pyrimidine biosynthesis is dysregulated when this protein is unable to assemble into inactive condensates. Taken together, our data suggest that HsUMPS self-assembles into biomolecular condensates as a means to store metabolic potential for the regulation of metabolic rates.

Analyzing the Reaction of Orotidine 5'-Phosphate Decarboxylase as a Way to Examine Some Key Catalytic Proposals.

This study analyzes the origin of enzyme catalysis by focusing on the reaction of orotidine 5'-phosphate decarboxylase (ODCase). This reaction involves an enormous catalytic effect of 23 kcal/mol that has been attributed to reactant state destabilization associated with the use of binding energy through the so-called Circe effect. However, our early studies and subsequent key experiments have shown that the presumed effect of the binding energy (namely, the strain exerted by a bond to a phosphate group) does not contribute to the catalysis. In this study, we perform quantitative empirical valence bond calculations that reproduce the catalytic effect of ODCase and the effect of removing the phosphate side chain. The calculations demonstrate that the effect of the phosphate is due to a change in reorganization energy and should not be described as an induced fit effect. Similarly, we show that the overall catalytic effect is due to electrostatic transition state stabilization, which again reflects the smaller reorganization energy in the enzyme than in water. We also elaborate on the problems with the induced fit proposal, including the fact that it does not serve to tell us what the actual origin of the action of the catalytic effect is. In addition to the above points, we use this paper to discuss misconceptions about the meaning of the preorganization effect, as well as other misunderstandings of what is being done in consistent calculations of enzyme catalysis.

Adenylate Kinase-Catalyzed Reactions of AMP in Pieces: Specificity for Catalysis at the Nucleoside Activator and Dianion Catalytic Sites.

The pressure to optimize the enzymatic rate acceleration for adenylate kinase (AK)-catalyzed phosphoryl transfer has led to the evolution of an induced-fit mechanism, where the binding energy from interactions between the protein and substrate adenosyl group is utilized to drive a protein conformational change that activates the enzyme for catalysis. The adenine group of adenosine contributes 11.8 kcal mol-1 to the total ≥14.7 kcal mol-1 adenosine stabilization of the transition state for AK-catalyzed phosphoryl transfer to AMP. The relative third-order rate constants for activation of adenylate kinase, by the C-5 truncated adenosine 1-(ß-d-erythrofuranosyl)adenine (EA), for catalysis of phosphoryl transfer from ATP to phosphite dianion (HP, kcat/KHPKAct = 260 M-2 s-1), fluorophosphate (47 M-2 s-1), and phosphate (9.6 M-2 s-1), show that substitution of -F for -H and of -OH for -H at HP results, respectively, in decreases in the reactivity of AK for catalysis of phosphoryl transfer due to polar and steric effects of the -F and -OH substituents. The addition of a 5'-CH2OH to the EA activator results in a 3.0 kcal mol-1 destabilization of the transition state for AK-activated phosphoryl transfer to HP due to a steric effect. This is smaller than the 8.3 kcal mol-1 steric effect of the 5'-CH2OH substituent at OMP on HP-activated OMPDC-catalyzed decarboxylation of 1-(ß-d-erythrofuranosyl)orotate. The 2'-OH ribosyl substituent shows significant interactions with the transition states for AK-catalyzed phosphoryl transfer from ATP to AMP and for adenosine-activated AK-catalyzed phosphoryl transfer from ATP to HP.

Uricase sensitizes hepatocellular carcinoma cells to 5-fluorouracil through uricase-uric acid-UMP synthase axis

Conventional chemotherapy plays a key role in hepatocellular carcinoma (HCC) treatment, however, with intrinsic or acquired chemoresistance being a major constraint. Here, we aimed to identify potential target to reverse such chemoresistance. In the present study, we found significant difference in uridine monophosphate synthetase (UMPS) expression between 5-FU resistant and sensitive HCC cell lines and the overexpression or downregulation of UMPS impacted 5-FU response in HCC cells. We further found that inhibition of UMPS activity with uric acid at concentration present in human plasma decreased the 5-FU sensitivity of HCC cells, while reduction of uric acid levels with uricase improved the 5-FU sensitivity of HCC cells as well as colorectal cancer cells. In vivo studies also suggested that modulation of uric acid levels did affect 5-FU sensitivity of tumors. These data indicated that UMPS was correlated with the 5-FU resistance in HCC cells and uricase sensitized cancer cells to 5-FU through uricase-uric acid-UMP synthase axis, which provided a potential strategy to improve the efficacy of 5-FU-based chemotherapy for human cancers. (AU)

Enabling Role of Ligand-Driven Conformational Changes in Enzyme Evolution.

Many enzymes that show a large specificity in binding the enzymatic transition state with a higher affinity than the substrate utilize substrate binding energy to drive protein conformational changes to form caged substrate complexes. These protein cages provide strong stabilization of enzymatic transition states. Using part of the substrate binding energy to drive the protein conformational change avoids a similar strong stabilization of the Michaelis complex and irreversible ligand binding. A seminal step in the development of modern enzyme catalysts was the evolution of enzymes that couple substrate binding to a conformational change. These include enzymes that function in glycolysis (triosephosphate isomerase), the biosynthesis of lipids (glycerol phosphate dehydrogenase), the hexose monophosphate shunt (6-phosphogluconate dehydrogenase), and the mevalonate pathway (isopentenyl diphosphate isomerase), catalyze the final step in the biosynthesis of pyrimidine nucleotides (orotidine monophosphate decarboxylase), and regulate the cellular levels of adenine nucleotides (adenylate kinase). The evolution of enzymes that undergo ligand-driven conformational changes to form active protein-substrate cages is proposed to proceed by selection of variants, in which the selected side chain substitutions destabilize a second protein conformer that shows compensating enhanced binding interactions with the substrate. The advantages inherent to enzymes that incorporate a conformational change into the catalytic cycle provide a strong driving force for the evolution of flexible protein folds such as the TIM barrel. The appearance of these folds represented a watershed event in enzyme evolution that enabled the rapid propagation of enzyme activities within enzyme superfamilies.

Uricase sensitizes hepatocellular carcinoma cells to 5-fluorouracil through uricase-uric acid-UMP synthase axis.

Conventional chemotherapy plays a key role in hepatocellular carcinoma (HCC) treatment, however, with intrinsic or acquired chemoresistance being a major constraint. Here, we aimed to identify potential target to reverse such chemoresistance. In the present study, we found significant difference in uridine monophosphate synthetase (UMPS) expression between 5-FU resistant and sensitive HCC cell lines and the overexpression or downregulation of UMPS impacted 5-FU response in HCC cells. We further found that inhibition of UMPS activity with uric acid at concentration present in human plasma decreased the 5-FU sensitivity of HCC cells, while reduction of uric acid levels with uricase improved the 5-FU sensitivity of HCC cells as well as colorectal cancer cells. In vivo studies also suggested that modulation of uric acid levels did affect 5-FU sensitivity of tumors. These data indicated that UMPS was correlated with the 5-FU resistance in HCC cells and uricase sensitized cancer cells to 5-FU through uricase-uric acid-UMP synthase axis, which provided a potential strategy to improve the efficacy of 5-FU-based chemotherapy for human cancers.

Protein-Ribofuranosyl Interactions Activate Orotidine 5'-Monophosphate Decarboxylase for Catalysis.

The role of a global, substrate-driven, enzyme conformational change in enabling the extraordinarily large rate acceleration for orotidine 5'-monophosphate decarboxylase (OMPDC)-catalyzed decarboxylation of orotidine 5'-monophosphate (OMP) is examined in experiments that focus on the interactions between OMPDC and the ribosyl hydroxyl groups of OMP. The D37 and T100' side chains of OMPDC interact, respectively, with the C-3' and C-2' hydroxyl groups of enzyme-bound OMP. D37G and T100'A substitutions result in 1.4 kcal/mol increases in the activation barrier ΔG⧧ for catalysis of decarboxylation of the phosphodianion-truncated substrate 1-(ß-d-erythrofuranosyl)orotic acid (EO) but result in larger 2.1-2.9 kcal/mol increases in ΔG⧧ for decarboxylation of OMP and for phosphite dianion-activated decarboxylation of EO. This shows that these substitutions reduce transition-state stabilization by the Q215, Y217, and R235 side chains at the dianion binding site. The D37G and T100'A substitutions result in <1.0 kcal/mol increases in ΔG⧧ for activation of OMPDC-catalyzed decarboxylation of the phosphoribofuranosyl-truncated substrate FO by phosphite dianions. Experiments to probe the effect of D37 and T100' substitutions on the kinetic parameters for d-glycerol 3-phosphate and d-erythritol 4-phosphate activators of OMPDC-catalyzed decarboxylation of FO show that ΔG⧧ for sugar phosphate-activated reactions is increased by ca. 2.5 kcal/mol for each -OH interaction eliminated by D37G or T100'A substitutions. We conclude that the interactions between the D37 and T100' side chains and ribosyl or ribosyl-like hydroxyl groups are utilized to activate OMPDC for catalysis of decarboxylation of OMP, EO, and FO.

A case-control study of a combination of single nucleotide polymorphisms and clinical parameters to predict clinically relevant toxicity associated with fluoropyrimidine and platinum-based chemotherapy in gastric cancer.

BACKGROUND: Fluoropyrimidine plus platinum chemotherapy remains the standard first line treatment for gastric cancer (GC). Guidelines exist for the clinical interpretation of four DPYD genotypes related to severe fluoropyrimidine toxicity within European populations. However, the frequency of these single nucleotide polymorphisms (SNPs) in the Latin American population is low (< 0.7%). No guidelines have been development for platinum. Herein, we present association between clinical factors and common SNPs in the development of grade 3-4 toxicity. METHODS: Retrospectively, 224 clinical records of GC patient were screened, of which 93 patients were incorporated into the study. Eleven SNPs with minor allelic frequency above 5% in GSTP1, ERCC2, ERCC1, TP53, UMPS, SHMT1, MTHFR, ABCC2 and DPYD were assessed. Association between patient clinical characteristics and toxicity was estimated using logistic regression models and classification algorithms. RESULTS: Reported grade ≤ 2 and 3-4 toxicities were 64.6% (61/93) and 34.4% (32/93) respectively. Selected DPYD SNPs were associated with higher toxicity (rs1801265; OR = 4.20; 95% CI = 1.70-10.95, p = 0.002), while others displayed a trend towards lower toxicity (rs1801159; OR = 0.45; 95% CI = 0.19-1.08; p = 0.071). Combination of paired SNPs demonstrated significant associations in DPYD (rs1801265), UMPS (rs1801019), ABCC2 (rs717620) and SHMT1 (rs1979277). Using multivariate logistic regression that combined age, sex, peri-operative chemotherapy, 5-FU regimen, the binary combination of the SNPs DPYD (rs1801265) + ABCC2 (rs717620), and DPYD (rs1801159) displayed the best predictive performance. A nomogram was constructed to assess the risk of developing overall toxicity. CONCLUSION: Pending further validation, this model could predict chemotherapy associated toxicity and improve GC patient quality of life.

Novel mutation in UMPS gene leads to false positive result of DUMPS (genetic disorder) in buffaloes: need for gene sequencing before confirming results of RFLP in new species.

Deficiency of uridine monophosphate synthase (DUMPS) is a lethal genetic disorder associated with early embryonic mortality. Murrah and Mehsana male buffaloes (n = 594) were screened for DUMPS by PCR-RFLP technique. A few Murrah buffalo male calves were found to be carriers of DUMPS in RFLP, which has not been reported earlier. On the Sanger sequencing, a novel A to G substitution mutation was identified in AvaI restriction recognition site of UMPS gene in buffaloes. This mutation hinders digestion of DNA by AvaI which leds to false positive results for DUMPS carrier in RFLP. The results indicated that genome sequencing must be performed before confirming results of RFLP in any new species. All the buffaloes that were tested had only wild-type genotype in exon 5 for DUMPS specific allele.

Development of a pyrF-based counterselectable system for targeted gene deletion in Streptomyces rimosus.

Streptomyces produces many valuable and important biomolecules with clinical and pharmaceutical applications. The development of simple and highly efficient gene editing tools for genetic modification of Streptomyces is highly desirable. In this study, we developed a screening system for targeted gene knockout using a uracil auxotrophic host (ΔpyrF) resistant to the highly toxic uracil analog of 5-fluoroorotic acid (5-FOA) converted by PyrF, and a non-replicative vector pKC1132-pyrF carrying the complemented pyrF gene coding for orotidine-5'-phosphate decarboxylase. The pyrF gene acts as a positive selection and counterselection marker for recombinants during genetic modifications. Single-crossover homologous integration mutants were selected on minimal medium without uracil by reintroducing pyrF along with pKC1132-pyrF into the genome of the mutant ΔpyrF at the targeted locus. Double-crossover recombinants were generated, from which the pyrF gene, plasmid backbone, and targeted gene were excised through homologous recombination exchange. These recombinants were rapidly screened by the counterselection agent, 5-FOA. We demonstrated the feasibility and advantage of using this pyrF-based screening system through deleting the otcR gene, which encodes the cluster-situated regulator that directly activates oxytetracycline biosynthesis in Streptomyces rimosus M4018. This system provides a new genetic tool for investigating the genetic characteristics of Streptomyces species.

Identification of critical genetic variants associated with metabolic phenotypes of the Japanese population.

We performed a metabolome genome-wide association study for the Japanese population in the prospective cohort study of Tohoku Medical Megabank. By combining whole-genome sequencing and nontarget metabolome analyses, we identified a large number of novel associations between genetic variants and plasma metabolites. Of the identified metabolite-associated genes, approximately half have already been shown to be involved in various diseases. We identified metabolite-associated genes involved in the metabolism of xenobiotics, some of which are from intestinal microorganisms, indicating that the identified genetic variants also markedly influence the interaction between the host and symbiotic bacteria. We also identified five associations that appeared to be female-specific. A number of rare variants that influence metabolite levels were also found, and combinations of common and rare variants influenced the metabolite levels more profoundly. These results support our contention that metabolic phenotyping provides important insights into how genetic and environmental factors provoke human diseases.

Inhibition of the de novo pyrimidine biosynthesis pathway limits ribosomal RNA transcription causing nucleolar stress in glioblastoma cells.

Glioblastoma is the most common and aggressive type of cancer in the brain; its poor prognosis is often marked by reoccurrence due to resistance to the chemotherapeutic agent temozolomide, which is triggered by an increase in the expression of DNA repair enzymes such as MGMT. The poor prognosis and limited therapeutic options led to studies targeted at understanding specific vulnerabilities of glioblastoma cells. Metabolic adaptations leading to increased synthesis of nucleotides by de novo biosynthesis pathways are emerging as key alterations driving glioblastoma growth. In this study, we show that enzymes necessary for the de novo biosynthesis of pyrimidines, DHODH and UMPS, are elevated in high grade gliomas and in glioblastoma cell lines. We demonstrate that DHODH's activity is necessary to maintain ribosomal DNA transcription (rDNA). Pharmacological inhibition of DHODH with the specific inhibitors brequinar or ML390 effectively depleted the pool of pyrimidines in glioblastoma cells grown in vitro and in vivo and impaired rDNA transcription, leading to nucleolar stress. Nucleolar stress was visualized by the aberrant redistribution of the transcription factor UBF and the nucleolar organizer nucleophosmin 1 (NPM1), as well as the stabilization of the transcription factor p53. Moreover, DHODH inhibition decreased the proliferation of glioblastoma cells, including temozolomide-resistant cells. Importantly, the addition of exogenous uridine, which reconstitutes the cellular pool of pyrimidine by the salvage pathway, to the culture media recovered the impaired rDNA transcription, nucleolar morphology, p53 levels, and proliferation of glioblastoma cells caused by the DHODH inhibitors. Our in vivo data indicate that while inhibition of DHODH caused a dramatic reduction in pyrimidines in tumor cells, it did not affect the overall pyrimidine levels in normal brain and liver tissues, suggesting that pyrimidine production by the salvage pathway may play an important role in maintaining these nucleotides in normal cells. Our study demonstrates that glioblastoma cells heavily rely on the de novo pyrimidine biosynthesis pathway to generate ribosomal RNA (rRNA) and thus, we identified an approach to inhibit ribosome production and consequently the proliferation of glioblastoma cells through the specific inhibition of the de novo pyrimidine biosynthesis pathway.

Generalized correlation-based dynamical network analysis: a new high-performance approach for identifying allosteric communications in molecular dynamics trajectories.

Molecular interactions are essential for regulation of cellular processes from the formation of multi-protein complexes to the allosteric activation of enzymes. Identifying the essential residues and molecular features that regulate such interactions is paramount for understanding the biochemical process in question, allowing for suppression of a reaction through drug interventions or optimization of a chemical process using bioengineered molecules. In order to identify important residues and information pathways within molecular complexes, the dynamical network analysis method was developed and has since been broadly applied in the literature. However, in the dawn of exascale computing, this method is frequently limited to relatively small biomolecular systems. In this work, we provide an evolution of the method, application, and interface. All data processing and analysis are conducted through Jupyter notebooks, providing automatic detection of important solvent and ion residues, an optimized and parallel generalized correlation implementation that is linear with respect to the number of nodes in the system, and subsequent community clustering, calculation of betweenness of contacts, and determination of optimal paths. Using the popular visualization program visual molecular dynamics (VMD), high-quality renderings of the networks over the biomolecular structures can be produced. Our new implementation was employed to investigate three different systems, with up to 2.5M atoms, namely, the OMP-decarboxylase, the leucyl-tRNA synthetase complexed with its cognate tRNA and adenylate, and respiratory complex I in a membrane environment. Our enhanced and updated protocol provides the community with an intuitive and interactive interface, which can be easily applied to large macromolecular complexes.

The development of an efficient RNAi system based on Agrobacterium-mediated transformation approach for studying functional genomics in medical fungus Wolfiporia cocos.

Genetic transformation methods reported for Wolfiporia cocos are limited. In this study, we describe an efficient RNA interference (RNAi) system based on Agrobacterium-mediated transformation approach in W. cocos for the first time. Actively growing mycelial plugs were used as recipients for transformation using endogenous orotidine-5'-phosphate decarboxylase gene (URA3) as both a selective marker and a silencing gene, under the control of the dual promoters of Legpd and Leactin from Lentinula edodes and the single promoter of Wcgpd from W. cocos, respectively. The results showed that both the two kinds of promoters effectively drive the expression of URA3 gene, and the URA3-silenced transformants could be selected on CYM medium containing 5'-fluoroorotic acid. In addition, silencing URA3 gene has no effect on the growth of W. cocos hyphae. The incomplete silencing of the URA3 locus was also observed in this study. This study will promote further study on the mechanism of substrate degradation, sclerotial formation, and biosynthesis network of pharmacological compounds in W. cocos.

Metabolic engineering generates a transgene-free safety switch for cell therapy.

Safeguard mechanisms can ameliorate the potential risks associated with cell therapies but currently rely on the introduction of transgenes. This limits their application owing to immunogenicity or transgene silencing. We aimed to create a control mechanism for human cells that is not mediated by a transgene. Using genome editing methods, we disrupt uridine monophosphate synthetase (UMPS) in the pyrimidine de novo synthesis pathway in cell lines, pluripotent cells and primary human T cells. We show that this makes proliferation dependent on external uridine and enables us to control cell growth by modulating the uridine supply, both in vitro and in vivo after transplantation in xenograft models. Additionally, disrupting this pathway creates resistance to 5-fluoroorotic acid, which enables positive selection of UMPS-knockout cells. We envision that this approach will add an additional level of safety to cell therapies and therefore enable the development of approaches with higher risks, especially those that are intended for limited treatment durations.

An efficient transformation system for Trichoderma atroviride using the pyr4 gene as a selectable marker.

The development of an efficient transformation system is essential to enrich the genetic understanding of Trichoderma atroviride. To acquire an additional homologous selectable marker, uracil auxotrophic mutants were generated. First, the pyr4 gene encoding OMP decarboxylase was replaced by the hph marker gene, encoding a hygromycin phosphotransferase. Then, uracil auxotrophs were employed to determine that 5 mM uracil restores their growth and conidia production, and 1 mg ml-1 is the lethal dose of 5-fluoroorotic acid in T. atroviride. Subsequently, uracil auxotrophic strains, free of a drug-selectable marker, were selected by 5-fluoroorotic acid resistance. Two different deletions in pyr4 were mapped in four auxotrophs, encoding a protein with frameshifts at the 310 and 335 amino acids in their COOH-terminal. Six auxotrophs did not have changes in the pyr4 ORF even though a specific cassette to delete the pyr4 was used, suggesting that 5-FOA could have mutagenic activity. The Ura-1 strain was selected as a genetic background to knock out the MAPKK Pbs2, MAPK Tmk3, and the blue light receptors Blr1/Blr2, using a short version of pyr4 as a homologous marker. The ∆tmk3 and ∆pbs2 mutants selected with pyr4 or hph marker were phenotypically identical, highly sensitive to different stressors, and affected in photoconidiation. The ∆blr1 and ∆blr2 mutants were not responsive to light, and complementation of uracil biosynthesis did not interfere in the expression of blu1, grg2, phr1, and env1 genes upregulated by blue light. Overall, uracil metabolism can be used as a tool for genetic manipulation in T. atroviride.

Direct Determination of Pseudouridine in RNA by Mass Spectrometry Coupled with Stable Isotope Labeling.

Pseudouridine (Ψ) is the only "mass-silent" nucleoside produced by post-transcriptional RNA modification. We developed a mass spectrometry (MS)-based technique coupled with in vivo deuterium (D) labeling of uridines for direct determination of Ψs in cellular RNA and applied it to the comprehensive analysis of post-transcriptional modifications in human ribosomal RNAs. The method utilizes human TK6/mouse FM3A cells deficient in uridine monophosphate synthase using a CRISPR-Cas9 technique to turn off de novo uridine synthesis and fully labels uridines with D at uracil positions 5 and 6 by cultivating the cells in a medium containing uridine-5,6-D2. The pseudouridylation reaction in those cells results in the exchange of the D at the C5 of uracil with hydrogen from solvent, which produces a -1 Da mass shift, thus allowing MS-based determination of RNA Ψs. We present here the experimental details of this method and show that it allows the identification of all Ψs in human major nuclear and nucleolar RNAs, including several previously unknown Ψs. Because the method allows direct determination of Ψs at the femtomole level of RNA, it will serve as a useful tool for structure/function studies of a wide variety of noncoding RNAs.