Computational Docking as a Tool in Guiding the Drug Design of Rutaecarpine Derivatives as Potential SARS-CoV-2 Inhibitors

COVID-19 continues to spread around the world. This is mainly because new variants of the SARS-CoV-2 virus emerge due to genomic mutations, evade the immune system and result in the effectiveness of current therapeutics being reduced. We previously established a series of detection platforms, comprising computational docking analysis, S-protein-based ELISA, pseudovirus entry, and 3CL protease activity assays, which allow us to screen a large library of phytochemicals from natural products and to determine their potential in blocking the entry of SARS-CoV-2. In this new screen, rutaecarpine (an alkaloid from Evodia rutaecarpa) was identified as exhibiting anti-SARS-CoV-2 activity. Therefore, we conducted multiple rounds of structure-activity-relationship (SAR) studies around this phytochemical and generated several rutaecarpine analogs that were subjected to in vitro evaluations. Among these derivatives, RU-75 and RU-184 displayed remarkable inhibitory activity when tested in the 3CL protease assay, S-protein-based ELISA, and pseudovirus entry assay (for both wild-type and omicron variants), and they attenuated the inflammatory response induced by SARS-CoV-2. Interestingly, RU-75 and RU-184 both appeared to be more potent than rutaecarpine itself, and this suggests that they might be considered as lead candidates for future pharmacological elaboration.


Introduction
The respiratory infection coronavirus disease 19 (COVID- 19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has had a catastrophic effect on all aspects of our lives on a global scale since the outbreak started in 2019.Various vaccines and oral treatments have been developed to treat COVID-19, and these have been employed worldwide to reduce the rate of infection and hence the course of the disease [1-3].Nevertheless, several variants of concern (VOCs) have emerged, given Greek letter monikers from alpha to omicron, and appear more transmittable and contagious than the initial wild-type virus.Indeed, several lines of evidence indicate that VOCs such as the omicron variant can trigger immune escape, which has resulted in the initial COVID-19 treatments that were developed having a reduced therapeutic effect [4,5].This suggests that there is still a need for the identification and/or development of novel antagonists that target the variants of SARS-CoV-2 to combat this wide-ranging disease.
The viral entry and replication mechanisms of SARS-CoV-2 have been well documented in previous reports [6][7][8].In brief, the spike (S)-protein from the virus binds to the angiotensin-converting enzyme 2 (ACE2) receptor on the host cells, and this leads to the virus fusing with the host cell membrane and subsequent invasion.The virus then replicates in the host cells via several key proteins, including 3CL protease, PL protease, and RNA polymerase [6,7].It is worth noting that, when compared with the original wild-type strain, VOCs contain mutations within the binding proteins, which results in them exhibiting distinct binding affinities, transmissibility, and eventual immune escape from antibodies.For example, approximately 37 residues are mutated within the S-protein of the omicron variant, and 15 of these (including K417N, T478K, and N501Y) occur in the receptor-binding domain (RBD) [8][9][10].In particular, it was reported that the K417N and E484K mutations were responsible for the immune escape of the omicron variant [9,10].Together, these proteins form a network to regulate the entry and growth of the virus in the host cells, thereafter resulting in somatic damage to the patient.Thus, the design of an inhibitor that effectively targets one or more protein(s) within this network would be extremely useful as it might lead to the development of potent therapeutics against SARS-CoV-2 and its variants [11].
In drug development, the "hit identification" and "hit-to-lead optimization" phases are crucial in discovering molecules with "drug-likeness", and they play a fundamental role in the design of new pharmaceuticals.In recent decades, fragment-based approaches have been widely utilized during these two phases [12].The concept of fragment-based drug discovery (FBDD) was first reported in 1981 by Jencks et al. [13], who proposed that compounds (while normally bound to their target proteins in a relatively weak manner) might also act as small fragments.In this way, they form efficacious interactions and thus allow for further optimization through structure-activity-relationship (SAR) studies, with the aim to grow the chemical scaffolds into highly active lead molecules [14,15].Interestingly, the structured-based drug design relies on an informative 3D structure between a ligand and a specific target.Molecular docking has become an increasingly popular and efficient tool in predicting and understanding the binding interactions between chemical molecules and biological targets, and this is beneficial to fragment-based drug design and development [16,17].
In our previous reports, we demonstrated that rutaecarpine (a pentacyclic indolopyridoquinazolinone alkaloid from Evodia rutaecarpa (Juss.)Benth) inhibited the formation of the S-protein-ACE2 complex and prevented the entry of SARS-CoV-2 into host cells [18].In addition, rutaecarpine was shown to inhibit the 3CL protease, which suggests that this alkaloid might help to suppress viral replication.Together, our data indicated that rutaecarpine might serve as a potent hit molecule in conducting SAR studies, to identify derivatives with higher efficacy.To guide our SAR investigations, we utilized computational docking stimulation to predict the binding affinities of newly designed rutaecarpine derivatives with the S-protein and 3CL protease.The derivatives were synthesized and evaluated with various in vitro assays, in which they appeared to be relatively potent.Here, we demonstrate that when compared with the parental rutaecarpine molecule, the RU-75 and RU-184 analogs have better efficacy in inhibiting the 3CL protease and preventing viral entry into host cells, as well as attenuating the inflammatory response induced by SARS-CoV-2.

In Silico-Guided Drug Design of Rutaecarpine Derivatives
The SAR study commenced with several rounds of in silico screening, where three of the rings of rutaecarpine (labelled A, B, C in Figure 1) were found to act as the binding core.Thus, they were predicted to be responsible for the binding of rutaecarpine to various protein targets, including the 3CL protease, as well as both the wild-type and omicron S-proteins.This allowed us to introduce various substituents in the D and E rings of rutaecarpine to develop active molecules with higher potency, including RU-75 and RU-184 (Figure 1).
As residues 438-506 of the S-protein have been revealed to be the RBD [8][9][10], this enabled us to establish a computational model for our in silico screening.Indeed, when the binding affinities between the newly designed analogs and the RBD were detected, RU-75 and RU-184 showed good interactions at the binding pocket by forming H-bonds with the Leu492 and Tyr449 residues, respectively (Figure 2).Interestingly, the required binding energies of these two molecules were both lower than that of rutaecarpine, with values of approximately −15.3 kJ/mol and −22.3 kJ/mol, respectively.This indicated that RU-75 and RU-184 might display more favorable binding activities at the binding site and might therefore have higher inhibitory activities than rutaecarpine.rutaecarpine to develop active molecules with higher potency, including RU-75 and RU-184 (Figure 1).As residues 438-506 of the S-protein have been revealed to be the RBD [8][9][10], this enabled us to establish a computational model for our in silico screening.Indeed, when the binding affinities between the newly designed analogs and the RBD were detected, RU-75 and RU-184 showed good interactions at the binding pocket by forming H-bonds with the Leu492 and Tyr449 residues, respectively (Figure 2).Interestingly, the required binding energies of these two molecules were both lower than that of rutaecarpine, with values of approximately −15.3 kJ/mol and −22.3 kJ/mol, respectively.This indicated that RU-75 and RU-184 might display more favorable binding activities at the binding site and might therefore have higher inhibitory activities than rutaecarpine.As residues 438-506 of the S-protein have been revealed to be the RBD [8][9][10], enabled us to establish a computational model for our in silico screening.Indeed, w the binding affinities between the newly designed analogs and the RBD were detec RU-75 and RU-184 showed good interactions at the binding pocket by forming H-bo with the Leu492 and Tyr449 residues, respectively (Figure 2).Interestingly, the requi binding energies of these two molecules were both lower than that of rutaecarpine, w values of approximately −15.3 kJ/mol and −22.3 kJ/mol, respectively.This indicated t RU-75 and RU-184 might display more favorable binding activities at the binding site a might therefore have higher inhibitory activities than rutaecarpine.To evaluate the inhibitory activities of the new analogs against VOCs, the S-protein of the omicron variant was obtained and employed in our docking model.It was observed that RU-75 and RU-184 both bound to the desired pocket of the RBD by forming pi-pi interactions with the Phe490 residue and H-bonds with the Tyr493 residue (Figure 3).As such, both molecules required lower binding energies at the pocket than rutaecarpine.This suggests that both derivatives are likely to inhibit the interaction between the S-protein of the omicron virus and ACE2.
To evaluate the inhibitory activities of the new analogs against VOCs, the S-protein of the omicron variant was obtained and employed in our docking model.It was observed that RU-75 and RU-184 both bound to the desired pocket of the RBD by forming pi-pi interactions with the Phe490 residue and H-bonds with the Tyr493 residue (Figure 3).As such, both molecules required lower binding energies at the pocket than rutaecarpine This suggests that both derivatives are likely to inhibit the interaction between the S-protein of the omicron virus and ACE2.The active site of the 3CL protease has previously been identified, and compounds binding to this site were reported to lead to protein suppression and a reduction in viral replication [19].In addition, we previously demonstrated, using an in vitro assay, that rutaecarpine could bind to the active site of the 3CL protease, expressing inhibitory effectiveness towards the protein [18].Here, somewhat intriguingly, it was found that RU-75 formed an extra H-bond with the His164 residue, while RU-184 established H-bonds with the Asn142, His164 and Glu166 residues at the desired pocket.As such, RU-75 and RU-184 were predicted to interact with the active site more favorably and required lower The active site of the 3CL protease has previously been identified, and compounds binding to this site were reported to lead to protein suppression and a reduction in viral replication [19].In addition, we previously demonstrated, using an in vitro assay, that rutaecarpine could bind to the active site of the 3CL protease, expressing inhibitory effectiveness towards the protein [18].Here, somewhat intriguingly, it was found that RU-75 formed an extra H-bond with the His164 residue, while RU-184 established H-bonds with the Asn142, His164 and Glu166 residues at the desired pocket.As such, RU-75 and RU-184 were predicted to interact with the active site more favorably and required lower energies than rutaecarpine (Figure 4).For this reason, both derivatives were subjected to organic synthesis prior to further in vitro investigations.
energies than rutaecarpine (Figure 4).For this reason, both derivatives were subje organic synthesis prior to further in vitro investigations.The synthesis pathways of RU-75 and RU-184 are shown in Figure 5. Specifically, compound 2 was obtained from amine 1 via ring closure in the presence of ethyl formate and POCl 3 .During the second ring closure, compounds 3 and 5 were both generated in an one-pot reaction.This is because both nitrogen atoms of intermediate 2 appeared to be equally active in the amide coupling reaction.Thus, good yields of both RU-75 and RU-184 were collected after stepwise nitro reduction and Buchwald-Hartwig amination [20].

In Vitro Evaluations of the Rutaecarpine Analogs
An S-protein-based ELISA was utilized to determine the binding activities of RU-75 and RU-184 to the S-protein-ACE2 complex.A standard inhibitor, provided by the supplier (calibrated to NIBSC code 20/136), was employed as a positive control (Figure 6A).Our data showed that RU-75 and RU-184 both displayed better inhibitory activities (with IC50 values of 1.03 µM and 0.78 µM, respectively) when disrupting the S-protein-ACE2 complex than RU (with an IC50 value of 1.26 µM) (Figure 6B).The maximal inhibition, triggered by RU-75 and RU-184, was at least 20% higher than that of RU.
Wild-type and omicron pseudovirus cell assays were established to evaluate whether RU-75 and RU-184 might prevent the viral entry of VOCs into host cells.First, an MTT assay was conducted to determine the potential toxicity of the analogs.Our data showed that, similar to rutaecarpine (RU), neither RU-75 nor RU-184 induced significant apoptosis, indicating that they had little impact on the cell viability (Figure S1).In addition, we found that RU-75 and RU-184 both exerted potent efficacy in blocking wild-type viral entry (Figure 7A).The IC50 values of RU, RU-75 and RU-184 were 1.34 µM, 0.46 µM and 0.21 µM, respectively.Thus, RU-75 and RU-184 were both more potent inhibitors than their parental compound.Indeed, the maximal inhibition of the analogs was almost double that of RU (Figure 7A).A similar scenario was observed when investigating the effect of RU, RU-75 and RU-184 on the entry of the omicron pseudovirus (Figure 7B).Here, the IC50 values of RU, RU-75 and RU-184 were 1.01 µM, 0.32 µM and 0.26 µM, respectively, indicating that once again the analogs were more potent inhibitors than the parent compound.

In Vitro Evaluations of the Rutaecarpine Analogs
An S-protein-based ELISA was utilized to determine the binding activities of RU-75 and RU-184 to the S-protein-ACE2 complex.A standard inhibitor, provided by the supplier (calibrated to NIBSC code 20/136), was employed as a positive control (Figure 6A).Our data showed that RU-75 and RU-184 both displayed better inhibitory activities (with IC 50 values of 1.03 µM and 0.78 µM, respectively) when disrupting the S-protein-ACE2 complex than RU (with an IC 50 value of 1.26 µM) (Figure 6B).The maximal inhibition, triggered by RU-75 and RU-184, was at least 20% higher than that of RU.
Wild-type and omicron pseudovirus cell assays were established to evaluate whether RU-75 and RU-184 might prevent the viral entry of VOCs into host cells.First, an MTT assay was conducted to determine the potential toxicity of the analogs.Our data showed that, similar to rutaecarpine (RU), neither RU-75 nor RU-184 induced significant apoptosis, indicating that they had little impact on the cell viability (Figure S1).In addition, we found that RU-75 and RU-184 both exerted potent efficacy in blocking wild-type viral entry (Figure 7A).The IC 50 values of RU, RU-75 and RU-184 were 1.34 µM, 0.46 µM and 0.21 µM, respectively.Thus, RU-75 and RU-184 were both more potent inhibitors than their parental compound.Indeed, the maximal inhibition of the analogs was almost double that of RU (Figure 7A).A similar scenario was observed when investigating the effect of RU, RU-75 and RU-184 on the entry of the omicron pseudovirus (Figure 7B).Here, the IC 50 values of RU, RU-75 and RU-184 were 1.01 µM, 0.32 µM and 0.26 µM, respectively, indicating that once again the analogs were more potent inhibitors than the parent compound.In our previous investigation, we revealed that RU attenuated the inflammatory response induced by SARS-CoV-2 and disrupted the viral entry [16].Here, we conducted a similar study with RU-75 and RU-184 and found that both of them significantly reduced the production of inflammatory cytokines (TNF-α, IL-6 and IL-1β) in the presence of pseudotyped SARS-CoV-2, when compared with the control group that was treated with the pseudotyped virus alone (Figure 8).In parallel, no significant inflammatory suppression was observed from these molecules when co-treated with lipopolysaccharides (LPS) instead of the pseudotyped SARS-CoV-2 (Figure S2).This suggests that RU-75 and RU-184 might disrupt the interaction between the S-protein and ACE2, leading to the inhibition of viral entry and a subsequent reduction in the inflammatory response.In our previous investigation, we revealed that RU attenuated the inflammatory response induced by SARS-CoV-2 and disrupted the viral entry [16].Here, we conducted a similar study with RU-75 and RU-184 and found that both of them significantly reduced the production of inflammatory cytokines (TNF-α, IL-6 and IL-1β) in the presence of pseudotyped SARS-CoV-2, when compared with the control group that was treated with the pseudotyped virus alone (Figure 8).In parallel, no significant inflammatory suppression was observed from these molecules when co-treated with lipopolysaccharides (LPS) instead of the pseudotyped SARS-CoV-2 (Figure S2).This suggests that RU-75 and RU-184 might disrupt the interaction between the S-protein and ACE2, leading to the inhibition of viral entry and a subsequent reduction in the inflammatory response.In parallel, a 3CL protease assay was carried out to determine the anti-replication effectiveness of both analogs.Indeed, both RU-75 and RU-184 displayed remarkable inhibitory activity towards the protease in a dose-dependent manner, with IC50 values of 1.31 µM and 1.03 µM, respectively (Figure 9).Thus, once again, RU-75 and RU-184 both appeared to be more potent than RU (1.60 µM).These findings support our observations on the computational docking modulation.Therefore, we suggest that the two derivatives might bind to the active site of the 3CL protease and, in this way, inhibit the action of this protein, resulting in the potential attenuation of viral replication.In parallel, a 3CL protease assay was carried out to determine the anti-replication effectiveness of both analogs.Indeed, both RU-75 and RU-184 displayed remarkable inhibitory activity towards the protease in a dose-dependent manner, with IC 50 values of 1.31 µM and 1.03 µM, respectively (Figure 9).Thus, once again, RU-75 and RU-184 both appeared to be more potent than RU (1.60 µM).These findings support our observations on the computational docking modulation.Therefore, we suggest that the two derivatives might bind to the active site of the 3CL protease and, in this way, inhibit the action of this protein, resulting in the potential attenuation of viral replication.

Discussion
As COVID-19 rapidly spread through communities, a vaccination program was rolled out on a global scale and effectively reduced the rate of infection [21].In parallel,

Discussion
As COVID-19 rapidly spread through communities, a vaccination program was rolled out on a global scale and effectively reduced the rate of infection [21].In parallel, several oral drugs were developed for clinical application that were shown to possess anti-viral efficacy.For example, the oral medicine molnupiravir was one of the major drugs that was used to treat SARS-CoV-2 patients.It was shown to be effective in attenuating the risk of hospitalization and death and provide relief from symptoms among patients with mild or moderate COVID-19 [22,23].Nevertheless, the emerging new SARS-CoV-2 VOCs displayed remarkable immune escape abilities and reduced the effectiveness of oral therapies [24,25].For example, paxlovid was initially identified as being an efficacious anti-COVID-19 treatment, but it was unable to induce fundamental viral clearance and to attenuate the mortality among patients who contracted the omicron variant [26,27].This indicates that more potent drugs are still required, specifically to tackle these VOCs.
In drug development, there are several crucial phases.These include target identification, target validation, hit identification, and hit-to-lead optimization, as well as pre-clinical and clinical evaluations [28].Drugs commonly fail in the clinic due to a lack of efficacy and safety issues; therefore, target identification and validation are becoming increasingly vital for drug discovery.Targets of interest include proteins, genes, and RNA, which enable the screening of hit molecules that express good to moderate activity for subsequent optimization.Fragment screening and high-throughput screening are two widely used strategies at this stage to yield hit candidates for subsequent expansion [28,29].The hit-to-lead phase of drug development typically commences with the determination of structure-activityrelationships to generate a large library of analogs and drive the hit molecules towards compounds with high drug-likeness.Prior to pre-clinical and clinical assessments, lead candidates are subjected to multiple rounds of in vitro assessments and in vivo studies, to identify those that possess superior efficacy and pharmaceutic properties [12,30].
Structural modifications based on rutaecarpine are well documented [31][32][33].For example, Kim et al. [31] conducted a series of SAR studies and found two rutaecarpine derivatives with improved inhibitory activity against topoisomerase I and II.In addition, Li et al. [32] reported that N-aryl/alkyl-substituted rutaecarpine analogs displayed better performance in upregulating ABCA1 activity when compared with the parent compound.In another comprehensive review, Son et al. [33] revealed that rutaecarpine analogs had more potent vasodilator activity than rutaecarpine.This inspired us to carry out intensive SAR investigations on rutaecarpine with the aim to enhance its anti-SARS-CoV-2 activity.
Previously, we established a series of detection assays.One of these was computational docking stimulation, which was found to be a relatively efficient in silico platform allowing us to screen a large library at a low cost.With these detection methods in hand, a number of rutaecarpine derivatives were designed and evaluated, and those displaying relatively good binding affinities to the targeted proteins were selected, synthesized, and evaluated.In this way, RU-75 and RU-184 (derived from rutaecarpine, a hit molecule that expressed good anti-SARS-CoV-2 activity) were generated via structural modification and multiple rounds of SAR studies.Both analogs were found to bind to the targeted proteins (i.e., the wild-type S-protein, omicron S-protein and 3CL protease) in our docking analysis.They also exhibited good potency in the S-protein-based ELISA and 3CL protease assay, as well as the wild-type and omicron pseudovirus entry tests.The two derivatives generally showed lower IC 50 values and appeared more effective than the parental phytochemical molecule-rutaecarpine.In general, the biological activity obtained from the in vitro investigations was in good agreement with the computational analysis, indicating that the molecular docking methodology was validated, and it efficiently guided the subsequent drug design and drug screening.Taken together, our investigations resulted in two potent lead candidates that might be subjected to follow-up pharmaceutical elaboration in the future, including in-depth pharmacokinetic and pharmacodynamic studies.

Cell Culture
HEK293T and RAW264.7 cells (American Type Culture Collection, Manassas, VA, USA) were cultured and maintained as described previously [34,35], except that the absorbance was measured at 490 nm.

Inhibiting SARS-CoV-2 Pseudovirus Entry
ACE2-overexpressing HEK293T cells were seeded into 48-well plates, 400 µL culture medium containing SARS-CoV-2 pseudovirus (100 µL) ± treatment reagents was applied, and the cells were incubated at 37 • C for 24 h.This medium was then replaced with fresh culture medium, and the cultures were allowed to recover for a further 48 h.Immediately prior to the start of the luciferase assay, the cultures were washed with PBS.An anti-SARS-CoV-2 neutralizing antibody (A19215, ABClonal, Woburn, MA, USA) was used as a positive control (1 µg/mL), whereas a solvent blank without the pseudovirus was used as a negative control.The inhibition percentage was determined according to the luciferase activity normalized to the activity of the enzyme without any treatment.

Luciferase Assay
The luciferase assay was conducted as previously described [35].The percentage inhibition of each sample was calculated as follows: Inhibition rate = (luciferase activity of the solvent blank − luciferase activity of the sample)/(luciferase activity of the solvent blank − luciferase activity of group without pseudovirus) × 100%.

Inflammation Test
RAW264.7 cells were incubated with a diluted concentration of 75% (v/v) wild-type SARS-CoV-2 virus or LPS (0.1 µg/mL), and then they were co-treated with various drugs for 24 h.Dexamethasone (DEX) at a concentration of 10 µM was used as a positive control.RNA was then extracted, and the levels of TNF-α, IL-6, and IL-1β mRNA were measured by RT-PCR.The RNA extraction and RT-PCR protocol were previously described [18].

S-Protein Inhibition
Spike (S)-protein inhibition was analyzed with the SARS-CoV-2 spike-ACE2 binding assay kit (ImmunoDiagnostics Ltd., Hong Kong, China), according to the manufacturer's instructions.The reaction was terminated by the addition of 2 M H 2 SO 4 , and the data were quantified using a microplate reader (FlexStation; Molecular Devices, San Jose, CA).The percentage of inhibition was calculated as follows: Percentage of inhibition = (PAvg − SAvg)/PAvg × 100%, where PAvg and SAvg are the mean OD values of the positive control and test sample, respectively.

3CL Protease Assay
Samples were tested for their ability to bind to 3CL protease on a fluorogenic substrate using the SensoLyte SARS-CoV-2 3CL protease assay kit (AnaSpec, San Jose, CA, USA), according to the manufacturer's instructions.When the 3CL protease became bound to the substrate, the fluorescence was captured using 360 nm excitation and 460 nm detection.The percentage of inhibition was calculated as follows: Percentage

Computational Docking Analysis
The chemical structure of each phytochemical was downloaded from Pubchem, and the protein structures were downloaded from the Protein Data Bank (https://www.rcsb.org/ (accessed on 28 March 2024).Virtual screening was performed with SEESAR (Version 12.0, https://www.biosolveit.de/(accessed on 28 March 2024).The docking analysis was performed as reported previously [16,36].

Conclusions
Although the global vaccine program has been widely deployed and several oral treatments have been developed, COVID-19 continues to be transmitted rapidly throughout communities due to the emerging variants of SARS-CoV-2, such as the omicron variant.Rutaecarpine was previously found to exert good anti-SARS-CoV-2 effectiveness and served as a hit molecule for structural modifications.Here, we employed computational docking to guide the drug design of rutaecarpine derivatives; these derivatives were subjected to drug screening for the following in vitro investigations.After multiple rounds of SAR studies, we identified that RU-75 and RU-184 exhibited better anti-viral efficacy than the parental molecule by disrupting the S-protein-ACE2 complex (including the wild-type and omicron variant of SARS-CoV-2), blocking viral entry and inhibiting the 3CL protease.Therefore, these two analogs might prove to be beneficial for subsequent pharmaceutical development as prospective new potent anti-COVID-19 treatments.

Figure 1 .
Figure 1.Design of rutaecarpine derivatives.The part of the scaffold shown in blue was identified as the binding core, responsible for the binding interactions, and rings D and E were the parts of the structure modified to yield the RU-75 and RU-184 analogs.

Figure 2 .
Figure 2. Docking studies to show the interaction of rutaecarpine or its derivatives (RU-75 and RU-184) with the wild-type S-protein.The information about the protein structure was downloaded

Figure 1 .
Figure 1.Design of rutaecarpine derivatives.The part of the scaffold shown in blue was identified as the binding core, responsible for the binding interactions, and rings D and E were the parts of the structure modified to yield the RU-75 and RU-184 analogs.

Figure 1 .
Figure 1.Design of rutaecarpine derivatives.The part of the scaffold shown in blue was identi as the binding core, responsible for the binding interactions, and rings D and E were the parts of structure modified to yield the RU-75 and RU-184 analogs.

Figure 2 .
Figure 2. Docking studies to show the interaction of rutaecarpine or its derivatives (RU-75 and 184) with the wild-type S-protein.The information about the protein structure was downloa

Figure 2 .
Figure 2. Docking studies to show the interaction of rutaecarpine or its derivatives (RU-75 and RU-184) with the wild-type S-protein.The information about the protein structure was downloaded from the Protein Data Bank (https://www.rcsb.org/,PDB ID: 6LZG, accessed on 28 March 2024), and the chemical structures were obtained from the Pubchem database (https://pubchem.ncbi.nlm.nih.gov/,accessed on 28 March 2024) and Chemdraw (Version 20.0).

Figure 3 .
Figure 3. Docking studies to show the interaction of rutaecarpine or its derivatives (RU-75 and RU-184) with the omicron S-protein.The protein structure was downloaded from the Protein Data Bank (https://www.rcsb.org/,PDB ID: 7T9L, accessed on 28 March 2024), and the chemical structures were obtained from the Pubchem database (https://pubchem.ncbi.nlm.nih.gov/,accessed on 28 March 2024) and Chemdraw (Version 20.0).

Figure 3 .
Figure 3. Docking studies to show the interaction of rutaecarpine or its derivatives (RU-75 and RU-184) with the omicron S-protein.The protein structure was downloaded from the Protein Data Bank (https://www.rcsb.org/,PDB ID: 7T9L, accessed on 28 March 2024), and the chemical structures were obtained from the Pubchem database (https://pubchem.ncbi.nlm.nih.gov/,accessed on 28 March 2024) and Chemdraw (Version 20.0).

Figure 4 .
Figure 4. Docking studies to show the interaction of rutaecarpine or its derivatives (RU-75 and RU-184) with the 3CL protease.The protein structure was downloaded from the Protein Data Bank (https://www.rcsb.org/,PDB ID: 6WNP, accessed on 28 March 2024), and the chemical structures were obtained from the Pubchem database (https://pubchem.ncbi.nlm.nih.gov/,accessed on 28 March 2024) and Chemdraw (Version 20.0).

Figure 6 .
Figure 6.In the S-protein-based ELISA test, RU-75 and RU-184 had a more profound inhibitory effect than RU.(A) A standard inhibitor, provided by supplier (calibrated to NIBSC code 20/136), was used as a positive control.(B) Rutaecarpine, RU-75 and RU-184 were detected at concentrations of 0.1, 1, 5, 10, 20 µM.The IC50 values are shown.The data indicate the mean + SD percentage of inhibition, compared to a no drug control (n = 4).

Figure 7 .
Figure 7. RU-75 and RU-184 block the viral entry of the wild-type (A) and omicron (B) pseudovirus.Wild-type and omicron antibodies (both used at 1 µg/mL) were used as positive controls for the wild-type and omicron viruses, respectively.RU, RU-75 and RU-184 were tested at concentrations

Figure 6 .
Figure 6.In the S-protein-based ELISA test, RU-75 and RU-184 had a more profound inhibitory effect than RU.(A) A standard inhibitor, provided by supplier (calibrated to NIBSC code 20/136), was used as a positive control.(B) Rutaecarpine, RU-75 and RU-184 were detected at concentrations of 0.1, 1, 5, 10, 20 µM.The IC 50 values are shown.The data indicate the mean ± SD percentage of inhibition, compared to a no drug control (n = 4).

Figure 6 .
Figure 6.In the S-protein-based ELISA test, RU-75 and RU-184 had a more profound inhibitory effect than RU.(A) A standard inhibitor, provided by supplier (calibrated to NIBSC code 20/136), was used as a positive control.(B) Rutaecarpine, RU-75 and RU-184 were detected at concentrations of 0.1, 1, 5, 10, 20 µM.The IC50 values are shown.The data indicate the mean + SD percentage of inhibition, compared to a no drug control (n = 4).

Figure 7 .
Figure 7. RU-75 and RU-184 block the viral entry of the wild-type (A) and omicron (B) pseudovirus.Wild-type and omicron antibodies (both used at 1 µg/mL) were used as positive controls for the wild-type and omicron viruses, respectively.RU, RU-75 and RU-184 were tested at concentrations

Figure 7 .
Figure 7. RU-75 and RU-184 block the viral entry of the wild-type (A) and omicron (B) pseudovirus.Wild-type and omicron antibodies (both used at 1 µg/mL) were used as positive controls for the wild-type and omicron viruses, respectively.RU, RU-75 and RU-184 were tested at concentrations of 0.1, 1, 5, 10, 20 µM.The IC 50 values are shown.The data indicate the mean ± SD percentage of inhibition, compared to a no drug control (n = 5).

Molecules 2024 ,
29, x FOR PEER REVIEW 8 of 15 of 0.1, 1, 5, 10, 20 µM.The IC50 values are shown.The data indicate the mean ± SD percentage of inhibition, compared to a no drug control (n = 5).

Figure 8 .
Figure 8. RU-75 and RU-184 attenuate the inflammatory response induced by the wild-type pseudovirus.The mRNA levels of TNF-α, IL-6 and IL-1β were determined by RT-PCR.The virus was used at a concentration of 75% (v/v), and dexamethasone (DEX) was utilized as the positive control (10 µM).The data show the mean ± SD (n = 4) mRNA fold change compared to a blank group (x basal), and the asterisks indicate statistically significant differences such that * p < 0.05, ** p < 0.01, *** p < 0.001 when compared with the blank.

Figure 8 .
Figure 8. RU-75 and RU-184 attenuate the inflammatory response induced by the wild-type pseudovirus.The mRNA levels of TNF-α, IL-6 and IL-1β were determined by RT-PCR.The virus was used at a concentration of 75% (v/v), and dexamethasone (DEX) was utilized as the positive control (10 µM).The data show the mean ± SD (n = 4) mRNA fold change compared to a blank group (x basal), and the asterisks indicate statistically significant differences such that * p < 0.05, ** p < 0.01, *** p < 0.001 when compared with the blank.

Molecules 2024 , 15 Figure 9 .
Figure 9. RU-75 and RU-184 suppress 3CL protease activity.GC376 (100 µM) was used as a positive control.RU-75 and RU-184 both inhibited the 3CL protease in a dose-dependent manner.The IC50 values are shown.The data indicate the mean ± SD percentage of inhibition, compared to a no drug control (n = 4).

Figure 9 .
Figure 9. RU-75 and RU-184 suppress 3CL protease activity.GC376 (100 µM) was used as a positive control.RU-75 and RU-184 both inhibited the 3CL protease in a dose-dependent manner.The IC 50 values are shown.The data indicate the mean ± SD percentage of inhibition, compared to a no drug control (n = 4).
of inhibition = (PAvg, b − SAvg, b)/PAvg, b × 100%, where PAvg, b and SAvg, b are the mean fluorescence of the positive control and test samples, respectively, subtracted from the mean fluorescence of a blank.