scholarly journals Enrichments of Ensemble Docking Strategy Based on the Bayesian Model

Author(s):  
Yu Lei ◽  
Sheng Guo ◽  
Yi Liu ◽  
Zhili Zuo

Abstract MotivationChallenges remained in structure-based drug discovery which include protein flexibility in binding site. Thus, concerning the flexibility of proteins, docking into an ensemble of rigid conformations (ensemble docking) have been proposed with incorporation into protein flexibility with expects that it could provide higher enrichments than rigid single receptor. Here we have developed the ensemble docking strategy by using Bayesian Model algorithms, and this method is validated by three proteins: BTK, JAK and PARP. The Bayesian Model was used to integrate independent docking runs of an ensemble of rigid crystal structures and MD simulations. ResultsThe structure of MD simulations outperforms the crystal structures in separating inhibitors from decoys in BTK and PARP. Further, the results demonstrated that the ensemble docking strategy has better performance than rigid single conformation.

2019 ◽  
Vol 35 (23) ◽  
pp. 4994-5002 ◽  
Author(s):  
Jiahua He ◽  
Huanyu Tao ◽  
Sheng-You Huang

AbstractMotivationGiven the importance of protein–ribonucleic acid (RNA) interactions in many biological processes, a variety of docking algorithms have been developed to predict the complex structure from individual protein and RNA partners in the past decade. However, due to the impact of molecular flexibility, the performance of current methods has hit a bottleneck in realistic unbound docking. Pushing the limit, we have proposed a protein-ensemble–RNA docking strategy to explicitly consider the protein flexibility in protein–RNA docking through an ensemble of multiple protein structures, which is referred to as MPRDock. Instead of taking conformations from MD simulations or experimental structures, we obtained the multiple structures of a protein by building models from its homologous templates in the Protein Data Bank (PDB).ResultsOur approach can not only avoid the reliability issue of structures from MD simulations but also circumvent the limited number of experimental structures for a target protein in the PDB. Tested on 68 unbound–bound and 18 unbound–unbound protein–RNA complexes, our MPRDock/DITScorePR considerably improved the docking performance and achieved a significantly higher success rate than single-protein rigid docking whether pseudo-unbound templates are included or not. Similar improvements were also observed when combining our ensemble docking strategy with other scoring functions. The present homology model-based ensemble docking approach will have a general application in molecular docking for other interactions.Availability and implementationhttp://huanglab.phys.hust.edu.cn/mprdock/Supplementary informationSupplementary data are available at Bioinformatics online.


2020 ◽  
Vol 9 (5) ◽  
pp. 1473 ◽  
Author(s):  
Umesh Kalathiya ◽  
Monikaben Padariya ◽  
Marcos Mayordomo ◽  
Małgorzata Lisowska ◽  
Judith Nicholson ◽  
...  

An important stage in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) life cycle is the binding of the spike (S) protein to the angiotensin converting enzyme-2 (ACE2) host cell receptor. Therefore, to explore conserved features in spike protein dynamics and to identify potentially novel regions for drugging, we measured spike protein variability derived from 791 viral genomes and studied its properties by molecular dynamics (MD) simulation. The findings indicated that S2 subunit (heptad-repeat 1 (HR1), central helix (CH), and connector domain (CD) domains) showed low variability, low fluctuations in MD, and displayed a trimer cavity. By contrast, the receptor binding domain (RBD) domain, which is typically targeted in drug discovery programs, exhibits more sequence variability and flexibility. Interpretations from MD simulations suggest that the monomer form of spike protein is in constant motion showing transitions between an “up” and “down” state. In addition, the trimer cavity may function as a “bouncing spring” that may facilitate the homotrimer spike protein interactions with the ACE2 receptor. The feasibility of the trimer cavity as a potential drug target was examined by structure based virtual screening. Several hits were identified that have already been validated or suggested to inhibit the SARS-CoV-2 virus in published cell models. In particular, the data suggest an action mechanism for molecules including Chitosan and macrolides such as the mTOR (mammalian target of Rapamycin) pathway inhibitor Rapamycin. These findings identify a novel small molecule binding-site formed by the spike protein oligomer, that might assist in future drug discovery programs aimed at targeting the coronavirus (CoV) family of viruses.


Molecules ◽  
2021 ◽  
Vol 26 (5) ◽  
pp. 1484
Author(s):  
Christopher Reinknecht ◽  
Anthony Riga ◽  
Jasmin Rivera ◽  
David A. Snyder

Proteins are molecular machines requiring flexibility to function. Crystallographic B-factors and Molecular Dynamics (MD) simulations both provide insights into protein flexibility on an atomic scale. Nuclear Magnetic Resonance (NMR) lacks a universally accepted analog of the B-factor. However, a lack of convergence in atomic coordinates in an NMR-based structure calculation also suggests atomic mobility. This paper describes a pattern in the coordinate uncertainties of backbone heavy atoms in NMR-derived structural “ensembles” first noted in the development of FindCore2 (previously called Expanded FindCore: DA Snyder, J Grullon, YJ Huang, R Tejero, GT Montelione, Proteins: Structure, Function, and Bioinformatics 82 (S2), 219–230) and demonstrates that this pattern exists in coordinate variances across MD trajectories but not in crystallographic B-factors. This either suggests that MD trajectories and NMR “ensembles” capture motional behavior of peptide bond units not captured by B-factors or indicates a deficiency common to force fields used in both NMR and MD calculations.


RSC Advances ◽  
2021 ◽  
Vol 11 (31) ◽  
pp. 18938-18944
Author(s):  
Jia-Hong Lei ◽  
Ling-Ling Ma ◽  
Jing-Hong Xian ◽  
Hai Chen ◽  
Jian-Jian Zhou ◽  
...  

Crystal structures of tubulin complexed with ELR510444 and parbendazole facilitate the design of novel colchicine binding site inhibitors.


2015 ◽  
Vol 43 (5) ◽  
pp. 1023-1032 ◽  
Author(s):  
Thomas Stockner ◽  
Anna Mullen ◽  
Fraser MacMillan

ABC transporters are primary active transporters found in all kingdoms of life. Human multidrug resistance transporter ABCB1, or P-glycoprotein, has an extremely broad substrate spectrum and confers resistance against chemotherapy drug treatment in cancer cells. The bacterial ABC transporter MsbA is a lipid A flippase and a homolog to the human ABCB1 transporter, with which it partially shares its substrate spectrum. Crystal structures of MsbA and ABCB1 have been solved in multiple conformations, providing a glimpse into the possible conformational changes the transporter could be going through during the transport cycle. Crystal structures are inherently static, while a dynamic picture of the transporter in motion is needed for a complete understanding of transporter function. Molecular dynamics (MD) simulations and electron paramagnetic resonance (EPR) spectroscopy can provide structural information on ABC transporters, but the strength of these two methods lies in the potential to characterise the dynamic regime of these transporters. Information from the two methods is quite complementary. MD simulations provide an all atom dynamic picture of the time evolution of the molecular system, though with a narrow time window. EPR spectroscopy can probe structural, environmental and dynamic properties of the transporter in several time regimes, but only through the attachment sites of an exogenous spin label. In this review the synergistic effects that can be achieved by combining the two methods are highlighted, and a brief methodological background is also presented.


Biochemistry ◽  
2009 ◽  
Vol 48 (5) ◽  
pp. 1144-1144
Author(s):  
Frédérick Faucher ◽  
Line Cantin ◽  
Van Luu-The ◽  
Fernand Labrie ◽  
Rock Breton

1995 ◽  
Vol 246 (2) ◽  
pp. 344-355 ◽  
Author(s):  
Vincent Mikol ◽  
Götz Baumann ◽  
Thomas H. Keller ◽  
Ute Manning ◽  
Mauro G.M. Zurini

Biomedicines ◽  
2021 ◽  
Vol 9 (3) ◽  
pp. 301
Author(s):  
Ming-Chun Chen ◽  
Yu-Chao Hsiao ◽  
Chun-Chun Chang ◽  
Sheng-Feng Pan ◽  
Chih-Wen Peng ◽  
...  

Congenital nephrogenic diabetes insipidus (CNDI) is a genetic disorder caused by mutations in arginine vasopressin receptor 2 (AVPR2) or aquaporin 2 genes, rendering collecting duct cells insensitive to the peptide hormone arginine vasopressin stimulation for water reabsorption. This study reports a first identified AVPR2 mutation in Taiwan and demonstrates our effort to understand the pathogenesis caused by applying computational structural analysis tools. The CNDI condition of an 8-month-old male patient was confirmed according to symptoms, family history, and DNA sequence analysis. The patient was identified to have a valine 279 deletion–mutation in the AVPR2 gene. Cellular experiments using mutant protein transfected cells revealed that mutated AVPR2 is expressed successfully in cells and localized on cell surfaces. We further analyzed the pathogenesis of the mutation at sub-molecular levels via long-term molecular dynamics (MD) simulations and structural analysis. The MD simulations showed while the structure of the extracellular ligand-binding domain remains unchanged, the mutation alters the direction of dynamic motion of AVPR2 transmembrane helix 6 toward the center of the G-protein binding site, obstructing the binding of G-protein, thus likely disabling downstream signaling. This study demonstrated that the computational approaches can be powerful tools for obtaining valuable information on the pathogenesis induced by mutations in G-protein-coupled receptors. These methods can also be helpful in providing clues on potential therapeutic strategies for CNDI.


Sign in / Sign up

Export Citation Format

Share Document