Use of 2D NMR, protein engineering, and molecular modeling to study the hapten-binding site of an antibody Fv fragment against 2-phenyloxazolone

Two-dimensional (2D) 1H NMR spectroscopy was used to study the hapten-binding site of a recombinant antibody Fv fragment expressed in Escherichia coli. Point mutations of residues in the CDR loops of the Fv fragment were designed in order to investigate their influence on hapten binding and to make site-specific assignments of aromatic NMR proton signals. Two tyrosines giving NOEs to the ligand 2-phenyloxazolone were identified, residue 33 in CDR1 of the heavy chain and residue 32 in CDR1 of the light chain. The benzyl portion of 2-phenyloxazolone is located between these two residues. The binding site is close to the surface of the Fv fragment. Comparison with a different anti-2-phenyloxazolone antibody, the crystal structure of which has recently been solved, shows that the general location of the hapten-binding site in both antibodies is similar. However, in the crystallographically solved antibody, the hapten is bound farther from the surface in a pocket created by a short CDR3 loop of the heavy chain. In the binding site identified in the Fv fragment studied in this report, this space is probably filled by the extra seven residues of the CDR3.     

Solution structure of the Grb2 SH2 domain complexed with a high-affinity inhibitor

The solution structure of the growth factor receptor-bound protein 2 (Grb2) SH2 domain complexed with a high-affinity inhibitor containing a non-phosphorus phosphate mimetic within a macrocyclic platform was determined by nuclear magnetic resonance (NMR) spectroscopy. Unambiguous assignments of the bound inhibitor and intermolecular NOEs between the Grb2 SH2 domain and the inhibitor was accomplished using perdeuterated Grb2 SH2 protein. The well-defined solution structure of the complex was obtained and compared to those by X-ray crystallography. Since the crystal structure of the Grb2 SH2 domain formed a domain-swapped dimer and several inhibitors were bound to a hinge region, there were appreciable differences between the solution and crystal structures. Based on the binding interactions between the inhibitor and the Grb2 SH2 domain in solution, we proposed a design of second-generation inhibitors that could be expected to have higher affinity.     

Hydrogen-1, carbon-13, and nitrogen-15 NMR spectroscopy of Anabaena 7120 flavodoxin: assignment of beta-sheet and flavin binding site resonances and analysis of protein-flavin interactions

Sequence-specific 1H and 13C NMR assignments have been made for residues that form the five-stranded parallel beta-sheet and the flavin mononucleotide (FMN) binding site of oxidized Anabaena 7120 flavodoxin. Interstrand nuclear Overhauser enhancements (NOEs) indicate that the beta-sheet arrangement is similar to that observed in the crystal structure of the 70% homologous long-chain flavodoxin from Anacystis nidulans [Smith et al. (1983) J. Mol. Biol. 165, 737-755]. A total of 62 NOEs were identified: 8 between protons of bound FMN, 29 between protons of the protein in the flavin binding site, and 25 between protons of bound FMN and protons of the protein. These constraints were used to determine the localized solution structure of the FMN binding site. The electronic environment and conformation of the protein-bound flavin isoalloxazine ring were investigated by determining 13C chemical shifts, one-bond 13C-13C and 15N-1H coupling constants, and three-bond 13C-1H coupling constants. The carbonyl edge of the flavin ring was found to be slightly polarized. The xylene ring was found to be nonplanar. Tyrosine 94, located adjacent to the flavin isoalloxazine ring, was shown to have a hindered aromatic ring flip rate.     

Solution-state molecular structure of apo and oleate-liganded liver fatty acid-binding protein

Rat liver fatty acid-binding protein (LFABP) is distinctive among intracellular lipid-binding proteins (iLBPs): more than one molecule of long-chain fatty acid and a variety of diverse ligands can be bound within its large cavity, and in vitro lipid transfer to model membranes follows a mechanism that is diffusion-controlled rather than mediated by protein-membrane collisions. Because the apoprotein has proven resistant to crystallization, nuclear magnetic resonance spectroscopy offers a unique route to functionally informative comparisons of molecular structure and dynamics for LFABP in free (apo) and liganded (holo) forms. We report herein the solution-state structures determined for apo-LFABP at pH 6.0 and for holoprotein liganded to two oleates at pH 7.0, as well as the structure of the complex including locations of the ligands. 1H, 13C, and 15N resonance assignments revealed very similar types and locations of secondary structural elements for apo- and holo-LFABP as judged from chemical shift indices. The solution-state tertiary structures of the proteins were derived with the CNS/ARIA computational protocol, using distance and angular restraints based on 1H-1H nuclear Overhauser effects (NOEs), hydrogen-bonding networks, 3J(HNHA) coupling constants, intermolecular NOEs, and residual dipolar (NH) couplings. The holo-LFABP solution-state conformation is in substantial agreement with a previously reported X-ray structure [Thompson, J., Winter, N., Terwey, D., Bratt, J., and Banaszak, L. (1997) The crystal structure of the liver fatty acid-binding protein. A complex with two bound oleates, J. Biol. Chem. 272, 7140-7150], including the typical beta-barrel capped by a helix-turn-helix portal. In the solution state, the internally bound oleate has the expected U-shaped conformation and is tethered electrostatically, but the extended portal ligand can adopt a range of conformations based on the computationally refined structures, in contrast to the single conformation observed in the crystal structure. The apo-LFABP also has a well-defined beta-barrel structural motif typical of other members of the iLBP protein family, but the portal region that is thought to facilitate ligand entry and exit exhibits conformational variability and an unusual "open cap" orientation with respect to the barrel. These structural results allow us to propose a model in which ligand binding to LFABP occurs through conformational fluctuations that adjust the helix-turn-helix motif to open or close the top of the beta-barrel, and solvent accessibility to the protein cavity favors diffusion-controlled ligand transport.     

Determining molecular binding sites on human serum albumin by displacement of oleic acid

An NMR method was developed for determining binding sites of small molecules on human serum albumin (HSA) by competitive displacement of (13)C-labeled oleic acid. This method is based on the observation that in the crystal structure of HSA complexed with oleic acid, two principal drug-binding sites, Sudlow's sites I (warfarin) and II (ibuprofen), are also occupied by fatty acids. In two-dimensional [(1)H,(13)C]heteronuclear single quantum coherence NMR spectra, seven distinct resonances were observed for the (13)C-methyl-labeled oleic acid as a result of its binding to HSA. Resonances corresponding to the major drug-binding sites were identified through competitive displacement of molecules that bind specifically to each site. Thus, binding of molecules to these sites can be followed by their displacement of oleic acids. Furthermore, the amount of bound ligand at each site can be determined from changes in resonance intensities. For molecules containing fluorine, binding results were further validated by direct observations of the bound ligands using (19)F NMR. Identifying the binding sites for drug molecules on HSA can aid in determining the structure-activity relationship of albumin binding and assist in the design of molecules with altered albumin binding.     

Dissection of the critical binding determinants of cellular retinoic acid binding protein II by mutagenesis and fluorescence binding assay

The binding of retinoic acid to mutants of Cellular Retinoic Acid Binding Protein II (CRABPII) was evaluated to better understand the importance of the direct protein/ligand interactions. The important role of Arg111 for the correct structure and function of the protein was verified and other residues that directly affect retinoic acid binding have been identified. Furthermore, retinoic acid binding to CRABPII mutants that lack all previously identified interacting amino acids was rescued by providing a carboxylic acid dimer partner in the form of a Glu residue. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

Analysis of the binding of 1,3-diacetylchloramphenicol to chloramphenicol acetyltransferase by isotope-edited 1H NMR and site-directed mutagenesis.

The binary complex of diacetylchloramphenicol and chloramphenicol acetyltransferase (CAT) has been studied by a combination of isotope-edited 1H NMR spectroscopy and site-directed mutagenesis. One-dimensional HMQC spectra of the complex between 1,3-[2-13C]diacetylchloramphenicol and the type III natural variant of CAT revealed the two methyl 1H signals arising from each 13C-labeled carbon atom in the acetyl groups of the bound ligand. Slow hydrolysis of the 3-acetyl group by the enzyme precluded further analysis of this binary complex. It was possible to slow down the rate of hydrolysis by use of the catalytically defective S148A mutant of CATIII (Lewendon et al., 1990); in the complex of diacetylchloramphenicol with S148A CATIII, the chemical shifts of the acetyl groups of the bound ligand were the same as in the wild-type complex. The acetyl signals were individually assigned by repeating the experiment using 1-[2-13C],3-[2-12C]diacetylchloramphenicol, where only one signal from the bound ligand was observed. A two-dimensional 1H, 1H NOESY experiment, with 13C(omega 2) half-filter, on the 1,3-[2-13C]diacetylchloramphenicol/S148A CATIII complex showed a number of intermolecular NOEs from each methyl group in the ligand to residues in the chloramphenicol binding site. The 3-acetyl group showed strong NOEs to two aromatic signals which were selected for assignment. The possibility that the NOEs originated from the aromatic protons of diacetylchloramphenicol itself was eliminated by assignment of the signals from enzyme-bound diacetylchloramphenicol and chloramphenicol using perdeuterated CATIII. Examination of the X-ray crystal structure of the chloramphenicol/CATIII binary complex indicated four plausible candidate aromatic residues: Y25, F33, F103, and F158.(ABSTRACT TRUNCATED AT 250 WORDS)     

Unbinding of retinoic acid from the retinoic acid receptor by random expulsion molecular dynamics

Unbinding pathways of retinoic acid (RA) bound to retinoic acid receptor (RAR) have been explored by the random expulsion molecular dynamics (REMD) method. Our results show that RA may exit the binding site of RAR through flexible regions close to the H1-H3 loop and beta-sheets, without displacing H12 from its agonist position. This result may explain kinetic differences between agonist and antagonist ligands observed for other nuclear receptors. The extended and flexible structure of RA initiated a methodological study in a simplified two-dimensional model system. The REMD force should in general be distributed to all atoms of the ligand to obtain the most unbiased results, but for a ligand which is tightly bound in the binding pocket through a strong electrostatic interaction, application of the REMD force on a single atom is preferred.     

Structure and dynamics in solution of the complex of Lactobacillus casei dihydrofolate reductase with the new lipophilic antifolate drug trimetrexate

We have determined the three-dimensional solution structure of the complex of Lactobacillus casei dihydrofolate reductase and the anticancer drug trimetrexate. Two thousand seventy distance, 345 dihedral angle, and 144 hydrogen bond restraints were obtained from analysis of multidimensional NMR spectra recorded for complexes containing 15N-labeled protein. Simulated annealing calculations produced a family of 22 structures fully consistent with the constraints. Several intermolecular protein-ligand NOEs were obtained by using a novel approach monitoring temperature effects of NOE signals resulting from dynamic processes in the bound ligand. At low temperature (5 degrees C) the trimethoxy ring of bound trimetrexate is flipping sufficiently slowly to give narrow signals in slow exchange, which give good NOE cross peaks. At higher temperature these broaden and their NOE cross peaks disappear thus allowing the signals in the lower-temperature spectrum to be identified as NOEs involving ligand protons. The binding site for trimetrexate is well defined and this was compared with the binding sites in related complexes formed with methotrexate and trimethoprim. No major conformational differences were detected between the different complexes. The 2,4-diaminopyrimidine-containing moieties in the three drugs bind essentially in the same binding pocket and the remaining parts of their molecules adapt their conformations such that they can make effective van der Waals interactions with essentially the same set of hydrophobic amino acids, the side-chain orientations and local conformations of which are not greatly changed in the different complexes (similar chi1 and chi2 values).     

NMR studies of the conformations of leghemoglobins from soybean and lupin

Phase-sensitive two-dimensional NMR methods have been used to obtain extensive proton resonance assignments for the carbon monoxide complexes of lupin leghemoglobins I and II and soybean leghemoglobin a. The assigned resonances provide information on the solution conformations of the proteins, particularly in the vicinity of the heme. The structure of the CO complex of lupin leghemoglobin II in solution is compared with the X-ray crystal structure of the cyanide complex by comparison of observed and calculated ring current shifts. The structures are generally very similar but significant differences are observed for the ligand contact residues, Phe30, His63 and Val67, and for the proximal His97 ligand. Certain residues are disordered and adopt two interconverting conformations in lupin leghemoglobin II in solution. The proximal heme pocket structure is closely conserved in the lupin leghemoglobins I and II but small differences in conformation in the distal heme pocket are apparent. Larger conformational differences are observed when comparisons are made with the CO complex of soybean leghemoglobin. Altered protein-heme packing is indicated on the proximal side of the heme and some conformational differences are evident in the distal heme pocket. The small conformational differences between the three leghemoglobins probably contribute to the known differences in their O2 and CO association and dissociation kinetics. The heme pocket conformations of the three leghemoglobins are more closely related to each other than to sperm whale myoglobin. The most notable differences between the leghemoglobins and myoglobin are: (a) reduced steric crowding of the ligand binding site in the leghemoglobins, (b) different orientations of the distal histidine, and (c) small but significant differences in proximal histidine coordination geometry. These changes probably contribute to the large differences in ligand binding kinetics between the leghemoglobins and myoglobin.     

1H NMR studies of porcine calbindin D9k in solution: sequential resonance assignment, secondary structure, and global fold

The 1H nuclear magnetic resonance (NMR) spectrum of Ca2+-saturated porcine calbindin D9k (78 amino acids, Mr 8800) has been assigned. Greater than 98% of the 1H resonances, including spin systems for each amino acid residue, have been identified by using an approach that integrates data from a wide range of two-dimensional scalar correlated NMR experiments [Chazin, Rance, & Wright (1988) J. Mol. Biol. 202, 603-626]. Due to the limited quantity of sample and conformational heterogeneity of the protein, two-dimensional nuclear Overhauser effect (NOE) experiments also played an essential role in the identification of spin systems. On the basis of the pattern of scalar connectivities, 43 of the 78 spin systems could be directly assigned to the appropriate residue type. This provided an ample basis for obtaining the sequence-specific resonance assignments. The elements of secondary structure are identified from sequential and medium-range NOEs, values of 3JNH alpha, and the location of slowly exchanging backbone amide protons. Four well-defined helices and a mini beta-sheet between the two calcium binding loops are present in solution. These elements of secondary structure and a few key long-range NOEs provided sufficient information to define the global fold of the protein in solution. Generally good agreement is found between the crystal structure of the minor A form of bovine calbindin D9k and the solution structure of intact porcine calbindin D9k. The only significant difference is a short one-turn helix in the loop between helices II and III in the bovine crystal structure, which is clearly absent in the porcine solution structure.     

Sequence-specific 1H NMR assignments and secondary structure in solution of Escherichia coli trp repressor

Sequence-specific 1H NMR assignments are reported for the active L-tryptophan-bound form of Escherichia coli trp repressor. The repressor is a symmetric dimer of 107 residues per monomer; thus at 25 kDa, this is the largest protein for which such detailed sequence-specific assignments have been made. At this molecular mass the broad line widths of the NMR resonances preclude the use of assignment methods based on 1H-1H scalar coupling. Our assignment strategy centers on two-dimensional nuclear Overhauser spectroscopy (NOESY) of a series of selectively deuterated repressor analogues. A new methodology was developed for analysis of the spectra on the basis of the effects of selective deuteration on cross-peak intensities in the NOESY spectra. A total of 90% of the backbone amide protons have been assigned, and 70% of the alpha and side-chain proton resonances are assigned. The local secondary structure was calculated from sequential and medium-range backbone NOEs with the double-iterated Kalman filter method [Altman, R. B., & Jardetzky, O. (1989) Methods Enzymol. 177, 218-246]. The secondary structure agrees with that of the crystal structure [Schevitz, R., Otwinowski, Z., Joachimiak, A., Lawson, C. L., & Sigler, P. B. (1985) Nature 317, 782], except that the solution state is somewhat more disordered in the DNA binding region and in the N-terminal region of the first alpha-helix. Since the repressor is a symmetric dimer, long-range intersubunit NOEs were distinguished from intrasubunit interactions by formation of heterodimers between two appropriate selectively deuterated proteins and comparison of the resulting NOESY spectrum with that of each selectively deuterated homodimer. Thus, from spectra of three heterodimers, long-range NOEs between eight pairs of residues were identified as intersubunit NOEs, and two additional long-range intrasubunits NOEs were assigned.     

Assignment of the 1H NMR resonances of protein residues in close proximity to the heme of the nitrophorins: similarities and differences among the four proteins from the saliva of the adult blood-sucking insect Rhodnius prolixus

The nuclear Overhauser effects (NOEs) observed between heme substituent protons and a small number of nearby protein side chain protons in the water-elimination Fourier transform NOE spectroscopy (WEFT-NOESY) spectra of high- and low-spin wild-type nitrophorin (NP) 2 and its ligand complexes have been analyzed and compared with those observed for the same complexes of wild-type NP3. These assignments were made on naturally abundant isotope samples, with the most useful protein side chains being those of Ile120, Leu122, and Leu132 for NP2 and NP3, and Thr121, Leu123, and Leu133 for NP1 and NP4. It is found that the NOEs observed are identical, with extremely similar protein side chain proton chemical shifts. This is strong evidence that the structure of NP3, for which no X-ray crystal structures are available, is essentially identical to that of NP2, at least near the heme binding pocket. Similarly, the NOEs observed between heme substituents and protein side chains for NP1 and NP4 also indicate that the structures of the protein having both A and B heme orientations are very similar to each other, as well as to the proteins with major B heme orientation of NP2 and NP3. These A and B connectivities can be seen, even though the two heme orientations have similar populations in NP1 and NP4, which complicates the analysis of the NOESY spectra. The histamine complex of wild-type NP2 shows significant shifts of the Leu132 side chain protons relative to all other ligand complexes of NP1-NP4 because of the perturbation of the structure near Leu132 caused by the histamine's side chain ammonium hydrogen bond to the Asp29 side chain carboxylate.     

Metal-ion binding and metal-ion induced folding of the adenine-sensing riboswitch aptamer domain

Divalent cations are important in the folding and stabilization of complex RNA structures. The adenine-sensing riboswitch controls the expression of mRNAs for proteins involved in purine metabolism by directly sensing intracellular adenine levels. Adenine binds with high affinity and specificity to the ligand binding or aptamer domain of the adenine-sensing riboswitch. The X-ray structure of this domain in complex with adenine revealed an intricate RNA-fold consisting of a three-helix junction stabilized by long-range base-pairing interactions and identified five binding sites for hexahydrated Mg²⁺-ions. Furthermore, a role for Mg²⁺-ions in the ligand-induced folding of this RNA was suggested. Here, we describe the interaction of divalent cations with the RNA–adenine complex in solution as studied by high-resolution NMR spectroscopy. Paramagnetic line broadening, chemical shift mapping and intermolecular nuclear Overhauser effects (NOEs) indicate the presence of at least three binding sites for divalent cations. Two of them are similar to those in the X-ray structure. The third site, which is important for the folding of this RNA, has not been observed previously. The ligand-free state of the RNA is conformationally heterogeneous and contains base-pairing patterns detrimental to ligand binding in the absence of Mg²⁺, but becomes partially pre-organized for ligand binding in the presence of Mg²⁺. Compared to the highly similar guanine-sensing riboswitch, the folding pathway for the adenine-sensing riboswitch aptamer domain is more complex and the influence of Mg²⁺ is more pronounced.     

A "structural" water molecule in the family of fatty acid binding proteins

A single water molecule (w135), buried within the structure of rat intestinal fatty acid binding protein (I-FABP), is investigated by NMR, molecular dynamics simulations, and analysis of known crystal structures. An ordered water molecule was found in structurally analogous position in 24 crystal structures of nine different members of the family of fatty acid binding proteins. There is a remarkable conservation of the local structure near the w135 binding site among different proteins from this family. NMR cross-relaxation measurements imply that w135 is present in the I-FABP:ANS (1-sulfonato-8-(1')anilinonaphthalene) complex in solution with the residence time of >300 ps. Mean-square positional fluctuations of w135 oxygen observed in MD simulations (0.18 and 0.13 A2) are comparable in magnitude to fluctuations exhibited by the backbone atoms and result from highly constrained binding pocket as revealed by Voronoi volumes (averages of 27.0 +/- 1.8 A3 and 24.7 +/- 2.2 A3 for the two simulations). Escape of w135 from its binding pocket was observed only in one MD simulation. The escape process was initiated by interactions with external water molecules and was accompanied by large deformations in beta-strands D and E. Immediately before the release, w135 assumed three distinct states that differ in hydrogen bonding topology and persisted for about 15 ps each. Computer simulations suggest that escape of w135 from the I-FABP matrix is primarily determined by conformational fluctuations of the protein backbone and interactions with external water molecules.     

High resolution crystal structures of unliganded and liganded human liver ACBP reveal a new mode of binding for the acyl-CoA ligand

The acyl-CoA binding protein (ACBP) is essential for the fatty acid metabolism, membrane structure, membrane fusion, and ceramide synthesis. Here high resolution crystal structures of human cytosolic liver ACBP, unliganded and liganded with a physiological ligand, myristoyl-CoA are described. The binding of the acyl-CoA molecule induces only few structural differences near the binding pocket. The crystal form of the liganded ACBP, which has two ACBP molecules in the asymmetric unit, shows that in human ACBP the same acyl-CoA binding pocket is present as previously described for the bovine and Plasmodium falciparum ACBP and the mode of binding of the 3'-phosphate-AMP moiety is conserved. Unexpectedly, in one of the acyl-CoA binding pockets the acyl moiety is bound in a reversed mode as compared with the bovine and P. falciparum structures. In this binding mode, the myristoyl-CoA molecule is fully ordered and bound across the two ACBP molecules of the crystallographic asymmetric unit: the 3'-phosphate-AMP moiety is bound in the binding pocket of one ACBP molecule and the acyl chain is bound in the pocket of the other ACBP molecule. The remaining binding pocket cavities of these two ACBP molecules are filled by other ligand fragments. This novel binding mode shows that the acyl moiety can flip out of its classical binding pocket and bind elsewhere, suggesting a mechanism for the acyl-CoA transfer between ACBP and the active site of a target enzyme. This mechanism is of possible relevance for the in vivo function of ACBP.     

Crystal structure of GGA2 VHS domain and its implication in plasticity in the ligand binding pocket

Golgi-localized, gamma-ear-containing, ARF binding (GGA) proteins regulate intracellular vesicle transport by recognizing sorting signals on the cargo surface in the initial step of the budding process. The VHS (VPS27, Hrs, and STAM) domain of GGA binds with the signal peptides. Here, a crystal structure of the VHS domain of GGA2 is reported at 2.2 A resolution, which permits a direct comparison with that of homologous proteins, GGA1 and GGA3. Significant structural difference is present in the loop between helices 6 and 7, which forms part of the ligand binding pocket. Intrinsic fluorescence spectroscopic study indicates that this loop undergoes a conformational change upon ligand binding. Thus, the current structure suggests that a conformational change induced by ligand binding occurs in this part of the ligand pocket.     

Direct photoaffinity labeling of cellular retinoic acid-binding protein I (CRABP-I) with all-trans-retinoic acid: identification of amino acids in the ligand binding site

Cellular retinoic acid-binding proteins I and II (CRABP-I and -II, respectively) are transport proteins for all-trans-retinoic acid (RA), an active metabolite of vitamin A (retinol), and have been reported to be directly involved in the metabolism of RA. In this study, direct photoaffinity labeling with [11,12-(3)H]RA was used to identify amino acids comprising the ligand binding site of CRABP-I. Photoaffinity labeling of CRABP-I with [(3)H]RA was light- and concentration-dependent and was protected by unlabeled RA and various retinoids, indicating that the labeling was directed to the RA-binding site. Photolabeled CRABP-I was hydrolyzed with endoproteinase Lys-C to yield radioactive peptides, which were separated by reversed-phase HPLC for analysis by Edman degradation peptide sequencing. This method identified five modified amino acids from five separate HPLC fractions: Trp7, Lys20, Arg29, Lys38, and Trp109. All five amino acids are located within one side of the "barrel" structure in the area indicated by the reported crystal structure as the ligand binding site. This is the first direct identification of specific amino acids in the RA-binding site of CRABPs by photoaffinity labeling. These results provide significant information about the ligand binding site of the CRABP-I molecule in solution.