Defining the Structural Basis for Allosteric Product Release from E. coli Dihydrofolate Reductase Using NMR Relaxation Dispersion

Molecular dynamics simulation on the role of CL5D in accelerating the product dissociation of SIRT6

SIRT6 is a member of the NAD+-dependent histone deacetylase family and is integral to maintaining genome stability and regulating metabolic transcription. SIRT6 transfers acetyl groups from the lysine side chains of protein substrates to the cofactor NAD+, generating nicotinamide, 2'-O-acyl-ADP-ribose (ADPr), and a deacetylated substrate. SIRT6 has been found to be activated by small molecule activators, such as CL5D. However, the process of dissociation of the SIRT6 product and the mechanism of activation by small molecule activators are unknown. In this work, we elucidated these activation mechanisms by performing extensive molecular dynamics simulations. The results of random acceleration molecular dynamics and umbrella sampling demonstrated that the dissociation sequence involves the exit of the deacetylated substrate first, followed by ADPr. The binding of CL5D does not alter the dissociation pathway of the products, but it increases the catalytic activity of SIRT6 by facilitating the dissociation of products within SIRT6. Our results suggest a mechanism of SIRT6 activation, which highlights the importance of product dissociation in enzyme catalysis. This result may help facilitate the development of new SIRT6 activators.

Hook loop dynamics engineering transcended the barrier of activity-stability trade-off and boosted the thermostability of enzymes

The improvement of enzyme thermostability often accompanies the decreased activity due to the loss of the key regions' flexibility. As a representative structure, unlocking the potential of loop dynamics will not only provide new ideas for stabilization strategies, but also help to deepen the understanding of the relationship between enzyme structural dynamics and function. In this study, a creative "hook loop dynamics engineering" (HLoD) strategy was successfully proposed for simultaneously improving the thermostability and maintaining activity of the model enzyme, Candida Antarctica lipase B. A small and smart mutant library involving five key residues located at the "hook loop" was meticulously identified and systematically investigated and thus yielded a five-point multiple mutant M1 (L147S/T244P/S250P/T256D/N292D), demonstrating a remarkable 7.0-fold increase in thermostability at 60 °C compared to the wild-type (WT). Furthermore, the activity of M1 remained comparable to that of WT, effectively transcending the barrier of activity-stability trade-off. Molecular dynamics simulations revealed that the precise regulation of hook loop dynamics via intermolecular interactions, such as salt bridges and hydrogen bonding, curbed the excessive flexibility of the pivotal regions α5 and α10 at high temperatures, thus driving the substantial enhancement of the thermostability of M1. Refining the dynamics of the flexible region via HLoD, which transcended the barrier of activity-stability trade-off, exhibited to be a robust and potentially universal strategy for designing enzymes with outstanding thermostability and activity.

Connecting Protein Millisecond Conformational Dynamics to Protein Thermal Stability

The stability of protein folded states is crucial for its function, yet the relationship with the protein sequence remains poorly understood. Prior studies have focused on the amino acid composition and thermodynamic couplings within a single folded conformation, overlooking the potential contribution of protein dynamics. Here, we address this gap by systematically analyzing the impact of alanine mutations in the C-terminal β-strand (β5) of ubiquitin, a model protein exhibiting millisecond timescale interconversion between two conformational states differing in the β5 position. Our findings unveil a negative correlation between millisecond dynamics and thermal stability, with alanine substitutions at seemingly flexible C-terminal residues significantly enhancing thermostability. Integrating spectroscopic and computational approaches, we demonstrate that the thermally unfolded state retains a substantial secondary structure but lacks β5 engagement, recapitulating the transition state for millisecond dynamics. Thus, alanine mutations that modulate the stabilities of the folded states with respect to the partially unfolded state impact both the dynamics and stability. Our findings underscore the importance of conformational dynamics with implications for protein engineering and design.

What is allosteric regulation? Exploring the exceptions that prove the rule!

"Allosteric" was first introduced to mean the other site (i.e., a site distinct from the active or orthosteric site), an adjective for "regulation" to imply a regulatory outcome resulting from ligand binding at another site. That original idea outlines a system with two ligand-binding events at two distinct locations on a macromolecule (originally a protein system), which defines a four-state energy cycle. An allosteric energy cycle provides a quantifiable allosteric coupling constant and focuses our attention on the unique properties of the four equilibrated protein complexes that constitute the energy cycle. Because many observed phenomena have been referenced as "allosteric regulation" in the literature, the goal of this work is to use literature examples to explore which systems are and are not consistent with the two-ligand thermodynamic energy cycle-based definition of allosteric regulation. We emphasize the need for consistent language so comparisons can be made among the ever-increasing number of allosteric systems. Building on the mutually exclusive natures of an energy cycle definition of allosteric regulation versus classic two-state models, we conclude our discussion by outlining how the often-proposed Rube-Goldberg-like mechanisms are likely inconsistent with an energy cycle definition of allosteric regulation.

Insights into the Structure of Invisible Conformations of Large Methyl Group Labeled Molecular Machines from High Pressure NMR

Most proteins are highly flexible and can adopt conformations that deviate from the energetically most favorable ground state. Structural information on these lowly populated, alternative conformations is often lacking, despite the functional importance of these states. Here, we study the pathway by which the Dcp1:Dcp2 mRNA decapping complex exchanges between an autoinhibited closed and an open conformation. We make use of methyl Carr-Purcell-Meiboom-Gill (CPMG) NMR relaxation dispersion (RD) experiments that report on the population of the sparsely populated open conformation as well as on the exchange rate between the two conformations. To obtain volumetric information on the open conformation as well as on the transition state structure we made use of RD measurements at elevated pressures. We found that the open Dcp1:Dcp2 conformation has a lower molecular volume than the closed conformation and that the transition state is close in volume to the closed state. In the presence of ATP the volume change upon opening of the complex increases and the volume of the transition state lies in-between the volumes of the closed and open state. These findings show that ATP has an effect on the volume changes that are associated with the opening-closing pathway of the complex. Our results highlight the strength of pressure dependent NMR methods to obtain insights into structural features of protein conformations that are not directly observable. As our work makes use of methyl groups as NMR probes we conclude that the applied methodology is also applicable to high molecular weight complexes.

DHFR Mutants Modulate Their Synchronized Dynamics with the Substrate by Shifting Hydrogen Bond Occupancies

Antibiotic resistance is a global health problem in which mutations occurring in functional proteins render drugs ineffective. The working mechanisms of the arising mutants are seldom apparent; a methodology to decipher these mechanisms systematically would render devising therapies to control the arising mutational pathways possible. Here we utilize C_α-C_β bond vector relaxations obtained from moderate length MD trajectories to determine conduits for functionality of the resistance conferring mutants of Escherichia coli dihydrofolate reductase. We find that the whole enzyme is synchronized to the motions of the substrate, irrespective of the mutation introducing gain-of-function or loss-of function. The total coordination of the motions suggests changes in the hydrogen bond dynamics with respect to the wild type as a possible route to determine and classify the mode-of-action of individual mutants. As a result, nine trimethoprim-resistant point mutations arising frequently in evolution experiments are categorized. One group of mutants that display the largest occurrence (L28R, W30G) work directly by modifying the dihydrofolate binding region. Conversely, W30R works indirectly by the formation of the E139-R30 salt bridge which releases energy resulting from tight binding by distorting the binding cavity. A third group (D27E, F153S, I94L) arising as single, resistance invoking mutants in evolution experiment trajectories allosterically and dynamically affects a hydrogen bonding motif formed at residues 59-69-71 which in turn modifies the binding site dynamics. The final group (I5F, A26T, R98P) consists of those mutants that have properties most similar to the wild type; these only appear after one of the other mutants is fixed on the protein structure and therefore display clear epistasis. Thus, we show that the binding event is governed by the entire enzyme dynamics while the binding site residues play gating roles. The adjustments made in the total enzyme in response to point mutations are what make quantifying and pinpointing their effect a hard problem. Here, we show that hydrogen bond dynamics recorded on sub-μs time scales provide the necessary fingerprints to decipher the various mechanisms at play.

Conserved conformational dynamics determine enzyme activity

Homologous enzymes often exhibit different catalytic rates despite a fully conserved active site. The canonical view is that an enzyme sequence defines its structure and function and, more recently, that intrinsic protein dynamics at different time scales enable and/or promote catalytic activity. Here, we show that, using the protein tyrosine phosphatase PTP1B, residues surrounding the PTP1B active site promote dynamically coordinated chemistry necessary for PTP1B function. However, residues distant to the active site also undergo distinct intermediate time scale dynamics and these dynamics are correlated with its catalytic activity and thus allow for different catalytic rates in this enzyme family. We identify these previously undetected motions using coevolutionary coupling analysis and nuclear magnetic resonance spectroscopy. Our findings strongly indicate that conserved dynamics drives the enzymatic activity of the PTP family. Characterization of these conserved dynamics allows for the identification of novel regulatory elements (therapeutic binding pockets) that can be leveraged for the control of enzymes.

Single-Molecule Sampling of Dihydrofolate Reductase Shows Kinetic Pauses and an Endosteric Effect Linked to Catalysis

The ability to sample multiple reactions on the same single enzyme is important to link rare intermediates with catalysis and to unravel the role of conformational changes. Despite decades of efforts, however, the single-molecule characterization of nonfluorogenic enzymes during multiple catalytic turnovers has been elusive. Here, we show that nanopore currents allow sampling the dynamic exchange between five structural intermediates during E. coli dihydrofolate reductase (DHFR) catalysis. We found that an endosteric effect promotes the binding of the substrate to the enzyme with a specific hierarchy. The chemical step then switched the enzyme from the closed to the occluded conformation, which in turn promotes the release of the reduced cofactor NADP⁺. Unexpectedly, only a few reactive complexes lead to catalysis. Furthermore, second-long catalytic pauses were observed, possibly reflecting an off-path conformation generated during the reaction. Finally, the free energy from multiple cofactor binding events were required to release the product and switch DHFR back to the reactive conformer. This catalytic fueled concerted mechanism is likely to have evolved to improve the catalytic efficiency of DHFR under the high concentrations of NADP⁺ in E. coli and might be a general feature for complex enzymatic reactions where the binding and release of the products must be tightly controlled.

Enhanced specificity mutations perturb allosteric signaling in CRISPR-Cas9

CRISPR-Cas9 (clustered regularly interspaced short palindromic repeat and associated Cas9 protein) is a molecular tool with transformative genome editing capabilities. At the molecular level, an intricate allosteric signaling is critical for DNA cleavage, but its role in the specificity enhancement of the Cas9 endonuclease is poorly understood. Here, multi-microsecond molecular dynamics is combined with solution NMR and graph theory-derived models to probe the allosteric role of key specificity-enhancing mutations. We show that mutations responsible for increasing the specificity of Cas9 alter the allosteric structure of the catalytic HNH domain, impacting the signal transmission from the DNA recognition region to the catalytic sites for cleavage. Specifically, the K855A mutation strongly disrupts the allosteric connectivity of the HNH domain, exerting the highest perturbation on the signaling transfer, while K810A and K848A result in more moderate effects on the allosteric communication. This differential perturbation of the allosteric signal correlates to the order of specificity enhancement (K855A > K848A ~ K810A) observed in biochemical studies, with the mutation achieving the highest specificity most strongly perturbing the signaling transfer. These findings suggest that alterations of the allosteric communication from DNA recognition to cleavage are critical to increasing the specificity of Cas9 and that allosteric hotspots can be targeted through mutational studies for improving the system's function.

Role of Active Site Loop Dynamics in Mediating Ligand Release from E. coli Dihydrofolate Reductase

Conformational fluctuations from ground-state to sparsely populated but functionally important excited states play a key role in enzyme catalysis. For Escherichia coli dihydrofolate reductase (DHFR), the release of the product tetrahydrofolate (THF) and oxidized cofactor NADP⁺ occurs through exchange between closed and occluded conformations of the Met20 loop. A "dynamic knockout" mutant of E. coli DHFR, where the E. coli sequence in the Met20 loop is replaced by the human sequence (N23PP/S148A), models human DHFR and is incapable of accessing the occluded conformation. ¹H and ¹⁵N CPMG relaxation dispersion analysis for the ternary product complex of the mutant enzyme with NADP⁺ and the product analogue 5,10-dideazatetrahydrofolate (ddTHF) (E:ddTHF:NADP⁺) reveals the mechanism by which NADP⁺ is released when the Met20 loop cannot undergo the closed-to-occluded conformational transition. Two excited states were observed: one related to a faster, relatively high-amplitude conformational fluctuation in areas near the active site, associated with the shuttling of the nicotinamide ring of the cofactor out of the active site, and the other to a slower process where ddTHF undergoes small-amplitude motions within the binding site that are consistent with disorder observed in a room-temperature X-ray crystal structure of the N23PP/S148A mutant protein. These motions likely arise due to steric conflict of the pterin ring of ddTHF with the ribose-nicotinamide moiety of NADP⁺ in the closed active site. These studies demonstrate that site-specific kinetic information from relaxation dispersion experiments can provide intimate details of the changes in catalytic mechanism that result from small changes in local amino acid sequence.

Two β-Lactamase Variants with Reduced Clavulanic Acid Inhibition Display Different Millisecond Dynamics

The β-lactamase of Mycobacterium tuberculosis, BlaC, is susceptible to inhibition by clavulanic acid. The ability of this enzyme to escape inhibition through mutation was probed using error-prone PCR combined with functional screening in Escherichia coli. The variant that was found to confer the most inhibitor resistance, K234R, as well as variant G132N that was found previously were characterized using X-ray crystallography and nuclear magnetic resonance (NMR) relaxation experiments to probe structural and dynamic properties. The G132N mutant exists in solution in two almost equally populated conformations that exchange with a rate of ca. 88 s-1. The conformational change affects a broad region of the enzyme. The crystal structure reveals that the Asn132 side chain forces the peptide bond between Ser104 and Ile105 in a cis-conformation. The crystal structure suggests multiple conformations for several side chains (e.g., Ser104 and Ser130) and a short loop (positions 214 to 216). In the K234R mutant, the active-site dynamics are significantly diminished with respect to the wild-type enzyme. These results show that multiple evolutionary routes are available to increase inhibitor resistance in BlaC and that active-site dynamics on the millisecond time scale are not required for catalytic function.

Altered expression of a quality control protease in E. coli reshapes the in vivo mutational landscape of a model enzyme

Protein mutational landscapes are shaped by the cellular environment, but key factors and their quantitative effects are often unknown. Here we show that Lon, a quality control protease naturally absent in common E. coli expression strains, drastically reshapes the mutational landscape of the metabolic enzyme dihydrofolate reductase (DHFR). Selection under conditions that resolve highly active mutants reveals that 23.3% of all single point mutations in DHFR are advantageous in the absence of Lon, but advantageous mutations are largely suppressed when Lon is reintroduced. Protein stability measurements demonstrate extensive activity-stability tradeoffs for the advantageous mutants and provide a mechanistic explanation for Lon’s widespread impact. Our findings suggest possibilities for tuning mutational landscapes by modulating the cellular environment, with implications for protein design and combatting antibiotic resistance.

Automated assignment of methyl NMR spectra from large proteins

As structural biology trends towards larger and more complex biomolecular targets, a detailed understanding of their interactions and underlying structures and dynamics is required. The development of methyl-TROSY has enabled NMR spectroscopy to provide atomic-resolution insight into the mechanisms of large molecular assemblies in solution. However, the applicability of methyl-TROSY has been hindered by the laborious and time-consuming resonance assignment process, typically performed with domain fragmentation, site-directed mutagenesis, and analysis of NOE data in the context of a crystal structure. In response, several structure-based automatic methyl assignment strategies have been developed over the past decade. Here, we present a comprehensive analysis of all available methods and compare their input data requirements, algorithmic strategies, and reported performance. In general, the methods fall into two categories: those that primarily rely on inter-methyl NOEs, and those that utilize methyl PRE- and PCS-based restraints. We discuss their advantages and limitations, and highlight the potential benefits from standardizing and combining different methods.

Protein Dynamics revealed by NMR Relaxation Methods

Structural biology often focuses primarily on three-dimensional structures of biological macromolecules, deposited in the Protein Data Bank (PDB). This resource is a remarkable entity for the world-wide scientific and medical communities, as well as the general public, as it is a growing translation into three-dimensional space of the vast information in genomic databases, e.g. GENBANK. There is, however, significantly more to understanding biological function than the three-dimensional coordinate space for ground-state structures of biomolecules. The vast array of biomolecules experiences natural dynamics, interconversion between multiple conformational states, and molecular recognition and allosteric events that play out on timescales ranging from picoseconds to seconds. This wide range of timescales demands ingenious and sophisticated experimental tools to sample and interpret these motions, thus enabling clearer insight into functional annotation of the PDB. NMR spectroscopy is unique in its ability to sample this range of timescales at atomic resolution and in physiologically relevant conditions using spin relaxation methods. The field is constantly expanding to provide new creative experiments, to yield more detailed coverage of timescales, and to broaden the power of interpretation and analysis methods. This review highlights the current state of the methodology and examines the extension of analysis tools for more complex experiments and dynamic models. The future for understanding protein dynamics is bright, and these extended tools bring greater compatibility with developments in computational molecular dynamics, all of which will further our understanding of biological molecular functions. These facets place NMR as a key component in integrated structural biology.

Sensitivity of Functional Loop Conformations on Long-Range Electrostatics: Implications for M20 Loop Dynamics in E. coli Dihydrofolate Reductase

In E. coli dihydrofolate reductase, unusual conformational motions of a functional M20 loop that interacts with substrate and coenzyme have been construed as evidence for dynamical effects in enzyme catalysis. By computing this loop's conformational free energies in the apoenzyme, we show that it is sensitive to the treatment of long-range electrostatic interactions and the solvation box size in modeling/simulations. These results provide important guidelines for computing reaction/binding free energy profiles of proteins with functional loops.

A Dynamic Switch in Inactive p38γ Leads to an Excited State on the Pathway to an Active Kinase

The inactive state of mitogen-activated protein kinases (MAPKs) adopts an open conformation while the active state exists in a compact form stabilized by phosphorylation. In the active state, eukaryotic kinases undergo breathing motions related to substrate binding and product release that have not previously been detected in the inactive state. However, docking interactions of partner proteins with inactive MAPK kinases exhibit allostery in binding of activating kinases. Interactions at a site distant from the activation loop are coupled to the configuration of the activation loop, suggesting that the inactive state may also undergo concerted dynamics. X-ray crystallographic studies of nonphosphorylated, inactive p38γ reveal differences in domain orientations and active site structure in the two molecules in the asymmetric unit. One molecule resembles an inactive kinase with an open active site. The second molecule has a rotation of the N-lobe that leads to partial compaction of the active site, resulting in a conformation that is intermediate between the inactive open state and the fully closed state of the activated kinase. Although the compact state of apo p38γ displays several of the features of the activated enzyme, it remains catalytically inert. In solution, the kinase fluctuates on a millisecond time scale between the open ground state and a weakly populated excited state that is similar in structure to the compact state observed in the crystal. The nuclear magnetic resonance and crystal structure data imply that interconversion between the open and compact states involves a molecular switch associated with the DFG loop.

Branch site bulge conformations in domain 6 determine functional sugar puckers in group II intron splicing

Although group II intron ribozymes are intensively studied the question how structural dynamics affects splicing catalysis has remained elusive. We report for the first time that the group II intron domain 6 exists in a secondary structure equilibrium between a single- and a two-nucleotide bulge conformation, which is directly linked to a switch between sugar puckers of the branch site adenosine. Our study determined a functional sugar pucker equilibrium between the transesterification active C2'-endo conformation of the branch site adenosine in the 1nt bulge and an inactive C3'-endo state in the 2nt bulge fold, allowing the group II intron to switch its activity from the branching to the exon ligation step. Our detailed NMR spectroscopic investigation identified magnesium (II) ions and the branching reaction as regulators of the equilibrium populations. The tuneable secondary structure/sugar pucker equilibrium supports a conformational selection mechanism to up- and downregulate catalytically active and inactive states of the branch site adenosine to orchestrate the multi-step splicing process. The conformational dynamics of group II intron domain 6 is also proposed to be a key aspect for the directionality selection in reversible splicing.

The crystal structure of a tetrahydrofolate-bound dihydrofolate reductase reveals the origin of slow product release

Dihydrofolate reductase (DHFR) catalyzes the stereospecific reduction of 7,8-dihydrofolate (FH2) to (6s)-5,6,7,8-tetrahydrofolate (FH4) via hydride transfer from NADPH. The consensus Escherichia coli DHFR mechanism involves conformational changes between closed and occluded states occurring during the rate-limiting product release step. Although the Protein Data Bank (PDB) contains over 250 DHFR structures, the FH4 complex structure responsible for rate-limiting product release is unknown. We report to our knowledge the first crystal structure of an E. coli. DHFR:FH4 complex at 1.03 Å resolution showing distinct stabilizing interactions absent in FH2 or related (6R)-5,10-dideaza-FH4 complexes. We discover the time course of decay of the co-purified endogenous FH4 during crystal growth, with conversion from FH4 to FH2 occurring in 2–3 days. We also determine another occluded complex structure of E. coli DHFR with a slow-onset nanomolar inhibitor that contrasts with the methotrexate complex, suggesting a plausible strategy for designing DHFR antibiotics by targeting FH4 product conformations.

Hongnan Cao et al. present the X-ray crystal structure of E. coli dihydrofolate reductase (DHFR) in complex with its reduced substrate, (6s)-5,6,7,8-tetrahydrofolate (FH4). This structure provides the first glimpse of the rate-limiting product release step of the DHFR mechanism and suggests a strategy for designing DHFR-targeting antibiotics.

Chemical Exchange

The phenomenon of chemical or conformational exchange in NMR spectroscopy has enabled detailed characterization of time-dependent aspects of biomolecular function, including folding, molecular recognition, allostery, and catalysis, on timescales from microsecond to second. Importantly, NMR methods based on a variety of spin relaxation parameters have been developed that provide quantitative information on interconversion kinetics, thermodynamic properties, and structural features of molecular states populated to a fraction of a percent at equilibrium and otherwise unobservable by other NMR approaches. The ongoing development of more sophisticated experimental techniques and the necessity to apply these methods to larger and more complex molecular systems engenders a corresponding need for theoretical advances describing such techniques and facilitating data analysis in applications. This review surveys current aspects of the theory of chemical exchange, as utilized in ZZ-exchange; Hahn and Carr-Purcell-Meiboom-Gill (CPMG) spin-echo; and R_1ρ, chemical exchange saturation transfer (CEST), and dark state saturation transfer (DEST) spin-locking experiments. The review emphasizes theoretical results for kinetic topologies with more than two interconverting states, both to obtain compact analytical forms suitable for data analysis and to establish conditions for distinguishability between alternative kinetic schemes.

High-pressure protein crystal structure analysis of Escherichia coli dihydrofolate reductase complexed with folate and NADP. Acta Crystallogr D Struct Biol 2018;74:895-905. [PMID: 30198899 PMCID: PMC6130465 DOI: 10.1107/s2059798318009397]

A high-pressure crystallographic study was conducted on Escherichia coli dihydrofolate reductase (ecDHFR) complexed with folate and NADP+ in crystal forms containing both the open and closed conformations of the M20 loop under high-pressure conditions of up to 800 MPa. At pressures between 270 and 500 MPa the crystal form containing the open conformation exhibited a phase transition from P21 to C2. Several structural changes in ecDHFR were observed at high pressure that were also accompanied by structural changes in the NADP+ cofactor and the hydration structure. In the crystal form with the closed conformation the M20 loop moved as the pressure changed, with accompanying conformational changes around the active site, including NADP+ and folate. These movements were consistent with the suggested hypothesis that movement of the M20 loop was necessary for ecDHFR to catalyze the reaction. In the crystal form with the open conformation the nicotinamide ring of the NADP+ cofactor undergoes a large flip as an intermediate step in the reaction, despite being in a crystalline state. Furthermore, observation of the water molecules between Arg57 and folate elucidated an early step in the substrate-binding pathway. These results demonstrate the possibility of using high-pressure protein crystallography as a method to capture high-energy substates or transient structures related to the protein reaction cycle.

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.