Functional architecture of vasopressin/oxytocin receptors

Maternal Aggression Driven by the Transient Mobilisation of a Dormant Hormone-Sensitive Circuit

Aggression, a sexually dimorphic behaviour, is prevalent in males and typically absent in virgin females. Following parturition, however, the transient expression of aggression in adult female mice protects pups from predators and infanticide by male conspecifics. While maternal hormones are known to elicit nursing, their potential role in maternal aggression remains elusive. Here, we show in mice that a molecularly defined subset of ventral premammillary (PMvDAT) neurons, instrumental for intermale aggression, switch from quiescence to a hyperexcitable state during lactation. We identify that the maternal hormones prolactin and oxytocin excite these cells through actions that include T-type Ca2+ channels. Optogenetic manipulation or genetic ablation of PMvDAT neurons profoundly affects maternal aggression, while activation of these neurons impairs the expression of non-aggression-related maternal behaviours. This work identifies a monomorphic neural substrate that can incorporate hormonal cues to enable the transient expression of a dormant behavioural program in lactating females.

[Oxytocin and its receptor: molecular and therapeutic approaches]

It is known since the fifties that oxytocin is a neurohormone synthesized in the brain and released in blood circulation to trigger uterus contraction during delivery. It is also involved in milk ejection during breast-feeding. Over the past 25 years, many other central and peripheral functions have been discovered, in particular for attachment between child and parents as well as between individuals and interaction between a human being and its social group. Over this period, we have studied the functional supramolecular architecture of the hormone bound to its receptor. This information was used to design pharmacological probes and drug candidates. This led to the discovery of the first non-peptide oxytocin receptor full agonist. This molecule, LIT-001, restores social interaction in an animal model of autism and paves the way for a treatment of this neurodevelopmental disorder.

Cryo-electron microscopy structure of the antidiuretic hormone arginine-vasopressin V2 receptor signaling complex

The antidiuretic hormone arginine-vasopressin (AVP) forms a signaling complex with the V2 receptor (V2R) and the G_s protein, promoting kidney water reabsorption. Molecular mechanisms underlying activation of this critical G protein-coupled receptor (GPCR) signaling system are still unknown. To fill this gap of knowledge, we report here the cryo-electron microscopy structure of the AVP-V2R-G_s complex. Single-particle analysis revealed the presence of three different states. The two best maps were combined with computational and nuclear magnetic resonance spectroscopy constraints to reconstruct two structures of the ternary complex. These structures differ in AVP and G_s binding modes. They reveal an original receptor-G_s interface in which the Gα_s subunit penetrates deep into the active V2R. The structures help to explain how V2R R137H or R137L/C variants can lead to two severe genetic diseases. Our study provides important structural insights into the function of this clinically relevant GPCR signaling complex.

Elucidating Solution Structures of Cyclic Peptides Using Molecular Dynamics Simulations

Protein-protein interactions are vital to biological processes, but the shape and size of their interfaces make them hard to target using small molecules. Cyclic peptides have shown promise as protein-protein interaction modulators, as they can bind protein surfaces with high affinity and specificity. Dozens of cyclic peptides are already FDA approved, and many more are in various stages of development as immunosuppressants, antibiotics, antivirals, or anticancer drugs. However, most cyclic peptide drugs so far have been natural products or derivatives thereof, with de novo design having proven challenging. A key obstacle is structural characterization: cyclic peptides frequently adopt multiple conformations in solution, which are difficult to resolve using techniques like NMR spectroscopy. The lack of solution structural information prevents a thorough understanding of cyclic peptides' sequence-structure-function relationship. Here we review recent development and application of molecular dynamics simulations with enhanced sampling to studying the solution structures of cyclic peptides. We describe novel computational methods capable of sampling cyclic peptides' conformational space and provide examples of computational studies that relate peptides' sequence and structure to biological activity. We demonstrate that molecular dynamics simulations have grown from an explanatory technique to a full-fledged tool for systematic studies at the forefront of cyclic peptide therapeutic design.

Genes, hormones, and circuits: an integrative approach to study the evolution of social behavior

Tremendous progress has been made in our understanding of the ultimate and proximate mechanisms underlying social behavior, yet an integrative evolutionary analysis of its underpinnings has been difficult. In this review, we propose that modern genomic approaches can facilitate such studies by integrating four approaches to brain and behavior studies: (1) animals face many challenges and opportunities that are ecologically and socially equivalent across species; (2) they respond with species-specific, yet quantifiable and comparable approach and avoidance behaviors; (3) these behaviors in turn are regulated by gene modules and neurochemical codes; and (4) these behaviors are governed by brain circuits such as the mesolimbic reward system and the social behavior network. For each approach, we discuss genomic and other studies that have shed light on various aspects of social behavior and its underpinnings and suggest promising avenues for future research into the evolution of neuroethological systems.

The impact of GPCR structures on pharmacology and structure-based drug design

After many years of effort, recent technical breakthroughs have enabled the X-ray crystal structures of three G-protein-coupled receptors (GPCRs) (beta1 and beta2 adrenergic and adenosine A(2a)) to be solved in addition to rhodopsin. GPCRs, like other membrane proteins, have lagged behind soluble drug targets such as kinases and proteases in the number of structures available and the level of understanding of these targets and their interaction with drugs. The availability of increasing numbers of structures of GPCRs is set to greatly increase our understanding of some of the key issues in GPCR biology. In particular, what constitutes the different receptor conformations that are involved in signalling and the molecular changes which occur upon receptor activation. How future GPCR structures might alter our views on areas such as agonist-directed signalling and allosteric regulation as well as dimerization is discussed. Knowledge of crystal structures in complex with small molecules will enable techniques in drug discovery and design, which have previously only been applied to soluble targets, to now be used for GPCR targets. These methods include structure-based drug design, virtual screening and fragment screening. This review considers how these methods have been used to address problems in drug discovery for kinase and protease targets and therefore how such methods are likely to impact GPCR drug discovery in the future.

VRQ397 (CRAVKY): a novel noncompetitive V2 receptor antagonist

Vasopressin type 2 receptor (V2R) exhibits mostly important properties for hydroosmotic equilibrium and, to a lesser extent, on vasomotricity. Drugs currently acting on this receptor are analogs of the natural neuropeptide, arginine vasopressin (AVP), and hence are competitive ligands. Peptides that reproduce specific sequences of a given receptor have lately been reported to interfere with its action, and if such molecules arise from regions remote from the binding site they would be anticipated to exhibit noncompetitive antagonism, but this has yet to be shown for V2R. Six peptides reproducing juxtamembranous regions of V2R were designed and screened; the most effective peptide, cravky (labeled VRQ397), was characterized. VRQ397 was potent (IC(50) = 0.69 +/- 0.25 nM) and fully effective in inhibiting V2R-dependent physiological function, specifically desmopressin-L-desamino-8-arginine-vasopressin (DDAVP)-induced cremasteric vasorelaxation; this physiological functional assay was utilized to avoid overlooking interference of specific signaling events. A dose-response profile revealed a noncompetitive property of VRQ397; correspondingly, VRQ397 bound specifically to V2R-expressing cells could not displace its natural ligand, AVP, but modulated AVP binding kinetics (dissociation rate). Specificity of VRQ397 was further confirmed by its inability to bind to homologous V1 and oxytocin receptors and its inefficacy to alter responses to stimulation of these receptors. VRQ397 exhibited pharmacological permissiveness on V2R-induced signals, as it inhibited DDAVP-induced PGI(2) generation but not that of cAMP or recruitment of beta-arrestin2. Consistent with in vitro and ex vivo effects as a V2R antagonist, VRQ397 displayed anticipated in vivo aquaretic efficacy. We hereby describe the discovery of a first potent noncompetitive antagonist of V2R, which exhibits functional selectivity, in line with properties of a negative allosteric modulator.

A vasopressin/oxytocin-related conopeptide with gamma-carboxyglutamate at position 8

Vasopressins and oxytocins are homologous, ubiquitous and multifunctional peptides present in animals. Conopressins are vasopressin/oxytocin-related peptides that have been found in the venom of cone snails, a genus of marine predatory molluscs that envenom their prey with a complex mixture of neuroactive peptides. In the present paper, we report the purification and characterization of a unique conopressin isolated from the venom of Conus villepinii, a vermivorous cone snail species from the western Atlantic Ocean. This novel peptide, designated gamma-conopressin-vil, has the sequence CLIQDCPgammaG* (gamma is gamma-carboxyglutamate and * is C-terminal amidation). The unique feature of this vasopressin/oxytocin-like peptide is that the eighth residue is gamma-carboxyglutamate instead of a neutral or basic residue; therefore it could not be directly classified into either the vasopressin or the oxytocin peptide families. Nano-NMR spectroscopy of the peptide isolated directly from the cone snails revealed that the native gamma-conopressin-vil undergoes structural changes in the presence of calcium. This suggests that the peptide binds calcium, and the calcium-binding process is mediated by the gamma-carboxyglutamate residue. However, the negatively charged residues in the sequence of gamma-conopressin-vil may mediate calcium binding by a novel mechanism not observed in other peptides of this family.

Systematic analysis of the entire second extracellular loop of the V(1a) vasopressin receptor: key residues, conserved throughout a G-protein-coupled receptor family, identified

The roles of extracellular residues of G-protein-coupled receptors (GPCRs) are not well defined compared with residues in transmembrane helices. Nevertheless, it has been established that extracellular domains of both peptide-GPCRs and amine-GPCRs incorporate functionally important residues. Extracellular loop 2 (ECL2) has attracted particular interest, because the x-ray structure of bovine rhodopsin revealed that ECL2 projects into the binding crevice within the transmembrane bundle. Our study provides the first comprehensive investigation into the role of the individual residues comprising the entire ECL2 domain of a small peptide-GPCR. Using the V(1a) vasopressin receptor, systematic substitution of all of the ECL2 residues by Ala generated 30 mutant receptors that were characterized pharmacologically. The majority of these mutant receptor constructs (24 in total) had essentially wild-type ligand binding and intracellular signaling characteristics, indicating that these residues are not critical for normal receptor function. However, four aromatic residues Phe(189), Trp(206), Phe(209), and Tyr(218) are important for agonist binding and receptor activation and are highly conserved throughout the neurohypophysial hormone subfamily of peptide-GPCRs. Located in the middle of ECL2, juxtaposed to the highly conserved disulfide bond, Trp(206) and Phe(209) project into the binding crevice. Indeed, Phe(209) is part of the Cys-X-X-X-Ar (where Ar is an aromatic residue) motif, which is well conserved in both peptide-GPCRs and amine-GPCRs. In contrast, Phe(189) and Tyr(218), located at the extreme ends of ECL2, may be important for determining the position of the ECL2 cap over the binding crevice. This study provides mechanistic insight into the roles of highly conserved ECL2 residues.

Charged extracellular residues, conserved throughout a G-protein-coupled receptor family, are required for ligand binding, receptor activation, and cell-surface expression

For G-protein-coupled receptors (GPCRs) in general, the roles of extracellular residues are not well defined compared with residues in transmembrane helices (TMs). Nevertheless, extracellular residues are important for various functions in both peptide-GPCRs and amine-GPCRs. In this study, the V(1a) vasopressin receptor was used to systematically investigate the role of extracellular charged residues that are highly conserved throughout a subfamily of peptide-GPCRs, using a combination of mutagenesis and molecular modeling. Of the 13 conserved charged residues identified in the extracellular loops (ECLs), Arg(116) (ECL1), Arg(125) (top of TMIII), and Asp(204) (ECL2) are important for agonist binding and/or receptor activation. Molecular modeling revealed that Arg(125) (and Lys(125)) stabilizes TMIII by interacting with lipid head groups. Charge reversal (Asp(125)) caused re-ordering of the lipids, altered helical packing, and increased solvent penetration of the TM bundle. Interestingly, a negative charge is excluded at this locus in peptide-GPCRs, whereas a positive charge is excluded in amine-GPCRs. This contrasting conserved charge may reflect differences in GPCR binding modes between peptides and amines, with amines needing to access a binding site crevice within the receptor TM bundle, whereas the binding site of peptide-GPCRs includes more extracellular domains. A conserved negative charge at residue 204 (ECL2), juxtaposed to the highly conserved disulfide bond, was essential for agonist binding and signaling. Asp(204) (and Glu(204)) establishes TMIII contacts required for maintaining the beta-hairpin fold of ECL2, which if broken (Ala(204) or Arg(204)) resulted in ECL2 unfolding and receptor dysfunction. This study provides mechanistic insight into the roles of conserved extracellular residues.

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.

A chemogenomic analysis of the transmembrane binding cavity of human G-protein-coupled receptors

The amino acid sequences of 369 human nonolfactory G-protein-coupled receptors (GPCRs) have been aligned at the seven transmembrane domain (TM) and used to extract the nature of 30 critical residues supposed--from the X-ray structure of bovine rhodopsin bound to retinal--to line the TM binding cavity of ground-state receptors. Interestingly, the clustering of human GPCRs from these 30 residues mirrors the recently described phylogenetic tree of full-sequence human GPCRs (Fredriksson et al., Mol Pharmacol 2003;63:1256-1272) with few exceptions. A TM cavity could be found for all investigated GPCRs with physicochemical properties matching that of their cognate ligands. The current approach allows a very fast comparison of most human GPCRs from the focused perspective of the predicted TM cavity and permits to easily detect key residues that drive ligand selectivity or promiscuity.

The N-Terminal Juxtamembrane Segment of the V1a Vasopressin Receptor Provides Two Independent Epitopes Required for High-Affinity Agonist Binding and Signaling

It is fundamentally important to define how agonist-receptor interaction differs from antagonist-receptor interaction. The V1a vasopressin receptor (V1aR) is a member of the neurohypophysial hormone subfamily of G protein-coupled receptors. Using alanine-scanning mutagenesis of the N-terminal juxtamembrane segment of the V1aR, we now establish that Glu54 (1.35) is critical for arginine vasopressin binding. The mutant [E54A]V1aR exhibited decreased arginine vasopressin affinity (1700-fold) and disrupted signaling, but antagonist binding was unaffected. Mutation of Glu54 had an almost identical pharmacological effect as mutation of Arg46, raising the possibility that agonist binding required a mutual interaction between Glu54 and Arg46. The role of these two charged residues was investigated by 1) substituting Glu54; 2) inserting additional Glu/Arg in transmembrane helix (TM) 1; 3) repositioning the Glu/Arg in TM1; and 4) characterizing the reciprocal mutant [R46E/E54R]V1aR. We conclude that 1) the positive/negative charges need to be precisely positioned in this N terminus/TM1 segment; and 2) Glu54 and Arg46 function independently, providing two discrete epitopes required for high-affinity agonist binding and signaling. This study explains why Glu and Arg, part of an -R(X3)L/V(X3)E(X3)L- motif, are conserved at these loci throughout this G protein-coupled receptor subfamily and provides molecular insight into key differences between agonist and antagonist binding requirements.

The 3D Structure of the Binding Pocket of the Human Oxytocin Receptor for Benzoxazine Antagonists, Determined by Molecular Docking, Scoring Functions and 3D-QSAR Methods

Molecular docking and 3D-QSAR studies were performed to determine the binding mode for a series of benzoxazine oxytocin antagonists taken from the literature. Structural hypotheses were generated by docking the most active molecule to the rigid receptor by means of AutoDock 3.05. The cluster analysis yielded seven possible binding conformations. These structures were refined by using constrained simulated annealing, and the further ligands were aligned in the refined receptor by molecular docking. A good correlation was found between the estimated deltaG(bind) and the pKi values for complex F. The Connolly-surface analysis, CoMFA and CoMSIA models q2(CoMFA) = 0.653, q2(CoMSA) = 0.630 and r2(pred,CoMFA) = 0.852 , r2(pred,CoMSIA) = 0.815) confirmed the scoring function results. The structural features of the receptor-ligand complex and the CoMFA and CoMSIA fields are in closely connected. These results suggest that receptor-ligand complex F is the most likely binding hypothesis for the studied benzoxazine analogs.

An arginyl in the N-terminus of the V1a vasopressin receptor is part of the conformational switch controlling activation by agonist

Defining how the agonist-receptor interaction differs from that of the antagonist-receptor and understanding the mechanisms of receptor activation are fundamental issues in cell signalling. The V1a vasopressin receptor (V1aR) is a member of a family of related G-protein coupled receptors that are activated by neurohypophysial peptide hormones, including vasopressin (AVP). It has recently been reported that an arginyl in the distal N-terminus of the V1aR is critical for binding agonists but not antagonists. To determine specific features required at this locus to support high affinity agonist binding and second messenger generation, Arg46 was substituted by all other 19 encoded amino acids. Our data establish that there is an absolute requirement for arginyl, as none of the [R46X]V1aR mutant constructs supported high affinity agonist binding and all 19 had defective signalling. In contrast, all of the mutant receptors possessed wildtype binding for both peptide and nonpeptide antagonists. The ratio of Ki to EC50, an indicator of efficacy, was increased for all substitutions. Consequently, although [R46X]V1aR constructs have a lower affinity for agonist, once AVP has bound all 19 are more likely than the wildtype V1aR to become activated. Therefore, in the wildtype V1aR, Arg46 constrains the inactive conformation of the receptor. On binding AVP this constraint is alleviated, promoting the transition to active V1aR. Our findings explain why arginyl is conserved at this locus throughout the evolutionary lineage of the neurohypophysial peptide hormone receptor family of G-protein coupled receptors.

Identification of the binding sites of the SR49059 nonpeptide antagonist into the V1a vasopressin receptor using sulfydryl-reactive ligands and cysteine mutants as chemical sensors

To identify the binding site of the human V1a vasopressin receptor for the selective nonpeptide antagonist SR49059, we have developed a site-directed irreversible labeling strategy that combines mutagenesis of the receptor and use of sulfydryl-reactive ligands. Based on a three-dimensional model of the antagonist docked into the receptor, hypothetical ligand-receptor interactions were investigated by replacing the residues potentially involved in the binding of the antagonist into cysteines and designing analogues of SR49059 derivatized with isothiocyanate or alpha-chloroacetamide moieties. The F225C, F308C, and K128C mutants of the V1a receptor were expressed in COS-7 or Chinese hamster ovary cells, and their pharmacological properties toward SR49059 and its sulfydryl-reactive analogues were analyzed. We demonstrated that treatment of the F225C mutant with the isothiocyanate-derivative compound led to dose-dependent inhibition of the residual binding of the radio-labeled antagonist [125I]HO-LVA. This inhibition is probably the consequence of a covalent irreversible chemical modification, which is only possible when close contacts and optimal orientations exist between reactive groups created both on the ligand and the receptor. This result validated the three-dimensional model hypothesis. Thus, we propose that residue Phe225, located in transmembrane domain V, directly participates in the binding of the V1a-selective nonpeptide antagonist SR49059. This conclusion is in complete agreement with all our previous data on the definition of the agonist/antagonist binding to members of the oxytocin/vasopressin receptor family.

Science review: Vasopressin and the cardiovascular system part 1--receptor physiology

Vasopressin is emerging as a rational therapy for vasodilatory shock states. Unlike other vasoconstrictor agents, vasopressin also has vasodilatory properties. The goal of the present review is to explore the vascular actions of vasopressin. In part 1 of the review we discuss structure, signaling pathways, and tissue distributions of the classic vasopressin receptors, namely V₁ vascular, V₂ renal, V₃ pituitary and oxytocin receptors, and the P₂ class of purinoreceptors. Knowledge of the function and distribution of vasopressin receptors is key to understanding the seemingly contradictory actions of vasopressin on the vascular system. In part 2 of the review we discuss the effects of vasopressin on vascular smooth muscle and the heart, and we summarize clinical studies of vasopressin in shock states.

Vasopressin increases milk flow and milk fat concentration in the goat

AIM

The aim of this study was to investigate whether physiological doses of vasopressin affect milk flow and milk composition in goats.

METHODS

Six goats were given vasopressin intravenously at doses of 0.1 ng kg(-1) min(-1) (a low dose, VL), 0.3 ng kg(-1) min(-1) (a medium dose, VM) and 1.0 ng kg(-1) min(-1) (a high dose, VH), oxytocin at 10 ng kg(-1) min(-1) and 0.15 m NaCl during 90 min at random.

RESULTS

Milk flow increased from 1.3 +/- 0.2 to 3.5 +/- 1.0 g min(-1) (VL), from 1.8 +/- 0.4 to 4.3 +/- 0.6 g min(-1) (VM) and from 1.5 +/- 0.2 to 7.9 +/- 0.9 g min(-1) (VH) during the first 30-min infusions (all P < 0.001), respectively, and then declined. The response to oxytocin corresponded to that of VH. Milk fat concentration increased from 2.1 +/- 0.1 to 3.3 +/- 0.4% during VH and from 2.8 +/- 0.6 to 4.1 +/- 0.6% during oxytocin (both P < 0.001), and tended to increase during VM and VL. The mammary blood flow of four goats was studied with ultrasound during VH and oxytocin infusion. Milk flow increased from 1.8 +/- 0.4 to 4.7 +/- 1.0 g min(-1) during VH and from 1.9 +/- 0.3 to 12.4 +/- 2.4 g min(-1) during oxytocin (both P < 0.001). Mammary blood flow decreased from 262 +/- 37 to 53 +/- 7 mL min(-1) (P < 0.001) during VH, but did not change during oxytocin.

CONCLUSION

We conclude that vasopressin increases milk flow and milk fat concentration and suggest that vasopressin, like oxytocin, acts on the myoepithelial cells.

Protein-based virtual screening of chemical databases. II. Are homology models of G-Protein Coupled Receptors suitable targets?

The aim of the current study is to investigate whether homology models of G-Protein-Coupled Receptors (GPCRs) that are based on bovine rhodopsin are reliable enough to be used for virtual screening of chemical databases. Starting from the recently described 2.8 A-resolution X-ray structure of bovine rhodopsin, homology models of an "antagonist-bound" form of three human GPCRs (dopamine D3 receptor, muscarinic M1 receptor, vasopressin V1a receptor) were constructed. The homology models were used to screen three-dimensional databases using three different docking programs (Dock, FlexX, Gold) in combination with seven scoring functions (ChemScore, Dock, FlexX, Fresno, Gold, Pmf, Score). Rhodopsin-based homology models turned out to be suitable, indeed, for virtual screening since known antagonists seeded in the test databases could be distinguished from randomly chosen molecules. However, such models are not accurate enough for retrieving known agonists. To generate receptor models better suited for agonist screening, we developed a new knowledge- and pharmacophore-based modeling procedure that might partly simulate the conformational changes occurring in the active site during receptor activation. Receptor coordinates generated by this new procedure are now suitable for agonist screening. We thus propose two alternative strategies for the virtual screening of GPCR ligands, relying on a different set of receptor coordinates (antagonist-bound and agonist-bound states).

Docking ligands to vasopressin and oxytocin receptors via genetic algorithm

The aim of the study was to computer-dock selected ligands to neurophyseal receptors in order to identify amino acid residues responsible for ligand-receptor interactions. To this aim, reliable oxytocin receptor (OTR) and arginine-vasopressin receptor (V1aR/V2R) models were built. The OTR-selective agonist [Thr4,Gly7]OT, the OTR-selective cyclohexapeptide antagonist L-366,948 and OT itself were docked via genetic algorithm to OTR, V1aR, and V2R and relaxed using a constrained simulated annealing protocol. For the analysis of receptor/ligand interactions a subset of initial conformations was chosen using energetic and steric criteria. All three ligands seem to prefer similar modes of binding to the receptors, manifested by repetitive residues of the receptors which directly interact with the ligands. Taking into account that many aspects of mechanisms of G protein-coupled receptor (GPCR) action are still unsolved, the results obtained with the docking simulations may propose future experimental research, especially in site-directed mutagenesis analysis and searching for key amino acid residues responsible for drug activities.

Molecular modeling of interactions of the non-peptide antagonist YM087 with the human vasopressin V1a, V2 receptors and with oxytocin receptors

The nonapeptide hormones arginine vasopressin (CYFQNCPRG-NH2, AVP) and oxytocin (CYIQNCPLG-NH2, OT), control many essential functions in mammals. Their main activities include the urine concentration (via stimulation of AVP V2 receptors, V2R, in the kidneys), blood pressure regulation (via stimulation of vascular V1a AVP receptors, V1aR), ACTH control (via stimulation of V1b receptors, V1bR, in the pituitary) and labor and lactation control (via stimulation of OT receptors, OTR, in the uterus and nipples, respectively). All four receptor subtypes belong to the GTP-binding (G) protein-coupled receptor (GPCR) family. This work consists of docking of YM087, a potent non-peptide V1aR and V2R - but not OTR - antagonist, into the receptor models based on relatively new theoretical templates of rhodopsin (RD) and opiate receptors, proposed by Mosberg et al. (Univ. of Michigan, Ann Arbor, USA). It is simultaneously demonstrated that this RD template satisfactorily compares with the first historical GPCR structure of bovine rhodopsin (Palczewski et al., 2000) and that homology-modeling of V2R, V1aR and OTR using opiate receptors as templates is rational, based on relatively high (20-60%) sequence homology among the set of 4 neurophyseal and 4 opiate receptors. YM087 was computer-docked to V1aR, V2R and OTR using the AutoDock (Olson et al., Scripps Research Institute, La Jolla, USA) and subsequently relaxed using restrained simulated annealing and molecular dynamics, as implemented in AMBER program (Kollman et al., University of California, San Francisco, USA). From about 80 diverse configurations, sampled for each of the three ligand/receptor systems, 3 best energy-relaxed complexes were selected for mutual comparisons. Similar docking modes were found for the YM087/V1aR and YM087/V2R complexes, diverse from those of the YM087/OTR complexes, in agreement with the molecular affinity data.