Molecular docking and molecular dynamics studies on β-lactamases and penicillin binding proteins

Discovery of putative inhibitors of human Pkd1 enzyme: Molecular docking, dynamics and simulation, QSAR, and MM/GBSA

Polycystic kidney disease is the most prevalent hereditary kidney disease globally and is mainly linked to the overexpression of a gene called PKD1. To date, there is no effective treatment available for polycystic kidney disease, and the practicing treatments only provide symptomatic relief. Discovery of the compounds targeting the PKD1 gene by inhibiting its expression under the disease condition could be crucial for effective drug development. In this study, a molecular docking and molecular dynamic simulation, QSAR, and MM/GBSA-based approaches were used to determine the putative inhibitors of the Pkd1 enzyme from a library of 1379 compounds. Initially, fourteen compounds were selected based on their binding affinities with the Pkd1 enzyme using MOE and AutoDock tools. The selected drugs were further investigated to explore their properties as drug candidates and the stability of their complex formation with the Pkd1 enzyme. Based on the physicochemical and ADMET (Absorption, Distribution, Metabolism, Excretion, and Toxicity) properties, and toxicity profiling, two compounds including olsalazine and diosmetin were selected for the downstream analysis as they demonstrated the best drug-likeness properties and highest binding affinity with Pkd1 in the docking experiment. Molecular dynamic simulation using Gromacs further confirmed the stability of olsalazine and diosmetin complexes with Pkd1 and establishing interaction through strong bonding with specific residues of protein. High biological activity and binding free energies of two complexes calculated using 3D QSAR and Schrodinger module, respectively further validated our results. Therefore, the molecular docking and dynamics simulation-based in-silico approach used in this study revealed olsalazine and diosmetin as potential drug candidates to combat polycystic kidney disease by targeting Pkd1 enzyme.

Dynamozones are the most obvious sign of the evolution of conformational dynamics in HIV-1 protease

Proteins are not static but are flexible molecules that can adopt many different conformations. The HIV-1 protease is an important target for the development of therapies to treat AIDS, due to its critical role in the viral life cycle. We investigated several dynamics studies on the HIV-1 protease families to illustrate the significance of examining the dynamic behaviors and molecular motions for an entire understanding of their dynamics-structure-function relationships. Using computer simulations and principal component analysis approaches, the dynamics data obtained revealed that: (i) The flap regions are the most obvious sign of the evolution of conformational dynamics in HIV-1 protease; (ii) There are dynamic structural regions in some proteins that contribute to the biological function and allostery of proteins via appropriate flexibility. These regions are a clear sign of the evolution of conformational dynamics of proteins, which we call dynamozones. The flap regions are one of the most important dynamozones members that are critical for HIV-1 protease function. Due to the existence of other members of dynamozones in different proteins, we propose to consider dynamozones as a footprint of the evolution of the conformational dynamics of proteins.

Ultrasonic-Assisted Synthesis of Heterocyclic Curcumin Analogs as Antidiabetic, Antibacterial, and Antioxidant Agents Combined with in vitro and in silico Studies

Background

Heterocyclic analogs of curcumin have a wide range of therapeutic potential and the ability to control the activity of a variety of metabolic enzymes.

Methods

1H-NMR and 13C-NMR spectroscopic techniques were used to determine the structures of synthesized compounds. The agar disc diffusion method and α-amylase inhibition assay were used to examine the antibacterial and anti-diabetic potential of the compounds against α-amylase enzyme inhibitory activity, respectively. DPPH-free radical scavenging and lipid peroxidation inhibition assays were used to assess the in vitro antioxidant potential.

Results and Discussion

In this work, nine heterocyclic analogs derived from curcumin precursors under ultrasonic irradiation were synthesized in excellent yields (81.4-93.7%) with improved reaction time. Results of antibacterial activities revealed that compounds 8, and 11 displayed mean inhibition zone of 13.00±0.57, and 19.66±00 mm, respectively, compared to amoxicillin (12.87±1.41 mm) at 500 μg/mL against E. coli, while compounds 8, 11 and 16 displayed mean inhibition zone of 17.67±0.57, 14.33±0.57 and 23.33±00 mm, respectively, compared to amoxicillin (13.75±1.83 mm) at 500 μg/mL against P. aeruginosa. Compound 11 displayed a mean inhibition zone of 11.33±0.57 mm compared to amoxicillin (10.75±1.83 mm) at 500 μg/mL against S. aureus. Compound 11 displayed higher binding affinities of -7.5 and -8.3 Kcal/mol with penicillin-binding proteins (PBPs) and β-lactamases producing bacterial strains, compared to amoxicillin (-7.2 and -7.9 Kcal/mol, respectively), these results are in good agreement with the in vitro antibacterial activities. In vitro antidiabetic potential on α-amylase enzyme revealed that compounds 11 (IC50=7.59 µg/mL) and 16 (IC50=4.08 µg/mL) have higher inhibitory activities than acarbose (IC50=8.0 µg/mL). Compound 8 showed promising antioxidant inhibition efficacy of DPPH (IC50 = 2.44 g/mL) compared to ascorbic acid (IC50=1.24 g/mL), while compound 16 revealed 89.9±20.42% inhibition of peroxide generation showing its potential in reducing the development of lipid peroxides. In silico molecular docking analysis, results are in good agreement with in vitro biological activity. In silico ADMET profiles suggested the adequate oral drug-likeness potential of the compounds without adverse effects.

Conclusion

According to our findings, both biological activities and in silico computational studies results demonstrated that compounds 8, 11, and 16 are promising α-amylase inhibitors and antibacterial agents against E. coli, P. aeruginosa, and S. aureus, whereas compound 8 was found to be a promising antioxidant agent.

Discovery of new thieno[2,3-d]pyrimidines as EGFR tyrosine kinase inhibitors for cancer treatment

Background: EGFR has been considered a vital molecular target in cancer management. Aim: The discovery of new thieno[2,3-d]pyrimidine derivatives as EGFR tyrosine kinase inhibitors. Methods: Nine derivatives were designed, synthesized and subjected to in vitro and in silico studies. Results: Compound 7a significantly inhibited the growth of HepG2 and PC3 cells for both EGFR wild-type and EGFRT790M. Compound 7a caused a significant apoptotic effect, arresting HepG2 cells' growth in the S and G2/M phases. Docking and molecular dynamics simulation studies confirmed the correct and stable binding modes of the synthesized compounds against the active sites. Conclusion: Compound 7a is a promising dual EGFR inhibitor for cancer treatment.

Discovery of new quinoline and isatine derivatives as potential VEGFR-2 inhibitors: design, synthesis, antiproliferative, docking and MD simulation studies

A new set of quinoline and isatine derivatives were synthesized as antiangiogenic VEGFR-2 inhibitors. On a biological level, the in vitro ability of the obtained candidates to inhibit VEGFR-2 was found to be strong with IC50 values in the range of 76.64-175.50 nM. To investigate the cytotoxicity and safety, all compounds were tested against a panel of four cancer cell lines (A549, Caco2, HepG2 and MDA) as well as two normal cell lines (Vero and WI-38). Interestingly, compound 12 exhibited noticeable cytotoxicity against A549, Caco2 and MDA with IC50 values of 5.40, 0.58 and 0.94 µM, respectively. These results were better and comparable to that of doxorubicin (0.70, 0.82 and 0.90 µM, respectively) with more than three folds higher selectivity index against the Caco2 cell lines. Compound 9 prevented the healing of the cancer cells at a low concentration. Also, the compound's potential to induce programmed cell death in Caco-2 was proved through the significant down regulating of the expression of Bcl2, Bcl-xl and Survivin in addition to the slight upregulation of the TGF-β gene. The cell cycle analysis indicated that compound 9 arrested the Caco-2 cells in the G2/M phase. Interestingly, the molecular docking studies against VEGFR-2 revealed the correct binding of the targeted compounds similar to sorafenib. Furthermore, MD experiments validated the binding of compound 12 with VEGFR-2 over 100 ns, as well as MM-PBSA analysis that confirmed the precise binding with optimum energy. Finally, ADMET analysis showed the general drug-likeness and confirmed the safety of the tested compounds.Communicated by Ramaswamy H. Sarma.

In silico drug screen reveals potential competitive MTHFR inhibitors for clinical repurposing

MTHFR (Methylenetetrahydrofolate reductase) is a pivotal enzyme involved in one-carbon metabolism, which is critical for the proliferation of cancer cells. In line with this, published literature showed that MTHFR knockdown caused impaired growth of multiple types of cancer cells. Moreover, higher MTHFR expression levels were linked to shorter overall survival in hepatocellular carcinoma, adrenocortical carcinoma, and low-grade glioma, bringing the need to design MTHFR inhibitors as a possible treatment option. No competitive inhibitors of MTHFR have been reported as of today. This study aimed to identify potential competitive MTHFR inhibitor candidates using an in silico drug screen. A total of 30470 molecules containing biogenic compounds, FDA-approved drugs, and those in clinical trials were screened against the catalytic pocket of MTHFR in the presence and absence of cofactors. Binding energy and ADMET analysis revealed that Vilanterol (β2-adrenergic agonist), Selexipag (prostacyclin receptor agonist), and Ramipril Diketopiperazine (ACE inhibitor) are potential competitive inhibitors of MTHFR. Molecular dynamics analyses and MM-PBSA calculations with these compounds particularly revealed the amino acids between 285-290 for ligand binding and highlighted Vilanterol as the strongest candidate for MTHFR inhibition. Our results could guide the development of novel MTHFR inhibitor compounds, which could be inspired by the drugs brought into the spotlight here. More importantly, these potential candidates could be quhickly tested as a repurposing strategy in pre-clinical and clinical studies of the cancers mentioned above.Communicated by Ramaswamy H. Sarma.

New quinoline and isatin derivatives as apoptotic VEGFR-2 inhibitors: design, synthesis, anti-proliferative activity, docking, ADMET, toxicity, and MD simulation studies

New quinoline and isatin derivatives having the main characteristics of VEGFR-2 inhibitors was synthesised. The antiproliferative effects of these compounds were estimated against A549, Caco-2, HepG2, and MDA-MB-231. Compounds 13 and 14 showed comparable activities with doxorubicin against the Caco-2 cells. These compounds strongly inhibited VEGFR-2 kinase activity. The cytotoxic activities were evaluated against Vero cells. Compound 7 showed the highest value of safety and selectivity. Cell migration assay displayed the ability of compound 7 to prevent healing and migration abilities in the cancer cells. Furthermore, compound 7 induced apoptosis in Caco-2 through the expressive down-regulation of the apoptotic genes, Bcl2, Bcl-xl, and Survivin, and the upregulation of the TGF gene. Molecular docking against VEGFR-2 emerged the interactions of the synthesised compounds in a similar way to sorafenib. Additionally, seven molecular dynamics simulations studies were applied and confirmed the stability of compound 13 in the active pocket of VEGFR-2 over 100 ns.

The Discovery of Potential SARS-CoV-2 Natural Inhibitors among 4924 African Metabolites Targeting the Papain-like Protease: A Multi-Phase In Silico Approach

Four compounds, hippacine, 4,2′-dihydroxy-4′-methoxychalcone, 2′,5′-dihydroxy-4-methoxychalcone, and wighteone, were selected from 4924 African natural metabolites as potential inhibitors against SARS-CoV-2 papain-like protease (PLpro, PDB ID: 3E9S). A multi-phased in silico approach was employed to select the most similar metabolites to the co-crystallized ligand (TTT) of the PLpro through molecular fingerprints and structural similarity studies. Followingly, to examine the binding of the selected metabolites with the PLpro (molecular docking. Further, to confirm this binding through molecular dynamics simulations. Finally, in silico ADMET and toxicity studies were carried out to prefer the most convenient compounds and their drug-likeness. The obtained results could be a weapon in the battle against COVID-19 via more in vitro and in vivo studies.

Discovery of New VEGFR-2 Inhibitors: Design, Synthesis, Anti-Proliferative Evaluation, Docking, and MD Simulation Studies

Four new nicotinamide-based derivatives were designed as antiangiogenic VEGFR-2 inhibitors. The congeners were synthesized possessing the pharmacophoric essential features to bind correctly with the VEGFR-2 active pocket. All members were evaluated for their cytotoxic and VEGFR-2 inhibitory potentialities. Compound 6 was the most potent showingIC₅₀ values of 9.3 ± 0.02 and 7.8 ± 0.025 µM against HCT-116 and HepG-2 cells, respectively, and IC₅₀ of 60.83 nM regarding VEGFR-2 enzyme inhibition. Compound 6 arrested the growth of HCT-116 cells at the pre-G1 and G2-M phases. Further, it induced both early and late apoptosis. Additionally, compound 6 caused a significant decrease in TNF-α and IL6 by 66.42% and 57.34%, respectively. The considered compounds had similar docking performances to that of sorafenib against the VEGFR-2 (PDB ID: 2OH4). The correct binding of compound 6 with VEGFR-2 was validated using MD simulations, and MM-GPSA calculations.

Structure-Based Virtual Screening, Docking, ADMET, Molecular Dynamics, and MM-PBSA Calculations for the Discovery of Potential Natural SARS-CoV-2 Helicase Inhibitors from the Traditional Chinese Medicine

Continuing our antecedent work against COVID-19, a set of 5956 compounds of traditional Chinese medicine have been virtually screened for their potential against SARS-CoV-2 helicase (PDB ID: 5RMM). Initially, a fingerprint study with VXG, the ligand of the target enzyme, disclosed the similarity of 187 compounds. Then, a molecular similarity study declared the most similar 40 compounds. Subsequently, molecular docking studies were carried out to examine the binding modes and energies. Then, the most appropriate 26 compounds were subjected to in silico ADMET and toxicity studies to select the most convenient inhibitors to be: (1R,2S)-ephedrine (57), (1R,2S)-norephedrine (59), 2-(4-(pyrrolidin-1-yl)phenyl)acetic acid (84), 1-phenylpropane-1,2-dione (195), 2-methoxycinnamic acid (246), 2-methoxybenzoic acid (364), (R)-2-((R)-5-oxopyrrolidin-3-yl)-2-phenylacetic acid (405), (Z)-6-(3-hydroxy-4-methoxystyryl)-4-methoxy-2H-pyran-2-one (533), 8-chloro-2-(2-phenylethyl)-5,6,7-trihydroxy-5,6,7,8-tetrahydrochromone (637), 3-((1R,2S)-2-(dimethylamino)-1-hydroxypropyl)phenol (818), (R)-2-ethyl-4-(1-hydroxy-2-(methylamino)ethyl)phenol (5159), and (R)-2-((1S,2S,5S)-2-benzyl-5-hydroxy-4-methylcyclohex-3-en-1-yl)propane-1,2-diol (5168). Among the selected 12 compounds, the metabolites, compound 533 showed the best docking scores. Interestingly, the MD simulation studies for compound 533, the one with the highest docking score, over 100 ns showed its correct binding to SARS-CoV-2 helicase with low energy and optimum dynamics. Finally, MM-PBSA studies showed that 533 bonded favorably to SARS-CoV-2 helicase with a free energy value of −83 kJ/mol. Further, the free energy decomposition study determined the essential amino acid residues that contributed favorably to the binding process. The obtained results give a huge hope to find a cure for COVID-19 through further in vitro and in vivo studies for the selected compounds.

Protein-Ligand Docking in the Machine-Learning Era

Molecular docking plays a significant role in early-stage drug discovery, from structure-based virtual screening (VS) to hit-to-lead optimization, and its capability and predictive power is critically dependent on the protein-ligand scoring function. In this review, we give a broad overview of recent scoring function development, as well as the docking-based applications in drug discovery. We outline the strategies and resources available for structure-based VS and discuss the assessment and development of classical and machine learning protein-ligand scoring functions. In particular, we highlight the recent progress of machine learning scoring function ranging from descriptor-based models to deep learning approaches. We also discuss the general workflow and docking protocols of structure-based VS, such as structure preparation, binding site detection, docking strategies, and post-docking filter/re-scoring, as well as a case study on the large-scale docking-based VS test on the LIT-PCBA data set.

The Inhibitory Potential of 2'-dihalo Ribonucleotides against HCV: Molecular Docking, Molecular Simulations, MM-BPSA, and DFT Studies

Sofosbuvir is the first approved direct-acting antiviral (DAA) agent that inhibits the HCV NS5B polymerase, resulting in chain termination. The molecular models of the 2′-dihalo ribonucleotides used were based on experimental biological studies of HCV polymerase inhibitors. They were modeled within HCV GT1a and GT1b to understand the structure–activity relationship (SAR) and the binding interaction of the halogen atoms at the active site of NS5B polymerase using different computational approaches. The outputs of the molecular docking studies indicated the correct binding mode of the tested compounds against the active sites in target receptors, exhibiting good binding free energies. Interestingly, the change in the substitution at the ribose sugar was found to produce a mild effect on the binding mode. In detail, increasing the hydrophobicity of the substituted moieties resulted in a better binding affinity. Furthermore, in silico ADMET investigation implied the general drug likeness of the examined derivatives. Specifically, good oral absorptions, no BBB penetration, and no CYP4502D6 inhibitions were expected. Likely, the in silico toxicity studies against several animal models showed no carcinogenicity and high predicted TD50 values. The DFT studies exhibited a bioisosteric effect between the substituents at the 2′-position and the possible steric clash between 2′-substituted nucleoside analogs and the active site in the target enzyme. Finally, compound 6 was subjected to several molecular dynamics (MD) simulations and MM-PBSA studies to examine the protein-ligand dynamic and energetic stability.

Isolation and In Silico SARS-CoV-2 Main Protease Inhibition Potential of Jusan Coumarin, a New Dicoumarin from Artemisia glauca

A new dicoumarin, jusan coumarin, (1), has been isolated from Artemisia glauca aerial parts. The chemical structure of jusan coumarin was estimated, by 1D, 2D NMR as well as HR-Ms spectroscopic methods, to be 7-hydroxy-6-methoxy-3-[(2-oxo-2H-chromen-6-yl)oxy]-2H-chromen-2-one. As the first time to be introduced in nature, its potential against SARS-CoV-2 has been estimated using various in silico methods. Molecular similarity and fingerprints experiments have been utilized for 1 against nine co-crystallized ligands of COVID-19 vital proteins. The results declared a great similarity between Jusan Coumarin and X77, the ligand of COVID-19 main protease (PDB ID: 6W63), Mpro. To authenticate the obtained outputs, a DFT experiment was achieved to confirm the similarity of X77 and 1. Consequently, 1 was docked against Mpro. The results clarified that 1 bonded in a correct way inside Mpro active site, with a binding energy of -18.45 kcal/mol. Furthermore, the ADMET and toxicity profiles of 1 were evaluated and showed the safety of 1 and its likeness to be a drug. Finally, to confirm the binding and understand the thermodynamic characters between 1 and Mpro, several molecular dynamics (MD) simulations studies have been administered. Additionally, the known coumarin derivative, 7-isopentenyloxycoumarin (2), has been isolated as well as β-sitosterol (3).

Truncated human angiotensin converting enzyme 2; a potential inhibitor of SARS-CoV-2 spike glycoprotein and potent COVID-19 therapeutic agent

The current pandemic of Covid-19 caused by SARS-CoV-2 is continued to spread globally and no potential drug or vaccine against it is available. Spike (S) glycoprotein is the structural protein of SARS-CoV-2 located on the envelope surface, involve in interaction with angiotensin converting enzyme 2 (ACE2), a cell surface receptor, followed by entry into the host cell. Thereby, blocking the S glycoprotein through potential inhibitor may interfere its interaction with ACE2 and impede its entry into the host cell. Here, we present a truncated version of human ACE2 (tACE2), comprising the N terminus region of the intact ACE2 from amino acid position 21-119, involved in binding with receptor binding domain (RBD) of SARS-CoV-2. We analyzed the in-silico potential of tACE2 to compete with intact ACE2 for binding with RBD. The protein-protein docking and molecular dynamic simulation showed that tACE2 has higher binding affinity for RBD and form more stabilized complex with RBD than the intact ACE2. Furthermore, prediction of tACE2 soluble expression in E. coli makes it a suitable candidate to be targeted for Covid-19 therapeutics. This is the first MD simulation based findings to provide a high affinity protein inhibitor for SARS-CoV-2 S glycoprotein, an important target for drug designing against this unprecedented challenge.Communicated by Ramaswamy H. Sarma.

Prevention of SARS-CoV-2 Proliferation with a Novel and Potent Main Protease Inhibitor by Docking, ADMET, MM-PBSA, and Molecular Dynamics Simulation

COVID-19 is the last disease caused by SARS-CoV-2 associated with a severe immune response and lung damage. The main protease (Mpro) has a vital role in SARS-CoV-2 proliferation. Moreover, humans lack homologous Mpro, which makes the Mpro a suitable drug target for the development of SARS-CoV-2 drugs. The purchasable L5000 library (Selleckchem Inc) includes 99,040 compounds that were used for virtual screening. After molecular docking and ADME studies, we selected a compound (WAY-604395) with a potent binding affinity to the Mpro active site and acceptable ADME properties compared to the reference drug (nelfinavir). Molecular dynamics (MD) simulation outcomes have proved that the Mpro-WAY604395 complex possesses a considerable value of flexibility, stability, compactness and binding energy. Our Molecular Mechanics Poisson–Boltzmann Surface Area (MM-PBSA) calculation demonstrates that WAY-604395 is more potent ([Formula: see text]272.19[Formula: see text]kcal mol[Formula: see text]) in comparison with nelfinavir ([Formula: see text]173.39[Formula: see text]kcal[Formula: see text]mol[Formula: see text]) against SARS-CoV-2 Mpro. In conclusion, we suggest that WAY-604395 has the potential for the treatment of SARS-CoV-2 by inhibition of the Mpro.

Transdermal delivery of a somatostatin receptor type 2 antagonist using microneedle patch technology for hypoglycemia prevention

Hypoglycemia is a serious and potentially fatal complication experienced by people with insulin-dependent diabetes. The complication is usually caused by insulin overdose, skipping meals, and/or excessive physical activities. In type 1 diabetes (T1D), on top of impaired pancreatic α-cells, excessive levels of somatostatin from δ-cells further inhibit glucagon secretion to counteract overdosed insulin. Herein, we aimed to develop a microneedle (MN) patch for transdermal delivery of a peptide (PRL-2903) that antagonizes somatostatin receptor type 2 (SSTR2) in α-cells. First, we investigated the efficacy of subcutaneously administered PRL-2903 and identified the optimal dose (i.e., the minimum effective dose) and treatment scheduling (i.e., the best administration time for hypoglycemia prevention) in a T1D rat model. We then designed an MN patch using a hyaluronic acid (HA)-based polymer. The possible effect of the polymer on stabilizing the native structure of PRL-2903 was studied by molecular dynamics (MD) simulations. The results showed that the HA-based polymer could stabilize the PRL-2903 structure by restricting water molecules, promoting intra-molecular H-bonding, and constraining torsional angles of important bonds. In vivo studies with an overdose insulin challenge revealed that the PRL-2903-loaded MN patch effectively increased the plasma glucagon level, restored the counter-regulation of blood glucose concentration, and prevented hypoglycemia. The proposed MN patch is the first demonstration of a transdermal microneedle patch designed to deliver an SSTR2 antagonist for the prevention of hypoglycemia. This counter-regulatory peptide delivery system may be applied alongside with insulin delivery systems to provide a more effective and safer treatment for people with insulin-dependent diabetes.

Designing a natural inhibitor against human kynurenine aminotransferase type II and a comparison with PF-04859989: a computational effort against schizophrenia

The Kynurenine aminotransferase II (KATII) enzyme has an essential role in L-kynurenine transmission to kynurenic acid (KYNA). High concentration of kynurenic acid associates with schizophrenia and some neurocognitive disorders. Decreasing KYNA production via inhibiting KATII would be an effective method for treating and understanding the related central nervous system (CNS) diseases. This study aimed to discover a potent inhibitor against human KATII (hKATII) in comparison with PF-04859989. We utilized the computational methods of molecular dynamics, virtual screening, docking, and binding free-energy calculations. Initially, the 58722 compounds from three drug libraries, including IBS library, DrugBank library, and Analyticon library, were obtained. At the next stage, these sets of compounds were screened by AutoDock Vina software, and a potent inhibitor (ZINC35466084) was selected. Following the screening, molecular dynamics simulations for both ZINC35466084 and PF-04859989 were performed by GROMACS software. MM-PBSA analysis showed that the amount of binding free energy for ZINC35466084 (-61.26 KJ mol^-1) is more potent than PF-04859989 (-43.14 KJ mol^-1). Furthermore, the ADME analysis results revealed that the pharmacokinetic parameters of ZINC35466084 are acceptable for human use. Eventually, our data demonstrated that ZINC35466084 is suitable for hKATII inhibition, and it is an appropriate candidate for further studies in the laboratory. Communicated by Ramaswamy H. Sarma.

Theoretical Evaluation of Novel Thermolysin Inhibitors from Bacillus thermoproteolyticus. Possible Antibacterial Agents

The search for new antibacterial agents that could decrease bacterial resistance is a subject in continuous development. Gram-negative and Gram-positive bacteria possess a group of metalloproteins belonging to the MEROPS peptidase (M4) family, which is the main virulence factor of these bacteria. In this work, we used the previous results of a computational biochemistry protocol of a series of ligands designed in silico using thermolysin as a model for the search of antihypertensive agents. Here, thermolysin from Bacillus thermoproteolyticus, a metalloprotein of the M4 family, was used to determine the most promising candidate as an antibacterial agent. Our results from docking, molecular dynamics simulation, molecular mechanics Poisson-Boltzmann (MM-PBSA) method, ligand efficiency, and ADME-Tox properties (Absorption, Distribution, Metabolism, Excretion, and Toxicity) indicate that the designed ligands were adequately oriented in the thermolysin active site. The Lig783, Lig2177, and Lig3444 compounds showed the best dynamic behavior; however, from the ADME-Tox calculated properties, Lig783 was selected as the unique antibacterial agent candidate amongst the designed ligands.

Discovery of Potential Chemical Probe as Inhibitors of CXCL12 Using Ligand-Based Virtual Screening and Molecular Dynamic Simulation

CXCL12 are small pro-inflammatory chemo-attractant cytokines that bind to a specific receptor CXCR4 with a role in angiogenesis, tumor progression, metastasis, and cell survival. Globally, cancer metastasis is a major cause of morbidity and mortality. In this study, we targeted CXCL12 rather than the chemokine receptor (CXCR4) because most of the drugs failed in clinical trials due to unmanageable toxicities. Until now, no FDA approved medication has been available against CXCL12. Therefore, we aimed to find new inhibitors for CXCL12 through virtual screening followed by molecular dynamics simulation. For virtual screening, active compounds against CXCL12 were taken as potent inhibitors and utilized in the generation of a pharmacophore model, followed by validation against different datasets. Ligand based virtual screening was performed on the ChEMBL and in-house databases, which resulted in successive elimination through the steps of pharmacophore-based and score-based screenings, and finally, sixteen compounds of various interactions with significant crucial amino acid residues were selected as virtual hits. Furthermore, the binding mode of these compounds were refined through molecular dynamic simulations. Moreover, the stability of protein complexes, Root Mean Square Deviation (RMSD), Root Mean Square Fluctuation (RMSF), and radius of gyration were analyzed, which led to the identification of three potent inhibitors of CXCL12 that may be pursued in the drug discovery process against cancer metastasis.

Truncated human angiotensin converting enzyme 2; a potential inhibitor of SARS-CoV-2 spike glycoprotein and potent COVID-19 therapeutic agent

The current pandemic of Covid-19 caused by SARS-CoV-2 is continued to spread globally and no potential drug or vaccine against it is available. Spike (S) glycoprotein is the structural protein of SARS-CoV-2 located on the envelope surface, involve in interaction with angiotensin converting enzyme 2 (ACE2), a cell surface receptor, followed by entry into the host cell. Thereby, blocking the S glycoprotein through potential inhibitor may interfere its interaction with ACE2 and impede its entry into the host cell. Here, we present a truncated version of human ACE2 (tACE2), comprising the N terminus region of the intact ACE2 from amino acid position 21-119, involved in binding with receptor binding domain (RBD) of SARS-CoV-2. We analyzed the in-silico potential of tACE2 to compete with intact ACE2 for binding with RBD. The protein-protein docking and molecular dynamic simulation showed that tACE2 has higher binding affinity for RBD and form more stabilized complex with RBD than the intact ACE2. Furthermore, prediction of tACE2 soluble expression in E. coli makes it a suitable candidate to be targeted for Covid-19 therapeutics. This is the first MD simulation based findings to provide a high affinity protein inhibitor for SARS-CoV-2 S glycoprotein, an important target for drug designing against this unprecedented challenge.Communicated by Ramaswamy H. Sarma.

Construction of an Electrochemical Receptor Sensor Based on Graphene/Thionine for the Sensitive Determination of β-Lactam Antibiotics Content in Milk

In antibiotics, β-lactam is one kind of major concern acknowledged as an unavoidable contaminant in milk. Thus, a facile and sensitive method is essential for rapid β-lactam antibiotics detection. In our work, a specific electrochemical receptor sensor based on the graphene/thionine (GO/TH) composite was established. The mechanism of the electrochemical receptor sensor was a direct competitive inhibition of the binding of horseradish peroxidase-labeled ampicillin (HRP-AMP) to the mutant BlaR-CTD protein by free β-lactam antibiotics. Then, horseradish peroxidase (HRP) catalyzed the hydrolysis of the substrate hydrogen peroxide (H₂O₂), which produced an electrochemical signal. Under optimal experimental conditions, this method could quantitatively detect cefquinome from 0.1 to 8 μg L⁻¹ and with the limit of detection (LOD) of 0.16 μg L⁻¹, much lower than the maximum residue limit (MRL) of 5 μg L⁻¹ set by the European Union. In addition, the LOD of spiked milk samples with cefalexin, cefquinoxime, cefotafur, penicillin G and ampicillin were 14.88 μg L⁻¹, 2.46 μg L⁻¹, 17.16 μg L⁻¹, 0.06 μg L⁻¹, 0.21 μg L⁻¹ and the limits of quantitation (LOQ) were 36.09 μg L⁻¹, 5.40 μg L⁻¹, 41.45 μg L⁻¹, 0.13 μg L⁻¹, 0.42 μg L⁻¹, respectively. The sensor showed a favorable recovery of 84.89–102.44%. Moreover, the electrochemical receptor sensor was successfully applied to assay β-lactam antibiotics in milk, which showed good correlation with the results obtained from liquid chromatography-tandem mass spectrometry (LC-MS/MS).

Solvent induced mononuclear and dinuclear mixed ligand Cu(ii) complex: structural diversity, supramolecular packing polymorphism and molecular docking studies

Phenolate bridged dinuclear and solvent induced mononuclear supramolecular isomers of Cu(ii) complex have been reported to explore the structural diversity and their antibacterial activity supported by molecular docking studies.

Comparative Molecular Field Analysis and Molecular Docking Studies on Quinolinone Derivatives Indicate Potential Hepatitis C Virus Inhibitors

Presently, there are no effective vaccines and anti-virals for the prevention and treatment of Hepatitis C virus infections and hence there is an urgent need to develop potent HCV inhibitors. In this study, we have carried out molecular docking, molecular dynamics and 3D-QSAR on heteroaryl 3-(1,1-dioxo-2H-(1,2,4)-benzothiadizin-3-yl)-4-hydroxy-2(1H)-quinolinone series using NS5B protein. Total of 41 quinolinone derivatives is used for molecular modeling study. The binding conformation and hydrogen bond interaction of the docked complexes were analyzed to model the inhibitors. We identified the molecule XXXV that had a higher affinity with NS5B. The molecular dynamics study confirmed the stability of the compound XXXV-NS5B complex. The developed CoMFA descriptors parameters, which were calculated using a test set of 13 compounds, were statistically significant. Our results will provide useful insights and lead to design potent anti-Hepatitis C virus molecules.

Identification of potential inhibitors for Klebsiella pneumoniae carbapenemase-3: a molecular docking and dynamics study

Klebsiella pneumoniae (K. pneumoniae) is a Gram-negative bacterium, which is a leading causal agent for nosocomial infections. Penicillin, cephalosporin and carbapenems along with the inhibitors such as tazobactam, sulbactam and clavulanic acid are prescribed for the treatment of K. pneumoniae infections. Prolonged exposure to β-lactam antibiotics leads to the development of resistance. The major reason for the β-lactam resistance in K. pneumoniae is the secretion of the enzyme K. pneumoniae carbapenemase (KPC). Secretion of KPC-2 and its variant KPC-3 by the K. pneumoniae strains causes resistance to both the substrate imipenem and the β-lactamase inhibitors. Hence, molecular docking and dynamics studies were carried out to analyze the resistance mechanism of KPC-2-imipenem and KPC-3-imipenem at the structural level. It reveals that KPC-3-imipenem has the highest c-score value of 4.03 with greater stability than the KPC-2-imipenem c-score value of 2.36. Greater the interaction between the substrate and the β-lactamase enzyme, higher the chances of hydrolysis of the substrate. Presently available β-lactamase inhibitors are also ineffective against KPC-3-expressing strains. This situation necessitates the need for development of novel and effective inhibitors for KPC-3. We have carried out the virtual screening process to identify more effective inhibitors for KPC-3, and this has resulted in ZINC48682523, ZINC50209041 and ZINC50420049 as the best binding energy compounds, having greater binding affinity and stability than KPC-3-tazobactam interactions. Our study provides a clear understanding of the mechanism of drug resistance and provides valuable inputs for the development of inhibitors against KPC-3 expressing K. pneumoniae. Communicated by Ramaswamy H. Sarma.

Identification of novel inhibitor against endonuclease subunit of Influenza pH1N1 polymerase: A combined molecular docking, molecular dynamics, MMPBSA, QMMM and ADME studies to combat influenza A viruses

The influenza H1N1 virus is the causative agent of the flu pandemic in the world. Due to the shortage of effective means of control, it is remained the serious threats to public and avian health. To battle the surge of viral outbreaks, new treatments are crucially needed. The viral RNA polymerase, which is responsible for transcription and replication of the RNA genome, is comprised of subunits PA, PB1 and PB2. PA has endonuclease activity and is a well known target for inhibitor and drug design. In the current study, we employed molecular docking, molecular dynamics (MD), MMPBSA, QMMM and ADME studies to find and propose an inhibitor among 11,873 structures against PA. Our molecular docking, MD, MMPBSA and QMMM studies showed that ZINC15340668 has ideal characteristics as a potent PA inhibitor, and can be used in experimental phase and further development. Also, ADME prediction demonstrated that all physico-chemical parameters are within the acceptable range defined for human use. Molecular mechanism based study revealed that upon inhibitor binding; the flexibility of PA backbone is increased. This observation demonstrates the plasticity of PA active site, and it should be noticed in drug design against PA Influenza A viruses. In the final phase of the study, the efficiency of our proposed hit was tested computationally against mutant drug resistant I38T_PA. Our results exhibited that the hit inhibits the I38T_PA in different manner with high potency.

Novel mutant of Escherichia coli asparaginase II to reduction of the glutaminase activity in treatment of acute lymphocytic leukemia by molecular dynamics simulations and QM-MM studies

L-Asparaginases (ASNase) belong to a family of amidohydrolases, have both asparaginase and glutaminase activity. Acute lymphocytic leukemia (ALL) is an outrageous disease worldwide. Bacterial ASNase has been used for the treatment of ALL. Glutaminase activity of enzyme causes some side effect and it is not essential for anticancer activity. The aim of this study was engineering of Escherichia coli asparaginase II to find a mutant with reduced glutaminase activity by molecular docking, molecular dynamics (MD) and QM-MM (Quantum mechanics molecular dynamics) simulations. Residues with low free energy of binding to Asn and high free binding energy to Gln were chosen for mutagenesis. Then, a mutant with higher glutaminase free binding energy was selected for further studies. Additionally, the MD simulation and QM-MM computation of wild type (WT) were employed and the selected mutated ASNase were analyzed and discussed. Our data showed that V27T is a good candidate to reduction the glutaminase activity, while has no remarkable effect on asparaginase activity of the enzyme. The simulation analysis revealed that V27T mutant is more stable than WT and mutant simulation was successful completely. QM-MM results confirmed the successfulness of our mutagenesis.

Ethyl 3-oxo-2-(2,5-dioxopyrrolidin-3-yl)butanoate Derivatives: Anthelmintic and Cytotoxic Potentials, Antimicrobial, and Docking Studies

Development of multidrug resistance (MDR) to antimicrobial, antiparasitic and chemotherapeutic agents is a global challenge for the scientific community. Despite of the emergence of MDR pathogens, the development of novel and more effective drugs is slow and scientist even speculate that we are going back the pre-antibiotic era. This work aims to study and evaluate the preliminary antibacterial, anthelmintic and cytotoxic potentials of ethyl 3-oxo-2-(2,5-dioxopyrrolidin-3-yl)butanoates. Among all of the four compounds, compound 2 has displayed remarkable potency with MIC values of 0.125, 0.083, 0.073, and 0.109 mg/ml against E. sakazakii, E. coli. S. aureus, and K. pneumonia, respectively. Compared to etoposide (LC50 9.8 μg/ml), the compounds demonstrated LC50 values from 280 to 765 μg/ml. For anthelmintic assay, three concentrations of each compound and standard drug were studied in determination of time of death of the two species. Excellent anthelmintic activity was observed by all four compounds against P. posthuma and A. galli better than standard albendazole. High GOLD fitness score data from docking analysis toward the targets represent better protein-ligand binding affinity and thus indicate a high propensity for all the active compounds to bind to the active site. The promising in-vitro antimicrobial, anthelmintic activity, and cytotoxicity data conclusively revealed that these compounds may serve as viable lead compounds for the treatment of bacterial and parasitic infections, and therefore, could help the medicinal chemists to design future chemotherapeutic agents to avoid rapid drug resistance.

Interactions of carbon nanotubes and/or graphene with manganese peroxidase during biodegradation of endocrine disruptors and triclosan

Molecular-level biodegradation processes of bisphenol A (BPA), nonylphenol (NP) and triclosan (TCS) mediated by manganese peroxidase (MnP) were investigated with and without single-walled carbon nanotube (SWCNT) and/or graphene (GRA). Although the incorporation of SWCNT, GRA or their combination (SWCNT+GRA) did not break up the complexes composed of manganese peroxidase (MnP) and these substrates, they had different effects on the native contacts between the substrates and MnP. GRA tended to decrease the overall stability of the binding between MnP and its substrates. SWCNT or SWCNT+GRA generally had a minor impact on the mean binding energy between MnP and its substrates. We detected some sensitive residues from MnP that were dramatically disturbed by the GRA, SWCNT or SWCNT+GRA. Nanomaterials changed the number and behavior of water molecules adjacent to both MnP and its substrates, which was not due to the destruction of H-bond network formed by sensitive regions and water molecules. The present results are useful for understanding the molecular basis of pollutant biodegradation affected by the nanomaterials in the environment, and are also helpful in assessing the risks of these materials to the environment.

MurB as a target in an alternative approach to tackle the Vibrio cholerae resistance using molecular docking and simulation study

Cholera is a serious threat to a large population in the under developed countries. Though oral rehydration therapy is the preferred choice of treatment, the use of antibiotics could reduce the microbial load in the case of severity. The use of antibiotics is also sought in the scenarios where there is problem with access to clean water. However, Vibrio cholera (V. cholerae) strains have developed resistance to antibiotics such as amoxicillin, ampicillin, chloramphenicol, doxycycline, erythromycin, and tetracycline. In this work, we have addressed the resistance issue by targeting MurB protein which is essential for the cell wall biosynthesis in V. cholerae. 20 Phytochemical compounds were chosen to screen the potential inhibitors against V. cholerae to avoid the complications faced by synthesis of small molecules. The molecular docking and dynamics study indicates that quercetin is the most potential and stable inhibitor of Murb.

Potential of Icariin Metabolites from Epimedium koreanum Nakai as Antidiabetic Therapeutic Agents

The therapeutic properties of Epimedium koreanum are presumed to be due to the flavonoid component icariin, which has been reported to have broad pharmacological potential and has demonstrated anti-diabetic, anti-Alzheimer’s disease, anti-tumor, and hepatoprotective activities. Considering these therapeutic properties of icariin, its deglycosylated icaritin and glycosylated flavonoids (icaeriside II, epimedin A, epimedin B, and epimedin C) were evaluated for their ability to inhibit protein tyrosine phosphatase 1B (PTP1B) and α-glucosidase. The results show that icaritin and icariside II exhibit potent inhibitory activities, with 50% inhibition concentration (IC₅₀) values of 11.59 ± 1.39 μM and 9.94 ± 0.15 μM against PTP1B and 74.42 ± 0.01 and 106.59 ± 0.44 μM against α-glucosidase, respectively. With the exceptions of icaritin and icariside II, glycosylated flavonoids did not exhibit any inhibitory effects in the two assays. Enzyme kinetics analyses revealed that icaritin and icariside II demonstrated noncompetitive-type inhibition against PTP1B, with inhibition constant (K_i) values of 11.41 and 11.66 μM, respectively. Moreover, molecular docking analysis confirmed that icaritin and icariside II both occupy the same site as allosteric ligand. Thus, the molecular docking simulation results were in close agreement with the experimental data with respect to inhibition activity. In conclusion, deglycosylated metabolites of icariin from E. koreanum might offer therapeutic potential for the treatment of type 2 diabetes mellitus.

Anti-Helicobacter pylori Activity of Isocoumarin Paepalantine: Morphological and Molecular Docking Analysis

The Helicobacterpylori bacterium is one of the main causes of chronic gastritis, peptic ulcers, and even gastric cancer. It affects an average of half of the world population. Its difficult eradication depends upon multi-drug therapy. Since its classification as a group 1 carcinogenic by International Agency for Research on Cancer (IARC), the importance of H. pylori eradication has obtained a novel meaning. There is considerable interest in alternative therapies for the eradication of H. pylori using compounds from a wide range of natural products. In the present study, we investigated the antibacterial property of the isocoumarin paepalantine against H. pylori and it exhibited significant anti-H. pylori activity at a minimum inhibitory concentration (MIC) of 128 μg/mL and at a minimum bactericidal concentration (MBC) of 256 μg/mL. The scanning electron microscopy (SEM) revealed significant morphological changes of the bacterial cell as a response to a sub-MIC of paepalantine, suggesting a penicillin-binding protein (PBP) inhibition. Computational studies were carried out in order to study binding modes for paepalantine in PBP binding sites, exploring the active and allosteric sites. The data from the present study indicates that paepalantine exhibits significant anti-H. pylori activity, most likely by inhibiting membrane protein synthesis.

The effect of β-glucan and its potential analog on the structure of Dectin-1 receptor

Dectin-1 is a recently discovered pattern-recognition receptor that plays an important role in antifungal innate immunity, which acts a specific receptor for β-glucan (BG). The present study, aimed at clarifying effect of BG and a new analog, maltotriose (MT) on Dectin-1 receptor. We implemented molecular docking of MT on Dectin-1 along with model-independent all-atom-molecular dynamics simulations. Simulations were carried out at three levels of complexity: (1) Apo-Dectin-1; (2) BG:Dectin-1; (3) MT:Dectin-1. All three system complexes were undergone stability check before showing a comparative analysis. A characteristic feature, noted for the MT:Dectin-1, is a shifting of loops (loop1 and loop2) orientation towards atoms of MT, a broad interaction suggested a robust and tight binding on comparison with BG:Dectin-1. Free energy estimation corroborated the observation, which furthermore, made a close agreement by revealing contribution of energy components of interacting residues. In addition, cluster analysis of complexes exhibit a smooth continuous transition to a new confirmation, represented by a series of clusters each having a longer lifetime. Principal component analysis revealed a broken pipe at binding site of BG:Dectin-1 during movement of atoms whereas in MT:Dectin-1 exhibited wide band and high amplitude motion of atoms in trajectory, was due to loop orientation toward MT. Observation was further shown by measuring distances and hydrogen binding calculation. Simulations of the BG:Dectin-1 and MT:Dectin-1 complex revealed first time the influence of BG and MT ligands. This study might extend the knowledge of the BG and MT interaction on Dectin-1 and proposed further potential bioassay of MT.

Synergistic antibacterial effects of herbal extracts and antibiotics on methicillin-resistant Staphylococcus aureus: A computational and experimental study

Antibiotic resistance has become a serious global concern, and the discovery of antimicrobial herbal constituents may provide valuable solutions to overcome the problem. In this study, the effects of therapies combining antibiotics and four medicinal herbs on methicillin-resistant Staphylococcus aureus (MRSA) were investigated. Specifically, the synergistic effects of Magnolia officinalis, Verbena officinalis, Momordica charantia, and Daphne genkwa in combination with oxacillin or gentamicin against methicillin-resistant (ATCC43300) and methicillin-susceptible (ATCC25923) S. aureus were examined. In vitro susceptibility and synergistic testing were performed to measure the minimum inhibitory concentration and fractional inhibitory concentration (FIC) index of the antibiotics and medicinal herbs against MRSA and methicillin-susceptible S. aureus. To identify the active constituents in producing these synergistic effects, in silico molecular docking was used to investigate the binding affinities of 139 constituents of the four herbs to the two common MRSA inhibitory targets, penicillin binding proteins 2a (PBP2a) and 4 (PBP4). The physicochemical and absorption, distribution, metabolism, and excretion properties and drug safety profiles of these compounds were also analyzed. D. genkwa extract potentiated the antibacterial effects of oxacillin against MRSA, as indicated by an FIC index value of 0.375. M. officinalis and V. officinalis produced partial synergistic effects when combined with oxacillin, whereas M. charantia was found to have no beneficial effects in inhibiting MRSA. Overall, tiliroside, pinoresinol, magnatriol B, and momorcharaside B were predicted to be PBP2a or PBP4 inhibitors with good drug-like properties. This study identifies compounds that deserve further investigation with the aim of developing therapeutic agents to modulate the effect of antibiotics on MRSA. Impact statement Antibiotic resistant is a well-known threat to global health and methicillin-resistant Staphylococcus aureus is one of the most significant ones. These resistant bacteria kill thousands of people every year and therefore a new effective antimicrobial treatment is necessary. This study identified the herbs and their associated bioactive ingredients that can potential the effects of current antibiotics. These herbs have long history of human usage in China and have well-defined monograph in the Chinese Pharmacopeia. These indicate their relatively high clinical safety and may have a quicker drug development process than that of a new novel antibiotic. Based on the results of this study, the authors will perform further in vitro and animal studies, aiming to accumulate significant data for the application of clinical trial.

Natural Inhibitors of HMG-CoA Reductase-An Insilico Approach Through Molecular Docking and Simulation Studies

Plant products have always been considered for many important metabolic disorders due to its abundant medicinal properties. Alarming adverse effects of overuse of statins has been reported for patients with dyslipidemia. This study was aimed to identify compounds with potent anti-dyslipidemic property from selected plants and analyze them for their efficiency in binding with HMG-CoA reductase, a key enzyme in lipid metabolism. The docking studies indicate rutin as the best compound that can inhibit HMG-CoA reductase as it had strong binding affinity to the enzyme. The molecular dynamics simulation studies confirmed the stability of the HMG-CoA reductase-rutin complex. RMSD, RMSF, Rg, H-bond results indicated that the HMG-CoA reductase-rutin complex is highly stable. Presently, statins are not preferred for individuals with pre-existing liver disease. Our study identified rutin as a promising lead compound which could be further developed into an anti-dyslipidemic molecule. Our results will be a good starting point for future experimental and clinical studies and if the results from such studies match international standards plant derived rutin might emerge as a good alternative to statins. J. Cell. Biochem. 118: 52-57, 2017. © 2016 Wiley Periodicals, Inc.