Engineering topics
Publications and source records attributed to Shanklin, John.
This project aims to study the microbial structure and potential functions of bacterial and fungal microbiomes in leaves, stems, roots, rhizospheres, and bulk soils of energy crops (oilcane) grown in greenhouses.
In eukaryotes, Target of Rapamycin (TOR), a conserved protein sensor kinase, integrates a diverse set of environmental cues, including growth factor signals, energy availability, and nutritional status, to direct cell growth. In plants, TOR is activated by light and sugars and regulates a wide range of cellular processes, including protein synthesis and metabolism. Fatty acid synthesis is key to membrane biogenesis that in turn, is required for cell growth. To elucidate the primary regulatory role(s) of TOR in lipid metabolism, we followed fatty acid and lipid changes in plants with altered TOR protein levels or activity for short durations, using Nicotiana benthamiana leaves, Arabidopsis seedlings and Brassica napus cell suspension cultures. Transient expression of TOR significantly elevated the levels of total fatty acids in Nicotiana benthamiana leaves, while treatment of Arabidopsis seedlings with Torin 2, a TOR specific inhibitor, for one day, caused significant reductions in fatty acids and membrane lipids. Similarly, incubating oil-producing Brassica napus suspension culture cells with Torin 2 for eight hours led to significant decreases in the levels of TFA and TAG. Taken together the results from three independent systems presented here establishes that TOR positively regulates lipid synthesis in plants, consistent with its role in animals. Furthermore, RNA-seq analysis of Torin 2-treated Arabidopsis seedlings showed that TOR promotes the upregulation of a number of genes involved in de novo fatty acid synthesis while downregulating genes involved in lipid turnover, which we propose as a mechanistic explanation for its promotion of lipid synthesis and accumulation.
Lipid remodeling plays a critical role in plant response to abiotic stress and metabolic perturbations. Key steps in this process involve modifications of phosphatidylcholine (PC) acyl chains mediated by lysophosphatidylcholine: acyl-CoA acyltransferases (LPCATs) and phosphatidylcholine: diacylglycerol cholinephosphotransferase (ROD1). To assess their importance in lipid homeostasis, we took advantage of the trigalactosyldiacylglycerol1 (tgd1) mutant that exhibits marked increases in fatty acid synthesis and fatty acid flux through PC due to a block in inter-organelle lipid trafficking. Here, we showed that the increased fatty acid synthesis in tgd1 is due to posttranslational activation of the plastidic acetyl-coenzyme A carboxylase. Genetic analysis showed that knockout of LPCAT1 and 2 resulted in a lethal phenotype in tgd1. In addition, plants homozygous for lpcat2 and heterozygous for lpcat1 in the tgd1 background showed reduced levels of PC and triacylglycerols (TAG) and alterations in their fatty acid profiles. We further showed that disruption of ROD1 in tgd1 resulted in changes in fatty acid composition of PC and TAG, decreased leaf TAG content and reduced seedling growth. Combined, our results reveal a critical role of LPCATs and ROD1 in maintaining cellular lipid homeostasis under conditions, in which fatty acid production largely exceeds the cellular demand for membrane lipid synthesis.
Summary Camelina ( Camelina sativa L.), a hexaploid member of the Brassicaceae family, is an emerging oilseed crop being developed to meet the increasing demand for plant oils as biofuel feedstocks. In other Brassicas, high oil content can be associated with a yellow seed phenotype, which is unknown for camelina. We sought to create yellow seed camelina using CRISPR/Cas9 technology to disrupt its Transparent Testa 8 (TT8) transcription factor genes and to evaluate the resulting seed phenotype. We identified three TT8 genes, one in each of the three camelina subgenomes, and obtained independent CsTT8 lines containing frameshift edits. Disruption of TT8 caused seed coat colour to change from brown to yellow reflecting their reduced flavonoid accumulation of up to 44%, and the loss of a well‐organized seed coat mucilage layer. Transcriptomic analysis of CsTT8 ‐edited seeds revealed significantly increased expression of the lipid‐related transcription factors LEC1 , LEC2 , FUS3 , and WRI1 and their downstream fatty acid synthesis‐related targets. These changes caused metabolic remodelling with increased fatty acid synthesis rates and corresponding increases in total fatty acid (TFA) accumulation from 32.4% to as high as 38.0% of seed weight, and TAG yield by more than 21% without significant changes in starch or protein levels compared to parental line. These data highlight the effectiveness of CRISPR in creating novel enhanced‐oil germplasm in camelina. The resulting lines may directly contribute to future net‐zero carbon energy production or be combined with other traits to produce desired lipid‐derived bioproducts at high yields.
SUCROSE-NON-FERMENTING1-RELATED PROTEIN KINASE1 (SnRK1), a central plant metabolic sensor kinase, phosphorylates its target proteins, triggering a global shift from anabolism to catabolism. Molecular modeling revealed that upon binding of KIN10 to GEMINIVIRUS REP-INTERACTING KINASE1 (GRIK1), KIN10’s activation T-loop reorients into GRIK1’s active site, enabling its phosphorylation and activation. Trehalose 6-phosphate (T6P) is a proxy for cellular sugar status and a potent inhibitor of SnRK1. T6P binds to KIN10, a SnRK1 catalytic subunit, weakening its affinity for GRIK1. Here, we investigate the molecular details of T6P inhibition of KIN10. Molecular dynamics simulations and in vitro phosphorylation assays identified and validated the T6P binding site on KIN10. Under high-sugar conditions, T6P binds to KIN10, blocking the reorientation of its activation loop and preventing its phosphorylation and activation by GRIK1. Under these conditions, SnRK1 maintains only basal activity levels, minimizing phosphorylation of its target proteins, thereby facilitating a general shift from catabolism to anabolism.
Brookhaven National Laboratory (BNL) was awarded a pilot project in FY22 under the U.S. Department of Energy (DOE) Office of Science Biopreparedness Research Virtual Environment (BRaVE) initiative, to define research priorities, needs, and requirements for a national virtual center devoted to the biosecurity of bioenergy crops. The mission of the proposed center, referred to as the National Virtual Biosecurity for Bioenergy Crop Center (NVBBCC), would be to provide the scientific basis and tools to detect, characterize, model, and mitigate biothreats to bioenergy crops. This function will be essential to ensure the projected increased US reliance over the next few decades on key plant-based energy products, such as biojet fuel. The NVBBCC is envisioned as a distributed, virtual center with multiple national laboratories at its core to maximize the use of existing unique facilities and expertise across the DOE complex. A major goal of the pilot project was to develop a roadmap for establishing NVBBCC through a series of meetings to gather community input. A total of about 150 individuals, drawn from DOE laboratories, the USDA, academia, NIH, DHS and the private sector participated in six planning meetings held in FY23. Four of the meetings were focused on specific research topics (disease detection, dispersion and disease propagation, biomolecular characterization of plant-pathogen interaction, and mitigation strategies). These four meetings were followed by a meeting that focused on computational needs to support collaborative, data-intensive research within a distributed center as well as workforce development. A final meeting focused on establishing and maintaining preparedness within NVBBCC to respond to an emerging disease within bioenergy crops and how it would collaborate and coordinate with USDA and DHS.
Development of sustainable and scalable technologies to convert lignocellulosic biomass to biofuels is critical to achieving carbon neutrality. The potential of transgenic bioenergy crops as a renewable source of sugars and lipids has been demonstrated at bench-scale. However, scaling up these processes is important for holistic analysis. Here proof-of-concept for chemical-free hydrothermal pretreatment of transgenic energycane-oilcane line L13 at an industrially relevant scale to recover vegetative lipids along with cellulosic sugars is presented. Pilot-scale processing of 97 kg of transgenic energycane-oilcane L13 stems and high solids pretreatment of bagasse enhanced the recovery of cellulosic glucose and xylose by 5-fold as compared to untreated bagasse and helped in the enrichment of vegetative lipids in the biomass residues which allowed its recovery at the end of the bioprocess. Palmitic and oleic acids were the predominant fatty acids (FAs) extracted from stems and leaves. The processing did not affect lipid composition. The efficiency of lipid recovery from untreated biomass was 75.9% which improved to 88.7% upon pretreatment. The vegetative tissues of transgenic energycane-oilcane L13 contained 0.42 metric tons/hectare of lipids. Processing vegetative tissues yielded 0.38 metric tons/hectare of lipids. This approaches an oil yield similar to soybean (global average 0.44 metric tons/hectare) and is almost twice as high as the oil yield from sugarcane engineered to hyperaccumulate lipids (0.20 metric tons/hectare). The study suggests that further optimization by state-of-the-art metabolic engineering and biomass processing can establish transgenic bioenergy crops for commercial drop-in fuel production.
The sensor kinase Sucrose Non-fermenting-1-Related Kinase 1 (SnRK1) plays a central role in energy and metabolic homeostasis. KIN10 is a major catalytic (a) kinase subunit of SnRK1 regulated by transcription, posttranslational modification, targeted protein degradation, and its subcellular localization. Geminivirus Rep Interacting Kinase 1 and 2 (GRIK1 and 2) are immediate upstream kinases of KIN10. In the transient protein expression assays carried out in Nicotiana benthamiana (N. benthamiana) leaves, GRIK1 not only phosphorylates KIN10 but also simultaneously initiates its degradation. Posttranslational GRIK-mediated KIN10 degradation is dependent on both GRIK kinase activity and phosphorylation of the KIN10 T-loop. KIN10 proteins are significantly enriched in the grik1-1 grik2-1 double mutant, consistent with the transient assays in N. benthamiana. Interestingly. Among the enriched KIN10 proteins from grik1-1 grik2-1, is a longer isoform, putatively derived by alternative splicing which is barely detectable in wild-type plants. The reduced stability of KIN10 upon phosphorylation and activation by GRIK represents a mechanism that enables the KIN10 activity to be rapidly reduced when the levels of intracellular sugar/energy are restored to their set point, representing an important homeostatic control that prevents a metabolic overreaction to low sugar conditions. Since GRIKs are activating kinases of KIN10, KIN10s in the grik1 grik2 double null mutant background remain un-phosphorylated, with only their basal level of activity, are more stable, and therefore increase in abundance, which also explains the longer isoform KIN10L which is a minor isoform in wild type is clearly detected in the grik1 grik2 double mutant.
Development of sustainable and scalable technologies to convert lignocellulosic biomass to biofuels is critical to achieving carbon neutrality. The potential of transgenic bioenergy crops as a renewable source of sugars and lipids has been demonstrated at bench-scale. However, scaling up these processes is important for holistic analysis. Here proof-of-concept for chemical-free hydrothermal pretreatment of transgenic energycane-oilcane line L13 at an industrially relevant scale to recover vegetative lipids along with cellulosic sugars is presented. Here, pilot-scale processing of 97 kg of transgenic energycane-oilcane L13 stems and high solids pretreatment of bagasse enhanced the recovery of cellulosic glucose and xylose by 5-fold as compared to untreated bagasse and helped in the enrichment of vegetative lipids in the biomass residues which allowed its recovery at the end of the bioprocess. Palmitic and oleic acids were the predominant fatty acids (FAs) extracted from stems and leaves. The processing did not affect lipid composition. The efficiency of lipid recovery from untreated biomass was 75.9% which improved to 88.7% upon pretreatment. The vegetative tissues of transgenic energycane-oilcane L13 contained 0.42 metric tons/hectare of lipids. Processing vegetative tissues yielded 0.38 metric tons/hectare of lipids. This approaches an oil yield similar to soybean (global average 0.44 metric tons/hectare) and is almost twice as high as the oil yield from sugarcane engineered to hyperaccumulate lipids (0.20 metric tons/hectare). The study suggests that further optimization by state-of-the-art metabolic engineering and biomass processing can establish transgenic bioenergy crops for commercial drop-in fuel production.
Summary Triacylglycerols (TAG), accumulate within lipid droplets (LD), predominantly surrounded by OLEOSINs (OLE), that protect TAG from hydrolysis. We tested the hypothesis that identifying and removing degradation signals from OLE would promote its abundance, preventing TAG degradation and enhancing TAG accumulation. We tested whether mutating potential ubiquitin‐conjugation sites in a previously reported improved Sesamum indicum OLE ( SiO ) variant, o3‐3 Cys‐OLE ( SiCO herein), would stabilize it and increase its lipogenic potential. SiCOv1 was created by replacing all five lysines in SiCO with arginines. Separately, six cysteine residues within SiCO were deleted to create SiCOv2. SiCOv1 and SiCOv2 mutations were combined to create SiCOv3. Transient expression of SiCOv3 in Nicotiana benthamiana increased TAG by two‐fold relative to SiCO. Constitutive expression of SiCOv3 or SiCOv5 , containing the five predominant TAG‐increasing mutations from SiCOv3 , in Arabidopsis along with mouse DGAT2 ( mD ) increased TAG accumulation by 54% in leaves and 13% in seeds compared with control lines coexpressing SiCO and mD. Lipid synthesis rates increased, consistent with an increase in lipid sink strength that sequesters newly synthesized TAG, thereby relieving the constitutive BADC‐dependent inhibition of ACCase reported for WT Arabidopsis. These OLE variants represent novel factors for potentially increasing TAG accumulation in a variety of oil crops.
Triacylglycerols (TAG), accumulate within lipid droplets (LD), predominantly surrounded by OLEOSINs (OLE), that protect TAG from hydrolysis. We tested the hypothesis that identifying and removing degradation signals from OLE would promote its abundance, preventing TAG degradation and enhancing TAG accumulation. We tested whether mutating potential ubiquitin-conjugation sites in a previously reported improved Sesamum indicum OLE ( SiO ) variant, o3-3 Cys-OLE ( SiCO herein), would stabilize it and increase its lipogenic potential. SiCOv1 was created by replacing all five lysines in SiCO with arginines. Separately, six cysteine residues within SiCO were deleted to create SiCOv2. SiCOv1 and SiCOv2 mutations were combined to create SiCOv3. Transient expression of SiCOv3 in Nicotiana benthamiana increased TAG by two-fold relative to SiCO. Constitutive expression of SiCOv3 or SiCOv5 , containing the five predominant TAG-increasing mutations from SiCOv3 , in Arabidopsis along with mouse DGAT2 ( mD ) increased TAG accumulation by 54% in leaves and 13% in seeds compared with control lines coexpressing SiCO and mD . Lipid synthesis rates increased, consistent with an increase in lipid sink strength that sequesters newly synthesized TAG, thereby relieving the constitutive BADC-dependent inhibition of ACCase reported for WT Arabidopsis. These OLE variants represent novel factors for potentially increasing TAG accumulation in a variety of oil crops.
Duckweeds have many practical applications, for example in human nutrition, as animal feed, in the production of bioplastics or vaccines, and phytoremediation (Acosta et al., 2021). Under most conditions, they reproduce asexually which provides genetically uniform material with predictable patterns of growth that make them ideal as sentinel organisms for phytotoxicity testing (Park et al., 2021). Asexual growth also results in high biomass production which makes duckweeds promising candidates as biofuel feedstocks (Acosta et al., 2021; Liang et al., 2023). In addition, duckweed species like Lemna minor and Spirodela polyrhiza are also reemerging as model organisms in plant biology as high-quality full genome assemblies and other genomic resources become available (Chang et al., 2016; Acosta et al., 2021). We argue that duckweed species are particularly of interest for the study of primary plant metabolism. Primary metabolism concerns the part of metabolism that is directly involved in the growth and development of plants, and which tends to be highly conserved among plant species. What makes duckweeds particularly attractive is that when grown on liquid media more precise control of physiological conditions can be attained relative to growth of plants in soil. Also, due to their relatively simple anatomical structure and asexual reproduction of fronds by budding, precise characterization of the physiological state under study is possible through one simple metric, i.e., the specific growth rate (rate of dry weight increase per existing dry weight), which can be incorporated relatively easily into metabolic models. This is not possible for land plants, such as Arabidopsis, because over the course of their life cycle, they go through multiple growth stages and phases of anatomical differentiation, which are much more complex to quantify. Furthermore, duckweeds can grow on organic substrates under heterotrophic or photomixotrophic conditions that facilitate isotope tracer studies. For example, in a previous study on duckweed by one of the authors, Lemna gibba (L). was grown on glucose with a position-specific 13 C-label that can be detected and resolved by Mass Spectrometry or Nuclear Magnetic Resonance spectrometry. Using this approach, the 13 C-label was traced into biomass compounds formed from glucose, particularly isoprenoid compounds. Some of the resulting labeling patterns were in apparent disagreement with predictions based on known metabolic pathways for the biosynthesis of isopentenyl pyrophosphate, the universal building block for isoprenoids, (Lichtenthaler et al., 1997). From this data it was deduced that isoprenoid compounds such as carotenoids and isoprenoid chains of phytol and plastoquinone, synthesized in the chloroplast, are produced via a previously unreported plant metabolic pathway, now known as the methylerythitol/deoxyxylulose-5-phosphate pathway (Lichtenthaler et al., 1997).
Abstract Metabolic engineering for hyperaccumulation of lipids in vegetative tissues of high biomass crops promises a step change in oil yields for the production of advanced biofuels. Energycane is the ideal feedstock for this approach due to its exceptional biomass production and persistence under marginal conditions. Here, we evaluated metabolically engineered energycane with constitutive expression of the lipogenic factors WRINKLED 1 ( WRI 1), DIACYLGLYCEROL ACYLTRANSFERASE 1 ( DGAT 1), and OLEOSIN 1 ( OLE 1) for the accumulation of triacylglycerol (TAG), total fatty acid (TFA), and biomass under field conditions at the University of Florida‐IFAS experiment station near Citra, Florida. TAG and TFA accumulation were highest in leaves (up to 9.9% and 12.9% of DW, respectively), followed by juice from crushed stems, stems, and roots. TAG and TFA accumulation increased up to harvest time and correlated highest with OLE 1 and DGAT 1 expression. Biomass dry weight, TAG, and TFA content differed greatly depending on DGAT 1 and OLE 1 expression in transgenic lines with similar WRI 1 expression. Biomass did not significantly differ between WT and line L2 with DAGT 1 and OLE 1 expressed at low levels and TAG and TFA accumulating to 12‐ and 1.6‐fold that of WT leaves, respectively. In contrast, line L13, with intron‐mediated enhancement of DGAT 1 expression, displayed a 245‐ to 330‐fold increase in TAG and a 4.75‐ to 6.45‐fold increase in TFA content compared with WT leaves and a biomass reduction of 52%. These results provide the basis for developing novel feedstocks for expanding plant lipid production and point to new prospects for advanced biofuels.
Raw data for publication: Viet Dang Cao, Baskaran Kannan, Guangbin Luo, Hui Liu, John Shanklin, and Fredy Altpeter. 2023. Triacylglycerol, total fatty acid and biomass accumulation of metabolically engineered energycane grown under field conditions. GCB-Bioenergy (accepted for publication).
Abstract Background Metabolic engineering for hyperaccumulation of lipids in vegetative tissues is a novel strategy for enhancing energy density and biofuel production from biomass crops. Energycane is a prime feedstock for this approach due to its high biomass production and resilience under marginal conditions. DIACYLGLYCEROL ACYLTRANSFERASE ( DGAT ) catalyzes the last and only committed step in the biosynthesis of triacylglycerol (TAG) and can be a rate-limiting enzyme for the production of TAG. Results In this study, we explored the effect of intron-mediated enhancement (IME) on the expression of DGAT 1 and resulting accumulation of TAG and total fatty acid (TFA) in leaf and stem tissues of energycane. To maximize lipid accumulation these evaluations were carried out by co-expressing the lipogenic transcription factor WRINKLED 1 ( WRI 1) and the TAG protect factor oleosin ( OLE 1). Including an intron in the codon-optimized TmDGAT 1 elevated the accumulation of its transcript in leaves by seven times on average based on 5 transgenic lines for each construct. Plants with WRI 1 (W), DGAT 1 with intron (Di), and OLE 1 (O) expression (WDiO) accumulated TAG up to a 3.85% of leaf dry weight (DW), a 192-fold increase compared to non-modified energycane (WT) and a 3.8-fold increase compared to the highest accumulation under the intron-less gene combination (WDO). This corresponded to TFA accumulation of up to 8.4% of leaf dry weight, a 2.8-fold or 6.1-fold increase compared to WDO or WT, respectively. Co-expression of WDiO resulted in stem accumulations of TAG up to 1.14% of DW or TFA up to 2.08% of DW that exceeded WT by 57-fold or 12-fold and WDO more than twofold, respectively. Constitutive expression of these lipogenic “push pull and protect” factors correlated with biomass reduction. Conclusions Intron-mediated enhancement (IME) of the expression of DGAT resulted in a step change in lipid accumulation of energycane and confirmed that under our experimental conditions it is rate limiting for lipid accumulation. IME should be applied to other lipogenic factors and metabolic engineering strategies. The findings from this study may be valuable in developing a high biomass feedstock for commercial production of lipids and advanced biofuels. Graphical abstract
Metabolic engineering for hyperaccumulation of lipids in vegetative tissues is a novel strategy for enhancing energy density and biofuel production from biomass crops. Energycane is a prime feedstock for this approach due to its high biomass production and resilience under marginal conditions. DIACYLGLYCEROL ACYLTRANSFERASE ( DGAT ) catalyzes the last and only committed step in the biosynthesis of triacylglycerol (TAG) and can be a rate-limiting enzyme for the production of TAG. In this study, we explored the effect of intron-mediated enhancement (IME) on the expression of DGAT1 and resulting accumulation of TAG and total fatty acid (TFA) in leaf and stem tissues of energycane. To maximize lipid accumulation these evaluations were carried out by co-expressing the lipogenic transcription factor WRINKLED1 ( WRI1 ) and the TAG protect factor oleosin ( OLE1 ). Including an intron in the codon-optimized TmDGAT1 elevated the accumulation of its transcript in leaves by seven times on average based on 5 transgenic lines for each construct. Plants with WRI1 (W), DGAT1 with intron (Di), and OLE1 (O) expression (WDiO) accumulated TAG up to a 3.85% of leaf dry weight (DW), a 192-fold increase compared to non-modified energycane (WT) and a 3.8-fold increase compared to the highest accumulation under the intron-less gene combination (WDO). This corresponded to TFA accumulation of up to 8.4% of leaf dry weight, a 2.8-fold or 6.1-fold increase compared to WDO or WT, respectively. Co-expression of WDiO resulted in stem accumulations of TAG up to 1.14% of DW or TFA up to 2.08% of DW that exceeded WT by 57-fold or 12-fold and WDO more than twofold, respectively. Constitutive expression of these lipogenic “push pull and protect” factors correlated with biomass reduction. Intron-mediated enhancement (IME) of the expression of DGAT resulted in a step change in lipid accumulation of energycane and confirmed that under our experimental conditions it is rate limiting for lipid accumulation. IME should be applied to other lipogenic factors and metabolic engineering strategies. The findings from this study may be valuable in developing a high biomass feedstock for commercial production of lipids and advanced biofuels.
Abstract Background Duckweeds are small, rapidly growing aquatic flowering plants. Due to their ability for biomass production at high rates they represent promising candidates for biofuel feedstocks. Duckweeds are also excellent model organisms because they can be maintained in well-defined liquid media, usually reproduce asexually, and because genomic resources are becoming increasingly available. To demonstrate the utility of duckweed for integrated metabolic studies, we examined the metabolic adaptation of growing Lemna gibba cultures to different nutritional conditions. Results To establish a framework for quantitative metabolic research in duckweeds we derived a central carbon metabolism network model of Lemna gibba based on its draft genome. Lemna gibba fronds were grown with nitrate or glutamine as nitrogen source. The two conditions were compared by quantification of growth kinetics, metabolite levels, transcript abundance, as well as by 13 C-metabolic flux analysis. While growing with glutamine, the fronds grew 1.4 times faster and accumulated more protein and less cell wall components compared to plants grown on nitrate. Characterization of photomixotrophic growth by 13 C-metabolic flux analysis showed that, under both metabolic growth conditions, the Calvin-Benson-Bassham cycle and the oxidative pentose-phosphate pathway are highly active, creating a futile cycle with net ATP consumption. Depending on the nitrogen source, substantial reorganization of fluxes around the tricarboxylic acid cycle took place, leading to differential formation of the biosynthetic precursors of the Asp and Gln families of proteinogenic amino acids. Despite the substantial reorganization of fluxes around the tricarboxylic acid cycle, flux changes could largely not be associated with changes in transcripts. Conclusions Through integrated analysis of growth rate, biomass composition, metabolite levels, and metabolic flux, we show that Lemna gibba is an excellent system for quantitative metabolic studies in plants. Our study showed that Lemna gibba adjusts to different nitrogen sources by reorganizing central metabolism. The observed disconnect between gene expression regulation and metabolism underscores the importance of metabolic flux analysis as a tool in such studies.
Mössbauer spectroscopy provides significant insights into the electronic structure and environment of the metal centers. In this work, we investigate the electronic structures of a set of nonheme diiron complexes by evaluating two key parameters pertaining to Mössbauer spectroscopy, namely, the isomer shift (δ) and quadrupole splitting (|ΔE Q |), using different levels of density functional theory (DFT). The diiron systems investigated here span diverse oxidation states, bridging motifs, and spin coupling patterns, which present a challenging case for theoretical predictions. We demonstrate that the combination of B97-D3/def2-TZVP is an efficient approach in modeling both the δ and |ΔEQ| values with high accuracy for the representative nonheme diiron complexes. We also show that δ is accurately predicted irrespective of the choice of approximate density functional while the |ΔE Q | is sensitive to the level of theory employed. Further investigation shows that the present methodology assessed using synthetic nonheme diiron complexes could be extended to nonheme diiron enzyme active sites, featuring both ferromagnetic and antiferromagnetic coupling between the iron centers.