A synthetic biology approach to self-regulatory recombinant protein production in Escherichia coli
© Dragosits et al; licensee BioMed Central Ltd. 2012
Received: 23 December 2011
Accepted: 30 March 2012
Published: 1 December 2012
Recombinant protein production is a process of great industrial interest, with products that range from pharmaceuticals to biofuels. Since high level production of recombinant protein imposes significant stress in the host organism, several methods have been developed over the years to optimize protein production. So far, these trial-and-error techniques have proved laborious and sensitive to process parameters, while there has been no attempt to address the problem by applying Synthetic Biology principles and methods, such as integration of standardized parts in novel synthetic circuits.
We present a novel self-regulatory protein production system that couples the control of recombinant protein production with a stress-induced, negative feedback mechanism. The synthetic circuit allows the down-regulation of recombinant protein expression through a stress-induced promoter. We used E. coli as the host organism, since it is widely used in recombinant processes. Our results show that the introduction of the self-regulatory circuit increases the soluble/insoluble ratio of recombinant protein at the expense of total protein yield. To further elucidate the dynamics of the system, we developed a computational model that is in agreement with the observed experimental data, and provides insight on the interplay between protein solubility and yield.
Our work introduces the idea of a self-regulatory circuit for recombinant protein products, and paves the way for processes with reduced external control or monitoring needs. It demonstrates that the library of standard biological parts serves as a valuable resource for initial synthetic blocks that needs to be further refined to be successfully applied in practical problems of biotechnological significance. Finally, the development of a predictive model in conjunction with experimental validation facilitates a better understanding of the underlying dynamics and can be used as a guide to experimental design.
KeywordsSynthetic biology Recombinant protein Self-regulatory Escherichia coli Stress promoter
Recombinant or heterologous protein production (RPP) is an important biotechnological process, with applications that range from catalysis (e.g. washing detergents) and therapeutic use (e.g. antibody production), to protein production for enzymatic characterization and crystallography. Production of human proteins in bacteria dates back to the production of the 14-codon somatostatin gene in Escherichia coli in 1977 . Since then, several hosts have been explored, including other prokaryote species , various yeast and fungal species , plant, insect, and mammalian cell lines . As there is no universally optimal host, the choice of host is based on various parameters (protein yield, production time, etc.) on a case-by-case basis.
One of the most critical parameters, especially for proteins of therapeutic interest, is the presence of post-translational modifications. Complex proteins might harbor disulfide bonds as well as complex glycan structures (e.g. antibodies and antibody fragments) that influence the 3D structure, serum stability and the protein effector functions . However, in the case of E. coli, many engineered strains and expression platforms were made available over the years , including strains that enable some complex post-translational modifications . These advances, together with its easy cultivation, fast growth, and well-studied physiology explain E. coli 's role as a major host for RPP, including sensitive applications such as therapeutic protein production .
Over-expression of recombinant proteins can lead to significant stress in the host cell, which in turn limits its capacity to function as a cell factory. First, it constitutes a general metabolic burden, as it is responsible for depletion of precursor metabolites . Recombinant proteins are usually produced in very high amounts (it is not unusual to comprise 30% or more of total cellular protein in the cell), which leads to a significant stress for the cell. The latter reacts by employing a heat-shock like response, which involves the induction of chaperones, foldases and proteases . Another potential drawback of RPP in bacterial cells is that many recombinant proteins form inclusion bodies (IBs) which represent insoluble protein aggregates and necessitates elaborate downstream processing including de/re-naturation methods . Recent work shows that IB-trapped proteins may actually be used directly , a result that contradicts the traditional thinking that IBs consist of misfolded, and thus inactive, protein. Still, in both cases further processing is needed to achieve soluble protein products . On the other hand, the formation of inclusion bodies can also represent a favorable factor as the formation of such insoluble protein greatly facilitates initial protein enrichment. As such, the protein of interest as well as the purpose of the recombinant product will determine the desirable approach.
Common techniques to optimize recombinant protein production in bacteria
Natural and engineered host strains can accommodate higher recombinant protein yields
Plasmid copy number
The choice of the plasmid backbone influences the production process through gene dosage
Inducer concentration influences transcription rate and therefore product formation/aggregation rate.
Different promoters can be considered. Weak/Strong, constitutive and inducible promoters.
Ribosome binding site (RBS)
Position and sequence of the RBS influences translational efficiency
mRNA stability and structure
mRNA turnover influences the production process as well as mRNA structure can influence ribosome binding and translational efficiency
Codon usage in the sequence of the recombinant gene greatly impacts translation efficiency
Temperature, oxygenation, pH and medium osmolarity impact on the production process
Optimization of the growth medium can lead to increase of the product yield and decrease of by-product formation
Heat shock protein co-overexpression and knockouts
Increased or decreased amount of several molecular chaperones, foldases and proteases influence protein yield and quality.
In order to address challenges such as balancing protein production and cellular stress, we developed a synthetic expression platform that enables the cell to shut down the RPP mechanism by itself, once stress signals are detected. For this implementation we used and created new standardized parts and developed a computational model to elucidate the dynamics of this protein expression system. This study illustrates the potential of synthetic biology to help traditional biotechnological fields by constructing customized circuits with desired behaviors from standardized parts.
Results and Discussion
Overview and parts selection
As shown in Figure 1B, we selected the pET expression system in combination with the E. coli C41 (DE3) strain, both widely used for laboratory scale protein expression . The C41 strain encodes the T7 polymerase under the control of the lac promoter, and we cloned the recombinant protein downstream of a strongly-regulated T7 promoter. Addition of IPTG (Isopropyl β-D-1-thiogalactopyranoside) to the medium triggers the expression of T7 polymerase, which in turn transcribes the recombinant gene. Next, we identified a stress-sensitive promoter that is induced in recombinant protein production settings based on available literature. The IbpAB operon, which encodes inclusion-body binding proteins A and B , is known to be significantly up-regulated during the expression of various recombinant products [9, 18–21], so we selected its promoter for the expression of a repressor protein. The TetR protein  was chosen as a repressor here due to the tight repression that it confers. To achieve transient repression for a relatively short time, we used a TetR protein with a degradation tag fusion  to decrease TetR protein half-life. We considered the application of repressor variant tagged for degradation necessary because we aimed at a relatively fast repressor turn-over. Previous research on synthetic gene circuits highlight, that otherwise such synthetic circuits suffer from a long response time [24, 25], which may render them unsuitable for recombinant protein production processes. Finally, green fluorescence protein (GFP) was used as the model recombinant protein due to its high yields and capacity for rapid, inexpensive screening. Furthermore, the GFP mutant used in this study is known to form inclusion bodies, even when using weaker expression systems , which makes it a suitable model for evaluating the effect of a negative feedback circuit on protein production and solubility. For GFP expression experiments, saturating IPTG concentrations (1 mM) were applied.
Integration of the Tet operator and stress promoter
Mutant libraries and dynamic range
The incorporation of the WT P IbpAB promoter led to high basal levels of promoter activity, significant decrease of intracellular GFP levels, and high variability in GFP expression. To increase the dynamic range of GFP expression that the circuit can operate and to investigate whether the variability and basal expression can be further reduced across cells, we created a P IbpAB mutant promoter library by using error-prone PCR that was characterized by using RFP as a reporter of promoter strength. A plasmid already containing the RFP gene was used to produce this mutant reporter library. Approximately 100 single clones of this library were analyzed with respect to RFP protein production and clones with reduced protein production with respect to the wild type P IbpAB promoter sequence were found. Sequencing elucidated the genetic basis of promoter variation, showing that 3 point mutations resulted in reduced expression in the 2 mutant promoters that were analyzed further. Furthermore, close proximity of the 3 nucleotide exchanges to the -35 region (promoter mutant m4_5) lead to stronger reduction of promoter activity (Additional file 1: Figure S2). We used RFP reporter constructs as described above in order to verify that the mutant promoters, despite lower basic activity, were still activated by recombinant protein stress which is typical for the IbpAB promoter (Figure 3C and 3D). Finally, those mutants were used to replace the original wild-type P IbpAB sequence in pNF_TetR. The application of these mutant promoter sequences lead to higher GFP levels and a lower clone-to-clone variability as compared with the negative feedback production system that uses the wildtype stress promoter for repressor production (Figure 3C).
For even higher flexibility, we engineered the ribosome binding site (RBS) to create an RBS library with mutants that have different (lower) levels of translational efficiency (Additional file 1: Table S3). Using lower translational efficiency RBS sequences in combination with promoter mutants for repressor expression, we were able to create a protein expression system that achieves comparable protein production levels with simultaneous induction of the stress-inducible system (Figure 3D). The delay that we observed in the production of the GFP is probably due to the basal activity of the IbpAB promoter.
Influence of the stress feedback system on protein solubility
Our analysis shows that the native expression system without feedback leads to the highest amount of total GFP, but at the same time with a high percentage of insoluble GFP in the cell (Figure 4B and 4C). Some variants with the feedback-based expression system were found to increase the soluble fraction, albeit at the expense of lower protein yield. The observed dynamics are consistent with what we would expect due to the trade-off between protein solubility and yield. The reduced protein yield is expected as the feedback mechanism is designed to reduce GFP expression upon stress. Interestingly, variability across clones was also observed here, with some clones retaining the protein yield and quality of the native expression system, while others favoring protein solubility over protein yield, for the wild type P IbpAB promoter (with strong RBS) (Additional file 1: Figure S1). The intracellular levels of TetR protein were analyzed by Western blot analysis and support data obtained throughout the study showing that different variants and combinations of stress promoter and RBS indeed resulted in different intracellular TetR levels, with similar variability as in the case of GFP expression (Figure 4D).
Effect of growth parameters on protein yield
The production of recombinant product was modeled by dividing the process into three main sub-components: (a) the induction of T7 polymerase production by IPTG, (b) the recombinant protein production via the T7 promoter and (c) the stress-induced expression of the TetR repressor. A system of delayed differential equations was used to capture the transcription and translation processes for T7 polymerase, GFP and TetR. In addition, the ratio between GFP fractions was also modeled and visualized. A detailed mathematical description of the model, together with all experimentally derived and estimated parameters, is given in the supplementary materials section. We performed sensitivity analysis of our model (Additional file 1: Figure S4) that identified two key processes that are crucial for its robustness. In addition, we performed cross-validation to avoid over-fitting and evaluate the generalization error of the model (Additional file 1: Table S9).
In our quest to create self-regulatory systems for recombinant protein production, we used an integrated synthetic biology approach to construct a synthetic circuit that limits recombinant protein production through stress-induced feedback. We validated the functionality of different variants of the synthetic circuit in their capacity to limit stressful protein production, and to increase the total soluble fraction. Since the protein yield was significantly lowered in the process, further investigation on promoter and repressor engineering for avoiding such loss would be welcomed. In addition, since different proteins lead to different levels of stress within the host cells , it would be interesting to test this approach with other recombinant protein species. Furthermore different inducer concentrations can be used to tune transcription rate and product formation, although there is evidence that inducer concentration does not necessarily influence the formation of active soluble protein .
The computational model that we constructed provided valid predictions on the system dynamics, and was useful as a first order guiding tool for our experimental design. An extended phenomenological description and inclusion of a larger set of measured experimental parameters would allow an increased predictive accuracy of the model, and it may help to test or generate alternative hypotheses regarding the the dynamics of inclusion body formation and their degradation rates. This study provides an example of how integration of computational, engineering and experimental methods, together with the synthetic biology concepts of parts standardization, can be applied to address biotechnological challenges from a new perspective. Altogether, this and similar future studies can be applied to guide the construction of robust auto-regulatory protein production systems.
Host strains and growth media
E. coli DH5α was used for all cloning procedures, whereas E. coli C41  was used as host expression strain. All cultivations were performed on LB medium and M9 minimal medium (0.4% w/v glucose) supplemented with antibiotics (carbenicillin 100 μg mL-1 and chloramphenicol 25 μg mL-1) when necessary and incubated at 37°C and 150 rpm on an orbital shaker. For experiments involving GFP expression and synthetic circuit characterizations, cultures were grown on M9 medium and induced with 1 mM IPTG (1 M 0.22 μm filtered sterile stock solution). Cultures were inoculated at an optical density (600 nm) of 0.1 and grown for 2 h before induction with IPTG. For experiments involving arabinose as inducer substance, arabinose was added to the growth medium in a concentration range of 0-1% (from a 20% w/v 0.22 μm filtered sterile stock solution). Samples were taken in regular intervals to monitor growth and product formation. All tests involved at least 3 biological replicates and were performed in 5 mL growth medium and 50 mL tubes unless stated otherwise.
Molecular cloning procedures
For standard molecular cloning restriction enzymes were purchased from New England Biolabs Inc. Taq polyermase (Qiagen) was used for screening and error prone PCR whereas other PCR reactions were performed using PfuTurbo proofreading polymerase (Stratagene). The native DNA sequence of P IbpAB was amplified from a genomic DNA isolation of E. coli MG1655. Subsequent error prone PCR was essentially performed as described previously . To achieve higher error rates during PCR, imbalanced dNTPs (0.2 mMdATP and GTP, 1 mMdCTP and dTTP) and increased MgCl2 (20 mM) and MnCl2 (0-0.5 mM) concentrations were applied. Modification of ribosome binding sites was performed by mismatch oligonucleotide PCR. Plasmid mini prep DNA was used as template for PCR with primers harboring the desired nucleotide mutations of the RBS. After PCR parental plasmid was digested with the methylation sensitive restriction enzyme DpnI and linear mutated plasmid DNA was transformed into E. coli DH5α. Plasmid DNA from single colonies was isolated and sequenced in order to select mutated sequences.
Primer sequences and further information on cloning procedures as well as standard parts and newly constructed DNA parts used in this study can be found in Supplementary Online Materials (Text and Additional file 1: Table S1, Table S2 and Table S3, respectively).
OD600 measurements and fluorescence measurements
Optical density, as a measure of cell mass, was determined on a Biophotometer (Biorad) and/or an Infinite M200 plate reader (Tecan). Fluorescence measurements of GFP were performed on an Infinite M200 microtiter-plate reader (Tecan) using the following settings: 37°C, excitation wavelength 485 nm, emission wavelength 535 nm and gain 75. GFP fluorescence values in the manuscript are given as GFP fluorescence per OD600.
Analysis of soluble and insoluble proteins and western blotting
Preparations of soluble and insoluble protein were performed according to the pET vector manual (Novagen, pET system manual, 11th edition, 2006) with small modifications. In short, cells were harvested by centrifugation, wet cell weight was determined and samples were immediately frozen at -20°C for at least 24 h. For protein analysis, samples were thawed to room temperature and re-suspended in cell lysis buffer. Lysozyme (Sigma) was added (60 kU per g wet biomass) and samples were incubated at room for 45 min on a rotary shaker. To complete cell lysis and decrease the viscosity of the solution, samples were treated in a Bioruptor sonication apparatus (Diagenode) at high power setting and a 30s sonication interval for 10 min. After centrifugation the supernatant containing soluble proteins was separated from insoluble cell fraction. The pellet was re-suspended in cell lysis buffer containing 1% SDS and kept as insoluble protein fraction. SDS PAGE analysis was performed using 12% BisTris polyacrylamide Novex gels (Invitrogen) according to standard protocols. Gels were stained with PageBlue staining solution (Thermo Scientific). One sample was loaded on all of the gels to account and correct for variations in gel staining. After scanning of gels, protein bands were quantified using GelAnalyzer v2010a http://www.gelanalyzer.com/. Intensity of gel bands was corrected for staining differences between individual gels and normalized to the wet cell weight of each sample to account for the different initial sample volumes.
Intracellular levels of TetR protein were determined by Western Blotting. Total soluble protein was separated by standard SDS-PAGE as described and blotted onto a PVDF membrane. After membrane blocking (TBS 0.1% Tween-20, 1% milk protein) at 4°C over night, TetR protein was detected using an anti-TetR polyclonal rabbit antibody serum (Thermo Scientific), a polyclonal anti-rabbit Peroxidase conjugate (Sigma Aldrich) and a BCL plus detection kit (GE Healthcare) using standard WB protocols. Western Blot signals were acquired on a Storm scanner (GE Healthcare) at 520 nm and analyzed using ImageQuant image analysis software (GE Healthcare). Normalization of TetR intensities on the blots was performed as described for GFP levels above.
The mathematical four replicates, where bars equations that incorporate the specific time required for particular events, such as transcription and translation. The equations were simulated through the "dde23" routine in MATLAB 7.10 (The MathWorks, Natick, MA). All results are displayed for each species over 12 hours of batch culture with induction at two hours. The system was allowed reach steady state prior to t = 0 hours for an extended period of time. During this pre-induction phase, the growth rate was set to μ 0 . Simulation results were obtained under induction by 1 mM of IPTG. In order to compare the model's single cell results, experimental data was converted from normalized fluorescence to molarity . Detailed transcription of the model can be found in Supplementary Material online.
This work was supported by NSF-CCF 1146926 and the UC Davis Opportunity Fund. We would like to thank Marc Facciotti, Karen McDonald, and the members of the Tagkopoulos Lab for the helpful discussions.
- Itakura K, Tadaaki H, Crea R, Riggs AD, Heyneker HL, Bolivar F, Boyer HW: Expression in Escherichia coli of a chemically synthesized gene for the hormone somatostatin. 1977. Biotechnology. 1992, 24: 84-91.Google Scholar
- Graumann K, Premstaller A: Manufacturing of recombinant therapeutic proteins in microbial systems. Biotechnol J. 2006, 1: 164-186. 10.1002/biot.200500051.View ArticleGoogle Scholar
- Gasser B, Mattanovich D: Antibody production with yeasts and filamentous fungi: on the road to large scale?. BiotechnolLett. 2007, 29: 201-212.Google Scholar
- Vaishnav P: Production of recombinant proteins by microbes and higher organisms. BiotechnolAdv. 2009, 27: 297-306.Google Scholar
- Hamilton S, Gerngross T: Glycosylation engineering in yeast: the advent of fully humanized yeast. CurrOpinBiotechnol. 2007, 18: 387-392.Google Scholar
- Baneyx F, Mujacic M: Recombinant protein folding and misfolding in Escherichia coli. Nat Biotechnol. 2004, 22: 1399-1408. 10.1038/nbt1029.View ArticleGoogle Scholar
- Wacker M, Linton D, Hitchen P, Nita-Lazar M, Haslam S, North S, Panico M, Morris H, Dell A, Wren B, Aebi M: N-linked glycosylation in Campylobacter jejuni and its functional transfer into E. coli. Science. 2002, 298: 1790-1793. 10.1126/science.298.5599.1790.View ArticleGoogle Scholar
- Ferrer-Miralles N, Domingo-Espín J, Corchero J, Vázquez E, Villaverde A: Microbial factories for recombinant pharmaceuticals. Microb Cell Fact. 2009, 8: 17-10.1186/1475-2859-8-17.View ArticleGoogle Scholar
- Gill R, Valdes J, Bentley W: A comparative study of global stress gene regulation in response to overexpression of recombinant proteins in Escherichia coli. MetabEng. 2000, 2: 178-189.Google Scholar
- Allen S, Polazzi J, Gierse J, Easton A: Two novel heat shock genes encoding proteins produced in response to heterologous protein expression in Escherichia coli. J Bacteriol. 1992, 174: 6938-6947.Google Scholar
- Martínez-Alonso M, González-Montalbán N, García-Fruitós E, Villaverde A: Learning about protein solubility from bacterial inclusion bodies. Microb Cell Fact. 2009, 8: 4-10.1186/1475-2859-8-4.View ArticleGoogle Scholar
- Nahálka J, Vikartovská A, Hrabárová E: A crosslinked inclusion body process for sialic acid synthesis. J Biotechnol. 2008, 134: 146-153. 10.1016/j.jbiotec.2008.01.014.View ArticleGoogle Scholar
- Schumann W: Production of recombinant proteins in Bacillus subtilis. AdvApplMicrobiol. 2007, 62: 137-189.Google Scholar
- Jana S, Deb J: Strategies for efficient production of heterologous proteins in Escherichia coli. ApplMicrobiolBiotechnol. 2005, 67: 289-298.Google Scholar
- Gasser B, Saloheimo M, Rinas U, Dragosits M, Rodríguez-Carmona E, Baumann K, Giuliani M, Parrilli E, Branduardi P, Lang C, et al: Protein folding and conformational stress in microbial cells producing recombinant proteins: a host comparative overview. Microb Cell Fact. 2008, 7: 11-10.1186/1475-2859-7-11.View ArticleGoogle Scholar
- Ramos CR, Abreu PA, Nascimento AL, Ho PL: A high-copy T7 Escherichia coli expression vector for the production of recombinant proteins with a minimal N-terminal His-tagged fusion peptide. Braz J Med Biol Res. 2004, 37: 1103-1109.View ArticleGoogle Scholar
- Chuang S, Burland V, Plunkett Gr, Daniels D, Blattner F: Sequence analysis of four new heat-shock genes constituting the hslTS/ibpAB and hslVU operons in Escherichia coli. Gene. 1993, 134: 1-6. 10.1016/0378-1119(93)90167-2.View ArticleGoogle Scholar
- Dürrschmid K, Reischer H, Schmidt-Heck W, Hrebicek T, Guthke R, Rizzi A, Bayer K: Monitoring of transcriptome and proteome profiles to investigate the cellular response of E. coli towards recombinant protein expression under defined chemostat conditions. J Biotechnol. 2008, 135: 34-44. 10.1016/j.jbiotec.2008.02.013.View ArticleGoogle Scholar
- Lesley S, Graziano J, Cho C, Knuth M, Klock H: Gene expression response to misfolded protein as a screen for soluble recombinant protein. Protein Eng. 2002, 15: 153-160. 10.1093/protein/15.2.153.View ArticleGoogle Scholar
- Choi J, Lee S, Lee S: Enhanced production of insulin-like growth factor I fusion protein in Escherichia coli by coexpression of the down-regulated genes identified by transcriptome profiling. Appl Environ Microbiol. 2003, 69: 4737-4742. 10.1128/AEM.69.8.4737-4742.2003.View ArticleGoogle Scholar
- Jürgen B, Lin H, Riemschneider S, Scharf C, Neubauer P, Schmid R, Hecker M, Schweder T: Monitoring of genes that respond to overproduction of an insoluble recombinant protein in Escherichia coli glucose-limited fed-batch fermentations. BiotechnolBioeng. 2000, 70: 217-224.Google Scholar
- Orth P, Schnappinger D, Hillen W, Saenger W, Hinrichs W: Structural basis of gene regulation by the tetracycline inducible Tet repressor-operator system. Nat StructBiol. 2000, 7: 215-219. 10.1038/73324.View ArticleGoogle Scholar
- Kim Y, Burton R, Burton B, Sauer R, Baker T: Dynamics of substrate denaturation and translocation by the ClpXP degradation machine. Mol Cell. 2000, 5: 639-648. 10.1016/S1097-2765(00)80243-9.View ArticleGoogle Scholar
- Stricker J, Cookson S, Bennett M, Mather W, Tsimring L, Hasty J: A fast, robust and tunable synthetic gene oscillator. Nature. 2008, 456: 516-519. 10.1038/nature07389.View ArticleGoogle Scholar
- Elowitz MB, Leibler S: A synthetic oscillatory network of transcriptional regulators. Nature. 2000, 403: 335-338. 10.1038/35002125.View ArticleGoogle Scholar
- Cormack BP, Valdivia RH, Falkow S: FACS-optimized mutants of the green fluorescent protein (GFP). Gene. 1996, 173: 33-38. 10.1016/0378-1119(95)00685-0.View ArticleGoogle Scholar
- Lee S, Chou H, Pfleger B, Newman J, Yoshikuni Y, Keasling J: Directed evolution of AraC for improved compatibility of arabinose- and lactose-inducible promoters. Appl Environ Microbiol. 2007, 73: 5711-5715. 10.1128/AEM.00791-07.View ArticleGoogle Scholar
- Hartinger D, Schwartz H, Hametner C, Schatzmayr G, Haltrich D, Moll WD: Enzyme characteristics of aminotransferase FumI of Sphingopyxis sp. MTA144 for deamination of hydrolyzed fumonisin B(1). ApplMicrobiolBiotechnol. 2011, 91 (3): 757-68.Google Scholar
- Miroux B, Walker JE: Over-production of proteins in Escherichia coli: mutant hosts that allow synthesis of some membrane proteins and globular proteins at high levels. J MolBiol. 1996, 260: 289-298.View ArticleGoogle Scholar
- Martin CT, Coleman JE: Kinetic analysis of T7 RNA polymerase-promoter interactions with small synthetic promoters. Biochemistry. 1987, 26: 2690-2696. 10.1021/bi00384a006.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.