Modeling Lignin Polymerization. I. Simulation Model of Dehydrogenation Polymers

Lignin is a heteropolymer that is thought to form in the cell wall by combinatorial radical coupling of monolignols. Here, we present a simulation model of in vitro lignin polymerization, based on the combinatorial coupling theory, which allows us to predict the reaction conditions controlling the primary structure of lignin polymers. Our model predicts two controlling factors for the β-O-4 content of syringyl-guaiacyl lignins: the supply rate of monolignols and the relative amount of supplied sinapyl alcohol monomers. We have analyzed the in silico degradability of the resulting lignin polymers by cutting the resulting lignin polymers at β-O-4 bonds. These are cleaved in analytical methods used to study lignin composition, namely thioacidolysis and derivatization followed by reductive cleavage, under pulping conditions, and in some lignocellulosic biomass pretreatments.Lignins are aromatic polymers that are predominantly present in secondarily thickened cell walls. These polymers make the cell wall rigid and impervious, allowing transport of water and nutrients through the vascular system and protecting plants against microbial invasion. Lignins are heterogeneous polymers derived from phenylpropanoid monomers, mainly the hydroxycinnamyl alcohols coniferyl alcohol (G-monomer) and sinapyl alcohol (S-monomer) and minor amounts of p-coumaryl alcohol (H-monomer). These monolignols differ in their degree of aromatic methoxylation (-OCH₃ group; Fig. 1). The resulting units in the lignin polymer are the guaiacyl (G), syringyl (S), and p-hydroxyphenyl (H) units. They are linked by a variety of chemical bonds (Fig. 2) that have different chemical properties (Boerjan et al., 2003; Ralph et al., 2004; Vanholme et al., 2008).Open in a separate windowFigure 1.Chemical structures of three monolignols. A, H-monomer (p-coumaryl alcohol). B, G-monomer (coniferyl alcohol). C, S-monomer (sinapyl alcohol). G- and S-monomers are considered in our simulations. The G-monomer is methoxylated (-OCH₃ group) on position 3, and the S-monomer is methoxylated on positions 3 and 5.Open in a separate windowFigure 2.Chemical structures resulting from the possible bonding between two monomers (A) or a monomer and the bindable end of an oligomer (B). X and Y in the monomers denote the absence (for a G-unit) or presence (for an S-unit) of a methoxyl group at position 5 (see Fig. 1). The red line indicates the bonds generated by couplings of the B position and B, 4, or 5 position.Lignification is the process by which monomers and/or oligomers are polymerized via radical coupling reactions and typically occurs after the polysaccharides have been laid down in the cell wall. Lignin composition varies among cell types and can even be different in individual cell wall layers (Ruel et al., 2009). Lignin composition is also influenced by environmental conditions; for example, lignin in compression wood is enriched in H-units (Timell, 1986). Hence, both developmental and environmental parameters influence the composition and thus the structure of the lignin polymer (Boerjan et al., 2003; Ralph et al., 2004).Lignin is one of the main negative factors in the conversion of lignocellulosic plant biomass into pulp and bioethanol (Lynd et al., 1991; Hill et al., 2006). In these processes, lignin needs to be degraded by chemical or mechanical processes that are expensive and often environmentally polluting. Hence, major research efforts are devoted toward understanding lignin biosynthesis and structure. It has already been shown that reducing lignin content and modifying its composition in transgenic plants can result in dramatic improvements in pulping efficiency (Pilate et al., 2002; Baucher et al., 2003; Huntley et al., 2003; Leplé et al., 2007) and in the conversion of biomass into bioethanol (Stewart et al., 2006; Chen and Dixon, 2007; Custers, 2009). These altered biomass properties are related to the alterations in lignin composition and structure in terms of the frequencies of the lignin units and the bond types connecting them and possibly also their interaction with hemicelluloses (Ralph et al., 2004; Ralph, 2006).To study the parameters that influence lignin structure, lignin polymerization has been mimicked in vitro by experiments with dehydrogenation polymers (DHPs; Terashima et al., 1995). Indeed, lignification can be mimicked by oxidizing monolignols using a peroxidase, such as horseradish peroxidase (HRP), and supplying its cofactor hydrogen peroxide, producing synthetic DHP lignins. Monolignol oxidation can also be achieved without enzymes (e.g. by using transition metal one-electron oxidants, such as copper acetate). Some of these biomimetic DHPs have been suggested to be better models for wood lignins than HRP-generated DHPs (Landucci, 2000).In DHP experiments, the monolignols are either added in bulk (Zulauf experiment) or dropwise (Zutropf experiment) to the reaction mixture, yielding lignin polymers with very different bond frequencies (Freudenberg, 1956). Zutropf experiments approach the in vivo formation of lignin, which depends on the slow introduction of monolignols into the wall matrix via diffusion to the site of incorporation (Hatfield and Vermerris, 2001). Because the exact reaction conditions are known, such in vitro experiments have provided insight into the lignification process in planta. In this way, numerous factors were shown to influence lignin structure, including the relative supply of the monolignols, the pH, the presence of polysaccharides, hydrogen peroxide concentrations, and cell wall matrix elements in general (Grabber et al., 2003; Vanholme et al., 2008).Computer simulations of lignin polymerization can help explain and predict lignin structure from low-level chemical kinetic factors, including subunit-coupling probabilities and monolignol synthesis rates. Such models are helpful in explaining the mechanism behind a range of controlling factors identified in the experimental work, including (1) the ratio of coniferyl versus sinapyl alcohol monolignols, (2) the monolignol supply rate, and (3) the abundance of alternative monomers present during lignin biosynthesis in mutants and transgenics. Thus, computer models will also help in suggesting new targets for controlled lignin biosynthesis.Here, we propose a simulation model of synthetic lignin polymerization that is based upon an emerging consensus from a variety of observations and derives from a series of previous models of lignin polymerization (Glasser and Glasser, 1974; Glasser et al., 1976; Jurasek, 1995; Roussel and Lim, 1995). Our model uses a symbolic grammar to describe a constructive dynamical system (Fontana, 1992) or a rule-based system (Feret et al., 2009) in which it is not necessary to define all possible products in advance. We assume that G- and S-monomers and newly formed oligomers couple in a well-mixed medium, depending on coupling rules and experimentally measured coupling probabilities. To develop the model, we have used information from DHP experiments rather than natural lignins, as they are formed in a well-mixed medium and their reaction conditions are well known (e.g. the influx rate of monomers). Using information from natural lignin would have further complicated our model, as the structures of natural lignin polymers are influenced by many factors, including the possible involvement of dirigent proteins (Davin and Lewis, 2005), steric hindrance by polysaccharides, spatiotemporal regulation, and modifications during isolation procedures (Boerjan et al., 2003; Ralph et al., 2004).Using our simulation models, we study how putative controlling factors of lignin primary structure, including the influx rate of monomers and the relative amount of S-monomers, affect in silico lignin synthesis, and we compare our predictions with in vitro experiments. To predict the degradability of lignins formed in our simulations, we apply an in silico thioacidolysis, which cleaves the polymers at their β-O-4 positions. This simulates the molecular action of two of the most used methods to analyze lignin composition, thioacidolysis (Lapierre, 1993; Baucher et al., 2003) and derivatization followed by reductive cleavage (Lu and Ralph, 1997). The G+S-monomer yield is often taken as a reflection of the fraction of units bound by β-O-4 bonds. Cleavage of β-O-4 bonds is also the most important reaction in kraft pulping of wood (Baucher et al., 2003). The model predicts from first principles (1) that DHP lignins formed under Zutropf conditions have a higher β-O-4 content than those formed under Zulauf conditions, (2) that DHP lignins formed with high S content have a higher β-O-4 content than those formed with high G content, and (3) that a higher β-O-4 content does not necessarily reduce the average length of lignin fragments generated during in silico thioacidolysis.     

Accelerated hand bone mineral density loss is associated with progressive joint damage in hands and feet in recent-onset rheumatoid arthritis

Introduction

To investigate whether accelerated hand bone mineral density (BMD) loss is associated with progressive joint damage in hands and feet in the first year of rheumatoid arthritis (RA) and whether it is an independent predictor of subsequent progressive total joint damage after 4 years.

Methods

In 256 recent-onset RA patients, baseline and 1-year hand BMD was measured in metacarpals 2-4 by digital X-ray radiogrammetry. Joint damage in hands and feet were scored in random order according to the Sharp-van der Heijde method at baseline and yearly up to 4 years.

Results

68% of the patients had accelerated hand BMD loss (>-0.003 g/cm²) in the first year of RA. Hand BMD loss was associated with progressive joint damage after 1 year both in hands and feet with odds ratios (OR) (95% confidence intervals [CI]) of 5.3 (1.3-20.9) and 3.1 (1.0-9.7). In univariate analysis, hand BMD loss in the first year was a predictor of subsequent progressive total joint damage after 4 years with an OR (95% CI) of 3.1 (1.3-7.6). Multivariate analysis showed that only progressive joint damage in the first year and anti-citrullinated protein antibody positivity were independent predictors of long-term progressive joint damage.

Conclusions

In the first year of RA, accelerated hand BMD loss is associated with progressive joint damage in both hands and feet. Hand BMD loss in the first year of recent-onset RA predicts subsequent progressive total joint damage, however not independent of progressive joint damage in the first year.     

The effect of immunomodulators on the immunogenicity of TNF-blocking therapeutic monoclonal antibodies: a review

Therapeutic monoclonal antibodies have revolutionized the treatment of various inflammatory diseases. Immunogenicity against these antibodies has been shown to be clinically important: it is associated with shorter response duration because of diminishing concentrations in the blood and with infusion reactions. Concomitant immunomodulators in the form of methotrexate or azathioprine reduced the immunogenicity of therapeutic antibodies in rheumatoid arthritis, Crohn disease, and juvenile idiopathic arthritis. The occurrence of adverse events does not increase when immunomodulators are added to therapeutic antibodies. The mechanism whereby methotrexate and azathioprine influence immunogenicity remains unclear. Evidence-based consensus on prescribing concomitant immunomodulators is needed.     

Pseudotyped recombinant adeno-associated viral vectors mediate efficient gene transfer into primary human CD34+ peripheral blood progenitor cells

Background and aimsBecause of their pluripotency, human CD34⁺ peripheral blood progenitor cells (PBPC) are targets of interest for the treatment of many acquired and inherited disorders using gene therapeutic approaches. Unfortunately, most current vector systems lack either sufficient transduction efficiency or an appropriate safety profile. Standard single-stranded recombinant adeno-associated virus 2 (AAV2)-based vectors offer an advantageous safety profile, yet lack the required efficiency in human PBPC.MethodsA panel of pseudotyped AAV vectors (designated AAV2/x, containing the vector genome of serotype 2 and capsid of serotype x, AAV2/1–AAV2/6) was screened on primary human granulocyte–colony-stimulating factor (G-CSF)-mobilized CD34⁺ PBPC to determine their gene transfer efficacy. Additionally, double-stranded self-complementary AAV (dsAAV) were used to determine possible second-strand synthesis limitations.ResultsAAV2/6 vectors proved to be the most efficient [12.8% (1.8–25.4%) transgene-expressing PBPC after a single transduction], being significantly more efficient (all P < 0.005) than the other vectors [AAV2/2, 2.0% (0.2–7.3%); AAV2/1, 1.3% (0.1–2.9%); others, <; 1% transgene-expressing PBPC]. In addition, the relevance of the single-to-double-strand conversion block in transduction of human PBPC could be shown using pseudotyped dsAAV vectors: for dsAAV2/2 [9.3% (8.3–20.3%); P < 0.001] and dsAAV2/6 [37.7% (23.6–61.0%); P < 0.001) significantly more PBPC expressed the transgene compared with their single-stranded counterparts; for dsAAV2/1, no significant increase could be observed.ConclusionsWe have shown that clinically relevant transduction efficiency levels using AAV-based vectors in human CD34⁺ PBPC are feasible, thereby offering an efficient alternative vector system for gene transfer into this important target cell population.     

Kinetic model of sucrose accumulation in maturing sugarcane culm tissue

Biochemically, it is not completely understood why or how commercial varieties of sugarcane (Saccharum officinarum) are able to accumulate sucrose in high concentrations. Such concentrations are obtained despite the presence of sucrose synthesis/breakdown cycles (futile cycling) in the culm of the storage parenchyma. Given the complexity of the process, kinetic modelling may help to elucidate the factors governing sucrose accumulation or direct the design of experimental optimisation strategies. This paper describes the extension of an existing model of sucrose accumulation (Rohwer, J.M., Botha, F.C., 2001. Analysis of sucrose accumulation in the sugar cane culm on the basis of in vitro kinetic data. Biochem. J. 358, 437-445) to account for isoforms of sucrose synthase and fructokinase, carbon partitioning towards fibre formation, and the glycolytic enzymes phosphofructokinase (PFK), pyrophosphate-dependent PFK and aldolase. Moreover, by including data on the maximal activity of the enzymes as measured in different internodes, a growth model was constructed that describes the metabolic behaviour as sugarcane parenchymal tissue matures from internodes 3-10. While there was some discrepancy between modelled and experimentally determined steady-state sucrose concentrations in the cytoplasm, steady-state fluxes showed a better fit. The model supports a hypothesis of vacuolar sucrose accumulation against a concentration gradient. A detailed metabolic control analysis of sucrose synthase showed that each isoform has a unique control profile. Fructose uptake by the cell and sucrose uptake by the vacuole had a negative control on the futile cycling of sucrose and a positive control on sucrose accumulation, while the control profile for neutral invertase was reversed. When the activities of these three enzymes were changed from their reference values, the effects on futile cycling and sucrose accumulation were amplified. The model can be run online at the JWS Online database (http://jjj.biochem.sun.ac.za/database/uys).     

Binding-site identification and genotypic profiling of hepatitis C virus polymerase inhibitors

The search for hepatitis C virus polymerase inhibitors has resulted in the identification of several nonnucleoside binding pockets. The shape and nature of these binding sites differ across and even within diverse hepatitis C virus genotypes. These differences confront antiviral drug discovery with the challenge of finding compounds that are capable of inhibition in variable binding pockets. To address this, we have established a hepatitis C virus mutant and genotypic recombinant polymerase panel as a means of guiding medicinal chemistry through the elucidation of the site of action of novel inhibitors and profiling against genotypes. Using a genotype 1b backbone, we demonstrate that the recombinant P495L, M423T, M414T, and S282T mutant enzymes can be used to identify the binding site of an acyl pyrrolidine analog. We assess the inhibitory activity of this analog and other nonnucleoside inhibitors with our panel of enzyme isolates generated from clinical sera representing genotypes 1a, 1b, 2a, 2b, 3a, 4a, 5a, and 6a.     

Elevated amount of Toll-like receptor 4 mRNA in bronchial epithelial cells is associated with airway inflammation in horses with recurrent airway obstruction

Recurrent airway obstruction (RAO) is characterized by neutrophilic airway inflammation and obstruction, and stabling of susceptible horses triggers acute disease exacerbations. Stable dust is rich in endotoxin, which is recognized by Toll-like receptor (TLR) 4. In human bronchial epithelium, TLR4 stimulation leads to elevation of interleukin (IL)-8 mRNA expression. The zinc finger protein A20 negatively regulates this pathway. We hypothesized that TLR4 and IL-8 mRNA and neutrophil numbers are elevated and that A20 mRNA is not increased in RAOs during stabling compared with controls and with RAOs on pasture. We measured the maximal change in pleural pressure (DeltaPpl(max)), determined inflammatory cell counts in bronchoalveolar lavage fluid (BAL), and quantified TLR4, IL-8, and A20 mRNA in bronchial epithelium by quantitative RT-PCR. We studied six horse pairs, each pair consisting of one RAO and one control horse. Each pair was studied when the RAO-affected horse had airway obstruction induced by stabling and after 7, 14, and 28 days on pasture. Stabling increased BAL neutrophils, DeltaPpl(max), and TLR4 (4.14-fold change) significantly in RAOs compared with controls and with RAOs on pasture. TLR4 correlated with IL-8 (R2 = 0.75). Whereas stabling increased IL-8 in all horses, A20 was unaffected. IL-8 was positively correlated with BAL neutrophils (R2 = 0.43) and negatively with A20 (R2 = 0.44) only in RAO-affected horses. Elevated TLR4 expression and lack of A20 upregulation in bronchial epithelial cells from RAO-affected horses may contribute to elevated IL-8 production, leading to exaggerated neutrophilic airway inflammation in response to inhalation of stable dust.     

The phylogeny and evolution of host choice in the Hippoboscoidea (Diptera) as reconstructed using four molecular markers

Hippoboscoidea is a superfamily of Diptera that contains the Glossinidae or tsetse flies, the Hippoboscidae or louse flies, and two families of bat flies, the Streblidae and the Nycteribiidae. We reconstruct the phylogenetic relationships within Hippoboscoidea using maximum parsimony and Bayesian methods based on nucleotide sequences from fragments of four genes: nuclear 28S ribosomal DNA and the CPSase domain of CAD, and mitochondrial 16S rDNA and cytochrome oxidase I. We recover monophyly for most of the presently recognized groups within Hippoboscoidea including the superfamily as a whole, the Hippoboscidae, the Nycteribiidae, the bat flies, and the Pupipara (=Hippoboscidae+Nycteribiidae+Streblidae), as well as several subfamilies within the constituent families. Streblidae appear to be paraphyletic. Our phylogenetic hypothesis is well supported and decisive in that most competing topological hypotheses for the Hippoboscoidea require significantly longer trees. We confirm a single shift from a free-living fly to a blood-feeding ectoparasite of vertebrates and demonstrate that at least two host shifts from mammals to birds have occurred. Wings have been repeatedly lost, but never regained. The hippoboscoid ancestor also evolved adenotrophic viviparity and our cladogram is consistent with a gradual reduction in the motility of the deposited final instar larvae from active burrowing in the soil to true pupiparity where adult females glue the puparium within the confines of bat roosts.     

Maize streak virus-resistant transgenic maize: a first for Africa

In this article, we report transgene-derived resistance in maize to the severe pathogen maize streak virus (MSV). The mutated MSV replication-associated protein gene that was used to transform maize showed stable expression to the fourth generation. Transgenic T₂ and T₃ plants displayed a significant delay in symptom development, a decrease in symptom severity and higher survival rates than non-transgenic plants after MSV challenge, as did a transgenic hybrid made by crossing T₂ Hi-II with the widely grown, commercial, highly MSV-susceptible, white maize genotype WM3. To the best of our knowledge, this is the first maize to be developed with transgenic MSV resistance and the first all-African-produced genetically modified crop plant.