Relaxation dispersion NMR spectroscopy for the study of protein allostery

Thermodynamic coupling of the tandem RRM domains of hnRNP A1 underlie its pleiotropic RNA binding functions

The functional properties of RNA binding proteins (RBPs) require allosteric regulation through interdomain communication. Despite the importance of allostery to biological regulation, only a few studies have been conducted to describe the biophysical nature by which interdomain communication manifests in RBPs. Here, we show for hnRNP A1 that interdomain communication is vital for the unique stability of its amino-terminal domain, which consists of two RNA recognition motifs (RRMs). These RRMs exhibit drastically different stability under pressure. RRM2 unfolds as an individual domain but remains stable when appended to RRM1. Variants that disrupt interdomain communication between the tandem RRMs show a significant decrease in stability. Carrying these mutations over to the full-length protein for in vivo experiments revealed that the mutations affected the ability of the disordered carboxyl-terminal domain to engage in protein-protein interactions and influenced the protein's RNA binding capacity. Collectively, this work reveals that thermodynamic coupling between the tandem RRMs of hnRNP A1 accounts for its allosteric regulatory functions.

Glycan Shape, Motions, and Interactions Explored by NMR Spectroscopy

Glycans in the form of oligosaccharides, polysaccharides, and glycoconjugates are ubiquitous in nature, and their structures range from linear assemblies to highly branched and decorated constructs. Solution state NMR spectroscopy facilitates elucidation of preferred conformations and shapes of the saccharides, motions, and dynamic aspects related to processes over time as well as the study of transient interactions with proteins. Identification of intermolecular networks at the atomic level of detail in recognition events by carbohydrate-binding proteins known as lectins, unraveling interactions with antibodies, and revealing substrate scope and action of glycosyl transferases employed for synthesis of oligo- and polysaccharides may efficiently be analyzed by NMR spectroscopy. By utilizing NMR active nuclei present in glycans and derivatives thereof, including isotopically enriched compounds, highly detailed information can be obtained by the experiments. Subsequent analysis may be aided by quantum chemical calculations of NMR parameters, machine learning-based methodologies and artificial intelligence. Interpretation of the results from NMR experiments can be complemented by extensive molecular dynamics simulations to obtain three-dimensional dynamic models, thereby clarifying molecular recognition processes involving the glycans.

Thermodynamic Coupling of the tandem RRM domains of hnRNP A1 underlie its Pleiotropic RNA Binding Functions

The functional properties of RNA-binding proteins (RBPs) require allosteric regulation through inter-domain communication. Despite the foundational importance of allostery to biological regulation, almost no studies have been conducted to describe the biophysical nature by which inter-domain communication manifests in RBPs. Here, we show through high-pressure studies with hnRNP A1 that inter-domain communication is vital for the unique stability of its N- terminal domain containing a tandem of RNA Recognition Motifs (RRMs). Despite high sequence similarity and nearly identical tertiary structures, the two RRMs exhibit drastically different stability under pressure. RRM2 unfolds completely under high-pressure as an individual domain, but when appended to RRM1, it remains stable. Variants in which inter-domain communication is disrupted between the tandem RRMs show a large decrease in stability under pressure. Carrying these mutations over to the full-length protein for in vivo experiments revealed that the mutations affected the ability of the disordered C-terminus to engage in protein-protein interactions and more importantly, they also influenced the RNA binding capacity. Collectively, this work reveals that thermodynamic coupling between the tandem RRMs of hnRNP A1 accounts for its allosteric regulatory functions.

Assessing the applicability of 19F labeled tryptophan residues to quantify protein dynamics

Nuclear magnetic resonance (NMR) spectroscopy is uniquely suited to study the dynamics of biomolecules in solution. Most NMR studies exploit the spins of proton, carbon and nitrogen isotopes, as these atoms are highly abundant in proteins and nucleic acids. As an alternative and complementary approach, fluorine atoms can be introduced into biomolecules at specific sites of interest. These labels can then be used as sensitive probes for biomolecular structure, dynamics or interactions. Here, we address if the replacement of tryptophan with 5-fluorotryptophan residues has an effect on the overall dynamics of proteins and if the introduced fluorine probe is able to accurately report on global exchange processes. For the four different model proteins (KIX, Dcp1, Dcp2 and DcpS) that we examined, we established that ¹⁵N CPMG relaxation dispersion or EXSY profiles are not affected by the 5-fluorotryptophan, indicating that this replacement of a proton with a fluorine has no effect on the protein motions. However, we found that the motions that the 5-fluorotryptophan reports on can be significantly faster than the backbone motions. This implies that care needs to be taken when interpreting fluorine relaxation data in terms of global protein motions. In summary, our results underscore the great potential of fluorine NMR methods, but also highlight potential pitfalls that need to be considered.

NMR applications to GPCR recognition by peptide ligands

Peptides form the largest group of ligands that modulate the activity of more than 120 different GPCRs. Among which linear disordered peptide ligands usually undergo significant conformational changes upon binding that is essential for receptor recognition and activation. Conformational selection and induced fit are the extreme mechanisms of coupled folding and binding that can be distinguished by analysis of binding pathways by methods that include NMR. However, the large size of GPCRs in membrane-mimetic environments limits NMR applications. In this review, we highlight advances in the field that can be adopted to address coupled folding and binding of peptide ligands to their cognate receptors.

Binding and Functional Folding (BFF): A Physiological Framework for Studying Biomolecular Interactions and Allostery

EF-hand Ca2+-binding proteins (CBPs), such as S100 proteins (S100s) and calmodulin (CaM), are signaling proteins that undergo conformational changes upon increasing intracellular Ca2+. Upon binding Ca2+, S100 proteins and CaM interact with protein targets and induce important biological responses. The Ca2+-binding affinity of CaM and most S100s in the absence of target is weak (CaKD > 1 μM). However, upon effector protein binding, the Ca2+ affinity of these proteins increases via heterotropic allostery (CaKD < 1 μM). Because of the high number and micromolar concentrations of EF-hand CBPs in a cell, at any given time, allostery is required physiologically, allowing for (i) proper Ca2+ homeostasis and (ii) strict maintenance of Ca2+-signaling within a narrow dynamic range of free Ca2+ ion concentrations, [Ca2+]free. In this review, mechanisms of allostery are coalesced into an empirical "binding and functional folding (BFF)" physiological framework. At the molecular level, folding (F), binding and folding (BF), and BFF events include all atoms in the biomolecular complex under study. The BFF framework is introduced with two straightforward BFF types for proteins (type 1, concerted; type 2, stepwise) and considers how homologous and nonhomologous amino acid residues of CBPs and their effector protein(s) evolved to provide allosteric tightening of Ca2+ and simultaneously determine how specific and relatively promiscuous CBP-target complexes form as both are needed for proper cellular function.

Dynamic Long-Range Interactions Influence Substrate Binding and Catalysis by Human Histidine Triad Nucleotide-Binding Proteins (HINTs), Key Regulators of Multiple Cellular Processes and Activators of Antiviral ProTides

Human histidine triad nucleotide-binding (hHINT) proteins catalyze nucleotide phosphoramidase and acyl-phosphatase reactions that are essential for the activation of antiviral proTides, such as Sofosbuvir and Remdesivir. hHINT1 and hHINT2 are highly homologous but exhibit disparate roles as regulators of opioid tolerance (hHINT1) and mitochondrial activity (hHINT2). NMR studies of hHINT1 reveal a pair of dynamic surface residues (Q62, E100), which gate a conserved water channel leading to the active site 13 Å away. hHINT2 crystal structures identify analogous residues (R99, D137) and water channel. hHINT1 Q62 variants significantly alter the steady-state k_cat and K_m for turnover of the fluorescent substrate (TpAd), while stopped-flow kinetics indicate that K_D also changes. hHINT2, like hHINT1, exhibits a burst phase of adenylation, monitored by fluorescent tryptamine release, prior to rate-limiting hydrolysis and nucleotide release. hHINT2 exhibits a much smaller burst-phase amplitude than hHINT1, which is further diminished in hHINT2 R99Q. Kinetic simulations suggest that amplitude variations can be accounted for by a variable fluorescent yield of the E·S complex from changes in the environment of bound TpAd. Isothermal titration calorimetry measurements of inhibitor binding show that these hHINT variants also alter the thermodynamic binding profile. We propose that these altered surface residues engender long-range dynamic changes that affect the orientation of bound ligands, altering the thermodynamic and kinetic characteristics of hHINT active site function. Thus, studies of the cellular roles and proTide activation potential by hHINTs should consider the importance of long-range interactions and possible protein binding surfaces far from the active site.

Binding Specificity of ASHH2 CW Domain Toward H3K4me1 Ligand Is Coupled to Its Structural Stability Through Its α1-Helix

The CW domain binds to histone tail modifications found in different protein families involved in epigenetic regulation and chromatin remodeling. CW domains recognize the methylation state of the fourth lysine on histone 3 and could, therefore, be viewed as a reader of epigenetic information. The specificity toward different methylation states such as me1, me2, or me3 depends on the particular CW subtype. For example, the CW domain of ASHH2 methyltransferase binds preferentially to H3K4me1, and MORC3 binds to both H3K4me2 and me3 modifications, while ZCWPW1 is more specific to H3K4me3. The structural basis for these preferential bindings is not well understood, and recent research suggests that a more complete picture will emerge if dynamical and energetic assessments are included in the analysis of interactions. This study uses fold assessment by NMR in combination with mutagenesis, ITC affinity measurements, and thermal denaturation studies to investigate possible couplings between ASHH2 CW selectivity toward H3K4me1 and the stabilization of the domain and loops implicated in binding. The key elements of the binding site—the two tryptophans and the α1-helix form and maintain the binding pocket— were perturbed by mutagenesis and investigated. Results show that the α1-helix maintains the overall stability of the fold via the I915 and L919 residues and that the correct binding consolidates the loops designated as η1 and η3, as well as the C-terminal. This consolidation is incomplete for H3K4me3 binding to CW, which experiences a decrease in overall thermal stability on binding. Loop mutations not directly involved in the binding site, nonetheless, affect the equilibrium positions of the key residues.

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).

A drug discovery toolbox for Nuclear Magnetic Resonance (NMR) characterization of ligands and their targets

Information about the structure, dynamics, and ligand-binding properties of biomolecules can be derived from Nuclear Magnetic Resonance (NMR) spectroscopy and provides valuable information for drug discovery. A multitude of experimental approaches provides a wealth of information that can be tailored to the system of interest. Methods to study the behavior of ligands upon target binding enable the identification of weak binders in a robust manner that is critical for the identification of truly novel binding interactions. This is particularly important for challenging targets. Observing the solution behavior of biomolecules yields information about their structure, dynamics, and interactions. This review describes the breadth of approaches that are available, many of which are under-utilized in a drug-discovery environment, and focuses on recent advances that continue to emerge.

Conformational exchange of the Zα domain of human RNA editing enzyme ADAR1 studied by NMR spectroscopy

Z-DNA binding proteins (ZBPs) play important roles in RNA editing, innate immune responses, and viral infections. Numerous studies have implicated a role for conformational motions during ZBPs binding upon DNA, but the quantitative intrinsic conformational exchanges of ZBP have not been elucidated. To understand the correlation between the biological function and dynamic feature of the Zα domains of human ADAR1 (hZα_ADAR1), we have performed the ¹⁵N backbone amide Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion experiments on the free hZα_ADAR1 at two different magnetic fields at 35 °C. The robust inter-dependence of parameters in the global fitting process using multi-magnetic field CPMG profiles allows us characterizing the dynamic properties of conformational changes in hZα_ADAR1. This study found that free hZα_ADAR1 exhibited the conformational exchange with a k_ex of 5784 s^-1 between the states "A" (89% population) and "B" (11% population). The different hydrophobic interactions among helices α1, α2, and α3 between these two states might correlate with efficient Z-DNA binding achieved by the hydrogen bonding interactions between its side-chains and the phosphate backbone of Z-DNA.

Hydrogen/Deuterium Exchange and Nuclear Magnetic Resonance Spectroscopy Reveal Dynamic Allostery on Multiple Time Scales in the Serine Protease Thrombin

A deeper understanding of how hydrogen/deuterium exchange mass spectrometry (HDX-MS) reveals allostery is important because HDX-MS can reveal allostery in systems that are not amenable to nuclear magnetic resonance (NMR) spectroscopy. We were able to study thrombin and its complex with thrombomodulin, an allosteric regulator, by both HDX-MS and NMR. In this Perspective, we compare and contrast the results from both experiments and from molecular dynamics simulations. NMR detects changes in the chemical environment around the protein backbone N-H bond vectors, providing residue-level information about the conformational exchange between distinct states. HDX-MS detects changes in amide proton solvent accessibility and H-bonding. Taking advantage of NMR relaxation dispersion measurements of the time scale of motions, we draw conclusions about the motions reflected in HDX-MS experiments. Both experiments detect allostery, but they reveal different components of the allosteric transition. The insights gained from integrating NMR and HDX-MS into thrombin dynamics enable a clearer interpretation of the evidence for allostery revealed by HDX-MS in larger protein complexes and assemblies that are not amenable to NMR.

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.

Systematic Approach to Find the Global Minimum of Relaxation Dispersion Data for Protein-Induced B-Z Transition of DNA

Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion spectroscopy is commonly used for quantifying conformational changes of protein in μs-to-ms timescale transitions. To elucidate the dynamics and mechanism of protein binding, parameters implementing CPMG relaxation dispersion results must be appropriately determined. Building an analytical model for multi-state transitions is particularly complex. In this study, we developed a new global search algorithm that incorporates a random search approach combined with a field-dependent global parameterization method. The robust inter-dependence of the parameters carrying out the global search for individual residues (GSIR) or the global search for total residues (GSTR) provides information on the global minimum of the conformational transition process of the Zα domain of human ADAR1 (hZαADAR1)-DNA complex. The global search results indicated that a α-helical segment of hZαADAR1 provided the main contribution to the three-state conformational changes of a hZαADAR1-DNA complex with a slow B-Z exchange process. The two global exchange rate constants, kex and kZB, were found to be 844 and 9.8 s-1, respectively, in agreement with two regimes of residue-dependent chemical shift differences-the "dominant oscillatory regime" and "semi-oscillatory regime". We anticipate that our global search approach will lead to the development of quantification methods for conformational changes not only in Z-DNA binding protein (ZBP) binding interactions but also in various protein binding processes.

The Picornavirus Precursor 3CD Has Different Conformational Dynamics Compared to 3Cpro and 3Dpol in Functionally Relevant Regions

Viruses have evolved numerous strategies to maximize the use of their limited genetic material, including proteolytic cleavage of polyproteins to yield products with different functions. The poliovirus polyprotein 3CD is involved in important protein-protein, protein-RNA and protein-lipid interactions in viral replication and infection. It is a precursor to the 3C protease and 3D RNA-dependent RNA polymerase, but has different protease specificity, is not an active polymerase, and participates in other interactions differently than its processed products. These functional differences are poorly explained by the known X-ray crystal structures. It has been proposed that functional differences might be due to differences in conformational dynamics between 3C, 3D and 3CD. To address this possibility, we conducted nuclear magnetic resonance spectroscopy experiments, including multiple quantum relaxation dispersion, chemical exchange saturation transfer and methyl spin-spin relaxation, to probe conformational dynamics across multiple timescales. Indeed, these studies identified differences in conformational dynamics in functionally important regions, including enzyme active sites, and RNA and lipid binding sites. Expansion of the conformational ensemble available to 3CD may allow it to perform additional functions not observed in 3C and 3D alone despite having nearly identical lowest-energy structures.

A Perspective on Biophysical Studies of Crystallin Aggregation and Implications for Cataract Formation

Lens crystallins are subject to various types of damage during their lifetime which triggers protein misfolding and aggregation, ultimately causing cataracts. There are several models for crystallin aggregation, but a comprehensive picture of the mechanism of cataract is still underway. The complex biomolecular interactions underlying crystallin aggregation have motivated major efforts to resolve the structural details and mechanism of aggregation using multiple biophysical techniques at different resolutions. Together, experimental and computational approaches identify and characterize both amyloidogenic and amorphous aggregates leading to an improved understanding of crystallin aggregation. A rigorous characterization of the aggregation-prone intermediates is crucial in cataract-mediated drug discovery. This Perspective summarizes recent biophysical studies on lens crystallin aggregation. We evaluate the outstanding challenges, future outlook, and rewards in this fertile field of research. With lessons learned from protein folding and multiple pathways of aggregation, we highlight the differences in the overall mechanisms of age-related and congenital cataracts. We expect that a correlation between the existing and developing biophysical techniques would provide a platform to study amyloid architecture in the eye lens and reduce the existing gaps in our understanding of crystallin biophysics.

Relationship between activity and stability: Design and characterization of stable variants of human frataxin

The relationships between conformational dynamics, stability and protein function are not obvious. Frataxin (FXN) is an essential protein that forms part of a supercomplex dedicated to the iron-sulfur (Fe-S) cluster assembly within the mitochondrial matrix. In humans, the loss of FXN expression or a decrease in its functionality results in Friedreich's Ataxia, a cardio-neurodegenerative disease. Recently, the way in which FXN interacts with the rest of the subunits of the supercomplex was uncovered. This opens a window to explore relationships between structural dynamics and function. In this study, we prepared a set of FXN variants spanning a broad range of conformational stabilities. Variants S160I, S160M and A204R were more stable than the wild-type and showed similar biological activity. Additionally, we prepared SILCAR, a variant that combines S160I, L203C and A204R mutations. SILCAR was 2.4 kcal mol-1 more stable and equally active. Some of the variants were significantly more resistant to proteolysis than the wild-type FXN. SILCAR showed the highest resistance, suggesting a more rigid structure. It was corroborated by means of molecular dynamics simulations. Relaxation dispersion NMR experiments comparing SILCAR and wild-type variants suggested similar internal motions in the microsecond to millisecond timescale. Instead, variant S157I showed higher denaturation resistance but a significant lower function, similarly to that observed for the FRDA variant N146K. We concluded that the contribution of particular side chains to the conformational stability of FXN might be highly subordinated to their impact on both the protein function and the stability of the functional supercomplex.

The Pierced Lasso Topology Leptin has a Bolt on Dynamic Domain Composed by the Disordered Loops I and III

Leptin is an important signaling hormone, mostly known for its role in energy expenditure and satiety. Furthermore, leptin plays a major role in other proteinopathies, such as cancer, marked hyperphagia, impaired immune function, and inflammation. In spite of its biological relevance in human health, there are no NMR resonance assignments of the human protein available, obscuring high-resolution characterization of the soluble protein and/or its conformational dynamics, suggested as being important for receptor interaction and biological activity. Here, we report the nearly complete backbone resonance assignments of human leptin. Chemical shift-based secondary structure prediction confirms that in solution leptin forms a four-helix bundle including a pierced lasso topology. The conformational dynamics, determined on several timescales, show that leptin is monomeric, has a rigid four-helix scaffold, and a dynamic domain, including a transiently formed helix. The dynamic domain is anchored to the helical scaffold by a secondary hydrophobic core, pinning down the long loops of leptin to the protein body, inducing motional restriction without a well-defined secondary or tertiary hydrogen bond stabilized structure. This dynamic region is well suited for and may be involved in functional allosteric dynamics upon receptor binding.

Dissecting the differential structural and dynamics features of CCL2 chemokine orthologs

Chemokines are a sub-group of cytokines that regulate the leukocyte migration. Monocyte chemoattractant protein-1 (MCP/CCL2) is one of the essential CC chemokine that regulates the migration of monocytes into inflamed tissues. It has been observed that the primary sequences of CCL2 orthologs among rodents and primates vary significantly at the C-terminal region. However, no structural details are available for the rodentia family CCL2 proteins. The current study unravelled the structural, dynamics and in-silico functional characteristics of murine CCL2 chemokine using a comprehensive set of NMR spectroscopy techniques and evolutionary approaches. The study unravelled that the N-terminal portion of the murine CCL2 forms a canonical CC chemokine dimer similar to that of human CCL2. However, unlike human CCL2, the murine ortholog exhibits extensive dynamics in the μs-ms timescales. The presence of C-terminal region of the murine CCL2 protein/rodentia family is highly glycosylated, completely disordered, and inhibits the folding of the structured CCL2 regions. Further, it has been observed that the glycosaminoglycan binding surfaces of these orthologs proteins are greatly differed. In a nut shell, this comparative study provided the role of molecular evolution in generating orthologous proteins with differential structural and dynamics characteristics to engage them in specific molecular interactions.

The Arabidopsis (ASHH2) CW domain binds monomethylated K4 of the histone H3 tail through conformational selection

Chromatin post-translational modifications are thought to be important for epigenetic effects on gene expression. Methylation of histone N-terminal tail lysine residues constitutes one of many such modifications, executed by families of histone lysine methyltransferase (HKMTase). One such protein is ASHH2 from the flowering plant Arabidopsis thaliana, equipped with the interaction domain, CW, and the HKMTase domain, SET. The CW domain of ASHH2 is a selective binder of monomethylation at lysine 4 on histone H3 (H3K4me1) and likely helps the enzyme dock correctly onto chromatin sites. The study of CW and related interaction domains has so far been emphasizing lock-key models, missing important aspects of histone-tail CW interactions. We here present an analysis of the ASHH2 CW-H3K4me1 complex using NMR and molecular dynamics, as well as mutation and affinity studies of flexible coils. β-augmentation and rearrangement of coils coincide with changes in the flexibility of the complex, in particular the η1, η3 and C-terminal coils, but also in the β1 and β2 strands and the C-terminal part of the ligand. Furthermore, we show that mutating residues with outlier dynamic behaviour affect the complex binding affinity despite these not being in direct contact with the ligand. Overall, the binding process is consistent with conformational selection. We propose that this binding mechanism presents an advantage when searching for the correct post-translational modification state among the highly modified and flexible histone tails, and also that the binding shifts the catalytic SET domain towards the nucleosome. DATABASES: Structural data are available in the PDB database under the accession code 6QXZ. Resonance assignments for CW42 in its apo- and holo-forms are available in the BMRB database under the accession code 27251.

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.

Chemical Exchange

The phenomenon of chemical or conformational exchange in NMR spectroscopy has enabled detailed characterization of time-dependent aspects of biomolecular function, including folding, molecular recognition, allostery, and catalysis, on timescales from microsecond to second. Importantly, NMR methods based on a variety of spin relaxation parameters have been developed that provide quantitative information on interconversion kinetics, thermodynamic properties, and structural features of molecular states populated to a fraction of a percent at equilibrium and otherwise unobservable by other NMR approaches. The ongoing development of more sophisticated experimental techniques and the necessity to apply these methods to larger and more complex molecular systems engenders a corresponding need for theoretical advances describing such techniques and facilitating data analysis in applications. This review surveys current aspects of the theory of chemical exchange, as utilized in ZZ-exchange; Hahn and Carr-Purcell-Meiboom-Gill (CPMG) spin-echo; and R_1ρ, chemical exchange saturation transfer (CEST), and dark state saturation transfer (DEST) spin-locking experiments. The review emphasizes theoretical results for kinetic topologies with more than two interconverting states, both to obtain compact analytical forms suitable for data analysis and to establish conditions for distinguishability between alternative kinetic schemes.

An Evaluation of the Potential of NMR Spectroscopy and Computational Modelling Methods to Inform Biopharmaceutical Formulations

Protein-based therapeutics are considered to be one of the most important classes of pharmaceuticals on the market. The growing need to prolong stability of high protein concentrations in liquid form has proven to be challenging. Therefore, significant effort is being made to design formulations which can enable the storage of these highly concentrated protein therapies for up to 2 years. Currently, the excipient selection approach involves empirical high-throughput screening, but does not reveal details on aggregation mechanisms or the molecular-level effects of the formulations under storage conditions. Computational modelling approaches have the potential to elucidate such mechanisms, and rapidly screen in silico prior to experimental testing. Nuclear Magnetic Resonance (NMR) spectroscopy can also provide complementary insights into excipient–protein interactions. This review will highlight the underpinning principles of molecular modelling and NMR spectroscopy. It will also discuss the advancements in the applications of computational and NMR approaches in investigating excipient–protein interactions.

Differential Conformational Dynamics Encoded by the Linker between Quasi RNA Recognition Motifs of Heterogeneous Nuclear Ribonucleoprotein H

Members of the heterogeneous nuclear ribonucleoprotein (hnRNP) F/H family are multipurpose RNA binding proteins that participate in most stages of RNA metabolism. Despite having similar RNA sequence preferences, hnRNP F/H proteins function in overlapping and, in some cases, distinct cellular processes. The domain organization of hnRNP F/H proteins is modular, consisting of N-terminal tandem quasi-RNA recognition motifs (F/HqRRM1,2) and a third C-terminal qRRM3 embedded between glycine-rich repeats. The tandem qRRMs are connected through a 10-residue linker, with several amino acids strictly conserved between hnRNP H and F. A significant difference occurs at position 105 of the linker, where hnRNP H contains a proline and hnRNP F an alanine. To investigate the influence of P105 on the conformational properties of hnRNP H, we probed the structural dynamics of its HqRRM1,2 domain with X-ray crystallography, NMR spectroscopy, and small-angle X-ray scattering. The collective results best describe that HqRRM1,2 exists in a conformational equilibrium between compact and extended structures. The compact structure displays an electropositive surface formed at the qRRM1-qRRM2 interface. Comparison of NMR relaxation parameters, including Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion, between HqRRM1,2 and FqRRM1,2 indicates that FqRRM1,2 primarily adopts a more extended and flexible conformation. Introducing the P105A mutation into HqRRM1,2 alters its conformational dynamics to favor an extended structure. Thus, our work demonstrates that the linker compositions confer different structural properties between hnRNP F/H family members that might contribute to their functional diversity.

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.

Magnesium Activates Microsecond Dynamics to Regulate Integrin-Collagen Recognition

Integrin receptors bind collagen via metal-mediated interactions that are modulated by magnesium (Mg²⁺) levels in the extracellular matrix. Nuclear magnetic resonance-based relaxation experiments, isothermal titration calorimetry, and adhesion assays reveal that Mg²⁺ functions as both a structural anchor and dynamic switch of the α₁β₁ integrin I domain (α₁I). Specifically, Mg²⁺ binding activates micro- to millisecond timescale motions of residues distal to the binding site, particularly those surrounding the salt bridge at helix 7 and near the metal ion-dependent adhesion site. Mutagenesis of these residues impacts α₁I functional activity, thereby suggesting that Mg-bound α₁I dynamics are important for collagen binding and consequent allosteric rearrangement of the low-affinity closed to high-affinity open conformation. We propose a multistep recognition mechanism for α₁I-Mg-collagen interactions involving both conformational selection and induced-fit processes. Our findings unravel the multifaceted role of Mg²⁺ in integrin-collagen recognition and assist in elucidating the molecular mechanisms by which metals regulate protein-protein interactions.

Insights on the conformational dynamics of human frataxin through modifications of loop-1

Human frataxin (FXN) is a highly conserved mitochondrial protein involved in iron homeostasis and activation of the iron-sulfur cluster assembly. FXN deficiency causes the neurodegenerative disease Friedreich's Ataxia. Here, we investigated the effect of alterations in loop-1, a stretch presumably essential for FXN function, on the conformational stability and dynamics of the native state. We generated four loop-1 variants, carrying substitutions, insertions and deletions. All of them were stable and well-folded proteins. Fast local motions (ps-ns) and slower long-range conformational dynamics (μs-ms) were altered in some mutants as judged by NMR. Particularly, loop-1 modifications impact on the dynamics of a distant region that includes residues from the β-sheet, helix α1 and the C-terminal. Remarkably, all the mutants retain the ability to activate cysteine desulfurase, even when two of them exhibit a strong decrease in iron binding, revealing a differential sensitivity of these functional features to loop-1 perturbation. Consequently, we found that even for a small and relatively rigid protein, engineering a loop segment enables to alter conformational dynamics through a long-range effect, preserving the native-state structure and important aspects of function.

A 13C-detected 15N double-quantum NMR experiment to probe arginine side-chain guanidinium 15Nη chemical shifts

Arginine side-chains are often key for enzyme catalysis, protein-ligand and protein-protein interactions. The importance of arginine stems from the ability of the terminal guanidinium group to form many key interactions, such as hydrogen bonds and salt bridges, as well as its perpetual positive charge. We present here an arginine ¹³C^ζ-detected NMR experiment in which a double-quantum coherence involving the two ¹⁵N^η nuclei is evolved during the indirect chemical shift evolution period. As the precession frequency of the double-quantum coherence is insensitive to exchange of the two ¹⁵N^η; this new approach is shown to eliminate the previously deleterious line broadenings of ¹⁵N^η resonances caused by the partially restricted rotation about the C^ζ-N^ε bond. Consequently, sharp and well-resolved ¹⁵N^η resonances can be observed. The utility of the presented method is demonstrated on the L99A mutant of the 19 kDa protein T4 lysozyme, where the measurement of small chemical shift perturbations, such as one-bond deuterium isotope shifts, of the arginine amine ¹⁵N^η nuclei becomes possible using the double-quantum experiment.

Solution NMR views of dynamical ordering of biomacromolecules

BACKGROUND

To understand the mechanisms related to the 'dynamical ordering' of macromolecules and biological systems, it is crucial to monitor, in detail, molecular interactions and their dynamics across multiple timescales. Solution nuclear magnetic resonance (NMR) spectroscopy is an ideal tool that can investigate biophysical events at the atomic level, in near-physiological buffer solutions, or even inside cells.

SCOPE OF REVIEW

In the past several decades, progress in solution NMR has significantly contributed to the elucidation of three-dimensional structures, the understanding of conformational motions, and the underlying thermodynamic and kinetic properties of biomacromolecules. This review discusses recent methodological development of NMR, their applications and some of the remaining challenges.

MAJOR CONCLUSIONS

Although a major drawback of NMR is its difficulty in studying the dynamical ordering of larger biomolecular systems, current technologies have achieved considerable success in the structural analysis of substantially large proteins and biomolecular complexes over 1MDa and have characterised a wide range of timescales across which biomolecular motion exists. While NMR is well suited to obtain local structure information in detail, it contributes valuable and unique information within hybrid approaches that combine complementary methodologies, including solution scattering and microscopic techniques.

GENERAL SIGNIFICANCE

For living systems, the dynamic assembly and disassembly of macromolecular complexes is of utmost importance for cellular homeostasis and, if dysregulated, implied in human disease. It is thus instructive for the advancement of the study of the dynamical ordering to discuss the potential possibilities of solution NMR spectroscopy and its applications. This article is part of a Special Issue entitled "Biophysical Exploration of Dynamical Ordering of Biomolecular Systems" edited by Dr. Koichi Kato.

NMR Methods to Study Dynamic Allostery

Nuclear magnetic resonance (NMR) spectroscopy provides a unique toolbox of experimental probes for studying dynamic processes on a wide range of timescales, ranging from picoseconds to milliseconds and beyond. Along with NMR hardware developments, recent methodological advancements have enabled the characterization of allosteric proteins at unprecedented detail, revealing intriguing aspects of allosteric mechanisms and increasing the proportion of the conformational ensemble that can be observed by experiment. Here, we present an overview of NMR spectroscopic methods for characterizing equilibrium fluctuations in free and bound states of allosteric proteins that have been most influential in the field. By combining NMR experimental approaches with molecular simulations, atomistic-level descriptions of the mechanisms by which allosteric phenomena take place are now within reach.