Flexible protein docking refinement using pose-dependent normal mode analysis

Modeling conformational changes in protein docking calculations is challenging. To make the calculations tractable, most current docking algorithms typically treat proteins as rigid bodies and use soft scoring functions that implicitly accommodate some degree of flexibility. Alternatively, ensembles of structures generated from molecular dynamics (MD) may be cross-docked. However, such combinatorial approaches can produce many thousands or even millions of docking poses, and require fast and sensitive scoring functions to distinguish them. Here, we present a novel approach called "EigenHex," which is based on normal mode analyses (NMAs) of a simple elastic network model of protein flexibility. We initially assume that the proteins to be docked are rigid, and we begin by performing conventional soft docking using the Hex polar Fourier correlation algorithm. We then apply a pose-dependent NMA to each of the top 1000 rigid body docking solutions, and we sample and re-score multiple perturbed docking conformations generated from linear combinations of up to 20 eigenvectors using a multi-threaded particle swarm optimization algorithm. When applied to the 63 "rigid body" targets of the Protein Docking Benchmark version 2.0, our results show that sampling and re-scoring from just one to three eigenvectors gives a modest but consistent improvement for these targets. Thus, pose-dependent NMA avoids the need to sample multiple eigenvectors and it offers a promising alternative to combinatorial cross-docking.     

Flexible protein-protein docking

Predicting the structure of protein-protein complexes using docking approaches is a difficult problem whose major challenges include identifying correct solutions, and properly dealing with molecular flexibility and conformational changes. Flexibility can be addressed at several levels: implicitly, by smoothing the protein surfaces or allowing some degree of interpenetration (soft docking) or by performing multiple docking runs from various conformations (cross or ensemble docking); or explicitly, by allowing sidechain and/or backbone flexibility. Although significant improvements have been achieved in the modeling of sidechains, methods for the explicit inclusion of backbone flexibility in docking are still being developed. A few novel approaches have emerged involving collective degrees of motion, multicopy representations and multibody docking, which should allow larger conformational changes to be modeled.     

Ras binding opens c-Raf to expose the docking site for mitogen-activated protein kinase kinase

A key signalling molecule, c-Raf, is situated downstream from Ras and upstream from the mitogen-activated protein kinase kinase (MEK). We studied the mechanism underlying the signal transduction from Ras to MEK by using probes based on the principle of fluorescence resonance energy transfer. In agreement with previous models, it was found that c-Raf adopted two conformations: open active and closed inactive. Ras binding induced the c-Raf transition from closed to open conformation, which enabled c-Raf to bind to MEK. In the presence of a cytosolic Ras mutant, c-Raf bound to, but failed to phosphorylate, MEK in the cytoplasm. In contrast, the cytosolic Ras mutant significantly enhanced MEK phosphorylation by a membrane-targeted c-Raf. These results demonstrated the essential role of Ras-induced conformational change in MEK activation by c-Raf.     

Molecular docking of balanol to dynamics snapshots of protein kinase A

Even if the structure of a receptor has been determined experimentally, it may not be a conformation to which a ligand would bind when induced fit effects are significant. Molecular docking using such a receptor structure may thus fail to recognize a ligand to which the receptor can bind with reasonable affinity. Here, we examine one way to alleviate this problem by using an ensemble of receptor conformations generated from a molecular dynamics simulation for molecular docking. Two molecular dynamics simulations were conducted to generate snapshots for protein kinase A: one with the ligand bound, the other without. The ligand, balanol, was then docked to conformations of the receptors presented by these trajectories. The Lamarckian genetic algorithm in Autodock [Goodsell et al. J Mol Recognit 1996;9(1):1-5; Morris et al. J Comput Chem 1998;19(14):1639-1662] was used in the docking. Three ligand models were used: rigid, flexible, and flexible with torsional potentials. When the snapshots were taken from the molecular dynamics simulation of the protein-ligand complex, the correct docking structure could be recovered easily by the docking algorithm in all cases. This was an easier case for challenging the docking algorithm because, by using the structure of the protein in a protein-ligand complex, one essentially assumed that the protein already had a pocket to which the ligand can fit well. However, when the snapshots were taken from the ligand-free protein simulation, which is more useful for a practical application when the structure of the protein-ligand complex is not known, several clusters of structures were found. Of the 10 docking runs for each snapshot, at least one structure was close to the correctly docked structure when the flexible-ligand models were used. We found that a useful way to identify the correctly docked structure was to locate the structure that appeared most frequently as the lowest energy structure in the docking experiments to different snapshots.     

Flexible relaxation of rigid-body docking solutions

Molecular Dynamics (MD) simulations have been performed on a set of rigid-body docking poses, carried out over 25 protein-protein complexes. The results show that fully flexible relaxation increases the fraction of native contacts (NC) by up to 70% for certain docking poses. The largest increase in the fraction of NC is observed for docking poses where anchor residues are able to sample their bound conformation. For each MD simulation, structural snap-shots were clustered and the centre of each cluster used as the MD-relaxed docking pose. A comparison between two energy-based scoring schemes, the first calculated for the MD-relaxed poses, the second for energy minimized poses, shows that the former are better in ranking complexes with large hydrophobic interfaces. Furthermore, complexes with large interfaces are generally ranked well, regardless of the type of relaxation method chosen, whereas complexes with small hydrophobic interfaces remain difficult to rank. In general, the results indicate that current force-fields are able to correctly describe direct intermolecular interactions between receptor and ligand molecules. However, these force-fields still fail in cases where protein-protein complexes are stabilized by subtle energy contributions.     

A lipid-regulated docking site on vinculin for protein kinase C.

During cell spreading, binding of actin-organizing proteins to acidic phospholipids and phosphorylation are important for localization and activity of these proteins at nascent cell-matrix adhesion sites. Here, we report on a transient interaction between the lipid-dependent protein kinase Calpha and vinculin, an early component of these sites, during spreading of HeLa cells on collagen. In vitro binding of protein kinase Calpha to vinculin tail was found dependent on free calcium and acidic phospholipids but independent of a functional kinase domain. The interaction was enhanced by conditions that favor the oligomerization of vinculin. Phosphorylation by protein kinase Calpha reached 1.5 mol of phosphate/mol of vinculin tail and required the C-terminal hydrophobic hairpin, a putative phosphatidylinositol 4,5-bisphosphate-binding site. Mass spectroscopy of peptides derived from in vitro phosphorylated vinculin tail identified phosphorylation of serines 1033 and 1045. Inhibition of C-terminal phospholipid binding at the vinculin tail by mutagenesis or deletion reduced the rate of phosphorylation to < or =50%. We suggest a possible mechanism whereby phospholipid-regulated conformational changes in vinculin may lead to exposure of a docking site for protein kinase Calpha and subsequent phosphorylation of vinculin and/or vinculin interaction partners, thereby affecting the formation of cell adhesion complexes.     

Flexible ligand docking using conformational ensembles.

Molecular docking algorithms suggest possible structures for molecular complexes. They are used to model biological function and to discover potential ligands. A present challenge for docking algorithms is the treatment of molecular flexibility. Here, the rigid body program, DOCK, is modified to allow it to rapidly fit multiple conformations of ligands. Conformations of a given molecule are pre-calculated in the same frame of reference, so that each conformer shares a common rigid fragment with all other conformations. The ligand conformers are then docked together, as an ensemble, into a receptor binding site. This takes advantage of the redundancy present in differing conformers of the same molecule. The algorithm was tested using three organic ligand protein systems and two protein-protein systems. Both the bound and unbound conformations of the receptors were used. The ligand ensemble method found conformations that resembled those determined in X-ray crystal structures (RMS values typically less than 1.5 A). To test the method's usefulness for inhibitor discovery, multi-compound and multi-conformer databases were screened for compounds known to bind to dihydrofolate reductase and compounds known to bind to thymidylate synthase. In both cases, known inhibitors and substrates were identified in conformations resembling those observed experimentally. The ligand ensemble method was 100-fold faster than docking a single conformation at a time and was able to screen a database of over 34 million conformations from 117,000 molecules in one to four CPU days on a workstation.     

Molecular docking of 1H-pyrazole derivatives to receptor tyrosine kinase and protein kinase for screening potential inhibitors

Tyrosine kinase receptor and protein kinases drawn much attention for the scientific fraternity in drug discovery due to its important role in different cancer, cardiovascular diseases and other hyper-proliferative disorders. Docking studies of pyrazole derivatives with tyrosine kinase and different serine/threonine protein kinases were employed by using flexible ligand docking approach of AutoDock 4.2. Among the molecules tested for docking study, 2-(4-chlorophenyl)-5-(3-(4-chlorophenyl)-5-methyl-1- phenyl-1H-pyrazol-4-yl)-1,3,4-thiadiazole (1b), 2-(4-methoxyphenyl)-5-(3-(4-methoxyphenyl)-5-methyl-1-phenyl-1H-pyrazol-4-yl)- 1,3,4-thiadiazole (1d) and 2-(4-chlorophenyl)-5-(3-(4-chlorophenyl)-5-methyl-1-phenyl-1H-pyrazol-4-yl)-1,3,4-thiadiazole (2b) revealed minimum binding energy of -10.09, -8.57 and -10.35 kJ/mol with VEGFR-2 (2QU5), Aurora A (2W1G) and CDK2 (2VTO) protein targets, respectively. These proteins are representatives of plausible models of interactions with different anticancer agents. All the ligands were docked deeply within the binding pocket region of all the three proteins, showing reasonable hydrogen bonds. The docking study results showed that these pyrazole derivatives are potential inhibitor of all the three protein targets; and also all these docked compounds have good inhibition constant, vdW + Hbond + desolv energy with best RMSD value.     

Role of mitogen-activated protein kinase (MAPK) docking sites on Staufen2 protein in dendritic mRNA transport

Although transport and subsequent translation of dendritic mRNA play an important role in neuronal synaptic plasticity, the underlying mechanisms for modulating dendritic mRNA transport are almost completely unknown. In this study, we identified and characterized an interaction between Staufen2 and mitogen-activated protein kinase (MAPK) with co-immunoprecipitation assays. Staufen2 utilized a docking (D) site to interact with ERK1/2; deleting the D-site decreased colocalization of Staufen2 with immunoreactive ERK1/2 in the cell body regions of cultured hippocampal neurons, and it reduced the amount of Staufen2-containing RNP complexes in the distal dendrites. In addition, the deletion completely abolished the depolarization-induced increase of Staufen2-containing RNP complexes. These results suggest that the MAPK pathway could modulate dendritic mRNA transport through its interaction with Staufen2.     

Flexible docking of DNA fragments and actinocin derivatives

We have applied molecular docking methods to systems containing nucleic acids as targets and biologically active substances as ligands. The complexes of DNA fragments and actinocin derivatives with different lengths of aminoalkyl side chains were obtained by molecular docking. It was observed that actinocin derivatives could form energetically favourable complexes with DNA both as intercalators and minor groove binders. It was shown that small changes in the binding energy (～1?kcal/mol) could result in complexes with substantially different structure. The complexes of actinocin derivatives and DNA fragments were stabilized by hydrogen bonding upon intercalation and minor groove binding. It was found that the change of solvent-accessible surface area upon binding of the actinocin derivative to DNA linear increased with the growth of methylene groups' number in ligand side chains. The solvation energy change upon binding of actinocin derivatives to DNA calculated by the WSAS method was favourable in the case of small uncharged ligands and unfavourable for positively charged ligands.     

Ensemble docking of multiple protein structures: considering protein structural variations in molecular docking

One approach to incorporate protein flexibility in molecular docking is the use of an ensemble consisting of multiple protein structures. Sequentially docking each ligand into a large number of protein structures is computationally too expensive to allow large-scale database screening. It is challenging to achieve a good balance between docking accuracy and computational efficiency. In this work, we have developed a fast, novel docking algorithm utilizing multiple protein structures, referred to as ensemble docking, to account for protein structural variations. The algorithm can simultaneously dock a ligand into an ensemble of protein structures and automatically select an optimal protein structure that best fits the ligand by optimizing both ligand coordinates and the conformational variable m, where m represents the m-th structure in the protein ensemble. The docking algorithm was validated on 10 protein ensembles containing 105 crystal structures and 87 ligands in terms of binding mode and energy score predictions. A success rate of 93% was obtained with the criterion of root-mean-square deviation <2.5 A if the top five orientations for each ligand were considered, comparable to that of sequential docking in which scores for individual docking are merged into one list by re-ranking, and significantly better than that of single rigid-receptor docking (75% on average). Similar trends were also observed in binding score predictions and enrichment tests of virtual database screening. The ensemble docking algorithm is computationally efficient, with a computational time comparable to that for docking a ligand into a single protein structure. In contrast, the computational time for the sequential docking method increases linearly with the number of protein structures in the ensemble. The algorithm was further evaluated using a more realistic ensemble in which the corresponding bound protein structures of inhibitors were excluded. The results show that ensemble docking successfully predicts the binding modes of the inhibitors, and discriminates the inhibitors from a set of noninhibitors with similar chemical properties. Although multiple experimental structures were used in the present work, our algorithm can be easily applied to multiple protein conformations generated by computational methods, and helps improve the efficiency of other existing multiple protein structure(MPS)-based methods to accommodate protein flexibility.     

Cooperative interactions of protein kinase C and cAMP-dependent protein kinase systems in human promyelocytic leukemia HL60 cells

Interactions of protein kinase C (PKC) and cAMP-dependent protein kinase (PKA) systems were investigated in HL60 cells. It was found that the differentiating effects of 12-O-tetradecanoylphorbol-13-acetate (TPA) were potentiated by dibutyryl cAMP (dbcAMP) or prostaglandin E2 (PGE2). In addition, dbcAMP or PGE2 inhibited TPA-induced binding of PKC to plasma membrane, leading to decreased protein phosphorylation, and promoted subsequent redistribution of enzyme to the nuclear membrane region. The findings are consistent with the hypothesis that PKC and PKA systems regulate cooperatively the phenotypical differentiation of leukemic cells.     

Use of docking peptides to design modular substrates with high efficiency for mitogen-activated protein kinase extracellular signal-regulated kinase.

The mitogen-activated protein kinase extracellular regulated kinase (ERK) plays a key role in the regulation of cellular proliferation. Mutations in the ERK cascade occur in 30% of malignant tumors. Thus understanding how the kinase identifies its cognate substrates as well as monitoring the activity of ERK is central to cancer research and therapeutic development. ERK binds to its protein targets, both downstream substrates and upstream activators, via a binding site distinct from the catalytic site of ERK. The substrate sequences that bind, or dock, to these sites on ERK influence the efficiency of phosphorylation. For this reason, simple peptide substrates containing only phosphorylation sequences typically possess low efficiencies for ERK. Appending short docking peptides derived from full-length protein substrates and activators of ERK to a phosphorylation sequence increased the affinity of ERK for the phosphorylation sequence by as much as 200-fold while only slightly diminishing the maximal velocity of the reaction. The efficiency of the phosphorylation reaction was increased by up to 150-fold, while the specificity of the substrate for ERK was preserved. Simple modular peptide substrates, which can be easily tailored to possess high phosphorylation efficiencies, will enhance our understanding of the regulation of ERK and provide a tool for the development of new kinase assays.     

A phosphoserine-regulated docking site in the protein kinase RSK2 that recruits and activates PDK1

The 90 kDa ribosomal S6 kinase-2 (RSK2) is a growth factor-stimulated protein kinase with two kinase domains. The C-terminal kinase of RSK2 is activated by ERK-type MAP kinases, leading to autophosphorylation of RSK2 at Ser386 in a hydrophobic motif. The N-terminal kinase is activated by 3-phosphoinositide-dependent protein kinase-1 (PDK1) through phosphorylation of Ser227, and phosphorylates the substrates of RSK. Here, we identify Ser386 in the hydrophobic motif of RSK2 as a phosphorylation-dependent docking site and activator of PDK1. Treatment of cells with growth factor induced recruitment of PDK1 to the Ser386-phosphorylated hydrophobic motif and phosphorylation of RSK2 at Ser227. A RSK2-S386K mutant showed no interaction with PDK1 or phosphorylation at Ser227. Interaction with Ser386-phosphorylated RSK2 induced autophosphorylation of PDK1. Addition of a synthetic phosphoSer386 peptide (RSK2(373-396)) increased PDK1 activity 6-fold in vitro. Finally, mutants of RSK2 and MSK1, a RSK-related kinase, with increased affinity for PDK1, were constitutively active in vivo and phosphorylated histone H3. Our results suggest a novel regulatory mechanism based on phosphoserine-mediated recruitment of PDK1 to RSK2, leading to coordinated phosphorylation and activation of PDK1 and RSK2.     

Flexible docking of peptides to class I major-histocompatibility-complex receptors

We present a new method for docking flexible peptides to class I Major-Histocompatibility-Complex (MHC) receptors. Docking is performed in two steps: (a) The charged terminal peptide residues are located by randomly distributing multiple copies of each in volumes of ≈ 150 ų at either end of the binding groove, and then minimizing the system energy using a modified multiple-copy search algorithm. This is followed by (b) bonstruction of the intervening chain using the multiple-copy bond-scaling-relaxation loop closure algorithm. In both steps, the copies tend to cluster and the size of the resulting clusters is proportional to the basin of attraction of the corresponding energy well. We show that native MHC-bound peptides have broad minima and, consequently, that misfolded, low-energy peptide conformations can be eliminated by restricting consideration to groups of peptides which cluster into broad minima. The accuracy of the method is assessed by comparing the predictions with crystallographic data for three different MHC peptide systems, at various degrees of stringency: (a) the extent to which we can determine side chain function (anchor vs. T-cell epitopes); (b) the extent to which we can determine the peptide-receptor orientation; and (c) the accuracy with which we can predict atomic coordinates. We find the method correct on (a) for 19 of the 22 non-Gly positions; the failures appearing to be a consequence of omitting solvation. Predictions related to (b) are also very encouraging, with the overall orientation of the predicted peptides being very similar to the crystal conformation, when measured by the hydrogen bonding pattern between the two. The degree of success in predicting atomic coordinates varied considerably, however, from 1.4 Å for the HLA-A2 peptide to 2.7 Å for the K^b peptide. The inaccuracy of the latter appears to reflect an incomplete target function, most likely the ommission of solvation. The calculations thus define the current limits of accuracy in docking flexible peptides to Class I receptors and identify the methodological improvements that must be made for the next advance in accuracy.     

Flexible protein-flexible ligand docking with disrupted velocity simulated annealing

By docking flexible balanol to a rigid model of protein kinase A (PKA), we found that a new simulated annealing protocol termed disrupted velocity simulated annealing (DIVE-SA) outperformed the replica-exchange method and the traditional simulated annealing method in identifying the correct docking pose. In this protocol, the atomic velocities were reassigned periodically to encourage the system to sample a large conformational space. We also found that scaling potential energy surface to reduce structural transition barriers could further facilitate docking. The DIVE-SA method was then evaluated on its ability to perform flexible ligand-flexible protein docking of three ligands (balanol, a balanol analog, and ATP) to PKA. To reduce computational time and to avoid possible unphysical structural changes resulting from the use of nonoptimal force fields, a soft restrain was applied to keep the root-mean-square-deviation (RMSD) between instantaneous protein structures and a chosen reference structure small. Because the restrain was applied to the overall RMSD rather than to individual atoms, a protein could still experience relatively large conformational changes during docking. To examine the impact of applying such a restrain on docking, we constructed two semi-flexible protein models by choosing two different crystal structures as reference. Both the balanol analog and ATP were able to dock to either one of these semi-flexible protein models. On the other hand, balanol could only dock well to one of them. Further analysis indicated that the restrain on the glycine-rich loop was too strong, preventing it to adjust its structure to accommodate balanol in the binding pocket of PKA. Removing the restrain on the glycine-rich loop resulted in much better docking poses. This finding demonstrates the important role that the flexibility of the glycine-rich loop play in accepting different ligands and should profitably not be restrained in molecular docking so that more diverse ligands can be studied.     

Two adjacent docking sites in the yeast Hog1 mitogen-activated protein (MAP) kinase differentially interact with the Pbs2 MAP kinase kinase and the Ptp2 protein tyrosine phosphatase

Functional interactions between a mitogen-activated protein kinase (MAPK) and its regulators require specific docking interactions. Here, we investigated the mechanism by which the yeast osmoregulatory Hog1 MAPK specifically interacts with its activator, the MAPK kinase Pbs2, and its major inactivator, the protein phosphatase Ptp2. We found, in the N-terminal noncatalytic region of Pbs2, a specific Hog1-binding domain, termed HBD-1. We also defined two adjacent Pbs2-binding sites in Hog1, namely, the common docking (CD) domain and Pbs2-binding domain 2 (PBD-2). The PBD-2 docking site appears to be sterically blocked in the intact Hog1 molecule, but its affinity to Pbs2 is apparent in shorter fragments of Hog1. Both the CD and the PBD-2 docking sites are required for the optimal activation of Hog1 by Pbs2, and in the absence of both sites, Hog1 cannot be activated by Pbs2. These data suggest that the initial interaction of Pbs2 with the CD site might induce a conformational change in Hog1 so that the PBD-2 site becomes accessible. The CD and PBD-2 docking sites are also involved in the specific interaction between Hog1 and Ptp2 and govern the dynamic dephosphorylation of activated Hog1. Thus, the CD and the PBD-2 docking sites play critical roles in both the activation and inactivation of Hog1.     

Functional characterization of peanut serine/threonine/tyrosine protein kinase: molecular docking and inhibition kinetics with tyrosine kinase inhibitors

Serine/threonine/tyrosine (STY) protein kinase from peanut is developmentally regulated and is induced by abiotic stresses. In addition, STY protein kinase activity is regulated by tyrosine phosphorylation. Kinetic mechanism of plant dual specificity protein kinases is not studied so far. Recombinant STY protein kinase occurs as a monomer in solution as shown by gel filtration chromatography. The relative phosphorylation rate of kinase against increasing enzyme concentrations follows a first-order kinetics indicating an intramolecular phosphorylation mechanism. Moreover, the active recombinant STY protein kinase could not transphosphorylate a kinase-deficient mutant of STY protein kinase. Molecular docking studies revealed that the tyrosine kinase inhibitors bind the protein kinase at the same region as ATP. STY protein kinase activity was inhibited by the tyrosine kinase inhibitors, and the inhibitor potency series against the recombinant STY protein kinase was tyrphostin > genistein > staurosporine. The inhibition constant (K(i)), and the IC(50) value of STY protein kinase for tyrosine kinase inhibitors with ATP and histone are discussed. All the inhibitors competed with ATP. Genistein was an uncompetitive inhibitor with histone, whereas staurosporine and tyrphostin were linear mixed type noncompetitive inhibitors with histone. Molecular docking and kinetic analysis revealed that Y148F mutant of the "ATP-binding loop" and Y297F mutant of the "activation loop" showed a dramatic increase in K(i) values for genistein and tyrphostin with respect to wild-type STY protein kinase. Data presented here provide the direct evidence on the mechanism of inhibition of plant protein kinases by tyrosine kinase inhibitors. This study also suggests that tyrosine kinase inhibitors may be useful in unraveling the plant tyrosine phosphorylation signaling cascades.