Synthesis, Anti-TMV Activities, and Action Mechanisms of a Novel Cytidine Peptide Compound.

Cytidine has a broad range of applications in the pharmaceutical field as an intermediate of antitumor or antiviral agent. Here, a series of new cytidine peptide compounds were synthesized using cytidine and Boc group-protected amino acids and analyzed for their antiviral activities against tobacco mosaic virus (TMV). Among these compounds, the structure of an effective antiviral cytidine peptide SN11 was characterized by 1H NMR, 13C NMR, and high-resolution mass spectrometer. The compound SN11 has a molecular formula of C15H22N6O8 and is named 2-amino-N-(2- ((1- (3,4-dihydroxy-5-(hydroxymethyl) tetrahydrofuran-2-yl) -2-oxo-1,2-dihydropyrimidin-4-yl) amino) -2-oxyethyl) amino). The protection, inactivation, and curation activities of SN11 at a concentration of 500 µg/mL against TMV in Nicotiana glutinosa were 82.6%, 84.2%, and 72.8%, respectively. SN11 also effectively suppressed the systemic transportation of a recombinant TMV carrying GFP reporter gene (p35S-30B:GFP) in Nicotiana benthamiana by reducing viral accumulation to 71.3% in the upper uninoculated leaves and inhibited the systemic infection of TMV in Nicotiana tabacum plants. Furthermore, the results of RNA-seq showed that compound SN11 induced differential expression of genes involved in the biogenesis and function of ribosome, plant hormone signal transduction, plant pathogen interaction, and chromatin. These results validate the antiviral mechanisms of the cytidine peptide compound and provide a theoretical basis for their potential application in the management of plant virus diseases.

Chemistry of installing epitranscriptomic 5-modified cytidines in RNA oligomers.

Studies of 5-hydroxymethylcytidine (hm5C), 5-formylcytidine (f5C) and 5-carboxycytidine (ca5C) modifications as products of the 5-methylcytidine (m5C) oxidative demethylation pathway in cellular mRNAs constitute an important element of the new epitranscriptomic field of research. The dynamic process of m5C conversion and final turnover to the parent cytidine is considered a post-transcriptional layer of gene-expression regulation. However, the regulatory mechanism associated with epitranscriptomic cytidine modifications remains largely unknown. Therefore, oligonucleotides containing m5C oxidation products are of great value for the next generation of biochemical, biophysical, and structural studies on their function, metabolism, and contribution to human diseases. Herein, we summarize the synthetic strategies developed for the incorporation of hm5C, f5C and ca5C into RNA oligomers by phosphoramidite chemistry, including post-synthetic C5-cytidine functionalization and enzymatic methods.

Investigating the Effect of Chemical Modifications on the Ribose Sugar Conformation, Watson-Crick Base Pairing, and Intrastrand Stacking Interactions: A Theoretical Approach.

Over the last few decades, chemically modified sugars have been incorporated into nucleic acid-based therapeutics to improve their pharmacological potential. Chemical modification can influence the sugar conformation, Watson-Crick hydrogen (W-C) bonding, and nucleobase stacking interactions, which play major roles in the structural integrity and dynamic properties of nucleic acid duplexes. In this study, we categorized 33 uridine (U*) and cytidine (C*) sugar modifications and calculated their sugar conformational parameters. We also calculated the Watson-Crick hydrogen bond energies of the modified RNA-type base pairs (U*:A and C*:G) using DFT and sSAPT0 methods. The W-C base pairing energy calculations suggested that the South-type modified sugar strengthens the C*:G base pair and weakens the U*:A base pair compared to the unmodified one. In contrast, the North-type sugar modifications form weaker C*:G base pair and marginally stronger U*:A base pair compared to the South-type modified sugars. Moreover, intrastrand base stacking energies were calculated for 15 modifications incorporated at the fourth position in 7-mer non-self-complementary RNA duplexes [(GCAU*GAC)2 and (GCAC*GAC)2], utilizing molecular dynamics simulation and quantum mechanical (DFT and sSAPT0) methods. The sugar modifications were found to have minimal effect on the intrastrand base-stacking interactions. However, the glycol nucleic acid modification disturbs the intrastrand base-stacking significantly, corroborating the experimental data.

Rational design of base, sugar and backbone modifications improves ADAR-mediated RNA editing.

AIMers are short, chemically modified oligonucleotides that induce A-to-I RNA editing through interaction with endogenous adenosine deaminases acting on RNA (ADAR) enzymes. Here, we describe the development of new AIMer designs with base, sugar and backbone modifications that improve RNA editing efficiency over our previous design. AIMers incorporating a novel pattern of backbone and 2' sugar modifications support enhanced editing efficiency across multiple sequences. Further efficiency gains were achieved through incorporation of an N-3-uridine (N3U), in place of cytidine (C), in the 'orphan base' position opposite the edit site. Molecular modeling suggests that N3U might enhance ADAR catalytic activity by stabilizing the AIMer-ADAR interaction and potentially reducing the energy required to flip the target base into the active site. Supporting this hypothesis, AIMers containing N3U consistently enhanced RNA editing over those containing C across multiple target sequences and multiple nearest neighbor sequence combinations. AIMers combining N3U and the novel pattern of 2' sugar chemistry and backbone modifications improved RNA editing both in vitro and in vivo. We provide detailed N3U synthesis methods and, for the first time, explore the impact of N3U and its analogs on ADAR-mediated RNA editing efficiency and targetable sequence space.

Discovery of A G-Quadruplex Unwinder That Unleashes the Translation of G-Quadruplex-Containing mRNA without Inducing DNA Damage.

To explore the mechanisms and therapeutic strategies for G-quadruplex (G4) mediated diseases, it is crucial to manipulate and intervene in intracellular G4 structures using small molecular tools. While hundreds of G4 stabilizers have been developed, there is a significant gap in the availability of G4 unwinding agents. Here, we propose a strategy to disrupt G-quadruplexes by forming G-C hydrogen bonds with chemically modified cytidine trimers. We validated a good G4 unwinder, the 2'-F cytidine trimer (2'-F C3). 2'-F C3 does not inhibit cell growth nor cause severe DNA damage at a concentration below 10 µM. Moreover, 2'-F C3 does not affect gene transcription nor RNA splicing, while it significantly enhances the translation of G4-containing mRNA and upregulates RNA splicing, RNA processing and cell cycle pathways. The discovery of this G4 unwinder provides a functional tool for the chemical modulation of G4s in living cells.

Crosslinking of base-modified RNAs by synthetic DYW-KP base editors implicates an enzymatic lysine as the nitrogen donor for U-to-C RNA editing.

Sequence-specific cytidine to uridine (C-to-U) and adenosine to inosine editing tools can alter RNA and DNA sequences and utilize a hydrolytic deamination mechanism requiring an active site zinc ion and a glutamate residue. In plant organelles, DYW-PG domain containing enzymes catalyze C-to-U edits through the canonical deamination mechanism. Proteins developed from consensus sequences of the related DYW-KP domain family catalyze what initially appeared to be uridine to cytidine (U-to-C) edits leading to this investigation into the U-to-C editing mechanism. The synthetic DYW-KP enzyme KP6 was found sufficient for C-to-U editing activity stimulated by the addition of carboxylic acids in vitro. Despite addition of putative amine/amide donors, U-to-C editing by KP6 could not be observed in vitro. C-to-U editing was found not to be concomitant with U-to-C editing, discounting a pyrimidine transaminase mechanism. RNAs containing base modifications were highly enriched in interphase fractions consistent with covalent crosslinks to KP6, KP2, and KP3 proteins. Mass spectrometry of purified KP2 and KP6 proteins revealed secondary peaks with mass shifts of 319 Da. A U-to-C crosslinking mechanism was projected to explain the link between crosslinking, RNA base changes, and the ∼319 Da mass. In this model, an enzymatic lysine attacks C4 of uridine to form a Schiff base RNA-protein conjugate. Sequenced RT-PCR products from the fern Ceratopteris richardii indicate U-to-C base edits do not preserve proteinaceous crosslinks in planta. Hydrolysis of a protonated Schiff base conjugate releasing cytidine is hypothesized to explain the completed pathway in plants.

Enhanced ac4C detection in RNA via chemical reduction and cDNA synthesis with modified dNTPs.

The functional analysis of epitranscriptomic modifications in RNA is constrained by a lack of methods that accurately capture their locations and levels. We previously demonstrated that the RNA modification N4-acetylcytidine (ac4C) can be mapped at base resolution through sodium borohydride reduction to tetrahydroacetylcytidine (tetrahydro-ac4C), followed by cDNA synthesis to misincorporate adenosine opposite reduced ac4C sites, culminating in C:T mismatches at acetylated cytidines (RedaC:T). However, this process is relatively inefficient, resulting in <20% C:T mismatches at a fully modified ac4C site in 18S rRNA. Considering that ac4C locations in other substrates including mRNA are unlikely to reach full penetrance, this method is not ideal for comprehensive mapping. Here, we introduce "RetraC:T" (reduction to tetrahydro-ac4C and reverse transcription with amino-dATP to induce C:T mismatches) as a method with enhanced ability to detect ac4C in cellular RNA. In brief, RNA is reduced through NaBH4 or the closely related reagent sodium cyanoborohydride (NaCNBH3) followed by cDNA synthesis in the presence of a modified DNA nucleotide, 2-amino-dATP, that preferentially binds to tetrahydro-ac4C. Incorporation of the modified dNTP substantially improved C:T mismatch rates, reaching stoichiometric detection of ac4C in 18S rRNA. Importantly, 2-amino-dATP did not result in truncated cDNA products nor increase mismatches at other locations. Thus, modified dNTPs are introduced as a new addition to the toolbox for detecting ac4C at base resolution.

A path from synthesis to emergency use authorization of molnupiravir as a COVID-19 therapy.

Coronaviruses are a group of enveloped viruses with non-segmented, single-stranded, and positive-sense RNA genomes. It belongs to the 'Coronaviridae family', responsible for various diseases, including the common cold, SARS, and MERS. The COVID-19 pandemic, which began in March 2020, has affected 209 countries, infected over a million people, and claimed over 50,000 lives. Significant efforts have been made by repurposing several approved drugs including antiviral, to combat the COVID-19 pandemic. Molnupiravir is found to be the first orally acting efficacious drug to treat COVID-19 cases. It was approved for medical use in the UK in November 2021 and other countries, including USFDA, which granted approval an emergency use authorization (EUA) for treating adults with mild to moderate COVID-19 patients. Considering the importance of molnupiravir, the present review deals with its various synthetic strategies, pharmacokinetics, bio-efficacy, toxicity, and safety profiles. The comprehensive information along with critical analysis will be very handy for a wide range of audience including medicinal chemists in the arena of antiviral drug discovery especially anti-viral drugs against any variant of COVID-19.

Structural insights into the decoding capability of isoleucine tRNAs with lysidine and agmatidine.

The anticodon modifications of transfer RNAs (tRNAs) finetune the codon recognition on the ribosome for accurate translation. Bacteria and archaea utilize the modified cytidines, lysidine (L) and agmatidine (agm2C), respectively, in the anticodon of tRNAIle to decipher AUA codon. L and agm2C contain long side chains with polar termini, but their functions remain elusive. Here we report the cryogenic electron microscopy structures of tRNAsIle recognizing the AUA codon on the ribosome. Both modifications interact with the third adenine of the codon via a unique C-A geometry. The side chains extend toward 3' direction of the mRNA, and the polar termini form hydrogen bonds with 2'-OH of the residue 3'-adjacent to the AUA codon. Biochemical analyses demonstrated that AUA decoding is facilitated by the additional interaction between the polar termini of the modified cytidines and 2'-OH of the fourth mRNA residue. We also visualized cyclic N6-threonylcarbamoyladenosine (ct6A), another tRNA modification, and revealed a molecular basis how ct6A contributes to efficient decoding.

Substrate specificity of Mycobacterium tuberculosis tRNA terminal nucleotidyltransferase toxin MenT3.

Mycobacterium tuberculosis transfer RNA (tRNA) terminal nucleotidyltransferase toxin, MenT3, incorporates nucleotides at the 3'-CCA end of tRNAs, blocking their aminoacylation and inhibiting protein synthesis. Here, we show that MenT3 most effectively adds CMPs to the 3'-CCA end of tRNA. The crystal structure of MenT3 in complex with CTP reveals a CTP-specific nucleotide-binding pocket. The 4-NH2 and the N3 and O2 atoms of cytosine in CTP form hydrogen bonds with the main-chain carbonyl oxygen of P120 and the side chain of R238, respectively. MenT3 expression in Escherichia coli selectively reduces the levels of seryl-tRNASers, indicating specific inactivation of tRNASers by MenT3. Consistently, MenT3 incorporates CMPs into tRNASer most efficiently, among the tested E. coli tRNA species. The longer variable loop unique to class II tRNASers is crucial for efficient CMP incorporation into tRNASer by MenT3. Replacing the variable loop of E. coli tRNAAla with the longer variable loop of M. tuberculosis tRNASer enables MenT3 to incorporate CMPs into the chimeric tRNAAla. The N-terminal positively charged region of MenT3 is required for CMP incorporation into tRNASer. A docking model of tRNA onto MenT3 suggests that an interaction between the N-terminal region and the longer variable loop of tRNASer facilitates tRNA substrate selection.

Photochemical Stability of 5-Methylcytidine Relative to Cytidine: Photophysical Insight for mRNA Therapeutic Applications.

5-Methylcytidine (5mCyd) has recently been investigated with renewed interest in its utilization in mRNA therapeutics. However, its photostability following exposure to electromagnetic radiation has been overlooked. This Letter compares the photostability and excited-state dynamics of 5mCyd with those of the canonical RNA nucleoside, cytidine (Cyd), using steady-state and femtosecond transient absorption spectroscopy under physiologic conditions. 5mCyd is shown to have a 5-fold higher fluorescence yield and a 5-fold longer 1ππ* excited-state decay lifetime. Importantly, however, the excited-state population in 5mCyd decays primarily by internal conversion, with a photodegradation rate 3 times smaller than that in Cyd. In Cyd, the population of a 1nπ* state with a lifetime of ca. 45 ps is implicated in the formation 6-hydroxycytidine and other photoproducts.

Study of protonated dimers of cytosine, cytidine, and deoxycytidine using survival yield method and quantum mechanics calculations.

RATIONALE: Cytosine and its conjugates are prone to form protonated, triply-bonded dimers. Therefore, the nucleic-acid cytosine-rich sequence forms the four-stranded noncanonical secondary structure known as the intercalated motif (i-motif). This process has resulted in studies on cytosine protonated dimers. This communication focuses on the protonated dimers of cytosine and its nucleoside using the survival yield (SY) method and quantum mechanics calculations. METHODS: To obtain the precursor ion fragmentation curve, the plot of SY against Ecomδ , the product ion spectra of the protonated dimers were obtained using a Waters/Micromass Q-TOF Premier mass spectrometer. Quantum mechanics calculations were performed using GAUSSIAN 16, and full geometry optimizations and energy calculations were performed within the density functional theory framework at B3LYP/6-31G(d,p). RESULTS: The precursor ion fragmentation curve allowed the rating of the gas-phase stabilities of the analyzed protonated dimers. Substitution of sugar moiety at N1 cytosine atom decreased the gas-phase stabilities of the protonated dimers. The deoxycytidine dimer was found to be more stable than the cytidine dimer and cytidine-deoxycytidine dimer. Quantum chemical calculations indicated that cytosine aminohydroxy tautomer may be involved in the formation of protonated cytosine-cytosine nucleoside dimers but not in the formation of cytosine dimers. CONCLUSIONS: The results obtained for nucleoside dimers indicated that the SY method may reflect the i-motif stabilities observed under physiological conditions. Therefore, the analysis of other protonated dimers of variously substituted cytosine-cytosine nucleoside using the SY method may be important to study the effect of cytosine substitution on the i-motif stabilities. Cytosine tautomer containing C2-OH… N(2H)-C4 moiety may be involved in the formation of protonated cytosine-cytosine nucleoside dimers.

Sequence-Specific Fluorescence Turn-On Sensing of RNA by DNA Probes Incorporating the Tricyclic Cytidine Analogue DEAtC.

Sequence-specific fluorescent probes for RNA are widely used in microscopy applications such as fluorescence in situ hybridization and a growing number of newer approaches to live-cell RNA imaging. The sequence specificity of most of these approaches relies on differential hybridization of the probe to the correct target. Competing sequences with only one or two base mismatches are prone to causing off-target recognition. Here, we report the sequence-specific fluorescent detection of model RNA targets using a tricyclic cytidine analogue DEAtC that is included as a surrogate for natural cytidine in DNA probe strands and that reports directly on Watson-Crick base pairing. The DEAtC-containing DNA oligonucleotide probes exhibit an average 8-fold increase in fluorescence intensity when hybridized to matched RNA with DEAtC base paired with G and little fluorescence turn-on when DEAtC is base paired with A. Duplex structure determination by NMR, time-resolved fluorescence studies, and Stern-Volmer quenching experiments suggest that the combination of greater π stacking and narrower grooves in the A-form DNA-RNA heteroduplex provides additional shielding and favorable electronic interactions between bases, explaining why DEAtC's fluorescence turn-on response to RNA targets is typically 3-fold greater than for DNA targets.

Design, Synthesis, and Evaluation of a Cross-Linked Oligonucleotide as the First Nanomolar Inhibitor of APOBEC3A.

Drug resistance is a major problem associated with anticancer chemo- and immunotherapies. Recent advances in the understanding of resistance mechanisms have revealed that enzymes of the APOBEC3 (A3) family contribute to the development of drug resistance in multiple cancers. A3 enzymes are polynucleotide cytidine deaminases that convert cytosine to uracil (C→U) in single-stranded DNA (ssDNA) and in this way protect humans against viruses and mobile retroelements. On the other hand, cancer cells use A3s, especially A3A and A3B, to mutate human DNA, and thus by increasing rates of evolution, cancer cells escape adaptive immune responses and resist drugs. However, as A3A and A3B are non-essential for primary metabolism, their inhibition opens up a strategy to augment existing anticancer therapies and suppress cancer evolution. To test our hypothesis that pre-shaped ssDNA mimicking the U-shape observed in ssDNA-A3 complexes can provide a better binder to A3 enzymes, a Cu(I)-catalyzed azide-alkyne cycloaddition was used to cross-link two distant modified nucleobases in ssDNA. The resultant cytosine-containing substrate, where the cytosine sits at the apex of the loop, was deaminated faster by the engineered C-terminal domain of A3B than a standard, linear substrate. The cross-linked ssDNA was converted into an A3 inhibitor by replacing the 2'-deoxycytidine in the preferred TCA substrate motif by 2'-deoxyzebularine, a known inhibitor of single nucleoside cytidine deaminases. This strategy yielded the first nanomolar inhibitor of engineered A3BCTD and wild-type A3A (Ki = 690 ± 140 and 360 ± 120 nM, respectively), providing a platform for further development of powerful A3 inhibitors.

Isoxazole Nucleosides as Building Blocks for a Plausible Proto-RNA.

The question of how RNA, as the principal carrier of genetic information evolved is fundamentally important for our understanding of the origin of life. The RNA molecule is far too complex to have formed in one evolutionary step, suggesting that ancestral proto-RNAs (first ancestor of RNA) may have existed, which evolved over time into the RNA of today. Here we show that isoxazole nucleosides, which are quickly formed from hydroxylamine, cyanoacetylene, urea and ribose, are plausible precursors for RNA. The isoxazole nucleoside can rearrange within an RNA-strand to give cytidine, which leads to an increase of pairing stability. If the proto-RNA contains a canonical seed-nucleoside with defined stereochemistry, the seed-nucleoside can control the configuration of the anomeric center that forms during the in-RNA transformation. The results demonstrate that RNA could have emerged from evolutionarily primitive precursor isoxazole ribosides after strand formation.

Protonation-Dependent Sequencing of 5-Formylcytidine in RNA.

Chemical modification of cytidine in noncoding RNAs plays a key role in regulating translation and disease. However, the distribution and dynamics of many of these modifications remain unknown due to a lack of sensitive site-specific sequencing technologies. Here, we report a protonation-dependent sequencing reaction for the detection of 5-formylcytidine (5fC) and 5-carboxycytidine (5caC) in RNA. First, we evaluate how protonation combined with electron-withdrawing substituents alters the molecular orbital energies and reduction of modified cytidine nucleosides, highlighting 5fC and 5caC as reactive species. Next, we apply this reaction to detect these modifications in synthetic oligonucleotides as well as endogenous human transfer RNA (tRNA). Finally, we demonstrate the utility of our method to characterize a patient-derived model of 5fC deficiency, where it enables facile monitoring of both pathogenic loss and exogenous rescue of NSUN3-dependent 5fC within the wobble base of human mitochondrial tRNAMet. These studies showcase the ability of protonation to enhance the reactivity and sensitive detection of 5fC in RNA and more broadly provide a molecular foundation for using optimized sequencing reactions to better understand the role of oxidized RNA cytidine residues in diseases.

A Systematic Review of the Global Intervention for SARS-CoV-2 Combating: From Drugs Repurposing to Molnupiravir Approval.

The rising outbreak of SARS-CoV-2 continues to unfold all over the world. The development of novel effective antiviral drugs to fight against SARS-CoV-2 is a time cost. As a result, some specific FDA-approved drugs have already been repurposed and authorized for COVID-19 treatment. The repurposed drugs used were either antiviral or non-antiviral drugs. Accordingly, the present review thoroughly focuses on the repurposing efficacy of these drugs including clinical trials experienced, the combination therapies used, the novel methods followed for treatment, and their future perspective. Therefore, drug repurposing was regarded as an effective avenue for COVID-19 treatment. Recently, molnupiravir is a prodrug antiviral medication that was approved in the United Kingdom in November 2021 for the treatment of COVID-19. On the other hand, PF-07321332 is an oral antiviral drug developed by Pfizer. For the treatment of COVID-19, the PF-07321332/ritonavir combination medication is used in Phase III studies and was marketed as Paxlovid. Herein, we represented the almost history of combating COVID-19 from repurposing to the recently available oral anti-SARS-CoV-2 candidates, as a new hope to end the current pandemic.

Computational Analysis of Molnupiravir.

In this work, we report in-depth computational studies of three plausible tautomeric forms, generated through the migration of two acidic protons of the N4-hydroxylcytosine fragment, of molnupiravir, which is emerging as an efficient drug to treat COVID-19. The DFT calculations were performed to verify the structure of these tautomers, as well as their electronic and optical properties. Molecular docking was applied to examine the influence of the structures of the keto-oxime, keto-hydroxylamine and hydroxyl-oxime tautomers on a series of the SARS-CoV-2 proteins. These tautomers exhibited the best affinity behavior (-9.90, -7.90, and -9.30 kcal/mol, respectively) towards RdRp-RTR and Nonstructural protein 3 (nsp3_range 207-379-MES).

Zebularine protects against blood-brain-barrier (BBB) disruption through increasing the expression of zona occludens-1 (ZO-1) and vascular endothelial (VE)-cadherin.

Blood-brain-barrier (BBB) disruption is an important pathological characteristic of ischemic stroke (IS) and mainly results from dysfunction of brain vascular endothelial cells and tight junctions. Zebularine is a novel inhibitor of DNA methyltransferase (DNMT). Here, we assessed its effects on BBB disruption in IS. Firstly, we reported that Zebularine maintained BBB integrity in middle cerebral artery occlusion (MCAO) mice by increasing the expressions of zona occludens-1 (ZO-1) and vascular endothelial (VE)-cadherin. Importantly, we found that Zebularine reduced the production of pro-inflammatory cytokines, attenuated brain edema, and improved neurological deficits. In in vitro experiments, the bEnd.3 brain endothelial cells were exposed to oxygen and glucose deprivation/reoxygenation (OGD/R), and the protective effects of Zebularine were assessed. Our findings demonstrated that Zebularine prevented OGD/R-induced cytotoxicity by reducing the release of lactate dehydrogenase (LDH). Additionally, Zebularine protected bEnd.3 cells against OGD/R-induced hyper-permeability and reduction of trans-endothelial electrical resistance (TEER). Notably, we found that treatment with Zebularine activated the Adenosine 5'-monophosphate (AMP)-activated protein kinase (AMPK) pathway by increasing the phosphorylation of adenosine monophosphate-activated protein kinase α (AMPKα). Blockage of AMPKα using its specific inhibitor compound C abolished the beneficial effects of Zebularine in mitigating endothelial hyper-permeability by reducing the expressions of ZO-1 and VE-cadherin. These findings suggest that the protective effects of Zebularine against OGD/R-induced endothelial hyper-permeability are mediated by the activation of AMPKα. In conclusion, our study sheds light on the potential application of Zebularine in the treatment of IS.

Potential COVID-19 Drug Candidates Based on Diazinyl-Thiazol-Imine Moieties: Synthesis and Greener Pastures Biological Study.

A novel series of 1-aryl-N-[4-phenyl-5-(arylazo)thiazol-2-yl)methanimines has been synthesized via the condensation of 2-amino-4-phenyl-5-arylazothiazole with various aromatic aldehydes. The synthesized imines were characterized by spectroscopic techniques, namely 1H and 13C-NMR, FTIR, MS, and Elemental Analysis. A molecular comparative docking study for 3a-f was calculated, with reference to two approved drugs, Molnupiravir and Remdesivir, using 7BQY (Mpro; PDB code 7BQY; resolution: 1.7 A°) under identical conditions. The binding scores against 7BQY were in the range of -7.7 to -8.7 kcal/mol for 3a-f. The high scores of the compounds indicated an enhanced binding affinity of the molecules to the receptor. This is due to the hydrophobic interactions and multi-hydrogen bonds between 3a-f ligands and the receptor's active amino acid residues. The main aim of using in silco molecular docking was to rank 3a-f with respect to the approved drugs, Molnupiravir and Remdesivir, using free energy methods as greener pastures. A further interesting comparison presented the laydown of the ligands before and after molecular docking. These results and other supporting statistical analyses suggested that ligands 3a-f deserve further investigation in the context of potential therapeutic agents for COVID-19. Free-cost, PASS, SwissADME, and Way2drug were used in this research paper to determine the possible biological activities and cytotoxicity of 3a-f.