Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Homology modelling”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7Linked to original sources
PubMed1998

Homology model of the quinohaemoprotein alcohol dehydrogenase from Comamonas testosteroni.

A molecular model of QH-ADH, the quinohaemoprotein alcohol dehydrogenase from Comamonas testosteroni, has been built by homology modelling. Sequence similarity of N-terminal residues 1-570 with the alpha-subunit of quinoprotein methanol dehydrogenases (MDHs) from Methylophilus methylotrophus W3A1 and Methylobacterium extorquens provided a basis for the design of the PQQ-binding domain of QH-ADH. Minimal sequence similarity with cytochrome c551 from Ectothiorhodospira halophila and cytochrome c5 from Azotobacter vinelandii has been used to model the C-terminal haem c-binding domain, residues 571-677, absent in MDHs. Distance constraints inferred from 19F-NMR relaxation studies of trifluoromethylphenylhydrazine-derivatized PQQ bound to QH-ADH apoenzyme as well as theoretical relations for optimal electron transfer have been employed to position the haem- and PQQ-binding domains relative to each other. The homology model obtained shows overall topological similarity with the crystal structure of cd1-nitrite reductase from Thiosphera pantotropha. The proposed model accounts for the following: (i) the site that is sensitive to in vivo proteolytic attack; (ii) the substrate specificity in comparison with MDHs; (iii) changes of the spectral properties of the haem c upon reconstitution of apo-enzyme with PQQ; (iv) electronic interaction between haem and PQQ; and (v) enantioselectivity in the conversion of a chiral sec alcohol.

Alcohol Dehydrogenase↗
PubMed2006

cDNA cloning, expression and homology modeling of a luciferase from the firefly Lampyroidea maculata.

The cDNA of a firefly luciferase from lantern mRNA of Lampyroidea maculata has been cloned, sequenced and functionally expressed. The cDNA has an open reading frame of 1647 bp and codes for a 548-residue-long polypeptide. Noteworthy, sequence comparison as well as homology modeling showed the highest degree of similarity with H. unmunsana and L. mingrelica luciferases, suggesting a close phylogenetic relationship despite the geographical distance separation. The deduced amino acid sequence of the luciferase gene of firefly L. maculata showed 93% identity to H. unmunsana. Superposition of the three-dimensional model of L. maculata luciferase (generated by homology modeling) and three dimensional structure of Photinus pyralis luciferase revealed that the spatial arrangements of Luciferin and ATP-binding residues are very similar. Putative signature of AMPbinding domain among the various firefly species and Lampyroidea maculata was compared and a striking similarity was found. Different motifs and sites have been identified in Lampyroidea maculata by sequence analysis. Expression and purification of luciferase from Lampyroidea maculata was carried out using Ni-NTA Sepharose. Bioluminescence emission spectrum was similar to Photinus pyralis luciferase.

Amino Acid Sequence↗
PubMed2001

Modularity and homology: modelling of the titin type I modules and their interfaces.

Titin is a giant muscle protein with a highly modular architecture consisting of multiple repeats of two sequence motifs, named type I and type II. Type I motifs are homologous to members of the fibronectin type 3 (Fn3) superfamily, one of the motifs most widespread in modular proteins. Fn3 domains are thought to mediate protein-protein interactions and to act as spacers. In titin, Fn3 modules are present in two different super-repeated patterns, likely to be involved in sarcomere assembly through interactions with A-band proteins. Here, we discuss results from homology modelling the whole family of Fn3 domains in titin. Homology modelling is a powerful tool that will play an increasingly important role in the post-genomic era. It is particularly useful for extending experimental structure determinations of parts of multidomain proteins that contain multiple copies of the same motif. The 3D structures of a representative titin type I domain and of other extracellular Fn3 modules were used as a template to model the structures of the 132 copies in titin. The resulting models suggest residues that contribute to the fold stability and allow us to distinguish these from residues likely to have functional importance. In particular, analysis of the models and mapping of the consensus sequence onto the 3D structure suggest putative surfaces of interaction with other proteins. From the structures of isolated modules and the pattern of conservation in the multiple alignment of the whole titin Ig and Fn3 families, it is possible to address the question of how tandem modules are assembled. Our predictions can be validated experimentally.

Amino Acid Motifs↗
PubMed2001

Competitive CYP2C9 inhibitors: enzyme inhibition studies, protein homology modeling, and three-dimensional quantitative structure-activity relationship analysis.

This study describes the generation of a three-dimensional quantitative structure activity relationship (3D-QSAR) model for 29 structurally diverse, competitive CYP2C9 inhibitors defined experimentally from an initial data set of 73 compounds. In parallel, a homology model for CYP2C9 using the rabbit CYP2C5 coordinates was built. For molecules with a known interaction mode with CYP2C9, this homology model, in combination with the docking program GOLD, was used to select conformers to use in the 3D-QSAR analysis. The remaining molecules were docked, and the GRID interaction energies for all conformers proposed by GOLD were calculated. This was followed by a principal component analysis (PCA) of the GRID energies for all conformers of all compounds. Based on the similarity in the PCA plot to the inhibitors with a known interaction mode, the conformer to be used in the 3D-QSAR analysis was selected. The compounds were randomly divided into two groups, the training data set (n = 21) to build the model and the external validation set (n = 8). The PLS (partial least-squares) analysis of the interaction energies against the K(i) values generated a model with r(2) = 0.947 and a cross-validation of q(2) = 0.730. The model was able to predict the entire external data set within 0.5 log units of the experimental K(i) values. The amino acids in the active site showed complementary features to the grid interaction energies in the 3D-QSAR model and were also in agreement with mutagenesis studies.

Animals↗
PubMed2005

Exploring blocker binding to a homology model of the open hERG K+ channel using docking and molecular dynamics methods.

Binding of blockers to the human voltage-gated hERG potassium channel is studied using a combination of homology modelling, automated docking calculations and molecular dynamics simulations, where binding affinities are evaluated using the linear interaction energy method. A homology model was constructed based on the available crystal structure of the bacterial KvAP channel and the affinities of a series of sertindole analogues predicted using this model. The calculations reproduce the relative binding affinities of these compounds very well and indicate that both polar interactions near the intracellular opening of the selectivity filter as well as hydrophobic complementarity in the region around F656 are important for blocker binding. These results are consistent with recent alanine scanning mutation experiments on the blocking of the hERG channel by other compounds.

Amino Acid Sequence↗
PubMed1998

New insulin-like proteins with atypical disulfide bond pattern characterized in Caenorhabditis elegans by comparative sequence analysis and homology modeling.

We have identified three new families of insulin homologs in Caenorhabditis elegans. In two of these families, concerted mutations suggest that an additional disulfide bond links B and A domains, and that the A-domain internal disulfide bond is substituted by a hydrophobic interaction. Homology modeling remarkably confirms these predictions and shows that despite this atypical disulfide bond pattern and the absence of C-like peptide, all these proteins may adopt the same fold as the insulin. Interestingly, whereas we identified 10 insulin-like peptides, only one insulin-like-receptor (daf-2) has been found. We propose that these insulin-related peptides may correspond to different activators or inhibitors of the daf-2 insulin-regulating pathway.

Amino Acid Sequence↗
PubMed2002

Homology modeling reveals the structural background of the striking difference in thermal stability between two related [NiFe]hydrogenases.

Hydrogenases are redox metalloenzymes in bacteria that catalyze the uptake or production of molecular hydrogen. Two homologous nickel-iron hydrogenases, HupSL and HydSL from the photosynthetic purple sulfur bacterium Thiocapsa roseopersicina, differ substantially in their thermal stabilities despite the high sequence similarity between them. The optimum temperature of HydSL activity is estimated to be at least 50 degrees C higher than that of HupSL. In this work, homology models of both proteins were constructed and analyzed for a number of structural properties. The comparison of the models reveals that the higher stability of HydSL can be attributed to increased inter-subunit electrostatic interactions: the homology models reliably predict that HydSL contains at least five more inter-subunit ion pairs than HupSL. The subunit interface of HydSL is more polar than that of HupSL, and it contains a few extra inter-subunit hydrogen bonds. A more optimized cavity system and amino acid replacements resulting in increased conformational rigidity may also contribute to the higher stability of HydSL. The results are in accord with the general observation that with increasing temperature, the role of electrostatic interactions in protein stability increases. Electronic supplementary material to this paper can be obtained by using the Springer Link server located at http://dx.doi.org/10.1007/s00894-001-0071-8.

Amino Acids↗
PubMedDate not supplied

Interactions of aryloxyphenoxypropionic acids with sensitive and resistant acetyl-coenzyme a carboxylase by homology modeling and molecular dynamic simulations.

Acetyl-coenzyme A carboxylase (ACCase) has been identified as one of the most important targets of herbicides. In the present study, we constructed homology models of the carboxyl-transferase (CT) domain of ACCase from sensitive and resistant foxtail and used these models as templates to study the molecular mechanism of herbicide resistance and stereochemistry-activity relationships of aryloxyphenoxypropionates (APPs). In the homology modeling structures, the dimer of the CT domain was formed by the side-to-side arrangement of the two monomers, in such a way that the N domain of one molecule is placed next to the C domain of the other. The dimeric association of sensitive foxtail CT was found to differ from that of resistant foxtail CT, and the spatial orientation of two key residues, Leu-695 and Ile-695, in these dimers also differed. The mutation of Ile to Leu may perturb the conformation of the dimeric interface, which may account for the molecular mechanism of herbicide resistance. Further docking analysis indicated that the binding model of high-active compounds is similar to that in the crystal structure of the enzyme-ligand complex. The different spatial orientations of ester groups of the isomers of APPs may explain the stereochemistry-activity relationship. Ser-698 formed a H-bonding interaction with all of the docked ligands, while Tyr-728 formed a pi-pi stacking interaction with some of the APPs. These findings may enhance our understanding of the molecular mechanism of herbicide resistance and stereochemistry-activity relationships, which may provide a new starting point for the identification of more potent inhibitors against both sensitive and resistant ACCase.

Acetyl-CoA Carboxylase↗
PubMed2003

Amino acid sequence and homology modeling of obtustatin, a novel non-RGD-containing short disintegrin isolated from the venom of Vipera lebetina obtusa.

Disintegrins represent a group of cysteine-rich peptides occurring in Crotalidae and Viperidae snake venoms, and are potent antagonists of several integrin receptors. A novel disintegrin, obtustatin, was isolated from the venom of the Vipera lebetina obtusa viper, and represents the first potent and selective inhibitor of the binding of integrin alpha(1)beta(1) to collagen IV. The primary structure of obtustatin contains 41 amino acids and is the shortest disintegrin described to date. Obtustatin shares the pattern of cysteines of other short disintegrins. However, in contrast to known short disintegrins, the integrin-binding loop of obtustatin is two residues shorter and does not express the classical RGD sequence. Using synthetic peptides, a KTS motif was identified as the integrin-binding sequence. A three-dimensional model of obtustatin, built by homology-modeling structure calculations using different templates and alignments, strongly indicates that the novel KTS motif may reside at the tip of a flexible loop.

Amino Acid Sequence↗
PubMed2005

Functional coverage of the human genome by existing structures, structural genomics targets, and homology models.

The bias in protein structure and function space resulting from experimental limitations and targeting of particular functional classes of proteins by structural biologists has long been recognized, but never continuously quantified. Using the Enzyme Commission and the Gene Ontology classifications as a reference frame, and integrating structure data from the Protein Data Bank (PDB), target sequences from the structural genomics projects, structure homology derived from the SUPERFAMILY database, and genome annotations from Ensembl and NCBI, we provide a quantified view, both at the domain and whole-protein levels, of the current and projected coverage of protein structure and function space relative to the human genome. Protein structures currently provide at least one domain that covers 37% of the functional classes identified in the genome; whole structure coverage exists for 25% of the genome. If all the structural genomics targets were solved (twice the current number of structures in the PDB), it is estimated that structures of one domain would cover 69% of the functional classes identified and complete structure coverage would be 44%. Homology models from existing experimental structures extend the 37% coverage to 56% of the genome as single domains and 25% to 31% for complete structures. Coverage from homology models is not evenly distributed by protein family, reflecting differing degrees of sequence and structure divergence within families. While these data provide coverage, conversely, they also systematically highlight functional classes of proteins for which structures should be determined. Current key functional families without structure representation are highlighted here; updated information on the "most wanted list" that should be solved is available on a weekly basis from http://function.rcsb.org:8080/pdb/function_distribution/index.html.

Databases, Protein↗
PubMed2003

Homology modeling of rat and human cytochrome P450 2D (CYP2D) isoforms and computational rationalization of experimental ligand-binding specificities.

The ligand-binding characteristics of rat and human CYP2D isoforms, i.e., rat CYP2D1-4 and human CYP2D6, were investigated by measuring IC(50) values of 11 known CYP2D6 ligands using 7-methoxy-4-(aminomethyl)coumarin (MAMC) as substrate. Like CYP2D6, all rat CYP2D isozymes catalyzed the O-demethylation of MAMC with K(m) and V(max) values ranging between 78 and 145 microM and 0.048 and 1.122 min(-1), respectively. To rationalize observed differences in the experimentally determined IC(50) values, homology models of the CYP2D isoforms were constructed. A homology model of CYP2D6 was generated on the basis of crystallized rabbit CYP2C5 and was validated on its ability to reproduce binding orientations corresponding to metabolic profiles of the substrates and to remain stable during unrestrained molecular dynamics simulations at 300 K. Twenty-two active site residues, sharing up to 59% sequence identity, were identified in the CYP2D binding pockets and included CYP2D6 residues Phe120, Glu216, and Asp301. Electrostatic potential calculations displayed large differences in the negative charge of the CYP2D active sites, which was consistent with observed differences in absolute IC(50) values. MD studies on the binding mode of sparteine, quinidine, and quinine in CYP2D2 and CYP2D6 furthermore concurred well with experimentally determined IC(50) values and metabolic profiles. The current study thus provides new insights into differences in the active site topology of the investigated CYP2D isoforms.

Alcohol Oxidoreductases↗
PubMed2004

Stochastic algorithm for kinase homology model construction.

A stochastic algorithm for constructing multiple loops in homology modeling of proteins is presented. The algorithm discards variable values in iterations based on a cost function and on statistical analysis of results. Values that remain are used for constructing an ensemble of best solutions. In test cases, the stochastic algorithm retains all the best solutions, compared to an exhaustive scan of the full set of conformations. Individual loops are constructed by adding dipeptide units. Dipeptide conformations are extracted from a database of proteins and their conformations include bond lengths and all angles. Single loops are constructed from both N- and C- terminals to the center, and loop closure is evaluated by a combination of penalties for the peptide closure and Miyazawa-Jernigan (MJ) [1]. residue-residue interactions with the rest of the protein. Large ensembles of each loop are clustered and re-evaluated with a refined [2]. energy term. The reduced, clustered set of each loop is then employed to construct simultaneously all the loops. The algorithm was applied to construct simultaneously six loops in c-Src kinase family proteins, incorporating a total of 37-40 residues. The best RMSD for reconstructing the loops is 1.45 A for Lck (structure 1QPE in the Protein Data Bank) and 2.54 A for human c-Src (structure 1FMK). The multiple loop conformations with lowest energy have higher RMSD values, of 2.06 A and 3.09 A, respectively. The average RMSD values for the first 1000 conformations are 3.00 A and 3.46 A, respectively. Models for the "open" structures of c-Src and of Jak-2 were constructed on the basis of 1QPE. The Jak-2 model is found to be more flexible in the loops region than its c-Src counterpart.

Algorithms↗
PubMed2003

Homology model of human corticosteroid binding globulin: a study of its steroid binding ability and a plausible mechanism of steroid hormone release at the site of inflammation.

Corticosteroid binding globulin (CBG) and thyroxin binding globulin (TBG) both belong to the same SERPIN superfamily of serine-proteinase inhibitors but in the course of evolution CBG has adapted to its new role as a transport agent of insoluble hormones. CBG binds corticosteroids in plasma, delivering them to sites of inflammation to modify the inflammatory response. CBG is an effective drug carrier for genetic manipulation, and hence there is immense biological interest in the location of the hormone binding site. The crystal structure of human CBG (hCBG) has not been determined, but sequence alignment with other SERPINs suggests that it conforms as a whole to the tertiary structure shared by the superfamily. Human CBG shares 52.15% and 55.50% sequence similarity with alpha1-antitrypsin and alpha1-antichymotrypsin, respectively. Multiple sequence alignment among the three sequences shows 73 conserved regions. The molecular structures of alpha1-antitrypsin and alpha1-antichymotrypsin, the archetype of the SERPIN superfamily, obtained by X-ray diffraction methods are used to develop a homology model of hCBG. Energy minimization was applied to the model to refine the structure further. The homology model of hCBG contains 371 residues (His13 to Val383 ). The secondary structure comprises 11 helices, 15 turns and 11 sheets. The putative corticosteroid binding region is found to exist in a pocket between beta-sheets S4, S10, S11 and alpha helix H10. Both cortisol and aldosterone are docked to the elongated hydrophobic ligand binding pocket with the polar residues at the two extremities. A difference accessible surface area (DASA) study revealed that cortisol binds with the native hCBG more tightly than aldosterone. Cleavage at the Val379-Met380 peptide bond causes a deformation of hCBG (also revealed through a DASA study). This deformation could probably trigger the release of the bound hormone. Figure Stereoscopic view of the ribbon diagram of hCBG complexed with cortisol. The bound cortisol is shown in space filling model in blue. Helices and sheets are shown in red and magenta respectively. Turns are shown in yellow.

Adrenal Cortex Hormones↗
PubMed2005

Recognizing and defining true Ras binding domains II: in silico prediction based on homology modelling and energy calculations.

Considering the large number of putative Ras effector proteins, it is highly desirable to develop computational methods to be able to identify true Ras binding molecules. Based on a limited sequence homology among members of the Ras association (RA) and Ras binding (RB) sub-domain families of the ubiquitin super-family, we have built structural homology models of Ras proteins in complex with different RA and RB domains, using the FOLD-X software. A critical step in our approach is to use different templates of Ras complexes, in order to account for the structural variation among the RA and RB domains. The homology models are validated by predicting the effect of mutating hot spot residues in the interface, and residues important for the specificity of interaction with different Ras proteins. The FOLD-X calculated energies of the best-modelled complexes are in good agreement with previously published experimental data and with new data reported here. Based on these results, we can establish energy thresholds above, or below which, we can predict with 96% confidence that a RA/RB domain will or will not interact with Ras. This study shows the importance of in depth structural analysis, high quality force-fields and modelling for correct prediction. Our work opens the possibility of genome-wide prediction for this protein family and for others, where there is enough structural information.

Amino Acid Sequence↗
PubMed1997

Diverse binding site structures revealed in homology models of polyreactive immunoglobulins.

We describe here computer-assisted homology models of the combining site structure of three polyreactive immunoglobulins. Template-based models of Fv (VL-VH) fragments were derived for the surface IgM expressed by the malignant CD5 positive B cells from three patients with chronic lymphocytic leukaemia (CLL). The conserved framework regions were constructed using crystal coordinates taken from highly homologous human variable domain structures (Pot and Hil). Complementarity determining regions (CDRs) were predicted by grafting loops, taken from known immunoglobulin structures, onto the Fv framework models. The CDR templates were chosen, where possible, to be of the same length and of high residue identity or similarity. LCDR1, 2 and 3 as well as HCDR1 and 2 for the Fv were constructed using this strategy. For HCDR3 prediction, a database containing the Cartesian coordinates of 30 of these loops was complied from unliganded antibody X-ray crystallographic structures and an HCDR3 of the same length as that of the B CLL Fv was selected as a template. In one case (Yar), the resulting HCDR3 model gave unfavourable interactions when incorporated into the Fv model. This HCDR3 was therefore modelled using an alternative strategy of construction of the loop stems, using a previously described HCDR3 conformation (Pot), followed by chain closure with a beta-turn. The template models were subjected to positional refinement using energy minimisation and molecular dynamics simulations (X-PLOR). An electrostatic surface description (GRASP) did not reveal a common structural feature within the binding sites of the three polyreactive Fv. Thus, polyreactive immunoglobulins may recognise similar and multiple antigens through a diverse array of binding site structures.

Amino Acid Sequence↗
PubMed2006

Expression and homology modeling of 2-aminobiphenyl-2,3-diol-1,2-dioxygenase from Pseudomonas stutzeri carbazole degradation pathway.

The enzyme 2'-aminobiphenyl-2,3-diol-1,2-dioxygenase (CarB), encoded by two genes (carBa and carBb), is an alpha(2)beta(2) heterotetramer that presents meta-cleavage activity toward the hydroxylated aromatic ring in the carbazole degradation pathway from petroleum-degrader bacteria Pseudomonas spp. The 1,082-base pair polymerase chain reaction product corresponding to carBaBb genes from Pseudomonas stutzeri ATCC 31258 was cloned by site-specific recombination and expressed in high levels in Escherichia coli BL21-SI with a histidine-tag and in native form. The CarB activity toward 2,3-dihydroxybiphenyl was similar for these two constructions. The alpha(2)beta(2)-heterotetrameric 3D model of CarB dioxygenase was proposed by homology modeling using the protocatechuate 4,5-dioxygenase (LigAB) structure as template. Accordingly, His12, His53, and Glu230 coordinate the Fe(II) in the catalytic site at the subunit CarBb. The model also indicates that His182 is the catalytic base responsible for deprotonating one of the hydroxyl group of the substrate by a hydrogen bond. The hydrophobic residues Trp257 and Phe258 in the CarB structure substituted the LigAB amino acid residues Ser269 and Asn270. These data could explain why the CarB was active for 2,3-dihydroxybiphenyl and not for protocatechuate.

Amino Acid Sequence↗
PubMed2005

Inhibitor-based validation of a homology model of the active-site of tripeptidyl peptidase II.

A homology model of the active site region of tripeptidyl peptidase II (TPP II) was constructed based on the crystal structures of four subtilisin-like templates. The resulting model was subsequently validated by judging expectations of the model versus observed activities for a broad set of prepared TPP II inhibitors. The structure-activity relationships observed for the prepared TPP II inhibitors correlated nicely with the structural details of the TPP II active site model, supporting the validity of this model and its usefulness for structure-based drug design and pharmacophore searching experiments.

Amino Acid Sequence↗
PubMed2006

Coverage of whole proteome by structural genomics observed through protein homology modeling database.

We have been developing FAMSBASE, a protein homology-modeling database of whole ORFs predicted from genome sequences. The latest update of FAMSBASE ( http://daisy.nagahama-i-bio.ac.jp/Famsbase/ ), which is based on the protein three-dimensional (3D) structures released by November 2003, contains modeled 3D structures for 368,724 open reading frames (ORFs) derived from genomes of 276 species, namely 17 archaebacterial, 130 eubacterial, 18 eukaryotic and 111 phage genomes. Those 276 genomes are predicted to have 734,193 ORFs in total and the current FAMSBASE contains protein 3D structure of approximately 50% of the ORF products. However, cases that a modeled 3D structure covers the whole part of an ORF product are rare. When portion of an ORF with 3D structure is compared in three kingdoms of life, in archaebacteria and eubacteria, approximately 60% of the ORFs have modeled 3D structures covering almost the entire amino acid sequences, however, the percentage falls to about 30% in eukaryotes. When annual differences in the number of ORFs with modeled 3D structure are calculated, the fraction of modeled 3D structures of soluble protein for archaebacteria is increased by 5%, and that for eubacteria by 7% in the last 3 years. Assuming that this rate would be maintained and that determination of 3D structures for predicted disordered regions is unattainable, whole soluble protein model structures of prokaryotes without the putative disordered regions will be in hand within 15 years. For eukaryotic proteins, they will be in hand within 25 years. The 3D structures we will have at those times are not the 3D structure of the entire proteins encoded in single ORFs, but the 3D structures of separate structural domains. Measuring or predicting spatial arrangements of structural domains in an ORF will then be a coming issue of structural genomics.

Amino Acid Sequence↗