A flexible docking scheme to explore the binding selectivity of PDZ domains

Modeling of protein binding site flexibility in molecular docking is still a challenging problem due to the large conformational space that needs sampling. Here, we propose a flexible receptor docking scheme: A dihedral restrained replica exchange molecular dynamics (REMD), where we incorporate the normal modes obtained by the Elastic Network Model (ENM) as dihedral restraints to speed up the search towards correct binding site conformations. To our knowledge, this is the first approach that uses ENM modes to bias REMD simulations towards binding induced fluctuations in docking studies. In our docking scheme, we first obtain the deformed structures of the unbound protein as initial conformations by moving along the binding fluctuation mode, and perform REMD using the ENM modes as dihedral restraints. Then, we generate an ensemble of multiple receptor conformations (MRCs) by clustering the lowest replica trajectory. Using ROSETTA LIGAND , we dock ligands to the clustered conformations to predict the binding pose and affinity. We apply this method to postsynaptic density‐95/Dlg/ZO‐1 (PDZ) domains; whose dynamics govern their binding specificity. Our approach produces the lowest energy bound complexes with an average ligand root mean square deviation of 0.36 Å. We further test our method on (i) homologs and (ii) mutant structures of PDZ where mutations alter the binding selectivity. In both cases, our approach succeeds to predict the correct pose and the affinity of binding peptides. Overall, with this approach, we generate an ensemble of MRCs that leads to predict the binding poses and specificities of a protein complex accurately.     

How to choose relevant multiple receptor conformations for virtual screening: a test case of Cdk2 and normal mode analysis

Better treatment of protein flexibility is essential in structure-based drug design projects such as virtual screening and protein-ligand docking. Diversity in ligand-binding mechanisms and receptor conformational changes makes it difficult to treat dynamic features of the receptor during the docking simulation. Thus, the use of pregenerated multiple receptor conformations is applied today in virtual screening studies. However, generation of a small relevant set of receptor conformations remains challenging. To address this problem, we propose a new protocol for the generation of multiple receptor conformations via normal mode analysis and for the selection of several receptor conformations suitable for docking/virtual screening. We validated this protocol on cyclin-dependent kinase 2, which possesses a binding site located at the interface between two subdomains and is known to undergo significant conformational changes in the active site region upon ligand binding. We believe that the suggested rules for the choice of suitable receptor conformations can be applied to other targets when dealing with in silico screening on flexible receptors.     

Optimized grid-based protein-protein docking as a global search tool followed by incorporating experimentally derivable restraints

Unbound protein docking, or the computational prediction of the structure of a protein complex from the structures of its separated components, is of importance but still challenging. A practical approach toward reliable results for unbound docking is to incorporate experimentally derived information with computation. To this end, truly systematic search of the global docking space is desirable. The fast Fourier transform (FFT) docking is a systematic search method with high computational efficiency. However, by using FFT to perform unbound docking, possible conformational changes upon binding must be treated implicitly. To better accommodate the implicit treatment of conformational flexibility, we develop a rational approach to optimize "softened" parameters for FFT docking. In connection with the increased "softness" of the parameters in this global search step, we use a revised rule to select candidate models from the search results. For complexes designated as of low and medium difficulty for unbound docking, these adaptations of the original FTDOCK program lead to substantial improvements of the global search results. Finally, we show that models resulted from FFT-based global search can be further filtered with restraints derivable from nuclear magnetic resonance (NMR) chemical shift perturbation or mutagenesis experiments, leading to a small set of models that can be feasibly refined and evaluated using computationally more expensive methods and that still include high-ranking near-native conformations.     

A hierarchical method for generating low-energy conformers of a protein-ligand complex

A novel dynamical protocol for finding the low-energy conformations of a protein-ligand complex is described. The energy functions examined consist of an empirical force field with four different dielectric screening models; the generalized Born/surface area model also is examined. Application of the method to three complexes of known crystal structure provides insights into the energy functions used for selecting low-energy docked conformations and into the structure of the binding-energy surface. Evidence is presented that the local energy minima of a ligand in a binding site are arranged in a hierarchical fashion. This observation motivates the construction of a hierarchical docking algorithm that substantially enriches the population of ligand conformations close to the crystal conformation. The algorithm is also adapted to permit docking into a flexible binding site and preliminary tests of this method are presented. Proteins 33:475–495, 1998. © 1998 Wiley-Liss, Inc.     

Induced fit in liver X receptor beta: a molecular dynamics-based investigation

Ligand induced fit phenomenon occurring at the ligand binding domain of the liver X receptor beta (LXRbeta) was investigated by means of molecular dynamics. Reliability of a 4-ns trajectory was tested from two distinct LXRbeta crystal complexes 1PQ6B/GW and 1PQ9B/T09 characterized by an open and a closed state of the pocket, respectively. Crossed complexes 1PQ6B/T09 and 1PQ9B/GW were then submitted to the same molecular dynamic conditions, which were able to recover LXRbeta conformations similar to the original crystallography data. Analysis of "open to closed" and "closed to open" conformational transitions pointed out the dynamic role of critical residues lining the ligand binding pocket involved in the local remodeling upon ligand binding (e.g., Phe271, Phe329, Phe340, Arg319, Glu281). Altogether, the present study indicates that the molecular dynamic protocol is a consistent approach for managing LXRbeta-related induced fit process. This protocol could therefore be used for refining ligand docking solutions of a structure-based design strategy.     

A genetic algorithm-based protocol for docking ensembles of small ligands using experimental restraints

Summary A genetic algorithm (GA) based method for docking ensembles of small, flexible ligands to receptor proteins using NMR-derived constraints is described. In this method, three translations and rotations of the ligand and the dihedral angles of the ligand are represented by binary strings and evolve under the genetic operators of cross-over, mutation, migration and selection. The fitness function for the selection process includes distance and dihedral restraints and a repulsive van der Waals term. The GA was applied to a three-atom model system as well as to the streptavidin-biotin complex using simulated intermolecular distance restraints. In both systems, the GA was able to obtain low-energy conformations when only a single binding site was simulated. Calculations were also performed using distance restraints from two distinct binding sites. In these simulations, the GA was able to obtain low-energy conformations corresponding to ligand molecules in each of the two sites. The inclusion of additional ligands in the ensemble did not result in an energetic benefit, confirming that only two ligand conformations were necessary to fulfill the distance restraints. This method allows for a direct investigation of the minimum number of ligand orientations necessary to fulfill experimental distance restraints, and simultaneously yields detailed structural information about each site.     

Flexible ligand docking to multiple receptor conformations: a practical alternative

State of the art docking algorithms predict an incorrect binding pose for about 50-70% of all ligands when only a single fixed receptor conformation is considered. In many more cases, lack of receptor flexibility results in meaningless ligand binding scores, even when the correct pose is obtained. Incorporating conformational rearrangements of the receptor binding pocket into predictions of both ligand binding pose and binding score is crucial for improving structure-based drug design and virtual ligand screening methodologies. However, direct modeling of protein binding site flexibility remains challenging because of the large conformational space that must be sampled, and difficulties remain in constructing a suitably accurate energy function. Here we show that using multiple fixed receptor conformations, either experimentally determined by crystallography or NMR, or computationally generated, is a practical shortcut that may improve docking calculations. In several cases, such an approach has led to experimentally validated predictions.     

Evaluation of site-directed spin labeling for characterizing protein-ligand complexes using simulated restraints

Simulation studies have been performed to evaluate the utility of site-directed spin labeling for determining the structures of protein-ligand complexes, given a known protein structure. Two protein-ligand complexes were used as model systems for these studies: a 1.9-A-resolution x-ray structure of a dihydrofolate reductase mutant complexed with methotrexate, and a 1.5-A-resolution x-ray structure of the V-Src tyrosine kinase SH2 domain complexed with a five-residue phosphopeptide. Nitroxide spin labels were modeled at five dihydrofolate reductase residue positions and at four SH2 domain residue positions. For both systems, after energy minimization, conformational ensembles of the spin-labeled residues were generated by simulated annealing while holding the remainder of the protein-ligand complex fixed. Effective distances, simulating those that could be obtained from (1)H-NMR relaxation measurements, were calculated between ligand protons and the spin labels. These were converted to restraints with several different levels of precision. Restrained simulated annealing calculations were then performed with the aim of reproducing target ligand-binding modes. The effects of incorporating a few supplementary short-range (< or =5.0 A) distance restraints were also examined. For the dihydrofolate reductase-methotrexate complex, the ligand-binding mode was reproduced reasonably well using relatively tight spin-label restraints, but methotrexate was poorly localized using loose spin-label restraints. Short-range and spin-label restraints proved to be complementary. For the SH2 domain-phosphopeptide complex without the short-range restraints, the peptide did not localize to the correct depth in the binding groove; nevertheless, the orientation and internal conformation of the peptide was reproduced moderately well. Use of the spin-label restraints in conjunction with the short-range restraints resulted in relatively well defined structural ensembles. These results indicate that restraints derived from site-directed spin labeling can contribute significantly to defining the orientations and conformations of bound ligands. Accurate ligand localization appears to require either a few supplementary short-range distance restraints, or relatively tight spin-label restraints, with at least one spin label positioned so that some of the restraints draw the ligand into the binding pocket in the latter case.     

Conformational Transitions upon Ligand Binding: Holo-Structure Prediction from Apo Conformations

Biological function of proteins is frequently associated with the formation of complexes with small-molecule ligands. Experimental structure determination of such complexes at atomic resolution, however, can be time-consuming and costly. Computational methods for structure prediction of protein/ligand complexes, particularly docking, are as yet restricted by their limited consideration of receptor flexibility, rendering them not applicable for predicting protein/ligand complexes if large conformational changes of the receptor upon ligand binding are involved. Accurate receptor models in the ligand-bound state (holo structures), however, are a prerequisite for successful structure-based drug design. Hence, if only an unbound (apo) structure is available distinct from the ligand-bound conformation, structure-based drug design is severely limited. We present a method to predict the structure of protein/ligand complexes based solely on the apo structure, the ligand and the radius of gyration of the holo structure. The method is applied to ten cases in which proteins undergo structural rearrangements of up to 7.1 Å backbone RMSD upon ligand binding. In all cases, receptor models within 1.6 Å backbone RMSD to the target were predicted and close-to-native ligand binding poses were obtained for 8 of 10 cases in the top-ranked complex models. A protocol is presented that is expected to enable structure modeling of protein/ligand complexes and structure-based drug design for cases where crystal structures of ligand-bound conformations are not available.     

Applying conformational selection theory to improve crossdocking efficiency in 3‐phosphoinositide dependent protein kinase‐1

The emerging picture of biomolecular recognition is that of conformational selection followed by induced‐fit. Conformational selection theory states that binding partners exist in various conformations in solution, with binding involving a “selection” between complementary conformers. In this study, we devise a docking protocol that mimics conformational selection in protein–ligand binding and demonstrate that it significantly enhances crossdocking accuracy over Glide's flexible docking protocol, which is widely used in the pharmaceutical industry. Our protocol uses a pregenerated conformational ensemble to simulate ligand flexibility. The ensemble was generated by thorough conformational sampling coupled with conformer minimization. The generated conformers were then rigidly docked in the active site of the protein along with a postdocking minimization step that allows limited induced fit effects to be modeled for the ligand. We illustrate the improved performance of our protocol through crossdocking of 31 ligands to cocomplexed proteins of the kinase 3‐phosphoinositide dependent protein kinase‐1 extracted from the crystal structures 1H1W (ATP bound), 1OKY (staurosporine bound) and 3QD0 (bound to a potent inhibitor). Consistent with conformational selection theory, the performance of our protocol was the best for crossdocking to the cognate protein bound to the natural ligand, ATP. Proteins 2014; 82:436–451. © 2013 Wiley Periodicals, Inc.     

Incorporating biochemical information and backbone flexibility in RosettaDock for CAPRI rounds 6-12

In CAPRI rounds 6-12, RosettaDock successfully predicted 2 of 5 unbound-unbound targets to medium accuracy. Improvement over the previous method was achieved with computational mutagenesis to select decoys that match the energetics of experimentally determined hot spots. In the case of Target 21, Orc1/Sir1, this resulted in a successful docking prediction where RosettaDock alone or with simple site constraints failed. Experimental information also helped limit the interacting region of TolB/Pal, producing a successful prediction of Target 26. In addition, we docked multiple loop conformations for Target 20, and we developed a novel flexible docking algorithm to simultaneously optimize backbone conformation and rigid-body orientation to generate a wide diversity of conformations for Target 24. Continued challenges included docking of homology targets that differ substantially from their template (sequence identity <50%) and accounting for large conformational changes upon binding. Despite a larger number of unbound-unbound and homology model binding targets, Rounds 6-12 reinforced that RosettaDock is a powerful algorithm for predicting bound complex structures, especially when combined with experimental data.     

Ligand-induced conformational capture of a synthetic tetracycline riboswitch revealed by pulse EPR

RNA aptamers are in vitro-selected binding domains that recognize their respective ligand with high affinity and specificity. They are characterized by complex three-dimensional conformations providing preformed binding pockets that undergo conformational changes upon ligand binding. Small molecule-binding aptamers have been exploited as synthetic riboswitches for conditional gene expression in various organisms. In the present study, double electron-electron resonance (DEER) spectroscopy combined with site-directed spin labeling was used to elucidate the conformational transition of a tetracycline aptamer upon ligand binding. Different sites were selected for post-synthetic introduction of either the (1-oxyl-2,2,5,5-tetramethylpyrroline-3-methyl) methanethiosulfonate by reaction with a 4-thiouridine modified RNA or of 4-isocyanato-2,6-tetramethylpiperidyl-N-oxid spin label by reaction with 2'-aminouridine modified RNA. The results of the DEER experiments indicate the presence of a thermodynamic equilibrium between two aptamer conformations in the free state and capture of one conformation upon tetracycline binding.     

Solvated protein–DNA docking using HADDOCK

Interfacial water molecules play an important role in many aspects of protein–DNA specificity and recognition. Yet they have been mostly neglected in the computational modeling of these complexes. We present here a solvated docking protocol that allows explicit inclusion of water molecules in the docking of protein–DNA complexes and demonstrate its feasibility on a benchmark of 30 high-resolution protein–DNA complexes containing crystallographically-determined water molecules at their interfaces. Our protocol is capable of reproducing the solvation pattern at the interface and recovers hydrogen-bonded water-mediated contacts in many of the benchmark cases. Solvated docking leads to an overall improvement in the quality of the generated protein–DNA models for cases with limited conformational change of the partners upon complex formation. The applicability of this approach is demonstrated on real cases by docking a representative set of 6 complexes using unbound protein coordinates, model-built DNA and knowledge-based restraints. As HADDOCK supports the inclusion of a variety of NMR restraints, solvated docking is also applicable for NMR-based structure calculations of protein–DNA complexes.     

Accuracy of binding mode prediction with a cascadic stochastic tunneling method

We investigate the accuracy of the binding modes predicted for 83 complexes of the high-resolution subset of the ASTEX/CCDC receptor-ligand database using the atomistic FlexScreen approach with a simple forcefield-based scoring function. The median RMS deviation between experimental and predicted binding mode was just 0.83 A. Over 80% of the ligands dock within 2 A of the experimental binding mode, for 60 complexes the docking protocol locates the correct binding mode in all of ten independent simulations. Most docking failures arise because (a) the experimental structure clashed in our forcefield and is thus unattainable in the docking process or (b) because the ligand is stabilized by crystal water.     

Enhanced sampling of protein conformational states for dynamic cross-docking within the protein-protein docking server SwarmDock

The formation of specific protein-protein interactions is often a key to a protein's function. During complex formation, each protein component will undergo a change in the conformational state, for some these changes are relatively small and reside primarily at the sidechain level; however, others may display notable backbone adjustments. One of the classic problems in the protein-docking field is to be able to a priori predict the extent of such conformational changes. In this work, we investigated three protocols to find the most suitable input structure conformations for cross-docking, including a robust sampling approach in normal mode space. Counterintuitively, knowledge of the theoretically best combination of normal modes for unbound-bound transitions does not always lead to the best results. We used a novel spatial partitioning library, Aether Engine (see Supplementary Materials ), to efficiently search the conformational states of 56 receptor/ligand pairs, including a recent CAPRI target, in a systematic manner and selected diverse conformations as input to our automated docking server, SwarmDock, a server that allows moderate conformational adjustments during the docking process. In essence, here we present a dynamic cross-docking protocol, which when benchmarked against the simpler approach of just docking the unbound components shows a 10% uplift in the quality of the top docking pose.     

A Fast Flexible Docking Method using an Incremental Construction Algorithm

We present an automatic method for docking organic ligands into protein binding sites. The method can be used in the design process of specific protein ligands. It combines an appropriate model of the physico-chemical properties of the docked molecules with efficient methods for sampling the conformational space of the ligand. If the ligand is flexible, it can adopt a large variety of different conformations. Each such minimum in conformational space presents a potential candidate for the conformation of the ligand in the complexed state. Our docking method samples the conformation space of the ligand on the basis of a discrete model and uses a tree-search technique for placing the ligand incrementally into the active site. For placing the first fragment of the ligand into the protein, we use hashing techniques adapted from computer vision. The incremental construction algorithm is based on a greedy strategy combined with efficient methods for overlap detection and for the search of new interactions. We present results on 19 complexes of which the binding geometry has been crystallographically determined. All considered ligands are docked in at most three minutes on a current workstation. The experimentally observed binding mode of the ligand is reproduced with 0.5 to 1.2 Å rms deviation. It is almost always found among the highest-ranking conformations computed.     

The relationship between the flexibility of proteins and their conformational states on forming protein-protein complexes with an application to protein-protein docking

We investigate the extent to which the conformational fluctuations of proteins in solution reflect the conformational changes that they undergo when they form binary protein-protein complexes. To do this, we study a set of 41 proteins that form such complexes and whose three-dimensional structures are known, both bound in the complex and unbound. We carry out molecular dynamics simulations of each protein, starting from the unbound structure, and analyze the resulting conformational fluctuations in trajectories of 5 ns in length, comparing with the structure in the complex. It is found that fluctuations take some parts of the molecules into regions of conformational space close to the bound state (or give information about it), but at no point in the simulation does each protein as whole sample the complete bound state. Subsequent use of conformations from a clustered MD ensemble in rigid-body docking is nevertheless partially successful when compared to docking the unbound conformations, as long as the unbound conformations are themselves included with the MD conformations and the whole globally rescored. For one key example where sub-domain motion is present, a ribonuclease inhibitor, principal components analysis of the MD was applied and was also able to produce conformations for docking that gave enhanced results compared to the unbound. The most significant finding is that core interface residues show a tendency to be less mobile (by size of fluctuation or entropy) than the rest of the surface even when the other binding partner is absent, and conversely the peripheral interface residues are more mobile. This surprising result, consistent across up to 40 of the 41 proteins, suggests different roles for these regions in protein recognition and binding, and suggests ways that docking algorithms could be improved by treating these regions differently in the docking process.     

Molecular modeling studies to predict the possible binding modes of endomorphin analogs in μ opioid receptor

The molecular docking of a series of endomorphin analog with the μ opioid receptor was performed. The successive molecular dynamics of several proposed ligand–receptor complexes inserted into the phospholipid bilayer were carried out to optimize the complex and explore the conformational changes. Meaningful differences of their binding modes were detected and the involvement of some essential residues in ligand binding was also identified. Our proposed ligand–receptor model is in good agreement with previous site-directed mutagenesis experiments.     

Combining NMR and small angle X-ray and neutron scattering in the structural analysis of a ternary protein-RNA complex

Many processes in the regulation of gene expression and signaling involve the formation of protein complexes involving multi-domain proteins. Individual domains that mediate protein-protein and protein-nucleic acid interactions are typically connected by flexible linkers, which contribute to conformational dynamics and enable the formation of complexes with distinct binding partners. Solution techniques are therefore required for structural analysis and to characterize potential conformational dynamics. Nuclear magnetic resonance spectroscopy (NMR) provides such information but often only sparse data are obtained with increasing molecular weight of the complexes. It is therefore beneficial to combine NMR data with additional structural restraints from complementary solution techniques. Small angle X-ray/neutron scattering (SAXS/SANS) data can be efficiently combined with NMR-derived information, either for validation or by providing additional restraints for structural analysis. Here, we show that the combination of SAXS and SANS data can help to refine structural models obtained from data-driven docking using HADDOCK based on sparse NMR data. The approach is demonstrated with the ternary protein-protein-RNA complex involving two RNA recognition motif (RRM) domains of Sex-lethal, the N-terminal cold shock domain of Upstream-to-N-Ras, and msl-2 mRNA. Based on chemical shift perturbations we have mapped protein-protein and protein-RNA interfaces and complemented this NMR-derived information with SAXS data, as well as SANS measurements on subunit-selectively deuterated samples of the ternary complex. Our results show that, while the use of SAXS data is beneficial, the additional combination with contrast variation in SANS data resolves remaining ambiguities and improves the docking based on chemical shift perturbations of the ternary protein-RNA complex.