Small molecule recognition of c-Src via the Imatinib-binding conformation

The cancer drug, Imatinib, is a selective Abl kinase inhibitor that does not inhibit the closely related kinase c-Src. This one drug and its ability to selectively inhibit Abl over c-Src has been a guiding principle in virtually all kinase drug discovery efforts in the last 15 years. A prominent hypothesis explaining the selectivity of Imatinib is that Abl has an intrinsic ability to adopt an inactive conformation (termed DFG-out), whereas c-Src appears to pay a high intrinsic energetic penalty for adopting this conformation, effectively excluding Imatinib from its ATP pocket. This explanation of the difference in binding affinity of Imatinib for Abl versus c-Src makes the striking prediction that it would not be possible to design an inhibitor that binds to the DFG-out conformation of c-Src with high affinity. We report the discovery of a series of such inhibitors. We use structure-activity relationships and X-ray crystallography to confirm our findings. These studies suggest that small molecules are capable of inducing the generally unfavorable DFG-out conformation in c-Src. Structural comparison between c-Src in complex with these inhibitors allows us to speculate on the differential selectivity of Imatinib for c-Src and Abl.     

Insights into the conformational switching mechanism of the human vascular endothelial growth factor receptor type 2 kinase domain

Human vascular endothelial growth factor receptor type 2 (h-VEFGR2) is a receptor tyrosine kinase involved in the angiogenesis process and regarded as an interesting target for the design of anticancer drugs. Its activation/inactivation mechanism is related to conformational changes in its cytoplasmatic kinase domain, involving first among all the αC-helix in N-lobe and the A-loop in C-lobe. Affinity of inhibitors for the active or inactive kinase form could dictate the open or closed conformation of the A-loop, thus making the different conformations of the kinase domain receptor (KDR) domain different drug targets in drug discovery. In this view, a detailed knowledge of the conformational landscape of KDR domain is of central relevance to rationalize the efficiency and selectivity of kinase inhibitors. Here, molecular dynamics simulations were used to gain insight into the conformational switching activity of the KDR domain and to identify intermediate conformations between the two limiting active and inactive conformations. Specific energy barriers have been selectively removed to induce, and hence highlight at the atomistic level, the regulation mechanism of the A-loop opening. The proposed strategy allowed to repeatedly observe the escape of the KDR domain from the DFG-out free energy basin and to identify rare intermediate conformations between the DFG-out and the DFG-in structures to be employed in a structure-based drug discovery process.     

Molecular dynamics simulation and free energy calculation studies of the binding mechanism of allosteric inhibitors with p38α MAP kinase

p38 MAP kinase is a promising target for anti-inflammatory treatment. The classical kinase inhibitors imatinib and sorafenib as well as BI-1 and BIRB-796 were reported to bind in the DFG-out form of human p38α, known as type II or allosteric kinase inhibitors. Although DFG-out conformation has attracted great interest in the design of type II kinase inhibitors, the structural requirements for binding and mechanism of stabilization of DFG-out conformation remain unclear. As allosteric inhibition is important to the selectivity of kinase inhibitor, herein the binding modes of imatinib, sorafenib, BI-1 and BIRB-796 to p38α were investigated by molecular dynamics simulation. Binding free energies were calculated by molecular mechanics/Poisson-Boltzmann surface area method. The predicted binding affinities can give a good explanation of the activity difference of the studied inhibitors. Furthermore, binding free energies decomposition analysis and further structural analysis indicate that the dominating effect of van der Waals interaction drives the binding process, and key residues, such as Lys53, Gly71, Leu75, Ile84, Thr106, Met109, Leu167, Asp168, and Phe169, play important roles by forming hydrogen bond, salt bridge, and hydrophobic interactions with the DFG-out conformation of p38α. Finally, we also conducted a detailed analysis of BI-1, imatinib, and sorafenib binding to p38α in comparison with BIRB-796 exploited for gaining potency as well as selectivity of p38 inhibitors. These results are expected to be useful for future rational design of novel type II p38 inhibitors.     

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.     

A molecular mechanics model for imatinib and imatinib:kinase binding

Imatinib is an important anticancer drug, which binds specifically to the Abl kinase and blocks its signalling activity. To model imatinib:protein interactions, we have developed a molecular mechanics force field for imatinib and four close analogues, which is consistent with the CHARMM force field for proteins and nucleic acids. Atomic charges and Lennard‐Jones parameters were derived from a supermolecule ab initio approach. We considered the ab initio energies and geometries of a probe water molecule interacting with imatinib fragments at 32 different positions. We considered both a neutral and a protonated imatinib. The final RMS deviation between the ab initio and force field energies, averaged over both forms, was 0.2 kcal/mol. The model also reproduces the ab initio geometry and flexibility of imatinib. To apply the force field to imatinib:Abl simulations, it is also necessary to determine the most likely imatinib protonation state when it binds to Abl. This was done using molecular dynamics free energy simulations, where imatinib is reversibly protonated during a series of MD simulations, both in solution and in complex with Abl. The simulations indicate that imatinib binds to Abl in its protonated, positively‐charged form. To help test the force field and the protonation prediction, we did MD free energy simulations that compare the Abl binding affinities of two imatinib analogs, obtaining good agreement with experiment. Finally, two new imatinib variants were considered, one of which is predicted to have improved Abl binding. This variant could be of interest as a potential drug. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Structure-guided optimization of small molecule c-Abl activators

c-Abl kinase is maintained in its normal inactive state in the cell through an assembled, compact conformation. We describe two chemical series that bind to the myristoyl site of the c-Abl kinase domain and stimulate c-Abl activation. We hypothesize that these molecules activate c-Abl either by blocking the C-terminal helix from adopting a bent conformation that is critical for the formation of the autoinhibited conformation or by simply providing no stabilizing interactions to the bent conformation of this helix. Structure-based molecular modeling guided the optimization of binding and activation of c-Abl of these two chemical series and led to the discovery of c-Abl activators with nanomolar potency. The small molecule c-Abl activators reported herein could be used as molecular tools to investigate the biological functions of c-Abl and therapeutic implications of its activation.     

Molecular Mechanics Investigation of an Adenine-Adenine Non-Canonical Pair Conformational Change

Conformational changes are important in RNA for binding and catalysis and understanding these changes is important for understanding how RNA functions. Computational techniques using all-atom molecular models can be used to characterize conformational changes in RNA. These techniques are applied to an RNA conformational change involving a single base pair within a nine base pair RNA duplex. The Adenine-Adenine (AA) non-canonical pair in the sequence 5'GGUGAAGGCU3' paired with 3'PCCGAAGCCG5', where P is Purine, undergoes conformational exchange between two conformations on the timescale of tens of microseconds, as demonstrated in a previous NMR solution structure [Chen, G., et al., Biochemistry, 2006. 45: 6889-903]. The more populated, major, conformation was estimated to be 0.5 to 1.3 kcal/mol more stable at 30 °C than the less populated, minor, conformation. Both conformations are trans-Hoogsteen/sugar edge pairs, where the interacting edges on the adenines change with the conformational change. Targeted Molecular Dynamics (TMD) and Nudged Elastic Band (NEB) were used to model the pathway between the major and minor conformations using the AMBER software package. The adenines were predicted to change conformation via intermediates in which they are stacked as opposed to hydrogen-bonded. The predicted pathways can be described by an improper dihedral angle reaction coordinate. Umbrella sampling along the reaction coordinate was performed to model the free energy profile for the conformational change using a total of 1800 ns of sampling. Although the barrier height between the major and minor conformations was reasonable, the free energy difference between the major and minor conformations was the opposite of that expected based on the NMR experiments. Variations in the force field applied did not improve the misrepresentation of the free energies of the major and minor conformations. As an alternative, the Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) approximation was applied to predict free energy differences between the two conformations using a total of 800 ns of sampling. MM-PBSA also incorrectly predicted the major conformation to be higher in free energy than the minor conformation.     

Temperature‐induced conformational transition of bovine serum albumin in neutral aqueous solution by reversed‐phase liquid chromatography

The temperature‐induced conformational transition of bovine serum albumin (BSA) in neutral aqueous solution was studied using intrinsic fluorescence emission spectrum, reversed‐phase liquid chromatography and sodium dodecyl sulfate–polyacrylamide gel electrophoresis, and the conformation transition thermodynamic parameters were determined in the temperature range 12–50 °C. The results showed that, in the temperature range 12–20 °C, BSA only existed in a single conformation state A, while in the temperature range 22–50 °C, it existed in two different conformation states: A and B. The percentage of conformation state A decreased while that of conformation state B increased with the increase in temperatures, and when temperature approached 50 °C conformation state B accounted for approximately 25% of all conformation states of BSA. In the conformational transition of BSA from conformation state A to conformation state B, the positive enthalpy change, entropy change and free energy changes demonstrated that the conformational transition was endothermic, nonspontaneous and mainly entropy‐driven. © 2013 The Authors. Biomedical Chromatography published by John Wiley & Sons, Ltd.     

Sequence, structure, and active site analyses of p38 MAP kinase: exploiting DFG-out conformation as a strategy to design new type II leads

A new knowledge, structure, and sequence based strategy involving the effective exploitation of the DFG-out conformation is delineated. A comprehensive analysis of the structure, sequence, cocrystals, and active sites of p38 MAP kinase crystal structures present in Protein Data Bank (PDB) and the FDA approved MAP kinase drugs has been done, and the information is used for the design of type II leads. The 98 crystal structures, 138 cocrystals, and 31 FDA drugs comprise of 7 different sequences of 2 organisms viz., Homo sapiens and Mus musculus differing in sequence length, constituting both homo- and heterochains. Multiple sequence alignment with ClustalW showed >95% sequence similarity with highly conserved domains and a high propensity for mutations in the activation loop. The bound ligands were extracted, and their interactions with DFG in and out conformations were studied. These cocrystals and FDA drugs were fragmented on the basis of their binding interactions and their affinity to ATP and allosteric sites. The fragment library thus generated contains 106 fragments with overlapping drug fragments. A blue print constituting three main parts viz., head (ATP region), linker (DFG region), and tail (allosteric region) has thus been formulated and used to design 64 type II p38 MAP kinase inhibitors. The above strategy has been employed to design potent type II p38 MAP kinase inhibitors, which are shown to be very promising.     

Precision of free energies calculated by molecular dynamics simulations of peptides in solution

Errors in free energies for molecular replacement and for conformation change of a small model peptide have been determined empirically by repeated simulations from different starting points. All calculations have been done using thermodynamic integration, in which the system's potential energy is coupled to a parameter λ, that is increased or decreased by a small amount at each step of the simulation. The effects of several factors that may alter the precision are evaluated. These factors include: the length of the simulation, the dependence of the potential energy on λ, the use of conformational restraints, and their magnitude and form. The methods used for restraint and conformational forcing are described in detail. The free energy change, calculated as the mean from several successive simulations with alternately increasing and decreasing λ, is found to be independent of the length of the simulations. As expected, longer simulations produce more precise results. The variation of the calculated free energies is found to consist of two parts, a random error and a systematic hysteresis, i.e., a dependence on the direction in which λ changes. The hysteresis varies as the inverse of the length of the simulation and the random error as the inverse square root The advantage of the use of a different (nonlinear) dependence of the attractive and repulsive parts of the nonbonded potential energy on the coupling parameter when “creating” particles in solution is found to be very large. This nonlinear coupling was found to be superior to the use of linear coupling and a nonlinear change of the coupling parameter with the simulation time. The hysteresis in conformational free energy calculations is found to increase markedly if too weak a forcing restraint is chosen. It is shown how to deconvolute the contribution of a torsional restraint from the dependence of the free energy on a torsion angle.     

The Confine-and-Release Method: Obtaining Correct Binding Free Energies in the Presence of Protein Conformational Change

Free energy calculations are increasingly being used to estimate absolute and relative binding free energies of ligands to proteins. However, computed free energies often appear to depend on the initial protein conformation, indicating incomplete sampling. This is especially true when proteins can change conformation on ligand binding, as free energies associated with these conformational changes are either ignored or assumed to be included by virtue of the sampling performed in the calculation. Here, we show that, in a model protein system (a designed binding site in T4 Lysozyme), conformational changes can make a difference of several kcal/mol in computed binding free energies, and that they are neglected in computed binding free energies if the system remains kinetically trapped in a particular metastable state on simulation timescales. We introduce a general "confine-and-release" framework for free energy calculations that accounts for these free energies of conformational change. We illustrate its use in this model system by demonstrating that an umbrella sampling protocol can obtain converged binding free energies that are independent of the starting protein structure and include these conformational change free energies.     

Discovery and characterization of a cell-permeable, small-molecule c-Abl kinase activator that binds to the myristoyl binding site

c-Abl kinase activity is regulated by a unique mechanism involving the formation of an autoinhibited conformation in which the N-terminal myristoyl group binds intramolecularly to the myristoyl binding site on the kinase domain and induces the bending of the αI helix that creates a docking surface for the SH2 domain. Here, we report a small-molecule c-Abl activator, DPH, that displays potent enzymatic and cellular activity in stimulating c-Abl activation. Structural analyses indicate that DPH binds to the myristoyl binding site and prevents the formation of the bent conformation of the αI helix through steric hindrance, a mode of action distinct from the previously identified allosteric c-Abl inhibitor, GNF-2, that also binds to the myristoyl binding site. DPH represents the first cell-permeable, small-molecule tool compound for c-Abl activation.     

Conformational elasticity theory of chain molecules

This paper develops a conformational elasticity theory of chain molecules, which is based on three key points: (i) the molecular model is the rotational isomeric state (RIS) model; (ii) the conformational distribution function of a chain molecule is described by a function of two variables, the end-to-end distance of a chain conformation and the energy of the conformation; (iii) the rule of changes in the chain conformational states during deformation is that a number of chain conformations would vanish. The ideal deformation behavior calculated by the theory shows that the change in chain conformations is physically able to make the upward curvature of the stress-strain curve at the large-scale deformation of natural rubber. With the theory, different deformation behaviors between polymers with different chemical structures can be described, the energy term of the stress in the deformations can be predicted, and for natural rubber the fraction of the energy term is around 13%, coinciding with the experi     

Conformational analysis of oligosaccharides in solution

A complete understanding of the role of carbohydrates in biological systems is to a large extent dependent on the information available about the equilibrium mixture and about the preferred conformation of the carbohydrate molecules in solution. The conformational analysis offers a tool which can determine all possible conformations which influence the solution behavior of carbohydrates. This paper attempts to survey the progress in the theoretical conformational analysis of saccharides in solution. The conformational analysis will be discussed in detail both with respect to the strategy for the investigation of conformational properties but also with regard to the quality of the method used for calculations of the energy of the isolated molecule and free energy of solvation. Finally, examples will be given to illustrate how the methods of conformational analysis can be used to estimate the solution behavior of cyclic model compounds of carbohydrates 2-methoxytetrahydropyran, monosaccharide D -glucopyranose, and two disaccharides; β-maltose and β-cellobiose.     

Intrinsic conformational preferences of C(alpha,alpha)-dibenzylglycine

The intrinsic conformational preferences of C (alpha,alpha)-dibenzylglycine, a symmetric alpha,alpha-dialkylated amino acid bearing two benzyl substituents on the alpha-carbon atom, have been determined using quantum chemical calculations at the B3LYP/6-31+G(d,p) level. A total of 46 minimum energy conformations were found for the N-acetyl- N'-methylamide derivative, even though only nine of them showed a relative energy lower than 5.0 kcal/mol. The latter involves C 7, C 5, and alpha' backbone conformations stabilized by intramolecular hydrogen bonds and/or N-H...pi interactions. Calculation of the conformational free energies in different environments (gas-phase, carbon tetrachloride, chloroform, methanol, and water solutions) indicates that four different minima (two C 5 and two C 7) are energetically accessible at room temperature in the gas phase, while in methanol and aqueous solutions one such minimum (C 5) becomes the only significant conformation. Comparison with results recently reported for C (alpha,alpha)-diphenylglycine indicates that substitution of phenyl side groups by benzyl enhances the conformational flexibility leading to (i) a reduction of the strain of the peptide backbone and (ii) alleviating the repulsive interactions between the pi electron density of the phenyl groups and the lone pairs of the carbonyl oxygen atoms.     

Treating entropy and conformational changes in implicit solvent simulations of small molecules

Implicit solvent models are increasingly popular for estimating aqueous solvation (hydration) free energies in molecular simulations and other applications. In many cases, parameters for these models are derived to reproduce experimental values for small molecule hydration free energies. Often, these hydration free energies are computed for a single solute conformation, neglecting solute conformational changes upon solvation. Here, we incorporate these effects using alchemical free energy methods. We find significant errors when hydration free energies are estimated using only a single solute conformation, even for relatively small, simple, rigid solutes. For example, we find conformational entropy (TDeltaS) changes of up to 2.3 kcal/mol upon hydration. Interestingly, these changes in conformational entropy correlate poorly (R2 = 0.03) with the number of rotatable bonds. The present study illustrates that implicit solvent modeling can be improved by eliminating the approximation that solutes are rigid.     

Structural basis for selective inhibition of Src family kinases by PP1.

BACKGROUND: Small-molecule inhibitors that can target individual kinases are powerful tools for use in signal transduction research. It is difficult to find such compounds because of the enormous number of protein kinases and the highly conserved nature of their catalytic domains. Recently, a novel, potent, Src family selective tyrosine kinase inhibitor was reported (PP1). Here, we study the structural basis for this inhibitor's specificity for Src family kinases. RESULTS: A single residue corresponding to Ile338 (v-Src numbering; Thr338 in c-Src) in Src family tyrosine kinases largely controls PP1's ability to inhibit protein kinases. Mutation of Ile338 to a larger residue such as methionine or phenylalanine in v-Src makes this inhibitor less potent. Conversely, mutation of Ile338 to alanine or glycine increases PP1's potency. PP1 can inhibit Ser/Thr kinases if the residue corresponding to Ile338 in v-Src is mutated to glycine. We have accurately predicted several non-Src family kinases that are moderately (IC(50) approximately 1 microM) inhibited by PP1, including c-Abl and the MAP kinase p38. CONCLUSIONS: Our mutagenesis studies of the ATP-binding site in both tyrosine kinases and Ser/Thr kinases explain why PP1 is a specific inhibitor of Src family tyrosine kinases. Determination of the structural basis of inhibitor specificity will aid in the design of more potent and more selective protein kinase inhibitors. The ability to desensitize a particular kinase to PP1 inhibition of residue 338 or conversely to sensitize a kinase to PP1 inhibition by mutation should provide a useful basis for chemical genetic studies of kinase signal transduction.     

Binding energies of tyrosine kinase inhibitors: Error assessment of computational methods for imatinib and nilotinib binding

The binding energies of imatinib and nilotinib to tyrosine kinase have been determined by quantum mechanical (QM) computations, and compared with literature binding energy studies using molecular mechanics (MM). The potential errors in the computational methods include these critical factors:

• •Errors in X-ray structures such as structural distortions and steric clashes give unrealistically high van der Waals energies, and erroneous binding energies.
• •MM optimization gives a very different configuration to the QM optimization for nilotinib, whereas the imatinib ion gives similar configurations
• •Solvation energies are a major component of the overall binding energy. The QM based solvent model (PCM/SMD) gives different values from those used in the implicit PBSA solvent MM models. A major error in inhibitor—kinase binding lies in the non-polar solvation terms.
• •Solvent transfer free energies and the required empirical solvent accessible surface area factors for nilotinib and imatinib ion to give the transfer free energies have been reverse calculated. These values differ from those used in the MM PBSA studies.
• •An intertwined desolvation—conformational binding selectivity process is a balance of thermodynamic desolvation and intramolecular conformational kinetic control.
• •The configurational entropies (TΔS) are minor error sources.

     

Folding of the 25 residue Abeta(12-36) peptide in TFE/water: temperature-dependent transition from a funneled free-energy landscape to a rugged one

The free-energy landscape of the Alzheimer beta-amyloid peptide Abeta(12-36) in a 40% (v/v) 2,2,2-trifluoroethanol (TFE)/water solution was determined by using multicanonical molecular dynamics simulations. Simulations using this enhanced conformational sampling technique were initiated from a random unfolded polypeptide conformation. Our simulations reliably folded the peptide to the experimental NMR structure, which consists of two linked helices. The shape of the free energy landscape for folding was found to be strongly dependent on temperature: Above 325 K, the overall shape was funnel-like, with the bottom of the funnel coinciding exactly with the NMR structure. Below 325 K, on the other hand, the landscape became increasingly rugged, with the emergence of new conformational clusters connected by low free-energy pathways. Finally, our simulations reveal that water and TFE solvate the polypeptide in different ways: The hydrogen bond formation between TFE and Abeta was enhanced with decreasing temperature, while that between water and Abeta was depressed.