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 91 records · Page 5Linked to original sources

Homology modeling of opioid receptor-ligand complexes using experimental constraints.

Opioid receptors interact with a variety of ligands, including endogenous peptides, opiates, and thousands of synthetic compounds with different structural scaffolds. In the absence of experimental structures of opioid receptors, theoretical modeling remains an important tool for structure-function analysis. The combination of experimental studies and modeling approaches allows development of realistic models of ligand-receptor complexes helpful for elucidation of the molecular determinants of ligand affinity and selectivity and for understanding mechanisms of functional agonism or antagonism. In this review we provide a brief critical assessment of the status of such theoretical modeling and describe some common problems and their possible solutions. Currently, there are no reliable theoretical methods to generate the models in a completely automatic fashion. Models of higher accuracy can be produced if homology modeling, based on the rhodopsin X-ray template, is supplemented by experimental structural constraints appropriate for the active or inactive receptor conformations, together with receptor-specific and ligand-specific interactions. The experimental constraints can be derived from mutagenesis and cross-linking studies, correlative replacements of ligand and receptor groups, and incorporation of metal binding sites between residues of receptors or receptors and ligands. This review focuses on the analysis of similarity and differences of the refined homology models of mu, delta, and kappa-opioid receptors in active and inactive states, emphasizing the molecular details of interaction of the receptors with some representative peptide and nonpeptide ligands, underlying the multiple modes of binding of small opiates, and the differences in binding modes of agonists and antagonists, and of peptides and alkaloids.

Analgesics, Opioid↗

Homology modeling and molecular dynamics simulation studies of an inward rectifier potassium channel.

A homology model has been generated for the pore-forming domain of Kir6.2, a component of an ATP-sensitive K channel, based on the x-ray structure of the bacterial channel KcsA. Analysis of the lipid-exposed and pore-lining surfaces of the model reveals them to be compatible with the known features of membrane proteins and Kir channels, respectively. The Kir6.2 homology model was used as the starting point for nanosecond-duration molecular dynamics simulations in a solvated phospholipid bilayer. The overall drift from the model structure was comparable to that seen for KcsA in previous similar simulations. Preliminary analysis of the interactions of the Kir6.2 channel model with K(+) ions and water molecules during these simulations suggests that concerted single-file motion of K(+) ions and water through the selectivity filter occurs. This is similar to such motion observed in simulations of KcsA. This suggests that a single-filing mechanism is conserved between different K channel structures and may be robust to changes in simulation details. Comparison of Kir6.2 and KcsA suggests some degree of flexibility in the filter, thus complicating models of ion selectivity based upon a rigid filter.

Amino Acid Sequence↗

High-quality homology models derived from NMR and X-ray structures of E. coli proteins YgdK and Suf E suggest that all members of the YgdK/Suf E protein family are enhancers of cysteine desulfurases.

The structural biology of proteins mediating iron-sulfur (Fe-S) cluster assembly is central for understanding several important biological processes. Here we present the NMR structure of the 16-kDa protein YgdK from Escherichia coli, which shares 35% sequence identity with the E. coli protein SufE. The SufE X-ray crystal structure was solved in parallel with the YdgK NMR structure in the Northeast Structural Genomics (NESG) consortium. Both proteins are (1) key components for Fe-S metabolism, (2) exhibit the same distinct fold, and (3) belong to a family of at least 70 prokaryotic and eukaryotic sequence homologs. Accurate homology models were calculated for the YgdK/SufE family based on YgdK NMR and SufE crystal structure. Both structural templates contributed equally, exemplifying synergy of NMR and X-ray crystallography. SufE acts as an enhancer of the cysteine desulfurase activity of SufS by SufE-SufS complex formation. A homology model of CsdA, a desulfurase encoded in the same operon as YgdK, was modeled using the X-ray structure of SufS as a template. Protein surface and electrostatic complementarities strongly suggest that YgdK and CsdA likewise form a functional two-component desulfurase complex. Moreover, structural features of YgdK and SufS, which can be linked to their interaction with desulfurases, are conserved in all homology models. It thus appears very likely that all members of the YgdK/SufE family act as enhancers of Suf-S-like desulfurases. The present study exemplifies that "refined" selection of two (or more) targets enables high-quality homology modeling of large protein families.

Animals↗

PERMOL: restraint-based protein homology modeling using DYANA or CNS.

SUMMARY: PERMOL is a new restraint-based program for homology modeling of proteins. Restraints are generated from the information contained in structures of homologous template proteins. Employing the restraints generated by PERMOL, three-dimensional structures are obtained using MD programs such as DYANA or CNS. In contrast to other programs PERMOL is mainly based on the use of dihedral angle information which is optimally suited to preserve the local secondary structure. The global arrangement of these elements is then facilitated by a small number of distance restraints. Using PERMOL homology, models of high quality are obtained. A key advantage of the proposed method is its flexibility, which allows the inclusion of data from other sources, such as experimental restraints and the use of modern molecular dynamics programs to calculate structures. AVAILABILITY: The software and a detailed manual are available free of charge (http://www.biologie.uni-regensburg.de/Biophysik/Kalbitzer/permol/permol.html)

Algorithms↗

Homology model and molecular dynamics simulation of carp ovum cystatin.

In this study, a homology model of carp ovum cystatin was constructed based on the crystal structure of chicken egg white cystatin. The results of amino acid sequence alignment indicate that these two proteins exhibit 36.11% of sequence identity. The resultant homology model reveals that carp ovum cystatin shares similar folds as chicken egg white cystatin, particularly in the conserved regions of Q48-V49-G52 and P98-W99 and the locations of two disulfide bonds, C67-C76 and C90-C110. However, the results of 1 ns molecular dynamics simulations show that carp ovum cystatin exhibits less structural integrity than chicken egg white cystatin in explicit water at 300 K. The relatively hydrophilic Met62 of carp ovum cystatin, corresponding to the hydrophobic Leu68 of human cystatin C and Ile66 of chicken egg white cystatin, may destabilize the hydrophobic core and form a dimeric structure more easily through domain swapping. A total of 16 positively charged residues are equally distributed on the surface of carp ovum cystatin, resulting in agglutination with the negatively charged spermatozoa via electrostatic interaction. Thus, carp ovum cystatin is considered to be important in preventing carp eggs from polyspermy.

Agglutination↗

A 3D structure model of integrin alpha 4 beta 1 complex: I. Construction of a homology model of beta 1 and ligand binding analysis.

It is well established that integrin alpha 4 beta 1 binds to the vascular cell adhesion molecule (VCAM) and fibronectin and plays an important role in signal transduction. Blocking the binding of VCAM to alpha 4 beta 1 is thought to be a way of controlling a number of disease processes. To better understand how various inhibitors might block the interaction of VCAM and fibronectin with alpha 4 beta 1, we began constructing a structure model for the integrin alpha 4 beta 1 complex. As the first step, we have built a homology model of the beta 1 subunit based on the I domain of the integrin CD11B subunit. The model, including a bound Mg(2+) ion, was optimized through a specially designed relaxation scheme involving restrained minimization and dynamics steps. The native ligand VCAM and two highly active small molecules (TBC772 and TBC3486) shown to inhibit binding of CS-1 and VCAM to alpha 4 beta 1 were docked into the active site of the refined model. Results from the binding analysis fit well with a pharmacophore model that was independently derived from active analog studies. A critical examination of residues in the binding site and analysis of docked ligands that are both potent and selective led to the proposal of a mechanism for beta 1/beta 7 ligand binding selectivity.

Amino Acid Sequence↗

Protein engineering tests of a homology model of Plasmodium falciparum lactate dehydrogenase.

This paper describes the testing of a homology model of Plasmodium falciparum lactate dehydrogenase (pfLDH) by protein engineering. The model had been validated in structural terms. It suggests explanations of the unusual properties of pfLDH (compared with all other LDHs). These unusual features are a lack of substrate inhibition, high activity with the synthetic coenzyme 3-acetylpyridine adenine dinucleotide (APAD+) and changes in residues at previously conserved positions. pfLDH shows several amino acid insertions and deletions in an alignment with protein sequences from all other known LDHs. The most notable is a five amino acid insertion into the active-site loop. In addition, a conserved serine at position 163 is replaced by leucine. The results showed that when the unique pfLDH structural features were engineered into Bacillus stearothermophilus lactate dehydrogenase, the thermophilic enzyme acquired the properties previously uniquely associated with the malarial enzyme. We conclude that the homology model of the malarial enzyme is adequate for the prediction of successful redesigns and, in the regions tested, is accurate.

Amino Acid Sequence↗

Molecular design of two sterol 14alpha-demethylase homology models and their interactions with the azole antifungals ketoconazole and bifonazole.

Sterol 14alpha-demethylase (CYP51) is one of the known major targets for azole antifungals. Therapeutic side effects of these antifungals are based on interactions of the azoles with the human analogue enzyme. This study describes for the first time a comparison of a human CYP51 (HU-CYP51) homology model with a homology model of the fungal CYP51 of Candida albicans (CA-CYP51). Both models are constructed by using the crystal structure of Mycobacterium tuberculosis MT-CYP51 (PDB code: 1EA1). The binding mode of the azole ketoconazole is investigated in molecular dynamics simulations with the GROMACS force field. The usage of special parameters for the iron azole complex binding is necessary to obtain the correct complex geometry in the active site of the enzyme models. Based on the dynamics simulations it is possible to explain the enantioselectivity of the human enzyme and also to predict the binding mode of the isomers of ketoconazole in the active site of the fungal model.

Amino Acid Sequence↗

Kinases, homology models, and high throughput docking.

With the many protein sequences coming from the genome sequencing projects, it is unlikely that we will ever have an atomic resolution structure of every relevant protein. With high throughput crystallography, however, we will soon have representative structures for the vast majority of protein families. Thus the drug discovery and design process will rely heavily on protein modeling to address issues such as designing combinatorial libraries for an entire class of targets and engineering genome-wide selectivity over a target class. In this study we assess the value of high throughput docking into homology models. To do this we dock a database of random compounds seeded with known inhibitors into homology models of six different kinases. In five of the six cases the known inhibitors were found to be enriched by factors of 4-5 in the top 5% of the overall scored and ranked compounds. Furthermore, in the same five cases the known inhibitors were found to be enriched by factors of 2-3 in the top 5% of the scored and ranked known kinase inhibitors, thus showing that the homology models can pick up some of the crucial selectivity information.

Catalytic Domain↗

Three dimensional structure of the soybean agglutinin Gal/GalNAc complexes by homology modeling.

Complexes of soybean agglutinin (SBA) with galactose (Gal) and N-acetyl galactosamine (GalNAc) have been modeled based on its homology to erythrina corallodendron (EcorL) lectin. The three dimensional structure of SBA-Gal modeled with homology techniques agrees well with SBA-(beta-LacNAc)2Gal-R complex determined by X-ray crystallographic techniques at the beta-sheet regions and the regions where Ca2+ and Mn2+ ions bind. However, significant deviations have been observed between the modeled and the X-ray structures, particularly at the loop regions where the polypeptide chain could not be unequivocally traced in the X-ray structure. The hydrogen bonding scheme, predicted from the homology model, shows that the invariant residues i.e. Asp, Gly, Asn, and aromatic residues (Phe) found in all other legume lectins, bind Gal, slightly in a different way than reported in X-ray structure of SBA-pentasaccharide complex. The higher binding affinity of GalNAc over Gal to SBA is due to additional hydrophobic interactions with Tyr107 rather than a hydrogen bond between N-acetamide group of the sugar and the side chain of Asp88 as suggested from X-ray crystal structure studies. Our modeling also suggest that the variation in the length of the loop D observed among galactose binding legume lectins may not have any effect on the binding of sugar at the monosaccharide specific site of the lectins. Soybean agglutinin (SBA) is a member of the leguminous family of lectins. They generally possess a single carbohydrate binding site, besides the tightly bound Ca2+ and Mn2+ ions which are required for their carbohydrate binding activity. They possess a high degree of sequence homology and about 50% of the amino acid residues are invariant. Some of these invariant amino acid residues are involved in the binding of sugar moieties and in metal ion coordination. X-ray crystallographic studies showed that their three-dimensional structures are very similar, though they differ in their carbohydrate binding specificity (1-6). Three of the invariant residues Asp, Gly, and Asn, besides an aromatic residue (Phe or Tyr), are involved in carbohydrate binding. Independent of their sugar specificity, these four residues in legume lectins provide the basic frame for the sugar to bind.

Acetylgalactosamine↗

A restraint molecular dynamics and simulated annealing approach for protein homology modeling utilizing mean angles.

BACKGROUND: We have developed the program PERMOL for semi-automated homology modeling of proteins. It is based on restrained molecular dynamics using a simulated annealing protocol in torsion angle space. As main restraints defining the optimal local geometry of the structure weighted mean dihedral angles and their standard deviations are used which are calculated with an algorithm described earlier by Doker et al. (1999, BBRC, 257, 348-350). The overall long-range contacts are established via a small number of distance restraints between atoms involved in hydrogen bonds and backbone atoms of conserved residues. Employing the restraints generated by PERMOL three-dimensional structures are obtained using standard molecular dynamics programs such as DYANA or CNS. RESULTS: To test this modeling approach it has been used for predicting the structure of the histidine-containing phosphocarrier protein HPr from E. coli and the structure of the human peroxisome proliferator activated receptor gamma (Ppar gamma). The divergence between the modeled HPr and the previously determined X-ray structure was comparable to the divergence between the X-ray structure and the published NMR structure. The modeled structure of Ppar gamma was also very close to the previously solved X-ray structure with an RMSD of 0.262 nm for the backbone atoms. CONCLUSION: In summary, we present a new method for homology modeling capable of producing high-quality structure models. An advantage of the method is that it can be used in combination with incomplete NMR data to obtain reasonable structure models in accordance with the experimental data.

Algorithms↗

CaspR: a web server for automated molecular replacement using homology modelling.

Molecular replacement (MR) is the method of choice for X-ray crystallography structure determination when structural homologues are available in the Protein Data Bank (PDB). Although the success rate of MR decreases sharply when the sequence similarity between template and target proteins drops below 35% identical residues, it has been found that screening for MR solutions with a large number of different homology models may still produce a suitable solution where the original template failed. Here we present the web tool CaspR, implementing such a strategy in an automated manner. On input of experimental diffraction data, of the corresponding target sequence and of one or several potential templates, CaspR executes an optimized molecular replacement procedure using a combination of well-established stand-alone software tools. The protocol of model building and screening begins with the generation of multiple structure-sequence alignments produced with T-COFFEE, followed by homology model building using MODELLER, molecular replacement with AMoRe and model refinement based on CNS. As a result, CaspR provides a progress report in the form of hierarchically organized summary sheets that describe the different stages of the computation with an increasing level of detail. For the 10 highest-scoring potential solutions, pre-refined structures are made available for download in PDB format. Results already obtained with CaspR and reported on the web server suggest that such a strategy significantly increases the fraction of protein structures which may be solved by MR. Moreover, even in situations where standard MR yields a solution, pre-refined homology models produced by CaspR significantly reduce the time-consuming refinement process. We expect this automated procedure to have a significant impact on the throughput of large-scale structural genomics projects. CaspR is freely available at http://igs-server.cnrs-mrs.fr/Caspr/.

Escherichia coli Proteins↗

Homology modeling of the structure of tobacco acetohydroxy acid synthase and examination of the active site by site-directed mutagenesis.

A reliable model of tobacco acetohydroxy acid synthase (AHAS) was obtained by homology modeling based on a yeast AHAS X-ray structure using the Swiss-Model server. Conserved residues at the dimer interface were identified, of which the functional roles of four residues, namely H142, E143, M489, and M542, were determined by site-directed mutagenesis. Eight mutants were successfully generated and purified, five of which (H142T, M489V, M542C, M542I, and M542V) were found to be inactive under various assay conditions. The H142K mutant was moderately altered in all kinetic parameters to a similar extent. In addition, the mutant was more thermo-labile than wild type enzyme. The E143A mutant increased the Km value more than 20-fold while other parameters were not significantly changed. All mutations carried out on residue M542 inactivated the enzyme. Though showing a single band on SDS-PAGE, the M542C mutant lost its native tertiary structure and was aggregated. Except M542C, each of the other mutants showed a secondary structure similar to that of wild type enzyme. Although all the inactive mutants were able to bind FAD, the mutants M489V and M542C showed a very low affinity for FAD. None of the active mutants constructed was strongly resistant to three tested herbicides. Taken together, the results suggest that the residues of H142, E143, M489, and M542 are essential for catalytic activity. Furthermore, it seems that H142 residue is involved in stabilizing the dimer interaction, while E143 residue may be involved in binding with substrate pyruvate. The data from the site-directed mutagenesis imply that the constructed homology model of tobacco AHAS is realistic.

Acetolactate Synthase↗

Structure/function analyses of human serum paraoxonase (HuPON1) mutants designed from a DFPase-like homology model.

Human serum paraoxonase (HuPON1) is a calcium-dependent enzyme that hydrolyzes esters, including organophosphates and lactones, and exhibits anti-atherogenic properties. A few amino acids have been shown to be essential for the enzyme's arylesterase and organophosphatase activities. Until very recently, a three-dimensional model was not available for HuPON1, so functional roles have not been assigned to those residues. Based on sequence-structure alignment studies, we have folded the amino acid sequence of HuPON1 onto the sixfold beta-propeller structure of squid diisopropylfluorophosphatase (DFPase). We tested the validity of this homology model by circular dichroism (CD) spectroscopy and site-directed mutagenesis. Consistent with predictions from the homology model, CD data indicated that the structural composition of purified HuPON1 consists mainly of beta-sheets. Mutants of HuPON1 were assayed for enzymatic activity against phenyl acetate and paraoxon. Substitution of residues predicted to be important for substrate binding (L69, H134, F222, and C284), calcium ion coordination (D54, N168, N224, and D269), and catalytic mechanism of HuPON1 (H285) led to enzyme inactivation. Mutants F222Y and H115W exhibited substrate-binding selectivity towards phenyl acetate and paraoxon, respectively. The homology model presented here is very similar to the recently obtained PON1 crystal structure, and has allowed identification of several residues within the enzyme active site.

Animals↗

Computational elucidation of the structural basis of ligand binding to the dopamine 3 receptor through docking and homology modeling.

The dopamine subtype 3 receptor (D3) is a promising therapeutic target for the treatment of cocaine addiction, schizophrenia, Parkinson's disease, and other disorders, but little is known about the binding of ligands to D3 at the atomic level. In the present study, binding of 29 known ligands to the D3 receptor was modeled computationally using four D3 receptor models which were obtained from homology modeling. The predicted binding models were validated with experimental data from site-directed mutagenesis, structure-activity relationship studies, and affinity labeling studies. Docking scores calculated for these 29 ligands correlate reasonably well with the experimentally determined binding affinities. A pharmacophore model is proposed that describes the binding of ligands at a single D3 receptor binding site and offers insights into the binding of structurally diverse D3 ligands to this receptor.

Binding Sites↗

Comparative molecular field analysis of substrates for an aryl sulfotransferase based on catalytic mechanism and protein homology modeling.

Comparative Molecular Field Analysis (CoMFA) methods were used to produce a 3D-QSAR model that correlated the catalytic efficiency of rat hepatic aryl sulfotransferase (AST) IV, expressed as log(k(cat)/K(m)), with the molecular structures of its substrates. A total of 35 substrate molecules were used to construct a CoMFA model that was evaluated on the basis of its leave-one-out cross-validated partial least-squares value (q(2)) and its ability to predict the activity of six additional substrates not used in the training set. The model was constructed using substrate conformations that favored (1) proton abstraction by the catalytic histidine residue, (2) an in-line sulfuryl-group transfer mechanism, and (3) constraints imposed by the residues lining the substrate binding pocket of a homology model of AST IV. This CoMFA model had a q(2) value of 0.691, and it successfully predicted the activities of the six molecules not used in the training set. A final CoMFA model was constructed using the same methodology but with molecules from both the training set and the test set. Its q(2) value was 0.701, and it had a non-cross-validated r(2) value of 0.922. The contour coefficient map generated by this CoMFA was overlaid on the amino acids in the substrate-binding pocket of the homology model of AST IV and found to show a good fit. Additionally external validation was obtained by using the CoMFA model to design substrates that show high activities. These results establish a methodology for prediction of the substrate specificity of this sulfotransferase based on CoMFA methods that are guided by both the homology model and the catalytic mechanism of the enzyme.

Animals↗

Structure of human annexin I: comparison of homology modelling and crystallographic experiment.

A model of domain II of annexin I has been built by homology modelling using an annexin V crystal structure as a template. The method used is based on that of Summers and Karplus (J Mol Biol (1989) 210, 785-811) and involves the calculation of torsion-angle rotational energy maps to position side chains. The RMS deviation of the backbone heavy atoms between the model and a crystal structure of annexin I is 1.1 A. Similarities and differences in the experimental and model-derived side-chain rotameric conformations and hydrogen-bonding interactions are examined. It is found that whereas many of the side chains are well positioned some of those placed using the 'entropy argument' in which the broadest of the available minima are preferred, are erroneous. The domain is subjected to molecular dynamics simulation in explicit solvent. The simulations are found to 'correct' some of the side-chain rotamer positions that were poorly placed in the homology modelling. Considerable helix instability is seen in the simulations, consistent with the requirement of domain interactions for the structural integrity of the protein.

Amino Acid Sequence↗

Construction of a model of the Candida albicans lanosterol 14-alpha-demethylase active site using the homology modelling technique.

On the basis of all hitherto known P450 X-ray structures and applying standard homology modelling procedures a three-dimensional model of the lanosterol-14 alpha-demethylase active site was constructed. The modelled active site nicely hosts the natural substrate lanosterol and the substrate-enzyme complex displayed stability in a 70 ps molecular dynamics simulation. The importance of Thr 122 of lanosterol 14 alpha-demethylase for hydrogen bond formation with the 3-hydroxyl group of lanosterol was found to be a characteristic feature of the interaction geometry.

Binding Sites↗