Stoichiometry of irreversible ligand binding to a one-dimensional lattice

In this paper we investigate the problem of irreversible binding of a ligand that covers several identical binding sites on a macromolecule with a one-dimensional lattice. Due to steric constraints, irreversible binding or binding with slow kinetics results in partial saturation of the binding sites thus impacting the stoichiometry of the interaction. Here we present a recursive formula to calculate the exact fraction of the occupied binding sites for a ligand and macromolecule of arbitrary lengths. We also provide an analytical result for the exact fraction of the occupied sites in case of an infinitely long lattice. We conclude with a simplified empirical formula for the exact fraction of the occupied sites in case of an infinitely long lattice.


Stoichiometry of irreversible ligand binding to a one-dimensional lattice Philipp O. Tsvetkov
In this paper we investigate the problem of irreversible binding of a ligand that covers several identical binding sites on a macromolecule with a one-dimensional lattice. Due to steric constraints, irreversible binding or binding with slow kinetics results in partial saturation of the binding sites thus impacting the stoichiometry of the interaction. Here we present a recursive formula to calculate the exact fraction of the occupied binding sites for a ligand and macromolecule of arbitrary lengths. We also provide an analytical result for the exact fraction of the occupied sites in case of an infinitely long lattice. We conclude with a simplified empirical formula for the exact fraction of the occupied sites in case of an infinitely long lattice.
One-dimensional lattice-like models of interaction are used to study the ligand binding to regular macromolecules containing long repetitive sequences of binding sites (binding lattice). Usually, such models are used to study DNA/RNA-ligand interactions 1 , but are also applicable for analyses of polysaccharide 2 and supramolecular polymers such as microfilaments and microtubules 3 . The equilibrium ligand binding to a one-dimensional lattice is usually described by the McGhee-von-Hippel equation 4 . Unfortunately, due to the slow kinetics of the interaction, a transient (pseudo-equilibrium) state could be misinterpreted as a true equilibrium. In extreme case of irreversible ligand binding, this problem can be described in terms of the so-called car parking problem. The likelihood of the problem occurring also increases with the increase in length of the acceptor lattice 5 . In this paper we investigate how the parking problem constraints influence the stoichiometry of ligand-macromolecule interactions in the case of the irreversible binding. We then derive a closed form solution for an extreme case where the length of the lattice increases infinitely.
If a ligand occupies only one site on the lattice and does not influence the other binding sites (non-cooperative binding), its interactions can be described by a simple model of binding with the stoichiometry 1:1. However, if the ligand occupies several neighboring binding sites on the macromolecule, the full saturation cannot be reached. Indeed, ligands binding to available sites at random can lead to formation of areas on the lattice where free sites cannot be occupied by the ligand molecule due to steric constraints (e.g., gap size is smaller than that of the ligand). For example, if the ligand occupies ℓ binding sites after saturating the regions with 1, …, ℓ − 1, the remaining free binding sites will stay unoccupied because they cannot accept ℓ-long ligand. The probability of the lattice becoming fully covered by ligands decreases dramatically with the decrease of the lattice length (A231580 in 8 ). This means that the stoichiometry of such interactions (the maximal ratio between bound binding sites and ligand) will always be higher than expected. In this study we do not consider models of the ligand binding to the lattice but determine the stoichiometry of such interactions through the ligand-to-binding-sites molar ratio after the saturation is reached.

Results
To solve this problem, we will consider a n-long lattice randomly filled with a ℓ-long ligand. Let a l,n be the expected fraction of sites occupied by the ℓ-long ligand at the n-long lattice. Obviously, if the length of the lattice is less than the length of the ligand then the fraction of occupied sites is equal to zero (a l,n = 0, forn < l), and when the lengths of the lattice and ligand are equal then all sites will be occupied a l,n = 1, forn = l . For a larger lattice lengths (n > l), the fraction of the occupied sites (   www.nature.com/scientificreports/ lines and simulation values (when a n-long lattice is randomly filled with ℓ-long ligand 1000 times) are shown by squares with error bars attached. Together with known a l,n for n < l this relation could be used to generate recurrence sequence A ℓ = {a l,0 , a l,1 , a l,2 , …} to find the fraction of sites occupied by a ℓ-long ligand binding to the lattice of any length. Let's solve this recurrence using generating functions approach. For this, we multiply the Eq. (1) by x n and sum for all n from ℓ to ∞: If we define A l (x) as ∞ n=0 a l,n x n and take into account that A ′ l (x) = ∞ n=1 na l,n x n−1 and that A ′′ l (x) = ∞ n=2 n(n − 1)a l,n x n−2 we can rewrite Eq. (2) as differential equation: It is easy to demonstrate that the following function is the solution of the obtained differential Eq. (3): where µ(x, l) = l−1 n=1 x n n and C 1 and C 2 are constants.

Taking into account that
= a l,n and that a l,n = 0 for n < ℓ, we can find both C 1 and C 2 and finally find A l (x): This analytical solution allows to calculate the fraction of the occupied sites by the ℓ-long ligand on the n-long lattice ( a l,n ) for any ℓ and n without using a recurrence sequence. For example, for 2-long ligand the fraction of occupied sites is equal to: where Ŵ(n) and Ŵ(n, z) are the Euler gamma function and the incomplete gamma function respectively. It also allows to find analytical solution for an extreme case when the lattice length approaches infinity: where H l is a harmonic number.

Discussion
The car parking problem could considerably impact the stoichiometry of the irreversible (in some cases also reversible) ligand binding to the lattice of binding sites. This problem was theoretically investigated about 80 years ago for a 2-long ligand in the context of intramolecular reaction between neighbouring subtituentes of vinyl polymers by Paul J. Flory 7 , who found the fraction of the occupied sites on a randomly filled infinite lattice to be equal to 1−e −2 ⋍ 0.865 (A219863 in OEIS 8 ). Further, in 1963 Gordon and Hillier extended this approach to ligands of any length but could not find an analytical solution for an infinite lattice 6 . In case of infinitely long ligand this problems turns into a continuous parking problem and have been solved by Renyi 9 . In this case the fraction of the occupied space called also the Rényi's parking constant is minimal and equal approximately to 0.748 (A050996 in OEIS 8 ). General formula (7) presented here allows to calculate the fraction of sites occupied on an infinite lattice by a fixed-length ligand (Fig. 2B). For example, for a ligand lengths equal to 3 and 4, this fraction is equal to ⋍ 0.824 (A307154 in OEIS 8 ) and ⋍ 0.804 (A307184 in OEIS 8 ) respectively. For ligand lengths equal to 2 or 3 this fraction can be expressed through closed formulae (1−e −2 and 3 √ π 2e 4 erfi(2) − erfi(1) , respectively), and for larger ligand lengths can be estimated numerically from (7). The dependence of a l,∞ as a function of the ligand length ℓ is presented in Fig. 2B. The Eq. (7) can be also approximated for any ℓ (Fig. 2B) with the following formula: This approximation is applicable for any ligand length in case of a very long lattice. To calculate the precise stoichiometry of the lattice of any length the recurrence Eq. (1) can be used. The difference in the number of sites occupied by a ℓ-long ligand in case of an equilibrium binding (at saturating concentrations) and irreversible (2) ∞ n=l n(n + 1 − l)a l,n x n = ∞ n=l lx n + ∞ n=l 2(n − l)a l,n−l x n + ∞ n=l (n − 1)(n − l)a l,n−1 x n  www.nature.com/scientificreports/ binding can be seen in Fig. 2D. For example, a 5-long ligand can occupy 6.3 binding sites on average. Therefore, finding an unexpectedly high stoichiometry (using isothermal titration calorimetry, for example) could be an indication of an irreversible binding or a state, close to pseudo-equilibrium formula presented here allows to find the correct stoichiometry.
In conclusion, here we present the exact analytical formula to calculate the fraction of occupied sites when an infinite lattice randomly filled with ligands of a fixed length (7). Additionally, we present an empiric expression that allows us to estimate this fraction with high precision.

Method
Fraction of occupied sites on n-length lattice after random fill with ℓ-length ligand was simulated using Wolfram Mathematica 10.1 software (https ://www.wolfr am.com/mathe matic a/) for 19 < n < 50 and 1 < ℓ < 11. Each simulation was run 100 times then average value and standard deviation were calculated for all values n and ℓ.