首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 8 毫秒
1.
    
A new approach to predicting the ligand-binding sites of proteins was developed, using protein-ligand docking computation. In this method, many compounds in a random library are docked onto the whole protein surface. We assumed that the true ligand-binding site would exhibit stronger affinity to the compounds in the random library than the other sites, even if the random library did not include the ligand corresponding to the true binding site. We also assumed that the affinity of the true ligand-binding site would be correlated to the docking scores of the compounds in the random library, if the ligand-binding site was correctly predicted. We call this method the molecular-docking binding-site finding (MolSite) method. The MolSite method was applied to 89 known protein-ligand complex structures extracted from the Protein Data Bank, and it predicted the correct binding sites with about 80-99% accuracy, when only the single top-ranked site was adopted. In addition, the average docking score was weakly correlated to the experimental protein-ligand binding free energy, with a correlation coefficient of 0.44.  相似文献   

2.
    
Fatty‐acid binding proteins (FABPs) are abundantly expressed proteins that bind a range of lipophilic molecules. They have been implicated in the import and intracellular distribution of their ligands and have been linked with metabolic and inflammatory responses in the cells in which they are expressed. Despite their high sequence identity, human intestinal FABP (hIFABP) and rat intestinal FABP (rIFABP) bind some ligands with different affinities. In order to address the structural basis of this differential binding, diffraction‐quality crystals have been obtained of hIFABP and rIFABP in complex with the fluorescent fatty‐acid analogue 11‐(dansylamino)undecanoic acid.  相似文献   

3.
    
We report on the synthesis, activity testing, docking, and quantum mechanical scoring of novel imidazo[1,2‐c]pyrimidin‐5(6H)‐one scaffold for cyclin‐dependent kinase 2 (CDK2) inhibition. A series of 26 compounds substituted with aromatic moieties at position 8 has been tested in in vitro enzyme assays and shown to inhibit CDK2. 2D structure‐activity relationships have ascertained that small substituents at position 8 (up to the size of naphtyl or methoxyphenyl) generally lead to single‐digit micromolar IC50 values, whereas bigger substituents (substituted biphenyls) decreased the compounds' activities. The binding modes of the compounds obtained using Glide docking have exhibited up to 2 hinge‐region hydrogen bonds to CDK2 and differed in the orientation of the inhibitor core and the placement of the 8‐substituents. Semiempirical quantum mechanics‐based scoring identified probable favourable binding modes, which will serve for future structure‐based design and synthetic optimization of substituents of the heterocyclic core. In summary, we have identified a novel core for CDK2 inhibition and will explore it further to increase the potencies of the compounds and also monitor selectivities against other protein kinases.  相似文献   

4.
    
The cyclic AMP receptor protein (CRP) from Escherichia coli regulates the expression of a large number of genes. In this work, CRP has been overexpressed, purified and digested by subtilisin and chymotrypsin. The fragments S‐CRP (digested by subtilisin) and CH‐CRP (digested by chymotrypsin) have been purified and crystallized. Crystals of S‐CRP diffracted to 2.0 Å resolution and belonged to space group P21, with unit‐cell parameters a = 59.7, b = 75.1, c = 128.3 Å, β = 91.5°. Crystals of CH‐CRP diffracted to 2.8 Å resolution and belonged to space group P222, with unit‐cell parameters a = 45.8, b = 60.9, c = 205.6 Å.  相似文献   

5.
Holm J  Hansen SI 《Bioscience reports》2002,22(3-4):455-463
Two molecular forms of the folate binding protein were isolated and purified from human milk by a combination of cation exchange- and affinity chromatography. One protein (27 kDa) was a cleavage product of the other 100 kDa protein as evidenced by N-terminal amino acid sequence homology and a reduction in the molecular size of the latter protein to 27 kDa after cleavage of its hydrophobic glycosylphosphatidylinositol tail by phosphatidylinositol-specific phospholipase C. High-affinity binding of [3H]folate was characterized by upward convex Scatchard plots and increasing ligand binding affinity with decreasing concentrations of both proteins. Downward convex Scatchard plots and binding affinities showing no dependence on the protein concentration were, however, observed in highly diluted solutions of both proteins. Radioligand binding was inhibited by folate analogs, and dissociation of radioligand was slow at pH 7.4 but rapid and complete at pH 5.0 and 3.5. Ligand binding quenched the tryptophan fluorescence of the 27 kDa protein suggesting that tryptophan is present at the binding site and/or ligand binding induces a conformation change that affects tryptophan environment in the protein. The 27 kDa protein representing soluble folate binding protein exhibited a greater affinity for ligand binding than the 100 kDa protein which possesses a hydrophobic tail identical to the one that anchors the folate receptor to the cell membrane.  相似文献   

6.
    
In the design of new enzymes and binding proteins, human intuition is often used to modify computationally designed amino acid sequences prior to experimental characterization. The manual sequence changes involve both reversions of amino acid mutations back to the identity present in the parent scaffold and the introduction of residues making additional interactions with the binding partner or backing up first shell interactions. Automation of this manual sequence refinement process would allow more systematic evaluation and considerably reduce the amount of human designer effort involved. Here we introduce a benchmark for evaluating the ability of automated methods to recapitulate the sequence changes made to computer‐generated models by human designers, and use it to assess alternative computational methods. We find the best performance for a greedy one‐position‐at‐a‐time optimization protocol that utilizes metrics (such as shape complementarity) and local refinement methods too computationally expensive for global Monte Carlo (MC) sequence optimization. This protocol should be broadly useful for improving the stability and function of designed binding proteins. Proteins 2014; 82:858–866. © 2013 Wiley Periodicals, Inc.  相似文献   

7.
    
Cavasotto CN  Orry AJ  Abagyan RA 《Proteins》2003,51(3):423-433
G-protein coupled receptors (GPCRs) are the largest family of cell-surface receptors involved in signal transmission. Drugs associated with GPCRs represent more than one fourth of the 100 top-selling drugs and are the targets of more than half of the current therapeutic agents on the market. Our methodology based on the internal coordinate mechanics (ICM) program can accurately identify the ligand-binding pocket in the currently available crystal structures of seven transmembrane (7TM) proteins [bacteriorhodopsin (BR) and bovine rhodopsin (bRho)]. The binding geometry of the ligand can be accurately predicted by ICM flexible docking with and without the loop regions, a useful finding for GPCR docking because the transmembrane regions are easier to model. We also demonstrate that the native ligand can be identified by flexible docking and scoring in 1.5% and 0.2% (for bRho and BR, respectively) of the best scoring compounds from two different types of compound database. The same procedure can be applied to the database of available chemicals to identify specific GPCR binders. Finally, we demonstrate that even if the sidechain positions in the bRho binding pocket are entirely wrong, their correct conformation can be fully restored with high accuracy (0.28 A) through the ICM global optimization with and without the ligand present. These binding site adjustments are critical for flexible docking of new ligands to known structures or for docking to GPCR homology models. The ICM docking method has the potential to be used to \"de-orphanize\" orphan GPCRs (oGPCRs) and to identify antagonists-agonists for GPCRs if an accurate model (experimentally and computationally validated) of the structure has been constructed or when future crystal structures are determined.  相似文献   

8.
    
We report the structural and biochemical characterization of a novel periplasmic ligand‐binding protein, Dret_0059, from Desulfohalobium retbaense DSM 5692, an organism isolated from Lake Retba, in Senegal. The structure of the protein consists of a unique combination of a periplasmic solute binding protein (SBP) domain at the N‐terminal and a tandem PAS‐like sensor domain at the C‐terminal region. SBP domains are found ubiquitously, and their best known function is in solute transport across membranes. PAS‐like sensor domains are commonly found in signal transduction proteins. These domains are widely observed as parts of many protein architectures and complexes but have not been observed previously within the same polypeptide chain. In the structure of Dret_0059, a ketoleucine moiety is bound to the SBP, whereas a cytosine molecule is bound in the distal PAS‐like domain of the tandem PAS‐like domain. Differential scanning flourimetry support the binding of ligands observed in the crystal structure. There is significant interaction between the SBP and tandem PAS‐like domains, and it is possible that the binding of one ligand could have an effect on the binding of the other. We uncovered three other proteins with this structural architecture in the non‐redundant sequence data base, and predict that they too bind the same substrates. The genomic context of this protein did not offer any clues for its function. We did not find any biological process in which the two observed ligands are coupled. The protein Dret_0059 could be involved in either signal transduction or solute transport.  相似文献   

9.
    
Braun P  Goldberg E  Negron C  von Jan M  Xu F  Nanda V  Koder RL  Noy D 《Proteins》2011,79(2):463-476
The cyclic tetrapyrroles, viz. chlorophylls (Chl), their bacterial analogs bacteriochlorophylls, and hemes are ubiquitous cofactors of biological catalysis that are involved in a multitude of reactions. One systematic approach for understanding how Nature achieves functional diversity with only this handful of cofactors is by designing de novo simple and robust protein scaffolds with heme and/or (bacterio)chlorophyll [(B)Chls]-binding sites. This strategy is currently mostly implemented for heme-binding proteins. To gain more insight into the factors that determine heme-/(B)Chl-binding selectivity, we explored the geometric parameters of (B)Chl-binding sites in a nonredundant subset of natural (B)Chl protein structures. Comparing our analysis to the study of a nonredundant database of heme-binding helical histidines by Negron et al. (Proteins 2009;74:400-416), we found a preference for the m-rotamer in (B)Chl-binding helical histidines, in contrast to the preferred t-rotamer in heme-binding helical histidines. This may be used for the design of specific heme- or (B)Chl-binding sites in water-soluble helical bundles, because the rotamer type defines the positioning of the bound cofactor with respect to the helix interface and thus the protein-binding site. Consensus sequences for (B)Chl binding were identified by combining a computational and database-derived approach and shown to be significantly different from the consensus sequences recommended by Negron et al. (Proteins 2009;74:400-416) for heme-binding helical proteins. The insights gained in this work on helix- (B)Chls-binding pockets provide useful guidelines for the construction of reasonable (B)Chl-binding protein templates that can be optimized by computational tools.  相似文献   

10.
    
Bile acids are generated in vivo from cholesterol in the liver, and they undergo an enterohepatic circulation involving the small intestine, liver, and kidney. To understand the molecular mechanism of this transportation, it is essential to gain insight into the three-dimensional (3D) structures of proteins involved in the bile acid recycling in free and complexed form and to compare them with homologous members of this protein family. Here we report the solution structure of the human ileal lipid-binding protein (ILBP) in free form and in complex with cholyltaurine. Both structures are compared with a previously published structure of the porcine ILBP-cholylglycine complex and with related lipid-binding proteins. Protein structures were determined in solution by using two-dimensional (2D)- and 3D-homo and heteronuclear NMR techniques, leading to an almost complete resonance assignment and a significant number of distance constraints for distance geometry and restrained molecular dynamics simulations. The identification of several intermolecular distance constraints unambiguously determines the cholyltaurine-binding site. The bile acid is deeply buried within ILBP with its flexible side-chain situated close to the fatty acid portal as entry region into the inner ILBP core. This binding mode differs significantly from the orientation of cholylglycine in porcine ILBP. A detailed analysis using the GRID/CPCA strategy reveals differences in favorable interactions between protein-binding sites and potential ligands. This characterization will allow for the rational design of potential inhibitors for this relevant system.  相似文献   

11.
    
Recent advances in modeling protein structures at the atomic level have made it possible to tackle \"de novo\" computational protein design. Most procedures are based on combinatorial optimization using a scoring function that estimates the folding free energy of a protein sequence on a given main-chain structure. However, the computation of the conformational entropy in the folded state is generally an intractable problem, and its contribution to the free energy is not properly evaluated. In this article, we propose a new automated protein design methodology that incorporates such conformational entropy based on statistical mechanics principles. We define the free energy of a protein sequence by the corresponding partition function over rotamer states. The free energy is written in variational form in a pairwise approximation and minimized using the Belief Propagation algorithm. In this way, a free energy is associated to each amino acid sequence: we use this insight to rescore the results obtained with a standard minimization method, with the energy as the cost function. Then, we set up a design method that directly uses the free energy as a cost function in combination with a stochastic search in the sequence space. We validate the methods on the design of three superficial sites of a small SH3 domain, and then apply them to the complete redesign of 27 proteins. Our results indicate that accounting for entropic contribution in the score function affects the outcome in a highly nontrivial way, and might improve current computational design techniques based on protein stability.  相似文献   

12.
    
The Lyme disease agent Borrelia burgdorferi, which is transmitted via a tick vector, is dependent on its tick and mammalian hosts for a number of essential nutrients. Like other bacterial diderms, it must transport these biochemicals from the extracellular milieu across two membranes, ultimately to the B. burgdorferi cytoplasm. In the current study, we established that a gene cluster comprising genes bb0215 through bb0218 is cotranscribed and is therefore an operon. Sequence analysis of these proteins suggested that they are the components of an ABC‐type transporter responsible for translocating phosphate anions from the B. burgdorferi periplasm to the cytoplasm. Biophysical experiments established that the putative ligand‐binding protein of this system, BbPstS (BB0215), binds to phosphate in solution. We determined the high‐resolution (1.3 Å) crystal structure of the protein in the absence of phosphate, revealing that the protein's fold is similar to other phosphate‐binding proteins, and residues that are implicated in phosphate binding in other such proteins are conserved in BbPstS. Taken together, the gene products of bb0215‐0218 function as a phosphate transporter for B. burgdorferi.  相似文献   

13.
14.
15.
  总被引:1,自引:1,他引:0  
Recent evidence suggests that the net effect of electrostatics is generally to destabilize protein binding due to large desolvation penalties. A novel method for computing ligand-charge distributions that optimize the tradeoff between ligand desolvation penalty and favorable interactions with a binding site has been applied to a model for barnase. The result is a ligand-charge distribution with a favorable electrostatic contribution to binding due, in part, to ligand point charges whose direct interaction with the binding site is unfavorable, but which make strong intra-molecular interactions that are uncloaked on binding and thus act to lessen the ligand desolvation penalty.  相似文献   

16.
    
Many proteins function by interacting with other small molecules (ligands). Identification of ligand‐binding sites (LBS) in proteins can therefore help to infer their molecular functions. A comprehensive comparison among local structures of LBSs was previously performed, in order to understand their relationships and to classify their structural motifs. However, similar exhaustive comparison among local surfaces of LBSs (patches) has never been performed, due to computational complexity. To enhance our understanding of LBSs, it is worth performing such comparisons among patches and classifying them based on similarities of their surface configurations and electrostatic potentials. In this study, we first developed a rapid method to compare two patches. We then clustered patches corresponding to the same PDB chemical component identifier for a ligand, and selected a representative patch from each cluster. We subsequently exhaustively as compared the representative patches and clustered them using similarity score, PatSim. Finally, the resultant PatSim scores were compared with similarities of atomic structures of the LBSs and those of the ligand‐binding protein sequences and functions. Consequently, we classified the patches into ~2000 well‐characterized clusters. We found that about 63% of these clusters are used in identical protein folds, although about 25% of the clusters are conserved in distantly related proteins and even in proteins with cross‐fold similarity. Furthermore, we showed that patches with higher PatSim score have potential to be involved in similar biological processes.  相似文献   

17.
    
Henzl MT  Tanner JJ  Tan A 《Proteins》2011,79(3):752-764
Birds express two β-parvalbumin isoforms, parvalbumin 3 and avian thymic hormone (ATH). Parvalbumin 3 from chicken (CPV3) is identical to rat β-parvalbumin (β-PV) at 75 of 108 residues. CPV3 displays intermediate Ca(2+) affinity--higher than that of rat β-parvalbumin, but lower than that of ATH. As in rat β-PV, the attenuation of affinity is associated primarily with the CD site (residues 41-70), rather than the EF site (residues 80-108). Structural data for rat α- and β-parvalbumins suggest that divalent ion affinity is correlated with the similarity of the unliganded and Ca(2+)-bound conformations. We herein present a comparison of the solution structures of Ca(2+)-free and Ca(2+)-bound CPV3. Although the structures are generally similar, the conformations of residues 47 to 50 differ markedly in the two protein forms. These residues are located in the C helix, proximal to the CD binding loop. In response to Ca(2+) removal, F47 experiences much greater solvent accessibility. The side-chain of R48 assumes a position between the C and D helices, adjacent to R69. Significantly, I49 adopts an interior position in the unliganded protein that allows association with the side-chain of L50. Concomitantly, the realignment of F66 and F70 facilitates their interaction with I49 and reduces their contact with residues in the N-terminal AB domain. This reorganization of the hydrophobic core, although less profound, is nevertheless reminiscent of that observed in rat β-PV. The results lend further support to the idea that Ca(2+) affinity correlates with the structural similarity of the apo- and bound parvalbumin conformations.  相似文献   

18.
    
Protein-design methodology can now generate models of protein structures and interfaces with computed energies in the range of those of naturally occurring structures. Comparison of the properties of native structures and complexes to isoenergetic design models can provide insight into the properties of the former that reflect selection pressure for factors beyond the energy of the native state. We report here that sidechains in native structures and interfaces are significantly more constrained than designed interfaces and structures with equal computed binding energy or stability, which may reflect selection against potentially deleterious non-native interactions.  相似文献   

19.
    
Allosteric HIV‐1 integrase (IN) inhibitors (ALLINIs) bind at the dimer interface of the IN catalytic core domain (CCD), and potently inhibit HIV‐1 by promoting aberrant, higher‐order IN multimerization. Little is known about the structural organization of the inhibitor‐induced IN multimers and important questions regarding how ALLINIs promote aberrant IN multimerization remain to be answered. On the basis of physical chemistry principles and from our analysis of experimental information, we propose that inhibitor‐induced multimerization is mediated by ALLINIs directly promoting inter‐subunit interactions between the CCD dimer and a C‐terminal domain (CTD) of another IN dimer. Guided by this hypothesis, we have built atomic models of inter‐subunit interfaces in IN multimers by incorporating information from hydrogen‐deuterium exchange (HDX) measurements to drive protein‐protein docking. We have also developed a novel free energy simulation method to estimate the effects of ALLINI binding on the association of the CCD and CTD. Using this structural and thermodynamic modeling approach, we show that multimer inter‐subunit interface models can account for several experimental observations about ALLINI‐induced multimerization, including large differences in the potencies of various ALLINIs, the mechanisms of resistance mutations, and the crucial role of solvent exposed R‐groups in the high potency of certain ALLINIs. Our study predicts that CTD residues Tyr226, Trp235 and Lys266 are involved in the aberrant multimer interfaces. The key finding of the study is that it suggests the possibility of ALLINIs facilitating inter‐subunit interactions between an external CTD and the CCD‐CCD dimer interface.  相似文献   

20.
    
Ruvinsky AM  Kozintsev AV 《Proteins》2006,62(1):202-208
We present two novel methods to predict native protein-ligand binding positions. Both methods identify the native binding position as the most probable position corresponding to a maximum of a probability distribution function (PDF) of possible binding positions in a protein active site. Possible binding positions are the origins of clusters composed, on the basis of root-mean square deviations (RMSD), from the multiple ligand positions determined by a docking algorithm. The difference between the methods lies in the ways the PDF is derived. To validate the suggested methods, we compare the averaged RMSD of the predicted ligand docked positions relative to the experimentally determined positions for a set of 135 PDB protein-ligand complexes. We demonstrate that the suggested methods improve docking accuracy by as much as 21-24% in comparison with a method that simply identifies the binding position as the energy top-scored ligand position.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号