Fragment-based screening using X-ray crystallography and NMR spectroscopy

Approaches which start from a study of the interaction of very simple molecules (fragments) with the protein target are proving to be valuable additions to drug design. Fragment-based screening allows the complementarity between a protein active site and drug-like molecules to be rapidly and effectively explored, using structural methods. Recent improvements in the intensities of laboratory X-ray sources permits the collection of greater amounts of high-quality diffraction data and have been matched by developments in automation, crystallisation and data analysis. Developments in NMR screening, including the use of cryogenically cooled NMR probes and (19)F-containing reporter molecules have expanded the scope of this technique, while increasing the availability of binding site and quantitative affinity data for the fragments. Application of these methods has led to a greater knowledge of the chemical variety, structural features and energetics of protein-fragment interactions. While fragment-based screening has already been shown to reduce the timescales of the drug discovery process, a more detailed characterisation of fragment screening hits can reveal unexpected similarities between fragment chemotypes and protein active sites leading to improved understanding of the pharmacophores and the re-use of this information against other protein targets.     

Structural analysis of lipocalin-type prostaglandin D synthase complexed with biliverdin by small-angle X-ray scattering and multi-dimensional NMR

Lipocalin-type prostaglandin D synthase (L-PGDS) acts as both a PGD₂ synthase and an extracellular transporter for small lipophilic molecules. From a series of biochemical studies, it has been found that L-PGDS has an ability to bind a variety of lipophilic ligands such as biliverdin, bilirubin and retinoids in vitro. Therefore, we considered that it is necessary to clarify the molecular structure of L-PGDS upon binding ligand in order to understand the physiological relevance of L-PGDS as a transporter protein. We investigated a molecular structure of L-PGDS/biliverdin complex by small-angle X-ray scattering (SAXS) and multi-dimensional NMR measurements, and characterized the binding mechanism in detail. SAXS measurements revealed that L-PGDS has a globular shape and becomes compact by 1.3 Å in radius of gyration on binding biliverdin. NMR experiments revealed that L-PGDS possessed an eight-stranded antiparallel β-barrel forming a central cavity. Upon the titration with biliverdin, some cross-peaks for residues surrounding the cavity and EF-loop and H2-helix above the β-barrel shifted, and the intensity of other cross-peaks decreased with signal broadenings in ¹H–¹⁵N heteronuclear single quantum coherence spectra. These results demonstrate that L-PGDS holds biliverdin within the β-barrel, and the conformation of the loop regions above the β-barrel changes upon binding biliverdin. Through such a conformational change, the whole molecule of L-PGDS becomes compact.     

Fibrinogen-catecholamine interaction as observed by NMR and Fourier transform infrared spectroscopy

In this work, the interactions between the main catecholamines-epinephrine and norepinephrine-and fibrinogen were investigated by NMR and Fourier transform infrared spectroscopies. The two hormones were found to interact with fibrinogen and to affect the protein secondary structure to a different extent. In particular, the protein selectively binds epinephrine at both the basal and stress concentrations, while it shows a weak nonspecific interaction with norepinephrine. The interaction with the stress level of epinephrine leads to drastic protein conformational changes, whereas norepinephrine does not affect fibrinogen secondary structure, even at stress concentration.     

Proteolytic E-cadherin activation followed by solution NMR and X-ray crystallography

Cellular adhesion by classical cadherins depends critically on the exact proteolytic removal of their N-terminal prosequences. In this combined solution NMR and X-ray crystallographic study, the consequences of propeptide cleavage of an epithelial cadherin construct (domains 1 and 2) were followed at atomic level. At low protein concentration, the N-terminal processing induces docking of the tryptophan-2 side-chain into a binding pocket on the same molecule. At high concentration, cleavage induces dimerization (KD=0.72 mM, k(off)=0.7 s(-1)) and concomitant intermolecular exchange of the betaA-strands and the tryptophan-2 side-chains. Thus, the cleavage represents the switch from a nonadhesive to the functional form of cadherin.     

Mirroring perfection: the structure of methylglyoxal synthase complexed with the competitive inhibitor 2-phosphoglycolate

The crystal structure of the transition-state analogue 2-phosphoglycolate (2PG) bound to methylglyoxal synthase (MGS) is presented at a resolution of 2.0 A. This structure is very similar to the previously determined structure of MGS complexed to formate and phosphate. Since 2PG is a competitive inhibitor of both MGS and triosephosphate isomerase (TIM), the carboxylate groups of each bound 2PG from this structure and the structure of 2PG bound to TIM were used to align and compare the active sites despite differences in their protein folds. The distances between the functional groups of Asp 71, His 98, His 19, and the carboxylate oxygens of the 2PG molecule in MGS are similar to the corresponding distances between the functional groups of Glu 165, His 95, Lys 13, and the carboxylate oxygens of the 2PG molecule in TIM. However, these spatial relationships are enantiomorphic to each other. Consistent with the known stereochemical data, the catalytic base Asp 71 is positioned on the opposite face of the 2PG-carboxylate plane as Glu 165 of TIM. Both His 98 of MGS and His 95 of TIM are in the plane of the carboxylate of 2PG, suggesting that these two residues are homologous in function. While His 19 of MGS and Lys 13 of TIM appear on the opposite face of the 2PG carboxylate plane, their relative location to the 2PG molecule is quite different, suggesting that they probably have different functions. Most remarkably, unlike the coplanar structure found in the 2PG molecule bound to TIM, the torsion angle around the C1-C2 bond of 2PG bound to MGS brings the phosphoryl moiety out of the molecule's carboxylate plane, facilitating elimination. Further, the superimposition of this structure with the structure of MGS bound to formate and phosphate suggests a model for the enzyme bound to the first transition state.     

Relative stability of protein structures determined by X-ray crystallography or NMR spectroscopy: a molecular dynamics simulation study

The relative stability of protein structures determined by either X-ray crystallography or nuclear magnetic resonance (NMR) spectroscopy has been investigated by using molecular dynamics simulation techniques. Published structures of 34 proteins containing between 50 and 100 residues have been evaluated. The proteins selected represent a mixture of secondary structure types including all alpha, all beta, and alpha/beta. The proteins selected do not contain cysteine-cysteine bridges. In addition, any crystallographic waters, metal ions, cofactors, or bound ligands were removed before the systems were simulated. The stability of the structures was evaluated by simulating, under identical conditions, each of the proteins for at least 5 ns in explicit solvent. It is found that not only do NMR-derived structures have, on average, higher internal strain than structures determined by X-ray crystallography but that a significant proportion of the structures are unstable and rapidly diverge in simulations.     

Unusual conformational changes in 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase as revealed by X-ray crystallography and NMR

The crystal structure of Escherichia coli 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) in complex with MgADP has been determined at 1.5-A resolution with a crystallographic R factor of 0.191. The solution structure of HPPK in complex with Mg(2+) and beta,gamma-methyleneadenosine 5'-triphosphate (MgAMPPCP) has been determined using a simulated annealing protocol with 3,523 experimental NMR restraints. The root mean square deviation of the ensemble of 20 refined conformers that represent the solution structure from the mean coordinate set derived from them is 0.74 +/- 0.26 A for all backbone atoms and 0.49 +/- 0.22 A when residues Pro(14), Pro(44)-Gln(50), and Arg(84)-Pro(91) are excluded. Binding of MgADP causes significant changes in the conformation and dynamical property of three loops of HPPK that are involved in catalysis. A dramatic, unusual conformational change is that loop 3 moves away from the active center significantly with some residues moving by >17 A. The binding of MgADP also stabilizes loop 1 and loop 3 but makes loop 2 more mobile. Very similar conformational and dynamical changes are observed in the NMR solution structure of HPPK.MgAMPPCP. The conformational and dynamical changes may play important roles in both substrate binding and product release in the catalytic cycle.     

Determination of absolute configurations by X-ray crystallography and 1H NMR anisotropy

To determine the absolute configurations of chiral compounds, many spectroscopic and diffraction methods have been developed. Among them, X-ray crystallographic Bijvoet method, CD exciton chirality method, and the combination of vibrational circular dichroism and quantum mechanical calculations are of nonempirical nature. On the other hand, X-ray crystallography using a chiral internal reference, and 1H NMR spectroscopy using chiral anisotropy reagents are relative and/or empirical methods. In addition to absolute configurational determinations, preparations of enantiopure compounds are strongly desired. As chiral reagents useful for both the preparation of enantiopure compounds by HPLC separation and the simultaneous determination of their absolute configurations, we have developed camphorsultam dichlorophthalic acid (CSDP acid) for X-ray crystallography and 2-methoxy-2-(1-naphthyl)propionic acid (MalphaNP acid) for 1H NMR spectroscopy. In this review, the principles and applications of these X-ray and NMR methods are explained using mostly our own data.     

The heme environment of mouse neuroglobin: histidine imidazole plane orientations obtained from solution NMR and EPR spectroscopy as compared with X-ray crystallography

The ¹H NMR chemical shifts of the heme methyl groups of the ferriheme complex of metneuroglobin (Du et al. in J. Am. Chem. Soc. 125:8080–8081, 2003) predict orientations of the axial histidine ligands (Shokhirev and Walker in J. Biol. Inorg. Chem. 3:581–594, 1998) that are not consistent with the X-ray data (Vallone et al. in Proteins Struct. Funct. Bioinf. 56:85–94, 2004), and the EPR spectrum (Vinck et al. in J. Am. Chem. Soc. 126:4516–4517, 2004) is only marginally consistent with these data. The reasons for these inconsistencies appear to be rooted in the high degree of aqueous solution exposure of the heme group and the fact that there are no strong hydrogen-bond acceptors for the histidine imidazole N–H protons provided by the protein. Similar inconsistencies may exist for other water-soluble heme proteins, and ¹H NMR spectroscopy provides a simple means to verify whether the solution structure of the heme center is the same as or different from that in the crystalline state.     

Combining on-line spectroscopy with synchrotron and X-ray free electron laser crystallography

With the recent advances in serial crystallography methods at both synchrotron and X-ray free electron laser sources, more details of intermediate or transient states of the catalytic reactions are being revealed structurally. These structural studies of reaction dynamics drive the need for on-line in crystallo spectroscopy methods to complement the crystallography experiment. The recent applications of combined spectroscopy and crystallography methods enable on-line determination of in crystallo reaction kinetics and structures of catalytic intermediates, sample integrity, and radiation-induced sample modifications, if any, as well as heterogeneity of crystals from different preparations or sample batches. This review describes different modes of spectroscopy that are combined with the crystallography experiment at both synchrotron and X-ray free-electron laser facilities, and the complementary information that each method can provide to facilitate the structural study of enzyme catalysis and protein dynamics.     

Conformational states of ADP ribosylation factor 1 complexed with different guanosine triphosphates as studied by 31P NMR spectroscopy

Guanine nucleotide binding proteins (GNB-proteins) play an essential role in cellular signaling, acting as molecular switches, cycling between the inactive, GDP-bound form and the active, GTP-bound form. It has been shown that conformational equilibria also exist within the active form of GNB-proteins between conformational states with different functional properties. Here we present (31)P NMR data on ADP ribosylation factor 1 (Arf1), a GNB-protein involved in Golgi traffic, promoting the coating of secretory vesicles. To investigate conformational equilibria in active Arf1, the wild type and switch I mutants complexed with GTP and a variety of commonly used GTP analogues, namely, GppCH(2)p, GppNHp, and GTPγS, were analyzed. To gain deeper insight into the conformational state of active Arf1, we titrated with Cu(2+)-cyclen and GdmCl and formed the complex with the Sec7 domain of nucleotide exchange factor ARNO and an effector GAT domain. In contrast to the related proteins Ras, Ral, Cdc42, and Ran, from (31)P NMR spectroscopic view, Arf1 exists predominantly in a single conformation independent of the GTP analogue used. This state seems to correspond to the so-called state 2(T) conformation, according to Ras nomenclature, which is interacting with the effector domain. The exchange of the highly conserved threonine in position 48 with alanine led to a shift of the equilibrium toward a conformational state with typical properties obtained for state 1(T) in Ras, such as interaction with guanine nucleotide exchange factors, a lower affinity for nucleoside triphosphates, and greater sensitivity to chaotropic agents. In active Arf1(wt), the effector interacting conformation is strongly favored. These intrinsic conformational equilibria of active GNB-proteins could be a fine-tuning mechanism of regulation and thereby an interesting target for the modulation of protein activity.     

Proline cis-trans isomers in calbindin D9k observed by X-ray crystallography.

In a structure of recombinant bovine calbindin D9k, determined crystallographically to 1.6 A resolution, a proline in mixed, approximately equally populated, cis and trans conformation is observed. Isomers of this kind have not been reported in structure determinations of calbindin D9k to 2.3 A resolution or in any other crystallographically determined protein structure. The cis-trans isomerization occurs at the peptide bond between Gly42 and Pro43, which is in agreement with results from two-dimensional 1H nuclear magnetic resonance spectroscopy experiments on solutions of calbindin D9k. Alternative backbone stretches have been modeled and refined by stereochemical restrained least-squares refinement for the segment Lys41 to Pro43. The final R-value was 0.188. The structural perturbations accompanying the cis-trans isomerization are found to be very localized. The largest positional differences are observed at residue Gly42, in which the alternative positions of the oxygen atom are 3.6 A apart.     

The use of solid-state NMR and X-ray crystallography as complementary tools for studying molecular recognition

Solid-state NMR is rapidly becoming available as a routine technique for studying the structure of crystalline or noncrystalline solids. This technique has an advantage over crystallography in that single crystals are not necessary, but it has the disadvantage that the information obtained does not produce a direct picture of the molecule and its environment. On the other hand, solid-state NMR can be done on mixtures, and it gives information about phase distribution in a manner similar to that of X-ray powder pattern analysis.Crystallographic effects such as polymorphism, multiple molecules per asymmetric unit, disorder and salvation can frequently be detected using NMR. Sometimes molecular point group symmetry can also be deduced based on the number of independent nuclei that are detected. The NMR method is sensitive to changes in the electronic structure of a molecule as sensed by the nuclei, and the effects are measured as changes in the isotropic chemical shift of individual nuclei.In this paper, we will give examples of the combined use of X-ray crystallography and ¹³CP/MAS (cross polarization/magic angle spinning) NMR for studying hostguest materials and cocrystals. We have learned how to use NMR to tell us about keto/enol composition in the solid state, to detect the presence of trapped solvent molecules, to detect hydrogen-bond formation and to evaluate molecular conformation and unusual packing pattern effects. We will also present a brief background of the ¹³CP/MAS NMR technique and three case studies in which solid-state NMR and X-ray crystallography are used together to understand materials' structures and properties     

X-ray crystallography and solution NMR spectroscopy characterization of heptakis(2,3-di-O-acetyl-6-bromo-6-deoxy)cyclomaltoheptaose

Heptakis(2,3-di-O-acetyl-6-bromo-6-deoxy)cyclomaltoheptaose has been characterized in aqueous solution by 1D and 2D NMR spectroscopy and in the solid state by X-ray crystallography. In methanol solution, the acetyl groups were found to interact with both inward and outward-pointing protons. This and the strong deshielding of the bridging carbons, relative to the nonacetylated precursor, indicate macrocyclic flexibility. In the crystalline state the macrocycle exists as a methanol complex. It exhibits elliptical distortion, all glucose residues been tilted with their primary side toward the cavity. The existing strain due to the congestion of 14 acetyl groups at the secondary site is relieved by two glucose rings acquiring the rarely observed skew-boat conformation, (0)S(2), by the increased tilting of two glucose residues, as well as by minor variations of the torsion angles of the acetyl groups. The seven bromine atoms are quite accessible to nucleophiles.     

Ligand orientation of human neuroglobin obtained from solution NMR and molecular dynamics simulation as compared with X-ray crystallography

Neuroglobin, a new member of hemoprotein family, can reversibly bind oxygen and take part in many biological processes such as enzymatic reaction, signal transduction and the mitochondria function. Different from myoglobin and hemoglobin, it has a hexacoordinated heme environment, with histidyl imidazole of proximal His⁹⁶(F8) and distal His⁶⁴(E7) directly bound to the metal ion. In the present work, solution ¹H NMR spectroscopy was employed to investigate the electronic structure of heme center of wild-type met-human neuroglobin. The resonances of heme protons and key residues in the heme pocket were assigned. Two heme orientations resulting from a 180° rotation about the α-γ-meso axis with a population ratio about 2:1 were observed. Then the ¹H NMR chemical shifts of the ferriheme methyl groups were used to predict orientations of the axial ligand. The obtained axial ligand plane angle φ is consistent with that from the molecular dynamics simulation but not with those from the crystal data. Compared with mouse neuroglobin, the obtained average ligand orientation of human neuroglobin reflects the changeability of heme environment for the Ngb family.     

Synthesis of new palladium(II) complexes containing tetradentate-nitrogen donor ligands: Combined structural studies by NMR spectroscopy and X-ray crystallography

Treatment of the tetradentate (NN′N′N) N-alkylaminopyrazole ligands 3,6-dimethyl-1,8-(3,5-dimethyl-1-pyrazolyl)-3,6-diazaoctane (ddad) and 1,4-bis[2-(3,5-dimethyl-1-pyrazolyl)ethyl]piperazine (bedp) with [PdCl₂(CH₃CN)₂] in a 1:1 M/L ratio in acetonitrile produces [Pd₂Cl₄(L)] and [PdCl₂(L)] (L = ddad and bedp). Treatment of the corresponding complex [PdCl₂(L)] (L = ddad, bedp) in the presence of AgBF₄ in CH₂Cl₂/methanol (2:1) or NaBF₄ in acetonitrile gives [Pd(L)](BF₄)₂. The Pd(II) complexes have been characterised by elemental analyses, conductivity measurements, IR and ¹H and ¹³C{¹H} NMR spectroscopies when possible. The X-ray structure of the complex [Pd(ddad)]Cl₂ · 3H₂O has been determined. The Pd(II) is coordinated to the ddad ligand by two nitrogen atoms of pyrazolyl groups and two nitrogen atoms of the amine groups, in a slightly distorted square-planar geometry.     

X-ray crystallography and the amino acid sequence of lysozyme

Certain parts of the amino acid sequence of hen's-egg lysozyme were reinvestigated. The parts in question are those in which X-ray-crystallographic data were used to decide between the two reported versions of the sequence that were based on chemical data. We conclude that the X-ray evidence, in the main, led to the correct choice between the sequences.     

Ring current effects in Adriamycin-flavin mononucleotide complexation as observed by 1H FT NMR spectroscopy

The addition of Adriamycin to a solution containing flavin mononucleotide (FMN) resulted in an upfield shift in the signals of the aromatic ring protons H(6,9) and the 8α, 7α methyl protons of FMN. The chemical shift of the H(6,9) and of the 8α and 7α methyl proton signals of FMN decreased from 7.92, 2.56 and 2.46 ppm, respectively, in the absence of Adriamycin to 7.61, 2.42 and 2.36 ppm, respectively, at 3 mM Adriamycin. Concomitant increases in the linewidth of aromatic and methyl proton siqnals of FMN were also observed. Variable temperature studies over the range of 5 to 43° showed an increase in the chemical shift of both the aromatic and aliphatic proton signals with increasing temperatures. These results suggest that FMN and Adriamycin form a complex via ring-ring stacking.