Conformation-activity relationships of cyclo-constrained micro/delta opioid agonists derived from the N-terminal tetrapeptide segment of dermorphin/deltorphin.

The N-terminal tetrapeptide segments of dermorphin (Tyr-D-Ala-Phe-Gly-Tyr-Pro-Ser-NH(2)) and deltorphin (Tyr-D-Ala-Phe-Asp/Glu-Val-Val-Gly-NH(2)) are agonists at the opioid receptors micro and delta, respectively. [D-Arg(2), Lys(4)]-dermorphin-(1-4) amide (Tyr-D-Arg-Phe-Lys-NH(2), DALDA) and [Dmt(1)]DALDA (where Dmt is 2',6'-dimethyltyrosine) are among the most potent and selective micro-agonists reported to date, both in vitro (having picomolar micro receptor affinity) and in vivo. In this communication, conformation-activity studies of the following four cyclic analogs of DALDA are presented and discussed: the lead peptide S(2),S(4)-cyclo (Tyr-D-Cys-Phe-Cys-NH(2)), constrained by means of an S(4.2)--S(4.4) disulfide between Cys(2) and Cys(4); its two cis and trans C(4.2)--C(4.4)-olefinic dicarba analogs, and the product of saturation of them both. They are potent nonselective or moderately micro-selective opioid agonists in vitro.They have been synthesized and tested earlier [Berezowska I, Chung NN, Lemieux C, Wilkes BC, and Schiller PW, Acta Biochim Polon 53, 2006, 73-76]. We have studied their conformations using NMR and molecular dynamics. With major conformational constraints imposed by the 11-membered ring spanning residues 2-4, they show well defined conformations of this ring, while the exocylic Tyr(1) and Phe(3) side chains still have significant conformational freedom. The more active and selective micro versus delta disulfide and saturated dicarba agonists seem to have in common: (i) their ring structures more flexilble than those of the other two and (ii) their ring structures similar to each other and more diverse than those in the other two. Given this and the small size of the peptides having confirmed bioactivity profiles, there is a chance that their conformations determined in solution approach receptor-bound conformations. Copyright (c) 2008 European Peptide Society and John Wiley & Sons, Ltd.     

Molecular dynamics study of 4-OH-phenylacetyl-D-Y(Me)FQNRPR-NH2 selectivity to V1a receptor

G protein-coupled receptors relay diverse extracellular signals into cells via a common mechanism, involving activation of cytosol G proteins. The mechanism underlies the actions of ~50% of all drugs. In this work, we focus on simulating three protein–ligand complexes of the neurohypophyseal hormone analog 4-OH-phenylacetyl-D-Y(Me)FQNRPR-NH₂ (I) with the human V1a, V2 and oxytocin receptors. The peptide I is a potent selective V1a receptor antagonist. To obtain relaxed models of the complexes, the following techniques were used: docking of I into the vasopressin V1a, V2 and oxytocin receptor models, optimization of the geometry of the resulting complexes and molecular dynamics in a fully hydrated 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine lipid bilayer. The results of the simulations allow us to draw some conclusions about the ligand selectivity to V1aR.     

Conformation of dioxopiperazines. III. Conformational-energy calculations and evaluation of specific forces contributing to the intramolecular dioxopiperazine-phenyl ring-ring folding in cyclo(glycyl-phenylalanyl)

Conformational energy calculations on cyclo(glycyl-L-phenylalanyl), c(Gly-Phe), were carried out by means of the semiempirical MO CNDO/2 method. They corroborated earlier experimental findings that in c(Gly-Phe) the conformations with an aromatic side-chain ring folded over the dioxopiperazine (DOP) ring contribute significantly to the overall equilibrium, provided that the empirical optimization of the DOP ring geometry, extracted from x-ray data, rather than that of the CNDO/2-derived data is used throughout the computations. Following these studies, more detailed calculations, based partly on CNDO/2-derived data and partly on experimental data, were carried out to clarify the question, which are the forces responsible for this ring-ring stacking? In contrast to early suggestions, it was found that the quadrupole-quadrupole and dispersion interactions mainly contribute to folded conformations of c(Gly-Phe). Some implications of this finding on devising force fields for molecular-mechanics calculations of peptides are briefly outlined.     

AMI and PM3 study of a low molecular weight structural mimic of hydrogen exchange within the catalytic center of aspartic proteases

Based on recent X-ray studies, a low molecular weight model of the active center of aspartic proteases is proposed. The model is small enough to enable unattended geometry optimizations (including search for saddle-points) by molecular orbital methods. It consists of two malonic acid molecules and a water molecule; there is a carboxylic dimer at one end and the water molecule is located between the carboxylate and the carboxyl group at the other. The latter structure reproduces the geometry of the catalytic center of the native enzyme penicillopepsin with a root-mean-square deviation of 0.46 Å for five O--O distances. The AMI and PM3 molecular orbital methods were used to study the H-bond exchange within the model. Both methods lead consistently to the following conclusions: Among 2 pairs of symmetry-equivalent stationary states of the catalytic center there are at least 4 symmetry-independent hydrogen-exchange pathways, and many more when including symmetry of the center. Energetics and geometry of all identified pathways are presented. In summary, they result in juggling all three active center protons (COON and HOH) among all five active center oxygens (COO^–, COOH and H₂O) providing the center with a high delocalisation with respect to the actual position of its anionic site and/or its protonation status. The relevance of the delocalisation of the acidic proton to the mechanism of enzymatic action is briefly discussed. Correspondence to: J. Ciarkowski     

The role of the L2 loop in the regulation and maintaining the proteolytic activity of HtrA (DegP) protein from Escherichia coli

The aim of this study was to characterize the role of particular elements of the regulatory loop L2 in the activation process and maintaining the proteolytic activity of HtrA (DegP) from Escherichia coli. We measured the effects of various mutations introduced to the L2 loop’s region (residues 228-238) on the stability of HtrA molecule and its proteolytic activity. We demonstrated that most mutations affected the activity of HtrA. In the case of the following substitutions: L229N, N235I, I238N, the proteolytic activity was undetectable. Thus, the majority of interactions mediated by the studied amino-acid residues seem to play important role in maintaining the active conformation. Formation of contacts between the apical parts (residues 231-234) of the L2 loops within the HtrA trimer, in particular the residues D232, was shown to play a crucial role in the activation process of HtrA. Stabilization of these intermolecular interactions by substitution of D232 with valine caused a stimulation of proteolytic activity whereas deletion of this region abolished the activity. Since the pathogenic E. coli strains require active HtrA for virulence, the apical part of L2 is of particular interest in terms of structure-based drug design for treatment E. coli infections.     

Deltorphin analogs restricted via a urea bridge: structure and opioid activity.

Eight cyclic heptapeptides related to the full sequence of deltorphin have been synthesized. The synthesis of linear peptides containing diamino acid residues in positions 2 and 4 was carried out on a 4-methylbenzhydrylamine resin. Depending on protection procedures, the N-protected peptide-resins or N-protected peptide amides with free amino groups in the side chains were obtained, which were subsequently treated with bis-(4-nitrophenyl)carbonate to form a urea unit. Opioid activities of the peptides were determined in the guinea pig ileum (GPI) and mouse vas deferens (MVD) assays. Several compounds showed high delta opioid agonist potency and high selectivity for delta receptors. The results were compared with those obtained earlier for respective 1-4 deltorphin analogs. The conformations of these peptides have been studied using 2D-NMR in H2O/D2O and molecular dynamics. We observed that the backbone rings had well defined conformations, while the Tyr and Phe side chains and the C-terminal tail had significant conformational freedom. The bioassay data and conformational parameters of these peptides were compared with those of previously described, corresponding 1-4 deltorphin analogs. This comparison permitted an assessment of the role of the C-terminal peptide segment in defining the conformation and receptor interaction of the N-terminal portion and provided insight into the relationship between the putative bioactive conformations and bioactivity.     

Interaction of class A G protein-coupled receptors with G proteins

A model for interaction of classA G protein-coupled receptor with the G protein G(alpha) subunit is proposed using the rhodopsin-transducin (RD/Gt) prototype. The model combines the resolved interactions/distances, essential in the active RD*/Gt system, with the structure of Gt(alpha) C-terminal peptide bound to RD* while stabilizing it. Assuming the interactions involve conserved parts of the partners, the model specifies the conserved Helix 2 non-polar X- - -X, Helix 3 DRY and Helix 7/8 NP- -Y- - F RD* motifs interacting with the Gt(alpha) C-terminal peptide, in compliance with the structure of the latter. A concomitant role of Gt(alpha) and Gt(gamma) C-termini in stabilizing RD* could possibly be resolved assuming a receptor dimer as requisite for G protein activation.     

A hypothesis for GPCR activation

Growing evidence that rhodopsin (RD) and related G protein-coupled receptors form functional dimers/oligomers, followed by direct proof (using atomic force microscopy) that in the retina disc membrane RD associates into a paracrystalline network of rows of dimers, need models of the RD-transducin (Gt) complex that would envision an optimal RD dimer/oligomer able to satisfy all well-documented interactions with Gt. Of the models proposed so far, only a few refer to RD dimers and only one of them proposes a complex of Gt with an RD oligomer (Filipek S, Krzyko KA, Fotiadis D, Liang Y, Saperstein DA, Engel, A, Palczewski K Photochem Photobiol Sci 3: 628–638, 2004). This paper puts forward a hypothesis on another arrangement of RD monomers into the reported network of rows of dimers. Arguments for the compatibility of this set-up with interactions and activation of RD in the complex with Gt, in particular, with the well-documented movement of transmembrane helix 6 and cytosolic loop 3, which is vital for RD activation, are provided and discussed.This revised version was published online in June 2005 with corrections to the acknowledgements.     

CNDO/2 quantum-mechanical calculations of the conformational flexibility of the diketopiperazine skeleton

Eight conformers typical of diketopiperazine (DKP) ring folding were chosen for analysis. Conformational energy calculations were carried out using the semiempirical quantummechanical CNDO/2 method. The results obtained confirm considerable flexibility of the DKP skeleton. As the degree of folding increases, twisted boat conformations with the nonplanar peptide bonds tend to be more stable, while more rigid regular boat conformations with planar peptide bonds appear to be less stable than a flat one. The CNDO/2 method was found to be reliable enough for conformational studies of cyclic peptide skeletons with cis-peptide bonds.