Defects in Glabrous 3 (GL3) functionality underlie the absence of trichomes in Brassica napus.

Previously, expression of the Arabidopsis thaliana GLABRA3 (GL3) induced trichome formation in Brassica napus. GL3 orthologues were examined from glabrous (B. oleracea), semi-glabrous (B. napus), moderately hirsute (B. rapa), and very hirsute (B. villosa) Brassica species. Ectopic expression of BnGL3, BrGL3 alleles, or BvGL3 induced trichome formation in glabrous B. napus with the effect on trichome number commensurate with density in the original accessions. Chimeric GL3 proteins in which the B. napus amino terminal region, which interacts with MYB proteins, or the middle region, which interacts with the WD40 protein TTG1, was exchanged with corresponding regions from A. thaliana were as stimulatory to trichome production as AtGL3. Exchange of the carboxy-terminal region containing a bHLH domain and an ACT domain did not alter the trichome stimulatory activity, although modeling of the ACT domain identified differences that could affect GL3 dimerization. B. napus A- and C-genomes orthologues differed in their abilities to form homo- and heterodimers. Modeling of the amino-terminal region revealed a conserved domain that may represent the MYB factor binding pocket. This region interacted with the MYB factors GL1, CPC, and TRY, as well as with JAZ8, which is involved in jasmonic acid-mediated regulation of MYC-like transcription factors. Protein interaction studies indicated that GL1 interaction with GL3 from B. napus and A. thaliana may underlie the difference in their respective abilities to induce trichome formation.

Cytonuclear interplay in auto- and allopolyploids: a multifaceted perspective from the Festuca-Lolium complex.

Restoring cytonuclear stoichiometry is necessary after whole-genome duplication (WGD) and interspecific/intergeneric hybridization in plants. We investigated this phenomenon in auto- and allopolyploids of the Festuca-Lolium complex providing insights into the mechanisms governing cytonuclear interactions in early polyploid and hybrid generations. Our study examined the main processes potentially involved in restoring the cytonuclear balance after WGD comparing diploids and new and well-established autopolyploids. We uncovered that both the number of chloroplasts and the number of chloroplast genome copies were significantly higher in the newly established autopolyploids and grew further in more established autopolyploids. The increase in the copy number of the chloroplast genome exceeded the rise in the number of chloroplasts and fully compensated for the doubling of the nuclear genome. In addition, changes in nuclear and organelle gene expression were insignificant. Allopolyploid Festuca × Lolium hybrids displayed potential structural conflicts in parental protein variants within the cytonuclear complexes. While biased maternal allele expression has been observed in numerous hybrids, our results suggest that its role in cytonuclear stabilization in the Festuca × Lolium hybrids is limited. This study provides insights into the restoration of the cytonuclear stoichiometry, yet it emphasizes the need for future research to explore post-transcriptional regulation and its impact on cytonuclear gene expression stoichiometry. Our findings may enhance the understanding of polyploid plant evolution, with broader implications for the study of cytonuclear interactions in diverse biological contexts.

Selective cross-linking of coinciding protein assemblies by in-gel cross-linking mass spectrometry.

Cross-linking mass spectrometry has developed into an important method to study protein structures and interactions. The in-solution cross-linking workflows involve time and sample consuming steps and do not provide sensible solutions for differentiating cross-links obtained from co-occurring protein oligomers, complexes, or conformers. Here we developed a cross-linking workflow combining blue native PAGE with in-gel cross-linking mass spectrometry (IGX-MS). This workflow circumvents steps, such as buffer exchange and cross-linker concentration optimization. Additionally, IGX-MS enables the parallel analysis of co-occurring protein complexes using only small amounts of sample. Another benefit of IGX-MS, demonstrated by experiments on GroEL and purified bovine heart mitochondria, is the substantial reduction of undesired over-length cross-links compared to in-solution cross-linking. We next used IGX-MS to investigate the complement components C5, C6, and their hetero-dimeric C5b6 complex. The obtained cross-links were used to generate a refined structural model of the complement component C6, resembling C6 in its inactivated state. This finding shows that IGX-MS can provide new insights into the initial stages of the terminal complement pathway.

GraphGPSM: a global scoring model for protein structure using graph neural networks.

The scoring models used for protein structure modeling and ranking are mainly divided into unified field and protein-specific scoring functions. Although protein structure prediction has made tremendous progress since CASP14, the modeling accuracy still cannot meet the requirements to a certain extent. Especially, accurate modeling of multi-domain and orphan proteins remains a challenge. Therefore, an accurate and efficient protein scoring model should be developed urgently to guide the protein structure folding or ranking through deep learning. In this work, we propose a protein structure global scoring model based on equivariant graph neural network (EGNN), named GraphGPSM, to guide protein structure modeling and ranking. We construct an EGNN architecture, and a message passing mechanism is designed to update and transmit information between nodes and edges of the graph. Finally, the global score of the protein model is output through a multilayer perceptron. Residue-level ultrafast shape recognition is used to describe the relationship between residues and the overall structure topology, and distance and direction encoded by Gaussian radial basis functions are designed to represent the overall topology of the protein backbone. These two features are combined with Rosetta energy terms, backbone dihedral angles and inter-residue distance and orientations to represent the protein model and embedded into the nodes and edges of the graph neural network. The experimental results on the CASP13, CASP14 and CAMEO test sets show that the scores of our developed GraphGPSM have a strong correlation with the TM-score of the models, which are significantly better than those of the unified field score function REF2015 and the state-of-the-art local lDDT-based scoring models ModFOLD8, ProQ3D and DeepAccNet, etc. The modeling experimental results on 484 test proteins demonstrate that GraphGPSM can greatly improve the modeling accuracy. GraphGPSM is further used to model 35 orphan proteins and 57 multi-domain proteins. The results show that the average TM-score of the models predicted by GraphGPSM is 13.2 and 7.1% higher than that of the models predicted by AlphaFold2. GraphGPSM also participates in CASP15 and achieves competitive performance in global accuracy estimation.

Multidimensional Cross-Linking and Real-Time Informatics for Multiprotein Interaction Studies.

Chemical cross-linking combined with mass spectrometry is a technique used to study protein structures and identify protein complexes. Traditionally, chemical cross-linkers contain two reactive groups, allowing them to covalently bond a pair of proximal residues, either within a protein or between two proteins. The output of a cross-linking experiment is a list of interacting site pairs that provide structural constraints for modeling of new structures and complexes. Due to the binary reactive nature of cross-linking reagents, only pairs of interacting sites can be directly observed, and assembly of higher-order structures typically requires prior knowledge of complex composition or iterative docking to produce a putative model. Here, we describe a new tetrameric cross-linker bearing four amine-reactive groups, allowing it to covalently link up to four proteins simultaneously and a real-time instrument method to facilitate the identification of these tetrameric cross-links. We applied this new cross-linker to isolated mitochondria and identified a number of higher-order cross-links in various OXPHOS complexes and ATP synthase, demonstrating its utility in characterizing complex interfaces. We also show that higher-order cross-links can be used to effectively filter models of large protein assemblies generated by using Alphafold. Higher-dimensional cross-linking provides a new avenue for characterizing multiple protein interfaces, even in complex samples such as intact mitochondria.

Invariant point message passing for protein side chain packing.

Protein side chain packing (PSCP) is a fundamental problem in the field of protein engineering, as high-confidence and low-energy conformations of amino acid side chains are crucial for understanding (and designing) protein folding, protein-protein interactions, and protein-ligand interactions. Traditional PSCP methods (such as the Rosetta Packer) often rely on a library of discrete side chain conformations, or rotamers, and a forcefield to guide the structure to low-energy conformations. Recently, deep learning (DL) based methods (such as DLPacker, AttnPacker, and DiffPack) have demonstrated state-of-the-art predictions and speed in the PSCP task. Building off the success of geometric graph neural networks for protein modeling, we present the Protein Invariant Point Packer (PIPPack) which effectively processes local structural and sequence information to produce realistic, idealized side chain coordinates using χ $$ \chi $$ -angle distribution predictions and geometry-aware invariant point message passing (IPMP). On a test set of ∼1400 high-quality protein chains, PIPPack is highly competitive with other state-of-the-art PSCP methods in rotamer recovery and per-residue RMSD but is significantly faster.

SpotLink enables sensitive and precise identification of site nonspecific cross-links at the proteome scale.

Nonspecific cross-linker can provide distance restraints between surface residues of any type, which could be used to investigate protein structure construction and protein-protein interaction (PPI). However, the vast number of potential combinations of cross-linked residues or sites obtained with such a cross-linker makes the data challenging to analyze, especially for the proteome-wide applications. Here, we developed SpotLink software for identifying site nonspecific cross-links at the proteome scale. Contributed by the dual pointer dynamic pruning algorithm and the quality control of cross-linking sites, SpotLink identified > 3000 cross-links from human cell samples within a short period of days. We demonstrated that SpotLink outperformed other approaches in terms of sensitivity and precision on the datasets of the simulated succinimidyl 4,4'-azipentanoate dataset and the condensin complexes with known structures. In addition, some valuable PPI were discovered in the datasets of the condensin complexes and the HeLa dataset, indicating the unique identification advantages of site nonspecific cross-linking. These findings reinforce the importance of SpotLink as a fundamental characteristic of site nonspecific cross-linking technologies.

Comparative studies of AlphaFold, RoseTTAFold and Modeller: a case study involving the use of G-protein-coupled receptors.

Neural network (NN)-based protein modeling methods have improved significantly in recent years. Although the overall accuracy of the two non-homology-based modeling methods, AlphaFold and RoseTTAFold, is outstanding, their performance for specific protein families has remained unexamined. G-protein-coupled receptor (GPCR) proteins are particularly interesting since they are involved in numerous pathways. This work directly compares the performance of these novel deep learning-based protein modeling methods for GPCRs with the most widely used template-based software-Modeller. We collected the experimentally determined structures of 73 GPCRs from the Protein Data Bank. The official AlphaFold repository and RoseTTAFold web service were used with default settings to predict five structures of each protein sequence. The predicted models were then aligned with the experimentally solved structures and evaluated by the root-mean-square deviation (RMSD) metric. If only looking at each program's top-scored structure, Modeller had the smallest average modeling RMSD of 2.17 Å, which is better than AlphaFold's 5.53 Å and RoseTTAFold's 6.28 Å, probably since Modeller already included many known structures as templates. However, the NN-based methods (AlphaFold and RoseTTAFold) outperformed Modeller in 21 and 15 out of the 73 cases with the top-scored model, respectively, where no good templates were available for Modeller. The larger RMSD values generated by the NN-based methods were primarily due to the differences in loop prediction compared to the crystal structures.

Regulation of Ras p21 and RalA GTPases activity by quinine in mammary epithelial cells.

Quinine, a bitter compound, can act as an agonist to activate the family of bitter taste G protein-coupled receptor family of proteins. Previous work from our laboratory has demonstrated that quinine causes activation of RalA, a Ras p21-related small G protein. Ral proteins can be activated directly or indirectly through an alternative pathway that requires Ras p21 activation resulting in the recruitment of RalGDS, a guanine nucleotide exchange factor for Ral. Using normal mammary epithelial (MCF-10A) and non-invasive mammary epithelial (MCF-7) cell lines, we investigated the effect of quinine in regulating Ras p21 and RalA activity. Results showed that in the presence of quinine, Ras p21 is activated in both MCF-10A and MCF-7 cells; however, RalA was inhibited in MCF-10A cells, and no effect was observed in the case of MCF-7 cells. MAP kinase, a downstream effector for Ras p21, was activated in both MCF-10A and MCF-7 cells. Western blot analysis confirmed the expression of RalGDS in MCF-10A cells and MCF-7 cells. The expression of RalGDS was higher in MCF-10A cells in comparison to the MCF-7 cells. Although RalGDS was detected in MCF-10A and MCF-7 cells, it did not result in RalA activation upon Ras p21 activation with quinine suggesting that the Ras p21-RalGDS-RalA pathway is not active in the MCF-10A cells. The inhibition of RalA activity in MCF-10A cells due to quinine could be as a result of a direct effect of this bitter compound on RalA. Protein modeling and ligand docking analysis demonstrated that quinine can interact with RalA through the R79 amino acid, which is located in the switch II region loop of the RalA protein. It is possible that quinine causes a conformational change that results in the inhibition of RalA activation even though RalGDS is present in the cell. More studies are needed to elucidate the mechanism(s) that regulate Ral activity in mammary epithelial cells.

Molecular genetics, neuroimaging outcomes, and structural analyses of novel and recurrent variants of WDR62 gene in two consanguineous Pakistani families with autosomal recessive primary microcephaly.

BACKGROUND: Autosomal recessive primary microcephaly (MCPH) is a rare neurodevelopmental and genetically heterogeneous disorder, characterized by small cranium size (> - 3 SD below mean) and often results in varying degree of intellectual disability. Thirty genes have been identified for the etiology of this disorder due to its clinical and genetic heterogeneity. METHODS AND RESULTS: Here, we report two consanguineous Pakistani families affected with MCPH exhibiting mutation in WDR62 gene. The investigation approach involved Next Generation Sequencing (NGS) gene panel sequencing coupled with linkage analysis followed by validation of identified variants through automated Sanger sequencing and Barcode-Tagged (BT) sequencing. The molecular genetic analysis revealed one novel splice site variant (NM_001083961.2(WDR62):c.1372-1del) in Family A and one known exonic variant NM_001083961.2(WDR62):c.3936dup (p.Val1313Argfs*18) in Family B. Magnetic Resonance Imaging (MRI) scans were also employed to gain insights into the structural architecture of affected individuals. Neurological assessments showed the reduced gyral and sulcal patterns along with normal corpus callosum in affected individuals harboring novel variant. In silico assessments of the identified variants were conducted using different tools to confirm the pathogenicity of these variants. Through In silico analyses, both variants were identified as disease causing and protein modeling of exonic variant indicates subtle conformational alterations in prophesied protein structure. CONCLUSION: This study identifies a novel variant (c.1372-1del) and a recurrent pathogenic variant c.3936dup (p.Val1313Argfs*18) in the WDR62 gene among the Pakistani population, expanding the mutation spectrum for MCPH. These findings emphasize the importance of genetic counseling and awareness to reduce consanguinity and address the burden of this disorder.

Molecular Prediction and Correlation of the Structure and Function of Universal Stress Protein A (UspA) from Salmonella Typhimurium.

Salmonella Typhimurium (ST) is a zoonotic pathogen that can cause gastroenteritis in humans when they consume contaminated food or water. When exposed to various stressors, both from living organisms (biotic) and the environment (abiotic), Salmonella Typhimurium produces Universal Stress Proteins (USPs). These proteins are gaining recognition for their crucial role in bacterial stress resistance and the ability to enter a prolonged state of growth arrest. Additionally, USPs exhibit diverse structures due to the fusion of the USP domain with different catalytic motifs, enabling them to participate in various reactions and cellular activities during stressful conditions. In this particular study, researchers cloned and analyzed the uspA gene obtained from poultry-derived strains of Salmonella Typhimurium. The gene comprises 435 base pairs, encoding a USP family protein consisting of 144 amino acids. Phylogenetic analysis demonstrated a close relationship between the uspA genes of Salmonella Typhimurium and those found in other bacterial species. We used molecular dynamics simulations and 3D structure prediction to ensure that the USPA protein was stable. Furthermore, we also carried out motif search and network analysis of protein-protein interactions. The findings from this study offer valuable insights for the development of inhibitors targeted against Salmonella Typhimurium.

Exploring the Structurally Conserved Regions and Functional Significance in Bacterial N-Terminal Nucleophile (Ntn) Amide-Hydrolases.

The initial adoption of penicillin as an antibiotic marked the start of exploring other compounds essential for pharmaceuticals, yet resistance to penicillins and their side effects has compromised their efficacy. The N-terminal nucleophile (Ntn) amide-hydrolases S45 family plays a key role in catalyzing amide bond hydrolysis in various compounds, including antibiotics like penicillin and cephalosporin. This study comprehensively analyzes the structural and functional traits of the bacterial N-terminal nucleophile (Ntn) amide-hydrolases S45 family, covering penicillin G acylases, cephalosporin acylases, and D-succinylase. Utilizing structural bioinformatics tools and sequence analysis, the investigation delineates structurally conserved regions (SCRs) and substrate binding site variations among these enzymes. Notably, sixteen SCRs crucial for substrate interaction are identified solely through sequence analysis, emphasizing the significance of sequence data in characterizing functionally relevant regions. These findings introduce a novel approach for identifying targets to enhance the biocatalytic properties of N-terminal nucleophile (Ntn) amide-hydrolases, while facilitating the development of more accurate three-dimensional models, particularly for enzymes lacking structural data. Overall, this research advances our understanding of structure-function relationships in bacterial N-terminal nucleophile (Ntn) amide-hydrolases, providing insights into strategies for optimizing their enzymatic capabilities.

AI-Driven Deep Learning Techniques in Protein Structure Prediction.

Protein structure prediction is important for understanding their function and behavior. This review study presents a comprehensive review of the computational models used in predicting protein structure. It covers the progression from established protein modeling to state-of-the-art artificial intelligence (AI) frameworks. The paper will start with a brief introduction to protein structures, protein modeling, and AI. The section on established protein modeling will discuss homology modeling, ab initio modeling, and threading. The next section is deep learning-based models. It introduces some state-of-the-art AI models, such as AlphaFold (AlphaFold, AlphaFold2, AlphaFold3), RoseTTAFold, ProteinBERT, etc. This section also discusses how AI techniques have been integrated into established frameworks like Swiss-Model, Rosetta, and I-TASSER. The model performance is compared using the rankings of CASP14 (Critical Assessment of Structure Prediction) and CASP15. CASP16 is ongoing, and its results are not included in this review. Continuous Automated Model EvaluatiOn (CAMEO) complements the biennial CASP experiment. Template modeling score (TM-score), global distance test total score (GDT_TS), and Local Distance Difference Test (lDDT) score are discussed too. This paper then acknowledges the ongoing difficulties in predicting protein structure and emphasizes the necessity of additional searches like dynamic protein behavior, conformational changes, and protein-protein interactions. In the application section, this paper introduces some applications in various fields like drug design, industry, education, and novel protein development. In summary, this paper provides a comprehensive overview of the latest advancements in established protein modeling and deep learning-based models for protein structure predictions. It emphasizes the significant advancements achieved by AI and identifies potential areas for further investigation.

In silico modeling human VPS13 proteins associated with donor and target membranes suggests lipid transfer mechanisms.

The VPS13 protein family constitutes a novel class of bridge-like lipid transferases. Autosomal recessive inheritance of mutations in VPS13 genes is associated with the development of neurodegenerative diseases in humans. Bioinformatic approaches previously recognized the domain architecture of these proteins. In this study, we model the first ever full-length structures of the four human homologs VPS13A, VPS13B, VPS13C, and VPS13D in association with model membranes, to investigate their lipid transfer ability and potential structural association with membrane leaflets. We analyze the evolutionary conservation and physicochemical properties of these proteins, focusing on conserved C-terminal amphipathic helices that disturb organelle surfaces and that, adjoined, resemble a traditional Venetian gondola. The gondola domains share significant structural homology with lipid droplet surface-binding proteins. We introduce in silico protein-membrane models displaying the mode of association of VPS13A, VPS13B, VPS13C, and VPS13D to donor and target membranes, and present potential models of action for protein-mediated lipid transfer.

Amphiphilic peptide Mastoparan-B induces conformational changes within the AdeB efflux pump, down-regulates adeB gene expression, and restores antibiotic susceptibility in an MDR strain of Acinetobacter baumannii.

Mastoparan B (MP-B) is an amphiphilic peptide with a potent antimicrobial activity against most Gram-negative bacteria. However, there is little information available on the inhibition of the Acinetobacter baumannii resistance-nodulation-cell-division (RND) efflux pump using this antimicrobial peptide. Here, we carried out a series of in-silico experiments to find the mechanisms underlying the anti-efflux activity of MP-B using a multi-drug resistant (MDR) strain of A. baumannii (AB). According to our findings, MP-B demonstrated a potent antibacterial activity against an MDR-AB (minimum inhibitory concentration [MIC] = 1 µg/mL) followed by a 20-fold reduction in the adeB gene expression in the presence of sub-MIC of this peptide. Using Groningen Machine for Chemicals Simulation (GROMACS) via PyMOL Graphical User Interface (GUI), (we observed that, the AdeB transporter had conserved helix-turn-helix regions and a tight pore rich in Phe and Ala residues. To understand how inhibition of the AdeB is achieved, we generated 20 apo-MP-B poses using the InterPep and SiteMap tools. The high-quality model was created by homology modeling and used for docking via AutoDock/Vina to identify the MP-B binding sites. We established that the most apo-MP-B formed H-bonds to the backbone of five amino acids in the Helix-5. As a result, the dihedral angles of the involved amino acids shift by 9.0-9.6 Ǻ, causing a change in the conformation of the AdeB protein. This led to helix conformation stereoisomerization and block the AdeB activity. MP-B presumably has dual mechanisms. (1) It blocks the AdeB transporter by changing its conformation. (2) MP-B influences the adeB gene expression by binding to G-protein which laterally controls efflux regulators like MarA, RamA, SoxS, and Rob proteins.

Cloning, Expression and Functional Characterization of a Novel α-Humulene Synthase, Responsible for the Formation of Sesquiterpene in Agarwood Originating from Aquilaria malaccensis.

This study describes the cloning, expression and functional characterization of α-humulene synthase, responsible for the formation of the key aromatic compound α-humulene in agarwood originating from Aquilaria malaccensis. The partial sesquiterpene synthase gene from the transcriptome data of A. malaccensis was utilized for full-length gene isolation via a 3' RACE PCR. The complete gene, denoted as AmDG2, has an open reading frame (ORF) of 1671 bp and encodes for a polypeptide of 556 amino acids. In silico analysis of the protein highlighted several conserved motifs typically found in terpene synthases such as Asp-rich substrate binding (DDxxD), metal-binding residues (NSE/DTE), and cytoplasmic ER retention (RxR) motifs at their respective sites. The AmDG2 was successfully expressed in the E. coli:pET-28a(+) expression vector whereby an expected band of about 64 kDa in size was detected in the SDS-PAGE gel. In vitro enzyme assay using substrate farnesyl pyrophosphate (FPP) revealed that AmDG2 gave rise to two sesquiterpenes: α-humulene (major) and ß-caryophyllene (minor), affirming its identity as α-humulene synthase. On the other hand, protein modeling performed using AlphaFold2 suggested that AmDG2 consists entirely of α-helices with short connecting loops and turns. Meanwhile, molecular docking via AutoDock Vina (Version 1.5.7) predicted that Asp307 and Asp311 act as catalytic residues in the α-humulene synthase. To our knowledge, this is the first comprehensive report on the cloning, expression and functional characterization of α-humulene synthase from agarwood originating from A. malaccensis species. These findings reveal a deeper understanding of the structure and functional properties of the α-humulene synthase and could be utilized for metabolic engineering work in the future.

Advanced computational analysis of CD40LG variants in atypical X-linked hyper-IgM syndrome.

X-Linked Hyper-IgM Syndrome is caused by pathogenic variants in CD40LG. Three patients with atypical clinical and immunological features were identified with variants in CD40LG requiring further characterization. Flow cytometry was used to evaluate CD40L protein expression and binding capacity to a surrogate receptor, CD40-muIg. Though functional anomalies were observed, there was still a lack of clarity regarding the underlying mechanism. We developed structural models for wild-type and the three variants of CD40L protein observed in these patients (p. Lys143Asn, Leu225Ser and Met36Arg) to evaluate structural alterations by molecular mechanic calculations, and assess protein movement by molecular dynamic simulations. These studies demonstrate that functional analysis of variants of unknown significance in CD40LG can be supplemented by advanced computational analysis in atypical clinical contexts. These studies in combination identify the deleterious effects of these variants and potential mechanisms for protein dysfunction.

Characterization of POLE c.1373A > T p.(Tyr458Phe), causing high cancer risk.

The cancer syndrome polymerase proofreading-associated polyposis results from germline mutations in the POLE and POLD1 genes. Mutations in the exonuclease domain of these genes are associated with hyper- and ultra-mutated tumors with a predominance of base substitutions resulting from faulty proofreading during DNA replication. When a new variant is identified by gene testing of POLE and POLD1, it is important to verify whether the variant is associated with PPAP or not, to guide genetic counseling of mutation carriers. In 2015, we reported the likely pathogenic (class 4) germline POLE c.1373A > T p.(Tyr458Phe) variant and we have now characterized this variant to verify that it is a class 5 pathogenic variant. For this purpose, we investigated (1) mutator phenotype in tumors from two carriers, (2) mutation frequency in cell-based mutagenesis assays, and (3) structural consequences based on protein modeling. Whole-exome sequencing of two tumors identified an ultra-mutator phenotype with a predominance of base substitutions, the majority of which are C > T. A SupF mutagenesis assay revealed increased mutation frequency in cells overexpressing the variant of interest as well as in isogenic cells encoding the variant. Moreover, exonuclease repair yeast-based assay supported defect in proofreading activity. Lastly, we present a homology model of human POLE to demonstrate structural consequences leading to pathogenic impact of the p.(Tyr458Phe) mutation. The three lines of evidence, taken together with updated co-segregation and previously published data, allow the germline variant POLE c.1373A > T p.(Tyr458Phe) to be reclassified as a class 5 variant. That means the variant is associated with PPAP.

A member of the OSCA/TMEM63 family of mechanosensitive calcium channels participates in cell wall integrity maintenance in Aspergillus nidulans.

The calF7 mutation in Aspergillus nidulans causes hypersensitivity to the cell wall compromising agents Calcofluor White (CFW) and Congo Red. In this research we demonstrate that the calF7 mutation resides in gene AN2880, encoding a predicted member of the OSCA/TMEM63 family of transmembrane glycoproteins. Those members of the family whose physiological functions have been investigated have been shown to act as mechanosensitive calcium transport channels. Deletion of AN2880 replicates the CFW hypersensitivity phenotype. Separately, we show that CFW hypersensitivity of calF deletion strains can be overcome by inclusion of elevated levels of extracellular calcium ions in the growth medium, and, correspondingly, wild type strains grown in media deficient in calcium ions are no longer resistant to CFW. These observations support a model in which accommodation to at least some forms of cell wall stress is mediated by a calcium ion signaling system in which the AN2880 gene product plays a role. The genetic lesion in calF7 is predicted to result in a glycine-to-arginine substitution at position 638 of the 945-residue CalF protein in a region of the RSN1_7TM domain that is highly conserved amongst filamentous fungi. Homology modeling predicts that the consequence of a G638R substitution is to structurally occlude the principal conductance pore in the protein. GFP-tagged wild type CalF localizes principally to the Spitzenkörper and the plasma membrane at growing tips and forming septa. However, both septation and hyphal morphology appear to be normal in calF7 and AN2880 deletion strains, indicating that any role played by CalF in normal hyphal growth and cytokinesis is dispensable.

A Novel De Novo Frameshift Pathogenic Variant in the FAM111B Resulting in Progressive Osseous Heteroplasia Phenotype.

INTRODUCTION: Progressive osseous heteroplasia (POH) is a rare, debilitating disorder characterized by heterotopic ossification in the skin and muscles, resulting in contractures of the joints and progressive loss of function. While 60-70% of the POH patients have paternally inherited, inactivating pathogenic variants in GNAS, the remaining 30-40% have no known etiology. FAM111B pathogenic variants, located on chromosome 11q12.1, cause POIKTMP (hereditary fibrosing poikiloderma with tendon contractures, myopathy, and pulmonary fibrosis), a very rare, autosomal-dominant disorder with high frequency of de novo missense pathogenic variants, which affects multiple tissues and organs, causing extensive fibrosis and muscle adiposis, though the exact mechanism is unknown. To our knowledge, there are no reports of FAM111B associated with POH. We describe the first case of POH phenotype associated with a novel de novo frameshift pathogenic variant in the FAM111B and present an analysis of the protein structure and function caused by this genomic disruption. CASE: A 15-year-old African-American male presented with generalized calcific nodules, progressive contractures, and muscle weakness leading to immobility, beginning at 6 years of age. Cutaneous examination showed generalized hard nodules varying from small to plaque-like ulcerated erupted skin lesions. Biochemical evaluation revealed 25(OH) vitamin D insufficiency (20 ng/mL), and normal levels of parathyroid hormone, FGF-23, alkaline phosphatase, calcium, and phosphorus. Skeletal survey radiographs and computed tomography (CT) of the chest, abdomen, and pelvis showed extensive soft tissue and muscle heterotopic ossifications involving shoulders, axillae, trunk, abdomen, pelvis, upper and lower extremities, in a clumped, conglomerate distribution within muscle, subcutaneous fat, and in some areas extending to the skin. There was no pulmonary fibrosis on the chest CT. The clinical and radiographic findings were most consistent with POH. A trio-clinical exome sequencing revealed a de novo heterozygous likely pathogenic variant in the FAM111B (OMIM # 615584) (c.1462delT [p.Cys488Valfs*21]). The resulted frameshift change in exon 4 replaced C-terminal region with 21 alternative amino acids. Multiple, previously reported disease-associated variants appear to localize within the trypsin-like cysteine/serine peptidase domain in which this variant occurs, supporting the functional significance of this region, though none have been previously reported to be associated with POH phenotype. Our 3D protein modeling showed obliteration of predicted protein folding and structure, and elimination of the zinc-binding domain, likely severely affecting protein function. CONCLUSION: This is the first case of POH phenotype associated with a novel de novo pathogenic frameshift variant in FAM111B. Whether the frameshift change in FAM111B predicts POH remains unclear. Further evaluations are necessary to fully elucidate this finding and the potential role and mechanism by which the FAM111B variants contributes to POH phenotype.