Docking-generated multiple ligand poses for bootstrapping bioactivity classifying Machine Learning: Repurposing covalent inhibitors for COVID-19-related TMPRSS2 as case study.
Ontology highlight
ABSTRACT:
SUBMITTER: Hatmal MM
PROVIDER: S-EPMC8373588 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:The identification of protein-ligand binding sites is critical to protein function annotation and drug discovery. The consensus algorithm COACH developed by us represents one of the most efficient approaches to protein-ligand binding sites prediction. One of the most commonly seen issues with the COACH prediction are the low quality of the predicted ligand-binding poses, which usually have severe steric clashes to the protein structure. Here, we present COACH-D, an enhanced version of COACH by utilizing molecular docking to refine the ligand-binding poses. The input to the COACH-D server is the amino acid sequence or the three-dimensional structure of a query protein. In addition, the users can also submit their own ligand of interest. For each job submission, the COACH algorithm is first used to predict the protein-ligand binding sites. The ligands from the users or the templates are then docked into the predicted binding pockets to build their complex structures. Blind tests show that the algorithm significantly outperforms other ligand-binding sites prediction methods. Benchmark tests show that the steric clashes between the ligand and the protein structures in the COACH models are reduced by 85% after molecular docking in COACH-D. The COACH-D server is freely available to all users at http://yanglab.nankai.edu.cn/COACH-D/.
| S-EPMC6030866 | biostudies-literature
Project description:A molecular dynamics-based protocol is proposed for finding and scoring protein-ligand binding poses. This protocol uses the recently developed reconnaissance metadynamics method, which employs a self-learning algorithm to construct a bias that pushes the system away from the kinetic traps where it would otherwise remain. The exploration of phase space with this algorithm is shown to be roughly six to eight times faster than unbiased molecular dynamics and is only limited by the time taken to diffuse about the surface of the protein. We apply this method to the well-studied trypsin-benzamidine system and show that we are able to refind all the poses obtained from a reference EADock blind docking calculation. These poses can be scored based on the length of time the system remains trapped in the pose. Alternatively, one can perform dimensionality reduction on the output trajectory and obtain a map of phase space that can be used in more expensive free-energy calculations.
| S-EPMC3325676 | biostudies-literature
Project description:Plasmid pLNBIV was used to overexpress the biosynthetic pathway of nucleoside-diphosphate (NDP)-activated l-digitoxose in the mithramycin producer Streptomyces argillaceus. This led to a "flooding" of the biosynthetic pathway of the antitumor drug mithramycin (MTM) with NDP-activated deoxysugars, which do not normally occur in the pathway, and consequently to the production of the four new mithramycin derivatives 1- 4 with altered saccharide patterns. Their structures reflect that NDP sugars produced by pLNBIV, namely, l-digitoxose and its biosynthetic intermediates, influenced the glycosyl transfer to positions B, D, and E, while positions A and C remained unaffected. All four new structures have unique, previously not found sugar decoration patterns, which arise from either overcoming the substrate specificity or inhibition of certain glycosyltransferases (GTs) of the MTM pathway with the foreign NDP sugars expressed by pLNBIV. An apoptosis TUNEL (=terminal deoxynucleotidyl transferase-mediated dUTP nick end labeling) assay revealed that compounds 1 (demycarosyl-3D-beta- d-digitoxosyl-MTM) and 3 (deoliosyl-3C-beta- d-mycarosyl-MTM) show improved activity (64.8 +/- 2% and 50.3 +/- 2.5% induction of apoptosis, respectively) against the estrogen receptor (ER)-positive human breast cancer cell line MCF-7 compared with the parent drug MTM (37.8 +/- 2.5% induction of apoptosis). In addition, compounds 1 and 4 (3A-deolivosyl-MTM) show significant effects on the ER-negative human breast cancer cell line MDA-231 (63.6 +/- 2% and 12.6 +/- 2.5% induction of apoptosis, respectively), which is not inhibited by the parent drug MTM itself (2.6 +/- 1.5% induction of apoptosis), but for which chemotherapeutic agents are urgently needed.
| S-EPMC2442402 | biostudies-literature
Project description:Molecular docking algorithms suggest possible structures for molecular complexes. They are used to model biological function and to discover potential ligands. A present challenge for docking algorithms is the treatment of molecular flexibility. Here, the rigid body program, DOCK, is modified to allow it to rapidly fit multiple conformations of ligands. Conformations of a given molecule are pre-calculated in the same frame of reference, so that each conformer shares a common rigid fragment with all other conformations. The ligand conformers are then docked together, as an ensemble, into a receptor binding site. This takes advantage of the redundancy present in differing conformers of the same molecule. The algorithm was tested using three organic ligand protein systems and two protein-protein systems. Both the bound and unbound conformations of the receptors were used. The ligand ensemble method found conformations that resembled those determined in X-ray crystal structures (RMS values typically less than 1.5 A). To test the method's usefulness for inhibitor discovery, multi-compound and multi-conformer databases were screened for compounds known to bind to dihydrofolate reductase and compounds known to bind to thymidylate synthase. In both cases, known inhibitors and substrates were identified in conformations resembling those observed experimentally. The ligand ensemble method was 100-fold faster than docking a single conformation at a time and was able to screen a database of over 34 million conformations from 117,000 molecules in one to four CPU days on a workstation.
| S-EPMC2143983 | biostudies-other
Project description:We describe an innovative protocol for ab initio prediction of ligand crystallographic binding poses and highly effective analysis of large datasets generated for protein-ligand dynamics. We include a procedure for setup and performance of distributed molecular dynamics simulations on cloud computing architectures, a model for efficient analysis of simulation data, and a metric for evaluation of model convergence. We give accurate binding pose predictions for five ligands ranging in affinity from 7 nM to > 200??M for the immunophilin protein FKBP12, for expedited results in cases where experimental structures are difficult to produce. Our approach goes beyond single, low energy ligand poses to give quantitative kinetic information that can inform protein engineering and ligand design.
| S-EPMC4302315 | biostudies-literature
Project description:Electrophilic peptides that form an irreversible covalent bond with their target have great potential for binding targets that have been previously considered undruggable. However, the discovery of such peptides remains a challenge. Here, we present Rosetta CovPepDock, a computational pipeline for peptide docking that incorporates covalent binding between the peptide and a receptor cysteine. We applied CovPepDock retrospectively to a dataset of 115 disulfide-bound peptides and a dataset of 54 electrophilic peptides. It produced a top-five scoring, near-native model, in 89% and 100% of the cases when docking from the native conformation, and 20% and 90% when docking from an extended peptide conformation, respectively. In addition, we developed a protocol for designing electrophilic peptide binders based on known non-covalent binders or protein-protein interfaces. We identified 7154 peptide candidates in the PDB for application of this protocol. As a proof-of-concept we validated the protocol on the non-covalent complex of 14-3-3σ and YAP1 phosphopeptide. The protocol identified seven highly potent and selective irreversible peptide binders. The predicted binding mode of one of the peptides was validated using X-ray crystallography. This case-study demonstrates the utility and impact of CovPepDock. It suggests that many new electrophilic peptide binders can be rapidly discovered, with significant potential as therapeutic molecules and chemical probes.
| S-EPMC8372624 | biostudies-literature
Project description:Here we present a virtual docking screen of 1648 commercially available covalent fragments, and identified covalent inhibitors of cysteine protease cathepsin L. These inhibitors did not inhibit closely related protease cathepsin B. Thus, we have established virtual docking of covalent fragments as an approach to discover covalent enzyme inhibitors.
| S-EPMC6319256 | biostudies-literature
Project description:In molecular docking, it is challenging to develop a scoring function that is accurate to conduct high-throughput screenings. Most scoring functions implemented in popular docking software packages were developed with many approximations for computational efficiency, which sacrifices the accuracy of prediction. With advanced technology and powerful computational hardware nowadays, it is feasible to use rigorous scoring functions, such as molecular mechanics/Poisson Boltzmann surface area (MM/PBSA) and molecular mechanics/generalized Born surface area (MM/GBSA) in molecular docking studies. Here, we systematically investigated the performance of MM/PBSA and MM/GBSA to identify the correct binding conformations and predict the binding free energies for 98 protein-ligand complexes. Comparison studies showed that MM/GBSA (69.4%) outperformed MM/PBSA (45.5%) and many popular scoring functions to identify the correct binding conformations. Moreover, we found that molecular dynamics simulations are necessary for some systems to identify the correct binding conformations. Based on our results, we proposed the guideline for MM/GBSA to predict the binding conformations. We then tested the performance of MM/GBSA and MM/PBSA to reproduce the binding free energies of the 98 protein-ligand complexes. The best prediction of MM/GBSA model with internal dielectric constant 2.0, produced a Spearman's correlation coefficient of 0.66, which is better than MM/PBSA (0.49) and almost all scoring functions used in molecular docking. In summary, MM/GBSA performs well for both binding pose predictions and binding free-energy estimations and is efficient to re-score the top-hit poses produced by other less-accurate scoring functions.
| S-EPMC3043139 | biostudies-literature
Project description:BackgroundProtein-protein docking, which aims to predict the structure of a protein-protein complex from its unbound components, remains an unresolved challenge in structural bioinformatics. An important step is the ranking of docked poses using a scoring function, for which many methods have been developed. There is a need to explore the differences and commonalities of these methods with each other, as well as with functions developed in the fields of molecular dynamics and homology modelling.ResultsWe present an evaluation of 115 scoring functions on an unbound docking decoy benchmark covering 118 complexes for which a near-native solution can be found, yielding top 10 success rates of up to 58%. Hierarchical clustering is performed, so as to group together functions which identify near-natives in similar subsets of complexes. Three set theoretic approaches are used to identify pairs of scoring functions capable of correctly scoring different complexes. This shows that functions in different clusters capture different aspects of binding and are likely to work together synergistically.ConclusionsAll functions designed specifically for docking perform well, indicating that functions are transferable between sampling methods. We also identify promising methods from the field of homology modelling. Further, differential success rates by docking difficulty and solution quality suggest a need for flexibility-dependent scoring. Investigating pairs of scoring functions, the set theoretic measures identify known scoring strategies as well as a number of novel approaches, indicating promising augmentations of traditional scoring methods. Such augmentation and parameter combination strategies are discussed in the context of the learning-to-rank paradigm.
| S-EPMC3850738 | biostudies-literature
Project description: Not available
| S-EPMC10964654 | biostudies-literature