A continuum mechanics model of the plant cell wall reveals interplay between enzyme action and cell wall structure

Abstract Plant cell growth is regulated through manipulation of the cell wall network, which consists of oriented cellulose microfibrils embedded within a ground matrix incorporating pectin and hemicellulose components. There remain many unknowns as to how this manipulation occurs. Experiments have shown that cellulose reorients in cell walls as the cell expands, while recent data suggest that growth is controlled by distinct collections of hemicellulose called biomechanical hotspots, which join the cellulose molecule together. The enzymes expansin and Cel12A have both been shown to induce growth of the cell wall; however, while Cel12A’s wall-loosening action leads to a reduction in the cell wall strength, expansin’s has been shown to increase the strength of the cell wall. In contrast, members of the XTH enzyme family hydrolyse hemicellulose but do not appear to cause wall creep. This experimentally observed behaviour still awaits a full explanation. We derive and analyse a mathematical model for the effective mechanical properties of the evolving cell wall network, incorporating cellulose microfibrils, which reorient with cell growth and are linked via biomechanical hotspots made up of regions of crosslinking hemicellulose. Assuming a visco-elastic response for the cell wall and using a continuum approach, we calculate the total stress resultant of the cell wall for a given overall growth rate. By changing appropriate parameters affecting breakage rate and viscous properties, we provide evidence for the biomechanical hotspot hypothesis and develop mechanistic understanding of the growth-inducing enzymes. Graphic Abstract


Introduction
Faced with climate change and population growth, humanity needs plants that can cope with extreme weather events, diseases and rising demands on yield.As a result, understanding plant growth is essential to secure a sustainable future [1,2].Modification and adaptations to the plant genomes can provide a solution to optimising plant development.However, there remain many intriguing mysteries behind the mechanisms underlying plant growth.
As plants grow, their cell walls need to be strong enough to resist bursting but weak enough to allow permanent expansion.Controlled cell wall growth is an amazing feat, where some cells can increase in volume by over 30,000 times as they move from the meristem into maturation [3].The driving force behind cell growth is the isotropic turgor pressure maintained by the uptake of water into the vacuole accompanied by an a e-mail: euan.smithers@slcu.cam.ac.uk (corresponding author) b e-mail: j.luo.5@bham.ac.uk c e-mail: r.j.dyson@bham.ac.uk increase in cytoplasmic volume.This pressure inflates the cell membrane whose expansion is then restricted by the cell wall (a balloon in a box) creating a great tension in the cell wall, equivalent to 100-1000 atmospheres of tensile stress [3].A stressed cell wall will deform elastically or plastically depending on the cell wall mechanical properties.Plastic deformation or growth begins when the mechanical load (turgor pressure) exceeds some critical value (yield threshold).This extension leads to thinning of the cell which can be balanced by the deposition of new wall material to maintain cell wall integrity.Since turgor acts in all directions, cell walls must be highly anisotropic to allow directional growth.One of the earliest models of plant cell growth in one direction was the Lockhart equation, which states that the relative elongation rate, α * , is proportional to the turgor pressure, P * , if it is above a yield threshold, Y * [4].This idea is expressed in the form: where L * is the length of the cell and φ * is the extensibility.If the pressure does not exceed the yield thresh-old, the cell does not grow as the turgor cannot overcome the strength of the wall.In this article, we use the word "strength" to imply the cell wall's ability to resist deformation.This formulation was a useful initial model, but our understanding of plant growth has now improved.Rather than changes in turgor pressure, growth is often controlled via careful manipulation of the cell wall, which results in changes to the yield threshold or the extensibility [5].General reviews on mechanical modelling of plant growth can be found in [6][7][8], while further information on the biological details can be found in [5,[9][10][11].
In this paper, we create a mathematical model to investigate the interplay between cell wall structure and enzyme action in order to understand experimentally observed behaviour.Enzyme action is a crucial process in cell wall growth; however, it is generally underresearched in the biological context with even less modelling undertaken [6,12].
The cell wall consists of three main components, cellulose (CMF), hemicellulose (e.g arabinoxylan or xyloglucan) and pectin [5,9,[13][14][15].The cell wall properties are carefully mediated via active control of the wall's mechanical structure (e.g. by enzymatic action or new material deposition), altering either the yield or the post-yield behaviour and ultimately affecting the direction and rate of growth.
On the microscale, bond breakage and polymer network rearrangement (wall loosening) result in the relaxation of wall stress, allowing for the viscous flow of the cell wall.Cell wall loosening can be mediated by the action of proteins or enzymes, such as expansins, xyloglucan endotransglucosylase/hydrolase (XTH), and pectin-modifying enzymes, and are regulated by the action of hormones (morphogens), such as auxin, gibberellins, and abscisic acid [3,10].
Cellulose are long and stiff molecules embedded in a matrix of hemicellulose and pectin.Cellulose are deposited in the cell wall in lamella layers at a variety of angles [16], typically perpendicular to the growth direction, circling the cell, where they reinforce the cell against radial expansion [9].Cellulose molecules are also responsible for resisting the majority of the cell wall tension [17].Experimental observations of the cell wall have found that cellulose molecules are transversely directed after deposition in the inner cell wall and reorient to a longitudinal direction as they move to the outside of the cell wall during growth [18][19][20][21][22], however this might not always occur [23].
How the rings of cellulose are connected is not precisely known [12].It was previously thought that cellulose molecules are joined together via a tethered network where the cellulose fibres run parallel to each other and are continuously joined together by hemicellulose, which form hydrogen bonds with the cellulose and peel off when the network is deformed [10].There are several problems with this theory.Simulations have revealed that assuming a tethered network structure results in a much weaker cell wall than experimentally observed [24].It was also found that some plants that lack xyloglucan (mutant forms of Arabidopsis and cel-ery) displayed only a small amount of growth reduction [14] implying the role of xyloglucan in cell wall strength could have been exaggerated.According to a finite element model featuring a network of cellulose molecules tethered together by hemicellulose via hydrogen bonds, a deformed network is not strong enough to withstand the strain caused by turgor [25].These studies present evidence that the tethered network model is not a feasible explanation as to how the cell wall retains integrity.These results emphasise the mechanical role of other molecules like pectin and suggest that when xyloglucan is present, it could be concentrated around a limited number of distinct biomechanical hotspots (hotspots) where cellulose molecules come into close proximity with one another [14].
These concentrated hotspots could allow for controlled extension of the cell wall where these distinct spots are selected to allow slippage [9,10,14,26].Pectin is also theorised to have a role to play in these hotspots [16].A model testing this theory considered a network of cellulose connected by hotspots represented as linear springs [27].The model hypothesises that a group of short xyloglucan strands is stiffer than a single long strand, and when combined with pectin, the cell wall can produce the requisite wall stiffness to oppose turgor.The hotspot hypothesis claims that a small amount of degradation of the hotspots could lead to the load being carried by pectin, which then enables the viscous flow of the cell wall, providing a possible mechanism for growth.
Some enzymes affect the cell walls and possibly interact with the hotspots.Note that in this article, we refer to "wall loosening" as the action that directly causes stress relaxation, creep and hence growth, and "wall softening" as a decrease in the Young's modulus [10].XTH has been shown to have a hydrolysing action on the cell wall where it can cut and rejoin xyloglucans; despite this observed effect, it strangely does not induce significant cell wall extension [3,10].Cel12A, an enzyme present in fungi, has been shown to cause wall loosening.It has been suggested that Cel12A targets the hotspots by performing hydrolysis at these sites and leads to a reduction in wall strength [10,28] and thus causes both wall loosening and wall softening [12].Modelling efforts have offered an explanation of why these two enzymes hydrolysis action's have different effects; using coarse-grained molecular dynamics it was found that cellulose is the main load-bearing component, which could be the reason why enzymes purely targeting xyloglucan are ineffective [17].Expansin is a pH-controlled wall-loosening protein [3].Some experiments have shown expansin action to induce growth [12,29], while in contrast, other tests have observed the cell walls withstand more force without bursting [10].Unlike Cel12A, Expansin action seems to loosen but not soften the wall [12], but the mechanism remains unknown.There has been no observed enzymatic action by expansin [30], so this effect could be due to force dissipation by α-expansin.It is hypothesised that expansin targets the hotspots as there is evidence that they act on cellulose-cellulose sites [12,26] where they may induce slippage of the fibres [30,31].
This paper aims to produce a proof-of-concept model to test the plausibility of the hotspots hypothesis as an explanation for observed enzyme behaviour.We focus on inner tissues cells with predominately transversely orientated CMF with negligible cell division, for example cells within the root elongation zone [32][33][34].We model the primary cell wall as a continuum incorporating crosslink dynamics (between cellulose and hemicellulose) and calculate the stress resultant when acted on by a prescribed growth rate.Using the model, we examine hypothesised expansin, Cel12A and XTH enzyme action, and investigate the consequences of cellulose reorientation.This model is designed to be a simplification of the system in order to test the feasibility behind proposed wall structure and protein mechanisms.The focus is therefore on the cellulose crosslinks and not on the pectin ground matrix, whose direct contribution (which was discussed in [35]) shall be neglected in the model, while its possible cellulose crosslinking will be incorporated.We begin with an explanation of the mathematical model of the cell wall in Sect. 2. This is followed by an analysis of the model outputs and a discussion of the implications of these findings in Sect. 4. We summarise the results and draw conclusions in Sect. 5.

Model formulation
The model comprises three distinct aspects of cell wall dynamics: the emergent macroscopic stress and CMF orientation evolution from the microscopic cell wall network (Sect.2.2.1), hotspot bond density evolution (Sect.2.2.2) and enzyme action (Sect.2.2.3).We will first detail the assumptions behind each aspect (Sect.2.1), before deriving the relevant governing equations (Sect.2.2).This model is based on the framework originally developed in [35], with significant differences in the treatment of material properties, network composition and fibre orientation evolution.We will simplify the governing equations through nondimensionalisation (Sect.2.3) and solve the resulting system (Sect.3).The solutions are further simplified via asymptotic reduction (Sect.3.2) providing insights into the principal components controlling cell wall behaviour.

Model assumptions and set-up
The crucial output of the model is an expression for the axial stress resultant Σ * of the cell wall (see Fig. 1).This stress resultant captures the strength of the cell wall generated by the underlying polymer network, as well as the relationship between turgor pressure and the growth rate; higher values of turgor increase the load on the cell wall, resulting in raised stress levels and possible further extension of the wall.The stress resultant is calculated by summing up the stress held by each component of the cell wall polymer network Fig. 1 An idealised segment of the cell wall.The two curved (light blue) surfaces are the two boundaries of the cell wall.The thick (purple) rods represent the cellulose molecules reorienting as they approach the outside of the cell wall.The groups of thin (orange) fibres represent the hotspots which become increasingly stretched as they traverse to the outside of the cell.Equation (2) imposes this movement and deformation in the x-y plane.Note that the figure represents a simplification of biological reality as per the model, where biomechanical hotspots are distributed throughout the domain according to some dynamical density at all points in the cell wall.We assume that the cell wall is an evolving continuum with embedded cellulose molecules that are initially orientated perpendicularly to the axial direction.As the cell wall elongates, with growth rate α * , all molecules are stretched, and the cellulose molecules are additionally reoriented.These cellulose molecules are connected by a mix of cell wall components including hotspots crosslinks as described in Sect. 1 where the hotspot bond density and the hemicellulose number functions are n * hot and n * , respectively.( * denotes dimensional quantities throughout.)All func- tions that represent the whole hotspot are denoted with a subscript hot ( hot ) where functions that represent a single fibre are without.This bond density is dependent on the energy held in the fibres, meaning that as they get increasingly stretched via cell wall extension, they become more likely to break.As we are examining the principal growth direction, forces are resolved in this direction.Therefore, only axial hemicellulose are included since mechanically they matter the most.Enzyme action will be modelled by changes in the bond density evolution parameters.The wall segment is initially unstressed at time t * = 0 and undergoes uniform stretching in such a way that the cell wall grows with a fixed growth rate α * .We now explain how this growth rate is imposed.
We introduce the coordinate system (x * , y * , z * ) with the origin positioned on the edge of the outer wall (Fig. 1).The cell wall is stretched in the axial x * direction with the y * , z * directions being perpendicular to the growth direction, such that y * points radially inwards, and z * tangential to the cell wall.We define θ to be the angle between the cellulose fibres and the z * axis.New wall material is deposited on the inner surface of the wall at y * = h * , moves through the wall at rate α * until it reaches the outer wall surface at y * = 0. We therefore model the growth via the incompressible flow field, meaning that dx * /dt * = α * x * and dy * /dt * = −α * y * .This flow describes the material moving through the flat x-y plane in a thin segment of the cell wall (Fig. 1 and see figure 2 and 3 in [35] for further details).Due to the modelling domain being thin, we assume the cellulose molecules do not bend radially, and since there is no stretching in the z * direction, the model simplifies to a 2D system.We now consider how a generic element of cell wall material changes in length as the cell wall stretches.Let L * (y * , t * ) be the length of such a segment lying in the x * direction.The material extends with rate and so segment is therefore under an axial strain of s * = L * L * 0 − 1 where L * 0 is the natural length of the material segment.If the initial length of a segment positioned at y * = y * i at the non-constant deposition time t * = τ * is L * i , the segment length evolves according to L * = L * i e α * (t * −τ * ) on the characteristic y * = y * i e −α * (t * −τ * ) (from Eq. 3).The evolving segment length is then related to its position in the cell wall by L * = L * i y * i y * , encoding both its extension and trajectory.As a result of this formulation, the system of equations described in Sect.2.2 will all be partial differential equations dependent on both time and space.
To consider the amount of wall material contained in the wall segment and its deposition, we assume that at the inner surface of the cell wall, material is deposited such that constant cell wall thickness is maintained.This new wall material then gets pushed through the wall towards the outer surface via the flow (2).We find that wall density remains uniform for all time, assuming appropriate initial conditions (see appendix A for details).When new fibres are deposited, they are assumed to have zero stress and oriented in the z direction as they are yet to be subjected to tension or the material flow.This assumption will be reflected in the boundary conditions in Sect.2.2.

Stress resultant from the composite cell wall
The total stress resultant from the composite cell wall in the axial direction when growing in the axial direction depends on the mechanical properties of each constitutive (and evolving) part: the CMF which are bound together by both the hotspots and other cell wall components, the hemicellulose within the hotspots and the surrounding pectin ground matrix.We assume a standard linear solid-like system, so the total stress σ * (taken from the stress components in the axial direction) is the sum of the stresses arising from the biomechanical hotspots σ * hot and the cellulose molecules σ * C (Fig. 2).We let σ * H be the stress held in a hemicellulose fibre which we assume to be characterised by a Maxwell element where the hemicellulose strain, s H , is equal to the wall strain, s = s H .The Maxwell assumption means that the hemicellulose strain is the sum of the elastic and viscous contributions (stretching a spring and viscous damper/dashpot in series), s H = s H,e + s H,v , with the hemicellulose stress being equal to its elastic and viscous stress, with the boundary and initial conditions of σ * H (h * , τ * ) = σ * H (y * , 0) = 0 where τ * is the time the fibre is deposited, h * the wall thickness, ν * H the hemicellulose Young's modulus and μ * the material constant of viscosity.As each hotspot is a compact collection of hemicellulose fibres, we assume that the stress arising from the hotspots, σ * hot , is then the sum of the hemicellulose stresses contained within it, meaning σ * hot = n * σ * H .Here it is assumed that the wall strain is uniformly distributed across all the fibres and that the deformation is smooth, i.e. the strain in the fibres is the same as the wall strain.It has been proposed that the strain distribution through the wall may be discontinuous, being smaller at the hotspots [11].This effect could be incorporated through extension of the fibres, for example s H = δs or s H = s−δs C where s C is the cellulose strain and δ < 1 and is some parameter to scale the amount of strain imposed onto the hotspots.This would, however, incorporate yet another unknown parameter while C is represented by a box that resists the strain elastically proportionally to the cellulose's angle θ.
The whole unit's stress is then dependant on the biomechanical hotspot density n * hot such that when it is zero, the only contribution is coming from cellulose crosslinked by other molecules.On the right a representation of the cell wall continuum consisting of the single units being progressively stretched as they move through the cell wall causing the individual hemicelluloses bonds and hotspots to break (red crosses) introducing little impact on the results.(It would scale the breakage rate down and the stress resultant up due to cellulose's contribution.)Thus, discontinuity of the strain is not considered.Importantly, even though these hotspots are being strained the same amount, this strain is being released by the viscous slippage of the dashpots; so these fibres are under less elastic strain than the wall.And as the fibre breakage rate is solely dependent on this elastic strain (see Sect. 2.2.2), some aspects of the discontinuity are already included in the model.
Assuming the cellulose molecules are elastic, the stress term for the CMF in the axial direction is then, where ν * C is the Young's modulus, θ the fibre angle, sin θ a modifier that changes the cellulose stress as they become increasingly reoriented to the axial direction, and s C the strain of cellulose fibre, such that where L * C,0 is the cellulose resting length.The cellulose fibre angle θ and length L * C are orientated and deformed by the flow u * .Letting a be the cellulose direction vector, we shall assume the cellulose network is uniform in the x * and z * -direction as they are in the plane of the wall, and thus, a = (sin θ, 0, cos θ).Assuming the cell wall is a composite material with a preferred direction under an incompressible, transversely isotropic viscous flow, the evolution of the fibre director field, a, can be described by [36], where ζ * = a • e * • a represents the strain rate in the direction of the matrix, with e * = ∇ * u * + ∇ * u * T /2.We now have an expression that describes how cellulose fibres convect, stretch and reorient as they move through the cell wall via the flow, u * .Equation ( 7) can be re-expressed in terms of θ using the definition of a meaning ζ * = α * sin 2 θ, and upon substituting in Eq. (7) with boundary and initial conditions, respectively, being θ(h * , τ * ) = θ(y * , 0) = θ 0 (note θ 0 = 0, else the solution is simply θ = 0 for all t * > 0).From the derivation of Eq. ( 7), we also have an expression that describes the evolution of the cellulose length L * C when stretched via the flow as with boundary and initial conditions being Finally, the total stress resultant in the cell wall is then the sum of the stress of all components when crosslinked in each "layer" of the cell wall, i.e. integrating over the thickness of the wall, which gives where ρ * is the density of cellulose and (1 + A * n * hot ) is the modification of the cellulose stress due to crosslinking, noting that cellulose can only contribute to the wall stress if they are connected to one another, otherwise they are just pulled apart..The cellulose modification term in (12) has two contributions: the first represents non-hotspot crosslinks and the second hotspots crosslinks, with A * controlling the magnitude of the contribution.The integral is the contribution from the CMF (first term with σ * C ) and the xyloglucan hotspots (second term with σ * hot ), with the xyloglucan term dependent on the hotspot density for the same reasons as cellulose; the final term is the pectin matrix contribution which is assumed to provide an extensional viscosity due to its properties [37,38], with Γ * being the stiffness density of the matrix.This concludes the description of the stress resultant, so we proceed to characterise the crosslinking dynamics.

Bond density evolution equation
We now introduce equations that describe the hotspots density and hemicellulose number, n * hot and n * respectively.Assuming that the CMF do not break before the hotspots, it is then the hotspots connections that rupture to allow slippage of the cell wall components.We assume that no new bonds are formed inside the cell wall.To address this potential limitation, we could have included a stochastic bond reformation term.However, this would be equivalent to uniformly decreasing all bond breakage rates, so it would introduce more unknown parameters without producing new effects on the results in Sect. 4. The hotspots and the hemicellulose fibres within them are advected through the cell wall according to a Smoluchowski equation [39], with both n * hot and n * undergoing energy-dependant breakage: where 0 , with the breakage rates defined as In Eqs. ( 15) and ( 16), k * b T * is the thermal energy (k * b being the Boltzmann's constant and T * the absolute temperature), k * 0,hot and k * 0 are the breakage rates when the fibres are unstressed, U * hot and U * the deformation free energies and β and γ are parameters controlling how strongly the energies affect the breakage rates.This formulation entails a direct relationship between the stress held in the fibres and the breakage rates and has been used successfully in other models of adhesion dynamics [40][41][42][43].To ensure that the crosslinks stretch a significant length before breaking [35], we take β, γ 1; alternatively, taking larger values of these parameters allows us to model increased breakage rates caused by, for example, heightened hydrolysis action.As hemicellulose bonds break, the hotspots become easier to rupture, and as a result, we take β2 = ( We also assume that the hotspots are stronger than the hemicellulose crosslinks, so β < γ. To calculate the free energy potential in a single fibre, we consider the force F * on a fibre as it moves along the trajectory x * (t * ) = (x * (t * ), y * (t * ), z * (t * )).We set y * (t) = z * (t) = 0 as the molecules are only being extended in the axial direction, so with s * H,e being the elastic extension.The work done by F * is then Since the bond breakage is dependent on elastic strain, the force is also elastic, meaning F * x = −κ * x * where κ * is the stiffness of the springs and is equal to the area multiplied by ν * H divided by the length of the cellulose molecules.Since temperature is constant, the free energy equals the potential/work done and so For the hotspot density potential, the force applied is F * x = −κ * hot x * for some stiffness κ * hot , which should be affected by the number of intact fibres in the hotspot, n * ; we therefore take κ * hot = n * κ * , as the springs are in parallel so stiffness is additive and hence and This completes the description of crosslinking, and we proceed to detail the implementation of enzyme action.

Enzyme action
We consider a simple model of enzyme action.Recall that the enzymes Cel12A and XTH perform hydrolysis/cutting action, where Cel12A targets the hotspots junctions (cutting both cellulose and hemicellulose) and XTH targets the hemicellulose fibres.We simulate the Cel12A and XTH actions by increasing β and γ, respectively, in Eqs. ( 15) and (16).As Cel12A also digests the cellulose molecules, it might also decrease the cellulose density ρ.As the focus is on the hotspots, the modelling of this effect will be left to the appendix E. Due to the lack of consensus on how expansin works, we try two different simple methods of enzyme action.As mentioned in Sect. 1, expansin may work by allowing slippage in the fibres [31].The first method entails decreasing the viscosity of the dashpots, i.e. decreasing their resistance to the flow.This effect is modelled by the equation , where E * is the expansin action rate, μ * 0 the initial viscosity and μ * 1 the target viscosity with μ * 1 < μ * 0 .We choose this form to keep the model as simple as possible while ensuring that μ * decreases and that μ * = 0.
The second method imposes expansin action by increasing the resting length of the springs in the system allowing for stress relaxation.We assume that the resting length growth rate depends on the strain of the spring, rather than the extension; this assumption avoids the resting length exceeding the actual length.Thus, with the boundary and initial conditions, L 0 (h * , τ * ) = L * 0 (y * , 0) = l * 0 .With the system of equations fully described, we now proceed to simplify them through nondimensionalisation.

Model analysis
To evaluate the effect of the hotspots and enzyme action on the stress resultant and hence the overall mechanical behaviour, we proceed to analyse the model as described in Sect.2.2.The equations may be solved via a (lengthy) semi-analytic method; we merely state the solutions in Sect.3.1 and leave full details to the Appendices B, C and D. Since breakage rates are assumed to be small, we also employ asymptotic expansions to determine the leading-order behaviour of the cell wall, as shown in Sect.3.2; as before full details are consigned to F for brevity.

Semi-analytical solutions
Throughout the analysis, we recognise two contrasting regions of cell wall material: the thinning region 0 ≤ y ≤ e −αt contains material already present at t = 0, while the expanding e −αt < y ≤ 1 region contains newly-deposited material which has been added to the y = 1 surface at some t > 0. Recall from Sect.2.1 that due to the fluid flow, fibres deposited at time t = τ (which is non-constant) lie on a characteristic y = y i e −α t, where t = t − τ .For the initially present material, the deposition time is τ = 0 and initial position is 0 ≤ y i ≤ 1 meaning the extension is described by L = y i e αt .For the later-deposited material, we have τ > 0 and y i = 1, and thus, the extension can be given in terms of its position in the wall as L = 1 y .Beginning with the fibre stress terms, in the absence of enzymatic action, L 0 and μ are constant, and hence, Eq. ( 25) may be solved using an integrating factor to find the axial hemicellulose stress, In contrast, when enzyme action is incorporated, L 0 and μ are no longer constant, and hence, the solution (34) is invalid.To solve for variable spring rest length, we first divide Eq. ( 32) by dy dt = −αy to reformulate the differential equation for L 0 in terms of y.This equation does not have an analytical solution, and hence, we proceed to derive a (forward) finite difference expression for L 0 , which can be used for solving Eq. ( 25) with an integrating factor, yielding where L 0 is treated as a function of L.
For expansin action on the viscosity, μ, we begin by solving Eq. ( 31) using separation of variables to get . Equation ( 25) is then solved to give To determine the evolving bond densities and number, we solve (30) for n using the method of characteristics and combining Eqs. ( 24) and (C18) yields with G(y) Similarly from Eq. ( 29), with Solving Eq. ( 28) gives the cellulose angle from the horizontal, θ = arctan (tan θ 0 L) , (42) and thus, the cellulose extension length is, The stress resultant is hence constructed by separating the domain into two regions ≤ y ≤ e −αt and e −αt < y ≤ 1, to give where σ H is given by Eqs. ( 34), (35) or (37) for the different expansin mechanisms.The small parameter denotes a cut-off for the outer region of the cell wall; this ensures that the cellulose stress does not go to infinity and that the molecules will eventually break rather than becoming infinitely long.The steady-state stress resultant as t → ∞ is where the contributions from the fibres present at t = 0 have now disappeared.

Leading-order cell wall behaviour
An asymptotic expansion simplifies the model and extracts the leading-order component of the stress resultant integral when expansin is neglected (i.e.taking σ H from Eq. 34).The reader can skip this section or directly go to Eq. ( 57) for the final form.Since both γ, β 1 by definition, we can expand the steady-state solution (45).
We first expand the integral G (Eq. 38) considering the respective regions γ y ≤ 1 and ≤ y γ.By bounding the integral for G, integrating by parts and neglecting higher-order terms, we find where α = α α+ω for notational simplicity.From Eq. ( 38), the asymptotic approximation of n is then, The expansion of G hot proceeds similarly to give Combining, we find Notice that the simplification of n hot is not dependant on n.To complete the approximations for all y, we find switch-over values where the asymptotic behaviour of n and n hot changes from one regime to another (denoted χ and χ hot for the hemicellulose and hotspots bond densities respectively).These are the y-values close to γ or β hot such that the derivatives of the dominant exponents in n and n hot exceed some threshold Q, causing n and n hot to rapidly decrease.Thus, Approximations of the trigonometric term describing how much cellulose contribute to the axial stress and the cellulose extension are also required.To balance simplicity and accuracy, the following expansions are taken to such an order that they achieve a meansquared error of ≈ 10 −3 with respect to the numerical solution, and the domain is split into three regions corresponding to tan θ0 y , ≈, 1.In the first region, tan θ 0 + δ < y ≤ 1, where δ is a small number of our choosing to optimise the expansion's overall accuracy, we Taylor expand sin and arctan in Eq. ( 45), where θ is described by Eq. ( 42).In the second and third regions, where tan θ0 2 < y ≤ tan θ 0 + δ and ≤ y ≤ tan θ0 2 , we use the expression sin arctan tan θ 0 and Taylor expand around y = tan θ 0 and y = 0, respectively.Hence, sin arctan tan(θ 0 ) 1 y By substituting Eq. ( 51) into Eq.( 43), integrating and Taylor expanding the exponential, we find the leading-order terms of L C to be, , with the q i s being integers; all b i , p i , q i and H are defined in appendix F.
The stress resultant consists of three terms each representing a physical effect.Firstly, the cellulose contribution independent of the hotspots connections is: where c i are constants dependent on tan θ 0 or (defined in Table 2).Secondly, the cellulose contribution when crosslinked by the hotspots is: The generic form of the constants d j,i is where w 1 , w 2 and w 3 are known constants often dependant on tan θ 0 , δ and χ hot (see Table 3).Lastly, the contribution from the hotspots is then  , cellulose and hemicellulose stiffness densities ratio In the model, we will be using the relevant parameters for the roots of Arabidopsis thaliana.We assume k0 = that the breakage rates of the hotspots and the hemicellulose fibres at rest are approximately the same.We set a1 = 0.1 in most simulations as this number does not affect the results apart from the expansin action simulations where we shall comment on the affect of this parameter and its implications in Sect. 4 Further simplifications to the stress resultant can be made by examining the cases for small and large strain rate, α 1 and α 1, respectively.When α 1, all terms of the form Z 1 α that appear in the expansions are negligible, since Z is less than one for all such expressions.When α 1, all the terms of the form Equation ( 57) gives the leading-order terms that control the dominant cell wall behaviour for small and large strain rates, neglecting the effect of expansin action.This expression is used to validate the numerical scheme presented in Sect.4, where we also discuss the biological implications.

Results and discussion
The complexity of the plant cell wall and its growth process present many modelling challenges.We have therefore created a simple model to focus our attention on the cell wall structure and its implications on possible enzyme action.
The system of Eqs. ( 38)-( 44) is solved for the bond number and densities n and n hot , the CMF angle θ, and the stress resultant Σ, subject to a choice of expansin action where σ H is determined by one of the Eqs.( 34)- (37).Where the steady-state stress resultant is required, we use (45) instead of (44).The asymptotically simplified expressions ( 53)-( 57) for the steadystate stress resultant are used where applicable.This study focusses on the cellulose crosslinking dynamics; hence, we neglect the effect of pectin on the cell wall dynamics by setting Γ = 0 throughout.All other fixed parameter values are listed in Table 1.All solutions are obtained using MATLAB; in particular, the integrals are computed numerically using the Legendre-Gauss Quadrature code by Greg von Winckel [44].
We begin the discussion of the results by first analysing the model's implications on the cell wall structure and yield threshold in Sect.4.1.Secondly, we explore possible enzyme mechanisms and the likelihood 1 Page 12 of 30 Eur.Phys.J. E (2024) 47 :1 Fig. 3 The effect on the stress resultant Σ ∞ plotted against the strain rate α of expansin acting with rate E on a the resting length, L0, which decreases the stress resultant and thus the yield threshold; or b the viscosity, μ, where the small (large) a1 value represents reduced (enhanced) cellulose crosslinking by the hotspots and hence decreased (increased) the stress resultant and yield threshold.Other parameters: β hot = 0.01, γ = 0.05, and others found in Table 1 of hotspots structures in the cell wall in Sect.4.2.We finish with some final remarks in Sect.4.3.

Implications for cell wall structure
Before delving into the results, we outline the logic as to how the stress resultant links to the cell wall yield threshold and growth rate using expansin's effect as an example.A fixed turgor pressure in a growing cell wall (constant strain rate α) will result in a specific cell wall tension and therefore a set stress resultant.When expansin acts on the hotspots resting length, L 0 , the same cell wall stress/turgor pressure produces different strain rates (growth rate) (Fig. 3a, solid line versus dashed line); Eq. ( 1) then implies that the yield threshold has changed.We know this from the plateauing effect; the extensibility determines the gradient, but it is the yield threshold that scales the plateau up and down as α 1.Further analysis of the stress-strainrate relationship and the crosslinks influence on the cellwall yield threshold can be found in [35].For the purpose of this article, it suffices to know that a decrease in stress resultant implies a decrease in the yield threshold.
From the collective results presented here, we observe some general trends.As the strain rate increases, the stress resultant plateaus (e.g.Fig. 3).This relationship arises due to the bond number and densities (n, n hot ) behaviour when α is changed, despite the inclusion of linear elasticity.Specifically, for low values of α, the bond density n hot remains low for much of the domain (see Fig. 6b).However, as α is increased, the bonds become increasingly loaded and remain intact for longer (increasing the overall stress held in the cell wall) before finally breaking.Eventually the hotspots reach a limit as to how much stress they can withstand without breaking, and the bond density drops to 0 when y ≈ 0; the amount of stress held across the cell wall has reached capacity (Fig. 6b).Reaching this capacity consequently limits further increases to the stress resultant, implying that the cell wall's yield threshold has a maximum, where further increasing the growth rate has limited effect.This result demonstrates that cell wall strength is capped by the mechanism of cellulose crosslinkage, despite the cellulose fibres being the major load-bearing component in the cell wall.
We now examine the effect on the stress resultant of cellulose reorientation during growth.For lower strain rates, we find that the stress resultant increases as time progresses (Fig. 7a) due to a combination of all the fibres being increasingly strained and the cellulose reorientating.For larger strain rates, α, the stress has very little variation over time because the large reorientation rate allows fibres to converge to the same distribution.For smaller times, the stress has a sharp maximum (Fig. 7a).This peak occurs because increasing the strain rate stretches and reorients the fibres uniformly, due to the equivalent effect of α and t on the fibres from the term e αt .On the other hand, fibres deposited at t > 0 are progressively stretched and rearranged as they move to the outside of the wall.Thus, for certain times/strain rates (between α = 10 − 30), the sum of the stress arising in this family of fibres is less than those present at t = 0.This raised stress resultant does not last as the crosslinks present at t = 0 eventually break, causing a rapid decrease in the cell wall strength and producing a sharp peak in the stress (Fig. 7a).The effect of cellulose reorientation can also be seen when the initial fibre angle, θ 0 , is increased (Fig. 7b).According to Eq. ( 42), for any fixed fibre position y, the CMF increasingly reorient as θ 0 increases.Consequently, increasing θ 0 also significantly increases the stress resultant (Fig. 7b).From these results, we conclude that cellulose reorientation leads to a higher Fig. 4 The effect on the hotspots density, n hot through the cell wall thickness y of expansin with rate E acting on a the resting length, L0; or b the viscosity, μ.Parameters: α = 5, β hot = 0.01, γ = 0.05, and others found in Table 1.In both cases, expansin has increased the bond densities and thus has not weakened the cell wall stress resultant in the plant cell wall and, subsequently, increases the cell wall yield threshold.The yield threshold increasing due to cellulose reorientation could provide a possible mechanism for growth slow-down and hence the beginning of the cell's secondary wall structure.This hypothesis has been previously put forward [19].Cellulose reorientation reducing axial growth has been previously reported in other mathematical models [53,54].
We now proceed to analyse the implications of the asymptotic reduction of the stress resultant.There is substantial agreement between the asymptotic expansions and the numerical results, for all three combinations of β hot and γ (Fig. 8).There is little difference in accuracy between the full asymptotic Eqs. ( 53)- (56) and the simplified expansion (57), demonstrating that the cell wall yield threshold's behaviour is dominated by the terms in the simplified equation.According to Eqs. ( 53) and (54), the constants c i and d j,i in the cellulose's contribution to the cell wall strength depend on tan θ 0 (Table 2), whereby increasing tan θ 0 increases the stress held in the cell wall (Fig. 7).Since a 2 < a 1 (see Table 1), we conclude that the main control of the cell wall yield threshold in the axial direction is the cellulose orientation.This result matches past simulation results [17] and the current consensus that the cellulose orientation controls the growth direction in most plant cells [21].
In Eq. ( 57) that characterises the cell wall stress behaviour, the hotspots density's contribution (the χ hot term) is controlled by b 3 .This coefficient determines the slope of the extension curve of the cellulose molecules, L C , as they rapidly extend close to y = 0 where they are most stressed.Thus, the effect on the yield threshold of breaking the hotspot crosslinks (changing β hot ) is predominantly actioned through loosening the smaller group of significantly stressed cellulose molecules.Equation ( 57) also demonstrates why the stress resultant plateaus (Fig. 3).Rearranging the equation gives As α → ∞, the first term on the first line remains constant.The second term can be recognised and reexpressed to be a limit of the form lim x→0 and therefore is a constant.The last term also converges to a constant (with α also converging to a constant) for the same reason.Collectively this means that as α increase the stress resultant will level out and become constant.We also note that for α 1, the stress resultant is not controlled by χ and χ hot at all.This implies for low levels of turgor pressure, enzyme action on the cell wall is ineffective at controlling the growth rate and that it is solely controlled by the orientation and concentration of fibres.

Enzyme action on the cell wall
We first discuss how our simulations can be directly compared to experimental results relating to enzyme action.Experiments have shown upon overexpressing or silencing expansins, the growth (through elonga- tion and division) of plant roots and leaves increase or decrease respectively in a variety of species [55,56], and root hair and leaf primordia initiation are also affected [57].Overexpression and knockout mutations were also performed on XTH proteins which made little difference to growth [58,59].Moreover, inactivating and then extending plant tissues while applying Cel12A or expansins permitted cell wall extension [10,28].The effect on growth and hence the strain rate of enzyme action can be seen in our model through the change in stress resultant.A decrease in the stress resultant will lead to an increased growth rate (strain rate) and vice versa as described in Sect.4.1.Additionally, from the extension experiments with multiple deformations (elastic and plastic) and examining the stress relaxation response, expansin was shown not to weaken the wall while Cel12A did.Additionally, XTH was found to have only minor effects on cell wall mechanics [58].This effect can be directly observed in the model by examining the bond density distributions; if there are fewer bonds intact, the wall integrity will be weaker.
We now examine different modes of enzyme action beginning with expansin.As previously mentioned in Sect. 1, expansin causes wall loosening without reducing the strength of the cell wall, potentially by allowing the slippage of fibres within the hotspot.We simulated expansin's hypothesised action via two methods: increasing the spring resting length and decreasing the viscosity of the dashpot by equations characterised in Sect.2.2.3 (Eqs.35 and 37).Simulations of expansin's theorised fibre slippage action on the resting length (Eq.35) reveal a decrease in the stress resultant (Fig. 3a), while ensuring the bond density does not reach zero (Fig. 4a).Similarly, expansin's slippage action on the viscosity of the dashpot (Eq.37) decreases the stress resultant provided hotspots cellulose crosslinking is small (Fig. 3b) and always increases the bond density (Fig. 4b).This bond density increase arises due to stress relaxation; both mechanisms reduce the elastic strain imposed on the fibres which are therefore less likely to break.Two opposing effects contribute to the overall stress resultant: stress relaxation from the hotspots, combined with the inclusion of more, increasingly stretched, CMF molecules towards the outer boundary of the cell wall.When the hotspotscellulose crosslinking is sufficiently large to counteract the hotspots stress relaxation, the stress resultant increases with expansin action due to the increased contribution of the CMF (Fig. 3b).Therefore, both mechanisms can decrease the axial stress resultant for all strain rates (Fig. 3) and are consequently effective at inducing growth.Moreover, the hotspots bond densities n hot increase for both enzyme action pathways (Fig. 4), increasing cell wall integrity and thus strengthening the cell wall.So, allowing the fibres to slide past one another (the dashpots) or relax (increasing resting lengths) loosens the fibre network and decreases the likelihood of the bonds breaking.This explains the observed experimental behaviour (as described in Sect. 1) and provides evidence for expansin's hypothesised slippage action.
In contrast, when the parameter a 1 is increased the hypothesised expansin action causes the yield threshold to increase (Fig. 3b) despite the loosening of the fibre network.This observation may be consistent with biological observations.It has been observed that plants lacking xyloglucan have only a small amount of growth reduction and alteration in cell wall structure [14,60].It has therefore been suggested that any xyloglucan present is concentrated in these compact biomechanical hotspots.Therefore, as a 1 represents the proportion of xyloglucan in the wall, we find that in order to recreate expansin's experimental behaviour, there must be a small amount (comparatively to the cellulose) of xyloglucan that is concentrated around the hotspots.
Continuing the investigation of hypothesised enzyme action, we now focus on XTH and Cel12A.As detailed in Sect. 1, XTH and Cel12A are theorised to hydrolyse hemicellulose and the hotspots respectively.Recall that we recreate this cutting action by increasing the breakage rate γ of the hemicellulose fibres to model XTH, and increasing the hotspot breakage rate β hot to model Cel12A action.By testing a range of values for both parameters, we discover that β hot significantly affects the stress resultant compared to γ (Fig. 5a), (Cel12A's cellulose digestion effect is investigated in appendix E where it simply increases its effectiveness on the stress resultant).These results agree with previous work, which also finds the yield threshold to be controlled by the rate of dissociation of tethering crosslinks [61].The effect of γ on the hotspot bond density, n hot , is also negligible (Fig. 5b).In contrast, an increased β hot leads to a considerable decrease in n hot (Fig. 6a).As the fibres move towards the outer y = 0 boundary, their lengths increase sharply (L = 1 y ), leading to significantly larger stress; the stress resultant is therefore highly sensitive to small changes in the transition location χ hot in the region 0 < y < 0.1 where n hot drops rapidly.Thus, increasing the breakage rate of the hemicellulose fibres (increasing γ) is ineffective when compared to targeting the hotspots themselves (increasing β hot ).
Simulated XTH action generates minimal reductions in the stress resultant and the hotspots density, implying that minimal growth is induced due to its limited wall-loosening ability.This effect occurs despite a natural assumption that the hydrolysis of bonds by some members of the XTH family could cause wall loosening [10,28].The model outputs offer two interpretations to explain this phenomenon: either the hotspots maintain their integrity when hemicellulose fibres are cut, or hydrolysing hemicellulose could be ineffective at breaking hotspots down if the hotspots are compact and inaccessible.Thus, the limited wall loosening is a consequence of the cell wall retaining its integrity, with minimal increase in cell wall growth as observed in experiments [10].In the asymptotic expansions (Eqs.( 57)), it is β hot , not γ (in the form of χ hot ), that controls the leading-order terms of cellulose's contribution.We therefore explicitly see that if XTH's ability to break  1 Fig. 6 The hotspot density, n hot , parameter dependence.a The effect of varying Cel12a's action on n hot with α = 10 demonstrating its effectiveness at decreasing the hotspot density.b The effect of varying α with β hot = 0.01 on n hot , both with γ = 0.05, showing the transitional behaviour of increasing α and the eventual plateauing under further increase.All other parameters can be found in Table 1 down the hotspots is ineffective, the observed experimental behaviour is reproduced.
By recreating Cel12A's hypothesised hotspot-cleaving action, we see that it weakens the cell walls by decreasing the hotspots density, which decreases the stress resultant and induces growth, matching biological experiments [10].This effect would not occur if a tethered network was assumed, in which case the model would only have one bond density that is not dependent on another (in a similar manner to [35]), leaving us unable to distinguish between Cell12A and XTH's action.Through analysis of the breakage rates and the model's assumption of the hotspots hypothesis, we have recreated experimentally observed enzyme behaviour.

Remarks
The results presented play into the hypothesis [27] first stated that growing plant cell walls seem to undergo a glass transition.This theory suggests that hemicellulose connecting the cellulose molecules is not the primary stress-bearing component (as a 1 is small) and holds the 1 Page 16 of 30 Eur.Phys.J. E (2024) 47 :1 Fig. 7 Changes in the stress resultant as we effect the cellulose fibre orientation by changing a the time point and b the cellulose angle at deposition, θ0 with β = 0.01 and γ = 0.05.Both graphs demonstrate that as the cellulose fibres increasingly reorient, they increase the stiffness and yield threshold in the axial direction.All other parameters can be found in Table 1 Fig. 8 Comparison between the axial stress resultant's numerical solution (solid line) and asymptotic approximations, Eqs. ( 53)-(56) (dashed line) and (57) (dotted line), at various values of β hot and γ, demonstrating the accuracy of the asymptotic expansions.All other parameters can be found in Table 1 wall in tension close to the yield threshold, such that when the stress relaxation imposed on the hotspots, by digestion or slippage occurring, it could be enough to cause the redistribution of stress to other cell wall components such as pectin (and its linkers with cellulose) and hence induce growth.Thus, we have a rapid change in viscosity with the extra load causing the pectin to 'melt' and flow.As the cell wall extends, other hotspots could become strained, hence raising the yield threshold and requiring further relaxation to continue growth.
The results in this article depend on the inclusion of the cellulose contributing to the axial stress.Even though, on average, the cellulose molecules are perpendicularly oriented to the growth axis [22], the findings show the importance of including cellulose reorientation or a distribution of angles in mathematical models of growth.
There is still further modelling work to be done.We have only examined axial growth, so it remains to be understood how enzyme action does not necessarily lead to radial growth and the consequences of cellulose reorientation, as it has been observed a reduction in anisotropy can lead to radial swelling [62].Expansin has been shown to induce growth as soon as it is applied, while Cel12A-induced growth is delayed after application [12].Our model cannot explain this time delay phenomenon, meaning we could benefit from further work on the crosslinking dynamics to understand this process.Many experiments are done on a tissue level, while our model examines only a small cell wall section.Therefore, there could be advantages to constructing a tissue-scale model for cells, including their interactions and the different cell layers to match with experiments.Models such as [17] also have great potential to aid our understanding but currently neglect growth and enzyme manipulation.The stress-strain relationship from this model could be used to inform constitutive laws and then incorporated into a larger-scale model with growth and wall modifications.

Conclusion
By deriving and analysing a mathematical model of the cell wall incorporating the biomechanical hotspot dynamics as proposed by [9], we have explained experimentally observed enzyme behaviour and thus provided insight into the cell wall structure and loosening mechanisms.We propose the following enzyme mechanisms in line with experimental observations (Fig. 9).XTH only cuts individual hemicellulose fibres and thus does not provide sufficient wall loosening to induce growth, possibly due to the hotspots remaining inaccessible or hemicellulose scissoring being insufficient to fully break down the hotspot.On the other hand, Cel12A can cleave the whole hotspot allowing the cell to grow while also weakening the cell wall.Expansin action causes fibre slippage, permitting stress relaxation, enabling cell wall extension without weakening the cell wall.We have therefore presented a mechanistic understanding of potential enzyme action.Finally, cellulose reorientation during cell wall extension leads to a decrease in the axial growth rate and its range orientation is the primary control for cell growth and the yield threshold.
Enzyme action on the cell wall is an under-researched area of plant growth often overlooked by modellers.The theory presented in this paper generates improved understanding of the fundamental mechanisms under-lying plant cell growth.We hope that it provides a key building block towards a unified framework for plant development.

(B2)
If we shift the time parameter so that t = t − τ where the fibre is deposited at time t = τ , we have that σH (1, 0) = 0.The above then becomes, Using an integrating factor and integrating both sides, we derive, where on substituting L = e α t (using boundary conditions L(0) = 1) we then have,

(B5)
The fibres follow the trajectories y = e −α t (from Sect. 2 using boundary condition y(0) = 1) and with y = e ω t − α ω = 1 L , we calculate the final form as For the expansin action on the spring resting length (Eq.( 32)), we were unable to find an analytic expression, as a result we use forward finite differences to solve the equation.We need this expression in terms of y so we first divide it by dy

(B8)
We can then use the calculated values of L0 in Eq. ( 32) to find the derivative of L0.We now rederive the hotspot stress, σH with non-constant resting length using the initial equation (Eq.25), shifting time again and the fact that and using the same method as above and Eq.B8 we calculate, (B10) which is the hemicellulose stress term for enzyme action on the resting length, L0.We repeat the same procedure for expansin action on the impedance, μ.By solving Eq. ( 31), we can obtain, (B11) . Substituting this into the nondimensional stress term (Eq.( 25)) and shifting time we have, where we again use an integrating factor and L = e α t to find Upon integrating we then get where we then make the substitution ŷ = e −α t in the integral.This implies that e −Et = y E α and dt = − 1 αy dy meaning

Appendix C: Solving the evolving bond densities
To determine the evolving bond number and densities we solve (30) for n, using the method of characteristics with dy dt = −αy, dx dt = αx, and ∂n ∂x = 0 meaning the left-hand side of Eq. (30)

(C21)
To find the hotspot bond density, n hot (Eq.( 29)), we repeat the same procedure using the total derivative and again using Eq. ( 24) we derive

Appendix D: Solving the cellulose angle equation
To find the angle cellulose makes to the horizontal, θ, we solve Eq. ( 28), by first re-expressing left-hand side as the total derivative to get dθ dt = α sin θ cos θ on dy dt = −αy.

(D30)
For the cellulose length we once again use the transformed co-ordinate t to get, using the separation of variables we get (D32) For our different domains we then get, and (D34)

Appendix E: Cel12A cellulose action
To model Cel12A's digestion action on cellulose, we assume it reduces the cellulose density ρ.For simplicity, this will take the form.
where C * is Cel12A's action rate and is proportional to β hot , meaning C * = C * β hot .We take this form to ensure that ρ * = 0 and to limit digestion.In the non-dimensional system, the above becomes which will have the same value as a2 shown in Table 1.Equation E36 can be solved to get 1 Page 20 of 30 Eur. Phys.J. E (2024) 47 :1 Fig. 10 Enzyme cutting action.The effect on the stress resultant of Cel12a with its additional cellulose digestion effect and XTH with α = 50, C = 100 and N = 10.Cel12a's action, β hot and its degradation, on the whole hotspot is much more effective at decreasing the stress resultant than XTH's action, γ, on the hemicelluloses.All other parameters can be found in Table 1 The steady-state stress result then becomes σH (y) dy + Γα.

(E38)
Figure 10 shows the effect of Cel12A additional cellulose digestion where it simply decreases the stress resultant (significantly) more when compared to Fig. 5 where it just affects the hotspots.

Appendix F: Asymptotic expansion
Beginning with the asymptotic expansion of G, we consider the two regions, γ y ≤ 1 and ≤ y γ.In considering the first case, we know that 1 z − z (F41) Using integration by parts and knowing that (F42) Using the fact that γ is small, we can ignore the first term, and if we then repeatedly integrate by parts, we can get the sequence (F44) From Eq. ( 38), the asymptotic approximation of n is then,

(F45)
Writing k0 = k 0 k 0,hot , we now calculate the asymptotic expansion of G hot .To find the approximation, we begin in the region γ y ≤ 1 and β hot y ≤ 1, where upon using the expansion of n in this same region (Eq.( 47

(F48)
Lets consider any term in the series, 1 − z k0 α i , with i ∈ N/{0}.We want to find the biggest term in this bracket.We begin by assuming that 1 is the biggest term and using the fact that k0 α > 0 and 0 < z ≤ 1 we have, So 1 is indeed the bigger term meaning 1−z k0 α ≤ 1.Now we need to show that the 1st term in the integral (Eq.(F48)) is greater than for each i ∈ N/{0}.From the workings above, we know the ith term is less than, so on assuming the first term, 1 is bigger than the other terms for any i, we have α2 and since β hot z 1, we now know 1 is indeed the biggest term in the sequence.Therefore, we can ignore the higher-order terms in Eq. (F48) to get G hot (y) ≈

(F50)
Deriving the expansion for < y β hot and < y γ we have that, (F52) The case β hot < y < γ is a hybrid of the two methods.You can ignore the exponential term in a similar way to (F51) and follow the same derivation for G in Eq. (F40) to derive G hot = ln 1 y .Collecting the results together we get,

(F53)
Notice that the simplification of n hot is not dependant on n.
We need to know exactly where the different asymptotic expansions in the different regions can be applied.Starting with n, the first boundary region, < y β hot , applies when n has rapidly decreased.So if we find the value for y for when this happens, we can get out a critical value for when the different expansions can be used.To find this value, we use the derivative of the dominate term of the exponent of n which is k0 α exp γ 2 α y 2 and find the value when it exceeds a certain number, as when this rapidly increases, n rapidly decreases.This value is then, where Q is a threshold value of our choice.The threshold value for n hot , similarly is, .

(F55)
To be able to solve the limit stress resultant, we need to also approximate the trigonometric term present, sin arctan tan(θ0) (F66) The approximation of the trigonometric term is then defined by three different expansions in three different regions, such that, sin arctan tan(θ0) Where tan θ 0 2 and tan θ0 + δ are interval values of our choosing to ensure accuracy.We can use this expression for the cellulose angle to find the leading terms that dominate the celluloses length's, LC behaviour.Substituting in (51), we get where ∼ is used for simplicity to represent the integrand that appears in the previous line for the same integral region.Upon evaluating our integral, we find, (F72) When we ignore the higher-order terms and expanding Taylor expanding the exponential we get the final form as, where b1 = e B1 (1) , ) , and (F76) We have all the needed approximations to calculate the stress resultant analytically.We start of with the cellulose contribution independent of the biomechanical hotspots such that,  (F78) which we find to be where all the ci values can be found in Table 2.
The second term is the contribution from the cellulose fibres when crosslinked by the biomechanical hotspots which we calculate to be, where the values of the coefficients d1,i, . . .,,d2,i, . . .,, d3,i, . . .,, and d4,i, . . ., can be found in Table 3.The third term is the hemicellulose contribution and is found to be, where upon integration Note that for all three contributions to the stress resultant, Σ ∞ 1 , Σ ∞ 2 , and Σ ∞ 3 for some critical values of α some of their terms are logarithmic upon integration which we choose to ignore.

Fig. 2
Fig. 2 On the left a depiction of a single unit of the continuum model with the cellulose (box A) and hotspot (box B) contributions to the stress joined in parallel.The hotspot stress σ * hot arises from a collection of hemicelluloses represented as Maxwell elements, consisting of a spring and viscous damper/dashpot, crosslinked with bond number n * (where the red cross represents a bond breakage).The cellulose stress σ *C is represented by a box that resists the strain elastically proportionally to the cellulose's angle θ.

where b 0
, b 1 , b 2 and b 3 are constants and H(y) is a function of the form i p i y tan θ0 qi

Fig. 5
Fig. 5 Enzyme cutting action.a The effect on the stress resultant, Σ ∞ of Cel12a and XTH with α = 50.Cel12a's action, β hot , on the whole hotspot is much more effective at decreasing the stress resultant than XTH's action, γ, on

Fig. 9
Fig. 9 Hypothesised enzyme action.The purple rods represent cellulose molecules, and the red rods are the hemicellulose fibres in a hotspot.a XTH hydrolysis action only cutting a select amount of hemicellulose fibres.b) Cel12A hydrolysis action cutting the whole hotspot.c Expansin loosening action allowing the hotspot to be intact while allowing extension can once again use the upper bound of,

Table 2
Asymptotic expansion constants for Σ ∞