# Kinetic proofreading of chromatin remodeling: from gene activation to gene repression and back

Raghvendra Singh, Guillaume Brysbaert, Marc Lensink, Fabrizio Cleri, Ralf Blossey

► To cite this version:

Raghvendra Singh, Guillaume Brysbaert, Marc Lensink, Fabrizio Cleri, Ralf Blossey. Kinetic proofreading of chromatin remodeling: from gene activation to gene repression and back. AIMS Biophysics, 2015, 2 (4), pp.398-411. 10.3934/biophy.2015.4.398. hal-03285276

HAL Id: hal-03285276

https://hal.univ-lille.fr/hal-03285276

Submitted on 13 Jul 2021

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_1.png" /></figure>

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

## Review

# Kinetic proofreading of chromatin remodeling: from gene activation to gene repression and back

Raghvendra P. Singh¹,², Guillaume Brysbaert¹, Marc F. Lensink¹, Fabrizio Cleri², and Ralf Blossey¹*

¹ University of Lille, CNRS UMR 8576 Unité de Glycobiologie Structurale et Fonctionnelle (UGSF) and FRABio FR3688, 59000 Lille, France

² University of Lille, CNRS UMR8520, Institut d'Electronique, de Microélectronique et de Nanotechnologie (IEMN), 59000 Lille, France

* Correspondence: ralf.blossey@univ-lille1.fr; Tel: +33-3-62-53-17-31; Fax: +33-3-20-43-65-55

**Abstract:** ATP-dependent chromatin remodeling is the active displacement of nucleosomes along or off DNA induced by chromatin remodeling complexes. This key process of gene regulation in eukaryote organisms has recently been argued to be controlled by a kinetic proofreading mechanism. In this paper we present a discussion of the current understanding of this process. We review the case of gene repression via heterochromatin formation by remodelers from the ISWI family and then discuss the activation of the IFN-β gene, where the displacement of the nucleosome is initiated by histone tail acetylations by the enzyme GCN5 which are required for the recruitment of SWI-SNF remodelers. We quantify the specificity of the acetylation step in the remodeling process by peptide docking simulations.

**Keywords:** chromatin; histone code; remodeler; bromodomain; nucleosome; histone

## 1. Introduction

A key problem in molecular biology is to understand how genes are switched "on" or "off". For prokaryotes (bacteria) and phages this question has been studied since the 1960's, culminating in the development of the operon model by Monod and Jacob [1]. A particularly well-studied model organism is the λ-phage, whose diverse regulatory mechanisms are beautifully explained in Mark Ptashne's book [2]. The main regulatory mechanism of gene control relies on the combinatorial and cooperative binding of transcription factors at regulatory sites which either block or favor the recruitment of RNA polymerase, the reader and translator of the genetic
code. Gene activation/inhibition is thus physically controlled by the energetics of binding of
the transcription factors, turning the theoretical description of transcription initiation into a
problem of statistical physics [3].

While the prokaryote mechanism of transcription factor binding remains valid in eukaryotes,
the presence of the chromatin fiber requires additional regulatory mechanisms, notably to con-
trol the positioning of nucleosomes along the chromatin fiber, relative to gene and regulatory
sequences on DNA. One may therefore talk of a step of ‘pre-initiation’ of transcription which
renders the fiber accessible or inaccessible to the molecular regulators that directly promote the
readout of the gene. A particular difference between eukaryotes and prokaryotes is the presence
of nucleosomes and the associated numerous post-translational modifications on the nucleosomal
histone tails (and cores) whose combinatorial presence has been linked to a potential ‘histone
code’ [4], and further, the existence of a dedicated machinery of ‘chromatin remodeling’ en-
zymes which actively, i.e. under ATP-consumption, displace nucleosomes on and also from the
chromatin fiber [5].

Recently, there has been an attempt to relate specific histone modifications and chromatin
remodeling in a mechanistic picture by invoking a kinetic proofreading scenario [6, 7, 8]. Other
kinetic scenarios have also appeared in the literature [9, 10, 11, 12, 13]. The advantage of the
explicit kinetic proofreading scenario is that it specifically couples a recognition step of the
chromatin remodeler, which is controlled by the free energy of binding or dissociation, with
a kinetic step, which is a non-equilbrium reaction and consumes ATP. This latter step, in any
kinetic proofreading scenario as introduced by Hopfield [14] and Ninio [15], is the one that confers
a high specificity to the reaction. In the case of mRNA translation, such a specificity is required
in order to guarantee the quality of the protein product. In dynamic intermediates, as it is the
case in chromatin remodeling, such a high degree of specificity is not needed. Nevertheless, as
will be argued here, the mechanism is essential in promoting gene activation or repression in
eukaryotes, as it allows a specificity increase by several orders of magnitude, a feat which cannot
be achieved by free energy-dependent processes alone.

Kinetic proofreading in chromatin remodeling has so far seen experimental support and more
detailed theoretical analysis in the case of a remodeler from the ISWI family, ACF [7, 8]. ACF
displaces nucleosomes to form arrays, and hence specifically acts to repress or ‘shut down’ genes.
The activation of the remodeler is controlled by the unmodified histone tails of the histones H4,
while for gene activation, e.g. by the RSC remodeler, specific histone tail modifications must be
present [16].

In this paper, we first review the argument underlying the kinetic proofreading scheme in
chromatin remodeling as it has been developed in [6, 7, 8]. In fact, the original proposal of ref.
[6] was concerned with the case of gene activation, while the first experimental case, ISWI/ACF,
to which the scenario has been applied concerns a case of gene repression. After briefly reviewing
this case, we address the original situation of gene activation and discuss, as a novel application of
the kinetic proofreading scenario, the case of the IFN-beta interferon gene which had earlier been
studied experimentally in great detail by D. Thanos and collaborators [17, 18, 19, 20, 21, 22].
## 2. The Kinetic Proofreading Scenario of Chromatin Remodeling

Kinetic proofreading as a biochemical process for error-correction was first proposed by Hopfield [14] and Ninio [15]. A particularly lucid introduction to this mechanism can be found in Uri Alon's book [23]. Here, we consider directly the scenario as relevant to the regulation of nucleosome position due to chromatin remodeling. In this context the scenario presupposes only two experimental facts. The first is that histone amino acids can bear different chemical groups (post-translational modifications) which can be read by chromatin remodeler domains. 'Reading' the histone tail state refers to the binding of a remodeler recognition domain with a rate $k'$, and an unbinding with a rate $k$. We take the rate of unbinding as specific: the rate of binding $k'$ is determined by the frequency of molecular collisions, while the stability of the binding between the remodeler recognition domain and the histone tail decides on the rate of unbinding $k$. The second fact is that the active engagement of the remodeler in displacing the nucleosome is irreversible and consumes ATP; irreversible thus means that in order to undo the induced motion, additional ATP consumption is required.

If we denote the remodeler by $R$, the nucleosome by $N$ and the remodeler-nucleosome complex by $I$, the activated complex by $I^*$ and the mobile complex by $M$, the proofreading scenario is given by the following reaction scheme in which $m'$ is the ATP-dependent activation rate, $t$ the translocation rate of the nucleosome remodeler complex, and $\ell'$ the dissociation rate of the activated remodeler complex from the nucleosome [8]

$$R + N \underset{k}{\overset{k'}{\rightleftharpoons}} I \rightarrow_{m'} I^* \rightarrow_t M \\ \downarrow \ell' \\ R + N \qquad (1)$$

The reaction scheme is also depicted schematically in Figure 1 where we have omitted the translocation reaction $t$ which occurs repeatedly if the remodeler acts in a processive manner. We note that, of course, the rate scheme we adopt is a very simplified one. This is particularly relevant for e.g. in the case of ISWI, as the remodeler carries both the recognition domain and an accessory domain which binds with the DNA linker; this last feature is not retained in the model. Whenever a back-reaction is missing in the reaction scheme it is assumed as irreversible. The reaction scheme can easily be rewritten in the form of the rate equations for concentrations indicated by brackets $[\cdot]$,

$$\frac{d[I]}{dt} = k'[R][N] - (m' + k)[I], \qquad (2)$$

$$\frac{d[I]^*}{dt} = m'[I] - (\ell' + t)[I^*]. \qquad (3)$$

From these equations the ratio

$$f \equiv \frac{[I^*]}{[R][N]} \qquad (4)$$

can be calculated when stationarity of the reaction is assumed, in this way characterizing its discriminatory capacity. One finds

$$f \equiv \frac{m'k'}{(\ell + t)(m' + k)}. \qquad (5)$$
<figure><img src="image_3.png" /><figcaption>Figure 1. Schematic representation of the remodeling reaction under proof-reading. The red circle is the nucleosome, the histone tail (green) emerges from it and carries a modification (yellow triangle). The remodeler is drawn as a blue ellipse.</figcaption></figure>

The expression for the error-rate, the key result in kinetic proofreading, follows from the comparison of two reactions ('correct' and 'incorrect' or 'favored' and 'disfavored') described by ratios $f_1$ and $f_2$, which spell out as

$$F = \frac{m'_1 k'_1 (k_2 + m'_2) (\ell'_2 + p_2)}{m'_2 k'_2 (k_1 + m'_1) (\ell'_1 + p_1)} \quad (6)$$

In the original scenario, the presence of a transcription factor was also included in the reaction scheme [6]. A rough estimate for the error factor $F$ yielded a value of ≈ 400, which, as we will see below, is very close to the value found from experiments by G. Narlikar and collaborators [27, 28].

### 3. Gene Repression: Remodelers from the ISWI Family

The ISWI family of chromatin remodelers is one of the major known groups of remodeling enzymes in eukaryotes [24]. Figure 2 displays a schematic drawing of the different domains of the remodeling complex. The ISWI ATPase contains two major domains one of which is the motor domain. In addition to consuming ATP it is also sensitive to the interaction with the unmarked histone H4 tail, and to linker DNA on both sides of the nucleosomes, however not to DNA sequence. Two accessory motifs, termed AutoN and NegC, regulate this interaction in a competitive way [25]. AutoN contains a basic patch whose amino acid composition mirrors that of the unmarked histone tail. The presence of the histone tail is therefore required to unblock
<figure><img src="image_4.png" /><figcaption>Figure 2. Schematic representation of the domains of the ISWI remodeler. The ATPase is the motor unit which interacts with the unmodified histone tail H4. The neighbouring regions AutoN and NegC modulate this interaction. Hand-Sant-Slide denotes domains involved in regulating the interaction of the remodeler with DNA.</figcaption></figure>

the functioning of the ATPase. The second major domain, Hand-Sant-Slide (HSS) confers the recognition of the DNA linker and thus steers the action of the activated complex.

G.J. Narlikar and collaborators have studied the regulation of a family member of ISWI remodelers, ACF, which acts as a dimer on a nucleosome [26, 27, 28]. This allowed to determine several of the parameters in the Hopfield formula. In the notation used by Narlikar we can rewrite the equation as

$$F \equiv \frac{k_{I,1}k'_{1}(k_{2} + k_{I,2})(k_{off,2} + k_{tr,2})}{k_{I,2}k'_{2}(k_{1} + k_{I,1})(k_{off,1} + k_{tr,1})} \quad (7)$$

where the following parameters were determined by her work: $k_{I,1} = 20$/min; $k_{I,2} = 1$/min; $k_{off,1} = 8$/min, $k_{off,2} = 160$/min; $k_{tr,1} = 80$/min, $k_{tr,2} = 80$/min. For the remaining parameters we assume $k'_{1} = k'_{2}$ and $k_{i} = l'_{i} = k_{off,i}$. When these numbers are used, an error factor of the correct substrate of $F \approx 313$ is obtained. This estimate is indeed of the same order of magnitude as the one put forward in [6]. More recently, we have extended the ISWI/ACF scenario to include the AutoN and NegC motifs into the kinetic proofreading scheme [29]. As it is common in kinetic proofreading scenarii, the addition self-regulatory mechanism of ISWI helps to increase the specificity of kinetic proofreading. A precise quantitative estimate of this effect is currently not available.

## 4. Gene Activation: the IFN-β Gene

We now return to the case of gene activation for which the proofreading scenario was first proposed in [6] and in which the presence of transcription factors had been included. In fact, as we will see below for the IFN-β gene, the reality of gene activation is more complex. Figure 3 shows a schematic representation of the build-up of the IFN-β gene and its regulatory regions. The first key element is an enhancer region which is positioned in a nucleosome free region (NFR) in between two positioned nucleosomes (see Figure 3 a). The downstream nucleosome partially occludes the TATA-box transcription factor binding site upstream of the gene. In order to pre-initiate transcription of this gene, the downstream nucleosome must be remodeled in order
<figure><img src="image_5.png" /><figcaption>Figure 3. The regulatory region of the IFN-β gene. a) An enhancer element is located upstream of the gene in a nucleosome-free region (NFR). The TATA-box of the gene is partially occluded by a nucleosome. b) After remodeling, the nucleosome is moved downstream, allowing the general transcription factor TFIID to bind. Drawn after ref. [17, 19]; numbers in both graphs denote the position relative to the IFN-β transcription start site on the DNA sequence *before* remodeling, as the nucleosome has no precise position afterwards.</figcaption></figure>

to shift it further downstream to allow for the generic transcription factor TFIID to bind.

How can the IFN-β gene go from the inactive state shown in Figure 3 a), to the transcription initiation-ready state shown in Figure 3 b)? The first step is the assembly of the enhanceosome which consists of a transcription factor complex which contains a set of three different transcription factors engaged in triggering the inflammation response following a viral infection, among them the well-known NF-κB factor [17]. The enhanceosome then recruits a histone-tail modification writer, GCN5, which modifies the histone tails of the nucleosome. On histone H4, the lysine residue K8 is acetylated while on H3 K9 and K14 are acetylated [19]. Acetylation of H4-K8 allows the recruitment of a remodeling complex, SWI-SNF, which senses the modification via its bromodomains. The remodeler then shifts the nucleosome downstream. Acetylation of H3-K4/K9 recruits the general transcription factor TFIID and hence the IFN-β gene is ready for transcription. While Thanos *et al.* established the time series of these events, they did not quantify the process by their corresponding free energies of binding or dissociation. At present it is therefore not useful to cast the process into rate equations. Qualitatively, we distinguish between three components:

i) the build-up of the enhanceosome is the analogue to transcription factor recruitment in prokaryotes [3] in which, as in prokaryotes, DNA looping also plays a relevant role. For the IFN-β gene, the free energies controlling the build-up of the transcription factor complex could
<figure><img src="image_6.png" /><figcaption>Figure 4. Ribbon-style representation of bromodomains. α-helices are shown as flat ribbons, loops as thin wires. (A) Yeast GCN5-BrD; (B) structure-based alignment of yeast (purple) and human (green) GCN5-BrD with an RMSD of 0.8; small differences in the ZA- and BC-loops are observed; (C) human GCN5.</figcaption></figure>

in principle be studied in a similar manner, at least *in vitro*. This is nevertheless a complex problem as the binding of an extended transcription-factor complex will certainly modify the elastic properties of the chromatin fiber and hence affects the overall state of the fiber [30, 31].

ii) the recruitment of the remodeler via the histone-tail state and its subsequent action can be measured experimentally, as was shown in the case of ISWI/ACF [7]. For the remodeler RSC, a remodeler related to SWI/SNF, *in-vitro* measurements have been performed earlier which allows to estimate the specificity of the binding to the tail and the activation effect the remodeler has on the nucleosome [16].

iii) The remaining key intermediate step is the recruitment of the histone tail writer GCN5, which occurs via protein-protein interactions with the enhanceosome, and its action on the histone tails. In the following we will focus on this last aspect in order to understand the specificity inherent to this step.

**5. Free Energies of Histone Tail Peptides in GCN5 Bromodomains**

The recognition of (acetylated) histone tails is brought about by bromodomains which provide key protein interfaces in gene expression mechanisms in eukaryotes [32, 33]. GCN5 bromodomains have been studied in yeast and human [34, 35]. The GCN5 protein is ~439 amino acids long and contains a C-terminal bromodomain (BrD) which is composed of ~115 amino acids and features the general characteristics of BrDs, i.e. a left-handed four-helical bundle (Z, A, B and C helices), and two ZA- and BC-loop regions with two (sometimes also three) small single-turn helices in the ZA-loop region. The active site of the GCN5-BrD is comparatively narrow. Figure 4 shows a comparison of the GCN5 bromodomain structures for the two species (crystallographic data: PDB entries 1E61 for yGCN5 and 1F68 for hGCN5). In order to quantify the specificity of the histone tail modifications we have compared different bromodomains and modified histone
tails by combining molecular dynamics (MD), peptide docking and umbrella sampling [36]. In the first step we equilibrate a solvated bromodomain with a bound histone tail peptide carrying a specific acetylation state on the lysine residues, with the binding characterized by the binding free energy. We then pull out the peptide of the bound complex and measure the dissociation free energy from steered MD and umbrella sampling. The expectation is that the dissociation free energy of the preferred modification must be lower than that of an unfavorably modified tail, as our selected docking process mimics the release of the enzyme *after the modification has already been placed*, and not its binding to a still unmodified tail.

For the docking calculations of the peptides in the bromodomains we built histone tail peptides of 15 amino acids with a central acetyl-marked lysine flanked by 7 amino acids on each side in the form (7aa - acK - 7aa). Five different acetylation marks taken at different positions along each of the histone tails H3 (lysine K4, K9, K14, K18 and K27) and H4 (lysine K5, K8, K12, K16 and K20) were selected. The corresponding 15aa peptides were generated using the Tinker modeling package [37]. The lysine residues were mutated into acetyl-lysines with CHIMERA [38], minimized using the steepest-descent algorithm and the Gromacs 43a2 force field [39].

In the first MD-step we took each of the BrD models and removed water, ions and any other ligands. The BrD was then put back into a simple point charge (SPC) water box, then equilibrated and relaxed. Box volume was 7.8 nm³ and the minimum solute-solvent distance was 1 Å with the assumption of a normal charged system at pH 7. To counter the net charge of the BrDs with the biological salt concentration of 0.15 mM, the system was neutralized using 32 sodium and 34 chloride ions which were added according to the potential gradient of the simulated system. Following minimization, equilibration was performed with the backbone of the BrDs partially restrained for 1 ns. Subsequently, unrestrained MD was initiated and the time-dependent evolution of trajectories recorded for 5 ns for further analysis. During this production run, bond lengths were constrained using LINCS [40] and SETTLE was used for water molecules [41]. The time step of the simulation was kept at 2 fs and the simulation was performed in the [NPT] ensemble, using Berendsen algorithms to impose constant P = 1 bar and T = 300 K. The van der Waals cutoff was set to 12 Å. The long-range electrostatic forces were treated with the Particle Mesh Ewald (PME) method. Each of the 20 MD-optimized peptides were docked to the respective yeast or human BrD's using the AutoDock4.2 [42]. The central part of the peptide, i.e. the acK plus the two neighboring amino acids to both sides were treated as a flexible part, while the terminal flanking amino acids (5 + 5) were treated as a rigid part of the peptide. For the BrD only the active site residues with 5 Å radius (corresponding to the complete ZA-loop and BC-loop regions) were treated as flexible, while the rest of the amino acids were treated as rigid in order to reduce computing time. We employed a knowledge-based docking approach using the Genetic Algorithm scoring function with a grid size of 1 nm³, and 25 million energy evaluations per grid. We assigned the maximum number of torsions allowed in AutoDock, 32 and distributed these around the central amino acid. We also used AutoDock/VINA to go beyond the limit on the number of torsions to include all possible torsions [43]. The ten best binding configurations for each of the BrD-acK pair were selected and subjected to MD, steered-MD [44, 45] and umbrella sampling simulations [46, 47] to obtain the potential of mean force (PMF) for in total 20 complexes (2 BrD × 5 H3 and 5 H4 tail peptides) from which the free energies can be obtained using the WHAM routine for each of the 20 complexes.
<table><caption>Table 1. Free energies of binding (FEB) and dissociation free energies (DFE) for the histone tails H4 and H3 to the bromodomain GCN5 for yeast and human. Free energies of binding compare both acetylated (ac) and non-acetylated (nac) 15-mers at the indicated residue, while DFE-values are given for the acetylated residue only.</caption><thead><tr><th>ht</th><th>acK</th><th>Sequence</th><th>FEB ac/nac (kCal/mol)</th><th>DFE (kCal/mol)</th></tr></thead><tbody><tr><td rowspan="5">yH4</td><td>H4K5</td><td>GGSSGRGKGGKGLGK</td><td>-6.6/-5.4</td><td>24.10</td></tr><tr><td>H4K8</td><td>SGRGKGGKGLGKGGA</td><td>-6.2/-6.1</td><td>16.37</td></tr><tr><td>H4K12</td><td>KGGKGLGKGGAKRHR</td><td>-8.3/-5.6</td><td>24.37</td></tr><tr><td>H4K16</td><td>GLGKGGAKRHRKILR</td><td>-8.2/-6.5</td><td>24.94</td></tr><tr><td>H4K20</td><td>GGAKRHRKILRDNIQ</td><td>-8.3/-5.4</td><td>27.04</td></tr><tr><td rowspan="5">hH4</td><td>H4K5</td><td>GGSSGRGKGGKGLGK</td><td>-6.2/-5.3</td><td>23.10</td></tr><tr><td>H4K8</td><td>SGRGKGGKGLGKGGA</td><td>-6.1/-6.7</td><td>15.80</td></tr><tr><td>H4K12</td><td>KGGKGLGKGGAKRH</td><td>-8.1/-5.4</td><td>25.89</td></tr><tr><td>H4K16</td><td>GLGKGGAKRHRKVLR</td><td>-7.9/-6.0</td><td>24.23</td></tr><tr><td>H4K20</td><td>GGAKRHRKVLRDNIQ</td><td>-8.3/-6.0</td><td>28.61</td></tr><tr><td rowspan="5">yH3</td><td>H3K4</td><td>STGGARTKQTARKST</td><td>-6.8/-5.1</td><td>24.46</td></tr><tr><td>H3K9</td><td>RTKQTARKSTGGKAP</td><td>-6.9/-6.6</td><td>17.53</td></tr><tr><td>H3K14</td><td>ARKSTGGKAPRKQLA</td><td>-8.4/-6.2</td><td>30.54</td></tr><tr><td>H3K18</td><td>TGGKAPRKQLASKAA</td><td>-6.6/-6.2</td><td>24.51</td></tr><tr><td>H3K27</td><td>LASKAARKSAPSTGG</td><td>-6.5/-6.2</td><td>24.21</td></tr><tr><td rowspan="5">hH3</td><td>H3K4</td><td>STGGARTKQTARKST</td><td>-6.3/-5.9</td><td>26.32</td></tr><tr><td>H3K9</td><td>RTKQTARKSTGGKAP</td><td>-6.5/-5.1</td><td>16.46</td></tr><tr><td>H3K14</td><td>ARKSTGGKAPRKQLA</td><td>-5.8/-5.8</td><td>29.80</td></tr><tr><td>H3K18</td><td>TGGKAPRKQLASKAA</td><td>-6.7/-6.2</td><td>22.05</td></tr><tr><td>H3K27</td><td>LASKAARKSAPSTGG</td><td>-6.5/-6.3</td><td>23.79</td></tr></tbody></table>

## 6. Results and Discussion

We report the results of the pulling experiments in both graphical and tabular form. Table 1 lists all data for human and yeast on H3 and H4 tails together with the amino acid sequence that has been employed in the simulations. In Figure 5 we show the dissociation free energies for the H3 and H4 tails and both GCN5 domains, with the residue positions arranged linearly at integer values along the abscissa, starting from the end of the tail. Values of the dissociation free energy do not vary significantly between species but display similar trends for position. For the H4 tail, a distinctly lower value of the acetylated residue K8 is clearly discernible. This result clearly shows that the positioning of the acetyl-mark on residue K8 is associated with a preferred release of the bromodomain from the tail. The case of H3 is distinct from the previous case in that we can see that for residue K9 the DFE is equally lowered (as for the case of K8 in H4) while the value for acetylation of K14 is in fact increased. At present the physical interpretation of these values is not as easy as for H4K8 since we have not studied the dissociation free energy in the case where both modifications have already been set. As far as the comparison with
<figure><img src="image_7.png" /><figcaption>Figure 5. Dissociation free energies (in kCal/mol) as determined from the pulling protocol, as taken from Table I. On the abscissa the modified histone amino acids are placed on arbitrary integer positions; the drawn line is a guide to the eye. a) H4; b) H3. The yellow (lighter) curve is the yeast data.</figcaption></figure>

experiment is concerned, the scenario we propose has so far not been tested experimentally. Experiments in the field so far typically test the enzymatic activity of the histone writer based on mass spectrometry approaches or other bulk measurements (see e.g. [48] and references therein), while our approach is in the spirit of a (simulational) force experiment, providing by design access to specific modifications on individual peptides. Simulations retaining molecular details of the complex formation between the histone writer and the tail are meanwhile also available [49].

## 7. Conclusions

To conclude, in this paper we have first briefly reviewed the basic idea of the kinetic proof-reading scenario of chromatin remodeling which was originally proposed with transcriptional activation in mind [6]. As it happens, G. Narlikar and collaborators proposed a related scenario for the repressive ISWI-type remodeler ACF, which was therefore the first case to which this idea could be applied quantitatively. After reviewing these results, we turn back to transcriptional activation and propose that the IFN-β gene as an interesting case for the application of the scenario as for this gene the sequence of regulatory events has been well-studied, however not in a quantitative fashion. We therefore have addressed one key step in the process which is the placement of the histone tail modifications, as step which was ignored in [6]. We determined dissociation free energies of the bromodomain of the histone code 'writer' GCN5 for both yeast and human with molecular pulling simulations, in the idealized case of employing representatively constructed tail peptides. Our results are indicative that a full quantitative characterization of the kinetic proofreading scheme in the case of gene activation is possible in principle, if not by experiment, then at least by simulation. As an example for the latter can serve the recent work by Teif et al. [30], in which the binding of HP1 to the chromatin fiber is considered, taking into account the folding state of the chromatin fiber.
Acknowledgement

R.P.S. gratefully thanks the President of the University Lille I for a collaborative, three-year
PhD grant.

Conflict of Interest

All authors declare that there are no conflicts of interest.

References

1. Jacob F, Monod J (1961) Genetic regulatory mechanisms in the synthesis of proteins. *J Mol Biol* 3: 318-356.

2. Ptashne M (2004) *A Genetic Switch: Phage Lambda Revisited* (3rd. ed.), Cold Spring Harbor Laboratory Press, USA.

3. Ackers GK, Johnson AD, Shea MA (1982) Quantitative model for gene regulation by lambda phage repressor. *Proc Natl Acad Sci (USA)* 79: 1129-1133.

4. Strahl BD, Allis CD (2000) The language of covalent histone modifications. *Nature* 403: 41-45.

5. Narlikar GJ, Sundaramoorthy R, Owen-Hughes T (2013) Mechanisms and Functions of ATP-Dependent Chromatin-Remodeling Enzymes. *Cell* 154: 440-503.

6. Blossey R, Schiessel H (2008) Kinetic proofreading of gene activation by chromatin remodeling. *HFSP Journal* 2: 167-170.

7. Narlikar GJ (2010) A proposal for kinetic proof reading by ISWI family chromatin remodeling motors. *Curr Opin Chem Biol* 14: 660-665.

8. Blossey R, Schiessel (2011) Kinetic Proofreading in Chromatin Remodeling: The Case of ISWI/ACF. *Biophys J* 101: L30-L32.

9. Herschlag D, Johnson FB (1993) Synergism in Transcriptional Activation: A Kinetic View. *Genes Development* 7: 173-179.

10. Cosgrove MS, Boeke JD, Wolberger C (2004) Regulated Nucleosome Mobility and the Histone Code. *Nature Struct Mol Biol* 11: 1037-1043.

11. Rippe K, Schrader A, Riede P, Strohner R, Lehmann E, Längst G (2007) DNA sequence-and conformation-directed positioning of nucleosomes by chromatin-remodeling complexes *Proc Natl Acad Sci (USA)* 104: 15635-15640.

12. Erdel F, Schubert T, Marth C, Längst G, Rippe K (2010) Human ISWI chromatin-remodeling complexes sample nucleosomes via transient binding reactions and become immobilized at active sites. *Proc Natl Acad Sci (USA)* 107: 19873-19878.

13. Michel D (2011) Hierarchical cooperativity mediated by chromatin remodeling; the model of the MMTV transcription regulation. *J Theoret Biol* 287: 74-81.
14. Hopfield JJ (1974) Kinetic Proofreading: A New Mechanism for Reducing Errors in Biosynthetic Processes Requiring High Specificity. *Proc Natl Acad Sci (USA)* 71: 4135-4139.

15. Ninio J (1975) Kinetic amplification of enzyme discrimination. *Biochimie* 57: 587-595.

16. Ferreira H, Flaus A, Owen-Hughes T (2007) Histone modifications influence the action of Snf2 family remodelling enzymes by different mechanisms. *J Mol Biol* 374: 563-579.

17. Agalioti T, Lomvardes S, Parekh B, Yie J, Maniatis T, Thanos D (2000) Ordered Recruitment of Chromatin Modifying and Basal Factors to the IFN-b promoter. *Cell* 103: 667-678.

18. Lomvardas S, Thanos D (2001) Nucleosome sliding via TBP DNA binding in vivo. *Cell* 106: 685-696.

19. Agalioti T, Chen G, Thanos D (2002) Deciphering the Transcriptional histone acetylation code for a human gene. *Cell* 111: 381-392.

20. Koutroubas G, Merika M, Thanos D (2008) Bypassing the requirements for epigenetic modifications in gene transcription by increasing enhancer strength. *Mol Cell Biol* 28: 926-938.

21. Nolis IK, McKay DJ, Mantouvalou E, Lomvardas S, Merika M, Thanos D (2009) Transcription factors mediate long range enhancer-promoter interactions. *Proc Natl Acad Sci (USA)* 106: 20222-20227.

22. Ford E, Thanos D (2010) The transcriptional code of human IFN-beta gene expression. *Biochim Biophys Acta* 1799: 328-336.

23. Alon U (2006) *An Introduction to Systems Biology: Design Principles of Biological Circuits*, Chapman & Hall/CRC, United Kingdom.

24. Clapier CR, Cairns BR (2009) The biology of chromatin remodeling complexes. *Annu Rev Biochem* 78: 273-304.

25. Clapier CR, Cairns BR (2012) Regulation of ISWI involves inhibitory modules antagonized by nucleosomal epitopes. *Nature* 492: 280-285.

26. Yang JG, Madrid TS, Sevastopoulos E, Narlikar GJ (2006) The chromatin-remodeling enzyme ACF is an ATP-dependent DNA length sensor that regulates nucleosome spacing. *Nat Struct Mol Biol* 13: 1078-1083.

27. Racki LR, Yang JG, Nariman N, Partensky PD, Acevedo A, Purcell TJ, Cooke R, Cheng Y, Narlikar GJ (2009) The chromatin remodeller ACF acts as a dimeric motor to space nucleosomes. *Nature* 462: 1016-1021.

28. Blosser TR, Yang JG, Stone MD, Narlikar GJ, Zhuang X (2009) Dynamics of nucleosome remodeling by individual ACF complexes. *Nature* 462: 1022-1028.

29. Brysbaert G, Lensink MF, Blossey R (2015) Regulatory motifs on ISWI chromatin remodelers: molecular mechanisms and kinetic proofreading. *J Phys Condens Matter* 27: 064108.

30. Teif VB, Kepper N, Yserentant K, Wiedemann G, Rippe K (2015) Affinity, stoichiometry and cooperativity of heterochromatin protein 1 (HP1) binding to nucleosomal arrays, *J Phys Condens Matter* 27: 064110.

31. Lesne A, Foray N, Cathala G, Forné T, Wong H, Victor JM (2015) Chromatin fiber allostery and the epigenetic code. *J Phys Condens Matter* 27: 064114.
32. Filippakopoulos P, Picaud S, Mangos M, Keates T, Lambert JP, Barsyte-Lovejoy D, Felletar I, Volkmer R, Mller S, Pawson T, Gingras AC, Arrowsmith CH, Knapp S. (2012) Histone recognition and large-scale structural analysis of the human bromodomain family. *Cell* 149: 214-231.

33. Josling G, Selvarajah SA, Petter M, Duffy MF (2012) The Role of Bromodomain Proteins in Regulating Gene Expression. *Genes* 3: 320-343.

34. Owen DJ, Ornaghi P, Yang JC, Lowe N, Evans PR, Ballario P, Neuhaus D, Filetici P, Travers AA (2000) The structural basis for the recognition of acetylated histone H4 by the bromodomain of histone acetyltransferase Gcn5p. *EMBO J* 19: 6141-6149.

35. Hudson BP, Martinez-Yamout MA, Dyson HJ, Wright PE (2000) Solution structure and acetyl-lysine binding activity of the GCN5 bromodomain. *J Mol Biol* 304: 355-370.

36. Singh RP (2013) Computer Modeling of the Application of Mechanical Forces to Biomolecules. Thesis University of Lille.

37. Ponder JW, Richards FM (1987) An efficient Newton-like method for molecular mechanics energy minimization of large molecules. *J Comput Chemistry* 8: 1016-1024.

38. Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, Ferrin TE (2004) UCSF Chimera—a visualization system for exploratory research and analysis. *J Comput Chem* 25: 1605-1612.

39. Van der Spoel D, Lindahl E, Hess B, Groenhof G, Mark AE, Berendsen HJ (2005) GRO-MACS: fast, flexible, and free. *J Comput Chem* 26: 1701-1718.

40. Hess HB, Berendsen HJC, Fraaije JGEM (1997) LINCS: A linear constraint solver for molecular simulations. *J Comp Chem* 18: 1463-1472.

41. Miyamoto S, Kollman PA (1992) Settle: An analytical version of the SHAKE and RATTLE algorithm for rigid water models. *J Comput Chem* 13: 952-962.

42. Morris GM, Huey R, Lindstrom W, Sanner MF, Belew RK, Goodsell DS, Olson AJ AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility. *J Comput Chem* 30: 2785-2761.

43. Trott O, Olson AJ (2010) AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. *J Comput Chem* 31: 455-461.

44. Isralewitz B, Baudry J, Gullingsrud J, Kosztin D, Schulten K (2001) Steered molecular dynamics investigations of protein function. *J Mol Graph Model* 19: 13-25.

45. Isralewitz B, Gao M, Schulten K (2001) Steered molecular dynamics and mechanical functions of proteins. *Curr Opin Struct Biol* 11: 224-230.

46. Torrie GM, Valleau JP (1977) Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling. *J Comp Phys* 23: 187-199.

47. Hummer G, Szabo A (2010) Free energy reconstruction from nonequilibrium single-molecule pulling experiments. *Proc Natl Acad Sci (USA)* 107: 21441-21446.
48. Cieniewicz AM, Moreland L, Ringel AE, Mackintosh SG, Raman A, Gilbert TM, Wolberger C, Tackett AJ, Taverna SD (2014) The bromodomain of Gcn5 regulates site specificity of lysine acetylation on histone H3. *Mol Cell Proteomics* 13: 2896-2910.

49. Jiang J, Lu J, Lu D, Liang Z, Li L, Ouyang S, Kong X, Jiang H, Shen B, Luo C (2012) Investigation of the Acetylation Mechanism by GCN5 Histone Acetyltransferase. *PloS One* 7: e36660.

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

AIMS Press

©2015, Raghvendra P. Singh *et al.*, licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)