Sequence-specific 1H NMR assignments and solution structure of bovine pancreatic polypeptide

Targeting type I DED interactions at the DED filament serves as a sensitive switch for cell fate decisions

Activation of procaspase-8 in the death effector domain (DED) filaments of the death-inducing signaling complex (DISC) is a key step in apoptosis. In this study, a rationally designed cell-penetrating peptide, DEDid, was engineered to mimic the h2b helical region of procaspase-8-DED2 containing a highly conservative FL motif. Furthermore, mutations were introduced into the DEDid binding site of the procaspase-8 type I interface. Additionally, our data suggest that DEDid targets other type I DED interactions such as those of FADD. Both approaches of blocking type I DED interactions inhibited CD95L-induced DISC assembly, caspase activation and apoptosis. We showed that inhibition of procaspase-8 type I interactions by mutations not only diminished procaspase-8 recruitment to the DISC but also destabilized the FADD core of DED filaments. Taken together, this study offers insights to develop strategies to target DED proteins, which may be considered in diseases associated with cell death and inflammation.

Illuminating the neuropeptide Y4 receptor and its ligand pancreatic polypeptide from a structural, functional, and therapeutic perspective

The neuropeptide Y4 receptor (Y4R), a rhodopsin-like G protein-coupled receptor (GPCR) and the hormone pancreatic polypeptide (PP) are members of the neuropeptide Y family consisting of four receptors (Y1R, Y2R, Y4R, Y5R) and three highly homologous peptide ligands (neuropeptide Y, peptide YY, PP). In this family, the Y4R is of particular interest as it is the only subtype with high affinity to PP over NPY. The Y4R, as a mediator of PP signaling, has a pivotal role in appetite regulation and energy homeostasis, offering potential avenues for the treatment of metabolic disorders such as obesity. PP as anorexigenic peptide is released postprandial from the pancreas in response to food intake, induces satiety signals and contributes to hamper excessive food intake. Moreover, this system was also described to be associated with different types of cancer: overexpression of Y4R have been found in human adenocarcinoma cells, while elevated levels of PP are related to the development of pancreatic endocrine tumors. The pharmacological relevance of the Y4R advanced the search for potent and selective ligands for this receptor subtype, which will be significantly progressed through the elucidation of the active state PP-Y4R cryo-EM structure. This review summarizes the development of novel PP-derived ligands, like Obinepitide as dual Y2R/Y4R agonist in clinical trials or UR-AK86c as small hexapeptide agonist with picomolar affinity, as well as the first allosteric modulators that selectively target the Y4R, e.g. VU0506013 as potent Y4R positive allosteric modulator or (S)-VU0637120 as allosteric antagonist. Here, we provide valuable insights into the complex physiological functions of the Y4R and PP and the pharmacological relevance of the system in appetite regulation to open up new avenues for the development of tool compounds for targeted therapies with potential applications in metabolic disorders.

Structural basis for Y2 receptor-mediated neuropeptide Y and peptide YY signaling

Neuropeptide Y (NPY) and its receptors are expressed in various human tissues including the brain where they regulate appetite and emotion. Upon NPY stimulation, the neuropeptide Y1 and Y2 receptors (Y₁R and Y₂R, respectively) activate G_I signaling, but their physiological responses to food intake are different. In addition, deletion of the two N-terminal amino acids of peptide YY (PYY(3-36)), the endogenous form found in circulation, can stimulate Y₂R but not Y₁R, suggesting that Y₁R and Y₂R may have distinct ligand-binding modes. Here, we report the cryo-electron microscopy structures of the PYY(3-36)‒Y₂R‒G_i and NPY‒Y₂R‒G_i complexes. Using cell-based assays, molecular dynamics simulations, and structural analysis, we revealed the molecular basis of the exclusive binding of PYY(3-36) to Y₂R. Furthermore, we demonstrated that Y₂R favors G protein signaling over β-arrestin signaling upon activation, whereas Y₁R does not show a preference between these two pathways.

Crystal structures of human neuropeptide Y (NPY) and peptide YY (PYY)

Neuropeptide Y (NPY), peptide YY (PYY) and pancreatic polypeptide (PP) form the evolutionarily conserved pancreatic polypeptide family. While the fold is widely utilized in nature, crystal structures remain elusive, particularly for the human forms, with only the structure of a distant avian form of PP reported. Here we utilize a crystallization chaperone (antibody Fab fragment), specifically recognizing the amidated peptide termini, to solve the structures of human NPY and human PYY. Intriguingly, and despite limited sequence identity (~50%), the structure of human PYY closely resembles that of avian PP, highlighting the broad structural conservation of the fold throughout evolution. Specifically, the PYY structure is characterized by a C-terminal amidated α-helix, preceded by a backfolded poly-proline N-terminus, with the termini in close proximity to each other. In contrast, in the structure of human NPY the N-terminal component is disordered, while the helical component of the peptide is observed in a four-helix bundle type arrangement, consistent with a propensity for multimerization suggested by NMR studies.

Arrestin-dependent internalization of rhodopsin-like G protein-coupled receptors

The internalization of G protein-coupled receptors (GPCRs) is an important mechanism regulating the signal strength and limiting the opportunity of receptor activation. Based on the importance of GPCRs, the detailed knowledge about the regulation of signal transduction is crucial. Here, current knowledge about the agonist-induced, arrestin-dependent internalization process of rhodopsin-like GPCRs is reviewed. Arrestins are conserved molecules that act as key players within the internalization process of many GPCRs. Based on highly conserved structural characteristics within the rhodopsin-like GPCRs, the identification of arrestin interaction sites in model systems can be compared and used for the investigation of internalization processes of other receptors. The increasing understanding of this essential regulation mechanism of receptors can be used for drug development targeting rhodopsin-like GPCRs. Here, we focus on the neuropeptide Y receptor family, as these receptors transmit various physiological processes such as food intake, energy homeostasis, and regulation of emotional behavior, and are further involved in pathophysiological processes like cancer, obesity and mood disorders. Hence, this receptor family represents an interesting target for the development of novel therapeutics requiring the understanding of the regulatory mechanisms influencing receptor mediated signaling.

Artificial Enzymes for Diels-Alder Reactions

The Diels-Alder (DA) reaction is a cycloaddition of a conjugated diene and an alkene (dienophile) leading to the formation of a cyclohexene derivative through a concerted mechanism. As DA reactions generally proceed with a high degree of regio- and stereoselectivity, they are widely used in synthetic organic chemistry. Considering eco-conscious public and governmental movements, efforts are now directed towards the development of synthetic processes that meet environmental concerns. Artificial enzymes, which can be developed to catalyze abiotic reactions, appear to be important synthetic tools in the synthetic biology field. This review describes the different strategies used to develop protein-based artificial enzymes for DA reactions, including for in cellulo approaches.

Stabilizing and Understanding a Miniprotein by Rational Redesign

Miniproteins reduce the complexity of the protein-folding problem allowing systematic studies of contributions to protein folding and stabilization. Here, we describe the rational redesign of a miniprotein, PPα, comprising a polyproline II helix, a loop, and an α helix. The redesign provides a de novo framework for interrogating noncovalent interactions. Optimized PPα has significantly improved thermal stability with a midpoint unfolding temperature (T_M) of 51 °C. Its nuclear magnetic resonance structure indicates a density of stabilizing noncovalent interactions that is higher than that of the parent peptide, specifically an increased number of CH-π interactions. In part, we attribute this to improved long-range electrostatic interactions between the two helical elements. We probe further sequence-stability relationships in the miniprotein through a series of rational mutations.

Structural Capacitance in Protein Evolution and Human Diseases

Canonical mechanisms of protein evolution include the duplication and diversification of pre-existing folds through genetic alterations that include point mutations, insertions, deletions, and copy number amplifications, as well as post-translational modifications that modify processes such as folding efficiency and cellular localization. Following a survey of the human mutation database, we have identified an additional mechanism that we term "structural capacitance," which results in the de novo generation of microstructure in previously disordered regions. We suggest that the potential for structural capacitance confers select proteins with the capacity to evolve over rapid timescales, facilitating saltatory evolution as opposed to gradualistic canonical Darwinian mechanisms. Our results implicate the elements of protein microstructure generated by this distinct mechanism in the pathogenesis of a wide variety of human diseases. The benefits of rapidly furnishing the potential for evolutionary change conferred by structural capacitance are consequently counterbalanced by this accompanying risk. The phenomenon of structural capacitance has implications ranging from the ancestral diversification of protein folds to the engineering of synthetic proteins with enhanced evolvability.

Adaptive Steered Molecular Dynamics of the Long-Distance Unfolding of Neuropeptide Y

Neuropeptide Y (NPY) has been found to adopt two stable conformations in vivo: (1) a monomeric form called the PP-fold in which a polyproline tail is folded onto an α-helix via a β-turn and (2) a dimeric form of the unfolded proteins in which the α-helices interact with each other via side chains. The transition pathway and rates between the two conformations remain unknown and are important to the nature of the binding of the protein. Toward addressing this question, the present work suggests that the unfolding of the PP-fold is too slow to play a role in NPY monomeric binding unless the receptor catalyzes it to do so. Specifically, the dynamics and structural changes of the unfolding of a monomeric NPY protein have been investigated in this work. Temperature accelerated molecular dynamics (MD) simulations at 500 K under constant (N,V,E) conditions suggests a hinge-like unraveling of the tail rather than a random unfolding. The free energetics of the proposed unfolding pathway have been described using an adaptive steered MD (SMD) approach at various temperatures. This approach generalizes the use of Jarzynski's equality through a series of stages that allows for better convergence along nonlinear and long-distance pathways. Results acquired using this approach provide a potential of mean force (PMF) with narrower error bars and are consistent with some of the earlier reports on the qualitative behavior of NPY binding.

Non-specific binding and general cross-reactivity of Y receptor agonists are correlated and should importantly depend on their acidic sectors

Non-specific binding of Y receptor agonists to intact CHO cells, and to CHO cell or rat brain particulates, is much greater for human neuropeptide Y (hNPY) compared to porcine peptide Y (pPYY), and especially relative to human pancreatic polypeptide (hPP). This binding of hNPY is reduced by alkali cations in preference to non-ionic chaotrope urea, while the much lower non-specific binding of pPYY is more sensitive to urea. The difference could mainly be due to the 10-16 stretch in 36-residue Y agonists (residues 8-14 in N-terminally clipped 34-peptides), located in the sector that contains all acidic residues of physiological Y agonists. Anionic pairs containing aspartate in the 10-16 zone could be principally responsible for non-specific attachments, but may also aid the receptor site binding. Two such pairs are found in hNPY, one in pPYY, and none in hPP. The hydroxyl amino acid residue at position 13 in mammalian PYY and PP molecules could lower conformational plasticity and the non-selective binding via intrachain hydrogen bonding. The acidity of this tract could also be important in agonist selectivity of the Y receptor subtypes. The differences point to an evolutionary reduction of promiscuous protein binding from NPY to PP, and should also be important for Y agonist selectivity within NPY receptor group, and correlate with partial agonism and out-of group cross-reactivity with other receptors.

Discriminating the native structure from decoys using scoring functions based on the residue packing in globular proteins

BACKGROUND

Setting the rules for the identification of a stable conformation of a protein is of utmost importance for the efficient generation of structures in computer simulation. For structure prediction, a considerable number of possible models are generated from which the best model has to be selected.

RESULTS

Two scoring functions, Rs and Rp, based on the consideration of packing of residues, which indicate if the conformation of an amino acid sequence is native-like, are presented. These are defined using the solvent accessible surface area (ASA) and the partner number (PN) (other residues that are within 4.5 A) of a particular residue. The two functions evaluate the deviation from the average packing properties (ASA or PN) of all residues in a polypeptide chain corresponding to a model of its three-dimensional structure. While simple in concept and computationally less intensive, both the functions are at least as efficient as any other energy functions in discriminating the native structure from decoys in a large number of standard decoy sets, as well as on models submitted for the targets of CASP7. Rs appears to be slightly more effective than Rp, as determined by the number of times the native structure possesses the minimum value for the function and its separation from the average value for the decoys.

CONCLUSION

Two parameters, Rs and Rp, are discussed that can very efficiently recognize the native fold for a sequence from an ensemble of decoy structures. Unlike many other algorithms that rely on the use of composite scoring function, these are based on a single parameter, viz., the accessible surface area (or the number of residues in contact), but still able to capture the essential attribute of the native fold.

Towards understanding the free and receptor bound conformation of neuropeptide Y by fluorescence resonance energy transfer studies

Despite a considerable sequence identity of the three mammalian hormones of the neuropeptide Y family, namely neuropeptide Y, peptide YY and pancreatic polypeptide, their structure in solution is described to be different. A so-called pancreatic polypeptide-fold has been identified for pancreatic polypeptide, whereas the structure of the N-terminal segment of neuropeptide Y is unknown. This element is important for the binding of neuropeptide Y to two of its relevant receptors, Y(1) and Y(5), but not to the Y(2) receptor subtype. In this study now, three doubly fluorescent-labeled analogs of neuropeptide Y have been synthesized that still bind to the Y(5) receptor with high affinity to investigate the conformation in solution and, for the first time, to probe the conformational changes upon binding of the ligand to its receptor in cell membrane preparations. The results obtained from the fluorescence resonance energy transfer investigations clearly show considerable differences in transfer efficiency that depend both on the solvent as well as on the peptide concentration. However, the studies do not support a pancreatic polypeptide-like folding of neuropeptide Y in the presence of membranes that express the human Y(5) receptor subtype.

Enantioselective artificial metalloenzymes based on a bovine pancreatic polypeptide scaffold

Site creation: Enantioselective artificial metalloenzymes have been created by grafting a new active site onto bovine pancreatic polypeptide through the introduction of an amino acid capable of coordinating a copper(II) ion. This hybrid catalyst gave good enantioselectivities in the Diels-Alder and Michael addition reactions in water (see scheme) and displayed a very high substrate selectivity.

Recognition of Neurohormones of the NPY Family by Their Receptors

In this review a structural approach developed to answer the question whether hormones from the neuropeptide Y (NPY) family are recognized directly from solution or from the membrane-bound state is described. The chosen strategy is built onto a comparison of a set of peptides with well-known pharmacology and investigates whether similarities of structures of pharmacologically related peptides are higher in solution or in the membrane-bound state. Moreover, we have established the membrane-association mode of these peptides and contributed to our understanding of the structural features of these hormones both when placed in bulk solution and when bound to membranes. As a result we propose a receptor recognition pathway that includes initial association with the membrane and requires the peptides to come off the membrane to diffuse into the binding pocket of the receptor. This review also presents methodology recently developed by us to simulate the structural transition the peptides undergo when diffusing from bulk solution onto the membrane.

Pyridone dipeptide backbone scan to elucidate structural properties of a flexible peptide segment

Whereas the C-terminal fragment of neuropeptide Y (NPY) has been structurally well-defined both in solution and as membrane-bound, detailed structural information regarding the proline-rich N-terminus is still missing. The systematic variation of each position by a conformationally constrained pyridone dipeptide building block within the amino terminal segment of NPY leads to a systematic receptor subtype selectivity of the neuropeptide. Thereby, the systematic dipeptide scan proved superior to the traditional L-Ala scan because it showed how to modify the N-terminus in order to obtain increasingly more Y1 or Y5 receptor selective ligands. NMR and CD spectroscopic analyses were used to characterize the stepwise rigidification of the N-terminus of NPY when up to three dipeptide building blocks were incorporated by solid-phase peptide synthesis. The pyridone dipeptide increases the hydrophobicity of the amino terminus of NPY, and this allows the tuning of the membrane affinity of NPY. The amphiphilic C-terminal helix of 3-fold-substituted NPY thus becomes visible by selective line broadening in the (1)H NMR. Accordingly, we could structurally characterize protein segments that are too flexible for other methods.

Recognition of GPCRs by Peptide Ligands and Membrane Compartments theory: Structural Studies of Endogenous Peptide Hormones in Membrane Environment

One of the largest family of cell surface proteins, G-protein coupled receptors (GPCRs) regulate virtually all known physiological processes in mammals. With seven transmembrane segments, they respond to diverse range of extracellular stimuli and represent a major class of drug targets. Peptidergic GPCRs use endogenous peptides as ligands. To understand the mechanism of GPCR activation and rational drug design, knowledge of three-dimensional structure of receptor–ligand complex is important. The endogenous peptide hormones are often short, flexible and completely disordered in aqueous solution. According to “Membrane Compartments Theory”, the flexible peptide binds to the membrane in the first step before it recognizes its receptor and the membrane-induced conformation is postulated to bind to the receptor in the second step. Structures of several peptide hormones have been determined in membrane-mimetic medium. In these studies, micelles, reverse micelles and bicelles have been used to mimic the cell membrane environment. Recently, conformations of two peptide hormones have also been studied in receptor-bound form. Membrane environment induces stable secondary structures in flexible peptide ligands and membrane-induced peptide structures have been correlated with their bioactivity. Results of site-directed mutagenesis, spectroscopy and other experimental studies along with the conformations determined in membrane medium have been used to interpret the role of individual residues in the peptide ligand. Structural differences of membrane-bound peptides that belong to the same family but differ in selectivity are likely to explain the mechanism of receptor selectivity and specificity of the ligands. Knowledge of peptide 3D structures in membrane environment has potential applications in rational drug design.

Re-evaluation of receptor-ligand interactions of the human neuropeptide Y receptor Y1: a site-directed mutagenesis study

Interactions of the human NPY (neuropeptide Y) receptor Y1 with the two endogenous agonists NPY and peptide YY and two non-peptide antagonists were investigated using site-directed mutagenesis at 17 positions. The present study was triggered by contradictions among previously published reports and conclusions that seemed inconsistent with sequence comparisons across species and receptor subtypes. Our results show that Asp287, at the border between TM (transmembrane) region 6 and EL3 (extracellular loop 3) influences peptide binding, while two aspartic residues in EL2 do not, in agreement with some previous studies but in disagreement with others. A hydrophobic pocket of the Y1 receptor consisting of Tyr100 (TM2), Phe286 (TM6) and His298 (EL3) has been proposed to interact with the amidated C-terminus of NPY, a theory that is unsupported by sequence comparisons between Y1, Y2 and Y5. Nevertheless, our results confirm that these amino acid residues are critical for peptide binding, but probably interact with NPY differently than proposed previously. Studies with the Y1-selective antagonist SR120819A identified a new site of interaction at Asn116 in TM3. Position Phe173 in TM4 is also important for binding of this antagonist. In contrast with previous reports, we found that Phe173 is not crucial for the binding of BIBP3226, another selective Y1 receptor antagonist. Also, we found that position Thr212 (TM5) is important for binding of both antagonists. Our mutagenesis results and our three-dimensional model of the receptor based on the high-resolution structure of bovine rhodopsin suggest new interactions for agonist as well as antagonist binding to the Y1 receptor.

Are Hormones from the Neuropeptide Y Family Recognized by Their Receptors from the Membrane-Bound State?

Hormones and many other neurotransmitters, growth factors, odorant molecules, and light all present stimuli for a class of membrane-anchored receptors called G protein-coupled receptors (GPCRs). The GPCRs are the largest family of cell-surface receptors involved in signal transduction. About 1% of all known genes of Drosophila and more than 5% of the genes of Caenorhabditis elegans encode GPCRs. In addition, more than 50% of current therapeutic agents on the market target these receptors. When the enormous biological and pharmaceutical importance of these receptors is considered, it is surprising how little is known about the mechanism with which these receptors recognize their natural ligands. In this review we present a structural approach, utilizing techniques of high-resolution NMR spectroscopy, to address the question of whether peptides from the neuropeptide Y family of neurohormones are recognized directly from solution or from the membrane-bound state. In our studies we discovered that the structures of the membrane-bound species are better correlated to the pharmacological properties of these peptides than the solution structures are. These findings are supported by the observation that many biophysical properties of these peptides seem to be optimized for membrane binding. We finally present a scenario of possible events during receptor recognition.

Backbone Thioester Exchange: A New Approach to Evaluating Higher Order Structural Stability in Polypeptides

An amide bond has been replaced by a thioester in bovine pancreatic polypeptide (bPP) to allow rapid and reversible (dynamic) exchange of the alpha-helical segment with other thiols in solution. We have begun to study the higher order structural stability of bPP by measuring the equilibrium constant of the "backbone thioester exchange" (BTE) reaction. The extent to which the equilibrium (KBTE) favors one set of peptides over the other, which can be easily measured, can be directly correlated to the energy gained from favorable noncovalent interactions that occur between peptide segments on either side of the thioester bond (Kfold).

Structural Similarities of Micelle-bound Peptide YY (PYY) and Neuropeptide Y (NPY) are Related to their Affinity Profiles at the Y Receptors

Here, we investigate the structure of porcine peptide YY (pPYY) both when unligated in solution at pH 4.2 and when bound to dodecylphosphocholine (DPC) micelles at pH 5.5. pPYY in solution displays the PP-fold, with the N-terminal segment being back-folded onto the C-terminal alpha-helix, which extends from residue 17 to 31. In contrast to the solution structure of Keire et al. published in the year 2000 the C-terminal helix does not display a kink around residue 23-25. The root mean square deviation (RMSD) for backbone atoms of the NMR ensemble of conformers to the mean structure is 0.99(+/-0.35) Angstrom for residues 14-31. The back-fold is supported by values of 0.60+/-0.1 for the (15)N(1)H-NOE and by generalized order parameters S(2) of 0.74+/-0.1 for residues 5-31 which indicate that the peptide is folded in that segment. We have additionally used DPC micelles as a membrane model and determined the structure of pPYY when bound to it. Therein, an alpha-helix occurs in the segment comprising residues 17-31 and the N terminus freely diffuses in solution. The hydrophobic side of the amphipathic helix forms the micelle-binding interface and hydrophobic side-chains extend into the micelle interior. A significant stabilization of helical conformation occurs in the C-terminal pentapeptide, which is important for receptor binding. The latter is supported by positive values of the heteronuclear NOE in that segment (0.52+/-0.1 compared to 0.08+/-0.4 for the unligated form) and by values of S(2) of 0.6+/-0.2 (versus 0.38+/-0.2 for the unligated form). The structures of micelle-bound pPYY and pNPY are much more similar than those of pPYY and bPP with pairwise RMSDs of 1.23(+/-0.21)A or 3.21(+/-0.39) Angstrom, respectively. In contrast to the conformational similarities in the DPC-bound state their structures in solution are very different. In fact pPYY is more similar to bPP, which with its strong preference for the Y(4) receptor displays a completely different binding profile. Considering the high degree of sequence homology of pNPY and pPYY (>80%) and the fact, that their binding affinities at all receptor subtypes are high and, more importantly, rather similar, it is much more likely that PYY and NPY are recognized by the Y receptors from the membrane-bound state. As a consequence of the latter the PP-fold is not important for recognition of PYY or NPY at the Y receptors. To our knowledge this work provides for the first time strong arguments derived from structural data that support a membrane-bound receptor recognition pathway.

Analysis of pancreatic polypeptide cDNA from the house musk shrew, Suncus murinus, suggests a phylogenetically closer relationship with humans than for other small laboratory animal species

Pancreatic polypeptide was isolated and sequenced from endocrine cells of the pancreas from an insectivore, the house musk shrew, Suncus murinus. The primary sequence was APLEPAYPGD(10)NATPEQMAQY(20)AAELRKYINM(30)VTRPRYamide. This is the first polypeptide hormone to be characterised from this species and is typical of the primary sequences of pancreatic polypeptide of other animals, being a C-terminal-amidated peptide with 36 residues. Comparison with several vertebrate sequences shows that it has more in common with the human form than do the forms from common laboratory animals such as rabbits, rats, mice and guinea-pigs.

Bovine pancreatic polypeptide (bPP) undergoes significant changes in conformation and dynamics upon binding to DPC micelles

The pancreatic polypeptide (PP), a 36-residue, C-terminally amidated polypeptide hormone is a member of the neuropeptide Y (NPY) family. Here, we have studied the structure and dynamics of bovine pancreatic polypeptide (bPP) when bound to DPC-micelles as a membrane-mimicking model as well as the dynamics of bPP in solution. The comparison of structure and dynamics of bPP in both states reveals remarkable differences. The overall correlation time of 5.08ns derived from the 15N relaxation data proves unambiguously that bPP in solution exists as a dimer. Therein, intermolecular as well as intramolecular hydrophobic interactions from residues of both the amphiphilic helix and of the back-folded N terminus contribute to the stability of the PP fold. The overall rigidity is well-reflected in positive values for the heteronuclear NOE for residues 4-34. The membrane-bound species displays a partitioning into a more flexible N-terminal region and a well-defined alpha-helical region comprising residues 17-31. The average RMSD value for residues 17-31 is 0.22(+/-0.09)A. The flexibility of the N terminus is compatible with negative values of the heteronuclear NOE observed for the N-terminal residues 4-12 and low values of the generalized order parameter S(2). The membrane-peptide interface was investigated by micelle-integrating spin-labels and H,2H exchange measurements. It is formed by those residues which make contacts between the C-terminal alpha-helix and the polyproline helix. In contrast to pNPY, also residues from the N terminus display spatial proximity to the membrane interface. Furthermore, the orientation of the C terminus, that presumably contains residues involved in receptor binding, is different in the two environments. We speculate that this pre-positioning of residues could be an important requirement for receptor activation. Moreover, we doubt that the PP fold is of functional relevance for binding at the Y(4) receptor.

Distinguishing native conformations of proteins from decoys with an effective free energy estimator based on the OPLS all-atom force field and the Surface Generalized Born solvent model

Protein decoy data sets provide a benchmark for testing scoring functions designed for fold recognition and protein homology modeling problems. It is commonly believed that statistical potentials based on reduced atomic models are better able to discriminate native-like from misfolded decoys than scoring functions based on more detailed molecular mechanics models. Recent benchmark tests on small data sets, however, suggest otherwise. In this work, we report the results of extensive decoy detection tests using an effective free energy function based on the OPLS all-atom (OPLS-AA) force field and the Surface Generalized Born (SGB) model for the solvent electrostatic effects. The OPLS-AA/SGB effective free energy is used as a scoring function to detect native protein folds among a total of 48,832 decoys for 32 different proteins from Park and Levitt's 4-state-reduced, Levitt's local-minima, Baker's ROSETTA all-atom, and Skolnick's decoy sets. Solvent electrostatic effects are included through the Surface Generalized Born (SGB) model. All structures are locally minimized without restraints. From an analysis of the individual energy components of the OPLS-AA/SGB energy function for the native and the best-ranked decoy, it is determined that a balance of the terms of the potential is responsible for the minimized energies that most successfully distinguish the native from the misfolded conformations. Different combinations of individual energy terms provide less discrimination than the total energy. The results are consistent with observations that all-atom molecular potentials coupled with intermediate level solvent dielectric models are competitive with knowledge-based potentials for decoy detection and protein modeling problems such as fold recognition and homology modeling.

The NMR-derived conformation of neuropeptide F from Moniezia expansa

The solution structure of neuropeptide F (NPF), from the flatworm (platyhelminthes) Moniezia expansa, has been determined by (1)H NMR spectroscopy at 800 MHz in 60%/40% CD(3)OH/H(2)O. NPF is the most abundant neuropeptide in platyhelminthes. The secondary structure of NPF contains an alpha helix from residues Lys(14) to Ile(31), while the N terminus, consisting of residues Pro(-2) to Asn(13), and the C-terminus, consisting of residues Gly(32) to Phe(36), are in a random conformation. The structure was calculated by a simulated annealing protocol, and the conformational data are compared to the porcine neuropeptide Y (NPY), a peptide hormone and neurotransmitter. The exact function of NPF is unknown, but structural similarity with porcine NPY indicates that its mode of action is similar. These structural data can serve as a starting point in the design of new antiparasitic drugs.

Structure and receptor binding of PYY analogs

Differences in the structure of PYY and two important analogs, PYY [3-36] and [Pro34]PYY, are evaluated. Y-receptor subtype ligand binding data are used in conjunction with structural data to develop a model for receptor subtype selective agonists. For PYY it is proposed that potent binding to Y1, Y4 and Y5 receptors requires the juxtaposition of the two termini while Y2 binding only requires the C-terminal helix. Further experiments that delineate between primary and tertiary structure contributions for receptor binding and activation are required to support the hypothesis that tertiary structure is stable enough to influence the expression of PYY's bioactivity.

Pancreatic polypeptide receptors: affinity, sodium sensitivity and stability of agonist binding

Cloned rat, human and guinea-pig Y4 pancreatic polypeptide (PP) receptors expressed in Chinese hamster ovary (CHO) cells, as well as the rabbit Y4-like PP receptor, show a selective sensitivity to Na+ over K+ ion in PP attachment, but little sensitivity to Na+ in dissociation of bound PP peptides. Agonist binding to Y4 receptors of intact CHO cells also shows much greater sensitivity to Na+ over K+, and a tenacious attachment of the bound agonist. Binding sensitivity to K+ is greatly enhanced upon receptor solubilization. Pancreatic polypeptide sites also show large sensitivity to modulators of Na+ transport such as N5-substituted amilorides and to RFamides, as different from Y1 or Y2 receptors. Thus, PP binding is modulated by cation-induced changes in site environment (with selectivity for Na+) and ultimately results in a blocking attachment. This would support receptor operation in the presence of ion gradients, as well as prolonged agonist-delimited signaling activity (which can include partial antagonism). Also, this could point to an evolutionary adaptation enabling small numbers of PP receptors to perform extensive metabolic tasks in response to low agonist signals.

Design, synthesis and characterisation of a peptide with oxaloacetate decarboxylase activity

The synthesis of Oxaldie-3, a synthetic 31-residue peptide with oxaloacetate decarboxylase activity, is described. Biophysical characterisation by gel filtration, CD and NMR spectroscopy indicated that the peptide adopted a folded structure in solution. Oxaldie-3 was an efficient catalyst at concentrations as low as 2 microM, 100-fold lower than the previously described Oxaldie-2, which relied on aggregating alpha-helices for activity. Oxaldie-3 speeded decarboxylation by more than three orders of magnitude relative to simple amines.

Ab initio protein structure prediction using physicochemical potentials and a simplified off-lattice model

This study describes a computational method for ab inito protein structure prediction. Protein conformation has been modeled by using six optimized backbone torsion angles and fixed side chains approximating rotationally averaged real side chains. The approximations aim to keep complexity of the structure description to a minimum without seriously compromising the accuracy of the structural representation. An evolutionary Monte Carlo algorithm has been developed to search through this restricted conformational space to locate low-energy protein structures. A simple physicochemical force field has been developed to assess the energies of different conformations within this structural description. The corresponding residue interaction energies are based on hydrophobic, hydrophilic, steric, and hydrogen-bonding potentials. The search procedure has been used to locate native energy minima from primary sequence alone. The 3-D structures of polypeptides up to 38 residues with both beta and alpha secondary structural elements have been accurately predicted. The search procedure has been found to be highly efficient and follows an energetically and structurally plausible pathway to locate native populations. The simple force field described in the study has been compared with a more complex all-atom model and been found to be similarly effective in predicting the structures of proposed independent folding units. Proteins 2001;43:186-202.

Primary structures of PYY, [Pro(34)]PYY, and PYY-(3-36) confer different conformations and receptor selectivity

We synthesized PYY-(1-36) (nonselective between Y(1) and Y(2) receptor subtype agonists), [Pro(34)]PYY (selective for Y(1)), and PYY-(3-36) (selective for Y(2)) to determine whether solution conformation plays a role in receptor subtype selectivity. The three peptides exhibited the expected specificities in displacing labeled PYY-(1-36) from cells transfected with Y(1) receptors (dissociation constants = 0.42, 0.21, and 1,050 nM, respectively) and from cells transfected with Y(2) receptors (dissociation constants = 0.03, 710, and 0.11 nM, respectively) for PYY-(1-36), [Pro(34)]PYY, and PYY-(3-36). Sedimentation equilibrium analyses revealed that the three PYY analogs were 80-90% monomer at the concentrations used for the subsequent circular dichroism (CD) and (1)H-nuclear magnetic resonance (NMR) studies. CD analysis measured helicities for PYY-(1-36), [Pro(34)]PYY, and PYY-(3-36) of 42%, 31%, and 24%, suggesting distinct differences in secondary structure. The backbone (1)H-NMR resonances of the three peptides further substantiated marked conformational differences. These patterns support the hypothesis that Y(1) and Y(2) receptor subtype binding affinities depend on the secondary and tertiary solution state structures of PYY and its analogs.

2-36[K4,RYYSA(19-23)]PP a novel Y5-receptor preferring ligand with strong stimulatory effect on food intake

Members of the neuropeptide Y (NPY) family regulate many physiological processes via interaction with at least four functional, pharmacologically distinct Y-receptors. However, selective antagonists developed for several subtypes have not been useful in defining particular Y-receptor functions in vivo. To identify critical residues within members of the NPY family required for Y-receptor subtype-selectivity we have determined the contribution of each residue within NPY to receptor binding by replacing them with L-alanine. In a second study, chimeric peptides where single or stretches of residues were interchanged between members of the NPY family were generated and tested in radioligand binding studies. Overall, substituted alanine analogues exhibited similar orders of affinities at each Y-receptor subtype with no obvious subtype-selectivity. Residues of particular interest are Leu30 which exhibited selectivity for the Y4-receptor, whereas Asp16 does not appear to play any role in ligand binding. Several chimeric peptides, e.g., [K4]pancreatic polypeptide ([K4]PP) and [RYYSA(19-23)]PP clearly showed higher affinity at the Y4 and Y5 subtypes compared to the Y1 and Y2 subtypes. In addition, the transfer of a proline residue from position 14 to 13 in peptide YY decreases its affinity at the Y1-, Y4- and Y5-receptors but is unchanged at the Y2 subtype. Combining these results, and with the help of molecular modelling, second generation chimeras were designed. The most significant improvement was achieved in chimera 2-36[K4,RYYSA(19-23)]PP where the affinity for the Y5 subtype increased by ninefold over that from NPY. Several of these compounds were also tested for their ability to stimulate food intake in a rat model. Interestingly, again 2-36[K4,RYYSA(19-23)]PP showed the most dramatic effect with a major increase on food intake over a range of doses compared to NPY suggesting a possible synergistic effect of several Y-receptors on feeding behaviour.

Molecular dynamics of NPY Y1 receptor activation

A three-dimensional model of the human neuropeptide Y(NPY)Y1 receptor (hY1) 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-N2-(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 hY1 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.

Neuropeptide families and their receptors: evolutionary perspectives

Examination of families of neuropeptides and their receptors can provide information about phyletic relationships and evolutionary processes. Within an individual a given signal molecule may serve many diverse functions, mediated via subtypes of the receptor which may be coupled to their transduction mechanisms in different ways. The rate of evolution of a peptide may reflect or be reflected in the rate of evolution of its receptor. For example, in the neuropeptide Y (NPY) family, pancreatic polypeptide (PP) shows significant structural diversity, while NPY is highly conserved. Molecular forms of a given subtype of NPY receptor that is selectively activated by NPY (Y1 or Y2 or Y5) are also highly conserved, but the subtype that is primarily activated by PP (Y4), shows remarkable diversity. Also, between receptor subtypes there can be remarkable diversity. This is evident in several neuropeptide families, where a neuropeptide sequence is highly conserved across a wide range of species but where the receptor homology of subtypes with species tends to be much lower than homology between species. For example, human and rat vasopressin are identical, but the human V(1)- or V(2)-vasopressin receptors are approximately 80% homologous with rat V(1)- or V(2)-receptors, but within humans or rats the V(1)-receptor is less than 50% homologous with the V(2)-receptor. Furthermore, duplication of an ancestral gene is thought to have led to the co-presence in eutherian mammals of oxytocin and vasopressin, which have maintained a close structural similarity, yet in many species the oxytocin receptor is only 30 to 50% homologous with vasopressin receptors. Thus it appears that there has been greater evolutionary pressure to conserve the signal molecule, than to conserve the structure of the receptor. Evaluation of the evolution of neuropeptides and their receptors may be useful in determining phyletic relationships. Traditional classification places the guinea pig as a hystricomorph rodent within the same order (Rodentia) as the muriform or myomorph rat and mouse. However, molecular analyses of polypeptides have led to the suggestion that guinea pigs belong to a distinct order. Analysis of several neuropeptide sequences and the Y4 receptor supports this view. In general terms for both neuropeptides and receptors, sequence homology reflects phylogeny and taxonomy as based on morphological features. Within the oxytocin/vasopressin family in which peptides and receptors have been characterised in invertebrate representatives as well as fish and amphibia in addition to mammals, the molecular diversity correlates well with evolutionary diversity.

Aspects of the molecular structure and dynamics of neuropeptide Y

Human neuropeptide Y (hNPY) and the Q34-->P34 mutant (P34-hNPY) have been characterized by CD spectroscopy. hNPY self-associates in aqueous solution with a dimerization constant in the micromolar range. The self-association correlates with an increase in secondary-structure content which was studied as a function of concentration, temperature and pH. The effects of temperature were measured in water (5-84 degrees C) and in ethanediol/water (2 : 1) (-90 degrees to +90 degrees C). A single-residue mutation, Q34-->P34, affects the pH, thermal and self-association properties of NPY. The CD results are correlated with photochemically induced dynamic nuclear polarization NMR experiments which show that the tyrosines at the interface between two monomer units present limited accessibility to a photoreactive dye. An equilibrium state is described, involving a PP-fold monomer form and a handshake dimer form, that accommodates the physicochemical properties of NPY.

A descriptive analysis of populations of three-dimensional structures calculated from primary sequences of proteins by OSIRIS

Among different ab initio approaches to calculate 3D-structures of proteins out of primary sequences, a few are using restricted dihedral spaces and empirical equations of energy as is OSIRIS. All those approaches were calibrated on a few proteins or fragments of proteins. To optimize the calculation over a larger diversity of structures, we need first to define for each sequence what are good conditions of calculations in order to choose a consensus procedure fitting most 3D-structures best. This requires objective classification of calculated 3D-structures. In this work, populations of avian and bovine pancreatic polypeptides (APP, BPP) and of calcium-binding protein (CaBP) are obtained by varying the rate of the angular dynamics of the second step of OSIRIS. Then, 3D-structures are clustered using a nonhierarchical method, SICLA, using rmsd as a distance parameter. A good clustering was obtained for four subpopulations of APP, BPP and CaBP. Each subpopulation was characterized by its barycenter, relative frequency and dispersion. For the three alpha-helix proteins, after the step 1 of OSIRIS, most secondary structures were correct but molecules have a few atomic contacts. Step 2, i.e., the angular dynamics, resolves those atomic contacts and clustering demonstrates that it generates subpopulations of topological conformers as the barycenter topologies show.

Prediction of the three-dimensional structure of proteins using the electrostatic screening model and hierarchic condensation

We describe a method for predicting the three-dimensional (3-D) structure of proteins from their sequence alone. The method is based on the electrostatic screening model for the stability of the protein main-chain conformation. The free energy of a protein as a function of its conformation is obtained from the potentials of mean force analysis of high-resolution x-ray protein structures. The free energy function is simple and contains only 44 fitted coefficients. The minimization of the free energy is performed by the torsion space Monte Carlo procedure using the concept of hierarchic condensation. The Monte Carlo minimization procedure is applied to predict the secondary, super-secondary, and native 3-D structures of 12 proteins with 28-110 amino acids. The 3-D structures of the majority of local secondary and super-secondary structures are predicted accurately. This result suggests that control in forming the native-like local structure is distributed along the entire protein sequence. The native 3-D structure is predicted correctly for 3 of 12 proteins composed mainly from the alpha-helices. The method fails to predict the native 3-D structure of proteins with a predominantly beta secondary structure. We suggest that the hierarchic condensation is not an appropriate procedure for simulating the folding of proteins made up primarily from beta-strands. The method has been proved accurate in predicting the local secondary and super-secondary structures in the blind ab initio 3-D prediction experiment.

Conformation of Vespa basalis mastoparan-B in trifluoroethanol-containing aqueous solution

Mastoparan-B, a tetradecapeptide isolated from the venom of the hornet Vespa basalis, belongs to the mastoparan analogs of vespid venom with the lysine residues common for all mastoparan family toxins at positions 4, 11 and 12. Here we use 1H-NMR spectroscopy and hybrid distance geometry-simulated annealing calculation to investigate its three-dimensional structure in trifluoroethanol-containing aqueous solution. The calculated structure shows that residues 3-14 adopt an amphiphilic alpha-helical structure in which the residues with hydrophilic side chains (i.e. Lys-4, Ser-5, Ser-8, Lys-11, Lys-12) are located on one side and the residues with hydrophobic side chains (i.e. Leu-3, Ile-6, Trp-9, Ala-10, Val-13, Leu-14) located on the other side of the molecule. The overall structural features a very similar to the conformation of mastoparan-X reconstituted in vesicle [Wakamatsu et al. (1992) Biochemistry 31, 5654-5660] in spite of the substitutions made for eight residues with distinctly different hydrophobicity. These substitutions lead to a larger hydrophobic moment for the alpha-helical segment and further mobilized N-terminal. This study will help reveal the conformational significance of mastoparan toxins with respect to their potency and activity in G protein regulation.

Simulating the folding of small proteins by use of the local minimum energy and the free solvation energy yields native-like structures

Assuming that the protein primary sequence contains all information required to fold a protein into its native tertiary structure, we propose a new computational approach to protein folding by distributing the total energy of the macromolecular system along the torsional axes. We further derive a new semiempirical equation to calculate the total energy of a macromolecular system including its free energy of solvation. The energy of solvation makes an important contribution to the stability of biological structures. The segregation of hydrophilic and hydrophobic domains is essential for the formation of micelles, lipid bilayers, and biological membranes, and it is also important for protein folding. The free energy of solvation consists of two components: one derived from interactions between the atoms of the protein, and the second resulting from interactions between the protein and the solvent. The latter component is expressed as a function of the fractional area of protein atoms accessible to the solvent. The protein-folding procedure described in this article consists of two successive steps: a theoretical transition from an ideal alpha helix to an ideal beta sheet is first imposed on the protein conformation, in order to calculate an initial secondary structure. The most stable secondary structure is built from a combination of the lowest energy structures calculated for each amino acid during this transition. An angular molecular dynamics step is then applied to this secondary structure. In this computational step, the total energy of the system consisting of the sum of the torsional energy, the van der Waals energy, the electrostatic energy, and the solvation energy is minimized. This process yields 3-D structures of minimal total energy that are considered to be the most probable native-like structures for the protein. This method therefore requires no prior hypothesis about either the secondary or the tertiary structure of the protein and restricts the input of data to its sequence. The validity of the results is tested by comparing the crystalline and computed structures of four proteins, i.e., the avian and bovine pancreatic polypeptide (36 residues each), uteroglobin (70 residues), and the calcium-binding protein (75 residues); the C alpha-C alpha maps show significant homologies and the position of secondary structure domains; that of the alpha helices is particularly close.

Structure of the presynaptic (Y2) receptor-specific neuropeptide Y analog ANA-NPY

Neuropeptide Y analog ANA-NPY or [Leu-17, Gln-19, Ala-20, Ala-23, Leu-28, Leu-31]NPY(13-36)-amide binds to postjunctional or Y1 receptors to raise blood pressure and to prejunctional or Y2 receptors to inhibit neurotransmitter release. ANA-NPY affects Y2 receptors in the same way as intact NPY but exhibits far less potent effects on Y1 receptors. The structure of ANA-NPY was examined using two-dimensional proton nuclear magnetic resonance spectroscopy. Complete assignment of all backbone and side chain hydrogens was accomplished with totally correlated spectroscopy (TOCSY) experiments providing through-bond 1H-1H connectivities, and nuclear Overhauser effect spectroscopy (NOESY), providing the through-space and sequential backbone connectivities. The tertiary solution structure of the peptide was performed using distance geometry and dynamic simulated annealing. ANA-NPY exhibits a helical structure with strong amphipathic character with a bend around Glu-24 indicating that the C-terminal segment 25-35 forms a single alpha-helical motif.

Conformational properties of the proline region of porcine neuropeptide Y by CD and 1H-nmr spectroscopy

We synthesized porcine neuropeptide Y (pNPY) N-terminal fragments by solid-phase synthesis techniques and analyzed them for solution conformational properties by CD and 1H-nmr spectroscopy. The analogues pNPY1-9 and pNPY1-14 displayed CD spectra indicative of random structures and showed no evidence for induced alpha-helical structures in trifluoroethanol (TFE) up to 50%. However, the CD spectra of pNPY1-9 suggested a conformational shift in tetrahydrofuran. Although in aqueous solution the CD spectra of pNPY1-21 indicated random structures with induction of only a small percentage of alpha-helix in aqueous TFE, pNPY1-25 displayed 13% alpha-helical structure in aqueous solution that increased to 40 and 41% by the addition of TFE and methanol, respectively. The nmr spectra of pNPY1-9 and the proline region of pNPY1-25 indicated extended structures with all-trans conformers at Pro5 and Pro8 for pNPY1-9 and at Pro5, Pro8, and Pro13 for pNPY1-25; in each case the Tyr1-Pro2 amide bond was in both cis and trans conformations. However, observed nuclear Overhauser effect correlations and HN exchange experiments indicated an alpha-helical segment in pNPY1-25 initiated by Pro13 and extending from residues 14 to 25. Thus, the N-terminal polyproline region of NPY has no propensity to fold into a regular secondary structure, although Pro13 is a helix initiator, a result consistent with the proposed role of this amino acid in the NPY structural model.

Structural characterizations of neuropeptide tyrosine (NPY) and its agonist analog [Ahx5-17]NPY by NMR and molecular modeling

The structures of human NPY and of its centrally truncated agonist analog [Ahx5-17]NPY have been investigated in DMSO-d6 by two-dimensional NMR and by molecular modeling. For both peptides, a complete resonance assignment was achieved and a large number (more than 200) of inter-residue NOE connectivities were observed, including long-range connectivities between the N- and C-terminal ends of the chain. Molecular models were calculated using NOE constraints by distance geometry, simulated annealing and conjugate gradient energy minimization. The results indicate that both peptides are folded in the center of their chain, NPY adopting the hairpin shape, whereas the central portion of [Ahx5-17]NPY is characterized by relatively large loops. In contrast to previous models, practically no alpha-helical structure exists for these peptides under our conditions, but two beta-turns are found in NPY and one in [Ahx5-17]NPY. The proximity of the terminal ends could be the determinant factor for their activity.