Quantification of natural abundance NMR data differentiates the solution behavior of monoclonal antibodies and their fragments.

Biotherapeutics are an important class of molecules for the treatment of a wide range of diseases. They include low molecular weight peptides, highly engineered protein scaffolds and monoclonal antibodies. During their discovery and development, assessments of the biophysical attributes is critical to understanding the solution behavior of therapeutic proteins and for de-risking liabilities. Thus, methods that can quantify, characterize, and provide a basis to inform risks and drive the selection of more optimal antibody and alternative scaffolds are needed. Nuclear Magnetic Resonance (NMR) spectroscopy is a technique that provides a means to probe antibody and antibody-like molecules in solution, at atomic resolution, under any formulated conditions. Here, all samples were profiled at natural abundance requiring no isotope enrichment. We present a numerical approach that quantitates two-dimensional methyl spectra. The approach was tested with a reference dataset that contained different types of antibody and antibody-like molecules. This dataset was processed through a procedure we call a Random Sampling of NMR Peaks for Covariance Analysis. This analysis revealed that the first two components were well correlated with the hydrodynamic radius of the molecules included in the reference set. Higher-order principal components were also linked to dynamic features between different tethered antibody-like molecules and contributed to decisions around candidate selection. The reference set provides a basis to characterize molecules with unknown solution behavior and is sensitive to the behavior of a molecule formulated under different conditions. The approach is independent of protein design, scaffold, formulation and provides a facile method to quantify solution behavior.

EGCG binds intrinsically disordered N-terminal domain of p53 and disrupts p53-MDM2 interaction.

Epigallocatechin gallate (EGCG) from green tea can induce apoptosis in cancerous cells, but the underlying molecular mechanisms remain poorly understood. Using SPR and NMR, here we report a direct, µM interaction between EGCG and the tumor suppressor p53 (KD = 1.6 ± 1.4 µM), with the disordered N-terminal domain (NTD) identified as the major binding site (KD = 4 ± 2 µM). Large scale atomistic simulations (>100 µs), SAXS and AUC demonstrate that EGCG-NTD interaction is dynamic and EGCG causes the emergence of a subpopulation of compact bound conformations. The EGCG-p53 interaction disrupts p53 interaction with its regulatory E3 ligase MDM2 and inhibits ubiquitination of p53 by MDM2 in an in vitro ubiquitination assay, likely stabilizing p53 for anti-tumor activity. Our work provides insights into the mechanisms for EGCG's anticancer activity and identifies p53 NTD as a target for cancer drug discovery through dynamic interactions with small molecules.

Solution structure and functional investigation of human guanylate kinase reveals allosteric networking and a crucial role for the enzyme in cancer.

Human guanylate kinase (hGMPK) is the only known enzyme responsible for cellular GDP production, making it essential for cellular viability and proliferation. Moreover, hGMPK has been assigned a critical role in metabolic activation of antiviral and antineoplastic nucleoside-analog prodrugs. Given that hGMPK is indispensable for producing the nucleotide building blocks of DNA, RNA, and cGMP and that cancer cells possess elevated GTP levels, it is surprising that a detailed structural and functional characterization of hGMPK is lacking. Here, we present the first high-resolution structure of hGMPK in the apo form, determined with NMR spectroscopy. The structure revealed that hGMPK consists of three distinct regions designated as the LID, GMP-binding (GMP-BD), and CORE domains and is in an open configuration that is nucleotide binding-competent. We also demonstrate that nonsynonymous single-nucleotide variants (nsSNVs) of the hGMPK CORE domain distant from the nucleotide-binding site of this domain modulate enzymatic activity without significantly affecting hGMPK's structure. Finally, we show that knocking down the hGMPK gene in lung adenocarcinoma cell lines decreases cellular viability, proliferation, and clonogenic potential while not altering the proliferation of immortalized, noncancerous human peripheral airway cells. Taken together, our results provide an important step toward establishing hGMPK as a potential biomolecular target, from both an orthosteric (ligand-binding sites) and allosteric (location of CORE domain-located nsSNVs) standpoint.

1H, 13C and 15N resonance assignment of human guanylate kinase.

Human guanylate kinase (hGMPK) is a critical enzyme that, in addition to phosphorylating its physiological substrate (d)GMP, catalyzes the second phosphorylation step in the conversion of anti-viral and anti-cancer nucleoside analogs to their corresponding active nucleoside analog triphosphates. Until now, a high-resolution structure of hGMPK is unavailable and thus, we studied free hGMPK by NMR and assigned the chemical shift resonances of backbone and side chain 1H, 13C, and 15N nuclei as a first step towards the enzyme's structural and mechanistic analysis with atomic resolution.

Solution NMR views of dynamical ordering of biomacromolecules.

BACKGROUND: To understand the mechanisms related to the 'dynamical ordering' of macromolecules and biological systems, it is crucial to monitor, in detail, molecular interactions and their dynamics across multiple timescales. Solution nuclear magnetic resonance (NMR) spectroscopy is an ideal tool that can investigate biophysical events at the atomic level, in near-physiological buffer solutions, or even inside cells. SCOPE OF REVIEW: In the past several decades, progress in solution NMR has significantly contributed to the elucidation of three-dimensional structures, the understanding of conformational motions, and the underlying thermodynamic and kinetic properties of biomacromolecules. This review discusses recent methodological development of NMR, their applications and some of the remaining challenges. MAJOR CONCLUSIONS: Although a major drawback of NMR is its difficulty in studying the dynamical ordering of larger biomolecular systems, current technologies have achieved considerable success in the structural analysis of substantially large proteins and biomolecular complexes over 1MDa and have characterised a wide range of timescales across which biomolecular motion exists. While NMR is well suited to obtain local structure information in detail, it contributes valuable and unique information within hybrid approaches that combine complementary methodologies, including solution scattering and microscopic techniques. GENERAL SIGNIFICANCE: For living systems, the dynamic assembly and disassembly of macromolecular complexes is of utmost importance for cellular homeostasis and, if dysregulated, implied in human disease. It is thus instructive for the advancement of the study of the dynamical ordering to discuss the potential possibilities of solution NMR spectroscopy and its applications. This article is part of a Special Issue entitled "Biophysical Exploration of Dynamical Ordering of Biomolecular Systems" edited by Dr. Koichi Kato.

Simultaneous determination of fast and slow dynamics in molecules using extreme CPMG relaxation dispersion experiments.

Molecular dynamics play a significant role in how molecules perform their function. A critical method that provides information on dynamics, at the atomic level, is NMR-based relaxation dispersion (RD) experiments. RD experiments have been utilized for understanding multiple biological processes occurring at micro-to-millisecond time, such as enzyme catalysis, molecular recognition, ligand binding and protein folding. Here, we applied the recently developed high-power RD concept to the Carr-Purcell-Meiboom-Gill sequence (extreme CPMG; E-CPMG) for the simultaneous detection of fast and slow dynamics. Using a fast folding protein, gpW, we have shown that previously inaccessible kinetics can be accessed with the improved precision and efficiency of the measurement by using this experiment.

Deciphering Conformational Changes Associated with the Maturation of Thrombin Anion Binding Exosite I.

Thrombin participates in procoagulation, anticoagulation, and platelet activation. This enzyme contains anion binding exosites, ABE I and ABE II, which attract regulatory biomolecules. As prothrombin is activated to thrombin, pro-ABE I is converted into mature ABE I. Unexpectedly, certain ligands can bind to pro-ABE I specifically. Moreover, knowledge of changes in conformation and affinity that occur at the individual residue level as pro-ABE I is converted to ABE I is lacking. Such changes are transient and were not captured by crystallography. Therefore, we employed nuclear magnetic resonance (NMR) titrations to monitor development of ABE I using peptides based on protease-activated receptor 3 (PAR3). Proton line broadening NMR revealed that PAR3 (44-56) and more weakly binding PAR3G (44-56) could already interact with pro-ABE I on prothrombin. 1H-15N heteronuclear single-quantum coherence NMR titrations were then used to probe binding of individual 15N-labeled PAR3G residues (F47, E48, L52, and D54). PAR3G E48 and D54 could interact electrostatically with prothrombin and tightened upon thrombin maturation. The higher affinity for PAR3G D54 suggests the region surrounding thrombin R77a is better oriented to bind D54 than the interaction between PAR3G E48 and thrombin R75. Aromatic PAR3G F47 and aliphatic L52 both reported on significant changes in the chemical environment upon conversion of prothrombin to thrombin. The ABE I region surrounding the 30s loop was more affected than the hydrophobic pocket (F34, L65, and I82). Our NMR titrations demonstrate that PAR3 residues document structural rearrangements occurring during exosite maturation that are missed by reported X-ray crystal structures.

A Small Molecule Causes a Population Shift in the Conformational Landscape of an Intrinsically Disordered Protein.

Intrinsically disordered proteins (IDPs) have roles in myriad biological processes and numerous human diseases. However, kinetic and amplitude information regarding their ground-state conformational fluctuations has remained elusive. We demonstrate using nuclear magnetic resonance (NMR)-based relaxation dispersion that the D2 domain of p27Kip1, a prototypical IDP, samples multiple discrete, rapidly exchanging conformational states. By combining NMR with mutagenesis and small-angle X-ray scattering (SAXS), we show that these states involve aromatic residue clustering through long-range hydrophobic interactions. Theoretical studies have proposed that small molecules bind promiscuously to IDPs, causing expansion of their conformational landscapes. However, on the basis of previous NMR-based screening results, we show here that compound binding only shifts the populations of states that existed within the ground state of apo p27-D2 without changing the barriers between states. Our results provide atomic resolution insight into how a small molecule binds an IDP and emphasize the need to examine motions on the low microsecond time scale when probing these types of interactions.

Recent advances in measuring the kinetics of biomolecules by NMR relaxation dispersion spectroscopy.

Protein function can be modulated or dictated by the amplitude and timescale of biomolecular motion, therefore it is imperative to study protein dynamics. Nuclear Magnetic Resonance (NMR) spectroscopy is a powerful technique capable of studying timescales of motion that range from those faster than molecular reorientation on the picosecond timescale to those that occur in real-time. Across this entire regime, NMR observables can report on the amplitude of atomic motion, and the kinetics of atomic motion can be ascertained with a wide variety of experimental techniques from real-time to milliseconds and several nanoseconds to picoseconds. Still a four orders of magnitude window between several nanoseconds and tens of microseconds has remained elusive. Here, we highlight new relaxation dispersion NMR techniques that serve to cover this "hidden-time" window up to hundreds of nanoseconds that achieve atomic resolution while studying the molecule under physiological conditions.

Kinetics of the Antibody Recognition Site in the Third IgG-Binding Domain of Protein G.

Protein dynamics occurring on a wide range of timescales play a crucial role in governing protein function. Particularly, motions between the globular rotational correlation time (τc ) and 40 µs (supra-τc window), strongly influence molecular recognition. This supra-τc window was previously hidden, owing to a lack of experimental methods. Recently, we have developed a high-power relaxation dispersion (RD) experiment for measuring kinetics as fast as 4 µs. For the first time, this method, performed under super-cooled conditions, enabled us to detect a global motion in the first ß-turn of the third IgG-binding domain of protein G (GB3), which was extrapolated to 371±115 ns at 310 K. Furthermore, the same residues show the plasticity in the model-free residual dipolar coupling (RDC) order parameters and in an ensemble encoding the supra-τc dynamics. This ß-turn is involved in antibody binding, exhibiting the potential link of the observed supra-τc motion with molecular recognition.

High-power (1)H composite pulse decoupling provides artifact free exchange-mediated saturation transfer (EST) experiments.

Exchange-mediated saturation transfer (EST) provides critical information regarding dynamics of molecules. In typical applications EST is studied by either scanning a wide range of (15)N chemical shift offsets where the applied (15)N irradiation field strength is on the order of hundreds of Hertz or, scanning a narrow range of (15)N chemical shift offsets where the applied (15)N irradiation field-strength is on the order of tens of Hertz during the EST period. The (1)H decoupling during the EST delay is critical as incomplete decoupling causes broadening of the EST profile, which could possibly result in inaccuracies of the extracted kinetic parameters and transverse relaxation rates. Currently two different (1)H decoupling schemes have been employed, intermittently applied 180° pulses and composite-pulse-decoupling (CPD), for situations where a wide range, or narrow range of (15)N chemical shift offsets are scanned, respectively. We show that high-power CPD provides artifact free EST experiments, which can be universally implemented regardless of the offset range or irradiation field-strengths.

Allosteric switch regulates protein-protein binding through collective motion.

Many biological processes depend on allosteric communication between different parts of a protein, but the role of internal protein motion in propagating signals through the structure remains largely unknown. Through an experimental and computational analysis of the ground state dynamics in ubiquitin, we identify a collective global motion that is specifically linked to a conformational switch distant from the binding interface. This allosteric coupling is also present in crystal structures and is found to facilitate multispecificity, particularly binding to the ubiquitin-specific protease (USP) family of deubiquitinases. The collective motion that enables this allosteric communication does not affect binding through localized changes but, instead, depends on expansion and contraction of the entire protein domain. The characterization of these collective motions represents a promising avenue for finding and manipulating allosteric networks.

Nucleophosmin integrates within the nucleolus via multi-modal interactions with proteins displaying R-rich linear motifs and rRNA.

The nucleolus is a membrane-less organelle formed through liquid-liquid phase separation of its components from the surrounding nucleoplasm. Here, we show that nucleophosmin (NPM1) integrates within the nucleolus via a multi-modal mechanism involving multivalent interactions with proteins containing arginine-rich linear motifs (R-motifs) and ribosomal RNA (rRNA). Importantly, these R-motifs are found in canonical nucleolar localization signals. Based on a novel combination of biophysical approaches, we propose a model for the molecular organization within liquid-like droplets formed by the N-terminal domain of NPM1 and R-motif peptides, thus providing insights into the structural organization of the nucleolus. We identify multivalency of acidic tracts and folded nucleic acid binding domains, mediated by N-terminal domain oligomerization, as structural features required for phase separation of NPM1 with other nucleolar components in vitro and for localization within mammalian nucleoli. We propose that one mechanism of nucleolar localization involves phase separation of proteins within the nucleolus.

Monitoring Ligand-Induced Protein Ordering in Drug Discovery.

While the gene for p53 is mutated in many human cancers causing loss of function, many others maintain a wild-type gene but exhibit reduced p53 tumor suppressor activity through overexpression of the negative regulators, Mdm2 and/or MdmX. For the latter mechanism of loss of function, the activity of endogenous p53 can be restored through inhibition of Mdm2 or MdmX with small molecules. We previously reported a series of compounds based upon the Nutlin-3 chemical scaffold that bind to both MdmX and Mdm2 [Vara, B. A. et al. (2014) Organocatalytic, diastereo- and enantioselective synthesis of nonsymmetric cis-stilbene diamines: A platform for the preparation of single-enantiomer cis-imidazolines for protein-protein inhibition. J. Org. Chem. 79, 6913-6938]. Here we present the first solution structures based on data from NMR spectroscopy for MdmX in complex with four of these compounds and compare them with the MdmX:p53 complex. A p53-derived peptide binds with high affinity (Kd value of 150nM) and causes the formation of an extensive network of hydrogen bonds within MdmX; this constitutes the induction of order within MdmX through ligand binding. In contrast, the compounds bind more weakly (Kd values from 600nM to 12µM) and induce an incomplete hydrogen bond network within MdmX. Despite relatively weak binding, the four compounds activated p53 and induced p21(Cip1) expression in retinoblastoma cell lines that overexpress MdmX, suggesting that they specifically target MdmX and/or Mdm2. Our results document structure-activity relationships for lead-like small molecules targeting MdmX and suggest a strategy for their further optimization in the future by using NMR spectroscopy to monitor small-molecule-induced protein order as manifested through hydrogen bond formation.

Discovery of Small Molecules that Inhibit the Disordered Protein, p27(Kip1).

Disordered proteins are highly prevalent in biological systems, they control myriad signaling and regulatory processes, and their levels and/or cellular localization are often altered in human disease. In contrast to folded proteins, disordered proteins, due to conformational heterogeneity and dynamics, are not considered viable drug targets. We challenged this paradigm by identifying through NMR-based screening small molecules that bound specifically, albeit weakly, to the disordered cell cycle regulator, p27(Kip1) (p27). Two groups of molecules bound to sites created by transient clusters of aromatic residues within p27. Conserved chemical features within these two groups of small molecules exhibited complementarity to their binding sites within p27, establishing structure-activity relationships for small molecule:disordered protein interactions. Finally, one compound counteracted the Cdk2/cyclin A inhibitory function of p27 in vitro, providing proof-of-principle that small molecules can inhibit the function of a disordered protein (p27) through sequestration in a conformation incapable of folding and binding to a natural regulatory target (Cdk2/cyclin A).

Sampling of Glycan-Bound Conformers by the Anti-HIV Lectin Oscillatoria agardhii agglutinin in the Absence of Sugar.

Lectins from different sources have been shown to interfere with HIV infection by binding to the sugars of viral-envelope glycoproteins. Three-dimensional atomic structures of a number of HIV-inactivating lectins have been determined, both as free proteins and in glycan-bound forms. However, details on the mechanism of recognition and binding to sugars are elusive. Herein we focus on the anti-HIV lectin OAA from Oscillatoria agardhii: We show that in the absence of sugars in solution, both the sugar-free and sugar-bound protein conformations that were observed in the X-ray crystal structures exist as conformational substates. Our results suggest that glycan recognition occurs by conformational selection within the ground state; this model differs from the popular "excited-state" model. Our findings provide further insight into molecular recognition of the major receptor on the HIV virus by OAA. These details can potentially be used for the optimization and/or development of preventive anti-HIV therapeutics.

Population shuffling of protein conformations.

Motions play a vital role in the functions of many proteins. Discrete conformational transitions to excited states, happening on timescales of hundreds of microseconds, have been extensively characterized. On the other hand, the dynamics of the ground state are widely unexplored. Newly developed high-power relaxation dispersion experiments allow the detection of motions up to a one-digit microsecond timescale. These experiments showed that side chains in the hydrophobic core as well as at protein-protein interaction surfaces of both ubiquitin and the third immunoglobulin binding domain of protein G move on the microsecond timescale. Both proteins exhibit plasticity to this microsecond motion through redistribution of the populations of their side-chain rotamers, which interconvert on the picosecond to nanosecond timescale, making it likely that this "population shuffling" process is a general mechanism.

A designed conformational shift to control protein binding specificity.

In a conformational selection scenario, manipulating the populations of binding-competent states should be expected to affect protein binding. We demonstrate how in silico designed point mutations within the core of ubiquitin, remote from the binding interface, change the binding specificity by shifting the conformational equilibrium of the ground-state ensemble between open and closed substates that have a similar population in the wild-type protein. Binding affinities determined by NMR titration experiments agree with the predictions, thereby showing that, indeed, a shift in the conformational equilibrium enables us to alter ubiquitin's binding specificity and hence its function. Thus, we present a novel route towards designing specific binding by a conformational shift through exploiting the fact that conformational selection depends on the concentration of binding-competent substates.

ORIUM: optimized RDC-based Iterative and Unified Model-free analysis.

Residual dipolar couplings (RDCs) are NMR parameters that provide both structural and dynamic information concerning inter-nuclear vectors, such as N-H(N) and Cα-Hα bonds within the protein backbone. Two approaches for extracting this information from RDCs are the model free analysis (MFA) (Meiler et al. in J Am Chem Soc 123:6098-6107, 2001; Peti et al. in J Am Chem Soc 124:5822-5833, 2002) and the direct interpretation of dipolar couplings (DIDCs) (Tolman in J Am Chem Soc 124:12020-12030, 2002). Both methods have been incorporated into iterative schemes, namely the self-consistent RDC based MFA (SCRM) (Lakomek et al. in J Biomol NMR 41:139-155, 2008) and iterative DIDC (Yao et al. in J Phys Chem B 112:6045-6056, 2008), with the goal of removing the influence of structural noise in the MFA and DIDC formulations. Here, we report a new iterative procedure entitled Optimized RDC-based Iterative and Unified Model-free analysis (ORIUM). ORIUM unifies theoretical concepts developed in the MFA, SCRM, and DIDC methods to construct a computationally less demanding approach to determine these structural and dynamic parameters. In all schemes, dynamic averaging reduces the actual magnitude of the alignment tensors complicating the determination of the absolute values for the generalized order parameters. To readdress this scaling issue that has been previously investigated (Lakomek et al. in J Biomol NMR 41:139-155, 2008; Salmon et al. in Angew Chem Int Edit 48:4154-4157, 2009), a new method is presented using only RDC data to establish a lower bound on protein motion, bypassing the requirement of Lipari-Szabo order parameters. ORIUM and the new scaling procedure are applied to the proteins ubiquitin and the third immunoglobulin domain of protein G (GB3). Our results indicate good agreement with the SCRM and iterative DIDC approaches and signify the general applicability of ORIUM and the proposed scaling for the extraction of inter-nuclear vector structural and dynamic content.

Measuring dynamic and kinetic information in the previously inaccessible supra-τ(c) window of nanoseconds to microseconds by solution NMR spectroscopy.

Nuclear Magnetic Resonance (NMR) spectroscopy is a powerful tool that has enabled experimentalists to characterize molecular dynamics and kinetics spanning a wide range of time-scales from picoseconds to days. This review focuses on addressing the previously inaccessible supra-tc window (defined as τ(c) < supra-τ(c) < 40 µs; in which tc is the overall tumbling time of a molecule) from the perspective of local inter-nuclear vector dynamics extracted from residual dipolar couplings (RDCs) and from the perspective of conformational exchange captured by relaxation dispersion measurements (RD). The goal of the first section is to present a detailed analysis of how to extract protein dynamics encoded in RDCs and how to relate this information to protein functionality within the previously inaccessible supra-τ(c) window. In the second section, the current state of the art for RD is analyzed, as well as the considerable progress toward pushing the sensitivity of RD further into the supra-τ(c) scale by up to a factor of two (motion up to 25 µs). From the data obtained with these techniques and methodology, the importance of the supra-τ(c) scale for protein function and molecular recognition is becoming increasingly clearer as the connection between motion on the supra-τ(c) scale and protein functionality from the experimental side is further strengthened with results from molecular dynamics simulations.