Full Text

Turn on search term navigation

ARTICLE

Received 26 Aug 2014 | Accepted 17 Mar 2015 | Published 22 Apr 2015

Switchable proteins that can be regulated through exogenous or endogenous inputs have a broad range of biotechnological and biomedical applications. Here we describe the design of switchable enzymes based on an ensemble allosteric model. First, we insert an enzyme domain into an effector-binding domain such that both domains remain functionally intact. Second, we induce the fusion to behave as a switch through the introduction of conditional conformational exibility designed to increase the conformational entropy of the enzyme domain in a temperature- or pH-dependent fashion. We conrm the switching behaviour in vitro and in vivo. Structural and thermodynamic studies support the hypothesis that switching result from an increase in conformational entropy of the enzyme domain in the absence of effector. These results support the ensemble model of allostery and embody a strategy for the design of protein switches.

DOI: 10.1038/ncomms7968

Design of protein switches based on an ensemble model of allostery

Jay H. Choi1, Abigail H. Laurent1, Vincent J. Hilser2 & Marc Ostermeier1

1 Department of Chemical and Biomolecular Engineering, Johns Hopkins University, 3400N Charles Street, Baltimore, Maryland 21218, USA. 2 Department of Biology, Johns Hopkins University, 3400N Charles Street, Baltimore, Maryland 21218, USA. Correspondence and requests for materials should be addressed to M.O. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 6:6968 | DOI: 10.1038/ncomms7968 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms7968

The ability to design allosteric proteins or protein switches has many applications from basic science, to biosensors, to selective protein therapeutics1. While a handful of switch

design strategies have met with some degree of success2,3, the design of switches for an application prescribed input signal (for example, effector and environmental conditions) and desired output function (for example, binding, catalysis, uorescence and so on) remains a challenge. Previously, we pursued a directed evolution approach to switch design in which the switches are built from existing protein domains with the prerequisite effector-binding (input) and functional (output) domains46. We constructed libraries that explore the sequence space comprising insertions of one domain into the other, with or without circular permutation of the insert domain. We identied switches from these libraries via suitable screens or genetic selections. These studies taught us that switches are rare among the domain fusions that retain the output functionality. A key question is what do the switches possess that non-switch fusions lack? While biochemical and biophysical studies of our switches have provided some insights into this question5,710, we decided to pursue a complementary approach by designing switches based on an ensemble allosteric model (EAM), wherein sensitivity can be tuned by remodelling the energy landscape11.

The allosteric regulation of protein function is central to essentially all aspects of cellular function. The phenomenological and structural models that have dominated the view of allostery have provided many insights1215. However, recent studies have illuminated the importance of protein dynamics in allostery11,1625 and conrmed the predictions of Cooper and Dryden26 that allostery can manifest in the absence of signicant changes in the ground state structure of the protein17,22,27. The importance of conformational degeneracy in allostery is encompassed in the conformational EAM2830. This model possesses the same characteristic of any ensemble model in that it posits that all possible conformations (or states) of a protein are populated according to their respective energies, and allosteric effectors remodel this landscape. The EAM2830, however, also provides a vehicle for rationalizing the relationship between binding and conformational changes in each domain, and thus provides insight into the role of conformational degeneracy

(that is, conformational entropy) in determining the overall allosteric effect. Allostery manifests, for example, from having the effector shift the ensemble from a set of states for which the ensemble-weighted average is less active to a set of states for which the ensemble-weighted average is more active by changing the relative stabilities of the different states. One way to achieve this is to have the ensemble in the absence of effector populating many states (especially less active ones) and have binding of the effector preferentially stabilize active states. This role for conformational entropy in allostery is central to the design of allosteric proteins described herein.

In nature, protein exibility has shown to be important for protein function31. In some cases, cold-adapted enzymes incorporate more glycine residues, which modulate functional activity by altering the exibility of the native state or increasing the population of states that do not bind the substrate in the absence of that substrate3234. Similarly, the substitution of glycines at sites distant from the active site of adenylate kinase affects its activity by selectively increasing the probability of locally unfolded states (and thus altering conformational entropy) without perturbing the structural properties of the native fold19. Based on these studies, we hypothesized that the appropriate introduction of conformational exibility into non-switch fusion proteins could impose differential effects on the structural uctuations on the output domain in the presence and absence of effector, and thereby confer switch properties to the fusion.

Here, we design fusion proteins to possess switching behaviours via the introduction of linkers designed to increase the conformational entropy of the enzyme domain in a temperature- or pH-dependent fashion. We conrm the switching behavior with cellular and enzymatic assays, and we present structural and thermodynamic studies that support our hypothesis. The switches we create have multilevel control, requiring both a particular environmental condition (temperature or pH) and ligand for activation.

ResultsDesign of switches. We have previously shown how protein switches can be constructed by fusing a ligand-binding domain

Inactive effector

binding domain

Active effecotor

binding domain

Energy

Inactive enzyme

domain

Active enzyme

domain

Fusion

protein

Representative non-allosteric ensembles Representative allosteric ensembles

Effector Effector

Mutations

that

increase

entropy

Energy

Figure 1 | Tuning an allosteric ensemble through mutations that promote conformational entropy. (a) The fusion proteins comprise two domains connected by two linkers (since one domain is inserted into the other). Each domain can be active or inactive (together with or separately from the other domain) resulting in four categories of representative states as shown. (b) Non-allosteric fusions lack allostery because states with an inactive enzyme domain are much less stable than states with an active enzyme domain regardless of the presence of ligand. (c) Allostery can be introduced though the introduction of conformational exibility through mutation, where the exibility increases the probability of the inactive states. Thus, the desired mutations remodel the energy landscape to increase the conformational entropy of the enzyme domain in the absence of ligand. The effector remodels the landscape by increasing the stability of states that bind the effector and thus their population in the ensemble, resulting in an effector-dependent increase in the observed enzyme activity. Although the schematic depicts the creation of positive allostery, negative allostery can be established in the same manner, as can be illustrated by switching the active and inactive domains in this gure (see Hilser and Thompson29).

2 NATURE COMMUNICATIONS | 6:6968 | DOI: 10.1038/ncomms7968 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms7968 ARTICLE

and an enzyme such that the ligand acts as an effector regulating enzyme activity4,6. However, fusions that will exhibit switching are difcult to predict. Fusions that retain at least some enzymatic activity rarely possess effector regulation. We sought to develop a strategy to convert such non-switch fusion proteins into switches via control of local conformational entropy using engineered linkers possessing conditional conformational exibility.

We cast the problem in terms of an EAM11,28,29 (Fig. 1). We hypothesized that most enzymatically active fusions lack allostery because the enzymatically active states of the fusion are considerably more stable than the inactive states even in the absence of the effector (Fig. 1b). Thus, enzymatically active states of the fusion dominate the distribution of states with or without effector. For such a scenario, allostery could be established by remodelling the energy landscape such that the active states in the absence of ligand have similar or lower stabilities than the inactive states in the absence of ligand (Fig. 1c). We hypothesized that such an increase in conformational entropy in the enzyme domain could be achieved by the proper introduction of conformational exibility in the protein. Since conformational exibility can be externally controlled (for example, through temperature), the resulting switches would be expected to have multilevel control, requiring the induction of conformational exibility for the switching property.

We chose to test this hypothesis for creating allostery using c4, a fusion protein between maltose-binding protein (MBP) and TEM1 b-lactamase (BLA), in which BLA is inserted into MBP in place of residues 317 and 318 of MBP, and is connected to residue 319 of MBP via a linker with the sequence DKT (ASP-LYS-THR)9

(Fig. 2a and Supplementary Table 1). c4 is not an allosteric protein and the c4 gene does not confer to Escherichia coli a maltose-dependent resistance to ampicillin (Amp)10. We reasoned that the introduction of conformational exibility in the linker region of fusion proteins was most likely to cause the emergence of allostery in this fusion protein since the linker lies at the interface between the two domains (Fig. 2b). Furthermore, the introduction of changes in the linker minimized the chance of directly disrupting enzymatic activity or effector afnity.

We chose two methods to introduce exibility: the substitution of glycines in the linker and the substitution of a peptide linker that undergoes pH-dependent conformational changes between a-helix and a random coil. Both of these methods afforded us non-mutagenic methods to control conformational entropy. Thus, we could interrogate the relationship between linker conformational exibility and switching by changing the temperature for the glycine linker constructs and by changing pH for the pH-responsive linker constructs. This ability to control conformational entropy would result in multi-input switches, requiring the correct temperature or pH to be turned on by the effector.

pH-induced conformational conversion of (EAQA)n peptides. Whereas the introduction of glycines was straightforward, the introduction of a peptide that would undergo pH-dependent conversion from an a-helix to a random coil required further design considerations. In previous studies, an amphipathic peptide was designed to interact with uncharged bilayers in a pH-dependent fashion through pH-induced random coil-a-helical transition35,36. We applied these design principles and chose an EAQA repeat peptide as a putative pH-responsive peptide (Fig. 3a,b). We designed the EAQA peptide to form an a-helix under acidic conditions. We positioned glutamate residues to align on one side of the helix so that they would destabilize the helix through charge repulsion at pH 7.5.

We evaluated our design using structural modelling, molecular dynamic (MD) simulations and spectroscopic characterization. We conducted a 100-ns MD simulation using an a-helix model of the (EAQA)3 peptide as a starting structure under three pH conditions: 3, 5 and 7 (Fig. 3c,d). Analysis of the positional root mean squared deviations of the backbone relative to the starting structure indicated that while the helix conformation at pH 3 remained relatively stable, the conformational conversion to random coil progressed as pH increased.

We synthesized a 30-residue EAQA peptide, WA-(EAQA)7, that included a tryptophan at the N terminus and monitored the variation of its helical content by scanning the circular dichroism (CD) and tryptophan uorescence spectra in solution at pH4.57.5 (Fig. 3e,f). Ellipticity increased and uorescence intensity decreased with increasing pH. Both ndings are indicative of the desired a-helix to random coil structural change. We concluded that (EAQA)n peptides undergo a pH-dependent conversion from an a-helix to a random coil as designed.

Designed proteins confer switch phenotype. We constructed a series of variants of c4 in which one to three glycines were substituted into the DKT linker (Fig. 2a). In two constructs, glycines were also substituted at the other fusion point between the BLA and MBP domains. In addition, we created variants of c4 in which two to four repeats of the EAQA peptide replaced the DKT linker.

We assayed for switching activity by measuring the minimum inhibitory concentration (MIC) for Amp for cells expressing these proteins in the presence and absence of maltose. We have previously shown that genes encoding MBP-BLA switches confer

c4:

c4-1G:

MBP 1-316 GGT

BLA 24-286 MBP 319-370

c4-2G:

c4-3G:

MBP 1-316 GKT

BLA 24-286 MBP 319-370

MBP 1-316 GGG

BLA 24-286 MBP 319-370

c4-4G:

MBP 1-314 GGG

BLA 24-286 MBP 319-370

c4-5G: GG

MBP 1-315 GGG

BLA 24-286 MBP 319-370

c4-2hx:

MBP 1-316 DKT

BLA 24-286 MBP 319-370

c4-3hx:

MBP 1-316 (EAQA)2

BLA 24-286 MBP 319-370

c4-4hx:

MBP 1-316 (EAQA)3

BLA 24-286 MBP 319-370

MBP 1-316 (EAQA)4

BLA 24-286 MBP 319-370

MBP 1-316

MBP 319-370

Linker: DKT

MBP 316

BLA 24

BLA 24-286

Figure 2 | Sequences and structural models of c4 and c4 variants.(a) Graphical illustration of the sequences of c4 and c4 variants indicating the MBP domain (red), BLA domain (blue), engineered glycine linker (yellow) and EAQA linker (green). (b) Structural model of the c4 fusion protein comprises MBP (red) and TEM-1 BLA (blue) domains with the DKT linker and fusion sites indicated by arrows.

NATURE COMMUNICATIONS | 6:6968 | DOI: 10.1038/ncomms7968 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms7968

Gln

Glu

a b

Gln

Glu

30 25 20 15 10

5 0

3.53.02.52.01.51.00.5 0

Trp

(+ Mal / Mal)

MIC Ampratio

(+ Mal / Mal)

35 30 25 20 15 10

5 0

pH 6 pH 7 pH 7.5 pH 8

Ala

Glu

MIC Ampratio

Isoelectric point

3.67

Gln

MBP

18 C 37 C

c4 4hx

c4-3hx

2hx

3hx

Gln

c d

3.58

3.51

Ratio of initial velocity

(+ Mal / Mal)

18 C 37 C

pH 6 pH 8

Glu Glu

Ratio of initial velocity

(+ Mal / Mal)

3.53.02.52.01.51.00.5 0

Glu

Trp

Glu

c d

86420 100

c4-3G

r.m.s.d. ()

pH 7

pH 5

pH 3

10 50 100

Time (ns)

Figure 4 | Introduction of conformational exibility results in cells with a switching phenotype and allosteric proteins. (a) Maltose-dependent fold increase in Amp resistance for cells expressing the indicated proteins at 18 and 37 C. (b) Maltose-dependent fold increase in Amp resistance for cells expressing the indicated proteins grown in media with the indicated initial pH. The MICAmp ratio ( maltose/ maltose) was calculated

using the median MICAmp value from three independent trials. (c) Maltose-dependent fold increase in initial rate of nitrocen hydrolysis enzymatic activity of c4 and c4-3G at 18 and 37 C. (d) Maltose-dependent fold increase in initial rates of nitrocen hydrolysis enzymatic activity of c4 and c4-3hx at pH 6 and 8, respectively. Average rates were calculated from three independent trials.

Time (ns)

e f

([afii9835])

(deg cm2 per dmol)

500 0.80.60.40.2 0

500

Fluorescence

intensity (AU)

1,000

200

220

240

350

325

300

400

375

425

Wavelength (nm) Wavelength (nm)

Figure 3 | Conformational conversion of the EAQA peptide. (a) a-Helical model of the WA-(EAQA)n peptide structure, (b) with isoelectric points and dimensions for two- to four-coiled a-helical structures: (EAQA)24 and a

cross-sectional view, showing side-chain residues. (c) MD simulation of the a-helical model of (EAQA)3 peptide with increasing pH: pH 3 (black),pH 5 (blue) and pH 7 (red). Backbone r.m.s.d. relative to the initial structure as a function of simulation time at 300 K and (d) nal structures after 10, 50 and 100-ns simulation. (e) CD spectra and (f) uorescence emission spectra of WA-(EAQA)7 at pH 4.5 (solid line), pH 5.5 (dashed line),pH 6.5 (dash-dotted line) and pH 7.5 (dotted lines).

to E. coli cells a switching phenotype, that is, maltose-dependent resistance to Amp5,9,10. This phenotype arises from one or both of the following two mechanisms: maltose binding increases the specic Amp hydrolysis activity of the switch or maltose binding increases the cellular abundance of the switch, resulting in an increase in Amp hydrolysis activity in the cell9,10. Although the latter mechanism may not be classically thought of as allostery, it can be encompassed by a broader denition, as maltose modulates the properties of the protein in the cell.

The switching phenotype is quantied by the fold increase in MIC afforded by the presence of maltose. For the glycine variants, we tested for a switching phenotype at both 18 and 37 C. According to our hypothesis, switching should preferentially emerge at 37 C. The progressive introduction of glycine residues (up to four glycines) monotonically increased the MICAmp ratio to 32-fold at 37 C, whereas no switching phenotype was observed at 18 C (Fig. 4a and Supplementary Tables 2 and 3). For the EAQA repeat variants, we tested for a switching phenotype at 37 C on media at different initial pHs ranging from 6.0 to 8.0. The fusion proteins will experience these pH differences because they are expressed in the periplasm of E. coli. According to our hypothesis, switching should preferentially emerge at higher pH.

We observed that the switching phenotype monotonically increased with pH and with the number of EAQA repeats, plateauing at three repeats and pH 8.0 with a 16-fold maltose-dependent switching phenotype (Fig. 4b and Supplementary Table 4). In both sets of proteins, switching emerged as predicted by the model of Fig. 1 through a reduction in MIC in the absence of maltose.

Designed proteins are multi-input allosteric enzymes. We puried c4 and c4-3G, and measured the rate of enzymatic hydrolysis of the BLA colorimetric substrate nitrocen in the presence and absence of maltose at 18 and 37 C (Fig. 4c and Supplementary Table 5). c4-3Gs catalytic activity depended on maltose at 37 C but not at 18 C. At 37 C, the substitution of glycines in the linker of c4 established allostery by lowering the catalytic activity in the absence of maltose by 46% but decreasing the activity in the presence of maltose by only 12%, such that maltose increased catalytic activity twofold (2.050.15). At 18 C, a temperature at which conformational exibility of the linker is reduced, the introduction of glycines in the linker of c4 resulted in similar o10% decrease in catalytic activity in the presence and absence of maltose and no switching. The induction of switching by shifting the temperature to 37 C was reversible, as shifting the temperature back to 18 C resulted in the loss of the switching. Thus, the maltose-controlled switching property of c4-3G arises in temperature-dependent fashion, and c4-3G is a multi-input, allosteric switch.

We compared the enzymatic activity of c4 and c4-3hx in the presence and absence of maltose at pH 6 and 8, respectively (Fig. 4d and Supplementary Table 6). Only c4-3hx exhibited a maltose-dependent catalytic activity only at pH 8. At pH 6, the substitution of the pH-responsive (EAQA)3 linker into c4 resulted in similar small increases in activity in the presence and absence of maltose and no switching. However, at pH 8, the introduction of the (EAQA)3 linker reduced catalytic activity in the absence of

4 NATURE COMMUNICATIONS | 6:6968 | DOI: 10.1038/ncomms7968 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms7968 ARTICLE

maltose more than in the presence of maltose resulting in an approximately threefold (2.740.23) maltose-dependent increase in enzyme activity. The induction of switching by shifting the pH to 8 was reversible, as shifting the pH back to 6 by dialysis resulted in the loss of switching. Thus, c4-3hx is a multi-input, allosteric switch that requires the correct pH to exhibit switching.

BLA domain experiences increased local structural entropy. We next sought evidence that the switching properties of c4-3G and c4-3hx arose because the linker caused a change in the conformational entropy of the BLA domain in a temperature-and pH-dependent fashion, respectively. We constructed structural models of c4, c4-3G and c4-3hx, the latter both with the (EAQA)3 peptide in helix and random coil form (Supplementary

Fig. 1). We analysed the effect of the new linkers on the local conformational entropy and stability using the ensemble-based structural thermodynamics algorithm COREX/BEST using these models. The algorithm predicted that substitution of the three glycines into c4 causes a signicant decrease in stability near the fusion sites and an overall decrease in stability in the BLA domain (Fig. 5a). Conversion of the (EAQA)3 peptide from helix form to random coil was predicted to reduced stability as well, but only in the (EAQA)3 linker (Fig. 5b).

We next acquired [15N and 1H] TROSY-HSQC spectra (in the presence of deuteration) of c4 and c4-3G in the presence and absence of maltose at 18 and 37 C, respectively. We compared these eight spectra with published chemical shift data of MBP37 (BMRB 4,986 ligand free and BMRB 4,987 containing maltotriose) and BLA38 (BMRB 6,357) as previously described7. The c4 and c4-3G spectra were consistent with the expectation that the individual domain structures were substantially conserved from MBP and BLA. As a quantication of this conservation of structure, we determined a peak count by counting peaks in the proteins spectrum that appeared at the same location in the published peak assignments of BLA and MBP within a tolerance of 0.2 p.p.m. on the 15N axis and 0.06 p.p.m. on the 1H axis. About 55 peaks in the spectra of c4 and c4-3G did not have corresponding peaks in spectra of MBP or BLA, indicating that the residue may experience a different chemical environment than in the parental protein. Peaks present in spectra of the parental proteins but missing from the fusion proteins can result from the residue being in a different chemical environment, from peak broadening, resulting from the residue sampling a broader array of environments, and experimental limitations.

Of the 634 residues of c4 and c4-3G, 579 residues had published peak assignments in the spectra of BLA38 and MBP37.

The peak counts of c4 and c4-3G (relative to those 579 potentially observable peaks) were very similar across all eight experiments (6369%) with one exception (Fig. 6a, Supplementary Figs 24 and Supplementary Table 7). The peak count of c4-3G in the absence of maltose at 37 C was only 43%. We compared the c4 and c4-3G spectra resonances at 37 C in the absence of maltose by overlaying the two spectra (Fig. 6b). Since the spectra for such a large protein is crowded near the centre, we focused our attention on resonances around the periphery (Fig. 6b) to reduce

No maltose Maltose

100

Mapped residues (%)

100

18 C

37 C 18 C

37 C

c4 c4-3G

534 529 385

31 436

473

252

105

110

115

120

125

130

243

474

558

152

201

410

465 413 492

472 564

572

249

437

537

543 549

54 417

516

367

364

15 N (p.p.m.)

468

391 511 346 545

428 507

440 357

518

533

538

359

459 430 427

343

541 457 349

362

542

478

366 182

377

339

250 432

181

461

525

469

425

171

387 172

443

341

10 9

1H (p.p.m.)

8 7 5

5 0 5 10 15 20 600

G(kcal mol1 )

5 0 5 10 15 20

G(kcal mol1 )

GGG linker

Figure 6 | Structural validation of protein switches. (a) Percentage of residues mapped from the 15N-TROSY-HSQC spectra of c4 and c4-3G onto the spectra of MBP37 and BLA38. (b) An overlay of the c4 (black) and c4-3G (red) 15N-TROSY-HSQC spectra at 37 C in the absence of maltose.

The labels represent the peripheral peaks in the c4 spectrum that are not present in the c4-3G spectrum. These missing peaks were manually examined and attributed to residues of MBP (underlined italic) and BLA (regular text) domains. The central region of the spectra that was excluded from examination is shaded in grey. The corresponding missing residues are (c) mapped onto the MBP and BLA domain of the model structure of c4-3G and highlighted in red. Green indicates residues with peripheral resonances present in the spectra of both c4 and c4-3G.

100

0 600

500

400

300

200

100

200

300

400

500

MBP

BLA

MBP

BLA

Residue Residue

Figure 5 | Predicted residue stability using COREX/BEST. The predicted residue stabilities of (a) c4 (black) and c4-3G (red) and (b) c4-3hx in acidic (black) and basic (red) conditions are shown. These stabilities are based on the ratio of the summed probabilities of the states in the ensemble in which a particular residue is in a folded conformation to the summed probability of the states in which a residue is in an unfolded conformation. Arrowheads indicate fusion/linker sites.

MBP domain BLA domain

NATURE COMMUNICATIONS | 6:6968 | DOI: 10.1038/ncomms7968 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms7968

the frequency of misattributing resonances in the fusion protein in our mapping process. Although, some of these peripheral resonances may be misattributed in our process, the good match between resonances of c4 and the set of resonances of MBP and BLA suggests that most are correctly mapped, and that we can draw conclusions from the overall statistics. Of the 189 peripheral resonances mapped onto the sequence of c4, 107 were attribute to the MBP domain and 82 to the BLA domain. Seventy one of these 189 peaks were absent in the c4-3G spectrum (Fig. 6b). The BLA domain of c4-3G had much higher fraction of the absent peaks than the MBP domain (71 and 12%, respectively) (Fig. 6c). In addition, peaks throughout the entire BLA domain displayed broadening and low signal intensity for c4-3G at 37 C in the absence of maltose (Fig. 6b and Supplementary Fig. 2b). Although this analysis is a selective examination of only residues with resonances in the periphery, this result is consistent with the BLA domain of c4-3G experiencing increased conformational entropy and exibility only in the absence of maltose and only at 37 C. Although similar selective broadening in fusion proteins might result from transient interactions between the domains or anisotropic tumbling, such effects would have to occur selectively in c4-3G at 37 C in the absence of maltose for it to explain the lower peak count and resonance broadening under these conditions. Thus, these nuclear magnetic resonance (NMR) experiments provide good evidence that the introduction of conformational exibility in the linker caused switch behavior through an increase in the conformational entropy of the BLA domain; however, more extensive NMR experiments will be necessary to solidify this interpretation.

We attempted to acquire NMR spectra of c4 and c4-3hx at 37 C at pH 6 and pH 8; however, the proteins precipitated in solution. We were able to obtain spectra at 25 C (Supplementary Fig. 5). Most identiable residues of c4 were also observed in c4-3hx. This indicates the preservation of structural integrity despite the replacement of the DKT linker with the relatively long and bulky (EAQA)3 linker. No signicant differences between the spectra of c4 and c4-3hx were apparent. Consistent with this result, we did not observe a switching phenotype when cells expressing c4-3hx were incubated at 25 C (Supplementary Table 8), unlike what we observed at 37 C (Fig. 4b).

Protein unfolding studies further supported the hypothesis that the introduction of conformational exibility in the linker induced conformational heterogeneity in the BLA domain in the absence of maltose. In temperature-induced unfolding studies of c4 in the absence of maltose, the BLA domain unfolded reversibly at a lower temperature than the MBP domain with a transition midpoint (Tm) slightly lower that reported for wild-type BLA39 (Fig. 7a and Supplementary Table 9). The addition of maltose shifted the unfolding of BLA to a higher temperature but

rendered the unfolding irreversible. The substitution of the three glycines in c4 to create c4-3G decreased the Tm of the BLA domain in the absence of maltose but not in the presence of maltose (Fig. 7a and Supplementary Table 9). Similar to c4, the unfolding of the BLA domain of c4-3G was reversible only in the absence of maltose. Differential scanning calorimetry (DSC) experiments on the BLA domain in the absence of maltose revealed that the introduction of glycines in the linker caused the unfolding of the BLA domain to shift from a two-state, cooperative process to a non-two-state process with three possible unfolding domains (Fig. 8a), including one with a Tm4.5 C lower than that of c4 (Supplementary Table 10).

We also conducted unfolding studies on c4 and c4-3hx at pH 6

and 8 (Fig. 7b,c and Supplementary Table 9). Introduction of the pH-responsive linker into c4 had opposite effects on stability at the two pH values. Under conditions where the linker is more structurally ordered (pH 6.0), the introduction of the linker stabilized the protein, while at pH 8.0, when the linker is less ordered, the pH-responsive linker destabilized the protein. For c4-3hx, shifting the pH from 6.0 to 8.0 caused a much larger decrease in Tm in the absence of maltose ( 16.30.3 C) than in

the presence ( 5.41.6 C). DSC experiments at pH 8 on the

BLA domain of c4 and c4-3hx in the absence of maltose mirrored that of c4 and c4-3G described above. The introduction of the linker caused a transition from two-state unfolding to non-two-state unfolding with unfolding domains with considerably lower Tm values than that of c4 (Fig. 8b and Supplementary Table 10), suggesting that multiple conformational states likely exist but with signicantly lower stability.

DiscussionThe EAM predicts that allostery can be established by remodelling the conformational ensemble through the introduction of conformational entropy in the absence of effector. We tested this prediction through the introduction of conformational exibility in the linker region between a fusion of MBP and BLA. The resulting allosteric behavior of the proteins and the switch phenotype conferred to cells by the corresponding genes depended on the degree of conformational exibility in the linker as modulated by linker composition or by environmental conditions (that is, temperature and pH). These results support our hypothesis in numerous ways on how the introduction of conformational exibility can create allosteric switches. First, the introduction of conformational exibility resulted in switch properties in vivo and in vitro. Second, switching depended on conformational exibility, as reducing the temperature (for the glycine-substituted switches) or inducing helix formation by reducing the pH (for EAQA switches) eliminated the switching

a b c

1.00.80.60.40.2 0

c4 versus c4-3G

pH 7.2

c4 pH 6

versus pH 8

c4-3hx pH 6

versus pH 8

Fraction folded or

relative activity

Fraction folded

Fraction folded or

relative activity

Temperature (C)

Temperature (C) Temperature (C)

Figure 7 | Temperature unfolding of conformational switches. (a) Temperature unfolding transitions of c4 (solid red), c4 with 5 mM maltose (solid blue), c4-3G (dotted pink) and c4-3G with 5 mM maltose (dotted light blue) as assessed by the far-ultraviolet CD spectra at wavelength 222 nm. Open circles represent relative enzymatic activity of c4 as a function of temperature in the absence of maltose. Temperature unfolding transitions of (b) c4 and (c) c4-3hx in pH 6 (dotted green), pH 6 with 5 mM maltose (solid green), pH 8 (dotted dark green) and pH 8 with 5 mM maltose (solid dark green) as assessed by the far-ultraviolet CD spectra at wavelength 222 nm. Open circles represent relative enzymatic activity of c4-3hx as a function of temperature in the absence of maltose.

6 NATURE COMMUNICATIONS | 6:6968 | DOI: 10.1038/ncomms7968 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms7968 ARTICLE

Cp (kcal mol1 C1 )

c4-3hx

c4-3G

10 20 30 40 50 60 20 30 40 50 60

Temperature (C)

40 30 20 10

0 10

Temperature (C)

Figure 8 | Thermodynamic validation of protein switches. DSC thermograms in the absence of maltose for (a) c4 and c4-3G at pH 7.2, and (b) c4 and c4-3hx at pH 8. The thermograms shown correspond to the BLA domain and are representative of three independent experiments.

property. Third, switching arose as predicted by the model of Fig. 1 through a reduction in enzyme activity in the absence of maltose and not through an increase in enzyme activity in the presence of maltose. Fourth, the degree of switching monotonically increases with the degree of conformational exibility introduced through the number of glycines or through the degree to which the pH-responsive peptide is expected to exhibit random coil behavior as controlled by pH. Fifth, biophysical studies support the hypothesis that linker exibility induced an increase in the conformational heterogeneity of the BLA domain, which presumably allowed less-active states to be populated more frequently.

Since we created allosteric proteins using two very different methods of introducing exibility applied to the same fusion protein, the molecular details of how conformational exibility is induced is not of primary concern. We propose that our approach potentially represents a general strategy for the introduction of allostery in non-allosteric proteins with the prerequisite effector-binding and functional properties. We note that the switches we created are multi-input switches, requiring a particular environmental condition (temperature or pH) and the presence of an effector (maltose) for activation.

The 23-fold allosteric effect observed in vitro in our engineered proteins is comparable to the magnitude of effect in many natural allosteric enzymes. For example, E. coli aspartate transcarbamoylase, a classic example of an allosteric enzyme, exhibits only B4-fold changes in catalytic activity in the presence of saturating ATP and CTP40. Furthermore, the biological effect our switches provide to cells (up to 32-fold effect on the resistance phenotype) is certainly large enough to be relevant to cell engineering.

We designed our proteins to have their switching property activated at a certain temperature or at a certain pH by the induction of disorder. Such activation has analogues in natural proteins. The yeast protein Hsp26 and the E. coli protein HdeA are chaperones that function to help the cell alleviate the negative effects of high temperature and low pH, respectively. Hsp26 is activated by an increase in the temperature41, and this activation is characterized by a decrease in ordered secondary structure42. HdeA shifts from its inactive, well-folded dimer to its active, partially unfolded monomer upon exposure to low pH43.

Our approach in inducing allostery is counterintuitive from a mechanical/structural view of allostery. From that viewpoint, increased linker exibility might be expected to decouple or dampen the transmittal of the signal from one domain to the other. However, theoretical28,29 and empirical44,45 studies have shown that allostery can be a property of proteins with intrinsically disordered domains. Thus, a structured linker is neither necessary nor sufcient for signal transmission. Whereas the aforementioned studies show that proteins with intrinsically disordered domains can possess allosteric properties, our

experiments demonstrate something different. Our experiments demonstrate that allostery can emerge from the introduction of short, unstructured protein segments at the interface between two domains. In so doing, our work adds to the growing appreciation of how intrinsically disordered regions can contribute to protein function.

Methods

Materials. All chemicals used were purchased from Sigma (St Louis, MO) unless otherwise noted. NMR tubes were purchased from Shigemi (Allison Park, PA). Nitrocen was purchased from EMD Millipore (Billerica, MA). BL21-competent cells were purchased from Agilent Technologies. Plasmids pDIMC8-c4h and pDIMC8-MBP that contain the c4 and malE genes, respectively, were previously described9,10. Plasmid pDIMC8-c4h was used as a template plasmid DNA to make linker variants using the PCR-based QuikChange mutagenesis procedure (Agilent Technologies). DNA sequencing of 2,0132,061 bp regions containing the targeted codon conrmed mutant clones (Genewiz). The WA-(EAQA)7 peptide (sequence:

WAEAQAEAQAEAQAEAQAEAQAEAQAEAQA) was custom synthesized with 98% purity by Peptide 2.0.synthesis, and peptide ends were modied by N-terminal acetylation and C-terminal amidation.

Molecular modeling and simulation. The molecular models of c4 and the c4 variants were constructed with the amino-acid sequence and coordinatesof MBP (PDB ID: 1OMP) and BLA (PDB ID: 1M40), and the distance constraints of the engineered MBP-BLA fusion protein, RG13 (PDB ID: 4DXC), using MODELLER software46. The sequence of residues 1316 and 319370 of MBP and residues 26290 of BLA were aligned with the sequence of RG13, and the initial model of c4 was built using RG13 as a template. Since there was no structural information about the linker between the MBP and BLA domains, intermolecular distance between two domain obtained from RG13 were used as constraints to construct the initial models, and then ab initio loop optimization methods were used to model the linkers between residues 313 and 323 and 575 and 585. The linker-optimized model was then used as a template to optimize the model ofc4 by repeating the steps described above. For c4-3G model, the amino-acid residues 580582 were replaced with three glycine residues. For c4-3hx model, the residues 580582 were replaced with a three-coiled a-helix and random coil structure model of the (EAQA)3 peptide. The three-coiled a-helix model of EAQA (3hx) was built using UCSF Chimera47, and the random coil model was the nal structure of the (EAQA)3 model after 100 ns of MD simulation at pH 7 (described below). The model-building step was then repeated again with the (EAQA)3 inserted models as templates to optimize the intermolecular distance between domains. The side-chain positions of the c4 and c4 variant models were subsequently optimized by using SCWRL 4.1 (ref. 48), and the nal models were optimized by energy minimization using the GROMACS 4.1.3 package49.

COREX/BEST50 was used for statistical thermodynamic analysis of the c4 and c4 variant structure models at 37 C. The model structures were used as templates to generate 4105 partially unfolded states using Monte Carlo sampling strategy with the assumption that folded regions maintain native-like geometry and unfolded regions are modelled as fully solvated16. Folding units were dened by assigning sequential residues to different folding units (11 residues per unit), and different states were enumerated by folding and unfolding each unit. Folding units were varied by sliding the boundaries by one-residue increments and repeating the enumeration. The free energy of each state was determined using a surface area-based parameterization of the energetic. Each of c4 and c4 variant model was then characterized with the residue stability constant and relative free energy, which were calculated from the ratio of the summed probabilities of the states in the enumerated protein ensembles in which a particular residue is in a folded conformation to the summed probability of the states in which a residue was in an unfolded conformation16,50.

MD simulations of EAQA peptide. All simulations of (EAQA)3 peptide were performed with the GROMACS software package49 using the GROMMOS 43a3 force eld51. The model structures used in the simulations were the simple a-helix model of (EAQA)3 peptide. Simple EAQA models were solvated individually in octahedron boxes lled with water molecules52. A single-point charge water model was used for the solvent molecules in the simulation53. Sodium ions were used to electroneutralize the system. Solutes, solvent and counterions were coupled independently to reference temperature baths at 300 K (ref. 54), and the pressure was maintained by coupling the system weakly to an external pressure bath at1 atm. Bond lengths were constrained by the LINCS procedure55 and non-bonded interactions were evaluated using twin-range cutoff of 0.8 and 1.4 nm for LennardJones and Coulomb potentials. The long-range electrostatic interactions beyond the cutoff were treated with the generalized reaction eld model by using a dielectric constant of 54. The integration time step was set to 0.002 ps and the trajectory coordinates and energies were stored at 0.5 ps intervals. The analysis was performed using the built-in programs of GROMACS software package.

To emulate the different pH environment, glutamate residues and the C terminus

NATURE COMMUNICATIONS | 6:6968 | DOI: 10.1038/ncomms7968 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms7968

were protonated accordingly. Analysis was performed using the built-in programs of the GROMACS software package49.

Expression and purication. Proteins were expressed and puried from BL21 cells grown in M9 minimal media containing ammonium chloride (18.5 mM), thiamine (30 mM), MgSO4 (2 mM), CaCl2 (100 mM), fructose (50 mM), 2% glycerol (w/v) and chloramphenicol (50 mg ml 1). For production of each protein, 1 lof minimal media was inoculated with 10 ml overnight culture of BL21 cells, harbouring the pDIMC8 plasmid encoding the proteins, and shaken at 37 C until the OD600 reached 0.7. The culture was induced with 1 mM isopropyl-b-D-thiogalactopyranoside (IPTG) and incubated at 25 C for another 4850 h. For NMR study, proteins were prepared from M9/D2O minimal media containing25 mM deuterated glucose-d7 (25 mM) and15NH4Cl (18.5 mM) as described in the previous study7. For production of labelled proteins in 1 l minimal media, a 10-ml overnight culture was prepared in lysogeny broth (LB)-rich media at 37 C. The cells were then pelleted, washed and resuspended in 1 ml of D2O. The resuspension was diluted to 10 ml of labelled M9/D2O minimal media and grown overnight at 37 C. The overnight culture was diluted to 100 ml with labelled M9/D2O minimal media and shaken at 37 C. Once the OD600 reached between 0.7, the cells were transferred to 2 l ask and diluted again to 1 l of labelled M9/D2O minimal media. The 1-l culture was incubated at 37 C until the OD600 reached 0.7, at which time IPTG was added to the culture to a nal concentration of 1 mM to induce the proteins of interest. The 1 l culture was switched to a 25-C shaker for 72 h.

After expression, the cells were harvested and resuspended in 10 ml buffer (50 mM sodium phosphate and 150 mM NaCl, pH 7.2) and lysed using a French press. Cell lysates were centrifuged at 20,000g for 1 h. The soluble proteins were puried using the HisTrap (GE Healthcare Life Sciences) and size-exclusion column in the AKTA FPLC purier system (GE Healthcare Life Sciences). Proteins were dialysed against 15 l of a dialysis buffer (25 mM phosphate, pH 7.2) and concentrated using Amicon centrifugal lter (Millipore). For the NMR study, NaN3 was added to 0.05%. Approximately 12 mg of puried proteins were obtained with 495% purity, estimated by Coomassie blue staining of SDSPAGE gels. Protein concentration was determined using a NanoDrop spectrophotometer (Thermo Scientic) by measuring ultraviolet absorbance at 280 nm.

Nitrocen hydrolysis enzymatic assay. Enzymatic assays were performed as previously described9,10 Proteins were added to 100 mM sodium phosphate buffer, pH 7.2, to the nal concentration of 15 nM, and were incubated in the absence and presence of 5 mM maltose for 530 min at 18 and 37 C. For pH experiment, c4 and c4-3hx were incubated in pH-adjusted phosphate buffer (pH 6 and 8)for 30 min before the incubation with maltose. Nitrocen was added to the nal concentration of 100 mM. The absorbance at the wavelength of 486 nm was recorded by using a Varian Cary 50 UV-Visible Spectrophotometer (Agilent Technologies) and the initial rate of nitrocen hydrolysis reaction was measured by Cary WinUV software package.

MICAmp assay. For c4 and the c4 variants, the MICAmp assay was performed as previously described9,10. BL21 cells were transformed with pDIMC8 plasmids containing genes that encode the protein variants. For each overnight culture, 5 ml of tryptone broth (10 g l 1 tryptone and 10 g l 1 NaCl) was inoculated by picking a single colony and incubated in a 37-C shaker for 1618 h. The optical density at 600 nm (OD600nm) was measured. For the MICAmp assay, B1 106 colony-

forming units (based on OD600nm) from the overnight culture were added to 1 or 5 ml of tryptone medium containing chloramphenicol (50 mg ml 1), IPTG(300 mM), Amp (08,192 mg ml 1), and either the absence or presence of 5 mM maltose in each culture tube. For temperature-dependent MICAmp assay, 5 ml cultures were incubated in 18 and 37 C shaker incubators for 36 and 18 h, respectively, and the OD600nm of each culture was measured. For pH-dependent

MICAmp assay, pH of the tryptone medium was adjusted with HCl and NaOH, and 1 ml cultures in 96-well assay block (Corning, Tewksbury, MA) were incubated in 37 C shaker incubators for 18 h. No growth was observed at pH 5. The MICAmp

was dened as the lowest concentration of Amp at which the OD600nm was o5% of

the OD600nm in the absence of Amp.

CD and uorescence spectroscopy. CD and uorescence emission spectra of the WA-(EAQA)7 peptide was scanned at 23 C as previously described35. The WA-(EAQA)7 peptide was dissolved in 50 mM TES at pH 7.5 (0.5 mg ml 1). The stock solution was diluted in a 5-mM 2-[[1,3-dihydroxy-2-(hydroxymethyl)

propan-2-yl]amino]ethanesulfonic acid (TES)/100 mM KCl, with varying pH by NaOH and HCl, and incubated for 30 min before the scan. The far-ultraviolet CD spectra were recorded in a Jasco J-710 CD Spectropolarimeter (Jasco, Easton, MD), averaged over three scans at a spectral bandwidth of 1 nm with 2 s integration time using a quartz optical cell with an optical pathlength of 1 mm (Fisher Scientic). Fluorescence emission spectra were recorded in a PTI Quantamaster30 Fluorescence Spectrouorometer (Photon Technology International) with excitation at 280 nm and spectral bandwidths of 2 nm (excitation) and 4 nm (emission) using quartz cuvette (pathlength of 10 mm). The nal concentration was 0.05 mg ml 1 for CD spectra and 0.005 mg ml 1 for uorescence emission spectra.

For c4 and c4 variants, thermal denaturation was examined by monitoring the far-ultraviolet CD spectra. Protein samples were incubated in the absence and presence of 5 mM maltose for 1 h at 4 C in 50 mM phosphate buffer (pH 7.2) before thermal denaturation. To obtain the CD spectra at varying pH, the proteins were prepared in pH-adjusted 50 mM phosphate buffer (pH 6 and 8) before the incubation. The nal protein concentration was 5 mM. CD absorbance at 222 nm was collected from 20 to 90 C at a temperature ramp rate of 1 C min 1 and a bandwidth of 1 nm. All spectra were buffer corrected. The data were normalized and tted to a two-state and three-state unfolding model using Prism analysis software (GraphPad, La Jolla, CA).

Differential scanning calorimetry. For the calorimetric study of protein unfolding, the protein samples were dialysed overnight against a buffer containing 50 mM sodium phosphate and 5% w/v glycerol at pH 7.2 and 8, which was also used in the reference cell and for baseline determination. The nal protein concentrations were 1 mg ml 1. Heat absorbance was collected from 20 to 55 C in a MicroCal VP-DSC instrument at a temperature ramp rate of 1 C min 1.

Thermograms were baseline corrected, normalized and analysed according to two-state and non-two-state models with a single and multiple transitions, in which Tm, DHcal and DHvH, and DCp of individual transitions were tted independently using the MicroCal Origin software (Microcal, Northhampton, MA).

NMR spectral acquisition. Acquisition of TROSY-HSQC spectra was conducted as previously described7. D2O was added to 10% to a 200-mM sample of c4 and c4 variants. The 350 ml samples were then loaded into a Shigemi tube. Data were acquired on an 800-MHz (1H) Varian INOVA spectrometer equipped with a cryogenic probe and z axis gradients at appropriate temperatures (18, 25 and37 C) in the presence and absence of 5 mM maltose. A sensitivity-enhanced TROSY sequence56 was used, with 16 scans/FID and an interscan delay of 1.9 s. Spectral widths of 16 kHz (1H) and 2,836 Hz (15N) were used, with an acquisition time of 60 ms in both 1H and 15N dimensions. The total data collection time was3.25 h per spectrum. Quadrature detection in the 15N dimension was achieved using the echoantiecho method57. The data were apodized with 60 phase-shifted cosine-squared bell window functions in both dimensions and zero-lled to a nal digital resolution of 4 Hz pt 1 in 1H and 1.4 Hz pt 1 in 15N. Data were processed and analysed using NMRPipe and NMRView software (One Moon Scientic,

Westeld, NJ).

References

1. Ostermeier, M. Engineering allosteric protein switches by domain insertion. Protein Eng. Des. Sel. 18, 359364 (2005).

2. Stratton, M. M. & Loh, S. N. Converting a protein into a switch for biosensing and functional regulation. Protein Sci. 20, 1929 (2011).

3. Vallee-Belisle, A. & Plaxco, K. W. Structure-switching biosensors: inspired by Nature. Curr. Opin. Struct. Biol. 20, 518526 (2010).

4. Wright, C. M., Wright, R. C., Eshleman, J. R. & Ostermeier, M. A protein therapeutic modality founded on molecular regulation. Proc. Natl Acad. Sci. USA 108, 1620616211 (2011).

5. Guntas, G., Mitchell, S. F. & Ostermeier, M. A molecular switch created by in vitro recombination of nonhomologous genes. Chem. Biol. 11, 14831487 (2004).

6. Guntas, G., Mansell, T. J., Kim, J. R. & Ostermeier, M. Directed evolution of protein switches and their application to the creation of ligand-binding proteins. Proc. Natl Acad. Sci. USA 102, 1122411229 (2005).

7. Wright, C. M., Majumdar, A., Tolman, J. R. & Ostermeier, M. NMR characterization of an engineered domain fusion between maltose binding protein and TEM1 beta-lactamase provides insight into its structure and allosteric mechanism. Proteins 78, 14231430 (2009).

8. Kim, J. R. & Ostermeier, M. Modulation of effector afnity by hinge region mutations also modulates switching activity in an engineered allosteric TEM1 beta-lactamase switch. Arch. Biochem. Biophys. 446, 4451 (2006).

9. Heins, R. A., Choi, J. H., Sohka, T. & Ostermeier, M. In vitro recombination of non-homologous genes can result in gene fusions that confer a switching phenotype to cells. PLoS One 6, e27302 (2011).

10. Choi, J. H., San, A. & Ostermeier, M. Non-allosteric enzyme switches possess larger effector-induced changes in thermodynamic stability than their non-switch analogs. Protein Sci. 22, 475485 (2013).

11. Motlagh, H. N., Wrabl, J. O., Li, J. & Hilser, V. J. The ensemble nature of allostery. Nature 508, 331339 (2014).

12. Cui, Q. & Karplus, M. Allostery and cooperativity revisited. Protein Sci. 17, 12951307 (2008).

13. Monod, J., Wyman, J. & Changeux, J. P. On the nature of allosteric transitions: a plausible model. J. Mol. Biol. 12, 88118 (1965).

14. Koshland, Jr. D. E. Enzyme exibility and enzyme action. J. Cell. Comp. Physiol. 54, 245258 (1959).

15. Changeux, J. P. Allostery and the Monod-Wyman-Changeux model after 50 years. Annu. Rev. Biophys. 41, 103133 (2012).

8 NATURE COMMUNICATIONS | 6:6968 | DOI: 10.1038/ncomms7968 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms7968 ARTICLE

16. Hilser, V. J., Garcia-Moreno, E. B., Oas, T. G., Kapp, G. & Whitten, S. T.

A statistical thermodynamic model of the protein ensemble. Chem. Rev. 106, 15451558 (2006).

17. Popovych, N., Sun, S., Ebright, R. H. & Kalodimos, C. G. Dynamically driven protein allostery. Nat. Struct. Mol. Biol. 13, 831838 (2006).

18. Boehr, D. D., Nussinov, R. & Wright, P. E. The role of dynamic conformational ensembles in biomolecular recognition. Nat. Chem. Biol. 5, 789796(2009).

19. Schrank, T. P., Bolen, D. W. & Hilser, V. J. Rational modulation of conformational uctuations in adenylate kinase reveals a local unfolding mechanism for allostery and functional adaptation in proteins. Proc. Natl Acad. Sci. USA 106, 1698416989 (2009).

20. Tzeng, S. R. & Kalodimos, C. G. Dynamic activation of an allosteric regulatory protein. Nature 462, 368372 (2009).

21. Wrabl, J. O. et al. The role of protein conformational uctuations in allostery, function, and evolution. Biophys. Chem. 159, 129141 (2011).

22. Tzeng, S. R. & Kalodimos, C. G. Protein dynamics and allostery: an NMR view. Curr. Opin. Struct. Biol. 21, 6267 (2011).

23. Tzeng, S. R. & Kalodimos, C. G. Protein activity regulation by conformational entropy. Nature 488, 236240 (2012).

24. Hilser, V. J. Structural biology: signalling from disordered proteins. Nature 498,

308310 (2013).

25. Tzeng, S. R. & Kalodimos, C. G. Allosteric inhibition through suppression of transient conformational states. Nat. Chem. Biol. 9, 462465 (2013).

26. Cooper, A. & Dryden, D. T. Allostery without conformational change. A plausible model. Eur. Biophys. J. 11, 103109 (1984).

27. Petit, C. M., Zhang, J., Sapienza, P. J., Fuentes, E. J. & Lee, A. L. Hidden dynamic allostery in a PDZ domain. Proc. Natl Acad. Sci. USA 106, 1824918254 (2009).

28. Motlagh, H. N. & Hilser, V. J. Agonism/antagonism switching in allosteric ensembles. Proc. Natl Acad. Sci. USA 109, 41344139 (2012).

29. Hilser, V. J. & Thompson, E. B. Intrinsic disorder as a mechanism to optimize allosteric coupling in proteins. Proc. Natl Acad. Sci. USA 104, 83118315 (2007).

30. Hilser, V. J., Wrabl, J. O. & Motlagh, H. N. Structural and energetic basis of allostery. Annu. Rev. Biophys. 41, 585609 (2012).

31. Wand, A. J. Dynamic activation of protein function: a view emerging from NMR spectroscopy. Nat. Struct. Biol. 8, 926931 (2001).

32. Feller, G., dAmico, D. & Gerday, C. Thermodynamic stability of a cold-active alpha-amylase from the Antarctic bacterium Alteromonas haloplanctis. Biochemistry 38, 46134619 (1999).

33. Holland, L. Z., McFall-Ngai, M. & Somero, G. N. Evolution of lactate dehydrogenase-A homologs of barracuda shes (genus Sphyraena) from different thermal environments: differences in kinetic properties and thermal stability are due to amino acid substitutions outside the active site. Biochemistry 36, 32073215 (1997).

34. Fields, P. A. & Somero, G. N. Hot spots in cold adaptation: localized increases in conformational exibility in lactate dehydrogenase A4 orthologs of Antarctic notothenioid shes. Proc. Natl Acad. Sci. USA 95, 1147611481 (1998).

35. Subbarao, N. K., Parente, R. A., Szoka, Jr. F. C., Nadasdi, L. & Pongracz, K. pH-dependent bilayer destabilization by an amphipathic peptide. Biochemistry 26, 29642972 (1987).

36. Li, W., Nicol, F. & Szoka, Jr. F. C. GALA: a designed synthetic pH-responsive amphipathic peptide with applications in drug and gene delivery. Adv. Drug Deliv. Rev. 56, 967985 (2004).

37. Evenas, J. et al. Ligand-induced structural changes to maltodextrin-binding protein as studied by solution NMR spectroscopy. J. Mol. Biol. 309, 961974 (2001).

38. Savard, P. Y. & Gagne, S. M. Backbone dynamics of TEM-1 determined by NMR: evidence for a highly ordered protein. Biochemistry 45, 1141411424 (2006).

39. Kather, I., Jakob, R. P., Dobbek, H. & Schmid, F. X. Increased folding stability of TEM-1 beta-lactamase by in vitro selection. J. Mol. Biol. 383, 238251 (2008).

40. West, J. M., Tsuruta, H. & Kantrowitz, E. R. A uorescent probe-labeled Escherichia coli aspartate transcarbamoylase that monitors the allosteric conformational state. J. Biol. Chem. 279, 945951 (2004).

41. Haslbeck, M. et al. Hsp26: a temperature-regulated chaperone. EMBO J. 18, 67446751 (1999).

42. Franzmann, T. M., Wuhr, M., Richter, K., Walter, S. & Buchner, J.

The activation mechanism of Hsp26 does not require dissociation of the oligomer. J. Mol. Biol. 350, 10831093 (2005).

43. Hong, W. et al. Periplasmic protein HdeA exhibits chaperone-like activity exclusively within stomach pH range by transforming into disordered conformation. J. Biol. Chem. 280, 2702927034 (2005).

44. Ferreon, A. C., Ferreon, J. C., Wright, P. E. & Deniz, A. A. Modulation of allostery by protein intrinsic disorder. Nature 498, 390394 (2013).

45. Garcia-Pino, A. et al. Allostery and intrinsic disorder mediate transcription regulation by conditional cooperativity. Cell 142, 101111 (2010).

46. Sali, A. & Blundell, T. L. Comparative protein modelling by satisfaction of spatial restraints. J. Mol. Biol. 234, 779815 (1993).

47. Pettersen, E. F. et al. UCSF Chimera--a visualization system for exploratory research and analysis. J. Comput. Chem. 25, 16051612 (2004).

48. Krivov, G. G., Shapovalov, M. V. & Dunbrack, Jr. R. L. Improved prediction of protein side-chain conformations with SCWRL4. Proteins 77, 778795 (2009).

49. Van Der Spoel, D. et al. GROMACS: fast, exible, and free. J. Comput. Chem. 26, 17011718 (2005).

50. Vertrees, J., Barritt, P., Whitten, S. & Hilser, V. J. COREX/BEST server: a web browser-based program that calculates regional stability variations within protein structures. Bioinformatics 21, 33183319 (2005).

51. Daura, X., Mark, A. E. & Van Gunsteren, W. F. Parametrization of aliphatic CHn united atoms of GROMOS96 force eld. J. Comput. Chem. 19, 535547 (1998).

52. Eisenberg, D. & McLachlan, A. D. Solvation energy in protein folding and binding. Nature 319, 199203 (1986).

53. Berendsen, H. J., Postma, J. P., van Gunsteren, W. F., DiNola, A. & Haak, J. R. in Intermolecular Forces. (ed. Pullman, B.) 331342 (Reidel Publishing Company, 1981).

54. Berendsen, H.J.C.P., J. P. M., van Gunsteren, W. F., Dinola, A. & Haak, J. R. Molecular dynamics with coupling to an external bath. J. Chem. Phys. 81, 36843690 (1984).

55. Hess, B., Bekker, H., Berendsen, H. J. C. & Fraaije, J. G. E. M. LINCS: a linear constraint solver for molecular simulations. J. Comput. Chem. 18, 14631472 (1997).

56. Rance, M., Loria, J. P. & Palmer, A. G. r. Sensitivity improvement of transverse relaxation-optimized spectroscopy. J. Magn. Reson. 136, 92101 (1999).

57. Davis, A. L., Keeler, J., Laue, E. D. & Moskau, D. Experiments for recording pureabsorption heteronuclear correlation spectra using pulsed eld gradients.J. Magn. Reson. 98, 207216 (1992).

Acknowledgements

We thank Francisco M. Goycoolea and Juliette Lecomte for their inspiring ideas and

valuable discussions, James O. Wrabl for assistance with COREX/BEST analysis, Ananya

Majumdar for NMR spectral acquisition, Katherine Tripp for assistance with DSC,

Angela San for preliminary MIC assay experiments and Robert Clay Wright for com

ments on the manuscript. This work was supported by funding from the Defense Threat

Reduction Agency (HDTRA1-09-1-0016), the National Institute of General Medicine at

the National Institutes of Health (R01 GM066972), the National Science Foundation

(MCB-0919377) to M.O., from the National Institute of General Medicine at the

National Institutes of Health (R01 GM63747) and the National Science Foundation

(MCB1330211) to V.J.H.

Author contributions

J.H.C. and M.O. designed the study; J.H.C. performed the experiments; A.H.L. assisted

J.H.C. with the NMR experiments; J.H.C., V.J.H. and M.O. performed data analysis and

interpretation; and J.H.C., V.H. and M.O. wrote the paper.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/

http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing nancial interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/

http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Choi, J.H. et al. Design of protein switches based on an

ensemble model of allostery. Nat. Commun. 6:6968 doi: 10.1038/ncomms7968 (2015).

NATURE COMMUNICATIONS | 6:6968 | DOI: 10.1038/ncomms7968 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 9

Word count: 9698

Show less

Abstract

Translate

Switchable proteins that can be regulated through exogenous or endogenous inputs have a broad range of biotechnological and biomedical applications. Here we describe the design of switchable enzymes based on an ensemble allosteric model. First, we insert an enzyme domain into an effector-binding domain such that both domains remain functionally intact. Second, we induce the fusion to behave as a switch through the introduction of conditional conformational flexibility designed to increase the conformational entropy of the enzyme domain in a temperature- or pH-dependent fashion. We confirm the switching behaviour in vitro and in vivo. Structural and thermodynamic studies support the hypothesis that switching result from an increase in conformational entropy of the enzyme domain in the absence of effector. These results support the ensemble model of allostery and embody a strategy for the design of protein switches.

Details

Title

Design of protein switches based on an ensemble model of allostery

Author

Choi, Jay H; Laurent, Abigail H; Hilser, Vincent J; Ostermeier, Marc

Pages

6968

Publication year

2015

Publication date

Apr 2015

Publisher

Nature Publishing Group

e-ISSN

20411723

Source type

Scholarly Journal

Language of publication

English

DOI

https://doi.org/10.1038/ncomms7968

ProQuest document ID

1674733156

Design of protein switches based on an ensemble model of allostery

Jump to:

Full Text

Abstract

Details

Suggested sources