Specificity of DNA-basic polypeptide interactions. III. Sequential and random copolymers of lysine and aromatic amino acids.

Sequential and random lysine copolymers containing various amounts of different aromatic amino acids were synthesized. The sequential copolypeptides exhibited strong dependence of yield and degree of polymerization on the amino acid sequence of the repeating unit. To elucidate the specific contributions of aromatic side chains to the interaction of these copolymers with DNA, direct-mixed complexes were studied by thermal denaturation and CD. The melting behaviour of peptide-bound DNA was found to be strongly affected by amino acid composition and sequence. The contribution of the different aromatic amino acids to thermal stability decreased in the order: polylysine > [Lys, Tyr]_n > [Lys,Phe]_n > [Lys,(OMe)Tyr]_n. The CD spectrum of DNA was altered by random copolymers, whereas sequential copolymers exhibited no changes. The influence of the random copolymers on the CD spectrum of DNA decreased in the series: polylysine > [Lys,Phe]_n > [Lys,(OMe)Tyr]_n > [Lys,Tyr]_n. The contribution of the different aromatic amino acids to thermal stability is interpreted as stacking tendencies toward denatured and, in the case of Tyr, H-bond formation with native DNA. The differences found for the random and the sequential polypeptides can best be explained by assuming a cooperative action of rather small peptide segments.     

Reinvestigation of the aromatic side-chains in the basic pancreatic trypsin inhibitor by heteronuclear two-dimensional nuclear magnetic resonance

The 1H nuclear magnetic resonance (n.m.r.) assignments for the aromatic spin systems of the four tyrosines and four phenylalanines in the basic pancreatic trypsin inhibitor (BPTI) were reinvestigated using novel 13C-1H heteronuclear two-dimensional experiments. Resonance lines which are degenerate in homonuclear 1H n.m.r. spectra could thus be resolved. Based on this new evidence the previous assignments for Phe22 and Phe33 had to be corrected. This affects the earlier conclusions on aromatic ring flips in BPTI in that Phe22 is rotating rapidly on the n.m.r. time scale at 36 degrees C, rather than being immobilized up to 80 degrees C.     

Dynamics of the aromatic amino acid residues in the globular conformation of the basic pancreatic trypsin inhibitor (BPTI). I. 1H NMR studies.

The basic pancreatic trypsin inhibitor (BPTI) was investigated by high resolution 1H NMR techniques at 360 MHz. Observation of the amide proton resonances of the polypeptide backbone showed that the globular conformation of BPTI determined by X-ray studies in single crystals is maintained in aqueous solution over the temperature range from 4 degrees to 87 degrees. NMR studies over this temperature range of the aromatic amino acid residues of BPTI. i.e. 4 tyrosines and 4 phenylalanines, led to complete assignments of all the aromatic spin systems in the protein. From this, information was obtained on the rotational motions about the C beta--Cv bond axis of the aromatic rings in the globular form of PBTI. At 25 degrees, two tyrosine rings and one phenylalanine ring are rotating rapidly on the NMR time scale. For the other rings the transitions from slow to rapid rotational motions were investigated at variable temperatures and energy barriers for these intramolecular rate processes determined. The studies of the tyrosine resonances had been described in detail in a previous publication. The present paper describes the identification of the phenylalanine resonances and comments on some technical aspects which might be of quite general interest for the analysis of highly resolved 1H NMR spectra of proteins. Data for the tyrosines and the phenylalanines are compiled in three tables, i.e. the pK alpha-values for the tyrosines, the NMR parameters for all eight aromatics, and the parameters delta G not equal to, and, where available, delta H not equal to and delta S not equal to for the rotational motions of the rings.     

Peptide hairpins with strand segments containing alpha- and beta-amino acid residues: cross-strand aromatic interactions of facing Phe residues

The incporation of beta-amino acid residues into the strand segments of designed beta-hairpin leads to the formation of polar sheets, since in the case of beta-peptide strands, all adjacent carbonyl groups point in one direction and the amide groups orient in the opposite direction. The conformational analysis of two designed peptide hairpins composed of alpha/beta-hybrid segments are described: Boc-Leu-betaPhe-Val-(D)-Pro-Gly-Leu-betaPhe-Val-OMe (1) and Boc-betaLeu-Phe-betaVal-D-Pro-Gly-betaLeu-Phe-betaVal-OMe (2). A 500-MHz 1H-NMR (nuclear magnetic resonance) analysis in methanol supports a significant population of hairpin conformations in both peptides. Diagnostic nuclear Overhauser effects (NOEs) are observed in both cases. X-ray diffraction studies on single crystals of peptide 1 reveal a beta-hairpin conformation in both the molecules, which constitute the crystallographic asymmetric unit. Three cross-strand hydrogen bonds and a nucleating type II' beta-turn at the D-Pro-Gly segment are observed in the two independent molecules. In peptide 1, the betaPhe residues at positions 2 and 7 occur at the nonhydrogen-bonding position, with the benzyl side chains pointing on opposite faces of the beta-sheet. The observed aromatic centroid-to-centroid distances are 8.92 A (molecule A) and 8.94 A (molecule B). In peptide 2, the aromatic rings must occupy facing positions in antiparallel strands, in the NMR-derived structure.Peptide 1 yields a normal "hairpin-like" CD spectrum in methanol with a minimum at 224 nm. The CD spectrum of peptide 2 reveals a negative band at 234 nm and a positive band at 221 nm, suggestive of an exciton split doublet. Modeling of the facing Phe side chains at the hydrogen-bonding position of a canonical beta-hairpin suggests that interring separation is approximately 4.78 A for the gauche+ gauche- (g+ g-) rotamer. A previously reported peptide beta-hairpin composed of only alpha-amino acids, Boc-Leu-Phe-Val-D-Pro-Gly-Leu-Phe-Val-OMe also exhibited an anomalous far-UV (ultraviolet) CD (circular dichroism) spectrum, which was interpreted in terms of interactions between facing aromatic chromophores, Phe 2 and Phe 7 (C. Zhao, P. L. Polavarapu, C. Das, and P. Balaram, Journal of the American Chemical Society, 2000, Vol 122, pp. 8228-8231).     

Aromatic side-chain contributions to the far ultraviolet circular dichroism of peptides and proteins

The rotational strength of the L_a transition in phenylalanine and tyrosine side chains has been calculated for dipeptides with various backbone and side-chain conformations. Similar calculations have also been performed for tripeptides in the β-turn conformation with aromatic residues at the corners of the turn. The interaction of the aromatic ring with neighboring peptides generates rotational strengths in the L_a transition of the order of 0.1 Debye-Bohr magneton. When the preferred backbone and side-chain conformations are considered, it is found that the most probable conformations have positive L_a bonds. This result accounts for the observation that the N-acyl amino acid amides of L -Tyr and L Phe have positive L_a bands. It also suggests that, although other interactions may affect the numerical value and even the sign, there will be a significant positive contribution to the rotational strength of aromatic residues in globular proteins from nearest-neighbor interactions. Calculations on proteins of known conformation at the nearest-neighbor level confirm the tendency toward positive L_a contributions for Phe and Tyr residues. This contribution can be of the order of 10% of the observed CD even in proteins with rather strong amide contributions. In some proteins, such as the gene 5 protein from bacteriophage fd and many snake-venom toxins, side-chain contributions from Tyr and Trp residues manifest themselves as positive CD bands in the 225–250-nm region. The magnitude of the nearest-neighbor contributions and the trend toward positive contributions are consistent with the observation of such CD bands in globular proteins. No special stacking interaction among aromatic side chains needs to be invoked.     

Characterization of a molten globule state of bovine carbonic anhydrase III: loss of asymmetrical environment of the aromatic residues has a profound effect on both the near- and far-UV CD spectrum

Bovine muscle carbonic anhydrase (isoenzyme III; BCAIII) exhibited a three-state unfolding process at equilibrium upon denaturation in guanidine hydrochloride (GuHCl). The stable folding intermediate appeared to be of molten globule type. The stability towards GuHCl in terms of mid-point concentrations of denaturation were very similar for BCAIII and human CAII (HCAII). It was further demonstrated that the aromatic amino acid residues contributed significantly to the circular dichroism (CD) spectrum in the far-UV wavelength region during the native-->molten globule state transition. Thus, the ellipiticity change at 218 nm was shown to monitor the loss of tertiary interactions of aromatic side chains at the first unfolding transition as well as the rupture of secondary structure at the second unfolding transition. Similar aromatic contributions to the far-UV CD spectrum, but with varying magnitudes, were also noted for BCAII and HCAII, further emphasizing that interference of aromatic residues should not be neglected at wavelengths that normally are assigned to secondary structural changes.     

Accuracy of protein secondary structure determination from circular dichroism spectra based on immunoglobulin examples

Strong contribution of the aromatic amino acid side chain chromophores to the far-UV circular dichroism (CD) spectra substantially distorts a relatively weak CD signal originating from beta sheet, the main type of immunoglobulin secondary structure. In this study we compared the secondary structure calculated from the far-UV CD spectra with the X-ray data for three antibody Fab fragments. Calculations were performed with three different algorithms, using two sets of reference proteins. Low standard deviations between all six estimates indicate stable mathematical solutions. Despite pronounced differences in the shape and amplitude of the CD spectra, we found a strong correlation between CD and X-ray data in the secondary structure for every protein studied. The number and average length of the secondary structure elements estimated from the CD spectra closely resemble those of the X-ray data. Agreement between spectroscopic and crystallographic results demonstrates that modern methods of secondary structure calculation are resilient to distortions of the far-UV CD spectra of immunoglobulins caused by aromatic side chain chromophores.     

Interaction of angiotensin II with the C-terminal 300-320 fragment of the rat angiotensin II receptor AT1a monitored by NMR

Interaction between angiotensin II (Ang II) and the fragment peptide 300-320 (fCT300-320) of the rat angiotensin II receptor AT1a was demonstrated by relaxation measurements, NOE effects, chemical shift variations, and CD measurements. The correlation times modulating dipolar interactions for the bound and free forms of Ang II were estimated by the ratio of the nonselective and single-selective longitudinal relaxation rates. The intermolecular NOEs observed in NOESY spectra between HN protons of 9Lys(fCT) and 6His(ang), 10Phe(fCT) and 8Phe(ang), HN proton of 3Tyr(fCT) and Halpha of 4Tyr(ang), 5Phe(fCT)Hdelta and Halpha of 4Tyr(ang) indicated that Ang II aromatic residues are directly involved in the interaction, as also verified by relaxation data. Some fCT300-320 backbone features were inferred by the CSI method and CD experiments revealing that the presence of Ang II enhances the existential probability of helical conformations in the fCT fragment. Restrained molecular dynamics using the simulated annealing protocol was performed with intermolecular NOEs as constraints, imposing an alpha-helix backbone structure to fCT300-320 fragment. In the built model, one strongly preferred interaction was found that allows intermolecular stacking between aromatic rings and forces the peptide to wrap around the 6Leu side chain of the receptor fragment.     

Role of interactions involving C-terminal nonpolar residues of hirudin in the formation of the thrombin-hirudin complex

The role of interactions involving C-terminal nonpolar residues of hirudin in the formation of the thrombin-hirudin complex has been investigated by site-directed mutagenesis. The residues Phe56, Pro60, and Tyr63 of hirudin were replaced by a number of different amino acids, and the kinetics of the inhibition of thrombin by the mutant proteins were determined. Phe56 could be replaced by aromatic amino acids without significant loss in binding energy. While substitution of Phe56 by alanine decreased the binding energy (delta G degrees b by only 1.9 kJ mol-1, replacement of this residue by amino acids with branched side chains caused larger decreases in delta G degrees b. For example, the mutant Phe56----Val displayed a decrease in delta G degrees b of 10.5 kJ mol-1. Substitution of Pro60 by alanine or glycine resulted in a decrease in delta G degrees b of about 6 kJ mol-1. Tyr63 could be replaced by phenylalanine without any loss in binding energy, and replacement of this residue by alanine caused a decrease of 2.2 kJ mol-1 in delta G degrees b. Substitution of Tyr63 by residues with branched side chains resulted in smaller decreases in delta G degrees b than those seen with the corresponding substitutions of Phe56; for example, the mutant Tyr63----Val showed a decrease in binding energy of 5.1 kJ mol-1. The effects of the mutations are discussed in terms of the crystal structure of the thrombin-hirudin complex.     

Identification of a residue critical for maintaining the functional conformation of BPTI

The effects of amino acid replacements on the backbone dynamics of bovine pancreatic trypsin inhibitor (BPTI) were examined using 15N NMR relaxation experiments. Previous studies have shown that backbone amide groups within the trypsin-binding region of the wild-type protein undergo conformational exchange processes on the micros time scale, and that replacement of Tyr35 with Gly greatly increases the number of backbone atoms involved in such motions. In order to determine whether these mutational effects are specific to the replacement of this residue with Gly, six additional replacements were examined in the present study. In two of these, Tyr35 was replaced with either Ala or Leu, and the other four were single replacements of Tyr23, Phe33, Asn43 or Asn44, all of which are highly buried in the native structure and conserved in homologous proteins. The Y35A and Y35L mutants displayed dynamic properties very similar to those of the Y35G mutant, with the backbone segments including residues 10-19 and 32-44 undergoing motions revealed by enhanced 15N transverse relaxation rates. On the other hand, the Y23L, N43G and N44A substitutions caused almost no detectable changes in backbone dynamics, on either the ns-ps or ms-micros time scales, even though each of these replacements significantly destabilizes the native conformation. Replacement of Phe33 with Leu caused intermediate effects, with several residues that have previously been implicated in motions in the wild-type protein displaying enhanced transverse relaxation rates. These results demonstrate that destabilizing amino acid replacements can be accommodated in a native protein with dramatically different effects on conformational dynamics and that Tyr35 plays a particularly important role in defining the conformation of the trypsin-binding site of BPTI.     

Solution structure of a Kunitz-type chymotrypsin inhibitor isolated from the elapid snake Bungarus fasciatus

Bungarus fasciatus fraction IX (BF9), a chymotrypsin inhibitor, consists of 65 amino acid residues with three disulfide bridges. It was isolated from the snake venom of B. fasciatus by ion-exchange chromatography and belongs to the bovine pancreatic trypsin inhibitor (BPTI)-like superfamily. It showed a dissociation constant of 5.8 x 10(-8) m with alpha-chymotrypsin as measured by a BIAcore binding assay system. The isothermal titration calorimetry revealed a 1:1 binding stoichiometry between this inhibitor and chymotrypsin and apparently no binding with trypsin. We further used CD and NMR to determine the solution structure of this venom-derived chymotrypsin inhibitor. The three-dimensional NMR solution structures of BF9 were determined on the basis of 582 restraints by simulated annealing and energy minimization calculations. The final set of 10 NMR structures was well defined, with average root mean square deviations of 0.47 A for the backbone atoms in the secondary structure regions and 0.86 A for residues The side chains of Phe(23), Tyr(24), Tyr(25), Phe(35), and Phe(47) exhibited many long-range nuclear Overhauser effects and were the principal components of the hydrophobic core in BF9. To gain insight into the structure-function relationships among proteins in the BPTI-like superfamily, we compared the three-dimensional structure of BF9 with three BPTI-like proteins that possess distinct biological functions. These proteins possessed similar secondary structure elements, but the loop regions and beta-turn were different from one another. Based on residues at the functional site of each protein, we suggest that the flexibility, rigidity, and variations of the amino acid residues in both the loop and beta-turn regions are related to their biological functions.     

Application of SAIL phenylalanine and tyrosine with alternative isotope-labeling patterns for protein structure determination

The extensive collection of NOE constraint data involving the aromatic ring signals is essential for accurate protein structure determination, although it is often hampered in practice by the pervasive signal overlapping and tight spin couplings for aromatic rings. We have prepared various types of stereo-array isotope labeled phenylalanines (ε- and ζ-SAIL Phe) and tyrosine (ε-SAIL Tyr) to overcome these problems (Torizawa et al. 2005), and proven that these SAIL amino acids provide dramatic spectral simplification and sensitivity enhancement for the aromatic ring NMR signals. In addition to these SAIL aromatic amino acids, we recently synthesized δ-SAIL Phe and δ-SAIL Tyr, which allow us to observe and assign δ-¹³C/¹H signals very efficiently. Each of the various types of SAIL Phe and SAIL Tyr yields well-resolved resonances for the δ-, ε- or ζ-¹³C/¹H signals, respectively, which can readily be assigned by simple and robust pulse sequences. Since the δ-, ε-, and ζ-proton signals of Phe/Tyr residues give rise to complementary NOE constraints, the concomitant use of various types of SAIL-Phe and SAIL-Tyr would generate more accurate protein structures, as compared to those obtained by using conventional uniformly ¹³C, ¹⁵N-double labeled proteins. We illustrated this with the case of an 18.2 kDa protein, Escherichia coli peptidyl-prolyl cis-trans isomerase b (EPPIb), and concluded that the combined use of ζ-SAIL Phe and ε-SAIL Tyr would be practically the best choice for protein structural determinations.     

Probing aromatic, hydrophobic, and steric effects on the self-assembly of an amyloid-β fragment peptide

Aromatic amino acids have been shown to promote self-assembly of amyloid peptides, although the basis for this amyloid-inducing behavior is not understood. We adopted the amyloid-β 16-22 peptide (Aβ(16-22), Ac-KLVFFAE-NH(2)) as a model to study the role of aromatic amino acids in peptide self-assembly. Aβ(16-22) contains two consecutive Phe residues (19 and 20) in which Phe 19 side chains form interstrand contacts in fibrils while Phe 20 side chains interact with the side chain of Va l18. The kinetic and thermodynamic effect of varying the hydrophobicity and aromaticity at positions 19 and 20 by mutation with Ala, Tyr, cyclohexylalanine (Cha), and pentafluorophenylalanine (F(5)-Phe) (order of hydrophobicity is Ala < Tyr < Phe < F(5)-Phe < Cha) was characterized. Ala and Tyr position 19 variants failed to undergo fibril formation at the peptide concentrations studied, but Cha and F(5)-Phe variants self-assembled at dramatically enhanced rates relative to wild-type. Cha mutation was thermodynamically stabilizing at position 20 (ΔΔG = -0.2 kcal mol(-1) relative to wild-type) and destabilizing at position 19 (ΔΔG = +0.2 kcal mol(-1)). Conversely, F(5)-Phe mutations were strongly stabilizing at both positions (ΔΔG = -1.3 kcal mol(-1) at 19, ΔΔG = -0.9 kcal mol(-1) at 20). The double Cha and F(5)-Phe mutants showed that the thermodynamic effects were additive (ΔΔG = 0 kcal mol(-1) for Cha 19,20 and -2.1 kcal mol(-1) for F(5)-Phe 19,20). These results indicate that sequence hydrophobicity alone does not dictate amyloid potential, but that aromatic, hydrophobic, and steric considerations collectively influence fibril formation.     

Structural basis for the unusual specificity of Escherichia coli aminopeptidase N

Aminopeptidase N from Escherichia coli is a M1 class aminopeptidase with the active-site region related to that of thermolysin. The enzyme has unusual specificity, cleaving adjacent to the large, nonpolar amino acids Phe and Tyr but also cleaving next to the polar residues Lys and Arg. To try to understand the structural basis for this pattern of hydrolysis, the structure of the enzyme was determined in complex with the amino acids L-arginine, L-lysine, L-phenylalanine, L-tryptophan, and L-tyrosine. These amino acids all bind with their backbone atoms close to the active-site zinc ion and their side chain occupying the S1 subsite. This subsite is in the form of a cylinder, about 10 A in cross-section and 12 A in length. The bottom of the cylinder includes the zinc ion and a number of polar side chains that make multiple hydrogen-bonding and other interactions with the alpha-amino group and the alpha-carboxylate of the bound amino acid. The walls of the S1 cylinder are hydrophobic and accommodate the nonpolar or largely nonpolar side chains of Phe and Tyr. The top of the cylinder is polar in character and includes bound water molecules. The epsilon-amino group of the bound lysine side chain and the guanidinium group of arginine both make multiple hydrogen bonds to this part of the S1 site. At the same time, the hydrocarbon part of the lysine and arginine side chains is accommodated within the nonpolar walls of the S1 cylinder. This combination of hydrophobic and hydrophilic binding surfaces explains the ability of ePepN to cleave Lys, Arg, Phe, and Tyr. Another favored substrate has Ala at the P1 position. The short, nonpolar side chain of this residue can clearly be bound within the hydrophobic part of the S1 cylinder, but the reason for its facile hydrolysis remains uncertain.     

Several properties of complexes of polynucleotides with aromatic amino acid amides of oligonucleotides

Residues of Phe, Tyr and Trp in the complexes of their oligonucleotide amidates and polynucleotides of A-U of G-C nucleotide composition are most likely localized in the minor groove of the Watson--Crick part of the triple helix where they interact with bases but do not intercalate into the helix. Formation of the complexes is accompanied with a change in the relative localization of amino acids and bases. The major geometrical parameters of the triple helices of the complexes are not changed by the residues of aromatic amino acids (according to CD data). A slight violation of stacking interactions between bases is observed along with an increase of the cooperativity of melting of the complexes of A-U composition (according to UV absorption data). The effect of the residues of aromatic amino acids on the stability of triple helices is determined by the nucleotide composition of the latter, i.e. complexes of A-U composition are destabilized with the Phe, Tyr and Trp residues, whereas the Trp residue does not affect the stability of the complexes of G-C composition. The hydrophobic character of aromatic amino acids and their different affinity for bases of different structure seem to account for this difference in stability. The dependence of the thermal stability of RNH-dp(An).2poly(U)-complexes on the structure of the amide radical (residues of glycin, aromatic amino acids, alkyl- and arylalkyl amines) testifies the ability of the radical to "regulate" the interaction between the oligomer and the complementary polynucleotide. This capacity for "regulation" is not observed in the system of G-C composition.     

Local control of peptide conformation: Stabilization of cis proline peptide bonds by aromatic proline interactions

In the native state of proteins there is a marked tendency for an aromatic amino acid to precede a cis proline. There are also significant differences between the three aromatic amino acids with Tyr exhibiting a noticeably higher propensity than Phe or Trp to precede a cis proline residue. In order to study the role that local interactions play in these conformation preferences, a set of tetrapeptides of the general sequence acetyl-Gly-X-Pro-Gly-carboxamide (GXPG), where X = Tyr, Phe, Trp, Ala, or cyclohexyl alanine, were synthesized and studied by nmr. Analysis of the nmr data shows that none of the peptides adopt a specific backbone structure. Ring current shifts, the equilibrium constants, the Van't Hoff enthalpy, and the measured rate of cis-trans isomerization all indicate that the cis proline conformer is stabilized by favorable interactions between the aromatic ring and the proline residue. Analysis of the side chain conformation of the aromatic residue and analysis of the chemical shifts of the pyrrolidine ring protons shows that the aromatic side chain adopts a preferred conformation in the cis form. The distribution of rotamers and the effect of an aromatic residue on the cis-trans equilibrium indicate that the preferred conformation is populated to approximately 62% for the Phe containing peptide, 67% for the Tyr containing peptide, and between 75 and 80% for the Trp containing peptide. The interaction is unaffected by the addition of 8M urea. These local interactions favor an aromatic residue immediately preceding a cis proline, but they cannot explain the relative propensities for Phe-Pro, Tyr-Pro, and Trp-Pro cis peptide bonds observed in the native state of proteins. In the model peptides the percentage of the cis proline conformer is 21% GYPG while it is 17% for GFPG. This difference is considerably smaller than the almost three to one preponderance observed for cis Tyr-Pro peptide bonds vs cis Phe-Pro peptide bonds in the protein database. © 1998 John Wiley & Sons, Inc. Biopoly 45: 381–394, 1998     

Topography of the prostaglandin endoperoxide H2 synthase-2 in membranes

The topology of association of the monotopic protein cyclooxygenase-2 (COX-2) with membranes has been examined using EPR spectroscopy of spin-labeled recombinant human COX-2. Twenty-four mutants, each containing a single free cysteine substituted for an amino acid in the COX-2 membrane binding domain were expressed using the baculovirus system and purified, then conjugated with a nitroxide spin label and reconstituted into liposomes. Determining the relative accessibility of the nitroxide-tagged amino acid side chains for the solubilized COX-2 mutants, or COX-2 reconstituted into liposomes to nonpolar (oxygen) and polar (NiEDDA or CrOx) paramagnetic reagents allowed us to map the topology of COX-2 interaction with the lipid bilayer. When spin-labeled COX-2 was reconstituted into liposomes, EPR power saturation curves showed that side chains for all but two of the 24 mutants tested had limited accessibility to both polar and nonpolar paramagnetic relaxation agents, indicating that COX-2 associates primarily with the interfacial membrane region near the glycerol backbone and phospholipid head groups. Two amino acids, Phe(66) and Leu(67), were readily accessible to the non-polar relaxation agent oxygen, and thus likely inserted into the hydrophobic core of the lipid bilayer. However these residues are co-linear with amino acids in the interfacial region, so their extension into the hydrophobic core must be relatively shallow. EPR and structural data suggest that membrane interaction of COX-2 is also aided by partitioning of 4 aromatic amino acids, Phe(59), Phe(66), Tyr(76), and Phe(84) to the interfacial region, and by the electrostatic interactions of two basic amino acids, Arg(62) and Lys(64), with the phospholipid head groups.     

Tuning β-sheet peptide self-assembly and hydrogelation behavior by modification of sequence hydrophobicity and aromaticity

Peptide self-assembly leading to cross-β amyloid structures is a widely studied phenomenon because of its role in amyloid pathology and the exploitation of amyloid as a functional biomaterial. The self-assembly process is governed by hydrogen bonding, hydrophobic, aromatic π-π, and electrostatic Coulombic interactions. A role for aromatic π-π interactions in peptide self-assembly leading to amyloid has been proposed, but the relative contributions of π-π versus general hydrophobic interactions in these processes are poorly understood. The Ac-(XKXK)(2)-NH(2) peptide was used to study the contributions of aromatic and hydrophobic interactions to peptide self-assembly. Position X was globally replaced by valine (Val), isoleucine (Ile), phenylalanine (Phe), pentafluorophenylalanine (F(5)-Phe), and cyclohexylalanine (Cha). At low pH, these peptides remain monomeric because of repulsion of charged lysine (Lys) residues. Increasing the solvent ionic strength to shield repulsive charge-charge interactions between protonated Lys residues facilitated cross-β fibril formation. It was generally found that as peptide hydrophobicity increased, the required ionic strength to induce self-assembly decreased. At [NaCl] ranging from 0 to 1000 mM, the Val sequence failed to assemble. Assembly of the Phe sequence commenced at 700 mM NaCl and at 300 mM NaCl for the less hydrophobic Ile variant, even though it displayed a mixture of random coil and β-sheet secondary structures over all NaCl concentrations. β-Sheet formation for F(5)-Phe and Cha sequences was observed at only 20 and 60 mM NaCl, respectively. Whereas self-assembly propensity generally correlated to peptide hydrophobicity and not aromatic character the presence of aromatic amino acids imparted unique properties to fibrils derived from these peptides. Nonaromatic peptides formed fibrils of 3-15 nm in diameter, whereas aromatic peptides formed nanotape or nanoribbon architectures of 3-7 nm widths. In addition, all peptides formed fibrillar hydrogels at sufficient peptide concentrations, but nonaromatic peptides formed weak gels, whereas aromatic peptides formed rigid gels. These findings clarify the influence of aromatic amino acids on peptide self-assembly processes and illuminate design principles for the inclusion of aromatic amino acids in amyloid-derived biomaterials.     

Molecular dynamics study of a complex between the human histocompatibility antigen HLA-A2 and the IMP58-66 nonapeptide from influenza virus matrix protein.

The structure of the influenza-virus-matrix-protein (IMP) 58-66 nonapeptide, bound to the major-histocompatibility-complex-encoded human leukocyte antigen (HLA) A2 protein was studied by molecular dynamics simulation. Starting from the extra electron density map of peptides co-crystallized with HLA-A2, the nonapeptide IMP58-66 was docked residue by residue in the protein binding cleft. The complex was simulated for 100 ps in a shell of 1372 water molecules. The averaged simulated HLA-A2 conformation was found to be similar to the crystal structure (0.182 nm RMS deviation, for the backbone atoms of the alpha 1-alpha 2 domain). Nine out of the 14 hydrogen bonds observed in the antigen-binding site were reproduced in the simulation. The IMP58-66 peptide exhibits an extended conformation with kinks at positions 3 and 5. The side chains of residues 2, 3 and 9 develop van der Waals' interactions with hydrophobic pockets of HLA-A2, corresponding to polymorphic residues of the major-histocompatibility-complex-encoded proteins. Both the N-terminus and C-terminus of the nonapeptide were anchored in the antigen-binding groove by hydrogen bonds with conserved amino acids. The N-terminus was more flexible and contacts four HLA-A2 conserved tyrosines (Tyr7, Tyr59, Tyr159 and Tyr171) and Glu63 by direct or water-relayed hydrogen bonds. Water intercalation occurred only around the N-terminus of the peptide, the C-terminal carboxylate forming strong hydrogen bonds with polar residues (Tyr84 and Thr143) and a salt bridge with Lys146 all over the molecular dynamics simulation. This model is fully compatible with the recently published crystal structure of the HLA-B27 protein, complexed by a mixture of self nonapeptides.