Hydrogen bond formation may enhance RDC-based discrimination of enantiomers.

The distinction of enantiomers based on residual anisotropic parameters obtained by alignment in chiral poly-γ-benzyl-L-glutamate (PBLG) is among the strongest in high-resolution NMR spectroscopy. However, large variations in enantiodifferentiation among different solutes are frequently observed. One hypothesis is that the formation of hydrogen bonds between solute and PBLG is important for the distinction of enantiomers. With a small set of three almost spherical enantiomeric pairs, for which 1DCH residual dipolar couplings are measured, we address this issue in a systematic way: borneol contains a single functional group that can act as a hydrogen bond donor, camphor has a single group that may act as a hydrogen bond acceptor, and quinuclidinol can act as both hydrogen bond donor and acceptor. The results are unambiguous: although camphor shows low enantiodifferentiation with PBLG and alignment that can be predicted well by the purely steric TRAMITE approach, the distinction of enantiomers for the other enantiomeric pairs is significantly higher with alignment properties that must involve a specific interaction in addition to steric alignment.

Intrinsically Disordered Regions in the Transcription Factor MYC:MAX Modulate DNA Binding via Intramolecular Interactions.

The basic helix-loop-helix leucine zipper (bHLH-LZ) transcription factor (TF) MYC is in large part an intrinsically disordered oncoprotein. In complex with its obligate heterodimerization partner MAX, MYC preferentially binds E-Box DNA sequences (CANNTG). At promoters containing these sequence motifs, MYC controls fundamental cellular processes such as cell cycle progression, metabolism, and apoptosis. A vast network of proteins in turn regulates MYC function via intermolecular interactions. In this work, we establish another layer of MYC regulation by intramolecular interactions. We used nuclear magnetic resonance (NMR) spectroscopy to identify and map multiple binding sites for the C-terminal MYC:MAX DNA-binding domain (DBD) on the intrinsically disordered regions (IDRs) in the MYC N-terminus. We find that these binding events in trans are driven by electrostatic attraction, that they have distinct affinities, and that they are competitive with DNA binding. Thereby, we observe the strongest effects for the N-terminal MYC box 0 (Mb0), a conserved motif involved in MYC transactivation and target gene induction. We prepared recombinant full-length MYC:MAX complex and demonstrate that the interactions identified in this work are also relevant in cis, i.e., as intramolecular interactions. These findings are supported by surface plasmon resonance (SPR) experiments, which revealed that intramolecular IDR:DBD interactions in MYC decelerate the association of MYC:MAX complexes to DNA. Our work offers new insights into how bHLH-LZ TFs are regulated by intramolecular interactions, which open up new possibilities for drug discovery.

Conformations of Macrocyclic Peptides Sampled by Nuclear Magnetic Resonance: Models for Cell-Permeability.

The biological activities and pharmacological properties of peptides and peptide mimetics are determined by their conformational states. Therefore, a detailed understanding of the conformational landscape is crucial for rational drug design. Nuclear magnetic resonance (NMR) is the only method for structure determination in solution. However, it remains challenging to determine the structures of peptides using NMR because of very weak nuclear Overhauser effects (NOEs), the semiquantitative nature of the rotating frame Overhauser effect (ROE), and the low number of NOEs/ROEs in N-methylated peptides. In this study, we introduce a new approach to investigating the structures of modified macrocyclic peptides. We utilize exact NOEs (eNOEs) in viscous solvent mixtures to replicate various cellular environments. eNOEs provide detailed structural information for highly dynamic modified peptides. Structures of high precision were obtained for cyclosporin A, with a backbone atom rmsd of 0.10 Å. Distinct conformational states in different environments were identified for omphalotin A (OmphA), a fungal nematotoxic and multiple backbone N-methylated macrocyclic peptides. A model for cell-permeation is presented for OmphA, based on its structures in polar, apolar, and mixed polarity solvents. During the transition from a polar to an apolar environment, OmphA undergoes a rearrangement of its H-bonding network, accompanied by a cis to trans isomerization of the ω torsion angle within a type VIa ß-turn. We hypothesize that the kinetics of these conformational transitions play a crucial role in determining the membrane-permeation capabilities of OmphA.

Affinity measurement of strong ligands with NMR spectroscopy: Limitations and ways to overcome them.

NMR spectroscopy is currently extensively used in binding assays for hit identification, but its use in dissociation constant determination is more limited when compared to other biophysical techniques, in particular for tight binders. Although NMR is quite suitable for measuring the binding strength of weak to medium affinity ligands with dissociation constant KD > 1 µM, it has some limitations in the determination of the binding strength of tight binders (KD < 1 µM). A theoretical analysis of the binding affinity determination of strong ligands using different types of NMR experiments is provided and practical guidelines are given for overcoming the limitations and for the proper set-up of the experiments. Some approaches require reagents with unique properties or highly specialized equipment, while others can be applied quite generally. We describe all approaches in detail, but give higher emphasis to the more general methods, like competition experiments, where we include actual experimental data and discuss the practical aspects.

A hybrid structure determination approach to investigate the druggability of the nucleocapsid protein of SARS-CoV-2.

The pandemic caused by SARS-CoV-2 has called for concerted efforts to generate new insights into the biology of betacoronaviruses to inform drug screening and development. Here, we establish a workflow to determine the RNA recognition and druggability of the nucleocapsid N-protein of SARS-CoV-2, a highly abundant protein crucial for the viral life cycle. We use a synergistic method that combines NMR spectroscopy and protein-RNA cross-linking coupled to mass spectrometry to quickly determine the RNA binding of two RNA recognition domains of the N-protein. Finally, we explore the druggability of these domains by performing an NMR fragment screening. This workflow identified small molecule chemotypes that bind to RNA binding interfaces and that have promising properties for further fragment expansion and drug development.

Nanobody GPS by PCS: An Efficient New NMR Analysis Method for G Protein Coupled Receptors and Other Large Proteins.

NMR chemical shift changes can report on the functional dynamics of biomacromolecules in solution with sizes >1 MDa. However, their interpretation requires chemical shift assignments to individual nuclei, which for large molecules often can only be obtained by tedious point mutations that may interfere with function. We present here an efficient pseudocontact shift NMR method to assign biomacromolecules using bound antibodies tagged with lanthanoid DOTA chelators. The stability of the antibody allows positioning the DOTA tag at many surface sites, providing triangulation of the macromolecule nuclei at distances >60 Å. The method provides complete assignments of valine and tyrosine 1H-15N resonances of the ß1-adrenergic receptor in various functional forms. The detected chemical shift changes reveal strong forces exerted onto the backbone of transmembrane helix 3 during signal transmission, which are absorbed by its electronic structure. The assignment method is applicable to any soluble biomacromolecule for which suitable complementary binders exist.

The Disordered MAX N-terminus Modulates DNA Binding of the Transcription Factor MYC:MAX.

The intrinsically disordered protein MYC belongs to the family of basic helix-loop-helix leucine zipper (bHLH-LZ) transcription factors (TFs). In complex with its cognate binding partner MAX, MYC preferentially binds to E-Box promotor sequences where it controls fundamental cellular processes such as cell cycle progression, metabolism, and apoptosis. Intramolecular regulation of MYC:MAX has not yet been investigated in detail. In this work, we use Nuclear Magnetic Resonance (NMR) spectroscopy to identify and map interactions between the disordered MAX N-terminus and the MYC:MAX DNA binding domain (DBD). We find that this binding event is mainly driven by electrostatic interactions and that it is competitive with DNA binding. Using NMR spectroscopy and Surface Plasmon Resonance (SPR), we demonstrate that the MAX N-terminus serves to accelerate DNA binding kinetics of MYC:MAX and MAX:MAX dimers, while it simultaneously provides specificity for E-Box DNA. We also establish that these effects are further enhanced by Casein Kinase 2-mediated phosphorylation of two serine residues in the MAX N-terminus. Our work provides new insights how bHLH-LZ TFs are regulated by intramolecular interactions between disordered regions and the folded DNA binding domain.

Biosynthetic potential of the global ocean microbiome.

Natural microbial communities are phylogenetically and metabolically diverse. In addition to underexplored organismal groups1, this diversity encompasses a rich discovery potential for ecologically and biotechnologically relevant enzymes and biochemical compounds2,3. However, studying this diversity to identify genomic pathways for the synthesis of such compounds4 and assigning them to their respective hosts remains challenging. The biosynthetic potential of microorganisms in the open ocean remains largely uncharted owing to limitations in the analysis of genome-resolved data at the global scale. Here we investigated the diversity and novelty of biosynthetic gene clusters in the ocean by integrating around 10,000 microbial genomes from cultivated and single cells with more than 25,000 newly reconstructed draft genomes from more than 1,000 seawater samples. These efforts revealed approximately 40,000 putative mostly new biosynthetic gene clusters, several of which were found in previously unsuspected phylogenetic groups. Among these groups, we identified a lineage rich in biosynthetic gene clusters ('Candidatus Eudoremicrobiaceae') that belongs to an uncultivated bacterial phylum and includes some of the most biosynthetically diverse microorganisms in this environment. From these, we characterized the phospeptin and pythonamide pathways, revealing cases of unusual bioactive compound structure and enzymology, respectively. Together, this research demonstrates how microbiomics-driven strategies can enable the investigation of previously undescribed enzymes and natural products in underexplored microbial groups and environments.

Relaxation optimized double acquisition (RODA) as an alternative for virtual decoupling of NMR spectra.

We introduce an alternative way for spin-state selection, RODA, which yields higher sensitivity for spin systems exhibiting a TROSY effect. With RODA, the TROSY component of a doublet is recorded twice using a double acquisition scheme. RODA works by simple addition of consecutive NMR signals, and does not require any special processing. Thus, this pulse sequence element can seamlessly be integrated into existing experiments. We demonstrate the broad applicability of RODA with several systems exhibiting a TROSY effect on 15N-1H, 19F-13C or 1H-13C moieties. Further, we show that virtual decoupling with increased sensitivity is possible in a single double acquisition experiment in situations as encountered with dissolution DNP.

The molecular basis of FimT-mediated DNA uptake during bacterial natural transformation.

Naturally competent bacteria encode sophisticated protein machinery for the uptake and translocation of exogenous DNA into the cell. If this DNA is integrated into the bacterial genome, the bacterium is said to be naturally transformed. Most competent bacterial species utilise type IV pili for the initial DNA uptake step. These proteinaceous cell-surface structures are composed of thousands of pilus subunits (pilins), designated as major or minor according to their relative abundance in the pilus. Here, we show that the minor pilin FimT plays an important role in the natural transformation of Legionella pneumophila. We use NMR spectroscopy, in vitro DNA binding assays and in vivo transformation assays to understand the molecular basis of FimT's role in this process. FimT binds to DNA via an electropositive patch, rich in arginines, several of which are well-conserved and located in a conformationally flexible C-terminal tail. FimT orthologues from other Gammaproteobacteria share the ability to bind to DNA. Our results suggest that FimT plays an important role in DNA uptake in a wide range of competent species.

The trans-to-cis proline isomerization in E. coli Trx folding is accelerated by trans prolines.

The conserved fold of thioredoxin (Trx)-like thiol/disulfide oxidoreductases contains an invariant cis-proline residue (P76 in Escherichia coli Trx) that is essential for Trx function and that is responsible for the folding rate-limiting step. E. coli Trx contains four additional prolines, which are all in the trans conformation in the native state. Notably, a recent study revealed that replacement of all four trans prolines in Trx by alanines (Trx variant Trx1P) further slowed the rate-limiting step 25-fold, indicating that one or several of the four trans prolines accelerate the trans-to-cis transition of P76 in Trx wild-type (wt). Here, we characterized the folding kinetics of Trx variants containing cisP76 and one or several of the natural trans prolines of Trx wt with NMR spectroscopy. First, we demonstrate that the isomerization reaction in Trx1P is a pure two-state transition between two distinct tertiary structures, in which all observed NMR resonances changes follow the same first-order kinetics. Moreover, we show that trans-P68 is the critical residue responsible for the faster folding of wt Trx relative to the single-proline (P76) variant Trx1P, as the two-proline variant Trx2P(P76P68) already folds seven times faster than Trx1P. trans-P34 also accelerates trans-to-cis isomerization of P76, albeit to a smaller extent. Overall, the results demonstrate that trans prolines can significantly modulate the kinetics of rate-limiting trans-to-cis proline isomerization in protein folding. Finally, we discuss possible mechanisms of acceleration and the potential significance of a protein-internal folding acceleration mechanism for Trx in a living cell.

GPCR Activation States Induced by Nanobodies and Mini-G Proteins Compared by NMR Spectroscopy.

In this work, we examine methyl nuclear magnetic resonance (NMR) spectra of the methionine ε-[13CH3] labelled thermostabilized ß1 adrenergic receptor from turkey in association with a variety of different effectors, including mini-Gs and nanobody 60 (Nb60), which have not been previously studied in complex with ß1 adrenergic receptor (ß1AR) by NMR. Complexes with pindolol and Nb60 induce highly similar inactive states of the receptor, closely resembling the resting state conformational ensemble. We show that, upon binding of mini-Gs or nanobody 80 (Nb80), large allosteric changes throughout the receptor take place. The conformation of tß1AR stabilized by the native-like mini-Gs protein is highly similar to the conformation induced by the currently used surrogate Nb80. Interestingly, in both cases residual dynamics are present, which were not observed in the resting states. Finally, we reproduce a pharmaceutically relevant situation, where an antagonist abolishes the interaction of the receptor with the mini-G protein in a competitive manner, validating the functional integrity of our preparation. The presented system is therefore well suited for reproducing the individual steps of the activation cycle of a G protein-coupled receptor (GPCR) in vitro and serves as a basis for functional and pharmacological characterizations of more native-like systems in the future.

Enabling NMR Studies of High Molecular Weight Systems Without the Need for Deuteration: The XL-ALSOFAST Experiment with Delayed Decoupling.

Current biological research increasingly focusses on large human proteins and their complexes. Such proteins are difficult to study by NMR spectroscopy because they often can only be produced in higher eukaryotic expression systems, where deuteration is hardly feasible. Here, we present the XL-ALSOFAST-[13 C,1 H]-HMQC experiment with much improved sensitivity for fully protonated high molecular weight proteins. For the tested systems ranging from 100 to 240 kDa in size, 3-fold higher sensitivity was obtained on average for fast relaxing signals compared to current state-of-the-art experiments. In the XL-ALSOFAST approach, non-observed magnetisation is optimally exploited and transverse relaxation is minimized by the newly introduced concept of delayed decoupling. The combination of high sensitivity and superior artefact suppression makes it ideal for studying inherently unstable membrane proteins or for analysing therapeutic antibodies at natural 13 C abundance. The XL-ALSOFAST and delayed decoupling will therefore expand the range of biomolecular systems accessible to NMR spectroscopy.

Determination of Configuration and Conformation of a Reserpine Derivative with Seven Stereogenic Centers Using Molecular Dynamics with RDC-Derived Tensorial Constraints*.

NMR-based determination of the configuration of complex molecules containing many stereocenters is often not possible using traditional NOE data and coupling patterns. Making use of residual dipolar couplings (RDCs), we were able to determine the relative configuration of a natural product containing seven stereocenters, including a chiral amine lacking direct RDC data. To identify the correct relative configuration out of 32 possible ones, experimental RDCs were used in three different approaches for data interpretation: by fitting experimental data based singular value decomposition (SVD) using a single alignment tensor and either (i) a single conformer or (ii) multiple conformers, or alternatively (iii) using molecular dynamics simulations with tensorial orientational constraints (MDOC). Even though in all three approaches one and the same configuration could be selected and clear discrimination between possible configurations was achieved, the experimental data was not fully satisfied by the methods based on single tensor approaches. While these two approaches are faster, only MDOC is able to fully reproduce experimental results, as the obtained conformational ensemble adequately covers the conformational space necessary to describe the molecule with inherent flexibility.

Efficient affinity ranking of fluorinated ligands by 19F NMR: CSAR and FastCSAR.

Fluorine NMR has recently gained high popularity in drug discovery as it allows efficient and sensitive screening of large numbers of ligands. However, the positive hits found in screening must subsequently be ranked according to their affinity in order to prioritize them for follow-up chemistry. Unfortunately, the primary read-out from the screening experiments, namely the increased relaxation rate upon binding, is not proportional to the affinity of the ligand, as it is polluted by effects such as exchange broadening. Here we present the method CSAR (Chemical Shift-anisotropy-based Affinity Ranking) for reliable ranking of fluorinated ligands by NMR, without the need of isotope labeled protein, titrations or setting up a reporter format. Our strategy is to produce relaxation data that is directly proportional to the binding affinity. This is achieved by removing all other contributions to relaxation as follows: (i) exchange effects are efficiently suppressed by using high power spin lock pulses, (ii) dipolar relaxation effects are approximately subtracted by measuring at two different magnetic fields and (iii) differences in chemical shift anisotropy are normalized using calculated values. A similar ranking can be obtained with the simplified approach FastCSAR that relies on a measurement of a single relaxation experiment at high field (preferably > 600 MHz). An affinity ranking obtained in this simple way will enable prioritizing ligands and thus improve the efficiency of fragment-based drug design.

Comprehensive and High-Throughput Exploration of Chemical Space Using Broadband 19 F†NMR-Based Screening.

Fragment-based lead discovery has become a fundamental approach to identify ligands that efficiently interact with disease-relevant targets. Among the numerous screening techniques, fluorine-detected NMR has gained popularity owing to its high sensitivity, robustness, and ease of use. To effectively explore chemical space, a universal NMR experiment, a rationally designed fragment library, and a sample composition optimized for a maximal number of compounds and minimal measurement time are required. Here, we introduce a comprehensive method that enabled the efficient assembly of a high-quality and diverse library containing nearly 4000 fragments and screening for target-specific binders within days. At the core of the approach is a novel broadband relaxation-edited NMR experiment that covers the entire chemical shift range of drug-like 19 F motifs in a single measurement. Our approach facilitates the identification of diverse binders and the fast ligandability assessment of new targets.

Structural basis and mechanism for metallochaperone-assisted assembly of the CuA center in cytochrome oxidase.

The mechanisms underlying the biogenesis of the structurally unique, binuclear Cu1.5+•Cu1.5+ redox center (CuA) on subunit II (CoxB) of cytochrome oxidases have been a long-standing mystery. Here, we reconstituted the CoxB•CuA center in vitro from apo-CoxB and the holo-forms of the copper transfer chaperones ScoI and PcuC. A previously unknown, highly stable ScoI•Cu2+•CoxB complex was shown to be rapidly formed as the first intermediate in the pathway. Moreover, our structural data revealed that PcuC has two copper-binding sites, one each for Cu1+ and Cu2+, and that only PcuC•Cu1+•Cu2+ can release CoxB•Cu2+ from the ScoI•Cu2+•CoxB complex. The CoxB•CuA center was then formed quantitatively by transfer of Cu1+ from a second equivalent of PcuC•Cu1+•Cu2+ to CoxB•Cu2+. This metalation pathway is consistent with all available in vivo data and identifies the sources of the Cu ions required for CuA center formation and the order of their delivery to CoxB.

Assessing molecular interactions with biophysical methods using the validation cross.

There are numerous methods for studying molecular interactions. However, each method gives rise to false negative- or false positive binding results, stemming from artifacts of the scientific equipment or from shortcomings of the experimental format. To validate an initial positive binding result, additional methods need to be applied to cover the shortcomings of the primary experiment. The aim of such a validation procedure is to exclude as many artifacts as possible to confirm that there is a true molecular interaction that meets the standards for publishing or is worth investing considerable resources for follow-up activities in a drug discovery project. To simplify this validation process, a graphical scheme - the validation cross - can be used. This simple graphic is a powerful tool for identifying blind spots of a binding hypothesis, for selecting the most informative combination of methods to reveal artifacts and, in general, for understanding more thoroughly the nature of a validation process. The concept of the validation cross was originally introduced for the validation of protein-ligand interactions by NMR in drug discovery. Here, an attempt is made to expand the concept to further biophysical methods and to generalize it for binary molecular interactions.

Inhibition of K-RAS4B by a Unique Mechanism of Action: Stabilizing Membrane-Dependent Occlusion of the Effector-Binding Site.

KRAS is frequently mutated in several of the most lethal types of cancer; however, the KRAS protein has proven a challenging drug target. K-RAS4B must be localized to the plasma membrane by prenylation to activate oncogenic signaling, thus we endeavored to target the protein-membrane interface with small-molecule compounds. While all reported lead compounds have low affinity for KRAS in solution, the potency of Cmpd2 was strongly enhanced when prenylated K-RAS4B is associated with a lipid bilayer. We have elucidated a unique mechanism of action of Cmpd2, which simultaneously engages a shallow pocket on KRAS and associates with the lipid bilayer, thereby stabilizing KRAS in an orientation in which the membrane occludes its effector-binding site, reducing RAF binding and impairing activation of RAF. Furthermore, enrichment of Cmpd2 on the bilayer enhances potency by promoting interaction with KRAS. This insight reveals a novel approach to developing inhibitors of membrane-associated proteins.

Cost-effective large-scale expression of proteins for NMR studies.

We present protocols for high-level expression of isotope-labelled proteins in E. coli in cost-effective ways. This includes production of large amounts of unlabeled proteins and 13C-methyl methionine labeling in rich media, where yields of up to a gram of soluble protein per liter of culture are reached. Procedures for uniform isotope labeling of 2H, 13C and 15N using auto-induction or isopropyl-ß-D-1-thiogalactopyranoside-induction are described, with primary focus on minimal isotope consumption and high reproducibility of protein expression. These protocols are based on high cell-density fermentation, but the key procedures are easily transferred to shake flask cultures.