Dependence of amino acid side chain 13C shifts on dihedral angle: application to conformational analysis

Structural basis of substrate recognition and allosteric activation of the proapoptotic mitochondrial HtrA2 protease

The mitochondrial serine protease HtrA2 is a human homolog of the Escherichia coli Deg-proteins exhibiting chaperone and proteolytic roles. HtrA2 is involved in both apoptotic regulation via its ability to degrade inhibitor-of-apoptosis proteins (IAPs), as well as in cellular maintenance as part of the cellular protein quality control machinery, by preventing the possible toxic accumulation of aggregated proteins. In this study, we use advanced solution NMR spectroscopy methods combined with biophysical characterization and biochemical assays to elucidate the crucial role of the substrate recognizing PDZ domain. This domain regulates the protease activity of HtrA2 by triggering an intricate allosteric network involving the regulatory loops of the protease domain. We further show that divalent metal ions can both positively and negatively modulate the activity of HtrA2, leading to a refined model of HtrA2 regulation within the apoptotic pathway.

Probing amino acid side chains of the integral membrane protein PagP by solution NMR: Side chain immobilization facilitates association of secondary structures

Solution NMR spectroscopy of large protein systems is hampered by rapid signal decay, so most multidimensional studies focus on long-lived 1H-13C magnetization in methyl groups and/or backbone amide 1H-15N magnetization in an otherwise perdeuterated environment. Herein we demonstrate that it is possible to biosynthetically incorporate additional 1H-12C groups that possess long-lived magnetization using cost-effective partially deuterated or unlabeled amino acid precursors added to Escherichia coli growth media. This approach is applied to the outer membrane enzyme PagP in membrane-mimetic dodecylphosphocholine micelles. We were able to obtain chemical shift assignments for a majority of side chain 1H positions in PagP using nuclear Overhauser enhancements (NOEs) to connect them to previously assigned backbone 1H-15N groups and newly assigned 1H-13C methyl groups. Side chain methyl-to-aromatic NOEs were particularly important for confirming that the amphipathic α-helix of PagP packs against its eight-stranded β-barrel, as indicated by previous X-ray crystal structures. Interestingly, aromatic NOEs suggest that some aromatic residues in PagP that are buried in the membrane bilayer are highly mobile in the micellar environment, like Phe138 and Phe159. In contrast, Tyr87 in the middle of the bilayer is quite rigid, held in place by a hydrogen bonded network extending to the surface that resembles a classic catalytic triad: Tyr87-His67-Asp61. This hydrogen bonded arrangement of residues is not known to have any catalytic activity, but we postulate that its role is to immobilize Tyr87 to facilitate packing of the amphipathic α-helix against the β-barrel.

Structure and dynamics of the mitochondrial DNA-compaction factor Abf2 from S. cerevisiae

Mitochondria are essential organelles that produce most of the energy via the oxidative phosphorylation (OXPHOS) system in all eukaryotic cells. Several essential subunits of the OXPHOS system are encoded by the mitochondrial genome (mtDNA) despite its small size. Defects in mtDNA maintenance and expression can lead to severe OXPHOS deficiencies and are amongst the most common causes of mitochondrial disease. The mtDNA is packaged as nucleoprotein structures, referred to as nucleoids. The conserved mitochondrial proteins, ARS-binding factor 2 (Abf2) in the budding yeast Saccharomyces cerevisiae (S. cerevisiae) and mitochondrial transcription factor A (TFAM) in mammals, are nucleoid associated proteins (NAPs) acting as condensing factors needed for packaging and maintenance of the mtDNA. Interestingly, gene knockout of Abf2 leads, in yeast, to the loss of mtDNA and respiratory growth, providing evidence for its crucial role. On a structural level, the condensing factors usually contain two DNA binding domains called high-mobility group boxes (HMG boxes). The co-operating mechanical activities of these domains are crucial in establishing the nucleoid architecture by bending the DNA strands. Here we used advanced solution NMR spectroscopy methods to characterize the dynamical properties of Abf2 together with its interaction with DNA. We find that the two HMG-domains react notably different to external cues like temperature and salt, indicating distinct functional properties. Biophysical characterizations show the pronounced difference of these domains upon DNA-binding, suggesting a refined picture of the Abf2 functional cycle.

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.

Conformational selection guides β-arrestin recruitment at a biased G protein-coupled receptor

G protein-coupled receptors (GPCRs) recruit β-arrestins to coordinate diverse cellular processes, but the structural dynamics driving this process are poorly understood. Atypical chemokine receptors (ACKRs) are intrinsically biased GPCRs that engage β-arrestins but not G proteins, making them a model system for investigating the structural basis of β-arrestin recruitment. Here, we performed nuclear magnetic resonance (NMR) experiments on ¹³CH₃-ε-methionine-labeled ACKR3, revealing that β-arrestin recruitment is associated with conformational exchange at key regions of the extracellular ligand-binding pocket and intracellular β-arrestin-coupling region. NMR studies of ACKR3 mutants defective in β-arrestin recruitment identified an allosteric hub in the receptor core that coordinates transitions among heterogeneously populated and selected conformational states. Our data suggest that conformational selection guides β-arrestin recruitment by tuning receptor dynamics at intracellular and extracellular regions.

Methodological advancements for characterising protein side chains by NMR spectroscopy

The surface of proteins is covered by side chains of polar amino acids that are imperative for modulating protein functionality through the formation of noncovalent intermolecular interactions. However, despite their tremendous importance, the unique structures of protein side chains require tailored approaches for investigation by nuclear magnetic resonance spectroscopy and so have traditionally been understudied compared with the protein backbone. Here, we review substantial recent methodological advancements within nuclear magnetic resonance spectroscopy to address this issue. Specifically, we consider advancements that provide new insight into methyl-bearing side chains, show the potential of using non-natural amino acids and reveal the actions of charged side chains. Combined, the new methods promise unprecedented characterisations of side chains that will further elucidate protein function.

Targeting the cryptic sites: NMR-based strategy to improve protein druggability by controlling the conformational equilibrium

Cryptic ligand binding sites, which are not evident in the unligated structures, are beneficial in tackling with difficult but attractive drug targets, such as protein-protein interactions (PPIs). However, cryptic sites have thus far not been rationally pursued in the early stages of drug development. Here, we demonstrated by nuclear magnetic resonance that the cryptic site in Bcl-xL exists in a conformational equilibrium between the open and closed conformations under the unligated condition. While the fraction of the open conformation in the unligated wild-type Bcl-xL is estimated to be low, F143W mutation that is distal from the ligand binding site can substantially elevate the population. The F143W mutant showed a higher hit rate in a phage-display peptide screening, and the hit peptide bound to the cryptic site of the wild-type Bcl-xL. Therefore, by controlling the conformational equilibrium in the cryptic site, the opportunity to identify a PPI inhibitor could be improved.

Determining isoleucine side-chain rotamer-sampling in proteins from 13C chemical shift

A framework is presented to derive the conformational sampling of isoleucine side chains from nuclear magnetic resonance ¹³C chemical shifts.

Chemical shifts are often the only nuclear magnetic resonance parameter that can be obtained for challenging macromolecular systems. Here we present a framework to derive the conformational sampling of isoleucine side chains from ¹³C chemical shifts and demonstrate that side-chain conformations in a low-populated folding intermediate can be determined.

Rotamer Jumps, Proton Exchange, and Amine Inversion Dynamics of Dimethylated Lysine Residues in Proteins Resolved by pH-Dependent 1H and 13C NMR Relaxation Dispersion

Post-translational methylation of lysine side chains is of great importance for protein regulation, including epigenetic control. Here, we present specific ¹³CHD₂ labeling of dimethylated lysines as a sensitive probe of the structure, interactions, and dynamics of these groups, and outline a theoretical and experimental framework for analyzing their conformational dynamics using ¹H and ¹³C CPMG relaxation dispersion experiments. Dimethylated lysine side chains in calcium-loaded calmodulin show a marked pH dependence of their Carr-Purcell-Meiboom-Gill (CPMG) dispersion profiles, indicating complex exchange behavior. Combined analysis of ¹H and ¹³C CPMG relaxation dispersions requires consideration of 12-state correlated exchange of the two methyl groups due to circular three-state rotamer jumps around the Cε-Nζ axis combined with proton exchange and amine inversion. Taking into account a number of fundamental constraints, the exchange model can be reduced to include only three fitted parameters, namely, the geometric average of the rotamer-jump rate constants, the rate constant of deprotonation of Nζ, and the chemical shift difference between the trans and gauge positions of the ¹³C or ¹H nuclei. The pH dependence indicates that protonation of the end group dramatically slows down rotamer exchange for some lysine residues, whereas deprotonation leads to rapid amine inversion coupled with rotamer scrambling. The observed variation among residues in their exchange behavior appears to depend on the structural environment of the side chain. Understanding this type of exchange process is critical to correctly interpreting NMR spectra of methylated lysine side chains. The exchange model presented here forms the basis for studying the structure and dynamics of epigenetically modified lysine side chains and perturbations caused by changes in pH or interactions with target proteins.

Transitions in DNA polymerase β μs-ms dynamics related to substrate binding and catalysis

DNA polymerase β (pol β) plays a central role in the DNA base excision repair pathway and also serves as an important model polymerase. Dynamic characterization of pol β from methyl-TROSY 13C-1H multiple quantum CPMG relaxation dispersion experiments of Ile and Met sidechains and previous backbone relaxation dispersion measurements, reveals transitions in μs-ms dynamics in response to highly variable substrates. Recognition of a 1-nt-gapped DNA substrate is accompanied by significant backbone and sidechain motion in the lyase domain and the DNA binding subdomain of the polymerase domain, that may help to facilitate binding of the apoenzyme to the segments of the DNA upstream and downstream from the gap. Backbone μs-ms motion largely disappears after formation of the pol β-DNA complex, giving rise to an increase in uncoupled μs-ms sidechain motion throughout the enzyme. Formation of an abortive ternary complex using a non-hydrolyzable dNTP results in sidechain motions that fit to a single exchange process localized to the catalytic subdomain, suggesting that this motion may play a role in catalysis.

Computational insights on agonist and antagonist mechanisms of estrogen receptor α induced by bisphenol A analogues

Structural analogues of bisphenol A (BPA) have become widely used as alternatives in BPA-free products. Most toxicological investigations have focused on the estrogenic activities of these analogues, which have been considered as potential environmental estrogens. However, recent studies revealed that certain BPA analogues could dramatically inhibit the proliferation of breast cancer cells, and exhibited strong anti-estrogenic effects compared with the antagonist 4-hydroxytamoxifen (OHT). Thus, we adopted computational models combining molecular dynamics simulations and binding free energy calculations to explore the underlying molecular basis of BPA analogues binding to estrogen receptor α (ERα). We also evaluated ligand-induced structural rearrangements of ERα at the atomic level. Conformational analyses showed that induced-fit H-bonding recognition by Thr347 was an important factor distinguishing antagonist from agonist BPA analogues. Moreover, antagonists of BPA analogues could indirectly induce the structural reposition of key helix 12 and produce an antagonistic conformation of ERα. Compared with OHT, the binding affinity of BPA analogues is stronger for antagonists than agonists. Taken together, we therefore propose computational indicators for screening of anti-estrogenic activities of BPA analogues, which may be beneficial for predicting the estrogenic or anti-estrogenic effects of BPA alternatives.

Formation of β-Lactoglobulin Aggregates from Quite, Unfolded Conformations upon Heat Activation

In presence of calcium ions, β-lactoglobulin (BLG) unfolds and subsequently aggregates after heating. This process has important pharmaceutical and agroalimentary applications. Nowadays, the molecular mechanism of unfolding and BLG aggregation, and the role of calcium in the mechanism, is poorly understood. Actually, in most studies, data have been acquired at room temperature, after heating and after aggregation, which makes it difficult to establish a clear causal-temporal relation between calcium binding, heat, and aggregation. Thus, the goal of the present study is to get accurate, nanoscale data about the molecular events leading to BLG unfolding and calcium-dependent aggregation. The molecular transformation of BLG during heating has been investigated, using the NMR pulse field gradient technique, operating in a high field (900 MHz). Thanks to this technique, the molecular conformation of newly formed unfolded BLG molecules can be distinguished in a large pool of native ones. The present work shows that BLG at neutral pH at 65 °C displays fast, cooperative-like unfolding, in which no long-lived intermediary state (as a molten globule one) is detected, before aggregation. These data also indicate that calcium ions bind unfolded BLG in specific sites which might be a necessary feature to form the aggregate. Finally, these data also provide an NMR-based methodology to monitor the rate of protein unfolding using NMR.

Solution NMR spectroscopy of GPCRs: Residue-specific labeling strategies with a focus on 13C-methyl methionine labeling of the atypical chemokine receptor ACKR3

The past decade has witnessed remarkable progress in the determination of G protein-coupled receptor (GPCR) structures, profoundly expanding our understanding of how GPCRs recognize ligands, become activated, and interact with intracellular signaling components. In recent years, numerous studies have used solution nuclear magnetic resonance (NMR) spectroscopy to investigate GPCRs, providing fundamental insights into GPCR conformational changes, allostery, dynamics, and other facets of GPCR function are challenging to study using other structural techniques. Despite these advantages, NMR-based studies of GPCRs are few relative to the number of published structures, due in part to the challenges and limitations of NMR for the characterization of large membrane proteins. Several studies have circumvented these challenges using a variety of isotopic labeling strategies, including side chain derivatization and metabolic incorporation of NMR-active nuclei. In this chapter, we provide an overview of different isotopic labeling strategies and describe an in-depth protocol for the expression, purification, and NMR studies of the chemokine GPCR atypical chemokine receptor 3 (ACKR3) via ¹³CH₃-methionine incorporation. The goal of this chapter is to provide a resource to the GPCR community for those interested in pursuing NMR studies of GPCRs.

GPCR drug discovery: integrating solution NMR data with crystal and cryo-EM structures

The 826 G protein-coupled receptors (GPCRs) in the human proteome regulate key physiological processes and thus have long been attractive drug targets. With the crystal structures of more than 50 different human GPCRs determined over the past decade, an initial platform for structure-based rational design has been established for drugs that target GPCRs, which is currently being augmented with cryo-electron microscopy (cryo-EM) structures of higher-order GPCR complexes. Nuclear magnetic resonance (NMR) spectroscopy in solution is one of the key approaches for expanding this platform with dynamic features, which can be accessed at physiological temperature and with minimal modification of the wild-type GPCR covalent structures. Here, we review strategies for the use of advanced biochemistry and NMR techniques with GPCRs, survey projects in which crystal or cryo-EM structures have been complemented with NMR investigations and discuss the impact of this integrative approach on GPCR biology and drug discovery.

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.

Phosphorylation-induced conformation of β2-adrenoceptor related to arrestin recruitment revealed by NMR

The C-terminal region of G-protein-coupled receptors (GPCRs), stimulated by agonist binding, is phosphorylated by GPCR kinases, and the phosphorylated GPCRs bind to arrestin, leading to the cellular responses. To understand the mechanism underlying the formation of the phosphorylated GPCR-arrestin complex, we performed NMR analyses of the phosphorylated β₂-adrenoceptor (β₂AR) and the phosphorylated β₂AR–β-arrestin 1 complex, in the lipid bilayers of nanodisc. Here we show that the phosphorylated C-terminal region adheres to either the intracellular side of the transmembrane region or lipids, and that the phosphorylation of the C-terminal region allosterically alters the conformation around M215^5.54 and M279^6.41, located on transemembrane helices 5 and 6, respectively. In addition, we found that the conformation induced by the phosphorylation is similar to that corresponding to the β-arrestin-bound state. The phosphorylation-induced structures revealed in this study propose a conserved structural motif of GPCRs that enables β-arrestin to recognize dozens of GPCRs.

Upon stimulation by agonist binding, the C-terminal regions of G-protein-coupled receptors (GPCRs) become phosphorylated by GPCR kinases, and phosphorylated GPCRs bind arrestin. Here the authors give structural insights into the phosphorylation induced conformational changes in GPCRs by performing NMR studies with the β2-adrenoceptor.

Pressure dependence of side chain 13C chemical shifts in model peptides Ac-Gly-Gly-Xxx-Ala-NH2

For evaluating the pressure responses of folded as well as intrinsically unfolded proteins detectable by NMR spectroscopy the availability of data from well-defined model systems is indispensable. In this work we report the pressure dependence of ¹³C chemical shifts of the side chain atoms in the protected tetrapeptides Ac-Gly-Gly-Xxx-Ala-NH₂ (Xxx, one of the 20 canonical amino acids). Contrary to expectation the chemical shifts of a number of nuclei have a nonlinear dependence on pressure in the range from 0.1 to 200 MPa. The size of the polynomial pressure coefficients B ₁ and B ₂ is dependent on the type of atom and amino acid studied. For H^N, N and C^α the first order pressure coefficient B ₁ is also correlated to the chemical shift at atmospheric pressure. The first and second order pressure coefficients of a given type of carbon atom show significant linear correlations suggesting that the NMR observable pressure effects in the different amino acids have at least partly the same physical cause. In line with this observation the magnitude of the second order coefficients of nuclei being direct neighbors in the chemical structure also are weakly correlated. The downfield shifts of the methyl resonances suggest that gauche conformers of the side chains are not preferred with pressure. The valine and leucine methyl groups in the model peptides were assigned using stereospecifically ¹³C enriched amino acids with the pro-R carbons downfield shifted relative to the pro-S carbons.

Recent developments in solution nuclear magnetic resonance (NMR)-based molecular biology

Visualizing post-translational modifications, conformations, and interaction surfaces of protein structures at atomic resolution underpins the development of novel therapeutics to combat disease. As computational resources expand, in silico calculations coupled with experimentally derived structures and functional assays have led to an explosion in structure-based drug design (SBDD) with several compounds in clinical trials. It is increasingly clear that "hidden" transition-state structures along activation trajectories can be harnessed to develop novel classes of allosteric inhibitors. The goal of this mini-review is to empower the clinical researcher with a general knowledge of the strengths and weaknesses of nuclear magnetic resonance (NMR) spectroscopy in molecular medicine. Although NMR can determine protein structures at atomic resolution, its unrivaled strength lies in sensing subtle changes in a nuclei's chemical environment as a result of intrinsic conformational dynamics, solution conditions, and binding interactions. These can be recorded at atomic resolution, without explicit structure determination, and then incorporated with static structures or molecular dynamics simulations to produce a complete biological picture.

The application of tailor-made force fields and molecular dynamics for NMR crystallography: a case study of free base cocaine

Motional averaging has been proven to be significant in predicting the chemical shifts in ab initio solid-state NMR calculations, and the applicability of motional averaging with molecular dynamics has been shown to depend on the accuracy of the molecular mechanical force field. The performance of a fully automatically generated tailor-made force field (TMFF) for the dynamic aspects of NMR crystallography is evaluated and compared with existing benchmarks, including static dispersion-corrected density functional theory calculations and the COMPASS force field. The crystal structure of free base cocaine is used as an example. The results reveal that, even though the TMFF outperforms the COMPASS force field for representing the energies and conformations of predicted structures, it does not give significant improvement in the accuracy of NMR calculations. Further studies should direct more attention to anisotropic chemical shifts and development of the method of solid-state NMR calculations.

Dynamic regulation of GDP binding to G proteins revealed by magnetic field-dependent NMR relaxation analyses

Heterotrimeric guanine-nucleotide-binding proteins (G proteins) serve as molecular switches in signalling pathways, by coupling the activation of cell surface receptors to intracellular responses. Mutations in the G protein α-subunit (Gα) that accelerate guanosine diphosphate (GDP) dissociation cause hyperactivation of the downstream effector proteins, leading to oncogenesis. However, the structural mechanism of the accelerated GDP dissociation has remained unclear. Here, we use magnetic field-dependent nuclear magnetic resonance relaxation analyses to investigate the structural and dynamic properties of GDP bound Gα on a microsecond timescale. We show that Gα rapidly exchanges between a ground-state conformation, which tightly binds to GDP and an excited conformation with reduced GDP affinity. The oncogenic D150N mutation accelerates GDP dissociation by shifting the equilibrium towards the excited conformation.

Heterotrimeric guanine-nucleotide-binding proteins (G proteins) act as molecular switches. Here the authors use NMR relaxation analyses, which reveal the dynamics of G protein alpha subunit binding to GDP on a microsecond timescale.

Metal Ion Binding at the Catalytic Site Induces Widely Distributed Changes in a Sequence Specific Protein-DNA Complex

Metal ion cofactors can alter the energetics and specificity of sequence specific protein–DNA interactions, but it is unknown if the underlying effects on structure and dynamics are local or dispersed throughout the protein–DNA complex. This work uses EcoRV endonuclease as a model, and catalytically inactive lanthanide ions, which replace the Mg²⁺ cofactor. Nuclear magnetic resonance (NMR) titrations indicate that four Lu³⁺ or two La³⁺ cations bind, and two new crystal structures confirm that Lu³⁺ binding is confined to the active sites. NMR spectra show that the metal-free EcoRV complex with cognate (GATATC) DNA is structurally distinct from the nonspecific complex, and that metal ion binding sites are not assembled in the nonspecific complex. NMR chemical shift perturbations were determined for ¹H–¹⁵N amide resonances, for ¹H–¹³C Ile-δ-CH₃ resonances, and for stereospecifically assigned Leu-δ-CH₃ and Val-γ-CH₃ resonances. Many chemical shifts throughout the cognate complex are unperturbed, so metal binding does not induce major conformational changes. However, some large perturbations of amide and side chain methyl resonances occur as far as 34 Å from the metal ions. Concerted changes in specific residues imply that local effects of metal binding are propagated via a β-sheet and an α-helix. Both amide and methyl resonance perturbations indicate changes in the interface between subunits of the EcoRV homodimer. Bound metal ions also affect amide hydrogen exchange rates for distant residues, including a distant subdomain that contacts DNA phosphates and promotes DNA bending, showing that metal ions in the active sites, which relieve electrostatic repulsion between protein and DNA, cause changes in slow dynamics throughout the complex.

NMR Methods to Study Dynamic Allostery

Nuclear magnetic resonance (NMR) spectroscopy provides a unique toolbox of experimental probes for studying dynamic processes on a wide range of timescales, ranging from picoseconds to milliseconds and beyond. Along with NMR hardware developments, recent methodological advancements have enabled the characterization of allosteric proteins at unprecedented detail, revealing intriguing aspects of allosteric mechanisms and increasing the proportion of the conformational ensemble that can be observed by experiment. Here, we present an overview of NMR spectroscopic methods for characterizing equilibrium fluctuations in free and bound states of allosteric proteins that have been most influential in the field. By combining NMR experimental approaches with molecular simulations, atomistic-level descriptions of the mechanisms by which allosteric phenomena take place are now within reach.

NMR Method for Characterizing Microsecond-to-Millisecond Chemical Exchanges Utilizing Differential Multiple-Quantum Relaxation in High Molecular Weight Proteins

Chemical exchange processes of proteins on the order of microseconds (μs) to milliseconds (ms) play critical roles in biological functions. Developments in methyl-transverse relaxation optimized spectroscopy (methyl-TROSY), which observes the slowly relaxing multiple quantum (MQ) coherences, have enabled the studies of biologically important large proteins. However, the analyses of μs to ms chemical exchange processes based on the methyl-TROSY principle are still challenging, because the interpretation of the chemical exchange contributions to the MQ relaxation profiles is complicated, as significant chemical shift differences occur in both (1)H and (13)C nuclei. Here, we report a new methyl-based NMR method for characterizing chemical exchanges, utilizing differential MQ relaxation rates and a heteronuclear double resonance pulse technique. The method enables quantitative evaluations of the chemical exchange processes, in which significant chemical shift differences exist in both the (1)H and (13)C nuclei. The versatility of the method is demonstrated with the application to KirBac1.1, with an apparent molecular mass of 200 kDa.

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.

Side Chain Conformational Distributions of a Small Protein Derived from Model-Free Analysis of a Large Set of Residual Dipolar Couplings

Accurate quantitative measurement of structural dispersion in proteins remains a prime challenge to both X-ray crystallography and NMR spectroscopy. Here we use a model-free approach based on measurement of many residual dipolar couplings (RDCs) in differentially orienting aqueous liquid crystalline solutions to obtain the side chain χ₁ distribution sampled by each residue in solution. Applied to the small well-ordered model protein GB3, our approach reveals that the RDC data are compatible with a single narrow distribution of side chain χ₁ angles for only about 40% of the residues. For more than half of the residues, populations greater than 10% for a second rotamer are observed, and four residues require sampling of three rotameric states to fit the RDC data. In virtually all cases, sampled χ₁ values are found to center closely around ideal g^–, g⁺ and t rotameric angles, even though no rotamer restraint is used when deriving the sampled angles. The root-mean-square difference between experimental ³J_HαHβ couplings and those predicted by the Haasnoot-parametrized, motion-adjusted Karplus equation reduces from 2.05 to 0.75 Hz when using the new rotamer analysis instead of the 1.1-Å X-ray structure as input for the dihedral angles. A comparison between observed and predicted ³J_HαHβ values suggests that the root-mean-square amplitude of χ₁ angle fluctuations within a given rotamer well is ca. 20°. The quantitatively defined side chain rotamer equilibria obtained from our study set new benchmarks for evaluating improved molecular dynamics force fields, and also will enable further development of quantitative relations between side chain chemical shift and structure.

Equilibrium transitions between side-chain conformations in leucine and isoleucine

Despite recent improvements in computational methods for protein design, we still lack a quantitative, predictive understanding of the intrinsic probabilities for amino acids to adopt particular side-chain conformations. Surprisingly, this question has remained unsettled for many years, in part because of inconsistent results from different experimental approaches. To explicitly determine the relative populations of different side-chain dihedral angles, we performed all-atom hard-sphere Langevin Dynamics simulations of leucine (Leu) and isoleucine (Ile) dipeptide mimetics with stereo-chemical constraints and repulsive-only steric interactions between non-bonded atoms. We determine the relative populations of the different χ(1) and χ(2) dihedral angle combinations as a function of the backbone dihedral angles ϕ and ψ. We also propose, and test, a mechanism for inter-conversion between the different side-chain conformations. Specifically, we discover that some of the transitions between side-chain dihedral angle combinations are very frequent, whereas others are orders of magnitude less frequent, because they require rare coordinated motions to avoid steric clashes. For example, to transition between different values of χ(2), the Leu side-chain bond angles κ(1) and κ(2) must increase, whereas to transition in χ(1), the Ile bond angles λ(1) and λ(2) must increase. These results emphasize the importance of computational approaches in stimulating further experimental studies of the conformations of side-chains in proteins. Moreover, our studies emphasize the power of simple steric models to inform our understanding of protein structure, dynamics, and design.

Dynamic multidrug recognition by multidrug transcriptional repressor LmrR

LmrR is a multidrug transcriptional repressor that controls the expression of a major multidrug transporter, LmrCD, in Lactococcus lactis. However, the molecular mechanism by which LmrR binds to structurally unrelated compounds and is released from the promoter region remains largely unknown. Here, we structurally and dynamically characterized LmrR in the apo, compound-bound and promoter-bound states. The compound-binding site of LmrR exhibits ps–μs dynamics in the apo state, and compound ligation shifts the preexisting conformational equilibrium to varying extents to achieve multidrug recognition. Meanwhile, the compound binding induces redistribution of ps–ns dynamics to the allosteric sites, which entropically favors the high-affinity recognition. Furthermore, the reciprocal compound/promoter binding by LmrR is achieved by the incompatible conformational ensembles between the compound- and promoter-bound states. Collectively, the data show how LmrR can dynamically exert its functions through promiscuous multi-target interactions, in a manner that cannot be understood by a static structural view.

Functional dynamics of cell surface membrane proteins

Cell surface receptors are integral membrane proteins that receive external stimuli, and transmit signals across plasma membranes. In the conventional view of receptor activation, ligand binding to the extracellular side of the receptor induces conformational changes, which convert the structure of the receptor into an active conformation. However, recent NMR studies of cell surface membrane proteins have revealed that their structures are more dynamic than previously envisioned, and they fluctuate between multiple conformations in an equilibrium on various timescales. In addition, NMR analyses, along with biochemical and cell biological experiments indicated that such dynamical properties are critical for the proper functions of the receptors. In this review, we will describe several NMR studies that revealed direct linkage between the structural dynamics and the functions of the cell surface membrane proteins, such as G-protein coupled receptors (GPCRs), ion channels, membrane transporters, and cell adhesion molecules.

Selective unfolding of one Ribonuclease H domain of HIV reverse transcriptase is linked to homodimer formation

HIV-1 reverse transcriptase (RT), a critical enzyme of the HIV life cycle and an important drug target, undergoes complex and largely uncharacterized conformational rearrangements that underlie its asymmetric folding, dimerization and subunit-selective ribonuclease H domain (RH) proteolysis. In the present article we have used a combination of NMR spectroscopy, small angle X-ray scattering and X-ray crystallography to characterize the p51 and p66 monomers and the conformational maturation of the p66/p66′ homodimer. The p66 monomer exists as a loosely structured molecule in which the fingers/palm/connection, thumb and RH substructures are connected by flexible (disordered) linking segments. The initially observed homodimer is asymmetric and includes two fully folded RH domains, while exhibiting other conformational features similar to that of the RT heterodimer. The RH′ domain of the p66′ subunit undergoes selective unfolding with time constant ∼6.5 h, consistent with destabilization due to residue transfer to the polymerase′ domain on the p66′ subunit. A simultaneous increase in the intensity of resonances near the random coil positions is characterized by a similar time constant. Consistent with the residue transfer hypothesis, a construct of the isolated RH domain lacking the two N-terminal residues is shown to exhibit reduced stability. These results demonstrate that RH′ unfolding is coupled to homodimer formation.

Side chain conformational averaging in human dihydrofolate reductase

The three-dimensional structures of the dihydrofolate reductase enzymes from Escherichia coli (ecDHFR or ecE) and Homo sapiens (hDHFR or hE) are very similar, despite a rather low level of sequence identity. Whereas the active site loops of ecDHFR undergo major conformational rearrangements during progression through the reaction cycle, hDHFR remains fixed in a closed loop conformation in all of its catalytic intermediates. To elucidate the structural and dynamic differences between the human and E. coli enzymes, we conducted a comprehensive analysis of side chain flexibility and dynamics in complexes of hDHFR that represent intermediates in the major catalytic cycle. Nuclear magnetic resonance relaxation dispersion experiments show that, in marked contrast to the functionally important motions that feature prominently in the catalytic intermediates of ecDHFR, millisecond time scale fluctuations cannot be detected for hDHFR side chains. Ligand flux in hDHFR is thought to be mediated by conformational changes between a hinge-open state when the substrate/product-binding pocket is vacant and a hinge-closed state when this pocket is occupied. Comparison of X-ray structures of hinge-open and hinge-closed states shows that helix αF changes position by sliding between the two states. Analysis of χ₁ rotamer populations derived from measurements of ³J_CγCO and ³J_CγN couplings indicates that many of the side chains that contact helix αF exhibit rotamer averaging that may facilitate the conformational change. The χ₁ rotamer adopted by the Phe31 side chain depends upon whether the active site contains the substrate or product. In the holoenzyme (the binary complex of hDHFR with reduced nicotinamide adenine dinucleotide phosphate), a combination of hinge opening and a change in the Phe31 χ₁ rotamer opens the active site to facilitate entry of the substrate. Overall, the data suggest that, unlike ecDHFR, hDHFR requires minimal backbone conformational rearrangement as it proceeds through its enzymatic cycle, but that ligand flux is brokered by more subtle conformational changes that depend on the side chain motions of critical residues.

Ab Initio Calculations of Possible γ-Gauche Effects in the 13C-NMR for Methine and Carbonyl Carbons in Precise Polyethylene Acrylic Acid Copolymers

The impacts of local polymer chain conformations on the methine and carbonyl ¹³C-NMR chemical shifts for polyethylene acrylic acid p(E-AA) copolymers were predicted using ab initio methods. Using small molecular cluster models, the magnitude and sign of the γ-gauche torsional angle effect, along with the impact of local tetrahedral structure distortions near the carbonyl group, on the ¹³C-NMR chemical shifts were determined. These ¹³C-NMR chemical shift variations were compared to the experimental trends observed for precise p(E-AA) copolymers as a function acid group spacing and degree of zinc-neutralization in the corresponding p(E-AA) ionomers. These ab initio calculations address the future ability of ¹³C-NMR chemical shift variations to provide information about the local chain conformations in p(E-AA) copolymer materials.

Allosteric communication in the KIX domain proceeds through dynamic repacking of the hydrophobic core

The KIX domain of the transcriptional coactivator CREB binding protein (CBP) co-operatively mediates interactions between transcription factors. Binding of the transcription factor mixed-lineage leukemia (MLL) induces the formation of a low-populated conformer of KIX that resembles the conformation of the KIX domain in the presence of a second transcription factor molecule. NMR spin relaxation studies have previously shown that allosteric coupling proceeds through a network of hydrophobic core residues that bridge the two binding sites. Here we describe high-resolution NMR solution structures of the binary complex of KIX with MLL and the ternary complex of KIX formed with MLL and phosphorylated kinase inducible domain of CREB (pKID) as a second ligand. We show that binding of pKID to the binary complex of KIX with MLL is accompanied by a defined repacking of the allosteric network in the hydrophobic core of the protein. Rotamer populations derived from methyl group ¹³C chemical shifts reveal a dynamic contribution to the repacking process that is not captured by the structural coordinates and exemplify the dynamic nature of allosteric communication in the KIX domain.

Contribution of high-energy conformations to NMR chemical shifts, a DFT-BOMD study

This paper highlights the relevance of including the high-energy conformational states sampled by Born-Oppenheimer molecular dynamics (BOMD) in the calculation of time-averaged NMR chemical shifts. Our case study is the very flexible glycerol molecule that undergoes interconversion between conformers in a nonrandom way. Along the sequence of structures from one backbone conformer to another, transition states have been identified. The three (13)C NMR chemical shifts of the molecule were estimated by averaging their calculated values over a large set of BOMD snapshots. The simulation time needed to obtain a good agreement with the two signals present in the experimental spectrum is shown to be dependent on the atomic orbital basis set used for the dynamics, with a necessary longer trajectory for the most extended basis sets. The large structural deformations with respect to the optimized conformer geometries that occur along the dynamics are related to a kinetically driven conformer distribution. Calculated conformer type populations are in good agreement with experimental gas phase microwave results.

Efficacy of the β₂-adrenergic receptor is determined by conformational equilibrium in the transmembrane region

Many drugs that target G-protein-coupled receptors (GPCRs) induce or inhibit their signal transduction with different strengths, which affect their therapeutic properties. However, the mechanism underlying the differences in the signalling levels is still not clear, although several structures of GPCRs complexed with ligands determined by X-ray crystallography are available. Here we utilized NMR to monitor the signals from the methionine residue at position 82 in neutral antagonist- and partial agonist-bound states of β₂-adrenergic receptor (β₂AR), which are correlated with the conformational changes of the transmembrane regions upon activation. We show that this residue exists in a conformational equilibrium between the inverse agonist-bound states and the full agonist-bound state, and the population of the latter reflects the signal transduction level in each ligand-bound state. These findings provide insights into the multi-level signalling of β₂AR and other GPCRs, including the basal activity, and the mechanism of signal transduction mediated by GPCRs.

Many drugs exist that target the β-adrenergic receptor, but they have different efficacies. Kofuku et al. use NMR to show that methionine 82 in the transmembrane domain undergoes conformational changes depending on whether agonists or inverse agonists are bound, explaining the differential drug efficacy.

Solid-state NMR study reveals collagen I structural modifications of amino acid side chains upon fibrillogenesis

In vivo, collagen I, the major structural protein in human body, is found assembled into fibrils. In the present work, we study a high concentrated collagen sample in its soluble, fibrillar, and denatured states using one and two dimensional {(1)H}-(13)C solid-state NMR spectroscopy. We interpret (13)C chemical shift variations in terms of dihedral angle conformation changes. Our data show that fibrillogenesis increases the side chain and backbone structural complexity. Nevertheless, only three to five rotameric equilibria are found for each amino acid residue, indicating a relatively low structural heterogeneity of collagen upon fibrillogenesis. Using side chain statistical data, we calculate equilibrium constants for a great number of amino acid residues. Moreover, based on a (13)C quantitative spectrum, we estimate the percentage of residues implicated in each equilibrium. Our data indicate that fibril formation greatly affects hydroxyproline and proline prolyl pucker ring conformation. Finally, we discuss the implication of these structural data and propose a model in which the attractive force of fibrillogenesis comes from a structural reorganization of 10 to 15% of the amino acids. These results allow us to further understand the self-assembling process and fibrillar structure of collagen.

Improved chemical shift prediction by Rosetta conformational sampling

Chemical shift frequencies represent a time-average of all the conformational states populated by a protein. Thus, chemical shift prediction programs based on sequence and database analysis yield higher accuracy for rigid rather than flexible protein segments. Here we show that the prediction accuracy can be significantly improved by averaging over an ensemble of structures, predicted solely from amino acid sequence with the Rosetta program. This approach to chemical shift and structure prediction has the potential to be useful for guiding resonance assignments, especially in solid-state NMR structural studies of membrane proteins in proteoliposomes.

Metal and ligand binding to the HIV-RNase H active site are remotely monitored by Ile556

HIV-1 reverse transcriptase (RT) contains a C-terminal ribonuclease H (RH) domain on its p66 subunit that can be expressed as a stable, although inactive protein. Recent studies of several RH enzymes demonstrate that substrate binding plays a major role in the creation of the active site. In the absence of substrate, the C-terminal helix E of the RT RNase H domain is dynamic, characterized by severe exchange broadening of its backbone amide resonances, so that the solution characterization of this region of the protein has been limited. Nuclear magnetic resonance studies of 13C-labeled RH as a function of experimental conditions reveal that the δ1 methyl resonance of Ile556, located in a short, random coil segment following helix E, experiences a large 13C shift corresponding to a conformational change of Ile556 that results from packing of helix E against the central β-sheet. This shift provides a useful basis for monitoring the effects of various ligands on active site formation. Additionally, we report that the RNase H complexes formed with one or both divalent ions can be individually observed and characterized using diamagnetic Zn2+ as a substitute for Mg2+. Ordering of helix E results specifically from the interaction with the lower affinity binding to the A divalent ion site.

A strong 13C chemical shift signature provides the coordination mode of histidines in zinc-binding proteins

Zinc is the second most abundant metal ion incorporated in proteins, and is in many cases a crucial component of protein three-dimensional structures. Zinc ions are frequently coordinated by cysteine and histidine residues. Whereas cysteines bind to zinc via their unique S(γ) atom, histidines can coordinate zinc with two different coordination modes, either N(δ1) or N(ε2) is coordinating the zinc ion. The determination of this coordination mode is crucial for the accurate structure determination of a histidine-containing zinc-binding site by NMR. NMR chemical shifts contain a vast amount of information on local electronic and structural environments and surprisingly their utilization for the determination of the coordination mode of zinc-ligated histidines has been limited so far to (15)N nuclei. In the present report, we observed that the (13)C chemical shifts of aromatic carbons in zinc-ligated histidines represent a reliable signature of their coordination mode. Using a statistical analysis of (13)C chemical shifts, we show that (13)C(δ2) chemical shift is sensitive to the histidine coordination mode and that the chemical shift difference δ{(13)C(ε1)} - δ{(13)C(δ2)} provides a reference-independent marker of this coordination mode. The present approach allows the direct determination of the coordination mode of zinc-ligated histidines even with non-isotopically enriched protein samples and without any prior structural information.

Is a malleable protein necessarily highly dynamic? The hydrophobic core of the nuclear coactivator binding domain is well ordered

The nuclear coactivator binding domain of CREB binding protein folds into remarkably different structures in complex with different ligands. To understand the mechanism of the structural adaptability in the nuclear coactivator binding domain (NCBD), we have compared the dynamics of the hydrophobic core of NCBD in the ligand-free state and in a well-folded complex with the ligand activator for thyroid hormone and retinoid receptors using multiple NMR methods including methyl chemical shifts, coupling constants, and methyl order parameters. From all NMR measures, the aliphatic side chains in the hydrophobic core are slightly more dynamic in the free protein than in the complex, but have mobility comparable to the hydrophobic cores of average folded proteins. Urea titration monitored by NMR reveals that all parts of the protein, including the side-chain packing in the hydrophobic core, denatures in a single cooperative process. The molten globule characteristics of NCBD are thus restricted to a slowly fluctuating tertiary structure. Consequently, the conformational plasticity of the protein is most likely related to its low overall stability rather than an intrinsically flexible protein structure. The well-defined structure supports a model of molecular recognition dominated by conformational selection, whereas only minor structural adjustments are necessary after the association.

The interplay between transient α-helix formation and side chain rotamer distributions in disordered proteins probed by methyl chemical shifts

The peptide backbones of disordered proteins are routinely characterized by NMR with respect to transient structure and dynamics. Little experimental information is, however, available about the side chain conformations and how structure in the backbone affects the side chains. Methyl chemical shifts can in principle report the conformations of aliphatic side chains in disordered proteins and in order to examine this two model systems were chosen: the acid denatured state of acyl-CoA binding protein (ACBP) and the intrinsically disordered activation domain of the activator for thyroid hormone and retinoid receptors (ACTR). We find that small differences in the methyl carbon chemical shifts due to the γ-gauche effect may provide information about the side chain rotamer distributions. However, the effects of neighboring residues on the methyl group chemical shifts obscure the direct observation of γ-gauche effect. To overcome this, we reference the chemical shifts to those in a more disordered state resulting in residue specific random coil chemical shifts. The (13)C secondary chemical shifts of the methyl groups of valine, leucine, and isoleucine show sequence specific effects, which allow a quantitative analysis of the ensemble of χ(2)-angles of especially leucine residues in disordered proteins. The changes in the rotamer distributions upon denaturation correlate to the changes upon helix induction by the co-solvent trifluoroethanol, suggesting that the side chain conformers are directly or indirectly related to formation of transient α-helices.

Mechanisms of allosteric gene regulation by NMR quantification of microsecond-millisecond protein dynamics

The trp RNA-binding attenuation protein (TRAP) is a paradigmatic allosteric protein that regulates the tryptophan biosynthetic genes associated with the trp operon in bacilli. The ring-shaped 11-mer TRAP is activated for recognition of a specific trp-mRNA target by binding up to 11 tryptophan molecules. To characterize the mechanisms of tryptophan-induced TRAP activation, we have performed methyl relaxation dispersion (MRD) nuclear magnetic resonance (NMR) experiments that probe the time-dependent structure of TRAP in the microsecond-to-millisecond "chemical exchange" time window. We find significant side chain flexibility localized to the RNA and tryptophan binding sites of the apo protein and that these dynamics are dramatically reduced upon ligand binding. Analysis of the MRD NMR data provides insights into the structural nature of transiently populated conformations sampled in solution by apo TRAP. The MRD data are inconsistent with global two-state exchange, indicating that conformational sampling in apo TRAP is asynchronous. These findings imply a temporally heterogeneous population of structures that are incompatible with RNA binding and substantiate the study of TRAP as a paradigm for probing and understanding essential dynamics in allosteric, regulatory proteins.

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.