A 13C labeling strategy reveals a range of aromatic side chain motion in calmodulin

Improved spectral resolution of [13C,1H]-HSQC spectra of aromatic amino acid residues in proteins produced by cell-free synthesis from inexpensive 13C-labelled precursors

Cell-free protein synthesis using eCells allows production of amino acids from inexpensive 13C-labelled precursors. We show that the metabolic pathway converting pyruvate, glucose and erythrose into aromatic amino acids is maintained in eCells. Judicious choice of 13C-labelled starting material leads to proteins, where the sidechains of aromatic amino acids display [13C,1H]-HSQC cross-peaks free of one-bond 13C-13C couplings. Selective 13C-labelling of tyrosine and phenylalanine residues is achieved simply by using different compositions of the reaction buffers.

Cell-free synthesis of proteins with selectively 13C-labelled methyl groups from inexpensive precursors

The novel eCell system maintains the activity of the entire repertoire of metabolic Escherichia coli enzymes in cell-free protein synthesis. We show that this can be harnessed to produce proteins with selectively 13 C-labelled amino acids from inexpensive 13 C-labelled precursors. The system is demonstrated with selective 13 C labelling of methyl groups in the proteins ubiquitin and peptidyl-prolyl cis-trans isomerase B. Starting from 3-13 C-pyruvate, 13 C-HSQC cross-peaks are obtained devoid of one-bond 13 C-13 C scalar couplings. Starting from 2-13 C-methyl-acetolactate, single methyl groups of valine and leucine are labelled. Labelling efficiencies are 70 % or higher, and the method allows us to produce perdeuterated proteins with protonated methyl groups in a residue-selective manner. The system uses the isotope-labelled precursors sparingly and is readily scalable.

NMR Studies of Aromatic Ring Flips to Probe Conformational Fluctuations in Proteins

Aromatic residues form a significant part of the protein core, where they make tight interactions with multiple surrounding side chains. Despite the dense packing of internal side chains, the aromatic rings of phenylalanine and tyrosine residues undergo 180° rotations, or flips, which are mediated by transient and large-scale "breathing" motions that generate sufficient void volume around the aromatic ring. Forty years after the seminal work by Wagner and Wüthrich, NMR studies of aromatic ring flips are now undergoing a renaissance as a powerful means of probing fundamental dynamic properties of proteins. Recent developments of improved NMR methods and isotope labeling schemes have enabled a number of advances in addressing the mechanisms and energetics of aromatic ring flips. The nature of the transition states associated with ring flips can be described by thermodynamic activation parameters, including the activation enthalpy, activation entropy, activation volume, and also the isothermal volume compressibility of activation. Consequently, it is of great interest to study how ring flip rate constants and activation parameters might vary with protein structure and external conditions like temperature and pressure. The field is beginning to gather such data for aromatic residues in a variety of environments, ranging from surface exposed to buried. In the future, the combination of solution and solid-state NMR spectroscopy together with molecular dynamics simulations and other computational approaches is likely to provide detailed information about the coupled dynamics of aromatic rings and neighboring residues. In this Perspective, we highlight recent developments and provide an outlook toward the future.

Aromatic ring flips in differently packed ubiquitin protein crystals from MAS NMR and MD

Probing the dynamics of aromatic side chains provides important insights into the behavior of a protein because flips of aromatic rings in a protein's hydrophobic core report on breathing motion involving a large part of the protein. Inherently invisible to crystallography, aromatic motions have been primarily studied by solution NMR. The question how packing of proteins in crystals affects ring flips has, thus, remained largely unexplored. Here we apply magic-angle spinning NMR, advanced phenylalanine 1H-13C/2H isotope labeling and MD simulation to a protein in three different crystal packing environments to shed light onto possible impact of packing on ring flips. The flips of the two Phe residues in ubiquitin, both surface exposed, appear remarkably conserved in the different crystal forms, even though the intermolecular packing is quite different: Phe4 flips on a ca. 10-20 ns time scale, and Phe45 are broadened in all crystals, presumably due to µs motion. Our findings suggest that intramolecular influences are more important for ring flips than intermolecular (packing) effects.

Characterizing Fast Conformational Exchange of Aromatic Rings Using Residual Dipolar Couplings: Distinguishing Jumplike Flips from Other Exchange Mechanisms

Aromatic ring flips are a hallmark of protein dynamics. They are experimentally studied by NMR spectroscopy, where recent advances have led to improved characterization across a wide range of time scales. Results on different proteins have been interpreted as continuous diffusive ring rotations or jumplike flips, leading to diverging views of the protein interior as being fluidlike or solidlike, respectively. It is challenging to distinguish between these mechanisms and other types of conformational exchange because chemical-shift-mediated line broadening provides only conclusive evidence for ring flips only if the system can be moved from the slow- to intermediate/fast-exchange regime. Moreover, whenever the chemical shift difference between the two symmetry-related sites is close to zero, it is not generally possible to determine the exchange time scale. Here we resolve these issues by measuring residual dipolar coupling (RDC)-mediated exchange contributions using NMR relaxation dispersion experiments on proteins dissolved in dilute liquid crystalline media. Excellent agreement is found between the experimental difference in RDC between the two symmetry-related sites and the value calculated from high-resolution X-ray structures, demonstrating that dynamics measured for F52 in the B1 domain of protein G reports on distinct, jumplike flips rather than other types of conformational exchange.

Measurement of 1Hα transverse relaxation rates in proteins: application to solvent PREs

It has recently been demonstrated that accurate near surface electrostatic potentials can be calculated for proteins from solvent paramagnetic relaxation enhancements (PREs) of amide protons measured using spin labels of similar structures but different charges (Yu et al. in Proc Natl Acad Sci 118(25):e2104020118, 2021). Here we develop methodology for extending such measurements to intrinsically disordered proteins at neutral pH where amide spectra are of very poor quality. Under these conditions it is shown that accurate PRE values can be measured using the haCONHA experiment that has been modified for recording ¹H^α transverse relaxation rates. The optimal pulse scheme includes a spin-lock relaxation element for suppression of homonuclear scalar coupled evolution for all ¹H^α protons, except those derived from Ser and Thr residues, and minimizes the radiation damping field from water magnetization that would otherwise increase measured relaxation rates. The robustness of the experiment is verified by developing a second approach using a band selective adiabatic decoupling scheme for suppression of scalar coupling modulations during ¹H^α relaxation and showing that the measured PRE values from the two methods are in excellent agreement. The near surface electrostatic potential of a 103-residue construct comprising the C-terminal intrinsically disordered region of the RNA-binding protein CAPRIN1 is obtained at pH 5.5 using both ¹H^N and ¹H^α-based relaxation rates, and at pH 7.4 where only ¹H^α rates can be quantified, with very good agreement between potentials obtained under all experimental conditions.

1H R1ρ relaxation dispersion experiments in aromatic side chains

Aromatic side chains are attractive probes of protein dynamic, since they are often key residues in enzyme active sites and protein binding sites. Dynamic processes on microsecond to millisecond timescales can be studied by relaxation dispersion experiments that attenuate conformational exchange contributions to the transverse relaxation rate by varying the refocusing frequency of applied radio-frequency fields implemented as either CPMG pulse trains or continuous spin-lock periods. Here we present an aromatic ¹H R_1ρ relaxation dispersion experiment enabling studies of two to three times faster exchange processes than achievable by existing experiments for aromatic side chains. We show that site-specific isotope labeling schemes generating isolated ¹H-¹³C spin pairs with vicinal ²H-¹²C moieties are necessary to avoid anomalous relaxation dispersion profiles caused by Hartmann-Hahn matching due to the ³J_HH couplings and limited chemical shift differences among ¹H spins in phenylalanine, tyrosine and the six-ring moiety of tryptophan. This labeling pattern is sufficient in that remote protons do not cause additional complications. We validated the approach by measuring ring-flip kinetics in the small protein GB1. The determined rate constants, k_flip, agree well with previous results from ¹³C R_1ρ relaxation dispersion experiments, and yield ¹H chemical shift differences between the two sides of the ring in good agreement with values measured under slow-exchange conditions. The aromatic¹H R_1ρ relaxation dispersion experiment in combination with the site-selective ¹H-¹³C/²H-¹²C labeling scheme enable measurement of exchange rates up to k_ex = 2k_flip = 80,000 s^-1, and serve as a useful complement to previously developed ¹³C-based methods.

Transition-State Compressibility and Activation Volume of Transient Protein Conformational Fluctuations

Proteins are dynamic entities that intermittently depart from their ground-state structures and undergo conformational transitions as a critical part of their functions. Central to understanding such transitions are the structural rearrangements along the connecting pathway, where the transition state plays a special role. Using NMR relaxation at variable temperature and pressure to measure aromatic ring flips inside a protein core, we obtain information on the structure and thermodynamics of the transition state. We show that the isothermal compressibility coefficient of the transition state is similar to that of short-chain hydrocarbon liquids, implying extensive local unfolding of the protein. Our results further indicate that the required local volume expansions of the protein can occur not only with a net positive activation volume of the protein, as expected from previous studies, but also with zero activation volume by compaction of remote void volume, when averaged over the ensemble of states.

Modulating Enzyme Function via Dynamic Allostery within Biliverdin Reductase B

The biliverdin reductase B (BLVRB) class of enzymes catalyze the NADPH-dependent reduction of multiple flavin substrates and are emerging as critical players in cellular redox regulation. However, the role of dynamics and allostery have not been addressed, prompting studies here that have revealed a position 15 Å away from the active site within human BLVRB (T164) that is inherently dynamic and can be mutated to control global micro-millisecond motions and function. By comparing the inherent dynamics through nuclear magnetic resonance (NMR) relaxation approaches of evolutionarily distinct BLVRB homologues and by applying our previously developed Relaxation And Single Site Multiple Mutations (RASSMM) approach that monitors both the functional and dynamic effects of multiple mutations to the single T164 site, we have discovered that the most dramatic mutagenic effects coincide with evolutionary changes and these modulate coenzyme binding. Thus, evolutionarily changing sites distal to the active site serve as dynamic "dials" to globally modulate motions and function. Despite the distal dynamic and functional coupling modulated by this site, micro-millisecond motions span an order of magnitude in their apparent kinetic rates of motions. Thus, global dynamics within BLVRB are a collection of partially coupled motions tied to catalytic function.

Protein conformational entropy is not slaved to water

Conformational entropy can be an important element of the thermodynamics of protein functions such as the binding of ligands. The observed role for conformational entropy in modulating molecular recognition by proteins is in opposition to an often-invoked theory for the interaction of protein molecules with solvent water. The "solvent slaving" model predicts that protein motion is strongly coupled to various aspects of water such as bulk solvent viscosity and local hydration shell dynamics. Changes in conformational entropy are manifested in alterations of fast internal side chain motion that is detectable by NMR relaxation. We show here that the fast-internal side chain dynamics of several proteins are unaffected by changes to the hydration layer and bulk water. These observations indicate that the participation of conformational entropy in protein function is not dictated by the interaction of protein molecules and solvent water under the range of conditions normally encountered.

Site-selective 1H/2H labeling enables artifact-free 1H CPMG relaxation dispersion experiments in aromatic side chains

Aromatic side chains are often key residues in enzyme active sites and protein binding sites, making them attractive probes of protein dynamics on the millisecond timescale. Such dynamic processes can be studied by aromatic ¹³C or ¹H CPMG relaxation dispersion experiments. Aromatic ¹H CPMG relaxation dispersion experiments in phenylalanine, tyrosine and the six-ring moiety of tryptophan, however, are affected by ³J ¹H-¹H couplings which are causing anomalous relaxation dispersion profiles. Here we show that this problem can be addressed by site-selective ¹H/²H labeling of the aromatic side chains and that artifact-free relaxation dispersion profiles can be acquired. The method has been further validated by measuring folding-unfolding kinetics of the small protein GB1. The determined rate constants and populations agree well with previous results from ¹³C CPMG relaxation dispersion experiments. Furthermore, the CPMG-derived chemical shift differences between the folded and unfolded states are in excellent agreement with those obtained directly from the spectra. In summary, site-selective ¹H/²H labeling enables artifact-free aromatic ¹H CPMG relaxation dispersion experiments in phenylalanine and the six-ring moiety of tryptophan, thereby extending the available methods for studying millisecond dynamics in aromatic protein side chains.

Emerging solution NMR methods to illuminate the structural and dynamic properties of proteins

The first recognition of protein breathing was more than 50 years ago. Today, we are able to detect the multitude of interaction modes, structural polymorphisms, and binding-induced changes in protein structure that direct function. Solution-state NMR spectroscopy has proved to be a powerful technique, not only to obtain high-resolution structures of proteins, but also to provide unique insights into the functional dynamics of proteins. Here, we summarize recent technical landmarks in solution NMR that have enabled characterization of key biological macromolecular systems. These methods have been fundamental to atomic resolution structure determination and quantitative analysis of dynamics over a wide range of time scales by NMR. The ability of NMR to detect lowly populated protein conformations and transiently formed complexes plays a critical role in its ability to elucidate functionally important structural features of proteins and their dynamics.

Structure and Dynamics of Stacking Interactions in an Antibody Binding Site

For years, antibodies (Abs) have been used as a paradigm for understanding how protein structure contributes to molecular recognition. However, with the ability to evolve Abs that recognize specific chromophores, they also have great potential as models for how protein dynamics contribute to molecular recognition. We previously raised murine Abs to different chromophores and, with the use of three-pulse photon echo peak shift spectroscopy, demonstrated that the immune system is capable of producing Abs with widely varying flexibility. We now report the characterization of the complexes formed between two Abs, 5D11 and 10A6, and the chromophoric ligand that they were evolved to recognize, 8-methoxypyrene-1,3,6-trisulfonic acid (MPTS). The sequences of the Ab genes indicate that they evolved from a common precursor. We also used a variety of spectroscopic methods to probe the photophysics and dynamics of the Ab-MPTS complexes and found that they are similar to each other but distinct from previously characterized anti-MPTS Abs. Structural studies revealed that this difference likely results from a unique mode of binding in which MPTS is sandwiched between the side chain of Phe^H98, which interacts with the chromophore via T-stacking, and the side chain of Trp^L91, which interacts with the chromophore via parallel stacking. The T-stacking interaction appears to mediate relaxation on the picosecond time scale, while the parallel stacking appears to mediate relaxation on an ultrafast, femtosecond time scale, which dominates the response. The anti-MPTS Abs thus not only demonstrate the simultaneous use of the two limiting modes of stacking for molecular recognition, but also provide a unique opportunity to characterize how dynamics might contribute to molecular recognition. Both types of stacking are common in proteins and protein complexes where they may similarly contribute to dynamics and molecular recognition.

Aromatic Ring Dynamics, Thermal Activation, and Transient Conformations of a 468 kDa Enzyme by Specific 1H-13C Labeling and Fast Magic-Angle Spinning NMR

Aromatic residues are located at structurally important sites of many proteins. Probing their interactions and dynamics can provide important functional insight but is challenging in large proteins. Here, we introduce approaches to characterize the dynamics of phenylalanine residues using ¹H-detected fast magic-angle spinning (MAS) NMR combined with a tailored isotope-labeling scheme. Our approach yields isolated two-spin systems that are ideally suited for artifact-free dynamics measurements, and allows probing motions effectively without molecular weight limitations. The application to the TET2 enzyme assembly of ∼0.5 MDa size, the currently largest protein assigned by MAS NMR, provides insights into motions occurring on a wide range of time scales (picoseconds to milliseconds). We quantitatively probe ring-flip motions and show the temperature dependence by MAS NMR measurements down to 100 K. Interestingly, favorable line widths are observed down to 100 K, with potential implications for DNP NMR. Furthermore, we report the first ¹³C R_1ρ MAS NMR relaxation-dispersion measurements and detect structural excursions occurring on a microsecond time scale in the entry pore to the catalytic chamber and at a trimer interface that was proposed as the exit pore. We show that the labeling scheme with deuteration at ca. 50 kHz MAS provides superior resolution compared to 100 kHz MAS experiments with protonated, uniformly ¹³C-labeled samples.

Characterizing Protein Hydration Dynamics Using Solution NMR Spectroscopy

Protein hydration is a critical aspect of protein stability, folding, and function and yet remains difficult to characterize experimentally. Solution NMR offers a route to a site-resolved view of the dynamics of protein-water interactions through the nuclear Overhauser effects between hydration water and the protein in the laboratory (NOE) and rotating (ROE) frames of reference. However, several artifacts and limitations including contaminating contributions from bulk water potentially plague this general approach and the corruption of measured NOEs and ROEs by hydrogen exchange-relayed magnetization. Fortunately, encapsulation of single protein molecules within the water core of a reverse micelle overcomes these limitations. The main advantages are the suppression hydrogen exchange and elimination of bulk water. Here we detail guidelines for the preparation solutions of encapsulated proteins that are suitable for characterization by NOE and ROE spectroscopy. Emphasis is placed on understanding the contribution of detected NOE intensity arising from magnetization relayed by hydrogen exchange. Various aspects of fitting obtained NOE, selectively decoupled NOE, and ROE time courses are illustrated.

Conformational exchange of aromatic side chains by 1H CPMG relaxation dispersion

Aromatic side chains are attractive probes of protein dynamics on the millisecond time scale, because they are often key residues in enzyme active sites and protein binding sites. Further they allow to study specific processes, like histidine tautomerization and ring flips. Till now such processes have been studied by aromatic ¹³C CPMG relaxation dispersion experiments. Here we investigate the possibility of aromatic ¹H CPMG relaxation dispersion experiments as a complementary method. Artifact-free dispersions are possible on uniformly ¹H and ¹³C labeled samples for histidine δ2 and ε1, as well as for tryptophan δ1. The method has been validated by measuring fast folding-unfolding kinetics of the small protein CspB under native conditions. The determined rate constants and populations agree well with previous results from ¹³C CPMG relaxation dispersion experiments. The CPMG-derived chemical shift differences between the folded and unfolded states are in good agreement with those obtained directly from the spectra. In contrast, the ¹H relaxation dispersion profiles in phenylalanine, tyrosine and the six-ring moiety of tryptophan, display anomalous behavior caused by ³J ¹H-¹H couplings and, if present, strong ¹³C-¹³C couplings. Therefore they require site-selective ¹H/²H and, in case of strong couplings, ¹³C/¹²C labeling. In summary, aromatic ¹H CPMG relaxation dispersion experiments work on certain positions (His δ2, His ε1 and Trp δ1) in uniformly labeled samples, while other positions require site-selective isotope labeling.

Optimal Isotope Labeling of Aromatic Amino Acid Side Chains for NMR Studies of Protein Dynamics

Aromatic side chains in proteins are often directly evolved in stabilizing the hydrophobic core, protein binding, or enzymatic activity. They are also responsible for specific local dynamic processes, such as histidine tautomerization or ring flips. Despite their importance, they are often not targeted directly by NMR spectroscopy, because of spectroscopic complications and challenges. This chapter addresses state-of-the-art site-selective ¹³C-labeling methods for aromatic side chains, and describes how they solve several of the spectroscopic issues. A special emphasis is put on thereby enabled protein dynamics experiments of aromatic side chains.

Late metabolic precursors for selective aromatic residue labeling

In recent years, we developed a toolbox of heavy isotope containing compounds, which serve as metabolic amino acid precursors in the E. coli-based overexpression of aromatic residue labeled proteins. Our labeling techniques show excellent results both in terms of selectivity and isotope incorporation levels. They are additionally distinguished by low sample production costs and meet the economic demands to further implement protein NMR spectroscopy as a routinely used method in drug development processes. Different isotopologues allow for the assembly of optimized protein samples, which fulfill the requirements of various NMR experiments to elucidate protein structures, analyze conformational dynamics, or probe interaction surfaces. In the present article, we want to summarize the precursors we developed so far and give examples of their special value in the probing of protein-ligand interaction.

Cyanylated Cysteine Reports Site-Specific Changes at Protein-Protein-Binding Interfaces Without Perturbation

To investigate the cyanylated cysteine vibrational probe group’s ability to report on binding-induced changes along a protein–protein interface, the probe group was incorporated at several sites in a peptide of the calmodulin (CaM)-binding domain of skeletal muscle myosin light chain kinase. Isothermal titration calorimetry was used to determine the binding thermodynamics between calmodulin and each peptide. For all probe positions, the binding affinity was nearly identical to that of the unlabeled peptide. The CN stretching infrared band was collected for each peptide free in solution and bound to calmodulin. Binding-induced shifts in the IR spectral frequencies were correlated with estimated solvent accessibility based on molecular dynamics simulations. This work generally suggests (1) that site-specific incorporation of this vibrational probe group does not cause major perturbations to its local structural environment and (2) that this small probe group might be used quite broadly to map dynamic protein-binding interfaces. However, site-specific perturbations due to artificial labeling groups can be somewhat unpredictable and should be evaluated on a site-by-site basis through complementary measurements. A fully quantitative, simulation-based interpretation of the rich probe IR spectra is still needed but appears to be possible given recent advances in simulation techniques.

Measuring Entropy in Molecular Recognition by Proteins

Molecular recognition by proteins is fundamental to the molecular basis of biology. Dissection of the thermodynamic landscape governing protein-ligand interactions has proven difficult because determination of various entropic contributions is quite challenging. Nuclear magnetic resonance relaxation measurements, theory, and simulations suggest that conformational entropy can be accessed through a dynamical proxy. Here, we review the relationship between measures of fast side-chain motion and the underlying conformational entropy. The dynamical proxy reveals that the contribution of conformational entropy can range from highly favorable to highly unfavorable and demonstrates the potential of this key thermodynamic variable to modulate protein-ligand interactions. The dynamical so-called entropy meter also refines the role of solvent entropy and directly determines the loss in rotational-translational entropy that occurs upon formation of high-affinity complexes. The ability to quantify the roles of entropy through an entropy meter based on measurable dynamical properties promises to highlight its role in protein function.

Anthranilic acid, the new player in the ensemble of aromatic residue labeling precursor compounds

The application of metabolic precursors for selective stable isotope labeling of aromatic residues in cell-based protein overexpression has already resulted in numerous NMR probes to study the structural and dynamic characteristics of proteins. With anthranilic acid, we present the structurally simplest precursor for exclusive tryptophan side chain labeling. A synthetic route to ¹³C, ²H isotopologues allows the installation of isolated ¹³C-¹H spin systems in the indole ring of tryptophan, representing a versatile tool to investigate side chain motion using relaxation-based experiments without the loss of magnetization due to strong ¹J_CC and weaker ²J_CH scalar couplings, as well as dipolar interactions with remote hydrogens. In this article, we want to introduce this novel precursor in the context of hitherto existing techniques of in vivo aromatic residue labeling.

Site-selective 13C labeling of histidine and tryptophan using ribose

Experimental studies on protein dynamics at atomic resolution by NMR-spectroscopy in solution require isolated ¹H-X spin pairs. This is the default scenario in standard ¹H-¹⁵N backbone experiments. Side chain dynamic experiments, which allow to study specific local processes like proton-transfer, or tautomerization, require isolated ¹H-¹³C sites which must be produced by site-selective ¹³C labeling. In the most general way this is achieved by using site-selectively ¹³C-enriched glucose as the carbon source in bacterial expression systems. Here we systematically investigate the use of site-selectively ¹³C-enriched ribose as a suitable precursor for ¹³C labeled histidines and tryptophans. The ¹³C incorporation in nearly all sites of all 20 amino acids was quantified and compared to glucose based labeling. In general the ribose approach results in more selective labeling. 1-¹³C ribose exclusively labels His δ2 and Trp δ1 in aromatic side chains and helps to resolve possible overlap problems. The incorporation yield is however only 37% in total and 72% compared to yields of 2-¹³C glucose. A combined approach of 1-¹³C ribose and 2-¹³C glucose maximizes ¹³C incorporation to 75% in total and 150% compared to 2-¹³C glucose only. Further histidine positions β, α and CO become significantly labeled at around 50% in total by 3-, 4- or 5-¹³C ribose. Interestingly backbone CO of Gly, Ala, Cys, Ser, Val, Phe and Tyr are labeled at 40-50% in total with 3-¹³C ribose, compared to 5% and below for 1-¹³C and 2-¹³C glucose. Using ribose instead of glucose as a source for site-selective ¹³C labeling enables a very selective labeling of certain positions and thereby expanding the toolbox for customized isotope labeling of amino-acids.

NMR Studies of Large Proteins

Recent breakthroughs in isotope-labeling and pulse sequence techniques have enabled the NMR characterization of large protein systems with molecular masses of hundreds of kilodaltons. NMR studies of a great variety of large proteins have provided unique insights into the binding, dynamic, and allosteric mechanisms. Here we present a brief summary of these developments by highlighting few cases that exemplify the uniqueness of NMR in providing atomic resolution information into key dynamic processes and structures of protein complexes with high degree of flexibility.

Highly Selective Stable Isotope Labeling of Histidine Residues by Using a Novel Precursor in E. coli-Based Overexpression Systems

The importance of NMR spectroscopy in unraveling the structural and dynamic properties of proteins is ever-expanding owing to progress in experimental techniques, hardware development, and novel labeling approaches. Multiple sophisticated methods of aliphatic residue labeling can be found in the literature, whereas the selective incorporation of NMR active isotopes into other amino acids still holds the potential for improvement. In order to close this methodological gap, we present a novel metabolic precursor for cell-based protein overexpression to assemble ¹³ C/² H isotope patterns in the peptide backbone, as well as in side chain positions of a mechanistically distinguished histidine residue.

Entropy in molecular recognition by proteins

Molecular recognition by proteins is fundamental to molecular biology. Dissection of the thermodynamic energy terms governing protein-ligand interactions has proven difficult, with determination of entropic contributions being particularly elusive. NMR relaxation measurements have suggested that changes in protein conformational entropy can be quantitatively obtained through a dynamical proxy, but the generality of this relationship has not been shown. Twenty-eight protein-ligand complexes are used to show a quantitative relationship between measures of fast side-chain motion and the underlying conformational entropy. We find that the contribution of conformational entropy can range from favorable to unfavorable, which demonstrates the potential of this thermodynamic variable to modulate protein-ligand interactions. For about one-quarter of these complexes, the absence of conformational entropy would render the resulting affinity biologically meaningless. The dynamical proxy for conformational entropy or "entropy meter" also allows for refinement of the contributions of solvent entropy and the loss in rotational-translational entropy accompanying formation of high-affinity complexes. Furthermore, structure-based application of the approach can also provide insight into long-lived specific water-protein interactions that escape the generic treatments of solvent entropy based simply on changes in accessible surface area. These results provide a comprehensive and unified view of the general role of entropy in high-affinity molecular recognition by proteins.

Probing conformational dynamics in biomolecules via chemical exchange saturation transfer: a primer

Although Chemical Exchange Saturation Transfer (CEST) type NMR experiments have been used to study chemical exchange processes in molecules since the early 1960s, there has been renewed interest in the past several years in using this approach to study biomolecular conformational dynamics. The methodology is particularly powerful for the study of sparsely populated, transiently formed conformers that are recalcitrant to investigation using traditional biophysical tools, and it is complementary to relaxation dispersion and magnetization transfer experiments that have traditionally been used to study chemical exchange processes. Here we discuss the concepts behind the CEST experiment, focusing on practical aspects as well, we review some of the pulse sequences that have been developed to characterize protein and RNA conformational dynamics, and we discuss a number of examples where the CEST methodology has provided important insights into the role of dynamics in biomolecular function.

Site-selective 13C labeling of proteins using erythrose

NMR-spectroscopy enables unique experimental studies on protein dynamics at atomic resolution. In order to obtain a full atom view on protein dynamics, and to study specific local processes like ring-flips, proton-transfer, or tautomerization, one has to perform studies on amino-acid side chains. A key requirement for these studies is site-selective labeling with ¹³C and/or ¹H, which is achieved in the most general way by using site-selectively ¹³C-enriched glucose (1- and 2-¹³C) as the carbon source in bacterial expression systems. Using this strategy, multiple sites in side chains, including aromatics, become site-selectively labeled and suitable for relaxation studies. Here we systematically investigate the use of site-selectively ¹³C-enriched erythrose (1-, 2-, 3- and 4-¹³C) as a suitable precursor for ¹³C labeled aromatic side chains. We quantify ¹³C incorporation in nearly all sites in all 20 amino acids and compare the results to glucose based labeling. In general the erythrose approach results in more selective labeling. While there is only a minor gain for phenylalanine and tyrosine side-chains, the ¹³C incorporation level for tryptophan is at least doubled. Additionally, the Phe ζ and Trp η2 positions become labeled. In the aliphatic side chains, labeling using erythrose yields isolated ¹³C labels for certain positions, like Ile β and His β, making these sites suitable for dynamics studies. Using erythrose instead of glucose as a source for site-selective ¹³C labeling enables unique or superior labeling for certain positions and is thereby expanding the toolbox for customized isotope labeling of amino-acid side-chains.

Bacterial production of site specific 13C labeled phenylalanine and methodology for high level incorporation into bacterially expressed recombinant proteins

Nuclear magnetic resonance spectroscopy studies of ever larger systems have benefited from many different forms of isotope labeling, in particular, site specific isotopic labeling. Site specific ¹³C labeling of methyl groups has become an established means of probing systems not amenable to traditional methodology. However useful, methyl reporter sites can be limited in number and/or location. Therefore, new complementary site specific isotope labeling strategies are valuable. Aromatic amino acids make excellent probes since they are often found at important interaction interfaces and play significant structural roles. Aromatic side chains have many of the same advantages as methyl containing amino acids including distinct ¹³C chemical shifts and multiple magnetically equivalent ¹H positions. Herein we report economical bacterial production and one-step purification of phenylalanine with ¹³C incorporation at the Cα, Cγ and Cε positions, resulting in two isolated ¹H-¹³C spin systems. We also present methodology to maximize incorporation of phenylalanine into recombinantly overexpressed proteins in bacteria and demonstrate compatibility with ILV-methyl labeling. Inexpensive, site specific isotope labeled phenylalanine adds another dimension to biomolecular NMR, opening new avenues of study.

Dynamics of Hydrophobic Core Phenylalanine Residues Probed by Solid-State Deuteron NMR

We conducted a detailed investigation of the dynamics of two phenylalanine side chains in the hydrophobic core of the villin headpiece subdomain protein (HP36) in the hydrated powder state over the 298-80 K temperature range. Our main tools were static deuteron NMR measurements of longitudinal relaxation and line shapes supplemented with computational modeling. The temperature dependence of the relaxation times reveals the presence of two main mechanisms that can be attributed to the ring-flips, dominating at high temperatures, and small-angle fluctuations, dominating at low temperatures. The relaxation is nonexponential at all temperatures with the extent of nonexponentiality increasing from higher to lower temperatures. This behavior suggests a distribution of conformers with unique values of activation energies. The central values of the activation energies for the ring-flipping motions are among the smallest reported for aromatic residues in peptides and proteins and point to a very mobile hydrophobic core. The analysis of the widths of the distributions, in combination with the earlier results on the dynamics of flanking methyl groups (Vugmeyster et al. J. Phys. Chem. B 2013, 117, 6129-6137), suggests that the hydrophobic core undergoes slow concerted fluctuations. There is a pronounced effect of dehydration on the ring-flipping motions, which shifts the distribution toward more rigid conformers. The crossover temperature between the regions of dominance of the small-angle fluctuations and ring-flips shifts from 195 K in the hydrated protein to 278 K in the dry one. This result points to the role of solvent in softening the core and highlights aromatic residues as markers of the protein dynamical transitions.

A Semiautomated Assignment Protocol for Methyl Group Side Chains in Large Proteins

The developments of biosynthetic specific labeling strategies for side-chain methyl groups have allowed structural and dynamic characterization of very large proteins and protein complexes. However, the assignment of the methyl-group resonances remains an Achilles' heel for NMR, as the experiments designed to correlate side chains to the protein backbone become rather insensitive with the increase of the transverse relaxation rates. In this chapter, we outline a semiempirical approach to assign the resonances of methyl-group side chains in large proteins. This method requires a crystal structure or an NMR ensemble of conformers as an input, together with NMR data sets such as nuclear Overhauser effects (NOEs) and paramagnetic relaxation enhancements (PREs), to be implemented in a computational protocol that provides a probabilistic assignment of methyl-group resonances. As an example, we report the protocol used in our laboratory to assign the side chains of the 42-kDa catalytic subunit of the cAMP-dependent protein kinase A. Although we emphasize the labeling of isoleucine, leucine, and valine residues, this method is applicable to other methyl group side chains such as those of alanine, methionine, and threonine, as well as reductively methylated cysteine side chains.

Novel approaches in selective tryptophan isotope labeling by using Escherichia coli overexpression media

NMR-based investigations of large protein complexes require optimized isotopic labeling schemes. We report new methods to introduce stable isotopes into tryptophan residues; these are fine-tuned to the requirements of the particular protein NMR experiment. Selective backbone labeling was performed by using a new α-ketoacid precursor as an additive in cell-based overexpression media. Additionally, we developed synthetic routes to certain isotopologues of indole with (13)C-(1)H spin systems surrounded by (12)C and (2)H. The corresponding proteins, overexpressed in the presence of these precursor compounds, can be effectively analyzed for conformational changes in tryptophan residues in response to external stimuli, such as interaction with other proteins or small molecules.

A sharp thermal transition of fast aromatic-ring dynamics in ubiquitin

Aromatic amino acid side chains have a rich role within proteins and are often central to their structure and function. Suitable isotopic-labelling strategies enable studies of sub-nanosecond aromatic-ring dynamics using solution NMR relaxation methods. Surprisingly, it was found that the three aromatic side chains in human ubiquitin show a sharp thermal dynamical transition at approximately 312 K. Hydrostatic pressure has little effect on the low-temperature behavior, but somewhat decreases the amplitude of motion in the high-temperature regime. Therefore, below the transition temperature, ring motion is largely librational. Above this temperature, a complete ring-rotation process that is fully consistent with a continuous diffusion not requiring the transient creation of a large activated free volume occurs. Molecular dynamics simulations qualitatively corroborate this view and reinforce the notion that the dynamical character of the protein interior has much more liquid-alkane-like properties than previously appreciated.

Dynamic and thermodynamic response of the Ras protein Cdc42Hs upon association with the effector domain of PAK3

Human cell division cycle protein 42 (Cdc42Hs) is a small, Rho-type guanosine triphosphatase involved in multiple cellular processes through its interactions with downstream effectors. The binding domain of one such effector, the actin cytoskeleton-regulating p21-activated kinase 3, is known as PBD46. Nitrogen-15 backbone and carbon-13 methyl NMR relaxation was measured to investigate the dynamical changes in activated GMPPCP·Cdc42Hs upon PBD46 binding. Changes in internal motion of the Cdc42Hs, as revealed by methyl axis order parameters, were observed not only near the Cdc42Hs-PBD46 interface but also in remote sites on the Cdc42Hs molecule. The binding-induced changes in side-chain dynamics propagate along the long axis of Cdc42Hs away from the site of PBD46 binding with sharp distance dependence. Overall, the binding of the PBD46 effector domain on the dynamics of methyl-bearing side chains of Cdc42Hs results in a modest rigidification, which is estimated to correspond to an unfavorable change in conformational entropy of approximately -10kcalmol(-1) at 298K. A cluster of methyl probes closest to the nucleotide-binding pocket of Cdc42Hs becomes more rigid upon binding of PBD46 and is proposed to slow the catalytic hydrolysis of the γ phosphate moiety. An additional cluster of methyl probes surrounding the guanine ring becomes more flexible on binding of PBD46, presumably facilitating nucleotide exchange mediated by a guanosine exchange factor. In addition, the Rho insert helix, which is located at a site remote from the PBD46 binding interface, shows a significant dynamic response to PBD46 binding.

Off-resonance rotating-frame relaxation dispersion experiment for 13C in aromatic side chains using L-optimized TROSY-selection

Protein dynamics on the microsecond-millisecond time scales often play a critical role in biological function. NMR relaxation dispersion experiments are powerful approaches for investigating biologically relevant dynamics with site-specific resolution, as shown by a growing number of publications on enzyme catalysis, protein folding, ligand binding, and allostery. To date, the majority of studies has probed the backbone amides or side-chain methyl groups, while experiments targeting other sites have been used more sparingly. Aromatic side chains are useful probes of protein dynamics, because they are over-represented in protein binding interfaces, have important catalytic roles in enzymes, and form a sizable part of the protein interior. Here we present an off-resonance R 1ρ experiment for measuring microsecond to millisecond conformational exchange of aromatic side chains in selectively (13)C labeled proteins by means of longitudinal- and transverse-relaxation optimization. Using selective excitation and inversion of the narrow component of the (13)C doublet, the experiment achieves significant sensitivity enhancement in terms of both signal intensity and the fractional contribution from exchange to transverse relaxation; additional signal enhancement is achieved by optimizing the longitudinal relaxation recovery of the covalently attached (1)H spins. We validated the L-TROSY-selected R 1ρ experiment by measuring exchange parameters for Y23 in bovine pancreatic trypsin inhibitor at a temperature of 328 K, where the ring flip is in the fast exchange regime with a mean waiting time between flips of 320 μs. The determined chemical shift difference matches perfectly with that measured from the NMR spectrum at lower temperatures, where separate peaks are observed for the two sites. We further show that potentially complicating effects of strong scalar coupling between protons (Weininger et al. in J Phys Chem B 117: 9241-9247, 2013b) can be accounted for using a simple expression, and provide recommendations for data acquisition when the studied system exhibits this behavior. The present method extends the repertoire of relaxation methods tailored for aromatic side chains by enabling studies of faster processes and improved control over artifacts due to strong coupling.

Optimized reverse micelle surfactant system for high-resolution NMR spectroscopy of encapsulated proteins and nucleic acids dissolved in low viscosity fluids

An optimized reverse micelle surfactant system has been developed for solution nuclear magnetic resonance studies of encapsulated proteins and nucleic acids dissolved in low viscosity fluids. Comprising the nonionic 1-decanoyl-rac-glycerol and the zwitterionic lauryldimethylamine-N-oxide (10MAG/LDAO), this mixture is shown to efficiently encapsulate a diverse set of proteins and nucleic acids. Chemical shift analyses of these systems show that high structural fidelity is achieved upon encapsulation. The 10MAG/LDAO surfactant system reduces the molecular reorientation time for encapsulated macromolecules larger than ∼20 kDa leading to improved overall NMR performance. The 10MAG/LDAO system can also be used for solution NMR studies of lipid-modified proteins. New and efficient strategies for optimization of encapsulation conditions are described. 10MAG/LDAO performs well in both the low viscosity pentane and ultralow viscosity liquid ethane and therefore will serve as a general surfactant system for initiating solution NMR studies of proteins and nucleic acids.

Evidence against the "Y-T coupling" mechanism of activation in the response regulator NtrC

The dominant theory on the mechanism of response regulators activation in two-component bacterial signaling systems is the "Y-T coupling" mechanism, wherein the χ1 rotameric state of a highly conserved aromatic residue correlates with the activation of the protein via structural rearrangements coupled to a conserved tyrosine. In this paper, we present evidence that, in the receiver domain of the response regulator nitrogen regulatory protein C (NtrC(R)), the interconversion of this tyrosine (Y101) between its rotameric states is actually faster than the rate of inactive/active conversion and is not correlated to the activation process. Data gathered from NMR relaxation dispersion experiments show that a subset of residues surrounding the conserved tyrosine sense a process that is occurring at a faster rate than the inactive/active conformational transition. We show that this process is related to χ1 rotamer exchange of Y101 and that mutation of this aromatic residue to a leucine eliminated this second faster process without affecting activation. Computational simulations of NtrC(R) in its active conformation further demonstrate that the rotameric state of Y101 is uncorrelated with the global conformational transition during activation. Moreover, the tyrosine does not appear to be involved in the stabilization of the active form upon phosphorylation and is not essential in propagating the signal downstream for ATPase activity of the central domain. Our data provide experimental evidence against the generally accepted "Y-T coupling" mechanism of activation in NtrC(R).

Synthesis of aromatic 13C/2H-α-ketoacid precursors to be used in selective phenylalanine and tyrosine protein labelling

A synthetic concept to access various ¹³C/²H patterns of phenylalanine and tyrosine precursors for cell-based protein overexpression systems is presented. The target compounds feature isolated ¹³C–¹H spin systems to be used in NMR probing of protein structure and dynamics.

α-Ketoacids as precursors for phenylalanine and tyrosine labelling in cell-based protein overexpression

(13)C-α-ketoacid metabolic precursors of phenylalanine and tyrosine effectively enter the metabolism of a protein overexpressing E. coli strain to label Phe- and Tyr-residues devoid of any cross-labelling. The methodology gives access to highly selective labelling patterns as valuable tools in protein NMR spectroscopy without the need of (15)N-chiral amino acid synthesis using organic chemistry.

Independent valine and leucine isotope labeling in Escherichia coli protein overexpression systems

The addition of labeled α-ketoisovalerate to the growth medium of a protein-expressing host organism has evolved into a versatile tool to achieve concomitant incorporation of specific isotopes into valine- and leucine- residues. The resulting target proteins represent excellent probes for protein NMR analysis. However, as the sidechain resonances of these residues emerge in a narrow spectral range, signal overlap represents a severe limitation in the case of high-molecular-weight NMR probes. We present a protocol to eliminate leucine labeling by supplying the medium with unlabeled α-ketoisocaproate. The resulting spectra of a model protein exclusively feature valine signals of increased intensity, confirming the method to be a first example of independent valine and leucine labeling employing α-ketoacid precursor compounds.