- Research article
- Open Access
Dynamics of receptor and protein transducer homodimerisation
© Vera et al; licensee BioMed Central Ltd. 2008
- Received: 19 December 2007
- Accepted: 31 October 2008
- Published: 31 October 2008
Signalling pathways are complex systems in which not only simple monomeric molecules interact, but also more complex structures that include constitutive or induced protein assemblies. In particular, the hetero-and homo-dimerisation of proteins is a commonly encountered motif in signalling pathways. Several authors have suggested in recent times that dimerisation relates to a series of physical and biological outcomes used by the cell in the regulation of signal transduction.
In this paper we investigate the role of homodimerisation in receptor-protein transducer interactions. Towards this end, mathematical modelling is used to analyse the features of such kind of interactions and to predict the behaviour of the system under different experimental conditions. A kinetic model in which the interaction between homodimers provokes a dual mechanism of activation (single and double protein transducer activation at the same time) is proposed. In addition, we analyse under which conditions the use of a power-law representation for the system is useful. Furthermore, we investigate the dynamical consequences of this dual mechanism and compare the performance of the system in different simulated experimental conditions.
The analysis of our mathematical model suggests that in receptor-protein interacting systems with dual mechanism there may be a shift between double and single activation in a way that intense double protein transducer activation could initiate and dominate the signal in the short term (getting a fast intense signal), while single protein activation could control the system in the medium and long term (when input signal is weaker and decreases slowly). Our investigation suggests that homodimerisation and oligomerisation are mechanisms used to enhance and regulate the dynamic properties of the initial steps in signalling pathways.
- Kinetic Order
- Protein Transducer
- Dual Mechanism
- Inactive Protein
- Double Activation
The processing of information in living cells is carried out by signal transduction pathways . Through the binding of external ligands to extracellular receptors, the cell can receive signals from its environment and transfer information into the cell. This information flow is regulated, amplified or modulated by different feedback mechanisms and interactions with other pathways (crosstalk). Moreover, signalling pathways are complex systems in which not only simple monomeric molecules interact but also more complex structures that include constitutive or induced protein assemblies [2–4]. In particular, the hetero- and homo-dimerisation of proteins is a commonly encountered motif in signalling pathways.
In Klemm  the role of dimerisation as a regulatory mechanism in signal transduction is analysed and discussed. Dimerisation is defined as an interaction producing a protein-protein complex composed of two subunits, either identical (homodimerisation) or non-identical (heterodimerisation). The authors argue that dimerisation relates to a series of physical and biological outcomes used by the cell in the regulation of signal transduction. The biophysical outcomes referred to facilitation of proximity and orientation in protein interaction, differential regulation through heterodimerisation, emergence of spatio-temporal boundaries, enhanced specificity and regulation of monomer-to-dimer transitions. The role of homodimeric receptors in the activation and dimerisation of intermediate proteins and transcription factors has been already described in the literature [6–8].
A well-studied example are the JAK/STAT signalling pathways . In case of the JAK2/STAT5 signalling pathway, the Epo receptor is a preformed inactive dimer in the plasma membrane [10, 11]. The binding of Epo results in the activation of the JAK2 kinase and subsequent phosphorylation of the cytosolic domain of each Epo receptor monomeric subunit. STAT5 proteins bind to the tyrosine phosphorylated Epo receptor and gets phosphorylation. Afterwards, they dimerised and translocate to the nucleus. The spatial conformation of the receptor as a dimer seems to indicate that each activated Epo receptor monomer could phosphorylate simultaneously at least one STAT5 molecule. The correspondence between the existence of a homodimer activated receptor and the activation of a homodimer transduction protein suggests a possible complex underlying molecular mechanism for the activation and dimerisation process. Similar behaviour has been shown in JAK/STAT pathways  and other signalling pathways [13–18], suggesting that this homodimer-homodimer interaction could constitute a more general pattern in cell signalling systems.
The purpose of this work is to use mathematical modelling to suggest mechanisms of interaction by which this homodimer-homodimer interaction can occur. In our work, two mathematical modelling frameworks are used and compared. We furthermore investigate the dynamical consequences of the interaction mechanisms suggested and propose general features of an experimental design to discriminate between the different mechanisms.
In this paper we support the hypothesis that in the transduction of signals via homodimeric proteins, the dimeric nature of the receptor plays an essential role in the fast response of the biological system. The surface density of many plasma membrane receptors tends to be very low . For a low density of receptors at the plasma membrane, active mechanisms for the homodimerisation of the cytosolic interacting proteins are required to boost the intracellular response of the system to external stimuli. Towards this end, we propose that the dimeric structure of the receptor allows not only the simple activation of a monomer binding protein per receptor but also a simultaneous coordinated activation of two monomers of the binding protein, one at each subunit of the receptor. The case in which we consider the single interaction between one subunit of the receptor and a monomer of the protein transducer, we call "single protein activation process" with respect to the number of receptor-protein interactions. In contrast, two subunits of the same receptor interacting simultaneously with two different monomers are referred to as a "double protein activation process". In the following we are going to assume that higher-order processes, allowing for simultaneous activation of two or more protein monomers at the same receptor monomer, are not significant. Several studies with the PDGF receptor in which all tyrosine phosphorylation sites were mutated individually or in combination showed that the signalling was unchanged when redundant sites were knocked out, suggesting that this assumption is in principle meaningful .
where γ j are rate coefficients, g jk kinetic orders and σ ij the stoichiometric coefficients. The stoichiometric coefficient σ ij is positive for products and negative for reactants and describes how many molecules of X i are converted in the considered reaction. The rate coefficients contain information about the physical properties of the reaction, like activation energies and internal states . If we assume that environmental conditions, like temperature, pressure and pH-value do not vary over time, γ j is a rate constant. The interpretation of the kinetic order g jk depends on the chosen level of description for the biochemical network. When (1) is derived within the framework of statistical physics and the complete reaction mechanism is known and considered, the coefficient g jk has a clear mechanistic meaning and defines the number of molecules of species X i involved on the considered reaction. The kinetic order has in this case a positive integer value. In contrast, if we aggregate elementary reaction steps into a combined power-law expression, with the same form of (1), the interpretation of the coefficient g jk changes and they can take non-integer values [22–25]. This latter case is investigated in the next section.
For the purpose of our investigation, we consider the generic case of a diffusive signalling protein (for simplicity, named protein, P, in our discussion), which is activated and dimerised after the interaction with a homodimeric activated receptor (named receptor, R*). In order to focus the investigation on general design principles underlying the homodimer receptor-homodimer transducer protein class of interactions, we simplify the mathematical modelling of the system, deliberately neglecting details that account for differences in the class. The procedure is similar to the one used in .
Equations (2a) and (2b) describe the independent monomer activation, whereas (2c) describes the formation of the dimer in the cytosol. In addition to its concentration, each reaction is determined by a rate constant k i , encapsulating physical information about the underlying biochemical reaction. For simplification purposes we apply additional assumptions. Firstly, since in our simplified model there is no production or additional recruitment of protein transducer P, there is an intrinsic conservation for P and P* in the model. In addition, we assume that the receptor activity remains constant as consequence of a constant external signal . The dynamic control of receptor activity/concentration and the recruitment/recycling of protein transducers  are important features in signalling systems that must be considered when modelling specific signalling pathways, but can be neglected here for the purpose of generalisation. Finally, we furthermore assume that in the class of interactions investigated, the intermediate state R*P cannot go down a monomeric activation route .
where the first term describes the double protein activation process and the second the dimer formation. In contrast to (2a) and (2b), the molecularity of protein P in the process of phosphorylation is now two instead of one. The process of dimerisation has the same structure as in the single protein activation process.
Because we focus on the consequences of the combination of single protein activation and double protein activation mechanism, we neglect some internal details of the double mechanisms and merge it into a single expression (For more details [see Additional file 1]). The temporal changes of the activated protein concentration are given now as
then the double protein activation is dominant. The apparent kinetic order tends to the limiting value of g2 = 2. Between the limiting cases discussed above, both activation mechanisms contribute to the dynamics and the system is not simply single or double protein activation. An elementary, positive and integer, kinetic order cannot be assigned to the overall reaction and an apparent positive non-integer kinetic order occurs. The origin of this apparent kinetic order as well as its dependence on the kinetic constants are discussed in the section below.
An experiment measuring the concentration of the protein, active and inactive receptors and/or protein activated dimers, cannot distinguish between active monomers which are produced in a single or double protein activation process. To distinguish the different activation mechanisms, an indirect method is required, another possibility is to investigate the structure of the protein. However, the detection of two binding sites at the same receptor, does not guarantee that both sites are used simultaneously nor that this is necessarily an effective activation process.
In order to investigate the possibility of two different and simultaneously acting phosphorylation mechanisms further, we aggregate both processes introduced in the above section into a single contribution. The order of the receptor R in the activation process is one, as in the mechanistic model (8). However, the contribution of protein P is a combination of the considered single and double activation processes. Since both activation mechanisms described by k1 and k3 are indistinguishable, the estimation of these kinetic parameters generates identifiability issues. If we consider a situation in which both processes cannot be distinguished, a feasible way to reproduce this complexity is the use of a power-law representation, allowing for non-integer kinetic orders. As discussed in , power-law models are useful for modelling cellular signalling when the exact reaction mechanism is unknown or if experimental data are not sufficient. For an example in cell signalling we refer to : in that work, quantitative time course data were used to identify a power-law model. If we used a power-law term, the production term (phosphorylation of P) of the new rate equation has the following structure:
V(R*, P) ≈ γ1·R·P g .
The apparent kinetic order g, which describes the role of the protein, has now a non-integer value between one and two. The specific value of g depends on the rate constants k1, k3 and the concentration of the inactive protein P as shown in the subsection "Mechanistic modelling". We note that this apparent kinetic order is only a mathematical analogue to the kinetic order of elementary reactions. It does not provide detailed information about the underlying mechanisms and hence, it cannot be interpreted in the same way as the kinetic order of an elementary reaction scheme. In addition to the apparent kinetic order g, we also introduce an apparent kinetic constant γ1 in the power-law representation (9). Due to the aggregation of production terms in (9), this apparent kinetic constant does not coincide with the kinetic constant in the mechanistic rate law (8). Nevertheless, it depends on the kinetic constants k1 and k3 but also on the ratio of efficiency of single and double protein activation processes.
In the general case and under the assumption that we do not know the exact relation between the underlying (indistinguishable) processes (encoded by the ratio k3/k1) the power-law representation allows us to reproduce essential dynamical properties of the system, regardless whether the ratio is low, medium or high [see Additional File 1]. A power-law model can therefore be used to elucidate the nature of the underlying mechanism. An estimated value for the kinetic order g near one would suggest that the single protein activation mechanism is the principal or unique contributor to the activation process while a value near two would indicate the relevance of the double protein activation process. An intermediate value would indicate that the mechanism has a dual nature as described in the previous section.
Dynamical Consequences of the Dual Mechanism
The comparison of the three situations shows clear differences between the different scenarios. If the system allows only the single protein activation process for P, the response of the system in both mutants would be indistinguishable (Figure 5[a]). In contrast, if the system allows the dual mechanism, single and double protein activation, the response is different for Mut1 and Mut2 (Figure 5[b]). Moreover, on the assumption of an intermediate value for the ratio k3 P/k1 in this case, the differences in the response between wild-type and mutant cells are more significant (Figure 5[b]). Finally, if the activation mechanism is a double protein activation process, only Mut1 should produce a significant signal after stimulation (Figure 5[c]).
A system, where the stimulation with the same input signal produces identical response in both mutants, would not present the dual mechanism of activation but only the single activation. A system where the responses for the same input signal are differentiated in both mutants would present the dual mechanism of activation. Finally, a system in which the stimulation of the system produces a signal for the Mut1 but not for the Mut2 (where any double protein activation process is intentionally blocked) is a system where only the double protein activation process is possible.
In case of k3 > k1, the concentration of inactivate protein for which both signal rates contribute with an identical amount of signal (k1Peq = k3(Peq)2) is defined by the ratio between rate the constants k1/k3 (Figure 6[b]). The higher the value of this ratio is, the more reduced the effect of double activation in the dynamics of the system. On the other hand, reduced values of the ratio imply that double activation dominates even at much reduced concentrations of inactive protein. When the value of k3 increases the position of Peq is shifted to smaller values of P.
In order to establish the domain of values for the rate constants in which the mechanism is effectively dual, we fix an interval of feasible values for the parameter k3 with respect to k1. The minimum value for k3 is such that double protein activation rate contributes a 10% of the single activation signal for the maximum amount of inactive protein: rDA(P = 1) = 0.1·rSA(P = 1.0). Below this value for k3, the contribution of double activation (rDA) is negligible for any concentration of P. The maximum value of k3 is such that the double protein activation is dominant for values of P higher than 10% of the total amount: Peq = k1/k3 = 0.1. Lower values for Peq implies that single protein activation process does not contribute at any significant concentration of P. Thus, we define the interval of values for k3 in which the mechanism is dual like k1/10 ≤ k3 ≤ 10k1. In this interval of values, the system presents a dual mechanism, with dominance of the single protein activation for low concentration of inactive protein and dominance of double protein activation for high values of P.
One of the dynamical consequences of this dual mechanism relates to the behaviour of the system during transient stimulation. Under certain circumstances there could be a switch between both mechanisms, with a dominance of the double protein activation in the short initial time (when the amount of inactive protein is high) and a take over the dynamics by the single protein activation during the medium-long term (when the amount of inactive protein is smaller).
In this paper, the possibility of a dual mechanism of activation for homodimer proteins by homodimer receptors is analysed. We propose a mechanistic model in which this peculiar interaction between homodimers provokes a dual mechanism of activation, which is simultaneously single and double protein activation process. Under the stated assumptions, the dual mechanism could appear essentially as either single or double activation, depending on the ratio between the rate constants of both process and the available amount of inactive protein. For a certain interval of values for these coefficients, the process can have an intermediate behaviour, which cannot be reduced to the single or double activation mechanism.
We further analysed a phenomenological power-law representation, which is able to adequately simulate the response of the system. We demonstrated that power-law representation is useful for cases where experimental data are available, but which are not sufficient to distinguish and characterise the reaction mechanisms (single or double activation). The power-law model, together with quantitative experimental data, is an alternative to investigate the structure of a given pathway where the estimated value for the kinetic order is used to decide whether the process is mainly of single activation (g ≈ 1), double activation (g ≈ 2) or dual (intermediate values).
Finally, the consequences of this dual mechanism on the dynamics were investigated through the simulation of different experimental conditions. Towards this end, we simulated the behaviour of the system for two mutants: a mutant with a reduced number of receptors, and a second mutant in which one of the subunits of the receptors is blocked for the activation of the studied protein. Only in cases in which the system has the dual mechanism of activation both mutants would induce a different response, while a pure double protein activation system would not produce a significant response to the input signal, when one of the subunits of the homodimer receptors is constitutively blocked for the activation of the studied protein. Furthermore, our analysis suggests that in receptor-protein interacting systems with dual mechanism there may be an active switch between double and single activation in a way that intense double protein activation could initiate and dominate the signal in the short term (getting a fast intense signal), while single protein activation could control the system in the medium and long term (when input signal is weaker and decreases slowly).
Open questions that require attention relate to the underlying reasons that justify the dual activation mechanism. Mathematical modelling could be used to investigate why and how this dual mechanism is used by the cell to improve the performance of monomer-to-dimer transitions . Moreover, there is increasing evidence for feedback inhibitor proteins that bind to activated receptors and compete with the true effectors for binding, e.g. CIS1, and Mig6 [32–34]. The expression of these feedback inhibitors is usually induced by the activated receptor. Thus, there may be an implicit switch between first and second order processes over the timecourse of stimulation dependent on feedback inhibitor levels that could be investigated in the next future.
The authors thank the collaboration of Taesoo Kwon and Ulf Schmitz elaborating a list of published homod-imer receptor-homodimer transducer protein interactions used in initial versions of this work. This work was supported by the European Commission 6th Framework program and as part of the COSBICS project under contract LSHG-CT-2004-512060 http://www.sbi.uni-rostock.de/cosbics and by the German Federal Ministry of Education and Research (BMBF) through the projects BaCell-SysMo (0313978F) and COSMIC-SysMo (0313981D).
- Downward J: The ins and outs of signalling. Nature. 2001, 411 (6839): 759-762. 10.1038/35081138View ArticlePubMedGoogle Scholar
- Jones N: Transcriptional regulation by dimerization: two sides to an incestuous relationship. Cell. 1990, 61: 9-11. 10.1016/0092-8674(90)90207-UView ArticlePubMedGoogle Scholar
- Jones S, Thornton J: Protein-protein interactions: a review of protein dimer structures. Prog Biophys Mol Biol. 1995, 63: 31-65. 10.1016/0079-6107(94)00008-WView ArticlePubMedGoogle Scholar
- Pereira-Leal J, Levy E, Kamp C, Teichmann S: Evolution of protein complexes by duplication of homomeric interactions. Gen Biol. 2007, 8 (4): R51-10.1186/gb-2007-8-4-r51.View ArticleGoogle Scholar
- Klemm J, Schreiber S, Crabtree G: Dimerization as a Regulatory Mechanism in Signal Transduction. Annu Rev Immunol. 1998, 16: 569-592. 10.1146/annurev.immunol.16.1.569View ArticlePubMedGoogle Scholar
- Reich Z, Boniface J, Lyons D, Borochov N, Wachtel E, Davis M: Ligand-specific oligomerization of T-cell receptor molecules. Nature. 1997, 387: 617-620. 10.1038/42500View ArticlePubMedGoogle Scholar
- Schlessinger J: Signal transduction by allosteric receptor oligomerization. Trends Biol Sci. 1988, 13: 443-447. 10.1016/0968-0004(88)90219-8.View ArticleGoogle Scholar
- Heldin C: Dimerization of cell surface receptors in signal transduction. Cell. 1995, 80: 213-223. 10.1016/0092-8674(95)90404-2View ArticlePubMedGoogle Scholar
- Ferrajoli A, Faderl S, Ravandi F, Estrov Z: The JAK-STAT pathway: a therapeutic target in hemato-logical malignancies. Curr Cancer Drug Targets. 2006, 6 (8): 671-679. 10.2174/156800906779010227View ArticlePubMedGoogle Scholar
- Livnah O, Stura E, Middleton S, Jollife L, Wilson I: Crystallographic evidence for preformed dimers of erythropoietin receptor before ligand activation. Science. 1999, 283 (5404): 987-990. 10.1126/science.283.5404.987View ArticlePubMedGoogle Scholar
- Remy I, Wilson I, Michnick S: Erythropoietin receptor activation by a ligand-induced conformation change. Science. 1999, 283: 990-993. 10.1126/science.283.5404.990View ArticlePubMedGoogle Scholar
- Shuai K, Horvath C, Huang L, Qureshi S, Cowburn D, Darnell J: Interferon activation of the transcription factor Stat91 involves dimerization through SH2-phosphotyrosyl peptide interactions. Cell. 1994, 76: 821-828. 10.1016/0092-8674(94)90357-3View ArticlePubMedGoogle Scholar
- Sakurai H, Nigam S: In vitro branching tubulogenesis: Implications for developmental and cystic disorders, nephron number, renal repair, and nephron engineering. Kidney International. 1998, 54: 14-26. 10.1046/j.1523-1755.1998.00969.xView ArticlePubMedGoogle Scholar
- Jain N, Zhang T, Kee W, Li W, Cao X: Protein kinase C delta associates with and phosphorylates Stat3 in an interleukin-6-dependent manner. J Biol Chem. 1999, 274 (34): 24392-400. 10.1074/jbc.274.34.24392View ArticlePubMedGoogle Scholar
- Moustakas A, Souchelnytskyi S, Heldin C: Smad regulation in TGF-β signal transduction. Journal of Cell Science. 2001, 114: 4359-4369.PubMedGoogle Scholar
- Stamper C, Zhang Y, Tobin J, Erbe D, Ikemizu S, Davis S, Stahl M, Seehra J, Somers W, Mosyak L: Crystal structure of the B7-1/CTLA-4 complex that inhibits human immune responses. Nature. 2001, 410: 608-611. 10.1038/35069118View ArticlePubMedGoogle Scholar
- Shao H, Cheng H, Cook R, Tweardy D: Identification and characterization of signal transducer and activator of transcription 3 recruitment sites within the epidermal growth factor receptor. Cancer Res. 2003, 63 (14): 3923-30.PubMedGoogle Scholar
- Schulze W, Deng L, Mann M: Phosphotyrosine interactome of the ErbB-receptor kinase family. Mol Sys Biol. 2005, 1: 2005.0008-Google Scholar
- Winter S, Howard T, Ware R: Regulation of Expression of the Human Erythropoietin Receptor Gene. Blood Cell Mol Dis. 1996, 22 (3): 214-224. 10.1006/bcmd.1996.0102.View ArticleGoogle Scholar
- Fambrough D, McClure K, Kazlauskas A, Lander E: Diverse signaling pathways activated by growth factor receptors induce broadly overlapping, rather than independent, sets of genes. Cell. 1999, 97: 727-741. 10.1016/S0092-8674(00)80785-0View ArticlePubMedGoogle Scholar
- Hänggi P, Talkner P, Borkovec M: Reaction-rate theory: fifty years after Kramers. Rev Mod Phys. 1990, 62 (2): 251-341. 10.1103/RevModPhys.62.251.View ArticleGoogle Scholar
- Savageau M: Biochemical systems analysis: I. Some mathematical properties of the rate law for the component enzymatic reactions. J Theor Biol. 1969, 25: 365-369. 10.1016/S0022-5193(69)80026-3View ArticlePubMedGoogle Scholar
- Savageau M: Biochemical systems analysis: II. Steady state solutions for an n-poll system using a power-law approximation. J Theor Biol. 1969, 25: 370-379. 10.1016/S0022-5193(69)80027-5View ArticlePubMedGoogle Scholar
- Savageau M: Biochemical systems analysis: III. Dynamic solutions using a power-law approximation. J Theor Biol. 1970, 26: 215-226. 10.1016/S0022-5193(70)80013-3View ArticlePubMedGoogle Scholar
- Vera J, Balsa-Canto E, Wellstead P, Banga J, Wolkenhauer O: Power-Law Models of Signal Transduction Pathways. Cell Signal. 2007, 19 (7): 1531-41. 10.1016/j.cellsig.2007.01.029View ArticlePubMedGoogle Scholar
- Barkai N, Leibler S: Robustness in simple biochemical networks. Nature. 1997, 387 (6636): 913-917. 10.1038/43199View ArticlePubMedGoogle Scholar
- Lauffenburger D, Linderman J: Receptors: Models for binding, trafficing, and signalling. 1993, Oxford University PressGoogle Scholar
- Swameye I, Müller T, Timmer J, Sandra O, Klingmüller U: Identification of nucleocytoplasmic cycling as a remote sensor in cellular signaling by data-based modeling. Proc Natl Acad Sci USA. 2003, 100 (3): 1028-1033. 10.1073/pnas.0237333100PubMed CentralView ArticlePubMedGoogle Scholar
- Vera J, Bachmann PAJ, Becker V, Hormiga J, Torres Darias N, Timmer J, Klingmuller U, Wolkenhauer O: A systems biology approach to analyse amplification in the JAK2-STAT5 signalling pathway. BMC Systems Biology. 2008, 2: 38- 10.1186/1752-0509-2-38PubMed CentralView ArticlePubMedGoogle Scholar
- Hernandez-Bermejo B, Fairen V, Sorribas A: Power-law modeling based on least-squares criteria: consequences for system analysis and simulation. Math Biosci. 2000, 167: 87-107. 10.1016/S0025-5564(00)00039-0View ArticlePubMedGoogle Scholar
- Behrmann I, Janzen C, Gerhartz C, Schmitz-Van de Leur H, Hermanns H, Heesel B, Graeve L, Horn F, Tav-ernier J, Heinrich C: A single STAT recruitment module in a chimeric cytokine receptor complex is sufficient for STAT activation. J Biol Chem. 1997, 272 (8): 5269-5274. 10.1074/jbc.272.8.5269View ArticlePubMedGoogle Scholar
- Yoshimura A: Negative regulation of cytokine signaling. Clin Rev Allergy Immunol. 2005, 28 (3): 205-220. 10.1385/CRIAI:28:3:205View ArticlePubMedGoogle Scholar
- Mason J, Morrison D, Basson M, Licht J: Sprouty proteins: multifaceted negative-feedback regulators of receptor tyrosine kinase signaling. Trends Cell Biol. 2006, 16: 45-54. 10.1016/j.tcb.2005.11.004View ArticlePubMedGoogle Scholar
- Zhang Y, Woude Vande G: Mig-6, signal transduction, stress response and cancer. Cell Cycle. 2007, 6 (5): 507-513.View ArticlePubMedGoogle 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.