Jerzy Ciarkowski

Exploration of the conformational space of oxytocin and arginine-vasopressin using the electrostatically driven Monte Carlo and molecular dynamics methods

Conformational analysis of the neurohypophyseal hormones oxytocin (OT) and arginine-vasopressin (AVP) has been carried out using two different computational approaches and three force fields, namely by the Electrostatically Driven Monte Carlo (EDMC) method, with the Empirical Conformational Energy Program for Peptides (ECEPP/3) force field or with the ECEPP/3 force field plus a hydration-shell model, and by simulated-annealing molecular dynamics with the Consistent Valence Force Field (CVFF). The low-energy conformations obtained for both hormones were classified using the minimal-tree clustering algorithm and characterized according to the locations of beta-turns in the cyclic moieties. Calculations with the CVFF force field located conformations with a beta-turn at residues 3 and 4 as the lowest energy ones both for OT and for AVP. In the ECEPP/3 force field the lowest energy conformation of OT contained a beta-turn at residues 2 and 3, conformations with this location of the turn being higher in energy for AVP. The latter difference can be attributed to the difference in the size of the side chain in position 3 of the sequences: the bulkier phenylalanine residue of AVP in combination with the bulky Tyr2 residue hinders the formation of a turn at residues 2 and 3. Conformations of OT and AVP with a turn at residues 3,4 were in the best agreement with the x-ray structures of deaminooxytocin and pressinoic acid (the cyclic moiety of vasopressin), respectively, and with the nmr-derived distance constraints. Generally, the low-energy conformations obtained with the hydration-shell model were in a better agreement with the experimental data than the conformations calculated in vacuo. It was found, however, that the obtained low-energy conformations do not satisfy all of the nmr-derived distance constraints and the nuclear Overhauser effect pattern observed in nmr studies can be fully explained only by assuming a dynamic equilibrium between conformations with beta-turns at residues 2,3, 3,4, and 4,5. The low-energy structures of OT with a beta-turn at residues 2,3 have the disulfide ring conformations close to the model proposed recently for a potent bicyclic antagonist of OT [M. D. Shenderovich et al. (1994) Polish Journal of Chemistry, Vol. 25, pp. 921-927], although the native hormone differs from the bicyclic analogue by the conformation of the C-terminal tripeptide. This finding confirms the hypothesis of different receptor-bound conformations of agonists and antagonists of OT.

Oxytocin receptors and cholesterol: interaction and regulation

Cholesterol affects the ligand binding function of the oxytocin receptor in a highly specific manner. While the structurally‐related cholecystokinin receptor shows a strong correlation between the membrane fluidity and its binding function, the oxytocin receptor behaves differently. A stringent and unique requirement of the affinity state of the oxytocin receptor for structural features of the sterol molecule has been found. The molecular requirements differ both from those postulated for sterol—phospholipid interactions and from those known to be necessary for the activity of other proteins. Employing a new detergent‐free subcellular fractionation protocol, a two‐fold enrichment of the oxytocin receptors (10–15% of total receptors) has been detected in the cholesterol‐rich, caveolin‐containing membrane domains of the plasma membrane. While most of the properties of the oxytocin receptors were indistinguishable in cholesterol‐poor versus cholesterol‐rich membrane compartments, high‐affinity oxytocin receptors localised in caveolin‐enriched low‐density membranes showed about a 3‐fold higher stability against thermal denaturation at 37°C compared with the oxytocin receptors localised in high‐density membranes. Moreover, addition of cholesterol to the cholesterol‐poor high‐density membranes fully protected the oxytocin receptors against thermal denaturation and partially rescued high‐affinity oxytocin binding. Although the membrane fluidity of the caveolin‐enriched domains was lower than that in the high‐density membranes, there was no correlation between the stability of oxytocin receptors and the fluidity level of the membrane domains. Finally, in a molecular modelling approach a putative cholesterol binding motif on the extracellular surface of the oxytocin receptor was found.

Molecular modeling of the human vasopressin V2 receptor/agonist complex.

The V2 vasopressin renal receptor (V2R), which controls antidiuresis in mammals, is a member of the large family of heptahelical transmembrane (7TM) G protein-coupled receptors (GPCRs). Using the automated GPCR modeling facility available via Internet (http://expasy.hcuge.ch/swissmod/SWISS-MODEL.html) for construction of the 7TM domain in accord with the bovine rhodopsin (RD) footprint, and the SYBYL software for addition of the intra- and extracellular domains, the human V2R was modeled. The structure was further refined and its conformational variability tested by the use of a version of the Constrained Simulated Annealing (CSA) protocol developed in this laboratory. An inspection of the resulting structure reveals that the V2R (likewise any GPCR modeled this way) is much thicker and accordingly forms a more spacious TM cavity than most of the hitherto modeled GPCR constructs do, typically based on the structure of bacteriorhodopsin (BRD). Moreover, in this model the 7TM helices are arranged differently than they are in any BRD-based model. Thus, the topology and geometry of the TM cavity, potentially capable of receiving ligands, is in this model quite different than it is in the earlier models. In the subsequent step, two ligands, the native [arginine8]vasopressin (AVP) and the selective agonist [d-arginine8]vasopressin (DAVP) were inserted, each in two topologically non-equivalent ways, into the TM cavity and the resulting structures were equilibrated and their conformational variabilities tested using CSA as above. The best docking was selected and justified upon consideration of ligand-receptor interactions and structure-activity data. Finally, the amino acid residues were indicated, mainly in TM helices 3-7, as potentially important in both AVP and DAVP docking. Among those Cys112, Val115-Lys116, Gln119, Met123 in helix 3; Glu174 in helix 4; Val206, Ala210, Val213-Phe214 in helix 5; Trp284, Phe287-Phe288, Gln291 in helix 6; and Phe307, Leu310, Ala314 and Asn317 in helix 7 appeared to be the most important ones. Many of these residues are invariant for either the GPCR superfamily or the neurophyseal (vasopressin V2R, V1aR and V1bR and oxytocin OR) subfamily of receptors. Moreover, some of the equivalent residues in V1aR have already been found critical for the ligand affinity [Mouillac et al., J. Biol. Chem, 270 (1995) 25771].

Temperature-induced conformational changes within the regulatory loops L1 -L2-LA of the HtrA heat-shock protease from Escherichia coli

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.

Binding modes of a new epoxysuccinyl–peptide inhibitor of cysteine proteases. Where and how do cysteine proteases express their selectivity?

Papain from Carica papaya, an easily available cysteine protease, is the best-studied representative of this family of enzymes. The three dimensional structure of papain is very similar to that of other cysteine proteases of either plant (actinidin, caricain, papaya protease IV) or animal (cathepsins B, K, L, H) origin. As abnormalities in the activities of mammalian cysteine proteases accompany a variety of diseases, there has been a long-lasting interest in the development of potent and selective inhibitors for these enzymes. A covalent inhibitor of cysteine proteases, designed as a combination of epoxysuccinyl and peptide moieties, has been modeled in the catalytic pocket of papain. A number of its configurations have been generated and relaxed by constrained simulated annealing-molecular dynamics in water. A clear conformational variability of this inhibitor is discussed in the context of a conspicuous conformational diversity observed earlier in several solid-state structures of other complexes between cysteine proteases and covalent inhibitors. The catalytic pockets S2 and even more so S3, as defined by the pioneering studies on the papain-ZPACK, papain-E64c and papain-leupeptin complexes, appear elusive in view of the evident flexibility of the present inhibitor and in confrontation with the obvious conformational scatter seen in other examples. This predicts limited chances for the development of selective structure-based inhibitors of thiol proteases, designed to exploit the minute differences in the catalytic pockets of various members of this family. A simultaneous comparison of the three published proenzyme structures suggests the enzymes prosegment binding loop-prosegment interface as a new potential target for selective inhibitors of papain-related thiol proteases.

Origin of the ring-ring interaction in cyclic dipeptides incorporating an aromatic amino acid

Abstract Cyclic dipeptides comprising an aromatic amino acid residue demonstrate the energetically most favored conformation in which an aromatic ring stacks over the 2,5-dioxopiperazine ring. The communication is the first report providing a more detailed insight into the mechanism of this stacking.

The LA Loop as an Important Regulatory Element of the HtrA (DegP) Protease from Escherichia coli STRUCTURAL AND FUNCTIONAL STUDIES

Background: An understanding of the HtrA protease activation mechanism is incomplete with respect to its LA regulatory loop. Results: A theoretical model of the LA structure is provided and experimentally verified. Conclusion: LA intersubunit contacts strongly contribute to the stabilization of the inactive HtrA. Significance: This is the first report that simultaneously offers a theoretical three-dimensional structure of LA and its biophysical and functional properties. 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.

Synthesis, biological activity and solution structure of new analogues of the antimicrobial Gramicidin S.

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

Prediction of conformation of rat galanin in the presence and absence of water with the use of monte Carlo methods and the ECEPP/3 force field

The conformation of the 29-residue rat galanin neuropeptide was studied using the Monte Carlo with energy minimization (MCM) and electrostatically driven Monte Carlo (EDMC) methods. According to a previously elaborated procedure, the polypeptide chain was first treated in a united-residue approximation, in order to enable extensive exploration of the conformational space to be carried out (with the use of MCM), Then the low-energy united-residue conformations were converted to the all-atom representations, and EDMC simulations were carried out for the all-atom polypeptide chains, using the ECEPP/3 force field with hydration included. In order to estimate the effect of environment on galanin conformation, the low-energy conformations obtained as a result of these simulations were taken as starting structures for further EDMC runs that did not include hydration. The lowest-energy conformation obtained in aqueous solution calculations had a nonhelical N-terminal part packed against the nonpolar face of a residual helix that extended from Pro13 toward the C-terminus. One next lowest-energy structure was a nearly-all-helical conformation, but with a markedly higher energy. In contrast, all of the low-energy conformations in the absence of water were all-helical differing only by the extent to which the helix was kinked around Pro13. These results are in qualitative agreement with the available NMR and CD data of galanin in aqueous and nonaqueous solvents.

Intra- and intersubunit changes accompanying thermal activation of the HtrA2(Omi) protease homotrimer.

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.