Report

Trapping a transition state in a computationally designed protein bottle

See allHide authors and affiliations

Science  20 Feb 2015:
Vol. 347, Issue 6224, pp. 863-867
DOI: 10.1126/science.aaa2424

A transition state holds a pose

The transition state of a chemical transformation is inherently fleeting because the structure is high in energy. Nonetheless, Pearson et al. trapped a classical example of a bond rotation transition state using a modified protein (see the Perspective by Romney and Miller). The biphenyl molecule passes through an energy maximum when its rings rotate through a parallel position. A pocket within the editing domain of threonyl–transfer RNA synthetase was modified to stabilize parallel biphenyl rings, allowing further characterization of this normally transient structure.

Science, this issue p. 863; see also p. 829

Abstract

The fleeting lifetimes of the transition states (TSs) of chemical reactions make determination of their three-dimensional structures by diffraction methods a challenge. Here, we used packing interactions within the core of a protein to stabilize the planar TS conformation for rotation around the central carbon-carbon bond of biphenyl so that it could be directly observed by x-ray crystallography. The computational protein design software Rosetta was used to design a pocket within threonyl-transfer RNA synthetase from the thermophile Pyrococcus abyssi that forms complementary van der Waals interactions with a planar biphenyl. This latter moiety was introduced biosynthetically as the side chain of the noncanonical amino acid p-biphenylalanine. Through iterative rounds of computational design and structural analysis, we identified a protein in which the side chain of p-biphenylalanine is trapped in the energetically disfavored, coplanar conformation of the TS of the bond rotation reaction.

The direct observation of the transition state (TS) of a chemical reaction requires highly sensitive spectroscopic techniques with a temporal resolution on the order of 10−13 to 10−14 s. This milestone was made possible with the advent of ultrafast laser spectroscopy, which allowed the direct observation of the transient species formed as reactants cross an energy barrier to products (1, 2). However, determination of the three-dimensional (3D) structure of a TS requires ultrafast electron or x-ray diffraction techniques with similar temporal resolution and sensitivity (3, 4). An alternative approach is to increase the lifetime of the transition state configuration by trapping it in a free energy (ΔG) well so that it can be directly observed by more conventional spectroscopic methods. In theory, the well-packed interior of a protein could act as a “programmable” solvent to selectively stabilize a TS configuration and thereby make it kinetically persistent and directly observable by conventional x-ray diffraction methods. Indeed, enzymes catalyze reactions by virtue of their ability to selectively stabilize the rate-limiting TS relative to the ground-state reactants (59). Here, we show the redesign of the interior of the protein threonyl-tRNA synthetase from the thermophile Pyrococcus abyssi (10) to create a microenvironment that has high complementarity to the planar TS configuration for rotation about the central C–C bond of biphenyl (11). Determination of the 3D structure of the designed protein at 2.05 Å resolution by x-ray crystallography revealed a planar biphenyl TS configuration stabilized by van der Waals interactions with side chains in the protein core.

Biphenyl bond rotation is a well-studied reaction, both experimentally and theoretically. In the gas phase, the ground state of biphenyl is twisted, with a dihedral angle of ~45°. Rotation around the central C­–C bond connecting the two phenyl rings in biphenyl is estimated to have an energy barrier of 5.8 ± 2.1 and 6.7 ± 2.1 kJ/mol around the 0° (planar) and 90° TSs, respectively, as determined experimentally by electron diffraction studies (Fig. 1A) (12). The energy barriers estimated using Raman data (13, 14) and ultraviolet absorption spectroscopy (15) are in general agreement with the results from the diffraction studies. Theoretical calculations show that these barriers result from opposing steric and electronic effects of the two phenyl rings but give a range of values depending on the method used (16). This simple reaction (which racemizes chiral biphenyls) provides an ideal system to pose the question whether side chain packing interactions in a folded protein can be exploited to make a transition state configuration kinetically persistent so that it can be directly observed by x-ray crystallography.

Fig. 1 Strategy to capture a transition state.

(A) The orientations and relative changes in energy (ΔGrel) between the ground (Φ = ~45°) and transition states (Φ = 0° and 90°) of rotation around the central bond of biphenyl are shown. (B) An overview of the computational design process used to generate candidate designs is shown.

To introduce a biphenyl moiety into a protein core, we used the genetically encoded noncanonical amino acid p-biphenylalanine (BiPhe). This amino acid can be introduced site-specifically into a protein in Escherichia coli in good yield in response to the amber nonsense codon TAG, with an orthogonal amber suppressor tRNA/aminoacyl-tRNA synthetase (aaRS) pair. The aaRS was evolved from the Methanococcus jannaschi tyrosyl-tRNA synthetase to be selective for BiPhe and not incorporate any of the 20 canonical amino acids (17, 18). The rotational barrier around the biphenyl C–C bond in BiPhe is not expected to be substantially different from biphenyl itself (19).

To identify an appropriate host protein in which to construct a site complementary to the planar biphenyl TS geometry, we first curated a set of ~2300 proteins of known structures from thermophilic organisms [a full description of the computational design process is found in the supplementary materials (20)]. Proteins with high thermostabilities have more negative ΔGs of folding and are in general more tolerant to mutations; both characteristics are likely to be beneficial when attempting to stabilize a TS configuration through packing interactions in a protein core (21, 22). RosettaMatch (23) was then used to identify residues in the core of each protein scaffold where BiPhe could be substituted without creating unfavorable steric interactions with the protein backbone (Fig. 1B), and make π-stacking interactions (24) with native Trp, Phe, or Tyr residues. Because the RosettaMatch calculations were carried out with a planar model of the BiPhe, these interactions should stabilize the biphenyl side chain in the desired planar TS conformation. To ensure that the substituted BiPhe side chain was sufficiently buried within a particular protein’s core, initial hits were first filtered on the basis of the change in solvent-accessible surface area (ΔSASA) (25) that occurs when BiPhe is removed from the substituted site. A ΔSASA cutoff of 0.9 (meaning the BiPhe was >90% buried) removed 35% of the initial matches from further consideration.

RosettaDesign (26) was then used to optimize the identities of residues surrounding the BiPhe (excluding residues already participating in π-stacking interactions) such that they packed tightly against the planar BiPhe side chain but did not create unfavorable hydrogen-bonding interactions (Fig. 1B). Candidate designs were ranked by shape complementarity (SC) (27) between the designed residues and the BiPhe (Fig. 1B). Ultimately, four designs (BIF_1 to BIF_4) with high SC values were chosen for experimental characterization (Table 1). The computer-modeled proteins containing the BiPhe were reverse-translated and optimized for expression in E. coli.

Table 1 First-round computational designs.

Designed protein names, parent scaffolds, and mutations made are listed. Single-letter abbreviations for the amino acid residues are as follows: A, Ala; C, Cys; D, Asp; E, Glu; F, Phe; G, Gly; H, His; I, Ile; K, Lys; L, Leu; M, Met; N, Asn; P, Pro; Q, Gln; R, Arg; S, Ser; T, Thr; V, Val; W, Trp; and Y, Tyr.

View this table:

To incorporate BiPhe into these four different protein scaffolds, we used the orthogonal amber suppressor tRNA/BiPheRS pair encoded on the dual-plasmid expression system (pUltra) (28). One plasmid contained the tRNA/BiPheRS pair specific for BiPhe, and the other contained the synthetic mutant gene of interest fused to a C-terminal hexa-histidine purification tag. A TAG amber nonsense codon was introduced at the desired site to encode BiPhe (17). These plasmids were cotransformed into E. coli Bl21(DE3), and protein expression was carried out in the presence of 1 mM BiPhe. Proteins were purified from cell lysate via Ni2+-affinity chromatography followed by size-exclusion chromatography. Mass spectrometric analysis of BIF_1 to BIF_4 indicated successful incorporation of BiPhe in all cases (fig. S1), and SDS–polyacrylamide gel electrophoresis indicated that the purified proteins were of suitable purity (>95%) for crystallographic analysis.

All four designed proteins were subjected to an initial crystallographic screen based on the conditions used to crystallize the wild-type protein (20, 2932). Crystals were obtained only for BIF_1, which has as its parent scaffold threonyl-tRNA synthetase from the thermophile P. abyssi, but they were needles and not suitable for x-ray crystallography. An automated high-throughput crystallization system was then used to identify new conditions (0.1 M sodium citrate, 15% polyethylene glycol 6000, pH = 5.5) in which large crystals of BIF_1 grew. The structure of BIF_1 was solved to 1.8 Å resolution by x-ray crystallography using molecular replacement with the parent scaffold (PDB ID 1y2q) serving as a search model. Density for BiPhe was clearly observed in a 2FoFc map (Fig. 2A), and the torsion angle between the two phenyl rings was determined to be ~28° (Fig. 2A). This value represented a rotation of ~17° toward the desired planar conformation relative to the dihedral found at the energetic minimum but remained far from the desired value of 0°. To force the BiPhe torsion angle closer to planarity, a second round of computational design was undertaken on the basis of a detailed analysis of the structure of BIF_1.

Fig. 2 X-ray crystallographic analysis of BIF_1.

(A) The x-ray crystal structure of BIF_1 is shown in yellow; BiPhe and surrounding residues are shown in sticks. Rings A and B are those closest to and farthest from the protein backbone, respectively. Electron density around the BiPhe side chain is shown as a 2FoFc map contoured to 2σ. (B) A comparison of the design model (gray) to the structure (yellow) of BIF_1 is shown; BiPhe, Trp42, and Trp81 are shown in sticks. A loop corresponding to residues 81 to 89 of the parent scaffold is shown in red. Missing density in the structure corresponding to residues 83 to 86 of BIF_1 is shown as a dashed red line. (C) A comparison of the structure of BIF_1 (yellow) to the design model (gray) is shown. BiPhe, Trp42, and Trp81 are shown in sticks. An arrow indicates rotation about χ2 in the structure relative to the design.

Globally, the structure of BIF_1 matched the design model quite well (Fig. 2B; root mean square deviation to the design model of ~1.3 Å over all atoms), although differences are apparent in the vicinity of the BiPhe residue. In the model, Trp42 and Trp81 both form edgewise interactions with the BiPhe side chain. However, in the crystal structure, the indole ring of Trp42 rotates such that it packs lengthwise against the BiPhe side chain in an orientation that would clash with the side-chain orientation of Trp81 predicted by the design model (Fig. 2C). To avoid this unfavorable steric interaction, a substantial displacement of the loop consisting of residues 81 to 89 likely occurs, as evidenced by the lack of electron density for residues 83 to 86 in the crystal structure (Fig. 2B). In an attempt to return the disordered loop to its native position, we independently mutated each Trp to Phe (the wild-type residue at both positions) in silico. Unconstrained repacking and minimization calculations in the context of each mutation showed that Trp42Phe increased SC and scored slightly better than the original design. As a result, this mutation was made standard for the remainder of the computational redesign.

A second focus of the redesign effort was to identify point mutations in residues packing against the biphenyl rings to force the side chain into the desired planar conformation. In the structure of BIF_1, the biphenyl side chain is rotated ~5° about χ1 and 25° about χ2, relative to its placement in the design model (Fig. 2C). Although the phenyl ring closest to the backbone (ring A) is out of plane with respect to the model, the ring farthest from the backbone (ring B) is essentially in the plane of the design model (Fig. 2, A and C). Thus, it appears that rotation about χ2 is the predominant determinant of the 28° deviation from planarity observed in the crystal structure relative to the design model. Because ring A is bounded on one face by the protein backbone, we believed it would be difficult to identify a mutant that would adjust this ring in the desired direction (Fig. 2C). In contrast, ring B is flanked by both Ala79 and Tyr123, suggesting that mutagenesis of one or both of these residues could potentially planarize the BiPhe dihedral angle (Fig. 2A). Analysis of the structure of BIF_1 suggested that the phenyl ring of Tyr123 likely prevents ring B from rotating into the plane of ring A (Fig. 2A); thus, we mutated Tyr123 in silico to the smaller residues Ala and Val. Concurrently, Ala79 was mutated to the bulkier residues Cys, Ser, Thr, and Val. Rosetta was then used to analyze these potential sequence alternatives, again by carrying out unconstrained repacking and minimization calculations. The energies for all mutants tested fell within ~6 Rosetta energy units (REU) of one another and gave SC values that differed at most by ~5%. The tight distribution of values for both metrics suggested that no clear preference exists for one mutant over another, so a series of four mutants was examined experimentally.

All four mutants expressed well and afforded diffraction-quality crystals. We solved the structures of BIF_1.1, BIF_1.2, BIF_ 1.3, and BIF_1.4 to 2.10, 2.50, 2.36, and 2.10 Å resolution, respectively, again by molecular replacement. In all cases, the Trp42Phe mutation returned the displaced loop to its native position (fig. S2). However, a distribution of BiPhe torsion angles between 35° and 15° was observed among the four structures (Table 2). In the majority of these cases, the change in the χ2 angle relative to the original design remained near the value of 25° observed in the initial BIF_1 structure. This result suggested that mutations to Ala79 and Trp123 have the desired effect of rotating ring B without affecting the absolute orientation of ring A. Unfortunately, substitution of Tyr123 with the beta-branched Val and the opposing Ala79 with Ser (BIF_1.1) had the effect of rotating ring B even farther out of plane (35°) than in the original structure (Fig. 3A). This undesired rotation was partially remedied in BIF_1.2 (Φ = 21°) (Fig. 3B) by substituting Ala79 with a bulkier Val residue and further corrected in BIF_1.3 and BIF_1.4 (Φ = 15° and 20°, respectively) (Fig. 3, C and D) by replacing the opposing Tyr123 with a smaller Ala residue and Ala79 with Ser (BIF_1.3) or Val (BIF_1.4). This analysis suggested that ring A could potentially be rotated into a coplanar geometry by further increasing the size of the amino acid at position 79 with an Ala79Ile mutation while maintaining Phe42 and the Tyr123Ala mutation. The additional methyl group of the isoleucine should force the side of ring A to rotate further in the desired direction.

Table 2 Second- and third-round crystallographic analysis.

Second-round mutant identities, biphenylalanine dihedral angle, and x-ray crystal structure resolution are listed. Dihedrals listed are averages of those measured on each side of the biphenyl ring.

View this table:
Fig. 3 A comparison of the crystal structures of BIF_1.1 to BIF_1.4.

(A to D) Crystal structures of second-round mutants BIF_1.1 to BIF_1.4 are shown. The side chains of BiPhe, and those at positions 79 and 123, are shown in sticks. Electron density from a 2FoFc map contoured to 1.5σ [(A to C)] and 2.0σ (D) is shown for the aforementioned residues. The measured dihedral angle between the two biphenyl rings is shown beneath the biphenyl side chain in each case.

We next generated the corresponding BIF_0 mutant (S8A, I11BiPhe, Y79I, F81W, K121I, and F123A), purified the protein, and solved its crystal structure to 2.05 Å resolution (Fig. 4). Analysis of the electron density showed that the two phenyl rings of BiPhe are coplanar, which matches the configuration of the TS for the bond rotation reaction. The structure of BIF_0 shows that, in addition to adding steric bulk beneath ring A, the V79I mutation also forces the side chain of Phe77 to adopt a different rotamer than was observed in BIF_1.4, which has the effect of further rotating ring A into the plane of ring B (Fig. 4). The mutations introduced into BIF_0 do not appear to substantially affect the thermal stability of the protein. The melting temperature of this mutant, as determined by differential scanning calorimetry, was ~110°C, consistent with the 3D structure of BIF_0, which shows that the protein core is well packed.

Fig. 4 X-ray crystal structure of BIF_0.

(A) The crystal structure of BIF_0 is shown in blue, and BIF_1.3 is shown in gray. The BiPhe side chain and surrounding residues are shown in sticks. (B) Packing interactions between the designed protein BIF_0 and the BiPhe side are highlighted with space-filling representations of the interacting residues. (C and D) The structure of BIF_0 is shown, highlighting the BiPhe side chain; views from the front and side are shown. The BiPhe side chain and surrounding residues are shown in sticks, and 2FoFc maps are contoured to 2σ in each case.

We have shown by iterative computational design, mutagenesis, and protein structure determination that one can design a protein core that stabilizes a simple conformational transition state to such a degree that one can determine its 3D x-ray crystal structure. However, we should note that the biphenyl energy landscape corresponds to a substructure within the protein relative to the energetics of the global protein conformational ensemble. A similar strategy was recently employed to directly observe catalyst-substrate interactions through x-ray crystallographic analysis (33). The results described here may not be all that surprising given that enzymes typically stabilize a rate-limiting TS by 8 to 12 kcal/mol. Nonetheless, these experiments underscore the ability of proteins to fold into defined 3D structures in which van der Waals, hydrogen-bonding, and electrostatic interactions can be controlled with exquisite precision.

Supplementary Materials

www.sciencemag.org/content/347/6224/863/suppl/DC1

Materials and Methods

Figs. S1 and S2

Tables S1 and S2

References (3442)

References and Notes

  1. Materials and methods are available as supplementary materials on Science Online.
  2. Acknowledgments: The authors thank N. P. King and P.-S. Huang for helpful discussions. D.B and J.H.M. were supported by the Defense Threat Reduction Agency (HDTRA1-11-1-0041). J.H.M. was supported by National Institute of General Medical Science of the National Institutes of Health under award F32GM099210. P.G.S. acknowledges support by the National Institutes of Health under award 2 R01 GM097206-05. The content is solely the responsibility of the authors and does not represent the official views of the National Institutes of Health. Structures of BIF_1, BIF_1.1 to BIF_1.4, and BIF_0 have been deposited in the Protein Data Bank under accession numbers 4S02, 4S0J, 4S0L, 4S0I, 4S0K, and 4S03.
View Abstract

Navigate This Article