Direct observation of small molecule activator binding to single PR65 protein

Direct observation of small molecule activator binding to single PR65 protein

Introduction

PR65, the scaffolding subunit of the serine/threonine protein phosphatase 2A (PP2A) family, provides a platform for the assembly of this heterotrimeric complex1. PP2A plays a critical role in regulating cell proliferation, maintaining cellular homeostasis and orchestrating cellular signaling pathways2,3. Dysregulation of PP2A has been linked to various cancers and chronic conditions such as Alzheimer’s disease and chronic obstructive pulmonary disease, thus presenting itself as a target for therapeutic intervention4,5,6,7,8. Reactivation of PP2A through a small-molecule activator of PP2A (SMAP), such as ATUX-8385, has shown to decrease hepatoblastoma and neuroblastoma cell proliferation in vitro and in animal models9. The aforementioned SMAPs are of the tricyclic sulfonamide class and bind directly to PR65, the subunit A of the trimeric PP2A. PR65 itself comprises of 15 α-helical HEAT repeats. The HEAT repeats of PR65 form an elongated curved structure, with one helix from each repeating forming the outer, convex surface, and the other forming the inner, concave surface onto which the regulatory (B) and catalytic (C) subunits assemble. These SMAPs have been reported to drive conformation changes that subsequently promote heterotrimer formation10. Mechanistically, the ligand-induced conformation change is analogous to mechanical extension or compression10,11. The structural impacts of ATUX-8385 on PR65 upon binding allow for a deeper understanding of the activation mechanism. Furthermore, for safety, efficacy and duration of response, it is important to understand the binding kinetics of any small molecule drug candidate.

While there have been a variety of techniques employed to study the structural changes and binding kinetics of small molecules with proteins, here we focus on nanoaperture optical tweezers (NOTs). NOTs allow for studying single proteins and their interactions without the need for modifications, such as labels or tethers12,13,14. Here we use a double-nanohole (DNH) NOT for generating a large gradient force capable of holding onto nanoparticles15,16. More generally, shaped nanoapertures have been used in the optical tweezer field for the study of small biomolecules and other nanoparticles16,17,18,19,20,21,22,23,24,25,26,27. The first report of trapping of a single protein demonstrated the technique’s sensitivity by observing the conformational changes in bovine serum albumin (BSA) via changes in light transmission intensity as measured by an avalanche photodiode (APD)12. Integration of NOT with a nanopore further supports single protein trapping25,28,29 and has been used to probe the interaction of antibodies with immunotherapy-relevant ligands30.

Conformational changes were suggested from the NOT signal changes for two different types of proteins: βamylase and heat-shock protein 90 (HSP90)31. Iron ion loading and unloading of a single ferritin was observed by seeing structural changes in real-time using this technique32. Recently, we have used NOTs to study the impact of point mutation on the conformational behavior of PR6533, showing correlation between the extension of the protein observed in the optical scattering signal of the NOT and those predicted from molecular dynamics (MD) simulations.

The goal of the present study is to use NOTs to probe the interaction of a single PR65 protein with the clinically-relevant experimental small molecule ATUX-8385. We performed binding assays under equilibrium conditions and observed biomolecular dynamics in real-time. By tracking ligand-induced conformational changes and subsequently discerning the optical scattering signal into binary states, bound and unbound, dissociation and association rate constants can be determined through the “bound” time and “unbound” time. This is similar to single-molecule Förster resonance energy transfer (FRET), however, the ligand concentration is not limited to the typical nanomolar range as we are not hindered by background fluorescence noise by looking directly at the polarizability change of a single protein34.

Previously, NOTs have been used to quantify established protein-drug interactions inhuman serum albumin with tolbutamide and with phenytoin, obtaining dissociation constants within the previously reported range of values13. NOTs were also used to study strong interactions where the binding was essentially irreversible on the timescale of observation, e.g., biotin with streptavidin and acetylsalicylic acid with cyclooxygenase 235. Here we apply and extend the technique to a much less understood binding pair, to understand the binding kinetics and to observe the conformational changes induced upon binding.

In addition to NOTs, we use several techniques to investigate the binding characteristics and effects of SMAP ATUX-8385 to PR65. These include all-atom MD simulations, nano differential scanning fluorimetry (nanoDSF), nuclear magnetic resonance (NMR) and fluorescence polarization anisotropy. In terms of molecular simulations, we explored the effect of ATUX-8385 binding on the structure and dynamics of PR65 by docking simulations followed by microseconds long MD runs, as well as microseconds simulations for apo PR65. Our docking simulations revealed a binding region that encompasses the inner helices of the HEAT repeats 4 and 5 (4i and 5i) and the outer helices of HEAT repeats 5 and 6 (5o and 6o). ATUX-8385 showed stable binding at this site throughout the 704 ns long triplicate MD runs of MD simulation, totaling 2.112 μs in combined trajectory length. ATUX-8385 binding stabilized an extended conformation in two runs and a compact form in a third run. The latter transitioned to the extended form towards the end of the 704 ns simulation, indicating the stabilizing effect of ATUX-8385 bound to PR65, and its tendency to favor the extended conformation. NanoDSF showed that the presence of ATUX-8385 stabilized PR65 through increase in the melting temperature. NMR confirmed the binding of ATUX-8385. Fluorescence polarization anisotropy experiments provided an ensemble dissociation constant in close agreement with the single molecule dissociation constant.

Results

Structural impact of ATUX-8385 binding: NOT experimental results

A home-built optical trapping setup and DNHs fabricated colloidally were used for the experiments as shown in Fig. 1. The trapping event was characterized as the change in transmission intensity due to dielectric loading and increase in noise due to Brownian motion of the protein in the trap16,36,37, see Fig. 1c. A single jump in transmission level, as depicted in Fig. 1c, typically represents an individual protein trapping event unless a protein exists naturally as a homo-oligomer. In the absence of other proteins, PR65 exists as a monomer and the trapping events observed were of single proteins. By contrast, multiple protein trapping events show up clearly in the signal, as has been reported in the literature for BSA and egg white protein38,39. For our study, we focused on the difference between the trapped signal (transmission intensity picked up by the APD after the trapping event) of apo-PR65 in 5% DMSO and PR65 with ATUX-8385 in 5% DMSO. A total of 8 measurements were collected at 20 μM of ligand and a protein concentration of 10 μM.

Fig. 1: Overview of NOTs trapping.
Direct observation of small molecule activator binding to single PR65 protein

a Depiction of a simplified nanoaperture optical tweezer. b Cartoon representation, not to scale, of laser beam focused onto a DNH aperture with PR65 trapped between the DNH cusp. c Sample trapping signal of PR65 protein without SMAP, pink region shows when the protein has entered the trap.

Full size image

A 100 s segment of the trapped signal of PR65 and PR65 with ATUX-8385 from separate trapping events is shown in Fig. 2. A 3 Hz low pass filter applied to the raw signal of PR65 renders a comparatively homogeneous signal when compared to the signal of PR65 in the presence of ATUX-8385. In the presence of SMAP, the intensity of transmission fluctuates between two levels, showing a conformational change of the protein40. The intensity changes correlate with a change in the polarizability of the protein, hence, refractive index change of the protein16,41. The molecular weight of PR65 is 65 kDa, which is reflected in its name. In contrast, ATUX-8385 has a molecular weight of approximately 500 Da, making it about two orders of magnitude smaller than PR65. In a previous work, the change in the standard deviation of the mean was linearly related to the protein size in the NOT system37. Based on this, we would expect a .7% change in standard deviation of the signal approximately with the additional mass, but in the experiment it approximately doubles (Fig. 2b). Furthermore, our all-atom MD simulations of both the apo and ATUX-8385-bound PR65, presented in the following section, show a change in the conformation upon binding, which shifts the polarizability41.

Fig. 2: Time series signal of the trapped PR65 with (pink) and without the presence (blue) of ATUX-8385 shows the interaction between protein PR65 and SMAP.
figure 2

a The SMAP concentration used was 20 μM while the protein concentration was fixed at 10 μM. Transmission signal from different DNH is shown. Lighter signal in the background is the raw signal while the solid lines are from the 3 Hz low pass filtered signal. PDF of the filtered signals are shown to the right of the time series signal. b PDF of the filtered signal fit with Gaussian functions further indicates the presence of a bound and unbound state. For the bound state, the Gaussian curve occupies an area(AB) of 0.76 with a standard deviation(σB) of 0.004 V. Likewise, the unbound state Gaussian curve has parameters: AU = 0.22, and σU = 0.002.

Full size image

The probability density function (PDF) of the filtered signal is shown to the right of the time series signal. The 20 μM ATUX-8385 histogram was fit with two Gaussians for the bound and unbound states. The larger scattering signal voltage came from increased scattering (larger polarizability) and showed that binding of ATUX-8385 extends the protein. Extension of protein was also seen through the increase in the standard deviations of the Gaussian fits. Higher noise from the extended form of PR65 has been reported previously for point mutations33. While higher polarizability leads to a stiffer trap, it has a stronger effect on the scattering, as shown in past works37.

To determine whether the extended form corresponded to the bound or unbound state, the experiment was repeated at a concentration of 5 μM ATUX-8385 SMAP. At this reduced concentration, a lower frequency of the extended form was observed (see Supplementary Fig. SI-1), suggesting that the higher transmission level corresponds to the bound state. Further support for this conclusion is provided by molecular simulation results discussed in the following section.

Molecular simulations

We performed molecular simulations to compare experimental and simulated structural impacts. Docking simulations for ATUX-8385 revealed a novel binding site. As no structural data on the binding site/poses of the SMAP molecule ATUX-8385 on PR65 exist, we conducted a series of docking simulations on PR65 structure to determine a binding pose for ATUX-8385. Based on hydroxyl radical footprinting experiments, the K194-L198 region of PR65 was identified as a putative SMAP-binding region for a different tricyclic molecule, with K194, E197, and L198 as potential interaction sites10. Using this information, we performed guided docking simulations of ATUX-8385 onto PR65.

First, we used the Rosetta Ligand Docking Protocol 342 in the ROSIE server43 to find the optimal binding pose of ATUX-8385 near E197. The compact conformation of PR65 (PDB: 6NTS, chain A) was used for docking simulations. The selected conformation (Supplementary Fig. SI-2a, ATUX-8385 in yellow) binds to 4i, 5i, 5o and 6o, and has an interaction energy of −10.5 kcal/mol. All interaction energies are calculated using PRODIGY-LIG44.

We then conducted additional docking simulations within 40Å; radial distance from E197, so as to explore alternative docking poses and potential binding to PR65 residues near E101, where a related SMAP, DT-061, was reported to bind and stabilize the PP2A trimer45. Using AutoDock Vina 1.2.046, docking simulations were carried out for both compact (PDB: 6NTS, chain A) and extended (PDB: 1B3U, chain A) forms of PR65. The selected conformations of ATUX-8385 with PR65 in the compact form (Supplementary Fig. SI-2a, ATUX-8385 in red) and extended form (Supplementary Fig. SI-2b, ATUX-8385 in red) led to interaction energies of −6.6 kcal/mol and −9.0 kcal/mol. Interestingly, the same binding site (near 4i and 5i) was shared between the compact and extended forms, while the extended form was significantly more favorable. We note that the coordinating residues D106 and R105 (among others) are adjacent to the DT-061 binding site in PP2A.

Since the conformation found using ROSIE was energetically more favorable, we proceeded with extensive MD simulations using this pose. MD simulations initiated from the compact PR65 structure showed that ATUX-8385 binding favored an extended conformation for PR65. Two sets of MD simulations were performed: one for apo PR65 and the other for the ATUX-8385 docked form of PR65. Each set was composed of three separate runs of 704 ns in length, totaling 4.224 μs of MD trajectory across both sets. ATUX-8385 remained in its docked position throughout all three runs, further validating this new binding site as a stable interaction point for ATUX-8385 on the PR65 monomer. ATUX-8385 binding impacted the structural preferences of PR65. While the apo PR65 runs fluctuated between compact and extended forms throughout the simulations, ATUX-8385 binding stabilized the extended PR65 form rapidly within 100 ns in two runs (Fig. 3). Interestingly, in the other run, the compact form was maintained for the first 500 ns before transitioning to an extended form. Collectively, these observations demonstrate that ATUX-8385 bound to the site formed by 4i,5i,5o and 6o stabilises the PR65 and favors the extended conformation. It should be noted that the fluctuations observed in the simulations are at a much faster timescale and do not correspond to the fluctuations seen in the experiment. The Boltzmann average binding affinity for this pose (Fig. 3a), calculated over 2.112 μs of MD trajectory, was −9.6 kcal/mol.

Fig. 3: Distribution of end-to-end distances in apo and SMAP-bound PR65 from MD simulations.
figure 3

a Representative extended PR65 conformation stabilized by ATUX-8385, obtained from the third run. Distance between N30 and F578 atoms was used to define the structural form of PR65. b The end-to-end distance time evolution for runs 1–3 shown in orange, purple, and green. Dashed lines depict the end-to-end distances observed in the compact (PDB: 6NTS) and extended (PDB: 1B3U) structures of PR65. c Violin plots depicting the distribution of end-to-end distances for PR65 during the simulations under apo and ATUX-8385-bound states. For each distribution, mean values are indicated, and standard deviations are shown with thick bars.

Full size image

Observation of protein binding kinetics: NOTs single molecule measurements

The 3 Hz low-pass filter was applied to the trapped signal of a single trapping event to track the mean movement of the transmission signal, filtering out higher-frequency noise while preserving the overall trend. The choice of filtering frequency can affect the results. Testing at 10 Hz showed that filtering was insufficient, allowing spurious noise artifacts to impact the data. At 1 Hz, however, the signal was overly smoothed and distorted. We selected 3 Hz as a compromise, balancing distortion and noise artifacts. This filtered signal was then segmented into approximately 30s sections for easier processing, given the large file size. These sections were fit into step functions through a Hidden Markov Method (HMM) for identification of the two states (Fig. 4a), see Methods for details. For each trapping event, ~100 s of data was collected, segmented and analyzed. Cumulative distribution functions (CDFs) of the residence time in the bound and unbound state were fit with a single exponential function to determine association (kon) and dissociation rate (koff) constants, Fig. 4b. The good fit to a single exponential function suggests a 1:1 binding model between the protein and ligand. Individual CDF fits, association and dissociation constants can be found in Supplementary Section 3. Across the measurements made (n = 8), the signal showed a good consistency. Per 100 seconds, 111 ± 16 binding events occurred. The bound transmission intensity normalized to unbound transmission intensity was (101 ± 0.19)%. Table 1 summarizes the values for all measurements made at the 20 μM ATUX-8385 SMAP concentration.

Fig. 4: Quantification of single molecule binding kinetics.
figure 4

a Step function fit to a 30 s segment of trapped filtered signal. Higher scattering state corresponds to the protein-ligand complex. b CDFs of the residence time in bound/unbound states are fitted to a single exponential for the extraction of binding constants. A total of 94 events occurred in the 110 s analyzed. Coefficients of determination for bound and unbound CDF fitting are 0.9885 and 0.991 respectively. c Cartoon depiction of the 1:1 binding model.

Full size image
Table 1 Summary of interaction parameters between PR65 and SMAP ATUX-8385 found from single molecule NOT measurements (n = 8)
Full size table

Additional characterization of SMAP binding to PR65

We used a range of biophysical techniques to confirm in vitro binding of SMAP ATUX-8385 to PR65 and to determine the binding affinity. The relatively low solubility of the small molecule (only soluble up to 100 μM in 10% DMSO, see Supplementary Fig. SI-5) precluded the use of standard techniques such as Isothermal Titration Calorimetry (ITC), but we collected orthogonal data from three different biophysical methods to confirm that ATUX-8385 binds to PR65.

We first explored NanoDSF, a label-free, fluorescence-based technique that can be used to analyze the binding of small molecule ligands to proteins. It measures changes in the intrinsic tryptophan and tyrosine fluorescence of proteins when they undergo thermal denaturation. A shift in the melting temperature, Tm, of the protein can be expected upon small molecule binding, as small molecules can locally or globally stabilize or destabilize the protein structure. Thermal denaturation of 2 μM PR65 was measured in 10% DMSO versus in the presence of 100 μM ATUX-8385, 10% DMSO. A positive shift in the melting temperature of PR65 was observed in the presence of ATUX-8385 (Tm PR65 = 52.7 ± 0.1 °C, Tm PR65+ATUX = 53.5 ± 0.1 °C), indicating binding and slight stabilization of the protein structure (Fig. 5a, b). The relatively low observed shift in Tm suggested a binding affinity in the micromolar range.

Fig. 5: Biophysical characterization of small molecule ATUX-8385 binding with protein PR65.
figure 5

a NanoDSF traces of the thermal denaturation of PR65 in the absence (green traces) and in the presence of ATUX-8385 (orange traces). The data shown are the first derivative of the ratio of fluorescence intensities read at 350 nm and 330 nm (dFIR (350 nm/330 nm)). The global minimum corresponds to the melting temperature of the protein, Tm. A shift towards higher Tm values indicates an increase in stability induced upon SMAP binding. b Extracted Tm values from the NanoDSF traces indicate an upwards shift in the presence of ATUX-8385(Tm PR65 = 52.7 ± 0.1 °C, Tm PR65+ATUX-8385 = 53.5 ± 0.1 °C, N = 7, p < 0.001 via an unpaired t-test). The NanoDSF experiments were performed with a Prometheus NanoDSF instrument (NanoTemper Technologies). 2 μM of PR65 in PBS, 2 mM DTT, was incubated either with 10% DMSO, or 100 μM ATUX-8385 in a final concentration of 10% DMSO, and thermal denaturation was performed from 20 °C to 90 °C with a 1 °C/min rate. c 19F NMR recorded at 298 K on ATUX-8385 indicates binding of the small molecule to PR65. The NMR spectra show a drop in signal intensity upon binding of the small molecule, while the change in chemical shift is small as indicative of relatively weak binding. 19F 1D CPMG experiments were performed for 100 μM ATUX-8385 in 10% DMSO (blue: transverse period 2 ms; cyan: transverse period 102 ms) or for 100 μM ATUX-8385 in 10% DMSO with 5 μM PR65 (red: transverse period 2 ms; magenta: transverse period 102 ms) on a Bruker Avance III 600 MHz (19F 564 MHz). d Fluorescence polarization experiments show ATUX-8385 binding on PR65 with a KD in the low micromolar range. Shown are the fluorescence polarization values (mP) for 2.5 μM, 5 μM, and 10 μM ATUX-8385 upon PR65 titration, N = 3. One-site fitting of the data (see methods for equation) gives KD = 9.4 ± 1.4 μM.

Full size image

To further establish the binding of ATUX-8385 to PR65, we used 19F-NMR, which probes Fluorine-19 on the ligand (ATUX-8385). The interaction between ATUX-8385 and PR65 was detected by measuring the signal intensity and the 19F relaxation rate of the small molecule in the presence and absence of PR65. The 19F NMR spectra were collected for 100 μM ATUX-8385 in 10% DMSO, and for 100 μM ATUX-8385 in 10% DMSO with 5 μM PR65. A singlet signal was detected at -57.46 ppm for ATUX-8385, corresponding to the -CF3 group of the molecule. Only a very small change in signal position was observed upon addition of PR65. An increase in the relaxation rate R2 upon addition of PR65 (R2ATUX = 6.93 ± 0.05, R2ATUX+PR65 = 1.85 ± 0.05), revealed an increase in the rotational correlation time, indicative of binding of ATUX-8385 to PR65 (Fig. 5c).

To quantify the binding affinity of ATUX-8385 for PR65, we exploited the fact that ATUX-8385 is a fluorescent molecule. ATUX-8385, when excited at 320 ± 10 nm wavelength light, emits at 340–550 nm, with a maximum fluorescence intensity at 370 nm (Fig. SI-5). Fluorescence polarization (FP) is, therefore, a useful technique that can be used in this system to calculate the binding affinity. FP is based on the principle that the rotational motion of a fluorescent molecule, in this case ATUX-8385, affects the polarization of emitted light. When ATUX-8385 was bound to PR65, its rotational motion was reduced, resulting in higher polarization values. FP was performed at 40 °C to match the optical trapping experiments. PR65 was titrated in concentrations ranging from 1 μM to 80 μM to samples containing 2.5 μM, 5 μM, and 10 μM of ATUX-8385 (10% DMSO) (Fig. 5d). By fitting the data to a one-site binding curve, the equilibrium dissociation constant, KD, was determined to be 9.4 ± 1.4 μM.

Discussion

As PP2A holds a central role in many biological processes, there have been clinical investigations into PP2A directed therapeutics for preventing dysregulation in the enzyme. Here we explored a SMAP (ATUX-8385) of tricyclic sulfonamide class that binds directly to the scaffold unit of PP2A, PR65, which has been shown to reactivate PP2A in a previous study9. We employed several approaches to gain novel insights into the effects of ATUX-8385 binding on the structure and dynamics of PR65, and the binding kinetics of PR65 and ATUX-8385.

Through docking simulations, a novel binding site, comprising residues at 4i, 5i, 5o, and 6o, for ATUX-8385 was identified. Subsequent extensive MD simulations confirmed the stability of this site and demonstrated that ATUX-8385 impacts the structural dynamics of PR65. In our simulations, ATUX-8385 binding was observed to favor the extended form. Specifically, binding of ATUX-8385 stabilized the extended form in two out of three simulations and the third run converged to an extended form towards the end of the simulation.

The docking of ATUX-8385 to the 4i and 5i shows promise, as indicated by its vicinity to known SMAP binding site. These novel binding sites differ from previously identified sites and suggest a potentially unique mechanism of action. These findings warrant further analysis in future work, including extensive MD simulations and experimental validation, to better understand the binding dynamics and functional implications of ATUX-8385 on PR65.

We produced two separate measurements of dissociation constants in this work: single molecule and ensemble. The dissociation constant found using the NOT approach comes from a single measurement on a single protein at a single SMAP concentration, while the ensemble measurement utilized the fluorescent property of ATUX-8385 to perform fluorescence polarization experiments. The ensemble measurement required repeated measurements at different PR65 concentrations to establish the dissociation constant. The KD obtained from both measurements are in good agreement with each other, and characterized ATUX-8385 as a molecule of weak affinity to bind to PR65. Using NOT, the observation of complex formation and dissociation in real time allows for determination of kon and koff. The change in the scattering signal with complex formation is consistent with the elongation of the protein upon binding, which provides biophysical data about the interaction. This stabilization in the elongated form was also confirmed in MD simulations. Our past work quantified this elongation for various point mutations of PR6533.

Since the NOT technique depends on detecting light scattering changes induced by the activator upon binding, its applicability may be limited to binding pairs that produce a measurable scattering change, as not all interactions may result in detectable shifts. However, due to the same principle of the technique, the single-molecule binding assay is not restricted to the typical picomolar- or nanomolar range for techniques such as total internal reflection microscopy, smFRET. Fluorescence-based single molecule techniques are often confined to a low concentration barrier, as high concentration of fluorescent molecules results in an increase in background signal that in turn obscures binding events34,47. This is especially limiting as many biological processes are transient and have a dissociation constant in the micromolar range. Therefore, the NOT technique has an appropriate niche to study kinetics that can be observed on a timescale ranging from a microsecond to hundreds of seconds. For longer observation, a stable microscope setup is required to avoid drift over extended duration.

The NOT has an interesting feature of being able to tune the temperature by adjusting the trapping laser power40,48,49,50. In the present work, we limited the power to achieve temperatures near physiological values; however, future studies can be undertaken to explore the impact of temperature on the binding kinetics. We did attempt a few temperature values by varying the laser power and saw insignificant variation in the dissociation constant over 4 °C.

In this study, we explored a SMAP (ATUX-8385) of tricyclic sulfonamide class that binds to PR65 and has been shown to reactivate PP2A in a previous study aimed at cancer therapeutics9. We used NOT, a single molecule technique, to determine the dissociation constant of 13.6 ± 2.5 μM, agreeing with fluorescence anisotropy measurements while not requiring a titration: the experiment was performed on a single protein molecule and was highly repeatable. The NOT approach showed increased light scattering upon binding, which is consistent with an extension of the protein upon binding, as found in MD simulations. The MD simulations also elucidated the binding site. Therefore, the NOT, combined with all-atom MD simulations, provided information on the conformational changes taking place upon SMAP binding. Such mechanistic data are useful for clinical certification of SMAPs.

Label-free single particle analysis techniques are changing the landscape of biomolecule investigations, with recent successful commercialization of iScat51,52, nanopores53,54 and tethering-based force-measurement55,56 techniques. The NOT approach allows for valuable data extraction on small molecule binding for drug discovery applications, not only quantifying the binding affinity, but also giving insight into the impact of binding on the structure of the protein. It also has the advantage of working on a single molecule, not requiring additional steps like titrations or labeling, using less sample and allowing for tuning the local temperature at the protein by changing the laser power. Therefore, we believe that there is potential for this technique to be widely adopted in drug discovery and other research areas.

Methods

Protein expression and purification

GST-tagged PR65 was expressed in E. coli as described previously57. In brief, plasmid encoding PR65 was transformed into chemically competent C41 E. coli (Komander laboratory, MRC-LMB, Cambridge). Cultures were grown at 37 °C in 2xYT medium containing ampicillin (50 μg/ml) until an OD600 of 0.6–0.8 was reached. Protein expression was induced with 500 μM isopropyl-β-d-thiogalactopyranoside (IPTG) (Generon) at 25 °C overnight. Cells were harvested by centrifugation at 4000 × g for 10 min at 4 °C before resuspending in lysis buffer [50 mM tris-HCl (pH 7.5), 500 mM NaCl, 2 mM dithiothreitol (DTT)] supplemented with EDTA-free protease inhibitor cocktail (Sigma-Aldrich) and deoxyribonuclease (DNase) I (Sigma-Aldrich). The cells were lysed by passing the suspension two to three times through an Emulsiflex-C5 (AVESTIN) at pressures of 10,000–15,000 psi. Soluble protein was separated from cell debris and other insoluble fractions by centrifugation at 35,000 × g for 35 min at 4 °C. The soluble protein fraction was applied to glutathione resin [Amintra Affinity, EGTA (0.5 g/liter), 2 mM DTT], the GST tag was cleaved with thrombin, and PR65 was eluted using a gravity column. After washing the column, the protein was subsequently eluted using a 20 × column-volume salt gradient from 0 to 1 M NaCl. MonoQ fractions containing the protein were concentrated before application to a HiLoad 26/600 Superdex 200 pg (GE Healthcare) equilibrated in phosphate-buffered saline (pH 7.4) and 2 mM DTT. Samples were analyzed by SDS-PAGE (polyacrylamide gel electrophoresis) comparing the lysed, flowthrough, and eluted fractions. The identity of the protein was confirmed via MALDI mass spectrometry (Department of Chemistry, University of Cambridge).

DNH nanoaperture fabrication

The DNH fabrication follows our past approaches15,58. A microscope slide with dimensions 75 × 50 × 1 mm was divided into three roughly equal pieces using a diamond scriber. The slides were sonicated in an ethanol bath for 10 min, then rinsed with acetone and isopropyl alcohol before being dried with nitrogen. To prepare the polystyrene spheres (PS) solution, 7.5 μL of ThermoFisher 300 nm PS was mixed into 1 mL of ethanol. Next, 10 μL of the PS-ethanol solution was deposited onto each slide in a zigzag pattern and left to dry overnight. The slides were plasma etched using a Harrick PDC plasma cleaner for 135 s to fuse randomly distributed dimers and reduce the cusp size as a function of etching time. A 7 nm layer of titanium was sputtered onto the slides to act as an adhesive layer for gold. Subsequently, a 70 nm layer of gold was sputtered onto the slides. To remove the PS, the slides were sonicated for 7 min. For verification purposes, the slides were then imaged using a SEM Hitachi S-4800 (Supplementary Fig. SI-6).

Sample preparation for trapping experiments

The protein solutions comprised of phosphate-buffered saline solution with 0.01% (w/v) concentration, 2 mM of Dithiothreitol and 5% (v/v) Dimethylsulfoxide (DMSO). DMSO was used as a solvent for ATUX-8385. PR65 concentration was 10 μM while the ATUX-8385 concentration in solution was 20 μM. It is important to note that upon laser illumination, ohmic losses in the metal lead to local heating at the nanoaperture. Studies have shown that the temperature increase ranges between 0.64 K/mW to 4.25 K/mW48,49. The temperature increase is experimentally shown to be dependent on the intensity of the incident laser, yet independent of aperture size, DNH cusp, and laser polarization49. The local temperature increase at the aperture was 0.64 K/mW39,48. Unless otherwise stated, trapping experiments were conducted at an incident power that raised the temperature of the trapped protein to approximately 35 °C. A 50 × 2 × 0.13 mm coverslip was rinsed with acetone and ethanol, then dried with nitrogen gas. A microwell was then formed by attaching an image spacer (Grace Bio-Laboratories GBL-654008-100EA) onto the coverslip. 9.5 μL of protein solution was then micro-pipetted into the microwell. The DNH slide, with the gold side facing the solution, was placed and sealed onto the image spacer.

Trapping setup: inverted microscope

A collimated laser beam exiting a fiber passes through a linear polarizer and then a half-wave plate, which adjusts the polarization. The beam fully fills the back aperture of a 100 ×, 1.25 NA Nikon E-plan oil immersion objective. From there, the laser passes through the objective, sample, and condenser (10 × objective, Nikon MRP70100). After the condenser, the beam is focused onto the APD. During trapping experiments, the laser is slightly under-focused to maximize transmission to the APD.

NOT data acquisition and analysis

Data was acquired at a sampling rate of 100 KHz using an Advantech USB-4711A data acquisition unit and analyzed in Python and vbFRET. vbFRET was used for HMM fitting of the data for identifying state transitions. As the program was designed for single-molecule FRET data, input to the program must contain two columns of data corresponding to the donor and acceptor intensity59. The FRET efficiency is calculated as ({I}_{{rm{Eff}}}=frac{{I}_{A}}{{I}_{D}+{I}_{A}}). Given that we want to fit the trapping transmission data (IT) as IEff, we must have the equivalence of IEff = IT. This is achieved by setting IT = IA and ID = 1 − IA. HMM fitting was applied with the setting of the number of maximum states as 2 and 20 fitting attempts.

After completing the HMM fitting, the fit was verified visually and refined using a custom Python script that corrects rapid state transitions in the time series data. The script identifies binary state values from the HMM fit and applies a 0.05 s duration threshold to detect potential spurious transitions likely caused by noise rather than actual state changes. For each transition detected within or below the threshold, the script compares the mean signal value during the brief state to the midpoint between the high and low states. If the mean value deviates by less than 10% from this midpoint, the transition is considered noise, and the state is reverted to the previous one. This correction process helps to remove spurious transitions, preserving both event duration and the integrity of the fit in the corrected time series. Refer to Supplementary Fig. SI-7 for the results of pre- and post-fitting improvements. Residence times are then collected after a quality check of the fitting.

Once all bound and unbound residence times are collected, we plot residence times as a CDF. This is then fit to the equation of (1-exp left(-frac{x}{tau }right)) using scipy.optimize.curve_fit, where τ is the time constant. Using equations given by the law of mass action, we calculated the binding rate constants as well as dissociation constant. We do not make the assumption of [L]total ≈ [L]free as this would introduce error to kon and KD as ligand depletion would remain unaccounted. Using time constants τunbound and τbound from the CDF fit, we use quadratic Eq. (5) to calculate for the free concentration of SMAP as well as τon and KD.

The protein concentration [P] and ligand concentration [L] bind as:

$$[P]+[L]mathop{rightleftharpoons }limits_{{k}_{{rm{off}}}}^{{{k}_{{rm{on}}}}}[PL]$$
(1)

The on binding rate is given by:

$${k}_{{rm{on}}}=frac{1}{{tau }_{{rm{on}}}cdot {[L]}_{{rm{free}}}}$$
(2)

and the off binding rate is given by:

$${k}_{{rm{off}}}=frac{1}{{tau }_{{rm{off}}}}$$
(3)

From these, the dissociation constant can be determined as:

$${K}_{D}=frac{{tau }_{{rm{on}}}cdot {[L]}_{free}}{{tau }_{{rm{off}}}}=frac{[P][L]}{[PL]}={tau }^{* }({[L]}_{t}-[PL])=frac{({[P]}_{t}-[PL])({[L]}_{t}-[PL])}{[PL]}$$
(4)

where we have exploited the relation:

$$left({[PL]}^{2}+{tau }^{* }{[PL]}^{2}right)-[PL]cdot left({[P]}_{t}+{[L]}_{t}+{tau }^{* }{[L]}_{t}right)+{[P]}_{t}cdot {[L]}_{t}=0$$
(5)

The Gibbs free energy of association is calculated using this equation.

$$Delta {G}^{circ },=RTln {K}_{{rm{D}}}$$
(6)

Docking simulations

ATUX-8385 was docked to PR65 using ROSIE43 and AutoDock Vina60. Before initiating the docking simulations using AutoDock Vina, hydrogens were added to both protein and ligand, and the structures were saved in the PDBQT file format utilizing AutoDock Tools (ADT) 1.5.646. Subsequently, partial charges for all ligands were calculated using the Gasteiger method within ADT. The added hydrogens in ligand were manually inspected. Following this verification, a 3-dimensional grid box was constructed with varying sizes ranging from 20 Å to 40 Å, maintaining a grid spacing of 1 Å. The center of the initial grid box was set based on the ATUX-8385 docked structure obtained from ROSIE.

System preparation for MD simulations

The structure of monomeric PR65 was derived from the DT-061 bound trimeric PP2A structure (PDB: 6NTS45) with a resolution of 3.63 Å, in which PR65 adopts a compact form. The coordinates for ATUX-8385 docked at 4i, 5i, 5o, and 6o were obtained from our docking simulations. Each of apo and ATUX-8385 bound PR65 structures was solvated in a water box containing explicit TIP3P61 water molecules, with a box size 144 Å in all directions. To neutralize the system and set the ion concentration to 150 mM NaCl, Na+ and Cl- ions were added. The system size was approximately 286,900 atoms. All system preparation steps were conducted using VMD62.

Details of MD simulations

MD simulations were conducted using the NAMD 363 software, employing the CHARMM36 all-atom additive protein force field64. A 2 fs time step was used throughout the simulations. Temperature was maintained at 310 K via Langevin dynamics, utilizing a damping coefficient of 1 ps-1, and pressure was held at 1 atm using the Langevin Nosé–Hoover method with an oscillation period of 100 fs and a damping time scale of 50 fs. Van der Waals interactions were calculated with a cut-off distance of 12 Å, while the particle-mesh Ewald method was applied for long-range electrostatic interactions. Two rounds of system minimization and equilibration were executed before each production run. Initially, the protein structure was kept fixed and the system was to 10,000 minimization steps, followed by a 1 ns of equilibration. This first round of minimization-equilibration was designed to equilibrate the solvent around the protein. Subsequently, we performed a second round of minimization-equilibration, in which the system underwent an additional 10,000-step minimization without any restrictions on protein structure and dynamics, followed by a 2 ns of equilibration, applying harmonic constraints (k = 1 kcal/mol/Å2) only on the Cα atoms. After these preparatory simulations, we removed all constraints and initiated the production runs. Apo PR65 simulations were extended to a total duration of 704 ns from our previous study33, where each simulation was 654 ns in length. In contrast, all system preparation and production runs for ATUX-8385 bound PR65 were entirely performed in this study.

Nano-differential scanning calorimetry (nanoDSF)

PR65 samples (2 μM) in PBS, 2 mM DTT were mixed with ATUX-8385 to a final concentration of 100 μM, 10% DMSO and in a total volume of 20 μL. The mixtures were allowed to equilibrate at room temperature for 10 min. High-sensitivity capillaries (NanoTemper) were filled with the equilibrated PR65-SMAP mixtures using capillary forces, ensuring no air bubbles were present. The capillaries were loaded into the Prometheus NT.48, and a thermal ramp from 20 °C to 90 °C at a rate of 1 °C per minute was applied. Intrinsic tryptophan and tyrosine fluorescence at 330 nm and 350 nm was continuously monitored. The melting temperature (Tm) was determined from the first derivative of the fluorescence ratio (I350/I330) with respect to temperature. The Tm corresponds to the inflection point of the melting curve.

Fluorescence polarization (FP)

A fixed concentration of ATUX-8385 (2.5 μM, 5 μM, and 10 μM) was mixed with varying concentrations of the protein (0–80 μM) in a total volume of 20 μL in each well of a 384-well microplate (Greiner), to a final DMSO concentration of 10%. The plate was incubated at 40 °C for 30 min to ensure binding equilibrium. Fluorescence polarization was measured at 40 °C using an excitation wavelength of 295 ± 10 nm and an emission wavelength of 360 ± 20 nm, on a CLARIOStar (BMG Labtech) measuring the light in parallel and perpendicular planes relative to the excitation plane. Wells containing only ATUX-8385 (without PR65) and wells containing only PR65 were used as controls to determine background fluorescence and polarization values. Fluorescence polarization (P) was calculated using the equation:

$$P=frac{{I}_{parallel }-{I}_{perp }}{{I}_{parallel }+{I}_{perp }}$$
(7)

where I and I are the intensities of fluorescence parallel and perpendicular to the excitation plane, respectively. Non-linear regression analysis was performed using a one-site binding model to fit the data and determine the dissociation constant KD, using the following equation:

$$P=frac{{P}_{max }cdot {C}_{{rm{PR65}}}}{{k}_{D}+{C}_{{rm{PR65}}}}$$
(8)

where P is the polarization, ({P}_{max }) is the maximum polarization at the plateau, CPR65 is the concentration of PR65 (μM) and KD is the equilibrium dissociation constant.

19F Nuclear magnetic resonance (NMR)

NMR samples containing 5% D2O were prepared with ATUX-8385 added to a final concentration of 100 μM in 10% DMSO. NMR spectra were recorded at 298 K on a 600 MHz (1H) Bruker Avance III spectrometer equipped with a 5 mm QCI HFCN/z cryoprobe (19F 564 MHz) (Department of Biochemistry, University of Cambridge). 1D 19F NMR data were obtained with 256 scans for ATUX-8385 in the presence of PR65 and 128 scans without the protein and recorded as 1D CPMG experiments (the separation between neighboring 180 pulses was 100 μs). To assess binding of ATUX-8385 the intensities of pairwise experiments recorded with 2 ms and 102 ms transverse delay time were compared and the intensity ratios were converted into transverse relaxation rate R2 values.

Related Articles

Iron homeostasis and ferroptosis in muscle diseases and disorders: mechanisms and therapeutic prospects

The muscular system plays a critical role in the human body by governing skeletal movement, cardiovascular function, and the activities of digestive organs. Additionally, muscle tissues serve an endocrine function by secreting myogenic cytokines, thereby regulating metabolism throughout the entire body. Maintaining muscle function requires iron homeostasis. Recent studies suggest that disruptions in iron metabolism and ferroptosis, a form of iron-dependent cell death, are essential contributors to the progression of a wide range of muscle diseases and disorders, including sarcopenia, cardiomyopathy, and amyotrophic lateral sclerosis. Thus, a comprehensive overview of the mechanisms regulating iron metabolism and ferroptosis in these conditions is crucial for identifying potential therapeutic targets and developing new strategies for disease treatment and/or prevention. This review aims to summarize recent advances in understanding the molecular mechanisms underlying ferroptosis in the context of muscle injury, as well as associated muscle diseases and disorders. Moreover, we discuss potential targets within the ferroptosis pathway and possible strategies for managing muscle disorders. Finally, we shed new light on current limitations and future prospects for therapeutic interventions targeting ferroptosis.

Type 2 immunity in allergic diseases

Significant advancements have been made in understanding the cellular and molecular mechanisms of type 2 immunity in allergic diseases such as asthma, allergic rhinitis, chronic rhinosinusitis, eosinophilic esophagitis (EoE), food and drug allergies, and atopic dermatitis (AD). Type 2 immunity has evolved to protect against parasitic diseases and toxins, plays a role in the expulsion of parasites and larvae from inner tissues to the lumen and outside the body, maintains microbe-rich skin and mucosal epithelial barriers and counterbalances the type 1 immune response and its destructive effects. During the development of a type 2 immune response, an innate immune response initiates starting from epithelial cells and innate lymphoid cells (ILCs), including dendritic cells and macrophages, and translates to adaptive T and B-cell immunity, particularly IgE antibody production. Eosinophils, mast cells and basophils have effects on effector functions. Cytokines from ILC2s and CD4+ helper type 2 (Th2) cells, CD8 + T cells, and NK-T cells, along with myeloid cells, including IL-4, IL-5, IL-9, and IL-13, initiate and sustain allergic inflammation via T cell cells, eosinophils, and ILC2s; promote IgE class switching; and open the epithelial barrier. Epithelial cell activation, alarmin release and barrier dysfunction are key in the development of not only allergic diseases but also many other systemic diseases. Recent biologics targeting the pathways and effector functions of IL4/IL13, IL-5, and IgE have shown promising results for almost all ages, although some patients with severe allergic diseases do not respond to these therapies, highlighting the unmet need for a more detailed and personalized approach.

Enhancer reprogramming: critical roles in cancer and promising therapeutic strategies

Transcriptional dysregulation is a hallmark of cancer initiation and progression, driven by genetic and epigenetic alterations. Enhancer reprogramming has emerged as a pivotal driver of carcinogenesis, with cancer cells often relying on aberrant transcriptional programs. The advent of high-throughput sequencing technologies has provided critical insights into enhancer reprogramming events and their role in malignancy. While targeting enhancers presents a promising therapeutic strategy, significant challenges remain. These include the off-target effects of enhancer-targeting technologies, the complexity and redundancy of enhancer networks, and the dynamic nature of enhancer reprogramming, which may contribute to therapeutic resistance. This review comprehensively encapsulates the structural attributes of enhancers, delineates the mechanisms underlying their dysregulation in malignant transformation, and evaluates the therapeutic opportunities and limitations associated with targeting enhancers in cancer.

Local structural flexibility drives oligomorphism in computationally designed protein assemblies

Many naturally occurring protein assemblies have dynamic structures that allow them to perform specialized functions. Although computational methods for designing novel self-assembling proteins have advanced substantially over the past decade, they primarily focus on designing static structures. Here we characterize three distinct computationally designed protein assemblies that exhibit unanticipated structural diversity arising from flexibility in their subunits. Cryo-EM single-particle reconstructions and native mass spectrometry reveal two distinct architectures for two assemblies, while six cryo-EM reconstructions for the third likely represent a subset of its solution-phase structures. Structural modeling and molecular dynamics simulations indicate that constrained flexibility within the subunits of each assembly promotes a defined range of architectures rather than nonspecific aggregation. Redesigning the flexible region in one building block rescues the intended monomorphic assembly. These findings highlight structural flexibility as a powerful design principle, enabling exploration of new structural and functional spaces in protein assembly design.

Targeting of TAMs: can we be more clever than cancer cells?

With increasing incidence and geography, cancer is one of the leading causes of death, reduced quality of life and disability worldwide. Principal progress in the development of new anticancer therapies, in improving the efficiency of immunotherapeutic tools, and in the personification of conventional therapies needs to consider cancer-specific and patient-specific programming of innate immunity. Intratumoral TAMs and their precursors, resident macrophages and monocytes, are principal regulators of tumor progression and therapy resistance. Our review summarizes the accumulated evidence for the subpopulations of TAMs and their increasing number of biomarkers, indicating their predictive value for the clinical parameters of carcinogenesis and therapy resistance, with a focus on solid cancers of non-infectious etiology. We present the state-of-the-art knowledge about the tumor-supporting functions of TAMs at all stages of tumor progression and highlight biomarkers, recently identified by single-cell and spatial analytical methods, that discriminate between tumor-promoting and tumor-inhibiting TAMs, where both subtypes express a combination of prototype M1 and M2 genes. Our review focuses on novel mechanisms involved in the crosstalk among epigenetic, signaling, transcriptional and metabolic pathways in TAMs. Particular attention has been given to the recently identified link between cancer cell metabolism and the epigenetic programming of TAMs by histone lactylation, which can be responsible for the unlimited protumoral programming of TAMs. Finally, we explain how TAMs interfere with currently used anticancer therapeutics and summarize the most advanced data from clinical trials, which we divide into four categories: inhibition of TAM survival and differentiation, inhibition of monocyte/TAM recruitment into tumors, functional reprogramming of TAMs, and genetic enhancement of macrophages.

Responses

Your email address will not be published. Required fields are marked *