Crystal structures of eukaryotic translation initiation factor 5A from Methanococcus jannaschii at 1.8 Å resolution

Eukaryotic translation initiation factor 5A (eIF-5A) is a ubiquitous protein found in all eukaryotic cells. The protein is closely associated with cell proliferation in the G1–S stage of the cell cycle. Recent findings show that the eIF-5A proteins are highly expressed in tumor cells and act as a cofactor of the Rev protein in HIV-1-infected cells. The mature eIF is the only protein known to have the unusual amino acid hypusine, a post-translationally modified lysine. The crystal structure of eIF-5A from Methanococcus jannaschii (MJ eIF-5A) has been determined at 1.9 Å and 1.8 Å resolution in two crystal forms by using the multiple isomorphous replacement method and the multiwavelength anomalous diffraction method for the first crystal form and the molecular replacement method for the second crystal form. The structure consists of two folding domains, one of which is similar to the oligonucleotide-binding domain found in the prokaryotic cold shock protein and the translation initiation factor IF1 despite the absence of any significant sequence similarities. The 12 highly conserved amino acid residues found among eIF-5As include the hypusine site and form a long protruding loop at one end of the elongated molecule.

Proposed mechanism for stability of proteins to evolutionary mutations

It is shown that the sequence-ordering tendencies induced by design into different fast-folding, thermally stable native structures interfere. This interference results in a type of quasiorthogonality between optimal native structures, which divides sequence space into fast-folding, thermally stable families surrounded by slow-folding, low stability shells. A concrete example of this effect is provided by using a simple α carbon type model in which a complete correspondence is established between sequence and structure. It is speculated that gaps can occur in the space of protein-like sequences separating the sequence families and resulting in a mechanism for stability and diversity of protein sequence information.

Functional expression of the Wilson disease protein reveals mislocalization and impaired copper-dependent trafficking of the common H1069Q mutation

Wilson disease is an autosomal recessive disorder of hepatic copper metabolism caused by mutations in a gene encoding a copper-transporting P-type ATPase. To elucidate the function of the Wilson protein, wild-type and mutant Wilson cDNAs were expressed in a Menkes copper transporter-deficient mottled fibroblast cell line defective in copper export. Expression of the wild-type cDNA demonstrated trans-Golgi network localization and copper-dependent trafficking of the Wilson protein identical to previous observations for the endogenously expressed protein in hepatocytes. Furthermore, expression of the Wilson cDNA rescued the mottled phenotype as evidenced by a reduction in copper accumulation and restoration of cell viability. In contrast, expression of an H1069Q mutant Wilson cDNA did not rescue the mottled phenotype, and immunofluorescence studies showed that this mutant Wilson protein was localized in the endoplasmic reticulum. Consistent with these findings, pulse–chase analysis demonstrated a 5-fold decrease in the half-life of the H1069Q mutant as compared with the wild-type protein. Maintenance of these transfected cell lines at 28°C resulted in localization of the H1069Q protein in the trans-Golgi network, suggesting that a temperature-sensitive defect in protein folding followed by degradation constitutes the molecular basis of Wilson disease in patients harboring the H1069Q mutation. Taken together, these studies describe a tractable expression system for elucidating the function and localization of the copper-transporting ATPases in mammalian cells and provide compelling evidence that the Wilson protein can functionally substitute for the Menkes protein, supporting the concept that these proteins use common biochemical mechanisms to effect cellular copper homeostasis.

Mutations in the DnaK chaperone affecting interaction with the DnaJ cochaperone

Hsp70 chaperones assist protein folding by ATP-controlled cycles of substrate binding and release. ATP hydrolysis is the rate-limiting step of the ATPase cycle that causes locking in of substrates into the substrate-binding cavity of Hsp70. This key step is strongly stimulated by DnaJ cochaperones. We show for the Escherichia coli Hsp70 homolog, DnaK, that stimulation by DnaJ requires the linked ATPase and substrate-binding domains of DnaK. Functional interaction with DnaJ is affected by mutations in an exposed channel located in the ATPase domain of DnaK. It is proposed that binding to this channel, possibly involving the J-domain, allows DnaJ to couple substrate binding with ATP hydrolysis by DnaK. Evolutionary conservation of the channel and the J-domain suggests conservation of the mechanism of action of DnaJ proteins.

Mode matches and their locations in the hydrophobic free energy sequences of peptide ligands and their receptor eigenfunctions

Patterns in sequences of amino acid hydrophobic free energies predict secondary structures in proteins. In protein folding, matches in hydrophobic free energy statistical wavelengths appear to contribute to selective aggregation of secondary structures in “hydrophobic zippers.” In a similar setting, the use of Fourier analysis to characterize the dominant statistical wavelengths of peptide ligands’ and receptor proteins’ hydrophobic modes to predict such matches has been limited by the aliasing and end effects of short peptide lengths, as well as the broad-band, mode multiplicity of many of their frequency (power) spectra. In addition, the sequence locations of the matching modes are lost in this transformation. We make new use of three techniques to address these difficulties: (i) eigenfunction construction from the linear decomposition of the lagged covariance matrices of the ligands and receptors as hydrophobic free energy sequences; (ii) maximum entropy, complex poles power spectra, which select the dominant modes of the hydrophobic free energy sequences or their eigenfunctions; and (iii) discrete, best bases, trigonometric wavelet transformations, which confirm the dominant spectral frequencies of the eigenfunctions and locate them as (absolute valued) moduli in the peptide or receptor sequence. The leading eigenfunction of the covariance matrix of a transmembrane receptor sequence locates the same transmembrane segments seen in n-block-averaged hydropathy plots while leaving the remaining hydrophobic modes unsmoothed and available for further analyses as secondary eigenfunctions. In these receptor eigenfunctions, we find a set of statistical wavelength matches between peptide ligands and their G-protein and tyrosine kinase coupled receptors, ranging across examples from 13.10 amino acids in acid fibroblast growth factor to 2.18 residues in corticotropin releasing factor. We find that the wavelet-located receptor modes in the extracellular loops are compatible with studies of receptor chimeric exchanges and point mutations. A nonbinding corticotropin-releasing factor receptor mutant is shown to have lost the signatory mode common to the normal receptor and its ligand. Hydrophobic free energy eigenfunctions and their transformations offer new quantitative physical homologies in database searches for peptide-receptor matches.

Rapid evolution of the DNA-binding site in LAGLIDADG homing endonucleases

Sequence analysis of chloroplast and mitochondrial large subunit rRNA genes from over 75 green algae disclosed 28 new group I intron-encoded proteins carrying a single LAGLIDADG motif. These putative homing endonucleases form four subfamilies of homologous enzymes, with the members of each subfamily being encoded by introns sharing the same insertion site. We showed that four divergent endonucleases from the I-CreI subfamily cleave the same DNA substrates. Mapping of the 66 amino acids that are conserved among the members of this subfamily on the 3-dimensional structure of I-CreI bound to its recognition sequence revealed that these residues participate in protein folding, homodimerization, DNA recognition and catalysis. Surprisingly, only seven of the 21 I-CreI amino acids interacting with DNA are conserved, suggesting that I-CreI and its homologs use different subsets of residues to recognize the same DNA sequence. Our sequence comparison of all 45 single-LAGLIDADG proteins identified so far suggests that these proteins share related structures and that there is a weak pressure in each subfamily to maintain identical protein–DNA contacts. The high sequence variability we observed in the DNA-binding site of homologous LAGLIDADG endonucleases provides insight into how these proteins evolve new DNA specificity.

The transcriptional program of a human B cell line in response to Myc

The proto-oncogene c-myc (myc) encodes a transcription factor (Myc) that promotes growth, proliferation and apoptosis. Myc has been suggested to induce these effects by induction/repression of downstream genes. Here we report the identification of potential Myc target genes in a human B cell line that grows and proliferates depending on conditional myc expression. Oligonucleotide microarrays were applied to identify downstream genes of Myc at the level of cytoplasmic mRNA. In addition, we identified potential Myc target genes in nuclear run-on experiments by changes in their transcription rate. The identified genes belong to gene classes whose products are involved in amino acid/protein synthesis, lipid metabolism, protein turnover/folding, nucleotide/DNA synthesis, transport, nucleolus function/RNA binding, transcription and splicing, oxidative stress and signal transduction. The identified targets support our current view that myc acts as a master gene for growth control and increases transcription of a large variety of genes.

CKAAPs DB: a conserved key amino acid positions database

The Conserved Key Amino Acid Positions DataBase (CKAAPs DB) provides access to an analysis of structurally similar proteins with dissimilar sequences where key residues within a common fold are identified. The derivation and significance of CKAAPs starting from pairwise structure alignments is described fully in Reddy et al. [Reddy,B.V.B., Li,W.W., Shindyalov,I.N. and Bourne,P.E. (2000) Proteins, in press]. The CKAAPs identified from this theoretical analysis are provided to experimentalists and theoreticians for potential use in protein engineering and modeling. It has been suggested that CKAAPs may be crucial features for protein folding, structural stability and function. Over 170 substructures, as defined by the Combinatorial Extension (CE) database, which are found in approximately 3000 representative polypeptide chains have been analyzed and are available in the CKAAPs DB. CKAAPs DB also provides CKAAPs of the representative set of proteins derived from the CE and FSSP databases. Thus the database contains over 5000 representative poly­peptide chains, covering all known structures in the PDB. A web interface to a relational database permits fast retrieval of structure-sequence alignments, CKAAPs and associated statistics. Users may query by PDB ID, protein name, function and Enzyme Classification number. Users may also submit protein alignments of their own to obtain CKAAPs. An interface to display CKAAPs on each structure from a web browser is also being implemented. CKAAPs DB is maintained by the San Diego Supercomputer Center and accessible at the URL http://ckaaps.sdsc.edu.

Pseudoknots in prion protein mRNAs confirmed by comparative sequence analysis and pattern searching

The human prion gene contains five copies of a 24 nt repeat that is highly conserved among species. An analysis of folding free energies of the human prion mRNA, in particular in the repeat region, suggested biased codon selection and the presence of RNA patterns. In particular, pseudoknots, similar to the one predicted by Wills in the human prion mRNA, were identified in the repeat region of all available prion mRNAs available in GenBank, but not those of birds and the red slider turtle. An alignment of these mRNAs, which share low sequence homology, shows several co-variations that maintain the pseudoknot pattern. The presence of pseudoknots in yeast Sup35p and Rnq1 suggests acquisition in the prokaryotic era. Computer generated three-dimensional structures of the human prion pseudoknot highlight protein and RNA interaction domains, which suggest a possible effect in prion protein translation. The role of pseudoknots in prion diseases is discussed as individuals with extra copies of the 24 nt repeat develop the familial form of Creutzfeldt–Jakob disease.

Ribosomal protein S7 from Escherichia coli uses the same determinants to bind 16S ribosomal RNA and its messenger RNA

Ribosomal protein S7 from Escherichia coli binds to the lower half of the 3′ major domain of 16S rRNA and initiates its folding. It also binds to its own mRNA, the str mRNA, and represses its translation. Using filter binding assays, we show in this study that the same mutations that interfere with S7 binding to 16S rRNA also weaken its affinity for its mRNA. This suggests that the same protein regions are responsible for mRNA and rRNA binding affinities, and that S7 recognizes identical sequence elements within the two RNA targets, although they have dissimilar secondary structures. Overexpression of S7 is known to inhibit bacterial growth. This phenotypic growth defect was relieved in cells overexpressing S7 mutants that bind poorly the str mRNA, confirming that growth impairment is controlled by the binding of S7 to its mRNA. Interestingly, a mutant with a short deletion at the C-terminus of S7 was more detrimental to cell growth than wild-type S7. This suggests that the C-terminal portion of S7 plays an important role in ribosome function, which is perturbed by the deletion.

Telomerase inhibitors based on quadruplex ligands selected by a fluorescence assay

The reactivation of telomerase activity in most cancer cells supports the concept that telomerase is a relevant target in oncology, and telomerase inhibitors have been proposed as new potential anticancer agents. The telomeric G-rich single-stranded DNA can adopt in vitro an intramolecular quadruplex structure, which has been shown to inhibit telomerase activity. We used a fluorescence assay to identify molecules that stabilize G-quadruplexes. Intramolecular folding of an oligonucleotide with four repeats of the human telomeric sequence into a G-quadruplex structure led to fluorescence excitation energy transfer between a donor (fluorescein) and an acceptor (tetramethylrhodamine) covalently attached to the 5′ and 3′ ends of the oligonucleotide, respectively. The melting of the G-quadruplex was monitored in the presence of putative G-quadruplex-binding molecules by measuring the fluorescence emission of the donor. A series of compounds (pentacyclic crescent-shaped dibenzophenanthroline derivatives) was shown to increase the melting temperature of the G-quadruplex by 2–20°C at 1 μM dye concentration. This increase in Tm value was well correlated with an increase in the efficiency of telomerase inhibition in vitro. The best telomerase inhibitor showed an IC50 value of 28 nM in a standard telomerase repeat amplification protocol assay. Fluorescence energy transfer can thus be used to reveal the formation of four-stranded DNA structures, and its stabilization by quadruplex-binding agents, in an effort to discover new potent telomerase inhibitors.

Synechocystis HSP17 is an amphitropic protein that stabilizes heat-stressed membranes and binds denatured proteins for subsequent chaperone-mediated refolding

The small heat shock proteins (sHSPs) are ubiquitous stress proteins proposed to act as molecular chaperones to prevent irreversible protein denaturation. We characterized the chaperone activity of Synechocystis HSP17 and found that it has not only protein-protective activity, but also a previously unrecognized ability to stabilize lipid membranes. Like other sHSPs, recombinant Synechocystis HSP17 formed stable complexes with denatured malate dehydrogenase and served as a reservoir for the unfolded substrate, transferring it to the DnaK/DnaJ/GrpE and GroEL/ES chaperone network for subsequent refolding. Large unilamellar vesicles made of synthetic and cyanobacterial lipids were found to modulate this refolding process. Investigation of HSP17-lipid interactions revealed a preference for the liquid crystalline phase and resulted in an elevated physical order in model lipid membranes. Direct evidence for the participation of HSP17 in the control of thylakoid membrane physical state in vivo was gained by examining an hsp17− deletion mutant compared with the isogenic wild-type hsp17+ revertant Synechocystis cells. We suggest that, together with GroEL, HSP17 behaves as an amphitropic protein and plays a dual role. Depending on its membrane or cytosolic location, it may function as a “membrane stabilizing factor” as well as a member of a multichaperone protein-folding network. Membrane association of sHSPs could antagonize the heat-induced hyperfluidization of specific membrane domains and thereby serve to preserve structural and functional integrity of biomembranes.

Point mutations and sequence variability in proteins: Redistributions of preexisting populations

Here we study the effect of point mutations in proteins on the redistributions of the conformational substates. We show that regardless of the location of a mutation in the protein structure and of its type, the observed movements of the backbone recur largely at the same positions in the structures. Despite the different interactions that are disrupted and formed by the residue substitution, not only are the conformations very similar, but the regions that move are also the same, regardless of their sequential or spatial distance from the mutation. This observation leads us to conclude that, apart from some extreme cases, the details of the interactions are not critically important in determining the protein conformation or in specifying which parts of the protein would be more prone to take on different local conformations in response to changes in the sequence. This finding further illustrates why proteins manifest a robustness toward many mutational events. This nonuniform distribution of the conformer population is consistently observed in a variety of protein structural types. Topology is critically important in determining folding pathways, kinetics, building block cutting, and anatomy trees. Here we show that topology is also very important in determining which regions of the protein structure will respond to sequence changes, regardless of the sequential or spatial location of the mutation.

Ligand-independent assembly of recombinant human CD1 by using oxidative refolding chromatography

CD1 is an MHC class I-like antigen-presenting molecule consisting of a heavy chain and β2-microglobulin light chain. The in vitro refolding of synthetic MHC class I molecules has always required the presence of ligand. We report here the use of a folding method using an immobilized chaperone fragment, a protein disulphide isomerase, and a peptidyl-prolyl cis-trans isomerase (oxidative refolding chromatography) for the fast and efficient assembly of ligand-free and ligand-associated CD1a and CD1b, starting with material synthesized in Escherichia coli. The results suggest that “empty” MHC class I-like molecules can assemble and remain stable at physiological temperatures in the absence of ligand. The use of oxidative refolding chromatography thus is extended to encompass complex multisubunit proteins and specifically to members of the extensive, functionally diverse and important immunoglobulin supergene family of proteins, including those for which a ligand has yet to be identified.

Requirement of the Lec35 Gene for All Known Classes of Monosaccharide-P-Dolichol-dependent Glycosyltransferase Reactions in Mammals

The Lec35 gene product (Lec35p) is required for utilization of the mannose donor mannose-P-dolichol (MPD) in synthesis of both lipid-linked oligosaccharides (LLOs) and glycosylphosphatidylinositols, which are important for functions such as protein folding and membrane anchoring, respectively. The hamster Lec35 gene is shown to encode the previously identified cDNA SL15, which corrects the Lec35 mutant phenotype and predicts a novel endoplasmic reticulum membrane protein. The mutant hamster alleles Lec35.1 and Lec35.2 are characterized, and the human Lec35 gene (mannose-P-dolichol utilization defect 1) was mapped to 17p12-13. To determine whether Lec35p was required only for MPD-dependent mannosylation of LLO and glycosylphosphatidylinositol intermediates, two additional lipid-mediated reactions were investigated: MPD-dependent C-mannosylation of tryptophanyl residues, and glucose-P-dolichol (GPD)-dependent glucosylation of LLO. Both were found to require Lec35p. In addition, the SL15-encoded protein was selective for MPD compared with GPD, suggesting that an additional GPD-selective Lec35 gene product remains to be identified. The predicted amino acid sequence of Lec35p does not suggest an obvious function or mechanism. By testing the water-soluble MPD analog mannose-β-1-P-citronellol in an in vitro system in which the MPD utilization defect was preserved by permeabilization with streptolysin-O, it was determined that Lec35p is not directly required for the enzymatic transfer of mannose from the donor to the acceptor substrate. These results show that Lec35p has an essential role for all known classes of monosaccharide-P-dolichol-dependent reactions in mammals. The in vitro data suggest that Lec35p controls an aspect of MPD orientation in the endoplasmic reticulum membrane that is crucial for its activity as a donor substrate.