In Silico Study of Coumarins and Quinolines Derivatives as Potent Inhibitors of SARS-CoV-2 Main Protease

Furanocoumarins promote proteasomal degradation of viral HBx protein and down-regulate cccDNA transcription and replication of hepatitis B virus

Nucleot(s)ide analogues, the current antiviral treatments against chronic hepatitis B (CHB) infection, are non-curative due to their inability to eliminate covalently closed circular DNA (cccDNA) from the infected hepatocytes. Preclinical studies have shown that coumarin derivatives can effectively reduce the HBV DNA replication. We evaluated the antiviral efficacy of thirty new coumarin derivatives in cell culture models for studying HBV. Furanocoumarins Fc-20 and Fc-31 suppressed the levels of pre-genomic RNA as well as cccDNA, and reduced the secretion of virions, HBsAg and HBeAg. The antiviral efficacies of Fc-20 and Fc31 improved further when used in combination with the hepatitis B antiviral drug Entecavir. There was a marked reduction in the intracellular HBx level in the presence of these furanocoumarins due to proteasomal degradation resulting in the down-regulation of HBx-dependent viral genes. Importantly, both Fc-20 and Fc-31 were non-cytotoxic to cells even at high concentrations. Further, our molecular docking studies confirmed a moderate to high affinity interaction between furanocoumarins and viral HBx via residues Ala3, Arg26 and Lys140. These data suggest that furanocoumarins could be developed as a new therapeutic for CHB infection.

Synthesis, SARS-CoV-2 main protease inhibition, molecular docking and in silico ADME studies of furanochromene-quinoline hydrazone derivatives

Seven furanochromene-quinoline derivatives containing a hydrazone linker were synthesized by condensing a furanochromene hydrazide with quinoline 2-, 3-, 4-, 5-, 6-, and 8-carbaldehydes, including 8-hydroxyquinoline-2-carbaldehye. Structure-activity correlations were investigated to determine the influence of the location of the hydrazone linker on the quinoline unit on SARS-CoV-2 Mpro enzyme inhibition. The 3-, 5-, 6- and 8-substituted derivatives showed moderate inhibition of SARS-CoV-2 Mpro with IC50 values ranging from 16 to 44 μM. Additionally, all of the derivatives showed strong interaction with the SARS-CoV-2 Mpro substrate binding pocket, with docking energy scores ranging from -8.0 to -8.5 kcal/mol. These values are comparable to that of N3 peptide (-8.1 kcal/mol) and more favorable than GC-373 (-7.6 kcal/mol) and ML-188 (-7.5 kcal/mol), all of which are known SARS-CoV-2 Mpro inhibitors. Furthermore, in silico absorption, distribution, metabolism, and excretion (ADME) profiles indicate that the derivatives have good drug-likeness properties. Overall, this study highlights the potential of the furanochromene-quinoline hydrazone scaffold as a SARS-CoV-2 Mpro inhibitor.

In silico investigation of falcipain-2 inhibition by hybrid benzimidazole-thiosemicarbazone antiplasmodial agents: A molecular docking, molecular dynamics simulation, and kinetics study

The emergence of artemisinin-resistant variants of Plasmodium falciparum necessitates the urgent search for novel antimalarial drugs. In this regard, an in silico study to screen antimalarial drug candidates from a series of benzimidazole-thiosemicarbazone hybrid molecules with interesting antiplasmodial properties and explore their falcipain-2 (FP2) inhibitory potentials has been undertaken herein. FP2 is a key cysteine protease that degrades hemoglobin in Plasmodium falciparum and is an important biomolecular target in the development of antimalarial drugs. Pharmacokinetic properties, ADMET profiles, MM/GBSA-based binding free energies, reaction mechanisms, and associated barrier heights have been investigated. DFT, molecular dynamics simulation, molecular docking, and ONIOM methods were used. From the results obtained, four 4N-substituted derivatives of the hybrid molecule (E)-2-(1-(5-chloro-1H-benzo[d]imidazol-2-yl)ethylidene)hydrazine-1-carbothioamide (1A) denoted 1B, 1C, 1D, and 1E are drug-like and promising inhibitors of FP2, exhibiting remarkably small inhibitory constants (5.94 × 10-14 - 2.59 × 10-04 n M) and favorable binding free energies (-30.32 to -17.17 kcal/mol). Moreover, the ONIOM results have revealed that 1B and possibly 1C and 1D may act as covalent inhibitors of FP2. The rate-determining step of the thermodynamically favorable covalent binding mechanism occurs across a surmountable barrier height of 24.18 kcal/mol in water and 28.42 kcal/mol in diethyl ether. Our findings are useful for further experimental investigations on the antimalarial activities of the hybrid molecules studied.

2,4-Dihydroxycinnamic acid as spike ACE2 inhibitor and apigenin as RdRp inhibitor in Nimbamritadi Panchatiktam Kashayam against COVID-19: an in silico and in vitro approach

Nimbamritadi Panchatiktam Kashayam (NPK) is an ayurvedic formulation composed of ingredients with potent anti-viral activities. We studied the interaction energy of 144 phytoconstituents present in NPK against spike receptor-binding domain (RBD) complexed with ACE2 protein (PDB ID: 6LZG) and RNA-dependent RNA polymerase protein (PDB ID: 7BTF) using Biovia Drug Discovery studio. The result indicated that 2,4-hydroxycinnamic acid exerts more significant binding affinities (28.43 kcal/mol) than Umifenovir (21.24 kcal/mol) against spike ACE2. Apigenin exhibited the highest binding affinities (54.63 kcal/mol) compared with Remdesivir (24.52 kcal/mol) against RdRp. An in vitro analysis showed a reduction in the number of lentiviral particles on transfected HEK293T-hACE2 cells as assessed by pseudovirus inhibition assay. At the same time, the tested compounds showed non-toxic up to 100 µg/ml in normal cells by MTT assay. The study highlights the plausible clinical utility of this traditional medicine against SARS CoV2.

Structure Activity Relationship and Molecular Docking of Some Quinazolines Bearing Sulfamerazine Moiety as New 3CLpro, cPLA2, sPLA2 Inhibitors

The current work was conducted to synthesize several novel anti-inflammatory quinazolines having sulfamerazine moieties as new 3CLpro, cPLA2, and sPLA2 inhibitors. The thioureido derivative 3 was formed when compound 2 was treated with sulfamerazine. Also, compound 3 was reacted with NH2-NH2 in ethanol to produce the N-aminoquinazoline derivative. Additionally, derivative 4 was reacted with 4-hydroxy-3-methoxybenzaldehyde, ethyl chloroacetate, and/or diethyl oxalate to produce quinazoline derivatives 5, 6, and 12, respectively. The results of the pharmacological study indicated that the synthesized 4-6 and 12 derivatives showed good 3CLpro, cPLA2, and sPLA2 inhibitory activity. The IC50 values of the target compounds 4-6, and 12 against the SARS-CoV-2 main protease were 2.012, 3.68, 1.18, and 5.47 µM, respectively, whereas those of baicalein and ivermectin were 1.72 and 42.39 µM, respectively. The IC50 values of the target compounds 4-6, and 12 against sPLA2 were 2.84, 2.73, 1.016, and 4.45 µM, respectively, whereas those of baicalein and ivermectin were 0.89 and 109.6 µM, respectively. The IC50 values of the target compounds 4-6, and 12 against cPLA2 were 1.44, 2.08, 0.5, and 2.39 µM, respectively, whereas those of baicalein and ivermectin were 3.88 and 138.0 µM, respectively. Also, incubation of lung cells with LPS plus derivatives 4-6, and 12 caused a significant decrease in levels of sPLA2, cPLA2, IL-8, TNF-α, and NO. The inhibitory activity of the synthesized compounds was more pronounced compared to baicalein and ivermectin. In contrast to ivermectin and baicalein, bioinformatics investigations were carried out to establish the possible binding interactions between the newly synthesized compounds 2-6 and 12 and the active site of 3CLpro. Docking simulations were utilized to identify the binding affinity and binding mode of compounds 2-6 and 12 with the active sites of 3CLpro, sPLA2, and cPLA2 enzymes. Our findings demonstrated that all compounds had outstanding binding affinities, especially with the key amino acids of the target enzymes. These findings imply that compound 6 is a potential lead for the development of more effective SARS-CoV-2 Mpro inhibitors and anti-COVID-19 quinazoline derivative-based drugs. Compound 6 was shown to have more antiviral activity than baicalein and against 3CLpro. Furthermore, the IC50 value of ivermectin against the SARS-CoV-2 main protease was revealed to be 42.39 µM, indicating that it has low effectiveness.

A Drug Repurposing Approach to Identify Therapeutics by Screening Pathogen Box Exploiting SARS-CoV-2 Main Protease

Coronavirus disease-19 (COVID-19) is caused by severe acute respiratory syndrome coronavirus -2 (SARS-CoV-2) and is responsible for a higher degree of morbidity and mortality worldwide. There is a smaller number of approved therapeutics available to target the SARS-CoV-2 virus, and the virus is evolving at a fast pace. So, there is a continuous need for new therapeutics to combat COVID-19. The main protease (Mpro ) enzyme of SARS-CoV-2 is essential for replication and transcription of the viral genome, thus could be a potent target for the treatment of COVID-19. In the present study, we performed an in-silico screening analysis of 400 diverse bioactive inhibitors with proven antibacterial and antiviral properties against Mpro drug target. Ten compounds showed a higher binding affinity for Mpro than the reference compound (N3), with desired physicochemical properties. Furthermore, in-depth docking and superimposition revealed that three compounds (MMV1782211, MMV1782220, and MMV1578574) are actively interacting with the catalytic domain of Mpro . In addition, the molecular dynamics simulation study showed a solid and stable interaction of MMV178221-Mpro complex compared to the other two molecules (MMV1782220, and MMV1578574). In line with this observation, MM/PBSA free energy calculation also demonstrated the highest binding free energy of -115.8 kJ/mol for MMV178221-Mpro compound. In conclusion, the present in silico analysis revealed MMV1782211 as a possible and potent molecule to target the Mpro and must be explored in vitro and in vivo to combat the COVID-19.

Multi-targeted molecular docking, pharmacokinetics, and drug-likeness evaluation of coumarin based compounds targeting proteins involved in development of COVID-19

COVID-19 is a progressing pandemic of coronavirus disease-2019, which had drowned the whole world in a deep sorrow sea. Uncountable deaths were extending the list of deaths every single day. The present research was aimed to study the multi-target interaction of coumarins against COVID-19 using molecular docking analysis. The structure of coumarin compounds was checked for ADME and Lipinski rule of five by using SwissADME, an online tool. SARS-CoV-2 proteins such as RdRp, PLpro, Mpro and spike protein were collected from the Protein Data Bank. The molecular docking study was performed in the PyRx tool, and the molecular interactions were visualised by Discovery Studio Visualizer. All the coumarin compounds used in the study were obeyed Lipinski’s rule of 5 without any violations. All the three designed derivatives of phenprocoumon, hymecromone, and psoralen were showed high binding affinity and prominent interactions with the drug target. The presence of –OH groups in the compound, His41, a catalytic dyad in Mpro, number of and the distance of hydrogen bond interactions with SARS-CoV-2 targets was accountable for the high binding attractions. The modified drug structures possess better binding efficacy towards at least three targets compared to their parent compounds. Further, molecular dynamic studies can be suggested to find the ligand–protein complex stability. The present study outcome reveals that the designed coumarins can be synthesised and examined as a potent inhibitory drug of SARS-CoV-2.

Efficient Synthesis of Fluorescent Coumarins and Phosphorous-Containing Coumarin-Type Heterocycles via Palladium Catalyzed Cross-Coupling Reactions

Quantum-chemical calculations on the spectral properties of some aryl substituted 3-phosphonocoumarins were performed, and the effect of the substituents in the aryl moiety was evaluated. The structures possessing promising fluorescent properties were successfully synthesized via Suzuki and Sonogashira cross-coupling. The synthetic protocol was also applied for the phosphorous chemoisomer of 3-phosphonocoumarin, 1,2-benzoxaphosphorin, and their carboxylate analogues. The optical properties of the arylated and alkynylated products were experimentally determined. The obtained quantum-chemical and experimental results give the possibility for a fine tuning of the optical properties of phosphorous-containing coumarin systems by altering the substituent at its C-6 position.

New Chemicals Suppressing SARS-CoV-2 Replication in Cell Culture

Candidates to being inhibitors of the main protease (Mpro) of SARS-CoV-2 were selected from the database of Voronezh State University using molecular modeling. The database contained approximately 19,000 compounds represented by more than 41,000 ligand conformers. These ligands were docked into Mpro using the SOL docking program. For one thousand ligands with best values of the SOL score, the protein–ligand binding enthalpy was calculated by the PM7 quantum-chemical method with the COSMO solvent model. Using the SOL score and the calculated protein–ligand binding enthalpies, eighteen compounds were selected for the experiments. Several of these inhibitors suppressed the replication of the coronavirus in cell culture, and we used the best three among them in the search for chemical analogs. Selection among analogs using the same procedure followed by experiments led to identification of seven inhibitors of the SARS-CoV-2 replication in cell culture with EC50 values at the micromolar level. The identified inhibitors belong to three chemical classes. The three inhibitors, 4,4-dimethyldithioquinoline derivatives, inhibit SARS-CoV-2 replication in Vero E6 cell culture just as effectively as the best published non-covalent inhibitors, and show low cytotoxicity. These results open up a possibility to develop antiviral drugs against the SARS-CoV-2 coronavirus.

Methodology-Centered Review of Molecular Modeling, Simulation, and Prediction of SARS-CoV-2

Despite tremendous efforts in the past two years, our understanding of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), virus-host interactions, immune response, virulence, transmission, and evolution is still very limited. This limitation calls for further in-depth investigation. Computational studies have become an indispensable component in combating coronavirus disease 2019 (COVID-19) due to their low cost, their efficiency, and the fact that they are free from safety and ethical constraints. Additionally, the mechanism that governs the global evolution and transmission of SARS-CoV-2 cannot be revealed from individual experiments and was discovered by integrating genotyping of massive viral sequences, biophysical modeling of protein-protein interactions, deep mutational data, deep learning, and advanced mathematics. There exists a tsunami of literature on the molecular modeling, simulations, and predictions of SARS-CoV-2 and related developments of drugs, vaccines, antibodies, and diagnostics. To provide readers with a quick update about this literature, we present a comprehensive and systematic methodology-centered review. Aspects such as molecular biophysics, bioinformatics, cheminformatics, machine learning, and mathematics are discussed. This review will be beneficial to researchers who are looking for ways to contribute to SARS-CoV-2 studies and those who are interested in the status of the field.

Pharmacokinetics of Nafamostat, a Potent Serine Protease Inhibitor, by a Novel LC-MS/MS Analysis

Nafamostat, a synthetic serine protease inhibitor, has been used for the treatment of inflammatory diseases such as pancreatitis. Recently, an increasing number of studies have shown the promising antiviral effects of nafamostat for the treatment of coronavirus disease-19 (COVID-19). This study aimed to develop a novel liquid chromatography–tandem mass spectrometry (LC-MS/MS) analysis and to characterize the pharmacokinetics of nafamostat in rats. Nafamostat in the rat plasma was extracted by solid phase extraction, and ¹³C₆-nafamostat was used as an internal standard. The quantification limit of nafamostat in the rat plasma was 0.5 ng/mL. The LC-MS/MS method was fully validated and applied to characterize the pharmacokinetics of nafamostat in rats. Following intravenous injection (2 mg/kg), nafamostat in the plasma showed a multiexponential decline with an average elimination half-life (t_1/2) of 1.39 h. Following oral administration of nafamostat solutions (20 mg/kg) in 10% dimethyl sulfoxide (DMSO) and in 10% DMSO with 10% Tween 80, nafamostat was rapidly absorbed, and the average oral bioavailability was 0.95% and 1.59%, respectively. The LC-MS/MS method and the pharmacokinetic information of nafamostat could be helpful for the further preclinical and clinical studies of nafamostat.

Insights into the binding and covalent inhibition mechanism of PF-07321332 to SARS-CoV-2 Mpro

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been causing the COVID-19 pandemic, resulting in several million deaths being reported. Numerous investigations have been carried out to discover a compound that can inhibit the biological activity of the SARS-CoV-2 main protease, which is an enzyme related to the viral replication. Among these, PF-07321332 (Nirmatrelvir) is currently under clinical trials for COVID-19 therapy. Therefore, in this work, atomistic and electronic simulations were performed to unravel the binding and covalent inhibition mechanism of the compound to M^pro. Initially, 5 μs of steered-molecular dynamics simulations were carried out to evaluate the ligand-binding process to SARS-CoV-2 M^pro. The successfully generated bound state between the two molecules showed the important role of the PF-07321332 pyrrolidinyl group and the residues Glu166 and Gln189 in the ligand-binding process. Moreover, from the MD-refined structure, quantum mechanics/molecular mechanics (QM/MM) calculations were carried out to unravel the reaction mechanism for the formation of the thioimidate product from SARS-CoV-2 M^pro and the PF-07321332 inhibitor. We found that the catalytic triad Cys145–His41–Asp187 of SARS-CoV-2 M^pro plays an important role in the activation of the PF-07321332 covalent inhibitor, which renders the deprotonation of Cys145 and, thus, facilitates further reaction. Our results are definitely beneficial for a better understanding of the inhibition mechanism and designing new effective inhibitors for SARS-CoV-2 M^pro.

The catalytic triad Cys145–His41–Asp187 of SARS-CoV-2 M^pro plays an important role in the activation of the PF-07321332 covalent inhibitor.

Molecular docking and dynamic simulations of Cefixime, Etoposide and Nebrodenside A against the pathogenic proteins of SARS-CoV-2

The catastrophe of the coronavirus continues from one part of the world to another, and hardly a country is left without its devastations. Millions of people were infected and several hundred thousand died of the COVID-19 pandemic across the world. There is no clear targeted drug therapy available for the treatment of the patients. The discovery of vaccines is not enough to curtail its spread and disastrous implications. An instantly qualifying approach is needed to utilize the current drugs and isolated compounds. The purpose of this work is to determine potent inhibitors against the target proteins of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). For this purpose, molecular docking study of pathogenic spike glycoproteins (S), nucleocapsid phosphoprotein (N), an envelope protein (E), two drugs i.e., cefixime, etoposide, and a previously isolated compound nebrodenside A is performed. Promising results were obtained via complimentary analysis of molecular dynamics (MD) simulations performed for the complexes of three proteins with etoposide drug. Minimum values were recorded for the docking scores and binding energies of the complexes. These results were further supported by the RMSD, RMSF data for the stability of proteins and ligands. Additionally, ligand properties and ligand-protein contacts were also explained with histograms of every simulation trajectory. The computational studies confirmed that cefixime, etoposide, and nebrodenoside A can be used as potent inhibitors of COVID-19. Nevertheless, additional experimental investigations and validation of the selected candidates are mandatory to confirm their applicability for clinical trials.

Searching and designing potential inhibitors for SARS-CoV-2 Mpro from natural sources using atomistic and deep-learning calculations

The spread of severe acute respiratory syndrome coronavirus 2 novel coronavirus (SARS-CoV-2) worldwide has caused the coronavirus disease 2019 (COVID-19) pandemic. A hundred million people were infected, resulting in several millions of death worldwide. In order to prevent viral replication, scientists have been aiming to prevent the biological activity of the SARS-CoV-2 main protease (3CL pro or Mpro). In this work, we demonstrate that using a reasonable combination of deep-learning calculations and atomistic simulations could lead to a new approach for developing SARS-CoV-2 main protease (Mpro) inhibitors. Initially, the binding affinities of the natural compounds to SARS-CoV-2 Mpro were estimated via atomistic simulations. The compound tomatine, thevetine, and tribuloside could bind to SARS-CoV-2 Mpro with nanomolar/high-nanomolar affinities. Secondly, the deep-learning (DL) calculations were performed to chemically alter the top-lead natural compounds to improve ligand-binding affinity. The obtained results were then validated by free energy calculations using atomistic simulations. The outcome of the research will probably boost COVID-19 therapy.

The hybrid DeepFrag/atomistic simulation approach could lead to a new scheme for developing SARS-CoV-2 3CLpro/Mpro inhibitors.

Migration of Excitation Energy in Furocoumarins

The migration of excitation energy of a number of psoralen compounds has been studied. For this, the methods of induced absorption spectroscopy, stationary electron spectroscopy, fluorescence and phosphorescence, as well as quantum chemistry were used. A comparative photostability of psoralen was achieved by exposure to a XeCl excilamp irradiation (emission wavelength λ_em = 308 nm) with parameters Δλ = 5-10 nm, W_peak = 18 mW/cm², p = 8.1 J/cm³, f = 200 kHz, pulse duration 1 μs. It was found that the singlet-triplet transition played a major role in the migration of excitation energy into triplet states. Among all tested compounds, substances with an OCH₃-group in the structure have the strongest effect on the spectral-luminescent characteristics.

Non-conventional interactions of N3 inhibitor with the main protease of SARS-CoV and SARS-CoV-2

The extensive spread of COVID-19 in every continent shows that SARS-CoV-2 virus has a higher transmission rate than SARS-CoV virus which emerged in 2002. This results in a global pandemic that is difficult to control. In this investigation, we analyze the interaction of N3 inhibitor and the main protease of SARS-CoV and SARS-CoV-2 by quantum chemistry calculations. Non-covalent interactions involved in these systems were studied using a model of 469 atoms. Density Functional Theory and Quantum Theory of Atoms in Molecules calculations lead us to the conclusion that non-conventional hydrogen bonds are important to describe attractive interactions in these complexes. The energy of these non-conventional hydrogen bonds represents more than a half of the estimated interaction energy for non-covalent contacts. This means that hydrogen bonds are crucial to correctly describe the bonds between inhibitors and the main proteases. These results could be useful for the design of new drugs, since non-covalent interactions are related to possible mechanisms of action of molecules used against these viruses.

Main interactions of dopamine and risperidone with the dopamine D2 receptor

Psychosis is one of the psychiatric disorders that is controlled by dopaminergic drugs such as antipsychotics that have affinity for the dopamine D2 receptor (DRD2). In this investigation we perform quantum chemical calculations of two molecules [dopamine and risperidone] within a large cavity of DRD2 that represents the binding site of the receptor. Dopamine is an endogenous neurotransmitter and risperidone is a second-generation antipsychotic. Non-covalent interactions of dopamine and risperidone with DRD2 are analyzed using the Quantum Theory of Atoms in Molecules (QTAIM) and the Non-Covalent Interaction index (NCI). The QTAIM results show that these molecules strongly interact with the receptor. There are 22 non-covalent interactions for dopamine and 54 for risperidone. The electron density evaluated at each critical binding point is small in both systems but it is higher for dopamine than for risperidone, indicating that the interactions of DRD2 with the first are stronger than with the second molecule. However, the binding energy is higher for risperidone (-72.6 kcal mol-1) than for dopamine (-22.8 kcal mol-1). Thus, the strength of the binding energy is due to the number of contacts rather than the strength of the interactions themselves. This could be related to the ability of risperidone to block DRD2 and may explain the efficacy of this drug for controlling the symptoms of schizophrenia, but likewise its secondary effects.

Theobroma cacao L. compounds: Theoretical study and molecular modeling as inhibitors of main SARS-CoV-2 protease

Cocoa beans contain antioxidant molecules with the potential to inhibit type 2 coronavirus (SARS-CoV-2), which causes a severe acute respiratory syndrome (COVID-19). In particular, protease. Therefore, using in silico tests, 30 molecules obtained from cocoa were evaluated. Using molecular docking and quantum mechanics calculations, the chemical properties and binding efficiency of each ligand was evaluated, which allowed the selection of 5 compounds of this series. The ability of amentoflavone, isorhoifolin, nicotiflorin, naringin and rutin to bind to the main viral protease was studied by means of free energy calculations and structural analysis performed from molecular dynamics simulations of the enzyme/inhibitor complex. Isorhoifolin and rutin stand out, presenting a more negative binding ΔG than the reference inhibitor N-[(5-methylisoxazol-3-yl)carbonyl]alanyl-l-valyl-N~1~-((1R,2Z)−4-(benzyloxy)−4-oxo-1-{[(3R)−2-oxopyrrolidin-3-yl]methyl}but-2-enyl)-L-leucinamide (N3). These results are consistent with high affinities of these molecules for the major SARS-CoV-2. The results presented in this paper are a solid starting point for future in vitro and in vivo experiments aiming to validate these molecules and /or test similar substances as inhibitors of SARS-CoV-2 protease.