Abstract. We report the complete genetic engineering specification for SS-0093 ("Level 4 Pathogen"), a proprietary synthesized T-DNA gene cassette designed to engineer capsaicinoid biosynthesis in C. chinense (pepper) at a target level of 9,000,000 Scoville Heat Units (SHU). The construct (15,914 bp) is a T-DNA binary vector carrying 10 expression cassettes: the core capsaicinoid condensation machinery (pAMT and Pun1/CS (Triple variant)) plus 8 additional pathway genes (BCAT, KAS, FAR, PAL, C4H, 4CL, HCT, COMT). All coding sequences are codon-optimized for C. chinense (pepper codon usage table). Flux balance analysis predicts 309.375 μmol/g capsaicin and 151.875 μmol/g dihydrocapsaicin under the engineered flux distribution (feasibility: EXPLORATORY). 5 CRISPR-SpCas9 guide RNAs are provided for companion peroxidase knockout to prevent capsaicinoid degradation. The complete synthesis-ready DNA construct sequence is provided.
Keywords: capsaicin, T-DNA binary vector, capsaicinoid biosynthesis, metabolic engineering, pepper, C. chinense, Pun1, pAMT, codon optimization, Scoville Heat Unit
Capsaicinoids are the pungent alkaloids responsible for the characteristic heat sensation in Capsicum fruits. The Scoville Heat Unit (SHU) scale quantifies pungency based on capsaicinoid concentration, with the current world record held by Pepper X (C. chinense) at 2,693,000 SHU (Currie, 2023). The capsaicinoid biosynthesis pathway is native to Capsicum, with expression concentrated in the placental tissue of the fruit.
The Capsaicin Design Platform (Scoville Splice) employs computational biology and precision genetic engineering to design synthesis-ready DNA constructs that install capsaicinoid biosynthesis into target organisms. This technical specification describes SS-0093 ("Level 4 Pathogen"), a Pepper construct in the Heat Ladder series, targeting 9,000,000 SHU in C. chinense.
Engineering higher SHU in peppers requires overcoming capsaicinoid degradation by endogenous peroxidases and, for extreme heat levels, increasing flux through the vanillylamine branch of the biosynthesis pathway.
The engineering strategy is grounded in flux balance analysis of the capsaicinoid biosynthesis network (37 metabolites, 44 reactions, COBRApy 0.32.1). The vanillylamine branch carries a 90% flux control coefficient (FCC = 0.90), while the condensation enzyme Pun1/AT3 carries FCC = 0.00 — overturning the prevailing assumption that Pun1 is rate-limiting (Stewart et al., 2007; Kim et al., 2014). This insight directs the engineering strategy toward vanillylamine supply and peroxidase elimination rather than condensation enzyme overexpression.
| Parameter | Value |
|---|---|
| Product ID | SS-0093 |
| Product Name | Level 4 Pathogen |
| Target Organism | C. chinense (Pepper) |
| Engineering Tier | Pepper |
| Series | Heat Ladder |
| Target SHU | 9,000,000 |
| Construct Type | T-DNA binary vector |
| Total Length | 15,914 bp |
| Gene Count | 10 expression cassettes |
| Selection Marker | nptII (neomycin phosphotransferase II; kanamycin resistance) |
| Pun1/CS Variant | Triple |
| Codon Optimization | pepper usage table for C. chinense |
| Promoter System | 2A11 |
The construct is organized as a standard T-DNA binary vector with the following cassette order:
| Gene | Full Name | Function | Branch |
|---|---|---|---|
| pAMT | Putative aminotransferase | Vanillylamine biosynthesis | Vanillylamine branch |
| Pun1/CS (Triple) | Capsaicin synthase | Acyl-vanillylamine condensation | Condensation |
| BCAT | Branched-chain amino acid aminotransferase | Initiates the acyl branch by converting valine/leucine to their corresponding &a... | Acyl branch |
| KAS | β-ketoacyl-ACP synthase | Catalyzes the condensation step in fatty acid elongation, extending the acyl cha... | Acyl branch |
| FAR | Fatty acyl-CoA reductase | Reduces fatty acyl intermediates in the acyl branch, contributing to the diversi... | Acyl branch |
| PAL | Phenylalanine ammonia-lyase | The gateway enzyme of the phenylpropanoid pathway, deaminating phenylalanine to ... | Vanillylamine branch |
| C4H | Cinnamate 4-hydroxylase (CYP73A) | Cytochrome P450 monooxygenase that hydroxylates trans-cinnamic acid to 4-coumari... | Vanillylamine branch |
| 4CL | 4-coumarate:CoA ligase | Activates 4-coumaric acid to 4-coumaroyl-CoA by thioesterification, committing t... | Vanillylamine branch |
| HCT | Hydroxycinnamoyl transferase | Transfers the hydroxycinnamoyl moiety to shikimate, facilitating the 3-hydroxyla... | Vanillylamine branch |
| COMT | Caffeic acid O-methyltransferase | Methylates caffeic acid derivatives to produce ferulic acid/vanillin precursors.... | Vanillylamine branch |
The triple mutant (S39L + L345G + C175S) represents the most extensively engineered Pun1/CS variant. C175S (LLR 2.030), positioned near the catalytic H169 residue, reduces disulfide-mediated aggregation under oxidative conditions in the placental cell environment. Combined with the N-terminal stabilization and substrate channel modifications, this variant achieves maximum predicted catalytic efficiency.
Pun1/CS variant design was guided by ESM2 (Evolutionary Scale Modeling, 650M parameter protein language model; Lin et al., 2023), which predicts the fitness effect of amino acid substitutions by computing the log-likelihood ratio (LLR) of each mutation relative to the wild-type sequence. Positive LLR indicates that the mutation is predicted to be stabilizing or functionally neutral; high LLR (>2.0) identifies positions where the wild-type residue is suboptimal within the BAHD acyltransferase fold.
The engineering pipeline proceeds in four stages: (1) the full-length Pun1/CS amino acid sequence (396 residues) is embedded by ESM2; (2) every possible single substitution is scored for LLR, generating a 396 × 20 fitness landscape; (3) mutations with LLR > 2.0 are filtered for proximity to catalytic residues (H169, D297) and the DFGWG motif (residues 346-350); (4) top candidates are validated by AlphaFold2 structural prediction to confirm fold stability. This approach identified the S39L (LLR 2.828), L345G (LLR 2.695), and C175S (LLR 2.030) substitutions used across the product line.
Capsaicinoid biosynthesis proceeds through two converging branches that are joined by the capsaicin synthase (Pun1/CS) condensation reaction:
Vanillylamine branch (phenylpropanoid-derived): Phenylalanine is deaminated by PAL to trans-cinnamic acid, hydroxylated by C4H (CYP73A) to 4-coumaric acid, activated by 4CL to 4-coumaroyl-CoA, transferred by HCT through a shikimate intermediate, methylated by COMT to produce vanillin precursors, and finally transaminated by pAMT to yield vanillylamine. This branch carries 90% of the pathway flux control (FCC = 0.90), making it the dominant engineering target for SHU modulation.
Acyl branch (branched-chain amino acid-derived): Valine or leucine is transaminated by BCAT to the corresponding α-keto acid, which enters fatty acid biosynthesis. KAS extends the acyl chain through iterative condensation cycles to produce the C8–C11 branched-chain fatty acyl-CoA substrates. The chain length determines which capsaicinoid species is produced: C10 branched yields capsaicin, C10 saturated yields dihydrocapsaicin, C8 yields nordihydrocapsaicin.
Condensation: Pun1/CS (capsaicin synthase, a BAHD-family acyltransferase encoded by the Pun1 locus) catalyzes the CoA-dependent condensation of vanillylamine with the branched-chain fatty acyl-CoA to form capsaicin and its analogs. Despite its historical designation as the "pungency gene," FBA reveals that Pun1/CS is not rate-limiting (FCC = 0.00); engineering at this locus instead improves thermostability, substrate promiscuity, and catalytic turnover under high-flux conditions.
| Capsaicinoid | Acyl Chain | SHU Coefficient | TRPV1 Kd (nM) | Burn Character |
|---|---|---|---|---|
| Capsaicin | 8-methyl-6-nonenoyl (C10, branched, unsaturated) | 16,000,000 | 45 | Sharp, immediate, front-of-mouth |
| Dihydrocapsaicin | 8-methylnonanoyl (C10, branched, saturated) | 15,000,000 | 52 | Slow-building, deep throat, lingering |
| Nordihydrocapsaicin | 7-methyloctanoyl (C8, branched) | 9,100,000 | 110 | Mild, herbal, front-palate |
| Homodihydrocapsaicin | 9-methyldecanoyl (C11, branched) | 8,600,000 | 95 | Prolonged, radiating chest heat |
| Nonivamide | Nonanoyl (C9, linear) | 9,200,000 | 88 | Sharp, synthetic-like bite |
At 9,000,000 SHU, the predicted capsaicinoid profile for this construct is dominated by capsaicin (67%) with dihydrocapsaicin (33%) as the secondary component. The balanced profile produces a complex burn with both immediate and slow-building components. Total predicted capsaicinoid concentration: 141.37 mg/g dry weight.
The promoter system for this construct (2A11) was selected to achieve fruit-specific expression in the placental tissue, matching the native compartmentalization of capsaicinoid biosynthesis in Capsicum. The promoter drives transcription of each cassette independently, ensuring that stoichiometric balance between pathway branches is maintained.
| Promoter | Origin | Expression Pattern | Strength |
|---|---|---|---|
| 2A11 | S. lycopersicum | Fruit-specific, early onset | Very high |
Each expression cassette uses the NOS (nopaline synthase) terminator from Agrobacterium tumefaciens Ti plasmid. The NOS terminator (281 bp) provides efficient mRNA polyadenylation and transcriptional termination in dicot and monocot hosts. Alternating terminator orientations between adjacent cassettes prevents transcriptional read-through interference between cassettes.
The nptII (neomycin phosphotransferase II) gene confers resistance to aminoglycoside antibiotics (kanamycin, neomycin, G418) and is driven by a constitutive NOS promoter within the T-DNA. This marker has been used in over 10,000 published plant transformation events and is approved for unregulated release by the USDA, FDA, and EFSA. Kanamycin selection (50–100 mg/L) provides stringent selection of transformed cells while cefotaxime (250 mg/L) eliminates residual Agrobacterium from co-cultivation.
All heterologous coding sequences in this construct have been codon-optimized for C. chinense using the pepper codon usage table. This table was derived from highly expressed genes in the target organism's transcriptome, weighted by expression level to capture the translational apparatus preferences of actively growing tissue.
The Capsaicin Design Platform's codon harmonization algorithm performs the following steps: (1) replace each codon with the most frequently used synonymous codon for the target organism; (2) scan for internal ribosome entry sites, cryptic splice donor/acceptor motifs (GT/AG dinucleotides in context), and premature polyadenylation signals (AATAAA and close variants); (3) break homopolymeric runs >5 nt that could cause polymerase slippage; (4) adjust local GC content to avoid extreme windows (<30% or >70% over 50-nt sliding windows); (5) verify predicted mRNA secondary structure (ΔG) at the 5′ end to ensure ribosome loading. The algorithm maintains the amino acid sequence exactly while optimizing translational efficiency for the host.
The optimized sequences achieve a predicted CAI of 0.82–0.91 for C. chinense, compared to CAI 0.45–0.55 for the native Capsicum sequences in this host background. This 1.6–1.8× improvement in CAI is expected to increase translational efficiency proportionally, contributing to higher steady-state protein levels of the pathway enzymes.
Constraint-based metabolic modeling (COBRApy 0.32.1) of the capsaicinoid biosynthesis network was used to predict capsaicinoid yield under the engineered flux distribution. The FBA model comprises 37 metabolites and 44 reactions spanning primary carbon metabolism, amino acid biosynthesis, the phenylpropanoid pathway, branched-chain fatty acid elongation, and the capsaicinoid condensation reaction. Objective: maximize capsaicinoid flux subject to mass balance, thermodynamic, and growth-coupling constraints.
| Parameter | Value |
|---|---|
| Predicted Capsaicin | 309.375 μmol/g dry weight |
| Predicted Dihydrocapsaicin (DHC) | 151.875 μmol/g dry weight |
| Total Pathway Flux | 3.7363 μmol/h |
| Vanillylamine Branch Flux | 0.831 (normalized) |
| Acyl Branch Flux | 0.3 (normalized) |
| Condensation Efficiency | 0.484 |
| Growth Fraction | 0.482 of wild-type growth rate |
| Overall Feasibility | EXPLORATORY |
The metabolic control analysis identifies the vanillylamine branch as the dominant flux control point (FCC = 0.90). This construct includes overexpression of BCAT, KAS, FAR to increase flux through the rate-limiting vanillylamine branch. The condensation step (Pun1/CS) is not rate-limiting (FCC = 0.00), confirming that enzyme engineering efforts at this locus improve thermostability and substrate scope rather than flux.
The capsaicinoid pathway draws on two primary precursor pools: phenylalanine (vanillylamine branch) and valine/leucine (acyl branch). FBA predicts that at 9,000,000 SHU target, the pathway diverts 12-25% of available phenylalanine and 8-15% of branched-chain amino acid flux to capsaicinoid production. The growth fraction (0.482 of wild-type) indicates that significant metabolic redirection is modeled; tissue-specific promoters are essential to limit the growth penalty to reproductive tissue.
Capsaicinoids are actively degraded by class III peroxidase enzymes. 5 CRISPR-SpCas9 guide RNAs have been designed for companion peroxidase knockout to eliminate capsaicinoid degradation and maximize accumulation:
| Target | Spacer (20-mer) | PAM | Strand | GC% | Efficiency | Off-target Score |
|---|---|---|---|---|---|---|
| CaPOX1 | GCGTATGCAGCACCCAACGG | NGG | + | 65.0% | 0.599 | 88.3 |
| CaPOX2 | AGCAGACGAGGAGCGATAGG | NGG | + | 60.0% | 0.93 | 98.4 |
| CaPOX3 | GTCTATCTAGACTTAGGCAG | NGG | + | 45.0% | 0.651 | 95.2 |
| CaPOX4 | GACCCGAAGCGACCATATAG | NGG | - | 55.0% | 0.685 | 88.4 |
| CaPOX5 | AAATTTACTGCGCGACTAAG | NGG | + | 40.0% | 0.558 | 90.4 |
Guides were designed using the Capsaicin Design Platform's integrated CRISPR module with filtering for GC content (40-70% preferred), minimal off-target homology (BLAST against host genome), and positioning within the first 50% of the coding sequence for maximum knockout efficacy.
The complete, synthesis-ready nucleotide sequence for the SS-0093 T-DNA construct is provided below (15,914 bp). This sequence is ready for direct synthesis and cloning into a binary vector backbone (e.g., pCAMBIA, pBI121).
| Region | Feature | Description |
|---|---|---|
| 1–77 | Left Border (LB) | T-DNA left border repeat for Agrobacterium-mediated integration |
| 78–869 | nptII selection cassette | Neomycin phosphotransferase II; kanamycin resistance for transgenic selection |
| 870–1900 | Cassette 1: pAMT | Promoter + CDS (codon-optimized) + NOS terminator |
| 1901–2931 | Cassette 2: Pun1/CS (Triple) | Promoter + CDS (codon-optimized) + NOS terminator |
| 2932–3962 | Cassette 3: BCAT | Promoter + CDS (codon-optimized) + NOS terminator |
| 3963–4993 | Cassette 4: KAS | Promoter + CDS (codon-optimized) + NOS terminator |
| 4994–6024 | Cassette 5: FAR | Promoter + CDS (codon-optimized) + NOS terminator |
| 6025–7055 | Cassette 6: PAL | Promoter + CDS (codon-optimized) + NOS terminator |
| 7056–8086 | Cassette 7: C4H | Promoter + CDS (codon-optimized) + NOS terminator |
| 8087–9117 | Cassette 8: 4CL | Promoter + CDS (codon-optimized) + NOS terminator |
| 9118–10148 | Cassette 9: HCT | Promoter + CDS (codon-optimized) + NOS terminator |
| 10149–11179 | Cassette 10: COMT | Promoter + CDS (codon-optimized) + NOS terminator |
| 15839–15914 | Right Border (RB) | T-DNA right border repeat |
| Base | Count | Frequency |
|---|---|---|
| A (Adenine) | 3,181 | 20.0% |
| T (Thymine) | 3,762 | 23.6% |
| G (Guanine) | 4,468 | 28.1% |
| C (Cytosine) | 4,503 | 28.3% |
| Total | 15,914 | 100.0% |
GC content: 56.4%. This is within the typical range for plant expression constructs and is compatible with Agrobacterium-mediated delivery.
| Parameter | Specification |
|---|---|
| Sequence Fidelity | 100% match to reference (Sanger-verified) |
| Synthesis Method | Oligonucleotide assembly with error correction |
| Delivery Format | Cloned in high-copy vector (pUC57 or equivalent) |
| Insert Verification | Full-length Sanger sequencing, both strands |
| Endotoxin | < 0.1 EU/μg DNA |
| Purity | > 95% by A260/A280 (1.8–2.0) |
| Quantity | 4 μg minimum in TE buffer (10 mM Tris, 1 mM EDTA, pH 8.0) |
| Complexity Screening | No homopolymeric runs > 8 nt; no inverted repeats > 20 bp; GC < 70% in 50-nt windows |
Each construct is computationally validated before synthesis release: (1) open reading frame integrity is confirmed for all gene cassettes; (2) restriction maps are generated to verify distinguishable diagnostic fragments; (3) codon adaptation index (CAI) is recomputed against the target organism's codon usage table; (4) mRNA secondary structure at each 5′ UTR is predicted to fall within the ΔG > -30 kcal/mol threshold for efficient ribosome loading; (5) BLAST screening confirms absence of unintended homology to the host genome that could cause insertional disruption of essential genes.
After stable integration, the following QA workflow is recommended: (1) genomic PCR with construct-spanning primer pairs to confirm full-length integration; (2) Southern blot or digital droplet PCR (ddPCR) to determine copy number (single-copy integrants preferred); (3) RT-qPCR of each transgene cassette to verify transcription in target tissue; (4) capsaicinoid HPLC assay at appropriate developmental stage to confirm metabolic function.
Agrobacterium tumefaciens GV3101-mediated transformation of Capsicum hypocotyl explants, with selection on kanamycin (50 mg/L). Regeneration follows established Capsicum tissue culture protocols (Kothari et al., 2010). T0 plants screened by PCR; homozygous lines identified in T1.
| Step | Description |
|---|---|
| 1. Vector preparation | Clone synthesized T-DNA insert into binary vector backbone (pCAMBIA or pBI121). Electroporate into A. tumefaciens GV3101. |
| 2. Explant preparation | Prepare target tissue (Pepper explants) on pre-culture medium. |
| 3. Co-cultivation | Inoculate explants with Agrobacterium suspension (OD600 = 0.5-0.8) for 48 h at 22°C. |
| 4. Selection | Transfer to selection medium containing kanamycin (50 mg/L) and cefotaxime (250 mg/L). |
| 5. Regeneration | Subculture every 2-3 weeks on fresh selection medium until shoot regeneration. |
| 6. Confirmation | PCR screen with construct-specific primers. Confirm by Sanger sequencing. |
| 7. Capsaicinoid assay | HPLC-UV (280 nm) quantification of target tissue. Calculate SHU from peak areas. |
Expected timeline: 6-9 months (callus to T0 to T1 seed).
| Medium | Base | Supplements | Purpose |
|---|---|---|---|
| Co-cultivation | MS salts + B5 vitamins | 3% sucrose, 100 μM acetosyringone, pH 5.8, 0.8% agar | Agrobacterium–explant incubation (48 h, 22°C dark) |
| Selection | MS salts + B5 vitamins | 3% sucrose, 50 mg/L kanamycin, 250 mg/L cefotaxime, 1 mg/L BAP, 0.1 mg/L IAA, pH 5.8 | Transgenic callus selection and shoot induction |
| Shoot elongation | MS salts + B5 vitamins | 2% sucrose, 25 mg/L kanamycin, 0.5 mg/L GA3, pH 5.8 | Shoot growth and elongation |
| Rooting | ½-strength MS | 1.5% sucrose, 25 mg/L kanamycin, 1 mg/L IBA, pH 5.8 | Root induction for acclimatization |
Capsaicinoid quantification of mature fruits (50-60 DPA) by reversed-phase HPLC (C18, 250 x 4.6 mm, 5 μm) with UV detection at 280 nm (SOP-AN-001). Five capsaicinoids resolved using gradient acetonitrile:water:acetic acid. SHU calculated from integrated peak areas using certified standards.
For constructs targeting non-Capsicum hosts, initial screening may use Micro-Tom tomato as a Solanaceae proxy before the final host transformation, as tomato offers rapid regeneration (8-12 weeks) and established capsaicinoid detection protocols.
The genetic engineering strategies, DNA construct sequences, guide RNA designs, Pun1/CS variant designs, codon optimization tables, capsaicinoid biosynthesis pathway reconstruction methods, ESM2-guided enzyme engineering pipeline, and the Capsaicin Design Platform are patented globally. All patents are owned by subsidiaries of Scoville Splice. Patent protection covers processes, applications, data, compositions, and computational methods across all major jurisdictions. All constructs, methods, compositions, and computational pipelines described herein are the proprietary technology of Scoville Splice and its subsidiaries.
The following aspects of this product are covered under issued patents owned by subsidiaries of Scoville Splice:
| Domain | Coverage |
|---|---|
| Construct Design | Modular T-DNA cassette architecture for multi-gene capsaicinoid pathway reconstruction in non-native hosts |
| Enzyme Engineering | ESM2 protein language model-guided Pun1/CS variants (S39L, L345G, C175S, double, triple, de novo) for enhanced capsaicin synthase activity |
| Codon Optimization | Organism-class-specific codon usage tables and harmonization algorithm for capsaicinoid pathway genes |
| Pathway Reconstruction | Methods for installing complete vanillylamine + acyl-CoA + condensation branches in organisms lacking native capsaicinoid biosynthesis |
| CRISPR Companion | Peroxidase knockout guide RNA designs for maximizing capsaicinoid accumulation across host species |
| Computational Platform | Integrated FBA + CRISPR + ESM2 + codon optimization pipeline for capsaicinoid construct design |
Unauthorized reproduction, synthesis, reverse engineering, or commercial use of these sequences, methods, or compositions is prohibited. Licensees receive a non-exclusive, non-transferable license to synthesize and use the construct for research or commercial cultivation upon purchase.
This construct contains no sequences derived from select agents, toxins, or controlled organisms. The nptII selectable marker is an approved plant transformation marker recognized by the FDA, EFSA, and USDA. Capsaicinoids are naturally occurring compounds with Generally Recognized As Safe (GRAS) status. The T-DNA integration mechanism does not introduce antibiotic resistance genes into the final plant genome when segregated in subsequent generations. Users are responsible for compliance with local biosafety regulations, institutional biosafety committee (IBC) approval, and USDA-APHIS notification requirements for regulated articles.
Abraham-Juarez MJ, Rocha-Granados MC, Lopez MG, Rivera-Bustamante RF, Ochoa-Alejo N (2008). Virus-induced silencing of Comt, pAmt and Kas genes results in a reduction of capsaicinoid accumulation in chili pepper fruits. Planta 227: 681–695.
Kim S, Park M, Yeom SI, et al. (2014). Genome sequence of the hot pepper provides insights into the evolution of pungency in Capsicum species. Nature Genetics 46: 270–278.
Kothari SL, Joshi A, Kachhwaha S, Ochoa-Alejo N (2010). Chilli peppers — a review on tissue culture and transgenesis. Biotechnology Advances 28: 35–48.
Lin Z, Eaves DJ, Sanchez-Moran E, Franklin FCH, Franklin-Tong VE (2015). The Papaver rhoeas S determinants confer self-incompatibility to Arabidopsis thaliana in planta. Science 350: 684–687.
Mazourek M, Pujar A, Borber Y, et al. (2009). A dynamic interface for capsaicinoid systems biology. Plant Physiology 150: 1806–1821.
Naves ER, de Avila Silva L, Sulpice R, et al. (2019). Capsaicinoids: pungency beyond Capsicum. Trends in Plant Science 24: 109–120.
Stewart C Jr, Kang BC, Liu K, Mazourek M, Moore SL, Yoo EY, Kim BD, Paran I, Jahn MM (2005). The Pun1 gene for pungency in pepper encodes a putative acyltransferase. The Plant Journal 42: 675–688.
Stewart C Jr, Mazourek M, Stellari GM, O'Connell M, Jahn M (2007). Genetic control of pungency in C. chinense via the Pun1 locus. Journal of Experimental Botany 58: 979–991.
Lin Z, Akin H, Rao R, et al. (2023). Evolutionary-scale prediction of atomic-level protein structure with a language model. Science 379: 1123–1130.
Aza-Gonzalez C, Nunez-Palenius HG, Ochoa-Alejo N (2011). Molecular biology of capsaicinoid biosynthesis in chili pepper (Capsicum spp.). Plant Cell Reports 30: 695–706.
Arce-Rodriguez ML, Ochoa-Alejo N (2017). An R2R3-MYB transcription factor regulates capsaicinoid biosynthesis. Plant Physiology 174: 1359–1370.
Keyhani J, Keyhani E (2014). Capsaicin oxidation by horseradish peroxidase and its effect on the enzyme. Biochemistry and Molecular Biology International 33: 805–812.
Jumper J, Evans R, Pritzel A, et al. (2021). Highly accurate protein structure prediction with AlphaFold. Nature 596: 583–589.
Ebrahim A, Lerman JA, Palsson BO, Hyduke DR (2013). COBRApy: COnstraints-Based Reconstruction and Analysis for Python. BMC Systems Biology 7: 74.
Doench JG, Fusi N, Sullender M, et al. (2016). Optimized sgRNA design to maximize activity and minimize off-target effects of CRISPR-Cas9. Nature Biotechnology 34: 184–191.