A computational study of the protein-ligand interactions in CDK2 inhibitors: using quantum mechanics/molecular mechanics interaction energy as a predictor of the biological activity

We report a combined quantum mechanics/molecular mechanics (QM/MM) study to determine the protein-ligand interaction energy between CDK2 (cyclin-dependent kinase 2) and five inhibitors with the N(2)-substituted 6-cyclohexyl-methoxy-purine scaffold. The computational results in this work show that the QM/MM interaction energy is strongly correlated to the biological activity and can be used as a predictor, at least within a family of substrates. A detailed analysis of the protein-ligand structures obtained from molecular dynamics simulations shows specific interactions within the active site that, in some cases, have not been reported before to our knowledge. The computed interaction energy gauges the strength of protein-ligand interactions. Finally, energy decomposition and multiple regression analyses were performed to check the contribution of the electrostatic and van der Waals energies to the total interaction energy and to show the capabilities of the computational model to identify new potent inhibitors.     

An insight into the inhibitory selectivity of 4-(Pyrazol- 4-yl)-pyrimidines to CDK4 over CDK2

Revealing selectivity mechanism of cyclin-dependent kinases (CDK) and their inhibitors is an important issue to develop potential anticancer drugs. The substituted 4-(Pyrazol-4-yl)-pyrimidines are potent inhibitors of CDK4 but not of the highly homologous CDK2. In order to reveal the inhibitory selectivity of these inhibitors to CDK4 over CDK2, we select one of substituted 4-(Pyrazol-4-yl)-pyrimidines as a representative (marked as A1 hereunder) and perform molecular docking, molecular dynamics simulations and binding free energy analysis for CDK4/A1 and CDK2/A1, respectively. The electrostatic and van der Waals (vdW) interactions of the A1 inhibitor with CDK4/CDK2 are discussed. The computed binding free energies based on the MM-PBSA method are consistent with experimental bioactivity ranking of A1 inhibitor to CDK4/CDK2. On the other hand, the conformational characteristics of CDK2 and CDK4 induced by A1 inhibitor are analysed and revealed. Results demonstrate that the vdW interactions considerably contribute to binding of CDK4/CDK2 with A1 inhibitor and are similar in size. The hydrogen bonding between A1 inhibitor and CDK4/CDK2 is considerably favourable to the binding, in which the hydrogen bond between the NH group of the pyrazole group of A1 and the residue Asp158 of CDK4 plays a crucial role in inhibitory selectivity of A1 inhibitor to CDK4 over CDK2. The electrostatic interaction energy differences between the corresponding residues of CDK4/A1 and CDK2/A1 confirm the above inference. The conformational changes of CDK2 and CDK4 induced by A1 inhibitor influence the selectivity of A1 inhibitor to CDK4/CDK2.     

A comparative study of trypsin specificity based on QM/MM molecular dynamics simulation and QM/MM GBSA calculation

Hydrogen bonding and polar interactions play a key role in identification of protein-inhibitor binding specificity. Quantum mechanics/molecular mechanics molecular dynamics (QM/MM MD) simulations combined with DFT and semi-empirical Hamiltonian (AM1d, RM1, PM3, and PM6) methods were performed to study the hydrogen bonding and polar interactions of two inhibitors BEN and BEN1 with trypsin. The results show that the accuracy of treating the hydrogen bonding and polar interactions using QM/MM MD simulation of PM6 can reach the one obtained by the DFT QM/MM MD simulation. Quantum mechanics/molecular mechanics generalized Born surface area (QM/MM-GBSA) method was applied to calculate binding affinities of inhibitors to trypsin and the results suggest that the accuracy of binding affinity prediction can be significantly affected by the accurate treatment of the hydrogen bonding and polar interactions. In addition, the calculated results also reveal the binding specificity of trypsin: (1) the amidinium groups of two inhibitors generate favorable salt bridge interaction with Asp189 and form hydrogen bonding interactions with Ser190 and Gly214, (2) the phenyl of inhibitors can produce favorable van der Waals interactions with the residues His58, Cys191, Gln192, Trp211, Gly212, and Cys215. This systematic and comparative study can provide guidance for the choice of QM/MM MD methods and the designs of new potent inhibitors targeting trypsin.     

Study of the inhibition of cyclin-dependent kinases with roscovitine and indirubin-3′-oxime from molecular dynamics simulations

Molecular dynamics simulations were performed to elucidate the interactions of CDK2 and CDK5 complexes with three inhibitors: R-roscovitine, S-roscovitine, and indirubin-3′-oxime. The preference of the two complexes for R-roscovitine over the S enantiomer, as reported by the experiment, was also found by the simulations. More importantly, the simulations showed that the cause of the stronger affinity for the R enantiomer is the presence of an important hydrogen bond between R-roscovitine and the kinases not found with S-roscovitine. The simulations also showed two amino acid mutations in the active site of CDK5/R-roscovitine that favor binding-enhanced electrostatic contributions, making the inhibitor more effective for CDK5 than for CDK2. This suggests that the effectiveness of roscovitine-like inhibitors can be improved by enhancing their electrostatic interaction with the kinases. Finally, molecular mechanics–Possion–Boltzmann/surface area calculations of the CDK5/indirubin-3′-oxime system in both water-excluded and water-included environments gave significantly different electrostatic contributions to the binding. The simulations detected the displacement of a water molecule in the active site of the water-included CDK/indirubin-3′-oxime system. This resulted in a more conserved binding pattern than the water-excluded structure. Hence, in the design of new indirubin-like inhibitors, it is important to include the water molecule in the analysis. Figure Hydrogen bonding networks at the active sites of both CDK5/R-roscotivine (light grey) and CDK2/R-roscovitine (black).     

Comparison of end-point continuum-solvation methods for the calculation of protein-ligand binding free energies

We have compared the predictions of ligand‐binding affinities from several methods based on end‐point molecular dynamics simulations and continuum solvation, that is, methods related to MM/PBSA (molecular mechanics combined with Poisson–Boltzmann and surface area solvation). Two continuum‐solvation models were considered, viz., the Poisson–Boltzmann (PB) and generalised Born (GB) approaches. The nonelectrostatic energies were also obtained in two different ways, viz., either from the sum of the bonded, van der Waals, nonpolar solvation energies, and entropy terms (as in MM/PBSA), or from the scaled protein–ligand van der Waals interaction energy (as in the linear interaction energy approach, LIE). Three different approaches to calculate electrostatic energies were tested, viz., the sum of electrostatic interaction energies and polar solvation energies, obtained either from a single simulation of the complex or from three independent simulations of the complex, the free protein, and the free ligand, or the linear‐response approximation (LRA). Moreover, we investigated the effect of scaling the electrostatic interactions by an effective internal dielectric constant of the protein (?_int). All these methods were tested on the binding of seven biotin analogues to avidin and nine 3‐amidinobenzyl‐1H‐indole‐2‐carboxamide inhibitors to factor Xa. For avidin, the best results were obtained with a combination of the LIE nonelectrostatic energies with the MM+GB electrostatic energies from a single simulation, using ?_int = 4. For fXa, standard MM/GBSA, based on one simulation and using ?_int = 4–10 gave the best result. The optimum internal dielectric constant seems to be slightly higher with PB than with GB solvation. © Proteins 2012; © 2012 Wiley Periodicals, Inc.     

Synthesis, NMR and conformational studies of fucoidan fragments. VI. Fragments, content of alpha-(1--->2)-bound fucobioside unit

A series of selectively sulfated di- and trisaccharide derivatives corresponding to the potential fragments of fucoidans with a (1-->2)-alpha-bound fucobioside unit were synthesized and studied by 1H and 13C NMR spectroscopy. NOE experiments and molecular modeling were used for a conformational analysis of the compounds synthesized. In the case of disaccharides, the experimental NOE values were found to agree with those obtained using modeling with the use of density functional theory (DFT) and differ from those resulting from modeling by the molecular mechanics MM3 force field. Trisaccharide fragments partially or completely sulfated in position 4 turned out to be correctly described by both MM3 force field and DFT computation. The English version of the paper: Russian Journal of Bioorganic Chemistry, 2004, vol. 30, no. 2; see also http://www.maik.ru.     

Understanding the conformational flexibility and electrostatic properties of curcumin in the active site of rhAChE via molecular docking,molecular dynamics,and charge density analysis

Acetylcholinesterase (AChE) is an important enzyme responsible for Alzheimer’s disease, as per report, keto-enol form of curcumin inhibits this enzyme. The present study aims to understand the binding mechanism of keto-enol curcumin with the recombinant human Acetylcholinesterase (rhAChE) from its conformational flexibility, intermolecular interactions, charge density distribution, and the electrostatic properties at the active site of rhAChE. To accomplish this, a molecular docking analysis of curcumin with the rhAChE was performed, which gives the structure and conformation of curcumin in the active site of rhAChE. Further, the charge density distribution and the electrostatic properties of curcumin molecule (lifted from the active site of rhAChE) were determined from the high level density functional theory (DFT) calculations coupled with the charge density analysis. On the other hand, the curcumin molecule was optimized (gas phase) using DFT method and further, the structure and charge density analysis were also carried out. On comparing the conformation, charge density distribution and the electrostatic potential of the active site form of curcumin with the corresponding gas phase form reveals that the above said properties are significantly altered when curcumin is present in the active site of rhAChE. The conformational stability and the interaction of curcumin in the active site are also studied using molecular dynamics simulation, which shows a large variation in the conformational geometry of curcumin as well as the intermolecular interactions.     

Calculation of chromophore excited state energy shifts in response to molecular dynamics of pigment-protein complexes

The absorption and energy transfer properties of photosynthetic pigments are strongly influenced by their local environment or “site.” Local electrostatic fields vary in time with protein and chromophore molecular movement and thus transiently influence the excited state transition properties of individual chromophores. Site-specific information is experimentally inaccessible in many light-harvesting pigment–proteins due to multiple chromophores with overlapping spectra. Full quantum mechanical calculations of each chromophores excited state properties are too computationally demanding to efficiently calculate the changing excitation energies along a molecular dynamics trajectory in a pigment–protein complex. A simplified calculation of electrostatic interactions with each chromophores ground to excited state transition, the so-called charge density coupling (CDC) for site energy, CDC, has previously been developed to address this problem. We compared CDC to more rigorous quantum chemical calculations to determine its accuracy in computing excited state energy shifts and their fluctuations within a molecular dynamics simulation of the bacteriochlorophyll containing light-harvesting Fenna–Mathews–Olson (FMO) protein. In most cases CDC calculations differed from quantum mechanical (QM) calculations in predicting both excited state energy and its fluctuations. The discrepancies arose from the inability of CDC to account for the differing effects of charge on ground and excited state electron orbitals. Results of our study show that QM calculations are indispensible for site energy computations and the quantification of contributions from different parts of the system to the overall site energy shift. We suggest an extension of QM/MM methodology of site energy shift calculations capable of accounting for long-range electrostatic potential contributions from the whole system, including solvent and ions.     

Understanding microscopic binding of macrophage migration inhibitory factor with phenolic hydrazones by molecular docking, molecular dynamics simulations and free energy calculations

Macrophage migration inhibitory factor (MIF), an immunoregulatory protein, is a potential target for a number of inflammatory diseases. In the current work, the interactions between MIF and a series of phenolic hydrazones were studied by molecular docking, molecular dynamics (MD) simulations, binding free energy calculations, and binding energy decomposition analysis to determine the structural requirement for achieving favorable biological activity of phenolic hydrazones. First, molecular docking was used to predict the binding modes of inhibitors in the binding site of MIF. The good correlation between the predicted docking scores and the experimental activities shows that the binding conformations of the inhibitors in the active site of MIF are well predicted. Moreover, our results suggest that the flexibility of MIF is essential in ligand binding process. Then, MD simulations and MM/GBSA free energy calculations were employed to determine the dynamic binding process and compare the binding modes of the inhibitors with different activities. The predicted binding free energies given by MM/GBSA are not well correlated with the experimental activities for the two subsets of the inhibitors; however, for each subset, a good correlation between the predicted binding free energies and the experimental activities is achieved. The MM/GBSA free energy decomposition analysis highlights the importance of hydrophobic residues for the MIF binding of the studied inhibitors. Based on the essential factors for MIF-inhibitor interactions derived from the theoretical predictions, some derivatives were designed and the higher inhibitory activities of several candidates were confirmed by molecular docking studies. The structural insights obtained from our study are useful for designing potent inhibitors of MIF.     

Theoretical modeling and experimental studies on N-n-Decyl-2-oxo-5-nitro-1-benzylidene-methylamine

The Schiff base compound, N-n-Decyl-2-oxo-5-nitro-1-benzylidene-methylamine, has been -synthesized and characterized by IR, electronic spectroscopy, and X-ray single-crystal determination. Molecular geometry from X-ray experiment of the title compound in the ground state have been compared using the Hartree-Fock (HF) and density functional method (B3LYP) with 6-31G(d) basis set. Calculated results show that density functional theory (DFT) at B3LYP/6-31G(d) level can well reproduce the structure of the title compound. To investigate the solvent effect for the atomic charge distributions of the title compound, self-consistent reaction field theory with Onsager reaction field model was used. In addition, DFT calculations of the title compound, molecular electrostatic potential and thermodynamic properties were performed at B3LYP/6-31G(d) level of theory.     

Conformational analysis of complex oligosaccharides: the CICADA approach to the uromodulin O-glycans

Uromodulin is the pregnancy-associated Tamm-Horsfall glycoprotein, with the enhanced ability to inhibit T-cell proliferation. Pregnancy-associated structural changes mainly occur in the O-glycosylation of this glycoprotein. These include up to 12 glycan structures, made up of an unusual core type 2 sequence terminated with one, two, or three sialyl Lewis(x) sequences; this type of O-glycans could serve as E- and P-selectin ligands. The present work focuses on the most complex one; a tetradecamer made up of a type 2 core carrying three sialyl Lewis(x) branches. Five different monosaccharides are assembled by 14 glycosidic linkages. The conformational behavior of the constituting disaccharide segments was evaluated using the flexible residue procedure of the MM3 molecular mechanics procedure. For each disaccharide, the adiabatic energy surface, along with the local energy minima were established. All these results were used for the generation, prior to complete optimization of the tetradecamer. This was followed by a complete exploration of conformational hyperspace throughout the use of the single coordinate method as implemented in the CICADA program. Despite the potential flexibility of the tetradecasaccharide, only four conformational families occur, accounting for more than 95% of the total low energy conformations. For each family, the molecular properties (electrostatic, lipophilicity, and hydrogen potential) were studied. The shape of the tetradecasaccharide is best described as a flat ribbon, flanked by three branches having terminal sialyl residues. Two of the branches interact through nonbonded interactions, bringing further energy stabilization, and limiting the conformational flexibility of the sialyl residues. Only one branch maintains the original conformational features of sialyl Lewis(x). This O-glycan can be seen as a fascinating example of 'dendrimeric' structure, where the spatial arrangement of three S-Le(x) epitopes may favor its complementary 'presentations' for the interactions with E- and P-selectins.     

Palbociclib (PD 0332991) Interaction with Kinases. Theoretical and Comparative Molecular Docking Study

The optimized geometry of palbociclib, (PD 0332991) (8-cyclopentyl-6-ethanoyl-5-methyl-2-(5-(piperazin-1-yl)pyridin-2-ylamino)pyrido[2,3-d]pyrimidin-7(8H)-one), electrostatic potential map, molecular orbitals were calculated using the density functional theory. The geometry was used in a molecular docking study of palbociclib-kinase complexes, results could be explained by the charge of the nitrogen and oxygen atoms within the palbociclib. Energy gap of HOMO-LUMO surfaces, could help to explain the reactivity of the ligand and the hydrogen bonding with three different kinases, two of CDK6 and one of CDK4 type. Docking results are similar and complementary with literature reports using molecular dynamics, were hydrogen bonding was obtained and analyzed. The promiscuity of three kinases with palbociclib was detected by the docking results, thus, palbociclib could be used in other types of cancer besides myeloid leukemia. Some similarities are found with CDK4/CDK6 kinases which allow us to determine that palbociclib could be used to control other resistant inhibitor types of cancer.     

Inhibitor binding to active and inactive CDK2: the crystal structure of CDK2-cyclin A/indirubin-5-sulphonate

BACKGROUND: Cyclin-dependent kinase 2 (CDK2) is an important target for structure-based design of antitumor agents. Monomeric CDK2 is inactive. Activation requires rearrangements to key structural elements of the enzyme's active site, which accompany cyclin binding and phosphorylation. To assess the validity of using monomeric CDK2 as a model for the active kinase in structure-based drug design, we have solved the structure of the inhibitor indirubin-5-sulphonate (E226) complexed with phospho-CDK2-cyclin A and compared it with the structure of E226 bound to inactive, monomeric CDK2. RESULTS: Activation of monomeric CDK2 leads to a rotation of its N-terminal domain relative to the C-terminal lobe. The accompanying change in position of E226 follows that of the N-terminal domain, and its interactions with residues forming part of the adenine binding pocket are conserved. The environment of the ATP-ribose site, not explored by E226, is significantly different in the binary complex compared to the monomeric complex due to movement of the glycine loop. Conformational changes also result in subtle differences in hydrogen bonding and electrostatic interactions between E226's sulphonate and CDK2's phosphate binding site. Affinities calculated by LUDI for the interaction of E226 with active or inactive CDK2 differ by a factor of approximately ten. CONCLUSIONS: The accuracy of monomeric CDK2 as an inhibitor design template is restricted to the adenine binding site. The general flexibility observed for the glycine loop and subtle changes to the phosphate binding site suggest a need to study interactions between inhibitors and active CDK2 in structure-based drug design programs.     

Conformational preferences of 1,4,7-trithiacyclononane: a molecular mechanics and density functional theory study

Conformational preferences of 1,4,7-trithiacyclononane were studied using a highly efficient sampling technique based on local nonstochastic deformations and the MM2(91) force field. The results show that conformers that the molecule adopts in the crystal state were found to be low-energy conformers (LECs) within 5 kcal mol(-1) of the global minimum. A conformation with C1 symmetry was the global minimum and the C3 and C2 conformations were calculated to be 0.03 and 1.78 kcal mol(-1) higher in energy, respectively. The structures were further minimized using Density Functional Theory (DFT) calculations with two different functionals. The C2 and the C1 conformations were found to be LECs with the C3 conformation more than 4.0 kcal mol(-1) above the global minimum. The relative energies and structural ordering obtained using the BP86 functional are in agreement with the previously reported relative energies calculated using second-order Moller-Plesset (MP2) ab initio calculations. With the energy ordering being dependent on the molecular mechanics force field used, the approach of MM-->DFT (searching exhaustively the available conformational space at the MM level followed by generating the energy ordering through DFT calculations) appears to be appropriate for thiacrown ethers.     

De novo design of quinazoline derivatives as CDK2 inhibitors: 3D-QSAR,molecular fragment replacement and Volsurf predictions

Cyclin-dependent kinase 2 (CDK2) has appeared as an important drug target over the years with a multitude of therapeutic potentials. To design compounds with enhanced inhibitory potencies against CDK2, 3D-QSAR and molecular fragment replacement studies were performed on the pyrazolo[4,3-h]quinazoline derivatives, a class of potent CDK2 inhibitors. The contours of 3D-QSAR model revealed important structural features of the inhibitors related to the active site of CDK2. Based on the pyrazolo[4,3-h]quinazoline core, the different substituents at three important points were replaced with diverse molecular fragments. The compounds resulting from fragments assembly with pyrazolo[4,3-h]quinazoline core were then scored with the robust 3D-QSAR model. Furthermore, the absorption, distribution, metabolism and excretion properties of these compounds were predicted by Volsurf to eliminate inappropriate compounds. Thirty-one new potential compounds were finally obtained. These results initiated us to further optimise and design new potential inhibitors.     

Reducing CDK4/6-p16(INK4a) interface. Computational alanine scanning of a peptide bound to CDK6 protein

The tumor suppressor gene p16INK4a is commonly found altered in numerous and different types of cancer. The encoded protein arrests cell cycle in G1 phase by binding to CDK4 and CDK6, inhibiting their kinase function. In 1995, a 20-residue peptide, extracted from p16INK4a protein sequence, was discovered that retains the cell cycle inhibition properties of the endogenous tumor suppressor. However, its structure has not been determined yet. In this article, the features of a theoretical structure of the peptide bound to CDK6 are reported. The complex was modeled from CDK6-p16INK4a X-ray crystal structure and through molecular dynamics. Final structure was assessed by comparing computed binding free energy changes, when single-alanine substitutions were brought about on the peptide, to experimental data. Better concordance was obtained when including a high level of solvation effects. Solute-solvent vdW energy and electrostatic energy between solute and first shells of water, computed through a force field and considering explicit waters, were also to be included to achieve reasonably good concordance between theoretical and experimental data.