Eric D. Watt

Characterization of Enzyme Motions by Solution NMR Relaxation Dispersion

In many enzymes, conformational changes that occur along the reaction coordinate can pose a bottleneck to the rate of conversion of substrates to products. Characterization of these rate-limiting protein motions is essential for obtaining a full understanding of enzyme-catalyzed reactions. Solution NMR experiments such as the Carr-Purcell-Meiboom-Gill (CPMG) spin-echo or off-resonance R 1rho pulse sequences enable quantitation of protein motions in the time range of microseconds to milliseconds. These experiments allow characterization of the conformational exchange rate constant, k ex, the equilibrium populations of the relevant conformations, and the chemical shift differences (Deltaomega) between the conformations. The CPMG experiments were applied to the backbone N-H positions of ribonuclease A (RNase A). To probe the role of dynamic processes in the catalytic cycle of RNase A, stable mimics of the apo enzyme (E), enzyme-substrate (ES) complex, and enzyme-product (EP) complex were formed. The results indicate that the ligand has relatively little influence on the kinetics of motion, which occurs at 1700 s (-1) and is the same as both k cat, and the product dissociation rate constant. Instead, the effect of ligand is to stabilize one of the pre-existing conformations. Thus, these NMR experiments indicate that the conformational change in RNase A is ligand-stabilized and does not appear to be ligand-induced. Further evidence for the coupling of motion and enzyme function comes from the similar solvent deuterium kinetic isotope effect on k ex derived from the NMR measurements and k cat from enzyme kinetic studies. This isotope effect of approximately 2 depends linearly on solvent deuterium content suggesting the involvement of a single proton in RNase A motion and function. Moreover, mutation of His48 to alanine eliminates motion in RNase A and decreases the catalytic turnover rate indicating the involvement of His48, which is far from the active site, in coupling motion and function. For the enzyme triosephosphate isomerase (TIM), the opening and closing motion of a highly conserved active site loop (loop 6) has been implicated in many studies to play an important role in the catalytic cycle of the enzyme. Off-resonance R 1rho experiments were performed on TIM, and results were obtained for amino acid residues in the N-terminal (Val167), and C-terminal (Lys174, Thr177) portions of loop 6. The results indicate that all three loop residues move between the open and closed conformation at about 10,000 s (-1), which is the same as the catalytic rate constant. The O (eta) atom of Tyr208 provides a hydrogen bond to stabilize the closed form of loop 6 by interacting with the amide nitrogen of Ala176; these atoms are outside of hydrogen bonding distance in the open form of the enzyme. Mutation of Tyr208 to phenylalanine results in significant loss of catalytic activity but does not appear to alter the kex value of the N-terminal part of loop 6. Instead, removal of this hydrogen bond appears to result in an increase in the equilibrium population of the open conformer of loop 6, thereby resulting in a loss of activity through a shift in the conformational equilibrium of loop 6. Solution NMR relaxation dispersion experiments are powerful experimental tools that can elucidate protein motions with atomic resolution and can provide insight into the role of these motions in biological function.

The mechanism of rate-limiting motions in enzyme function

The ability to use conformational flexibility is a hallmark of enzyme function. Here we show that protein motions and catalytic activity in a RNase are coupled and display identical solvent isotope effects. Solution NMR relaxation experiments identify a cluster of residues, some distant from the active site, that are integral to this motion. These studies implicate a single residue, histidine-48, as the key modulator in coupling protein motion with enzyme function. Mutation of H48 to alanine results in loss of protein motion in the isotope-sensitive region of the enzyme. In addition, kcat decreases for this mutant and the kinetic solvent isotope effect on kcat, which was 2.0 in WT, is near unity in H48A. Despite being located 18 Å from the enzyme active site, H48 is essential in coordinating the motions involved in the rate-limiting enzymatic step. These studies have identified, of ≈160 potential exchangeable protons, a single site that is integral in the rate-limiting step in RNase A enzyme function.

Hydrogen exchange of monomeric α‐synuclein shows unfolded structure persists at physiological temperature and is independent of molecular crowding in Escherichia coli

Amide proton NMR signals from the N‐terminal domain of monomeric α‐synuclein (αS) are lost when the sample temperature is raised from 10°C to 35°C at pH 7.4. Although the temperature‐induced effects have been attributed to conformational exchange caused by an increase in α‐helix structure, we show that the loss of signals is due to fast amide proton exchange. At low ionic strength, hydrogen exchange rates are faster for the N‐terminal segment of αS than for the acidic C‐terminal domain. When the salt concentration is raised to 300 mM, exchange rates increase throughout the protein and become similar for the N‐ and C‐terminal domains. This indicates that the enhanced protection of amide protons from the C‐terminal domain at low salt is electrostatic in nature. Cα chemical shift data point to <10% residual α‐helix structure at 10°C and 35°C. Conformational exchange contributions to R2 are negligible at both temperatures. In contrast to the situation in vitro, the majority of amide protons are observed at 37°C in 1H‐15N HSQC spectra of αS encapsulated within living Escherichia coli cells. Our finding that temperature effects on αS NMR spectra can be explained by hydrogen exchange obviates the need to invoke special cellular factors. The retention of signals is likely due to slowed hydrogen exchange caused by the lowered intracellular pH of high‐density E. coli cultures. Taken together, our results emphasize that αS remains predominantly unfolded at physiological temperature and pH—an important conclusion for mechanistic models of the association of αS with membranes and fibrils.

Characterizing complex dynamics in the transactivation response element apical loop and motional correlations with the bulge by NMR, molecular dynamics, and mutagenesis.

The HIV-1 transactivation response element (TAR) RNA binds a variety of proteins and is a target for developing anti-HIV therapies. TAR has two primary binding sites: a UCU bulge and a CUGGGA apical loop. We used NMR residual dipolar couplings, carbon spin relaxation (R(1) and R(2)), and relaxation dispersion (R(1rho)) in conjunction with molecular dynamics and mutagenesis to characterize the dynamics of the TAR apical loop and investigate previously proposed long-range interactions with the distant bulge. Replacement of the wild-type apical loop with a UUCG loop did not significantly affect the structural dynamics at the bulge, indicating that the apical loop and the bulge act largely as independent dynamical recognition centers. The apical loop undergoes complex dynamics at multiple timescales that are likely important for adaptive recognition: U31 and G33 undergo limited motions, G32 is highly flexible at picosecond-nanosecond timescales, and G34 and C30 form a dynamic Watson-Crick basepair in which G34 and A35 undergo a slow (approximately 30 mus) likely concerted looping in and out motion, with A35 also undergoing large amplitude motions at picosecond-nanosecond timescales. Our study highlights the power of combining NMR, molecular dynamics, and mutagenesis in characterizing RNA dynamics.

Tiered High-Throughput Screening Approach to Identify Thyroperoxidase Inhibitors Within the ToxCast Phase I and II Chemical Libraries

High-throughput screening for potential thyroid-disrupting chemicals requires a system of assays to capture multiple molecular-initiating events (MIEs) that converge on perturbed thyroid hormone (TH) homeostasis. Screening for MIEs specific to TH-disrupting pathways is limited in the U.S. Environmental Protection Agency ToxCast screening assay portfolio. To fill 1 critical screening gap, the Amplex UltraRed-thyroperoxidase (AUR-TPO) assay was developed to identify chemicals that inhibit TPO, as decreased TPO activity reduces TH synthesis. The ToxCast phase I and II chemical libraries, comprised of 1074 unique chemicals, were initially screened using a single, high concentration to identify potential TPO inhibitors. Chemicals positive in the single-concentration screen were retested in concentration-response. Due to high false-positive rates typically observed with loss-of-signal assays such as AUR-TPO, we also employed 2 additional assays in parallel to identify possible sources of nonspecific assay signal loss, enabling stratification of roughly 300 putative TPO inhibitors based upon selective AUR-TPO activity. A cell-free luciferase inhibition assay was used to identify nonspecific enzyme inhibition among the putative TPO inhibitors, and a cytotoxicity assay using a human cell line was used to estimate the cellular tolerance limit. Additionally, the TPO inhibition activities of 150 chemicals were compared between the AUR-TPO and an orthogonal peroxidase oxidation assay using guaiacol as a substrate to confirm the activity profiles of putative TPO inhibitors. This effort represents the most extensive TPO inhibition screening campaign to date and illustrates a tiered screening approach that focuses resources, maximizes assay throughput, and reduces animal use.

Development and Validation of a Computational Model for Androgen Receptor Activity.

Testing thousands of chemicals to identify potential androgen receptor (AR) agonists or antagonists would cost millions of dollars and take decades to complete using current validated methods. High-throughput in vitro screening (HTS) and computational toxicology approaches can more rapidly and inexpensively identify potential androgen-active chemicals. We integrated 11 HTS ToxCast/Tox21 in vitro assays into a computational network model to distinguish true AR pathway activity from technology-specific assay interference. The in vitro HTS assays probed perturbations of the AR pathway at multiple points (receptor binding, coregulator recruitment, gene transcription, and protein production) and multiple cell types. Confirmatory in vitro antagonist assay data and cytotoxicity information were used as additional flags for potential nonspecific activity. Validating such alternative testing strategies requires high-quality reference data. We compiled 158 putative androgen-active and -inactive chemicals from a combination of international test method validation efforts and semiautomated systematic literature reviews. Detailed in vitro assay information and results were compiled into a single database using a standardized ontology. Reference chemical concentrations that activated or inhibited AR pathway activity were identified to establish a range of potencies with reproducible reference chemical results. Comparison with existing Tier 1 AR binding data from the U.S. EPA Endocrine Disruptor Screening Program revealed that the model identified binders at relevant test concentrations (<100 μM) and was more sensitive to antagonist activity. The AR pathway model based on the ToxCast/Tox21 assays had balanced accuracies of 95.2% for agonist (n = 29) and 97.5% for antagonist (n = 28) reference chemicals. Out of 1855 chemicals screened in the AR pathway model, 220 chemicals demonstrated AR agonist or antagonist activity and an additional 174 chemicals were predicted to have potential weak AR pathway activity.

Recent advances in solid-state nuclear magnetic resonance techniques to quantify biomolecular dynamics.

The emergence of solid-state NMR (SSNMR) techniques over the last decade has greatly expanded the range of accessible biomolecules amenable to structural study.1-2 These methods were developed and first applied to a variety of microcrystalline model proteins including BPTI,3 SH3,4-5 ubiquitin,6-7 kaliotoxin,8 and GB1.9-11 These studies demonstrated that individual 13C, 15N and 1H sites throughout the entire protein could be uniquely resolved and assigned for purposes of high-resolution structure determination. Whereas traditional solid-state NMR approaches required site-specific isotopic labeling, the new developments enable entire proteins to be examined by resolving signals in multiple dimensions. These capabilities also permit the measurement of dynamic parameters through correlation spectroscopy, which increases the throughput, sensitivity, reliability and reproducibility of such measurements. A major advantage of SSNMR is that a large range of sample conditions and types can be examined, including many that are not accessible to solution NMR methods, such as membrane proteins and fibrils in addition to microcrystals. Recent examples include fibrils of the HET-s prion,12 high molecular weight oligomers of αB crystallin,13 the tetrameric assembly of the M2 peptide with bound amantadine drug in physiologically relevant bilayers,14 the trimeric membrane protein YadA,15 the photoreceptor proteorhodopsin,16 fibrils of the core domain of transthyretin,17 and beta-amyloid plaques derived from the brains of Alzheimers disease patients.18 Such studies have opened up new avenues for exploration of fundamental biophysics and biochemistry, and yielded insights into clinically relevant events in membranes and aggregated states not accessible to solution NMR or X-ray crystallography. Moreover, SSNMR experiments are applicable over a larger range of temperatures than can be accessed by solution NMR19 and thereby enable a continuum of experimental conditions from physiological temperatures to the cryogenic temperatures at which most crystal structures are solved. Another advantage of NMR methods in general is that they are non-perturbing and do not require covalent modification of the protein. Thus, the dynamic information procured from such experiments has the potential to accurately represent the native protein behavior. This review focuses on developments over the last five years, regarding methods and applications of dynamics studies by solid-state NMR over timescales ranging from nanoseconds to days.

Reengineering Rate-Limiting, Millisecond Enzyme Motions by Introduction of an Unnatural Amino Acid

Rate-limiting millisecond motions in wild-type (WT) Ribonuclease A (RNase A) are modulated by histidine 48. Here, we incorporate an unnatural amino acid, thia-methylimidazole, at this site (H48C-4MI) to investigate the effects of a single residue on protein motions over multiple timescales and on enzyme catalytic turnover. Molecular dynamics simulations reveal that H48C-4MI retains some crucial WT-like hydrogen bonding interactions but the extent of protein-wide correlated motions in the nanosecond regime is decreased relative to WT. NMR Carr-Purcell-Meiboom-Gill relaxation dispersion experiments demonstrate that millisecond conformational motions in H48C-4MI are present over a similar pH range compared to WT. Furthermore, incorporation of this nonnatural amino acid allows retention of WT-like catalytic activity over the full pH range. These studies demonstrate that the complexity of the protein energy landscape during the catalytic cycle can be maintained using unnatural amino acids, which may prove useful in enzyme design efforts.

Experimental Protein Structure Verification by Scoring with a Single, Unassigned NMR Spectrum

Standard methods for de novo protein structure determination by nuclear magnetic resonance (NMR) require time-consuming data collection and interpretation efforts. Here we present a qualitatively distinct and novel approach, called Comparative, Objective Measurement of Protein Architectures by Scoring Shifts (COMPASS), which identifies the best structures from a set of structural models by numerical comparison with a single, unassigned 2D (13)C-(13)C NMR spectrum containing backbone and side-chain aliphatic signals. COMPASS does not require resonance assignments. It is particularly well suited for interpretation of magic-angle spinning solid-state NMR spectra, but also applicable to solution NMR spectra. We demonstrate COMPASS with experimental data from four proteins--GB1, ubiquitin, DsbA, and the extracellular domain of human tissue factor--and with reconstructed spectra from 11 additional proteins. For all these proteins, with molecular mass up to 25 kDa, COMPASS distinguished the correct fold, most often within 1.5 Å root-mean-square deviation of the reference structure.

On selecting a minimal set of in vitro assays to reliably determine estrogen agonist activity

ABSTRACT The US EPA is charged with screening chemicals for their ability to be endocrine disruptors through interaction with the estrogen, androgen and thyroid axes. The agency is exploring the use of high‐throughput in vitro assays to use in the Endocrine Disruptor Screening Program (EDSP), potentially as replacements for lower‐throughput in vitro and in vivo tests. The first replacement is an integrated computational and experimental model for estrogen receptor (ER) activity, to be used as an alternative to the EDSP Tier 1 in vitro ER binding and transactivation assays and the in vivo uterotrophic bioassay. The ER agonist model uses a set of 16 in vitro assays that incorporate multiple technologies and cell lines and probe multiple points in the ER pathway. Here, we demonstrate that subsets of assays with as few as 4 assays can predict the activity of all 1811 chemicals tested with accuracy equivalent to that of the full 16‐assay model. The prediction accuracy against reference chemicals is higher than that of the full chemical set, partly because the larger set contains many chemicals that can cause a variety of types of assay interference There are multiple accurate assay subsets, allowing flexibility in the construction of a multiplexed assay battery. We also discuss the issue of challenging chemicals, i.e. those that can give false positive results in certain assays, and could hence be more problematic when only a few assays are used. Graphical abstract Figure. No Caption available. HighlightsIn vitro assays can be used to prioritize chemicals for more detailed testing.A 16‐assay estrogen receptor (ER) prioritization battery has been developed.We show that smaller subsets of assays can be used instead of the full 16‐assay set.The subset assays still need to probe multiple points in the ER signaling pathway.