3D Protein Models Research Articles

Congenital disorder of glycosylation with defective fucosylation 2 (FCSK gene defect): The third report in the literature with a mild phenotype.

Congenital disorders of glycosylation (CDG) are a group of heterogeneous disorders caused by abnormal lipid or protein glycosylation. Variants in the FCSK gene have been reported to cause CDG. Defective FCSK-induced CDG (FCSK-CDG) has only been reported previously in three unrelated children. In this study, we genetically and clinically examined a 3-year-old proband with resolved infantile spasms and normal development. Standard whole-exome sequencing (WES) and Sanger sequencing were performed to identify the functional impact of the variant. WES revealed a rare biallelic missense variant (c.3013G>C; p.Val1005Leu) in FCSK. RT-qPCR showed a significant depletion in FCSK gene expression in the affected individual. Western blotting revealed reduced FCSK expression at the protein level compared to that in the control. Furthermore, 3D protein modeling suggested changes in the secondary structure, which might affect the overall FCSK protein function. This study broadens the mutation and phenotypic spectrum of FCSK-associated developmental disorders.

Genome-wide identification, comprehensive characterization of transcription factors, cis-regulatory elements, protein homology, and protein interaction network of DREB gene family in Solanum lycopersicum.

Tomato is a drought-sensitive crop which has high susceptibility to adverse climatic changes. Dehydration-responsive element-binding (DREB) are significant plant transcription factors that have a vital role in regulating plant abiotic stress tolerance by networking with DRE/CRT cis-regulatory elements in response to stresses. In this study, bioinformatics analysis was performed to conduct the genome-wide identification and characterization of DREB genes and promoter elements in Solanum lycopersicum. In genome-wide coverage, 58 SlDREB genes were discovered on 12 chromosomes that justified the criteria of the presence of AP2 domain as conserved motifs. Intron-exon organization and motif analysis showed consistency with phylogenetic analysis and confirmed the absence of the A3 class, thus dividing the SlDREB genes into five categories. Gene expansion was observed through tandem duplication and segmental duplication gene events in SlDREB genes. Ka/Ks values were calculated in ortholog pairs that indicated divergence time and occurrence of purification selection during the evolutionary period. Synteny analysis demonstrated that 32 out of 58 and 47 out of 58 SlDREB genes were orthologs to Arabidopsis and Solanum tuberosum, respectively. Subcellular localization predicted that SlDREB genes were present in the nucleus and performed primary functions in DNA binding to regulate the transcriptional processes according to gene ontology. Cis-acting regulatory element analysis revealed the presence of 103 motifs in 2.5-kbp upstream promoter sequences of 58 SlDREB genes. Five representative SlDREB proteins were selected from the resultant DREB subgroups for 3D protein modeling through the Phyre2 server. All models confirmed about 90% residues in the favorable region through Ramachandran plot analysis. Moreover, active catalytic sites and occurrence in disorder regions indicated the structural and functional flexibility of SlDREB proteins. Protein association networks through STRING software suggested the potential interactors that belong to different gene families and are involved in regulating similar functional and biological processes. Transcriptome data analysis has revealed that the SlDREB gene family is engaged in defense response against drought and heat stress conditions in tomato. Overall, this comprehensive research reveals the identification and characterization of SlDREB genes that provide potential knowledge for improving abiotic stress tolerance in tomato.

Characteristics of Allan-Herndon-Dudley Syndrome in Chinese children: Identification of two novel pathogenic variants of the SLC16A2 gene.

The aim of this study was to identify causative variants associated with Allan-Herndon-Dudley syndrome (AHDS) in two unrelated Chinese families, and to determine their potential pathogenicity. We also summarized the core clinical symptoms of AHDS by reviewing the related literature. Genomic DNA was isolated from the peripheral blood of AHDS patients and their family members. Whole exome sequencing (WES) was performed on the proband from each family to identify the candidate variants. Subsequently, Sanger sequencing was used to verify the identified candidate variants and to assess co-segregation among the available family members. In silico prediction combined with 3D protein modeling was conducted to predict the functional effects of the variants on the encoded protein. Two novel hemizygous variants of SLC16A2, c.1111_1112insGTCTTGT (Gly375fs*6) and c.942delA (Val315fs*28), were detected in two patients. We compared the clinical symptoms of the patients with all patients with AHDS reported in China and those reported in the literature. While both our patients presented symptoms mostly consistent with AHDS, Patient 1 had no abnormal brain structure and thyroid function, and yet showed other symptoms including lactic aciduria, conjunctival hyperemia, vomiting, laryngeal stridor, low immunoglobulin and iron levels. This study expands the mutation spectrum of AHDS and has clinical value for variant-based prenatal and postnatal screening for this condition. Doctors often have difficulty identifying AHDS by using clinical symptoms. WES can help to identify specific disorder when diagnosis cannot be made based on symptoms alone.

The in silico interaction analysis of CARMIL1 protein-containing leucine-rich repeat (LRR) regions with interleukin-1 receptor-associated kinase 1 (IRAK1) protein and LLR peptide

Objectives: Capping protein Arp2/3 and myosin-I linker protein 1 (CARMIL1) encoded by the CARMIL, is a major, multidomain, membrane-linked protein regulating actin assembly; however, its function in inflammatory signaling is not fully elucidated. The leucine-rich repeat (LRR) region of CARMIL1 has been associated with interleukin (IL)-1 receptor-associated kinase (IRAK) in fibroblasts by many methods including tandem mass tag mass spectrometry, immunoprecipitation, and CRISPR-Cas9. This study, therefore, set out to assess the interaction of CARMIL1 with each IRAK1 protein and a novel LRR peptide. Methods: The molecular docking techniques were employed to compare the binding modes and affinities of the 3D structure of CARMIL1 each of LRR peptides and IRAK1 protein. 3D structure model of CARMIL1 protein and LRR peptide was predicted through Robetta tool considering the structures and function of these proteins. Results: As an overall conclusion of docking, the LRR peptide was observed to contact the residues in the LRR 1-2 of the human CARMIL1, whereas the IRAK1 protein was to interact with the residues in the LRR 1, 2, and 10 regions of the human CARMIL1. Conclusions: Our computational results suggest that LRRs in CARMIL1 are involved in the formation of protein-peptide binding interfaces with its structural conformation.

Characterization of the conserved regions of E1A protein from human adenovirus for reinforcement of cytotoxic T lymphocytes responses to the all genogroups causes ocular manifestation through an in silico approach.

Adenovirus species B, C, D, and E are the most common causes of ocular manifestations caused by adenoviruses. FDA-approved treatment agents for adenovirus infections are not available. Cell-mediated immunity is the major protective mechanism versus human adenoviruses (HAdVs) infection and T cells specific for peptide epitopes from nonstructural proteins can prevent adenoviral dissemination. E1A CR2 region of HAdVs Epitopes predicted for reinforcing cytotoxic T lymphocytes (CTLs) in the EKC patients. Among human adenoviruses E1 protein, four distinct E1A regions had a significantly higher level of homology than the rest of E1A protein. E1A protein inhibits IFN signal transduction. Epitope-based vaccines are designed to have flexible and simple methods to synthesize a vaccine, using an adjuvant to trigger fast immune responses. CTL epitopes were applied to create a multiepitope vaccine. Conserve region1 (CR1) and CR3 have less antigenicity compared to CR2. Additionally, CR3 in HAdV-D8 contains three toxic areas. CR4 similar to the two regions CR1 and CR3 do not show acceptable antigenic properties. Bioinformatics' tools were used to predict, refine and validate the 3D structure of the construct. Effective binding was predicted by protein-protein docking of the epitope vaccine with MHC-I molecules and revealed the safety and efficacy of the predicted vaccine construct. In silico analysis show that rising levels of cytotoxic CD8 + T cells, TH1 cells, macrophages, and neutrophils are linked to IFN-dominant TH1-type responses, which are detected in putative immune individuals. Combined with 3D protein modeling, this study predicted the epitopes of E1A CR2 protein in HAdVs.

Cloning and in silico characterization of an abiotic stress-inducible U-box domain-containing protein gene GsPUB8 from Glycine soja

The ubiquitination pathway is involved in the posttranslational modification of cellular proteins. However, the role of E3 ubiquitin ligase family proteins under abiotic stress conditions remains unclear, particularly in soybean. The core objective of the current study was to isolate and functionally characterize the GsPUB8 protein gene from wild soybean (Glycine soja) by using a homologous cloning method to investigate its abiotic stress responses. The GsPUB8 is a 40,562 Da molecular weight protein with 373 amino acid residues. The sequence alignment revealed the presence of U-box domain while the phylogenetic analysis showed an abundance of PUB8 proteins in both monocot and dicot plants. Analysis of gene structure predicted the absence of introns along with the presence of one exon. Furthermore, the activity of the GsPUB8 protein was anticipated in the plasma membrane and its expression was persuaded with NaCl, ABA, PEG6000, and NaHCO3 treatments with considerably higher manifestation in roots than leaves although, expressed in both vegetative and reproductive parts of G. soja. GsPUB8 protein showed 54% and 32% sequence identity to U-box domain containing 8 and 12 proteins from Arabidopsis thaliana and Oryza sativa subsp. japonica, respectively. GsPUB8 exhibited relatively higher expression under saline and drought stress particularly in roots. Whereas, the 3D model of GsPUB8 protein was generated using the SWISS-MODEL. This study can be used to manipulate the GsPUB8 protein or GsPUB8 gene for transformation purposes and its functional characterization under abiotic stress conditions.

Phylogenetic analysis and characterization of arsenic (As) transforming bacterial marker proteins following isolation of As-tolerant indigenous bacteria.

Marker proteins play a significant role in bacterial arsenic (As) transformation. Phylogenetic analysis and three-dimensional (3D) characteristics of As transforming bacterial marker proteins guide the evolutionary origin and As transforming potential of the species. Indeed, As-tolerant bacteria also show a significant level of As transformation. Hence, characterization of As transforming bacterial marker proteins, isolation of As transforming bacteria, and proper integration of the findings may guide to elucidate how bacteria transform As. Therefore, phylogenetic analysis and 3D characterization of As transforming bacterial marker protein following isolation of potential indigenous As-tolerant indigenous bacteria were done to explore the mechanism of bacterial As transformation. Phylogenetic analysis of ten As transforming marker proteins (arsA, arsB, arsC, arsD, arsR, aioA, arrA, aioB,acr1, andacr3) in 20 potential bacterial genomes (except19 for the acr3) were studied. Some bacterial genomes featured up to five marker proteins, and therefore, 3D characteristics of the marker proteins were analyzed in those genomes having three-to-five marker proteins. In phylogeny, species in close clades represent their phylogenetic resemblances and may have similar functions. P. aeruginosa, E. coli, and K. pneumonia were found to be more effective due to having the highest number (five) of marker proteins. In 3D protein modeling, most of the marker proteins were found to be active. Among 19 indigenous bacterial isolates, multiple isolates showed tolerance up to 50mM As(III) and 250mM As(V), which may potentially transform a significant quantities of As. Hence, integration of the results of phylogenetic analysis, 3D protein characteristics, and As tolerance in the bacterial isolates could guide to explore the mechanism of how bacteria transform As at cellular and molecular levels.

Homozygous Missense Variant in the N-Terminal Region of ANK3 Gene Is Associated with Developmental Delay, Seizures, Speech Abnormality, and Aggressive Behavior

Introduction: Intellectual disability (ID) is a lifelong disability that affects an individual’s learning capacity and adaptive behavior. Such individuals depend on their families for day-to-day survival and pose a significant challenge to the healthcare system, especially in developing countries. ID is a heterogeneous condition, and genetic studies are essential to unravel the underlying cellular pathway for brain development and functioning. Methods: Here we studied a female index patient, born to a consanguineous Pakistani couple, showing clinical symptoms of ID, ataxia, hypotonia, developmental delay, seizures, speech abnormality, and aggressive behavior. Whole exome sequencing (WES) coupled with Sanger sequencing was performed for molecular diagnosis. Further, 3D protein modeling was performed to see the effect of variant on protein structure. Results: WES identified a novel homozygous missense variant (c.178T>C; p.Tyr60His) in the ANK3 gene. In silico analysis and 3-dimensional (3D) protein modeling supports the deleterious impact of this variant on the encoding protein, which compromises the protein’s overall structure and function. Conclusion: Our finding supports the clinical and genetic diversity of the ANK3 gene as a plausible candidate gene for ID syndrome. Intelligence is a complex polygenic human trait, and understanding molecular and biological pathways involved in learning and memory can solve the complex puzzle of how cognition develops. Intellectual disability (ID) is defined as a deficit in an individual’s learning and adaptive behavior at an early age of onset [American Psychiatric Association, 2013]. It is one of the major medical, and cognitive disorders with a prevalence of 1–3% in the population worldwide [Leonard and Wen, 2002]. ID often exists with other disabling mental conditions such as autism, attention deficit hyperactivity disorder, epilepsy, schizophrenia, bipolar disorder, or depression. Almost half of the cases appear to have a genetic explanation that ranges from cytogenetically visible abnormalities to monogenic defects [Flint, 2001; Ropers, 2010; Tucker-Drob et al., 2013]. Intellectual disability is a genetically heterogeneous condition, and more than 700 genes have been identified to cause ID alone or as a part of the syndrome. Research in X-linked ID has identified more than 100 disease-causing genes on the X chromosome that play a role in cognition; however, research into autosomal causes of ID is still ongoing [Vissers et al., 2016].

Oceanimonas sp. BPMS22-derived protein protease inhibitor induces anti-leishmanial immune responses through macrophage M2 to M1 repolarization

The present study aimed to validate the potential of a novel serine protein protease inhibitor (PPI), purified from marine Oceanimonas sp. BPMS22, induced M2 to M1 repolarization of the macrophages to treat visceral leishmaniasis (VL). Peptide mass fingerprint of the purified trypsin digested PPI peptide was obtained using matrix-assisted laser desorption ionization-time of flight combined with tandem mass spectrometry (MALDI-TOF MS/MS) and the sequence was used to construct a 3D protein model by homology modelling. The IC50 of PPI were 25.28 ± 1.675 μg/mL and 0.415 ± 0.015 μg/mL against promastigotes and intracellular amastigotes, respectively, indicating the host-directed therapy using PPI. The PPI enhanced the effector molecule i.e., nitric oxide (NO), and dampened the arginase activity in a dose-dependent manner. In vitro studies revealed that the BPMS22-derived PPI significantly (p < 0.05) decreased the mRNA expressions of M2 markers (FIZZ-1, YM-1, CD206, Arg-1) and increased the mRNA expressions of M1 markers (iNOS, IL-1β, IL-12) in rIL-4 + rIL-10 induced M2 macrophages. Interestingly, the BPMS22-derived PPI also significantly (p < 0.05) decreased the FIZZ-1, YM-1, CD206, and Arg-1; significantly (p < 0.05) increased iNOS, IL-12, and IFN-γ mRNA expression in L. donovani -infected murine macrophages, alongside the decreased parasite load in it. Hence, PPI has the potential to repolarize the cytokines (rIL-4 + rIL-10) pre-stimulated and L. donovani-infected M2 macrophages to M1 phenotype in vitro. A decrease in parasite burden after treatment with PPI indicated the acceleration of the parasite killing by enhancing the macrophage effector functions. Further, in vivo PPI treatment reduced hepatic and splenic Leishman donovan units (LDU) up to 93.34 % and 87.63 %, respectively. This was followed by a surge in pro-inflammatory cytokines and dampening anti-inflammatory cytokines (p < 0.01), which exhibited anti-VL immunity. These observations might open new perspectives on PPI in macrophage repolarization to treat VL.

Autosomal recessive congenital ichthyosis caused by a pathogenic missense variant in CLDN1.

Autosomal recessive congenital ichthyosis (ARCI) refers to a large and genetically heterogenous group of non-syndromic disorders of cornification featuring diffuse scaling. Ichthyosis, leukocyte vacuoles, alopecia, and sclerosing cholangitis (ILVASC) syndrome is a rare autosomal recessive syndromic form of ichthyosis. The disease usually results from premature termination codon-causing pathogenic variants in CLDN1 encoding CLAUDIN-1 (CLDN1). We used whole exome sequencing (WES), Sanger sequencing, 3D protein modeling, Western blotting, and immunofluorescence confocal microscopy to delineate the genetic basis of ichthyosis in two siblings with ichthyosis but no other ectodermal abnormalities. One of the two siblings underwent liver transplantation in early childhood due to biliary atresia. Both patients were found to carry a homozygous missense pathogenic variant, c.242G>A (p.Arg81His), in CLDN1. The variant resulted in decreased CLDN1 expression in patient skin. 3D protein modeling predicted that p.Arg81His induces deleterious conformational changes. Accordingly, HaCaT cells transfected with a construct expressing the mutant CLDN1 cDNA featured decreased levels and mislocation of CLDN1 as compared with cells expressing the wildtype cDNA. In conclusion, we describe the first pathogenic missense variant in CLDN1 shown to result in ARCI.

Prediction and Validation of 3-Dimensional Structure of Rice OsTHIC Abiotic Stress Responsive Gene

Rice is an important cereal crop believed to have been cultivated for hundred years. The ancestry and evolution of this plant and its diverse cultivars remain contentious. Thus, determining the cultivar's molecular mechanism of abiotic stress-responsive genes using a bioinformatics system continues to be an area to investigate. For better knowledge of the mechanisms of the abiotic stress-responsive gene from rice, a novel stress-related gene named OsTHIC was selected for this study. The OsTHIC gene accessions were collected from NCBI, then predict its secondary structure. I-TASSER server was used to predict the 3-dimensional [3D] model of the gene protein and validated using ERRAT and Ramachandran plot. The OsTHIC gene appeared to be highly conserved and structurally functional. The secondary structure of the OsTHIC indicated a high percentage of a random coil [47.35%], followed by an alpha helix [45.05%], then an extended strand [10.61%]. The prediction using I-TASSER modeller produced structures based on the protein sequence. The protein 3D model showed different model qualities using the two-refinement software which ERRAT Plot model refinement having the best score at 97.266%.

Identification of six novel variants from nine Chinese families with hypophosphatemic rickets

BackgroundHypophosphatemic rickets (HR) is a rare genetic disorder associated with renal phosphate wasting and characterized by bone defects. Inactivating mutations in the phosphate regulating endopeptidase homolog X‑linked gene (PHEX) account for most cases of HR. The aim of this study was to identify causative variants in nine unrelated Chinese families associated with HR, and to determine potential pathogenicity of the identified variants.MethodsGenomic DNA was isolated from the peripheral blood of HR patients and their healthy relatives, followed by next-generation sequencing and/or Sanger sequencing. In silico prediction combined with conservation analysis was performed to assess the effects of the variants, and 3D protein modeling was conducted to predict the functional effects on the encoded protein.ResultsAll HR patients recruited in this study displayed bone deformities and tooth agenesis, as well as reduced serum phosphate levels and elevated urine phosphate levels. Nine PHEX variants were identified in eight families, including four novel variants (c.1661_1726del, c.980A > G, c.1078A > T, and c.1017_1051dup). Of the nine identified PHEX variants, five caused a truncated protein, two caused an altered amino acid, and the other two were the canonical splicing variants. Novel variants c.1336G > A and c.1364 T > C in SLC34A3 were also found in one family. Conservation analysis showed that all the amino acids corresponding to the missense variants were highly conserved. In silico analysis and 3D protein structure modeling confirmed the pathogenicity of these variants.ConclusionsThis study identified four novel variants in PHEX and two novel variants in SLC34A3 in a Chinese cohort with HR. Our findings highlight the dominant role of PHEX in HR, and expand the genotypic and phenotypic spectra of this disorder.

Mitochondrial iron transporter (MIT) gene in potato (Solanum tuberosum): comparative bioinformatics, physiological and expression analyses in response to drought and salinity.

Mitochondrial iron transporter (MIT) genes are essential for mitochondrial acquisition/import of iron and vital to proper functioning of mitochondria. Unlike other organisms, research on the MITs in plants is limited. The present study provides comparative bioinformatics assays for the potato MIT gene (StMIT) as well as gene expression analyses. The phylogenetic analyses revealed monocots-dicot divergence in MIT proteins and it was also found clade specific motif diversity. In addition, docking analyses indicated that Asp172 and Gly100 residues to be identified as the closest residues binding to ferrous iron. The percentage of structure overlap of the StMIT 3D protein model with Arabidopsis, maize and rice MIT proteins was found between 80.18% and 85.71%. The transcript analyses exhibited that the expression of StMIT was triggered under drought and salinity stresses. The findings of the present study would provide valuable leads for further studies targeting specifically the MIT gene and generally the plant iron metabolism.

Gain and loss of TASK3 channel function and its regulation by novel variation cause KCNK9 imprinting syndrome

BackgroundGenomics enables individualized diagnosis and treatment, but large challenges remain to functionally interpret rare variants. To date, only one causative variant has been described for KCNK9 imprinting syndrome (KIS). The genotypic and phenotypic spectrum of KIS has yet to be described and the precise mechanism of disease fully understood.MethodsThis study discovers mechanisms underlying KCNK9 imprinting syndrome (KIS) by describing 15 novel KCNK9 alterations from 47 KIS-affected individuals. We use clinical genetics and computer-assisted facial phenotyping to describe the phenotypic spectrum of KIS. We then interrogate the functional effects of the variants in the encoded TASK3 channel using sequence-based analysis, 3D molecular mechanic and dynamic protein modeling, and in vitro electrophysiological and functional methodologies.ResultsWe describe the broader genetic and phenotypic variability for KIS in a cohort of individuals identifying an additional mutational hotspot at p.Arg131 and demonstrating the common features of this neurodevelopmental disorder to include motor and speech delay, intellectual disability, early feeding difficulties, muscular hypotonia, behavioral abnormalities, and dysmorphic features. The computational protein modeling and in vitro electrophysiological studies discover variability of the impact of KCNK9 variants on TASK3 channel function identifying variants causing gain and others causing loss of conductance. The most consistent functional impact of KCNK9 genetic variants, however, was altered channel regulation.ConclusionsThis study extends our understanding of KIS mechanisms demonstrating its complex etiology including gain and loss of channel function and consistent loss of channel regulation. These data are rapidly applicable to diagnostic strategies, as KIS is not identifiable from clinical features alone and thus should be molecularly diagnosed. Furthermore, our data suggests unique therapeutic strategies may be needed to address the specific functional consequences of KCNK9 variation on channel function and regulation.

FAR1/FHY3 Transcription Factors Positively Regulate the Salt and Temperature Stress Responses in Eucalyptus grandis.

FAR-RED ELONGATED HYPOCOTYLS3 (FHY3) and its homolog FAR-RED IMPAIRED RESPONSE1 (FAR1), which play pivotal roles in plant growth and development, are essential for the photo-induced phyA nuclear accumulation and subsequent photoreaction. The FAR1/FHY3 family has been systematically characterized in some plants, but not in Eucalyptus grandis. In this study, genome-wide identification of FAR1/FHY3 genes in E. grandis was performed using bioinformatic methods. The gene structures, chromosomal locations, the encoded protein characteristics, 3D models, phylogenetic relationships, and promoter cis-elements were analyzed with this gene family. A total of 33 FAR1/FHY3 genes were identified in E. grandis, which were divided into three groups based on their phylogenetic relationships. A total of 21 pairs of duplicated repeats were identified by homology analysis. Gene expression analysis showed that most FAR1/FHY3 genes were differentially expressed in a spatial-specific manner. Gene expression analysis also showed that FAR1/FHY3 genes responded to salt and temperature stresses. These results and observation will enhance our understanding of the evolution and function of the FAR1/FHY3 genes in E. grandis and facilitate further studies on the molecular mechanism of the FAR1/FHY3 gene family in growth and development regulations, especially in response to salt and temperature.

Representation of insulin synthesis and construction of a three-dimensional insulin model as didactic tools in basic education

ABSTRACT Proteins are essential for cell function, and their biological performance depends on three-dimensional (3D) structures. Insulin is made up of 51 amino acids and is a critical hormone involved in glycaemia regulation. In basic education, concepts involving proteins (synthesis, structure, interactions, and function) are abstract and difficult for students to understand. The use of protein 3D models that assist in visualising biomolecular structures can be an important didactic tool to understand the complexity of protein structure. Here, we propose a practical activity for the construction of a 3D model of insulin, as a tool for teaching protein synthesis and structure to students in the 9th-grade of basic education. The activity involved three classes of one hour each. 62 students role-played the insulin mRNA translation and represented the primary structure through balloons tied to a string. They also built an 3D insulin model using Styrofoam balls and wire. The activities proved to be effective, as there was a better understanding by students regarding proteins and their characteristics. Therefore, practical classes for teaching proteins can be used in basic education to improve learning of abstract and sub-microscopic concepts.

Impact of protein conformational diversity on AlphaFold predictions.

After the outstanding breakthrough of AlphaFold in predicting protein 3D models, new questions appeared and remain unanswered. The ensemble nature of proteins, for example, challenges the structural prediction methods because the models should represent a set of conformers instead of single structures. The evolutionary and structural features captured by effective deep learning techniques may unveil the information to generate several diverse conformations from a single sequence. Here, we address the performance of AlphaFold2 predictions obtained through ColabFold under this ensemble paradigm. Using a curated collection of apo-holo pairs of conformers, we found that AlphaFold2 predicts the holo form of a protein in ∼70% of the cases, being unable to reproduce the observed conformational diversity with the same error for both conformers. More importantly, we found that AlphaFold2's performance worsens with the increasing conformational diversity of the studied protein. This impairment is related to the heterogeneity in the degree of conformational diversity found between different members of the homologous family of the protein under study. Finally, we found that main-chain flexibility associated with apo-holo pairs of conformers negatively correlates with the predicted local model quality score plDDT, indicating that plDDT values in a single 3D model could be used to infer local conformational changes linked to ligand binding transitions. Data and code used in this manuscript are publicly available at https://gitlab.com/sbgunq/publications/af2confdiv-oct2021. Supplementary data are available at Bioinformatics online.

De novo variants in ATP2B1 lead to neurodevelopmental delay

Calcium (Ca2+) is a universal second messenger involved in synaptogenesis and cell survival; consequently, its regulation is important for neurons. ATPase plasma membrane Ca2+ transporting 1 (ATP2B1) belongs to the family of ATP-driven calmodulin-dependent Ca2+ pumps that participate in the regulation of intracellular free Ca2+. Here, we clinically describe a cohort of 12 unrelated individuals with variants in ATP2B1 and an overlapping phenotype of mild to moderate global development delay. Additional common symptoms include autism, seizures, and distal limb abnormalities. Nine probands harbor missense variants, seven of which were in specific functional domains, and three individuals have nonsense variants. 3D structural protein modeling suggested that the variants have a destabilizing effect on the protein. We performed Ca2+ imaging after introducing all nine missense variants in transfected HEK293 cells and showed that all variants lead to a significant decrease in Ca2+ export capacity compared with the wild-type construct, thus proving their pathogenicity. Furthermore, we observed for the same variant set an incorrect intracellular localization of ATP2B1. The genetic findings and the overlapping phenotype of the probands as well as the functional analyses imply that de novo variants in ATP2B1 lead to a monogenic form of neurodevelopmental disorder.

Feasibility of Follow-Up Studies and Reclassification in Spinocerebellar Ataxia Gene Variants of Unknown Significance.

Spinocerebellar ataxia (SCA) is a heterogeneous group of neurodegenerative disorders with autosomal dominant inheritance. Genetic testing for SCA leads to diagnosis, prognosis and risk assessment for patients and their family members. While advances in sequencing and computing technologies have provided researchers with a rapid expansion in the genetic test content that can be used to unravel the genetic causes that underlie diseases, the large number of variants with unknown significance (VUSes) detected represent challenges. To minimize the proportion of VUSes, follow-up studies are needed to aid in their reclassification as either (likely) pathogenic or (likely) benign variants. In this study, we addressed the challenge of prioritizing VUSes for follow-up using (a combination of) variant segregation studies, 3D protein modeling, in vitro splicing assays and functional assays. Of the 39 VUSes prioritized for further analysis, 13 were eligible for follow up. We were able to reclassify 4 of these VUSes to LP, increasing the molecular diagnostic yield by 1.1%. Reclassification of VUSes remains difficult due to limited possibilities for performing variant segregation studies in the classification process and the limited availability of routine functional tests.

Identification of variants in ACAN and PAPSS2 leading to spondyloepi(meta)physeal dysplasias in four Chinese families.

BackgroundSpondyloepi(meta)physeal dysplasias (SE[M]D) are a group of inherited skeletal disorders that mainly affect bone and cartilage, and next‐generation sequencing has aided the detection of genetic defects of such diseases. In this study, we aimed to identify causative variants in four Chinese families associated with SE(M)D.MethodsWe recruited four unrelated Chinese families all displaying short stature and growth retardation. Clinical manifestations and X‐ray imaging were recorded for all patients. Candidate variants were identified by whole‐exome sequencing (WES) and verified by Sanger sequencing. Pathogenicity was assessed by conservation analysis, 3D protein modeling and in silico prediction, and was confirmed according to American College of Medical Genetics and Genomics.ResultsThree novel SE(M)D‐related variants c.1090dupG, c.7168 T > G, and c.2947G > C in ACAN, and one reported variant c.712C > T in PAPSS2 were identified. Among them, c.1090dupG in ACAN and c.712C > T in PAPSS2 caused truncated protein and the other two variants led to amino acid alterations. Conservation analysis revealed sites of the two missense variants were highly conserved, and bioinformatic findings confirmed their pathogenicity. 3D modeling of mutant protein encoded by c.7168 T > G(p.Trp2390Gly) in ACAN proved the structural alteration in protein level.ConclusionOur data suggested ACAN is a common pathogenic gene of SE(M)D. This study enriched the genetic background of skeletal dysplasias, and expanded the mutation spectra of ACAN and PAPSS2.