Theory, practice, and applications of paramagnetic relaxation enhancement for the characterization of transient low-population states of biological macromolecules and their complexes

Excellent Fe(II) Binding Tag in Protein Paramagnetic NMR Spectroscopy

Most first series of transition metal ions have one or more unpaired electrons and show great variations in the paramagnetic property. The magnetic anisotropy of some transition metal ions, including Co(II), as well as lanthanide ions [Ln(III)], has been well examined in proteins by NMR. In contrast, few examples of Fe(II) complexes reporting the magnetic anisotropy were analyzed in proteins, except for the ones containing a heme motif or iron-sulfur clusters. Here, we showed that [2,2':6',2″-terpyridine]-6,6″-dicarboxylic acid (TDA) is an excellent iron-binding ligand. It forms a stable iron complex in aqueous solution and demonstrates distinct iron-binding properties. TDA forms a 1:1 stable complex with both Fe(II) and Fe(III), but Fe(II) presents a high-spin state in the complex. The TDA moiety can be site-specifically attached to a protein, and its protein conjugate generates sizable pseudocontact shifts (PCSs) in complex with Fe(II), which are larger than those of commonly used metal binding tags. In contrast, the protein-TDA-Fe(III) complex produces negligible paramagnetic relaxation enhancement (PRE) effects in the protein signals, indicating a low-spin state of Fe(III) in the protein-TDA complex. The high stability of the protein-TDA-Fe(II) complex allows one to measure accurate PCSs in cell lysate even in the presence of other transition metal ions and an excess of GSH.

Dynamic conformations of the P. furiosus MR-DNA complex link Mre11 nuclease activity to DNA-stimulated Rad50 ATP hydrolysis

The MRE11-RAD50-NBS1/Xrs2 (MRN/X) protein complex has essential roles in the repair of damaged DNA. The current understanding of the conformational landscape of the core MR complex comes from various structural studies. However, given the heterogeneous nature of these structures, we suspect that several conformational states may still be unaccounted for. Here, we use methyl-based NMR experiments on P. furiosus MR to determine an ensemble of distinct conformations of MR bound to DNA, consistent with the highly dynamic nature of the MR-DNA complex. Interrogation of these structures via in vitro activity assays on MR mutants reveal an unexpected, striking correlation between the nuclease activity of Mre11 and the magnitude of DNA-stimulated ATP hydrolysis by Rad50. Together, the structures and activity data support a model for MR activity where DNA-stimulated ATP hydrolysis unlocks Rad50 to provide access to the Mre11 active sites and further demonstrate how a heterogeneous ensemble of conformations can be used to coordinate various functions to direct biological outcomes. By elucidating the dynamic conformations of the DNA-bound MR complex, this work lays the foundation for future studies aimed at further characterizing this landscape and dissecting its role in the molecular mechanism of DNA repair and genome stability.

Sizes, conformational fluctuations, and SAXS profiles for intrinsically disordered proteins

The preponderance of intrinsically disordered proteins (IDPs) in the eukaryotic proteome, and their ability to interact with each other, and with folded proteins, RNA, and DNA for functional purposes, have made it important to quantitatively characterize their biophysical properties. Toward this end, we developed the transferable self-organized polymer (SOP-IDP) model to calculate the properties of several IDPs. The values of the radius of gyration (R g $$ {R}_g $$ ) obtained from SOP-IDP simulations are in excellent agreement (correlation coefficient of 0.96) with those estimated from SAXS experiments. For AP180 and Epsin, the predicted values of the hydrodynamic radii (R h s $$ {R}_h\mathrm{s} $$ ) are in nearly quantitative agreement with those from fluorescence correlation spectroscopy (FCS) experiments. Strikingly, the calculated SAXS profiles for 36 IDPs are also nearly superimposable on the experimental profiles. The dependence ofR g $$ {R}_g $$ and the mean end-to-end distance (R ee $$ {R}_{ee} $$ ) on chain length,N $$ N $$ , follows Flory's scaling law,R α ≈ a α N 0.588 $$ {R}_{\alpha}\approx {a}_{\alpha }{N}^{0.588} $$ (α = g , $$ \alpha =g, $$ ande $$ e $$ ), suggesting that globally IDPs behave as synthetic polymers in a good solvent. This finding depends on the solvent quality, which can be altered by changing variables such as pH and salt concentration. The values ofa g $$ {a}_g $$ anda e $$ {a}_e $$ are 0.20 and 0.48 nm, respectively. Surprisingly, finite size corrections to scaling, expected on theoretical grounds, are negligible forR g $$ {R}_g $$ andR ee $$ {R}_{ee} $$ . In contrast, only by accounting for the finite sizes of the IDPs, the dependence of experimentally measurableR h $$ {R}_h $$ onN $$ N $$ can be quantitatively explained usingν = 0.588 $$ \nu =0.588 $$ . Although Flory scaling law captures the estimates forR g $$ {R}_g $$ ,R ee $$ {R}_{ee} $$ , andR h $$ {R}_h $$ accurately, the spread of the simulated data around the theoretical curve is suggestive of of sequence-specific features that emerge through a fine-grained analysis of the conformational ensembles using hierarchical clustering. Typically, the ensemble of conformations partitions into three distinct clusters, having different equilibrium populations and structural properties. Without any further readjustments to the parameters of the SOP-IDP model, we also obtained nearly quantitative agreement with paramagnetic relaxation enhancement (PRE) measurements for α-synuclein. The transferable SOP-IDP model sets the stage for several applications, including the study of phase separation in IDPs and interactions with nucleic acids.

Structure-Based Experimental Datasets for Benchmarking Protein Simulation Force Fields [Article v0.1]

This review article provides an overview of structurally oriented experimental datasets that can be used to benchmark protein force fields, focusing on data generated by nuclear magnetic resonance (NMR) spectroscopy and room temperature (RT) protein crystallography. We discuss what the observables are, what they tell us about structure and dynamics, what makes them useful for assessing force field accuracy, and how they can be connected to molecular dynamics simulations carried out using the force field one wishes to benchmark. We also touch on statistical issues that arise when comparing simulations with experiment. We hope this article will be particularly useful to computational researchers and trainees who develop, benchmark, or use protein force fields for molecular simulations.

Probing Biomolecular Interactions with Paramagnetic Nuclear Magnetic Resonance Spectroscopy

Recent advances in computational methods like AlphaFold have transformed structural biology, enabling accurate modeling of protein complexes and driving applications in drug discovery and protein engineering. However, predicting the structure of systems involving weak, transient, or dynamic interactions, or of complexes with disordered regions, remains challenging. Nuclear Magnetic Resonance (NMR) spectroscopy offers atomic-level insights into biomolecular complexes, even in weakly interacting and dynamic systems. Paramagnetic NMR, in particular, provides long-range structural restraints, easily exceeding distances over 25 Å, making it ideal for studying large protein complexes. Advances in chemical tools for introducing paramagnetic tags into proteins, combined with progress in electron paramagnetic resonance (EPR) spectroscopy, have enhanced the method's utility. This perspective article discusses paramagnetic NMR approaches for analyzing biomolecular complexes in solution and in the solid state, emphasizing quantities like pseudocontact shifts, residual dipolar couplings, and paramagnetic relaxation enhancements. Additionally, dynamic nuclear polarization offers a promising method to amplify NMR signals of large complexes, even in complex environments. The integration of AlphaFold protein structure prediction with paramagnetic NMR holds great potential for advancing our understanding of biomolecular interactions.

Terminal spin labeling of xylotriose strongly affects interactions in the active site of xylanase BcX

Paramagnetic probes provide long-range distance information and report on minor conformations of biomacromolecules. However, it is important to realize that any probe can affect the system of interest. Here, we report on the effects of attaching a small nitroxide spin label [TEMPO, (2,2,6,6-tetramethylpiperidin-1-yl)oxyl] to xylotriose, a substrate of the enzyme xylanase from Bacillus circulans (BcX). BcX has a long and narrow active site cleft accommodating six xylose units and a secondary binding site on its surface. The aim of the study was to probe the interactions of the substrate with the enzyme using paramagnetic relaxation enhancements (PREs). Binding of the substrate to the surface exposed secondary binding site resulted in strong and localized PREs, indicative of well-defined binding. The xylotriose with diamagnetic control tag was still able to bind the active site cleft, though the rate of exchange was reduced relative to that of untagged xylotriose. The substrate with the paramagnetic TEMPO was not able to bind inside the active site cleft. Also, additional interactions on another surface location showed differences between the paramagnetic substrate and the diamagnetic control, despite the minimal chemical differences between TEMPO modified xylotriose and its reduced, diamagnetic counterpart. Our findings underscore the sensitivity of BcX substrate binding to minor substrate modifications. This study serves as a reminder that any probe, including the attachment of a small paramagnetic group, can affect the behavior of the system under investigation. Even the chemical difference between a paramagnetic tag and its diamagnetic control can result in differences in the molecular interactions.

Quantitative Interpretation of Transverse Spin Relaxation by Translational Diffusion in Liquids

Intermolecular spin relaxation by translational diffusion of spin pairs has been widely used to study the properties of biomolecules in liquids. Notably, solvent paramagnetic relaxation enhancement (sPRE) arising from paramagnetic cosolutes has gained significant attention for various applications in structural biology, including the structural refinement of intrinsically disordered proteins, the elucidation of the molecular mechanisms driving cosolute-induced protein denaturation, and the characterization of residue-specific effective near-surface electrostatic potentials (ENS). Furthermore, sPRE has been extensively applied in magnetic resonance imaging (MRI), where paramagnetic ions, such as Gd(III)-based ions, are used as contrast agents. Among these applications, the transverse sPRE rate (Γ2) has predominantly been interpreted empirically as being proportional to the average interspin distance ⟨r-6⟩norm. In this study, we present a rigorous theoretical interpretation of Γ2 for spherically symmetric intermolecular potentials, demonstrating that it is proportional to ⟨r-4⟩norm. We provide an explicit formula for calculating ⟨r-4⟩norm without any adjustable parameters, offering valuable insights into the interaction potential independent of the type or strength of interactions. It has broad applicability, including the precise interpretation of the relaxation properties of the MRI contrast agents and calculation of the ENS.

Biophysical investigation of the molecular interaction between minichromosome maintenance protein 6 and Bloom syndrome helicase

The minichromosome maintenance protein (MCM) complex and Bloom syndrome helicase (BLM) are crucial components in DNA replication and cell division. MCM, a hexameric helicase that unwinds double-stranded DNA, serves as an important diagnostic and prognostic biomarker for cancer cells and a target for anticancer drug development. BLM, associated with G-quadruplex structures, is another key helicase in maintaining genomic stability. In this study, we investigate the interaction between MCM6 and BLM at the atomic level, as their expression levels are highly correlated in various cancer types, with elevated levels indicating poor prognosis. To elucidate the molecular basis of MCM6/BLM interaction, we employed fluorescence polarization anisotropy analysis, NMR chemical shifts perturbation analysis (CSP), and paramagnetic relaxation enhancement (PRE) experiments. MCM6 binding domain (MBD) C and D exhibit similar binding affinities to MCM6 winged-helix domain (WHD). However, significant CSPs with MBD-D and PRE experiments suggested that MBD-D is closer to MCM6 WHD than MBD-C. Despite both proteins containing numerous negatively charged residues, hydrophobic interactions govern the association between MCM6 WHD and BLM MBD-D. This biophysical characterization of the MCM6/BLM interaction provides new insights into their functional relationship and challenges existing models. Our findings reveal that MCM6 binds BLM at a different site than its other known partner chromatin licensing and DNA replication factor. Understanding these protein-protein interactions at the molecular level may contribute to the development of novel anticancer therapies targeting the MCM6/BLM interaction.

Interaction of Cyclosporin C with Dy3+ Ions in Acetonitrile and in Complex with Dodecylphosphocholine Micelles Determined by NMR Spectroscopy

The spectral characteristics of cyclosporin C (CsC) with the addition of Dy3+ ions in acetonitrile (CD3CN) and CsC with Dy3+ incorporated into dodecylphosphocholine (DPC) micelle in deuterated water were investigated by high-resolution NMR spectroscopy. The study was focused on the interaction between Dy3+ ions and CsC molecules in different environments. Using a combination of one-dimensional and two-dimensional NMR techniques, we obtained information on the spatial features of the peptide molecule and the interaction between CsC and the metal ion. The non-uniform effect of the metal ion on different NMR signals of CsC was observed. The paramagnetic attenuation parameter was calculated for the amide, alpha, and beta protons of CsC upon the addition of Dy3+. The metal ion was found to interact with the polar part of the DPC micelle, and the ion also has a significant effect on the NMR signals of amino acid residues from Sar3 to d-Ala8. This pattern is reproduced in both environments studied here and also agrees with earlier investigations of the CsA-Dy3+ complex.

Molecular insights into the interaction between a disordered protein and a folded RNA

Intrinsically disordered protein regions (IDRs) are well established as contributors to intermolecular interactions and the formation of biomolecular condensates. In particular, RNA-binding proteins (RBPs) often harbor IDRs in addition to folded RNA-binding domains that contribute to RBP function. To understand the dynamic interactions of an IDR-RNA complex, we characterized the RNA-binding features of a small (68 residues), positively charged IDR-containing protein, Small ERDK-Rich Factor (SERF). At high concentrations, SERF and RNA undergo charge-driven associative phase separation to form a protein- and RNA-rich dense phase. A key advantage of this model system is that this threshold for demixing is sufficiently high that we could use solution-state biophysical methods to interrogate the stoichiometric complexes of SERF with RNA in the one-phase regime. Herein, we describe our comprehensive characterization of SERF alone and in complex with a small fragment of the HIV-1 Trans-Activation Response (TAR) RNA with complementary biophysical methods and molecular simulations. We find that this binding event is not accompanied by the acquisition of structure by either molecule; however, we see evidence for a modest global compaction of the SERF ensemble when bound to RNA. This behavior likely reflects attenuated charge repulsion within SERF via binding to the polyanionic RNA and provides a rationale for the higher-order assembly of SERF in the context of RNA. We envision that the SERF-RNA system will lower the barrier to accessing the details that support IDR-RNA interactions and likewise deepen our understanding of the role of IDR-RNA contacts in complex formation and liquid-liquid phase separation.

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.

Likely Overstabilization of Charge-Charge Interactions in CHARMM36m(w): A Case for a99SB-disp Water

Recent years have witnessed drastic improvements in general-purpose explicit solvent protein force fields, partially driven by the need to study intrinsically disordered proteins (IDPs), and yet the state-of-the-art force fields such as CHARMM36m (c36m) and a99SB-disp still provide different performances in simulating disordered protein states, where c36m has a bias toward overcompaction for large IDPs. Here, we examine the performance of c36m and a99SB-disp in describing the stabilities of a set of 46 amino acid backbone and side chain pairs in various configurations. The free energy results show that c36m systematically predicts stronger interactions compared to a99SB-disp by an average of 0.2 kcal/mol for nonpolar pairs, 0.6 kcal/mol for polar pairs, and 0.8 kcal/mol for salt bridges. The most severe overstabilization in c36m is observed for charged pairs involving the Arg and Glu side chains by up to 2.9 kcal/mol. Importantly, the systematic overstabilization of c36m is only marginally alleviated by c36mw, an ad hoc patch to c36m that increases the dispersion interactions between TIP3P hydrogens and protein atoms. Guided by free energy decomposition, we evaluated if revising the charges alone could alleviate the severe overstabilization of salt bridges of c36m(w) vs a99SB-disp. The results suggested that the direct modification of protein-water interactions is also necessary. Toward this end, we proposed a tentative modification to c36m, referred to as c36mrb-disp, which combines modified Arg side chain charges, retuned backbone hydrogen bonding strength, and the a99SB-disp water model. The modified force field successfully reproduces the secondary structures of several intrinsically disordered peptides and proteins including (AAQAA)3, GB1p, and p53 transactivation domain, while maintaining the stability of a set of folded proteins. This work provides a set of useful systems for benchmarking and optimizing protein force fields and highlights the importance of balancing protein-protein and protein-water electrostatic interactions for accurately describing both folded and disordered proteins.

Multivalent interactions of the disordered regions of XLF and XRCC4 foster robust cellular NHEJ and drive the formation of ligation-boosting condensates in vitro

In mammalian cells, DNA double-strand breaks are predominantly repaired by non-homologous end joining (NHEJ). During repair, the Ku70-Ku80 heterodimer (Ku), X-ray repair cross complementing 4 (XRCC4) in complex with DNA ligase 4 (X4L4) and XRCC4-like factor (XLF) form a flexible scaffold that holds the broken DNA ends together. Insights into the architectural organization of the NHEJ scaffold and its regulation by the DNA-dependent protein kinase catalytic subunit (DNA-PKcs) were recently obtained by single-particle cryo-electron microscopy analysis. However, several regions, especially the C-terminal regions (CTRs) of the XRCC4 and XLF scaffolding proteins, have largely remained unresolved in experimental structures, which hampers the understanding of their functions. Here we used magnetic resonance techniques and biochemical assays to comprehensively characterize the interactions and dynamics of the XRCC4 and XLF CTRs at residue resolution. We show that the CTRs of XRCC4 and XLF are intrinsically disordered and form a network of multivalent heterotypic and homotypic interactions that promotes robust cellular NHEJ activity. Importantly, we demonstrate that the multivalent interactions of these CTRs lead to the formation of XLF and X4L4 condensates in vitro, which can recruit relevant effectors and critically stimulate DNA end ligation. Our work highlights the role of disordered regions in the mechanism and dynamics of NHEJ and lays the groundwork for the investigation of NHEJ protein disorder and its associated condensates inside cells with implications in cancer biology, immunology and the development of genome-editing strategies.

HEXIM1 homodimer binds two sites on 7SK RNA to release autoinhibition for P-TEFb inactivation

Hexim proteins are RNA-dependent regulators whose main target is 7SK long non-coding RNA, a major regulator of eukaryotic mRNA transcription. 7SK RNPs control available intracellular concentrations of the kinase P-TEFb (Cdk9-CyclinT1/2) by sequestering it in an inactive form. Active P-TEFb phosphorylates NELF, DSIF, and the RNA polymerase II CTD to transition it from promoter-proximal pausing to productive elongation. P-TEFb associates with 7SK RNP via Hexim, which directly binds 7SK RNA. However, free Hexim is in an autoinhibited state that cannot inactivate P-TEFb, and how Hexim autoinhibition is released by 7SK remains unknown. Here, we show that one Hexim1 homodimer binds two sites on linear 7SK RNA in a manner that exposes the Cdk9 binding sites, which are otherwise masked within the autoinhibited dimer. These results provide mechanistic insights into Hexim-RNA specificity and explain how P-TEFb can be effectively regulated to respond to changing levels of transcriptional signaling.

Human CSTF2 RNA Recognition Motif Domain Binds to a U-Rich RNA Sequence through a Multistep Binding Process

The RNA recognition motif (RRM) is a conserved and ubiquitous RNA-binding domain that plays essential roles in mRNA splicing, polyadenylation, transport, and stability. RRM domains exhibit remarkable diversity in binding partners, interacting with various sequences of single- and double-stranded RNA, despite their small size and compact fold. During pre-mRNA cleavage and polyadenylation, the RRM domain from CSTF2 recognizes U- or G/U-rich RNA sequences downstream from the cleavage and polyadenylation site to regulate the process. Given the importance of alternative cleavage and polyadenylation in increasing the diversity of mRNAs, the exact mechanism of binding of RNA to the RRM of CSTF2 remains unclear, particularly in the absence of a structure of this RRM bound to a native RNA substrate. Here, we performed a series of NMR titration and spin relaxation experiments, which were complemented by paramagnetic relaxation enhancement measurements and rigid-body docking, to characterize the interactions of the CSTF2 RRM with a U-rich ligand. Our results reveal a multistep binding process involving differences in ps-ns time scale dynamics and potential structural changes, particularly in the C-terminalα-helix. These results provide insights into how the CSTF2 RRM domain binds to U-rich RNA ligands and offer a greater understanding for the molecular basis of the regulation of pre-mRNA cleavage and polyadenylation.

Utilizing Molecular Dynamics Simulations, Machine Learning, Cryo-EM, and NMR Spectroscopy to Predict and Validate Protein Dynamics

Protein dynamics play a crucial role in biological function, encompassing motions ranging from atomic vibrations to large-scale conformational changes. Recent advancements in experimental techniques, computational methods, and artificial intelligence have revolutionized our understanding of protein dynamics. Nuclear magnetic resonance spectroscopy provides atomic-resolution insights, while molecular dynamics simulations offer detailed trajectories of protein motions. Computational methods applied to X-ray crystallography and cryo-electron microscopy (cryo-EM) have enabled the exploration of protein dynamics, capturing conformational ensembles that were previously unattainable. The integration of machine learning, exemplified by AlphaFold2, has accelerated structure prediction and dynamics analysis. These approaches have revealed the importance of protein dynamics in allosteric regulation, enzyme catalysis, and intrinsically disordered proteins. The shift towards ensemble representations of protein structures and the application of single-molecule techniques have further enhanced our ability to capture the dynamic nature of proteins. Understanding protein dynamics is essential for elucidating biological mechanisms, designing drugs, and developing novel biocatalysts, marking a significant paradigm shift in structural biology and drug discovery.

Gadolinium-Based NMR Spin Relaxation Measurements of Near-Surface Electrostatic Potentials of Biomolecules

NMR spectroscopy is an important tool for the measurement of the electrostatic properties of biomolecules. To this point, paramagnetic relaxation enhancements (PREs) of 1H nuclei arising from nitroxide cosolutes in biomolecular solutions have been used to measure effective near-surface electrostatic potentials (ϕENS) of proteins and nucleic acids. Here, we present a gadolinium (Gd)-based NMR method, exploiting Gd chelates with different net charges, for measuring ϕENS values and demonstrate its utility through applications to a number of biomolecular systems. The use of Gd-based cosolutes offers several advantages over nitroxides for ϕENS measurements. First, unlike nitroxide compounds, Gd chelates enable electrostatic potential measurements on oxidation-sensitive proteins that require reducing agents. Second, the large electron spin quantum number of Gd (7/2) results in notably larger PREs for Gd chelates when used at the same concentrations as nitroxide radicals. Thus, it is possible to measure ϕENS values exclusively from + and - charged compounds even for highly charged biomolecules, avoiding the use of neutral cosolutes that, as we further establish here, limits the accuracy of the measured electrostatic potentials. In addition, the smaller concentrations of cosolutes required minimize potential binding to sites on macromolecules. Fourth, the closer proximity of the paramagnetic center and charged groups within Gd chelates, in comparison to the corresponding nitroxide compounds, enables more accurate predictions of ϕENS potentials for cross-validation of the experimental results. The Gd-based method described here, thus, broadens the applicability of studies of biomolecular electrostatics using solution NMR spectroscopy.

Structural details of helix-mediated TDP-43 C-terminal domain multimerization

The primarily disordered C-terminal domain (CTD) of TAR DNA binding protein-43 (TDP-43), a key nuclear protein in RNA metabolism, forms neuronal inclusions in several neurodegenerative diseases. A conserved region (CR, spanning residues 319-341) in CTD forms transient helix-helix contacts important for its higher-order oligomerization and function that are disrupted by ALS-associated mutations. However, the structural details of CR assembly and the explanation for several ALS-associated variants' impact on phase separation and function remain unclear due to challenges in analyzing the dynamic association of TDP-43 CTD using traditional structural biology approaches. By employing an integrative approach, combining biophysical experiments, biochemical assays, AlphaFold2-Multimer (AF2-Multimer), and atomistic simulations, we generated structural models of helical oligomerization of TDP-43 CR. Using NMR, we first established that the native state of TDP-43 CR under physiological conditions is α-helical. Next, alanine scanning mutagenesis revealed that while hydrophobic residues in the CR are important for CR assembly, phase separation and TDP-43 nuclear retention function, polar residues down regulate these processes. Finally, pairing AF2-Multimer modeling with AAMD simulations indicated that dynamic, oligomeric assemblies of TDP-43 that are stabilized by a methionine-rich core with specific contributions from a tryptophan/leucine pair. In conclusion, our results advance the structural understanding of the mechanisms driving TDP-43 function and provide a window into the initial stages of its conversion into pathogenic aggregates.

Resolution in cryogenic solid state NMR: Challenges and solutions

NMR at cryogenic temperatures has the potential to provide rich site-specific details regarding biopolymer structure, function, and mechanistic intermediates. Broad spectral lines compared with room temperature NMR can sometimes present practical challenges. A number of hypotheses regarding the origins of line broadening are explored. One frequently considered explanation is the presence of inhomogeneous conformational distributions. Possibly these arise when the facile characteristic motions that occur near room temperature become dramatically slower or "frozen out" at temperatures below the solvent phase change. Recent studies of low temperature spectra harness the distributions in properties in these low temperature spectra to uncover information regarding the conformational ensembles that drive biological function.

Molecular insights into the interaction between a disordered protein and a folded RNA

Intrinsically disordered protein regions (IDRs) are well-established as contributors to intermolecular interactions and the formation of biomolecular condensates. In particular, RNA-binding proteins (RBPs) often harbor IDRs in addition to folded RNA-binding domains that contribute to RBP function. To understand the dynamic interactions of an IDR-RNA complex, we characterized the RNA-binding features of a small (68 residues), positively charged IDR-containing protein, SERF. At high concentrations, SERF and RNA undergo charge-driven associative phase separation to form a protein- and RNA-rich dense phase. A key advantage of this model system is that this threshold for demixing is sufficiently high that we could use solution-state biophysical methods to interrogate the stoichiometric complexes of SERF with RNA in the one-phase regime. Herein, we describe our comprehensive characterization of SERF alone and in complex with a small fragment of the HIV-1 TAR RNA (TAR) with complementary biophysical methods and molecular simulations. We find that this binding event is not accompanied by the acquisition of structure by either molecule; however, we see evidence for a modest global compaction of the SERF ensemble when bound to RNA. This behavior likely reflects attenuated charge repulsion within SERF via binding to the polyanionic RNA and provides a rationale for the higher-order assembly of SERF in the context of RNA. We envision that the SERF-RNA system will lower the barrier to accessing the details that support IDR-RNA interactions and likewise deepen our understanding of the role of IDR-RNA contacts in complex formation and liquid-liquid phase separation.

Structural basis for the regulation of plant transcription factor WRKY33 by the VQ protein SIB1

The WRKY transcription factors play essential roles in a variety of plant signaling pathways associated with biotic and abiotic stress response. The transcriptional activity of many WRKY members are regulated by a class of intrinsically disordered VQ proteins. While it is known that VQ proteins interact with the WRKY DNA-binding domains (DBDs), also termed as the WRKY domains, structural information regarding VQ-WRKY interaction is lacking and the regulation mechanism remains unknown. Herein we report a solution NMR study of the interaction between Arabidopsis WRKY33 and its regulatory VQ protein partner SIB1. We uncover a SIB1 minimal sequence neccessary for forming a stable complex with WRKY33 DBD, which comprises not only the consensus "FxxhVQxhTG" VQ motif but also its preceding region. We demonstrate that the βN-strand and the extended βN-β1 loop of WRKY33 DBD form the SIB1 docking site, and build a structural model of the complex based on the NMR paramagnetic relaxation enhancement and mutagenesis data. Based on this model, we further identify a cluster of positively-charged residues in the N-terminal region of SIB1 to be essential for the formation of a SIB1-WRKY33-DNA ternary complex. These results provide a framework for the mechanism of SIB1-enhanced WRKY33 transcriptional activity.

Probing protein aggregation through spectroscopic insights and multimodal approaches: A comprehensive review for counteracting neurodegenerative disorders

Aberrant accumulation of protein misfolding can cause aggregation and fibrillation and is one of the primary characteristic features of neurodegenerative diseases. Because they are disordered, misfolded, and aggregated proteins pose a significant setback in drug designing. The structural study of intermediate steps in these kinds of aggregated proteins will allow us to determine the conformational changes as well as the probable pathways encompassing various neurodegenerative disorders. The analysis of protein aggregates involved in neurodegenerative diseases relies on a diverse toolkit of biophysical techniques, encompassing both morphological and non-morphological methods. Additionally, Thioflavin T (ThT) assays and Circular Dichroism (CD) spectroscopy facilitate investigations into aggregation kinetics and secondary structure alterations. The collective application of these biophysical techniques empowers researchers to comprehensively unravel the intricate nature of protein aggregates associated with neurodegeneration. Furthermore, the topics covered in this review have summed up a handful of well-established techniques used for the structural analysis of protein aggregation. This multifaceted approach advances our fundamental understanding of the underlying mechanisms driving neurodegenerative diseases and informs potential therapeutic strategies.

Experimental methods to study the structure and dynamics of intrinsically disordered regions in proteins

Eukaryotic proteins often feature long stretches of amino acids that lack a well-defined three-dimensional structure and are referred to as intrinsically disordered proteins (IDPs) or regions (IDRs). Although these proteins challenge conventional structure-function paradigms, they play vital roles in cellular processes. Recent progress in experimental techniques, such as NMR spectroscopy, single molecule FRET, high speed AFM and SAXS, have provided valuable insights into the biophysical basis of IDP function. This review discusses the advancements made in these techniques particularly for the study of disordered regions in proteins. In NMR spectroscopy new strategies such as 13C detection, non-uniform sampling, segmental isotope labeling, and rapid data acquisition methods address the challenges posed by spectral overcrowding and low stability of IDPs. The importance of various NMR parameters, including chemical shifts, hydrogen exchange rates, and relaxation measurements, to reveal transient secondary structures within IDRs and IDPs are presented. Given the high flexibility of IDPs, the review outlines NMR methods for assessing their dynamics at both fast (ps-ns) and slow (μs-ms) timescales. IDPs exert their functions through interactions with other molecules such as proteins, DNA, or RNA. NMR-based titration experiments yield insights into the thermodynamics and kinetics of these interactions. Detailed study of IDPs requires multiple experimental techniques, and thus, several methods are described for studying disordered proteins, highlighting their respective advantages and limitations. The potential for integrating these complementary techniques, each offering unique perspectives, is explored to achieve a comprehensive understanding of IDPs.

Analyzing paramagnetic NMR data on target DNA search by proteins using a discrete-state kinetic model for translocation

Before reaching their targets, sequence-specific DNA-binding proteins nonspecifically bind to DNA through electrostatic interactions and stochastically change their locations on DNA. Investigations into the dynamics of DNA-scanning by proteins are nontrivial due to the simultaneous presence of multiple translocation mechanisms and many sites for the protein to nonspecifically bind to DNA. Nuclear magnetic resonance (NMR) spectroscopy can provide information about the target DNA search processes at an atomic level. Paramagnetic relaxation enhancement (PRE) is particularly useful to study how the proteins scan DNA in the search process. Previously, relatively simple two-state or three-state exchange models were used to explain PRE data reflecting the target search process. In this work, using more realistic discrete-state stochastic kinetics models embedded into an NMR master equation, we analyzed the PRE data for the HoxD9 homeodomain interacting with DNA. The kinetic models that incorporate sliding, dissociation, association, and intersegment transfer can reproduce the PRE profiles observed at some different ionic strengths. The analysis confirms the previous interpretation of the PRE data and shows that the protein's probability distribution among nonspecific sites is nonuniform during the target DNA search process.

Membrane-Driven Dimerization of the Peripheral Membrane Protein KRAS: Implications for Downstream Signaling

Transient homo-dimerization of the RAS GTPase at the plasma membrane has been shown to promote the mitogen-activated protein kinase (MAPK) signaling pathway essential for cell proliferation and oncogenesis. To date, numerous crystallographic studies have focused on the well-defined GTPase domains of RAS isoforms, which lack the disordered C-terminal membrane anchor, thus providing limited structural insight into membrane-bound RAS molecules. Recently, lipid-bilayer nanodisc platforms and paramagnetic relaxation enhancement (PRE) analyses have revealed several distinct structures of the membrane-anchored homodimers of KRAS, an isoform that is most frequently mutated in human cancers. The KRAS dimerization interface is highly plastic and altered by biologically relevant conditions, including oncogenic mutations, the nucleotide states of the protein, and the lipid composition. Notably, PRE-derived structures of KRAS homodimers on the membrane substantially differ in terms of the relative orientation of the protomers at an "α-α" dimer interface comprising two α4-α5 regions. This interface plasticity along with the altered orientations of KRAS on the membrane impact the accessibility of KRAS to downstream effectors and regulatory proteins. Further, nanodisc platforms used to drive KRAS dimerization can be used to screen potential anticancer drugs that target membrane-bound RAS dimers and probe their structural mechanism of action.

Paramagnetic NMR to study iron sulfur proteins: 13C detected experiments illuminate the vicinity of the metal center

The robustness of NMR coherence transfer in proximity of a paramagnetic center depends on the relaxation properties of the nuclei involved. In the case of Iron-Sulfur Proteins, different pulse schemes or different parameter sets often provide complementary results. Tailored versions of HCACO and CACO experiments significantly increase the number of observed Cα/C' connectivities in highly paramagnetic systems, by recovering many resonances that were lost due to paramagnetic relaxation. Optimized 13C direct detected experiments can significantly extend the available assignments, improving the overall knowledge of these systems. The different relaxation properties of Cα and C' nuclei are exploited in CACO vs COCA experiments and the complementarity of the two experiments is used to obtain structural information. The two [Fe2S2]+ clusters containing NEET protein CISD3 and the one [Fe4S4]2+ cluster containing HiPIP protein PioC have been taken as model systems. We show that tailored experiments contribute to decrease the blind sphere around the cluster, to extend resonance assignment of cluster bound cysteine residues and to retrieve details on the topology of the iron-bound ligand residues.

Localization of the ibuprofen molecule in model lipid membranes revealed by spin-label-enhanced NMR relaxation

Non-steroidal anti-inflammatory drugs (NSAIDs) have antipyretic, anti-inflammatory and analgesic effects, and can be used in the treatment of various diseases. These drugs have also a number of side effects, which may be related to their interaction with lipid membranes. In this study, we use the spin-labeled NSAID ibuprofen (ibuprofen-SL) as a relaxation enhancer to study its interaction with model lipid membranes employing liquid-state 1H NMR at 500 MHz. The high magnetic moment of unpaired electron in the spin label made it possible to reduce the concentration of the studied drug in the membrane to tenths of a mole percent. As model membranes, unilamellar POPC liposomes and bicelles consisting of a 2:1 mixture of DHPC:DMPC or DHPC:POPC lipids were used. An increase in the rate of proton spin-lattice relaxation, T1-1, selectively detected for protons at different positions in the lipid molecule, showed that ibuprofen-SL is localized in the hydrophobic part of the lipid bilayer. As the concentration of ibuprofen-SL increases to 0.5 mol%, the distribution of positions of ibuprofen-SL across the bilayer becomes wider. In the presence of 20 mol% of cholesterol, ibuprofen-SL is displaced from the core of the membrane to a region closer to the head group of the bilayer. This displacement was also confirmed by the NMR NOESY experiment conducted with unlabeled ibuprofen. For bilayers containing unsaturated POPC lipids, the distribution of ibuprofen positions across the bilayer becomes narrower compared to the presence of saturated DMPC lipids.

Unexpected dynamics in femtomolar complexes of binding proteins with peptides

Ultra-tight binding is usually observed for proteins associating with rigidified molecules. Previously, we demonstrated that femtomolar binders derived from the Armadillo repeat proteins (ArmRPs) can be designed to interact very tightly with fully flexible peptides. Here we show for ArmRPs with four and seven sequence-identical internal repeats that the peptide-ArmRP complexes display conformational dynamics. These dynamics stem from transient breakages of individual protein-residue contacts that are unrelated to overall unbinding. The labile contacts involve electrostatic interactions. We speculate that these dynamics allow attaining very high binding affinities, since they reduce entropic losses. Importantly, only NMR techniques can pick up these local events by directly detecting conformational exchange processes without complications from changes in solvent entropy. Furthermore, we demonstrate that the interaction surface of the repeat protein regularizes upon peptide binding to become more compatible with the peptide geometry. These results provide novel design principles for ultra-tight binders.

Observation of pH-Dependent Residual Structure in the Pmel17 Repeat Domain and the Implication for Its Amyloid Formation

The varying conformational states of amyloid-forming protein monomers can determine their fibrillation outcome. In this study, we utilize solution NMR and the paramagnetic relaxation enhancement (PRE) effect to observe monomer properties of the repeat domain (RPT) from a human functional amyloid, premelanosomal protein, Pmel17. After excision from the full-length protein, RPT can self-assemble into amyloid fibrils, functioning as a scaffold for melanin deposition. Here, we report possible conformational states of the short RPT (sRPT) isoform, which has been demonstrated to be a fibrillation nucleator. NMR experiments were performed to determine conformational differences in sRPT by comparing aggregation-prone vs nonaggregating solution conditions. We observed significant chemical shift perturbations localized to residues near the C-terminus, demonstrating that the local chemical environment of the amyloid core region is highly sensitive to changes in pH. Next, we introduced cysteine point mutations for the covalent attachment of PRE ligands to sRPT to facilitate the observation of intramolecular interactions. We also utilized solvent PRE molecules with opposing charges to measure changes in the electrostatic potential of sRPT in different pH environments. These observed PRE effects offer insight into initial molecular events that might promote intermolecular interactions, which can trigger fibrillation. Taken together, our results show that sRPT monomers adopt a conformation inconsistent with a fully random coil at neutral pH and undergo conformational changes at lower pH values. These observations highlight regulatory mechanisms via organelle-associated pH conditions that can affect the fibrillation activity of proteins like RPT.

The internal structure of gadolinium and perfluorocarbon-loaded polymer nanoparticles affects 19F MRI relaxation times

19F magnetic resonance imaging (19F MRI) is an emerging technique for quantitative imaging in novel therapies, such as cellular therapies and theranostic nanocarriers. Nanocarriers loaded with liquid perfluorocarbon (PFC) typically have a (single) core-shell structure with PFC in the core due to the poor miscibility of PFC with organic and inorganic solvents. Paramagnetic relaxation enhancement acts only at a distance of a few angstroms. Thus, efficient modulation of the 19F signal is possible only with fluorophilic PFC-soluble chelates. However, these chelates cannot interact with the surrounding environment and they might result in image artifacts. Conversely, chelates bound to the nanoparticle shell typically have a minimal effect on the 19F signal and a strong impact on the aqueous environment. We show that the confinement of PFC in biodegradable polymeric nanoparticles (NPs) with a multicore structure enables the modulation of longitudinal (T1) and transverse (T2) 19F relaxation, as well as proton (1H) signals, using non-fluorophilic paramagnetic chelates. We compared multicore NPs versus a conventional single core structure, where the PFC is encapsulated in the core(s) and the chelate in the surrounding polymeric matrix. This modulated relaxation also makes multicore NPs sensitive to various acidic pH environments, while preserving their stability. This effect was not observed with single core nanocapsules (NCs). Importantly, paramagnetic chelates affected both T1 and T219F relaxation in multicore NPs, but not in single core NCs. Both relaxation times of the 19F nucleus were enhanced with an increasing concentration of the paramagnetic chelate. Moreover, as the polymeric matrix remained water permeable, proton enhancement additionally was observed in MRI.

Directed evolution unlocks oxygen reactivity for a nicotine-degrading flavoenzyme

The flavoenzyme nicotine oxidoreductase (NicA2) is a promising injectable treatment to aid in the cessation of smoking, a behavior responsible for one in ten deaths worldwide. NicA2 acts by degrading nicotine in the bloodstream before it reaches the brain. Clinical use of NicA2 is limited by its poor catalytic activity in the absence of its natural electron acceptor CycN. Without CycN, NicA2 is instead oxidized slowly by dioxygen (O2), necessitating unfeasibly large doses in a therapeutic setting. Here, we report a genetic selection strategy that directly links CycN-independent activity of NicA2 to growth of Pseudomonas putida S16. This selection enabled us to evolve NicA2 variants with substantial improvement in their rate of oxidation by O2. The encoded mutations cluster around a putative O2 tunnel, increasing flexibility and accessibility to O2 in this region. These mutations further confer desirable clinical properties. A variant form of NicA2 is tenfold more effective than the wild type at degrading nicotine in the bloodstream of rats.

Developments in small-angle X-ray scattering (SAXS) for characterizing the structure of surfactant-macromolecule interactions and their complex

The surfactant-macromolecule interactions (SMI) are one of the most critical topics for scientific research and industrial application. Small-angle X-ray scattering (SAXS) is a powerful tool for comprehensively studying the structural and conformational features of macromolecules at a size ranging from Angstroms to hundreds of nanometers with a time-resolve in milliseconds scale. The SAXS integrative techniques have emerged for comprehensively analyzing the SMI and the structure of their complex in solution. Here, the various types of emerging interactions of surfactant with macromolecules, such as protein, lipid, nuclear acid, polysaccharide and virus, etc. have been systematically reviewed. Additionally, the principle of SAXS and theoretical models of SAXS for describing the structure of SMI as well as their complex has been summarized. Moreover, the recent developments in the applications of SAXS for charactering the structure of SMI have been also highlighted. Prospectively, the capacity to complement artificial intelligence (AI) in the structure prediction of biological macromolecules and the high-throughput bioinformatics sequencing data make SAXS integrative structural techniques expected to be the primary methodology for illuminating the self-assembling dynamics and nanoscale structure of SMI. As advances in the field continue, we look forward to proliferating uses of SAXS based upon its abilities to robustly produce mechanistic insights for biology and medicine.

Nuclear magnetic resonance/single molecule fluorescence combinations to study dynamic protein systems

Many proteins require different structural states or conformations for function, and intrinsically disordered proteins, i.e. proteins without stable three-dimensional structure, are certainly an extreme. Single molecule fluorescence and nuclear magnetic resonance (NMR) spectroscopy are both exceptionally well suited to decipher and describe these states and their interconversion. Different time scales, from picoseconds to several milliseconds, can be addressed by both techniques. The length scales probed and the sample requirements (e.g. concentration, molecular weight, sample complexity) are, however, vastly different, making NMR and single molecule fluorescence an excellent combination for integrated studies. Here, we review recently undertaken approaches for the combined use of NMR and single molecule fluorescence to study protein dynamics.

Direct measurements of biomolecular electrostatics through experiments

Biomolecular electrostatics has been a subject of computational investigations based on 3D structures. This situation is changing because emerging experimental tools allow us to quantitatively investigate biomolecular electrostatics without any use of structure information. Now, electrostatic potentials around biomolecules can directly be measured for many residues simultaneously by nuclear magnetic resonance (NMR) spectroscopy. This NMR method can be used to study electrostatic aspects of various processes, including macromolecular association and liquid-liquid phase separation. Applications to structurally flexible biomolecules such as intrinsically disordered proteins are particularly useful. The new tools also facilitate examination of theoretical models and methods for biomolecular electrostatics.

19F electron nuclear double resonance (ENDOR) spectroscopy for distance measurements using trityl spin labels in DNA duplexes

The combination of fluorine labeling and pulsed electron-nuclear double resonance (ENDOR) is emerging as a powerful technique for obtaining structural information about proteins and nucleic acids. In this work, we explored the capability of Mims 19F ENDOR experiments on reporting intermolecular distances in trityl- and 19F-labeled DNA duplexes at three electron paramagnetic resonance (EPR) frequencies (34, 94, and 263 GHz). For spin labeling, we used the hydrophobic Finland trityl radical and hydrophilic OX063 trityl radical. Fluorine labels were introduced into two positions of a DNA oligonucleotide. The results indicated that hyperfine splittings visible in the ENDOR spectra are consistent with the most populated interspin distances between 19F and the trityl radical predicted from molecular dynamic (MD) simulations. Moreover, for some cases, ENDOR spectral simulations based on MD results were able to reproduce the conformational distribution reflected in the experimental ENDOR line broadening. Additionally, MD simulations provided more detailed information about the melting of terminal base pairs of the oligonucleotides and about the configuration of the trityls relative to a DNA end.

The translational repressor Glorund uses interchangeable RNA recognition domains to recognize Drosophila nanos

The Drosophila melanogaster protein Glorund (Glo) represses nanos (nos) translation and uses its quasi-RNA recognition motifs (qRRMs) to recognize both G-tract and structured UA-rich motifs within the nos translational control element (TCE). We showed previously that each of the three qRRMs is multifunctional, capable of binding to G-tract and UA-rich motifs, yet if and how the qRRMs combine to recognize the nos TCE remained unclear. Here we determined solution structures of a nos TCEI_III RNA containing the G-tract and UA-rich motifs. The RNA structure demonstrated that a single qRRM is physically incapable of recognizing both RNA elements simultaneously. In vivo experiments further indicated that any two qRRMs are sufficient to repress nos translation. We probed interactions of Glo qRRMs with TCEI_III RNA using NMR paramagnetic relaxation experiments. Our in vitro and in vivo data support a model whereby tandem Glo qRRMs are indeed multifunctional and interchangeable for recognition of TCE G-tract or UA-rich motifs. This study illustrates how multiple RNA recognition modules within an RNA-binding protein may combine to diversify the RNAs that are recognized and regulated.

Lipid acyl chain protrusion induced by the influenza virus hemagglutinin fusion peptide detected by NMR paramagnetic relaxation enhancement

The glycoprotein spikes of membrane-enveloped viruses include a subunit that catalyzes fusion (joining) of the viral and target cell membranes. For influenza virus, this is subunit 2 of hemagglutinin which has a ∼ 20-residue N-terminal fusion peptide (Fp) region that binds target membrane. An outstanding question is whether there are associated membrane changes important for fusion. Several computational studies have found increased "protrusion" of lipid acyl chains near Fp, i.e. one or more chain carbons are closer to the aqueous region than the headgroup phosphorus. Protrusion may accelerate initial joining of outer leaflets of the two membranes into a stalk intermediate. In this study, higher protrusion probability in membrane with vs. without Fp is convincingly detected by larger Mn2+-associated increases in chain 13C NMR transverse relaxation rates (Γ2's). Data analysis provides a ratio Γ2,neighbor/Γ2,distant for lipids neighboring vs. more distant from the Fp. The calculated ratio depends on the number of Fp-neighboring lipids and the experimentally-derived range of 4 to 24 matches the range of increased protrusion probabilities from different simulations. For samples either with or without Fp, the Γ2 values are well-fitted by an exponential decay as the 13C site moves closer to the chain terminus. The decays correlate with free-energy of protrusion proportional to the number of protruded -CH2 groups, with free energy per -CH2 of ∼0.25 kBT. The NMR data support one major fusion role of the Fp to be much greater protrusion of lipid chains, with highest protrusion probability for chain regions closest to the headgroups.

Native Top-Down Mass Spectrometry Uncovers Two Distinct Binding Motifs of a Functional Neomycin-Sensing Riboswitch Aptamer

Understanding how ligands bind to ribonucleic acids (RNA) is important for understanding RNA recognition in biological processes and drug development. Here, we have studied neomycin B binding to neomycin-sensing riboswitch aptamer constructs by native top-down mass spectrometry (MS) using electrospray ionization (ESI) and collisionally activated dissociation (CAD). Our MS data for a 27 nt aptamer construct reveal the binding site and ligand interactions, in excellent agreement with the structure derived from nuclear magnetic resonance (NMR) studies. Strikingly, for an extended 40 nt aptamer construct, which represents the sequence with the highest regulatory factor for riboswitch function, we identified two binding motifs for neomycin B binding, one corresponding to the bulge-loop motif of the 27 nt construct and the other one in the minor groove of the lower stem, which according to the MS data are equally populated. By replacing a noncanonical with a canonical base pair in the lower stem of the 40 nt aptamer, we can reduce binding to the minor groove motif from ∼50 to ∼30%. Conversely, the introduction of a CUG/CUG motif in the lower stem shifts the binding equilibrium in favor of minor groove binding. The MS data reveal site-specific and stoichiometry-resolved information on aminoglycoside binding to RNA that is not directly accessible by other methods and underscore the role of noncanonical base pairs in RNA recognition by aminoglycosides.

Multivalent interactions of the disordered regions of XLF and XRCC4 foster robust cellular NHEJ and drive the formation of ligation-boosting condensates in vitro

In mammalian cells, DNA double-strand breaks are predominantly repaired by non-homologous end joining (NHEJ). During repair, the Ku70/80 heterodimer (Ku), XRCC4 in complex with DNA Ligase 4 (X4L4), and XLF form a flexible scaffold that holds the broken DNA ends together. Insights into the architectural organization of the NHEJ scaffold and its regulation by the DNA-dependent protein kinase catalytic subunit (DNA-PKcs) have recently been obtained by single-particle cryo-electron microscopy analysis. However, several regions, especially the C-terminal regions (CTRs) of the XRCC4 and XLF scaffolding proteins, have largely remained unresolved in experimental structures, which hampers the understanding of their functions. Here, we used magnetic resonance techniques and biochemical assays to comprehensively characterize the interactions and dynamics of the XRCC4 and XLF CTRs at atomic resolution. We show that the CTRs of XRCC4 and XLF are intrinsically disordered and form a network of multivalent heterotypic and homotypic interactions that promotes robust cellular NHEJ activity. Importantly, we demonstrate that the multivalent interactions of these CTRs led to the formation of XLF and X4L4 condensates in vitro which can recruit relevant effectors and critically stimulate DNA end ligation. Our work highlights the role of disordered regions in the mechanism and dynamics of NHEJ and lays the groundwork for the investigation of NHEJ protein disorder and its associated condensates inside cells with implications in cancer biology, immunology and the development of genome editing strategies.

Selective polypeptide ligand binding to the extracellular surface of the transmembrane domains of the class B GPCRs GLP-1R and GCGR

Crystal and cryo-EM structures of the glucagon-like peptide-1 receptor (GLP-1R) and glucagon receptor (GCGR) bound with their peptide ligands have been obtained with full-length constructs, indicating that the extracellular domain (ECD) is indispensable for specific ligand binding. This article complements these data with studies of ligand recognition of the two receptors in solution. Paramagnetic NMR relaxation enhancement measurements using dual labeling with fluorine-19 probes on the receptor and nitroxide spin labels on the peptide ligands provided new insights. The glucagon-like peptide-1 (GLP-1) was found to interact with GLP-1R by selective binding to the extracellular surface. The ligand selectivity toward the extracellular surface of the receptor was preserved in the transmembrane domain (TMD) devoid of the ECD. The dual labeling approach further provided evidence of cross-reactivity of GLP-1R and GCGR with glucagon and GLP-1, respectively, which is of interest in the context of medical treatments using combinations of the two polypeptides.

Conformational tuning improves the stability of spirocyclic nitroxides with long paramagnetic relaxation times

Nitroxides are widely used as probes and polarization transfer agents in spectroscopy and imaging. These applications require high stability towards reducing biological environments, as well as beneficial relaxation properties. While the latter is provided by spirocyclic groups on the nitroxide scaffold, such systems are not in themselves robust under reducing conditions. In this work, we introduce a strategy for stability enhancement through conformational tuning, where incorporating additional substituents on the nitroxide ring effects a shift towards highly stable closed spirocyclic conformations, as indicated by X-ray crystallography and density functional theory (DFT) calculations. Closed spirocyclohexyl nitroxides exhibit dramatically improved stability towards reduction by ascorbate, while maintaining long relaxation times in electron paramagnetic resonance (EPR) spectroscopy. These findings have important implications for the future design of new nitroxide-based spin labels and imaging agents.

Quantitative Paramagnetic NMR-Based Analysis of Protein Orientational Dynamics on Membranes: Dissecting the KRas4B-Membrane Interactions

Peripheral membrane proteins can adopt distinct orientations on the surfaces of lipid bilayers that are often short-lived and challenging to characterize by conventional experimental methods. Here we describe a robust approach for mapping protein orientational landscapes through quantitative interpretation of paramagnetic relaxation enhancement (PRE) data arising from membrane mimetics with spin-labeled lipids. Theoretical analysis, followed by experimental verification, reveals insights into the distinct properties of the PRE observables that are generally distorted in the case of stably membrane-anchored proteins. To suppress the artifacts, we demonstrate that undistorted Γ₂ values can be obtained via transient membrane anchoring, based on which a computational framework is established for deriving accurate orientational ensembles obeying Boltzmann statistics. Application of the approach to KRas4B, a classical peripheral membrane protein whose orientations are critical for its functions and drug design, reveals four distinct orientational states that are close but not identical to those reported previously. Similar orientations are also found for a truncated KRas4B without the hypervariable region (HVR) that can sample a broader range of orientations, suggesting a confinement role of the HVR geometrically prohibiting severe tilting. Comparison of the KRas4B Γ₂ rates measured using nanodiscs containing different types of anionic lipids reveals identical Γ₂ patterns for the G-domain but different ones for the HVR, indicating only the latter is able to selectively interact with anionic lipids.

Recent Advances in NMR Protein Structure Prediction with ROSETTA

Nuclear magnetic resonance (NMR) spectroscopy is a powerful method for studying the structure and dynamics of proteins in their native state. For high-resolution NMR structure determination, the collection of a rich restraint dataset is necessary. This can be difficult to achieve for proteins with high molecular weight or a complex architecture. Computational modeling techniques can complement sparse NMR datasets (<1 restraint per residue) with additional structural information to elucidate protein structures in these difficult cases. The Rosetta software for protein structure modeling and design is used by structural biologists for structure determination tasks in which limited experimental data is available. This review gives an overview of the computational protocols available in the Rosetta framework for modeling protein structures from NMR data. We explain the computational algorithms used for the integration of different NMR data types in Rosetta. We also highlight new developments, including modeling tools for data from paramagnetic NMR and hydrogen-deuterium exchange, as well as chemical shifts in CS-Rosetta. Furthermore, strategies are discussed to complement and improve structure predictions made by the current state-of-the-art AlphaFold2 program using NMR-guided Rosetta modeling.

Function and dynamics of the intrinsically disordered carboxyl terminus of β2 adrenergic receptor

Advances in structural biology have provided important mechanistic insights into signaling by the transmembrane core of G-protein coupled receptors (GPCRs); however, much less is known about intrinsically disordered regions such as the carboxyl terminus (CT), which is highly flexible and not visible in GPCR structures. The β2 adrenergic receptor's (β2AR) 71 amino acid CT is a substrate for GPCR kinases and binds β-arrestins to regulate signaling. Here we show that the β2AR CT directly inhibits basal and agonist-stimulated signaling in cell lines lacking β-arrestins. Combining single-molecule fluorescence resonance energy transfer (FRET), NMR spectroscopy, and molecular dynamics simulations, we reveal that the negatively charged β2AR-CT serves as an autoinhibitory factor via interacting with the positively charged cytoplasmic surface of the receptor to limit access to G-proteins. The stability of this interaction is influenced by agonists and allosteric modulators, emphasizing that the CT plays important role in allosterically regulating GPCR activation.

Interpreting Transient Interactions of Intrinsically Disordered Proteins

The flexible nature of intrinsically disordered proteins (IDPs) gives rise to a conformational ensemble with a diverse set of conformations. The simplest way to describe this ensemble is through a homopolymer model without any specific interactions. However, there has been growing evidence that the conformational properties of IDPs and their relevant functions can be affected by transient interactions between specific and even nonlocal pairs of amino acids. Interpreting these interactions from experimental methods, each of which is most sensitive to a different distance regime referred to as probing length, remains a challenging and unsolved problem. Here, we first show that transient interactions can be realized between short fragments of charged amino acids by generating conformational ensembles using model disordered peptides and coarse-grained simulations. Using these ensembles, we investigate how sensitive different types of experimental measurements are to the presence of transient interactions. We find methods with shorter probing lengths to be more appropriate for detecting these transient interactions, but one experimental method is not sufficient due to the existence of other weak interactions typically seen in IDPs. Finally, we develop an adjusted polymer model with an additional short-distance peak which can robustly reproduce the distance distribution function from two experimental measurements with complementary short and long probing lengths. This new model can suggest whether a homopolymer model is insufficient for describing a specific IDP and meets the challenge of quantitatively identifying specific, transient interactions from a background of nonspecific, weak interactions.

Peptides that Mimic RS repeats modulate phase separation of SRSF1, revealing a reliance on combined stacking and electrostatic interactions

Phase separation plays crucial roles in both sustaining cellular function and perpetuating disease states. Despite extensive studies, our understanding of this process is hindered by low solubility of phase-separating proteins. One example of this is found in SR and SR-related proteins. These proteins are characterized by domains rich in arginine and serine (RS domains), which are essential to alternative splicing and in vivo phase separation. However, they are also responsible for a low solubility that has made these proteins difficult to study for decades. Here, we solubilize the founding member of the SR family, SRSF1, by introducing a peptide mimicking RS repeats as a co-solute. We find that this RS-mimic peptide forms interactions similar to those of the protein's RS domain. Both interact with a combination of surface-exposed aromatic residues and acidic residues on SRSF1's RNA Recognition Motifs (RRMs) through electrostatic and cation-pi interactions. Analysis of RRM domains from human SR proteins indicates that these sites are conserved across the protein family. In addition to opening an avenue to previously unavailable proteins, our work provides insight into how SR proteins phase separate and participate in nuclear speckles.

Peptidase Activation by a Leader Peptide-Bound RiPP Recognition Element

The RiPP precursor recognition element (RRE) is a conserved domain found in many prokaryotic ribosomally synthesized and post-translationally modified peptide (RiPP) biosynthetic gene clusters (BGCs). RREs bind with high specificity and affinity to a recognition sequence within the N-terminal leader region of RiPP precursor peptides. Lasso peptide biosynthesis involves an RRE-dependent leader peptidase, which is discretely encoded or fused to the RRE as a di-domain protein. Here we leveraged thousands of predicted BGCs to define the RRE:leader peptidase interaction through evolutionary covariance analysis. Each interacting domain contributes a three-stranded β-sheet to form a hydrophobic β-sandwich-like interface. The bioinformatics-guided predictions were experimentally confirmed using proteins from discrete and fused lasso peptide BGC architectures. Support for the domain-domain interface derived from chemical shift perturbation, paramagnetic relaxation enhancement experiments, and rapid variant activity screening using cell-free biosynthesis. Further validation of selected variants was performed with purified proteins. We developed a p-nitroanilide-based leader peptidase assay to illuminate the role of RRE domains. Our data show that RRE domains play a dual function. RRE domains deliver the precursor peptide to the leader peptidase, and the rate is saturable as expected for a substrate. RRE domains also partially compose the elusive S2 proteolytic pocket that binds the penultimate threonine of lasso leader peptides. Because the RRE domain is required to form the active site, leader peptidase activity is greatly diminished when the RRE domain is supplied at substoichiometric levels. Full proteolytic activation requires RRE engagement with the recognition sequence-containing portion of the leader peptide. Together, our observations define a new mechanism for protease activity regulation.

Interaction of the bacterial division regulator MinE with lipid bicelles studied by NMR spectroscopy

The bacterial MinE and MinD division regulatory proteins form a standing wave enabling MinC, which binds MinD, to inhibit FtsZ polymerization everywhere except at the midcell, thereby assuring correct positioning of the cytokinetic septum and even distribution of contents to daughter cells. The MinE dimer undergoes major structural rearrangements between a resting six-stranded state present in the cytoplasm, a membrane-bound state, and a four-stranded active state bound to MinD on the membrane, but it is unclear which MinE motifs interact with the membrane in these different states. Using NMR, we probe the structure and global dynamics of MinE bound to disc-shaped lipid bicelles. In the bicelle-bound state, helix α1 no longer sits on top of the six-stranded β-sheet, losing any contact with the protein core, but interacts directly with the bicelle surface; the structure of the protein core remains unperturbed and also interacts with the bicelle surface via helix α2. Binding may involve a previously identified excited state of free MinE in which helix α1 is disordered, thereby allowing it to target the membrane surface. Helix α1 and the protein core undergo nanosecond rigid body motions of differing amplitudes in the plane of the bicelle surface. Global dynamics on the sub-millisecond time scale between a ground state and a sparsely populated excited state are also observed and may represent a very early intermediate on the transition path between the resting six-stranded and active four-stranded conformations. In summary, our results provide insights into MinE structural rearrangements important during bacterial cell division.

Autoinhibitory structure of preligand association state implicates a new strategy to attain effective DR5 receptor activation

Members of the tumor necrosis factor receptor superfamily (TNFRSF) are important therapeutic targets that can be activated to induce death of cancer cells or stimulate proliferation of immune cells. Although it has long been implicated that these receptors assemble preligand associated states that are required for dominant interference in human disease, such states have so far eluded structural characterization. Here, we find that the ectodomain of death receptor 5 (DR5-ECD), a representative member of TNFRSF, can specifically self-associate when anchored to lipid bilayer, and we report this self-association structure determined by nuclear magnetic resonance (NMR). Unexpectedly, two non-overlapping interaction interfaces are identified that could propagate to higher-order clusters. Structure-guided mutagenesis indicates that the observed preligand association structure is represented on DR5-expressing cells. The DR5 preligand association serves an autoinhibitory role as single-domain antibodies (sdAbs) that partially dissociate the preligand cluster can sensitize the receptor to its ligand TRAIL and even induce substantial receptor signaling in the absence of TRAIL. Unlike most agonistic antibodies that require multivalent binding to aggregate receptors for activation, these agonistic sdAbs are monovalent and act specifically on an oligomeric, autoinhibitory configuration of the receptor. Our data indicate that receptors such as DR5 can form structurally defined preclusters incompatible with signaling and that true agonists should disrupt the preligand cluster while converting it to signaling-productive cluster. This mechanism enhances our understanding of a long-standing question in TNFRSF signaling and suggests a new opportunity for developing agonistic molecules by targeting receptor preligand clustering.

Cosolvent Dimethyl Sulfoxide Influences Protein-Ligand Binding Kinetics via Solvent Viscosity Effects: Revealing the Success Rate of Complex Formation Following Diffusive Protein-Ligand Encounter

Protein-ligand-exchange kinetics determines the duration of biochemical signals and consequently plays an important role in drug design. Binding studies commonly require solubilization of designed ligands in solvents such as dimethyl sulfoxide (DMSO), resulting in residual amounts of DMSO following titration of solubilized ligands into aqueous protein samples. Therefore, it is critical to establish whether DMSO influences protein-ligand binding. Here, we address the general and indirect effect of DMSO on protein-ligand binding caused by solvent viscosity, which is strongly dependent on the relative concentrations of DMSO and water. As a model system, we studied the binding of a drug-like ligand to the carbohydrate recognition domain of galectin-3 in the presence of variable amounts of DMSO. We used isothermal titration calorimetry to characterize binding thermodynamics and ¹⁵N NMR relaxation to monitor kinetics. The binding enthalpy is not affected, but we observe a subtle trend of increasingly unfavorable entropy of binding, and consequently decreased affinity, with increasing DMSO concentration. The increasing concentration of DMSO results in a reduced association rate of binding, while the dissociation rate is less affected. The observed association rate is inversely proportional to the viscosity of the DMSO-water mixture, as expected from theory, but significantly reduced from the diffusion-controlled limit. By comparing the viscosity dependence of the observed association rate with that of the theoretical diffusion-controlled association rate, we estimate the success rate of productive complex formation following an initial encounter of proteins and ligands, showing that only one out of several hundred binding "attempts" are successful.