Structure-based screening and validation of bioactive compounds as Zika virus methyltransferase (MTase) inhibitors through first-principle density functional theory, classical molecular simulation and QM/MM affinity estimation

Role of Alkannin in the Therapeutic Targeting of Protein-Tyrosine Phosphatase 1B and Aldose Reductase in Type 2 Diabetes: An In Silico and In Vitro Evaluation

Alkannin is a plant-derived naphthoquinone that is isolated from the Boraginaceae family plants. In our previous studies, we found that shikonin, which is the R-enantiomer of alkannin, has potent antidiabetic activity by inhibiting the action of the aldose reductase (AR) enzyme and the protein-tyrosine phosphatase 1B (PTP1B). Therefore, in this study, we aim to explore the antidiabetic effect of alkannin targeting PTP1B and AR by employing in silico and in vitro techniques. For in silico, we used different parameters such as ADMET analysis, molecular docking, MD simulation, Root Mean Square Deviation (RMSD), protein-ligand mapping, and free binding energy calculation. The in vitro evaluation was done by assessing the inhibitory activity and enzyme kinetics of PTP1B and AR inhibition by alkannin. The in silico studies indicate that alkannin possesses favorable pharmacological properties and possesses strong binding affinity for diabetes target proteins. Hydrogen bonds (Val297, Ala299, Leu300, and Ser302) and hydrophobic interactions (Trp20, Val47, Tyr48, Trp79, Trp111, Phe122, Trp219, Val297, Cys298, Ala299, Leu300, and Leu301) are established by the compound, which potentially improves specificity and aids in the stabilization of the protein-ligand complex. The results from in vitro studies show a potent dose-dependent PTP1B inhibitory activity with an IC50 value of 19.47 μM, and toward AR it was estimated at 22.77 μM. Thus, from the results it is concluded that a low IC50 value of alkannin for both PTP1B and AR along with favorable pharmacological properties and optimal intra-molecular interactions indicates its utilization as a potential drug candidate for the management of diabetes and its end complications.

Identification of Novel Non-Nucleoside Inhibitors of Zika Virus NS5 Protein Targeting MTase Activity

Zika virus (ZIKV) is a positive-sense single-stranded virus member of the Flaviviridae family. Among other arboviruses, ZIKV can cause neurological disorders such as Guillain Barré syndrome, and it can have congenital neurological manifestations and affect fertility. ZIKV nonstructural protein 5 (NS5) is essential for viral replication and limiting host immune detection. Herein, we performed virtual screening to identify novel small-molecule inhibitors of the ZIKV NS5 methyltransferase (MTase) domain. Compounds were tested against the MTases of both ZIKV and DENV, demonstrating good inhibitory activities against ZIKV MTase. Extensive molecular dynamic studies conducted on the series led us to identify other derivatives with improved activity against the MTase and limiting ZIKV infection with an increased selectivity index. Preliminary pharmacokinetic parameters have been determined, revealing excellent stability over time. Preliminary in vivo toxicity studies demonstrated that the hit compound 17 is well tolerated after acute administration. Our results provide the basis for further optimization studies on novel non-nucleoside MTase inhibitors.

Application of in-silico drug discovery techniques to discover a novel hit for target-specific inhibition of SARS-CoV-2 Mpro's revealed allosteric binding with MAO-B receptor: A theoretical study to find a cure for post-covid neurological disorder

Several studies have revealed that SARS-CoV-2 damages brain function and produces significant neurological disability. The SARS-CoV-2 coronavirus, which causes COVID-19, may infect the heart, kidneys, and brain. Recent research suggests that monoamine oxidase B (MAO-B) may be involved in metabolomics variations in delirium-prone individuals and severe SARS-CoV-2 infection. In light of this situation, we have employed a variety of computational to develop suitable QSAR model using PyDescriptor and genetic algorithm-multilinear regression (GA-MLR) models (R2 = 0.800-793, Q2LOO = 0.734-0.727, and so on) on the data set of 106 molecules whose anti-SARS-CoV-2 activity was empirically determined. QSAR models generated follow OECD standards and are predictive. QSAR model descriptors were also observed in x-ray-resolved structures. After developing a QSAR model, we did a QSAR-based virtual screening on an in-house database of 200 compounds and found a potential hit molecule. The new hit's docking score (-8.208 kcal/mol) and PIC50 (7.85 M) demonstrated a significant affinity for SARS-CoV-2's main protease. Based on post-covid neurodegenerative episodes in Alzheimer's and Parkinson's-like disorders and MAO-B's role in neurodegeneration, the initially disclosed hit for the SARS-CoV-2 main protease was repurposed against the MAO-B receptor using receptor-based molecular docking, which yielded a docking score of -12.0 kcal/mol. This shows that the compound that inhibits SARS-CoV-2's primary protease may bind allosterically to the MAO-B receptor. We then did molecular dynamic simulations and MMGBSA tests to confirm molecular docking analyses and quantify binding free energy. The drug-receptor complex was stable during the 150-ns MD simulation. The first computational effort to show in-silico inhibition of SARS-CoV-2 Mpro and allosteric interaction of novel inhibitors with MAO-B in post-covid neurodegenerative symptoms and other disorders. The current study seeks a novel compound that inhibits SAR's COV-2 Mpro and perhaps binds MAO-B allosterically. Thus, this study will enable scientists design a new SARS-CoV-2 Mpro that inhibits the MAO-B receptor to treat post-covid neurological illness.

Inhibition of Mycobacterium tuberculosis resuscitation-promoting factor B (RpfB) by microbially derived natural compounds: a computational study

The majority of the world population (around 25%) has latent Mycobacterium tuberculosis (Mtb) infection, among which only 5-10% of individuals develop active tuberculosis (TB), and 90-95% continue to have latent tuberculosis infection. This makes it the biggest global health concern. It has been reported that the resuscitation-promoting factor B (RpfB) is an exciting potential target for tuberculosis drug discovery due to its significant role in the reactivation of latent TB infection to an active infection. Several attempts have been made to investigate potential inhibitors against RpfB utilizing in-silico approaches. The present study also utilized a computational approach to investigate microbially derived natural compounds against the Mtb RpfB protein which is a very cost-effective This evaluation used structure-based virtual screening (SBVS), drug-likeness profiling, molecular docking, molecular dynamics simulation, and free-binding energy calculations. Six potential natural compounds, viz. Cyclizidine I, Boremexin C, Xenocoumacin 2, PM-94128, Cutinostatin B, and (+)1-O-demethylvariecolorquinone A were selected, which displayed a potential binding affinity between -52.39 and -60.87 Kcal/mol MMGBSA score and docking energy between -7.307 Kcal/mol to -6.972 Kcal/mol. All the complexes showed acceptable stability (<2.7 Å RMSD) during 100 ns MD simulation time except the RpfB protein-xenocoumacin 2 complex. This result exhibited that the selected compounds have high efficiency in inhibiting the Mtb RpfB and can be taken into account for additional in vitro and in vivo experimental validation.Communicated by Ramaswamy H. Sarma.

Active site determination of novel plant versatile peroxidase extracted from Citrus sinensis and bioconversion of β-naphthol

A ligninolytic peroxidase called versatile peroxidase, VP, (EC 1.11.1.16) is an iron-containing metalloenzyme. The most distinctive feature of this enzyme is its composite molecular framework, which combines lignin peroxidase's capacity to oxidize compounds with high-redox potential with manganese peroxidase's capacity to oxidize Mn2+ to Mn3+. In this study, we have extracted amino acid sequences from the Citrus sinensis source and subjected them to various computation tools to visualize the insight secondary and 3D structure, physicochemical properties, and validation of the structure which have not been studied so far to further investigate the catalytic efficiency and effectiveness of VP. The binding energies of HEME and HEME C (HEC) ligands with produced PDB (6rqf.1. A) have been also assessed, analyzed, and confirmed utilizing AutoDock. Binding energies were calculated using the AutoDock and validated by MD simulation using SCHRODINGER DESMOND. Most stable confirmation was achieved through a protein-ligand interaction study. Bio-technological use of VP in the biotransformation of β-naphthol has also been studied. The findings in the current study will have a substantial impact on proteomics, biochemistry, biotechnology, and possible uses of versatile peroxidase in the bio-remediation of different toxic organic compounds.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03758-x.

Selected phytochemicals of Momordica charantia L. as potential anti-DENV-2 through the docking, DFT and molecular dynamic simulation

Dengue fever is now one of the major global health concerns particularly for tropical and sub-tropical countries. However, there has been no FDA approved medication to treat dengue fever. Researchers are looking into DENV NS5 RdRp protease as a potential therapeutic target for discovering effective anti-dengue agents. The aim of this study to discover dengue virus inhibitor from a set of five compounds from Momordica charantia L. using a series of in-silico approaches. The compounds were docked into the active area of the DENV-2 NS5 RdRp protease to obtain the hit compounds. The successful compounds underwent additional testing for a study on drug-likeness similarity. Our study obtained Momordicoside-I as a lead compound which was further exposed to the Cytochrome P450 (CYP450) toxicity analysis to determine the toxicity based on docking scores and drug-likeness studies. Moreover, DFT studies were carried out to calculate the thermodynamic, molecular orbital and electrostatic potential properties for the lead compound. Moreover, the lead compound was next subjected to molecular dynamic simulation for 200 ns in order to confirm the stability of the docked complex and the binding posture discovered during docking experiment. Overall, the lead compound has demonstrated good medication like qualities, non-toxicity, and significant binding affinity towards the DENV-2 RdRp enzyme.Communicated by Ramaswamy H. Sarma.

Computational Prediction of 3,5-Diaryl-1H-Pyrazole and spiropyrazolines derivatives as potential acetylcholinesterase inhibitors for alzheimer disease treatment by 3D-QSAR, molecular docking, molecular dynamics simulation, and ADME-Tox

The efficacy of 40 synthesized variants of 3,5-diaryl-1H-pyrazole and spiropyrazoline' derivatives as acetylcholinesterase inhibitors is verified using a quantitative three-dimensional structure-activity relationship (3D-QSAR) by comparative molecular field analysis (CoMFA) and molecular similarity index analysis (CoMSIA) models. In this research, different field models proved that CoMSIA/SE model is the best model with high predictive power compared to several models (Qved2 = O.65; R2 = 0.980; R2test = 0.727). Also, contour maps produced by CoMSIA/SE model have been employed to prove the key structural needs of the activity. Consequently, six new compounds have been generated. Among these compounds, M4 and M5 were the most active but remained toxic and had poor absorption capacities. While the M1, M2, M3 and M6 remained highly active while respecting ADMET's characteristics. Molecular docking results showed compound M2 better with acetylcholinesterase than compound 22. The interactions are classical hydrogen bonding with residues TYR:124, TYR:72, and SER:293, which play a critical role in the biological activity as AChE inhibitors. MD results confirmed the docking results and showed that compound M2 had satisfactory stability with (ΔGbinding = -151.225 KJ/mol) in the active site of AChE receptor compared with compound 22 (ΔGbinding = -133.375 KJ/mol). In addition, both compounds had good stability regarding RMSD, Rg, and RMSF. The previous results show that the newly designed compound M2 is more active in the active site of AChE receptor than compound 22.Communicated by Ramaswamy H. Sarma.

Identification of promising methionine aminopeptidase enzyme inhibitors: A combine study of comprehensive virtual screening and dynamics simulation study

Methionine aminopeptidase (MetAP) enzymes play a critical role in bacterial cell survival by cleaving formyl-methionine initiators at N-terminal of nascent protein, a process which is vital in proper protein folding. This makes MetAP an attractive and novel antibacterial target to unveil promising antibiotics. In this study, the crystal structure of R. prowazekii MetAP was used in structure-based virtual screening of drug libraries such as Asinex antibacterial library and Comprehensive Marine Natural Products Database (CMNPD) to identify promising lead molecules against the enzyme. This shortlisted three drug molecules; BDE-25098678, BDE-30686468 and BDD_25351157 as most potent leads that showed strong binding to the MetAP enzyme. The static docked conformation of the compounds to the MetAP was reevaluated in molecular dynamics simulation studies. The analysis observed the docked complexes as stable structure with no major local or global deviations noticed. These findings suggest the formation of strong intermolecular docked complexes, which showed stable dynamics and atomic level interactions network. The binding free energy analysis predicted net MMGBSA energy of complexes as: BDE-25098678 (-73.41 kcal/mol), BDE-30686468 (-59.93 kcal/mol), and BDD_25351157 (-75.39 kcal/mol). In case of MMPBSA, the complexes net binding energy was as; BDE-25098678 (-77.47 kcal/mol), BDE-30686468 (-69.47 kcal/mol), and BDD_25351157 (-75.6 kcal/mol). Further, the compounds were predicted to follow the famous Lipinski rule of five and have non-toxic, non-carcinogenic and non-mutagenic profile. The screened compounds might be used in experimental test to highlight the real anti- R. prowazekii MetAP activity.

Recent advances in the development of methyltransferase (MTase) inhibitors against (re)emerging arboviruses diseases dengue and Zika

Emerging and/or re-emerging viral diseases such as dengue and Zika are a worldwide concern. Therefore, new antiviral therapeutics are necessary. In this sense, a non-structural protein with methyltransferase (MTase) activity is an attractive drug target because it plays a crucial role in dengue and Zika virus replication. Different drug strategies such as virtual screening, molecular docking, and molecular dynamics have identified new inhibitors that bind on the MTase active site. Therefore, in this review, we analyze MTase inhibitors, including S-adenosyl-L-methionine (SAM), S-adenosyl-l-homocysteine (SAH) and guanosine-5'-triphosphate (GTP) analogs, nitrogen-containing heterocycles (pyrimidine, adenosine, and pyridine), urea derivatives, and natural products. Advances in the design of MTase inhibitors could lead to the optimization of a possible single or broad-spectrum antiviral drug against dengue and Zika virus.

Synthesis, biological evaluation, and computational studies of some novel quinazoline derivatives as anticancer agents

A series of quinazolinone derivatives (7a-7h) were synthesized as antiproliferative agents. All compounds, were synthesized through three steps method and structurally evaluated by FTIR, ¹H-NMR, ¹³CNMR and Mass spectroscopy. Their cytotoxic activities were assessed using MTT protocol against three humans cancerous (MCF-7, A549 and 5637) and normal (MRC-5) cell lines. In addition, molecular docking and simulation studies of the synthesized compounds were performed to assessment their orientation, interaction mode against EGFR as plausible mechanism of quinazoline compounds as anticancer agents. The synthesized compounds mostly showed moderate activity against the three studied cell lines. They also indicated an appropriate selectivity against tumorigenic and non-tumorigenic cell line. The molecular docking results also confirmed biological activity. Most of the compounds fulfilled Lipinski rule. Collectively, these compounds with further modification can be considered as potent antiproliferative agents.

Three-dimensional quantitative structural-activity relationship and molecular dynamics study of multivariate substituted 4-oxyquinazoline HDAC6 inhibitors

3D-QSAR models were established by collecting 46 multivariate-substituted 4-oxyquinazoline HDAC6 inhibitors. The relationship of molecular structure and inhibitory activity was studied by comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA). The results showed the models established by CoMFA (q² = 0.590, r² = 0.965) and CoMSIA (q² = 0.594, r² = 0.931) had good prediction ability. At the same time, 3D-QSAR models met the internal verification, external verification and AD test. Ten new compounds were designed based on CoMFA and CoMSIA contour maps and their pharmacokinetic/toxic properties (ADME/T) were evaluated. It was found that most compounds have well safety profile and pharmacokinetic property. Then, we explored the interaction between HDAC6 and compounds by molecular docking. The results showed that the binding mode of the new compounds with HDAC6 was the same as the template compound 46, and the hydrogen bond and hydrophobic bond played a vital role in the binding process. Molecular dynamics simulation results showed that residues Ser531, His574 and Tyr745 played key roles in the binding process. All newly designed compounds had lower energy gap and binding energy than compound 46 according to DFT analysis and free energy analysis. This study provided a theoretical reference for designing compounds of higher activity and a new idea for the development of novel HDAC6 inhibitors.

Review on development of potential inhibitors of SARS-CoV-2 main protease (MPro)

Background

The etiological agent for the coronavirus illness outbreak in 2019-2020 is a novel coronavirus known as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) (COVID-19), whereas coronavirus disease pandemic of 2019 (COVID-19) has compelled the implementation of novel therapeutic options.

Main body of the abstract

There are currently no targeted therapeutic medicines for this condition, and effective treatment options are quite restricted; however, new therapeutic candidates targeting the viral replication cycle are being investigated. The primary protease of the severe acute respiratory syndrome coronavirus 2 virus is a major target for therapeutic development (M^Pro). Severe acute respiratory syndrome coronavirus 2, severe acute respiratory syndrome coronavirus, and Middle East respiratory syndrome coronavirus (MERS-CoV) all seem to have a structurally conserved substrate-binding domain that can be used to develop novel protease inhibitors.

Short conclusion

With the recent publication of the X-ray crystal structure of the severe acute respiratory syndrome coronavirus 2 Mm, virtual and in vitro screening investigations to find M^Pro inhibitors are fast progressing. The focus of this review is on recent advancements in the quest for small-molecule inhibitors of the severe acute respiratory syndrome coronavirus 2 main protease.

Efficient synthesis of 1,3-naphtoxazine derivatives using reusable magnetic catalyst (GO-Fe3O4-Ti(IV)): anticonvulsant evaluation and computational studies

A series of 2-aryl/alkyl-2,3-dihydro-1H-naphtho[1,2-e][1,3]oxazines (S₁-S₁₁) were synthesized with an eco-friendly and recoverable nanocatalyst (GO-Fe₃O₄-Ti^(IV)) as an efficient magnetic composite. The new nanocatalyst was characterized by FT-IR, XRD and, EDS analysis. A conformable procedure, easy to work up and having a short reaction time with high yields are some advantages of this method. The new catalyst is also thermal-stable, reusable and, environment-friendly. The chemical structures of the synthesized 1,3-oxazine compounds were confirmed by comparing their melting points with those reported in literature. Then, the anticonvulsant activity of these compounds was assessed by the intraperitoneal pentylenetetrazole test (ipPTZ). Compounds S₁₀ and S₁₁ displayed considerable activity against chemically-induced seizure tests. The molecular simulation was also done to achieve their binding affinities as γ-aminobutyric acid A (GABA-A) receptor agonists as an assumptive mechanism of their anticonvulsant action. The result of molecular studies represented strongly matched with biological activity. Molecular docking simulation of the potent compound (S₁₀) and diazepam as the positive control was performed and some critical residues like Thr262, Asn265, Met286, Phe289, and Val290 were identified. Based on the anticonvulsant results and also in silico ADME predictions, S₁₁ can be to become a potential drug candidate as an anticonvulsant agent.

A review on structural genomics approach applied for drug discovery against three vector-borne viral diseases: Dengue, Chikungunya and Zika

Structural genomics involves the advent of three-dimensional structures of the genome encoded proteins through various techniques available. Numerous structural genomics research groups have been developed across the globe and they contribute enormously to the identification of three-dimensional structures of various proteins. In this review, we have discussed the applications of the structural genomics approach towards the discovery of potential lead-like molecules against the genomic drug targets of three vector-borne diseases, namely, Dengue, Chikungunya and Zika. Currently, all these three diseases are associated with the most important global public health problems and significant economic burden in tropical countries. Structural genomics has accelerated the identification of novel drug targets and inhibitors for the treatment of these diseases. We start with the current development status of the drug targets and antiviral drugs against these three diseases and conclude by describing challenges that need to be addressed to overcome the shortcomings in the process of drug discovery.

A multi-targeted approach to identify potential flavonoids against three targets in the SARS-CoV-2 life cycle

The advent and persistence of the Severe Acute Respiratory Syndrome Coronavirus - 2 (SARS-CoV-2)-induced Coronavirus Disease (COVID-19) pandemic since December 2019 has created the largest public health emergency in over a century. Despite the administration of multiple vaccines across the globe, there continues to be a lack of approved efficacious non-prophylactic interventions for the disease. Flavonoids are a class of phytochemicals with historically established antiviral, anti-inflammatory and antioxidative properties that are effective against cancers, type 2 diabetes mellitus, and even other human coronaviruses. To identify the most promising bioactive flavonoids against the SARS-CoV-2, this article screened a virtual library of 46 bioactive flavonoids against three promising targets in the SARS-CoV-2 life cycle: human TMPRSS2 protein, 3CLpro, and PLpro. By examining the effects of glycosylation and other structural-activity relationships, the presence of sugar moiety in flavonoids significantly reduces its binding energy. It increases the solubility of flavonoids leading to reduced toxicity and higher bioavailability. Through protein-ligand contact profiling, it was concluded that naringin formed more hydrogen bonds with TMPRSS2 and 3CLpro. In contrast, hesperidin formed a more significant number of hydrogen bonds with PLpro. These observations were complimented by the 100 ns molecular dynamics simulation and binding free energy analysis, which showed a considerable stability of docked bioflavonoids in the active site of SARS-CoV-2 target proteins. Finally, the binding affinity and stability of the selected docked complexes were compared with the reference ligands (camostat for TMPRSS2, GC376 for 3CLpro, and GRL0617 for PLpro) that strongly inhibit their respective SARS-COV-2 targets. Overall analysis revealed that the selected flavonoids could be potential therapeutic agents against SARS-CoV-2. Naringin showed better affinity and stability for TMPRSS2 and 3CLpro, whereas hesperidin showed a better binding relationship and stability for PLpro.

De novo design of novel protease inhibitor candidates in the treatment of SARS-CoV-2 using deep learning, docking, and molecular dynamic simulations

The main protease of SARS-CoV-2 is a critical target for the design and development of antiviral drugs. 2.5 M compounds were used in this study to train an LSTM generative network via transfer learning in order to identify the four best candidates capable of inhibiting the main proteases in SARS-CoV-2. The network was fine-tuned over ten generations, with each generation resulting in higher binding affinity scores. The binding affinities and interactions between the selected candidates and the SARS-CoV-2 main protease are predicted using a molecular docking simulation using AutoDock Vina. The compounds selected have a strong interaction with the key MET 165 and Cys145 residues. Molecular dynamics (MD) simulations were run for 150ns to validate the docking results on the top four ligands. Additionally, root-mean-square deviation (RMSD), root-mean-square fluctuation (RMSF), and hydrogen bond analysis strongly support these findings. Furthermore, the MM-PBSA free energy calculations revealed that these chemical molecules have stable and favorable energies, resulting in a strong binding with Mpro's binding site. This study's extensive computational and statistical analyses indicate that the selected candidates may be used as potential inhibitors against the SARS-CoV-2 in-silico environment. However, additional in-vitro, in-vivo, and clinical trials are required to demonstrate their true efficacy.

Mechanistic insights into the Japanese encephalitis virus RNA dependent RNA polymerase protein inhibition by bioflavonoids from Azadirachta indica

Japanese encephalitis (JE) virus is a flavivirus causing encephalitis causing neurological damage. RNA-dependent-RNA-polymerase (RdRp) is responsible for genome replication making it excellent anti-viral target. In this study, the crystal structure of JE RdRp (jRdRp) and bioflavonoids reported in Azadirachta indica were retrieved from specific databases. Structure-based virtual screening was employed using MTiOpenScreen server and top four compounds selected with the most negative docking scores. Conformations were redocked using AutoDock Vina; these complexes showed mechanistic interactions with Arg⁴⁷⁴, Gly⁶⁰⁵, Asp⁶⁶⁸, and Trp⁸⁰⁰ residues in the active site of jRdRp, i.e., guanosine-5′-triphosphate. Furthermore, 100 ns classical molecular dynamics simulation and binding free energy calculation showed stability of docked bioflavonoids in the active jRdRp pocket and significant contribution of van-der-Waals interactions for docked complex stability during simulation. Therefore, this study predicted the anti-viral activity of Gedunin, Nimbolide, Ohchinin acetate, and Kulactone against jRdRp and can be considered for further antiviral drug development.

The identification of novel inhibitors of human angiotensin-converting enzyme 2 and main protease of Sars-Cov-2: A combination of in silico methods for treatment of COVID-19

The angiotensin-converting enzyme 2 (ACE2) and main protease (MPro), are the putative drug candidates for the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). In this study, we performed 3D-QSAR pharmacophore modeling and screened 1,264,479 ligands gathered from Pubchem and Zinc databases. Following the calculation of the ADMET properties, molecular docking was carried out. Moreover, the de novo ligand design was performed. MD simulation was then applied to survey the behavior of the ligand-protein complexes. Furthermore, MMPBSA has utilized to re-estimate the binding affinities. Then, a free energy landscape was used to find the most stable conformation of the complexes. Finally, the hybrid QM-MM method was carried out for the precise calculation of the energies. The Hypo1 pharmacophore model was selected as the best model. Our docking results indicate that the compounds ZINC12562757 and 112,260,215 were the best potential inhibitors of the ACE2 and MPro, respectively. Furthermore, the Evo_1 compound enjoys the highest docking energy among the designed de novo ligands. Results of RMSD, RMSF, H-bond, and DSSP analyses have demonstrated that the lead compounds preserve the stability of the complexes' conformation during the MD simulation. MMPBSA data confirmed the molecular docking results. The results of QM-MM showed that Evo_1 has a stronger potential for specific inhibition of MPro, as compared to the 112,260,215 compound.

In silico structural elucidation of the rabies RNA-dependent RNA polymerase (RdRp) toward the identification of potential rabies virus inhibitors

The rabies virus (RABV) is a non-segmented, negative single-stranded RNA virus which causes acute infection of the central nervous system in humans. Once symptoms appear, the result is nearly always death, and to date, post-exposure prophylaxis (PEP) is the only treatment applicable only immediately after an exposure. Previous studies have identified viral RNA-dependent RNA polymerase (RdRp) as a potential drug target due to its significant role in viral replication and transcription. Herein we generated an energy-minimized homology model of RABIES-RdRp and used it for virtual screening against 2045 NCI Diversity Set III library. The best five ligand-RdRp complexes were picked for further energy minimization via molecular dynamics (MDs) where the complex with ligand Z01690699 shows a minimum score characterized with stable hydrogen bonds and hydrophobic interactions with the catalytic site residues. Our study identified an important ligand for development of remedial approach for treatment of rabies infection.

Structure-Based Identification of Natural Products as SARS-CoV-2 Mpro Antagonist from Echinacea angustifolia Using Computational Approaches

Coronavirus disease-19 (COVID-19) pandemic, caused by the novel SARS-CoV-2 virus, continues to be a global threat. The number of cases and deaths will remain escalating due to the lack of effective therapeutic agents. Several studies have established the importance of the viral main protease (Mpro) in the replication of SARS-CoV-2 which makes it an attractive target for antiviral drug development, including pharmaceutical repurposing and other medicinal chemistry approaches. Identification of natural products with considerable inhibitory potential against SARS-CoV-2 could be beneficial as a rapid and potent alternative with drug-likeness by comparison to de novo antiviral drug discovery approaches. Thereof, we carried out the structure-based screening of natural products from Echinacea-angustifolia, commonly used to prevent cold and other microbial respiratory infections, targeting SARS-CoV-2 Mpro. Four natural products namely, Echinacoside, Quercetagetin 7-glucoside, Levan N, Inulin from chicory, and 1,3-Dicaffeoylquinic acid, revealed significant docking energy (>-10 kcal/mol) in the SARS-CoV-2 Mpro catalytic pocket via substantial intermolecular contacts formation against co-crystallized ligand (<-4 kcal/mol). Furthermore, the docked poses of SARS-CoV-2 Mpro with selected natural products showed conformational stability through molecular dynamics. Exploring the end-point net binding energy exhibited substantial contribution of Coulomb and van der Waals interactions to the stability of respective docked conformations. These results advocated the natural products from Echinacea angustifolia for further experimental studies with an elevated probability to discover the potent SARS-CoV-2 Mpro antagonist with higher affinity and drug-likeness.

Computational and In Vitro Investigation of (-)-Epicatechin and Proanthocyanidin B2 as Inhibitors of Human Matrix Metalloproteinase 1

Matrix metalloproteinases 1 (MMP-1) energetically triggers the enzymatic proteolysis of extracellular matrix collagenase (ECM), resulting in progressive skin aging. Natural flavonoids are well known for their antioxidant properties and have been evaluated for inhibition of matrix metalloproteins in human. Recently, (-)-epicatechin and proanthocyanidin B2 were reported as essential flavanols from various natural reservoirs as potential anti-inflammatory and free radical scavengers. However, their molecular interactions and inhibitory potential against MMP-1 are not yet well studied. In this study, sequential absorption, distribution, metabolism, and excretion (ADME) profiling, quantum mechanics calculations, and molecular docking simulations by extra precision Glide protocol predicted the drug-likeness of (-)-epicatechin (−7.862 kcal/mol) and proanthocyanidin B2 (−8.145 kcal/mol) with the least reactivity and substantial binding affinity in the catalytic pocket of human MMP-1 by comparison to reference bioactive compound epigallocatechin gallate (−6.488 kcal/mol). These flavanols in docked complexes with MMP-1 were further studied by 500 ns molecular dynamics simulations that revealed substantial stability and intermolecular interactions, viz. hydrogen and ionic interactions, with essential residues, i.e., His218, Glu219, His222, and His228, in the active pocket of MMP-1. In addition, binding free energy calculations using the Molecular Mechanics Generalized Born Surface Area (MM/GBSA) method suggested the significant role of Coulomb interactions and van der Waals forces in the stability of respective docked MMP-1-flavonol complexes by comparison to MMP-1-epigallocatechin gallate; these observations were further supported by MMP-1 inhibition assay using zymography. Altogether with computational and MMP-1–zymography results, our findings support (-)-epicatechin as a comparatively strong inhibitor of human MMP-1 with considerable drug-likeness against proanthocyanidin B2 in reference to epigallocatechin gallate.

Computational insights into tetracyclines as inhibitors against SARS-CoV-2 Mpro via combinatorial molecular simulation calculations

The COVID-19 pandemic raised by SARS-CoV-2 is a public health emergency. However, lack of antiviral drugs and vaccine against human coronaviruses demands a concerted approach to challenge the SARS-CoV-2 infection. Under limited resource and urgency, combinatorial computational approaches to identify the potential inhibitor from known drugs could be applied against risen COVID-19 pandemic. Thereof, this study attempted to purpose the potent inhibitors from the approved drug pool against SARS-CoV-2 main protease (M^pro). To circumvent the issue of lead compound from available drugs as antivirals, antibiotics with broad spectrum of viral activity, i.e. doxycycline, tetracycline, demeclocycline, and minocycline were chosen for molecular simulation analysis against native ligand N3 inhibitor in SARS-CoV-2 M^pro crystal structure. Molecular docking simulation predicted the docking score >−7 kcal/mol with significant intermolecular interaction at the catalytic dyad (His41 and Cys145) and other essential substrate binding residues of SARS-CoV-2 M^pro. The best ligand conformations were further studied for complex stability and intermolecular interaction profiling with respect to time under 100 ns classical molecular dynamics simulation, established the significant stability and interactions of selected antibiotics by comparison to N3 inhibitor. Based on combinatorial molecular simulation analysis, doxycycline and minocycline were selected as potent inhibitor against SARS-CoV-2 M^pro which can used in combinational therapy against SARS-CoV-2 infection.