<figure><img src="image_1.png" /></figure>

HAL
open science

# Stoichiometry of irreversible ligand binding to a one-dimensional lattice

Philipp A Tsvetkov

<figure><img src="image_2.png" /></figure>

## To cite this version:

Philipp A Tsvetkov. Stoichiometry of irreversible ligand binding to a one-dimensional lattice. Scientific Reports, 2020, 10 (1), pp.21308. 10.1038/s41598-020-77896-0. hal-03104103

HAL Id: hal-03104103

https://amu.hal.science/hal-03104103

Submitted on 19 Jun 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_3.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_4.png" /></figure>

OPEN

# 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¹, but are also applicable for analyses of polysaccharide² and supramolecular polymers such as microfilaments and microtubules³. The equilibrium ligand binding to a one-dimensional lattice is usually described by the McGhee-von-Hippel equation⁴. 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⁵. 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 $\ell$ binding sites after saturating the regions with 1, ..., $\ell-1$, the remaining free binding sites will stay unoccupied because they cannot accept $\ell$-long ligand. The probability of the lattice becoming fully covered by ligands decreases dramatically with the decrease of the lattice length (A231580 in⁸). 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 $\ell$-long ligand. Let $a_{l,n}$ be the expected fraction of sites occupied by the $\ell$-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 (Fig. 1) could be found by modifying the difference equation used by Gordon⁶ for the number of vacant sites, which is in turn a generalization of the simplest case of 2-long ligand⁷:

CNRS, UMR 7051, INP, Inst Neurophysiopathol, Fac Médecine, Aix Marseille Univ, Marseille, France. email: philipp.tsvetkov@univ-amu.fr
<figure><img src="image_5.png" /><figcaption>**Figure 1.** The random filling of n-length lattice by ℓ-length segments. First segment (shown in light blue) has an equal probability to occupy one of the n − ℓ + 1 position shown in figure. The expected number of occupied sites after the complete filling of the lattice for each configuration will be equal to ℓ plus the expected number of occupied sites for lattice on the right and on the left from the first segment.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Figure 2.** (A) Fraction of occupied sites on n-length lattice for different values of ligand length (ℓ). Theoretical curves are shown in solid lines and numerical simulation values (see “Method” section for details) in squares with error bars. (B) The fraction of occupied sites on infinite lattice as a function of ligand length is shown in black squares and best fit with Eq. (8) as red curve. (C) Fraction of occupied sites of n-length lattice for different values of ligand length (ℓ) from 2 to 10. (D) Stoichiometry of ℓ-length ligand binding on infinite lattice in case of equilibrium (black dots) and irreversible (red dots) intercation.</figcaption></figure>

$$a_{l,n} = \frac{l + 2(n-l)a_{l,n-l} + (n-1)(n-l)a_{l,n-1}}{n(n+1-l)}. \qquad (1)$$

The dependences between the fraction of the occupied sites and the length of the lattice (n) is presented in Fig. 2A,C for different values of ligand length (ℓ). Theoretical curves derived from Eq. (1) are shown by solid
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^l = \{a_{l,0}, a_{l,1}, a_{l,2}, \dots\}$ to find the fraction of sites occupied by a $\ell$-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 $\ell$ to $\infty$:

$$ \sum_{n=l}^{\infty} n(n+1-l)a_{l,n}x^n = \sum_{n=l}^{\infty} lx^n + \sum_{n=l}^{\infty} 2(n-l)a_{l,n-l}x^n + \sum_{n=l}^{\infty} (n-1)(n-l)a_{l,n-1}x^n \quad (2) $$

If we define $A_l(x)$ as $\sum_{n=0}^{\infty} a_{l,n}x^n$ and take into account that $A'_l(x) = \sum_{n=1}^{\infty} na_{l,n}x^{n-1}$ and that $A''_l(x) = \sum_{n=2}^{\infty} n(n-1)a_{l,n}x^{n-2}$ we can rewrite Eq. (2) as differential equation:

$$ (x^2 A''_l(x) - (2-l)x A'_l(x))(x-1) + 2x^{l+1} A'_l(x) + \frac{lx^l}{1-x} = 0 \quad (3) $$

It is easy to demonstrate that the following function is the solution of the obtained differential Eq. (3):

$$ C_1 + \int_1^x \left( \frac{e^{-2\mu(t,l)} t^{l-2}}{(t-1)^2} \left( C_2 + \int_1^t le^{2\mu(k,l)} dk \right) \right) dt, \quad (4) $$

where $\mu(x, l) = \sum_{n=1}^{l-1} \frac{x^n}{n}$ and $C_1$ and $C_2$ are constants.

Taking into account that $\frac{A_l^{(n)}(0)}{n!} = a_{l,n}$ and that $a_{l,n} = 0$ for $n < \ell$, we can find both $C_1$ and $C_2$ and finally find $A_l(x)$:

$$ A_l(x) = \int_0^x \left( \frac{e^{-2\mu(t,l)} t^{l-2}}{(t-1)^2} \int_0^t le^{2\mu(k,l)} dk \right) dt. \quad (5) $$

This analytical solution allows to calculate the fraction of the occupied sites by the $\ell$-long ligand on the n-long lattice ($a_{l,n}$) for any $\ell$ and n without using a recurrence sequence. For example, for 2-long ligand the fraction of occupied sites is equal to:

$$ a_{2,n} = \frac{(-2)^{1+n} e^2 + \Gamma(3+n, -2)}{n e^2 \Gamma(2+n)}, \quad (6) $$

where $\Gamma(n)$ and $\Gamma(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:

$$ a_{l,\infty} = le^{-2\mu(1,l)} \int_0^1 e^{2\mu(t,l)} dt, \quad \text{or} \quad (7) $$

$$ a_{l,\infty} = le^{-2H_{l-1}} \int_0^1 e^{2\mu(t,l)} dt, \quad (7) $$

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 substituents of vinyl polymers by Paul J. Flory⁷, who found the fraction of the occupied sites on a randomly filled infinite lattice to be equal to $1-e^{-2 \approx 0.865}$ (A219863 in OEIS⁸). 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⁶. In case of infinitely long ligand this problems turns into a continuous parking problem and have been solved by Renyi⁹. 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⁸). 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⁸) and ≈ 0.804 (A307184 in OEIS⁸) respectively. For ligand lengths equal to 2 or 3 this fraction can be expressed through closed formulae ($1-e^{-2}$ and $\frac{3\sqrt{\pi}}{2e^4}$ ($erfi(2) - erfi(1)$), respectively), and for larger ligand lengths can be estimated numerically from (7). The dependence of $a_{l,\infty}$ as a function of the ligand length $\ell$ is presented in Fig. 2B. The Eq. (7) can be also approximated for any $\ell$ (Fig. 2B) with the following formula:

$$ a_{l,\infty} \approx 0.748 + \frac{0.228}{l} \pm 10^{-3}. \quad (8) $$

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 $\ell$-long ligand in case of an equilibrium binding (at saturating concentrations) and irreversible
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 l-length ligand was simulated using Wolf-ram Mathematica 10.1 software (https://www.wolfram.com/mathematica/) for 19<n<50 and 1<l<11. Each simulation was run 100 times then average value and standard deviation were calculated for all values n and l.

Received: 29 June 2020; Accepted: 23 October 2020

Published online: 04 December 2020

## References

1. Latt, S. A. & Sober, H. A. Protein-nucleic acid interactions. II. Oligopeptide-polyribonucleotide binding studies. *Biochemistry* **6**, 3293-3306 (1967).

2. Olson, S. T., Halvorson, H. R. & Björk, I. Quantitative characterization of the thrombin-heparin interaction. Discrimination between specific and nonspecific binding models. *J. Biol. Chem.* **266**, 6342-6352 (1991).

3. Tsvetkov, P. O., Makarov, A. A., Malesinski, S., Peyrot, V. & Devred, F. New insights into tau-microtubules interaction revealed by isothermal titration calorimetry. *Biochimie* **94**, 916-919 (2012).

4. McGhee, J. D. & von Hippel, P. H. Theoretical aspects of DNA-protein interactions: Co-operative and non-co-operative binding of large ligands to a one-dimensional homogeneous lattice. *J. Mol. Biol.* **86**, 469-489 (1974).

5. Munro, P. D., Jackson, C. M. & Winzor, D. J. On the need to consider kinetic as well as thermodynamic consequences of the parking problem in quantitative studies of nonspecific binding between proteins and linear polymer chains. *Biophys. Chem.* **71**, 185-198 (1998).

6. Gordon, M. & Hillier, I. H. Statistics of random placement, subject to restrictions, on a linear lattice. *J. Chem. Phys.* **38**, 1376-1380 (1963).

7. Flory, P. J. Intramolecular reaction between neighboring substituents of vinyl polymers. *J. Am. Chem. Soc.* **61**(6), 1518-1521 (1939).

8. Sloane, N. J. A. *The On-Line Encyclopedia of Integer Sequences* http://oeis.org/.

9. Rényi, A. On a one-dimensional problem concerning random place filling. *Magyar Tud Akad Mat Intézeténik Kozleményei* **3**, 109-127 (1958).

## Author contributions

P.O.T. wrote the main manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

Correspondence and requests for materials should be addressed to P.O.T.

Reprints and permissions information is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_7.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2020