15N Chemical Shifts Research Articles

NMR assignment of the conserved bacterial DNA replication protein DnaA domain IV

Chromosomal replication is a ubiquitous and essential cellular process. In bacteria, the master replication initiator DnaA plays a key role in promoting an open complex at the origin (oriC) and recruiting helicase in a tightly regulated process. The C-terminal domain IV specifically recognises consensus sequences of double-stranded DNA in oriC, termed DnaA-boxes, thereby facilitating the initial engagement of DnaA to oriC. Here, we report the 13Cβ and backbone 1H, 15N, and 13C chemical shift assignments of soluble DnaA domain IV from Bacillus subtilis at pH 7.6 and 298 K.

Imino chemical shift assignments of tRNAAsp, tRNAVal and tRNAPhe from Escherichia coli

Transfer RNAs (tRNAs) are an essential component of the protein synthesis machinery. In order to accomplish their cellular functions, tRNAs go through a highly controlled biogenesis process leading to the production of correctly folded tRNAs. tRNAs in solution adopt the characteristic L-shape form, a stable tertiary conformation imperative for the cellular stability of tRNAs, their thermotolerance, their interaction with protein and RNA complexes and their activity in the translation process. The introduction of post-transcriptional modifications by modification enzymes, the global conformation of tRNAs, and their cellular stability are highly interconnected. We aim to further investigate this existing link by monitoring the maturation of bacterial tRNAs in E. coli extracts using NMR. Here, we report on the 1H, 15N chemical shift assignment of the imino groups and some amino groups of unmodified and modified E. coli tRNAAsp, tRNAVal and tRNAPhe, which are essential for characterizing their maturation process using NMR spectroscopy.

The essential synergy of MD simulation and NMR in understanding amorphous drug forms.

Molecular dynamics (MD) simulations and chemical shifts from machine learning are used to predict 15N, 13C and 1H chemical shifts for the amorphous form of the drug irbesartan. The local environments are observed to be highly dynamic well below the glass transition, and averaging over the dynamics is essential to understanding the observed NMR shifts. Predicted linewidths are about 2 ppm narrower than observed experimentally, which is hypothesised to largely result from susceptibility effects. Previously observed differences in the 13C shifts associated with the two tetrazole tautomers can be rationalised in terms of differing conformational dynamics associated with the presence of an intramolecular interaction in one tautomer. 1H shifts associated with hydrogen bonding can also be rationalised in terms of differing average frequencies of transient hydrogen bonding interactions.

Solid-state NMR assignment of α-synuclein polymorph prepared from helical intermediate

Synucleinopathies are neurodegenerative diseases characterized by the accumulation of α-synuclein protein aggregates in the neurons and glial cells. Both ex vivo and in vitro α-synuclein fibrils tend to show polymorphism. Polymorphism results in structure variations among fibrils originating from a single polypeptide/protein. The polymorphs usually have different biophysical, biochemical and pathogenic properties. The various pathologies of a single disease might be associated with distinct polymorphs. Similarly, in the case of different synucleinopathies, each condition might be associated with a different polymorph. Fibril formation is a nucleation-dependent process involving the formation of transient and heterogeneous intermediates from monomers. Polymorphs are believed to arise from heterogeneous oligomer populations because of distinct selection mechanisms in different conditions. To test this hypothesis, we isolated and incubated different intermediates during in vitro fibrillization of α-synuclein to form different polymorphs. Here, we report 13C and 15N chemical shifts and the secondary structure of fibrils prepared from the helical intermediate using solid-state nuclear magnetic spectroscopy.

Including the Ensemble of Unstructured Conformations in the Analysis of Protein's Native State by High-Pressure NMR Spectroscopy.

The analysis of pressure induced changes in the chemical shift of proteins allows statements on structural fluctuations proteins exhibit at ambient pressure. The inherent issue of separating general pressure effects from structural related effects on the pressure dependence of chemical shifts has so far been addressed by considering the characteristics of random coil peptides on increasing pressure. In this work, chemically and pressure denatured states of the cold shock protein B from Bacillus subtilis (BsCspB) have been assigned in 2D 1H-15N HSQC NMR spectra and their dependence on increasing hydrostatic pressure has been evaluated. The pressure denatured polypeptide chain has been used to separate general from structural related effects on 1H and 15N chemical shifts of native BsCspB and the implications on the interpretation of pressure induced changes in the chemical shift regarding the structure of BsCspB are discussed. It has been found that the ensemble of unstructured conformations of BsCspB shows different responses to increasing pressure than random coil peptides do. Thus, the approach used for considering the general effects that arise when hydrostatic pressure increases changes the structural conclusions that are drawn from high pressure NMR spectroscopic experiments that rely on the analysis of chemical shifts.

Analyse des Nativen Zustands von Proteinen durch Hochdruck NMR‐Spektroskopie unter Berücksichtigung des Ensembles Unstrukturierter Konformationen

AbstractDie Analyse druckabhängiger Änderungen der chemischen Verschiebung von Proteinen erlaubt es, Aussagen über strukturelle Fluktuationen zu machen, die diese bereits bei Umgebungsdruck aufweisen. Das Problem, allgemeine Effekte, die mit einer Erhöhung des hydrostatischen Drucks einhergehen, von strukturell bedingten Effekten auf die chemische Verschiebung zu trennen, wurde bisher durch das Einbeziehen der Eigenschaften unstrukturierter Modellpeptide angegangen. In dieser Arbeit wurden chemische Verschiebungen in 2D 1H‐15N‐HSQC Spektren von chemisch als auch druckdenaturierten Zuständen des Kälteschockproteins B aus Bacillus Subtilis (BsCspB) zugeordnet und ihre Abhängigkeit von Änderungen des hydrostatischen Drucks ausgewertet. Die durch hydrostatischen Druck denaturierte Polypeptidkette wurde verwendet, um allgemeine Effekte von strukturbedingten Effekten auf 1H und 15N chemische Verschiebungen zu trennen. Daraus resultierende Interpretationen der druckabhängigen Änderungen der chemischen Verschiebungen werden im strukturellen Kontext diskutiert. Es wurde festgestellt, dass das Ensemble unstrukturierter Konformationen von BsCspB anders auf eine Erhöhung des Drucks reagiert, als es unstrukturierte Modellpeptide tun. Daraus kann geschlussfolgert werden, dass der Ansatz, mit dem man allgemeine Effekte, die durch eine Erhöhung des Drucks hervorgerufen werden, abschätzt, die strukturellen Rückschlüsse ändert, die anhand einer Analyse von chemischen Verschiebungen getroffen werden.

Solution NMR chemical shift assignment of apo and molybdate-bound ModA at two pHs.

ModA is a soluble periplasmic molybdate-binding protein found in most gram-negative bacteria. It is part of the ABC transporter complex ModABC that moves molybdenum into the cytoplasm, to be used by enzymes that carry out various redox reactions. Since there is no clear analog for ModA in humans, this protein could be a good target for antibacterial drug design. Backbone 1H, 13C and 15N chemical shifts of apo and molybdate-bound ModA from E. coli were assigned at pHs 6.0 and 4.5. In addition, side chain atoms were assigned for apo ModA at pH 6.0. When comparing apo and molybdate-bound ModA at pH 6.0, large chemical shift perturbations are observed, not only in areas near the bound metal, but also in regions that are distant from the metal-binding site. Given the significant conformational change between apo and holo ModA, we might expect the large chemical shift changes to be more widespread; however, since they are limited to specific regions, the residues with large perturbations may reveal allosteric sites that could ultimately be important for the design of antibiotics that target ModA.

Raman, ATRIR and NMR spectra, conformational stability, vibrational assignments, normal coordinates analysis and DFT calculations of hydrochlorothiazide drug

The Raman (4000-50 cm−1), attenuated total reflectance (ATR) infrared (4000-400 cm−1) spectra of Hydrochlorothiazide (C7H8ClN3O4S2; HCTZ) have been recorded in the solid state in addition to LCMS in the dimer and trimer ranges. The 1H, 13C and 15N NMR spectra were also recorded in addition to DEPT (45, 90 and 135) of HCTZ sample dissolved DMSO‑d6. Initially, three conformers resulting from the rotation of -SO2NH2 moiety around C9-S11 were proposed, namely gauche, trans and cis conformers. And then Density Functional Theory (DFT) calculations for an isolated molecule were carried out employing both B3LYP and P3PW91 methods at 6-31G(d), 6-311+G(d), 6–311++G(d,p) basis set. The cis isomer is a transition state owing to the prediction of an imaginary wavenumber, while the gauche being more stable than trans by ∼438–538 cm−1 (∼1.25–1.54 kcal) employing the afore-mentioned computations. Moreover, we carried out a potential surface scan for the gauche to determine SO2NH2 barriers to internal rotations and to ensure the optimal gauche is the most stable conformer. Aided by normal coordinate analysis (NCA) in addition to the atomic displacements in Cartesian coordinates (ADCC), confident vibrational assignment of the observed infrared and Raman bands are proposed for all fundamentals of the gauche rotamer. Nevertheless, force constants (FCs) in internal coordinates were also estimated using the above-mentioned methods and basis sets as well. The 1H, 13C and 15N chemical shifts were also correlated to their corresponding atoms. The results are reported herein and compared with HCTZ and its cocrystals whenever appropriate.

Backbone 1H, 13C, and 15N chemical shift assignments for human SERF2.

Human small EDRK-rich factor protein SERF2 is a cellular driver of protein amyloid formation, a process that has been linked to neurodegenerative diseases including Alzheimer's and Parkinson's disease. SERF2 is a 59 amino acid protein, highly charged, and well conserved whose structure and physiological function is unclear. SERF family proteins including human SERF2 have shown a tendency to form fuzzy complexes with misfolded proteins such as α-Synuclein which has been linked to Parkinson's disease. SERF family proteins have been recently identified to bind nucleic acids, but the binding mechanism(s) remain enigmatic. Here, using multidimensional solution NMR, we report the 1H, 15N, and 13C chemical shift assignments (~ 86% of backbone resonance assignments) for human SERF2. TALOS-N predicted secondary structure of SERF2 showed three very short helices (3-4 residues long) in the N-terminal region of the protein and a long helix in the C-terminal region spanning residues 37-46 which is consistent with the helical content indicated by circular dichroism spectroscopy. Paramagnetic relaxation enhancement NMR analysis revealed that a short C-terminal region E53-K55 is in the proximity of the N-terminus. Having the backbone assignment of SERF2 allowed us to probe its interaction with α-Synuclein and to identify the residues in SERF2 binding interfaces that likely promote α-Synuclein aggregation.

Predicted and Experimental NMR Chemical Shifts at Variable Temperatures: The Effect of Protein Conformational Dynamics.

NMR chemical shifts provide a sensitive probe of protein structure and dynamics but remain challenging to predict and interpret. We examine the effect of protein conformational distributions on 15N chemical shifts for dihydrofolate reductase (DHFR), comparing QM/MM predicted shifts with experimental shifts in solution as well as frozen distributions. Representative snapshots from MD trajectories exhibit variation in predicted 15N chemical shifts of up to 25 ppm. The average over the fluctuations is in significantly better agreement with room temperature solution experimental values than the prediction for any single optimal conformations. Meanwhile, solid-state NMR (SSNMR) measurements of frozen solutions at 105 K exhibit broad lines whose widths agree well with the widths of distributions of predicted shifts for samples from the trajectory. The backbone torsion angle ψi-1 varies over 60° on the picosecond time scale, compensated by φi. These fluctuations can explain much of the shift variation.

The Scope of the Applicability of Non-relativistic DFT Calculations of NMR Chemical Shifts in Pyridine-Metal Complexes for Applied Applications.

Heavy metals are toxic, but it is impossible to stop using them. Considering the variety of molecular systems in which they can be present, the multicomponent nature and disorder of the structure of such systems, one of the most effective methods for studying them is NMR spectroscopy. This determines the need to calculate NMR chemical shifts for expected model systems. For elements beyond the third row of the periodic table, corrections for relativistic effects are necessary when calculating NMR parameters. Such corrections may be necessary even for light atoms due to the shielding effect of a neighboring heavy atom. This work examines the extent to which non-relativistic DFT calculations are able to reproduce experimental 15N and 113Cd NMR chemical shift tensors in pyridine-metal coordination complexes. It is shown that while for the calculation of 15N NMR chemical shift tensors there is no real need to consider relativistic corrections, for 113Cd, on the contrary, none of the tested calculation methods could reproduce the experimentally obtained tensor to any extent correctly.

The interplay of density functional selection and crystal structure for accurate NMR chemical shift predictions.

Ab initio chemical shift prediction plays a central role in nuclear magnetic resonance (NMR) crystallography, and the accuracy with which chemical shifts can be predicted relative to experiment impacts the confidence with which structures can be assigned. For organic crystals, periodic density functional theory calculations with the gauge-including projector augmented wave (GIPAW) approximation and the PBE functional are widely used at present. Many previous studies have examined how using more advanced density functionals can increase the accuracy of predicted chemical shifts relative to experiment, but nearly all of those studies employed crystal structures that were optimized with generalized-gradient approximation (GGA) functionals. Here, we investigate how the accuracy of the predicted chemical shifts in organic crystals is affected by replacing GGA-level PBE-D3(BJ) crystal geometries with more accurate hybrid functional PBE0-D3(BJ) ones. Based on benchmark data sets containing 132 13C and 35 15N chemical shifts, plus case studies on testosterone, acetaminophen, and phenobarbital, we find that switching from GGA-level geometries and chemical shifts to hybrid-functional ones reduces 13C and 15N chemical shift errors by ∼40-60% versus experiment. However, most of the improvement stems from the use of the hybrid functional for the chemical shift calculations, rather than from the refined geometries. In addition, even with the improved geometries, we find that double-hybrid functionals still do not systematically increase chemical shift agreement with experiment beyond what hybrid functionals provide. In the end, these results suggest that the combination of GGA-level crystal structures and hybrid-functional chemical shifts represents a particularly cost-effective combination for NMR crystallography in organic systems.

NMR analysis of structural geometry and molecular dynamics in perovskite-type N(CH3)4CdBr3 crystal near high-temperature phase transition.

The NMR chemical shifts, linewidths, spin-lattice relaxation times in the rotating system T1ρ, and spin-lattice relaxation times in the laboratory system T1 were evaluated for the perovskite-type N(CH3)4CdBr3 crystal, aiming to understand the changes in the structural geometry and molecular dynamics from phase I to phase II. From the temperature-dependence of the 1H, 13C, 14N, and 113Cd NMR chemical shifts, the structural geometry underwent a continuous change, without anomalous changes around (TC = 390 K). However, the linewidths in phase I were narrower than those in phase II, indicating that the motional averaging effects were caused by the rapid rotation of the N(CH3)4 group. Sudden changes in T1 and T1ρ were observed near TC, for which the activation energy Ea in phase I was approximately 12 times larger than that in phase II; the small Ea values in phase II indicate a large degree of freedom for the methyl group and CdBr6 octahedra, whereas the large Ea in phase I was primarily attributed to the overall N(CH3)4 and the 113Cd in the CdBr6 groups. Consequently, the phase transition mechanisms of N(CH3)4CdBr3 are related to reorientation of the N(CH3)4 group and the arrangement of the CdBr6 groups.

Distinguishing Different Hydrogen-Bonded Helices in Proteins by Efficient 1H-Detected Three-Dimensional Solid-State NMR.

Helical structures in proteins include not only α-helices but also 310 and π helices. These secondary structures differ in the registry of the C═O···H-N hydrogen bonds, which are i to i + 4 for α-helices, i to i + 3 for 310 helices, and i to i + 5 for π-helices. The standard NMR observable of protein secondary structures are chemical shifts, which are, however, insensitive to the precise type of helices. Here, we introduce a three-dimensional (3D) 1H-detected experiment that measures and assigns CO-HN cross-peaks to distinguish the different types of hydrogen-bonded helices. This hCOhNH experiment combines efficient cross-polarization from CO to HN with 13C, 15N, and 1H chemical shift correlation to detect the relative proximities of the COi-Hi+jN spin pairs. We demonstrate this experiment on the membrane-bound transmembrane domain of the SARS-CoV-2 envelope (E) protein (ETM). We show that the C-terminal five residues of ETM form a 310-helix, whereas the rest of the transmembrane domain have COi-Hi+4N hydrogen bonds that are characteristic of α-helices. This result confirms the recent high-resolution solid-state NMR structure of the open state of ETM, which was solved in the absence of explicit hydrogen-bonding restraints. This C-terminal 310 helix may facilitate proton and calcium conduction across the hydrophobic gate of the channel. This hCOhNH experiment is generally applicable and can be used to distinguish not only different types of helices but also different types of β-strands and other hydrogen-bonded conformations in proteins.

15N NMR study of azo-hydrazo tautomerism in some industrially produced azo dyes

15N NMR data (1J (15N, 1H) coupling constants and 15N chemical shifts) were used to characterize quantitatively azo-hydrazo tautomerism in some industrially produced azo dyes. Both 1J (15N, 1H) coupling constants and 15N chemical shifts were obtained at the natural abundance level of 15N using inverse detection of resonances. All three possible types of azo dye existence were found: (i) dyes existing completely in azo forms, (ii) dyes existing completely in hydrazo forms or (iii) dyes forming azo-hydrazo tautomeric equilibria.

Automated Fragmentation Quantum Mechanical Calculation of 15N and 13C Chemical Shifts in a Membrane Protein.

In this work, we developed an accurate and cost-effective automated fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) method to calculate the chemical shifts of 15N and 13C of membrane proteins. The convergence of the AF-QM/MM method was tested using Krokinobacter eikastus rhodopsin 2 as a test case. When the distance threshold of the QM region is equal to or larger than 4.0 Å, the results of the AF-QM/MM calculations are close to convergence. In addition, the effects of selected density functionals, basis sets, and local chemical environment of target atoms on the chemical shift calculations were systematically investigated. Our results demonstrate that the predicted chemical shifts are more accurate when important environmental factors including cross-protomer interactions, lipid molecules, and solvent molecules are taken into consideration, especially for the 15N chemical shift prediction. Furthermore, with the presence of sodium ions in the environment, the chemical shift of residues, retinal, and retinal Schiff base are affected, which is consistent with the results of the solid-state nuclear magnetic resonance (NMR) experiment. Upon comparing the performance of various density functionals (namely, B3LYP, B3PW91, M06-2X, M06-L, mPW1PW91, OB95, and OPBE), the results show that mPW1PW91 is a suitable functional for the 15N and 13C chemical shift prediction of the membrane proteins. Meanwhile, we find that the improved accuracy of the 13Cβ chemical shift calculations can be achieved by the employment of the triple-ζ basis set. However, the employment of the triple-ζ basis set does not improve the accuracy of the 15N and 13Cα chemical shift calculations nor does the addition of a diffuse function improve the overall prediction accuracy of the chemical shifts. Our study also underscores that the AF-QM/MM method has significant advantages in predicting the chemical shifts of key ligands and nonstandard residues in membrane proteins than most widely used empirical models; therefore, it could be an accurate computational tool for chemical shift calculations on various types of biological systems.

Revisiting the influence of pH on 1JCαH and chemical shifts of glycine and alanine short oligopeptides.

The pH dependence of several NMR parameters of glycine and alanine short oligopeptides has been reported previously in different studies. Here we have thoroughly examined, summarized and demonstrated the dependences of 1H, 13C and 15N chemical shifts and protonation states of amino acids using two-dimensional NMR experiments. Nevertheless, 1JCαH one bond spin-spin coupling constants are more informative and convenient for determination of the position and protonation state of glycine and alanine residue in the oligopeptide chain. In particular, for various oligopeptides (up to six residues), it was shown that the pH dependence of 1JCαH of N-terminal glycine and alanine residues is larger than that of C-terminal groups, and in backbone residues, it is not influenced by pH and only slightly depends on the position of the amino acid residue in the chain.

1H, 13C, and 15N NMR chemical shift assignment of LytM N-terminal domain (residues 26–184)

Antibiotic resistance is a growing problem and a global threat for modern healthcare. New approaches complementing the traditional antibiotic drugs are urgently needed to secure the ability to treat bacterial infections also in the future. Among the promising alternatives are bacteriolytic enzymes, such as the cell wall degrading peptidoglycan hydrolases. Staphylococcus aureus LytM, a Zn2+-dependent glycyl-glycine endopeptidase of the M23 family, is one of the peptidoglycan hydrolases. It has a specificity towards staphylococcal peptidoglycan, making it an interesting target for antimicrobial studies. LytM hydrolyses the cell wall of S. aureus, a common pathogen with multi-resistant strains that are difficult to treat, such as the methicillin-resistant S. aureus, MRSA. Here we report the 1H, 15N and 13C chemical shift assignments of S. aureus LytM N-terminal domain and linker region, residues 26–184. These resonance assignments can provide the basis for further studies such as elucidation of structure and interactions.

NMR assignment and dynamics of the dimeric form of soluble C-terminal domain major ampullate spidroin 2 from Latrodectus hesperus.

Spider dragline silk has attracted great interest due to its outstanding mechanical properties, which exceed those of man-made synthetic materials. Dragline silk, which is composed of at least major ampullate spider silk protein 1 and 2 (MaSp1 and MaSp2), contains a long repetitive domain flanked by N-terminal and C-terminal domains (NTD and CTD). Despite the small size of the CTD, this domain plays a crucial role as a molecular switch that regulates and directs spider silk self-assembly. In this study, we report the 1H, 13C, and 15N chemical shift assignments of the Latrodectus hesperus MaSp2 CTD in dimeric form at pH 7. Our solution NMR data demonstrated that this protein contains five helix regions connected by a flexible linker.

Charge transport in protic ionic liquids: Effect of protonation state in 1-methylimidazolium – acetate/trifluoroacetate mixtures

Despite the potential to overcome slow charge transport in ionic liquids by decoupling charge transport from mass transport in protic ionic liquids (PILs), the average proton distribution and the contribution of proton conduction to the overall charge transport in PILs has not been fully established. To elucidate the predominant underlying molecular species in the PILs and the molecular dynamics of the PILs, we study the dipolar response and conductivity using dielectric relaxation (DR) spectroscopy and the nucleis chemical environment and molecular mobility using nuclear magnetic resonance (NMR) spectroscopy of prototypical PILs based on 1–methylimidazolium as Brønsted base. To elucidate the effect of the acidity of the Brønsted acid we gradually vary the composition of the PILs from 1–methylimidazoliumacetate [MimH][AcO] to 1–methylimidazoliumtrifluoroacetate [MimH][TFA]. Both, 15N-NMR chemical shifts and the sample conductivities suggest that in neat [MimH][AcO] electro-neutral molecular species dominate, consistent with earlier findings. The dipolar response of these molecular species as detected via the DR spectra is stronger than what one would expect based on the molecular dipole moments estimated from density functional theory calculations, which points to the formation of an extended hydrogen-bonded network. Upon addition of [MimH][TFA], this dipolar response is reduced as the degree of protonation of 1–methylimidazolium gradually increases, which is apparent from the conductance and 15N chemical shifts. Interestingly, diffusion ordered NMR spectroscopy indicates protons to become more mobile than the constituting acid and base molecules, consistent with Grotthuss-like proton transport. To vary the immediate environment of the PILs, we study solutions of the PIL mixtures in aprotic dimethylformamide (DMF) and protic methanol (MeOH) as two limiting cases. We find that [MimH][AcO] exclusively forms electroneutral species in DMF, as evident from the negligible conductivity, while with increasing TFA content the fraction of charged molecular species increases. Although these ionic species give rise to a finite conductivity, demonstrating the presence of mobile ions, a large fraction of these ions is bound in [MimH][TFA] ion-pairs, as detected in the DR spectra. Yet, temperature dependent 1H-NMR linewidths of the acidic proton signals suggest proton exchange to occur on the NMR timescale, presumably between [MimH]+ and [Mim]. Conversely, in MeOH we find no evidence for the presence of ion-pairs, while the high conductivity of these solutions shows that MeOH can efficiently solvate ions and may participate in proton transfer. Together, our results highlight the sensitivity of protonation equilibria to the immediate environment and of proton dynamics to these protonation equilibria. Our results show that tuning these equilibria using PIL mixtures can provide routes to tailor the contribution of proton transport to charge transport in PILs.