Chemical shifts in amino acids, peptides, and proteins: from quantum chemistry to drug design

Deep Learning Protocol for Predicting Full-Spectrum Infrared and Raman Spectra of Polypeptides and Proteins Using All-Atom Models

Infrared (IR) spectroscopy and Raman spectroscopy are powerful tools for probing protein and peptide structures due to their capability to provide molecular fingerprints. As a popular spectral simulation method, the quantum chemistry (QC) calculation is usually hampered by the high computational cost and low efficiency. In this study, we developed a comprehensive data set of IR and Raman spectra for amino acids, dipeptides, and tripeptides. Using this data set, we applied transfer learning with DetaNet (a deep equivariant tensor attention network) to simulate full-spectrum IR and Raman spectra for large polypeptides and proteins. We have demonstrated that the transfer-learned DetaNet (TL-DetaNet) model successfully simulated the vibrational spectra of proteins with thousands of atoms, far exceeding traditional QC limitations. Additionally, TL-DetaNet achieved an efficiency that was 103-105 times greater than that of QC methods. This work highlights the importance of data sets in machine learning and positions transfer learning as a transformative tool for large-scale biomolecular simulations, marking a substantial advancement in protein vibrational spectroscopy.

Tailoring solid-state DNP methods to the study of α-synuclein LLPS

Dynamic Nuclear Polarization (DNP) is a technique that leverages the quantum sensing capability of electron spins to enhance the sensitivity of nuclear magnetic resonance (NMR) signals, especially for insensitive samples. Glassing agents play a crucial role in the DNP process by facilitating the transfer of polarization from the unpaired electron spins to the nuclear spins along with cryoprotection of biomolecules. DNPjuice comprising of glycerol-d8/D2O/H2O has been extensively used for this purpose over the past two decades. Polyethylene glycol (PEG), also used as a cryoprotectant, is often used as a crowding agent in experimental setups to mimic cellular conditions, particularly the invitro preparation of liquid-liquid phase separated (LLPS) condensates. In this study, we investigate the efficacy of PEG as an alternative to glycerol in the DNP juice, critical for signal enhancement. The modified DNP matrix leads to high DNP enhancement which enables direct study of LLPS condensates by solid-state DNP methods without adding any external constituents. An indirect advantage of employing PEG is that the PEG signals appear at ∼72.5 ppm and are relatively well-separated from the aliphatic region of the protein spectra. Large cross-effect DNP enhancement is attained for 13C-glycine by employing the PEG-water mixture as a glassing agent and ASYMPOL-POK as the state-of-art polarizing agent, without any deuteration. The DNP enhancement and the buildup rates are similar to results obtained with DNP juice, conforming to that PEG serves as a good candidate for both inducing crowding and glassing agent in the study of LLPS.

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.

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. Prediction of shifts, and therefore interpretation of shifts, particularly for the frequently measured amidic 15 N sites, remains a tall challenge. We demonstrate that protein 15 N chemical shift prediction from QM/MM predictions can be improved if conformational variation is included via MD sampling, focusing on the antibiotic target, E. coli Dihydrofolate reductase (DHFR). Variations of up to 25 ppm in predicted 15 N chemical shifts are observed over the trajectory. For solution shifts the average of fluctuations on the low picosecond timescale results in a superior prediction to a single optimal conformation. For low temperature solid state measurements, the histogram of predicted shifts for locally minimized snapshots with specific solvent arrangements sampled from the trajectory explains the heterogeneous linewidths; in other words, the conformations and associated solvent are 'frozen out' at low temperatures and result in inhomogeneously broadened NMR peaks. We identified conformational degrees of freedom that contribute to chemical shift variation. Backbone torsion angles show high amplitude fluctuations during the trajectory on the low picosecond timescale. For a number of residues, including I60, ψ varies by up to 60º within a conformational basin during the MD simulations, despite the fact that I60 (and other sites studied) are in a secondary structure element and remain well folded during the trajectory. Fluctuations in ψ appear to be compensated by other degrees of freedom in the protein, including φ of the succeeding residue, resulting in "rocking" of the amide plane with changes in hydrogen bonding interactions. Good agreement for both room temperature and low temperature NMR spectra provides strong support for the specific approach to conformational averaging of computed chemical shifts.

Atomic Resolution Insights into pH Shift Induced Deprotonation Events in LS-Shaped Aβ(1-42) Amyloid Fibrils

Alzheimer's disease is a neurodegenerative disorder associated with the deposition of misfolded aggregates of the amyloid-β protein (Aβ). Aβ(1-42) is one of the most aggregation-prone components in senile plaques of AD patients. We demonstrated that relatively homogeneous Aβ(1-42) fibrils with one predominant fold visible in solid-state NMR spectra can be obtained at acidic pH. The structure of these fibrils differs remarkably from some other polymorphs obtained at neutral pH. In particular, the entire N-terminal region is part of the rigid fibril core. Here, we investigate the effects of a pH shift on the stability and the fold of these fibrils at higher pH values. Fibril bundling at neutral pH values renders cryo-EM studies impractical, but solid-state NMR spectroscopy, molecular dynamics simulations, and biophysical methods provide residue-specific structural information under these conditions. The LS-fold of the Aβ(1-42) fibrils does not change over the complete pH range from pH 2 to pH 7; in particular, the N-terminus remains part of the fibril core. We observe changes in the protonation state of charged residues starting from pH 5 on a residue-specific level. The deprotonation of the C-terminal carboxyl group of A42 in the intermolecular salt bridge with D1 and K28 is slow on the NMR time scale, with a local pK_a of 5.4, and local conformations of the involved residues are affected by deprotonation of A42. Thus, we demonstrate that this fibril form is stable at physiological pH values.

Quantum Chemistry Calculations for Metabolomics

A primary goal of metabolomics studies is to fully characterize the small-molecule composition of complex biological and environmental samples. However, despite advances in analytical technologies over the past two decades, the majority of small molecules in complex samples are not readily identifiable due to the immense structural and chemical diversity present within the metabolome. Current gold-standard identification methods rely on reference libraries built using authentic chemical materials ("standards"), which are not available for most molecules. Computational quantum chemistry methods, which can be used to calculate chemical properties that are then measured by analytical platforms, offer an alternative route for building reference libraries, i.e., in silico libraries for "standards-free" identification. In this review, we cover the major roadblocks currently facing metabolomics and discuss applications where quantum chemistry calculations offer a solution. Several successful examples for nuclear magnetic resonance spectroscopy, ion mobility spectrometry, infrared spectroscopy, and mass spectrometry methods are reviewed. Finally, we consider current best practices, sources of error, and provide an outlook for quantum chemistry calculations in metabolomics studies. We expect this review will inspire researchers in the field of small-molecule identification to accelerate adoption of in silico methods for generation of reference libraries and to add quantum chemistry calculations as another tool at their disposal to characterize complex samples.

Automated fragmentation quantum mechanical calculation of 13C and 1H chemical shifts in molecular crystals

In this work, the automated fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) approach was applied to calculate the ¹³C and ¹H nuclear magnetic resonance (NMR) chemical shifts in molecular crystals. Two benchmark sets of molecular crystals were selected to calculate the NMR chemical shifts. Systematic investigation was conducted to examine the convergence of AF-QM/MM calculations and the impact of various density functionals with different basis sets on the NMR chemical shift prediction. The result demonstrates that the calculated NMR chemical shifts are close to convergence when the distance threshold for the QM region is larger than 3.5 Å. For ¹³C chemical shift calculations, the mPW1PW91 functional is the best density functional among the functionals chosen in this study (namely, B3LYP, B3PW91, M06-2X, M06-L, mPW1PW91, OB98, and OPBE), while the OB98 functional is more suitable for the ¹H NMR chemical shift prediction of molecular crystals. Moreover, with the B3LYP functional, at least a triple-ζ basis set should be utilized to accurately reproduce the experimental ¹³C and ¹H chemical shifts. The employment of diffuse basis functions will further improve the accuracy for ¹³C chemical shift calculations, but not for the ¹H chemical shift prediction. We further proposed a fragmentation scheme of dividing the central molecule into smaller fragments. By comparing with the results of the fragmentation scheme using the entire central molecule as the core region, the AF-QM/MM calculations with the fragmented central molecule can not only achieve accurate results but also reduce the computational cost. Therefore, the AF-QM/MM approach is capable of predicting the ¹³C and ¹H NMR chemical shifts for molecular crystals accurately and effectively, and could be utilized for dealing with more complex periodic systems such as macromolecular polymers and biomacromolecules. The AF-QM/MM program for molecular crystals is available at https://github.com/shiman1995/NMR.

Accurate and cost-effective NMR chemical shift predictions for proteins using a molecules-in-molecules fragmentation-based method

We have developed an efficient protocol using our two-layer Molecules-in-Molecules (MIM2) fragmentation-based quantum chemical method for the prediction of NMR chemical shifts of large biomolecules. To investigate the performance of our fragmentation approach and demonstrate its applicability, MIM-NMR calculations are first calibrated on a test set of six proteins. The MIM2-NMR method yields a mean absolute deviation (MAD) from unfragmented full molecule calculations of 0.01 ppm for 1H and 0.06 ppm for 13C chemical shifts. Thus, the errors from fragmentation are only about 3% of our target accuracy of ∼0.3 ppm for 1H and 2-3 ppm for 13C chemical shifts. To compare with experimental chemical shifts, a standard protocol is first derived using two smaller proteins 2LHY (176 atoms) and 2LI1 (146 atoms) for obtaining an appropriate protein structure for NMR chemical shift calculations. The effect of the solvent environment on the calculated NMR chemical shifts is incorporated through implicit, explicit, or explicit-implicit solvation models. The expensive first solvation shell calculations are replaced by a micro-solvation model in which only the immediate interaction between the protein and the explicit solvation environment is considered. A single explicit water molecule for each amine and amide proton is found to be sufficient to yield accurate results for 1H chemical shifts. The 1H and 13C NMR chemical shifts calculated using our protocol give excellent agreement with experiments for two larger proteins, 2MC5 (the helical part with 265 atoms) and 3UMK (33 residue slice with 547 atoms). Overall, our target accuracy of ∼0.3 ppm for 1H and ∼2-3 ppm for 13C has been achieved for the larger proteins. The proposed MIM-NMR method is accurate and computationally cost-effective and should be applicable to study a wide range of large proteins.

Determination of accurate backbone chemical shift tensors in microcrystalline proteins by integrating MAS NMR and QM/MM

Chemical shifts are highly sensitive probes of local conformation and overall structure. Both isotropic shifts and chemical shift tensors are readily accessible from NMR experiments but their quantum mechanical calculations remain challenging. In this work, we report and compare accurately measured and calculated 15NH and 13Cα chemical shift tensors in proteins, using the microcrystalline agglutinin from Oscillatoria agardhii (OAA). Experimental 13Cα and 15NH chemical tensors were obtained by solid-state NMR spectroscopy, employing tailored recoupling sequences, and for their quantum mechanics/molecular mechanics (QM/MM) calculations different sets of functionals were evaluated. We show that 13Cα chemical shift tensors are primarily determined by backbone dihedral angles and dynamics, while 15NH tensors mainly depend on local electrostatic contributions from solvation and hydrogen bonding. In addition, the influence of including crystallographic waters, the molecular mechanics geometry optimization protocol, and the level of theory on the accuracy of the calculated chemical shift tensors is discussed. Specifically, the power of QM/MM calculations in accurately predicting the unusually upfield shifted 1HN G26 and G93 resonances is highlighted. Our integrated approach is expected to benefit structure refinement of proteins and protein assemblies.

Evaluating electronic structure methods for accurate calculation of 19 F chemical shifts in fluorinated amino acids

The ability of electronic structure methods (11 density functionals, HF, and MP2 calculations; two basis sets and two solvation models) to accurately calculate the ¹⁹ F chemical shifts of 31 structures of fluorinated amino acids and analogues with known experimental ¹⁹ F NMR spectra has been evaluated. For this task, BHandHLYP, ωB97X, and Hartree-Fock with scaling factors (provided within) are most accurate. Additionally, the accuracy of methods to calculate relative changes in fluorine shielding across 23 sets of structural variants, such as zwitterionic amino acids versus side chains only, was also determined. This latter criterion may be a better indicator of reliable methods for the ultimate goal of assigning and interpreting chemical shifts of fluorinated amino acids in proteins. It was found that MP2 and M062X calculations most accurately assess changes in shielding among analogues. These results serve as a guide for computational developments to calculate ¹⁹ F chemical shifts in biomolecular environments. © 2017 Wiley Periodicals, Inc.

Unraveling the meaning of chemical shifts in protein NMR

Chemical shifts are among the most informative parameters in protein NMR. They provide wealth of information about protein secondary and tertiary structure, protein flexibility, and protein-ligand binding. In this report, we review the progress in interpreting and utilizing protein chemical shifts that has occurred over the past 25years, with a particular focus on the large body of work arising from our group and other Canadian NMR laboratories. More specifically, this review focuses on describing, assessing, and providing some historical context for various chemical shift-based methods to: (1) determine protein secondary and super-secondary structure; (2) derive protein torsion angles; (3) assess protein flexibility; (4) predict residue accessible surface area; (5) refine 3D protein structures; (6) determine 3D protein structures and (7) characterize intrinsically disordered proteins. This review also briefly covers some of the methods that we previously developed to predict chemical shifts from 3D protein structures and/or protein sequence data. It is hoped that this review will help to increase awareness of the considerable utility of NMR chemical shifts in structural biology and facilitate more widespread adoption of chemical-shift based methods by the NMR spectroscopists, structural biologists, protein biophysicists, and biochemists worldwide. This article is part of a Special Issue entitled: Biophysics in Canada, edited by Lewis Kay, John Baenziger, Albert Berghuis and Peter Tieleman.

Computing the Diamagnetic Susceptibility and Diamagnetic Anisotropy of Membrane Proteins from Structural Subunits

The behavior of large, complex molecules in the presence of magnetic fields is experimentally challenging to measure and computationally intensive to predict. This work proposes a novel, mixed-methods approach for efficiently computing the principal magnetic susceptibilities and diamagnetic anisotropy of membrane proteins. The hierarchical primary (amino acid), secondary (α helical and β sheet), and tertiary (α helix and β barrel) structure of transmembrane proteins enables analysis of a complex molecule using discrete subunits of varying size and resolution. The proposed method converts the magnetic susceptibility tensor for all protein subunits to a unit coordinate system and sums them to build the magnetic susceptibility tensor for the membrane protein. Using this approach, we calculate the diamagnetic anisotropy for all transmembrane proteins of known structure and investigate the effect of different subunit resolutions on the resulting predictions of diamagnetic anisotropy. We demonstrate that amino acid residues with aromatic side groups exhibit higher diamagnetic anisotropies. On average, high percentages of aromatic amino acid subunits, a β barrel tertiary structure, and a small volume are correlated with high volumetric diamagnetic anisotropy. Finally, we demonstrate that accounting for the spatial position of the residues with respect to one another is critical to accurately computing the magnetic properties of the complex protein molecule.

On the calculation of second-order magnetic properties using subsystem approaches in a relativistic framework

The implementation of second-order magnetic properties in a frozen density embedding scheme in a four component relativistic framework is outlined and applied to model H₂X–H₂O systems (X = Se, Te, Po).

Automated assignment of NMR chemical shifts based on a known structure and 4D spectra

Apart from their central role during 3D structure determination of proteins the backbone chemical shift assignment is the basis for a number of applications, like chemical shift perturbation mapping and studies on the dynamics of proteins. This assignment is not a trivial task even if a 3D protein structure is known and needs almost as much effort as the assignment for structure prediction if performed manually. We present here a new algorithm based solely on 4D [(1)H,(15)N]-HSQC-NOESY-[(1)H,(15)N]-HSQC spectra which is able to assign a large percentage of chemical shifts (73-82 %) unambiguously, demonstrated with proteins up to a size of 250 residues. For the remaining residues, a small number of possible assignments is filtered out. This is done by comparing distances in the 3D structure to restraints obtained from the peak volumes in the 4D spectrum. Using dead-end elimination, assignments are removed in which at least one of the restraints is violated. Including additional information from chemical shift predictions, a complete unambiguous assignment was obtained for Ubiquitin and 95 % of the residues were correctly assigned in the 251 residue-long N-terminal domain of enzyme I. The program including source code is available at https://github.com/thomasexner/4Dassign .

Demystifying fluorine chemical shifts: electronic structure calculations address origins of seemingly anomalous (19)F-NMR spectra of fluorohistidine isomers and analogues

Fluorine NMR spectroscopy is a powerful tool for studying biomolecular structure, dynamics, and ligand binding, yet the origins of (19)F chemical shifts are not well understood. Herein, we use electronic structure calculations to describe the changes in (19)F chemical shifts of 2F- and 4F-histidine/(5-methyl)-imidazole upon acid titration. While the protonation of the 2F species results in a deshielded chemical shift, protonation of the 4F isomer results in an opposite, shielded chemical shift. The deshielding of 2F-histidine/(5-methyl)-imidazole upon protonation can be rationalized by concomitant decreases in charge density on fluorine and a reduced dipole moment. These correlations do not hold for 4F-histidine/(5-methyl)-imidazole, however. Molecular orbital calculations reveal that for the 4F species, there are no lone pair electrons on the fluorine until protonation. Analysis of a series of 4F-imidazole analogues, all with delocalized fluorine electron density, indicates that the deshielding of (19)F chemical shifts through substituent effects correlates with increased C-F bond polarity. In summary, the delocalization of fluorine electrons in the neutral 4F species, with gain of a lone pair upon protonation may help explain the difficulty in developing a predictive framework for fluorine chemical shifts. Ideas debated by chemists over 40 years ago, regarding fluorine's complex electronic effects, are shown to have relevance for understanding and predicting fluorine NMR spectra.

Quantum calculation of protein NMR chemical shifts based on the automated fragmentation method

The performance of quantum mechanical methods on the calculation of protein NMR chemical shifts is reviewed based on the recently developed automatic fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) approach. By using the Poisson-Boltzmann (PB) model and first solvation water molecules, the influence of solvent effect is also discussed. Benefiting from the fragmentation algorithm, the AF-QM/MM approach is computationally efficient, linear-scaling with a low pre-factor, and thus can be applied to routinely calculate the ab initio NMR chemical shifts for proteins of any size. The results calculated using Density Functional Theory (DFT) show that when the solvent effect is included, this method can accurately reproduce the experimental ¹H NMR chemical shifts, while the ¹³C NMR chemical shifts are less affected by the solvent. However, although the inclusion of solvent effect shows significant improvement for ¹⁵N chemical shifts, the calculated values still have large deviations from the experimental observations. Our study further demonstrates that AF-QM/MM calculated results accurately reflect the dependence of ¹³C(α) NMR chemical shifts on the secondary structure of proteins, and the calculated ¹H chemical shift can be utilized to discriminate the native structure of proteins from decoys.

Composite-180° pulse-based symmetry sequences to recouple proton chemical shift anisotropy tensors under ultrafast MAS solid-state NMR spectroscopy

There is considerable interest in the measurement of proton ((1)H) chemical shift anisotropy (CSA) tensors to obtain deeper insights into H-bonding interactions which find numerous applications in chemical and biological systems. However, the presence of strong (1)H/(1)H dipolar interaction makes it difficult to determine small size (1)H CSAs from the homogeneously broadened NMR spectra. Previously reported pulse sequences for (1)H CSA recoupling are prone to the effects of radio frequency field (B1) inhomogeneity. In the present work we have carried out a systematic study using both numerical and experimental approaches to evaluate γ-encoded radio frequency (RF) pulse sequences based on R-symmetries that recouple (1)H CSA in the indirect dimension of a 2D (1)H/(1)H anisotropic/isotropic chemical shift correlation experiment under ultrafast magic angle spinning (MAS) frequencies. The spectral resolution and sensitivity can be significantly improved in both frequency dimensions of the 2D (1)H/(1)H correlation spectrum without decoupling (1)H/(1)H dipolar couplings but by using ultrafast MAS rates up to 70 kHz. We successfully demonstrate that with a reasonable RF field requirement (<200 kHz) a set of symmetry-based recoupling sequences, with a series of phase-alternating 270°0-90°180 composite-180° pulses, are more robust in combating B1 inhomogeneity effects. In addition, our results show that the new pulse sequences render remarkable (1)H CSA recoupling efficiency and undistorted CSA lineshapes. Experimental results on citric acid and malonic acid comparing the efficiencies of these newly developed pulse sequences with that of previously reported CSA recoupling pulse sequences are also reported under ultrafast MAS conditions.

Accurate ab initio prediction of NMR chemical shifts of nucleic acids and nucleic acids/protein complexes

NMR chemical shift predictions based on empirical methods are nowadays indispensable tools during resonance assignment and 3D structure calculation of proteins. However, owing to the very limited statistical data basis, such methods are still in their infancy in the field of nucleic acids, especially when non-canonical structures and nucleic acid complexes are considered. Here, we present an ab initio approach for predicting proton chemical shifts of arbitrary nucleic acid structures based on state-of-the-art fragment-based quantum chemical calculations. We tested our prediction method on a diverse set of nucleic acid structures including double-stranded DNA, hairpins, DNA/protein complexes and chemically-modified DNA. Overall, our quantum chemical calculations yield highly/very accurate predictions with mean absolute deviations of 0.3–0.6 ppm and correlation coefficients (r²) usually above 0.9. This will allow for identifying misassignments and validating 3D structures. Furthermore, our calculations reveal that chemical shifts of protons involved in hydrogen bonding are predicted significantly less accurately. This is in part caused by insufficient inclusion of solvation effects. However, it also points toward shortcomings of current force fields used for structure determination of nucleic acids. Our quantum chemical calculations could therefore provide input for force field optimization.

Cytochrome-P450-cytochrome-b5 interaction in a membrane environment changes 15N chemical shift anisotropy tensors

It has been well realized that the dependence of chemical shift anisotropy (CSA) tensors on the amino acid sequence, secondary structure, dynamics, and electrostatic interactions can be utilized in the structural and dynamic studies of proteins by NMR spectroscopy. In addition, CSA tensors could also be utilized to measure the structural interactions between proteins in a protein-protein complex. To this end, we report the experimentally measured backbone amide-(15)N CSA tensors for a membrane-bound 16.7 kDa full-length rabbit cytochrome-b5 (cytb5), in complexation with a 55.8 kDa microsomal rabbit cytochrome P450 2B4 (cytP4502B4). The (15)N-CSAs, determined using the (15)N CSA/(15)N-(1)H dipolar coupling transverse cross-correlated rates, for free cytb5 are compared with those for the cytb5 bound to cytP4502B4. An overall increase in backbone amide-(15)N transverse cross-correlated rates for the cytb5 residues in the cytb5-cytP450 complex is observed as compared to the free cytb5 residues. Due to fast spin-spin relaxation (T2) and subsequent broadening of the signals in the complex, we could measure amide-(15)N CSAs only for 48 residues of cytb5 as compared to 84 residues of free cytb5. We observed a change in (15)N CSA for most residues of cytb5 in the complex, as compared to free cytb5, suggesting a dynamic interaction between the oppositely charged surfaces of anionic cytb5 and cationic cytP450. The mean values of (15)N CSA determined for residues in helical, sheet, and turn regions of cytb5 in the complex are -184.5, -146.8, and -146.2 ppm, respectively, with an overall average value of -165.5 ppm (excluding the values from residues in more flexible termini). The measured CSA value for residues in helical conformation is slightly larger as compared to previously reported values. This may be attributed to the paramagnetic effect from Fe(III) of the heme in cytb5, which is similar to our previously reported values for the free cytb5.

NMR crystallography of enzyme active sites: probing chemically detailed, three-dimensional structure in tryptophan synthase

NMR crystallography--the synergistic combination of X-ray diffraction, solid-state NMR spectroscopy, and computational chemistry--offers unprecedented insight into three-dimensional, chemically detailed structure. Initially, researchers used NMR crystallography to refine diffraction data from organic and inorganic solids. Now we are applying this technique to explore active sites in biomolecules, where it reveals chemically rich detail concerning the interactions between enzyme site residues and the reacting substrate. Researchers cannot achieve this level of detail from X-ray, NMR,or computational methodologies in isolation. For example, typical X-ray crystal structures (1.5-2.5 Å resolution) of enzyme-bound intermediates identify possible hydrogen-bonding interactions between site residues and substrate but do not directly identify the protonation states. Solid-state NMR can provide chemical shifts for selected atoms of enzyme-substrate complexes, but without a larger structural framework in which to interpret them only empirical correlations with local chemical structure are possible. Ab initio calculations and molecular mechanics can build models for enzymatic processes, but they rely on researcher-specified chemical details. Together, however, X-ray diffraction, solid-state NMR spectroscopy, and computational chemistry can provide consistent and testable models for structure and function of enzyme active sites: X-ray crystallography provides a coarse framework upon which scientists can develop models of the active site using computational chemistry; they can then distinguish these models by comparing calculated NMR chemical shifts with the results of solid-state NMR spectroscopy experiments. Conceptually, each technique is a puzzle piece offering a generous view of the big picture. Only when correctly pieced together, however, can they reveal the big picture at the highest possible resolution. In this Account, we detail our first steps in the development of NMR crystallography applied to enzyme catalysis. We begin with a brief introduction to NMR crystallography and then define the process that we have employed to probe the active site in the β-subunit of tryptophan synthase with unprecedented atomic-level resolution. This approach has resulted in a novel structural hypothesis for the protonation state of the quinonoid intermediate in tryptophan synthase and its surprising role in directing the next step in the catalysis of L-Trp formation.

Allostery and substrate channeling in the tryptophan synthase bienzyme complex: evidence for two subunit conformations and four quaternary states

The allosteric regulation of substrate channeling in tryptophan synthase involves ligand-mediated allosteric signaling that switches the α- and β-subunits between open (low activity) and closed (high activity) conformations. This switching prevents the escape of the common intermediate, indole, and synchronizes the α- and β-catalytic cycles. (19)F NMR studies of bound α-site substrate analogues, N-(4'-trifluoromethoxybenzoyl)-2-aminoethyl phosphate (F6) and N-(4'-trifluoromethoxybenzenesulfonyl)-2-aminoethyl phosphate (F9), were found to be sensitive NMR probes of β-subunit conformation. Both the internal and external aldimine F6 complexes gave a single bound peak at the same chemical shift, while α-aminoacrylate and quinonoid F6 complexes all gave a different bound peak shifted by +1.07 ppm. The F9 complexes exhibited similar behavior, but with a corresponding shift of -0.12 ppm. X-ray crystal structures show the F6 and F9 CF3 groups located at the α-β subunit interface and report changes in both the ligand conformation and the surrounding protein microenvironment. Ab initio computational modeling suggests that the change in (19)F chemical shift results primarily from changes in the α-site ligand conformation. Structures of α-aminoacrylate F6 and F9 complexes and quinonoid F6 and F9 complexes show the α- and β-subunits have closed conformations wherein access of ligands into the α- and β-sites from solution is blocked. Internal and external aldimine structures show the α- and β-subunits with closed and open global conformations, respectively. These results establish that β-subunits exist in two global conformational states, designated open, where the β-sites are freely accessible to substrates, and closed, where the β-site portal into solution is blocked. Switching between these conformations is critically important for the αβ-catalytic cycle.

Automated Fragmentation QM/MM Calculation of Amide Proton Chemical Shifts in Proteins with Explicit Solvent Model

We have performed a density functional theory (DFT) calculation of the amide proton NMR chemical shift in proteins using a recently developed automated fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) approach. Systematic investigation was carried out to examine the influence of explicit solvent molecules, cooperative hydrogen bonding effects, density functionals, size of the basis sets, and the local geometry of proteins on calculated chemical shifts. Our result demonstrates that the predicted amide proton ((1)HN) NMR chemical shift in explicit solvent shows remarkable improvement over that calculated with the implicit solvation model. The cooperative hydrogen bonding effect is also shown to improve the accuracy of (1)HN chemical shifts. Furthermore, we found that the OPBE exchange-correlation functional is the best density functional for the prediction of protein (1)HN chemical shifts among a selective set of DFT methods (namely, B3LYP, B3PW91, M062X, M06L, mPW1PW91, OB98, OPBE), and the locally dense basis set of 6-311++G**/4-31G* is shown to be sufficient for (1)HN chemical shift calculation. By taking ensemble averaging into account, (1)HN chemical shifts calculated by the AF-QM/MM approach can be used to validate the performance of various force fields. Our study underscores that the electronic polarization of protein is of critical importance to stabilizing hydrogen bonding, and the AF-QM/MM method is able to describe the local chemical environment in proteins more accurately than most widely used empirical models.

Site-specific perturbations of alpha-synuclein fibril structure by the Parkinson's disease associated mutations A53T and E46K

Parkinson's disease (PD) is pathologically characterized by the presence of Lewy bodies (LBs) in dopaminergic neurons of the substantia nigra. These intracellular inclusions are largely composed of misfolded α-synuclein (AS), a neuronal protein that is abundant in the vertebrate brain. Point mutations in AS are associated with rare, early-onset forms of PD, although aggregation of the wild-type (WT) protein is observed in the more common sporadic forms of the disease. Here, we employed multidimensional solid-state NMR experiments to assess A53T and E46K mutant fibrils, in comparison to our recent description of WT AS fibrils. We made de novo chemical shift assignments for the mutants, and used these chemical shifts to empirically determine secondary structures. We observe significant perturbations in secondary structure throughout the fibril core for the E46K fibril, while the A53T fibril exhibits more localized perturbations near the mutation site. Overall, these results demonstrate that the secondary structure of A53T has some small differences from the WT and the secondary structure of E46K has significant differences, which may alter the overall structural arrangement of the fibrils.

Analysis of the contributions of ring current and electric field effects to the chemical shifts of RNA bases

Ring current and electric field effects can considerably influence NMR chemical shifts in biomolecules. Understanding such effects is particularly important for the development of accurate mappings between chemical shifts and the structures of nucleic acids. In this work, we first analyzed the Pople and the Haigh-Mallion models in terms of their ability to describe nitrogen base conjugated ring effects. We then created a database (DiBaseRNA) of three-dimensional arrangements of RNA base pairs from X-ray structures, calculated the corresponding chemical shifts via a hybrid density functional theory approach and used the results to parametrize the ring current and electric field effects in RNA bases. Next, we studied the coupling of the electric field and ring current effects for different inter-ring arrangements found in RNA bases using linear model fitting, with joint electric field and ring current, as well as only electric field and only ring current approximations. Taken together, our results provide a characterization of the interdependence of ring current and electric field geometric factors, which is shown to be especially important for the chemical shifts of non-hydrogen atoms in RNA bases.

Quantum chemical calculations of amide-15N chemical shift anisotropy tensors for a membrane-bound cytochrome-b5

There is considerable interest in determining amide-(15)N chemical shift anisotropy (CSA) tensors from biomolecules and understanding their variation for structural and dynamics studies using solution and solid-state NMR spectroscopy and also by quantum chemical calculations. Due to the difficulties associated with the measurement of CSA tensors from membrane proteins, NMR-based structural studies heavily relied on the CSA tensors determined from model systems, typically single crystals of model peptides. In the present study, the principal components of backbone amide-(15)N CSA tensors have been determined using density functional theory for a 16.7 kDa membrane-bound paramagnetic heme containing protein, cytochrome-b(5) (cytb(5)). All the calculations were performed by taking residues within 5 Å distance from the backbone amide-(15)N nucleus of interest. The calculated amide-(15)N CSA spans agree less well with our solution NMR data determined for an effective internuclear distance r(N-H) = 1.023 Å and a constant angle β = 18° that the least shielded component (δ(11)) makes with the N-H bond. The variation of amide-(15)N CSA span obtained using quantum chemical calculations is found to be smaller than that obtained from solution NMR measurements, whereas the trends of the variations are found to be in close agreement. We believe that the results reported in this study will be useful in studying the structure and dynamics of membrane proteins and heme-containing proteins, and also membrane-bound protein-protein complexes such as cytochromes-b5-P450.

Toward the Quantum Chemical Calculation of NMR Chemical Shifts of Proteins. 3. Conformational Sampling and Explicit Solvents Model

Fragment-based quantum chemical calculations are able to accurately calculate NMR chemical shifts even for very large molecules like proteins. But even with systematic optimization of the level of theory and basis sets as well as the use of implicit solvents models, some nuclei like polar protons and nitrogens suffer from poor predictions. Two properties of the real system, strongly influencing the experimental chemical shifts but almost always neglected in the calculations, will be discussed here in great detail: (1) conformational averaging and (2) interactions with first-shell solvent molecules. Classical molecular dynamics simulations in explicit water were carried out for obtaining a representative ensemble including the arrangement of neighboring solvent molecules, which was then subjected to quantum chemical calculations. We could demonstrate with the small test system N-methyl acetamide (NMA) that the calculated chemical shifts show immense variations of up to 6 ppm and 50 ppm for protons and nitrogens, respectively, depending on the snapshot taken from a classical molecular dynamics simulation. Applying the same approach to the HA2 domain of the influenza virus glycoprotein hemagglutinin, a 32-amino-acid-long polypeptide, and comparing averaged values to the experiment, chemical shifts of nonpolar protons and carbon atoms in proteins were calculated with unprecedented accuracy. Additionally, the mean absolute error could be reduced by a factor of 2.43 for polar protons, and reasonable correlations were obtained for nitrogen and carbonyl carbon in contrast to all other studies published so far.

Assessing and refining molecular dynamics simulations of proteins with nuclear magnetic resonance data

The sophistication of the force fields, algorithms and hardware used for molecular dynamics (MD) simulations of proteins is continuously increasing. No matter how advanced the methodology, however, it is essential to evaluate the appropriateness of the structures sampled in a simulation by comparison with quantitative experimental data. Solution nuclear magnetic resonance (NMR) data are particularly useful for checking the quality of protein simulations, as they provide both structural and dynamic information on a variety of temporal and spatial scales. Here, various features and implications of using NMR data to validate and bias MD simulations are outlined, including an overview of the different types of NMR data that report directly on structural properties and of relevant simulation techniques. The focus throughout is on how to properly account for conformational averaging, particularly within the context of the assumptions inherent in the relationships that link NMR data to structural properties.

Determination of 15N chemical shift anisotropy from a membrane-bound protein by NMR spectroscopy

Chemical shift anisotropy (CSA) tensors are essential in the structural and dynamic studies of proteins using NMR spectroscopy. Results from relaxation studies in biomolecular solution and solid-state NMR experiments on aligned samples are routinely interpreted using well-characterized CSA tensors determined from model compounds. Since CSA tensors, particularly the (15)N CSA, highly depend on a number of parameters including secondary structure, electrostatic interaction, and the amino acid sequence, there is a need for accurately determined CSA tensors from proteins. In this study, we report the backbone amide-(15)N CSA tensors for a 16.7-kDa membrane-bound and paramagnetic-heme containing protein, rabbit Cytochrome b(5) (cytb(5)), determined using the (15)N CSA/(15)N-(1)H dipolar transverse cross-correlation rates. The mean values of (15)N CSA determined for residues in helical, sheet, and turn regions are -187.9, -166.0, and -161.1 ppm, respectively, with an overall average value of -171.7 ppm. While the average CSA value determined from this study is in good agreement with previous solution NMR experiments on small globular proteins, the CSA value determined for residues in helical conformation is slightly larger, which may be attributed to the paramagnetic effect from Fe(III) of the heme unit in cytb(5). However, like in previous solution NMR studies, the CSA values reported in this study are larger than the values measured from solid-state NMR experiments. We believe that the CSA parameters reported in this study will be useful in determining the structure, dynamics, and orientation of proteins, including membrane proteins, using NMR spectroscopy.

Conformational distributions of denatured and unstructured proteins are similar to those of 20 × 20 blocked dipeptides

Understanding intrinsic conformational preferences of amino-acids in unfolded proteins is important for elucidating the underlying principles of their stability and re-folding on biological timescales. Here, to investigate the neighbor interaction effects on the conformational propensities of amino-acids, we carried out (1)H NMR experiments for a comprehensive set of blocked dipeptides and measured the scalar coupling constants between alpha protons and amide protons as well as their chemical shifts. Detailed inspection of these NMR properties shows that, irrespective of amino-acid side-chain properties, the distributions of the measured coupling constants and chemical shifts of the dipeptides are comparatively narrow, indicating small variances of their conformation distributions. They are further compared with those of blocked amino-acids (Ac-X-NHMe), oligopeptides (Ac-GGXGG-NH(2)), and native (lysozyme), denatured (lysozyme and outer membrane protein X from Escherichia coli), unstructured (Domain 2 of the protein 5A of Hepatitis C virus), and intrinsically disordered (hNlg3cyt: intracellular domain of human NL3) proteins. These comparative investigations suggest that the conformational preferences and local solvation environments of the blocked dipeptides are quite similar to not only those of other short oligopeptides but also those of denatured and natively unfolded proteins.

Toward the Quantum Chemical Calculation of NMR Chemical Shifts of Proteins. 2. Level of Theory, Basis Set, and Solvents Model Dependence

It has been demonstrated that the fragmentation scheme of our adjustable density matrix assembler (ADMA) approach for the quantum chemical calculations of very large systems is well-suited to calculate NMR chemical shifts of proteins [ Frank et al. Proteins2011, 79, 2189-2202 ]. The systematic investigation performed here on the influences of the level of theory, basis set size, inclusion or exclusion of an implicit solvent model, and the use of partial charges to describe additional parts of the macromolecule on the accuracy of NMR chemical shifts demonstrates that using a valence triple-ζ basis set leads to large improvement compared to the results given in the previous publication. Additionally, moving from the B3LYP to the mPW1PW91 density functional and including partial charges and implicit solvents gave the best results with mean absolute errors of 0.44 ppm for hydrogen atoms excluding H(N) atoms and between 1.53 and 3.44 ppm for carbon atoms depending on the size and also on the accuracy of the protein structure. Polar hydrogen and nitrogen atoms are more difficult to predict. For the first, explicit hydrogen bonds to the solvents need to be included and, for the latter, going beyond DFT to post-Hartree-Fock methods like MP2 is probably required. Even if empirical methods like SHIFTX+ show similar performance, our calculations give for the first time very reliable chemical shifts that can also be used for complexes of proteins with small-molecule ligands or DNA/RNA. Therefore, taking advantage of its ab initio nature, our approach opens new fields of application that would otherwise be largely inaccessible due to insufficient availability of data for empirical parametrization.

Mutant protein A30P α-synuclein adopts wild-type fibril structure, despite slower fibrillation kinetics

α-Synuclein (AS) is associated with both sporadic and familial forms of Parkinson disease (PD). In sporadic disease, wild-type AS fibrillates and accumulates as Lewy bodies within dopaminergic neurons of the substantia nigra. The accumulation of misfolded AS is associated with the death of these neurons, which underlies many of the clinical features of PD. In addition, a rare missense mutation in AS, A30P, is associated with highly penetrant, autosomal dominant PD, although the pathogenic mechanism is unclear. A30P AS fibrillates more slowly than the wild-type (WT) protein in vitro and has been reported to preferentially adopt a soluble, protofibrillar conformation. This has led to speculation that A30P forms aggregates that are distinct in structure compared with wild-type AS. Here, we perform a detailed comparison of the chemical shifts and secondary structures of these fibrillar species, based upon our recent characterization of full-length WT fibrils. We have assigned A30P AS fibril chemical shifts de novo and used them to determine its secondary structure empirically. Our results illustrate that although A30P forms fibrils more slowly than WT in vitro, the chemical shifts and secondary structure of the resultant fibrils are in high agreement, demonstrating a conserved β-sheet core.

Structure, dynamics, lipid binding, and physiological relevance of the putative GTPase-binding domain of Dictyostelium formin C

Dictyostelium Formin C (ForC) is involved in the regulation of local actin cytoskeleton reorganization (e.g. during cellular adhesion or migration). ForC contains formin homology 2 and 3 (FH2 and -3) domains and an N-terminal putative GTPase-binding domain (GBD) but lacks a canonical FH1 region. To better understand the role of the GBD, its structure, dynamics, lipid-binding properties, and cellular functions were analyzed by NMR and CD spectroscopy and by in vivo fluorescence microscopy. Moreover, the program CS-Rosetta was tested for the structure prediction based on chemical shift data only. The ForC GBD adopts an ubiquitin-like α/β-roll fold with an unusually long loop between β-strands 1 and 2. Based on the lipid-binding data, the presence of DPC micelles induces the formation of α-helical secondary structure and a rearrangement of the tertiary structure. Lipid-binding studies with a mutant protein and a peptide suggest that the β1-β2 loop is not relevant for these conformational changes. Whereas small amounts of negatively charged phosphoinositides (1,2-dioctanoyl-sn-glycero-3-(phosphoinositol 4,5-bisphosphate) and 1,2-dihexanoyl-sn-glycero-3-(phosphoinositol 3,4,5-trisphosphate)) lower the micelle concentration necessary to induce the observed spectral changes, other negatively charged phospholipids (1,2-dihexanoyl-sn-glycero-3-(phospho-L-serine) and 1,2-dihexanoyl-sn-glycero-3-phospho-(1'-rac-glycerol)) had no such effect. Interestingly, bicelles and micelles composed of diacylphosphocholines had no effect on the GBD structure. Our data suggest a model in which part of the large positively charged surface area of the GBD mediates localization to specific membrane patches, thereby regulating interactions with signaling proteins. Our cellular localization studies show that both the GBD and the FH3 domain are required for ForC targeting to cell-cell contacts and early phagocytic cups and macropinosomes.

Toward the quantum chemical calculation of nuclear magnetic resonance chemical shifts of proteins

Despite the many protein structures solved successfully by nuclear magnetic resonance (NMR) spectroscopy, quality control of NMR structures is still by far not as well established and standardized as in crystallography. Therefore, there is still the need for new, independent, and unbiased evaluation tools to identify problematic parts and in the best case also to give guidelines that how to fix them. We present here, quantum chemical calculations of NMR chemical shifts for many proteins based on our fragment-based quantum chemical method: the adjustable density matrix assembler (ADMA). These results show that (13)C chemical shifts of reasonable accuracy can be obtained that can already provide a powerful measure for the structure validation. (1)H and even more (15)N chemical shifts deviate more strongly from experiment due to the insufficient treatment of solvent effects and conformational averaging.

Analysis of organic compounds in an urban wastewater treatment plant effluent

In the present work we fractioned the effluent water from an urban sewage treatment plant (USTP) of Ribarroja (Valencia, Spain) using the conventional protocol based on DAX8 and XAD4 resins. The fractions were analyzed by elemental analysis, FT-IR, 1H-NMR, COSY-NMR, HSQC-NMR, FAB+-MS and also by derivatization with bis(trimethylsilyl)trifluoroacetamide with 10% of trimethylchlorosilane. The four fractions obtained have common spectroscopic features and individual compounds indicating that the fractioning procedure is inefficient at separating different families of compounds. Gas chromatography/mass spectrometry (GC-MS) analysis of the derivatized fractions allowed identification of many individual compounds. The main classes of organic compounds present in the effluent are saccharides, amino acids, fatty acids, hydroxyacids, aromatic compounds and steroids. Also we were able to identify in the effluent the emerging pollutants paracetamol and ketoprofen - two best-selling antiinflammatory drugs used in humans.

Protegrin-1 orientation and physicochemical properties in membrane bilayers studied by potential of mean force calculations

Protegrin-1 (PG-1) belongs to the family of antimicrobial peptides. It interacts specifically with the membrane of a pathogen and kills the pathogen by releasing its cellular contents. To fully understand the energetics governing the orientation of PG-1 in different membrane environments and its effects on the physicochemical properties of the peptide and membrane bilayers, we have performed the potential of mean force (PMF) calculations as a function of its tilt angle at four distinct rotation angles in explicit membranes composed of either DLPC (1,2-dilauroylphosphatidylcholine) or POPC (1-palmitoyl-2-oleoylphosphatidylcholine) lipid molecules. The resulting PMFs in explicit lipid bilayers were then used to search for the optimal hydrophobic thickness of the EEF1/IMM1 implicit membrane model in which a two-dimensional PMF in the tilt and rotation space was calculated. The PMFs in explicit membrane systems clearly reveal that the energetically favorable tilt angle is affected by both the membrane hydrophobic thickness and the PG-1 rotation angle. Local thinning of the membrane around PG-1 is observed upon PG-1 tilting. The thinning is caused by both hydrophobic mismatch and arginine-lipid head group interactions. The two-dimensional PMF in the implicit membrane is in good accordance with those from the explicit membrane simulations. The ensemble-averaged Val16 (15)N and (13)CO chemical shifts weighted by the two-dimensional PMF agree fairly well with the experimental values, suggesting the importance of peptide dynamics in calculating such ensemble properties for direct comparison with experimental observables.

Solid-state NMR techniques for the structural determination of amyloid fibrils

This review discusses the solid-state NMR techniques developed for the study of amyloid fibrils. Literature up to the end of 2010 has been surveyed and the materials are organized according to five categories, viz. homonuclear dipolar recoupling and polarization transfer via J-coupling, heteronuclear dipolar recoupling, correlation spectroscopy, recoupling of chemical shift anisotropy, and tensor correlation. Our emphasis is on the NMR techniques and their practical aspects. The biological implications of the results obtained for amyloid fibrils are only briefly discussed. Our main objective is to showcase the power of NMR in the study of biological unoriented solids.

Radiofrequency quadrupolar NMR stark spectroscopy: steady state response calibration and tensorial mapping

Radiofrequency electric (E) fields oscillating at twice the usual NMR frequency (2ω(0)) can induce double-quantum transitions in quadrupolar nuclei, an NMR Stark effect. Characterization of such is of interest to aid understanding of electrostatic effects in NMR spectra. Calibration of Stark responses to an applied electric field may also be used to assess native fields within molecules and materials. We present high-field (14.1 T), room-temperature NMR experiments to calibrate the 2ω(0) Stark response in crystalline GaAs. This system presents an important test of current techniques and conditions, as historical studies at low field (500-900 mT) and low temperature (77 K) provide a basis for comparison. Our measurements of steady state response reveal the quadrupolar Stark tuning rate for (69)Ga in this material. The value, β(Q) = (11.5 ± 0.1) × 10(12) m(-1), is 3.6 times larger than the most-reliable prior result. In the process, we also uncovered a previously unobserved double-quantum steady state coherence. It appears as a completely separable dispersive signal component in quadrature-detected presaturation spectra versus offset from 2ω(0). The new component may eventually afford an independent route to calibrating β(Q). Finally, we demonstrated exceptional agreement with theory of the orientation-dependent Stark response for rotation of the sample relative to B(0) over a range of 90° and for E-field amplitudes from 30-180 V/cm.

Assignment strategies for large proteins by magic-angle spinning NMR: the 21-kDa disulfide-bond-forming enzyme DsbA

We present strategies for chemical shift assignments of large proteins by magic-angle spinning solid-state NMR, using the 21-kDa disulfide-bond-forming enzyme DsbA as prototype. Previous studies have demonstrated that complete de novo assignments are possible for proteins up to approximately 17 kDa, and partial assignments have been performed for several larger proteins. Here we show that combinations of isotopic labeling strategies, high field correlation spectroscopy, and three-dimensional (3D) and four-dimensional (4D) backbone correlation experiments yield highly confident assignments for more than 90% of backbone resonances in DsbA. Samples were prepared as nanocrystalline precipitates by a dialysis procedure, resulting in heterogeneous linewidths below 0.2 ppm. Thus, high magnetic fields, selective decoupling pulse sequences, and sparse isotopic labeling all improved spectral resolution. Assignments by amino acid type were facilitated by particular combinations of pulse sequences and isotopic labeling; for example, transferred echo double resonance experiments enhanced sensitivity for Pro and Gly residues; [2-(13)C]glycerol labeling clarified Val, Ile, and Leu assignments; in-phase anti-phase correlation spectra enabled interpretation of otherwise crowded Glx/Asx side-chain regions; and 3D NCACX experiments on [2-(13)C]glycerol samples provided unique sets of aromatic (Phe, Tyr, and Trp) correlations. Together with high-sensitivity CANCOCA 4D experiments and CANCOCX 3D experiments, unambiguous backbone walks could be performed throughout the majority of the sequence. At 189 residues, DsbA represents the largest monomeric unit for which essentially complete solid-state NMR assignments have so far been achieved. These results will facilitate studies of nanocrystalline DsbA structure and dynamics and will enable analysis of its 41-kDa covalent complex with the membrane protein DsbB, for which we demonstrate a high-resolution two-dimensional (13)C-(13)C spectrum.