Search

Human HtrA1-4 proteins belong to the HtrA family of evolutionarily conserved serine proteases and function as important modulators of many physiological processes, including maintenance of mitochondrial homeostasis, cell signaling and apoptosis. Disturbances in their action are linked to severe diseases, including oncogenesis and neurodegeneration. The HtrA1-4 proteins share structural and functional features of other members of the HtrA protein family, however there are several significant differences in structural architecture and mechanisms of action which makes each of them unique. Our goal is to present recent studies regarding human HtrAs. We focus on their physiological functions, structure and regulation, and describe current models of activation mechanisms. Knowledge of molecular basis of the human HtrAs' action is a subject of great interest; it is crucial for understanding their relevance in cellular physiology and pathogenesis as well as for using them as targets in future therapies of diseases such as neurodegenerative disorders and cancer., (Copyright © 2017 Elsevier Inc. All rights reserved.)

HtrA2(Omi) protease controls protein quality in mitochondria and plays a major role in apoptosis. Its HtrA2S306A mutant (with the catalytic serine routinely disabled for an X-ray study to avoid self-degradation) is a homotrimer whose subunits contain the serine protease domain (PD) and the regulatory PDZ domain. In the inactive state, a tight interdomain interface limits penetration of both PDZ-activating ligands and PD substrates into their respective target sites. We successfully crystalized HtrA2V226K/S306A, whose active counterpart HtrA2V226K has had higher proteolytic activity, suggesting higher propensity to opening the PD-PDZ interface than that of the wild type HtrA2. Yet, the crystal structure revealed the HtrA2V226K/S306A architecture typical of the inactive protein. To get a consistent interpretation of crystallographic data in the light of kinetic results, we employed molecular dynamics (MD). V325D inactivating mutant was used as a reference. Our simulations demonstrated that upon binding of a specific peptide ligand NH2-GWTMFWV-COOH, the PDZ domains open more dynamically in the wild type protease compared to the V226K mutant, whereas the movement is not observed in the V325D mutant. The movement relies on a PDZ vs. PD rotation which opens the PD-PDZ interface in a lid-like (budding flower-like in trimer) fashion. The noncovalent hinges A and B are provided by two clusters of interfacing residues, harboring V325D and V226K in the C- and N-terminal PD barrels, respectively. The opening of the subunit interfaces progresses in a sequential manner during the 50 ns MD simulation. In the systems without the ligand only minor PDZ shifts relative to PD are observed, but the interface does not open. Further activation-associated events, e.g. PDZ-L3 positional swap seen in any active HtrA protein (vs. HtrA2), were not observed. In summary, this study provides hints on the mechanism of activation of wtHtrA2, the dynamics of the inactive HtrA2V325D, but does not allow to explain an increased activity of HtrA2V226K., Competing Interests: The authors have declared that no competing interests exist.

Deamination of vasopressin (AVP) enhances its antidiuretic activity. Moreover, introduction of D-Arg8 instead of its L enantiomer in deamino-vasopressin (dAVP) results in an extremely potent and selective antidiuretic agonist - desmopressin (dDAVP). In this study we describe the synthesis, pharmacological properties and structures of these two potent antidiuretic agonists, and their inverso analogs. The structures of the peptides are studied in micellar and liposomic models of cell membrane using CD spectroscopy. Additionally, three-dimensional structures in mixed anionic-zwitterionic micelles are obtained using NMR spectroscopy supported by molecular dynamics simulations. Our conformational studies have shown that desmopressin in a membrane mimicking environment adopts one of the characteristic for vasopressin-like peptides β-turn - in position 3,4. Furthermore, dDAVP shows the tendency to create a β-turn in the Cys6-Gly9 C-tail, considered to be important for the antidiuretic activity, and also some tendency to adopt a 5,6 β-turn. In desmopressin, in contrast to the native vasopressin, deamino-vasopressin and [D-Arg8]-vasopressin (DAVP), the Arg8 side chain, crucial for the pressor and antidiuretic activities, is very well exposed for interaction with the receptor, whereas Gly9, crucial for the pressor and uterotonic activities, is situated together with the C-terminal amide group very close to the tocin ring. The arrangements of the Gln4 and Asn5 side chains, being crucial for OT activity, also differ in desmopressin as compared to those of AVP, dAVP and DAVP. These differences in arrangement of the important for activities side chains are likely to explain extremely potent and selective antidiuretic activities of desmopressin. © 2016 Wiley Periodicals, Inc. Biopolymers (Pept Sci) 106: 245-259, 2016., (© 2016 Wiley Periodicals, Inc.)

HtrA2(Omi) protease is involved in the maintenance of mitochondrial homeostasis and stimulation of apoptosis as well as in development of cancer and neurodegenerative disorders. The protein is a homotrimer whose subunits comprise serine protease domain (PD) and PDZ regulatory domain. In the basal, inactive state, a tight interdomain interface limits access both to the PDZ peptide (carboxylate) binding site and to the PD catalytic center. The molecular mechanism of activation is not well understood. To further the knowledge of HtrA2 thermal activation we monitored the dynamics of the PDZ-PD interactions during temperature increase using tryptophan-induced quenching (TrIQ) method. The TrIQ results suggested that during activation the PDZ domain changed its position versus PD inside a subunit, including a prominent change affecting the L3 regulatory loop of PD, and also changed its interactions with the PD of the adjacent subunit (PD*), specifically with its L1* regulatory loop containing the active site serine. The α5 helix of PDZ was involved in both, the intra- and intersubunit changes of interactions and thus seems to play an important role in HtrA2 activation. The amino acid substitutions designed to decrease the PDZ interactions with the PD or PD* promoted protease activity at a wide range of temperatures, which supports the conclusions based on the TrIQ analysis. The model presented in this work describes PDZ movement in relation to PD and PD*, resulting in an increased access to the peptide binding and active sites, and conformational changes of the L3 and L1* loops., (Copyright © 2015 Elsevier B.V. All rights reserved.)

We describe the synthesis, pharmacological properties, and structures of antidiuretic agonists, arginine vasopressin (AVP) and [D-Arg(8)]-vasopressin (DAVP), and their inverso analogues. The structures of the peptides are studied based on micellar and liposomic models of cell membranes using CD spectroscopy. Additionally, three-dimensional structures in mixed anionic-zwitterionic micelles are obtained using NMR spectroscopy and molecular dynamics simulations. NMR data have shown that AVP and DAVP tend to adopt typical of vasopressin-like peptides β-turns: in the 2-5 and 3-6 fragments. The inverso-analogues also adopt β-turns in the 3-6 fragments. For this reason, their inactivity seems to be due to the difference in side chains orientations of Tyr(2), Phe(3), and Arg(8), important for interactions with the receptors. Again, the potent antidiuretic activity of DAVP can be explained by CD data suggesting differences in mutual arrangement of the aromatic side chains of Tyr(2) and Phe(3) in this peptide in liposomes rather than of native AVP. In the presence of liposomes, the smallest conformational changes of the peptides are noticed with DPPC and the largest with DPPG liposomes. This suggests that electrostatic interactions are crucial for the peptide-membrane interactions. We obtained similar, probably active, conformations of the antidiuretic agonists in the mixed DPC/SDS micelles (5:1) and in the mixed DPPC/DPPG (7:3) liposomes. Thus it can be speculated that the anionic-zwitterionic liposomes as well as the anionic-zwitterionic micelles, mimicking the eukaryotic cell membrane environment, partially restrict conformational freedom of the peptides and probably induce conformations resembling those of biologically relevant ones.

Serum amyloid A (SAA) is a multifunctional acute-phase protein whose concentration in serum increases markedly following a number of chronic inflammatory and neoplastic diseases. Prolonged high SAA level may give rise to reactive systemic amyloid A (AA) amyloidosis, where the N-terminal segment of SAA is deposited as amyloid fibrils. Besides, recently, well-documented association of SAA with high-density lipoprotein or glycosaminoglycans, in particular heparin/heparin sulfate (HS), and specific interaction between SAA and human cystatin C (hCC), the ubiquitous inhibitor of cysteine proteases, was proved. Using a combination of selective proteolytic excision and high-resolution mass spectrometry, a hCC binding site in the SAA sequence was determined as SAA(86-104). The role of this SAA C-terminal fragment as a ligand-binding locus is still not clear. It was postulated important in native SAA folding and in pathogenesis of AA amyloidosis. In the search of conformational details of this SAA fragment, we did its structure and affinity studies, including its selected double/triple Pro → Ala variants. Our results clearly show that the SAA(86-104) 19-peptide has rather unordered structure with bends in its C-terminal part, which is consistent with the previous results relating to the whole protein. The results of affinity chromatography, fluorescent ELISA-like test, CD and NMR studies point to an importance of proline residues on structure of SAA(86-104). Conformational details of SAA fragment, responsible for hCC binding, may help to understand the objective of hCC-SAA complex formation and its importance for pathogenesis of reactive amyloid A amyloidosis., (Copyright © 2015 John Wiley & Sons, Ltd.)

Two glycosylated peptides have been studied using NMR spectroscopy supported by molecular modeling. Peptide I is an oxytocin (OT) analogue in which glutamine 4 was replaced by serine with attached α-d-mannose through the oxygen β atom, whereas peptide II is a lysine-vasopressin (LVP) analogue with lysine 8 side chain modified by the attachment of glucuronic acid through an amide bond. Both peptides exhibit very weak uterotonic effect and are less susceptible to proteolytic degradation than the mother hormones. Additionally, peptide II reveals very weak pressor and antidiuretic activities. Our results have shown that the conformational preferences of glycosylated analogues are highly similar to those of their respective mother hormones. OT glycosylated analogue (I) exhibits a 3,4 β-turn characteristic of OT-like peptides, and vasopressin-glycosylated analogue (II) exhibits β-turns typical of vasopressin-like peptides. Therefore, the lack of binding of the glycosylated analogues to the receptors can be attributed to a steric interference between the carbohydrate moieties and the receptors. We also consider this to be the reason of the very low activity of the analyzed glycopeptides. We expect that results from these studies will be helpful in designing new OT-like and vasopressin-like drugs., (Copyright © 2014 European Peptide Society and John Wiley & Sons, Ltd.)

Bacterial HtrAs are serine proteases engaged in extracytoplasmic protein quality control and are required for the virulence of several pathogenic species. The proteolytic activity of HtrA (DegP) from Escherichia coli, a model prokaryotic HtrA, is stimulated by stressful conditions; the regulation of this process is mediated by the LA, LD, L1, L2, and L3 loops. The precise mechanism of action of the LA loop is not known due to a lack of data concerning its three-dimensional structure as well as its mode of interaction with other regulatory elements. To address these issues we generated a theoretical model of the three-dimensional structure of the LA loop as per the resting state of HtrA and subsequently verified its correctness experimentally. We identified intra- and intersubunit contacts that formed with the LA loops; these played an important role in maintaining HtrA in its inactive conformation. The most significant proved to be the hydrophobic interactions connecting the LA loops of the hexamer and polar contacts between the LA' (the LA loop on an opposite subunit) and L1 loops on opposite subunits. Disturbance of these interactions caused the stimulation of HtrA proteolytic activity. We also demonstrated that LA loops contribute to the preservation of the integrity of the HtrA oligomer and to the stability of the monomer. The model presented in this work explains the regulatory role of the LA loop well; it should also be applicable to numerous Enterobacteriaceae pathogenic species as the amino acid sequences of the members of this bacterial family are highly conserved., (© 2014 by The American Society for Biochemistry and Molecular Biology, Inc.)

HtrA2(Omi), belonging to the high-temperature requirement A (HtrA) family of stress proteins, is involved in the maintenance of mitochondrial homeostasis and in the stimulation of apoptosis, as well as in cancer and neurodegenerative disorders. The protein comprises a serine protease domain and a postsynaptic density of 95 kDa, disk large, and zonula occludens 1 (PDZ) regulatory domain and functions both as a protease and a chaperone. Based on the crystal structure of the HtrA2 inactive trimer, it has been proposed that PDZ domains restrict substrate access to the protease domain and that during protease activation there is a significant conformational change at the PDZ-protease interface, which removes the inhibitory effect of PDZ from the active site. The crystal structure of the HtrA2 active form is not available yet. HtrA2 activity markedly increases with temperature. To understand the molecular basis of this increase in activity, we monitored the temperature-induced structural changes using a set of single-Trp HtrA2 mutants with Trps located at the PDZ-protease interface. The accessibility of each Trp to aqueous medium was assessed by fluorescence quenching, and these results, in combination with mean fluorescence lifetimes and wavelength emission maxima, indicate that upon an increase in temperature the HtrA2 structure relaxes, the PDZ-protease interface becomes more exposed to the solvent, and significant conformational changes involving both domains occur at and above 30 °C. This conclusion correlates well with temperature-dependent changes of HtrA2 proteolytic activity and the effect of amino acid substitutions (V226K and R432L) located at the domain interface, on HtrA2 activity. Our results experimentally support the model of HtrA2 activation and provide an insight into the mechanism of temperature-induced changes in HtrA2 structure.

Two variants of NMR-based conformational analyses of flexible peptides are compared using two examples meeting the formula Tyr-D-Daa-Phe-Daa-NH₂ (Daa=diamino acid): 1 combining D-Dab² (α,γ-diaminobutyryl) with Lys⁴, and 2 -D-Dap² (α,β-diaminopropionyl) with Orn⁴. The ω-amino groups of D-Daa² and Daa⁴ are coupled with C=O into the urea, restraining 1 and 2 with 16- and 14-membered rings and leading to potent and impotent μ/δ opioid peptides, respectively. To the current task, we took from an earlier work (Filip et al, J. Pept. Sci. 11 (2005) 347-352) the NMR NOE- and J-data in H₂O/D₂O; and the selection of the ensembles of 1 and 2, 822 and 788 conformational families, respectively, obtained by using the EDMC/ECEPP3 method. Here, we generated ensembles of 1 and 2 using AMBER molecular dynamics in explicit water to eventually selected 686 and 761 conformers for 1 and 2, respectively. We did numbers of fits for both types of the conformational ensembles of 1 and 2 to their NOE- and J-data using a common method i.e. maximum entropy approach (Groth et al, J. Biomol. NMR 15 (1999) 315-330). Both types of the well structurally diversified ensembles fit to quite different equilibria in regressions to common experimental NOE- and J-restraints using maximum entropy principle, which is a disappointing message. Intriguing is startlingly small standard deviation in J-couplings: σ(JNHαH) ≈ 0.01 Hz for LES-MD/AMBER ensemble, contrary to σ(JNHαH) = 0.8 - 1.1 Hz for the EDMC/ECEPP ensemble, over the wide range of entropy, i.e. relatively insensitive to it. A similar feature is not the case when comparing σ(NOE) in both methods. Hence, at minute entropy contributions, it follows that J does or does not transpose "overfitted" into the final σ(J) in the AMBER or ECEPP ensemble, respectively. Could this be an effect of softness of the AMBER flexible-valence force field compared to ECEPP rigid-geometry, and its effect on ensemble sampling? We do not know an answer., (Copyright © 2011 Elsevier Inc. All rights reserved.)

Gramicidin S (GS) is a cyclo-decapeptide antibiotic isolated from Bacillus brevis. The structural studies have shown that GS forms a two-stranded antiparallel β-sheet imposed by two II' β-turns. Despite its wide Gram+ and Gram- antimicrobial spectrum, GS is useless in therapy because of its high hemotoxicity in humans. It was found, however, that the analogues of GS-14 (GS with 14 amino acid residues) attained a better antimicrobial selectivity when their amphipatic moments were perturbed. In this study, we report effects of similar perturbations imposed on GS cyclo-decapeptide analogues. Having solved their structures by NMR/molecular dynamics and having tested their activities/selectivities, we have concluded that the idea of perturbation of the amphipatic moment does not work for GS-10_0 analogues. An innovative approach to the synthesis of head-to-tail cyclopeptides was used., (Copyright © 2010 European Peptide Society and John Wiley & Sons, Ltd.)

In this paper, we use NMR spectroscopy and molecular modeling to examine four new vasopressin analogs modified with alpha-2-indanylglycine (Igl) at position 2, [L-Igl(2)]AVP (I), [D-Igl(2)]AVP (II), [Mpa(1),L-Igl(2)]AVP (III) and [Mpa(1),D-Igl(2)]AVP (IV), embedded in a sodium dodecyl sulfate (SDS) micelle. All the analogs display antiuterotonic activity. In addition, the analogs with D-Igl reveal antipressor properties. Each analog exhibits the tendency to adopt beta-turns at positions 2, 3 and/or 3, 4, which is characteristic of oxytocin-like peptides. Mutual arrangement of aromatic residues at positions 2 and 3 has been found to be crucial for binding antagonists with the OT and V(1a) receptors. The orientation of the Gln(4) side chain seems to be important for the V(1a) receptor affinity. In each of the peptides studied, the Gln(4) side chain is folded back over the ring moiety. However, it lies on the opposite face of the tocin moiety in analogs with L and D enantiomers of Igl., (2010 Elsevier B.V. All rights reserved.)

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., (2010 Elsevier Inc. All rights reserved.)

The present investigation was undertaken to characterize mechanism of thermal activation of serine protease HtrA (DegP) from Escherichia coli. We monitored the temperature-induced structural changes within the regulatory loops L1, L2 and LA using a set of single-Trp HtrA mutants. The accessibility of each Trp residue to aqueous medium at temperature range 25-45 degrees C was assessed by steady-state fluorescence quenching using acrylamide and these results in combination with mean fluorescence lifetimes (tau) and wavelength emission maxima (lambda(em)max) were correlated with the induction of the HtrA proteolytic activity. Generally the temperature shift caused better exposure of Trps to the quencher; although, each of the loops was affected differently. The LA loop seemed to be the most prone to temperature-induced conformational changes and a significant opening of its structure was observed even at the lowest temperatures tested (25-30 degrees C). To the contrary, the L1 loop, containing the active site serine, remained relatively unchanged up to 40 degrees C. The L2 loop was the most exposed element and showed the most pronounced changes at temperatures exceeding 35 degrees C. Summing up, the HtrA structure appears to open gradually, parallel to the gradual increase of its proteolytic activity.

The role of the internal water molecules in vasopressin and oxytocin receptors has been investigated via molecular dynamics simulations in hydrated membrane model. Several water molecules have been identified within the binding pockets of receptors, where they interact with the conserved residues. In all unliganded receptors, the water molecules bound to the highly conserved D2.50 cluster have been observed. It has been proposed which water molecules may significantly contribute to the stability of overall receptor structure. In receptor-ligand complexes the water molecules are involved in the receptor-ligand interactions by forming water-mediated hydrogen bonds at their contact surface.

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 flexible 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.)

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., (Copyright (c) 2008 European Peptide Society and John Wiley & Sons, Ltd.)

The peptidomimetic Z-Arg-Leu-Arg-Agly-Ile-Val-OMe (where Agly means alpha-aza-glycyl, -NHNHCO-) is the strongest (K(i) = 480 pM) and the most selective inhibitor of cathepsin B to date, being approximately 2310 times as active to cathepsin B as to cathepsin K. In this paper we introduce the peptide and seek to rationalize its structure-activity relationships using molecular dynamics (MD) and NMR. It is shown that the -Agly-moiety restrains the peptide backbone to a bent shape, contrary to its parent peptide (with Gly in position 4), having its backbone extended and flexible. This fold is maintained in the plug covalently bound to the cathepsin B Cys29, in compliance with similar bends already observed in two other azapeptides attached to the active sites of cathepsin B. The MD simulation of the Z-Arg-Leu-Arg-Agly approximately cathepsin B complex suggests that, contrary to other potent inhibitors of cathepsin B, the current double Arg(1)/Arg(3) inhibitor, while maintaining the fold is able to form a unique ion cluster involving both Arg residues on the inhibitor part and two acidic Glu171 and Glu245 on the cathepsin B part, thus enhancing the affinity and subsequently the inhibiting power and selectivity of Z-Arg-Leu-Arg-Agly-Ile-Val-OMe to the observed extreme extent., (Copyright (c) 2007 European Peptide Society and John Wiley & Sons, Ltd.)

In this study, four cyclic vasopressin (CYFQNCPRG-NH(2), AVP) analogues substituted at positions 2 and 3 with four combinations of enantiomers of N-methylphenylalanine have been investigated. Three-dimensional structures of analogues have been formerly determined using NMR spectroscopy in dimethyl sulfoxide. Three-dimensional models of the vasopressin and oxytocin receptors were constructed by combining the multiple sequence alignment and the RD crystal structure as a template. The analogues have been docked into the receptor using the AutoDock program. The relaxation of the receptor-ligand complexes using energy minimization, followed by the constrained simulated annealing protocols (CSA), has been performed. The receptor-bound conformations of the investigated analogues have been proposed. We concluded that the N-methylated residues at positions 2 and 3 act as a structural restraint, determining the conformation of analogues, their location inside the receptor cavity, and mutual arrangement of the aromatic side chains. The conserved polar residues constitute the handles keeping the biologically active analogues inside the binding cavity. The Arg(8)-D(2.65) salt bridge might be responsible for analogue-selective binding in OTR and V1aR versus V2R, where the positively charged K(2.65) 100 is present at the equivalent position.

The vasopressin V2 receptor (V2R) belongs to the Class A G protein-coupled receptors (GPCRs). V2R is expressed in the renal collecting duct (CD), where it mediates the antidiuretic action of the neurohypophyseal hormone arginine vasopressin (CYFQNCPRG-NH2, AVP). Desmopressin ([1-deamino, 8-D]AVP, dDAVP) is strong selective V2R agonist with negligible pressor and uterotonic activity. In this paper, the interactions responsible for binding of dDAVP to vasopressin V2 receptor versus vasopressin V1a and oxytocin receptors has been examined. Three-dimensional activated models of the receptors were constructed using the multiple sequence alignment and the complex of activated rhodopsin with Gt(alpha) C-terminal peptide of transducin MII-Gt(alpha) (338-350) prototype (Slusarz, R.; Ciarkowski, J. Acta Biochim Pol 2004 51, 129-136) as a template. The 1-ns unconstrained molecular dynamics (MD) of receptor-dDAVP complexes immersed in the fully hydrated 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine (POPC) membrane model was conducted in an Amber 7.0 force field. Highly conserved transmembrane residues have been proposed as being responsible for V2R activation and G protein coupling. Molecular mechanism of the dDAVP binding has been suggested. The internal water molecules involved in an intricate network of the hydrogen bonds inside the receptor cavity have been identified and their role in the stabilization of the agonist-bound state proposed., (Copyright 2005 Wiley Periodicals, Inc.)

The neurohypophyseal hormone oxytocin (CYIQNCPLG-NH(2), OT) is involved in the control of labor, secretion of milk and many social and behavioral functions via interaction with its receptors (OTR) located in the uterus, mammary glands and peripheral tissues, respectively. In this paper we propose the interactions responsible for OT binding and selectivity to OTR versus vasopressin ([F3,R8]OT, AVP) receptors: V1aR and V2R, all three belonging to the Class A G protein-coupled receptors (GPCRs). Three-dimensional models of the activated receptors were constructed using a multiple sequence alignment and the activated rhodopsin-transducin (MII-Gt) prototype [Slusarz and Ciarkowski, 2004] as a template. The 1 ns unconstrained molecular dynamics (MD) of three pairs of receptor-OT complexes (two complexes per each receptor) immersed in the fully hydrated 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine (POPC) lipid bilayer was conducted in the AMBER 7.0 force field. The relaxed models of ligand-receptor complexes were used to identify the putative binding sites of OT. The stabilizing interactions with conserved Gln residues in all complexes were identified. The nonconserved hydrophobic residues were proposed as responsible for OTR-OT selectivity and ligand recognition. These results provide guidelines for experimental site-directed mutagenesis and if confirmed, they may be helpful in designing new selective OT analogs with both agonistic or antagonistic properties., (Copyright (c) 2005 European Peptide Society and John Wiley & Sons, Ltd.)

Vasopressin (CYFQNCPRG-NH(2), AVP) is a semicyclic endogenous peptide, which exerts a variety of biological effects in mammals. The main physiological roles of AVP are the regulation of water balance and the control of blood pressure and adrenocorticotropin hormone (ACTH) secretion, mediated via three different subtypes of vasopressin receptors: V1a, V1b and V2 receptors (V1aR, V1bR and V2R, respectively). They are the members of the class A, G-protein-coupled receptors (GPCRs). AVP also modulates several behavioral and social functions. In this study, the interactions responsible for AVP binding to vasopressin V1a and V2 receptors versus the closely related oxytocin ([I3,L8]AVP, OT) receptor (OTR) have been investigated. Three-dimensional models of the activated receptors were constructed using multiple sequence alignment, followed by homology modeling using the complex of activated rhodopsin with Gt(alpha) C-terminal peptide of transducin MII-Gt(338-350) prototype as a template. AVP was docked into the receptor-G(alpha) systems. The three lowest-energy pairs of receptor-AVP-G(alpha) (two complexes per each receptor) were selected. The 1-ns unconstrained molecular dynamics (MD) of complexes embedded into the fully hydrated 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine (POPC) lipid bilayer was conducted in the AMBER 7.0 force field. Six relaxed receptor-AVP-G(alpha) models were obtained. The residues responsible for AVP binding to vasopressin receptors have been identified and a different mechanism of AVP binding to V2R than to V1aR has been proposed., (Copyright (c) 2005 European Peptide Society and John Wiley & Sons, Ltd.)

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, Krzyśko 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.

A model for interaction of class A 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.

In this review the conformational studies of natural enkephalins (H-Tyr-Gly-Gly-Phe-Met-OH; the [Met(5)]enkephalin and H-Tyr-Gly-Gly-Phe-Leu-OH; the [Leu(5)]enkephalin), their acyclic and cyclic analogues, including those carried out in our laboratory, performed by experimental and theoretical methods and their combination, are described. Emphasis is given on the role of conformational constraints introduced by cyclization on activity at the micro and delta opioid receptors. Comparison of the conformations of cyclic enkephalin analogues with high delta-receptor activity with those of potent rigid non-peptide delta-receptor agonists indicates that the proximity of the aromatic side chains in positions 1 and 4 as well as the N-terminal amino group is desirable for the activity at the delta opioid receptors; early conformational studies also suggested that spatial separation of the aromatic side chains and rigidity of the cyclic backbone is desirable for micro-receptor activity. The results of our recent conformational studies performed with the use of fluorescence and NMR spectroscopy as well as theoretical calculations indicate, however, that these structural features are not necessary for activity at the micro opioid receptors. Methods applied to the determination of the conformation of flexible peptides, such as Nuclear Magnetic Resonance (NMR), fluorescence spectroscopy, and theoretical conformational analysis are also discussed briefly.

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 approximately 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-NH2 (I) with the human V1a, V2 and oxytocin receptors. The peptideI is a potent selective V1a receptor antagonist. To obtain relaxed models of the complexes, the following techniques were used: docking ofI 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.

The neurohypophyseal nonapeptide hormone oxytocin (OT) is the strongest uterotonic substance known and is responsible for the initiation of labor. Conversely, oxytocin antagonists blocking uterine OT receptor can suppress uterus contraction. In this paper we describe a computer simulated docking pertinent to affinity of an oxytocin antagonist atosiban towards OT receptor, versus vasopressin V1a and V2 receptors.

Solution structures of two analogues of vasopressin with an amino acid sequence of c[Mpa1-Tyr2-Phe3-Gln4-Asn5-Cys6]-Xaa7-Arg8-Gly9-NH2 (Xaa = Sar [I] or MeAla [II]) were established using ROE and the 3J(HNH alpha) couplings. Each of the peptides was found to exist in two stable isomers, pertaining to the cis or trans status of the Cys6-Xaa7 peptide bond, thus giving rise to four study objects. Two methods for studies of the conformational properties of the structures were compared. In the first, the algorithm consisted of three steps: (i) An Electrostatically Driven Monte-Carlo (EDMC) search for low-energy conformations. (ii) Simulations of the NOESY spectra and the vicinal couplings for these conformations. (iii) Determination of the statistical weights of the conformations with the ANALYZE package, so as to meet the best fit of the averaged NOE intensities and couplings to the experimental data. In the second method, the distance constraints and the torsion angles were used as the usual constraints in the Distance Geometry and Simulated Annealing algorithms. The flexibility of the pressin ring and the C-terminus was characterized by a large number of families of conformations. The presence of the beta-turn at position 4,5 was confirmed for all low energy structures found. The use of the EDMC method for the elaboration of the NMR data for small flexible peptides yielded an adequate description of their conformational diversity and is the method of choice for the analysis of their solution structures.