Matthew A.
Molski
a,
Jessica L.
Goodman
b,
Fang-Chieh
Chou
c,
David
Baker
*d,
Rhiju
Das
*c and
Alanna
Schepartz
*ae
aDepartment of Chemistry, Yale University, New Haven, CT 06511, USA. E-mail: matthew.molski@yale.edu; Fax: +1 203-432-3486; Tel: +1 203-432-8276
bWhitehead Institute for Biomedical Research, Nine Cambridge Center, Cambridge, MA 02142, USA. E-mail: jgoodman@wi.mit.edu; Fax: +1 617-258-7226; Tel: +1 617-258-5184
cDepartment of Biochemistry, Stanford University, Stanford, CA 94305, USA. E-mail: fcchou@stanford.edu; rhiju@stanford.edu; Fax: +1 650-723-5976; Tel: +1 650-723-7310
dDepartment of Biochemistry, University of Washington, Seattle, WA 98195, USA. E-mail: dabaker@u.washington.edu; Fax: +1 206-685-1792; Tel: +1 206-543-1295
eDepartments of Chemistry and of Molecular, Cellular and Developmental Biology, Yale University, New Haven, CT 06511, USA. E-mail: alanna.schepartz@yale.edu; Fax: +1 203-432-3486; Tel: +1 203-432-5094
First published on 21st September 2012
Natural biopolymers fold with fidelity, burying diverse side chains into well-packed cores and protecting their backbones from solvent. Certain β-peptide oligomers assemble into bundles of defined octameric stoichiometry that resemble natural proteins in many respects. These β-peptide bundles are thermostable, fold cooperatively, exchange interior amide N–H protons slowly, exclude hydrophobic dyes, and can be characterized at high resolution using X-ray crystallography – just like many proteins found in nature. But unlike natural proteins, all octameric β-peptide bundles contain a sequence-uniform hydrophobic core composed of 32 leucine side chains. Here we apply rational design principles, including the Rosetta computational design methodology, to introduce sequence diversity into the bundle core while retaining the characteristic β-peptide bundle fold. Using circular dichroism spectroscopy and analytical ultracentrifugation, we confirmed the prediction that an octameric bundle still assembles upon a major remodelling of its core: the mutation of sixteen core β-homo-leucine side chains into sixteen β-homo-phenylalanine side chains. Nevertheless, the bundle containing a partially β-homo-phenylalanine core poorly protects interior amide protons from exchange, suggesting molten-globule-like properties. We further improve stability by the incorporation of eight β-homo-pentafluorophenyalanine side chains, giving an assembly with amide protection factors comparable to prior well-structured bundles. By demonstrating that their cores tolerate significant sequence variation, the β-peptide bundles reported here represent a starting point for the “bottom-up” construction of β-peptide assemblies possessing both structure and sophisticated function.
In this work, we applied the Rosetta software package29,30 to predict β-peptide sequences that could effectively recapitulate the structurally characterized β-peptide bundle core using a mixture of leucine and non-leucine side chains. One such sequence (Acid-1YFF), containing an equal number of core β-homo-phenylalanine and β-homo-leucine residues, assembles into a 314-helical, relatively thermostable, octameric bundle. Despite this stability, Acid-1YFF displayed two properties associated with a molten globule state: rapid amide NH exchange and conformational heterogeneity as judged by NMR, indicating the potential for improved side chain organization within the core. It is well known that mixtures of phenylalanine and penta-fluorophenylalanine can improve protein stability when introduced into an otherwise all-hydrocarbon protein core.31–40 To test whether the thermodynamic stability and structural uniqueness of Acid-1YFF could be improved by fluorocarbon substitution, we synthesized an analogue (Acid-1YFF★) of Acid-1YFF containing β-pentafluoro-homo-phenylalanine (F5βPhe) at position 8. Acid-1YFF★ displays improved folding properties, resulting in a more conformationally distinct and stable core as judged by NMR. Acid-1YFF and Acid-1YFF★ are the first β-peptide bundles containing mixed sequence hydrophobic cores, suggesting further use of Rosetta and rational design principles as tools to remodel novel β-peptide bundles.
Fig. 1 β-peptide bundles studied in this work. (A) Helical net diagram and sequences of Base-1F, Acid-1F, Zwit-1F, Acid-1Y, Acid-1YFF, and Acid-1YFF★. F represents β-homo-phenylalanine; F★ represents β-homo-pentafluorophenylalanine. Colors distinguish side chains on the aromatic (β-hY-containing), hydrophobic (β-hL containing) and salt bridge (β-hO- and β-hD-containing) faces. Also shown are ln Ka values characterizing each octameric assembly, as determined previously or in this work by SE-AU. (B) Ribbon representation of the crystal structure of Acid-1Y highlighting the packing of the leucine side chains, along with computationally predicted structure of Acid-1YFF and a color-coded guide to the locations of F5β-hPhe (green) in Acid-1YFF★. |
Next we searched for alternative side chains that could be accommodated in place of some or all of the β-homo-leucines within the bundle interior. While the fixed-backbone conformational search described above gave the all-leucine core as the optimal solution, a small void observed in the bundle center suggested that larger side chains might be tolerated at positions 5, 8 and/or 11 with minor backbone adjustments. We modeled the effects of introducing between four and eight alternative, proteinogenic side chains at these positions within Acid-1Y, and optimized the backbone torsion angles with a constraint potential tethering the angles to within ∼10° of their crystallographic values. Rosetta calculations predicted variant bundles would generally be poorer in energy than the starting bundle. However, a few were predicted to possess better hydrophobic packing of non-polar side-chains, as assessed by the sum of van der Waals and solvation energies. In particular, a variant containing β-homo-phenylalanine residues at positions 5 and 8 (Acid-1YFF, Fig. 1A) gave more favorable hydrophobic packing energy by ∼1 kT/monomer. Acid-1YFF contains 16 β-homo-leucine to β-homo-phenylalanine substitutions, placing an additional 48 carbon atoms into the bundle core. Nevertheless, the backbone atoms of the Acid-1YFF★ bundle were shifted by less than 0.8 Å RMSD from those in the starting bundle (Fig. 1B), suggesting that the bulky aromatic side chains in the interior could be accommodated without disrupting the octamer.
The β-peptide monomer Acid-1YFF was prepared using solid phase, microwave-assisted methods, purified to homogeneity by HPLC, and characterized initially using wavelength-dependent circular dichroism (CD) spectroscopy. As predicted, Acid-1YFF underwent a concentration-dependent increase in 314-helical structure (as judged by the molar residue ellipticity at 209 nm, MRE209)48 between 12 and 200 μM (Fig. 2a), consistent with an equilibrium between a partially structured monomer and a folded oligomer. A plot of MRE209vs. [Acid-1YFF] was first fit to a monomer–octamer equilibrium with ln Ka = 66.9 ± 0.5, suggesting that oligomerization of Acid-1YFF was less favorable than that of the Acid-1Y bundle (ln Ka = 82.5 ± 1.8).4 The fit, however, was imperfect [P = 2 × 10−10; see Table S3 and Fig. S1†], with a closer agreement at higher concentrations ([Acid-1YFF] > 50 μM) than at lower concentrations ([Acid-1YFF] < 50 μM). A plot of MRE209vs. [Acid-1YFF] was subsequently fit to alternative, three-state models containing either a dimeric (1–2–8) or tetrameric (1–4–8) intermediate (Fig. 2a). The resulting association constants for the 1–2–8 model were ln Ka1 = 18.4 ± 4.1 and ln Ka2 = 79.2 ± 9.8, and those for the 1–4–8 model were ln Ka1 = 36.8 ± 3.2 and ln Ka2 = 73.7 ± 11.2. Both three-state models fit the CD data substantially better than the original two-state monomer–octamer model (1–8), albeit with greater uncertainties in equilibrium constants due to the added fit parameters. Both fits gave predictions within the error of the data (Table S3 and Fig. S1†). Temperature-dependent CD studies also supported formation of a relatively stable Acid-1YFF bundle; the TM of a 200 μM solution of Acid-1YFF was 52 °C (Fig. 2b); this value is lower than that of Acid-1Y, whose TM was 82 °C at 150 μM.4 Taken together, these CD data established the relative stability of the Acid-1YFF bundle, but could not precisely define the stoichiometry of the putative intermediate or the difference in stability between Acid-1Y and Acid-1YFF.
Fig. 2 Self-association of designed β-peptide bundles. Circular dichroism spectra of Acid-1YFF (a and b) and Acid-1YFF★ (c and d) as a function of concentration (a and c) and temperature (b and d). Plots of MRE209 as a function of [β-peptide] were fit to a monomer–dimer–octamer (1–2–8) equilibrium (dotted line), a monomer–tetramer–octamer (1–4–8) equilibrium (dashed line) or monomer–octamer equilibrium (solid line). Inset: wavelength-dependent CD spectra of Acid-1YFF and Acid-1YFF★ (MRE in units of 103 deg cm2 dmol−1). The TM is defined as the maximum of a plot of δMRE209·δT−1versus temperature. |
We therefore turned to sedimentation equilibrium analytical ultracentrifugation to more precisely characterize the stoichiometry and relative stability of the Acid-1YFF bundle. Sedimentation of Acid-1YFF at concentrations of 20, 80 and 200 μM was monitored at four speeds (36000, 42000, 50000 and 60000 RPM). The AU data was fit to both two-state monomer-n-mer equilibrium models as well as three-state models proceeding through a dimer or tetramer intermediate (1–2–8 or 1–4–8, respectively). Poor fits with high RMSD values and larger and more systematic residuals were observed when n in the two-state model was set to any value other than 8 between 2 and 10 (see Fig. S2†). Both the 1–2–8 and the 1–4–8 model fit the AU data better than the two state 1–8 model (P < 10−6; F-test, see ESI†), in agreement with the CD analysis, and the 1–4–8 model produced the best fit. The ln Ka values calculated from the monomer–tetramer–octamer fit were ln Ka1 = 27.6 ± 0.2 and ln Ka2 = 68.2 ± 0.3, and the position of the monomer–octamer equilibrium (ln Ka2, 68.2 ± 0.3) agrees within error with the ln Ka2 value determined by CD (73.7 ± 11.2). Taken together, the CD and AU data suggest that Acid-1YFF assembly is a three state process involving a tetrameric intermediate towards a final octameric assembly, and that the oligomerization of Acid-1YFF is less favorable than that of Acid-1Y (ln Ka = 82.5 ± 1.8).4
Many de novo designed proteins exist as molten globules, and thus we hypothesized that the lower thermodynamic stability of the Acid-1YFF bundle compared to Acid-1Y might signal the presence of an undefined or heterogeneous hydrophobic core.49 Molten globules often bind and increase the fluorescence of dyes such as 1-anilino-8-naphthalenesulfonate (ANS) by factors ΔF > 100.49 By contrast, well-folded or unfolded proteins do not provide favorable ANS binding sites, and elicit little or no change in ANS fluorescence (ΔF < 10). All octameric β-peptide bundles reported thus far, including the Acid-1Y bundle, behave like well-folded proteins, causing minimal (<2-fold) changes in ANS fluorescence. Like previously characterized β-peptide bundles, even at concentrations as high as 200 μM (70% octamer), Acid-1YFF had little or no effect on the fluorescence of ANS (ΔF = 2) (Fig. S4†). This analysis suggests that the Acid-1YFF β-peptide bundle provides only limited access of solvent to its hydrophobic core.
Greater insight into the differences between the Acid-1YFF bundle relative to Acid-1Y was revealed by NMR experiments that monitored the rate of amide NH hydrogen/deuterium exchange. Unlike all previously characterized β-peptide bundles,4–6,9 the Acid-1YFF spectrum revealed no slowly exchanging amide NH protons; no amide NH resonances were visible as soon as 15 min after addition of D2O at 25 °C (Fig. 3a). In addition, relative to the NMR spectrum of Acid-1Y,4 the NMR spectrum of Acid-1YFF showed significant line broadening in the adjacent aromatic region. These data suggest that the aromatic sub-core of the Acid-1YFF bundle, although octameric and relatively thermostable, possesses significant conformational heterogeneity on the NMR time scale.
Fig. 3 Selective substitution of β-homo-pentafluorophenylalanine improves β-peptide bundle stability as judged by NMR. (a) The amide NH regions of the Acid-1YFF and (b) Acid-1YFF★ β-peptide bundles after dissolution in D2O. In the case of Acid-1YFF, no amide NH resonances are visible within 15 min after addition of D2O. In the case of Acid-1YFF★, the amide signals persist for more than 1 h. (c) Integration of the resonances indicated in (B) normalized to the aromatic reference peak and fit to exponential decays to determine average exchange rate constants. |
Fluorocarbon side chains are more hydrophobic than their hydrocarbon counterparts50 and mixtures of phenylalanine and penta-fluorophenylalanine can improve protein stability when introduced into an otherwise all-hydrocarbon protein core.31–40 In certain cases, favorable face-to-face interactions between phenylalanine and pentafluorophenylalanine side chains account for improved stability,51–54 while in others steric and/or hydrophobic effects have been invoked.33,55 To evaluate whether the thermodynamic stability and structural uniqueness of Acid-1YFF could be improved by fluorocarbon substitution, we synthesized an analogue of Acid-1YFF containing β-pentafluoro-homo-phenylalanine (F5βPhe) at position 8 (Acid-1YFF★). Previous work has shown that β-peptide bundles containing β-homo-hexafluoroleucine at position 8 contain a discrete fluorous sub-domain and retain the characteristic β-peptide bundle fold.8
The Acid-1YFF★ sequence was synthesized and its assembly characterized by CD, SE-AU, ANS binding, and NMR to evaluate the effects of β-homo-pentafluorophenylalanine substitution on β-peptide bundle structure and stability. As expected, the CD spectrum of Acid-1YFF★ was characterized by concentration-dependent changes in MRE212, and the data was first fit to a monomer–octamer equilibrium with a ln Ka of 73.0 ± 0.5. Again the fits deviated from the experimental data in the lower concentration range ([Acid-1YFF★] < 25 μM (Fig. S1†) and could be improved by including a dimer or tetramer intermediate along the folding pathway. Fitting to a monomer–dimer–octamer equilibrium (1–2–8) resulted in ln Ka values of 18.4 ± 4.0 (ln Ka1) and 82.9 ± 2.3 (ln Ka2) (Fig. 2c), and fitting to a monomer–tetramer–octamer equilibrium (1–4–8) yielded ln Ka values of 41.4 ± 4.6 (ln Ka1) and 84.7 ± 3.0 (ln Ka2). Both fits agreed with the CD data within error (Fig. S1 and Table S3†). Although measured with large errors, the ln Ka2 values suggested that the Acid-1YFF★ octamer stability was greater than the initial Acid-1YFF design and now comparable to the starting Acid-1Y bundle with the all-leucine core (ln Ka of 82.5 ± 1.8). The temperature-dependence of the Acid-1YFF★ CD signal was consistent with improved stability; the TM of Acid-1YFF★ at 200 μM was 67 °C (Fig. 2d), a significantly higher value than the TM of 52 °C for Acid-1YFF at the same concentration.
We turned again to SE-AU to more precisely and accurately determine the relative association constants of the Acid-1YFF and Acid-1YFF★ bundles. The equilibrium sedimentation of Acid-1YFF★, performed at 20, 80 and 200 μM, was fit to two-state and three-state equilibrium models as described previously. Once again, the data fit best to a three-state model with a tetramer intermediate (1–4–8 model), although the 1–2–8 model was also an improvement over a two-state monomer–octamer equilibrium (P-value < 10−6 for both three state models according to the F-test; see ESI†). The ln Ka values calculated from the 1–2–8 fit were ln Ka1 = 11.6 ± 0.4 and ln Ka2 = 79.7 ± 1.3, and those from the 1–4–8 fit were ln Ka1 = 29.7 ± 0.2 and ln Ka2 = 72.6 ± 0.2 (Fig. S3†). Comparison to the AU analysis of the initial design indicates that, irrespective of model, the Acid-1YFF★ bundle is significantly more stable than the Acid-1YFF bundle.
The effects of perfluoro substitution on the aromatic β-peptide bundle core were studied further by monitoring the effect of Acid-1YFF★ on the intrinsic fluorescence of ANS (see ESI†). In previous work, designed α-amino acid proteins containing fluorous cores induced larger increases in ANS fluorescence than analogous proteins containing hydrocarbon cores; for example, the perfluorinated helical bundle α4-F2 increased ANS fluorescence by 5-fold.56 This value is greater than that observed in the presence of α4 (ΔF = 3), but less than that seen in the presence of classic molten globule states such as lactalbumin (ΔF ≥ 135).57,58 Whereas the Acid-1YFF bundle, like the Acid-1Y bundle, led to little or no increase in ANS fluorescence intensity (ΔF = 2-fold at [Acid-1YFF] = 200 μM), at an equivalent concentration the Acid-1YFF★ bundle caused ANS fluorescence to increase >19-fold (Fig. S4†). Interestingly, the ΔF of ANS observed in the presence of Acid-1YFF★ is comparable to that observed in the presence of the Zwit-8L*, in which eight leucine side chains were substituted by hexafluoroleucine (ΔF = 15), and whose structure is known at atomic detail.8 Thus, it may be that the increased fluorescence of ANS in the presence of the Acid-1YFF★ and Zwit-8L* results from an intrinsic affinity of ANS for fluorous domains, and not necessarily from a molten core.
Finally, H/D exchange NMR experiments gave the most incisive assessment of whether introduction of β-homo-pentafluorophenylalanine into the Acid-1YFF β-peptide bundle core would impart greater stability and structural uniqueness. Whereas the amide NH resonances in the Acid-1YFF bundle exchanged completely within 15 min, those of Acid-1YFF★ were better resolved and exchanged over the course of hours (Fig. 3b and c). The average rate of H/D exchange for two regions of the Acid-1YFF★ spectrum each fit a first-order decay function with a rate constant (kex) of 2.2 × 10−4 s−1. By comparison, the exchange rate constant for the random coil model poly-β-alanine (β3-homo-glycine) (krc) is 4.4 s−1. The ratio of these values (krc/kex), often defined as a protection factor P, is 2.0 × 104, a value comparable to that of previously reported β-peptide bundles (P between 9 × 103 and 6 × 104),5,6 including Acid-1Y (P = 6 × 104).4 This equivalence of the P values calculated for Acid-1YFF★ underscores the effect that the F5βPhe substitution has on the stability of the diversified hydrophobic core. The octamer bundle containing a mixed core of β-homo-leucine, β-homo-phenylalanine, and β-homo-pentafluorophenylalanine is characterized by monomer–octamer equilibrium constants and NMR amide hydrogen exchange rates that are comparable to previous bundles containing homogeneous cores.
Footnote |
† Electronic supplementary information (ESI) available: β-Peptide design and synthesis and characterization using CD, SE-AU, and NMR. See DOI: 10.1039/c2sc21117c |
This journal is © The Royal Society of Chemistry 2013 |