High-throughput NMR-based screening with competition binding experiments

The Achilles heel of ligand-based NMR screening methods is their failure to detect high-affinity ligands and molecules that bind covalently to the receptor. We have developed a novel approach for performing high-throughput screening with NMR spectroscopy that overcomes this limitation. The method also permits detection of potential high-affinity molecules that are only marginally soluble, thus significantly enlarging the diversity of compounds amenable to NMR screening. The techniques developed utilize transverse and/or selective longitudinal relaxation parameters in combination with competition binding experiments. Mathematical expressions are derived for proper setup of the NMR experiments and for extracting an approximate value of the binding constant for the identified ligand from a single-point measurement. With this approach it is possible to screen thousands of compounds in a short period of time against protein or DNA and RNA fragments. The methodology can also be applied for screening plant and fungi extracts.     

Application of NMR screening techniques for observing ligand binding with a protein receptor

Water ligand observed via gradient spectroscopy (WaterLOGSY), saturation transfer difference and NOE pumping NMR techniques were used to identify ligand binding with a receptor. Although these experiments were originally designed to observe ligands in complexes, their application is limited by the affinity of ligands towards target molecules. Here the improved WaterLOGSY pulse sequence was developed by incorporating the double pulsed field gradient spin-echo and gradient-tailored excitation WATERGATE sequences. The efficiency of these ligand-observed NMR screening techniques was investigated using the ribonuclease T1-inhibitor system.     

Nuclear magnetic resonance relaxometry as a spectroscopic probe of the coordination sphere of a paramagnetic metal bound to a humic acid mixture

Protons on water molecules are strongly affected by paramagnetic ions. Since the acid-base properties of water facilitate rapid proton exchange, a single proton nuclear magnetic resonance (NMR) signal is seen in aqueous solutions of paramagnetic ions. Proton relaxation times are significantly affected by paramagnetic species and the readily detectable single signal serves as a powerful amplifier of the information contained concerning the protons in the paramagnetic environment. Where water molecules coordinated to free paramagnetic ions and to metal complexes of ligands that form non-labile (on the NMR time scale) complexes, the effects on water in the two environments can be distinguished. This can provide information on the nature of the ligand binding sites. The example of Cu²⁺ bound to the Laurentian humic acid mixture reported here using convenient low field NMR relaxometers shows that the information can enrich our understanding of complexation and speciation in the presence of complex mixture ligands characteristic of natural water systems. In this case, the data underline the role of aggregation and conformation in defining the complexation sites.     

Impact of planarity of unfused aromatic molecules on G-quadruplex binding: learning from isaindigotone derivatives

G-quadruplex structures are a new class of attractive targets for DNA-interactive anticancer agents. The primary building block of this structure is the G-quartet, which is composed of four coplanar guanines and serves as the major binding site for small molecules. NMR studies and molecular dynamics simulations have suggested that the planarity of G-quartet surface has been highly dynamic in solution. To better investigate how the planarity of unfused aromatic ligand impacts on its quadruplex binding properties, a variety of planarity controllable isaindigotone derivatives were designed and synthesized. The interaction of G-quadruplex DNA with these designed ligands was systematically explored using a series of biophysical studies. The FRET-melting, SPR, and CD spectroscopy results showed that reducing the planarity of their unfused aromatic core resulted in their decreased binding affinity and stabilization ability for G-quadruplex. NMR studies also suggested that these compounds could stack on the G-quartet surface. Such results are in parallel with subsequent molecular modeling studies. A detailed binding energy analysis indicated that van der Waals energy (ΔE(vdw)) and entropy (TΔS) are responsible for their decreased quadruplex binding and stabilization effect. All these results provided insight information about how quadruplex recognition could be controlled by adjusting the planarity of ligands, which shed light on further development of unfused aromatic molecules as optimal G-quadruplex binding ligands.     

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.     

19F NMR isotropic chemical shift for efficient screening of fluorinated fragments which are racemates and/or display multiple conformers

Fluorine ligand‐based NMR spectroscopy is now an established method for performing binding screening against a macromolecular target. Typically, the transverse relaxation rate of the fluorine signals is monitored in the absence and presence of the target. However, useful structural information can sometimes be obtained from the analysis of the fluorine isotropic chemical shift. This is particularly relevant for molecules that are racemates and/or display multiple conformers. The large difference in fluorine isotropic chemical shift between free and bound state deriving mainly from the breaking and/or making of intramolecular and/or intermolecular hydrogen bonds allows the detection of very weak affinity ligands. According to our experimental results, racemates should always be included in the generation of the fluorinated fragment libraries. The selection or the availability of only one of the enantiomers for the fluorinated screening library could result in missing relevant chemical scaffold motifs.     

Binding interaction of [Re(H2O)3(CO)3]+ with the DNA fragment d(CpGpG)

Insights into the interaction of the [Re(H2O)3(CO)3]+ complex (1) with the DNA fragment d(CpGpG) have been obtained by one- (1D) and two-dimensional (2D) NMR spectroscopy. The H8 resonances of the single major [Re(H2O)d(CpGpG)(CO)3]- adduct (2) exhibit pH-independent chemical shift changes attributable to metal N7 binding. The structure of this adduct has been characterized by molecular modeling studies based on 1D and 2D NMR data. In solution, 2 shows the presence of two N7-coordinated guanine moieties in a head-to-head (HH) orientation as evidenced by G2H8/G3H8 cross-peaks in the 1H-1H NOESY NMR spectrum. The presence of the 5'-bridging phosphodiester appears to stabilize the HH1 L conformer, as was previously described for related Pt and Rh complexes.     

Design of new bidentate ligands constructed of two Hoechst 33258 units for discrimination of the length of two A3T3 binding motifs

[structure: see text] The aim of this study is to develop bidentate minor-groove binders that bind the double binding motifs cooperatively. The new bidentate ligands (1) have been designed by connecting two Hoechst 33258 units with a polyether linker for cooperative binding with two remote A3T3 sites of DNA. The linker is introduced to the benzimidazole ring so that it is located at the convex side of the Hoechst unit. DNA binding affinity of the ligands was evaluated by measuring surface plasmon resonance (SPR), circular dichroism, and fluorescence spectra. Interestingly, the bidentate ligands (1) did not show affinity to DNA1 with a single A3T3 motif but showed selective affinity to DNA2 with two A3T3 motifs. The Long Bis-H (1L) having a long polyether linker showed specific binding to DNA2(6) with two A3T3 motifs separated by six nonbinding base pairs. The Long Bis-H (1L) has also shown specific binding to the three-way junction DNA4 with two A3T3 motifs. This study has demonstrated that DNA with double binding motifs can be selectively recognized by the newly designed bidentate ligands.     

In silico footprinting of ligands binding to the minor groove of DNA

The sequence selectivity of small molecules binding to the minor groove of DNA can be predicted by "in silico footprinting". Any potential ligand can be docked in the minor groove and then moved along it using simple simulation techniques. By applying a simple scoring function to the trajectory after energy minimization, the preferred binding site can be identified. We show application to all known noncovalent binding modes, namely 1:1 ligand:DNA binding (including hairpin ligands) and 2:1 side-by-side binding, with various DNA base pair sequences and show excellent agreement with experimental results from X-ray crystallography, NMR, and gel-based footprinting.     

meso-[{Ru(phen)2}2(mu-HAT)]4+: a high-affinity DNA hairpin probe {HAT = 1,4,5,8,9,12-hexaazatriphenylene; phen = 1,10-phenanthroline}

1H NMR spectroscopy and fluorescent intercalator displacement (FID) assays have been used to investigate the DNA-binding abilities of two series of dinuclear polypyridyl ruthenium(II) complexes of the form [{Ru(L)2}2(mu-BL)]4+ {L = 2,2'-bipyridine (bpy), 4,4'-dimethyl-2,2'-bipyridine (Me2bpy), 1,10-phenanthroline (phen), or 4,7-dimethyl-1,10-phenanthroline (Me2phen); BL = 2,2'-bipyrimidine (bpm) or 1,4,5,8,9,12-hexaazatriphenylene (HAT)}. Preliminary FID surveys of these metal complexes against a variety of different oligonucleotides revealed that those complexes based upon the HAT bridging ligand induced greater fluorescence decreases in dye-bound DNA than did their bpm-bridged counterparts, suggesting a higher binding affinity by the HAT-bridged species. Furthermore, the greatest fluorescence decreases were typically observed in an oligonucleotide featuring a six-base hairpin loop. The apparent binding affinity of the metal complexes was also found to be a function of the stereochemistry and identity of the terminal ligands of the complex. The meso (DeltaLambda) stereoisomer generally induced greater fluorescence decreases than did either enantiomer (DeltaDelta or LambdaLambda), phen-based terminal ligands performed better than bpy-based terminal ligands, and those terminal ligands with methyl substituents demonstrated stronger apparent binding than did their non-methylated analogues. NMR experiments on meso-[{Ru(phen)2}2(mu-HAT)]4+ and meso-[{Ru(Me2phen)2}2(mu-HAT)]4+ demonstrated that both complexes bound with high affinity to the six-base hairpin oligonucleotide at the stem-loop interface and provided evidence to support stronger binding by the methylated species. meso-[{Ru(phen)2}2(mu-HAT)]4+ was found to bind poorly to duplex DNA and smaller four-base hairpin loops in FID and NMR experiments, whereas FID data suggest that the methylated analogue binds relatively strongly to most oligonucleotide sequences (the four- and six-base hairpins in particular). These results demonstrate that binding affinity can come at the expense of selectivity, with meso-[{Ru(phen)2}2(mu-HAT)]4+ proving to be an efficient compromise between the two as a high-affinity DNA hairpin probe.     

Heteronuclear NMR spectroscopy for lysine NH(3) groups in proteins: unique effect of water exchange on (15)N transverse relaxation

In this paper, we present a series of heteronuclear NMR experiments for the direct observation and characterization of lysine NH3 groups in proteins. In the context of the HoxD9 homeodomain bound specifically to DNA we were able to directly observe three cross-peaks, arising from lysine NH3 groups, with 15N chemical shifts around approximately 33 ppm at pH 5.8 and 35 degrees C. Measurement of water-exchange rates and various types of 15N transverse relaxation rates for these NH3 groups, reveals that rapid water exchange dominates the 15N relaxation for antiphase coherence with respect to 1H through scalar relaxation of the second kind. As a consequence of this phenomenon, 15N line shapes of NH3 signals in a conventional 1H-15N heteronuclear single quantum coherence (HSQC) correlation experiment are much broader than those of backbone amide groups. A 2D 1H-15N correlation experiment that exclusively observes in-phase 15N transverse coherence (termed HISQC for heteronuclear in-phase single quantum coherence spectroscopy) is independent of scalar relaxation in the t(1) (15N) time domain and as a result exhibits strikingly sharper 15N line shapes and higher intensities for NH3 cross-peaks than either HSQC or heteronuclear multiple quantum coherence (HMQC) correlation experiments. Coherence transfer through the relatively small J-coupling between 15Nzeta and 13Cepsilon (4.7-5.0 Hz) can be achieved with high efficiency by maintaining in-phase 15N coherence owing to its slow relaxation. With the use of a suite of triple resonance experiments based on the same design principles as the HISQC, all the NH3 cross-peaks observed in the HISQC spectrum could be assigned to lysines that directly interact with DNA phosphate groups. Selective observation of functional NH3 groups is feasible because of hydrogen bonding or salt bridges that protect them from rapid water exchange. Finally, we consider the potential use of lysine NH3 groups as an alternative probe for larger systems as illustrated by data obtained on the 128-kDa enzyme I dimer.     

Differential Epitope Mapping by STD NMR Spectroscopy To Reveal the Nature of Protein–Ligand Contacts

Saturation transfer difference (STD) NMR spectroscopy is extensively used to obtain epitope maps of ligands binding to protein receptors, thereby revealing structural details of the interaction, which is key to direct lead optimization efforts in drug discovery. However, it does not give information about the nature of the amino acids surrounding the ligand in the binding pocket. Herein, we report the development of the novel method differential epitope mapping by STD NMR (DEEP‐STD NMR) for identifying the type of protein residues contacting the ligand. The method produces differential epitope maps through 1) differential frequency STD NMR and/or 2) differential solvent (D₂O/H₂O) STD NMR experiments. The two approaches provide different complementary information on the binding pocket. We demonstrate that DEEP‐STD NMR can be used to readily obtain pharmacophore information on the protein. Furthermore, if the 3D structure of the protein is known, this information also helps in orienting the ligand in the binding pocket.     

Platinum phenanthroimidazole complexes as G-quadruplex DNA selective binders

Complexes that bind and stabilize G-quadruplex DNA structures are of significant interest due to their potential to inhibit telomerase and halt tumor cell proliferation. We here report the synthesis of the first Pt(II) G-quadruplex selective molecules, containing pi-extended phenanthroimidazole ligands. Binding studies of these complexes with duplex and quadruplex d(T(4)G(4)T(4))(4) DNA were performed. Intercalation to duplex DNA was established through UV/Vis titration, CD spectroscopy, and thermal denaturation studies. Significantly stronger binding affinity of these phenanthroimidazole Pt(II) complexes to G-quadruplex DNA was observed by UV/Vis spectroscopy and competitive equilibrium dialysis studies. Observed binding constants to quadruplex DNA were nearly two orders of magnitude greater than for duplex DNA. Circular dichroism studies show that an increase in pi-surface leads to a significant increase in the thermal stability of the Pt(II)/quadruplex DNA complex. The match in the pi-surface of these phenanthroimidazole Pt(II) complexes with quadruplex DNA was further substantiated by molecular modeling studies. Numerous favorable pi-stacking interactions with the large aromatic surface of the intermolecular G-quadruplex, and unforeseen hydrogen bonds between the ancillary ethylenediamine ligands and the quadruplex phosphate backbone are predicted. Thus, both biological and computational studies suggest that coupling the square-planar geometry of Pt(II) with pi-extended ligands results in a simple and modular method to create effective G-quadruplex selective binders, which can be readily optimized for use in telomerase-based antitumor therapy.     

Polymer-tethered ligand-receptor interactions between surfaces II

In this paper, we analyze the thermodynamics of interactions between surfaces mediated by polymer-tethered ligand-receptor binding. From statistical thermodynamics calculations, we obtain an effective two-dimensional binding constant reflecting contributions from the microscopic binding affinity as well as from the conformation entropy of the polymer tethers. The total interaction is a result of the interplay between attractive binding and repulsion due to confinement of the polymer chains. We illustrate the differences between three scenarios when the binding molecules are (1) immobile, (2) mobile with a fixed density, and (3) mobile with a fixed chemical potential (connected to a reservoir), which correspond to different biological or experimental situations. The key features of interactions, including the range of adhesion (onset of binding) and the equilibrium separation, can be obtained from scaling analysis and are verified by numerical solutions. In addition, we also extend our method of treating the quenched case with immobile ligands and receptors developed in a previous paper [Martin, J. I.; Zhang, C.-Z.; Wang, Z.-G. J. Polym. Sci., Part B: Polym. Phys. 2006, 44, 2621-2637] as a density expansion in terms of both ligand and receptor densities. Finally, we examine several cases having ligand-receptor pairs with different tether lengths and binding affinities, and/or nonbinding linear polymers as steric repellers. Such systems can exhibit complex interactions such as a double-well potential, or a bound state with an adjustable barrier (due to the repellers), which have both biological and bioengineering relevance.     

Use of 19F NMR spectroscopy to screen chemical libraries for ligands that bind to proteins

Identification of compounds from chemical libraries that bind to macromolecules by use of NMR spectroscopy has gained increasing importance during recent years. A simple methodology based on (19)F NMR spectroscopy for the screening of ligands that bind to proteins, which also provides qualitative information about relative binding strengths and the presence of multiple binding sites, is presented here. A library of fluorinated compounds was assembled and investigated for binding to the two bacterial chaperones PapD and FimC, and also to human serum albumin (HSA). It was found that library members which are bound to a target protein could be identified directly from line broadening and/or induced chemical shifts in a single, one-dimensional (19)F NMR spectrum. The results obtained for binding to PapD using (19)F NMR spectroscopy agreed well with independent studies based on surface plasmon resonance, providing support for the versatility and accuracy of the technique. When the library was titrated to a solution of PapD chemical shift and linewidth changes were observed with increasing ligand concentration, which indicated the presence of several binding sites on PapD and enabled the assessment of relative binding strengths for the different ligands. Screening by (19)F NMR spectroscopy should thus be a valuable addition to existing NMR techniques for evaluation of chemical libraries in bioorganic and medicinal chemistry.     

Thermodynamics of binding of 2-methoxy-3-isopropylpyrazine and 2-methoxy-3-isobutylpyrazine to the major urinary protein

In the present study we examine the thermodynamics of binding of two related pyrazine-derived ligands to the major urinary protein, MUP-I, using a combination of isothermal titration calorimetry (ITC), X-ray crystallography, and NMR backbone (15)N and methyl side-chain (2)H relaxation measurements. Global thermodynamics data derived from ITC indicate that binding is driven by favorable enthalpic contributions, rather than the classical entropy-driven hydrophobic effect. Unfavorable entropic contributions from the protein backbone and side-chain residues in the vicinity of the binding pocket are partially offset by favorable entropic contributions at adjacent positions, suggesting a "conformational relay" mechanism whereby increased rigidity of residues on ligand binding are accompanied by increased conformational freedom of side chains in adjacent positions. The principal driving force governing ligand affinity and specificity can be attributed to solvent-driven enthalpic effects from desolvation 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.     

Analysis of hairpin polyamide complexes having DNA binding sites in close proximity

The binding of two hairpin polyamide ligands at adjacent sites on DNA has been studied using NMR spectroscopy. The ligands ImPyPy-gamma-PyPyPy-Gly-Dp and Ac-ImPyPy-gamma-PyPyPy-Gly-Dp were studied binding to oligomers containing one or two matched binding sites: 5'-XGTTA-3' and 5'-TAACXNGTTA-3', where X is G, C, or A and N = 0, 1 or 2. At these sites the C-terminal ring shows an equilibrium between normal and inverted conformations. Better binding was observed with the ligand running 5' to 3' along the contacted strand than in the opposite direction. Complexes of DNAs with two binding sites indicated that at least one spacing base pair was required, and that the identity of this base pair was not critical. Binding with 5' to 3' contact is again preferred. Demonstrated binding at adjacent sites indicates that it may be possible to engineer cooperative binding for enhanced specificity or affinity.     

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.