Role of hydrogen bond networks and dynamics in positive and negative cooperative stabilization of a protein

Cooperativity mediated through hydrogen bond networks in yeast iso-1-cytochrome c was studied using a thermodynamic triple mutant cycle. Three known stabilizing mutations, Asn 26 to His, Asn 52 to Ile, and Tyr 67 to Phe, were used to construct the triple mutant cycle. The side chain of His 26, a wild-type residue, forms two hydrogen bonds that bridge two substructures of the wild-type protein, and Tyr 67 and Asn 52 are part of an extensive buried hydrogen bond network. The stabilities of all variants in the triple mutant cycle were determined by guanidine hydrochloride denaturation methods and used to determine the pairwise, Delta(2)G(int), and triple interaction energies. His 26 and Ile 52 interact cooperatively (Delta(2)G(int) is 1-2 kcal/mol), whereas the two other pairs of mutations interact anticooperatively (Delta(2)G(int) is -0.5 to -1.5 kcal/mol). Previously reported structural data for iso-1-cytochrome c variants containing these mutations show that changes in the strength of the His 26 to Glu 44 hydrogen bond, apparently caused by changes in main chain dynamics, provide a mechanism for the long distance (His 26 to Phe 67 and His 26 to Ile 52) propagation of pairwise interaction energies. Opposing changes in the strength of the His 26 to Glu 44 hydrogen bond caused by the N52I and Y67F mutations generate a negative triple interaction energy (-0.9 +/-0.7 kcal/mol) that combined with cancellation of cooperative and anticooperative pairwise interactions produce apparent additivity for the stabilizing effects of the single mutations in the triple mutant variant.     

Investigation of the PDZ domain ligand binding site using chemically modified peptides

Several chemically modified analogues to a tightly binding ligand for the second PDZ domain of MAGI-3 were synthesized and evaluated for their ability to compete with native peptide ligands. N-methyl scanning of the ligand backbone amides revealed the energetically important hydrogen bonds between the ligand backbone and the PDZ domain. Analogues to the ligand's conserved threonine/serine(-2) residue, involved in a side chain to side chain hydrogen bond with a conserved histidine in the PDZ domain, revealed that the interaction is highly sensitive to the steric structure around the hydroxyl group of this residue. Analogues of the ligand carboxy terminus revealed that the full hydrogen bond network of the GLGF loop is important in ligand binding.     

Analysis of the structure and stability of omega loop A replacements in yeast iso-1-cytochrome c.

Omega (omega)-loop A, residues 18-32 in wild-type yeast iso-1-cytochrome c, has been deleted and replaced with loop sequences from three other cytochromes c and one from esterase. Yeast expressing a partial loop deletion do not contain perceptible amounts of holoprotein as measured by low-temperature spectroscopy and cannot grow on nonfermentable media. Strains expressing loop replacement mutations accumulate holoprotein in vivo, but the protein function varies depending on the sequence and length of the replacement loop; in vivo expression levels do not correlate with their thermal denaturation temperatures. In vitro spectroscopic studies of the loop replacement proteins indicate that all fold into a native-like cytochrome c conformation, but are less stable than the wild-type protein. Decreases in thermal stability are caused by perturbation of loop C backbone in one case and a slight reorganization of the protein hydrophobic core in another case, rather than rearrangement of the loop A backbone. A single-site mutation in one of the replacement mutants designed to relieve inefficient hydrophobic core packing caused by the new loop recovers some, but not all, of the lost stability.     

The structure and function of omega loop A replacements in cytochrome c.

The structural and functional consequences of replacing omega-loop A (residues 18-32) in yeast iso-1-cytochrome c with the corresponding loop of Rhodospirillum rubrum cytochrome c2 have been examined. The three-dimensional structure of this loop replacement mutant RepA2 cytochrome c, and a second mutant RepA2(Val 20) cytochrome c in which residue 20 was back substituted to valine, were determined using X-ray diffraction techniques. A change in the molecular packing is evident in the RepA2 mutant protein, which has a phenylalanine at position 20, a residue considerably larger than the valine found in wild-type yeast iso-1-cytochrome c. The side chain of Phe 20 is redirected toward the molecular surface, altering the packing of this region of omega-loop A with the hydrophobic core of the protein. In the RepA2(Val 20) structure, omega-loop A contains a valine at position 20, which restores the original wild-type packing arrangement of the hydrophobic core. Also, as a result of omega-loop A replacement, residue 26 is changed from a histidine to asparagine, which results in displacements of the main-chain atoms near residue 44 to which residue 26 is hydrogen bonded. In vivo studies of the growth rate of the mutant strains on nonfermentable media indicate that the RepA2(Val 20) cytochrome c behaves much like the wild-type yeast iso-1 protein, whereas the stability and function of the RepA2 cytochrome c showed a temperature dependence. The midpoint reduction potential measured by cyclic voltammetry of the RepA2 mutant is 271 mV at 25 degrees C. This is 19 mV less than the wild-type and RepA2(Val 20) proteins (290 mV) and may result from disruption of the hydrophobic packing in the heme pocket and increased mobility of omega-loop A in RepA2 cytochrome c. The temperature dependence of the reduction potential is also greatly enhanced in the RepA2 protein.     

Rapid intrachain binding of histidine-26 and histidine-33 to heme in unfolded ferrocytochrome C.

Time-resolved spectroscopic studies of unfolded horse iron(II) cytochrome c have suggested that the imidazole side chains of His26 and His33 bind transiently to the heme iron on microsecond time scales, after photodissociation of a carbon monoxide ligand from the heme. Our studies of four variants of cytochrome c (horse wild type, horse H33N, horse H33N/H26Q, and tuna wild type), unfolded in guanidine hydrochloride at pH 6.5, demonstrate that these side chains are responsible for the observed microsecond spectral changes. As His33 and then His26 are eliminated from the horse wild-type sequence, transient optical absorption spectra show systematic suppression of a rapid (approximately 10-100 micros) Soret absorbance change that follows photolysis of CO. Transient binding of these histidine side chains to the heme therefore generates one of the fast kinetic phases observed in previous photochemically triggered spectroscopic studies of dynamics in unfolded iron(II) cytochrome c. Furthermore, both His33 and His26 appear to contribute to a similar extent in these early kinetics. Thus, the stiffness of the polypeptide chain creates a deviation from Gaussian chain behavior by impeding, although not preventing, the formation of short (<10 peptide bonds) intrachain loops around the heme group.     

Crystal structure of the oxidized cytochrome c(2) from Blastochloris viridis

The crystal structure of the oxidized cytochrome c(2) from Blastochloris (formerly Rhodopseudomonas) viridis was determined at 1.9 A resolution. Structural comparison with the reduced form revealed significant structural changes according to the oxidation state of the heme iron. Slight perturbation of the polypeptide chain backbone was observed, and the secondary structure and the hydrogen patterns between main-chain atoms were retained. The oxidation state-dependent conformational shifts were localized in the vicinity of the methionine ligand side and the propionate group of the heme. The conserved segment of the polypeptide chain in cytochrome c and cytochrome c(2) exhibited some degree of mobility, interacting with the heme iron atom by the hydrogen bond network. These results indicate that the movement of the internal water molecule conserved in various c-type cytochromes drives the adjustments of side-chain atoms of nearby residue, and the segmental temperature factor changes along the polypeptide chain.     

Rupture of the hydrogen bond linking two Omega-loops induces the molten globule state at neutral pH in cytochrome c

His26Tyr and His33Tyr mutants were obtained from the Cys102Thr variant of yeast iso-1-cytochrome c. Spectroscopic studies show that a mutation at position 26 at pH 7.0 enhances flexibility of the peptide, alters the heme pocket region and the axial coordination to heme-iron, and reduces protein stability. The His26Tyr mutant shows properties typical of the molten globule. Further, formation of an axially misligated minor low spin species occurs with partial displacement of Met80, the axial ligand of the heme-iron in the native protein. The pK(a) determined for the alkaline transition of this mutant is 7.48 (+/- 0.05), approximately 0.5 lower than that of the wild-type protein. Hence, the alkaline conformer is populated at pH 7.0, and the sixth ligand of the misligated species is proposed to be a lysine. Furthermore, a reduction in catalytic activity indicates that the functional properties are altered. The results suggest that the structural and functional changes observed in the His26Tyr mutant are because the mutation frees the two Omega-loops that, in the native protein, are linked by the hydrogen bond between His26 and Glu44. Hence, one may infer that the His26-Glu44 hydrogen bond is essential for the rigidity and stability of the native protein. In its absence, the heightened flexibility of the peptide fold results in conversion of the macromolecule to a molten globule state, even at neutral pH. Ligand exchange at the sixth coordination position of the heme-iron(III) observed as the minor species (i.e., the alkaline conformer) is therefore induced by a long-range effect. This result is of interest since mutations reported to date, which stabilize the alkaline conformer, all occur in the loop including Met80. By contrast, only very minor spectroscopic (and, thus, structural) changes are observed for the His33Tyr mutant. This suggests that His33 does not form intramolecular bonds considered important for the protein structure and stability, and is consistent with the high variability of residues at position 33 in cytochromes c.     

An intersubunit hydrogen bond in the nicotinic acetylcholine receptor that contributes to channel gating

The muscle nicotinic acetylcholine receptor is a large, allosteric, ligand-gated ion channel with the subunit composition alpha2betagammadelta. Although much is now known about the structure of the binding site, relatively little is understood about how the binding event is communicated to the channel gate, causing the pore to open. Here we identify a key hydrogen bond near the binding site that is involved in the gating pathway. Using mutant cycle analysis with the novel unnatural residue alpha-hydroxyserine, we find that the backbone N-H of alphaSer-191 in loop C makes a hydrogen bond to an anionic side chain of the complementary subunit upon agonist binding. However, the anionic partner is not the glutamate predicted by the crystal structures of the homologous acetylcholine-binding protein. Instead, the hydrogen-bonding partner is the extensively researched aspartate gammaAsp-174/deltaAsp-180, which had originally been identified as a key binding residue for cationic agonists.     

The structural basis of a conserved P2 threonine in canonical serine proteinase inhibitors

Bowman-Birk inhibitors (BBIs) are a well-studied family of canonical inhibitor proteins of serine proteinases. In nature, the active region of BBIs possesses a highly conserved Thr at the P2 position. The importance of this residue has been reemphasized by synthetic BBI reactive site loop proteinomimetics. In particular, this residue was exclusively identified for active chymotrypsin inhibitors selected from a BBI template-assisted combinatorial peptide library. A further kinetic analysis of 26 P2 variant peptides revealed that Thr provides both optimal binding affinity and optimal resistance against enzymatic turnover by chymotrypsin. Herein, we report the (1)H-NMR spectroscopic study of a 5-membered sub-set of these reactive site loop peptides representing a stepwise elimination of the Thr side-chain functionalities and inversion of its side-chain chirality. The P2 Thr variant adopts a three-dimensional structure that closely mimics the one of the corresponding region of the complete protein. This validates the use of this template for the investigation of structure-function relationships. While the overall backbone geometry is similar in all studied variants, conformational changes induced by the modification of the P2 side chain have now been identified and provide a rational explanation of the kinetically observed functional differences. Eliminating the gamma-methyl group has little structural effect, whereas the elimination of the gamma-oxygen atom or the inversion of the side-chain chirality results in characteristic changes to the intramolecular hydrogen bond network. We conclude that the transannular hydrogen bond between the P2 Thr side-chain hydroxyl and the P5' backbone amide is an important conformational constraint and directs the hydrophobic contact of the P2 Thr side chain with the enzyme surface in a functionally optimal geometry, both in the proteinomimetic and the native protein. In at least four canonical inhibitor protein families similar structural arrangements for a conserved P2 Thr have been observed, which suggests an analogous functional role. Substitutions at P2 of the proteinomimetic also affect the conformational balance between cis and trans isomers at a distant Pro-Pro motif (P3'-P4'). Presented with a mixture of cis/trans isomers chymotrypsin appears to interact preferably with the conformer that retains the cis-P3' Pro-trans-P4' Pro geometry found in the parent BBI protein.     

Catalytic and structural role of a metal-free histidine residue in bovine Cu-Zn superoxide dismutase

Cu-Zn superoxide dismutase (SOD) contains a conserved, metal-free His residue at an opening of the backbone beta-barrel in addition to six Cu- and/or Zn-bound His residues in the active site. We examined the protonation and hydrogen bonding state of the metal-free His residue (His41) in bovine SOD by UV Raman spectroscopy. Analysis of the His Raman intensity at 1406 cm(-1) in a D2O solution has shown that His41 has a pKa of 9.4, consistent with the NMR and X-ray structures at acidic to neutral pH, in which two imidazole nitrogen atoms of cationic His41 are hydrogen bonded to the main chain C=O groups of Thr37 and His118. Upon deprotonation of His41 at pH 9.4, the Thr37-His41-His118 hydrogen bond bridge breaks on the His118 side and SOD loses 70% of its activity. Concomitantly, hydrogen-deuterium exchange is accelerated for amide groups of beta-strands, indicating an increased conformational fluctuation of the beta-barrel. Thr37 and His41 are in direct contact with Leu36, whose hydrophobic side chain closes off the opening of the beta-barrel, while His118 is indirectly connected to Arg141 that assists the docking of superoxide to Cu. These Raman findings strongly suggest that the His41-mediated hydrogen bond bridge plays a crucial role in keeping the protein structure suitable for highly efficient catalytic reactions. The catalytic and structural role of His41 is consistent with the observation that the mutation of His43 in human SOD (equivalent to His41 in bovine SOD) to Arg largely reduces the dismutase activity and the protein structural stability.     

Solution conformation of the His-47 to Ala-47 mutant of Pseudomonas stutzeri ZoBell ferrocytochrome c-551

In the cytochrome c-551 family, the heme 17-propionate caboxylate group is always hydrogen bonded to an invariant Trp-56 and conserved residues (His and Arg mainly, Lys occasionally) at position 47. The mutation of His-47 to Ala-47 for Pseudomas stutzeri ZoBell cytochrome c-551 removes this otherwise invariant hydrogen bond. The solution structure of ferrous H47A has been solved based on NMR-derived constraints. Results indicate that the mutant has very similar main chain folding compared to wild-type. However, less efficient packing of residues in the mutant surrounding the heme propionates leads to more solvent exposure for both propionate groups, which may account for decreased stability of the mutant. The mutant has a reduction potential different from wild-type, and furthermore, the pH dependence of this potential is not the same as for wild-type. The structure of the mutant suggests that these changes are related to the loss of the residue-47 propionate hydrogen bond and the loss of charge on the side chain of residue 47.     

The importance of anchorage in determining a strained protein loop conformation.

We examine the role of the conformational restriction imposed by constrained ends of a protein loop on the determination of a strained loop conformation. The Lys 116-Pro 117 peptide bond of staphylococcal nuclease A exists in equilibrium between the cis and trans isomers. The folded protein favors the strained cis isomer with an occupancy of 90%. This peptide bond is contained in a solvent-exposed, flexible loop of residues 112-117 whose ends are anchored by Val 111 and Asn 118. Asn 118 is constrained by 2 side-chain hydrogen bonds. We investigate the importance of this constraint by replacing Asn 118 with aspartate, alanine, and glycine. We found that removing 1 or more of the hydrogen bonds observed in Asn 118 stabilizes the trans configuration over the cis configuration. By protonating the Asp 118 side chain of N118D through decreased pH, the hydrogen bonding character of Asp 118 approached that of Asn 118 in nuclease A, and the cis configuration was stabilized relative to the trans configuration. These data suggest that the rigid anchoring of the loop end is important in establishing the strained cis conformation. The segment of residues 112-117 in nuclease A provides a promising model system for study of the basic principles that determine polypeptide conformations. Such studies could be useful in the rational design or redesign of protein molecules.     

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.     

The interaction of borate ions with cytochrome c surface sites: a molecular dynamics study.

Ionic interactions of cytochrome c play an important role in the electron transfer process. Molecular dynamics simulations of the binding of borate ion, which serves as a model ion, at three different cytochrome c surface sites are performed. This work is motivated by previous NMR studies of cytochrome c in borate solution, which indicate the existence of two types of binding sites, a slow exchange site and a fast exchange site. These two types of binding behavior were observed in the dynamic simulations, offering a molecular interpretation of "loose" and "tight" binding. At the "loose" binding sites (near Lys25/Lys27 and Lys55/Lys73) the ion forms two to three hydrogen bonds to the nearest lysine residue. This binding is transient on the time scale of the simulation, demonstrating the feasibility of fast exchange. At the "tight" binding site (near Lys13/Lys86), on the other hand, the ion becomes integrated into the protein hydrogen bond network and remains there for the duration of the simulation (exemplifying slow exchange). Binding simulations of the ion at the "tight" site of H26Q mutant cytochrome c also showed integration of the ion into the protein's hydrogen bond network. However, this integration differs in details from the binding of the ion to the native protein, in agreement with previous NMR observations.     

Chemical-scale studies on the role of a conserved aspartate in preorganizing the agonist binding site of the nicotinic acetylcholine receptor

The nicotinic acetylcholine receptor and related Cys-loop receptors are ligand-gated ion channels that mediate fast synaptic transmission throughout the central and peripheral nervous system. A highly conserved aspartate residue (D89) that is near the agonist binding site but does not directly contact the ligand plays a critical part in receptor function. Here we probe the role of D89 using unnatural amino acid mutagenesis coupled with electrophysiology. Homology modeling implicates several hydrogen bonds involving D89. We find that no single hydrogen bond is essential to proper receptor function. Apparently, the side chain of D89 establishes a redundant network of hydrogen bonds; these bonds preorganize the agonist binding site by positioning a critical tryptophan residue that directly contacts the ligand. Earlier studies of the D89N mutant led to the proposal that a negative charge at this position is essential for receptor function. However, we find that receptors with neutral side chains at position 89 can function well, if the side chain is less perturbing than the amide of asparagine (nitro or keto groups allow function) or if a compensating backbone mutation is introduced to relieve unfavorable electrostatics.     

Alanine-scanning of the 50's loop in the Clostridium beijerinckii flavodoxin: evaluation of additivity and the importance of interactions provided by the main chain in the modulation of the oxidation-reduction potentials.

The four-residue reverse turn -Met56-Gly-Asp-Glu59- in the Clostridium beijerinckii flavodoxin provides the majority of the critical interactions with the isoalloxazine ring of the flavin mononucleotide (FMN) cofactor that contribute to the binding and the differential stabilization of its three redox states. Direct side chain contacts include the sulfur-ring interaction of Met56, which primarily influences the oxidized and hydroquinone states, and the hydrogen bond by Glu59 with the N3H, which directly (and indirectly through its "anchoring" function) influences all three states to various extents. Involving a novel redox-dependent conformational change, the hydrogen bond formed between the carbonyl group of Gly57 and the N5H of the reduced cofactor strongly influences the stability of the semiquinone state. In this study, the sequential elimination of all side chain interactions in various combinations through a systematic alanine-scanning mutagenesis approach was conducted to more completely understand the functional inter-relationships as well as any synergistic interactions that might occur within the loop. In general, additive effects for each side chain on the midpoint potentials for both couples were observed except for the hydroquinone state where some degree of nonadditivity was noted in multiple mutants involving Glu59. The study concluded with the generation of the triple mutant -Ala56-Gly-Ala-Ala59- in which all side chain interactions are removed. Gly57 was left unchanged because of its critical conformational contribution. Remarkably, this mutant retained the ability to bind the FMN and to thermodynamically stabilize the semiquinone state despite the absence of all side chain interactions. Collectively, these observations emphasize the overriding importance of the main chain interactions with the N5H of the FMN and the associated redox-dependent conformational change in this loop and leaves little doubt as to its role in the thermodynamic stabilization of the neutral semiquinone state of the FMN cofactor.     

Protein dynamics: imidazole binding to class I C-type cytochromes.

The oxidized cytochrome c(2) from the purple phototrophic bacteria, Rhodobacter sphaeroides and Rhodobacter capsulatus, bind the neutral species of imidazole (K(a) = 1440 +/- 40 M(-1)) 50 times more strongly than does horse mitochondrial cytochrome c (K(a) = 30 +/- 1 M(-1)). The kinetics of imidazole binding are consistent with a change in rate-limiting step at high ligand concentrations for all three proteins. This is attributed to a conformational change leading to breakage of the iron-methionine bond which precedes imidazole binding. The three-dimensional structure of the Rb. sphaeroides cytochrome c(2) imidazole complex (Axelrod et al., Acta Crystalogr. D50, 596-602) supports the view that the conformational changes are essentially localized to approximately seven residues on either side of the ligated methionine and there is a hydrogen bond between the Phe 102 carbonyl, an internal water, and the bound imidazole. Insertions and deletions in this region of cytochrome c(2), the presence of a proline near the methionine, and the smaller size of the dynamic region of horse cytochrome c suggest that the stabilizing hydrogen bond is not present in horse cytochrome c, hence, the dramatic difference in affinity for imidazole. The kinetics of ligand binding do not correlate with either the strength of the iron-methionine bond as measured by the pK of the 695-nm absorption band or the overall stability of the cytochromes studied. However, the very similar imidazole binding properties of the two cytochromes c(2) indicate that the Rb. sphaeroides cytochrome c(2)-imidazole complex structure is an excellent model for the corresponding Rb. capsulatus cytochrome c(2) complex. It is notable that the movement of the peptide chain in the vicinity of the ligated methionine has been preserved throughout evolution and suggests a role in the function of c-type cytochromes.     

Evaluation of cooperative interactions between substructures of iso-1-cytochrome c using double mutant cycles

A double mutant cycle has been used to evaluate interaction energies between the global stabilizer mutation asparagine 52 --> isoleucine (N52I) in iso-1-cytochrome c and mutations producing single surface histidines at positions 26, 33, 39, 54, 73, 89, and 100. These histidine mutation sites are distributed through the four cooperative folding units of cytochrome c. The double mutant cycle starts with the iso-1-cytochrome c variant AcTM, a variant with no surface histidines and with asparagine at position 52. Isoleucine is added singly at position 52, AcTMI52 variant, as are the surface histidines, AcHX variants, where X indicates the histidine sequence position. The double mutant variants, AcHXI52, provide the remaining corner of the double mutant cycle. The stabilities of all variants were determined by guanidine hydrochloride denaturation and interaction energies were calculated between position 52 and each histidine site. Six of the seven double mutants show additive (AcH33I52, AcH39I52, AcH54I52, AcH89I52, and AcH100I52) stability effects or weak interaction energies (AcH73I52) of the histidine mutations and the N52I mutation, consistent with cooperative effects on protein folding and stability being sparsely distributed through the protein structure. The AcH26I52 variant shows a strong favorable interaction energy, 2.0 +/- 0.5 kcal/mol, between the N52I mutation in one substructure and the addition of His 26 to an adjacent substructure. The data are consistent with an entropic stabilization of the intersubstructure hydrogen bond between His 26 and Glu 44 by the Ile 52 mutation.     

Molecular Modeling of the NPY Binding Site on the Y1 Receptor

A three-dimensional model of the neuropeptide Y (NPY) - rat Y₁ (rY₁) receptor complex and of the NPY _13-36 - rY₁ receptor complex was constructed by molecular modeling based on the electron density projection map of rhodopsin and on site-directed mutagenesis studies of neuropeptide receptors. In order to further guide the modeling, the nucleotide sequences encoding Trp287, Cys295 and His297 in the third extracellular loop of the rY₁ receptor, were altered by site-directed mutagenesis experiments. Single-point mutated receptors were expressed in COS-7 cells, and tested for their ability to bind radio labelled NPY (³H-NPY). Mutations of Trp287 and His297 completely abolished binding of ³H-NPY. The Cys295Ser mutation only slightly decreased the binding of ³H-NPY, suggesting that the involvement of Cys295 in a disulphide bond is not essential for maintaining the correct three-dimensional structure of the binding site for NPY. Molecular dynamics simulations of NPY-rY₁ receptor interactions suggested that Asp199, Asp103 and Asp286 in the receptor interact, respectively, with Lys4, Arg33 and Arg35 of NPY. The simulations also suggested that His297 acts as a hydrogen acceptor from Arg35 in NPY, and that Tyr1 of NPY interacts with a binding pocket on the receptor formed by Asn115, Asp286, Trp287 and His297. Tyr36 in NPY interacted both with Thr41 and Tyr99 via hydrogen bonds, and also with Asn296, His297 and Phe301. The present study suggests that amino acid residues at the extracellular end of the transmembrane helices and in the extracellular loops are strongly involved in binding to NPY and NPY_13-36.Electronic Supplementary Material available.     

Major histocompatibility complex class II binding characteristics of peptoid-peptide hybrids

The major histocompatibility complex (MHC) class II binding requirements for solvent-exposed peptide residues were systematically studied using amino acid and peptoid substitutions. In a peptoid residue, the side chain is present on the backbone nitrogen atom as opposed to the alpha-carbon atom in an amino acid residue. To investigate the effect of this side chain shifting on MHC binding, three amino acids in the central part of the peptide sticking out of the binding groove were replaced by corresponding peptoid residues. Two peptoid-peptide hybrids showed large affinity decreases in the MHC-peptide binding assay. To investigate this affinity loss, the individual contributions to MHC binding affinity of the side chain (position), the putative hydrogen bond, and the flexibility were dissected. We conclude that the side chain position as well as the backbone nitrogen atom hydrogen bonding features of solvent-exposed residues in the peptide can be important for MHC binding affinity.