Exploring Multi-Subsite Binding Pockets in Proteins: DEEP-STD NMR Fingerprinting and Molecular Dynamics Unveil a Cryptic Subsite at the GM1 Binding Pocket of Cholera Toxin B

Ligand-based NMR techniques to study protein–ligand interactions are potent tools in drug design. Saturation transfer difference (STD) NMR spectroscopy stands out as one of the most versatile techniques, allowing screening of fragments libraries and providing structural information on binding modes. Recently, it has been shown that a multi-frequency STD NMR approach, differential epitope mapping (DEEP)-STD NMR, can provide additional information on the orientation of small ligands within the binding pocket. Here, the approach is extended to a so-called DEEP-STD NMR fingerprinting technique to explore the binding subsites of cholera toxin subunit B (CTB). To that aim, the synthesis of a set of new ligands is presented, which have been subject to a thorough study of their interactions with CTB by weak affinity chromatography (WAC) and NMR spectroscopy. Remarkably, the combination of DEEP-STD NMR fingerprinting and Hamiltonian replica exchange molecular dynamics has proved to be an excellent approach to explore the geometry, flexibility, and ligand occupancy of multi-subsite binding pockets. In the particular case of CTB, it allowed the existence of a hitherto unknown binding subsite adjacent to the GM1 binding pocket to be revealed, paving the way to the design of novel leads for inhibition of this relevant toxin.     

Probing the binding of propranolol enantiomers to alpha1-acid glycoprotein with ligand-detected NMR experiments

Mapping the interactions of a small molecule ligand with a protein can provide information important for biochemical studies and for drug design and development. This information can be determined using the ligand-detected (1)H NMR experiments T(1rho)-NOESY, diffusion, and saturation transfer difference (STD). This work compares the results of these experiments and examines their ability to distinguish the binding epitopes of propranolol enantiomers with alpha 1-acid glycoprotein (AGP). The epitope maps for the propranolol enantiomers are fairly similar, as expected from their similar binding affinities; however, the STD epitope maps provide unique insights into the different orientations of the enantiomers with respect to the AGP binding pocket. Our results suggest that it is best to consider the data provided by several NMR epitope mapping experiments in drawing conclusions about ligand-protein binding interactions.     

Investigation of the binding of a carbohydrate-mimetic peptide to its complementary anticarbohydrate antibody by STD-NMR spectroscopy and molecular-dynamics simulations

Saturation transfer difference (STD)‐NMR spectroscopy was used to probe experimentally the bioactive solution conformation of the carbohydrate mimic MDWNMHAA 1 of the O‐polysaccharide of Shigella flexneri Y when bound to its complementary antibody, mAb SYA/J6. Molecular dynamics simulations using the ZymeCAD? Molecular Dynamics platform were also undertaken to give a more accurate picture of the conformational flexibility and the possibilities for bound ligand conformations. The ligand topology, or the dynamic epitope, was mapped with the CORCEMA‐ST (COmplete Relaxation and Conformational Exchange Matrix Analysis of Saturation Transfer) program that calculates a total matrix analysis of relaxation and exchange effects to generate predicted STD‐NMR intensities from simulation. The comparison of these predicted STD enhancements with experimental data was used to select a representative binding mode. A protocol that employed theoretical STD effects calculated at snapshots during the entire course of a molecular dynamics (MD) trajectory of the peptide bound to the Fv portion of the antibody, and not the averaged atomic positions of receptor–ligand complexes, was also examined. In addition, the R factor was calculated on the basis of STD (fit) to avoid T1 bias, and an effective R factor, R_eff, was defined such that if the calculated STD (fit) for proton k was within error of the experimental STD (fit) for proton k, then that calculated STD (fit) for proton k was not included in the calculation of the R factor. This protocol was effective in deriving the antibody‐bound solution conformation of the peptide which also differed from the bound conformation determined by X‐ray crystallography; however, several discrepancies between experimental and calculated STD (fit) values were observed. The bound conformation was therefore further refined with a simulated annealing refinement protocol known as STD‐NMR intensity‐restrained CORCEMA optimization (SICO) to give a more accurate representation of the bound peptide epitope. Further optimization was required in this case, but a satisfactory correlation between experimental and calculated STD values was obtained. Attempts were also made to obtain STD enhancements with a synthetic pentasaccharide hapten, corresponding to the O‐polysaccharide, while bound to the antibody. However, unfavorable kinetics of binding in this system prevented sufficient STD build‐up, which, in turn, hindered a rigorous analysis via full STD build‐up curves.     

Ligand–Receptor Binding Affinities from Saturation Transfer Difference (STD) NMR Spectroscopy: The Binding Isotherm of STD Initial Growth Rates

The direct evaluation of dissociation constants (K_D) from the variation of saturation transfer difference (STD) NMR spectroscopy values with the receptor–ligand ratio is not feasible due to the complex dependence of STD intensities on the spectral properties of the observed signals. Indirect evaluation, by competition experiments, allows the determination of K_D, as long as a ligand of known affinity is available for the protein under study. Herein, we present a novel protocol based on STD NMR spectroscopy for the direct measurements of receptor–ligand dissociation constants (K_D) from single‐ligand titration experiments. The influence of several experimental factors on STD values has been studied in detail, confirming the marked impact on standard determinations of protein–ligand affinities by STD NMR spectroscopy. These factors, namely, STD saturation time, ligand residence time in the complex, and the intensity of the signal, affect the accumulation of saturation in the free ligand by processes closely related to fast protein–ligand rebinding and longitudinal relaxation of the ligand signals. The proposed method avoids the dependence of the magnitudes of ligand STD signals at a given saturation time on spurious factors by constructing the binding isotherms using the initial growth rates of the STD amplification factors, in a similar way to the use of NOE growing rates to estimate cross relaxation rates for distance evaluations. Herein, it is demonstrated that the effects of these factors are cancelled out by analyzing the protein–ligand association curve using STD values at the limit of zero saturation time, when virtually no ligand rebinding or relaxation takes place. The approach is validated for two well‐studied protein–ligand systems: the binding of the saccharides GlcNAc and GlcNAcβ1,4GlcNAc (chitobiose) to the wheat germ agglutinin (WGA) lectin, and the interaction of the amino acid L ‐tryptophan to bovine serum albumin (BSA). In all cases, the experimental K_D measured under different experimental conditions converged to the thermodynamic values. The proposed protocol allows accurate determinations of protein–ligand dissociation constants, extending the applicability of the STD NMR spectroscopy for affinity measurements, which is of particular relevance for those proteins for which a ligand of known affinity is not available.     

Refinement of the conformation of UDP-galactose bound to galactosyltransferase using the STD NMR intensity-restrained CORCEMA optimization

The STD NMR technique has originally been described as a tool for screening large compound libraries to identify the lead compounds that are specific to target proteins of interest. The application of this technique in the qualitative epitope mapping of ligands weakly binding to proteins, virus capsid shells, and nucleic acids has also been described. Here we describe the application of the STD NMR intensity-restrained CORCEMA optimization (SICO) procedure for refining the bound conformation of UDP-galactose in galactosyltransferase complex using STD NMR intensities recorded at 500 MHz as the experimental constraints. A comparison of the SICO structure for the bound UDP-galactose in solution with that in the crystal structure for this complex shows some differences in ligand torsion angles and V253 side-chain orientation in the protein. This work describes the first application of an STD NMR intensity-restrained CORCEMA optimization procedure for refining the torsion angles of a bound ligand structure. This method is likely to be useful in structure-based drug design programs since most initial lead compounds generally exhibit weak affinity (millimolar to micromolar) to target proteins of pharmaceutical interest, and the bound conformation of these lead compounds in the protein binding pocket can be determined by the CORCEMA-ST refinement.     

Unwinding of the C‐Terminal Residues of Neuropeptide Y is critical for Y2 Receptor Binding and Activation

Despite recent breakthroughs in the structural characterization of G‐protein‐coupled receptors (GPCRs), there is only sparse data on how GPCRs recognize larger peptide ligands. NMR spectroscopy, molecular modeling, and double‐cycle mutagenesis studies were integrated to obtain a structural model of the peptide hormone neuropeptide Y (NPY) bound to its human G‐protein‐coupled Y₂ receptor (Y₂R). Solid‐state NMR measurements of specific isotope‐labeled NPY in complex with in vitro folded Y₂R reconstituted into phospholipid bicelles provided the bioactive structure of the peptide. Guided by solution NMR experiments, it could be shown that the ligand is tethered to the second extracellular loop by hydrophobic contacts. The C‐terminal α‐helix of NPY, which is formed in a membrane environment in the absence of the receptor, is unwound starting at T³² to provide optimal contacts in a deep binding pocket within the transmembrane bundle of the Y₂R.     

Virus-ligand interactions: identification and characterization of ligand binding by NMR spectroscopy

We demonstrate the detection and characterization of ligand binding to viruses via NMR. To illustrate the methodology, the interaction of an antiviral compound with human rhinovirus serotype 2 (HRV2) was investigated. Specific interaction of a capsid-binding inhibitor and native HRV2 was monitored utilizing saturation transfer difference (STD) NMR. STD NMR experiments at atomic resolution allowed those regions of the ligand that are involved in the interaction with the virus to be determined. The approach allows for (i) the fast and robust assessment of binding, (ii) the determination of the ligand binding epitope at atomic resolution without the necessity to crystallize virus-ligand complexes, and (iii) the reuse of the virus in subsequent assays. This methodology enables one to easily identify binding of drugs, peptides, and receptor or antibody fragments to the viral capsid.     

A Combination of Spin Diffusion Methods for the Determination of Protein–Ligand Complex Structural Ensembles

Structure‐based drug design (SBDD) is a powerful and widely used approach to optimize affinity of drug candidates. With the recently introduced INPHARMA method, the binding mode of small molecules to their protein target can be characterized even if no spectroscopic information about the protein is known. Here, we show that the combination of the spin‐diffusion‐based NMR methods INPHARMA, trNOE, and STD results in an accurate scoring function for docking modes and therefore determination of protein–ligand complex structures. Applications are shown on the model system protein kinase A and the drug targets glycogen phosphorylase and soluble epoxide hydrolase (sEH). Multiplexing of several ligands improves the reliability of the scoring function further. The new score allows in the case of sEH detecting two binding modes of the ligand in its binding site, which was corroborated by X‐ray analysis.     

An NMR Method To Pinpoint Supramolecular Ligand Binding to Basic Residues on Proteins

Targeting protein surfaces involved in protein–protein interactions by using supramolecular chemistry is a rapidly growing field. NMR spectroscopy is the method of choice to map ligand‐binding sites with single‐residue resolution by amide chemical shift perturbation and line broadening. However, large aromatic ligands affect NMR signals over a greater distance, and the binding site cannot be determined unambiguously by relying on backbone signals only. We herein employed Lys‐ and Arg‐specific H2(C)N NMR experiments to directly observe the side‐chain atoms in close contact with the ligand, for which the largest changes in the NMR signals are expected. The binding of Lys‐ and Arg‐specific supramolecular tweezers and a calixarene to two model proteins was studied. The H2(C)N spectra track the terminal CH₂ groups of all Lys and Arg residues, revealing significant differences in their binding kinetics and chemical shift perturbation, and can be used to clearly pinpoint the order of ligand binding.     

Determination of the binding specificity of the 12S subunit of the transcarboxylase by saturation transfer difference NMR

In this study we present the characterization of the interaction of biotin and methylmalonyl-CoA (MMCoA) with the carboxyltransferase subunit (12S) from the transcarboxylase (TC) from Propionibacterium shermanii. This biotin dependent multienzyme complex catalyses the transfer of carbon dioxide from methylmalonyl-CoA (MMCoA) to pyruvate. The Saturation Transfer Difference NMR (STD) technique was performed to determine the binding epitope from biotin and MMCoA to the 12S subunit. We could show by titrations during STD experiments that biotin and MMCoA bind cooperatively in one binding pocket.     

Group epitope mapping by saturation transfer difference NMR to identify segments of a ligand in direct contact with a protein receptor.

A protocol based on saturation transfer difference (STD) NMR spectra was developed to characterize the binding interactions at an atom level, termed group epitope mapping (GEM). As an example we chose the well-studied system of galactose binding to the 120-kDa lectin Ricinus communis agglutinin I (RCA(120)). As ligands we used methyl beta-D-galactoside and a biantennary decasaccharide. Analysis of the saturation transfer effects of methyl beta-D-galactoside showed that the H2, H3, and H4 protons are saturated to the highest degree, giving evidence of their close proximity to protons of the RCA(120) lectin. The direct interaction of the lectin with this region of the galactose is in excellent agreement with results obtained from the analysis of the binding specificities of many chemically modified galactose derivatives (Bhattacharyya, L.; Brewer, C. F. Eur. J. Biochem. 1988, 176, 207-212). This new NMR technique can identify the binding epitope of even complex ligands very quickly, which is a great improvement over time-consuming chemical modifications. Efficient GEM benefits from a relatively high off rate of the ligand and a large excess of the ligand over the receptor. Even for a ligand like the biantennary decasaccharide with micromolar binding affinity, the binding epitopes could easily be mapped to the terminal beta-D-Gal-(1-4)-beta-D-GlcNAc (beta-D-GlcNAc = N-acetyl-D-glucosamine) residues located at the nonreducing end of the two carbohydrate chains. The binding contribution of the terminal galactose residue is stronger than those of the penultimate GlcNAc residues. We could show that the GlcNAc residues bind "edge-on" with the region from H2 to H4, making contact with the protein. Analysis of STD NMR experiments performed under competitive conditions proved that the two saccharides studied bind at the same receptor site, thereby ruling out unspecific binding.     

Coarse‐grained molecular dynamics simulations of protein–ligand binding

Coarse‐grained molecular dynamics (CGMD) simulations with the MARTINI force field were performed to reproduce the protein–ligand binding processes. We chose two protein–ligand systems, the levansucrase–sugar (glucose or sucrose), and LinB–1,2‐dichloroethane systems, as target systems that differ in terms of the size and shape of the ligand‐binding pocket and the physicochemical properties of the pocket and the ligand. Spatial distributions of the Coarse‐grained (CG) ligand molecules revealed potential ligand‐binding sites on the protein surfaces other than the real ligand‐binding sites. The ligands bound most strongly to the real ligand‐binding sites. The binding and unbinding rate constants obtained from the CGMD simulation of the levansucrase–sucrose system were approximately 10 times greater than the experimental values; this is mainly due to faster diffusion of the CG ligand in the CG water model. We could obtain dissociation constants close to the experimental values for both systems. Analysis of the ligand fluxes demonstrated that the CG ligand molecules entered the ligand‐binding pockets through specific pathways. The ligands tended to move through grooves on the protein surface. Thus, the CGMD simulations produced reasonable results for the two different systems overall and are useful for studying the protein–ligand binding processes. © 2014 Wiley Periodicals, Inc.     

Synthesis and optimization of peptidomimetics as HIV entry inhibitors against the receptor protein CD4 using STD NMR and ligand docking

We recently described the design and synthesis of a novel CD4 binding peptidomimetic as a potential HIV entry inhibitor with a KD value of approximately 35 microM and a high proteolytic stability [A. T. Neffe and B. Meyer, Angew. Chem., Int. Ed., 2004, 43, 2937-2940]. Based on saturation transfer difference (STD) NMR analyses and docking studies of peptidomimetics we now report the rational design, synthesis, and binding properties of 11 compounds with improved binding affinity. Surface plasmon resonance (SPR) resulted in a KD = 10 microM for the best peptidomimetic XI, whose binding affinity is confirmed by STD NMR (KD = 9 microM). The STD NMR determined binding epitope of the ligand indicates a very similar binding mode as that of the lead structure. The binding studies provide structure activity relationships and demonstrate the utility of this approach.     

Affinity Screening Using Competitive Binding with Fluorine‐19 Hyperpolarized Ligands

Fluorine‐19 NMR and hyperpolarization form a powerful combination for drug screening. Under a competitive equilibrium with a selected fluorinated reporter ligand, the dissociation constant (K_D) of other ligands of interest is measurable using a single‐scan Carr–Purcell–Meiboom–Gill (CPMG) experiment, without the need for a titration. This method is demonstrated by characterizing the binding of three ligands with different affinities for the serine protease trypsin. Monte Carlo simulations show that the highest accuracy is obtained when about one‐half of the bound reporter ligand is displaced in the binding competition. Such conditions can be achieved over a wide range of affinities, allowing for rapid screening of non‐fluorinated compounds when a single fluorinated ligand for the binding pocket of interest is known.     

Free enthalpies of replacing water molecules in protein binding pockets

Water molecules in the binding pocket of a protein and their role in ligand binding have increasingly raised interest in recent years. Displacement of such water molecules by ligand atoms can be either favourable or unfavourable for ligand binding depending on the change in free enthalpy. In this study, we investigate the displacement of water molecules by an apolar probe in the binding pocket of two proteins, cyclin-dependent kinase 2 and tRNA-guanine transglycosylase, using the method of enveloping distribution sampling (EDS) to obtain free enthalpy differences. In both cases, a ligand core is placed inside the respective pocket and the remaining water molecules are converted to apolar probes, both individually and in pairs. The free enthalpy difference between a water molecule and a CH₃ group at the same location in the pocket in comparison to their presence in bulk solution calculated from EDS molecular dynamics simulations corresponds to the binding free enthalpy of CH₃ at this location. From the free enthalpy difference and the enthalpy difference, the entropic contribution of the displacement can be obtained too. The overlay of the resulting occupancy volumes of the water molecules with crystal structures of analogous ligands shows qualitative correlation between experimentally measured inhibition constants and the calculated free enthalpy differences. Thus, such an EDS analysis of the water molecules in the binding pocket may give valuable insight for potency optimization in drug design.     

Epitope mapping of ligand-receptor interactions by diffusion NMR

A novel method based on diffusion NMR for the epitope mapping of ligand binding is presented. The intermolecular NOE builds up during a long diffusion period and creates a deviation from the linearity. The ligand proton nearest the protein generates the strongest NOE from protein during the diffusion period and has the largest deviation. Therefore, this diffusion artifact can be used to characterize the ligand binding epitope. The concept was investigated using dihydrofolate reductase (DHFR) and its ligand trimethoprim (TMP), and the epitope map of TMP on DHFR generated with this method is in excellent agreement with the structural and dynamic studies by crystallography and NMR, as well as the medicinal chemistry results.     

Determination of the binding epitope of lidocaine with AGP: minimizing the effects of nonspecific binding in saturation transfer difference experiments

The ligand epitope map and the effect of nonspecific binding is assessed for lidocaine binding to α₁-acid glycoprotein using the saturation transfer difference nuclear magnetic resonance experiment performed as a function of the ligand/protein ratio. The experimental design tested two different approaches for preparing solutions with various ligand/protein ratios; holding the protein concentration constant and increasing the ligand concentration; and holding the ligand concentration constant while decreasing the protein concentration. Nonspecific binding effects were more prevalent in experiments in which the ligand concentration was increased, although spectra with higher signal-to-noise ratios were obtained under these conditions. The epitope map determined for achiral lidocaine is compared with previously determined results for the (R)- and (S)-enantiomers of propranolol. The weaker binding affinity of lidocaine may be partially attributed to steric hindrance by the lidocaine N-ethyl groups which may prevent close contact of the lidocaine amine with the negatively charged amino acids at the apex of the protein binding pocket.     

Characterization of Ligand Binding by Saturation Transfer Difference NMR Spectroscopy

Fast identification of binding activity directly from mixtures of potential ligands is possible with the NMR method described, which is based on saturation transfer to molecules in direct contact to a protein. In addition, the ligand's binding epitope is easily identified. High sensitivity and ease of use are the principal advantages of this method. The picture shows the normal 1D NMR spectrum of a mixture and the spectrum obtained by applying the STD method, which exclusively shows signals from molecules with binding affinity.     

Competition STD NMR for the detection of high-affinity ligands and NMR-based screening

The reported competition STD NMR method combines saturation transfer difference (STD) NMR with competition binding experiments to allow the detection of high-affinity ligands that undergo slow chemical exchange on the NMR time-scale. With this technique, the presence of a competing high-affinity ligand in the compound mixture can be detected by the disappearance or reduction of the STD signals of a low-affinity indicator ligand. This is demonstrated on a BACE1 (beta-site amyloid precursor protein cleaving enzyme 1) protein-inhibitor system. This method can also be used to derive an approximate value, or a lower limit, for the dissociation constant of the potential ligand based on the reduction of the signal intensity of the STD indicator, which is illustrated on an HSA (human serum albumin) model system. This leads to important applications of the competition STD NMR method for lead discovery: it can be used (i) for compound library screening against a broad range of drug targets to identify both high- and low-affinity ligands and (ii) to rank order analogs rapidly and derive structure-activity relationships, which are used to optimize these NMR hits into viable drug leads.     

Pose scoring by NMR

Recently, we have developed a fast approach to calculate NMR chemical shifts using the divide and conquer method at the semiempirical level. To demonstrate the utility of this approach for characterizing protein-ligand interactions, we used the deviation of calculated chemical shift perturbations from experiment to determine the orientation of a ligand (GPI-1046) in the binding pocket of the FK506 binding protein (FKBP12). Moreover, we were able to select the native state of the ligand from a collection of decoy poses. A key hydrogen bond between O1 and HN in Ile56 was also identified. Our results suggest that ligand-induced chemical shift perturbations can be used to refine protein/ligand structures.