Probing Protein-Ligand Methyl-π Interaction Geometries through Chemical Shift Measurements of Selectively Labeled Methyl Groups.

Fragment-based drug design is heavily dependent on the optimization of initial low-affinity binders. Herein we introduce an approach that uses selective labeling of methyl groups in leucine and isoleucine side chains to directly probe methyl-π contacts, one of the most prominent forms of interaction between proteins and small molecules. Using simple NMR chemical shift perturbation experiments with selected BRD4-BD1 binders, we find good agreement with a commonly used model of the ring-current effect as well as the overall interaction geometries extracted from the Protein Data Bank. By combining both interaction geometries and chemical shift calculations as fit quality criteria, we can position dummy aromatic rings into an AlphaFold model of the protein of interest. The proposed method can therefore provide medicinal chemists with important information about binding geometries of small molecules in fast and iterative matter, even in the absence of high-resolution experimental structures.

Ligand 1 H NMR Chemical Shifts as Accurate Reporters for Protein-Ligand Binding Interfaces in Solution.

The availability of high-resolution 3D structural information is crucial for investigating guest-host systems across a wide range of fields. In the context of drug discovery, the information is routinely used to establish and validate structure-activity relationships, grow initial hits from screening campaigns, and to guide molecular docking. For the generation of protein-ligand complex structural information, X-ray crystallography is the experimental method of choice, however, with limited information on protein flexibility. An experimentally verified structural model of the binding interface in the native solution-state would support medicinal chemists in their molecular design decisions. Here we demonstrate that protein-bound ligand 1 H NMR chemical shifts are highly sensitive and accurate probes for the immediate chemical environment of protein-ligand interfaces. By comparing the experimental ligand 1 H chemical shift values with those computed from the X-ray structure using quantum mechanics methodology, we identify significant disagreements for parts of the ligand between the two experimental techniques. We show that quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) ensembles can be used to refine initial X-ray co-crystal structures resulting in a better agreement with experimental 1 H ligand chemical shift values. Overall, our findings highlight the usefulness of ligand 1 H NMR chemical shift information in combination with a QM/MM MD workflow for generating protein-ligand ensembles that accurately reproduce solution structural data.

Hyperphosphorylation of Human Osteopontin and Its Impact on Structural Dynamics and Molecular Recognition.

Protein phosphorylation is an abundant post-translational modification (PTM) and an essential modulator of protein functionality in living cells. Intrinsically disordered proteins (IDPs) are particular targets of PTM protein kinases due to their involvement in fundamental protein interaction networks. Despite their dynamic nature, IDPs are far from having random-coil conformations but exhibit significant structural heterogeneity. Changes in the molecular environment, most prominently in the form of PTM via phosphorylation, can modulate these structural features. Therefore, how phosphorylation events can alter conformational ensembles of IDPs and their interactions with binding partners is of great interest. Here we study the effects of hyperphosphorylation on the IDP osteopontin (OPN), an extracellular target of the Fam20C kinase. We report a full characterization of the phosphorylation sites of OPN using a combined nuclear magnetic resonance/mass spectrometry approach and provide evidence for an increase in the local flexibility of highly phosphorylated regions and the ensuing overall structural elongation. Our study emphasizes the simultaneous importance of electrostatic and hydrophobic interactions in the formation of compact substates in IDPs and their relevance for molecular recognition events.

PI by NMR: Probing CH-π Interactions in Protein-Ligand Complexes by NMR Spectroscopy.

While CH-π interactions with target proteins are crucial determinants for the affinity of arguably every drug molecule, no method exists to directly measure the strength of individual CH-π interactions in drug-protein complexes. Herein, we present a fast and reliable methodology called PI (π interactions) by NMR, which can differentiate the strength of protein-ligand CH-π interactions in solution. By combining selective amino-acid side-chain labeling with 1 H-13 C NMR, we are able to identify specific protein protons of side-chains engaged in CH-π interactions with aromatic ring systems of a ligand, based solely on 1 H chemical-shift values of the interacting protein aromatic ring protons. The information encoded in the chemical shifts induced by such interactions serves as a proxy for the strength of each individual CH-π interaction. PI by NMR changes the paradigm by which chemists can optimize the potency of drug candidates: direct determination of individual π interactions rather than averaged measures of all interactions.

Modulation of Correlated Segment Fluctuations in IDPs upon Complex Formation as an Allosteric Regulatory Mechanism.

Molecular recognition of and by intrinsically disordered proteins (IDPs) is an intriguing and still largely elusive phenomenon. Typically, protein recognition involving IDPs requires either folding upon binding or, alternatively, the formation of "fuzzy complexes." Here we show via correlation analyses of paramagnetic relaxation enhancement data unprecedented and striking alterations of the concerted fluctuations within the conformational ensemble of IDPs upon ligand binding. We study the binding of α-synuclein to calmodulin, a ubiquitous calcium-binding protein, and the binding of the extracellular matrix IDP osteopontin to heparin, a mimic of the extracellular matrix ligand hyaluronic acid. In both cases, binding leads to reduction of correlated long-range motions in these two IDPs and thus indicates a loosening of structural compaction upon binding. Most importantly, however, the simultaneous presence of correlated and anti-correlated fluctuations in IDPs suggests the prevalence of "energetic frustration" and provides an explanation for the puzzling observation of disordered allostery in IDPs.

Late metabolic precursors for selective aromatic residue labeling.

In recent years, we developed a toolbox of heavy isotope containing compounds, which serve as metabolic amino acid precursors in the E. coli-based overexpression of aromatic residue labeled proteins. Our labeling techniques show excellent results both in terms of selectivity and isotope incorporation levels. They are additionally distinguished by low sample production costs and meet the economic demands to further implement protein NMR spectroscopy as a routinely used method in drug development processes. Different isotopologues allow for the assembly of optimized protein samples, which fulfill the requirements of various NMR experiments to elucidate protein structures, analyze conformational dynamics, or probe interaction surfaces. In the present article, we want to summarize the precursors we developed so far and give examples of their special value in the probing of protein-ligand interaction.

1H, 13C, 15N resonance assignment of human YAP 50-171 fragment.

Yes associated protein (YAP) is an intrinsically disordered protein that plays a major role in the Hippo pathway, regulating organ size, cell proliferation, apoptosis, and is associated with cancer development. Therefore, the binding between YAP and TEAD is an interesting target for cancer therapy. The TEAD binding domain of YAP was mapped to protein residues 50-171. To obtain further structural insights into this 12 kDa segment of YAP, we report a backbone and a partial sidechain assignment of recombinant YAP 50-171.

Highly Selective Stable Isotope Labeling of Histidine Residues by Using a Novel Precursor in E.†coli-Based Overexpression Systems.

The importance of NMR spectroscopy in unraveling the structural and dynamic properties of proteins is ever-expanding owing to progress in experimental techniques, hardware development, and novel labeling approaches. Multiple sophisticated methods of aliphatic residue labeling can be found in the literature, whereas the selective incorporation of NMR active isotopes into other amino acids still holds the potential for improvement. In order to close this methodological gap, we present a novel metabolic precursor for cell-based protein overexpression to assemble 13 C/2 H isotope patterns in the peptide backbone, as well as in side chain positions of a mechanistically distinguished histidine residue.

NMR probing and visualization of correlated structural fluctuations in intrinsically disordered proteins.

A novel statistical analysis of paramagnetic relaxation enhancement (PRE) and paramagnetic relaxation interference (PRI) based nuclear magnetic resonance (NMR) data is proposed based on the computation of correlation matrices. The technique is demonstrated with an example of the intrinsically disordered proteins (IDPs) osteopontin (OPN) and brain acid soluble protein 1 (BASP1). The correlation analysis visualizes in detail the subtleties of conformational averaging in IDPs and highlights the presence of correlated structural fluctuations of individual sub-domains in IDPs.

Arginine: Its pKa value revisited.

Using complementary approaches of potentiometry and NMR spectroscopy, we have determined that the equilibrium acid dissociation constant (pKa value) of the arginine guanidinium group is 13.8 ± 0.1. This is substantially higher than that of ∼ 12 often used in structure-based electrostatics calculations and cited in biochemistry textbooks. The revised intrinsic pKa value helps explains why arginine side chains in proteins are always predominantly charged, even at pH values as great as 10. The high pKa value also reinforces the observation that arginine side chains are invariably protonated under physiological conditions of near neutral pH. This occurs even when the guanidinium moiety is buried in a hydrophobic micro-environment, such as that inside a protein or a lipid membrane, thought to be incompatible with the presence of a charged group.

(1)H, (15)N, (13)C resonance assignment of human osteopontin.

Osteopontin (OPN) is a 33.7 kDa intrinsically disordered protein and a member of the SIBLING family of proteins. OPN is bearing a signal peptide for secretion into the extracellular space, where it exerts its main physiological function, the control of calcium biomineralization. It is often involved in tumorigenic processes influencing proliferation, migration and survival, as well as the adhesive properties of cancer cells via CD44 and integrin signaling pathways. Here we report the nearly complete NMR chemical shift assignment of recombinant human osteopontin.

pH-dependent random coil (1)H, (13)C, and (15)N chemical shifts of the ionizable amino acids: a guide for protein pK a measurements.

The pK a values and charge states of ionizable residues in polypeptides and proteins are frequently determined via NMR-monitored pH titrations. To aid the interpretation of the resulting titration data, we have measured the pH-dependent chemical shifts of nearly all the (1)H, (13)C, and (15)N nuclei in the seven common ionizable amino acids (X = Asp, Glu, His, Cys, Tyr, Lys, and Arg) within the context of a blocked tripeptide, acetyl-Gly-X-Gly-amide. Alanine amide and N-acetyl alanine were used as models of the N- and C-termini, respectively. Together, this study provides an essentially complete set of pH-dependent intra-residue and nearest-neighbor reference chemical shifts to help guide protein pK a measurements. These data should also facilitate pH-dependent corrections in algorithms used to predict the chemical shifts of random coil polypeptides. In parallel, deuterium isotope shifts for the side chain (15)N nuclei of His, Lys, and Arg in their positively-charged and neutral states were also measured. Along with previously published results for Asp, Glu, Cys, and Tyr, these deuterium isotope shifts can provide complementary experimental evidence for defining the ionization states of protein residues.

Compensatory adaptations of structural dynamics in an intrinsically disordered protein complex.

Intrinsically disordered proteins (IDPs) play crucial roles in protein interaction networks and in this context frequently constitute important hubs and interfaces. Here we show by a combination of NMR and EPR spectroscopy that the binding of the cytokine osteopontin (OPN) to its natural ligand, heparin, is accompanied by thermodynamically compensating structural adaptations. The core segment of OPN expands upon binding. This "unfolding-upon-binding" is governed primarily through electrostatic interactions between heparin and charged patches along the protein backbone and compensates for entropic penalties due to heparin-OPN binding. It is shown how structural unfolding compensates for entropic losses through ligand binding in IDPs and elucidates the interplay between structure and thermodynamics of rapid substrate-binding and -release events in IDP interaction networks.

Cooperative unfolding of compact conformations of the intrinsically disordered protein osteopontin.

Intrinsically disordered proteins (IDPs) constitute a class of biologically active proteins that lack defined tertiary and often secondary structure. The IDP Osteopontin (OPN), a cytokine involved in metastasis of several types of cancer, is shown to simultaneously sample extended, random coil-like conformations and stable, cooperatively folded conformations. By a combination of two magnetic resonance methods, electron paramagnetic resonance and nuclear magnetic resonance spectroscopy, we demonstrate that the OPN ensemble exhibits not only characteristics of an extended and flexible polypeptide, as expected for an IDP, but also simultaneously those of globular proteins, in particular sigmoidal structural denaturation profiles. Both types of states, extended and cooperatively folded, are populated simultaneously by OPN in its apo state. The heterogeneity of the structural properties of IDPs is thus shown to even involve cooperative folding and unfolding events.

The metastasis-associated extracellular matrix protein osteopontin forms transient structure in ligand interaction sites.

Osteopontin (OPN) is an acidic hydrophilic glycophosphoprotein that was first identified as a major sialoprotein in bones. It functions as a cell attachment protein displaying a RGD cell adhesion sequence and as a cytokine that signals through integrin and CD44 cell adhesion molecules. OPN is also implicated in human tumor progression and cell invasion. OPN has intrinsic transforming activity, and elevated OPN levels promote metastasis. OPN gene expression is also strongly activated in avian fibroblasts simultaneously transformed by the v-myc and v-mil(raf) oncogenes. Here we have investigated the solution structure of a 220-amino acid recombinant OPN protein by an integrated structural biology approach employing bioinformatic sequence analysis, multidimensional nuclear magnetic resonance spectroscopy, synchrotron radiation circular dichroism spectroscopy, and small-angle X-ray scattering. These studies suggest that OPN is an intrinsically unstructured protein in solution. Although OPN does not fold into a single defined structure, its conformational flexibility significantly deviates from random coil-like behavior. OPN comprises distinct local secondary structure elements with reduced conformational flexibility and substantially populates a compact subspace displaying distinct tertiary contacts. These compacted regions of OPN encompass the binding sites for α(V)ß(III) integrin and heparin. The conformational flexibility combined with the modular architecture of OPN may represent an important structural prerequisite for its functional diversity.