Conformational and dynamic considerations in the design of peptide hormone analogs

Efforts to understand the chemical-physical basis for peptide hormone and neurotransmitter action requires integration of conformational parameters and biological properties. Since most peptide hormones are conformationally flexible, the question arises as to which of the manifold of conformations is of biological significance. In molecular terms, it is necessary to carefully distinguish chemical-physical features important to binding (the binding message) from those involved in transduction (the biological activity message). One approach to this involves the design, synthesis, and conformational analysis of semirigid hormone analogs. The distinction between binding and transduction can best be examined by evaluation of full biological profiles of partial agonists, antagonists, and analogs with prolonged biological activity. Using this multidisciplinary approach, we have prepared several semirigid [Pen¹]-oxytocin antagonist analogs and evaluated their conformational properties and biological activities. Specific conformational features can be related to inhibitory activities in several cases. On the basis of structure–activity relationships and conformational considerations, we have designed a series of conformationally restricted cyclic and acyclic analogs of the linear peptide α-melanotropin. Some of these peptides have exceptionally prolonged in vivo activity (weeks), and others exhibit superagonist potency (10,000 times the native hormone). We have evidence that potency and prolonged activity have different structural and conformational requirements. It is suggested that potency is primarily a function of receptor recognition (the binding message), whereas prolonged activity is related to transduction (the biological activity message).     

Reversible scaling of dihedral angle barriers during molecular dynamics to improve structure prediction of cyclic peptides.

Peptide cyclization or chemical cross-linking has frequently been used to restrict the conformational freedom of a peptide, for example, to enhance its capacity for selective binding to a target receptor molecule. Structure prediction of cyclic peptides is important to evaluate possible conformations prior to synthesis. Because of the conformational constraints imposed by cyclization low energy conformations of cyclic peptides can be separated by large energy barriers. In order to improve the conformational search properties of molecular dynamics (MD) simulations a potential scaling method has been designed. The approach consists of several consecutive MD simulations with a specific lowering of dihedral energy barriers and reduced nonbonded interactions between atoms separated by three atoms followed by gradually scaling the potential until the original barriers are reached. Application to four cyclic penta- and hexa-peptide test cases and a protein loop of known structure indicates that the potential scaling method is more efficient and faster in locating low energy conformations than standard MD simulations. Combined with a generalized Born implicit solvation model the low energy cyclic peptide conformations and the loop structure are in good agreement with experiment. Applications in the presence of explicit water molecules during the simulations showed also improved convergence to structures close to experiment compared with regular MD.     

NMR and molecular modeling characterization of RGD containing peptides.

The tripeptide sequence arginine-glycine-aspartic acid (RGD) has been shown to be the key recognition segment in numerous cell adhesion proteins. The solution conformation and dynamics in DMSO-d6 of the cyclic pentapeptides, [formula: see text], a potent fibrinogen receptor antagonist, and [formula: see text], a weak fibrinogen receptor antagonist, have been characterized by nuclear magnetic resonance (NMR) spectroscopy and molecular modeling. 1H-1H distance constraints derived from two-dimensional NOE spectroscopy and torsional angle constraints obtained from 3JNH-H alpha coupling constants, combined with computer-assisted modeling using conformational searching algorithms and energy minimization have allowed several low energy conformations of the peptides to be determined. Low temperature studies in combination with molecular dynamics simulations suggest that each peptide does not exist in a single, well-defined conformation, but as an equilibrating mixture of conformers in fast exchange on the NMR timescale. The experimental results can be fit by considering pairs of low energy conformers. Despite this inherent flexibility, distinct conformational preferences were found which may be related to the biological activity of the peptides.     

Accessible conformations of melanin-concentrating hormone: a molecular dynamics approach

Molecular dynamics simulations have been used to search for the accessible conformations of the melanin-concentrating hormone (MCH). The studies have been performed on native MCH and two of its peptide fragments, a cyclic MCH(5-14) fragment and a linear MCH(5-14) fragment. An analysis of the molecular dynamics trajectories of the three peptides indicates that two regions of the peptide have characteristic conformational properties that may be important for the biological activity. One is a region around Gly8, which is conformationally mobile, and the other is around Pro13, which shows unusual rigidity. The molecular dynamics simulation results are discussed in terms of backbone structural features like beta turns, side-chain interactions, and orientations of the disulfide bridge. The results of this analysis are used to suggest new analogues that will modify the conformational features of the peptide and further define the conformational requirements for activity. Finally, the results are related to nmr studies of the peptide and reveal agreements between the experimental nuclear Overhauser effect constraints and some of the accessible conformations obtained from the simulation.     

Role of the conformational element in peptide-receptor interactions. Studies with cyclic opioid peptide analogs

Biological activity profiles of three different families of cyclic opioid peptide analogs are presented. It is illustrated that conformational constraints introduced through peptide cyclizations can have drastic effects on receptor affinity, selectivity and 'efficacy' ('intrinsic activity'). Conformational studies of cyclic opioid peptides by various physico-chemical techniques have been initiated and have already produced insight into the conformational requirements of the various opioid receptor types. On the basis of the results obtained, conformational restriction of opioid peptides may represent a first promising step towards the goal of developing peptide mimetics.     

Cyclic peptides — Small and big and their conformational aspects

Cyclic peptides form an interesting class of compounds for study by conformational analysis, by virtue of their unique conformational features and biological properties. The small cyclic peptides having 3-6 peptide units in their ring, show a variety of conformational characteristics such as occurrence ofcis peptide units, flexibility of peptide dimension and variety in hydrogen bonding. The different possible conformations of cyclic tri- and hexa-peptides are given and certain specific conformational features are discussed for cyclic tetra and pentapeptides. For higher cyclic peptides, the hydrogen bonding requirement for stability of the backbone of the ring, is seen to be kept to a minimum. These various features and their significance are examined and discussed in the light of energy minimization studies and analysis of available experimental data.     

Improvement of biological activity and proteolytic stability of peptides by coupling with a cyclic peptide

The cyclic moiety of an endothelin antagonist peptide RES-701-1, composed of 10 amino acids with an amide bond between alpha-NH(2) of Gly1 and beta-COOH of Asp9, was coupled to some biologically active peptides aiming to improve their activities and stabilities against proteolytic degradation. Coupling of the cyclic peptide to the N-terminal of RGD-peptides, maximally 4-fold improvement of in vitro activity compared to the original peptide has been achieved. Coupling of it to protein farnesyltransferase inhibiting peptides resulted to improve in vitro activity maximally 3-fold. These peptides coupled with the cyclic peptide also showed enhanced stability against some typical proteases. These results indicate that this cyclic peptide can stabilize the conformations of the peptides coupled to its C-terminus. Coupling of our cyclic peptide is anticipated to be a novel conformational stabilizing method for biologically active peptides, results to improve their activity and stability.     

Modeling an active conformation for linear peptides and design of a competitive inhibitor for HMG-CoA reductase

This study presents an approach that can be used to search for lead peptide candidates, including unconstrained structures in a recognized sequence. This approach was performed using the design of a competitive inhibitor for 3-hydroxy-3-methylglutaryl CoA reductase (HMGR). In a previous design for constrained peptides, a head-to-tail cyclic structure of peptide was used as a model of linear analog in searches for lead peptides with a structure close to an active conformation. Analysis of the conformational space occupied by the peptides suggests that an analogical approach can be applied for finding a lead peptide with an unconstrained structure in a recognized sequence via modeling a cycle using fixed residues of the peptide backbone. Using the space obtained by an analysis of the bioactive conformations of statins, eight cyclic peptides were selected for a peptide library based on the YVAE sequence as a recognized motif. For each cycle, the four models were assessed according to the design criterion ("V" parameter) applied for constrained peptides. Three cyclic peptides (FGYVAE, FPYVAE, and FFYVAE) were selected as lead cycles from the library. The linear FGYVAE peptide (IC(50) = 0.4 microM) showed a 1200-fold increase the inhibitory activity compared to the first isolated LPYP peptide (IC(50) = 484 microM) from soybean. Experimental analysis of the modeled peptide structures confirms the appropriateness of the proposed approach for the modeling of active conformations of peptides.     

Structure-activity relationship and metabolic stability studies of backbone cyclization and N-methylation of melanocortin peptides

Backbone cyclization (BC) and N-methylation have been shown to enhance the activity and/or selectivity of biologically active peptides and improve metabolic stability and intestinal permeability. In this study, we describe the synthesis, structure-activity relationship (SAR) and intestinal metabolic stability of a backbone cyclic peptide library, BL3020, based on the linear alpha-Melanocyte stimulating hormone analog Phe-D-Phe-Arg-Trp-Gly. The drug lead, BL3020-1, selected from the BL3020 library (compound 1) has been shown to inhibit weight gain in mice following oral administration. Another member of the BL3020 library, BL3020-17, showed improved biological activity towards the mMC4R, in comparison to BL3020-1, although neither were selective for MC4R or MC5R. N-methylation, which restrains conformational freedom while increasing metabolic stability beyond that which is imparted by BC, was used to find analogs with increased selectivity. N-methylated backbone cyclic libraries were synthesized based on the BL3020 library. SAR studies showed that all the N-methylated backbone cyclic peptides demonstrated reduced biological activity and selectivity for all the analyzed receptors. N-methylation of active backbone cyclic peptides destabilized the active conformation or stabilized an inactive conformation, rendering the peptides biologically inactive. N-methylation of backbone cyclic peptides maintained stability to degradation by intestinal enzymes.     

Conformational studies of antimetastatic laminin-1 derived peptides in different solvent systems, using solution NMR spectroscopy.

Due to its critical role in cancer progression, interactions between laminin-1 and the 67 kDa Laminin-Binding Protein (the 67 kDa LBP) have been the focus of a number of structural and biological studies. As laminin-1 is such a large and complex molecule, research interests have turned to the investigation of bioactive peptides derived from binding domains of laminin-1. Two peptides of interest, CDPGYIGSR (peptide 11) and YIGSR, both derived from the beta1 chain of laminin-1, have been shown to block invasion of basement membranes by tumor cells. Substituting the C-terminal arginine to lysine, a conservative substitution, results in a loss of peptide antimetastatic activity. This difference in bioactivity has been attributed, based on numerous modeling studies of free peptide conformations, to structural differences between YIGSR and YIGSK. Yet the nature of the 'active' free peptide backbone conformation has been a matter of debate and controversy. In order to test the validity of the structural modeling claims, we have undertaken detailed conformational studies of the two laminin-1 derived peptides YIGSR and CDPGYIGSR along with the biologically inactive YIGSK analog by two-dimensional solution 1H NMR spectroscopy in three different solvent systems. Herein we report that although both the active (YIGSR, CDPGYIGSR) and the inactive (YIGSK) peptides can adopt several closely related conformations in solution, the two peptides share similar conformational preferences, and there are no significant structural differences between the active and inactive peptides, contrary to previously reported modeling data. We conclude that the basis of the peptide biological activity, in contrast to published models, cannot be attributed to well-defined structural preferences of the free peptides. We infer that the difference in bioactivity observed between YIGSR and YIGSK originates primarily from the chemical nature of the arginine versus lysine sidechain substitution, rather than being due to a structural change in the free peptide conformations.     

An approach to the problem of the structuro-functional organization of natural peptides

Theory and computational scheme of three-dimensional structure and dynamic conformational properties of naturally occurring peptides are proposed basing on a known amino acid sequence. The diverse biological activity of a low-molecular peptide is shown to arise from a restricted number of preferable spatial structures which may occur under physiological conditions. Each particular function of an oligopeptide is connected to a definite spatial structure, belonging to the set of low-energy conformations from one biological activity of a peptide shift of the conformational equilibrium caused by a change of environmental conditions. This shift is provided for by specific intramolecular interactions, alternative in their nature, which stabilize a particular structure. An approach is suggested which enables to construct a synthetic analog with the predetermined physiologically active conformation, prior to all chemical and biological tests.     

Solution conformational study of nociceptin an(d its 1-13 and 1-11 fragments using circular dichroism and two-dimensional NMR in conjunction with theoretical conformational analysis.

Conformational studies of nociceptin (NC-NH2), its fully active fragment, NC(1-13)-NH2, and two significantly less potent fragments, NC(1-13)-OH and NC(1-11)-OH, were conducted in water and TFE solutions by the employment of circular dichroism, and in DMSO-d6 by 2DNMR spectroscopy in conjunction with theoretical conformational analysis. The conformations of all thepeptides studied were calculated taking two approaches. The first assumes multiconformational equilibrium of the peptide studied, which is characterized by a set of conformations (and their statistical weight values)obtained from a global conformational analysis using three methods: the electrostatically driven Monte-Carlo (EDMC) with the ECEPP/3 force field, the simulated annealing (SA) protocols in the AMBER and CHARMM force fields. The second approach incorporates the interproton distance and dihedral angle constraints into the starting conformation. Calculations were performed using the distance geometry and SA protocol in the CHARMM force field implemented in the X-PLOR program. The CD experiments indicated that for the active peptides, hydrophobic solvents induced a significantly higher (compared with those remaining)content order, probably a helical structure. Unfortunately, as a result of the conformational flexibility of thepeptides, the analysis of conformations obtained with both approaches and different force fields did not alllow the selection of any structural elements of the NC peptides that might be connected with their bioactivity. The only common element found in most conformations of the active peptides was a helical character of fragment 8-13, which allowed the side chains of basic amino acid residues to be exposed to the outside of the molecule and probably to interact with the ORL1 receptor.     

Conformational and dynamic considerations in peptide structure-function studies

Most small peptide hormones and neurotransmitters are highly flexible, conformationally labile molecules in aqueous and other environments. Thus efforts to determine the relationships between conformational properties of these peptides in aqueous and other solvents and their biological activities at membrane receptors have been difficult and of limited success. One approach which may provide a more rational basis for conformation-activity relationships is the design of conformationally restricted, semi-rigid analogs of the native peptides which still possess high potency and/or antagonist properties. In addition to the increased likelihood that the conformational properties determined for these derivatives in aqueous or other solvent environments will have biological relevance, such analogs are likely to have higher specificity for particular receptors, greater in vivo stability, and perhaps even oral activity. The application of this approach to the design of highly potent oxytocin antagonists is discussed with particular emphasis on the conformational and dynamic properties which appear to differentiate agonist and antagonist analogs. The results of these studies are briefly compared with similar studies with somatostatin, angiotensin, bradykinin, α-melanotropin and enkephalin, and discussed in terms of likely further developments.     

Conformation-activity relationship of peptide T and new pseudocyclic hexapeptide analogs.

Peptide T (ASTTTNYT), a segment corresponding to residues 185-192 of gp120, the coat protein of HIV, has several important biological properties in vitro that have stimulated the search for simpler and possibly more active analogs. We have previously shown that pseudocyclic hexapeptide analogs containing the central residues of peptide T retain considerable chemotactic activity. We have now extended the design of this type of analogs to peptides containing different aromatic residues and/or Ser in lieu of Thr. The complex conformation-activity relationship of these analogs called for a reexamination of the basic conformational tendencies of peptide T itself. Here, we present an exhaustive NMR conformational study of peptide T in different media. Peptide T assumes a gamma-turn in aqueous mixtures of ethylene glycol, a type-IV beta-turn conformation in aqueous mixtures of DMF, and a type-II beta-turn conformation in aqueous mixtures of DMSO. The preferred conformations for the analogs were derived from modeling, starting from the preferred conformations of peptide T. The best models derived from the gamma-turn conformation of peptide T are those of peptides XII (DSNYSR), XIII (ETNYTK) and XVI (ESNYSR). The best models derived from the type-IV beta-turn conformation of peptide T are those of peptides XIV (KTTNYE) and XV (DSSNYR). No low-energy models could be derived starting from the type-II beta-turn conformation of peptide T. The analogs with the most favored conformations are also the most active in the chemotactic test.     

The Relationship Between Structure and Activity Among Opioid Peptides

Summary Since the discovery and isolation of the endogenous opioid peptides Leu- and Met-enkephalin, structural studies have been focused on deducing the bioactive conformation of the peptide ligands. Theoretically, linear peptides can have many different backbone conformations, yet early, X-ray studies on enkephalin and its analogues showed only two different backbone conformations: extended and single β-bend. More recent reports include a third conformation for Leu-enkephalin and constrained opioid peptides from two ‘new’ classes (i.e. cyclic and ‘allaromatic’ peptides). In this report the relationship between solid-state X-ray structure and opioid peptide activity is examined. The N-terminal amine nitrogen and the two aromatic rings have previously been identified as structural features important to the biological activity of opioid peptides. From X-ray studies we find that the distances between the centroids of the aromatic rings, and between the N-terminal amine nitrogen and the centroid of the phenylalanine ring, vary over a large range. There is a discernible relationship, however, between the separation of the two rings and their orientation that correlates with activity.     

NMR study of cyclic peptides with renin inhibitor activity

Several cyclic analogues of renin inhibitors, based on Glu-D-Phe-Lys motif have been investigated by NMR spectroscopy and molecular dynamics calculations (MD). The 15 membered macrocycle, resulting from Glu and Lys side-chain cyclization, exhibits conformational preference. The structural evidence from NMR shows the presence of hydrogen bond between Lys NH and Glu side-chain carbonyl, resulting in a 10 membered pseudo beta-turn-like structure. The structure of the cyclic moiety is similar in all the peptides, which takes at least two conformations around Calpha-Cbeta in Glu side chain. The restrained MD calculations further support such observations and show that the macrocycle is fairly rigid, with two conformations about the Glu Calpha-Cbeta bond. The linear peptide appendages, which are essential for activity in cyclic peptides, show an extended structure in the beta-region of Ramchandran plot. These calculations also demonstrate that for the most active peptide, two major conformers each exist about the Calpha-CO bond of the Lys, D-Trp and Leu residues. In this peptide, the cyclic moiety presents a negatively charged surface formed due to the carbonyl oxygens, which are thus available to form hydrogen bonds with the receptor. The linear fragment presents further binding sites with a surface which has the hydrophobic side chains of D-Trp, Leu and D-Met on one side and carbonyls on the other side.     

Photomodulation of conformational states. IV. Integrin-binding RGD-peptides with (4-aminomethyl)phenylazobenzoic acid as backbone constituent

In previous studies we have investigated octapeptides backbone-cyclized by (4-amino)phenyl azobenzoic acid (APB) or (4-aminomethyl)phenylazobenzoic acid (AMPB) and containing the active-site sequence Cys-Ala-Thr-Cys-Asp from the thioredoxin reductase. The conformational and redox properties of these peptides were strongly dependent on the isomeric state of the azobenzene chromophore. Using the same approach we were successful in constructing photoresponsive ligands for alphavbeta3 integrin containing the Arg-Gly-Asp (RGD) sequence as binding motif. For achieving maximal conformational restriction of the peptide a reduced ring size compared to our previous azobenzene peptides was employed in the cyclic peptide c[Asp-D-Phe-Val-AMPB-Lys-Ala-Arg-Gly-]. Conformational properties of the trans and cis isomers of this peptide in solution were investigated by CD and NMR and were found to differ markedly from the thioredoxin derived azobenzene peptides. In a second peptide, c[Asp-D-Phe-Val-Lys-AMPB-Ala-Arg-Gly-], shifting the position of the chromophore lead to a marked decrease in affinity. With the availability of the x-ray structure of a cyclic RGD-pentapeptide bound to alphavbeta3 integrin (PDB entry 1L5G) modeling of possible bound conformations for trans and cis isomers of both azobenzene peptides was possible. Notably, both peptides in either isomeric form share the same overall conformation in the bound state according to our molecular dynamics simulations.     

Conformation/activity studies of rationally designed potent anti-adhesive RGD peptides.

The Arg-Gly-Asp (RGD) sequence is a universal cell-recognition site of various extracellular proteins that interact with integrin cell-surface receptors. In order to design low-molecular-mass RGD protein antagonists, the determination of the biologically active conformation is a prerequisite. We present a method that yields detailed insight into the steric factors which govern the binding of the ligands to their receptors by systematically scanning the conformational space accessible for the tripeptide sequence RGD. The investigation is based on the conformationally controlled design of homodetic cyclic oligopeptides and their structural determination, coupled with biological assays. For this purpose, a whole set of cyclic pentapeptides and hexapeptides has been synthesized and their three-dimensional structures in solution analyzed by modern two-dimensional NMR techniques in combination with restrained and free molecular dynamics simulations. Their biological activity was compared with that of linear GRGDS in inhibition assays of tumor cell adhesion to laminin P1 and vitronectin substrates. An up to 100-fold, and in part selective, increase in activity was observed for two cyclic pentapeptides. Most other peptides showed a decreased activity which, however, was useful to correlate activity with rather small variations in conformation. Detailed comparative studies of the systematically designed conformations and the corresponding anti-adhesive activities offer an access to lead structures for a rational indirect drug design of peptide and peptidomimetic pharmaceuticals with strong interfering activity for integrin-mediated cell-cell and cell-matrix interactions.     

Conformations and pharmacophores of cyclic RGD containing peptides which selectively bind integrin alpha(v)beta3.

This paper reports a detailed conformational characterization in solution by 1H-NMR in H2O and DMSO-d6 and molecular modeling simulations of cyclic peptides containing the RGDDV pharmacophore and the RGDY(Me)R pharmacophore. These two pentapeptide sequences when properly constrained in cyclic peptides are low to sub-nanomolar inhibitors of integrin alpha(v)beta3. The peptides containing the RGDDY(Me)R sequence bind potently to integrin alphaIIb3 as well. The conformations found in H2O and in DMSO-d6 solutions are valuable for the design of peptidomimetics of these two pharmacophores. The structure-activity relationships of the RGDDV and RGDY(Me)R pharmacophores within cyclic peptides are discussed. Specifically, the orientation of surface-accessible chemical features on the ligand, such as hydrophobic, positive and negative ionizable groups, which are considered to be responsible for the desired biological activity, is focused on.