Computational structural proteomics.
Explore the source record for details and available documents.
Animals↗
Biomedical subjects
R Abagyan
Publications and source records attributed to R Abagyan.
In silico discovery of novel retinoic acid receptor agonist structures.
BACKGROUND: Several Retinoic Acid Receptors (RAR) agonists have therapeutic activity against a variety of cancer types; however, unacceptable toxicity profiles have hindered the development of drugs. RAR agonists presenting novel structural and chemical features could therefore open new avenues for the discovery of leads against breast, lung and prostate cancer or leukemia. RESULTS: We have analysed the induced fit of the active site residues upon binding of a known ligand. The derived binding site models were used to dock over 150,000 molecules in silico (or virtually) to the structure of the receptor with the Internal Coordinates Mechanics (ICM) program. Thirty ligand candidates were tested in vitro. CONCLUSIONS: Two novel agonists resulting from the predicted receptor model were active at 50 nM. One of them displays novel structural features which may translate into the development of new ligands for cancer therapy.
Algorithms↗
Rapid boundary element solvation electrostatics calculations in folding simulations: successful folding of a 23-residue peptide.
Solvation effects play a profound role in the energetics of protein folding. While a continuum dielectric model of solvation may provide a sufficiently accurate estimate of the solvation effects, until now this model was too computationally expensive and unstable for folding simulations. Here we proposed a fast yet accurate and robust implementation of the boundary element solution of the Poisson equation, the REBEL algorithm. Using our earlier double-energy scheme, we included for the first time the mathematically rigorous continuous REBEL solvation term in our Biased Probability Monte Carlo (BPMC) simulations of the peptide folding. The free energy of a 23-residue beta beta alpha-peptide was then globally optimized with and without the solvation electrostatics contribution. An ensemble of beta beta alpha conformations was found at and near the global minimum of the energy function with the REBEL electrostatic solvation term. Much poorer correspondence to the native solution structure was found in the "control" simulations with a traditional method to account for solvation via a distance-dependent dielectric constant. Each simulation took less than 40 h (21 h without electrostatic solvation calculation) on a single Alpha 677 MHz CPU and involved more than 40,000 solvation energy evaluations. This work demonstrates for the first time that such a simulation can be performed in a realistic time frame. The proposed procedure may eliminate the energy evaluation accuracy bottleneck in folding simulations.
Algorithms↗
High-throughput docking for lead generation.
Recent improvements in flexible docking technology may lead to a bigger role for computational methods in lead discovery. Although fast and accurate computational prediction of binding affinities for an arbitrary molecule is still beyond the limits of current methods, the docking and screening procedures can select small sets of likely lead candidates from large libraries of either commercially or synthetically available compounds.
Binding Sites↗
Aminoglycoside binding in the major groove of duplex RNA: the thermodynamic and electrostatic forces that govern recognition.
We use a combination of spectroscopic, calorimetric, viscometric and computer modeling techniques to characterize the binding of the aminoglycoside antibiotic, tobramycin, to the polymeric RNA duplex, poly(rI).poly(rC), which exhibits the characteristic A-type conformation that is conserved among natural and synthetic double-helical RNA sequences. Our results reveal the following significant features: (i) CD-detected binding of tobramycin to poly(rI).poly(rC) reveals an apparent site size of four base-pairs per bound drug molecule; (ii) tobramycin binding enhances the thermal stability of the host poly(rI).poly(rC) duplex, the extent of which decreases upon increasing in Na(+) concentration and/or pH conditions; (iii) the enthalpy of tobramycin- poly(rI).poly(rC) complexation increases with increasing pH conditions, an observation consistent with binding-induced protonation of one or more drug amino groups; (iv) the affinity of tobramycin for poly(rI).poly(rC) is sensitive to both pH and Na(+) concentration, with increases in pH and/or Na(+) concentration resulting in a concomitant reduction in binding affinity. The salt dependence of the tobramycin binding affinity reveals that the drug binds to the host RNA duplex as trication. (v) The thermodynamic driving force for tobramycin- poly(rI).poly(rC) complexation depends on pH conditions. Specifically, at pH< or =6.0, tobramycin binding is entropy driven, but is enthalpy driven at pH > 6.0. (vi) Viscometric data reveal non-intercalative binding properties when tobramycin complexes with poly(rI).poly(rC), consistent with a major groove-directed mode of binding. These data also are consistent with a binding-induced reduction in the apparent molecular length of the host RNA duplex. (vii) Computer modeling studies reveal a tobramycin-poly(rI). poly(rC) complex in which the drug fits snugly at the base of the RNA major groove and is stabilized, at least in part, by an array of hydrogen bonding interactions with both base and backbone atoms of the host RNA. These studies also demonstrate an inability of tobramycin to form a stable low-energy complex with the minor groove of the poly(rI).poly(rC) duplex. In the aggregate, our results suggest that tobramycin-RNA recognition is dictated and controlled by a broad range of factors that include electrostatic interactions, hydrogen bonding interactions, drug protonation reactions, and binding-induced alterations in the structure of the host RNA. These modulatory effects on tobramycin-RNA complexation are discussed in terms of their potential importance for the selective recognition of specific RNA structural motifs, such as asymmetric internal loops or hairpin loop-stem junctions, by aminoglycoside antibiotics and their derivatives.
Anti-Bacterial Agents↗
Expression of the adenovirus receptor and its interaction with the fiber knob.
The coxsackievirus group B (CVB) and adenovirus (Ad) receptor (HCVADR, formerly HCAR) is a cell surface protein with two immunoglobulin-like regions (IG1 and IG2) that serves as a receptor for two structurally unrelated viruses. We have established the tissue distribution of the receptor in the rodent by immunohistochemistry and show that the receptor is broadly expressed during embryonic development in the central and peripheral nervous systems and in several types of epithelial cells. The tissue distribution is more restricted in the adult but remains high mainly in epithelial cells. Using site-directed mutagenesis, based on computer modeling of the IG1 region, Ad5 binding could be inhibited but CVB attachment was unaffected. A double amino acid substitution in a three-stranded anti-parallel beta sheet that may form a face of the receptor completely inhibited Ad5 binding. Therefore, we conclude that the molecular interactions critical for Ad5 binding to HCVADR do not overlap with those of CVB3. In fact a specific antibody interfering with only CVB binding recognizes the IG2 domain in the receptor, suggesting that the CVB interacts with this region or an overlap between the IG1 and the IG2 regions.
3T3 Cells↗
Rational discovery of novel nuclear hormone receptor antagonists.
Nuclear hormone receptors (NRs) are potential targets for therapeutic approaches to many clinical conditions, including cancer, diabetes, and neurological diseases. The crystal structure of the ligand binding domain of agonist-bound NRs enables the design of compounds with agonist activity. However, with the exception of the human estrogen receptor-alpha, the lack of antagonist-bound "inactive" receptor structures hinders the rational design of receptor antagonists. In this study, we present a strategy for designing such antagonists. We constructed a model of the inactive conformation of human retinoic acid receptor-alpha by using information derived from antagonist-bound estrogen receptor-alpha and applied a computer-based virtual screening algorithm to identify retinoic acid receptor antagonists. Thus, the currently available crystal structures of NRs may be used for the rational design of antagonists, which could lead to the development of novel drugs for a variety of diseases.
Algorithms↗
Estimating local backbone structural deviation in homology models.
After the atomic coordinates themselves, the most important data in a homology model are the spatial reliability estimates associated with each of the atoms (atom annotation). Recent blind homology modeling predictions have demonstrated that principally correct sequence-structure alignments are achievable to sequence identities as low as 25% [Martin, A.C., MacArthur, M.W., Thornton, J.M., 1997. Assessment of comparative modeling in CASP2. Proteins Suppl(1), 14-28]. The locations and extent of spatial deviations in the backbone between correctly aligned homologous protein structures remained very poorly estimated however, and these errors were the cause of errant loop predictions [Abagyan, R., Batalov, S., Cardozo, T., Totrov, M., Webber, J., Zhou, Y., 1997. Homology modeling with internal coordinate mechanics: deformation zone mapping and improvements of models via conformational search. Proteins Suppl(1), 29-37]. In order to derive accurate measures for local backbone deviations, we made a systematic study of static local backbone deviations between homologous pairs of protein structures. We found that 'through space' proximity to gaps and chain termini, local three-dimensional 'density', three-dimensional environment conservation, and B-factor of the template contribute to local deviations in the backbone in addition to local sequence identity. Based on these finding, we have identified the meaningful ranges of values within which each of these parameters correlates with static local backbone deviation and produced a combined scoring function to greatly improve the estimation of local backbone deviations. The optimized function has more than twice the accuracy of local sequence identity or B-factor alone and was validated in a recent blind structure prediction experiment. This method may be used to evaluate the utility of a preliminary homology model for a particular biological investigation (e.g. drug design) or to provide an improved starting point for molecular mechanics loop prediction methods.
Computer Simulation↗
Genetic linkage mapping of zebrafish genes and ESTs.
Genetic screens in zebrafish (Danio rerio) have isolated mutations in hundreds of genes essential for vertebrate development, physiology, and behavior. We have constructed a genetic linkage map that will facilitate the identification of candidate genes for these mutations and allow comparisons among the genomes of zebrafish and other vertebrates. On this map, we have localized 771 zebrafish genes and expressed sequence tags (ESTs) by scoring single-stranded conformational polymorphisms (SSCPs) in a meiotic mapping panel. Of these sequences, 642 represent previously unmapped genes and ESTs. The mapping panel was comprised of 42 homozygous diploid individuals produced by heat shock treatment of haploid embryos at the one-cell stage (HS diploids). This "doubled haploid" strategy combines the advantages of mapping in haploid and standard diploid systems, because heat shock diploid individuals have only one allele at each locus and can survive to adulthood, enabling a relatively large quantity of genomic DNA to be prepared from each individual in the mapping panel. To integrate this map with others, we also scored 593 previously mapped simple-sequence length polymorphisms (SSLPs) in the mapping panel. This map will accelerate the molecular analysis of zebrafish mutations and facilitate comparative analysis of vertebrate genomes.
Animals↗
Antibodies against thrombospondin-related anonymous protein do not inhibit Plasmodium sporozoite infectivity in vivo.
Thrombospondin-related anonymous protein (TRAP), a candidate malaria vaccine antigen, is required for Plasmodium sporozoite gliding motility and cell invasion. For the first time, the ability of antibodies against TRAP to inhibit sporozoite infectivity in vivo is evaluated in detail. TRAP contains an A-domain, a well-characterized adhesive motif found in integrins. We modeled here a three-dimensional structure of the TRAP A-domain of Plasmodium yoelii and located regions surrounding the MIDAS (metal ion-dependent adhesion site), the presumed business end of the domain. Mice were immunized with constructs containing these A-domain regions but were not protected from sporozoite challenge. Furthermore, monoclonal and rabbit polyclonal antibodies against the A-domain, the conserved N terminus, and the repeat region of TRAP had no effect on the gliding motility or sporozoite infectivity to mice. TRAP is located in micronemes, secretory organelles of apicomplexan parasites. Accordingly, the antibodies tested here stained cytoplasmic TRAP brightly by immunofluorescence. However, very little TRAP could be detected on the surface of sporozoites. In contrast, a dramatic relocalization of TRAP onto the parasite surface occurred when sporozoites were treated with calcium ionophore. This likely mimics the release of TRAP from micronemes when a sporozoite contacts its target cell in vivo. Contact with hepatoma cells in culture also appeared to induce the release of TRAP onto the surface of sporozoites. If large amounts of TRAP are released in close proximity to its cellular receptor(s), effective competitive inhibition by antibodies may be difficult to achieve.
Amino Acid Sequence↗
Soft docking an L and a D peptide to an anticholera toxin antibody using internal coordinate mechanics.
BACKGROUND: The tremendous increase in sequential and structural information is a challenge for computer-assisted modelling to predict the binding modes of interacting biomolecules. One important area is the structural understanding of protein-peptide interactions, information that is increasingly important for the design of biologically active compounds. RESULTS: We predicted the three-dimensional structure of a complex between the monoclonal antibody TE33 and its cholera-toxin-derived peptide epitope VPGSQHID. Using the internal coordinate mechanics (ICM) method of flexible docking, the bound conformation of the initially extended peptide epitope to the antibody crystal or modelled structure reproduced the known binding conformation to a root mean square deviation of between 1.9 A and 3.1 A. The predicted complexes are in good agreement with binding data obtained from substitutional analyses in which each epitope residue is replaced by all other amino acids. Furthermore, a de novo prediction of the recently discovered TE33-binding D peptide dwGsqhydp (single-letter amino acid code where D amino acids are represented by lower-case letters) explains results obtained from binding studies with 172 peptide analogues. CONCLUSIONS: Despite the difficulties arising from the huge conformational space of a peptide, this approach allowed the prediction of the correct binding orientation and the majority of essential binding features of a peptide-antibody complex.
Antibodies, Monoclonal↗
Identification of the individual residues that determine human CD59 species selective activity.
Formation of the cytolytic membrane attack complex of complement on host cells is inhibited by the membrane-bound glycoprotein, CD59. The inhibitory activity of CD59 is species restricted, and human CD59 is not effective against rat complement. Previous functional analysis of chimeric human/rat CD59 proteins indicated that the residues responsible for the species selective function of human CD59 map to a region contained between positions 40 and 66 in the primary structure. By comparative analysis of rat and human CD59 models and by mutational analysis of candidate residues, we now identify the individual residues within the 40-66 region that confer species selective function on human CD59. All nonconserved residues within the 40-66 sequence were substituted from human to rat residues in a series of chimeric human/rat CD59 mutant proteins. Functional analysis revealed that the individual human to rat residue substitutions F47A, T51L, R55E, and K65Q each produced a mutant human CD59 protein with enhanced rat complement inhibitory activity with the single F47A substitution having the most significant effect. Interestingly, the side chains of the residues at positions 47, 51, and 55 are all located on the short single helix (residues 47-55) of CD59 and form an exposed continuous strip parallel to the helix axis. A single human CD59 mutant protein containing rat residue substitutions at all three helix residues produced a protein with species selective activity comparable to that of rat CD59. We further found that synthetic peptides spanning the human CD59 helix sequence were able to inhibit the binding of human CD59 to human C8, but had little effect on the binding of rat CD59 to rat C8.
Amino Acid Sequence↗
NRIF3 is a novel coactivator mediating functional specificity of nuclear hormone receptors.
Many nuclear receptors are capable of recognizing similar DNA elements. The molecular event(s) underlying the functional specificities of these receptors (in regulating the expression of their native target genes) is a very important issue that remains poorly understood. Here we report the cloning and analysis of a novel nuclear receptor coactivator (designated NRIF3) that exhibits a distinct receptor specificity. Fluorescence microscopy shows that NRIF3 localizes to the cell nucleus. The yeast two-hybrid and/or in vitro binding assays indicated that NRIF3 specifically interacts with the thyroid hormone receptor (TR) and retinoid X receptor (RXR) in a ligand-dependent fashion but does not bind to the retinoic acid receptor, vitamin D receptor, progesterone receptor, glucocorticoid receptor, or estrogen receptor. Functional experiments showed that NRIF3 significantly potentiates TR- and RXR-mediated transactivation in vivo but has little effect on other examined nuclear receptors. Domain and mutagenesis analyses indicated that a novel C-terminal domain in NRIF3 plays an essential role in its specific interaction with liganded TR and RXR while the N-terminal LXXLL motif plays a minor role in allowing optimum interaction. Computer modeling and subsequent experimental analysis suggested that the C-terminal domain of NRIF3 directly mediates interaction with liganded receptors through an LXXIL (a variant of the canonical LXXLL) module while the other part of the NRIF3 protein may still play a role in conferring its receptor specificity. Identification of a coactivator with such a unique receptor specificity may provide new insight into the molecular mechanism(s) of receptor-mediated transcriptional activation as well as the functional specificities of nuclear receptors.
Amino Acid Motifs↗
A genetic linkage map for zebrafish: comparative analysis and localization of genes and expressed sequences.
Genetic screens in zebrafish (Danio rerio) have isolated mutations in hundreds of genes with essential functions. To facilitate the identification of candidate genes for these mutations, we have genetically mapped 104 genes and expressed sequence tags by scoring single-strand conformational polymorphisms in a panel of haploid siblings. To integrate this map with existing genetic maps, we also scored 275 previously mapped genes, microsatellites, and sequence-tagged sites in the same haploid panel. Systematic phylogenetic analysis defined likely mammalian orthologs of mapped zebrafish genes, and comparison of map positions in zebrafish and mammals identified significant conservation of synteny. This comparative analysis also identified pairs of zebrafish genes that appear to be orthologous to single mammalian genes, suggesting that these genes arose in a genome duplication that occurred in the teleost lineage after the divergence of fish and mammal ancestors. This comparative map analysis will be useful in predicting the locations of zebrafish genes from mammalian gene maps and in understanding the evolution of the vertebrate genome.
Animals↗
Identification and analysis of PH domain-containing targets of phosphatidylinositol 3-kinase using a novel in vivo assay in yeast.
Phosphatidylinositol 3-kinase (PI3K) mediates a variety of cellular responses by generating PtdIns(3,4)P2 and PtdIns(3,4,5)P3. These 3-phosphoinositides then function directly as second messengers to activate downstream signaling molecules by binding pleckstrin homology (PH) domains in these signaling molecules. We have established a novel assay in the yeast Saccharomyces cerevisiae to identify proteins that bind PtdIns(3,4)P2 and PtdIns(3,4,5)P3 in vivo which we have called TOPIS (Targets of PI3K Identification System). The assay uses a plasma membrane-targeted Ras to complement a temperature-sensitive CDC25 Ras exchange factor in yeast. Coexpression of PI3K and a fusion protein of activated Ras joined to a PH domain known to bind PtdIns(3,4)P2 (AKT) or PtdIns(3,4,5)P3 (BTK) rescues yeast growth at the non-permissive temperature of 37 degreesC. Using this assay, we have identified several amino acids in the beta1-beta2 region of PH domains that are critical for high affinity binding to PtdIns(3,4)P2 and/or PtdIns(3,4,5)P3, and we have proposed a structural model for how these PH domains might bind PI3K products with high affinity. From these data, we derived a consensus sequence which predicts high-affinity binding to PtdIns(3, 4)P2 and/or PtdIns(3,4,5)P3, and we have identified several new PH domain-containing proteins that bind PI3K products, including Gab1, Dos, myosinX, and Sbf1. Use of this assay to screen for novel cDNAs which rescue yeast at the non-permissive temperature should provide a powerful approach for uncovering additional targets of PI3K.
Amino Acid Sequence↗
Novel anchorage of GluR2/3 to the postsynaptic density by the AMPA receptor-binding protein ABP.
We report the cloning of alpha-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid (AMPA) receptor-binding protein (ABP), a postsynaptic density (PSD) protein related to glutamate receptor-interacting protein (GRIP) with two sets of three PDZ domains, which binds the GluR2/3 AMPA receptor subunits. ABP exhibits widespread CNS expression and is found at the postsynaptic membrane. We show that the protein interactions of the ABP/GRIP family differ from the PSD-95 family, which binds N-methyl-D-aspartate (NMDA) receptors. ABP binds to the GluR2/3 C-terminal VKI-COOH motif via class II hydrophobic PDZ interactions, distinct from the class I PSD-95-NMDA receptor interaction. ABP and GRIP also form homo- and heteromultimers through PDZ-PDZ interactions but do not bind PSD-95. We suggest that the ABP/GRIP and PSD-95 families form distinct scaffolds that anchor, respectively, AMPA and NMDA receptors.
Amino Acid Sequence↗
Energy strain in three-dimensional protein structures.
BACKGROUND: Steric strain in protein three-dimensional structures is related to unfavorable inter-atomic interactions. The steric strain may be a result of packing or functional requirements, or may indicate an error in the coordinates of a structure. Detailed energy functions are, however, usually considered too noisy for error detection. RESULTS: After a short energy refinement, a full-atom, detailed energy function becomes a sensitive indicator of errors. The statistics of the energy distribution of amino acid residues in high-resolution crystal structures, represented by models with idealized covalent geometry, were calculated. The interaction energy of each residue with the whole protein structure and with the solvent was considered. Normalized deviations of amino acid residue energies from their average values were used for detecting energy-strained and, therefore, potentially incorrect fragments of a polypeptide chain. Protein three-dimensional structures of different origin (X-ray crystallography, NMR spectroscopy, theoretical models and deliberately misfolded decoys) were compared. Examples of the applications to loop and homology modeling are provided. CONCLUSIONS: Elevated levels of energy strain may point at a problematic fragment in a protein three-dimensional structure of either experimental or theoretical origin. The approach may be useful in model building and refinement, modeling by homology, protein design, folding calculations, and protein structure analysis.
Amino Acid Sequence↗
How and why phosphotyrosine-containing peptides bind to the SH2 and PTB domains.
BACKGROUND: Specific recognition of phosphotyrosine-containing protein segments by Src homology 2 (SH2) and phosphotyrosine-binding (PTB) domains plays an important role in intracellular signal transduction. Although many SH2/PTB-domain-containing receptor-peptide complex structures have been solved, little has been done to study the problem computationally. Prediction of the binding geometry and the binding constant of any peptide-protein pair is an extremely important problem. RESULTS: A procedure to predict binding energies of phosphotyrosine-containing peptides with SH2/PTB domains was developed. The average deviation between experimentally measured binding energies and theoretical evaluations was 1.8 kcal/mol. Binding states of unphosphorylated peptides were also predicted reasonably well. Ab initio predictions of binding geometry of fully flexible peptides correctly identified conformations of two pentapeptides and a hexapeptide complexed with a v-Src SH2 domain receptor with root mean square deviations (rmsds) of 0.3 A, 1.2 A and 1.5 A, respectively. CONCLUSIONS: The binding energies of phosphotyrosine-containing complexes can be effectively predicted using the procedure developed here. It was also possible to predict the bound conformations of flexible short peptides correctly from random starting conformations.
Animals↗