Targeting the cell wall of Mycobacterium tuberculosis: a molecular modeling investigation of the interaction of imipenem and meropenem with L,D-transpeptidase 2

The single crystal X-ray structure of the extracellular portion of the L,D-transpeptidase (ex-Ldt_Mt2 – residues 120–408) enzyme was recently reported. It was observed that imipenem and meropenem inhibit activity of this enzyme, responsible for generating L,D-transpeptide linkages in the peptidoglycan layer of Mycobacterium tuberculosis. Imipenem is more active and isothermal titration calorimetry experiments revealed that meropenem is subjected to an entropy penalty upon binding to the enzyme. Herein, we report a molecular modeling approach to obtain a molecular view of the inhibitor/enzyme interactions. The average binding free energies for nine commercially available inhibitors were calculated using MM/GBSA and Solvation Interaction Energy (SIE) approaches and the calculated energies corresponded well with the available experimentally observed results. The method reproduces the same order of binding energies as experimentally observed for imipenem and meropenem. We have also demonstrated that SIE is a reasonably accurate and cost-effective free energy method, which can be used to predict carbapenem affinities for this enzyme. A theoretical explanation was offered for the experimental entropy penalty observed for meropenem, creating optimism that this computational model can serve as a potential computational model for other researchers in the field.     

Interactions of selected indole derivatives with COX-2 and their in silico structure modifications towards the development of novel NSAIDs

Cyclooxygenase-2 (COX-2) is an important enzyme responsible for the formation of potent inflammatory mediators like prostaglandins, prostacyclin and thromboxane. Hence, inhibition of COX-2 is one of the best ways to control the inflammation. Non-steroidal anti-inflammatory drugs can control inflammation by inhibiting Cyclooxygenase. Selective inhibition of COX-2 is preferable over the inhibition of COX-1 because of the fewer adverse effects produced. Molecular modeling and docking of 134 selected indole compounds were done against COX-2. The pharmacophore-based in silico structural modifications of the best scored compounds were carried out in order to enhance the binding affinity and selectivity. The modification resulted in derivatives with better binding energies than that of known COX-2 inhibitors. The four best derivatives in terms of the binding energies were selected and their binding stabilities were studied by molecular dynamics simulation methods.     

A profound computational study to prioritize the natural compound inhibitors against the P. falciparum orotidine-5-monophosphate decarboxylase enzyme

Abstract

In the current contribution, a multicomplex-based pharmacophore modeling approach was employed on the structural proteome of Plasmodium falciparum orotidine-5-monophosphate decarboxylase enzyme (PfOMPDC). Among the constructed pharmacophore models, the representative hypotheses were selected as the primary filter to screen the molecules with the complementary features responsible for showing inhibition. Thereafter, auxiliary evaluations were performed on the screened candidates via drug-likeness and molecular docking studies. Subsequently, the stability of the docked protein-ligand complexes was scrutinized by employing molecular dynamics simulations and molecular mechanics-Poisson Boltzmann surface area based free binding energy calculations. The stability the docked candidates was compared with the highly active crystallized inhibitor (3S9Y-FNU) to seek more potential candidates. All the docked molecules displayed stable dynamic behavior and high binding free energy in comparison to 3S9Y-FNU. The employed workflow resulted in the retrieval of five drug-like candidates with diverse scaffolds that may show inhibitory activity against PfOMPDC and could be further used as the novel scaffold to develop novel antimalarials.

Communicated by Ramaswamy H. Sarma     

The mechanism of solvent effect on the positional selectivity of Candida antarctica lipase B during 1,3-diolein synthesis by esterification

We investigated the influence of solvent on the positional selectivity of Novozym 435 which was the immobilized Candida antarctica lipase B (CALB) during the esterification of oleic acid with glycerol for 1,3-diolein preparation previously. Herein, molecular modeling was used to elucidate the underlying mechanism of the solvent effect on the positional selectivity of the enzyme. The results showed that the binding energy of sn-1 hydroxyl of glycerol molecular with CALB became higher, and the binding energy of sn-2 hydroxyl of glycerol molecular with CALB became lower along with the increase of the solvent log P. It was demonstrated that, increasing log P of the solvent, the enzyme selectivity to sn-1 hydroxyl of glycerol molecular grew weaker, and the selectivity to sn-2 hydroxyl of glycerol molecular grew stronger.     

Discover potential inhibitors for PFKFB3 using 3D-QSAR,virtual screening,molecular docking and molecular dynamics simulation

Abstract

The 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase-3 (PFKFB3) is a master regulator of glycolysis in cancer cells by synthesizing fructose-2,6-bisphosphate (F-2,6-BP), a potent allosteric activator of phosphofructokinase-1 (PFK-1), which is a rate-limiting enzyme of glycolysis. PFKFB3 is an attractive target for cancer treatment. It is valuable to discover promising inhibitors by using 3D-QSAR pharmacophore modeling, virtual screening, molecular docking and molecular dynamics simulation. Twenty molecules with known activity were used to build 3D-QSAR pharmacophore models. The best pharmacophore model was ADHR called Hypo1, which had the highest correlation value of 0.98 and the lowest RMSD of 0.82. Then, the Hypo1 was validated by cost value method, test set method and decoy set validation method. Next, the Hypo1 combined with Lipinski's rule of five and ADMET properties were employed to screen databases including Asinex and Specs, total of 1,048,159 molecules. The hits retrieved from screening were docked into protein by different procedures including HTVS, SP and XP. Finally, nine molecules were picked out as potential PFKFB3 inhibitors. The stability of PFKFB3-lead complexes was verified by 40?ns molecular dynamics simulation. The binding free energy and the energy contribution of per residue to the binding energy were calculated by MM-PBSA based on molecular dynamics simulation.     

Computational investigations of gram-negative bacteria phosphopantetheine adenylyltransferase inhibitors using 3D-QSAR,molecular docking and molecular dynamic simulations

Abstract

Phosphopantetheine adenylyltransferase (PPAT) has been recognized as a promising target to develop novel antimicrobial agents, which is a hexameric enzyme that catalyzes the penultimate step in coenzyme A biosynthesis. In this work, molecular modeling study was performed with a series of PPAT inhibitors using molecular docking, three-dimensional qualitative structure-activity relationship (3D-QSAR) and molecular dynamic (MD) simulations to reveal the structural determinants for their bioactivities. Molecular docking study was applied to understand the binding mode of PPAT with its inhibitors. Subsequently, 3D-QSAR model was constructed to find the features required for different substituents on the scaffolds. For the best comparative molecular field analysis (CoMFA) model, the Q² and R² values of which were calculated as 0.702 and 0.989, while they were calculated as 0.767 and 0.983 for the best comparative molecular similarity index analysis model. The statistical data verified the significance and accuracy of our 3D-QSAR models. Furthermore, MD simulations were carried out to evaluate the stability of the receptor–ligand contacts in physiological conditions, and the results were consistent with molecular docking studies and 3D-QSAR contour map analysis. Binding free energy was calculated with molecular mechanics generalized born surface area approach, the result of which coincided well with bioactivities and demonstrated that van der Waals accounted for the largest portion. Overall, our study provided a valuable insight for further research work on the recognition of potent PPAT inhibitors.

Communicated by Ramaswamy H. Sarma     

Validation of potential inhibitors for SrtA against Bacillus anthracis by combined approach of ligand-based and molecular dynamics simulation

The development of SrtA inhibitors targeting the biothreat organism namely Bacillus anthracis was achieved by the combined approach of pharmacophore modeling, binding interactions, electron transferring capacity, ADME, and Molecular dynamics studies. In this study, experimentally reported Ba-SrtA inhibitors (pyridazinone and pyrazolethione derivatives) were considered for the development of enhanced pharmacophoric model. The obtained AAAHR hypothesis was a pure theoretical concept that accounts for common molecular interaction network present in experimentally active pyridazinone and pyrazolethione derivatives. Pharmacophore-based screening of AAAHR hypothesis provides several new compounds, and those compounds were treated with four phases of docking protocols with combined Glide-QPLD docking approach. In this approach, scoring and charge accuracy variations were seen to be dominated by QM/MM approach through the allocation of partial charges. Finally, we reported the best compounds from binding db, Chembridge db, and Toslab based on scoring values, energy parameters, electron transfer reaction, ADME, and cell adhesion inhibition activity. The dynamic state of interaction and binding energy assess that new compounds are more active inside the binding pocket and these compounds on experimental validations will survive as better inhibitors for targeting the cell adhesion mechanism of Ba-SrtA.     

Virtual screening of small molecules databases for discovery of novel PARP-1 inhibitors: combination of in silico and in vitro studies

Poly(ADP-ribose) polymerase-1 (PARP-1) enzyme has critical roles in DNA replication repair and recombination. Thus, PARP-1 inhibitors play an important role in the cancer therapy. In the current study, we have performed combination of in silico and in vitro studies in order to discover novel inhibitors against PARP-1 target. Structure-based virtual screening was carried out for an available small molecules database. A total of 257,951 ligands from Otava database were screened at the binding pocket of PARP-1 using high-throughput virtual screening techniques. Filtered structures based on predicted binding energy results were then used in more sophisticated molecular docking simulations (i.e. Glide/standard precision, Glide/XP, induced fit docking – IFD, and quantum mechanics polarized ligand docking – QPLD). Potential high binding affinity compounds that are predicted by molecular simulations were then tested by in vitro methods. Computationally proposed compounds as PARP-1 inhibitors (Otava Compound Codes: 7111620047 and 7119980926) were confirmed by in vitro studies. In vitro results showed that compounds 7111620047 and 7119980926 have IC₅₀ values of 0.56 and 63 μM against PARP-1 target, respectively. The molecular mechanism analysis, free energy perturbation calculations using long multiple molecular dynamics simulations for the discovered compounds which showed high binding affinity against PARP-1 enzyme, as well as structure-based pharmacophore development (E-pharmacophore) studies were also studied.     

Molecular docking,synthesis, and biological evaluation of naphthoquinone as potential novel scaffold for H5N1 neuraminidase inhibition

A series of dimeric naphthoquinones containing natural 2-hydroxy-1-4-naphthoquinone moiety was designed, synthesized, and evaluated against neuraminidase of H5N1 virus. p-hydroxy derivatives showed higher inhibition when compared to p-halogenated compounds. Molecular docking studies conducted with H5N1 neuraminidase clearly demonstrated different binding modes of the most active compound onto the open and closed conformations of loop 150. The results thus provide not only evidences of a novel scaffold evaluated as inhibitor, but also a rational explanation involving molecular modeling and the role of loop 150 in the binding.     

Insights into the structure and inhibition of Giardia intestinalis arginine deiminase: homology modeling,docking, and molecular dynamics studies

Giardia intestinalis arginine deiminase (GiADI) is an important metabolic enzyme involved in the energy production and defense of this protozoan parasite. The lack of this enzyme in the human host makes GiADI an attractive target for drug design against G. intestinalis. One approach in the design of inhibitors of GiADI could be computer-assisted studies of its crystal structure, such as docking; however, the required crystallographic structure of the enzyme still remains unresolved. Because of its relevance, in this work, we present a three-dimensional structure of GiADI obtained from its amino acid sequence using the homology modeling approximation. Furthermore, we present an approximation of the most stable dimeric structure of GiADI identified through molecular dynamics simulation studies. An in silico analysis of druggability using the structure of GiADI was carried out in order to know if it is a good target for design and optimization of selective inhibitors. Potential GiADI inhibitors were identified by docking of a set of 3196 commercial and 19 in-house benzimidazole derivatives, and molecular dynamics simulation studies were used to evaluate the stability of the ligand–enzyme complexes.     

Molecular modeling of human lanosterol 14α-demethylase complexes with substrates and their derivatives

Lanosterol 14α-demethylase (CYP51A1) is a key enzyme in sterol biosynthesis. In humans, this enzyme is involved in the cholesterol biosynthesis pathway. The majority of antifungal drugs are aimed at the inhibition of CYP51 in fungi. To elucidate the molecular mechanisms of highly specific protein-ligand recognition, we have developed a full-atomic model of human CYP51A1 and performed docking of natural substrates and their derivatives to the active site of the enzyme. The parameters of the binding enthalpy of substrates, intermediates, and final products of the reaction of 14α-demethylation were estimated using the MMPB(GB)SA algorithm. Dynamic properties and conformational changes of the protein globule upon binding of the ligand near the active site have been investigated by the molecular dynamics method. Our studies reveal that hydroxylated intermediate reaction products have a greater affinity than the initial substrates, which facilitates the multistage reaction without accumulation of intermediate products. The contribution to the free energy of steroid ligand binding of 30 amino acids forming the substrate-binding region of CYP51A1, as well as the influence of their substitutions to alanine on the stability of the protein molecule, has been clarified using alanine scanning modeling. We demonstrate that the most serious weakening of the binding is observed in the case of substitutions Y137A, F145A, V149A, I383A, and R388A. The results of molecular modeling are in agreement with the data obtained by analysis of primary sequences of representatives of the CYP51 family.     

Computational analysis of spiro-oxindole inhibitors of the MDM2-p53 interaction: insights and selection of novel inhibitors

Since MDM2 is an inhibitor of the p53 tumor suppressor, disrupting the MDM2-p53 interaction is a promising approach for cancer therapy. Here, we used molecular dynamics simulations followed by free energy decomposition analysis to study conformational changes in MDM2 induced by three known spiro-oxindole inhibitors. Analysis of individual energy terms suggests that van der Waals and electrostatic interactions explain much of the binding affinities of these inhibitors. Binding free energies calculated for the three inhibitors using the molecular mechanics-generalized Born surface area model were consistent with experimental data, suggesting the validity of this approach. Based on this structure-function analysis, several novel spiro-oxindole derivatives were selected and evaluated for their ability to block the MDM2-p53 interaction in vitro. These results suggest that combining in silico and experimental techniques can provide insights into the structure-function relationships of MDM2 inhibitors and guide the rational design of anticancer drugs targeting the MDM2-p53 interaction.     

An optimized MM/PBSA virtual screening approach applied to an HIV-1 gp41 fusion peptide inhibitor

VIRus Inhibitory Peptide (VIRIP), a 20 amino acid peptide, binds to the fusion peptide (FP) of human immunodeficiency virus type 1 (HIV-1) gp41 and blocks viral entry. VIRIP derivatives with improved antiviral activity have been developed, and one of those derivatives has recently proven effective and safe in a phase 1/2 clinical trial. Here, molecular dynamics were executed in combination with molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) free energy calculations to explore the binding interaction between VIRIP derivatives and gp41 FP. A promising correlation between antiviral activity and simulated binding free energy was established thanks to restriction of the flexibility of the peptides, inclusion of configurational entropy calculations, and the use of multiple internal dielectric constants for the MM/PBSA calculations depending on the amino acid sequence. Based on these results, a virtual screening experiment was carried out to design VIRIP analogs with further improved antiretroviral activity. A selection of peptides was tested for inhibitory activity and several VIRIP derivatives were identified with significantly enhanced activity compared to the reference peptides. The results demonstrate that computational modeling strategies using an adapted MM/PBSA methodology improve the accuracy of binding free energy calculations of peptide complexes compared to the classic MM/PBSA protocol. As such, this virtual screening approach generated HIV-1 gp41 FP inhibitors with improved antiviral activity that could be useful for future clinical applications.     

Probing the binding mechanism of mercaptoguanine derivatives as inhibitors of HPPK by docking and molecular dynamics simulations

6-Hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) is a promising antimicrobial target involved in the folate biosynthesis pathway. Although, the results from crystallographic studies of HPPK have attracted a great interest in the design of novel HPPK inhibitors, the mechanism of action of HPPK due to inhibitor binding remains questionable. Recently, mercaptoguanine derivatives were reported to inhibit the pyrophosphoryl transfer mechanism of Staphylococcus aureus HPPK (SaHPPK). The present study is an attempt to understand the SaHPPK-inhibitors binding mechanism and to highlight the key residues that possibly involve in the complex formation. To decipher these questions, we used the state-of-the-art advanced insilico approach such as molecular docking, molecular dynamics (MD), molecular mechanics-generalized Born surface area approach. Domain cross correlation and principle component analysis were applied to the snapshots obtained from MD revealed that the compounds with high binding affinity stabilize the conformational dynamics of SaHPPK. The binding free energy estimation showed that the van der Waals and electrostatic interactions played a vital role for the binding mechanism. Additionally, the predicted binding free energy was in good agreement with the experimental values (R² = .78). Moreover, the free energy decomposition on per-residue confirms the key residues that significantly contribute to the complex formation. These results are expected to be useful for rational design of novel SaHPPK inhibitors.     

Molecular interaction of ctDNA and HSA with sulfadiazine sodium by multispectroscopic methods and molecular modeling

Interactions of sulfadiazine sodium (SD‐Na) with calf thymus DNA (ctDNA) and human serum albumin (HSA) were studied using fluorescence spectroscopy, UV absorption spectroscopy and molecular modeling. The fluorescence experiments showed that the processes were static quenching. The results of UV spectra and molecular modeling of the interaction between SD‐Na and ctDNA indicated that the binding mode might be groove binding. In addition, the interaction of SD‐Na with HSA under simulative physiological conditions was also investigated. The binding constants (K) and the number of binding sites (n) at different temperatures (292, 302, 312 K) were 5.23 × 10³ L/mol, 2.18; 4.50 × 10³ L/mol, 2.35; and 4.08 × 10³ L/mol, 2.47, respectively. Thermodynamic parameters including enthalpy change (ΔH) and entropy change (ΔS) were calculated, the results suggesting that hydrophobic force played a very important role in SD‐Na binding to HSA, which was in good agreement with the molecular modeling study. Moreover, the effect of SD‐Na on the conformation of HSA was analyzed using three‐dimensional fluorescence spectra. Copyright © 2013 John Wiley & Sons, Ltd.     

Ugi Bis-Amide Derivatives as Tyrosinase Inhibitor; Synthesis,Biology Assessment,and in Silico Analysis

Herein, a straightforward synthetic strategy mediated by Ugi reaction was developed to synthesize novel series of compounds as tyrosinase inhibitors. The structures of all compounds were confirmed by FT-IR, ¹H-NMR, ¹³C-NMR, and CHNOS techniques. The tyrosinase inhibitory activities of all synthesized derivatives 5a – m were determined against mushroom tyrosinase and it was found that derivative 5c possesses the best inhibition with an IC₅₀ value of 69.53±0.042 μM compared to the rest of the synthesized derivatives. Structure–activity relationships (SARs) showed that the presence of 4-MeO or 4-NO₂ at the R² position plays a key role in tyrosinase inhibitory activities. The enzyme kinetics studies showed that compound 5c is an noncompetitive inhibitor. For in silico study, the allosteric site detection was first applied to find the appropriate binding site and then molecular docking and molecular dynamic studies were performed to reveal the position and interactions of 5c as the most potent inhibitor within the tyrosinase active site. The results showed that 5c bind well with the proposed binding site and formed a stable complex with the target protein.     

Insight into the structural requirements of thiophene-3-carbonitriles-based MurF inhibitors by 3D-QSAR,molecular docking and molecular dynamics study

The discovery of clinically relevant inhibitors against MurF enzyme has proven to be a challenging task. In order to get further insight into the structural features required for the MurF inhibitory activity, we performed pharmacophore and atom-based three-dimensional quantitative structure–activity relationship studies for novel thiophene-3-carbonitriles based MurF inhibitors. The five-feature pharmacophore model was generated using 48 inhibitors having IC₅₀ values ranging from 0.18 to 663?μm. The best-fitted model showed a higher coefficient of determination (R²?=?0.978), cross-validation coefficient (Q²?=?0.8835) and Pearson coefficient (0.9406) at four component partial least-squares factor. The model was validated with external data set and enrichment study. The effectiveness of the docking protocol was validated by docking the co-crystallized ligand into the catalytic pocket of MurF enzyme. Further, binding free energy calculated by the molecular mechanics generalized Born surface area approach showed that van der Waals and non-polar solvation energy terms are the main contributors to ligand binding in the active site of MurF enzyme. A 10-ns molecular dynamic simulation was performed to confirm the stability of the 3ZM6-ligand complex. Four new molecules are also designed as potent MurF inhibitors. These results provide insights regarding the development of novel MurF inhibitors with better binding affinity.     

A combination of pharmacophore modeling,atom-based 3D-QSAR,molecular docking and molecular dynamics simulation studies on PDE4 enzyme inhibitors

Phosphodiesterases 4 enzyme is an attractive target for the design of anti-inflammatory and bronchodilator agents. In the present study, pharmacophore and atom-based 3D-QSAR studies were carried out for pyrazolopyridine and quinoline derivatives using Schrödinger suite 2014-3. A four-point pharmacophore model was developed using 74 molecules having pIC₅₀ ranging from 10.1 to 4.5. The best four feature model consists of one hydrogen bond acceptor, two aromatic rings, and one hydrophobic group. The pharmacophore hypothesis yielded a statistically significant 3D-QSAR model, with a high correlation coefficient (R^2?=?.9949), cross validation coefficient (Q^2?=?.7291), and Pearson-r (.9107) at six component partial least square factor. The external validation indicated that our QSAR model possessed high predictive power with R² value of .88. The generated model was further validated by enrichment studies using the decoy test. Molecular docking, free energy calculation, and molecular dynamics (MD) simulation studies have been performed to explore the putative binding modes of these ligands. A 10-ns MD simulation confirmed the docking results of both stability of the 1XMU–ligand complex and the presumed active conformation. Outcomes of the present study provide insight in designing novel molecules with better PDE4 inhibitory activity.     

Insight into the key active sites of afChiA1 based on molecular dynamics simulations and free energy calculations

In this study, the binding of the enzyme chitinase A1 (afChiA1) from the plant-type Aspergillus fumigatus with four potent inhibitors, allosamidin (ASM), acetazolamide (AZM), 8-chloro-theophylline (CTP) and kinetin (KIT) is investigated by molecular docking, molecular dynamics simulation and binding free energy calculation. The results reveal that the electrostatic interactions play an important role in the stabilisation of the binding of afChiA1 with inhibitors. Based on the binding energy of afChiA1-ligands, the key residues (Gln37 and Trp312) in the active binding pocket of the complex systems are confirmed by molecular mechanics/Poisson–Boltzmann surface area method, and the active inhibitors, ASM and AZM, both could form strong interaction with Gln37 and Trp312, and the non-active ligands, CTP and KIT, could not interact with these two residues, which is consistent with the result of experimental report. Then, it is identified that Gln37 and Trp312 should be one of the important active site residues of afChiA1.