Self-Assembly of Human Profilin-1 Detected by ... - ACS Publications

Jan 4, 2017 - In particular, Carr–Purcell–Meiboom–Gill (CPMG) NMR experiments(17, 18) have been used to characterize conformationally excited st...
0 downloads 0 Views 5MB Size
Subscriber access provided by UNIV OF CAMBRIDGE

Article

Self-Assembly of Human Profilin-1 Detected by CPMG NMR Spectroscopy Enrico Rennella, Ashok Sekhar, and Lewis E. Kay Biochemistry, Just Accepted Manuscript • DOI: 10.1021/acs.biochem.6b01263 • Publication Date (Web): 04 Jan 2017 Downloaded from http://pubs.acs.org on January 7, 2017

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 38

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

Self-Assembly of Human Profilin-1 Detected by CPMG NMR Spectroscopy Enrico Rennella1, Ashok Sekhar1, and Lewis E. Kay1,2

1

Departments of Molecular Genetics, Biochemistry and Chemistry, The University of Toronto, Toronto,

Ontario, M5S 1A8, Canada 2

Hospital for Sick Children, Program in Molecular Structure and Function, 555 University Avenue,

Toronto, Ontario, M5G 1X8, Canada

Email: [email protected]

ACS Paragon Plus Environment

1

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 2 of 38

Abstract Protein oligomerization in the cell has important implications in both health and disease and an understanding of the mechanisms by which proteins can self-associate is, therefore, of critical interest. Initial stages of the oligomerization process can be hard to detect, as they often involve the formation of sparsely populated and transient states that are difficult to characterize by standard biophysical approaches. Using relaxation dispersion NMR spectroscopy we study the oligomerization of human profilin-1, a protein that regulates the polymerization of actin. We show that in solution and at millimolar concentrations profilin-1 is predominantly monomeric. However, fits of concentration dependent relaxation data are consistent with the formation of a higher order oligomer that is generated via a multi-step process. Together with crystallographic data on profilin-2, a homologue of the protein studied here, our results suggest that profilin-1 forms a sparsely populated tetrameric conformer in solution.

ACS Paragon Plus Environment

2

Page 3 of 38

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

Introduction

Self-association is a ubiquitous feature of many proteins1 and it estimated that in E. coli approximately 65% of the total protein content forms homo-oligomers2. An evolutionary advantage of such complexes is that they can be produced from just a single small gene, reducing the impact of errors during translation relative to large proteins that involve only a single polypeptide chain. Moreover, homooligomeric receptors and enzymes can present multiple copies of identical binding sites. These can be arranged in specific topologies, with ligand binding occurring in a cooperative manner2 and with a strength that can vary depending on the oligomeric state, so that affinities can be modulated by protein concentrations3,4. In many cases, a shift in the equilibrium towards homo-oligomerization controls function, for example by activating/inhibiting enzymatic activity5,6 or by triggering/blocking binding to other proteins7. Polymerization of proteins can also lead to misfunction and to disease, especially in cases of neurodegeneration8–10 Several techniques are available for the characterization of protein self-association.11 These include, among others, small-angle X-ray scattering,12 mass spectrometry,13 electron cryomicroscopy14,15 and ultracentrifugation16. However, under conditions where self-association is weak, so that only a small fraction of an oligomeric state is present, the association process becomes difficult to study, and the sparsely populated oligomers so formed are challenging to characterize in detail. In many cases only average properties are extracted that essentially reflect those of the monomeric state (although scattering techniques are sensitive to small fractions of large complexes). Nuclear magnetic resonance (NMR) spectroscopy has emerged as a valuable method for the study of sparsely populated and transiently formed states of macromolecules. In particular, Carr-Purcell-Meiboom-Gill (CPMG) NMR experiments17,18 have been used to characterize conformationally excited states involved in protein folding19–21 in enzyme function,22–25 in allostery26,27 and in oligomerization21,28–31. Herein we apply CPMG methods to investigate the oligomeric states of human profilin-1, PFN1, ACS Paragon Plus Environment

3

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 4 of 38

and the mechanism by which it oligomerizes. Profilins are evolutionary conserved proteins that bind monomeric

actin

with

sub-micromolar

dissociation

constants,32

and

thus

regulate

actin

polymerization.33–35 The profilin fold consists of a 5-stranded antiparallel β-sheet, with 2 helices on one face (N-terminal and C-terminal helices α1 and α4, respectively) and a mixture of a number of small α and β elements on the other face (helices α2 and α3, strands β3 and β4), Fig. 1. Profilins can inhibit growth of actin filaments by sequestering actin monomers,33,36 using an interface for actin binding comprising α3, β4, β5 and the C-terminal part of α4, or promote polymerization by catalyzing the exchange of ADP to ATP in actin37. Several signalling as well as actin-binding proteins associate with profilins through a specific pocket, located between α1 and α4, that binds proline-rich regions,38 Fig. 1. Profilins can also bind to poly-phosphoinositides,39 although the site of interaction remains to be defined. It has been shown that post-translational modifications of profilin can regulate binding in vivo, with the phosphorylation of Tyr129 and Ser138 or Tyr140 increasing affinities for actin40 and poly-Lproline41, respectively.

Figure 1. Ribbon representation of PFN1 shown as an isolated monomer (left, PDB ID 1PFL42) or in the ternary complex with actin (blue) and a proline-rich peptide (cyan) (right, PDB ID 2PAV43). The sidechains of Leu123 and His134 are shown as white spheres.

ACS Paragon Plus Environment

4

Page 5 of 38

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

In addition to their role in actin biology, profilins have also been found to localize to transcriptionally active loci,44 possibly regulating the activity of transcription factors containing proline-rich regions45 and/or playing a role in pre-mRNA processing46. They play an inhibitory role in multiple types of carcinomas47–49 and, recently, mutations in the PFN1 gene were identified in cases of familial amyotrophic lateral sclerosis (ALS)50 where PFN1 aggregates sequester TDP-4351, a key protein in the pathogenesis of ALS52. Profilins analyzed under non-reducing conditions have been found to exist as monomers, dimers and tetramers53–55 and a crystal structure of human profilin-2 (PFN2) shows it to be tetrameric in the crystal lattice56. Given the biological importance of PFN1 it is of considerable interest to elucidate the nature of its oligomeric states in solution and to obtain insight into the oligomerization process in general. Herein, we show that for the millimolar (mM) protein concentrations considered in this study PFN1 is largely monomeric. It does, however, form sparsely populated oligomers that are not observable in NMR spectra but that can be characterized by protein concentration dependent CPMG relaxation dispersion NMR experiments. Key to our analysis has been an approach that separates the effects of chemical exchange from contributions to transverse relaxation rates that arise from viscosity effects and potentially non-specific, transient protein aggregation that often also increase with protein concentration. Together with available structural data on PFN2 our NMR results strongly suggest the formation of a sparsely populated tetrameric PFN1 conformer in solution.

ACS Paragon Plus Environment

5

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 6 of 38

Materials and Methods Protein production. BL21 E. Coli cells were transformed with a pET-29b plasmid containing the PFN1 gene linked to an N-terminal HisTag (GenScript). Cells were grown at 37 °C in M9 medium containing 15

NH4Cl and 12C-glucose (3 g/L) as the sole nitrogen and carbon sources until OD ~ 0.7, then induced

overnight at 25 °C. In the case where 2H,15N PFN1 was produced, the M9 medium was D2O based and 2

H-glucose (3 g/L) was used as the carbon source. Cells were harvested, resuspended in 25 mM

phosphate buffer, 150 mM NaCl, 25 mM imidazole, 2 mM DTT, pH 6.4, and sonicated. Both WT and L123R/H134R PFN1 mutant proteins localized to the soluble fraction. The protein was purified by nickel affinity chromatography, before removing the HisTag by TEV proteolysis overnight in 25 mM HEPES buffer, 150 mM NaCl, 1 mM EDTA, 2 mM DTT, pH 7. The final purification step consisted of gel-filtration using a Superdex 75 column (GE Healthcare Life Sciences) in 12 mM phosphate buffer, 150 mM NaCl, 2 mM DTT, pH 6.4. NMR samples (ranging in concentration from 0.5 – 3 mM) were dissolved in buffer comprised of 50 mM BisTris, 5 mM TCEP, 1 mM EDTA, pH 6.4, 10% D2O.

NMR experiments. All experiments were recorded at 10 °C, using either a 14.1 T Bruker AVANCE III HD spectrometer equipped with a cryogenically cooled x,y,z gradient probe (15N-PFN1, in-phase CPMG recorded at 5 protein concentrations) or both a 14.1 T spectrometer (Bruker) and an 18.8 T Varian spectrometer, the latter equipped with a room temperature probe (2H,15N-PFN1, in-phase, TROSY, anti-TROSY data sets at 2 protein concentrations). In-phase, TROSY (TR) and anti-TROSY (AT)

15

N CPMG experiments were acquired as previously described.58,63 Constant-time relaxation

intervals57, Trelax, varied between 20-40 ms, where shorter Trelax values (20 ms) were used for the most concentrated WT samples to compensate for the large R2 rates. 26 νCPMG values up to 1 kHz were used in all the experiments, with the CPMG refocusing pulses applied at a field strength of γB1/2π = 5.8 kHz. For the in-phase experiment, 1H continuous wave decoupling was applied during Trelax using a γB1/2π = 16.7 kHz field. To keep the amount of heating constant in these experiments the 1H continuous wave ACS Paragon Plus Environment

6

Page 7 of 38

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

element was applied during the recycle delay for reference experiments recorded with Trelax = 0 s. A recycle delay of 2.5 s was used in all CPMG experiments. Lorentz-to-Gauss transformation was applied to both direct and indirect acquisition dimensions before Fourier transform, and signal intensities were quantified using the multidimensional lineshape fitting routine nlinLS in the NMRPipe software package73. In total, relaxation dispersion data recorded on wild-type PFN1 required 13.5 days of acquisition time corresponding to (66 h, 33 h, 17 h, 17 h, 17 h) for in-phase 15N data (14.1 T) obtained at protein concentrations of (0.5 mM, 1 mM, 1.3 mM, 2 mM, 3 mM); (29 h, 57 h) for in-phase and TR/AT profiles, respectively, for both protein concentrations of 1 mM and 1.8 mM (18.8T). In addition, 6 in-phase

15

N dispersion profiles were recorded in total for the two mutants (L123R and H134R

PFN1), each of duration 17 hours, for a total of 4.2 days. Data fitting. CPMG data for 6 selected residues of wild-type PFN1 (those with ∆R2,eff > 2 s-1 at [PFN1] = 1 mM and free of spectral overlap, corresponding to Gly3, Gly121, Leu123, Ile124, Asn125 and Arg137) were fit to a model of chemical exchange via the Bloch-McConnell equations64,  = − 





where the explicit elements of the Liouvillian 

[1]



depend on the model of exchange and the

subscripts K and J are the number of oligomeric states and the size of the Liouvillian (JxJ) that is

required for each exchanging state, l, respectively. When K = 2 or 3,  ⋅ =



0 [2] 

 ⋅ =  0 [3]

0

0 + −  0 

0

0

0 + 

− − −

⊗  , +

− −

(K=2)

+

− −

+

⊗  ,

(K=3)

where  is the J-dimensional identity matrix. For the in-phase CPMG experiment (J=3) ACS Paragon Plus Environment

7

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

  = −



0 0



0

Page 8 of 38

0

[4]



so that the explicit Liouvillian describing exchange between two states (K=2) is given by

 ⋅

+    −  0 = −  0   0 [5]

where



,



0 0  + 0 0 −

+ 0 0 − 0



and



0 0  + − 0

− 

0

0

+ 0

− 

0 0

0 0 +

   ,   

are longitudinal and transverse spin relaxation rates and the chemical shift

(rad/sec) for a nucleus in state l, and  = [NAx, NAy, NAz, NBx, NBy, NBz]+ (see Eq [1]) where the

superscript ‘+’ denotes transpose.

For TROSY and anti-TROSY based CPMG experiments the equations are slightly more complicated because both in-phase and anti-phase operators are needed to describe the evolution of magnetization for each oligomeric state l (J=6). Thus Eq. [4] is replaced by74 ,  −   0  =    − ⋅  0 



[6]

where



and



,

0



0

0 0



,

− ⋅ 0



0

0 0



,





0

,

0

0  0    , 0  0   , 

are longitudinal and transverse cross-correlated relaxation rates for an 15N spin of an

amide pair in state l, JHN is the one-bond 1H-15N scalar coupling constant and the i/a subscripts denote in-phase / anti-phase coherences. Thus, the appropriate Liouvillian for exchange between 2 and 3 states

is  ⋅ (12x12 matrix) and  ⋅ , (18x18 matrix) which can be constructed from Eqs. [2] and [3], respectively, with the sub-Liouvillian for each state l as defined by Eq. [6]. In the case of exchange

ACS Paragon Plus Environment

8

Page 9 of 38

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

between two states (J=6) the magnetization vector for each state l is given by l = [Nlx, Nly, Nlz, 2HzNlx,

2HzNly, 2HzNlz]+ and  = [NAx, NAy, NAz, 2HzNAx, 2HzNAy, 2HzNAz, NBx, NBy, NBz, 2HzNBx, 2HzNBy, 2HzNBz ]+ where subscripts A, B denote the pair of exchanging states.

Because concentration dependent CPMG relation dispersion profiles are fit that depend on the overall tumbling of the protein (that, in turn, increases with protein concentration due to the increase of solvent viscosity and, possibly, due to non-specific protein-protein interactions) intrinsic transverses relaxation rate , 

!

1#$ =

, 

in all states l were corrected using the relation 0

$⋅

,

[7]

are obtained from linear fits of the R2 vs [PFN1] profiles for each

where the per-residue factors

residue in the L123R mutant (Fig. S4). Values of

are equal to one at zero protein concentration

and range between 1.4-1.8, depending on internal protein dynamics, at [PFN1] = 3 mM.



rates for

the monomeric state were measured in a separate experiment and were not corrected for concentration in fits of CPMG data. Eq. [7] assumes that transverse relaxation rates are proportional to the overall tumbling time, which is a good approximation for a 15 kDa protein at 10 oC. In simultaneous fits of TROSY, anti-TROSY and in-phase

perdeuterated samples, all transverse relaxation terms, calculated from the product,

,



,

15

N CPMG profiles recorded on

, , and

 for

each state, l, were

$, which was a floating variable in the fits (S2,l is the square of

the backbone order parameter in state l and

is the correlation time in that state). Additional

contributions to relaxation from 1H spin flips RHl z ,ext were included for anti-phase relaxation terms, ,

to

and ,



,

,

=

$ + RHl z ,ext where RHl z ,ext was assumed to be linearly proportional

$. As in the analysis of 15N in-phase CPMG data sets, the change in sample viscosity with

protein concentration was included using equations analogous to Eq. [7].

Different models of chemical exchange were considered, as indicated in Table 1. Exchange data

ACS Paragon Plus Environment

9

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

were fit to pseudo 2-state or 3-state models,



Page 10 of 38



or [

,





,

#, respectively,

and the extracted rates kij, i,j ∈ {A,B,C} converted into k±1,k±2 via relations in Table 1 (see Supporting Information for further details). It is, of course, not surprising that fits to multi-step models would be more time consuming than for those based on two-state exchange. In Fig. S6 it is shown that robust estimates of exchange parameters and chemical shift differences can be readily obtained from fits of simulated data generated k1   → A model. In fits to multi-step models, such as 2 using the 4 A ←  4 k −1

⇌ '

 ,2

('





)

()

*,

we have

initially used a grid-search focusing on dispersion data for G121 that has the largest of the experimental profiles examined. A 3-dimensional space (K1,k-1,k-2) was explored, with K1=k1/k-1; k-1 and k-2 were sampled logarithmically from 10 to 108 s-1, while K1 was varied linearly between 0.4 and 20.4 M-1 (that 2 gives p A values ranging from 0.2% to 10% for a monomer-dimer model). All fits for which χ red < 1.5 2

2 were clustered. Parameters corresponding to the lowest χ red were, therefore, used as starting values for

a final optimization where the four microscopic rates and chemical shift values were allowed to float.

Table 1: Chemical exchange models considered in the analysis of PFN1 CPMG profiles. Model

[2]

=



[2]

= 2⋅



⋅ ! #,

*

[2]

= 4⋅



⋅ ! # ,

Eqn.

⇌ '

('

(k1 and k-1 are allowed to ‘float’ with concentration) 2

⇌ '

('

3 2



'

('

4

⇌ '

('

⇌ '

,

('



+



[2]

⇌ )

()



[3]

Rates

=

,

= 3⋅ = 2⋅





⋅ ! # , ⋅ ! #,

+

=

= =

=

+

+ +

+ ,

=





ACS Paragon Plus Environment

10

Page 11 of 38

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

2



'

('

 ,2



⇌ )

()

! #, *

[3]

+

= ⋅  

+ ,

= 2 ⋅  ⋅ ! #, !  #, = +,

=

= = 0,



⋅!

+ ,

 #,

=0

= ⋅

=2⋅

 





All fits were carried out using the chemex package (https://github.com/gbouvignies/chemex), minimizing a χ2 target function, as described previously.75 Specific chemex modules for the various oligomerization models considered in this paper are available upon request.

Proline-rich peptide. A peptide, GAGGGPPPAPPLPAAQ (MW = 1354.54), corresponding to the proline-rich PFN1 binding site of human VASP43, was purchased from ThermoFisher. A complex of 2.0 mM WT PFN1 was generated with 3.7 mM peptide; based on the published affinity constant (≤ 84 1kHz µM43), the estimated amount of PFN1 bound to the peptide is over 95%. Expected R2,eff values for the

1:1 complex have been extrapolated from rates measured for the L123R mutant, that take into account viscosity effects, Eq. [7]. Values obtained in this manner were further multiplied by a factor of 1.09 = ( MWPFN1 + MWpeptide ) / MWPFN1, to account for the size increase upon binding.

Results Experimental Approach In what follows we have used CPMG-based experiments to study the oligomerization process of PFN1. In this approach signal intensities derived from the ground state of the protein are measured as a function of the repetition rate of chemical shift refocusing pulses, νCPMG, that are applied during a delay of constant duration, Trelax. These intensities, I(νCPMG), are subsequently converted to effective transverse relaxation rates, R2,eff = (1/Trelax) ln[I0/I(νCPMG)], where I0 is the reference signal intensity when Trelax = 0 s57. Random fluctuations between exchanging states on the millisecond (ms) timescale lead to a stochastic modulation of chemical shifts and for small νCPMG values (i.e., only a few refocusing pulses are applied) dephasing of the magnetization occurs, resulting in a reduction of the net ACS Paragon Plus Environment

11

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 38

signal and an increase in R2,eff values. In contrast, when refocusing pulses are applied more frequently (large νCPMG) the fluctuations lead to less dephasing, the signal is thus affected less by the exchange process and R2,eff values become smaller. The resulting dispersion profile, R2,eff vs νCPMG, can be fit to appropriate models to extract the kinetics and thermodynamics of the exchange process as well as the chemical shifts of the interconverting low populated state(s), so long as the effective two-site exchange rate lies between several hundred to several thousand per second and the population of the rare state is 0.5% or larger.17,18 A large variety of CPMG experiments have been developed to measure additional structural parameters in rare states such as bond-vector orientations58–61 as well as amplitudes of bond vector motions62. Taken together the data so obtained facilitates the calculation of atomic-resolution structural models of excited state conformers, such as protein folding and misfolding intermediates19,20,28. To date relatively few studies have, however, focused on oligomerization processes more complex than the monomer-dimer variety, that forms the basis of the present work.

PFN1 forms oligomers on the ms timescale 15

N CPMG relaxation dispersion experiments58,63 have been performed on samples of PFN1 at different

protein concentrations (Fig. 2). An initial analysis of the raw CPMG data can be performed simply by calculating the sizes of the dispersion profiles, ∆R2,eff = R2,eff(νCPMG,low) - R2,eff(νCPMG,high) where

νCPMG,low and νCPMG,high are apparent R2 values recorded at the lowest and highest CPMG frequencies, respectively (~ 30 Hz and 1 kHz in this study). Non-zero values of ∆R2,eff delineate regions of the protein where exchange processes are occurring.

15

N ∆R2,eff values for PFN1 acquired at different

protein concentrations indicate the presence of at least two distinct conformational exchange processes occurring on the ms timescale, Fig. 2A. One of them, involving the loop between strands β4 and β5, corresponding to residues Asp76-Met86, does not depend on protein concentration in a significant way and hence is not related to protein oligomerization. It will not be discussed further. The second process is localized mainly to the C-terminal helix α4, but also includes the N-terminal helix α1 and Glu83 in the ACS Paragon Plus Environment

12

Page 13 of 38

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

loop between strands β4 and β5. 15N ∆R2,eff values for this process depend on protein concentration, Fig. 2A, consistent with self-assembly of PFN1.

Figure 2. Analysis of in-phase concentrations, 10oC.

15

N CPMG data acquired on a sample of

15

N WT PFN1 at several protein

30 Hz 1kHz 1kHz ∆R2,eff = R2,eff − R2,eff and R2,eff are plotted at two different concentrations in panels A

and B, respectively. The position of Glu83, the only residue in loop β4-β5 that clearly displays concentrationdependent ∆R2,eff values, is indicated by the dashed vertical line in panel A. (C) 15N-CPMG profiles for Leu123 are shown at three different protein concentrations, clearly establishing that both the baseline and the sizes of the dispersion profiles increase with concentration. (D) Crystal structure of the tetrameric state of human profilin-2 (PDB ID: 1D1J56). Regions colored in red in the yellow chain indicate residues with ∆R2,eff > 1 s-1 at 1 mM ACS Paragon Plus Environment

13

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 38

protein concentration, 600 MHz; the backbone 15N atoms of Leu123 and His134 are shown as white spheres.

Fig. 2B shows a plot of apparent R2 values measured at the highest νCPMG frequency of 1 kHz, 1kHz 1kHz , where the exchange process of interest is quenched (i.e., flat dispersions, see below). R2,eff rates R2,eff

are sensitive to the rotational diffusion of the protein and their increase with concentration is thus consistent with an oligomerization process. However, the effective solution viscosity also grows with 1kHz protein concentration, leading to increased R2,eff values, and it is not possible to separate these two

effects from this data alone. A representative

15

N CPMG profile acquired as a function of protein

concentration is shown for Leu123, Fig. 2C, clearly establishing that both the size of the dispersions and the apparent R2 values at 1 kHz increase as a function of protein concentration. Residues with concentration dependent CPMG profiles (red) are localized to regions defined as points of contact between monomers in a previously published tetrameric crystal structure of human PFN256, Fig. 2D, a protein that is a close homologue to PFN1 (sequence identity of 62%). Notably, helix 1 and the Cterminus of helix 4 (yellow molecule, Fig. 2D) form one of the interfaces with an adjacent molecule (blue), while the N-terminus of helix 4 and the loop between strands β4 and β5 contact a second molecule (green). The fact that both isoforms of the same protein share similar oligomerization properties may reflect a functional role for the self-assembly of profilins.

Mutation of either L123 or H134 inhibits oligomerization of PFN1 In order to gain insight into the nature of the concentration dependent exchange process that is quantified by our CPMG relaxation dispersion experiments we exploited the available X-ray structure of the PFN2 tetramer showing that Leu123 and His134 sit in two distinct dimer-dimer interfaces of the tetramer56 (Fig, 2D, Fig. S1). These residues were modified (one at a time) to Arg since the burial of positively charged Arg residues in the mutant structures was predicted to be unfavourable, potentially stabilizing the monomeric structure. Both L123R and H134R mutations have little effect on the tertiary

ACS Paragon Plus Environment

14

Page 15 of 38

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

structure of the protein as judged by 1H,15N HSQC spectra (Fig. S2). 15N CPMG experiments acquired for L123R-PFN1 (H134R-PFN1) at a number of protein concentrations have been analyzed and, in contrast to the wild-type (WT) protein, ∆R2,eff rates in the terminal helices are close to zero and independent of concentration (Fig. 3 and Fig. S3). Values of ∆R2,eff ~ 0 s-1 for both mutants very likely reflect an impaired ability to oligomerize or, less likely, a shift of the exchange timescale outside the CPMG window (i.e., exchange lifetime not within the ms range). If a decrease in the kinetics of exchange is the reason for the lack of concentration dependent dispersions then it might be expected that additional resonances would be present in spectra, yet these are not observed. Conversely, an increase in kinetics that would place the exchange rate in the fast regime would lead to small shifts in peak positions, that again is not consistent with experiment. Although we cannot unequivocally exclude either of these two possibilities, we can estimate upper limits for the residual populations of the oligomers, pB, in each of the two cases. Assuming that peaks at a threshold of three times the noise can be detected in high sensitivity 1H,15N HSQC spectra an upper bound for pB of 1.5% is calculated if slow exchange between oligomeric states is operative. Conversely, assuming that a shift of 0.05 ppm (15N) can be detected then pB is estimated to be less than 5% in the fast exchange regime. These bounds hold for both mutants at a protein concentration of 2.5 mM.

ACS Paragon Plus Environment

15

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 3. Analysis of in-phase

15

Page 16 of 38

N CPMG dispersion profiles for U-15N L123R and H134R PFN1 at several

1kHz protein concentrations, 10oC, as for WT PFN1 in Figure 2. (A,B) ∆R2,eff and R2,eff for L123R PFN1, plotted at

two protein concentrations, with 15N CPMG profiles for Arg123 of L123R PFN1 at three concentrations shown in 1kHz panel C. (D) Effect of protein concentration on average R2 values ( R2,eff ) for the most rigid regions of the

protein (85 residues); the standard errors are less than 0.25 s-1 for wild-type PFN1 and less than 0.15 s-1 for both mutants.

1kHz Notably, and in a manner analogous to the WT protein, R2,eff values for amide 15N nuclei from

ACS Paragon Plus Environment

16

Page 17 of 38

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

either of the two mutants increase with protein concentration (Fig. 3B), as do the CPMG profiles that remain flat over the entire νCPMG range, illustrated for Arg123 of L123R-PFN1 in Fig. 3C. However, the 1kHz increase in R2,eff as a function of protein concentration for both L123R and H134R mutants is less than

what is observed for WT PFN1, Fig. 3D, reflecting the fact that the oligomerization process observed for the wild-type protein is inhibited in both mutants. As shown in Fig. 3D, the relationship between 1kHz and [PFN1] is linear for both variants of PFN1 considered, where rates have been averaged for R2,eff

approximately 80 residues that are associated with rigid portions of the protein. For these residues 1kHz rates report primarily on the overall tumbling of the molecule that is affected by solvent viscosity R2,eff

and, possibly, by non-specific protein-protein interactions. As expected, the slopes are smaller in the case of residues that reside in flexible regions where the amplitudes of amide bond vector motions increase (Fig. S4). Oligomerization of PFN1 proceeds through multiple steps The concentration dependent CPMG profiles from 6 residues of WT PFN1 (largest dispersions, see Materials and Methods), were globally fit to the Bloch-McConnell equations that describe the exchange of magnetization during the CPMG interval64. In addition to the ‘traditional’

15

N-based experiments

where in-phase magnetization profiles are measured63 (5 protein concentrations ranging from 0.5 to 3.0 mM), we also recorded both TROSY- and anti-TROSY CPMG experiments58 (2 protein concentrations, 1.0 and 1.8 mM) where the coherences of interest are slowly and rapidly relaxing

15

N elements,

respectively, that correspond to nitrogen transverse magnetization coupled to β and α amide proton spin states. Unlike the in-phase data set, where protonated samples were used for measurements, a perdeuterated sample was employed in the TROSY-/anti-TROSY-based experiments to minimize 1H spin flips. Initially, we considered a simple kinetic model,



, where the populated and sparse

states, A and B, respectively, interconvert with an effective exchange rate kex = kAB+kBA. We recognize

ACS Paragon Plus Environment

17

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 38

that this simplest of models is in some sense a ‘virtual scheme’ because it does not explicitly take into account oligomerization, yet as we show below, it provides a starting framework by which to begin to make choices about different oligomerization models. The interconversion process could be well fit by this pseudo first-order exchange model (reduced χ2red ~ 1), Figs. 4A-C, that also included the contribution to the intrinsic R2 rates of the exchanging states from the viscosity vs protein concentration dependence (see Materials and Methods). Values of kAB and kBA were allowed to ‘float’ with concentration. Notably, kAB and pB, increase with [PFN1] in a reasonably linear manner, while kBA is essentially independent of protein concentration. Such a dependence would be expected for a dimerization reaction, yet in that case kAB(pB) would extrapolate to 0 s-1(0%) for [PFN1]=0 mM and this is not what is observed. Indeed, the utility of the pseudo first-order model considered here is that it shows clearly that a dimerization model cannot explain the concentration dependence of the dispersion profiles. However, it provides no physical insight into what the oligomerization process might be. This can only be obtained by considering more complex models, as described below.

Figure 4. Fits of concentration dependent 15N CPMG datasets for WT PFN1 (10oC) using a pseudo first-order

ACS Paragon Plus Environment

18

Page 19 of 38

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

model of exchange. The CPMG profiles for three residues are shown in panels (A-C), with experimental data (best global fit) indicated by points (solid lines). In (A) and (B) dispersion profiles are shown for Gly3 and Leu123 at [PFN1] = 3, 2, 1.3, 1 and 0.5 mM, from the top to bottom (light gray to black), while in (C) TROSY- (TR, light-gray at the bottom), in-phase (CW, dark-gray in the middle) and anti-TROSY (AT, black at the top) 15N CPMG profiles are illustrated for Asn125, recorded on a U-[2H,15N]-PFN1 sample at [PFN1] = 1 mM (closed circles) and 1.8 mM (open diamonds). Fitted rates and populations are plotted in panels (D-F) as a function [PFN1]. The solid lines in (D-F) connect the fitted parameters as a function of concentration and are intended to guide the eye.

Having established that a model with higher order kinetics is necessary to describe the CPMG data, we initially considered several simple oligomerization schemes involving the interconversion between two states. These include a monomer-dimer (Fig. 5A) scheme where the relatively poor fit (χ2reduced = 3.3) confirms the results of the simple



model, discussed above. We have also

examined more complex models such as those involving monomer-trimer (Fig. S5) and monomertetramer (Fig. 5B) interconversions. These models were, however, not able to explain the recorded CPMG profiles, even in cases where the effect of viscosity on the microscopic rates was included (not shown). It is worth emphasizing that the poor fits are not a reflection of lack of convergence. Figure S6 shows correlations of extracted exchange parameters from fits of simulated data generated from the two k1   → A model, where excellent agreement is obtained with the input parameters. In contrast state 4 A ←  4 k −1

to fits involving two-state models, a more complex model involving consecutive oligomerization steps, !2



'

('

 ,2



⇌ )

()

* #,

Fig. 5C, was able to fit the data reasonably well (χ2reduced = 1.6) with the

fractional populations of A2 and A4, pA and pA , decreasing to 0% in the limit that [PFN1] = 0 mM, as 2

4

would be expected. In these fits the value of k-1 was reasonably well defined (~ 500 s-1, Fig. S7A) as were the equilibrium constants and hence populations (Fig. S7B), but other rate constants were found to be correlated and thus not well determined. Other models also gave reasonable fits, such as a scheme

ACS Paragon Plus Environment

19

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 38

where dimers combine with monomers to form trimers (Fig. S5). Thus, while it is not possible to unambiguously determine the underlying kinetic model, it is clear that schemes with more than a single reaction step are required to fit the data, strongly supporting the notion that oligomerization of PFN1 proceeds through a series of consecutive binding events.

Figure 5. (A-C) Fits of CPMG data for WT PFN1, 10oC, using several models of oligomerization, as indicated. The models used to fit the data are shown on the left side, along with the reduced χ2. Shown also are the populations of states A2 and A4 as a function of [PFN1] (defined as p A = n[ An ] / [PT ] where PT is the total

protein concentration) derived from the fit in (C), !2

⇌ '

('

n

 ,2





)

()

*# ,

along with the errors in the

populations (shaded region), as obtained via a boot-strap analysis65,66. Concentrations of PFN1 used in the CPMG analysis are as in the legend to Figure 4.

Binding of poly-proline substrate shifts equilibrium to monomer The X-ray structure of PFN2 shows that the surface sequestered upon oligomerization overlaps with the ACS Paragon Plus Environment

20

Page 21 of 38

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

binding pockets for both actin and poly-proline and the µM binding affinities that have been reported previously67 indicate that binding of either partner would thus prevent oligomer formation. If a similar oligomeric structure to that observed in the solid state for PFN2 is adopted by PFN1 in solution then binding of a proline rich peptide to PFN1 would be expected to shift the equilibrium to the monomeric state, that would lead to quenching of exchange contributions to transverse relaxation rates. Fig. 6A plots ∆R2,eff values for WT PFN1 (2.0 mM) measured in the absence (black) and in the presence (grey) of a two-fold excess of proline-rich peptide. Binding of the peptide reduces ∆R2,eff values in the N- and C-terminal helices, that contain reporter residues of oligomerization, to close to 0 s-1, consistent with 1kHz inhibition of PFN1 oligomerization. In addition, R2,eff values for residues in the sample with peptide

1kHz are consistent with PFN1 binding to peptide as a monomer (Fig. 6B). The offset between R2,eff rates for

PFN1 with or without peptide in Fig. 6B reflects the fact that the oligomerization process observed for the wild-type protein is inhibited in the presence of the proline-rich substrate.

ACS Paragon Plus Environment

21

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 38

1kHz Figure 6. Addition of a proline-rich peptide eliminates chemical exchange in WT PFN1. ∆R2,eff (A) and R2,eff (B)

are plotted vs residue based on measurements recorded on a 2 mM PFN1 sample with (grey) or without (black) 1kHz 3.7 mM peptide (10oC). Calculated R2,eff values for the 1:1 complex from extrapolation of rates measured for

the L123R mutant to 2 mM protein concentration and taking into account the increased molecular weight of the complex are shown with a solid grey line in (B) (see Materials and Methods). Dashed lines in (A) and (B) are meant to guide the eye.

Discussion Protein oligomerization plays a critically important role in both cellular function and misfunction.1,2,9 For example, many of the vital processes in the cell as related to homeostasis are carried out by oligomeric molecular machines, where the multi-subunit nature of the structure is clearly important functionally.68 Oligomeric conformers are also implicated in disease, in particular in neurodegeneration, where excursions from a populated monomeric state to higher order structures can lead to the formation

ACS Paragon Plus Environment

22

Page 23 of 38

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

of toxic species.10 The continued development of NMR spin-relaxation experiments provides an avenue for investigating oligomerization at atomic resolution, generating mechanistic and physical chemical insights into this important process. Herein we have characterized homo-oligomerization of the actin binding protein PFN1 by analyzing

15

N CPMG relaxation dispersion profiles of the protein as a function of concentration. An

assumption in most analyses of CPMG relaxation dispersion data is that the intrinsic relaxation rates of spins in the interconverting states are identical. This is clearly incorrect in general, and in the present case pathologic. It is also important that the concentration dependence of the intrinsic relaxation rates be taken into account properly. Here we have kept intrinsic rates as fitting parameters, subject to the constraint that they vary with protein concentration on a per residue basis as measured for L123RPFN1 that does not show evidence of chemical exchange (i.e., flat dispersion profiles) but whose tumbling is sensitive to the increased viscosity and potentially to non-specific protein-protein interactions that accompany higher PFN1 protein concentrations (Fig. S4, see Materials and Methods). Figures 3C and D indicate that these effects can be considerable. In the present study we show that PFN1 undergoes a conformational exchange process involving the interconversion between an abundant monomer species and a sparsely populated oligomeric structure in solution. Notably, the dispersion data cannot be well fit to a simple oligomerization model,



, but rather must involve at least a pair of association steps. 15N in-

phase, TROSY, and anti-TROSY CPMG data can be reasonably well fit to a monomer-dimer-tetramer model although this is not the only scheme that is consistent with the dispersion profiles. We have attempted to use our concentration dependent NMR spin relaxation measurements to provide an estimate of the effective increase in particle size upon oligomerization in a model independent manner. Using a simple



scheme, average values of R2 for state B can be obtained

(A is the monomer) that, in principle, might provide insight into the size of B. Detailed simulations

ACS Paragon Plus Environment

23

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 38

establish that such an approach works well for simple interconversion processes of the form,



,

but not for more complex, multi-process schemes. For example, consider a 3-step process such as ⇌



, where B and C are sparse states. Since dispersion profiles are critically sensitive to

chemical shift differences between spins in each of the states and if the majority of these happen to be between A and B then effective R2 values obtained for the rare state

′ from fits to the simple ′ ⇌



scheme are skewed to those in B. Thus, the extracted relaxation rates depend on chemical shift differences, so that, in general, reliable rates cannot be obtained. We have further attempted to use a ‘generic’ 3 state model to obtain information about molecular size but the rates so obtained are not well defined. Despite the lack of ‘proof’ from our experimental data there are, nevertheless, several lines of evidence to suggest that the sparse PFN1 state characterized in the present study is predominantly tetrameric. First, an X-ray structure of a related profilin, PFN2, that shares a high degree of sequence homology to PFN1, is tetrameric in the crystal56 and significant dispersion profiles are observed precisely in regions that form contacts between adjacent monomers in the solid state, including terminal helices α1 and α4 and Glu83 in the loop connecting β4 with β5. The fact that other residues within this loop do not show significant concentration dependent dispersion profiles may reflect the fact that only small chemical shift differences arise from the oligomerization process in this region and/or that the concentration independent exchange event in this region (Fig. 3) masks the one of interest in this study. We have shown that in solution the oligomeric structures can be disrupted by either of two mutants, L123R and H134R, which are at the interface between the monomers in the Xray model, supporting the notion that the oligomeric states characterized by NMR and X-ray methods are structurally similar. Second, the results from the present study establish that binding of ligands, such as poly-proline rich peptides, compete with PFN1 oligomerization. Dispersion profiles recorded in the presence of saturating amounts of a poly-proline peptide were flat, consistent with a significant

ACS Paragon Plus Environment

24

Page 25 of 38

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

decrease in the population of the sparse oligomeric state. This is, again, in keeping with X-ray studies that establish that binding of such peptides would be expected to disrupt the interface between protomers in the reported tetrameric X-ray structure of PFN2 by preventing contacts between the Nterminal helix of one protomer and the C-terminus of the C-terminal helix from an adjacent molecule (Figure S1C). The X-ray structure also predicts that actin binding would lead to a similar result, although we have not tested this here due to the difficulties in working with monomeric actin. The fact that the PFN1 oligomer is only sparsely populated at mM concentrations ( p A ~7% for 4

[PFN1]=3mM) implies low association constants (Fig. S7); assuming a monomer-dimer-tetramer model a value of pA + pA ~ 50% is predicted when [PFN1] = 12 mM. Given the weak association it is 2

4

reasonable to ask whether the tetramer observed in solution is physiologically relevant. Profilins are amongst the most highly abundant proteins in the cytoplasim with an estimated concentration of 20 µM - 100 µM,69 and it is known that the local concentration of this protein is significantly enriched in proximity to proline-rich binding partners70 or in highly dynamic microfilament structures associated with cellular membranes in motile cells71,72. Notably, one calculation suggests that profilins can be concentrated to approximately 1 mM at polyproline sites inside cells.70 In addition, rapid modulation of microfilament architecture in response to external signals is thought to require elevated concentrations of profilins.72 It may be that at sites of high profilin concentration the protein is stored, at least partially, as homo-oligomers. Notably, the potential to oligomerize is not restricted to PFN1 or PFN2. Structures of a pair of additional profilins from plants showed dimerization contacts involving one of the oligomerization interfaces that has been identified in this work, involving helices α1 and α473. In summary, we have demonstrated that relaxation dispersion NMR approaches can be applied to study protein oligomerization and have the potential to discriminate between alternative models for the initial steps of the oligomerization pathway. In combination with results from other studies it is possible to establish the likely oligomerization state of sparse conformers that are formed. This has

ACS Paragon Plus Environment

25

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 26 of 38

important implications in characterizing the energy landscapes of other aggregation prone protein systems, where protein association can lead to the formation of toxic species that have been implicated in a number of debilitating neurological diseases10.

Acknowledgements This work was supported by grants from the CIHR. L.E.K holds a Canada Research Chair in Biochemistry.

Supporting Information Figures of (i) putative PFN1 tetramer and stabilizing interfaces, (ii) HSQC spectra of WT, L123R and 1kHz 1kHz rates for H134R, (iv) R2,eff values as a function of protein H134R PFN1, (iii) ∆R2,eff and R2,eff

concentration, (v) Fits of CPMG data to additional exchange models, (vi) Demonstration of robust fitting of two-state exchange data, (vii) Reduced χ2 surface for k-1 from monomer-dimer-trimer model as well as contour plot of k-1/k1 vs k-2/k2,. A derivation of effective 3-site exchange rates for the monomer-dimer-tetramer exchange model is provided as well. This information is available free of charge via the internet at http://pubs.acs.org.

ACS Paragon Plus Environment

26

Page 27 of 38

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) Ali, M. H., and Imperiali, B. (2005) Protein oligomerization: how and why. Bioorg. Med. Chem. 13, 5013–5020. (2) Goodsell, D. S., and Olson, A. J. (2000) Structural symmetry and protein function. Annu. Rev. Biophys. Biomol. Struct. 29, 105–153. (3) Swanson, C. J., Ritt, M., Wang, W., Lang, M. J., Narayan, A., Tesmer, J. J., Westfall, M., and Sivaramakrishnan, S. (2014) Conserved modular domains team up to latch-open active protein kinase Cα. J. Biol. Chem. 289, 17812–17829. (4) Nooren, I. M. A., and Thornton, J. M. (2003) Structural characterisation and functional significance of transient protein-protein interactions. J. Mol. Biol. 325, 991–1018. (5) Sarabipour, S., and Hristova, K. (2016) Mechanism of FGF receptor dimerization and activation. Nat. Commun. 7, 10262. (6) Brooks, A. J., Dai, W., O’Mara, M. L., Abankwa, D., Chhabra, Y., Pelekanos, R. A., Gardon, O., Tunny, K. A., Blucher, K. M., Morton, C. J., Parker, M. W., Sierecki, E., Gambin, Y., Gomez, G. A., Alexandrov, K., Wilson, I. A., Doxastakis, M., Mark, A. E., and Waters, M. J. (2014) Mechanism of activation of protein kinase JAK2 by the growth hormone receptor. Science 344, 1249783. (7) Diaz, J. L., Oltersdorf, T., Horne, W., McConnell, M., Wilson, G., Weeks, S., Garcia, T., and Fritz, L. C. (1997) A common binding site mediates heterodimerization and homodimerization of Bcl-2 family members. J. Biol. Chem. 272, 11350–11355. (8) Haass, C., and Selkoe, D. J. (2007) Soluble protein oligomers in neurodegeneration: lessons from the Alzheimer’s amyloid beta-peptide. Nat. Rev. Mol. Cell Biol. 8, 101–112. (9) Protein Misfolding, Aggregation and Conformational | Vladimir N. Uversky | Springer. (10) Chiti, F., and Dobson, C. M. (2009) Amyloid formation by globular proteins under native conditions. Nat. Chem. Biol. 5, 15–22. (11) Gell, D. A., Grant, R. P., and Mackay, J. P. (2012) The detection and quantitation of protein oligomerization. Adv. Exp. Med. Biol. 747, 19–41. (12) Blobel, J., Bernadó, P., Svergun, D. I., Tauler, R., and Pons, M. (2009) Low-resolution structures of transient protein-protein complexes using small-angle X-ray scattering. J. Am. Chem. Soc. 131, 4378–4386. (13) Mehmood, S., Allison, T. M., and Robinson, C. V. (2015) Mass spectrometry of protein complexes: from origins to applications. Annu. Rev. Phys. Chem. 66, 453–474. (14) Zhang, H., Griggs, A., Rochet, J.-C., and Stanciu, L. A. (2013) In vitro study of α-synuclein protofibrils by cryo-EM suggests a Cu(2+)-dependent aggregation pathway. Biophys. J. 104, 2706– 2713. (15) Chen, S. W., Drakulic, S., Deas, E., Ouberai, M., Aprile, F. A., Arranz, R., Ness, S., Roodveldt, C., Guilliams, T., De-Genst, E. J., Klenerman, D., Wood, N. W., Knowles, T. P. J., Alfonso, C., Rivas, G., Abramov, A. Y., Valpuesta, J. M., Dobson, C. M., and Cremades, N. (2015) Structural characterization of toxic oligomers that are kinetically trapped during α-synuclein fibril formation. Proc. Natl. Acad. Sci. U. S. A. 112, E1994-2003. (16) Patel, T. R., Winzor, D. J., and Scott, D. J. (2016) Analytical ultracentrifugation: A versatile tool for the characterisation of macromolecular complexes in solution. Methods San Diego Calif 95, 55–61. (17) Palmer, A. G., Kroenke, C. D., and Loria, J. P. (2001) Nuclear magnetic resonance methods for quantifying microsecond-to-millisecond motions in biological macromolecules. Methods Enzymol. 339, 204–238. (18) Korzhnev, D. M., Salvatella, X., Vendruscolo, M., Di Nardo, A. A., Davidson, A. R., Dobson, C. M., and Kay, L. E. (2004) Low-populated folding intermediates of Fyn SH3 characterized by relaxation dispersion NMR. Nature 430, 586–590. (19) Neudecker, P., Robustelli, P., Cavalli, A., Walsh, P., Lundström, P., Zarrine-Afsar, A., Sharpe, S., Vendruscolo, M., and Kay, L. E. (2012) Structure of an intermediate state in protein folding and ACS Paragon Plus Environment

27

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 28 of 38

aggregation. Science 336, 362–366. (20) Korzhnev, D. M., Religa, T. L., Banachewicz, W., Fersht, A. R., and Kay, L. E. (2010) A transient and low-populated protein-folding intermediate at atomic resolution. Science 329, 1312–1316. (21) Hill, R. B., Bracken, C., DeGrado, W. F., and Palmer, A. G. (2000) Molecular Motions and Protein Folding:  Characterization of the Backbone Dynamics and Folding Equilibrium of α2D Using 13C NMR Spin Relaxation. J. Am. Chem. Soc. 122, 11610–11619. (22) Henzler-Wildman, K. A., Thai, V., Lei, M., Ott, M., Wolf-Watz, M., Fenn, T., Pozharski, E., Wilson, M. A., Petsko, G. A., Karplus, M., Hübner, C. G., and Kern, D. (2007) Intrinsic motions along an enzymatic reaction trajectory. Nature 450, 838–844. (23) Boehr, D. D., McElheny, D., Dyson, H. J., and Wright, P. E. (2006) The dynamic energy landscape of dihydrofolate reductase catalysis. Science 313, 1638–1642. (24) Whittier, S. K., Hengge, A. C., and Loria, J. P. (2013) Conformational motions regulate phosphoryl transfer in related protein tyrosine phosphatases. Science 341, 899–903. (25) Ishima, R., Louis, J. M., and Torchia, D. A. (2001) Characterization of two hydrophobic methyl clusters in HIV-1 protease by NMR spin relaxation in solution. J. Mol. Biol. 305, 515–521. (26) Shi, L., and Kay, L. E. (2014) Tracing an allosteric pathway regulating the activity of the HslV protease. Proc. Natl. Acad. Sci. U. S. A. 111, 2140–2145. (27) Tzeng, S.-R., and Kalodimos, C. G. (2013) Allosteric inhibition through suppression of transient conformational states. Nat. Chem. Biol. 9, 462–465. (28) Sekhar, A., Rumfeldt, J. A. O., Broom, H. R., Doyle, C. M., Bouvignies, G., Meiering, E. M., and Kay, L. E. (2015) Thermal fluctuations of immature SOD1 lead to separate folding and misfolding pathways. eLife 4, e07296. (29) Sekhar, A., Bain, A. D., Rumfeldt, J. A. O., Meiering, E. M., and Kay, L. E. (2016) Evolution of magnetization due to asymmetric dimerization: theoretical considerations and application to aberrant oligomers formed by apoSOD1(2SH). Phys. Chem. Chem. Phys. PCCP 18, 5720–5728. (30) Jee, J., Ishima, R., and Gronenborn, A. M. (2008) Characterization of specific protein association by 15N CPMG relaxation dispersion NMR: the GB1(A34F) monomer-dimer equilibrium. J. Phys. Chem. B 112, 6008–6012. (31) Korchuganov, D. S., Nolde, S. B., Reibarkh, M. Y., Orekhov, V. Y., Schulga, A. A., Ermolyuk, Y. S., Kirpichnikov, M. P., and Arseniev, A. S. (2001) NMR study of monomer-dimer equilibrium of barstar in solution. J. Am. Chem. Soc. 123, 2068–2069. (32) Xue, B., and Robinson, R. C. (2013) Guardians of the actin monomer. Eur. J. Cell Biol. 92, 316– 332. (33) Carlsson, L., Nyström, L. E., Sundkvist, I., Markey, F., and Lindberg, U. (1977) Actin polymerizability is influenced by profilin, a low molecular weight protein in non-muscle cells. J. Mol. Biol. 115, 465–483. (34) Karlsson, R., and Lindberg, U. (2007) Profilin, an Essential Control Element for Actin Polymerization, in Actin-Monomer-Binding Proteins, pp 29–44. Springer New York. (35) Witke, W. (2004) The role of profilin complexes in cell motility and other cellular processes. Trends Cell Biol. 14, 461–469. (36) Markey, F., Larsson, H., Weber, K., and Lindberg, U. (1982) Nucleation of actin polymerization from profilactin. Opposite effects of different nuclei. Biochim. Biophys. Acta 704, 43–51. (37) Mockrin, S. C., and Korn, E. D. (1980) Acanthamoeba profilin interacts with G-actin to increase the rate of exchange of actin-bound adenosine 5’-triphosphate. Biochemistry (Mosc.) 19, 5359–5362. (38) Mahoney, N. M., Janmey, P. A., and Almo, S. C. (1997) Structure of the profilin-poly-L-proline complex involved in morphogenesis and cytoskeletal regulation. Nat. Struct. Biol. 4, 953–960. (39) Goldschmidt-Clermont, P. J., Machesky, L. M., Baldassare, J. J., and Pollard, T. D. (1990) The actin-binding protein profilin binds to PIP2 and inhibits its hydrolysis by phospholipase C. Science 247, 1575–1578. ACS Paragon Plus Environment

28

Page 29 of 38

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

(40) Fan, Y., Arif, A., Gong, Y., Jia, J., Eswarappa, S. M., Willard, B., Horowitz, A., Graham, L. M., Penn, M. S., and Fox, P. L. (2012) Stimulus-dependent phosphorylation of profilin-1 in angiogenesis. Nat. Cell Biol. 14, 1046–1056. (41) Sathish, K., Padma, B., Munugalavadla, V., Bhargavi, V., Radhika, K. V. N., Wasia, R., Sairam, M., and Singh, S. S. (2004) Phosphorylation of profilin regulates its interaction with actin and poly (Lproline). Cell. Signal. 16, 589–596. (42) Metzler, W. J., Farmer, B. T., Constantine, K. L., Friedrichs, M. S., Lavoie, T., and Mueller, L. (1995) Refined solution structure of human profilin I. Protein Sci. Publ. Protein Soc. 4, 450–459. (43) Ferron, F., Rebowski, G., Lee, S. H., and Dominguez, R. (2007) Structural basis for the recruitment of profilin-actin complexes during filament elongation by Ena/VASP. EMBO J. 26, 4597– 4606. (44) Söderberg, E., Hessle, V., von Euler, A., and Visa, N. (2012) Profilin is associated with transcriptionally active genes. Nucl. Austin Tex 3, 290–299. (45) Lederer, M., Jockusch, B. M., and Rothkegel, M. (2005) Profilin regulates the activity of p42POP, a novel Myb-related transcription factor. J. Cell Sci. 118, 331–341. (46) Skare, P., Kreivi, J.-P., Bergström, A., and Karlsson, R. (2003) Profilin I colocalizes with speckles and Cajal bodies: a possible role in pre-mRNA splicing. Exp. Cell Res. 286, 12–21. (47) Yao, W., Ji, S., Qin, Y., Yang, J., Xu, J., Zhang, B., Xu, W., Liu, J., Shi, S., Liu, L., Liu, C., Long, J., Ni, Q., Li, M., and Yu, X. (2014) Profilin-1 suppresses tumorigenicity in pancreatic cancer through regulation of the SIRT3-HIF1α axis. Mol. Cancer 13, 187. (48) Janke, J., Schlüter, K., Jandrig, B., Theile, M., Kölble, K., Arnold, W., Grinstein, E., Schwartz, A., Estevéz-Schwarz, L., Schlag, P. M., Jockusch, B. M., and Scherneck, S. (2000) Suppression of tumorigenicity in breast cancer cells by the microfilament protein profilin 1. J. Exp. Med. 191, 1675– 1686. (49) Zoidakis, J., Makridakis, M., Zerefos, P. G., Bitsika, V., Esteban, S., Frantzi, M., Stravodimos, K., Anagnou, N. P., Roubelakis, M. G., Sanchez-Carbayo, M., and Vlahou, A. (2012) Profilin 1 is a potential biomarker for bladder cancer aggressiveness. Mol. Cell. Proteomics MCP 11, M111.009449. (50) Wu, C.-H., Fallini, C., Ticozzi, N., Keagle, P. J., Sapp, P. C., Piotrowska, K., Lowe, P., Koppers, M., McKenna-Yasek, D., Baron, D. M., Kost, J. E., Gonzalez-Perez, P., Fox, A. D., Adams, J., Taroni, F., Tiloca, C., Leclerc, A. L., Chafe, S. C., Mangroo, D., Moore, M. J., Zitzewitz, J. A., Xu, Z.-S., van den Berg, L. H., Glass, J. D., Siciliano, G., Cirulli, E. T., Goldstein, D. B., Salachas, F., Meininger, V., Rossoll, W., Ratti, A., Gellera, C., Bosco, D. A., Bassell, G. J., Silani, V., Drory, V. E., Brown, R. H., and Landers, J. E. (2012) Mutations in the profilin 1 gene cause familial amyotrophic lateral sclerosis. Nature 488, 499–503. (51) Tanaka, Y., Nonaka, T., Suzuki, G., Kametani, F., and Hasegawa, M. (2016) Gain-of-function profilin 1 mutations linked to familial amyotrophic lateral sclerosis cause seed-dependent intracellular TDP-43 aggregation. Hum. Mol. Genet. 25, 1420–1433. (52) Arai, T., Hasegawa, M., Akiyama, H., Ikeda, K., Nonaka, T., Mori, H., Mann, D., Tsuchiya, K., Yoshida, M., Hashizume, Y., and Oda, T. (2006) TDP-43 is a component of ubiquitin-positive taunegative inclusions in frontotemporal lobar degeneration and amyotrophic lateral sclerosis. Biochem. Biophys. Res. Commun. 351, 602–611. (53) Babich, M., Foti, L. R., Sykaluk, L. L., and Clark, C. R. (1996) Profilin forms tetramers that bind to G-actin. Biochem. Biophys. Res. Commun. 218, 125–131. (54) Wopfner, N., Willeroidee, M., Hebenstreit, D., van Ree, R., Aalbers, M., Briza, P., Thalhamer, J., Ebner, C., Richter, K., and Ferreira, F. (2002) Molecular and immunological characterization of profilin from mugwort pollen. Biol. Chem. 383, 1779–1789. (55) Korupolu, R. V., Achary, M. S., Aneesa, F., Sathish, K., Wasia, R., Sairam, M., Nagarajaram, H. A., and Singh, S. S. (2009) Profilin oligomerization and its effect on poly (L-proline) binding and phosphorylation. Int. J. Biol. Macromol. 45, 265–273. ACS Paragon Plus Environment

29

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 30 of 38

(56) Nodelman, I. M., Bowman, G. D., Lindberg, U., and Schutt, C. E. (1999) X-ray structure determination of human profilin II: A comparative structural analysis of human profilins. J. Mol. Biol. 294, 1271–1285. (57) Mulder, F. A., Skrynnikov, N. R., Hon, B., Dahlquist, F. W., and Kay, L. E. (2001) Measurement of slow (micros-ms) time scale dynamics in protein side chains by (15)N relaxation dispersion NMR spectroscopy: application to Asn and Gln residues in a cavity mutant of T4 lysozyme. J. Am. Chem. Soc. 123, 967–975. (58) Vallurupalli, P., Hansen, D. F., Stollar, E., Meirovitch, E., and Kay, L. E. (2007) Measurement of bond vector orientations in invisible excited states of proteins. Proc. Natl. Acad. Sci. U. S. A. 104, 18473–18477. (59) Baldwin, A. J., Hansen, D. F., Vallurupalli, P., and Kay, L. E. (2009) Measurement of methyl axis orientations in invisible, excited states of proteins by relaxation dispersion NMR spectroscopy. J. Am. Chem. Soc. 131, 11939–11948. (60) Igumenova, T. I., Brath, U., Akke, M., and Palmer, A. G. (2007) Characterization of chemical exchange using residual dipolar coupling. J. Am. Chem. Soc. 129, 13396–13397. (61) Vallurupalli, P., Hansen, D. F., and Kay, L. E. (2008) Probing structure in invisible protein states with anisotropic NMR chemical shifts. J. Am. Chem. Soc. 130, 2734–2735. (62) Hansen, D. F., Vallurupalli, P., and Kay, L. E. (2009) Measurement of methyl group motional parameters of invisible, excited protein states by NMR spectroscopy. J. Am. Chem. Soc. 131, 12745– 12754. (63) Hansen, D. F., Vallurupalli, P., and Kay, L. E. (2008) An improved 15N relaxation dispersion experiment for the measurement of millisecond time-scale dynamics in proteins. J. Phys. Chem. B 112, 5898–5904. (64) McConnell, H. M. (1958) Reaction Rates by Nuclear Magnetic Resonance. J. Chem. Phys. 28, 430–431. (65) Press, W. H. (2007) Numerical Recipes 3rd Edition: The Art of Scientific Computing. Cambridge University Press. (66) Choy, W.-Y., Zhou, Z., Bai, Y., and Kay, L. E. (2005) An 15N NMR spin relaxation dispersion study of the folding of a pair of engineered mutants of apocytochrome b562. J. Am. Chem. Soc. 127, 5066–5072. (67) Lu, J., and Pollard, T. D. (2001) Profilin binding to poly-L-proline and actin monomers along with ability to catalyze actin nucleotide exchange is required for viability of fission yeast. Mol. Biol. Cell 12, 1161–1175. (68) Goldberg, A. L. (2003) Protein degradation and protection against misfolded or damaged proteins. Nature 426, 895–899. (69) Finkel, T., Theriot, J. A., Dise, K. R., Tomaselli, G. F., and Goldschmidt-Clermont, P. J. (1994) Dynamic actin structures stabilized by profilin. Proc. Natl. Acad. Sci. U. S. A. 91, 1510–1514. (70) Kang, F., Laine, R. O., Bubb, M. R., Southwick, F. S., and Purich, D. L. (1997) Profilin interacts with the Gly-Pro-Pro-Pro-Pro-Pro sequences of vasodilator-stimulated phosphoprotein (VASP): implications for actin-based Listeria motility. Biochemistry (Mosc.) 36, 8384–8392. (71) Buss, F., Temm-Grove, C., Henning, S., and Jockusch, B. M. (1992) Distribution of profilin in fibroblasts correlates with the presence of highly dynamic actin filaments. Cell Motil. Cytoskeleton 22, 51–61. (72) Mayboroda, O., Schlüter, K., and Jockusch, B. M. (1997) Differential colocalization of profilin with microfilaments in PtK2 cells. Cell Motil. Cytoskeleton 37, 166–177. (73) Delaglio, F., Grzesiek, S., Vuister, G. W., Zhu, G., Pfeifer, J., and Bax, A. (1995) NMRPipe: a multidimensional spectral processing system based on UNIX pipes. J. Biomol. NMR 6, 277–293. (74) Helgstrand, M., Härd, T., and Allard, P. Simulations of NMR pulse sequences during equilibrium and non-equilibrium chemical exchange. J. Biomol. NMR 18, 49–63. ACS Paragon Plus Environment

30

Page 31 of 38

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

(75) Bouvignies, G., and Kay, L. E. (2012) Measurement of proton chemical shifts in invisible states of slowly exchanging protein systems by chemical exchange saturation transfer. J. Phys. Chem. B 116, 14311–14317.

ACS Paragon Plus Environment

31

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 38

TOC Figure

ACS Paragon Plus Environment

32

Paragon Page Biochemistry 33Plus of 38Environm

Biochemistry Page 34 of 38

1 2 3 Paragon Plus Environm 4 5 6

Page Biochemistry 35 of 38

1 2 3 Paragon Plus Environm 4 5 6

S Paragon Plus36 Environm Biochemistry Page of 38

Page 37 of 38

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

Biochemistry

ACS Paragon Plus Environment

Biochemistry Page 38 of 38

Paragon Plus Environm 1 2