Modeling of human prokineticin receptors: interactions with novel small-molecule binders and potential off-target drugs

Therapeutic Potential of Targeting Prokineticin Receptors in Diseases

The prokineticins (PKs) were discovered approximately 20 years ago as small peptides inducing gut contractility. Today, they are established as angiogenic, anorectic, and proinflammatory cytokines, chemokines, hormones, and neuropeptides involved in variety of physiologic and pathophysiological pathways. Their altered expression or mutations implicated in several diseases make them a potential biomarker. Their G-protein coupled receptors, PKR1 and PKR2, have divergent roles that can be therapeutic target for treatment of cardiovascular, metabolic, and neural diseases as well as pain and cancer. This article reviews and summarizes our current knowledge of PK family functions from development of heart and brain to regulation of homeostasis in health and diseases. Finally, the review summarizes the established roles of the endogenous peptides, synthetic peptides and the selective ligands of PKR1 and PKR2, and nonpeptide orthostatic and allosteric modulator of the receptors in preclinical disease models. The present review emphasizes the ambiguous aspects and gaps in our knowledge of functions of PKR ligands and elucidates future perspectives for PK research. SIGNIFICANCE STATEMENT: This review provides an in-depth view of the prokineticin family and PK receptors that can be active without their endogenous ligand and exhibits "constitutive" activity in diseases. Their non- peptide ligands display promising effects in several preclinical disease models. PKs can be the diagnostic biomarker of several diseases. A thorough understanding of the role of prokineticin family and their receptor types in health and diseases is critical to develop novel therapeutic strategies with safety concerns.

Non-Peptide Agonists and Antagonists of the Prokineticin Receptors

The prokineticin family comprises a group of secreted peptides that can be classified as chemokines based on their structural features and chemotactic and immunomodulatory functions. Prokineticins (PKs) bind with high affinity to two G protein-coupled receptors (GPCRs). Prokineticin receptor 1 (PKR1) and prokineticin receptor 2 (PKR2) are involved in a variety of physiological functions such as angiogenesis and neurogenesis, hematopoiesis, the control of hypothalamic hormone secretion, the regulation of circadian rhythm and the modulation of complex behaviors such as feeding and drinking. Dysregulation of the system leads to an inflammatory process that is the substrate for many pathological conditions such as cancer, pain, neuroinflammation and neurodegenerative diseases such as Alzheimer's and Parkinson's disease. The use of PKR's antagonists reduces PK2/PKRs upregulation triggered by various inflammatory processes, suggesting that a pharmacological blockade of PKRs may be a successful strategy to treat inflammatory/neuroinflammatory diseases, at least in rodents. Under certain circumstances, the PK system exhibits protective/neuroprotective effects, so PKR agonists have also been developed to modulate the prokineticin system.

Anosmin 1 N-terminal domains modulate prokineticin receptor 2 activation by prokineticin 2

The X-linked form of Kallmann syndrome (KS), characterized by hypogonadotropic hypogonadism and anosmia, is due to mutations in the ANOS1 gene that encodes for the extracellular matrix (ECM) protein anosmin 1. Prokineticins (PKs) exert their biological functions through the activation of the G protein-coupled receptors (GPCRs) prokineticin receptor 1 and 2 (PKR1, 2), and mutations in the PK2 and PKR2 genes are involved in the pathogenesis of KS. We have previously shown interaction between PKR2 and anosmin 1 in vitro. In the current report we present evidence of the modulation of PK2/PKR2 activity by anosmin 1, since this protein is able to enhance the activation of the ERK1/2 (extracellular signal-regulated kinase 1/2) pathway elicited by PK2 through PKR2. We also show that the N-terminal region of anosmin 1, capable of binding to the PK2-binding domain of PKR2, seems to be responsible for this effect. The whey acidic protein domain (WAP) is necessary for this modulatory activity, although data from GST pull-down (glutathione-S-transferase) and analysis of the N267K mutation in the fibronectin type III domain 1 (FnIII.1) suggest the cysteine-rich (CR) and the FnIII.1 domains could assist the WAP domain both in the binding to PKR2 and in the modulation of the activation of the receptor by PK2. Our data support the idea of a modulatory role of anosmin 1 in the biological effects controlled by the PK2/PKR2 system.

Structural and Functional Diversity of Animal Toxins Interacting With GPCRs

Peptide toxins from venoms have undergone a long evolutionary process allowing host defense or prey capture and making them highly selective and potent for their target. This has resulted in the emergence of a large panel of toxins from a wide diversity of species, with varied structures and multiple associated biological functions. In this way, animal toxins constitute an inexhaustible reservoir of druggable molecules due to their interesting pharmacological properties. One of the most interesting classes of therapeutic targets is the G-protein coupled receptors (GPCRs). GPCRs represent the largest family of membrane receptors in mammals with approximately 800 different members. They are involved in almost all biological functions and are the target of almost 30% of drugs currently on the market. Given the interest of GPCRs in the therapeutic field, the study of toxins that can interact with and modulate their activity with the purpose of drug development is of particular importance. The present review focuses on toxins targeting GPCRs, including peptide-interacting receptors or aminergic receptors, with a particular focus on structural aspects and, when relevant, on potential medical applications. The toxins described here exhibit a great diversity in size, from 10 to 80 amino acids long, in disulfide bridges, from none to five, and belong to a large panel of structural scaffolds. Particular toxin structures developed here include inhibitory cystine knot (ICK), three-finger fold, and Kunitz-type toxins. We summarize current knowledge on the structural and functional diversity of toxins interacting with GPCRs, concerning first the agonist-mimicking toxins that act as endogenous agonists targeting the corresponding receptor, and second the toxins that differ structurally from natural agonists and which display agonist, antagonist, or allosteric properties.

Identification and Characterization of a New Splicing Variant of Prokineticin 2

Prokineticin 2 (PROK2) is a secreted bioactive peptide that regulates a variety of biological responses via two GPCRs, the prokineticin receptors (PROKRs). The aim of this study was to characterize a new alternatively spliced product of the prok2 gene consisting of four exons. The 40-amino acid peptide, designated PROK2C, is encoded by exon 1 and exon 4, and its expression was detected in the hippocampus and spinal cord of mice. PROK2C was expressed in a heterologous system, Pichia pastoris, and its binding specificity to the amino-terminal regions of PROKR1 and PROKR2 was investigated by GST pull-down experiments. In addition, the introduction of the unnatural amino acid p-benzoyl-L-phenylalanine using amber codon suppression technology demonstrated the role of tryptophan at position 212 of PROKR2 for PROK2C binding by photoactivatable cross-linking. The functional significance of this new isoform was determined in vivo by nociceptive experiments, which showed that PROK2C elicits strong sensitization of peripheral nociceptors to painful stimuli. In order to analyze the induction of PROK2C signal transduction, STAT3 and ERK phosphorylation levels were determined in mammalian CHO cells expressing PROKR1 and PROKR2. Our data show by in vivo and in vitro experiments that PROK2C can bind and activate both prokineticin receptors.

Prokineticin-Receptor Network: Mechanisms of Regulation

Prokineticins are a new class of chemokine-like peptides that bind their G protein-coupled receptors, PKR1 and PKR2, and promote chemotaxis and the production of pro-inflammatory cytokines following tissue injury or infection. This review summarizes the major cellular and biochemical mechanisms of prokineticins pathway regulation that, like other chemokines, include: genetic polymorphisms; mRNA splice modulation; expression regulation at transcriptional and post-transcriptional levels; prokineticins interactions with cell-surface glycosaminoglycans; PKRs degradation, localization, post-translational modifications and oligomerization; alternative signaling responses; binding to pharmacological inhibitors. Understanding these mechanisms, which together exert substantial biochemical control and greatly enhance the complexity of the prokineticin-receptor network, leads to novel opportunities for therapeutic intervention. In this way, besides targeting prokineticins or their receptors directly, it could be possible to indirectly influence their activity by modulating their expression and localization or blocking the downstream signaling pathways.

Analysis of role of aromatic residues in extracellular loop 2 of Prokineticin receptor 2 in ligand binding probed with genetically encoded photo-crosslinkers

Prokineticin 2 (PK2) and Prokineticin 2 beta (PK2β), products of alternative splicing of pk2 gene, are chemokine-like proteins. While PK2 mediates its biological activities by signaling with the same efficiency through two homologous G protein coupled receptors, prokineticin receptor 1 (PKR1) and prokineticin receptor 2 (PKR2), PK2β is able to bind specifically PKR1. Extracellular loop 2 (ECL2) of chemokine receptors is a part of a transmembrane (TM) ligand binding site. In the ECL2 of PKR2 is present, as well as in all chemokine receptors, an aromatic residue cluster, involving tryptophan 212 localized four residues after an ECL2 conserved cysteine, and Phenylalanine 198 located in the top of TM 4. In this work, the photoactivatable unnatural amino acid p-benzoyl-L-phenylalanine is incorporated by amber codon suppression technology into PKR2 in position 212. Experiments of photoactivatable cross-linking demonstrated the role of tryptophan in position 212 for binding the ligand contacting Tryptophan in position 24. We also analyzed the role of Phenylalanine 198 in the specificity of PKRs binding. The comparison of TM-bundle binding sites between PKR1 and PKR2 revealed that they are completely conserved except for one residue: valine 207 in human PKR1, which is phenylalanine 198 in human PKR2. The F198V mutation in PKR2 permits to obtain a receptor able to bind more efficiently PK2β, a ligand highly specific for PKR1.

Discovery of New Phosphoinositide 3-kinase Delta (PI3Kδ) Inhibitors via Virtual Screening using Crystallography-derived Pharmacophore Modelling and QSAR Analysis

BACKGROUND

PI3Kδ is predominantly expressed in hematopoietic cells and participates in the activation of leukocytes. PI3Kδ inhibition is a promising approach for treating inflammatory diseases and leukocyte malignancies. Accordingly, we decided to model PI3Kδ binding.

METHODS

Seventeen PI3Kδ crystallographic complexes were used to extract 94 pharmacophore models. QSAR modelling was subsequently used to select the superior pharmacophore(s) that best explain bioactivity variation within a list of 79 diverse inhibitors (i.e., upon combination with other physicochemical descriptors).

RESULTS

The best QSAR model (r2 = 0.71, r2 LOO = 0.70, r2 press against external testing list of 15 compounds = 0.80) included a single crystallographic pharmacophore of optimal explanatory qualities. The resulting pharmacophore and QSAR model were used to screen the National Cancer Institute (NCI) database for new PI3Kδ inhibitors. Two hits showed low micromolar IC50 values.

CONCLUSION

Crystallography-based pharmacophores were successfully combined with QSAR analysis for the identification of novel PI3Kδ inhibitors.

Designing a multi-epitope based vaccine to combat Kaposi Sarcoma utilizing immunoinformatics approach

Kaposi’s sarcoma-associated herpesvirus (KSHV) responsible for causing Kaposi sarcoma (KS), an opportunistic angioproliferative neoplasm is emerging rapidly. Despite this there is no permanent cure for this disease. The present study was aimed to design a multi-epitope based vaccine targeting the major glycoproteins of KSHV which plays an important role in the virus entry. After the application of rigorous immunoinformatics analysis and several immune filters, the multi-epitope vaccine was constructed, consisting of CD4, CD8 and IFN-γ inducing epitopes. Several physiochemical characteristics, allergenicity and antigenicity of the multi-epitope vaccine were analyzed in order to ensure its safety and immunogenicity. Further, the binding affinity and stability of the vaccine with Toll like receptor -9 (TLR-9) was analyzed by molecular docking and dynamics simulation studies. In addition, an in silico cloning was performed to ensure the expression and translation efficiency of the vaccine, utilizing pET-28a (+) vector. Such T-cell-based immunotherapies which leverage this mechanism could prove their potential against cancer. Further, the authors propose to test the present findings in the lab settings to ensure the safety, immunogenicity and efficacy of the presented vaccine which may help in controlling KSHV infection.

The Prokineticins: Neuromodulators and Mediators of Inflammation and Myeloid Cell-Dependent Angiogenesis

The mammalian prokineticins family comprises two conserved proteins, EG-VEGF/PROK1 and Bv8/PROK2, and their two highly related G protein-coupled receptors, PKR1 and PKR2. This signaling system has been linked to several important biological functions, including gastrointestinal tract motility, regulation of circadian rhythms, neurogenesis, angiogenesis and cancer progression, hematopoiesis, and nociception. Mutations in PKR2 or Bv8/PROK2 have been associated with Kallmann syndrome, a developmental disorder characterized by defective olfactory bulb neurogenesis, impaired development of gonadotropin-releasing hormone neurons, and infertility. Also, Bv8/PROK2 is strongly upregulated in neutrophils and other inflammatory cells in response to granulocyte-colony stimulating factor or other myeloid growth factors and functions as a pronociceptive mediator in inflamed tissues as well as a regulator of myeloid cell-dependent tumor angiogenesis. Bv8/PROK2 has been also implicated in neuropathic pain. Anti-Bv8/PROK2 antibodies or small molecule PKR inhibitors ameliorate pain arising from tissue injury and inhibit angiogenesis and inflammation associated with tumors or some autoimmune disorders.

PK2β ligand, a splice variant of prokineticin 2, is able to modulate and drive signaling through PKR1 receptor

Prokineticin-2 (PK2) is a secreted bioactive peptide that signals through two GPCRs, the prokineticin receptors (PKRs), and regulates a variety of biological processes including angiogenesis, immunity and nociception. The PK2 primary transcript has two alternative splice variants, PK2 and PK2L (a Long form) which is cleaved in an active peptide, named PK2β that preferentially binds to PKR1 receptor. The aim of this study was to characterize the PK2β. Using different Saccharomyces cerevisiae strains, we examined the specificity of PKR1 and PKR2 G-protein coupling following PK2β binding. Data obtained in yeast confirmed that PK2 binds both receptors, inducing a comparable response throughout a promiscuous coupling of G protein subtypes. Conversely, we demonstrated, for the first time, that PK2β preferentially binding to PKR1, activates a signaling cascade that not depends on Gα_i/o coupling. The binding specificity of PK2β for PKR1 was evaluated by the analysis of PKR mutant in yeast and GST pull-down experiments, suggesting an important role of PKR1 amino-terminal region. We also evaluated the ability of PK2β to differentially activate PKR1 and/or PKR2 by in vivo nociceptive experiments and we showed that PK2β induces intense sensitization of peripheral nociceptors to painful stimuli through the activation of PKR1. To analyze PK2β-induced signal transduction, we demonstrated the inability of PK2β to induce STAT3 protein phosphorylation in organotypic primary explants from mice Dorsal Root Ganglion (DRG), an important pain station. The control of the concentration ratio between PK2β and PK2 could be one of the keys to allow the specificity of the cell response of prokineticin signaling pathway.

Discovery of FIXa inhibitors by combination of pharmacophore modeling, molecular docking, and 3D-QSAR modeling

Human Coagulation Factor IXa (FIXa), specifically inhibited at the initiation stage of the blood coagulation cascade, is an excellent target for developing selective and safe anticoagulants. To explore this inhibitory mechanism, 86 FIXa inhibitors were selected to generate pharmacophore models and subsequently SAR models. Both best pharmacophore model and ROC curve were built through the Receptor-Ligand Pharmacophore Generation module. CoMFA model based on molecular docking and PLS factor analysis methods were developed. Model propagations values are q² = 0.709, r² = 0.949, and r²_pred = 0.905. The satisfactory q² value of 0.609, r² value of 0.962, and r²_pred value of 0.819 for CoMSIA indicated that the CoMFA and CoMSIA models are both available to predict the inhibitory activity on FIXa. On the basis of pharmacophore modeling, molecular docking, and 3D-QSAR modeling screening, six molecules are screened as potential FIXa inhibitors.

Identification of novel selective MMP-9 inhibitors as potential anti-metastatic lead using structure-based hierarchical virtual screening and molecular dynamics simulation

MMP-9 is an attractive target for the development of new anticancer drugs. In the current study, pharmacophore modeling was employed using two highly active and selective gelatinase inhibitors obtained from two X-ray crystal structures (PDB IDs: and ) to identify novel selective MMP-9 inhibitors. The derived model was refined manually and also validated by the GH scoring method. The refined pharmacophore model, ADRR, was able to retrieve 86% of actives with a GH score of 0.774, indicating that the model was capable of retrieving the active MMP-9 inhibitors. ADRR was used to screen 2 838 166 unique structures. Hit filtration was carried out using a fitness score >1.5 and drug-likeness properties. Hierarchical clustering generates 33 clusters based on diversity. A total of 33 molecules were obtained and these molecules were taken for cross-docking studies with 5 subtype MMPs. Among 33 tested, 2 molecules, P10A-0000088030 (Lig-1) and P10A-0001383812 (Lig-2), were found to have the highest docking scores (-8.59 kcal mol(-1) and -8.27 kcal mol(-1)) towards MMP-9 compared with the other MMPs. Further MM-GBSA analysis was performed for two hits with 5 subtype MMPs to reveal the essential features that contribute to selectivity. The results showed that van der Waals contributions play a central role in determining the selectivity of MMP-9 inhibitors. Molecular dynamics studies were carried out for total time of 330 ns to assess the stability of ligands at the active site. MD analysis showed that binding of Lig-1 with MMP-9 is stable compared to that with Lig-2. Hence, we suggest the Lig-1 compound as a good lead in designing novel potent inhibitors of MMP-9.

A Novel Drug-Mouse Phenotypic Similarity Method Detects Molecular Determinants of Drug Effects

The molecular mechanisms that translate drug treatment into beneficial and unwanted effects are largely unknown. We present here a novel approach to detect gene-drug and gene-side effect associations based on the phenotypic similarity of drugs and single gene perturbations in mice that account for the polypharmacological property of drugs. We scored the phenotypic similarity of human side effect profiles of 1,667 small molecules and biologicals to profiles of phenotypic traits of 5,384 mouse genes. The benchmarking with known relationships revealed a strong enrichment of physical and indirect drug-target connections, causative drug target-side effect links as well as gene-drug links involved in pharmacogenetic associations among phenotypically similar gene-drug pairs. The validation by in vitro assays and the experimental verification of an unknown connection between oxandrolone and prokineticin receptor 2 reinforces the ability of this method to provide new molecular insights underlying drug treatment. Thus, this approach may aid in the proposal of novel and personalized treatments.

In order to avoid unwanted effects of current drug interventions, it is necessary to expand the knowledge of the molecular mechanisms related to drug action. Side effects offer insight into drug action, as for example similar side effects of unrelated drugs can be caused by their common off-targets. Moreover, the phenotypes of systematic single gene perturbation screenings in mice strongly contribute to the comprehension of gene function. Here, we present a novel approach that detects molecular interactions of drugs based on the phenotypic similarity of drugs and mouse models. The method is benchmarked with diverse data sets including drug-target interactions as well as gene-drug links of pharmacogenetic associations and validated by in vitro assays.

Design and synthesis of novel protein kinase R (PKR) inhibitors

Protein kinase RNA-activated (PKR) plays an important role in a broad range of intracellular regulatory mechanisms and in the pathophysiology of many human diseases, including microbial and viral infections, cancer, diabetes and neurodegenerative disorders. Recently, several potent PKR inhibitors have been synthesized. However, the enzyme's multifunctional character and a multitude of PKR downstream targets have prevented the successful transformation of such inhibitors into effective drugs. Thus, the need for additional PKR inhibitors remains. With the help of computer-aided drug-discovery tools, we designed and synthesized potential PKR inhibitors. Indeed, two compounds were found to inhibit recombinant PKR in pharmacologically relevant concentrations. One compound, 6-amino-3-methyl-2-oxo-N-phenyl-2,3-dihydro-1H-benzo[d]imidazole-1-carboxamide, also showed anti-apoptotic properties. The novel molecules diversify the existing pool of PKR inhibitors and provide a basis for the future development of compounds based on PKR signal transduction mechanism.

Discovery and cardioprotective effects of the first non-Peptide agonists of the G protein-coupled prokineticin receptor-1

Prokineticins are angiogenic hormones that activate two G protein-coupled receptors: PKR1 and PKR2. PKR1 has emerged as a critical mediator of cardiovascular homeostasis and cardioprotection. Identification of non-peptide PKR1 agonists that contribute to myocardial repair and collateral vessel growth hold promises for treatment of heart diseases. Through a combination of in silico studies, medicinal chemistry, and pharmacological profiling approaches, we designed, synthesized, and characterized the first PKR1 agonists, demonstrating their cardioprotective activity against myocardial infarction (MI) in mice. Based on high throughput docking protocol, 250,000 compounds were computationally screened for putative PKR1 agonistic activity, using a homology model, and 10 virtual hits were pharmacologically evaluated. One hit internalizes PKR1, increases calcium release and activates ERK and Akt kinases. Among the 30 derivatives of the hit compound, the most potent derivative, IS20, was confirmed for its selectivity and specificity through genetic gain- and loss-of-function of PKR1. Importantly, IS20 prevented cardiac lesion formation and improved cardiac function after MI in mice, promoting proliferation of cardiac progenitor cells and neovasculogenesis. The preclinical investigation of the first PKR1 agonists provides a novel approach to promote cardiac neovasculogenesis after MI.

Cross talk between Leishmania donovani CpG DNA and Toll-like receptor 9: an immunoinformatics approach

The precise and potential contribution of Toll-like receptors (TLRs) signaling pathways in fighting parasitic infections of Leishmania spp., an intracellular protozoan parasite, has gained significant attention during the last decades. Although it is well established that TLR9 recognizes CpG motifs in microbial genomes, the specificity of the CpG DNA pattern of Leishmania parasite interacting with endosomal TLR9 is still unknown. Hence in our study to identify the CpG DNA pattern of Leishmania donovani acting as ligand for TLR9, consecutive homology searches were performed using known CpG ODN 2216 as initial template until a consistent CpG pattern in L. donovani was found. A reliable model of TLR9 ectodomains (ECDs) as well as CpG DNA patterns was predicted to develop the 3D structural complexes of TLR9 ECD-CpG DNA utilizing molecular modeling and docking approaches. The results revealed the preferential specificity of L. donovani CpG DNA to TLR9 compared to control ODN and other CpG patterns. The interface between TLR9 and L. donovani CpG DNA was also found to be geometrically complementary with the LRR11 region of TLR9, acting as the critical region for ligand recognition. The L. donovani CpG pattern identified can be employed to derive a platform for development of an innate immunomodulatory agent for deadly disease.

The bitter pill: clinical drugs that activate the human bitter taste receptor TAS2R14

Bitter taste receptors (TAS2Rs) mediate aversive response to toxic food, which is often bitter. These G-protein-coupled receptors are also expressed in extraoral tissues, and emerge as novel targets for therapeutic indications such as asthma and infection. Our goal was to identify ligands of the broadly tuned TAS2R14 among clinical drugs. Molecular properties of known human bitter taste receptor TAS2R14 agonists were incorporated into pharmacophore- and shape-based models and used to computationally predict additional ligands. Predictions were tested by calcium imaging of TAS2R14-transfected HEK293 cells. In vitro testing of the virtual screening predictions resulted in 30-80% success rates, and 15 clinical drugs were found to activate the TAS2R14. hERG potassium channel, which is predominantly expressed in the heart, emerged as a common off-target of bitter drugs. Despite immense chemical diversity of known TAS2R14 ligands, novel ligands and previously unknown polypharmacology of drugs were unraveled by in vitro screening of computational predictions. This enables rational repurposing of traditional and standard drugs for bitter taste signaling modulation for therapeutic indications.

Structure-activity relationships of FMRF-NH2 peptides demonstrate A role for the conserved C terminus and unique N-terminal extension in modulating cardiac contractility

FMRF-NH2 peptides which contain a conserved, identical C-terminal tetrapeptide but unique N terminus modulate cardiac contractility; yet, little is known about the mechanisms involved in signaling. Here, the structure-activity relationships (SARs) of the Drosophila melanogaster FMRF-NH2 peptides, PDNFMRF-NH2, SDNFMRF-NH2, DPKQDFMRF-NH2, SPKQDFMRF-NH2, and TPAEDFMRF-NH2, which bind FMRFa-R, were investigated. The hypothesis tested was the C-terminal tetrapeptide FMRF-NH2, particularly F1, makes extensive, strong ligand-receptor contacts, yet the unique N terminus influences docking and activity. To test this hypothesis, docking, binding, and bioactivity of the C-terminal tetrapeptide and analogs, and the FMRF-NH2 peptides were compared. Results for FMRF-NH2 and analogs were consistent with the hypothesis; F1 made extensive, strong ligand-receptor contacts with FMRFa-R; Y → F (YMRF-NH2) retained binding, yet A → F (AMRF-NH2) did not. These findings reflected amino acid physicochemical properties; the bulky, aromatic residues F and Y formed strong pi-stacking and hydrophobic contacts to anchor the ligand, interactions which could not be maintained in diversity or number by the small, aliphatic A. The FMRF-NH2 peptides modulated heart rate in larva, pupa, and adult distinctly, representative of the contact sites influenced by their unique N-terminal structures. Based on physicochemical properties, the peptides each docked to FMRFa-R with one best pose, except FMRF-NH2 which docked with two equally favorable poses, consistent with the N terminus influencing docking to define specific ligand-receptor contacts. Furthermore, SDNAMRF-NH2 was designed and, despite lacking the aromatic properties of one F, it binds FMRFa-R and demonstrated a unique SAR, consistent with the N terminus influencing docking and conferring binding and activity; thus, supporting our hypothesis.

Targeting angiogenesis in the pathological ovary

The ovary is a key tissue in the study of physiological neo-vascularisation in the adult and its study has highlighted important molecules involved in the regulation of angiogenesis in vivo. These include vascular endothelial growth factor, delta-like ligand 4, thrombospondin-1, prokineticin-1 and prostaglandin E2. Targeting these molecular pathways has therapeutic potential and their manipulation has an increasing preclinical and clinical role in the management of the pathological ovary. Targeting angiogenic pathways has utility in the promotion of ovarian angiogenesis to improve tissue and follicle survival and function as well as the prevention and management of ovarian hyperstimulation syndrome. There is a theoretical possibility that targeting angiogenesis may improve the function of the polycystic ovary and a real role for targeting angiogenesis in ovarian cancer.

Prokineticin receptor variants (PKR1-I379V and PKR2-V331M) are protective genotypes in human early pregnancy

Endocrine gland-derived vascular endothelial growth factor (EG-VEGF) and its receptor genes (PROKR1(PKR1) andPROKR2(PKR2)) play an important role in human early pregnancy. We have previously shown thatPROKR1andPROKR2are associated with recurrent miscarriage (RM) using the tag-SNP method. In this study, we aimed to identifyPROKR1andPROKR2variants in idiopathic RM patients by genotyping of the entire coding regions. Peripheral blood DNA samples of 100 RM women and 100 controls were subjected to sequence the entire exons ofPROKR1andPROKR2. Significant non-synonymous variant genotypes present in the original 200 samples were further confirmed in the extended samples of 144 RM patients and 153 controls. Genetic variants that were over- or under-represented in the patients were ectopically expressed in HEK293 and JAR cells to investigate their effects on intracellular calcium influx, cell proliferation, cell invasion, cell–cell adhesion, and tube organization. We found that the allele and genotype frequencies ofPROKR1(I379V) andPROKR2(V331M) were significantly increased in the normal control groups compared with idiopathic RM women (P<0.05).PROKR1(I379V) andPROKR2(V331M) decreased intracellular calcium influx but increased cell invasiveness (P<0.05), whereas cell proliferation, cell–cell adhesion, and tube organization were not significantly affected. In conclusion,PROKR1(I379V) andPROKR2(V331M) variants conferred lower risk for RM and may play protective roles in early pregnancy by altering calcium signaling and facilitating cell invasiveness.

A combination of receptor-based pharmacophore modeling & QM techniques for identification of human chymase inhibitors

Inhibition of chymase is likely to divulge therapeutic ways for the treatment of cardiovascular diseases, and fibrotic disorders. To find novel and potent chymase inhibitors and to provide a new idea for drug design, we used both ligand-based and structure-based methods to perform the virtual screening(VS) of commercially available databases. Different pharmacophore models generated from various crystal structures of enzyme may depict diverse inhibitor binding modes. Therefore, multiple pharmacophore-based approach is applied in this study. X-ray crystallographic data of chymase in complex with different inhibitors were used to generate four structure-based pharmacophore models. One ligand-based pharmacophore model was also developed from experimentally known inhibitors. After successful validation, all pharmacophore models were employed in database screening to retrieve hits with novel chemical scaffolds. Drug-like hit compounds were subjected to molecular docking using GOLD and AutoDock. Finally four structurally diverse compounds with high GOLD score and binding affinity for several crystal structures of chymase were selected as final hits. Identification of final hits by three different pharmacophore models necessitates the use of multiple pharmacophore-based approach in VS process. Quantum mechanical calculation is also conducted for analysis of electrostatic characteristics of compounds which illustrates their significant role in driving the inhibitor to adopt a suitable bioactive conformation oriented in the active site of enzyme. In general, this study is used as example to illustrate how multiple pharmacophore approach can be useful in identifying structurally diverse hits which may bind to all possible bioactive conformations available in the active site of enzyme. The strategy used in the current study could be appropriate to design drugs for other enzymes as well.

Toll-like receptor 9 interaction with CpG ODN--an in silico analysis approach

Background

Toll-like receptor 9 (TLR9) recognises unmethylated CpG DNA and activates a signalling cascade, leading to the production of inflammatory cytokines such as TNF-α, IL-1, IL-6 and IL-12 via the adaptor protein MyD88. However, the specific sequence and structural requirements of the CpG DNA for the recognition of and binding to TLR9 are unknown. Moreover, the 3D structures of TLR9 and the TLR9-ODN complex have not been determined. In this study, we propose a reliable model of the interaction of the TLR9 ECD with CpG ODN using bioinformatics tools.

Results

The three-dimensional structures of two TLR9 ECD-CpG ODN complexes were constructed using a homology modelling and docking strategy. Based on the models of these complexes, the TLR9 ECD-CpG ODN interaction patterns were calculated. The results showed that the interface between the human TLR9 and the CpG ODN molecule is geometrically complementary. The computed molecular interactions indicated that LRR11 is the main region of TLR9 that binds to CpG ODN and that five positively charged residues within LRR11 are involved in the binding of the TLR9 ECD to the CpG ODN. Observations in the close-up view of these interactions indicated that these five positively charged residues contribute differently to the binding region within the TLR9 ECD-CpG ODN complex. 337Arg and 338Lys reside in the binding sites of ODN, forming hydrogen bonds and direct contacts with the CpG ODN, whereas 347Lys, 348Arg, and 353His do not directly contact the CpG ODN. These results are in agreement with previously reported experimental data.

Conclusion

In this study, we present two structural models for the human and mouse TLR9 ECD in a complex with CpG ODN. Some features predicted by this model are consistent with previously reported experimental data. This complex model may lead to a better understanding of the function of TLR9 and its interaction with CpG ODN and will improve our understanding of TLR9-ligand interaction in general.