Browse
Submit Data
Databases
API
Help

Dataset Information

26 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Assessing the performance of docking scoring function, FEP, MM-GBSA, and QM/MM-GBSA approaches on a series of PLK1 inhibitors.

ABSTRACT: Over-expressed polo-like kinases 1, a key regulator of cell mitosis, is associated with carcinogenesis and poor prognosis. It is very necessary to develop a reliable computational affinity prediction protocol targeting PLK1. In this study, the performance of different docking scoring function, free energy perturbation, MM-GBSA and QM/MM-GBSA were evaluated. The ranking capability of FEP is the best with rs = 0.854. However, the rs obtained from MM-GBSA can reach 0.767, which requires only about one-eighth of the simulation time of FEP. As for the sampling method, single long molecular dynamics (SLMD) surpass the multiple short molecular dynamics (MSMD) in ranking of the 20 congeneric compounds by about 0.1 in rs. In addition, ligands treated by QM can significantly improve the ranking performance. As for the docking scoring functions, a force field-based scoring function is more suitable for ranking congeneric compounds.

SUBMITTER: Pu C

PROVIDER: S-EPMC6072424 | biostudies-literature | 2017 Jul

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Publications

Assessing the performance of docking scoring function, FEP, MM-GBSA, and QM/MM-GBSA approaches on a series of PLK1 inhibitors.

Pu Chunlan C Yan Guoyi G Shi Jianyou J Li Rui R

MedChemComm 20170522 7

Over-expressed polo-like kinases 1, a key regulator of cell mitosis, is associated with carcinogenesis and poor prognosis. It is very necessary to develop a reliable computational affinity prediction protocol targeting PLK1. In this study, the performance of different docking scoring function, free energy perturbation, MM-GBSA and QM/MM-GBSA were evaluated. The ranking capability of FEP is the best with rs = 0.854. However, the rs obtained from MM-GBSA can rea ...[more]

PMID: 30108856

Similar Datasets

Comparison of radii sets, entropy, QM methods, and sampling on MM-PBSA, MM-GBSA, and QM/MM-GBSA ligand binding energies of F. tularensis enoyl-ACP reductase (FabI).

Project description:To validate a method for predicting the binding affinities of FabI inhibitors, three implicit solvent methods, MM-PBSA, MM-GBSA, and QM/MM-GBSA were carefully compared using 16 benzimidazole inhibitors in complex with Francisella tularensis FabI. The data suggests that the prediction results are sensitive to radii sets, GB methods, QM Hamiltonians, sampling protocols, and simulation length, if only one simulation trajectory is used for each ligand. In this case, QM/MM-GBSA using 6 ns MD simulation trajectories together with GB(neck2) , PM3, and the mbondi2 radii set, generate the closest agreement with experimental values (r(2) = 0.88). However, if the three implicit solvent methods are averaged from six 1 ns MD simulations for each ligand (called "multiple independent sampling"), the prediction results are relatively insensitive to all the tested parameters. Moreover, MM/GBSA together with GB(HCT) and mbondi, using 600 frames extracted evenly from six 0.25 ns MD simulations, can also provide accurate prediction to experimental values (r(2) = 0.84). Therefore, the multiple independent sampling method can be more efficient than a single, long simulation method. Since future scaffold expansions may significantly change the benzimidazole's physiochemical properties (charges, etc.) and possibly binding modes, which may affect the sensitivities of various parameters, the relatively insensitive "multiple independent sampling method" may avoid the need of an entirely new validation study. Moreover, due to large fluctuating entropy values, (QM/)MM-P(G)BSA were limited to inhibitors' relative affinity prediction, but not the absolute affinity. The developed protocol will support an ongoing benzimidazole lead optimization program.

| S-EPMC4688044 | biostudies-literature

Binding-affinity predictions of HSP90 in the D3R Grand Challenge 2015 with docking, MM/GBSA, QM/MM, and free-energy simulations.

Project description:We have estimated the binding affinity of three sets of ligands of the heat-shock protein 90 in the D3R grand challenge blind test competition. We have employed four different methods, based on five different crystal structures: first, we docked the ligands to the proteins with induced-fit docking with the Glide software and calculated binding affinities with three energy functions. Second, the docked structures were minimised in a continuum solvent and binding affinities were calculated with the MM/GBSA method (molecular mechanics combined with generalised Born and solvent-accessible surface area solvation). Third, the docked structures were re-optimised by combined quantum mechanics and molecular mechanics (QM/MM) calculations. Then, interaction energies were calculated with quantum mechanical calculations employing 970-1160 atoms in a continuum solvent, combined with energy corrections for dispersion, zero-point energy and entropy, ligand distortion, ligand solvation, and an increase of the basis set to quadruple-zeta quality. Fourth, relative binding affinities were estimated by free-energy simulations, using the multi-state Bennett acceptance-ratio approach. Unfortunately, the results were varying and rather poor, with only one calculation giving a correlation to the experimental affinities larger than 0.7, and with no consistent difference in the quality of the predictions from the various methods. For one set of ligands, the results could be strongly improved (after experimental data were revealed) if it was recognised that one of the ligands displaced one or two water molecules. For the other two sets, the problem is probably that the ligands bind in different modes than in the crystal structures employed or that the conformation of the ligand-binding site or the whole protein changes.

| S-EPMC5078160 | biostudies-literature

Delineating Pixantrone Maleate's adroit activity against cervical cancer proteins through multitargeted docking-based MM\GBSA, QM-DFT and MD simulation.

Project description:Cervical cancer poses a substantial worldwide health challenge, especially in low- and middle-income nations, caused by high-risk types of human papillomavirus. It accounted for a significant percentage of cancer-related deaths among women, particularly in areas with limited healthcare resources, necessitating innovative therapeutic approaches, and single-targeted studies have produced significant results, with a considerable chance of developing resistance. Therefore, the multitargeted studies can work as a beacon of hope. This study is focused on performing the multitargeted molecular docking of FDA-approved drugs with the three crucial proteins TBK1, DNA polymerase epsilon, and integrin α-V β-8 of cervical cancer. The docking studies using multisampling algorithms HTVS, SP, and XP reveal Pixantrone Maleate (DB06193) as a multitargeted inhibitor with docking scores of -8.147, -8.206 and -7.31 Kcal/mol and pose filtration with MM\GBSA computations with scores -40.55, -33.67, and -37.64 Kcal/mol. We also have performed QM-based DFT and pharmacokinetics studies of the compound and compared it with the standard values, which results in the compound being entirely suitable against cervical cancer proteins. The interaction fingerprints have revealed that PHE, VAL, SER and ALA are the residues among most interactions. We also explore the stability of the multitargeted potential of Pixantrone Maleate through 100ns MD simulations and investigate the RMSD, RMSF and intermolecular interactions between all three proteins-ligand complexes. All computational studies favour Pixantrone Maleate as a multitargeted inhibitor of the TBK1, DNA polymerase epsilon, and integrin α-V β-8 and can be validated experimentally before use.

| S-EPMC10723688 | biostudies-literature

Delineated 3-1-BenCarMethInYlPro-Phosphonic Acid's Adroit Activity against Lung Cancer through Multitargeted Docking, MM\GBSA, QM-DFT and Multiscale Simulations.

Project description:Lung cancer is a pervasive and challenging disease with limited treatment options, with global health challenges often present with complex molecular profiles necessitating the exploration of innovative therapeutic strategies. Single-target drugs have shown limited success due to the heterogeneity of this disease. Multitargeted drug designing is imperative to combat this complexity by simultaneously targeting multiple target proteins and pathways, which can enhance treatment efficacy and overcome resistance by addressing the dynamic nature of the disease and stopping tumour growth and spread. In this study, we performed the molecular docking studies of Drug Bank compounds with a multitargeted approach against crucial proteins of lung cancer such as heat shock protein 5 (BIP/GRP78) ATPase, myosin 9B RhoGAP, EYA2 phosphatase inhibitor, RSK4 N-terminal kinase, and collapsin response mediator protein-1 (CRMP-1) using HTVS, SP with XP algorithms, and poses were filtered using MM\GBSA which identified [3-(1-Benzyl-3-Carbamoylmethyl-2-Methyl-1h-Indol-5-Yloxy)-Propyl-]-Phosphonic Acid (3-1-BenCarMethIn YlPro-Phosphonic Acid) (DB02504) as multitargeted drug candidate with docking and MM\GBSA score ranges from -5.83 to -10.66 and -7.56 to -50.14 Kcal/mol, respectively. Further, the pharmacokinetic and QM-based DFT studies have shown complete acceptance results, and interaction fingerprinting reveals that ILE, GLY, VAL, TYR, LEU, and GLN were among the most interacting residues. The 100 ns MD simulation in the SPC water model with NPT ensemble showed stable performance with deviation and fluctuations <2 Å with huge interactions, making it a promising multitargeted drug candidate; however, experimental studies are needed before use.

| S-EPMC10779231 | biostudies-literature

Insights into the Interaction between Polyphenols and β-Lactoglobulin through Molecular Docking, MD Simulation, and QM/MM Approaches.

Project description:In this work, we have explored the interaction of three different polyphenols with the food protein β-lactoglobulin. Antioxidant activities of polyphenols are influenced by complexation with the protein. However, studies have shown that polyphenols after complexation with the protein can be more beneficial due to enhanced antioxidant activities. We have carried out molecular docking, molecular dynamics (MD) simulation, and quantum mechanics/molecular mechanics (QM/MM) studies on the three different protein-polyphenol complexes. We have found from molecular docking studies that apigenin binds in the internal cavity, luteolin binds at the mouth of the cavity, and eriodictyol binds outside the cavity of the protein. Docking studies have also provided binding free energy and inhibition constant values that showed that eriodictyol and apigenin exhibit better binding interactions with the protein than luteolin. For eriodictyol and luteolin, van der Waals, hydrophobic, and hydrogen bonding interactions are the main interacting forces, whereas for apigenin, hydrophobic and van der Waals interactions play major roles. We have calculated the root mean square deviation (RMSD), root mean square fluctuations (RMSF), solvent-accessible surface area (SASA), interaction energies, and hydrogen bonds of the protein-polyphenol complexes. Results show that the protein-eriodictyol complex is more stable than the other complexes. We have performed ONIOM calculations to study the antioxidant properties of the polyphenols. We have found that apigenin and luteolin act as better antioxidants than eriodictyol does on complexation with the protein, which is consistent with the results obtained from MD simulations.

| S-EPMC9280950 | biostudies-literature

A Mixed QM/MM Scoring Function to Predict Protein-Ligand Binding Affinity.

Project description:Computational methods for predicting protein-ligand binding free energy continue to be popular as a potential cost-cutting method in the drug discovery process. However, accurate predictions are often difficult to make as estimates must be made for certain electronic and entropic terms in conventional force field based scoring functions. Mixed quantum mechanics/molecular mechanics (QM/MM) methods allow electronic effects for a small region of the protein to be calculated, treating the remaining atoms as a fixed charge background for the active site. Such a semi-empirical QM/MM scoring function has been implemented in AMBER using DivCon and tested on a set of 23 metalloprotein-ligand complexes, where QM/MM methods provide a particular advantage in the modeling of the metal ion. The binding affinity of this set of proteins can be calculated with an R(2) of 0.64 and a standard deviation of 1.88 kcal/mol without fitting and 0.71 and a standard deviation of 1.69 kcal/mol with fitted weighting of the individual scoring terms. In this study we explore using various methods to calculate terms in the binding free energy equation, including entropy estimates and minimization standards. From these studies we found that using the rotational bond estimate to ligand entropy results in a reasonable R(2) of 0.63 without fitting. We also found that using the ESCF energy of the proteins without minimization resulted in an R(2) of 0.57, when using the rotatable bond entropy estimate.

| S-EPMC3017370 | biostudies-literature

A Comprehensive Docking and MM/GBSA Rescoring Study of Ligand Recognition upon Binding Antithrombin.

Project description:BackgroundA high-throughput virtual screening pipeline has been extended from single energetically minimized structure Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) rescoring to ensemble-average MM/GBSA rescoring. The correlation coefficient (R2) of calculated and experimental binding free energies for a series of antithrombin ligands has been improved from 0.36 to 0.69 when switching from the single-structure MM/GBSA rescoring to ensemble-average one. The electrostatic interactions in both solute and solvent are identified to play an important role in determining the binding free energy after the decomposition of the calculated binding free energy. The increasing negative charge of the compounds provides a more favorable electrostatic energy change but creates a higher penalty for the solvation free energy. Such a penalty is compensated by the electrostatic energy change, which results in a better binding affinity. A highly hydrophobic ligand is determined by the docking program to be a non-specific binder.ResultsOur results have demonstrated that it is very important to keep a few top poses for rescoring, if the binding is non-specific or the binding mode is not well determined by the docking calculation.

| S-EPMC5403970 | biostudies-literature

Study of the Differential Activity of Thrombin Inhibitors Using Docking, QSAR, Molecular Dynamics, and MM-GBSA.

Project description:Non-peptidic thrombin inhibitors (TIs; 177 compounds) with diverse groups at motifs P1 (such as oxyguanidine, amidinohydrazone, amidine, amidinopiperidine), P2 (such as cyanofluorophenylacetamide, 2-(2-chloro-6-fluorophenyl)acetamide), and P3 (such as phenylethyl, arylsulfonate groups) were studied using molecular modeling to analyze their interactions with S1, S2, and S3 subsites of the thrombin binding site. Firstly, a protocol combining docking and three dimensional quantitative structure-activity relationship was performed. We described the orientations and preferred active conformations of the studied inhibitors, and derived a predictive CoMSIA model including steric, donor hydrogen bond, and acceptor hydrogen bond fields. Secondly, the dynamic behaviors of some selected TIs (compounds 26, 133, 147, 149, 162, and 177 in this manuscript) that contain different molecular features and different activities were analyzed by creating the solvated models and using molecular dynamics (MD) simulations. We used the conformational structures derived from MD to accomplish binding free energetic calculations using MM-GBSA. With this analysis, we theorized about the effect of van der Waals contacts, electrostatic interactions and solvation in the potency of TIs. In general, the contents reported in this article help to understand the physical and chemical characteristics of thrombin-inhibitor complexes.

| S-EPMC4657979 | biostudies-literature

Assessing the performance of MM/PBSA and MM/GBSA methods. 8. Predicting binding free energies and poses of protein-RNA complexes.

Project description:Molecular docking provides a computationally efficient way to predict the atomic structural details of protein-RNA interactions (PRI), but accurate prediction of the three-dimensional structures and binding affinities for PRI is still notoriously difficult, partly due to the unreliability of the existing scoring functions for PRI. MM/PBSA and MM/GBSA are more theoretically rigorous than most scoring functions for protein-RNA docking, but their prediction performance for protein-RNA systems remains unclear. Here, we systemically evaluated the capability of MM/PBSA and MM/GBSA to predict the binding affinities and recognize the near-native binding structures for protein-RNA systems with different solvent models and interior dielectric constants (εin). For predicting the binding affinities, the predictions given by MM/GBSA based on the minimized structures in explicit solvent and the GBGBn1 model with εin = 2 yielded the highest correlation with the experimental data. Moreover, the MM/GBSA calculations based on the minimized structures in implicit solvent and the GBGBn1 model distinguished the near-native binding structures within the top 10 decoys for 117 out of the 148 protein-RNA systems (79.1%). This performance is better than all docking scoring functions studied here. Therefore, the MM/GBSA rescoring is an efficient way to improve the prediction capability of scoring functions for protein-RNA systems.

| S-EPMC6097651 | biostudies-literature

Probing Origin of Binding Difference of inhibitors to MDM2 and MDMX by Polarizable Molecular Dynamics Simulation and QM/MM-GBSA Calculation.

Project description:Binding abilities of current inhibitors to MDMX are weaker than to MDM2. Polarizable molecular dynamics simulations (MD) followed by Quantum mechanics/molecular mechanics generalized Born surface area (QM//MM-GBSA) calculations were performed to investigate the binding difference of inhibitors to MDM2 and MDMX. The predicted binding free energies not only agree well with the experimental results, but also show that the decrease in van der Walls interactions of inhibitors with MDMX relative to MDM2 is a main factor of weaker bindings of inhibitors to MDMX. The analyses of dihedral angles based on MD trajectories suggest that the closed conformation formed by the residues M53 and Y99 in MDMX leads to a potential steric clash with inhibitors and prevents inhibitors from arriving in the deep of MDMX binding cleft, which reduces the van der Waals contacts of inhibitors with M53, V92, P95 and L98. The calculated results using the residue-based free energy decomposition method further prove that the interaction strength of inhibitors with M53, V92, P95 and L98 from MDMX are obviously reduced compared to MDM2. We expect that this study can provide significant theoretical guidance for designs of potent dual inhibitors to block the p53-MDM2/MDMX interactions.

| S-EPMC4663504 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

Tweets

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data