Computational structure‐based drug design: Predicting target flexibility

The role of molecular modeling in drug design has experienced a significant revamp in the last decade. The increase in computational resources and molecular models, along with software developments, is finally introducing a competitive advantage in early phases of drug discovery. Medium and small companies with strong focus on computational chemistry are being created, some of them having introduced important leads in drug design pipelines. An important source for this success is the extraordinary development of faster and more efficient techniques for describing flexibility in three‐dimensional structural molecular modeling. At different levels, from docking techniques to atomistic molecular dynamics, conformational sampling between receptor and drug results in improved predictions, such as screening enrichment, discovery of transient cavities, etc. In this review article we perform an extensive analysis of these modeling techniques, dividing them into high and low throughput, and emphasizing in their application to drug design studies. We finalize the review with a section describing our Monte Carlo method, PELE, recently highlighted as an outstanding advance in an international blind competition and industrial benchmarks.We acknowledge the BSC-CRG-IRB Joint Research Program in Computational Biology. This work was supported by a grant from the Spanish Government CTQ2016-79138-R.J.I. acknowledges support from SVP-2014-068797, awarded by the Spanish Government.Peer ReviewedPostprint (author's final draft

Sampling of conformational ensemble for virtual screening using molecular dynamics simulations and normal mode analysis

Aim: Molecular dynamics simulations and normal mode analysis are well-established approaches to generate receptor conformational ensembles (RCEs) for ligand docking and virtual screening. Here, we report new fast molecular dynamics-based and normal mode analysis-based protocols combined with conformational pocket classifications to efficiently generate RCEs. Materials \& methods: We assessed our protocols on two well-characterized protein targets showing local active site flexibility, dihydrofolate reductase and large collective movements, CDK2. The performance of the RCEs was validated by distinguishing known ligands of dihydrofolate reductase and CDK2 among a dataset of diverse chemical decoys. Results \& discussion: Our results show that different simulation protocols can be efficient for generation of RCEs depending on different kind of protein flexibility

In silico assessment of potential druggable pockets on the surface of α1-Antitrypsin conformers

The search for druggable pockets on the surface of a protein is often performed on a single conformer, treated as a rigid body. Transient druggable pockets may be missed in this approach. Here, we describe a methodology for systematic in silico analysis of surface clefts across multiple conformers of the metastable protein α1-antitrypsin (A1AT). Pathological mutations disturb the conformational landscape of A1AT, triggering polymerisation that leads to emphysema and hepatic cirrhosis. Computational screens for small molecule inhibitors of polymerisation have generally focused on one major druggable site visible in all crystal structures of native A1AT. In an alternative approach, we scan all surface clefts observed in crystal structures of A1AT and in 100 computationally produced conformers, mimicking the native solution ensemble. We assess the persistence, variability and druggability of these pockets. Finally, we employ molecular docking using publicly available libraries of small molecules to explore scaffold preferences for each site. Our approach identifies a number of novel target sites for drug design. In particular one transient site shows favourable characteristics for druggability due to high enclosure and hydrophobicity. Hits against this and other druggable sites achieve docking scores corresponding to a Kd in the µM–nM range, comparing favourably with a recently identified promising lead. Preliminary ThermoFluor studies support the docking predictions. In conclusion, our strategy shows considerable promise compared with the conventional single pocket/single conformer approach to in silico screening. Our best-scoring ligands warrant further experimental investigation

Langerin-Heparin Interaction: Two Binding Sites for Small and Large Ligands as revealed by a combination of NMR Spectroscopy and Cross-Linking Mapping Experiments

Langerin is a C-type lectin present on Langerhans cells that mediates capture of pathogens in a carbohydrate-dependent manner, leading to subsequent internalization and elimination in the cellular organelles called Birbeck granules. This mechanism mediated by langerin was shown to constitute a natural barrier for HIV-1 particle transmission. Besides interacting specifically with high mannose and fucosylated neutral carbohydrate structures, langerin has the ability to bind sulfated carbohydrate ligands as 6-sulfated galactosides in the Ca2+ dependent binding site. Very recently langerin was demonstrated to interact with sulfated glycosaminoglycans (GAGs), in a Ca2+ independent way, resulting in the proposal of a new binding site for GAGs. Based on those results, we have conducted a structural study of the interactions of small heparin (HEP) like oligosaccharides with langerin in solution. Heparin-bead cross-linking experiments, an approach specifically designed to identify HEP/HS binding sites in proteins were first carried out and experimentally validated the previously proposed model for the interaction of Lg ECD with 6 kDa HEP. High-resolution NMR studies of a set of 8 synthetic HEP-like trisaccharides harboring different sulfation patterns demonstrated that all of them bound to langerin in a Ca2+ dependent way. The binding epitopes were determined by STD NMR and the bound conformations by transferred NOESY experiments. These experimental data were combined with docking and molecular dynamics and resulted in the proposal of a binding mode characterized by the coordination of calcium by the two equatorial hydroxyl groups OH3 and OH4 at the non-reducing end. The binding also includes the carboxylate group at the adjacent iduronate residue. Such epitope is shared by all the 8 ligands, explaining the absence of any impact on binding from their differences in substitution pattern. Finally, in contrast to the small trisaccharides, we demonstrated that a longer HEP-like hexasaccharide, bearing an additional O-sulfate group at the non-reducing end, which precludes binding to the Ca2+ site, interacts with langerin in the previously identified Ca2+ independent binding site

The benefits of in silico modeling to identify possible small-molecule drugs and their off-target interactions

Accepted for publication in a future issue of Future Medicinal Chemistry.The research into the use of small molecules as drugs continues to be a key driver in the development of molecular databases, computer-aided drug design software and collaborative platforms. The evolution of computational approaches is driven by the essential criteria that a drug molecule has to fulfill, from the affinity to targets to minimal side effects while having adequate absorption, distribution, metabolism, and excretion (ADME) properties. A combination of ligand- and structure-based drug development approaches is already used to obtain consensus predictions of small molecule activities and their off-target interactions. Further integration of these methods into easy-to-use workflows informed by systems biology could realize the full potential of available data in the drug discovery and reduce the attrition of drug candidates.Peer reviewe

A New Multi-Objective Approach for Molecular Docking Based on RMSD and Binding Energy

Ligand-protein docking is an optimization problem based on predicting the position of a ligand with the lowest binding energy in the active site of the receptor. Molecular docking problems are traditionally tackled with single-objective, as well as with multi-objective approaches, to minimize the binding energy. In this paper, we propose a novel multi-objective formulation that considers: the Root Mean Square Deviation (RMSD) difference in the coordinates of ligands and the binding (intermolecular) energy, as two objectives to evaluate the quality of the ligand-protein interactions. To determine the kind of Pareto front approximations that can be obtained, we have selected a set of representative multi-objective algorithms such as NSGA-II, SMPSO, GDE3, and MOEA/D. Their performances have been assessed by applying two main quality indicators intended to measure convergence and diversity of the fronts. In addition, a comparison with LGA, a reference single-objective evolutionary algorithm for molecular docking (AutoDock) is carried out. In general, SMPSO shows the best overall results in terms of energy and RMSD (value lower than 2A for successful docking results). This new multi-objective approach shows an improvement over the ligand-protein docking predictions that could be promising in in silico docking studies to select new anticancer compounds for therapeutic targets that are multidrug resistant.Universidad de Málaga. Campus de Excelencia Internacional Andalucía Tech