Carbanionic intermediates play a central role in the catalytic transformations of amino acids performed by pyridoxal-5′-phosphate (PLP)-dependent enzymes. Here, we make use of NMR crystallography—the synergistic combination of solid-state nuclear magnetic resonance, X-ray crystallography, and computational chemistry—to interrogate a carbanionic/quinonoid intermediate analogue in the β-subunit active site of the PLP-requiring enzyme tryptophan synthase. The solid-state NMR chemical shifts of the PLP pyridine ring nitrogen and additional sites, coupled with first-principles computational models, allow a detailed model of protonation states for ionizable groups on the cofactor, substrates, and nearby catalytic residues to be established. Most significantly, we find that a deprotonated pyridine nitrogen on PLP precludes formation of a true quinonoid species and that there is an equilibrium between the phenolic and protonated Schiff base tautomeric forms of this intermediate. Natural bond orbital analysis indicates that the latter builds up negative charge at the substrate Cα and positive charge at C4′ of the cofactor, consistent with its role as the catalytic tautomer. These findings support the hypothesis that the specificity for β-elimination/replacement versus transamination is dictated in part by the protonation states of ionizable groups on PLP and the reacting substrates and underscore the essential role that NMR crystallography can play in characterizing both chemical structure and dynamics within functioning enzyme active sites.

The acid–base chemistry that drives catalysis in pyridoxal-5′-phosphate (PLP)-dependent enzymes has been the subject of intense interest and investigation since the initial identification of PLP’s role as a coenzyme in this extensive class of enzymes. It was first proposed over 50 years ago that the initial step in the catalytic cycle is facilitated by a protonated Schiff base form of the holoenzyme in which the linking lysine ε-imine nitrogen, which covalently binds the coenzyme, is protonated. Here we provide the first 15N NMR chemical shift measurements of such a Schiff base linkage in the resting holoenzyme form, the internal aldimine state of tryptophan synthase. Double-resonance experiments confirm the assignment of the Schiff base nitrogen, and additional 13C, 15N, and 31P chemical shift measurements of sites on the PLP coenzyme allow a detailed model of coenzyme protonation states to be established.

Ab initio chemical shielding calculations greatly facilitate the interpretation of nuclear magnetic resonance (NMR) chemical shifts in biological systems, but the large sizes of these systems requires approximations in the chemical models used to represent them. Achieving good convergence in the predicted chemical shieldings is necessary before one can unravel how other complex structural and dynamical factors affect the NMR measurements. Here, we investigate how to balance trade-offs between using a better basis set or a larger cluster model for predicting the chemical shieldings of the substrates in two representative examples of protein-substrate systems involving different domains in tryptophan synthase: the N-(4′-trifluoromethoxybenzoyl)-2-aminoethyl phosphate (F9) ligand which binds in the α active site, and the 2-aminophenol (2AP) quinonoid intermediate formed in the β active site. We first demonstrate that a chemically intuitive three-layer, locally dense basis model that uses a large basis on the substrate, a medium triple-zeta basis to describe its hydrogen-bonding partners and/or surrounding van derWaals cavity, and a crude basis set for more distant atoms provides chemical shieldings in good agreement with much more expensive large basis calculations. Second, long-range quantum mechanical interactions are important, and one can accurately estimate them as a small-basis correction to larger-basis calculations on a smaller cluster. The combination of these approaches enables one to perform density functional theory NMR chemical shift calculations in protein systems that are well-converged with respect to both basis set and cluster size.

Intrinsically disordered protein domains are not only found in soluble proteins but can also be part of large protein complexes or protein aggregates. For example, several amyloid fibrils have intrinsically disordered domains framing a rigid β-sheet rich core. These disordered domains can often be observed using solution NMR methods in combination with modest magic angle spinning and without perdeuteration. But how can these regions be detected using solution NMR methods when they are part of a fibril that is not tumbling isotropically in solution? Here we addressed this question by investigating the dynamic C-terminus of huntingtin exon-1 (HTTex1) fibrils that are important in Huntington’s Disease. We assigned the most dynamic regions of the C-terminus of three HTTex1 variants. Based on this assignment, we measured site-specific residual chemical shifts, peak intensities, and R1, R′2 and R1ρ 15N relaxation rates. In addition, we determined the residual 1H-15N dipolar couplings of this region. Our results show that the dipolar couplings are averaged to a very high degree resulting in an order parameter that is essentially zero. Together, our data show that the C-terminus of HTTex1 is intrinsically disordered and undergoes motions in the high picosecond to low nanosecond range.

Four new X-ray structures of tryptophan synthase (TS) crystallized with varying numbers of the amphipathic N-(4′-trifluoromethoxybenzoyl)-2-aminoethyl phosphate (F6) molecule are presented. These structures show one of the F6 ligands threaded into the tunnel from the β-site and reveal a distinct hydrophobic region. Over this expanse, the interactions between F6 and the tunnel are primarily nonpolar, while the F6 phosphoryl group fits into a polar pocket of the β-subunit active site. Further examination of TS structures reveals that one portion of the tunnel (T1) binds clusters of water molecules, whereas waters are not observed in the nonpolar F6 binding region of the tunnel (T2). MD simulation of another TS structure with an unobstructed tunnel also indicates the T2 region of the tunnel excludes water, consistent with a dewetted state that presents a significant barrier to the transfer of water into the closed β-site. We conclude that hydrophobic molecules can freely diffuse between the α- and β-sites via the tunnel, while water does not. We propose that exclusion of water serves to inhibit reaction of water with the α-aminoacrylate intermediate to form ammonium ion and pyruvate, a deleterious side reaction in the αβ-catalytic cycle. Finally, while most TS structures show βPhe280 partially blocking the tunnel between the α- and β-sites, new structures show an open tunnel, suggesting the flexibility of the βPhe280 side chain. Flexible docking studies and MD simulations confirm the dynamic behavior of βPhe280 allows unhindered transfer of indole through the tunnel, therefore excluding a gating role for this residue.