Molecular modeling of the oxytocin receptor/bioligand interactions

Oxytocin is a nonapeptide hormone (CYIQNCPLG-NH2, OT), controlling labor and lactation in mammalian females, via interactions with specific cellular membrane receptors (OTRs). The native hormone is cyclized via a 1-6 disulfide and its receptor belongs to the GTP-binding (G) protein-coupled receptor (GPCR) family, also known as heptahelical transmembrane (7TM) or serpentine receptors. Using a technique combining multiple sequence alignments with available experimental constraints, a reliable OTR model was built. Subsequently, the OTR complexes with a selective agonist [Thr4,Gly7]OT, a selective cyclohexapeptide antagonist L-366,948 and oxytocin itself were modeled and relaxed using a constrained simulated annealing (CSA) protocol. All three ligands seem to prefer similar modes of binding to the receptor, manifested by repeating receptor residues which directly interact with the ligands. Those involved in the three complexes are putative helices: TM3: R113, K116, Q119, M123; TM4: Q171, and TM5: I201 and T205. Most of them are the equivalent residues/positions to those found in our earlier studies, regarding related vasopressin V2 receptor/bioligand interactions.     

Theoretical studies of binding modes of two covalent inhibitors of cysteine proteases

Physiological and pathological roles of cysteine proteases make them important targets for inhibitor development. Although highly potent inhibitors of this group of enzymes are known, their major drawback is a lack of sufficient specificity. Two cysteine protease covalent inhibitors, viz. (i) Z-RL-deoxo-V-peptide-epoxysuccinyl hybrid, and (ii) Z-RLVG-methyl-, have been developed and modeled in the catalytic pocket of papain, an archetypal thiol protease. A number of configurations have been generated and relaxed for each system using the AMBER force field. The catalytic pockets S3 and S4 appear rather elusive in view of the observed inhibitors' flexibility. This suggest rather limited chances for the development of selective structure-based inhibitors of thiol proteases, designed to exploit differences in the structure of catalytic pockets of various members of this family.     

Vasopressin V2 Receptor/Bioligand Interactions

We predict some essential interactions between the V2 vasopressin renal receptor (V2R) and its agonists [Arg⁸]vasopressin (AVP) and [D-Arg⁸]vasopressin (DAVP), and the non-peptide antagonist OPC-31260. V2R controls antidiuresis and belongs to the superfamily of heptahelical transmembrane (7TM) G-protein-coupled receptors (GPCRs). The receptor was built, the ligands were docked and the structures relaxed using advanced molecular modeling techniques. Docked agonists and antagonists appear to prefer similar V2R compartments. A number of receptor amino acid residues are indicated, mainly in the TM3–TM7 helices, as potentially important in ligand binding. Many of these residues are invariant for either the GPCR superfamily or the subfamily of related (vasopressin V2R, V1aR and V1bR and oxytocin OR) receptors. Moreover, some of the equivalent residues in V1aR have already been found critical for ligand affinity [Mouillac et al., J. Biol. Chem., 270 (1995) 25771].     

Structural studies of the C‐terminal 19‐peptide of serum amyloid A and its Pro→Ala variants interacting with human cystatin C

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.     

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.     

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.     

Structural studies of cysteine proteases and their inhibitors.

Cysteine proteases (CPs) are responsible for many biochemical processes occurring in living organisms and they have been implicated in the development and progression of several diseases that involve abnormal protein turnover. The activity of CPs is regulated among others by their specific inhibitors: cystatins. The main aim of this review is to discuss the structure-activity relationships of cysteine proteases and cystatins, as well as of some synthetic inhibitors of cysteine proteases structurally based on the binding fragments of cystatins.     

Molecular modelling study of the role of cholesterol in the stimulation of the oxytocin receptor.

Cholesterol, an integral component of membranes in Eucaryota, is a modifier of membrane properties. In vivo studies have demonstrated that cholesterol can also modulate activities of some G protein-coupled receptors (GPCRs), which are integral membrane proteins. This can result either from an effect of cholesterol on the membrane fluidity or from specific interactions of the membrane cholesterol with the receptor, as recently demonstrated for the cholecystokinin type beta (CCKRbeta) or the oxytocin receptor (OTR). Using molecular modelling, we studied conformational preferences of cholesterol and several of its analogues. Subsequently, we simulated the distributions of their preferred conformations around the surface of OTR, CCKRbeta and a chimeric oxytocin/cholecystokinin receptor. Consequently, we suggest residues on the surface of OTR which are potentially significant in the OTR/cholesterol interaction.