1
|
Flachsenberg F, Ehrt C, Gutermuth T, Rarey M. Redocking the PDB. J Chem Inf Model 2024; 64:219-237. [PMID: 38108627 DOI: 10.1021/acs.jcim.3c01573] [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: 12/19/2023]
Abstract
Molecular docking is a standard technique in structure-based drug design (SBDD). It aims to predict the 3D structure of a small molecule in the binding site of a receptor (often a protein). Despite being a common technique, it often necessitates multiple tools and involves manual steps. Here, we present the JAMDA preprocessing and docking workflow that is easy to use and allows fully automated docking. We evaluate the JAMDA docking workflow on binding sites extracted from the complete PDB and derive key factors determining JAMDA's docking performance. With that, we try to remove most of the bias due to manual intervention and provide a realistic estimate of the redocking performance of our JAMDA preprocessing and docking workflow for any PDB structure. On this large PDBScan22 data set, our JAMDA workflow finds a pose with an RMSD of at most 2 Å to the crystal ligand on the top rank for 30.1% of the structures. When applying objective structure quality filters to the PDBScan22 data set, the success rate increases to 61.8%. Given the prepared structures from the JAMDA preprocessing pipeline, both JAMDA and the widely used AutoDock Vina perform comparably on this filtered data set (the PDBScan22-HQ data set).
Collapse
Affiliation(s)
- Florian Flachsenberg
- Universität Hamburg, ZBH - Center for Bioinformatics, Bundesstraße 43, 20146 Hamburg, Germany
| | - Christiane Ehrt
- Universität Hamburg, ZBH - Center for Bioinformatics, Bundesstraße 43, 20146 Hamburg, Germany
| | - Torben Gutermuth
- Universität Hamburg, ZBH - Center for Bioinformatics, Bundesstraße 43, 20146 Hamburg, Germany
| | - Matthias Rarey
- Universität Hamburg, ZBH - Center for Bioinformatics, Bundesstraße 43, 20146 Hamburg, Germany
| |
Collapse
|
2
|
Penner P, Guba W, Schmidt R, Meyder A, Stahl M, Rarey M. The Torsion Library: Semiautomated Improvement of Torsion Rules with SMARTScompare. J Chem Inf Model 2022; 62:1644-1653. [PMID: 35318851 DOI: 10.1021/acs.jcim.2c00043] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The Torsion Library is a collection of torsion motifs associated with angle distributions, derived from crystallographic databases. It is used in strain assessment, conformer generation, and geometry optimization. A hierarchical structure of expert curated SMARTS defines the chemical environments of rotatable bonds and associates these with preferred angles. SMARTS can be very complex and full of implications, which make them difficult to maintain manually. Recent developments in automatically comparing SMARTS patterns can be applied to the Torsion Library to ensure its correctness. We specifically discuss the implementation and the limits of such a procedure in the context of torsion motifs and show several examples of how the Torsion Library benefits from this. All automated changes are validated manually and then shown to have an effect on the angle distributions by correcting matching behavior. The corrected Torsion Library itself is available including both PDB as well as CSD histograms in the Supporting Information and can be used to evaluate rotatable bonds at https://torsions.zbh.uni-hamburg.de.
Collapse
Affiliation(s)
- Patrick Penner
- Universität Hamburg,ZBH - Center for Bioinformatics, Bundesstraße 43, 20146 Hamburg, Germany
| | - Wolfgang Guba
- Roche Pharma Research & Early Development, Roche Innovation Center Basel, F. Hoffmann-La Roche Ltd., CH-4070 Basel, Switzerland
| | - Robert Schmidt
- Universität Hamburg,ZBH - Center for Bioinformatics, Bundesstraße 43, 20146 Hamburg, Germany
| | - Agnes Meyder
- Universität Hamburg,ZBH - Center for Bioinformatics, Bundesstraße 43, 20146 Hamburg, Germany
| | - Martin Stahl
- Roche Pharma Research & Early Development, Roche Innovation Center Basel, F. Hoffmann-La Roche Ltd., CH-4070 Basel, Switzerland
| | - Matthias Rarey
- Universität Hamburg,ZBH - Center for Bioinformatics, Bundesstraße 43, 20146 Hamburg, Germany
| |
Collapse
|
3
|
Schöning-Stierand K, Diedrich K, Fährrolfes R, Flachsenberg F, Meyder A, Nittinger E, Steinegger R, Rarey M. ProteinsPlus: interactive analysis of protein-ligand binding interfaces. Nucleic Acids Res 2020; 48:W48-W53. [PMID: 32297936 PMCID: PMC7319454 DOI: 10.1093/nar/gkaa235] [Citation(s) in RCA: 121] [Impact Index Per Article: 30.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2020] [Revised: 03/19/2020] [Accepted: 04/14/2020] [Indexed: 01/22/2023] Open
Abstract
Due to the increasing amount of publicly available protein structures searching, enriching and investigating these data still poses a challenging task. The ProteinsPlus web service (https://proteins.plus) offers a broad range of tools addressing these challenges. The web interface to the tool collection focusing on protein–ligand interactions has been geared towards easy and intuitive access to a large variety of functionality for life scientists. Since our last publication, the ProteinsPlus web service has been extended by additional services as well as it has undergone substantial infrastructural improvements. A keyword search functionality was added on the start page of ProteinsPlus enabling users to work on structures without knowing their PDB code. The tool collection has been augmented by three tools: StructureProfiler validates ligands and active sites using selection criteria of well-established protein–ligand benchmark data sets, WarPP places water molecules in the ligand binding sites of a protein, and METALizer calculates, predicts and scores coordination geometries of metal ions based on surrounding complex atoms. Additionally, all tools provided by ProteinsPlus are available through a REST service enabling the automated integration in structure processing and modeling pipelines.
Collapse
Affiliation(s)
| | - Konrad Diedrich
- Universität Hamburg, ZBH - Center for Bioinformatics (ZBH), 20146 Hamburg, Germany
| | - Rainer Fährrolfes
- Universität Hamburg, ZBH - Center for Bioinformatics (ZBH), 20146 Hamburg, Germany
| | - Florian Flachsenberg
- Universität Hamburg, ZBH - Center for Bioinformatics (ZBH), 20146 Hamburg, Germany
| | - Agnes Meyder
- Universität Hamburg, ZBH - Center for Bioinformatics (ZBH), 20146 Hamburg, Germany
| | - Eva Nittinger
- Universität Hamburg, ZBH - Center for Bioinformatics (ZBH), 20146 Hamburg, Germany
| | - Ruben Steinegger
- Universität Hamburg, ZBH - Center for Bioinformatics (ZBH), 20146 Hamburg, Germany
| | - Matthias Rarey
- Universität Hamburg, ZBH - Center for Bioinformatics (ZBH), 20146 Hamburg, Germany
| |
Collapse
|
4
|
Friedrich NO, Flachsenberg F, Meyder A, Sommer K, Kirchmair J, Rarey M. Conformator: A Novel Method for the Generation of Conformer Ensembles. J Chem Inf Model 2019; 59:731-742. [PMID: 30747530 DOI: 10.1021/acs.jcim.8b00704] [Citation(s) in RCA: 44] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Computer-aided drug design methods such as docking, pharmacophore searching, 3D database searching, and the creation of 3D-QSAR models need conformational ensembles to handle the flexibility of small molecules. Here, we present Conformator, an accurate and effective knowledge-based algorithm for generating conformer ensembles. With 99.9% of all test molecules processed, Conformator stands out by its robustness with respect to input formats, molecular geometries, and the handling of macrocycles. With an extended set of rules for sampling torsion angles, a novel algorithm for macrocycle conformer generation, and a new clustering algorithm for the assembly of conformer ensembles, Conformator reaches a median minimum root-mean-square deviation (measured between protein-bound ligand conformations and ensembles of a maximum of 250 conformers) of 0.47 Å with no significant difference to the highest-ranked commercial algorithm OMEGA and significantly higher accuracy than seven free algorithms, including the RDKit DG algorithm. Conformator is freely available for noncommercial use and academic research.
Collapse
Affiliation(s)
- Nils-Ole Friedrich
- Center for Bioinformatics , Universität Hamburg , Bundesstrasse 43 , 20146 Hamburg , Germany
| | - Florian Flachsenberg
- Center for Bioinformatics , Universität Hamburg , Bundesstrasse 43 , 20146 Hamburg , Germany
| | - Agnes Meyder
- Center for Bioinformatics , Universität Hamburg , Bundesstrasse 43 , 20146 Hamburg , Germany
| | - Kai Sommer
- Center for Bioinformatics , Universität Hamburg , Bundesstrasse 43 , 20146 Hamburg , Germany
| | - Johannes Kirchmair
- Center for Bioinformatics , Universität Hamburg , Bundesstrasse 43 , 20146 Hamburg , Germany.,Department of Chemistry , University of Bergen , N-5020 Bergen , Norway.,Computational Biology Unit (CBU) , University of Bergen , N-5020 Bergen , Norway
| | - Matthias Rarey
- Center for Bioinformatics , Universität Hamburg , Bundesstrasse 43 , 20146 Hamburg , Germany
| |
Collapse
|