Minimal and optimal mechanisms for GroE-mediated protein folding

We have analyzed the effects of different components of the GroE chaperonin system on protein folding by using a nonpermissive substrate (i.e., one that has very low spontaneous refolding yield) for which rate data can be acquired. In the absence of GroES and nucleotides, the rate of GroEL-mediated refolding of heat- and DTT-denatured mitochondrial malate dehydrogenase was extremely low, but some three times higher than the spontaneous rate. This GroEL-mediated rate was increased 17-fold by saturating concentrations of ATP, 11-fold by ADP and GroES, and 465-fold by ATP and GroES. Optimal refolding activity was observed when the dissociation of GroES from the chaperonin complex was dramatically reduced. Although GroEL minichaperones were able to bind denatured mitochondrial malate dehydrogenase, they were ineffective in enhancing the refolding rate. The spectrum of mechanisms for GroE-mediated protein folding depends on the nature of the substrate. The minimal mechanism for permissive substrates (i.e., having significant yields of spontaneous refolding), requires only binding to the apical domain of GroEL. Slow folding rates of nonpermissive substrates are limited by the transitions between high- and low-affinity states of GroEL alone. The optimal mechanism, which requires holoGroEL, physiological amounts of GroES, and ATP hydrolysis, is necessary for the chaperonin-mediated folding of nonpermissive substrates at physiologically relevant rates under conditions in which retention of bound GroES prevents the premature release of aggregation-prone folding intermediates from the chaperonin complex. The different mechanisms are described in terms of the structural features of mini- and holo-chaperones.

Correct protein folding in glycerol

Water is the natural medium for protein folding, which is also used in all in vitro studies. In the present work, we posed, and answered affirmatively, a question of whether it is possible to fold correctly a typical protein in a nonaqueous solvent. To this end, unfolded and reduced hen egg-white lysozyme was refolded and reoxidized in glycerol containing varying amounts of water. The unfolded/reduced enzyme was found to regain spontaneously substantial catalytic activity even in the nearly anhydrous solvent; for example, the refolding yield in 99% glycerol was still some one-third of that in pure water, and one-half of that was regained even in 99.8% glycerol. The less than full recovery of the enzymatic activity in glycerol is, as in water, because of competing protein aggregation during the refolding. Lysozyme reoxidation in glycerol was successfully mediated by two dissimilar oxidizing systems, and the refolding yield was markedly affected by the pH of the last aqueous solution before the transfer into glycerol. No recovery of the lysozyme activity was observed when the refolding/reoxidation reaction was carried out in the denaturing solvent dimethyl sulfoxide. This study paves the way for a systematic investigation of the solvent effect on protein folding and demonstrates that water is not a unique milieu for this process.

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.

Structural determinants of oxidative folding in proteins

A method for determining the kinetic fate of structured disulfide species (i.e., whether they are preferentially oxidized or reshuffle back to an unstructured disulfide species) is introduced. The method relies on the sensitivity of unstructured disulfide species to low concentrations of reducing agents. Because a structured des species that preferentially reshuffles generally first rearranges to an unstructured species, a small concentration of reduced DTT (e.g., 260 μM) suffices to distinguish on-pathway intermediates from dead-end species. We apply this method to the oxidative folding of bovine pancreatic ribonuclease A (RNase A) and show that des[40–95] and des[65–72] are productive intermediates, whereas des[26–84] and des[58–110] are metastable dead-end species that preferentially reshuffle. The key factor in determining the kinetic fate of these des species is the relative accessibility of both their thiol groups and disulfide bonds. Productive intermediates tend to be disulfide-secure, meaning that their structural fluctuations preferentially expose their thiol groups, while keeping their disulfide bonds buried. By contrast, dead-end species tend to be disulfide-insecure, in that their structural fluctuations expose their disulfide bonds in concert with their thiol groups. This distinction leads to four generic types of oxidative folding pathways. We combine these results with those of earlier studies to suggest a general three-stage model of oxidative folding of RNase A and other single-domain proteins with multiple disulfide bonds.

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.

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.

Formation of a GNRA tetraloop in P5abc can disrupt an interdomain interaction in the Tetrahymena group I ribozyme

The secondary structure of a truncated P5abc subdomain (tP5abc, a 56-nucleotide RNA) of the Tetrahymena thermophila group I intron ribozyme changes when its tertiary structure forms. We have now used heteronuclear NMR spectroscopy to determine its conformation in solution. The tP5abc RNA that contains only secondary structure is extended compared with the tertiary folded form; both forms coexist in slow chemical exchange (the interconversion rate constant is slower than 1 s−1) in the presence of magnesium. Kinetic experiments have shown that tertiary folding of the P5abc subdomain is one of the earliest folding transitions in the group I intron ribozyme, and that it leads to a metastable misfolded intermediate. Previous mutagenesis studies suggest that formation of the extended P5abc structure described here destabilize a misfolded intermediate. This study shows that the P5abc RNA subdomain containing a GNRA tetraloop in P5c (in contrast to the five-nucleotide loop P5c in the tertiary folded ribozyme) can disrupt the base-paired interdomain (P14) interaction between P5c and P2.

Crystal structure of the ectodomain of Methuselah, a Drosophila G protein-coupled receptor associated with extended lifespan

The Drosophila mutant methuselah (mth) was identified from a screen for single gene mutations that extended average lifespan. Mth mutants have a 35% increase in average lifespan and increased resistance to several forms of stress, including heat, starvation, and oxidative damage. The protein affected by this mutation is related to G protein-coupled receptors of the secretin receptor family. Mth, like secretin receptor family members, has a large N-terminal ectodomain, which may constitute the ligand binding site. Here we report the 2.3-Å resolution crystal structure of the Mth extracellular region, revealing a folding topology in which three primarily β-structure-containing domains meet to form a shallow interdomain groove containing a solvent-exposed tryptophan that may represent a ligand binding site. The Mth structure is analyzed in relation to predicted Mth homologs and potential ligand binding features.

Monitoring the structure of Escherichia coli RNase P RNA in the presence of various divalent metal ions

Lead(II)-induced cleavage can be used as a tool to probe conformational changes in RNA. In this report, we have investigated the conformation of M1 RNA, the catalytic subunit of Escherichia coli RNase P, by studying the lead(II)-induced cleavage pattern in the presence of various divalent metal ions. Our data suggest that the overall conformation of M1 RNA is very similar in the presence of Mg2+, Mn2+, Ca2+, Sr2+ and Ba2+, while it is changed compared to the Mg2+-induced conformation in the presence of other divalent metal ions, Cd2+ for example. We also observed that correct folding of some M1 RNA domains is promoted by Pb2+, while folding of other domain(s) requires the additional presence of other divalent metal ions, cobalt(III) hexamine or spermidine. Based on the suppression of Pb2+ cleavage at increasing concentrations of various divalent metal ions, our findings suggest that different divalent metal ions bind with different affinities to M1 RNA as well as to an RNase P hairpin–loop substrate and yeast tRNAPhe. We suggest that this approach can be used to obtain information about the relative binding strength for different divalent metal ions to RNA in general, as well as to specific RNA divalent metal ion binding sites. Of those studied in this report, Mn2+ is generally among the strongest RNA binders.

Hydrophobic moments of protein structures: Spatially profiling the distribution

It is generally accepted that globular proteins fold with a hydrophobic core and a hydrophilic exterior. Might the spatial distribution of amino acid hydrophobicity exhibit common features? The hydrophobic profile detailing this distribution from the protein interior to exterior has been examined for 30 relatively diverse structures obtained from the Protein Data Bank, for 3 proteins of the 30S ribosomal subunit, and for a simple set of 14 decoys. A second-order hydrophobic moment has provided a simple measure of the spatial variation. Shapes of the calculated spatial profiles of all native structures have been found to be comparable. Consequently, profile shapes as well as particular profile features should assist in validating predicted protein structures and in discriminating between different protein-folding pathways. The spatial profiles of the 14 decoys are clearly distinguished from the profiles of their native structures.

Non-Arrhenius kinetics for the loop closure of a DNA hairpin

Intramolecular chain diffusion is an elementary process in the conformational fluctuations of the DNA hairpin-loop. We have studied the temperature and viscosity dependence of a model DNA hairpin-loop by FRET (fluorescence resonance energy transfer) fluctuation spectroscopy (FRETfs). Apparent thermodynamic parameters were obtained by analyzing the correlation amplitude through a two-state model and are consistent with steady-state fluorescence measurements. The kinetics of closing the loop show non-Arrhenius behavior, in agreement with theoretical prediction and other experimental measurements on peptide folding. The fluctuation rates show a fractional power dependence (β = 0.83) on the solution viscosity. A much slower intrachain diffusion coefficient in comparison to that of polypeptides was derived based on the first passage time theory of SSS [Szabo, A., Schulten, K. & Schulten, Z. (1980) J. Chem. Phys. 72, 4350–4357], suggesting that intrachain interactions, especially stacking interaction in the loop, might increase the roughness of the free energy surface of the DNA hairpin-loop.

Three-dimensional domain swapping in p13suc1 occurs in the unfolded state and is controlled by conserved proline residues

p13suc1 has two native states, a monomer and a domain-swapped dimer. We show that their folding pathways are connected by the denatured state, which introduces a kinetic barrier between monomer and dimer under native conditions. The barrier is lowered under conditions that speed up unfolding, thereby allowing, to our knowledge for the first time, a quantitative dissection of the energetics of domain swapping. The monomer–dimer equilibrium is controlled by two conserved prolines in the hinge loop that connects the exchanging domains. These two residues exploit backbone strain to specifically direct dimer formation while preventing higher-order oligomerization. Thus, the loop acts as a loaded molecular spring that releases tension in the monomer by adopting its alternative conformation in the dimer. There is an excellent correlation between domain swapping and aggregation, suggesting they share a common mechanism. These insights have allowed us to redesign the domain-swapping propensity of suc1 from a fully monomeric to a fully dimeric protein.

Evidence in support of a four transmembrane-pore-transmembrane topology model for the Arabidopsis thaliana Na+/K+ translocating AtHKT1 protein, a member of the superfamily of K+ transporters

The Arabidopsis thaliana AtHKT1 protein, a Na+/K+ transporter, is capable of mediating inward Na+ currents in Xenopus laevis oocytes and K+ uptake in Escherichia coli. HKT1 proteins are members of a superfamily of K+ transporters. These proteins have been proposed to contain eight transmembrane segments and four pore-forming regions arranged in a mode similar to that of a K+ channel tetramer. However, computer analysis of the AtHKT1 sequence identified eleven potential transmembrane segments. We have investigated the membrane topology of AtHKT1 with three different techniques. First, a gene fusion alkaline phosphatase study in E. coli clearly defined the topology of the N-terminal and middle region of AtHKT1, but the model for membrane folding of the C-terminal region had to be refined. Second, with a reticulocyte-lysate supplemented with dog-pancreas microsomes, we demonstrated that N-glycosylation occurs at position 429 of AtHKT1. An engineered unglycosylated protein variant, N429Q, mediated Na+ currents in X. laevis oocytes with the same characteristics as the wild-type protein, indicating that N-glycosylation is not essential for the functional expression and membrane targeting of AtHKT1. Five potential glycosylation sites were introduced into the N429Q. Their pattern of glycosylation supported the model based on the E. coli-alkaline phosphatase data. Third, immunocytochemical experiments with FLAG-tagged AtHKT1 in HEK293 cells revealed that the N and C termini of AtHKT1, and the regions containing residues 135–142 and 377–384, face the cytosol, whereas the region of residues 55–62 is exposed to the outside. Taken together, our results show that AtHKT1 contains eight transmembrane-spanning segments.