Kevin M.
Burridge‡
,
Ben A.
Shurina‡
,
Caleb T.
Kozuszek
,
Ryan F.
Parnell
,
Jonathan S.
Montgomery
,
Jamie L.
VanPelt
,
Nicholas M.
Daman
,
Robert M.
McCarrick
,
Theresa A.
Ramelot
,
Dominik
Konkolewicz
* and
Richard C.
Page
*
Department of Chemistry and Biochemistry, Miami University, 651 E High St., Oxford, OH 45056, USA. E-mail: d.konkolewicz@MiamiOH.edu; pagerc@MiamiOH.edu
First published on 3rd June 2020
Rational design of protein–polymer bioconjugates is hindered by limited experimental data and mechanistic understanding on interactions between the two. In this communication, nuclear magnetic resonance (NMR) paramagnetic relaxation enhancement (PRE) reports on distances between paramagnetic spin labels and NMR active nuclei, informing on the conformation of conjugated polymers. 1H/15N-heteronuclear single quantum coherence (HSQC) NMR spectra were collected for ubiquitin (Ub) modified with block copolymers incorporating spin labels at different positions along their backbone. The resultant PRE data show that the conjugated polymers have conformations biased towards the nonpolar β-sheet face of Ub, rather than behaving as if in solution. The bioconjugates are stabilized against denaturation by guanidine-hydrochloride, as measured by circular dichroism (CD), and this stabilization is attributed to the interaction between the protein and conjugated polymer.
Within the field of protein polymer conjugates there has been much seemingly conflicting data produced, without a satisfying unifying model. Intuitively, soluble proteins with a predominantly hydrophilic surface should prefer interactions with polar polymers featuring complementary functional groups. However, a recent publication showed that proteins such as glucose-oxidase can be modified with polystyrene to the point of becoming globular nanoparticles, with improved or retained activity.11 Recent work on surface-immobilization of enzymes has shown that lipase catalytic performance can be dramatically improved by attachment to hydrophilic surface-grafted polymers.12 In a subsequent publication, by the same group, it was found that there is a trade-off in this system – more protein–polymer attachments resulted in higher stability but lower activity due to decreased motion in the folded state of the protein.13 This experimental result conflicts with recent work on soluble chymotrypsin conjugates, where MD models showed more motion in the conjugate with the lowest enzymatic activity, a manifestation of a measured decrease in the substrate binding constant, Km.14
The degree, specificity, location, and timescale of protein–polymer associations are all open questions. It has been found that even small proteins such as insulin can have improved or decreased activity depending on the location of the modification. In all cases, insulin was stabilized against heat degradation.15 If the polymer is strongly interacting with the insulin surface, it would be expected to inhibit receptor binding no matter where it is attached, while the stabilization observed in these conjugates would intuitively be caused by some sort of physical association.
The understanding of interactions between proteins and polymers, and the proximities and conformations of polymers within bioconjugates have relied on molecular dynamics (MD) simulations, small-angle neutron scattering, and a single crystal structure. Although MD gives hypotheses on interactions and polymer conformations,14,16–20 the results are limited by the accuracy of the force fields used and an absence of data from physical samples. Neutron scattering data, while providing useful insights on polymer conformation and protein dynamics,21–23 are limited by: the quality of models and fits; specialized instrumentation; and resultant coarse-grained information. While the single existing crystal structure was determined from data collected on a protein–polymer bioconjugate, the polymer was not resolved, likely due to polymer flexibility.24 Other isolated coarse-grained experiments probing polymer behavior have been published, such as a 15A resolution Cryo-EM structure25 of a PEG-bottlebrush-modified virus-like particle (a privileged system with high symmetry), and motional dynamics probed by electron paramagnetic resonance (EPR) on the surface of polymer-modified self-assembling peptide amphiphiles.26
Finally, most structural information has been on PEG-conjugates, which constitute a subset of protein–polymer bioconjugates. Thus, experiments that report on conformations and proximities within protein–polymer bioconjugates are greatly needed.
There is no consensus on the impact of polymer–protein non-covalent associations and conformation on stability.16,17,27,28 One model proposes that protein–polymer associations are stabilizing through non-covalent bonds,17 others suggest protein–polymer associations are destabilizing by disrupting the hydration shell,28 or that the polymer acts as a chaperone.16 The dearth of experimental data and conflicting models limit the mechanistic understanding of protein–polymer hybrids and rational bioconjugate engineering. Often distinct metrics of a bioconjugate's ability to tolerate environmental challenges are used, such as stability measured by loss of activity or function in the bioconjugate upon exposure to a challenge, versus structural stability where the bioconjugate's fold or structure is lost upon exposure to a challenge. The various metrics of stability and lack of experimental data on conformations, interactions and proximities of the synthetic and biological components limit rational design of bioconjugates for enhanced performance. To address the limited information on the interplay of biological and synthetic components of a bioconjugate, this work applies paramagnetic relaxation enhancement (PRE) to bioconjugates, giving information on the conformation of these bioconjugates. PRE is a nuclear magnetic resonance (NMR) technique, providing up to 30 Å distance information.29,30 In PRE, paramagnetic species increase the relaxation rate of NMR-active nuclei,31 broadening resonances that either diminish (15–30 Å) or eliminate (0–15 Å) when the paramagnetic center approaches the NMR active nucleus (Scheme 1B). Through the PRE method, the proximity of a spin-labeled polymer to NMR active nuclei on a protein can be evaluated, creating a platform to correlate protein–polymer interactions and proximities with bioconjugate performance. This can open the door to rational bioconjugate design and engineering.
Molecular weight and chain length data for each block are presented in Tables 1 and S1,† with number average molecular weight by NMR determined using eqn S1.† Four units of spin label precursor N-(2,2,6,6-tetramethyl-4-piperidinyl)-2-propenamide (TMPA) were targeted to account for (1) Poisson-like distribution of monomers incorporated during polymerization33 and (2) the non-quantitative yield of TEMPO in post-polymerization oxidation.
The dispersity (Mw/Mn) and length of the first blocks are critical as these parameters determine the distance between the spin label and Ub when conjugated. In general, the number average molecular weights (Mn) of each block agreed well with the theoretical values, where the Mn was determined using nuclear magnetic resonance (NMR) to calculate the number of DMAm and TMPA repeating units per polymer on average. Molar mass dispersities were determined by size exclusion chromatography (SEC) and we obtained Mw/Mn < 1.12 for the first block of all three block copolymers (Table 1), and Mw/Mn < 1.35 for all polymers overall.
Block polymer | 0-Block | A-Block | B-Block | C-Block |
---|---|---|---|---|
1st block (DMAm) | 34 | 5 | 15 | 25 |
2nd block (TMPA) | 0 | 4 | 4 | 4 |
3rd block (DMAm) | 0 | 25 | 15 | 5 |
M n,theory-block1 | 3600 | 710 | 1700 | 2700 |
M n,NMR-block1 | 3900 | 780 | 1900 | 3100 |
M w/Mn,SEC-block1 | 1.10 | 1.12 | 1.12 | 1.12 |
rTEMPO (Å) | — | 13 ± 3 | 21 ± 4 | 29 ± 6 |
M n,theory-overall | 3600 | 4000 | 4000 | 4000 |
M n,NMR-overall | 3900 | 4700 | 4400 | 4900 |
M w/Mn,SEC-overall | 1.10 | 1.32 | 1.26 | 1.25 |
The chain lengths (Table 1) place the TEMPO spin label distances (rTEMPO) in key PRE-effect ranges of approximately rTEMPO = 13 ± 3 Å for A-block, rTEMPO = 21 ± 4 Å for B-block, and rTEMPO = 29 ± 6 Å for C-block, based on end-to-end distances (eqn S2–S8†). These values were estimated based on the end-to-end distance estimated from the weight average molecular weight (Mw) and the dispersity (Mw/Mn) of the poly(DMAm) block between the attachment site and the spin labeled block. This is because the spin label will be placed precisely one repeat unit after this first poly(DMAm) block. The dispersity of the first block would impact the placement of the spin labels most significantly, since the first spin label is found exactly 1 unit after the first block synthesized, therefore uncertainty in the placement of the spin label is relatively minimal, since the first poly(DMAm) blocks have dispersities in the range of 1.10–1.12.34
To account for the impact of molar mass dispersity on the placement of the spin label, the standard deviation (σM) in molecular weights of the first poly(DMAm) block was estimated for each of the A, B and C-block polymer.35 The uncertainty in the distance between the protein's surface and the spin label was found by estimating the end-to-end distance of the first poly(DMAm) block based: Mw, giving the mean value of rTEMPO; Mw + σM, giving the mean squared end-to-end distance of a polymer of molecular weight Mw + σM, which is denoted as 〈r(Mw + σM)2〉1/2; and finally Mw − σM, giving the mean squared end-to-end distance of a polymer of molecular weight Mw − σM: 〈r(Mw − σM)2〉1/2. The absolute difference between rTEMPO and 〈r(Mw + σM)2〉1/2 and the absolute distance between rTEMPO and 〈r(Mw − σM)2〉1/2 is calculated. The larger of |rTEMPO − 〈r(Mw + σM)2〉1/2| and |rTEMPO − 〈r(Mw − σM)2〉1/2| was used as the uncertainty in the placement of the spin label.
Oxidation by mCPBA gave TEMPO spin labeled polymers (Fig. S1†). A control DMAm homopolymer (0-block) of the same length as the block copolymers was also synthesized. This polymer was used to gauge the effect on the Ub spectrum from increased viscosity due to free polymer and structural perturbations from conjugation (Fig. 1b). Each polymer was conjugated to Ub using N-(3-dimethylaminopropyl)-N′-ethylcarbodiimide hydrochloride/N-hydroxysuccinimide (EDC/NHS) targeting amine groups within Ub, such as the N-terminus and lysine side-chain ε-amines as indicated in Scheme 1.36 Conjugation was confirmed using polyacrylamide gel electrophoresis (PAGE) for the 0-block conjugate (Fig. 1a) and all spin labeled conjugates (Fig. S2†).
To determine the location of polymers on Ub, conjugation with a 5 unit chain of DMAm, was used. This DMAm polymer gives peptide conjugates that ionize efficiently in MALDI-TOF-MS after trypsin digestion (Fig. S4†). Peaks resulting from polymer-modified species had masses differing by 99 m/z and isotopic signatures from natural-abundance versus15N-containing species (Fig. S4b†). Masses in Table S2† are the highest-intensity peak in each MALDI-TOF-MS series. Predicted mass matches are within 0.25 m/z of the observed mass. The MALDI-TOF-MS data indicate that conjugation can occur at K11, K48, and K63, however it is possible that not every protein will have precisely 3 conjugations. Instead there could be some variation where some proteins have a smaller number of attached polymers, still at one or two of the residues: K11, K48, and K63. Minor peaks observed in the MS data correspond to conjugation at K33, K6, or the N-terminus, indicating that a small number of proteins could also have additional polymers attached. The MALDI-TOF-MS data agree with a model developed by Russell et al.,37 which predicts K11, K48, and K63 to be fast-reacting, while all other amines are slower.
1H/15N-heteronuclear single quantum coherence (HSQC) spectroscopy was performed on all conjugates and native Ub. This NMR technique gives the proton and nitrogen chemical shift for each N–H bonded group in the protein, averaged over all molecules and conformations. Fig. S5† compares the 1H/15N-HSQC NMR spectrum of the native ubiquitin in aqueous media, to the protein in the presence of 40 equivalents of the B-block spin labeled polymer. The data in Fig. S5† indicate that although there are some small changes in the 1H/15N-HSQC spectrum of ubiquitin, the spectrum is overall preserved even in the presence of large excesses of spin labeled polymer. Therefore, unattached spin labelled polymers in solution only cause minimal perturbations in the protein's NMR spectrum, indicating that the polymer in solution does not have any specific affinity for the protein.
Fig. 1b overlays the NMR spectra of Ub conjugated with the 0-block polymer and native Ub in the presence of 40 equivalents of B-block polymer. Although the 1H/15N-HSQC of the 0-block Ub conjugate exhibits some differences compared to native Ub, due to heterogeneity in conjugation density, chain lengths, and conjugation sites, the similarity between the spectra suggests that perturbations in Ub structure caused by polymer conjugation are minor. However, due to the lack of spin label on the conjugated 0-block synthetic polymer, the data in Fig. 1b do not provide information on potential interactions, proximities, or conformations for the attached polymer and the protein.
To ascertain the possible conformation of the polymer with respect to the protein in the bioconjugates, PRE of protein residues by spin-labeled polymer was probed. If the spin label is within 15 Å of an amide proton, almost complete broadening and signal loss is observed for the corresponding resonance in the 1H/15N-HSQC. When the spin label is at a distance of 15–30 Å, partial broadening and reduced signal intensity occurs. If the spin label is greater than 30 Å from an amide proton, no change in signal intensity is observed (Scheme 1B).31 In this way, the proximities between the spin labeled polymer and the protein's individual amino acid residues can be evaluated, averaged over all bioconjugates in solution. However, it is important to note that substantial perturbations in the protein's NMR spectrum through PRE require the spin labeled polymer to have a dominant conformation in close proximity to the protein's surface. Fig. S5† shows that unconjugated polymers do not lead to notable PREs, suggesting that the poly(DMAm) in solution, without covalent bonding to the protein, is unlikely to have specific affinities or conformations on the protein's surface.
In A-block conjugate, the TEMPO is estimated to be within the total-signal-loss distance rTEMPO = 13 ± 3 Å (Fig. 2) for at least some residues (even if the polymer behaves in its solution form as an unperturbed random coil), and most residues are completely broadened, with only four protein peaks present. B-Block conjugate puts the TEMPO further from the protein with the expected distance in the weak PRE range with rTEMPO = 21 ± 4 Å, however, 47% of residues are completely broadened. Finally, 31% of C-block residues are completely broadened even though the TEMPO is expected to be rTEMPO = 29 ± 6 Å away from the nearest protein residue – almost double the distance that could enable complete PRE broadening. These data suggest the polymer, on average, adopts conformations that place TEMPO near the protein surface, much closer than solution behavior of DMAm would predict. NMR gives information of the ensemble averaged environment experienced by the NMR active nuclei. In this way, the 1H/15N HSQC NMR allows the environment of each of the protein's residue to be probed. Most of the residues experiencing PRE, which is attributed to close proximity of the NMR active residue to the polymer's spin label, located in the β-sheet region. This is most clear in B-block (Fig. 2), where there is a band of PRE-relaxed residues across the face of the β-sheet, but the helical region is mostly unaffected. The β-sheet region harbors surface-exposed hydrophobic and hydrogen-bond donating residues, providing complementary functional groups to those in DMAm. This information, combined with the gradual decrease in the number of PRE-relaxed residues moving from A-block to B-block to C-block leads to the following model: DMAm in the first block of the copolymer lays across or near the beta-sheet face of Ub, where its carbonyl groups may interact with hydrogen-bond donors and shield nonpolar residues from water. The end of the chain is still free to coil as normal. This puts TEMPO closer to the surface than expected, and explains the decrease in the number of PRE-relaxed resonances in C-block. If the polymer was not adopting conformations biased closer to the protein surface, only very weak or negligible PRE effects would be observed in C-block. In contrast, if the polymer was adopting a wider array of conformations that allowed for potential interactions across the entire Ub surface and were not biased to a particular face of Ub, nearly complete PRE would be expected for all residues of ubiquitin and there would likely be more substantial chemical-shift perturbations in the zero-block spectrum (Fig. 1).
To evaluate the impact of conjugation, stability against denaturation by guanidine hydrochloride (G-HCl) was performed (Fig. 3). Stability was monitored by circular dichroism (CD) on native Ub, Ub modified with 2–3 C-block polymers and Ub highly modified with 4–6 C-block polymers (SDS-PAGE in Fig. S3†). The number of polymers is estimated based on predicted reactivities37 and conjugation sites observed by MALDI-TOF (Table S2†). CD gives the mean density of secondary structural elements (α-helices, β-strands, and coils) in a given system. Monitoring the CD signal as a function of a chemical denaturant enables the averaged structural stability to be evaluated as the protein or conjugate faces an increasingly harsh environment. While native ubiquitin is completely denatured at 6 M G-HCl, both bioconjugates contain folded populations even at 7.2 M G-HCl. Based on a 4-parameter logistic symmetric sigmoidal fit, the inflection point for native ubiquitin is 4.0 M; C-block modified conjugate is 4.4 M; and C-block highly modified conjugate is >9.3 M. We attribute stabilization to protein–polymer interactions resulting from the conjugated polymer adopting conformations biasing its position to be near the Ub β-sheet face, conformations consistent with the PREs, occurring predominantly at residues located in the Ub β-sheet regions. Fig. S6† compares the 1H/15N-HSQC of C-block modified and highly modified species. The highly modified conjugate has many more PREs (Fig. S6†), suggesting that an increase in potential polymer–protein interactions are correlated with enhanced protein stability.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d0sc02200d |
‡ Equal contribution. |
This journal is © The Royal Society of Chemistry 2020 |