Automated NMR resonance assignments and structure determination using a minimal set of 4D spectra

Elucidation of the Mechanisms of Inter-domain Coupling in the Monomeric State of Enzyme I by High-pressure NMR

The catalytic cycle of Enzyme I (EI), a phosphotransferase enzyme responsible for converting phosphoenolpyruvate (PEP) into pyruvate, is characterized by a series of local and global conformational rearrangements. This multistep process includes a monomer-to-dimer transition, followed by an open-to-closed rearrangement of the dimeric complex upon PEP binding. In the present study, we investigate the thermodynamics of EI dimerization using a range of high-pressure solution NMR techniques complemented by SAXS experiments. 1H-15N TROSY and 1H-13C methyl TROSY NMR spectra combined with 15N relaxation measurements revealed that a native-like engineered variant of full-length EI fully dissociates into stable monomeric state above 1.5 kbar. Conformational ensembles of EI monomeric state were generated via a recently developed protocol combining coarse-grained molecular simulations with experimental backbone residual dipolar coupling measurements. Analysis of the structural ensembles provided detailed insights into the molecular mechanisms driving formation of the catalytically competent dimeric state, and reveals that each step of EI catalytical cycle is associated with a significant reduction in either inter- or intra-domain conformational entropy. Altogether, this study completes a large body work conducted by our group on EI and establishes a comprehensive structural and dynamical description of the catalytic cycle of this prototypical multidomain, oligomeric enzyme.

Conformational plasticity of RAS Q61 family of neoepitopes results in distinct features for targeted recognition

The conformational landscapes of peptide/human leucocyte antigen (pHLA) protein complexes encompassing tumor neoantigens provide a rationale for target selection towards autologous T cell, vaccine, and antibody-based therapeutic modalities. Here, using complementary biophysical and computational methods, we characterize recurrent RAS55-64 Q61 neoepitopes presented by the common HLA-A*01:01 allotype. We integrate sparse NMR restraints with Rosetta docking to determine the solution structure of NRASQ61K/HLA-A*01:01, which enables modeling of other common RAS55-64 neoepitopes. Hydrogen/deuterium exchange mass spectrometry experiments alongside molecular dynamics simulations reveal differences in solvent accessibility and conformational plasticity across a panel of common Q61 neoepitopes that are relevant for recognition by immunoreceptors. Finally, we predict binding and provide structural models of NRASQ61K antigens spanning the entire HLA allelic landscape, together with in vitro validation for HLA-A*01:191, HLA-B*15:01, and HLA-C*08:02. Our work provides a basis to delineate the solution surface features and immunogenicity of clinically relevant neoepitope/HLA targets for cancer therapy.

Recognition and coacervation of G-quadruplexes by a multifunctional disordered region in RECQ4 helicase

Biomolecular polyelectrolyte complexes can be formed between oppositely charged intrinsically disordered regions (IDRs) of proteins or between IDRs and nucleic acids. Highly charged IDRs are abundant in the nucleus, yet few have been functionally characterized. Here, we show that a positively charged IDR within the human ATP-dependent DNA helicase Q4 (RECQ4) forms coacervates with G-quadruplexes (G4s). We describe a three-step model of charge-driven coacervation by integrating equilibrium and kinetic binding data in a global numerical model. The oppositely charged IDR and G4 molecules form a complex in the solution that follows a rapid nucleation-growth mechanism leading to a dynamic equilibrium between dilute and condensed phases. We also discover a physical interaction with Replication Protein A (RPA) and demonstrate that the IDR can switch between the two extremes of the structural continuum of complexes. The structural, kinetic, and thermodynamic profile of its interactions revealed a dynamic disordered complex with nucleic acids and a static ordered complex with RPA protein. The two mutually exclusive binding modes suggest a regulatory role for the IDR in RECQ4 function by enabling molecular handoffs. Our study extends the functional repertoire of IDRs and demonstrates a role of polyelectrolyte complexes involved in G4 binding.

NMR-Chemical-Shift-Driven Protocol Reveals the Cofactor-Bound, Complete Structure of Dynamic Intermediates of the Catalytic Cycle of Oncogenic KRAS G12C Protein and the Significance of the Mg2+ Ion

In this work, catalytically significant states of the oncogenic G12C variant of KRAS, those of Mg2+-free and Mg2+-bound GDP-loaded forms, have been determined using CS-Rosetta software and NMR-data-driven molecular dynamics simulations. There are several Mg2+-bound G12C KRAS/GDP structures deposited in the Protein Data Bank (PDB), so this system was used as a reference, while the structure of the Mg2+-free but GDP-bound state of the RAS cycle has not been determined previously. Due to the high flexibility of the Switch-I and Switch-II regions, which also happen to be the catalytically most significant segments, only chemical shift information could be collected for the most important regions of both systems. CS-Rosetta was used to derive an "NMR ensemble" based on the measured chemical shifts, which, however, did not contain the nonprotein components of the complex. We developed a torsional restraint set for backbone torsions based on the CS-Rosetta ensembles for MD simulations, overriding the force-field-based parametrization in the presence of the reinserted cofactors. This protocol (csdMD) resulted in complete models for both systems that also retained the structural features and heterogeneity defined by the measured chemical shifts and allowed a detailed comparison of the Mg2+-bound and Mg2+-free states of G12C KRAS/GDP.

Solution Structure Ensembles of the Open and Closed Forms of the ∼130 kDa Enzyme I via AlphaFold Modeling, Coarse Grained Simulations, and NMR

Large-scale interdomain rearrangements are essential to protein function, governing the activity of large enzymes and molecular machineries. Yet, obtaining an atomic-resolution understanding of how the relative domain positioning is affected by external stimuli is a hard task in modern structural biology. Here, we show that combining structural modeling by AlphaFold2 with coarse-grained molecular dynamics simulations and NMR residual dipolar coupling data is sufficient to characterize the spatial domain organization of bacterial enzyme I (EI), a ∼130 kDa multidomain oligomeric protein that undergoes large-scale conformational changes during its catalytic cycle. In particular, we solve conformational ensembles for EI at two different experimental temperatures and demonstrate that a lower temperature favors sampling of the catalytically competent closed state of the enzyme. These results suggest a role for conformational entropy in the activation of EI and demonstrate the ability of our protocol to detect and characterize the effect of external stimuli (such as mutations, ligand binding, and post-translational modifications) on the interdomain organization of multidomain proteins. We expect the ensemble refinement protocol described here to be easily transferrable to the investigation of the structure and dynamics of other uncharted multidomain systems and have assembled a Google Colab page (https://potoyangroup.github.io/Seq2Ensemble/) to facilitate implementation of the presented methodology elsewhere.

Linear discriminant analysis reveals hidden patterns in NMR chemical shifts of intrinsically disordered proteins

NMR spectroscopy is key in the study of intrinsically disordered proteins (IDPs). Yet, even the first step in such an analysis-the assignment of observed resonances to particular nuclei-is often problematic due to low peak dispersion in the spectra of IDPs. We show that the assignment process can be aided by finding "hidden" chemical shift patterns specific to the amino acid residue types. We find such patterns in the training data from the Biological Magnetic Resonance Bank using linear discriminant analysis, and then use them to classify spin systems in an α-synuclein sample prepared by us. We describe two situations in which the procedure can greatly facilitate the analysis of NMR spectra. The first involves the mapping of spin systems chains onto the protein sequence, which is part of the assignment procedure-a prerequisite for any NMR-based protein analysis. In the second, the method supports assignment transfer between similar samples. We conducted experiments to demonstrate these cases, and both times the majority of spin systems could be unambiguously assigned to the correct residue types.

Key aspects of the past 30 years of protein design

Proteins are the workhorse of life. They are the building infrastructure of living systems; they are the most efficient molecular machines known, and their enzymatic activity is still unmatched in versatility by any artificial system. Perhaps proteins' most remarkable feature is their modularity. The large amount of information required to specify each protein's function is analogically encoded with an alphabet of just ∼20 letters. The protein folding problem is how to encode all such information in a sequence of 20 letters. In this review, we go through the last 30 years of research to summarize the state of the art and highlight some applications related to fundamental problems of protein evolution.

A Single Point Mutation Controls the Rate of Interconversion Between the g + and g - Rotamers of the Histidine 189 χ2 Angle That Activates Bacterial Enzyme I for Catalysis

Enzyme I (EI) of the bacterial phosphotransferase system (PTS) is a master regulator of bacterial metabolism and a promising target for development of a new class of broad-spectrum antibiotics. The catalytic activity of EI is mediated by several intradomain, interdomain, and intersubunit conformational equilibria. Therefore, in addition to its relevance as a drug target, EI is also a good model for investigating the dynamics/function relationship in multidomain, oligomeric proteins. Here, we use solution NMR and protein design to investigate how the conformational dynamics occurring within the N-terminal domain (EIN) affect the activity of EI. We show that the rotameric g⁺-to-g⁻ transition of the active site residue His¹⁸⁹ χ2 angle is decoupled from the state A-to-state B transition that describes a ∼90° rigid-body rearrangement of the EIN subdomains upon transition of the full-length enzyme to its catalytically competent closed form. In addition, we engineered EIN constructs with modulated conformational dynamics by hybridizing EIN from mesophilic and thermophilic species, and used these chimeras to assess the effect of increased or decreased active site flexibility on the enzymatic activity of EI. Our results indicate that the rate of the autophosphorylation reaction catalyzed by EI is independent from the kinetics of the g⁺-to-g⁻ rotameric transition that exposes the phosphorylation site on EIN to the incoming phosphoryl group. In addition, our work provides an example of how engineering of hybrid mesophilic/thermophilic chimeras can assist investigations of the dynamics/function relationship in proteins, therefore opening new possibilities in biophysics.

Distinct EH domains of the endocytic TPLATE complex confer lipid and protein binding

Clathrin-mediated endocytosis (CME) is the gatekeeper of the plasma membrane. In contrast to animals and yeasts, CME in plants depends on the TPLATE complex (TPC), an evolutionary ancient adaptor complex. However, the mechanistic contribution of the individual TPC subunits to plant CME remains elusive. In this study, we used a multidisciplinary approach to elucidate the structural and functional roles of the evolutionary conserved N-terminal Eps15 homology (EH) domains of the TPC subunit AtEH1/Pan1. By integrating high-resolution structural information obtained by X-ray crystallography and NMR spectroscopy with all-atom molecular dynamics simulations, we provide structural insight into the function of both EH domains. Both domains bind phosphatidic acid with a different strength, and only the second domain binds phosphatidylinositol 4,5-bisphosphate. Unbiased peptidome profiling by mass-spectrometry revealed that the first EH domain preferentially interacts with the double N-terminal NPF motif of a previously unidentified TPC interactor, the integral membrane protein Secretory Carrier Membrane Protein 5 (SCAMP5). Furthermore, we show that AtEH/Pan1 proteins control the internalization of SCAMP5 via this double NPF peptide interaction motif. Collectively, our structural and functional studies reveal distinct but complementary roles of the EH domains of AtEH/Pan1 in plant CME and connect the internalization of SCAMP5 to the TPLATE complex.

Natural Products, the Fourth Industrial Revolution, and the Quintuple Helix

The profound interconnectedness of the sciences and technologies embodied in the Fourth Industrial Revolution is discussed in terms of the global role of natural products, and how that interplays with the development of sustainable and climate-conscious practices of cyberecoethnopharmacolomics within the Quintuple Helix for the promotion of a healthier planet and society.

Highly Efficient Determination of Complex NMR Multiplet Structures in Inhomogeneous Magnetic Fields

Proton-proton scalar (J) coupling plays an important role in disentangling molecular structures and spatial conformations. But it is challenging to extract J coupling networks from congested ¹H NMR spectra, especially in inhomogeneous magnetic fields. Herein, we propose a general liquid NMR protocol, named HR-G-SERF, to implement highly efficient determination of individual J couplings and corresponding coupling networks via simultaneously suppressing effects of spectral congestions and magnetic field inhomogeneity. This method records full-resolved 2D absorption-mode spectra to deliver great convenience for multipet analyses on complex samples. More meaningfully, it is capable of disentangling multiplet structures of biological samples, that is, grape sarcocarp, despite of its heterogeneous semisolid state and extensive compositions. In addition, a modification, named AH-G-SERF, is developed to compress experimental acquisition and subsequently improve unit-time SNR, while maintaining satisfactory spectral performance. This accelerated variant may further boost the applicability for rapid NMR detections and afford the possibility of adopting hyperpolarized substances to enhance the overall sensitivity. Therefore, this study provides a promising tool for molecular structure elucidations and composition analyses in chemistry, biochemistry, and metabonomics among others.

Phosphorylation-induced changes in the PDZ domain of Dishevelled 3

The PDZ domain of Dishevelled 3 protein belongs to a highly abundant protein recognition motif which typically binds short C-terminal peptides. The affinity of the PDZ towards the peptides could be fine-tuned by a variety of post-translation modifications including phosphorylation. However, how phosphorylations affect the PDZ structure and its interactions with ligands remains elusive. Combining molecular dynamics simulations, NMR titration, and biological experiments, we explored the role of previously reported phosphorylation sites and their mimetics in the Dishevelled PDZ domain. Our observations suggest three major roles for phosphorylations: (1) acting as an on/off PDZ binding switch, (2) allosterically affecting the binding groove, and (3) influencing the secondary binding site. Our simulations indicated that mimetics had similar but weaker effects, and the effects of distinct sites were non-additive. This study provides insight into the Dishevelled regulation by PDZ phosphorylation. Furthermore, the observed effects could be used to elucidate the regulation mechanisms in other PDZ domains.

Macromolecular modeling and design in Rosetta: recent methods and frameworks

The Rosetta software for macromolecular modeling, docking and design is extensively used in laboratories worldwide. During two decades of development by a community of laboratories at more than 60 institutions, Rosetta has been continuously refactored and extended. Its advantages are its performance and interoperability between broad modeling capabilities. Here we review tools developed in the last 5 years, including over 80 methods. We discuss improvements to the score function, user interfaces and usability. Rosetta is available at http://www.rosettacommons.org.

A Practical Perspective on the Roles of Solution NMR Spectroscopy in Drug Discovery

Solution nuclear magnetic resonance (NMR) spectroscopy is a powerful tool to study structures and dynamics of biomolecules under physiological conditions. As there are numerous NMR-derived methods applicable to probe protein–ligand interactions, NMR has been widely utilized in drug discovery, especially in such steps as hit identification and lead optimization. NMR is frequently used to locate ligand-binding sites on a target protein and to determine ligand binding modes. NMR spectroscopy is also a unique tool in fragment-based drug design (FBDD), as it is able to investigate target-ligand interactions with diverse binding affinities. NMR spectroscopy is able to identify fragments that bind weakly to a target, making it valuable for identifying hits targeting undruggable sites. In this review, we summarize the roles of solution NMR spectroscopy in drug discovery. We describe some methods that are used in identifying fragments, understanding the mechanism of action for a ligand, and monitoring the conformational changes of a target induced by ligand binding. A number of studies have proven that ¹⁹F-NMR is very powerful in screening fragments and detecting protein conformational changes. In-cell NMR will also play important roles in drug discovery by elucidating protein-ligand interactions in living cells.

Hybrid Thermophilic/Mesophilic Enzymes Reveal a Role for Conformational Disorder in Regulation of Bacterial Enzyme I

Conformational disorder is emerging as an important feature of biopolymers, regulating a vast array of cellular functions, including signaling, phase separation, and enzyme catalysis. Here we combine NMR, crystallography, computer simulations, protein engineering, and functional assays to investigate the role played by conformational heterogeneity in determining the activity of the C-terminal domain of bacterial Enzyme I (EIC). In particular, we design chimeric proteins by hybridizing EIC from thermophilic and mesophilic organisms, and we characterize the resulting constructs for structure, dynamics, and biological function. We show that EIC exists as a mixture of active and inactive conformations and that functional regulation is achieved by tuning the thermodynamic balance between active and inactive states. Interestingly, we also present a hybrid thermophilic/mesophilic enzyme that is thermostable and more active than the wild-type thermophilic enzyme, suggesting that hybridizing thermophilic and mesophilic proteins is a valid strategy to engineer thermostable enzymes with significant low-temperature activity.

Comparative phosphorylation map of Dishevelled 3 links phospho-signatures to biological outputs

Background

Dishevelled (DVL) is an essential component of the Wnt signaling cascades. Function of DVL is controlled by phosphorylation but the molecular details are missing. DVL3 contains 131 serines and threonines whose phosphorylation generates complex barcodes underlying diverse DVL3 functions. In order to dissect the role of DVL phosphorylation we analyzed the phosphorylation of human DVL3 induced by previously reported (CK1ε, NEK2, PLK1, CK2α, RIPK4, PKCδ) and newly identified (TTBK2, Aurora A) DVL kinases.

Methods

Shotgun proteomics including TiO₂ enrichment of phosphorylated peptides followed by liquid chromatography tandem mass spectrometry on immunoprecipitates from HEK293T cells was used to identify and quantify phosphorylation of DVL3 protein induced by 8 kinases. Functional characterization was performed by in-cell analysis of phospho-mimicking/non-phosphorylatable DVL3 mutants and supported by FRET assays and NMR spectroscopy.

Results

We used quantitative mass spectrometry and calculated site occupancies and quantified phosphorylation of > 80 residues. Functional validation demonstrated the importance of CK1ε-induced phosphorylation of S268 and S311 for Wnt-3a-induced β-catenin activation. S630–643 cluster phosphorylation by CK1, NEK2 or TTBK2 is essential for even subcellular distribution of DVL3 when induced by CK1 and TTBK2 but not by NEK2. Further investigation showed that NEK2 utilizes a different mechanism to promote even localization of DVL3. NEK2 triggered phosphorylation of PDZ domain at S263 and S280 prevents binding of DVL C-terminus to PDZ and promotes an open conformation of DVL3 that is more prone to even subcellular localization.

Conclusions

We identify unique phosphorylation barcodes associated with DVL function. Our data provide an example of functional synergy between phosphorylation in structured domains and unstructured IDRs that together dictate the biological outcome.

High-Resolution Probing of Heterogeneous Samples by Spatially Selective Pure Shift NMR Spectroscopy

Liquid NMR spectroscopy generally encounters two major challenges for high-resolution measurements of heterogeneous samples, namely, magnetic field inhomogeneity caused by spatial variations in magnetic susceptibility and spectral congestion induced by crowded NMR resonances. In this study, we demonstrate a spatially selective pure shift NMR approach for high-resolution probing of heterogeneous samples by suppressing effects of field inhomogeneity and J coupling simultaneously. A Fourier phase encoding strategy is proposed and implemented for spatially selective pure shift experiments to enhance signal intensity and further boost the applicability. The spatially selective pure shift method can serve as an effective tool for high-resolution probing of heterogeneous samples, thus presenting interesting prospects for extensive applications in the fields of chemistry, physics, biology, and food science.

N-terminal fusion of the N-terminal domain of bacterial enzyme I facilitates recombinant expression and purification of the human RNA demethylases FTO and Alkbh5

Various fusion tags are commonly employed to increase the heterologous expression and solubility of aggregation-prone proteins within Escherichia coli. Herein, we present a protocol for efficient recombinant expression and purification of the human RNA demethylases Alkbh5 and FTO. Our method incorporates a novel fusion tag (the N-terminal domain of bacterial enzyme I, EIN) that dramatically increases the solubility of its fusion partner and is promptly removed upon digestion with a protease. The presented protocol allows for the production of mg amounts of Alkbh5 and FTO in 1L of both rich and minimal media. We developed a liquid chromatography-mass spectrometry (LC-MS)-based assay to confirm that both proteins are enzymatically active. Furthermore, the LC-MS method developed here is applicable to other members of the AlkB family of Fe(II)/α-ketoglutarate-dependent dioxygenases. The superior protein yield, afforded by our expression and purification method, will facilitate biochemical investigations into the biological function of the human RNA demethylases and endorse employment of EIN as a broadly applicable fusion tag for recombinant expression projects.

Resonance assignment of the 128 kDa enzyme I dimer from Thermoanaerobacter tengcongensis

Enzyme I (EI) of the bacterial phosphotransferase system (PTS) utilizes phosphoenolpyruvate (PEP) as a source of energy in order to transport sugars across the cellular membrane. PEP binding to EI initiates a phosphorylation cascade that regulates a variety of essential pathways in the metabolism of bacterial cells. Given its central role in controlling bacterial metabolism, EI has been often suggested as a good target for antimicrobial research. Here, we report the ¹H_N, ¹⁵N, ¹³C', ¹H_methyl, and ¹³C_methyl chemical shifts of the 128 kDa homodimer EI from the thermophile Thermoanaerobacter tengcongensis. In total 79% of the expected backbone amide correlations and 80% of the expected methyl TROSY peaks from U-[²H, ¹³C, ¹⁵N], Ile^δ1-[¹³CH₃], Val-Leu-[¹³CH₃/¹²CD₃] labeled EI were assigned. The reported assignments will enable future structural studies aimed at illuminating the fundamental mechanisms governing long-range interdomain communication in EI and at indicating new therapeutic strategies to combat bacterial infections.

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

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

Automated, context-free assignment of asymmetric rotor microwave spectra

We present a new algorithm, Robust Automated Assignment of Rigid Rotors (RAARR), for assigning rotational spectra of asymmetric tops. The RAARR algorithm can automatically assign experimental spectra under a broad range of conditions, including spectra comprised of multiple mixture components, in ≲100 s. The RAARR algorithm exploits constraints placed by the conservation of energy to find sets of connected lines in an unassigned spectrum. The highly constrained structure of these sets eliminates all but a handful of plausible assignments for a given set, greatly reducing the number of potential assignments that must be evaluated. We successfully apply our algorithm to automatically assign 15 experimental spectra, including 5 previously unassigned species, without prior estimation of molecular rotational constants. In 9 of the 15 cases, the RAARR algorithm successfully assigns two or more mixture components.

Dishevelled-3 conformation dynamics analyzed by FRET-based biosensors reveals a key role of casein kinase 1

Dishevelled (DVL) is the key component of the Wnt signaling pathway. Currently, DVL conformational dynamics under native conditions is unknown. To overcome this limitation, we develop the Fluorescein Arsenical Hairpin Binder- (FlAsH-) based FRET in vivo approach to study DVL conformation in living cells. Using this single-cell FRET approach, we demonstrate that (i) Wnt ligands induce open DVL conformation, (ii) DVL variants that are predominantly open, show more even subcellular localization and more efficient membrane recruitment by Frizzled (FZD) and (iii) Casein kinase 1 ɛ (CK1ɛ) has a key regulatory function in DVL conformational dynamics. In silico modeling and in vitro biophysical methods explain how CK1ɛ-specific phosphorylation events control DVL conformations via modulation of the PDZ domain and its interaction with DVL C-terminus. In summary, our study describes an experimental tool for DVL conformational sampling in living cells and elucidates the essential regulatory role of CK1ɛ in DVL conformational dynamics.

De novo design of a non-local β-sheet protein with high stability and accuracy

β-sheet proteins carry out critical functions in biology, and hence are attractive scaffolds for computational protein design. Despite this potential, de novo design of all-β-sheet proteins from first principles lags far behind the design of all-α or mixed-αβ domains owing to their non-local nature and the tendency of exposed β-strand edges to aggregate. Through study of loops connecting unpaired β-strands (β-arches), we have identified a series of structural relationships between loop geometry, side chain directionality and β-strand length that arise from hydrogen bonding and packing constraints on regular β-sheet structures. We use these rules to de novo design jellyroll structures with double-stranded β-helices formed by eight antiparallel β-strands. The nuclear magnetic resonance structure of a hyperthermostable design closely matched the computational model, demonstrating accurate control over the β-sheet structure and loop geometry. Our results open the door to the design of a broad range of non-local β-sheet protein structures.

CS-ROSETTA

Chemical Shift-Rosetta (CS-Rosetta) is an automated method that employs NMR chemical shifts to model protein structures de novo. In this chapter, we introduce the terminology and central concepts of CS-Rosetta. We describe the architecture and functionality of automatic NOESY assignment (AutoNOE) and structure determination protocols (Abrelax and RASREC) within the CS-Rosetta framework. We further demonstrate how CS-Rosetta can discriminate near-native structures against a large conformational search space using restraints obtained from NMR data, and/or sequence and structure homology. We highlight how CS-Rosetta can be combined with alternative automated approaches to (i) model oligomeric systems and (ii) create NMR-based structure determination pipeline. To show its practical applicability, we emphasize on the computational requirements and performance of CS-Rosetta for protein targets of varying molecular weight and complexity. Finally, we discuss the current Python interface, which enables easy execution of protocols for rapid and accurate high-resolution structure determination.

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.