Exploring free energy landscapes of large conformational changes: molecular dynamics with excited normal modes

Sampling of conformational ensemble for virtual screening using molecular dynamics simulations and normal mode analysis

Aim: Molecular dynamics simulations and normal mode analysis are well-established approaches to generate receptor conformational ensembles (RCEs) for ligand docking and virtual screening. Here, we report new fast molecular dynamics-based and normal mode analysis-based protocols combined with conformational pocket classifications to efficiently generate RCEs. Materials & Methods: We assessed our protocols on two well-characterized protein targets showing local active site flexibility, dihydrofolate reductase and large collective movements, CDK2. The performance of the RCEs was validated by distinguishing known ligands of dihydrofolate reductase and CDK2 among a dataset of diverse chemical decoys. Results & discussion: Our results show that different simulation protocols can be efficient for generation of RCEs depending on different kind of protein flexibility.[Formula: see text]

Conformational Equilibrium of CDK/Cyclin Complexes by Molecular Dynamics with Excited Normal Modes

Cyclin-dependent kinases (CDKs) and their associated regulatory cyclins are central for timely regulation of cell-cycle progression. They constitute attractive pharmacological targets for development of anticancer therapeutics, since they are frequently deregulated in human cancers and contribute to sustained, uncontrolled tumor proliferation. Characterization of their structural/dynamic features is essential to gain in-depth insight into structure-activity relationships. In addition, the identification of druggable pockets or key intermediate conformations yields potential targets for the development of novel classes of inhibitors. Structural studies of CDK2/cyclin A have provided a wealth of information concerning monomeric/heterodimeric forms of this kinase. There is, however, much less structural information for other CDK/cyclin complexes, including CDK4/cyclin D1, which displays an alternative (open) position of the cyclin partner relative to CDK, contrasting with the closed CDK2/cyclin A conformation. In this study, we carried out normal-mode analysis and enhanced sampling simulations with our recently developed method, molecular dynamics with excited normal modes, to understand the conformational equilibrium on these complexes. Interestingly, the lowest-frequency normal mode computed for each complex described the transition between the open and closed conformations. Exploration of these motions with an explicit-solvent representation using molecular dynamics with excited normal modes confirmed that the closed conformation is the most stable for the CDK2/cyclin A complex, in agreement with their experimentally available structures. On the other hand, we clearly show that an open↔closed equilibrium may exist in CDK4/cyclin D1, with closed conformations resembling that captured for CDK2/cyclin A. Such conformational preferences may result from the distinct distributions of frustrated contacts in each complex. Using the same approach, the putative roles of the Thr(160) phosphoryl group and the T-loop conformation were investigated. These results provide a dynamic view of CDKs revealing intermediate conformations not yet characterized for CDK members other than CDK2, which will be useful for the design of inhibitors targeting critical conformational transitions.