Analysis of tachykinin-binding site interactions using NMR and energy calculation data of potent cyclic analogues of substance P

Solution structure of amphibian tachykinin Uperolein bound to DPC micelles

Uperolein, a physalaemin-like endecapeptide, has been shown to be selective for Neurokinin 1 receptor. As a first step towards understanding the structure-activity relationship, we report the membrane-induced structure of Uperolein with the aid of circular dichroism and 2D (1)H NMR spectroscopy. Sequence-specific resonance assignments of protons have been made using correlation spectroscopy (TOCSY, DQF-COSY) and NOESY spectroscopy. The interproton distance constraints and dihedral angle constraints have been utilized to generate a family of structures using torsion angle molecular dynamics within program DYANA. The conformational range of the peptide revealed by NMR and CD studies has been analysed in terms of characteristic secondary features. Analysis of NMR data indicates that the global fold of Uperolein can be explained in terms of equilibrium between 3(10)-helix and alpha-helix from residues 5 to 11. An extended highly flexible N-terminus displays some degree of order and a possible turn structure. A comparison between the structures of Uperolein and Substance P, a prototype and endogenous Neurokinin 1 receptor agonist, indicates several common features in the distribution of hydrophobic and hydrophilic residues. Both the peptides show an amphiphilic character towards the middle region. The similarities suggest that the molecules interact with the receptor in an analogous manner.

Conformational analysis of the C-terminal Gly-Leu-Met-NH2 tripeptide of substance P bound to the NK-1 receptor

We examined the effect of simultaneously incorporating proline or proline-amino acid chimeras in positions 9, 10, and/or 11 of substance P, on the affinity for the two NK-1 binding sites and on second-messenger activation. Because these 3-substituted prolines constrain not only the (phi,psi) values of the peptide backbone, but also the chi space of the amino acid side chain, we were able to gather data on the structural requirements for high-affinity binding to the NK-1 receptor. We were able to confirm that this C-terminal component is crucial and that it should adopt an extended conformation close to a polyproline II structure when bound to the receptor. The partial additivity of these constraints, more specifically, for the NK-1M site, suggests that the peptide backbone flexibility around the hinge-point residue Gly9 is essential to subtly position crucial side chains.

Membrane-induced structure of the mammalian tachykinin neuropeptide gamma

Neuropeptide gamma (NPgamma) is a neurokinin-2 (NK-2) receptor selective agonist, which plays an important role in mediation of asthma and elicits a wide range of biological responses like bronchoconstriction, vasodepression and regulation of endocrine functions. The structure determination of this peptide agonist is important in understanding the molecular basis of peptide ligand recognition by the receptor and for rational drug design. In the present study we report the solution structure of NPgamma characterized by circular dichroism (CD) spectropolarimetry and 2D (1)H NMR spectroscopy in both aqueous and membrane mimetic solvents. Effect of calcium ions on the conformation of NPgamma was also studied using CD spectropolarimetry. Sequence-specific resonance assignments of protons have been made with the aid of correlation spectroscopy experiments and nuclear Overhauser effect spectroscopy experiments. The distance constraints obtained from the NMR data have been utilized to generate a family of structures, which have been refined using restrained energy minimization and dynamics. These data show that in water NPgamma prefers to be in an extended chain conformation whereas a helical conformation is induced in the central core and the C-terminal region of the peptide (K13-M21) in the presence of perdeuterated dodecylphosphocholine micelles, a membrane model system. A type II' beta turn from H9 to R11 precedes the helical core in the C-terminus of NPgamma. N-terminus of NPgamma also displays some degree of order and a possible turn structure. Conformation adopted by NPgamma in presence of lipid micelles represents a structural motif typical of NK-2 selective agonists and is similar to that observed for Neurokinin A in hydrophobic environment. The observed conformational features have been correlated to the binding ability and biological activity of NPgamma.

Membrane-Induced Structure of Scyliorhinin I: A Dual NK1/NK2 Agonist

Scyliorhinin I, a linear decapeptide, is the only known tachykinin that shows high affinity for both NK-1 and NK-2 binding sites and low affinity for NK-3 binding sites. As a first step to understand the structure-activity relationship, we report the membrane-induced structure of scyliorhinin I with the aid of circular dichroism and 2D-(1)H NMR spectroscopy. Sequence specific resonance assignments of protons have been made from correlation spectroscopy (TOCSY, DQF-COSY) and NOESY spectroscopy. The interproton distance constraints and dihedral angle constraints have been utilized to generate a family of structures using DYANA. The superimposition of 20 final structures has been reported with backbone pairwise root mean-square deviation of 0.38 +/- 0.19 A. The results show that scyliorhinin I exists in a random coil state in aqueous environments, whereas helical conformation is induced toward the C-terminal region of the peptide (D4-M10) in the presence of dodecyl phosphocholine micelles. Analysis of NMR data is suggestive of the presence of a 3(10)-helix that is in equilibrium with an alpha-helix in this region from residue 4 to 10. An extended highly flexible N-terminus of scyliorhinin I displays some degree of order and a possible turn structure. Observed conformational features have been compared with respect to that of substance P and neurokinin A, which are endogenous agonists of NK-1 and NK-2 receptors, respectively.

Three-dimensional structure of the mammalian tachykinin peptide neurokinin A bound to lipid micelles

The solution structure of NKA, a decapeptide of mammalian origin, has been characterized by CD spectropolarimetry and 2D proton nuclear magnetic resonance (2D 1H-NMR) spectroscopy in both aqueous and membrane mimetic solvents. Unambiguous NMR assignments of protons have been made with the aid of correlation spectroscopy (DQF-COSY and TOCSY) experiments and nuclear Overhauser effect spectroscopy (NOESY and ROESY) experiments. The distance constraints obtained from the NMR data have been utilized to generate a family of structures, which have been refined using restrained energy minimization and dynamics. These data show that in water NKA prefers to be in an extended chain conformation whereas a helical conformation is induced in the central core and the C-terminal region (D4-M10) of the peptide in the presence of perdeuterated dodecylphosphocholine (DPC) micelles, a membrane model system. Though less defined the N-terminus also displays some degree of order and a possible turn structure. The conformation adopted by NKA in the presence of DPC micelles represents a structural motif typical of neurokinin-2 selective agonists and is similar to that reported for eledoisin in hydrophobic environment.

Solution structure of the tachykinin peptide eledoisin

Both the aqueous and the lipid-induced structure of eledoisin, an undecapeptide of mollusk origin, have been studied by two-dimensional proton nuclear magnetic resonance spectroscopy and distance geometry calculations. Unambiguous nuclear magnetic resonance assignments of protons have been made with the aid of correlation spectroscopy experiments and nuclear Overhauser effect spectroscopy experiments. The distance constraints obtained from the nuclear magnetic resonance data have been utilized in a distance geometry algorithm to generate a family of structures, which have been refined using restrained energy minimization and dynamics. These data show that, while in water and dimethyl sulfoxide, eledoisin prefers to be in an extended chain conformation, whereas in the presence of perdeuterated dodecylphosphocholine micelles, a membrane model system, helical conformation is induced in the central core and C-terminal region (K4-M11) of the peptide. N terminus, though less defined, also displays some degree of order and a possible turn structure. The conformation adopted by eledoisin in the presence of dodecylphosphocholine micelles is similar to the structural motif typical of neurokinin-2 selective agonists and with that reported for kassinin in hydrophobic environment.

Lipid induced conformation of the tachykinin peptide Kassinin

Both the aqueous and lipid-induced structure of Kassinin, a dodecapeptide of amphibian origin, has been studied by two-dimensional proton nuclear magnetic resonance (2D 1H-NMR) spectroscopy and distance geometry calculations. Unambiguous NMR assignments of protons have been made with the aid of correlation spectroscopy (DQF-COSY and TOCSY) experiments and nuclear Overhauser effect spectroscopy (NOESY and ROESY) experiments. The distance constraints obtained from the NMR data have been utilized in a distance geometry algorithm to generate a family of structures, which have been refined using restrained energy minimization and dynamics. These data show that, while in water Kassinin prefers to be in an extended chain conformation, in the presence of perdeuterated dodecylphosphocholine (DPC) micelles, a membrane model system, helical conformation is induced in the central core and C-terminal region (K4-M12) of the peptide. N-terminus though less defined also displays some degree of order and a possible turn structure. The conformation adopted by Kassinin in the presence of DPC micelles is consistent with the structural motif typical of neurokinin-1 selective agonists and with that reported for Eledoisin in hydrophobic environment.

NMR and molecular dynamics studies of tachykinins: conformation of substance P fragment 4-11

The conformation of the C-terminal octapeptide fragment of Substance P (SP4-11, Pro-Gln-Gln-Phe-Phe-Gly-Leu-Met-NH2) has been investigated by 2D-NMR and MD methods. The octapeptide exists in a blend of conformations. The molecule seems to shuttle between conformations with gamma-bends either at Phe5 or Gly6 or Gln3 or Leu7 and between a nearly extended structure.

Lipid-induced conformation of substance P

Both the aqueous and the lipid-induced structure of a representative and widely studied tachykinin, substance P, has been investigated by two-dimensional proton nuclear magnetic resonance (2D 1H-NMR) spectroscopy and distance geometry calculations. Unambiguous NMR assignments of protons have been made with the aid of correlation spectroscopy (COSY and TOCSY) experiments and Overhauser enhancement spectroscopy (ROESY and NOESY; experiments. The NMR data obtained were utilized in a distance geometry algorithm to generate a family of structures which were further refined using restrained energy minimization. These data show that, while in water substance P appears to favour an extended chain conformation, in the presence of perdeuterated dodecylphosphocholine (DPC) micelles as membrane model system an amphiphilic helical conformation is induced in the mid-region (Q5-Q8) of substance P. The conformation adopted by substance P in the presence of DPC micelles yields a structural motif typical of neurokinin-1 selective ligands, as proposed by Convert and coworkers (O. Convert et al., Neuropeptides 19, 259-270 (1991)).

Topographic analysis of the S7 binding subsite of the tachykinin neurokinin-1 receptor

Conformationally and configurationally restricted rotameric probes of phenylalanine have been incorporated in the sequence of substance P (SP)-Arg-Pro-Lys-Pro-Gln-Gln-Phe-Phe-Gly-Leu-Met-NH2-for analyzing the binding pockets of Phe7 (S7) and Phe8 (S8), in the neurokinin-1 receptor. These analogues of phenylalanine are (2S. 3R)- and (2S, 3S)-indanylglycines, E- and Z-alpha, beta-dehydrophenylalanines, and 2(S)-alpha, beta-cyclopropylphenylalanines [delta E Phe. delta Z Phe. inverted delta E2 (S) Phe, and inverted delta Z 2 (S) Phe]. Binding data obtained with either conformationally (Ing diastereoisomers) or configurationally (delta E Phe, delta Z Phe) probes have unveiled large differences in the binding potencies of these rotameric probes. With the support of nmr data and energy calculations done on these SP-substituted analogues, we attempt to answer questions inherent to such study. First, none of these six probes prevents the formation of bioactive conformation(s) of the backbone of SP. Second, both diastereoisomers (S, S) and (S, R) of indanylglycine preferentially adopt, in the sequence of SP, the gauche (-) and trans side-chain orientations, respectively, as previously postulated from energy calculations with model peptides. However, in solution, the difference in energy between these rotamers included in the sequence of SP, compared to model peptides, is small since the other rotamer can be detected in [(2S, 3R)Ing7]SP. Finally, from this study we can hypothesize that the large variations observed in the affinities of Phe7 substituted analogues of SP must come from steric hindrance in the S7 binding site, which drastically restricts the space filling around the C alpha-C beta bond of residue 7.

260 ps molecular dynamics simulation of substance P with hydrated dimyristoyl phosphatidyl choline bilayer

We present here results on 260 pico seconds (ps) molecular dynamics (MD) simulation of substance P (SP) in hydrated bilayer of dimyristoyl phosphatidyl choline (DMPC) (39 molecules of DMPC with 776 water molecules). 260 ps MD simulation has been carried out in 0.001 ps time interval with united atom force field, using AMBER 4.0 package. Non bonded pair list was updated every 20 cycles using 12.5 Angstrom cut off distance. Analysis of MD data is done using our package ANALMD. The obtained models are presented using graphics package RASMOL. All simulations, analysis of MD data and graphics is done on INDIGO-2, R-4400 extreme graphics work station. Our results show no systematic change in order parameter, but reduction in transfraction of the chain torsional angles, compared to our earlier results on MD simulation on hydrated DMPC bilayer without SP. C-terminal and central peptide residues adopt partial helical conformation. Helix type as classified on the basis of H-bonds is between alpha and 3(10). The peptide backbone shows flexibility during heating runs. Later, it is stabilized and there was not much change in the spatial position of the backbone. Lipid matrix serves the role of immobilization of the peptide backbone in a preferred conformation.

Computational molecular modeling in peptide drug design

The review concentrates on practical applications of computer molecular modeling in peptide drug design. The examples of the predictions (successful or not) made by computational modeling before synthesis of peptide analogs, not the explanations provided after synthesis and biological testing of peptides, are discussed. The review spans over 20 years of predictions made by computer molecular modeling for bradykinin, angiotensin, thyrotropin-releasing factor, tuftsin, substance P, CCK-related peptides, luliberin, alpha-melanotropin and opioid peptides. The described examples are discussed in terms of finding the optimal way to use computer modeling for peptide design. The step-by-step 'technology' of peptide design is outlined in detail.

NMR and molecular modeling investigations of the neuropeptide substance P in the presence of 15 mM sodium dodecyl sulfate micelles

To better understand the structural basis of the biological activity of the neuropeptide substance P SP; (Arg-Pro-Lys-Pro-Gln-Gln-Phe-Phe-Gly-Leu-Met-NH2), two-dimensional nmr spectroscopy experiments and simulated annealing calculations were used to investigate the conformation adopted in the presence of the membrane model system sodium dodecyl sulfate. It was determined that SP in the presence of SDS micelles undergoes a conformational equilibrium between an alpha- and a 3(10)-helix involving the midregion (Pro4-Gln5-Gln6-Phe7-Phe8) of the peptide. The C-terminus adopts an extended conformation while the N-terminus remains quite flexible. The conformation adopted by SP in the presence of SDS micelles yields a structure that is consistent with the model of a neurokinin-1 selective ligand proposed by Convert.

Binding of substance P to monolayers and vesicles made of phosphatidylcholine and/or phosphatidylserine

Analyses of interactions between substance P (SP) and phospholipids were performed by combined surface pressure and surface potential measurements in monolayers and by 13C-NMR experiments on liposomes. This study was carried out using synthetic SP molecules: [1-13C-Gly9]SP and [1-13C-Gly2]SP. Injection of SP into the aqueous subphase led to an expansion of phosphatidylcholine (PtdCho) or phosphatidylserine (PtdSer) monolayer surface area. An apparent association constant of SP for PtdSer was estimated to be around 10(6)-10(-7) M-1. The surface potential delta V/n varied linearly with the molecular area whereas the variation of surface pressure was biphasic, suggesting that at least two binding states contributed to the monolayer expansion. These two states Si (SP is inserted into the bilayer) and Ss (SP is stuck on the surface) were observed on vesicular membranes by 13C-NMR. The kinetic of interconversion between these two states can be estimated by NMR, the Ss state being the stablest one. No perpendicular insertion of SP into these vesicular preparations seemed to occur, as previously postulated. However, SP might form aggregates in contact with these model systems, leading to a loss of permeability of the lipid vesicles.

Influence of the replacement of amino acid by its D-enantiomer in the sequence of substance P. 2. Conformational analysis by NMR and energy calculations

The D-enantiomer of residues 2, 4, 5, 6, 7, 8, 10 and 11 was introduced in the sequence of Substance P: Arg-Pro-Lys-Pro-Gln-Gln-Phe-Phe-Gly-Leu-Met-NH1. The achiral glycine was replaced by a D-Ala residue. The conformations of the D-substituted analogues were analysed by NMR and molecular energy calculations. Introduction of a D-amino acid in the address sequence of SP (1 to 5) distorted the helical structure of [D-Pro2]SP and [D-Pro4]SP. A D-glutamine in position 5 hampered the formation of an helix, the core of [D-Gln5]SP was stabilized by the presence of two beta-turns. The exact fitting of both Phe7 and Phe8 in the binding pocket can be achieved by either an alpha- or a 3(10) helix or two beta-turns types I and II'. Replacement of an amino acid in the message sequence, 6 to 11, drastically decreased the potencies of the corresponding analogues, different conformational modifications were observed. [D-Gln6]SP presented two beta-turns, however, the types of beta-turns should orientate the side-chains in such a way that [D-Gln6]SP did not fit in the binding site. The conformations of [D-Phe7]SP and [D-Phe8]SP were completely changed, a more or less extended structure being observed. Modifications in the Gly-Leu-Met-NH2 sequence did not affect the helical structure of the core of [D-Ala9]SP, [D-Leu10]SP and [D-Met11]SP, but decreased the percentage of extended structure of the C-terminal tripeptide.

Selective agonists of tachykinin binding sites

Three types of binding sites for the mammalian tachykinins, ie Substance P (SP) Neurokinin A (NKA) and Neurokinin B (NKB), have been found in both the central and peripheral nervous systems. Substance P binds to the NK-1 subclass of binding site while NKA and NKB are less selective endogenous ligands, which preferentially interact with the NK-2 and NK-3 subclasses of binding sites, respectively. Complementary strategies, including 3-dimensional structure analysis by NMR spectroscopy and structure-activity relationship led to the design of selective agonists of these binding sites. [Pro9] SP, [Pro10] SP and the cyclic analogues [Cys3,6, Tyr8, Pro9] SP and [Cys3,6, Tyr8, Pro10] SP are selective NK-1 agonists. [Lys5] NKA(4-10) is a water soluble NK-2 potent agonist. Finally, [Pro7] NKB, which completely discriminates NK-2 and NK-3 binding sites, is a water-soluble NK-3 selective agonist.