Ligand binding characteristics of two molecular forms,one equipped with a hydrophobic glycosyl phosphatidylinositol tail,of the folate binding protein purified from human milk

Two molecular forms of the folate binding protein were isolated and purified from human milk by a combination of cation exchange- and affinity chromatography. One protein (27 kDa) was a cleavage product of the other 100 kDa protein as evidenced by N-terminal amino acid sequence homology and a reduction in the molecular size of the latter protein to 27 kDa after cleavage of its hydrophobic glycosylphosphatidylinositol tail by phosphatidylinositol-specific phospholipase C. High-affinity binding of [³H]folate was characterized by upward convex Scatchard plots and increasing ligand binding affinity with decreasing concentrations of both proteins. Downward convex Scatchard plots and binding affinities showing no dependence on the protein concentration were, however, observed in highly diluted solutions of both proteins. Radioligand binding was inhibited by folate analogs, and dissociation of radioligand was slow at pH 7.4 but rapid and complete at pH 5.0 and 3.5. Ligand binding quenched the tryptophan fluorescence of the 27 kDa protein suggesting that tryptophan is present at the binding site and/or ligand binding induces a conformation change that affects tryptophan environment in the protein. The 27 kDa protein representing soluble folate binding protein exhibited a greater affinity for ligand binding than the 100 kDa protein which possesses a hydrophobic tail identical to the one that anchors the folate receptor to the cell membrane.     

Disaccharide binding to galectin-1: free energy calculations and molecular recognition mechanism

Galectin-1, a member of the conserved family of carbohydrate-binding proteins with affinity for β-galactosides, is a key modulator of diverse cell functions such as immune response and regulation. The binding affinity and specificity of galectin-1 for eight different β-galactosyl terminal disaccharides was studied using molecular-dynamics simulations in which the ligand was pulled away from the binding site using a mechanical force. We present what we believe to be a novel procedure, based on combinations of multistep trajectories, that was used to estimate the binding free energy (ΔG) of each disaccharide. The computed binding free energy differences show excellent correlation with experimental values determined previously. The small differences in affinity among the disaccharides are the result of an exquisite balance between the strengths of the galectin-sugar H-bonds and the H-bonds the protein and the disaccharides make with the solvent. Analysis of the free energies along the reaction coordinate shows that disaccharide unbinding/binding presents no energetic barrier and, therefore, is diffusion-limited. In addition, the calculations revealed that as the ligand is undocked from the binding site, breaking of protein-disaccharide H-bonds takes place in stages with intermediate states in which the interactions are bridged by water molecules.     

A combination of cation exchange and ligand-affinity chromatography for purification of two molecular species of the folate binding protein in human milk,one equipped with a hydrophobic glycosyl phosphatidylinositol tail: characterization of hydrophobicity and electrical charge

Cation exchange chromatography combined with ligand (methotrexate) affinity chromatography on a column desorbed with a pH-gradient was used for separation and large scale purification of two folate binding proteins in human milk. One of the proteins, which had a molecular size of 27 kDa on gel filtration and eluted from the affinity column at pH 5-6 was a cleavage product of a 100 kDa protein eluted at pH 3-4 as evidenced by identical N-terminal amino acid sequences and a reduction in the molecular size of the latter protein to 27 kDa after cleavage of its hydrophobic glycosylphosphatidyl-inositol tail that inserts into Triton X-100 micelles. Chromatofocusing showed that both proteins possessed multiple isoelectric points within the pH range 7-9. The 100 kDa protein exhibited a high affinity to hydrophobic interaction chromatographic gels, whereas this was only the case with unliganded forms of the 27 kDa protein indicative of a decrease in the hydrophobicity of the protein after ligand binding.     

Conformational Destabilization of Immunoglobulin G Increases the Low pH Binding Affinity with the Neonatal Fc Receptor

Crystallographic evidence suggests that the pH-dependent affinity of IgG molecules for the neonatal Fc receptor (FcRn) receptor primarily arises from salt bridges involving IgG histidine residues, resulting in moderate affinity at mildly acidic conditions. However, this view does not explain the diversity in affinity found in IgG variants, such as the YTE mutant (M252Y,S254T,T256E), which increases affinity to FcRn by up to 10×. Here we compare hydrogen exchange measurements at pH 7.0 and pH 5.5 with and without FcRn bound with surface plasmon resonance estimates of dissociation constants and FcRn affinity chromatography. The combination of experimental results demonstrates that differences between an IgG and its cognate YTE mutant vary with their pH-sensitive dynamics prior to binding FcRn. The conformational dynamics of these two molecules are nearly indistinguishable upon binding FcRn. We present evidence that pH-induced destabilization in the CH2/3 domain interface of IgG increases binding affinity by breaking intramolecular H-bonds and increases side-chain adaptability in sites that form intermolecular contacts with FcRn. Our results provide new insights into the mechanism of pH-dependent affinity in IgG-FcRn interactions and exemplify the important and often ignored role of intrinsic conformational dynamics in a protein ligand, to dictate affinity for biologically important receptors.     

Properties of purified folate-binding proteins from chronic myelogenous leukemia cells

Folate-binding protein(s) from chronic myelogenous leukemia cells have been purified using acid dialysis, ammonium sulfate fractionation and affinity chromatography. The purified preparation which migrates as a single band on disc electrophoresis could be separated by DEAE agarose chromatography into two folate-binding proteins (binders I and II) which bind molar equivalents of folic acid. One binder (I) eluted from DEAE at 1 mM sodium phosphate, pH 6.0, and the other (II) at 100 mM sodium phosphate, pH 7.4. Analysis of the purified mixture, which contained more than 90% binder II, by sedimentation equilibrium centrifugation indicated a homogeneous protein with a calculated molecular weight of 44000. Antiserum raised against the purified mixture gave a single precipitin line by immunodiffusion against a preparation of partially purified cell lysate.Hydrolysis of the more acidic binder (II) with neuraminidase converted it to a weakly acidic protein similar to binder I suggesting that these binders are glycoproteins which differ in sialic acid content. With isoelectric focusing, the binding of folic acid would be demonstrated at pH 6.7, 7.3, 7.8 and 8.2 for binder I, and at pH 5.1, 5.8 and 6.5 for binder II. Binders I and II had equally high affinity for folic acid and dihydroflate, lower affinity of N⁵-methyl-tetrahydrofolate, and no apparent affinity for N⁵-formytetrahydrofolate or methotrexate.     

Reverse and conventional chemical ecology approaches for the development of oviposition attractants for Culex mosquitoes

Synthetic mosquito oviposition attractants are sorely needed for surveillance and control programs for Culex species, which are major vectors of pathogens causing various human diseases, including filariasis, encephalitis, and West Nile encephalomyelitis. We employed novel and conventional chemical ecology approaches to identify potential attractants, which were demonstrated in field tests to be effective for monitoring populations of Cx. p. quinquefasciatus in human dwellings. Immunohistochemistry studies showed that an odorant-binding protein from this species, CquiOBP1, is expressed in trichoid sensilla on the antennae, including short, sharp-tipped trichoid sensilla type, which house an olfactory receptor neuron sensitive to a previously identified mosquito oviposition pheromone (MOP), 6-acetoxy-5-hexadecanolide. CquiOBP1 exists in monomeric and dimeric forms. Monomeric CquiOBP1 bound MOP in a pH-dependent manner, with a change in secondary structure apparently related to the loss of binding at low pH. The pheromone antipode showed higher affinity than the natural stereoisomer. By using both CquiOBP1 as a molecular target in binding assays and gas chromatography-electroantennographic detection (GC-EAD), we identified nonanal, trimethylamine (TMA), and skatole as test compounds. Extensive field evaluations in Recife, Brazil, a region with high populations of Cx. p. quinquefasciatus, showed that a combination of TMA (0.9 microg/l) and nonanal (0.15 ng/microl) is equivalent in attraction to the currently used infusion-based lure, and superior in that the offensive smell of infusions was eliminated in the newly developed synthetic mixture.     

Insights on pH-dependent conformational changes of mosquito odorant binding proteins by molecular dynamics simulations

Chemical recognition plays an important role for the survival and reproduction of many insect species. Odorant binding proteins (OBPs) are the primary components of the insect olfactory mechanism and have been documented to play an important role in the host-seeking mechanism of mosquitoes. They are “transport proteins” believed to transport odorant molecules from the external environment to their respective membrane targets, the olfactory receptors. The mechanism by which this transport occurs in mosquitoes remains a conundrum in this field. Nevertheless, OBPs have proved to be amenable to conformational changes mediated by a pH change in other insect species. In this paper, the effect of pH on the conformational flexibility of mosquito OBPs is assessed computationally using molecular dynamics simulations of a mosquito OBP “CquiOBP1” bound to its pheromone 3OG (PDB ID: 3OGN). Conformational twist of a loop, driven by a set of well-characterized changes in intramolecular interactions of the loop, is demonstrated. The concomitant (i) closure of what is believed to be the entrance of the binding pocket, (ii) expansion of what could be an exit site, and (iii) migration of the ligand towards this putative exit site provide preliminary insights into the mechanism of ligand binding and release of these proteins in mosquitoes. The correlation of our results with previous experimental observations based on NMR studies help us provide a cardinal illustration on one of the probable dynamics and mechanism by which certain mosquito OBPs could deliver their ligand to their membrane-bound receptors at specific pH conditions.     

The crystal structure of an odorant binding protein from Anopheles gambiae: evidence for a common ligand release mechanism

The Anopheles gambiae mosquito is the main vector of malaria transmission in sub-Saharan Africa. We present here a 1.5A crystal structure of AgamOBP1, an odorant binding protein (OBP) from the A. gambiae mosquito. The protein crystallized as a dimer with a unique binding pocket consisting of a continuous tunnel running through both subunits of the dimer and occupied by a PEG molecule. We demonstrate that AgamOBP1 undergoes a pH dependent conformational change that is associated with reduced ligand binding. A predominance of acid-labile hydrogen bonds involving the C-terminal loop suggests a mechanism in which a drop in pH causes C-terminal loop to open, leaving the binding tunnel solvent exposed, thereby lowering binding affinity for ligand. Because proteins from two distantly related insects also undergo a pH dependent conformational change involving the C-terminus that is associated with reduced ligand affinity, our results suggest a common mechanism for OBP activity.     

Investigation of silent information regulator 1 (Sirt1) agonists from Traditional Chinese Medicine

Silent information regulator 1 (Sirt1), a class III nicotinamide adenine dinucleotide dependent histone deacetylases, is important in cardioprotection, neuroprotection, metabolic disease, calorie restriction, and diseases associated with aging. Traditional Chinese Medicine (TCM) compounds from TCM Database@Taiwan (http://tcm.cmu.edu.tw/) were employed for screening potent Sirt1 agonists, and molecular dynamics (MD) simulation was implemented to simulate ligand optimum docking poses and protein structure under dynamic conditions. TCM compounds such as (S)-tryptophan-betaxanthin, 5-O-feruloylquinic acid, and RosA exhibited good binding affinity across different computational methods, and their drug-like potential were validated by MD simulation. Docking poses indicate that the carboxylic group of the three candidates generated H-bonds with residues in the protein chain from Ser441 to Lys444 and formed H-bond, π–cation interactions, or hydrophobic contacts with Phe297 and key active residue, His363. During MD, stable π–cation interactions with residues Phe273 or Arg274 were formed by (S)-tryptophan-betaxanthin and RosA. All candidates were anchored to His363 by stable π- or H-bonds. Hence, we propose (S)-tryptophan-betaxanthin, 5-O-feruloylquinic acid, and RosA as potential lead compounds that can be further tested in drug development process for diseases associated with aging

An animated interactive 3D complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:28     

Dynamic conformational equilibria in the physiological function of the Bombyx mori pheromone-binding protein

The Bombyx mori pheromone-binding protein (BmorPBP) undergoes a pH-dependent conformational transition from a form at basic pH, which contains an open cavity suitable for ligand binding (BmorPBP^B), to a form at pH 4.5, where this cavity is occupied by an additional helix (BmorPBP^A). This helix α7 is formed by the C-terminal dodecapeptide 131-142, which is flexibly disordered on the protein surface in BmorPBP^B and in its complex with the pheromone bombykol. Previous work showed that the ligand-binding cavity cannot accommodate both bombykol and helix α7. Here we further investigated mechanistic aspects of the physiologically crucial ejection of the ligand at lower pH values by solution NMR studies of the variant protein BmorPBP(1-128), where the C-terminal helix-forming tetradecapeptide is removed. The NMR structure of the truncated protein at pH 6.5 corresponds closely to BmorPBP^B. At pH 4.5, BmorPBP(1-128) maintains a B-type structure that is in a slow equilibrium, on the NMR chemical shift timescale, with a low-pH conformation for which a discrete set of ¹⁵N-¹H correlation peaks is NMR unobservable. The full NMR spectrum was recovered upon readjusting the pH of the protein solution to 6.5. These data reveal dual roles for the C-terminal tetradecapeptide of BmorPBP in the mechanism of reversible pheromone binding and transport, where it governs dynamic equilibria between two locally different protein conformations at acidic pH and competes with the ligand for binding to the interior cavity.     

Sulfate-binding protein dislikes protonated oxyacids. A molecular explanation

We have determined the effect of pH on the binding affinities of the conjugate bases of four different tetrahedral oxyacids to the sulfate-binding protein. The equilibrium dissociation constants of the binding of sulfate (Kd = 0.12 microM) and selenate (Kd = 5 microM) were found to be pH independent over the range pH 5 to pH 8.1, whereas chromate binding exhibited a pH dependence that is approximately attributable to the pK2 of the chromic acid. Phosphate was bound with an affinity five orders of magnitude weaker than that of sulfate. In light of the highly refined 2 A structure of the complex of the sulfate-binding protein with sulfate, and considering the protonation state and net charge of the various oxyacids, we conclude that the pH dependence of chromate binding and the extremely low affinity of phosphate are attributable mainly to a lack of hydrogen bond acceptors in the binding site. These studies demonstrate that the sulfate-binding site is stringently designed to bind tightly tetrahedral, fully ionized, oxyacid dianions. The presence of a donatable proton on the ligand reduces binding energy by approximately 7 kcal/mol.     

Multiple functions of aromatic-carbohydrate interactions in a processive cellulase examined with molecular simulation

Proteins employ aromatic residues for carbohydrate binding in a wide range of biological functions. Glycoside hydrolases, which are ubiquitous in nature, typically exhibit tunnels, clefts, or pockets lined with aromatic residues for processing carbohydrates. Mutation of these aromatic residues often results in significant activity differences on insoluble and soluble substrates. However, the thermodynamic basis and molecular level role of these aromatic residues remain unknown. Here, we calculate the relative ligand binding free energy by mutating tryptophans in the Trichoderma reesei family 6 cellulase (Cel6A) to alanine. Removal of aromatic residues near the catalytic site has little impact on the ligand binding free energy, suggesting that aromatic residues immediately upstream of the active site are not directly involved in binding, but play a role in the glucopyranose ring distortion necessary for catalysis. Removal of aromatic residues at the entrance and exit of the Cel6A tunnel, however, dramatically impacts the binding affinity, suggesting that these residues play a role in chain acquisition and product stabilization, respectively. The roles suggested from differences in binding affinity are confirmed by molecular dynamics and normal mode analysis. Surprisingly, our results illustrate that aromatic-carbohydrate interactions vary dramatically depending on the position in the enzyme tunnel. As aromatic-carbohydrate interactions are present in all carbohydrate-active enzymes, these results have implications for understanding protein structure-function relationships in carbohydrate metabolism and recognition, carbon turnover in nature, and protein engineering strategies for biomass utilization. Generally, these results suggest that nature employs aromatic-carbohydrate interactions with a wide range of binding affinities for diverse functions.     

Histidine-mediated pH-sensitive regulation of M-ficolin:GlcNAc binding activity in innate immunity examined by molecular dynamics simulations

Background

M-ficolin, a pathogen recognition molecule in the innate immune system, binds sugar residues including N-acetyl-D-glucosamine (GlcNAc), which is displayed on invading microbes and on apoptotic cells. The cis and trans Asp282-Cys283 peptide bond in the M-ficolin, which was found to occur at neutral and acidic pH in crystal structures, has been suggested to represent binding and non-binding activity, respectively. A detailed understanding of the pH-dependent conformational changes in M-ficolin and pH-mediated discrimination mechanism of GlcNAc-binding activity are crucial to both immune-surveillance and clearance of apoptotic cells.

Methodology/Principal Findings

By immunodetection analysis, we found that the pH-sensitive binding of GlcNAc is regulated by a conformational equilibrium between the active and inactive states of M-ficolin. We performed constant pH molecular dynamics (MD) simulation at a series of pH values to explore the pH effect on the cis-trans isomerization of the Asp282-Cys283 peptide bond in the M-ficolin fibrinogen-like domain (FBG). Analysis of the hydrogen bond occupancy of wild type FBG compared with three His mutants (H251A, H284A and H297A) corroborates that His284 is indispensible for pH-dependent binding. H251A formed new but weaker hydrogen bonds with GlcNAc. His297, unlike the other two His mutants, is more dependent on the solution pH and also contributes to cis-trans isomerization of the Asp282-Cys283 peptide bond in weak basic solution.

Conclusions/Significance

Constant pH MD simulation indicated that the cis active isomer of Asp282-Cys283 peptide bond was predominant around neutral pH while the trans bond gradually prevailed towards acidic environment. The protonation of His284 was found to be associated with the trans-to-cis isomerization of Asp282-Cys283 peptide bond which dominantly regulates the GlcNAc binding. Our MD simulation approach provides an insight into the pH-sensitive proteins and hence, ligand binding activity.     

Theoretical investigation on the diatomic ligand migration process and ligand binding properties in non‐O2‐binding H‐NOX domain

The Nostoc sp (Ns) H‐NOX (heme‐nitric oxide or OXygen‐binding) domain shares 35% sequence identity with soluble guanylate cyclase (sGC) and exhibits similar ligand binding property with the sGC. Previously, our molecular dynamic (MD) simulation work identified that there exists a Y‐shaped tunnel system hosted in the Ns H‐NOX interior, which servers for ligand migration. The tunnels were then confirmed by Winter et al. [PNAS 2011;108(43):E 881–889] recently using x‐ray crystallography with xenon pressured conditions. In this work, to further investigate how the protein matrix of Ns H‐NOX modulates the ligand migration process and how the distal residue composition affects the ligand binding prosperities, the free energy profiles for nitric oxide (NO), carbon monooxide (CO), and O₂ migration are explored using the steered MDs simulation and the ligand binding energies are calculated using QM/MM schemes. The potential of mean force profiles suggest that the longer branch of the tunnel would be the most favorable route for NO migration and a second NO trapping site other than the distal heme pocket along this route in the Ns H‐NOX was identified. On the contrary, CO and O₂ would prefer to diffuse via the shorter branch of the tunnel. The QM/MM (quantum mechanics/molecular mechanics) calculations suggest that the hydrophobic distal pocket of Ns H‐NOX would provide an approximately vacuum environment and the ligand discrimination would be determined by the intrinsic binding properties of the diatomic gas ligand to the heme group. Proteins 2013; 81:1363–1376. © 2013 Wiley Periodicals, Inc.     

The preparation and properties of folate-binding protein from cow''s milk.

An improved affinity-chromatographic method for the preparation of folate-binding protein from cow's milk is described. Under dissociating conditions the protein appeared homogeneous in the ultracentrifuge, with a molecular weight of 35 000 +/- 1500, but it was heterogeneous on electrophoresis and ion-exchange chromatography and evidently consisted of several glycoproteins with similar molecular weights that all bound folic acid. Overall, the protein contained a high proportion of half-cystine (18 residues/molecule) and 10.3% of carbohydrate. At saturation it bound approx. 1 mol of folate/mol of protein at pH 7.2. Equilibrium-dialysis measurements of the binding of folic acid and 5-methyltetrahydrofolate to the purified protein gave non-linear Scatchard plots, the shapes of which depended on pH. The results were interpreted in terms of ligand binding to a polymerizing system in which the affinity of ligand for monomer was greater than its affinity for polymer. When the protein concentration was similar to that in cow's milk, dissociation constants (Kd) for folate and 5-methyltetrahydrofolate were 3 nM and 5 nM respectively at pH 7.2 and 37 degrees C, whereas Kd for the binding of folate to monomer was about 50 pM. The properties of the binding protein are discussed in relation to its possible role in folate absorption in the gut.     

Lupinus luteus pathogenesis-related protein as a reservoir for cytokinin

Plant pathogenesis-related (PR) proteins of class 10 (PR-10) are small and cytosolic. The main feature of their three-dimensional structure is a large cavity between a seven-stranded antiparallel β-sheet and a long C-terminal α-helix. Although PR-10 proteins are abundant in plants, their physiological role remains unknown. Recent data have indicated ligand binding as their possible biological function. The article describes the structure of a complex between a classic PR-10 protein (yellow lupine LlPR-10.2B) and the plant hormone, trans-zeatin. Previously, trans-zeatin binding has been reported in a structurally related cytokinin-specific binding protein, which has a distant sequence relation with classic PR-10 proteins. In the present 1.35 Å resolution crystallographic model, three perfectly ordered zeatin molecules are found in the binding cavity of the protein. The fact that three zeatin molecules are bound by the protein when only a fourfold molar excess of the ligand was used indicates an unusual type of affinity for this ligand and suggests that LlPR-10.2B, and perhaps other PR-10 proteins as well, acts as a reservoir of cytokinin molecules in the aqueous environment of the cell.     

The molecular mechanisms of interactions between bioactive peptides and angiotensin-converting enzyme

The ability of milk protein derived Ile-Pro-Ala (IPA), Phe-Pro (FP) and Gly-Lys-Pro (GKP) peptides to inhibit angiotensin I-converting enzyme (ACE), a protein with an important role in blood-pressure regulation, were verified in vitro and in vivo. This work elucidates the modes and molecular mechanisms of the interaction of IPA, FP and GKP with ACE, including mechanisms that bind the peptides to the cofactor Zn²⁺. It was observed that the best docking poses obtained for IPA, FP and GKP were at the ACE catalytic site with very similar modes of interaction, including the interaction with Zn²⁺. The interactions, including H-bonds, hydrophobic, hydrophilic, and electrostatic interactions, as well as the interaction with Zn²⁺, were responsible for the binding between the bioactive peptides and ACE.     

A structural and mutagenic blueprint for molecular recognition of strychnine and d-tubocurarine by different cys-loop receptors

Cys-loop receptors (CLR) are pentameric ligand-gated ion channels that mediate fast excitatory or inhibitory transmission in the nervous system. Strychnine and d-tubocurarine (d-TC) are neurotoxins that have been highly instrumental in decades of research on glycine receptors (GlyR) and nicotinic acetylcholine receptors (nAChR), respectively. In this study we addressed the question how the molecular recognition of strychnine and d-TC occurs with high affinity and yet low specificity towards diverse CLR family members. X-ray crystal structures of the complexes with AChBP, a well-described structural homolog of the extracellular domain of the nAChRs, revealed that strychnine and d-TC adopt multiple occupancies and different ligand orientations, stabilizing the homopentameric protein in an asymmetric state. This introduces a new level of structural diversity in CLRs. Unlike protein and peptide neurotoxins, strychnine and d-TC form a limited number of contacts in the binding pocket of AChBP, offering an explanation for their low selectivity. Based on the ligand interactions observed in strychnine- and d-TC-AChBP complexes we performed alanine-scanning mutagenesis in the binding pocket of the human α1 GlyR and α7 nAChR and showed the functional relevance of these residues in conferring high potency of strychnine and d-TC, respectively. Our results demonstrate that a limited number of ligand interactions in the binding pocket together with an energetic stabilization of the extracellular domain are key to the poor selective recognition of strychnine and d-TC by CLRs as diverse as the GlyR, nAChR, and 5-HT(3)R.     

Discrete molecular dynamics distinguishes nativelike binding poses from decoys in difficult targets

Virtual screening is one of the major tools used in computer-aided drug discovery. In structure-based virtual screening, the scoring function is critical to identifying the correct docking pose and accurately predicting the binding affinities of compounds. However, the performance of existing scoring functions has been shown to be uneven for different targets, and some important drug targets have proven especially challenging. In these targets, scoring functions cannot accurately identify the native or near-native binding pose of the ligand from among decoy poses, which affects both the accuracy of the binding affinity prediction and the ability of virtual screening to identify true binders in chemical libraries. Here, we present an approach to discriminating native poses from decoys in difficult targets for which several scoring functions failed to correctly identify the native pose. Our approach employs Discrete Molecular Dynamics simulations to incorporate protein-ligand dynamics and the entropic effects of binding. We analyze a collection of poses generated by docking and find that the residence time of the ligand in the native and nativelike binding poses is distinctly longer than that in decoy poses. This finding suggests that molecular simulations offer a unique approach to distinguishing the native (or nativelike) binding pose from decoy poses that cannot be distinguished using scoring functions that evaluate static structures. The success of our method emphasizes the importance of protein-ligand dynamics in the accurate determination of the binding pose, an aspect that is not addressed in typical docking and scoring protocols.