Fast and accurate fitting of relaxation dispersion data using the flexible software package GLOVE

An integrated approach of NMR experiments and MD simulations visualizes structural dynamics of a cyclic multi-domain protein

Cyclization can stabilize the structure of proteins, as previously demonstrated in single-domain proteins. Although Lys48-linked polyubiquitin, a multi-domain protein, is also known to be cyclized in human cells, the structural effects of cyclization remain unclear. Here, we examined the impact of cyclization on the structural stability and dynamics of cyclic Lys48-linked diubiquitin (Ub2 ). As expected, cyclization increased the thermal stability of Ub2 and its resistance to proteolytic digestion, indicating that cyclization stabilized the structure of Ub2 . Furthermore, cyclization repressed the interdomain motion in Ub2 , but cyclic Ub2 still exhibited microsecond conformational exchange in NMR relaxation dispersion experiments. A series of long coarse-grained (CG) MD simulations visualized how cyclization slowed down the intrinsic nanosecond open-closed domain motion of Ub2 to microseconds. Thus, CG-MD analysis helped to explain the unexpected NMR relaxation results, thereby facilitating characterization of the structural stabilization of cyclic Ub2 .

Solution structure of the HOIL-1L NZF domain reveals a conformational switch regulating linear ubiquitin affinity

Attachment of polyubiquitin (poly-Ub) chains to proteins is a major posttranslational modification in eukaryotes. Linear ubiquitin chain assembly complex, consisting of HOIP (HOIL-1-interacting protein), HOIL-1L (heme-oxidized IRP2 Ub ligase 1), and SHARPIN (Shank-associated RH domain-interacting protein), specifically synthesizes "head-to-tail" poly-Ub chains, which are linked via the N-terminal methionine α-amino and C-terminal carboxylate of adjacent Ub units and are thus commonly called "linear" poly-Ub chains. Linear ubiquitin chain assembly complex-assembled linear poly-Ub chains play key roles in immune signaling and suppression of cell death and have been associated with immune diseases and cancer; HOIL-1L is one of the proteins known to selectively bind linear poly-Ub via its Npl4 zinc finger (NZF) domain. Although the structure of the bound form of the HOIL-1L NZF domain with linear di-Ub is known, several aspects of the recognition specificity remain unexplained. Here, we show using NMR and orthogonal biophysical methods, how the NZF domain evolves from a free to the specific linear di-Ub-bound state while rejecting other potential Ub species after weak initial binding. The solution structure of the free NZF domain revealed changes in conformational stability upon linear Ub binding, and interactions between the NZF core and tail revealed conserved electrostatic contacts, which were sensitive to charge modulation at a reported phosphorylation site: threonine-207. Phosphomimetic mutations reduced linear Ub affinity by weakening the integrity of the linear di-Ub-bound conformation. The described molecular determinants of linear di-Ub binding provide insight into the dynamic aspects of the Ub code and the NZF domain's role in full-length HOIL-1L.

Characterization of a transitionally occupied state and thermal unfolding of domain 1.1 of σ A factor of RNA polymerase from Bacillus subtilis

σ factors are essential parts of bacterial RNA polymerase (RNAP) as they allow to recognize promotor sequences and initiate transcription. Domain 1.1 of vegetative σ factors occupies the primary channel of RNAP and also prevents binding of the σ factor to promoter DNA alone. Here, we show that domain 1.1 of Bacillus subtilisσ A exists in more structurally distinct variants in dynamic equilibrium. The major conformation at room temperature is represented by a previously reported well-folded structure solved by nuclear magnetic resonance (NMR), but 4% of the protein molecules are present in a less thermodynamically favorable state. We show that this population increases with temperature and we predict its significant elevation at higher but still biologically relevant temperatures. We characterized the minor state of the domain 1.1 using specialized methods of NMR. We found that, in contrast to the major state, the detected minor state is partially unfolded. Its propensity to form secondary structure elements is especially decreased for the first and third α helices, while the second α helix and β strand close to the C-terminus are more stable. We also analyzed thermal unfolding of the domain 1.1 and performed functional experiments with full lengthσ A and its shortened version lacking domain 1.1 (σ A _ Δ 1.1 ). The results revealed that while full lengthσ A increases transcription activity of RNAP with increasing temperature, transcription withσ A _ Δ 1.1 remains constant. In summary, this study reveals conformational dynamics of domain 1.1 and provides a basis for studies of its interaction with RNAP and effects on transcription regulation.

Machine Learning Generation of Dynamic Protein Conformational Ensembles

Machine learning has achieved remarkable success across a broad range of scientific and engineering disciplines, particularly its use for predicting native protein structures from sequence information alone. However, biomolecules are inherently dynamic, and there is a pressing need for accurate predictions of dynamic structural ensembles across multiple functional levels. These problems range from the relatively well-defined task of predicting conformational dynamics around the native state of a protein, which traditional molecular dynamics (MD) simulations are particularly adept at handling, to generating large-scale conformational transitions connecting distinct functional states of structured proteins or numerous marginally stable states within the dynamic ensembles of intrinsically disordered proteins. Machine learning has been increasingly applied to learn low-dimensional representations of protein conformational spaces, which can then be used to drive additional MD sampling or directly generate novel conformations. These methods promise to greatly reduce the computational cost of generating dynamic protein ensembles, compared to traditional MD simulations. In this review, we examine recent progress in machine learning approaches towards generative modeling of dynamic protein ensembles and emphasize the crucial importance of integrating advances in machine learning, structural data, and physical principles to achieve these ambitious goals.

Transient On- and Off-Pathway Protein Folding Intermediate States Characterized with NMR Relaxation Dispersion

The earliest events in the folding of a protein are in general poorly understood. We used NMR R2 relaxation dispersion experiments to study transient local collapse events in the unfolded-state (U) conformational ensemble of apomyoglobin (apoMb). Local residual secondary structure (seen in regions corresponding to the A, D, E, and H helices of the folded protein) is largely unchanged over the pH range of 2.3-2.75, yet a significant pH-dependent increase in the conformational exchange contribution to the R2 relaxation rate (Rex) indicates that transient intramolecular contacts occur on a microsecond to millisecond time scale at pH 2.75. A comparison of 15N and 13CO relaxation dispersion data at pH 2.75 for residues in the A, B, G, and H regions, which participate in the earliest folding intermediates, indicates that chain collapse and secondary structure formation are rapid and concomitant. Increasingly stabilizing conditions (lower temperature, higher pH) result in the observation of a relaxation dispersion in the C, CD, and E regions of the protein, which are known to fold at later stages. Mutation of Trp14 in the A-helix region to Ala eliminates conformational exchange throughout the protein, and the mutation of hydrophobic residues in other regions results in the selective inhibition of conformational exchange in the B, G, or H regions. The R2 dispersion data for WT apoMb at pH 2.75 and 10 °C are best fit to a four-state model ABGH ⇆ AGH ⇆ U ⇆ ABCD that includes on-pathway (AGH and ABGH) and off-pathway (ABCD) transiently folded states, both of which are required to explain the behavior of the mutant proteins. The off-pathway intermediate is destabilized at higher temperatures. Our analysis provides insights into the earliest stages of apoMb folding where the collapsing polypeptide chain samples both productive and nonproductive states with stabilized secondary structure.

Structural and dynamic studies of DNA recognition by NF-κB p50 RHR homodimer using methyl NMR spectroscopy

Protein dynamics involving higher-energy sparsely populated conformational substates are frequently critical for protein function. This study describes the dynamics of the homodimer (p50)2 of the p50 Rel homology region (RHR) of the transcription factor NF-κB, using 13C relaxation dispersion experiments with specifically (13C, 1H)-labeled methyl groups of Ile (δ), Leu and Val. Free (p50)2 is highly dynamic in solution, showing μs-ms relaxation dispersion consistent with exchange between the ground state and higher energy substates. These fluctuations propagate from the DNA-binding loops through the core of the domain. The motions are damped in the presence of κB DNA, but the NMR spectra of the DNA complexes reveal multiple local conformations of the p50 RHR homodimer bound to certain κB DNA sequences. Varying the length and sequence of κB DNA revealed two factors that promote a single bound conformation for the complex: the length of the κB site in the duplex and a symmetrical sequence of guanine nucleotides at both ends of the recognition motif. The dynamic nature of the DNA-binding loops, together with the multiple bound conformations of p50 RHR with certain κB sites, is consistent with variations in the transcriptional activity of the p50 homodimer with different κB sequences.

Characteristic H3 N-tail dynamics in the nucleosome core particle, nucleosome, and chromatosome

The nucleosome core particle (NCP) comprises a histone octamer, wrapped around by ∼146-bp DNA, while the nucleosome additionally contains linker DNA. We previously showed that, in the nucleosome, H4 N-tail acetylation enhances H3 N-tail acetylation by altering their mutual dynamics. Here, we have evaluated the roles of linker DNA and/or linker histone on H3 N-tail dynamics and acetylation by using the NCP and the chromatosome (i.e., linker histone H1.4-bound nucleosome). In contrast to the nucleosome, H3 N-tail acetylation and dynamics are greatly suppressed in the NCP regardless of H4 N-tail acetylation because the H3 N-tail is strongly bound between two DNA gyres. In the chromatosome, the asymmetric H3 N-tail adopts two conformations: one contacts two DNA gyres, as in the NCP; and one contacts linker DNA, as in the nucleosome. However, the rate of H3 N-tail acetylation is similar in the chromatosome and nucleosome. Thus, linker DNA and linker histone both regulate H3-tail dynamics and acetylation.

Structural Insights into Methylated DNA Recognition by the Methyl-CpG Binding Domain of MBD6 from Arabidopsis thaliana

Cytosine methylation is an epigenetic modification essential for formation of mature heterochromatin, gene silencing, and genomic stability. In plants, methylation occurs not only at cytosine bases in CpG but also in CpHpG and CpHpH contexts, where H denotes A, T, or C. Methyl-CpG binding domain (MBD) proteins, which recognize symmetrical methyl-CpG dinucleotides and act as gene repressors in mammalian cells, are also present in plant cells, although their structural and functional properties still remain poorly understood. To fill this gap, in this study, we determined the solution structure of the MBD domain of the MBD6 protein from Arabidopsis thaliana and investigated its binding properties to methylated DNA by binding assays and an in-depth NMR spectroscopic analysis. The AtMBD6 MBD domain folds into a canonical MBD structure in line with its binding specificity toward methyl-CpG and possesses a DNA binding interface similar to mammalian MBD domains. Intriguingly, however, the binding affinity of the AtMBD6 MBD domain toward methyl-CpG-containing DNA was found to be much lower than that of known mammalian MBD domains. The main difference arises from the absence of positively charged residues in AtMBD6 that supposedly interact with the DNA backbone as seen in mammalian MBD/methyl-CpG-containing DNA complexes. Taken together, we have established a structural basis for methyl-CpG recognition by AtMBD6 to develop a deeper understanding how MBD proteins work as mediators of epigenetic signals in plant cells.

NMR Relaxation Dispersion Methods for the Structural and Dynamic Analysis of Quickly Interconverting, Low-Populated Conformational Substates

Most biomolecular processes involve proteins shuttling among different conformational states, particularly from highly populated ground states to the lowly populated excited states that determine the interconversion rates and biological function, and which are invisible to most structural biology techniques. These structural transitions are rare and relatively fast: happen in the millisecond-microsecond timescale (ms-μs). NMR spectroscopy can access these timescales via relaxation dispersion techniques (RD-NMR). The exchange parameters extracted from RD-NMR experiments provide pivotal information on these otherwise invisible states that reports on key properties of the high free energy, reactive regions of the protein's energy landscape, including the mechanisms of folding/unfolding and of the interconversion between active and inactive states. Here, we describe a simple, step-by-step protocol to carry out RD-NMR experiments on proteins to detect the existence of such conformational substates and characterize their structural properties (chemical shifts).

Rigorous analysis of the interaction between proteins and low water-solubility drugs by qNMR-aided NMR titration experiments

Drugs are designed and validated based on physicochemical data on their interactions with target proteins. For low water-solubility drugs, however, quantitative analysis is practically impossible without accurate estimation of precipitation. Here we combined quantitative NMR with NMR titration experiments to rigorously quantify the interaction of the low water-solubility drug pimecrolimus with its target protein FKBP12. Notably, the dissociation constants estimated with and without consideration of precipitation differed by more than tenfold. Moreover, the method enabled us to quantitate the FKBP12-pimecrolimus interaction even under a crowded condition established using the protein crowder BSA. Notably, the FKBP12-pimecrolimus interaction was slightly hampered under the crowded environment, which is explained by transient association of BSA with the drug molecules. Collectively, the described method will contribute to both quantifying the binding properties of low water-solubility drugs and to elucidating the drug behavior in complex crowded solutions including living cells.

Role of Active Site Loop Dynamics in Mediating Ligand Release from E. coli Dihydrofolate Reductase

Conformational fluctuations from ground-state to sparsely populated but functionally important excited states play a key role in enzyme catalysis. For Escherichia coli dihydrofolate reductase (DHFR), the release of the product tetrahydrofolate (THF) and oxidized cofactor NADP⁺ occurs through exchange between closed and occluded conformations of the Met20 loop. A "dynamic knockout" mutant of E. coli DHFR, where the E. coli sequence in the Met20 loop is replaced by the human sequence (N23PP/S148A), models human DHFR and is incapable of accessing the occluded conformation. ¹H and ¹⁵N CPMG relaxation dispersion analysis for the ternary product complex of the mutant enzyme with NADP⁺ and the product analogue 5,10-dideazatetrahydrofolate (ddTHF) (E:ddTHF:NADP⁺) reveals the mechanism by which NADP⁺ is released when the Met20 loop cannot undergo the closed-to-occluded conformational transition. Two excited states were observed: one related to a faster, relatively high-amplitude conformational fluctuation in areas near the active site, associated with the shuttling of the nicotinamide ring of the cofactor out of the active site, and the other to a slower process where ddTHF undergoes small-amplitude motions within the binding site that are consistent with disorder observed in a room-temperature X-ray crystal structure of the N23PP/S148A mutant protein. These motions likely arise due to steric conflict of the pterin ring of ddTHF with the ribose-nicotinamide moiety of NADP⁺ in the closed active site. These studies demonstrate that site-specific kinetic information from relaxation dispersion experiments can provide intimate details of the changes in catalytic mechanism that result from small changes in local amino acid sequence.

Kinetic Constraints in the Specific Interaction between Phosphorylated Ubiquitin and Proteasomal Shuttle Factors

Ubiquitin (Ub) specifically interacts with the Ub-associating domain (UBA) in a proteasomal shuttle factor, while the latter is involved in either proteasomal targeting or self-assembly coacervation. PINK1 phosphorylates Ub at S65 and makes Ub alternate between C-terminally relaxed (pUb_RL) and retracted conformations (pUb_RT). Using NMR spectroscopy, we show that pUb_RL but not pUb_RT preferentially interacts with the UBA from two proteasomal shuttle factors Ubqln2 and Rad23A. Yet discriminatorily, Ubqln2-UBA binds to pUb more tightly than Rad23A does and selectively enriches pUb_RL upon complex formation. Further, we determine the solution structure of the complex between Ubqln2-UBA and pUb_RL and uncover the thermodynamic basis for the stronger interaction. NMR kinetics analysis at different timescales further suggests an indued-fit binding mechanism for pUb-UBA interaction. Notably, at a relatively low saturation level, the dissociation rate of the UBA-pUb_RL complex is comparable with the exchange rate between pUb_RL and pUb_RT. Thus, a kinetic constraint would dictate the interaction between Ub and UBA, thus fine-tuning the functional state of the proteasomal shuttle factors.

Substitution of a Surface-Exposed Residue Involved in an Allosteric Network Enhances Tryptophan Synthase Function in Cells

Networks of noncovalent amino acid interactions propagate allosteric signals throughout proteins. Tryptophan synthase (TS) is an allosterically controlled bienzyme in which the indole product of the alpha subunit (αTS) is transferred through a 25 Å hydrophobic tunnel to the active site of the beta subunit (βTS). Previous nuclear magnetic resonance and molecular dynamics simulations identified allosteric networks in αTS important for its function. We show here that substitution of a distant, surface-exposed network residue in αTS enhances tryptophan production, not by activating αTS function, but through dynamically controlling the opening of the indole channel and stimulating βTS activity. While stimulation is modest, the substitution also enhances cell growth in a tryptophan-auxotrophic strain of Escherichia coli compared to complementation with wild-type αTS, emphasizing the biological importance of the network. Surface-exposed networks provide new opportunities in allosteric drug design and protein engineering, and hint at potential information conduits through which the functions of a metabolon or even larger proteome might be coordinated and regulated.

Molecular recognition and deubiquitination of cyclic K48-linked ubiquitin chains by OTUB1

Conjugation of K48-linked ubiquitin chains to intracellular proteins mainly functions as a signal for proteasomal degradation. The conjugating enzyme E2-25K synthesizes not only canonical (noncyclic) but also cyclic K48-linked ubiquitin chains. Although the cyclic conformation is expected to repress molecular recognition by ubiquitin binding proteins due to restricting the flexibility of the ubiquitin subunits in a chain, multiple proteins are reported to associate with cyclic ubiquitin chains similar to noncyclic chains. However, the molecular mechanism of how cyclic ubiquitin chains are recognized remains unclear. Here we investigated the effect of cyclization on ubiquitin-chain cleavage and molecular recognition by a K48-linkage specific deubiquitinating enzyme OTUB1 for cyclic diubiquitin by NMR spectroscopic analyses. Compared to noncyclic diubiquitin, we observed slow but unambiguously detectable cleavage of cyclic diubiquitin to monoubiquitin by OTUB1. Intriguingly, upon ubiquitin chain cleavage, cyclic diubiquitin appeared to alter its "autoinhibited" conformation to an incompletely but partially accessible conformation, induced by interaction with OTUB1 via the ubiquitin-subunit specific recognition patches and adjacent surfaces. These data imply that cyclic ubiquitin chains may exist stably in cells in spite of the presence of deubiquitinating enzymes and that these chains can be recognized by intracellular proteins in a manner distinct from that of noncyclic ubiquitin chains.

Serine protease dynamics revealed by NMR analysis of the thrombin-thrombomodulin complex

Serine proteases catalyze a multi-step covalent catalytic mechanism of peptide bond cleavage. It has long been assumed that serine proteases including thrombin carry-out catalysis without significant conformational rearrangement of their stable two-β-barrel structure. We present nuclear magnetic resonance (NMR) and hydrogen deuterium exchange mass spectrometry (HDX-MS) experiments on the thrombin-thrombomodulin (TM) complex. Thrombin promotes procoagulative fibrinogen cleavage when fibrinogen engages both the anion binding exosite 1 (ABE1) and the active site. It is thought that TM promotes cleavage of protein C by engaging ABE1 in a similar manner as fibrinogen. Thus, the thrombin-TM complex may represent the catalytically active, ABE1-engaged thrombin. Compared to apo- and active site inhibited-thrombin, we show that thrombin-TM has reduced μs-ms dynamics in the substrate binding (S1) pocket consistent with its known acceleration of protein C binding. Thrombin-TM has increased μs-ms dynamics in a β-strand connecting the TM binding site to the catalytic aspartate. Finally, thrombin-TM had doublet peaks indicative of dynamics that are slow on the NMR timescale in residues along the interface between the two β-barrels. Such dynamics may be responsible for facilitating the N-terminal product release and water molecule entry that are required for hydrolysis of the acyl-enzyme intermediate.

Transient Diffusive Interactions with a Protein Crowder Affect Aggregation Processes of Superoxide Dismutase 1 β-Barrel

Aggregate formation of superoxide dismutase 1 (SOD1) inside motor neurons is known as a major factor in onset of amyotrophic lateral sclerosis. The thermodynamic stability of the SOD1 β-barrel has been shown to decrease in crowded environments such as inside a cell, but it remains unclear how the thermodynamics of crowding-induced protein destabilization relate to SOD1 aggregation. Here we have examined the effects of a protein crowder, lysozyme, on fibril aggregate formation of the SOD1 β-barrel. We found that aggregate formation of SOD1 is decelerated even in mildly crowded solutions. Intriguingly, transient diffusive interactions with lysozyme do not significantly affect the static structure of the SOD1 β-barrel but stabilize an alternative excited "invisible" state. The net effect of crowding is to favor species off the aggregation pathway, thereby explaining the decelerated aggregation in the crowded environment. Our observations suggest that the intracellular environment may have a similar negative (inhibitory) effect on fibril formation of other amyloidogenic proteins in living cells. Deciphering how crowded intracellular environments affect aggregation and fibril formation of such disease-associated proteins will probably become central in understanding the exact role of aggregation in the etiology of these enigmatic diseases.

Structural Dynamic Heterogeneity of Polyubiquitin Subunits Affects Phosphorylation Susceptibility

Polyubiquitin is a multifunctional protein tag formed by the covalent conjugation of ubiquitin molecules. Due to the high rigidity of the ubiquitin fold, the ubiquitin moieties in a polyubiquitin chain appear to be structurally equivalent to each other. It is therefore unclear how a specific ubiquitin moiety in a chain may be preferentially recognized by some proteins, such as the kinase PINK1. Here we show that there is structural dynamic heterogeneity in the two ubiquitin moieties of K48-linked diubiquitin by NMR spectroscopic analyses. Our analyses capture subunit-asymmetric structural fluctuations that are not directly related to the closed-to-open transition of the two ubiquitin moieties in diubiquitin. Strikingly, these newly identified heterogeneous structural fluctuations may be linked to an increase in susceptibility to phosphorylation by PINK1. Coupled with the fact that there are almost no differences in static tertiary structure among ubiquitin moieties in a chain, the observed subunit-specific structural fluctuations may be an important factor that distinguishes individual ubiquitin moieties in a chain, thereby aiding both efficiency and specificity in post-translational modifications.

Structural dynamics of double-stranded DNA with epigenome modification

Modification of cytosine plays an important role in epigenetic regulation of gene expression and genome stability. Cytosine is converted to 5-methylcytosine (5mC) by DNA methyltransferase; in turn, 5mC may be oxidized to 5-hydroxymethylcytosine (5hmC) by ten-eleven translocation enzyme. The structural flexibility of DNA is known to affect the binding of proteins to methylated DNA. Here, we have carried out a semi-quantitative analysis of the dynamics of double-stranded DNA (dsDNA) containing various epigenetic modifications by combining data from imino 1H exchange and imino 1H R1ρ relaxation dispersion NMR experiments in a complementary way. Using this approach, we characterized the base-opening (kopen) and base-closing (kclose) rates, facilitating a comparison of the base-opening and -closing process of dsDNA containing cytosine in different states of epigenetic modification. A particularly striking result is the increase in the kopen rate of hemi-methylated dsDNA 5mC/C relative to unmodified or fully methylated dsDNA, indicating that the Watson-Crick base pairs undergo selective destabilization in 5mC/C. Collectively, our findings imply that the epigenetic modulation of cytosine dynamics in dsDNA mediates destabilization of the GC Watson-Crick base pair to allow base-flipping in living cells.

RING NMR dynamics: software for analysis of multiple NMR relaxation experiments

Molecular motions are fundamental to the existence of life, and NMR spectroscopy remains one of the most useful and powerful methods to measure their rates and molecular characteristics. Multiple experimental methods are available for measuring the NMR relaxation properties and these can require different methods for extracting model parameters. We present here a new software application, RING NMR Dynamics, that is designed to support analysis of multiple relaxation types. The initial release of RING NMR Dynamics supports the analysis of exponential decay experiments such as T₁ and T₂, as well as CEST and R₂ and R_1ρ relaxation dispersion. The software runs on multiple operating systems in both a command line mode and a user-friendly GUI that allows visualizing and simulating relaxation data. Interaction with another program, NMRFx Analyst, allows drilling down from the derived relaxation parameters to the raw spectral data.

Distinct conformational dynamics and allosteric networks in alpha tryptophan synthase during active catalysis

Experimental observations of enzymes under active turnover conditions have brought new insight into the role of protein motions and allosteric networks in catalysis. Many of these studies characterize enzymes under dynamic chemical equilibrium conditions, in which the enzyme is actively catalyzing both the forward and reverse reactions during data acquisition. We have previously analyzed conformational dynamics and allosteric networks of the alpha subunit of tryptophan synthase under such conditions using NMR. We have proposed that this working state represents a four to one ratio of the enzyme bound with the indole-3-glycerol phosphate substrate (E:IGP) to the enzyme bound with the products indole and glyceraldehyde-3-phosphate (E:indole:G3P). Here, we analyze the inactive D60N variant to deconvolute the contributions of the substrate- and products-bound states to the working state. While the D60N substitution itself induces small structural and dynamic changes, the D60N E:IGP and E:indole:G3P states cannot entirely account for the conformational dynamics and allosteric networks present in the working state. The act of chemical bond breakage and/or formation, or possibly the generation of an intermediate, may alter the structure and dynamics present in the working state. As the enzyme transitions from the substrate-bound to the products-bound state, millisecond conformational exchange processes are quenched and new allosteric connections are made between the alpha active site and the surface which interfaces with the beta subunit. The structural ordering of the enzyme and these new allosteric connections may be important in coordinating the channeling of the indole product into the beta subunit.

Modeling of Hidden Structures Using Sparse Chemical Shift Data from NMR Relaxation Dispersion

NMR relaxation dispersion measurements report on conformational changes occurring on the μs-ms timescale. Chemical shift information derived from relaxation dispersion can be used to generate structural models of weakly populated alternative conformational states. Current methods to obtain such models rely on determining the signs of chemical shift changes between the conformational states, which are difficult to obtain in many situations. Here, we use a "sample and select" method to generate relevant structural models of alternative conformations of the C-terminal-associated region of Escherichia coli dihydrofolate reductase (DHFR), using only unsigned chemical shift changes for backbone amides and carbonyls (¹H, ¹⁵N, and ¹³C'). We find that CS-Rosetta sampling with unsigned chemical shift changes generates a diversity of structures that are sufficient to characterize a minor conformational state of the C-terminal region of DHFR. The excited state differs from the ground state by a change in secondary structure, consistent with previous predictions from chemical shift hypersurfaces and validated by the x-ray structure of a partially humanized mutant of E. coli DHFR (N23PP/G51PEKN). The results demonstrate that the combination of fragment modeling with sparse chemical shift data can determine the structure of an alternative conformation of DHFR sampled on the μs-ms timescale. Such methods will be useful for characterizing alternative states, which can potentially be used for in silico drug screening, as well as contributing to understanding the role of minor states in biology and molecular evolution.

Algebraic expressions for Carr-Purcell-Meiboom-Gill relaxation dispersion for N-site chemical exchange

The Carr-Purcell-Meiboom-Gill (CPMG) NMR relaxation dispersion experiment measures the effective relaxation rate constant during a train of spin-echo pulse sequence elements as a function of the echo time. The CPMG experiment is a powerful method for characterizing chemical and conformational dynamic processes, termed chemical and conformational exchange, on μs-ms time scales, comparable to the experimentally accessible echo times. Approximate theoretical expressions for the effective relaxation rate constant for N-site chemical exchange have been reported (H. Koss, M. Rance, and A. G. Palmer, Biochemistry 57, 4753-4763 (2018)). Expressions for the effective relaxation rate constant have been improved by using the Cayley-Hamilton theorem to obtain simple and accurate approximations of the average Liouvillian for the CPMG experiment. The improved accuracy of the results allows efficient analyses of experimental data. In addition, the relationship is clarified between the approach of Koss and coworkers and that of Jen (1978).

Detection of key sites of dimer dissociation and unfolding initiation during activation of acid-stress chaperone HdeA at low pH

HdeA is a small acid-stress chaperone protein with a unique activity profile. At physiological pH, it forms a folded, but inactive, dimer. Below pH 3.0, HdeA unfolds and dissociates into disordered monomers, utilizing exposed hydrophobic patches to bind other unfolded proteins and prevent their irreversible aggregation. In this way, HdeA has a key role in helping pathogenic bacteria survive our acidic stomach and colonize our intestines, facilitating the spread of dysentery. Despite numerous publications on the topic, there remain questions about the mechanism by which HdeA unfolding and activation are triggered. Previous studies usually assessed HdeA unfolding over pH increments that are too far apart to gain fine detail of the process of unfolding and dimer dissociation, and often employed techniques that prevented thorough evaluation of specific regions of the protein. We used a variety of heteronuclear NMR experiments to investigate changes to backbone and side chain structure and dynamics of HdeA at four pHs between 3.0 and 2.0. We found that the long loop in the dimer interface is an early site of initiation of dimer dissociation, and that a molecular "clasp" near the disulfide bond is broken open at low pH as part, or as a trigger, of unfolding; this process also results in the separation of C-terminal helices and exposure of key hydrophobic client binding sites. Our results highlight important regions of HdeA that may have previously been overlooked because they lie too close to the disulfide bond or are thought to be too dynamic in the folded state to influence unfolding processes.

Pinpoint analysis of a protein in slow exchange using F1F2-selective ZZ-exchange spectroscopy: assignment and kinetic analysis

ZZ-exchange spectroscopy is widely used to study slow exchange processes in biomolecules, especially determination of exchange rates and assignment of minor peaks. However, if the exchange cross peaks overlap or the populations are skewed, kinetic analysis is hindered. In order to analyze slow exchange protein dynamics under such conditions, here we have developed a new method by combining ZZ-exchange and F₁F₂-selective NMR spectroscopy. We demonstrate the utility of this method by examining the monomer-dimer transition of the ubiquitin-associated domain of p62, successfully assigning the minor (monomeric) peaks and obtaining the exchange rates, which cannot be achieved by ZZ-exchange alone.

Protein Dynamics revealed by NMR Relaxation Methods

Structural biology often focuses primarily on three-dimensional structures of biological macromolecules, deposited in the Protein Data Bank (PDB). This resource is a remarkable entity for the world-wide scientific and medical communities, as well as the general public, as it is a growing translation into three-dimensional space of the vast information in genomic databases, e.g. GENBANK. There is, however, significantly more to understanding biological function than the three-dimensional coordinate space for ground-state structures of biomolecules. The vast array of biomolecules experiences natural dynamics, interconversion between multiple conformational states, and molecular recognition and allosteric events that play out on timescales ranging from picoseconds to seconds. This wide range of timescales demands ingenious and sophisticated experimental tools to sample and interpret these motions, thus enabling clearer insight into functional annotation of the PDB. NMR spectroscopy is unique in its ability to sample this range of timescales at atomic resolution and in physiologically relevant conditions using spin relaxation methods. The field is constantly expanding to provide new creative experiments, to yield more detailed coverage of timescales, and to broaden the power of interpretation and analysis methods. This review highlights the current state of the methodology and examines the extension of analysis tools for more complex experiments and dynamic models. The future for understanding protein dynamics is bright, and these extended tools bring greater compatibility with developments in computational molecular dynamics, all of which will further our understanding of biological molecular functions. These facets place NMR as a key component in integrated structural biology.

Methyl TROSY spectroscopy: A versatile NMR approach to study challenging biological systems

A major goal in structural biology is to unravel how molecular machines function in detail. To that end, solution-state NMR spectroscopy is ideally suited as it is able to study biological assemblies in a near natural environment. Based on methyl TROSY methods, it is now possible to record high-quality data on complexes that are far over 100 kDa in molecular weight. In this review, we discuss the theoretical background of methyl TROSY spectroscopy, the information that can be extracted from methyl TROSY spectra and approaches that can be used to assign methyl resonances in large complexes. In addition, we touch upon insights that have been obtained for a number of challenging biological systems, including the 20S proteasome, the RNA exosome, molecular chaperones and G-protein-coupled receptors. We anticipate that methyl TROSY methods will be increasingly important in modern structural biology approaches, where information regarding static structures is complemented with insights into conformational changes and dynamic intermolecular interactions.

Determining Binding Kinetics of Intrinsically Disordered Proteins by NMR Spectroscopy

The unique structural flexibility of intrinsically disordered proteins (IDPs) is central to their diverse functions in cellular processes. Protein-protein interactions involving IDPs are frequently transient and dynamic in nature. Nuclear magnetic resonance (NMR) spectroscopy is an especially powerful tool for characterizing the structural propensities, dynamics, and interactions of IDPs. Here we describe applications of the Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion experiment in combination with NMR titrations to characterize the kinetics and mechanisms of interactions between intrinsically disordered proteins and their targets. We illustrate the method with reference to interactions between the activation domain of the human T-cell leukemia virus type-I (HTLV-1) basic leucine zipper protein (HBZ) and its cellular binding partner, the KIX domain of the transcriptional coactivator CBP.

A Dynamic Switch in Inactive p38γ Leads to an Excited State on the Pathway to an Active Kinase

The inactive state of mitogen-activated protein kinases (MAPKs) adopts an open conformation while the active state exists in a compact form stabilized by phosphorylation. In the active state, eukaryotic kinases undergo breathing motions related to substrate binding and product release that have not previously been detected in the inactive state. However, docking interactions of partner proteins with inactive MAPK kinases exhibit allostery in binding of activating kinases. Interactions at a site distant from the activation loop are coupled to the configuration of the activation loop, suggesting that the inactive state may also undergo concerted dynamics. X-ray crystallographic studies of nonphosphorylated, inactive p38γ reveal differences in domain orientations and active site structure in the two molecules in the asymmetric unit. One molecule resembles an inactive kinase with an open active site. The second molecule has a rotation of the N-lobe that leads to partial compaction of the active site, resulting in a conformation that is intermediate between the inactive open state and the fully closed state of the activated kinase. Although the compact state of apo p38γ displays several of the features of the activated enzyme, it remains catalytically inert. In solution, the kinase fluctuates on a millisecond time scale between the open ground state and a weakly populated excited state that is similar in structure to the compact state observed in the crystal. The nuclear magnetic resonance and crystal structure data imply that interconversion between the open and compact states involves a molecular switch associated with the DFG loop.

Structural and thermodynamic basis for the recognition of the substrate-binding cleft on hen egg lysozyme by a single-domain antibody

Single-domain antibodies (VHHs or nanobodies), developed from heavy chain-only antibodies of camelids, are gaining attention as next-generation therapeutic agents. Despite their small size, the high affinity and specificity displayed by VHHs for antigen molecules rival those of IgGs. How such small antibodies achieve that level of performance? Structural studies have revealed that VHHs tend to recognize concave surfaces of their antigens with high shape-complementarity. However, the energetic contribution of individual residues located at the binding interface has not been addressed in detail, obscuring the actual mechanism by which VHHs target the concave surfaces of proteins. Herein, we show that a VHH specific for hen egg lysozyme, D3-L11, not only displayed the characteristic binding of VHHs to a concave region of the surface of the antigen, but also exhibited a distribution of energetic hot-spots like those of IgGs and conventional protein-protein complexes. The highly preorganized and energetically compact interface of D3-L11 recognizes the concave epitope with high shape complementarity by the classical lock-and-key mechanism. Our results shed light on the fundamental basis by which a particular VHH accommodate to the concave surface of an antigens with high affinity in a specific manner, enriching the mechanistic landscape of VHHs.

Conformational flexibility of adenine riboswitch aptamer in apo and bound states using NMR and an X-ray free electron laser

Riboswitches are structured cis-regulators mainly found in the untranslated regions of messenger RNA. The aptamer domain of a riboswitch serves as a sensor for its ligand, the binding of which triggers conformational changes that regulate the behavior of its expression platform. As a model system for understanding riboswitch structures and functions, the add adenine riboswitch has been studied extensively. However, there is a need for further investigation of the conformational dynamics of the aptamer in light of the recent real-time crystallographic study at room temperature (RT) using an X-ray free electron laser (XFEL) and femtosecond X-ray crystallography (SFX). Herein, we investigate the conformational motions of the add adenine riboswitch aptamer domain, in the presence or absence of adenine, using nuclear magnetic resonance relaxation measurements and analysis of RT atomic displacement factors (B-factors). In the absence of ligand, the P1 duplex undergoes a fast exchange where the overall molecule exhibits a motion at k_ex ~ 319 s^-1, based on imino signals. In the presence of ligand, the P1 duplex adopts a highly ordered conformation, with k_ex~ 83 s^-1, similar to the global motion of the molecule, excluding the loops and binding pocket, at 84 s^-1. The µs-ms motions in both the apo and bound states are consistent with RT B-factors. Reduced spatial atomic fluctuation, ~ 50%, in P1 upon ligand binding coincides with significantly attenuated temporal dynamic exchanges. The binding pocket is structured in the absence or presence of ligand, as evidenced by relatively low and similar RT B-factors. Therefore, despite the dramatic rearrangement of the binding pocket, those residues exhibit similar spatial thermal fluctuation before and after binding.

Overview of Relaxation Dispersion NMR Spectroscopy to Study Protein Dynamics and Protein-Ligand Interactions

Proteins and nucleic acids are central to all biological processes. NMR spectroscopy has proven to be excellent for studying the dynamics of these macromolecules over various timescales. Relaxation rates and heteronuclear nuclear Overhauser-effect values can resolve motion on pico- to nanosecond timescales, residual dipolar couplings provide information on submicro- to millisecond timescales, and even slower dynamics over seconds to hours can be resolved by hydrogen-exchange experiments. Relaxation dispersion experiments are especially valuable because they resolve motion on micro- to millisecond timescales, encompassing biomolecular motions associated with ligand binding, enzymatic catalysis, and domain-domain opening. These experiments provide structural, kinetic, and thermodynamic information on "invisible" excited conformational states. Relaxation dispersion can be applied not only to single biomolecules but also to protein-ligand complexes to study the kinetics and thermodynamics of association and dissociation. We review recent developments in relaxation dispersion methodology, outline the R_1ρ relaxation dispersion experiment, and discuss application to biomolecular interactions. © 2018 by John Wiley & Sons, Inc.

Millisecond Timescale Motions Connect Amino Acid Interaction Networks in Alpha Tryptophan Synthase

Tryptophan synthase is a model system for understanding allosteric regulation within enzyme complexes. Amino acid interaction networks were previously delineated in the isolated alpha subunit (αTS) in the absence of the beta subunit (βTS). The amino acid interaction networks were different between the ligand-free enzyme and the enzyme actively catalyzing turnover. Previous X-ray crystallography studies indicated only minor localized changes when ligands bind αTS, and so, structural changes alone could not explain the changes to the amino acid interaction networks. We hypothesized that the network changes could instead be related to changes in conformational dynamics. As such, we conducted nuclear magnetic resonance relaxation studies on different substrate- and products-bound complexes of αTS. Specifically, we collected ¹⁵N R₂ relaxation dispersion data that reports on microsecond-to-millisecond timescale motion of backbone amide groups. These experiments indicated that there are conformational exchange events throughout αTS. Substrate and product binding change specific motional pathways throughout the enzyme, and these pathways connect the previously identified network residues. These pathways reach the αTS/βTS binding interface, suggesting that the identified dynamic networks may also be important for communication with the βTS subunit.

Characterization of Dynamic IDP Complexes by NMR Spectroscopy

NMR spectroscopy has proven to be a key method for studying intrinsically disordered proteins (IDPs). Nonetheless, traditional NMR methods developed for solving structures of ordered protein complexes are insufficient for the full characterization of dynamic IDP complexes, where the energy landscape is broader and more rugged. Furthermore, due to their high sensitivity to environmental changes, NMR studies of IDP complexes must be conducted with extra care and the observed NMR parameters thoroughly evaluated to enable disentanglement of binding events from ensemble distribution changes. In this chapter, written for the non-NMR expert, we start out by outlining sample preparation for IDP complexes, guide through the recording and evaluation of diagnostic ¹H,¹⁵N-HSQC spectra, and delineate more sophisticated NMR strategies to follow for the particular type of complex. The most relevant experiments are then described in terms of aims, needs, pitfalls, analysis, and expected outcomes, with references to recent examples.

Intra- and inter-protein couplings of backbone motions underlie protein thiol-disulfide exchange cascade

The thioredoxin (Trx)-coupled arsenate reductase (ArsC) is a family of enzymes that catalyzes the reduction of arsenate to arsenite in the arsenic detoxification pathway. The catalytic cycle involves a series of relayed intramolecular and intermolecular thiol-disulfide exchange reactions. Structures at different reaction stages have been determined, suggesting significant conformational fluctuations along the reaction pathway. Herein, we use two state-of-the-art NMR methods, the chemical exchange saturation transfer (CEST) and the CPMG-based relaxation dispersion (CPMG RD) experiments, to probe the conformational dynamics of B. subtilis ArsC in all reaction stages, namely the enzymatic active reduced state, the intra-molecular C10-C82 disulfide-bonded intermediate state, the inactive oxidized state, and the inter-molecular disulfide-bonded protein complex with Trx. Our results reveal highly rugged energy landscapes in the active reduced state, and suggest global collective motions in both the C10-C82 disulfide-bonded intermediate and the mixed-disulfide Trx-ArsC complex.

Residue selective 15N CEST and CPMG experiments for studies of millisecond timescale protein dynamics

Proteins are intrinsically dynamic molecules and undergo exchanges among multiple conformations to perform biological functions. The CPMG relaxation dispersion and CEST experiments are two important solution NMR techniques for characterizing the conformational exchange processes on the millisecond timescale. Traditional pseudo 3D ¹⁵N CEST and CPMG experiments have certain limitations in their applications. For example, both experiments have low sensitivity for broadened resonances, and the process of optimizing sample conditions and experimental parameters are often time consuming. To overcome these limitations, we herein present a new set of residue selective ¹⁵N CEST and CPMG pulse sequences by employing the Hartmann-Hahn cross-polarization transfer of magnetization in both 1D and 2D schemes. Combined with frequency labeling in the indirect dimension using only a small number of increments, the pulse sequences in the 2D scheme can be applied on resonances in overlapped regions of the ¹H-¹⁵N HSQC spectrum. The pulse sequences were further applied on several proteins, demonstrating their advantages over the traditional CEST and CPMG experiments under specific circumstances.

NMR Measurements Reveal the Structural Basis of Transthyretin Destabilization by Pathogenic Mutations

Inherited mutations of transthyretin (TTR) destabilize its structure, leading to aggregation and familial amyloid disease. Although numerous crystal structures of wild-type (WT) and mutant TTRs have been determined, they have failed to yield a comprehensive structural explanation for destabilization by pathogenic mutations. To identify structural and dynamic variations that are not readily observed in the crystal structures, we used NMR to study WT TTR and three kinetically and/or thermodynamically destabilized pathogenic variants (V30M, L55P, and V122I). Sequence-corrected chemical shifts reveal important structural differences between WT and mutant TTR. The L55P mutation linked to aggressive early onset cardiomyopathy and polyneuropathy induces substantial structural perturbations in both the DAGH and CBEF β-sheets, whereas the V30M polyneuropathy-linked substitution perturbs primarily the CBEF sheet. In both variants, the structural perturbations propagate across the entire width of the β-sheets from the site of mutation. Structural changes caused by the V122I cardiomyopathy-associated mutation are restricted to the immediate vicinity of the mutation site, directly perturbing the subunit interfaces. NMR relaxation dispersion measurements show that WT TTR and the three pathogenic variants undergo millisecond time scale conformational fluctuations to populate a common excited state with an altered structure in the subunit interfaces. The excited state is most highly populated in L55P. The combined application of chemical shift analysis and relaxation dispersion to these pathogenic variants reveals differences in ground state structure and in the population of a transient excited state that potentially facilitates tetramer dissociation, providing new insights into the molecular mechanism by which mutations promote TTR amyloidosis.

General Expressions for Carr-Purcell-Meiboom-Gill Relaxation Dispersion for N-Site Chemical Exchange

The Carr-Purcell-Meiboom-Gill (CPMG) nuclear magnetic resonance experiment is widely used to characterize chemical exchange phenomena in biological macromolecules. Theoretical expressions for the nuclear spin relaxation rate constant for two-site chemical exchange during CPMG pulse trains valid for all time scales are well-known as are descriptions of N-site exchange in the fast limit. We have obtained theoretical expressions for N-site exchange outside of the fast limit by using approximations to an average Liouvillian describing the decay of magnetization during a CPMG pulse train. We obtain general expressions for CPMG experiments for any N-site scheme and all experimentally accessible time scales. For sufficiently slow chemical exchange, we obtain closed-form expressions for the relaxation rate constant and a general characteristic polynomial for arbitrary kinetic schemes. Furthermore, we highlight features that qualitatively characterize CPMG curves obtained for various N-site kinetic topologies, quantitatively characterize CPMG curves obtained from systems in various N-site exchange situations, and test distinguishability of kinetic models.

Hydrogen-Deuterium Exchange Profiles of Polyubiquitin Fibrils

Ubiquitin and its polymeric forms are conjugated to intracellular proteins to regulate diverse intracellular processes. Intriguingly, polyubiquitin has also been identified as a component of pathological protein aggregates associated with Alzheimer’s disease and other neurodegenerative disorders. We recently found that polyubiquitin can form amyloid-like fibrils, and that these fibrillar aggregates can be degraded by macroautophagy. Although the structural properties appear to function in recognition of the fibrils, no structural information on polyubiquitin fibrils has been reported so far. Here, we identify the core of M1-linked diubiquitin fibrils from hydrogen-deuterium exchange experiments using solution nuclear magnetic resonance (NMR) spectroscopy. Intriguingly, intrinsically flexible regions became highly solvent-protected in the fibril structure. These results indicate that polyubiquitin fibrils are formed by inter-molecular interactions between relatively flexible structural components, including the loops and edges of secondary structure elements.

Isolation and characterization of a minimal building block of polyubiquitin fibrils

As a posttranslational modifier, polyubiquitin is involved in the regulation of diverse intracellular processes; however, it is also found in pathological protein aggregates associated with Alzheimer's disease and other neurodegenerative disorders. We previously observed that various types of polyubiquitin can form amyloid-like fibrils; however, the structural properties of these polyubiquitin fibrils have not been examined at an atomic level. Here we demonstrate that a soluble intermediate species can be extracted from disulfide-conjugated diubiquitin fibrils after cleaving the disulfide bonds in the fibrils. This newly discovered molecule is structurally and physicochemically distinguishable from native ubiquitin. In addition, it is thermodynamically metastable, as demonstrated by real-time NMR measurements. Collectively, our results suggest that the fibril-derived molecule is a minimal building block of polyubiquitin fibrils that reflects their structural and physicochemical properties.

Slow Dynamics of Tryptophan-Water Networks in Proteins

Water has a profound effect on the dynamics of biomolecules and governs many biological processes, leading to the concept that function is slaved to solvent dynamics within and surrounding the biomolecule. Protein conformational changes on μs-ms time scales are frequently associated with protein function, but little is known about the behavior of protein-bound water on these time scales. Here we have used NMR relaxation dispersion measurements to probe the tryptophan indoles in the enzyme dihydrofolate reductase (DHFR). We find that during structural changes on the μs-ms time scale, large chemical shift changes are often observed for the NH proton on the indole ring, while relatively smaller chemical shift changes are observed for the ring nitrogen atom. Comparison with experimental chemical shifts and density functional theory-based chemical shift predictions show that during the structural change the tryptophan indole NHs remain bound to water, but the geometry of the protein-bound water networks changes. These results establish that relaxation dispersion measurements can indirectly probe water dynamics and indicate that water can influence, or be influenced by, protein conformational changes on the μs-ms time scale. Our data show that structurally conserved bound water molecules can play a critical role in transmitting information between functionally important regions of the protein and provide evidence that internal protein motions can be coupled through the mediation of hydrogen-bonded water bound in the protein structure.

Elucidation of potential sites for antibody engineering by fluctuation editing

Target-specific monoclonal antibodies can be routinely acquired, but the sequences of naturally acquired antibodies are not always affinity-matured and methods that increase antigen affinity are desirable. Most biophysical studies have focused on the complementary determining region (CDR), which directly contacts the antigen; however, it remains difficult to increase the affinity as much as desired. While strategies to alter the CDR to increase antibody affinity are abundant, those that target non-CDR regions are scarce. Here we describe a new method, designated fluctuation editing, which identifies potential mutation sites and engineers a high-affinity antibody based on conformational fluctuations observed by NMR relaxation dispersion. Our data show that relaxation dispersion detects important fluctuating residues that are not located in the CDR and that increase antigen–antibody affinity by point mutation. The affinity-increased mutants are shown to fluctuate less in their free form and to form a more packed structure in their antigen-bound form.

Ligand-Mediated Folding of the OmpA Periplasmic Domain from Acinetobacter baumannii

The periplasmic domain of OmpA from Acinetobacter baumannii (AbOmpA-PD) binds to diaminopimelate and anchors the outer membrane to the peptidoglycan layer in the cell wall. Although the crystal structure of AbOmpA-PD with its ligands has been reported, the mechanism of ligand-mediated folding of AbOmpA remains elusive. Here, we report that in vitro refolded apo-AbOmpA-PD in the absence of ligand exists as a mixture of two partially folded forms in solution: mostly unfolded (apo-state I) and hololike (apo-state II) states. Binding of the diaminopimelate or glycine ligand induced complete folding of AbOmpA-PD. The apo-state I was highly flexible and contained some secondary structural elements, whereas the apo-state II closely resembled the holo-state in terms of both structure and backbone dynamics, except for the ligand-binding region. ¹⁵N-relaxation-dispersion analyses for apo-state II revealed substantial motion on a millisecond timescale of residues in the H3 helix near the ligand-binding site, with this motion disappearing upon ligand binding. These results provide an insight into the ligand-mediated folding mechanism of AbOmpA-PD in solution.

Defining the Structural Basis for Allosteric Product Release from E. coli Dihydrofolate Reductase Using NMR Relaxation Dispersion

The rate-determining step in the catalytic cycle of E. coli dihydrofolate reductase is tetrahydrofolate (THF) product release, which can occur via an allosteric or an intrinsic pathway. The allosteric pathway, which becomes accessible when the reduced cofactor NADPH is bound, involves transient sampling of a higher energy conformational state, greatly increasing the product dissociation rate as compared to the intrinsic pathway that obtains when NADPH is absent. Although the kinetics of this process are known, the enzyme structure and the THF product conformation in the transiently formed excited state remain elusive. Here, we use side-chain proton NMR relaxation dispersion measurements, X-ray crystallography, and structure-based chemical shift predictions to explore the structural basis of allosteric product release. In the excited state of the E:THF:NADPH product release complex, the reduced nicotinamide ring of the cofactor transiently enters the active site where it displaces the pterin ring of the THF product. The p-aminobenzoyl-l-glutamate tail of THF remains weakly bound in a widened binding cleft. Thus, through transient entry of the nicotinamide ring into the active site, the NADPH cofactor remodels the enzyme structure and the conformation of the THF to form a weakly populated excited state that is poised for rapid product release.

Dynamics of the Extended String-Like Interaction of TFIIE with the p62 Subunit of TFIIH

General transcription factor II E (TFIIE) contains an acid-rich region (residues 378-393) in its α-subunit, comprising 13 acidic and two hydrophobic (Phe387 and Val390) residues. Upon binding to the p62 subunit of TFIIH, the acidic region adopts an extended string-like structure on the basic groove of the pleckstrin homology domain (PHD) of p62, and inserts Phe387 and Val390 into two shallow pockets in the groove. Here, we have examined the dynamics of this interaction by NMR and molecular dynamics (MD) simulations. Although alanine substitution of Phe387 and/or Val390 greatly reduced binding to PHD, the binding mode of the mutants was similar to that of the wild-type, as judged by the chemical-shift changes of the PHD. NMR relaxation dispersion profiles of the interaction exhibited large amplitudes for residues in the C-terminal half-string in the acidic region (Phe387, Glu388, Val390, Ala391, and Asp392), indicating a two-site binding mode: one corresponding to the final complex structure, and one to an off-pathway minor complex. To probe the off-pathway complex structure, an atomically detailed free-energy landscape of the binding mode was computed by all-atom multicanonical MD. The most thermodynamically stable cluster corresponded to the final complex structure. One of the next stable clusters was the off-pathway structure cluster, showing the reversed orientation of the C-terminal half-string on the PHD groove, as compared with the final structure. MD calculations elucidated that the C-terminal half-acidic-string forms encounter complexes mainly around the positive groove region with nearly two different orientations of the string, parallel and antiparallel to the final structure. Interestingly, the most encountered complexes exhibit a parallel-like orientation, suggesting that the string has a tendency to bind around the groove in the proper orientation with the aid of Phe387 and/or Val390 to proceed smoothly to the final complex structure.

F 1 F 2-selective NMR spectroscopy

Fourier transform NMR spectroscopy has provided unprecedented insight into the structure, interaction and dynamic motion of proteins and nucleic acids. Conventional biomolecular NMR relies on the acquisition of three-dimensional and four-dimensional (4D) data matrices to establish correlations between chemical shifts in the frequency domains F ₁, F ₂, F ₃ and F ₁, F ₂, F ₃, F ₄ respectively. While rich in information, these datasets require a substantial amount of acquisition time, are visually highly unintuitive, require expert knowledge to process, and sample dark and bright regions of the frequency domains equally. Here, we present an alternative approach to obtain multidimensional chemical shift correlations for biomolecules. This strategy focuses on one narrow frequency range, F ₁ F ₂, at a time and records the resulting F ₃ F ₄ correlation spectrum by two-dimensional NMR. As a result, only regions of the frequency domain that contain signals in F ₁ F ₂ ("bright regions") are sampled. F ₁ F ₂ selection is achieved by Hartmann-Hahn cross-polarization using weak radio frequency fields. This approach reveals information equivalent to that of a conventional 4D experiment, while the dimensional reduction may shorten the total acquisition time and simplifies spectral processing, interpretation and comparative analysis. Potential applicability of the F ₁ F ₂-selective approach is illustrated by de novo assignment, structural and dynamics studies of ubiquitin and fatty-acid binding protein 4 (FABP4). Further extension of this concept may spawn new selective NMR experiments to aid studies of site-specific structural dynamics, protein-protein interactions and allosteric modulation of protein structure.

Practical considerations for investigation of protein conformational dynamics by 15N R 1ρ relaxation dispersion

It is becoming increasingly apparent that proteins are not static entities and that their function often critically depends on accurate sampling of multiple conformational states in aqueous solution. Accordingly, the development of methods to study conformational states in proteins beyond their ground-state structure ("excited states") has crucial biophysical importance. Here we investigate experimental schemes for optimally probing chemical exchange processes in proteins on the micro- to millisecond timescale by ¹⁵N R _1ρ relaxation dispersion. The schemes use selective Hartmann-Hahn cross-polarization (CP) transfer for excitation, and derive peak integrals from 1D NMR spectra (Korzhnev et al. in J Am Chem Soc 127:713-721, 2005; Hansen et al. in J Am Chem Soc 131:3818-3819, 2009). Simulation and experiment collectively show that in such CP-based schemes care has to be taken to achieve accurate suppression of undesired off-resonance coherences, when using weak spin-lock fields. This then (i) ensures that relaxation dispersion profiles in the absence of chemical exchange are flat, and (ii) facilitates extraction of relaxation dispersion profiles in crowded regions of the spectrum. Further improvement in the quality of the experimental data is achieved by recording the free-induction decays in an interleaved manner and including a heating-compensation element. The reported considerations will particularly benefit the use of CP-based R _1ρ relaxation dispersion to analyze conformational exchange processes in larger proteins, where resonance line overlap becomes the main limiting factor.

Multi-probe relaxation dispersion measurements increase sensitivity to protein dynamics

Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion measurements are a valuable tool for the characterization of structural transitions on the micro-millisecond timescale. While the measurement of (15)N relaxation dispersion is now routine, the measurements with alternative nuclei remain limited. Here we report (15)N as well as (1)H R2 relaxation dispersion measurements of the N23PP/S148A "dynamic knockout" mutant of dihydrofolate reductase. The (1)H dispersion measurements are complementary to (15)N data as many additional residues are observed to have dispersive behavior for the (1)H nucleus. Simultaneous fitting of the dispersion profiles for the two nuclei increases the accuracy of exchange parameters determined for individual residues and clustered groups of residues. The different sensitivity of the two nuclei to changes in backbone torsional angles, ring currents, and hydrogen bonding effects provides important insights into the nature of the structural changes that take place during the exchange process. We observe clear evidence of direct and indirect hydrogen bond effects for the (15)N and (1)H chemical shift changes in the active-site, modulation of ring current shielding in the CD-loop and backbone torsional changes in a cluster of residues associated with the C-terminus. This work demonstrates the power of combined (1)H and (15)N probes for the study of backbone dynamics on the micro-millisecond timescale though the analysis of chemical shift changes.