Accurate measurement of methyl 13C chemical shifts by solid-state NMR for the determination of protein side chain conformation: the influenza a M2 transmembrane peptide as an example

Profiling Enzyme Activity of l-Asparaginase II by NMR-Based Methyl Fingerprinting at Natural Abundance

l-asparaginase II (MW 135 kDa) from E. coli is an FDA-approved protein drug used for the treatment of childhood leukemia. Despite its long history as a chemotherapeutic, the structural basis of enzyme action, in solution, remains widely contested. In this work, methyl-based 2D [¹H-¹³C]-heteronuclear single-quantum correlation (HSQC) NMR, at natural abundance, has been used to profile the enzymatic activity of the commercially available enzyme drug. The [¹H-¹³C]-HSQC NMR spectra of the protein reveal the role of a flexible loop segment in the activity of the enzyme, in solution. Addition of asparagine to the protein results in distinct conformational changes of the loop that could be signatures of intermediates formed in the catalytic reaction. To this end, an isothermal titration calorimetry (ITC)-based assay has been developed to measure the enzymatic reaction enthalpy, as a marker for its activity. Combining both ITC and NMR, it was shown that the disruption of the protein conformation can result in the loss of function. The scope, robustness, and validity of the loop fingerprints in relation to enzyme activity have been tested under different solution conditions. Overall, our results indicate that 2D NMR can be used reliably to gauge the structure-function of this enzyme, bypassing the need to label the protein. Such natural abundant NMR methods can be potentially extended to probe the structure-function aspects of high-molecular-weight protein therapeutics (glycosylated protein drugs, enzymes, therapeutic monoclonal antibodies, antibody-drug conjugates, and Fc-fusion proteins), where (a) flexible loops are required for their function and (b) isotope labeling may not be straightforward.

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.

Hydration and Dynamics of Full-Length Tau Amyloid Fibrils Investigated by Solid-State Nuclear Magnetic Resonance

The microtubule-associated protein tau aggregates into distinct neurofibrillary tangles in brains afflicted with multiple neurodegenerative diseases such as Alzheimer's disease and corticobasal degeneration (CBD). The mechanism of tau misfolding and aggregation is poorly understood. Determining the structure, dynamics, and water accessibility of tau filaments may provide insight into the pathway of tau misfolding. Here, we investigate the hydration and dynamics of the β-sheet core of heparin-fibrillized 0N4R tau using solid-state nuclear magnetic resonance spectroscopy. This β-sheet core consists of the second and third microtubule-binding repeats, R2 and R3, respectively, which form a hairpin. Water-edited two-dimensional (2D) 13C-13C and 15N-13C correlation spectra show that most residues in R2 and R3 domains have low water accessibility, indicating that this hairpin is surrounded by other proteinaceous segments. However, a small number of residues, especially S285 and S316, are well hydrated compared to other Ser and Thr residues, suggesting that there is a small water channel in the middle of the hairpin. To probe whether water accessibility correlates with protein dynamics, we measured the backbone N-H dipolar couplings of the β-sheet core. Interestingly, residues in the fourth microtubule-binding repeat, R4, show rigid-limit N-H dipolar couplings, even though this domain exhibits weaker intensities in the 2D 15N-13C correlation spectra. These results suggest that the R4 domain participates in cross-β hydrogen bonding in some of the subunits but exhibits dynamic disorder in other subunits. Taken together, these hydration and dynamics data indicate that the R2-R3 hairpin of 0N4R tau is shielded from water by other proteinaceous segments on the exterior but contains a small water pore in the interior. This structural topology has various similarities with the CBD tau fibril structure but also shows specific differences. The disorder of the R4 domain and the presence of a small water channel in the heparin-fibrillized 4R tau have implications for the structure of tau fibrils in diseased brains.

Chemical shift-based methods in NMR structure determination

Chemical shifts are highly sensitive probes harnessed by NMR spectroscopists and structural biologists as conformational parameters to characterize a range of biological molecules. Traditionally, assignment of chemical shifts has been a labor-intensive process requiring numerous samples and a suite of multidimensional experiments. Over the past two decades, the development of complementary computational approaches has bolstered the analysis, interpretation and utilization of chemical shifts for elucidation of high resolution protein and nucleic acid structures. Here, we review the development and application of chemical shift-based methods for structure determination with a focus on ab initio fragment assembly, comparative modeling, oligomeric systems, and automated assignment methods. Throughout our discussion, we point out practical uses, as well as advantages and caveats, of using chemical shifts in structure modeling. We additionally highlight (i) hybrid methods that employ chemical shifts with other types of NMR restraints (residual dipolar couplings, paramagnetic relaxation enhancements and pseudocontact shifts) that allow for improved accuracy and resolution of generated 3D structures, (ii) the utilization of chemical shifts to model the structures of sparsely populated excited states, and (iii) modeling of sidechain conformations. Finally, we briefly discuss the advantages of contemporary methods that employ sparse NMR data recorded using site-specific isotope labeling schemes for chemical shift-driven structure determination of larger molecules. With this review, we aim to emphasize the accessibility and versatility of chemical shifts for structure determination of challenging biological systems, and to point out emerging areas of development that lead us towards the next generation of tools.

CryoEM Structure Refinement by Integrating NMR Chemical Shifts with Molecular Dynamics Simulations

Single particle cryoEM has emerged as a powerful method for structure determination of proteins and complexes, complementing X-ray crystallography and NMR spectroscopy. Yet, for many systems, the resolution of cryoEM density map has been limited to 4-6 Å, which only allows for resolving bulky amino acids side chains, thus hindering accurate model building from the density map. On the other hand, experimental chemical shifts (CS) from solution and solid state MAS NMR spectra provide atomic level data for each amino acid within a molecule or a complex; however, structure determination of large complexes and assemblies based on NMR data alone remains challenging. Here, we present a novel integrated strategy to combine the highly complementary experimental data from cryoEM and NMR computationally by molecular dynamics simulations to derive an atomistic model, which is not attainable by either approach alone. We use the HIV-1 capsid protein (CA) C-terminal domain as well as the large capsid assembly to demonstrate the feasibility of this approach, termed NMR CS-biased cryoEM structure refinement.

Access to side-chain carbon information in deuterated solids under fast MAS through non-rotor-synchronized mixing

MOCCA provides comprehensive solid-state NMR side chain carbon correlations despite perdeuteration and fast sample rotation, thereby inducing minimal power dissipation.

Magic angle spinning NMR of viruses

Viruses, relatively simple pathogens, are able to replicate in many living organisms and to adapt to various environments. Conventional atomic-resolution structural biology techniques, X-ray crystallography and solution NMR spectroscopy provided abundant information on the structures of individual proteins and nucleic acids comprising viruses; however, viral assemblies are not amenable to analysis by these techniques because of their large size, insolubility, and inherent lack of long-range order. In this article, we review the recent advances in magic angle spinning NMR spectroscopy that enabled atomic-resolution analysis of structure and dynamics of large viral systems and give examples of several exciting case studies.

Practical use of chemical shift databases for protein solid-state NMR: 2D chemical shift maps and amino-acid assignment with secondary-structure information

We introduce a Python-based program that utilizes the large database of (13)C and (15)N chemical shifts in the Biological Magnetic Resonance Bank to rapidly predict the amino acid type and secondary structure from correlated chemical shifts. The program, called PACSYlite Unified Query (PLUQ), is designed to help assign peaks obtained from 2D (13)C-(13)C, (15)N-(13)C, or 3D (15)N-(13)C-(13)C magic-angle-spinning correlation spectra. We show secondary-structure specific 2D (13)C-(13)C correlation maps of all twenty amino acids, constructed from a chemical shift database of 262,209 residues. The maps reveal interesting conformation-dependent chemical shift distributions and facilitate searching of correlation peaks during amino-acid type assignment. Based on these correlations, PLUQ outputs the most likely amino acid types and the associated secondary structures from inputs of experimental chemical shifts. We test the assignment accuracy using four high-quality protein structures. Based on only the Cα and Cβ chemical shifts, the highest-ranked PLUQ assignments were 40-60 % correct in both the amino-acid type and the secondary structure. For three input chemical shifts (CO-Cα-Cβ or N-Cα-Cβ), the first-ranked assignments were correct for 60 % of the residues, while within the top three predictions, the correct assignments were found for 80 % of the residues. PLUQ and the chemical shift maps are expected to be useful at the first stage of sequential assignment, for combination with automated sequential assignment programs, and for highly disordered proteins for which secondary structure analysis is the main goal of structure determination.

A straightforward method for stereospecific assignment of val and leu prochiral methyl groups by solid-state NMR: Scrambling in the [2-13C]Glucose labeling scheme

The unambiguous stereospecific assignment of the prochiral methyl groups in Val and Leu plays an important role in the structural investigation of proteins by NMR. Here, we present a straightforward method for their stereospecific solid-state NMR assignment based on [2-(13)C]Glucose ([2-(13)C]Glc) as the sole carbon source during protein expression. The approach is fundamentally based on the stereo-selective biosynthetic pathway of Val and Leu, and the co-presence of [2-(13)C]pyruvate produced mainly by glycolysis and [3-(13)C]/[1,3-(13)C]pyruvate most probably formed through scrambling in the pentose phosphate pathway. As a consequence, the isotope spin pairs (13)Cβ-(13)Cγ2 and (13)Cα-(13)Cγ1 in Val, and (13)Cγ-(13)Cδ2 and (13)Cβ-(13)Cδ1 in Leu are obtained. The approach is successfully demonstrated with the stereospecific assignment of the methyl groups of Val and Leu of type 3 secretion system PrgI needles and microcrystalline ubiquitin.

Structural basis for proton conduction and inhibition by the influenza M2 protein

The influenza M2 protein forms an acid-activated and drug-sensitive proton channel in the virus envelope that is important for the virus lifecycle. The functional properties and high-resolution structures of this proton channel have been extensively studied to understand the mechanisms of proton conduction and drug inhibition. We review biochemical and electrophysiological studies of M2 and discuss how high-resolution structures have transformed our understanding of this proton channel. Comparison of structures obtained in different membrane-mimetic solvents and under different pH using X-ray crystallography, solution NMR, and solid-state NMR spectroscopy revealed how the M2 structure depends on the environment and showed that the pharmacologically relevant drug-binding site lies in the transmembrane (TM) pore. Competing models of proton conduction have been evaluated using biochemical experiments, high-resolution structural methods, and computational modeling. These results are converging to a model in which a histidine residue in the TM domain mediates proton relay with water, aided by microsecond conformational dynamics of the imidazole ring. These mechanistic insights are guiding the design of new inhibitors that target drug-resistant M2 variants and may be relevant for other proton channels.

Intramolecular 1H-13C distance measurement in uniformly 13C, 15N labeled peptides by solid-state NMR

A (1)H-(13)C frequency-selective REDOR (FS-REDOR) experiment is developed for measuring intramolecular (1)H-(13)C distances in uniformly (13)C, (15)N-labeled molecules. Theory and simulations show that the experiment removes the interfering homonuclear (1)H-(1)H, (13)C-(13)C and heteronuclear (1)H-(15)N, (13)C-(15)N dipolar interactions while retaining the desired heteronuclear (1)H-(13)C dipolar interaction. Our results indicate that this technique, combined with the numerical fitting, can be used to measure a (1)H-(13)C distance up to 5Å. We also demonstrate that the measured intramolecular (1)H-(13)C distances are useful to determine dihedral angles in proteins.

Structured regions of α-synuclein fibrils include the early-onset Parkinson's disease mutation sites

α-Synuclein (AS) fibrils are the major component of Lewy bodies, the pathological hallmark of Parkinson's disease (PD). Here, we use results from an extensive investigation employing solid-state NMR to present a detailed structural characterization and conformational dynamics quantification of full-length AS fibrils. Our results show that the core extends with a repeated structural motif. This result disagrees with the previously proposed fold of AS fibrils obtained with limited solid-state NMR data. Additionally, our results demonstrate that the three single point mutations associated with early-onset PD-A30P, E46K and A53T-are located in structured regions. We find that E46K and A53T mutations, located in rigid β-strands of the wild-type fibrils, are associated with major and minor structural perturbations, respectively.

Structure-based prediction of methyl chemical shifts in proteins

Protein methyl groups have recently been the subject of much attention in NMR spectroscopy because of the opportunities that they provide to obtain information about the structure and dynamics of proteins and protein complexes. With the advent of selective labeling schemes, methyl groups are particularly interesting in the context of chemical shift based protein structure determination, an approach that to date has exploited primarily the mapping between protein structures and backbone chemical shifts. In order to extend the scope of chemical shifts for structure determination, we present here the CH3Shift method of performing structure-based predictions of methyl chemical shifts. The terms considered in the predictions take account of ring current, magnetic anisotropy, electric field, rotameric type, and dihedral angle effects, which are considered in conjunction with polynomial functions of interatomic distances. We show that the CH3Shift method achieves an accuracy in the predictions that ranges from 0.133 to 0.198 ppm for (1)H chemical shifts for Ala, Thr, Val, Leu and Ile methyl groups. We illustrate the use of the method by assessing the accuracy of side-chain structures in structural ensembles representing the dynamics of proteins.

Free energy calculations on the two drug binding sites in the M2 proton channel

Two alternative binding sites of adamantane-type drugs in the influenza A M2 channel have been suggested, one with the drug binding inside the channel pore and the other with four drug molecule S-binding to the C-terminal surface of the transmembrane domain. Recent computational and experimental studies have suggested that the pore binding site is more energetically favorable but the external surface binding site may also exist. Nonetheless, which drug binding site leads to channel inhibition in vivo and how drug-resistant mutations affect these sites are not completely understood. We applied molecular dynamics simulations and potential of mean force calculations to examine the structures and the free energies associated with these putative drug binding sites in an M2-lipid bilayer system. We found that, at biological pH (~7.4), the pore binding site is more thermodynamically favorable than the surface binding site by ~7 kcal/mol and, hence, would lead to more stable drug binding and channel inhibition. This result is in excellent agreement with several recent studies. More importantly, a novel finding of ours is that binding to the channel pore requires overcoming a much higher energy barrier of ~10 kcal/mol than binding to the C-terminal channel surface, indicating that the latter site is more kinetically favorable. Our study is the first computational work that provides both kinetic and thermodynamic energy information on these drug binding sites. Our results provide a theoretical framework to interpret and reconcile existing and often conflicting results regarding these two binding sites, thus helping to expand our understanding of M2-drug binding, and may help guide the design and screening of novel drugs to combat the virus.

Protonation, tautomerization, and rotameric structure of histidine: a comprehensive study by magic-angle-spinning solid-state NMR

Histidine structure and chemistry lie at the heart of many enzyme active sites, ion channels, and metalloproteins. While solid-state NMR spectroscopy has been used to study histidine chemical shifts, the full pH dependence of the complete panel of (15)N, (13)C, and (1)H chemical shifts and the sensitivity of these chemical shifts to tautomeric structure have not been reported. Here we use magic-angle-spinning solid-state NMR spectroscopy to determine the (15)N, (13)C, and (1)H chemical shifts of histidine from pH 4.5 to 11. Two-dimensional homonuclear and heteronuclear correlation spectra indicate that these chemical shifts depend sensitively on the protonation state and tautomeric structure. The chemical shifts of the rare π tautomer were observed for the first time, at the most basic pH used. Intra- and intermolecular hydrogen bonding between the imidazole nitrogens and the histidine backbone or water was detected, and N-H bond length measurements indicated the strength of the hydrogen bond. We also demonstrate the accurate measurement of the histidine side-chain torsion angles χ(1) and χ(2) through backbone-side chain (13)C-(15)N distances; the resulting torsion angles were within 4° of the crystal structure values. These results provide a comprehensive set of benchmark values for NMR parameters of histidine over a wide pH range and should facilitate the study of functionally important histidines in proteins.

Detection of a transient intermediate in a rapid protein folding process by solid-state nuclear magnetic resonance

We describe the use of solid-state NMR spectroscopy to characterize a partially folded state of the 35-residue helical protein HP35 created by rapid freeze-quenching from a thermally unfolded state on the 10-20 micros time scale. Two-dimensional solid-state (13)C NMR spectra of (13)C-labeled HP35 in frozen glycerol/water solution exhibit two sets of signals, one corresponding to strongly unfolded protein molecules and the other to an ensemble of molecules having native helical secondary structure but incomplete tertiary structure. The NMR data indicate that secondary structure forms within the freeze-quenching time scale but that full folding involves a slower phase of structural annealing. The approximately 5 micros folding time observed in earlier studies of HP35 by time-resolved optical techniques may not represent the time scale for full folding.

Structure of the amantadine binding site of influenza M2 proton channels in lipid bilayers

The M2 protein of influenza A virus is a membrane-spanning tetrameric proton channel targeted by the antiviral drugs amantadine and rimantadine. Resistance to these drugs has compromised their effectiveness against many influenza strains, including pandemic H1N1. A recent crystal structure of M2(22-46) showed electron densities attributed to a single amantadine in the amino-terminal half of the pore, indicating a physical occlusion mechanism for inhibition. However, a solution NMR structure of M2(18-60) showed four rimantadines bound to the carboxy-terminal lipid-facing surface of the helices, suggesting an allosteric mechanism. Here we show by solid-state NMR spectroscopy that two amantadine-binding sites exist in M2 in phospholipid bilayers. The high-affinity site, occupied by a single amantadine, is located in the N-terminal channel lumen, surrounded by residues mutated in amantadine-resistant viruses. Quantification of the protein-amantadine distances resulted in a 0.3 A-resolution structure of the high-affinity binding site. The second, low-affinity, site was observed on the C-terminal protein surface, but only when the drug reaches high concentrations in the bilayer. The orientation and dynamics of the drug are distinct in the two sites, as shown by (2)H NMR. These results indicate that amantadine physically occludes the M2 channel, thus paving the way for developing new antiviral drugs against influenza viruses. The study demonstrates the ability of solid-state NMR to elucidate small-molecule interactions with membrane proteins and determine high-resolution structures of their complexes.

Structure and function of the influenza A M2 proton channel

The M2 protein of influenza A viruses forms a tetrameric pH-activated proton-selective channel that is targeted by the amantadine class of antiviral drugs. Its ion channel function has been extensively studied by electrophysiology and mutagenesis; however, the molecular mechanism of proton transport is still elusive, and the mechanism of inhibition by amantadine is controversial. We review the functional data on proton channel activity, molecular dynamics simulations of the proton conduction mechanism, and high-resolution structural and dynamical information of this membrane protein in lipid bilayers and lipid-mimetic detergents. These studies indicate that elucidation of the structural basis of M2 channel activity and inhibition requires thorough examination of the complex dynamics and conformational plasticity of the protein in different lipid bilayers and lipid-mimetic environments.

Mapping the dynamics of ligand reorganization via 13CH3 and 13CH2 relaxation dispersion at natural abundance

Flexible ligands pose challenges to standard structure-activity studies since they frequently reorganize their conformations upon protein binding and catalysis. Here, we demonstrate the utility of side chain (13)C relaxation dispersion measurements to identify and quantify the conformational dynamics that drive this reorganization. The dispersion measurements probe methylene (13)CH(2) and methyl (13)CH(3) groups; the latter are highly prevalent side chain moieties in known drugs. Combining these side chain studies with existing backbone dispersion studies enables a comprehensive investigation of mus-ms conformational dynamics related to binding and catalysis. We perform these measurements at natural (13)C abundance, in congruence with common pharmaceutical research settings. We illustrate these methods through a study of the interaction of a phosphopeptide ligand with the peptidyl-prolyl isomerase, Pin1. The results illuminate the side-chain moieties that undergo conformational readjustments upon complex formation. In particular, we find evidence that multiple exchange processes influence the side chain dispersion profiles. Collectively, our studies illustrate how side-chain relaxation dispersion can shed light on ligand conformational transitions required for activity, and thereby suggest strategies for its optimization.

Structural basis for the function and inhibition of an influenza virus proton channel

The M2 protein from influenza A virus is a pH-activated proton channel that mediates acidification of the interior of viral particles entrapped in endosomes. M2 is the target of the anti-influenza drugs amantadine and rimantadine; recently, resistance to these drugs in humans, birds and pigs has reached more than 90% (ref. 1). Here we describe the crystal structure of the transmembrane-spanning region of the homotetrameric protein in the presence and absence of the channel-blocking drug amantadine. pH-dependent structural changes occur near a set of conserved His and Trp residues that are involved in proton gating. The drug-binding site is lined by residues that are mutated in amantadine-resistant viruses. Binding of amantadine physically occludes the pore, and might also perturb the pK(a) of the critical His residue. The structure provides a starting point for solving the problem of resistance to M2-channel blockers.