New Views of Functionally Dynamic Proteins by Solution NMR Spectroscopy

Atomic structures of excited state A-T Hoogsteen base pairs in duplex DNA by combining NMR relaxation dispersion, mutagenesis, and chemical shift calculations

NMR relaxation dispersion studies indicate that in canonical duplex DNA, Watson-Crick base pairs (bps) exist in dynamic equilibrium with short-lived low abundance excited state Hoogsteen bps. N1-methylated adenine (m¹A) and guanine (m¹G) are naturally occurring forms of damage that stabilize Hoogsteen bps in duplex DNA. NMR dynamic ensembles of DNA duplexes with m¹A-T Hoogsteen bps reveal significant changes in sugar pucker and backbone angles in and around the Hoogsteen bp, as well as kinking of the duplex towards the major groove. Whether these structural changes also occur upon forming excited state Hoogsteen bps in unmodified duplexes remains to be established because prior relaxation dispersion probes provided limited information regarding the sugar-backbone conformation. Here, we demonstrate measurements of C3' and C4' spin relaxation in the rotating frame (R1ρ) in uniformly ¹³C/¹⁵N labeled DNA as sensitive probes of the sugar-backbone conformation in DNA excited states. The chemical shifts, combined with structure-based predictions using an automated fragmentation quantum mechanics/molecular mechanics method, show that the dynamic ensemble of DNA duplexes containing m¹A-T Hoogsteen bps accurately model the excited state Hoogsteen conformation in two different sequence contexts. Formation of excited state A-T Hoogsteen bps is accompanied by changes in sugar-backbone conformation that allow the flipped syn adenine to form hydrogen-bonds with its partner thymine and this in turn results in overall kinking of the DNA toward the major groove. Results support the assignment of Hoogsteen bps as the excited state observed in canonical duplex DNA, provide an atomic view of DNA dynamics linked to formation of Hoogsteen bps, and lay the groundwork for a potentially general strategy for solving structures of nucleic acid excited states.

Solution NMR Studies of Anesthetic Interactions with Ion Channels

NMR spectroscopy is one of the major tools to provide atomic resolution protein structural information. It has been used to elucidate the molecular details of interactions between anesthetics and ion channels, to identify anesthetic binding sites, and to characterize channel dynamics and changes introduced by anesthetics. In this chapter, we present solution NMR methods essential for investigating interactions between ion channels and general anesthetics, including both volatile and intravenous anesthetics. Case studies are provided with a focus on pentameric ligand-gated ion channels and the voltage-gated sodium channel NaChBac.

Probing the excited-state chemical shifts and exchange parameters by nitrogen-decoupled amide proton chemical exchange saturation transfer (HNdec-CEST)

CEST-NMR spectroscopy is a powerful tool for probing the conformational dynamics of macromolecules. We present a HN^dec-CEST experiment that simplifies the relaxation matrix, reduces fitting parameters, and enhances signal resolution. Importantly, fitting of HN^dec-CEST profiles enables robust extraction of exchange rates as well as excited-state chemical shifts and populations.

Confirmation for conformational selection

NMR studies settle part of a long-standing debate about the mechanism used by the Hsp70 chaperone to recognize substrates.

Slow Dynamics of Tryptophan-Water Networks in Proteins

Water has a profound effect on the dynamics of biomolecules and governs many biological processes, leading to the concept that function is slaved to solvent dynamics within and surrounding the biomolecule. Protein conformational changes on μs-ms time scales are frequently associated with protein function, but little is known about the behavior of protein-bound water on these time scales. Here we have used NMR relaxation dispersion measurements to probe the tryptophan indoles in the enzyme dihydrofolate reductase (DHFR). We find that during structural changes on the μs-ms time scale, large chemical shift changes are often observed for the NH proton on the indole ring, while relatively smaller chemical shift changes are observed for the ring nitrogen atom. Comparison with experimental chemical shifts and density functional theory-based chemical shift predictions show that during the structural change the tryptophan indole NHs remain bound to water, but the geometry of the protein-bound water networks changes. These results establish that relaxation dispersion measurements can indirectly probe water dynamics and indicate that water can influence, or be influenced by, protein conformational changes on the μs-ms time scale. Our data show that structurally conserved bound water molecules can play a critical role in transmitting information between functionally important regions of the protein and provide evidence that internal protein motions can be coupled through the mediation of hydrogen-bonded water bound in the protein structure.

Enhanced spectral density mapping through combined multiple-field deuterium 13CH2D methyl spin relaxation NMR spectroscopy

Quadrupolar relaxation of ²H (D) nuclear spins is a powerful probe of conformational dynamics in biological macromolecules. Deuterium relaxation rate constants are determined by the spectral density function for reorientation of the C-D bond vector at zero, single-quantum, and double-quantum ²H frequencies. In the present work, ²H relaxation rate constants were measured for an E. coli ribonuclease H [U-²H, ¹⁵N] ILV-[¹³CH₂D] sample using 400, 500, 800, and 900 MHz NMR spectrometers and analyzed by three approaches to determine spectral density values. First, data recorded at each static magnetic field were analyzed independently. Second, data recorded at 400 and 800 MHz were analyzed jointly and data recorded at other fields were analyzed independently. Third, data recorded at 400 and 500 MHz were interpolated to 450 MHz, and the resulting two pairs of data, corresponding to 400 MHz/800 MHz and 450 MHz/900 MHz, were analyzed jointly. The second and third approaches rely on the identity between the double quantum frequency at the lower field and the single quantum frequency at the higher field. Spectral density values for 32 of the 48 resolvable ILV methyl resonances were fit by the Lipari-Szabo model-free formalism and used to validate the three methods. The three spectral density mapping methods performed equally well in cross validation with data recorded at 700 MHz. However, the third method yielded approximately 10-15% more precise estimates of model-free parameters and consequently provides a general strategy for analysis of ²H spin relaxation data in biological macromolecules.

15N CEST data and traditional model-free analysis capture fast internal dynamics of DJ-1

Previous studies have shown that relaxation parameters and fast protein dynamics can be quickly elucidated from ¹⁵N-CEST experiments [1]. Longitudinal R₁ and transverse R₂ values were reliably derived from fitting of CEST profiles. Herein we show that ¹⁵N-CEST experiments and traditional modelfree analysis provide the internal dynamics of three states of human protein DJ-1 at physiological temperature. The chemical exchange profiles show the absence of a minor state conformation and, in conjunction with ¹H-¹⁵N NOEs, show increased mobility. R₁ and R₂ values remained relatively unchanged at the three naturally occurring oxidation states of DJ-1, but exhibit striking NOE differences. The NOE data was, therefore, essential in determining the internal motions of the DJ-1 proteins. To the authors' knowledge, we present the first study that combines ¹⁵N CEST data with traditional model-free analyses in the study of a biological system and affirm that more 'lean' model-free approaches should be used cautiously.

An invisible ubiquitin conformation is required for efficient phosphorylation by PINK1

The Ser/Thr protein kinase PINK1 phosphorylates the well-folded, globular protein ubiquitin (Ub) at a relatively protected site, Ser65. We previously showed that Ser65 phosphorylation results in a conformational change in which Ub adopts a dynamic equilibrium between the known, common Ub conformation and a distinct, second conformation wherein the last β-strand is retracted to extend the Ser65 loop and shorten the C-terminal tail. We show using chemical exchange saturation transfer (CEST) nuclear magnetic resonance experiments that a similar, C-terminally retracted (Ub-CR) conformation also exists at low population in wild-type Ub. Point mutations in the moving β5 and neighbouring β-strands shift the Ub/Ub-CR equilibrium. This enabled functional studies of the two states, and we show that while the Ub-CR conformation is defective for conjugation, it demonstrates improved binding to PINK1 through its extended Ser65 loop, and is a superior PINK1 substrate. Together our data suggest that PINK1 utilises a lowly populated yet more suitable Ub-CR conformation of Ub for efficient phosphorylation. Our findings could be relevant for many kinases that phosphorylate residues in folded protein domains.

Comparison of NMR and crystal structures of membrane proteins and computational refinement to improve model quality

Membrane proteins are challenging to study and restraints for structure determination are typically sparse or of low resolution because the membrane environment that surrounds them leads to a variety of experimental challenges. When membrane protein structures are determined by different techniques in different environments, a natural question is "which structure is most biologically relevant?" Towards answering this question, we compiled a dataset of membrane proteins with known structures determined by both solution NMR and X-ray crystallography. By investigating differences between the structures, we found that RMSDs between crystal and NMR structures are below 5 Å in the membrane region, NMR ensembles have a higher convergence in the membrane region, crystal structures typically have a straighter transmembrane region, have higher stereo-chemical correctness, and are more tightly packed. After quantifying these differences, we used high-resolution refinement of the NMR structures to mitigate them, which paves the way for identifying and improving the structural quality of membrane proteins.

Forbidden Coherence Transfer of 19F Nuclei to Quantitatively Measure the Dynamics of a CF₃-Containing Ligand in Receptor-Bound States

The dynamic property of a ligand in the receptor-bound state is an important metric to characterize the interactions in the ligand–receptor interface, and the development of an experimental strategy to quantify the amplitude of motions in the bound state is of importance to introduce the dynamic aspect into structure-guided drug development (SGDD). Fluorine modifications are frequently introduced at the hit-to-lead optimization stage to enhance the binding potency and other characteristics of a ligand. However, the effects of fluorine modifications are generally difficult to predict, owing to the pleiotropic nature of the interactions. In this study, we report an NMR-based approach to experimentally evaluate the local dynamics of trifluoromethyl (CF₃)-containing ligands in the receptor-bound states. For this purpose, the forbidden coherence transfer (FCT) analysis, which has been used to study the dynamics of methyl moieties in proteins, was extended to the ¹⁹F nuclei of CF₃-containing ligands. By applying this CF₃–FCT analysis to a model interaction system consisting of a ligand, AST-487, and a receptor, p38α, we successfully quantified the amplitude of the CF₃ dynamics in the p38α-bound state. The strategy would bring the CF₃-containing ligands within the scope of dynamic SGDD to improve the affinity and specificity for the drug-target receptors.

A Unique and Simple Approach to Improve Sensitivity in 15N-NMR Relaxation Measurements for NH₃⁺ Groups: Application to a Protein-DNA Complex

NMR spectroscopy is a powerful tool for research on protein dynamics. In the past decade, there has been significant progress in the development of NMR methods for studying charged side chains. In particular, NMR methods for lysine side-chain NH₃⁺ groups have been proven to be powerful for investigating the dynamics of hydrogen bonds or ion pairs that play important roles in biological processes. However, relatively low sensitivity has been a major practical issue in NMR experiments on NH₃⁺ groups. In this paper, we present a unique and simple approach to improve sensitivity in ¹⁵N relaxation measurements for NH₃⁺ groups. In this approach, the efficiency of coherence transfers for the desired components are maximized, whereas undesired anti-phase or multi-spin order components are purged through pulse schemes and rapid relaxation. For lysine side-chain NH₃⁺ groups of a protein-DNA complex, we compared the data obtained with the previous and new pulse sequences under the same conditions and confirmed that the ¹⁵N relaxation parameters were consistent for these datasets. While retaining accuracy in measuring ¹⁵N relaxation, our new pulse sequences for NH₃⁺ groups allowed an 82% increase in detection sensitivity of ¹⁵N longitudinal and transverse relaxation measurements.

Complementary uses of small angle X-ray scattering and X-ray crystallography

Most proteins function within networks and, therefore, protein interactions are central to protein function. Although stable macromolecular machines have been extensively studied, dynamic protein interactions remain poorly understood. Small-angle X-ray scattering probes the size, shape and dynamics of proteins in solution at low resolution and can be used to study samples in a large range of molecular weights. Therefore, it has emerged as a powerful technique to study the structure and dynamics of biomolecular systems and bridge fragmented information obtained using high-resolution techniques. Here we review how small-angle X-ray scattering can be combined with other structural biology techniques to study protein dynamics. This article is part of a Special Issue entitled: Biophysics in Canada, edited by Lewis Kay, John Baenziger, Albert Berghuis and Peter Tieleman.

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.

Crowding in Cellular Environments at an Atomistic Level from Computer Simulations

The effects of crowding in biological environments on biomolecular structure, dynamics, and function remain not well understood. Computer simulations of atomistic models of concentrated peptide and protein systems at different levels of complexity are beginning to provide new insights. Crowding, weak interactions with other macromolecules and metabolites, and altered solvent properties within cellular environments appear to remodel the energy landscape of peptides and proteins in significant ways including the possibility of native state destabilization. Crowding is also seen to affect dynamic properties, both conformational dynamics and diffusional properties of macromolecules. Recent simulations that address these questions are reviewed here and discussed in the context of relevant experiments.

The case for defined protein folding pathways

We consider the differences between the many-pathway protein folding model derived from theoretical energy landscape considerations and the defined-pathway model derived from experiment. A basic tenet of the energy landscape model is that proteins fold through many heterogeneous pathways by way of amino acid-level dynamics biased toward selecting native-like interactions. The many pathways imagined in the model are not observed in the structure-formation stage of folding by experiments that would have found them, but they have now been detected and characterized for one protein in the initial prenucleation stage. Analysis presented here shows that these many microscopic trajectories are not distinct in any functionally significant way, and they have neither the structural information nor the biased energetics needed to select native vs. nonnative interactions during folding. The opposed defined-pathway model stems from experimental results that show that proteins are assemblies of small cooperative units called foldons and that a number of proteins fold in a reproducible pathway one foldon unit at a time. Thus, the same foldon interactions that encode the native structure of any given protein also naturally encode its particular foldon-based folding pathway, and they collectively sum to produce the energy bias toward native interactions that is necessary for efficient folding. Available information suggests that quantized native structure and stepwise folding coevolved in ancient repeat proteins and were retained as a functional pair due to their utility for solving the difficult protein folding problem.

Applications of NMR and computational methodologies to study protein dynamics

Overwhelming evidence now illustrates the defining role of atomic-scale protein flexibility in biological events such as allostery, cell signaling, and enzyme catalysis. Over the years, spin relaxation nuclear magnetic resonance (NMR) has provided significant insights on the structural motions occurring on multiple time frames over the course of a protein life span. The present review article aims to illustrate to the broader community how this technique continues to shape many areas of protein science and engineering, in addition to being an indispensable tool for studying atomic-scale motions and functional characterization. Continuing developments in underlying NMR technology alongside software and hardware developments for complementary computational approaches now enable methodologies to routinely provide spatial directionality and structural representations traditionally harder to achieve solely using NMR spectroscopy. In addition to its well-established role in structural elucidation, we present recent examples that illustrate the combined power of selective isotope labeling, relaxation dispersion experiments, chemical shift analyses, and computational approaches for the characterization of conformational sub-states in proteins and enzymes.

Conformational Dynamics and Allostery in E2:E3 Interactions Drive Ubiquitination: gp78 and Ube2g2

Conformational dynamics plays a fundamental role in molecular recognition and activity in enzymes. The ubiquitin-conjugating enzyme (E2) Ube2g2 functions with the ubiquitin ligase (E3) gp78 to assemble poly-ubiquitin chains on target substrates. Two domains in gp78, RING and G2BR, bind to two distant regions of Ube2g2, and activate it for ubiquitin (Ub) transfer. G2BR increases the affinity between the RING and Ube2g2 by 50-fold, while the RING catalyzes the transfer of Ub from the Ube2g2∼Ub conjugate. How G2BR and RING activate Ube2g2 is unclear. In this work, conformational dynamics in Ube2g2 revealed a clear correlation of binding G2BR and RING with the sequential progression toward Ub transfer. The interrelationship of the existence and exchange between ground and excited states leads to a dynamic energy landscape model, in which redistribution of populations contributes to allostery and activation. These findings provide insight into gp78's modulation of conformational exchange in Ube2g2 to stimulate ubiquitination.

The Michaelis Complex of Arginine Kinase Samples the Transition State at a Frequency That Matches the Catalytic Rate

Arginine kinase (AK), which is a member of the phosphagen kinase family, serves as a model system for studying the structural and dynamic determinants of biomolecular enzyme catalysis of all major states involved of the enzymatic cycle. These states are the apo state (substrate free), the Michaelis complex analogue AK:Arg:Mg·AMPPNP (MCA), a product complex analogue AK:pAIE:Mg·ADP (PCA), and the transition state analogue AK:Arg:Mg·ADP:NO₃^- (TSA). The conformational dynamics of these states have been studied by NMR relaxation dispersion measurements of the methyl groups of the Ile, Leu, and Val residues at two static magnetic fields. Although all states undergo significant amounts of μs-ms time scale dynamics, only the MCA samples a dominant excited state that resembles the TSA, as evidenced by the strong correlation between the relaxation dispersion derived chemical shift differences Δω and the equilibrium chemical shift differences Δδ of these states. The average lifetime of the MCA is 36 ms and the free energy difference to the TSA-like form is 8.5 kJ/mol. It is shown that the conformational energy landscape of the Michaelis complex analogue is shaped in a way that at room temperature it channels passage to the transition state, thereby determining the rate-limiting step of the phosphorylation reaction of arginine. Conversely, relaxation dispersion experiments of the TSA reveal that it samples the structures of the Michaelis complex analogue or the apo state as its dominant excited state. This reciprocal behavior shows that the free energy of the TSA, with all ligands bound, is lower by only about 8.9 kJ/mol than that of the Michaelis or apo complex conformations with the TSA ligands present.

Quantitative FRET studies and integrative modeling unravel the structure and dynamics of biomolecular systems

Förster Resonance Energy Transfer (FRET) combined with single-molecule spectroscopy probes macromolecular structure and dynamics and identifies coexisting conformational states. We review recent methodological developments in integrative structural modeling by satisfying spatial restraints on networks of FRET pairs (hybrid-FRET). We discuss procedures to incorporate prior structural knowledge and to obtain optimal distance networks. Finally, a workflow for hybrid-FRET is presented that automates integrative structural modeling and experiment planning to put hybrid-FRET on rails. To test this workflow, we simulate realistic single-molecule experiments and resolve three protein conformers, exchanging at 30μs and 10ms, with accuracies of 1-3Å RMSD versus the target structure. Incorporation of data from other spectroscopies and imaging is also discussed.

Mapping insulin non-covalent interactions with natural polysaccharides by hydrogen/deuterium exchange mass spectrometry

RATIONALE

Drug development efforts involving therapeutic peptides or proteins strongly lead optimization of drug delivery, drug stability, solubility and functionality. The key feature of controlled drug delivery is the use of biocompatible polymers able to interact via non-covalent bonds with an active principle through multiple functional groups. Here amide hydrogen/deuterium exchange (HDX) mass spectrometry was employed to localize insulin dynamics induced by interactions with three natural polysaccharides, i.e. chitosan (CH), sodium alginate (ALG) and chondroitin sulfate (CS).

METHODS

LTQ-Orbitap continuous-labelling mass spectra were collected by diluting insulin stock solution (10 mM in 0.1% formic acid) to a final concentration of 0.1 mM in D₂ O containing 1 mM deuterated ammonium acetate (final pH .6) (insulin:polysaccharide ratio 1:2, w/w). For peptide mapping, deuterated samples were quenched after 0.5, 30, 60, 120 minutes exchange by adding HCl (pH ) and digested with pepsin before LC-MS/MS analysis.

RESULTS

Differences in the insulin backbone dynamics in the presence of the three polysaccharides were highlighted by monitoring peptic peptides at different time points. No significant differences were observed in the presence of CH, whereas the negatively charged ALG and CS were able to induce significant conformational variations at the B-chain level resulting in more protection against H/D exchange. The A-chain interacted only with CS reducing the protein mobility on a long time scale (120 min). HDX data evidenced heterogeneous insulin dynamics in the presence of ALG and CS.

CONCLUSIONS

The studies reported here demonstrated the capabilities of mass spectrometry techniques and HDX methods to obtain useful information toward the flexibility and the behavior of native insulin in the presence of natural polysaccharides, and could provide insights to study the behavior of pharmaceutical formulations. Copyright © 2016 John Wiley & Sons, Ltd.

Coexistence of multiple minor states of fatty acid binding protein and their functional relevance

Proteins are dynamic over a wide range of timescales, but determining the number of distinct dynamic processes and identifying functionally relevant dynamics are still challenging. Here we present the study on human intestinal fatty acid binding protein (hIFABP) using a novel analysis of ¹⁵N relaxation dispersion (RD) and chemical shift saturation transfer (CEST) experiments. Through combined analysis of the two types of experiments, we found that hIFABP exists in a four-state equilibrium in which three minor states interconvert directly with the major state. According to conversion rates from the major “closed” state to minor states, these minor states are irrelevant to the function of fatty acid transport. Based on chemical shifts of the minor states which could not be determined from RD data alone but were extracted from a combined analysis of RD and CEST data, we found that all the minor states are native-like. This conclusion is further supported by hydrogen-deuterium exchange experiments. Direct conversions between the native state and native-like intermediate states may suggest parallel multitrack unfolding/folding pathways of hIFABP. Moreover, hydrogen-deuterium exchange data indicate the existence of another locally unfolded minor state that is relevant to the fatty acid entry process.

High-power (1)H composite pulse decoupling provides artifact free exchange-mediated saturation transfer (EST) experiments

Exchange-mediated saturation transfer (EST) provides critical information regarding dynamics of molecules. In typical applications EST is studied by either scanning a wide range of (15)N chemical shift offsets where the applied (15)N irradiation field strength is on the order of hundreds of Hertz or, scanning a narrow range of (15)N chemical shift offsets where the applied (15)N irradiation field-strength is on the order of tens of Hertz during the EST period. The (1)H decoupling during the EST delay is critical as incomplete decoupling causes broadening of the EST profile, which could possibly result in inaccuracies of the extracted kinetic parameters and transverse relaxation rates. Currently two different (1)H decoupling schemes have been employed, intermittently applied 180° pulses and composite-pulse-decoupling (CPD), for situations where a wide range, or narrow range of (15)N chemical shift offsets are scanned, respectively. We show that high-power CPD provides artifact free EST experiments, which can be universally implemented regardless of the offset range or irradiation field-strengths.

Intrinsic local destabilization of the C-terminus predisposes integrin α1 I domain to a conformational switch induced by collagen binding

Integrin-collagen interactions play a critical role in a myriad of cellular functions that include immune response, and cell development and differentiation, yet their mechanism of binding is poorly understood. There is increasing evidence that conformational flexibility assumes a central role in the molecular mechanisms of protein-protein interactions and here we employ NMR hydrogen-deuterium exchange (HDX) experiments to explore the impact of slower timescale dynamic events. To gain insight into the mechanisms underlying collagen-induced conformational switches, we have undertaken a comparative study between the wild type integrin α1 I and a gain-of-function E317A mutant. NMR HDX results suggest a relationship between regions exhibiting a reduced local stability in the unbound I domain and those that undergo significant conformational changes upon binding. Specifically, the αC and α7 helices within the C-terminus are at the center of such major perturbations and present reduced local stabilities in the unbound state relative to other structural elements. Complementary isothermal titration calorimetry experiments have been performed to derive complete thermodynamic binding profiles for association of the collagen-like triple-helical peptide with wild type α1 I and E317A mutant. The differential energetics observed for E317A are consistent with the HDX experiments and support a model in which intrinsically destabilized regions predispose conformational rearrangement in the integrin I domain. This study highlights the importance of exploring different timescales to delineate allosteric and binding events.

Conformational Sub-states and Populations in Enzyme Catalysis

Enzyme function involves substrate and cofactor binding, precise positioning of reactants in the active site, chemical turnover, and release of products. In addition to formation of crucial structural interactions between enzyme and substrate(s), coordinated motions within the enzyme-substrate complex allow reaction to proceed at a much faster rate, compared to the reaction in solution and in the absence of enzyme. An increasing number of enzyme systems show the presence of conserved protein motions that are important for function. A wide variety of motions are naturally sampled (over femtosecond to millisecond time-scales) as the enzyme complex moves along the energetic landscape, driven by temperature and dynamical events from the surrounding environment. Areas of low energy along the landscape form conformational sub-states, which show higher conformational populations than surrounding areas. A small number of these protein conformational sub-states contain functionally important structural and dynamical features, which assist the enzyme mechanism along the catalytic cycle. Identification and characterization of these higher-energy (also called excited) sub-states and the associated populations are challenging, as these sub-states are very short-lived and therefore rarely populated. Specialized techniques based on computer simulations, theoretical modeling, and nuclear magnetic resonance have been developed for quantitative characterization of these sub-states and populations. This chapter discusses these techniques and provides examples of their applications to enzyme systems.

Visualizing chaperone-assisted protein folding

Challenges in determining the structures of heterogeneous and dynamic protein complexes have greatly hampered past efforts to obtain a mechanistic understanding of many important biological processes. One such process is chaperone-assisted protein folding. Obtaining structural ensembles of chaperone-substrate complexes would ultimately reveal how chaperones help proteins fold into their native state. To address this problem, we devised a new structural biology approach based on X-ray crystallography, termed residual electron and anomalous density (READ). READ enabled us to visualize even sparsely populated conformations of the substrate protein immunity protein 7 (Im7) in complex with the Escherichia coli chaperone Spy, and to capture a series of snapshots depicting the various folding states of Im7 bound to Spy. The ensemble shows that Spy-associated Im7 samples conformations ranging from unfolded to partially folded to native-like states and reveals how a substrate can explore its folding landscape while being bound to a chaperone.

STARD6 on steroids: solution structure, multiple timescale backbone dynamics and ligand binding mechanism

START domain proteins are conserved α/β helix-grip fold that play a role in the non-vesicular and intracellular transport of lipids and sterols. The mechanism and conformational changes permitting the entry of the ligand into their buried binding sites is not well understood. Moreover, their functions and the identification of cognate ligands is still an active area of research. Here, we report the solution structure of STARD6 and the characterization of its backbone dynamics on multiple time-scales through ¹⁵N spin-relaxation and amide exchange studies. We reveal for the first time the presence of concerted fluctuations in the Ω₁ loop and the C-terminal helix on the microsecond-millisecond time-scale that allows for the opening of the binding site and ligand entry. We also report that STARD6 binds specifically testosterone. Our work represents a milestone for the study of ligand binding mechanism by other START domains and the elucidation of the biological function of STARD6.

Biophysical and computational methods to analyze amino acid interaction networks in proteins

Globular proteins are held together by interacting networks of amino acid residues. A number of different structural and computational methods have been developed to interrogate these amino acid networks. In this review, we describe some of these methods, including analyses of X-ray crystallographic data and structures, computer simulations, NMR data, and covariation among protein sequences, and indicate the critical insights that such methods provide into protein function. This information can be leveraged towards the design of new allosteric drugs, and the engineering of new protein function and protein regulation strategies.

Chemo-enzymatic labeling for rapid assignment of RNA molecules

Even though Nuclear Magnetic Resonance (NMR) spectroscopy is one of the few techniques capable of determining atomic resolution structures of RNA, it is constrained by two major problems of chemical shift overlap of resonances and rapid signal loss due to line broadening. Emerging tools to tackle these problems include synthesis of atom specifically labeled or chemically modified nucleotides. Herein we review the synthesis of these nucleotides, the design and production of appropriate RNA samples, and the application and analysis of the NMR experiments that take advantage of these labels.

Optimization of Conformational Dynamics in an Epistatic Evolutionary Trajectory

The understanding of protein evolution depends on the ability to relate the impact of mutations on molecular traits to organismal fitness. Biological activity and robustness have been regarded as important features in shaping protein evolutionary landscapes. Conformational dynamics, which is essential for protein function, has received little attention in the context of evolutionary analyses. Here we employ NMR spectroscopy, the chief experimental tool to describe protein dynamics at atomic level in solution at room temperature, to study the intrinsic dynamic features of a metallo- Β: -lactamase enzyme and three variants identified during a directed evolution experiment that led to an expanded substrate profile. We show that conformational dynamics in the catalytically relevant microsecond to millisecond timescale is optimized along the favored evolutionary trajectory. In addition, we observe that the effects of mutations on dynamics are epistatic. Mutation Gly262Ser introduces slow dynamics on several residues that surround the active site when introduced in the wild-type enzyme. Mutation Asn70Ser removes the slow dynamics observed for few residues of the wild-type enzyme, but increases the number of residues that undergo slow dynamics when introduced in the Gly262Ser mutant. These effects on dynamics correlate with the epistatic interaction between these two mutations on the bacterial phenotype. These findings indicate that conformational dynamics is an evolvable trait, and that proteins endowed with more dynamic active sites also display a larger potential for promoting evolution.