Nmrglue: an open source Python package for the analysis of multidimensional NMR data

Photodynamic treatment increases the lifespan and oxidative stress resistance of Caenorhabditis elegans

Photodynamic therapy is a noninvasive treatment in which specific photosensitizers and light are used to produce high amounts of reactive oxygen species (ROS), which can be employed for targeted tissue destruction in cancer treatment or antimicrobial therapy. However, it remains unknown whether lower amounts of ROS produced by mild photodynamic therapy increase lifespan and stress resistance at the organism level. Here, we introduce a novel photodynamic treatment (PDTr) that uses 20 μM hypericin, a photosensitizer that originates from Hypericum perforatum, and orange light (590 nm, 5.4 W/m2, 1 min) to induce intracellular ROS formation (ROS), thereby resulting in lifespan extension and improved stress resistance in C. elegans. The PDTr-induced increase in longevity was abrogated by N-acetyl cysteine, suggesting the hormetic response was driven by prooxidative mechanisms. PDTr activated the translocation of SKN-1/NRF-2 and DAF-16/FOXO, leading to elevated expression of downstream oxidative stress-responsive genes, including ctl-1, gst-4, and sod-3. In summary, our findings suggest a novel PDTr method that extends the lifespan of C. elegans under both normal and oxidative stress conditions through the activation of SKN-1 and DAF-16 via the involvement of many antioxidant genes.

Structure and dynamics of the proton-selective histidine and the gating tryptophan in an inward rectifying hybrid influenza B and A virus M2 proton channel

The M2 proteins of influenza A and B viruses form acid-activated proton channels that are essential for the virus lifecycle. Proton selectivity is achieved by a transmembrane (TM) histidine whereas gating is achieved by a tryptophan residue. Although this functional apparatus is conserved between AM2 and BM2 channels, AM2 conducts protons exclusively inward whereas BM2 conducts protons in either direction depending on the pH gradient. Previous studies showed that in AM2, mutations of D44 abolished inward rectification of AM2, suggesting that the tryptophan gate is destabilized. To elucidate how charged residues C-terminal to the tryptophan regulates channel gating, here we investigate the structure and dynamics of H19 and W23 in a BM2 mutant, GDR-BM2, in which three BM2 residues are mutated to the corresponding AM2 residues, S16G, G26D and H27R. Whole-cell electrophysiological data show that GDR-BM2 conducts protons with inward rectification, identical to wild-type (WT) AM2 but different from WT-BM2. Solid-state NMR 15N and 13C spectra of H19 indicate that the mutant BM2 channel contains higher populations of cationic histidine and neutral τ tautomers compared to WT-BM2 at acidic pH. Moreover, 19F NMR spectra of 5-19F-labeled W23 resolve three peaks at acidic pH, suggesting three tryptophan sidechain conformations. Comparison of these spectra with the tryptophan spectra of other M2 peptides suggests that these indole sidechain conformations arise from interactions with the C-terminal charged residues and with the N-terminal cationic histidine. Taken together, these solid-state NMR data show that inward rectification in M2 proton channels is accomplished by tryptophan interactions with charged residues on both its C-terminal and N-terminal sides. Gating of these M2 proton channels is thus accomplished by a multi-residue complex with finely tuned electrostatic and aromatic interactions.

Tracking Metabolite Variations during the Degradation of Vegetables in Rice Bran Bed with Intact-State Nuclear Magnetic Resonance Spectroscopy

Fermentation-a process of compound degradation by microorganisms-is a traditional food processing method utilized worldwide for the long-term preservation of fresh foods. In recent years, fermented foods have gained attention as health foods. Fermentation increases the nutritional value of ingredients, producing complex flavors and aromas. To identify unknown components in fermented foods, it is necessary to analyze compounds and conditions nondestructively and comprehensively. We performed intact-state nuclear magnetic resonance (NMR) spectroscopy using intermolecular single quantum coherence (iSQC) to detect the degradation of vegetables directly and nondestructively. We used two types of vegetables and a rice bran bed (nukazuke), which is used for traditional vegetable fermentation in Japan. Major metabolites such as saccharides, organic acids, and amino acids were identified in iSQC-sliced spectra. Comparing NMR signal intensities during degradation revealed the transition of metabolites characteristic of lactic acid fermentation. A pathway-based network analysis showed pathways involved in amino acid metabolism and lactic acid fermentation. Our analytical approach with intact-state NMR spectroscopy using iSQC demonstrated that it may be effective in other experimental systems, allowing for the evaluation of phenomena that have been conventionally overlooked in their true state.

Pathological mutations in the phospholamban cytoplasmic region affect its topology and dynamics modulating the extent of SERCA inhibition

Phospholamban (PLN) is a 52 amino acid regulin that allosterically modulates the activity of the sarco(endo)plasmic reticulum Ca2+-ATPase (SERCA) in the heart muscle. In its unphosphorylated form, PLN binds SERCA within its transmembrane (TM) domains, approximately 20 Å away from the Ca2+ binding site, reducing SERCA's apparent Ca2+ affinity (pKCa) and decreasing cardiac contractility. During the enzymatic cycle, the inhibitory TM domain of PLN remains anchored to SERCA, whereas its cytoplasmic region transiently binds the ATPase's headpiece. Phosphorylation of PLN at Ser16 by protein kinase A increases the affinity of its cytoplasmic domain to SERCA, weakening the TM interactions with the ATPase, reversing its inhibitory function, and augmenting muscle contractility. How the structural changes caused by pathological mutations in the PLN cytoplasmic region are transmitted to its inhibitory TM domain is still unclear. Using solid-state NMR spectroscopy and activity assays, we analyzed the structural and functional effects of a series of mutations and their phosphorylated forms located in the PLN cytoplasmic region and linked to dilated cardiomyopathy. We found that these missense mutations affect the overall topology and dynamics of PLN and ultimately modulate its inhibitory potency. Also, the changes in the TM tilt angle and cytoplasmic dynamics of PLN caused by these mutations correlate well with the extent of SERCA inhibition. Our study unveils new molecular determinants for designing variants of PLN that outcompete endogenous PLN to regulate SERCA in a tunable manner.

Supporting the assignment of NMR spectra with variable-temperature experiments

Nuclear magnetic resonance (NMR) spectroscopy is one of the most powerful tools in analytical chemistry. An important step in the analysis of NMR data is the assignment of resonance frequencies to the corresponding atoms in the molecule being investigated. The traditional approach considers the spectrum's characteristic parameters: chemical shift values, internuclear couplings, and peak intensities. In this paper, we show how to support the process of assigning a series of spectra of similar organic compounds by using temperature coefficients, that is, the rates of change in chemical shift values associated with given changes in temperature.

Improved analysis of NMR chemical shift perturbations through an error estimation method

In solution NMR, chemical shift perturbation (CSP) experiments are widely employed to study intermolecular interactions. However, excluding the nonsignificant peak shift is difficult because little is known about errors in CSP. Here, to address this issue, we introduce a method for estimating errors in CSP based on the noise level. First, we developed a technique that involves line shape fitting to estimate errors in peak position via Monte Carlo simulations. Second, this technique was applied to estimate errors in CSP. In intermolecular interaction analysis of VAP-A with SNX2, error estimation of CSP enabled the evaluation of small but significant changes in peak position and yielded detailed insights that are unattainable with conventional CSP analysis. Third, this technique was successfully applied to estimate errors in residual dipolar couplings. In conclusion, our error estimation method improves CSP analysis by excluding the nonsignificant peak shift.

Quantitative analysis of the slow exchange process by 19F NMR in the presence of scalar and dipolar couplings: applications to the ribose 2'-19F probe in nucleic acids

Solution NMR spectroscopy is a particularly powerful technique for characterizing the functional dynamics of biomolecules, which is typically achieved through the quantitative characterization of chemical exchange processes via the measurement of spin relaxation rates. In addition to the conventional nuclei such as 15N and 13C, which are abundant in biomolecules, fluorine-19 (19F) has recently garnered attention and is being widely used as a site-specific spin probe. While 19F offers the advantages of high sensitivity and low background, it can be susceptible to artifacts in quantitative relaxation analyses due to a multitude of dipolar and scalar coupling interactions with nearby 1H spins. In this study, we focused on the ribose 2'-19F spin probe in nucleic acids and investigated the effects of 1H-19F spin interactions on the quantitative characterization of slow exchange processes on the millisecond time scale. We demonstrated that the 1H-19F dipolar coupling can significantly affect the interpretation of 19F chemical exchange saturation transfer (CEST) experiments when 1H decoupling is applied, while the 1H-19F interactions have a lesser impact on Carr-Purcell-Meiboom-Gill relaxation dispersion applications. We also proposed a modified CEST scheme to alleviate these artifacts along with experimental verifications on self-complementary RNA systems. The theoretical framework presented in this study can be widely applied to various 19F spin systems where 1H-19F interactions are operative, further expanding the utility of 19F relaxation-based NMR experiments.

Solution-state methyl NMR spectroscopy of large non-deuterated proteins enabled by deep neural networks

Methyl-TROSY nuclear magnetic resonance (NMR) spectroscopy is a powerful technique for characterising large biomolecules in solution. However, preparing samples for these experiments is demanding and entails deuteration, limiting its use. Here we demonstrate that NMR spectra recorded on protonated, uniformly 13C labelled samples can be processed using deep neural networks to yield spectra that are of similar quality to typical deuterated methyl-TROSY spectra, potentially providing information for proteins that cannot be produced in bacterial systems. We validate the methodology experimentally on three proteins with molecular weights in the range 42-360 kDa. We further demonstrate the applicability of our methodology to 3D NOESY spectra of Escherichia coli Malate Synthase G (81 kDa), where observed NOE cross-peaks are in good agreement with the available structure. The method represents an advance in the field of using deep learning to analyse complex magnetic resonance data and could have an impact on the study of large biomolecules in years to come.

Improved detection of magnetic interactions in proteins based on long-lived coherences

Living systems rely on molecular building blocks with low structural symmetry. Therefore, constituent amino acids and nucleotides yield short-lived nuclear magnetic responses to electromagnetic radiation. Magnetic signals are at the basis of molecular imaging, structure determination and interaction studies. In solution state, as the molecular weight of analytes increases, coherences with long lifetimes are needed to yield advantageous through-space magnetisation transfers. Interactions between magnetic nuclei can only be detected provided the lifetimes of spin order are sufficient. In J-coupled pairs of nuclei, long-lived coherences (LLC's) connect states with different spin-permutation symmetry. Here in, we show sustained LLC's in protein Lysozyme, weighing 14.3 kDa, with lifetimes twice as long as those of classical magnetisation for the aliphatic protons of glycine residues. We found for the first time that, in a protein of significant molecular weight, LLC's yield substantial through-space magnetisation transfers: spin-order transfer stemming from LLC's overcame transfers from classical coherences by factors > 2. Furthermore, in agreement with theory, the permutation symmetry of LLC-based transfers allows mapping interacting atoms in the protein structure with respect to the molecular plane of glycine residues in a stereospecific manner. These findings can extend the scope of liquid-state high-resolution biomolecular spectroscopy.

Streamlining the automated discovery of porous organic cages

Self-assembly through dynamic covalent chemistry (DCC) can yield a range of multi-component organic assemblies. The reversibility and dynamic nature of DCC has made prediction of reaction outcome particularly difficult and thus slows the discovery rate of new organic materials. In addition, traditional experimental processes are time-consuming and often rely on serendipity. Here, we present a streamlined hybrid workflow that combines automated high-throughput experimentation, automated data analysis, and computational modelling, to accelerate the discovery process of one particular subclass of molecular organic materials, porous organic cages. We demonstrate how the design and implementation of this workflow aids in the identification of organic cages with desirable properties. The curation of a precursor library of 55 tri- and di-topic aldehyde and amine precursors enabled the experimental screening of 366 imine condensation reactions experimentally, and 1464 hypothetical organic cage outcomes to be computationally modelled. From the screen, 225 cages were identified experimentally using mass spectrometry, 54 of which were cleanly formed as a single topology as determined by both turbidity measurements and 1H NMR spectroscopy. Integration of these characterisation methods into a fully automated Python pipeline, named cagey, led to over a 350-fold decrease in the time required for data analysis. This work highlights the advantages of combining automated synthesis, characterisation, and analysis, for large-scale data curation towards an accessible data-driven materials discovery approach.

NMR characterization of RNA binding property of the DEAD-box RNA helicase DDX3X and its implications for helicase activity

The DEAD-box RNA helicase (DDX) plays a central role in many aspects of RNA metabolism by remodeling the defined structure of RNA molecules. While a number of structural studies have revealed the atomistic details of the interaction between DDX and RNA ligands, the molecular mechanism of how this molecule unwinds a structured RNA into an unstructured single-stranded RNA (ssRNA) has largely remained elusive. This is due to challenges in structurally characterizing the unwinding intermediate state and the lack of thermodynamic details underlying this process. In this study, we use solution nuclear magnetic resonance (NMR) spectroscopy to characterize the interaction of human DDX3X, a member of the DDX family, with various RNA ligands. Our results show that the inherent binding affinity of DDX3X for ssRNA is significantly higher than that for structured RNA elements. This preferential binding, accompanied by the formation of a domain-closed conformation in complex with ssRNA, effectively stabilizes the denatured ssRNA state and thus underlies the unwinding activity of DDX3X. Our results provide a thermodynamic and structural basis for the DDX function, whereby DDX can recognize and remodel a distinct set of structured RNAs to participate in a wide range of physiological processes.

Transmembrane conformation of the envelope protein of an alpha coronavirus, NL63

The envelope (E) proteins of coronaviruses (CoVs) form cation-conducting channels that are associated with the pathogenicity of these viruses. To date, high-resolution structural information about these viroporins is limited to the SARS-CoV E protein. To broaden our structural knowledge of other members of this family of viroporins, we now investigate the conformation of the E protein of the human coronavirus (hCoV), NL63. Using two- and three-dimensional magic-angle-spinning NMR, we have measured 13 C and 15 N chemical shifts of the transmembrane domain of E (ETM), which yielded backbone (ϕ, ψ) torsion angles. We further measured the water accessibility of NL63 ETM at neutral pH versus acidic pH in the presence of Ca2+ ions. These data show that NL63 ETM adopts a regular α-helical conformation that is unaffected by pH and the N-terminal ectodomain. Interestingly, the water accessibility of NL63 ETM increases only modestly at acidic pH in the presence of Ca2+ compared to neutral pH, in contrast to SARS ETM, which becomes much more hydrated at acidic pH. This difference suggests a structural basis for the weaker channel conductance of α-CoV compared to β-CoV E proteins. The weaker E channel activity may in turn contribute to the reduced virulence of hCoV-NL63 compared to SARS-CoV viruses.

Exploring Host-Guest Interactions within a 600 kDa DegP Protease Cage Complex Using Hydrodynamics Measurements and Methyl-TROSY NMR

The DegP protease-chaperone operates within the periplasm of Gram-negative bacteria, where it assists in the regulation of protein homeostasis, promotes virulence, and is essential to survival under stress. To carry out these tasks, DegP forms a network of preorganized apo oligomers that facilitate the capture of substrates within distributions of cage-like complexes which expand to encapsulate clients of various sizes. Although the architectures of DegP cage complexes are well understood, little is known about the structures, dynamics, and interactions of client proteins within DegP cages and the relationship between client structural dynamics and function. Here, we probe host-guest interactions within a 600 kDa DegP cage complex throughout the DegP activation cycle using a model α-helical client protein through a combination of hydrodynamics measurements, methyl-transverse relaxation optimized spectroscopy-based solution nuclear magnetic resonance studies, and proteolytic activity assays. We find that in the presence of the client, DegP cages assemble cooperatively with few intermediates. Our data further show that the N-terminal half of the bound client, which projects into the interior of the cages, is predominantly unfolded and flexible, and exchanges between multiple conformational states over a wide range of time scales. Finally, we show that a concerted structural transition of the protease domains of DegP occurs upon client engagement, leading to activation. Together, our findings support a model of DegP as a highly cooperative and dynamic molecular machine that stabilizes unfolded states of clients, primarily via interactions with their C-termini, giving rise to efficient cleavage.

Conformational Preferences at the Glycosidic Linkage of Saccharides in Solution as Deduced from NMR Experiments and MD Simulations: Comparison to Crystal Structures

Glycans are central to information content and regulation in biological systems. These carbohydrate molecules are active either as oligo- or polysaccharides, often in the form of glycoconjugates. The monosaccharide entities are joined by glycosidic linkages and stereochemical arrangements are of utmost importance in determining conformation and flexibility of saccharides. The conformational preferences and population distributions at the glycosidic torsion angles φ and ψ have been investigated for O-methyl glycosides of three disaccharides where the substitution takes place at a secondary alcohol, viz., in α-l-Fucp-(1→3)-β-d-Glcp-OMe, α-l-Fucp-(1→3)-α-d-Galp-OMe and α-d-Glcp-(1→4)-α-d-Galp-OMe, corresponding to disaccharide structural elements present in bacterial polysaccharides. Stereochemical differences at or adjacent to the glycosidic linkage were explored by solution state NMR spectroscopy using one-dimensional 1 H,1 H-NOESY NMR experiments to obtain transglycosidic proton-proton distances and one- and two-dimensional heteronuclear NMR experiments to obtain 3 JCH transglycosidic coupling constants related to torsion angles φ and ψ. Computed effective proton-proton distances from molecular dynamics (MD) simulations showed excellent agreement to experimentally derived distances for the α-(1→3)-linked disaccharides and revealed that for the bimodal distribution at the ψ torsion angle for the α-(1→4)-linked disaccharide experiment and simulation were at variance with each other, calling for further force field developments. The MD simulations disclosed a highly intricate inter-residue hydrogen bonding pattern for the α-(1→4)-linked disaccharide, including a nonconventional hydrogen bond between H5' in the glucosyl residue and O3 in the galactosyl residue, supported by a large downfield 1 H NMR chemical shift displacement compared to α-d-Glcp-OMe. Comparison of population distributions of the glycosidic torsion angles φ and ψ in the disaccharide entities to those of corresponding crystal structures highlighted the potential importance of solvation on the preferred conformation.

Solid-state NMR MAS CryoProbe enables structural studies of human blood protein vitronectin bound to hydroxyapatite

The low sensitivity of nuclear magnetic resonance (NMR) is a major bottleneck for studying biomolecular structures of complex biomolecular assemblies. Cryogenically cooled probe technology overcomes the sensitivity limitations enabling NMR applications to challenging biomolecular systems. Here we describe solid-state NMR studies of the human blood protein vitronectin (Vn) bound to hydroxyapatite (HAP), the mineralized form of calcium phosphate, using a CryoProbe designed for magic angle spinning (MAS) experiments. Vn is a major blood protein that regulates many different physiological and pathological processes. The high sensitivity of the CryoProbe enabled us to acquire three-dimensional solid-state NMR spectra for sequential assignment and characterization of site-specific water-protein interactions that provide initial insights into the organization of the Vn-HAP complex. Vn associates with HAP in various pathological settings, including macular degeneration eyes and Alzheimer's disease brains. The ability to probe these assemblies at atomic detail paves the way for understanding their formation.

Electron and Spin Delocalization in [Co6 Se8 (PEt3 )6 ]0/+1 Superatoms

Molecular clusters can function as nanoscale atoms/superatoms, assembling into superatomic solids, a new class of solid-state materials with designable properties through modifications on superatoms. To explore possibilities on diversifying building blocks, here we thoroughly studied one representative superatom, Co6 Se8 (PEt3 )6 . We probed its structural, electronic, and magnetic properties and revealed its detailed electronic structure as valence electrons delocalize over inorganic [Co6 Se8 ] core while ligands function as an insulated shell. 59 Co SSNMR measurements on the core and 31 P, 13 C on the ligands show that the neutral Co6 Se8 (PEt3 )6 is diamagnetic and symmetric, with all ligands magnetically equivalent. Quantum computations cross-validate NMR results and reveal degenerate delocalized HOMO orbitals, indicating aromaticity. Ligand substitution keeps the inorganic core nearly intact. After losing one electron, the unpaired electron in [Co6 Se8 (PEt3 )6 ]+1 is delocalized, causing paramagnetism and a delocalized electron spin. Notably, this feature of electron/spin delocalization over a large cluster is attractive for special single-electron devices.

Characterization of Levan Fructan Produced by a Gluconobacter japonicus Strain Isolated from a Sugarcane Processing Facility

During raw sugarcane processing, a significant portion of lost sucrose is attributable to microbial degradation. Sucrose consumption by many bacteria is also linked to the production of exopolysaccharides (EPS) such as dextrans and fructans. These resulting EPS cause operational challenges during raw sugar manufacturing. Here, we report the characterization of EPS from a fructan-forming Gluconobacter japonicus bacterium that we previously isolated from a Louisiana sugarcane factory. The genome sequencing revealed the presence of two encoded levansucrase genes, lsrA and lsrB. One levansucrase, LsrB, was detected in the secreted protein fraction of G. japonicus LASM12 by QTOF LC-MS. The spotting assays indicated that G. japonicus produces EPS using sucrose and raffinose as substrates. The G. japonicus EPS correlated with levan fructan commercial standards by 1H-NMR, and with the characteristic carbohydrate fingerprint region for FTIR spectra, confirming that the G. japonicus EPS is levan fructan. The glycosyl composition and glycosyl linkage analysis revealed a linear β-2,6-fructofuranosyl polysaccharide with occasional (5.7%) β-2,1-fructofuranosyl branches. The gel permeation chromatography of the levan fructan EPS showed two main peaks at 4.5 kDa and 8 kDa and a very minor peak at 500 kDa. G. japonicus was identified as a producer of levan fructan. These findings will be useful for future studies aimed at evaluating the impact of levan fructans on sugar crop processing, which have been historically underestimated in industry.

Hydrogen bonding in duplex DNA probed by DNP enhanced solid-state NMR N-H bond length measurements

Numerous biological processes and mechanisms depend on details of base pairing and hydrogen bonding in DNA. Hydrogen bonds are challenging to quantify by X-ray crystallography and cryo-EM due to difficulty of visualizing hydrogen atom locations but can be probed with site specificity by NMR spectroscopy in solution and the solid state with the latter particularly suited to large, slowly tumbling DNA complexes. Recently, we showed that low-temperature dynamic nuclear polarization (DNP) enhanced solid-state NMR is a valuable tool for distinguishing Hoogsteen base pairs (bps) from canonical Watson-Crick bps in various DNA systems under native-like conditions. Here, using a model 12-mer DNA duplex containing two central adenine-thymine (A-T) bps in either Watson-Crick or Hoogsteen confirmation, we demonstrate DNP solid-state NMR measurements of thymine N3-H3 bond lengths, which are sensitive to details of N-H···N hydrogen bonding and permit hydrogen bonds for the two bp conformers to be systematically compared within the same DNA sequence context. For this DNA duplex, effectively identical TN3-H3 bond lengths of 1.055 ± 0.011 Å and 1.060 ± 0.011 Å were found for Watson-Crick A-T and Hoogsteen A (syn)-T base pairs, respectively, relative to a reference amide bond length of 1.015 ± 0.010 Å determined for N-acetyl-valine under comparable experimental conditions. Considering that prior quantum chemical calculations which account for zero-point motions predict a somewhat longer effective peptide N-H bond length of 1.041 Å, in agreement with solution and solid-state NMR studies of peptides and proteins at ambient temperature, to facilitate direct comparisons with these earlier studies TN3-H3 bond lengths for the DNA samples can be readily scaled appropriately to yield 1.083 Å and 1.087 Å for Watson-Crick A-T and Hoogsteen A (syn)-T bps, respectively, relative to the 1.041 Å reference peptide N-H bond length. Remarkably, in the context of the model DNA duplex, these results indicate that there are no significant differences in N-H···N A-T hydrogen bonds between Watson-Crick and Hoogsteen bp conformers. More generally, high precision measurements of N-H bond lengths by low-temperature DNP solid-state NMR based methods are expected to facilitate detailed comparative analysis of hydrogen bonding for a range of DNA complexes and base pairing environments.

Automated High-Order Shimming for Neuroimaging Studies

B0 inhomogeneity presents a significant challenge in MRI and MR spectroscopy, particularly at high-field strengths, leading to image distortion, signal loss, and spectral broadening. Existing high-order shimming methods can alleviate these issues but often require time-consuming and subjective manual selection of regions of interest (ROIs). To address this, we proposed an automated high-order shimming (autoHOS) method, incorporating deep-learning-based brain extraction and image-based high-order shimming. This approach performs automated real-time brain extraction to define the ROI of the field map to be used in the shimming algorithm. The shimming performance of autoHOS was assessed through in vivo echo-planar imaging (EPI) and spectroscopic studies at both 3T and 7T field strengths. AutoHOS outperforms linear shimming and manual high-order shimming, enhancing both the image and spectral quality by reducing the EPI image distortion and narrowing the MRS spectral lineshapes. Therefore, autoHOS demonstrated a significant improvement in correcting B0 inhomogeneity while eliminating the need for additional user interaction.

MagMet: A fully automated web server for targeted nuclear magnetic resonance metabolomics of plasma and serum

Nuclear magnetic resonance (NMR) spectral analysis of biofluids can be a time-consuming process, requiring the expertise of a trained operator. With NMR becoming increasingly popular in the field of metabolomics, there is a growing need to change this paradigm and to automate the process. Here we introduce MagMet, an online web server, that automates the processing and quantification of 1D 1 H NMR spectra from biofluids-specifically, human serum/plasma metabolites, including those associated with inborn errors of metabolism (IEM). MagMet uses a highly efficient data processing procedure that performs automatic Fourier Transformation, phase correction, baseline optimization, chemical shift referencing, water signal removal, and peak picking/peak alignment. MagMet then uses the peak positions, linewidth information, and J-couplings from its own specially prepared standard metabolite reference spectral NMR library of 85 serum/plasma compounds to identify and quantify compounds from experimentally acquired NMR spectra of serum/plasma. MagMet employs linewidth adjustment for more consistent quantification of metabolites from higher field instruments and incorporates a highly efficient data processing procedure for more rapid and accurate detection and quantification of metabolites. This optimized algorithm allows the MagMet webserver to quickly detect and quantify 58 serum/plasma metabolites in 2.6 min per spectrum (when processing a dataset of 50-100 spectra). MagMet's performance was also assessed using spectra collected from defined mixtures (simulating other biofluids), with >100 previously measured plasma spectra, and from spiked serum/plasma samples simulating known IEMs. In all cases, MagMet performed with precision and accuracy matching the performance of human spectral profiling experts. MagMet is available at http://magmet.ca.

5D solid-state NMR spectroscopy for facilitated resonance assignment

1H-detected solid-state NMR spectroscopy has been becoming increasingly popular for the characterization of protein structure, dynamics, and function. Recently, we showed that higher-dimensionality solid-state NMR spectroscopy can aid resonance assignments in large micro-crystalline protein targets to combat ambiguity (Klein et al., Proc. Natl. Acad. Sci. U.S.A. 2022). However, assignments represent both, a time-limiting factor and one of the major practical disadvantages within solid-state NMR studies compared to other structural-biology techniques from a very general perspective. Here, we show that 5D solid-state NMR spectroscopy is not only justified for high-molecular-weight targets but will also be a realistic and practicable method to streamline resonance assignment in small to medium-sized protein targets, which such methodology might not have been expected to be of advantage for. Using a combination of non-uniform sampling and the signal separating algorithm for spectral reconstruction on a deuterated and proton back-exchanged micro-crystalline protein at fast magic-angle spinning, direct amide-to-amide correlations in five dimensions are obtained with competitive sensitivity compatible with common hardware and measurement time commitments. The self-sufficient backbone walks enable efficient assignment with very high confidence and can be combined with higher-dimensionality sidechain-to-backbone correlations from protonated preparations into minimal sets of experiments to be acquired for simultaneous backbone and sidechain assignment. The strategies present themselves as potent alternatives for efficient assignment compared to the traditional assignment approaches in 3D, avoiding user misassignments derived from ambiguity or loss of overview and facilitating automation. This will ease future access to NMR-based characterization for the typical solid-state NMR targets at fast MAS.

Automated identification and quantification of metabolites in human fecal extracts by nuclear magnetic resonance spectroscopy

We report the development of a software program, called MagMet-F, that automates the processing and quantification of 1D 1 H NMR of human fecal extracts. To optimize the program, we identified 82 potential fecal metabolites using 1D 1 H NMR of six human fecal extracts using manual profiling and a literature review of known fecal metabolites. We acquired pure versions of those metabolites and then acquired their 1D 1 H NMR spectra at 700 MHz to generate a fecal metabolite spectral library for MagMet-F. The fitting of these metabolites by MagMet-F was iteratively optimized to replicate manual profiling. We validated MagMet-F's automated profiling using a test set of six fecal extracts. It correctly identified 80% of the compounds and quantified those within <20% of the values determined by manual profiling using Chenomx. We also compared MagMet-F's profiling performance to two other open-access NMR profiling tools, Bayesil and Batman. MagMet-F outperformed both. Bayesil repeatedly overestimated metabolite concentrations by 10% to 40% while Batman was unable to properly quantify any compounds and took 10-20× longer. We have implemented MagMet-F as a freely accessible web server to enable automated, fast and convenient 1D 1 H NMR spectral profiling of fecal samples. MagMet-F is available at https://www.magmet.ca.

Optimal-Control-Based Cβ Chemical Shift Encoding for Efficient Signal Assignment of Solid Proteins

Fast magic-angle spinning (MAS) solid-state NMR spectroscopy is a powerful tool for gaining structural and dynamic information on solid proteins. To access such information site-specifically, the signal assignment process is unavoidable. In the assignment process, Cα and Cβ chemical shifts are of paramount importance in identifying the type of amino acid residues. Conventionally, however, recording the Cβ chemical shift of solid proteins with relatively short transverse relaxation time is often limited by the long delay required for the magnetization transfer to Cβ spins and its evolution, that is, by the sensitivity drop. In this article, we propose a new method that encodes the Cβ chemical shifts onto the intensities of the scalar-coupled Cα signals by combining an optimal control-based spin manipulation pulse and a spin-state filter. This reduces the total required transverse evolution to less than half of that for the previously proposed method, opening up the concept of the Cβ-encoding nearest-neighbor NMR, for the first time, to solid proteins. Also, the total measurement time was shorter than that required for the explicit Cβ shift evolution. We demonstrate the sequential signal assignment for microcrystalline protein GB1, and then discuss the prospects for more challenging proteins.

Deep-Learning-Based Mixture Identification for Nuclear Magnetic Resonance Spectroscopy Applied to Plant Flavors

Nuclear magnetic resonance (NMR) is a crucial technique for analyzing mixtures consisting of small molecules, providing non-destructive, fast, reproducible, and unbiased benefits. However, it is challenging to perform mixture identification because of the offset of chemical shifts and peak overlaps that often exist in mixtures such as plant flavors. Here, we propose a deep-learning-based mixture identification method (DeepMID) that can be used to identify plant flavors (mixtures) in a formulated flavor (mixture consisting of several plant flavors) without the need to know the specific components in the plant flavors. A pseudo-Siamese convolutional neural network (pSCNN) and a spatial pyramid pooling (SPP) layer were used to solve the problems due to their high accuracy and robustness. The DeepMID model is trained, validated, and tested on an augmented data set containing 50,000 pairs of formulated and plant flavors. We demonstrate that DeepMID can achieve excellent prediction results in the augmented test set: ACC = 99.58%, TPR = 99.48%, FPR = 0.32%; and two experimentally obtained data sets: one shows ACC = 97.60%, TPR = 92.81%, FPR = 0.78% and the other shows ACC = 92.31%, TPR = 80.00%, FPR = 0.00%. In conclusion, DeepMID is a reliable method for identifying plant flavors in formulated flavors based on NMR spectroscopy, which can assist researchers in accelerating the design of flavor formulations.

NMR chemical shift prediction and structural elucidation of linker-containing oligo- and polysaccharides using the computer program CASPER

Carbohydrate structures containing alkyl groups as aglycones are useful for investigating enzyme activity and glycan-protein interactions. Moreover, linker-containing oligosaccharides with a spacer group are commonly used to print glycan microarrays or to prepare protein-conjugates as vaccine candidates. The structural accuracy of these synthesized glycans are essential for interpretation of results from biological experiments in which the compounds have been used and NMR spectroscopy can unravel and confirm their structures. An approach for efficient 1H and 13C NMR chemical shift assignments employed a parallel NOAH-10 measurement followed by NMR spin-simulation to refine the 1H NMR chemical shifts, as exemplified for a disaccharide with an azidoethyl group as an aglycone, the NMR chemical shifts of which have been used to enhance the quality of CASPER (http://www.casper.organ.su.se/casper/). The CASPER program has been further developed to aid characterization of linker-containing oligo- and polysaccharides, either by chemical shift prediction for comparison to experimental NMR data or as structural investigation of synthesized glycans based on acquired unassigned NMR data. The ability of CASPER to elucidate structures of linker-containing oligosaccharides is demonstrated and comparisons to assigned or unassigned NMR data show the utility of CASPER in supporting a proposed oligosaccharide structure. Prediction of NMR chemical shifts of an oligosaccharide, corresponding to the repeating unit of an O-antigen polysaccharide, having a linker as an aglycone and a non-natural substituent derivative thereof are presented to exemplify the diversity of structures handled. Furthermore, NMR chemical shift predictions of synthesized polysaccharides, corresponding to bacterial polysaccharides, containing a linker are described showing that in addition to oligosaccharide structures also polysaccharide structures having an aglycone spacer group can be analyzed by CASPER.

Sustainable and cost-effective MAS DNP-NMR at 30 K with cryogenic sample exchange

We report here instrumental developments to achieve sustainable, cost-effective cryogenic Helium sample spinning in order to conduct dynamic nuclear polarisation (DNP) and solid-state NMR (ssNMR) at ultra-low temperatures (<30 K). More specifically, we describe an efficient closed-loop helium system composed of a powerful heat exchanger (95% efficient), a single cryocooler, and a single helium compressor to power the sample spinning and cooling. The system is integrated with a newly designed triple-channel NMR probe that minimizes thermal losses without compromising the radio frequency (RF) performance and spinning stability (±0.05%). The probe is equipped with an innovative cryogenic sample exchange system that allows swapping samples in minutes without introducing impurities in the closeloop system. We report that significant gain in sensitivity can be obtained at 30-40 K on large micro-crystalline molecules with unfavorable relaxation timescales, making them difficult or impossible to polarize at 100 K. We also report rotor-synchronized 2D experiments to demonstrate the stability of the system.

Atomic structure of the open SARS-CoV-2 E viroporin

The envelope (E) protein of the SARS-CoV-2 virus forms cation-conducting channels in the endoplasmic reticulum Golgi intermediate compartment (ERGIC) of infected cells. The calcium channel activity of E is associated with the inflammatory responses of COVID-19. Using solid-state NMR (ssNMR) spectroscopy, we have determined the open-state structure of E's transmembrane domain (ETM) in lipid bilayers. Compared to the closed state, open ETM has an expansive water-filled amino-terminal chamber capped by key glutamate and threonine residues, a loose phenylalanine aromatic belt in the middle, and a constricted polar carboxyl-terminal pore filled with an arginine and a threonine residue. This structure gives insights into how protons and calcium ions are selected by ETM and how they permeate across the hydrophobic gate of this viroporin.

RNA structure and multiple weak interactions balance the interplay between RNA binding and phase separation of SARS-CoV-2 nucleocapsid

The nucleocapsid (N) protein of SARS-CoV-2 binds viral RNA, condensing it inside the virion, and phase separating with RNA to form liquid-liquid condensates. There is little consensus on what differentiates sequence-independent N-RNA interactions in the virion or in liquid droplets from those with specific genomic RNA (gRNA) motifs necessary for viral function inside infected cells. To identify the RNA structures and the N domains responsible for specific interactions and phase separation, we use the first 1,000 nt of viral RNA and short RNA segments designed as models for single-stranded and paired RNA. Binding affinities estimated from fluorescence anisotropy of these RNAs to the two-folded domains of N (the NTD and CTD) and comparison to full-length N demonstrate that the NTD binds preferentially to single-stranded RNA, and while it is the primary RNA-binding site, it is not essential to phase separation. Nuclear magnetic resonance spectroscopy identifies two RNA-binding sites on the NTD: a previously characterized site and an additional although weaker RNA-binding face that becomes prominent when binding to the primary site is weak, such as with dsRNA or a binding-impaired mutant. Phase separation assays of nucleocapsid domains with double-stranded and single-stranded RNA structures support a model where multiple weak interactions, such as with the CTD or the NTD's secondary face promote phase separation, while strong, specific interactions do not. These studies indicate that both strong and multivalent weak N-RNA interactions underlie the multifunctional abilities of N.

On the effects of quadrupolar relaxation in Earth's field NMR spectra

There is growing interest in using low-field magnetic resonance experiments for routine chemical characterization. Earth's field NMR is one such technique that can garner structural information and enable sample differentiation with low cost and highly portable designs. The resulting NMR spectra are primarily influenced by J-couplings, resulting in so-called J-coupled spectra (JCS). Many small molecules include atoms with NMR-active nuclei that are quadrupolar either at natural abundance or are often isotopically enriched (e.g.,2H, 6Li, 11B, 14N, 17O, etc.) where the effects of quadrupolar J-couplings and relaxation on JCS of strongly- and weakly-coupled spin systems have not been explored to date. Herein, using a set of seven fluoropyridine samples with unique substitution and J-couplings, we demonstrate that the 14N relaxation rates can induce drastic line-broadening in the JCS. This includes a previously unexplored unique line broadening mechanism enabled by strongly coupled spins at low-field. Numerical simulations are used to model and refine the magnitudes and signs of J-couplings, as well as indirectly determine the 14N relaxation rates in a single 1D experiment that has a higher fidelity than observed in high-field NMR experiments.

Discovery of a selective and biologically active low-molecular weight antagonist of human interleukin-1β

Human interleukin-1β (hIL-1β) is a pro-inflammatory cytokine involved in many diseases. While hIL-1β directed antibodies have shown clinical benefit, an orally available low-molecular weight antagonist is still elusive, limiting the applications of hIL-1β-directed therapies. Here we describe the discovery of a low-molecular weight hIL-1β antagonist that blocks the interaction with the IL-1R1 receptor. Starting from a low affinity fragment-based screening hit 1, structure-based optimization resulted in a compound (S)-2 that binds and antagonizes hIL-1β with single-digit micromolar activity in biophysical, biochemical, and cellular assays. X-ray analysis reveals an allosteric mode of action that involves a hitherto unknown binding site in hIL-1β encompassing two loops involved in hIL-1R1/hIL-1β interactions. We show that residues of this binding site are part of a conformationally excited state of the mature cytokine. The compound antagonizes hIL-1β function in cells, including primary human fibroblasts, demonstrating the relevance of this discovery for future development of hIL-1β directed therapeutics.

Pharmacological perturbation of the phase-separating protein SMNDC1

SMNDC1 is a Tudor domain protein that recognizes di-methylated arginines and controls gene expression as an essential splicing factor. Here, we study the specific contributions of the SMNDC1 Tudor domain to protein-protein interactions, subcellular localization, and molecular function. To perturb the protein function in cells, we develop small molecule inhibitors targeting the dimethylarginine binding pocket of the SMNDC1 Tudor domain. We find that SMNDC1 localizes to phase-separated membraneless organelles that partially overlap with nuclear speckles. This condensation behavior is driven by the unstructured C-terminal region of SMNDC1, depends on RNA interaction and can be recapitulated in vitro. Inhibitors of the protein's Tudor domain drastically alter protein-protein interactions and subcellular localization, causing splicing changes for SMNDC1-dependent genes. These compounds will enable further pharmacological studies on the role of SMNDC1 in the regulation of nuclear condensates, gene regulation and cell identity.

Continuous flow synthesis of pyridinium salts accelerated by multi-objective Bayesian optimization with active learning

We report a human-in-the-loop implementation of the multi-objective experimental design via a Bayesian optimization platform (EDBO+) towards the optimization of butylpyridinium bromide synthesis under continuous flow conditions. The algorithm simultaneously optimized reaction yield and production rate (or space-time yield) and generated a well defined Pareto front. The versatility of EDBO+ was demonstrated by expanding the reaction space mid-campaign by increasing the upper temperature limit. Incorporation of continuous flow techniques enabled improved control over reaction parameters compared to common batch chemistry processes, while providing a route towards future automated syntheses and improved scalability. To that end, we applied the open-source Python module, nmrglue, for semi-automated nuclear magnetic resonance (NMR) spectroscopy analysis, and compared the acquired outputs against those obtained through manual processing methods from spectra collected on both low-field (60 MHz) and high-field (400 MHz) NMR spectrometers. The EDBO+ based model was retrained with these four different datasets and the resulting Pareto front predictions provided insight into the effect of data analysis on model predictions. Finally, quaternization of poly(4-vinylpyridine) with bromobutane illustrated the extension of continuous flow chemistry to synthesize functional materials.

Strategies for acquisition of resonance assignment spectra of highly dynamic membrane proteins: a GPCR case study

In protein nuclear magnetic resonance (NMR), chemical shift assignment provides a wealth of information. However, acquisition of high-quality solid-state NMR spectra depends on protein-specific dynamics. For membrane proteins, bilayer heterogeneity further complicates this observation. Since the efficiency of cross-polarization transfer is strongly entwined with protein dynamics, optimal temperatures for spectral sensitivity and resolution will depend not only on inherent protein dynamics, but temperature-dependent phase properties of the bilayer environment. We acquired 1-, 2-, and 3D homo- and heteronuclear experiments of the chemokine receptor CCR3 in a 7:3 phosphatidylcholine:cholesterol lipid environment. 1D direct polarization, cross polarization (CP), and T2' experiments indicate sample temperatures below - 25 °C facilitate higher CP enhancement and longer-lived transverse relaxation times. T1rho experiments indicate intermediate timescales are minimized below a sample temperature of - 20 °C. 2D DCP NCA experiments indicated optimal CP efficiency and resolution at a sample temperature of - 30 °C, corroborated by linewidth analysis in 3D NCACX at - 30 °C compared to - 5 °C. This optimal temperature is concluded to be directly related the lipid phase transition, measured to be between - 20 and 15 °C based on rINEPT signal of all-trans and trans-gauche lipid acyl conformations. Our results have critical implications in acquisition of SSNMR membrane protein assignment spectra, as we hypothesize that different lipid compositions with different phase transition properties influence protein dynamics and therefore the optimal acquisition temperature.

Airway and systemic biomarkers of health effects after short-term exposure to indoor ultrafine particles from cooking and candles - A randomized controlled double-blind crossover study among mild asthmatic subjects

BACKGROUND

There is insufficient knowledge about the systemic health effects of exposure to fine (PM_2.5) and ultrafine particles emitted from typical indoor sources, including cooking and candlelight burning. We examined whether short-term exposure to emissions from cooking and burning candles cause inflammatory changes in young individuals with mild asthma. Thirty-six non-smoking asthmatics participated in a randomized controlled double-blind crossover study attending three exposure sessions (mean PM_2.5 µg/m³_; polycyclic aromatic hydrocarbons ng/m³): (a) air mixed with emissions from cooking (96.1; 1.1), (b) air mixed with emissions from candles (89.8; 10), and (c) clean filtered air (5.8; 1.0). Emissions were generated in an adjacent chamber and let into a full-scale exposure chamber where participants were exposed for five hours. Several biomarkers were assessed in relation to airway and systemic inflammatory changes; the primary outcomes of interest were surfactant Protein-A (SP-A) and albumin in droplets in exhaled air - novel biomarkers for changes in the surfactant composition of small airways. Secondary outcomes included cytokines in nasal lavage, cytokines, C-reactive protein (CRP), epithelial progenitor cells (EPCs), genotoxicity, gene expression related to DNA-repair, oxidative stress, and inflammation, as well as metabolites in blood. Samples were collected before exposure start, right after exposure and the next morning.

RESULTS

SP-A in droplets in exhaled air showed stable concentrations following candle exposure, while concentrations decreased following cooking and clean air exposure. Albumin in droplets in exhaled air increased following exposure to cooking and candles compared to clean air exposure, although not significant. Oxidatively damaged DNA and concentrations of some lipids and lipoproteins in the blood increased significantly following exposure to cooking. We found no or weak associations between cooking and candle exposure and systemic inflammation biomarkers including cytokines, CRP, and EPCs.

CONCLUSIONS

Cooking and candle emissions induced effects on some of the examined health-related biomarkers, while no effect was observed in others; Oxidatively damaged DNA and concentrations of lipids and lipoproteins were increased in blood after exposure to cooking, while both cooking and candle emissions slightly affected the small airways including the primary outcomes SP-A and albumin. We found only weak associations between the exposures and systemic inflammatory biomarkers. Together, the results show the existence of mild inflammation following cooking and candle exposure.

Direct Observation of Membrane-Associated H-Ras in the Native Cellular Environment by In-Cell 19F-NMR Spectroscopy

Ras acts as a molecular switch to control intracellular signaling on the plasma membrane (PM). Elucidating how Ras associates with PM in the native cellular environment is crucial for understanding its control mechanism. Here, we used in-cell nuclear magnetic resonance (NMR) spectroscopy combined with site-specific ¹⁹F-labeling to explore the membrane-associated states of H-Ras in living cells. The site-specific incorporation of p-trifluoromethoxyphenylalanine (OCF₃Phe) at three different sites of H-Ras, i.e., Tyr32 in switch I, Tyr96 interacting with switch II, and Tyr157 on helix α5, allowed the characterization of their conformational states depending on the nucleotide-bound states and an oncogenic mutational state. Exogenously delivered ¹⁹F-labeled H-Ras protein containing a C-terminal hypervariable region was assimilated via endogenous membrane-trafficking, enabling proper association with the cell membrane compartments. Despite poor sensitivity of the in-cell NMR spectra of membrane-associated H-Ras, the Bayesian spectral deconvolution identified distinct signal components on three ¹⁹F-labeled sites, thus offering the conformational multiplicity of H-Ras on the PM. Our study may be helpful in elucidating the atomic-scale picture of membrane-associated proteins in living cells.

A Toolbox for Glutamine Use in Dissolution Dynamic Nuclear Polarization: from Enzymatic Reaction Monitoring to the Study of Cellular Metabolic Pathways and Imaging

Glutamine is under scrutiny regarding its metabolic deregulation linked to energetic reprogramming in cancer cells. Many analytical techniques have been used to better understand the impact of the metabolism of amino acids on biological processes, however only a few are suited to work with complex samples. Here, we report the use of a general dissolution dynamic nuclear polarization (D-DNP) formulation using an unexpensive radical as a multipurpose tool to study glutamine, with insights from enzymatic modelling to complex metabolic networks and fast imaging. First, hyperpolarized [5-13 C] glutamine is used as molecular probe to study the kinetic action of two enzymes: L-asparaginase that has been used as an anti-metabolic treatment for cancer, and glutaminase. These results are also compared with those acquired with another hyperpolarized amino acid, [1,4-13 C] asparagine. Second, we explored the use of hyperpolarized (HP) substrates to probe metabolic pathways by monitoring metabolic profiles arising from hyperpolarized glutamine in E. coli extracts. Finally, a highly concentrated sample formulation is proposed for the purpose of fast imaging applications. We think that this approach can be extended to formulate other amino acids as well as other metabolites and provide complementary insights into the analysis of metabolic networks.

Sensitivity considerations on denoising series of spectra by singular value decomposition

When acquiring series of spectra ( T 1 , T 2 , CP buildup curves, etc.) on samples with poor SNR, we are usually faced with choosing between taking a few points with a large number of scans to maximize the SNR or more points with a smaller number of scans to maximize the information content. In this Letter, we show how low-rank decomposition can be used to denoise a series of spectra, reducing the trade-off between the number of scans and the number of experiments.

Early events in amyloid-β self-assembly probed by time-resolved solid state NMR and light scattering

Self-assembly of amyloid-β peptides leads to oligomers, protofibrils, and fibrils that are likely instigators of neurodegeneration in Alzheimer's disease. We report results of time-resolved solid state nuclear magnetic resonance (ssNMR) and light scattering experiments on 40-residue amyloid-β (Aβ40) that provide structural information for oligomers that form on time scales from 0.7 ms to 1.0 h after initiation of self-assembly by a rapid pH drop. Low-temperature ssNMR spectra of freeze-trapped intermediates indicate that β-strand conformations within and contacts between the two main hydrophobic segments of Aβ40 develop within 1 ms, while light scattering data imply a primarily monomeric state up to 5 ms. Intermolecular contacts involving residues 18 and 33 develop within 0.5 s, at which time Aβ40 is approximately octameric. These contacts argue against β-sheet organizations resembling those found previously in protofibrils and fibrils. Only minor changes in the Aβ40 conformational distribution are detected as larger assemblies develop.

Nuclear spin diffusion under fast magic-angle spinning in solid-state NMR

Solid-state nuclear spin diffusion is the coherent and reversible process through which spin order is transferred via dipolar couplings. With the recent increases in magic-angle spinning (MAS) frequencies and magnetic fields becoming routinely applied in solid-state nuclear magnetic resonance, understanding how the increased 1H resolution obtained affects spin diffusion is necessary for interpretation of several common experiments. To investigate the coherent contributions to spin diffusion with fast MAS, we have developed a low-order correlation in Liouville space model based on the work of Dumez et al. (J. Chem. Phys. 33, 224501, 2010). Specifically, we introduce a new method for basis set selection, which accounts for the resonance-offset dependence at fast MAS. Furthermore, we consider the necessity of including chemical shift, both isotropic and anisotropic, in the modeling of spin diffusion. Using this model, we explore how different experimental factors change the nature of spin diffusion. Then, we show case studies to exemplify the issues that arise in using spin diffusion techniques at fast spinning. We show that the efficiency of polarization transfer via spin diffusion occurring within a deuterated and 100% back-exchanged protein sample at 60 kHz MAS is almost entirely dependent on resonance offset. We additionally identify temperature-dependent magnetization transfer in beta-aspartyl L-alanine, which could be explained by the influence of an incoherent relaxation-based nuclear Overhauser effect.

Liquid-Liquid Phase Separation Modifies the Dynamic Properties of Intrinsically Disordered Proteins

Liquid-liquid phase separation of flexible biomolecules has been identified as a ubiquitous phenomenon underlying the formation of membraneless organelles that harbor a multitude of essential cellular processes. We use nuclear magnetic resonance (NMR) spectroscopy to compare the dynamic properties of an intrinsically disordered protein (measles virus N_TAIL) in the dilute and dense phases at atomic resolution. By measuring ¹⁵N NMR relaxation at different magnetic field strengths, we are able to characterize the dynamics of the protein in dilute and crowded conditions and to compare the amplitude and timescale of the different motional modes to those present in the membraneless organelle. Although the local backbone conformational sampling appears to be largely retained, dynamics occurring on all detectable timescales, including librational, backbone dihedral angle dynamics and segmental, chainlike motions, are considerably slowed down. Their relative amplitudes are also drastically modified, with slower, chain-like motions dominating the dynamic profile. In order to provide additional mechanistic insight, we performed extensive molecular dynamics simulations of the protein under self-crowding conditions at concentrations comparable to those found in the dense liquid phase. Simulation broadly reproduces the impact of formation of the condensed phase on both the free energy landscape and the kinetic interconversion between states. In particular, the experimentally observed reduction in the amplitude of the fastest component of backbone dynamics correlates with higher levels of intermolecular contacts or entanglement observed in simulations, reducing the conformational space available to this mode under strongly self-crowding conditions.

Elucidating dynamic anaerobe metabolism with HRMAS 13C NMR and genome-scale modeling

Anaerobic microbial metabolism drives critical functions within global ecosystems, host-microbiota interactions, and industrial applications, yet remains ill-defined. Here we advance a versatile approach to elaborate cellular metabolism in obligate anaerobes using the pathogen Clostridioides difficile, an amino acid and carbohydrate-fermenting Clostridia. High-resolution magic angle spinning nuclear magnetic resonance (NMR) spectroscopy of C. difficile, grown with fermentable 13C substrates, informed dynamic flux balance analysis (dFBA) of the pathogen's genome-scale metabolism. Analyses identified dynamic recruitment of oxidative and supporting reductive pathways, with integration of high-flux amino acid and glycolytic metabolism at alanine's biosynthesis to support efficient energy generation, nitrogen handling and biomass generation. Model predictions informed an approach leveraging the sensitivity of 13C NMR spectroscopy to simultaneously track cellular carbon and nitrogen flow from [U-13C]glucose and [15N]leucine, confirming the formation of [13C,15N]alanine. Findings identify metabolic strategies used by C. difficile to support its rapid colonization and expansion in gut ecosystems.

Three segment ligation of a 104 kDa multi-domain protein by SrtA and OaAEP1

NMR spectroscopy is an excellent tool for studying protein structure and dynamics which provides a deeper understanding of biological function. As the size of the biomolecule of interest increases, it can become advantageous to dilute the number of observed signals in the NMR spectrum to decrease spectral overlap and increase resolution. One way to limit the number of resonances in the NMR data is by selectively labeling a smaller domain within the larger macromolecule, a process called segmental isotopic labeling. Many examples of segmental isotopic labeling have been described where two segments of a protein are ligated together by chemical or enzymatic means, but there are far fewer descriptions of a three or more segment ligation reaction. Herein, we describe an enzymatic segmental labeling scheme that combines the widely used Sortase A and more recently described OaAEP1 for a two site ligation strategy. In preparation to study proposed long-range allostery in the 104 kDa DNA damage repair protein Rad50, we ligated side-chain methyl group labeled Zn Hook domain between two long segments of otherwise unlabeled P.furiosus Rad50. Enzymatic activity data demonstrated that the scars resulting from the ligation reactions did not affect Rad50 function within the Mre11-Rad50 DNA double strand break repair complex. Finally, methyl-based NMR spectroscopy confirmed the formation of the full-length ligated protein. Our strategy highlights the strengths of OaAEP1 for segmental labeling, namely faster reaction times and a smaller recognition sequence, and provides a straightforward template for using these two enzymes in multisite segmental labeling reactions.

Asynchronising five-fold symmetry sequence for better homonuclear polarisation transfer in magic-angle-spinning solid-state NMR

Recoupling, decoupling, and multidimensional correlation experiments in magic-angle-spinning (MAS) solid-state NMR can be designed by exploiting the symmetry of internal spin interactions. One such scheme, namely, C5₂¹, and its supercycled version SPC5₂¹, notated as a five-fold symmetry sequence, is widely used for double-quantum dipole-dipole recoupling. Such schemes are generally rotor synchronised by design. We demonstrate an asynchronous implementation of the SPC5₂¹ sequence leading to higher double-quantum homonuclear polarisation transfer efficiency compared to the normal synchronous implementation. Rotor-synchronisation is broken in two different ways: lengthening the duration of one of the pulses, denoted as pulse-width variation (PWV), and mismatching the MAS frequency denoted as MAS variation (MASV). The application of this asynchronous sequence is shown on three different samples, namely, U-¹³C-alanine and 1,4-¹³C-labelled ammonium phthalate which include ¹³C_α-¹³C_β, ¹³C_α-¹³C_o, and ¹³C_o-¹³C_o spin systems, and adenosine 5'- triphosphate disodium salt trihydrate (ATP⋅3H₂O). We show that the asynchronous version performs better for spin pairs with small dipole-dipole couplings and large chemical-shift anisotropies, for example, ¹³C_o-¹³C_o. Simulations and experiments are shown to corroborate the results.

Practical considerations for the measurement of near-surface electrostatics based on solvent paramagnetic relaxation enhancements

Electrostatic interactions can play important roles in regulating various biological processes. Quantifying surface electrostatics of biomolecules is, therefore, of significant interest. Recent advances in solution NMR spectroscopy have enabled site-specific measurements of de novo near-surface electrostatic potentials (ϕ_ENS) based on a comparison of solvent paramagnetic relaxation enhancements generated from differently charged paramagnetic co-solutes with similar structures. Although the NMR-derived near-surface electrostatic potentials have been shown to agree with theoretical calculations in the context of folded proteins and nucleic acids, such benchmark comparisons may not always be possible, particularly in cases where high-resolution structural models are lacking, such as in the study of intrinsically disordered proteins. Cross-validation of ϕ_ENS potentials can be achieved by comparing values obtained using three pairs of paramagnetic co-solutes, each with a different net charge. Notably we have found cases where agreement of ϕ_ENS potentials between the three pairs is poor and herein we investigate the source of this discrepancy in some detail. We show that for the systems considered here ϕ_ENS potentials obtained from cationic and anionic co-solutes are accurate and that the use of paramagnetic co-solutes with different structures can be a viable option for validation, although the optimal choice of paramagnetic compounds depends on the system of interest.

Transition Metal Dissolution Mechanisms and Impacts on Electronic Conductivity in Composite LiNi0.5Mn1.5O4 Cathode Films

The high-voltage LiNi0.5Mn1.5O4 (LNMO) spinel cathode material offers high energy density storage capabilities without the use of costly Co that is prevalent in other Li-ion battery chemistries (e.g., LiNixMnyCozO2 (NMC)). Unfortunately, LNMO-containing batteries suffer from poor cycling performance because of the intrinsically coupled processes of electrolyte oxidation and transition metal dissolution that occurs at high voltage. In this work, we use operando electron paramagnetic resonance (EPR) and nuclear magnetic resonance (NMR) spectroscopies to demonstrate that transition metal dissolution in LNMO is tightly coupled to HF formation (and thus, electrolyte oxidation reactions as detected with operando and in situ solution NMR), indicative of an acid-driven disproportionation reaction that occurs during delithiation (i.e., battery charging). Leveraging the temporal resolution (s-min) of magnetic resonance, we find that the LNMO particles accelerate the rate of LiPF6 decomposition and subsequent Mn2+ dissolution, possibly due to the acidic nature of terminal Mn-OH groups. X-ray photoemission electron microscopy (XPEEM) provides surface-sensitive and localized X-ray absorption spectroscopy (XAS) measurements, in addition to X-ray photoelectron spectroscopy (XPS), that indicate disproportionation is enabled by surface reconstruction upon charging, which leads to surface Mn3+ sites on the LNMO particle surface that can disproportionate into Mn2+(dissolved) and Mn4+(s). During discharge of the battery, we observe high quantities of metal fluorides (in particular, MnF2) in the cathode electrolyte interphase (CEI) on LNMO as well as the conductive carbon additives in the composite. Electronic conductivity measurements indicate that the MnF2 decreases film conductivity by threefold compared to LiF, suggesting that this CEI component may impede both the ionic and electronic properties of the cathode. Ultimately, to prevent transition metal dissolution and the associated side reactions in spinel-type cathodes (particularly those that operate at high voltages like LNMO), the use of electrolytes that offer improved anodic stability and prevent acid byproducts will likely be necessary.

The Cytoplasmic Domain of the SARS-CoV-2 Envelope Protein Assembles into a β-Sheet Bundle in Lipid Bilayers

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) envelope (E) protein forms a pentameric ion channel in the lipid membrane of the endoplasmic reticulum Golgi intermediate compartment (ERGIC) of the infected cell. The cytoplasmic domain of E interacts with host proteins to cause virus pathogenicity and may also mediate virus assembly and budding. To understand the structural basis of these functions, here we investigate the conformation and dynamics of an E protein construct (residues 8-65) that encompasses the transmembrane domain and the majority of the cytoplasmic domain using solid-state NMR. ¹³C and ¹⁵N chemical shifts indicate that the cytoplasmic domain adopts a β-sheet-rich conformation that contains three β-strands separated by turns. The five subunits associate into an umbrella-shaped bundle that is attached to the transmembrane helices by a disordered loop. Water-edited NMR spectra indicate that the third β-strand at the C terminus of the protein is well hydrated, indicating that it is at the surface of the β-bundle. The structure of the cytoplasmic domain cannot be uniquely determined from the inter-residue correlations obtained here due to ambiguities in distinguishing intermolecular and intramolecular contacts for a compact pentameric assembly of this small domain. Instead, we present four structural topologies that are consistent with the measured inter-residue contacts. These data indicate that the cytoplasmic domain of the SARS-CoV-2 E protein has a strong propensity to adopt β-sheet conformations when the protein is present at high concentrations in lipid bilayers. The equilibrium between the β-strand conformation and the previously reported α-helical conformation may underlie the multiple functions of E in the host cell and in the virion.

Histone H3 core domain in chromatin with different DNA linker lengths studied by 1H-Detected solid-state NMR spectroscopy

Chromatin, a dynamic protein-DNA complex that regulates eukaryotic genome accessibility and essential functions, is composed of nucleosomes connected by linker DNA with each nucleosome consisting of DNA wrapped around an octamer of histones H2A, H2B, H3 and H4. Magic angle spinning solid-state nuclear magnetic resonance (NMR) spectroscopy can yield unique insights into histone structure and dynamics in condensed nucleosomes and nucleosome arrays representative of chromatin at physiological concentrations. Recently we used J-coupling-based solid-state NMR methods to investigate with residue-specific resolution the conformational dynamics of histone H3 N-terminal tails in 16-mer nucleosome arrays containing 15, 30 or 60 bp DNA linkers. Here, we probe the H3 core domain in the 16-mer arrays as a function of DNA linker length via dipolar coupling-based ¹H-detected solid-state NMR techniques. Specifically, we established nearly complete assignments of backbone chemical shifts for H3 core residues in arrays with 15-60 bp DNA linkers reconstituted with ²H,¹³C,¹⁵N-labeled H3. Overall, these chemical shifts were similar irrespective of the DNA linker length indicating no major changes in H3 core conformation. Notably, however, multiple residues at the H3-nucleosomal DNA interface in arrays with 15 bp DNA linkers exhibited relatively pronounced differences in chemical shifts and line broadening compared to arrays with 30 and 60 bp linkers. These findings are consistent with increased heterogeneity in nucleosome packing and structural strain within arrays containing short DNA linkers that likely leads to side-chains of these interfacial residues experiencing alternate conformations or shifts in their rotamer populations relative to arrays with the longer DNA linkers.

NMR spectrum reconstruction as a pattern recognition problem

A new deep neural network based on the WaveNet architecture (WNN) is presented, which is designed to grasp specific patterns in the NMR spectra. When trained at a fixed non-uniform sampling (NUS) schedule, the WNN benefits from pattern recognition of the corresponding point spread function (PSF) pattern produced by each spectral peak resulting in the highest quality and robust reconstruction of the NUS spectra as demonstrated in simulations and exemplified in this work on 2D ¹H-¹⁵N correlation spectra of three representative globular proteins with different sizes: Ubiquitin (8.6 kDa), Azurin (14 kDa), and Malt1 (44 kDa). The pattern recognition by WNN is also demonstrated for successful virtual homo-decoupling in a 2D methyl ¹H-¹³C - HMQC spectrum of MALT1. We demonstrate using WNN that prior knowledge about the NUS schedule, which so far was not been fully exploited, can be used for designing new powerful NMR processing techniques that surpass the existing algorithmic methods.

Less is more: A simple methyl-TROSY based pulse scheme offers improved sensitivity in applications to high molecular weight complexes

The HMQC pulse sequence and variants thereof have been exploited in studies of high molecular weight protein complexes, taking advantage of the fact that fast and slow relaxing magnetization components are sequestered along two distinct magnetization transfer pathways. Despite the simplicity of the HMQC scheme an even shorter version can be designed, based on elimination of the terminal refocusing period, as a further means of increasing signal. Here we present such an experiment, and show that significant sensitivity gains, in some cases by factors of two or more, are realized in studies of proteins varying in molecular masses from 100 kDa to 1 MDa.