Molecular modelling of beta turns in a cyclic melanotropin

Aza-scanning of the Potent Melanocortin Receptor Agonist Ac-His-d-Phe-Arg-Trp-NH2

The melanocortin pathway is an important participant in the regulation of skin pigmentation, steroidogenesis, obesity, energy homeostasis, and exocrine gland function. Melanocortin agonists contain the putative sequence 'His-Phe-Arg-Trp', which has been designated as the 'message' sequence for melanocortin peptides, and this sequence has been hypothesized to adopt a bioactive reverse turn conformation. Exploring the relationship between its structure and biological activity, we report the synthesis and evaluation of seven aza-analogs of the potent melanocortin receptor agonist Ac-His-D-Phe-Arg-Trp-NH2. Aza-amino acids, in which the alpha-carbon was replaced by nitrogen, were inserted along the peptide sequence to probe the importance of local configuration and turn conformation on the biology of this tetrapeptide. Although systematic substitution of aza-amino acids for the D-Phe and Arg residues led to a significant loss of activity relative to the parent peptide for all melanocortin receptor subtypes examined, substitution of aza-amino acids at the C-terminal Trp residue gave analogs equipotent to the parent peptide. In summary, the aza-scan has demonstrated that recognition of this tetrapeptide by the melanocortin receptors is particularly sensitive to modifications of configuration and conformation at the D-Phe and Arg residues versus the Trp amino acid. In light of aza-residues imparting resistance from enzymatic degradation, C-terminal aza-amino acid analogs may be used to design new peptide mimics with enhanced metabolic stability.

Structure and function of the potent cyclic and linear melanocortin analogues

The MC3R and MC4R proteins comprise two melanocortin receptor subtypes that are involved in obesity, with each protein displaying a unique mechanism of action. To enable the design of a selective drug candidate, the solution structures of four peptidyl analogues of the melanocyte stimulating hormones, NDP-MSH, NDP-MSH(4-10) and two cyclic forms ([C5,C10]NDP-MSH(5-10), [C5,C10]NDP-MSH(5-11)), were characterized by two-dimensional nuclear magnetic resonance (NMR) spectroscopy and simulated annealing calculations. Using data from c-AMP assays in combination with structural analysis of melanocortin receptor/ligand models, we conclude that a lysine residue at the C-terminus of the His-Phe-Arg-Trp core sequence of melanocortin hormone is an important determinant for receptor selectivity in the both cyclic and linear MSH analogues. Our results suggest that side-chain orientation and charge-charge interactions with the ligand molecule play critical roles in receptor selectivity, whereas the overall backbone conformation or turn type contributes mainly to receptor binding.

Conformational analysis of tripeptide Ac-Lys-Pro-Val-NH2, COOH-terminal sequence of alpha-MSH

Alpha-melanocyte stimulating hormone (alpha-MSH) is an endogenous linear tridecapeptide which interacts with the melanocortin receptors (MC1-R to MC5-R) to mediate its biological effects. Antipyretic and anti-inflammatory activities of alpha-MSH are due to the COOH-terminal peptide sequence, Lys-Pro-Val (alpha-MSH[11-13]). This tripeptide might be useful as a therapeutic agent in the control of fever and inflammatory reactions. With this aim, a theoretical conformational study of the tripeptide has been carried out using molecular dynamics. The obtained conformational space has been classified into families according to the letter-code convention to partition the phi-psi map. The lowest energy conformations of each family were used as templates to design six models of conformationally constrained nonpeptide analogues.

Modeling of alpha-MSH conformations with implicit solvent

A conformational search for the most probable structures of the hormone alpha-MSH in aqueous solution was performed in order to help determine the structural features necessary for biological activity. The free-energy surface was modeled using methods from integral equation theory, and high-temperature molecular dynamics was used to enhance conformational sampling. Families of low free-energy structures have been found. The minimum energy structure shows a stable beta-turn conformation in the putative message region that is stabilized by a salt bridge between Glu5 and Lys11. The orientation of the side chains reflects the amphiphilic nature of the peptide, and a close interaction between the side chains of the His6, Phe7 and Trp9 was observed. Several structural features observed in the minimum energy structure agree well with experimental results. The conformational features led to a hypothesis of a receptor-hormone interaction model in which the hydrophobic side chains of Phe7 and Trp9 interact with the transmembrane portion of the human melanocortin (MC1) receptor. Also, the positively charged side chain of Arg8 and the imidazole side chain of His6 may interact with the negatively charged portions of the receptor which may even be on the receptor's extracellular loops.

Structure and dynamics of alpha-MSH using DRISM integral equation theory and stochastic dynamics

The structural and dynamical features of the hormone alpha-MSH in solution have been examined over a 100 ns time scale by using free energy molecular mechanics models at room temperature. The free energy surface has been modeled using methods from integral equation theory and the dynamics by the Langevin equation. A modification of the accessible surface area friction drag model was used to calculate the atomic friction coefficients. The molecule shows a stable beta-turn conformation in the message region and a close interaction between the side chains of His6, Phe7, and Trp9. A salt bridge between Glu5 and Arg8 was found not to be a preferred interaction, whereas a Glu5 and Lys11 salt bridge was not sampled, presumably due to relatively high free energy barriers. The message region was more conformationally rigid than the N-terminal region. Several structural features observed here agree well with experimental results. The conformational features suggest a receptor-hormone interaction model where the hydrophobic side chains of Phe7 and Trp9 interact with the transmembrane portion of the MC1 receptor. Also, the positively charged side chain of Arg8 and the imidazole side chain of His6 may interact with the negatively charged portions of the receptor which may even be on the receptor's extracellular loops.

NMR and quenched molecular dynamics studies of superpotent linear and cyclic alpha-melanotropins

Conformational searching, computer simulations, synthesis and NMR are used on a variety of alpha melanocyte-stimulating hormone (alpha-MSH) analogues to understand the physical characteristics required for biological potency. Peptides I (Ac-[Nle4,Asp5,D-Phe7,Lys10]alpha-MSH(4-10)-NH2), II (Ac-c[Nle4,Asp5,D-Phe7,Lys10]alpha-MSH(4-10)-NH2) and III (Ac-[Nle4,Asp5,D-Phe7,Dap10]alpha-MSH(4-10)-NH2 all show very similar conformational properties (backbone and side-chain torsional angles), and all display high biological potencies. The modeling results for these compounds are supported by the NMR data. Peptide IV (Ac-c[Nle4,Asp5,D-Phe7,Dap10]alpha-MSH(4-10)-NH2) appears to have a markedly different conformation and has decreased biological potency.

Comparative conformational studies on cyclic hexapeptides corresponding to message sequence His-Phe-Arg-Trp of alpha-melanotropin by NMR

Solution conformation of cyclo(Gly1-His2-Phe3-Arg4-Trp5-Gly6) and its D-Phe analog corresponding to the message sequence [Gly-alpha-MSH5-10] of alpha-MSH has been studied by 1D and 2D proton magnetic resonance spectroscopy in dimethyl sulfoxide (DMSO)-d6 solution and in a DMSO-d6/H2O cryoprotective mixture. The NMR data for both the analogs in solution at 300 K cannot be interpreted based on a single ordered conformation, as evidenced by the broadening of only -NH resonances as well as the temperature coefficients of the amide protons. An analysis of the nuclear Overhauser effect (NOE) cross-peaks in conjunction with temperature coefficient data indicates an equilibrium of multiple conformers with a substantial population of particular conformational states at least in the D-analog. The molecular dynamics simulations without and with NOE constraints also reveal numerous low-energy conformers with two gamma-turns, a gamma-turn and a beta-turn, two beta-turns, etc. for both the analogs. The observed NMR spectra can be rationalized by a dynamic equilibrium of conformers characterized by a gamma-bend at Gly6, two gamma-bends at Phe3 and Gly6 and a conformer with a single beta-turn and a gamma-bend for the L-Phe analog. On the other hand, a conformation with two fused beta-turns around the two tetrads His2-D-Phe3-Arg4-Trp5 and Trp5-Gly6-Gly1-His2 dominates the equilibrium mixture for the D-Phe analog. For the D-Phe analog, the experimentally observed average conformation is corroborated by molecular dynamics simulations as well as by studies in cryoprotective solvent.