1
|
Clemente CM, Prieto JM, Martí M. Unlocking Precision Docking for Metalloproteins. J Chem Inf Model 2024; 64:1581-1592. [PMID: 38373276 DOI: 10.1021/acs.jcim.3c01853] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/21/2024]
Abstract
Metalloproteins play a fundamental role in molecular biology, contributing to various biological processes. However, the discovery of high-affinity ligands targeting metalloproteins has been delayed due, in part, to a lack of suitable tools and data. Molecular docking, a widely used technique for virtual screening of small-molecule ligand interactions with proteins, often faces challenges when applied to metalloproteins due to the particular nature of the ligand metal bond. To address these limitations associated with docking metalloproteins, we introduce a knowledge-driven docking approach known as "metalloprotein bias docking" (MBD), which extends the AutoDock Bias technique. We assembled a comprehensive data set of metalloprotein-ligand complexes from 15 different metalloprotein families, encompassing Ca, Co, Fe, Mg, Mn, and Zn metal ions. Subsequently, we conducted a performance analysis of our MBD method and compared it to the conventional docking (CD) program AutoDock4, applied to various metalloprotein targets within our data set. Our results demonstrate that MBD outperforms CD, significantly enhancing accuracy, selectivity, and precision in ligand pose prediction. Additionally, we observed a positive correlation between our predicted ligand free energies and the corresponding experimental values. These findings underscore the potential of MBD as a valuable tool for the effective exploration of metalloprotein-ligand interactions.
Collapse
Affiliation(s)
- Camila M Clemente
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires (FCEyN-UBA) e Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN) CONICET, Pabellón 2 de Ciudad Universitaria, Ciudad de Buenos Aires C1428EHA, Argentina
| | - Juan M Prieto
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires (FCEyN-UBA) e Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN) CONICET, Pabellón 2 de Ciudad Universitaria, Ciudad de Buenos Aires C1428EHA, Argentina
| | - Marcelo Martí
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires (FCEyN-UBA) e Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN) CONICET, Pabellón 2 de Ciudad Universitaria, Ciudad de Buenos Aires C1428EHA, Argentina
| |
Collapse
|
2
|
Kudo G, Yanagisawa K, Yoshino R, Hirokawa T. AAp-MSMD: Amino Acid Preference Mapping on Protein-Protein Interaction Surfaces Using Mixed-Solvent Molecular Dynamics. J Chem Inf Model 2023; 63:7768-7777. [PMID: 38085669 PMCID: PMC10751795 DOI: 10.1021/acs.jcim.3c01677] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2023] [Revised: 11/28/2023] [Accepted: 11/29/2023] [Indexed: 12/26/2023]
Abstract
Peptides have attracted much attention recently owing to their well-balanced properties as drugs against protein-protein interaction (PPI) surfaces. Molecular simulation-based predictions of binding sites and amino acid residues with high affinity to PPI surfaces are expected to accelerate the design of peptide drugs. Mixed-solvent molecular dynamics (MSMD), which adds probe molecules or fragments of functional groups as solutes to the hydration model, detects the binding hotspots and cryptic sites induced by small molecules. The detection results vary depending on the type of probe molecule; thus, they provide important information for drug design. For rational peptide drug design using MSMD, we proposed MSMD with amino acid residue probes, named amino acid probe-based MSMD (AAp-MSMD), to detect hotspots and identify favorable amino acid types on protein surfaces to which peptide drugs bind. We assessed our method in terms of hotspot detection at the amino acid probe level and binding free energy prediction with amino acid probes at the PPI site for the complex structure that formed the PPI. In hotspot detection, the max-spatial probability distribution map (max-PMAP) obtained from AAp-MSMD detected the PPI site, to which each type of amino acid can bind favorably. In the binding free energy prediction using amino acid probes, ΔGFE obtained from AAp-MSMD roughly estimated the experimental binding affinities from the structure-activity relationship. AAp-MSMD, with amino acid probes, provides estimated binding sites and favorable amino acid types at the PPI site of a target protein.
Collapse
Affiliation(s)
- Genki Kudo
- Physics
Department, Graduate School of Pure and Applied Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba 305-8571, Ibaraki Japan
| | - Keisuke Yanagisawa
- Department
of Computer Science, School of Computing, Tokyo Institute of Technology, 2-12-1 Ookayama, Meguro 152-8550, Tokyo Japan
- Middle
Molecule IT-based Drug Discovery Laboratory, Tokyo Institute of Technology, 2-12-1 Ookayama, Meguro 152-8550, Tokyo Japan
| | - Ryunosuke Yoshino
- Faculty
of Medicine, University of Tsukuba, 1-1-1 Tennodai, Tsukuba 305-8575, Ibaraki Japan
- Transborder
Medical Research Center, University of Tsukuba, 1-1-1 Tennodai, Tsukuba 305-8577, Ibaraki Japan
| | - Takatsugu Hirokawa
- Faculty
of Medicine, University of Tsukuba, 1-1-1 Tennodai, Tsukuba 305-8575, Ibaraki Japan
- Transborder
Medical Research Center, University of Tsukuba, 1-1-1 Tennodai, Tsukuba 305-8577, Ibaraki Japan
| |
Collapse
|
3
|
Moreno-Cabezuelo JÁ, Del Carmen Muñoz-Marín M, López-Lozano A, Athayde D, Simón-García A, Díez J, Archer M, Issoglio FM, García-Fernández JM. Production, homology modeling and mutagenesis studies on GlcH glucose transporter from Prochlorococcus sp. strain SS120. BIOCHIMICA ET BIOPHYSICA ACTA. BIOENERGETICS 2023; 1864:148954. [PMID: 36563737 DOI: 10.1016/j.bbabio.2022.148954] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/01/2022] [Revised: 12/13/2022] [Accepted: 12/14/2022] [Indexed: 12/24/2022]
Abstract
The marine cyanobacterium Prochlorococcus is one of the main primary producers on Earth, which can take up glucose by using the high affinity, multiphasic transporter GlcH. We report here the overexpression of glcH from Prochlorococcus marinus strain SS120 in Escherichia coli. Modeling studies of GlcH using the homologous MelB melibiose transporter from Salmonella enterica serovar Typhimurium showed high conservation at the overall fold. We observed that an important structural interaction, mediated by a strong hydrogen bond between D8 and R141, is conserved in Prochlorococcus, although the corresponding amino acids in MelB from Salmonella are different. Biased docking studies suggested that when glucose reaches the pocket of the transporter and interacts with D8 and R141, the hydrogen bond network in which these residues are involved could be disrupted, favoring a conformational change with the subsequent translocation of the glucose molecule towards the cytoplasmic region of the pmGlcH structure. Based on these theoretical predictions and on the conservation of N117 and W348 in other MelB structures, D8, N117, R141 and W348 were mutated to glycine residues. Their key role in glucose transport was evaluated by glucose uptake assays. N117G and W348G mutations led to 17 % decrease in glucose uptake, while D8G and R141G decreased the glucose transport by 66 % and 92 % respectively. Overall, our studies provide insights into the Prochlorococcus 3D-structure of GlcH, paving the way for further analysis to understand the features which are involved in the high affinity and multiphasic kinetics of this transporter.
Collapse
Affiliation(s)
- José Ángel Moreno-Cabezuelo
- Departamento de Bioquímica y Biología Molecular, Campus de Excelencia Internacional Agroalimentario CeiA3, Universidad de Córdoba, Córdoba, Spain
| | - María Del Carmen Muñoz-Marín
- Departamento de Bioquímica y Biología Molecular, Campus de Excelencia Internacional Agroalimentario CeiA3, Universidad de Córdoba, Córdoba, Spain
| | - Antonio López-Lozano
- Departamento de Bioquímica y Biología Molecular, Campus de Excelencia Internacional Agroalimentario CeiA3, Universidad de Córdoba, Córdoba, Spain
| | - Diogo Athayde
- Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa (ITQB NOVA), 2780-157 Oeiras, Portugal
| | - Ana Simón-García
- Departamento de Bioquímica y Biología Molecular, Campus de Excelencia Internacional Agroalimentario CeiA3, Universidad de Córdoba, Córdoba, Spain
| | - Jesús Díez
- Departamento de Bioquímica y Biología Molecular, Campus de Excelencia Internacional Agroalimentario CeiA3, Universidad de Córdoba, Córdoba, Spain
| | - Margarida Archer
- Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa (ITQB NOVA), 2780-157 Oeiras, Portugal
| | - Federico M Issoglio
- Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa (ITQB NOVA), 2780-157 Oeiras, Portugal; CONICET-Universidad de Buenos Aires, Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN), Buenos Aires, Argentina.
| | - José Manuel García-Fernández
- Departamento de Bioquímica y Biología Molecular, Campus de Excelencia Internacional Agroalimentario CeiA3, Universidad de Córdoba, Córdoba, Spain.
| |
Collapse
|
4
|
Mayol GF, Defelipe LA, Arcon JP, Turjanski AG, Marti MA. Solvent Sites Improve Docking Performance of Protein–Protein Complexes and Protein–Protein Interface-Targeted Drugs. J Chem Inf Model 2022; 62:3577-3588. [DOI: 10.1021/acs.jcim.2c00264] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Gonzalo F. Mayol
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires (FCEyN-UBA) e Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN) CONICET, Pabellòn 2 de Ciudad Universitaria, Ciudad de Buenos Aires C1428EHA, Argentina
| | - Lucas A. Defelipe
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires (FCEyN-UBA) e Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN) CONICET, Pabellòn 2 de Ciudad Universitaria, Ciudad de Buenos Aires C1428EHA, Argentina
- European Molecular Biology Laboratory - Hamburg Unit, Notkestrasse 85, Hamburg 22607, Germany
| | - Juan Pablo Arcon
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires (FCEyN-UBA) e Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN) CONICET, Pabellòn 2 de Ciudad Universitaria, Ciudad de Buenos Aires C1428EHA, Argentina
- Institute for Research in Biomedicine (IRB), 08028 Barcelona, Spain
- The Barcelona Institute of Science and Technology, 08036 Barcelona, Spain
| | - Adrian G. Turjanski
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires (FCEyN-UBA) e Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN) CONICET, Pabellòn 2 de Ciudad Universitaria, Ciudad de Buenos Aires C1428EHA, Argentina
| | - Marcelo A. Marti
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires (FCEyN-UBA) e Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN) CONICET, Pabellòn 2 de Ciudad Universitaria, Ciudad de Buenos Aires C1428EHA, Argentina
| |
Collapse
|
5
|
Burastero O, Defelipe LA, Gola G, Tateosian NL, Lopez ED, Martinena CB, Arcon JP, Traian MD, Wetzler DE, Bento I, Barril X, Ramirez J, Marti MA, Garcia-Alai MM, Turjanski AG. Cosolvent Sites-Based Discovery of Mycobacterium Tuberculosis Protein Kinase G Inhibitors. J Med Chem 2022; 65:9691-9705. [PMID: 35737472 PMCID: PMC9344462 DOI: 10.1021/acs.jmedchem.1c02012] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
![]()
Computer-aided
drug discovery methods play a major role in the
development of therapeutically important small molecules, but their
performance needs to be improved. Molecular dynamics simulations in
mixed solvents are useful in understanding protein–ligand recognition
and improving molecular docking predictions. In this work, we used
ethanol as a cosolvent to find relevant interactions for ligands toward
protein kinase G, an essential protein of Mycobacterium
tuberculosis (Mtb).
We validated the hot spots by screening a database of fragment-like
compounds and another one of known kinase inhibitors. Next, we performed
a pharmacophore-guided docking simulation and found three low micromolar
inhibitors, including one with a novel chemical scaffold that we expanded
to four derivative compounds. Binding affinities were characterized
by intrinsic fluorescence quenching assays, isothermal titration calorimetry,
and the analysis of melting curves. The predicted binding mode was
confirmed by X-ray crystallography. Finally, the compounds significantly
inhibited the viability of Mtb in infected
THP-1 macrophages.
Collapse
Affiliation(s)
- Osvaldo Burastero
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pabellón 2, Buenos Aires C1428EGA, Argentina.,Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN), Universidad de Buenos Aires, Ciudad Universitaria, Pabellón 2, Buenos Aires C1428EGA, Argentina.,European Molecular Biology Laboratory Hamburg, Notkestrasse 85, Hamburg D-22607, Germany
| | - Lucas A Defelipe
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pabellón 2, Buenos Aires C1428EGA, Argentina.,Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN), Universidad de Buenos Aires, Ciudad Universitaria, Pabellón 2, Buenos Aires C1428EGA, Argentina.,European Molecular Biology Laboratory Hamburg, Notkestrasse 85, Hamburg D-22607, Germany
| | - Gabriel Gola
- Departamento de Química Orgánica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pabellón 2, Buenos Aires C1428EGA, Argentina.,Unidad de Microanálisis y Métodos Físicos Aplicados a Química Orgánica (UMYMFOR), Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires. CONICET, Buenos Aires C1428EGA, Argentina
| | - Nancy L Tateosian
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pabellón 2, Buenos Aires C1428EGA, Argentina.,Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN), Universidad de Buenos Aires, Ciudad Universitaria, Pabellón 2, Buenos Aires C1428EGA, Argentina
| | - Elias D Lopez
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pabellón 2, Buenos Aires C1428EGA, Argentina.,Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN), Universidad de Buenos Aires, Ciudad Universitaria, Pabellón 2, Buenos Aires C1428EGA, Argentina
| | - Camila Belen Martinena
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pabellón 2, Buenos Aires C1428EGA, Argentina.,Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN), Universidad de Buenos Aires, Ciudad Universitaria, Pabellón 2, Buenos Aires C1428EGA, Argentina
| | - Juan Pablo Arcon
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pabellón 2, Buenos Aires C1428EGA, Argentina.,Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN), Universidad de Buenos Aires, Ciudad Universitaria, Pabellón 2, Buenos Aires C1428EGA, Argentina
| | - Martín Dodes Traian
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pabellón 2, Buenos Aires C1428EGA, Argentina.,Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN), Universidad de Buenos Aires, Ciudad Universitaria, Pabellón 2, Buenos Aires C1428EGA, Argentina
| | - Diana E Wetzler
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pabellón 2, Buenos Aires C1428EGA, Argentina.,Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN), Universidad de Buenos Aires, Ciudad Universitaria, Pabellón 2, Buenos Aires C1428EGA, Argentina
| | - Isabel Bento
- European Molecular Biology Laboratory Hamburg, Notkestrasse 85, Hamburg D-22607, Germany
| | - Xavier Barril
- Catalan Institution for Research and Advanced Studies (ICREA), Passeig Lluís Companys 23, Barcelona 08010, Spain.,Faculty of Pharmacy and Institute of Biomedicine (IBUB), University of Barcelona, Av.Joan XXIII 27-31, Barcelona 08028, Spain
| | - Javier Ramirez
- Departamento de Química Orgánica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pabellón 2, Buenos Aires C1428EGA, Argentina.,Unidad de Microanálisis y Métodos Físicos Aplicados a Química Orgánica (UMYMFOR), Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires. CONICET, Buenos Aires C1428EGA, Argentina
| | - Marcelo A Marti
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pabellón 2, Buenos Aires C1428EGA, Argentina.,Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN), Universidad de Buenos Aires, Ciudad Universitaria, Pabellón 2, Buenos Aires C1428EGA, Argentina
| | - Maria M Garcia-Alai
- European Molecular Biology Laboratory Hamburg, Notkestrasse 85, Hamburg D-22607, Germany
| | - Adrián G Turjanski
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pabellón 2, Buenos Aires C1428EGA, Argentina.,Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN), Universidad de Buenos Aires, Ciudad Universitaria, Pabellón 2, Buenos Aires C1428EGA, Argentina
| |
Collapse
|
6
|
Computational Design of Inhibitors Targeting the Catalytic β Subunit of Escherichia coli FOF1-ATP Synthase. Antibiotics (Basel) 2022; 11:antibiotics11050557. [PMID: 35625201 PMCID: PMC9138118 DOI: 10.3390/antibiotics11050557] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2022] [Revised: 03/22/2022] [Accepted: 03/25/2022] [Indexed: 01/27/2023] Open
Abstract
With the uncontrolled growth of multidrug-resistant bacteria, there is an urgent need to search for new therapeutic targets, to develop drugs with novel modes of bactericidal action. FoF1-ATP synthase plays a crucial role in bacterial bioenergetic processes, and it has emerged as an attractive antimicrobial target, validated by the pharmaceutical approval of an inhibitor to treat multidrug-resistant tuberculosis. In this work, we aimed to design, through two types of in silico strategies, new allosteric inhibitors of the ATP synthase, by targeting the catalytic β subunit, a centerpiece in communication between rotor subunits and catalytic sites, to drive the rotary mechanism. As a model system, we used the F1 sector of Escherichia coli, a bacterium included in the priority list of multidrug-resistant pathogens. Drug-like molecules and an IF1-derived peptide, designed through molecular dynamics simulations and sequence mining approaches, respectively, exhibited in vitro micromolar inhibitor potency against F1. An analysis of bacterial and Mammalia sequences of the key structural helix-turn-turn motif of the C-terminal domain of the β subunit revealed highly and moderately conserved positions that could be exploited for the development of new species-specific allosteric inhibitors. To our knowledge, these inhibitors are the first binders computationally designed against the catalytic subunit of FOF1-ATP synthase.
Collapse
|
7
|
Alvarez-Garcia D, Schmidtke P, Cubero E, Barril X. Extracting Atomic Contributions to Binding Free Energy Using Molecular Dynamics Simulations with Mixed Solvents (MDmix). Curr Drug Discov Technol 2022; 19:62-68. [PMID: 34951392 PMCID: PMC9906626 DOI: 10.2174/1570163819666211223162829] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2021] [Revised: 09/28/2021] [Accepted: 10/05/2021] [Indexed: 11/22/2022]
Abstract
BACKGROUND Mixed solvents MD (MDmix) simulations have proved to be a useful and increasingly accepted technique with several applications in structure-based drug discovery. One of the assumptions behind the methodology is the transferability of free energy values from the simulated cosolvent molecules to larger drug-like molecules. However, the binding free energy maps (ΔGbind) calculated for the different moieties of the cosolvent molecules (e.g. a hydroxyl map for the ethanol) are largely influenced by the rest of the solvent molecule and do not reflect the intrinsic affinity of the moiety in question. As such, they are hardly transferable to different molecules. METHOD To achieve transferable energies, we present here a method for decomposing the molecular binding free energy into accurate atomic contributions. RESULT We demonstrate with two qualitative visual examples how the corrected energy maps better match known binding hotspots and how they can reveal hidden hotspots with actual drug design potential. CONCLUSION Atomic decomposition of binding free energies derived from MDmix simulations provides transferable and quantitative binding free energy maps.
Collapse
Affiliation(s)
- Daniel Alvarez-Garcia
- Gain Therapeutics, Parc Cientific de Barcelona, Baldiri Reixac 10, 08029 Barcelona, Spain
| | - Peter Schmidtke
- Facultat de Farmacia, Universitat de Barcelona, Av. Joan XXIII 27-31, 08028 Barcelona, Spain;,Current address: Discngine, 79 Avenue Ledru Rollin, 75012 Paris, France;
| | - Elena Cubero
- Gain Therapeutics, Parc Cientific de Barcelona, Baldiri Reixac 10, 08029 Barcelona, Spain
| | - Xavier Barril
- Gain Therapeutics, Parc Cientific de Barcelona, Baldiri Reixac 10, 08029 Barcelona, Spain;,Facultat de Farmacia, Universitat de Barcelona, Av. Joan XXIII 27-31, 08028 Barcelona, Spain;,Catalan Institution for Research and Advanced Studies (ICREA), Passeig Lluis Companys 23, 08010 Barcelona, Spain,Address correspondence to this author at the Gain Therapeutics, Parc Cientific de Barcelona, Baldiri Reixac 10, 08029 Barcelona, Spain; E-mail:
| |
Collapse
|
8
|
Burastero O, Cabrera M, Lopez ED, Defelipe LA, Arcon JP, Durán R, Marti MA, Turjanski AG. Specificity and Reactivity of Mycobacterium tuberculosis Serine/Threonine Kinases PknG and PknB. J Chem Inf Model 2022; 62:1723-1733. [PMID: 35319884 DOI: 10.1021/acs.jcim.1c01358] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
Mycobacterium tuberculosis (Mtb), the causative agent of Tuberculosis, has 11 eukaryotic-like serine/threonine protein kinases, which play essential roles in cell growth, signal transduction, and pathogenesis. Protein kinase G (PknG) regulates the carbon and nitrogen metabolism by phosphorylation of the glycogen accumulation regulator (GarA) protein at Thr21. Protein kinase B (PknB) is involved in cell wall synthesis and cell shape, as well as phosphorylates GarA but at Thr22. While PknG seems to be constitutively activated and recognition of GarA requires phosphorylation in its unstructured tail, PknB activation is triggered by phosphorylation of its activation loop, which allows binding of the forkhead-associated domain of GarA. In the present work, we used molecular dynamics and quantum-mechanics/molecular mechanics simulations of the catalytically competent complex and kinase activity assays to understand PknG/PknB specificity and reactivity toward GarA. Two hydrophobic residues in GarA, Val24 and Phe25, seem essential for PknG binding and allow specificity for Thr21 phosphorylation. On the other hand, phosphorylated residues in PknB bind Arg26 in GarA and regulate its specificity for Thr22. We also provide a detailed analysis of the free energy profile for the phospho-transfer reaction and show why PknG has a constitutively active conformation not requiring priming phosphorylation in contrast to PknB. Our results provide new insights into these two key enzymes relevant for Mtb and the mechanisms of serine/threonine phosphorylation in bacteria.
Collapse
Affiliation(s)
- Osvaldo Burastero
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pabellón 2, C1428EGA Buenos Aires, Argentina.,Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN), Universidad de Buenos Aires, Ciudad Universitaria, Pabellón 2, C1428EGA Buenos Aires, Argentina
| | - Marisol Cabrera
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pabellón 2, C1428EGA Buenos Aires, Argentina.,Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN), Universidad de Buenos Aires, Ciudad Universitaria, Pabellón 2, C1428EGA Buenos Aires, Argentina
| | - Elias D Lopez
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pabellón 2, C1428EGA Buenos Aires, Argentina.,Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN), Universidad de Buenos Aires, Ciudad Universitaria, Pabellón 2, C1428EGA Buenos Aires, Argentina
| | - Lucas A Defelipe
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pabellón 2, C1428EGA Buenos Aires, Argentina.,Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN), Universidad de Buenos Aires, Ciudad Universitaria, Pabellón 2, C1428EGA Buenos Aires, Argentina
| | - Juan Pablo Arcon
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pabellón 2, C1428EGA Buenos Aires, Argentina.,Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN), Universidad de Buenos Aires, Ciudad Universitaria, Pabellón 2, C1428EGA Buenos Aires, Argentina
| | - Rosario Durán
- Institut Pasteur de Montevideo, 11400 Montevideo, Uruguay.,Instituto de Investigaciones BiológicasClemente Estable, 11600 Montevideo, Uruguay
| | - Marcelo A Marti
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pabellón 2, C1428EGA Buenos Aires, Argentina.,Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN), Universidad de Buenos Aires, Ciudad Universitaria, Pabellón 2, C1428EGA Buenos Aires, Argentina
| | - Adrian G Turjanski
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pabellón 2, C1428EGA Buenos Aires, Argentina.,Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN), Universidad de Buenos Aires, Ciudad Universitaria, Pabellón 2, C1428EGA Buenos Aires, Argentina
| |
Collapse
|
9
|
Goel H, Hazel A, Yu W, Jo S, MacKerell AD. Application of Site-Identification by Ligand Competitive Saturation in Computer-Aided Drug Design. NEW J CHEM 2022; 46:919-932. [PMID: 35210743 PMCID: PMC8863107 DOI: 10.1039/d1nj04028f] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
Site Identification by Ligand Competitive Saturation (SILCS) is a molecular simulation approach that uses diverse small solutes in aqueous solution to obtain functional group affinity patterns of a protein or other macromolecule. This involves employing a combined Grand Canonical Monte Carlo (GCMC)-molecular dynamics (MD) method to sample the full 3D space of the protein, including deep binding pockets and interior cavities from which functional group free energy maps (FragMaps) are obtained. The information content in the maps, which include contributions from protein flexibilty and both protein and functional group desolvation contributions, can be used in many aspects of the drug discovery process. These include identification of novel ligand binding pockets, including allosteric sites, pharmacophore modeling, prediction of relative protein-ligand binding affinities for database screening and lead optimization efforts, evaluation of protein-protein interactions as well as in the formulation of biologics-based drugs including monoclonal antibodies. The present article summarizes the various tools developed in the context of the SILCS methodology and their utility in computer-aided drug design (CADD) applications, showing how the SILCS toolset can improve the drug-development process on a number of fronts with respect to both accuracy and throughput representing a new avenue of CADD applications.
Collapse
Affiliation(s)
- Himanshu Goel
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, 20, Penn St. Baltimore, Maryland 21201, United States
| | - Anthony Hazel
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, 20, Penn St. Baltimore, Maryland 21201, United States
| | - Wenbo Yu
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, 20, Penn St. Baltimore, Maryland 21201, United States
| | - Sunhwan Jo
- SilcsBio LLC, 1100 Wicomico St. Suite 323, Baltimore, MD, 21230, United States
| | - Alexander D. MacKerell
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, 20, Penn St. Baltimore, Maryland 21201, United States., SilcsBio LLC, 1100 Wicomico St. Suite 323, Baltimore, MD, 21230, United States.,, Tel: 410-706-7442, Fax: 410-706-5017
| |
Collapse
|
10
|
Szabó PB, Sabanés Zariquiey F, Nogueira JJ. Cosolvent and Dynamic Effects in Binding Pocket Search by Docking Simulations. J Chem Inf Model 2021; 61:5508-5523. [PMID: 34730967 PMCID: PMC8659376 DOI: 10.1021/acs.jcim.1c00924] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2021] [Indexed: 11/30/2022]
Abstract
The lack of conformational sampling in virtual screening projects can lead to inefficient results because many of the potential drugs may not be able to bind to the target protein during the static docking simulations. Here, we performed ensemble docking for around 2000 United States Food and Drug Administration (FDA)-approved drugs with the RNA-dependent RNA polymerase (RdRp) protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) as a target. The representative protein structures were generated by clustering classical molecular dynamics trajectories, which were evolved using three solvent scenarios, namely, pure water, benzene/water and phenol/water mixtures. The introduction of dynamic effects in the theoretical model showed improvement in docking results in terms of the number of strong binders and binding sites in the protein. Some of the discovered pockets were found only for the cosolvent simulations, where the nonpolar probes induced local conformational changes in the protein that lead to the opening of transient pockets. In addition, the selection of the ligands based on a combination of the binding free energy and binding free energy gap between the best two poses for each ligand provided more suitable binders than the selection of ligands based solely on one of the criteria. The application of cosolvent molecular dynamics to enhance the sampling of the configurational space is expected to improve the efficacy of virtual screening campaigns of future drug discovery projects.
Collapse
Affiliation(s)
- P. Bernát Szabó
- Department
of Chemistry, KU Leuven, Celestijnenlaan 200F, 3001 Leuven, Belgium
- Department
of Chemistry, Universidad Autónoma
de Madrid, Calle Francisco Tomás y Valiente, 7, 28049 Madrid, Spain
| | | | - Juan J. Nogueira
- Department
of Chemistry, Universidad Autónoma
de Madrid, Calle Francisco Tomás y Valiente, 7, 28049 Madrid, Spain
- IADCHEM,
Institute for Advanced Research in Chemistry, Universidad Autónoma de Madrid, Calle Francisco Tomás y Valiente, 7, 28049 Madrid, Spain
| |
Collapse
|
11
|
Yanagisawa K, Moriwaki Y, Terada T, Shimizu K. EXPRORER: Rational Cosolvent Set Construction Method for Cosolvent Molecular Dynamics Using Large-Scale Computation. J Chem Inf Model 2021; 61:2744-2753. [PMID: 34061535 DOI: 10.1021/acs.jcim.1c00134] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Cosolvent molecular dynamics (CMD) simulations involve an MD simulation of a protein in the presence of explicit water molecules mixed with cosolvent molecules to perform hotspot detection, binding site identification, and binding energy estimation, while other existing methods (e.g., MixMD, SILCS, and MDmix) utilize small molecules that represent functional groups of compounds. However, the cosolvent selections employed in these methods differ and there are only a few cosolvents that are commonly used in these methods. In this study, we proposed a systematic method for constructing a set of cosolvents for drug discovery, termed the EXtended PRObes set construction by REpresentative Retrieval (EXPRORER). First, we extracted typical substructures from FDA-approved drugs, generated 138 cosolvent structures, and for each cosolvent molecule, we conducted CMD simulations to generate a spatial probability distribution map of cosolvent atoms (PMAP). Analyses of PMAP similarity revealed that a cosolvent pair with a PMAP similarity greater than 0.70-0.75 shared similar structural features. We present a method for the construction of a cosolvent subset that satisfies a similarity threshold for all cosolvents, and we tested the constructed sets for four proteins. To our knowledge, this is the first study to include a systematic proposal for cosolvent set construction, and thus, the EXPRORER cosolvents will provide deeper insights into ligand binding sites of various proteins.
Collapse
Affiliation(s)
- Keisuke Yanagisawa
- Department of Computer Science, School of Computing, Tokyo Institute of Technology, Tokyo 152-8550, Japan.,Department of Biotechnology, Graduate School of Agricultural and Life Sciences, The University of Tokyo, Tokyo 113-8657, Japan
| | - Yoshitaka Moriwaki
- Agricultural Bioinformatics Research Unit, Graduate School of Agricultural and Life Sciences, The University of Tokyo, Tokyo 113-8657, Japan
| | - Tohru Terada
- Department of Biotechnology, Graduate School of Agricultural and Life Sciences, The University of Tokyo, Tokyo 113-8657, Japan.,Collaborative Research Institute for Innovative Microbiology, The University of Tokyo, Tokyo 113-8657, Japan.,Agricultural Bioinformatics Research Unit, Graduate School of Agricultural and Life Sciences, The University of Tokyo, Tokyo 113-8657, Japan
| | - Kentaro Shimizu
- Department of Biotechnology, Graduate School of Agricultural and Life Sciences, The University of Tokyo, Tokyo 113-8657, Japan.,Collaborative Research Institute for Innovative Microbiology, The University of Tokyo, Tokyo 113-8657, Japan.,Agricultural Bioinformatics Research Unit, Graduate School of Agricultural and Life Sciences, The University of Tokyo, Tokyo 113-8657, Japan
| |
Collapse
|
12
|
Arcon JP, Turjanski AG, Martí MA, Forli S. Biased Docking for Protein-Ligand Pose Prediction. Methods Mol Biol 2021; 2266:39-72. [PMID: 33759120 PMCID: PMC10708986 DOI: 10.1007/978-1-0716-1209-5_3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2023]
Abstract
The interaction between a protein and its ligands is one of the basic and most important processes in biological chemistry. Docking methods aim to predict the molecular 3D structure of protein-ligand complexes starting from coordinates of the protein and the ligand separately. They are widely used in both industry and academia, especially in the context of drug development projects. AutoDock4 is one of the most popular docking tools and, as for any docking method, its performance is highly system dependent. Knowledge about specific protein-ligand interactions on a particular target can be used to successfully overcome this limitation. Here, we describe how to apply the AutoDock Bias protocol, a simple and elegant strategy that allows users to incorporate target-specific information through a modified scoring function that biases the ligand structure towards those poses (or conformations) that establish selected interactions. We discuss two examples using different bias sources. In the first, we show how to steer dockings towards interactions derived from crystal structures of the receptor with different ligands; in the second example, we define and apply hydrophobic biases derived from Molecular Dynamics simulations in mixed solvents. Finally, we discuss general concepts of biased docking, its performance in pose prediction, and virtual screening campaigns as well as other potential applications.
Collapse
Affiliation(s)
- Juan Pablo Arcon
- Departamento de Química Biológica e IQUIBICEN-UBA/CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Buenos Aires, Argentina.
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, Barcelona, Spain.
| | - Adrián G Turjanski
- Departamento de Química Biológica e IQUIBICEN-UBA/CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Buenos Aires, Argentina
| | - Marcelo A Martí
- Departamento de Química Biológica e IQUIBICEN-UBA/CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Buenos Aires, Argentina
| | - Stefano Forli
- Department of Integrative Structural and Computational Biology, Scripps Research, La Jolla, CA, USA.
| |
Collapse
|
13
|
Wakefield AE, Yueh C, Beglov D, Castilho MS, Kozakov D, Keserű GM, Whitty A, Vajda S. Benchmark Sets for Binding Hot Spot Identification in Fragment-Based Ligand Discovery. J Chem Inf Model 2020; 60:6612-6623. [PMID: 33291870 DOI: 10.1021/acs.jcim.0c00877] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Binding hot spots are regions of proteins that, due to their potentially high contribution to the binding free energy, have high propensity to bind small molecules. We present benchmark sets for testing computational methods for the identification of binding hot spots with emphasis on fragment-based ligand discovery. Each protein structure in the set binds a fragment, which is extended into larger ligands in other structures without substantial change in its binding mode. Structures of the same proteins without any bound ligand are also collected to form an unbound benchmark. We also discuss a set developed by Astex Pharmaceuticals for the validation of hot and warm spots for fragment binding. The set is based on the assumption that a fragment that occurs in diverse ligands in the same subpocket identifies a binding hot spot. Since this set includes only ligand-bound proteins, we added a set with unbound structures. All four sets were tested using FTMap, a computational analogue of fragment screening experiments to form a baseline for testing other prediction methods, and differences among the sets are discussed.
Collapse
Affiliation(s)
- Amanda E Wakefield
- Department of Biomedical Engineering, Boston University, Boston, Massachusetts 02215, United States.,Department of Chemistry, Boston University, Boston, Massachusetts 02215, United States
| | - Christine Yueh
- Acpharis Inc., Holliston, Massachusetts 01746, United States
| | - Dmitri Beglov
- Acpharis Inc., Holliston, Massachusetts 01746, United States
| | - Marcelo S Castilho
- Faculdade de Farmácia da Universidade Federal da Bahia, Bahia, Salvador, BA 40170-115, Brazil
| | - Dima Kozakov
- Department of Applied Mathematics and Statistics, Stony Brook University, Stony Brook, New York 11794, United States.,Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York 11794, United States
| | - György M Keserű
- Medicinal Chemistry Research Group, Research Centre for Natural Sciences, Magyar tudósok krt. 2, 1117 Budapest, Hungary
| | - Adrian Whitty
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, United States
| | - Sandor Vajda
- Department of Biomedical Engineering, Boston University, Boston, Massachusetts 02215, United States.,Department of Chemistry, Boston University, Boston, Massachusetts 02215, United States
| |
Collapse
|
14
|
Yoshidome T, Ikeguchi M, Ohta M. Comprehensive 3D-RISM analysis of the hydration of small molecule binding sites in ligand-free protein structures. J Comput Chem 2020; 41:2406-2419. [PMID: 32815201 PMCID: PMC7540010 DOI: 10.1002/jcc.26406] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2020] [Revised: 07/28/2020] [Accepted: 08/05/2020] [Indexed: 12/25/2022]
Abstract
Hydration is a critical factor in the ligand binding process. Herein, to examine the hydration states of ligand binding sites, the three‐dimensional distribution function for the water oxygen site, gO(r), is computed for 3,706 ligand‐free protein structures based on the corresponding small molecule–protein complexes using the 3D‐RISM theory. For crystallographic waters (CWs) close to the ligand, gO(r) reveals that several CWs are stabilized by interaction networks formed between the ligand, CW, and protein. Based on the gO(r) for the crystallographic binding pose of the ligand, hydrogen bond interactions are dominant in the highly hydrated regions while weak interactions such as CH‐O are dominant in the moderately hydrated regions. The polar heteroatoms of the ligand occupy the highly hydrated and moderately hydrated regions in the crystallographic (correct) and wrongly docked (incorrect) poses, respectively. Thus, the gO(r) of polar heteroatoms may be used to distinguish the correct binding poses.
Collapse
Affiliation(s)
- Takashi Yoshidome
- Department of Applied Physics, Graduate School of Engineering, Tohoku University, Sendai, Japan
| | - Mitsunori Ikeguchi
- Drug Development Data Intelligence Platform Group, Medical Science Innovation Hub Program, Cluster of Science, Technology and Innovation Hub, RIKEN, Yokohama, Japan.,Graduate School of Medical Life Science, Yokohama City University, Yokohama, Japan
| | - Masateru Ohta
- Drug Development Data Intelligence Platform Group, Medical Science Innovation Hub Program, Cluster of Science, Technology and Innovation Hub, RIKEN, Yokohama, Japan
| |
Collapse
|
15
|
Schaller D, Šribar D, Noonan T, Deng L, Nguyen TN, Pach S, Machalz D, Bermudez M, Wolber G. Next generation 3D pharmacophore modeling. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2020. [DOI: 10.1002/wcms.1468] [Citation(s) in RCA: 57] [Impact Index Per Article: 14.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- David Schaller
- Pharmaceutical and Medicinal Chemistry Freie Universität Berlin Berlin Germany
| | - Dora Šribar
- Pharmaceutical and Medicinal Chemistry Freie Universität Berlin Berlin Germany
| | - Theresa Noonan
- Pharmaceutical and Medicinal Chemistry Freie Universität Berlin Berlin Germany
| | - Lihua Deng
- Pharmaceutical and Medicinal Chemistry Freie Universität Berlin Berlin Germany
| | - Trung Ngoc Nguyen
- Pharmaceutical and Medicinal Chemistry Freie Universität Berlin Berlin Germany
| | - Szymon Pach
- Pharmaceutical and Medicinal Chemistry Freie Universität Berlin Berlin Germany
| | - David Machalz
- Pharmaceutical and Medicinal Chemistry Freie Universität Berlin Berlin Germany
| | - Marcel Bermudez
- Pharmaceutical and Medicinal Chemistry Freie Universität Berlin Berlin Germany
| | - Gerhard Wolber
- Pharmaceutical and Medicinal Chemistry Freie Universität Berlin Berlin Germany
| |
Collapse
|
16
|
Li J, Liu W, Song Y, Xia J. Improved method of structure-based virtual screening based on ensemble learning. RSC Adv 2020; 10:7609-7618. [PMID: 35492172 PMCID: PMC9049841 DOI: 10.1039/c9ra09211k] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2019] [Accepted: 01/10/2020] [Indexed: 01/19/2023] Open
Abstract
Virtual screening has become a successful alternative and complementary technique to experimental high-throughput screening technologies for drug design. Since the scoring function of docking software cannot predict binding affinity accurately, how to improve the hit rate remains a common issue in structure-based virtual screening. This paper proposed a target-specific virtual screening method based on ensemble learning named ENS-VS. In this method, protein–ligand interaction energy terms and structure vectors of the ligands were used as a combination descriptor. Support vector machine, decision tree and Fisher linear discriminant classifiers were integrated into ENS-VS for predicting the activity of the compounds. The results showed that the enrichment factor (EF) 1% of ENS-VS was 6 times higher than that of Autodock vina. Compared with the newest virtual screening method SIEVE-Score, the mean EF 1% and AUC of ENS-VS (mean EF 1% = 52.77, AUC = 0.982) were statistically significantly higher than those of SIEVE-Score (mean EF 1% = 42.64, AUC = 0.912) on DUD-E datasets; and the mean EF 1% and AUC of ENS-VS (mean EF 1% = 29.73, AUC = 0.793) were also higher than those of SIEVE-Score (mean EF 1% = 25.56, AUC = 0.765) on eight DEKOIS datasets. ENS-VS also showed significant improvements compared with other similar research. The source code is available at https://github.com/eddyblue/ENS-VS. Virtual screening has become a successful alternative and complementary technique to experimental high-throughput screening technologies for drug design. This paper proposed a target-specific virtual screening method based on ensemble learning named ENS-VS.![]()
Collapse
Affiliation(s)
- Jin Li
- College of Computer and Information Science, Southwest University Chongqing 400715 China.,Key Laboratory of Medical Electrophysiology of Ministry of Education, Medical Electrophysiological Key Laboratory of Sichuan Province, Institute of Cardiovascular Research, School of Medical Information and Engineering, Southwest Medical University Luzhou 646000 China
| | - WeiChao Liu
- Key Laboratory of Medical Electrophysiology of Ministry of Education, Medical Electrophysiological Key Laboratory of Sichuan Province, Institute of Cardiovascular Research, School of Medical Information and Engineering, Southwest Medical University Luzhou 646000 China
| | | | - JiYi Xia
- Key Laboratory of Medical Electrophysiology of Ministry of Education, Medical Electrophysiological Key Laboratory of Sichuan Province, Institute of Cardiovascular Research, School of Medical Information and Engineering, Southwest Medical University Luzhou 646000 China
| |
Collapse
|
17
|
Modenutti CP, Capurro JIB, Di Lella S, Martí MA. The Structural Biology of Galectin-Ligand Recognition: Current Advances in Modeling Tools, Protein Engineering, and Inhibitor Design. Front Chem 2019; 7:823. [PMID: 31850312 PMCID: PMC6902271 DOI: 10.3389/fchem.2019.00823] [Citation(s) in RCA: 70] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2019] [Accepted: 11/12/2019] [Indexed: 12/25/2022] Open
Abstract
Galectins (formerly known as “S-type lectins”) are a subfamily of soluble proteins that typically bind β-galactoside carbohydrates with high specificity. They are present in many forms of life, from nematodes and fungi to animals, where they perform a wide range of functions. Particularly in humans, different types of galectins have been described differing not only in their tissue expression but also in their cellular location, oligomerization, fold architecture and carbohydrate-binding affinity. This distinct yet sometimes overlapping distributions and physicochemical attributes make them responsible for a wide variety of both intra- and extracellular functions, including tremendous importance in immunity and disease. In this review, we aim to provide a general description of galectins most important structural features, with a special focus on the molecular determinants of their carbohydrate-recognition ability. For that purpose, we structurally compare the human galectins, in light of recent mutagenesis studies and novel X-ray structures. We also offer a detailed description on how to use the solvent structure surrounding the protein as a tool to get better predictions of galectin-carbohydrate complexes, with a potential application to the rational design of glycomimetic inhibitory compounds. Finally, using Gal-1 and Gal-3 as paramount examples, we review a series of recent advances in the development of engineered galectins and galectin inhibitors, aiming to dissect the structure-activity relationship through the description of their interaction at the molecular level.
Collapse
Affiliation(s)
- Carlos P Modenutti
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Buenos Aires, Argentina.,Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN), CONICET, Buenos Aires, Argentina
| | - Juan I Blanco Capurro
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Buenos Aires, Argentina.,Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN), CONICET, Buenos Aires, Argentina
| | - Santiago Di Lella
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Buenos Aires, Argentina.,Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN), CONICET, Buenos Aires, Argentina
| | - Marcelo A Martí
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Buenos Aires, Argentina.,Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN), CONICET, Buenos Aires, Argentina
| |
Collapse
|