Mendelian Disease Associations Reveal Novel Insights into Inflammatory Bowel Disease

Background

Monogenic diseases have been shown to contribute to complex disease risk and may hold new insights into the underlying biological mechanism of Inflammatory Bowel Disease (IBD).

Methods

We analyzed Mendelian disease associations with IBD using over 55 million patients from the Optum's deidentified electronic health records dataset database. Using the significant Mendelian diseases, we performed pathway enrichment analysis and constructed a model using gene expression datasets to differentiate Crohn's disease (CD), ulcerative colitis (UC), and healthy patient samples.

Results

We found 50 Mendelian diseases were significantly associated with IBD, with 40 being significantly associated with both CD and UC. Our results for CD replicated those from previous studies. Pathways that were enriched consisted of mainly immune and metabolic processes with a focus on tolerance and oxidative stress. Our 3-way classifier for UC, CD, and healthy samples yielded an accuracy of 72%.

Conclusions

Mendelian diseases that are significantly associated with IBD may reveal novel insights into the genetic architecture of IBD.

Reverse translation of adverse event reports paves the way for de-risking preclinical off-targets

The Food and Drug Administration Adverse Event Reporting System (FAERS) remains the primary source for post-marketing pharmacovigilance. The system is largely un-curated, unstandardized, and lacks a method for linking drugs to the chemical structures of their active ingredients, increasing noise and artefactual trends. To address these problems, we mapped drugs to their ingredients and used natural language processing to classify and correlate drug events. Our analysis exposed key idiosyncrasies in FAERS, for example reports of thalidomide causing a deadly ADR when used against myeloma, a likely result of the disease itself; multiplications of the same report, unjustifiably increasing its importance; correlation of reported ADRs with public events, regulatory announcements, and with publications. Comparing the pharmacological, pharmacokinetic, and clinical ADR profiles of methylphenidate, aripiprazole, and risperidone, and of kinase drugs targeting the VEGF receptor, demonstrates how underlying molecular mechanisms can emerge from ADR co-analysis. The precautions and methods we describe may enable investigators to avoid confounding chemistry-based associations and reporting biases in FAERS, and illustrate how comparative analysis of ADRs can reveal underlying mechanisms.

DOI:http://dx.doi.org/10.7554/eLife.25818.001

New treatments are tested in clinical trials before they are licensed for use in patients, but until the drugs are available for prescribing it’s not always possible to identify every side effect. When the drugs enter the clinic, they might be prescribed to patients with multiple medical conditions, or combined with other treatments. The drugs may also be taken for longer periods of time than tested in trials. It is therefore common for new adverse reactions to emerge after a drug is in widespread use.

The FDA Adverse Event Reporting System (FAERS) is a surveillance system used in the United States for reporting drug side effects after new treatments have been licensed. Healthcare professionals and patients can submit reports to the database, logging the adverse drug reactions that they have experienced.

FAERS currently contains over 8.5 million entries, and is growing all the time. However, Maciejewski et al. show that the database has several shortcomings that are reducing its usefulness. For instance, on average any given drug will have 16 different names in the system; this makes it challenging to group all of the reported side effects so that trends and patterns can be correctly seen.

To address this first problem, Maciejewski et al. grouped together drugs according to their active ingredients, rather than their name. This made it much easier to account for subsequent, and more crucial conflating factors such as multiple reports for the same adverse event and patient, or cases where adverse reactions were confused with the diseases that the drugs are trying to treat. For example, diabetes was listed as a side effect for drugs used to treat diabetes.

Building on this cleaned-up dataset, Maciejewski et al. monitored how adverse event signals evolve over time and uncovered biases that were hard to see otherwise. For example, side-effects were reported more often when drugs were in the news. More strikingly, this bias affected not only the drug in question, but also other drugs that acted in the same way or on the same molecular target.

The computational method developed by Maciejewski et al. allows the data in FAERS to be combined and corrected, making easier to evaluate the safety of different medicines. The link between adverse side effects and the molecular targets of the drug, via the ingredient’s chemical structure, furthermore makes it possible to analyze such clinical data reliably by using chemical and genetic information. In the future, this method could also help to identify previously unknown side effects and the biological mechanisms behind them. This could help researchers to develop new drugs with improved side effect profiles.

DOI:http://dx.doi.org/10.7554/eLife.25818.002

Implications of Dynamic Occupancy, Binding Kinetics, and Channel Gating Kinetics for hERG Blocker Safety Assessment and Mitigation

Blockade of the hERG potassium channel prolongs the ventricular action potential (AP) and QT interval, and triggers early after depolarizations (EADs) and torsade de pointes (TdP) arrhythmia. Opinions differ as to the causal relationship between hERG blockade and TdP, the relative weighting of other contributing factors, definitive metrics of preclinical proarrhythmicity, and the true safety margin in humans. Here, we have used in silico techniques to characterize the effects of channel gating and binding kinetics on hERG occupancy, and of blockade on the human ventricular AP. Gating effects differ for compounds that are sterically compatible with closed channels (becoming trapped in deactivated channels) versus those that are incompatible with the closed/closing state, and expelled during deactivation. Occupancies of trappable blockers build to equilibrium levels, whereas those of non-trappable blockers build and decay during each AP cycle. Occupancies of ~83% (non-trappable) versus ~63% (trappable) of open/inactive channels caused EADs in our AP simulations. Overall, we conclude that hERG occupancy at therapeutic exposure levels may be tolerated for nontrappable, but not trappable blockers capable of building to the proarrhythmic occupancy level. Furthermore, the widely used Redfern safety index may be biased toward trappable blockers, overestimating the exposure-IC50 separation in nontrappable cases.

[Metameric macular and papular skin mastocytosis]

BACKGROUND

Mastocytosis is characterised by the presence of abnormal quantities of mastocytes in one or more organs. Although it occurs in systemic forms of mastocytosis, isolated skin involvement is the predominant presentation, particularly in children, in the form of more or less extensive though non-systematic lesions. Herein, we report a case of maculopapular cutaneous mastocytosis that is unusual in terms of its metameric topography.

PATIENTS AND METHODS

A 16-year-old youth presented with an erythematous maculopapular rash of 18 months' duration and involving pruritic inflammatory episodes strictly localised in segment T8 to the left. The skin biopsy showed a significant increase in the number of dermal mastocytes (CD117+). No KIT mutations were found in the skin lesions nor in the unimpaired skin of the opposite side. Further investigations ruled out systemic mastocytis.

DISCUSSION

Herein, we report a case of cutaneous mastocytosis that is unusual in terms of its metameric disposition. There have been only two previous reports of segmental cutaneous mastocytis. The two pathological hypotheses involved precessional dermatitis that renders the skin surface susceptible to homing, and somatic mosaicism (type 1) with local mastocyte proliferation.

Public Domain HTS Fingerprints: Design and Evaluation of Compound Bioactivity Profiles from PubChem's Bioassay Repository

Molecular profiling efforts aim at characterizing the biological actions of small molecules by screening them in hundreds of different biochemical and/or cell-based assays. Together, these assays yield a rich data landscape of target-based and phenotypic effects of the tested compounds. However, submitting an entire compound library to a molecular profiling panel can easily become cost-prohibitive. Here, we make use of historical screening assays to create comprehensive bioactivity profiles for more than 300 000 small molecules. These bioactivity profiles, termed PubChem high-throughput screening fingerprints (PubChem HTSFPs), report small molecule activities in 243 different PubChem bioassays. Although the assays originate from originally independently pursued drug or probe discovery projects, we demonstrate their value as molecular signatures when used in combination. We use these PubChem HTSFPs as molecular descriptors in hit expansion experiments for 33 different targets and phenotypes, showing that, on average, they lead to 27 times as many hits in a set of 1000 chosen molecules as a random screening subset of the same size (average ROC score: 0.82). Moreover, we demonstrate that PubChem HTSFPs retrieve hits that are structurally diverse and distinct from active compounds retrieved by chemical similarity-based hit expansion methods. PubChem HTSFPs are made freely available for the chemical biology research community.

Experimental design strategy: weak reinforcement leads to increased hit rates and enhanced chemical diversity

High Throughput Screening (HTS) is a common approach in life sciences to discover chemical matter that modulates a biological target or phenotype. However, low assay throughput, reagents cost, or a flowchart that can deal with only a limited number of hits may impair screening large numbers of compounds. In this case, a subset of compounds is assayed, and in silico models are utilized to aid in iterative screening design, usually to expand around the found hits and enrich subsequent rounds for relevant chemical matter. However, this may lead to an overly narrow focus, and the diversity of compounds sampled in subsequent iterations may suffer. Active learning has been recently successfully applied in drug discovery with the goal of sampling diverse chemical space to improve model performance. Here we introduce a robust and straightforward iterative screening protocol based on naı̈ve Bayes models. Instead of following up on the compounds with the highest scores in the in silico model, we pursue compounds with very low but positive values. This includes unique chemotypes of weakly active compounds that enhance the applicability domain of the model and increase the cumulative hit rates. We show in a retrospective application to 81 Novartis assays that this protocol leads to consistently higher compound and scaffold hit rates compared to a standard expansion around hits or an active learning approach. We recommend using the weak reinforcement strategy introduced herein for iterative screening workflows.

The role of peroxisome-proliferator-activating receptor gamma agonists: rosiglitazone and 15-deoxy-delta12,14-prostaglandin J2 in chronic experimental cyclosporine A-induced nephrotoxicity

Cyclosporine A(CsA) is an immunosuppressor frequently used in the transplant surgery and in the treatment of autoimmune diseases. The therapeutic benefits of CsA are often limited by it's main side effect-nephrotoxicity. Mechanisms of chronic CsA- induced renal damage include: activation of renin-angiotensin-aldosterone system, upregulation of transforming growth factor beta (TGF-β), oxidative stress. This study was undertaken to investigate the protective effect of the peroxisome-proliferator-activated receptors gamma (PPARs-γ) agonists: rosiglitazone and 15-deoxy-Δ12,14-prostaglandin J2 (PGDJ2), against CsA-induced kidney injury in male Wistar rats. CsA was administered subcutaneously at a dose of 15 mg/kg/day for 28 days. Both PPAR-γ agonists were given for 28 days 0.5 hour before the administration of CsA. Rosiglitazone was administered orally at a dose of 8 mg/kg/day and PGDJ2 was given intraperitoneally at a dose of 30 μg/kg/day. CsA induced renal failure was evidenced by increased serum levels of urea, uric acid and creatinine. Serum concentrations of GSH and GSSG, lipid peroxidation products as well as NAD+/NADH, NADP+/NADPH and ADP/ATP ratios showed, that CsA induced oxidative stress and evoked an imbalanced red-ox state in the kidney. Light and electron microscope studies showed degenerative changes within renal tubules with damage to their mitochondria, interstitial fibrosis and arteriolopathy. Immunohistochemical expression of profibrotic TGF-β was assessed. The biochemical and morphological changes induced by CsA were limited by administration of both rosiglitazone and PGDJ2. Ultrastructural examination of renal tubular epithelial cells showed marked improvement within mitochondria. Our results indicate that both PPAR-γ agonists used in the experiment may play an important role in protecting against CsA-induced damage in the kidney.

Translation of off-target effects: prediction of ADRs by integrated experimental and computational approach

Adverse drug reactions (ADRs) are associated with most drugs, often discovered late in drug development and sometimes only during extended course of clinical use.

Decoding the components of dynamics in three-domain proteins

In this study, we examine the feasibility and limitations of describing the motional behavior of three-domain proteins in which the domains are linearly connected. In addition to attempting the determination of the internal and overall reorientational correlation times, we investigate the existence of correlations in the motions between the three domains. Since in linearly arranged three-domain proteins, there are typically no experimental data that can directly report on motional correlation between the first and the third domain, we address this question by dynamics simulations. Two limiting cases occur: (1) for weak repulsive potentials and (2) when strong repulsive potentials are applied between sequential domains. The motions of the terminal domains become correlated in the case of strong interdomain repulsive potentials when these potentials do not allow the angle between the sequential domains to be smaller than about 60°. Using the model-free (MF) and extended MF formalisms of Lipari and Szabo, we find that the motional behavior can be separated into two components; the first component represents the concerted overall motion of the three domains, and the second describes the independent component of the motion of each individual domain. We find that this division of the motional behavior of the protein is maintained only when their timescales are distinct and can be made when the angles between sequential domains remain between 60° and 160°. In this work, we identify and quantify interdomain motional correlations.

Using mutagenesis and structural biology to map the binding site for the Plasmodium falciparum merozoite protein PfRh4 on the human immune adherence receptor

To survive and replicate within the human host, malaria parasites must invade erythrocytes. Invasion can be mediated by the P. falciparum reticulocyte-binding homologue protein 4 (PfRh4) on the merozoite surface interacting with complement receptor type 1 (CR1, CD35) on the erythrocyte membrane. The PfRh4 attachment site lies within the three N-terminal complement control protein modules (CCPs 1-3) of CR1, which intriguingly also accommodate binding and regulatory sites for the key complement activation-specific proteolytic products, C3b and C4b. One of these regulatory activities is decay-accelerating activity. Although PfRh4 does not impact C3b/C4b binding, it does inhibit this convertase disassociating capability. Here, we have employed ELISA, co-immunoprecipitation, and surface plasmon resonance to demonstrate that CCP 1 contains all the critical residues for PfRh4 interaction. We fine mapped by homologous substitution mutagenesis the PfRh4-binding site on CCP 1 and visualized it with a solution structure of CCPs 1-3 derived by NMR and small angle x-ray scattering. We cross-validated these results by creating an artificial PfRh4-binding site through substitution of putative PfRh4-interacting residues from CCP 1 into their homologous positions within CCP 8; strikingly, this engineered binding site had an ∼30-fold higher affinity for PfRh4 than the native one in CCP 1. These experiments define a candidate site on CR1 by which P. falciparum merozoites gain access to human erythrocytes in a non-sialic acid-dependent pathway of merozoite invasion.

Solution structure of CCP modules 10-12 illuminates functional architecture of the complement regulator, factor H

The 155-kDa plasma glycoprotein factor H (FH), which consists of 20 complement control protein (CCP) modules, protects self-tissue but not foreign organisms from damage by the complement cascade. Protection is achieved by selective engagement of FH, via CCPs 1-4, CCPs 6-8 and CCPs 19-20, with polyanion-rich host surfaces that bear covalently attached, activation-specific, fragments of complement component C3. The role of intervening CCPs 9-18 in this process is obscured by lack of structural knowledge. We have concatenated new high-resolution solution structures of overlapping recombinant CCP pairs, 10-11 and 11-12, to form a three-dimensional structure of CCPs 10-12 and validated it by small-angle X-ray scattering of the recombinant triple-module fragment. Superimposing CCP 12 of this 10-12 structure with CCP 12 from the previously solved CCP 12-13 structure yielded an S-shaped structure for CCPs 10-13 in which modules are tilted by 80-110° with respect to immediate neighbors, but the bend between CCPs 10 and 11 is counter to the arc traced by CCPs 11-13. Including this four-CCP structure in interpretation of scattering data for the longer recombinant segments, CCPs 10-15 and 8-15, implied flexible attachment of CCPs 8 and 9 to CCP 10 but compact and intimate arrangements of CCP 14 with CCPs 12, 13 and 15. Taken together with difficulties in recombinant production of module pairs 13-14 and 14-15, the aberrant structure of CCP 13 and the variability of 13-14 linker sequences among orthologues, a structural dependency of CCP 14 on its neighbors is suggested; this has implications for the FH mechanism.

New analogs of the clinical complement inhibitor compstatin with subnanomolar affinity and enhanced pharmacokinetic properties

Therapeutic modulation of the complement system has become increasingly important in line with the growing recognition of the role of complement in numerous diseases. Compstatin, a peptidic inhibitor that acts at the central level of the complement cascade, is currently in clinical evaluation but routes to improve its efficacy have not yet been fully explored. Here, we report improvements in both the inhibitory potency and pharmacokinetic parameters of compstatin that broaden its clinical applications. Selective modification of the compstatin N-terminus with non-proteinogenic amino acids resulted in the first analogue with subnanomolar binding affinity (KD=0.5nM) and other similarly potent derivatives with improved solubility in clinically relevant solvents. Detailed structure-activity relationship studies based on biophysical and computational methods revealed key structural determinants for the observed improvements. Importantly, pharmacokinetic evaluation in non-human primates revealed target-driven elimination kinetics with plasma half-life values exceeding expectations for peptidic drugs (close to 12h). This successful optimization strategy is expected to pave the way for systemic administration of compstatin in a range of clinical conditions.

Estimation of interdomain flexibility of N-terminus of factor H using residual dipolar couplings

Characterization of segmental flexibility is needed to understand the biological mechanisms of the very large category of functionally diverse proteins, exemplified by the regulators of complement activation, that consist of numerous compact modules or domains linked by short, potentially flexible, sequences of amino acid residues. The use of NMR-derived residual dipolar couplings (RDCs), in magnetically aligned media, to evaluate interdomain motion is established but only for two-domain proteins. We focused on the three N-terminal domains (called CCPs or SCRs) of the important complement regulator, human factor H (i.e., FH1-3). These domains cooperate to facilitate cleavage of the key complement activation-specific protein fragment, C3b, forming iC3b that no longer participates in the complement cascade. We refined a three-dimensional solution structure of recombinant FH1-3 based on nuclear Overhauser effects and RDCs. We then employed a rudimentary series of RDC data sets, collected in media containing magnetically aligned bicelles (disklike particles formed from phospholipids) under three different conditions, to estimate interdomain motions. This circumvents a requirement of previous approaches for technically difficult collection of five independent RDC data sets. More than 80% of conformers of this predominantly extended three-domain molecule exhibit flexions of <40°. Such segmental flexibility (together with the local dynamics of the hypervariable loop within domain 3) could facilitate recognition of C3b via initial anchoring and eventual reorganization of modules to the conformation captured in the previously solved crystal structure of a C3b:FH1-4 complex.