8The light-dependent reactions of photosynthesis take place in the plant chloroplast thylakoid 9 membrane, a complex three-dimensional structure divided into the stacked grana and 10 unstacked stromal lamellae domains. Plants regulate the macro-organization of 11 photosynthetic complexes within the thylakoid membrane to adapt to changing environmental 12 conditions and avoid oxidative stress. One such mechanism is the state transition which 13 regulates photosynthetic light harvesting and electron transfer. State transitions are driven by 14 changes in the phosphorylation of light harvesting antenna complex II (LHCII), which cause 15 a decrease in grana diameter and stacking, a decreased energetic connectivity between 16 photosystem II (PSII) reaction centres and an increase in the relative LHCII antenna size of 17 photosystem I (PSI) compared to PSII. Phosphorylation is believed to drive these changes by 18 weakening the intra-membrane lateral PSII-LHCII and LHCII-LHCII interactions and the 19 inter-membrane stacking interactions between these complexes, while simultaneously 20 increasing the affinity of LHCII for PSI. We investigated the relative roles and contributions 21 of these three types of interaction to state transitions using a lattice-based model of the 22 thylakoid membrane based on existing structural data, developing a novel algorithm to 23 simulate protein complex dynamics. Monte Carlo simulations revealed that state transitions 24 are unlikely to lead to a large-scale migration of LHCII from the grana to the stromal 25 lamellae. Instead, the increased light harvesting capacity of PSI is largely due to the more 26 efficient recruitment of LHCII already residing in the stromal lamellae into PSI-LHCII 27 supercomplexes upon its phosphorylation. Likewise, the increased light harvesting capacity 28 of PSII upon dephosphorylation was found to be driven by a more efficient recruitment of 29 LHCII already residing in the grana into functional PSII-LHCII clusters, primarily driven by 30 lateral interactions. 31 32 2 33 Statement of significance 34 For photosynthesis to operate at maximum efficiency the activity of the light-driven 35 chlorophyll-protein complexes, photosystems I and II (PSI and PSII) must be fine-tuned to 36 environmental conditions. Plants achieve this balance through a regulatory mechanism 37 known as the state transition, which modulates the relative light-harvesting antenna size and 38 therefore excitation rate of each photosystem. State transitions are driven by changes in the 39 extent of the phosphorylation of light harvesting complex II (LHCII), which modulate the 40 interactions between PSI, PSII and LHCII. Here we developed a novel algorithm to simulate 41 protein complex dynamics and then ran Monte Carlo simulations to understand how these 42 interactions cooperate to affect the organization of the photosynthetic membrane and bring 43 about state transitions. 44 45 Here, the integral membrane protein complexes; light-harvesting complex II (LHCII), 48 photosystem II (PS...