Structure-based design of an aromatic helical
foldamer-protein interface

Lingfei Wang, Céline Douat, Johannes Sigl, Post Sai Reddy, Lucile Fischer, Béatrice Langlois d'Estaintot, Zhiwei Liu, Vojislava Pophristic, Yuwei Yang, Yingkai Zhang, et al.

► To cite this version:

Lingfei Wang, Céline Douat, Johannes Sigl, Post Sai Reddy, Lucile Fischer, et al. Structure-based design of an aromatic helical foldamer-protein interface. Chemical Science, 2025, 16 (27), pp.12385-12396. 10.1039/D5SC01826A. hal-05396621

HAL Id: hal-05396621

https://hal.science/hal-05396621v1

Submitted on 3 Dec 2025

**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
EDGE ARTICLE

View Article Online

View Journal | View Issue

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

Check for updates

Cite this: *Chem. Sci.*, 2025, 16, 12385

All publication charges for this article have been paid for by the Royal Society of Chemistry

# Structure-based design of an aromatic helical foldamer-protein interface†

Lingfei Wang,ª Céline Douat,ª ID a Johannes Sigl, ID a Post Sai Reddy,ab Lucile Fischer, b Béatrice Langlois d'Estaintot, b Zhiwei Liu, c Vojislava Pophristic, c Yuwei Yang, d Yingkai Zhang ID d and Ivan Huc ID *a

The starting point of this study is the solid state structure of a complex between human carbonic anhydrase II (HCAII) and a helically folded tetradecaamide aromatic foldamer with a nanomolar HCAII ligand appended at the N terminus of the helix. In this complex, the foldamer is achiral but its handedness is biased by diastereoselective interaction with the protein. Computational analysis of the HCAII surface and inspection of the initial solid state structure led to the suggestion of main chain and side chain modifications of the foldamer helix that would result in an extension of the foldamer protein interface as well as in absolute helix handedness control. Molecular dynamics simulations validated several of these suggested modifications as potentially resulting in favorable foldamer-protein contacts. Five new Fmoc-protected amino acid building blocks bearing new biogenic-like side chains were synthesized. Nine new tetradecaamide sequences with or without the appended HCAII ligand were synthesized on solid phase and purified by RP-HPLC. The solid state structures of four of these sequences in complex with HCAII were obtained and validated the main design principles: (i) side chains can be predictably introduced at precise positions of the foldamer surface to create new contacts with the protein; (ii) side chains modifications do not alter main chain behavior and can be implemented independent from each other; (iii) some main chain units derived from quinoline-, pyridine-, or benzene-based δ-amino acids are largely interchangeable without altering the overall helix curvature in the context of a complex with a protein. An assessment of the K_D values required the adaptation of an existing fluorescence competition assay and suggested that the side chain and main chain modifications introduced in the new sequences did not result in significant improvement of the affinity of the foldamers to HCA.

Received 7th March 2025

Accepted 30th May 2025

DOI: 10.1039/d5sc01826a

rsc.li/chemical-science

## Introduction

Aromatic oligoamides represent a large class of compounds that can be used to recognize proteins and nucleic acids and that may interfere with their functions in multiple ways.¹ They comprise natural products such as distamycin,² cytobactamids,³ and albicidin,<fnref n="4" /> drug molecules such as suramin that has been crystallized bound to numerous proteins,⁵ rod-like oligomers many of which have been developed as α-helix mimetics,⁶ and oligomers that adopt helically folded conformations.⁷⁻¹⁰ We have been interested in the latter because their relatively large size offers the possibility to cover a large

surface area of a protein target, which is relevant to protein-protein and protein-nucleic acid interactions, two types of interactions that are difficult to inhibit with small molecules.¹¹ Helical aromatic oligoamide foldamers (AOFs) and in particular those derived from 8-amino-2-quinolinecarboxylic acid (Fig. 1) also possess the advantage that their conformations are very stable in particular in protic solvents,¹² and that synthetic

<figure><img src="image_3.png" /><figcaption>Fig. 1 Chemical formula (left) of a hexadecaamide of 8-amino-2-quinolinecarboxylic acid bearing a protein ligand (Lig) at the N terminus and biogenic-like side chains (R groups) in position 4. Schematic representation (right) of the helical structure of the hexadecaamide with some R groups interacting with a protein surface to which the ligand is also bound.</figcaption></figure>

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

<footnote>ªDepartment Pharmazie, Ludwig-Maximilians-Universität München, Butenandstr. 5-13, 81377 München, Germany. E-mail: ivan.huc@cup.lmu.de</footnote>

<footnote>ᵇCNRS, Bordeaux INP, CBMN, UMR5248, Univ. Bordeaux, F-33600 Pessac, France</footnote>

<footnote>Department of Chemistry & Biochemistry, Rowan University, 201 Mullica Hill Road, 08028 Glassboro, New Jersey, USA</footnote>

<footnote n="4">⁴Department of Chemistry, New York University, 10003 New York, New York, USA</footnote>

<footnote>†Electronic supplementary information (ESI) available: Experimental protocols and characterisation of new compounds. See DOI: https://doi.org/10.1039/d5sc01826a</footnote>
<figure><img src="image_5.png" /></figure>

methods exist to introduce various biogenic-like side chains at their periphery.¹³

The potential of helical AOFs to interfere with protein function has been highlighted in the context of amyloid fibers⁷ and DNA-binding proteins.⁸ For the latter, AOFs that specifically mimic the shape of charge distribution of DNA have been developed. In contrast, methods are still missing to design *ab initio* a helical AOF protein binder that does not mimic an already known binding epitope, for example through the introduction of biogenic-like side chains complementary to the protein surface. To assess the potential of some helical AOFs to interact with a given protein surface, we introduced a tethering approach where a covalent or non-covalent linkage confines the AOF at the surface of the protein (Fig. 1).⁸,⁹ Tethering between two small molecules or between a small molecule and its protein target is a common approach in the context of drug research to compensate for initially weak binding. It lies at the heart of linker design in fragment-based approaches,¹⁴ including in the context of template-assisted strategies,¹⁵ and of covalent ligands.¹⁶ Tethering of AOFs to a protein target was used to detect foldamer-protein interactions upon observing a preferred handedness in an achiral oligomer.⁸,⁹ Taken alone, the achiral AOF exists as a racemic mixture of right-handed (P) and left-handed (M) enantiomeric conformers. If either helix (M or P) interacts better than the other with the protein surface, the resulting change of proportion leads to an induced circular dichroism (CD) signal. An AOF CD signal is easy to detect because AOFs absorb above 350 nm, in regions where proteins are transparent. For these studies, we have used human carbonic anhydrase II (HCAII) as a model system.⁸,⁹ This protein was initially selected because it is robust and commercially available as an easy-to-handle freeze-dried solid,

i.e. it can be supplied to a chemistry laboratory without requiring recombinant protein expression. Furthermore, commercial HCAII was found to readily crystallize and small nanomolar ligands for this enzyme that can conveniently be appended to foldamers exist. Thus, strong helix handedness induction was observed with several helical AOFs linked to the protein *via* a nanomolar ligand. Subsequently, solid state structures of such complexes were obtained that confirmed the preferred P helix handedness and informed about foldamer protein contacts.⁸,⁹

In the case or tetradecamide sequence 1, the solid state structure showed a large contact surface area between the foldamer and HCAII (Fig. 2).⁹ However, it has not yet been shown whether such a structure could serve to further design the foldamer to extend its contact with the protein surface. Furthermore, although the AOFs are known to be rigid, it remained to be demonstrated whether side chain and main chain modifications could be implemented without altering their overall structure, a task difficult to achieve with *e.g.* a peptide or an aliphatic peptidic foldamer.¹⁷ Here, we show that, with the help of computational tools, the solid state structure of the complex 1·HCAII can be used as a starting point to place side chains at defined positions in space to further elaborate the foldamer-protein interface. We validate that the foldamer structure remains independent of side chain variations and even some main chain variations. Although the changes implemented have not resulted in significant changes in the dissociation constant of the complexes, the results further validate the concept that AOF helices can serve as reliable scaffolds to display biogenic-like side chains at the surface of a protein. The results also show that the HCAII system had some shortcomings alongside the advantages mentioned above.

<figure><img src="image_6.png" /><figcaption>Fig. 2 (a)–(d) Show different views of the solid state structure of the complex between HCAII and AOF 1.⁹ Only the surface of the protein is shown. The foldamer is shown in stick representation. Pyridine units (P) are shown in green and the HCAII ligand is shown in gold. In (b), two residues of interest are shown in purple. In (c) three carbon atoms of three P residues are shown in purple balls. These carbons would belong to the additional benzenic ring when implementing a P→Q mutation at these positions. Hydrogen atoms have been omitted for clarity. (e) Structural formula of 1.</figcaption></figure>
<figure><img src="image_7.png" /></figure>

# Results and discussion

## Synthesis and general design principles

Oligoamide sequence 1 consists of eight $Q^{Xxx}$ δ-amino acid monomers presenting different biogenic-like side chains in position 4, and of six P residues (Fig. 2e and 3). P residues bring the same contribution to helix curvatures as $Q^{Xxx}$ but they are more flexible. Their initial role was to make helix handedness dynamics fast enough to be practically monitored, e.g. in the course of minutes to hours,¹⁸ for example when helix handedness bias takes place upon binding of 1 to HCAII.⁹ Unexpectedly, P residues were found to be directly involved in

foldamer-protein contacts in the solid structure of 1·HCAII (Fig. 2b). Their role thus extends to that of interacting units despite the fact that they carry no biogenic-like side chains.

In order to extend foldamer-HCAII contacts, a number of $Q^{Xxx}$ monomers were considered bearing side chains in position 4, 5 or 6 (Fig. 3). $B^{Xxx}$ monomers were also involved. $B^{Xxx}$ monomers are δ-amino acids as well and may thus bring a contribution to helix curvature similar to that of Q and P.¹⁹ In addition, they carry a stereogenic center that has been shown to quantitatively bias helix handedness.²⁰ Sequences comprising $B^{Xxx}$ monomers are thus designed to be one handed even when they are not bound to HCAII. The S configuration of $B^{Xxx}$

<figure><img src="image_8.png" /><figcaption>**Fig. 3** Structural formulas of P, $Q^{Xxx}$, and $B^{Xxx}$ δ-amino acids and of N-terminal Lig and Deg functional groups. Carbon atoms in position 4, 5, and 6 of $Q^{Xxx}$ are indicated and carry R, $R^1$, and $R^2$ side chains, respectively. The Xxx three letter code used for the side chains is sometimes inspired by the three letter code of α-amino acids bearing similar side chains, even when they may not exactly match. The three letter code also indicates when the side chain is in position 5 or 6 of the quinoline ring. Sequences 1–21 are defined with the letter code used in this study. To facilitate residue identification, the colors of the highlighted residues match with the side chain colors in the adjacent boxes. Sequence numbers with a star (4\*-14\*) indicate sequences that were investigated in molecular dynamics simulations only and were not synthesized.</figcaption></figure>
<figure><img src="image_9.png" /></figure>

monomers is intended to favor *P* helicity. Among the monomers used in this study (Fig. 3), some were previously described, some were used only in computations, and five – *B*<sup>Ida</sup>, *B*<sup>Gpr</sup>, *Q*<sup>Gly</sup>, *Q*<sup>5Ph</sup> and *Q*<sup>6gp</sup> – were newly synthesized (Fig. S1†). The preparation of these building blocks in a form suitable for solid phase synthesis is presented in detail in the ESI.† All monomers were produced with a free carboxylic acid and an Fmoc-protected main chain amine. In addition, the side chains of *B*<sup>Gpr</sup> and *Q*<sup>6gp</sup> were protected with Boc groups. Typically, side chain installation involved Sonogashira or Suzuki cross-coupling reactions on a bromoaryl precursor.<sup>13</sup> The synthesis of *Q*<sup>5In</sup> is also described in the ESI.† However, although this monomer is a stable compound in its protected form (with a Boc group on the indole ring), sequences that contain it degrade over time after deprotection, apparently through an oxidation process, and could not be investigated. A representative example of this degradation is shown in Fig. S2.† It appears that oxidation is mediated by the foldamer and possibly involves light. For instance, we found that biotinylated sequences must be protected from air and light for storage to avoid the formation of biotin oxide.<sup>10a</sup>

Sequences **1–3** and **15–21** (Fig. 3) were synthesized on solid phase using an established *in situ* acid chloride activation protocol for the coupling steps.<sup>21</sup> An improved procedure for the on-resin introduction of the HCAII ligand at the N-terminus of the helix using a urea linkage was also developed. Sequences were purified by RP-HPLC after TFA-mediated resin cleavage and side chain deprotection.

Sequences **4**\*-\* **14**\* were investigated in computational studies but not synthesized. Sequences **2–3**, **4**\* and **5**\* were the focus of a first phase of our investigation. They are analogues of **1** in which P10 is replaced by different, more rigid Q10, residues. In the case of **3**, it could be verified experimentally that helix handedness bias upon binding to HCAII takes place despite the added rigidity, albeit significantly slower than with **1** (Fig. S3†). This first phase led to the installation of a *Q*<sup>5Ph</sup>10 residue in sequence **3**, instead of P10 in **1**. *Q*<sup>5Ph</sup>10 was conserved in all subsequently synthesized sequences except in **15**. In a second phase, residue variations in positions 6, 7, 11, 12 and 13 were assessed computationally in sequences **6**\*-\* **14**\* (see next section) and a selection of these variations was experimentally implemented in **15–20**. Sequence **21** is an analogue of **3** lacking the N-terminal ligand.

Including chiral *B*<sup>XXx</sup> units to favor *P* helix handedness was desirable, for example to avoid conformational changes in the course of a *K*<sub>D</sub> value determination. For this purpose, residues 6 and 11 were chosen as possible locations. This choice was based on the observation that the side chains in position 4 of the quinoline rings of *Q*<sup>ASP6</sup> and *Q*<sup>Ace11</sup> of **1** lie far from the HCAII surface (Fig. 2b). Removing these side chains and part of the pyridine ring of *Q* when performing a *Q*→*B* mutation should not alter the HCAII–foldamer interactions observed in the solid state. In contrast, the carbon atoms in position 5 and 6 of the quinoline rings of *Q*<sup>ASP6</sup> and *Q*<sup>Ace11</sup> in **1** seem better oriented to introduce a side chain that may interact with the HCAII surface and a *B*<sup>XXx</sup> monomer may offer a similar side chain presentation.

Before implementing *Q*→*B* mutations, another design feature had to be considered. While chiral *B* units have been shown to quantitatively bias helix handedness in the context of (*Q*)<sub>n</sub> oligomers,<sup>20</sup> this has not been validated when the helix also contains multiple more flexible *P* units, as in **1**. Indeed, partial handedness bias has occasionally been observed when *Q* monomers are mixed with other monomers.<sup>22</sup> To mitigate the risk that the chiral *B*-containing sequences would not be quantitatively one handed, we replaced some *P* units by *Q* in the vicinity of the *B*6 and *B*11 monomers. With their additional fused benzenic ring, *Q* monomers are bulkier than *P*. The structure of the **1**·HCAII complex showed that this extra bulk could be accommodated without generating clashes in P7, P10 and P12, but not in P5 and P8 (Fig. 2c). Chiral *B*-containing sequences **15–20** therefore contain at least one and sometimes up to three *Q* monomers at residues 7, 10 or 12.

The one-handed nature of the new chiral foldamers could be verified by <sup>1</sup>H NMR spectroscopy through the observation of a single set of signals. On top of ensuring quantitative handedness bias, the additional *Q* residues also resulted in slow helix handedness inversion in water. In case handedness bias was incomplete when the foldamer was first dissolved in water, e.g. for RP-HPLC purification, it may no longer proceed to completion. This pitfall is easily detected by the observation of two distinct sets of signals on the <sup>1</sup>H NMR spectra, corresponding to *P* and *M* diastereomeric conformers. To solve this problem, one can dissolve and incubate the compound in an organic solvent such as DMF, where helix handedness inversion takes place faster,<sup>12</sup> before evaporating and redissolving in water.

## Computational design

**Protein surface analysis.** The potential of the HCAII surface for interacting with biogenic-like residues was assessed with AlphaSpace, a computational analysis tool designed for fragment-centric topographical mapping.<sup>23</sup> The assessment proceeded in two phases. In a first phase, the surface in the vicinity of the HCAII active site was analysed, leading to the identification of potential binding pockets Po1–Po5 (Fig. 4b). Po1 has the highest ligandability (highest Bscore) and corresponds to the HCAII active site where HCAII ligands usually bind.<sup>24</sup> In the **1**·HCAII complex, Po1, Po2, Po3 and Po5 are filled by the N-terminal ligand and the helix backbone, as indicated by the color patches in Fig. 4a, leaving essentially no space to add functionalities on the foldamer helix to further enhance contacts with the protein surface. In contrast, Po4 was identified as a sizeable (<sup>158</sup>Å<sup>3</sup>) cavity nearby P10. Since a P10Q mutation appeared to be feasible without causing steric clashes (Fig. 2c), various side chains were docked in Po4 while being connected to the C5 carbon of the quinoline ring of Q10. All 274 side chains of the Swiss amino acid database<sup>25</sup> were tested. In each case, the amino acid was replaced by the quinoline residue and Autodock Vina<sup>26</sup> was used to score interactions between the side chain in position 5 and Po4 (Fig. S4†). *Q*<sup>5Ph</sup>, *Q*<sup>5Pa</sup>, and *Q*<sup>5In</sup> were selected as having a sufficiently low estimated Δ*G* and as being at the same time synthetically accessible. As presented in detail below, subsequent computational steps, synthesis and
<figure><img src="image_10.png" /></figure>

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

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

<table><thead><tr><th>Pocket Label</th><th>Pocket Color</th><th>Space (Å<sup>3</sup>)</th><th>Bscore (kcal/mol)</th><th>Nearby residues</th></tr></thead><tbody><tr><td>Po1</td><td>Green</td><td>379</td><td>-7.5</td><td>Lig</td></tr><tr><td>Po2</td><td>Yellow</td><td>153</td><td>-1.6</td><td>Lig</td></tr><tr><td>Po3</td><td>Teal</td><td>64</td><td>-2.6</td><td>5, 7, 8</td></tr><tr><td>Po4</td><td>Pink</td><td>158</td><td>-3.5</td><td>10</td></tr><tr><td>Po5</td><td>Orange</td><td>35</td><td>-1.1</td><td>3</td></tr></tbody></table>

<table><caption>Fig. 4 (a) Structural formula of 3 including Q and P residue numbering. The bonds shown in red highlight parts of the molecule in direct contact with the HCAII surface in the solid state structure of the 3·HCAII complex (Fig. 5). The bonds shown in green highlight parts of the foldamer involved in intercomplex contacts in the crystal lattice of the solid state structure of the 3·HCAII complex. The color patches indicate the pockets near the foldamer main chain or side chains in the solid state structure of the complex. Pockets are colored and numbered as in (b) and (c). (b) Pocket analysis of the surface of HCAII restricted to the vicinity of the ligand binding site and the contact area with the foldamer helix in the solid state structure of the 1·HCAII complex (Fig. 2). (c) Pocket analysis of the surface of HCAII restricted to the vicinity of the contact area with the foldamer helix in the solid structure of the 3·HCAII complex (Fig. 5). In (b) and (c), pockets have been assessed in terms of their volume and their ligandability (BScore). Pockets Po10-Po14 were not included in this study.</caption><thead><tr><th>Pocket Label</th><th>Pocket Color</th><th>Space (Å<sup>3</sup>)</th><th>Bscore (kcal/mol)</th><th>Nearby Residues</th></tr></thead><tbody><tr><td>Po6</td><td>Rosy</td><td>254</td><td>-4.6</td><td>6</td></tr><tr><td>Po7</td><td>Aqua</td><td>127</td><td>-2.8</td><td>11</td></tr><tr><td>Po8</td><td>Blue</td><td>57</td><td>-1.0</td><td>11, 13</td></tr><tr><td>Po9</td><td>Light pink</td><td>83</td><td>-2.2</td><td>10</td></tr></tbody></table>

structural analysis eventually delivered the solid state structure of the 3·HCAII complex where Po4 is indeed filled by the phenethyl side chains of Q⁵Ph¹⁰.

The second surface analysis was performed on the 3·HCAII complex in order to identify potentially ligandable sites in the vicinity of the foldamer helix where foldamer-protein contacts may be extended through the addition of foldamer side chains (Fig. 4c). This analysis led to the identification of pockets Po6-Po9. Po9 consists of the space left in Po4 that is not occupied by the side chain of Q⁵Phe¹⁰, hinting at the possibility that this side chain may be further elaborated (Fig. S5†). However, this option was not explored as many of the suggested side chains were synthetically challenging. We focused instead on Po6-Po8 which all lie on the same side of the foldamer helix, and may potentially be reached with side chains on residues Q6, P8, Q11, and Q13 (Fig. 4a). As explained above for Q10, side chains that both had a reasonable docking score and appeared to be

synthetically accessible were kept for subsequent investigations (Fig. S6 and S7†). For Q6 and Q11, we have mentioned above that side chains in position 4 of the quinoline do not establish contacts with the protein surface (Fig. 2b) and that these positions were considered for the introduction of chiral B residues to control helix handedness. Instead, the HCAII surface analysis suggested the side chains in position 6 of the quinoline ring might establish contacts with the protein. This eventually led to the mutation of Q⁵Phe⁶ into guanidinium-containing B⁵Phe⁶ or B⁵Phe⁶ in sequences 9*, 10*, 13*, 14*, 18, 19, as well as indane-containing residue B¹⁰Phe⁶ in sequences 8*, 12*, 20. Similarly, mutation of Q⁵Phe⁶ to guanidinium-containing residues Q⁵Phe⁶, Q⁶Phe⁶, and B⁵Phe⁶ was implemented in sequences 7*-14* and 16-19.

Possible modifications of Q⁵Phe⁶ were inspired by a salt bridge between this residue and Lys24 of HCAII observed in one of the solid state structures presented below. To better reach

<figure><img src="image_13.png" /></figure>
<figure><img src="image_14.png" /></figure>

this Lys24, residues Q<sup>But</sup>13, Q<sup>Bph</sup>13 and Q<sup>5Bu</sup>13 were considered in sequences 7*-9* and 11*-13*. It should be pointed that, given the extensive HCAII surface that the foldamer helix covers, opportunities for mutations and for the creation of new foldamer-protein contacts were too numerous to be considered at the same time. For instance, pockets Po10-Po14 were not investigated (Fig. 4c).

**Molecular dynamics (MD) simulations.** Prior to investing time and resources in the preparation of new momomers and new sequences suggested by the HCAII surface analysis, the effect of foldamer modifications in sequences 3 and 4*-14* on interactions with HCAII were evaluated using MD simulations in explicit water using the AMBER22 package.<sup>27,28</sup> The initial HCAII and foldamer structures and positions were based on structure alignments to the 3·HCAII solid state structure. The ff14SB force field<sup>29</sup> was used for α-amino acid residues. The general AMBER force field (GAFF),<sup>30</sup> with improved torsional parameters for arylamides;<sup>31</sup> was used for the foldamer (see ESI† for details). One additional simulation in the presence of 125 mM NaCl was performed on 3·HCAII. It resulted in minor changes such as slightly larger fluctuations, deviations, and a reduction of the occurrence of salt bridges. To inspect the interactions between foldamer and HCAII, we carried out a combination of structure visualization, calculations of root mean square displacements (RMSD) of protein and foldamer backbone atoms with respect to the solid state structure (Fig. S8†), as well as analysis of specific residue-to-residue distances. With the exceptions of sequences 7* and 11*,‡ all backbone RMSDs stayed within 3 Å of the solid state structure. Furthermore, sequence pairs with similar sets of mutations (7* and 11*, 8* and 12*, 9* and 13*, 10* and 14*) exhibit similar RMSDs that differ from those of other pairs, with these differences beginning to emerge after approximately 100 ns. Partial dissociation of sequence 7* was observed after 200 ns (Fig. S8†). The observed RMSDs, along with subsequent residue-to-residue distance analyses, support that conformational sampling from the 500 ns trajectory is sufficient for examining the interactions between HCAII and foldamers in the vicinity of the ligand binding site.

A first aspect concerns the rationale that led to selecting sequences 3 and 4*-14*. Because the number of possible single mutations was large, these were not investigated individually. Most sequences, carry two, three or four simultaneous side chain modifications with respect to 1, at positions 6, 10, 11 and 13. This way, all but two side chain modifications were examined in at least two distinct MD simulations. The first essential result is that the various side chain and, sometimes, main chain modifications mainly depend on where they are implemented, and generally do not depend from one another. When a mutation is performed at a given position, the behavior of the new residue tends not to vary with mutations at other positions. This is a major advantage for making predictions and sharply contrasts with aliphatic peptides where local modifications may impact global behavior.<sup>17</sup> The consistent behavior of each new residue regardless of other sequence modifications also suggests that no mutation led to a major steric clash that would disturb the whole structure. In addition, owing to the

independent behavior of the side chains, we could perform an analysis per interaction site/pocket, instead of an analysis per sequence. The results are presented in Fig. S9-S14† and a representative example is shown in Fig. 5.

The MD simulations reflected the strong interaction between the ligand and HCAII (Fig. S9-S10†). In all simulations, the bond between the HCAII-bound Zn<sup>2+</sup> ion and the ligand sulfonamide group, tight contacts between the two aryl groups of the ligand and pockets Po1 and Po2, and contacts between Q3 and P5 with HCAII in Po2, all remained well in place (Fig. 2a, 4a and b). Concerning the mutations of P10 implemented to fill Po4, Q<sup>5Ph</sup>10 was found to form stable hydrophobic contacts with Pro137, Leu203, Glu204 and Cys205, as highlighted by the histograms of distance shown in Fig. 5. Note that sequences 4* and 5* show some deviations in these histograms because their different Q<sup>5In</sup>10 and Q<sup>5Pa</sup>10 residues establish distinct contacts. The indole side chain of Q<sup>5In</sup>10 in sequence 4* lies closer to Pro137 (Fig. S12†). This residue was synthesized but sequences containing it had stability issues that hampered experimental investigations. The benzamidinium side chain of Q<sup>5Pa</sup>10 in sequence 5* appeared to be too large for pocket Po4 and its position fluctuated more (Fig. S12†). It was not considered further and Q<sup>5Ph</sup>10 was conserved in all subsequent experiments.

The interactions between side chains on residues 6, 11 and 13 and pockets Po6-Po8 (Fig. S13†) can be summarized as follows. Overall, the distance histograms show larger variations than for the contacts in pockets Po1-Po4. The salt bridge between Q<sup>Asp</sup>13 and Lys24 seen in one solid state structure (see below) was absent or present in small percentage (3% to 41%) of the time along the trajectories in aqueous solution. Using other negatively charged side residues Q<sup>But</sup>13, Q<sup>Bph</sup>13, and Q<sup>5Bu</sup>13 made little difference. These were therefore not tested experimentally. In the case of residue 6, the hydrophobic side chain of B<sup>Ida</sup>6 or the cationic side chains of B<sup>Gpe</sup>6 or B<sup>Gpr</sup>6 could potentially form contacts with Phe20 and Asp19, respectively, within pocket 6. Some of these residues were subsequently synthesized and implemented in sequences 19–20. Finally, the benefit of cationic residues in position 11 to fill pocket Po8 was not clear. Salt bridges were established only during small fractions of simulation time. With sequence 10*, a possible exception to the independent role of the side chains was observed with an apparent positive cooperative effect of the guanidinium-containing side chains of B<sup>Gpe</sup>6 and Q<sup>6gp</sup>11 (Fig. S13†). Specifically, salt bridging (59% of simulation time for Q<sup>6gp</sup>11-Asp19 and 95% for B<sup>Gpe</sup>6-Asp19) and hydrogen bonding (68% of simulation time for Q<sup>6gp</sup>11-Asp19) are more frequent in 10*·HCAII than with 7*, 8*, 9* and 14* in which either B<sup>Gpe</sup>6 or Q<sup>6gp</sup>11 is present (0 to 28% for Q<sup>6gp</sup>11-Asp19, 74% for B<sup>Gpe</sup>6-Asp19 and 0–11% H-bond for Q<sup>6gp</sup>11-Asp19).

MD simulations revealed an additional, unplanned, favorable foldamer-protein contact between Q<sup>Hyd</sup>7 and Gln135 (Fig. S14†). As mentioned above, P7Q<sup>Hyd</sup> and P12Q<sup>Hyd</sup> mutations were introduced to rigidify the helix and mitigate the risk that chiral B residue may not quantitatively bias helix handedness. Q<sup>Hyd</sup> residues were selected for that purpose because their relatively acidic yet small hydroxy side chain would not decrease
<figure><img src="image_15.png" /><figcaption>Fig. 5 (a) MD simulations of complexes between HCAII and foldamers 3 and 4*-14* showing the occupancy of pocket Po4 by the side chain of Q10, that is, Q⁵Ph¹⁰ for all sequences but 4* (Q⁵In¹⁰) and 5* (Q⁵Pa¹⁰). The histograms show the probability, through the entire simulation time, of the distance between the center of mass (COM) of the phenyl ring of Q⁵Ph¹⁰ or Q⁵Pa or of the pyrrole ring of Q⁵In, and HCAII residues Cys205 (position of Cβ), Glu204 (COM of Cβ and Cγ), Leu203 (COM of all side chain C atoms), and Pro137 (COM of Cβ, Cγ, and Cδ). Two distinct simulations were performed with sequence 3. The simulation marked with # included 125 mM NaCl. Arrows highlight the different positioning of Q⁵In¹⁰ and Q⁵Pa¹⁰ in the pocket. The histograms show weaker probabilities for sequences 7* and 11* due to strong deviations from the initial structure in these two cases (the Q10 residues are most of the time at distances >12 Å from Po4). (b) Snapshot from the MD simulation of 7*. HCAII showing the relevant residues. In this snapshot, one can also spot a transient salt bridge between the carboxylate side chain of Q⁵But¹³ and the Arg27 (dashed yellow lines). This salt bridge does not occur frequently.</figcaption></figure>

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

foldamer water solubility. Nevertheless, MD simulations suggest that the side chain of QHyd7 can also hydrogen bond to Gln135.

## Structure elucidation

Crystallization was attempted for all foldamer-HCAII complexes. In the case of 2·HCAII, 3·HCAII, 16·HCAII and 20·HCAII, single crystals suitable for X-ray diffraction analysis were obtained (Fig. S15†) and the solid state structures were elucidated in the P2₁2₁2₁₂ space group at a resolution of 1.4, 2.1, 1.6, and 2.1 Å, respectively (Fig. S16 and S17†). For the four structures, crystallization conditions were similar to that of 2·HCAII, and so were the unit cells and packing arrangements. Some parts of the foldamer molecule are involved in intercomplex contacts in the crystal lattice, including the side chain of QAcce11 (Fig. 4a and S18†), these contacts are all conserved in foldamers 1, 2, 3, 16 and 20. In retrospect, we hypothesized that the unsuccessful crystallization of the complexes with 17, 18 and 19 may be assigned to the mutation of QAcce11 into BGpr11 or Q⁶gP11 in these compounds.

The structure of 2·HCAII validated that a P10QGly mutation of the foldamer could be performed without any steric apparent clash between the protein and QGly10 or any alteration of the foldamer helix shape (Fig. 6a and b). This structure also revealed a salt bridge between the carboxylate of the QAsp13 and residue Lys24 (Fig. 7a). In the structure of 1·HCAII, the side chain of Lys24 was only partly visible in the electron density map and the salt bridge was overlooked. As mentioned above, MD simulations suggested that this salt bridge is not stable in aqueous solution and is not convincingly stabilized when using anionic side chains longer than in QAsp13, or placed in position

5 of the quinoline ring, or having a dianionic phosphonate group, as in QBut, QBph, and Q⁵Bu.

The structure of 3·HCAII then validated that the phenethyl side chain of Q⁵Ph¹⁰ filled Po4 as predicted by computations (Fig. 6c and d). The methylene carbon atom linked to the quinoline ring was found at 3.9 Å from a methyl group of Leu203, and one carbon atom of the phenyl ring lies within 3.5 Å from the nitrogen atom of Pro137. The structure of 16·HCAII confirmed the position of Q⁵Ph¹⁰ found in 3·HCAII and validated the double mutation QAsp6BGly and P7QHyd intended to introduce helix handedness bias (through BGly6) and to make the helix more rigid (through QHyd7). The proximity of QHyd7 and Gln135 (Fig. 6e) makes the hydrogen bonding observed in MD simulations plausible. In the solid state, the Gln135 amide NH₂ hydrogen bonds to the main chain carbonyl of QHyd7. In addition, the proximity between the primary amide of Gln135 and the hydroxy side chain of QHyd7 likely favors contacts with the latter as well, be it in a protonated or deprotonated state. Finally, the structure of 20·HCAII validated that the indane side chain of BIda6 filled Po6 again as predicted by computations, establishing contacts with Phe20 (Fig. 6f). In this structure, Lys24 was again visible in the electron density map, but in a conformation where hydrogen bonding to QAsp13 is not established (Fig. 7b), different from the structure of 2·HCAII. Of note, in all structures, QAsp13 is involved in a salt bridge with a lysine (Lys80) belonging to another HCAII molecule of the crystal lattice as part of the intercomplex contacts (Fig. S18†). This probably influences, that is, competes with the formation of the salt bridge with Lys24 in the solid state.

Altogether, the solid state structures validate the predictions made by computations. They demonstrate the equivalent contribution to main chain helix curvature of P, Q and B
<figure><img src="image_17.png" /></figure>

<figure><img src="image_18.png" /><figcaption>**Fig. 6** Solid state structures of HCAII in complex with: (a) 1, (b) 2, (c) 3, (d and e) 16, (f) 20. In (a)–(c), the protein is shown as a white soft surface and the foldamer is shown in stick representation in gray except the residue in position 10 colored in purple. In (d)–(f), the protein is shown in gray ribbon representation except relevant amino acids which are in space filling representation. The foldamer is shown in blue stick representation except the residue in position 10 in (d), the residue in position 7 in (e) and the residue in position 6 in (f), which are colored in purple. Hydrogen bonds are indicated as yellow dashed lines. Pockets Po4 and Po6 are defined in Fig. 4.</figcaption></figure>

<figure><img src="image_19.png" /><figcaption>**Fig. 7** Solid state structures of HCAII in complex with: (a) 2, (b) 20. The protein is shown in gray ribbon representation except P21 and Lys24, which are in space filling representation. The foldamer is shown in blue stick representation except the residue in position 13, which is colored in purple. Hydrogen bonds (here within a salt bridge) are indicated as yellow dashed lines.</figcaption></figure>

monomers in the context of a foldamer-protein contact area. They also demonstrate that the positions of the foldamer side chains and the type of interactions they may engage with the protein are predictable.

## Binding studies

We set out to measure the binding affinities of *P*-helical, chiral *B*-containing sequences **15–20** for HCAII to assess the extent to which they reflected the changes introduced in the foldamers. Sequences **1–3** are potentially problematic as they exist as

a racemic mixture of *M*- and *P*-helical conformers that must have different $K_D$ values and whose proportion evolve with time upon binding to HCAII, hence the focus on **15–20**. This assessment proved challenging. Simple ligands such as **22** (Fig. 8a) – the fragment of **15–20** that fills pockets Po1 and Po2 of HCAII – bind in the low nM range. Getting accurate $K_D$ values to comment on potentially small effects for such strong binding is delicate. Furthermore, we have shown that appending a foldamer on **22** has one major consequence: both the association and dissociation kinetics are slowed down by almost two orders of magnitude.⁸ᵃ Unlike with classical small molecule HCAII ligands, dissociation becomes so slow that techniques such as surface plasmon resonance (SPR) or biolayer interferometry (BLI) no longer deliver reliable results.

We turned to a recently published assay that exploits the quenching of the fluorescence of nanomolar ligand **23** upon binding to HCAII (Fig. 8a).³² The low $K_D$ value of **23** makes it suitable to perform competition (displacement) assays with compounds binding in the same concentration range. However, some optimization of the assay was required to perform experiments with foldamer-containing ligands because the foldamers absorb both at the excitation (373 nm) and emission (400–450 nm) wavelengths of **23**. Performing classical direct displacement titrations where a foldamer is added to a solution containing **23** and HCAII would be complicated by variable inner filter effects. Instead, we performed titrations in which aliquots of an HCAII solution, typically 10 µM, were added to a solution already containing a foldamer (50 nM) and **23** (50 nM). The foldamer and **23** were also present at the same concentrations in the HCAII solution. This way, the concentrations of fluorophore and foldamer were kept constant and only the ratio of HCAII was varied. A representative titration is
<figure><img src="image_20.png" /></figure>

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

<figure><img src="image_22.png" /><figcaption>**Fig. 8** (a) Structural formula of HCAII ligand **22** and of fluorescence probe **23** used in the competition assay. (b) Changes in fluorescence spectra (380–600 nm) upon HCAII titration of sequence **19**. [19] = [23] = 50 nM; (c) Experimental (■) and calculated values using a 1:1 binding isotherm (–) of fluorescence intensity of sequences **15–20** and **22** titrated with HCAII. [competing species] = [23] = 50 nM. Note that the curve fitting is shown at one wavelength (419 nm) as an illustration but that the $K_D$ values (Table 1) were calculated by simultaneously fitting data recorded in the 380–600 nm range. Blue arrows point to a systematic error discussed in the text.</figcaption></figure>

shown in Fig. 8b and the corresponding $K_D$ values are shown in Table 1. With this assay, the $K_D$ value for simple ligand **22** was 10 nM compared to 5 nM previously measured by SPR with HCAII immobilized on the SPR chip under slightly different buffer conditions.⁸ᵃucht For the titrations of the foldamers, one may note a slight but persistent systematic deviation in the curve fitting of the data by the 1:1 binding isotherm (blue arrows in Fig. 8). It appears that the fluorescence drops upon adding HCAII, i.e. that ligand **23** binds to HCAII, in presence of the foldamer faster than predicted by the binding model. A better fit can be artificially obtained by inputting HCAII concentrations 20% larger than those measured or foldamer concentrations 20% lower than those measured in the curve fitting program. Using a 2:1 binding model (two foldamers per

<table><caption>**Table 1** Dissociation constants of the complexes formed with HCAII determined by the fluorescence competition assay</caption><thead><tr><th>HCAII binder<sup>a</sup></th><th><i>K</i><sub>D</sub><sup>b</sup> (nM)</th></tr></thead><tbody><tr><td>22 (reference ligand)</td><td>10</td></tr><tr><td>15</td><td>1.5</td></tr><tr><td>16 (with Q<sup>5</sup>Phe<sup>10</sup>)</td><td>5.2</td></tr><tr><td>17 (with Q<sup>5</sup>Phe<sup>10</sup>, B<sup>Gpr11</sup>)</td><td>10</td></tr><tr><td>18 (with B<sup>Gpr6</sup>, Q<sup>5</sup>Phe<sup>10</sup>, B<sup>Gpr11</sup>)</td><td>7.2</td></tr><tr><td>19 (with B<sup>Gpr6</sup>, Q<sup>5</sup>Phe<sup>10</sup>, B<sup>6Gp11</sup>)</td><td>9.3</td></tr><tr><td>20 (with B<sup>Ida6</sup>, Q<sup>5</sup>Phe<sup>10</sup>)</td><td>7.4</td></tr></tbody></table>

<sup>a</sup> Some remarkable features are indicated in parenthesis. <sup>b</sup> Values were found to be repeatable within ±15% in duplicate experiments.

HCAII) also gave a better fit. The results are nevertheless presented as observed because the modifications that lead to a better fit have no experimental justification. Furthermore, these modifications do not change the trend and conclusions presented below.

In order to confirm these data, we tried to develop an alternate competition assay using BLI. A new biotinylated HCAII ligand **24** (Fig. 9a) was synthesized which, after immobilization on streptavidin sensors allowed for an accurate $K_D$ determination of its association with HCAII (Fig. 9). Immobilized **24** may in principle act as a reporter of the concentrations of free HCAII in solution. However, in this case as well, the kinetics were slow, a steady state regime was not reached. A calibration curve could in principle be produced by intercepting a value on the sensorgrams after a fixed amount of time instead of waiting until a steady state is reached. However, this proved not to be accurate enough to reliably determine the free HCAII concentration in solution.

Coming back to the $K_D$ values measured with the fluorescence competition assay, it appears that sequence **15**, with no added side chain in positions 6, 10 and 11, is the best binder and that all others bind similarly. These results should

<figure><img src="image_23.png" /><figcaption>**Fig. 9** (a) Structural formula of biotinylated HCAII ligand **24**. (b) BLI sensorgrams (black solid lines) of the titration of **24** immobilized on the streptavidin sensors at different HCAII concentrations. Calculated curves based on a 1:1 binding isotherm (colored dashed lines) fit with the measured values, yielding $K_D$ = 38.2 nM.</figcaption></figure>
<figure><img src="image_24.png" /></figure>

nevertheless be taken with caution. The foldamers suffer from low water solubility and a contribution from foldamer aggregation cannot be excluded. For instance, when the fluorescence titrations were performed at higher concentrations (e.g. [23] = 100 nM, [foldamer] = 200 nM, Fig. S19†), the apparent $K_D$ were higher than in Table 1, consistent with an effect of aggregation that reduces the effective foldamer concentration for binding to HCAII. The values in Table 1 could therefore also reflect that sequences **16–20** aggregate more than **15**, e.g. because of their hydrophobic Q⁵Ph¹⁰ residue.

## Discussion

It appears that none of the additional side chain combination of **16–20** result in a strong enhancement of their affinity for HCAII. It should be pointed that the structure-based design is intended to stabilize the complex, that is, to slow down complex dissociation. The effect of the additional side chains on the kinetics of complex formation, e.g. potentially slowing it down, remains unknown and is not taken into account in the computations. Early studies on HCAII ligands had shown that higher affinity correlated with faster complex formation rather than slower complex dissociation.³³ Another early study also reported the lack of effect of extending an HCAII ligand in the search for secondary binding sites, a result comparable to ours in an approach conceptually similar, albeit with much smaller molecules.³⁴ Finally, it may be that the very architecture of foldamers **15–20** makes it difficult for side chain modification to result in strong effects. These compounds consist of small nanomolar ligand to which is appended a much larger foldamer that has inherently no affinity for HCAII even in the low micromolar range – no induced CD is observed at 35 µM with **21** which lacks an HCAII ligand (Fig. S3†). In other words, the starting affinity of the foldamers for HCAII is too low to hope that a few modifications will bring it to an interesting range of $K_D$ value. In this respect, the choice of HCAII as a model system was perhaps not ideal. HCAII is a therapeutically relevant target and transmembrane isoforms HCAIX and HCAXII are overexpressed in some cancers and identified as potential targets as well.²⁴,²⁸ Nevertheless, we selected HCAII as a model system mainly for its robustness, easy overexpression, good crystal growth ability and the availability of simple nanomolar ligands that could act as tethers. To our knowledge, the vicinity of the HCAII active site is not involved in protein-protein interactions and deprived of any hotspot that may facilitate foldamer binding.

## Conclusions

In summary, starting from the crystal structure of the complex between HCAII and tetradecaamide foldamer 1, we have used computational tools to identify main chain and side chain modifications that may result in an extended foldamer-protein interface. New monomers and sequences incorporating these monomers were synthesized and several solid state structures of complexes with HCAII validated the design principles. We find that Q, B, and P main chain variations are interchangeable also

in the context of a foldamer-protein complex. We also find that side chains may generally be introduced independently from one another, a result that was consistent in both MD simulations and solid state structures. This behavior is in sharp contrast with that of peptides and peptidic foldamers in which a local change, e.g. a side chain modification, may also result in a different behavior of the main chain.¹⁷ The robustness of aromatic foldamer helices should therefore represent a good starting point for the structure-based design of protein binders that cover large protein surface areas. Nevertheless, the modifications explored in this study did not result into stronger associations nor did they deliver foldamers that would bind HCAII without a ligand or a covalent tether to mediate the interactions. The grand challenge of the *ab initio* design of an aromatic foldamer to bind a given protein surface remains unmet. In the meantime, other studies have revealed the potential of some aromatic helical foldamers to mimic α-helices¹³ᵃ or DNA double helices.¹⁰ Solid state structures have been obtained of complexes between chromosomal protein Sac7d and a DNA-mimic foldamer¹⁰ᵈ and between a fragment of ubiquitin ligase E6AP and a foldamer-peptide macrocycle.³⁵ In these complexes, no ligand or covalent tethering are involved. Furthermore, reliable $K_D$ determination methods are available. These structure thus represent new candidates to apply the structure-based design principles validated here. Steps in these directions are being made and will be reported in due course.

## Data availability

The data supporting this article have been included as part of the ESI.† Crystallographic data for 2·HCAII, 3·HCAII, 16·HCAII and 20·HCAII has been deposited at the Protein Data bank under accession numbers 9GAM, 9GAK, 9GAJ, 9HGB, respectively, and can be obtained from https://doi.org/10.2210/pdbXXXX/pdb where XXXX is the accession number.

## Author contributions

Except for the first and last authors, the author list groups the authors by institution and does not reflect a ranking of their contributions. LW, PSR and CD performed solution phase and solid phase syntheses. LW performed protein expression. BLE, LW and JS performed crystal growth. LF, JS and TG carried out crystallographic structure elucidation. LW and CD performed fluorescence and BLI titrations. ZL and YY performed computational studies. VP, YZ and IH supervised the research. IH, LW, CD and ZL wrote the manuscript. All authors reviewed and edited the manuscript and approved its final version.

## Conflicts of interest

There are no conflicts to declare.

## Acknowledgements

This work was supported by the Deutsche Forschungsgemeinschaft (DFG) through project CRC1309-C7
<figure><img src="image_25.png" /></figure>

(project ID 325871075) to I. H. Y. Y. and Y. Z. would like to acknowledge the support of the *National Institutes of Health* (NIH, R35GM127040 to Y. Z.). The work has benefited from the facilities and expertise of IECB Biophysical and Structural Chemistry platform (BPCS), CNRS UAR3033, Inserm US001, Univ. Bordeaux. We thank P. Mateus for valuable advice for the implementation of the fluorescence competition assay and Dr V. Morozov for providing some HCAII. We also thank J. Buratto, M. Soler-Lopez (ESFR beamline ID30B), M. Bowler (ESRF beamline ID30A-1), and M. Savko (SOLEIL beamline Proxima-2A) for assistance during X-ray data collection.

## Notes and references

‡ Sequences 7 and 11 are very similar – they differ only from one methylene group in the side chain of Q11 – so the fact that in both cases the MD trajectories significantly deviate from initial solid state coordinates is probably not a coincidence (Fig. S8). The exact reason is not clear but it appears that the combination of an anionic side chain on Q6 and a cationic side chain on Q11 may lead to this behavior. These two side chains are five units apart, i.e. exactly two helix turns, and thus only about 7 Å from each other. Salt bridges between them are observed especially in the case of 7 (40% of the simulation time) which has a longer cationic side chain, and most significant deviation from the solid state structure (Fig. S8).
§ For SPR measurements: 2: 98 DMSO/aqueous phosphate saline buffer PBS at pH 7.4 (vol/vol) at 25 °C, with PBS = 10 mM Na₂HPO₄, 1.8 mM KH₂PO₄, 2.7 mM KCl and 137 mM NaCl. For fluorescence measurements: 50 mM HEPES (4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid) at pH 7.2.

1. T. Seedorf, A. Kirschning and D. Solga, *Chem.-Eur. J.*, 2021, **27**, 7321.

2. (a) P. G. Baraldi, M. Del Carmen Núnez, A. Espinosa and R. Romagnoli, *Curr. Top. Med. Chem.*, 2004, **4**, 231; (b) Y. Hiraku, S. Oikawa and S. Kawanishi, *Nucleic Acids Res.*, 2002, **2**, 95; (c) A. Paul, P. Guo, D. W. Boykin and W. D. Wilson, *Molecules*, 2019, **24**, 1.

3. S. Baumann, J. Herrmann, R. Raju, H. Steinmetz, K. I. Mohr, S. Huttel, K. Harmrolfs, M. Stadler and R. Muller, *Angew. Chem., Int. Ed.*, 2014, **53**, 14605.

4. (a) E. Michalczyk, K. Hommernick, I. Behroz, M. Kulike, Z. Pakosz-Stepien, L. Mazurek, M. Seidel, M. Kunert, K. Santos, H. von Moeller, B. Loll, J. B. Weston, A. Mainz, J. G. Heddle, R. D. Sussmuth and D. Ghilarov, *Nat. Catal.*, 2023, **6**, 52; (b) S. M. Hashimi, *J. Antibiot.*, 2019, **72**, 785.

5. (a) J. Zeelen, M. van Straaten, J. Verdi, A. Hempelmann, H. Hashemi, K. Perez, P. D. Jeffrey, S. Halg, N. Wiedemar, P. Maser, F. N. Papavasiliou and C. E. Stebbins, *Nat. Microbiol.*, 2021, **6**, 392; (b) G. H. M. Salvador, T. R. Dreyer, A. A. S. Gomes, W. L. G. Cavalcante, J. I. Dos Santos, C. A. Gandin, M. de Oliveira Neto, M. Gallacci and M. R. M. Fontes, *Sci. Rep.*, 2018, **8**, 10317; (c) L. Jiao, S. Ouyang, M. Liang, F. Niu, N. Shaw, W. Wu, W. Ding, C. Jin, Y. Peng, Y. Zhu, F. Zhang, T. Wang, C. Li, X. Zuo, C. H. Luan, D. Li and Z. J. Liu, *J. Virol.*, 2013, **87**, 6829; (d) M. T. Murakami, E. Z. Arruda, P. A. Melo, A. B. Martinez, S. Calil-Elias, M. A. Tomaz, B. Lomonte, J. M. Gutierrez and R. K. Arni, *J. Mol. Biol.*, 2005, **350**, 416; (e) G. H. M. Salvador, T. R. Dreyer, W. L. Cavalcante, F. F. Matioli, J. I. Dos Santos, A. Velazquez-Campoy, M. Gallacci and M. R. Fontes, *Acta Crystallogr.*, 2015, **71**, 2066; (f) E. Mastrangelo, M. Pezzullo, D. Tarantino, R. Petazzi, F. Germani, D. Kramer, I. Robel, J. Rohayem, M. Bolognesi and M. Milani, *J. Mol. Biol.*, 2012, **419**, 198.

6. (a) I. Saraogi, J. A. Hebda, J. Becerril, L. A. Estroff, A. D. Miranker and A. D. Hamilton, *Angew. Chem., Int. Ed.*, 2010, **49**, 736; (b) T. Flack, C. Romain, A. J. P. White, P. R. Haycock and A. Barnard, *Org. Lett.*, 2019, **21**, 4433; (c) R. A. Dohoney, J. A. Joseph, C. Baysah, A. G. Thomas, A. Siwakoti, T. D. Ball and S. Kumar, *ACS Chem. Biol.*, 2023, **18**, 1510; (d) I. Arrata, C. M. Grison, H. M. Coubrough, P. Prabhakaran, M. A. Little, D. C. Tomlinson, M. E. Webb and A. J. Wilson, *Org. Biomol. Chem.*, 2019, **17**, 3861; (e) G. M. Burslem, H. F. Kyle, A. L. Breeze, T. A. Edwards, A. Nelson, S. L. Warriner and A. J. Wilson, *Chembiochem*, 2014, **15**, 1083; (f) S. Kumar and A. D. Hamilton, *J. Am. Chem. Soc.*, 2017, **139**, 5744; (g) J. P. Plante, T. Burnley, B. Malkova, M. E. Webb, S. L. Warriner, T. A. Edwards and A. J. Wilson, *Chem. Commun.*, 2009, **34**, 5091; (h) N. H. Stillman, J. A. Joseph, J. Ahmed, C. Zuwu Baysah, R. A. Dohoney, T. D. Ball, A. G. Thomas, T. C. Fitch, C. M. Donnelly and S. Kumar, *Nat. Commun.*, 2024, **15**, 3658.

7. (a) S. Kumar, M. Birol, D. E. Schlamadinger, S. P. Wojcik, E. Rhoades and A. D. Miranker, *Nat. Commun.*, 2016, **7**, 11412; (b) J. Ahmed, T. C. Fitch, C. M. Donnelly, J. A. Joseph, T. D. Ball, M. M. Bassil, A. Son, C. Zhang, A. Ledreux, S. Horowitz, Y. Qin, D. Paredes and S. Kumar, *Nat. Commun.*, 2022, **13**, 2273.

8. (a) J. Buratto, C. Colombo, M. Stupfel, S. J. Dawson, C. Dolain, B. Langlois d'Estaintot, L. Fischer, T. Granier, M. Laguerre, B. Gallois and I. Huc, *Angew. Chem., Int. Ed.*, 2014, **53**, 883; (b) M. Jewginski, T. Granier, B. Langlois d'Estaintot, L. Fischer, C. D. Mackereth and I. Huc, *J. Am. Chem. Soc.*, 2017, **139**, 2928; (c) M. Vallade, M. Jewginski, L. Fischer, J. Buratto, K. Bathany, J.-M. Schmitter, M. Stupfel, F. Godde, C. D. Mackereth and I. Huc, *Bioconjugate Chem.*, 2019, **30**, 54.

9. P. S. Reddy, B. Langlois d'Estaintot, T. Granier, C. D. Mackereth, L. Fischer and I. Huc, *Chem.-Eur. J.*, 2019, **25**, 11042.

10. (a) K. Ziach, C. Chollet, V. Parissi, P. Prabhakaran, M. Marchivie, V. Corvaglia, P. P. Bose, K. Laxmi-Reddy, F. Godde, J.-M. Schmitter, S. Chaignepain, P. Pourquier and I. Huc, *Nat. Chem.*, 2018, **10**, 511; (b) V. Corvaglia, D. Carbajo, P. Prabhakaran, K. Ziach, P. K. Mandal, V. D. Santos, C. Legeay, R. Vogel, V. Parissi, P. Pourquier and I. Huc, *Nucleic Acids Res.*, 2019, **47**, 5511; (c) V. Kleene, V. Corvaglia, E. Chacin, I. Forne, D. B. Konrad, P. Khosravani, C. Douat, C. F. Kurat, I. Huc and A. Imhof, *Nucleic Acids Res.*, 2023, **51**, 9629; (d) D. Deepak, J. Wu, V. Corvaglia, L. Allmendinger, M. Scheckenbach, P. Tinnefeld and I. Huc, *Angew. Chem., Int. Ed.*, 2025, **64**, e202422958.

11. (a) L.-G. Milroy, T. N. Grossmann, S. Hennig, L. Brunsveld and C. Ottmann, *Chem. Rev.*, 2014, **114**, 4695; (b) H. Lu, Q. Zhou, J. He, Z. Jiang, C. Peng, R. Tong and J. Shi, *Signal Transduct. Target. Ther.*, 2020, **5**, 213; (c) J. M. Ostrem, U. Peters, M. L. Sos, J. A. Wells and K. M. Shokat, *Nature*, 2066; (f) E. Mastrangelo, M. Pezzullo, D. Tarantino, R. Petazzi, F. Germani, D. Kramer, I. Robel, J. Rohayem, M. Bolognesi and M. Milani, *J. Mol. Biol.*, 2012, **419**, 198.
2013, **503**, 548; (d) C. V. Pagba, A. K. Gupta, A. K. Naji, D. van der Hoeven, K. Churion, X. Liang, J. Jakubec, M. Hook, Y. Zuo, M. Martinez de Kraatz, J. A. Frost and A. A. Gorf, *ACS Bio. Med. Chem. Au.*, 2022, **2**, 617.

12 T. Qi, V. Maurizot, H. Noguchi, T. Charoenraks, B. Kauffmann, M. Takafuji, H. Ihara and I. Huc, *Chem. Commun.*, 2012, **48**, 6337.

13 (a) M. Zwillinger, P. S. Reddy, B. Wicher, P. K. Mandal, M. Csékei, L. Fischer, A. Kotschy and I. Huc, *Chem.-Eur. J.*, 2020, **26**, 17366; (b) M. Zwillinger, P. Sőregi, F. Sanchez, C. Douat, M. Csékei, I. Huc and A. Kotschy, *J. Org. Chem.*, 2025, **90**, 3043; (c) X. Hu, S. J. Dawson, P. K. Mandal, X. de Hatten, B. Baptiste and I. Huc, *Chem. Sci.*, 2017, **8**, 3741.

14 (a) O. Ichihara, J. Barker, R. J. Law and M. Whittaker, *Mol. Inform.*, 2011, **30**, 298; (b) J. M. Mattheisen, C. Limberakis, R. B. Ruggeri, M. S. Dowling, C. W. am Ende, E. Ceraudo, T. Huber, C. L. McClendon and T. P. Sakmar, *J. Am. Chem. Soc.*, 2023, **145**, 11173; (c) A. Heikal, Y. Nakatani, W. Jiao, C. Wilson, D. Rennison, M. R. Weimar, E. J. Parker, M. A. Brimble and G. M. Cook, *Bioorg. Med. Chem. Lett.*, 2018, **28**, 2239; (d) X. Wu, Y. Zhang, S. Liu, C. Liu, G. Tang, X. Cao, X. Lei and J. Peng, *Bioorg. Chem.*, 2022, **127**, 105921; (e) B. C. Doak, R. S. Norton and M. J. Scanlon, *Pharmacol. Ther.*, 2016, **167**, 28.

15 (a) M. F. Schmidt and J. Rademann, *Trends Biotechnol.*, 2009, **27**, 512; (b) M. Mondal and A. K. H. Hirsch, *Chem. Soc. Rev.*, 2015, **44**, 2455; (c) P. Frei, R. Hevey and B. Ernst, *Chem.-Eur. J.*, 2019, **25**, 60.

16 W. Lu, M. Kostic, T. Zhang, J. Che, M. P. Patricelli, L. H. Jones, E. T. Chouchani and N. S. Gray, *RSC Chem. Biol.*, 2021, **2**, 354.

17 (a) M. E. Perrin, B. Li, J. Mbianda, M. Bakail, C. André, G. Moal, P. Legrand, V. Ropars, C. Douat, F. Ochsenbein and G. Guichard, *Chem. Commun.*, 2023, **59**, 8696; (b) J. Mbianda, M. Bakail, C. André, G. Moal, M. E. Perrin, G. Pinna, R. Guerois, F. Becher, P. Legrand, S. Traoré, C. Douat, G. Guichard and F. Ochsenbein, *Sci. Adv.*, 2021, **7**, eabd9153; (c) H. Wang, R. S. Dawber, P. Zhang, M. Walko, A. J. Wilson and X. Wang, *Chem. Sci.*, 2021, **17**, 5977.

18 M. Vallade, P. Sai Reddy, L. Fischer and I. Huc, *Eur. J. Org Chem.*, 2018, 5489.

19 D. Bindl, P. K. Mandal and I. Huc, *Chem.-Eur. J.*, 2022, **28**, e202200538.

20 D. Bindl, E. Heinemann, P. K. Mandal and I. Huc, *Chem. Commun.*, 2021, **57**, 5662.

21 (a) S. Dengler, P. K. Mandal, L. Allmendinger, C. Douat and I. Huc, *Chem. Sci.*, 2021, **12**, 11004; (b) V. Corvaglia, F. Sanchez, F. S. Menke, C. Douat and I. Huc, *Chem.-Eur. J.*, 2023, **29**, e202300898.

22 V. Corvaglia, J. Wu, D. Deepak, M. Loos and I. Huc, *Chem.-Eur. J.*, 2024, **28**, e202303650.

23 (a) D. Rooklin, C. Wang, J. Katigbak, P. S. Arora and Y. Zhang, *J. Chem. Inf. Model.*, 2015, **55**, 1585; (b) J. Katigbak, H. Li, D. Rooklin and Y. Zhang, *J. Chem. Inf. Model.*, 2020, **60**, 1494.

24 (a) T. Gokcen, I. Gulcin, T. Ozturk and A. C. Goren, *J. Enzyme Inhib. Med. Chem.*, 2016, **31**, 180; (b) A. Maresca, C. Temperini, L. Pochet, B. Masereel, A. Scozzafava and C. T. Supuran, *J. Med. Chem.*, 2010, 53-335; (c) C. T. Supuran, *Nat. Rev. Drug Discovery*, 2008, 7, 168.

25 D. Gfeller, O. Michielin and V. Zoete, *Nucleic Acids Res.*, 2013, **41**, D327.

26 O. Trott and A. J. Olson, *J. Comput. Chem.*, 2010, **31**, 455-461.

27 D. A. Case, H. M. Aktulga, K. Belfon, D. S. Cerutti, G. A. Cisneros, V. W. D. Cruz eiro, N. Forouzesh, T. J. Giese, A. W. Götz, H. Gohlke, S. Izadi, K. Kasavajhala, M. C. Kaymak, E. King, T. Kurtzman, T.-S. Lee, P. Li, J. Liu, T. Luchko, R. Luo, M. Manathunga, M. R. Machado, H. M. Nguyen, K. A. O'Hearn, A. V. Onufriev, F. Pan, S. Pantano, R. Qi, A. Rahnamoun, A. Risheh, S. Schott-Verdugo, A. Shajan, J. Swails, J. Wang, H. Wei, X. Wu, Y. Wu, S. Zhang, S. Zhao, Q. Zhu, T. E. Cheatham III, D. R. Roe, A. Roitberg, C. Simmerling, D. M. York, M. C. Nagan and K. M. Merz Jr, *J. Chem. Inf. Model.*, 2023, **63**, 6183.

28 T. O. Wambo, L. Y. Chen, S. F. McHardy and A. T. Tsin, *Biophys. Chem.*, 2016, **214-215**, 54.

29 J. A. Maier, C. Martinez, K. Kasavajhala, L. Wickstrom, K. E. Hauser and C. Simmerling, *J. Chem. Theory Comput.*, 2015, **11**, 3696.

30 J. Wang, R. M. Wolf, J. W. Caldwell, P. A. Kollman and D. A. Case, *J. Comput. Chem.*, 2004, **25**, 1157.

31 Z. Liu, A. M. Abramyan and V. Pophristic, *New J. Chem.*, 2015, **39**, 3229.

32 P. Koutnik, E. G. Shcherbakova, S. Gozem, M. G. Caglayan, T. Minami and P. Anzenbacher, *Chem*, 2017, **2**, 27.

33 R. W. King and A. S. V. Burgen, *Proc. R. Soc. Lond. B Biol. Sci.*, 1976, **193**, 107.

34 A. Jain, S. G. Huang and G. M. Whitesides, *J. Am. Chem. Soc.*, 1994, **116**, 5057.

35 S. Dengler, R. T. Howard, V. Morozov, C. Tsiamantas, W. Huang, Z. Liu, C. Dobrzanski, V. Pophristic, S. Brameyer, C. Douat, H. Suga and I. Huc, *Angew. Chem., Int. Ed.*, 2023, **62**, e202308408.