PubMed Health⌕ Search

Biomedical subjects

Jerzy Ciarkowski

Publications and source records attributed to Jerzy Ciarkowski.

10 recordsLinked to original sources

Molecular docking-based study of vasopressin analogues modified at positions 2 and 3 with N-methylphenylalanine: influence on receptor-bound conformations and interactions with vasopressin and oxytocin receptors.

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.

Amino Acid Sequence↗

Investigation of mechanism of desmopressin binding in vasopressin V2 receptor versus vasopressin V1a and oxytocin receptors: molecular dynamics simulation of the agonist-bound state in the membrane-aqueous system.

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.

Amino Acid Motifs↗

Molecular dynamics simulation of human neurohypophyseal hormone receptors complexed with oxytocin-modeling of an activated state.

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.

Amino Acid Sequence↗

Analysis of interactions responsible for vasopressin binding to human neurohypophyseal hormone receptors-molecular dynamics study of the activated receptor-vasopressin-G(alpha) systems.

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.

Amino Acid Sequence↗

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, 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.

Dimerization↗

Use of NMR and fluorescence spectroscopy as well as theoretical conformational analysis in conformation-activity studies of cyclic enkephalin analogues.

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.

Enkephalins↗

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

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.

Amino Acid Substitution↗

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

Computer Simulation↗

Molecular modeling of the neurohypophyseal receptor/atosiban complexes.

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.

Amino Acids↗

Conformational solution studies of (Sar7)desamino- and (MeAla7)desamino-vasopressin analogues using NMR spectroscopy.

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.

Algorithms↗