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 217 records · Page 12Linked to original sources
PubMed1993

Tertiary structure of calcineurin B by homology modeling.

The crystal structure of the calcium-binding protein calmodulin is used to model the immunologically important calcineurin subunit B. The rough structure is produced by computer-aided homology modeling. Refinement of this using molecular dynamics leads to a suggested structure which appears to satisfy reasonable hydrophilicity and hydrogen-bonding criteria. In the absence of a crystal structure, the model may prove useful in modeling of its interactions with the phosphatase catalytic subunit calcineurin A, and help to explain the calcium modulation of this protein.

Amino Acid Sequence↗
PubMed2003

Structure of Plasmodium vivax dihydrofolate reductase determined by homology modeling and molecular dynamics refinement.

The structure of Plasmodium vivax dihydrofolate reductase (PvDHFR), a potentially important target for antimalarial chemotherapy, was determined by means of homology modeling and molecular dynamics refinement. The structure proved to be consistent with DHFRs of known crystal structure. The comparison of the complexes of the antifolate inhibitor pyrimethamine bound at the active sites of PvDHFR and PfDHFR, the related enzyme from Plasmodium falciparum, prospected the possibility of using structure-based drug design to develop inhibitors that are effective against both malarial enzymes. This study constitutes a first step toward understanding of the antifolate-PvDHFR molecular interactions and possible rationalization of resistance in vivax malaria.

Amino Acid Sequence↗
PubMed2006

Homology modeling and molecular dynamics study of West Nile virus NS3 protease: a molecular basis for the catalytic activity increased by the NS2B cofactor.

The West Nile virus (WNV) NS3 serine protease, which plays an important role in assembly of infective virion, is an attractive target for anti-WNV drug development. Cofactors NS2B and NS4A increase the catalytic activity of NS3 in dengue virus and Hepatitis C virus, respectively. Recent studies on the WNV-NS3 characterize the catalytically active form of NS3 by tethering the 40-residue cofactor NS2B. It is suggested that NS2B is essential for the NS3 activity in WNV, while there is no information of the WNV-NS3-related crystal structure. To understand the role of NS2B/substrate in the NS3 catalytic activity, we built a series of models: WNV-NS3 and WNV-NS3-NS2B and WNV-NS3-NS2B-substrate using homology modeling and molecular modeling techniques. Molecular dynamics (MD) simulations were performed for 2.75 ns on each model, to investigate the structural stabilization and catalytic triad motion of the WNV NS3 protease with and without NS2B/substrate. The simulations show that the NS3 rearrangement occurs upon the NS2B binding, resulting in the stable D75-OD1...H51-NH hydrogen bonding. After the substrate binds to the NS3-NS2B active site, the NS3 protease becomes more stable, and the catalytic triad is formed. These results provide a structural basis for the activation and stabilization of the enzyme by its cofactor and substrate.

Amino Acid Sequence↗
PubMed2001

Homology modeling using multiple molecular dynamics simulations and docking studies of the human androgen receptor ligand binding domain bound to testosterone and nonsteroidal ligands.

To facilitate the rational design of novel and more potent androgen receptor ligands, three-dimensional models for the human androgen receptor ligand binding domain bound to testosterone have been developed. These models of the androgen receptor were based on the crystal structure of the highly homologous human progesterone receptor ligand binding domain. The homology modeled androgen receptor was refined using unrestrained multiple molecular dynamics simulations in explicit solvent. Key H-bonding partners with the 17-hydroxy group and 3-keto group of testosterone are Asn705 and Thr877, and Gln711 and Arg752, respectively. These models show the presence of a unique unoccupied cavity within the androgen receptor binding pocket which may be valuable in the development of novel selective androgen receptor ligands. A qualitative analysis of amino acid mutations within the hAR binding pocket that affect ligand binding are consistent with these androgen receptor models. In addition to testosterone, the binding modes of several hydroxyflutamide-like nonsteroidal ligands for the androgen receptor are investigated using flexible docking with FlexX followed by refinement of the initial complexes with molecular dynamics simulations. These docking studies indicate that Asn705 is an important determinant in binding hydroxyflutamide and its derivatives by participating in H-bond interactions with the alpha-hydroxy moiety of these ligands. In addition, the nitro functionality mimics the 3-keto group of the natural ligand testosterone and is involved in H-bonding interactions with Gln711 and Arg752. From these docking studies, we suggest a mechanism for the enantioselective binding of chiral hydroxyflutamide derivatives and expand upon the previously reported structure-activity relationship for hydroxyflutamide and its derivatives.

Amino Acid Substitution↗
PubMed1994

Homology modelling of annexin I: implicit solvation improves side-chain prediction and combination of evaluation criteria allows recognition of different types of conformational error.

Annexin I homology models were built from the annexin V crystal structure. Three methods for side-chain prediction were tested based on molecular mechanics conformational search, the use of a rotamer database, or a combination of these two methods. We showed that rotamer-based methods were more efficient and that molecular mechanics energy minimizations, prior to rotamer selection, did not afford clearly improved predictions. Models built in vacuo and with an implicit solvation term were compared with the annexin I crystal structure which became available during the course of this study. The analysis of solvation energies, root mean square deviations, chi 1 angles and hydrogen bonds showed that models built with implicit solvation were of better quality. In annexin V, repeat III displays A-B and D-E loop conformations quite different from other repeats. Since the sequence differences suggest that repeat III in annexin I might present a conformation similar to other repeats, two annexin I models with different repeat III conformations were built and compared to determine whether the correct conformation could have been predicted. We show that using a combination of evaluation criteria, it is possible to discriminate unequivocally between the native and the incorrect fold, stressing that only one criterion should not be used to evaluate protein structures.

Amino Acid Sequence↗
PubMed2004

Splice variants: a homology modeling approach.

Splice variants play an important role within the cell in both increasing the proteome diversity and in cellular function. Splice variants are also associated with disease states and may play a role in their etiology. Information about splice variants has, until now, mostly been derived from the primary transcript or through cellular studies. In this study information from the transcript and other studies is combined with tertiary structure information derived from homology models. Through this method we have determined that it is possible to effectively model splice variants. Forty models of splice variants for fourteen proteins were produced. Analysis of the models shows that deletions produce superior model validation values. Additions to sequences where there is little homology become increasingly difficult to model with increasing sequence length. Many of the splicing events are associated with post-translational modification either in the N-terminal region by changing the signal peptide or by affecting the number or availability of glycosylation sites. Often the alternative exon combinations are associated with loss or gain of whole structural units, as opposed to just changing small loop regions. Losing part of the secondary structure may destabilize neighboring parts of the same secondary structure. Detailed analysis is given of four biomedically relevant proteins (Beta-site Amyloid Precursor Protein Cleaving enzyme (BACE), Interleukin-4, Frataxin and Hereditary hemochromatosis protein) and their associated splice variant models. The visualization of these possible structures provides new insights about their functionality and the possible etiology of associated diseases.

Alternative Splicing↗
PubMed2003

Homology modelling of CYP2A6 based on the CYP2C5 crystallographic template: enzyme-substrate interactions and QSARs for binding affinity and inhibition.

The results of homology modelling of the human P450 enzyme CYP2A6, based on the CYP2C5 crystallographic template structure are reported. A substantial number of selective substrates of the CYP2A6 enzyme fit the putative active site in a manner that is consistent with their known metabolites. Moreover, the evidence from site-directed mutagenesis experiments is in accordance with the current model, particularly in relation to complementary amino acid contacts within the haem environment. The binding of substrates is rationalized in terms of QSAR analyses and from a consideration of the contributory factors affecting the binding affinity. The latter approach appears to represent a highly correlated (R=0.99) method for estimating the relative strength of enzyme-substrate binding within CYP2A6-selective compounds, albeit within a fairly limited dataset of substrates.

Amino Acid Sequence↗
PubMedDate not supplied

Amino acid similarity matrix for homology modeling derived from structural alignment and optimized by the Monte Carlo method.

In this paper, we obtained a similarity matrix for homology modeling based on the structure of proteins in a structural alignment. The alignment procedure was executed within dynamic programming generally used in alignment methods. An initial matrix derived from the structural alignment was optimized by the Markov chain Monte Carlo method at low temperature to fit its sequence alignment to the structural alignment. Structural alignment was performed on the basis of the superposition of C alpha atoms for two protein structures. The objective function in the Monte Carlo procedure was defined by entropy in the information theory, allowing us to show that the amino acid similarity matrix aligned accurately. When compared with the structural alignment, the average number of incorrect amino acid residues in the sequence alignment was 22.6 for all residues and about 3.7 for residues in structurally conserved regions. The alignment with our matrix was more similar to structural alignment than to sequence alignments using other amino acid substitution matrices.

Algorithms↗
PubMed1996

Testing homology modeling on mutant proteins: predicting structural and thermodynamic effects in the Ala98-->Val mutants of T4 lysozyme.

BACKGROUND: Current approaches to homology modeling predict how amino acid substitutions will alter a protein's structure, primarily by modeling sidechain conformations upon essentially immobile backbone frameworks. However, recent crystal structures of T4 lysozyme mutants reveal significant shifts of the mainchain and other potentially serious problems for sidechain rotamer-based modeling. This paper evaluates the accuracy of structural and thermodynamic predictions from two common sidechain modeling approaches to measure errors caused by the fixed-backbone approximation. RESULTS: Tested on a series of T4 lysozyme mutants, this sidechain rotamer library approach did not handle mainchain shifts well, correctly predicting the sidechain conformations of only two of six mutants. By contrast, allowing sidechains to move more flexibly appeared to compensate for the rigidity of the mainchain and gave reasonably accurate coordinate predictions (rms errors of 0.5-1.0 A for each mutated sidechain), better on average than 90% of possible conformations. The calculated packing energies correlated well with experimental stabilities (r2 = 0.81) and correctly captured the cooperative interactions of several neighboring mutations. CONCLUSIONS: Mutant modeling can be relatively accurate despite the fixed-backbone approximation. Mainchain shifts (0.2-0.5 A) cause increased sidechain coordinate errors of 0.1-0.8 A, torsional errors of 10-30 degrees, and exaggerated strain energy for overpacked mutants, compared with the same calculations performed with the correct mutant backbones.

Bacteriophage T4↗
PubMed1995

Testing homology modeling on mutant proteins: predicting structural and thermodynamic effects in the Ala98-->Val mutants of T4 lysozyme.

Backgound. Current approaches to homology modeling predict how amino acid substitutions will alter a protein's structure, primarily by modeling sidechain conformations upon essentially immobile backbone frameworks. However, recent crystal structures of T4 lysozyme mutants reveal significant shifts of the mainchain and other potentially serious problems for sidechain rotamer-based modeling. This paper evaluates the accuracy of structural and thermodynamic predictions from two common sidechain modeling approaches to measure errors caused by the fixed-backbone approximation. Results. Tested on a series of T4 lysozyme mutants, this sidechain rotamer library approach did not handle mainchain shifts well, correctly predicting the sidechain conformations of only two of six mutants. By contrast, allowing sidechains to move more flexibly appeared to compensate for the rigidity of the mainchain and gave reasonably accurate coordinate predictions (rms errors of 0.5-1.0 Å for each mutated sidechain), better on average than 90% of possible conformations. The calculated packing energies correlated well with experimental stabilities (r2=0.81) and correctly captured the cooperative interactions of several neighboring mutations. Conclusion. Mutant modeling can be relatively accurate despite the fixed-backbone approximation. Mainchain shifts (0.2-0.5 Å) cause increased sidechain coordinate errors of 0.1-0.8 Å, torsional errors of 10-30°, and exaggerated strain energy for overpacked mutants, compared with the same calculations performed with the correct mutant backbones.

Journal Article↗
PubMed1998

Homology models of the C domains of blood coagulation factors V and VIII: a proposed membrane binding mode for FV and FVIII C2 domains.

We present homology models of the C domains of coagulation factors V (FV) and VIII (FVIII). Using a threading approach, we identified the binding domain of galactose oxidase as an appropriate template for each C domain. The C1 and C2 domains of FV associate to form an elongated cylinder of 80A long and 30A diameter. The folding unit is a beta-sandwich with a long axis of 40A and a diameter of 30A. The current model allows us to propose a membrane binding mode for the C2 domains of FV and FVIII with three major characteristics: 1) solvent-exposed hydrophobic side chains from three loops at one end of the beta-sandwich are buried in the hydrophobic layer of the outer phospholipid leaflet; 2) a crown of positively charged residues is located in the polar zone of the phospholipid head groups; and 3) the long axis of the beta-sandwich of the C2 domain is perpendicular to the plane of the membrane. This proposal satisfies experimentally observed characteristics of membrane binding for the C2 domain and the light chain of FVa.

Amino Acid Sequence↗
PubMed2005

Ligand docking in the gastric H+/K+-ATPase: homology modeling of reversible inhibitor binding sites.

Using the recent high-resolution X-ray structures determined for the Ca2+-ATPase, we have generated two homology models of the gastric H+/K+-ATPase reflecting the E1 and E2 conformations adopted by P-type ATPases in their catalytic cycle. In regimes where the in situ solid-state NMR-determined structure for 1,2,3-trimethyl-8-(pentafluorophenylmethoxy)imidazo[1,2-a]pyridinium iodide (TMPFPIP), a reversible inhibitor of the gastric H+/K+-ATPase, was retained in its predefined conformation and was allowed full torsional flexibility in docking, the ligands localized to discrete binding volumes in the E1 model and to a single central binding space, together with secondary peripheral locations, in the E2 conformation. The results of these binding studies are in good agreement with current site-directed mutagenesis data and support the suggestion that the binding site is proximal to the loop between TM5 and TM6 and TM8, the transmembrane (TM) region considered important for cation translocation. Furthermore, the results of the simulation with the flexible ligand complement the solid-state NMR structural constraints of this inhibitor when bound in situ to the protein.

Binding Sites↗
PubMed1995

A homology model of the three-dimensional structure of human O6-alkylguanine-DNA alkyltransferase based on the crystal structure of the C-terminal domain of the Ada protein from Escherichia coli.

O6-Alkylguanine-DNA alkyltransferase (EC 2.1.1.63) repairs O6-alkylguanine lesions in DNA. A homology model of the human protein (hAT) was built, based on the crystal structure of the C-terminal domain of the Ada protein, which carries out a similar repair in Escherichia coli. Sequence alignments of known O6-alkylguanine-DNA alkyltransferases were used to aid the model building using QUANTA and CHARMm software. Despite low homology in the N-terminal half (hAT residues 1-85), a well-defined topology over this region in Ada permitted successful modelling. The C-terminal half of hAT (residues 92-207) was modelled almost entirely by residue-for-residue superposition onto the Ada structure up to residue hAT175. The model was solvated to a residue radius of 8.0 A [corrected] and then minimized using CHARMm. This structural model was used to rationalize findings from site-directed mutagenesis experiments on hAT, to make further predictions on the relationship between structure and function for the alkyltransferase family of proteins, and to explain the specificity towards known small-molecule inhibitors of the protein.

Amino Acid Sequence↗
PubMed1994

Homology modelling of the dihydrofolate reductase-thymidylate synthase bifunctional enzyme of Leishmania major, a potential target for rational drug design in leishmaniasis.

A homology model of the dimeric Leishmania major dihydrofolate reductase-thymidylate synthase bifunctional enzyme was built based upon the available crystallographic structures of the dihydrofolate reductase and thymidylate synthase monofunctional enzymes. The corresponding substrates, cofactors and inhibitors of the monofunctional enzymes were modelled in the active-sites of the homodimeric model of the bifunctional enzyme. The small number of residues in the interdomain region of the L. major sequence, compared to the other DHFR-TS sequences, imposed severe constraints on the relative locations of the two functional domains. The DHFR and TS domains in the malarial, trypanosomal and toxoplasmia structures are likely to have comparable locations. The leishmanial model reveals features which may be of use in designing novel antifolates for leishmaniasis and protozoal parasites in general.

Amino Acid Sequence↗
PubMed1998

Evidence for subsites in the galectins involved in sugar binding at the nonreducing end of the central galactose of oligosaccharide ligands: sequence analysis, homology modeling and mutagenesis studies of hamster galectin-3.

A model of the carbohydrate recognition domain CRD, residues 111-245, of hamster galectin-3 has been made using homology modeling and dynamics minimization methods. The model is based on the known x-ray structures of bovine galectin-1 and human galectin-2. The oligosaccharides NeuNAc-alpha2,3-Gal-beta1,4-Glc and GalNAc-alpha1, 3-[Fuc-alpha1,2]-Gal-beta1,4-Glc, known to be specific high-affinity ligands for galectin-3, as well as lactose recognized by all galectins were docked in the galectin-3 CRD model structure and a minimized binding conformation found in each case. These studies indicate a putative extended carbohydrate-binding subsite in the hamster galectin-3 involving Arg139, Glu230, and Ser232 for NeuNAc-alpha2,3-; Arg139 and Glu160 for fucose-alpha1,2-; and Arg139 and Ile141 for GalNAc-alpha1,3- substituents on the primary galactose. Each of these positions is variable within the whole galectin family. Two of these residues, Arg139 and Ser232, were selected for mutagenesis to probe their importance in this newly identified putative subsite. Residue 139 adopts main-chain dihedral angles characteristic of an isolated bridge structural feature, while residue 232 is the C-terminal residue of beta-strand-11, and is followed immediately by an inverse gamma-turn. A systematic series of mutant proteins have been prepared to represent the residue variation present in the aligned sequences of galectins-1, -2, and -3. Minimized docked models were generated for each mutant in complex with NeuNAc-alpha2,3-Gal-beta1,4-Glc, GalNAc-alpha1, 3-[Fuc-alpha1,2]-Gal-beta1,4- Glc, and Gal-beta1,4-Glc. Correlation of the computed protein-carbohydrate interaction energies for each lectin-oligosaccharide pair with the experimentally determined binding affinities for fetuin and asialofetuin or the relative potencies of lactose and sialyllactose in inhibiting binding to asiolofetuin is consistent with the postulated key importance of Arg139 in recognition of the extended sialylated ligand.

Amino Acid Sequence↗
PubMed2002

Homology modeling of the cellulose-binding domain of a pollen allergen from rye grass: structural basis for the cellulose recognition and associated allergenic properties.

A three-dimensional model of the cellulose-binding domain of the rye-grass pollen allergen Lol pI built by homology modeling is proposed as a structural scaffold for expansins and other expansin-related proteins. A groove and an extended strip of aromatic and polar residues presumably account for the cellulose-binding properties of the protein domain. Two of the four predicted T-cell epitopes readily exposed on the surface of the cellulose-binding domain match with previously reported IgE-binding regions. A close structural relationship occurs between the cellulose-binding and allergenic properties.

Allergens↗
PubMed2006

C-(4,5,6-trimethoxyindan-1-yl)methanamine: a mescaline analogue designed using a homology model of the 5-HT2A receptor.

A conformationally restricted analogue of mescaline, C-(4,5,6-trimethoxyindan-1-yl)-methanamine, was designed using a 5-HT(2A) receptor homology model. The compound possessed 3-fold higher affinity and potency than and efficacy equal to that of mescaline at the 5-HT(2A) receptor. The new analogue substituted fully for LSD in drug discrimination studies and was 5-fold more potent than mescaline. Resolution of this analogue into its enantiomers corroborated the docking experiments, showing the R-(+) isomer to have higher affinity and potency and to have efficacy similar to that of mescaline at the 5-HT(2A) receptor.

Animals↗
PubMed1995

Structural basis for the extreme thermostability of D-glyceraldehyde-3-phosphate dehydrogenase from Thermotoga maritima: analysis based on homology modelling.

D-Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) from a hyperthermophilic eubacterium, Thermotoga maritima, is remarkably heat stable (Tm = 109 degrees C). In this work, we have applied homology modelling to predict the 3-D structure of Th.maritima GAPDH to reveal the structural basis of thermostability. Three known GAPDH structures were used as reference proteins. First, the rough model of one subunit was constructed using the identified structurally conserved and variable regions of the reference proteins. The holoenzyme was assembled from four subunits and the NAD molecules. The structure was refined by energy minimization and molecular dynamics simulated annealing. No errors were detected in the refined model using the 3-D profile method. The model was compared with the structure of Bacillus stearothermophilus GAPDH to identify structural details underlying the increased thermostability. In all, 12 extra ion pairs per subunit were found at the protein surface. This seems to be the most important factor responsible for thermostability. Differences in the non-specific interactions, including hydration effects, were also found. Minor changes were detected in the secondary structure. The model predicts that a slight increase in alpha-helical propensities and helix-dipole interactions also contribute to increased stability, but to a lesser degree.

Amino Acid Sequence↗