Paul J. Sapienza

Hidden dynamic allostery in a PDZ domain

Structure–function relationships in proteins are predicated on the spatial proximity of noncovalently interacting groups of atoms. Thus, structural elements located away from a proteins active site are typically presumed to serve a stabilizing or scaffolding role for the larger structure. Here we report a functional role for a distal structural element in a PDZ domain, even though it is not required to maintain PDZ structure. The third PDZ domain from PSD-95/SAP90 (PDZ3) has an unusual additional third alpha helix (α3) that packs in contiguous fashion against the globular domain. Although α3 lies outside the active site and does not make direct contact with C-terminal peptide ligand, removal of α3 reduces ligand affinity by 21-fold. Further investigation revealed that the difference in binding free energies between the full-length and truncated constructs is predominantly entropic in nature and that without α3, picosecond-nanosecond side-chain dynamics are enhanced throughout the domain, as determined by 2H methyl NMR relaxation. Thus, the distal modulation of binding function appears to occur via a delocalized conformational entropy mechanism. Without removal of α3 and characterization of side-chain dynamics, this dynamic allostery would have gone unnoticed. Moreover, what appeared at first to be an artificial modification of PDZ3 has been corroborated by experimentally verified phosphorylation of α3, revealing a tangible biological mechanism for this novel regulatory scheme. This hidden dynamic allostery raises the possibility of as-yet unidentified or untapped allosteric regulation in this PDZ domain and is a very clear example of function arising from dynamics rather than from structure.

Using NMR to study fast dynamics in proteins: methods and applications

Proteins exist not as singular structures with precise coordinates, but rather as fluctuating bodies that move rapidly through an enormous number of conformational substates. These dynamics have important implications for understanding protein function and for structure-based drug design. NMR spectroscopy is particularly well suited to characterize the dynamics of proteins and other molecules in solution at atomic resolution. Here, NMR relaxation methods for characterizing thermal motions on the picosecond-nanosecond (ps-ns) timescale are reviewed. Motion on this timescale can be conveniently captured by the Lipari-Szabo order parameter, S², a bond-specific measure of restriction of motion. Approaches for determining order parameters are discussed, as are recent examples from the literature that link ps-ns dynamics with conformational entropy, allostery, and protein function in general.

Multi-timescale dynamics study of FKBP12 along the rapamycin-mTOR binding coordinate.

Drugs can affect function in proteins by modulating their flexibility. Despite this possibility, there are very few studies on how drug binding affects the dynamics of target macromolecules. FKBP12 (FK506 binding protein 12) is a prolyl cis-trans isomerase and a drug target. The immunosuppressant drug rapamycin exerts its therapeutic effect by serving as an adaptor molecule between FKBP12 and the cell proliferation regulator mTOR (mammalian target of rapamycin). To understand the role of dynamics in rapamycin-based immunosuppression and to gain insight into the role of dynamics in the assembly of supramolecular complexes, we used (15)N, (13)C, and (2)H NMR spin relaxation to characterize FKBP12 along the binding coordinate that leads to cell cycle arrest. We show that sequential addition of rapamycin and mTOR leads to incremental rigidification of the FKBP12 backbone on the picosecond-nanosecond timescale. Both binding events lead to perturbation of main-chain and side-chain dynamics at sites distal to the binding interfaces, suggesting tight coupling interactions dispersed throughout the FKBP12-rapamycin interface. Binding of the first molecule, rapamycin, quenches microsecond-millisecond motions of the FKBP12 80s loop. This loop provides much of the surface buried at the protein-protein interface of the ternary complex, leading us to assert that preorganization upon rapamycin binding facilitates binding of the second molecule, mTOR. Widespread microsecond-millisecond motions of the backbone persist in the drug-bound enzyme, and we provide evidence that these slow motions represent coupled dynamics of the enzyme and isomerization of the bound drug. Finally, the pattern of microsecond-millisecond dynamics reported here in the rapamycin complex is dramatically different from the pattern in the complex with the structurally related drug FK506. This raises the important question of how two complexes that are highly isomorphic based on high-resolution static models have such different flexibilities in solution.

Structure and Dynamics of the G121V Dihydrofolate Reductase Mutant: Lessons from a Transition-State Inhibitor Complex

It is well known that enzyme flexibility is critical for function. This is due to the observation that the rates of intramolecular enzyme motions are often matched to the rates of intermolecular events such as substrate binding and product release. Beyond this role in progression through the reaction cycle, it has been suggested that enzyme dynamics may also promote the chemical step itself. Dihydrofolate reductase (DHFR) is a model enzyme for which dynamics have been proposed to aid in both substrate flux and catalysis. The G121V mutant of DHFR is a well studied form that exhibits a severe reduction in the rate of hydride transfer yet there remains dispute as to whether this defect is caused by altered structure, dynamics, or both. Here we address this by presenting an NMR study of the G121V mutant bound to reduced cofactor and the transition state inhibitor, methotrexate. NMR chemical shift markers demonstrate that this form predominantly adopts the closed conformation thereby allowing us to provide the first glimpse into the dynamics of a catalytically relevant complex. Based on 15N and 2H NMR spin relaxation, we find that the mutant complex has modest changes in ps-ns flexibility with most affected residues residing in the distal adenosine binding domain rather than the active site. Thus, aberrant ps-ns dynamics are likely not the main contributor to the decreased catalytic rate. The most dramatic effect of the mutation involves changes in µs-ms dynamics of the F-G and Met20 loops. Whereas loop motion is quenched in the wild type transition state inhibitor complex, the F-G and Met20 loops undergo excursions from the closed conformation in the mutant complex. These excursions serve to decrease the population of conformers having the correct active site configuration, thus providing an explanation for the G121V catalytic defect.

Mg2+ binds to the surface of thymidylate synthase and affects hydride transfer at the interior active site.

Thymidylate synthase (TSase) produces the sole intracellular de novo source of thymidine (i.e., the DNA base T) and thus is a common target for antibiotic and anticancer drugs. Mg(2+) has been reported to affect TSase activity, but the mechanism of this interaction has not been investigated. Here we show that Mg(2+) binds to the surface of Escherichia coli TSase and affects the kinetics of hydride transfer at the interior active site (16 Å away). Examination of the crystal structures identifies a Mg(2+) near the glutamyl moiety of the folate cofactor, providing the first structural evidence for Mg(2+) binding to TSase. The kinetics and NMR relaxation experiments suggest that the weak binding of Mg(2+) to the protein surface stabilizes the closed conformation of the ternary enzyme complex and reduces the entropy of activation on the hydride transfer step. Mg(2+) accelerates the hydride transfer by ~7-fold but does not affect the magnitude or temperature dependence of the intrinsic kinetic isotope effect. These results suggest that Mg(2+) facilitates the protein motions that bring the hydride donor and acceptor together, but it does not change the tunneling ready state of the hydride transfer. These findings highlight how variations in cellular Mg(2+) concentration can modulate enzyme activity through long-range interactions in the protein, rather than binding at the active site. The interaction of Mg(2+) with the glutamyl tail of the folate cofactor and nonconserved residues of bacterial TSase may assist in designing antifolates with polyglutamyl substitutes as species-specific antibiotic drugs.

Native State Volume Fluctuations in Proteins as a Mechanism for Dynamic Allostery

Allostery enables tight regulation of protein function in the cellular environment. Although existing models of allostery are firmly rooted in the current structure-function paradigm, the mechanistic basis for allostery in the absence of structural change remains unclear. In this study, we show that a typical globular protein is able to undergo significant changes in volume under native conditions while exhibiting no additional changes in protein structure. These native state volume fluctuations were found to correlate with changes in internal motions that were previously recognized as a source of allosteric entropy. This finding offers a novel mechanistic basis for allostery in the absence of canonical structural change. The unexpected observation that function can be derived from expanded, low density protein states has broad implications for our understanding of allostery and suggests that the general concept of the native state be expanded to allow for more variable physical dimensions with looser packing.

The Effect of Protein Mass Modulation on Human Dihydrofolate Reductase.

Dihydrofolate reductase (DHFR) from Escherichia coli has long served as a model enzyme with which to elucidate possible links between protein dynamics and the catalyzed reaction. Such physical properties of its human counterpart have not been rigorously studied so far, but recent computer-based simulations suggest that these two DHFRs differ significantly in how closely coupled the protein dynamics and the catalyzed C-H → C hydride transfer step are. To test this prediction, two contemporary probes for studying the effect of protein dynamics on catalysis were combined here: temperature dependence of intrinsic kinetic isotope effects (KIEs), which are sensitive to the physical nature of the chemical step, and protein mass modulation, which slows down fast dynamics (femto- to picosecond time scale) throughout the protein. The intrinsic H/T KIEs of human DHFR, like those of E. coli DHFR, are shown to be temperature-independent in the range from 5 to 45 °C, indicating fast sampling of donor and acceptor distances (DADs) at the reactions transition state (or tunneling ready state, TRS). Mass modulation of these enzymes through isotopic labeling with (13)C, (15)N, and (2)H at nonexchangeable hydrogens yields an 11% heavier enzyme. The additional mass has no effect on the intrinsic KIEs of the human enzyme. This finding indicates that the mass modulation of the human DHFR affects neither DAD distribution nor the DADs conformational sampling dynamics. Furthermore, reduction in the enzymatic turnover number and the dissociation rate constant for the product indicate that the isotopic substitution affects kinetic steps that are not the catalyzed C-H → C hydride transfer. The findings are discussed in terms of fast dynamics and their role in catalysis, the comparison of calculations and experiments, and the interpretation of isotopically modulated heavy enzymes in general.

Chemical shift imprint of intersubunit communication in a symmetric homodimer.

Significance Regulatory signals between protein subunits depend on communication between sequential binding events. How such long-range communication, or allostery, operates at the subdomain level has been elusive, especially for homooligomeric proteins. To address this problem, homodimers of thymidylate synthase were generated for optimized study of individual protomers of the singly bound dimer. Mixed 15N-labeled dimers were created with a single functional active site, allowing site-specific, protomer-specific chemical shifts to report on step-wise binding effects. Long-range intersubunit communication was observed although this communication was apparent only in the second ligand-binding step, in which changes were in the first ligand-bound region. Visualization of up to four peaks for each residue amide provides a unique way to assess the allosteric mechanism. Allosteric communication is critical for protein function and cellular homeostasis, and it can be exploited as a strategy for drug design. However, unlike many protein–ligand interactions, the structural basis for the long-range communication that underlies allostery is not well understood. This lack of understanding is most evident in the case of classical allostery, in which a binding event in one protomer is sensed by a second symmetric protomer. A primary reason why study of interdomain signaling is challenging in oligomeric proteins is the difficulty in characterizing intermediate, singly bound species. Here, we use an NMR approach to isolate and characterize a singly ligated state (“lig1”) of a homodimeric enzyme that is otherwise obscured by rapid exchange with apo and saturated forms. Mixed labeled dimers were prepared that simultaneously permit full population of the lig1 state and isotopic labeling of either protomer. Direct visualization of peaks from lig1 yielded site-specific ligand-state multiplets that provide a convenient format for assessing mechanisms of intersubunit communication from a variety of NMR measurements. We demonstrate this approach on thymidylate synthase from Escherichia coli, a homodimeric enzyme known to be half-the-sites reactive. Resolving the dUMP1 state shows that active site communication occurs not upon the first dUMP binding, but upon the second. Surprisingly, for many sites, dUMP1 peaks are found beyond the limits set by apo and dUMP2 peaks, indicating that binding the first dUMP pushes the enzyme ensemble to further conformational extremes than the apo or saturated forms. The approach used here should be generally applicable to homodimers.

Thermodynamic and NMR assessment of ligand cooperativity and intersubunit communication in symmetric dimers: application to thymidylate synthase

Thymidylate synthase (TS) is a homodimeric enzyme with evidence for negative regulation of one protomer while the other protomer acts on substrate, so called half-the-sites reactivity. The mechanisms by which multisubunit allosteric proteins communicate between protomers is not well understood, and the simplicity of dimeric systems has advantages for observing conformational and dynamic processes that functionally connect distance-separated active sites. This review considers progress in overcoming the inherent challenges of accurate thermodynamic and atomic-resolution characterization of interprotomer communication mechanisms in symmetric protein dimers, with TS used as an example. Isothermal titration calorimetry (ITC) is used to measure ligand binding cooperativity, even in cases where the two binding enthalpies are similar, and NMR spectroscopy is used to detect site-specific changes occurring in the two protomers. The NMR approach makes use of mixed-labeled dimers, enabling protomer-specific detection of signals in the singly ligated state. The rich informational content of the NMR signals from the singly ligated state, relative to the apo and saturated states, requires new considerations that do not arise in simple cases of 1:1 protein-ligand interactions.