Chemical exchange in biomacromolecules: past, present, and future

Activation dynamics of ubiquitin-specific protease 7

Ubiquitin-specific protease 7 (USP7) is a deubiquitinating enzyme that plays a crucial role in cellular processes, including the maintenance of genome stability and regulation of antiviral and immune responses. Its dysfunction is linked to various cancers and neurodevelopmental disorders such as Hao-Fountain syndrome. Unlike other USP-family enzymes, the triad of catalytic residues in USP7 adopts an inactive conformation and undergoes rearrangement into the active state upon substrate binding. Despite its potential importance for regulating the enzyme's activity, the dynamics of USP7 have not been explored. In this study, we combine advanced CPMG NMR relaxation dispersion measurements with the analysis of enzyme kinetics to investigate the conformational dynamics of USP7 in solution and its role in enzyme activation. Our results suggest that apo-USP7 exists in a dynamic equilibrium, transiently switching between inactive and low-populated active conformations, indicating that enzyme activation can occur spontaneously, even in the absence of a substrate. Furthermore, we show that the Hao-Fountain syndrome-associated variant G392D enhances the conformational dynamics of the enzyme, leading to a significant increase in its catalytic activity. This study captures the sparsely populated, "invisible" active conformation of USP7 and demonstrates how changes in enzyme dynamics can contribute to activity, offering broader insights into enzyme function and disease mechanisms.

Unraveling structural transitions and kinetics along the fold-switching pathway of the RfaH C-terminal domain using exchange-based NMR

The bacterial transcriptional regulator RfaH comprises structurally and functionally distinct N- (NTD) and C- (CTD) terminal domains. The latter switches from a helical hairpin packed against the NTD to a five-stranded β-roll upon displacement by RNA polymerase binding. Here, we use exchange-based NMR to probe fold-switching intermediates sampled by the isolated CTD. In addition to the predominant (~76 to 77%), semistable β-roll conformation (state A), we identify four structurally and kinetically distinct states: A', B, B', and B″. State B is NMR observable with an occupancy of ~23%, exchanges slowly (τex ~ 300 ms) with the major A species, and comprises a largely unfolded ensemble with transient occupancy of helical (α5*) and β-hairpin (β1*/β2*) elements. Backbone chemical shift-based structure predictions using the program CS-ROSETTA suggest that the two transient structural elements within the B state may interact with one another to form a semicompact structure. A' (~0.35%) is an off-pathway state that exchanges rapidly (τex ~ 1 ms) with state A and likely entails a minor localized conformational change in the β1/β2 loop. State B' (~0.3%) exchanges rapidly (τex ~ 1.2 ms) with state B and exhibits downfield 15N backbone shifts (relative to B) in the α5* region indicative of reduced helicity. Finally state B″ (~0.05%) exchanges rapidly (τex ~ 0.8 to 1 ms) with either B' (linear model) or B (branched model), displays significant differences in absolute 15N chemical shift from states B and B', and likely represents a further intermediate with increased helicity along the fold-switching pathway.

Autonomous adaptive optimization of NMR experimental conditions for precise inference of minor conformational states of proteins based on chemical exchange saturation transfer

In scientific experiments where measurement sensitivity is a major limiting factor, the optimization of experimental conditions, such as measurement parameters, is essential to maximize the information obtained per unit time and the number of experiments performed. When optimization in advance is not possible because of limited prior knowledge of the system, autonomous, adaptive optimization must be implemented during the experiment. One approach to this involves sequential Bayesian optimal experimental design, which adopts mutual information as the utility function to be maximized. In this study, we applied this optimization method to the chemical exchange saturation transfer (CEST) experiment in nuclear magnetic resonance (NMR) spectroscopy, which is used to study minor but functionally important invisible states of certain molecules, such as proteins. Adaptive optimization was utilized because prior knowledge of minor states is limited. To this end, we developed an adaptive optimization system of 15N-CEST experimental conditions for proteins using Markov chain Monte Carlo (MCMC) to calculate the posterior distribution and utility function. To ensure the completion of MCMC computations within a reasonable period with sufficient precision, we developed a second-order approximation of the CEST forward model. Both simulations and actual measurements using the FF domain of the HYPA/FBP11 protein with the A39G mutation demonstrated that the adaptive method outperformed the conventional one in terms of estimation precision of minor-state parameters based on equal numbers of measurements. Because the algorithm used for the evaluation of the utility function is independent of the type of experiment, the proposed method can be applied to various spectroscopic measurements in addition to NMR, if the forward model or its approximation can be calculated sufficiently quickly.

Advances in Isotope Labeling for Solution Nucleic Acid Nuclear Magnetic Resonance Spectroscopy

The availability of structural biology methods for nucleic acid still lags behind that of proteins, as evidenced by the smaller number of structures (DNA: 2513, RNA: 1899, nucleic acid-protein complexes: 13 842, protein: 196 887) deposited in the protein database. The skewed ratio of nucleic acid structures, relative to proteins (≈1:50), is inverted with respect to the cellular output of RNA and proteins in higher organisms (≈50:1). While nuclear magnetic resonance (NMR) is an attractive biophysical tool capable of bridging this gap at the molecular level, the conformational flexibility, line broadening, and low chemical shift dispersion of nucleic acids have made the NMR method challenging, especially for structures larger than 35 nucleotides. The incorporation of NMR-active isotopes is a f strategy to combat these problems. Significant strides made to push the size limits of nucleic acid structures solved by NMR using chemoenzymatic 13C- methyl and aromatic 15N- and 19F-13C-labeling are reviewed and challenges and opportunities are evaluated. Combining these isotopic labeling patterns with superior NMR spectroscopic properties, and new DNA/RNA synthesis methods (palindrome-nicking-dependent amplification and segmental labeling and site-specific modifications by template-directed tension), may stimulate advances in NMR studies of large DNA/RNA and their complexes with important biological functions.

Quantitative and Systematic NMR Measurements of Sequence-Dependent A-T Hoogsteen Dynamics in the DNA Double Helix

The dynamic properties of DNA depend on the sequence, providing an important source of sequence-specificity in biochemical reactions. However, comprehensively measuring how these dynamics vary with sequence is challenging, especially when they involve lowly populated and short-lived conformational states. Using 1H CEST supplemented by targeted 13C R1ρ NMR experiments, we quantitatively measured Watson-Crick to Hoogsteen dynamics for an A-T base pair in 13 trinucleotide sequence contexts. The Hoogsteen population and exchange rate varied 4-fold and 16-fold, respectively, and were dependent on both the 3'- and 5'-neighbors but only weakly dependent on monovalent ion concentration (25 versus 100 mM NaCl) and pH (6.8 versus 8.0). Flexible TA and CA dinucleotide steps exhibited the highest Hoogsteen populations, and their kinetics rates strongly depended on the 3'-neighbor. In contrast, the stiffer AA and GA steps had the lowest Hoogsteen population, and their kinetics were weakly dependent on the 3'-neighbor. The Hoogsteen lifetime was especially short when G-C neighbors flanked the A-T base pair. Our results uncover a unique conformational basis for sequence-specificity in the DNA double helix and establish the utility of NMR to quantitatively and comprehensively measure sequence-dependent DNA dynamics.

Navigating the challenges of engineering composite specialized metabolite pathways in plants

Plants are a valuable source of diverse specialized metabolites with numerous applications. However, these compounds are often produced in limited quantities, particularly under unfavorable ecological conditions. To achieve sufficient levels of target metabolites, alternative strategies such as pathway engineering in heterologous systems like microbes (e.g., bacteria and fungi) or cell-free systems can be employed. Another approach is plant engineering, which aims to either enhance the native production in the original plant or reconstruct the target pathway in a model plant system. Although increasing metabolite production in the native plant is a promising strategy, these source plants are often exotic and pose significant challenges for genetic manipulation. Effective pathway engineering requires comprehensive prior knowledge of the genes and enzymes involved, as well as the precursor, intermediate, branching, and final metabolites. Thus, a thorough elucidation of the biosynthetic pathway is closely linked to successful metabolic engineering in host or model systems. In this review, we highlight recent advances in strategies for biosynthetic pathway elucidation and metabolic engineering. We focus on efforts to engineer complex, multi-step pathways that require the expression of at least eight genes for transient and three genes for stable transformation. Reports on the engineering of complex pathways in stably transformed plants remain relatively scarce. We discuss the major hurdles in pathway elucidation and strategies for overcoming them, followed by an overview of achievements, challenges, and solutions in pathway reconstitution through metabolic engineering. Recent advances including computer-based predictions offer valuable platforms for the sustainable production of specialized metabolites in plants.

17O NMR relaxation measurements for investigation of molecular dynamics in static solids using sodium nitrate as a model compound

17O NMR methods are emerging as a powerful tool for determination of structure and dynamics in materials and biological solids. We present experimental and theoretical frameworks for measurements of 17O NMR relaxation times in static solids focusing on the excitation of the central transition of the 17O spin 5/2 system. We employ 17O-enriched NaNO3 as a model compound, in which the nitrate oxygen atoms undergo 3-fold jumps. Rotating frame (T1ρ), transverse (T2) and longitudinal (T1) relaxation times as well as line shapes were measured for the central transition in the 280 to 195 K temperature range at 14.1 and 18.8 T field strengths. We conduct experimental and theoretical comparison between different relaxation methods and demonstrate the advantage of combining data from multiple relaxation time and line shape measurements to obtain a more accurate determination of the dynamics as compared to either of the techniques alone. The computational framework for relaxation of spin 5/2 nuclei is developed using the numerical integration of the Liouville - von Neumann equation.

Decoding chaperone complexes: Insights from NMR spectroscopy

Molecular chaperones play a key role in protein homeostasis by preventing misfolding and aggregation, assisting in proper protein folding, and sometimes even disaggregating formed aggregates. Chaperones achieve this through a range of transient weak protein-protein interactions, which are difficult to study using traditional structural and biophysical techniques. Nuclear magnetic resonance (NMR) spectroscopy, however, is well-suited for studying such dynamic states and interactions. A wide range of NMR experiments have been particularly valuable in understanding the mechanisms of chaperone function, as they can characterize disordered protein structures, detect weak and nonspecific interactions involving sparsely populated states, and probe the conformational dynamics of proteins and their complexes. Recent advances in NMR have significantly enhanced our knowledge of chaperone mechanisms, especially chaperone-client interactions, despite the inherent challenges posed by the flexibility and complexity of these systems. In this review, we highlight contributions of NMR to the chaperone field, focusing on the work carried out in our laboratory, which have provided insights into how chaperones maintain function within the cellular environment and interact with various protein substrates.

Conformational penalties: New insights into nucleic acid recognition

The energy cost accompanying changes in the structures of nucleic acids when they bind partner molecules is a significant but underappreciated thermodynamic contribution to binding affinity and specificity. This review highlights recent advances in measuring conformational penalties and determining their contribution to the recognition, folding, and regulatory activities of nucleic acids. Notable progress includes methods for measuring and structurally characterizing lowly populated conformational states, obtaining ensemble information in high throughput, for large macromolecular assemblies, and in complex cellular environments. Additionally, quantitative and predictive thermodynamic models have been developed that relate conformational penalties to nucleic acid-protein association and cellular activity. These studies underscore the crucial role of conformational penalties in nucleic acid recognition.

Solid-State NMR Analysis of the Dynamics of Cofactors: Comparison of Heliobacterial and Purple Bacterial Reaction Centers

Photosynthetic reaction centers (RCs) serve as natural engines converting solar energy to chemical energy. Understanding the principles of efficient charge separation and light-induced electron transfer (ET) between the chlorophyll-type pigments might guide the synthesis for artificial photosynthetic systems. We present detailed insight into the dynamics at the atomic level using solid-state NMR techniques applied to the RCs of Heliobacillus (Hb.) mobilis (HbRCs) and the purple bacterium Rhodobacter (R.) sphaeroides (PbRCs). It is assumed that heliobacteria were among the first phototrophic organisms; therefore, their RC can be regarded as ancient. They are constructed homodimerically with perfect C2 symmetry, enabling ET over both branches of cofactors. Modern RCs of R. sphaeroides wild-type (WT) have higher redox power and are functionally highly asymmetric. The dynamics of the cofactors in both RCs has been explored using nuclear hyperpolarization, induced by the solid-state photochemically induced dynamic nuclear polarization (photo-CIDNP) effect. Based on the individual incorporation of 13C positions of the cofactors (through supplementation by 13C-δ-aminolevulinic acid), photo-CIDNP magic-angle spinning (MAS) NMR experiments provide access to the local dynamics of the cofactors along the ET path over a broad range of time scales. Theoretical analysis of the dynamic deformation of these macrocycles is also discussed in terms of function. The dynamics observed in HbRCs appears to be correlated to ET. The cofactors in PbRC are significantly less dynamic than those in the HbRC. Relevance for efficiency and redox properties are discussed.

A molten globule ensemble primes Arf1-GDP for the nucleotide switch

The adenosine di-phosphate (ADP) ribosylation factor (Arf) small guanosine tri-phosphate (GTP)ases function as molecular switches to activate signaling cascades that control membrane organization in eukaryotic cells. In Arf1, the GDP/GTP switch does not occur spontaneously but requires guanine nucleotide exchange factors (GEFs) and membranes. Exchange involves massive conformational changes, including disruption of the core β-sheet. The mechanisms by which this energetically costly switch occurs remain to be elucidated. To probe the switch mechanism, we coupled pressure perturbation with nuclear magnetic resonance (NMR), Fourier Transform infra-red spectroscopy (FTIR), small-angle X-ray scattering (SAXS), fluorescence, and computation. Pressure induced the formation of a classical molten globule (MG) ensemble. Pressure also favored the GDP to GTP transition, providing strong support for the notion that the MG ensemble plays a functional role in the nucleotide switch. We propose that the MG ensemble allows for switching without the requirement for complete unfolding and may be recognized by GEFs. An MG-based switching mechanism could constitute a pervasive feature in Arfs and Arf-like GTPases, and more generally, the evolutionarily related (Ras-like small GTPases) Rags and Gα GTPases.

Insights into the A-C Mismatch Conformational Ensemble in Duplex DNA and its Role in Genetic Processes through a Structure-based Review

Knowing the conformational ensembles formed by mismatches is crucial for understanding how they are generated and repaired and how they contribute to genomic instability. Here, we review structural and energetic studies of the A-C mismatch in duplex DNA and use the information to identify critical conformational states in its ensemble and their significance in genetic processes. In the 1970s, Topal and Fresco proposed the A-C wobble stabilized by two hydrogen bonds, one requiring protonation of adenine-N1. Subsequent NMR and X-ray crystallography studies showed that the protonated A-C wobble was in dynamic equilibrium with a neutral inverted wobble. The mismatch was shown to destabilize duplex DNA in a sequence- and pH-dependent manner by 2.4-3.8 kcal/mol and to have an apparent pKa ranging between 7.2 and 7.7. The A-C mismatch conformational repertoire expanded as structures were determined for damaged and protein-bound DNA. These structures included Watson-Crick-like conformations forming through tautomerization of the bases that drive replication errors, the reverse wobble forming through rotation of the entire nucleotide proposed to increase the fidelity of DNA replication, and the Hoogsteen base-pair forming through the flipping of the adenine base which explained the unusual specificity of DNA polymerases that bypass DNA damage. Thus, the A-C mismatch ensemble encompasses various conformational states that can be selectively stabilized in response to environmental changes such as pH shifts, intermolecular interactions, and chemical modifications, and these adaptations facilitate critical biological processes. This review also highlights the utility of existing 3D structures to build ensemble models for nucleic acid motifs.

19F NMR Reveals the Dynamics of Substrate Binding and Lid Closure for Iodotyrosine Deiodinase as a Complement to Steady-State Kinetics and Crystallography

Active site lids are common features of enzymes and typically undergo conformational changes upon substrate binding to promote catalysis. Iodotyrosine deiodinase is no exception and contains a lid segment in all of its homologues from human to bacteria. The solution-state dynamics of the lid have now been characterized using 19F NMR spectroscopy with a CF3-labeled enzyme and CF3O-labeled ligands. From two-dimensional 19F-19F NMR exchange spectroscopy, interconversion rates between the free and bound states of a CF3O-substituted tyrosine (45 ± 10 s-1) and the protein label (40 ± 3 s-1) are very similar and suggest a correlation between ligand binding and conformational reorganization of the lid. Both occur at rates that are ∼100-fold faster than turnover, and therefore these steps do not limit catalysis. A simple CF3O-labeled phenol also binds to the active site and induces a conformational change in the lid segment that was not previously detectable by crystallography. Exchange rates of the ligand (130 ± 20 s-1) and protein (98 ± 8 s-1) in this example are faster than those above but remain self-consistent to affirm a correlation between ordering of the lid and binding of the ligand. Both ligands also protect the protein from limited proteolysis, as expected from their ability to stabilize a compact lid structure. However, the minimal turnover of simple phenol substrates indicates that such stabilization may be necessary but is not sufficient for efficient catalysis.

A Transient, Excited Species of the Autoinhibited α-State of the Bacterial Transcription Factor RfaH May Represent an Early Intermediate on the Fold-Switching Pathway

RfaH is a two-domain transcription factor in which the C-terminal domain switches fold from an α-helical hairpin to a β-roll upon binding the ops-paused RNA polymerase. To ascertain the presence of a sparsely populated excited state that may prime the autoinhibited resting state of RfaH for binding ops-paused RNA polymerase, we carried out a series of NMR-based exchange experiments to probe for conformational exchange on the millisecond time scale. Quantitative analysis of these data reveals exchange between major ground (∼95%) and sparsely populated excited (∼5%) states with an exchange lifetime of ∼3 ms involving residues at the interface between the N-terminal and C-terminal domains formed by the β3/β4 hairpin and helix α3 of the N-terminal domain and helices α4 and α5 of the C-terminal domain. The largest 15N backbone chemical shift differences are associated with the β3/β4 hairpin, leading us to suggest that the excited state may involve a rigid body lateral displacement/rotation away from the C-terminal domain to adopt a position similar to that seen in the active RNA polymerase-bound state. Such a rigid body reorientation would result in a reduction in the interface between the N- and C-terminal domains with the possible introduction of a cavity or cavities. This hypothesis is supported by the observation that the population of the excited species and the exchange rate of interconversion between ground and excited states are reduced at a high (2.5 kbar) pressure. Mechanistic implications for fold switching of the C-terminal domain in the context of RNA polymerase binding are discussed.

Quantitative Characterization of Chain-Flipping of Acyl Carrier Protein of Escherichia coli Using Chemical Exchange NMR

The acyl carrier protein of Escherichia coli, termed AcpP, is a prototypical example of type II fatty acid synthase systems found in many bacteria. It serves as a central hub by accepting diverse acyl moieties (4-18 carbons) and shuttling them between its multiple enzymatic partners to generate fatty acids. Prior structures of acyl-AcpPs established that thioester-linked acyl cargos are sequestered within AcpP's hydrophobic lumen. In contrast, structures of enzyme-bound acyl-AcpPs showed translocation of AcpP-tethered acyl chains into the active sites of enzymes. The mechanistic underpinnings of this conformational interplay, termed chain-flipping, are unclear. Here, using heteronuclear NMR spectroscopy, we reveal that AcpP-tethered acyl chains (6-10 carbons) spontaneously adopt lowly populated solvent-exposed conformations. To this end, we devised a new strategy to replace AcpP's thioester linkages with 15N-labeled amide bonds, which facilitated direct "visualization" of these excited states using NMR chemical exchange saturation transfer and relaxation dispersion measurements. Global fitting of the corresponding data yielded kinetic rate constants of the underlying equilibrium and populations and lifetimes of solvent-exposed states. The latter were influenced by acyl chain composition and ranged from milliseconds to submilliseconds for chains containing six, eight, and ten carbons, owing to their variable interactions with AcpP's hydrophobic core. Although transient, the exposure of AcpP-tethered acyl chains to the solvent may allow relevant enzymes to gain access to its active thioester, and the enzyme-induced selection of this conformation will culminate in the production of fatty acids.

Quantitative and systematic NMR measurements of sequence-dependent A-T Hoogsteen dynamics uncovers unique conformational specificity in the DNA double helix

The propensities to form lowly-populated short-lived conformations of DNA could vary with sequence, providing an important source of sequence-specificity in biochemical reactions. However, comprehensively measuring how these dynamics vary with sequence is challenging. Using 1H CEST and 13C R 1 ρ NMR, we measured Watson-Crick to Hoogsteen dynamics for an A-T base pair in thirteen trinucleotide sequence contexts. The Hoogsteen population and exchange rate varied 4-fold and 16-fold, respectively, and were dependent on both the 3'- and 5'-neighbors but only weakly dependent on monovalent ion concentration (25 versus 100 mM NaCl) and pH (6.8 versus 8.0). Flexible TA and CA dinucleotide steps exhibited the highest Hoogsteen populations, and their kinetics rates strongly depended on the 3'-neighbor. In contrast, the stiffer AA and GA steps had the lowest Hoogsteen population, and their kinetics were weakly dependent on the 3'-neighbor. The Hoogsteen lifetime was especially short when G-C neighbors flanked the A-T base pair. The Hoogsteen dynamics had a distinct sequence-dependence compared to duplex stability and minor groove width. Thus, our results uncover a unique source of sequence-specificity hidden within the DNA double helix in the form of A-T Hoogsteen dynamics and establish the utility of 1H CEST to quantitively measure sequence-dependent DNA dynamics.

Direct Measurement of 8OG Syn-Anti Flips in Mutagenic 8OG·A and Long-Range Damage-Dependent Hoogsteen Breathing Dynamics Using 1H CEST NMR

Elucidating how damage impacts DNA dynamics is essential for understanding the mechanisms of damage recognition and repair. Many DNA lesions alter their propensities to form low-populated and short-lived conformational states. However, NMR methods to measure these dynamics require isotopic enrichment, which is difficult for damaged nucleotides. Here, we demonstrate the utility of the 1H chemical exchange saturation transfer (CEST) NMR experiment in measuring the dynamics of oxidatively damaged 8-oxoguanine (8OG) in the mutagenic 8OGsyn·Aanti mismatch. Using 8OG-H7 as an NMR probe of the damaged base, we directly measured 8OG syn-anti flips to form a lowly populated (pop. ∼ 5%) and short-lived (lifetime ∼50 ms) nonmutagenic 8OGanti·Aanti. These exchange parameters were in quantitative agreement with values from 13C off-resonance R1ρ and CEST on the labeled partner adenine. The Watson-Crick-like 8OGsyn·Aanti mismatch also rescued the kinetics of Hoogsteen motions at distant A-T base pairs, which the G·A mismatch had slowed down. The results lend further support for 8OGanti·Aanti as a minor conformational state of 8OG·A, reveal that 8OG damage can impact Hoogsteen dynamics at a distance, and demonstrate the utility of 1H CEST for measuring damage-dependent dynamics in unlabeled DNA.

Ligand-induced protein transition state stabilization switches the binding pathway from conformational selection to induced fit

Protein-ligand complex formation is fundamental to biological function. A central question is whether proteins spontaneously adopt binding-competent conformations to which ligands bind conformational selection (CS) or whether ligands induce the binding-competent conformation induced fit (IF). Here, we resolve the CS and IF binding pathways by characterizing protein conformational dynamics over a wide range of ligand concentrations using NMR relaxation dispersion. We determined the relative flux through the two pathways using a four-state binding model that includes both CS and IF. Experiments conducted without ligand show that galectin-3 exchanges between the ground-state conformation and a high-energy conformation similar to the ligand-bound conformation, demonstrating that CS is a plausible pathway. Near-identical crystal structures of the apo and ligand-bound states suggest that the high-energy conformation in solution corresponds to the apo crystal structure. Stepwise additions of the ligand lactose induce progressive changes in the relaxation dispersions that we fit collectively to the four-state model, yielding all microscopic rate constants and binding affinities. The ligand affinity is higher for the bound-like conformation than for the ground state, as expected for CS. Nonetheless, the IF pathway contributes greater than 70% of the total flux even at low ligand concentrations. The higher flux through the IF pathway is explained by considerably higher rates of exchange between the two protein conformations in the ligand-associated state. Thus, the ligand acts to decrease the activation barrier between protein conformations in a manner reciprocal to enzymatic transition-state stabilization of reactions involving ligand transformation.

Optimization of TROSY- and anti-TROSY-based 15N CPMG relaxation dispersion experiments through phase cycling

CPMG relaxation dispersion studies of biomolecular dynamics on the μs-ms timescale can provide detailed kinetic, thermodynamic, and structural insights into function. Frequently, the 15N spin serves as the probe of choice, as uniform incorporation of the 15N isotope is facile and cost-effective, and the interpretation of the resulting data is often relatively straightforward. In conventional CPMG relaxation dispersion experiments the application of CPMG pulses with constant radiofrequency (RF) phase can lead to artifactual dispersion profiles that result from off-resonance effects, RF field inhomogeneity, and pulse miscalibration. The development of CPMG experiments with the [0013]-phase cycle has significantly reduced the impact of pulse imperfections over a greater bandwidth of frequency offsets in comparison to constant phase experiments. Application of 15N-TROSY-based CPMG schemes to studies of the dynamics of large molecules is necessary for high sensitivity, yet the correct incorporation of the [0013]-phase cycle is non-trivial. Here we present TROSY- and anti-TROSY-based 15N CPMG experiments with the [0013]-phase cycling scheme and demonstrate, through comprehensive numerical simulations and experimental validation, enhanced resistance to pulse imperfections relative to traditional schemes utilizing constant phase CPMG pulses. Notably, exchange parameters derived from the new experiments are in good agreement with those obtained using other, more established, 15N-based CPMG approaches.

Carbon-detected deuterium solid-state NMR rotating frame relaxation measurements for protein methyl groups under magic angle spinning

Deuterium rotating frame solid-state NMR relaxation measurements (2H R1ρ) are important tools in quantitative studies of molecular dynamics. We demonstrate how 2H to 13C cross-polarization (CP) approaches under 10-40 kHz magic angle spinning rates can be combined with the 2H R1ρ blocks to allow for extension of deuterium rotating frame relaxation studies to methyl groups in biomolecules. This extension permits detection on the 13C nuclei and, hence, for the achievement of site-specific resolution. The measurements are demonstrated using a nine-residue low complexity peptide with the sequence GGKGMGFGL, in which a single selective -13CD3 label is placed at the methionine residue. Carbon-detected measurements are compared with the deuterium direct-detection results, which allows for fine-tuning of experimental approaches. In particular, we show how the adiabatic respiration CP scheme and the double adiabatic sweep on the 2H and 13C channels can be combined with the 2H R1ρ relaxation rates measurement. Off-resonance 2H R1ρ measurements are investigated in addition to the on-resonance condition, as they extent the range of effective spin-locking field.

Persistence of Methionine Side Chain Mobility at Low Temperatures in a Nine-Residue Low Complexity Peptide, as Probed by 2 H Solid-State NMR

Methionine side chains are flexible entities which play important roles in defining hydrophobic interfaces. We utilize deuterium static solid-state NMR to assess rotameric inter-conversions and other dynamic modes of the methionine in the context of a nine-residue random-coil peptide (RC9) with the low-complexity sequence GGKGMGFGL. The measurements in the temperature range of 313 to 161 K demonstrate that the rotameric interconversions in the hydrated solid powder state persist to temperatures below 200 K. Removal of solvation significantly reduces the rate of the rotameric motions. We employed 2 H NMR line shape analysis, longitudinal and rotation frame relaxation, and chemical exchange saturation transfer methods and found that the combination of multiple techniques creates a significantly more refined model in comparison with a single technique. Further, we compare the most essential features of the dynamics in RC9 to two different methionine-containing systems, characterized previously. Namely, the M35 of hydrated amyloid-β1-40 in the three-fold symmetric polymorph as well as Fluorenylmethyloxycarbonyl (FMOC)-methionine amino acid with the bulky hydrophobic group. The comparison suggests that the driving force for the enhanced methionine side chain mobility in RC9 is the thermodynamic factor stemming from distributions of rotameric populations, rather than the increase in the rate constant.

Measuring Hydroxyl Exchange Rates in Glycans Using a Synergistic Combination of Saturation Transfer and Relaxation Dispersion NMR

Molecular recognition events mediated by glycans play pivotal roles in controlling the fate of diverse biological processes such as cellular communication and the immune response. The affinity of glycans for their target receptors is governed primarily by the hydrogen bonds formed by hydroxyl groups decorating the glycan surface. Hydroxyl exchange rate constants are therefore vital parameters that report on glycan structure and dynamics. Here we present a strategy for characterizing hydroxyl hydrogen/deuterium (H/D) exchange in glycans that employs a synergistic combination of 13C chemical exchange saturation transfer (CEST) and Carr-Purcell-Meiboom-Gill relaxation dispersion (CPMG) NMR methods. We show that the combination of CEST and CPMG experiments facilitates the sensitive detection of the small (∼0.1 ppm) two-bond deuterium isotope shift on a 13C nucleus when the attached hydroxyl group fluctuates between protonated and deuterated states. This shift is leveraged for measuring site-specific kinetic H/D exchange rate constants as well as thermodynamic free energies of isotope fractionation. The CEST and CPMG modules are integrated with a selective J-cross-polarization scheme that provides the flexibility for rapid characterization of H/D exchange at a specific hydroxyl site. Moreover, our approach enables the precise isothermal measurement of hydroxyl exchange rate constants without the need for cumbersome isotope labeling. The H/D exchange rate constants of three different glycans assessed using this method highlight its potential for detecting transient intra- and intermolecular hydrogen bonds. In addition, the trends in H/D exchange rate constants establish site-specific steric accessibility as a key determinant of solvent exchange dynamics in glycans.

Direct Measurement of 8OG syn-anti Flips in Mutagenic 8OG•A and Long-Range Damage-Dependent Hoogsteen Breathing Dynamics Using 1H CEST NMR

Elucidating how damage impacts DNA dynamics is essential for understanding the mechanisms of damage recognition and repair. Many DNA lesions alter the propensities to form lowly-populated and short-lived conformational states. However, NMR methods to measure these dynamics require isotopic enrichment, which is difficult for damaged nucleotides. Here, we demonstrate the utility of the 1H chemical exchange saturation transfer (CEST) NMR experiment in measuring the dynamics of oxidatively damaged 8-oxoguanine (8OG) in the mutagenic 8OGsyn•Aanti mismatch. Using 8OG-H7 as an NMR probe of the damaged base, we directly measured 8OG syn-anti flips to form a lowly-populated (pop. ~ 5%) and short-lived (lifetime ~ 50 ms) non-mutagenic 8OGanti•Aanti. These exchange parameters were in quantitative agreement with values from 13C off-resonance R1ρ and CEST on a labeled partner adenine. The Watson-Crick-like 8OGsyn•Aanti mismatch also rescued the kinetics of Hoogsteen motions at distance A-T base pairs, which the G•A mismatch had slowed down. The results lend further support for 8OGanti•Aanti as a minor conformational state of 8OG•A, reveal that 8OG damage can impact Hoogsteen dynamics at a distance, and demonstrate the utility of 1H CEST for measuring damage-dependent dynamics in unlabeled DNA.

Z-Form Adoption of Nucleic Acid is a Multi-Step Process Which Proceeds through a Melted Intermediate

The left-handed Z-conformation of nucleic acids can be adopted by both DNA and RNA when bound by Zα domains found within a variety of innate immune response proteins. Zα domains stabilize this higher-energy conformation by making specific interactions with the unique geometry of Z-DNA/Z-RNA. However, the mechanism by which a right-handed helix contorts to become left-handed in the presence of proteins, including the intermediate steps involved, is poorly understood. Through a combination of nuclear magnetic resonance (NMR) and other biophysical measurements, we have determined that in the absence of Zα, under low salt conditions at room temperature, d(CpG) and r(CpG) constructs show no observable evidence of transient Z-conformations greater than 0.5% on either the intermediate or slow NMR time scales. At higher temperatures, we observed a transient unfolded intermediate. The ease of melting a nucleic acid duplex correlates with Z-form adoption rates in the presence of Zα. The largest contributing factor to the activation energies of Z-form adoption as calculated by Arrhenius plots is the ease of flipping the sugar pucker, as required for Z-DNA and Z-RNA. Together, these data validate the previously proposed "zipper model" for Z-form adoption in the presence of Zα. Overall, Z-conformations are more likely to be adopted by double-stranded DNA and RNA regions flanked by less stable regions and by RNAs experiencing torsional/mechanical stress.

An RNA excited conformational state at atomic resolution

Sparse and short-lived excited RNA conformational states are essential players in cell physiology, disease, and therapeutic development, yet determining their 3D structures remains challenging. Combining mutagenesis, NMR spectroscopy, and computational modeling, we determined the 3D structural ensemble formed by a short-lived (lifetime ~2.1 ms) lowly-populated (~0.4%) conformational state in HIV-1 TAR RNA. Through a strand register shift, the excited conformational state completely remodels the 3D structure of the ground state (RMSD from the ground state = 7.2 ± 0.9 Å), forming a surprisingly more ordered conformational ensemble rich in non-canonical mismatches. The structure impedes the formation of the motifs recognized by Tat and the super elongation complex, explaining why this alternative TAR conformation cannot activate HIV-1 transcription. The ability to determine the 3D structures of fleeting RNA states using the presented methodology holds great promise for our understanding of RNA biology, disease mechanisms, and the development of RNA-targeting therapeutics.

Representing structures of the multiple conformational states of proteins

Biomolecules exhibit dynamic behavior that single-state models of their structures cannot fully capture. We review some recent advances for investigating multiple conformations of biomolecules, including experimental methods, molecular dynamics simulations, and machine learning. We also address the challenges associated with representing single- and multiple-state models in data archives, with a particular focus on NMR structures. Establishing standardized representations and annotations will facilitate effective communication and understanding of these complex models to the broader scientific community.

Biomolecular NMR spectroscopy in the era of artificial intelligence

Biomolecular nuclear magnetic resonance (NMR) spectroscopy and artificial intelligence (AI) have a burgeoning synergy. Deep learning-based structural predictors have forever changed structural biology, yet these tools currently face limitations in accurately characterizing protein dynamics, allostery, and conformational heterogeneity. We begin by highlighting the unique abilities of biomolecular NMR spectroscopy to complement AI-based structural predictions toward addressing these knowledge gaps. We then highlight the direct integration of deep learning approaches into biomolecular NMR methods. AI-based tools can dramatically improve the acquisition and analysis of NMR spectra, enhancing the accuracy and reliability of NMR measurements, thus streamlining experimental processes. Additionally, deep learning enables the development of novel types of NMR experiments that were previously unattainable, expanding the scope and potential of biomolecular NMR spectroscopy. Ultimately, a combination of AI and NMR promises to further revolutionize structural biology on several levels, advance our understanding of complex biomolecular systems, and accelerate drug discovery efforts.

Evolution of Cd2+ and Cu+ binding in Helix pomatia metallothioneins

Metallothioneins (MTs) are small proteins present in all kingdoms of life. Their high cysteine content enables them to bind metal ions, such as Zn2+, Cd2+, and Cu+, providing means for detoxification and metal homeostasis. Three MT isoforms with distinct metal binding preferences are present in the Roman Snail Helix pomatia. Here, we use nuclear magnetic resonance (NMR) to follow the evolution of Cd2+ and Cu+ binding from the reconstructed ancestral Stylommatophora MT to the three H. pomatia MT (HpMT) isoforms. Information obtained from [15N,1H]-HSQC spectra and T2 relaxation times are combined to describe the conformational stability of the MT-metal complexes. A well-behaved MT-metal complex adopts a unique structure and does not undergo additional conformational exchange. The ancestor to all three HpMTs forms conformationally stable Cd2+ complexes and closely resembles the Cd2+-specific HpCdMT isoform, suggesting a role in Cd2+ detoxification for the ancestral protein. All Cu+-MT complexes, including the Cu+-specific HpCuMT isoform, undergo a considerable amount of conformational exchange. The unspecific HpCd/CuMT and the Cu+-specific HpCuMT isoforms form Cu+ complexes with comparable characteristics. It is possible to follow how Cd2+ and Cu+ binding changed throughout evolution. Interestingly, Cu+ binding improved independently in the lineages leading to the unspecific and the Cu+-specific HpMT isoforms. C-terminal domains are generally less capable of coordinating the non-cognate metal ion than N-terminal domains, indicating a higher level of specialization of the C-domain. Our findings provide new insights into snail MT evolution, helping to understand the interplay between biological function and structural features toward a comprehensive understanding of metal preference.

Discerning conformational dynamics and binding kinetics of GPCRs by 19F NMR

19F NMR provides a way of monitoring conformational dynamics of G-protein coupled receptors (GPCRs) from the perspective of an ensemble. While X-ray crystallography provides exquisitely resolved high-resolution structures of specific states, it generally does not recapitulate the true ensemble of functional states. Fluorine (19F) NMR provides a highly sensitive spectroscopic window into the conformational ensemble, generally permitting the direct quantification of resolvable states. Moreover, straightforward T1- and T2-based relaxation experiments allow for the study of fluctuations within a given state and exchange between states, on timescales spanning nanoseconds to seconds. Conveniently, most biological systems are free of fluorine. Thus, via fluorinated amino acid analogues or thiol-reactive fluorinated tags, F or CF3 reporters can be site specifically incorporated into proteins of interest. In this review, fluorine labeling protocols and 19F NMR experiments will be presented, from the perspective of small molecule NMR (i.e. drug or small molecule interactions with receptors) or macromolecular NMR (i.e. conformational dynamics of receptors and receptor-G-protein complexes).

Probing excited state 1Hα chemical shifts in intrinsically disordered proteins with a triple resonance-based CEST experiment: Application to a disorder-to-order switch

Over 40% of eukaryotic proteomes and 15% of bacterial proteomes are predicted to be intrinsically disordered based on their amino acid sequence. Intrinsically disordered proteins (IDPs) exist as heterogeneous ensembles of interconverting conformations and pose a challenge to the structure-function paradigm by apparently functioning without possessing stable structural elements. IDPs play a prominent role in biological processes involving extensive intermolecular interaction networks and their inherently dynamic nature facilitates their promiscuous interaction with multiple structurally diverse partner molecules. NMR spectroscopy has made pivotal contributions to our understanding of IDPs because of its unique ability to characterize heterogeneity at atomic resolution. NMR methods such as Chemical Exchange Saturation Transfer (CEST) and relaxation dispersion have enabled the detection of 'invisible' excited states in biomolecules which are transiently and sparsely populated, yet central for function. Here, we develop a 1Hα CEST pulse sequence which overcomes the resonance overlap problem in the 1Hα-13Cα plane of IDPs by taking advantage of the superior resolution in the 1H-15N correlation spectrum. In this sequence, magnetization is transferred after 1H CEST using a triple resonance coherence transfer pathway from 1Hα (i) to 1HN(i + 1) during which the 15N(t1) and 1HN(t2) are frequency labelled. This approach is integrated with spin state-selective CEST for eliminating spurious dips in CEST profiles resulting from dipolar cross-relaxation. We apply this sequence to determine the excited state 1Hα chemical shifts of the intrinsically disordered DNA binding domain (CytRN) of the bacterial cytidine repressor (CytR), which transiently acquires a functional globally folded conformation. The structure of the excited state, calculated using 1Hα chemical shifts in conjunction with other excited state NMR restraints, is a three-helix bundle incorporating a helix-turn-helix motif that is vital for binding DNA.

ATP induces folding of ALS-causing C71G-hPFN1 and nascent hSOD1

ALS-causing C71G-hPFN1 coexists in both folded and unfolded states, while nascent hSOD1 is unfolded. So far, the mechanisms underlying their ALS-triggering potential remain enigmatic. Here we show by NMR that ATP completely converts C71G-hPFN1 into the folded state at a 1:2 ratio, while inducing nascent hSOD1 into two co-existing states at a 1:8 ratio. Surprisingly, the inducing capacity of ATP comes from its triphosphate, but free triphosphate triggers aggregation. The inducing capacity ranks as: ATP = ATPP = PPP > ADP = AMP-PNP = AMP-PCP = PP, while AMP, adenosine, P, and NaCl show no conversion. Mechanistically, ATP and triphosphate appear to enhance the intrinsic folding capacity encoded in the sequences, as unveiled by comparing conformations and dynamics of ATP- and Zn2+-induced hSOD1 folded states. Our study provides a mechanism for the finding that some single-cell organisms employ polyphosphates as primordial chaperones, and sheds light on the enigma of age-related onset of familial ALS and risk increase of neurodegenerative diseases.

Studying micro to millisecond protein dynamics using simple amide 15N CEST experiments supplemented with major-state R2 and visible peak-position constraints

Over the last decade amide 15N CEST experiments have emerged as a popular tool to study protein dynamics that involves exchange between a 'visible' major state and sparsely populated 'invisible' minor states. Although initially introduced to study exchange between states that are in slow exchange with each other (typical exchange rates of, 10 to 400 s-1), they are now used to study interconversion between states on the intermediate to fast exchange timescale while still using low to moderate (5 to 350 Hz) 'saturating' B1 fields. The 15N CEST experiment is very sensitive to exchange as the exchange delay TEX can be quite long (~0.5 s) allowing for a large number of exchange events to occur making it a very powerful tool to detect minor sates populated ([Formula: see text]) to as low as 1%. When systems are in fast exchange and the 15N CEST data has to be described using a model that contains exchange, the exchange parameters are often poorly defined because the [Formula: see text] versus [Formula: see text] and [Formula: see text] versus exchange rate ([Formula: see text]) plots can be quite flat with shallow or no minima and the analysis of such 15N CEST data can lead to wrong estimates of the exchange parameters due to the presence of 'spurious' minima. Here we show that the inclusion of experimentally derived constraints on the intrinsic transverse relaxation rates and the inclusion of visible state peak-positions during the analysis of amide 15N CEST data acquired with moderate B1 values (~50 to ~350 Hz) results in convincing minima in the [Formula: see text] versus [Formula: see text] and the [Formula: see text] versus [Formula: see text] plots even when exchange occurs on the 100 μs timescale. The utility of this strategy is demonstrated on the fast-folding Bacillus stearothermophilus peripheral subunit binding domain that folds with a rate constant ~104 s-1. Here the analysis of 15N CEST data alone results in [Formula: see text] versus [Formula: see text] and [Formula: see text] versus [Formula: see text] plots that contain shallow minima, but the inclusion of visible-state peak positions and restraints on the intrinsic transverse relaxation rates of both states during the analysis of the 15N CEST data results in pronounced minima in the [Formula: see text] versus [Formula: see text] and [Formula: see text] versus [Formula: see text] plots and precise exchange parameters even in the fast exchange regime ([Formula: see text]~5). Using this strategy we find that the folding rate constant of PSBD is invariant (~10,500 s-1) from 33.2 to 42.9 °C while the unfolding rates (~70 to ~500 s-1) and unfolded state populations (~0.7 to ~4.3%) increase with temperature. The results presented here show that protein dynamics occurring on the 10 to 104 s-1 timescale can be studied using amide 15N CEST experiments.

Deuteron off-resonance rotating frame relaxation for the characterization of slow motions in rotating and static solid-state proteins

We demonstrate the feasibility of deuterium solid-state NMR off-resonance rotating frame relaxation measurements for studies of slow motions in biomolecular solids. The pulse sequence, which includes adiabatic pulses for magnetization alignment, is illustrated for static and magic-angle spinning conditions away from rotary resonances. We apply the measurements for three systems with selective deuterium labels at methyl groups: a) a model compound, Fluorenylmethyloxycarbonyl methionine-D3 amino acid, for which the principles of the measurements and corresponding motional modeling based on rotameric interconversions are demonstrated; b) amyloid-β1-40 fibrils labeled at a single alanine methyl group located in the disordered N-terminal domain. This system has been extensively studied in prior work and here serves as a test of the method for complex biological systems. The essential features of the dynamics consist of large-scale rearrangements of the disordered N-terminal domain and the conformational exchange between the free and bound forms of the domain, the latter one due to transient interactions with the structured core of the fibrils. and c) a 15-residue helical peptide which belongs to the predicted α-helical domain near the N-terminus of apolipoprotein B. The peptide is solvated with triolein and incorporates a selectively labeled leucine methyl groups. The method permits model refinement, indicating rotameric interconversions with a distribution of rate constants.

Synergism between x-ray crystallography and NMR residual dipolar couplings in characterizing protein dynamics

The important role of structural dynamics in protein function is widely recognized. Thermal or B-factors and their anisotropy, seen in x-ray analysis of protein structures, report on the presence of atomic coordinate heterogeneity that can be attributed to motion. However, their quantitative evaluation in terms of protein dynamics by x-ray ensemble refinement remains challenging. NMR spectroscopy provides quantitative information on the amplitudes and time scales of motional processes. Unfortunately, with a few exceptions, the NMR data do not provide direct insights into the atomic details of dynamic trajectories. Residual dipolar couplings, measured by solution NMR, are very precise parameters reporting on the time-averaged bond-vector orientations and may offer the opportunity to derive correctly weighted dynamic ensembles of structures for cases where multiple high-resolution x-ray structures are available. Applications to the SARS-CoV-2 main protease, Mpro, and ubiquitin highlight this complementarity of NMR and crystallography for quantitative assessment of internal motions.

SPICY: a method for single scan rotating frame relaxometry

T _1ρ is an NMR relaxation mode that is sensitive to low frequency molecular motions, making it an especially valuable tool in biomolecular research. Here, we introduce a new method, SPICY, for measuring T_1ρ relaxation times. In contrast to conventional T_1ρ experiments, in which the sequence is repeated many times to determine the T_1ρ time, the SPICY sequence allows determination of T_1ρ within a single scan, shortening the experiment time remarkably. We demonstrate the method using ¹H T_1ρ relaxation dispersion experiments. Additionally, we combine the sequence with spatial encoding to produce 1D images in a single scan. We show that T_1ρ relaxation times obtained using the single scan approach are in good agreement with those obtained using the traditional experiments.

Robust automated backbone triple resonance NMR assignments of proteins using Bayesian-based simulated annealing

Assignment of resonances of nuclear magnetic resonance (NMR) spectra to specific atoms within a protein remains a labor-intensive and challenging task. Automation of the assignment process often remains a bottleneck in the exploitation of solution NMR spectroscopy for the study of protein structure-dynamics-function relationships. We present an approach to the assignment of backbone triple resonance spectra of proteins. A Bayesian statistical analysis of predicted and observed chemical shifts is used in conjunction with inter-spin connectivities provided by triple resonance spectroscopy to calculate a pseudo-energy potential that drives a simulated annealing search for the most optimal set of resonance assignments. Termed Bayesian Assisted Assignments by Simulated Annealing (BARASA), a C++ program implementation is tested against systems ranging in size to over 450 amino acids including examples of intrinsically disordered proteins. BARASA is fast, robust, accommodates incomplete and incorrect information, and outperforms current algorithms - especially in cases of sparse data and is sufficiently fast to allow for real-time evaluation during data acquisition.

Slow conformational changes in the rigid and highly stable chymotrypsin inhibitor 2

Slow conformational changes are often directly linked to protein function. It is however less clear how such processes may perturb the overall folding stability of a protein. We previously found that the stabilizing double mutant L49I/I57V in the small protein chymotrypsin inhibitor 2 from barley led to distributed increased nanosecond and faster dynamics. Here we asked what effects the L49I and I57V substitutions, either individually or together, have on the slow conformational dynamics of CI2. We used ¹⁵ N CPMG spin relaxation dispersion experiments to measure the kinetics, thermodynamics and structural changes associated with slow conformational change in CI2. These changes result in an excited state that is populated to 4.3% at 1 °C. As the temperature is increased the population of the excited state decreases. Structural changes in the excited state are associated with residues that interact with water molecules that have well defined positions and are found at these positions in all crystal structures of CI2. The substitutions in CI2 have only little effect on the structure of the excited state whereas the stability of the excited state to some extent follows the stability of the main state. The minor state is thus most populated for the most stable CI2 variant and least populated for the least stable variant. We hypothesize that the interactions between the substituted residues and the well-ordered water molecules links subtle structural changes around the substituted residues to the region in the protein that experience slow conformational changes. This article is protected by copyright. All rights reserved.

Nine-residue low-complexity disordered peptide as a model system, an NMR/CD study

Disordered proteins and protein segments can be crucial for biological function. In this work we present a detailed biophysical characterization of the low-complexity nine-residue peptide with the sequence GGKGMGFGL. Based on proton solution NMR chemical shifts, circular dichroism measurements, as well as the analysis of concentration dependence of NMR linewidth, proton longitudinal relaxation times, hydrogen-deuterium exchange measurements, and ¹⁵ N rotating frame NMR relaxation measurements, we conclude that the peptide is fully disordered and monomeric in solution. The peptide will serve as a model system for future structural and dynamics studies of biologically relevant disordered peptides in solution and solid states.

Phosphates form spectroscopically dark state assemblies in common aqueous solutions

Phosphates and polyphosphates play ubiquitous roles in biology as integral structural components of cell membranes and bone, or as vehicles of energy storage via adenosine triphosphate and phosphocreatine. The solution phase space of phosphate species appears more complex than previously known. We present nuclear magnetic resonance (NMR) and cryogenic transmission electron microscopy (cryo-TEM) experiments that suggest phosphate species including orthophosphates, pyrophosphates, and adenosine phosphates associate into dynamic assemblies in dilute solutions that are spectroscopically "dark." Cryo-TEM provides visual evidence of the formation of spherical assemblies tens of nanometers in size, while NMR indicates that a majority population of phosphates remain as unassociated ions in exchange with spectroscopically invisible assemblies. The formation of these assemblies is reversibly and entropically driven by the partial dehydration of phosphate groups, as verified by diffusion-ordered spectroscopy (DOSY), indicating a thermodynamic state of assembly held together by multivalent interactions between the phosphates. Molecular dynamics simulations further corroborate that orthophosphates readily cluster in aqueous solutions. This study presents the surprising discovery that phosphate-containing molecules, ubiquitously present in the biological milieu, can readily form dynamic assemblies under a wide range of commonly used solution conditions, highlighting a hitherto unreported property of phosphate's native state in biological solutions.

A study on multi-exponential inversion of nuclear magnetic resonance relaxation data using deep learning

Nuclear magnetic resonance (NMR) is a powerful tool for formation evaluation in the oil industry to determine parameters, such as pore structure, fluid saturation, and permeability of porous materials, which are critical to reservoir engineering. The inversion of the measured relaxation data is an ill-posed problem and may lead to deviations of inversion results, which may degrade the accuracy of further data analysis and evaluation. This paper proposes a deep learning method for multi-exponential inversion of NMR relaxation data to improve accuracy. Simulated NMR data are first constructed using a priori knowledge based on the signal parameters and Gaussian distribution. These data are then used to train the neural network designed to consider noise characteristics, signal decay characteristics, signal energy variations, and non-negative features of the T₂ spectra. With the validation from simulated data, the models introduced by multi-scale convolutional neural network (CNN) and attention mechanism outperform other approaches in terms of denoising and T₂ inversion. Finally, NMR measurements of rock cores are used to compare the effectiveness of the attention multi-scale convolutional neural network (ATT-CNN) model in practical applications. The results demonstrate that the proposed method based on deep learning has better performance than the regularization method.

Probing Substrate-Loaded Carrier Proteins by Nuclear Magnetic Resonance

Carrier proteins (CPs) are central actors in nonribosomal peptide synthetases (NRPSs) as they interact with all catalytic domains, and because they covalently hold the substrates and intermediates leading to the final product. Thus, how CPs and their partner domains recognize and engage with each other as a function of CP cargos is paramount to understanding and engineering NRPSs. However, rapid hydrolysis of the labile thioester bonds holding substrates challenges molecular and biophysical studies to determine the molecular mechanisms of domain recognition. In this chapter, we describe a protocol to counteract hydrolysis and study loaded carrier proteins at the atomic level with nuclear magnetic resonance (NMR) spectroscopy. The method relies on loading CPs in situ, with adenylation domains in the NMR tube, to reach substrate-loaded CPs at steady state. We describe controls and experimental readouts necessary to assess the integrity of the sample and maintain loading on CPs. Our approach provides a basis to conduct subsequent NMR experiments and obtain kinetic, thermodynamic, dynamic, and structural parameters of substrate-loaded CPs alone or in the presence of other domains.

High Pressure CPMG and CEST Reveal That Cavity Position Dictates Distinct Dynamic Disorder in the PP32 Repeat Protein

Given the central role of conformational dynamics in protein function, it is essential to characterize the time scales and structures associated with these transitions. High pressure (HP) perturbation favors transitions to excited states because they typically occupy a smaller molar volume, thus facilitating characterization of conformational dynamics. Repeat proteins, with their straightforward architecture, provide good models for probing the sequence dependence of protein conformational dynamics. Investigations of chemical exchange by 15N CPMG relaxation dispersion analysis revealed that introduction of a cavity via substitution of isoleucine 7 by alanine in the N-terminal capping motif of the pp32 leucine-rich repeat protein leads to pressure-dependent conformational exchange detected on the 500 μs-2 ms CPMG time scale. Exchange amplitude decreased from the N- to C-terminus, revealing a gradient of conformational exchange across the protein. In contrast, introduction of a cavity in the central core of pp32 via the L60A mutation led to pressure-induced exchange on a slower (>2 ms) time scale detected by 15N-CEST analysis. Excited state 15N chemical shifts indicated that in the excited state detected by HP CEST, the N-terminal region is mostly unfolded, while the core retains native-like structure. These HP chemical exchange measurements reveal that cavity position dictates exchange on distinct time scales, highlighting the subtle, yet central role of sequence in determining protein conformational dynamics.

Charge neutralization of the active site glutamates does not limit substrate binding and transport by small multidrug resistance transporter EmrE

EmrE, a small multidrug resistance transporter from Escherichia coli, confers broad-spectrum resistance to polyaromatic cations and quaternary ammonium compounds. Previous transport assays demonstrate that EmrE transports a +1 and a +2 substrate with the same stoichiometry of two protons:one cationic substrate. This suggests that EmrE substrate binding capacity is limited to neutralization of the two essential glutamates, E14_A and E14_B (one from each subunit in the antiparallel homodimer), in the primary binding site. Here, we explicitly test this hypothesis, since EmrE has repeatedly broken expectations for membrane protein structure and transport mechanism. We previously showed that EmrE can bind a +1 cationic substrate and proton simultaneously, with cationic substrate strongly associated with one E14 residue, whereas the other remains accessible to bind and transport a proton. Here, we demonstrate that EmrE can bind a +2 cation substrate and a proton simultaneously using NMR pH titrations of EmrE saturated with divalent substrates, for a net +1 charge in the transport pore. Furthermore, we find that EmrE can alternate access and transport a +2 substrate and proton at the same time. Together, these results lead us to conclude that E14 charge neutralization does not limit the binding and transport capacity of EmrE.

A "Steady-State" Relaxation Dispersion Nuclear Magnetic Resonance Experiment for Studies of Chemical Exchange in Degenerate 1H Transitions of Methyl Groups

Degenerate spin-systems consisting of magnetically equivalent nuclear spins, such as a ¹H₃ spin-system in selectively ¹³CH₃-labeled proteins, present considerable challenges for the design of Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion NMR experiments to characterize chemical exchange on the micro-to-millisecond time-scale. Several approaches have been previously proposed for the elimination of deleterious artifacts observed in methyl ¹H CPMG relaxation dispersion profiles obtained for (¹³C)¹H₃ groups. We describe an alternative, experimentally simple solution and design a "steady-state" methyl ¹H CPMG scheme, where 90° or acute-angle (<90°) ¹H radiofrequency pulses are applied after each CPMG echo in-phase with methyl ¹H magnetization, resulting in the establishment of a "steady-state" for effective rates of magnetization decay. A simple computational procedure for quantitative analysis of the "steady-state" CPMG relaxation dispersion profiles is developed. The "steady-state" CPMG methodology is applied to two protein systems where exchange between major and minor species occurs in different regimes on the chemical shift time-scale.

Optimizing frequency sampling in CEST experiments

For the past decade chemical exchange saturation transfer (CEST) experiments have been successfully applied to study exchange processes in biomolecules involving sparsely populated, transiently formed conformers. Initial implementations focused on extensive sampling of the CEST frequency domain, requiring significant measurement times. Here we show that the lengthy sampling schemes often used are not optimal and that reduced frequency sampling schedules can be developed without a priori knowledge of the exchange parameters, that only depend on the chosen B₁ field, and, to a lesser extent, on the intrinsic transverse relaxation rates of ground state spins. The reduced sampling approach described here can be used synergistically with other methods for reducing measurement times such as those that excite multiple frequencies in the CEST dimension simultaneously, or make use of non-uniform sampling of indirectly detected time domains, to further decrease measurement times. The proposed approach is validated by analysis of simulated and experimental datasets.

Ion influence on surface water dynamics and proton exchange at protein surfaces - A unified model for transverse and longitudinal NMR relaxation dispersion

In all biologically relevant media, proteins interact in the presence of surrounding ions, and such interactions are water-mediated. Water molecules play a crucial role in the restructuring of proteins in solution and indeed in their biological activity. Surface water dynamics and proton exchange at protein surfaces is investigated here using NMR relaxometry, for two well-known globular proteins, lysozyme and bovine serum albumin, with particular attention to the role of surface ions. We present a unified model of surface water dynamics and proton exchange, accounting simultaneously for the observed longitudinal and transverse relaxation rates. The most notable effect of salt (0.1 M) concerns the slow surface water dynamics, related to rare water molecules embedded in energy wells on the protein surface. This response is protein-specific. On the other hand, the proton exchange time between labile protein-protons and water-protons at the protein surface seems to be very similar for the two proteins and is insensitive to the addition of salts at the concentration studied.

Combining NMR, SAXS and SANS to characterize the structure and dynamics of protein complexes

Understanding the structure and dynamics of biological macromolecules is essential to decipher the molecular mechanisms that underlie cellular functions. The description of structure and conformational dynamics often requires the integration of complementary techniques. In this review, we highlight the utility of combining nuclear magnetic resonance (NMR) spectroscopy with small angle scattering (SAS) to characterize these challenging biomolecular systems. NMR can assess the structure and conformational dynamics of multidomain proteins, RNAs and biomolecular complexes. It can efficiently provide information on interaction surfaces, long-distance restraints and relative domain orientations at residue-level resolution. Such information can be readily combined with high-resolution structural data available on subcomponents of biomolecular assemblies. Moreover, NMR is a powerful tool to characterize the dynamics of biomolecules on a wide range of timescales, from nanoseconds to seconds. On the other hand, SAS approaches provide global information on the size and shape of biomolecules and on the ensemble of all conformations present in solution. Therefore, NMR and SAS provide complementary data that are uniquely suited to investigate dynamic biomolecular assemblies. Here, we briefly review the type of data that can be obtained by both techniques and describe different approaches that can be used to combine them to characterize biomolecular assemblies. We then provide guidelines on which experiments are best suited depending on the type of system studied, ranging from fully rigid complexes, dynamic structures that interconvert between defined conformations and systems with very high structural heterogeneity.

Retrospective study for the universal applicability of the residue-based linear free energy relationship in the two-state exchange of protein molecules

Multiprobe measurements, such as NMR and hydrogen exchange studies, can provide the equilibrium constant, K, and rate constants for forward and backward processes, k and k′, of the two-state structural changes of a polypeptide on a per-residue basis. We previously found a linear relationship between log K and log k and between log K and log k′ for the topological exchange of a 27-residue bioactive peptide. To test the general applicability of the residue-based linear free energy relationship (rbLEFR), we performed a literature search to collect residue-specific K, k, and k′ values in various exchange processes, including folding-unfolding equilibrium, coupled folding and binding of intrinsically disordered peptides, and structural fluctuations of folded proteins. The good linearity in a substantial number of the log–log plots proved that the rbLFER holds for the structural changes in a wide variety of protein-related phenomena. Among the successful cases, the hydrogen exchange study of apomyoglobin folding intermediates is particularly interesting. We found that the residues that deviated from the linear relationship corresponded to the α-helix, for which transient translocation had been identified by other experiments. Thus, the rbLFER is useful for studying the structures and energetics of the dynamic states of protein molecules.