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.

Abstract 611: Blocking CD40-TRAF6 Signaling is a Novel Therapeutic Target in Obesity-Associated Insulin Resistance

The immune system plays an instrumental role in obesity and insulin resistance. Here we unravel the role of the co-stimulatory molecule, CD40, and its signaling intermediates, TNF-Receptor-Associated-Factors (TRAFs), in diet-induced obesity (DIO). Although not exhibiting increased weight gain, male CD40-/- mice in DIO displayed worsened insulin resistance, as compared to wild type mice. This was associated with excessive inflammation of adipose tissue (AT), characterized by increased accumulation of CD8+ T cells and M1 macrophages, and enhanced hepatosteatosis. Mice with deficient CD40-TRAF2/3/5 signaling in MHCII+ cells exhibited a similar phenotype in DIO as CD40-/- mice. In contrast, mice with deficient CD40-TRAF6 signaling in MHCII+ cells displayed no insulin resistance, and showed a reduction in both AT inflammation and hepatosteatosis in DIO. To prove the therapeutic potential of inhibition of CD40-TRAF6 in obesity, DIO mice were treated with a small-molecule inhibitor that we designed to specifically block CD40-TRAF6 interactions; this improved insulin sensitivity, reduced AT inflammation and decreased hepatosteatosis. Our study reveals that the CD40-TRAF2/3/5 signaling pathway in MHCII+ cells protects against AT inflammation and metabolic complications associated with obesity, whereas CD40-TRAF6 interactions in MHCII+ cells aggravate these complications. Inhibition of CD40-TRAF6 signaling by our newly developed compound may provide a novel therapeutic option in obesity-associated insulin resistance.

Collapse

Application of a cell-based protease assay for testing inhibitors of picornavirus 3C proteases

Proteolytical cleavage of the picornaviral polyprotein is essential for viral replication. Therefore, viral proteases are attractive targets for anti-viral therapy. Most assays available for testing proteolytical activity of proteases are performed in vitro, using heterologously expressed proteases and peptide substrates. To deal with the disadvantages associated with in vitro assays, we modified a cell-based protease assay for picornavirus proteases. The assay is based on the induction of expression of a firefly luciferase reporter by a chimeric transcription factor in which the viral protease and cleavage sites are inserted between the GAL4 binding domain and the VP16 activation domain. Firefly luciferase expression is dependent on cleavage of the transcription factor by the viral protease. This biosafe assay enables testing the effect of compounds on protease activity in cells while circumventing the need for infection. We designed the assay for 3C proteases (3C(pro)) of various enteroviruses as well as of viruses of several other picornavirus genera, and show that the assay is amenable for use in a high-throughput setting. Furthermore, we show that the spectrum of activity of 3C(pro) inhibitor AG7088 (rupintrivir) not only encompasses enterovirus 3C(pro) but also 3C(pro) of foot-and-mouth disease virus (FMDV), an aphthovirus. In contrary, AG7404 (compound 1), an analogue of AG7088, had no effect on FMDV 3C(pro) activity, for which we provide a structural explanation.

X-linked Charcot-Marie-Tooth disease, Arts syndrome, and prelingual non-syndromic deafness form a disease continuum: evidence from a family with a novel PRPS1 mutation

Background

X-linked Charcot-Marie-Tooth disease type 5 (CMTX5), Arts syndrome, and non-syndromic sensorineural deafness (DFN2) are allelic syndromes, caused by reduced activity of phosphoribosylpyrophosphate synthetase 1 (PRS-I) due to loss-of-function mutations in PRPS1. As only few families have been described, knowledge about the relation between these syndromes, the phenotypic spectrum in patients and female carriers, and the relation to underlying PRS-I activity is limited.

Methods

We investigated a family with a novel PRPS1 mutation (c.830A > C, p.Gln277Pro) by extensive phenotyping, MRI, and genetic and enzymatic tests.

Results

The male index subject presented with an overlap of CMTX5 and Arts syndrome features, whereas his sister presented with prelingual DFN2. Both showed mild parietal and cerebellar atrophy on MRI. Enzymatically, PRS-I activity was undetectable in the index subject, reduced in his less affected sister, and normal in his unaffected mother.

Conclusions

Our findings demonstrate that CMTX5, Arts syndrome and DFN2 are phenotypic clusters on an intrafamilial continuum, including overlapping phenotypes even within individuals. The respective phenotypic presentation seems to be determined by the exact PRPS1 mutation and the residual enzyme activity, the latter being largely influenced by the degree of skewed X-inactivation. Finally, our findings show that brain atrophy might be more common in PRPS1-disorders than previously thought.

A 3-base pair deletion, c.9711_9713del, in DMD results in intellectual disability without muscular dystrophy

We have identified a deletion of 3 base pairs in the dystrophin gene (DMD), c.9711_9713del, in a family with nonspecific X-linked intellectual disability (ID) by sequencing of the exons of 86 known X-linked ID genes. This in-frame deletion results in the deletion of a single-amino-acid residue, Leu3238, in the brain-specific isoform Dp71 of dystrophin. Linkage analysis supported causality as the mutation was present in the 7.6 cM linkage interval on Xp22.11-Xp21.1 with a maximum positive LOD score of 2.41 (MRX85 locus). Molecular modeling predicts that the p.(Leu3238del) deletion results in the destabilization of the C-terminal domain of dystrophin and hence reduces the ability to interact with β-dystroglycan. Correspondingly, Dp71 protein levels in lymphoblastoid cells from the index patient are 6.7-fold lower than those in control cell lines (P=0.08). Subsequent determination of the creatine kinase levels in blood of the index patient showed a mild but significant elevation in serum creatine kinase, which is in line with impaired dystrophin function. In conclusion, we have identified the first DMD mutation in Dp71 that results in ID without muscular dystrophy.

Abstract 14: Small Molecule Inhibitors of the CD40-TRAF6 Interaction Reduce Atherosclerosis by Inducing Hypo-inflammatory Myeloid Cells

The co-stimulatory CD40-CD40L dyad is well known for its pro-inflammatory role in atherosclerosis. Inhibition of CD40(L) in hyperlipidemic mice drastically reduces atherosclerosis. However, long term inhibition of CD40(L) results in immunosuppression and/or thromboembolic events. Therefore more targeted inhibition of the CD40L-CD40 dyad is required. To elicit intracellular signalling upon activation, CD40 needs to recruit adaptor proteins: the tumour necrosis factor receptor-associated factors (TRAFs). Using mice with a mutation in the CD40-TRAF binding site, we previously showed that CD40-TRAF6, and not CD40-TRAF2/3/5, are pivotal in atherosclerosis.

To discover small drug-like molecules that inhibit the CD40-TRAF6 interaction an in silico approach of virtual ligand screening was used; this resulted in the identification of 40.000 compounds. Analysis of the top 800 compounds identified 7 compounds that dose-dependently reduced NFκB activation and IL1β and IL6 expression in RAW cells and CD40-stimulated bone marrow-derived macrophages, respectively. Surface plasmon resonance experiments confirmed direct binding of the compounds to the TRAF6 C-domain.

To analyze the effects of the top two compounds on atherosclerosis, Apoe -/- mice were daily treated with compounds 6877002 or 6860766 for 6 weeks. Compound treatment reduced atherosclerosis in the aortic arch by 47.1% (6877002) and 66.8% (6860766). The number of plaque monocytes/macrophages was decreased by 41.4% (6877002) and 53.0% (6860766), and granulocytes by 62.5% (6877002) and 41.5% (6860766). Intravital microscopy showed that the compounds reduced monocyte recruitment to the endothelium by 40.1% (6877002) and 51.2% (6860766), neutrophil adhesion was reduced by 40.2% (6877002) and 51.2% (6860766). In accordance, expression of CCL2-CCR2 and CCL5-CCR5 was remarkably reduced in compound treated macrophages. Moreover, compound treatment abolished CD40-induced expression of TNFα, IL1β, IL6, IL10, and IL12.

These results establish the importance and possibilities of long-term therapeutic inhibition of CD40-TRAF6 interactions in atherosclerosis and other inflammatory diseases and harbor the potential to overcome the current limitations of CD40(L) blocking therapies.

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.

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.

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.

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.

In silico identification of potential cholestasis-inducing agents via modeling of Na(+)-dependent taurocholate cotransporting polypeptide substrate specificity

Na(+)-dependent taurocholate cotransporting polypeptide (NTCP, SLC10A1) is the main transporter facilitating the hepatic uptake of bile acids from the circulation. Consequently, the interaction of xenobiotics, including therapeutic drugs, with the bile acid binding pocket of NTCP could lead to impairment of hepatic bile acid uptake. We pursued a 3D-pharmacophore approach to model the NTCP substrate and inhibitor specificity and investigated whether it is possible to identify compounds with intrinsic NTCP inhibitory properties. Based on known endogenous NTCP substrates, a 3D-pharmacophore model was built, which was subsequently used to screen two virtual libraries together containing the structures of 10 million compounds. Studies with Chinese hamster ovary cells overexpressing human NTCP, human hepatocytes, ex vivo perfused rat livers, and bile duct-cannulated rats were conducted to validate the activity of the virtual screening hits. Modeling yielded a 3D-pharmacophore, consisting of two hydrogen bond acceptors and three hydrophobic features. Six out of 10 structurally diverse compounds selected in the first virtual screening procedure significantly inhibited taurocholate uptake in the NTCP overexpressing cells. For the most potent inhibitor identified, an anthraquinone derivative, this finding was confirmed in human hepatocytes and perfused rat livers. Subsequent structure and activity relationship studies with analogs of this derivative indicated that an appropriate distance between hydrogen bond acceptor features and presence of one or two negative charges appear critical for a successful NTCP interaction. In conclusion, pharmacophore modeling was successfully used to identify compounds that inhibit NTCP. Our approach represents an important first step toward the in silico flagging of potential cholestasis-inducing molecules.

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.

Structure based hypothesis of a mitochondrial ribosome rescue mechanism

Background

mtRF1 is a vertebrate mitochondrial protein with an unknown function that arose from a duplication of the mitochondrial release factor mtRF1a. To elucidate the function of mtRF1, we determined the positions that are conserved among mtRF1 sequences but that are different in their mtRF1a paralogs. We subsequently modeled the 3D structure of mtRF1a and mtRF1 bound to the ribosome, highlighting the structural implications of these differences to derive a hypothesis for the function of mtRF1.

Results

Our model predicts, in agreement with the experimental data, that the 3D structure of mtRF1a allows it to recognize the stop codons UAA and UAG in the A-site of the ribosome. In contrast, we show that mtRF1 likely can only bind the ribosome when the A-site is devoid of mRNA. Furthermore, while mtRF1a will adopt its catalytic conformation, in which it functions as a peptidyl-tRNA hydrolase in the ribosome, only upon binding of a stop codon in the A-site, mtRF1 appears specifically adapted to assume this extended, peptidyl-tRNA hydrolyzing conformation in the absence of mRNA in the A-site.

Conclusions

We predict that mtRF1 specifically recognizes ribosomes with an empty A-site and is able to function as a peptidyl-tRNA hydrolase in those situations. Stalled ribosomes with empty A-sites that still contain a tRNA bound to a peptide chain can result from the translation of truncated, stop-codon less mRNAs. We hypothesize that mtRF1 recycles such stalled ribosomes, performing a function that is analogous to that of tmRNA in bacteria.

Reviewers

This article was reviewed by Dr. Eugene Koonin, Prof. Knud H. Nierhaus (nominated by Dr. Sarah Teichmann) and Dr. Shamil Sunyaev.

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.

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.

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.

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.

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.

Two divergent leptin paralogues in zebrafish (Danio rerio) that originate early in teleostean evolution

We describe duplicate leptin genes in zebrafish (Danio rerio) that share merely 24% amino acid identity with each other and only 18% with human leptin. We were also able to retrieve a second leptin gene in medaka (Oryzias latipes). The presence of duplicate leptin genes in these two distantly related teleosts suggests that duplicate leptin genes are a common feature of teleostean fishes. Despite low primary sequence conservation, we are confident in assigning orthology between mammalian and zebrafish leptins for several reasons. First, both zebrafish leptins share their characteristic gene structure and display key features of conserved synteny with mammalian leptin genes. Secondly, the cysteine residues that make up leptin's single disulphide bridge are equally spaced in mammalian and zebrafish leptins and are unique among all members of the class-I helical cytokine family. Thirdly, the zebrafish leptins cluster with other fish leptins and mammalian leptins in phylogenetic analysis, supported by high bootstrap values. Within the leptin cluster, leptin-b forms a separate clade with the leptin-b orthologue from medaka. Finally, our prediction of the tertiary structures shows that both leptins conform to the typical four alpha-helix bundle structure of the class-I alpha-helical cytokines. The zebrafish leptins are differentially expressed; the liver shows high leptin-a expression (in concordance with what we observed for carp leptins), while leptin-b is expressed at much lower levels, which are downregulated further upon fasting. The finding of duplicate leptin genes in teleosts adds to our understanding of the evolution of leptin physiology in the early vertebrate lineage.

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.

Arts syndrome is caused by loss-of-function mutations in PRPS1

Arts syndrome is an X-linked disorder characterized by mental retardation, early-onset hypotonia, ataxia, delayed motor development, hearing impairment, and optic atrophy. Linkage analysis in a Dutch family and an Australian family suggested that the candidate gene maps to Xq22.1-q24. Oligonucleotide microarray expression profiling of fibroblasts from two probands of the Dutch family revealed reduced expression levels of the phosphoribosyl pyrophosphate synthetase 1 gene (PRPS1). Subsequent sequencing of PRPS1 led to the identification of two different missense mutations, c.455T-->C (p.L152P) in the Dutch family and c.398A-->C (p.Q133P) in the Australian family. Both mutations result in a loss of phosphoribosyl pyrophosphate synthetase 1 activity, as was shown in silico by molecular modeling and was shown in vitro by phosphoribosyl pyrophosphate synthetase activity assays in erythrocytes and fibroblasts from patients. This is in contrast to the gain-of-function mutations in PRPS1 that were identified previously in PRPS-related gout. The loss-of-function mutations of PRPS1 likely result in impaired purine biosynthesis, which is supported by the undetectable hypoxanthine in urine and the reduced uric acid levels in serum from patients. To replenish low levels of purines, treatment with S-adenosylmethionine theoretically could have therapeutic efficacy, and a clinical trial involving the two affected Australian brothers is currently underway.

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.

Negative Constraints Underlie the ErbB Specificity of Epidermal Growth Factor-like Ligands

Epidermal growth factor (EGF)-like growth factors bind their ErbB receptors in a highly selective manner, but the molecular basis for this specificity is poorly understood. We have previously shown that certain residues in human EGF (Ser(2)-Asp(3)) and TGFalpha (Glu(26)) are not essential for their binding to ErbB1 but prevent binding to ErbB3 and ErbB4. In the present study, we have used a phage display approach to affinity-optimize the C-terminal linear region of EGF-like growth factors for binding to each ErbB receptor and thereby shown that Arg(45) in EGF impairs binding to both ErbB3 and ErbB4. By omitting all these so-called negative constraints from EGF, we designed a ligand designated panerbin that binds ErbB1, ErbB3, and ErbB4 with similarly high affinity as their wild-type ligands. Homology models, based on the known crystal structure of TGFalpha-bound ErbB1, showed that panerbin is able to bind ErbB1, ErbB3, and ErbB4 in a highly similar manner with respect to position and number of interaction sites. Upon in silico introduction of the experimentally known negative constraints into panerbin, we found that Arg(45) induced local charge repulsion and Glu(26) induced steric hindrance in a receptor-specific manner, whereas Ser(2)-Asp(3) impaired binding due to a disordered conformation. Furthermore, radiolabeled panerbin was used to quantify the level of all three receptors on human breast cancer cells in a single radioreceptor assay. It is concluded that the ErbB specificity of EGF-like growth factors primarily results from the presence of a limited number of residues that impair the unintended interaction with other ErbB receptors.

Making optimal use of empirical energy functions: force-field parameterization in crystal space

Today's energy functions are not able yet to distinguish reliably between correct and almost correct protein models. Improving these near-native models is currently a major bottle-neck in homology modeling or experimental structure determination at low resolution. Increasingly accurate energy functions are required to complete the "last mile of the protein folding problem," for example during a molecular dynamics simulation. We present a new approach to reach this goal. For 50 high resolution X-ray structures, the complete unit cell was reconstructed, including disordered water molecules, counter ions, and hydrogen atoms. Simulations were then run at the pH at which the crystal was solved, while force-field parameters were iteratively adjusted so that the damage done to the structures was minimal. Starting with initial parameters from the AMBER force field, the optimization procedure converged at a new force field called YAMBER (Yet Another Model Building and Energy Refinement force field), which is shown to do significantly less damage to X-ray structures, often move homology models in the right direction, and occasionally make them look like experimental structures. Application of YAMBER during the CASP5 structure prediction experiment yielded a model for target 176 that was ranked first among 150 submissions. Due to its compatibility with the well-established AMBER format, YAMBER can be used by almost any molecular dynamics program. The parameters are freely available from www.yasara.org/yamber.

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.

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.

Three-dimensional biomolecular structures provide an invaluable source of biologically relevant information. To be able to learn the most of the wealth of information that these structures can provide us, it is of great importance that the quality and accuracy of the protein structure models deposited in the Protein Data Bank are as high as possible. In this work, the authors describe an analysis that illustrates that this is unfortunately not the case for many protein structures solved using nuclear magnetic resonance spectroscopy. They present an example in which two strikingly different models describing the same protein are analyzed using commonly available structure validation tools, and the results of this analysis show one of the two models to be incorrect. Subsequently, using a large set of recently determined structures, the authors demonstrate that unfortunately this example does not stand on its own. The analyses and examples clearly illustrate that relying solely on the experimental data to evaluate structural quality can provide a false sense of correctness and the combination of multiple sophisticated structure validation tools is required to detect the presence of errors in protein nuclear magnetic resonance structures.

ZNF674: a new kruppel-associated box-containing zinc-finger gene involved in nonsyndromic X-linked mental retardation

Array-based comparative genomic hybridization has proven to be successful in the identification of genetic defects in disorders involving mental retardation. Here, we studied a patient with learning disabilities, retinal dystrophy, and short stature. The family history was suggestive of an X-linked contiguous gene syndrome. Hybridization of full-coverage X-chromosomal bacterial artificial chromosome arrays revealed a deletion of ~1 Mb in Xp11.3, which harbors RP2, SLC9A7, CHST7, and two hypothetical zinc-finger genes, ZNF673 and ZNF674. These genes were analyzed in 28 families with nonsyndromic X-linked mental retardation (XLMR) that show linkage to Xp11.3; the analysis revealed a nonsense mutation, p.E118X, in the coding sequence of ZNF674 in one family. This mutation is predicted to result in a truncated protein containing the Kruppel-associated box domains but lacking the zinc-finger domains, which are crucial for DNA binding. We characterized the complete ZNF674 gene structure and subsequently tested an additional 306 patients with XLMR for mutations by direct sequencing. Two amino acid substitutions, p.T343M and p.P412L, were identified that were not found in unaffected individuals. The proline at position 412 is conserved between species and is predicted by molecular modeling to reduce the DNA-binding properties of ZNF674. The p.T343M transition is probably a polymorphism, because the homologous ZNF674 gene in chimpanzee has a methionine at that position. ZNF674 belongs to a cluster of seven highly related zinc-finger genes in Xp11, two of which (ZNF41 and ZNF81) were implicated previously in XLMR. Identification of ZNF674 as the third XLMR gene in this cluster may indicate a common role for these zinc-finger genes that is crucial to human cognitive functioning.

A novel D458V mutation in the SANS PDZ binding motif causes atypical Usher syndrome

Homozygosity mapping and linkage analysis in a Turkish family with autosomal recessive prelingual sensorineural hearing loss revealed a 15-cM critical region at 17q25.1-25.3 flanked by the polymorphic markers D17S1807 and D17S1806. The maximum two-point lod score was 4.07 at theta=0.0 for the marker D17S801. The linkage interval contains the Usher syndrome 1G gene (USH1G) that is mutated in patients with Usher syndrome (USH) type 1g and encodes the SANS protein. Mutation analysis of USH1G led to the identification of a homozygous missense mutation D458V at the -3 position of the PDZ binding motif of SANS. This mutation was also present homozygously in one out of 64 additional families from Turkey with autosomal recessive nonsyndromic hearing loss and heterozygously in one out of 498 control chromosomes. By molecular modeling, we provide evidence that this mutation impairs the interaction of SANS with harmonin. Ophthalmologic examination and vestibular evaluation of patients from both families revealed mild retinitis pigmentosa and normal vestibular function. These results suggest that these patients suffer from atypical USH.

Definition of a new information-based per-residue quality parameter

For biomolecular NMR structures typically only a poor correspondence is observed between statistics derived from the experimental input data and structural quality indicators obtained from the structure ensembles. Here, we investigate the relationship between the amount of available NMR data and structure quality. By generating datasets with a predetermined information content and evaluating the quality of the resulting structure ensembles we show that there is, in contrast to previous findings, a linear relation between the information contained in experimental data and structural quality. From this relation, a new quality parameter is derived that provides direct insight, on a per-residue basis, into the extent to which structural quality is governed by the experimental input data.

Solution structure of the second PDZ domain of the neuronal adaptor X11alpha and its interaction with the C-terminal peptide of the human copper chaperone for superoxide dismutase

Protection against reactive oxygen species is provided by the copper containing enzyme superoxide dismutase 1 (SOD1). The copper chaperone CCS is responsible for copper insertion into apo-SOD1. This role is impaired by an interaction between the second PDZ domain (PDZ2alpha) of the neuronal adaptor protein X11alpha and the third domain of CCS (McLoughlin et al. (2001) J. Biol. Chem., 276, 9303-9307). The solution structure of the PDZ2alpha domain has been determined and the interaction with peptides derived from CCS has been explored. PDZ2alpha binds to the last four amino acids of the CCS protein (PAHL) with a dissociation constant of 91 +/- 2 microM. Peptide variants have been used to map the interaction areas on PDZ2alpha for each amino acid, showing an important role for the C-terminal leucine, in line with canonical PDZ-peptide interactions.

RECOORD: A recalculated coordinate database of 500+ proteins from the PDB using restraints from the BioMagResBank

State-of-the-art methods based on CNS and CYANA were used to recalculate the nuclear magnetic resonance (NMR) solution structures of 500+ proteins for which coordinates and NMR restraints are available from the Protein Data Bank. Curated restraints were obtained from the BioMagResBank FRED database. Although the original NMR structures were determined by various methods, they all were recalculated by CNS and CYANA and refined subsequently by restrained molecular dynamics (CNS) in a hydrated environment. We present an extensive analysis of the results, in terms of various quality indicators generated by PROCHECK and WHAT_CHECK. On average, the quality indicators for packing and Ramachandran appearance moved one standard deviation closer to the mean of the reference database. The structural quality of the recalculated structures is discussed in relation to various parameters, including number of restraints per residue, NOE completeness and positional root mean square deviation (RMSD). Correlations between pairs of these quality indicators were generally low; for example, there is a weak correlation between the number of restraints per residue and the Ramachandran appearance according to WHAT_CHECK (r = 0.31). The set of recalculated coordinates constitutes a unified database of protein structures in which potential user- and software-dependent biases have been kept as small as possible. The database can be used by the structural biology community for further development of calculation protocols, validation tools, structure-based statistical approaches and modeling. The RECOORD database of recalculated structures is publicly available from http://www.ebi.ac.uk/msd/recoord.

DRESS: a database of REfined solution NMR structures

Several studies have shown that biomolecular NMR structures are often of lower quality when compared to crystal structures, and consequently they are often excluded from structural analyses. We present a publicly available database of re-refined NMR structures, exhibiting significantly improved quality. This database (available at http://www.cmbi.kun.nl/dress/) presents a uniformly refined and validated set of structural models that improves the value of these NMR structures as input for experimental and theoretical studies in many fields of research.

A Closed Binding Pocket and Global Destabilization Modify the Binding Properties of an Alternatively Spliced Form of the Second PDZ Domain of PTP-BL

PTP-BL is a large phosphatase that is implicated in cellular processes as diverse as cytokinesis, actin-cytoskeletal rearrangement, and apoptosis. Five PDZ domains mediate its cellular role by binding to the C termini of target proteins, forming multiprotein complexes. The second PDZ domain (PDZ2) binds to the C termini of the tumor suppressor protein APC and the LIM domain-containing protein RIL; however, in one splice variant, PDZ2as, a 5 residue insertion abrogates this binding. The insert causes distinct structural and dynamical changes in the alternatively spliced PDZ2: enlarging the L1 loop between beta2 and beta3, both lengthening and changing the orientation of the alpha2 helix, giving the base of the binding pocket less flexibility to accommodate ligands, and destabilizing the entire domain. These changes render the binding pocket incapable of binding C termini, possibly having implications in the functional role of PTP-BL.

Quantitative Evaluation of Experimental NMR Restraints

Nuclear Overhauser effect (NOE) data are an indispensable source of structural information in biomolecular structure determination by NMR spectroscopy. The number and type of experimental restraints used in the structure calculation and the RMS deviation of the restraints are usually reported. We present a new method for quantifying the information contained in the experimental NMR restraints. The method is based on a description of the structure in distance space and concepts derived from information theory. It allows for an objective description of the amount of available experimental information, which we show to be related to the positional uncertainty of the NMR ensemble. The measure of information presented is not affected by redundancy in the experimental restraints. Using various examples, we show that the method successfully identifies the crucial restraints in a structure determination: those restraints that are both important and unique. Finally, we demonstrate that the method can detect a wider range of redundancy in experimental datasets when compared to currently available methods. Because our method describes the quantitative evaluation of experimental NMR restraints, we propose the acronym QUEEN.

The precision of NMR structure ensembles revisited

Biomolecular structures provide the basis for many studies in research areas such as structure-based drug design and homology modeling. In order to use molecular coordinates it is important that they are reliable in terms of accurate description of the experimental data and in terms of the overall and local geometry. Besides these primary quality criteria an indication is needed for the uncertainty in the atomic coordinates that may arise from the dynamic behavior of the considered molecules as well as from experimental- and computational procedures. In contrast to the crystallographic B-factor, a good measure for the uncertainty in NMR-derived atomic coordinates is still not available. It has become clear in recent years that the widely used atomic Root Mean Square Deviation (RMSD), which is a measure for the precision of the data, overestimates the accuracy of NMR structure ensembles and therefore is a problematic measure for the uncertainty in the atomic coordinates. In this study we report a method that yields a more realistic estimate of the uncertainty in the atomic coordinates by maximizing the RMSD of an ensemble of structures, while maintaining the accordance with the experimentally derived data. The results indicate that the RMSD of most NMR structure ensembles can be significantly increased compromising neither geometric quality nor NMR data. This maximized RMSD therefore seems a better estimate of the true uncertainty in the atomic coordinates.