Molecular Mechanics Poisson-Boltzmann Surface 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.

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.

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

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

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.

Structural requirements and interaction mechanisms of ACE inhibitory peptides: molecular simulation and thermodynamics studies on LAPYK and its modified peptides

The understanding of the structural requirements and the intermolecular-interaction mechanism are important for discovering potent angiotensin-converting enzyme (ACE) inhibitory peptides. In this study, we modified an egg-white derived peptide, LAPYK, using the amino acids with different properties to produce the LAPYK-modified peptides. The ACE inhibitory activities of the modified peptides were determined to explore the structural requirements of ACE inhibitory peptides (ACEIPs). Molecular simulation and isothermal titration calorimetry analysis were used to investigate interactions between the peptides and ACE. We found that hydrophobicity and the amino acids with ring structures were beneficial for the ACE inhibitory activities of the peptides. The results of the molecular mechanics poisson boltzmann surface area (MMPBSA) binding free energy calculations indicated that the polar solvation free energy (ΔGpolar) of the charged peptides (LAPYK, LAPYE) were unfavorable for binding to ACE. On the other hand, the results of isothermal titration calorimetry analyses suggested that the enthalpy-driven ACE-peptide interactions were more favorable than the entropy-driven ACE-peptide interaction counterparts.

Molecular Dynamics and MM-PBSA Analysis of the SARS-CoV-2 Gamma Variant in Complex with the hACE-2 Receptor.

The SARS-CoV-2 virus, since its appearance in 2019, has caused millions of cases and deaths. To date, there is no effective treatment or a vaccine that is fully protective. Despite the efforts made by governments and health institutions around the globe to control its propagation, the evolution of the virus has accelerated, diverging into hundreds of variants. However, not all of them are variants of concern (VoC’s). VoC’s have appeared in different regions and throughout the two years of the pandemic they have spread around the world. Specifically, in South America, the gamma variant (previously known as P.1) appeared in early 2021, bringing with it a second wave of infections. This variant contains the N501Y, E484K and K417T mutations in the receptor binding domain (RBD) of the spike protein. Although these mutations have been described experimentally, there is still no clarity regarding their role in the stabilization of the complex with the human angiotensin converting enzyme 2 (hACE-2) receptor. In this article we dissect the influence of mutations on the interaction with the hACE-2 receptor using molecular dynamics and estimations of binding affinity through a screened version of the molecular mechanics Poisson Boltzmann surface area (MM-PBSA) and interaction entropy. Our results indicate that mutations E484K and K417T compensate each other in terms of binding affinity, while the mutation N501Y promotes a more convoluted effect. This effect consists in the adoption of a cis configuration in the backbone of residue Y495 within the RBD, which in turn promotes polar interactions with the hACE-2 receptor. These results not only correlate with experimental observations and complement previous knowledge, but also expose new features associated with the specific contribution of concerned mutations. Additionally, we propose a recipe to assess the residue-specific contribution to the interaction entropy.

Discovery of potential competitive inhibitors against With-No-Lysine kinase 1 for treating hypertension by virtual screening, inverse pharmacophore-based lead optimization, and molecular dynamics simulations

ABSTRACT The With-No-Lysine (WNK) has received attention because of its involvement in hypertension. Genetic mutation in the genes of WNK, leading to its overexpression, has been reported in Familial Hyperkalaemic Hypertension, and thus WNK is considered a potential drug target. Herein, we have performed a high-throughput virtual screening of ~11,000 compounds, mainly the natural phytochemical compounds and kinase inhibitory libraries, to find potential competitive inhibitors against WNK1. Initially, candidates with a docking score of ~ −10.0 kcal/mol or less were selected to further screen their good pharmacological properties by applying absorption, distribution, metabolism, excretion, and toxicity (ADMET). Finally, six docked compounds bearing appreciable binding affinities and WNK1 selectivity were complimented with 500 ns long all-atom molecular dynamic simulations. Subsequently, the MMPBSA scheme (Molecular Mechanics Poisson Boltzmann Surface Area) suggested three phytochemical compounds, C00000947, C00020451, and C00005031, with favourable binding affinity against WNK1. Among them, C00000947 acts as the most potent competitive inhibitor of WNK1. Further, inverse pharmacophore-based lead optimization of the C00000947 leads to one potential compound, meciadanol, which shows better binding affinity and specificity than C00000947 towards WNK1, which may be further exploited to develop effective therapeutics against WNK1-associated hypertension after in vitro and in vivo validation.

Pharmacophore screening to identify natural origin compounds to target RNA-dependent RNA polymerase (RdRp) of SARS-CoV2.

Several existing drugs have gained initial consideration due to their therapeutic characteristics against COVID-19 (Corona Virus Disease 2019). Hydroxychloroquine (HCQ) was proposed as possible therapy for shortening the duration of COVID-19, but soon after this, it was discarded. Similarly, known antiviral compounds were also proposed and investigated to treat COVID-19. We report a pharmacophore screening using essential chemical groups derived from HCQ and known antivirals to search a natural compound chemical space. Molecular docking of HCQ under physiological condition with spike protein, 3C-like protease (3CLpro), and RNA-dependent RNA polymerase (RdRp) of SARS-CoV2 showed − 8.52 kcal/mole binding score with RdRp, while the other two proteins showed relatively weaker binding affinity. Docked complex of RdRp-HCQ is further examined using 100 ns molecular dynamic simulation. Docking and simulation study confirmed active chemical moieties of HCQ, treated as 6-point pharmacophore to screen ZINC natural compound database. Pharmacophore screening resulted in the identification of potent hit molecule [(3S,3aR,6R,6aS)-3-(5-phenylsulfanyltetrazol-1-yl)-2,3,3a,5,6,6a-hexahydrofuro[3,2-b]furan-6-yl]N-naphthalen-ylcarbamate from natural compound library. Additionally, a set of antiviral compounds with similar chemical scaffolds are also used to design a separate ligand-based pharmacophore screening. Antiviral pharmacophore screening produced a potent hit 4-[(1,5-dimethyl-3-oxo-2-phenylpyrazol-4-yl)-(2-hydroxyphenyl)methyl]-1,5-dimethyl-2-phenylpyrazol-3-one containing essential moieties that showed affinity towards RdRp. Further, both these screened compounds are docked (− 8.69 and − 8.86 kcal/mol) and simulated with RdRp protein for 100 ns in explicit solvent medium. They bind at the active site of RdRp and form direct/indirect interaction with ASP618, ASP760, and ASP761 catalytic residues of the protein. Successively, their molecular mechanics Poisson Boltzmann surface area (MMPBSA) binding energies are calculated over the simulation trajectory to determine the dynamic atomistic interaction details. Overall, this study proposes two key natural chemical moieties: (a) tetrazol and (b) phenylpyrazol that can be investigated as a potential chemical group to design inhibitors against SARS-CoV2 RdRp.Graphical abstract Supplementary InformationThe online version contains supplementary material available at 10.1007/s11030-021-10358-5.

Impact of Cd(II) on the stability of human uracil DNA glycosylase enzyme; an implication of molecular dynamics trajectories on stability analysis

Uracil DNA glycosylase is a key enzyme that identifies and removes damaged bases from DNA in the base excision repair pathway. Experimentalists have identified the possibility of Cd(II) reducing the activity of human uracil DNA glycosylase (hUNG) by binding with the enzyme replacing the catalytic water molecule. The present study focus on the stability variation of the enzyme in the presence and absence of Cd(II) and confirms the reported results with the stability analysis done using molecular dynamic (MD) simulation trajectories. The CavityPlus web server identified seven cavities for the free enzyme as possible binding sites and a cavity containing the active site of the enzyme as the best binding cavity for a ligand. Based on the CavityPlus results and the previously reported work, a free hUNG system and two systems of the enzyme with Cd(II); one with Cd(II) replacing the catalytic water molecule in the active site of the enzyme and the other replacing a non-catalytic water molecule in the active site were generated for the simulation. The simulation trajectories were used for the structural stability analysis of the enzyme in all three systems. The binding free energy of the Cd(II) with the enzyme was calculated using molecular mechanics Poisson Boltzmann surface area method. The results showed that the enzyme achieves comparatively high stability with the removal of catalytic water of the enzyme by Cd(II). Therefore, this supports the previously reported idea that Cd(II) replaces catalytic water molecules and affects enzyme activity.

Recent Developments in Free Energy Calculations for Drug Discovery.

The grand challenge in structure-based drug design is achieving accurate prediction of binding free energies. Molecular dynamics (MD) simulations enable modeling of conformational changes critical to the binding process, leading to calculation of thermodynamic quantities involved in estimation of binding affinities. With recent advancements in computing capability and predictive accuracy, MD based virtual screening has progressed from the domain of theoretical attempts to real application in drug development. Approaches including the Molecular Mechanics Poisson Boltzmann Surface Area (MM-PBSA), Linear Interaction Energy (LIE), and alchemical methods have been broadly applied to model molecular recognition for drug discovery and lead optimization. Here we review the varied methodology of these approaches, developments enhancing simulation efficiency and reliability, remaining challenges hindering predictive performance, and applications to problems in the fields of medicine and biochemistry.

In silico binding profile characterization of SARS-CoV-2 spike protein and its mutants bound to human ACE2 receptor.

Severe acute respiratory syndrome coronavirus (SARS-CoV-2), a novel coronavirus, has brought an unprecedented pandemic to the world and affected over 64 million people. The virus infects human using its spike glycoprotein mediated by a crucial area, receptor-binding domain (RBD), to bind to the human ACE2 (hACE2) receptor. Mutations on RBD have been observed in different countries and classified into nine types: A435S, D364Y, G476S, N354D/D364Y, R408I, V341I, V367F, V483A and W436R. Employing molecular dynamics (MD) simulation, we investigated dynamics and structures of the complexes of the prototype and mutant types of SARS-CoV-2 spike RBDs and hACE2. We then probed binding free energies of the prototype and mutant types of RBD with hACE2 protein by using an end-point molecular mechanics Poisson Boltzmann surface area (MM-PBSA) method. According to the result of MM-PBSA binding free energy calculations, we found that V367F and N354D/D364Y mutant types showed enhanced binding affinities with hACE2 compared to the prototype. Our computational protocols were validated by the successful prediction of relative binding free energies between prototype and three mutants: N354D/D364Y, V367F and W436R. Thus, this study provides a reliable computational protocol to fast assess the existing and emerging RBD mutations. More importantly, the binding hotspots identified by using the molecular mechanics generalized Born surface area (MM-GBSA) free energy decomposition approach can guide the rational design of small molecule drugs or vaccines free of drug resistance, to interfere with or eradicate spike-hACE2 binding.

In silico drug designing for COVID-19: an approach of high-throughput virtual screening, molecular, and essential dynamics simulations

Severe acute respiratory syndrome-coronavirus2 (SARS-CoV2), a new coronavirus has emerged in Wuhan city of China, December 2019 causing pneumonia named Coronavirus disease-19 (COVID-19), which has spread to the entire world. By January 2021, number of confirmed cumulative cases crossed ∼104 million worldwide. Till date, no effective treatment or drug is available for this virus. Availability of X-ray structures of SARS-CoV2 main protease (Mpro) provides the potential opportunity for structure-based drug designing. Here, we have made an attempt to do computational drug design by targeting main protease of SARS-CoV2. High-throughput virtual screening of million molecules and natural compounds databases were performed followed by docking. After that, the protein-ligand complexes were optimized and rescoring of binding energies were accomplished through molecular dynamics simulation and Molecular mechanics Poisson Boltzmann surface area approaches, respectively. In addition, conformational effect of various ligands on protein was also examined through essential dynamics simulation. Three compounds namely ZINC14732869, ZINC19774413, and ZINC19774479 were finally filtered that displayed better binding affinities than N3 (known) inhibitor and formed conformationally stable complexes. Hence, the current study features the potential novel inhibitors against main protease of SARS-CoV2 which might provide an effective therapeutic strategy against COVID-19. Communicated by Ramaswamy H. Sarma

A molecular docking study of potential inhibitors and repurposed drugs against SARS-CoV-2 main protease enzyme

COVID-19 has affected millions of people. Although many drugs are in use to combat disease, there is not any sufficient treatment yet. Having critical role in propagation of the novel coronavirus (SARS-CoV-2) works Main Protease up into a significant drug target. We have performed a molecular docking study to define possible inhibitor candidates against SARS-CoV-2 Main Protease enzyme. Besides docking Remdesivir, Ribavirin, Chloroquine and 28 other antiviral inhibitors (totally 31 inhibitors) to Main Protease enzyme, we have also performed a molecular docking study of 2177 ligands, which are used against Main Protease for the first time by using molecular docking program Autodock4. All ligands were successfully docked into Main Protease enzyme binding site. Among all ligands, EY16 coded ligand which previously used as EBNA1-DNA binding blocker candidate showed the best score for Main Protease with a binding free energy of −10.83 ​kcal/mol which was also lower than re-docking score of N3 ligand (−10.72 ​kcal/mol) contained in crystal structure of Main Protease. After analyzing the docking modes and docking scores we have found that our ligands have better binding free energy values than the inhibitors in use of treatment. We believe that further studies such as molecular dynamics or Molecular Mechanic Poisson Boltzmann Surface Area studies can make contribution that is more exhaustive to the docking results.