scholarly journals Combining Molecular Dynamic Information and an Aspherical-Atom Data Bank in the Evaluation of the Electrostatic Interaction Energy in Multimeric Protein-Ligand Complex: A Case Study for HIV-1 Protease

Molecules ◽  
2021 ◽  
Vol 26 (13) ◽  
pp. 3872
Author(s):  
Prashant Kumar ◽  
Paulina Maria Dominiak

Computational analysis of protein–ligand interactions is of crucial importance for drug discovery. Assessment of ligand binding energy allows us to have a glimpse of the potential of a small organic molecule to be a ligand to the binding site of a protein target. Available scoring functions, such as in docking programs, all rely on equations that sum each type of protein–ligand interactions in order to predict the binding affinity. Most of the scoring functions consider electrostatic interactions involving the protein and the ligand. Electrostatic interactions constitute one of the most important part of total interactions between macromolecules. Unlike dispersion forces, they are highly directional and therefore dominate the nature of molecular packing in crystals and in biological complexes and contribute significantly to differences in inhibition strength among related enzyme inhibitors. In this study, complexes of HIV-1 protease with inhibitor molecules (JE-2147 and darunavir) were analyzed by using charge densities from the transferable aspherical-atom University at Buffalo Databank (UBDB). Moreover, we analyzed the electrostatic interaction energy for an ensemble of structures, using molecular dynamic simulations to highlight the main features of electrostatic interactions important for binding affinity.

2021 ◽  
Author(s):  
Prashant Kumar ◽  
Paulina Dominiak

<div> <div> <div> <p>Computational analysis of protein-ligand interactions is of crucial importance for drug discovery. Assessment of ligand binding energy allows us to have a glimpse on the potential of a small organic molecule to be a ligand to the binding site of a protein target. Available scoring functions such as in docking programs, we could say that they all rely on equations that sum each type of protein-ligand interactions to model the binding affinity. Most of the scoring functions consider electrostatic interactions involving the protein and the ligand. Electrostatic interactions contribute one of the most important part of total interaction energies between macromolecules, unlike dispersion forces they are highly directional and therefore dominate the nature of molecular packing in crystals and in biological complexes and contribute significantly to differences in inhibition strength among related enzyme inhibitors. In this paper, complexes of HIV-1 protease with inhibitor molecules (JE-2147 and Darunavir) have been analysed using charge densities from a transferable aspherical-atom data bank. Moreover, we analyse the electrostatic interaction energy for an ensemble of structures using molecular dynamic simulation to highlight the main features related to the importance of this interaction for binding affinity. </p> </div> </div> </div>


2021 ◽  
Author(s):  
Prashant Kumar ◽  
Paulina Dominiak

<div> <div> <div> <p>Computational analysis of protein-ligand interactions is of crucial importance for drug discovery. Assessment of ligand binding energy allows us to have a glimpse on the potential of a small organic molecule to be a ligand to the binding site of a protein target. Available scoring functions such as in docking programs, we could say that they all rely on equations that sum each type of protein-ligand interactions to model the binding affinity. Most of the scoring functions consider electrostatic interactions involving the protein and the ligand. Electrostatic interactions contribute one of the most important part of total interaction energies between macromolecules, unlike dispersion forces they are highly directional and therefore dominate the nature of molecular packing in crystals and in biological complexes and contribute significantly to differences in inhibition strength among related enzyme inhibitors. In this paper, complexes of HIV-1 protease with inhibitor molecules (JE-2147 and Darunavir) have been analysed using charge densities from a transferable aspherical-atom data bank. Moreover, we analyse the electrostatic interaction energy for an ensemble of structures using molecular dynamic simulation to highlight the main features related to the importance of this interaction for binding affinity. </p> </div> </div> </div>


2019 ◽  
Vol 13 (1) ◽  
pp. 40-49 ◽  
Author(s):  
Pedro Fong ◽  
Hong-Kong Wong

Background: DNA has been a pharmacological target for different types of treatment, such as antibiotics and chemotherapy agents, and is still a potential target in many drug discovery processes. However, most docking and scoring approaches were parameterised for protein-ligand interactions; their suitability for modelling DNA-ligand interactions is uncertain. Objective: This study investigated the performance of four scoring functions on DNA-ligand complexes. Material & Methods: Here, we explored the ability of four docking protocols and scoring functions to discriminate the native pose of 33 DNA-ligand complexes over a compiled set of 200 decoys for each DNA-ligand complexes. The four approaches were the AutoDock, ASP@GOLD, ChemScore@GOLD and GoldScore@GOLD. Results: Our results indicate that AutoDock performed the best when predicting binding mode and that ChemScore@GOLD achieved the best discriminative power. Rescoring of AutoDock-generated decoys with ChemScore@GOLD further enhanced their individual discriminative powers. All four approaches have no discriminative power in some DNA-ligand complexes, including both minor groove binders and intercalators. Conclusion: This study suggests that the evaluation for each DNA-ligand complex should be performed in order to obtain meaningful results for any drug discovery processes. Rescoring with different scoring functions can improve discriminative power.


2020 ◽  
Vol 17 (2) ◽  
pp. 233-247
Author(s):  
Krishna A. Gajjar ◽  
Anuradha K. Gajjar

Background: Pharmacophore mapping and molecular docking can be synergistically integrated to improve the drug design and discovery process. A rational strategy, combiphore approach, derived from the combined study of Structure and Ligand based pharmacophore has been described to identify novel GPR40 modulators. Methods: DISCOtech module from Discovery studio was used for the generation of the Structure and Ligand based pharmacophore models which gave hydrophobic aromatic, ring aromatic and negative ionizable as essential pharmacophoric features. The generated models were validated by screening active and inactive datasets, GH scoring and ROC curve analysis. The best model was exposed as a 3D query to screen the hits from databases like GLASS (GPCR-Ligand Association), GPCR SARfari and Mini-Maybridge. Various filters were applied to retrieve the hit molecules having good drug-like properties. A known protein structure of hGPR40 (pdb: 4PHU) having TAK-875 as ligand complex was used to perform the molecular docking studies; using SYBYL-X 1.2 software. Results and Conclusion: Clustering both the models gave RMSD of 0.89. Therefore, the present approach explored the maximum features by combining both ligand and structure based pharmacophore models. A common structural motif as identified in combiphore for GPR40 modulation consists of the para-substituted phenyl propionic acid scaffold. Therefore, the combiphore approach, whereby maximum structural information (from both ligand and biological protein) is explored, gives maximum insights into the plausible protein-ligand interactions and provides potential lead candidates as exemplified in this study.


2021 ◽  
Vol 13 (1) ◽  
Author(s):  
Xujun Zhang ◽  
Chao Shen ◽  
Xueying Guo ◽  
Zhe Wang ◽  
Gaoqi Weng ◽  
...  

AbstractVirtual screening (VS) based on molecular docking has emerged as one of the mainstream technologies of drug discovery due to its low cost and high efficiency. However, the scoring functions (SFs) implemented in most docking programs are not always accurate enough and how to improve their prediction accuracy is still a big challenge. Here, we propose an integrated platform called ASFP, a web server for the development of customized SFs for structure-based VS. There are three main modules in ASFP: (1) the descriptor generation module that can generate up to 3437 descriptors for the modelling of protein–ligand interactions; (2) the AI-based SF construction module that can establish target-specific SFs based on the pre-generated descriptors through three machine learning (ML) techniques; (3) the online prediction module that provides some well-constructed target-specific SFs for VS and an additional generic SF for binding affinity prediction. Our methodology has been validated on several benchmark datasets. The target-specific SFs can achieve an average ROC AUC of 0.973 towards 32 targets and the generic SF can achieve the Pearson correlation coefficient of 0.81 on the PDBbind version 2016 core set. To sum up, the ASFP server is a powerful tool for structure-based VS.


2017 ◽  
Vol 74 (1) ◽  
pp. 108-112 ◽  
Author(s):  
Urszula Anna Budniak ◽  
Paulina Maria Dominiak

Isoguanine, an analogue of guanine, is of intrinsic interest as a noncanonical nucleobase. The crystal structure of isoguaninium chloride (systematic name: 6-amino-2-oxo-1H,7H-purin-3-ium chloride), C5H6N5O+·Cl−, has been determined by single-crystal X-ray diffraction. Structure analysis was supported by electrostatic interaction energy (E es) calculations based on charge density reconstructed with the UBDB databank. In the structure, two kinds of molecular tapes are observed, one parallel to (010) and the other parallel to (50\overline{4}). The tapes are formed by dimers of isoguaninium cations interacting with chloride anions. E es analysis indicates that cations in one kind of tape are oriented so as to minimize repulsive electrostatic interactions.


2021 ◽  
Vol 35 (08) ◽  
pp. 2130002
Author(s):  
Connor J. Morris ◽  
Dennis Della Corte

Molecular docking and molecular dynamics (MD) are powerful tools used to investigate protein-ligand interactions. Molecular docking programs predict the binding pose and affinity of a protein-ligand complex, while MD can be used to incorporate flexibility into docking calculations and gain further information on the kinetics and stability of the protein-ligand bond. This review covers state-of-the-art methods of using molecular docking and MD to explore protein-ligand interactions, with emphasis on application to drug discovery. We also call for further research on combining common molecular docking and MD methods.


2020 ◽  
Vol 21 (15) ◽  
pp. 5183 ◽  
Author(s):  
Eric D. Boittier ◽  
Yat Yin Tang ◽  
McKenna E. Buckley ◽  
Zachariah P. Schuurs ◽  
Derek J. Richard ◽  
...  

A promising protein target for computational drug development, the human cluster of differentiation 38 (CD38), plays a crucial role in many physiological and pathological processes, primarily through the upstream regulation of factors that control cytoplasmic Ca2+ concentrations. Recently, a small-molecule inhibitor of CD38 was shown to slow down pathways relating to aging and DNA damage. We examined the performance of seven docking programs for their ability to model protein-ligand interactions with CD38. A test set of twelve CD38 crystal structures, containing crystallized biologically relevant substrates, were used to assess pose prediction. The rankings for each program based on the median RMSD between the native and predicted were Vina, AD4 > PLANTS, Gold, Glide, Molegro > rDock. Forty-two compounds with known affinities were docked to assess the accuracy of the programs at affinity/ranking predictions. The rankings based on scoring power were: Vina, PLANTS > Glide, Gold > Molegro >> AutoDock 4 >> rDock. Out of the top four performing programs, Glide had the only scoring function that did not appear to show bias towards overpredicting the affinity of the ligand-based on its size. Factors that affect the reliability of pose prediction and scoring are discussed. General limitations and known biases of scoring functions are examined, aided in part by using molecular fingerprints and Random Forest classifiers. This machine learning approach may be used to systematically diagnose molecular features that are correlated with poor scoring accuracy.


2020 ◽  
Vol 21 (3) ◽  
pp. 893 ◽  
Author(s):  
Nurbubu T. Moldogazieva ◽  
Daria S. Ostroverkhova ◽  
Nikolai N. Kuzmich ◽  
Vladimir V. Kadochnikov ◽  
Alexander A. Terentiev ◽  
...  

Alpha-fetoprotein (AFP) is a major embryo- and tumor-associated protein capable of binding and transporting a variety of hydrophobic ligands, including estrogens. AFP has been shown to inhibit estrogen receptor (ER)-positive tumor growth, which can be attributed to its estrogen-binding ability. Despite AFP having long been investigated, its three-dimensional (3D) structure has not been experimentally resolved and molecular mechanisms underlying AFP–ligand interaction remains obscure. In our study, we constructed a homology-based 3D model of human AFP (HAFP) with the purpose of molecular docking of ERα ligands, three agonists (17β-estradiol, estrone and diethylstilbestrol), and three antagonists (tamoxifen, afimoxifene and endoxifen) into the obtained structure. Based on the ligand-docked scoring functions, we identified three putative estrogen- and antiestrogen-binding sites with different ligand binding affinities. Two high-affinity binding sites were located (i) in a tunnel formed within HAFP subdomains IB and IIA and (ii) on the opposite side of the molecule in a groove originating from a cavity formed between domains I and III, while (iii) the third low-affinity binding site was found at the bottom of the cavity. Here, 100 ns molecular dynamics (MD) simulation allowed us to study their geometries and showed that HAFP–estrogen interactions were caused by van der Waals forces, while both hydrophobic and electrostatic interactions were almost equally involved in HAFP–antiestrogen binding. Molecular mechanics/Generalized Born surface area (MM/GBSA) rescoring method exploited for estimation of binding free energies (ΔGbind) showed that antiestrogens have higher affinities to HAFP as compared to estrogens. We performed in silico point substitutions of amino acid residues to confirm their roles in HAFP–ligand interactions and showed that Thr132, Leu138, His170, Phe172, Ser217, Gln221, His266, His316, Lys453, and Asp478 residues, along with two disulfide bonds (Cys224–Cys270 and Cys269–Cys277), have key roles in both HAFP–estrogen and HAFP–antiestrogen binding. Data obtained in our study contribute to understanding mechanisms underlying protein–ligand interactions and anticancer therapy strategies based on ERα-binding ligands.


Sign in / Sign up

Export Citation Format

Share Document