Cytidine Acetylation Across the Tree of Life.

Acetylation plays a critical role in regulating eukaryotic transcription via the modification of histones. Beyond this well-documented function, a less explored biological frontier is the potential for acetylation to modify and regulate the function of RNA molecules themselves. N4-Acetylcytdine (ac4C) is a minor RNA nucleobase conserved across all three domains of life (archaea, bacteria, and eukarya), a conservation that suggests a fundamental role in biological processes. Unlike many RNA modifications that are controlled by large enzyme families, almost all organisms catalyze ac4C using a homologue of human Nat10, an essential disease-associated acetyltransferase enzyme.A critical step in defining the fundamental functions of RNA modifications has been the development of methods for their sensitive and specific detection. This Account describes recent progress enabling the use of chemical sequencing reactions to map and quantify ac4C with single-nucleotide resolution in RNA. To orient readers, we first provide historical background of the discovery of ac4C and the enzymes that catalyze its formation. Next, we describe mechanistic experiments that led to the development of first- and second-generation sequencing reactions able to determine ac4C's position in a polynucleotide by exploiting the nucleobase's selective susceptibility to reduction by hydride donors. A notable feature of this chemistry, which may serve as a prototype for nucleotide resolution RNA modification sequencing reactions more broadly, is its ability to drive a penetrant and detectable gain of signal specifically at ac4C sites. Emphasizing practical applications, we present how this optimized chemistry can be integrated into experimental workflows capable of sensitive, transcriptome-wide analysis. Such readouts can be applied to quantitatively define the ac4C landscape across the tree of life. For example, in human cell lines and yeast, this method has uncovered that ac4C is highly selective, predominantly occupying dominant sites within rRNA (rRNA) and tRNA (tRNA). By contrast, when we extend these analyses to thermophilic archaea they identify the potential for much more prevalent patterns of cytidine acetylation, leading to the discovery of a role for this modification in adaptation to environmental stress. Nucleotide resolution analyses of ac4C have also allowed for the determination of structure-activity relationships required for short nucleolar RNA (snoRNA)-catalyzed ac4C deposition and the discovery of organisms with unexpectedly divergent tRNA and rRNA acetylation signatures. Finally, we share how these studies have shaped our approach to evaluating novel ac4C sites reported in the literature and highlight unanswered questions and new directions that set the stage for future research in the field.

Multiplexed LC-MS/MS quantification of salivary RNA modifications in periodontitis.

OBJECTIVE: To analyse the salivary epitranscriptomic profiles as periodontitis biomarkers using multiplexed mass spectrometry (MS). BACKGROUND: The field of epitranscriptomics, which relates to RNA chemical modifications, opens new perspectives in the discovery of diagnostic biomarkers, especially in periodontitis. Recently, the modified ribonucleoside N6-methyladenosine (m6A) was revealed as a crucial player in the etiopathogenesis of periodontitis. However, no epitranscriptomic biomarker has been identified in saliva to date. MATERIALS AND METHODS: Twenty-four saliva samples were collected from periodontitis patients (n = 16) and from control subjects (n = 8). Periodontitis patients were stratified according to stage and grade. Salivary nucleosides were directly extracted and, in parallel, salivary RNA was digested into its constituent nucleosides. Nucleoside samples were then quantified by multiplexed MS. RESULTS: Twenty-seven free nucleosides were detected and an overlapping set of 12 nucleotides were detected in digested RNA. Among the free nucleosides, cytidine and three other modified nucleosides (inosine, queuosine and m6Am) were significantly altered in periodontitis patients. In digested RNA, only uridine was significantly higher in periodontitis patients. Importantly there was no correlation between free salivary nucleoside levels and the levels of those same nucleotides in digested salivary RNA, except for cytidine, m5C and uridine. This statement implies that the two detection methods are complementary. CONCLUSION: The high specificity and sensitivity of MS allowed the detection and quantification of multiple nucleosides from RNA and free nucleosides in saliva. Some ribonucleosides appear to be promising biomarkers of periodontitis. Our analytic pipeline opens new perspectives for diagnostic periodontitis biomarkers.

NAT10: An RNA cytidine transferase regulates fatty acid metabolism in cancer cells.

BACKGROUND: N-4 cytidine acetylation (ac4C) is an epitranscriptomics modification catalyzed by N-acetyltransferase 10 (NAT10); important for cellular mRNA stability, rRNA biogenesis, cell proliferation and epithelial to mesenchymal transition (EMT). However, whether other crucial pathways are regulated by NAT10-dependent ac4C modification in cancer cells remains unclear. Therefore, in this study, we explored the impact of NAT10 depletion in cancer cells using unbiased RNA-seq. METHODS: High-throughput sequencing of knockdown NAT10 in cancer cells was conducted to identify enriched pathways. Acetylated RNA immunoprecipitation-seq (acRIP-seq) and RIP-PCR were used to map and determine ac4C levels of RNA. Exogenous palmitate uptake assay was conducted to assess NAT10 knockdown cancer cells using Oil Red O staining and lipid content analysis. Gas-chromatography-tandem mass spectroscopy (GC/MS) was used to perform untargeted lipidomics. RESULTS: High-throughput sequencing of NAT10 knockdown in cancer cells revealed fatty acid (FA) metabolism as the top enriched pathway through the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis in differentially downregulated genes. FA metabolic genes such as ELOLV6, ACSL1, ACSL3, ACSL4, ACADSB and ACAT1 were shown to be stabilised via NAT10-dependent ac4C RNA acetylation. Additionally, NAT10 depletion was shown to significantly reduce the levels of overall lipid content, triglycerides and total cholesterol. Further, NAT10 depletion in palmitate-loaded cancer cells showed decrease in ac4C levels across the RNA transcripts of FA metabolic genes. In untargeted lipidomics, 496 out of 2 279 lipids were statistically significant in NAT10 depleted cancer cells, of which pathways associated with FA metabolism are the most enriched. CONCLUSIONS: Conclusively, our results provide novel insights into the impact of NAT10-mediated ac4C modification as a crucial regulatory factor during FA metabolism and showed the benefit of targeting NAT10 for cancer treatment.

Nucleotide resolution profiling of m3C RNA modification by HAC-seq.

Cellular RNAs are subject to a myriad of different chemical modifications that play important roles in controlling RNA expression and function. Dysregulation of certain RNA modifications, the so-called 'epitranscriptome', contributes to human disease. One limitation in studying the functional, physiological, and pathological roles of the epitranscriptome is the availability of methods for the precise mapping of individual RNA modifications throughout the transcriptome. 3-Methylcytidine (m3C) modification of certain tRNAs is well established and was also recently detected in mRNA. However, methods for the specific mapping of m3C throughout the transcriptome are lacking. Here, we developed a m3C-specific technique, Hydrazine-Aniline Cleavage sequencing (HAC-seq), to profile the m3C methylome at single-nucleotide resolution. We applied HAC-seq to analyze ribosomal RNA (rRNA)-depleted total RNAs in human cells. We found that tRNAs are the predominant m3C-modified RNA species, with 17 m3C modification sites on 11 cytoplasmic and 2 mitochondrial tRNA isoacceptors in MCF7 cells. We found no evidence for m3C-modification of mRNA or other non-coding RNAs at comparable levels to tRNAs in these cells. HAC-seq provides a novel method for the unbiased, transcriptome-wide identification of m3C RNA modification at single-nucleotide resolution, and could be widely applied to reveal the m3C methylome in different cells and tissues.

Visualizing RNA Cytidine Acetyltransferase Activity by Northern Blotting.

Cytidine acetyltransferases are an emerging class of nucleic-acid-modifying enzymes responsible for the establishment of N4 -acetylcytidine (ac4C) in RNA. In contrast to histone acetyltransferases, whose activity is commonly studied by western blotting, relatively few methods exist for quickly assessing the activity of cytidine acetyltransferases from a biological sample of interest or the distribution of ac4C across different RNA species. In this protocol, we describe a method for analysis of cellular cytidine acetyltransferase activity using dot- and immuno-northern-blotting-based detection. © 2020 Wiley Periodicals LLC. Basic Protocol 1: Detection of N4 -Acetylcytidine in RNA by dot blotting Basic Protocol 2: Visualizing N4 -Acetylcytidine Distribution in RNA by northern blotting.

Chemical labeling - Assisted mass spectrometry analysis for sensitive detection of cytidine dual modifications in RNA of mammals.

RNA molecules carry diverse modifications that exert important influences in many cellular processes. In addition to the single modification occurring in either nucleobase or 2' hydroxyl of ribose in RNA, some dual modifications occur in both the nucleobase and 2' hydroxyl of ribose in RNA. 2'-O-methyl-5-methylcytidine (m5Cm), the dual modifications of cytidine, was first discovered from the tRNA of archaea. Recent studies identified that 2'-O-methyl-5-hydroxymethylcytidine (hm5Cm) and 2'-O-methyl-5-formylcytidine (f5Cm) were present in the anticodon of cytoplasmic tRNA of mammals. Similar to the series of single modification of cytidines of 5-methylcytosine (m5C), 5-hydroxymethylcytidine (hm5C), 5-formylcytidine (f5C), and 5-carboxylcytidine (ca5C) in nucleic acids, the dual modifications of m5Cm, hm5Cm, f5Cm and 2'-O-methyl-5-carboxylcytidine (ca5Cm) may also constitute the series of cytidine modifications in mammals. However, it is normally challenging to detect these modifications because of their low endogenous levels. Here, we established a method by chemical labeling-assisted liquid chromatography - electrospray ionization - tandem mass spectrometry (LC-ESI-MS/MS) analysis for the sensitive and simultaneous determination of all these four cytidine dual modifications, i.e., m5Cm, hm5Cm, f5Cm and ca5Cm. Three different labeling reagents (2-bromo-1-(3,4-dimeth oxyphenyl)-ethanone, BDMOPE; 2-bromo-1-(4-methoxyphenyl)-ethanone, BMOPE; 2-bromo-1-(4-diethylaminophenyl)-ethanone, BDEPE) were used for the chemical labeling. The results showed that the detection sensitivities of m5Cm, hm5Cm, f5Cm and ca5Cm increased up to 462 folds after chemical labeling. With the developed method, we achieved the simultaneous detection of m5Cm, hm5Cm and f5Cm in RNA of mammals. In addition, we found these cytidine dual modifications mainly exist in small RNA (<200 nt) and barely detected in other types of RNA. Moreover, we found that the levels of m5Cm in RNA of human lung carcinoma tissues significantly increased, while hm5Cm and f5Cm significantly decreased compared to tumor adjacent normal tissues. The significant changes of m5Cm, hm5Cm and f5Cm levels may serve as indicator for the detection and prognosis of lung cancer.

Intensive optimization and evaluation of global DNA methylation quantification using LC-MS/MS.

DNA methylation is a typical epigenetic phenomenon. Numerous methods for detecting global DNA methylation levels have been developed, among which LC-MS/MS has emerged as an excellent method from the viewpoint of sensitivity, reproducibility, and cost. However, LC-MS/MS methods have limitations due to a lack of the stability and the standardization required for a laboratory assay. The present study aimed to establish a robust assay that guarantees highly accurate measurements of global DNA methylation levels. There are at least three facets of the developed method. The first is discovery of the solvent conditions to minimize sodium adducts. The second is improvement of separation of nucleosides by LC using the columns that had not been used in previous similar studies. The third is success in reduction of the uncertainty of the measurement results, which was achieved by the calibration using the ratio of mdC but not the absolute amount in the presence of internal standards. These facets represent the advantage over methods reported previously. Our developed method enables quantification of DNA methylation with a short time length (8 min) for one analysis as well as with the high reproducibility of measurements that is represented by the inter-day CV% being less than 5%. In addition, data obtained from measuring global DNA methylation levels in cultured cell lines, with or without pharmacological demethylation, support its use for biomedical research. This assay is expected to allow us to conduct initial screening of epigenetic alterations or aberration in a variety of cells.

Identification and Quantification of Modified Nucleosides in Saccharomyces cerevisiae mRNAs.

Post-transcriptional modifications to messenger RNAs (mRNAs) have the potential to alter the biological function of this important class of biomolecules. The study of mRNA modifications is a rapidly emerging field, and the full complement of chemical modifications in mRNAs is not yet established. We sought to identify and quantify the modifications present in yeast mRNAs using an ultra-high performance liquid chromatography tandem mass spectrometry method to detect 40 nucleoside variations in parallel. We observe six modified nucleosides with high confidence in highly purified mRNA samples (N7-methylguanosine, N6-methyladenosine, 2'-O-methylguanosine, 2'-O-methylcytidine, N4-acetylcytidine, and 5-formylcytidine) and identify the yeast protein responsible for N4-acetylcytidine incorporation in mRNAs (Rra1). In addition, we find that mRNA modification levels change in response to heat shock, glucose starvation, and/or oxidative stress. This work expands the repertoire of potential chemical modifications in mRNAs and highlights the value of integrating mass spectrometry tools in the mRNA modification discovery and characterization pipeline.

Determination of five nucleosides by LC-MS/MS and the application of the method to quantify N6 -methyladenosine level in liver messenger ribonucleic acid of an acetaminophen-induced hepatotoxicity mouse model.

Ribonucleic acid N6 -methyladenosine methylation plays an important role in a variety of biological processes and diseases. Acetaminophen-induced hepatotoxicity is one of the major challenges faced by clinicians. To date, the link between N6 -methyladenosine and acetaminophen-induced hepatotoxicity has not been studied. In this study, a simple ultra high performance liquid chromatography with tandem mass spectrometry method was developed for the simultaneous determination of five nucleosides (adenosine, uridine, cytidine, guanosine, and N6 -methyladenosine) in messenger ribonucleic acid. After enzymatic digestion of messenger ribonucleic acid, the nucleosides sample was separated on an Acquity UPLC column with gradient elution using methanol and 0.02% formic acid water, and detected by a Qtrap 4500 mass spectrometer with an electrospray ionization mode. The method was validated over the concentration ranges of 4-800 ng/mL for adenosine, uridine, cytidine, and guanosine and 0.1-20 ng/mL for N6 -methyladenosine. It was successfully applied to the determination of N6 -methyladenosine levels in liver messenger ribonucleic acid in an acetaminophen-induced hepatotoxicity mouse model and a control group. This study offers a method for the determination of nucleoside contents in epigenetic studies and constitutes the first step toward the investigation of ribonucleic acid methylation in acetaminophen-induced hepatotoxicity, which will facilitate the elucidation of its mechanism.

Nucleotide resolution sequencing of N4-acetylcytidine in RNA.

Posttranscriptional modifications of RNA represent an emerging class of regulatory elements in human biology. Improved methods for studying how these elements are controlled and where they occur has the potential to transform our understanding of gene expression in development and disease. Here we describe a chemical method for nucleotide resolution sequencing of N4-acetylcytidine (ac4C), a highly conserved modified nucleobase whose formation is catalyzed by the essential cytidine acetyltransferase enzyme NAT10. This approach enables the sensitive, PCR-amplifiable detection of individual ac4C sites from nanograms of unfractionated cellular RNA. The sensitive and quantitative nature of this assay provides a powerful tool to understand how cytidine acetylation is targeted, profile RNA acetyltransferase dynamics, and validate the sites and stoichiometry of ac4C in novel RNA species.

Determination of genomic N3-methylthymidine in human cancer cells treated with nitrosamines using capillary electrophoresis with laser-induced fluorescence.

Methylating substances alter DNA by forming N3-methylthymidine (N3mT), a mutagenic base modification. To develop a sensitive analytical method for the detection of N3mT in DNA based on capillary electrophoresis with laser-induced fluorescence detection (CE-LIF), we synthesized the N3mT-3'-phosphate as a chemical standard. The limit of detection was 1.9 amol of N3mT, which corresponds to one molecule of N3mT per 1000 normal nucleotides or 0.1%. With this method, we demonstrated that the carcinogenic nitrosamine N'-nitrosonornicotine (NNN) induced N3mT in the human lung cancer cell line A549. Treatment with NNN also caused an elevated degree of 5-hydroxymethylcytidine (5hmdC) in DNA, while the methylation degree (i.e. 5-methylcytidine; 5mdC) stayed constant. According to our data, NNN could, via yet unknown mechanisms, play a role in the formation of N3mT as well as 5hmdC. In this study we have developed a new sensitive analytical method using CE-LIF for the simultaneous detection of the three DNA modifications, 5mdC, 5hmdC and N3mT.

AlkAniline-Seq: Profiling of m7 G and m3 C RNA Modifications at Single Nucleotide Resolution.

RNA modifications play essential roles in gene expression regulation. Only seven out of >150 known RNA modifications are detectable transcriptome-wide by deep sequencing. Here we describe a new principle of RNAseq library preparation, which relies on a chemistry based positive enrichment of reads in the resulting libraries, and therefore leads to unprecedented signal-to-noise ratios. The proposed approach eschews conventional RNA sequencing chemistry and rather exploits the generation of abasic sites and subsequent aniline cleavage. The newly generated 5'-phosphates are used as unique entry for ligation of an adapter in library preparation. This positive selection, embodied in the AlkAniline-Seq, enables a deep sequencing-based technology for the simultaneous detection of 7-methylguanosine (m7 G) and 3-methylcytidine (m3 C) in RNA at single nucleotide resolution. As a proof-of-concept, we used AlkAniline-Seq to comprehensively validate known m7 G and m3 C sites in bacterial, yeast, and human cytoplasmic and mitochondrial tRNAs and rRNAs, as well as for identifying previously unmapped positions.

A Chemical Signature for Cytidine Acetylation in RNA.

N4-acetylcytidine (ac4C) is a highly conserved modified RNA nucleobase whose formation is catalyzed by the disease-associated N-acetyltransferase 10 (NAT10). Here we report a sensitive chemical method to localize ac4C in RNA. Specifically, we characterize the susceptibility of ac4C to borohydride-based reduction and show this reaction can cause introduction of noncognate base pairs during reverse transcription (RT). Combining borohydride-dependent misincorporation with ac4C's known base-sensitivity provides a unique chemical signature for this modified nucleobase. We show this unique reactivity can be used to quantitatively analyze cellular RNA acetylation, study adapters responsible for ac4C targeting, and probe the timing of RNA acetylation during ribosome biogenesis. Overall, our studies provide a chemical foundation for defining an expanding landscape of cytidine acetyltransferase activity and its impact on biology and disease.

Identification of ENTPD8 and cytidine in pancreatic cancer by metabolomic and transcriptomic conjoint analysis.

To identify metabolic pathways that were perturbed in pancreatic cancer (PC), we investigated gene-metabolite networks by integration of metabolomic and transcriptomic. In this research, we undertook the metabolomic study of 43 paired human PC samples, aiming to identify key metabolic alterations in PC. We also carried out in vitro experiments to validate that the key metabolite cytidine and its related gene ENTPD8 played an important role in PC cell proliferation. We screened out 13 metabolites differentially expressed in PC tissue (PCT) by liquid chromatography/mass spectrometry analysis on 34 metabolites, and the partial least square discrimination analysis results revealed that 9 metabolites among them were remarkably altered in PCT compared to adjacent noncancerous tissue (variable importance in projection >1, P < .05). Among the 9 metabolites, 7 might be potential biomarkers. The most significantly enriched metabolic pathway was pyrimidine metabolism. We analyzed 351 differentially expressed genes from The Cancer Genome Atlas and intersected them with Kyoto Encyclopedia of Genes and Genomes metabolic pathways. We found that ENTPD8 had a gene-metabolite association with cytidine in the CTP dephosphorylation pathway. We verified by in vitro experiments that the CTP dephosphorylation pathway was changed in PCT compared with adjacent noncancerous tissue. ENTPD8 was downregulated in PCT, causing a reduction in cytidine formation and hence weakened CTP dephosphorylation in pyrimidine metabolism.

Osmium Tag for Post-transcriptionally Modified RNA.

5-Methylcytidine (m5 C) and 5-methyluridine (m5 U) are highly abundant post-transcriptionally modified nucleotides that are observed in various natural RNAs. Such nucleotides were labeled through a chemical approach, as both underwent oxidation at the C5=C6 double bond, leading to the formation of osmium-bipyridine complexes, which could be identified by mass spectrometry. This osmium tag made it possible to distinguished m5 C and m5 U from their isomers, 2'-O-methylcytidine and 2'-O-methyluridine, respectively. Queuosine and 2-methylthio-N6 -isopentenyladenosine in tRNA were also tagged through complex formation.

Surface enhanced Raman detection of the colon cancer biomarker cytidine by using magnetized nanoparticles of the type Fe3O4/Au/Ag.

Cytidine is regarded as an early marker of colon cancer. The authors describe a surface enhanced Raman scattering (SERS) technique to detect trace levels of cytidine in urine. The Raman band at 784 cm-1 can be acquired best. Compared to earlier methods, an improvement in detection sensitivity by a factor of 6.2 × 105 is achieved by using a magnetically induced method in which cytidine is captured in the vicinity of the SERS hot spots of the type Fe3O4/Au/Ag. Cytidine can be quantified at 1 nM levels by this method which is simple and reliable. Graphical Abstract Clusters consisting of magnetite (Fe3O4) nanoparticles, gold nanoparticles and silver nanoparticles were prepared and used in a SERS based method for detection of cytidine in urine by using magnetic improvement. The lowest detectable concentration of cytidine are at the nM level.

Agrocin 108 is a 5'-cytidine nucleotide bacteriocin containing a carbocyclic phosphoryl-ascorbate group.

Agrocin 108 is a 3'-O-ß-D-xylopyranosyl-cytidine-5'-O-phosphodiester of an ascorbate-carbocyclic cyclopentenone analogue, with bacteriocin-like properties. This bacteriocin exhibits orders of magnitude greater than the inhibition zone diameter towards the indicator strain than either ampicillin or streptomycin. It has been isolated from cultures of Rhizobium rhizogenes strain K108. The structure of the agrocin 108 without detail, has been previously published. We now report a detailed structure elucidation, including the hitherto undetermined residual 5'-phospho-diester fragment by a combination of 1D and 2D NMR studies at various pH values in H2O/D2O, high resolution MS, pKa determination, and chemical degradation.

Novel ribonuclease activity of cusativin from Cucumis sativus for mapping nucleoside modifications in RNA.

A recombinant ribonuclease, cusativin, was characterized for its cytidine-specific cleavage ability of RNA to map chemical modifications. Following purification of native cusativin protein as described before (Rojo et al. Planta 194:328, 17), partial amino acid sequencing was carried out to identify the corresponding protein coding gene in cucumber genome. Cloning and heterologous expression of the identified gene in Escherichia coli resulted in successful production of active protein as a C-terminal His-tag fusion protein. The ribonuclease activity and cleavage specificity of the fusion protein were confirmed with a variety of tRNA isoacceptors and total tRNA. Characterization of cusativin digestion products by ion-pairing reverse-phase liquid chromatography coupled with mass spectrometry (IP-RP-LC-MS) analysis revealed cleavage of CpA, CpG, and CpU phosphodiester bonds at the 3'-terminus of cytidine under optimal digestion conditions. Ribose methylation or acetylation of cytosine inhibited RNA cleavage. The CpC phosphodiester bond was also resistant to cusativin-mediated RNA cleavage; a feature to our knowledge has not been reported for other nucleobase-specific ribonucleases. Here, we demonstrate the analytical utility of such a novel feature for obtaining high-sequence coverage and accurate mapping of modified residues in substrate RNAs. Graphical abstract Cytidine-specific novel ribonuclease activity of cusativin.

Label-free, direct localization and relative quantitation of the RNA nucleobase methylations m6A, m5C, m3U, and m5U by top-down mass spectrometry.

Nucleobase methylations are ubiquitous posttranscriptional modifications of ribonucleic acids (RNA) that can substantially increase the structural diversity of RNA in a highly dynamic fashion with implications for gene expression and human disease. However, high throughput, deep sequencing does not generally provide information on posttranscriptional modifications (PTMs). A promising alternative approach for the characterization of PTMs, i.e. their identification, localization, and relative quantitation, is top-down mass spectrometry (MS). In this study, we have investigated how specific nucleobase methylations affect RNA ionization in electrospray ionization (ESI), and backbone cleavage in collisionally activated dissociation (CAD) and electron detachment dissociation (EDD). For this purpose, we have developed two new approaches for the characterization of RNA methylations in mixtures of either isomers of RNA or nonisomeric RNA forms. Fragment ions from dissociation experiments were analyzed to identify the modification type, to localize the modification sites, and to reveal the site-specific, relative extent of modification for each site.

Intelligent peak deconvolution through in-depth study of the data matrix from liquid chromatography coupled with a photo-diode array detector applied to pharmaceutical analysis.

Multivariate curve resolution-alternating least squares (MCR-ALS) method was investigated for its potential to accelerate pharmaceutical research and development. The fast and efficient separation of complex mixtures consisting of multiple components, including impurities as well as major drug substances, remains a challenging application for liquid chromatography in the field of pharmaceutical analysis. In this paper we suggest an integrated analysis algorithm functioning on a matrix of data generated from HPLC coupled with photo-diode array detector (HPLC-PDA) and consisting of the mathematical program for the developed multivariate curve resolution method using an expectation maximization (EM) algorithm with a bidirectional exponentially modified Gaussian (BEMG) model function as a constraint for chromatograms and numerous PDA spectra aligned with time axis. The algorithm provided less than ±1.0% error between true and separated peak area values at resolution (Rs) of 0.6 using simulation data for a three-component mixture with an elution order of a/b/c with similarity (a/b)=0.8410, (b/c)=0.9123 and (a/c)=0.9809 of spectra at peak apex. This software concept provides fast and robust separation analysis even when method development efforts fail to achieve complete separation of the target peaks. Additionally, this approach is potentially applicable to peak deconvolution, allowing quantitative analysis of co-eluted compounds having exactly the same molecular weight. This is complementary to the use of LC-MS to perform quantitative analysis on co-eluted compounds using selected ions to differentiate the proportion of response attributable to each compound.