In silico identification and analysis of the binding site for aminocoumarin type inhibitors in the C-terminal domain of Hsp90

Hsp90 contains two Nucleotide Binding Sites (NBS): one each in its N-terminal domain (NTD) and C-terminal domain (CTD), respectively. Previously we used computational techniques to locate a nucleotide-binding site in the CTD. Nucleotide binding at this site stabilized the structurally labile region within this domain, thus providing a rationale for increased resistance to thermal denaturation and proteolysis. A scan for ligand-binding sites in CTD revealed four potential sites with the requisite volume to accommodate aminocoumarins and -derived inhibitors. Only one of these reproducibly formed docked complexes with inhibitors and showed excellent interactions with residues lining the site. Fortuitously, it was identical to the aforementioned nucleotide-binding site thus providing an explanation for the reported direct competition between inhibitors and nucleotides. Further studies with carefully chosen inhibitors and some inactive analogues provided an explanation for the known Structure-Activity Relationships (SAR) of aminocoumarin and -derived inhibitors. We also performed similar studies of the NTD to discern the reason(s) for its inability to bind aminocoumarins, given the family resemblance to prokaryotic Top-IV and Gyr-B. Our studies permitted the identification of the putative inhibitor binding site in the CTD, an explanation for increased resistance to thermal denaturation and proteolysis upon inhibitor binding as well as direct competition with ATP.

Molecular cloning, characterization and expression analysis of heat shock protein 90 in albino northern snakehead Channa argus

The great albino northern snakehead Channa argus is habitual to only the Sichuan Jialing Rivers System in China, making its introduction difficult to other riverine systems. Here, we characterized heat shock protein 90 (AcaHSP90) and probed its molecular responses toward the environmental stressors that C. argus can face during its introduction and breeding in the other southern latitudes of China. To serve the purpose, cDNA encoding of AcaHSP90 were cloned and characterized in albino C. argus. The cDNA was 2752bps that contained an open reading frame (ORF), encoding a 726-amino-acid polypeptide of 83.35kDa (theoretical isoelectric point [pI]: 4.89). Genomic DNA analysis showed that the AcaHSP90 gene consisted of 7 introns, five conserved amino acid blocks and other motifs or domains. The AcaHSP90 structure was highly similar with the other known HSP90s except those identified in the bacteria. The expression profiles of AcaHSP90 gene in albino C. argus were also investigated after experimentally exposed to different temperature stresses (8.5, 26 and 37°C) and infected with Edwardsiella tarda (strain NO. DL1476) at different time intervals (0, 6, 12, 24, 36, 48, 72h). In addition, the AcaHSP90 expression in different tissues of albino C. argus were also analyzed. The quantitative real-time PCR and western blot analysis revealed tissue-specific AcaHSP90 expressions in control group, and expressions were significantly stimulated in the brain, heart, kidney, liver, muscle and spleen after the heat shock (37°C), while showed no significant difference after the cold treatment (8.5°C). The mRNA levels of AcaHSP90 were also significantly upregulated in the spleen and muscle at 12h and in the kidney at 12 and 48h post pathogen injections. In a nut shell, these novel results showed tissue-specific responses of AcaHSP90 and indicated that this heat shock protein might also be sensitive to pathogen infection, but closely related to the thermal resistance in albino C. argus.

In silico identification and computational analysis of the nucleotide binding site in the C-terminal domain of Hsp90

Hsp90 contains two distinct Nucleotide Binding Sites (NBS), in its N-terminal domain (NTD) and C-terminal domain (CTD), respectively. The NTD site belongs to the GHKL super-family of ATPases and has been the subject of extensive characterization. However, a structure of the nucleotide-bound form of CTD is still unavailable. In this study molecular modeling was employed to incorporate experimental data using partial constructs of the CTD, from work published by many research groups, onto existing structural models of its apo- form. Our attempts to locate potential nucleotide ligand-binding sites or cavities yielded one major candidate-a structurally unconventional site-exhibiting the requisite shape and volume for accommodation of tri-phosphate nucleotides. Its structure was refined by molecular dynamics (MD)-based techniques. We reproducibly docked the Mg²⁺ complexed form of ATP, GTP, CTP, TTP and UTP to this putative NBS. These docking simulations and calculated ligand-binding scores are in general agreement with published data about experimentally measured binding to the CTD. The overall pattern of interactions between residues lining the site and docked nucleotides is conserved and broadly similar to that of other nucleotide-binding sites. Our docking simulations suggest that nucleotide binding stabilizes the only structurally labile region, thereby providing a rationale for the increased resistance to thermal denaturation and proteolysis. The docked nucleotides do not intrude onto the surface of residues involved in dimerization or chaperoning. Our molecular modeling permitted recognition of larger structural changes in the nucleotide-bound CTD dimer, including stabilization of helix-2 in both chains and intra- and inter- chain interactions between three residues (I613, Q617, R620).

Molecular paleoscience: systems biology from the past

Experimental paleomolecular biology, paleobiochemistry, and paleogenetics are closely related emerging fields that infer the sequences of ancient genes and proteins from now-extinct organisms, and then resurrect them for study in the laboratory. The goal of paleogenetics is to use information from natural history to solve the conundrum of modern genomics: How can we understand deeply the function of biomolecular structures uncovered and described by modern chemical biology? Reviewed here are the first 20 cases where biomolecular resurrections have been achieved. These show how paleogenetics can lead to an understanding of the function of biomolecules, analyze changing function, and put meaning to genomic sequences, all in ways that are not possible with traditional molecular biological studies.

Structural models for the protein family characterized by gamete surface protein Pfs230 of Plasmodium falciparum

Ps230 is the largest representative of a 10-member family of proteins found in all Plasmodium species. The family is defined by partially conserved, cysteine-rich double domains that are approximately 350 aa in length and have one to three predicted disulfide bridges in each half. In Plasmodium falciparum, the most dangerous human malaria, Pf12 is the smallest member of the family, comprising just one double domain. Pfs230, with 7 double domains, and Pfs48/45 and Pfs47, with 1.5 double domains each, are found on the gamete surfaces and are thus potential candidates for a transmission-blocking vaccine. Fold prediction analyses of the double domains in Pfs230 reveal structural resemblance to SAG1 (surface antigen 1), a surface protein with a double beta-sandwich structure from another apicomplexan parasite, Toxoplasma gondii. Template-directed modeling onto SAG1 clearly establishes the structural link between SAG1 and Pfs230 and produces positions for the cysteines that accord with the disulfide-bonding arrangement predicted for the Pfs230 family in earlier work. A highly clustered region of polymorphisms within the second double domain in Pfs230 maps to one side of the sandwich surface. This observation suggests that this region may be functional and reinforces the validity of these molecular models for the core domains of the Pfs230 family of proteins.

Contributions of residue pairing to beta-sheet formation: conservation and covariation of amino acid residue pairs on antiparallel beta-strands

In an effort to better understand beta-sheet assembly, we have investigated the evolutionary behavior of neighboring residues on adjacent antiparallel beta-strands. Residue pairs were classified according to solvent exposure as well as by whether their backbone NH and C==O groups are hydrogen bonded. The conservation and covariation of 19,241 pairs in 219 sequence alignments was analyzed. Buried pairs were found to be the most conserved, while stronger covariation was detected in the solvent-exposed pairs. However, residues on neighboring strands showed a degree of conservation and covariation similar to that of well-separated residues on the same strand, suggesting that evolutionary pressure to maintain complementarity between pairs on neighboring strands is weak. Moreover, in spite of the preference of certain amino acid pairs to occupy neighboring positions on adjacent strands, such favored pairs are neither more strongly mutually conserved nor covary more strongly than pairs of the same type in non-interacting positions. Although the beta-sheet pairs did not show outstanding evolutionary coupling, in many protein families significant conservation and covariation patterns were detected for some of the residue pairs. Overall, the weak evolutionary conservation and covariation of the beta-sheet pairs indicates that sheet structure is unlikely to be dictated by specific side-chain interactions.

Analysis of covariation in an SH3 domain sequence alignment: applications in tertiary contact prediction and the design of compensating hydrophobic core substitutions

We have analyzed sequence covariation in an alignment of 266 non-redundant SH3 domain sequences using chi-squared statistical methods. Artifactual covariations arising from close evolutionary relationships among certain sequence subgroups were eliminated using empirically derived sequence diversity thresholds. This covariation detection method was able to predict residue-residue contacts (side-chain centres of mass within 8 A) in the structure of the SH3 domain with an accuracy of 85 %, which is greater than that achieved in many previous covariation studies. In examining the positions involved most frequently in covariations, we discovered a dramatic over-representation of a subset of five hydrophobic core positions. This covariation information was used to design second and third site substitutions that could compensate for highly destabilizing hydrophobic core substitutions in the Fyn SH3 domain, thus providing experimental data to validate the covariation analysis. The testing of our covariation detection method on 15 other alignments showed that the accuracy of contact prediction is highly variable depending on which sequence alignment is used, and useful levels of prediction accuracy were obtained with only approximately one-third of alignments. The results presented here provide insight into the difficulties inherent in covariation analysis, and suggest that it may have limited usefulness in tertiary structure prediction. On the other hand, our ability to use covariation analysis to design stabilizing combinations of hydrophobic core substitutions attests to its potential utility for gaining deeper insight into the stability determinants and functional mechanisms of proteins with known three-dimensional structures.

Ligand interactions in the adenosine nucleotide-binding domain of the Hsp90 chaperone, GRP94. II. Ligand-mediated activation of GRP94 molecular chaperone and peptide binding activity

The N-terminal domain of eukaryotic Hsp90 proteins contains a conserved adenosine nucleotide binding pocket that also serves as the binding site for the Hsp90 inhibitors geldanamycin and radicicol. Although this domain is essential for Hsp90 function, the molecular basis for adenosine nucleotide-dependent regulation of GRP94, the endoplasmic reticulum paralog of Hsp90, remains to be established. We report that bis-ANS (1,1'-bis(4-anilino-5-napthalenesulfonic acid), an environment sensitive fluorophore known to interact with nucleotide-binding domains, binds to the adenosine nucleotide-binding domain of GRP94 and thereby activates its molecular chaperone and peptide binding activities. bis-ANS was observed to elicit a tertiary conformational change in GRP94 similar to that occurring upon heat shock, which also activates GRP94 function. bis-ANS activation of GRP94 function was efficiently blocked by radicicol, an established inhibitory ligand for the adenosine nucleotide binding pocket. Confirmation of the N-terminal nucleotide binding pocket as the bis-ANS-binding site was obtained following covalent incorporation of bis-ANS into GRP94, trypsinolysis, and sequencing of bis-ANS-labeled limit digestion products. These data identify a ligand dependent regulation of GRP94 function and suggest a model whereby GRP94 function is regulated through a ligand-dependent conversion of GRP94 from an inactive to an active conformation.

Functional inferences from reconstructed evolutionary biology involving rectified databases--an evolutionarily grounded approach to functional genomics

If bioinformatics tools are constructed to reproduce the natural, evolutionary history of the biosphere, they offer powerful approaches to some of the most difficult tasks in genomics, including the organization and retrieval of sequence data, the updating of massive genomic databases, the detection of database error, the assignment of introns, the prediction of protein conformation from protein sequences, the detection of distant homologs, the assignment of function to open reading frames, the identification of biochemical pathways from genomic data, and the construction of a comprehensive model correlating the history of biomolecules with the history of planet Earth.

The importance of ATP binding and hydrolysis by hsp90 in formation and function of protein heterocomplexes

The chaperone hsp90 is capable of binding and hydrolyzing ATP. Using information on a related ATPase, DNA gyrase B, we selected three conserved residues in hsp90's ATP-binding domain for mutation. Two of these mutations eliminate nucleotide binding, while the third retains nucleotide binding but is apparently deficient in ATP hydrolysis. We first analyzed how these mutations affect hsp90's binding to the co-chaperones p23 and Hop, and to the hydrophobic resin, phenyl-Sepharose. These experiments showed that ATP's effects, specifically, increased affinity for p23 and decreased affinity for Hop and phenyl-Sepharose, are brought on by ATP binding alone. We also tested the ability of hsp90 mutants to assist hsp70, hsp40, and Hop in the refolding of denatured firefly luciferase. While hsp90 is capable of participating in this process in a nucleotide-independent manner, the ability to hydrolyze ATP markedly potentiates hsp90's effect. Finally, we assembled progesterone receptor heterocomplexes with hsp70, hsp40, Hop, p23, and wild type or mutant hsp90. While neither ATP binding nor hydrolysis was necessary to bind hsp90 to the receptor, mature complexes containing p23 and capable of hormone binding were only obtained with wild type hsp90.

Protein structure prediction by threading. Why it works and why it does not

We developed a novel Monte Carlo threading algorithm which allows gaps and insertions both in the template structure and threaded sequence. The algorithm is able to find the optimal sequence-structure alignment and sample suboptimal alignments. Using our algorithm we performed sequence-structure alignments for a number of examples for three protein folds (ubiquitin, immunoglobulin and globin) using both "ideal" set of potentials (optimized to provide the best Z-score for a given protein) and more realistic knowledge-based potentials. Two physically different scenarios emerged. If a template structure is similar to the native one (within 2 A RMS), then (i) the optimal threading alignment is correct and robust with respect to deviations of the potential from the "ideal" one; (ii) suboptimal alignments are very similar to the optimal one; (iii) as Monte Carlo temperature decreases a sharp cooperative transition to the optimal alignment is observed. In contrast, if the template structure is only moderately close to the native structure (RMS greater than 3.5 A), then (i) the optimal alignment changes dramatically when an "ideal" potential is substituted by the real one; (ii) the structures of suboptimal alignments are very different from the optimal one, reducing the reliability of the alignment; (iii) the transition to the apparently optimal alignment is non-cooperative. In the intermediate cases when the RMS between the template and the native conformations is in the range between 2 A and 3.5 A, the success of threading alignment may depend on the quality of potentials used. These results are rationalized in terms of a threading free energy landscape. Possible ways to overcome the fundamental limitations of threading are discussed briefly.

Structure prediction in a post-genomic environment: a secondary and tertiary structural model for the initiation factor 5A family

Two predictions have been prepared for the fold of initiation factor 5A (IF5A) starting from a set of homologous sequences. In the first, a secondary structural model was predicted for the protein in 1994, when only eleven homologs (and no eubacterial homologs) had been sequenced. The second was made recently, after genome projects had generated a total of 33 sequences for the protein family from species of all three kingdoms of life. With the second set of sequences, but not with the first, it was possible to predict that the N-terminal domain of the protein folds in a possibly open beta-barrel/sandwich core structure, with a short helix capping one side of the barrel. We place the pair of predictions in the public domain before an experimental structure is known. This example illustrates the impact of genome sequencing projects on structure prediction from sequence alignments.

The 90-kDa molecular chaperone family: structure, function, and clinical applications. A comprehensive review

The 90-kDa molecular chaperone family (which comprises, among other proteins, the 90-kDa heat-shock protein, hsp90 and the 94-kDa glucose-regulated protein, grp94, major molecular chaperones of the cytosol and of the endoplasmic reticulum, respectively) has become an increasingly active subject of research in the past couple of years. These ubiquitous, well-conserved proteins account for 1-2% of all cellular proteins in most cells. However, their precise function is still far from being elucidated. Their involvement in the aetiology of several autoimmune diseases, in various infections, in recognition of malignant cells, and in antigen-presentation already demonstrates the essential role they likely will play in clinical practice of the next decade. The present review summarizes our current knowledge about the cellular functions, expression, and clinical implications of the 90-kDa molecular chaperone family and some approaches for future research.

Systems for categorizing functions of gene products

As the sequencing of the total DNA of many organisms continues, attention is turning next to the interpretation of the function of all of the genes and gene products, with the aim of learning the full meaning of the entire genetic blueprint of a sequenced organism in concrete terms. To set the stage for accomplishing this, systems need to be constructed for the expression of the functions of gene products in systematic yet rich ways. As much as possible, the same systems should be applicable to all organisms, so that when comparability exists among organisms the connections will become clear.

Fold assembly of small proteins using monte carlo simulations driven by restraints derived from multiple sequence alignments

The feasibility of predicting the global fold of small proteins by incorporating predicted secondary and tertiary restraints into ab initio folding simulations has been demonstrated on a test set comprised of 20 non-homologous proteins, of which one was a blind prediction of target 42 in the recent CASP2 contest. These proteins contain from 37 to 100 residues and represent all secondary structural classes and a representative variety of global topologies. Secondary structure restraints are provided by the PHD secondary structure prediction algorithm that incorporates multiple sequence information. Predicted tertiary restraints are derived from multiple sequence alignments via a two-step process. First, seed side-chain contacts are identified from correlated mutation analysis, and then a threading-based algorithm is used to expand the number of these seed contacts. A lattice-based reduced protein model and a folding algorithm designed to incorporate these predicted restraints is described. Depending upon fold complexity, it is possible to assemble native-like topologies whose coordinate root-mean-square deviation from native is between 3.0 A and 6.5 A. The requisite level of accuracy in side-chain contact map prediction can be roughly 25% on average, provided that about 60% of the contact predictions are correct within +/-1 residue and 95% of the predictions are correct within +/-4 residues. Precision in tertiary contact prediction is more critical than absolute accuracy. Furthermore, only a subset of the tertiary contacts, on the order of 25% of the total, is sufficient for successful topology assembly. Overall, this study suggests that the use of restraints derived from multiple sequence alignments combined with a fold assembly algorithm holds considerable promise for the prediction of the global topology of small proteins.

The amino-terminal domain of heat shock protein 90 (hsp90) that binds geldanamycin is an ATP/ADP switch domain that regulates hsp90 conformation

Many functions of the chaperone, heat shock protein 90 (hsp90), are inhibited by the drug geldanamycin that specifically binds hsp90. We have studied an amino-terminal domain of hsp90 whose crystal structure has recently been solved and determined to contain a geldanamycin-binding site. We demonstrate that, in solution, drug binding is exclusive to this domain. This domain also binds ATP linked to Sepharose through the gamma-phosphate. Binding is specific for ATP and ADP and is inhibited by geldanamycin. Mutation of four glycine residues within two proposed ATP binding motifs diminishes both geldanamycin binding and the ATP-dependent conversion of hsp90 to a conformation capable of binding the co-chaperone p23. Since p23 binding requires regions outside the 1-221 domain of hsp90, these results indicate a common site for nucleotides and geldanamycin that regulates the conformation of other hsp90 domains.