Interaction mechanisms of a melatonergic inhibitor in the melatonin synthesis pathway

Identifying plant-derived antiviral alkaloids as dual inhibitors of SARS-CoV-2 main protease and spike glycoprotein through computational screening

COVID-19 is currently considered the ninth-deadliest pandemic, spreading through direct or indirect contact with infected individuals. It has imposed a consistent strain on both the financial and healthcare resources of many countries. To address this challenge, there is a pressing need for the development of new potential therapeutic agents for the treatment of this disease. To identify potential antiviral agents as novel dual inhibitors of SARS-CoV-2, we retrieved 404 alkaloids from 12 selected medicinal antiviral plants and virtually screened them against the renowned catalytic sites and favorable interacting residues of two essential proteins of SARS-CoV-2, namely, the main protease and spike glycoprotein. Based on docking scores, 12 metabolites with dual inhibitory potential were subjected to drug-likeness, bioactivity scores, and drug-like ability analyses. These analyses included the ligand-receptor stability and interactions at the potential active sites of target proteins, which were analyzed and confirmed through molecular dynamic simulations of the three lead metabolites. We also conducted a detailed binding free energy analysis of pivotal SARS-CoV-2 protein inhibitors using molecular mechanics techniques to reveal their interaction dynamics and stability. Overall, our results demonstrated that 12 alkaloids, namely, adouetine Y, evodiamide C, ergosine, hayatinine, (+)-homoaromoline, isatithioetherin C, N,alpha-L-rhamnopyranosyl vincosamide, pelosine, reserpine, toddalidimerine, toddayanis, and zanthocadinanine, are shortlisted as metabolites based on their interactions with target proteins. All 12 lead metabolites exhibited a higher unbound fraction and therefore greater distribution compared with the standards. Particularly, adouetine Y demonstrated high docking scores but exhibited a nonspontaneous binding profile. In contrast, ergosine and evodiamide C showed favorable binding interactions and superior stability in molecular dynamics simulations. Ergosine demonstrated exceptional performance in several key pharmaceutical metrics. Pharmacokinetic evaluations revealed that ergosine exhibited pronounced bioactivity, good absorption, and optimal bioavailability. Additionally, it was predicted not to cause skin sensitivity and was found to be non-hepatotoxic. Importantly, ergosine and evodiamide C emerged as superior drug candidates for dual inhibition of SARS-CoV-2 due to their strong binding affinity and drug-like ability, comparable to known inhibitors like N3 and molnupiravir. This study is limited by its in silico nature and demands the need for future in vitro and in vivo studies to confirm these findings.

Efficacy assessment of novel methanimine derivatives in chronic constriction injury-induced neuropathic model: An in-vivo, ex-vivo and In-Silico approach

The multicomponent etiology, complex clinical implications, dose-based side effect and degree of pain mitigation associated with the current pharmacological therapy is incapable in complete resolution of chronic neuropathic pain patients which necessitates the perpetual requirement of novel medication therapy. Therefore, this study explored the ameliorative aptitude of two novel methanimine imitative like (E)-N-(4-nitrobenzylidene)-4‑chloro-2-iodobenzamine (KB 09) and (E)-N-(4-methylbenzylidene)-4‑chloro-2-iodobenzamine (KB 10) in chronic constriction injury (CCI) of sciatic nerve induced neuropathic pain in rat model. Standard behavioral tests like dynamic and static allodynia, cold, thermal and mechanical hyperalgesia along with rotarod activity were performed at various experimental days like 0, 3, 7, 14 and 21. Enzyme linked immunosorbent assay (ELISA) on spinal tissue and antioxidant assays on sciatic nerve were executed accompanied by molecular docking and simulation studies. Prolonged ligation of sciatic nerve expressively induced hyperalgesia as well as allodynia in rats. KB 09 and KB 10 substantially attenuated the CCI elicited hyperalgesia and allodynia. They significantly reduced the biomarkers of pain and inflammation like Interleukin-6 (IL-6) and tumor necrosis factor-alpha (TNF-α) in ELISA and while enhanced the GSH, SOD and CAT and diminished the MDA levels during antioxidant assays. KB 09 displayed -9.62 kcal/mol with TNF-α and -7.68 kcal/mol binding energy with IL-6 whereas KB 10 exhibited binding energy of -8.20 kcal/mol with IL-6 while -11.68 kcal/mol with TNF-α and hence both trial compounds ensured stable interaction with IL-6 and TNF-α during computational analysis. The results advocated that both methanimine derivatives might be novel candidates for attenuation of CCI-induced neuropathic pain prospects via anti-nociceptive, anti-inflammatory and antioxidant mechanisms.

A hierarchical approach towards identification of novel inhibitors against L, D-transpeptidase YcbB as an anti-bacterial therapeutic target

The bacterial cell wall, being a vital component for cell viability, is regarded as a promising drug target. The L, D-Transpeptidase YcbB enzyme has been implicated for a significant role in cell wall polymers cross linking during typhoid toxin release, β-lactam resistance and outer membrane defect rescue. These observations have been recorded in different bacterial pathogens such as Salmonella Typhimurium, Citrobacter rodentium, and Salmonella typhi. In this work, we have shown structure based virtual screening of diverse natural and synthetic drug libraries against the enzyme and revealed three compounds as LAS_32135590, LAS_34036730 and LAS-51380924. These compounds showed highly stable energies and the findings are very competitive with the control molecule ((1RG or (4 R,5S)-3-({(3S,5S)-5-[(3-carboxyphenyl)carbamoyl]pyrrolidin-3-yl}sulfanyl)-5-[(1S,2R)-1-formyl-2-hydroxypropyl]-4-methyl-4,5-dihydro-1H-pyrrole-2-carboxylic acid or ertapenem)) used. Compared to control (which has binding energy score of -11.63 kcal/mol), the compounds showed better binding energy. The binding energy score of LAS_32135590, LAS_34036730 and LAS-51380924 is -12.63 kcal/mol, -12.22 kcal/mol and -12.10 kcal/mol, respectively. Further, the docked snapshot of the lead compounds and control were investigated for stability under time dependent dynamics environment. All the three leads complex and control system showed significant equilibrium (mean RMSD < 3 Å) both in term of intermolecular docked conformation and binding interactions network. Further validation on the complex's stability was acquired from the end-state MMPB/GBSA analysis that observed greater contribution from van der Waals forces and electrostatic energy while less contribution was noticed from solvation part. The compounds were also showed good drug-likeness and are non-toxic and non-mutagenic. In short, the compounds can be used in experimental testing's and might be subjected to structure modification to get better results.Communicated by Ramaswamy H. Sarma.

Structural Dynamics of P-Rex1 Complexed with Natural Leads Establishes the Protein as an Attractive Target for Therapeutics to Suppress Cancer Metastasis

Phosphatidylinositol 3,4,5-trisphosphate- (PIP3-) dependent Rac exchanger 1 (P-Rex1) functions as Rho guanine nucleotide exchange factor and is activated by synergistic activity of Gβγ and PIP3 of the heterotrimeric G protein. P-Rex1 activates Rac GTPases for regulating cell invasion and migration and promotes metastasis in several human cancers including breast, prostate, and skin cancer. The protein is a promising therapeutic target because of its multifunction roles in human cancers. Herein, the present study attempts to identify selective P-Rex1 natural inhibitors by targeting PIP3-binding pocket using large-size multiple natural molecule libraries. Each library was filtered subsequently in FAF-Drugs4 based on Lipinski's rule of five (RO5), toxicity, and filter pan assay interference compounds (PAINS). The output hits were virtually screened at the PIP3-binding pocket through PyRx AutoDock Vina and cross-checked by GOLD. The best binders at the PIP3-binding pocket were prioritized using a comparative analysis of the docking scores. Top-ranked two compounds with high GOLD fitness score (>80) and lowest AutoDock binding energy (< -12.7 kcal/mol) were complexed and deciphered for molecular dynamics along with control-P-Rex1 complex to validate compound binding conformation and disclosed binding interaction pattern. Both the systems were seen in good equilibrium, and along the simulation time, the compounds are in strong contact with the P-Rex1 PIP3-binding site. Hydrogen bonding analysis towards simulation end identified the formation of 16 and 22 short- and long-distance hydrogen bonds with different percent of occupancy to the PIP3 residues for compound I and compound 2, respectively. Radial distribution function (RDF) analysis of the key hydrogen bonds between the compound and the PIP3 residues demonstrated a strong affinity of the compounds to the mentioned PIP3 pocket. Additionally, MMGB/PBSA energies were performed that confirmed the dominance of Van der Waals energy in complex formation along with favorable contribution from hydrogen bonding. These findings were also cross-validated by a more robust WaterSwap binding energy predictor, and the results are in good agreement with a strong binding affinity of the compounds for the protein. Lastly, the key contribution of residues in interaction with the compounds was understood by binding free energy decomposition and alanine scanning methods. In short, the results of this study suggest that P-Rex1 is a good druggable target to suppress cancer metastasis; therefore, the screened druglike molecules of this study need in vitro and in vivo anti-P-Rex1 validation and may serve as potent leads to fight cancer.

Biological Evaluation and Computational Studies of Methoxy-flavones from Newly Isolated Radioresistant Micromonospora aurantiaca Strain TMC-15

This study aims to determine UV-B resistance and to investigate computational analysis and antioxidant potential of methoxy-flavones of Micromonospora aurantiaca TMC-15 isolated from Thal Desert, Pakistan. The cellular extract was purified through solid-phase extraction and UV-Vis spectrum analysis indicated absorption peaks at λmax 250 nm, 343 nm, and 380 nm that revealed the presence of methoxy-flavones named eupatilin and 5-hydroxyauranetin. The flavones were evaluated for their antioxidant as well as protein and lipid peroxidation inhibition potential using di(phenyl)-(2,4,6-trinitrophenyl) iminoazanium (DPPH), 2,4-dinitrophenyl hydrazine (DNPH), and thiobarbituric acid reactive substances (TBARS) assays, respectively. The methoxy-flavones were further studied for their docking affinity and interaction dynamics to determine their structural and energetic properties at the atomic level. The antioxidant potential, protein, and lipid oxidation inhibition and DNA damage preventive abilities were correlated as predicted by computational analysis. The eupatilin and 5-hydroxyauranetin binding potential to their targeted proteins 1N8Q and 1OG5 is - 4.1 and - 7.5 kcal/mol, respectively. Moreover, the eupatiline and 5-hydroxyauranetin complexes illustrate van der Waals contacts and strong hydrogen bonds to their respective enzymes target. Both in vitro studies and computational analysis results revealed that methoxy-flavones of Micromonospora aurantiaca TMC-15 can be used against radiation-mediated oxidative damages due to its kosmotrophic nature. The demonstration of good antioxidant activities not only protect DNA but also protein and lipid oxidation and therefore could be a good candidate in radioprotective drugs and as sunscreen due to its kosmotropic nature.

A cheminformatics-biophysics correlate to identify promising lead molecules against matrix metalloproteinase-2 (MMP-2) enzyme: A promising anti-cancer target

Matrix metalloproteinase-2 (MMP-2) is an endopeptidase enzyme that is devoted to extracellular matrix proteins degradation. The enzyme is warranted as promising drugs target for different light threating diseases such as arthritis, cancer and fibrosis. Herein, in this study, three drug molecules: CMNPD8322, CMNPD8320, and CMNPD8318 were filtered as high affinity binding compounds with binding energy score of -9.75 kcal/mol, -9.11 kcal/mol, -9.05 kcal/mol, respectively. The control binding energy score was -9.01 kcal/mol. The compounds docked deeply inside the pocket interacting with S1 pocket residues. The docked complexes dynamics in real time at cellular environment was then done to decipher the stable binding conformation and intermolecular interactions network. The compounds complexes achieved very stable dynamics with root mean square deviation (RMSD) with mean value of around 2-3 Å compared to control complex that showed higher fluctuations of 5 Å. The simulation trajectories frames based binding free energy demonstrated all the compounds-MMP-2 complexes reported highly stable energy, particularly the van der Waals energy dominate the overall net energy. Similarly, the complexes revalidation of WaterSwap based energies also disclosed the complexes highly stable in term docked conformation. Also, the compounds illustrated the compounds favorable pharmacokinetics and were non-toxic and non-mutagenic. Thus, the compounds might be used thorough experimental assays to confirm compounds selective biological potency against MMP-2 enzyme.

Role of hinge motion and ATP dynamics in factors for inversion stimulation FIS protein deduced while targeting drug resistant Orientia tsutsugamushi

Orientia tsutsugamushi, the causative agent of scrub typhus has been found resistant to various classes of antibiotics such as penicillins, gentamycin and cephalosporins. Review of current literature suggests that the prevalence of scrub typhus has increased globally. Therefore, the current study has aimed at exploring the genome of O. tsutsugamushi to identify potential drug target proteins that can be used for developing novel antibiotics against the pathogen. Subtractive proteomics approach has revealed FIS as a potential drug target protein involved in two component system (TCS), a signaling pathway crucial for bacteria to survive and adjust in changing environmental conditions. Molecular docking studies have revealed compound-356 (CHEMBRIDGE-10040641-3710.356) as a potential inhibitor in both chains A and B of the FIS protein. Simulation results suggest that the docked complex has remained stable and compact throughout the 200 ns run. Significant conformational changes including the hinge motion was observed in the DNA binding domain. Furthermore, the presence of salt bridge between GLU910 and ARG417, rearrangement of interaction residues and displacement of ATP in the central AAA + domain upon binding to the inhibitor were also observed playing a role in stabilizing the protein structure.

Repurposing FIASMAs against Acid Sphingomyelinase for COVID-19: A Computational Molecular Docking and Dynamic Simulation Approach

Over the past few years, COVID-19 has caused widespread suffering worldwide. There is great research potential in this domain and it is also necessary. The main objective of this study was to identify potential inhibitors against acid sphingomyelinase (ASM) in order to prevent coronavirus infection. Experimental studies revealed that SARS-CoV-2 causes activation of the acid sphingomyelinase/ceramide pathway, which in turn facilitates the viral entry into the cells. The objective was to inhibit acid sphingomyelinase activity in order to prevent the cells from SARS-CoV-2 infection. Previous studies have reported functional inhibitors against ASM (FIASMAs). These inhibitors can be exploited to block the entry of SARS-CoV-2 into the cells. To achieve our objective, a drug library containing 257 functional inhibitors of ASM was constructed. Computational molecular docking was applied to dock the library against the target protein (PDB: 5I81). The potential binding site of the target protein was identified through structural alignment with the known binding pocket of a protein with a similar function. AutoDock Vina was used to carry out the docking steps. The docking results were analyzed and the inhibitors were screened based on their binding affinity scores and ADME properties. Among the 257 functional inhibitors, Dutasteride, Cepharanthine, and Zafirlukast presented the lowest binding affinity scores of −9.7, −9.6, and −9.5 kcal/mol, respectively. Furthermore, computational ADME analysis of these results revealed Cepharanthine and Zafirlukast to have non-toxic properties. To further validate these findings, the top two inhibitors in complex with the target protein were subjected to molecular dynamic simulations at 100 ns. The molecular interactions and stability of these compounds revealed that these inhibitors could be a promising tool for inhibiting SARS-CoV-2 infection.

Effect of Rumex dentatus on Gastrointestinal Protection and Toxicology in Rodents via Investigating H+/K+-ATPase, Calcium Channels, and PDE Mediated Signaling

This present study aims to delineate Rumex dentatus crude extract (Rd.Cr), n-Hexane, ethyl acetate, aqueous fractions (Rd.n-Hex, Rd.ETAC, and Rd.Aq), and emodin for antidiarrheal, antisecretory effects, anti-spasmodic, gastrointestinal transient time, anti-H. pylori, antiulcer effects, and toxicology. Plant extracts attributed dose-dependent protection against castor oil-induced diarrhea and dose-dependently inhibited intestinal fluid secretions in mice. They decreased the distance transverse by charcoal in the gastrointestinal transit model in rats. In rabbit jejunum preparations, it causes a concentration-dependent relaxation of both spontaneous and K⁺ (80 mM)-induced contraction, Rd.n-Hex and verapamil were relatively potent against K⁺-induced contractions and shifted the Ca²⁺ concentration-response curves (CRCs) to the right, Rd.Cr and Rd.ETAC shifted the isoprenaline-induced inhibitory CRCs to the left, showing potentiating effect similar to papaverine. Rd.n-Hex showed anti-H. pylori effect. Extracts and emodin also show an inhibitory effect against H⁺/K⁺-ATPase. Rumex dentatus showed a gastroprotective and antioxidant effect. Histopathological evaluation showed improvement in cellular architecture and decrease in the expression of inflammatory markers such as cyclooxygenase (COX2), tumor necrosis factor (TNF-α), and phosphorylated nuclear factor kappa B (p-NFƙB), validated through immunohistochemistry, ELISA, and western blot techniques. In RT-PCR, it decreases H⁺/K⁺-ATPase mRNA levels. Rumex dentatus was analyzed for certain safety aspects and exhibited a relative safety profile as no impairment was observed in kidneys, heart, liver, and brain further assisted by biochemical and hematological analysis. Docking studies revealed that emodin against H⁺/K⁺-ATPase pump and voltage gated L-type calcium channel showed E-value of −7.9 and −7.4 kcal/mol, respectively. MD simulations and molecular mechanics Poisson Boltzmann surface area and molecular mechanics Generalized Born surface area MMPBSA/GBSA findings are consistent with the in-vitro, in-vivo, and docking results. In conclusion, Rumex dentatus extracts and its phytoconstituent could be considered a potent antioxidant and anti-inflammatory drug candidates that possess anti-diarrheal, anti-secretary, antispasmodic, anti-H. pylori, and anti-ulcer potential. Toxicity studies were done according to OECD standards 425. It belongs to group 5 (LD50 > 2000 mg/kg), which suggests that it is in the lower toxicity class.

Synthesis, virtual screening and computational approach of a quinoxaline derivative as potent anti-HIV agent targeting the reverse transcriptase enzyme

Infection by the human immunodeficiency virus still represents a continuous serious concern and a global threat to human health. Due to the appearance of multi-resistant virus strains and the serious adverse side effects of the antiretroviral therapy administered, there is an urgent need for the development of new treatment agents that are more active, less toxic, and with increased tolerability to mutations. Quinoxaline derivatives are a class of heterocyclic compounds with a wide range of organic and remedial applications. In addition, they are known to significantly inhibit HIV reverse transcriptase (RT) and HIV replication in cell cultures. For these reasons, we are investigating the synthesis and computational studies of quinoxaline derivatives with a focus on their effects on the HIV RT enzyme, and we present here the structure of one such molecule, methyl 2-[(2E)-3-oxo-1,2,3,4-tetrahydroquinoalin-2-ylidene] acetate, which was confirmed by X-ray diffraction studies. In the crystal, N-H···O and C-H···O hydrogen bonds form ribbons whose mean planes are inclined to (111) by 25.69(8)°. The ribbons are formed into stacks by C-H···π(ring) interactions and π-stacking interactions between carbonyl groups. The Hirshfeld surface map allows us to understand the nature of interactions in the contribution to crystal packing. A density functional theory (DFT) calculation was performed to optimize the geometrical parameters and then they were compared with the solid-state phase. The molecular electrostatic potential map displays reactive sites on the surface, which are responsible for intermolecular interaction in the chemical species. Computational molecular docking, in addition to molecular dynamics simulations and MMGB/PBSA binding energy techniques, was used to assess the affinity of the molecule for the HIV reverse transcriptase enzyme. The new quinoxaline derivative is more powerful in terms of binding affinity and binding conformation stability with the HIV reverse transcriptase enzyme, which suggests the molecule is a good candidate for further biological optimization.Communicated by Ramaswamy H. Sarma.

Virtual Screening and Hit Selection of Natural Compounds as Acetylcholinesterase Inhibitors

Acetylcholinesterase (AChE) is one of the classical targets in the treatment of Alzheimer's disease (AD). Inhibition of AChE slows down the hydrolysis of acetycholine and increases choline levels, improving the cognitive function. The achieved success of plant-based natural drugs acting as AChE inhibitors, such as galantamine (GAL) from Galanthus genus and huperzine A from Huperzia serrate (approved drug in China), in the treatment of AD, and the fact that natural compounds (NCs) are considered as safer and less toxic compared to synthetic drugs, led us to screen the available NCs (almost 150,000) in the ZINC12 database for AChE inhibitory activity. The compounds were screened virtually by molecular docking, filtered for suitable ADME properties, and 32 ligands from 23 structural groups were selected. The stability of the complexes was estimated via 1 μs molecular dynamics simulation. Ten compounds formed stable complexes with the enzyme and had a vendor and a reasonable price per mg. They were tested for AChE inhibitory and antioxidant activity. Five compounds showed weak AChE inhibition and three of them exhibited high antioxidant activity.

Integrated computer-aided drug design and biophysical simulation approaches to determine natural anti-bacterial compounds for Acinetobacter baumannii

Acinetobacter baumannii is a nosocomial bacterial pathogen and is responsible for a wide range of diseases including pneumonia, necrotizing fasciitis, meningitis, and sepsis. The enzyme 5-enolpyruvylshikimate-3-phosphate (EPSP) synthase (encoded by aroA gene) in ESKAPE pathogens catalyzes the sixth step of shikimate pathway. The shikimate pathway is an attractive drug targets pathway as it is present in bacteria but absent in humans. As EPSP is essential for the A. baumannii growth and needed during the infection process, therefore it was used as a drug target herein for high-throughput screening of a comprehensive marine natural products database (CMNPD). The objective was to identify natural molecules that fit best at the substrate binding pocket of the enzyme and interact with functionally critical residues. Comparative assessment of the docking scores allowed selection of three compounds namely CMNPD31561, CMNPD28986, and CMNPD28985 as best binding molecules. The molecules established a balanced network of hydrophobic and hydrophilic interactions, and the binding pose remained in equilibrium throughout the length of molecular simulation time. Radial distribution function (RDF) analysis projected key residues from enzyme active pocket which actively engaged the inhibitors. Further validation is performed through binding free energies estimation that affirms very low delta energy of <-22 kcal/mol in MM-GBSA method and <-12 kcal/mol in MM-PBSA method. Lastly, the most important active site residues were mutated and their ligand binding potential was re-investigated. The molecules also possess good druglike properties and better pharmacokinetics. Together, these findings suggest the potential biological potency of the leads and thus can be used by experimentalists in vivo and in vitro studies.

Discovery of Potential Antiviral Compounds against Hendra Virus by Targeting Its Receptor-Binding Protein (G) Using Computational Approaches

Hendra virus (HeV) belongs to the paramyxoviridae family of viruses which is associated with the respiratory distress, neurological illness, and potential fatality of the affected individuals. So far, no competitive approved therapeutic substance is available for HeV. For that reason, the current research work was conducted to propose some novel compounds, by adopting a Computer Aided Drug Discovery approach, which could be used to combat HeV. The G attachment Glycoprotein (Ggp) of HeV was selected to achieve the primary objective of this study, as this protein makes the entry of HeV possible in the host cells. Briefly, a library of 6000 antiviral compounds was screened for potential drug-like properties, followed by the molecular docking of short-listed compounds with the Protein Data Bank (PDB) structure of Ggp. Docked complexes of top two hits, having maximum binding affinities with the active sites of Ggp, were further considered for molecular dynamic simulations of 200 ns to elucidate the results of molecular docking analysis. MD simulations and Molecular Mechanics Energies combined with the Generalized Born and Surface Area (MMGBSA) or Poisson-Boltzmann and Surface Area (MMPBSA) revealed that both docked complexes are stable in nature. Furthermore, the same methodology was used between lead compounds and HeV Ggp in complex with its functional receptor in human, Ephrin-B2. Surprisingly, no major differences were found in the results, which demonstrates that our identified compounds can also perform their action even when the Ggp is attached to the Ephrin-B2 ligand. Therefore, in light of all of these results, we strongly suggest that compounds (S)-5-(benzylcarbamoyl)-1-(2-(4-methyl-2-phenylpiperazin-1-yl)-2-oxoethyl)-6-oxo-3,6-dihydropyridin-1-ium-3-ide and 5-(cyclohexylcarbamoyl)-1-(2-((2-(3-fluorophenyl)-2-methylpropyl)amino)-2-oxoethyl)-6-oxo-3,6-dihydropyridin-1-ium-3-ide could be considered as potential therapeutic agents against HeV; however, further in vitro and in vivo experiments are required to validate this study.

Determination of Novel Anti-Cancer Agents by Targeting OGG1 Enzyme Using Integrated Bioinformatics Methods

The 8-oxoguanine DNA glycosylase (OGG1) enzyme is a key DNA glycosylase mediating the excision of 7,8-dihydro-8-oxoguanine (8-oxoG) from DNA molecule to the start base excision repair pathway. The OGG1 glycosylase function depletion has been seen to obstruct pathological conditions such as inflammation, A3 T-cell lymphoblastic acute leukemia growth, and neurodegenerative diseases, thus warranting OGG1 as an attractive anti-cancer enzyme. Herein, we employed several drug libraries intending to screen non-toxic inhibitory molecules against the active pocket of the enzyme that achieved stable binding mode in dynamics. Two anti-cancer compounds ([O-]C1=C(CC2=CC=CC=C2)SC(=[N+]1CC(=O)NC3=NC=C(CC4=CC=CC=C4)S3)S and CCCN(CCC)[S]-(=O)(=O)C1=CC=C(C=C1)C(=O)NNC2=NC3=CC=C(Br)C=C3C(=N2)C4=CC=CC=C4) from Selleckchem.com were identified to occupy the active pocket of OGG1 and bind with greater affinity than Control TH5487. The binding affinity of Top-1 is -11.6 kcal/mol while that of Top-2 is -10.7 kcal/mol in contrast to TH5487 Control (-9 kcal/mol). During molecular dynamic simulations versus time, the said compounds are tightly held by the enzyme with no minor structural deviations reported except flexible loops in particular those present at the N and C-terminal. Both the compounds produced extensive hydrophobic interactions with the enzyme along with stable hydrogen bonding. The docking and molecular dynamics simulations predictions were further validated by molecular mechanics with generalized Born and surface area solvation (MM/GBSA) and Poisson Boltzmann surface area (MM/PBSA), and WaterSwap binding energies that validated strong binding of the compounds to the enzyme. The MM/GBSA binding free energy for Top-1 complex is -28.10 kcal/mol, Top-2 complex is -50.14 kcal/mol) and Control is -46.91 kcal/mol while MM/PBSA value for Top-1, Top-2 and Control is -23.38 kcal/mol, -35.29 kcal/mol and -38.20 kcal/mol, respectively. Computational pharmacokinetics support good druglike candidacy of the compounds with acceptable profile of pharmacokinetics and very little toxicity. All these findings support the notion that the compounds can be used in experiments to test their anti-cancer activities.

Molecular Insights into Binding Mode and Interactions of Structure-Based Virtually Screened Inhibitors for Pseudomonas aeruginosa Multiple Virulence Factor Regulator (MvfR)

Multi-drug resistance (MDR) bacterial pathogens pose a threat to global health and warrant the discovery of new therapeutic molecules, particularly those that can neutralize their virulence and stop the evolution of new resistant mechanisms. The superbug nosocomial pathogen, Pseudomonas aeruginosa, uses a multiple virulence factor regulator (MvfR) to regulate the expression of multiple virulence proteins during acute and persistent infections. The present study targeted MvfR with the intention of designing novel anti-virulent compounds, which will function in two ways: first, they will block the virulence and pathogenesis P. aeruginosa by disrupting the quorum-sensing network of the bacteria, and second, they will stop the evolution of new resistant mechanisms. A structure-based virtual screening (SBVS) method was used to screen druglike compounds from the Asinex antibacterial library (~5968 molecules) and the comprehensive marine natural products database (CMNPD) (~32 thousand compounds), against the ligand-binding domain (LBD) of MvfR, to identify molecules that show high binding potential for the relevant pocket. In this way, two compounds were identified: Top-1 (4-((carbamoyloxy)methyl)-10,10-dihydroxy-2,6-diiminiodecahydropyrrolo[1,2-c]purin-9-yl sulfate) and Top-2 (10,10-dihydroxy-2,6-diiminio-4-(((sulfonatocarbamoyl)oxy)methyl)decahydropyrrolo[1,2-c]purin-9-yl sulfate), in contrast to the co-crystallized M64 control. Both of the screened leads were found to show deep pocket binding and interactions with several key residues through a network of hydrophobic and hydrophilic interactions. The docking results were validated by a long run of 200 ns of molecular dynamics simulation and MM-PB/GBSA binding free energies. All of these analyses confirmed the presence of strong complex formation and rigorous intermolecular interactions. An additional analysis of normal mode entropy and a WaterSwap assay were also performed to complement the aforementioned studies. Lastly, the compounds were found to show an acceptable range of pharmacokinetic properties, making both compounds potential candidates for further experimental studies to decipher their real biological potency.

Mechanistic evaluation of a novel cyclohexenone derivative's functionality against nociception and inflammation: An in-vitro, in-vivo and in-silico approach

The synthesis of a novel cyclohexanone derivative (CHD; Ethyl 6-(4-metohxyphenyl)-2-oxo-4-phenylcyclohexe-3-enecarboxylate) was described and the subsequent aim was to perform an in vitro, in vivo and in silico pharmacological evaluation as a putative anti-nociceptive and anti-inflammatory agent in mice. Initial in vitro studies revealed that CHD inhibited both cyclooxygenase-2 (COX-2) and 5-lipoxygenase (5-LOX) enzymes and it also reduced mRNA expression of COX-2 and the pro-inflammatory cytokines TNF-α and IL-1β. It was then shown that CHD dose dependently inhibited chemically induced tonic nociception in the abdominal constriction assay and also phasic thermal nociception (i.e. anti-nociception) in the hot plate and tail immersion tests in comparison with aspirin and tramadol respectively. The thermal test outcomes indicated a possible moderate centrally mediated anti-nociception which, in the case of the hot plate test, was pentylenetetrazole (PTZ) and naloxone reversible, implicating GABAergic and opioidergic mechanisms. CHD was also effective against both the neurogenic and inflammatory mediator phases induced in the formalin test and it also disclosed anti-inflammatory activity against the phlogistic agents, carrageenan, serotonin, histamine and xylene compared with standard drugs in edema volume tests. In silico studies indicated that CHD possessed preferential affinity for GABAA, opioid and COX-2 target sites and this was supported by molecular dynamic simulations where computation of free energy of binding also favored the formation of stable complexes with these sites. These findings suggest that CHD has prospective anti-nociceptive and anti-inflammatory properties, probably mediated through GABAergic and opioidergic interactions supplemented by COX-2 and 5-LOX enzyme inhibition in addition to reducing pro-inflammatory cytokine expression. CHD may therefore possess potentially beneficial therapeutic effectiveness in the management of inflammation and pain.

Proteomic profiling the molecular signatures of plectranthoic acid in prostate cancer cells

Among cancers, prostate cancer (PCa) is frequently detected solid tumor and a growing problem for the male population, globally. Newer treatment modalities with specific targets are required for management. Plant-derived agents/drugs have historically been useful in cancer therapeutics. Natural metabolite i.e. plectranthoic acid (PA), inhibits the proliferation of PCa cells and has potent anti-cancer potential. Herein, we aim to identify the molecular signatures of PA. Proteins from control and PA-treated PCa cells were analysed using high-throughput labeled free proteomics approach. Data was processed with the SIEVE software and thoroughly analysed by using Ingenuity pathway analysis (IPA) and PANTHER. A total of 98 unique peptides, showing >2 fold change, were identified. Results indicated that PA modulates oncogenic pro-survival and pro-apoptotic signaling pathways in PCa cells. mTOR was the major canonical pathway targeted by PA, the inhibition of which was likely to induce PA mediated apoptosis. Moreover, PA interacts with the rapamycin binding domain of mTOR, demonstrated by the molecular dynamic (MD) simulation and binding free energy calculations. Furthermore, the biological process moderated by PA with a high percentage was a metabolic process. Taken together, PA appears to have pleiotropic effects, as it modulates multiple key signaling pathways, supporting the potential usefulness. SIGNIFICANCE: Studies on the mechanism of action of therapeutic agents are crucial for drug development. These studies support the selection of a therapeutic agent, appropriate models of its efficacy, and designing of further experiments. Furthermore, information on mechanism of action may suggest strategies for combination therapies. In this regard Proteomics provide the platform for comprehensive understanding of the molecular action mechanisms of newly discovered therapeutic agents. Current research highlights the new insights into mode of action of novel therapeutic metabolite i.e. Plectranthoic acid (PA). Using labeled free proteomics approach we extracted the underlying mechanisms for the anticancer activity of PA using prostate cancer model. The result of the study will pay the way for further investigations on this potent natural compound in different cancers and will provide a root for its development as a lead.

Bacillus species; a potential source of anti-SARS-CoV-2 main protease inhibitors

The COVID-19 being a preconized global pandemic by the World Health Organization needs persuasive immediate research for possible medications. The present study was carried out with a specific aim to computationally evaluate and identify compounds derived from Bacillus species as the plausible inhibitors against 3-chymotrypsin-like main protease (3CLpro) or main protease (M^Pro), which is a key enzyme in the life-cycle of coronavirus. The compounds were isolated from the crude extracts of Bacillus species. Among the isolated compounds, novel inhibitory leads were identified using in silico techniques. Molecular docking revealed that stigmasterol (-8.3 kcal/mol), chondrillasterol (-7.9 kcal/mol) and hexadecnoic acid (-6.9 kcal/mol)) among others bind in the substrate-binding pocket and also interacted with the catalytic dyad of the 3-CLpro. Further evaluation using 50 ns molecular dynamic simulation and MMPB-GBSA indicated that among the top three docking hits, hexadecanoic acid was found to be the most promising anti-COVID-19 lead against the main protease. Hexadecanoic acid might serve as a potent anti-SARS-CoV-2 compound to combat COVID-19, however, in vitro and in vivo validation and optimization is needed.

Communicated by Ramaswamy H. Sarma

Identification of potential inhibitors of Zika virus NS5 RNA-dependent RNA polymerase through virtual screening and molecular dynamic simulations

Zika virus (ZIKV) is one of the mosquito borne flavivirus with several outbreaks in past few years in tropical and subtropical regions. The non-structural proteins of flaviviruses are suitable active targets for inhibitory drugs due to their role in pathogenicity. In ZIKV, the non-structural protein 5 (NS5) RNA-Dependent RNA polymerase replicates its genome. Here we have performed virtual screening to identify suitable ligands that can potentially halt the ZIKV NS5 RNA dependent RNA polymerase (RdRp). During this process, we searched and screened a library of ligands against ZIKV NS5 RdRp. The selected ligands with significant binding energy and ligand-receptor interactions were further processed. Among the selected docked conformations, top five was further optimized at atomic level using molecular dynamic simulations followed by binding free energy calculations. The interactions of ligands with the target structure of ZIKV RdRp revealed that they form strong bonds within the active sites of the receptor molecule. The efficacy of these drugs against ZIKV can be further analyzed through in-vitro and in-vivo studies.

Immunoinformatics characterization of SARS-CoV-2 spike glycoprotein for prioritization of epitope based multivalent peptide vaccine

The COVID-19 pandemic caused by SARS-CoV-2 is a public health emergency of international concern and thus calling for the development of effective and safe therapeutics and prophylactics particularly a vaccine to protect against the infection. SARS-CoV-2 spike glycoprotein is an attractive candidate for a vaccine, antibodies, and inhibitors development because of the many roles it plays in attachment, fusion and entry into the host cell. In the present investigation, we characterized the SARS-CoV-2 spike glycoprotein by immunoinformatics techniques to put forward potential B and T cell epitopes, followed by the use of epitopes in construction of a multi-epitope peptide vaccine construct (MEPVC). The MEPVC revealed robust host immune system simulation with high production of immunoglobulins, cytokines and interleukins. Stable conformation of the MEPVC with a representative innate immune TLR3 receptor was observed involving strong hydrophobic and hydrophilic chemical interactions, along with enhanced contribution from salt-bridges towards inter-molecular stability. Molecular dynamics simulation in aqueous milieu aided further in interpreting strong affinity of the MEPVC for TLR3. This stability is the attribute of several vital residues from both TLR3 and MEPVC as shown by radial distribution function (RDF) and a novel axial frequency distribution (AFD) analytical tool. Comprehensive binding free energies estimation was provided at the end that concluded major domination by electrostatic and minor from van der Waals. Summing all, the designed MEPVC has tremendous potential of providing protective immunity against COVID-19 and thus could be considered in experimental studies.

Deciphering the Inhibition Mechanism of under Trial Hsp90 Inhibitors and Their Analogues: A Comparative Molecular Dynamics Simulation

Heat shock protein 90 (Hsp90) performs functions in cellular activities together with other signaling pathways. Hsp90 is evolutionarily conserved and universally articulated as a human cancer-causing agent involved in lung cancer and breast cancer followed by colon and rectum cancers. It has emerged as an effective drug candidate, and inhibition may affect several signaling pathways associated with cancer spread. Therefore, in-silico approaches, molecular docking, molecular dynamics simulation, and binding free energy calculations were applied to create insights into the inhibition mechanism against Hsp90 to identify new cancer therapeutic drugs. Top-docked Hsp90-inhibitor complexes with their analogues were selected as the best complexes based on the GOLD fitness score and orientation. The significant interaction of Hsp90 inhibitors and their analogues were observed to be bound with active site residues as well as residing within the same cavity region. System stability factors RMSD, RMSF, beta-factor, and radius of gyration were analyzed for top-docked complexes and ensure strong binding interaction between inhibitors and the Hsp90 cavity. Cavity bound inhibitors were found to retain consistent hydrogen bonding during the simulation. The radial distribution function (RDF) illustrated that interacting active site residues drive the binding and stability of the inhibitors. Similarly, the axial frequency distribution, which is an indigenously developed analytical tool, produced noteworthy knowledge of the hydrogen-bonding pattern. Results yielded new insights into the design of cancer therapeutic drugs against Hsp90. This finding suggests that under trial Hsp90 inhibitors MPC-3100 could be a potential starting point into the development of potential anticancer agents with the possibility of future directions for the improvement of early existing Hsp90 inhibitors CNF-2024 and SNX-5422 as an anticancer agent.

Insight into natural inhibitors and bridging docking to dynamic simulation against sugar Isomerase (SIS) domain protein

The pathogen Legionella longbeachae is a causative agent of legionellosis. The antibiotic resistance is the major problem of this modern world. Thus, selective pressure warrants the need for identification of newer drug target. In current study, subtractive proteomics approach screen out SIS (sugar isomerase) domain protein as an attractive receptor molecule for rational drug design. This protein is involved in lipopolysaccharide biosynthesis and catalyzes the isomerization of sedoheptulose 7-phosphate in D-glycero-D-manno-heptose 7-phosphate. Molecular docking revealed compound 1 (2-(6-(N,N-dimethyl sulfamoyl)pipridin-4-yl)pyrazin-2-yl)imidazol-3-ium-1-ide) as the potent inhibitor having GOLD fitness score of 69. The complex is affirmed by half-site effect via simulation analysis. Complex stability was investigated via several approaches that follows dynamic simulation and binding energies. Trajectory analysis revealed slight change in ring positioning of inhibitor inside the active pocket during 130 ns (nanosecond). Interestingly, it was affirmed via binding interactions' density distribution. Hence, radial distribution function (RDF) inferred that SER55 and SER83 are the major residues that take part in hydrogen bonding and complex stability. Furthermore, an indigenously developed method axial frequency distribution (AFD) has revealed that ligand moved closer to the active site with both the residues SER55 and SER83 binding to the ligand. The phenomena was observed via rotating motion with respect to receptor center cavity. Thus, inhibitor movement towards allosteric site was observed at the end of simulations. Finally, binding free energy calculations by MMPB/GBSA predicts high compound affinity for the complex. Hence, findings from the current study will aid in the novel drug discovery and future experimental studies. Graphical abstract.

Molecular docking, simulation and MM-PBSA studies of nigella sativa compounds: a computational quest to identify potential natural antiviral for COVID-19 treatment

Nigella sativa or black seed is used as a medicinal plant around the globe. Oil and seeds have a long tradition of folklore use in various medicinal and food systems. The conventional therapeutic use of Nigella sativa, in different ways, has been reported in several studies to treat different diseases including influenza, headache, hypertension, diabetes, inflammation, eczema, fever, cough, asthma, bronchitis, and fever. Based on previously reported potential therapeutic uses of N. sativa compounds, and keeping in mind the dire need of time for the development of potent antiviral, a combined docking, ADMET properties calculation, molecular dynamics, and MM-PBSA approaches were applied in the current study to check the therapeutic potentials of N. sativa chief constituents against COVID-19. Among the studied compounds, we found that dithymoquinone (DTQ), with binding affinity of −8.6 kcal/mol compared to a positive control (chloroquine, −7.2 kcal/mol) , has the high potential of binding at SARS-CoV-2:ACE2 interface and thus could be predicted as a plausible inhibitor to disrupt viral-host interactions. Molecular dynamics simulation of 100 ns well complemented binding affinity of the compound and revealed strong stability of DTQ at the docked site. Additionally, MM-PBSA also affirms the docking results. Compound DTQ of the present study, if validated in wet lab experiments, could be used to treat COVID-19 and could serve as a lead in the future for development of more effective natural antivirals against COVID-19.

Communicated by Ramaswamy H. Sarma

Conformational transition of Acinetobacter baumannii KdsC enzyme and the role of magnesium in binding: An insight from comparative molecular dynamics simulation and its implications in novel antibiotics design

The 3-deoxy-d-manno-octulosonate 8-phosphate phosphatase (KdsC) catalyzes the hydrolysis of 3-deoxy-d-manno-octulosonate 8-phosphate (KDO 8-P) to 3-deoxy-d-manno-octulosonate (KDO) and inorganic phosphate in KDO biosynthesis pathway of Gram-negative bacteria lipopolysaccharide (LPS) hydrophobic lipid-A core. The essentiality of KDO for bacterial cell viability presents the possibility of its targeting to develop broad-spectrum antibacterial agents. In this study, a receptor based virtually screening method was put forward to identify novel lead inhibitory molecules for KdsC enzyme. Dynamics evaluation in solution revealed three complexes: Asinex-1197, Asinex-1705, and Asinex-1710 from Asinex antibacterial library as highly stable, involving conformational transition from open to close upon lead molecules binding and eloquent role of active pocket magnesium towards inhibitors binding and movements. Interconversion of local secondary structure elements in sequence region of Asp192-Asp208 covering motif β-turn, β-hairpin, and β-sheets is seen recurrently that could be in all likelihood of the pressure excreted on this region during closing conformation event or magnesium driven inhibitor adjustments. The binding free energy estimation predicted gas phase energy for all the three complexes dominating with major contribution from van der Waals energy (in case of Asinex-1705 and Asinex-1710) and balanced contributions of both electrostatic and van der Waals (in case of Asinex-1197). Key residues-scanning shortlisted Leu45, Asp185, Gy188, Arg231, and Lys255 as vital in the interaction network of magnesium and inhibitors at the binding site. Their crucial roles in net binding energy were reaffirmed via in silico site directed alanine scanning method. The filtered hits might be useful to further scaffolds addition and structural optimization to yield high affinity binders of KdsC enzyme, whose inhibition, in turn, will disrupt the outer membrane synthesis.

Comparative subtractive proteomics based ranking for antibiotic targets against the dirtiest superbug: Acinetobacter baumannii

Multidrug-resistant Acinetobacter baumannii is indeed to be the most successful nosocomial pathogen responsible for myriad infections in modern health care system. Computational methodologies based on genomics and proteomics proved to be powerful tools for providing substantial information about different aspects of A. baumannii biology that made it possible to design new approaches for treating multi, extensive and total drug resistant isolates of A. baumannii. In this current approach, 35 completely annotated proteomes of A. bauamnnii were filtered through a comprehensive subtractive proteomics pipeline for broad-spectrum drug candidates. In total, 10 proteins (KdsA, KdsB, LpxA, LpxC, LpxD, GpsE, PhoB, UvrY, KdpE and OmpR) could serve as ideal candidates for designing novel antibiotics. The work was extended with KdsA enzyme for structure information, prediction of intrinsic disorders, active site details, and structure based virtual screening of library containing natural product-like scaffolds. Most of the enzyme structure has fixed three-dimensional conformation. The selection of inhibitor for KdsA enzyme was based on druglikeness, pharmacokinetics and docking scores. Compound-4636 (5-((3-chloro-5-methyl-2-phenyl-2,3-dihydro-1H-pyrazol-4-yl)methoxy)-2-(((1-hydroxy-4-methylpentan-2-yl)amino)methyl)phenol) was revealed as the most potent inhibitor against A. baumannii KdsA enzyme having Gold fitness score of 77.68 and Autodock binding energy of -6.2 kcal/mol. The inhibitor completely follows Lipinski rule of five, Ghose rule, and Egan rule. Molecular dynamics simulation for KdsA and KdsA-4636 complex was performed for 100 ns to unveil what conformational changes the enzyme underwent in the absence and presence of the inhibitor, respectively. The average root means square deviation (RMSD) for both systems was found 3.5 Å, which signifies stable structure of the enzyme in both bounded and unbounded states. Absolute binding energy using Molecular Mechanics-Generalized Born Surface Area (MM-GBSA) reflected high affinity and vigorous interactions of the inhibitor with enzyme active residues. Findings of the current study could open up new avenues for experimentalists to design new potent antibiotics by targeting the targets screened in this study.

Molecular dynamics simulation studies of novel β-lactamase inhibitor

New Delhi Metallo-β-Lactamase-1 (NDM-1) has drawn great attention due to its diverse antibiotic resistant activity. It can hydrolyze almost all clinically available β-lactam antibiotics. To inhibit the activity of NDM-1 a new strategy is proposed using computational methods. Molecular dynamics (MD) simulations are used to analyze the molecular interactions between selected inhibitor candidates and NDM-1 structure. The enzyme-ligand complex is subject to binding free energy calculations using MM(PB/GB)SA methods. The role of each residue of the active site contributing in ligand binding affinity is explored using energy decomposition analysis. Furthermore, a hydrogen bonding network between ligand and enzyme active site is observed and key residues are identified ensuring that the ligand stays inside the active site and maintains its movement towards the active site pocket. A production run of 150ns is carried out and results are analyzed using root mean square deviation (RMSD), root mean square fluctuation (RMSF), and radius of gyration (Rg) to explain the stability of enzyme ligand complex. Important active site residue e.g. PHE70, VAL73, TRP93, HIS122, GLN123, ASP124, HIS189, LYS216, CYS208, LYS211, ALA215, HIS250, and SER251 were observed to be involved in ligand attachemet inside the active site pocket, hence depicting its inhibitor potential. Hydrogen bonds involved in structural stability are analyzed through radial distribution function (RDF) and contribution of important residues involved in ligand movement is explained using a novel analytical tool, axial frequency distribution (AFD) to observe the role of important hydrogen bonding partners between ligand atoms and active site residues.