Merridee A. Wouters

Mutation of the LUNATIC FRINGE Gene in Humans Causes Spondylocostal Dysostosis with a Severe Vertebral Phenotype

The spondylocostal dysostoses (SCDs) are a heterogeneous group of vertebral malsegmentation disorders that arise during embryonic development by a disruption of somitogenesis. Previously, we had identified two genes that cause a subset of autosomal recessive forms of this disease: DLL3 (SCD1) and MESP2 (SCD2). These genes are important components of the Notch signaling pathway, which has multiple roles in development and disease. Here, we have used a candidate-gene approach to identify a mutation in a third Notch pathway gene, LUNATIC FRINGE (LFNG), in a family with autosomal recessive SCD. LFNG encodes a glycosyltransferase that modifies the Notch family of cell-surface receptors, a key step in the regulation of this signaling pathway. A missense mutation was identified in a highly conserved phenylalanine close to the active site of the enzyme. Functional analysis revealed that the mutant LFNG was not localized to the correct compartment of the cell, was unable to modulate Notch signaling in a cell-based assay, and was enzymatically inactive. This represents the first known mutation in the human LFNG gene and reinforces the hypothesis that proper regulation of the Notch signaling pathway is an absolute requirement for the correct patterning of the axial skeleton.

Analysis of protein sequence and interaction data for candidate disease gene prediction

Linkage analysis is a successful procedure to associate diseases with specific genomic regions. These regions are often large, containing hundreds of genes, which make experimental methods employed to identify the disease gene arduous and expensive. We present two methods to prioritize candidates for further experimental study: Common Pathway Scanning (CPS) and Common Module Profiling (CMP). CPS is based on the assumption that common phenotypes are associated with dysfunction in proteins that participate in the same complex or pathway. CPS applies network data derived from protein–protein interaction (PPI) and pathway databases to identify relationships between genes. CMP identifies likely candidates using a domain-dependent sequence similarity approach, based on the hypothesis that disruption of genes of similar function will lead to the same phenotype. Both algorithms use two forms of input data: known disease genes or multiple disease loci. When using known disease genes as input, our combined methods have a sensitivity of 0.52 and a specificity of 0.97 and reduce the candidate list by 13-fold. Using multiple loci, our methods successfully identify disease genes for all benchmark diseases with a sensitivity of 0.84 and a specificity of 0.63. Our combined approach prioritizes good candidates and will accelerate the disease gene discovery process.

Mutated MESP2 Causes Spondylocostal Dysostosis in Humans

Spondylocostal dysostosis (SCD) is a term given to a heterogeneous group of disorders characterized by abnormal vertebral segmentation (AVS). We have previously identified mutations in the Delta-like 3 (DLL3) gene as a major cause of autosomal recessive spondylocostal dysostosis. DLL3 encodes a ligand for the Notch receptor and, when mutated, defective somitogenesis occurs resulting in a consistent and distinctive pattern of AVS affecting the entire spine. From our study cohort of cases of AVS, we have identified individuals and families with abnormal segmentation of the entire spine but no mutations in DLL3, and, in some of these, linkage to the DLL3 locus at 19q13.1 has been excluded. Within this group, the radiological phenotype differs mildly from that of DLL3 mutation-positive SCD and is variable, suggesting further heterogeneity. Using a genomewide scanning strategy in one consanguineous family with two affected children, we demonstrated linkage to 15q21.3-15q26.1 and furthermore identified a 4-bp duplication mutation in the human MESP2 gene that codes for a basic helix-loop-helix transcription factor. No MESP2 mutations were found in a further 7 patients with related radiological phenotypes in whom abnormal segmentation affected all vertebrae, nor in a further 12 patients with diverse phenotypes.

Mechanism of allosteric regulation of transglutaminase 2 by GTP

Allosteric regulation is a fundamental mechanism of biological control. Here, we investigated the allosteric mechanism by which GTP inhibits cross-linking activity of transglutaminase 2 (TG2), a multifunctional protein, with postulated roles in receptor signaling, extracellular matrix assembly, and apoptosis. Our findings indicate that at least two components are involved in functionally coupling the allosteric site and active center of TG2, namely (i) GTP binding to mask a conformationally destabilizing switch residue, Arg-579, and to facilitate interdomain interactions that promote adoption of a compact, catalytically inactive conformation and (ii) stabilization of the inactive conformation by an uncommon H bond between a cysteine (Cys-277, an active center residue) and a tyrosine (Tyr-516, a residue located on a loop of the β-barrel 1 domain that harbors the GTP-binding site). Although not essential for GTP-mediated inhibition of cross-linking, this H bond enhances the rate of formation of the inactive conformer.

Evolution of distinct EGF domains with specific functions

EGF domains are extracellular protein modules cross‐linked by three intradomain disulfides. Past studies suggest the existence of two types of EGF domain with three‐disulfides, human EGF‐like (hEGF) domains and complement C1r‐like (cEGF) domains, but to date no functional information has been related to the two different types, and they are not differentiated in sequence or structure databases. We have developed new sequence patterns based on the different C‐termini to search specifically for the two types of EGF domains in sequence databases. The exhibited sensitivity and specificity of the new pattern‐based method represents a significant advancement over the currently available sequence detection techniques. We re‐annotated EGF sequences in the latest release of Swiss‐Prot looking for functional relationships that might correlate with EGF type. We show that important post‐translational modifications of three‐disulfide EGFs, including unusual forms of glycosylation and post‐translational proteolytic processing, are dependent on EGF subtype. For example, EGF domains that are shed from the cell surface and mediate intercellular signaling are all hEGFs, as are all human EGF receptor family ligands. Additional experimental data suggest that functional specialization has accompanied subtype divergence. Based on our structural analysis of EGF domains with three‐disulfide bonds and comparison to laminin and integrin‐like EGF domains with an additional inter‐domain disulfide, we propose that these hEGF and cEGF domains may have arisen from a four‐disulfide ancestor by selective loss of different cysteine residues.

Potential role of glutathione in evolution of thiol-based redox signaling sites in proteins

Cysteine is susceptible to a variety of modifications by reactive oxygen and nitrogen oxide species, including glutathionylation; and when two cysteines are involved, disulfide formation. Glutathione-cysteine adducts may be removed from proteins by glutaredoxin, whereas disulfides may be reduced by thioredoxin. Glutaredoxin is homologous to the disulfide-reducing thioredoxin and shares similar binding modes of the protein substrate. The evolution of these systems is not well characterized. When a single Cys is present in a protein, conjugation of the redox buffer glutathione may induce conformational changes, resulting in a simple redox switch that effects a signaling cascade. If a second cysteine is introduced into the sequence, the potential for disulfide formation exists. In favorable protein contexts, a bistable redox switch may be formed. Because of glutaredoxins similarities to thioredoxin, the mutated protein may be immediately exapted into the thioredoxin-dependent redox cycle upon addition of the second cysteine. Here we searched for examples of protein substrates where the number of redox-active cysteine residues has changed throughout evolution. We focused on cross-strand disulfides (CSDs), the most common type of forbidden disulfide. We searched for proteins where the CSD is present, absent and also found as a single cysteine in protein orthologs. Three different proteins were selected for detailed study—CD4, ERO1, and AKT. We created phylogenetic trees, examining when the CSD residues were mutated during protein evolution. We posit that the primordial cysteine is likely to be the cysteine of the CSD which undergoes nucleophilic attack by thioredoxin. Thus, a redox-active disulfide may be introduced into a protein structure by stepwise mutation of two residues in the native sequence to Cys. By extension, evolutionary acquisition of structural disulfides in proteins can potentially occur via transition through a redox-active disulfide state.

A despecialization step underlying evolution of a family of serine proteases.

In the trypsin superfamily of serine proteases, non-trypsin-like primary specificities have arisen in only two monophyletic descendent subbranches. We have recreated an ancestor to one of these subbranches (granzyme) using phylogenetic inference, gene synthesis, and protein expression. This ancestor has two unusual properties. First, it has broad primary specificity encompassing the entire repertoire of novel primary specificities found in its descendents. Second, unlike extant members that have narrow primary specificities, the ancestor exhibits tolerance to mutational changes in primary specificity-conferring residues-that is, structural plasticity. Molecular modeling and mutagenesis studies indicate that these unusual properties are due to a particularly wide substrate binding pocket. These two crucial properties of the ancestor not only distinguish it from its extant descendents but also from the trypsin-like proteases that preceded it. This indicates that a despecialization step, characterized by broad specificity and structural plasticity, underlies evolution of new primary specificities in this protease superfamily.

Epistatic Effects of Potassium Channel Variation on Cardiac Repolarization and Atrial Fibrillation Risk

OBJECTIVES The aim of this study was to evaluate the role of cardiac K(+) channel gene variants in families with atrial fibrillation (AF). BACKGROUND The K(+) channels play a major role in atrial repolarization but single mutations in cardiac K(+) channel genes are infrequently present in AF families. The collective effect of background K(+) channel variants of varying prevalence and effect size on the atrial substrate for AF is largely unexplored. METHODS Genes encoding the major cardiac K(+) channels were resequenced in 80 AF probands. Nonsynonymous coding sequence variants identified in AF probands were evaluated in 240 control subjects. Novel variants were characterized using patch-clamp techniques and in silico modeling was performed using the Courtemanche atrial cell model. RESULTS Nineteen nonsynonymous variants in 9 genes were found, including 11 rare variants. Rare variants were more frequent in AF probands (18.8% vs. 4.2%, p < 0.001), and the mean number of variants was greater (0.21 vs. 0.04, p < 0.001). The majority of K(+) channel variants individually had modest functional effects. Modeling simulations to evaluate combinations of K(+) channel variants of varying population frequency indicated that simultaneous small perturbations of multiple current densities had nonlinear interactions and could result in substantial (>30 ms) shortening or lengthening of action potential duration as well as increased dispersion of repolarization. CONCLUSIONS Families with AF show an excess of rare functional K(+) channel gene variants of varying phenotypic effect size that may contribute to an atrial arrhythmogenic substrate. Atrial cell modeling is a useful tool to assess epistatic interactions between multiple variants.

General mutation databases: analysis and review

Databases of mutations causing Mendelian disease play a crucial role in research, diagnostic and genetic health care and can play a role in life and death decisions. These databases are thus heavily used, but only gene or locus specific databases have been previously reviewed for completeness, accuracy, currency and utility. We have performed a review of the various general mutation databases that derive their data from the published literature and locus specific databases. Only two—the Human Gene Mutation Database (HGMD) and Online Mendelian Inheritance in Man (OMIM)—had useful numbers of mutations. Comparison of a number of characteristics of these databases indicated substantial inconsistencies between the two databases that included absent genes and missing mutations. This situation strengthens the case for gene specific curation of mutations and the need for an overall plan for collection, curation, storage and release of mutation data.

Evolutionary specialization of a tryptophan indole group for transition-state stabilization by eukaryotic transglutaminases

Covalent posttranslational protein modifications by eukaryotic transglutaminases proceed by a kinetic pathway of acylation and deacylation. Ammonia is released as the acylenzyme is formed, whereas the cross-linked product is released later in the deacylation step. Superposition of the active sites of transglutaminase type 2 (TG2) and the structurally related cysteine protease, papain, indicates that in the formation of tetrahedral intermediates, the backbone nitrogen of the catalytic Cys-277 and the Nε1 nitrogen of Trp-241 of TG2 could contribute to transition-state stabilization. The importance of this Trp-241 side chain was demonstrated by examining the kinetics of dansylcadaverine incorporation into a model peptide. Although substitution of the Trp-241 side chain with Ala or Gly had only a small effect on the Michaelis constant Km (1.5-fold increase), it caused a >300-fold lowering of the catalytic rate constant kcat. The wild-type and mutant TG2-catalyzed release of ammonia showed kinetics similar to the kinetics for the formation of cross-linked product, indicating that transition-state stabilization in the acylation step was rate-limiting. In papain, a Gln residue is at the position of TG2-Trp-241. The conservation of Trp-241 in all eukaryotic transglutaminases and the finding that W241Q-TG2 had a much lower kcat than wild-type enzyme suggest evolutionary specialization in the use of the indole group. This notion is further supported by the observation that transition-state-stabilizing side chains of Tyr and His that operate in some serine and metalloproteases only partially substituted for Trp.