James A. R. Dalton

An evaluation of automated homology modelling methods at low target–template sequence similarity

MOTIVATION There are two main areas of difficulty in homology modelling that are particularly important when sequence identity between target and template falls below 50%: sequence alignment and loop building. These problems become magnified with automatic modelling processes, as there is no human input to correct mistakes. As such we have benchmarked several stand-alone strategies that could be implemented in a workflow for automated high-throughput homology modelling. These include three new sequence-structure alignment programs: 3D-Coffee, Staccato and SAlign, plus five homology modelling programs and their respective loop building methods: Builder, Nest, Modeller, SegMod/ENCAD and Swiss-Model. The SABmark database provided 123 targets with at least five templates from the same SCOP family and sequence identities </=50%. RESULTS When using Modeller as the common modelling program, 3D-Coffee outperforms Staccato and SAlign using both multiple templates and the best single template, and across the sequence identity range 20-50%. The mean model RMSD generated from 3D-Coffee using multiple templates is 15 and 28% (or using single templates, 3 and 13%) better than those generated by Staccato and Salign, respectively. 3D-Coffee gives equivalent modelling accuracy from multiple and single templates, but Staccato and SAlign are more successful with single templates, their quality deteriorating as additional lower sequence identity templates are added. Evaluating the different homology modelling programs, on average Modeller performs marginally better in overall modelling than the others tested. However, on average Nest produces the best loops with an 8% improvement by mean RMSD compared to the loops generated by Builder.

An allosteric modulator to control endogenous G protein-coupled receptors with light

Controlling drug activity with light offers the possibility of enhancing pharmacological selectivity with spatial and temporal regulation, thus enabling highly localized therapeutic effects and precise dosing patterns. Here we report on the development and characterization of what is to our knowledge the first photoswitchable allosteric modulator of a G protein-coupled receptor. Alloswitch-1 is selective for the metabotropic glutamate receptor mGlu5 and enables the optical control of endogenous mGlu5 receptors.

TERT Promoter Mutations Are Predictive of Aggressive Clinical Behavior in Patients with Spitzoid Melanocytic Neoplasms

Spitzoid neoplasms constitute a morphologically distinct category of melanocytic tumors, encompassing Spitz nevus (benign), atypical Spitz tumor (intermediate malignant potential), and spitzoid melanoma (fully malignant). Currently, no reliable histopathological criteria or molecular marker is known to distinguish borderline from overtly malignant neoplasms. Because TERT promoter (TERT-p) mutations are common in inherently aggressive cutaneous conventional melanoma, we sought to evaluate their prognostic significance in spitzoid neoplasms. We analyzed tumors labeled as atypical Spitz tumor or spitzoid melanoma from 56 patients with available follow-up data for the association of TERT-p mutations, biallelic CDKN2A deletion, biallelic PTEN deletion, kinase fusions, BRAF/NRAS mutations, nodal status, and histopathological parameters with risk of hematogenous metastasis. Four patients died of disseminated disease and 52 patients were alive and disease free without extranodal metastasis (median follow-up, 32.5 months). We found TERT-p mutations in samples from the 4 patients who developed hematogenous metastasis but in none of tumors from patients who had favorable outcomes. Presence of TERT-p mutations was the most significant predictor of haematogenous dissemination (P < 0.0001) among variables analyzed. We conclude that TERT-p mutations identify a clinically high-risk subset of patients with spitzoid tumors. Application of TERT-p mutational assays for risk stratification in the clinic requires large-scale validation.

Calculation of helix packing angles in protein structures

UNLABELLED Software is presented for the calculation of packing angles and geometry of helical secondary structure elements in protein structures. AVAILABILITY C language source code and documentation is available from http://www.bioinformatics.leeds.ac.uk.

Quantifying conformational changes in GPCRs: glimpse of a common functional mechanism

BackgroundG-protein-coupled receptors (GPCRs) are important drug targets and a better understanding of their molecular mechanisms would be desirable. The crystallization rate of GPCRs has accelerated in recent years as techniques have become more sophisticated, particularly with respect to Class A GPCRs interacting with G-proteins. These developments have made it possible for a quantitative analysis of GPCR geometrical features and binding-site conformations, including a statistical comparison between Class A GPCRs in active (agonist-bound) and inactive (antagonist-bound) states.ResultsHere we implement algorithms for the analysis of interhelical angles, distances, interactions and binding-site volumes in the transmembrane domains of 25 Class A GPCRs (7 active and 18 inactive). Two interhelical angles change in a statistically significant way between average inactive and active states: TM3-TM6 (by -9°) and TM6-TM7 (by +12°). A third interhelical angle: TM5-TM6 shows a trend, changing by -9°. In the transition from inactive to active states, average van der Waals interactions between TM3 and TM7 significantly increase as the average distance between them decreases by >2 Å. Average H-bonding between TM3 and TM6 decreases but is seemingly compensated by an increase in H-bonding between TM5 and TM6. In five Class A GPCRs, crystallized in both active and inactive states, increased H-bonding of agonists to TM6 and TM7, relative to antagonists, is observed. These protein-agonist interactions likely favour a change in the TM6-TM7 angle, which creates a narrowing in the binding pocket of activated receptors and an average ~200 Å3 reduction in volume.ConclusionsIn terms of similar conformational changes and agonist binding pattern, Class A GPCRs appear to share a common mechanism of activation, which can be exploited in future drug development.

The landscape of fusion transcripts in spitzoid melanoma and biologically indeterminate spitzoid tumors by RNA sequencing

Kinase activation by chromosomal translocations is a common mechanism that drives tumorigenesis in spitzoid neoplasms. To explore the landscape of fusion transcripts in these tumors, we performed whole-transcriptome sequencing using formalin-fixed, paraffin-embedded (FFPE) tissues in malignant or biologically indeterminate spitzoid tumors from 7 patients (age 2–14 years). RNA sequence libraries enriched for coding regions were prepared and the sequencing was analyzed by a novel assembly-based algorithm designed for detecting complex fusions. In addition, tumor samples were screened for hotspot TERT promoter mutations, and telomerase expression was assessed by TERT mRNA in situ hybridization (ISH). Two patients had widespread metastasis and subsequently died of disease, and 5 patients had a benign clinical course on limited follow-up (mean: 30 months). RNA sequencing and TERT mRNA ISH were successful in six tumors and unsuccessful in one disseminating tumor because of low RNA quality. RNA sequencing identified a kinase fusion in five of the six sequenced tumors: TPM3–NTRK1 (2 tumors), complex rearrangements involving TPM3, ALK, and IL6R (1 tumor), BAIAP2L1–BRAF (1 tumor), and EML4–BRAF (1 disseminating tumor). All predicted chimeric transcripts were expressed at high levels and contained the intact kinase domain. In addition, two tumors each contained a second fusion gene, ARID1B–SNX9 or PTPRZ1–NFAM1. The detected chimeric genes were validated by home-brew break-apart or fusion fluorescence in situ hybridization (FISH). The two disseminating tumors each harbored the TERT promoter –124C>T (Chr 5:1,295,228 hg19 coordinate) mutation, whereas the remaining five tumors retained the wild-type gene. The presence of the –124C>T mutation correlated with telomerase expression by TERT mRNA ISH. In summary, we demonstrated complex fusion transcripts and novel partner genes for BRAF by RNA sequencing of FFPE samples. The diversity of gene fusions demonstrated by RNA sequencing defines the molecular heterogeneity of spitzoid neoplasms.

Illuminating Phenylazopyridines To Photoswitch Metabotropic Glutamate Receptors: From the Flask to the Animals

Phenylazopyridines are photoisomerizable compounds with high potential to control biological functions with light. We have obtained a series of phenylazopyridines with light dependent activity as negative allosteric modulators (NAM) of metabotropic glutamate receptor subtype 5 (mGlu5). Here we describe the factors needed to achieve an operational molecular photoisomerization and its effective translation into in vitro and in vivo receptor photoswitching, which includes zebrafish larva motility and the regulation of the antinociceptive effects in mice. The combination of light and some specific phenylazopyridine ligands displays atypical pharmacological profiles, including light-dependent receptor overactivation, which can be observed both in vitro and in vivo. Remarkably, the localized administration of light and a photoswitchable compound in the peripheral tissues of rodents or in the brain amygdalae results in an illumination-dependent analgesic effect. The results reveal a robust translation of the phenylazopyridine photoisomerization to a precise photoregulation of biological activity.

Computational Analysis of Negative and Positive Allosteric Modulator Binding and Function in Metabotropic Glutamate Receptor 5 (In)Activation

Metabotropic glutamate receptors (mGluRs) are high-profile G-protein coupled receptors drug targets because of their involvement in several neurological disease states, and mGluR5 in particular is a subtype whose controlled allosteric modulation, both positive and negative, can potentially be useful for the treatment of schizophrenia and relief of chronic pain, respectively. Here we model mGluR5 with a collection of positive and negative allosteric modulators (PAMs and NAMs) in both active and inactive receptor states, in a manner that is consistent with experimental information, using a specialized protocol that includes homology to increase docking accuracy, and receptor relaxation to generate an individual induced fit with each allosteric modulator. Results implicate two residues in particular for NAM and PAM function: NAM interaction with W785 for receptor inactivation, and NAM/PAM H-bonding with S809 for receptor (in)activation. Models suggest the orientation of the H-bond between allosteric modulator and S809, controlled by PAM/NAM chemistry, influences the position of TM7, which in turn influences the shape of the allosteric site, and potentially the receptor state. NAM-bound and PAM-bound mGluR5 models also reveal that although PAMs and NAMs bind in the same pocket and share similar binding modes, they have distinct effects on the conformation of the receptor. Our models, together with the identification of a possible activation mechanism, may be useful in the rational design of new allosteric modulators for mGluR5.

A double effect molecular switch leads to a novel potent negative allosteric modulator of metabotropic glutamate receptor 5

Compounds that modulate the function of G-protein-coupled receptors (GPCRs) by binding to their allosteric sites are of potential interest for the treatment of multiple CNS and non-CNS disorders. Allosteric ligands can act either as positive (PAM), negative (NAM), or silent (SAM) receptor modulators and have numerous advantages over classic orthosteric compounds, including improved GPCR-subtype selectivity; the capacity to adapt to physiological conditions; and better safety profiles. Despite these benefits, allosteric modulators are difficult to design and optimize and are often prone to “molecular switching”: a structural phenomenon by which very subtle chemical variations in the ligand result in unexpected changes in selectivity profiles or pharmacology, changing PAMs to NAMs or vice versa. Here, we report the discovery of a nanomolar and subtype selective NAM of metabotropic glutamate receptor 5 (mGlu5) through a targeted “double effect molecular switch” of a potent mGlu4 PAM, and suggests a promising approach towards the discovery of novel mGluR allosteric modulators.

Shining Light on an mGlu5 Photoswitchable NAM: A Theoretical Perspective

Metabotropic glutamate receptors (mGluRs) are important drug targets because of their involvement in several neurological diseases. Among mGluRs, mGlu5 is a particularly high-profile target because its positive or negative allosteric modulation can potentially treat schizophrenia or anxiety and chronic pain, respectively. Here, we computationally and experimentally probe the functional binding of a novel photoswitchable mGlu5 NAM, termed alloswitch-1, which loses its NAM functionality under violet light. We show alloswitch-1 binds deep in the allosteric pocket in a similar fashion to mavoglurant, the co-crystallized NAM in the mGlu5 transmembrane domain crystal structure. Alloswitch-1, like NAM 2-Methyl-6-(phenylethynyl)pyridine (MPEP), is significantly affected by P655M mutation deep in the allosteric pocket, eradicating its functionality. In MD simulations, we show alloswitch-1 and MPEP stabilize the co-crystallized water molecule located at the bottom of the allosteric site that is seemingly characteristic of the inactive receptor state. Furthermore, both NAMs form H-bonds with S809 on helix 7, which may constitute an important stabilizing interaction for NAM-induced mGlu5 inactivation. Alloswitch-1, through isomerization of its amide group from trans to cis is able to form an additional interaction with N747 on helix 5. This may be an important interaction for amide-containing mGlu5 NAMs, helping to stabilize their binding in a potentially unusual cis-amide state. Simulated conformational switching of alloswitch-1 in silico suggests photoisomerization of its azo group from trans to cis may be possible within the allosteric pocket. However, photoexcited alloswitch-1 binds in an unstable fashion, breaking H-bonds with the protein and destabilizing the co-crystallized water molecule. This suggests photoswitching may have destabilizing effects on mGlu5 binding and functionality.