A quantitative, high-throughput screen for protein stability

In proteomic research, it is often necessary to screen a large number of polypeptides for the presence of stable structure. Described here is a technique (referred to as SUPREX, stability of unpurified proteins from rates of H/D exchange) for measuring the stability of proteins in a rapid, high-throughput fashion. The method uses hydrogen exchange to estimate the stability of microgram quantities of unpurified protein extracts by using matrix-assisted laser desorption/ionization MS. The stabilities of maltose binding protein and monomeric λ repressor variants determined by SUPREX agree well with stability data obtained from conventional CD denaturation of purified protein. The method also can detect the change in stability caused by the binding of maltose to maltose binding protein. The results demonstrate the precision of the method over a wide range of stabilities.

Site-specific conformational determination in thermal unfolding studies of helical peptides using vibrational circular dichroism with isotopic substitution

Understanding the detailed mechanism of protein folding requires dynamic, site-specific stereochemical information. The short time response of vibrational spectroscopies allows evaluation of the distribution of populations in rapid equilibrium as the peptide unfolds. Spectral shifts associated with isotopic labels along with local stereochemical sensitivity of vibrational circular dichroism (VCD) allow determination of the segment sequence of unfolding. For a series of alanine-rich peptides that form α-helices in aqueous solution, we used isotopic labeling and VCD to demonstrate that the α-helix noncooperatively unwinds from the ends with increasing temperature. For these blocked peptides, the C-terminal is frayed at 5°C. Ab initio level theoretical simulations of the IR and VCD band shapes are used to analyze the spectra and to confirm the conformation of the labeled components. The VCD signals associated with the labeled residues are amplified by coupling to the nonlabeled parts of the molecule. Thus small labeled segments are detectable and stereochemically defined in moderately large peptides in this report of site-specific peptide VCD conformational analysis.

A meanfield approach to the thermodynamics of a protein-solvent system with application to the oligomerization of the tumor suppressor p53

The thermodynamic stability and oligomerization status of the tumor suppressor p53 tetramerization domain have been studied experimentally and theoretically. A series of hydrophilic mutations at Met-340 and Leu-344 of human p53 were designed to disrupt the hydrophobic dimer–dimer interface of the tetrameric oligomerization domain of p53 (residues 325–355). Meanfield calculations of the free energy of the solvated mutants as a function of interdimer distance were compared with experimental data on the thermal stability and oligomeric state (tetramer, dimer, or equilibrium mixture of both) of each mutant. The calculations predicted a decreasing stability and oligomeric state for the following amino acids at residue 340: Met (tetramer) > Ser Asp, His, Gln, > Glu, Lys (dimer), whereas the experimental results showed the following order: Met (tetramer) > Ser > Gln > His, Lys > Asp, Glu (dimers). For residue 344, the calculated trend was Leu (tetramer) > Ala > Arg, Gln, Lys (dimer), and the experimental trend was Leu (tetramer) > Ala, Arg, Gln, Lys (dimer). The discrepancy for the lysine side chain at residue 340 is attributed to the dual nature of lysine, both hydrophobic and charged. The incorrect prediction of stability of the mutant with Asp at residue 340 is attributed to the fact that within the meanfield approach, we use the wild-type backbone configuration for all mutants, but low melting temperatures suggest a softening of the α-helices at the dimer–dimer interface. Overall, this initial application of meanfield theory toward a protein-solvent system is encouraging for the application of the theoretical model to more complex systems.

Unbinding forces of single antibody-antigen complexes correlate with their thermal dissociation rates

Point mutants of three unrelated antifluorescein antibodies were constructed to obtain nine different single-chain Fv fragments, whose on-rates, off-rates, and equilibrium binding affinities were determined in solution. Additionally, activation energies for unbinding were estimated from the temperature dependence of the off-rate in solution. Loading rate-dependent unbinding forces were determined for single molecules by atomic force microscopy, which extrapolated at zero force to a value close to the off-rate measured in solution, without any indication for multiple transition states. The measured unbinding forces of all nine mutants correlated well with the off-rate in solution, but not with the temperature dependence of the reaction, indicating that the same transition state must be crossed in spontaneous and forced unbinding and that the unbinding path under load cannot be too different from the one at zero force. The distance of the transition state from the ground state along the unbinding pathway is directly proportional to the barrier height, regardless of the details of the binding site, which most likely reflects the elasticity of the protein in the unbinding process. Atomic force microscopy thus can be a valuable tool for the characterization of solution properties of protein-ligand systems at the single molecule level, predicting relative off-rates, potentially of great value for combinatorial chemistry and biology.

Functional and physiological consequences of genetic variation at phosphoglucose isomerase: Heat shock protein expression is related to enzyme genotype in a montane beetle

Allele frequency variation at the phosphoglucose isomerase (PGI) locus in Californian populations of the beetle Chrysomela aeneicollis suggests that PGI may be undergoing natural selection. We quantified (i) apparent Michaelis-Menten constant (Km) of fructose 6-phosphate at different temperatures and (ii) thermal stability for three common PGI genotypes (1–1, 1–4, and 4–4). We also measured air temperature (Ta) and beetle body temperature (Tb) in three montane drainages in the Sierra Nevada, California. Finally, we measured 70-kDa heat shock protein (Hsp70) expression in field-collected and laboratory-acclimated beetles. We found that PGI allele 1 predominated in the northernmost drainage, Rock Creek (RC), which was also significantly cooler than the southernmost drainage, Big Pine Creek (BPC), where PGI allele 4 predominated. Allele frequencies and air temperatures were intermediate in the middle drainage, Bishop Creek (BC). Differences among genotypes in Km (1–1 > 1–4 > 4–4) and thermal stability (4–4 > 1–4 > 1–1) followed a pattern consistent with temperature adaptation. In nature, Tb was closely related to Ta. Hsp70 expression in adult beetles decreased with elevation and differed among drainages (BPC > BC > RC). After laboratory acclimation (8 days, 20°C day, 4°C night) and heat shock (4 h, 28–36°C), Hsp70 expression was greater for RC than BPC beetles. In RC, field-collected beetles homozygous for PGI 1–1 had higher Hsp70 levels than heterozygotes or a 4–4 homozygote. These results reveal functional and physiological differences among PGI genotypes, which suggest that montane populations of this beetle are locally adapted to temperature.

The domain structure of protoporphyrinogen oxidase, the molecular target of diphenyl ether-type herbicides

Protoporphyrinogen oxidase (EC 1–3-3–4), the 60-kDa membrane-bound flavoenzyme that catalyzes the final reaction of the common branch of the heme and chlorophyll biosynthesis pathways in plants, is the molecular target of diphenyl ether-type herbicides. It is highly resistant to proteases (trypsin, endoproteinase Glu-C, or carboxypeptidases A, B, and Y), because the protein is folded into an extremely compact form. Trypsin maps of the native purified and membrane-bound yeast protoporphyrinogen oxidase show that this basic enzyme (pI > 8.5) was cleaved at a single site under nondenaturing conditions, generating two peptides with relative molecular masses of 30,000 and 35,000. The endoproteinase Glu-C also cleaved the protein into two peptides with similar masses, and there was no additional cleavage site under mild denaturing conditions. N-terminal peptide sequence analysis of the proteolytic (trypsin and endoproteinase Glu-C) peptides showed that both cleavage sites were located in putative connecting loop between the N-terminal domain (25 kDa) with the βαβ ADP-binding fold and the C-terminal domain (35 kDa), which possibly is involved in the binding of the isoalloxazine moiety of the FAD cofactor. The peptides remained strongly associated and fully active with the Km for protoporphyrinogen and the Ki for various inhibitors, diphenyl-ethers, or diphenyleneiodonium derivatives, identical to those measured for the native enzyme. However, the enzyme activity of the peptides was much more susceptible to thermal denaturation than that of the native protein. Only the C-terminal domain of protoporphyrinogen oxidase was labeled specifically in active site-directed photoaffinity-labeling experiments. Trypsin may have caused intramolecular transfer of the labeled group to reactive components of the N-terminal domain, resulting in nonspecific labeling. We suggest that the active site of protoporphyrinogen oxidase is in the C-terminal domain of the protein, at the interface between the C- and N-terminal domains.

Interaction of the histone (H3-H4)2 tetramer of the nucleosome with positively supercoiled DNA minicircles: Potential flipping of the protein from a left- to a right-handed superhelical form.

We have studied the ability of the histone (H3-H4)2 tetramer, the central part of the nucleosome of eukaryotic chromatin, to form particles on DNA minicircles of negative and positive superhelicities, and the effect of relaxing these particles with topoisomerase I. The results show that even modest positive torsional stress from the DNA, and in particular that generated by DNA thermal fluctuations, can trigger a major, reversible change in the conformation of the particle. Neither a large excess of naked DNA, nor a crosslink between the two H3s prevented the transition from one form to the other. This suggested that during the transition, the histones neither dissociated from the DNA nor were even significantly reshuffled. Moreover, the particles reconstituted on negatively and positively supercoiled minicircles look similar under electron microscopy. These data agree best with a transition involving a switch of the wrapped DNA from a left- to a right-handed superhelix. It is further proposed, based on the left-handed overall superhelical conformation of the tetramer within the octamer [Arents, G., Burlingame, R. W., Wang, B. C., Love, W. E. & Moudrianakis, E. N. (1991) Proc. Natl.Acad. Sci. USA 88, 10148-10152] that this change in DNA topology is mediated by a similar change in the topology of the tetramer itself, which may occur through a rotation (or a localized deformation) of the two H3-H4 dimers about their H3-H3 interface. Potential implications of this model for nucleosome dynamics in vivo are discussed.

Thermal unfolding of human high-density apolipoprotein A-1: implications for a lipid-free molten globular state.

Apolipoprotein A-1 (apoA-1) in complex with high-density lipoprotein is critically involved in the transport and metabolism of cholesterol and in the pathogenesis of atherosclerosis. We reexamined the thermal unfolding of lipid-free apoA-1 in low-salt solution at pH approximately 7, by using differential scanning calorimetry and circular dichroism. At protein concentrations <5 mg/ml, thermal unfolding of apoA-1 is resolved as an extended peak (25 degrees C-90 degrees C) that can be largely accounted for by a single reversible non-two-state transition with midpoint Tm 57 +/- 1 degree C, calorimetric enthalpy deltaH(Tm)= 200 +/- 20 kcal/mol (1 kcal = 4.18 kJ), van't Hoff enthalpy deltaHv(Tm) approximately 32.5 kcal/mol, and cooperativity deltaHv(Tm)/deltaH(Tm) approximately 0.16. The enthalpy deltaH(Tm) can be accounted for by melting of the alpha-helical structure that is inferred by CD to constitute approximately 60% of apoA-1 amino acids. Farand near-UV CD spectra reveal noncoincident melting of the secondary and tertiary structural elements and indicate a well-defined secondary structure but a largely melted tertiary structure for apoA-1 at approximately 37 degrees C and pH 7. This suggests a molten globular-like state for lipid-free apoA-1 under near-physiological conditions. Our results suggest that in vivo lipid binding by apoA-1 may be mediated via the molten globular apolipoprotein state in plasma.

Molecular mechanism of protein-retinal coupling in bacteriorhodopsin.

Bacteriorhodopsin is a membrane protein that functions as a light-driven proton pump. Each cycle of proton transport is initiated by the light-induced isomerization of retinal from the all-trans to 13-cis configuration and is completed by the protein-driven reisomerization of retinal to the all-trans configuration. Previous studies have shown that replacement of Leu-93, a residue in close proximity to the 13-methyl group of retinal, by alanine, resulted in a 250-fold increase in the time required to complete each photocycle. Here, we show that the kinetic defect in the photocycle of the Leu-93-->Ala mutant occurs at a stage after the completion of proton transport and can be overcome in the presence of strong background illumination. Time-resolved retinal-extraction experiments demonstrate the continued presence of a 13-cis intermediate in the photocycle of the Leu-93-->Ala mutant well after the completion of proton release and uptake. These results indicate that retinal reisomerization is kinetically the rate-limiting step in the photocycle of this mutant and that the slow thermal reisomerization can be bypassed by the absorption of a second photon. The effects observed for the Leu-93-->Ala mutant are not observed upon replacement of any other residue in van der Waals contact with retinal or upon replacement of Leu-93 by valine. We conclude that the contact between Leu-93 and the 13-methyl group of retinal plays a key role in controlling the rate of protein conformational changes associated with retinal reisomerization and return of the protein to the initial state.

Tertiary structure of an amyloid immunoglobulin light chain protein: a proposed model for amyloid fibril formation.

An immunoglobulin light chain protein was isolated from the urine of an individual (BRE) with systemic amyloidosis. Complete amino acid sequence of the variable region of the light chain (VL) protein established it as a kappa I, which when compared with other kappa I amyloid associated proteins had unique residues, including Ile-34, Leu-40, and Tyr-71. To study the tertiary structure, BRE VL was expressed in Escherichia coli by using a PCR product amplified from the patient BRE's bone marrow DNA. The PCR product was ligated into pCZ11, a thermal-inducible replication vector. Recombinant BRE VL was isolated, purified to homogeneity, and crystallized by using ammonium sulfate as the precipitant. Two crystal forms were obtained. In crystal form I the BRE VL kappa domain crystallizes as a dimer with unit cell constants isomorphous to previously published kappa protein structures. Comparison with a nonamyloid VL kappa domain from patient REI, identified significant differences in position of residues in the hypervariable segments plus variations in framework region (FR) segments 40-46 (FR2) and 66-67 (FR3). In addition, positional differences can be seen along the two types of local diads, corresponding to the monomer-monomer and dimer-dimer interfaces. From the packing diagram, a model for the amyloid light chain (AL) fibril is proposed based on a pseudohexagonal spiral structure with a rise of approximately the width of two dimers per 360 degree turn. This spiral structure could be consistent with the dimensions of amyloid fibrils as determined by electron microscopy.

Stringent chemical and thermal regulation of recombinant gene expression by vaccinia virus vectors in mammalian cells.

We developed a stringently regulated expression system for mammalian cells that uses (i) the RNA polymerase, phi 10 promoter, and T phi transcriptional terminator of bacteriophage T7; (ii) the lac repressor, lac operator, rho-independent transcriptional terminators and the gpt gene of Escherichia coli; (iii) the RNA translational enhancer of encephalomyocarditis virus; and (iv) the genetic background of vaccinia virus. In cells infected with the recombinant vaccinia virus, reporter beta-galactosidase synthesis was not detected in the absence of inducer. An induction of at least 10,000- to 20,000-fold occurred upon addition of isopropyl beta-D-thiogalactopyranoside or by temperature elevation from 30 to 37 degrees C using a temperature-sensitive lac repressor. Regulated synthesis of the secreted and highly glycosylated human immunodeficiency virus 1 envelope protein gp120 was also demonstrated. Yields of both proteins were approximately 2 mg per 10(8) cells in 24 hr. Plasmid transfer vectors for cloning and expression of complete or incomplete open reading frames in recombinant vaccinia viruses are described.

Study of zinc protein ligands in a halophilic enzyme

Glucose dehydrogenase (EC 1.1.1.47) from the halophilic Archaeon Haloferax mediterranei belongs to the medium-chain alcohol dehydrogenase superfamily and requires a zinc ion for catalysis. The zinc ion is coordinated by a histidine, a water molecule and two other ligands from the protein or the substrate, which vary during the catalytic cycle of the enzyme. In many enzymes of this superfamily one of the zinc ligands is commonly cysteine, which is replaced by an aspartate residue at position 38 in the halophilic enzyme. This change has been only observed in glucose dehydrogenases from extremely halophilic microorganisms belonging to the Archaea Domain. This paper describes biochemical studies and structural comparisons to analyze the role of sequence differences between thermophilic and halophilic glucose dehydrogenases which contain a zinc ion within the protein surrounded by three ligands. Whilst the catalytic activity of the D38C GlcDH mutant is reduced, its thermal stability is enhanced, consistent with the greater structural similarity between this mutant and the homologous thermophilic enzyme from Thermoplasma acidophilum.

Mechanisms of Thermal Adaptation Revealed From the Genomes of the Antarctic Archaea Methanogenium frigidum and Methanococcoides burtonii

We generated draft genome sequences for two cold-adapted Archaea, Methanogenium frigidum and Methanococcoides burtonii, to identify genotypic characteristics that distinguish them from Archaea with a higher optimal growth temperature (OGT). Comparative genomics revealed trends in amino acid and tRNA composition, and structural features of proteins. Proteins from the cold-adapted Archaea are characterized by a higher content of noncharged polar amino acids, particularly Gin and Thr and a lower content of hydrophobic amino acids, particularly Leu. Sequence data from nine methanogen genomes (OGT 15degrees-98degreesC) were used to generate IIII modeled protein structures. Analysis of the models from the cold-adapted Archaea showed a strong tendency in the solvent-accessible area for more Gin, Thr, and hydrophobic residues and fewer charged residues. A cold shock domain (CSD) protein (CspA homolog) was identified in M. frigidum, two hypothetical proteins with CSD-folds in M. burtonii, and a unique winged helix DNA-binding domain protein in M. burtonii. This suggests that these types of nucleic acid binding proteins have a critical role in cold-adapted Archaea. Structural analysis of tRNA sequences from the Archaea indicated that GC content is the major factor influencing tRNA stability in hyperthermophiles, but not in the psychrophiles, mesophiles or moderate thermophiles. Below an OGT of 60degreesC, the GC content in tRNA was largely unchanged, indicating that any requirement for flexibility of tRNA in psychrophiles is mediated by other means. This is the first time that comparisons have been performed with genome data from Archaea spanning the growth temperature extremes. from psychrophiles to hyperthermophiles

Scleractinian corals with photoprotective host pigments are hypersensitive to thermal bleaching

Recent episodes of mass coral bleaching, the loss of symbiotic dinoflagellates or photosynthetic pigment from hermatypic corals, have been triggered by elevated sea temperatures. Photosynthetic irradiance is an important secondary factor. Host based pigments (pocilloporins or Green Fluorescent Protein homologues) have been proposed to reduce the impact of elevated temperature by shading the dinoflagellate symbionts of corals, thereby reducing light stress. This study investigates this phenomenon in the reef-building coral Acropora aspera from Heron Island Research Station (Great Barrier Reef, Australia), which occurs as 3 distinct colour morphs. Experimental data showed that the host pigments are photoprotective at normal temperatures or

Prediction of protein B-factor profiles

The polypeptide backbones and side chains of proteins are constantly moving due to thermal motion and the kinetic energy of the atoms. The B-factors of protein crystal structures reflect the fluctuation of atoms about their average positions and provide important information about protein dynamics. Computational approaches to predict thermal motion are useful for analyzing the dynamic properties of proteins with unknown structures. In this article, we utilize a novel support vector regression (SVR) approach to predict the B-factor distribution (B-factor profile) of a protein from its sequence. We explore schemes for encoding sequences and various settings for the parameters used in SVR. Based on a large dataset of high-resolution proteins, our method predicts the B-factor distribution with a Pearson correlation coefficient (CC) of 0.53. In addition, our method predicts the B-factor profile with a CC of at least 0.56 for more than half of the proteins. Our method also performs well for classifying residues (rigid vs. flexible). For almost all predicted B-factor thresholds, prediction accuracies (percent of correctly predicted residues) are greater than 70%. These results exceed the best results of other sequence-based prediction methods. (C) 2005 Wiley-Liss, Inc.