Structural and Dynamic Analyses of Pathogenic Variants in PIK3R1 Reveal a Shared Mechanism Associated among Cancer, Undergrowth, and Overgrowth Syndromes.

The PI3K enzymes modify phospholipids to regulate cell growth and differentiation. Somatic variants in PI3K are recurrent in cancer and drive a proliferative phenotype. Somatic mosaicism of PIK3R1 and PIK3CA are associated with vascular anomalies and overgrowth syndromes. Germline PIK3R1 variants are associated with varying phenotypes, including immunodeficiency or facial dysmorphism with growth delay, lipoatrophy, and insulin resistance associated with SHORT syndrome. There has been limited study of the molecular mechanism to unify our understanding of how variants in PIK3R1 drive both undergrowth and overgrowth phenotypes. Thus, we compiled genomic variants from cancer and rare vascular anomalies and sought to interpret their effects using an unbiased physics-based simulation approach for the protein complex. We applied molecular dynamics simulations to mechanistically understand how genetic variants affect PIK3R1 and its interactions with PIK3CA. Notably, iSH2 genetic variants associated with undergrowth destabilize molecular interactions with the PIK3CA receptor binding domain in simulations, which is expected to decrease activity. On the other hand, overgrowth and cancer variants lead to loss of inhibitory interactions in simulations, which is expected to increase activity. We find that all disease variants display dysfunctions on either structural characteristics or intermolecular interaction energy. Thus, this comprehensive characterization of novel mosaic somatic variants associated with two opposing phenotypes has mechanistic importance and biomedical relevance and may aid in future therapeutic developments.

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.

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.

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.

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.