Kinetic and Thermodynamic Analyses of ... - ACS Publications

Oct 21, 2016 - Department of Basic Medical Sciences, Institute of Medical Science, University of Tokyo, Shirokanedai, Minato-ku, Tokyo 108-8639,. Japa...
2 downloads 0 Views 1MB Size
Subscriber access provided by RYERSON UNIV

Article

Kinetic and thermodynamic analyses of interaction between a high-affinity RNA aptamer and its target protein Ryo Amano, Kenta Takada, Yoichiro Tanaka, Yoshikazu Nakamura, Gota Kawai, Tomoko Kozu, and Taiichi Sakamoto Biochemistry, Just Accepted Manuscript • DOI: 10.1021/acs.biochem.6b00748 • Publication Date (Web): 21 Oct 2016 Downloaded from http://pubs.acs.org on October 22, 2016

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Biochemistry is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 49

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biochemistry

Kinetic and Thermodynamic Analyses of Interaction between a High-Affinity RNA Aptamer and its Target Protein Ryo Amano,† Kenta Takada,† Yoichiro Tanaka,‡ Yoshikazu Nakamura,§,ǁ Gota Kawai,† Tomoko Kozu,# and Taiichi Sakamoto*,†



Department of Life and Environmental Sciences, Faculty of Engineering, Chiba Institute of

Technology, 2-17-1 Tsudanuma, Narashino, Chiba 275-0016, Japan ‡

Facility for RI Research and Education, Instrumental Analysis Center, Yokohama National

University, 79-5 Tokiwadai, Hodogaya-ku, Yokohama 240-8501, Japan §

Department of Basic Medical Sciences, Institute of Medical Science, University of Tokyo,

Shirokanedai, Minato-ku, Tokyo 108-8639, Japan ǁ

Ribomic Inc., 3-16-13 Shirokanedai, Minato-ku, Tokyo 108-0071, Japan

#

Research Institute for Clinical Oncology, Saitama Cancer Center, Ina, Saitama 362-0806, Japan

1 ACS Paragon Plus Environment

Biochemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Keywords SELEX, RNA aptamer, kinetics, thermodynamics

2 ACS Paragon Plus Environment

Page 2 of 49

Page 3 of 49

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biochemistry

ABSTRACT: AML1 (RUNX1) protein is an essential transcription factor involved in the development of hematopoietic cells. Several genetic aberrations that disrupt the function of AML1 have been frequently observed in human leukemia. AML1 contains a DNA-binding domain known as the Runt domain (RD), which recognizes the RD-binding double-stranded DNA element of target genes. In this study, we identified high-affinity RNA aptamers that bind to RD by SELEX. The binding assay using surface plasmon resonance indicated that a shortened aptamer retained the ability to bind to RD when 1 M potassium acetate was used. Thermodynamic study using isothermal titration calorimetry (ITC) showed that the aptamer–RD interaction is driven by a large enthalpy change, and its unfavorable entropy change is compensated by a favorable enthalpy change. Furthermore, binding heat capacity change was identified from the ITC data at various temperatures. The aptamer binding showed a large negative heat capacity change, which suggests that a large apolar surface is buried upon such binding. Thus, we proposed that the aptamer binds to RD with long-range electrostatic force in the early stage of the association and then changes its conformation and recognizes a large surface area of RD. These findings on the biophysics of aptamer binding should be useful for

3 ACS Paragon Plus Environment

Biochemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

understanding the mechanism of RNA–protein interaction and optimizing and modifying RNA aptamers.

4 ACS Paragon Plus Environment

Page 4 of 49

Page 5 of 49

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biochemistry

Aptamers are short single-stranded nucleic acid molecules that are selected in vitro from large random sequence libraries based on their high affinity for target molecules by a process known as Systematic Evolution of Ligands by EXponential enrichment (SELEX).1–5 The target molecules range from small molecules, such as nucleotides, cofactors, and amino acids, followed by peptides, polysaccharides, and proteins, to complex structures, such as whole cells, viruses, and single-celled organisms. Similar to antibodies, aptamers show high affinity and specificity against target molecules as well as additional advantages, including their small size and ease of chemical production.6–8 Thus, aptamers are promising therapeutic agents. However, the details of the biophysical aspects of high-affinity binding of aptamers to target molecules are still unknown. The transcription factor AML1, also known as RUNX1, is one of the key regulators of hematopoiesis and is involved in regulating the transcription of a range of blood cell-specific genes.9–11 The region evolutionarily most conserved in the AML1 protein is the DNA-binding domain, referred to as the Runt domain (RD), which is responsible for binding to a specific DNA sequence.12,13 The tertiary structure of RD has been investigated by X-ray crystallography and nuclear magnetic resonance (NMR) spectroscopy.13–16 These studies revealed that RD recognizes

5 ACS Paragon Plus Environment

Biochemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

the DNA consensus sequence “YGYGGTY” (Y = pyrimidine) using two loop-containing regions and a C-terminal tail, which clamp around the sugar–phosphate backbone. The C-terminal tail and one of the loop regions interact with the major groove. Three guanines of the consensus sequence are recognized by three arginine residues that form hydrogen bonds in the major groove. The other loop region interacts with the minor groove. Complex formation is further supported by several interactions between the side chains and backbone of RD and the sugar– phosphate backbone of DNA. In addition to its role in the regulation of normal hematopoiesis, the AML1 gene is a protooncogene. One of the most common chromosomal translocations in acute myeloid leukemia, t(8;21), generates the oncogenic fusion gene AML1-ETO, which is generally considered to act as a transcriptional repressor by interacting with other proteins.17–29 RNA aptamers that bind to RD were previously studied regarding their potential utility in the diagnosis and treatment of AML1related diseases.30–32 We have already obtained an RNA aptamer (Apt1-S) that shows higher affinity (Kd = 0.99 ± 0.02 nM) than the Runt-binding double-stranded DNA element (RDE, Kd = 9.6 ± 0.2 nM).31,32

6 ACS Paragon Plus Environment

Page 6 of 49

Page 7 of 49

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biochemistry

Nonspecific RNA binding to proteins is assumed to be due mainly to electrostatic interactions between positively charged proteins and the negatively charged phosphate backbone of RNA. Therefore, in this study, to remove the nonspecific binding of RNAs to RD, we performed SELEX under a high salt concentration and isolated RNA aptamers targeting RD with a higher affinity than previously reported aptamers. To reveal the kinetic and thermodynamic aspects of the high affinity of the selected aptamer, surface plasmon resonance (SPR) and isothermal titration calorimetry (ITC) measurements were performed. Our data show that the aptamer binds to a larger surface area of RD than that of RDE.

EXPERIMENTAL PROCEDURES Protein expression and purification. The coding sequence of an AML1 N-terminal fragment (amino acids 51–188, referred to as the Runt domain) was cloned into the pQE80 vector (QIAGEN). Recombinant hexahistidine-tagged AML1–RD was expressed in Escherichia coli strain BL21 (DE3) (Novagen). The crude extract from a 1-L culture was applied to a 5-mL His trap HP column (GE Healthcare) equilibrated with buffer A [20 mM Tris–HCl (pH 8.0), 500 mM

7 ACS Paragon Plus Environment

Biochemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 49

NaCl, and 1 mM DTT] and 5 mM imidazole. A Protease Inhibitor Cocktail Set, EDTA-Free (Merck), was added to E. coli lysate. After washing with buffer A and 5 mM imidazole, the protein was eluted using a 5–500 mM imidazole gradient, followed by dialysis against buffer B [20 mM sodium phosphate (pH 6.5), 50 mM NaCl, and 1 mM DTT]. The sample was then applied to a 5-mL Hi trap SP column (GE Healthcare) equilibrated with buffer B and eluted using a 0–1 M NaCl gradient. The purified protein was then dialyzed against buffer C [20 mM sodium phosphate (pH 6.5), 2 mM magnesium acetate, 300 mM potassium acetate, 50% glycerol, and 1 mM DTT] and stored at −25°C. The concentration of RD was determined based on its molecular absorption coefficient at 280 nm. SELEX.

A

DNA

template

[5′-CTCTCATGTCGGCCGTTA-40N-

CGTCCATTGTGTCCCTATAGTGATCGTATTA-3′; T7 promoter sequence is underlined, and 40N

represents

40

nucleotides

(nt)

of

random

TAATACGACTCACTATAGGGACACAATGGACG-3′),

sequence], and

primer

primer

1 2

(5′(5′-

CTCTCATGTCGGCCGTTA-3′) were synthesized and purchased from Hokkaido System Science Co., Ltd. An initial DNA library was constructed from 400 pmol synthetic DNA

8 ACS Paragon Plus Environment

Page 9 of 49

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biochemistry

template. Thus, a pool of, in theory, 2.4 × 1014 different RNA molecules randomized over 40 nucleotides was prepared by in vitro transcription using T7 RNA polymerase. Before the selection process, the purified RNAs were heat denatured at 95°C for 5 min, refolded by snap cooling to 4°C, and incubated in binding buffer [20 mM sodium phosphate (pH 6.5), 2 mM magnesium acetate, 500 mM potassium acetate, 5% glycerol, 0.05% Triton-X100, and 5 mM βmercaptoethanol] for 10 min at 37°C. To reduce undesired, nonspecific adsorption of RNAs to RD or resins, 2–5 mg/mL tRNA (Invitrogen) was added to the folded RNA pool. Two types of resins, Ni-agarose (QIAGEN) and Talon metal affinity (Clontech) resins, were alternatively used for the affinity capture of the recombinant hexahistidine-tagged RD. Negative selection was performed before and after positive selection. For negative selection, the folded RNA pools were transferred to 3 µL of Ni agarose and Talon metal affinity resins without protein immobilization, followed by incubation for 30 min at room temperature. Then, the supernatant containing unbound RNAs was collected and subjected to positive selection. For positive selection, the RNA pools were transferred to 3 µL of the indicated resin immobilized with RD, followed by incubation for 30 min at room temperature. After the samples had been washed for the indicated

9 ACS Paragon Plus Environment

Biochemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 49

times with 400 µL of binding buffer, tightly bound RNAs were eluted from the resin via the addition of 200 µL of binding buffer containing 500 mM imidazole, subjected to phenol– chloroform extraction, and precipitated using ethanol. The recovered RNAs were reverse transcribed using ReverTra Ace reverse transcriptase (TOYOBO) and amplified by PCR, employing KOD-Plus-DNA polymerase (TOYOBO), with primers 1 and 2, followed by T7 transcription. After eight rounds of SELEX, cDNAs were cloned into the pGEM-T Easy vector (Promega) and sequenced. Details of the conditions for selection in each round are listed in Table S1. Aptamer preparation. S4, S4-S, and Apt1-S were synthesized by in vitro transcription. The template of S4-S was amplified from a cloning vector containing S4 by PCR. PCR products were purified by phenol–chloroform extraction and precipitated using ethanol. The template of Apt1-S was purchased from Hokkaido System Science Co., Ltd. All RNAs were purified by denaturing PAGE and their concentrations were determined based on their molecular absorption coefficients at 260 nm. All RNAs were annealed by heating at 95°C for 5 min, followed by snap cooling on ice.

10 ACS Paragon Plus Environment

Page 11 of 49

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biochemistry

Surface plasmon resonance (SPR) assays. SPR assays were performed as previously described using a BIAcore X instrument (GE Healthcare).31 The 5′-biotinylated dT16 oligomer was immobilized to flow cell 2 on the surface of a streptavidin sensor chip (GE Healthcare) at a level of approximately 800 Resonance Units. The sensor chip was washed and equilibrated in SPR buffer [20 mM sodium phosphate (pH 6.5), 2 mM magnesium acetate, 0.1% Tween 20, and 1 mM DTT] containing 300 mM or 1 M potassium acetate at 298 K. The 3′-A16-tagged RNAs or 3′-A16-tagged RDE (5′-GTCGTTTGCGGTTTGGGGAAAAAAAAAAAAAAAAA-3′ or 5′TCCCCAAACCGCAAACGAC-3′) dissolved in SPR buffer was immobilized to approximately 100 Resonance Units in flow cell 2 at a flow rate of 90 µL/min. RDE and the 40N random RNA pool were used as positive and negative controls, respectively. RD (0.625–2400 nM) or Rev peptide (62.5–2000 nM) in SPR buffer were injected into flow cells 1–2 of the sensor chip for 90 s and dissociated for 180 s. To regenerate the sensor chip, the bound material was completely removed by injecting 4 M urea. The signal of flow cell 2 was subtracted from that of flow cell 1 to eliminate nonspecific interactions. The sensorgrams were analyzed using BIA evaluation software (GE Healthcare). A Langmuir (1:1) binding model was used to analyze the association

11 ACS Paragon Plus Environment

Biochemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 49

rate constant kon and the dissociation rate constant koff (M–1 s−1 and s−1, respectively). The dissociation constant Kd was also determined as the ratio of koff and kon as follows: Kd = koff/kon, and is represented by the mean ± standard error of three independent measurements. Isothermal titration calorimetry (ITC) experiments. ITC experiments were performed at 288, 293, 298, or 303 K with a Microcal iTC200 (Malvern). S4-S, RDE (5′GTCGTTTGCGGTTTGGGGA-3′/5′-TCCCCAAACCGCAAACGAC-3′),

and

RD

were

dissolved in ITC buffer [20 mM sodium phosphate (pH 6.5), 2 mM magnesium acetate, 300 mM potassium acetate, 10% glycerol, and 1 mM DTT]. S4-S or RDE solution in an injection syringe (100 µM) was injected into RD solutions in a cell (10 µM). The volume of each injection was 2 µL, except that of the first injection was 0.4 µL. All measurements involved 19 injections at 150s intervals, a syringe rotational speed of 1250 rpm, and a reference power of 5.0 µcal s−1. To subtract the heat of dilution of the samples, the oligonucleotide samples were injected into the buffer alone. Each area of the heat pulse was integrated, and the resulting values were plotted as a function of the molar ratio. The stoichiometry of binding N, the dissociation constant Kd, and the enthalpy change ∆H were obtained according to single-site-binding model curve fitting provided

12 ACS Paragon Plus Environment

Page 13 of 49

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biochemistry

by the instrument software Origin 7.0 (Microcal), although the c value for ITC was higher than 1000.33-35 The Gibbs free energy change ∆G and the entropy change ∆S were calculated based on Kd and ∆H. Each thermodynamic parameter is represented by the mean ± standard error of three independent measurements. The binding heat capacity change ∆Cp is obtained by measuring the binding enthalpy change as a function of temperature.36 Negative heat capacity change is usually associated with specific recognition and burial of an apolar surface area.

RESULTS Selection and identification of RNA aptamers that bind to RD. We isolated RNA aptamers against RD by SELEX using an initial RNA pool of 73-nt RNA with 40-nt random sequences, referred to as the 40N RNA pool. To select tightly bound aptamers, we performed selection under a high salt concentration and progressively increased the selection stringency by increasing the number of washes and concentration of tRNA as the competitor and decreasing the target concentration for positive selection. Following eight rounds of SELEX, variants of the RNA pool were cloned and sequenced. A total of 36 clones were randomly selected, and 11

13 ACS Paragon Plus Environment

Biochemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 49

independent sequences, S1–S11, were isolated (Table 1). The binding affinity of these 11 sequences to RD was analyzed by SPR, and the dissociation constants (Kd) of all of these clones were lower than 0.3 nM. The Kd of S1, which is most frequently observed, was 0.27 ± 0.02 nM, and that of S4, which shows the highest binding affinity, was 0.044 ± 0.002 nM. Secondary structure prediction and truncation of S4. The secondary structure of S4, which contains a multibranched loop, was predicted using the vs_subopt program (http://www.rna.it-chiba.ac.jp/~vsfold/vs_subopt/)37,38 (Figure 1A). Based on the predicted secondary structure of S4, single-stranded primer binding sequences were truncated and a stem (5′-GGA/CCGA-3′) was added for effective transcription by T7 RNA polymerase and for stabilization of the terminal stem structure, which resulted in a length of 52 nucleotides (Figure 1B). The SPR analysis showed that the Kd of S4-S was 0.034 ± 0.004 nM, which is approximately 300 and 30 times lower than that of RDE (Kd = 10 ± 1 nM) and Apt1-S (Kd = 1.1 ± 0.1 nM), respectively (Figures 1, 2, and Table 2). The difference in the RD-binding affinities between S4-S, RDE, and Ap1-S is mainly attributed to the lower dissociation rate (koff) of S4S/RD interaction than those of the RDE–RD and Apt1-S–RD interactions.

14 ACS Paragon Plus Environment

Page 15 of 49

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biochemistry

To confirm the specificity of S4-S binding to RD, the binding of S4-S, RDE, and Apt1-S to Rev peptide derived from the arginine-rich motif (ARM) of HIV Rev protein as a negative control was examined by SPR (Figure S1 and Table S2). It is well known that the Rev peptide binds to the Rev-response element (RRE) RNA with high affinity (Kd = 6.8 nM at 293 K; ITC experiment) and specificity.39,40 S4-S, RDE, and Apt1-S exhibited low affinity in the micromolar range against Rev peptide because of non-specific binding. It was suggested that S4-S binds to RD with high affinity and specificity. Salt concentration dependence of S4-S binding to RD. The salt concentration dependence of S4-S binding to RD was investigated because S4 was selected under the highly stringent conditions with a high salt concentration. Although the binding affinity of S4-S to RD was decreased as the potassium acetate concentration was increased to 1 M, S4-S retained its binding affinity (Kd = 3.5 ± 0.4 nM) (Table 2 and Figure 2D). On the other hand, the binding of RDE to RD was too weak to be detected using 1 M potassium acetate (Table 2 and Figure 2E). Apt1-S also markedly diminished the binding affinity (Kd = 2800 ± 200 nM) (Table 2 and Figure

15 ACS Paragon Plus Environment

Biochemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 49

2F). The loss of affinity of S4-S at a high salt concentration was because of a decrease in kon without a change in koff. Thermodynamic analysis of S4-S for RD binding. In addition to kinetic analysis using SPR, the thermodynamics of the interactions between S4-S and RD was analyzed by ITC. An exothermic heat pulse was observed after each injection of S4-S or RDE into RD (Figure 3 and Figure S2). Although the Kd values of S4-S involve error due to inaccuracy of curve fitting (c > 1000),33-35 the thermodynamic parameters obtained at 288, 293, 298, and 303 K are shown in Table 3. The stoichiometry N for RD binding to S4-S, RDE, or Apt1-S is almost 1; the difference may depend on the precision of S4-S, RDE, Apt1-S and/or RD concentrations or the proportion of inactive RD in the solution. Comparison of the Kd values indicates that S4-S retains almost the same affinity to RD in the range of 288−303 K, whereas its affinity to RDE declines with increasing temperature. Furthermore, S4-S binds to RD with an approximately 10-fold higher affinity (Kd = 5 ± 3 nM at 298 K) than its affinity to RDE (Kd = 66 ± 5 nM at 298 K). The observed change in enthalpy (∆H) for S4-S binding (∆H = −56 ± 1 kcal mol−1 at 298 K) is greater than that for RDE binding (∆H = −25.5 ± 0.4 kcal mol−1 at 298 K). However, the change in

16 ACS Paragon Plus Environment

Page 17 of 49

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biochemistry

entropy (T∆S) for S4-S binding (−T∆S = −45 ± 1 kcal mol−1 at 298 K) is smaller than that for RDE binding (−T∆S = −15.7 ± 0.4 kcal mol−1 at 298 K). These results indicate that S4-S binding to RD is contributed to by higher enthalpic interaction than its binding to RDE, and its unfavorable T∆S is compensated for by the favorable ∆H. To elucidate the heat capacity change (∆Cp) for RD binding to S4-S and RDE, ∆H values versus temperatures were plotted (Figure 4). The heat capacity change plot revealed that the binding of S4-S or RDE to RD becomes increasingly exothermic with increasing temperature. The ∆Cp values were calculated from the linear fit of ∆H versus temperature data in the 288−303 K range. Both S4-S–RD and RDE–RD interactions exhibit negative heat capacity change in this temperature range, and the ∆Cp value of S4-S–RD interaction (−0.82 ± 0.14 kcal mol−1 K−1) is lower than that of RDE–RD interaction (−0.52 ± 0.05 kcal mol−1 K−1). The T∆S values versus temperatures are plotted for S4-S and RDE, revealing that the binding entropy change for both S4-S and RDE becomes increasingly unfavorable with increasing temperature, akin to the binding enthalpy change (Figure S3).

17 ACS Paragon Plus Environment

Biochemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 49

DISCUSSION

In this work, we obtained high-affinity RNA aptamers against RD under highly stringent conditions (Table 1). Among the selected aptamers, S4 showed the highest binding affinity (Kd = 0.044 ± 0.002 nM), whereas the Kd of the most frequently observed S1 was 0.27 ± 0.02 nM. It was previously reported that the most frequently observed aptamer showed fast association rather than high affinity.41 However, the association of S1 with RD was slower than that of other aptamers, except for S8, in this case. Thus, the reason why the most frequently observed S1 did not show the highest binding affinity could include bias in the efficiency of transcription or reverse-transcription PCR in the SELEX scheme.42 The sequence of S1 may be suitable for amplification in the SELEX scheme. We succeeded in developing S4-S, which is a shortened variant of S4 and shows higher binding affinity than the parent S4. According to the computational prediction of secondary structure, S4-S comprises a multibranched loop, two hairpin stem-loops, and a distal stem. It was suggested that the truncated 3′- and 5′-regions of S4 are dispensable for RD binding and the stabilization of the distal stem structure is favorable for RD binding.

18 ACS Paragon Plus Environment

Page 19 of 49

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biochemistry

It has often been observed that the Kd values obtained by SPR are smaller than those obtained by ITC.43,44 Although the reason for this difference is not clearly understood, it may be caused by the following two differences. First, the solution components are different among the two experiments. The solution used for ITC contained 10% glycerol to prevent RD from precipitation, whereas no glycerol was included in the SPR solution. Thus, glycerol at the interface between S4-S and RD may prevent interaction between them. Second, basic physical conditions are different among the two methods. SPR involves solid-phase measurement of RD binding to S4-S immobilized on a sensor chip, whereas ITC involves solution-phase measurement of RD binding to free S4-S. Therefore, it is assumed that the entropy change upon association of RD with an immobilized S4-S at SPR is smaller than that of RD with a free S4-S in solution at ITC. 45,46 These differences may further lead to a difference in Kd values observed by SPR and ITC. SPR analyses revealed that S4-S binds to RD with faster association and slower dissociation compared with RDE and Apt1-S (Table 2). Although the binding affinity of RDE to RD was completely diminished and that of Apt1-S was extremely reduced when 1 M potassium acetate

19 ACS Paragon Plus Environment

Biochemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 49

was used, S4-S retained its binding activity under the same conditions. Selection under a high salt concentration would result in retention of the binding activity. Furthermore, kinetic analysis of S4-S binding shows that koff did not change and kon decreased by increasing salt concentration, whereas kinetic analysis of Apt1-S binding shows that koff increased and kon decreased. Ions counteract the electrostatic interactions between the negatively charged phosphate backbone of S4-S and the basic residues in RD; this may lead to a decrease in kon.47–50 The result that koff was not changed could be explained by good geometric complementarity between the interacting surfaces S4-S and RD. In some cases of the aptamer that binds to its target protein, the crystal structure revealed that the aptamer binds to its target protein with good shape complementarity and a wide range of contours on protein surfaces.51–54 Thus, we are attempting a structural study of the complex for S4-S and RD to reveal the complementarity of their shape. Figure 5 shows that the partitioning of the contributions of ∆H and T∆S to S4-S/RD interaction and RDE–RD interaction markedly differs. S4-S binding to RD relies on a large favorable ∆H offset by a strongly unfavorable T∆S, whereas RDE binding to RD displays relatively less favorable ∆H and less unfavorable T∆S. The thermodynamic parameters provide

20 ACS Paragon Plus Environment

Page 21 of 49

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biochemistry

information about the interaction between biomolecules.36,40,43–47,55–62 A more favorable ∆H for S4-S–RD interaction suggests more favorable interactions than the RDE–RD interactions; this may be caused by specific binding with good surface complementarity formed by optimal hydrogen bonding and van der Waals interactions. The highly unfavorable T∆S in S4-S binding to RD could arise from decreased solvation entropy and/or decreased conformational entropy in the bound state. Furthermore, the larger negative ∆Cp of S4-S binding than that of RDE binding suggests that a larger apolar surface area is buried upon complex formation for S4-S binding than that for RDE binding (Table 3 and Figure 4). Thus, it was suggested that the unfavorable entropy change is due to conformational restraint in the bound form. These ITC findings are consistent with the fact that the dissociation of S4-S from RD occurs very slowly, as observed by SPR. In summary, SPR and ITC experiments have shown the binding kinetics and thermodynamics of S4-S binding to RD. From the obtained results, it was proposed that in comparison with RDE, S4-S binds to RD with stronger long-range electrostatic force in the early stage of the association and then S4-S changes its conformation and recognizes the larger surface area of RD by optimal hydrogen bonding, van der Waals contact, and/or hydrophobic interaction. Recently, a

21 ACS Paragon Plus Environment

Biochemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 49

combination of structural, kinetic, and thermodynamic analyses has been used for the effective screening and rational optimization of lead chemical compounds and antibody molecules48,61,62. However, optimization and modification of aptamers have been performed by trial and error. Thus, these biophysical analyses may provide some valuable information for the optimization and modification of aptamers. Furthermore, these approaches should be useful for understanding the mechanism underlying various RNA–protein interactions.

22 ACS Paragon Plus Environment

Page 23 of 49

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biochemistry

ASSOCIATED CONTENT Supporting Information

In vitro selection conditions (Table S1), Kinetic parameters of RNA aptamers or RDE binding to Rev peptide (Table S2), SPR analysis of RNA aptamers or RDE binding to Rev peptide (Figure S1), ITC profiles at various temperatures of S4-S/RD interaction and RDE/RD interaction (Figure S2), Temperature dependence of entropy for S4-S binding (Figure S3). These materials are available online free of charge at http://pubs.acs.org.

AUTHOR INFORMATION

Corresponding Author

*Telephone: +81-47-478-0317; Fax: +81-47-478-0317;

Email: [email protected]

Funding

23 ACS Paragon Plus Environment

Biochemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 49

This study was supported by research grants of Innovative Area, Structural Cell Biology (Number 23121528), and the Strategic Research Foundation Grant-aided Project for Private Universities (Number S1101001) from The Ministry of Education, Sports, Culture, Science and Technology (MEXT) of Japan.

Notes

The authors declare no competing financial interests.

ACKNOWLEDGMENT The authors thank Prof. Masato Katahira and Prof. Takashi Nagata for helpful discussions. We also thank Ms. Masako Hirose for technical support with ITC measurements.

ABBREVIATIONS SELEX, systematic evolution of ligands by exponential enrichment; AML1, acute myeloid leukemia 1; RD, Runt domain; RDE, Runt-binding double-stranded DNA element; SPR, surface plasmon resonance; ITC, isothermal titration calorimetry

24 ACS Paragon Plus Environment

Page 25 of 49

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biochemistry

REFERENCES (1) Ellington, A. D., and Szostak, J. (1990) In vitro selection of RNA molecules that bind specific ligands. Nature 346, 818–822. (2) Robertson, D. L., and Joyce, G. F. (1990) Selection in vitro of an RNA enzyme that specifically cleaves single-stranded-DNA. Nature 344, 467–468. (3) Tuerk, C., and Gold, L. (1990) Systematic evolution of ligands by exponential enrichment – RNA ligands to bacteriophage-T4 DNA-polymerase. Science 249, 505–510. (4) Stoltenburg, R., Reinemann, C., and Strehlitz, B. (2007) SELEX-A (r)evolutionary method to generate high-affinity nucleic acid ligands. Biomol Eng. 24, 381–403. (5) Szeitner, Z., András, J., Gyurcsányi, R. E., and Mészáros, T. (2014) Is less more? Lessons from aptamer selection strategies. J Pharm Biomed Anal. 101, 58–65. (6) Klussmann, S. (2006) The Aptamer Handbook: Functional Oligonucleotides and Their Applications, Wiley-VCH, Weinheim, Germany.

25 ACS Paragon Plus Environment

Biochemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(7) Que-Gewirth, N. S., and Sullenger, B. A. (2007) Gene therapy progress and prospects: RNA aptamaer. Gene Therapy 14, 283–291. (8) Keefe, A. D., Pai, S., and Ellington, A. D. (2010) Aptamers as therapeutics. Drug discovery 9, 537–548. (9) Mikhail, F. M., Sinha, K. K., Saunthararajah, Y., and Nucifora, G. (2006) Normal and transforming functions of RUNX1: a perspective. J. Cell Physiol. 207, 582–593. (10) Speck, N. A. (2001) Core binding factor and its role in normal hematopoietic development. Curr. Opin. Hematol. 8, 192–196. (11) de Bruijn, M. F., and Speck, N. A. (2004) Core-binding factors in hematopoiesis and immune function. Oncogene. 23, 4238–4248. (12) Collins, A., Littman, D. R., and Taniuchi, I. (2009) RUNX proteins in transcription factor networks that regulate T-cell lineage choice. Nat Rev Immunol 9, 106–115. (13) Bravo, J., Li, Z., Speck,N.A., and Warren, A. J. (2001) The leukemia-associated AML1 (Runx1)-CBF beta complex functions as a DNA-induced molecular clamp. Nat. Struct. Biol. 8, 371–378.

26 ACS Paragon Plus Environment

Page 26 of 49

Page 27 of 49

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biochemistry

(14) Nagata, T., Gupta, V., Sorce, D., Kim, W.Y., Sali, A., Chait, B.T., Shigesada, K., Ito, Y., and Werner, M.H. (1999) Immunoglobulin motif DNA recognition and heterodimerization of the PEBP2/CBF Runt domain. Nat. Struct. Biol. 6, 615–619. (15) Tahirov, T.H., Inoue-Bungo, T., Morii, H., Fujikawa, A., Sasaki, M., Kimura, K., Shiina, M., Sato, K., Kumasaka, T., Yamamoto, M., Ishii, S., and Ogata, K. (2001) Structural analyses of DNA recognition by the AML1/Runx-1 Runt domain and its allosteric control by CBFbeta. Cell 104, 755–767. (16) Bäckström, S., Wolf-Watz, M., Grundström, C., Härd, T., Grundström, T., and Sauer, U. H. (2002) The RUNX1 Runt domain at 1.25A resolution: a structural switch and specifically bound chloride ions modulate DNA binding. J. Mol. Biol. 322, 259-272. (17) Gelmetti, V., Zhang, J., Fanelli, M., Minucci, S., Pelicci, P.G., and Lazar, M. A. (1998) Aberrant recruitment of the nuclear receptor corepressor-histone deacetylase complex by the acute myeloid leukemia fusion partner ETO. Mol. Cell. Biol. 18, 7185–7191. (18) Lutterbach, B., Westendorf, J. J., Linggi, B., Patten, A., Moniwa, M., Davie, J. R., Huynh, K. D., Bardwell, V. J., Lavinsky, R. M., Rosenfeld, M. G., Glass, C., Seto, E., and Hiebert, S. W.

27 ACS Paragon Plus Environment

Biochemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(1998) ETO, a target of t(8;21) in acute leukemia, interacts with the N-CoR and mSin3 corepressors. Mol. Cell. Biol. 18, 7176–7184. (19) Wang, J., Hoshino, T., Redner, R. L., Kajigaya, S., and Liu, J. M. (1998) ETO, fusion partner in t(8;21) acute myeloid leukemia, represses transcription by interaction with the human N-CoR/mSin3/HDAC1 complex. Proc. Natl Acad. Sci. USA 95, 10860–10865. (20) Westendorf, J. J., Yamamoto, C. M., Lenny, N., Downing, J. R., Selsted, M. E., and Hiebert, S. W. (1998) The t(8;21) fusion product, AML-1–ETO, associates with C/EBP-α, inhibits C/EBP-α-dependent transcription, and blocks granulocytic differentiation. Mol. Cell. Biol. 18, 322–333. (21) Mao, S., Frank, R. C., Zhang, J., Miyazaki, Y., and Nimer, S. D. (1999) Functional and physical interactions between AML1 proteins and an ETS protein,MEF: implications for the pathogenesis of t(8;21)-positive leukemias. Mol. Cell. Biol. 19, 3635–3644. (22) Elagib, K. E., Racke, F. K., Mogass, M., Khetawat, R., Delehanty, L. L., and Goldfarb, A. N. (2003) RUNX1 and GATA-1 coexpression and cooperation in megakaryocytic differentiation, Blood 101. 4333–4341.

28 ACS Paragon Plus Environment

Page 28 of 49

Page 29 of 49

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biochemistry

(23) Zhang, J., Kalkum, M., Yamamura, S., Chait, B. T., and Roeder, R. G. (2004) E protein silencing by the leukemogenic AML1-ETO fusion protein. Science 305, 1286–1289. (24) Gardini, A., Cesaroni, M., Luzi, L., Okumura, A. J., Biggs, J. R., Minardi, S. P., Venturini, E., Zhang, D. E., Pelicci, P. G., and Alcalay, M. (2008) AML1/ETO oncoprotein is directed to AML1 binding regions and co-localizes with AML1 and HEB on its targets. PLoS Genet. 4, e1000275. (25) Guo, C., Hu, Q., Yan, C., and Zhang, J. (2009) Multivalent binding of the ETO corepressor to E proteins facilitates dual repression controls targeting chromatin and the basal transcription machinery. Mol. Cell. Biol. 29, 2644–2657. (26) Martens, J. H., Mandoli, A., Simmer, F., Wierenga, B. J., Saeed, S., Singh, A. A., Altucci, L., Vellenga, E., and Stunnenberg, H. G. (2012) ERG and FLI1 binding sites demarcate targets for aberrant epigenetic regulation by AML1-ETO in acute myeloid leukemia. Blood 120, 4038– 4048. (27) Wang, L., Gural, A., Sun, X. J., Zhao, X., Perna, F., Huang, G., Hatlen, M. A., Vu, L., Liu, F., Xu, H., Asai, T., Xu, H., Deblasio, T., Menendez, S., Voza, F., Jiang, Y., Cole, P. A., Zhang,

29 ACS Paragon Plus Environment

Biochemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

J., Melnick, A., Roeder, R. G., and Nimer, S. D. (2011) The leukemogenicity of AML1-ETO is dependent on site-specific lysine acetylation. Science 333, 765–769. (28) Shia, W. J., Okumura, A. J., Yan, M., Sarkeshik, A., Lo, M. C., Matsuura, S., Komeno, Y., Zhao, X., Nimer, S. D., Yates, J. R. 3rd., and Zhang, D. E. (2012) PRMT1 interacts with AML1ETO to promote its transcriptional activation and progenitor cell proliferative potential. Blood 119, 4953–4962. (29) Sun, X., Wang, Z., Wang, L., Jiang, Y., Kost, N., Soong, T. D., Chen, W., Tang, Z., Nakadai, T., Elemento, O., Fischle, W., Melnick, A., Patel, D. J., Nimer, S. D., and Roeder, R. G. (2013) A stable transcription factor complex nucleated by oligomeric AML1–ETO controls leukaemogenesis. Nature 500, 93–97. (30) Barton, J. L., Bunka, D. H. J., Knowling, S. E., Lefevre, P., Warren, A. J., Bonifer, C., and Stockley, P. G. (2009) Characterization of RNA aptamers that disrupt the RUNX1CBFbeta/DNA complex. Nucleic Acids Res 37, 6818–6830.

30 ACS Paragon Plus Environment

Page 30 of 49

Page 31 of 49

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biochemistry

(31) Fukunaga, J., Tanaka, Y., Nomura, Y., Amano, R., Tanaka, T., Nakamura, Y., Kawai, G., Sakamoto, T., and Kozu, T. (2013) The Runt domain of AML1 (RUNX1) binds a sequenceconserved RNA motif that mimics a DNA element. RNA 19, 927–936. (32) Nomura, Y., Tanaka, Y., Fukunaga, J., Fujiwara, K., Chiba, M., Iibuchi, H., Tanaka, T., Nakamura, Y., Kawai, G., Kozu, T., and Sakamoto, T. (2013) Solution structure of a DNA mimicking motif of an RNA aptamer against transcription factor AML1 Runt domain. J. Biochem. 154, 513–519. (33) Krainer, G., and Keller, S. (2015) Single-experiment displacement assay for quantifying high-affinity binding by isothermal titration calorimetry. Methods 76, 116–123. (34) Tellinghuisen, J. (2005) Optimizing Experimental Parameters in Isothermal Titration Calorimetry. J Phys Chem B. 109, 20027–20035. (35) Velázquez, C. A., and Freire, E. (2005) ITC in the post-genomic era...? Priceless. Biophys Chem. 115, 115–124.

31 ACS Paragon Plus Environment

Biochemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 32 of 49

(36) Sokolosk, J. E., Dombrowski, S. E., and Bevilacqua, P. C. (2012) Thermodynamics of Ligand Binding to a Heterogeneous RNA Population in the Malachite Green Aptamer. Biochemistry 51, 565–572. (37) Dawson, W., Fujiwara, K., Futamura, Y., Yamamoto, K., and Kawai, G. (2006) A method for finding optimal RNA secondary structures using a new entropy model (vsfold). Nucleosides Nucleotides Nucl. Acids 25, 171–189. (38) Dawson, W., Takai, T., Ito, I., Shimizu, K., and Kawai, G. (2014) A new entropy model for RNA: part III, Is the folding free energy landscape of RNA funnel shaped? Journal of Nucleic Acids Investigation 5, 2652. (39) Battiste, J. L., Mao, H., Rao, N. S., Tan, R., Muhandiram, D. R., Kay, L. E., Frankel, A. D., and Williamson, J. R. (1996) Alpha helix-RNA major groove recognition in an HIV-1 rev peptide-RRE RNA complex. Science 273, 1547–1551. (40) Jayaraman, B., Mavor, D., Gross, J. D., and Frankel, A. D. (2015) Thermodynamics of Rev−RNA Interactions in HIV‑1 Rev−RRE Assembly. Biochemistry 54, 6545–6554.

32 ACS Paragon Plus Environment

Page 33 of 49

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biochemistry

(41) Ruigrok, V. J., van Duijn, E., Barendregt, A., Dyer, K., Tainer, J. A., Stoltenburg, R., Strehlitz, B., Levisson, M., Smidt, H., and van der Oost, J. (2012) Kinetic and stoichiometric characterisation of streptavidin-binding aptamers. Chembiochem.13, 829–836. (42) Tsuji, S., Hirabayashi, N., Kato, S., Akitomi, J., Egashira, H., Tanaka, T., Waga, I., and Ohtsu, T. (2009) Effective isolation of RNA aptamer through suppression of PCR bias. Biochem Biophys Res Commun. 386, 223–226.

(43) Kamatari, Y. O., Ohta, S., Inoshima, Y., Oda, M., Maruno, T., Kobayashi, Y., and Ishiguro, N. (2014), Identification and characterization of a multispecific monoclonal antibody G2 against chicken prion protein. Protein Sci. 23, 1050–1059 (44) Potty, A. S., Kourentzi, K., Fang, H., Jackson, G. W., Zhang, X., Legge, G. B., and Willson R. C. (2009) Biophysical characterization of DNA aptamer interactions with vascular endothelial factor. Biopolymers 91, 145–156. (45) Mounier, C. M., Hackeng, T. M., Schaeffer, F., Faure, G., Bon, C., and Griffin, J. H. (1998) Inhibition of prothrombinase by human secretory phospholipase A2 involves binding to factor Xa.

33 ACS Paragon Plus Environment

Biochemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 34 of 49

J Biol Chem 273, 23764–23772. (46) Mouratou, B., Schaeffer, F., Guilvout, I., Tello-Manigne, D., Pugsley, A. P., Alzari, P. M., and Pecorari, F. (2007) Remodeling a DNA-binding protein as a specific in vivo inhibitor of bacterial secretin PulD. Proc Natl Acad Sci U S A. 104, 17983–17988. (47) Oda, M., and Nakamura, H. (2000) Thermodynamic and kinetic analyses for understanding sequence-specific DNA recognition. Genes to Cells 5, 319–326. (48) Fukunaga, A., and Tsumoto, K. (2013) Improving the affinity of an antibody for its antigen via long-range electrostatic interactions. Protein Eng. Des. Sel. 26, 773–780. (49) Tang, Q., Su, X., and Loh, K. P. (2007) Surface plasmon resonance spectroscopy study of interfacial binding of thrombin to antithrombin DNA aptamers. J Colloid Interface Sci. 315, 99– 106. (50) Su, X., Lin, C. Y., O'Shea, S. J., Teh, H. F., Peh, W. Y., and Thomsen, J. S. (2006) Combinational application of surface plasmon resonance spectroscopy and quartz crystal microbalance for studying nuclear hormone receptor-response element interactions. Anal Chem. 78, 5552–5558.

34 ACS Paragon Plus Environment

Page 35 of 49

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biochemistry

(51) Miyakawa, S., Nomura, Y., Sakamoto, T., Yamaguchi, Y., Kato, K., Yamazaki, Y., and Nakamura, Y. (2008) Structural and molecular basis for hyperspecificity of RNA aptamer to human immunoglobulin G. RNA 14, 1154–1163. (52) Nomura, Y., Sugiyama, S., Sakamoto, T., Miyakawa, S., Adachi, H., Takano, K., Murakami, S., Inoue, T., Mori, Y., Nakamura, Y., and Matsumura, H. (2010) Conformational plasticity of RNA for target recognition as revealed by the 2.15 A crystal structure of a human IgG-aptamer complex. Nucleic Acids Res. 38, 7791–7799. (53) Kato, K., Ikeda, H., Miyakawa, S., Futakawa, S., Nonaka, Y., Fujiwara, M., Okudaira, S., Kano, K., Aoki, J., Morita, J., Ishitani, R., Nishimasu, H., Nakamura, Y., and Nureki, O. (2016) Structural basis for specific inhibition of Autotaxin by a DNA aptamer. Nat Struct Mol Biol. 23, 395–401. (54) Gelinas, A. D., Davies, D. R., and Janjic, N. (2016) Embracing proteins: structural themes in aptamer-protein complexes. Curr Opin Struct Biol. 36, 122–132.

35 ACS Paragon Plus Environment

Biochemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 36 of 49

(55) Freisz, S., Mezher, J., Hafirassou, L., Wolff, P., Nominé, Y., Romier, C., Dumas, P., and Ennifar, E. (2012) Sequence and structure requirements for specific recognition of HIV-1 TAR and DIS RNA by the HIV-1 Vif protein. RNA Biol. 9, 966–977. (56) Bec, G., Meyer, B., Gerard, M. A., Steger, J., Fauster, K., Wolff, P., Burnouf, D., Micura, R., Dumas, P., and Ennifar, E. (2013) Thermodynamics of HIV-1 reverse transcriptase in action elucidates the mechanism of action of non-nucleoside inhibitors. J Am Chem Soc. 135, 9743– 9752. (57) Nagatoishi, S., Tanaka, Y., Kudou, M., and Tsumoto, K. (2009) The interaction of hyperthermophilic TATA-box binding protein with single-stranded DNA is entropically favorable and exhibits a large negative heat capacity change at high salt concentration. Mol Biosyst. 5, 957–961. (58) Da Costa, J. B., Andreiev, A. I., and Dieckmann, T. (2013) Thermodynamics and kinetics of adaptive binding in the malachite green RNA aptamer. Biochemistry 52, 6575–6583. (59) Lin, P. H., Chen, R. H., Lee, C. H., Chang, Y., Chen, C. S., and Chen, W. Y. (2011) Studies of the binding mechanism between aptamers and thrombin by circular dichroism, surface

36 ACS Paragon Plus Environment

Page 37 of 49

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biochemistry

plasmon resonance and isothermal titration calorimetry. Colloids Surf B Biointerfaces. 88, 552– 558. (60) Potty, A. S., Kourentzi, K., Fang, H., Schuck, P., and Willson, R. C. (2011) Biophysical characterization of DNA and RNA aptamer interactions with hen egg lysozyme. Int J Biol Macromol. 48, 392–397. (61) Freire, E. (2009) A thermodynamic approach to the affinity optimization of drug candidates. Chem Biol Drug Des. 74, 468–472. (62) Ennifar, E., Aslam, M. W., Strasser, P., Hoffmann, G., Dumas, P., and van Delft, F. L. (2013) Structure-guided discovery of a novel aminoglycoside conjugate targeting HIV-1 RNA viral genome. ACS Chem Biol. 8, 2509–2517.

37 ACS Paragon Plus Environment

Biochemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48

Page 38 of 49

Table 1. Sequences of RD-binding RNA aptamers and their kinetic parameters a Clone ID S1

S2

S3

S4

S5

S6

S7

S8

S9

Sequence (5′ to 3′) b

kon

koff

(×106 M−1s−1)

(×10−3 s−1)

10/36

8.5 ± 0.4

2.3 ± 0.1

0.27 ± 0.02

8/36

9.6 ± 0.2

1.1 ± 0.1

0.11 ± 0.01

6/36

9.7 ± 0.3

2.1 ± 0.1

0.21 ± 0.01

3/36

9.1 ± 0.1

0.40 ± 0.02

0.044 ± 0.002

2/36

9.3 ± 0.4

1.1 ± 0.1

0.11 ± 0.13

2/36

13 ± 2

1.5 ± 0.1

0.12 ± 0.09

1/36

11.6 ± 0.2

0.7 ± 0.3

0.056 ± 0.003

1/36

6.8 ± 0.3

0.53 ± 0.8

0.077 ± 0.004

1/36

13.1 ± 0.3

0.82 ± 0.2

0.063 ± 0.002

Frequency

gggacacaauggacgUGUCGGCCCUGCCGUGUAACGC UGGCGCGGGAUGUUCUCCuaacggccgacaugagag gggacacaauggacgGCCCAGCCACCUAGAGCGAGCG CGCAAUGGAGACCCAUUGuaacggccgacaugagag gggacacaauggacgUGUCGGCCCUGCCGUGUAACGC UGGCGCGGGACUUCUCCuaacggccgacaugagag gggacacaauggacgGCCCUGCCACGAUAGCGGCGCG GGAAGUAAAGUAUACACCuaacggccgacaugagag gggacacaauggacgUGUCGGCCCUGCCGUGUAAUGC UGGCGCGGGACGUUCUCCuaacggccgacaugagag gggacacaauggacgGUCAGCCACCACUGUGCGGCGA GCGGAAGCACACCGUCCGuaacggccgacaugagag gggacacaauggacgGCCCUGCCACCUAGAGCGAGCG CGCAAUGGAGACCCAUUGuaacggccgacaugagag gggacacaauggacgGCCCUGCCACGAAGGCGGCGCG CAGGCUACCCGCACCUGuaacggccgacaugagag gggacacaauggacgAUGCCGGCCCUGCCACACCAAU GCGGCGCGGUCAAUAGACuaacggccgacaugagag

38 ACS Paragon Plus Environment

Kd (nM)

Page 39 of 49

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48

Biochemistry

S10

S11

a

gggacacaauggacgGCCCAGCCACCUAGUGCGAGCG CGCAAUGGAGACCCAUUGuaacggccgacaugagag gggacacaauggacgACGCCGGCCCUGCCACACCAAU GCGGCGCGGUCAAUAGACuaacggccgacaugagag

1/36

10.7 ± 0.5

0.97 ± 0.4

0.090 ± 0.001

1/36

13.7 ± 0.2

1.27 ± 0.03

0.093 ± 0.002

All SPR experiments were performed in 20 mM sodium phosphate (pH 6.5) with 2 mM magnesium acetate, 300 mM

potassium acetate, 0.1% Tween 20, and 1 mM DTT at 298 K. The kon, koff, and Kd values represent the mean ± SE (n = 3).

b

Lower-case letters indicate nucleotides derived from the primer-binding sites.

39 ACS Paragon Plus Environment

Biochemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 40 of 49

Table 2. Kinetics of the binding of S4-S, RDE, or Apt1-S to RD at 0.3 M or 1 M potassium acetate a

0.3 M potassium acetate

1 M potassium acetate

kon

koff

Kd

kon

koff

Kd

(×106 M−1s−1)

(×10−3 s−1)

(nM)

(×106 M−1s−1)

(×10−3 s−1)

(nM)

S4-S

10.7 ± 0.3

0.37 ± 0.06

0.034 ± 0.004

0.13 ± 0.02

0.47 ± 0.02

3.5 ± 0.4

RDE

6.4 ± 0.8

65 ± 3

10 ± 1

Apt1-S

6.7 ± 0.2

7.6 ± 0.6

1.1 ± 0.1

a

N.D.b 0.10 ± 0.01

280 ± 10

2800 ± 200

All experiments were performed in 20 mM sodium phosphate (pH 6.5) with 2 mM magnesium

acetate, 300 mM potassium acetate, or 1 M potassium acetate, 0.1% Tween 20, and 1 mM DTT at 298 K. The kon, koff, and Kd values represent the mean ± SE (n = 3).

b

Not determined for no significant binding.

40 ACS Paragon Plus Environment

Page 41 of 49

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48

Biochemistry

Table 3. Thermodynamics of the binding of S4-S or RDE to RD at 0.3 M potassium acetate a

S4-S

RDE

a

Temperature (K)

N

Kd (nM)

∆G

∆H

(kcal mol−1)

288

0.96 ± 0.03

3± 2

293

0.89 ± 0.03

298

(kcal mol−1)

(kcal mol−1 K−1)

−11.2 ± 0.3

−45.4 ± 0.2

−34.2 ± 0.2

−0.82 ± 0.14

5±3

−11.2 ± 0.4

−51.4 ± 0.4

−40.2 ± 0.8

0.87 ± 0.02

5±3

−11.4 ± 0.4

−56 ± 1

−45 ± 1

303

0.90 ± 0.03

6±3

−11.4 ± 0.3

−57.5 ± 0.8

−46 ± 1

288

1.02 ± 0.04

20 ± 4

−10.2 ± 0.1

−21.1± 0.4

−11.0 ± 0.3

293

0.93 ± 0.02

28 ± 2

−10.13 ± 0.03

−23.69 ± 0.02

−13.56± 0.04

298

0.99 ± 0.08

66 ± 5

−9.82 ± 0.05

−25.5 ± 0.4

−15.7 ± 0.4

303

0.87 ± 0.08

130 ± 30

−9.6 ± 0.2

−29 ± 1

−20 ± 1

−0.52 ± 0.05

All experiments were performed in 20 mM sodium phosphate (pH 6.5) with 2 mM magnesium acetate, 300 mM potassium

acetate, 10% glycerol, and 1 mM DTT. The N, Kd, ∆G, ∆H, and T∆S values represent the mean ± SE (n = 3). b

∆Cpb

T∆S (kcal mol−1)

The ∆Cp value was obtained between 288 K and 303 K.

ACS Paragon Plus Environment

Biochemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 42 of 49

Figure Legends Figure 1. The secondary structures of S4, S4-S, RDE, and Apt1-S. The secondary structures of S4 (A) and S4-S (B) were predicted by vs_subopt. Lower-case letters indicate nucleotides derived from the primer-binding sites. The secondary structures of RDE (C) and Apt1-S (D). The important bases for RD binding are highlighted in bold.

Figure 2. Effect of salt concentration on S4-S binding to RD. SPR sensorgrams of S4-S/RD interaction at RD concentrations of 0.625, 1.25, 2.5, 5, 10, and 20 nM with 0.3 M potassium acetate (A) and at RD concentrations of 37.5, 75, 150, 300, 600, and 1200 nM with 1 M potassium acetate (D). SPR sensorgrams of RDE/RD interaction at RD concentrations of 2.5, 5, 10, 20, 40, and 80 nM with 0.3 M potassium acetate (B) and at an RD concentration of 1000 nM with 1 M potassium acetate (E). SPR sensorgrams of Apt1-S/RD interaction at RD concentrations of 2.5, 5, 10, 20, 30, and 40 nM with 0.3 M potassium acetate (C) and at RD concentrations of 150, 300, 600, 900, 1200, and 2400 nM with 1 M potassium acetate (F).

Figure 3. ITC data for S4-S/RD interaction and RDE/RD interaction. Titrations of RD with S4-S (A) and with RDE (B) were carried out at 298 K. Top panels show ITC traces and bottom panels show integrated heat values. Data were fitted using the one set of sites binding model.

Figure 4. Temperature dependence of enthalpy for S4-S binding and RDE binding to RD. ∆H values a given temperature for S4-S (open) and RDE (filled). Error bars show the standard error of the mean. These plots provide ∆Cp values for S4-S/RD and RDE/RD interactions.

ACS Paragon Plus Environment

42

Page 43 of 49

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biochemistry

Figure 5. Graphical view of enthalpy−entropy partition. Comparison of thermodynamic parameters for S4-S/RD interaction (gray bars) and RDE/RD interaction (black bars) at 0.3 M potassium acetate at 298 K indicates different partitioning of free energy into enthalpy and entropy terms.

ACS Paragon Plus Environment

43

Biochemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 44 of 49

For Table of Contents Use Only Title “Kinetic and thermodynamic analyses of interaction between a high-affinity RNA aptamer and its target protein” Authors Ryo Amano, Kenta Takada, Yoichiro Tanaka, Yoshikazu Nakamura, Gota Kawai, Tomoko Kozu, and Taiichi Sakamoto

ACS Paragon Plus Environment

44

Page 45 of 49

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biochemistry

Figure 1. The secondary structures of S4, S4-S, RDE, and Apt1-S. The secondary structures of S4 (A) and S4-S (B) were predicted by vs_subopt. Lower-case letters indicate nucleotides derived from the primerbinding sites. The secondary structures of RDE (C) and Apt1-S (D). The important bases for RD binding are highlighted in bold. 294x414mm (300 x 300 DPI)

ACS Paragon Plus Environment

Biochemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 2. Effect of salt concentration on S4-S binding to RD. SPR sensorgrams of S4-S/RD interaction at RD concentrations of 0.625, 1.25, 2.5, 5, 10, and 20 nM with 0.3 M potassium acetate (A) and at RD concentrations of 37.5, 75, 150, 300, 600, and 1200 nM with 1 M potassium acetate (D). SPR sensorgrams of RDE/RD interaction at RD concentrations of 2.5, 5, 10, 20, 40, and 80 nM with 0.3 M potassium acetate (B) and at an RD concentration of 1000 nM with 1 M potassium acetate (E). SPR sensorgrams of Apt1-S/RD interaction at RD concentrations of 2.5, 5, 10, 20, 30, and 40 nM with 0.3 M potassium acetate (C) and at RD concentrations of 150, 300, 600, 900, 1200, and 2400 nM with 1 M potassium acetate (F). 412x373mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 46 of 49

Page 47 of 49

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biochemistry

Figure 3. ITC data for S4-S/RD interaction and RDE/RD interaction. Titrations of RD with S4-S (A) and with RDE (B) were carried out at 298 K. Top panels show ITC traces and bottom panels show integrated heat values. Data were fitted using the one set of sites binding model. 367x279mm (300 x 300 DPI)

ACS Paragon Plus Environment

Biochemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 4. Temperature dependence of enthalpy for S4-S binding and RDE binding to RD. ∆H values a given temperature for S4-S (open) and RDE (filled). Error bars show the standard error of the mean. These plots provide ∆Cp values for S4-S/RD and RDE/RD interactions. 294x216mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 48 of 49

Page 49 of 49

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biochemistry

Figure 5. Graphical view of enthalpy−entropy partition. Comparison of thermodynamic parameters for S4S/RD interaction (gray bars) and RDE/RD interaction (black bars) at 0.3 M potassium acetate at 298 K indicates different partitioning of free energy into enthalpy and entropy terms. 311x213mm (300 x 300 DPI)

ACS Paragon Plus Environment