Physical reasons for the unusual α-helix stabilization afforded by charged or neutral polar residues in alanine-rich peptides

We have carried out conformational energy calculations on alanine-based copolymers with the sequence Ac-AAAAAXAAAA-NH2 in water, where X stands for lysine or glutamine, to identify the underlying source of stability of alanine-based polypeptides containing charged or highly soluble polar residues in the absence of charge–charge interactions. The results indicate that ionizable or neutral polar residues introduced into the sequence to make them soluble sequester the water away from the CO and NH groups of the backbone, thereby enabling them to form internal hydrogen bonds. This solvation effect dictates the conformational preference and, hence, modifies the conformational propensity of alanine residues. Even though we carried out simulations for specific amino acid sequences, our results provide an understanding of some of the basic principles that govern the process of folding of these short sequences independently of the kind of residues introduced to make them soluble. In addition, we have investigated through simulations the effect of the bulk dielectric constant on the conformational preferences of these peptides. Extensive conformational Monte Carlo searches on terminally blocked 10-mer and 16-mer homopolymers of alanine in the absence of salt were carried out assuming values for the dielectric constant of the solvent ɛ of 80, 40, and 2. Our simulations show a clear tendency of these oligopeptides to augment the α-helix content as the bulk dielectric constant of the solvent is lowered. This behavior is due mainly to a loss of exposure of the CO and NH groups to the aqueous solvent. Experimental evidence indicates that the helical propensity of the amino acids in water shows a dramatic increase on addition of certain alcohols, such us trifluoroethanol. Our results provide a possible explanation of the mechanism by which alcohol/water mixtures affect the free energy of helical alanine oligopeptides relative to nonhelical ones.

An amino acid as a cofactor for a catalytic polynucleotide

Natural ribozymes require metal ion cofactors that aid both in structural folding and in chemical catalysis. In contrast, many protein enzymes produce dramatic rate enhancements using only the chemical groups that are supplied by their constituent amino acids. This fact is widely viewed as the most important feature that makes protein a superior polymer for the construction of biological catalysts. Herein we report the in vitro selection of a catalytic DNA that uses histidine as an active component for an RNA cleavage reaction. An optimized deoxyribozyme from this selection requires l-histidine or a closely related analog to catalyze RNA phosphoester cleavage, producing a rate enhancement of ≈1-million-fold over the rate of substrate cleavage in the absence of enzyme. Kinetic analysis indicates that a DNA–histidine complex may perform a reaction that is analogous to the first step of the proposed catalytic mechanism of RNase A, in which the imidazole group of histidine serves as a general base catalyst. Similarly, ribozymes of the “RNA world” may have used amino acids and other small organic cofactors to expand their otherwise limited catalytic potential.

Role of the J-domain in the cooperation of Hsp40 with Hsp70

The Escherichia coli Hsp40 DnaJ and Hsp70 DnaK cooperate in the binding of proteins at intermediate stages of folding, assembly, and translocation across membranes. Binding of protein substrates to the DnaK C-terminal domain is controlled by ATP binding and hydrolysis in the N-terminal ATPase domain. The interaction of DnaJ with DnaK is mediated at least in part by the highly conserved N-terminal J-domain of DnaJ that includes residues 2–75. Heteronuclear NMR experiments with uniformly 15N-enriched DnaJ2–75 indicate that the chemical environment of residues located in helix II and the flanking loops is perturbed on interaction with DnaK or a truncated DnaK molecule, DnaK2–388. NMR signals corresponding to these residues broaden and exhibit changes in chemical shifts in the presence of DnaK(MgADP). Addition of MgATP largely reversed the broadening, indicating that NMR signals of DnaJ2–75 respond to ATP-dependent changes in DnaK. The J-domain interaction is localized to the ATPase domain of DnaK and is likely to be dominated by electrostatic interactions. The results suggest that the J-domain tethers DnaK to DnaJ-bound substrates, which DnaK then binds with its C-terminal peptide-binding domain.

Dynamical principles in biological processes: A model of charge migration in proteins and DNA

The generalized master equations (GMEs) that contain multiple time scales have been derived quantum mechanically. The GME method has then been applied to a model of charge migration in proteins that invokes the hole hopping between local amino acid sites driven by the torsional motions of the floppy backbones. This model is then applied to analyze the experimental results for sequence-dependent long-range hole transport in DNA reported by Meggers et al. [Meggers, E., Michel-Beyerle, M. E., & Giese, B. (1998) J. Am. Chem. Soc. 120, 12950–12955]. The model has also been applied to analyze the experimental results of femtosecond dynamics of DNA-mediated electron transfer reported by Zewail and co-workers [Wan, C., Fiebig, T., Kelley, S. O., Treadway, C. R., Barton, J. K. & Zewail, A. H. (1999) Proc. Natl. Acad. Sci. USA 96, 6014–6019]. The initial events in the dynamics of protein folding have begun to attract attention. The GME obtained in this paper will be applicable to this problem.

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.