Mechanics Poisson-Boltzmann Surface Area Research Articles

Cefiderocol susceptibility endows hope in treating carbapenem-resistant Pseudomonas aeruginosa: insights from in vitro and in silico evidence.

'High-risk' hypermutable clones of Pseudomonas aeruginosa disseminating extensive drug-resistance (XDR) have raised global health concerns with escalating mortality rates in immunocompromised patients. Mutations in conventional drug-targets under antibiotic stress necessitate structural understanding to formulate sustainable therapeutics. In the present study, the major β-lactam antibiotic target, penicillin-binding protein-3 (PBP3) with mutations F533L and T91A, were identified in carbapenemase-positive P. aeruginosa isolates (n = 6) using whole genome sequencing. Antibiotic susceptibility tests showed susceptibility to cefiderocol (MIC ≤ 4 μg ml-1) despite pan-β-lactam resistance in the isolates. Both the mutations reduced local intra-chain interactions in PBP3 that marginally increased the local flexibility (∼1%) in the structures to affect antibiotic-interactions. Molecular dynamics simulations confirmed the overall stability of the PBP3 mutants through root-mean square deviations, radius of gyration, solvent-accessibility and density curves, which favored their selection. Docking studies unveiled that the mutations in PBP3 elicited unfavorable stereochemical clashes with the conventional antibiotics thereby increasing their inhibition constants (IC) up to ∼50 fold. It was deciphered that cefiderocol retained its susceptibility despite mutations in PBP3, due to its higher average binding affinity (ΔG: -8.2 ± 0.4 kcal mol-1) towards multiple PBP-targets and lower average binding affinity (ΔG: -6.7 ± 0.7 kcal mol-1) to β-lactamases than the other β-lactam antibiotics. The molecular dynamics simulations and molecular mechanics Poisson Boltzmann surface area calculations further indicated energetically favorable binding for cefiderocol with PBP3 proteins. The study gave structural insight into emerging non-polar amino acid substitutions in PBP3 causing XDR and recommends prioritizing available antibiotics based on multi-target affinities to overcome challenges imposed by target-protein mutations.

Biophysical and structural characterization of tetramethrin serum protein complex and its toxicological implications.

Tetramethrin (TMT) is a commonly used insecticide and has a carcinogenic and neurodegenerative effect on humans. The binding mechanism and toxicological implications of TMT to human serum albumin (HSA) were examined in this study employing a combination of biophysical and computational methods indicating moderate binding affinity and potential hepato and renal toxicity. Fluorescence quenching experiments showed that TMT binds to HSA with a moderate affinity, and the binding process was spontaneous and predominantly enthalpy-driven. Circular dichroism spectroscopy revealed that TMT binding did not induce any significant conformational changes in HSA, resulting in no changes in its alpha-helix content. The binding site and modalities of TMT interactions with HSA as computed by molecular docking and molecular dynamics simulations revealed that it binds to Sudlow site II of HSA via hydrophobic interactions through its dimethylcyclopropane carboxylate methyl propanyl group. The structural dynamics of TMT induce proper fit into the binding site creating increased and stabilizing interactions. Additionally, molecular mechanics-Poisson Boltzmann surface area calculations also indicated that non-polar and van der Waals were found to be the major contributors to the high binding free energy of the complex. Quantum mechanics (QM) revealed the conformational energies of the binding confirmation and the degree of deviation from the global minimum energy conformation of TMT. The results of this study provide a comprehensive understanding of the binding mechanism of TMT with HSA, which is important for evaluating the toxicity of this insecticide in humans.

Computational study of the piperidine and FtsZ interaction in Salmonella Typhi: implications for disrupting cell division machinery

FtsZ, a bacterial cell division protein, is essential for assembling the contractile Z-ring crucial in bacterial cytokinesis. Consequently, inhibiting FtsZ could impede proto-filaments, disrupting FtsZ and other associated proteins vital for cell division machinery. Conduct an in-silico drug interaction study to identify novel drug candidates that inhibit the FtsZ protein, aiming to prevent Multi-Drug Resistant (MDR) Salmonella Typhi. Data mining was performed based on piperidine compounds, which were subsequently screened for safe pharmacokinetic profiles. Compounds that met favorable drug-likeness criteria underwent virtual screening against the FtsZ drug target. Two compounds were chosen for molecular docking and molecular dynamic simulation to verify the binding affinity and stability between the target protein and the potential compounds. The 400 isoforms of piperidine analogues were curated, among them potent compound ZINC000000005416 found to possess high binding affinity (-8.49 kcal/mol) and low dissociation constant (0.597 µM). The highest binding affinity shown by ZINC000000005416 was validated by hydrogen bonds, hydrophobic interaction, and salt bridges with the functional domain of the cell division regulatory protein. Docking profiles, when correlated with molecular dynamic simulation (MDS) depicted stable trajectories and compatible conformational changes in the FtsZ-ZINC000000005416 complex. The stable simulated trajectories were validated through free-energy calculations using the Molecular Mechanics-Poisson Boltzmann Surface Area (MM/PBSA) module. Low energy conformations, although the simulation trajectory confirmed the stable ZINC000000005416-FtsZ interaction, which encouraged experimental validations. This study encourages further exploration of the compound ZINC000000005416 as a drug candidate inhibiting FtsZ protein against MDR Salmonella Typhi. Communicated by Ramaswamy H. Sarma

In Silico Investigation on the Molecular Behavior and Structural Stability of the Rosette Nanotubes as the Drug Vehicles for Paclitaxel, an Anti-Cancer Drug

Most anticancer drugs affect healthy cells in addition to cancer cells, causing severe side effects. Targeted delivery by nano-based drug delivery systems (NDDS) can reduce these severe side effects while maintaining therapeutic efficacy. This work introduced rosette nanotube (RNT) as a potential drug vehicle for paclitaxel (PTX) due to its self-assembling property, biocompatibility, amphiphilicity, and low toxicity. Molecular dynamics (MD) simulations aided with molecular mechanics Poisson Boltzmann surface area (MMPBSA) analysis are used here to investigate the molecular behavior and the loading energetics of each type of RNT (K1, xK1, and iEt-xK1) with PTX. Analysis showed that the most probable configuration of PTX is on either end of each RNT. The binding free energies (−117.74 to −69.29 kJ/mol) when PTX is closer to one end were stronger than when it is in the inner channel (−53.51 to −40.88 kJ/mol). The latter alludes to the encapsulation of the PTX by each RNT. Thus, loading is possible by encapsulation during the self-assembly process given the favorable estimated binding free energies. Based on the results, RNT has potential as a drug vehicle for PTX, which warrants further investigation.

Computational Screening of Neuropilin-1 Unveils Novel Potential Anti-SARS-CoV-2 Therapeutics.

Neuropilin 1 (NRP-1) inhibition has shown promise in reducing the infectivity of severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2) and preventing the virus entry into nerve tissues, thereby mitigating neurological symptoms in COVID-19 patients. In this study, we employed virtual screening, including molecular docking, Molecular Dynamics (MD) simulation, and Molecular Mechanics-Poisson Boltzmann Surface Area (MM-PBSA) calculations, to identify potential NRP-1 inhibitors. From a compendium of 1930 drug-like natural compounds, we identified five potential leads: CNP0435132, CNP0435311, CNP0424372, CNP0429647, and CNP0427474, displaying robust binding energies of -8.2, -8.1, -10.7, -8.2, and -8.2 kcal/mol, respectively. These compounds demonstrated interactions with critical residues Tyr297, Trp301, Thr316, Asp320, Ser346, Thr349, and Tyr353 located within the b1 subdomain of NRP-1. Furthermore, MD simulations and MM-PBSA calculations affirmed the stability of the complexes formed, with average root mean square deviation, radius of gyration, and solvent accessible surface area values of 0.118 nm, 1.516 nm, and 88.667 nm2 , respectively. Notably, these lead compounds were estimated to penetrate the blood-brain barrier and displayed antiviral properties, with Pa values ranging from 0.414 to 0.779. The antagonistic effects of these lead compounds merit further investigation, as they hold the potential to serve as foundational scaffolds for the development of innovative therapeutics aimed at reducing the neuroinfectivity of SARS-CoV-2.

Molecular level investigation for identifying potential inhibitors against thymidylate kinase of monkeypox through in silico approaches.

The need for more advanced and effective monkeypox (Mpox) treatments has become evident with numerous Mpox virus (MPXV) outbreaks. Over the years, interest has increased in developing targeted medicines that are efficient, safe, and precise while avoiding adverse effects. Here, we screened 32409 compounds against thymidylate kinase (TMPK), an emerging target for Mpox treatment. We studied their pharmacological characteristics and analyzed those through all-atom molecular dynamics simulations followed by molecular mechanics Poisson Boltzmann surface area (MM-PBSA) based free energy calculations. According to our findings, the leads CID40777874 and CID28960001 had the highest binding affinities towards TMPK with ΔGbind of -8.04 and -5.58 kcal/mol, respectively, which outperformed our control drug cidofovir (ΔGbind = -2.92 kcal/mol) in terms of binding favourability. Additionally, we observed crucial TMPK dynamics brought on by ligand-binding and identified key residues such as Phe68 and Tyr101 as the critical points of the protein-ligand interaction. The DCCM analysis revealed the role of ligand binding in stabilizing TMPK's binding region, as indicated by residual correlation motions. Moreover, the PSN analysis revealed that the interaction with ligand induces changes in residual network properties, enhancing the stability of complexes. We successfully identified novel compounds that may serve as potential building blocks for constructing contemporary antivirals against MPXV and highlighted the molecular mechanisms underlying their binding with TMPK. Overall, our findings will play a significant role in advancing the development of new therapies against Mpox and facilitating a comprehensive understanding of their interaction patterns.Communicated by Ramaswamy H. Sarma.

Unlocking potential inhibitors for Bruton's tyrosine kinase through in-silico drug repurposing strategies

Bruton's tyrosine kinase (BTK) is a non-receptor protein kinase that plays a crucial role in various biological processes, including immune system function and cancer development. Therefore, inhibition of BTK has been proposed as a therapeutic strategy for various complex diseases. In this study, we aimed to identify potential inhibitors of BTK by using a drug repurposing approach. To identify potential inhibitors, we performed a molecular docking-based virtual screening using a library of repurposed drugs from DrugBank. We then used various filtrations followed by molecular dynamics (MD) simulations, principal component analysis (PCA), and Molecular Mechanics Poisson Boltzmann Surface Area (MM-PBSA) analysis to further evaluate the binding interactions and stability of the top-ranking compounds. Molecular docking-based virtual screening approach identified several repurposed drugs as potential BTK inhibitors, including Eltrombopag and Alectinib, which have already been approved for human use. All-atom MD simulations provided insights into the binding interactions and stability of the identified compounds, which will be helpful for further experimental validation and optimization. Overall, our study demonstrates that drug repurposing is a promising approach to identify potential inhibitors of BTK and highlights the importance of computational methods in drug discovery.

In silico and in vitro investigation of dual targeting Prima-1MET as precision therapeutic against lungs cancer

This study emphasizes the explorations of binding of Prima-1MET with two targets, p53 a tumor suppressor protein, and tyrosine kinase of epidermal growth factor receptor. In silico investigations reveal that Prima-1MET showed robust binding with both targets. Molecular docking simulations demonstrated the binding affinity of Prima-1MET with p53 and tyrosine kinase was found to be −38.601 kJ/mol and −38.976 kJ/mol. In addition, the stability of Prima-1MET was explored by molecular dynamics simulation. Prima-1MET attains stability in the binding site of the respective protein till the simulation period is over. Moreover, the free binding energy ΔG bind was calculated by the molecular mechanics Poisson Boltzmann surface area method. The ΔG bind of Prima-1MET with tyrosine kinase was found to be −58.585 ± 0.327 kJ/mol and with p53 it was −35.910 ± 0.335 kJ/mol. Next, cytotoxicity of the Prima-1MET was evaluated using multiple cancer cell lines and the IC50 value were ranging between 4.5 and 30 μM. The cell death was identified by apoptosis assay. Further, the p53 and tyrosine kinase expression was monitored using immunofluorescence techniques, it was found Prima-1MET induces the expression of p53 protein and mimics the level of tyrosine kinase oncogenic target. Also, reactive oxygen species (ROS) and membrane potential activity of Prima-1MET was evaluated by using a lung cancer cell line. A significant decrease in intracellular ROS was observed and resulted in disruption of mitochondrial transmembrane potential. This study uncovers the underlying mechanism of Prima-1MET and could be helpful to design further leads against lung cancers. Communicated by Ramaswamy H. Sarma

“Identification of alkaloid compounds as potent inhibitors of Mycobacterium tuberculosis NadD using computational strategies”

Mycobacterium tuberculosis is leading cause of death worldwide. NAD participates in a host of redox reactions in energy landscape of organisms. Several studies implicate surrogate energy pathways involving NAD pools as important in survival of active as well as dormant mycobacteria. One of the NAD metabolic pathway enzyme, nicotinate mononucleotide adenylyltransferase (NadD) is indispensable in mycobacterial NAD metabolism and is perceived as an attractive drug target in pathogen. In this study, we have employed in silico screening, simulation and MM-PBSA strategies to identify potentially important alkaloid compounds against mycobacterial NadD for structure-based inhibitor development. We have performed an exhaustive structure-based virtual screening of an alkaloid library, ADMET, DFT profiling followed by Molecular Dynamics (MD) simulation, and Molecular Mechanics-Poisson Boltzmann Surface Area (MM-PBSA) calculation to identify 10 compounds which exhibit favourable drug like properties and interactions. Interaction energies of these 10 alkaloid molecules range between −190 kJ/mol and -250 kJ/mol. These compounds could be promising starting point in the development of selective inhibitors against Mycobacterium tuberculosis.

In silico screening of phytochemicals from Dissotis rotundifolia against Plasmodium falciparum Dihydrofolate Reductase

BackgroundMalaria remains a major health concern in developing countries with high morbidity and mortality, especially in pregnant women and infants. A major obstacle to the treatment of malaria is a low effectiveness and an increase resistance of the parasite to antimalarial drugs. As a result, there is an ongoing demand for new and potent antimalarial drugs. Medicinal plants remain a potential source for the development of new antimalarial drugs. Amongst them is Dissotis rotundifolia is an ethnomedical important plant used in West Africa to treat malaria. PurposeThis study aimed at identifying new potential antifolates by virtually screening phytochemicals characterized from the whole plant methanolic extract of D. rotundifolia against Plasmodium falciparum Dihydrofolate Reductase (PfDHFR). MethodsLC-ESI-Q-TOF-MS analysis was employed to identify the phytochemicals present in the whole plant methanolic extract of D. rotundifolia. These phytochemicals were docked against the catalytic site of PfDHFR. The docking protocol was evaluated using the Area Under the Curve (AUC) of a Receiver Operating Characteristic (ROC) curve. The binding mechanisms and the drug-likeness of the phytochemicals were characterized. A 100 ns Molecular Dynamics (MD) simulation and Molecular Mechanics-Poisson Boltzmann Surface Area (MM-PBSA) calculations were utilized to analyze the stability, the energy decomposition per residue and the binding free energy of the potential leads. ResultsTwenty nine phytochemicals were characterized and docked against PfDHFR. Dimethylmatairesinol, flavodic acid, sakuranetin, and sesartemin were identified as potential leads with binding affinities of -8.4, -8.9, -8.6, and -8.9 kcal/mol respectively, greater than a stringent threshold of -8.0 kcal/mol. The potential leads also interacted hydrophobically with critical residue Phe58. A novel critical residue, Leu46 was identified to be crucial in the catalytic activity of PfDHFR. The potential leads were also predicted to be anti-protozoal with a probability of active (Pa) value ranging from 0.319 to 0.537. ConclusionThis study elucidates the potential inhibition of PfDHFR by dimethylmatairesinol, flavodic acid, sakuranetin and sesartemin present in D. rotundifolia. These compounds are druglike, do not violate Lipinski's rule of five, have a high binding affinity to PfDHFR, and interact with crucial residues involved in the catalytic activity PfDHFR. Dimethylmatairesinol, flavodic acid, sakuranetin and sesartemin could therefore be further investigated and developed as new antifolate drugs for malaria.

2D-QSAR, molecular docking and MD simulation based virtual screening of the herbal molecules against Alzheimer’s disorder: an approach to predict CNS activity

Acetylcholinesterase (AChE) is one of the key enzyme targets that have been used clinically for the management of Alzheimer’s Disorder (AD). Numerous reports in the literature predict and demonstrate in-vitro, and in-silico anticholinergic activity of herbal molecules, however, majority of them failed to find clinical application. To address these issues, we developed a 2D-QSAR model that could efficiently predict the AChE inhibitory activity of herbal molecules along with predicting their potential to cross the blood-brain-barrier (BBB) to exert their beneficial effects during AD. Virtual screening of the herbal molecules was performed and amentoflavone, asiaticoside, astaxanthin, bahouside, biapigenin, glycyrrhizin, hyperforin, hypericin, and tocopherol were predicted as the most promising herbal molecules for inhibiting AChE. Results were validated through molecular docking, atomistic molecular dynamics simulations and Molecular mechanics-Poisson Boltzmann surface area (MM-PBSA) studies against human AChE (PDB ID: 4EY7). To determine whether or not these molecules can cross BBB to inhibit AChE within the central nervous system (CNS) for being beneficial for the management of AD, we determined a CNS Multi-parameter Optimization (MPO) score, which was found in the range of 1 to 3.76. Overall, the best results were observed for amentoflavone and our results demonstrated a PIC50 value of 7.377 nM, molecular docking score of −11.5 kcal/mol, and CNS MPO score of 3.76. In conclusion, we successfully developed a reliable and efficient 2D-QSAR model and predicted amentoflavone to be the most promising molecule that could inhibit human AChE enzyme within the CNS and could prove beneficial for the management of AD. Communicated by Ramaswamy H. Sarma

Development and test of highly accurate endpoint free energy methods. 2: Prediction of logarithm of n-octanol-water partition coefficient (logP) for druglike molecules using MM-PBSA method.

The logarithm of n-octanol-water partition coefficient (logP) is frequently used as an indicator of lipophilicity in drug discovery, which has substantial impacts on the absorption, distribution, metabolism, excretion, and toxicity of a drug candidate. Considering that the experimental measurement of the property is costly and time-consuming, it is of great importance to develop reliable prediction models for logP. In this study, we developed a transfer free energy-based logP prediction model-FElogP. FElogP is based on the simple principle that logP is determined by the free energy change of transferring a molecule from water to n-octanol. The underlying physical method to calculate transfer free energy is the molecular mechanics-Poisson Boltzmann surface area (MM-PBSA), thus this method is named as free energy-based logP (FElogP). The superiority of FElogP model was validated by a large set of 707 structurally diverse molecules in the ZINC database for which the measurement was of high quality. Encouragingly, FElogP outperformed several commonly-used QSPR or machine learning-based logP models, as well as some continuum solvation model-based methods. The root-mean-square error (RMSE) and Pearson correlation coefficient(R) between the predicted and measured values are 0.91 log units and 0.71, respectively, while the runner-up, the logP model implemented in OpenBabel had an RMSE of 1.13 log units and R of 0.67. Given the fact that FElogP was not parameterized against experimental logP directly, its excellent performance is likely to be expanded to arbitrary organic molecules covered by the general AMBER force fields.

In vitro contraceptive activities, molecular docking, molecular dynamics, MM-PBSA, non-covalent interaction and DFT studies of bioactive compounds from Aegle marmelos. Linn., leaves.

Introduction: Bioactive molecules from natural sources having contraceptive properties were excellent alternatives for modern hormonal contraceptives. Researchers around the world were working on identifying contraceptive leads targeting the male reproductive system rather than the usual female contraceptives. The lack of proper understanding on male contraceptive protein drug targets leads to insufficient evidence on activities of identified contraceptive compounds. The proteins specific to the male reproductive system and involved in sperm-egg fusion will be an excellent drug target to identify the male non-hormonal, reversible contraceptive leads. Inhibiting sperm hyaluronidase activity by natural non-hormonal compounds will lead to reversible and non-hormonal male contraception. The Aegle marmelos Linn. is one such important medicinal plant with valuable phytocompounds, used traditionally as a potential contraceptive measure. The in vivo experiments on leaf extracts of Aegle marmelos. Linn containing terpenes, sterols, and alkaloids shows prominent contraceptive activities. Moreover, this study explores the potential ability of the leaf extract on inhibiting the sperm hyaluronidase action with additional molecular details on the interaction between sperm hyaluronidases and three phytocompounds such as aegeline, marmin, and marminol. Material and methods: The in vitro hyaluronidase inhibition assay and Computer Assisted Sperm Analysis (CASA) were used to evaluate the male contraceptive properties of the Aegle marmelos Linn. leaf extract. To identify the interaction profile of aegeline, marmin, and marmenol on sperm cell hyaluronidases the in-silico methods such as molecular docking, Non-Covalent Interaction analysis, Molecular dynamics, and Molecular Mechanics Poisson Boltzmann Surface Area were used. Results and discussion: The results of in vitro hyaluronidase inhibition assay and Computer Assisted Sperm Analysis shows the inhibition of hyaluronidase enzymatic activity and reduced sperm activities in the presence of leaf extracts. After incubation with leaf extracts for about 30 minutes time intervals show, the motility drops from progressive to non-progressive and ended up with complete immotile in 100 μg/ml concentration of leaf extract. The results of molecular docking, Non-Covalent Interaction analysis, Molecular dynamics, and Molecular mechanics Poisson Boltzmann Surface Area show that the phytocompounds marmin, and aegeline have the potential ability to inhibit sperm hyaluronidase.

Identification of potent EGFR-TKD inhibitors from NPACT database through combined computational approaches

Cancer is the world’s second leading cause of death, and there are no approved herbal therapies. The epidermal growth factor receptor tyrosine kinase (EGFR-TK) receptor is a transmembrane protein with eight domains that is found in almost every cancer type and plays an important role in abnormal cell cellular function and causes malignant outcomes. The current study aimed to virtually screen phytochemicals from the NPACT database against EGFR-TKD and also to identify potential inhibitors of this transmembrane protein among plant candidates for anticancer drug development. The docking scores of the chosen phytochemicals were compared with the control (erlotinib). Kurarinone, (2S)-2-methoxykurarnione, and Sophoraflavanone-G exhibited a stronger binding affinity of −18.102 kcal/mol, −14.243 kcal/mol, and −13.759 kcal/mol than erlotinib −12.783 kcal/mol. Moreover, several online search engines were used to predict ADME and toxicity. The drug-likeness of selected phytochemicals was higher than the reference (erlotinib). A 100 ns molecular dynamic (MD) simulation was also applied to the docked conformations to examine the stability and molecular mechanics of protein-ligand interactions. Furthermore, the calculated molecular mechanics Poisson Boltzmann surface area energy of (2S)-2-methoxykurarnione was found to be −129.555 ± 0.512 kJ/mol, which approximately corresponds to the free energy of the reference molecule −130.595 ± 0.908 kJ/mol. We identify phytoconstituents present in Sophora flavescens from the NPACT database, providing key insights into tyrosine kinase inhibition and may serve as better chemotherapeutic agents. Experimental validation is required to determine the anti-EGFR potency of the potent lead molecules discussed in this study. Communicated by Ramaswamy H. Sarma

In silico identification of potential PvFKBP35 inhibitors from Entadrophragma angolense Limonoids extracts as antimalarial agents

Plasmodium species, which are spread by female Anopheles mosquitoes, are responsible for malaria. Out of the five major Plasmodium species, Plasmodium falciparum and Plasmodium vivax are the most deadly and invasive species responsible for 99.7% and 75% of malaria cases in Africa and America respectively. Despite the invasive nature of malaria, the Plasmodium parasite continues to develop resistance to current drugs. It is therefore imperative to come up with new therapeutics to combat malaria. Previous studies have reported that Limonoids from the Meliaceae family possess antimalarial properties. This study therefore aims at employing computational approaches to identify potential antimalarial Limonoids by targeting PvFKBP35. PvFKBP35 has been reported to be a suitable target for antimalarial therapeutics as it is involved in various physiological activities including transcription, protein stability and folding. Molecular docking, Molecular Dynamics simulation and Molecular Mechanics-Poisson Boltzmann Surface Area calculation were employed to identify the potential leads. Sixteen [16] Limonoids extracted from the bark of the stem of Entadrophragma angolense were virtually screened against PvFKPB35. The top hit compounds were subjected to 500 ns Molecular Dynamics simulation and Molecular Mechanics – Poisson Boltzmann Surface Area calculations to examine their stability and free binding energy. Two potential leads, compounds 1 and 11 with binding energies −6.3 and −5.4 kcal/mol respectively were identified. The potential leads in complexed with PvFKBP35 had an average root mean square deviation of 1.18 ± 0.19 Å and 3.12 ± 0.60 Å, indicating their stability. Solvent Accessible Surface Area was utilized to predict the penetrative ability of the compounds into the binding pocket. Average Solvent Accessible Surface Area values of 327.88 ± 47.54 A2, 402.18 ± 39.81 A2 were obtained for compounds 1 and 11 respectively. ADMET estimations of compounds 1 and 11 predicted them to be druglike and do not violate Lipinski's rule of five. Compounds 1 and 11 need be tested in vitro to validate their antimalarial activity although they were predicted to be antiprotozoal with Pa values 0.207 and 0.162. These compounds can then serve as the scaffold for the design of novel antimalarial therapeutics.

Molecular docking and dynamics studies show: Phytochemicals from Papaya leaves extracts as potential inhibitors of SARS–CoV–2 proteins targets and TNF–alpha and alpha thrombin human targets for combating COVID-19

<abstract> <p>Tackling COVID-19 requires halting virus proliferation and reducing viral complications in humans. Papaya leaf extract (PLE) is well known for its ability to inhibit numerous viral replications <italic>in vitro</italic> and <italic>in vivo</italic> and reduce viral complications in humans such as thrombocytopenia and cytokine storm. The goal of this research is to evaluate the possible use of papaya leaf extract as a multifaceted antiviral and potential therapy for COVID-19 using an in-silico docking followed by a 100 ns molecular dynamics simulation (MDS) approach. The targeted proteins are the SARS-CoV-2's proteins such as the nucleocapsid, main protease (MPro), RNA-dependent RNA polymerase (RdRP), spike protein (Wuhan, Delta, and Omicron variants) and human TNF-alpha and alpha-thrombin protein targets. Several compounds from PLE such as protodioscin, clitorin, glycyrrhizic acid, manghaslin, kaempferol–3–(2g–glucosylrutinoside), rutin, isoquercetrin and acacic acid were found to exhibit strong binding to these targets. The free energies of binding (Autodock) with protodioscin, the best PLE compound for nucleocapsid, main protease (MPro), RdRP and spike protein were –13.83, –13.19, –11.62 and –10.77 (Omicron), kcal/mol, respectively, while the TNF-alpha and alpha-thrombin binding free energies were –13.64 and –13.50 kcal/mol, respectively. The calculated inhibition constants for protodioscin were in the nanomolar and picomolar range at 216.34, 27.07, 73.28, and 99.93 pM, respectively, whilst RdRp and spike protein (Omicron) were in the nanomolar range at 3.02 and 12.84 nM, respectively. Protodioscin interacted with key residues of all protein targets. The binding affinity poses were confirmed by molecular dynamics simulation. Analysis of the binding affinities calculated employing the molecular mechanics-Poisson–Boltzmann surface area (MM-PBSA) shows favorable interaction between protodioscin, and all targets based on total binding-free energies corroborating the Autodock's docking results. In conclusion, compounds from PLE, especially protodioscin have good potentials in combating COVID-19.</p> </abstract>

Identification of potent COVID-19 main protease inhibitors by loading of favipiravir on Mg12O12 and Zn12O12 nanoclusters: an in silico strategy for COVID-19 treatment

Pandemic new severe acute respiratory syndrome coronavirus (SARS-CoV-2) virus has increased throughout the world. There is no effective treatment against this virus until now. Since its appearance in Wuhan, China in December 2019, SARS-CoV-2 becomes the largest challenge the world is opposite today, including the discovery of an antiviral drug for this virus. Several viral proteins have been prioritized as SARS-CoV-2 antiviral drug targets, among them the papain-like protease (PLpro) and the main protease (Mpro). Inhibition of these proteases would target viral replication, viral maturation and suppression of host innate immune responses. Potential candidates have been identified to show inhibitory effects against Mpro, both in biochemical assays and viral replication in cells. There are different molecules such as lopinavir and favipiravir considerably inhibit the activity of Mpro in vitro. Different studies have shown that structurally improved favipiravir and other similar compounds can inhibit SARS-CoV-2 main protease. In this work, we study the interactions between favipiravir with Mg12O12 and Zn12O12 nanoclusters by density functional theory (DFT) and quantum mechanics atoms in molecules (QMAIM) methods to summarize the ability to load favipiravir onto Mg12O12 and Zn12O12 nanoclusters. Favipiravir-Mg12O12 and favipiravir-Zn12O12 lowest structures complexes were chosen to dock inside the SARS-CoV-2 main protease by molecular docking study. The molecular docking analysis revealed that the binding affinity of Mg12O12 and Zn12O12 nanoclusters inside the Mpro receptor is larger than that of favipiravir. Also, the loading of favipiravir on the surface of Mg12O12 and Zn12O12 nanoclusters increased the binding affinity against the Mpro receptor. Subsequently, 100 ns molecular dynamics simulation of the favipiravir-Mg12O12, and favipiravir-Zn12O12 docked inside the Mpro complexes established that favipiravir-Mg12O12, forms the most stable complex with the Mpro. Further molecular mechanics Poisson Boltzmann surface area (MMPBSA) analyses using the MD trajectories also demonstrated the higher binding affinity of favipiravir-Mg12O12 inside the Mpro. In summary, this study demonstrates a new way to characterize leads for novel anti-viral drugs against SARS-CoV-2, by improving the drug ability of favipiravir via loading it on Mg12O12 and Zn12O12 nanoclusters. Communicated by Ramaswamy H. Sarma

Finding inhibitors and deciphering inhibitor-induced conformational plasticity in the Janus kinase via multiscale simulations

ABSTRACT The Janus kinase (JAK) is a master regulator of the JAK/STAT pathway. Dysregulation of this signalling cascade causes neuroinflammation and autoimmune disorders. Therefore, JAKs have been characterized as an attractive target for developing anti-inflammatory drugs. Nowadays, designing efficient, effective, and specific targeted therapeutics without being cytotoxic has gained interest. We performed the virtual screening of natural products in combination with pharmacological analyses. Subsequently, we performed molecular dynamics simulations to study the stability of the ligand-bound complexes and ligand-induced inactive conformations. Notably, inactive kinases display remarkable conformational plasticity; however, ligand-induced molecular mechanisms of these conformations are still poorly understood. Herein, we performed a free energy landscape analysis to explore the conformational plasticity of the JAK1 kinase. Leonurine, STOCK1N-68642, STOCK1N-82656, and STOCK1N-85809 bound JAK1 exhibited a smooth transition from an active (αC-in) to a completely inactive conformation (αC-out). Ligand binding induces disorders in the αC-helix. Molecular mechanics Poisson Boltzmann surface area (MM/PBSA) calculation suggested three phytochemicals, namely STOCK1N-68642, Epicatechin, and STOCK1N-98615, have higher binding affinity compared to other ligand molecules. The ligand-induced conformational plasticity revealed by our simulations differs significantly from the available crystal structures, which might help in designing allosteric drugs.

Pharmacological Basis of Rumex hastatus D. Don in Gastrointestinal Diseases with Focusing Effects on H+/K+-ATPase, Calcium Channels Inhibition and PDE Mediated Signaling: Toxicological Evaluation on Vital Organs.

This present study aimed to delineate Rumex hastatus D. Don crude extract (Rh.Cr), n-Hexane, ethyl acetate, aqueous fractions (Rh.n-Hex, Rh.ETAC, Rh.Aq) and rutin for antidiarrheal, antisecretory effects, anti-spasmodic, gastrointestinal transient time, anti H. pylori, antiulcer effects, and toxicology. The preliminary phytochemical analysis of Rumex hastatus showed different phytoconstituents and shows different peaks in GC-MC chromatogram. Rumex hastatus crude extract (Rh.Cr), fractions, and rutin attributed dose-dependent (50–300 mg/kg) protection (0–100%) against castor oil-induced diarrhea and dose-dependently inhibited intestinal fluid secretions in mice. They decreased the distance traversed by charcoal in the gastrointestinal transit model in rats. In rabbit jejunum preparations, Rh.Cr and Rh.ETAC caused a concentration-dependent relaxation of both spontaneous and K+ (80 mM)-induced contractions at a similar concentration range, whereas Rh.n-Hex, rutin, and verapamil were relatively potent against K+-induced contractions and shifted the Ca2+ concentration–response curves (CRCs) to the right, Rh.Cr (0.3–1 mg/mL) and Rh.ETAC (0.1–0.3 mg/mL) shifted the isoprenaline-induced inhibitory CRCs to the left. Rh.n-Hex, Rh.ETAC and rutin showed anti-H. pylori effect, also shows an inhibitory effect against H+/K+-ATPase. Rumex hastatus showed gastroprotective and antioxidant effects. Histopathological evaluation showed improvement in cellular architecture and a decrease in the expression of inflammatory markers such as, cyclooxygenase (COX-2), tumor necrosis factor (TN,F-α) and phosphorylated nuclear factor kappa B (p-NFƙB), validated through immunohistochemistry and ELISA techniques. In RT-PCR it decreases H+/K+-ATPase mRNA levels. Rumex hastatus was found to be safe to consume up to a dose of 2000 mg/kg in a comprehensive toxicity profile. Docking studies revealed that rutin against H+/K+-ATPase pump and voltage-gated L-type calcium channel showed E-values of −8.7 and −9.4 Kcal/mol, respectively. MD simulations 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.

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 Ca2+ 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.