Mechanistic Insights into the Ligand-Induced Unfolding of an RNA G-Quadruplex.

The cationic porphyrin TMPyP4 is a well-established DNA G-quadruplex (G4) binding ligand that can stabilize different topologies via multiple binding modes. However, TMPyP4 can have both a stabilizing and destabilizing effect on RNA G4 structures. The structural mechanisms that mediate RNA G4 unfolding remain unknown. Here, we report on the TMPyP4-induced RNA G4 unfolding mechanism studied by well-tempered metadynamics (WT-MetaD) with supporting biophysical experiments. The simulations predict a two-state mechanism of TMPyP4 interaction via a groove-bound and a top-face-bound conformation. The dynamics of TMPyP4 stacking on the top tetrad disrupts Hoogsteen H-bonds between guanine bases, resulting in the consecutive TMPyP4 intercalation from top-to-bottom G-tetrads. The results reveal a striking correlation between computational and experimental approaches and validate WT-MetaD simulations as a powerful tool for studying RNA G4-ligand interactions.

Mutations in TFAP2B and previously unimplicated genes of the BMP, Wnt, and Hedgehog pathways in syndromic craniosynostosis.

Craniosynostosis (CS) is a frequent congenital anomaly featuring the premature fusion of 1 or more sutures of the cranial vault. Syndromic cases, featuring additional congenital anomalies, make up 15% of CS. While many genes underlying syndromic CS have been identified, the cause of many syndromic cases remains unknown. We performed exome sequencing of 12 syndromic CS cases and their parents, in whom previous genetic evaluations were unrevealing. Damaging de novo or transmitted loss of function (LOF) mutations were found in 8 genes that are highly intolerant to LOF mutation (P = 4.0 × 10-8); additionally, a rare damaging mutation in SOX11, which has a lower level of intolerance, was identified. Four probands had rare damaging mutations (2 de novo) in TFAP2B, a transcription factor that orchestrates neural crest cell migration and differentiation; this mutation burden is highly significant (P = 8.2 × 10-12). Three probands had rare damaging mutations in GLI2, SOX11, or GPC4, which function in the Hedgehog, BMP, and Wnt signaling pathways; other genes in these pathways have previously been implicated in syndromic CS. Similarly, damaging de novo mutations were identified in genes encoding the chromatin modifier KAT6A, and CTNNA1, encoding catenin α-1. These findings establish TFAP2B as a CS gene, have implications for assessing risk to subsequent children in these families, and provide evidence implicating other genes in syndromic CS. This high yield indicates the value of performing exome sequencing of syndromic CS patients when sequencing of known disease loci is unrevealing.

Dynamic Profiling of ß-Coronavirus 3CL Mpro Protease Ligand-Binding Sites.

ß-coronavirus (CoVs) alone has been responsible for three major global outbreaks in the 21st century. The current crisis has led to an urgent requirement to develop therapeutics. Even though a number of vaccines are available, alternative strategies targeting essential viral components are required as a backup against the emergence of lethal viral variants. One such target is the main protease (Mpro) that plays an indispensable role in viral replication. The availability of over 270 Mpro X-ray structures in complex with inhibitors provides unique insights into ligand-protein interactions. Herein, we provide a comprehensive comparison of all nonredundant ligand-binding sites available for SARS-CoV2, SARS-CoV, and MERS-CoV Mpro. Extensive adaptive sampling has been used to investigate structural conservation of ligand-binding sites using Markov state models (MSMs) and compare conformational dynamics employing convolutional variational auto-encoder-based deep learning. Our results indicate that not all ligand-binding sites are dynamically conserved despite high sequence and structural conservation across ß-CoV homologs. This highlights the complexity in targeting all three Mpro enzymes with a single pan inhibitor.

Selective recognition of c-MYC Pu22 G-quadruplex by a fluorescent probe.

Nucleic acid mimics of fluorescent proteins can be valuable tools to locate and image functional biomolecules in cells. Stacking between the internal G-quartet, formed in the mimics, and the exogenous fluorophore probes constitutes the basis for fluorescence emission. The precision of recognition depends upon probes selectively targeting the specific G-quadruplex in the mimics. However, the design of probes recognizing a G-quadruplex with high selectivity in vitro and in vivo remains a challenge. Through structure-based screening and optimization, we identified a light-up fluorescent probe, 9CI that selectively recognizes c-MYC Pu22 G-quadruplex both in vitro and ex vivo. Upon binding, the biocompatible probe emits both blue and green fluorescence with the excitation at 405 nm. With 9CI and c-MYC Pu22 G-quadruplex complex as the fluorescent response core, a DNA mimic of fluorescent proteins was constructed, which succeeded in locating a functional aptamer on the cellular periphery. The recognition mechanism analysis suggested the high selectivity and strong fluorescence response was attributed to the entire recognition process consisting of the kinetic match, dynamic interaction, and the final stacking. This study implies both the single stacking state and the dynamic recognition process are crucial for designing fluorescent probes or ligands with high selectivity for a specific G-quadruplex structure.

Phosphorylation of Histone Deacetylase 8: Structural and Mechanistic Analysis of the Phosphomimetic S39E Mutant.

Histone deacetylase (HDAC) enzymes that catalyze removal of acetyl-lysine post-translational modifications are frequently post-translationally modified. HDAC8 is phosphorylated within the deacetylase domain at conserved residue serine 39, which leads to decreased catalytic activity. HDAC8 phosphorylation at S39 is unique in its location and function and may represent a novel mode of deacetylation regulation. To better understand the impact of phosphorylation of HDAC8 on enzyme structure and function, we performed crystallographic, kinetic, and molecular dynamics studies of the S39E HDAC8 phosphomimetic mutant. This mutation decreases the level of deacetylation of peptides derived from acetylated nuclear and cytoplasmic proteins. However, the magnitude of the effect depends on the peptide sequence and the identity of the active site metal ion [Zn(II) vs Fe(II)], with the value of kcat/KM for the mutant decreasing 9- to >200-fold compared to that of wild-type HDAC8. Furthermore, the dissociation rate constant of the active site metal ion increases by ∼10-fold. S39E HDAC8 was crystallized in complex with the inhibitor Droxinostat, revealing that phosphorylation of S39, as mimicked by the glutamate side chain, perturbs local structure through distortion of the L1 loop. Molecular dynamics simulations of both S39E and phosphorylated S39 HDAC8 demonstrate that the perturbation of the L1 loop likely occurs because of the lost hydrogen bond between D29 and S39. Furthermore, the S39 perturbation causes structural changes that propagate through the protein scaffolding to influence function in the active site. These data demonstrate that phosphorylation plays an important regulatory role for HDAC8 by affecting ligand binding, catalytic efficiency, and substrate selectivity.

Defect in phosphoinositide signalling through a homozygous variant in PLCB3 causes a new form of spondylometaphyseal dysplasia with corneal dystrophy.

BACKGROUND: Bone dysplasias are a large group of disorders affecting the growth and structure of the skeletal system. METHODS: In the present study, we report the clinical and molecular delineation of a new form of syndromic autosomal recessive spondylometaphyseal dysplasia (SMD) in two Emirati first cousins. They displayed postnatal growth deficiency causing profound limb shortening with proximal and distal segments involvement, narrow chest, radiological abnormalities involving the spine, pelvis and metaphyses, corneal clouding and intellectual disability. Whole genome homozygosity mapping localised the genetic cause to 11q12.1-q13.1, a region spanning 19.32 Mb with ~490 genes. Using whole exome sequencing, we identified four novel homozygous variants within the shared block of homozygosity. Pathogenic variants in genes involved in phospholipid metabolism, such as PLCB4 and PCYT1A, are known to cause bone dysplasia with or without eye anomalies, which led us to select PLCB3 as a strong candidate. This gene encodes phospholipase C ß 3, an enzyme that converts phosphatidylinositol 4,5 bisphosphate (PIP2) to inositol 1,4,5 triphosphate (IP3) and diacylglycerol. RESULTS: The identified variant (c.2632G>T) substitutes a serine for a highly conserved alanine within the Ha2' element of the proximal C-terminal domain. This disrupts binding of the Ha2' element to the catalytic core and destabilises PLCB3. Here we show that this hypomorphic variant leads to elevated levels of PIP2 in patient fibroblasts, causing disorganisation of the F-actin cytoskeleton. CONCLUSIONS: Our results connect a homozygous loss of function variant in PLCB3 with a new SMD associated with corneal dystrophy and developmental delay (SMDCD).

Zoledronic acid and bone cellular respiration.

Phosphorescence O2 analyzer was used to measure calvarial bone cellular respiration (cellular mitochondrial O2 consumption) in Taylor Outbred mice in the presence and absence of zoledronic acid. This potent bisphosphonate inhibits osteoclast-mediated calcium resorption, and its effects on bone respiration have not been previously investigated. The change of O2 concentration with time was measured in closed vials containing phosphate-buffered saline (PBS), 5 mM glucose and 5-25 mg calvarial bone fragments, and it was complex for t = 0-30 h. Cyanide (specific inhibitor of cytochrome oxidase) halted O2 consumption, confirming the oxidation occurred in the respiratory chain. Initial rate of respiration was estimated from the zero-order plots d[O2]/dt for t = 0-4 h. For untreated specimens, the rate (mean ± SD) was 2.0 ± 1.2 µM O2 h-1 mg-1 (n = 6). This value was 7-10 times lower than that of other murine organs, but similar to that reported for rat and Guinea pig calvaria (averaging, 2.7 nmol O2 h-1 mg-1). The corresponding rate in the presence of 10-100 µM zoledronic acid was 2.7 ± 0.7 µM O2 h-1 mg-1 (n = 11), p = 0.216. The first-order plots ln ([O2] t  ÷ [O2] t=0) versus time for t = 0-30 h were also used to compare treated and untreated specimens. The rate (h-1 mg-1 103) for specimens incubated in PBS without glucose was 1.3 ± 0.6 (n = 3, p = 0.007), in PBS + glucose it was 10.7 ± 6.9 (n = 10), in PBS + glucose + 10 µM zoledronic acid it was 12.1 ± 6.7 (n = 10, p = 0.579), in PBS + glucose + 20 µM zoledronic acid it was 12.9 ± 3.3 (n = 9, p = 0.356), and in PBS + glucose + 100 µM zoledronic acid it was 13.7 ± 7.7 (n = 9, p = 0.447). Thus, exposure to high-doses of zoledronic acid over several hours imposed a statistically insignificant increase in calvarial bone cellular respiration.

Folding of guanine quadruplex molecules-funnel-like mechanism or kinetic partitioning? An overview from MD simulation studies.

BACKGROUND: Guanine quadruplexes (GQs) play vital roles in many cellular processes and are of much interest as drug targets. In contrast to the availability of many structural studies, there is still limited knowledge on GQ folding. SCOPE OF REVIEW: We review recent molecular dynamics (MD) simulation studies of the folding of GQs, with an emphasis paid to the human telomeric DNA GQ. We explain the basic principles and limitations of all types of MD methods used to study unfolding and folding in a way accessible to non-specialists. We discuss the potential role of G-hairpin, G-triplex and alternative GQ intermediates in the folding process. We argue that, in general, folding of GQs is fundamentally different from funneled folding of small fast-folding proteins, and can be best described by a kinetic partitioning (KP) mechanism. KP is a competition between at least two (but often many) well-separated and structurally different conformational ensembles. MAJOR CONCLUSIONS: The KP mechanism is the only plausible way to explain experiments reporting long time-scales of GQ folding and the existence of long-lived sub-states. A significant part of the natural partitioning of the free energy landscape of GQs comes from the ability of the GQ-forming sequences to populate a large number of syn-anti patterns in their G-tracts. The extreme complexity of the KP of GQs typically prevents an appropriate description of the folding landscape using just a few order parameters or collective variables. GENERAL SIGNIFICANCE: We reconcile available computational and experimental studies of GQ folding and formulate basic principles characterizing GQ folding landscapes. This article is part of a Special Issue entitled "G-quadruplex" Guest Editor: Dr. Concetta Giancola and Dr. Daniela Montesarchio.

Selectivity of major isoquinoline alkaloids from Chelidonium majus towards telomeric G-quadruplex: A study using a transition-FRET (t-FRET) assay.

BACKGROUND: Natural bioproducts are invaluable resources in drug discovery. Isoquinoline alkaloids of Chelidonium majus constitute a structurally diverse family of natural products that are of great interest, one of them being their selectivity for human telomeric G-quadruplex structure and telomerase inhibition. METHODS: The study focuses on the mechanism of telomerase inhibition by stabilization of telomeric G-quadruplex structures by berberine, chelerythrine, chelidonine, sanguinarine and papaverine. Telomerase activity and mRNA levels of hTERT were estimated using quantitative telomere repeat amplification protocol (q-TRAP) and qPCR, in MCF-7 cells treated with different groups of alkaloids. The selectivity of the main isoquinoline alkaloids of Chelidonium majus towards telomeric G-quadruplex forming sequences were explored using a sensitive modified thermal FRET-melting measurement in the presence of the complementary oligonucleotide CT22. We assessed and monitored G-quadruplex topologies using circular dichroism (CD) methods, and compared spectra to previously well-characterized motifs, either alone or in the presence of the alkaloids. Molecular modeling was performed to rationalize ligand binding to the G-quadruplex structure. RESULTS: The results highlight strong inhibitory effects of chelerythrine, sanguinarine and berberine on telomerase activity, most likely through substrate sequestration. These isoquinoline alkaloids interacted strongly with telomeric sequence G-quadruplex. In comparison, chelidonine and papaverine had no significant interaction with the telomeric quadruplex, while they strongly inhibited telomerase at transcription level of hTERT. Altogether, all of the studied alkaloids showed various levels and mechanisms of telomerase inhibition. CONCLUSIONS: We report on a comparative study of anti-telomerase activity of the isoquinoline alkaloids of Chelidonium majus. Chelerythrine was most effective in inhibiting telomerase activity by substrate sequesteration through G-quadruplex stabilization. GENERAL SIGNIFICANCE: Understanding structural and molecular mechanisms of anti-cancer agents can help in developing new and more potent drugs with fewer side effects. Isoquinolines are the most biologically active agents from Chelidonium majus, which have shown to be telomeric G-quadruplex stabilizers and potent telomerase inhibitors.

Extended molecular dynamics of a c-kit promoter quadruplex.

The 22-mer c-kit promoter sequence folds into a parallel-stranded quadruplex with a unique structure, which has been elucidated by crystallographic and NMR methods and shows a high degree of structural conservation. We have carried out a series of extended (up to 10 µs long, ∼50 µs in total) molecular dynamics simulations to explore conformational stability and loop dynamics of this quadruplex. Unfolding no-salt simulations are consistent with a multi-pathway model of quadruplex folding and identify the single-nucleotide propeller loops as the most fragile part of the quadruplex. Thus, formation of propeller loops represents a peculiar atomistic aspect of quadruplex folding. Unbiased simulations reveal µs-scale transitions in the loops, which emphasizes the need for extended simulations in studies of quadruplex loops. We identify ion binding in the loops which may contribute to quadruplex stability. The long lateral-propeller loop is internally very stable but extensively fluctuates as a rigid entity. It creates a size-adaptable cleft between the loop and the stem, which can facilitate ligand binding. The stability gain by forming the internal network of GA base pairs and stacks of this loop may be dictating which of the many possible quadruplex topologies is observed in the ground state by this promoter quadruplex.

Conformations of Human Telomeric G-Quadruplex Studied Using a Nucleotide-Independent Nitroxide Label.

Guanine-rich oligonucleotides can form a unique G-quadruplex (GQ) structure with stacking units of four guanine bases organized in a plane through Hoogsteen bonding. GQ structures have been detected in vivo and shown to exert their roles in maintaining genome integrity and regulating gene expression. Understanding GQ conformation is important for understanding its inherent biological role and for devising strategies to control and manipulate functions based on targeting GQ. Although a number of biophysical methods have been used to investigate structure and dynamics of GQs, our understanding is far from complete. As such, this work explores the use of the site-directed spin labeling technique, complemented by molecular dynamics simulations, for investigating GQ conformations. A nucleotide-independent nitroxide label (R5), which has been previously applied for probing conformations of noncoding RNA and DNA duplexes, is attached to multiple sites in a 22-nucleotide DNA strand derived from the human telomeric sequence (hTel-22) that is known to form GQ. The R5 labels are shown to minimally impact GQ folding, and inter-R5 distances measured using double electron-electron resonance spectroscopy are shown to adequately distinguish the different topological conformations of hTel-22 and report variations in their occupancies in response to changes of the environment variables such as salt, crowding agent, and small molecule ligand. The work demonstrates that the R5 label is able to probe GQ conformation and establishes the base for using R5 to study more complex sequences, such as those that may potentially form multimeric GQs in long telomeric repeats.

Structure-phenotype correlations of human CYP21A2 mutations in congenital adrenal hyperplasia.

Mutations in the cytochrome p450 (CYP)21A2 gene, which encodes the enzyme steroid 21-hydroxylase, cause the majority of cases in congenital adrenal hyperplasia, an autosomal recessive disorder. To date, more than 100 CYP21A2 mutations have been reported. These mutations can be associated either with severe salt-wasting or simple virilizing phenotypes or with milder nonclassical phenotypes. Not all CYP21A2 mutations have, however, been characterized biochemically, and the clinical consequences of these mutations remain unknown. Using the crystal structure of its bovine homolog as a template, we have constructed a humanized model of CYP21A2 to provide comprehensive structural explanations for the clinical manifestations caused by each of the known disease-causing missense mutations in CYP21A2. Mutations that affect membrane anchoring, disrupt heme and/or substrate binding, or impair stability of CYP21A2 cause complete loss of function and salt-wasting disease. In contrast, mutations altering the transmembrane region or conserved hydrophobic patches cause up to a 98% reduction in enzyme activity and simple virilizing disease. Mild nonclassical disease can result from interference in oxidoreductase interactions, salt-bridge and hydrogen-bonding networks, and nonconserved hydrophobic clusters. A simple in silico evaluation of previously uncharacterized gene mutations could, thus, potentially help predict the often diverse phenotypes of a monogenic disorder.

Dominant Driving Forces in Human Telomere Quadruplex Binding-Induced Structural Alterations.

Recently various pathways of human telomere (ht) DNA folding into G-quadruplexes and of ligand binding to these structures have been proposed. However, the key issue as to the nature of forces driving the folding and recognition processes remains unanswered. In this study, structural changes of 22-mer ht-DNA fragment (Tel22), induced by binding of ions (K(+), Na(+)) and specific bisquinolinium ligands, were monitored by calorimetric and spectroscopic methods and by gel electrophoresis. Using the global model analysis of a wide variety of experimental data, we were able to characterize the thermodynamic forces that govern the formation of stable Tel22 G-quadruplexes, folding intermediates, and ligand-quadruplex complexes, and then predict Tel22 behavior in aqueous solutions as a function of temperature, salt concentration, and ligand concentration. On the basis of the above, we believe that our work sets the framework for better understanding the heterogeneity of ht-DNA folding and binding pathways, and its structural polymorphism.

Conformational dynamics of the human propeller telomeric DNA quadruplex on a microsecond time scale.

The human telomeric DNA sequence with four repeats can fold into a parallel-stranded propeller-type topology. NMR structures solved under molecular crowding experiments correlate with the crystal structures found with crystal-packing interactions that are effectively equivalent to molecular crowding. This topology has been used for rationalization of ligand design and occurs experimentally in a number of complexes with a diversity of ligands, at least in the crystalline state. Although G-quartet stems have been well characterized, the interactions of the TTA loop with the G-quartets are much less defined. To better understand the conformational variability and structural dynamics of the propeller-type topology, we performed molecular dynamics simulations in explicit solvent up to 1.5 µs. The analysis provides a detailed atomistic account of the dynamic nature of the TTA loops highlighting their interactions with the G-quartets including formation of an A:A base pair, triad, pentad and hexad. The results present a threshold in quadruplex simulations, with regards to understanding the flexible nature of the sugar-phosphate backbone in formation of unusual architecture within the topology. Furthermore, this study stresses the importance of simulation time in sampling conformational space for this topology.

Asparagine synthetase deficiency detected by whole exome sequencing causes congenital microcephaly, epileptic encephalopathy and psychomotor delay.

Deficiency of Asparagine Synthetase (ASNSD, MIM 615574) is a very rare autosomal recessive disorder presenting with some brain abnormalities. Affected individuals have congenital microcephaly and progressive encephalopathy associated with severe intellectual disability and intractable seizures. The loss of function of the asparagine synthetase (ASNS, EC 6.3.5.4), particularly in the brain, is the major cause of this particular congenital microcephaly. In this study, we clinically evaluated an affected child from a consanguineous Emirati family presenting with congenital microcephaly and epileptic encephalopathy. In addition, whole-exome sequencing revealed a novel homozygous substitution mutation (c.1193A > C) in the ASNS gene. This mutation resulted in the substitution of highly conserved tyrosine residue by cysteine (p.Y398C). Molecular modeling analysis predicts hypomorphic and damaging effects of this mutation on the protein structure and altering its enzymatic activity. Therefore, we conclude that the loss of ASNS function is most likely the cause of this condition in the studied family. This report brings the number of reported families with this very rare disorder to five and the number of pathogenic mutations in the ASNS gene to four. This finding extends the ASNS pathogenic mutations spectrum and highlights the utility of whole-exome sequencing in elucidation the causes of rare recessive disorders that are heterogeneous and/or overlap with other conditions.

Protectiveness of NAM-based hazard assessment - which testing scope is required?

Hazard assessment (HA) requires toxicity tests to allow deriving protective points of departure (PoDs) for risk assessment irrespective of a compound's mode of action (MoA). The scope of in vitro test batteries (ivTB) thereby necessitated for systemic toxicity is still unclear. We explored the protectiveness regarding systemic toxicity of an ivTB with a scope, which was guided by previous findings from rodent studies, where examining six main targets, including liver and kidney, was sufficient to predict the guideline scope-based PoD with high probability. The ivTB comprises human in vitro models representing liver, kidney, lung and the neuronal system covering transcriptome, mitochondrial dysfunction and neuronal outgrowth. Additionally, 32 CALUX®- and 10 HepG2 BAC-GFP reporters cover a broad range of disturbance mechanisms. Eight compounds were chosen for causing adverse effects such as immunotoxicity or anemia in vivo, i.e., effects not directly covered by assays in the ivTB. PoDs derived from the ivTB and from oral repeated dose studies in rodents were extrapolated to maximum unbound plasma concentrations for comparison. The ivTB-based PoDs were one to five orders of magnitude lower than in vivo PoDs for six of eight compounds, implying that they were protective. The extent of in vitro response varied across test compounds. Especially for hematotoxic substances, the ivTB showed either no response or only cytotoxicity. Assays better capturing this type of hazard would be needed to complement the ivTB. This study highlights the potentially broad applicability of ivTBs for deriving protective PoDs of compounds with unknown MoA.

Animal tests are used to determine which amount of a chemical is toxic ('threshold of toxicity') and which organs are affected. In principle, the threshold can also be derived solely from tests with cultured cells. However, only a limited number of cell types can practically be tested, so one challenge is to determine how many and which types shall be tested. In animal studies, only few organs including liver and kidney are regularly among those most sensitively affected. We explored whether a cell-based test battery representing these sensitive organs and covering important mechanisms of toxicity can be used to derive protective human thresholds. To challenge this approach, eight chemicals were tested that primarily cause effects in organs not directly represented in our test battery. Results provided protective thresholds for most of the investigated compounds and gave indications how to further improve the approach towards a full-fledged replacement for animal tests.

Maternal Undernutrition Induces Cell Signalling and Metabolic Dysfunction in Undifferentiated Mouse Embryonic Stem Cells.

Peri-conceptional environment can induce permanent changes in embryo phenotype which alter development and associate with later disease susceptibility. Thus, mouse maternal low protein diet (LPD) fed exclusively during preimplantation is sufficient to lead to cardiovascular, metabolic and neurological dysfunction in adult offspring. Embryonic stem cell (ESC) lines were generated from LPD and control NPD C57BL/6 blastocysts and characterised by transcriptomics, metabolomics, bioinformatics and molecular/cellular studies to assess early potential mechanisms in dietary environmental programming. Previously, we showed these lines retain cellular and epigenetic characteristics of LPD and NPD embryos after several passages. Here, three main changes were identified in LPD ESC lines. First, their derivation capacity was reduced but pluripotency marker expression was similar to controls. Second, LPD lines had impaired Mitogen-activated protein kinase (MAPK) pathway with altered gene expression of several regulators (e.g., Maff, Rassf1, JunD), reduced ERK1/2 signalling capacity and poorer cell survival characteristics which may contribute to reduced derivation. Third, LPD lines had impaired glucose metabolism comprising reduced upstream enzyme expression (e.g., Gpi, Mpi) and accumulation of metabolites (e.g., glucose-6-P, fructose-6-P) above the phosphofructokinase (PFK) gateway with PFK enzyme activity reduced. ESC lines may therefore permit investigation of peri-conceptional programming mechanisms with reduced need for animal experimentation.

The Identification of Human Translational Biomarkers of Neuropathic Pain and Cross-Species Validation Using an Animal Model.

Neuropathic pain is a common chronic condition, which remains poorly understood. Many patients receiving treatment continue to experience severe pain, due to limited diagnostic/treatment management programmes. The development of objective clinical diagnostic/treatment strategies requires identification of robust biomarkers of neuropathic pain. To this end, we looked to identify biomarkers of chronic neuropathic pain by assessing gene expression profiles in an animal model of neuropathic pain, and differential gene expression in patients to determine the potential translatability. We demonstrated cross-species validation of several genes including those identified through bioinformatic analysis by assessing their expression in blood samples from neuropathic pain patients, according to conservative assessments of significance measured using Bonferroni-corrected p-values. These include CASP5 (p = 0.00226), CASP8 (p = 0.00587), CASP9 (p = 2.09 × 10-9), FPR2 (p = 0.00278), SH3BGRL3 (p = 0.00633), and TMEM88 (p = 0.00038). A ROC analysis revealed several combinations of genes to show high levels of discriminatory power in the comparison of neuropathic pain patients and control participants, of which the combination SH3BGRL3, TMEM88, and CASP9 achieved the highest level (AUROC = 0.923). The CASP9 gene was found to be common in five combinations of three genes revealing the highest levels of discriminatory power. In contrast, the gene combination PLAC8, ROMO1, and A3GALT2 showed the highest levels of discriminatory power in the comparison of neuropathic pain and nociceptive pain (AUROC = 0.919), when patients were grouped by S-LANSS scores. Molecules that demonstrate an active role in neuropathic pain have the potential to be developed into a biological measure for objective diagnostic tests, or as novel drug targets for improved pain management.

SimRFlow: An R-based workflow for automated high-throughput PBPK simulation with the Simcyp® simulator.

SimRFlow is a high-throughput physiologically based pharmacokinetic (PBPK) modelling tool which uses Certara's Simcyp® simulator. The workflow is comprised of three main modules: 1) a Data Collection module for automated curation of physicochemical (from ChEMBL and the Norman Suspect List databases) and experimental data (i.e.: clearance, plasma-protein binding, and blood-to-plasma ratio, from httk-R package databases), 2) a Simulation module which activates the Simcyp® simulator and runs Monte Carlo simulations on virtual subjects using the curated data, and 3) a Data Visualisation module for understanding the simulated compound-specific profiles and predictions. SimRFlow has three administration routes (oral, intravenous, dermal) and allows users to change some simulation parameters including the number of subjects, simulation duration, and dosing. Users are only expected to provide a file of the compounds they wish to simulate, and in return the workflow provides summary statistics, concentration-time profiles of various tissue types, and a database file (containing in-depth results) for each simulated compound. This is presented within a guided and easy-to-use R Shiny interface which provides many plotting options for the visualisation of concentration-time profiles, parameter distributions, trends between the different parameters, as well as comparison of predicted parameters across all batch-simulated compounds. The in-built R functions can be assembled in user-customised scripts which allows for the modification of the workflow for different purposes. SimRFlow proves to be a time-efficient tool for simulating a large number of compounds without any manual curation of physicochemical or experimental data necessary to run Simcyp® simulations.

The Identification of Blood Biomarkers of Chronic Neuropathic Pain by Comparative Transcriptomics.

In this study, we recruited 50 chronic pain (neuropathic and nociceptive) and 43 pain-free controls to identify specific blood biomarkers of chronic neuropathic pain (CNP). Affymetrix microarray was carried out on a subset of samples selected 10 CNP and 10 pain-free control participants. The most significant genes were cross-validated using the entire dataset by quantitative real-time PCR (qRT-PCR). In comparative analysis of controls and CNP patients, WLS (P = 4.80 × 10-7), CHPT1 (P = 7.74 × 10-7) and CASP5 (P = 2.30 × 10-5) were highly significant, whilst FGFBP2 (P = 0.00162), STAT1 (P = 0.00223), FCRL6 (P = 0.00335), MYC (P = 0.00335), XCL2 (P = 0.0144) and GZMA (P = 0.0168) were significant in all CNP patients. A three-arm comparative analysis was also carried out with control as the reference group and CNP samples differentiated into two groups of high and low S-LANSS score using a cut-off of 12. STAT1, XCL2 and GZMA were not significant but KIR3DL2 (P = 0.00838), SH2D1B (P = 0.00295) and CXCR31 (P = 0.0136) were significant in CNP high S-LANSS group (S-LANSS score > 12), along with WLS (P = 8.40 × 10-5), CHPT1 (P = 7.89 × 10-4), CASP5 (P = 0.00393), FGFBP2 (P = 8.70 × 10-4) and FCRL6 (P = 0.00199), suggesting involvement of immune pathways in CNP mechanisms. None of the genes was significant in CNP samples with low (< 12) S-LANSS score. The area under the receiver operating characteristic (AUROC) analysis showed that combination of MYC, STAT1, TLR4, CASP5 and WLS gene expression could be potentially used as a biomarker signature of CNP (AUROC - 0.852, (0.773, 0.931 95% CI)).