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 199 records · Page 11Linked to original sources

Homology models of the cannabinoid CB1 and CB2 receptors. A docking analysis study.

The 3D models of both CB1 and CB2 human receptors have been established by homology modeling using as template the X-ray structure of bovine Rhodopsin (code pdb: 1F88) a G-protein-coupled receptor (GPCR). A recursive approach comprising sequence alignment and model building was used to build both models, followed by the refinement of non-conserved regions. The cannabinoid system has been studied by means of docking techniques, using the 3D models of both CB1 and CB2 and well known reference inverse agonist/antagonist compounds. An approach based on the flexibility of the structures has been used to model the receptor-ligand complexes. The structural effects of ligand binding were studied and analyzed on the basis of hydrogen bond interactions, and binding energy calculations. Potential interaction sites of the receptor were determined from analysis of the difference accessible surface area (DASA) study of the protein with and without ligand.

Computer Simulation↗

Homology model of human retinoic acid metabolising enzyme cytochrome P450 26A1 (CYP26A1): active site architecture and ligand binding.

Homology models of cytochrome P450 RA1 (CYP26A1) were constructed using three human P450 structures, CYP2C8, CYP2C9 and CYP3A4 as templates for the model building. Using MOE software the lowest energy CYP26A1 model was then assessed for stereochemical quality and side chain environment. Further active site optimisation of the CYP26A1 model built using the CYP3A4 template was performed by molecular dynamics to generate a final CYP26A1 model. The natural substrate, all-trans-retinoic acid (atRA), and inhibitor R 15866, were docked into the model allowing further validation of the active site architecture. Using the docking studies structurally and functionally important residues were identified with subsequent characterisation of secondary structure. Multiple hydrophobic interactions, including the side chains of TRP112, PHE299, PHE222, PHE84, PHE374 and PRO371, are important for binding of atRA and R115866. Additional hydrogen bonding interactions were noted as follows: atRA-- C==O of the atRA carboxylate group and ARG86; R115866--benzothiazole nitrogen and the backbone NH of SER115.

Alcohol Dehydrogenase↗

Homology model building of the HMG-1 box structural domain.

Nucleoproteins belonging to the HMG-1/2 family possess homologous domains approximately 75 amino acids in length. These domains, termed HMG-1 boxes, are highly structured, compact, and mediate the interaction between HMG-1 box-containing proteins and DNA in a variety of biological contexts. Homology model building experiments on HMG-1 box sequences 'threaded' through the 1H-NMR structure of an HMG-1 box from rat indicate that the domain does not have rigid sequence requirements for its formation. Energy calculations indicate that the structure of all HMG-1 box domains is stabilized primarily through hydrophobic interactions. We have found structural relationships in the absence of statistically significant sequence similarity, identifying several candidate proteins which could possibly assume the same three-dimensional conformation as the rat HMG-1 box motif. The threading technique provides a method by which significant structural similarities in a diverse protein family can be efficiently detected, and the 'structural alignment' derived by this method provides a rational basis through which phylogenetic relationships and the precise sites of interaction between HMG-1 box proteins and DNA can be deduced.

Amino Acid Sequence↗

Homology modelling and docking studies on Varicella Zoster Virus Thymidine kinase.

Thymidine kinase (TK) is the key enzyme in antiviral and suicide gene therapies. While herpes simplex virus type 1 thymidine kinase has been widely studied and crystallised less is known on Varicella Zoster Virus thymidine kinase (VZV TK) and its three-dimensional structure. In this paper we report the model of the three-dimensional structure of VZV TK resulting from a homology modelling study. Subsequent docking studies of the natural substrate deoxythymidine (dT) and known antiviral drugs were performed and shaded new light on the binding characteristics of the enzyme.

Acyclovir↗

Accuracy of a structural homology model for a class II histocompatibility protein, HLA-DR1: comparison to the crystal structure.

Structural homology modeling is used to test the accuracy by which a Class I major histocompatibility complex (MHC) could be used to model a Class II MHC. The crystal structure of HLA-aw68 served as a reference molecule to model HLA-DR1. The resulting model was compared to the recently released crystal structure by Brown et al. (Nature, Vol. 364, p. 33-39 (1993)). The overall tertiary structure motif (two alpha-helices and a beta-sheet forming a peptide binding cleft) was maintained. However, significant deviations in the secondary structure elements were found between the model and the DR1 crystal structure. These deviations were consistent with the differences between Class I and Class II crystal structures. In regions where the model and DR1 crystals structures are most similar, side chain orientations are also similar. Specific peptide-MHC interactions are discussed and compared with the crystal structure results.

Amino Acid Sequence↗

Homology modelling by distance geometry.

BACKGROUND: Unknown protein structures can be predicted from known structures (the scaffolds) with sequences sufficiently homologous to that of the target, based on the observation that similar sequences usually adopt the same fold. When structural equivalences between residues in the scaffold and target proteins are expressed in terms of conserved interatomic distances, the resulting 'distance geometry' representation provides an elegant mechanism for simultaneous restraint satisfaction and bias-free conformation space exploration. RESULTS: We present a homology modelling algorithm based on distance geometry that relies on the gradual projection of simple model chain coordinates into Euclidean spaces with decreasing dimensionality. The similarity between the unknown target structure and the scaffold proteins with known structures was described by mapping secondary structure assignments and specific distance restraints between C alpha atoms onto the model through a multiple alignment. This information was complemented by additional restraints derived from stereochemical considerations and other general aspects of protein structure such as hydrophobic core formation or the absence of tangled mainchains. CONCLUSIONS: The method was capable of quickly locating the correct fold even from an alignment with modest average conservation indicating that it could serve as a fast tool for obtaining correct low-resolution starting conformations for detailed refinement.

Amino Acid Sequence↗

Alignment algorithm for homology modeling and threading.

A DNA/protein sequence comparison is a popular computational tool for molecular biologists. Finding a good alignment implies an evolutionary and/or functional relationship between proteins or genomic loci. Sequential similarity between two proteins indicates their structural resemblance, providing a practical approach for structural modeling, when structure of one of these proteins is known. The first step in the homology modeling is a construction of an accurate sequence alignment. The commonly used alignment algorithms do not provide an adequate treatment of the structurally mismatched residues in locally dissimilar regions. We propose a simple modification of the existing alignment algorithm which treats these regions properly and demonstrate how this modification improves sequence alignments in real proteins.

Algorithms↗

Using multiple structure alignments, fast model building, and energetic analysis in fold recognition and homology modeling.

We participated in the fold recognition and homology sections of CASP5 using primarily in-house software. The central feature of our structure prediction strategy involved the ability to generate good sequence-to-structure alignments and to quickly transform them into models that could be evaluated both with energy-based methods and manually. The in-house tools we used include: a) HMAP (Hybrid Multidimensional Alignment Profile)-a profile-to-profile alignment method that is derived from sequence-enhanced multiple structure alignments in core regions, and sequence motifs in non-structurally conserved regions. b) NEST-a fast model building program that applies an "artificial evolution" algorithm to construct a model from a given template and alignment. c) GRASP2-a new structure and alignment visualization program incorporating multiple structure superposition and domain database scanning modules. These methods were combined with model evaluation based on all atom and simplified physical-chemical energy functions. All of these methods were under development during CASP5 and consequently a great deal of manual analysis was carried out at each stage of the prediction process. This interactive model building procedure has several advantages and suggests important ways in which our and other methods can be improved, examples of which are provided.

Algorithms↗

Homology models of the Yersinia pseudotuberculosis and Yersinia pestis general porins and comparative analysis of their functional and antigenic regions.

The amino acid sequences of the Yersinia pseudotuberculosis porin (YPS) and Y. pestis porin (YPT) have recently deduced but their three-dimensional structures were not known. These sequences were analyzed using the servers 3D-PSSM and PredPort. The YPS and YPT porins were shown to have a high degree of identity (above 50%) in primary and secondary structures. The three-dimensional models of the Yersinia pseudotuberculosis porin (YPS) and Y. pestis porin (YPT) were obtained using the homology modeling approach, SWISS-MODEL Protein Modeling Server and 3-D structure of PhoE porin from E. coli as template. The superposition of the Calpha-atoms of the monomers of the Yersinia porins and PhoE porin gave a root mean square deviations of 0.47 A and 0.43 A for YPS and YPT respectively. Yersinia porins were found to be very similar in their three-dimensional structure to other non-specific enterobacterial porins, having the same features of overall fold and disposition of loop L3. The intrinsic structures of the monomer pores of YPS and YPT were investigated and their conductances were predicted with the program HOLE. The good correspondence between the theoretical and experimental magnitudes of YPS conductance was found. The Yersinia porins were determined to be unusual in containing the substitution, Glu replaced by Val, in a highly conserved pentapeptide (Pro-Glu-Phe-Gly-Gly-Asp), located in the loop L3 tip that disturbs the functionally important cluster of the acidic amino acids in the constriction site. Comparative analysis of structural organization of YPS and E. coli OmpF porin in the regions involved in subunit association and pore lumen was performed. The YPS porin functional properties were predicted. The differences between these porins in polar interactions playing a significant role in stabilization of the porin trimers were found and discussed in term of the variations in trimer stability. The Yersinia porins were shown to have the highest degree of the structural similarity. The differences between the porins were observed in their external loops. Their loops L6 and loops L8 showed 71.4 and 52.9% of sequence identity, respectively. The arrangement of charged residues clustered in the channel external vestibule of these porins was found to be also different suggesting the possible differences in their functional properties. The surface exposed regions of Yersinia porins involved in their potential sequential antigenic determinants were compared. The structural basis of their cross reactivity and antigenic differences is discussed.

Amino Acid Sequence↗

A three-dimensional structure of Plasmodium falciparum serine hydroxymethyltransferase in complex with glycine and 5-formyl-tetrahydrofolate. Homology modeling and molecular dynamics.

Cytosolic Plasmodium falciparum serine hydroxymethyltransferase (pfSHMT) is a potential target for antimalarial chemotherapy. Contrasting with the other enzymes involved in the parasite folate cycle, little information is available about this enzyme, and its crystallographic structure is unknown yet. In this paper, we propose a theoretical low-resolution 3D model for pfSHMT in complex with glycine and 5-formyl tetrahydrofolate (5-FTHF) based on homology modeling by multiple alignment followed by intensive optimization, validation and dynamics simulations in water. Comparison between the active sites of our model and that of crystallographic Human SHMT (hSHMT) revealed key differences that could be useful for the design of new selective inhibitors of pfSHMT.

Amino Acid Sequence↗

Use of homology modeling in conjunction with site-directed mutagenesis for analysis of structure-function relationships of mammalian cytochromes P450.

In recent years, homology modeling has become an important tool to study cytochrome P450 function, especially in conjunction with experimental approaches such as site-directed mutagenesis. Molecular models of mammalian P450s can be constructed based on crystal structures of four bacterial enzymes, P450cam, P450 BM-3, P450terp and P450eryF, using molecular replacement or consensus methods. In a model built by molecular replacement, the coordinates are copied from those of a given template protein, while consensus methods utilize more then one protein as a template and are based on distance geometry calculations. The models can be used to identify or confirm key residues, evaluate enzyme-substrate interactions and explain changes in protein stability and/or regio- and stereospecificity of substrate oxidation upon residue substitution by site-directed mutagenesis. P450 models have also been utilized to analyze binding of inhibitors or activators, as well as alterations in inhibition and activation due to residue replacement.

Animals↗

Homology models of two isozymes of manganese peroxidase: prediction of a Mn(II) binding site.

The three-dimensional structures of two isozymes of manganese peroxidase (MnP) have been predicted from homology modeling using lignin peroxidase as a template. Although highly homologous, MnP differs from LiP by the requirement of Mn(II) as an intermediate in its oxidation of substrates. The Mn(II) site is absent in LiP and unique to the MnP family of peroxidases. The model structures were used to identify the unique Mn(II) binding sites, to determine to what extent they were conserved in the two isozymes, and to provide insight into why this site is absent in LiP. For each isozyme of MnP, three candidate Mn(II) binding sites were identified. Energy optimizations of the three possible Mn(II) enzyme complexes allowed the selection of the most favorable Mn(II) binding site as one with the most anionic oxygen moieties best configured to act as ligands for the Mn(II). At the preferred site, the Mn(II) is coordinated to the carboxyl oxygens of Glu-35, Glu-39, and Asp-179, and a propionate group of the heme. The predicted Mn(II) binding site is conserved in both isozymes. Comparison between the residues at this site in MnP and the corresponding residues in LiP shows that two of the three anionic residues in MnP are replaced by neutral residues in LiP, explaining why LiP does not bind Mn(II).

Amino Acid Sequence↗

Structure of the human glycogen-associated protein phosphatase 1 regulatory subunit hGM: homology modeling revealed an (alpha/beta)8-barrel-like fold in the multidomain protein.

Protein phosphatase 1 (PP1) is widely distributed among tissues and species and acts as a regulator of many important cellular processes. By targeting the catalytic part of PP1 (PP1C) toward particular loci and substrates, regulatory subunits constitute key elements conferring specificity to the holoenzyme. Here, we report the identification of an (alpha/beta)8-barrel-like structure within the N-ter stretch of the human PP1 regulatory subunit hGM, which is part of the family of diverse proteins associated with glycogen metabolism. Protein homology modeling gave rise to a three-dimensional (3D) model for the 381 N-ter residue stretch of hGM, based on sequence similarity with Streptomyces olivochromogenes xylose isomerase, identified by using FASTA. The alignment was subsequently extended by using hydrophobic cluster analysis. The homology-derived model includes the putative glycogen binding area located within the 142-230 domain of hGM as well as a structural characterization of the PP1C interacting domain (segment 51-67). Refinement of the latter by molecular dynamics afforded a topology that is in agreement with previous X-ray studies (Egloff et al., 1997). Finite difference Poisson-Boltzmann calculations performed on the interacting domains of PP1C and hGM confirm the complementarity of the local electrostatic potentials of the two partners. This work highlights the presence of a conserved fold among distant species (mammalian, Caenorhabditis elegans, yeast) and, thus, emphasizes the involvement of PP1 in crucial basic cellular functions.

Amino Acid Motifs↗

Homology modelling of human CYP2E1 based on the CYP2C5 crystal structure: investigation of enzyme-substrate and enzyme-inhibitor interactions.

The construction of a homology model of human cytochrome P450 2E1 (CYP2E1) is reported, based on the CYP2C5 crystallographic template. A relatively high degree of primary sequence homology (identity=59%), as expected for proteins of the same CYP family, ensured a straightforward generation of the 3-dimensional model due to relatively few deletions and insertions of amino acid residues with respect to the CYP2C5 crystal structure. Probing the CYP2E1 model with typical substrates of the enzyme showed a good agreement with experimental information in the form of positions of metabolism for substrates, and with site-directed mutagenesis data on certain residues. Furthermore, quantitative relationships between substrate binding affinity and various structural parameters associated with the substrate molecules facilitated the formulation of a procedure for estimating relative binding energy and, consequently, K(m) or K(D) values towards the CYP2E1 enzyme. This method has been based on a consideration of the active site interactions between substrates and key amino acid residues lining the haem pocket, together with compound lipophilicity data from partition coefficients.

Amino Acid Sequence↗

A model for human cytochrome P450 2D6 based on homology modeling and NMR studies of substrate binding.

The cytochrome P450 responsible for the debrisoquine/sparteine polymorphism (P450 2D6) has been produced in large quantities by expression of a modified cDNA in baculovirus. A polyhistidine extension was incorporated at the C-terminus of the expressed protein, which, after purification of the protein on a nickel-agarose column, could be removed proteolytically by treatment with thrombin. Purified yields of P450 2D6 were 2.4 mg from 700 mL of cell culture. The protein had a greater than 90% heme content and was fully active, having no residual absorbance at 420 nm in the reduced CO complex. The quantities produced allowed direct study of the interaction of the substrate codeine with the enzyme by paramagnetic relaxation effects on the NMR spectrum of the substrate. Distances between the heme iron atom and substrate protons were calculated from these experiments, and the orientation of the substrate in the binding pocket was determined. This showed that codeine was bound with the methoxy group of the molecule closest to the heme iron (iron-methyl proton distance of 3.1 +/- 0.1 A), consistent with the observed O-demethylation to morphine. A model of the complex Of P450 2D6 with codeine was built from a multiple sequence and structure alignment of the known crystal structures for P450s, incorporating the experimental constraints derived from the NMR studies. This showed that the overall fold Of P450 2D6 is more similar to that of P450 BM3 than to either P450 cam or P450 terp. Codeine binds to P450 2D6 so that the methoxy group is directly above the A ring of the heme, while the basic nitrogen interacts with the carboxylate of aspartate 301.

Amino Acid Sequence↗

Homology modeling and docking mechanism of the mercaptosuccinate and methotrexate to P. falciparum 1-Cys peroxiredoxin: a preliminary molecular study.

A three-dimensional (3-D) model of 1-Cys peroxiredoxin from P. falciparum (Pf-Prx) has been constructed by homology modeling. The model building was based on a structural alignment with the human 1-Cys peroxiredoxin ray structure. First, mercaptosuccinate was docked by Molecular and Quantum Mechanics at the active site in both isozymes, evidencing the role of different residues in the ligand-protein interaction. Stable conformation of the inhibitor in the active site was obtained from the conformational analysis by molecular dynamics. Next, The complex was reoptimized by semiempirical molecular orbital AM1 method. Conformational and frontier orbitals analyses of the ligand-protein complex were carried out in an attempt to obtain structural insight into the inhibition mechanism. Finally, the docking study of the methotrexate (MTX), an anticancer drug also used as an antimalarial inhibitor, into the modes binding site was performed. From the resulting stable complex structure, it was found that the glutamate ring of MTX fits the active site with high complementarity. The glutamate ring formed two hydrogen bonds to the imidazol group of His41 and the amino groups of Arg129. The side-chain of glutamate was in close proximity to the sulfur atom of the catalytic residue, Cys47. This binding mode suggests a possible inhibition mechanism, whereby the cysteine residue is covered with the glutamate ring of the MTX inhibitor, forming an enzyme-ligand adduct. In addition, the higher interaction energies and the molecular orbitals localization between the Pf-Prx active site and the inhibitors alluded to the probable binding sites of the ligand nucleophilic ring.

Amino Acid Sequence↗

Design of novel N-(2,4-dioxo-1,2,3,4-tetrahydro-thieno[3,2-d]pyrimidin-7-yl)-guanidines as thymidine phosphorylase inhibitors, and flexible docking to a homology model.

A novel class of thymidine phosphorylase (TP) inhibitors has been designed based on analogy to the enzyme substrate as well as known inhibitors. Flexible docking studies, using a homology model of human TP, of the designed N-(2,4-dioxo-1,2,3,4-tetrahydro-thieno[3,2-d]pyrimidin-7-yl)-guanidines as well as their synthetic precursors provide insight into the observed experimental trends in binding affinity.

Binding Sites↗

Structural study of homeodomain protein-DNA complexes using a homology modeling approach.

The homeodomain is a conserved protein motif that binds to DNA and plays a central role in gene regulation. We use homeodomain as a model system to study the specific interactions between protein and DNA in a complex. Following the fundamental concept of homology modeling, we have developed an algorithm for predicting structures of both protein and DNA using the known structure of a similar complex as the template. The accuracies of the algorithm in predicting the complex structures are evaluated when two of the homeodomain protein-DNA complexes with known structures (antennapedia and MATalpha2) are selected as test systems. This algorithm allows structural studies of homeodomain binds to DNA with different sequences.

Amino Acid Sequence↗