In silico screening of anti-inflammatory compounds from Lichen by targeting cyclooxygenase-2

Multi-Target In-Silico modeling strategies to discover novel angiotensin converting enzyme and neprilysin dual inhibitors

Cardiovascular diseases, including heart failure, stroke, and hypertension, affect 608 million people worldwide and cause 32% of deaths. Combination therapy is required in 60% of patients, involving concurrent Renin-Angiotensin-Aldosterone-System (RAAS) and Neprilysin inhibition. This study introduces a novel multi-target in-silico modeling technique (mt-QSAR) to evaluate the inhibitory potential against Neprilysin and Angiotensin-converting enzymes. Using both linear (GA-LDA) and non-linear (RF) algorithms, mt-QSAR classification models were developed using 983 chemicals to predict inhibitory effects on Neprilysin and Angiotensin-converting enzymes. The Box-Jenkins method, feature selection method, and machine learning algorithms were employed to obtain the most predictive model with ~ 90% overall accuracy. Additionally, the study employed virtual screening of designed scaffolds (Chalcone and its analogues, 1,3-Thiazole, 1,3,4-Thiadiazole) applying developed mt-QSAR models and molecular docking. The identified virtual hits underwent successive filtration steps, incorporating assessments of drug-likeness, ADMET profiles, and synthetic accessibility tools. Finally, Molecular dynamic simulations were then used to identify and rank the most favourable compounds. The data acquired from this study may provide crucial direction for the identification of new multi-targeted cardiovascular inhibitors.

Computational identification of new TKI as potential noncovalent reversible EGFRL858R/T790M inhibitors: VHTS, molecular docking, DFT study and molecular dynamic simulation

The mutations concerned with non-small cell lung cancer involving epidermal growth factor receptor of tyrosine kinase family have primarily targeted. In this study, we employed a scalable high-throughput virtual screening (HTVS) framework and a targeted compound library of over 50.000 Erlotinib-derived compounds as noncovalent reversible EGFRL858R/T790M inhibitors. Our HTVS work flow leverages include HTVS, SP (Standard Precision) and XP (Extra Precision) docking protocol along with its relative binding free energy calculation, cluster analysis study and ADMET properties. Then we used multiple ns-time scale molecular dynamics (MD) simulations and density functional theory (DFT) precise calculation techniques to elucidate how the bound ligand interact with the complexes conformational states involving motions both proximal and distal to the binding site. Based on glide score and protein-ligand interactions, the highest scoring molecule was selected for molecular dynamic simulation providing a complete insight into the conformational stability. A hyperfine analysis of DFT based refinement strategy highly supported their stability by strong intermolecular interactions. Together, our results demonstrate that the virtually screened top retained molecules present the best moieties introduced to Erlotinib. They exhibit interesting pharmacokinetic properties that can act as potent antitumor drug candidates than the lead compound drug and in some extent tackling the drug resistance problem which offer a springboard for further therapeutic experiments and applications.Communicated by Ramaswamy H. Sarma.

Identification of Pyrazole Derivatives of Usnic Acid as Novel Inhibitor of SARS-CoV-2 Main Protease Through Virtual Screening Approaches

The infection produced by the SARS-CoV-2 virus remains a significant health crisis worldwide. The lack of specific medications for COVID-19 necessitates a concerted effort to find the much-desired therapies for this condition. The main protease (Mpro) of SARS-CoV-2 is a promising target, vital for virus replication and transcription. In this study, fifty pyrazole derivatives were tested for their pharmacokinetics and drugability, resulting in eight hit compounds. Subsequent molecular docking simulations on SARS-CoV-2 main protease afforded two lead compounds with strong affinity at the active site. Additionally, the molecular dynamics (MD) simulations of lead compounds (17 and 39), along with binding free energy calculations, were accomplished to validate the stability of the docked complexes and the binding poses achieved in docking experiments. Based on these findings, compound 17 and 39, with their favorable projected pharmacokinetics and pharmacological characteristics, are the proposed potential antiviral candidates which require further investigation to be used as anti-SARS-CoV-2 medication.

In silico Molecular Docking Approach to Identify Potential Antihypertensive Compounds from Ajuga integrifolia Buch.-Ham. Ex D. Don (Armagusa)

Background

Ajuga integrifolia (Armagusa) is used as a decoction to treat high blood pressure and diabetes, widely in Ethiopia. Specific compounds for anti-hypertension activity were not identified so far. This study aims to provide a scientific basis for the therapeutic use of A. integrifolia as an antihypertension agent.

Methods

In silico studies were used to evaluate the antihypertensive components of A. integrifolia. Flavonoids identified using HPLC analysis and iridoid glycosides isolated from A. integrifolia in this study and those isolated from synonyms (A. remota and A. bractosa) were considered in the molecular docking study. Interactions were studied by using Autodock vina (1.2) on PyRx 0.8 and visualizing in 2D and 3D using ligPlot+ and Discovery studio software. Activities like vasoprotection and druglikeness properties were predicted using online servers.

Results

Flavonoids such as quercetin, myricetin, and rutin were identified and quantified by HPLC analysis from different extracts of A. integrifolia. Reptoside and 8-O-acetylharpgide isolated from the aerial part of A. integrifolia. The binding energies of all 17 candidates considered in this study range from -10.2 kcal/mol to -7.5 kcal/mol and are lower than enalapril (reference drug: -5.9 kcal/mol). The binding energies, in most case, constitute hydrogen bonding. Biological activity predicted using PASS test also showed that the flavonoids have more probability of activity than the iridoid glycosides. Druglikeness properties of the candidate molecules showed that most follow the Lipinski rule of five with few violations.

Conclusion

Lower binding energies involving hydrogen bonding and predicted activities concerning hypertension confirm the traditional use of the aerial part of the medicinal plant concerned. Flavonoids: rutin, myricetin, quercetin, and kaempferol take the leading role in the antihypertensive activity of the aerial part of A. integrifolia. The iridoid glycosides studied are almost similar in their effect on their antihypertensive activity and still better than the reference drug.

Unveiling the ESR1 Conformational Stability and Screening Potent Inhibitors for Breast Cancer Treatment

BACKGROUND

The current study recognizes the significance of estrogen receptor alpha (ERα) as a member of the nuclear receptor protein family, which holds a central role in the pathophysiology of breast cancer. ERα serves as a valuable prognostic marker, with its established relevance in predicting disease outcomes and treatment responses.

METHODS

In this study, computational methods are utilized to search for suitable drug-like compounds that demonstrate analogous ligand binding kinetics to ERα.

RESULTS

Docking-based simulation screened out the top 5 compounds - ZINC13377936, NCI35753, ZINC35465238, ZINC14726791, and NCI663569 against the targeted protein. Further, their dynamics studies reveal that the compounds ZINC13377936 and NCI35753 exhibit the highest binding stability and affinity.

CONCLUSION

Anticipating the competitive inhibition of ERα protein expression in breast cancer, we envision that both ZINC13377936 and NCI35753 compounds hold substantial promise as potential therapeutic agents. These candidates warrant thorough consideration for rigorous In vitro and In vivo evaluations within the context of clinical trials. The findings from this current investigation carry significant implications for the advancement of future diagnostic and therapeutic approaches for breast cancer.

Pleiotropic Potential of Evernia prunastri Extracts and Their Main Compounds Evernic Acid and Atranorin: In Vitro and In Silico Studies

Evernia prunastri is a lichen widely distributed in the Northern Hemisphere. Its biological properties still need to be discovered. Therefore, our paper focuses on studies of E. prunastri extracts, including its main metabolites evernic acid (EA) or atranorin (ATR). Phytochemical profiles using chromatographic analysis were confirmed. The antioxidant activity was evaluated using in vitro chemical tests and in vitro enzymatic cells-free tests, namely superoxide dismutase (SOD), glutathione peroxidase (GPx), glutathione reductase (GR), and catalase (CAT). The anti-inflammatory potential using cyclooxygenase-2 (COX-2) and hyaluronidase were determined. The neuroprotective potential using acetylcholinesterase, (AChE), butyrylcholinesterase (BChE), and tyrosinase (Tyr) was estimated. The hypoglycemic activity was also confirmed (α-glucosidase). Principal component analysis was performed to determine the relationship between the biological activity of extracts. The inhibitory effect of EA and ATR on COX-2 AChE, BChE, Tyr, and α-glucosidase was evaluated using molecular docking techniques and confirmed for EA and ATR (besides α-glucosidase). The penetration of EA and ATR from extracts through the blood-brain barrier was confirmed using the parallel artificial membrane permeability assay blood-brain barrier test. In conclusion, depending on chemical surroundings and the concentration, the E. prunastri extracts, EA or ATR, showed attractive pleiotropic properties, which should be further investigated.

Biosynthetic gene cluster synteny: Orthologous polyketide synthases in Hypogymnia physodes, Hypogymnia tubulosa, and Parmelia sulcata

Lichens are symbiotic associations consisting of a photobiont (algae or cyanobacteria) and a mycobiont (fungus), which together generate a variety of unique secondary metabolites. To access this biosynthetic potential for biotechnological applications, deeper insights into the biosynthetic pathways and corresponding gene clusters are necessary. Here, we provide a comparative view of the biosynthetic gene clusters of three lichen mycobionts derived from Hypogymnia physodes, Hypogymnia tubulosa, and Parmelia sulcata. In addition, we present a high-quality PacBio metagenome of Parmelia sulcata, from which we extracted the mycobiont bin containing 214 biosynthetic gene clusters. Most biosynthetic gene clusters in these genomes were associated with T1PKSs, followed by NRPSs and terpenes. This study focused on biosynthetic gene clusters related to polyketide synthesis. Based on ketosynthase homology, we identified nine highly syntenic clusters present in all three species. Among the four clusters belonging to nonreducing PKSs, two are putatively linked to lichen substances derived from orsellinic acid (orcinol depsides and depsidones, e.g., lecanoric acid, physodic acid, lobaric acid), one to compounds derived from methylated forms of orsellinic acid (beta orcinol depsides, e.g., atranorin), and one to melanins. Five clusters with orthologs in all three species are linked to reducing PKSs. Our study contributes to sorting and dereplicating the vast PKS diversity found in lichenized fungi. High-quality sequences of biosynthetic gene clusters of these three common species provide a foundation for further exploration into biotechnological applications and the molecular evolution of lichen substances.

Biosynthetic Potential of Hypogymnia Holobionts: Insights into Secondary Metabolite Pathways

Lichens are symbiotic associations consisting of a photobiont (algae or cyanobacteria) and a mycobiont (fungus). They are known to produce a variety of unique secondary metabolites. To access this biosynthetic potential for biotechnological applications, deeper insights into the biosynthetic pathways and corresponding gene clusters are necessary. Here we provide a comprehensive view of the biosynthetic gene clusters of all organisms comprising a lichen thallus: fungi, green algae, and bacteria. We present two high-quality PacBio metagenomes, in which we identified a total of 460 biosynthetic gene clusters. Lichen mycobionts yielded 73-114 clusters, other lichen associated ascomycetes 8-40, green algae of the genus Trebouxia 14-19, and lichen-associated bacteria 101-105 clusters. The mycobionts contained mainly T1PKSs, followed by NRPSs, and terpenes; Trebouxia reads harbored mainly clusters linked to terpenes, followed by NRPSs and T3PKSs. Other lichen-associated ascomycetes and bacteria contained a mix of diverse biosynthetic gene clusters. In this study, we identified for the first time the biosynthetic gene clusters of entire lichen holobionts. The yet untapped biosynthetic potential of two species of the genus Hypogymnia is made accessible for further research.

Tridepsides as potential bioactives: a review on their chemistry and the global distribution of their lichenic and non-lichenic natural sources

Tridepsides, as fully oxidized polyketides, have been known to exist in lichens for more than a century. Recent studies have showed that these possible defensive lichenochemicals possess various biological activities. Also, a candidate biosynthetic gene cluster was recently reported for gyrophoric acid (GA), an important tridepside. The present study focused on biosynthesis, natural sources, biological activities, and bioanalytical methods of tridepside molecules. Our survey shows that, so far, lichenic tridepsides have been reported from 37 families, 111 genera, and 526 species of lichen. Because many of their species contain tridepsides, the families Parmeliaceae, Lobariaceae, and Peltigeraceae can be considered critical lichenic sources of tridepsides. Furthermore, several species of Hypotrachyna in Parmeliaceae family showed lichenic tridepsides, suggesting that this genus is a viable source of tridepsides. This research also explored tridepsides from non-lichenic sources, such as non-lichenized fungi, lichenicolous fungi, endophytes, parasites, and liverworts, which offer substantial potential as biotechnological sources to produce tridepsides, which are produced in small amounts in lichen thalli. Two lichenic tridepsides have also been detected in non-lichenic sources: GA and tenuiorin (TE). Additionally, no significant correlation was found between tridepside biosynthesis and geographical distribution patterns for several potentially tridepside-producing lichens. We further showed that GA is the most studied tridepside with various reported biological activities, including anticancer, wound healing, photoprotection, anti-aging, antioxidant, cardiovascular effect, DNA interaction, anti-diabetes, anti-Alzheimer's, anti-bacterial, and antifungal. Last but not least, this study provides an overview of some bioanalytical methods used to analyze tridepsides over the past few years.

Molecular Docking of the Terpenes in Gorgonian Corals to COX-2 and iNOS Enzymes as Anti-Inflammatory

Background:

Because the inflammatory pathway is triggered by the enzymes cyclooxygenase- 2 (COX-2) and inducible nitric oxide synthase (iNOS), inhibitors, such as nonsteroidal anti-inflammatory drugs (NSAIDs), are needed, although these have side effects. Therefore, the discovery and development of natural medicine as a lead compound are needed. The gorgonian corals have been reported to contain cyclic diterpenes with anti-inflammatory activities. The specific anti-inflammatory inhibitor potential has not been reported regarding these secondary metabolites, whether in COX-2 or iNOS. Thus, the in silico method is the right alternative.

Objective:

This study aimed to determine the potency of fifteen terpenes of the various gorgonian corals to COX-2 and iNOS enzymes as an anti-inflammatory.

Methods:

Molecular docking was performed using ChemDraw Ultra 12.0, Chem3D Pro 12.0, Biovia Discovery Studio 2016 Client®, Autodock Tools 4.2, prediction pharmacokinetics (Pre-ADMET), and oral administration (Lipinski rule of five).

Results:

Potential terpenes based on ΔG (kcal/mol) and Ki (nM) to COX-2 were gyrosanol B (-10,32; 27,15), gyrosanol A (-10,20; 33,57), echinolabdane A (-9,81; 64,76). Only nine terpenes were specific to COX-2 active sites, while for iNOS were palmonine F (-7.76; 2070), briarenol C (-7.55; 2910), and all test compounds binding to the iNOS active sites. Pre-ADMET prediction obtained that HIA was very excellent (70–100%), Caco-2 had moderate permeability (4–70 nm sec-1), and PPB had strong binding (> 90%). Eight terpenes qualified for the Lipinski rule of five.

Conclusion:

iNOS was a specific target for terpenes based on the free energy of binding (ΔG).

Repurposing of FDA approved drugs against Salmonella enteric serovar Typhi by targeting dihydrofolate reductase: an in silico study

Drug-resistant Salmonella enteric serovar Typhi (S. Typhi) poses a vital public health issue. To overcome drug resistance issues, the development of effective drugs with novel mechanism(s) of action is required. In this regard, drug repurposing is a viable alternative approach to find novel drugs to overcome drug resistance. Therefore, a FDA-approved-drug-library containing 1930 drugs was analyzed against the dihydrofolate reductase (DHFR) of S. Typhi using deep learning regression algorithms. Initially, a total of 500 compounds were screened, followed by rescreening by molecular docking. Further, from screened compounds by molecular docking, the top eight compounds were subjected to molecular dynamics (MD) simulation. Analysis of MD simulation resulted in four potential compounds, namely; Duvelisib, Amenamevir, Lifitegrast and Nilotinib against the DHFR enzyme. During the MD simulation, these four drugs achieved good stability during the 100 ns trajectory period at 300 K. Further to know the insights of the complex's stability, we calculated RMSF, RG, SASA and interaction energy for the last 60 ns trajectory period because all complexes showed the stability after 40 ns trajectory period. MM-PBSA analysis of the last 10 ns of MD trajectories showed the stability of the complexes. From our results, we conclude that these drugs can also be useful for treating typhoid fever and can inhibit S. Typhi by interfering with the function of the DHFR enzyme. Communicated by Ramaswamy H. Sarma.

In vitro and in silico screening of novel typhoid drugs from endangered herb (Enicostema axillare)

In the present study, antibacterial activity of swertiamarin from Enicostema axillare (Lam) was checked against three different human gram-negative pathogens namely Salmonella typhi, Klebsiella pneumoniae and Shigella flexneri. Minimum inhibitory concentration assay revealed low dose and efficient activity of swertiamarin on the above said pathogens. Though swertiamarin is a well-studied and characterized compound, there is no experimental proof available for its antibacterial activity. To gain more insight about the antibacterial efficiency of swertiamarin against typhoid causing S. typhi, a comparative molecular docking of S. typhi OmpF (3NSG) was performed with swertiamarin and other typhoid drugs available in the market which exposed better activity strength of swertiamarin compared with that of the other drugs. Further, molecular dynamics of S. typhi OmpF-swertiamarin shows good flexibility and stability at 100 ns. The outcome of this work will definitely provide an idea of using very low dose of swertiamarin as a potent and promising drug against typhoid fever. Communicated by Ramaswamy H. Sarma.

Mollusc-Derived Brominated Indoles for the Selective Inhibition of Cyclooxygenase: A Computational Expedition

Inflammation plays an important role in different chronic diseases. Brominated indoles derived from the Australian marine mollusk Dicathais orbita (D. orbita) are of interest for their anti-inflammatory properties. This study evaluates the binding mechanism and potentiality of several brominated indoles (tyrindoxyl sulfate, tyrindoleninone, 6-bromoisatin, and 6,6′-dibromoindirubin) against inflammatory mediators cyclooxygenases-1/2 (COX-1/2) using molecular docking, followed by molecular dynamics simulation, along with physicochemical, drug-likeness, pharmacokinetic (pk), and toxicokinetic (tk) properties. Molecular docking identified that these indole compounds are anchored, with the main amino acid residues, positioned in the binding pocket of the COX-1/2, required for selective inhibition. Moreover, the molecular dynamics simulation based on root mean square deviation (RMSD), radius of gyration (Rg), solvent accessible surface area (SASA), and root mean square fluctuation (RMSF) analyses showed that these natural brominated molecules transit rapidly to a progressive constant configuration during binding with COX-1/2 and seem to accomplish a consistent dynamic behavior by maintaining conformational stability and compactness. The results were comparable to the Food and Drug Administration (FDA)-approved selective COX inhibitor, aspirin. Furthermore, the free energy of binding for the compounds assessed by molecular mechanics–Poisson–Boltzmann surface area (MM–PBSA) confirmed the binding capacity of indoles towards COX-1/2, with suitable binding energy values except for the polar precursor tyrindoxyl sulfate (with COX-1). The physicochemical and drug-likeness analysis showed zero violations of Lipinski’s rule, and the compounds are predicted to have excellent pharmacokinetic profiles. These indoles are projected to be non-mutagenic and free from hepatotoxicity, with no inhibition of human ether-a-go–go gene (hERG) I inhibitors, and the oral acute toxicity LD₅₀ in rats is predicted to be similar or lower than aspirin. Overall, this work has identified a plausible mechanism for selective COX inhibition by natural marine indoles as potential therapeutic candidates for the mitigation of inflammation.

An Insilico evaluation of phytocompounds from Albizia amara and Phyla nodiflora as cyclooxygenase-2 enzyme inhibitors

PURPOSE

The enzyme Cyclooxygenases (COX-1 and COX-2) catalyze the formation of prostaglandin, a mediator of the inflammatory pathway. Inflammation related pathological conditions may be alleviated by targeting the Cox enzymes.COX-2 inhibitors that are currently available in the market causes undesirable side effects. Our present study focuses on the in-silico inhibition of COX -2 enzyme by the phytocompounds from Albizia amara and Phyla nodiflora.

METHODS

The phytochemicals present in Albizia amara and Phyla nodiflora were analyzed for their COX-2 inhibition potential. Eight compounds from Albizia amara and eleven compounds from Phyla nodiflora obtained from GC-MS analysis was used for the current study. Molecular docking was performed using AutoDock vina. The crystal structure of COX-2 (PDB ID: 5IKR) was obtained from Protein data bank. PyMol was used to remove any solvent, organic and inorganic molecules. Energy minimization of the protein was carried out using SPDBV software. Geometrical optimizations of the ligands were performed using Avogadro software. Celecoxib was used as the positive control. ADMET properties of the compounds were analyzed using SwissADME and ProtoxII online servers. Molecular mechanics/generalized born surface area (MM/GBSA) calculations were performed to evaluate the binding efficiency. Molecular dynamics of the protein and protein-ligand complex was studied for about 100 ns using Desmond package of Schrodinger suite.

RESULTS

Among the eighteen compounds, Squalene present in both the plants showed a better binding energy of -7.7 kcal/mol, when compare to other phytocompounds present in the extract. The control celecoxib showed a binding energy of about - 9.4 kcal/mol. The toxicity and ADMET properties of squalene indicated that it is non-toxic and followed Lipinski's rule. Molecular Dynamics (MD) analysis showed that the binding of squalene to the enzyme was stable.

CONCLUSION

Squalene could potentially inhibit COX2 and o wing to its properties, squalene can be formulated in gels/creams and could be possibly used for external edema and inflammation.

Virtual screening-driven drug discovery of SARS-CoV2 enzyme inhibitors targeting viral attachment, replication, post-translational modification and host immunity evasion infection mechanisms

The novel coronavirus SARS-CoV2, the causative agent of the pandemic disease COVID-19, emerged in December 2019 forcing lockdown of communities in many countries. The absence of specific drugs and vaccines, the rapid transmission of the virus, and the increasing number of deaths worldwide necessitated the discovery of new substances for anti-COVID-19 drug development. With the aid of bioinformatics and computational modelling, ninety seven antiviral secondary metabolites from fungi were docked onto five SARS-CoV2 enzymes involved in viral attachment, replication, post-translational modification, and host immunity evasion infection mechanisms followed by molecular dynamics simulation and in silico ADMET prediction (absorption, distribution, metabolism, excretion and toxicity) of the hit compounds. Thus, three fumiquinazoline alkaloids scedapin C (15), quinadoline B (19) and norquinadoline A (20), the polyketide isochaetochromin D1 (8), and the terpenoid 11a-dehydroxyisoterreulactone A (11) exhibited high binding affinities on the target proteins, papain-like protease (PLpro), chymotrypsin-like protease (3CLpro), RNA-directed RNA polymerase (RdRp), non-structural protein 15 (nsp15), and the spike binding domain to GRP78. Molecular dynamics simulation was performed to optimize the interaction and investigate the stability of the top-scoring ligands in complex with the five target proteins. All tested complexes were found to have dynamic stability. Of the five top-scoring metabolites, quinadoline B (19) was predicted to confer favorable ADMET values, high gastrointestinal absorptive probability and poor blood-brain barrier crossing capacities.Communicated by Ramaswamy H. Sarma.

Berberis lyceum and Fumaria indica: in vitro cytotoxicity, antioxidant activity, and in silico screening of their selected phytochemicals as novel hepatitis C virus nonstructural protein 5A inhibitors

Berberis lyceum and Fumaria indica are two Pakistani indigenous herbal medicines used to treat liver infections, including hepatitis C virus (HCV). This study aimed to evaluate the cytotoxicity, and antioxidant activity of these plant extracts and computationally screen their selected phytoconstituents as HCV NS5A inhibitors. The viability of HepG2 cells was assessed 24 h and 48 h post-treatment using colorimetric and dye exclusion methods. Antioxidant properties were examined by the 2,2-diphenyl-1-picrylhydrazyl (DPPH), reducing power, and total antioxidant capacity assays. Seventeen known phytochemicals identified from each plant were docked into the active binding site of HCV NS5A protein. The top hit ligands were analyzed for their druglikeness properties and the indices of absorption, distribution, metabolism, elimination, and toxicity (ADMET). The results showed that both plant extracts were non-toxic (CC₅₀ > 200 µg/ml). The IC₅₀ values of DPPH-radical scavenging activity were 51.02 ± 0.94 and 62.91 ± 1.85 µg/ml for B. lyceum and F. indica, respectively. They also exhibited reducing power and total antioxidant capacity.The phytochemicals were identified as potent HCV NS5A inhibitors with good druglikeness and ADMET properties. Six of the docked phytochemicals exhibited higher binding scores (-17.9 to -19.2 kcal/mol) with HCV NS5A protein than the standard drug, daclatasvir (-17.2 kcal/mol). Molecular dynamics (MD) simulation confirmed the stability of two compounds, berbamine and paprafumine at 100 ns with active site of HCV NS5A protein. The identified compounds through molecular docking and MD simulation could have potential as HCV NS5A inhibitor after further validation. Communicated by Ramaswamy H. Sarma.

In-silico screening of plant-derived antivirals against main protease, 3CLpro and endoribonuclease, NSP15 proteins of SARS-CoV-2

Novel Coronavirus or SARS-CoV-2 outbreak has developed a pandemic condition all over the world. The virus is highly infectious and spreads by human to human local transmission mode. Till date, there is no vaccination or drugs been approved for the treatment by the World Health Organisation. Henceforth, the discovery of the potential drugs is an urgent and utmost requirement for the medical fraternity. Since, the side effects of plant-derived compounds will be lower compared to synthetic/chemical drugs. The Main protease (3CL^pro or NSP5) and endoribonuclease (NSP15) proteins are necessity for viral replication and its survival in the host cell. In the present study, in-silico approach of drug development was used to search for potential antiviral plant-derived compounds as inhibitors against SARS-CoV-2 replication proteins. Eight plant-derived compounds of which the antiviral activity was known and available, and two reported drugs against SARS-CoV-2 selected for the molecular docking analysis. The docking results suggested that bisdemethoxycurcumin, demethoxycurcumin, scutellarin, quercetin and myricetin showed least binding energy, i.e., greater than −6.5 Kcal/mol against 3CL^pro and endoribonuclease of SARS-CoV-2. Further studies of ADME-Tox and bioavailability of drugs were also performed that exhibited efficient parameters of drug likeness. Molecular dynamics simulation calculations were performed for the most negative binding affinity of the compound to evaluate the dynamic behavior,and stability of protein-ligand complex. Our findings suggest that these compounds could be potential inhibitors of SARS‐CoV-2 main protease and endoribonuclease. However, further in-vitro and pre-clinical experiments would validate the potential inhibitors of SARS‐CoV‐2 proteins.

Bacterial protein azurin and derived peptides as potential anti-SARS-CoV-2 agents: insights from molecular docking and molecular dynamics simulations

The current pandemic SARS-CoV-2 has wreaked havoc in the world, and neither drugs nor vaccine is available for the treatment of this disease. Thus, there is an immediate need for novel therapeutics that can combat this deadly infection. In this study, we report the therapeutic assessment of azurin and its peptides: p18 and p28 against the viral structural S-protein and non-structural 3CL^pro and PL^pro proteins. Among the analyzed complexes, azurin docked relatively well with the S2 domain of S-protein compared to the other viral proteins. The derived peptide p18 bound to the active site domain of the PL^pro protein; however, in other complexes, lesser interactions were recorded. The second azurin derived peptide p28, fared the best among the docked proteins. p28 interacted with all the three viral proteins and the host ACE-2 receptor by forming several electrostatic and hydrogen bonds with the S-protein, 3CL^pro, and PL^pro. MD simulations indicated that p28 exhibited a strong affinity to S-protein and ACE-2 receptor, indicating a possibility of p28 as a protein-protein interaction inhibitor. Our data suggest that the p28 has potential as an anti-SARS-CoV-2 agent and can be further exploited to establish its validity in the treatment of current and future SARS-CoV crisis.Communicated by Ramaswamy H. Sarma.

Structure-based screening of novel lichen compounds against SARS Coronavirus main protease (Mpro) as potentials inhibitors of COVID-19

The outbreak of SARS-CoV-2 and deaths caused by it all over the world have imposed great concern on the scientific community to develop potential drugs to combat Coronavirus disease-19 (COVID-19). In this regard, lichen metabolites may offer a vast reservoir for the discovery of antiviral drug candidates. Therefore, to find novel compounds against COVID-19, we created a library of 412 lichen compounds and subjected to virtual screening against the SARS-CoV-2 Main protease (Mpro). All the ligands were virtually screened, and 27 compounds were found to have high affinity with Mpro. These compounds were assessed for drug-likeness analysis where two compounds were found to fit well for redocking studies. Molecular docking, drug-likeness, X-Score, and toxicity analysis resulting in two lichen compounds, Calycin and Rhizocarpic acid with Mpro-inhibiting activity. These compounds were finally subjected to molecular dynamics simulation to compare the dynamics behavior and stability of the Mpro after ligand binding. The binding energy was calculated by MM-PBSA method to determine the intermolecular protein-ligand interactions. Our results showed that two compounds; Calycin and Rhizocarpic acid had the binding free energy of - 42.42 kJ mol/1 and - 57.85 kJ mol/1 respectively as compared to reference X77 (- 91.78 kJ mol/1). We concluded that Calycin and Rhizocarpic acid show considerable structural and pharmacological properties and they can be used as hit compounds to develop potential antiviral agents against SARS-CoV-2. These lichen compounds may be a suitable candidate for further experimental analysis.

Quantitative structure activity relationship and molecular simulations for the exploration of natural potent VEGFR-2 inhibitors: an in silico anti-angiogenic study

VEGFR-2 has recently become an eye-catching molecular target for the novel therapeutic designs against cancer for its well known role in persuading angiogenesis in tumor cells. The current study set sights on the exploration of novel potent natural compound targeting VEGFR-2 via computational ligand-based modeling and database screening followed by binding pattern analysis, reactivity site prediction and MD simulation studies. The known 53 VEGFR-2 inhibitors (with IC50 ranging from 0.7 nM to 9700 nM) were headed for development of Ligand based pharmacophore model using 3 D QSAR pharmacophore generation module of DS Client. Training set inhibitors (23 compounds) were exploited to create pharmacophore model based on their chemical features. The model was validated through 30 test set inhibitors and exploited further for screening of 62,082 natural compounds from InterBioscreen natural compound database. Screened compounds further went through Drug-Likeliness study, ADMET prediction, Binding pattern analysis, In silico prediction of reactivity sites, Biological activity spectra prediction, pan assay interference compound identification and MD simulation analysis. Out of 5 screened compounds, Compound A and Compound B exhibited highest binding energy judged against the standard drug "Sorafenib". On further conducting reactivity site prediction, BAS prediction, and pan assay interference compound identification, Compound B exhibited better result which was carried forward for MD simulation study for 50 ns. MD simulation results suggested that Compound B exhibited more stable binding to the active site of VEGFR-2 without causing any conformational changes in protein-ligand complex. Thereby, the investigation proposes Compound B to hold potent antiangiogenic potential targeting VEGFR-2. [Formula: see text] Communicated by Ramaswamy H. Sarma.

Molecular docking and molecular dynamics simulation approach to screen natural compounds for inhibition of Xanthomonas oryzae pv. Oryzae by targeting peptide deformylase

Xanthomonas oryzae pv. Oryzae (Xoo) causes bacterial leaf blight (BLB) of rice which results in a huge loss in production. Many chemicals are used to control BLB disease. However, these chemicals are toxic to the environments, animals and human beings. Thus, there is a demand to discover potential and safe natural pesticides to manage BLB disease successfully. Therefore, we screened a library of phytochemicals of different plants having antibacterial activity by targeting Peptide Deformylase (PDF) of Xoo using in silico techniques. A library of 318 phytochemicals was prepared and subjected to rigid and flexible molecular docking against PDF followed by molecular dynamics simulation and free energy analysis of protein-ligand complexes. The results of virtual screening showed that 14 compounds from different plants have good binding energy as compare to reference molecule (3 R)-2,3-dihydro[1,3] thiazolo [3,2 a]benzimidazol-3-ol) (-7.7 kcal mol^-1). Out of 14 hit compounds, eight compounds that were selected based on binding energy were analyzed by Molecular dynamic (MD) simulation. Analysis of MD simulation revealed that eight compounds namely; Bisdemethoxycurcumin, Rosmarinic acid, Piperanine, Dihydropiperlonguminine, Piperdardine, Dihydrocurcumin and Lonhumosides B achieved good stability during the 80 ns MD simulation at 300 K in term of the RMSD. Further, we calculated RMSF, RG, SASA, and interaction energy after 40 ns due to showing the stability of complexes. From our results, we conclude that these natural compounds could inhibit Xoo by targeting PDF receptor and can be used as potential bactericidal candidates against BLB disease of rice against Xoo and other bacteria. Communicated by Ramaswamy H. Sarma.