Saad Raza

Binding mode analysis, dynamic simulation and binding free energy calculations of the MurF ligase from Acinetobacter baumannii

MurF ligase catalyzes the final cytoplasmic step of bacterial peptidoglycan biosynthesis and, as such, is a validated target for therapeutic intervention. Herein, we performed molecular docking to identify putative inhibitors of Acinetobacter baumannii MurF (AbMurF). Based on comparative docking analysis, compound 114 (ethyl pyridine substituted 3-cyanothiophene) was predicted to potentially be the most active ligand. Computational pharmacokinetic characterization of drug-likeness of the compound showed it to fulfil all the parameters of Muegge and the MDDR rule. A molecular dynamic simulation of 114 indicated the complex to be stable on the basis of an average root mean square deviation (RMSD) value of 2.09Å for the ligand. The stability of the complex was further supported by root mean square fluctuation (RMSF), beta factor and radius of gyration values. Analyzing the complex using radial distribution function (RDF) and a novel analytical tool termed the axial frequency distribution (AFD) illustrated that after simulation the ligand is positioned in close vicinity of the protein active site where Thr42 and Asp43 participate in hydrogen bonding and stabilization of the complex. Binding free energy calculations based on the Poisson-Boltzmann or Generalized-Born Surface Area Continuum Solvation (MM(PB/GB)SA) method indicated the van der Waals contribution to the overall binding energy of the complex to be dominant along with electrostatic contributions involving the hot spot amino acids from the protein active site. The present results indicate that the screened compound 114 may act as a parent structure for designing potent derivatives against AbMurF in specific and MurF of other bacterial pathogens in general.

A one-pot multicomponent facile synthesis of dihydropyrimidin-2(1H)-thione derivatives using triphenylgermane as a catalyst and its binding pattern validation

A series of substituted dihydropyrimidin-2(1H)-thione derivatives (1–8) have been synthesized using a facile and modified procedure with triphenylgermyl propionate as a catalyst. In comparison with the classical Biginelli reaction, this new protocol has the advantages of excellent yield and shorter reaction times. The synthesized compounds have been characterized by various spectroscopic techniques such as FT-IR, multinuclear (1H/13C) NMR spectroscopy and single crystal XRD analysis. Molecular docking studies were performed to identify the probable binding modes of potent inhibitors in the active site of the enzymes human topoisomerase II alpha (4FM9) and Helicobacter pylori urease (1E9Y). Compound 3 was found to be the most potent inhibitor according to the molecular docking scores and molecular dynamic simulations which suggests it can be further processed as a lead molecule to interpret the pharmacological properties of these compounds.

Molecular dynamics simulation studies of novel β-lactamase inhibitor

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

Binding pattern analysis and structural insight into the inhibition mechanism of Sterol 24-C methyltransferase by docking and molecular dynamics approach.

Sterol 24-C methyltransferase (SMT) plays a major role during the production of steroids, especially in the biosynthesis of ergosterol, which is the major membrane sterol in leishmania parasite, and the etiological basis of leishmaniasis. Mechanism-based inactivators bind irreversibly to SMT and interfere with its activity to provide leads for the design of antileishmanial inhibitors. In this study, computational methods are used for studying enzyme–inhibitor interactions. fifty-seven mechanism-based inactivators are docked using 3 docking/scoring approaches (FRED, GoldScore, and ChemScore). A consensus is generated from the results of different scoring functions which are also validated with already reported experimental values. The most active compound thus obtained is subjected to molecular dynamics simulation of length 20 ns. Stability of simulation is analyzed through root-mean-square deviation, beta factor (B-factor), and radius of gyration (Rg). Hydrogen bonds and their involvement in the structural stability of the enzyme are evaluated through radial distribution function. Newly developed application of axial frequency distribution that determines three-particle correlation on frequency distributions before and after simulation has provided a clear evidence for the movement of the inhibitor into active pocket of the enzyme. Results yielded strong interaction between enzyme and the inhibitor throughout the simulation. Binding of the inhibitor with enzyme has stabilized the enzyme structure; thus, the inhibitor has the potential to become a lead compound.

Structure and dynamics studies of sterol 24-C-methyltransferase with mechanism based inactivators for the disruption of ergosterol biosynthesis

The enzyme sterol 24-C-methyltransferase (SMT) belongs to the family of transferases, specifically to the one-carbon transferring methyltransferases. SMT has been found playing a major role during the production of steroids, especially for the biosynthesis of ergosterol, which is the major membrane sterol in leishmania parasites, causing leishmaniasis. However, SMT and ergosterol are not found in mammals, so, an extensive study has been carried out over the susceptible SMT protein, which is found to be highly conserved among all the Leishmania species and holds a significant anti-leishmanial drug target. To date, there is no computational data available for SMT, due to its highly unexplored profile. In this work, a complete set of structural attributes have been examined through the available computational procedures, along with an attempt to characterize the most capable modeling server available. The exploration ranges from physicochemical characterization, pairwise alignment, secondary structure prediction, to active site detection. With this information, a docking study was carried out to find the compound that best binds into the active site. Moreover, molecular dynamics simulation was conducted to examine the stability of the homology modeled protein and the ligand–enzyme complex. The results indicate that the ligand–enzyme complex is more stable.

Towards novel inhibitors against KdsB: A highly specific and selective broad-spectrum bacterial enzyme

KdsB (3-deoxy-manno-octulosonate cytidylyltransferase) is a highly specific and selective bacterial enzyme that catalyzes KDO (3-Deoxy-D-mano-oct-2-ulosonic acid) activation in KDO biosynthesis pathway. Failure in KDO biosynthesis causes accumulation of lipid A in the bacterial outer membrane that leads to cell growth arrest. This study reports a combinatorial approach comprising virtual screening of natural drugs library, molecular docking, computational pharmacokinetics, molecular dynamics simulation, and binding free energy calculations for the identification of potent lead compounds against the said enzyme. Virtual screening demonstrated 1460 druglike compounds in a total of 4800, while molecular docking illustrated Ser13, Arg14, and Asp236 as the anchor amino acids for recognizing and binding the inhibitors. Functional details of the enzyme in complex with the best characterized compound-226 were explored through two hundred nanoseconds of MD simulation. The ligand after initial adjustments jumps into the active cavity, followed by the deep cavity, and ultimately backward rotating movement toward the initial docked site of the pocket. During the entire simulation period, Asp236 remained in contact with the ligand and can be considered as a major catalytic residue of the enzyme. Radial distribution function confirmed that toward the end of the simulation, strengthening of ligand-receptor occurred with ligand and enzyme active residues in close proximity. Binding free energy calculations via MM(PB/GB)SA and Waterswap reaction coordinates, demonstrated the high affinity of the compound for enzyme active site residues. These findings can provide new avenues for designing potent compounds against notorious bacterial pathogens.

Improvement in thermostability of xylanase from Geobacillus thermodenitrificans C5 by site directed mutagenesis

Enzymes activity and stability at extreme temperature can be intensified by regularly applying protein engineering. In the present study, two amino acids were perceived to mark the temperature dependability of xylanase from Geobacillus thermodenitrificans C5. Six mutants of G. thermodenitrificans C5 were built through site-directed mutagenesis by interchanging the residue with proline and glutamic acid (R81P, H82E, W185P, D186E, double mutant W185P/D186E and triple mutant H82E/W185P/D186E). Both mutant and wild type enzymes were quantified in host E. coli BL21. In comparison to wild type, the temperature was enhanced by 4 °C, 5 °C and 11 °C in H82E, W185P/D186E and H82E/W185P/D186E mutant models, respectively. The mutant H82E and the combined substitutions (H82E/W185P/D186E) showed the most pronounced shifts in their half-lives for thermal inactivation. Half-life was increased 13 times at 60 °C, 15 times at 65 °C, 9 times at 70 °C and 5 times at 75 °C by H82E/W185P/D186E mutant. Mutations in xylanase enzyme causes rigidification of essential chain and filling of groove that leads to stabilization of mutants and finally resulted into enhancement in their thermostability.

Structure modeling and hybrid virtual screening study of Alzheimer’s associated protease kallikrein 8 for the identification of novel inhibitors

Kallikrein 8 (KLK8) is an extra cellular serine protease which is responsible for nerve growth, and degeneration and nervous plasticity have been associated with Alzheimer’s disease. In silico analysis of KLK8 has not been performed until now. This study is aimed at molecular modeling and lead identification for a potent inhibitor. The hybrid of ligand- and structure-based virtual screening was applied on commercially available compounds. Ligand similarity search was employed, followed by ligand docking protocol. Compound’s potential for their activity was deduced from docking scores and their interactions with the active site. For active permeation of compounds through blood brain barrier, their molecular properties were checked with previously reported compounds. The compound that showed the most potential according to criteria was 1-(3,5-difluorophenyl)-5-hydroxy-7-(4-hydroxy-3,5-dimethoxyphenyl)-6,7-dihydro-1H-pyrrolo[3,2-b] pyridine-3-carboxylic acid (ZINC 61720639).

AFD: an application for bi-molecular interaction using axial frequency distribution

Conformational flexibility and generalized structural features are responsible for specific phenomena existing in biological pathways. With advancements in computational chemistry, novel approaches and new methods are required to compare the dynamic nature of biomolecules, which are crucial not only to address dynamic functional relationships but also to gain detailed insights into the disturbance and positional fluctuation responsible for functional shifts. Keeping this in mind, axial frequency distribution (AFD) has been developed, designed, and implemented. AFD can profoundly represent distribution and density of ligand atom around a particular atom or set of atoms. It enabled us to obtain an explanation of local movements and rotations, which are not significantly highlighted by any other structural and dynamical parameters. AFD can be implemented on biological models representing ligand and protein interactions. It shows a comprehensive view of the binding pattern of ligand by exploring the distribution of atoms relative to the x-y plane of the system. By taking a relative centroid on protein or ligand, molecular interactions like hydrogen bonds, van der Waals, polar or ionic interaction can be analyzed to cater the ligand movement, stabilization or flexibility with respect to the protein. The AFD graph resulted in the residual depiction of bi-molecular interaction in gradient form which can yield specific information depending upon the system of interest.

Moleculer dynamics simulaiton revealed reciever domain of Acinetobacter baumannii BfmR enzyme as the hot spot for future antibiotics designing

Abstract Acinetobacter baumannii is an alarming nosocomial pathogen that is resistant to multiple drugs. The pathogen is forefront of scientific attention because of high mortality and morbidity found for its complications in the past decade. As a consequence, identification of novel drug candidates and subsequent designing of novel chemical scaffolds is an imperative need of time. In the present study, we used a recently reported structure of BfmR enzyme and performed structure based virtual screening, MD simulation and binding free energies calculations. MD simulation revealed a profound movement of the best-characterized inhibitor towards the α4-β5-α5 face of the enzyme receiver domain, thus indicating its high affinity for this site compared to phosphorylation. Furthermore, it was observed that the enzyme and enzyme-inhibitor complex have high structure stability with mean RMSD of 1.2 and 1.1 Å, respectively. Binding free energy calculations for the complex unraveled high stability with MMGBSA score of −26.21 kcal/mol and MMPBSA score of −1.47 kcal/mol. Van der Waal energy was found highly favorable with value of −30.25 kcal/mol and dominated significantly the overall binding energy. Furthermore, a novel WaterSwap assay was used to circumvent the limitations of MMGB/PBSA that complements the inhibitor affinity for enzyme active pocket as depicted by the low convergence of Bennett, TI and FEP algorithms. Results yielded from this study will not only give insight into the phenomena of inhibitor movement towards the enzyme receiver domain, but will also provide a useful baseline for designing derivatives with improved biological and pharmacokinetics profiles. Communicated by Ramaswamy H. Sarma