Osteopathia striata with cranial sclerosis as a cancer predisposition syndrome: The first report of neuroblastoma and review of all cancers in OSCS.

Osteopathia Striata with Cranial Sclerosis (OSCS) is a rare genetic condition primarily characterized by metaphyseal striations of long bones, bone sclerosis, macrocephaly, and other congenital anomalies. It is caused by pathogenic variants in AMER1, a tumor suppressor and a WNT signaling repressor gene with key roles in tissue regeneration, neurodevelopment, tumorigenesis, and other developmental processes. While somatic AMER1 pathogenic variants have frequently been identified in several tumor types (e.g., Wilms tumor and colorectal cancer), whether OSCS (i.e., with AMER1 germline variants) is a tumor predisposition syndrome is not clear, with only nine cases reported with tumors. We here report the first case of neuroblastoma diagnosed in a male child with OSCS, review all previously reported tumors diagnosed in individuals with OSCS, and discuss potential tumorigenic mechanisms of AMER1. Our report adds to the accumulating evidence suggesting OSCS is a tumor predisposition condition, highlighting the importance of maintaining a high index of suspicion for the associated tumors when evaluating patients with OSCS. Importantly, Wilms tumor stands out as the most commonly observed tumor in OSCS patients, underscoring the need for regular surveillance.

A novel FAM20C mutation causes a rare form of neonatal lethal Raine syndrome.

Raine syndrome is a rare, autosomal recessive, osteosclerotic bone dysplasia due to pathogenic variants in FAM20C. The clinical phenotype is characterized by generalized osteosclerosis affecting all bones, cerebral calcifications, and craniofacial dysmorphism. Most cases present during the neonatal period with early lethality due to pulmonary hypoplasia and respiratory compromise while only few affected individuals have been reported to survive into adulthood. FAM20C is a ubiquitously expressed protein kinase that contains five functional domains including a catalytic domain, a binding pocket for FAM20A and three distinct N-glycosylation sites. We report a newborn infant with a history of prenatal onset fractures, generalized osteosclerosis, and craniofacial dysmorphism and early lethality. The clinical presentation was highly suggestive of Raine syndrome. A homozygous, novel missense variant in exon 5 of FAM20C (c.1007T>G; p.Met336Arg) was identified by targeted Sanger sequencing. Following in silico analysis and mapping of the variant on a three-dimensional (3D) model of FAM20C it is predicted to be deleterious and to affect N-glycosylation, protein folding, and subsequent secretion of FAM20C. In addition, we reviewed all published FAM20C mutations and observed that most pathogenic variants affect functional regions within the protein establishing evidence for an emerging genotype-phenotype correlation.

Spectrum of neurodevelopmental disease associated with the GNAO1 guanosine triphosphate-binding region.

OBJECTIVE: To characterize the phenotypic spectrum associated with GNAO1 variants and establish genotype-protein structure-phenotype relationships. METHODS: We evaluated the phenotypes of 14 patients with GNAO1 variants, analyzed their variants for potential pathogenicity, and mapped them, along with those in the literature, on a three-dimensional structural protein model. RESULTS: The 14 patients in our cohort, including one sibling pair, had 13 distinct, heterozygous GNAO1 variants classified as pathogenic or likely pathogenic. We attributed the same variant in two siblings to parental mosaicism. Patients initially presented with seizures beginning in the first 3 months of life (8/14), developmental delay (4/14), hypotonia (1/14), or movement disorder (1/14). All patients had hypotonia and developmental delay ranging from mild to severe. Nine had epilepsy, and nine had movement disorders, including dystonia, ataxia, chorea, and dyskinesia. The 13 GNAO1 variants in our patients are predicted to result in amino acid substitutions or deletions in the GNAO1 guanosine triphosphate (GTP)-binding region, analogous to those in previous publications. Patients with variants affecting amino acids 207-221 had only movement disorder and hypotonia. Patients with variants affecting the C-terminal region had the mildest phenotypes. SIGNIFICANCE: GNAO1 encephalopathy most frequently presents with seizures beginning in the first 3 months of life. Concurrent movement disorders are also a prominent feature in the spectrum of GNAO1 encephalopathy. All variants affected the GTP-binding domain of GNAO1, highlighting the importance of this region for G-protein signaling and neurodevelopment.

ExoLocator--an online view into genetic makeup of vertebrate proteins.

ExoLocator (http://exolocator.eopsf.org) collects in a single place information needed for comparative analysis of protein-coding exons from vertebrate species. The main source of data--the genomic sequences, and the existing exon and homology annotation--is the ENSEMBL database of completed vertebrate genomes. To these, ExoLocator adds the search for ostensibly missing exons in orthologous protein pairs across species, using an extensive computational pipeline to narrow down the search region for the candidate exons and find a suitable template in the other species, as well as state-of-the-art implementations of pairwise alignment algorithms. The resulting complements of exons are organized in a way currently unique to ExoLocator: multiple sequence alignments, both on the nucleotide and on the peptide levels, clearly indicating the exon boundaries. The alignments can be inspected in the web-embedded viewer, downloaded or used on the spot to produce an estimate of conservation within orthologous sets, or functional divergence across paralogues.

Active site coupling in PDE:PKA complexes promotes resetting of mammalian cAMP signaling.

Cyclic 3'5' adenosine monophosphate (cAMP)-dependent-protein kinase (PKA) signaling is a fundamental regulatory pathway for mediating cellular responses to hormonal stimuli. The pathway is activated by high-affinity association of cAMP with the regulatory subunit of PKA and signal termination is achieved upon cAMP dissociation from PKA. Although steps in the activation phase are well understood, little is known on how signal termination/resetting occurs. Due to the high affinity of cAMP to PKA (KD ∼ low nM), bound cAMP does not readily dissociate from PKA, thus begging the question of how tightly bound cAMP is released from PKA to reset its signaling state to respond to subsequent stimuli. It has been recently shown that phosphodiesterases (PDEs) can catalyze dissociation of bound cAMP and thereby play an active role in cAMP signal desensitization/termination. This is achieved through direct interactions with the regulatory subunit of PKA, thereby facilitating cAMP dissociation and hydrolysis. In this study, we have mapped direct interactions between a specific cyclic nucleotide phosphodiesterase (PDE8A) and a PKA regulatory subunit (RIα isoform) in mammalian cAMP signaling, by a combination of amide hydrogen/deuterium exchange mass spectrometry, peptide array, and computational docking. The interaction interface of the PDE8A:RIα complex, probed by peptide array and hydrogen/deuterium exchange mass spectrometry, brings together regions spanning the phosphodiesterase active site and cAMP-binding sites of RIα. Computational docking combined with amide hydrogen/deuterium exchange mass spectrometry provided a model for parallel dissociation of bound cAMP from the two tandem cAMP-binding domains of RIα. Active site coupling suggests a role for substrate channeling in the PDE-dependent dissociation and hydrolysis of cAMP bound to PKA. This is the first instance, to our knowledge, of PDEs directly interacting with a cAMP-receptor protein in a mammalian system, and highlights an entirely new class of binding partners for RIα. This study also highlights applications of structural mass spectrometry combined with computational docking for mapping dynamics in transient signaling protein complexes. Together, these results present a novel and critical role for phosphodiesterases in moderating local concentrations of cAMP in microdomains and signal resetting.

Cube-DB: detection of functional divergence in human protein families.

Cube-DB is a database of pre-evaluated results for detection of functional divergence in human/vertebrate protein families. The analysis is organized around the nomenclature associated with the human proteins, but based on all currently available vertebrate genomes. Using full genomes enables us, through a mutual-best-hit strategy, to construct comparable taxonomical samples for all paralogues under consideration. Functional specialization is scored on the residue level according to two models of behavior after divergence: heterotachy and homotachy. In the first case, the positions on the protein sequence are scored highly if they are conserved in the reference group of orthologs, and overlap poorly with the residue type choice in the paralogs groups (such positions will also be termed functional determinants). The second model additionally requires conservation within each group of paralogs (functional discriminants). The scoring functions are phylogeny independent, but sensitive to the residue type similarity. The results are presented as a table of per-residue scores, and mapped onto related structure (when available) via browser-embedded visualization tool. They can also be downloaded as a spreadsheet table, and sessions for two additional molecular visualization tools. The database interface is available at http://epsf.bmad.bii.a-star.edu.sg/cube/db/html/home.html.

WLS-dependent secretion of WNT3A requires Ser209 acylation and vacuolar acidification.

Wnt proteins are secreted post-translationally modified proteins that signal locally to regulate development and proliferation. The production of bioactive Wnts requires a number of dedicated factors in the secreting cell whose coordinated functions are not fully understood. A screen for small molecules identified inhibitors of vacuolar acidification as potent inhibitors of Wnt secretion. Inhibition of the V-ATPase or disruption of vacuolar pH gradients by diverse drugs potently inhibited Wnt/ß-catenin signaling both in cultured human cells and in vivo, and impaired Wnt-regulated convergent extension movements in Xenopus embryos. WNT secretion requires its binding to the carrier protein wntless (WLS); we find that WLS is ER-resident in human cells and WNT3A binding to WLS requires PORCN-dependent lipid modification of WNT3A at serine 209. Inhibition of vacuolar acidification results in accumulation of the WNT3A-WLS complex both in cells and at the plasma membrane. Modeling predictions suggest that WLS has a lipid-binding ß-barrel that is similar to the lipocalin-family fold. We propose that WLS binds Wnts in part through a lipid-binding domain, and that vacuolar acidification is required to release palmitoylated WNT3A from WLS in secretory vesicles, possibly to facilitate transfer of WNT3A to a soluble carrier protein.

Distinct faces of the Ku heterodimer mediate DNA repair and telomeric functions.

The Ku heterodimer, comprised of Ku70 and Ku80 subunits, is a conserved complex involved in nonhomologous end-joining (NHEJ). However, it also functions in maintenance of telomeres, chromosome termini normally resistant to end-joining events. To elucidate the spatial organization of these functions, we rationally guided Ku mutagenesis in yeast with real-valued evolutionary trace (rvET). This revealed two ancestrally related alpha-helices: one on the Ku70 surface that is required in yeast for NHEJ, and a second on the Ku80 surface that is required in yeast for telomeric heterochromatin formation. When bound to a DNA end, the surface containing the NHEJ-specific Ku70 helix is oriented toward the DNA terminus, whereas the surface containing the telomeric function-specific Ku80 helix faces inward, toward telomeric chromatin, when bound to a telomere. We propose a 'two-face' model for Ku and that divergent evolution of these faces allowed Ku's dual role in NHEJ and telomere maintenance.

deconSTRUCT: general purpose protein database search on the substructure level.

deconSTRUCT webserver offers an interface to a protein database search engine, usable for a general purpose detection of similar protein (sub)structures. Initially, it deconstructs the query structure into its secondary structure elements (SSEs) and reassembles the match to the target by requiring a (tunable) degree of similarity in the direction and sequential order of SSEs. Hierarchical organization and judicious use of the information about protein structure enables deconSTRUCT to achieve the sensitivity and specificity of the established search engines at orders of magnitude increased speed, without tying up irretrievably the substructure information in the form of a hash. In a post-processing step, a match on the level of the backbone atoms is constructed. The results presented to the user consist of the list of the matched SSEs, the transformation matrix for rigid superposition of the structures and several ways of visualization, both downloadable and implemented as a web-browser plug-in. The server is available at http://epsf.bmad.bii.a-star.edu.sg/struct_server.html.

Reduced representation of protein structure: implications on efficiency and scope of detection of structural similarity.

BACKGROUND: Computational comparison of two protein structures is the starting point of many methods that build on existing knowledge, such as structure modeling (including modeling of protein complexes and conformational changes), molecular replacement, or annotation by structural similarity. In a commonly used strategy, significant effort is invested in matching two sets of atoms. In a complementary approach, a global descriptor is assigned to the overall structure, thus losing track of the substructures within. RESULTS: Using a small set of geometric features, we define a reduced representation of protein structure, together with an optimizing function for matching two representations, to provide a pre-filtering stage in a database search. We show that, in a straightforward implementation, the representation performs well in terms of resolution in the space of protein structures, and its ability to make new predictions. CONCLUSIONS: Perhaps unexpectedly, a substantial discriminating power already exists at the level of main features of protein structure, such as directions of secondary structural elements, possibly constrained by their sequential order. This can be used toward efficient comparison of protein (sub)structures, allowing for various degrees of conformational flexibility within the compared pair, which in turn can be used for modeling by homology of protein structure and dynamics.

pFlexAna: detecting conformational changes in remotely related proteins.

The pFlexAna (protein flexibility analyzer) web server detects and displays conformational changes in remotely related proteins, without relying on sequence homology. To do so, it first applies a reliable statistical test to align core protein fragments that are structurally similar and then clusters these aligned fragment pairs into 'super-alignments', according to the similarity of geometric transformations that align them. The result is that the dominant conformational changes occur between the clusters, while the smaller conformational changes occur within a cluster. pFlexAna is available at http://bigbird.comp.nus.edu.sg/pfa2/.

Cancer-associated mutations in the ribosomal protein L5 gene dysregulate the HDM2/p53-mediated ribosome biogenesis checkpoint.

Perturbations in ribosome biogenesis have been associated with cancer. Such aberrations activate p53 through the RPL5/RPL11/5S rRNA complex-mediated inhibition of HDM2. Studies using animal models have suggested that this signaling pathway might constitute an important anticancer barrier. To gain a deeper insight into this issue in humans, here we analyze somatic mutations in RPL5 and RPL11 coding regions, reported in The Cancer Genome Atlas and International Cancer Genome Consortium databases. Using a combined computational and statistical approach, complemented by a range of biochemical and functional analyses in human cancer cell models, we demonstrate the existence of several mechanisms by which RPL5 mutations may impair wild-type p53 upregulation and ribosome biogenesis. Unexpectedly, the same approach provides only modest evidence for a similar role of RPL11, suggesting that RPL5 represents a preferred target during human tumorigenesis in cancers with wild-type p53. Furthermore, we find that several functional cancer-associated RPL5 somatic mutations occur as rare germline variants in general population. Our results shed light on the so-far enigmatic role of cancer-associated mutations in genes encoding ribosomal proteins, with implications for our understanding of the tumor suppressive role of the RPL5/RPL11/5S rRNA complex in human malignancies.

The C-Terminal Extension Unique to the Long Isoform of the Shelterin Component TIN2 Enhances Its Interaction with TRF2 in a Phosphorylation- and Dyskeratosis Congenita Cluster-Dependent Fashion.

TIN2 is central to the shelterin complex, linking the telomeric proteins TRF1 and TRF2 with TPP1/POT1. Mutations in TINF2, which encodes TIN2, that are found in dyskeratosis congenita (DC) result in very short telomeres and cluster in a region shared by the two TIN2 isoforms, TIN2S (short) and TIN2L (long). Here we show that TIN2L, but not TIN2S, is phosphorylated. TRF2 interacts more with TIN2L than TIN2S, and both the DC cluster and phosphorylation promote this enhanced interaction. The binding of TIN2L, but not TIN2S, is affected by TRF2-F120, which is also required for TRF2's interaction with end processing factors such as Apollo. Conversely, TRF1 interacts more with TIN2S than with TIN2L. A DC-associated mutation further reduces TIN2L-TRF1, but not TIN2S-TRF1, interaction. Cells overexpressing TIN2L or phosphomimetic TIN2L are permissive to telomere elongation, whereas cells overexpressing TIN2S or phosphodead TIN2L are not. Telomere lengths are unchanged in cell lines in which TIN2L expression has been eliminated by clustered regularly interspaced short palindromic repeat (CRISPR)/Cas9-mediated mutation. These results indicate that TIN2 isoforms are biochemically and functionally distinguishable and that shelterin composition could be fundamentally altered in patients with TINF2 mutations.

Utility of rapid whole-exome sequencing in the diagnosis of Niemann-Pick disease type C presenting with fetal hydrops and acute liver failure.

Rapid whole-exome sequencing (rWES) is used in critically ill newborn infants to inform about diagnosis, clinical management, and prognosis. Here we report a male newborn infant with hydrops, pancytopenia, and acute liver failure who was listed for liver transplantation. Given the acuity of the presentation, the procedure-related morbidity and mortality, and lack of diagnosis, we used rWES in the proband and both parents with a turnaround time of 10 business days. rWES returned one maternally inherited, likely pathogenic and one paternally inherited, likely pathogenic variant in NPC1, suggestive of a diagnosis of Niemann-Pick disease type C (NPC). Interestingly, a diagnosis of NPC was entertained prior to rWES, but deemed unlikely in light of absent cholesterol storage on liver biopsy and near-normal oxysterol levels in dried blood. The diagnosis of NPC was confirmed on filipin stain in fibroblasts demonstrating defective cholesterol trafficking. NPC is a slowly progressive neurodegenerative disorder that may also affect the liver with overall poor prognosis. It was decided to take the infant off the transplant list and transfer to palliative care, where he died after 4 wk. This case highlights the utility of rWES in an acute clinical setting for several domains of precision medicine including (1) diagnosis, (2) prognosis and outcome, (3) management and therapy, and (4) utilization of resources.

A defect in the inner kinetochore protein CENPT causes a new syndrome of severe growth failure.

Primordial growth failure has been linked to defects in the biology of cell division and replication. The complex processes involved in microtubule spindle formation, organization and function have emerged as a dominant patho-mechanism in these conditions. The majority of reported disease genes encode for centrosome and centriole proteins, leaving kinetochore proteins by which the spindle apparatus interacts with the chromosomes largely unaccounted for. We report a novel disease gene encoding the constitutive inner kinetochore member CENPT, which is involved in kinetochore targeting and assembly, resulting in severe growth failure in two siblings of a consanguineous family. We herein present studies on the molecular and cellular mechanisms that explain how genetic mutations in this gene lead to primordial growth failure. In both, affected human cell lines and a zebrafish knock-down model of Cenpt, we observed aberrations in cell division with abnormal accumulation of micronuclei and of nuclei with increased DNA content arising from incomplete and/or irregular chromosomal segregation. Our studies underscore the critical importance of kinetochore function for overall body growth and provide new insight into the cellular mechanisms implicated in the spectrum of these severe growth disorders.

Rank information: a structure-independent measure of evolutionary trace quality that improves identification of protein functional sites.

Protein functional sites are key targets for drug design and protein engineering, but their large-scale experimental characterization remains difficult. The evolutionary trace (ET) is a computational approach to this problem that has been useful in a variety of case studies, but its proteomic scale application is partially hindered because automated retrieval of input sequences from databases often includes some with errors that degrade functional site identification. To recognize and purge these sequences, this study introduces a novel and structure-free measure of ET quality called rank information (RI). It is shown that RI decreases in response to errors in sequences, alignments, or functional classifications. Conversely, an automated procedure to increase RI by selectively removing sequences improves functional site identification so as to nearly match manually curated traces in kinases and in a test set of 79 diverse proteins. Thus we conclude that RI partially reflects the evolutionary consistency of sequence, structure, and function. In practice, as the size of the proteome continues to grow exponentially, it provides a novel and structure-free measure of ET quality that increases its accuracy for large-scale automated annotation of protein functional sites.

Role of transmembrane domain/transmembrane domain interfaces of P-glycoprotein (ABCB1) in solute transport. Convergent information from photoaffinity labeling, site directed mutagenesis and in silico importance prediction.

Human P-glycoprotein (P-gp, ABCB1) plays an important role in the development of resistance to anticancer therapy. This ABC-transporter (ATP-binding cassette transporter) intercepts drugs at the level of the plasma membrane and effluxes them before they are able to reach their intracellular target structures. Inhibition of P-gp by low molecular weight compounds has been advocated as a concept for resensitization of cells to anticancer agents and several clinical studies in oncological patients have advanced to phase III. Even more importantly, P-glycoprotein also represents an antitarget. Its expression in cells lining the intestinal tract, the canalicular side of hepatocytes, renal tubuli and the blood brain barrier lead to interference with pharmacokinetics of compounds that are recognized as pump substrates. An early prediction of ADMET (Absorption-Distribution-Metabolism-Excretion-Toxicity) properties is important during drug development, since interference of a compound with P-gp might compromise its future development into a drug. Despite considerable efforts, the mechanism by which P-gp binds and transports its solutes remains unclear. Generation of homology models of the protein allowed integration of data obtained by photoaffinity labeling, in silico prediction of functional importance by evolutionary tracing and site directed mutagenesis. An integral view of data indicates that these three lines of evidence converge to indicate two pseudosymmetric P-gp drug binding pockets located at the two transmembrane domain interfaces.

An accurate, sensitive, and scalable method to identify functional sites in protein structures.

Functional sites determine the activity and interactions of proteins and as such constitute the targets of most drugs. However, the exponential growth of sequence and structure data far exceeds the ability of experimental techniques to identify their locations and key amino acids. To fill this gap we developed a computational Evolutionary Trace method that ranks the evolutionary importance of amino acids in protein sequences. Studies show that the best-ranked residues form fewer and larger structural clusters than expected by chance and overlap with functional sites, but until now the significance of this overlap has remained qualitative. Here, we use 86 diverse protein structures, including 20 determined by the structural genomics initiative, to show that this overlap is a recurrent and statistically significant feature. An automated ET correctly identifies seven of ten functional sites by the least favorable statistical measure, and nine of ten by the most favorable one. These results quantitatively demonstrate that a large fraction of functional sites in the proteome may be accurately identified from sequence and structure. This should help focus structure-function studies, rational drug design, protein engineering, and functional annotation to the relevant regions of a protein.

Cube - an online tool for comparison and contrasting of protein sequences.

When comparing sequences of similar proteins, two kinds of questions can be asked, and the related two kinds of inference made. First, one may ask to what degree they are similar, and then, how they differ. In the first case one may tentatively conclude that the conserved elements common to all sequences are of central and common importance to the protein's function. In the latter case the regions of specialization may be discriminative of the function or binding partners across subfamilies of related proteins. Experimental efforts - mutagenesis or pharmacological intervention - can then be pointed in either direction, depending on the context of the study. Cube simplifies this process for users that already have their favorite sets of sequences, and helps them collate the information by visualization of the conservation and specialization scores on the sequence and on the structure, and by spreadsheet tabulation. All information can be visualized on the spot, or downloaded for reference and later inspection. SERVER HOMEPAGE: http://eopsf.org/cube.

TRF2 interaction with Ku heterotetramerization interface gives insight into c-NHEJ prevention at human telomeres.

Telomeres are protected from nonhomologous end-joining (NHEJ) to avoid deleterious chromosome fusions, yet they associate with the Ku heterodimer that is principal in the classical NHEJ (c-NHEJ) pathway. T-loops have been proposed to inhibit Ku's association with telomeric ends, thus inhibiting c-NHEJ; however, deficiencies in the t-loop model suggest additional mechanisms are in effect. We demonstrate that TRF2 interacts with Ku at telomeres and via residues in Ku70 helix 5 (α5), which are vital for NHEJ. We show that Ku's interaction with a TRF2 mutant that induces telomeric fusions is significantly impaired. Additionally, we demonstrate that Ku70 α5 is required for Ku self-association in live cells, which can bridge DNA ends. Together, these findings lead us to propose a model in which telomeres are directly protected from c-NHEJ via TRF2 impeding Ku's ability to synapse telomere ends.