@TOME-2: a new pipeline for comparative modeling of protein-ligand complexes

@TOME 3.0: Interfacing Protein Structure Modeling and Ligand Docking

Knowledge of protein-ligand complexes is essential for efficient drug design. Virtual docking can bring important information on putative complexes but it is still far from being simultaneously fast and accurate. Receptors are flexible and adapt to the incoming small molecules while docking is highly sensitive to small conformational deviations. Conformation ensemble is providing a mean to simulate protein flexibility. However, modeling multiple protein structures for many targets is seldom connected to ligand screening in an efficient and straightforward manner. @TOME-3 is an updated version of our former pipeline @TOME-2, in which protein structure modeling is now directly interfaced with flexible ligand docking. Sequence-sequence profile comparisons identify suitable PDB templates for structure modeling and ligands from these templates are used to deduce binding sites to be screened. In addition, bound ligand can be used as pharmacophoric restraint during the virtual docking. The latter is performed by PLANTS while the docking poses are analysed through multiple chemoinformatics functions. This unique combination of tools allows rapid and efficient ligand docking on multiple receptor conformations in parallel. @TOME-3 is freely available on the web at https://atome.cbs.cnrs.fr.

The evolutionary novelty of insect defensins: from bacterial killing to toxin neutralization

Insect host defense comprises two complementary dimensions, microbial killing-mediated resistance and microbial toxin neutralization-mediated resilience, both jointly providing protection against pathogen infections. Insect defensins are a class of effectors of innate immunity primarily responsible for resistance to Gram-positive bacteria. Here, we report a newly originated gene from an ancestral defensin via genetic deletion following gene duplication in Drosophila virilis, which confers an enhanced resilience to Gram-positive bacterial infection. This gene encodes an 18-mer arginine-rich peptide (termed DvirARP) with differences from its parent gene in its pattern of expression, structure and function. DvirARP specifically expresses in D. virilis female adults with a constitutive manner. It adopts a novel fold with a 310 helix and a two CXC motif-containing loop stabilized by two disulfide bridges. DvirARP exhibits no activity on the majority of microorganisms tested and only a weak activity against two Gram-positive bacteria. DvirARP knockout flies are viable and have no obvious defect in reproductivity but they are more susceptible to the DvirARP-resistant Staphylococcus aureus infection than the wild type files, which can be attributable to its ability in neutralization of the S. aureus secreted toxins. Phylogenetic distribution analysis reveals that DvirARP is restrictedly present in the Drosophila subgenus, but independent deletion variations also occur in defensins from the Sophophora subgenus, in support of the evolvability of this class of immune effectors. Our work illustrates for the first time how a duplicate resistance-mediated gene evolves an ability to increase the resilience of a subset of Drosophila species against bacterial infection.

Design, synthesis and preliminary biological evaluation of rivastigmine-INDY hybrids as multitarget ligands against Alzheimer's disease by targeting butyrylcholinesterase and DYRK1A/CLK1 kinases

Based on a multitarget approach implementing rivastigmine-INDY hybrids 1, we identified a set of pseudo-irreversible carbamate-type inhibitors of eqBuChE that, after carbamate transfer at the active site serine residue, released the corresponding INDY analogues 2 endowed with hDYRK1A/hCLK1 kinases inhibitory properties. A SAR study and molecular docking investigation of both series of compounds 1 and 2 revealed that appropriate structural modifications at the carbamate moiety and at the N-appendage of the benzothiazole core led to potent and selective eqBuChE inhibitors with IC50 up to 27 nM and potent hDYRK1A and hCLK1 inhibitors with IC50 up to 106 nM and 17 nM respectively. Pleasingly, identification of the matched pair of compounds 1b/2b with a good balance between inhibition of eqBuChE and hDYRK1A/hCLK1 kinases (IC50 = 68 nM and IC50 = 529/54 nM, respectively) further validated our multitarget approach based on a sequential mechanism of action. In addition, target compound 1b exhibited a suitable ADMET profile, including good brain permeability and high stability in PBS, encouraging further biological investigation as a drug candidate.

ORFans in Mitochondrial Genomes of Marine Polychaete Polydora

Most characterized metazoan mitochondrial genomes are compact and encode a small set of proteins that are essential for oxidative phosphorylation, as well as rRNA and tRNA for their expression. However, in rare cases, invertebrate taxa have additional open reading frames (ORFs) in their mtDNA sequences. Here, we sequenced and analyzed the mitochondrial genome of a polychaete worm, Polydora cf. ciliata, part of whose life cycle takes place in low-oxygen conditions. In the mitogenome, we found three "ORFan" regions (544, 1,060, and 427 bp) that have no resemblance to any standard metazoan mtDNA gene but lack stop codons in one of the reading frames. Similar regions are found in the mitochondrial genomes of three other Polydora species and Bocardiella hamata. All five species share the same gene order in their mitogenomes, which differ from that of other known Spionidae mitogenomes. By analyzing the ORFan sequences, we found that they are under purifying selection pressure and contain conservative regions. The codon adaptation indices (CAIs) of the ORFan genes were in the same range of values as the CAI of conventional protein-coding genes in corresponding mitochondrial genomes. The analysis of the P. cf. ciliata mitochondrial transcriptome showed that ORFan-544, ORFan-427, and a portion of the ORFan-1060 are transcribed. Together, this suggests that ORFan-544 and ORFan-427 encode functional proteins. It is likely that the ORFans originated when the Polydora/Bocardiella species complex separated from the rest of the Spionidae, and this event coincided with massive gene rearrangements in their mitochondrial genomes and tRNA-Met duplication.

Synthesis and structure-activity relationship studies of original cyclic diadenosine derivatives as nanomolar inhibitors of NAD kinase from pathogenic bacteria

Nicotinamide adenine dinucleotide kinases (NAD kinases) are essential and ubiquitous enzymes involved in the production of NADP(H) which is an essential cofactor in many metabolic pathways. Targeting NAD kinase (NADK), a rate limiting enzyme of NADP biosynthesis pathway, represents a new promising approach to treat bacterial infections. Previously, we have produced the first NADK inhibitor active against staphylococcal infection. From this linear di-adenosine derivative, namely NKI1, we designed macrocyclic analogues. Here, we describe the synthesis and evaluation of an original series of cyclic diadenosine derivatives as NADK inhibitors of two pathogenic bacteria, Listeria monocytogenes and Staphylococcus aureus. The nature and length of the link between the two adenosine units were examined leading to sub-micromolar inhibitors of NADK1 from L. monocytogenes, including its most potent in vitro inhibitor reported so far (with a 300-fold improvement compared to NKI1).

Structural and Functional Characterization of CreFH1, the Frataxin Homolog from Chlamydomonas reinhardtii

Frataxin plays a key role in cellular iron homeostasis of different organisms. It has been implicated in iron storage, detoxification, delivery for Fe-S cluster assembly and heme biosynthesis. However, its specific role in iron metabolism remains unclear, especially in photosynthetic organisms. To gain insight into the role and properties of frataxin in algae, we identified the gene CreFH1, which codes for the frataxin homolog from Chlamydomonas reinhardtii. We performed the cloning, expression and biochemical characterization of CreFH1. This protein has a predicted mitochondrial transit peptide and a significant structural similarity to other members of the frataxin family. In addition, CreFH1 was able to form a dimer in vitro, and this effect was increased by the addition of Cu²⁺ and also attenuated the Fenton reaction in the presence of a mixture of Fe²⁺ and H₂O₂. Bacterial cells with overexpression of CreFH1 showed increased growth in the presence of different metals, such as Fe, Cu, Zn and Ni and H₂O₂. Thus, results indicated that CreFH1 is a functional protein that shows some distinctive features compared to its more well-known counterparts, and would play an important role in response to oxidative stress in C. reinhardtii.

UNC45A deficiency causes microvillus inclusion disease-like phenotype by impairing myosin VB-dependent apical trafficking

Variants in the UNC45A cochaperone have been recently associated with a syndrome combining diarrhea, cholestasis, deafness, and bone fragility. Yet the mechanism underlying intestinal failure in UNC45A deficiency remains unclear. Here, biallelic variants in UNC45A were identified by next-generation sequencing in 6 patients with congenital diarrhea. Corroborating in silico prediction, variants either abolished UNC45A expression or altered protein conformation. Myosin VB was identified by mass spectrometry as client of the UNC45A chaperone and was found misfolded in UNC45AKO Caco-2 cells. In keeping with impaired myosin VB function, UNC45AKO Caco-2 cells showed abnormal epithelial morphogenesis that was restored by full-length UNC45A, but not by mutant alleles. Patients and UNC45AKO 3D organoids displayed altered luminal development and microvillus inclusions, while 2D cultures revealed Rab11 and apical transporter mislocalization as well as sparse and disorganized microvilli. All those features resembled the subcellular abnormalities observed in duodenal biopsies from patients with microvillus inclusion disease. Finally, microvillus inclusions and shortened microvilli were evidenced in enterocytes from unc45a-deficient zebrafish. Taken together, our results provide evidence that UNC45A plays an essential role in epithelial morphogenesis through its cochaperone function of myosin VB and that UNC45A loss causes a variant of microvillus inclusion disease.

Protein targeting to starch 1, a functional protein of starch biosynthesis in wheat (Triticum aestivum L.)

TaPTST1, a wheat homolog of AtPTST1 containing CBM can interact with GBSSI and regulate starch metabolism in wheat endosperm. In cereal endosperm, native starch comprising amylose and amylopectin is synthesized by the coordinated activities of several pathway enzymes. Amylose in starch influences its physio-chemical properties resulting in several human health benefits. The Granule-Bound Starch Synthase I (GBSSI) is the most abundant starch-associated protein. GBSSI lacks dedicated Carbohydrate-binding module (CBM). Previously, Protein Targeting To Starch 1 (PTST1) was identified as a crucial protein for the localization of GBSSI to the starch granules in Arabidopsis. The function of its homologous protein in the wheat endosperm is not known. In this study, TaPTST1, an AtPTST1 homolog, containing a CBM and a coiled-coil domain was identified in wheat. Protein-coding nucleotide sequence of TaPTST1 from Indian wheat variety 'C 306' was cloned and characterized. Homology modelling and molecular docking suggested the potential interaction of TaPTST1 with glucans and GBSSI. The TaPTST1 expression was higher in wheat grain than the other tissues, suggesting a grain-specific function. In vitro binding assays demonstrated different binding affinities of TaPTST1 for native starch, amylose, and amylopectin. Furthermore, the immunoaffinity pull-down assay revealed that TaPTST1 directly interacts with GBSSI, and the interaction is mediated by a coiled-coil domain. The direct protein-protein interaction was further confirmed by bimolecular fluorescence complementation assay (BiFC) in planta. Based on our findings we postulate a functional role for TaPTST1 in starch metabolism by targeting GBSSI to starch granules in wheat endosperm.

CBM20CP, a novel functional protein of starch metabolism in green algae

Ostreococcus tauri is a picoalga that contains a small and compact genome, which resembles that of higher plants in the multiplicity of enzymes involved in starch synthesis (ADP-glucose pyrophosphorylase, ADPGlc PPase; granule bound starch synthase, GBSS; starch synthases, SSI, SSII, SSIII; and starch branching enzyme, SBE, between others), except starch synthase IV (SSIV). Although its genome is fully sequenced, there are still many genes and proteins to which no function was assigned. Here, we identify the OT_ostta06g01880 gene that encodes CBM20CP, a plastidial protein which contains a central carbohydrate binding domain of the CBM20 family, and a coiled coil domain at the C-terminus that lacks catalytic activity. We demonstrate that CBM20CP has the ability to bind starch, amylose and amylopectin with different affinities. Furthermore, this protein interacts with OsttaSSIII-B, increasing its binding to starch granules, its catalytic efficiency and promoting granule growth. The results allow us to postulate a functional role for CBM20CP in starch metabolism in green algae. KEY MESSAGE: CBM20CP, a plastidial protein that has a modular structure but lacks catalytic activity, regulates the synthesis of starch in Ostreococcus tauri.

Nitric oxide synthases from photosynthetic organisms improve growth and confer nitrosative stress tolerance in E. coli. Insights on the pterin cofactor

Nitric oxide synthase (NOS) catalyzes NO formation from the substrate l-arginine (Arg). Previously, NOS with distinct biochemical properties were characterized from two photosynthetic microorganisms, the unicellular algae Ostreococcus tauri (OtNOS) and the cyanobacteria Synechococcus PCC 7335 (SyNOS). In this work we studied the effect of recombinant OtNOS and SyNOS expressed under IPTG-induced promoter in E. coli, a bacterium that lacks NOS. Results show that OtNOS and SyNOS expression promote E. coli growth in a nutrient replete medium and allow to better metabolize Arg as N source. In LB medium, OtNOS induces the expression of the NO dioxygenase hmp in E. coli, in accordance with high NO levels visualized with the probe DAF-FM DA. In contrast, SyNOS expression does not induce hmp and show a slight increase of NO production compared to OtNOS. NOS expression reduces ROS production and increases viability of E. coli cultures growing in LB. A strong nitrosative stress provoked by the addition of 1 mM of the NO donors sodium nitroprusside (SNP) and nitrosoglutathione (GSNO) inhibits bacterial growth rate. Under these conditions, the expression of OtNOS or SyNOS counteracts NO donor toxicity restoring bacterial growth. Finally, using bioinformatic tools and ligand docking analyses, we postulate that tetrahydromonapterin (MH4), an endogenous pterin found in E. coli, could act as cofactor required for NOS catalytic activity. Our findings could be useful for the development of biotechnological applications using NOS expression to improve growth in NOS-lacking bacteria.

Quantitative high-throughput phenotypic screening for environmental estrogens using the E-Morph Screening Assay in combination with in silico predictions

BACKGROUND

Exposure to environmental chemicals that interfere with normal estrogen function can lead to adverse health effects, including cancer. High-throughput screening (HTS) approaches facilitate the efficient identification and characterization of such substances.

OBJECTIVES

We recently described the development of the E-Morph Assay, which measures changes at adherens junctions as a clinically-relevant phenotypic readout for estrogen receptor (ER) alpha signaling activity. Here, we describe its further development and application for automated robotic HTS.

METHODS

Using the advanced E-Morph Screening Assay, we screened a substance library comprising 430 toxicologically-relevant industrial chemicals, biocides, and plant protection products to identify novel substances with estrogenic activities. Based on the primary screening data and the publicly available ToxCast dataset, we performed an insilico similarity search to identify further substances with potential estrogenic activity for follow-up hit expansion screening, and built seven insilico ER models using the conformal prediction (CP) framework to evaluate the HTS results.

RESULTS

The primary and hit confirmation screens identified 27 'known' estrogenic substances with potencies correlating very well with the published ToxCast ER Agonist Score (r=+0.95). We additionally detected potential 'novel' estrogenic activities for 10 primary hit substances and for another nine out of 20 structurally similar substances from insilico predictions and follow-up hit expansion screening. The concordance of the E-Morph Screening Assay with the ToxCast ER reference data and the generated CP ER models was 71% and 73%, respectively, with a high predictivity for ER active substances of up to 87%, which is particularly important for regulatory purposes.

DISCUSSION

These data provide a proof-of-concept for the combination of in vitro HTS approaches with insilico methods (similarity search, CP models) for efficient analysis of large substance libraries in order to prioritize substances with potential estrogenic activity for subsequent testing against higher tier human endpoints.

Retinitis Punctata Albescens and RLBP1-Allied Phenotypes

Purpose

To identify relevant criteria for gene therapy based on clinical and genetic characteristics of rod–cone dystrophy associated with RLBP1 pathogenic variants in a large cohort comprising children and adults.

Design

Retrospective cohort study.

Participants

Patients with pathogenic variants in RLBP1 registered in a single French reference center specialized in inherited retinal dystrophies.

Methods

Clinical, multimodal imaging, and genetic findings were reviewed.

Main Outcome Measures

Age of onset; visual acuity; ellipsoid line length; nasal, temporal, and foveal retinal thickness; and pathogenic variants and related phenotypes, including Newfoundland rod–cone and Bothnia dystrophies (NFRCDs), were reappraised.

Results

Twenty-one patients (15 families) were included. The most frequent form was NFRCD with 12 patients (8 families) homozygous for the recurrent deletion of exons 7 through 9 in RLBP1 and 5 patients (4 families) with biallelic protein-truncating variants (2 novel: p.Gln16∗ and p.Tyr251∗). A novel combination of the p.Arg234Trp Bothnia variant with a nonsense variant in trans led to Bothnia dystrophy in 2 sisters. One proband carrying the p.Met266Lys Bothnia variant and in trans p.Arg121Trp and a second, with the p.Arg9Cys and p.Tyr111∗ combination, both demonstrated mild retinitis punctata albescens. Independently of genotype, all patients showed a visual acuity of worse than 20/200, an ellipsoid line width of less than 1000 μm, and a mean foveal thickness of less than 130 to 150 μm, with loss of both the interdigitation and ellipsoid lines.

Conclusions

The eligibility for RLBP1 gene therapy first should be determined according to the biallelic variant combination using a robust classification as proposed herein. An ellipsoid line width of more than 1200 μm and a central thickness of more than 130 to 150 μm with detectable ellipsoid and interdigitation lines should be 2 prerequisite imaging indicators for gene therapy.

Characterization of SdGA, a cold-adapted glucoamylase from Saccharophagus degradans

We investigated the structural and functional properties of SdGA, a glucoamylase (GA) from Saccharophagus degradans, a marine bacterium which degrades different complex polysaccharides at high rate. SdGA is composed mainly by a N-terminal GH15_N domain linked to a C-terminal catalytic domain (CD) found in the GH15 family of glycosylhydrolases with an overall structure similar to other bacterial GAs. The protein was expressed in Escherichia coli cells, purified and its biochemical properties were investigated. Although SdGA has a maximum activity at 39 °C and pH 6.0, it also shows high activity in a wide range, from low to mild temperatures, like cold-adapted enzymes. Furthermore, SdGA has a higher content of flexible residues and a larger CD due to various amino acid insertions compared to other thermostable GAs. We propose that this novel SdGA, is a cold-adapted enzyme that might be suitable for use in different industrial processes that require enzymes which act at low or medium temperatures.

Development of Novel Potential Pleiotropic Compounds of Interest in Alzheimer's Disease Treatment through Rigidification Strategy

The development of Multi-Target Directed Ligand is of clear interest for the treatment of multifactorial pathology such as Alzheimer’s disease (AD). In this context, acetylcholinesterase (AChE) inhibitors have been modulated in order to generate novel pleiotropic compounds targeting a second protein of therapeutic interest in AD. Among them, donecopride was the first example of a dual acetylcholinesterase inhibitor and 5-HT₄ receptor agonist. In order to explore the structural diversity around this preclinical candidate we have explored the preparation of novel constrained analogs through late-stage rigidification strategy. A series of phenylpyrazoles was prepared in a late-stage functionalization process and all compounds were evaluated in vitro towards AChE and 5-HTRs. A docking study was performed in order to better explain the observed SAR towards AChE, 5-HT₄R and 5-HT₆R and this study led to the description of novel ligand targeting both AChE and 5-HT₆R.

New 3-Aryl-2-(2-thienyl)acrylonitriles with High Activity Against Hepatoma Cells

New 2-(thien-2-yl)-acrylonitriles with putative kinase inhibitory activity were prepared and tested for their antineoplastic efficacy in hepatoma models. Four out of the 14 derivatives were shown to inhibit hepatoma cell proliferation at (sub-)micromolar concentrations with IC50 values below that of the clinically relevant multikinase inhibitor sorafenib, which served as a reference. Colony formation assays as well as primary in vivo examinations of hepatoma tumors grown on the chorioallantoic membrane of fertilized chicken eggs (CAM assay) confirmed the excellent antineoplastic efficacy of the new derivatives. Their mode of action included an induction of apoptotic capsase-3 activity, while no contribution of unspecific cytotoxic effects was observed in LDH-release measurements. Kinase profiling of cancer relevant protein kinases identified the two 3-aryl-2-(thien-2-yl)acrylonitrile derivatives 1b and 1c as (multi-)kinase inhibitors with a preferential activity against the VEGFR-2 tyrosine kinase. Additional bioinformatic analysis of the VEGFR-2 binding modes by docking and molecular dynamics calculations supported the experimental findings and indicated that the hydroxy group of 1c might be crucial for its distinct inhibitory potency against VEGFR-2. Forthcoming studies will further unveil the underlying mode of action of the promising new derivatives as well as their suitability as an urgently needed novel approach in HCC treatment.

Structural analysis of alternate sigma factor ComX with RpoC, RpoB and its cognate CIN promoter reveals a distinctive promoter melting mechanism

Alternate sigma factors play a major role in the survival of pathogenic bacteria such as Streptococcus pyogenes in adverse environment conditions. Stress induced sigma factors mediate gene expression under conditions of pathogenesis, dormancy and unusual environmental cues. In the present work, ComX, an alternate sigma factor from S. pyogenes has been characterized. The structures of ComX, RpoB β subunit and RpoC β' subunit of RNA polymerase have been predicted using comparative and homology modelling respectively and validated. Attempts have been made to study RpoB-RpoC-ComX complex interactions with Double Strand (DS) and Single Strand (SS) promoter regions. Stability of these complexes and the promoter melting mechanism have been analysed using Molecular Dynamic (MD) simulations. This study suggests that ComX, although identifies promoter analogous to the alternate sigma factor SigH of M. tuberculosis, follows a distinctive promoter flip out mechanism.Communicated by Ramaswamy H. Sarma.

An oligomeric switch controls the Mrr-induced SOS response in E. coli

Mrr from Escherichia coli K12 is a type IV restriction endonuclease whose role is to recognize and cleave foreign methylated DNA. Beyond this protective role, Mrr can inflict chromosomal DNA damage that elicits the SOS response in the host cell upon heterologous expression of specific methyltransferases such as M.HhaII, or after exposure to high pressure (HP). Activation of Mrr in response to these perturbations involves an oligomeric switch that dissociates inactive homo-tetramers into active dimers. Here we used scanning number and brightness (sN&B) analysis to determine in vivo the stoichiometry of a constitutively active Mrr mutant predicted to be dimeric and examine other GFP-Mrr mutants compromised in their response to either M.HhaII activity or HP shock. We also observed in vitro the direct pressure-induced tetramer dissociation by HP fluorescence correlation spectroscopy of purified GFP-Mrr. To shed light on the linkages between subunit interactions and activity of Mrr and its variants, we built a structural model of the full-length tetramer bound to DNA. Similar to functionally related endonucleases, the conserved DNA cleavage domain would be sequestered by the DNA recognition domain in the Mrr inactive tetramer, dissociating into an enzymatically active dimer upon interaction with multiple DNA sites.

Profilage in silico des inhibiteurs de protéine kinases

Les protéine kinases ont été rapidement identifiées comme favorisant l’apparition de cancers, à travers leur implication dans la régulation du développement et du cycle cellulaire. Il y a une vingtaine d’années, la mise sur le marché des premiers traitements par inhibiteur de protéine kinase, ouvrait la voie vers de nouvelles solutions médicamenteuses plus ciblées contre le cancer. Depuis, nombreuses sont les données structurales et fonctionnelles acquises sur ces cibles thérapeutiques. Les techniques informatiques ont elles aussi évolué, notamment les méthodes d’apprentissage automatique. En tirant parti de la grande quantité d’informations disponibles aujourd’hui, ces méthodes devraient permettre prochainement la prédiction fine de l’interaction d’un inhibiteur donné avec chaque protéine kinase humaine et donc, à terme, la construction d’outils de profilage de leurs inhibiteurs spécifiques. Cette approche intégrative devrait aider la découverte de solutions thérapeutiques anti-cancéreuses plus efficaces et plus sûres.

Iron-Sulfur Cluster Complex Assembly in the Mitochondria of Arabidopsis thaliana

In plants, the cysteine desulfurase (AtNFS1) and frataxin (AtFH) are involved in the formation of Fe-S groups in mitochondria, specifically, in Fe and sulfur loading onto scaffold proteins, and the subsequent formation of the mature Fe-S cluster. We found that the small mitochondrial chaperone, AtISD11, and AtFH are positive regulators for AtNFS1 activity in Arabidopsis. Moreover, when the three proteins were incubated together, a stronger attenuation of the Fenton reaction was observed compared to that observed with AtFH alone. Using pull-down assays, we found that these three proteins physically interact, and sequence alignment and docking studies showed that several amino acid residues reported as critical for the interaction of their human homologous are conserved. Our results suggest that AtFH, AtNFS1 and AtISD11 form a multiprotein complex that could be involved in different stages of the iron-sulfur cluster (ISC) pathway in plant mitochondria.

Hybrid MM/CG Webserver: Automatic Set Up of Molecular Mechanics/Coarse-Grained Simulations for Human G Protein-Coupled Receptor/Ligand Complexes

Hybrid Molecular Mechanics/Coarse-Grained (MM/CG) simulations help predict ligand poses in human G protein-coupled receptors (hGPCRs), the most important protein superfamily for pharmacological applications. This approach allows the description of the ligand, the binding cavity, and the surrounding water molecules at atomistic resolution, while coarse-graining the rest of the receptor. Here, we present the Hybrid MM/CG Webserver (mmcg.grs.kfa-juelich.de) that automatizes and speeds up the MM/CG simulation setup of hGPCR/ligand complexes. Initial structures for such complexes can be easily and efficiently generated with other webservers. The Hybrid MM/CG server also allows for equilibration of the systems, either fully automatically or interactively. The results are visualized online (using both interactive 3D visualizations and analysis plots), helping the user identify possible issues and modify the setup parameters accordingly. Furthermore, the prepared system can be downloaded and the simulation continued locally.

Refined template selection and combination algorithm significantly improves template-based modeling accuracy

In contrast to ab-initio protein modeling methodologies, comparative modeling is considered as the most popular and reliable algorithm to model protein structure. However, the selection of the best set of templates is still a major challenge. An effective template-ranking algorithm is developed to efficiently select only the reliable hits for predicting the protein structures. The algorithm employs the pairwise as well as multiple sequence alignments of template hits to rank and select the best possible set of templates. It captures several key sequences and structural information of template hits and converts into scores to effectively rank them. This selected set of templates is used to model a target. Modeling accuracy of the algorithm is tested and evaluated on TBM-HA domain containing CASP8, CASP9 and CASP10 targets. On an average, this template ranking and selection algorithm improves GDT-TS, GDT-HA and TM_Score by 3.531, 4.814 and 0.022, respectively. Further, it has been shown that the inclusion of structurally similar templates with ample conformational diversity is crucial for the modeling algorithm to maximally as well as reliably span the target sequence and construct its near-native model. The optimal model sampling also holds the key to predict the best possible target structure.

Probing Protein Interaction Networks by Combining MS-Based Proteomics and Structural Data Integration

Protein-protein interactions play a major role in the molecular machinery of life, and various techniques such as AP-MS are dedicated to their identification. However, those techniques return lists of proteins devoid of organizational structure, not detailing which proteins interact with which others. Proposing a hierarchical view of the interactions between the members of the flat list becomes highly tedious for large data sets when done by hand. To help hierarchize this data, we introduce a new bioinformatics protocol that integrates information of the multimeric protein 3D structures available in the Protein Data Bank using remote homology detection, as well as information related to Short Linear Motifs and interaction data from the BioGRID. We illustrate on two unrelated use-cases of different complexity how our approach can be useful to decipher the network of interactions hidden in the list of input proteins, and how it provides added value compared to state-of-the-art resources such as Interactome3D or STRING. Particularly, we show the added value of using homology detection to distinguish between orthologs and paralogs, and to distinguish between core obligate and more facultative interactions. We also demonstrate the potential of considering interactions occurring through Short Linear Motifs.

Virtual screening web servers: designing chemical probes and drug candidates in the cyberspace

The interplay between life sciences and advancing technology drives a continuous cycle of chemical data growth; these data are most often stored in open or partially open databases. In parallel, many different types of algorithms are being developed to manipulate these chemical objects and associated bioactivity data. Virtual screening methods are among the most popular computational approaches in pharmaceutical research. Today, user-friendly web-based tools are available to help scientists perform virtual screening experiments. This article provides an overview of internet resources enabling and supporting chemical biology and early drug discovery with a main emphasis on web servers dedicated to virtual ligand screening and small-molecule docking. This survey first introduces some key concepts and then presents recent and easily accessible virtual screening and related target-fishing tools as well as briefly discusses case studies enabled by some of these web services. Notwithstanding further improvements, already available web-based tools not only contribute to the design of bioactive molecules and assist drug repositioning but also help to generate new ideas and explore different hypotheses in a timely fashion while contributing to teaching in the field of drug development.

Towards accurate high-throughput ligand affinity prediction by exploiting structural ensembles, docking metrics and ligand similarity

MOTIVATION

Nowadays, virtual screening (VS) plays a major role in the process of drug development. Nonetheless, an accurate estimation of binding affinities, which is crucial at all stages, is not trivial and may require target-specific fine-tuning. Furthermore, drug design also requires improved predictions for putative secondary targets among which is Estrogen Receptor alpha (ERα).

RESULTS

VS based on combinations of Structure-Based VS (SBVS) and Ligand-Based VS (LBVS) is gaining momentum to improve VS performances. In this study, we propose an integrated approach using ligand docking on multiple structural ensembles to reflect receptor flexibility. Then, we investigate the impact of the two different types of features (structure-based and ligand molecular descriptors) on affinity predictions using a random forest algorithm. We find that ligand-based features have lower predictive power (rP = 0.69, R2 = 0.47) than structure-based features (rP = 0.78, R2 = 0.60). Their combination maintains high accuracy (rP = 0.73, R2 = 0.50) on the internal test set, but it shows superior robustness on external datasets. Further improvement and extending the training dataset to include xenobiotics, leads to a novel high-throughput affinity prediction method for ERα ligands (rP = 0.85, R2 = 0.71). The presented prediction tool is provided to the community as a dedicated satellite of the @TOME server in which one can upload a ligand dataset in mol2 format and get ligand docked and affinity predicted.

AVAILABILITY AND IMPLEMENTATION

http://edmon.cbs.cnrs.fr.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

ClusPro LigTBM: Automated Template-based Small Molecule Docking

The template-based approach has been essential for achieving high-quality models in the recent rounds of blind protein-protein docking competition CAPRI (Critical Assessment of Predicted Interactions). However, few such automated methods exist for protein-small molecule docking. In this paper, we present an algorithm for template-based docking of small molecules. It searches for known complexes with ligands that have partial coverage of the target ligand, performs conformational sampling and template-guided energy refinement to produce a variety of possible poses, and then scores the refined poses. The algorithm is available as the automated ClusPro LigTBM server. It allows the user to specify the target protein as a PDB file and the ligand as a SMILES string. The server then searches for templates and uses them for docking, presenting the user with top-scoring poses and their confidence scores. The method is tested on the Astex Diverse benchmark, as well as on the targets from the last round of the D3R (Drug Design Data Resource) Grand Challenge. The server is publicly available as part of the ClusPro docking server suite at https://ligtbm.cluspro.org/.

Binding site matching in rational drug design: algorithms and applications

Interactions between proteins and small molecules are critical for biological functions. These interactions often occur in small cavities within protein structures, known as ligand-binding pockets. Understanding the physicochemical qualities of binding pockets is essential to improve not only our basic knowledge of biological systems, but also drug development procedures. In order to quantify similarities among pockets in terms of their geometries and chemical properties, either bound ligands can be compared to one another or binding sites can be matched directly. Both perspectives routinely take advantage of computational methods including various techniques to represent and compare small molecules as well as local protein structures. In this review, we survey 12 tools widely used to match pockets. These methods are divided into five categories based on the algorithm implemented to construct binding-site alignments. In addition to the comprehensive analysis of their algorithms, test sets and the performance of each method are described. We also discuss general pharmacological applications of computational pocket matching in drug repurposing, polypharmacology and side effects. Reflecting on the importance of these techniques in drug discovery, in the end, we elaborate on the development of more accurate meta-predictors, the incorporation of protein flexibility and the integration of powerful artificial intelligence technologies such as deep learning.

In Silico Predictions of Endocrine Disruptors Properties

Endocrine-disrupting chemicals (EDCs) are a broad class of molecules present in our environment that are suspected to cause adverse effects in the endocrine system by interfering with the synthesis, transport, degradation, or action of endogenous ligands. The characterization of the harmful interaction between environmental compounds and their potential cellular targets and the development of robust in vivo, in vitro, and in silico screening methods are important for assessment of the toxic potential of large numbers of chemicals. In this context, computer-aided technologies that will allow for activity prediction of endocrine disruptors and environmental risk assessments are being developed. These technologies must be able to cope with diverse data and connect chemistry at the atomic level with the biological activity at the cellular, organ, and organism levels. Quantitative structure-activity relationship methods became popular for toxicity issues. They correlate the chemical structure of compounds with biological activity through a number of molecular descriptors (e.g., molecular weight and parameters to account for hydrophobicity, topology, or electronic properties). Chemical structure analysis is a first step; however, modeling intermolecular interactions and cellular behavior will also be essential. The increasing number of three-dimensional crystal structures of EDCs' targets has provided a wealth of structural information that can be used to predict their interactions with EDCs using docking and scoring procedures. In the present review, we have described the various computer-assisted approaches that use ligands and targets properties to predict endocrine disruptor activities.

Small-Molecule Inhibitors of Cyclophilins Block Opening of the Mitochondrial Permeability Transition Pore and Protect Mice From Hepatic Ischemia/Reperfusion Injury

BACKGROUND & AIMS

Hepatic ischemia/reperfusion injury is a complication of liver surgery that involves mitochondrial dysfunction resulting from mitochondrial permeability transition pore (mPTP) opening. Cyclophilin D (PPIF or CypD) is a peptidyl-prolyl cis-trans isomerase that regulates mPTP opening in the inner mitochondrial membrane. We investigated whether and how recently created small-molecule inhibitors of CypD prevent opening of the mPTP in hepatocytes and the resulting effects in cell models and livers of mice undergoing ischemia/reperfusion injury.

METHODS

We measured the activity of 9 small-molecule inhibitors of cyclophilins in an assay of CypD activity. The effects of the small-molecule CypD inhibitors or vehicle on mPTP opening were assessed by measuring mitochondrial swelling and calcium retention in isolated liver mitochondria from C57BL/6J (wild-type) and Ppif-/- (CypD knockout) mice and in primary mouse and human hepatocytes by fluorescence microscopy. We induced ischemia/reperfusion injury in livers of mice given a small-molecule CypD inhibitor or vehicle before and during reperfusion and collected samples of blood and liver for histologic analysis.

RESULTS

The compounds inhibited peptidyl-prolyl isomerase activity (half maximal inhibitory concentration values, 0.2-16.2 μmol/L) and, as a result, calcium-induced mitochondrial swelling, by preventing mPTP opening (half maximal inhibitory concentration values, 1.4-132 μmol/L) in a concentration-dependent manner. The most potent inhibitor (C31) bound CypD with high affinity and inhibited swelling in mitochondria from livers of wild-type and Ppif-/- mice (indicating an additional, CypD-independent effect on mPTP opening) and in primary human and mouse hepatocytes. Administration of C31 in mice with ischemia/reperfusion injury before and during reperfusion restored hepatic calcium retention capacity and oxidative phosphorylation parameters and reduced liver damage compared with vehicle.

CONCLUSIONS

Recently created small-molecule inhibitors of CypD reduced calcium-induced swelling in mitochondria from mouse and human liver tissues. Administration of these compounds to mice during ischemia/reperfusion restored hepatic calcium retention capacity and oxidative phosphorylation parameters and reduced liver damage. These compounds might be developed to protect patients from ischemia/reperfusion injury after liver surgery or for other hepatic or nonhepatic disorders related to abnormal mPTP opening.

Binding mode of Pyridoclax to myeloid cell leukemia-1 (Mcl-1) revealed by nuclear magnetic resonance spectroscopy, docking and molecular dynamics approaches

Myeloid cell leukemia-1 (Mcl-1) is an anti-apoptotic member of the Bcl-2 family proteins. Its amplification is one of the most frequent genetic aberrations found in human cancers. Pyridoclax, a promising BH3 mimetic inhibitor, interacts directly with Mcl-1 and induces massive apoptosis at a concentration of 15 µM in combination with anti-Bcl-x_L strategies in chemo-resistant ovarian cancer cell lines. In this study, a combined experimental and theoretical approach was used to investigate the binding mode of Pyridoclax to Mcl-1. The representative poses generated from dynamics simulations compared with NMR data revealed: (i) Pyridoclax bound to P1 and P2 pockets of Mcl-1 BH3 binding groove through its styryl and methyl groups establishing mainly hydrophobic contacts, (ii) one of the ending pyridines interacts through electrostatic interaction with K234 side chain, a negatively charged residue present only in this position in Mcl-1. Communicated by Ramaswamy H. Sarma.

The ellagitannin metabolite urolithin C is a glucose-dependent regulator of insulin secretion through activation of L-type calcium channels

BACKGROUND AND PURPOSE

The pharmacology of polyphenol metabolites on beta-cell function is largely undetermined. We sought to identify polyphenol metabolites that enhance the insulin-secreting function of beta-cells and to explore the underlying mechanisms.

EXPERIMENTAL APPROACH

INS-1 beta-cells and rat isolated islets of Langerhans or perfused pancreas preparations were used for insulin secretion experiments. Molecular modelling, intracellular Ca2+ monitoring, and whole-cell patch-clamp recordings were used for mechanistic studies.

KEY RESULTS

Among a set of polyphenol metabolites, we found that exposure of INS-1 beta-cells to urolithins A and C enhanced glucose-stimulated insulin secretion. We further characterized the activity of urolithin C and its pharmacological mechanism. Urolithin C glucose-dependently enhanced insulin secretion in isolated islets of Langerhans and perfused pancreas preparations. In the latter, enhancement was reversible when glucose was lowered from a stimulating to a non-stimulating concentration. Molecular modelling suggested that urolithin C could dock into the Cav 1.2 L-type Ca2+ channel. Calcium monitoring indicated that urolithin C had no effect on basal intracellular Ca2+ but enhanced depolarization-induced increase in intracellular Ca2+ in INS-1 cells and dispersed cells isolated from islets. Electrophysiology studies indicated that urolithin C dose-dependently enhanced the L-type Ca2+ current for levels of depolarization above threshold and shifted its voltage-dependent activation towards more negative potentials in INS-1 cells.

CONCLUSION AND IMPLICATIONS

Urolithin C is a glucose-dependent activator of insulin secretion acting by facilitating L-type Ca2+ channel opening and Ca2+ influx into pancreatic beta-cells. Our work paves the way for the design of polyphenol metabolite-inspired compounds aimed at ameliorating beta-cell function.

In silico chemical library screening and experimental validation of novel compounds with potential varroacide activities

The mite Varroa destructor is an ectoparasite and has been identified as a major cause of worldwide honey bee colony losses. The use of yearly treatments for the control of varroosis is the most common answer to prevent collapses of honey bee colonies due to the mite. However, the number of effective acaricides is small and the mite tends to become resistant to these few active molecules. In this study, we have been looking for a new original varroacide treatment inhibiting selectively Varroa destructor AChE (vdAChE) with respect to Apis mellifera AChE (amAChE). To do this an original drug design methodology was used applying virtual screening of the CERMN chemolibrary, starting from a vdAChE homology sequence model. By combining the in silico screening with in vitro experiments, two promising compounds were found. In vitro tests of AChE inhibition for both species have confirmed good selectivity toward the mite vdAChE. Moreover, an in vivo protocol was performed and highlighted a varroacide activity without acute consequences on honey bee survival. The two compounds discovered have the potential to become new drug leads for the development of new treatments against the mite varroa. The method described here clearly shows the potential of a drug-design approach to develop new solutions to safeguard honey bee health.

Comprehensive classification of the plant non-specific lipid transfer protein superfamily towards its sequence-structure-function analysis

Background

Non-specific Lipid Transfer Proteins (nsLTPs) are widely distributed in the plant kingdom and constitute a superfamily of related proteins. Several hundreds of different nsLTP sequences—and counting—have been characterized so far, but their biological functions remain unclear. It has been clear for years that they present a certain interest for agronomic and nutritional issues. Deciphering their functions means collecting and analyzing a variety of data from gene sequence to protein structure, from cellular localization to the physiological role. As a huge and growing number of new protein sequences are available nowadays, extracting meaningful knowledge from sequence–structure–function relationships calls for the development of new tools and approaches. As nsLTPs show high evolutionary divergence, but a conserved common right handed superhelix structural fold, and as they are involved in a large number of key roles in plant development and defense, they are a stimulating case study for validating such an approach.

Methods

In this study, we comprehensively investigated 797 nsLTP protein sequences, including a phylogenetic analysis on canonical protein sequences, three-dimensional structure modeling and functional annotation using several well-established bioinformatics programs. Additionally, two integrative methodologies using original tools were developed. The first was a new method for the detection of (i) conserved amino acid residues involved in structure stabilization and (ii) residues potentially involved in ligand interaction. The second was a structure–function classification based on the evolutionary trace display method using a new tree visualization interface. We also present a new tool for visualizing phylogenetic trees.

Results

Following this new protocol, an updated classification of the nsLTP superfamily was established and a new functional hypothesis for key residues is suggested. Lastly, this work allows a better representation of the diversity of plant nsLTPs in terms of sequence, structure and function.

Novel and Selective TLR7 Antagonists among the Imidazo[1,2-a]pyrazines, Imidazo[1,5-a]quinoxalines, and Pyrazolo[1,5-a]quinoxalines Series

The Toll-like receptors (TLRs) 7 and 8 play an important role in the immune system activation, and their agonists may therefore serve as promising candidate vaccine adjuvants. However, the chronic immune activation by excessive TLR stimulation is a hallmark of several clinically important infectious and autoimmune diseases, which warrants the search for TLR antagonists. In this study, we have synthesized and characterized a variety of compounds belonging to three heterocyclic chemical series: imidazo[1,2-a]pyrazine, imidazo[1,5-a]quinoxaline, and pyrazolo[1,5-a]quinoxaline. These compounds have been tested for their TLR7 or TLR8 agonistic and antagonistic activities. Several of them are shown to be selective TLR7 antagonists without any TLR7 or TLR8 agonistic activity. The selectivity was confirmed by a comparative ligand-docking study in TLR7 antagonist pocket. Two compounds of the pyrazolo[1,5-a]quinoxaline series (10a and 10b) are potent selective TLR7 antagonists and may be considered as promising starting points for the development of new therapeutic agents.

A Tiny Change Makes a Big Difference in the Anti-Parasitic Activities of an HDAC Inhibitor

We previously synthesized an hydroxamate derivative (N-hydroxy-4-[2-(3-methoxyphenyl)acetamido]benzamide) named 363 with potent anti-Toxoplasma gondii activity and histone deacetylase inhibitor (HDACi) effects. Here we show that 1-N-hydroxy-4-N-[(2-methoxyphenyl)methyl]benzene-1,4-dicarboxamide, a 363 isomer, does not have antiparasitic potency and has a 13-fold decrease in HDACi activity. The in silico modeling of T. gondii HDACs of the type II strain discloses identity varying from 25% to 62% on more than 250 residues for S8EP32_TOXG and A0A125YPH4_TOXGM. We observed a high conservation degree with the human HDAC2 (53% and 64% identity, respectively) and a moderate one with the human HDAC8 (30–40%). Two other TgHDACs, S8F6L4_TOXGM and S8GEI3_TOXGM, were identified as displaying a higher similarity with some bacterial orthologs (~35%) than with the human enzymes (~25%). The docking in parallel of the two compounds on the models generated allowed us to gain insights on the docking of these hydroxamate derivatives that guide their specificity and potency against T. gondii histone deacetylase. This information would constitute the rationale from which more specific derivatives can be synthetized.

Insights into Substrate and Inhibitor Selectivity among Human GLUT Transporters through Comparative Modeling and Molecular Docking

The solute carrier 2 family is composed of 14 transporters, which are members of the major facilitator superfamily. Despite their high physiological importance, there are still many open questions concerning their function and specificity, and in some cases, their physiological substrate is still unknown. To understand the determinants of the substrate and inhibitor specificity, we modeled all human glucose transport carriers (GLUTs) and simulated their interaction with known ligands. Comparative modeling was performed with the @TOME-2 pipeline, employing multiple templates and providing an ensemble of models for each GLUT. We analyzed models in both outward-occluded and inward-open conformations, to compare exofacial and endofacial binding sites throughout the family and understand differences in susceptibility of GLUTs to the inhibitor cytochalasin B. Finally, we employed molecular docking and bioinformatics to identify residues likely critical for recognition of myo-inositol by GLUT13 and urate by GLUT9. These results provide insights into the molecular basis for the specificity for these substrates. In addition, we suggested a potential recognition site of glucosamine by GLUT11 to be evaluated in future experiments.

First-in-class allosteric inhibitors of bacterial IMPDHs

Inosine-5'-monophosphate dehydrogenase (IMPDH) is an essential enzyme in many bacterial pathogens and is considered as a potential drug target for the development of new antibacterial agents. Our recent work has revealed the crucial role of one of the two structural domains (i.e. Bateman domain) in the regulation of the quaternary structure and enzymatic activity of bacterial IMPDHs. Thus, we have screened chemical libraries to search for compounds targeting the Bateman domain and identified first in-class allosteric inhibitors of a bacterial IMPDH. These inhibitors were shown to counteract the activation by the natural positive effector, MgATP, and to block the enzyme in its apo conformation (low affinity for IMP). Our structural studies demonstrate the versatility of the Bateman domain to accommodate totally unrelated chemical scaffolds and pave the way for the development of allosteric inhibitors, an avenue little explored until now.

Structural basis for two efficient modes of agropinic acid opine import into the bacterial pathogen Agrobacterium tumefaciens

Agrobacterium tumefaciens pathogens genetically modify their host plants to drive the synthesis of opines in plant tumors. The mannityl-opine family encompasses mannopine, mannopinic acid, agropine and agropinic acid. These opines serve as nutrients and are imported into bacteria via periplasmic-binding proteins (PBPs) in association with ABC transporters. Structural and affinity data on agropine and agropinic acid opines bound to PBPs are currently lacking. Here, we investigated the molecular basis of AgtB and AgaA, proposed as the specific PBP for agropine and agropinic acid import, respectively. Using genetic approaches and affinity measurements, we identified AgtB and its transporter as responsible for agropine uptake in agropine-assimilating agrobacteria. Nonetheless, we showed that AgtB binds agropinic acid with a higher affinity than agropine, and we structurally characterized the agropinic acid-binding mode through three crystal structures at 1.4, 1.74 and 1.9 Å resolution. In the crystallization time course, obtaining a crystal structure of AgtB with agropine was unsuccessful due to the spontaneous lactamization of agropine into agropinic acid. AgaA binds agropinic acid only with a similar affinity in nanomolar range as AgtB. The structure of AgaA bound to agropinic acid at 1.65 Å resolution defines a different agropinic acid-binding signature. Our work highlights the structural and functional characteristics of two efficient agropinic acid assimilation pathways, of which one is also involved in agropine assimilation.

In Silico Target Druggability Assessment: From Structural to Systemic Approaches

This chapter will focus on today's in silico direct and indirect approaches to assess therapeutic target druggability. The direct approach tries to infer from the 3D structure the capacity of the target protein to bind small molecule in order to modulate its biological function. Algorithms to recognize and characterize the quality of the ligand interaction sites whether within buried protein cavities or within large protein-protein interface will be reviewed in the first part of the paper. In the case a ligand-binding site is already identified, indirect aspects of target druggability can be assessed. These indirect approaches focus first on target promiscuity and the potential difficulties in developing specific drugs. It is based on large-scale comparison of protein-binding sites. The second aspect concerns the capacity of the target to induce resistant pathway once it is inhibited or activated by a drug. The emergence of drug-resistant pathways can be assessed through systemic analysis of biological networks implementing metabolism and/or cell regulation signaling.

Characterization of the Anti-Hepatitis C Virus Activity of New Nonpeptidic Small-Molecule Cyclophilin Inhibitors with the Potential for Broad Anti-Flaviviridae Activity

Although members of the Flaviviridae display high incidence, morbidity, and mortality rates, the development of specific antiviral drugs for each virus is unlikely. Cyclophilins, a family of host peptidyl-prolyl cis-trans isomerases (PPIases), play a pivotal role in the life cycles of many viruses and therefore represent an attractive target for broad-spectrum antiviral development. We report here the pangenotypic anti-hepatitis C virus (HCV) activity of a small-molecule cyclophilin inhibitor (SMCypI). Mechanistic and modeling studies revealed that the SMCypI bound to cyclophilin A in competition with cyclosporine (CsA), inhibited its PPIase activity, and disrupted the CypA-nonstructural protein 5A (NS5A) interaction. Resistance selection showed that the lead SMCypI hardly selected amino acid substitutions conferring low-level or no resistance in vitro. Interestingly, the SMCypI selected D320E and Y321H substitutions, located in domain II of the NS5A protein. These substitutions were previously associated with low-level resistance to cyclophilin inhibitors such as alisporivir. Finally, the SMCypI inhibited the replication of other members of the Flaviviridae family with higher 50% effective concentrations (EC₅₀s) than for HCV. Thus, because of its chemical plasticity and simplicity of synthesis, our new family of SMCypIs represents a promising new class of drugs with the potential for broad-spectrum anti-Flaviviridae activity as well as an invaluable tool to explore the role of cyclophilins in viral life cycles.

MinOmics, an Integrative and Immersive Tool for Multi-Omics Analysis

Proteomic and transcriptomic technologies resulted in massive biological datasets, their interpretation requiring sophisticated computational strategies. Efficient and intuitive real-time analysis remains challenging. We use proteomic data on 1417 proteins of the green microalga Chlamydomonas reinhardtii to investigate physicochemical parameters governing selectivity of three cysteine-based redox post translational modifications (PTM): glutathionylation (SSG), nitrosylation (SNO) and disulphide bonds (SS) reduced by thioredoxins. We aim to understand underlying molecular mechanisms and structural determinants through integration of redox proteome data from gene- to structural level. Our interactive visual analytics approach on an 8.3 m² display wall of 25 MPixel resolution features stereoscopic three dimensions (3D) representation performed by UnityMol WebGL. Virtual reality headsets complement the range of usage configurations for fully immersive tasks. Our experiments confirm that fast access to a rich cross-linked database is necessary for immersive analysis of structural data. We emphasize the possibility to display complex data structures and relationships in 3D, intrinsic to molecular structure visualization, but less common for omics-network analysis. Our setup is powered by MinOmics, an integrated analysis pipeline and visualization framework dedicated to multi-omics analysis. MinOmics integrates data from various sources into a materialized physical repository. We evaluate its performance, a design criterion for the framework.

HadD, a novel fatty acid synthase type II protein, is essential for alpha- and epoxy-mycolic acid biosynthesis and mycobacterial fitness

Mycolic acids (MAs) have a strategic location within the mycobacterial envelope, deeply influencing its architecture and permeability, and play a determinant role in the pathogenicity of mycobacteria. The fatty acid synthase type II (FAS-II) multienzyme system is involved in their biosynthesis. A combination of pull-downs and proteomics analyses led to the discovery of a mycobacterial protein, HadD, displaying highly specific interactions with the dehydratase HadAB of FAS-II. In vitro activity assays and homology modeling showed that HadD is, like HadAB, a hot dog folded (R)-specific hydratase/dehydratase. A hadD knockout mutant of Mycobacterium smegmatis produced only the medium-size alpha’-MAs. Data strongly suggest that HadD is involved in building the third meromycolic segment during the late FAS-II elongation cycles, leading to the synthesis of the full-size alpha- and epoxy-MAs. The change in the envelope composition induced by hadD inactivation strongly altered the bacterial fitness and capacities to aggregate, assemble into colonies or biofilms and spread by sliding motility, and conferred a hypersensitivity to the firstline antimycobacterial drug rifampicin. This showed that the cell surface properties and the envelope integrity were greatly affected. With the alarmingly increasing case number of nontuberculous mycobacterial diseases, HadD appears as an attractive target for drug development.

Synthesis of aminophenylhydroxamate and aminobenzylhydroxamate derivatives and in vitro screening for antiparasitic and histone deacetylase inhibitory activity

A series of aminophenylhydroxamates and aminobenzylhydroxamates were synthesized and screened for their antiparasitic activity against Leishmania, Trypanosoma, and Toxoplasma. Their anti-histone deacetylase (HDAC) potency was determined. Moderate to no antileishmanial or antitrypanosomal activity was found (IC50 > 10 μM) that contrast with the highly efficient anti-Toxoplasma activity (IC50 < 1.0 μM) of these compounds. The antiparasitic activity of the synthetized compounds correlates well with their HDAC inhibitory activity. The best-performing compound (named 363) express a high anti-HDAC6 inhibitory activity (IC50 of 0.045 ± 0.015 μM) a moderate cytotoxicity and a high anti-Toxoplasma activity in the range of known anti-Toxoplasma compounds (IC50 of 0.35-2.25 μM). The calculated selectivity index (10-300 using different human cell lines) of the compound 363 makes it a lead compound for the future development of anti-Toxoplasma molecules.

Dimerization of the Pragmin Pseudo-Kinase Regulates Protein Tyrosine Phosphorylation

The pseudo-kinase and signaling protein Pragmin has been linked to cancer by regulating protein tyrosine phosphorylation via unknown mechanisms. Here we present the crystal structure of the Pragmin 906-1,368 amino acid C terminus, which encompasses its kinase domain. We show that Pragmin contains a classical protein-kinase fold devoid of catalytic activity, despite a conserved catalytic lysine (K997). By proteomics, we discovered that this pseudo-kinase uses the tyrosine kinase CSK to induce protein tyrosine phosphorylation in human cells. Interestingly, the protein-kinase domain is flanked by N- and C-terminal extensions forming an original dimerization domain that regulates Pragmin self-association and stimulates CSK activity. A1329E mutation in the C-terminal extension destabilizes Pragmin dimerization and reduces CSK activation. These results reveal a dimerization mechanism by which a pseudo-kinase can induce protein tyrosine phosphorylation. Further sequence-structure analysis identified an additional member (C19orf35) of the superfamily of dimeric Pragmin/SgK269/PEAK1 pseudo-kinases.

Culex pipiens crossing type diversity is governed by an amplified and polymorphic operon of Wolbachia

Culex pipiens mosquitoes are infected with Wolbachia (wPip) that cause an important diversity of cytoplasmic incompatibilities (CIs). Functional transgenic studies have implicated the cidA-cidB operon from wPip and its homolog in wMel in CI between infected Drosophila males and uninfected females. However, the genetic basis of the CI diversity induced by different Wolbachia strains was unknown. We show here that the remarkable diversity of CI in the C. pipiens complex is due to the presence, in all tested wPip genomes, of several copies of the cidA-cidB operon, which undergoes diversification through recombination events. In 183 isofemale lines of C. pipiens collected worldwide, specific variations of the cidA-cidB gene repertoires are found to match crossing types. The diversification of cidA-cidB is consistent with the hypothesis of a toxin–antitoxin system in which the gene cidB co-diversifies with the gene cidA, particularly in putative domains of reciprocal interactions.

Wolbachia causes cytoplasmic incompatibility (CI) between mosquitoes infected with different strains, but the genetic basis of observed CI diversity is unknown. Here, Bonneau et al. sequence Wolbachia from over 100 Culex pipiens lines and show that crossing types match variations of the toxin-antitoxin cidA-cidB genes.

Understanding the Role of Intrinsic Disorder of Viral Proteins in the Oncogenicity of Different Types of HPV

Intrinsic disorder is very important in the biological function of several proteins, and is directly linked to their foldability during interaction with their targets. There is a close relationship between the intrinsically disordered proteins and the process of carcinogenesis involving viral pathogens. Among these pathogens, we have highlighted the human papillomavirus (HPV) in this study. HPV is currently among the most common sexually transmitted infections, besides being the cause of several types of cancer. HPVs are divided into two groups, called high- and low-risk, based on their oncogenic potential. The high-risk HPV E6 protein has been the target of much research, in seeking treatments against HPV, due to its direct involvement in the process of cell cycle control. To understand the role of intrinsic disorder of the viral proteins in the oncogenic potential of different HPV types, the structural characteristics of intrinsically disordered regions of high and low-risk HPV E6 proteins were analyzed. In silico analyses of primary sequences, prediction of tertiary structures, and analyses of molecular dynamics allowed the observation of the behavior of such disordered regions in these proteins, thereby proving a direct relationship of structural variation with the degree of oncogenicity of HPVs. The results obtained may contribute to the development of new therapies, targeting the E6 oncoprotein, for the treatment of HPV-associated diseases.

Studying Kv Channels Function using Computational Methods

In recent years, molecular modeling techniques, combined with MD simulations, provided significant insights on voltage-gated (Kv) potassium channels intrinsic properties. Among the success stories are the highlight of molecular level details of the effects of mutations, the unraveling of several metastable intermediate states, and the influence of a particular lipid, PIP₂, in the stability and the modulation of Kv channel function. These computational studies offered a detailed view that could not have been reached through experimental studies alone. With the increase of cross disciplinary studies, numerous experiments provided validation of these computational results, which endows an increase in the reliability of molecular modeling for the study of Kv channels. This chapter offers a description of the main techniques used to model Kv channels at the atomistic level.

Identification of two frataxin isoforms in Zea mays: Structural and functional studies

Frataxin is a ubiquitous protein that plays a role in Fe-S cluster biosynthesis and iron and heme metabolism, although its molecular functions are not entirely clear. In non-photosynthetic eukaryotes, frataxin is encoded by a single gene, and the protein localizes to mitochondria. Here we report the presence of two functional frataxin isoforms in Zea mays, ZmFH-1 and ZmFH-2. We confirmed our previous findings regarding plant frataxins: both proteins have dual localization in mitochondria and chloroplasts. Physiological, biochemical and biophysical studies show some differences in the expression pattern, protection against oxidants and in the aggregation state of both isoforms, suggesting that the two frataxin homologs would play similar but not identical roles in plant cell metabolism. In addition, two specific features of plant frataxins were evidenced: their ability to form dimers and their tendency to undergo conformational change under oxygen exposure.