Binding to δ and μ opioid receptors by deltorphin I/II analogues modified at the Phe and Asp/Glu side chains: a report of 32 new analogues and a QSAR study

The synthesis and binding affinities of 32 X³Gly⁴ dual-substitution analogues of the natural opioid heptapeptides deltorphin I and II are reported. A multiple regression QSAR analysis was performed using those results along with literature data for the X³Asp⁴ and Phe³X⁴ side chain analogues. Fitting to a three-term potential well model with hydrophobic and van der Waals attraction terms and a steric repulsion term indicates that the δ and μ receptor sites for binding the residue three side chain are similar, and that the binding interaction is primarily van der Waals and secondarily hydrophobic. Further analysis indicates that both sites are more constrained with respect to side chain length than width or thickness, and the μ site appears to be somewhat larger. A binding model consistent with these findings pictures the native third residues Phe ring laying on a step notched out of the receptor surface, pointing toward the back (riser) of the step, and sandwiched between the receptor and ligand. However, the binding sites for the residue four side chains are quite different on δ and μ receptors. Binding to the δ site appears to involve both electrostatic attraction (probably to a partial positive charge) and van der Waals attraction, but not necessarily hydrogen bonding, and more constraint with respect to side chain length than width or thickness. In contrast, there is no evidence for any kind of binding attraction between the side chain of residue four and the μ site, which acts more as steric repulsion site, as though the space that is a pocket on the δ receptor is filled in on the μ receptor. A regression model based only on steric repulsion by van der Waals bulk and/or the effective bulk of a hydration layer accounts for over 80% of the residue four related variation in μ affinity.

Abstract

Thirty-two new X³Gly⁴ analogues of deltrophin I/II opioid peptides are described. A QSAR study of the X³Gly⁴, X³Asp⁴, and Phe³X⁴ analogue series using a potential well model reveals the roles of hydrophobic, van der Waals, electrostatic, hydrogen bonding and steric interactions in δ and μ receptor binding of X³ and X⁴ side chains.     

Refinement of the structure of Escherichia coli-derived rat intestinal fatty acid binding protein with bound oleate to 1.75-A resolution. Correlation with the structures of the apoprotein and the protein with bound palmitate.

The structure of rat intestinal fatty acid binding protein (I-FABP) with bound oleate (C18:1) has been refined with x-ray diffraction data to a resolution of 1.75 A. The protein contains 10 anti-parallel beta strands composed of 99 residues and 2 short helices of 14 residues. Oleate is located in the interior of the protein in a bent conformation with C1-C12 more ordered than C13-C18. Two of the eight ordered waters in I-FABP:oleate are part of a hydrogen bond network that includes the carboxylate of oleate, the guanidinium group of Arg106, the nitrogen of the indole group of Trp82, and the side chain of Gln115. Most of the methylenes of bound oleate reside in a crevice formed by hydrophobic and aromatic side chains. Tyr70 and Tyr117 envelop the acyl chain from C3 to C8 forming contacts with both the convex and concave faces of its van der Waals surface. The hydroxyls of each phenolic side chain hydrogen bond to ordered water molecules. Two ordered waters make van der Waals contact with the concave face of the bound fatty acid. The omega-terminal methyl of oleate is oriented so that it points toward the center of the benzene of Phe55 allowing it to form van der Waals interactions with its component methylenes. Comparison of the structure of I-FABP:oleate with a recently refined 1.19-A model of apoI-FABP and an earlier 2.0-A model of I-FABP:palmitate revealed a remarkable degree of similarity in the positions of their main chain and side chain atoms and in the conformations of the bound oleate and palmitate. The principal differences were confined to a few discrete regions of the protein. The helical domain, the type I turn between beta strands C and D, and the ring of Phe55 together form a solvent-accessible portal to the interior of the protein. They are repositioned in I-FABP:oleate (and I-FABP:palmitate) so that the binding cavity is even more accessible to solvent and its volume is increased. The side chain of Phe55 which shows discrete disorder in the apoprotein functions as an omega-terminal "sensing device": moving progressively outward toward the surface as the chain length of the bound fatty acid increases by 2 methylenes. Tyr70 and Tyr117 which also show discrete disorder in the apoprotein structure due to rotation around their C alpha-C beta bonds, are stabilized in a single, well ordered position in the holoproteins.(ABSTRACT TRUNCATED AT 400 WORDS)     

High resolution X-ray studies of mammalian intestinal and muscle fatty acid-binding proteins provide an opportunity for defining the chemical nature of fatty acid: protein interactions

The structure ofE. coli-derived rat intestinal fatty acid-binding protein has recently been refined to 1.2 Å without bound fatty acid and to 2.0 Å and 1.75 Å with bound hexadecanoate (palmitate) and 9Z-octadecenoate (oleate), respectively. The structure ofE. coli-derived human muscle fatty acid-binding protein has also been solved to 2.1 Å with a C16 bacterial fatty acid. Both proteins contain 10 anti-parallel -strands in a+1, +1, +1... motif. The strands are arranged in two -pleated sheets that are orthogonally oriented. In each case, the fatty acid is enclosed by the -sheets and is bound to the proteins by feeble forces. These feeble forces consist of (i) a hydrogen bonding network between the fatty acid's carboxylate group, ordered solvent, and side chains of polar/ionizable amino acid residues; (ii) van der Waals contacts between the methylene chain of the fatty acid and the side chain atoms of hydrophobic and aromatic residues; (iii) van der Waals interactions between the methyl and the component methenyls of the phenyl side chain of a Phe which serves as an adjustable terminal sensor situated over a surface opening or portal connecting interior and exterior solvent; and (iv) van der Waals contacts between methylenes of the alkyl chain and oxygens of ordered waters that have been located inside the binding cavity. These waters are positioned over one face of the ligand and are held in place by hydrogen bonding with one another and with the side chains of protein's polar and ionizable residues. Binding of the fatty acid ligand is associated with minimal adjustments of the positions of main chain or side chain atoms. However, acquisition of ligand is associated with removal of ordered interior solvent suggesting that the free energy of dehydration of the binding site may be as important for the energy of the binding reaction as the free energy of stabilization of the fatty acid: protein complex.     

Analysis of the interactions of sulfur‐containing amino acids in membrane proteins

The interactions of Met and Cys with other amino acid side chains have received little attention, in contrast to aromatic–aromatic, aromatic–aliphatic or/and aliphatic–aliphatic interactions. Precisely, these are the only amino acids that contain a sulfur atom, which is highly polarizable and, thus, likely to participate in strong Van der Waals interactions. Analysis of the interactions present in membrane protein crystal structures, together with the characterization of their strength in small‐molecule model systems at the ab‐initio level, predicts that Met–Met interactions are stronger than Met–Cys ≈ Met–Phe ≈ Cys–Phe interactions, stronger than Phe–Phe ≈ Phe–Leu interactions, stronger than the Met–Leu interaction, and stronger than Leu–Leu ≈ Cys–Leu interactions. These results show that sulfur‐containing amino acids form stronger interactions than aromatic or aliphatic amino acids. Thus, these amino acids may provide additional driving forces for maintaining the 3D structure of membrane proteins and may provide functional specificity.     

Carbonic anhydrase inhibitors: Valdecoxib binds to a different active site region of the human isoform II as compared to the structurally related cyclooxygenase II "selective" inhibitor celecoxib

The high resolution X-ray crystal structure of the adduct of human carbonic anhydrase (CA, EC 4.2.1.1) isoform II (hCA II) with the clinically used painkiller valdecoxib, acting as a potent CA II and cyclooxygenase-2 (COX-2) inhibitor, is reported. The ionized sulfonamide moiety of valdecoxib is coordinated to the catalytic Zn(II) ion with a tetrahedral geometry. The phenyl-isoxazole moiety of the inhibitor fills the active site channel and interacts with the side chains of Gln92, Val121, Leu198, Thr200, and Pro202. Its 3-phenyl group is located into a hydrophobic pocket, simultaneously establishing van der Waals interactions with the aliphatic side chain of various hydrophobic residues (Val135, Ile91, Val121, Leu198, and Leu141) and a strong offset face-to-face stacking interaction with the aromatic ring of Phe131 (the chi1 angle of which is rotated about 90 degrees with respect to what was observed in the structure of the native enzyme and those of other sulfonamide complexes). Celecoxib, a structurally related COX-2 inhibitor for which the X-ray crystal structure was reported earlier, binds in a completely different manner to hCA II as compared to valdecoxib. Celecoxib completely fills the entire CA II active site, with its trifluoromethyl group in the hydrophobic part of the active site and the p-tolyl moiety in the hydrophilic one, not establishing any interaction with Phe131. In contrast to celecoxib, valdecoxib was rotated about 90 degrees around the chemical bond connecting the benzensulfonamide and the substituted isoxazole ring allowing for these multiple favorable interactions. These different binding modes allow for the further drug design of various CA inhibitors belonging to the benzenesulfonamide class.     

Side-chain flexibility in proteins upon ligand binding

Ligand binding may involve a wide range of structural changes in the receptor protein, from hinge movement of entire domains to small side-chain rearrangements in the binding pocket residues. The analysis of side chain flexibility gives insights valuable to improve docking algorithms and can provide an index of amino-acid side-chain flexibility potentially useful in molecular biology and protein engineering studies. In this study we analyzed side-chain rearrangements upon ligand binding. We constructed two non-redundant databases (980 and 353 entries) of "paired" protein structures in complexed (holo-protein) and uncomplexed (apo-protein) forms from the PDB macromolecular structural database. The number and identity of binding pocket residues that undergo side-chain conformational changes were determined. We show that, in general, only a small number of residues in the pocket undergo such changes (e.g., approximately 85% of cases show changes in three residues or less). The flexibility scale has the following order: Lys > Arg, Gln, Met > Glu, Ile, Leu > Asn, Thr, Val, Tyr, Ser, His, Asp > Cys, Trp, Phe; thus, Lys side chains in binding pockets flex 25 times more often then do the Phe side chains. Normalizing for the number of flexible dihedral bonds in each amino acid attenuates the scale somewhat, however, the clear trend of large, polar amino acids being more flexible in the pocket than aromatic ones remains. We found no correlation between backbone movement of a residue upon ligand binding and the flexibility of its side chain. These results are relevant to 1. Reduction of search space in docking algorithms by inclusion of side-chain flexibility for a limited number of binding pocket residues; and 2. Utilization of the amino acid flexibility scale in protein engineering studies to alter the flexibility of binding pockets.     

The impact of <Emphasis Type="Italic">Trichoderma reesei</Emphasis> Cel7A carbohydrate binding domain mutations on its binding to a cellulose surface: a molecular dynamics free energy study

A critical role of the Family 7 cellobiohydrolase (Cel7A) carbohydrate binding domain (CBD) is to bind to a cellulose surface and increase the enzyme concentration on the surface. Several residues of Trichoderma reesei Cel7A CBD, including Y5, N29, Y31, Y32 and Q34, contribute to cellulose binding, as revealed by early experimental studies. To investigate the interactions between these important residues and cellulose, we applied a thermodynamic integration method to calculate the cellulose–Cel7A CBD binding free energy changes caused by Y5A, N29A, Y31A, Y32A and Q34A mutations. The experimental binding trend was successfully predicted, proving the effectiveness of the complex model. For the two polar residue mutants N29A and Q34A, the changes in the electrostatics are comparable to those of van der Waals, while for three Y to A mutants, the free energy differences mainly come from van der Waals interactions. However, in both cases, the electrostatics dominates the interactions between individual residues and cellulose. The side chains of these residues are rigidified after the complex is formed. The binding free energy changes for the two mutants Y5W and Y31W were also determined, and for these the van der Waals interaction was strengthened but the electrostatics was weakened.     

The 2.3-A resolution structure of the maltose- or maltodextrin-binding protein, a primary receptor of bacterial active transport and chemotaxis

The three-dimensional structure of the maltose- or maltodextrin-binding protein (Mr = 40,622) with bound maltose has been obtained by crystallographic analysis at 2.8-A resolution. The structure, which has been partially refined at 2.3 A, is ellipsoidal with overall dimensions of 30 x 40 x 65 A and divided into two distinct globular domains by a deep groove. Although each domain is built from two peptide segments from the amino- and carboxyl-terminal halves, both domains exhibit similar supersecondary structure, consisting of a central beta-pleated sheet flanked on both sides with two or three parallel alpha-helices. The groove, which has a depth of 18 A and a base of about 9 x 18 A, contains the maltodextrin-binding site. We have previously observed the same general features in the well-refined structures of six other periplasmic receptors with specificities for L-arabinose, D-galactose/D-glucose, sulfate, phosphate, leucine/isoleucine/valine, and leucine. The bound maltose is buried in the groove and almost completely inaccessible to the bulk solvent. The groove is heavily populated by polar and aromatic groups many of which are involved in extensive hydrogen-bonding and van der Waals interactions with the maltose. All the disaccharide hydroxyl groups, which form a peripheral polar surface approximately in the plane of the sugar rings, are tied in a total of 11 direct hydrogen bonds with six charged side chains, one Trp side chain, and one peptide backbone NH, and five indirect hydrogen bonds via water molecules. The maltose is wedged between four aromatic side chains. The resulting stacking of these aromatic residues on the faces of the glucosyl units provides a majority of the van der Waals contacts in the complex. The nonreducing glucosyl unit of the maltose is involved in approximately twice as many hydrogen bonds and van der Waals contacts as the glucosyl unit at the reducing end. The binding protein-maltose complex shows the best example of the extensive use of polar and aromatic residues in binding oligosaccharides. The tertiary structure of the maltodextrin-binding protein, along with the results of genetic studies by a number of investigators, has also enabled us for the first time to map the different regions on the surface of the protein involved in the interactions with the membrane-bound protein components necessary for transport of and chemotaxis toward maltodextrins. These sites permit distinction of the "open cleft" (without bound sugar) and closed (with bound sugar) conformations of the binding protein by the chemotactic signal transducer with which the maltodextrin-binding protein interacts.(ABSTRACT TRUNCATED AT 250 WORDS)     

Thermodynamics of interactions between amino acid side chains: experimental differentiation of aromatic-aromatic, aromatic-aliphatic, and aliphatic-aliphatic side-chain interactions in water

A stationary phase for high-pressure liquid chromatography has been prepared by derivatizing microparticulate silica gel with functionality mimicking the side chain of isoleucine. The chromatographic retentions of a series of hydrophobic and amphiphilic amino acid analytes on this stationary phase (Ile MSP) using an aqueous mobile phase were measured as a function of temperature from 273 K to 323 K. Observed temperature dependencies are consistent with a constant change in heat capacity, DeltaC degrees P, upon binding of the analyte to the stationary phase. The curvatures of plots of retention data versus temperature (related to the magnitude of DeltaC degrees P) are distinctly different for retention of aromatic and aliphatic analytes, with retention of aliphatic analytes Val, Ile, and Leu exhibiting the characteristic signature of the hydrophobic effect, i.e., a large negative DeltaC degrees P upon desolvation from water and a maximum of retention around room temperature. Retention of aromatic analytes (Trp, Phe, and Tyr) involves smaller heat capacity changes and pronounced negative enthalpies of interaction with the stationary phase. Estimates of DeltaC degrees P for the interactions of analyte side chains with the Ile side chain were obtained by fitting the temperature dependence of retention to an expression derived from thermodynamic considerations and chromatographic theory. Similar estimates were made for interactions with the Phe side chain, using previously published data for a phenylalanine mimic stationary phase (Phe MSP) (. Protein Sci. 1:786-795). As with the Ile MSP, the retentions of aliphatic analytes show temperature dependencies markedly different from those of aromatic analytes. Data from both phases indicate that a realistic differentiation can be made between the interactions of various types of amino acid side chains tested (i.e., aliphatic/aliphatic, aliphatic/aromatic, and aromatic/aromatic) by comparison of the corresponding thermodynamic functions for pairwise interactions. The retention of leucine on the Phe MSP and that of phenylalanine on the Ile MSP showed similar DeltaC degrees P values, suggesting that the aromatic-aliphatic interaction is reasonably independent of the residue attached to the stationary phase. This result is consistent with a one-to-one interaction and suggests a simple way to estimate the column-dependent phase factor, making it possible to compare entropies and free energies of interaction obtained using different MSPs. The possibilities for using MSP-derived interaction potentials in folding simulations are discussed.     

A novel site-directed mutant of myoglobin with an unusually high O2 affinity and low autooxidation rate.

Mutants of sperm whale myoglobin were constructed at position 29 (B10 in helix notation) to examine the effects of distal pocket size on the rates of ligand binding and autooxidation. Leu29 was replaced with Ala, Val, and Phe using the synthetic gene and Escherichia coli expression system of Springer and Sligar (Springer, B. A., and Sligar, S. G. (1987) Proc. Natl. Acad. Sci. U. S. A. 84, 8961-8965). Structures of the ferric forms of Val29 and Phe29, and the oxy form of Phe29 myoglobin were determined to 1.7 A by x-ray crystallography. The ferric mutant proteins are remarkably isomorphous with the wild type protein except in the immediate vicinity of residue 29. Thus, the protein structure in the distal pocket of myoglobin can accommodate either a large "hole" (i.e. Ala or Val) or a large side chain (i.e. Phe) at position 29 without perturbation of tertiary structure. Phe29 oxymyoglobin is also identical to the native oxy protein in terms of overall structure and interactions between the bound O2 and His64, Val68, Phe43, and Ile107. The distance between the nearest side chain atom of residue 29 and the second atom of the bound oxygen molecule is 3.2 A in the Phe29 protein and 4.9 A in native myoglobin. The equilibrium constants for O2 binding to Ala29, Val29, and Leu29 (native) myoglobin are the same, approximately 1.0 x 10(6) M-1 at 20 degrees C, whereas that for the Phe29 protein is markedly greater, 15 x 10(6) M-1. This increase in affinity is due primarily to a 10-fold decrease in the O2 dissociation rate constant for the Phe29 mutant and appears to be the result of stabilizing interactions between the negative portion of the bound O2 dipole and the partially positive edge of the phenyl ring. Increasing the size of residue 29 causes large decreases in the rate of autooxidation of myoglobin: k(ox) = 0.24, 0.23, 0.055, and 0.005 h-1 for Ala29, Val29, Leu29 (native), and Phe29 myoglobin, respectively, in air at 37 degrees C. Thus, the Leu29----Phe mutation produces a reduced protein that is remarkably stable and is expressed in E. coli as 100% MbO2. The selective pressure to conserve Leu29 at the B10 position probably represents a compromise between reducing the rate of autooxidation and maintaining a large enough O2 dissociation rate constant to allow rapid oxygen release during respiration.     

Replacements in a conserved leucine cluster in the hydrophobic heme pocket of cytochrome c.

A cluster of highly conserved leucine side chains from residues 9, 68, 85, 94, and 98 is located in the hydrophobic heme pocket of cytochrome c. The contributions of two of these, Leu 85 and Leu 94, have been studied using a protein structure-function-mutagenesis approach to probe their roles in the maintenance of overall structural integrity and electron transfer activity. Structural studies of the L85C, L85F, L85M, and L94S mutant proteins show that, in each case, the overall fold of cytochrome c is retained, but that localized conformational shifts are required to accommodate the introduced side chains. In particular, the side chains of Cys 85 and Phe 85 form energetically favorable interactions with Phe 82, whereas Met 85 takes on a more remote conformation to prevent an unfavorable interaction with the phenyl ring of Phe 82. In the case of the L94S mutant protein, the new polar group introduced is found to form hydrogen bonds to nearby carbonyl groups. In all proteins with substitutions at Leu 85, the hydrophobic nature of the heme pocket is preserved and no significant decrease in heme reduction potential is observed. Despite earlier predictions that Leu 85 is an important determinant in cytochrome c electron transfer partner complexation, our studies show this is unlikely to be the case because the considerable surface contour perturbations made by substitutions at this residue do not correspondingly translate into significant changes in electron transfer rates. For the L94S mutant protein, the substitution of a polar hydroxyl group directly into the hydrophobic heme pocket has a larger effect on heme reduction potential, but this is mitigated by two factors. First, the side chain of Ser 94 is rotated away from the heme group and, second, the side chain of Leu 98 shifts into a portion of the new space available, partially shielding the heme group. The Leu 94 Ser substitution does not perturb the highly conserved interface formed by the nearly perpendicular packing of the N- and C-terminal helices of cytochrome c, ruling this out as the cause of this mutant protein becoming thermally labile and having a lower functional activity. Our results show these effects are most likely attributable to disruption of the heme pocket region. Much of the ability of cytochrome c to absorb the introduction of mutations at Leu 85 and Leu 94 appears to be a consequence of the conformational flexibility afforded by the leucine cluster in this region as well as the presence of a nearby internal cavity.(ABSTRACT TRUNCATED AT 400 WORDS)     

Role of amino acid hydrophobicity, aromaticity, and molecular volume on IAPP(20-29) amyloid self-assembly

Aromatic amino acids strongly promote cross-β amyloid formation; whether the amyloidogenicity of aromatic residues is due to high hydrophobicity and β-sheet propensity or formation of stabilizing π-π interactions has been debated. To clarify the role of aromatic residues on amyloid formation, the islet amyloid polypeptide 20-29 fragment [IAPP(20-29)], which contains a single aromatic residue (Phe 23), was adopted as a model. The side chain of residue 23 does not self-associate in cross-β fibrils of IAPP(20-29) (Nielsen et al., Angew Chem Int Ed 2009;48:2118-2121), allowing investigation of the amyloidogenicity of aromatic amino acids in a context where direct π-π interactions do not occur. We prepared variants of IAPP(20-29) in which Tyr, Leu, Phe, pentafluorophenylalanine (F5-Phe), Trp, cyclohexylalanine (Cha), α-naphthylalanine (1-Nap), or β-naphthylalanine (2-Nap) (in order of increasing peptide hydrophobicity) were incorporated at position 23 (SNNXGAILSS-NH2), and the kinetic and thermodynamic effects of these mutations on cross-β self-assembly were assessed. The Tyr, Leu, and Trp 23 variants failed to readily self-assemble at concentrations up to 1.5 mM, while the Cha 23 mutant fibrillized with attenuated kinetics and similar thermodynamic stability relative to the wild-type Phe 23 peptide. Conversely, the F5-Phe, 1-Nap, and 2-Nap 23 variants self-assembled at enhanced rates, forming fibrils with greater thermodynamic stability than the wild-type peptide. These results indicate that the high amyloidogenicity of aromatic amino acids is a function of hydrophobicity, β-sheet propensity, and planar geometry and not the ability to form stabilizing or directing π-π bonds.     

X-ray crystallography, CD and kinetic studies revealed the essence of the abnormal behaviors of the cytochrome b5 Phe35-->Tyr mutant.

Conserved phenylalanine 35 is one of the hydrophobic patch residues on the surface of cytochrome b5 (cyt b5). This patch is partially exposed on the surface of cyt b5 while its buried face is in direct van der Waals' contact with heme b. Residues Phe35 and Phe/Tyr74 also form an aromatic channel with His39, which is one of the axial ligands of heme b. By site-directed mutagenesis we have produced three mutants of cyt b5: Phe35-->Tyr, Phe35-->Leu, and Phe35-->His. We found that of these three mutants, the Phe35-->Tyr mutant displays abnormal properties. The redox potential of the Phe35-->Tyr mutant is 66 mV more negative than that of the wild-type cyt b5 and the oxidized Phe35-->Tyr mutant is more stable towards thermal and chemical denaturation than wild-type cyt b5. In this study we studied the most interesting mutant, Phe35-->Tyr, by X-ray crystallography, thermal denaturation, CD and kinetic studies of heme dissociation to explore the origin of its unusual behaviors. Analysis of crystal structure of the Phe35-->Tyr mutant shows that the overall structure of the mutant is basically the same as that of the wild-type protein. However, the introduction of a hydroxyl group in the heme pocket, and the increased van der Waals' and electrostatic interactions between the side chain of Tyr35 and the heme probably result in enhancement of stability of the Phe35-->Tyr mutant. The kinetic difference of the heme trapped by the heme pocket also supports this conclusion. The detailed conformational changes of the proteins in response to heat have been studied by CD for the first time, revealing the existence of the folding intermediate.     

NMR studies of 4-19F-phenylalanine-labeled intestinal fatty acid binding protein: evidence for conformational heterogeneity in the native state

(19)F-Nuclear magnetic resonance (NMR) studies have been carried out after incorporation of 4-(19)F-phenylalanine into the intestinal fatty acid binding protein (IFABP), a protein composed of two beta-sheets containing a large hydrophobic cavity into which ligands bind. NMR spectra have been obtained with both the ligand-free and ligand-bound (oleate) forms. There are 29 residues involved in van der Waals or hydrophobic interactions or both to form a U-shaped ligand binding pocket (Sacchettni J. C., Scapin G., Gopaul D., and Gordon J. I. (1992) J. Biol. Chem. 267, 23534-23545). The protein contains eight phenylalanines, and all are included in those residues that line the pocket. Peak assignments were made using site-specific incorporation of 4-(19)F-phenylalanine. Fluorine is a highly sensitive probe to monitor the conformation and dynamics of the side chains in native state. We find that chemical exchange in the binding pocket exists in the native apo- and holo-state. Of the eight phenylalanine residues, Phe2, Phe47, Phe62, Phe68, and Phe93 are arranged on one side of the binding pocket, and all exist in two conformations with Phe2, Phe47, and Phe62 showing exchange cross-peaks with minor conformation in (19)F-(19)F nuclear Overhauser effect (NOESY) spectra. The line widths of Phe68 and Phe93 are broader than those of other phenylalanine residues and can be deconvoluted into two peaks. Phe47, Phe62, Phe68, Phe93, and Trp82 have been proposed to be involved in the early stage of collapse (Ropson, I. J., and Frieden, C. (1992) Proc. Natl. Acad. Sci U.S.A. 89, 7222-7226), but a temperature study suggests that Phe47 behaves differently than other residues and may be more involved in a later stage of folding, for example, side chain stabilization. In the holo-form, Phe17 shows an extra exchange cross-peak in addition to those exchange cross-peaks observed in apo-form. Holo-IFABP exhibits broader line width than the apo-form, suggesting more flexibility of the binding cavity upon ligand binding.     

Critical contribution of aromatic rings to specific recognition of polyether rings. The case of ciguatoxin CTX3C-ABC and its specific antibody 1C49

To address how proteins recognize polyether toxin compounds, we focused on the interaction between the ABC ring compound of ciguatoxin 3C and its specific antibody, 1C49. Surface plasmon resonance analyses indicated that Escherichia coli-expressed variable domain fragments (Fv) of 1C49 had the high affinity constants and slow dissociation constants typical of antigen-antibody interactions. Linear van't Hoff analyses suggested that the interaction is enthalpy-driven. We resolved the crystal structure of 1C49 Fv bound to ABC ring compound of ciguatoxin 3C at a resolution of 1.7A. The binding pocket of the antibody had many aromatic rings and bound the antigen by shape complementarity typical of hapten-antibody interactions. Three hydrogen bonds and many van der Waals interactions were present. We mutated several residues of the antibody to Ala, and we used surface plasmon resonance to analyze the interactions between the mutated antibodies and the antigen. This analysis identified Tyr-91 and Trp-96 in the light chain as hot spots for the interaction, and other residues made incremental contributions by conferring enthalpic advantages and reducing the dissociation rate constant. Systematic mutation of Tyr-91 indicated that CH-pi and pi-pi interactions between the aromatic ring at this site and the antigen made substantial contributions to the association, and van der Waals interactions inhibited dissociation, suggesting that aromaticity and bulkiness are critical for the specific recognition of polyether compounds by proteins.     

Evidence for lipid packaging in the crystal structure of the GM2-activator complex with platelet activating factor

GM2-activator protein (GM2-AP) is a lipid transfer protein that has the ability to stimulate the enzymatic processing of gangliosides as well as T-cell activation through lipid presentation. Our previous X-ray crystallographic studies of GM2-AP have revealed a large lipid binding pocket as the central overall feature of the structure with non-protein electron density within this pocket suggesting bound lipid. To extend these studies, we present here the 2A crystal structure of GM2-AP complexed with platelet activating factor (PAF). PAF is a potent phosphoacylglycerol whose toxic patho-physiological effects can be inhibited by GM2-AP. The structure shows an ordered arrangement of two bound lipids and a fatty acid molecule. One PAF molecule binds in an extended conformation within the hydrophobic channel that has an open and closed conformation, and was seen to contain bound phospholipid in the low pH apo structure. The second molecule is submerged inside the pocket in a U-shaped conformation with its head group near the single polar residue S141. It was refined as lyso-PAF as it lacks electron density for the sn-2 acetate group. The alkyl chains of PAF interact through van der Waals' contacts, while the head groups bind in different environments with their phosphocholine moieties in contact with aromatic rings (Y137, F80). The structure has revealed further insights into the lipid binding properties of GM2-AP, suggesting an unexpected unique mode of lipid packaging that may explain the efficiency of GM2-AP in inhibiting the detrimental biological effects of PAF.     

Phenol red interacts with the protofibril-like oligomers of an amyloidogenic hexapeptide NFGAIL through both hydrophobic and aromatic contacts

Amyloid-associated diseases affect millions of people worldwide. Phenol red exhibits modest inhibition toward fibril formation of human Islet amyloid polypeptide (hIAPP) and its toxicity, which is associated with type II diabetes mellitus. However, the molecular level mechanisms of interactions remain elusive. The binding of phenol red molecules to the protofibrils of an amyloidogenic fragment (NFGAIL) of hIAPP has been investigated by molecular dynamics simulations with explicit solvent. The phenol red molecules were observed to bind primarily along either beta-sheet stacking or beta-strand directions. Through its three aromatic rings, the phenol red molecule preferentially interacted with the hydrophobic side chains of Phe, Leu, and Ile; and the polar sulfone and hydroxyl groups were mainly exposed in solvent. Thus, phenol red improves the solubility of the early protofibrils and represses further growth. Interestingly, there was no obvious preference toward the aromatic Phe residue in comparison to the hydrophobic Leu or Ile residues. The lack of binding along the hydrogen bond direction indicates that phenol red does not directly block the beta-sheet extension. Further free energy analysis suggested that a phenol red analog may potentially improve the binding affinity.     

Contributions of the interdomain loop, amino terminus, and subunit interface to the ligand-facilitated dimerization of neurophysin: crystal structures and mutation studies of bovine neurophysin-I

Current evidence indicates that the ligand-facilitated dimerization of neurophysin is mediated in part by dimerization-induced changes at the hormone binding site of the unliganded state that increase ligand affinity. To elucidate other contributory factors, we investigated the potential role of neurophysin's short interdomain loop (residues 55-59), particularly the effects of loop residue mutation and of deleting amino-terminal residues 1-6, which interact with the loop and adjacent residues 53-54. The neurophysin studied was bovine neurophysin-I, necessitating determination of the crystal structures of des 1-6 bovine neurophysin-I in unliganded and liganded dimeric states, as well as the structure of its liganded Q58V mutant, in which peptide was bound with unexpectedly increased affinity. Increases in dimerization constant associated with selected loop residue mutations and with deletion of residues 1-6, together with structural data, provided evidence that dimerization of unliganded neurophysin-I is constrained by hydrogen bonding of the side chains of Gln58, Ser56, and Gln55 and by amino terminus interactions, loss or alteration of these hydrogen bonds, and probable loss of amino terminus interactions, contributing to the increased dimerization of the liganded state. An additional intersubunit hydrogen bond from residue 81, present only in the liganded state, was demonstrated as the largest single effect of ligand binding directly on the subunit interface. Comparison of bovine neurophysins I and II indicates broadly similar mechanisms for both, with the exception in neurophysin II of the absence of Gln55 side chain hydrogen bonds in the unliganded state and a more firmly established loss of amino terminus interactions in the liganded state. Evidence is presented that loop status modulates dimerization via long-range effects on neurophysin conformation involving neighboring Phe22 as a key intermediary.