Dynamic mathematics software, such as GeoGebra, is a kind of subject-specific digital tool used for enabling users to create mathematical objects and operate them dynamically and interactively, which is very suitable for mathematics teaching and learning at all school levels, especially at the secondary school level. However, limited research has focused on how multiple influencing factors of secondary school teachers' usage behavior of dynamic mathematics software work together. Based on the unified theory of acceptance and use of technology (UTAUT) model, combined with the concept of self-efficacy, this study proposed a conceptual model used to analyze the factors influencing secondary school teachers' usage behavior of dynamic mathematics software. Valid questionnaire data were provided by 393 secondary school mathematics teachers in the Hunan province of China and analyzed using a partial least squares structural equation modeling (PLS-SEM) method. The results showed that social influence, performance expectancy and effort expectancy significantly and positively affected secondary school teachers' behavioral intentions of dynamic mathematics software, and social influence was the greatest influential factor. In the meantime, facilitating conditions, self-efficacy and behavioral intention had significant and positive effects on secondary school teachers' usage behavior of dynamic mathematics software, and facilitating conditions were the greatest influential factor. Results from the multi-group analysis indicated that gender and teaching experience did not have significant moderating effects on all relationships in the dynamic mathematics software usage conceptual model. However, major had a moderating effect on the relationship between self-efficacy and usage behavior, as well as the relationship between behavioral intention and usage behavior. In addition, training had a moderating effect on the relationship between social influence and behavioral intention. This study has made a significant contribution to the development of a conceptual model that could be used to explore how multiple factors affected secondary school teachers' usage behavior of dynamic mathematics software. It also benefits the government, schools and universities in enhancing teachers' digital teaching competencies.
Citation: Zhiqiang Yuan, Xi Deng, Tianzi Ding, Jing Liu, Qi Tan. Factors influencing secondary school teachers' usage behavior of dynamic mathematics software: A partial least squares structural equation modeling (PLS-SEM) method[J]. Electronic Research Archive, 2023, 31(9): 5649-5684. doi: 10.3934/era.2023287
[1] | Djo Hasan . Purinergic P2X7R expressed on regulatory T cells potentially links molecular mimicry to autoimmune responses. AIMS Allergy and Immunology, 2024, 8(2): 80-123. doi: 10.3934/Allergy.2024006 |
[2] | Rosa Molfetta, Beatrice Zitti, Angela Santoni, Rossella Paolini . Ubiquitin and ubiquitin-like modifiers modulate NK cell-mediated recognition and killing of damaged cells. AIMS Allergy and Immunology, 2017, 1(4): 164-180. doi: 10.3934/Allergy.2017.4.164 |
[3] | Ken S. Rosenthal, Daniel H. Zimmerman . J-LEAPS vaccines elicit antigen specific Th1 responses by promoting maturation of type 1 dendritic cells (DC1). AIMS Allergy and Immunology, 2017, 1(2): 89-100. doi: 10.3934/Allergy.2017.2.89 |
[4] | Waleed M. Hussein, Phil M. Choi, Cheng Zhang, Emma Sierecki, Wayne Johnston, Zhongfan Jia, Michael J. Monteiro, Mariusz Skwarczynski, Yann Gambin, Istvan Toth . Investigating the affinity of poly tert-butyl acrylate toward Toll-Like Receptor 2. AIMS Allergy and Immunology, 2018, 2(3): 141-147. doi: 10.3934/Allergy.2018.3.141 |
[5] | Stefano Regis, Fabio Caliendo, Alessandra Dondero, Francesca Bellora, Beatrice Casu, Cristina Bottino, Roberta Castriconi . Main NK cell receptors and their ligands: regulation by microRNAs. AIMS Allergy and Immunology, 2018, 2(2): 98-112. doi: 10.3934/Allergy.2018.2.98 |
[6] | Mohamed El Yaqoubi, Mouad Lahyaoui, Yousra Seqqat, Taoufiq Saffaj, Bouchaib Ihssane, Nabil Saffaj, Rachid Mamouni, Fouad Ouazzani Chahdi, Fahad M Alshabrmi, Alaa Abdulaziz Alnahari, Saad M. Howladar, Ammar AL-Farga, Youssef Kandri Rodi . Pyridopyrazine derivatives as highly selective histamine H4 receptor antagonist for the treatment of atopic dermatitis: QSAR modeling and molecular docking studies. AIMS Allergy and Immunology, 2024, 8(4): 303-323. doi: 10.3934/Allergy.2024019 |
[7] | Daniil Shevyrev, Valeriy Tereshchenko, Olesya Manova, Vladimir kozlov . Homeostatic proliferation as a physiological process and a risk factor for autoimmune pathology. AIMS Allergy and Immunology, 2021, 5(1): 18-32. doi: 10.3934/Allergy.2021002 |
[8] | Tuba Karakurt, Nurhan Kasap, Kübra Aslan, Hayrunnisa Bekis Bozkurt, Fatma Bal Cetinkaya, Gizem Uslu, Zafer Bicakci, Filiz Özen, Özlem Cavkaytar, Ahmet Eken, Mustafa Arga . Suspected ALPS with clinical and laboratory findings: Three patients—three different diagnoses. AIMS Allergy and Immunology, 2023, 7(4): 304-312. doi: 10.3934/Allergy.2023020 |
[9] | Tamara Tuuminen, Jouni Lohi . Immunological and toxicological effects of bad indoor air to cause Dampness and Mold Hypersensitivity Syndrome. AIMS Allergy and Immunology, 2018, 2(4): 190-204. doi: 10.3934/Allergy.2018.4.190 |
[10] | Katarzyna Nazimek . The complex functions of microRNA-150 in allergy, autoimmunity and immune tolerance. AIMS Allergy and Immunology, 2021, 5(4): 195-221. doi: 10.3934/Allergy.2021016 |
Dynamic mathematics software, such as GeoGebra, is a kind of subject-specific digital tool used for enabling users to create mathematical objects and operate them dynamically and interactively, which is very suitable for mathematics teaching and learning at all school levels, especially at the secondary school level. However, limited research has focused on how multiple influencing factors of secondary school teachers' usage behavior of dynamic mathematics software work together. Based on the unified theory of acceptance and use of technology (UTAUT) model, combined with the concept of self-efficacy, this study proposed a conceptual model used to analyze the factors influencing secondary school teachers' usage behavior of dynamic mathematics software. Valid questionnaire data were provided by 393 secondary school mathematics teachers in the Hunan province of China and analyzed using a partial least squares structural equation modeling (PLS-SEM) method. The results showed that social influence, performance expectancy and effort expectancy significantly and positively affected secondary school teachers' behavioral intentions of dynamic mathematics software, and social influence was the greatest influential factor. In the meantime, facilitating conditions, self-efficacy and behavioral intention had significant and positive effects on secondary school teachers' usage behavior of dynamic mathematics software, and facilitating conditions were the greatest influential factor. Results from the multi-group analysis indicated that gender and teaching experience did not have significant moderating effects on all relationships in the dynamic mathematics software usage conceptual model. However, major had a moderating effect on the relationship between self-efficacy and usage behavior, as well as the relationship between behavioral intention and usage behavior. In addition, training had a moderating effect on the relationship between social influence and behavioral intention. This study has made a significant contribution to the development of a conceptual model that could be used to explore how multiple factors affected secondary school teachers' usage behavior of dynamic mathematics software. It also benefits the government, schools and universities in enhancing teachers' digital teaching competencies.
In a West Nile Virus (WNV) infection, there is a substantial self damage component and this is probably due to the way that the virus infects two cell populations differently. This difference, which involves an larger upregulation of MHC-1 sites on the surface of non-dividing infected cells over dividing infected cells, is critical in establishing a self damage or collateral damage response. In [1], we develop a macro level model of the nonlinear interactions between two critical signaling agents that mediate the interaction between these two sets of cell populations. In the case of WNV infection, the two signals are the MHC-1 upregulation level of the cell and the free WNV antigen level. This macro model allowed us to predict a hosts health response to varying levels of initial virus dose. Hence, we could begin to understand the oscillations in collateral damage and host health that lead to the survival data we see in WNV infections. In this work, we examine more general immunopathological interactions and develop a self affinity threshold model based on a first principles analysis of the immunosynapse. This model can then be used to gain insight into inappropriate autoimmunity events. In addition, this work serves as the background for two additional investigations: a general cytokine signal model and how this influences the self affinity threshold equation and a discussion of how the immunosynapse network model implied manifold structure allows us to interpret self damage as a topological defect.
The methods used here are quantitative and based on extensions of the Cubic Ternary Complex Model (CTC) model as discussed in [2,3] and [4]. The modifications are used to develop a simple model of an immunosynapse and how an autoimmune interaction can occur.
The CTC model is quite involved and to use it and also to modify it requires a good understanding of all of its details. Therefore, we will begin with a tutorial on how the CTC model works and afterwords, when we can discuss alterations and extensions, they will make more sense. The general assumptions of the CTC model are:
1. The receptor has two distinct binding sites: an external site accessible to agonists and antagonists and an internal site accessible to G-proteins.
2. External ligands ( agonists and antagonists ) and G-proteins exist in separate phases and do not encounter each other.
3. Receptors exist in two states with respect to their ability to activate G-proteins and initiate biological responses: active and inactive.
4. The interactions of external ligands, G-proteins, and receptor activation states are assumed to be governed by the laws of mass action.
5. All possible two-way and three-way interactions between external ligands, G-proteins, and receptor activation states (subject to assumption two) are potentially significant and are represented by coefficients in the model.
Letting A denote the ligand, G, the G-protein and R the receptor, we have 8 possible receptor species: Ri, Ra, ARi, ARa, RiG, RaG, ARiG and ARaG where the subscripts i and a indicate inactive and active, respectively.
These eight possibilities define six sets of interactions which, in this case can be drawn as a cube. See [2] for much better figures! We have shown the faces of this cube in the equations below. It helps to think of a sugar cube. If you hold it in your hand you can touch its six faces. Drawing the cube and labeling all the necessary equations makes for a very cluttered image. But we can draw the cube and label the corners as in
Each face of this cube has four corners which you can read off from Figure 1a but it is easier to see what is going on by drawing the cube in an exploded view and explicitly writing out all the interactions for each face.
The six faces with interaction pathways are shown in the equations below.
[ARi→ARa↑↑Ri→RaFrontFace:GFree][ARiG→ARaG↑↑RiG→RAGBackFace:GInvolved] |
[RiG→RaG↑↑Ri→RaBottomFace:AFree][ARiG→ARaG↑↑ARi→ARaTopFace:AInvolved] |
[ARi→ARiG↑↑Ri→RiGLeftFace:RiInvolved][ARa→ARaG↑↑Ra→RaGRightFace:RaInvolved] |
Note these faces then define 12 interaction edges:
Ri↔Ra, Ri↔ARi, Ra↔ARa, ARi↔ARaRa↔RaG, Ri↔RiG, ARi↔ARiG, ARa↔ARaGRAG↔ARaG, RiG↔ARiG, RiG↔RaG, ARiG↔ARaG |
where the double edged arrows indicate these reactions are reversible. An easy way to see this as a cube is to write all the equations as follows. We place the Back, Top and Front face above the Left, Bottom and Right faces for convenience of exposition.
[ARiG→ARaG↑↑RiG→RaGBack][ARiG→ARaG↑↑ARi→ARaTop][ARi→ARa↑↑Ri→RaFront][ARi→ARiG↑↑Ri→RiGLeft][RiG→RaG↑↑Ri→RaBottom][ARa→ARaG↑↑Ra→RaGRight] |
You can see this visually by taking the information in the columns above and assembling a cube as shown in Figure 1b. We want to make this notation more compact. The faces can be organized using mnemonics.
• The right face is Ra, A, G interactions.
• The left face is Ri, A, G interactions.
• The back face is Ri, Ra, A, G interactions, i.e. interactions with the G-protein.
• The front face is Ri, Ra, A interactions; i.e. no G-protein present
• The bottom face is Ri, Ra, G interactions; i.e. no ligand present
• The top face is Ri, Ra, A, G interactions; i.e. ligand is present
Let's organize this in column form with four systems: the R, A, G, the R, A, H, the R, B, G and the R, B, H for two G-proteins, G and H and two ligands A and B. Here we let Front = F, Back = Bk, Left = L, Right = R, Top = T and Bottom = Bo.
[RRa,B,HLRi,B,HBkRi,Ra,B,HFRi,Ra,BBoRi,Ra,HTRi,Ra,B,HB,H]...↔...[RRa,B,GLRi,B,GBkRi,Ra,B,GFRi,Ra,BBoRi,Ra,GTRi,Ra,B,GB,G]....↔..[RRa,A,GLRi,A,GBkRi,Ra,A,GFRi,Ra,ABoRi,Ra,GTRi,Ra,A,GA,G]...↔...[RRa,A,HLRi,A,HBkRi,Ra,A,HFRi,Ra,ABoRi,Ra,HTRi,Ra,A,HA,H] |
Connections between these subsystems are made with the ↔ markers. The R, B, H system connects to the R, B, G along the face with no G-protein interaction, the R, B, G then connects to the R, A, G system on the face with no ligand interaction and the R, A, G system interacts with the R, A, H system on the face with no G-protein interaction.
Activation is approached this way. Each face of these interaction cubes can be thought of as providing two paths to a given goal state. We assume thermodynamic closure so that the reaction rate constants for the reactions on the cube face are compatible. Here, Kact is association constant for the conversion of the inactive receptor Ri to the active form Ra; KG is the constant for the conversion of the inactivated form Ri to RiG and thermodynamic closure refers to the fact that the corresponding constant for the conversion of the activated form Ra to RaG corresponds to a scaling by β. Similarly, the conversion of the inactivated receptor Ri to the activated receptor Ra is scaled by α when ligand is present. The ligand pathway for A is also shown below and the comments about association constants are similar. Hence, the discussed pathways are as follows:
[RiGβKact↔RaG↑KG↑βKGRiKact↔Ra] [ARiαKact↔ARa↑KA↑αKARiKact↔Ra] |
For the inactive pathway γ denotes how the ligand A binding alters the equilibrium constant for the Ri to RiG pathway. We once again use thermodynamic closure to close the top pathway.
[RiGγKA↔ARiG↑KG↑γKGRiKA↔ARi] |
How do we interpret γ? It measures how ligand binding effects the RiG to ARiG pathway-i.e. how ligand binding effects the G-protein coupling and also it measures how the addition of the G-protein effects the other pathway, ARi to ARiG. Now, at this point, the association constants that connect the three remaining pathways to ARaG have not been specified. The left, front and bottom faces of the interaction cube have been totally specified, but the right, the top and the back each have two edges free. We have
[ARiG?→ARaG↑γKA↑?RiGβKact→RaGBack][ARiG?→ARaG↑γKG↑?ARiαKact→ARaTop][ARiαKact→ARa↑KA↑αKARiKact→RaFront][ARiγKA→ARiG↑KG↑γKGRiKA→RiGLeft][RiGβKact→RaG↑KG↑βKGRiKact→RaBottom][ARa?→ARaG↑αKA↑?RaβKG→RaGRight] |
On the back, right and top face, there are two edges each unspecified. Note thermodynamic closure tells us that if we fix one edge on one of these faces, closure will fill in other edges. For example, specifying the RaG to ARaG path will fill in the rest of the right face and the rest of the back face. That will fill in the ARiG to ARaG path in the top face too and then closure fills in the rest of the top face. So specifying one edge, fills in the full cube of interactions. Hence, there are just the three free edges RaG to ARaG, ARa to ARaG and ARiG to ARaG; these are the three edges that lead to ARaG. For example, if we set RaG ↔ ARaG to have constant αγKA, by closure we fill in all pathways to find
[ARiGαβKact→ARaG↑γKA↑αγKARiGβKact→RaGBack][ARiGαβKact→ARaG↑γKG↑γβKGARiαKact→ARaTop][ARiαKact→ARa↑KA↑αKARiKact→RaFront][ARiγKA→ARiG↑KG↑γKGRiKA→RiGLeft][RiGβKact→RaG↑KG↑βKGRiKact→RaBottom][ARaγβKG→ARaG↑αKA↑αγKARaβKG→RaGRight] |
This gives the pathways
[RaG0.6in↑ αγKAARiGαβKact→ARaG0.6in↑ βγKGARa] |
To make things as general as possible, add an additional reaction parameter δ to these pathways:
[RaG0.6in↑ δαγKAARiGδαβKact→ARaG0.6in↑ δβγKGARa] |
Normally, the reactions Ra + G ↔ RaG (constant βKG) and ARa + G ↔ ARaG (constant γβKG) differ by a factor γ. If the addition of the ligand enhances the G-protein coupling, then we can model that by changing the factor to δγ where δ > 1. If the G-protein coupling is diminished by the ligand, we can use δ < 1 to indicate that. Finally, if δ = 1, there is no advantage or disadvantage to the G - protein coupling due to using the ligand. This moves us to the following diagram
[ARiGδαβKact→ARaG↑γKA↑δαγKARiGβKact→RaGBack][ARiGδαβKact→ARaG↑γKG↑γδβKGARiαKact→ARaTop][ARiαKact→ARa↑KA↑αKARiKact→RaFront][ARiγKA→ARiG↑KG↑γKGRiKA→RiGLeft][RiGβKact→RaG↑KG↑βKGRiKact→RaBottom][ARaδγβKG→ARaG↑αKA↑δαγKARaβKG→RaGRight] |
and in it, we can see how the factor δ shapes the interactions.
Now let's look at how to model affinity in these sorts of interactions.
We say the unbound receptor cloud here is the ensemble consisting of the four members Ri, Ra, RiG and RaG. There is then the ligand bound receptor cloud having members ARi, ARa, ARiG and ARaG. Thus, in the cubic ternary complex model, there are four ligand bound types in the ligand cloud. The total concentration of bound species is then RB
[RB]=[ARi]+[ARa]+[ARiG]+[ARiG]. |
The total concentration of unbound ligand is then denoted by RU and we know the total amount of ligand RT satisfies RT = RB + RU. Let KLA be the association constant for the reaction
A+RUKLA→RB. |
In other words, KLA is the equilibrium association constant for the overall conversion of unbound species into bound species. Standard equilibrium kinetics then says
[RiGβKact→RaG↑KG↑βKGRiKact→RaReceptors[RU]]KLA×[A]×[RU]=[RB]⇒[ARiGδαβKact→ARaG↑γKG↑δγβKGARiαKact→ARaLigands[RB]] |
There are four separate association constants here: KA, αKA, γKA and δαγKA. In [3], we find
KLA=w1KA+w2(αKA)+w3(γKA)+w4(δαγKA) |
where w1, w2, w3 and w4 are the relative frequencies of the various receptor species. Thus,
w1=[Ri][Ri]+[Ra]+[RiG]+[RaG], w2=[Ra][Ri]+[Ra]+[RiG]+[RaG]w3=[RiG][Ri]+[Ra]+[RiG]+[RaG], w4=[RaG][Ri]+[Ra]+[RiG]+[RaG] |
Using these weights, we find
KLA=KA([Ri]+α[Ra]+γ[RiG]+δαγ[RaG])[Ri]+[Ra]+[RiG]+[RaG] |
But Ra = KactRi, RiG = KGRiG, RaG = βKGKactRiG and so after rearranging we have
KLA=KA[Ri](1+αKact+γKG[G]+δαβγKGKact[Ri][G])[Ri](1+Kact+KG[G]+KGKact[G])=KA1+αKact+γKG[G]+δαβγKGKact[G]1+Kact+KG[G]+βKGKact[G] |
We call KLA the apparent affinity to distinguish it from the association constant KA which governs the addition of ligand from the receptor cloud. Note if there is no G-protein, we have
KLA=KA[Ri](1+αKact+γKG[G]+δαβγKGKact[Ri][G])[Ri](1+Kact+KG[G]+KGKact[G])=KA1+αKact+γKG[G]+δαβγKGKact[G]1+Kact+KG[G]+βKGKact[G] |
Note the apparent affinity now is a weighted sum of the constants for the RiG → ARiG and RaG → ARaG pathways in the ligand cloud. However, it is important to note that the weights themselves depend only on receptor cloud constants.
Eventually, we want to define avidity in an immunological setting, but let's start with a notion of efficacy. We will follow the exposition in [4] with some modifications. In general, a compound that binds to a receptor and after being bound induces a physiological response is called an agonist. The ability of a ligand to bind to a receptor is quantified by its apparent affinity which we have just discussed. The efficacy of the ligand is the measure of how effective it is at eliciting a physiological response. We know there are antagonists which bind to receptors but generate no physiological response and so we suspect that affinity and efficacy are independent parameters. Of course, the agent must be able to bind the receptor in order to talk about efficacy, but as long as the affinity is not zero, the affinity and efficacy of a ligand - receptor pair is not correlated.
Before the ligand A is introduced to the CTC model system, there is only one active species present; RaG. Following [4], we denote its concentration by [rag]. After the ligand is introduced, there are two active species, ARaG and RaG. The generation of active species can be thought of as a stimulus, S and so efficacy should be somehow related to the net stimulus ΔS = [ARaG] + [RaG] - [rag]. This means we call a ligand efficacious only if it changes the production of active species beyond what is observed in the absence of ligand. This kind of efficacy is called internal efficacy because it arises from the internal connections of the model. There is an additional source of efficacy we call external efficacy which comes from external connections. A useful diagram for this is as follows:
[RiGβKact→RaG↑KG↑βKGRiKact→RaReceptorCloud]KLA⇒[ARiGδαβKact→ARaG↑γKGInternalEfficacy↑δγβKGARiαKact→ARaLigandCloud]ExternalEfficacy⇒ |
We can derive an expression for internal efficacy in the CTC model. Recall an efficacious response for a ligand requires
[ARaG]+[RaG]TotalReceptorPostLigand−[rag]TotalReceptorPriorLigand>0 |
The total receptor concentration before the ligand is introduced is
RbeforeT=[ri]+[ra]+[rig]+[rag]=[ri]+Kact[ri]+KG[g][ri]+KactβKG[g][ri]=[ri](1+Kact+KG[g]+KactβKG[g])=[ri]a |
where a is the long term in the parenthesis. After the ligand is introduced, we have
RafterT=[Ri]+[Ra]+[RiG]+[RaG][ARi]+[ARa]+[ARiG]+[ARaG]=[Ri]+Kact[Ri]+KG[G][Ri]+KactβKG[G][Ri]+KA[A][Ri]+αKA[A]Kact[Ri]γKA[A]KG[G][Ri]+αγδKA[A]KactβKG[G][Ri]=[Ri](1+Kact+KG[G]+KactβKG[G]+)[Ri](KA[A]+αKA[A]Kact+γKA[A]KG[G]+αγδKA[A]KactβKG[G])=[Ri](a*+b) |
where a* = 1 + Kact + KG[G] + KactβKG[G] and b is the other term. Next, we note
[ARaG]+[RaG]=αγδKactβ[G]KA[A][Ri]+KactβKG[G][Ri]=KactβKG[G][Ri](αγδKA[A]+1) |
and [rag] = KactβKG[g][ri]. Thus,
ΔS=KactβKG[G][Ri](αγδKA[A]+1)RafterT−KactβKG[g][ri]RbeforeT=KactβKG[G](αγδKA[A]+1)a*+b−KactβKG[g]a=βKactKG([G](αγδKA[A]+1)a*+b−[g]a)=βKactKG[G](αγδKA[A]+1)a−[g](a*+b)(a*+b)a |
We want to find out when the efficacy is positive, so we only need look at the numerator of this large fraction. Plugging in for a*, a and b, and letting the numerator be N, we find
N=[G](αγδKA[A]+1)(1+Kact+KG[g]+KactβKG[g])−[g](1+Kact+KG[G]+KactβKG[G])−[g](KA[A]+αKA[A]Kact+γKA[A]KG[G]+αγδKA[A]KactβKG[G]) |
We can make a lot more sense out of this equation by starting to use vectors and dot products. We can rewrite N as
N=(αγδKA[A]+1)[1111]·[G][1KactKG[g]βKactKG[g]]−[1111]·[g][1KactKG[G]βKactKG[G]]−[1111]·[A]KA[g][1αKactγKG[G]αγβδKactKG[G]] |
Let
Zg=[1,Kact,KG[g],βKactKG[g]]', ZG=[1,Kact,KG[G],βKactKG[G]]'X=[1,α,γ,αγδ]' |
and then letting I be <1,1,1,1>', we can rewrite N as
N=(αγδKA[A]+1)[G]〈I,Zg〉−[g]〈I,ZG〉−KA[A][g]〈X,ZG〉=(αγδKA[A])[G]〈I,Zg〉+[G]〈I,Zg〉−[g]〈I,ZG〉−KA[A][g]〈X,ZG〉 |
This simplifies to
N=[1111]·[[G]−[g]Kact([G]−[g])00]+KA[A][1111]·[αγδ[G]−[g]αγδKact[G]−αKact[g]αγδKG[g][G]−γKG[g][G]αβγδKactKG[g][G]−αβγδKactKG[g][G]]=[1111]·[[G]−[g]Kact([G]−[g])00]+KA[A][1111]·[αγδ[G]−[g]αγδKact[G]−αKact[g]αγδKG[g][G]−γKG[g][G]0]=[1111]·[[G]−[g]Kact([G]−[g])00]+KA[A](αγδ[1111]·[G][1KactKG[g]0]−[1αγ0]·[g][1KactKG[G]0]) |
Now define the vectors ΔG, U and V as
ΔG=[[G]−[g]Kact([G]−[g]),0,0,]', U=αγδ[1,1,1,0]', V=[1,α,γ,0]' |
Then, N has the form
N=〈I,ΔG〉+KA[A]([G]〈U,Zg〉−[g]〈V,ZG〉) |
In [4], to simplify the analysis, it is assumed [G] = [g]. In that case, Zg = ZG and we have ΔG = 0 giving N = KA[A][G] < U - V, ZG >. The equation we derive is a slight extension of the one in [4] as we want to apply this more general form to our discussion of the immunosynapse. Since we want N > 0, we can take the more general equation and break it apart some. Look at the top two components of each vector dot product and the bottom two components separately and assume [G] = ξ[g]. Then, we have
N=(ξ−1)(1+Kact)[g]+KA[A][g](ξ[αγδαγδ]−[1α])·[1Kact]+KA[A][g](ξ[αγδαβδ]·[KG[g]0]−[γ0]·[ξKG[g]0]) |
Thus, we have
N=(ξ−1)(1+Kact)[g]+KA[A][g]((ξαγδ−1)+αKact(ξγδ−1))+KG[A][g]2γKAξ(αδ−1) |
Note if ξ = 1, we want (cancelling out the common [g]),
N=KA[A]((αγδ−1)+αKact(γδ−1))+KG[A][g]γKA(αδ−1)>0. |
We can write this in dot product form as N = KA[A][αγδ - 1, αγδ - α, αγδ - γ, 0]' · [1, Kact, KG[g], 0]' which is the form used in [4]. We can put our modified equation into this form as well. Again, we let [G] = ξ[g]. Then, we have
N=(ξ−1)(1+Kact)[g]+KA[A][g]〈[ξαγδξαγδξαγδ0],[1KactKG[g]0]〉−KA[A][g]〈[1αξγ0],[1KactKG[g]0]〉=(ξ−1)(1+Kact)[g]+KA[A]ξ[g](〈[αγδαγδαγδ0],[1KactKG[g]0]〉−〈[ξ−1ξ−1αγ0],[1KactKG[g]0]〉) |
If we let U = αγδ[1,1,1,0]' (as defined before), WG and Wg as WG = [1, Kact, KG[G], 0]' and Wg = [1, Kact, KG[g], 0]', the new vector Vξ = [ξ−1, ξ−1α, γ, 0]', we can rewrite this as
N=(ξ−1)(1+Kact)[g]+KA[A]ξ[g]〈U−Vξ,Wg〉 |
We now use the cubic ternary complex model to model the immunosynapse which forms when a T Cell binds to a peptide and MHC-I complex, pMHCI. We define the vector U - Vξ to be the internal efficacy vector, eI. This is based on the components that determine interactions in the ligand cloud. Then, we define Wg to be the receptor cloud vector which is based on the interactions that determine the relationships between the receptor species. Then, the dot product < U - Vξ, Wg > is called the realized efficacy value, eR. Note, the Cauchy - Schwartz inequality tells us that
eR=〈U−Vξ,Wg〉=‖U−Vξ‖‖Wg‖cos(φ) |
where φ is the angle between ligand space and receptor space. The external efficacy is then eE given by
eE=G-ProteinChange+ExternalScalingFactor×RealizedEfficacy=(ξ−1)(1+Kact)[g]+KA[A]ξ[g]‖U−Vξ‖‖Wg‖cos(φ) |
Note the concentration [g] can be factored out to give
eE=[g]((ξ−1)(1+Kact)+KA[A]ξ‖U−Vξ‖‖Wg‖cos(φ)) |
Since [g] is a scaling factor, the algebraic sign of the change in stimulus is determined by the sign of the discriminant 𝒟 given by
𝒟=(ξ−1)(1+Kact)+KA[A]ξ‖U−Vξ‖‖Wg‖cos(φ) |
Recall, previously we found N ( the numerator which determined algebraic sign ) was given by
N=[g]((ξ−1)(1+Kact)+KA[A]ξ〈U−Vξ,Wg〉) |
Hence, N is the same as the external efficacy.
However, to understand any of this, we need to fully understand the meaning the the parameters α, δ and γ; let's review the discussion of [2]. The term α is concerned with the effect the binding of the ligand A has on the activation of the receptor. Hence, α > 1 tells us that the ligand binding increases the production of the activated form on the receptor Ra. The term γ models the effect the binding of the ligand A has on the coupling of the G-protein G to the receptor. If γ > 1, then the binding of the ligand helps later coupling to the G-protein. Finally, the parameter δ is the most interesting. There are several cases to consider; in all of them, if δ > 1, this means that the coupling of G-protein is enhanced by the addition of the ligand, if δ = 1, the addition of the ligand did not effect the G-protein binding and if δ < 1, the coupling of G-protein is reduced.
We need to define our terms. Here the inactivated receptor Ri is TcR prior to adding the auxiliary proteins needed to prepare it for binding to the pMHC; the activated receptor Ra is the TcR after these preparations are completed. The ligand A is the pMHC and the G-protein G here is a collection of proteins that provide a scaffolding for the TcR-pMHC interaction.
• The effect a bound ligand has on whether or not a receptor will couple with a G-protein can be different if the receptor is inactive or active. Hence, the pathways Ra + G ↔ RaG and ARa + G ↔ ARaG can be independent from each other. Without the δ interaction, the association constants for these two pathways are βKG versus γβKG, but the δ option changes the second pathway's constant to δγβKG. Hence, if δ goes up here, there is enhanced coupling of the activated TcR and the scaffolding.
• We can also look at the pathways Ri + G ↔ RiG ( association constant KG) and ARi + G ↔ ARiG (association constant γKG). The δ option changes the second pathway's constant to δγKG. We interpret the δ choices the same way. If δ increases here, there is an enhancement of the coupling of the scaffolding between the inactivated TcR and the scaffolding.
• The pathways RiG + A ↔ ARiG (association constant γKA) and RaG + A ↔ ARaG (association constant αγKA). The δ option changes the second pathway's constant to δαγKA. Thus, the binding of pMHC to the inactivated TcR and scaffolding complex is enhanced.
• The pathways RiG ↔ RaG (association constant βKact) and ARiG ↔ ARaG (association constant αβKact). The δ option changes the second pathway's constant to δαβKact. This means that the coupling between the inactivated TcR and the activated TcR is enhanced if δ increases.
• The pathways ARi ↔ ARa ( association constant αKact) and ARiG ↔ ARaG (association constant αβKact). The δ option changes the second pathway's constant to δαβKact. Here, if δ increases, the activation of the receptor TcR enhances the binding of pMHC to the receptor.
• The pathways ARi + G ↔ ARiG (association constant γKG) and ARa + G ↔ ARaG (association constant γβKG). The δ option changes the second pathway's constant to δαβKact. Now, if δ increases, the interpretation is that the addition of the scaffolding to the inactive TcR binding with pMHC is enhanced.
In all cases, the parameter δ has to do with a synergism between the ligand and G-protein binding. Indeed, increasing δ in several cases can potentially encourage autoimmune interactions.
Now let's add a bit more detail for our model of TcR and pMHC binding in the context of the adaptive immune system. We assume the ligand A is the pMHC complex, the receptor R is the TcR complex and there are accessory proteins P1 and P2 to determine a scaffolding that is build up to facilitate the TcR and pMHC interaction that forms the immunosynapse. We have lumped this activity into the G-protein actions of the CTC model. We will therefore the G label with the new label Sc. To save space, we will let T be the TcR receptor and P, the pMHC complex. Let's redo our fundamental equations for this context.
[PTiScαβKact→PTaSc↑γKP↑αγKPTiScβKact→TaScBack][PTiScαβKact→PTaSc↑γKSc↑γβKScPTiαKact→PTaTop][PTiαKact→PTa↑KP↑αKPTiKact→TaFront][PTiγKP→PTiSc↑KSc↑γKScTiKP→TiScLeft][TiScβKact→TaSc↑KSc↑βKScTiKact→TaBottom][PTaγβKSc→PTaSc↑αKP↑αγKPTaβKSc→TaScRight] |
This gives the pathways
[TaSc↑αγKPPTiScαβKact→PTaSc↑βγKScPTa] |
which we make as general as possible by adding the δ synergy parameter as before:
[TaSc↑δαγKPPTiScδαβKact→PTaSc↑δβγKScPTa] |
Normally, the reactions Ta + Sc ↔ TaSc (constant βKSc) and PTa + Sc ↔ PTaSc (constant γβKSc) differ by a factor γ. If the pMHC binding enhances the scaffold formation, then we can model that by changing the factor to δγ where δ > 1. If the scaffold coupling is diminished by the ligand, we can use δ < 1 to indicate that. Finally, if δ = 1, there is no advantage or disadvantage to the scaffold coupling due to the pMHC binding. This moves us to the following diagram
[PTiScδαβKact→PTaSc↑γKP↑δαγKPTiScβKact→TaScBack][PTiScδαβKact→PTaSc↑γKSc↑δγβKScPTiαKact→PTaTop][PTiαKact→PTa↑KP↑αKPTiKact→TaFront][PTiγKP→PTiSc↑KSc↑γKScTiKP→TiScLeft][TiScβKact→TaSc↑KSc↑βKScTiKact→TaBottom][PTaδγβKSc→PTaSc↑αKP↑δαγKPTaβKSc→TaScRight] |
and in it, we can see how the factor δ shapes the interactions. The external efficacy eE here is given by
eE=Sc-Change+ExternalScalingFactor×RealizedEfficacy=(ξ−1)(1+Kact)[sc]+KP[P]ξ[sc][αγδ−ξ−1αγδ−αξ−1αγδ−γ0]·[1KactKSc[sc]0]=[sc]((ξ−1)(1+Kact)+KP[P]ξ〈U−Vξ,Wsc〉) |
where [sc] denotes the concentration of the scaffolding system Sc prior to the pMHC binding to the receptor cloud.
We can calculate the rate of change of the external efficacy eE with respect to a given parameter p in the usual way:
∂eE∂p=∂∂p{[sc]((ξ−1)(1+Kact)+KP[P]ξ〈U−Vξ,Wsc〉)}=∂[sc]∂p{(ξ−1)(1+Kact)+KP[P]ξ〈U−Vξ,Wsc〉}+[sc]KP[P]ξ∂∂p(〈U−Vξ,Wsc〉) |
where p is any parameter that is interesting to us. Now if [sc] does not depend on the parameter p, this becomes
∂eE∂p=[sc]KA[P]ξ〈∂(U−Vξ)∂p,Wsc〉 |
We can calculate these partials. Since,
U−Vξ=[αγδ−ξ−1,αγδ−αξ−1,αγδ−γ,0]' |
we find
(U−Vξ)α=[γδ,γδ−ξ−1,γδ,0]', (U−Vξ)δ=[αγ,αγ,αγ,0]'(U−Vξ)γ=[αδ,αδ,αδ−1,0]' |
Thus,
∂eE∂α=[sc]KP[P]ξ〈∂(U−Vξ)∂α,Wsc〉=[sc]KP[P]ξ[γδγδ−ξ−1γδ0]·[1KactKSc[g]0]=[sc]KP[P]ξ(γδ+(γδ−ξ−1)Kact+γδKSc[sc]) |
and
∂eE∂δ=[sc]KA[P]ξ〈∂(U−Vξ)∂δ,Wsc〉=[sc]KA[P]ξ[αγαγαδ0]·[1KactKSc[g]0]=[sc]KA[P]ξ(αγ+αγKact+αγKG[sc]) |
and
∂eE∂γ=[sc]KP[P]ξ〈∂(U−Vξ)∂γ,Wsc〉=[sc]KP[P]ξ[αδαδαδ−10]·[1KactKSc[sc]0]=[sc]KP[P]ξ(αδ+αδKact+(αδ−1)KSc[sc]) |
In the three dimensional ligand association constant space (α, δ, γ), we see
∇eE=[sc]KP[P]ξ(1+Kact+KSc[sc])[γδαγαδ]−[sc]KP[P][Kact0ξKG[sc]] |
which implies the <∇eE, u> gives us the directional derivative Du(eE) in any unit vector direction: i.e. ||u|| = 1 and for u = < Δα, Δδ, Δγ >', this implies Δα2 + Δδ2 + Δγ2 = 1. To get a sense of what this implies, look at a unit change in γ only and also assume αδ = 1. Then
• Δγ = 1, Δα = 0 and Δδ = 0. We have
Du(eE)=〈[sc]KP[P]ξ(1+Kact+KSc[sc])[γδαγαδ]−[sc]KP[P][Kact0ξKSc[sc]],[0,0,1]〉=[sc]KP[P]ξ(1+Kact)>0. |
We know γ measures the change in the Ri → RiG pathway or how the ligand A binding effects G-protein coupling. In the context of an immunosynapse model, the ligand A is the pMHC containing an epitope on the cell's surface and the receptor Ri is the inactivated TcR cell. Hence, a one dimensional alteration of γ only in the ligand association space induces an increase in the external efficacy eE. The positive increase in γ thus represents an increase in G - protein coupling to the inactive receptor due to the presence of the ligand A or pMHC. We can think of this as saying that as the TcR approaches the pMHC, there is a positive feedback loop which induces structural proteins to enhance the binding of the TcR to the pMHC. A similar analysis works for the case of the activated TcR also. Then an increase in γ enhances the G-protein (scaffolding) coupling between the pMHC and the activated TcR via another positive feedback loop.
• Δγ = -1, Δα = 0 and Δδ = 0. We have
Du(eE)=〈[sc]KP[P]ξ(1+Kact+KSc[sc])[γδαγαδ]−[sc]KP[P][Kact0ξKSc[sc]],[0,0,−1]〉=−[sc]KP[P]ξ(1+Kact)<0. |
implying a decrease in the external efficacy eE.
In general, if the directional derivative in ligand associative space gives rise to an increase in external efficacy there is an increased probability of a successful TCR and pMHC binding to create the full immunological synapse.
The modifications of the CTC model above now give us the tools to discuss how autoimmune events might occur.
There is much that is not known about T cell recognition algorithms in vivo. A current review of salient issues can be found in [5] and we will rely heavily on what is said in that review in what follows. Huppa et. al. distinguish three levels of complexity in the TcR-pMHC interaction: the kinetics and a structural basis for TcR-pMHC engagement, the TcR interactions with the structural protein scaffolding we have mentioned. His paper focuses on the TcR-CD3 complex and antigen recognition and the immunological synapse. New techniques have surfaced to help answer these questions. There has been a lot more crystallographic analysis of various TcR-pMHC complexes which has helped. Surface Plasmon Resonance (SPR) has allowed us to study these engagements in solution. The problem, of course, is that these bindings are constrained to be essentially 2D due to the ligand cloud and receptor cloud interaction at a planar interface. So there is a significant change in how many degrees of freedom are available which alters the measurement of what we call external efficacy. Other tools are total internal reflection (TIRF) microscopy, Förster resonance energy transfer (FRET) microscopy, superresolution fluorescence and so on. All have their problems and so there is no simple analytical tool to resolve our questions. However, these tools are letting us make a lot of progress.
It is now clear the external affinity between the TcR and a pMHC is very low ( 1-100 μM) which is caused by a slow association and a fast disassociation rate. Also, peptides in the MHC groove which differ by only a slight amount in peptide composition from an antigen can produce very different T cell responses. This is very strange and hard to explain; yet, that is our job! There is evidence that TcR engagement of pMHC binding occurs in two stages: the first involves MHC binding and the second handles peptide scanning and binding. In the CDR 1-3 regions within the TcR, it seems to work this way: the CDR 1 and CDR 2 loops bind mostly to the MHC while the more diverse and flexible CDR 3 loop engages the MHC embedded peptide. TcRs engage the MHC primarily on a diagonal. Vα fits with the MHC helix α2 for class I and the β1 for class 2 and the peptide N-terminus where as Vβ contacts the MHC α1 helix and the peptide C-terminus. Hence, the germline encoded CDR 1 and CDR 2 contact the α helices of the MHC and the CDR 3 regions bind to the central and most diverse part of the peptide in the MHC groove. There is a wide range of docking angles ( our φ in the external efficacy equation comes to mind!) of ±100 deg. But, it appears, there is a restricted range of angles that promote signaling.
Hence, there appears to be a germline bias for TcR recognition which suggest evolutionary pressures shape V-gene segments within the TcRα and TcRβ gene loci that allow them to bind to MHC regardless of the peptide cargo. The codon hypothesis posits that a given TcR equipped with a set of CDR 1 and CDR 2 loops can bind more than one MHC molecule. Binding can be achieved at a wide range of docking angles for MHC molecules specifically. The cost of this ability is probably a lack of binding strength which fits well with the generally low external affinity measured. Thus, there is significant TcR cross-reactivity. In more detail, each Vα and Vβ interacts with diverse MHC surface residues on top of the helices of different MHC molecules. Suppose this happens for three MHC molecules called MA, MB and MC. Our TcR has two keys: the CDR 1 and CDR 2 loops. Let MA have two structural regions or locks, CA1 and CA2, MB have locks CB1 and CB2 and finally MC have CC1 and CC2. Then we can get the bindings CDR 1 + CA1 and CDR 2 + CA2, CDR 1 + CB1 and CDR 2 + CB2 and CDR 1 + CC1 and CDR 2 + CC2. These interactions are then influenced by the more flexible CDR 3 interaction with the peptide in the MHC groove. One of these interactions is the low binding energy solution. Also, these interactions are influenced by the binding angle ( our φ). The TcR-MHC interaction is very short lived and if this reaction had high affinity, then it would not be very helpful in helping discriminate between subtle differences in peptide structure. The human TcR repertoire is about 108 TcRs at any given time and without cross-reactivity would not be enough to cover enough peptide variability for our protection. However, there is a price. If each TcR can interact profitably with multiple peptide loaded MHC cradles, there is a distinct possibility that a TcR could interact with a self MHC leading to autoimmunity issues.
As discussed in [6], it is now understood that α, β T cells have a broadly sensory role in detecting perturbations from self and normalcy.
Because the T cells must not only detect foreign antigens but also act on that information very rapidly, another important issue is what the threshold for an irreversible action will be.
The models for TcR mediated T cell activation have evolved over the years. First there was the monomer/ co-receptor model in which a single pMHC engages a TcR and also almost simultaneously engages CD4 or CD8 co-receptors. However, in general even high concentrations of monomeric pMHC fail to stimulate most T cells. Hence, the dimer of dimers model was developed. It was suggested that MHC dimerization could promote TcR dimerization and a following activation. But there is data that shows that even single pMHC complexes can trigger some T cells, so the pseudodimer model was developed. Here, a single pMHC could join up with an equivalent pMHC and with cross linking by CD4 or CD8 would provide the base unit for the T cell activation. The linking up with the CD4 or CD8 is why the model is called a pseudodimer, of course, as it is not really a dimer at all. Observations have shown that self or endogenous pMHC is strongly recruited to the immunological synapse by interactions with the TcR. The pseudodimer model was then modified in [6] to give a model of serial engagement. If an α or β TcR receptor on a CD4 T cell, TCR1, interacts with a pMHC having a foreign peptide or agonist pMHC complex, this interaction allows the recruitment of a second TcR, TCR2, through its associated CD4 molecule. TcR2 then interacts with a self pMHC. So the chain is TCR1 + agonist + TCR2 + self at this point. It is postulated the second TcR2 + self pMHC is very short lived. One possibility is that this short lived second complex stabilizes the first TCR1 + agonist interaction and allows the phosphorylation of TcR1's CD3 molecules. The second possibility is that the TcR1 + agonist acts like a catalyst and facilitates the addition of multiple short lived TcR + self complexes. Each additional TcR + self complex would be short lived but would cause additional CD3 phosphorylations on the T Cell. Hence, the chain would be
TcR1+agonist+TcR2+self A+TcR3+self B+… |
where the additional TcRs need not recruit the same self pMHC and so we call them self A and self B here. Hence, T cells can use endogenous pMHC-TcR interactions to augment T cell responsiveness. In [6], this discussion of serial engagement is restricted to CD4 T cell which works with MHC class II and it is said there is less evidence for such a recruitment of endogenous pMHC for the CD8 T cells. However, this could well be that the experimental means to see this recruitment is currently not known. As is said in [6], there is every reason to believe the activation sequence for CD4 and CD8 T cells is similar, so from now on, we will assume we can model the TcR and pMHC activation sequence in this way for both a CD4 and a CD8 T cell.
We characterize the recruited endogenous pMHC recruited in this process as follows. The initial foreign peptide TcR pMHC interaction event causes a movement of self pMHCs, possibly similar by some measure, to the foreign one on the cell surface to the binding site. It is hard to say what we mean by similar at this point. Does it mean the peptides lying in the MHC cradles are similar in some fashion so that the binding of the CDR 3 flexible loop to the peptide in the cradle is similar? Or is it possible for quite dissimilar peptides to initiate similar external affinity values? Or can any self pMHC be recruited to help strength the original TcR and agonist pMHC?
Following [5], we note that T cell antigen recognition and the formation of the immunological synapse (IS) are intimately connected.
The term “synapse” is now applied to any bimembranous junction between T cells and APCs...The common denominator of all immunological synapses is probably that they promote, with varying complexity, one -on-one conversations between two conjugated cells through engagement of cell-specific receptor-ligand pairs and in most cases, also through the directed secretion of cytokines and perforins.
It is known the formation and maintenance of the IS is dependent on the TcR and APC engagement and the TcR mediated signaling. The typical look of this interaction in terms of what we have been calling the scaffolding is shown in Figure 2.
There are many proteins that are involved in the construction of what we call the scaffolding Sc. These scaffolding proteins accumulate in the mature IS in a nested sequence of proteins which is called the bulls eye as shown in Figure 3. As the serial interactions we have listed as
TcR1+agonist+TcR2+self A+TcR3+self B+… |
continue, the bulls eye forms in the maturing IS. It must be noted though that this kind of synaptic binding is very different from binding in solution as the system is open instead of closed with energy moving in and out all the time. As [5] says
Molecular partners are preoriented and tied to flexible membrane surfaces, which can adjust their morphology in response to cell stimuli and acting forces. This should ... increase the on - rate of binding ... however, only if the two participating membranes are properly aligned. Moreover, molecular size could ... affect binding within the crowded confines of the immunological synapse.
For example, the protein CD45 is rather large and becomes pushed out of the central bulls eye because of its size. That is why CD45 is shown in the outer most rectangle in our bulls eye drawing. So it can not access the central part of the IS and so the equilibrium shifts towards the TcR-CD3 complex and its phosphorylation. During the IS formation, there is a huge flow of actin which imposes shear forces on the IS which destabilized weaker TcR-pMHC interactions. Hence, a more stable TcR-pMHC interaction is more probable if it followed the serial or chain interaction scheme which uses a stream of frequent bindings. It is theorized that these forces are an important part of the IS formation. Experiments have shown that the TcR-pMHC chains have a more vigorous response once it is attached to a glass supported lipid bilayer. Before such an attachment, TcR movement is in all three dimensions, but once there is an attachment, one degree of freedom in the movement is lost with a subsequent increase in the external affinity. Still, much is not known for sure. It is possible that LFA-1 when it binds to ICAM1 protects it from the shear stresses we mentioned earlier. And the CD2-CD48 interaction might help setup the T cell and APC membranes for an optimal TcR-pMHC interaction.
Also we should remember the T cell is like a large bag whose surface is studded with many protein complexes. The distribution of TcRs on the surface of a T cell appears to be organized as clusters. From [5], we have
a single quantum-dot labeled pMHC within the T cell APC interface caused the formation of TcR microclusters, in which hundreds of TcRs surrounded the antigenic pMHC.
We can now begin to use the affinity model developed above to do quantitative experiments.
It is illuminating to look at what happens during a binding event chain. In our model, the original TcR1 + agonist interaction has an external efficacy of
e1E=(ξ1−1)(1+K1act)[(sc)1]+K1P[P1]ξ1[(sc)1][α1γ1δ1−ξ−11α1γ1δ1−α1ξ−11α1γ1δ1−γ10]·[1K1actK1Sc[(sc)1]0]=(ξ1−1)(1+K1act)[(sc)1]+K1P[P1]ξ1[(sc)1]<U1−V1ξ1,W(sc1)> |
We would then have similar expressions for the external efficacy of the TcR2, TcR3 and so forth additional complexes that are being recruited catalytically. This suggests the external efficacy of this growing T cell activation due to a chain of N recruitments should be
𝒜=e1E+e2E+e3E+...+eNE |
where 𝒜 is the avidity of the interaction. The original e1E, although low, is enough to begin this complicated recruitment chain. The self pMHCs brought into this process could even be part of the thymus clonal expansion that results in a huge increase in the T cell population that targets the agonist.
Next, note this definition of avidity does not take into account time issues though, so next we will talk about that a bit. Let's assume a pMHC corresponding to an antigen binds to a TcR. There is a critical external efficacy which must be exceeded in order to the pMHC and TcR connection to take hold. Call this threshold e1*E where the superscript 1 refers to the fact that this peptide anchors the chain of additional pMHC interactions we have talked about already. Hence, for a valid interaction, we must have e1E>e1*E and moreover, the binding of the pMHC much last longer than a critical time interval we will label Δ1*t. If the time interval that this interaction occurs is denoted by Δ1t, we must have Δ1t > Δ1*t for a successful interaction. Each additional pMHC has an associated critical external efficacy level, ej*E, a critical time interval for interaction Δj*t and an actual time interval of contact Δjt. Hence, a better way to define avidity is
𝒜=e1EΔ1t+e2EΔ2t+e3EΔ3t+...+eNEΔNt |
The developing T cell-pMHC interface we call the immunological synapse sends out signals to its surroundings that enhance these interactions. Eventually, if the foreign peptide initiated first interaction exceeds the threshold affinity level needed for a sustained reaction, the resulting avidity increase causes a clonal expansion of the T cells that bind with the foreign antigen. It may also cause a clonal expansion of the pMHCs that are sufficiently similar to the foreign peptide. Hence, in addition to a greatly increased concentration of the T cells needed to remove the cells exhibiting the foreign peptide in an pMHC complex, there may also be a vast increase in the T cells bearing pMHCs that are similar to the pMHC with the foreign peptide. Now the foreign peptide bearing cells are lysed and within a short time removed. So how does this kind of event trigger autoimmunity in some circumstances?
If something did happen, there would be a catalyzing event of the form
𝒜s=e1EΔ1t+e2EΔ2t+e3EΔ3t+...+eNEΔNt |
where 𝒜s is large enough to begin an interaction chain culminating in the construction of the full immunosynapse and a corresponding lysis event. If we let e*E be the small avidity usually seen with a TcR-pMHC interaction of a T cell with a self protein, to start this chain we want a foreign peptide catalyst event whose avidity exceeds e*E. Hence, we need
e1E=(ξ1−1)(1+K1act)[(sc)1]+K1P[P1]ξ1[(sc)1][α1γ1δ1−ξ−11α1γ1δ1−α1ξ−11α1γ1δ1−γ1]·[1K1actK1Sc[(sc)1]]>e*E=(ξ*−1)(1+K*act)[(sc)*]+K*P[P*]ξ*[(sc)*][α*γ*δ*−ξ−1*α*γ*δ*−α*ξ−1*α*γ*δ*−γ*]·[1K*actK*Sc[(sc)*]] | (4.1) |
The response e1E may be less than the value required for an immunosynaptic event. If this case e1E<e*E. The directional derivative is largest if we move from the current point in ligand space to the new one determined by moving in the direction of Du where u=∇eE/‖∇eE‖. In general, any unit vector u which leads to a positive directional derivative is acceptable. This implies the directional derivative value at the base point must satisfy
〈∇*eE,[ΔαΔδΔγ]〉=〈[(sc)*]KP*[P*]ξ*(1+K*act)+K*Sc[sc*][γ*δ*α*γ*α*δ*],[ΔαΔδΔγ]〉−〈[sc*]K*P[P*][K*act0ξ*K*Sc[sc*]],[ΔαΔδΔγ]〉>0. | (4.2) |
Looking at
(ξ1−1)(1+K1act)[(sc)1]+K1P[P1]ξ1[(sc)1]<U1−V1ξ1,W(sc)1>>(ξ*−1)(1+K*act)[(sc)*]+K*P[P*]ξ*[(sc)*]<U*−V*ξ*,W(sc)*> | (4.3) |
or
(ξ1−1)(1+K1act)[(sc)1]+K1P[P1]ξ1[(sc)1]‖U1−V1ξ1‖‖W(sc)1‖cos(φ)>(ξ*−1)(1+K*act)[(sc)*]+K*P[P*]ξ*[(sc)*]‖U*−V*ξ*‖‖W(sc)*‖cos(φ*) | (4.4) |
Hence, if the parameter changes resulted in a decrease in φ from φ*, the avidity would increase. Evidently, under normal circumstances this does not occur, but perhaps if there are underlying genetic abnormalities in T cell generation and so on, this avidity increase is possible and there is an autoimmune effect.
Since this is so busy, let's try a thought experiment. First, let's assume the only parameters that matter are δ and ξ. Also, redefine all the units here so that all the other parameters have value 1. Then, we can write
e1E=(ξ1−1)(2)+ξ1[δ1−ξ−11δ1−ξ−11δ1−10]·[1110]>e*E=(ξ*−1)(2)+ξ*[δ*−ξ−1*δ*−ξ−1*δ*−10]·[1110] |
or
ξ1−3+3ξ1δ1>ξ*−3+3ξ*δ* |
Thus, to get over the threshold, we need
ξ1(1+3δ1)>ξ*(1+3δ*). | (4.5) |
Clearly, this occurs if ξ1 > ξ* and δ1 > δ* which is the easiest possibility. Now ξ1 is the multiplier for the concentration of the scaffolding proteins we have prior to the binding of the ligand-i.e. the pMHC. The term δ1 is related to the coupling of the ligand pMHC and the scaffolding. If δ1 > 1, this means the coupling of the pMHC enhances the scaffolding protein complexes we lump together as Sc. So these two terms are related. Since we assume ξ1 > ξ*, we know this means the concentration of the scaffolding complex is upregulated due to the pMHC binding over the threshold value for the self pMHC. Further, since δ1 > δ*, we know the binding of the this self pMHC enhances the coupling between the scaffolding a bit over the threshold amount for this self pMHC. How could this happen? Here is a possible scenario. There is an initial viral infection creates an inflammatory response which upregulates the scaffolding protein complexes around a self pMHC. This pushes both ξ1 and δ1 up and allows the the serial method of T cell activation to create a mature IS. Although simplified enormously, this thought experiment allows us to focus on what we think is really important here. What is the mechanism or process that allows the self pMHC and TcR interaction to exceed the putative threshold value? Once this happens, the formation of the IS is essentially assured and self damage ensues. This clearly does not happen except in people who are susceptible.
We can do better if we do a numerical analysis of this model. If we assume all of the rate constants are set to one by either convenience or a suitable relabeling of units, we have the prototypical affinity equation
eE=(ξ−1)(1+Kact)[(sc)]+KP[P]ξ[(sc)][αγδ−ξ−1αγδ−αξ−1αγδ−γ]·[1KactKSc[(sc)]]=(ξ−1)(1+Kact)[(sc)]+KP[P][(sc)](αγδξ(1+Kact+KSc[(sc)]))−1−α−γKSc[(sc)]ξ |
Let 1 + Kact = 2 and KP[p][sc] = 1, KSc[(sc)] = 1 and so 1 + Kact + KSc[(sc)] = 3. Then we have
eE=2(ξ−1)+3αγδξ−1−α−γξ |
Finally, let A = αδγξ. Then we have
eE=2(ξ−1)+3A−1−α−γξ |
We can perform some simulation studies with this equation and use it to describe a model of autoimmunity. Consider the function FindAffinities as shown in Appendix A.1. Our simulations calculate all the affinities associated with a fixed value of A = αδξγ. These parameters control scaffolding, TcR - pMHC binding, cytokine mediated positive feedback to strengthen the binding and so forth. For such a constant A, the affinities could be negative corresponding to no interaction between the TcR and the pMHC at all. The code is self explanatory. First, we compute all the affinities to find the minimum and maximum values, Am and AM, for this choice of A. Then scale the affinities into [0,1] using the transformation y ∈ [Am, AM] → (y - Am)/(AM - Am). Then for the choice Nbin divide [0,1] into Nbin-1 uniform pieces and count how many affinities land in each bin. From this, we can then calculate frequencies and the original affinity values associated with each bin.
It is easy to see that as A increases AM increases also. So a small A is equivalent to a low affinity for binding. In Figure 4a we see the affinity vs. frequency curve for A = .05. This is generated with
Note, the affinities range from −2.5 to −1.3 so for A = 0.05 so there is no binding event. If we increase A to A = .15, we see a slightly different graph in Figure 4b. Note the affinity spread has increased from [−2.5, −1.3] to [−2.3, −1.1] which shows the affinities are increasing although well below what is needed for binding event. If we increase to A = .75, we see the affinity range has now included some positive possibilities: the affinities belong to [−1.3, 0.07] so there is a small probability of a positive affinity with a very small probability of a binding event. We show this in Figure 5a.
Increasing to A = 1.5, the affinity range increase to [−0.04, 166] which means the probability of a binding event is greatly enhanced. We show this in Figure 5b. In a real situation, there will be a pool of T cells whose binding possibilities are modeled by a range of the parameters α, δ, ξ and γ determined from a distribution of A choices. The DrawSummedAffinities shown in Appendix A.2 allows us to draw affinity results from multiple families of P distributions. We can use this code to sum the affinities for A values in various ranges. Here we use A in [0.15, 0.35], [1.15, 1.35], [3.15, 5.35] and [7.15, 10.35].
This code generates the plot we see in Figure 6a. In Figure 6a, we see how the affinity families has affinities that increase from families with no possibility of a binding event to families with affinities as high as 15. If the formation on an immunosynapse requires affinities of say 60 or so, this can be achieved by positive feedback from earlier binding events to slowly increase the affinity. To see this graphically, consider another plot of affinity families as shown in Figure 6b. The probability of immunosynapse formation is quite high now.
We will denote the αδγξ Ligand - Receptor interaction model by LRαδγξ. The affinity/avidity obtained in this binding depends on many factors: this parameter space is actually nine dimensional and is comprised of realistic biological values of α, δ, γ, ξ, Kact, Kp, KSc, [P] and [sc]. The corresponding point cloud of data values lives in a topological space which whose structure is determined by the data. We suspect that once we move in parameter space from the maximum affinity value for a given A, that the remaining parameter values show a decrease in affinity which will look like part of a standard normal distribution. We can show this by looking at some computational results. The code in the function MaxAffinityCurve shown in Appendix A.3 computes the affinities for a select collection of parameter values for a given A and then returns two vectors: one is the set of all γξ values and other is the affinities. To find the maximum affinity curves for A = 1.0, 5.0, 10.0 and A = 30, we use the code below. In Figure 7, we show the results. Note as A increases, we see the right half of a typical normal distribution curve for the maximum affinity. As γxi moves away for the maximum, the affinities decrease. The LRαδγξ makes it clear that the maximum affinity value depends on much more that the epitope T cell binding achieved in the pMHC groove.
From Figure 7, we see that the probability of an effective immunosynapse creation increases as A increases. The affinities for the lower A values, by themselves, as not enough to initiate such an event. We can also see that for an autoimmune event to occur there must be a chain of T Cell and pMHC interactions that push diagonally upward in the γξ - affinity positive quadrant.
There are a number of very interesting papers that are relevant to a study of viral infections in the CNS and their possible involvement in autoimmune disease. For example, Multiple Sclerosis (MS) ( see the review [7]). may be due to a reaction against a virus and some believe the virus takes part in the immunopathology of the disease. There are other diseases which are similar such as Devic's disease, Marburg's variant, tumefactive MS, Balos's concentric sclerosis and acute disseminated encephalomyletis (ADEM) which all have interesting differences. From [7],
The most common hypothesis describes MS as an autoimmune disease induced by autoreactive T and B cells, which recognize myelin autoantigens and trigger inflammation in the CNS. The latter leads to the loss of myelin sheaths and CNS nervous conductivity and subsequently to the death of neurons.
It is known that autoreactive lymphocytes are present in all healthy people but it is not known why these autoreactive cells are triggered in MS patients. Again, from [7],
One possibility is that specific major histocompatibility complex (MHC) haplotype/autogen combinations together with specific genetic traits mainly of the immune system and additional triggering by infectious diseases initiates disease in MS patients.
There are a number of possibilities, We know there is a group of patients for which the primary problem is some sort of oligodendrocyte (ODC) degeneration which then triggers a subsequent inflammation that targets the mylenin sheathing. Four patterns of MS lesions were described.
1. Pattern I is primarily T cell based inflammatory type with few B cells and their complement.
2. Pattern II has a major antibody component so in addition to T cells there are also B cells, antibodies and complement.
3. Pattern III is dominated by apoptotic ODCs.
4. Pattern IV is characterized by a non apoptotic death of the ODCs.
Pattern I and II are associated with all forms of MS, Pattern III dominated in acute early MS. There is evidence that the ODC degenerative process could precede inflammation as well from a study of patients who died soon after a relapse. In a majority of them, Type III lesions were found but no leukocytes and mylenin degradation. So it is possible that in some cases of MS, the autoreactive lymphocyte scenario is secondary following an ODC problem. Of course, which of these scenarios is right determines the course of treatment. If the disease is primarily autoreactive then we would do one thing, but if the disease is primarily ODC related, we would need another approach. Hence, the animal model chosen to test the efficacy of new drugs depends on which of these scenarios is being followed. The primary animal model now being used is the Experimental autoimmune encephalomyletis (EAE) which is induced in mice. We know people who have been immunized against mylenin basic protein for rabies using material purified from rabbit brains have developed ADEM which resembles EAE. So the EAE mice model makes a lot of sense to pursue. EAE is induced in mice a variety of ways and generates Type I lesions and sometimes Type II. The methods at first involved injections of various mylenin proteins such as mylenin oligodendrocyte glycoprotein (MOG). However, new techniques have arisen using transgenesis. These used TcR transgenic models established in various mouse lines. These models showed spontaneous EAE and allowed new discoveries about tolerogenic mechanisms due to peripheral tolerance to T regulatory cells. A breakthrough occurred when the variant MOG35-55 peptide was discovered to by autoantigenic and an efficient EAE inducing peptide.
An unresolved question in MS is that the number of CD8 T cells might outnumber the number of CD4 T cells in the MS lesions. There are questions about whether the mylenin damage is primarily due to CD4 T cells or whether it can be due to CD8 T cells. A number of experiments have been performed which show CD8 T cells can cause mylenin damage without the CD4 T cells being involved, but there remain questions and as of the year of the review in [7], it has not been resolved. Nevertheless, we are taking the view that the CD8 T cell involvement is possible along with the CD4 T cell pathway. Experiments in the 1980's showed conclusively that T cells are critical for the induction of EAE. However, mysteries remained. At first, EAE was thought of as a Th1 autoimmune disease where Th1 was an influential line of T helper cells. The amount of IFN-γ released from periphery of these cells and the lesions seemed very important. But then it was found that mice lacking IFN-γ were not resistant to EAE and instead were hyper susceptible. Further, when mice lacked the materials to develop Th1, such as IL-12,they were also hyper susceptible to EAE. Eventually it was shown that IL-23 is essential for the development of EAE. It was then shown that IL-23 is needed to develop T helper cells that produce IL-17. These cells are now named TH17s and have been accepted as a new T helper population. Further work has shown that IL-23 may not be needed for the development of the TH17 sub population but actually helps stabilize them.
Recently, there have been a number of illuminating reviews. There is a good summary discussion of how immunomodulation might effect MS in [8]. The genetics of multiple sclerosis is explored in [9]. The creation of the correct pools of B and T Cells from naive ones is determined by the individual's genetics to some extent. Hence, it is possible there is a genetic component to the self affinity model equations. Further information about the interplay between bacteria and viral vectors for possible MS causal analysis can be found in examining how therapies to find infections helps ameliorate MS progression [10]. Further, the role of viruses in the onset of MS is not well understood and in [11], there is a careful review of what we think we know. There are also new tools on the horizon: engineered T Cells can be used to both probe immune response and potentially help with disease. A review of the use of synthetic biology tools used to expand T Cell capability is in [12]. Some of the engineered T Cells are called CAR-T Cells and their design and implementation are outlined in [13]. The models we have presented may be useful in delineating how to find the right engineered T Cell design for therapeutic use. Finally, a most interesting connection between the intestinal biome and autoimmunity in the brain is presented in [14]. The model discussed shows how such gut interactions might enable immune interactions to cross the blood brain barrier surrounding the brain. It is notable that a healthy gut biome probably includes hundreds if not thousands of coexisting bacterial species whose presence in the gut influences immune response. It is essential that we study abstract models of these systems as the number of independent components is beyond our ability to simulate in detail. The models we present here are a small step in this direction.
So the question is why does a T cell become pathogenic? A T cell type produces many different factors and cytokines. Further experiments have shown that the Th17 sub population is very plastic and these cells tend to lose their cytokines and transcription factors both in vitro and in vivo. They readily express IFN-γ under appropriate conditions.
Our model gives us some insight into a pathway to an autoimmune event in the case of a long term viral infection as described above. If a viral agent has an epitope similar to a self protein--for example, simple to a protein carried by a mylenin sheath cell--then a large pool of T Cells is maintained whose affinity for this self protein is close to a binding event. From the graphs of frequencies above, we can see if this pool is large enough, even a small probability of a positive binding can lead to a successful binding event.
If the initial binding event between the TcR and the pMHC holds for T seconds before the scaffolding proteins break bonds, there is a window of T seconds for another virally generated T Cell to bind with the pMHC containing the self protein peptide. The probability of another binding is much enhanced by having a large pool of T cells that bind weakly to the self protein pMHC. If this binding occurs before T seconds elapses, additional positive feedback occurs.
Each binding event which leads to positive feedback is effectively increasing A and from Figure 6b, we see the increase in A moves us to the right on the Figure. This increases the probability of immunosynapse formation and a subsequent lysis event. The avidity chain
𝒜=e1EΔ1t+e2EΔ2t+e3EΔ3t+...+eNEΔNt |
thus effectively is a chain for increasing A values leading to affinity spreads which are better at initiating immunosynapse formation. It is easy, however, for this chain to break. All it takes is one Δi to be too large and scaffolding proteins and other positive feedback events based on cytokine signalling will move backwards and the movement forward to immunosynapse formation will be compromised.
Effectively, a successful autoimmune event requires a smooth continuous path through α, δ, ξ and γ space to be successful. As a mathematical aside, this will not be possible in general in the end state desired lies in a separate component of the topological space which the α, δ, ξ and γ reside in.
Note our model of an affinity threshold computation depends on a wide variety of factors including geometry of the interaction between the T cell and the pMHC via the angle φ and ligand interactions of various sorts. The affinity threshold equation in conjunction with serial chains of T cell-pMHC bindings sheds insight into some autoimmune events by helping us think about the rise of a weakly bound T cell-pMHC affinity to the level of lysis. The affinity threshold equation is just the first step however, As discussed above, we also need a general model of cytokine signalling and a way to organize the immunosynapse responses in a network setting. To organize the immunosynapse responses, we also need to look at the cytokine signal data modeling in terms of manifolds and to apply ideas from topology.
This research has been supported by grant W911NF-17-1-0455 from the Army Research Office.
The author has no conflicts of interest in this paper.
[1] |
S. Timotheou, O. Miliou, Y. Dimitriadis, S. V. Sobrino, N. Giannoutsou, R. Cachia, et al., Impacts of digital technologies on education and factors influencing schools' digital capacity and transformation: A literature review, Educ. Inf. Technol., 28 (2023), 6695–6726. https://doi.org/10.1007/s10639-022-11431-8 doi: 10.1007/s10639-022-11431-8
![]() |
[2] |
P. Twining, J. Raffaghelli, P. Albion, D. Knezek, Moving education into the digital age: The contribution of teachers' professional development, J. Comput. Assisted Learn., 29 (2013), 426–437. https://doi.org/10.1111/jcal.12031 doi: 10.1111/jcal.12031
![]() |
[3] |
J. Engelbrecht, S. Llinares, M. C. Borba, Transformation of the mathematics classroom with the internet, ZDM-Math. Educ., 52 (2020), 825–841. https://doi.org/10.1007/s11858-020-01176-4 doi: 10.1007/s11858-020-01176-4
![]() |
[4] |
M. Beardsley, L. Albo, P. Aragon, D. Hernandez-Leo, Emergency education effects on teacher abilities and motivation to use digital technologies, Br. J. Educ. Technol., 52 (2021), 1455–1477. https://doi.org/10.1111/bjet.13101 doi: 10.1111/bjet.13101
![]() |
[5] |
S. J. Hegedus, L. Moreno-Armella, The emergence of mathematical structures, Educ. Stud. Math., 77 (2011), 369–388. https://doi.org/10.1007/s10649-010-9297-7 doi: 10.1007/s10649-010-9297-7
![]() |
[6] | C. Zhou, Children's cognitive abilities in the context of the digital transformation of basic education, Interact. Learn. Environ., 30 (2022), 1–13. |
[7] | MOE, Standards of Mathematics Curriculum for Compulsory Educaiton (2022 year version)(in Chinese), Beijing Normal University Publishing House, Beijing, China, 2022. |
[8] | NCTM, Principles and Standards for School Mathematics, NCTM, Reston, VA, 2000. |
[9] |
O. Birgin, F. Topuz, Effect of the GeoGebra software-supported collaborative learning environment on seventh grade student' geometry achievement, retention and attitudes, J. Educ. Res., 114 (2021), 474–494. https://doi.org/10.1080/00220671.2021.1983505 doi: 10.1080/00220671.2021.1983505
![]() |
[10] |
O. Birgin, K. U. Yazici, The effect of GeoGebra software-supported mathematics instruction on eighth-grade students' conceptual understanding and retention, J. Comput. Assisted Learn., 37 (2021), 925–939. https://doi.org/10.1111/jcal.12532 doi: 10.1111/jcal.12532
![]() |
[11] |
R. Bozic, D. Takaci, G. Stankov, Influence of dynamic software environment on students' achievement of learning functions with parameters, Interact. Learn. Environ., 29 (2019), 655–669. https://doi.org/10.1080/10494820.2019.1602842 doi: 10.1080/10494820.2019.1602842
![]() |
[12] |
E. Cekmez, Investigating the effect of computer-supported instruction on students' understanding of different representations of two-variable inequalities, Interact. Learn. Environ., 31 (2021). https://doi.org/10.1080/10494820.2021.1926288 doi: 10.1080/10494820.2021.1926288
![]() |
[13] |
M. Demir, Y. Zengin, The effect of a technology-enhanced collaborative learning environment on secondary school students' mathematical reasoning: A mixed method design, Educ. Inf. Technol., 28 (2023), 9855–9883. https://doi.org/10.1007/s10639-023-11587-x doi: 10.1007/s10639-023-11587-x
![]() |
[14] |
S. Radović, M. Radojičić, K. Veljković, M. Marić, Examining the effects of Geogebra applets on mathematics learning using interactive mathematics textbook, Interact. Learn. Environ., 28 (2020), 32–49. https://doi.org/10.1080/10494820.2018.1512001 doi: 10.1080/10494820.2018.1512001
![]() |
[15] |
M. S. Uwurukundo, J. F. Maniraho, M. T. Rwibasira, Effect of GeoGebra software on secondary school students' achievement in 3-D geometry, Educ. Inf. Technol., 27 (2022), 5749–5765. https://doi.org/10.1007/s10639-021-10852-1 doi: 10.1007/s10639-021-10852-1
![]() |
[16] |
H. Zulnaidi, E. Oktavika, R. Hidayat, Effect of use of GeoGebra on achievement of high school mathematics students, Educ. Inf. Technol., 25 (2020), 51–72. https://doi.org/10.1007/s10639-019-09899-y doi: 10.1007/s10639-019-09899-y
![]() |
[17] |
D. Dalby, M. Swan, Using digital technology to enhance formative assessment in mathematics classrooms, Br. J. Educ. Technol., 50 (2019), 832–845. https://doi.org/10.1111/bjet.12606 doi: 10.1111/bjet.12606
![]() |
[18] |
D. Juandi, Y. S. Kusumah, M. Tamur, K. S. Perbowo, T. T. Wijaya, A meta-analysis of GeoGebra software decade of assisted mathematics learning: What to learn and where to go?, Heliyon, 7 (2021), e06953. https://doi.org/10.1016/j.heliyon.2021.e06953 doi: 10.1016/j.heliyon.2021.e06953
![]() |
[19] |
J. S. Pyper, Learning about ourselves: A review of the mathematics teacher in the digital era, Can. J. Sci. Math. Technol. Educ., 17 (2017), 234–242. https://doi.org/10.1080/14926156.2017.1297509 doi: 10.1080/14926156.2017.1297509
![]() |
[20] |
A. Baccaglini-Frank, To tell a story, you need a protagonist: How dynamic interactive mediators can fulfill this role and foster explorative participation to mathematical discourse, Educ. Stud. Math., 106 (2021), 291–312. https://doi.org/10.1007/s10649-020-10009-w doi: 10.1007/s10649-020-10009-w
![]() |
[21] |
G. Bozkurt, K. Ruthven, Classroom-based professional expertise: A mathematics teacher's practice with technology, Educ. Stud. Math., 94 (2017), 309–328. https://doi.org/10.1007/s10649-016-9732-5 doi: 10.1007/s10649-016-9732-5
![]() |
[22] |
I. F. Rahmadi, Z. Lavicza, S. Arkün Kocadere, T. Houghton, M. Hohenwarter, The strengths and weaknesses of user-generated microgames for assisting learning, Educ. Inf. Technol., 27 (2022), 979–995. https://doi.org/10.1007/s10639-021-10635-8 doi: 10.1007/s10639-021-10635-8
![]() |
[23] |
F. Zhu, B. Xu, The role of dynamic geometry software in teacher-student interactions: Stories from three Chinese mathematics teachers, Sustainability, 15 (2023), 7660. https://doi.org/10.3390/su15097660 doi: 10.3390/su15097660
![]() |
[24] |
M. J. Bossé, E. S. Young, A. Bayaga, K. Lynch-Davis, A. DeMarte, C. Fountain, Cognitive processes in problem solving in a dynamic mathematics environment, Int. J. Math. Teach. Learn., 21 (2020), 174–196. https://doi.org/10.4256/ijmtl.v21i2.273 doi: 10.4256/ijmtl.v21i2.273
![]() |
[25] |
M. Cevikbas, G. Kaiser, A systematic review on task design in dynamic and interactive mathematics learning environments (DIMLEs), Mathematics, 9 (2021), 399. https://doi.org/10.3390/math9040399 doi: 10.3390/math9040399
![]() |
[26] |
D. Martinovic, Z. Karadag, Dynamic and interactive mathematics learning environments: The case of teaching the limit concept, Teach. Math. Its Appl., 31 (2012), 41–48. https://doi.org/10.1093/teamat/hrr029 doi: 10.1093/teamat/hrr029
![]() |
[27] |
Z. Yuan, J. Liu, X. Deng, T. Ding, T. T. Wijaya, Facilitating conditions as the biggest factor influencing elementary school teachers' usage behavior of dynamic mathematics software in China, Mathematics, 11 (2023), 1536. https://doi.org/10.3390/math11061536 doi: 10.3390/math11061536
![]() |
[28] |
S. S. Dogruer, D. Akyuz, Mathematical practices of eighth graders about 3D shapes in an argumentation, technology, and design-based classroom environment, Int. J. Sci. Math. Educ., 18 (2020), 1485–1505. https://doi.org/10.1007/s10763-019-10028-x doi: 10.1007/s10763-019-10028-x
![]() |
[29] |
D. Oner, Tracing the change in discourse in a collaborative dynamic geometry environment: From visual to more mathematical, Int. J. Comput Supported Collab. Learn., 11 (2016), 59–88. https://doi.org/10.1007/s11412-016-9227-5 doi: 10.1007/s11412-016-9227-5
![]() |
[30] |
A. H. Abdullah, N. S. Misrom, U. H. A. Kohar, M. H. Hamzah, Z. M. Ashari, D. F. Ali, et al., The effects of an inductive reasoning learning strategy assisted by the GeoGebra software on students' motivation for the functional graph Ⅱ topic, IEEE Access, 8 (2020), 143848–143861. https://doi.org/10.1109/ACCESS.2020.3014202 doi: 10.1109/ACCESS.2020.3014202
![]() |
[31] | P. Shadaan, K. E. Leong, Effectiveness of using GeoGebra on students' understanding in learning circles, Malays. Online J. Educ.Technol., 1 (2013), 1–11. |
[32] |
E. Cekmez, Using dynamic mathematics software to model a real-world phenomenon in the classroom, Interact. Learn. Environ., 28 (2020), 526–538. https://doi.rog/10.1080/10494820.2019.1674882 doi: 10.1080/10494820.2019.1674882
![]() |
[33] |
W. Daher, Middle school students' motivation in solving modelling activities with technology, EURASIA J. Math. Sci. Technol. Educ., 17 (2021), em1999. https://doi.org/10.29333/ejmste/11127 doi: 10.29333/ejmste/11127
![]() |
[34] | M. Hohenwarter, K. J. Fuchs, Combination of dynamic geometry, algebra and calculus in the software system GeoGebra, in Computer Algebra Systems and Dynamic Geometry Systems in Mathematics Teaching, Proceedings of Sprout-Slecting Conference (Sarvari, Cs. Hrsg.), Pecs, Hungary, (2005), 128–133. |
[35] |
E. Z. Hutkemri, Impact of using GeoGebra on students' conceptual and procedural knowledge of limit function, Mediterr. J. Soc. Sci., 5 (2014), 873–881. https://doi.org/10.5901/mjss.2014.v5n23p873 doi: 10.5901/mjss.2014.v5n23p873
![]() |
[36] |
E. Çekmez, B.Ö. Bülbül, An example of the use of dynamic mathematics software to create problem-solving environments that serve multiple purposes, Interact. Learn. Environ., 26 (2018), 654–663. https://doi.org/10.1080/10494820.2017.1385029 doi: 10.1080/10494820.2017.1385029
![]() |
[37] |
H. Jacinto, S. Carreira, Mathematical problem solving with technology: The techno-mathematical fluency of a student-with-GeoGebra, Int. J. Sci. Math. Educ., 15 (2017), 1115–1136. https://doi.org/10.1007/s10763-016-9728-8 doi: 10.1007/s10763-016-9728-8
![]() |
[38] |
A. Leung, A. M. S. Lee, Students' geometrical perception on a task-based dynamic geometry platform, Educ. Stud. Mathe., 82 (2013), 361–377. https://doi.org/10.1007/s10649-012-9433-7 doi: 10.1007/s10649-012-9433-7
![]() |
[39] |
J. Olsson, The contribution of reasoning to the utilization of feedback from software when solving mathematical problems, Int. J. Sci. Math. Educ., 16 (2018), 715–735. https://doi.org/10.1007/s10763-016-9795-x doi: 10.1007/s10763-016-9795-x
![]() |
[40] |
R. Ramirez-Ucles, J. F. Ruiz-Hidalgo, Reasoning, representing, and generalizing in geometric proof problems among 8th grade talented students, Mathematics, 10 (2022), 789. https://doi.org/10.3390/math10050789 doi: 10.3390/math10050789
![]() |
[41] | Q. Tan, Z. Yuan, A framework on mathematical problem solving with technology and its applications (in Chinese), J. Math. Educ., 30 (2021), 48–54 |
[42] |
K. K. Chan, S. W. Leung, Dynamic geometry software improves mathematical achievement: Systematic review and meta-analysis, J. Educ. Comput. Res., 51 (2014), 311–325. https://doi.org/10.2190/EC.51.3.c doi: 10.2190/EC.51.3.c
![]() |
[43] |
A. Sokolowski, Y. Li, V. Willson, The effects of using exploratory computerized environments in grades 1 to 8 mathematics: A meta-analysis of research, Int. J. STEM Educ., 2 (2015), 8. https://doi.org/10.1186/s40594-015-0022-z doi: 10.1186/s40594-015-0022-z
![]() |
[44] |
A. Dvir, M. Tabach, Learning extrema problems using a non-differential approach in a digital dynamic environment: The case of high-track yet low-achievers, ZDM-Math. Educ., 49 (2017), 785–798. https://doi.org/10.1007/s11858-017-0862-8 doi: 10.1007/s11858-017-0862-8
![]() |
[45] |
R. S. R. da Silva, L. M. Barbosa, M. C. Borba, A. L. A. Ferreira, The use of digital technology to estimate a value of Pi: Teachers' solutions on squaring the circle in a graduate course in Brazil, ZDM-Math. Educ., 53 (2021), 605–619. https://doi.org/10.1007/s11858-021-01246-1 doi: 10.1007/s11858-021-01246-1
![]() |
[46] |
M. Turgut, Reinventing geometric linear transformations in a dynamic geometry environment: Multimodal analysis of student reasoning, Int. J. Sci. Math. Educ., 20 (2022), 1203–1223. https://doi.org/10.1007/s10763-021-10185-y doi: 10.1007/s10763-021-10185-y
![]() |
[47] |
M. Turgut, Sense-making regarding matrix representation of geometric transformations in R2: A semiotic mediation perspective in a dynamic geometry environment, ZDM-Math. Educ., 51 (2019), 1199–1214. https://doi.org/10.1007/s11858-019-01032-0 doi: 10.1007/s11858-019-01032-0
![]() |
[48] |
Y. Zengin, Construction of proof of the fundamental theorem of calculus using dynamic mathematics software in the calculus classroom, Educ. Inf. Technol., 27 (2022), 2331–2366. https://doi.org/10.1007/s10639-021-10666-1 doi: 10.1007/s10639-021-10666-1
![]() |
[49] |
O. Birgin, H. Acar, The effect of computer-supported collaborative learning using GeoGebra software on 11th grade students' mathematics achievement in exponential and logarithmic functions, Int. J. Math. Educ. Sci. Technol., 53 (2022), 872–889. https://doi.org/10.1080/0020739X.2020.1788186 doi: 10.1080/0020739X.2020.1788186
![]() |
[50] |
M. Mthethwa, A. Bayaga, M. J. Bossé, D. Williams, GeoGebra for learning and teaching: A parallel investigation, South African J. Educ., 40 (2020), 1669. https://doi.org/10.15700/saje.v40n2a1669 doi: 10.15700/saje.v40n2a1669
![]() |
[51] |
K. Acikgul, R. Aslaner, Effects of Geogebra supported micro teaching applications and technological pedagogical content knowledge (TPACK) game practices on the TPACK levels of prospective teachers, Educ. Inf. Technol., 25 (2020), 2023–2047. https://doi.org/10.1007/s10639-019-10044-y doi: 10.1007/s10639-019-10044-y
![]() |
[52] |
O. Erkek, M. I. Bostan, Prospective middle school mathematics teachers' global argumentation structures, Int. J. Sci. Math. Educ., 17 (2019), 613–633. https://doi.org/10.1007/s10763-018-9884-0 doi: 10.1007/s10763-018-9884-0
![]() |
[53] |
S. B. Khoza, A. T. Biyela, Decolonising technological pedagogical content knowledge of first year mathematics students, Educ. Inf. Technol., 25 (2020), 2665–2679. https://doi.org/10.1007/s10639-019-10084-4 doi: 10.1007/s10639-019-10084-4
![]() |
[54] |
N. C. Verhoef, F. Coenders, J. M. Pieters, D. van Smaalen, D. O. Tall, Professional development through lesson study: Teaching the derivative using GeoGebra, Prof. Develop. Educ., 41 (2015), 109–126. https://doi.org/10.1080/19415257.2014.886285 doi: 10.1080/19415257.2014.886285
![]() |
[55] |
Y. Zengin, Development of mathematical connection skills in a dynamic learning environment, Educ. Inf. Technol., 24 (2019), 2175–2194. http://doi.org/10.1007/s10639-019-09870-x doi: 10.1007/s10639-019-09870-x
![]() |
[56] |
Y. Zengin, Effectiveness of a professional development course based on information and communication technologies on mathematics teachers' skills in designing technology-enhanced task, Educ. Inf. Technol., 28 (2023). https://doi.org/10.1007/s10639-023-11728-2 doi: 10.1007/s10639-023-11728-2
![]() |
[57] |
S. Gokce, P. Guner, Dynamics of GeoGebra ecosystem in mathematics education, Educ. Inf. Technol., 27 (2022), 5301–5323. https://doi.org/10.1007/s10639-021-10836-1 doi: 10.1007/s10639-021-10836-1
![]() |
[58] |
A. Yohannes, H. L. Chen, GeoGebra in mathematics education: A systematic review of journal articles published from 2010 to 2020, Interact. Learn. Environ., 29 (2021), 1–16. https://doi.org/10.1080/10494820.2021.2016861 doi: 10.1080/10494820.2021.2016861
![]() |
[59] |
M. D. Gurer, Examining technology acceptance of pre-service mathematics teachers in Turkey: A structural equation modeling approach, Educ. Inf. Technol., 26 (2021), 4709–4729. https://doi.org/10.1007/s10639-021-10493-4 doi: 10.1007/s10639-021-10493-4
![]() |
[60] |
K. F. Hollebrands, The role of a dynamic software program for geometry in the strategies high school mathematics students employ, J. Res. Math. Educ., 38 (2007), 164–192. https://doi.org/10.2307/30034955 doi: 10.2307/30034955
![]() |
[61] |
F. Ulusoy, İ. B. Turuş, The mathematical and technological nature of tasks containing the use of dynamic geometry software in middle and secondary school mathematics textbooks, Educ. Inf. Technol., 27 (2022), 11089–11113. https://doi.org/10.1007/s10639-022-11070-z doi: 10.1007/s10639-022-11070-z
![]() |
[62] |
K. Ruthven, R. Deaney, S. Hennessy, Using graphing software to teach about algebraic forms: A study of technology-supported practice in secondary-school mathematics, Educ. Stud. Math., 71 (2009), 279–297. https://doi.org/10.1007/s10649-008-9176-7 doi: 10.1007/s10649-008-9176-7
![]() |
[63] |
J. Hohenwarter, M. Hohenwarter, Z. Lavicza, Introducing dynamic mathematics software to second school teachers: The case of GeoGebra, J. Comput. Math. Sci. Teach., 28 (2008), 135–146. https://doi.org/10.13140/RG.2.2.15003.05921 doi: 10.13140/RG.2.2.15003.05921
![]() |
[64] |
R. Pierce, K. Stacey, Mapping pedagogical opportunities provided by mathematics analysis software, Int. J. Comput. Math. Learn., 15 (2010), 1–20. https://doi.org/10.1007/s10758-010-9158-6 doi: 10.1007/s10758-010-9158-6
![]() |
[65] |
R. Marrades, Á. Gutiérrez, Proofs produced by secondary school students learning geometry in a dynamic computer environment, Educ. Stud. Math., 44 (2000), 87–125. https://doi.org/10.1023/A:1012785106627 doi: 10.1023/A:1012785106627
![]() |
[66] |
G. G. Nagar, S. Hegedus, C. H. Orrill, High school mathematics teachers' discernment of invariant properties in a dynamic geometry environment, Educ. Stud. Math., 111 (2022), 127–145. https://doi.org/10.1007/s10649-022-10144-6 doi: 10.1007/s10649-022-10144-6
![]() |
[67] |
M. S. Alabdulaziz, S. M. Aldossary, S. A. Alyahya, H. M. Althubiti, The effectiveness of the GeoGebra programme in the development of academic achievement and survival of the learning impact of the mathematics among secondary stage students, Educ. Inf. Technol., 26 (2021), 2685–2713. https://doi.org/10.1007/s10639-020-10371-5 doi: 10.1007/s10639-020-10371-5
![]() |
[68] |
F. D. Davis, Perceived usefulness, perceived ease of use, and user acceptance of information technology, MIS Quarterly, 13 (1989), 319–340. https://doi.org/10.2307/249008 doi: 10.2307/249008
![]() |
[69] |
F. D. Davis, R. P. Bagozzi, P. R. Warshaw, User acceptance of computer technology: A comparison of two theoretical models, Manage. Sci., 35 (1989), 982–1003. https://doi.org/10.1287/mnsc.35.8.982 doi: 10.1287/mnsc.35.8.982
![]() |
[70] |
V. Venkatesh, F. D. Davis, A theoretical extension of the technology acceptance model: Four longitudinal field studies, Manage. Sci., 46 (2000), 186–204. https://doi.org/10.1287/mnsc.46.2.186.11926 doi: 10.1287/mnsc.46.2.186.11926
![]() |
[71] |
I. Ajzen, The theory of planned behavior, Organ. Behav. Hum. Dec. Processes, 50 (1991), 179–211. https://doi.org/10.1016/0749-5978(91)90020-T doi: 10.1016/0749-5978(91)90020-T
![]() |
[72] | E. M. Rogers, Diffusion of Innovations, The Free Press, New York, 1963. |
[73] |
V. Venkatesh, M.G. Morris, G. B. Davis, F. D. Davis, User acceptance of information technology: Toward a unified view, MIS Quarterly, 27 (2003), 425–478. https://doi.org/10.5555/2017197.2017202 doi: 10.5555/2017197.2017202
![]() |
[74] |
V. Venkatesh, J. Y. L. Thong, X. Xu, Consumer acceptance and use of information technology: Extending the unified theory of acceptance and use of technology, MIS Quarterly, 36 (2012), 157–178. https://doi.org/10.2307/41410412 doi: 10.2307/41410412
![]() |
[75] |
A. M. Al-Rahmi, A. Shamsuddin, E. Wahab, W. M. Al-Rahmi, U. Alturki, A. Aldraiweesh, et al., Integrating the role of UTAUT and TTF model to evaluate social media use for teaching and learning in higher education, Front. Public Health, 10 (2022), 905968. https://doi.org/10.3389/fpubh.2022.905968 doi: 10.3389/fpubh.2022.905968
![]() |
[76] |
L. Wan, S. Xie, A. Shu, Toward an understanding of university students' continued Intention to use MOOCs: When UTAUT model meets TTF model, SAGE Open, 10 (2020). https://doi.org/10.1177/2158244020941858 doi: 10.1177/2158244020941858
![]() |
[77] |
M. D. Williams, N. P. Rana, Y. K. Dwivedi, The unified theory of acceptance and use of technology (UTAUT): A literature review, J. Enterp. Inf. Manage., 28 (2015), 443–488. https://doi.org/10.1108/JEIM-09-2014-0088 doi: 10.1108/JEIM-09-2014-0088
![]() |
[78] |
C. M. Chao, Factors determining the behavioral intention to use mobile learning: An application and extension of the UTAUT model, Front. Psychol., 10 (2019), 1652. https://doi.org/10.3389/fpsyg.2019.01652 doi: 10.3389/fpsyg.2019.01652
![]() |
[79] |
A. Kundu, T. Bej, K. N. Dey, Investigating effects of self-efficacy and infrastructure on teachers' ICT use: An extension of UTAUT, Int. J. Web-Based Learn. Teach. Technol. (IJWLTT), 16 (2021), 21. https://doi.org/10.4018/IJWLTT.20211101.oa10 doi: 10.4018/IJWLTT.20211101.oa10
![]() |
[80] |
A. Clark-Wilson, O. Robutti, M. Thomas, Teaching with digital technology, ZDM-Math. Educ., 52 (2020), 1223–1242. https://doi.org/10.1007/s11858-020-01196-0 doi: 10.1007/s11858-020-01196-0
![]() |
[81] |
T. T. Wijaya, Y. Zhou, T. Houghton, R. Weinhandl, Z. Lavicza, F.D. Yusop, Factors affecting the use of digital mathematics textbooks in Indonesia, Mathematics, 10 (2022), 1808. https://doi.org/10.3390/math10111808 doi: 10.3390/math10111808
![]() |
[82] |
Ş. B. Tosuntaş, E. Karadağ, S. Orhan, The factors affecting acceptance and use of interactive whiteboard within the scope of FATIH project: A structural equation model based on the unified theory of acceptance and use of technology, Comput. Educ., 81 (2015), 169–178. https://doi.org/10.1016/j.compedu.2014.10.009 doi: 10.1016/j.compedu.2014.10.009
![]() |
[83] |
V. Venkatesh, H. Bala, Technology acceptance model 3 and a research agenda on interventions, Decis. Sci., 39 (2008), 273–315. https://doi.org/10.1111/j.1540-5915.2008.00192.x doi: 10.1111/j.1540-5915.2008.00192.x
![]() |
[84] |
M. A. O. Abbad, Using the UTAUT model to understand students' usage of e-learning systems in developing countries, Educ. Inf. Technol., 26 (2021), 7205–7224. https://doi.org/10.1007/s10639-021-10573-5 doi: 10.1007/s10639-021-10573-5
![]() |
[85] |
T. T. Wijaya, R. Weinhandl, Factors influencing students's continuous intentions for using micro-lectures in the post-COVID-19 period: A modification of the UTAUT-2 approach, Electronics, 11 (2022), 1924. https://doi.org/10.3390/electronics11131924 doi: 10.3390/electronics11131924
![]() |
[86] |
T. T. Wijaya, Y. Cao, R. Weinhandl, E. Yusron, Z. Lavicza, Applying the UTAUT model to understand factors affecting micro-lecture usage by mathematics teachers in China, Mathematics, 10 (2022), 1008. https://doi.org/10.3390/math10071008 doi: 10.3390/math10071008
![]() |
[87] |
T. Zhou, Y. Lu, B. Wang, Integrating TTF and UTAUT to explain mobile banking user adoption, Comput. Hum. Behav., 26 (2010), 760–767. https://doi.org/10.1016/j.chb.2010.01.013 doi: 10.1016/j.chb.2010.01.013
![]() |
[88] |
H. J. Lai, Investigating older adults' decisions to use mobile devices for learning, based on the unified theory of acceptance and use of technology, Int. Learn. Environ., 28 (2020), 890–901. https://doi.org/10.1080/10494820.2018.1546748 doi: 10.1080/10494820.2018.1546748
![]() |
[89] |
B. Šumak, A. Šorgo, The acceptance and use of interactive whiteboards among teachers: Differences in UTAUT determinants between pre- and post-adopters, Comput. Hum. Behav., 64 (2016), 602–620. https://doi.org/10.1016/j.chb.2016.07.037 doi: 10.1016/j.chb.2016.07.037
![]() |
[90] |
A. Gunasinghe, J. A. Hamid, A. Khatibi, S. M. F. Azam, The viability of UTAUT-3 in understanding the lecturer's acceptance and use of virtual learning environments, Int. J. Technol. Enhanced Learn., 12 (2020), 458–481. https://doi.org/10.1504/IJTEL.2020.110056 doi: 10.1504/IJTEL.2020.110056
![]() |
[91] |
J. Kim, K. S. S. Lee, Conceptual model to predict Filipino teachers' adoption of ICT-based instruction in class: Using the UTAUT model, Asia Pac. J. Educ., 42 (2022), 699–713. https://doi.org/10.1080/02188791.2020.1776213 doi: 10.1080/02188791.2020.1776213
![]() |
[92] |
A. Bandura, The explanatory and predictive scope of self-efficacy theory, J. Soc. Clin. Psychol., 4 (1986), 359–373. https://doi.org/10.1521/jscp.1986.4.3.359 doi: 10.1521/jscp.1986.4.3.359
![]() |
[93] |
J. Wu, H. Du, Toward a better understanding of behavioral intention and system usage constructs, Eur. J. Inf. Syst., 21 (2012), 680–698. https://doi.org/10.1057/ejis.2012.15 doi: 10.1057/ejis.2012.15
![]() |
[94] |
W. H. Loo, P. H. P. Yeow, S. C. Chong, User acceptance of Malaysian government multipurpose smartcard applications, Gov. Inf. Q., 26 (2009), 358–367. https://doi.org/10.1016/j.giq.2008.07.004 doi: 10.1016/j.giq.2008.07.004
![]() |
[95] |
S. Lakhal, H. Khechine, J. Mukamurera, Explaining persistence in online courses in higher education: A difference-in-differences analysis, Int. J. Educ. Technol. Higher Educ., 18 (2021), 19. https://doi.org/10.1186/s41239-021-00251-4 doi: 10.1186/s41239-021-00251-4
![]() |
[96] |
V. Marinković, A. Đorđević, Z. Kalinić, The moderating effects of gender on customer satisfaction and continuance intention in mobile commerce: A UTAUT-based perspective, Technol. Anal. Strategic Manage., 32 (2020), 306–318. https://doi.org/10.1080/09537325.2019.1655537 doi: 10.1080/09537325.2019.1655537
![]() |
[97] |
S. J. Jang, M. F. Tsai, Reasons for using or not using interactive whiteboards: Perspectives of Taiwanese elementary mathematics and science teachers, Australas. J. Educ. Technol., 28 (2012), 1451–1465. https://doi.org/10.14742/ajet.781 doi: 10.14742/ajet.781
![]() |
[98] |
S. Hu, K. Laxman, K. Lee, Exploring factors affecting academics' adoption of emerging mobile technologies: An extended UTAUT perspective, Educ. Inf. Technol., 25 (2020), 4615–4635. https://doi.org/10.1007/s10639-020-10171-x doi: 10.1007/s10639-020-10171-x
![]() |
[99] | Z. Yuan, S. Li, Developing prospective mathematics teachers' technological pedagogical content knowledge (TPACK): A case of normal distribution, in The 12th International Congress on Mathematical Education, (2012), 5804–5813. |
[100] |
C. Y. Lee, M. J. Chen, Developing a questionnaire on technology-integrated mathematics instruction: A case study of the AMA training course in Xinjiang and Taiwan, Br. J. Educ. Technol., 46 (2016), 1287–1303. https://doi.org/10.1111/bjet.12339 doi: 10.1111/bjet.12339
![]() |
[101] |
K. Sungur Gül, H. Ateş, An examination of the effect of technology-based STEM education training in the framework of technology acceptance model, Educ. Inf. Technol., 28 (2023), 8761–8787. https://doi.org/10.1007/s10639-022-11539-x doi: 10.1007/s10639-022-11539-x
![]() |
[102] | W. W. Chin, The partial least squares approach to structural equation modeling, in Mordern Methods for Business Research, Erlbaum, (1998), 295–358. |
[103] |
J. F. Hair, C. M. Ringle, M. Starstedt, PLS-SEM: Indeed a silver bullet, J. Mark. Theory Practice, 19 (2011), 139–151. https://doi.org/10.2753/MTP1069-6679190202 doi: 10.2753/MTP1069-6679190202
![]() |
[104] |
J. F. Hair, M. Sarstedt, C. M. Ringle, J. A. Mena, An assessment of the use of partial least squares structural equation modeling in marketing research, J. Acad. Mark Sci., 40 (2012), 414–433. https://doi.org/10.1007/s11747-011-0261-6 doi: 10.1007/s11747-011-0261-6
![]() |
[105] |
J. F. Hair, J. J. Risher, M. Sarstedt, C. M. Ringle, When to use and how to report the results of PLS-SEM, Eur. Business Rev., 31 (2019), 2–24. https://doi.org/10.1108/EBR-11-2018-0203 doi: 10.1108/EBR-11-2018-0203
![]() |
[106] | J. F. Hair, G. T. M. Hult, C.M. Ringle, M. Sarstedt, A Primer on Partial Least Squares Structural Equation Modeling (PLS-SEM), 3rd edition, Sage Publications, Thousand oaks, CA, 2022. |
[107] | J. F. Hair, M. Sarstedt, C. M. Ringle, S. P. Gudergan, Advanced Issues in Partial Least Squares Structural Equation Modeling, 2nd edition, SAGE, Thousand Oaks, CA, 2023. |
[108] |
H. M. Lin, M. H. Lee, J. C. Liang, H. Y. Chang, P. Huang, C. C. Tsai, A review of using partial least square structural equation modeling in e-learning research, Br. J. Educ. Technol., 51 (2020), 1354–1372. https://doi.org/10.1111/bjet.12890 doi: 10.1111/bjet.12890
![]() |
[109] | H. Wold, Partial least squares, in Encyclopedia of Statistical Sciences, (Eds. S. Kotz and N. L. Johnson), Wiley, (1985), 581–591. |
[110] | Z. Karadag, D. McDougall, GeoGebra as a cognitive tool: Where cognitive theories and technology meet, in Model-Centered Learning: Pathways to Mathematical Understanding Using GeoGebra, (Eds. L. Bu and R. Schoen), Sense Publishers, (2011), 169–181. |
[111] |
Z. A. Reis, S. Ozdemir, Using GeoGebra as an information technology tool: Parabola teaching, Proc. Soc. Behav. Sci., 9 (2010), 565–572. https://doi.org/10.1016/j.sbspro.2010.12.198 doi: 10.1016/j.sbspro.2010.12.198
![]() |
[112] |
R. Ziatdinov, J. R. Valles, Synthesis of modeling, visualization, and programming in GeoGebra as an effective approach for teaching and learning STEM topics, Mathematics, 10 (2022), 398. https://doi.org/10.3390/math10030398 doi: 10.3390/math10030398
![]() |
[113] |
M. Pittalis, Extending the technology acceptance model to evaluate teachers' intention to use dynamic geometry software in geometry teaching, Int. J. Math. Educ. Sci. Technol., 52 (2021), 1385–1404. https://doi.org/10.1080/0020739X.2020.1766139 doi: 10.1080/0020739X.2020.1766139
![]() |
[114] | Z. Yuan, M. Milner-Bolotin, A study of a TPACK-based subject-specific educational technology course and its implications: The case of "teaching mathematics and science through technology" course at the University of British Columbia (in Chinese), J. Math. Educ., 29 (2020), 23–28 |
[115] |
Z. Drezner, O. Turel, D. Zerom, A modified Kolmogorov–Smirnov test for normality, Commun. Stat. Simul. Comput., 39 (2010), 693–704. https://doi.org/10.1080/03610911003615816 doi: 10.1080/03610911003615816
![]() |
[116] | B. M. Byrne, Structural Equation Modeling With AMOS: Basic Concepts, Applications, and Programming, 3rd edition, Routledge, New York, 2016. |
[117] |
M. A. Moteri, M. Alojail, Factors influencing the supply chain management in e-health using UTAUT model, Electron. Res. Arch., 31 (2023), 2855–2877. https://doi.org/10.3934/era.2023144 doi: 10.3934/era.2023144
![]() |
[118] |
J. M. Santos-Jaén, A. León-Gómez, D. Ruiz-Palomo, F. García-Lopera, M. D. V. Martínez, Exploring information and communication technologies as driving forces in hotel SMEs performance: Influence of corporate social responsibility, Mathematics, 10 (2022), 3629. https://doi.org/10.3390/math10193629 doi: 10.3390/math10193629
![]() |
[119] |
M. D. Valls Martínez, P. A. Martín-Cervantes, A. M. Sánchez Pérez, M. D. Martínez Victoria, Learning mathematics of financial operations during the COVID-19 era: An assessment with partial least squares structural equation modeling, Mathematics, 9 (2021), 2120. https://doi.org/10.3390/math9172120 doi: 10.3390/math9172120
![]() |
[120] |
T. T. Wijaya, B. Yu, F. Xu, Z. Yuan, M. Mailizar, Analysis of factors affecting academic performance of mathematics education doctoral students: A structural equation modeling approach, Int. J. Environ. Res. Pub. Health, 20 (2023), 4518. https://doi.org/10.3390/ijerph20054518 doi: 10.3390/ijerph20054518
![]() |
[121] |
L. J. Cronbach, Coefficient alpha and the internal structure of tests, Psychometrika, 16 (1951), 297–334. https://doi.org/10.1007/BF02310555 doi: 10.1007/BF02310555
![]() |
[122] |
T. K. Dijkstra, J. Henseler, Consistent partial least squares path modeling, MIS Q., 39 (2015), 297–316. https://doi.org/10.25300/MISQ/2015/39.2.02 doi: 10.25300/MISQ/2015/39.2.02
![]() |
[123] |
C. E. Werts, R. L. Linn, K. G. Jöreskog, Intraclass reliability estimates: Testing structural assumptions, Educ. Psychol. Meas., 34 (1974), 25–33. https://doi.org/10.1177/001316447403400104 doi: 10.1177/001316447403400104
![]() |
[124] |
C. Fornell, D. F. Larcker, Evaluating structural equation models with unobservable variables and measurement error, J. Mark. Res., 18 (1981), 39–50. https://doi.org/10.1177/002224378101800104 doi: 10.1177/002224378101800104
![]() |
[125] |
J. Henseler, C. M. Ringle, M. Sarstedt, A new criterion for assessing discriminant validity in variance-based structural equation modeling, J. Acad. Mark. Sci., 43 (2015), 115–135. https://doi.org/10.1007/s11747-014-0403-8 doi: 10.1007/s11747-014-0403-8
![]() |
[126] |
J. Henseler, G. Hubona, P. A. Ray, Using PLS path modeling in new technology research: Updated guidelines, Ind. Manage. Data Syst., 116 (2016), 2–20. https://doi.org/10.1108/IMDS-09-2015-0382 doi: 10.1108/IMDS-09-2015-0382
![]() |
[127] | P. M. Bentler, EQS Structural Equations Program Manual, CA: Multivariate Software, Encino, CA, 1995. |
[128] |
P. M. Bentler, D. G. Bonett, Significance tests and goodness of fit in the analysis of covariance structures, Psychol. Bull., 88 (1980), 588–606. https://doi.org/10.1037/0033-2909.88.3.588 doi: 10.1037/0033-2909.88.3.588
![]() |
[129] |
F. Schuberth, M. E. Rademaker, J. Henseler, Assessing the overall fit of composite models estimated by partial least squares path modeling, Eur. J. Mark., 57 (2023), 1678–1702. https://doi.org/10.1108/EJM-08-2020-0586 doi: 10.1108/EJM-08-2020-0586
![]() |
[130] |
G. Shmueli, M. Sarstedt, J. F. Hair, J. H. Cheah, H. Ting, S. Vaithilingam, et al., Predictive model assessment in PLS-SEM: Guidelines for using PLSpredict, Eur. J. Mark., 53 (2019), 2322-2347. https://doi.org/10.1108/EJM-02-2019-0189 doi: 10.1108/EJM-02-2019-0189
![]() |
[131] | J. Henseler, C. M. Ringle, R. R. Sinkovics, The use of partial least squares path modeling in international marketing, in New Challenges to International Marketing, (Eds. R. R. Sinkovics and P. N. Ghauri), Emerald Group Publishing Limited, Bingley, (2009), 277–319. |
[132] | W. W. Chin, J. Dibbern, A permutation-based procedure for multi-group PLS analysis: Results of tests of differences on simulated data and a cross cultural analysis of the sourcing of information system services between Germany and the USA, (Eds. V. Esposito Vinzi, et al.), Springer, Berlin, (2010). |
[133] | M. Sarstedt, J. Henseler, C. M. Ringle, Multigroup analysis in partial least squares (PLS) path modeling: Alternative methods and empirical results, in Measurement and Research Methods in International Marketing, (Eds. M. Sarstedt, M. Schwaiger and C. R. Taylor), Emerald Group Publishing Limited, (2011), 195–218. |
[134] |
J. M. Bland, D. G. Altman, Multiple significance tests: The Bonferroni method, BMJ, 310 (1995), 170. https://doi.org/10.1136/bmj.310.6973.170 doi: 10.1136/bmj.310.6973.170
![]() |
[135] | J. F. Hair, W. Black, B. Babin, R. Anderson, Multivariate Data Analysis, 8th edition, Cengage Learning EMEA, Hampshire, United Kingdom, 2019. |