Decoupling of catalysis and transition state analog binding from mutations throughout a phosphatase revealed by high-throughput enzymology.

Using high-throughput microfluidic enzyme kinetics (HT-MEK), we measured over 9,000 inhibition curves detailing impacts of 1,004 single-site mutations throughout the alkaline phosphatase PafA on binding affinity for two transition state analogs (TSAs), vanadate and tungstate. As predicted by catalytic models invoking transition state complementary, mutations to active site and active-site-contacting residues had highly similar impacts on catalysis and TSA binding. Unexpectedly, most mutations to more distal residues that reduced catalysis had little or no impact on TSA binding and many even increased tungstate affinity. These disparate effects can be accounted for by a model in which distal mutations alter the enzyme's conformational landscape, increasing the occupancy of microstates that are catalytically less effective but better able to accommodate larger transition state analogs. In support of this ensemble model, glycine substitutions (rather than valine) were more likely to increase tungstate affinity (but not more likely to impact catalysis), presumably due to increased conformational flexibility that allows previously disfavored microstates to increase in occupancy. These results indicate that residues throughout an enzyme provide specificity for the transition state and discriminate against analogs that are larger only by tenths of an Ångström. Thus, engineering enzymes that rival the most powerful natural enzymes will likely require consideration of distal residues that shape the enzyme's conformational landscape and fine-tune active-site residues. Biologically, the evolution of extensive communication between the active site and remote residues to aid catalysis may have provided the foundation for allostery to make it a highly evolvable trait.

High-Throughput Affinity Measurements of Transcription Factor and DNA Mutations Reveal Affinity and Specificity Determinants.

Transcription factors (TFs) bind regulatory DNA to control gene expression, and mutations to either TFs or DNA can alter binding affinities to rewire regulatory networks and drive phenotypic variation. While studies have profiled energetic effects of DNA mutations extensively, we lack similar information for TF variants. Here, we present STAMMP (simultaneous transcription factor affinity measurements via microfluidic protein arrays), a high-throughput microfluidic platform enabling quantitative characterization of hundreds of TF variants simultaneously. Measured affinities for ∼210 mutants of a model yeast TF (Pho4) interacting with 9 oligonucleotides (>1,800 Kds) reveal that many combinations of mutations to poorly conserved TF residues and nucleotides flanking the core binding site alter but preserve physiological binding, providing a mechanism by which combinations of mutations in cis and trans could modulate TF binding to tune occupancies during evolution. Moreover, biochemical double-mutant cycles across the TF-DNA interface reveal molecular mechanisms driving recognition, linking sequence to function. A record of this paper's Transparent Peer Review process is included in the Supplemental Information.

Role of Polyubiquitin Chain Flexibility in the Catalytic Mechanism of Cullin-RING Ubiquitin Ligases.

Cullin-RING ubiquitin ligases are a diverse family of ubiquitin ligases that catalyze the synthesis of K48-linked polyubiquitin (polyUb) chains on a variety of substrates, ultimately leading to their degradation by the proteasome. The cullin-RING enzyme scaffold processively attaches a Ub molecule to the distal end of a growing chain up to lengths of eight Ub monomers. However, the molecular mechanism governing how chains of increasing size are built using a scaffold of largely fixed dimensions is not clear. We developed coarse-grained molecular dynamics simulations to describe the dependence of kcat for cullin-RING ligases on the length and flexibility of the K48-linked polyUb chain attached to the substrate protein, key factors that determine the rate of subsequent Ub attachment to the chain, and therefore, the ensuing biological outcomes of ubiquitination. The results suggest that a number of regulatory mechanisms may lead to variations in the rate of chain elongation for different cullin-RING ligases. Specifically, modulation of the distance between the target lysine and the phosphodegron sequence of the substrate, the distance between the substrate lysine and the active site cysteine of the Ub conjugation enzyme (E2) bound to the cullin-RING scaffold, and flexibility of the bound E2 can lead to significant differences in the processing of K48-linked chains on substrates, potentially leading to differences in biological outcomes.

An Open-Source, Programmable Pneumatic Setup for Operation and Automated Control of Single- and Multi-Layer Microfluidic Devices.

Microfluidic technologies have been used across diverse disciplines (e.g. high-throughput biological measurement, fluid physics, laboratory fluid manipulation) but widespread adoption has been limited in part due to the lack of openly disseminated resources that enable non-specialist labs to make and operate their own devices. Here, we report the open-source build of a pneumatic setup capable of operating both single and multilayer (Quake-style) microfluidic devices with programmable scripting automation. This setup can operate both simple and complex devices with 48 device valve control inputs and 18 sample inputs, with modular design for easy expansion, at a fraction of the cost of similar commercial solutions. We present a detailed step-by-step guide to building the pneumatic instrumentation, as well as instructions for custom device operation using our software, Geppetto, through an easy-to-use GUI for live on-chip valve actuation and a scripting system for experiment automation. We show robust valve actuation with near real-time software feedback and demonstrate use of the setup for high-throughput biochemical measurements on-chip. This open-source setup will enable specialists and novices alike to run microfluidic devices easily in their own laboratories.

Covalent Inhibition of Ubc13 Affects Ubiquitin Signaling and Reveals Active Site Elements Important for Targeting.

Ubc13 is an E2 ubiquitin conjugating enzyme that functions in nuclear DNA damage signaling and cytoplasmic NF-κB signaling. Here, we present the structures of complexes of Ubc13 with two inhibitors, NSC697923 and BAY 11-7082, which inhibit DNA damage and NF-κB signaling in human cells. NSC697923 and BAY 11-7082 both inhibit Ubc13 by covalent adduct formation through a Michael addition at the Ubc13 active site cysteine. The resulting adducts of both compounds exploit a binding groove unique to Ubc13. We developed a Ubc13 mutant which resists NSC697923 inhibition and, using this mutant, we show that the inhibition of cellular DNA damage and NF-κB signaling by NSC697923 is largely due to specific Ubc13 inhibition. We propose that unique structural features near the Ubc13 active site could provide a basis for the rational development and design of specific Ubc13 inhibitors.

Stochastic gate dynamics regulate the catalytic activity of ubiquitination enzymes.

Initiation of the DNA damage and innate immune responses is dependent upon the flow of chemical information through coupled protein-protein interaction networks and driven by the synthesis and recognition of Lys 63 linked polyubiquitin (polyUb) chains on adaptor proteins. The central chemical step in Lys 63-linked protein ubiquitination involves the reaction of a specific lysine on a target protein with Ub that is covalently attached as a thioester conjugate to the Ub conjugating enzyme (E2) Ubc13. The active site cysteine of Ubc13, and E2 enzymes in general, is buttressed by a flexible loop. The role of loop dynamics in catalysis was investigated by mutating the central and hinge residues to glycine. The loop dynamics were experimentally characterized through measurement of enzyme kinetics, main chain NMR relaxation, X-ray crystallographic studies, and in vivo studies in yeast. The experimental data were complemented by analysis of MD simulations of the dynamics and kinetics for the loop motion. The results show that fast pico- to nanosecond time scale active site loop fluctuations play a crucial role in regulating the catalytic activity of Ubc13 by functioning as a stochastic active site gate, which is characterized by precisely balanced rates of opening and closing. In vivo functional complementation assays in yeast demonstrate that defects within this regulatory mechanism can have profound biological consequences, given that Ubc13 is the only E2 dedicated to synthesizing Lys 63-linked polyUb chains.

Molecular basis for impaired DNA damage response function associated with the RAP80 ΔE81 defect.

Signal transduction within the DNA damage response is driven by the flux of protein-protein interaction cascades that ultimately recruit repair complexes to sites of damage. The protein RAP80 plays a central role in the damage response by targeting BRCA1/BRCA2 tumor suppressors to DNA damage foci through multivalent binding of Lys-63-linked polyubiquitin chains. Mutations within the high penetrance BRCA1/BRCA2 genes account for ∼20% of familial breast cancers. The genetic basis for the remaining cancers remains unknown, but may involve defects in binding partners for BRCA1 and BRCA2 that lead to impaired targeting to foci and a concomitant role in the pathogenesis of cancer. Recently, an in-frame deletion mutation (ΔE81) in a conserved region from the first ubiquitin interaction motif of RAP80 has been linked to an increase in chromosomal abnormalities. Using NMR spectroscopy, we demonstrate that the N-cap motif within the α-helix of the first ubiquitin interaction motif from ΔE81 undergoes a structural frameshift that leads to abolishment of multivalent binding of polyubiquitin chains. Loss of this single glutamate residue disrupts favorable electrostatic interactions between RAP80 and ubiquitin, establishing a plausible molecular basis for a potential predisposition to cancer unrelated to mutations within BRCA1/BRCA2 genes.

Accuracy and precision of protein-ligand interaction kinetics determined from chemical shift titrations.

NMR-monitored chemical shift titrations for the study of weak protein-ligand interactions represent a rich source of information regarding thermodynamic parameters such as dissociation constants (K ( D )) in the micro- to millimolar range, populations for the free and ligand-bound states, and the kinetics of interconversion between states, which are typically within the fast exchange regime on the NMR timescale. We recently developed two chemical shift titration methods wherein co-variation of the total protein and ligand concentrations gives increased precision for the K ( D ) value of a 1:1 protein-ligand interaction (Markin and Spyracopoulos in J Biomol NMR 53: 125-138, 2012). In this study, we demonstrate that classical line shape analysis applied to a single set of (1)H-(15)N 2D HSQC NMR spectra acquired using precise protein-ligand chemical shift titration methods we developed, produces accurate and precise kinetic parameters such as the off-rate (k ( off )). For experimentally determined kinetics in the fast exchange regime on the NMR timescale, k ( off ) ~ 3,000 s(-1) in this work, the accuracy of classical line shape analysis was determined to be better than 5 % by conducting quantum mechanical NMR simulations of the chemical shift titration methods with the magnetic resonance toolkit GAMMA. Using Monte Carlo simulations, the experimental precision for k ( off ) from line shape analysis of NMR spectra was determined to be 13 %, in agreement with the theoretical precision of 12 % from line shape analysis of the GAMMA simulations in the presence of noise and protein concentration errors. In addition, GAMMA simulations were employed to demonstrate that line shape analysis has the potential to provide reasonably accurate and precise k ( off ) values over a wide range, from 100 to 15,000 s(-1). The validity of line shape analysis for k ( off ) values approaching intermediate exchange (~100 s(-1)), may be facilitated by more accurate K ( D ) measurements from NMR-monitored chemical shift titrations, for which the dependence of K ( D ) on the chemical shift difference (Δω) between free and bound states is extrapolated to Δω = 0. The demonstrated accuracy and precision for k ( off ) will be valuable for the interpretation of biological kinetics in weakly interacting protein-protein networks, where a small change in the magnitude of the underlying kinetics of a given pathway may lead to large changes in the associated downstream signaling cascade.

Increased precision for analysis of protein-ligand dissociation constants determined from chemical shift titrations.

NMR is ideally suited for the analysis of protein-protein and protein ligand interactions with dissociation constants ranging from ~2 µM to ~1 mM, and with kinetics in the fast exchange regime on the NMR timescale. For the determination of dissociation constants (K ( D )) of 1:1 protein-protein or protein-ligand interactions using NMR, the protein and ligand concentrations must necessarily be similar in magnitude to the K ( D ), and nonlinear least squares analysis of chemical shift changes as a function of ligand concentration is employed to determine estimates for the parameters K ( D ) and the maximum chemical shift change (Δδ(max)). During a typical NMR titration, the initial protein concentration, [P (0)], is held nearly constant. For this condition, to determine the most accurate parameters for K ( D ) and Δδ(max) from nonlinear least squares analyses requires initial protein concentrations that are ~0.5 × K ( D ), and a maximum concentration for the ligand, or titrant, of ~10 × [P (0)]. From a practical standpoint, these requirements are often difficult to achieve. Using Monte Carlo simulations, we demonstrate that co-variation of the ligand and protein concentrations during a titration leads to an increase in the precision of the fitted K ( D ) and Δδ(max) values when [P (0)] > K ( D ). Importantly, judicious choice of protein and ligand concentrations for a given NMR titration, combined with nonlinear least squares analyses using two independent variables (ligand and protein concentrations) and two parameters (K ( D ) and Δδ(max)) is a straightforward approach to increasing the accuracy of measured dissociation constants for 1:1 protein-ligand interactions.

Catalytic proficiency of ubiquitin conjugation enzymes: balancing pK(a) suppression, entropy, and electrostatics.

Biological organisms orchestrate coordinated responses to external stimuli through temporal fluctuations in protein-protein interaction networks using molecular mechanisms such as the synthesis and recognition of polyubiquitin (polyUb) chains on signaling adaptor proteins. One of the pivotal chemical steps in ubiquitination involves reaction of a lysine amino group with a thioester group on an activated E2, or ubiquitin conjugation enzyme, to form an amide bond between Ub and a target protein. In this study, we demonstrate a nominal 14-fold range for the rate of the chemical step, k(cat), catalyzed by different E2 enzymes using non-steady-state, single-turnover assays. However, the observed range for k(cat) is as large as ∼100-fold for steady-state, single-turnover assays. Biochemical assays were used in combination with measurement of the underlying protein-protein interaction kinetics using NMR line-shape and ZZ-exchange analyses to determine the rate of polyUb chain synthesis catalyzed by the heterodimeric E2 enzyme Ubc13-Mms2. Modest variations in substrate affinity and k(cat) can achieve functional diversity in E2 mechanism, thereby influencing the biological outcomes of polyubiquitination. E2 enzymes achieve reaction rate enhancements through electrostatic effects such as suppression of substrate lysine pK(a) and stabilization of transition states by the preorganized, polar enzyme active site as well as the entropic effects of binding. Importantly, modestly proficient enzymes such as E2s maintain the ability to tune reaction rates; this may confer a biological advantage for achieving specificity in the diverse cellular roles for which these enzymes are involved.

Context-dependent remodeling of structure in two large protein fragments.

Protein folding involves the formation of secondary structural elements from the primary sequence and their association with tertiary assemblies. The relation of this primary sequence to a specific folded protein structure remains a central question in structural biology. An increasing body of evidence suggests that variations in homologous sequence ranging from point mutations to substantial insertions or deletions can yield stable proteins with markedly different folds. Here we report the structural characterization of domain IV (D4) and ΔD4 (polypeptides with 222 and 160 amino acids, respectively) that differ by virtue of an N-terminal deletion of 62 amino acids (28% of the overall D4 sequence). The high-resolution crystal structures of the monomeric D4 and the dimeric ΔD4 reveal substantially different folds despite an overall conservation of secondary structure. These structures show that the formation of tertiary structures, even in extended polypeptide sequences, can be highly context dependent, and they serve as a model for structural plasticity in protein isoforms.

Mechanism for recognition of polyubiquitin chains: balancing affinity through interplay between multivalent binding and dynamics.

RAP80 plays a key role in signal transduction in the DNA damage response by recruiting proteins to DNA damage foci by binding K63-polyubiquitin chains with two tandem ubiquitin-interacting motifs (tUIM). It is generally recognized that the typically weak interaction between ubiquitin (Ub) and various recognition motifs is intensified by themes such as tandem recognition motifs and Ub polymerization to achieve biological relevance. However, it remains an intricate problem to develop a detailed molecular mechanism to describe the process that leads to amplification of the Ub signal. A battery of solution-state NMR methods and molecular dynamics simulations were used to demonstrate that RAP80-tUIM employs mono- and multivalent interactions with polyUb chains to achieve enhanced affinity in comparison to monoUb interactions for signal amplification. The enhanced affinity is balanced by unfavorable entropic effects that include partial quenching of rapid reorientation between individual UIM domains and individual Ub domains in the bound state. For the RAP80-tUIM-polyUb interaction, increases in affinity with increasing chain length are a result of increased numbers of mono- and multivalent binding sites in the longer polyUb chains. The mono- and multivalent interactions are characterized by intrinsically weak binding and fast off-rates; these weak interactions with fast kinetics may be an important factor underlying the transient nature of protein-protein interactions that comprise DNA damage foci.

Dynamics of the RING domain from human TRAF6 by 15N NMR spectroscopy: implications for biological function.

Activation of transcription factor NF-kappaB requires Lys63-linked polyubiquitination of the E3 ubiquitin ligase TRAF6 via protein-protein interactions mediated by a RING domain. In this study, intra- and intermolecular chemical exchange processes of the TRAF6 RING domain were analyzed by (15)N NMR spectroscopy. Micro- to millisecond time scale motions were assessed through R 1, R 2, NOE, and cross-correlated relaxation measurements, and the kinetics of these motions were quantified with relaxation dispersion. The relaxation experiments indicate that the protein core is rigid, consistent with the functional requirement that RING domains form a binding scaffold for E2 ubiquitin conjugation enzymes. Chemical exchange is observed at the C-terminal end of the main alpha-helix of the RING domain. The C-terminal end of the main alpha-helix from the RING domain is involved in E2-E3 interactions, and modulation of slow motions for this region of the helix may be a general mechanism by which these interactions achieve ubiquitin transfer. Chemical shift mapping indicates that the TRAF6 RING domain does not self-associate in solution. Numerous RING domains are homo- or heterodimeric, and this is thought to be a functional necessity for recruitment of substrates for ubiquitination, or recruitment of multiple E2 enzymes for efficient substrate ubiquitination. However, lack of self-association for the RING domain from TRAF6, and the observation that the intact protein is a trimer, suggests that close association of RING domains within a homodimeric scaffold may not be a fundamental requirement for biological function.