共查询到20条相似文献，搜索用时 15 毫秒
1.
D. M. Lorber B. K. Shoichet 《Protein science : a publication of the Protein Society》1998,7(4):938-950
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. 相似文献
2.
Docking is a computational technique that places a small molecule (ligand) in the binding site of its macromolecular target (receptor) and estimates its binding affinity. This review addresses methodological developments that have occurred in the docking field in 2009, with a particular focus on the more difficult, and sometimes controversial, aspects of this promising computational discipline. These developments aim to address the main challenges of docking: receptor representation (such aspects as structural waters, side chain protonation, and, most of all, flexibility (from side chain rotation to domain movement)), ligand representation (protonation, tautomerism and stereoisomerism, and the effect of input conformation), as well as accounting for solvation and entropy of binding. This review is strongly focused on docking advances in the context of drug design, specifically in virtual screening and fragment-based drug design. 相似文献
3.
蛋白质分子间相互作用与识别是21世纪生命科学研究的前沿和热点.分子对接方法是研究这一课题有效的计算机模拟手段.通常,蛋白质-蛋白质分子对接包括四个阶段:搜索受体与配体分子间的结合模式,过滤对接结构以排除不合理的结合模式,优化结构,用精细的打分函数评价、排序对接模式并挑选近天然构象.结合国内外研究蛋白质-蛋白质分子对接方法进展和本研究小组的工作,对以上四个环节做了详细的综述.另外,还分析了目前存在的主要问题,并提出对未来工作的展望. 相似文献
4.
Ligand induced fit phenomenon occurring at the ligand binding domain of the liver X receptor beta (LXRbeta) was investigated by means of molecular dynamics. Reliability of a 4-ns trajectory was tested from two distinct LXRbeta crystal complexes 1PQ6B/GW and 1PQ9B/T09 characterized by an open and a closed state of the pocket, respectively. Crossed complexes 1PQ6B/T09 and 1PQ9B/GW were then submitted to the same molecular dynamic conditions, which were able to recover LXRbeta conformations similar to the original crystallography data. Analysis of \"open to closed\" and \"closed to open\" conformational transitions pointed out the dynamic role of critical residues lining the ligand binding pocket involved in the local remodeling upon ligand binding (e.g., Phe271, Phe329, Phe340, Arg319, Glu281). Altogether, the present study indicates that the molecular dynamic protocol is a consistent approach for managing LXRbeta-related induced fit process. This protocol could therefore be used for refining ligand docking solutions of a structure-based design strategy. 相似文献
5.
Cancer cells are resistant to apoptosis and this is one of the most obvious symptoms of cancer in humans. One of the most exciting strategies for treating cancer is to design regulators that increase cell death and stop cell growth. Members of the BCL-2 family of proteins play an important role in the regulation of apoptosis. In this study, an attempt was made to improve the performance of one of the anticancer drugs by designing new analogs of venetoclax (VNT). For this purpose, molecular docking studies were performed to determine the best binding state of VNT and its newly designed derivatives at the protein-binding site to estimate the binding energy. The best analog in terms of free energy was VNT-12 with the lowest energy (−12.15 kcal/mol). Finally, to investigate the inhibitory effect of the compounds on BCL-2 protein, molecular dynamics simulation was used, and by performing the relevant analyses during the simulation, it was observed that the newly designed ligand had better performance in inhibiting BCL-2 protein compared to VNT. 相似文献
6.
A simple method is described to perform docking of subtrates to proteins or probes to receptor molecules by a modification of molecular dynamics simulations. The method consists of a separation of the center-of-mass motion of the substrate from its internal and rotational motions, and a separate coupling to different thermal baths for both types of motion of the substrate and for the motion of the receptor. Thus the temperatures and the time constants of coupling to the baths can be arbitrarily varied for these three types of motion, allowing either a frozen or a flexible receptor and allowing control of search rate without disturbance of internal structure. In addition, an extra repulsive term between substrate and protein was applied to smooth the interaction. The method was applied to a model substrate docking onto a model surface, and to the docking of phosphocholine onto immunoglobulin McPC603, in both cases with a frozen receptor. Using transrational temperatures of the substrate in the range of 1300–1700 K and room temperature for the internal degrees of freedom of the substrate, an efficient nontrapping exploratory search (“helicopter view”) is obtained, which visits the correct binding sites. Low energy conformations can then be further investigated by separate search or by dynamic simulated annealing. In both cases the correct minima were identified. The possibility to work with flexible receptors is discussed. © 1994 Wiley-Liss, Inc. 相似文献
7.
The 1TW7 crystal structure of HIV-1 protease shows the flaps placed wider and more open than what is seen in other examples of the semi-open, apo form. It has been proposed that this might be experimental evidence of allosteric control, because crystal packing creates contacts to the \"elbow region\" of the protease, which may cause deformation of the flaps. Recent dynamics simulations have shown that the conformation seen in 1TW7 relaxes into the typical semi-open conformation in the absence of the crystal contacts, definitively showing that the crystal contacts cause the deformation (Layten et al., J Am Chem Soc 2006;128:13360-13361). However, this does not prove or disprove allosteric modulation at the elbow. In this study, we have conducted additional simulations, supplemented with experimental testing, to further probe the possibility of 1TW7 providing an example of allosteric control of the flap region. We show that the contacts are unstable and do not restrict the conformational sampling of the flaps. The deformation seen in the 1TW7 crystal structure is simply opportunistic crystal packing and not allosteric control. 相似文献
8.
A recent crystal structure of HIV-1 protease (HIVp) was the first to experimentally observe a ligand targeting an open-flap conformation. Researchers studying a symmetric pyrrolidine inhibitor found that two ligands cocrystallized with the protease, forcing an unusual configuration and unique crystallographic contacts. One molecule is centered in the traditional binding site (α pose) and the other binds between the flaps (β pose). The ligands stack against each other in a region termed the \"eye\" site. Ligands bound to the eye site should prevent flap closure, but it is unclear if the pyrrolidine inhibitors or the crystal packing are causing the open state. Molecular dynamics simulations were used to examine the solution-state behavior of three possible binding modes: the ternary complex of HIVp+αβ and the binary complexes, HIVp+α and HIVp+β. We show that HIVp+α is the most stable of the three states. During conformational sampling, α takes an asymmetric binding pose, with one naphthyl ring occupying the eye site and the other reoriented down to occupy positions seen with traditional inhibitors. This finding supports previous studies that reveal a requirement for asymmetric binding at the eye site. In fact, if the α pose is modified to splay both naphthyl rings across the binding site like traditional inhibitors, one ring consistently flips to occupy the eye site. Our simulations reveal that interactions to the eye site encourage a conformationally restrained state, and understanding those contacts may aid the design of ligands to specifically target alternate conformations of the protease. 相似文献
9.
10.
Computational biology methods are now firmly entrenched in the drug discovery process. These methods focus on modeling and
simulations of biological systems to complement and direct conventional experimental approaches. Two important branches of
computational biology include protein homology modeling and the computational biophysics method of molecular dynamics. Protein
modeling methods attempt to accurately predict three-dimensional (3D) structures of uncrystallized proteins for subsequent
structure-based drug design applications. Molecular dynamics methods aim to elucidate the molecular motions of the static
representations of crystallized protein structures. In this review we highlight recent novel methodologies in the field of
homology modeling and molecular dynamics. Selected drug discovery applications using these methods conclude the review. 相似文献
11.
Kazeem O. Sulaiman Temitope U. Kolapo Abdulmujeeb T. Onawole Md. Ataul Islam Rukayat O. Adegoke Suaibu O. Badmus 《Journal of biomolecular structure & dynamics》2019,37(12):3029-3040
Ebola virus (EBOV) is a lethal human pathogen with a risk of global spread of its zoonotic infections, and Ebolavirus Zaire specifically has the highest fatality rate amongst other species. There is a need for continuous effort towards having therapies, as a single licensed treatment to neutralize the EBOV is yet to come into reality. This present study virtually screened the MCULE database containing almost 36 million compounds against the structure of a Zaire Ebola viral protein (VP) 35 and a consensus scoring of both MCULE and CLCDDW docking programs remarked five compounds as potential hits. These compounds, with binding energies ranging from –7.9 to –8.9?kcal/mol, were assessed for predictions of their physicochemical and bioactivity properties, as well as absorption, distribution, metabolism, excretion, and toxicity (ADMET) criteria. The results of the 50?ns molecular dynamics simulations showed the presence of dynamic stability between ligand and protein complexes, and the structures remained significantly unchanged at the ligand-binding site throughout the simulation period. Both docking analysis and molecular dynamics simulation studies suggested strong binding affinity towards the receptor cavity and these selected compounds as potential inhibitors against the Zaire Ebola VP 35. With respect to inhibition constant values, bioavailability radar and other physicochemical properties, compound A (MCULE-1018045960-0-1) appeared to be the most promising hit compound. However, the ligand efficiency and ligand efficiency scale need improvement during optimization, and also validation via in vitro and in vivo studies are necessary to finally make a lead compound in treating Ebola virus diseases.
Communicated by Ramaswamy H. Sarma 相似文献
12.
13.
The aim of this study was to investigate the usefulness of structure-based virtual screening (VS) for focused library design in G protein-coupled receptors (GPCR) projects on the example of 5-HT(2c) agonists. We compared the performance of structure-based VS against two different homology models using FRED for docking and ScreenScore, FlexX, and PMF for rescoring with the results of 12 ligand-based similarity searches using four different query compounds and three different similarity metrics (Daylight, FTree, Phacir). The result of the similarity search showed much variation, from an enrichment factor up to 3.2 to worse than random, whereas the structure-based VS gave a more stable result with a constant enrichment factor around 2. Additionally, actives retrieved by the structure-based approach were more diverse than the actives among the top scorers of the similarity searches. Based on these results, we suggest basing a focused library design for a GPCR project on a combination of a ligand-based similarity search and structure-based docking. 相似文献
14.
Angiotensin II type 1 (AT(1)) receptor belongs to the super-family of G-protein-coupled receptors, and antagonists of the AT(1) receptor are effectively used in the treatment of hypertension. To understand the molecular interactions of these antagonists, such as losartan and telmisartan, with the AT(1) receptor, a homology model of the human AT(1) (hAT(1)) receptor with all connecting loops was constructed from the 2.6 A resolution crystal structure (PDB i.d., 1L9H) of bovine rhodopsin. The initial model generated by MODELLER was subjected to a stepwise ligand-supported model refinement. This protocol involved initial docking of non-peptide AT(1) antagonists in the putative binding site, followed by several rounds of iterative energy minimizations and molecular dynamics simulations. The final model was validated based on its correlation with several structure-activity relationships and site-directed mutagenesis data. The final model was also found to be in agreement with a previously reported AT(1) antagonist pharmacophore model. Docking studies were performed for a series of non-peptide AT(1) receptor antagonists in the active site of the final hAT(1) receptor model. The docking was able to identify key molecular interactions for all the AT(1) antagonists studied. Reasonable correlation was observed between the interaction energy values and the corresponding binding affinities of these ligands, providing further validation for the model. In addition, an extensive unrestrained molecular dynamics simulation showed that the docking-derived bound pose of telmisartan is energetically stable. Knowledge gained from the present studies can be used in structure-based drug design for developing novel ligands for the AT(1) receptor. 相似文献
15.
Knowing the ligand or peptide binding site in proteins is highly important to guide drug discovery, but experimental elucidation of the binding site is difficult. Therefore, various computational approaches have been developed to identify potential binding sites in protein structures. However, protein and ligand flexibility are often neglected in these methods due to efficiency considerations despite the recognition that protein–ligand interactions can be strongly affected by mutual structural adaptations. This is particularly true if the binding site is unknown, as the screening will typically be performed based on an unbound protein structure. Herein we present DynaBiS, a hierarchical sampling algorithm to identify flexible binding sites for a target ligand with explicit consideration of protein and ligand flexibility, inspired by our previously presented flexible docking algorithm DynaDock. DynaBiS applies soft-core potentials between the ligand and the protein, thereby allowing a certain protein–ligand overlap resulting in efficient sampling of conformational adaptation effects. We evaluated DynaBiS and other commonly used binding site identification algorithms against a diverse evaluation set consisting of 26 proteins featuring peptide as well as small ligand binding sites. We show that DynaBiS outperforms the other evaluated methods for the identification of protein binding sites for large and highly flexible ligands such as peptides, both with a holo or apo structure used as input. 相似文献
16.
Abigail L. Emtage Shailesh N. Mistry Peter M. Fischer Barrie Kellam Charles A. Laughton 《Journal of biomolecular structure & dynamics》2017,35(12):2604-2619
G protein-coupled receptors (GPCRs) are proteins of pharmaceutical importance, with over 30% of all drugs in clinical use targeting them. Increasing numbers of X-ray crystal (XRC) structures of GPCRs offer a wealth of data relating to ligand binding. For the β-adrenoceptors (β-ARs), XRC structures are available for human β2- and turkey β1-subtypes, in complexes with a range of ligands. While these structures provide insight into the origins of ligand structure-activity relationships (SARs), questions remain. The ligands in all published complexed XRC structures lack extensive substitution, with no obvious way the ligand-binding site can accommodate β1-AR-selective antagonists with extended side-chains para- to the common aryloxypropanolamine pharmacophore. Using standard computational docking tools with such ligands generally returns poses that fail to explain known SARs. Application of our Active Site Pressurisation modelling method to β-AR XRC structures and homology models, however, reveals a dynamic area in the ligand-binding pocket that, through minor changes in amino acid side chain orientations, opens a fissure between transmembrane helices H4 and H5, exposing intra-membrane space. This fissure, which we term the “keyhole”, is ideally located to accommodate extended moieties present in many high-affinity β1-AR-selective ligands, allowing the rest of the ligand structure to adopt a canonical pose in the orthosteric binding site. We propose the keyhole may be a feature of both β1- and β2-ARs, but that subtle structural differences exist between the two, contributing to subtype-selectivity. This has consequences for the rational design of future generations of subtype-selective ligands for these therapeutically important targets. 相似文献
17.
Hayes JM Skamnaki VT Archontis G Lamprakis C Sarrou J Bischler N Skaltsounis AL Zographos SE Oikonomakos NG 《Proteins》2011,79(3):703-719
With an aim toward glycogenolysis control in Type 2 diabetes, we have investigated via kinetic experiments and computation the potential of indirubin (IC?? > 50 μM), indirubin-3'-oxime (IC?? = 144 nM), KT5720 (K(i) = 18.4 nM) and staurosporine (K(i) = 0.37 nM) as phosphorylase kinase (PhKγtrnc) ATP-binding site inhibitors, with the latter two revealed as potent inhibitors in the low nM range. Because of lack of structural information, we have exploited information from homologous kinase complexes to direct in silico calculations (docking, molecular dynamics, and MMGBSA) to predict the binding characteristics of the four ligands. All inhibitors are predicted to bind in the same active site area as the ATP adenine ring, with binding dominated by hinge region hydrogen bonds to Asp104:O and Met106:O (all four ligands) and also Met106:NH (for the indirubins). The PhKγtrnc-staurosporine complex has the greatest number of receptor-ligand hydrogen bonds, while for the indirubin-3'-oxime and KT5720 complexes there is an important network of interchanging water molecules bridging inhibitor-enzyme contacts. The MM-GBSA results revealed the source of staurosporine's low nM potency to be favorable electrostatic interactions, while KT5720 has strong van der Waals contributions. KT5720 interacts with the greatest number of protein residues either by direct or 1-water bridged hydrogen bond interactions, and the potential for more selective PhK inhibition based on a KT5720 analogue has been established. Including receptor flexibility in Schr?dinger induced-fit docking calculations in most cases correctly predicted the binding modes as compared with the molecular dynamics structures; the algorithm was less effective when there were key structural waters bridging receptor-ligand contacts. 相似文献
18.
Elizabeth Yuriev Jessica Holien Paul A. Ramsland 《Journal of molecular recognition : JMR》2015,28(10):581-604
Molecular docking is a computational method for predicting the placement of ligands in the binding sites of their receptor(s). In this review, we discuss the methodological developments that occurred in the docking field in 2012 and 2013, with a particular focus on the more difficult aspects of this computational discipline. The main challenges and therefore focal points for developments in docking, covered in this review, are receptor flexibility, solvation, scoring, and virtual screening. We specifically deal with such aspects of molecular docking and its applications as selection criteria for constructing receptor ensembles, target dependence of scoring functions, integration of higher‐level theory into scoring, implicit and explicit handling of solvation in the binding process, and comparison and evaluation of docking and scoring methods. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
19.
The adaptation of forcefield-based scoring function to specific receptors remains an important challenge for in-silico drug discovery. Here we compare binding energies of forcefield-based scoring functions with models that are reparameterized on the basis of large-scale quantum calculations of the receptor. We compute binding energies of eleven ligands to the human estrogen receptor subtype alpha (ERalpha) and four ligands to the human retinoic acid receptor of isotype gamma (RARgamma). Using the FlexScreen all-atom receptor-ligand docking approach, we compare docking simulations parameterized by quantum-mechanical calculation of a large protein fragment with purely forcefield-based models. The use of receptor flexibility in the FlexScreen permits the treatment of all ligands in the same receptor model. We find a high correlation between the classical binding energy obtained in the docking simulation and quantum mechanical binding energies and a good correlation with experimental affinities R=0.81 for ERalpha and R=0.95 for RARgamma using the quantum derived scoring functions. A significant part of this improvement is retained, when only the receptor is treated with quantum-based parameters, while the ligands are parameterized with a purely classical model. 相似文献
20.
Dharmendra Kumar Yadav Surendra Kumar Mahesh Kumar Teli Mi-Hyun Kim 《Journal of cellular biochemistry》2020,121(7):3570-3583
In recent years, pharmacophore modeling and molecular docking approaches have been extensively used to characterize the structural requirements and explore the conformational space of a ligand in the binding pocket of the selected target protein. Herein, we report a pharmacophore modeling and molecular docking of 45 compounds comprising of the indole scaffold as vitamin D receptor (VDR) inhibitors. Based on the selected best hypothesis (DRRRR.61), an atom-based three-dimensional quantitative structure-activity relationships model was developed to rationalize the structural requirement of biological activity modulating components. The developed model predicted the binding affinity for the training set and test set with R2(training) = 0.8869 and R2(test) = 0.8139, respectively. Furthermore, molecular docking and dynamics simulation were performed to understand the underpinning of binding interaction and stability of selected VDR inhibitors in the binding pocket. In conclusion, the results presented here, in the form of functional and structural data, agreed well with the proposed pharmacophores and provide further insights into the development of novel VDR inhibitors with better activity. 相似文献