ABSTRACT
We have developed two ligand- and receptor-based computational approaches to study the physicochemical properties relevant to the biological activity of vasopressin V2 receptor (V2R) antagonist and eventually to predict the expected binding mode to V2R. The obtained quantitative structure activity relationship (QSAR) model showed a correlation of the antagonist activity with the hydration energy (EH2O), the polarizability (P), and the calculated partial charge on atom N7 (q6) of the common substructure. The first two descriptors showed a positive contribution to antagonist activity, while the third one had a negative contribution. V2R was modeled and further relaxed on a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocoline (POPC) membrane by molecular dynamics simulations. The receptor antagonist complexes were guessed by molecular docking, and the stability of the most relevant structures was also evaluated by molecular dynamics simulations. As a result, amino acid residues Q96, W99, F105, K116, F178, A194, F307, and M311 were identified with the probably most relevant antagonist-receptor interactions on the studied complexes. The proposed QSAR model could explain the molecular properties relevant to the antagonist activity. The contributions to the antagonist-receptor interaction appeared also in agreement with the binding mode of the complexes obtained by molecular docking and molecular dynamics. These models will be used in further studies to look for new V2R potential antagonist molecules.
Subject(s)
Antidiuretic Hormone Receptor Antagonists/chemistry , Chemical Phenomena , Models, Molecular , Receptors, Vasopressin/chemistry , Algorithms , Amino Acid Sequence , Antidiuretic Hormone Receptor Antagonists/pharmacology , Binding Sites , Cluster Analysis , Ligands , Molecular Conformation , Molecular Docking Simulation , Molecular Dynamics Simulation , Molecular Structure , Protein Binding , Quantitative Structure-Activity Relationship , Structure-Activity RelationshipABSTRACT
Ghrelin is a peptide hormone involved in multiple functions, including growth hormone release stimulation, food intake regulation, and metabolic and cytoprotective effect. A novel family of peptides with internal cycles was designed as ghrelin analogs and the biological activity of two of them (A228 and A233) was experimentally studied in-depth. In this work, an in silico strategy was developed for describing and assessing the binding modes of A228 and A233 to GHS-R1a (ghrelin receptor) comparing it with ghrelin and GHRP-6 peptides. Several reported structures of different G protein coupled receptors were used as templates, to obtain a good quality model of GHS-R1a. The best model was selected by preliminary molecular docking with ghrelin and GHRP-6. Docking was used to estimate peptide orientations in the binding site of the best model, observing a superposition of its N-terminal and its first aromatic residue. To test the complex stability in time, the C-terminal fragments of each peptide were added and the complexes were inserted a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) membrane, performing a molecular dynamic simulation for 100 ns using the CHARMM36 force field. Despite of the structural differences, the studied peptides share a common binding mode; the N-terminal interacts with E124 and the aromatic residue close to it, with the aromatic cluster (F279, F309, and F312). A preliminary pharmacophore model, consisting in a positive charged amine and an aromatic ring at an approximate distance of 0.79 nm, can be proposed. The results here described could represent a step forward in the efficient search of new ghrelin analogs.
