Andrew M.
Hartley‡
a,
Athraa J.
Zaki‡
b,
Adam R.
McGarrity‡
a,
Cecile
Robert-Ansart
a,
Andriy V.
Moskalenko
b,
Gareth F.
Jones
c,
Monica F.
Craciun
c,
Saverio
Russo
c,
Martin
Elliott
b,
J. Emyr
Macdonald
b and
D. Dafydd
Jones
*a
aSchool of Biosciences, Cardiff University, Cardiff, UK. E-mail: jonesdd@cardiff.ac.uk
bSchool of Physics and Astronomy, Cardiff University, Cardiff, Wales, UK
cCentre for Graphene Science, University of Exeter, Exeter, Devon, UK
First published on 31st March 2015
Post-translational modification (PTM) modulates and supplements protein functionality. In nature this high precision event requires specific motifs and/or associated modification machinery. To overcome the inherent complexity that hinders PTM's wider use, we have utilized a non-native biocompatible Click chemistry approach to site-specifically modify TEM β-lactamase that adds new functionality. In silico modelling was used to design TEM β-lactamase variants with the non-natural amino acid p-azido-L-phenylalanine (azF) placed at functionally strategic positions permitting residue-specific modification with alkyne adducts by exploiting strain-promoted azide–alkyne cycloaddition. Three designs were implemented so that the modification would: (i) inhibit TEM activity (Y105azF); (ii) restore activity compromised by the initial mutation (P174azF); (iii) facilitate assembly on pristine graphene (W165azF). A dibenzylcyclooctyne (DBCO) with amine functionality was enough to modulate enzymatic activity. Modification of TEMW165azF with a DBCO–pyrene adduct had little effect on activity despite the modification site being close to a key catalytic residue but allowed directed assembly of the enzyme on graphene, potentially facilitating the construction of protein-gated carbon transistor systems.
Covalent modification with non-natural adducts is traditionally achieved using chemistry inherent to the natural amino acid repertoire, mostly amine, carboxyl and thiol groups. The main problem is lack of specificity as such chemistry is normally distributed across the surface of a protein and is ubiquitous in the proteome. A more general and powerful approach is the introduction of new chemically reactive handles not present in the native 20 amino acid set through the use of an expanded genetic code6 (see ref. 7–10 for recent reviews). The unique reactivity of a chosen non-natural amino acid (nAA) together with the ability to select both the target protein and the residue within the target means that specificity comparable to natural PTM events can be achieved.10,11 The use of highly specialized caged nAAs has been very effective in controlling activity12–15 but is restricted to certain protein types and chemistry, and can add significant bulk to the amino acid side chain during the production and folding of the nascent polypeptide and its subsequent folded form.
The incorporation of phenyl azide chemistry into proteins through the use of p-azidophenylalanine (azF)17 is an attractive alternative. As well as being only one atom bigger than the natural amino acid tyrosine and having been successfully incorporated into a wide variety of proteins,7–9 the phenyl azide moiety opens up different routes to non-natural PTM (nnPTM): photochemical transformations and Click chemistry adduct addition.9 The use of phenyl azide photochemistry to control protein activity has recently been demonstrated.18–22 Azide–alkyne cycloaddition is fast becoming a useful approach for orthogonal biomolecule conjugation but has largely been used in a passive way, for example, to label proteins.7,23–25 Given the versatility in terms of the array of adducts available coupled with the inherent bioorthogonality and biocompatibility,26 it is surprising Click chemistry has not been used more extensively as a general direct modulator of protein activity. Additionally, useful adducts can be placed at strategic positions to expand and facilitate protein function in a manner akin to co-factors. This includes attachment of entities to facilitate interfacing and assembly with secondary systems or materials.27–31 As in many natural biomolecular assemblies, a defined and optimal protein–material interface is critical for maximal communication between the individual elements. Interfacing proteins with carbon sp2 materials such as graphene is gaining significant interest as it forms the basis for constructing hybrid bio-transistors in which events at even the single protein molecule level can be used to gate conductance through graphene.27,32
Here, we show that the activity and assembly onto pristine graphene of the antibiotic resistance protein TEM β-lactamase33,34 can be directly controlled by Click chemistry. By genetically encoding phenyl azide chemistry at designed positions in TEM and using Cu-free biocompatible strain promoted azide–alkyne cycloaddition (SPAAC; Fig. 1)35,36 different adducts can be attached (1, 2 and 3; Fig. 1) that either reduce or restore activity on modification, and define assembly of TEM on graphene.
Fig. 1 Modification of TEM β-lactamase. The structure of TEM β-lactamase (top left) with the inhibitor imipenem (yellow sticks to emphasise the substrate binding pocket) highlights the residues targeted for replacement with azF (red spheres) together with the key catalytic residues (cyan sticks) and the Ω-loop (green). The protein structures were generated using PyMol.16 A representation of the SPAAC reaction using the activated alkyne present in the core DBCO moiety is shown (top right). The R groups attached to the DBCO used in this study are shown (bottom half). |
Fig. 2 Molecular design models of TEM azF variants. (A) Model of TEMY105azF (yellow) aligned with wt TEM β-lactamase (grey; PDB 1btl), with the active site S70 shown as stick representation. (B) Alignment of the TEMP174azF (cyan) with wt TEM (grey; PDB 1btl). The region highlighted in red covers residues 173–175. (C) Interaction of AzF174 with neighbouring residues. |
Fig. 3 Enzyme kinetics of ampicillin hydrolysis by TEM β-lactamase variants pre- and post-Click modification. (A) Catalytic efficiency of each variant before and after modification with 1 (+1) or 2 (+2) derived from KM and kcat (ESI Table 1†). (B) The effect of modifying TEMP174azF with 1 (Click 1) or 2 (Click 2) on Michaelis–Menten kinetics. (C) Modelled changes to backbone structure of TEMP174azF (cyan) on modification with 1 (green) compared to wt TEM (grey) (top panel) and the configuration and local interactions (hydrogen bond in red) of 1 (bottom panel). |
As predicted, TEMY105azF was largely inhibited on modification with either 1 or 2 (Fig. 3) despite neither having any significant inherent inhibitory effect on wild-type TEM (ESI Table 1†). Both catalysis and substrate binding were disrupted as evident by the decrease in kcat and increase in KM (ESI Table 1†) for ampicillin resulting in an overall decrease in catalytic efficiency of ∼70% and 85% when modified with 1 or 2 respectively compared to unmodified protein. This is similar to the most deactivating site-directed mutants of Y105.40 Modelling of the covalent complex (vide infra) between TEMY105azF and 1 suggested that modification blocks access at one end of the substrate binding site, with the amine of 1 forming a hydrogen bond with the backbone carbonyl group of G238 (ESI Fig. 3†). The enzyme kinetics suggests that simple steric blocking of substrate binding is unlikely to be the sole mechanism of action as catalysis and substrate binding are equally affected, similar to that observed for classical mixed inhibition models. The relatively small difference in catalytic efficiency of TEMY105azF modified with either 1 or 2 indicates that the most significant effect is occurring close to the linkage site, as suggested by the covalent complex model (ESI Fig. 3†). However, the bulkier group of 2 is exerting a slightly greater effect with regards to kcat rather than KM (Fig. 3 and ESI Table 1†) suggesting that long range interactions made by the rhodamine dye moiety may be forcing the structure around residue 105 to adapt less catalytically proficient form. Currently, the separation of modified from unmodified protein has proved unsuccessful, so the observed activity may be contributed by unmodified protein only (estimated to be ∼20%; vide supra). Therefore, the inhibition levels observed here represent the lowest currently achievable and we cannot rule out modified species are essentially inactive.
P174 is another largely surface exposed residue (∼80 Å2 ≈ 60% relative surface accessibility) that is relatively distant from the active site (Fig. 1). P174 however contributes towards one of the turns that comprise the Ω loop (Fig. 1 and ESI Fig. 4†), a region critical for substrate binding, substrate specificity and catalysis.42,43 It was hypothesized in the design process that mutating P174 to a non-cyclic amino acid would result in local conformational changes around the turn region. Molecular modelling supported this hypothesis with the structure of the turn slightly shifted compared to wt TEM (Fig. 2b and ESI Fig. 5†). The side chain of azF174 lies approximately perpendicular to the native proline, with the azido moiety fitting into a shallow pocket and making a local polar interaction network with R42 and T265 (Fig. 2c). These new interactions may act as the driver of local conformation changes associated with the P174azF mutation, which in turn alters the local bonding network that subtly shifts the conformation of both the active site and substrate binding residues, including S70 (ESI Fig. 5†). Based on the TEMP174azF model it was further postulated that on modification with 1, the formation of the triazole link would break the interaction network with R42 and T265 so reconstituting the original loop structure and activating the protein. The amine group also has the potential to form new long range hydrogen bonds in spatially local regions (Fig. 3 and ESI Fig. 4†).
As suggested by the model, replacement of P174 with azF results in a significant change in activity. The overall catalytic efficiency of TEMP174azF was 60% lower compared to wt TEM (Fig. 3). The major effect was on kcat which was 3 fold lower compared to wt TEM (ESI Table 1†), supporting the evidence from the model that the likely effect of azF incorporation at residue 174 is to disrupt interactions of residues associated with catalysis. This was offset by a slightly increased ampicillin affinity.
TEMP174azF was accessible to Click modification, but efficiency was lower (estimated ∼33% using the absorbance properties of 2) compared to TEMY105azF. It is unknown how local protein microenvironment dictates SPAAC efficiency but dynamics, relative exposure to aqueous solvent and the character of shallow “pockets” have been suggested as possible determinants.23 The molecular model suggests that interaction of the azide moiety with other residues, which may play a role through, for example, altering the relative populations of resonance structures sampled by the azide moiety (R–NN+N−/R–N−–N+N), and/or its accessibility and available orientations to the incoming DBCO. The effect of the SPAAC nnPTM on the ability of TEMP174azF to hydrolyse ampicillin varied depending on the DBCO adduct. Modification with 2 exerted a similar but somewhat smaller inhibitory effect (∼40% drop) as observed for TEMT105azF (Fig. 3 and ESI Table 1†). Modification with 1 resulted in a significant increase in overall activity restoring apparent catalytic efficiency essentially to wild-type levels (Fig. 3). This should be considered as a lower estimate of activation as ∼70% of the protein may be unmodified (vide supra). With both adducts, the most significant contribution was the change in kcat of TEMP174azF compared to unmodified protein (ESI Table 1†); modification with 1 increased kcat by 210% while addition of 2 reduced kcat by almost half.
To understand how 1 exerts its beneficial effect on TEMP174azF, the nnPTM product was modelled. Using the TEMP174azF model as a starting point, the triazole linkage between azF and DBCO of 1 was generated in silico by parameterizing the azF–DBCO complex to calculate optimized geometries and the electrostatic potentials. The model was subjected to molecular dynamics for a total of 5 ns. The model of nnPTM product, termed TEMP174azF+1 suggested that the general structure was closer to wt TEM than the original TEMP174azF model (Fig. 3c and ESI Fig. 5†). The Ω-loop structure around residue 174 and the local bonding networks for TEMP174azF+1 were largely comparable with wt TEM. The orientation of residue 174 differs significantly to accommodate 1 (ESI Fig. 5†). The amine group of 1 bends back and appears to come within hydrogen bonding distance (∼2.6 Å) of the side-chain carbonyl group of the adjacent residue N175 (Fig. 3c). The loss of the azide group on formation of the trizole removes the interactions with R43 and T238. However, the reason for 1 acting as an activator and 2 as an inhibitor may be down to the reduced bulk and ability of 1 to form local H-bonds.
W165 is another partially surface exposed residue (∼100 Å2 ≈ 50% relative surface accessibility) that fulfils the required criteria (Fig. 1). It resides in the Ω-loop and splits two functionally important residues; E166 is a conserved catalytic residue involved in proton shuttling34 and R164 forms buried ionic interactions critical for stabilizing the Ω loop and thus the positioning of E166.45 The approach used to model TEMP174azF modified with 1 was applied to the W165azF mutation to ascertain the potential position of the adduct group and its effect on structure (Fig. 4a). Modelling suggested that the overall structure of the protein would be largely unperturbed and the adduct would point outwards into the solvent. Mutating W165 to azF did slightly reduce overall catalytic efficiency towards ampicillin (1.5 fold) but crucially, modification with 3 was feasible and had little effect on overall activity (ESI Fig. 6†).
Fig. 4 AFM imaging of pyrene-modified TEMW165azF on graphene. (A) Model of TEMW165azF (magenta) modified with 1 aligned with the crystal structure of wt TEM (grey) with the residues either side of residue 165 and the catalytic S70 highlighted. The surface representation of TEMW165azF modified with 1 (left panel) illustrates that the Click adduct points out and away from the protein. Repeated AFM imaging of the surface after the (B) sixth and (C) tenth scans shows that the proteins bind stably to the surface. The height analysis of a single TEMW165azF-3 molecule, which is selected in C (highlighted as *) is shown as an inlay to C. ESI Movie 1† provides a full trace for the imaged area above over 10 scans. |
By coupling 3 to TEMW165azF prior to surface assembly we can directly interface the protein with graphene in a defined manner without surface coating with pyrene, as has been previously used.27,46 Intermittent-contact mode AFM (tapping mode) imaging indicated that TEMW165azF modified with pyrene binds stably to graphene surfaces (Fig. 4b and c; ESI Fig. 7 & Movie 1†). The protein molecules are not significantly disturbed by multiple scans. The average apparent heights were 3 nm, which is close to the predicted height of the protein bound in the designed orientation (∼3 nm; ESI Fig. 7†). The average apparent lateral dimension was larger than predicted (10 nm versus 5 nm), as a result of tip convolution effects. In comparison, wt TEM did not stably bind to the surface, with tip contamination commonly observed after multiple scans (ESI Fig. 8†).
Footnotes |
† Electronic supplementary information (ESI) available: Detailed experimental methods, supplementary Fig. 1 to 11 and supplementary Tables 1 to 3. See DOI: 10.1039/c4sc03900a |
‡ AMH, AJZ and ARM contributed equally to this work. |
This journal is © The Royal Society of Chemistry 2015 |