Analysis of CDK2 active-site hydration: a method to design new inhibitors

Are Waters around RNA More than Just a Solvent? - An Insight from Molecular Dynamics Simulations

Hydrating water molecules are believed to be an inherent part of the RNA structure and have a considerable impact on RNA conformation. However, the magnitude and mechanism of the interplay between water molecules and the RNA structure are still poorly understood. In principle, such hydration effects can be studied by molecular dynamics (MD) simulations. In our recent MD studies, we observed that the choice of water model has a visible impact on the predicted structure and structural dynamics of RNA and, in particular, has a larger effect than type, parametrization, and concentration of the ions. Furthermore, the water model effect is sequence dependent and modulates the sequence dependence of A-RNA helical parameters. Clearly, the sensitivity of A-RNA structural dynamics to the water model parametrization is a rather spurious effect that complicates MD studies of RNA molecules. These results nevertheless suggest that the sequence dependence of the A-RNA structure, usually attributed to base stacking, might be driven by the structural dynamics of specific hydration. Here, we present a systematic MD study that aimed to (i) clarify the atomistic mechanism of the water model sensitivity and (ii) discover whether and to what extent specific hydration modulates the A-RNA structural variability. We carried out an extended set of MD simulations of canonical A-RNA duplexes with TIP3P, TIP4P/2005, TIP5P, and SPC/E explicit water models and found that different water models provided a different extent of water bridging between 2'-OH groups across the minor groove, which in turn influences their distance and consequently also inclination, roll, and slide parameters. Minor groove hydration is also responsible for the sequence dependence of these helical parameters. Our simulations suggest that TIP5P is not optimal for RNA simulations.

Protein kinase-inhibitor database: structural variability of and inhibitor interactions with the protein kinase P-loop

Structure-based drug design of protein-kinase inhibitors has been facilitated by availability of an enormous number of structures in the Protein Databank (PDB), systematic analyses of which can provide insight into the factors that govern ligand-protein kinase interactions and into the conformational variability of the protein kinases. In this study, a nonredundant database containing 755 unique, curated, and annotated PDB protein kinase-inhibitor complexes (each consisting of a single protein kinase chain, a ligand, and water molecules around the ligand) was created. With this dataset, analyses were performed of protein conformational variability and interactions of ligands with 11 P-loop residues. Analysis of ligand-protein interactions included ligand atom preference, ligand-protein hydrogen bonds, and the number and position of crystallographic water molecules around important P-loop residues. Analysis of variability in the conformation of the P-loop considered backbone and side-chain dihedral angles, and solvent accessible surface area (SASA). A distorted conformation of the P-loop was observed for some of the protein kinase structures. Lower SASA was observed for the hydrophobic residue in beta1 of several members of the AGC family of protein kinases. Our systematic studies were performed amino acid-by-amino acid, which is unusual for analyses of protein kinase-inhibitor complexes.

Explicit Water Models Affect the Specific Solvation and Dynamics of Unfolded Peptides While the Conformational Behavior and Flexibility of Folded Peptides Remain Intact

Conventional molecular dynamics simulations on 50 ns to 1 μs time scales were used to study the effects of explicit solvent models on the conformational behavior and solvation of two oligopeptide solutes: α-helical EK-peptide (14 amino acids) and a β-hairpin chignolin (10 amino acids). The widely used AMBER force fields (ff99, ff99SB, and ff03) were combined with four of the most commonly used explicit solvent models (TIP3P, TIP4P, TIP5P, and SPC/E). Significant differences in the specific solvation of chignolin among the studied water models were identified. Chignolin was highly solvated in TIP5P, whereas reduced specific solvation was found in the TIP4P, SPC/E, and TIP3P models for kinetic, thermodynamic, and both kinetic and thermodynamic reasons, respectively. The differences in specific solvation did not influence the dynamics of structured parts of the folded peptide. However, substantial differences between TIP5P and the other models were observed in the dynamics of unfolded chignolin, stability of salt bridges, and specific solvation of the backbone carbonyls of EK-peptide. Thus, we conclude that the choice of water model may affect the dynamics of flexible parts of proteins that are solvent-exposed. On the other hand, all water models should perform similarly for well-structured folded protein regions. The merits of the TIP3P model include its high and overestimated mobility, which accelerates simulation processes and thus effectively increases sampling.

Why pyridine containing pyrido[2,3-d]pyrimidin-7-ones selectively inhibit CDK4 than CDK2: insights from molecular dynamics simulation

Designing selective cyclin-dependent kinase 4 (CDK4) inhibitors is an area of intense research to develop potential anticancer drugs. The molecular basis governing the selective inhibition of CDK4 by lig17 (6-bromo-8-cyclopentyl-2-(5-piperazin-1-yl-pyridin-2-ylamino)-8H-pyrido[2,3-d]pyrimidin-7-one) has been investigated using molecular dynamics simulation. The positive charge on the ligand was determined to be an important contributor for CDK4 selectivity due to the electronegative nature of its active site. Similar studies on CDK2 indicated that Lys89 intrudes into the active site displacing the positive charge on lig17 away from the active center. This intrusion was observed to propel a drastic conformational change in lig17, weakening its binding interactions with the protein. The pyridine nitrogen (N(AR)) of lig17 was capable of interacting with His95 (CDK4) through hydrogen bonding. N(AR) also showed a strong tendency to mediate protein-ligand interactions through a bridged water molecule, only when bound to CDK4. The G-loop of CDK4 was observed to fluctuate extensively when complexed with lig17 and a novel "flipping-out" mechanism exhibited by Tyr17(CDK4/CDK4-17) is reported in this study. Although these proteins have similar folds, the results from principal component analysis (PCA) indicate that CDK4 and CDK2 follow an anti-correlated behavior towards the accessibility of the active site.

Flexibility of human cytochromes P450: molecular dynamics reveals differences between CYPs 3A4, 2C9, and 2A6, which correlate with their substrate preferences

Molecular dynamics (MD) simulations at normal and high temperature were used to study the flexibility and malleability of three microsomal cytochromes P450 (CYPs): CYP3A4, CYP2C9, and CYP2A6. Comparison of B-factors (describing the atomic fluctuations) between X-ray and MD data shows that the X-ray B-factors are significantly lower in the regions where the crystal contacts occur than for other regions. Consequently, the conclusions about CYP flexibility based solely on the X-ray data might be misleading. Comparison of flexibility patterns of the three CYPs enabled common features and variations in flexibility and malleability of the studied CYPs to be identified. The previously described pattern of flexibility in topological elements of microsomal CYPs (a rigid heme binding core, a malleable distal side and intermediately flexible proximal side) was confirmed. These topological features provide an important combination of high stereo- and regio-specificity (mediated by the relative rigidity in the neighborhood of the heme), together with high substrate promiscuity due to the more flexible active site and the malleability of the distal side. The data acquired here show that the malleability of the three studied CYPs correlates with their substrate specificity: CYP2A6 has a narrow substrate range and is the most rigid, CYP3A4 is the most promiscuous CYP known and is the most malleable, and CYP2C9 is intermediate in terms of both its substrate specificity and malleability. Thus, the malleability of CYPs is probably a major determinant of their substrate specificity.

Ligand-protein docking with water molecules

The presence of water molecules plays an important role in the accuracy of ligand-protein docking predictions. Comprehensive docking simulations have been performed on a large set of ligand-protein complexes whose crystal structures contain water molecules in their binding sites. Only those water molecules found in the immediate vicinity of both the ligand and the protein were considered. We have investigated whether prior optimization of the orientation of water molecules in either the presence or absence of the bound ligand has any effect on the accuracy of docking predictions. We have observed a statistically significant overall increase in accuracy when water molecules are included during docking simulations and have found this to be independent of the method of optimization of the orientation of water molecules. These results confirm the importance of including water molecules whenever possible in a ligand-protein docking simulation. Our findings also reveal that prior optimization of the orientation of water molecules, in the absence of any bound ligand, does not have a detrimental effect on the improved accuracy of ligand-protein docking. This is important, given the use of docking simulations to predict the binding modes of new ligands or drug molecules.

Novel structural features of CDK inhibition revealed by an ab initio computational method combined with dynamic simulations

The rational development of specific inhibitors for the approximately 500 protein kinases encoded in the human genome is impeded by a poor understanding of the structural basis for the activity and selectivity of small molecules that compete for ATP binding. Combining classical dynamic simulations with a novel ab initio computational approach linear-scalable to molecular interactions involving thousands of atoms, we have investigated the binding of five distinct inhibitors to the cyclin-dependent kinase CDK2. We report here that polarization and dynamic hydrogen bonding effects, so far undetected by crystallography, affect both their activity and selectivity. The effects arise from the specific solvation patterns of water molecules in the ATP binding pocket or the intermittent formation of hydrogen bonds during the dynamics of CDK/inhibitor interactions and explain the unexpectedly high potency of certain inhibitors such as 3-(3H-imidazol-4-ylmethylene)-5-methoxy-1,3-dihydro-indol-2-one (SU9516). The Lys89 residue in the ATP-binding pocket of CDK2 is observed to form temporary hydrogen bonds with the three most potent inhibitors. This residue is replaced in CDK4 by Thr89, whose shorter side-chain cannot form similar bonds, explaining the relative selectivity of the inhibitors for CDK2. Our results provide a generally applicable computational method for the analysis of biomolecular structures and reveal hitherto unrecognized features of the interaction between protein kinases and their inhibitors.

Optimized and parallelized implementation of the electronegativity equalization method and the atom-bond electronegativity equalization method

The most common way to calculate charge distribution in a molecule is ab initio quantum mechanics (QM). Some faster alternatives to QM have also been developed, the so-called "equalization methods" EEM and ABEEM, which are based on DFT. We have implemented and optimized the EEM and ABEEM methods and created the EEM SOLVER and ABEEM SOLVER programs. It has been found that the most time-consuming part of equalization methods is the reduction of the matrix belonging to the equation system generated by the method. Therefore, for both methods this part was replaced by the parallel algorithm WIRS and implemented within the PVM environment. The parallelized versions of the programs EEM SOLVER and ABEEM SOLVER showed promising results, especially on a single computer with several processors (compact PVM). The implemented programs are available through the Web page http://ncbr.chemi.muni.cz/~n19n/eem_abeem.

Different Mechanisms of CDK5 and CDK2 Activation as Revealed by CDK5/p25 and CDK2/Cyclin A Dynamics

A detailed analysis is presented of the dynamics of human CDK5 in complexes with the protein activator p25 and the purine-like inhibitor roscovitine. These and other findings related to the activation of CDK5 are critically reviewed from a molecular perspective. In addition, the results obtained on the behavior of CDK5 are compared with data on CDK2 to assess the differences and similarities between the two kinases in terms of (i) roscovitine binding, (ii) regulatory subunit association, (iii) conformational changes in the T-loop following CDK/regulatory subunit complex formation, and (iv) specificity in CDK/regulatory subunit recognition. An energy decomposition analysis, used for these purposes, revealed why the binding of p25 alone is sufficient to stabilize the extended active T-loop conformation of CDK5, whereas the equivalent conformational change in CDK2 requires both the binding of cyclin A and phosphorylation of the Thr(160) residue. The interaction energy of the CDK5 T-loop with p25 is about 26 kcal.mol(-1) greater than that of the CDK2 T-loop with cyclin A. The binding pattern between CDK5 and p25 was compared with that of CDK2/cyclin A to find specific regions involved in CDK/regulatory subunit recognition. The analyses performed revealed that the alphaNT-helix of cyclin A interacts with the alpha6-alpha7 loop and the alpha7 helix of CDK2, but these regions do not interact in the CDK5/p25 complex. Further differences between the CDK5/p25 and CDK2/cyclin A systems studied are discussed with respect to their specific functionality.

The effect of a tightly bound water molecule on scaffold diversity in the computer-aided de novo ligand design of CDK2 inhibitors

We have determined the effects that tightly bound water molecules have on the de novo design of cyclin-dependent kinase-2 (CDK2) ligands. In particular, we have analyzed the impact of a specific structural water molecule on the chemical diversity and binding mode of ligands generated through a de novo structure-based ligand generation method in the binding site of CDK2. The tightly bound water molecule modifies the size and shape of the binding site and we have found that it also imposed constraints on the observed binding modes of the generated ligands. This in turn had the indirect effect of reducing the chemical diversity of the underlying molecular scaffolds that were able to bind to the enzyme satisfactorily. [Figure: see text].

Study of a ligand complexed with Cdk2/Cdk4 by computer simulation

Cyclin-dependent kinases (Cdks) play important roles in the regulation of the cell cycle. Their inhibitors have entered clinical trials to treat cancer. Very recently, Davis et al. (Nat Struct Biol 9:745-749, 2002) have found a ligand NU6102, which has a high affinity with cyclin-dependent kinase 2 (K(i) = 6 nM) but a low affinity with cyclin-dependent kinase 4 (K(i) = 1,600 nM). To understand the selectivity, we use homology modeling, molecular docking, molecular dynamics and free-energy calculations to analyze the interactions. A rational 3D model of the Cdk4-NU6102 complex is built. Asp86 is a key residue that recognizes NU6102 more effectively with Cdk2 rather than Cdk4. Good binding free energies are obtained. Energetic analysis reveals that van der Waals interaction and nonpolar contributions to solvent are favorable in the formation of complexes and the sulfonamide group of the ligand plays a crucial role for binding selectivity between Cdk2 and Cdk4.

The mechanism of inhibition of the cyclin-dependent kinase-2 as revealed by the molecular dynamics study on the complex CDK2 with the peptide substrate HHASPRK

Molecular dynamics (MD) simulations were used to explain structural details of cyclin-dependent kinase-2 (CDK2) inhibition by phosphorylation at T14 and/or Y15 located in the glycine-rich loop (G-loop). Ten-nanosecond-long simulations of fully active CDK2 in a complex with a short peptide (HHASPRK) substrate and of CDK2 inhibited by phosphorylation of T14 and/or Y15 were produced. The inhibitory phosphorylations at T14 and/or Y15 show namely an ATP misalignment and a G-loop shift (~5 A) causing the opening of the substrate binding box. The biological functions of the G-loop and GxGxxG motif evolutionary conservation in protein kinases are discussed. The position of the ATP gamma-phosphate relative to the phosphorylation site (S/T) of the peptide substrate in the active CDK2 is described and compared with inhibited forms of CDK2. The MD results clearly provide an explanation previously not known as to why a basic residue (R/K) is preferred at the P(2) position in phosphorylated S/T peptide substrates.

Loop Flexibility and Solvent Dynamics as Determinants for the Selective Inhibition of Cyclin-Dependent Kinase 4: Comparative Molecular Dynamics Simulation Studies of CDK2 and CDK4

The design and discovery of selective cyclin-dependent kinase 4 (CDK4) inhibitors have been actively pursued in order to develop therapeutic cancer treatments. By means of a consecutive computational protocol involving homology modeling, docking experiments, and molecular dynamics simulations, we examine the characteristic structural and dynamic properties that distinguish CDK4 from CDK2 in its complexation with selective inhibitors. The results for all three CDK4-selective inhibitors under investigation show that the large-amplitude motion of a disordered loop of CDK4 is damped out in the presence of the inhibitors whereas their binding in the CDK2 active site has little effect on the loop flexibility. It is also found that the binding preference of CDK4- selective inhibitors for CDK4 over CDK2 stems from the reduced solvent accessibility in the active site of the former due to the formation of a stable hydrogen-bond triad by the Asp99, Arg101, and Thr102 side chains at the top of the active-site gorge. Besides the differences in loop flexibility and solvent accessibility, the dynamic stabilities of the hydrogen bonds between the inhibitors and the side chain of the lysine residue at the bottom of the active site also correlate well with the relative binding affinities of the inhibitors for the two CDKs. These results highlight the usefulness of this computational approach in evaluating the selectivity of a CDK inhibitor, and demonstrate the necessity of considering protein flexibility and solvent effects in designing new selective CDK4-selective inhibitors.