Dynamic multiscale metabolic network modeling of Chinese hamster ovary cell metabolism integrating N-linked glycosylation in industrial biopharmaceutical manufacturing.

Experimental and modeling work, described in this article, is focused on the metabolic pathway of Chinese hamster ovary (CHO) cells, which are the preferred expression system for monoclonal antibody protein production. CHO cells are one of the primary hosts for monoclonal antibodies production, which have extensive applications in multiple fields like biochemistry, biology and medicine. Here, an approach to explain cellular metabolism with in silico modeling of a microkinetic reaction network is presented and validated with unique experimental results. Experimental data of 25 different fed-batch bioprocesses included the variation of multiple process parameters, such as pH, agitation speed, oxygen and CO2 content, and dissolved oxygen. A total of 151 metabolites were involved in our proposed metabolic network, which consisted of 132 chemical reactions that describe the reaction pathways, and include 25 reactions describing N-glycosylation and additional reactions for the accumulation of the produced glycoforms. Additional eight reactions are considered for accumulation of the N-glycosylation products in the extracellular environment and one reaction to correlate cell degradation. The following pathways were considered: glycolysis, pentose phosphate pathway, nucleotide synthesis, tricarboxylic acid cycle, lipid synthesis, protein synthesis, biomass production, anaplerotic reactions, and membrane transport. With the applied modeling procedure, different operational scenarios and fed-batch techniques can be tested.

Metabolomic profiling of CHO fed-batch growth phases at 10, 100, and 1,000 L.

Established bioprocess monitoring is based on quick and reliable methods, including cell count and viability measurement, extracellular metabolite measurement, and the measurement of physicochemical qualities of the cultivation medium. These methods are sufficient for monitoring of process performance, but rarely give insight into the actual physiological states of the cell culture. However, understanding of the latter is essential for optimization of bioprocess development. Our study used LC-MS metabolomics as a tool for additional resolution of bioprocess monitoring and was designed at three bioreactors scales (10 L, 100 L, and 1,000 L) to gain insight into the basal metabolic states of the Chinese hamster ovary (CHO) cell culture during fed-batch. Metabolites characteristics of the four growth stages (early and late exponential phase, stationary phase, and the phase of decline) were identified by multivariate analysis. Enriched metabolic pathways were then established for each growth phase using the CHO metabolic network model. Biomass generation and nucleotide synthesis were enriched in early exponential phase, followed by increased protein production and imbalanced glutathione metabolism in late exponential phase. Glycolysis became downregulated in stationary phase and amino-acid metabolism increased. Phase of culture decline resulted in rise of oxidized glutathione and fatty acid concentrations. Intracellular metabolic profiles of the CHO fed-batch culture were also shown to be consistent with scale and thus demonstrate metabolomic profiling as an informative method to gain physiological insight into the cell culture states during bioprocess regardless of scale.

The Bark-Beetle-Associated Fungus, Endoconidiophora polonica, Utilizes the Phenolic Defense Compounds of Its Host as a Carbon Source.

Norway spruce (Picea abies) is periodically attacked by the bark beetle Ips typographus and its fungal associate, Endoconidiophora polonica, whose infection is thought to be required for successful beetle attack. Norway spruce produces terpenoid resins and phenolics in response to fungal and bark beetle invasion. However, how the fungal associate copes with these chemical defenses is still unclear. In this study, we investigated changes in the phenolic content of Norway spruce bark upon E. polonica infection and the biochemical factors mediating these changes. Although genes encoding the rate-limiting enzymes in Norway spruce stilbene and flavonoid biosynthesis were actively transcribed during fungal infection, there was a significant time-dependent decline of the corresponding metabolites in fungal lesions. In vitro feeding experiments with pure phenolics revealed that E. polonica transforms both stilbenes and flavonoids to muconoid-type ring-cleavage products, which are likely the first steps in the degradation of spruce defenses to substrates that can enter the tricarboxylic acid cycle. Four genes were identified in E. polonica that encode catechol dioxygenases carrying out these reactions. These enzymes catalyze the cleavage of phenolic rings with a vicinal dihydroxyl group to muconoid products accepting a wide range of Norway spruce-produced phenolics as substrates. The expression of these genes and E. polonica utilization of the most abundant spruce phenolics as carbon sources both correlated positively with fungal virulence in several strains. Thus, the pathways for the degradation of phenolic compounds in E. polonica, initiated by catechol dioxygenase action, are important to the infection, growth, and survival of this bark beetle-vectored fungus and may play a major role in the ability of I. typographus to colonize spruce trees.

Gene discovery for enzymes involved in limonene modification or utilization by the mountain pine beetle-associated pathogen Grosmannia clavigera.

To successfully colonize and eventually kill pine trees, Grosmannia clavigera (Gs cryptic species), the main fungal pathogen associated with the mountain pine beetle (Dendroctonus ponderosae), has developed multiple mechanisms to overcome host tree chemical defenses, of which terpenoids are a major component. In addition to a monoterpene efflux system mediated by a recently discovered ABC transporter, Gs has genes that are highly induced by monoterpenes and that encode enzymes that modify or utilize monoterpenes [especially (+)-limonene]. We showed that pine-inhabiting Ophiostomale fungi are tolerant to monoterpenes, but only a few, including Gs, are known to utilize monoterpenes as a carbon source. Gas chromatography-mass spectrometry (GC-MS) revealed that Gs can modify (+)-limonene through various oxygenation pathways, producing carvone, p-mentha-2,8-dienol, perillyl alcohol, and isopiperitenol. It can also degrade (+)-limonene through the C-1-oxygenated pathway, producing limonene-1,2-diol as the most abundant intermediate. Transcriptome sequencing (RNA-seq) data indicated that Gs may utilize limonene 1,2-diol through beta-oxidation and then valine and tricarboxylic acid (TCA) metabolic pathways. The data also suggested that at least two gene clusters, located in genome contigs 108 and 161, were highly induced by monoterpenes and may be involved in monoterpene degradation processes. Further, gene knockouts indicated that limonene degradation required two distinct Baeyer-Villiger monooxygenases (BVMOs), an epoxide hydrolase and an enoyl coenzyme A (enoyl-CoA) hydratase. Our work provides information on enzyme-mediated limonene utilization or modification and a more comprehensive understanding of the interaction between an economically important fungal pathogen and its host's defense chemicals.

The cytochromes P450 of Grosmannia clavigera: Genome organization, phylogeny, and expression in response to pine host chemicals.

Grosmannia clavigera is a fungal associate of the mountain pine beetle (Dendroctonus ponderosae) and a pathogen of lodgepole pine (Pinus contorta) that must overcome terpenoid oleoresin and phenolic defenses of host trees. G. clavigera responds to monoterpene influx with complementary mechanisms that include export and the use of these compounds as a carbon source. Cytochromes P450 (CYPs) may also be involved in the metabolism of host defense compounds. We have identified and phylogenetically classified G. clavigera CYPs (CYPome). We show that although the G. clavigera CYPome has contracted in evolution, certain CYP families have expanded by duplication. We analyzed RNA-seq data for CYP expression following treatment with terpenes and pine phloem extracts to identify CYPs potentially involved in detoxification of these pine defense compounds. We also used transcriptome analysis of G. clavigera grown on monoterpenes, triglycerides or oleic acid as a carbon source to identify up-regulated CYPs that may be involved in the utilization of these compounds to support fungal growth. Finally, we identify secondary metabolite biosynthetic gene clusters that contain CYPs, and CYPs in clusters that may be involved in conversion of host chemicals.

The versatility of the fungal cytochrome P450 monooxygenase system is instrumental in xenobiotic detoxification.

Cytochromes P450 (CYPs) catalyse diverse reactions and are key enzymes in fungal primary and secondary metabolism, and xenobiotic detoxification. CYP enzymatic properties and substrate specificity determine the reaction outcome. However, CYP-mediated reactions may also be influenced by their redox partners. Filamentous fungi with numerous CYPs often possess multiple microsomal redox partners, cytochrome P450 reductases (CPRs). In the plant pathogenic ascomycete Cochliobolus lunatus we recently identified two CPR paralogues, CPR1 and CPR2. Our objective was to functionally characterize two endogenous fungal cytochrome P450 systems and elucidate the putative physiological roles of CPR1 and CPR2. We reconstituted both CPRs with CYP53A15, or benzoate 4-hydroxylase from C. lunatus, which is crucial in the detoxification of phenolic plant defence compounds. Biochemical characterization using RP-HPLC shows that both redox partners support CYP activity, but with different product specificities. When reconstituted with CPR1, CYP53A15 converts benzoic acid to 4-hydroxybenzoic acid, and 3-methoxybenzoic acid to 3-hydroxybenzoic acid. However, when the redox partner is CPR2, both substrates are converted to 3,4-dihydroxybenzoic acid. Deletion mutants and gene expression in mycelia grown on media with inhibitors indicate that CPR1 is important in primary metabolism, whereas CPR2 plays a role in xenobiotic detoxification.

Virtual screening yields inhibitors of novel antifungal drug target, benzoate 4-monooxygenase.

Fungal CYP53 enzymes are highly conserved proteins, involved in phenolic detoxification, and have no homologues in higher eukaryotes, rendering them favorable drug targets. Aiming to discover novel CYP53 inhibitors, we employed two parallel virtual screening protocols and evaluated highest scoring hit compounds by analyzing the spectral binding interactions, by surveying the antifungal activity, and assessing the inhibition of catalytic activity. On the basis of combined results, we selected 3-methyl-4-(1H-pyrrol-1-yl)benzoic acid (compound 2) as the best candidate for hit-to-lead follow-up in the antifungal drug discovery process.

A comparative genomic analysis of calcium and proton signaling/homeostasis in Aspergillus species.

A large number of proteins involved in calcium and intracellular pH signaling and homeostasis have previously been discovered and characterized in Saccharomyces cerevisiae, Homo sapiens and Arabidopsis thaliana,but relatively few have been identified in Aspergillus species. The aim of this study was to identify proteins regulating the intracellular concentration of calcium ions and protons in Aspergillus spp. and compare these with other fungi. For Aspergillus spp. we identified 46, 97 and 105 putative Ca2+-permeable channels, cation/proton transporters and P-ATPases, respectively, the majority of them previously unknown. The subunits composing V-type H+ ATPase and F0F1 ATP synthase (F-type ATPase) from Aspergillus spp. were identified. The greater redundancy of Ca2+-permeable channels, cation/proton exchangers and P-ATPases in filamentous fungi (between 28 putative proteins from A. clavatus and 34 from A. oryzae)compared to that of S. cerevisiae (19 proteins) reflects a more complex cellular organization and filamentous growth form. On the other hand the complexity of V-type H+ ATPase and F0F1 ATP synthase in filamentous fungi is comparable to that in ascomycetous yeast species indicating that both ATPase complexes are a basic universal requirement of the fungal cell.

High diversity and complex evolution of fungal cytochrome P450 reductase: cytochrome P450 systems.

Cytochrome P450 reductase (CPR) is the redox partner of P450 monooxygenases, involved in primary and secondary metabolism of eukaryotes. Two novel CPR genes, sharing 34% amino acid identity, were found in the filamentous ascomycete Cochliobolus lunatus. Fungal genomes were searched for putative CPR enzymes. Phylogenetic analysis suggests that multiple independent CPR duplication events occurred in fungi, whereas P450-CPR fusion occurred before the diversification of Dikarya and Zygomycota. Additionally, losses of methionine synthase reductase were found in certain fungal taxa; a truncated form of this enzyme was conserved in Pezizomycotina. In fungi, high numbers of cytochrome P450 enzymes, multiple CPRs, and P450-CPR fusion proteins were associated with filamentous growth. Evolution of multiple CPR-like oxidoreductases in filamentous fungi might have been driven by the complexity of biochemical functions necessitated by their growth form, as opposed to yeast.

CYP53A15 of Cochliobolus lunatus, a target for natural antifungal compounds.

A novel cytochrome P450, CYP53A15, was identified in the pathogenic filamentous ascomycete Cochliobolus lunatus. The protein, classified into the CYP53 family, was capable of para hydroxylation of benzoate. Benzoate is a key intermediate in the metabolism of aromatic compounds in fungi and yet basically toxic to the organism. To guide functional analyses, protein structure was predicted by homology modeling. Since many naturally occurring antifungal phenolic compounds are structurally similar to CYP53A15 substrates, we tested their putative binding into the active site of CYP53A15. Some of these compounds inhibited CYP53A15. Increased antifungal activity was observed when tested in the presence of benzoate. Some results suggest that CYP53A15 O-demethylation activity is important in detoxification of other antifungal substances. With the design of potent inhibitors, CYP53 enzymes could serve as alternative antifungal drug targets.

High-quality whole-genome sequence of an abundant Holarctic odontocete, the harbour porpoise (Phocoena phocoena).

The harbour porpoise (Phocoena phocoena) is a highly mobile cetacean found across the Northern hemisphere. It occurs in coastal waters and inhabits basins that vary broadly in salinity, temperature and food availability. These diverse habitats could drive subtle differentiation among populations, but examination of this would be best conducted with a robust reference genome. Here, we report the first harbour porpoise genome, assembled de novo from an individual originating in the Kattegat Sea (Sweden). The genome is one of the most complete cetacean genomes currently available, with a total size of 2.39 Gb and 50% of the total length found in just 34 scaffolds. Using 122 of the longest scaffolds, we were able to show high levels of synteny with the genome of the domestic cattle (Bos taurus). Our draft annotation comprises 22,154 predicted genes, which we further annotated through matches to the NCBI nucleotide database, GO categorization and motif prediction. Within the predicted genes, we have confirmed the presence of >20 genes or gene families that have been associated with adaptive evolution in other cetaceans. Overall, this genome assembly and draft annotation represent a crucial addition to the genomic resources currently available for the study of porpoises and Phocoenidae evolution, phylogeny and conservation.

A genomic comparison of putative pathogenicity-related gene families in five members of the Ophiostomatales with different lifestyles.

Ophiostomatoid fungi are vectored by their bark-beetle associates and colonize different host tree species. To survive and proliferate in the host, they have evolved mechanisms for detoxification and elimination of host defence compounds, efficient nutrient sequestration, and, in pathogenic species, virulence towards plants. Here, we assembled a draft genome of the spruce pathogen Ophiostoma bicolor. For our comparative and phylogenetic analyses, we mined the genomes of closely related species (Ophiostoma piceae, Ophiostoma ulmi, Ophiostoma novo-ulmi, and Grosmannia clavigera). Our aim was to acquire a genomic and evolutionary perspective of gene families important in host colonization. Genome comparisons showed that both the nuclear and mitochondrial genomes in our assembly were largely complete. Our O. bicolor 25.3 Mbp draft genome had 10 018 predicted genes, 6041 proteins with gene ontology (GO) annotation, 269 carbohydrate-active enzymes (CAZymes), 559 peptidases and inhibitors, and 1373 genes likely involved in pathogen-host interactions. Phylogenetic analyses of selected protein families revealed core sets of cytochrome P450 genes, ABC transporters and backbone genes involved in secondary metabolite (SM) biosynthesis (polyketide synthases (PKS) and non-ribosomal synthases), and species-specific gene losses and duplications. Phylogenetic analyses of protein families of interest provided insight into evolutionary adaptations to host biochemistry in ophiostomatoid fungi.

Spatially Explicit Analysis of Genome-Wide SNPs Detects Subtle Population Structure in a Mobile Marine Mammal, the Harbor Porpoise.

The population structure of the highly mobile marine mammal, the harbor porpoise (Phocoena phocoena), in the Atlantic shelf waters follows a pattern of significant isolation-by-distance. The population structure of harbor porpoises from the Baltic Sea, which is connected with the North Sea through a series of basins separated by shallow underwater ridges, however, is more complex. Here, we investigated the population differentiation of harbor porpoises in European Seas with a special focus on the Baltic Sea and adjacent waters, using a population genomics approach. We used 2872 single nucleotide polymorphisms (SNPs), derived from double digest restriction-site associated DNA sequencing (ddRAD-seq), as well as 13 microsatellite loci and mitochondrial haplotypes for the same set of individuals. Spatial principal components analysis (sPCA), and Bayesian clustering on a subset of SNPs suggest three main groupings at the level of all studied regions: the Black Sea, the North Atlantic, and the Baltic Sea. Furthermore, we observed a distinct separation of the North Sea harbor porpoises from the Baltic Sea populations, and identified splits between porpoise populations within the Baltic Sea. We observed a notable distinction between the Belt Sea and the Inner Baltic Sea sub-regions. Improved delineation of harbor porpoise population assignments for the Baltic based on genomic evidence is important for conservation management of this endangered cetacean in threatened habitats, particularly in the Baltic Sea proper. In addition, we show that SNPs outperform microsatellite markers and demonstrate the utility of RAD-tags from a relatively small, opportunistically sampled cetacean sample set for population diversity and divergence analysis.

Oleic acid metabolism via a conserved cytochrome P450 system-mediated ω-hydroxylation in the bark beetle-associated fungus Grosmannia clavigera.

The bark beetle-associated fungus Grosmannia clavigera participates in the large-scale destruction of pine forests. In the tree, it must tolerate saturating levels of toxic conifer defense chemicals (e.g. monoterpenes). The fungus can metabolize some of these compounds through the ß-oxidation pathway and use them as a source of carbon. It also uses carbon from pine triglycerides, where oleic acid is the most common fatty acid. High levels of free fatty acids, however, are toxic and can cause additional stress during host colonization. Fatty acids induce expression of neighboring genes encoding a cytochrome P450 (CYP630B18) and its redox partner, cytochrome P450 reductase (CPR2). The aim of this work was to study the function of this novel P450 system. Using LC/MS, we biochemically characterized CYP630 as a highly specific oleic acid ω-hydroxylase. We explain oleic acid specificity using protein interaction modeling. Our results underscore the importance of ω-oxidation when the main ß-oxidation pathway may be overwhelmed by other substrates such as host terpenoid compounds. Because this CYP-CPR gene cluster is evolutionarily conserved, our work has implications for metabolism studies in other fungi.

Aegerolysins: structure, function, and putative biological role.

Aegerolysins, discovered in fungi, bacteria and plants, are highly similar proteins with interesting biological properties. Certain aegerolysins possess antitumoral, antiproliferative, and antibacterial activities. Further possible medicinal applications include their use in the prevention of atherosclerosis, or as vaccines. Additional biotechnological value of fungal aegerolysins lies in their involvement in development, which could improve cultivation of commercially important edible mushrooms. Besides, new insights on microheterogeneity of raft-like membrane domains could be gained by using aegerolysins as specific markers in cell and molecular biology. Although the exact function of aegerolysins in their producing organisms remains to be explained, they are biochemically well characterized all-beta structured proteins sharing the following common features: low isoelectric points, similar molecular weights (15-17 kDa), and stability in a wide pH range.