De novo coding variants in the AGO1 gene cause a neurodevelopmental disorder with intellectual disability

BACKGROUND

High-impact pathogenic variants in more than a thousand genes are involved in Mendelian forms of neurodevelopmental disorders (NDD).

METHODS

This study describes the molecular and clinical characterisation of 28 probands with NDD harbouring heterozygous AGO1 coding variants, occurring de novo for all those whose transmission could have been verified (26/28).

RESULTS

A total of 15 unique variants leading to amino acid changes or deletions were identified: 12 missense variants, two in-frame deletions of one codon, and one canonical splice variant leading to a deletion of two amino acid residues. Recurrently identified variants were present in several unrelated individuals: p.(Phe180del), p.(Leu190Pro), p.(Leu190Arg), p.(Gly199Ser), p.(Val254Ile) and p.(Glu376del). AGO1 encodes the Argonaute 1 protein, which functions in gene-silencing pathways mediated by small non-coding RNAs. Three-dimensional protein structure predictions suggest that these variants might alter the flexibility of the AGO1 linker domains, which likely would impair its function in mRNA processing. Affected individuals present with intellectual disability of varying severity, as well as speech and motor delay, autistic behaviour and additional behavioural manifestations.

CONCLUSION

Our study establishes that de novo coding variants in AGO1 are involved in a novel monogenic form of NDD, highly similar to the recently reported AGO2-related NDD.

Enhanced interpretation of 935 hotspot and non-hotspot RAS variants using evidence-based structural bioinformatics

In the current study, we report computational scores for advancing genomic interpretation of disease-associated genomic variation in members of the RAS family of genes. For this purpose, we applied 31 sequence- and 3D structure-based computational scores, chosen by their breadth of biophysical properties. We parametrized our data by assembling a numerically homogenized experimentally-derived dataset, which when use in our calculations reveal that computational scores using 3D structure highly correlate with experimental measures (e.g., GAP-mediated hydrolysis R_Spearman = 0.80 and RAF affinity R_spearman = 0.82), while sequence-based scores are discordant with this data. Performing all-against-all comparisons, we applied this parametrized modeling approach to the study of 935 RAS variants from 7 RAS genes, which led us to identify 4 groups of mutations according to distinct biochemical scores within each group. Each group was comprised of hotspot and non-hotspot KRAS variants, indicating that poorly characterized variants could functionally behave like pathogenic mutations. Combining computational scores using dimensionality reduction indicated that changes to local unfolding propensity associate with changes in enzyme activity by genomic variants. Hence, our systematic approach, combining methodologies from both clinical genomics and 3D structural bioinformatics, represents an expansion for interpreting genomic data, provides information of mechanistic value, and that is transferable to other proteins.

Computational modeling reveals key molecular properties and dynamic behavior of disruptor of telomeric silencing 1-like (DOT1L) and partnering complexes involved in leukemogenesis

Disruptor of telomeric silencing 1‐like (DOT1L) is the only non‐SET domain histone lysine methyltransferase (KMT) and writer of H3K79 methylation on nucleosomes marked by H2B ubiquitination. DOT1L has elicited significant attention because of its interaction or fusion with members of the AF protein family in blood cell biology and leukemogenic transformation. Here, our goal was to extend previous structural information by performing a robust molecular dynamic study of DOT1L and its leukemogenic partners combined with mutational analysis. We show that statically and dynamically, D161, G163, E186, and F223 make frequent time‐dependent interactions with SAM, while additional residues T139, K187, and N241 interact with SAM only under dynamics. Dynamics models reveal DOT1L, SAM, and H4 moving as one and show that more than twice the number of DOT1L residues interacts with these partners, relative to the static structure. Mutational analyses indicate that six of these residues are intolerant to substitution. We describe the dynamic behavior of DOT1L interacting with AF10 and AF9. Studies on the dynamics of a heterotrimeric complex of DOT1L1‐AF10 illuminated describe coordinated motions that impact the relative position of the DOT1L HMT domain to the nucleosome. The molecular motions of the DOT1L–AF9 complex are less extensive and highly dynamic, resembling a swivel‐like mechanics. Through molecular dynamics and mutational analysis, we extend the knowledge previous provided by static measurements. These results are important to consider when describing the biochemical properties of DOT1L, under normal and in disease conditions, as well as for the development of novel therapeutic agents.

Structural bioinformatics enhances mechanistic interpretation of genomic variation, demonstrated through the analyses of 935 distinct RAS family mutations

MOTIVATION

Protein-coding genetic alterations are frequently observed in Clinical Genetics, but the high yield of variants of uncertain significance remains a limitation in decision making. RAS-family GTPases are cancer drivers, but only 54 variants, across all family members, fall within well-known hotspots. However, extensive sequencing has identified 881 non-hotspot variants for which significance remains to be investigated.

RESULTS

Here, we evaluate 935 missense variants from seven RAS genes, observed in cancer, RASopathies and the healthy adult population. We characterized hotspot variants, previously studied experimentally, using 63 sequence- and 3D structure-based scores, chosen by their breadth of biophysical properties. Applying scores that display best correlation with experimental measures, we report new valuable mechanistic inferences for both hot-spot and non-hotspot variants. Moreover, we demonstrate that 3D scores have little-to-no correlation with those based on DNA sequence, which are commonly used in Clinical Genetics. Thus, combined, these new knowledge bear significant relevance.

AVAILABILITY AND IMPLEMENTATION

All genomic and 3D scores, and markdown for generating figures, are provided in our supplemental data.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Molecular mechanics and dynamic simulations of well-known Kabuki syndrome-associated KDM6A variants reveal putative mechanisms of dysfunction

BACKGROUND

Kabuki syndrome is a genetic disorder that affects several body systems and presents with variations in symptoms and severity. The syndrome is named for a common phenotype of faces resembling stage makeup used in a Japanese traditional theatrical art named kabuki. The most frequent cause of this syndrome is mutations in the H3K4 family of histone methyltransferases while a smaller percentage results from genetic alterations affecting the histone demethylase, KDM6A. Because of the rare presentation of the latter form of the disease, little is known about how missense changes in the KDM6A protein sequence impact protein function.

RESULTS

In this study, we use molecular mechanic and molecular dynamic simulations to enhance the annotation and mechanistic interpretation of the potential impact of eleven KDM6A missense variants found in Kabuki syndrome patients. These variants (N910S, D980V, S1025G, C1153R, C1153Y, P1195L, L1200F, Q1212R, Q1248R, R1255W, and R1351Q) are predicted to be pathogenic, likely pathogenic or of uncertain significance by sequence-based analysis. Here, we demonstrate, for the first time, that although Kabuki syndrome missense variants are found outside the functionally critical regions, they could affect overall function by significantly disrupting global and local conformation (C1153R, C1153Y, P1195L, L1200F, Q1212R, Q1248R, R1255W and R1351Q), chemical environment (C1153R, C1153Y, P1195L, L1200F, Q1212R, Q1248R, R1255W and R1351Q), and/or molecular dynamics of the catalytic domain (all variants). In addition, our approaches predict that many mutations, in particular C1153R, could allosterically disrupt the key enzymatic interactions of KDM6A.

CONCLUSIONS

Our study demonstrates that the KDM6A Kabuki syndrome variants may impair histone demethylase function through various mechanisms that include altered protein integrity, local environment, molecular interactions and protein dynamics. Molecular dynamics simulations of the wild type and the variants are critical to gain a better understanding of molecular dysfunction. This type of comprehensive structure- and MD-based analyses should help develop improved impact scoring systems to interpret the damaging effects of variants in this protein and other related proteins as well as provide detailed mechanistic insight that is not currently predictable from sequence alone.

Defining the genotypic and phenotypic spectrum of X-linked MSL3-related disorder

Purpose

We sought to delineate the genotypic and phenotypic spectrum of female and male individuals with X-linked, MSL3-related disorder (Basilicata–Akhtar syndrome).

Methods

Twenty-five individuals (15 males, 10 females) with causative variants in MSL3 were ascertained through exome or genome sequencing at ten different sequencing centers.

Results

We identified multiple variant types in MSL3 (ten nonsense, six frameshift, four splice site, three missense, one in-frame-deletion, one multi-exon deletion), most proven to be de novo, and clustering in the terminal eight exons suggesting that truncating variants in the first five exons might be compensated by an alternative MSL3 transcript. Three-dimensional modeling of missense and splice variants indicated that these have a deleterious effect. The main clinical findings comprised developmental delay and intellectual disability ranging from mild to severe. Autism spectrum disorder, muscle tone abnormalities, and macrocephaly were common as well as hearing impairment and gastrointestinal problems. Hypoplasia of the cerebellar vermis emerged as a consistent magnetic resonance image (MRI) finding. Females and males were equally affected. Using facial analysis technology, a recognizable facial gestalt was determined.

Conclusion

Our aggregated data illustrate the genotypic and phenotypic spectrum of X-linked, MSL3-related disorder (Basilicata–Akhtar syndrome). Our cohort improves the understanding of disease related morbidity and allows us to propose detailed surveillance guidelines for affected individuals.

Interaction between Mas1 and AT1RA contributes to enhancement of skeletal muscle angiogenesis by angiotensin-(1-7) in Dahl salt-sensitive rats

The heptapeptide angiotensin-(1-7) (Ang-(1-7)) is protective in the cardiovascular system through its induction of vasodilator production and angiogenesis. Despite acting antagonistically to the effects of elevated, pathophysiological levels of angiotensin II (AngII), recent evidence has identified convergent and beneficial effects of low levels of both Ang-(1-7) and AngII. Previous work identified the AngII receptor type I (AT1R) as a component of the protein complex formed when Ang-(1-7) binds its receptor, Mas1. Importantly, pharmacological blockade of AT1R did not alter the effects of Ang-(1-7). Here, we use a novel mutation of AT1RA in the Dahl salt-sensitive (SS) rat to test the hypothesis that interaction between Mas1 and AT1R contributes to proangiogenic Ang-(1-7) signaling. In a model of hind limb angiogenesis induced by electrical stimulation, we find that the restoration of skeletal muscle angiogenesis in SS rats by Ang-(1-7) infusion is impaired in AT1RA knockout rats. Enhancement of endothelial cell (EC) tube formation capacity by Ang-(1-7) is similarly blunted in AT1RA mutant ECs. Transcriptional changes elicited by Ang-(1-7) in SS rat ECs are altered in AT1RA mutant ECs, and tandem mass spectrometry-based proteomics demonstrate that the protein complex formed upon binding of Ang-(1-7) to Mas1 is altered in AT1RA mutant ECs. Together, these data support the hypothesis that interaction between AT1R and Mas1 contributes to proangiogenic Ang-(1-7) signaling.

Novel KLHL26 variant associated with a familial case of Ebstein's anomaly and left ventricular noncompaction

BACKGROUND

Ebstein's anomaly (EA) is a rare congenital heart disease of the tricuspid valve and right ventricle. Patients with EA often manifest with left ventricular noncompaction (LVNC), a cardiomyopathy. Despite implication of cardiac sarcomere genes in some cases, very little is understood regarding the genetic etiology of EA/LVNC. Our study describes a multigenerational family with at least 10 of 17 members affected by EA/LVNC.

METHODS

We performed echocardiography on all family members and conducted exome sequencing of six individuals. After identifying candidate variants using two different bioinformatic strategies, we confirmed segregation with phenotype using Sanger sequencing. We investigated structural implications of candidate variants using protein prediction models.

RESULTS

Exome sequencing analysis of four affected and two unaffected members identified a novel, rare, and damaging coding variant in the Kelch-like family member 26 (KLHL26) gene located on chromosome 19 at position 237 of the protein (GRCh37). This variant region was confirmed by Sanger sequencing in the remaining family members. KLHL26 (c.709C > T p.R237C) segregates only with EA/LVNC-affected individuals (FBAT p < .05). Investigating structural implications of the candidate variant using protein prediction models suggested that the KLHL26 variant disrupts electrostatic interactions when binding to part of the ubiquitin proteasome, specifically Cullin3 (CUL3), a component of E3 ubiquitin ligase.

CONCLUSION

In this familial case of EA/LVNC, we have identified a candidate gene variant, KLHL26 (p.R237C), which may have an important role in ubiquitin-mediated protein degradation during cardiac development.

Discovery, expression, cellular localization, and molecular properties of a novel, alternative spliced HP1γ isoform, lacking the chromoshadow domain

By reading the H3K9Me3 mark through their N-terminal chromodomain (CD), HP1 proteins play a significant role in cancer-associated processes, including cell proliferation, differentiation, chromosomal stability, and DNA repair. Here, we used a combination of bioinformatics-based methodologies, as well as experimentally-derived datasets, that reveal the existence of a novel short HP1γ (CBX3) isoform, named here sHP1γ, generated by alternative splicing of the CBX3 locus. The sHP1γ mRNA encodes a protein composed of 101 residues and lacks the C-terminal chromoshadow domain (CSD) that is required for dimerization and heterodimerization in the previously described 183 a. a HP1γ protein. Fold recognition, order-to-disorder calculations, threading, homology-based molecular modeling, docking, and molecular dynamic simulations show that the sHP1γ is comprised of a CD flanked by intrinsically disordered regions (IDRs) with an IDR-CD-IDR domain organization and likely retains the ability to bind to the H3K9Me3. Both qPCR analyses and mRNA-seq data derived from large-scale studies confirmed that sHP1γ mRNA is expressed in the majority of human tissues at approximately constant ratios with the chromoshadow domain containing isoform. However, sHP1γ mRNA levels appear to be dysregulated in different cancer types. Thus, our data supports the notion that, due to the existence of functionally different isoforms, the regulation of HP1γ-mediated functions is more complex than previously anticipated.

Novel biallelic variants in MSTO1 associated with mitochondrial myopathy

Mitochondrial disorders are caused by nuclear and mitochondrial pathogenic variants leading to defects in mitochondrial function and cellular respiration. Recently, the nuclear-encoded mitochondrial fusion gene MSTO1 (Misato 1) has been implicated in mitochondrial myopathy and ataxia. Here we report on a 30-yr-old man presenting with a maternally inherited NM_018116.3:c.651C>G, p.F217L missense variant as well as a paternally inherited arr[GRCh37] 1q22(155581773_155706887) × 1 deletion encompassing exons 7-14 of MSTO1 His phenotype included muscle weakness, hypotonia, early motor developmental delay, pectus excavatum, and scoliosis. Testing revealed elevated plasma creatine kinase, and electromyogram results were consistent with longstanding generalized myopathy. These phenotypic features overlap well with previously reported patients harboring biallelic MSTO1 variants. Additionally, our patient presents with dysphagia and restrictive lung disease, not previously reported for MSTO1-associated disorders. The majority of patients with disease-associated variants in MSTO1 present with biallelic variants suggesting autosomal recessive inheritance; however, one family has been reported with a single variant and presumed autosomal dominant inheritance. The pattern of inheritance we observed is consistent with the majority of previous reports suggesting an autosomal recessive disorder. We add to our knowledge of the syndrome caused by variants in MSTO1 and provide additional evidence supporting autosomal recessive inheritance. We also describe phenotypic features not reported in previous cases, although further research is needed to confirm they are associated with defects in MSTO1.

Genetic variants in DGAT1 cause diverse clinical presentations of malnutrition through a specific molecular mechanism

BACKGROUND

DGAT1, a gene encoding a protein involved in lipid metabolism, has been recently implicated in causing a rare nutritional and digestive disease presenting as Congenital Diarrheal Disorder (CDD). Genetic causes of malnutrition can be classified as metabolic disorders, caused by loss of a specific enzyme's function. However, disease driven by genetic variants in lipid metabolism genes is not well understood, and additional information is needed to better understand these effects.

METHODS

We gathered a multi-institutional cohort of undiagnosed patients with a constellation of phenotypes presenting as malnutrition and metal ion dysregulation. Clinical Whole Exome Sequencing (WES) was performed on four patients and their unaffected parents. We prioritized genetic variants based on multiple criteria including population allele frequency and presumed inheritance pattern, and identified a candidate gene. Computational modeling was used to investigate if the altered amino acids are likely to result in a dysfunctional enzyme.

RESULTS

We identified a multi-institutional cohort of patients presenting with malnutrition-like symptoms and likely pathogenic genomic variants within DGAT1. Multiple approaches were used to profile the effect these variants have on protein structure and function. Laboratory and nutritional intervention studies showed rapid and robust patient responses.

CONCLUSIONS

This report adds on to the database for existing mutations known within DGAT1, a gene recently implicated with CDD, and also expands its clinical spectrum. Identification of these DGAT1 mutations by WES has allowed for changes in the patients' nutritional rehabilitation, reversed growth failure and enabled them to be weaned off of total parenteral nutrition (TPN).

A case of Coffin-Siris syndrome with severe congenital heart disease and a novel SMARCA4 variant

Coffin–Siris syndrome (CSS) is a developmental disability, caused by genomic variants in the gene SMARCA4, in addition to other known genes, but the full spectrum of SMARCA4 variants that can cause CSS is unknown with 40% of cases not having molecular confirmation. In this report, we identify a patient with CSS, a severe cardiac phenotype, and a novel SMARCA4 variant. There is no experimental structure of human SMARCA4, so we use molecular modeling techniques to generate a structural model of human SMARCA4. We then map known SMARCA4 variants causative of CSS and our novel variant to the model. We use the resulting information to support the interpretation that the novel variant is causative of disease in our patient. Modeling demonstrates that the variant found in our patient is in a region of SMARCA4 associated with DNA binding, as are the other known pathogenic SMARCA4 variants mapped. Because of this structural information, we discuss how these variants may be disease-causing through a dominant negative effect of disrupting DNA binding.

Discovery of small molecule inhibitors of adenovirus by disrupting E3-19K/HLA-A2 interactions

The binding of the adenovirus (Ad) protein E3-19K with the human leukocyte antigen (HLA) plays an important role in Ad infections, which is the causative agent of a series of gastrointestinal, respiratory and ocular diseases. The objective of this research is to evaluate the essential interactions between E3-19K and HLA-A2 using the X-ray crystal structure of the E3-19K/HLA-A2 complex, and to identify small molecules that could potentially disrupt their binding. Computational methods, including molecular dynamic simulations, MM/GBSA calculations, and computational solvent mapping, were implemented to determine potential binding site(s) for small molecules. The previous experimentally determined hot spot residues, Q54 and E177 in HLA-A2, were also predicted to be the dominant residues for binding to E3-19K by our theoretical calculations. Several other residues were also found to play pivotal roles for the binding of E3-19K with HLA-A2. Residues adjacent to E177, including Q54 and several other residues theoretically predicted to be crucial in HLA-A2 were selected as a potential binding pocket to perform virtual screening with 1200 compounds from the Prestwick library. Seven hits were validated by surface plasmon resonance (SPR) as binders to HLA-A2 as a first step in identifying molecules that can perturb its association with the Ad E3-19K protein.