VX680 binding in Aurora A: π-π interactions involving the conserved aromatic amino acid of the flexible glycine-rich loop

The regulation of protein kinases requires flexibility, especially near the ATP binding site. The cancer drug target Aurora A is inhibited by the ATP site inhibitor VX680, and published crystal structures show two distinct conformations. In one, a refolded glycine-rich loop creates a stacked π-π interaction between the conserved aromatic residue of the glycine-rich loop hairpin turn (F144) and the inhibitor. This refolding, associated with binding to a peptide derived from the cofactor TPX2, is absent in the other structure. We use surface plasmon resonance to measure VX680 binding to native and mutant F144A Aurora A kinase domains, with and without the TPX2 peptide. Results show that the F144 aromatic side chain contributes 2 kcal/mol to the VX680 binding energy, independent of the TPX2 peptide. This indicates that distinct VX680 bound conformations of Aurora A cannot be simply correlated with TPX2 binding and that Aurora A retains flexibility when inhibitor-bound. Molecular dynamics simulations show that alternate geometries for the π-π interactions are feasible in the absence of the rigidifying packing interactions seen in the crystal lattice.     

Dynamic Equilibrium of the Aurora A Kinase Activation Loop Revealed by Single‐Molecule Spectroscopy

The conformation of the activation loop (T‐loop) of protein kinases underlies enzymatic activity and influences the binding of small‐molecule inhibitors. By using single‐molecule fluorescence spectroscopy, we have determined that phosphorylated Aurora A kinase is in dynamic equilibrium between a DFG‐in‐like active T‐loop conformation and a DFG‐out‐like inactive conformation, and have measured the rate constants of interconversion. Addition of the Aurora A activating protein TPX2 shifts the equilibrium towards an active T‐loop conformation whereas addition of the inhibitors MLN8054 and CD532 favors an inactive T‐loop. We show that Aurora A binds TPX2 and MLN8054 simultaneously and provide a new model for kinase conformational behavior. Our approach will enable conformation‐specific effects to be integrated into inhibitor discovery across the kinome, and we outline some immediate consequences for structure‐based drug discovery.     

In silico analysis of nonsynonymous single nucleotide polymorphisms of the human adiponectin receptor 2 (ADIPOR2) gene

Polymorphisms of the ADIPOR2 gene are frequently linked to a higher risk of developing diseases including obesity, type 2 diabetes and cardiovascular diseases. Though mutations of the ADIPOR2 gene are detrimental, there is a lack of comprehensive in silico analyses of the functional and structural impacts at the protein level. Considering the involvement of ADIPOR2 in glucose uptake and fatty acid oxidation, an in silico functional analysis was conducted to explore the possible association between genetic mutations and phenotypic variations. A genomic analysis of 82 nonsynonymous SNPs in ADIPOR2 was initiated using SIFT followed by the SNAP2, nsSNPAnalyzer, PolyPhen-2, SNPs&GO, FATHMM and PROVEAN servers. A total of 10 mutations (R126W, L160Q, L195P, F201S, L235R, L235P, L256R, Y328H, E334K and Q349H) were predicted to have deleterious effects on the ADIPOR2 protein and were therefore selected for further analysis. Theoretical models of the variants were generated by comparative modeling via MODELLER 9.16. A protein structural analysis of these amino acid variants was performed using SNPeffect, I-Mutant, ConSurf, Swiss-PDB Viewer and NetSurfP to explore their solvent accessibility, molecular dynamics and energy minimization calculations. In addition, FTSite was used to predict the ligand binding sites, while NetGlycate, NetPhos2.0, UbPerd and SUMOplot were used to predict post-translational modification sites. All of the variants showed increased free energy, though F201S exhibited the highest energy increase. The root mean square deviation values of the modeled mutants strongly indicated likely pathogenicity. Remarkably, three binding sites were detected on ADIPOR2, and two mutations at positions 328 and 201 were found in the first and second binding pockets, respectively. Interestingly, no mutations were found at the post-translational modification sites. These genetic variants can provide a better understanding of the wide range of disease susceptibility associated with ADIPOR2 and aid the development of new molecular diagnostic markers for these diseases. The findings may also facilitate the development of novel therapeutic elements for associated diseases.     

Molecular insight into the interaction between IFABP and PA by using MM-PBSA and alanine scanning methods

The molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method combined with alanine-scanning mutagenesis is a very important tool for rational drug design. In this study, molecular dynamics (MD) and MM-PBSA were applied to calculate the binding free energy between the rat intestinal fatty acid binding protein (IFABP) and palmitic acid (PA) to gain insight to the interaction details. Equally spaced snapshots along the trajectory were chosen to perform the binding free energy calculation, which yields a result highly consistent with experimental value with a deviation of 0.4 kcal/mol. Computational alanine scanning was performed on the same set of snapshots by mutating the residues in IFABP to alanine and recomputing the DeltaDeltaG(binding). By postprocessing a single trajectory of the wild-type complex, the average unsigned error of our calculated DeltaDeltaG(binding) is below 1.5 kcal/mol for most of the alanine mutations of the noncharged residues (67% in total). To further investigate some particular mutants, three additional dynamical simulations of IFABP Arg126Ala, Arg106Ala, and Arg106Gln mutants were conducted. Recalculated binding free energies are well consistent with the experimental data. Moreover, the ambiguous role of Arg106 caused by the free energy change of the opposite sign when it is mutated to alanine and glutamine respectively is clarified both structurally and energetically. Typically, this can be attributed to the partial electrostatic compensation mainly from Arg56 and the obvious entropy gain in Arg106Ala mutant while not in Arg106Gln mutant. The presented structural model of IFABP-PA complex could be used to guide future studies.     

An estimation method of binding free energy in terms of ABEEMσπ/MM and continuum electrostatics fused into LIE method

A method is proposed for the estimation of absolute binding free energy of interaction between proteins and ligands. The linear interaction energy method is combined with atom‐bond electronegativity equalization method at σπ level Force field (fused into molecular mechanics) and generalized Born continuum model calculation of electrostatic solvation for the estimation of the absolute free energy of binding. The parameters of this method are calibrated by using a training set of 24 HIV‐1 protease–inhibitor complexes (PDB entry 1AAQ). A correlation coefficient of 0.93 was obtained with a root mean square deviation of 0.70 kcal mol^?1. This approach is further tested on seven inhibitor and protease complexes, and it provides small root mean square deviation between the calculated binding free energy and experimental binding free energy without reparametrization. By comparing the radii of gyration and the hydrogen bond distances between ligand and protein of three training model molecules, the consistent comparison result of binding free energy is obtained. It proves that this method of calculating the binding free energy with appropriate structural analysis can be applied to quickly assess new inhibitors of HIV‐1 proteases. To test whether the parameters of this method can apply to other drug targets, we have validated this method for the drug target cyclooxygenase‐2. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2011     

Cytochrome P450CAM enzymatic catalysis cycle: a quantum mechanics/molecular mechanics study

The catalytic pathway of cytochrome P450cam is studied by means of a hybrid quantum mechanics/molecular mechanics method. Our results reveal an active role of the enzyme in the different catalytic steps. The protein initially controls the energy gap between the high- and low-spin states in the substrate binding process, allowing thermodynamic reduction by putidaredoxin reductase and molecular oxygen addition. A second electron reduction activates the delivery of protons to the active site through a selective interaction of Thr252 and the distal oxygen causing the O--O cleavage. Finally, the protein environment catalyzes the substrate hydrogen atom abstraction step with a remarkably low free energy barrier ( approximately 8 kcal/mol). Our results are consistent with the effect of mutations on the enzymatic efficacy and provide a satisfactory explanation for the experimental failure to trap the proposed catalytically competent species, a ferryl Fe(IV) heme.     

The molecular mechanism studies of chirality effect of PHA-739358 on Aurora kinase A by molecular dynamics simulation and free energy calculations

Aurora kinase family is one of the emerging targets in oncology drug discovery and several small molecules targeting aurora kinases have been discovered and evaluated under early phase I/II trials. Among them, PHA-739358 (compound 1r) is a 3-aminopyrazole derivative with strong activity against Aurora A under early phase II trial. Inhibitory potency of compound 1r (the benzylic substituent at the pro-R position) is 30 times over that of compound 1s (the benzylic substituent at the pro-S position). In present study, the mechanism of how different configurations influence the binding affinity was investigated using molecular dynamics (MD) simulations, free energy calculations and free energy decomposition analysis. The predicted binding free energies of these two complexes are consistent with the experimental data. The analysis of the individual energy terms indicates that although the van der Waals contribution is important for distinguishing the binding affinities of these two inhibitors, the electrostatic contribution plays a more crucial role in that. Moreover, it is observed that different configurations of the benzylic substituent could form different binding patterns with protein, thus leading to variant inhibitory potency of compounds 1r and 1s. The combination of different molecular modeling techniques is an efficient way to interpret the chirality effects of inhibitors and our work gives valuable information for the chiral drug design in the near future.     

Steered Molecular Dynamics for Investigating the Interactions Between Insulin Receptor Tyrosine Kinase (IRK) and Variants of Protein Tyrosine Phosphatase 1B (PTP1B)

The aim of this study is to use steered molecular dynamics to investigate the dissociation process between IRK and PTP1Bs for wild type and five mutants (consisting of p.D181E, p.D181A, p.Q262A, p.D181A-Y46F, and p.D181A-Q262A). The gained results are observed not only the unbinding mechanism of IRK-PTP1B complexes came from pulling force profile, number of hydrogen bonds, and interaction energy between IRK and PTP1Bs but also described PTP1B’s point mutations could variably change its binding affinity towards IRK. Additionally, the binding free energy calculated by Molecular Mechanics/Poisson-Boltzmann Surface Area (MM-PBSA) is also revealed that electrostatic energy and polar solvation energy mainly made up the binding free energy of PTP1B-IRK complexes.     

乙酰胆碱酯酶AChE与1,7-二氮杂咔唑抑制剂的作用机理的分子动力学模拟

通过分子对接、分子动力学(MD)模拟以及成键自由能分析方法,从原子水平上模拟研究了3种1,7-二氮杂咔唑衍生物(分别记为M1、M2和M3)与ACh E的结合模式及相互作用机理,分析和讨论了研究体系的静电相互作用和范德华相互作用(vd W)。用MM-PBSA方法计算的3种抑制剂与ACh E之间的结合自由能与抑制剂的实验生物活性数据(IC50值)相对应。分析结果表明,残基S286与抑制剂之间形成的氢键作用有利于抑制剂与ACh E之间的结合。范德华相互作用,尤其是抑制剂与关键残基W279和Y334的作用,对抑制剂与ACh E之间的结合自由能有较大的贡献,在区分抑制剂M1(或M2)和M3的生物活性上发挥着重要的作用。     

乙酰胆碱酯酶AChE与1,7-二氮杂咔唑抑制剂的作用机理的分子动力学模拟

通过分子对接、分子动力学（MD）模拟以及成键自由能分析方法,从原子水平上模拟研究了3种1,7-二氮杂咔唑衍生物（分别记为M1、M2和M3）与AChE的结合模式及相互作用机理,分析和讨论了研究体系的静电相互作用和范德华相互作用（vdW）。用MM-PBSA方法计算的3种抑制剂与AChE之间的结合自由能与抑制剂的实验生物活性数据（IC₅₀值）相对应。分析结果表明,残基S286与抑制剂之间形成的氢键作用有利于抑制剂与AChE之间的结合。范德华相互作用,尤其是抑制剂与关键残基W279和Y334的作用,对抑制剂与AChE之间的结合自由能有较大的贡献,在区分抑制剂M1（或M2）和M3的生物活性上发挥着重要的作用。     

Computational design of protein–ligand binding: Modifying the specificity of asparaginyl‐tRNA synthetase

A method for computational design of protein–ligand interactions is implemented and tested on the asparaginyl‐ and aspartyl‐tRNA synthetase enzymes (AsnRS, AspRS). The substrate specificity of these enzymes is crucial for the accurate translation of the genetic code. The method relies on a molecular mechanics energy function and a simple, continuum electrostatic, implicit solvent model. As test calculations, we first compute AspRS‐substrate binding free energy changes due to nine point mutations, for which experimental data are available; we also perform large‐scale redesign of the entire active site of each enzyme (40 amino acids) and compare to experimental sequences. We then apply the method to engineer an increased binding of aspartyl‐adenylate (AspAMP) into AsnRS. Mutants are obtained using several directed evolution protocols, where four or five amino acid positions in the active site are randomized. Promising mutants are subjected to molecular dynamics simulations; Poisson‐Boltzmann calculations provide an estimate of the corresponding, AspAMP, binding free energy changes, relative to the native AsnRS. Several of the mutants are predicted to have an inverted binding specificity, preferring to bind AspAMP rather than the natural substrate, AsnAMP. The computed binding affinities are significantly weaker than the native, AsnRS:AsnAMP affinity, and in most cases, the active site structure is significantly changed, compared to the native complex. This almost certainly precludes catalytic activity. One of the designed sequences has a higher affinity and more native‐like structure and may represent a valid candidate for Asp activity. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Fragment‐based quantum mechanical calculation of protein–protein binding affinities

The electrostatically embedded generalized molecular fractionation with conjugate caps (EE‐GMFCC) method has been successfully utilized for efficient linear‐scaling quantum mechanical (QM) calculation of protein energies. In this work, we applied the EE‐GMFCC method for calculation of binding affinity of Endonuclease colicin–immunity protein complex. The binding free energy changes between the wild‐type and mutants of the complex calculated by EE‐GMFCC are in good agreement with experimental results. The correlation coefficient (R) between the predicted binding energy changes and experimental values is 0.906 at the B3LYP/6‐31G*‐D level, based on the snapshot whose binding affinity is closest to the average result from the molecular mechanics/Poisson–Boltzmann surface area (MM/PBSA) calculation. The inclusion of the QM effects is important for accurate prediction of protein–protein binding affinities. Moreover, the self‐consistent calculation of PB solvation energy is required for accurate calculations of protein–protein binding free energies. This study demonstrates that the EE‐GMFCC method is capable of providing reliable prediction of relative binding affinities for protein–protein complexes. © 2018 Wiley Periodicals, Inc.     

新型二氟甲基磷酸类酪氨酸蛋白磷酸酯酶1B抑制剂的分子动力学模拟和结合自由能计算

通过分子对接建立了一系列含二氟甲基磷酸基团(DFMP)或二氟甲基硫酸基团(DFMS)的抑制剂与酪氨酸蛋白磷酸酯酶1B(PTP1B)的相互作用模式, 并通过1 ns的分子动力学模拟和molecular mechanics/generalized Born surface area (MM/GBSA)方法计算了其结合自由能. 计算获得的结合自由能排序和抑制剂与靶酶间结合能力排序一致; 通过基于主方程的自由能计算方法, 获得了抑制剂与靶酶残基间相互作用的信息, 这些信息显示DFMP/DFMS基团的负电荷中心与PTP1B的221位精氨酸正电荷中心之间的静电相互作用强弱决定了此类抑制剂的活性, 进一步的分析还显示位于DFMP/DFMS基团中的氟原子或其他具有适当原子半径的氢键供体原子会增进此类抑制剂与PTP1B活性位点的结合能力.     

Ferrocene‐substituted bis(ethynyl)anthracene compounds as anticancer agents

Three compounds ( 1 – 3 ) were synthesized, where ethynylferrocene is substituted at different positions of anthracene and anthraquinone, and their biological properties were investigated. Compounds 1 – 3 were characterized using NMR and mass spectroscopies and confirmed by their single‐crystal X‐ray structure. They were also characterized from electronic and photophysical properties. All three crystal structures were optimized using density functional theory calculations. The presence of C–H???π interactions in 1 – 3 leads to the formation of two‐ and three‐dimensional networks. The bioactivity of 1 – 3 was expressed by molecular docking with various cancerous proteins, which participated in progression of cancer. Compound 2 displayed the best interaction with cancer‐related Aurora A protein in terms of both binding energy (?10.61 kcal mol^?1) as well as inhibition constant (16.74 nM). The molecular docking result also coincides with cytotoxicity on cancer cell lines (A375, HeLa) and DNA/protein binding affinity.     

β-环糊精和甾类化合物的分子动力学模拟

运用分子动力学(molecular dynamics, MD)和MM-PBSA (molecular mechanics/Poisson Boltzmann surface area)相结合的方法预测了β-环糊精(cyclodextrin, CD)和甾类客体分子包结模式. 通过重原子均方根偏差(root mean square deviation, RMSD)分析可得, 两种包结模式下客体分子都可以和β-CD形成稳定的包结. 在MD轨迹采样基础上, 采用高效MM-PBSA方法计算了两种包结模式下的包结自由能. 计算结果显示, β-CD和三个甾类客体分子包结的主要驱动力为范德华相互作用, 而溶剂化能和熵变则不利于体系的包结. 进一步分析平均构象和包结自由能发现, 对于波尼松龙, D-up (D-ring up orientation)取向为优势包结模式; 而乙炔雌二醇和雌三醇的优势包结模式均为A-up (A-ring up orientation)取向. 通过比较β-CD和三个客体分子的理论包结自由能, 预测包结稳定性的次序为乙炔雌二醇＞雌三醇＞波尼松龙, 和实验结果相一致.     

Relative importance of secondary structure and solvent accessibility to the stability of protein mutants. A case study with amino acid properties and energetics on T4 and human lysozymes

Understanding the factors influencing the stability of protein mutants is an important task in molecular and computational biology. In this work, we have approached this problem by examining the relative importance of secondary structure and solvent accessibility of the mutant residue for understanding/predicting the stability of protein mutants. We have used hydrophobic, electrostatic and hydrogen bond free energy terms and nine unique physicochemical, energetic and conformational properties of amino acids in the present study and these parameters have been related with changes in thermal stability (DeltaTm) of all the single mutants of lysozymes based on single and multiple correlation coefficients. As expected the properties reflecting hydrophobicity and hydrophobic free energy play a major role to distinguish stabilizing and destabilizing mutants. The hydrophobic free energy due to carbon and nitrogen atoms distinguish the stability of coil and strand mutations to the accuracy of 100 and 90%, respectively. In agreement with previous results, the subgroup classification based on secondary structure and the information about its location in the structure yielded good relationship with the experimental DeltaTm. We revealed that the secondary structure information is equally or more important than solvent accessibility for understanding the stability of protein mutants. The comparison of amino acid properties with free-energy terms indicate that the energetic contribution explains the mutant stability better in coil region whereas the amino acid properties do better in strand region. Further, the combination of free energies with amino acid properties increased the correlation significantly. The present study demonstrates the importance of classifying the mutants based on secondary structure to the stability of proteins upon mutations.     

Engineering the Pseudomonas aeruginosa II lectin: designing mutants with changed affinity and specificity

This article focuses on designing mutations of the PA-IIL lectin from Pseudomonas aeruginosa that lead to change in specificity. Following the previous results revealing the importance of the amino acid triad 22–23–24 (so-called specificity-binding loop), saturation in silico mutagenesis was performed, with the intent of finding mutations that increase the lectin’s affinity and modify its specificity. For that purpose, a combination of docking, molecular dynamics and binding free energy calculation was used. The combination of methods revealed mutations that changed the performance of the wild-type lectin and its mutants to their preferred partners. The mutation at position 22 resulted in 85 % in inactivation of the binding site, and the mutation at 23 did not have strong effects thanks to the side chain being pointed away from the binding site. Molecular dynamics simulations followed by binding free energy calculation were performed on mutants with promising results from docking, and also at those where the amino acid at position 24 was replaced for bulkier or longer polar chain. The key mutants were also prepared in vitro and their binding properties determined by isothermal titration calorimetry. Combination of the used methods proved to be able to predict changes in the lectin performance and helped in explaining the data observed experimentally.     

Molecular dynamic and free energy studies of primary resistance mutations in HIV-1 protease-ritonavir complexes

To understand the basis of drug resistance of the HIV-1 protease, molecular dynamic (MD) and free energy calculations of the wild-type and three primary resistance mutants, V82F, I84V, and V82F/I84V, of HIV-1 protease complexed with ritonavir were carried out. Analysis of the MD trajectories revealed overall structures of the protein and the hydrogen bonding of the catalytic residues to ritonavir were similar in all four complexes. Substantial differences were also found near the catalytic binding domain, of which the double mutant complex has the greatest impact on conformational changes of the protein and the inhibitor. The tip of the HIV-1 protease flap of the double mutant has the greater degree of opening with respect to that of the others. Additionally, the phenyl ring of Phe82 moves away from the binding pocket S1', and the conformational change of ritonavir subsite P1' consequently affects the cavity size of the protein and the conformational energy of the inhibitor. Calculations of binding free energy using the solvent continuum model were able to reproduce the same trend of the experimental inhibition constant. The results show that the resistance mutants require hydrophobic residues to maintain the interactions in the binding pocket. Changes of the cavity volume correlate well with free energy penalties due to the mutation and are responsible for the loss of drug susceptibility.     

Molecular basis of drug resistance in aurora kinases

Aurora kinases have emerged as potential targets in cancer therapy, and several drugs are currently undergoing preclinical and clinical validation. Whether clinical resistance to these drugs can arise is unclear. We exploited a hypermutagenic cancer cell line to select mutations conferring resistance to a well-studied Aurora inhibitor, ZM447439. All resistant clones contained dominant point mutations in Aurora B. Three mutations map to residues in the ATP-binding pocket that are distinct from the "gatekeeper" residue. The mutants retain wild-type catalytic activity and were resistant to all of the Aurora inhibitors tested. Our studies predict that drug-resistant Aurora B mutants are likely to arise during clinical treatment. Furthermore, because the plasticity of the ATP-binding pocket renders Aurora B insensitive to multiple inhibitors, our observations indicate that the drug-resistant Aurora B mutants should be exploited as novel drug targets.     

Molecular mechanical study of halogen bonding in drug discovery

A halogen bond is a noncovalent bond between a halogen atom (X) and a Lewis base (Y). This type of bond is attributed to the anisotropic distribution of the charge density on the halogen atom, resulting in the formation of a positive cap (called the σ-hole) centered on the A-X axis. The current research is the first reported molecular mechanical study of halogen bonding, the positive region centered on the halogen atom was represented by an extra-point (EP) of charge. The correlation between the X-EP distance and the X…Y bond length was explored to determine the optimal position of the EP. A test set of 27 halogen-containing molecules complexed to various Lewis bases was studied using molecular mechanical potentials. The molecular mechanical minimized halogen bond lengths and binding energies were in good agreement with the corresponding quantum mechanical values. The EP inclusion on the halogen atom resulted in an improvement in the accuracy of the electrostatic-potential derived charges. The solvation free energies of halobenzene molecules relative to benzene were calculated with and without EP inclusion to assess the accuracy of the developed approach. Molecular mechanical study of halo derivatives of benzotriazole complexed to cyclin-dependent protein kinase 2 (CDK2) was performed, and MM-PB(GB)SA binding energies were calculated as a case study in finding potent halogenated inhibitors that can serve as antitumor drugs.