Pharmacophore-directed Homology Modeling and Molecular Dynamics Simulation of G Protein-coupled Receptor： Study of Possible Binding Modes of 5-HT2c Receptor Agonists

A new pharmacophore-based modeling procedure, including homology modeling, pharmacophore study, flexible molecular docking, and long-time molecular dynamics (MD) simulations, was employed to construct the structure of the human 5-HT_(2C) receptor and determine the characteristics of binding modes of 5-HT_(2C) receptor agonists. An agonist-receptor complex has been constructed based on homology modeling and a pharmacophore hypothesis model based on some high active compounds. Then MD simulations of the ligand-receptor complex in an explicit membrane environment were carried out. The conformation of the 5- HT_(2C) receptor during MD simulation was explored, and the stable binding modes of the studied agonist were determined. Flexible molecular docking of several structurally diverse agonists of the human 5-HT_(2C) receptor was carried out, and the general binding modes of these agonists were investigated. According to the models presented in this work and the results of Flexi-Dock, the involvement of the amino acid residues Asp134, Ser138, Ash210, Asn331, Tyr358, Ile131, Ser132, Val135, Thr139, Ile189, Val202, Val208, Leu209, Phe214, Val215, Gly218, Ser219, Phe223, Trp324, Phe327, and Phe328 in agonist recognition was studied. The obtained binding modes of the human 5-HT_(2C) receptor agonists have good agreement with the site-directed mutagenesis data and other studies.     

The unique extracellular disulfide loop of the glycine receptor is a principal ligand binding element.

A loop structure, formed by the putative disulfide bridging of Cys198 and Cys209, is a principal element of the ligand binding site in the glycine receptor (GlyR). Disruption of the loop's tertiary structure by Ser mutations of these Cys residues either prevented receptor assembly on the cell surface, or created receptors unable to be activated by agonists or to bind the competitive antagonist, strychnine. Mutation of residues Lys200, Tyr202 and Thr204 within this loop reduced agonist binding and channel activation sensitivities by up to 55-, 520- and 190-fold, respectively, without altering maximal current sizes, and mutations of Lys200 and Tyr202 abolished strychnine binding to the receptor. Removal of the hydroxyl moiety from Tyr202 by mutation to Phe profoundly reduced agonist sensitivity, whilst removal of the benzene ring abolished strychnine binding, thus demonstrating that Tyr202 is crucial for both agonist and antagonist binding to the GlyR. Tyr202 also influences receptor assembly on the cell surface, with only large chain substitutions (Phe, Leu and Arg, but not Thr, Ser and Ala) forming functional receptors. Our data demonstrate the presence of a second ligand binding site in the GlyR, consistent with the three-loop model of ligand binding to the ligand-gated ion channel superfamily.     

Probing the affinity and specificity of yeast alcohol dehydrogenase I for coenzymes.

Yeast (Saccharomyces cerevisiae) alcohol dehydrogenase I (SceADH) binds NAD+ and NADH less tightly and turns over substrates more rapidly than does horse (Equus caballus) liver alcohol dehydrogenase E isoenzyme (EcaADH), and neither enzyme uses NADP efficiently. Amino acid residues in the proposed adenylate binding pocket of SceADH were substituted in attempts to improve affinity for coenzymes or reactivity with NADP. Substitutions in SceADH (Gly202Ile or Ser246Ile) with the corresponding residues in the adenine binding site of the homologous EcaADH have modest effects on coenzyme binding and other kinetic constants, but the Ser246Ile substitution decreases turnover numbers by 350-fold. The Ser176Phe substitution (also near adenine site) significantly decreases affinity for coenzymes and turnover numbers. In the consensus nucleotide-binding betaalphabeta fold sequence, SceADH has two alanine residues (177-GAAGGLG-183) instead of the Leu200 in EcaADH (199-GLGGVG-204); the Ala178-Ala179 to Leu substitution significantly decreases affinity for coenzymes and turnover numbers. Some NADP-dependent enzymes have an Ala corresponding to Gly183 in SceADH; the Gly183Ala substitution significantly decreases affinity for coenzymes and turnover numbers. NADP-dependent enzymes usually have a neutral residue instead of the Asp (Asp201 in SceADH) that interacts with the hydroxyl groups of the adenosine ribose, along with a basic residue (at position 202 or 203) to stabilize the 2'-phosphate of NADP. The Gly203Arg change in SceADH does not significantly affect the kinetics. The Gly183Ala or Gly203Arg substitutions do not enable SceADH to use NADP+ as coenzyme. SceADH with the single Asp201Gly or double Asp201Gly:Gly203Arg substitutions have similar, low activity with NADP+. The results suggest that several of the amino acid residues participate in coenzyme binding and that conversion of specificity for coenzyme requires multiple substitutions.     

Identification of the critical residues of bradykinin receptor B1 for interaction with the kinins guided by site-directed mutagenesis and molecular modeling

We report the critical residues for the interaction of the kinins with human bradykinin receptor 1 (B1) using site-directed mutagenesis in conjunction with molecular modeling of the binding modes of the kinins in the homology model of the B1 receptor. Mutation of Lys118 in transmembrane (TM) helix 3, Ala270 in TM6, and Leu294 in TM7 causes a significant decrease in the affinity for the peptide agonists des-Arg10kallidin (KD) and des-Arg9BK but not the peptide antagonist des-Arg10Leu9KD. In contrast, mutations in TM2, TM3, TM6, and TM7 cause a significant decrease in the affinity for both the peptide agonists and the antagonist. These data indicate that the B1 bradykinin binding pocket for agonists and antagonists is similar, but the manners in which they interact with the receptor do not completely overlap. Therefore, there is a potential to influence the receptor's ligand selectivity.     

Identification of functional segments within the beta2I-domain of integrin alphaMbeta2

The alpha(M)beta(2) integrin plays an important role in leukocyte biology through its interactions with a diverse set of ligands. Efficient ligand binding requires the involvement of both the alpha(M) and beta(2) subunits. Past ligand binding studies have focused mainly on the alpha(M) subunit, with the beta(2) subunit being largely unexplored. Therefore, in this study we conducted homolog-scanning mutagenesis on the I-domain (residues 125-385) within the beta(2) subunit. We identified four noncontiguous sequences (Arg(144)-Lys(148), Gln(199)-Ala(203), Leu(225)-Leu(230), and Gly(305)-His(309)) that are critical for fibrinogen and C3bi binding to alpha(M)beta(2). Molecular modeling revealed that these four sequences reside within a narrow region on the surface of the beta(2)I-domain, in close proximity to three potential cation-binding sites. Among these sequences, Gln(199)-Ala(203), Leu(225)-Leu(230), and Gly(305)-His(309) are important for the binding of both ligands, whereas Arg(144)-Lys(148) is more critical for fibrinogen than for C3bi binding. These sequences within the beta(2)I-domain are directly involved in ligand binding, since 1) switching these segments to their corresponding beta(1) sequences destroyed ligand binding; 2) loss of function was not due to a nonspecific gross conformational change, since the defective alpha(M)beta(2) mutants reacted well with a panel of conformation-dependent mAbs; 3) mutation of these functional sequences did not effect Ca(2+) binding; and 4) synthetic peptides corresponding to sequences Gln(199)-Ala(203) and Gly(305)-His(309) blocked ligand binding to alpha(M)beta(2), and the peptides interacted directly with fibrinogen and C3bi. Given the similarity among all integrin beta subunits, our results may help us to understand the underlying mechanism of integrin-ligand interactions in general.     

Pyruvate carboxylase: isolation of the biotin-containing tryptic peptide and the determination of its primary sequency

The biotin-containing tryptic peptides of pyruvate carboxylase from sheep, chicken, and turkey liver mitochondria have been isolated and their primary structures determined. The amino acid sequences of the 19 residue peptides from chicken and turkey are identical and share a common sequence of 14 residues around biocytin with the 24-residue peptide isolated from sheep. The sequences obtained were: residue 1 → 11 Avian: Gly Ala Pro Leu Val Leu Ser Ala Met Biocytin Met Sheep: Gly Gln Pro Leu Val Leu Ser Ala Met Biocytin Met residues 12 → 19 or 24 Avian: Glu Thr Val Val Thr Ala Pro Arg Sheep: Glu Thr Val Val Thr Ser Pro Val Thr Glu Gly Val Arg A sensitive radiochemical assay for biotin was developed based on the tight binding of biotin by avidin. The ability of zinc sulfate to precipitate, without dissociating, the avidin-biotin complex provided a convenient procedure for separating free and bound biotin, and hence, for back-titrating a standard amount of avidin with [¹⁴C]biotin.     

The covalent structure of pig kidney fructose 1,6-bisphosphatase: sequence of the 60-residue NH2-terminal peptide produced by digestion with subtilisin

Digestion of the native pig kidney fructose 1,6-bisphosphatase tetramer with subtilisin cleaves each of the 35,000-molecular-weight subunits to yield two major fragments: the S-subunit (M_{r ca.} 29,000), and the S-peptide (M_r 6,500). The following amino acid sequence has been determined for the S peptide: AcThrAspGlnAlaAlaPheAspThrAsnIle Val ThrLeuThrArgPheValMetGluGlnGlyArgLysAla ArgGlyThrGlyGlu MetThrGlnLeuLeuAsnSerLeuCysThrAlaValLys AlaIleSerThrAla z.sbnd;ValArgLysAlaGlyIleAlaHisLeuTyrGlyIleAla. Comparison of this sequence with that of the NH₂-terminal 60 residues of the enzyme from rabbit liver (El-Dorry et al., 1977, Arch. Biochem. Biophys.182, 763) reveals strong homology with 52 identical positions and absolute identity in sequence from residues 26 to 60.Although subtilisin cleavage of fructose 1,6-bisphosphatase results in diminished sensitivity of the enzyme to AMP inhibition, we have found no AMP inhibition-related amino acid residues in the sequenced S-peptide. The loss of AMP sensitivity that occurs upon pyridoxal-P modification of the enzyme does not result in the modification of lysyl residues in the S-peptide. Neither photoaffinity labeling of fructose 1,6-bisphosphatase with 8-azido-AMP nor modification of the cysteinyl residue proximal to the AMP allosteric site resulted in the modification of residues located in the NH₂-terminal 60-amino acid peptide.     

Molecular dynamics and free energy studies on the Drosophila melanogaster and Leptinotarsa decemlineata ecdysone receptor complexed with agonists: Mechanism for binding and selectivity

The ecdysone receptor is a nuclear hormone receptor that plays a pivotal role in the insect metamorphosis and development. To address the molecular mechanisms of binding and selectivity, the interactions of two typical agonists Ponasterone A and 20-Hydroxyecdysone with Drosophila melanogaster (DME) and Leptinotarsa decemlineata ecdysone (LDE) receptors were investigated by homology modeling, molecular docking, molecular dynamic simulation, and thermodynamic analysis. We discover that 1) the L5-loop, L11-loop, and H12 helix for DME, L7-loop, and L11-loop for LDE are more flexible, which affect the global dynamics of the ligand-binding pocket, thus facilitating the ligand recognition of ecdysone receptor; 2) several key residues (Thr55/Thr37, Phe109/Phe91, Arg95/Arg77, Arg99/Arg81, Phe108/Leu90, and Ala110/Val92) are responsible for the binding of the proteins; 3) the binding-free energy is mainly contributed by the van der Waals forces as well as the electrostatic interactions of ligand and receptor; 4) the computed binding-free energy difference between DME-C1 and LDE-C1 is –4.65 kcal/mol, explains that C1 can form many more interactions with the DME; 5) residues Phe108/Leu90 and Ala110/Val92 have relatively position and orientation difference in the two receptors, accounting most likely for the ligand selectivity of ecdysone receptor from different orders of insects. This study underscores the expectation that different insect pests should be able to discriminate among compounds from different as yet undiscovered compounds, and the results firstly show a structural and functional relay between the agonists and receptors (DME and LDE), which can provide an avenue for the development of target-specific insecticides.

Communicated by Ramaswamy H. Sarma     

Importance of each residue within secretin for receptor binding and biological activity

Secretin is a linear 27-residue peptide hormone that stimulates pancreatic and biliary ductular bicarbonate and water secretion by acting at its family B G protein-coupled receptor. While, like other family members, the carboxyl-terminal region of secretin is most important for high affinity binding and its amino-terminal region is most important for receptor selectivity and receptor activation, determinants for these activities are distributed throughout the entire length of this peptide. In this work, we have systematically investigated changing each residue within secretin to alanine and evaluating the impact on receptor binding and biological activity. The residues most critical for receptor binding were His1, Asp3, Gly4, Phe6, Thr7, Ser8, Leu10, Asp15, Leu19, and Leu23. The residues most critical for biological activity included His1, Gly4, Thr7, Ser8, Glu9, Leu10, Leu19, Leu22, and Leu23, with Asp3, Phe6, Ser11, Leu13, Asp15, Leu26, and Val27 also contributing. While the importance of residues in positions analogous to His1, Asp3, Phe6, Thr7, and Leu23 is conserved for several closely related members of this family, Leu19 is uniquely important for secretin. We, therefore, have further studied this residue by molecular modeling and molecular dynamics simulations. Indeed, the molecular dynamics simulations showed that mutation of Leu19 to alanine was destabilizing, with this effect greater than that observed for the analogous position in the other close family members. This could reflect reduced contact with the receptor or an increase in the solvent-accessible surface area of the hydrophobic residues in the carboxyl terminus of secretin as bound to its receptor.     

A single residue (arg46) located within the N-terminus of the V1a vasopressin receptor is critical for binding vasopressin but not peptide or nonpeptide antagonists

A fundamental issue in molecular endocrinology is to define how agonist:receptor interaction differs from antagonist:receptor interaction. The vasopressin V1a receptor (V1aR) is a member of a subfamily of related G protein-coupled receptors that are activated by the hormone AVP or related peptides. The N-terminus of the V1aR has recently been shown to be critical for binding agonists but not antagonists. Using a combination of N-terminally truncated constructs and alanine-scanning mutagenesis, individual residues that provide these agonist-specific binding epitopes have now been identified in this study. Our data establish that a single residue, Arg46, is critical for AVP binding to the V1aR. Systematic substitution revealed that Arg was required at this locus and could not be substituted by Lys, Glu, Leu, or Ala. In contrast, antagonist binding (cyclic or linear, peptide or nonpeptide) was unaffected. Disruption of Arg46 also resulted in defective intracellular signaling. Arginine is conserved at this locus in all members of the neurohypophysial peptide hormone receptor family cloned to date, indicative of a fundamental role in receptor function. In addition to Arg46, the residues Leu42, Gly43, Asp45 form a patch contributing to AVP binding. This study provides molecular insight into the role of the V1aR N-terminus and key differences between agonist and antagonist binding requirements.     

Strategies for Positioning Fluorescent Probes and Crosslinkers on Formyl Peptide Ligands

Abstract

Chemoattractant receptors represent a major subset of the G-protein coupled receptor (GPCR) family. One of the best characterized, the N-formyl peptide receptor (FPR), participates in host defense responses of neutrophils. The features of the ligand which regulate its interaction with the FPR are well-known. By manipulating these features we have developed new ligands to probe structural and mechanistic aspects of the peptide-receptor interaction. Three ligand groups have been developed: 1) ligands containing a Lys residue located in positions 2 through 7 that can be conjugated to FITC (N-formyl-Met1-Lys2-Phe3-Phe4, N-formyl-Met1-Leu2-Lys3-Phe4, N-formyl-Met1-Leu2-Phe3-Lys4, N-formyl-Met1-Leu2-Phe3-Phe4-Lys5, N-formyl-nLeu1-Leu2-Phe3-nLeu4-Tyr5-Lys6 and N-formyl-Met1-Leu2-Phe3-Phe4-Gly5-Gly6-Lys7; 2) fluorescent pentapeptide ligands (N-formyl-Met-X-Phe-Phe-Lys(FITC) where X = Leu, Ala, Val or Gly); and 3) small crosslinking ligands where the photoaffinity crosslinker 4-azidosalicylic acid (ASA) was conjugated to Lys in positions 3 and 4 and p-benzoyl-phenylalanine (Bpa) was located in position 2 in N-formyl-Met1-Bpa2-Phe3-Tyr4. The peptides were characterized according to activity and affinity in human neutrophils and cell lines transfected with FPR. All of the peptides were agonists, with parallel affinity and activity. In the first group, the peptide activity decreases as Lys is placed closer to the N-formyl group and the activity is improved by 1–3 orders of magnitude by conjugation with FITC. In the second group, the dissociation rate of the peptide from the receptor increases as position 2 is replaced by aliphatic amino acids with smaller alkyl groups. In the third group, crosslinking ligands remain biologically active, display nM affinity and covalently label the FPR.     

Targeting of hepatocellular carcinoma with glypican‐3‐targeting peptide ligand

Hepatocellular carcinoma is a common malignancy. The carcinoma cells express glypican‐3 (GPC‐3) on the cell membrane. GPC‐3 is also expressed in melanoma cells. Therefore, GPC‐3 might be a potential target for tumor imaging or therapy. Here, proteomic mass spectrometry was used to identify peptides that target GPC‐3‐expressing tumors. A mammalian expression vector expressing a FLAG‐GPC‐3 fusion protein was cloned for immunoprecipitation. With the use of liposomes, the vector was transfected into HepG2 (HepG2/FLAG‐GPC‐3) and HEK 293 cells, and the transfected cell lines were selected with geneticin. HepG2/FLAG‐GPC‐3 cells were used for immunoprecipitation of FLAG‐GPC‐3 fusion protein. Seven peptide candidates (L1–L7) were selected for GPC‐3‐targeting ligands by mass spectrometric analysis. The L5 peptide with 14 amino acids (Arg‐Leu‐Asn‐Val‐Gly‐Gly‐Thr‐Tyr‐Phe‐Leu‐Thr‐Thr‐Arg‐Gln) showed selective binding to the GPC‐3‐expressing tumor cells, as did a shortened L5 peptide (L5‐2) with seven amino acids (Tyr‐Phe‐Leu‐Thr‐Thr‐Arg‐Gln). These peptide ligands have potential as targeting moieties to GPC‐3‐expressing tumors for diagnostic and/or therapeutic purposes. Copyright © 2011 European Peptide Society and John Wiley & Sons, Ltd.     

X-ray structure of papaya chitinase reveals the substrate binding mode of glycosyl hydrolase family 19 chitinases

The crystal structure of a chitinase from Carica papaya has been solved by the molecular replacement method and is reported to a resolution of 1.5 A. This enzyme belongs to family 19 of the glycosyl hydrolases. Crystals have been obtained in the presence of N-acetyl- d-glucosamine (GlcNAc) in the crystallization solution and two well-defined GlcNAc molecules have been identified in the catalytic cleft of the enzyme, at subsites -2 and +1. These GlcNAc moieties bind to the protein via an extensive network of interactions which also involves many hydrogen bonds mediated by water molecules, underlying their role in the catalytic mechanism. A complex of the enzyme with a tetra-GlcNAc molecule has been elaborated, using the experimental interactions observed for the bound GlcNAc saccharides. This model allows to define four major substrate interacting regions in the enzyme, comprising residues located around the catalytic Glu67 (His66 and Thr69), the short segment E89-R90 containing the second catalytic residue Glu89, the region 120-124 (residues Ser120, Trp121, Tyr123, and Asn124), and the alpha-helical segment 198-202 (residues Ile198, Asn199, Gly201, and Leu202). Water molecules from the crystal structure were introduced during the modeling procedure, allowing to pinpoint several additional residues involved in ligand binding that were not previously reported in studies of poly-GlcNAc/family 19 chitinase complexes. This work underlines the role played by water-mediated hydrogen bonding in substrate binding as well as in the catalytic mechanism of the GH family 19 chitinases. Finally, a new sequence motif for family 19 chitinases has been identified between residues Tyr111 and Tyr125.     

Identification of residues important for ligand binding of thromboxane A2 receptor in the second extracellular loop using the NMR experiment-guided mutagenesis approach

The second extracellular loop (eLP2) of the thromboxane A(2) receptor (TP) had been proposed to be involved in ligand binding. Through two-dimensional (1)H NMR experiments, the overall three-dimensional structure of a constrained synthetic peptide mimicking the eLP2 had been determined by our group (Ruan, K.-H., So, S.-P., Wu, J., Li, D., Huang, A., and Kung, J. (2001) Biochemistry 40, 275-280). To further identify the residues involved in ligand binding, a TP receptor antagonist, SQ29,548 was used to interact with the synthetic peptide. High resolution two-dimensional (1)H NMR experiments, NOESY, and TOCSY were performed for the peptide, SQ29,548, and peptide with SQ29,548, respectively. Through completed (1)H NMR assignment and by comparing the different spectra, extra peaks were observed on the NOESY spectrum of the peptide with SQ29,548, which implied the contacts between residues of eLP2 at Val(176), Leu(185), Thr(186), and Leu(187) with SQ29,548 at position H2, H7, and H8. Site-directed mutagenesis was used to confirm the possible ligand-binding sites on native human TP receptor. Each of the four residues was mutated to the residues either in the same group, with different structure or different charged. The mutated receptors were then tested for their ligand binding activity. The receptor with V176L mutant retained binding activity to SQ29,548. All other mutations resulted in decreased or lost binding activity to SQ29,548. These mutagenesis results supported the prediction from NMR experiments in which Val(176), Leu(185), Thr(186), and Leu(187) are the possible residues involved in ligand binding. This information facilitates the understanding of the molecular mechanism of thromboxane A(2) binding to the important receptor and its signal transduction.     

The amino acid sequence of human chorionic gonadotropin. The alpha subunit and beta subunit.

The amino acid sequences of both the alpha and beta subunits of human chorionic gonadotropin have been determined. The amino acid sequence of the alpha subunit is: Ala - Asp - Val - Gln - Asp - Cys - Pro - Glu - Cys-10 - Thr - Leu - Gln - Asp - Pro - Phe - Ser - Gln-20 - Pro - Gly - Ala - Pro - Ile - Leu - Gln - Cys - Met - Gly-30 - Cys - Cys - Phe - Ser - Arg - Ala - Tyr - Pro - Thr - Pro-40 - Leu - Arg - Ser - Lys - Lys - Thr - Met - Leu - Val - Gln-50 - Lys - Asn - Val - Thr - Ser - Glu - Ser - Thr - Cys - Cys-60 - Val - Ala - Lys - Ser - Thr - Asn - Arg - Val - Thr - Val-70 - Met - Gly - Gly - Phe - Lys - Val - Glu - Asn - His - Thr-80 - Ala - Cys - His - Cys - Ser - Thr - Cys - Tyr - Tyr - His-90 - Lys - Ser. Oligosaccharide side chains are attached at residues 52 and 78. In the preparations studied approximately 10 and 30% of the chains lack the initial 2 and 3 NH2-terminal residues, respectively. This sequence is almost identical with that of human luteinizing hormone (Sairam, M. R., Papkoff, H., and Li, C. H. (1972) Biochem. Biophys. Res. Commun. 48, 530-537). The amino acid sequence of the beta subunit is: Ser - Lys - Glu - Pro - Leu - Arg - Pro - Arg - Cys - Arg-10 - Pro - Ile - Asn - Ala - Thr - Leu - Ala - Val - Glu - Lys-20 - Glu - Gly - Cys - Pro - Val - Cys - Ile - Thr - Val - Asn-30 - Thr - Thr - Ile - Cys - Ala - Gly - Tyr - Cys - Pro - Thr-40 - Met - Thr - Arg - Val - Leu - Gln - Gly - Val - Leu - Pro-50 - Ala - Leu - Pro - Gin - Val - Val - Cys - Asn - Tyr - Arg-60 - Asp - Val - Arg - Phe - Glu - Ser - Ile - Arg - Leu - Pro-70 - Gly - Cys - Pro - Arg - Gly - Val - Asn - Pro - Val - Val-80 - Ser - Tyr - Ala - Val - Ala - Leu - Ser - Cys - Gln - Cys-90 - Ala - Leu - Cys - Arg - Arg - Ser - Thr - Thr - Asp - Cys-100 - Gly - Gly - Pro - Lys - Asp - His - Pro - Leu - Thr - Cys-110 - Asp - Asp - Pro - Arg - Phe - Gln - Asp - Ser - Ser - Ser - Ser - Lys - Ala - Pro - Pro - Pro - Ser - Leu - Pro - Ser-130 - Pro - Ser - Arg - Leu - Pro - Gly - Pro - Ser - Asp - Thr-140 - Pro - Ile - Leu - Pro - Gln. Oligosaccharide side chains are found at residues 13, 30, 121, 127, 132, and 138. The proteolytic enzyme, thrombin, which appears to cleave a limited number of arginyl bonds, proved helpful in the determination of the beta sequence.     

The primary structure of the β-subunit of protocatechuate 3,4-dioxygenase from Pseudomonas aeruginosa

The complete amino acid sequence of the β-subunit of protocatechuate 3,4-dioxygenase was determined. The β-subunit contained four methionine residues. Thus, five peptides were obtained after cleavage of the carboxymethylated β-subunit with cyanogen bromide, and were isolated on Sephadex G-75 column chromatography. The amino acid sequences of the cyanogen bromide peptides were established by characterization of the peptides obtained after digestion with trypsin, chymotrypsin, thermolysin, or Staphylococcus aureus protease. The major sequencing techniques used were automated and manual Edman degradations. The five cyanogen bromide peptides were aligned by means of the amino acid sequences of the peptides containing methionine purified from the tryptic hydrolysate of the carboxymethylated β-subunit. The amino acid sequence of all the 238 residues was as follows: ProAlaGlnAspAsnSerArgPheValIleArgAsp ArgAsnTrpHis ProLysAlaLeuThrPro-Asp — TyrLysThrSerIleAlaArg SerProArgGlnAla LeuValSerIleProGlnSer — IleSerGluThrThrGly ProAsnPheSerHisLeu GlyPheGlyAlaHisAsp-His — AspLeuLeuLeuAsnPheAsn AsnGlyGlyLeu ProIleGlyGluArgIle-Ile — ValAlaGlyArgValValAsp GlnTyrGlyLysPro ValProAsnThrLeuValGluMet — TrpGlnAlaAsnAla GlyGlyArgTyrArg HisLysAsnAspArgTyrLeuAlaPro — LeuAspProAsn PheGlyGlyValGly ArgCysLeuThrAspSerAspGlyTyrTyr — SerPheArg ThrIleLysProGlyPro TyrProTrpArgAsnGlyProAsnAsp — TrpArgProAla HisIleHisPheGlyIle SerGlyProSerIleAlaThr-Lys — LeuIleThrGlnLeuTyr PheGluGlyAspPro LeuIleProMetCysProIleVal — LysSerIleAlaAsn ProGluAlaValGlnGln LeuIleAlaLysLeuAspMetAsnAsn — AlaAsnProMet AsnCysLeuAlaTyr ArgPheAspIleValLeuArgGlyGlnArgLysThrHis PheGluAsnCys. The sequence published earlier in summary form (Iwaki et al., 1979, J. Biochem.86, 1159–1162) contained a few errors which are pointed out in this paper.     

Molecular dynamics of NPY Y1 receptor activation

A three-dimensional model of the human neuropeptide Y(NPY)Y₁ receptor (hY₁) was constructed, energy refined and used to simulate molecular receptor interactions of the peptide ligands NPY, [L31, P34]NPY, peptide YY (PYY) and pancreatic polypeptide (PP), and of the nonpeptide antagonist R-N²-(diphenylacetyl)-N-(4-hydroxyphenyl)methyl-argininamide (BIBP3226) and its S-enantiomer BIBP3435. The best complementarity in charges between the receptor and the peptides, and the best structural accordance with experimental studies, was obtained with amino acid 1–4 of the peptides interacting with Asp194, Asp200, Gln201, Phe202 and Trp288 in the receptor. Arg33 and Arg35 of the peptides formed salt bridges with Asp104 and Asp287, respectively, while Tyr36 interacted in a binding pocket formed by Phe41, Thr42, Tyr100, Asn297, His298 and Phe302. Calculated electrostatic potentials around NPY and hY₁ molecules indicated that ligand binding is initiated by electrostatic interactions between a highly positive region in the N- and C-terminal parts of the peptides, and a negative region in the extracellular receptor domains. Molecular dynamics simulations of NPY and BIBP3226 interactions with the receptor indicated rigid body motions of TMH5 and TMH6 upon NPY binding as mechanisms of receptor activation, and that BIBP3226 may act as an antagonist by constraining these motions.     

Parallel analysis of mutant human glucose 6-phosphate dehydrogenase in yeast using PCR colonies

We demonstrate a highly parallel strategy to analyze the impact of single nucleotide mutations on protein function. Using our method, it is possible to screen a population and quickly identify a subset of functionally interesting mutants. Our method utilizes a combination of yeast functional complementation, growth competition of mutant pools, and polymerase colonies. A defined mutant human glucose-6-phosphate-dehydrogenase library was constructed which contains all possible single nucleotide missense mutations in the eight-residue glucose-6-phosphate binding peptide of the enzyme. Mutant human enzymes were expressed in a zwf1 (gene encoding yeast homologue) deletion strain of Saccharomyces cerevisiae. Growth rates of the 54 mutant strains arising from this library were measured in parallel in conditions selective for active hG6PD. Several residues were identified which tolerated no mutations (Asp200, His201 and Lys205) and two (Ile199 and Leu203) tolerated several substitutions. Arg198, Tyr202, and Gly204 tolerated only 1-2 specific substitutions. Generalizing from the positions of tolerated and non-tolerated amino acid substitutions, hypotheses were generated about the functional role of specific residues, which could, potentially, be tested using higher resolution/lower throughput methods.     

Study of the binding residues between ANEPII and insect sodium channel receptor

The present study aimed at determining the functional characteristics of anti-neuroexcitation peptide II (ANEPII). The depressant insect toxin ANEPII from the Chinese scorpion Buthus martensii Karsch had an effect on insect sodium channels. Previous studies showed that scorpion depressant toxins induce insect flaccid paralysis upon binding to receptor site-4, so we tried to predict the functional residues involved using computational techniques. In this study, three-dimensional structure modeling of ANEPII and site-4 of the insect sodium channel were carried out by homology modeling, and these models were used as the starting point for nanosecond-duration molecular dynamics simulations. Docking studies of ANEPII in the sodium channel homology model were conducted, and likely ANEPII binding loci were investigated. Based on these analyses, the residues Tyr34, Trp36, Gly39, Leu40, Trp53, Asn58, Gly61 and Gly62 were predicted to interact with sodium channel receptor and to act as functional residues.     

Evolutionarily conserved residues at glucagon-like peptide-1 (GLP-1) receptor core confer ligand-induced receptor activation

Glucagon-like peptide-1 (GLP-1) and glucose-dependent insulinotropic polypeptide (GIP) play important roles in insulin secretion through their receptors, GLP1R and GIPR. Although GLP-1 and GIP are attractive candidates for treatment of type 2 diabetes and obesity, little is known regarding the molecular interaction of these peptides with the heptahelical core domain of their receptors. These core domains are important not only for specific ligand binding but also for ligand-induced receptor activation. Here, using chimeric and point-mutated GLP1R/GIPR, we determined that evolutionarily conserved amino acid residues such as Ile(196) at transmembrane helix 2, Leu(232) and Met(233) at extracellular loop 1, and Asn(302) at extracellular loop 2 of GLP1R are responsible for interaction with ligand and receptor activation. Application of chimeric GLP-1/GIP peptides together with molecular modeling suggests that His(1) of GLP-1 interacts with Asn(302) of GLP1R and that Thr(7) of GLP-1 has close contact with a binding pocket formed by Ile(196), Leu(232), and Met(233) of GLP1R. This study may provide critical clues for the development of peptide and/or nonpeptide agonists acting at GLP1R.