Discovery of small molecule CD40-TRAF6 inhibitors.

The CD154-CD40 receptor complex plays a pivotal role in several inflammatory pathways. Attempts to inhibit the formation of this complex have resulted in systemic side effects. Downstream inhibition of the CD40 signaling pathway therefore seems a better way to ameliorate inflammatory disease. To relay a signal, the CD40 receptor recruits adapter proteins called tumor necrosis factor receptor-associated factors (TRAFs). CD40-TRAF6 interactions are known to play an essential role in several inflammatory diseases. We used in silico, in vitro, and in vivo experiments to identify and characterize compounds that block CD40-TRAF6 interactions. We present in detail our drug docking and optimization pipeline and show how we used it to find lead compounds that reduce inflammation in models of peritonitis and sepsis. These compounds appear to be good leads for drug development, given the observed absence of side effects and their demonstrated efficacy for peritonitis and sepsis in mouse models.

A complex V ATP5A1 defect causes fatal neonatal mitochondrial encephalopathy.

Whole exome sequencing is a powerful tool to detect novel pathogenic mutations in patients with suspected mitochondrial disease. However, the interpretation of novel genetic variants is not always straightforward. Here, we present two siblings with a severe neonatal encephalopathy caused by complex V deficiency. The aim of this study was to uncover the underlying genetic defect using the combination of enzymatic testing and whole exome sequence analysis, and to provide evidence for causality by functional follow-up. Measurement of the oxygen consumption rate and enzyme analysis in fibroblasts were performed. Immunoblotting techniques were applied to study complex V assembly. The coding regions of the genome were analysed. Three-dimensional modelling was applied. Exome sequencing of the two siblings with complex V deficiency revealed a heterozygous mutation in the ATP5A1 gene, coding for complex V subunit α. The father carried the variant heterozygously. At the messenger RNA level, only the mutated allele was expressed in the patients, whereas the father expressed both the wild-type and the mutant allele. Gene expression data indicate that the maternal allele is not expressed, which is supported by the observation that the ATP5A1 expression levels in the patients and their mother are reduced to ∼50%. Complementation with wild-type ATP5A1 restored complex V in the patient fibroblasts, confirming pathogenicity of the defect. At the protein level, the mutation results in a disturbed interaction of the α-subunit with the ß-subunit of complex V, which interferes with the stability of the complex. This study demonstrates the important value of functional studies in the diagnostic work-up of mitochondrial patients, in order to guide genetic variant prioritization, and to validate gene defects.

Evolution and diversification of the organellar release factor family.

Translation termination is accomplished by proteins of the Class I release factor family (RF) that recognize stop codons and catalyze the ribosomal release of the newly synthesized peptide. Bacteria have two canonical RFs: RF1 recognizes UAA and UAG, RF2 recognizes UAA and UGA. Despite that these two release factor proteins are sufficient for de facto translation termination, the eukaryotic organellar RF protein family, which has evolved from bacterial release factors, has expanded considerably, comprising multiple subfamilies, most of which have not been functionally characterized or formally classified. Here, we integrate multiple sources of information to analyze the remarkable differentiation of the RF family among organelles. We document the origin, phylogenetic distribution and sequence structure features of the mitochondrial and plastidial release factors: mtRF1a, mtRF1, mtRF2a, mtRF2b, mtRF2c, ICT1, C12orf65, pRF1, and pRF2, and review published relevant experimental data. The canonical release factors (mtRF1a, mtRF2a, pRF1, and pRF2) and ICT1 are derived from bacterial ancestors, whereas the others have resulted from gene duplications of another release factor. These new RF family members have all lost one or more specific motifs relevant for bona fide release factor function but are mostly targeted to the same organelle as their ancestor. We also characterize the subset of canonical release factor proteins that bear nonclassical PxT/SPF tripeptide motifs and provide a molecular-model-based rationale for their retained ability to recognize stop codons. Finally, we analyze the coevolution of canonical RFs with the organellar genetic code. Although the RF presence in an organelle and its stop codon usage tend to coevolve, we find three taxa that encode an RF2 without using UGA stop codons, and one reverse scenario, where mamiellales green algae use UGA stop codons in their mitochondria without having a mitochondrial type RF2. For the latter, we put forward a "stop-codon reinvention" hypothesis that involves the retargeting of the plastid release factor to the mitochondrion.

PRPS1 mutations: four distinct syndromes and potential treatment.

Phosphoribosylpyrophosphate synthetases (PRSs) catalyze the first step of nucleotide synthesis. Nucleotides are central to cell function, being the building blocks of nucleic acids and serving as cofactors in cellular signaling and metabolism. With this in mind, it is remarkable that mutations in phosphoribosylpyrophosphate synthetase 1 (PRPS1), which is the most ubiquitously expressed gene of the three PRS genes, are compatible with life. Mutations described thus far in PRPS1 are all missense mutations that result in PRS-I superactivity or in variable levels of decreased activity, resulting in X-linked Charcot-Marie-Tooth disease-5 (CMTX5), Arts syndrome, and X-linked nonsyndromic sensorineural deafness (DFN2). Patients with PRS-I superactivity primarily present with uric acid overproduction, mental retardation, ataxia, hypotonia, and hearing impairment. Postlingual progressive hearing loss is found as an isolated feature in DFN2 patients. Patients with CMTX5 and Arts syndrome have peripheral neuropathy, including hearing impairment and optic atrophy. However, patients with Arts syndrome are more severely affected because they also have central neuropathy and an impaired immune system. The neurological phenotype in all four PRPS1-related disorders seems to result primarily from reduced levels of GTP and possibly other purine nucleotides including ATP, suggesting that these disorders belong to the same disease spectrum. Preliminary results of S-adenosylmethionine (SAM) supplementation in two Arts syndrome patients show improvement of their condition, indicating that SAM supplementation in the diet could alleviate some of the symptoms of patients with PRPS1 spectrum diseases by replenishing purine nucleotides (J.C., unpublished data).

Homozygosity mapping reveals PDE6C mutations in patients with early-onset cone photoreceptor disorders.

Cone photoreceptor disorders form a clinical spectrum of diseases that include progressive cone dystrophy (CD) and complete and incomplete achromatopsia (ACHM). The underlying disease mechanisms of autosomal recessive (ar)CD are largely unknown. Our aim was to identify causative genes for these disorders by genome-wide homozygosity mapping. We investigated 75 ACHM, 97 arCD, and 20 early-onset arCD probands and excluded the involvement of known genes for ACHM and arCD. Subsequently, we performed high-resolution SNP analysis and identified large homozygous regions spanning the PDE6C gene in one sibling pair with early-onset arCD and one sibling pair with incomplete ACHM. The PDE6C gene encodes the cone alpha subunit of cyclic guanosine monophosphate (cGMP) phosphodiesterase, which converts cGMP to 5'-GMP, and thereby plays an essential role in cone phototransduction. Sequence analysis of the coding region of PDE6C revealed homozygous missense mutations (p.R29W, p.Y323N) in both sibling pairs. Sequence analysis of 104 probands with arCD and 10 probands with ACHM revealed compound heterozygous PDE6C mutations in three complete ACHM patients from two families. One patient had a frameshift mutation and a splice defect; the other two had a splice defect and a missense variant (p.M455V). Cross-sectional retinal imaging via optical coherence tomography revealed a more pronounced absence of cone photoreceptors in patients with ACHM compared to patients with early-onset arCD. Our findings identify PDE6C as a gene for cone photoreceptor disorders and show that arCD and ACHM constitute genetically and clinically overlapping phenotypes.

CING: an integrated residue-based structure validation program suite.

We present a suite of programs, named CING for Common Interface for NMR Structure Generation that provides for a residue-based, integrated validation of the structural NMR ensemble in conjunction with the experimental restraints and other input data. External validation programs and new internal validation routines compare the NMR-derived models with empirical data, measured chemical shifts, distance- and dihedral restraints and the results are visualized in a dynamic Web 2.0 report. A red-orange-green score is used for residues and restraints to direct the user to those critiques that warrant further investigation. Overall green scores below ~20 % accompanied by red scores over ~50 % are strongly indicative of poorly modelled structures. The publically accessible, secure iCing webserver ( https://nmr.le.ac.uk ) allows individual users to upload the NMR data and run a CING validation analysis.

Flexible protein-ligand docking using the Fleksy protocol.

Considering protein plasticity is important in accurately predicting the three-dimensional geometry of protein-ligand complexes. Here, we present the first public release of our flexible docking tool Fleksy, which is able to consider both ligand and protein flexibility in the docking process. We describe the workflow and different features of the software and present its performance on two cross-docking benchmark datasets.

Enzyme-specific activation versus leaving group ability.

Enzyme-specific activation and the substrate mimetics strategy are effective ways to circumvent the limited substrate recognition often encountered in protease-catalyzed peptide synthesis. A key structural element in both approaches is the guanidinophenyl (OGp) ester, which enables important interactions for affinity and recognition by the enzyme--at least, this is usually the explanation given for its successful application. In this study we show that leaving group ability is of equal or even greater importance. To this end we used both experimental and computational methods: 1) synthesis of close analogues of OGp, and their evaluation in a dipeptide synthesis assay with trypsin, 2) molecular docking studies to provide insights into the binding mode, and 3) ab initio calculations to evaluate their electronic properties.

Papain-specific activating esters in aqueous dipeptide synthesis.

Enzymatic peptide synthesis has the potential to be a viable alternative for chemical peptide synthesis. Because of the increasing commercial interest in peptides, new and improved enzymatic synthesis methods are desirable. In recently developed enzymatic strategies such as substrate mimetic approaches and enzyme-specific activation, use of the guanidinophenyl ester (OGp) group has been shown to suffer from some drawbacks. OGp esters are sensitive to spontaneous chemical hydrolysis and the group is expensive to synthesize and therefore not suitable for large-scale applications. On the basis of earlier computational studies, we hypothesized that OGp might be replaceable by simpler ester groups to make the enzyme-specific activation approach to peptide bond formation more accessible. To this end, a set of potential activating esters (Z-Gly-Act) was designed, synthesized, and evaluated. Both the benzyl (OBn) and the dimethylaminophenyl (ODmap) esters gave promising results. For these esters, the scope of a model dipeptide synthesis reaction under aqueous conditions was investigated by varying the amino acid donor. The results were compared with those obtained from a previous study of Z-X(AA) -OGp esters. Computational docking analysis of the set of esters was performed in order to provide insight into the differences in the reactivities of all the potential activating esters. Finally, selected ODmap- and OBn-activated amino acids were applied in the synthesis of two biologically active dipeptides on preparative scales.

Gaussian interaction profile kernels for predicting drug-target interaction.

MOTIVATION: The in silico prediction of potential interactions between drugs and target proteins is of core importance for the identification of new drugs or novel targets for existing drugs. However, only a tiny portion of all drug-target pairs in current datasets are experimentally validated interactions. This motivates the need for developing computational methods that predict true interaction pairs with high accuracy. RESULTS: We show that a simple machine learning method that uses the drug-target network as the only source of information is capable of predicting true interaction pairs with high accuracy. Specifically, we introduce interaction profiles of drugs (and of targets) in a network, which are binary vectors specifying the presence or absence of interaction with every target (drug) in that network. We define a kernel on these profiles, called the Gaussian Interaction Profile (GIP) kernel, and use a simple classifier, (kernel) Regularized Least Squares (RLS), for prediction drug-target interactions. We test comparatively the effectiveness of RLS with the GIP kernel on four drug-target interaction networks used in previous studies. The proposed algorithm achieves area under the precision-recall curve (AUPR) up to 92.7, significantly improving over results of state-of-the-art methods. Moreover, we show that using also kernels based on chemical and genomic information further increases accuracy, with a neat improvement on small datasets. These results substantiate the relevance of the network topology (in the form of interaction profiles) as source of information for predicting drug-target interactions. AVAILABILITY: Software and Supplementary Material are available at http://cs.ru.nl/~tvanlaarhoven/drugtarget2011/. CONTACT: tvanlaarhoven@cs.ru.nl; elenam@cs.ru.nl. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Phosphoribosylpyrophosphate synthetase superactivity and recurrent infections is caused by a p.Val142Leu mutation in PRS-I.

We identified a novel missense mutation, c.424G>C (p.Val142Leu) in PRPS1 in a patient with uric acid overproduction without gout but with developmental delay, hypotonia, hearing loss, and recurrent respiratory infections. The uric acid overproduction accompanying this combination of symptoms suggests that the patient presented with phosphoribosylpyrophosphate (PRPP) synthetase superactivity, but recurrent infections have not been associated with superactivity until now. However, recurrent infections are a prominent feature of patients with Arts syndrome, which is caused by PRPS1 loss-of-function mutations, indicating that the patient reported here has an intermediate phenotype. Molecular modeling predicts that the p.Val142Leu change affects both allosteric sites that are involved in inhibition of PRPS1 and the ATP-binding site, which suggests that this substitution can result both in a gain-of-function and loss-of-function of PRPP synthetase. This finding is in line with the normal PRPP synthetase activity in fibroblasts and the absence of activity in erythrocytes of the present patient. We postulate that the overall effect of the p.Val142Leu change on protein activity is determined by the cell type, being a gain-of-function in proliferating cells and a loss-of-function in postmitotic cells. Our results show that missense mutations in PRPS1 can cause a continuous spectrum of features ranging from progressive non-syndromic postlingual hearing impairment to uric acid overproduction, neuropathy, and recurrent infections depending on the functional sites that are affected.

Basal laminar drusen caused by compound heterozygous variants in the CFH gene.

Age-related macular degeneration (AMD) is a multifactorial disease that is strongly associated with the Tyr402His variant in the complement factor H (CFH) gene. Drusen are hallmark lesions of AMD and consist of focal-inflammatory and/or immune-mediated depositions of extracellular material at the interface of the retinal pigment epithelium (RPE) and the Bruch membrane. We evaluated the role of CFH in 30 probands with early-onset drusen and identified heterozygous nonsense, missense, and splice variants in five families. The affected individuals all carried the Tyr402His AMD risk variant on the other allele. This supports an autosomal-recessive disease model in which individuals who carry a CFH mutation on one allele and the Tyr402His variant on the other allele develop drusen. Our findings strongly suggest that monogenic inheritance of CFH variants can result in basal laminar drusen in young adults, and this can progress to maculopathy and severe vision loss later in life.

Papain-catalyzed peptide bond formation: enzyme-specific activation with guanidinophenyl esters.

The substrate mimetics approach is a versatile method for small-scale enzymatic peptide-bond synthesis in aqueous systems. The protease-recognized amino acid side chain is incorporated in an ester leaving group, the substrate mimetic. This shift of the specific moiety enables the acceptance of amino acids and peptide sequences that are normally not recognized by the enzyme. The guanidinophenyl group (OGp), a known substrate mimetic for the serine proteases trypsin and chymotrypsin, has now been applied for the first time in combination with papain, a cheap and commercially available cysteine protease. To provide insight in the binding mode of various Z-X(AA)-OGp esters, computational docking studies were performed. The results strongly point at enzyme-specific activation of the OGp esters in papain through a novel mode of action, rather than their functioning as mimetics. Furthermore, the scope of a model dipeptide synthesis was investigated with respect to both the amino acid donor and the nucleophile. Molecular dynamics simulations were carried out to prioritize 22 natural and unnatural amino acid donors for synthesis. Experimental results correlate well with the predicted ranking and show that nearly all amino acids are accepted by papain.

Snooker: a structure-based pharmacophore generation tool applied to class A GPCRs.

G-protein coupled receptors (GPCRs) are important drug targets for various diseases and of major interest to pharmaceutical companies. The function of individual members of this protein family can be modulated by the binding of small molecules at the extracellular side of the structurally conserved transmembrane (TM) domain. Here, we present Snooker, a structure-based approach to generate pharmacophore hypotheses for compounds binding to this extracellular side of the TM domain. Snooker does not require knowledge of ligands, is therefore suitable for apo-proteins, and can be applied to all receptors of the GPCR protein family. The method comprises the construction of a homology model of the TM domains and prioritization of residues on the probability of being ligand binding. Subsequently, protein properties are converted to ligand space, and pharmacophore features are generated at positions where protein ligand interactions are likely. Using this semiautomated knowledge-driven bioinformatics approach we have created pharmacophore hypotheses for 15 different GPCRs from several different subfamilies. For the beta-2-adrenergic receptor we show that ligand poses predicted by Snooker pharmacophore hypotheses reproduce literature supported binding modes for ∼75% of compounds fulfilling pharmacophore constraints. All 15 pharmacophore hypotheses represent interactions with essential residues for ligand binding as observed in mutagenesis experiments and compound selections based on these hypotheses are shown to be target specific. For 8 out of 15 targets enrichment factors above 10-fold are observed in the top 0.5% ranked compounds in a virtual screen. Additionally, prospectively predicted ligand binding poses in the human dopamine D3 receptor based on Snooker pharmacophores were ranked among the best models in the community wide GPCR dock 2010.

Stress and innate immunity in carp: corticosteroid receptors and pro-inflammatory cytokines.

The stress hormone cortisol is deeply involved in immune regulation in all vertebrates. Common carp (Cyprinus carpio L.) express four corticoid receptors that may modulate immune responses: three glucocorticoid receptors (GR); GR1, with two splice variants (GR1a and GR1b), GR2 and a single mineralocorticoid receptor (MR). All receptors are expressed as of 4 days post-fertilization and may thus play a critical role in development and functioning of the adult immune system. Immune tissues and cells predominantly express mRNA for GRs compared to mRNA for the MR. Three-dimensional protein structure modeling predicts, and transfection assays confirm that alternative splicing of GR1 does not influence the capacity to induce transcription of effector genes. When tested for cortisol activation, GR2 is the most sensitive corticoid receptor in carp, followed by the MR and GR1a and GR1b. Lipopolysacharide (LPS) treatment of head kidney phagocytes quickly induces GR1 expression and inhibits GR2 expression. Cortisol treatment in vivo enhances GR1a and MR mRNA expression, but only mildly, and cortisol treatment in vitro does not affect receptor expression of phagocytes. Cortisol has no direct effect on the LPS-induced receptor profile. Therefore, an immune rather than a stress stimulus regulates GR expression. Cortisol administered at stress levels to phagocytes in vitro significantly inhibits LPS-induced expression of the pro-inflammatory cytokines tumor necrosis factor alpha (TNF-alpha) and interleukin-12 (IL-12) (subunit p35) and of inducible nitric oxide synthase (iNOS) expression. A physiologically differential function for GR1 and GR2 in the immune response of fish to infection is indicated.

Features of a twin-arginine signal peptide required for recognition by a Tat proofreading chaperone.

The twin-arginine translocation (Tat) system is a bacterial protein targeting pathway. Tat-targeted proteins display signal peptides containing a distinctive SRRxFLK 'twin-arginine' motif. The Escherichia coli trimethylamine N-oxide reductase (TorA) bears a bifunctional Tat signal peptide, which directs protein export and serves as a binding site for the TorD biosynthetic chaperone. Here, the physical interaction between TorD and the TorA signal peptide was investigated. A single substitution within the TorA signal peptide (L31Q) was sufficient to impair TorD binding. Screening of a random torD mutant library identified a variant TorD protein (Q7L) that displayed increased binding affinity for the TorA signal peptide.

A flexible approach to induced fit docking.

We present Fleksy, a new approach to consider both ligand and receptor flexibility in small molecule docking. Pivotal to our method is the use of a receptor ensemble to describe protein flexibility. To construct these ensembles, we use a backbone-dependent rotamer library and implement the concept of interaction sampling. The latter allows the evaluation of different orientations of ambivalent interaction partners. The docking stage consists of an ensemble-based soft-docking experiment using FlexX-Ensemble, followed by an effective flexible receptor-ligand complex optimization using Yasara. Fleksy produces a set of receptor-ligand complexes ranked using a consensus scoring function combining docking scores and force field energies. Averaged over three cross-docking datasets, containing 35 different receptor-ligand complexes in total, Fleksy reproduces the observed binding mode within 2.0 A for 78% of the complexes. This compares favorably to the rigid receptor FlexX program, which on average reaches a success rate of 44% for these datasets.

Traditional biomolecular structure determination by NMR spectroscopy allows for major errors.

One of the major goals of structural genomics projects is to determine the three-dimensional structure of representative members of as many different fold families as possible. Comparative modeling is expected to fill the remaining gaps by providing structural models of homologs of the experimentally determined proteins. However, for such an approach to be successful it is essential that the quality of the experimentally determined structures is adequate. In an attempt to build a homology model for the protein dynein light chain 2A (DLC2A) we found two potential templates, both experimentally determined nuclear magnetic resonance (NMR) structures originating from structural genomics efforts. Despite their high sequence identity (96%), the folds of the two structures are markedly different. This urged us to perform in-depth analyses of both structure ensembles and the deposited experimental data, the results of which clearly identify one of the two models as largely incorrect. Next, we analyzed the quality of a large set of recent NMR-derived structure ensembles originating from both structural genomics projects and individual structure determination groups. Unfortunately, a visual inspection of structures exhibiting lower quality scores than DLC2A reveals that the seriously flawed DLC2A structure is not an isolated incident. Overall, our results illustrate that the quality of NMR structures cannot be reliably evaluated using only traditional experimental input data and overall quality indicators as a reference and clearly demonstrate the urgent need for a tight integration of more sophisticated structure validation tools in NMR structure determination projects. In contrast to common methodologies where structures are typically evaluated as a whole, such tools should preferentially operate on a per-residue basis.

The presence of multiple and differentially regulated interleukin-12p40 genes in bony fishes signifies an expansion of the vertebrate heterodimeric cytokine family.

Interleukin-12 (IL-12) is the founding member of a rapidly growing family of heterodimeric cytokines. It consists of two subunits, designated p35 and p40 that together constitute a disulphide-linked heterodimeric cytokine. IL-12 is well known for its prominent role in both the early innate immune response and the skewing of the ensuing acquired immune response towards Th1. Here, we report the presence of IL-12p35 and three highly distinct IL-12p40 genes in common carp (Cyprinus carpio). The carp is a bony fish species genetically similar to the zebrafish, but its substantially larger body size facilitates immunological studies. A comparison of IL-12p35 genes of mammalian and non-mammalian species reveals the presence of a duplicated exon that is unique to the mammalian lineage. The organisation of the three carp IL-12p40 genes is similar to that of higher vertebrates. Phylogenetic analyses that include the p40-related subunits of other composite cytokines confirm the presence of three genuine IL-12p40 genes in carp and indicate that they are evolutionary ancient and possibly not restricted to bony fishes. The orthology of the different carp p40 subunits to mammalian IL-12p40 is further evident from the conservation of key residues involved in the formation of intra- and interchain disulphide bridges and the tight interlocking topology between p35 and p40. The expression of each of the carp IL-12p40 genes differs profoundly, constitutively as well as in response to in vitro stimulation of carp macrophages. Collectively, the presence of multiple and substantially different IL-12 genes signifies a considerable expansion of the vertebrate heterodimeric cytokine family.

Increased leptin expression in common Carp (Cyprinus carpio) after food intake but not after fasting or feeding to satiation.

Leptin is a key factor in the regulation of food intake and is an important factor in the pathophysiology of obesity. However, more than a decade after the discovery of leptin in mouse, information regarding leptin in any nonmammalian species is still scant. We report the identification of duplicate leptin genes in common carp (Cyprinus carpio). The unique gene structure, the conservation of both cysteines that form leptin's single disulfide bridge, and stable clustering in phylogenetic analyses substantiate the unambiguous orthology of mammalian and carp leptins, despite low amino acid identity. The liver is a major yet not the only site of leptin expression. However, neither 6 d nor 6 wk of fasting nor subsequent refeeding affected hepatic leptin expression, although the carp predictably shifted from carbohydrate to lipid metabolism. Animals that were fed to satiation grew twice as fast as controls; however, they did not show increased leptin expression at the termination of the study. Hepatic leptin expression did, however, display an acute and transient postprandial increase that follows the postprandial plasma glucose peak. In summary, leptin mRNA expression in carp changes acutely after food intake, but involvement of leptin in the long-term regulation of food intake and energy metabolism was not evident from fasting for days or weeks or long-term feeding to satiation. These are the first data on the regulation of leptin expression in any nonmammalian species.