Protein Dynamics revealed by NMR Relaxation Methods

The path to the G protein-coupled receptor structural landscape: Major milestones and future directions

G protein-coupled receptors (GPCRs) play a crucial role in cell function by transducing signals from the extracellular environment to the inside of the cell. They mediate the effects of various stimuli, including hormones, neurotransmitters, ions, photons, food tastants and odorants, and are renowned drug targets. Advancements in structural biology techniques, including X-ray crystallography and cryo-electron microscopy (cryo-EM), have driven the elucidation of an increasing number of GPCR structures. These structures reveal novel features that shed light on receptor activation, dimerization and oligomerization, dichotomy between orthosteric and allosteric modulation, and the intricate interactions underlying signal transduction, providing insights into diverse ligand-binding modes and signalling pathways. However, a substantial portion of the GPCR repertoire and their activation states remain structurally unexplored. Future efforts should prioritize capturing the full structural diversity of GPCRs across multiple dimensions. To do so, the integration of structural biology with biophysical and computational techniques will be essential. We describe in this review the progress of nuclear magnetic resonance (NMR) to examine GPCR plasticity and conformational dynamics, of atomic force microscopy (AFM) to explore the spatial-temporal dynamics and kinetic aspects of GPCRs, and the recent breakthroughs in artificial intelligence for protein structure prediction to characterize the structures of the entire GPCRome. In summary, the journey through GPCR structural biology provided in this review illustrates how far we have come in decoding these essential proteins architecture and function. Looking ahead, integrating cutting-edge biophysics and computational tools offers a path to navigating the GPCR structural landscape, ultimately advancing GPCR-based applications.

Parsing Dynamics of Protein Backbone NH and Side-Chain Methyl Groups using Molecular Dynamics Simulations

Experimental NMR spectroscopy and theoretical molecular dynamics (MD) simulations provide complementary insights into protein conformational dynamics and hence into biological function. The present work describes an extensive set of backbone NH and side-chain methyl group generalized order parameters for the Escherichia coli ribonuclease HI (RNH) enzyme derived from 2-μs microsecond MD simulations using the OPLS4 and AMBER-FF19SB force fields. The simulated generalized order parameters are compared with values derived from NMR 15N and 13CH2D spin relaxation measurements. The squares of the generalized order parameters, S2 for the N-H bond vector and Saxis2 for the methyl group symmetry axis, characterize the equilibrium distribution of vector orientations in a molecular frame of reference. Optimal agreement between simulated and experimental results was obtained by averaging S2 or Saxis2 calculated by dividing the simulated trajectories into 50 ns blocks (∼five times the rotational diffusion correlation time for RNH). With this procedure, the median absolute deviations (MAD) between experimental and simulated values of S2 and Saxis2 are 0.030 (NH) and 0.061 (CH3) for OPLS4 and 0.041 (NH) and 0.078 (CH3) for AMBER-FF19SB. The MAD between OPLS4 and AMBER-FF19SB are 0.021 (NH) and 0.072 (CH3). The generalized order parameters for the methyl group symmetry axis can be decomposed into contributions from backbone fluctuations, between-rotamer dihedral angle transitions, and within-rotamer dihedral angle fluctuations. Analysis of the simulation trajectories shows that (i) backbone and side chain conformational fluctuations exhibit little correlation and that (ii) fluctuations within rotamers are limited and highly uniform with values that depend on the number of dihedral angles considered. Low values of Saxis2, indicative of enhanced side-chain flexibility, result from between-rotamer transitions that can be enhanced by increased local backbone flexibility.

Uncovering allostery and regulation in SORCIN through molecular dynamics simulations

Soluble resistance-related calcium-binding protein or Sorcin is an allosteric, calcium-binding Penta-EF hand (PEF) family protein implicated in multi-drug resistant cancers. Sorcin is known to bind chemotherapeutic molecules such as Doxorubicin. This study uses in-silico molecular dynamics simulations to explore the dynamics and allosteric behavior of Sorcin in the context of Ca2+ uptake and Doxorubicin binding. The results show that Ca2+ binding induces large, but reversible conformational changes in the Sorcin structure which manifest as rigid body reorientations that preserve the local secondary structure. A reciprocal allosteric handshake centered around the EF5 hand is found to be key in Sorcin dimer formation and stabilization. Binding of Doxorubicin results in rearrangement of allosteric communities which disrupts long-range allosteric information transfer from the N-terminal domain to the middle lobe. However, this binding does not result in secondary structure destabilization. Sorcin does not appear to have a distinct Ca2+ activated mode of Doxorubicin binding.Communicated by Ramaswamy H. Sarma.

Combining NMR, SAXS and SANS to characterize the structure and dynamics of protein complexes

Understanding the structure and dynamics of biological macromolecules is essential to decipher the molecular mechanisms that underlie cellular functions. The description of structure and conformational dynamics often requires the integration of complementary techniques. In this review, we highlight the utility of combining nuclear magnetic resonance (NMR) spectroscopy with small angle scattering (SAS) to characterize these challenging biomolecular systems. NMR can assess the structure and conformational dynamics of multidomain proteins, RNAs and biomolecular complexes. It can efficiently provide information on interaction surfaces, long-distance restraints and relative domain orientations at residue-level resolution. Such information can be readily combined with high-resolution structural data available on subcomponents of biomolecular assemblies. Moreover, NMR is a powerful tool to characterize the dynamics of biomolecules on a wide range of timescales, from nanoseconds to seconds. On the other hand, SAS approaches provide global information on the size and shape of biomolecules and on the ensemble of all conformations present in solution. Therefore, NMR and SAS provide complementary data that are uniquely suited to investigate dynamic biomolecular assemblies. Here, we briefly review the type of data that can be obtained by both techniques and describe different approaches that can be used to combine them to characterize biomolecular assemblies. We then provide guidelines on which experiments are best suited depending on the type of system studied, ranging from fully rigid complexes, dynamic structures that interconvert between defined conformations and systems with very high structural heterogeneity.

The transient expression of recombinant proteins in plant cell packs facilitates stable isotope labelling for NMR spectroscopy

Nuclear magnetic resonance (NMR) spectroscopy can be used to determine the structure, dynamics and interactions of proteins. However, protein NMR requires stable isotope labelling for signal detection. The cells used for the production of recombinant proteins must therefore be grown in medium containing isotopically labelled substrates. Stable isotope labelling is well established in Escherichia coli, but bacteria are only suitable for the production of simple proteins without post-translational modifications. More complex proteins require eukaryotic production hosts, but their growth can be impaired by labelled media, thus reducing product yields and increasing costs. To address this limitation, we used media supplemented with isotope-labelled substrates to cultivate the tobacco-derived cell line BY-2, which was then cast into plant cell packs (PCPs) for the transient expression of a labelled version of the model protein GB1. Mass spectrometry confirmed the feasibility of isotope labelling with ¹⁵ N and ² H using this approach. The resulting NMR spectrum featured a signal dispersion comparable to recombinant GB1 produced in E. coli. PCPs therefore offer a rapid and cost-efficient alternative for the production of isotope-labelled proteins for NMR analysis, especially suitable for complex proteins that cannot be produced in microbial systems.

Effect of Ligands and Transducers on the Neurotensin Receptor 1 Conformational Ensemble

Using a discrete, intracellular ¹⁹F nuclear magnetic resonance (NMR) probe on transmembrane helix 6 of the neurotensin receptor 1 (NTS1), we aim to understand how ligands and transducers modulate the receptor's structural ensemble in a solution. For apo NTS1, ¹⁹F NMR spectra reveal an ensemble of at least three conformational substates (one inactive and two active-like) in equilibrium that exchange on the millisecond to second timescale. Dynamic NMR experiments reveal that these substates follow a linear three-site exchange process that is both thermodynamically and kinetically remodeled by orthosteric ligands. As previously observed in other G protein-coupled receptors (GPCRs), the full agonist is insufficient to completely stabilize the active-like state. The inactive substate is abolished upon coupling to β-arrestin-1 (βArr1) or the C-terminal helix of Gα_q, which comprises ≳60% of the GPCR/G protein interface surface area. Whereas βArr1 exclusively selects for pre-existing active-like substates, the Gα_q peptide induces a new substate. Both transducer molecules promote substantial line broadening of active-like states, suggesting contributions from additional microsecond to millisecond exchange processes. Together, our study suggests that (i) the NTS1 allosteric activation mechanism may be alternatively dominated by induced fit or conformational selection depending on the coupled transducer, and (ii) the available static structures do not represent the entire conformational ensemble observed in a solution.

Insights into the Cross Talk between Effector and Allosteric Lobes of KRAS from Methyl Conformational Dynamics

KRAS is the most frequently mutated RAS protein in cancer patients, and it is estimated that about 20% of the cancer patients in the United States carried mutant RAS proteins. To accelerate therapeutic development, structures and dynamics of RAS proteins had been extensively studied by various biophysical techniques for decades. Although 31P NMR studies revealed population equilibrium of the two major states in the active GMPPNP-bound form, more complex conformational dynamics in RAS proteins and oncogenic mutants subtly modulate the interactions with their downstream effectors. We established a set of customized NMR relaxation dispersion techniques to efficiently and systematically examine the ms-μs conformational dynamics of RAS proteins. This method allowed us to observe varying synchronized motions that connect the effector and allosteric lobes in KRAS. We demonstrated the role of conformational dynamics of KRAS in controlling its interaction with the Ras-binding domain of the downstream effector RAF1, the first kinase in the MAPK pathway. This allows one to explain, as well as to predict, the altered binding affinities of various KRAS mutants, which was neither previously reported nor apparent from the structural perspective.

Subsecond Time-Resolved Mass Spectrometry in Dynamic Structural Biology

Life at the molecular level is a dynamic world, where the key players-proteins, oligonucleotides, lipids, and carbohydrates-are in a perpetual state of structural flux, shifting rapidly between local minima on their conformational free energy landscapes. The techniques of classical structural biology, X-ray crystallography, structural NMR, and cryo-electron microscopy (cryo-EM), while capable of extraordinary structural resolution, are innately ill-suited to characterize biomolecules in their dynamically active states. Subsecond time-resolved mass spectrometry (MS) provides a unique window into the dynamic world of biological macromolecules, offering the capacity to directly monitor biochemical processes and conformational shifts with a structural dimension provided by the electrospray charge-state distribution, ion mobility, covalent labeling, or hydrogen-deuterium exchange. Over the past two decades, this suite of techniques has provided important insights into the inherently dynamic processes that drive function and pathogenesis in biological macromolecules, including (mis)folding, complexation, aggregation, ligand binding, and enzyme catalysis, among others. This Review provides a comprehensive account of subsecond time-resolved MS and the advances it has enabled in dynamic structural biology, with an emphasis on insights into the dynamic drivers of protein function.

Theoretical classification of exchange geometries from the perspective of NMR relaxation dispersion

NMR relaxation dispersion experiments have been widely applied to probe important conformational exchange of macro-molecules in many biological systems. The current improvements in computational techniques as well as the theoretical breakthroughs make the quantitative data analysis of complex exchange models possible. However, the topology of a given exchange model is also one of the main factors affecting the solution of Bloch-McConnell equation. The lack of a theoretical analysis of the exchange topologies at n-site exchange hinders further progress of such data analysis. Here, using graph theory, we reveal the topological complexity of n-site exchange and present all exchange models when n is less than 6. Furthermore, we introduce an alternative way, using machine learning, to select an exchange model based on a set of relaxation dispersion data without fitting them with every individual exchange model.

Structural Dynamics from NMR Relaxation by SRLS Analysis: Local Geometry, Potential Energy Landscapes, and Spectral Densities

We have developed the two-body coupled-rotator slowly relaxing local structure (SRLS) approach for elucidating protein dynamics by nuclear magnetic resonance (NMR) relaxation. The rotators are represented by diffusion tensors D₁ for overall protein tumbling and D₂ for locally ordered probe motion. D₁ and D₂ are coupled dynamically by a potential, u, typically given by linear combinations of the Wigner functions D₀₀² and (D₀₂² + D_0-2²). Until now, our SRLS analyses provided the tensors, D₁ and D₂, the potential, u, and the geometric link between SRLS and NMR. Here we enhance this description by also examining the SRLS spectral densities obtained by solving the SRLS Smoluchowski equation. In addition, we show that the form of u specified above complies with two NMR-detected potential energy landscapes representing preferential ordering along N-H or C^α-C^α. Pictorial illustrations thereof are provided. The extended SRLS analysis is applied to ¹⁵N-H relaxation from the carbohydrate recognition domain of galectin-3 (Gal3C) in complex with two diastereomeric ligands, S and R. We find that D₂ is isotropic with a principal value, D₂, of 10¹⁰ s^-1 on average, and it is faster in the strands β₃, β₅, and β₈. The potential, u, is strong (∼20 kT); it is slightly rhombic when N-H is the main ordering axis and highly rhombic when C^α-C^α is the main ordering axis. Gal3C-S exhibits primarily preferential ordering along C^α-C^α; Gal3C-R exhibits both types of ordering. The binding-associated polypeptide chain segment of Gal3C-S is homogeneous, whereas that of Gal3C-R is diversified, with regard to D₂ and ordering preference. We associate these features with the previously determined diminished binding constant of Gal3C-R in comparison with Gal3C-S. Thus, the present study enhances the SRLS analysis, in general, and provides new insights into the dynamic structure and binding properties of Gal3C-S and Gal3C-R, in particular.

Reductionism in the study of enzyme adaptation

One of the principal goals of comparative biology is the elucidation of mechanisms by which organisms adapt to different environments. The study of enzyme structure, function, and stability has contributed significantly to this effort, by revealing adaptation at a molecular level. Comparative biochemistry, including enzymology, necessarily pursues a reductionist approach in describing the function and structure of biomolecules, allowing more straightforward study of molecular systems by removing much of the complexity of their biological milieu. Although this reductionism has allowed a remarkable series of discoveries linking chemical processes to metabolism and to whole-organism function in the context of the environment, it also has the potential to mislead when careful consideration is not made of the simplifying assumptions inherent to such research. In this review, a brief history of the growth of enzymology, its reliance on a reductionist philosophy, and its contributions to our understanding of biological systems is given. Examples then are provided of research techniques, based on a reductionist approach, that have advanced our knowledge about enzyme adaptation to environmental stresses, including stability assays, enzyme kinetics, and the impact of solute composition on enzyme function. In each case, the benefits of the reductionist nature of the approach is emphasized, notable advances are described, but potential drawbacks due to inherent oversimplification of the study system are also identified.

The Role of Conformational Dynamics in the Recognition and Regulation of Ubiquitination

The ubiquitination pathway is central to many cell signaling and regulatory events. One of the intriguing aspects of the pathway is the combinatorial sophistication of substrate recognition and ubiquitin chain building determinations. The abundant structural and biological data portray several characteristic protein folds among E2 and E3 proteins, and the understanding of the combinatorial complexity that enables interaction with much of the human proteome is a major goal to developing targeted and selective manipulation of the pathway. With the commonality of some folds, there are likely other aspects that can provide differentiation and recognition. These aspects involve allosteric effects and conformational dynamics that can direct recognition and chain building processes. In this review, we will describe the current state of the knowledge for conformational dynamics across a wide timescale, address the limitations of present approaches, and illustrate the potential to make new advances in connecting dynamics with ubiquitination regulation.

Unravelling the Time Scale of Conformational Plasticity and Allostery in Glycan Recognition by Human Galectin-1

The interaction of human galectin-1 with a variety of oligosaccharides, from di-(N-acetyllactosamine) to tetra-saccharides (blood B type-II antigen) has been scrutinized by using a combined approach of different NMR experiments, molecular dynamics (MD) simulations, and isothermal titration calorimetry. Ligand- and receptor-based NMR experiments assisted by computational methods allowed proposing three-dimensional structures for the different complexes, which explained the lack of enthalpy gain when increasing the chemical complexity of the glycan. Interestingly, and independently of the glycan ligand, the entropy term does not oppose the binding event, a rather unusual feature for protein-sugar interactions. CLEANEX-PM and relaxation dispersion experiments revealed that sugar binding affected residues far from the binding site and described significant changes in the dynamics of the protein. In particular, motions in the microsecond-millisecond timescale in residues at the protein dimer interface were identified in the presence of high affinity ligands. The dynamic process was further explored by extensive MD simulations, which provided additional support for the existence of allostery in glycan recognition by human galectin-1.

Probing the Broad Time Scale and Heterogeneous Conformational Dynamics in the Catalytic Core of the Arf-GAP ASAP1 via Methyl Adiabatic Relaxation Dispersion

Methyl-TROSY is one of the most powerful NMR spectroscopic tools for studying structures and conformational dynamics of large protein complexes in solution. In studying conformational dynamics, side chains usually display heterogeneous dynamics, including collective and local motions, that can be difficult to detect and analyze by conventional relaxation dispersion (RD) approaches. The combination of NH-based heteronuclear adiabatic relaxation dispersion (HARD) experiments and a geometric approximation (geoHARD) has been shown to have several advantages over conventional RD in revealing conformational dynamics over a broad time scale. Here, we demonstrate a new technique that has been developed to detect both heterogeneous and wide time scale conformational dynamics in the hydrophobic interior of large macromolecules utilizing methyl-geoHARD. It is shown that methyl-geoHARD will be feasible at ultrahigh magnetic fields (>1 GHz), when this technology becomes available. For the ZA domain of Arf-GAP ASAP1, with a global correlational time of 24 ns at 15 °C, a wide range of conformational dynamics (exhibiting chemical exchange rates (kex) between 102 and 105 s-1) are observed in the methyl groups of isoleucine, leucine, and valine. The dynamics include collective and independent local motions. Furthermore, portions of the collective motions have been confirmed by single-quantum Carr-Purcell-Meiboom-Gill (SQ-CPMG) RD experiments; however, motions outside of the detectable CPMG window (400-8000 s-1) cannot be accurately determined by SQ-CPMG experiments. The methyl-geoHARD experiment allows the dissection of heterogeneous conformational dynamics and pinpoints important motions that, potentially, can be correlated with important biological functions and recognition.