Structural Insights into the Interaction of Clinically Relevant Phosphorothioate mRNA Cap Analogs with Translation Initiation Factor 4E Reveal Stabilization via Electrostatic Thio-Effect

mRNA-based therapies and vaccines constitute a disruptive technology with the potential to revolutionize modern medicine. Chemically modified 5′ cap structures have provided access to mRNAs with superior translational properties that could benefit the currently flourishing mRNA field. Prime examples of compounds that enhance mRNA properties are antireverse cap analog diastereomers that contain an O-to-S substitution within the β-phosphate (β-S-ARCA D1 and D2), where D1 is used in clinically investigated mRNA vaccines. The compounds were previously found to have high affinity for eukaryotic translation initiation factor 4E (eIF4E) and augment translation in vitro and in vivo. However, the molecular basis for the beneficial “thio-effect” remains unclear. Here, we employed multiple biophysical techniques and captured 11 cap analog-eIF4E crystallographic structures to investigate the consequences of the β-O-to-S or -Se substitution on the interaction with eIF4E. We determined the SP/RP configurations of β-S-ARCA and related compounds and obtained structural insights into the binding. Unexpectedly, in both stereoisomers, the β-S/Se atom occupies the same binding cavity between Lys162 and Arg157, indicating that the key driving force for complex stabilization is the interaction of negatively charged S/Se with positively charged amino acids. This was observed for all structural variants of the cap and required significantly different conformations of the triphosphate for each diastereomer. This finding explains why both β-S-ARCA diastereomers have higher affinity for eIF4E than unmodified caps. Binding affinities determined for di-, tri-, and oligonucleotide cap analogs suggested that the “thio-effect” was preserved in longer RNAs. Our observations broaden the understanding of thiophosphate biochemistry and enable the rational design of translationally active mRNAs and eIF4E-targeting drugs.


■ INTRODUCTION
Phosphorothioate (PS) modification of (oligo)nucleotides is one of the most frequently used synthetic approaches for investigating nucleic acid metabolism or for improving the pharmacokinetics of therapeutic or diagnostic molecules. 1−6 Many applications of thio-modified oligonucleotides, such as the elucidation of ion-binding sites on RNA 7 and identification of phosphate sites important for ligand or protein binding or for catalysis, 8−12 assume that this modification has little or no effect on the secondary structure of nucleic acids. However, it has been shown that incorporation of the PS moiety alters the negative charge distribution 13 and can significantly alter the conformation of nucleic acids, thus changing their biological activity. 14, 15 Furthermore, substitution of a nonbridging oxygen in the phosphate residue with sulfur usually generates an additional stereogenic center at the phosphorus atom yielding a pair of diastereomers, 16,17 which are often differentially recognized by specific proteins or enzymes. 18−24 Interestingly, a natural PS backbone modification was shown to be an epigenetic mark in bacterial DNA, 25−27 and recently it was also found in both prokaryotic and eukaryotic rRNAs. 28 The 7-methylguanosine (m 7 G) cap structure is present at the 5′ end of all eukaryotic mRNAs ( Figure 1A), protecting them from rapid degradation by 5′-exonucleases and regulating many important gene expression steps, including mRNA splicing, nuclear export, initiation of translation, and mRNA turnover. 29,30 Chemical modification of the phosphate moieties in the 5′ cap structure has emerged as a particularly appealing strategy to stabilize transcripts and improve their translational activity, 31,32 providing a new way to improve mRNA based therapeutics. 33−36 Recently, several clinical trials on the development of mRNA vaccines against SARS-CoV-2 gave very promising results. 37,38 We have previously reported a series of mRNA cap analogs carrying a PS moiety as reagents for mRNA preparation. 39 We have shown that the two analogs termed β-S-ARCA (ARCA stands for Anti-Reverse Cap Analog, meaning that the ribose of 7-methylguanosine was blocked by a methyl group to prevent priming of in vitro transcription at m 7 G, which produces inactive mRNA), carrying an O-to-S substitution within the β-phosphate of the triphosphate chain, have particularly beneficial biological properties. The analogs, differing in the absolute configurations at the β-phosphate ( Figure 1B), were labeled D1 and D2 according to their elution order from the RP HPLC column. β-S-ARCA D1 and D2 had, to different extents, decreased susceptibility to cleavage by Dcp2 decapping enzyme, and increased affinity for eukaryotic translation initiation factor 4E (eIF4E). Interestingly, isomer D1 had a 2-fold higher affinity for eIF4E than D2 and 4-fold higher affinity than the unmodified parent compound, 39 which resulted in higher translational efficiency in vitro and in cultured cells for mRNAs incorporating either of these analogs. 31 This indicates that the "thio-effect" resulting from a single-atom substitution is sufficient to evoke biologically relevant effects ( Figure 1C). mRNAs incorporating β-S-ARCAs also had elongated cellular half-lives and showed augmented protein expression. Because of these favorable properties, β-S-ARCA D1 has been applied in experimental mRNA-based anticancer vaccines 40,41 and is currently used as a capping reagent for mRNAs under clinical investigation in multiple trials. However, the R P /S P -absolute configurations of the β-modified PS cap analogs and the electronic and structural consequences of this modification have not been elucidated so far.
Here, we performed X-ray crystallography structural studies on a series of therapeutically relevant β-phosphorothioate and β-phosphoroselenoate mRNA 5′ cap analogs ( Figure  1D), 39,42,43 including β-S-ARCA isomers, to elucidate the molecular mechanism of their tight interaction with eIF4E as well as the differences in the recognition of diastereomers. 31,39 To this end, we determined 11 original X-ray crystal structures of eIF4E in complexes with ligands carrying the O to S or O to Se substitutions at the β-position of the triphosphate chain, which enabled us to define both the differences and the common features of their recognition. We also determined binding affinities of model compounds and RNAs incorporating them for eIF4E, using two independent methods to verify whether the observed "thio-effect" is an artifact observed only for dinucleotide cap analogs or is preserved in longer RNAs.

■ RESULTS AND DISCUSION
Chemical Modification of the Cap Has No Effect on eIF4E Protein Conformation. To explain the "thio-effect" on the molecular level, we determined the structures of β-S-ARCA (m 2 7,2 ′ -O Gpp S pG) D1 and D2 ( Figure 2) and related compounds ( Figure 3) in complex with murine eIF4E (meIF4E; residues 28− 217), which is 99% identical and 100% similar to the human protein.
All the structures were solved at 1.7−2.2 Å resolution, most of them in the P1 space group (Table S1) with four copies of the protein-cap complex in a unit cell, except for complexes with m 2 7,2 ′ -O Gpp S pA D1 and D2, which were solved in the C2 space group with two copies of ACS Chemical Biology pubs.acs.org/acschemicalbiology Articles the complex (the overall crystal packing was very similar in both space groups). We were able to find noncrystallographic symmetry operators between those copies and verified that their conformation and interaction patterns are consistent (except for m 7 Gpp S pApG D2, which is discussed later). Protein conformations in the structures of β-S-ARCA diastereomer complexes are very similar to that of murine eIF4E in complex with an unmodified cap (m 7 GpppG) reported earlier ( Figure 2A; RMSD 1.41 and 1.29 Å for D1 and D2, respectively, RMSD 0.61 Å between D1 and D2; all the atoms were included). 44 The Absolute Configuration of β-Phosphorus Affects the Cap-Binding Mode by eIF4E. Surprisingly, the ligand conformations adopted by β-S-ARCA D1 and D2 are quite different. In both diastereomers, 7-methylguanosine and γ-phosphate moieties have similar binding modes, similar to the previously observed structures of mono-and dinucleotide caps in complex with either murine or human eIF4E. 44,45 However, the position and conformation of βand αphosphates in β-S-ARCA D1 and D2 are markedly distinct from each other ( Figure 2A). β-S-ARCA D1 adopts a conformation similar to that determined previously for m 7 GpppG, 44 with the sulfur atom of the β-thiophosphate placed in a small cavity between a positively charged Lys162 and Arg157 (Figure 2A). β-S-ARCA D2 is bound in a different conformation, in which the β-phosphate appeared to rotate around the P β −O β,γ bond by ∼120°compared to D1, and αphosphate points outside of the binding pocket (Movie 1). Interestingly, for both diastereomers, sulfur atoms appear to be localized in the same position of the cap binding pocket, which suggests that the stereochemistry of β-phosphorus determines the arrangement of the triphosphate bridge. Unfortunately, the initial data did not enable us to clearly distinguish sulfur from oxygen based on the electron density; hence, an unambiguous assignment of the absolute configuration of the β-phosphorus atom required additional experiments. To this end, we cocrystallized eIF4E in complex with the corresponding phosphoroselenoate cap analogues (β-Se-ARCA = m 2 7,2 ′ -O Gpp Se pG) and collected diffraction data at the wavelength corresponding to the maximum absorption of selenium (0.97957 Å − 12.657 keV). 46 This enabled unambiguous determination of the position of selenium atoms in the complex.
The phosphoroselenoate cap analogs have biological properties similar to those of phosphorothioate analogs ( Table 1, Table S2;   42 so we expected that they should be closely related in structure and the molecular mechanism of binding to eIF4E. Indeed, the shapes of F o −F c electron densities for β-Se-ARCAs in the cap binding pocket are very similar to those observed for the corresponding β-S-ARCAs. The anomalous difference electron density maps revealed strong signals in both β-Se-ARCAs complexes ( Figure 2B), confirming the absolute configurations around the β-phosphorus center for D1 and D2 as R P and S P , respectively, and clearly positioning the selenium between the side chains of Lys162 and Arg157 regardless of the absolute configuration. The occupancy of this positively charged cavity by the sulfur/ selenium atom in β-S-ARCA D2 and β-Se-ARCA D2 is enabled by a significant conformational rearrangement of the 5′,5′-triphosphate chain, resulting in different positioning of the α-phosphate (and likely also a guanosine moiety, for which the electron density is unfortunately not visible in any of the analyzed structures) compared to m 7 GpppG or corresponding isomers D1.
The Nature of the "Thio-Effect" Is Mainly Electrostatic. Sulfur and selenium atoms form hydrogen bonds of similar strength to those with oxygen atoms. 47 Therefore, their strong tendency to occupy space between the two positively charged groups is most likely caused by the electrostatic attraction between opposite charges. This, in turn, suggests that the negative charge in β-phosphorothioate and β-phosphoroselenoate moieties is localized at the heteroatom, instead of being equally distributed between the heteroatom and nonbridging β-oxygen. This hypothesis is consistent with an extensive analysis based on NMR, X-ray, and vibrational spectroscopy evidence, 13 as well as ab initio calculations. 48 These calculations led to the conclusion that the most accurate representation of charge localization in phosphorothioate esters and anhydrides in solution is a resonance structure with a single P−S bond and the negative charge localized at the sulfur atom ( Figure 2C). We envisage that in the case of the studied cap analogues, the negatively charged S/Se atoms from the modified β-phosphate residue form a very strong contact with positively charged Arg157 and Lys162, which can be classified as a salt bridge. The presence of these strong stabilizing interactions involving S/Se atoms arguably explains why both the D1 and D2 isomers display increased affinity for eIF4E compared to m 7 GpppG. On the other hand, the conformational change observed for the triphosphate in D2 isomers entails different positioning of the α-phosphate in those complexes, which disrupts its interaction with Lys162 and is most likely responsible for the lower affinity to eIF4E observed for the D2 diastereomer as compared to D1. Comparison of the "z-score normalized" 49 B-factors of D1 and D2 ligands in complexes with eIF4E show that the D2 isomer is generally more flexible than the D1 isomer, particularly at the β-phosphate, which is consistent with the number of observed contacts with the protein.
To further examine the significance of charge distribution within the phosphorothioate residue in arranging the conformation of the 5′,5′-triphosphate bridge, we determined crystal structures of a diastereomeric pair of another β-phosphorothioate cap analog, in which the m 7 G ribose ring is constrained by a 2′,3′-O,Oisopropylidene substituent (iPr-m 7 Gpp S pG). 43 We have previously obtained iPr-m 7 Gpp S pG D1 and D2 as synthetically simpler alternatives to β-S-ARCA but found that iPr modification slightly destabilizes eIF4E/cap complexes, likely by affecting the conformational equilibrium of m 7 G pucker, which results in relatively lower translational efficiencies for capped mRNAs. 43 Interestingly, we determined that the D2 isomer of iPr-m 7 Gpp S pG has higher affinity for eIF4E than the D1 isomer, in contrast to what was observed for β-S-ARCA, suggesting that the conformation of the triphosphate bridge    35 n.d. no binding may also be affected by the bulky substituent. 43 Therefore, we were interested to obtain structural insights into iPr-m 7 Gpp S pG binding by eIF4E at the molecular level, to see whether it can provide a rationale for the observed biophysical and biological properties. As expected, the crystal structures of eIF4E/iPr-m 7 Gpp S pG complexes showed unnatural sugar puckering in m 7 G, resulting in the displacement of the γ-phosphate toward the C-terminal loop (residues 206−212) by ca. 2 and 3 Å for diastereomers D1 (R P ) and D2 (S P ), respectively ( Figure  3). Consequently, for the D1 (R P ) isomer, the positioning of the βphosphate moiety is also altered, resulting in suboptimal contact with Arg157 and Lys162, which is probably the main reason for its lower affinity for eIF4E as compared to the D2 (S P ) isomer. Despite the significant differences in the conformation of the triphosphate bridge between iPr-m 7 Gpp S pG diastereomers, the sulfur atoms are clearly attracted toward Lys162 and Arg157 side chains. This observation again confirms that the electrostatic interaction between basic amino acids and heteroatoms at the β-position of the triphosphate chain plays a dominant role in determining the arrangement of the cap analogues in complexes with eIF4E. Chemically Modified Cap Analogs Reveal a New Potentially Druggable Binding Site on the Surface of eIF4E. Notably, in the case of the iPr-m 7 Gpp S pG D2 (S P ) complex, the difference electron density map is very well-defined for the entire ligand molecule, revealing the binding mode of the second nucleoside ( Figure 3). This unprecedented conformation of the cap in complex with eIF4E is stabilized by several hydrogen bonds to Lys54 and Asn59, engaging the Watson−Crick edge of guanine (N 1 , N 2 , and O 6 ) as well as an edge-to-face contact between the nucleobase and the −CONH 2 group of Asn50 ( Figure S1A). Binding of the second base is accommodated by rotation of the Asn50 side chain, which exposes a small cavity formed by the N-terminal loop (residues 50−55). After a careful inspection of the refined structures of (S P )-m 2 7,2 ′ -O Gpp S pG and (S P )m 2 7,2 ′ -O Gpp Se pG complexes (D2 isomers), we observed some residual electron density peaks in a pocket formed by the N-terminal loop, which could not be satisfactorily modeled by water molecules (in contrast to the corresponding structures of R P isomers). The superposition of this map with the structure of the iPr-m 7 Gpp S pG D2 (S P ) complex shows a very good colocalization of an unmodeled electron density with the guanine part of isopropylidene cap analogs ( Figure S1B), suggesting that a similar binding mode may be relevant to m 2 7,2 ′ -O Gpp S pG and m 2 7,2 ′ -O Gpp Se pG D2 (S P ) complexes, but these conformers are apparently not fully occupied.
To verify whether the observed conformational differences between D1 and D2 isomers are more general phenomena related to the presence of sulfur (or selenium) atoms ("thio-effect") or are observed due to the particular affinity of guanine to interact with the N-terminal loop pocket, we cocrystallized eIF4E in complexes with another pair of β-PS cap analogs, m 2 7,2 ′ -O Gpp S pA D1 and D2, which contained adenosine instead of guanosine at the transcription start site (TSS or +1 nucleotide; Figure 1A,D). Adenine differs from guanine in the hydrogen bonding pattern at the Watson−Crick edge; therefore, it should not form similar contacts with Lys54 and Asn59. Indeed, in these structures, we did not observe any unmodeled electron density peaks in the N-terminal loop region nor the rotation of the Asn50 side chain, suggesting the absence of nucleobase in the pocket. However, the simulated annealing F o −F c omit map showed some residual electron density peaks in the C-terminal loop region of eIF4Em 2 7,2 ′ -O Gpp S pA D1 (R P ) complex. Importantly, the density does not colocalize with the adenosine portion of an unmodified m 7 GpppA cap in complex with the human eIF4E reported earlier ( Figure S2). 45 Rather, it suggests that the α-phosphate is rotated around P α -O α,β by ∼120°compared to m 7 GpppA, placing the second nucleoside in a different region of the loop. The 3′-OH group in the TSS (which is the site of attachment of the RNA chain) for both m 7 GpppA and m 2 7,2 ′ -O Gpp S pA D1 is buried deep in the cap binding pocket, which suggests that these conformations cannot be accommodated by fulllength mRNA. Nonetheless, the structures reveal new structural features in eIF4E that may be targeted for the design of small molecule translational inhibitors with anticancer activity. 50−52 Importantly, for both diastereomers of m 2 7,2 ′ -O Gpp S pA, the sulfur atom is localized between Arg157 and Lys162 side chains and the conformations of the triphosphate bridge were essentially identical to the conformations of corresponding β-S-ARCAs ( Figure 3). The absolute configurations of β-phosphorus atoms also match those of m 2 7,2 ′ -O Gpp S pG (namely, R P for the D1 isomer and S P for the D2 isomer).
Trinucleotide mRNA 5′ Cap Analogs Provide a More Detailed View of the Structure of the Translation Initiation Complex. We next aimed to verify whether the two distinct conformations observed for β-PS caps are also relevant to longer RNAs. To provide better model molecules to study the impact of the RNA chain on the cap-binding mode, we synthesized trinucleotide cap analogs (the third nucleotide was attached through a 3′-5′phosphodiester linkage, identical to the natural mRNA 5′ end structure; Figure 1D, compounds m 7 GpppApG, m 7 Gpp S pApG D1 and D2). We obtained high-resolution structures of complexes with m 7 GpppApG and both diastereomers of m 7 Gpp S pApG. Notably, these are the first experimentally determined models of eIF4E in complex with a portion of mRNA longer than two nucleotides. The overall conformations of the triphosphate bridges of trinucleotide analogs are virtually identical to those of their dinucleotide counterparts (m 7 GpppA and m 2 7,2 ′ -O Gpp S pA; Figure 3), including the position of the sulfur atom, which means that the "thio-effect" is also preserved for trinucleotides. The structures also allowed us to assign the absolute configurations of the β-phosphorus atom for both PS diastereoisomers as R P for the D1 isomer and S P for the D2 isomer, which are analogous to the configurations determined for dinucleotides. The electron density maps for m 7 GpppApG and ACS Chemical Biology pubs.acs.org/acschemicalbiology Articles m 7 Gpp S pApG D1 (R P ) are interpretable only for the 7-methylguanosine and triphosphate bridge portions, while the rest of the ligand is disordered with no apparent electron density for adenosine near the C-terminal loop. This suggests that the first transcribed nucleotide (TSS) and downstream nucleotides in mRNA are generally not involved in any specific contacts with eIF4E or that they form multiple weaker contacts that do not stabilize any particular conformation. This is in contrast to the eIF4E-m 7 GpppA X-ray structure, which showed an interaction between adenine and the Cterminal loop of the protein. 45 However, as mentioned above, the conformation of adenosine observed in dinucleotide complexes could not be adopted by trinucleotide or longer RNAs, since its 3′-OH group is buried in the C-terminal loop, which makes no room for phosphodiester bonds and downstream nucleotides. Therefore, we hypothesize that the presence of the mRNA chain attached to the cap might influence the conformation of the flexible loop by steric hindrance. To verify this, we recorded 1 H− 15 N HSQC spectra of eIF4E/m 7 GTP and eIF4E/m 7 GpppApG complexes and compared the chemical shifts of amide NH groups between those complexes. Indeed, large chemical shift perturbations (CSPs) were observed for resonances corresponding to residues 200−204, which form a hinge region of the C-terminal loop, indicating some conformational differences between the complexes (Figure 4, Figure S4). A similar effect was observed for trinucleotide cap analogs differing in the methylation status of adenosine at the TSS (m 7 GpppApG vs m 7 Gppp m6 A m pG). 53 Although the identity and methylation status of the TSS nucleotide have minor effects on the affinity of the cap for eIF4E, it might have some biological implications, since the Cterminal loop of eIF4E is phosphorylated at Ser209 in response to mitogenic stimuli and cytokines. 54 Comparison of the HSQC spectra of mono-and trinucleotide complexes also showed significant shift differences for some amides (residues G160, L117, and N118) located on the surface of eIF4E. These CSPs indicate the position of the TSS (+1) and +2 nucleotides in the eIF4E/m 7 GpppApG complex, and hence imply a possible path taken by the RNA chain in complex with the full-length transcript ( Figure 4). Surprisingly, in the cap binding pocket of one of the eIF4E/m 7 Gpp S pApG D2 (S P ) complex subunits (chain A), we observed a very well-defined electron density, which allowed us to model almost the entire ligand, excluding only the nucleobase of the third nucleotide (Figures 3 and 5). This particular conformation of the cap is stabilized by an intramolecular hydrogen bond between m 7 G ribose (3′-OH) and phosphodiester linkage and probably by the contact between ribose of the third nucleotide and aspartate of the symmetry-related protein molecule (the one that belongs to the adjacent crystallographic unit cell −1 + x, y, z), and as such may be only an artifact of crystal packing. Nonetheless, it provides the idea of a possible alignment of the mRNA 5′ end related to eIF4E in biologically relevant complexes. This is consistent with the NMR chemical shift perturbation experiment ( Figures 5 and S3A), showing that an unnatural alignment of the triphosphate bridge does not have to affect the path taken by the RNA chain along the eIF4E surface. Such a well-defined electron density corresponding to +1 and +2 nucleotides of the ligand is not observed in the other copy of the complex present in a crystallographic unit cella noncrystallographic symmetry (NCS)-related subunit (chain B), which differs only slightly (by ∼3.3 Å) in the position of the aspartate from analogous symmetry-related molecule (x, −1 + y, z). Instead, chain B revealed a new and interesting binding pattern. At a hydrogen bond distance from the ε-amine group of Lys162, we observed a large electron density (4.6 RMSD peak) in the F o −F c simulated annealing omit map, which can be attributed to the phosphorus atom of the phosphodiester linkage ( Figure S3B). Such an interaction might compensate for the lack of a hydrogen bond between Lys162 and αphosphate. It also allows the RNA chain to take the path along the βstrand core of eIF4E, similar to that observed for m 7 GpppApG by NMR CSP mapping (Figure 4).
The Conformation of the 5′,5′-Triphosphate Bridge Observed for S P (D2) Isomers May Be Biologically Relevant.
To obtain a more comprehensive image of β-PS modification in the context of cap-eIF4E interactions, we superimposed all the structures determined in this work (excluding isopropylidene analogs). This superposition revealed the existence of two distinct alignments of the triphosphate bridge, dependent solely on the β-phosphorus stereochemistry ( Figure 6). Comparing our structures with other X-ray structures of cap-bound translation initiation factor 4E deposited in the Protein Data Bank, it was found that the conformation characteristic of R P β-phosphorothioates is virtually identical to that of unmodified caps (m 7 GTP, m 7 GpppG, and m 7 GpppA) in complex with mammalian eIF4E proteins. Interestingly, the conformation ACS Chemical Biology pubs.acs.org/acschemicalbiology Articles characteristic of S P isomers closely resembles that of m 7 GpppG in complex with the parasite (Schistosoma mansoni, Sch) eIF4E protein (32% sequence identity and 51% similarity to human eIF4E). 56 Similar to the S P β-phosphorothioate cap analogs in mammalian eIF4E, the α-phosphate of m 7 GpppG in the Sch-eIF4E complex points away from the lysine amine group, directing the second nucleoside to the N-terminal loop ( Figure S5). It has been suggested that such conformation of the cap allows for the movement of the Cterminal loop of Sch-eIF4E toward the cap-binding pocket to clasp the ligand and consequently impacts the pathway taken by the mRNA chain in the translation initiation complex. 56 However, mammalian eIF4Es have considerably shorter C-terminal loops, which cannot fully clasp the ligand. The structure of the eIF4E/m 7 Gpp S pApG D2 (S P ) complex suggests that the RNA phosphodiester backbone is flexible enough to bypass the changes in alignment of the triphosphate bridge caused by PS modification.
The "Thio-Effect" Is Preserved in Longer RNAs. Unfortunately, we did not succeed in crystallizing longer capped RNA fragments in complex with eIF4E. Therefore, to verify whether the "thio-effect" is preserved in longer RNAs, we focused on analyzing the thermodynamic properties of various cap-eIF4E complexes. First, to corroborate the structural information with thermodynamic data, we determined the dissociation constants (K D ) for eIF4E-m 2 7,2 ′ -O Gpp S pA complexes by fluorescence quenching titration (FQT) and compared them to values previously determined for β-S-ARCA and β-Se-ARCA (Tables 1 and S2). Similar to guanine analogs, R P -β-PS modification in m 7 Gpp S pA D1 stabilizes the complex, yielding 4.0-fold lower K D compared to m 7 GpppA. Interestingly, the stabilization effect for S P modification in m 7 Gpp S pA D2 was slightly larger than that for the corresponding guanine analog (2.5-fold increase in binding compared to 1.8-fold increase for guanine analog), although no additional contacts were observed. Using the FQT method, we also determined the binding affinities of trinucleotide cap analogs, finding that the stabilizing "thio-effect" is preserved for those compounds, although the differences between stereoisomers are less pronounced. For trinucleotides, β-PS D1 increases binding affinity to eIF4E by ∼3.5fold, which is slightly less than for dinucleotides, while D2 increases binding affinity by 2.8-fold, which is slightly higher than that for dinucleotides. Finally, we prepared 35-nt fragments of RNA capped with ARCA or β-S-ARCA D1 by in vitro transcription to determine their K D values for eIF4E. The FQT method was not suitable for this purpose because of the large sample consumption and high inner filter effect occurring in the case of RNA. To overcome this issue, we developed a microscale thermophoresis (MST) competition assay for eIF4E ( Figure S6), based on a fluorescently labeled tightly binding probe [m 7 Gp 5 OC 3 (5)FAM, Supporting Information]. The assay enabled indirect comparison of eIF4E-cap dissociation constants for different ligands, based on the apparent binding constant values K d,app derived from fluorescent probe replacement experiments ( Figure  S6C). To validate the MST assay, we first determined the K d,app for ligands that had already been characterized by FQT. These two data sets (Table S2) Figure S6B). The MST assay independently confirmed that the stabilization effect of β-PS was observed for both di-and trinucleotides. Finally, we determined K d,app for β-S-ARCA D1-and ARCA-capped RNA. In this case, the β-PS modification also provided a statistically significant stabilizing effect, albeit less pronounced than in dinucleotides. Namely, RNA capped with β-S-ARCA D1 binds to eIF4E with an affinity 2.2-fold higher than ARCA-capped RNA. The thio-effect is then generally conserved from di-through tri-to oligonucleotide cap analogs, but its contribution to overall complex stability seems to be slightly diminished in the presence of the RNA chain.

■ CONCLUSIONS
In conclusion, a series of therapeutically relevant βphosphorothioate-modified analogs of mRNA 5′ ending in complex with eukaryotic translation initiation factor 4E was structurally characterized. On the basis of the crystal structures of phosphoroselenoate analogs and their anomalous diffraction, we were able to unambiguously determine the position of selenium atoms and thus assign the absolute configuration of the asymmetric phosphorus atom in each diastereomer as well as indirectly confirm assignments for phosphorothioate analogs. The local environment of sulfur or selenium formed by positively charged amino acids clearly suggests that the negative charge in the modified phosphate residues is mainly localized at the S/Se atoms. The ionic contact between S/Se and basic amino acids was identified as a major factor in complex stabilization and the driving force for the alignment of ligands in the cap binding pocket. As a consequence, phosphorothioate cap analog diastereomers adopted two different sets of conformations, one of which has never been observed before in any mammalian eIF4E complex. In this conformation, adopted preferentially by S P β-phosphorothioate cap analogs, the triphosphate bridge is bound by the protein only through β and γ phosphates, lacking interaction between α-phosphate and Lys162. This is the likely reason for the lower affinity of S P isomers in comparison to R P isomers, but the electrostatic interactions between a negatively charged sulfur atom and two basic amino acids (Arg157 and Lys162) provide higher stabilization energy than the similar interactions involving nonmodified α and β phosphate residues with a delocalized charge in m 7 GpppG. The missing interaction involving α-phosphate is compensated to some extent in RNAs capped with β-S-ARCA D2 (S P ), most likely by a contact between Lys162 and phosphodiester linkage, which is not present in dinucleotide analogs. On the other hand, the lower number of contacts involving the triphosphate bridge of S P diastereomers might facilitate dissociation of the transcript from eIF4E, which was suggested to be necessary for the loading of mRNA onto the 43S preinitiation complex and subsequent start codon scanning. 57,58 This conclusion is consistent with the polysomal distribution experiment reported earlier, which showed that mRNAs capped with m 2 7,2 ′ -O Gpp S pG D2 (S P ) are more efficiently recruited to polysomes in HC11 cells and are shifted to heavier polysomes than the analogous mRNAs capped with m 2 7,2 ′ -O GpppG. 31 The structural insight into stabilization of cap/eIF4E complexes by Figure 6. Superposition of eight β-S/β-Se cap-eIF4E structures analyzed in this work. The R P diastereoisomers are colored purple. The S P diastereoisomers are colored cyan, and the sulfur or selenium atoms are represented as orange spheres in all of the analogs. Only m 7 G and triphosphate moieties are shown for clarity. A single protein model is shown as a gray surface representation. The superposition reveals that S or Se atoms always occupy the same space, regardless of the absolute configuration around the phosphorus atom.
ACS Chemical Biology pubs.acs.org/acschemicalbiology Articles phosphorothioate modification presented here provides an explanation of the observed biological properties of clinically relevant β-S-ARCA-capped RNAs and may guide the design of biologically superior mRNA caps that are chemically modified within the phosphodiester linkage between the TSS (+1) and +2 nucleotides. It also reveals two new potentially druggable sites on the eIF4E surface, which might be attractive targets for the inhibition of eukaryotic translation: (i) the N-terminal loop (residues 50−55), occupied by guanine in some of the structures, and (ii) the groove in the vicinity of residues 117− 118, which is probably the path taken by the RNA chain in eIF4E/mRNA complexes.