Discovery of ZAP70 inhibitors by high-throughput docking into a conformation of its kinase domain generated by molecular dynamics

Very few selective inhibitors of the zeta-chain associated protein kinase 70 kDa (ZAP70) have been reported despite its importance in autoimmune diseases. Here, to induce a fit of the so-called gatekeeper residue (Met414) and hydrophobic pocket next to it, a potent Janus kinase 2 (JAK2) inhibitor was first docked into the ATP binding site of ZAP70 by structural alignment of the kinase domains. The resulting model of the complex between ZAP70 and the JAK2 inhibitor was then relaxed by an explicit solvent molecular dynamics simulation with restraints on the backbone atoms. High-throughput docking into the induced-fit conformation of ZAP70 generated by molecular dynamics has revealed 10 low-micromolar inhibitors which correspond to six distinct chemotypes. One of these ZAP70 inhibitors has an IC₅₀ of 110 nM for JAK2.     

Analysis of cyclosporin A and a set of analogs as inhibitors of a T. cruzi cyclophilin by docking and molecular dynamics

Cyclophilins (CyPs) are enzymes involved in protein folding. In Trypanosoma cruzi (T. cruzi), the most abundantly expressed CyP is the isoform TcCyP19. It has been shown that TcCyP19 is inhibited by the immunosuppressive drug cyclosporin A (CsA) and analogs, which also proved to have potent trypanosomicidal activity in vitro. In this work, we continue and expand a previous study on the molecular interactions of CsA, and a set of analogs modeled in complexes with TcCyP19. The modeled complexes were used to evaluate binding free energies by molecular dynamics (MD), applying the Linear Interaction Energy (LIE) method. In addition, putative binding sites were identified by molecular docking. In our analysis, the binding free energy calculations did not correlate with experimental data. The heterogeneity of the non-bonded energies and the variation in the pattern of hydrogen bonds suggest that the systems may not be suitable for the application of the LIE method. Further, the docking calculations identified two other putative binding sites with comparable scoring energies to the active site, a fact that may also explain the lack of correlation found. Kinetic experiments are needed to confirm or reject the multiple binding sites hypothesis. In the meantime, MD simulations at the alternative sites, employing other methods to compute binding free energies, might be successful at finding good correlations with the experimental data.     

New insights into the effect of mutations on affibody-Fc interaction,a molecular dynamics simulation approach

Staphylococcal protein A (SpA) domain B (the basis of affibody) has been widely used in affinity chromatography and found therapeutic applications against inflammatory diseases through targeting the Fc part of immunoglobulin G (IgG). We have performed extensive molecular dynamics simulation of 41 SpA mutants and compared their dynamics and conformations to wild type. The simulations revealed the molecular details of structural and dynamics changes that occurred due to introducing point mutations and helped to explain the SPR results. It was observed in some variants a point mutation caused extensive structural changes far from the mutation site, while an effect of some other mutations was limited to the site of the mutated residue. Also, the pattern of hydrogen bond networks and hydrophobic core arrangements were investigated. We figured out mutations that occurred at positions 128, 136, 150 and 153, affected two hydrophobic cores at the interface as well as mutations introduced at positions 129 and 154 interrupted two hydrogen bond networks of the interface, SPR data showed all of these mutations reduced binding affinity significantly. Overall, by scanning the SpA-Fc interface through the large numbers of introduced mutations, the new insights have been gained which would help to design high- affinity ligands of IgG.     

2D-QSAR study,molecular docking,and molecular dynamics simulation studies of interaction mechanism between inhibitors and transforming growth factor-beta receptor I (ALK5)

Transforming growth factor type 1 receptor (ALK5) is kinase associated with a wide variety of pathological processes, and inhibition of ALK5 is a good strategy to treat many kinds of cancer and fibrotic diseases. Recently, a series of compounds have been synthesized as ALK5 inhibitors. However, the study of their selectivity against other potential targets remains elusive. In this research, a data-set of ALK5 inhibitors were collected and studied based on the combination of 2D-QSAR, molecular docking and molecular dynamics simulation. The quality of QSAR models were assessed statistically by F, R², and R²_ADJ, proved to be credible. The cross-validations for the models (q²_LOO = 0.571 and 0.629, respectively) showed their robustness, while the external validations (r²_test = 0.703 and 0.764, respectively) showed their predictive power. Besides, the predicted binding free energy results calculated by MM/GBSA method were in accordance with the experimental data, and the van der Waals energy term was the factor that had the most significant impact on ligand binding. What is more, several important residues were found to significantly affect the binding affinity. Finally, based on our analyses above, a proposed series of molecules were designed.     

High throughput screening,docking, and molecular dynamics studies to identify potential inhibitors of human calcium/calmodulin-dependent protein kinase IV

Calcium/calmodulin-dependent protein kinase IV (CAMKIV) is associated with many diseases including cancer and neurodegenerative disorders and thus being considered as a potential drug target. Here, we have employed the knowledge of three-dimensional structure of CAMKIV to identify new inhibitors for possible therapeutic intervention. We have employed virtual high throughput screening of 12,500 natural compounds of Zinc database to screen the best possible inhibitors of CAMKIV. Subsequently, 40 compounds which showed significant docking scores (?11.6 to ?10.0?kcal/mol) were selected and further filtered through Lipinski rule and drug likeness parameter to get best inhibitors of CAMKIV. Docking results are indicating that ligands are binding to the hydrophobic cavity of the kinase domain of CAMKIV and forming a significant number of non-covalent interactions. Four compounds, ZINC02098378, ZINC12866674, ZINC04293413, and ZINC13403020, showing excellent binding affinity and drug likeness were subjected to molecular dynamics simulation to evaluate their mechanism of interaction and stability of protein-ligand complex. Our observations clearly suggesting that these selected ligands may be further employed for therapeutic intervention to address CAMKIV associated diseases.

Communicated by Ramaswamy H. Sarma     

Exploring the molecular basis for selective binding of Mycobacterium tuberculosis Asp kinase toward its natural substrates and feedback inhibitors: A docking and molecular dynamics study

Tuberculosis (TB) is still a major public health problem, compounded by the human immunodeficiency virus (HIV)-TB co-infection and recent emergence of multidrug-resistant (MDR) and extensively drug resistant (XDR)-TB. In this context, aspartokinase of mycobacterium tuberculosis has drawn attention for designing novel anti-TB drugs. Asp kinase is an enzyme responsible for the synthesis of 4-phospho-L-aspartate from L-aspartate and involved in the branched biosynthetic pathway leading to the synthesis of amino acids lysine, threonine, methionine and isoleucine. An intermediate of lysine biosynthetic branch, mesodiaminopimelate is also a component of the peptidoglycan which is a component of bacterial cell wall. To interfere with the production of all these amino acids and cell wall, it is possible to inhibit Asp kinase activity. This can be achieved using Asp kinase inhibitors. In order to design novel Asp kinase inhibitors as effective anti-TB drugs, it is necessary to have an understanding of the binding sites of Asp kinase. As no crystal structure of the enzyme has yet been published, we built a homology model of Asp kinase using the crystallized Asp kinase from M. Jannaschii, as template structures (2HMF and 3C1M). After the molecular dynamics refinement, the optimized homology model was assessed as a reliable structure by PROCHECK, ERRAT, WHAT-IF, PROSA2003 and VERIFY-3D. The results of molecular docking studies with natural substrates, products and feedback inhibitors are in agreement with the published data and showed that ACT domain plays an important role in binding to ligands. Based on the docking conformations, pharmacophore model can be developed by probing the common features of ligands. By analyzing the results, ACT domain architecture, certain key residues that are responsible for binding to feedback inhibitors and natural substrates were identified. This would be very helpful in understanding the blockade mechanism of Asp kinase and providing insights into rational design of novel Asp kinase inhibitors for M.tuberculosis.     

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.     

A multi-spectroscopic and molecular docking approach to investigate the interaction of antiviral drug oseltamivir with ct-DNA

The possible interaction between the antiviral drug oseltamivir and calf thymus DNA at physiological pH was studied by spectrophotometry, competitive spectrofluorimetry, differential pulse voltammogram (DPV), circular dichroism spectroscopy (CD), viscosity measurements, salt effect, and computational studies. Intercalation of oseltamivir between the base pairs of DNA was shown by a sharp increase in specific viscosity of DNA and a decrease of the peak current and a positive shift in differential pulse voltammogram. Competitive fluorescence experiments were performed using neutral red (NR) as a probe for the intercalation binding mode. The studies showed that oseltamivir is able to release the NR.     

Investigating the structural impacts of I64T and P311S mutations in APE1-DNA complex: a molecular dynamics approach

Background

Elucidating the molecular dynamic behavior of Protein-DNA complex upon mutation is crucial in current genomics. Molecular dynamics approach reveals the changes on incorporation of variants that dictate the structure and function of Protein-DNA complexes. Deleterious mutations in APE1 protein modify the physicochemical property of amino acids that affect the protein stability and dynamic behavior. Further, these mutations disrupt the binding sites and prohibit the protein to form complexes with its interacting DNA.

Principal Findings

In this study, we developed a rapid and cost-effective method to analyze variants in APE1 gene that are associated with disease susceptibility and evaluated their impacts on APE1-DNA complex dynamic behavior. Initially, two different in silico approaches were used to identify deleterious variants in APE1 gene. Deleterious scores that overlap in these approaches were taken in concern and based on it, two nsSNPs with IDs rs61730854 (I64T) and rs1803120 (P311S) were taken further for structural analysis.

Significance

Different parameters such as RMSD, RMSF, salt bridge, H-bonds and SASA applied in Molecular dynamic study reveals that predicted deleterious variants I64T and P311S alters the structure as well as affect the stability of APE1-DNA interacting functions. This study addresses such new methods for validating functional polymorphisms of human APE1 which is critically involved in causing deficit in repair capacity, which in turn leads to genetic instability and carcinogenesis.     

Predicting the impact of mutations on the specific activity of Bacillus thermocatenulatus lipase using a combined approach of docking and molecular dynamics

Lipases are important biocatalysts owing to their ability to catalyze diverse reactions with exceptional substrate specificities. A combined docking and molecular dynamics (MD) approach was applied to study the chain‐length selectivity of Bacillus thermocatenulatus lipase (BTL2) towards its natural substrates (triacylglycerols). A scoring function including electrostatic, van der Waals (vdW) and desolvation energies along with conformational entropy was developed to predict the impact of mutation. The native BTL2 and its 6 mutants (F17A, V175A, V175F, D176F, T178V and I320F) were experimentally analyzed to determine their specific activities towards tributyrin (C4) or tricaprylin (C8), which were used to test our approach. Our scoring methodology predicted the chain‐length selectivity of BTL2 with 85.7% (6/7) accuracy with a positive correlation between the calculated scores and the experimental activity values (r = 0.82, p = 0.0004). Additionally, the impact of mutation on activity was predicted with 75% (9/12) accuracy. The described study represents a fast and reliable approach to accurately predict the effect of mutations on the activity and selectivity of lipases and also of other enzymes. Copyright © 2016 John Wiley & Sons, Ltd.     

An integrated approach of network-based systems biology,molecular docking,and molecular dynamics approach to unravel the role of existing antiviral molecules against AIDS-associated cancer

A serious challenge in cancer treatment is to reposition the activity of various already known drug candidates against cancer. There is a need to rewrite and systematically analyze the detailed mechanistic aspect of cellular networks to gain insight into the novel role played by various molecules. Most Human Immunodeficiency Virus infection-associated cancers are caused by oncogenic viruses like Human Papilloma Viruses and Epstein–Bar Virus. As the onset of AIDS-associated cancers marks the severity of AIDS, there might be possible interconnections between the targets and mechanism of both the diseases. We have explored the possibility of certain antiviral compounds to act against major AIDS-associated cancers: Kaposi’s Sarcoma, Non-Hodgkin Lymphoma, and Cervical Cancer with the help of systems pharmacology approach that includes screening for targets and molecules through the construction of a series of drug–target and drug–target–diseases network. Two molecules (Calanolide A and Chaetochromin B) and the target “HRAS” were finally screened with the help of molecular docking and molecular dynamics simulation. The results provide novel antiviral molecules against HRAS target to treat AIDS defining cancers and an insight for understanding the pharmacological, therapeutic aspects of similar unexplored molecules against various cancers.     

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.     

Identification of curcumin derivatives as human glyoxalase I inhibitors: A combination of biological evaluation, molecular docking, 3D-QSAR and molecular dynamics simulation studies

Several recent developments suggest that the human glyoxalase I (GLO I) is a potential target for anti-tumor drug development. In present study, a series of curcumin derivatives with high inhibitory activity against human GLO I were discovered. Inhibition constant (K(i)) values of compounds 8, 9, 10, 11 and 13 to GLO I are 4.600μM, 2.600μM, 3.200μM, 3.600μM and 3.600μM, respectively. To elucidate the structural features of potent inhibitors, docking-based three-dimensional structure-activity relationship (3D-QSAR) analyses were performed. Satisfactory agreement between experiment and theory suggests that comparative molecular similarity index analysis (CoMSIA) modeling exhibit much better correlation and predictive power. The cross-validated q(2) value is 0.638 while no-validation r(2) value is 0.930. Integrated with docking-based 3D-QSAR CoMSIA modeling, molecular surface property (electrostatic and steric) mapping and molecular dynamics simulation, a set of receptor-ligand binding models and bio-affinity predictive models for rational design of more potent inhibitors of GLO I are established.     

Understanding the molecular interaction of human argonaute-2 and miR-20a complex: A molecular dynamics approach

Argonaute-2 (AGO2), a member of the Argonaute family, is the only member possessing catalytic and RNA silencing activity. In here, a molecular dynamics (MDs) simulation was performed using the crystal structure of human AGO2 protein complex with miR-20a. miR-20a is involved with various kind of biological process like heart and lung development, oncogenic process, etc. In precise, MD simulation was carried out with AGO2 protein complex with wild type, two mutant sites and four mutant sites in guided microRNA (miRNA). It has been noted that root-mean-square deviation (RMSD) of atomic positions of nucleic acid for wild type and two mutant sites guided miRNA has the same pattern of fluctuations, which stabilizes around 0.27 nm after 2 ns. Cα atom of AGO2 protein in the complex shows that this complex with wild type and two mutant site mutation duplex has a stable RMSD value after 20 ns, ranging between 0.14 and 0.21 nm. From the root-mean-square fluctuation (RMSF), we observed an increased pattern of fluctuations for the atoms of four mutant complex of AGO2-miR-20a complex. This increased RMSF of non-mutated nucleic acids is contributed by U-A bond breaking at the site of the nucleotide of U2 of guided miRNA, as observed from the duplex structure taken at different time steps of the simulation. Superimposed structure of the miRNA-mRNA duplex for the three complexes depicts that the three miRNA-mRNA duplexes are stable during the simulation. Current work demonstrates the possible correlations between the conformational changes of this AGO2-miR-20a duplex structure and the interactions of different atoms.     

Study on the interaction of fisetholz with BSA/HSA by multi-spectroscopic,cyclic voltammetric,and molecular docking technique

The interaction of fisetholz with bovine serum albumin (BSA) and human serum albumin (HSA) was investigated by multi-spectroscopic, cyclic voltammetric, and molecular docking technique. The results revealed that there was a static quenching of BSA/HSA induced by fisetholz. The binding constants (K_a) and binding sites (n) were calculated at different temperatures (293, 303, and 311?K). The enthalpy change (ΔH) were calculated to be –17.20?kJ mol^?1 (BSA) and –18.28?kJ mol^?1 (HSA) and the entropy change (ΔS) were calculated to be 35.41?J mol^?1 (BSA) and 24.02?J mol^?1 (HSA), respectively, which indicated that the interaction between fisetholz and BSA/HSA was mainly by electrostatic attraction. Based on displacement experiments using site probes, indomethacin and ibuprofen, the binding site of fisetholz to BSA/HSA was identified as sub-domain IIIA, which was further confirmed by molecular docking method. There was little effect of K⁺, Ca²⁺, Cu²⁺, Zn²⁺, and Fe³⁺ on fisetholz-BSA or fisetholz-HSA complex. The spectra of synchronous fluorescence, circular dichroism (CD) and Fourier transform infrared (FT-IR) all showed that fisetholz binding to BSA/HSA leads to secondary structures change of the two serum albumins. According to the Förster non-radiation energy transfer theory, the binding distance between fisetholz and BSA/HSA was 2.94/4.68?nm. The cyclic voltammetry as a supporting tool also indicated that fisetholz interacted with protein.

Communicated by Ramaswamy H. Sarma     

Structural requirements of pyrimidine, thienopyridine and ureido thiophene carboxamide-based inhibitors of the checkpoint kinase 1: QSAR, docking, molecular dynamics analysis

Our focus of current research is directed toward clarification of novel inhibitors (pyrazolo[1,5-a] pyrimidine (PP), thienopyridines (TP) and 2-ureido thiophene carboxamide (UTC) derivatives) targeting Checkpoint kinase 1 (CHK(1)), which is an oncology target of significant current interest. Our computational approaches include: (i) QSAR analysis was carried out on the computed steric/electrostatic/hydrophobic/hydrogen bond donor/hydrogen bond acceptor interactions with the pseudoreceptor surface, which yielded predictive models capable of explaining much of the variance of inhibitors. The resultant optimum QSAR/CoMFA models exhibited (N(training) = 51, N(test) = 16, R(cv) (2) = 0.47, R(pred) (2) = 0.7) for PP, (N(training) = 45, N(test) = 9, R(cv) (2) = 0.52, R(pred) (2) = 0.75) for TP and (N(training) = 58, N(test) = 15, R(cv) (2) = 0.67, R(pred) (2) = 0.88) for UTC. (ii) Molecular docking and molecular dynamics simulations experiments of the inhibitors into the binding site of CHK(1) aided the interpretation of the QSAR models and demonstrated the binding modes in the aspects of inhibitor's conformation, subsite interaction, and hydrogen bonding interactions, which indicated that a set of critical residues (Cys87, Glu91, Glu85, Ser147, Asp148, Glu17, Leu84 and Asn135) played a key role in the drug-target interactions. The obtained results in the present work will be fruitful for the design of new potent and selective inhibitors of CHK(1).     

Effects of the V82A and I54V mutations on the dynamics and ligand binding properties of HIV-1 protease

A major problem in the antiretroviral treatment of HIV-infections with protease-inhibitors is the emergence of resistance, resulting from the occurrence of distinct mutations within the protease molecule. In the present work we investigated the structural properties of a triple mutant (I54V-V82A-L90M) and a double mutant (V82A-L90M) that both confer strong resistance to ritonavir (RTV), but not to amprenavir (APV). For the unliganded double mutant protease molecular dynamics simulations revealed a contraction of the ligand binding pocket, which is enhanced by the I54V mutation. The observed displacement of backbone atoms of the 80s loops (residues 80–85 and 80’–85’ of the dimer) was found to primarily affect binding of the larger RTV molecule. The pocket contraction detected for the unbound protease upon mutation is also observed in the presence of APV, but not of RTV. As a consequence, the protein-ligand contacts lost upon the V82A mutation are restored by 80s loop motions for the APV-bound, but not for the RTV-bound form. RTV binding is therefore both hampered in the initial recognition step due to the poor fit of the bulky inhibitor into the small pocket of the mutant free protease and by the loss of protein-ligand interactions in the RTV-bound protease. The synergistic nature of both effects offers an explanation for the high level of resistance observed. These findings demonstrate that large inhibitors, which tightly bind to wild-type protease, may nevertheless be prone to the emergence of resistance in the presence of particular patterns of mutations. This information should be helpful for the design of novel and more effective drugs, e.g., by targeting different residues or by developing allosteric inhibitors that are capable of regulating protease dynamics.