Increased peptide deformylase activity for N-formylmethionine processing of proteins overexpressed in Escherichia coli: application to homogeneous rubredoxin production

Activity-based probes for rhomboid proteases discovered in a mass spectrometry-based assay

Rhomboid proteases are evolutionary conserved intramembrane serine proteases. Because of their emerging role in many important biological pathways, rhomboids are potential drug targets. Unfortunately, few chemical tools are available for their study. Here, we describe a mass spectrometry-based assay to measure rhomboid substrate cleavage and inhibition. We have identified isocoumarin inhibitors and developed activity-based probes for rhomboid proteases. The probes can distinguish between active and inactive rhomboids due to covalent, reversible binding of the active-site serine and stable modification of a histidine residue. Finally, the structure of an isocoumarin-based inhibitor with Escherichia coli rhomboid GlpG uncovers an unusual mode of binding at the active site and suggests that the interactions between the 3-substituent on the isocoumarin inhibitor and hydrophobic residues on the protease reflect S' subsite binding. Overall, these probes represent valuable tools for rhomboid study, and the structural insights may facilitate future inhibitor design.

Activation of BlaR1 protein of methicillin-resistant Staphylococcus aureus, its proteolytic processing, and recovery from induction of resistance

The fates of BlaI, the gene repressor protein for the bla operon, BlaR1, the β-lactam sensor/signal transducer, and PC1 β-lactamase in four strains of Staphylococcus aureus upon exposure to four different β-lactam antibiotics were investigated as a function of time. The genes for the three proteins are encoded by the bla operon, the functions of which afford inducible resistance to β-lactam antibiotics in S. aureus. BlaR1 protein is expressed at low copy number. Acylation of the sensor domain of BlaR1 by β-lactam antibiotics initiates signal transduction to the cytoplasmic domain, a zinc protease, which is activated and degrades BlaI. This proteolytic degradation derepresses transcription of all three genes, resulting in inducible resistance. These processes take place within minutes of exposure to the antibiotics. The BlaR1 protein was shown to undergo fragmentation in three S. aureus strains within the time frame relevant for manifestation of resistance and was below the detection threshold in the fourth. Two specific sites of fragmentation were identified, one cytoplasmic and the other in the sensor domain. This is proposed as a means for turnover, a process required for recovery from induction of resistance in S. aureus in the absence of the antibiotic challenge. In S. aureus not exposed to β-lactam antibiotics (i.e. not acylated by antibiotic) the same fragmentation of BlaR1 is still observed, including the shedding of the sensor domain, an observation that leads to the conclusion that the sites of proteolysis might have evolved to predispose the protein to degradation within a set period of time.

Characterization of unusual families of ATG8-like proteins and ATG12 in the protozoan parasite Leishmania major

Leishmania major possesses, apparently uniquely, four families of ATG8-like genes, designated ATG8, ATG8A, ATG8B and ATG8C, and 25 genes in total. L. major ATG8 and examples from the ATG8A, ATG8B and ATG8C families are able to complement a Saccharomyces cerevisiae ATG8-deficient strain, indicating functional conservation. Whereas ATG8 has been shown to form putative autophagosomes during differentiation and starvation of L. major, ATG8A primarily form puncta in response to starvation-suggesting a role for ATG8A in starvation-induced autophagy. Recombinant ATG8A was processed at the scissile glycine by recombinant ATG4.2 but not ATG4.1 cysteine peptidases of L. major and, consistent with this, ATG4.2-deficient L. major mutants were unable to process ATG8A and were less able to withstand starvation than wild-type cells. GFP-ATG8-containing puncta were less abundant in ATG4.2 overexpression lines, in which unlipidated ATG8 predominated, which is consistent with ATG4.2 being an ATG8-deconjugating enzyme as well as an ATG8A-processing enzyme. In contrast, recombinant ATG8, ATG8B and ATG8C were all processed by ATG4.1, but not by ATG4.2. ATG8B and ATG8C both have a distinct subcellular location close to the flagellar pocket, but the occurrence of the GFP-labeled puncta suggest that they do not have a role in autophagy. L. major genes encoding possible ATG5, ATG10 and ATG12 homologues were found to complement their respective S. cerevisiae mutants, and ATG12 localized in part to ATG8-containing puncta, suggestive of a functional ATG5-ATG12 conjugation pathway in the parasite. L. major ATG12 is unusual as it requires C-terminal processing by an as yet unidentified peptidase.

A billion-fold range in acidity for the solvent-exposed amides of Pyrococcus furiosus rubredoxin

The exchange rates of the static solvent-accessible amide hydrogens of Pyrococcus furiosus rubredoxin range from near the diffusion-limited rate to a billion-fold slower for the non-hydrogen-bonded Val 38 (eubacterial numbering). Hydrogen exchange directly monitors the kinetic acidity of the peptide nitrogen. Electrostatic solvation free energies were calculated by Poisson-Boltzmann methods for the individual peptide anions that form during the hydroxide-catalyzed exchange reaction to examine how well the predicted thermodynamic acidities match the experimentally determined kinetic acidities. With the exception of the Ile 12 amide, the differential exchange rate constant for each solvent-exposed amide proton that is not hydrogen bonded to a backbone carbonyl can be predicted within a factor of 6 (10 (0.78)) root-mean-square deviation (rmsd) using the CHARMM22 electrostatic parameter set and an internal dielectric value of 3. Under equivalent conditions, the PARSE parameter set yields a larger rmsd value of 1.28 pH units, while the AMBER parm99 parameter set resulted in a considerably poorer correlation. Either increasing the internal dielectric value to 4 or reducing it to a value of 2 significantly degrades the quality of the prediction. Assigning the excess charge of the peptide anion equally between the peptide nitrogen and the carbonyl oxygen also reduces the correlation to the experimental data. These continuum electrostatic calculations were further analyzed to characterize the specific structural elements that appear to be responsible for the wide range of peptide acidities observed for these solvent-exposed amides. The striking heterogeneity in the potential at sites along the protein-solvent interface should prove germane to the ongoing challenge of quantifying the contribution that electrostatic interactions make to the catalytic acceleration achieved by enzymes.

Processing of N-terminal unnatural amino acids in recombinant human interferon-beta in Escherichia coli

Incorporation of unnatural amino acids into recombinant proteins represents a powerful tool for protein engineering and protein therapeutic development. While the processing of the N-terminal methionine (Met) residues in proteins is well studied, the processing of unnatural amino acids used for replacing the N-terminal Met remains largely unknown. Here we report the effects of the penultimate residue (the residue after the initiator Met) on the processing of two unnatural amino acids, L-azidohomoalanine (AHA) and L-homopropargylglycine (HPG), at the N terminus of recombinant human interferon-beta in E. coli. We have identified specific amino acids at the penultimate position that can be used to efficiently retain or remove N-terminal AHA or HPG. Retention of N-terminal AHA or HPG can be achieved by choosing amino acids with large side chains (such as Gln, Glu, and His) at the penultimate position, while Ala can be selected for the removal of N-terminal AHA or HPG. Incomplete processing of N-terminal AHA and HPG (in terms of both deformylation and cleavage) was observed with Gly or Ser at the penultimate position.

NMR and X-ray analysis of structural additivity in metal binding site-swapped hybrids of rubredoxin

Background

Chimeric hybrids derived from the rubredoxins of Pyrococcus furiosus (Pf) and Clostridium pasteurianum (Cp) provide a robust system for the characterization of protein conformational stability and dynamics in a differential mode. Interchange of the seven nonconserved residues of the metal binding site between the Pf and Cp rubredoxins yields a complementary pair of hybrids, for which the sum of the thermodynamic stabilities is equal to the sum for the parental proteins. Furthermore, the increase in amide hydrogen exchange rates for the hyperthermophile-derived metal binding site hybrid is faithfully mirrored by a corresponding decrease for the complementary hybrid that is derived from the less thermostable rubredoxin, indicating a degree of additivity in the conformational fluctuations that underlie these exchange reactions.

Results

Initial NMR studies indicated that the structures of the two complementary hybrids closely resemble "cut-and-paste" models derived from the parental Pf and Cp rubredoxins. This protein system offers a robust opportunity to characterize differences in solution structure, permitting the quantitative NMR chemical shift and NOE peak intensity data to be analyzed without recourse to the conventional conversion of experimental NOE peak intensities into distance restraints. The intensities for 1573 of the 1652 well-resolved NOE crosspeaks from the hybrid rubredoxins were statistically indistinguishable from the intensities of the corresponding parental crosspeaks, to within the baseplane noise level of these high sensitivity data sets. The differences in intensity for the remaining 79 NOE crosspeaks were directly ascribable to localized dynamical processes. Subsequent X-ray analysis of the metal binding site-swapped hybrids, to resolution limits of 0.79 Å and 1.04 Å, demonstrated that the backbone and sidechain heavy atoms in the NMR-derived structures lie within the range of structural variability exhibited among the individual molecules in the crystallographic asymmetric unit (~0.3 Å), indicating consistency with the "cut-and-paste" structuring of the hybrid rubredoxins in both crystal and solution.

Conclusion

Each of the significant energetic interactions in the metal binding site-swapped hybrids appears to exhibit a 1-to-1 correspondence with the interactions present in the corresponding parental rubredoxin structure, thus providing a structural basis for the observed additivity in conformational stability and dynamics. The congruence of these X-ray and NMR experimental data offers additional support for the interpretation that the conventional treatment of NOE distance restraints contributes substantially to the systematic differences that are commonly reported between NMR- and X-ray-derived protein structures.

Tris is a non-innocent buffer during intein-mediated protein cleavage

Fusion protein purification systems based on self-cleavable protein splicing elements are well established nowadays and have the advantage of producing recombinant proteins with their native amino acid composition while abolishing the need of an additional proteolytic cleavage step for removal of a purification tag. However, a potential disadvantage is the concomitant generation of reactive thioester intermediates during the protein self-splicing process, which are prone to undergo side reactions yielding undesired adducts. We followed the formation of these adducts as well as ways to avoid them with electrospray ionization mass spectrometry using one of our target proteins, Triticum aestivum (wheat) E(c)-1, a plant metallothionein with the ability to bind a total of six zinc or cadmium ions in the form of metal-thiolate clusters. Our investigations show that one of the most commonly used buffer substances, tris(hydroxymethyl)aminomethane (Tris), has to be applied with caution in combination with the described purification system, as it can itself react with the thioester intermediate forming a yet unreported stable adduct. This makes Tris a so called non-innocent buffer during the protein isolation procedure. Additionally, the results presented open up an interesting possibility to directly couple the one-step purification strategy with selective carboxy-terminal protein or peptide modification, e.g. the addition of fluorophors or PEGylation of peptides. Unrelated to the purification system used, we further observed a high amount of N-formylmethionine in the mass spectra when the protein of interest was expressed in cadmium-supplemented growth media.

Production and comprehensive quality control of recombinant human Interleukin-1beta: a case study for a process development strategy

We describe an efficient strategy to produce high-quality proteins by using a single large IMAC chromatography column and enzymatic His-tag removal via the TAGZyme system in pilot scale. Numerous quality assays demonstrated a high purity of the final product, the human cytokine Interleukin-1beta (IL-1beta). The protein preparation was apparently free of host cell proteins, endotoxins, protease, and aggregates. The N-terminal amino acid sequence of IL-1beta was in full agreement with the natural mature form of IL-1beta. The homogeneity of the product was further shown by X-ray structure determination which confirmed the previously solved structure of the protein. We propose the applied workflow as a strategy for industrial production of protein-based biopharmaceuticals.

Normal carbon acid referencing for protein amide hydrogen exchange

Measurement of the exchange kinetics for amide hydrogens along the protein backbone continues to offer valuable insights into structural stability and conformational dynamics. Since such studies routinely compare samples that differ in solution conditions or mechanical handling, normalization of the relative exchange rates can present a potentially significant source of experimental uncertainty. The carbon acids 1,3-dimethylimidazolium cation and thiomethylacetonitrile exhibit base catalyzed exchange rates similar to those of the slowly exchanging amides, under conditions typical for protein studies. With 13C enrichment at the acidic carbon position to facilitate selective observation, such carbon acids offer practical internal calibration of exchange.

Spatial distribution of dielectric shielding in the interior of Pyrococcus furiosus rubredoxin as sampled in the subnanosecond timeframe by hydrogen exchange

Experimental pK values of ionizable sidechains provide the most direct test for models representing dielectric shielding within the interior of a protein. However, only the strongly shifted pK values are particularly useful for discriminating among models. NMR titration studies have usually found only one or two such shifted pK values in each protein, so that the fitting of the experimental data to a uniform internal dielectric (epsilon(int)) model is not well constrained. The observed variation among proteins for such epsilon(int) estimates may reflect nonuniformity of dielectric shielding within each protein interior or qualitative differences between individual proteins. The differential amide kinetic acidities for a series of metal-substituted rubredoxins are shown to be consistent with Poisson-Boltzmann predictions of dielectric shielding that is relatively uniform for all of the amides that are sensitive to the metal charge, a region which corresponds to roughly 1/3 of the internal volume. The effective epsilon(int) values near 6 that are found in this study are significantly lower than many such estimates derived from sidechain pK measurements. The differing timeframes in which dielectric relaxation can respond to the highly transient peptide anion as compared to the longer lived states of the charged sidechains offers an explanation for the lower apparent dielectric constant deduced from these measurements.

Residue cluster additivity of thermodynamic stability in the hydrophobic core of mesophile vs. hyperthermophile rubredoxins

The branched sidechain residues 24 and 33 in the hydrophobic core of rubredoxin differ between the Clostridium pasteurianum (Cp) and Pyrococcus furiosus (Pf) sequences. Their X-ray structures indicate that these two sidechains are in van der Waals contact with each other, while neither appears to significantly interact with the other nonconserved residues. The simultaneous interchange of residues 24 and 33 between the Cp and Pf rubredoxin sequences yield a complementary pair of hybrid proteins for which the sum of their thermodynamic stabilities equals that of the parental rubredoxins. The 1.2 kcal/mol change arising from this two residues interchange accounts for 21% of the differential thermodynamic stability between the mesophile and hyperthermophile proteins. The additional interchange of the sole nonconserved aromatic residue in the hydrophobic core yields a 0.78 kcal/mol deviation from thermodynamic additivity.

Enhanced thermal stability achieved without increased conformational rigidity at physiological temperatures: Spatial propagation of differential flexibility in rubredoxin hybrids

The extreme thermal stability of proteins from hyperthermophilic organisms is widely believed to arise from an increased conformational rigidity in the native state. In apparent contrast to this paradigm, both Pyrococcus furiosus (Pf) rubredoxin, the most thermostable protein characterized to date, and its Clostridium pasteurianum (Cp) mesophile homolog undergo a transient conformational opening of their multi-turn segments, which is more favorable in hyperthermophile proteins below room temperature. Substitution of the hyperthermophile multi-turn sequence into the mesophile protein sequence yields a hybrid, (14-33(Pf)) Cp, that exhibits a 12 degrees increase in its reversible thermal unfolding transition midpoint. Nuclear magnetic resonance (NMR) magnetization transfer-based hydrogen exchange was used to monitor backbone conformational dynamics in the subsecond time regime. Despite the substantially increased thermostability, flexibility throughout the entire main chain of the more thermostable hybrid is equal to or greater than that of the wild type mesophile rubredoxin near its normal growth temperature. In comparison to the identical core residues of the (14-33(Pf)) Cp rubredoxin hybrid, six spatially clustered residues in the parental mesophile protein exhibit a substantially larger temperature dependence of exchange. The exchange behavior of these six residues closely matches that observed in the multi-turn segment, consistent with a more extensive conformational process. These six core residues exhibit a much weaker temperature dependence of exchange in the (14-33(Pf)) Cp hybrid, similar to that observed for the multi-turn segment in its parental Pf rubredoxin. These results suggest that differential temperature dependence of flexibility can underlie variations in thermostability observed for mesophile versus hyperthermophile homologs.

Additivity in Both Thermodynamic Stability and Thermal Transition Temperature for Rubredoxin Chimeras via Hybrid Native Partitioning

Given any operational definition of pairwise interaction, the set of residues that differ between two structurally homologous proteins can be uniquely partitioned into subsets of clusters for which no such interactions occur between clusters. Although hybrid protein sequences that preserve such clustering are consistent with tertiary structures composed of only parental native-like interactions, the stability of such predicted structures will depend upon the physical robustness of the assumed interaction potential. A simple distance cutoff criterion was applied to the most thermostable protein known to predict such a seven-residue cluster in the metal binding site region of Pyrococcus furiosus rubredoxin and a mesophile homolog. Both conformational stability and thermal transition temperature measurements demonstrate that 39% of the differential stability arises from these seven residues.

Contribution of the multi-turn segment in the reversible thermal stability of hyperthermophile rubredoxin: NMR thermal chemical exchange analysis of sequence hybrids

Pyrococcus furiosus (Pf) rubredoxin is the most thermostable protein characterized to date. Reflecting the complications arising from irreversible denaturation of this protein, predictions of which structural regions confer differential thermal stability have utilized kinetic stability measurements, hydrogen exchange protection factors, long range hydrogen bond NMR spin couplings, and molecular dynamics simulations, and have primarily implicated the three-stranded beta-sheet and the adjacent metal binding site. Herein, NMR chemical exchange experiments demonstrate reversible two-state unfolding at the thermal transition temperature (T(m)) for hybrids of Pf and the mesophile Clostridium pasteurianum (Cp) rubredoxins which interchange residues 14-33, the so-called multi-turn segment. This complementary pair of hybrid rubredoxins exhibits largely additive incremental thermal stabilizations vs. the parental proteins. Both stabilization free energy measurements as well as incremental T(m) values indicate that a minimum of 37% of the total differential thermal stability resides in this multi-turn segment. Such a proportionality between DeltaDeltaG and incremental T(m) values is predicted for hybrid pairs exhibiting thermodynamic additivity in which the differential stability is predominantly enthalpic.

Absence of kinetic thermal stabilization in a hyperthermophile rubredoxin indicated by 40 microsecond folding in the presence of irreversible denaturation

The striking kinetic stability of many proteins derived from hyperthermophilic organisms has led to the proposal that such stability may result from a heightened activation barrier for unfolding independent of a corresponding increase in the thermodynamic stability. This in turn implies a corresponding retardation of the folding reaction. A commonly cited model for kinetic thermal stabilization is the rubredoxin from Pyrococcus furiosus (Pf), which exhibits an irreversible denaturation lifetime at 100 degrees C of nearly a week. Utilizing protein resonances shifted well outside of the random coil chemical shift envelope, nuclear magnetic resonance (NMR) chemical exchange measurements on Pf rubredoxin as well as on the mesophile Clostridium pasteurianum (Cp) rubredoxin demonstrate reversible thermal transition temperatures of 144 degrees C (137 degrees C for the N-terminal modified A2K variant) and 104 degrees C, respectively, with similar (un)folding rates of approximately 25,000 s(-1), only modestly slower than the diffusion controlled rate. The absence of a substantial activation barrier to rubredoxin folding as well as the similar folding kinetics of the mesophile protein indicate that kinetic stabilization has not been utilized by the hyperthermophile rubredoxin in achieving its extreme thermal stability. The two-state folding kinetics observed for Pf rubredoxin contradict a previous assertion of multiphasic folding based on hydrogen exchange data extrapolated to an estimated midpoint of transition temperature (T(m)) of nearly 200 degrees C. This discrepancy is resolved by the observation that the base-catalyzed hydrogen exchange of the model dipeptide (N-acetyl-L-cysteine-N-methylamide)4-Cd2+ is 23-fold slower than that of the free cysteine model dipeptide used to normalize the Pf rubredoxin hydrogen exchange data.