Interaction of 8-anilinonaphthalene-1-sulfonate with SARS-CoV-2 main protease and its application as a fluorescent probe for inhibitor identification

Enhanced targeting efficacy of baicalein analogues on the dimeric state of SARS-CoV-2 3CL protease compared to monomeric state

The COVID-19 pandemic caused by the novel coronavirus, SARS-CoV-2, has been a global threat affecting the entire world. It is a single-stranded RNA virus that belongs to the coronavirus family. In SARS-CoV2, the 3CL protease protein significantly contributes to viral replication and is responsible for viral polyprotein cleavage. These factors make 3CL protease a promising drug target to inhibit the growth of SARS-CoV-2. In this study, using in silico approaches, we have targeted the 3CL protease of SARS-CoV-2 to identify promising antiviral candidates for COVID-19 treatment. Here, 463 structural analogs of Baicalein compounds were collected initially, and by employing the quantitative structure-activity relationship (QSAR) technique on 76 antiviral compounds, screening was done against monomeric and dimeric versions of the target protein. Further, based on the molecular interaction studies and MD simulation, followed by validation of the obtained simulation trajectories using PCA and MM/PBSA calculation, it was observed that ligands showed better binding stability with dimeric proteins than monomeric proteins and can be used as suitable therapeutic candidates for SARS-CoV2 treatment. The MD simulation showed a favorable, robust outcome for the 46885476 when bound to the dimeric state. It matched the control in the number of hydrogen bonds and conformational stability. This molecule also directly impacted the catalytic dyads of the protein, suggesting potential inhibitory action. In addition, this study helps to accelerate the drug development process against SARS-CoV2 through the observed in-silico results, which need to be validated using clinical experiments in future studies.

Alpha and gamma mangostins inhibit wild-type B SARS-CoV-2 more effectively than the SARS-CoV-2 variants and the major target is unlikely the 3C-like protease

Background

Anti-SARS-CoV-2 and immunomodulatory drugs are important for treating clinically severe patients with respiratory distress symptoms. Alpha- and gamma-mangostins (AM and GM) were previously reported as potential 3C-like protease (3CLpro) and Angiotensin-converting enzyme receptor 2 (ACE2)-binding inhibitors in silico.

Objective

We aimed to evaluate two active compounds, AM and GM, from Garcinia mangostana for their antivirals against SARS-CoV-2 in live virus culture systems and their cytotoxicities using standard methods. Also, we aimed to prove whether 3CLpro and ACE2 neutralization were major targets and explored whether any additional targets existed.

Methods

We tested the translation and replication efficiencies of SARS-CoV-2 in the presence of AM and GM. Initial and subgenomic translations were evaluated by immunofluorescence of SARS-CoV-2 3CLpro and N expressions at 16 h after infection. The viral genome was quantified and compared with the untreated group. We also evaluated the efficacies and cytotoxicities of AM and GM against four strains of SARS-CoV-2 (wild-type B, B.1.167.2, B.1.36.16, and B.1.1.529) in Vero E6 cells. The potential targets were evaluated using cell-based anti-attachment, time-of-drug addition, in vitro 3CLpro activities, and ACE2-binding using a surrogated viral neutralization test (sVNT). Moreover, additional targets were explored using combinatorial network-based interactions and Chemical Similarity Ensemble Approach (SEA).

Results

AM and GM reduced SARS-CoV-2 3CLpro and N expressions, suggesting that initial and subgenomic translations were globally inhibited. AM and GM inhibited all strains of SARS-CoV-2 at EC50 of 0.70-3.05 μM, in which wild-type B was the most susceptible strain (EC50 0.70-0.79 μM). AM was slightly more efficient in the variants (EC50 0.88-2.41 μM), resulting in higher selectivity indices (SI 3.65-10.05), compared to the GM (EC50 0.94-3.05 μM, SI 1.66-5.40). GM appeared to be more toxic than AM in both Vero E6 and Calu-3 cells. Cell-based anti-attachment and time-of-addition suggested that the potential molecular target could be at the post-infection. 3CLpro activity and ACE2 binding were interfered with in a dose-dependent manner but were insufficient to be a major target. Combinatorial network-based interaction and chemical similarity ensemble approach (SEA) suggested that fatty acid synthase (FASN), which was critical for SARS-CoV-2 replication, could be a target of AM and GM.

Conclusion

AM and GM inhibited SARS-CoV-2 with the highest potency at the wild-type B and the lowest at the B.1.1.529. Multiple targets were expected to integratively inhibit viral replication in cell-based system.

FMO-guided design of darunavir analogs as HIV-1 protease inhibitors

The prevalence of HIV-1 infection continues to pose a significant global public health issue, highlighting the need for antiretroviral drugs that target viral proteins to reduce viral replication. One such target is HIV-1 protease (PR), responsible for cleaving viral polyproteins, leading to the maturation of viral proteins. While darunavir (DRV) is a potent HIV-1 PR inhibitor, drug resistance can arise due to mutations in HIV-1 PR. To address this issue, we developed a novel approach using the fragment molecular orbital (FMO) method and structure-based drug design to create DRV analogs. Using combinatorial programming, we generated novel analogs freely accessible via an on-the-cloud mode implemented in Google Colab, Combined Analog generator Tool (CAT). The designed analogs underwent cascade screening through molecular docking with HIV-1 PR wild-type and major mutations at the active site. Molecular dynamics (MD) simulations confirmed the assess ligand binding and susceptibility of screened designed analogs. Our findings indicate that the three designed analogs guided by FMO, 19-0-14-3, 19-8-10-0, and 19-8-14-3, are superior to DRV and have the potential to serve as efficient PR inhibitors. These findings demonstrate the effectiveness of our approach and its potential to be used in further studies for developing new antiretroviral drugs.

Targeting the receptor binding domain and heparan sulfate binding for antiviral drug development against SARS-CoV-2 variants

The emergence of SARS-CoV-2 variants diminished the efficacy of current antiviral drugs and vaccines. Hence, identifying highly conserved sequences and potentially druggable pockets for drug development was a promising strategy against SARS-CoV-2 variants. In viral infection, the receptor-binding domain (RBD) proteins are essential in binding to the host receptor. Others, Heparan sulfate (HS), widely distributed on the surface of host cells, is thought to play a central role in the viral infection cycle of SARS-CoV-2. Therefore, it might be a reasonable strategy for antiviral drug design to interfere with the RBD in the HS binding site. In this study, we used computational approaches to analyze multiple sequences of coronaviruses and reveal important information about the binding of HS to RBD in the SARS-CoV-2 spike protein. Our results showed that the potential hot-spots, including R454 and E471, in RBD, exhibited strong interactions in the HS-RBD binding region. Therefore, we screened different compounds in the natural product database towards these hot-spots to find potential antiviral candidates using LibDock, Autodock vina and furthermore applying the MD simulation in AMBER20. The results showed three potential natural compounds, including Acetoside (ACE), Hyperoside (HYP), and Isoquercitrin (ISO), had a strong affinity to the RBD. Our results demonstrate a feasible approach to identify potential antiviral agents by evaluating the binding interaction between viral glycoproteins and host receptors. The present study provided the applications of the structure-based computational approach for designing and developing of new antiviral drugs against SARS-CoV-2 variants.

A combination of structure-based virtual screening and experimental strategies to identify the potency of caffeic acid ester derivatives as SARS-CoV-2 3CLpro inhibitor from an in-house database

Drug development requires significant time and resources, and computer-aided drug discovery techniques that integrate chemical and biological spaces offer valuable tools for the process. This study focused on the field of COVID-19 therapeutics and aimed to identify new active non-covalent inhibitors for 3CLpro, a key protein target. By combining in silico and in vitro approaches, an in-house database was utilized to identify potential inhibitors. The drug-likeness criteria were considered to pre-filter 553 compounds from 12 groups of natural products. Using structure-based virtual screening, 296 compounds were identified that matched the chemical features of SARS-CoV-2 3CLpro peptidomimetic inhibitor pharmacophore models. Subsequent molecular docking resulted in 43 hits with high binding affinities. Among the hits, caffeic acid analogs showed significant interactions with the 3CLpro active site, indicating their potential as promising candidates. To further evaluate their efficacy, enzyme-based assays were conducted, revealing that two ester derivatives of caffeic acid (4k and 4l) exhibited more than a 30% reduction in 3CLpro activity. Overall, these findings suggest that the screening approach employed in this study holds promise for the discovery of novel anti-SARS-CoV-2 therapeutics. Furthermore, the methodology could be extended for optimization or retrospective evaluation to enhance molecular targeting and antiviral efficacy of potential drug candidates.

Docking and Electronic Structure of Rutin, Myricetin, and Baicalein Targeting 3CLpro

Understanding the role of 3CLpro protease for SARS-CoV-2 replication and knowing the potential of flavonoid molecules like rutin, myricetin, and baicalein against 3CLpro justify an investigation into their inhibition. This study investigates possible bonds and reactivity descriptors of rutin, myricetin, and baicalein through conformational and electronic properties. Density functional theory was used to determine possible interactions. Analyses were carried out through the molecular electrostatic potential, electron localization function, Fukui function descriptors based on frontier orbitals, and non-covalent interactions. A docking study was performed using a resolution of 1.55 Å for 3CLpro to analyze the interactions of rutin, myricetin, and baicalein. Scores of structures showed that rutin is the best ligand, followed by myricetin and baicalein. Docking studies showed that baicalein and rutin can establish effective interactions with residues of the catalytic dyad (Cys145 and His41), but just rutin forms a hydrogen bond. Myricetin, in turn, could not establish an effective interaction with Cys145. Baicalein interaction arose with active residues such as Arg188, Val186, Gln189, and Gln192. Interactions of rutin and myricetin with Arg188 and Gln189 were also found. A critical interaction was observed only for rutin with the hydroxyls of ring A with His41, and also for Cys145 with rings B and C, which is probably related to the highest score of rutin.

In Silico and In Vitro Potential of FDA-Approved Drugs for Antimalarial Drug Repurposing against Plasmodium Serine Hydroxymethyltransferases

Malaria has spread in many countries, with a 12% increase in deaths after the coronavirus disease 2019 pandemic. Malaria is one of the most concerning diseases in the Greater Mekong subregion, showing increased drug-resistant rates. Serine hydroxymethyltransferase (SHMT), a key enzyme in the deoxythymidylate synthesis pathway, has been identified as a promising antimalarial drug target due to its conserved folate binding pocket. This study used a molecular docking approach to screen 2509 US Food and Drug Administration (FDA)-approved drugs against seven Plasmodium SHMT structures. Eight compounds had significantly lower binding energies than the known SHMT inhibitors pyrazolopyran(+)-86, tetrahydrofolate, and antimalarial drugs, ranging from 4 to 10 kcal/mol. Inhibition assays testing the eight compounds against Plasmodium falciparum SHMT (PfSHMT) showed that amphotericin B was a competitive inhibitor of PfSHMT with a half-maximal inhibitory concentration (IC50) of 106 ± 1 μM. Therefore, a 500 ns molecular dynamics simulation of PfSHMT/PLS/amphotericin B was performed. The backbone root-mean-square deviation of the protein-ligand complex indicated the high complex stability during simulations, supported by its radius of gyration, hydrogen-bond interactions, and number of atom contacts. The appreciable binding affinity of amphotericin B for PfSHMT was indicated by their solvated interaction energy (-11.15 ± 0.09 kcal/mol) and supported by strong ligand-protein interactions (≥80% occurrences) with its essential residues (i.e., Y78, K151, N262, F266, and V365) predicted by pharmacophore modeling and per-residue decomposition free energy methods. Therefore, our findings identify a promising new PfSHMT inhibitor, albeit with less inhibitory activity, and suggest a core structure that differs from that of previous SHMT inhibitors, thus being a rational approach for novel antimalarial drug design.

Galangal-Cinnamon Spice Mixture Blocks the Coronavirus Infection Pathway through Inhibition of SARS-CoV-2 MPro, Three HCoV-229E Targets; Quantum-Chemical Calculations Support In Vitro Evaluation

Natural products such as domestic herbal drugs which are easily accessible and cost-effective can be used as a complementary treatment in mild and moderate COVID-19 cases. This study aimed to detect and describe the efficiency of phenolics detected in the galangal-cinnamon mixture in the inhibition of SARS-CoV-2's different protein targets. The potential antiviral effect of galangal-cinnamon aqueous extract (GCAE) against Low Pathogenic HCoV-229E was assessed using cytopathic effect inhibition assay and the crystal violet method. Low Pathogenic HCoV-229E was used as it is safer for in vitro laboratory experimentation and due to the conformation and the binding pockets similarity between HCoV-229E and SARS-CoV-2 MPro. The GCAE showed a significant antiviral effect against HCoV-229E (IC50 15.083 µg/mL). Twelve phenolic compounds were detected in the extract with ellagic, cinnamic, and gallic acids being the major identified phenolic acids, while rutin was the major identified flavonoid glycoside. Quantum-chemical calculations were made to find molecular properties using the DFT/B3LYP method with 6-311++G(2d,2p) basis set. Quantum-chemical values such as EHOMO, ELUMO, energy gap, ionization potential, chemical hardness, softness, and electronegativity values were calculated and discussed. Phenolic compounds detected by HPLC-DAD-UV in the GCAE were docked into the active site of 3 HCoV-229E targets (PDB IDs. 2ZU2, 6U7G, 7VN9, and 6WTT) to find the potential inhibitors that block the Coronavirus infection pathways from quantum and docking data for these compounds. There are good adaptations between the theoretical and experimental results showing that rutin has the highest activity against Low Pathogenic HCoV-229E in the GCAE extract.

Identification of Promising Sulfonamide Chalcones as Inhibitors of SARS-CoV-2 3CLpro through Structure-Based Virtual Screening and Experimental Approaches

3CLpro is a viable target for developing antiviral therapies against the coronavirus. With the urgent need to find new possible inhibitors, a structure-based virtual screening approach was developed. This study recognized 75 pharmacologically bioactive compounds from our in-house library of 1052 natural product-based compounds that satisfied drug-likeness criteria and exhibited good bioavailability and membrane permeability. Among these compounds, three promising sulfonamide chalcones were identified by combined theoretical and experimental approaches, with SWC423 being the most suitable representative compound due to its competitive inhibition and low cytotoxicity in Vero E6 cells (EC50 = 0.89 ± 0.32 μM; CC50 = 25.54 ± 1.38 μM; SI = 28.70). The binding and stability of SWC423 in the 3CLpro active site were investigated through all-atom molecular dynamics simulation and fragment molecular orbital calculation, indicating its potential as a 3CLpro inhibitor for further SARS-CoV-2 therapeutic research. These findings suggested that inhibiting 3CLpro with a sulfonamide chalcone such as SWC423 may pave the effective way for developing COVID-19 treatments.

Designing Potent α-Glucosidase Inhibitors: A Synthesis and QSAR Modeling Approach for Biscoumarin Derivatives

Nineteen biscoumarins were synthesized, well-characterized, and evaluated against α-glucosidases in vitro. Of these, six compounds (10, 12, 16, and 17-19) were newly synthesized and not previously reported in the chemical literature. The majority of the synthesized derivatives demonstrated significant inhibitory activity. A quantitative structure-activity relationship (QSAR) model was developed, revealing a strong correlation between the anti-α-glucosidase activity and selected molecular descriptors. Based on this model, two new compounds (18 and 19) were designed, which exhibited the strongest inhibition with IC₅₀ values of 0.62 and 1.21 μM, respectively, when compared to the positive control (acarbose) with an IC₅₀ value of 93.63 μM. Enzyme kinetic studies of compounds 18 and 19 revealed their competitive inhibition with K_i values of 3.93 and 1.80 μM, respectively. Computational studies demonstrated that compound 18 could be inserted into the original binding site (OBS) of α-glucosidase MAL12 and form multiple hydrophobic interactions with nearby amino acids, with the bromo group playing an essential role in enhancing the binding strength and stability at the OBS of the enzyme based on the quantum mechanical calculations using the fragment molecular orbital method. These findings provide valuable insights into the design of potent α-glucosidase inhibitors, which may have potential therapeutic applications in the treatment of diabetes and related diseases.

Tenofovir antiviral drug solubility enhancement with β-cyclodextrin inclusion complex and in silico study of potential inhibitor against SARS-CoV-2 main protease (Mpro)

Tenofovir (TFR) is an antiviral drug commonly used to fight against viral diseases infection due to its good potency and high genetic barrier to drug resistance. In physiological conditions, TFR is less water soluble, more unstable, and less permeable, limiting its effective therapeutic applications. In addition to their use in treating the Coronavirus disease 2019 (COVID-19), cyclodextrins (CDs) are also being used as a molecule to develop therapies for other diseases due to its enhance solubility and stability. This study is designed to synthesize and characterization of β-CD:TFR inclusion complex and its interaction against SARS-CoV-2 (M^Pro) protein (PDB ID;7cam). Several techniques were used to characterize the prepared β-CD:TFR inclusion complex, including UV-Visible, FT-IR, XRD, SEM, TGA, and DSC, which provided appropriate evidence to confirm the formation. A 1:1 stoichiometry was determined for β-CD:TFR inclusion complex in aqueous medium from UV-Visible absorption spectra by using the Benesi-Hildebrand method. Phase solubility studies proposed that β-CD enhanced the excellent solubility of TFR and the stability constant was obtained at 863 ± 32 M^-1. Moreover, the molecular docking confirmed the experimental results demonstrated the most desirable mode of TFR encapsulated into the β-CD nanocavity via hydrophobic interactions and possible hydrogen bonds. Moreover, TFR was validated in the β-CD:TFR inclusion complex as potential inhibitors against SARS-CoV-2 main protease (M^pro) receptors by using in silico methods. The enhanced solubility, stability, and antiviral activity against SARS-CoV-2 (M^Pro) suggest that β-CD:TFR inclusion complexes can be further used as feasible water-insoluble antiviral drug carriers in viral disease infection.

Quantum Chemical Interaction Analysis between SARS-CoV-2 Main Protease and Ensitrelvir Compared with Its Initial Screening Hit

A non-covalent oral drug targeting SARS-CoV-2 main protease (Mpro), ensitrelvir (Xocova), has been developed using structure-based drug design (SBDD). To elucidate the factors responsible for enhanced inhibitory activities from an in silico screening hit compound to ensitrelvir, we analyzed the interaction energies of the inhibitors with each residue of Mpro using fragment molecular orbital (FMO) calculations. This analysis reveals that functional group conversion for P1' and P1 parts in the inhibitors increases the strength of existing interactions with Mpro and also provides novel interactions for ensitrelvir; the associated changes in the conformation of Mpro induce further interactions for ensitrelvir in other parts, including hydrogen bonds, a halogen bond, and π-orbital interactions. Thus, we illuminate the promising strategies of SBDD for leading ensitrelvir to get higher activity against Mpro by elucidating microscopic interactions through FMO-based analysis. These detailed mechanism findings, including water cross-linkings, will help to design novel inhibitors in SBDD.

N-Containing α-Mangostin Analogs via Smiles Rearrangement as the Promising Cytotoxic, Antitrypanosomal, and SARS-CoV-2 Main Protease Inhibitory Agents

New N-containing xanthone analogs of α-mangostin were synthesized via one-pot Smiles rearrangement. Using cesium carbonate in the presence of 2-chloroacetamide and catalytic potassium iodide, α-mangostin (1) was subsequently transformed in three steps to provide ether 2, amide 3, and amine 4 in good yields at an optimum ratio of 1:3:3, respectively. The evaluation of the biological activities of α-mangostin and analogs 2-4 was described. Amine 4 showed promising cytotoxicity against the non-small-cell lung cancer H460 cell line fourfold more potent than that of cisplatin. Both compounds 3 and 4 possessed antitrypanosomal properties against Trypanosoma brucei rhodesiense at a potency threefold stronger than that of α-mangostin. Furthermore, ether 2 gave potent SARS-CoV-2 main protease inhibition by suppressing 3-chymotrypsinlike protease (3CL^pro) activity approximately threefold better than that of 1. Fragment molecular orbital method (FMO-RIMP2/PCM) indicated the improved binding interaction of 2 in the 3CL^pro active site regarding an additional ether moiety. Thus, the series of N-containing α-mangostin analogs prospectively enhance druglike properties based on isosteric replacement and would be further studied as potential biotically active chemical entries, particularly for anti-lung-cancer, antitrypanosomal, and anti-SARS-CoV-2 main protease applications.

SARS-CoV-2 main protease (3CLpro) interaction with acyclovir antiviral drug/methyl-β-cyclodextrin complex: Physiochemical characterization and molecular docking

During the current outbreak of the novel coronavirus disease 2019 (COVID-19), researchers have examined several antiviral drugs with the potential to inhibit the proliferation of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The antiviral drug acyclovir (AVR), which is used to treat COVID-19, in complex with methyl-β-cyclodextrin (Mβ-CD) was examined in the solution and solid phases. UV-visible and fluorescence spectroscopic analyses confirmed that the guest (AVR) was included inside the host (Mβ-CD) cavity. A solid inclusion complex of AVR was prepared by co-precipitation, physical mixing, kneading, and bath sonication methods at a 1:1 ratio of Mβ-CD:AVR. The prepared Mβ-CD:AVR inclusion complex was characterized using Fourier transform infrared spectroscopy (FTIR), X-ray powder diffraction (XRD), thermogravimetric analysis (TGA), differential scanning calorimetry (DSC), and scanning electron microscopy (SEM) analysis. Phase solubility studies indicated the Mβ-CD:AVR inclusion complex exhibited a higher stability constant and linear enhancement in AVR solubility with increasing Mβ-CD concentrations. In silico analysis of the Mβ-CD/AVR inclusion complex confirmed that AVR drugs show potential as inhibitors of SARS-CoV-2 3C-like protease (3CL^pro) receptors. Results obtained using the PatchDock and FireDock servers indicated that the most favorable docking ligand was Mβ-CD:AVR, which interacted with SARS-CoV-2 (3CL^Pro) protease inhibitors with high geometric shape complementarity scores (2522 and 5872) and atomic contact energy (-313.77 and -214.70 kcal mol^-1). Our results suggest that the Mβ-CD/AVR inclusion complex inhibits the main protease of SARS-CoV-2, although further wet-lab experiments are needed to verify these findings.

Plant Molecular Pharming and Plant-Derived Compounds towards Generation of Vaccines and Therapeutics against Coronaviruses

The current century has witnessed infections of pandemic proportions caused by Coronaviruses (CoV) including severe acute respiratory syndrome-related CoV (SARS-CoV), Middle East respiratory syndrome-related CoV (MERS-CoV) and the recently identified SARS-CoV2. Significantly, the SARS-CoV2 outbreak, declared a pandemic in early 2020, has wreaked devastation and imposed intense pressure on medical establishments world-wide in a short time period by spreading at a rapid pace, resulting in high morbidity and mortality. Therefore, there is a compelling need to combat and contain the CoV infections. The current review addresses the unique features of the molecular virology of major Coronaviruses that may be tractable towards antiviral targeting and design of novel preventative and therapeutic intervention strategies. Plant-derived vaccines, in particular oral vaccines, afford safer, effectual and low-cost avenues to develop antivirals and fast response vaccines, requiring minimal infrastructure and trained personnel for vaccine administration in developing countries. This review article discusses recent developments in the generation of plant-based vaccines, therapeutic/drug molecules, monoclonal antibodies and phytochemicals to preclude and combat infections caused by SARS-CoV, MERS-CoV and SARS-CoV-2 viruses. Efficacious plant-derived antivirals could contribute significantly to combating emerging and re-emerging pathogenic CoV infections and help stem the tide of any future pandemics.

Promising SARS-CoV-2 main protease inhibitor ligand-binding modes evaluated using LB-PaCS-MD/FMO

Parallel cascade selection molecular dynamics-based ligand binding-path sampling (LB-PaCS-MD) was combined with fragment molecular orbital (FMO) calculations to reveal the ligand path from an aqueous solution to the SARS-CoV-2 main protease (M^pro) active site and to customise a ligand-binding pocket suitable for delivering a potent inhibitor. Rubraxanthone exhibited mixed-inhibition antiviral activity against SARS-CoV-2 M^pro, relatively low cytotoxicity, and high cellular inhibition. However, the atomic inhibition mechanism remains ambiguous. LB-PaCS-MD/FMO is a hybrid ligand-binding evaluation method elucidating how rubraxanthone interacts with SARS-CoV-2 M^pro. In the first step, LB-PaCS-MD, which is regarded as a flexible docking, efficiently samples a set of ligand-binding pathways. After that, a reasonable docking pose of LB-PaCS-MD is evaluated by the FMO calculation to elucidate a set of protein-ligand interactions, enabling one to know the binding affinity of a specified ligand with respect to a target protein. A possible conformation was proposed for rubraxanthone binding to the SARS-CoV-2 M^pro active site, and allosteric inhibition was elucidated by combining blind docking with k-means clustering. The interaction profile, key binding residues, and considerable interaction were elucidated for rubraxanthone binding to both M^pro sites. Integrated LB-PaCS-MD/FMO provided a more reasonable complex structure for ligand binding at the SARS-CoV-2 M^pro active site, which is vital for discovering and designing antiviral drugs.

Identification of repurposing therapeutics toward SARS-CoV-2 main protease by virtual screening

SARS-CoV-2 causes the current global pandemic coronavirus disease 2019. Widely-available effective drugs could be a critical factor in halting the pandemic. The main protease (3CL^pro) plays a vital role in viral replication; therefore, it is of great interest to find inhibitors for this enzyme. We applied the combination of virtual screening based on molecular docking derived from the crystal structure of the peptidomimetic inhibitors (N3, 13b, and 11a), and experimental verification revealed FDA-approved drugs that could inhibit the 3CL^pro of SARS-CoV-2. Three drugs were selected using the binding energy criteria and subsequently performed the 3CL^pro inhibition by enzyme-based assay. In addition, six common drugs were also chosen to study the 3CL^pro inhibition. Among these compounds, lapatinib showed high efficiency of 3CL^pro inhibition (IC₅₀ value of 35 ± 1 μM and K_i of 23 ± 1 μM). The binding behavior of lapatinib against 3CL^pro was elucidated by molecular dynamics simulations. This drug could well bind with 3CL^pro residues in the five subsites S1’, S1, S2, S3, and S4. Moreover, lapatinib’s key chemical pharmacophore features toward SAR-CoV-2 3CL^pro shared important HBD and HBA with potent peptidomimetic inhibitors. The rational design of lapatinib was subsequently carried out using the obtained results. Our discovery provides an effective repurposed drug and its newly designed analogs to inhibit SARS-CoV-2 3CL^pro.

Semi-Synthesis of N-Aryl Amide Analogs of Piperine from Piper nigrum and Evaluation of Their Antitrypanosomal, Antimalarial, and Anti-SARS-CoV-2 Main Protease Activities

Piper nigrum, or black pepper, produces piperine, an alkaloid that has diverse pharmacological activities. In this study, N-aryl amide piperine analogs were prepared by semi-synthesis involving the saponification of piperine (1) to yield piperic acid (2) followed by esterification to obtain compounds 3, 4, and 5. The compounds were examined for their antitrypanosomal, antimalarial, and anti-SARS-CoV-2 main protease activities. The new 2,5-dimethoxy-substituted phenyl piperamide 5 exhibited the most robust biological activities with no cytotoxicity against mammalian cell lines, Vero and Vero E6, as compared to the other compounds in this series. Its half-maximal inhibitory concentration (IC₅₀) for antitrypanosomal activity against Trypanosoma brucei rhodesiense was 15.46 ± 3.09 μM, and its antimalarial activity against the 3D7 strain of Plasmodium falciparum was 24.55 ± 1.91 μM, which were fourfold and fivefold more potent, respectively, than the activities of piperine. Interestingly, compound 5 inhibited the activity of 3C-like main protease (3CL^Pro) toward anti-SARS-CoV-2 activity at the IC₅₀ of 106.9 ± 1.2 μM, which was threefold more potent than the activity of rutin. Docking and molecular dynamic simulation indicated that the potential binding of 5 in the 3CL^pro active site had the improved binding interaction and stability. Therefore, new aryl amide analogs of piperine 5 should be investigated further as a promising anti-infective agent against human African trypanosomiasis, malaria, and COVID-19.

Halogenated Baicalein as a Promising Antiviral Agent toward SARS-CoV-2 Main Protease

The coronavirus disease pandemic is a constant reminder that global citizens are in imminent danger of exposure to emerging infectious diseases. Therefore, developing a technique for inhibitor discovery is essential for effective drug design. Herein, we proposed fragment molecular orbital (FMO)-based virtual screening to predict the molecular binding energy of potential severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main protease inhibitors. The integration of quantum mechanical approaches and trajectory analysis from a microsecond molecular dynamics simulation was used to identify potential inhibitors. We identified brominated baicalein as a potent inhibitor of the SARS-CoV-2 main protease and confirmed its inhibitory activity in an in vitro assay. Brominated baicalein did not demonstrate significant toxicity in either in vitro or in vivo studies. The pair interaction energy from FMO-RIMP2/PCM and inhibitory constants based on the protease enzyme assay suggested that the brominated baicalein could be further developed into novel SARS-CoV-2 protease inhibitors.

Pharmacogenetics and Precision Medicine Approaches for the Improvement of COVID-19 Therapies

Many drugs are being administered to tackle coronavirus disease 2019 (COVID-19) pandemic situations without establishing clinical effectiveness or tailoring safety. A repurposing strategy might be more effective and successful if pharmacogenetic interventions are being considered in future clinical studies/trials. Although it is very unlikely that there are almost no pharmacogenetic data for COVID-19 drugs, however, from inferring the pharmacokinetic (PK)/pharmacodynamic(PD) properties and some pharmacogenetic evidence in other diseases/clinical conditions, it is highly likely that pharmacogenetic associations are also feasible in at least some COVID-19 drugs. We strongly mandate to undertake a pharmacogenetic assessment for at least these drug-gene pairs (atazanavir-UGT1A1, ABCB1, SLCO1B1, APOA5; efavirenz-CYP2B6; nevirapine-HLA, CYP2B6, ABCB1; lopinavir-SLCO1B3, ABCC2; ribavirin-SLC28A2; tocilizumab-FCGR3A; ivermectin-ABCB1; oseltamivir-CES1, ABCB1; clopidogrel-CYP2C19, ABCB1, warfarin-CYP2C9, VKORC1; non-steroidal anti-inflammatory drugs (NSAIDs)-CYP2C9) in COVID-19 patients for advancing precision medicine. Molecular docking and computational studies are promising to achieve new therapeutics against SARS-CoV-2 infection. The current situation in the discovery of anti-SARS-CoV-2 agents at four important targets from in silico studies has been described and summarized in this review. Although natural occurring compounds from different herbs against SARS-CoV-2 infection are favorable, however, accurate experimental investigation of these compounds is warranted to provide insightful information. Moreover, clinical considerations of drug-drug interactions (DDIs) and drug-herb interactions (DHIs) of the existing repurposed drugs along with pharmacogenetic (e.g., efavirenz and CYP2B6) and herbogenetic (e.g., andrographolide and CYP2C9) interventions, collectively called multifactorial drug-gene interactions (DGIs), may further accelerate the development of precision COVID-19 therapies in the real-world clinical settings.

Ligand Binding Path Sampling Based on Parallel Cascade Selection Molecular Dynamics: LB-PaCS-MD

Parallel cascade selection molecular dynamics (PaCS-MD) is a rare-event sampling method that generates transition pathways between a reactant and product. To sample the transition pathways, PaCS-MD repeats short-time MD simulations from important configurations as conformational resampling cycles. In this study, PaCS-MD was extended to sample ligand binding pathways toward a target protein, which is referred to as LB-PaCS-MD. In a ligand-concentrated environment, where multiple ligand copies are randomly arranged around the target protein, LB-PaCS-MD allows for the frequent sampling of ligand binding pathways. To select the important configurations, we specified the center of mass (COM) distance between each ligand and the relevant binding site of the target protein, where snapshots generated by the short-time MD simulations were ranked by their COM distance values. From each cycle, snapshots with smaller COM distance values were selected as the important configurations to be resampled using the short-time MD simulations. By repeating conformational resampling cycles, the COM distance values gradually decreased and converged to constants, meaning that a set of ligand binding pathways toward the target protein was sampled by LB-PaCS-MD. To demonstrate relative efficiency, LB-PaCS-MD was applied to several proteins, and their ligand binding pathways were sampled more frequently than conventional MD simulations.

Discovery of C-12 dithiocarbamate andrographolide analogues as inhibitors of SARS-CoV-2 main protease: In vitro and in silico studies

•

Andrographolide analogues were found to inhibit SARS-CoV-2 main protease.

•

The compounds 3k, 3l, 3m and 3t showed promising in vitro inhibitory activity.

•

Most of the candidates could bind well to the SARS-CoV-2 main protease active site.

A global crisis of coronavirus disease 2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has impacted millions of people’s lives throughout the world. In parallel to vaccine development, identifying potential antiviral agents against SARS-CoV-2 has become an urgent need to combat COVID-19. One of the most attractive drug targets for discovering anti-SARS-CoV-2 agents is the main protease (M^pro), which plays a pivotal role in the viral life cycle. This study aimed to elucidate a series of twenty-one 12-dithiocarbamate-14-deoxyandrographolide analogues as SARS-CoV-2 M^pro inhibitors using in vitro and in silico studies. These compounds were initially screened for the inhibitory activity toward SARS-CoV-2 M^pro by in vitro enzyme-based assay. We found that compounds 3k, 3l, 3m and 3t showed promising inhibitory activity against SARS-CoV-2 M^pro with >50% inhibition at 10 μM. Afterward, the binding mode of each compound in the active site of SARS-CoV-2 M^pro was explored by molecular docking. The optimum docked complexes were then chosen and subjected to molecular dynamic (MD) simulations. The MD results suggested that all studied complexes were stable along the simulation time, and most of the compounds could fit well with the SARS-CoV-2 M^pro active site, particularly at S1, S2 and S4 subsites. The per-residue decomposition free energy calculations indicated that the hot-spot residues essential for ligand binding were T25, H41, C44, S46, M49, C145, H163, M165, E166, L167, D187, R188, Q189 and T190. Therefore, the obtained information from the combined experimental and computational techniques could lead to further optimization of more specific and potent andrographolide analogues toward SARS-CoV-2 M^pro.

Monitoring protein conformational changes using fluorescent nanoantennas

Understanding the relationship between protein structural dynamics and function is crucial for both basic research and biotechnology. However, methods for studying the fast dynamics of structural changes are limited. Here, we introduce fluorescent nanoantennas as a spectroscopic technique to sense and report protein conformational changes through noncovalent dye-protein interactions. Using experiments and molecular simulations, we detect and characterize five distinct conformational states of intestinal alkaline phosphatase, including the transient enzyme-substrate complex. We also explored the universality of the nanoantenna strategy with another model protein, Protein G and its interaction with antibodies, and demonstrated a rapid screening strategy to identify efficient nanoantennas. These versatile nanoantennas can be used with diverse dyes to monitor small and large conformational changes, suggesting that they could be used to characterize diverse protein movements or in high-throughput screening applications.

Food phytochemicals, epigallocatechin gallate and myricetin, covalently bind to the active site of the coronavirus main protease in vitro

SARS-CoV-2 main protease is a possible target for protection against viral infection. This study examined the inhibitory effect of food phytochemicals on the main protease of SARS-CoV-2 by determining a cleaved product after chromatographic separation. First, 37 phytochemicals, including glycosides and metabolites, were screened at 20 µM; epigallocatechin gallate, myricetin, theaflavin, herbacetin, piceatannol, myricitrin, and isothiocyanates inhibited the enzyme in varying degrees. The IC₅₀ values were estimated from 0.4 to 33.3 µM against the 0.5-µM enzyme. The dose-dependent adduction of epigallocatechin gallate and myricetin was confirmed by quinone staining of protein blotted onto a membrane. The enzyme activity was decreased by increasing the concentration of the two phytochemicals, accompanied by increasing the respective adducted molecule estimated by intact mass spectrometry. Reduced glutathione canceled the formation of conjugate and the inhibitory effect of epigallocatechin gallate or myricetin on the enzyme, suggesting that the formation of the quinone moiety in the phytochemicals is critical for the inhibition. The covalent binding of epigallocatechin gallate or myricetin to the cysteine residue at the active site was confirmed by analyzing peptides from the chymotrypsin-digested main protease.

Special Features of COVID-19 in the FMODB: Fragment Molecular Orbital Calculations and Interaction Energy Analysis of SARS-CoV-2-Related Proteins

SARS-CoV-2 is the causative agent of coronavirus (known as COVID-19), the virus causing the current pandemic. There are ongoing research studies to develop effective therapeutics and vaccines against COVID-19 using various methods and many results have been published. The structure-based drug design of SARS-CoV-2-related proteins is promising, however, reliable information regarding the structural and intra- and intermolecular interactions is required. We have conducted studies based on the fragment molecular orbital (FMO) method for calculating the electronic structures of protein complexes and analyzing their quantitative molecular interactions. This enables us to extensively analyze the molecular interactions in residues or functional group units acting inside the protein complexes. Such precise interaction data are available in the FMO database (FMODB) (https://drugdesign.riken.jp/FMODB/). Since April 2020, we have performed several FMO calculations on the structures of SARS-CoV-2-related proteins registered in the Protein Data Bank. We have published the results of 681 structures, including three structural proteins and 11 nonstructural proteins, on the COVID-19 special page (as of June 8, 2021). In this paper, we describe the entire COVID-19 special page of the FMODB and discuss the calculation results for various proteins. These data not only aid the interpretation of experimentally determined structures but also the understanding of protein functions, which is useful for rational drug design for COVID-19.