Christopher M. Shuford

Membrane protein complexes catalyze both 4- and 3-hydroxylation of cinnamic acid derivatives in monolignol biosynthesis

The hydroxylation of 4- and 3-ring carbons of cinnamic acid derivatives during monolignol biosynthesis are key steps that determine the structure and properties of lignin. Individual enzymes have been thought to catalyze these reactions. In stem differentiating xylem (SDX) of Populus trichocarpa, two cinnamic acid 4-hydroxylases (PtrC4H1 and PtrC4H2) and a p-coumaroyl ester 3-hydroxylase (PtrC3H3) are the enzymes involved in these reactions. Here we present evidence that these hydroxylases interact, forming heterodimeric (PtrC4H1/C4H2, PtrC4H1/C3H3, and PtrC4H2/C3H3) and heterotrimeric (PtrC4H1/C4H2/C3H3) membrane protein complexes. Enzyme kinetics using yeast recombinant proteins demonstrated that the enzymatic efficiency (Vmax/km) for any of the complexes is 70–6,500 times greater than that of the individual proteins. The highest increase in efficiency was found for the PtrC4H1/C4H2/C3H3-mediated p-coumaroyl ester 3-hydroxylation. Affinity purification-quantitative mass spectrometry, bimolecular fluorescence complementation, chemical cross-linking, and reciprocal coimmunoprecipitation provide further evidence for these multiprotein complexes. The activities of the recombinant and SDX plant proteins demonstrate two protein-complex–mediated 3-hydroxylation paths in monolignol biosynthesis in P. trichocarpa SDX; one converts p-coumaric acid to caffeic acid and the other converts p-coumaroyl shikimic acid to caffeoyl shikimic acid. Cinnamic acid 4-hydroxylation is also mediated by the same protein complexes. These results provide direct evidence for functional involvement of membrane protein complexes in monolignol biosynthesis.

Complete Proteomic-Based Enzyme Reaction and Inhibition Kinetics Reveal How Monolignol Biosynthetic Enzyme Families Affect Metabolic Flux and Lignin in Populus trichocarpa

A proteomic-based predictive kinetic metabolic-flux model was developed for monolignol biosynthesis in Populus trichocarpa. Absolute quantities of all monolignol pathway proteins and 189 kinetic parameters were generated to construct the model, which was experimentally validated in transgenic P. trichocarpa and provides a comprehensive description of the monolignol biosynthetic pathway. We established a predictive kinetic metabolic-flux model for the 21 enzymes and 24 metabolites of the monolignol biosynthetic pathway using Populus trichocarpa secondary differentiating xylem. To establish this model, a comprehensive study was performed to obtain the reaction and inhibition kinetic parameters of all 21 enzymes based on functional recombinant proteins. A total of 104 Michaelis-Menten kinetic parameters and 85 inhibition kinetic parameters were derived from these enzymes. Through mass spectrometry, we obtained the absolute quantities of all 21 pathway enzymes in the secondary differentiating xylem. This extensive experimental data set, generated from a single tissue specialized in wood formation, was used to construct the predictive kinetic metabolic-flux model to provide a comprehensive mathematical description of the monolignol biosynthetic pathway. The model was validated using experimental data from transgenic P. trichocarpa plants. The model predicts how pathway enzymes affect lignin content and composition, explains a long-standing paradox regarding the regulation of monolignol subunit ratios in lignin, and reveals novel mechanisms involved in the regulation of lignin biosynthesis. This model provides an explanation of the effects of genetic and transgenic perturbations of the monolignol biosynthetic pathway in flowering plants.

Evaluation of Normalization Methods on GeLC-MS/MS Label-Free Spectral Counting Data to Correct for Variation during Proteomic Workflows

Normalization of spectral counts (SpCs) in label-free shotgun proteomic approaches is important to achieve reliable relative quantification. Three different SpC normalization methods, total spectral count (TSpC) normalization, normalized spectral abundance factor (NSAF) normalization, and normalization to selected proteins (NSP) were evaluated based on their ability to correct for day-to-day variation between gel-based sample preparation and chromatographic performance. Three spectral counting data sets obtained from the same biological conidia sample of the rice blast fungus Magnaporthe oryzae were analyzed by 1D gel and liquid chromatography-tandem mass spectrometry (GeLC-MS/MS). Equine myoglobin and chicken ovalbumin were spiked into the protein extracts prior to 1D-SDS- PAGE as internal protein standards for NSP. The correlation between SpCs of the same proteins across the different data sets was investigated. We report that TSpC normalization and NSAF normalization yielded almost ideal slopes of unity for normalized SpC versus average normalized SpC plots, while NSP did not afford effective corrections of the unnormalized data. Furthermore, when utilizing TSpC normalization prior to relative protein quantification, t-testing and fold-change revealed the cutoff limits for determining real biological change to be a function of the absolute number of SpCs. For instance, we observed the variance decreased as the number of SpCs increased, which resulted in a higher propensity for detecting statistically significant, yet artificial, change for highly abundant proteins. Thus, we suggest applying higher confidence level and lower fold-change cutoffs for proteins with higher SpCs, rather than using a single criterion for the entire data set. By choosing appropriate cutoff values to maintain a constant false positive rate across different protein levels (i.e., SpC levels), it is expected this will reduce the overall false negative rate, particularly for proteins with higher SpCs.

Coupling of a vented column with splitless nanoRPLC-ESI-MS for the improved separation and detection of brain natriuretic peptide-32 and its proteolytic peptides

The circulating concentration of a biomarker for congestive heart failure, Brain (B-type) Natriuretic Peptide (BNP-32), is measured using ELISA based assays in order to rapidly diagnose and monitor disease progression. The lack of molecular specificity afforded by these assays has recently come into question as emerging studies indicate there are potentially multiple heterogeneous forms of BNP in circulation with immunoreactive capabilities. In order to better understand the molecular biology of BNP-32 as it relates to congestive heart failure, it would thus be advantageous to use a detection platform such as Fourier transform ion cyclotron resonance mass spectrometry. This high resolving power mass spectrometer can provide unparalleled molecular specificity and can facilitate identification and characterization of the various molecular forms across all disease states. Unfortunately, BNP circulates at low concentrations (as low as 3fmol/mL). Thus, it will require a collaborative effort from a number of orthogonal front-end technologies to overcome the disconnect between the practical detection limits of this instrument platform and the physiological levels of BNP-32 and its alternative molecular forms. Herein, we begin optimization of these front-end techniques by first enhancing the conditions for online nanoLC-ESI-MS separations of BNP-32 and its proteolytic fragments. Through extensive analysis of various chromatographic parameters we determined that Michrom Magic C8 stationary phase used in conjunction with a continuous, vented column configuration provided advanced chromatographic performance for the nano-flow separations involving intact BNP-32 and its associated tryptic peptides. Furthermore, conditions for the tryptic digestion of BNP-32 were also studied. We demonstrate that the use of free cysteine as an alkylation quenching agent and a secondary digestion within the digestion scheme can provide targeted tryptic peptides with increased abundances. Combined, these data will serve to further augment the detection of BNP-32 by LC-MS.

Peptide production and decay rates affect the quantitative accuracy of protein cleavage isotope dilution mass spectrometry (PC-IDMS)

No consensus has been reached on the proper time to add stable-isotope labeled (SIL) peptides in protein cleavage isotope dilution mass spectrometry workflows. While quantifying 24 monolignol pathway enzymes in the xylem tissue of Populus trichocarpa, we compared the protein concentrations obtained when adding the SIL standard peptides concurrently with the enzyme or after quenching of the digestion (i.e. postdigestion) and observed discrepancies for nearly all tryptic peptides investigated. In some cases, greater than 30-fold differences were observed. To explain these differences and potentially correct for them, we developed a mathematical model based on pseudo-first-order kinetics to account for the dynamic production and decay (e.g. degradation and precipitation) of the native peptide targets in conjunction with the decay of the SIL peptide standards. A time course study of the digests confirmed the results predicted by the proposed model and revealed that the discrepancy between concurrent and postdigestion introduction of the SIL standards was related to differential decay experienced by the SIL peptide and the native peptide in each method. Given these results, we propose concurrent introduction of the SIL peptide is most appropriate, though not free from bias. Mathematical modeling of this method reveals that overestimation of protein quantities would still result when rapid peptide decay occurs and that this bias would be further exaggerated by slow proteolysis. We derive a simple equation to estimate the bias for each peptide based on the relative rates of production and decay. According to this equation, nearly half of the peptides evaluated here were estimated to have quantitative errors greater than 10% and in a few cases over 100%. We conclude that the instability of peptides can often significantly bias the protein quantities measured in protein cleavage isotope dilution mass spectrometry-based assays and suggest peptide stability be made a priority when selecting peptides to use for quantification.

Comprehensive Quantification of Monolignol-Pathway Enzymes in Populus trichocarpa by Protein Cleavage Isotope Dilution Mass Spectrometry

The economic value of wood/pulp from many tree species is largely dictated by the quantity and chemical properties of lignin, which is directly related to the composition and linkages of monolignols comprising the polymer. Although much is known regarding the monolignol biosynthetic pathway, our understanding is still deficient due to the lack of quantitative information at the proteomic level. We developed an assay based on protein cleavage isotope dilution mass spectrometry (PC-IDMS) for the determination of all potential, primary enzymes involved in the biosynthesis of monolignols and the peroxidases responsible for their polymerization to form lignin in the model tree species, Populus trichocarpa. Described is the identification of quantitative surrogate peptides through shotgun analysis of native and recombinant proteins, optimization of trypsin proteolysis using fractional factorial design of experiments, and development of a liquid chromatography-selected reaction monitoring method for specific detection of all targeted peptides. Of the 25 targeted enzymes, three were undetected in the normal xylem tissues, and all but two of the detectable species showed good day-to-day precision (CV < 10%). This represents the most comprehensive assay for quantification of proteins regulating monolignol biosynthesis and will lead to a better understanding of lignin formation at a systems level.

Systems Biology of Lignin Biosynthesis in Populus trichocarpa: Heteromeric 4-Coumaric Acid:Coenzyme A Ligase Protein Complex Formation, Regulation, and Numerical Modeling

This work shows that 4CL, an enzyme in monolignol biosynthesis, is found as a heterotetrameric complex of two isoforms in Populus trichocarpa. The activity of the heterotetramer can be described by a mathematical model that explains the effects of each isoform with mixtures of substrates and three types of inhibition, providing insights into the regulation of metabolic flux for this pathway. As a step toward predictive modeling of flux through the pathway of monolignol biosynthesis in stem differentiating xylem of Populus trichocarpa, we discovered that the two 4-coumaric acid:CoA ligase (4CL) isoforms, 4CL3 and 4CL5, interact in vivo and in vitro to form a heterotetrameric protein complex. This conclusion is based on laser microdissection, coimmunoprecipitation, chemical cross-linking, bimolecular fluorescence complementation, and mass spectrometry. The tetramer is composed of three subunits of 4CL3 and one of 4CL5. 4CL5 appears to have a regulatory role. This protein–protein interaction affects the direction and rate of metabolic flux for monolignol biosynthesis in P. trichocarpa. A mathematical model was developed for the behavior of 4CL3 and 4CL5 individually and in mixtures that form the enzyme complex. The model incorporates effects of mixtures of multiple hydroxycinnamic acid substrates, competitive inhibition, uncompetitive inhibition, and self-inhibition, along with characteristic of the substrates, the enzyme isoforms, and the tetrameric complex. Kinetic analysis of different ratios of the enzyme isoforms shows both inhibition and activation components, which are explained by the mathematical model and provide insight into the regulation of metabolic flux for monolignol biosynthesis by protein complex formation.

Monolignol Pathway 4-Coumaric Acid:Coenzyme A Ligases in Populus. trichocarpa: Novel Specificity, Metabolic Regulation, and Simulation of Coenzyme A Ligation Fluxes

Two 4-coumaric acid:CoA ligases regulate CoA flux with novel specificity. 4-Coumaric acid:coenzyme A ligase (4CL) is involved in monolignol biosynthesis for lignification in plant cell walls. It ligates coenzyme A (CoA) with hydroxycinnamic acids, such as 4-coumaric and caffeic acids, into hydroxycinnamoyl-CoA thioesters. The ligation ensures the activated state of the acid for reduction into monolignols. In Populus spp., it has long been thought that one monolignol-specific 4CL is involved. Here, we present evidence of two monolignol 4CLs, Ptr4CL3 and Ptr4CL5, in Populus trichocarpa. Ptr4CL3 is the ortholog of the monolignol 4CL reported for many other species. Ptr4CL5 is novel. The two Ptr4CLs exhibited distinct Michaelis-Menten kinetic properties. Inhibition kinetics demonstrated that hydroxycinnamic acid substrates are also inhibitors of 4CL and suggested that Ptr4CL5 is an allosteric enzyme. Experimentally validated flux simulation, incorporating reaction/inhibition kinetics, suggested two CoA ligation paths in vivo: one through 4-coumaric acid and the other through caffeic acid. We previously showed that a membrane protein complex mediated the 3-hydroxylation of 4-coumaric acid to caffeic acid. The demonstration here of two ligation paths requiring these acids supports this 3-hydroxylation function. Ptr4CL3 regulates both CoA ligation paths with similar efficiencies, whereas Ptr4CL5 regulates primarily the caffeic acid path. Both paths can be inhibited by caffeic acid. The Ptr4CL5-catalyzed caffeic acid metabolism, therefore, may also act to mitigate the inhibition by caffeic acid to maintain a proper ligation flux. A high level of caffeic acid was detected in stem-differentiating xylem of P. trichocarpa. Our results suggest that Ptr4CL5 and caffeic acid coordinately modulate the CoA ligation flux for monolignol biosynthesis.

Utilizing Spectral Counting To Quantitatively Characterize Tandem Removal of Abundant Proteins (TRAP) in Human Plasma

Biomarker discovery efforts in serum and plasma are greatly hindered by the presence of high abundance proteins that prevent the detection and quantification of less abundant, yet biologically significant, proteins. The most common method for addressing this problem is to specifically remove the few abundant proteins through immunoaffinity depletion/subtraction. Herein, we improved upon this method by utilizing multiple depletion columns in series, so as to increase the efficiency of the abundant protein removal and augment the detection/identification of less abundant plasma proteins. Spectral counting was utilized to make quantitative comparisons between undepleted plasma, plasma depleted with a single depletion column, and plasma depleted using two or three depletion columns in tandem. In the undepleted plasma only 29 lower abundance protein groups were identified with the top-scoring protein from each group having a median spectral count of 3, while in the plasma processed using a single HSA depletion column 61 such protein groups were identified with a median spectral count of 8. In comparison, 76 lesser abundant protein groups were identified with a median spectral count of 11.5 in the two column setup (i.e., HSA followed by MARS Hu14). However, in the ultimate depleted plasma sample, which was created using three depletion columns in tandem, the number of less abundant protein groups identified increase to 81 and the median spectral count for the top-scoring proteins from each group increased to 15 counts per protein. Moreover, exogenous B-type natriuretic peptide-32, which was added to the plasma as a detection benchmark at 12 μg/mL, was only detected in the plasma sample depleted using three depletion columns in tandem. Collectively, these data demonstrate that this method, tandem removal of abundant proteins or TRAP, provides superior removal efficiency compared to traditional applications and improves the depth of proteome coverage in plasma.

Improving limits of detection for B-type natriuretic peptide using PC-IDMS: An application of the ALiPHAT strategy

Hydrophobic tagging of biomolecules has been reported by our group and others to increase their ionization efficiency during electrospray ionization and facilitate their detection by mass spectrometry. As such, hydrophobic tagging should provide a viable method for augmenting MS-based quantification of low abundance proteins by decreasing their detection limits. Herein we have evaluated two commercial alkylation reagents and several newly synthesized hydrophobic alkylation reagents for their utility in quantifying B-type Natriuretic Peptide, a low abundance cardiac biomarker, by protein cleavage isotope dilution mass spectrometry. For the cysteine containing tryptic peptide evaluated, a approximately 3.5-fold decrease in the detection limit was observed for the best performing hydrophobic reagent, 2-iodo-N-octylacetamide, relative to the commonly used alkylation reagent, iodoacetamide. Additionally, we have evaluated the use of nonpolar surface areas as a metric for assessing the effectiveness of the alkylation reagents in improving ESI response.