Challenges of docking in large, flexible and promiscuous binding sites

After decades of work, the correct determination of the binding mode of a small molecule into a target protein is still a challenging problem, whose difficulty depends on: (i) the sizes of the binding site and the ligand; (ii) the flexibility of both interacting partners, and (iii) the differential solvation of bound and unbound partners. We have evaluated the performance of standard rigid(receptor)/flexible(ligand) docking approaches with respect to last-generation fully flexible docking methods to obtain reasonable poses in a very challenging case: soluble Epoxide Hydrolase (sEH), a flexible protein showing different binding sites. We found that full description of the flexibility of both protein and ligand and accurate description of solvation leads to significant improvement in the ability of docking to reproduce well known binding modes, and at the same time capture the intrinsic binding promiscuity of the protein.Peer ReviewedPostprint (author's final draft

Pushing the Limits of Computational Structure-Based Drug Design with a Cryo-EM Structure: The Ca2+ Channel α2δ-1 Subunit as a Test Case

Cryo-electron microscopy (cryo-EM) is emerging as a real alternative for structural elucidation. In spite of this, very few cryo-EM structures have been described so far as successful platforms for in silico drug design. Gabapentin and pregabalin are some of the most successful drugs in the treatment of epilepsy and neuropathic pain. Although both are in clinical use and are known to exert their effects by binding to the regulatory α2δ subunit of voltage gated calcium channels, their binding modes have never been characterized. We describe here the successful use of an exhaustive protein–ligand sampling algorithm on the α2δ-1 subunit of the recently published cryo-EM structure, with the goal of characterizing the ligand entry path and binding mode for gabapentin, pregabalin, and several other amino acidic α2δ-1 ligands. Our studies indicate that (i) all simulated drugs explore the same path for accessing the occluded binding site on the interior of the α2δ-1 subunit; (ii) they all roughly occupy the same pocket; (iii) the plasticity of the binding site allows the accommodation of a variety of amino acidic modulators, driven by the flexible “capping loop” delineated by residues Tyr426-Val435 and the floppy nature of Arg217; (iv) the predicted binding modes are in line with previously available mutagenesis data, confirming Arg217 as key for binding, with Asp428 and Asp467 highlighted as additional anchoring points for all amino acidic drugs. The study is one of the first proofs that latest-generation cryo-EM structures combined with exhaustive computational methods can be exploited in early drug discovery.The authors would like to thank Modesto Orozco for fruitful discussions and feedback on the manuscript. Suggestions by the reviewers are also gratefully acknowledged. Nostrum is supported by Fundacion Marcelino Botin (Mind the Gap) and CDTI (Neotec grant −EXP 00094141/SNEO-20161127). M.K. and R.S. would like to thank BSC and IRB for technical support. V.G. is supported by the CTQ2016-79138-R grant.Peer ReviewedPostprint (author's final draft

Challenges of docking in large, flexible and promiscuous binding sites

After decades of work, the correct determination of the binding mode of a small molecule into a target protein is still a challenging problem, whose difficulty depends on: i) the sizes of the binding site and the ligand; ii) the flexibility of both interacting partners, and iii) the differential solvation of bound and unbound partners. We have evaluated the performance of standard rigid(receptor)/flexible(ligand) docking approaches with respect to last-generation fully flexible docking methods to obtain reasonable poses in a very challenging case: soluble Epoxide Hydrolase (sEH), a flexible protein showing different binding sites. We found that full description of the flexibility of both protein and ligand and accurate description of solvation leads to significant improvement in the ability of docking to reproduce well known binding modes, and at the same time capture the intrinsic binding promiscuity of the protein

Unveiling Prolyl Oligopeptidase Ligand Migration by Comprehensive Computational Techniques

Prolyl oligopeptidase (POP) is a large 80 kDa protease, which cleaves oligopeptides at the C-terminal side of proline residues and constitutes an important pharmaceutical target. Despite the existence of several crystallographic structures, there is an open debate about migration (entrance and exit) pathways for ligands, and their coupling with protein dynamics. Recent studies have shown the capabilities of molecular dynamics and classical force fields in describing spontaneous binding events and nonbiased ligand migration pathways. Due to POP's size and to the buried nature of its active site, an exhaustive sampling by means of conventional long enough molecular dynamics trajectories is still a nearly impossible task. Such a level of sampling, however, is possible with the breakthrough protein energy landscape exploration technique. Here, we present an exhaustive sampling of POP with a known inhibitor, Z-pro-prolinal. In >3000 trajectories Z-pro-prolinal explores all the accessible surface area, showing multiple entrance events into the large internal cavity through the pore in the β-propeller domain. Moreover, we modeled a natural substrate binding and product release by predicting the entrance of an undecapeptide substrate, followed by manual active site cleavage and nonbiased exit of one of the products (a dipeptide). The product exit shows preference from a flexible 18-amino acid residues loop, pointing to an overall mechanism where entrance and exit occur in different sites.This study was supported by The European Research Council (2009-Adg25027-PELE) to V.G., MICIN-FEDER (BIO2013-40716R) to E.G., the Generalitat de Catalunya (XRB and 2014SGR-521) to E.G., and SEV-2011-00067 of Severo Ochoa Program, awarded by the Spanish Government to D.L.Peer ReviewedPostprint (published version

Active site-directed inhibitors of prolyl oligopeptidase abolishes its conformational dynamics

Deciphering conformational dynamics is crucial for understanding the biological functions of proteins and for designing compounds targeting them. In particular, providing an accurate description of microsecond-millisecond motions opens the opportunity for regulating protein-protein interactions (PPIs) by modulating the dynamics of one interacting partner. Here we analyzed the conformational dynamics of prolyl oligopeptidase (POP) and the effects of active-site-directed inhibitors on the dynamics. We used an integrated structural biology approach based on NMR spectroscopy and SAXS experiments complemented by MD simulations. We found that POP is in a slow equilibrium in solution between open and closed conformations, and that inhibitors effectively abolished this equilibrium by stabilizing the enzyme in the closed conformation

Challenges of docking in large, flexible and promiscuous binding sites

After decades of work, the correct determination of the binding mode of a small molecule into a target protein is still a challenging problem, whose difficulty depends on: (i) the sizes of the binding site and the ligand; (ii) the flexibility of both interacting partners, and (iii) the differential solvation of bound and unbound partners. We have evaluated the performance of standard rigid(receptor)/flexible(ligand) docking approaches with respect to last-generation fully flexible docking methods to obtain reasonable poses in a very challenging case: soluble Epoxide Hydrolase (sEH), a flexible protein showing different binding sites. We found that full description of the flexibility of both protein and ligand and accurate description of solvation leads to significant improvement in the ability of docking to reproduce well known binding modes, and at the same time capture the intrinsic binding promiscuity of the protein.Peer Reviewe

Challenges of docking in large, flexible and promiscuous binding sites

After decades of work, the correct determination of the binding mode of a small molecule into a target protein is still a challenging problem, whose difficulty depends on: i) the sizes of the binding site and the ligand; ii) the flexibility of both interacting partners, and iii) the differential solvation of bound and unbound partners. We have evaluated the performance of standard rigid(receptor)/flexible(ligand) docking approaches with respect to last-generation fully flexible docking methods to obtain reasonable poses in a very challenging case: soluble Epoxide Hydrolase (sEH), a flexible protein showing different binding sites. We found that full description of the flexibility of both protein and ligand and accurate description of solvation leads to significant improvement in the ability of docking to reproduce well known binding modes, and at the same time capture the intrinsic binding promiscuity of the protein

Challenges of docking in large, flexible and promiscuous binding sites

After decades of work, the correct determination of the binding mode of a small molecule into a target protein is still a challenging problem, whose difficulty depends on: i) the sizes of the binding site and the ligand; ii) the flexibility of both interacting partners, and iii) the differential solvation of bound and unbound partners. We have evaluated the performance of standard rigid(receptor)/flexible(ligand) docking approaches with respect to last-generation fully flexible docking methods to obtain reasonable poses in a very challenging case: soluble Epoxide Hydrolase (sEH), a flexible protein showing different binding sites. We found that full description of the flexibility of both protein and ligand and accurate description of solvation leads to significant improvement in the ability of docking to reproduce well known binding modes, and at the same time capture the intrinsic binding promiscuity of the protein

Unveiling Prolyl Oligopeptidase Ligand Migration by Comprehensive Computational Techniques

Prolyl oligopeptidase (POP) is a large 80 kDa protease, which cleaves oligopeptides at the C-terminal side of proline residues and constitutes an important pharmaceutical target. Despite the existence of several crystallographic structures, there is an open debate about migration (entrance and exit) pathways for ligands, and their coupling with protein dynamics. Recent studies have shown the capabilities of molecular dynamics and classical force fields in describing spontaneous binding events and nonbiased ligand migration pathways. Due to POP's size and to the buried nature of its active site, an exhaustive sampling by means of conventional long enough molecular dynamics trajectories is still a nearly impossible task. Such a level of sampling, however, is possible with the breakthrough protein energy landscape exploration technique. Here, we present an exhaustive sampling of POP with a known inhibitor, Z-pro-prolinal. In >3000 trajectories Z-pro-prolinal explores all the accessible surface area, showing multiple entrance events into the large internal cavity through the pore in the β-propeller domain. Moreover, we modeled a natural substrate binding and product release by predicting the entrance of an undecapeptide substrate, followed by manual active site cleavage and nonbiased exit of one of the products (a dipeptide). The product exit shows preference from a flexible 18-amino acid residues loop, pointing to an overall mechanism where entrance and exit occur in different sites.This study was supported by The European Research Council (2009-Adg25027-PELE) to V.G., MICIN-FEDER (BIO2013-40716R) to E.G., the Generalitat de Catalunya (XRB and 2014SGR-521) to E.G., and SEV-2011-00067 of Severo Ochoa Program, awarded by the Spanish Government to D.L.Peer Reviewe

Pushing the Limits of Computational Structure-Based Drug Design with a Cryo-EM Structure: The Ca2+ Channel α2δ-1 Subunit as a Test Case

Cryo-electron microscopy (cryo-EM) is emerging as a real alternative for structural elucidation. In spite of this, very few cryo-EM structures have been described so far as successful platforms for in silico drug design. Gabapentin and pregabalin are some of the most successful drugs in the treatment of epilepsy and neuropathic pain. Although both are in clinical use and are known to exert their effects by binding to the regulatory α2δ subunit of voltage gated calcium channels, their binding modes have never been characterized. We describe here the successful use of an exhaustive protein–ligand sampling algorithm on the α2δ-1 subunit of the recently published cryo-EM structure, with the goal of characterizing the ligand entry path and binding mode for gabapentin, pregabalin, and several other amino acidic α2δ-1 ligands. Our studies indicate that (i) all simulated drugs explore the same path for accessing the occluded binding site on the interior of the α2δ-1 subunit; (ii) they all roughly occupy the same pocket; (iii) the plasticity of the binding site allows the accommodation of a variety of amino acidic modulators, driven by the flexible “capping loop” delineated by residues Tyr426-Val435 and the floppy nature of Arg217; (iv) the predicted binding modes are in line with previously available mutagenesis data, confirming Arg217 as key for binding, with Asp428 and Asp467 highlighted as additional anchoring points for all amino acidic drugs. The study is one of the first proofs that latest-generation cryo-EM structures combined with exhaustive computational methods can be exploited in early drug discovery.The authors would like to thank Modesto Orozco for fruitful discussions and feedback on the manuscript. Suggestions by the reviewers are also gratefully acknowledged. Nostrum is supported by Fundacion Marcelino Botin (Mind the Gap) and CDTI (Neotec grant −EXP 00094141/SNEO-20161127). M.K. and R.S. would like to thank BSC and IRB for technical support. V.G. is supported by the CTQ2016-79138-R grant.Peer Reviewe