Spectrum of disease-causing mutations in protein secondary structures

Metabolic and hormonal profiling in polycystic ovarian syndrome: insights into INSR gene variations

BACKGROUND

Polycystic Ovary Syndrome (PCOS) is a hormonal disorder characterized by irregular periods, excess androgen levels, and polycystic ovaries, affecting many women of reproductive age.

METHODS AND RESULTS

This study employed statistical and molecular analyses to compare hormone and metabolic markers between PCOS patients and controls. Sanger sequencing identified two INSR gene variants linked to high insulin and pre-diabetic conditions. Statistically, no significant age differences were detected (p = 0.492) between the overall PCOS patient pool and controls. However, a substantial variation in Vitamin D levels was observed within PCOS patients compared to controls (p = 0.0006), suggesting an association with PCOS. Correlations between Vitamin D and insulin, as well as HbA1c levels (R2 = 0.141 and 0.143, respectively), suggest Vitamin D's potential impact on glycemic control. Significant differences were found in HbA1c (p < 0.0001), insulin (p < 0.0001), and LDL (p = 0.0004) levels between PCOS patients and controls, highlighting marked disparities in these metabolic markers. LH levels also showed a significant contrast (p < 0.0001), while progesterone levels displayed a notable difference (p = 0.007) between the two groups. Correlation analyses within PCOS patients demonstrated associations among LDL, HbA1c, and insulin, with no such correlations observed in control cases. Additionally, Sanger sequencing identified two INSR gene variants, c.3614C > T (p.Pro1205Leu) and c.3355C > T (p.Arg1119Trp), associated with high insulin, LH, and pre-diabetic conditions. These amino acid changes may trigger metabolic imbalances and hormonal irregularities, potentially contributing to the development of PCOS.

CONCLUSIONS

The findings highlight the multifaceted nature of PCOS, revealing significant metabolic, hormonal, and genetic differences compared to controls. These insights may inform tailored interventions and management strategies for the complex associations characteristic of PCOS.

Breast Cancer is Increased in Women with Primary Ovarian Insufficiency

CONTEXT

DNA damage/repair gene variants are associated with both primary ovarian insufficiency (POI) and cancer risk.

OBJECTIVE

We hypothesized that a subset of women with POI and family members would have increased risk for cancer.

DESIGN

Case-control population-based study using records from 1995-2022.

SETTING

Two major Utah academic healthcare systems serving 85% of the state.

SUBJECTS

Women with POI (n=613) were identified using ICD codes and reviewed for accuracy. Relatives were linked using the Utah Population Database.

INTERVENTION

Cancer diagnoses were identified using the Utah Cancer Registry.

MAIN OUTCOME MEASURES

The relative risk of cancer in women with POI and relatives was estimated by comparison to population rates. Whole genome sequencing was performed on a subset of women.

RESULTS

Breast cancer was increased in women with POI (OR [95%CI] 2.20 [1.30, 3.47]; p=0.0023) and there was a nominally significant increase in ovarian cancer. Probands with POI were 36.5±4.3 years and 59.5±12.7 years when diagnosed with POI and cancer, respectively. Causal and candidate gene variants for cancer and POI were identified.Among second-degree relatives of these women, there was an increased risk of breast (1.28 [1.08, 1.52]; p=0.0078) and colon cancer (1.50 [1.14, 1.94]; p=0.0036). Prostate cancer was increased in first- (1.64 [1.18, 2.23]; p=0.0026), second- (1.54 [1.32, 1.79]; p<0.001), and third-degree relatives (1.33 [1.20, 1.48]; p<0.001).

CONCLUSIONS

Data suggest common genetic risk for POI and reproductive cancers. Tools are needed to predict cancer risk in women with POI and potentially to counsel about risks of hormone replacement therapy.

In silico functional, structural and pathogenicity analysis of missense single nucleotide polymorphisms in human MCM6 gene

Single nucleotide polymorphisms (SNPs) are one of the most common determinants and potential biomarkers of human disease pathogenesis. SNPs could alter amino acid residues, leading to the loss of structural and functional integrity of the encoded protein. In humans, members of the minichromosome maintenance (MCM) family play a vital role in cell proliferation and have a significant impact on tumorigenesis. Among the MCM members, the molecular mechanism of how missense SNPs of minichromosome maintenance complex component 6 (MCM6) contribute to DNA replication and tumor pathogenesis is underexplored and needs to be elucidated. Hence, a series of sequence and structure-based computational tools were utilized to determine how mutations affect the corresponding MCM6 protein. From the dbSNP database, among 15,009 SNPs in the MCM6 gene, 642 missense SNPs (4.28%), 291 synonymous SNPs (1.94%), and 12,500 intron SNPs (83.28%) were observed. Out of the 642 missense SNPs, 33 were found to be deleterious during the SIFT analysis. Among these, 11 missense SNPs (I123S, R207C, R222C, L449F, V456M, D463G, H556Y, R602H, R633W, R658C, and P815T) were found as deleterious, probably damaging, affective and disease-associated. Then, I123S, R207C, R222C, V456M, D463G, R602H, R633W, and R658C missense SNPs were found to be highly harmful. Six missense SNPs (I123S, R207C, V456M, D463G, R602H, and R633W) had the potential to destabilize the corresponding protein as predicted by DynaMut2. Interestingly, five high-risk mutations (I123S, V456M, D463G, R602H, and R633W) were distributed in two domains (PF00493 and PF14551). During molecular dynamics simulations analysis, consistent fluctuation in RMSD and RMSF values, high Rg and hydrogen bonds in mutant proteins compared to wild-type revealed that these mutations might alter the protein structure and stability of the corresponding protein. Hence, the results from the analyses guide the exploration of the mechanism by which these missense SNPs of the MCM6 gene alter the structural integrity and functional properties of the protein, which could guide the identification of ways to minimize the harmful effects of these mutations in humans.

Rare genetic coding variants associated with age-related episodic memory decline implicate distinct memory pathologies in the hippocampus

Background

Approximately 40% of people aged 65 or older experience memory loss, particularly in episodic memory. Identifying the genetic basis of episodic memory decline is crucial for uncovering its underlying causes.

Methods

We investigated common and rare genetic variants associated with episodic memory decline in 742 (632 for rare variants) Ashkenazi Jewish individuals (mean age 75) from the LonGenity study. All-atom MD simulations were performed to uncover mechanistic insights underlying rare variants associated with episodic memory decline.

Results

In addition to the common polygenic risk of Alzheimer's Disease (AD), we identified and replicated rare variant association in ITSN1 and CRHR2 . Structural analyses revealed distinct memory pathologies mediated by interfacial rare coding variants such as impaired receptor activation of corticotropin releasing hormone and dysregulated L-serine synthesis.

Discussion

Our study uncovers novel risk loci for episodic memory decline. The identified underlying mechanisms point toward heterogeneous memory pathologies mediated by rare coding variants.

HIV-1 envelope facilitates the development of protease inhibitor resistance through acquiring mutations associated with viral entry and immune escape

Introduction

There is increasing evidence supporting a role for HIV-1 envelope in the development of Protease Inhibitor drug resistance, and a recent report from our group suggested that Env mutations co-evolve with Gag-Protease mutations in the pathway to Lopinavir resistance. In this study, we investigated the effect of co-evolving Env mutations on virus function and structure.

Methods

Co-receptor usage and n-linked glycosylation were investigated using Geno2Pheno as well as tools available at the Los Alamos sequence database. Molecular dynamics simulations were performed using Amber 18 and analyzed using Cpptraj, and molecular interactions were calculated using the Ring server.

Results

The results showed that under Protease Inhibitor drug selection pressure, the envelope gene modulates viral entry by protecting the virus from antibody recognition through the increased length and number of N-glycosylation sites observed in V1/V2 and to some extent V5. Furthermore, gp120 mutations appear to modulate viral entry through a switch to the CXCR4 coreceptor, induced by higher charge in the V3 region and specific mutations at the coreceptor binding sites. In gp41, S534A formed a hydrogen bond with L602 found in the disulfide loop region between the Heptad Repeat 1 and Heptad Repeat 2 domains and could negatively affect the association of gp120-gp41 during viral entry. Lastly, P724Q/S formed both intermolecular and intramolecular interactions with residues within the Kennedy loop, a known epitope.

Discussion

In conclusion, the results suggest that mutations in envelope during Protease Inhibitor treatment failure are related to immune escape and that S534A mutants could preferentially use the cell-to-cell route of infection.

Variant Classification Discordance: Contributing Factors and Predictive Models

An ever-growing catalog of human variants is hosted in the ClinVar database. In this database, submissions on a variant are combined into a multisubmitter record; and in the case of discordance in variant classification between submitters, the record is labeled as conflicting. The current study used ClinVar data to identify characteristics that would make variants more likely to be associated with the conflict class of variants. Furthermore, the Extreme Gradient Boosting algorithm was used to train classifier models to provide prediction of classification discordance for single submission variants in ClinVar database. Population allele frequency, the gene harboring the variant, variant type, consequence on protein, variant deleteriousness score, first submitter identity, and submission count were associated with conflict in variant classification. Using such features, the optimized classifier showed accuracy on the test set of 88% with the weighted average of precision, recall, and f1-score of 0.84, 0.88, and 0.85, respectively. There were pronounced associations between variant classification discordance and allele frequency, gene type, and the identity of the first submitter. The study provides the predicted discordance status for single-submitter variants deposited in ClinVar. This approach can be used to assess whether single-submitter variants are likely to be supported, or in conflict with, future entries; this knowledge may help laboratories with clinical variant assessment.

Biological importance of arginine: A comprehensive review of the roles in structure, disorder, and functionality of peptides and proteins

Arginine shows Jekyll and Hyde behavior in several respects. It participates in protein folding via ionic and H-bonds and cation-pi interactions; the charge and hydrophobicity of its side chain make it a disorder-promoting amino acid. Its methylation in histones; RNA binding proteins; chaperones regulates several cellular processes. The arginine-centric modifications are important in oncogenesis and as biomarkers in several cardiovascular diseases. The cross-links involving arginine in collagen and cornea are involved in pathogenesis of tissues but have also been useful in tissue engineering and wound-dressing materials. Arginine is a part of active site of several enzymes such as GTPases, peroxidases, and sulfotransferases. Its metabolic importance is obvious as it is involved in production of urea, NO, ornithine and citrulline. It can form unusual functional structures such as molecular tweezers in vitro and sprockets which engage DNA chains as part of histones in vivo. It has been used in design of cell-penetrating peptides as drugs. Arginine has been used as an excipient in both solid and injectable drug formulations; its role in suppressing opalescence due to liquid-liquid phase separation is particularly very promising. It has been known as a suppressor of protein aggregation during protein refolding. It has proved its usefulness in protein bioseparation processes like ion-exchange, hydrophobic and affinity chromatographies. Arginine is an amino acid, whose importance in biological sciences and biotechnology continues to grow in diverse ways.

DIS3 Variants are Associated With Primary Ovarian Insufficiency: Importance of Transcription/Translation in Oogenesis

CONTEXT

A genetic etiology accounts for the majority of unexplained primary ovarian insufficiency (POI).

OBJECTIVE

We hypothesized a genetic cause of POI for a sister pair with primary amenorrhea.

DESIGN

The study was an observational study. Subjects were recruited at an academic institution.

SUBJECTS

Subjects were sisters with primary amenorrhea caused by POI and their parents. Additional subjects included women with POI analyzed previously (n = 291). Controls were recruited for health in old age or were from the 1000 Genomes Project (total n = 233).

INTERVENTION

We performed whole exome sequencing, and data were analyzed using the Pedigree Variant Annotation, Analysis and Search Tool, which identifies genes harboring pathogenic variants in families. We performed functional studies in a Drosophila melanogaster model.

MAIN OUTCOME

Genes with rare pathogenic variants were identified.

RESULTS

The sisters carried compound heterozygous variants in DIS3. The sisters did not carry additional rare variants that were absent in publicly available datasets. DIS3 knockdown in the ovary of D. melanogaster resulted in lack of oocyte production and severe infertility.

CONCLUSIONS

Compound heterozygous variants in highly conserved amino acids in DIS3 and failure of oocyte production in a functional model suggest that mutations in DIS3 cause POI. DIS3 is a 3' to 5' exoribonuclease that is the catalytic subunit of the exosome involved in RNA degradation and metabolism in the nucleus. The findings provide further evidence that mutations in genes important for transcription and translation are associated with POI.

Development of a novel prediction model based on protein structure for identifying RPE65-associated inherited retinal disease (IRDs) of missense variants

Purpose

This study aimed to develop a prediction model to classify RPE65-mediated inherited retinal disease (IRDs) based on protein secondary structure and to analyze phenotype-protein structure correlations of RPE65 missense variants in a Chinese cohort.

Methods

Pathogenic or likely pathogenic missense variants of RPE65 were obtained from UniProt, ClinVar, and HGMD databases. The three-dimensional structure of RPE65 was retrieved from the Protein Data Bank (PDB) and modified with Pymol software. A novel prediction model was developed using LASSO regression and multivariate logistic regression to identify RPE65-associated IRDs. A total of 21 Chinese probands with RPE65 variants were collected to analyze phenotype-protein structure correlations of RPE65 missense variants.

Results

The study found that both pathogenic and population missense variants were associated with structural features of RPE65. Pathogenic variants were linked to sheet, β-sheet, strands, β-hairpins, Fe2+ (iron center), and active site cavity, while population variants were related to helix, loop, helices, and helix-helix interactions. The novel prediction model showed accuracy and confidence in predicting the disease type of RPE65 variants (AUC = 0.7531). The study identified 25 missense variants in Chinese patients, accounting for 72.4% of total mutations. A significant correlation was observed between clinical characteristics of RPE65-associated IRDs and changes in amino acid type, specifically for missense variants of F8 (H68Y, P419S).

Conclusion

The study developed a novel prediction model based on the protein structure of RPE65 and investigated phenotype-protein structure correlations of RPE65 missense variants in a Chinese cohort. The findings provide insights into the precise diagnosis of RPE65-mutated IRDs.

Structural analysis and conformational dynamics of SOCS1 gene mutations involved in diffuse large B-cell lymphoma

OBJECTIVES

The SOCS1 gene is frequently mutated in primary Diffuse Large B-Cell Lymphoma (DLBCL) patients and is associated with a reduced survival rate. Using various computational techniques, the current study aims to identify Single Nucleotide Polymorphisms (SNPs) in the SOCS1 gene that are associated with the mortality rate of DLBCL patients. This study also evaluates the effects of SNPs on the structural instability of the SOCS1 protein in DLBCL patient.

METHODS

The cBioPortal webserver was used for mutations and determining how the SNP mutations affect the SOCS1 protein with various algorithms (PolyPhen-2.0, Provean, PhD-SNPg, SNPs&GO, SIFT, FATHMM, Predict SNP and SNAP). Five webservers (I-Mutant 2.0, MUpro, mCSM, DUET and SDM) were used for protein instability and the conserved status and were also predicted through different tools (ConSurf, Expasy, SOMPA). Lastly, MD simulations were run on the two chosen mutations (S116N and V128G) using GROMACS 5.0.1 to study how the mutations change the structure of SOCS1.

RESULTS

Among the 93 SOCS1 mutations detected in DLBCL patients, nine mutations were found to have a detrimental effect (damaging/deleterious/pathogenic/altered) on the SOCS1 protein. All the nine selected mutations are in the conserved region and four are on the extended strand site, four on the random coil site and one on the alpha helix position of the secondary protein structure. After anticipating the structural effects of these nine mutations, two were chosen (S116N and V128G) based on mutational frequency, location within the protein, structural effect (primary, secondary and tertiary) on stability and conservation status within the SOCS1 protein. The simulation of a 50 ns time interval revealed that the Rg value of S116N (2.17 nm) is higher than that of WT (1.98 nm), indicating a loss of structural compactness. In the case of the RMSD value, this mutated type (V128G) shows more deviation (1.54 nm) in comparison to the wild-type (2.14 nm) and another mutant type (S116N) (2.12 nm). The average RMSF values of wild-type and mutant types (V128G and S116N) were 0.88 nm, 0.49 nm, and 0.93 nm, respectively. The RMSF result shows that the mutant V128G structure is more stable than the wild-type and mutant S116N structures.

CONCLUSION

Based on all these computational predictions, this study finds that certain mutations, particularly S116N, have a destabilising and robust effect on the SOCS1 protein. These results can be used to learn more about the importance of SOCS1 mutations in DLBCL patients and to develop new ways to treat DLBCL.

A spinal muscular atrophy modifier implicates the SMN protein in SNARE complex assembly at neuromuscular synapses

Reduced survival motor neuron (SMN) protein triggers the motor neuron disease, spinal muscular atrophy (SMA). Restoring SMN prevents disease, but it is not known how neuromuscular function is preserved. We used model mice to map and identify an Hspa8G470R synaptic chaperone variant, which suppressed SMA. Expression of the variant in the severely affected mutant mice increased lifespan >10-fold, improved motor performance, and mitigated neuromuscular pathology. Mechanistically, Hspa8G470R altered SMN2 splicing and simultaneously stimulated formation of a tripartite chaperone complex, critical for synaptic homeostasis, by augmenting its interaction with other complex members. Concomitantly, synaptic vesicular SNARE complex formation, which relies on chaperone activity for sustained neuromuscular synaptic transmission, was found perturbed in SMA mice and patient-derived motor neurons and was restored in modified mutants. Identification of the Hspa8G470R SMA modifier implicates SMN in SNARE complex assembly and casts new light on how deficiency of the ubiquitous protein causes motor neuron disease.

Missense MED12 variants in 22 males with intellectual disability: From nonspecific symptoms to complete syndromes

We describe the phenotype of 22 male patients (20 probands) carrying a hemizygous missense variant in MED12. The phenotypic spectrum is very broad ranging from nonspecific intellectual disability (ID) to the three well-known syndromes: Opitz-Kaveggia syndrome, Lujan-Fryns syndrome, or Ohdo syndrome. The identified variants were randomly distributed throughout the gene (p = 0.993, χ² test), but mostly outside the functional domains (p = 0.004; χ² test). Statistical analyses did not show a correlation between the MED12-related phenotypes and the locations of the variants (p = 0.295; Pearson correlation), nor the protein domain involved (p = 0.422; Pearson correlation). In conclusion, establishing a genotype-phenotype correlation in MED12-related diseases remains challenging. Therefore, we think that patients with a causative MED12 variant are currently underdiagnosed due to the broad patients' clinical presentations.

Variant predictions in congenital adrenal hyperplasia caused by mutations in CYP21A2

CYP21A2 deficiency represents 95% of congenital adrenal hyperplasia (CAH) cases, a group of genetic disorders that affect steroid biosynthesis. The genetic and functional analysis provide critical tools to elucidate complex CAH cases. One of the most accessible tools to infer the pathogenicity of new variants is in silico prediction. Here, we analyzed the performance of in silico prediction tools to categorize missense single nucleotide variants (SNVs) of CYP21A2. SNVs of CYP21A2 characterized in vitro by functional assays were selected to assess the performance of online single and meta predictors. SNVs were tested separately or in combination with the related phenotype (severe or mild CAH form). In total, 103 SNVs of CYP21A2 (90 pathogenic and 13 neutral) were used to test the performance of 13 single-predictors and four meta-predictors. All SNVs associated with the severe phenotypes were well categorized by all tools, with an accuracy of between 0.69 (PredictSNP2) and 0.97 (CADD), and Matthews’ correlation coefficient (MCC) between 0.49 (PoredicSNP2) and 0.90 (CADD). However, SNVs related to the mild phenotype had more variation, with the accuracy between 0.47 (S3Ds&amp;GO and MAPP) and 0.88 (CADD), and MCC between 0.18 (MAPP) and 0.71 (CADD). From our analysis, we identified four predictors of CYP21A2 variant pathogenicity with good performance, CADD, ConSurf, DANN, and PolyPhen2. These results can be used for future analysis to infer the impact of uncharacterized SNVs in CYP21A2.

Missense variants in ANKRD11 cause KBG syndrome by impairment of stability or transcriptional activity of the encoded protein

PURPOSE

Although haploinsufficiency of ANKRD11 is among the most common genetic causes of neurodevelopmental disorders, the role of rare ANKRD11 missense variation remains unclear. We characterized clinical, molecular, and functional spectra of ANKRD11 missense variants.

METHODS

We collected clinical information of individuals with ANKRD11 missense variants and evaluated phenotypic fit to KBG syndrome. We assessed pathogenicity of variants through in silico analyses and cell-based experiments.

RESULTS

We identified 20 unique, mostly de novo, ANKRD11 missense variants in 29 individuals, presenting with syndromic neurodevelopmental disorders similar to KBG syndrome caused by ANKRD11 protein truncating variants or 16q24.3 microdeletions. Missense variants significantly clustered in repression domain 2 at the ANKRD11 C-terminus. Of the 10 functionally studied missense variants, 6 reduced ANKRD11 stability. One variant caused decreased proteasome degradation and loss of ANKRD11 transcriptional activity.

CONCLUSION

Our study indicates that pathogenic heterozygous ANKRD11 missense variants cause the clinically recognizable KBG syndrome. Disrupted transrepression capacity and reduced protein stability each independently lead to ANKRD11 loss-of-function, consistent with haploinsufficiency. This highlights the diagnostic relevance of ANKRD11 missense variants, but also poses diagnostic challenges because the KBG-associated phenotype may be mild and inherited pathogenic ANKRD11 (missense) variants are increasingly observed, warranting stringent variant classification and careful phenotyping.

In-Silico Characterization of Estrogen Reactivating β-Glucuronidase Enzyme in GIT Associated Microbiota of Normal Human and Breast Cancer Patients

Estrogen circulating in blood has been proved to be a strong biomarker for breast cancer. A β-glucuronidase enzyme (GUS) from human gastrointestinal tract (GIT) microbiota including probiotics has significant involvement in enhancing the estrogen concentration in blood through deconjugation of glucuronidated estrogens. The present project has been designed to explore GIT microbiome-encoded GUS enzymes (GUSOME) repertoire in normal human and breast cancer patients. For this purpose, a total of nineteen GUS enzymes from human GIT microbes, i.e., seven from healthy and twelve from breast cancer patients have been focused on. Protein sequences of enzymes retrieved from UniProt database were subjected to ProtParam, CELLO2GO, SOPMA (secondary structure prediction method), PDBsum (Protein Database summaries), PHYRE2 (Protein Homology/AnalogY Recognition Engine), SAVES v6.0 (Structure Validation Server), MEME version 5.4.1 (Multiple Em for Motif Elicitation), Caver Web server v 1.1, Interproscan and Predicted Antigenic Peptides tool. Analysis revealed the number of amino acids, isoelectric point, extinction coefficient, instability index and aliphatic index of GUS enzymes in the range of 586−795, 4.91−8.92, 89,980−155,075, 25.88−40.93 and 71.01−88.10, respectively. Sub-cellular localization of enzyme was restricted to cytoplasm and inner-membrane in case of breast cancer patients’ bacteria as compared to periplasmic space, outer membrane and extracellular space in normal GIT bacteria. The 2-D structure analysis showed α helix, extended strand, β turn and random coil in the range of 27.42−22.66%, 22.04−25.91%, 5.39−8.30% and 41.75−47.70%, respectively. The druggability score was found to be 0.05−0.45 and 0.06−0.80 in normal and breast cancer patients GIT, respectively. The radius, length and curvature of catalytic sites were observed to be 1.1−2.8 Å, 1.4−15.9 Å and 0.65−1.4, respectively. Ten conserved protein motifs with p < 0.05 and width 25−50 were found. Antigenic propensity-associated sequences were 20−29. Present study findings hint about the use of the bacterial GUS enzymes against breast cancer tumors after modifications via site-directed mutagenesis of catalytic sites involved in the activation of estrogens and through destabilization of these enzymes.

PHEXL222P Mutation Increases Phex Expression in a New ENU Mouse Model for XLH Disease

PhexL222P mouse is a new ENU mouse model for XLH disease due to Leu to Pro amino acid modification at position 222. PhexL222P mouse is characterized by growth retardation, hypophosphatemia, hypocalcemia, reduced body bone length, and increased epiphyseal growth plate thickness and femur diameter despite the increase in PHEXL222P expression. Actually, PhexL222P mice show an increase in Fgf23, Dmp1, and Mepe and Slc34a1 (Na-Pi IIa cotransporter) mRNA expression similar to those observed in Hyp mice. Femoral osteocalcin and sclerostin and Slc34a1 do not show any significant variation in PhexL222P mice. Molecular dynamics simulations support the experimental data. P222 might locally break the E217-Q224 β-sheet, which in turn might disrupt inter-β-sheet interactions. We can thus expect local protein misfolding, which might be responsible for the experimentally observed PHEXL222P loss of function. This model could be a valuable addition to the existing XLH model for further comprehension of the disease occurrence and testing of new therapies.

Insight into the Conformational Transitions of Serine Acetyl Transferase Isoforms in E. histolytica: Implications for Structural and Functional Balance

Serine acetyl transferase (SAT) is one of the crucial enzymes in the cysteine biosynthetic pathway and an essential enzyme for the survival of Entamoeba histolytica, the causative agent of amoebiasis. E. histolytica expresses three isoforms of SAT, where SAT1 and SAT2 are inhibited by the final product cysteine, while SAT3 is not inhibited. SAT3 has a slightly elongated C-terminus compared to SAT1. To understand the stability and conformational transition between two secondary structures of proteins, we measured the effect of urea, a chemical denaturant, on two isoforms of SAT (SAT1 and SAT3) of E. histolytica. The effect of urea on the structure and stability of SAT1 and SAT3 was determined by measuring changes in their far-UV circular dichroism (CD), Trp fluorescence, and near-UV absorption spectra. The urea-induced normal transition curves suggested that the structural transition is reversible and follows a two-state process. Analysis of the urea-induced transition of all optical properties for the stability parameters ΔG _D° (Gibbs free energy change (ΔG _D) in the absence of urea), m (dependence of ΔG _D on urea concentration), and C _m (midpoint of urea transition) suggested that SAT1 is more stable than SAT3. Characterization of the end product of the urea-induced transition of both proteins by the far-UV CD and Trp-fluorescence and near-UV absorbance suggested that urea causes α-helix to β-sheet transition and burial of Trp residues, respectively. To support the in vitro findings, 100 ns molecular dynamics simulations (in silico study) were performed. Both the spectroscopic and molecular dynamics approaches clearly indicated that SAT1 is more stable than SAT3. SAT3 has evolved to escape the feedback inhibition to keep producing cysteine, but in the process, it compromises its structural stability relative to SAT1.

Toxicological Implications of Biological Heterogeneity

From a micro to macro scale of biological organization, macromolecular diversity and biological heterogeneity are fundamental properties of biological systems. Heterogeneity may result from genetic, epigenetic, and non-genetic characteristics (e.g., tissue microenvironment). Macromolecular diversity and biological heterogeneity are tolerated as long as the sustenance and propagation of life are not disrupted. They also provide the raw materials for microevolutionary changes that may help organisms adapt to new selection pressures arising from the environment. Sequence evolution, functional divergence, and positive selection of gene and promoter dosage play a major role in the evolution of life's diversity including complex metabolic networks, which is ultimately reflected in changes in the allele frequency over time. Robustness in evolvable biological systems is conferred by functional redundancy that is often created by macromolecular diversity and biological heterogeneity. The ability to investigate biological macromolecules at an increasingly finer level has uncovered a wealth of information in this regard. Therefore, the dynamics of biological complexity should be taken into consideration in biomedical research.

VARIDT 2.0: structural variability of drug transporter

The structural variability data of drug transporter (DT) are key for research on precision medicine and rational drug use. However, these valuable data are not sufficiently covered by the available databases. In this study, a major update of VARIDT (a database previously constructed to provide DTs' variability data) was thus described. First, the experimentally resolved structures of all DTs reported in the original VARIDT were discovered from PubMed and Protein Data Bank. Second, the structural variability data of each DT were collected by literature review, which included: (a) mutation-induced spatial variations in folded state, (b) difference among DT structures of human and model organisms, (c) outward/inward-facing DT conformations and (d) xenobiotics-driven alterations in the 3D complexes. Third, for those DTs without experimentally resolved structural variabilities, homology modeling was further applied as well-established protocol to enrich such valuable data. As a result, 145 mutation-induced spatial variations of 42 DTs, 1622 inter-species structures originating from 292 DTs, 118 outward/inward-facing conformations belonging to 59 DTs, and 822 xenobiotics-regulated structures in complex with 57 DTs were updated to VARIDT (https://idrblab.org/varidt/ and http://varidt.idrblab.net/). All in all, the newly collected structural variabilities will be indispensable for explaining drug sensitivity/selectivity, bridging preclinical research with clinical trial, revealing the mechanism underlying drug-drug interaction, and so on.

Investigation of Interferon Gamma Activity Using Bioinformatics Methods

Breast cancer grows from the breast tissue and is a severe health problem worldwide. Genetics is believed to be the primary cause of all cases of breast cancer via gene mutation. Bioinformatics methodology has been used to determine the sequences and structures of bioactive substances. This study aimed to analyze the function and structure of the Interferon Gamma (IFNγ) in healthy controls and patients with breast cancer using bioinformatics methods. Blood samples were collected from 75 patients with breast cancer and 25 healthy subjects as control samples. The results showed transition mutation (30%) and transversion mutation (70%) in patients with breast cancer. Moreover, missense mutations (84%) and silent mutations (16%) were detected by BLAST. In addition, the amino acid of the IFNγ protein consisting of alpha-helical, β-sheet, and coil of secondary structure was determined in this study using BioEdit. The results of the physicochemical properties of the IFNγ protein reflect the function, stability, molecular weight, isoelectric point, and instability index of the IFNγ protein using ProtParam. Moreover, the results of mutation affected the percentage of alpha-helix, β-turns, and coil in breast cancer patients compared to healthy groups with reference of NCBI using PSIpred program. Additionally, the PHYRE2 server and RasMol program showed a tertiary structure of the IFNγ protein in breast cancer patients. Furthermore, the STRING program revealed the poly IFNγ protein interacted with other proteins to perform its functions normally. From the recorded data in the current study, it was concluded that IFNγ is considered a marker for patients with breast cancer.

Mutational sources of trans-regulatory variation affecting gene expression in Saccharomyces cerevisiae

Heritable variation in a gene’s expression arises from mutations impacting cis- and trans-acting components of its regulatory network. Here, we investigate how trans-regulatory mutations are distributed within the genome and within a gene regulatory network by identifying and characterizing 69 mutations with trans-regulatory effects on expression of the same focal gene in Saccharomyces cerevisiae. Relative to 1766 mutations without effects on expression of this focal gene, we found that these trans-regulatory mutations were enriched in coding sequences of transcription factors previously predicted to regulate expression of the focal gene. However, over 90% of the trans-regulatory mutations identified mapped to other types of genes involved in diverse biological processes including chromatin state, metabolism, and signal transduction. These data show how genetic changes in diverse types of genes can impact a gene’s expression in trans, revealing properties of trans-regulatory mutations that provide the raw material for trans-regulatory variation segregating within natural populations.

Clinical, Biochemical, and Genetic Heterogeneity in Glutaric Aciduria Type II Patients

The variants of electron transfer flavoprotein (ETFA, ETFB) and ETF dehydrogenase (ETFDH) are the leading cause of glutaric aciduria type II (GA-II). In this study, we identified 13 patients harboring six variants of two genes associated with GA-II. Out of the six variants, four were missense, and two were frameshift mutations. A missense variant (ETFDH:p.Gln269His) was observed in a homozygous state in nine patients. Among nine patients, three had experienced metabolic crises with recurrent vomiting, abdominal pain, and nausea. In one patient with persistent metabolic acidosis, hypoglycemia, and a high anion gap, the ETFDH:p.Gly472Arg, and ETFB:p.Pro94Thrfs*8 variants were identified in a homozygous, and heterozygous state, respectively. A missense variant ETFDH:p.Ser442Leu was detected in a homozygous state in one patient with metabolic acidosis, hypoglycemia, hyperammonemia and liver dysfunction. The ETFDH:p.Arg41Leu, and ETFB:p.Ile346Phefs*19 variants were observed in a homozygous state in one patient each. Both these variants have not been reported so far. In silico approaches were used to evaluate the pathogenicity and structural changes linked with these six variants. Overall, the results indicate the importance of a newborn screening program and genetic investigations for patients with GA-II. Moreover, careful interpretation and correlation of variants of uncertain significance with clinical and biochemical findings are needed to confirm the pathogenicity of such variants.

Comprehensive Evaluation of PAXgene Fixation on Oral Cancer Tissues Using Routine Histology, Immunohistochemistry, and FTIR Microspectroscopy

The choice of tissue fixation is critical for preserving the morphology and biochemical information of tissues. Fragile oral tissues with lower tensile strength are challenging to process for histological applications as they are prone to processing damage, such as tissue tear, wrinkling, and tissue fall-off from slides. This leads to loss of morphological information and unnecessary delay in experimentation. In this study, we have characterized the new PAXgene tissue fixation system on oral buccal mucosal tissue of cancerous and normal pathology for routine histological and immunohistochemical applications. We aimed to minimize the processing damage of tissues and improve the quality of histological experiments. We also examined the preservation of biomolecules by PAXgene fixation using FTIR microspectroscopy. Our results demonstrate that the PAXgene-fixed tissues showed significantly less tissue fall-off from slides. Hematoxylin and Eosin staining showed comparable morphology between formalin-fixed and PAXgene-fixed tissues. Good quality and slightly superior immunostaining for cancer-associated proteins p53 and CK5/6 were observed in PAXgene-fixed tissues without antigen retrieval than formalin-fixed tissues. Further, FTIR measurements revealed superior preservation of glycogen, fatty acids, and amide III protein secondary structures in PAXgene-fixed tissues. Overall, we present the first comprehensive evaluation of the PAXgene tissue fixation system in oral tissues. This study concludes that the PAXgene tissue fixation system can be applied to oral tissues to perform diagnostic molecular pathology experiments without compromising the quality of the morphology or biochemistry of biomolecules.

R869C mutation in molecular motor KIF17 gene is involved in dementia with Lewy bodies

INTRODUCTION THE

GBA-N370S mutation is one of the most frequent risk factors for dementia with Lewy bodies (DLB) and Parkinson's disease (PD). We looked for genetic variations that contribute to the outcome in N370S-carriers, whether PD or DLB.

METHODS

Whole-genome sequencing of 95 Ashkenazi-N370S-carriers affected with either DLB (n = 19) or PD (n = 76) was performed, and 564 genes related to dementia and PD analyzed.

RESULTS

We identified enrichment of linked alleles in PINK1 locus in DLB patients (false discovery rate P = .0412). Haplotype analysis delineated 1.8 Mb interval encompassing 29 genes and 87 unique variants, of them, KIF17-R869C received the highest functional prediction score (Combined Annotation Dependent Depletion = 34). Its frequency was significantly higher in 26 DLB-N370S-carriers compared to 140 PD-N370S-carriers (odds ratio [OR] = 33.4 P = .001, and OR = 70.2 when only heterozygotes were included).

DISCUSSION

Because KIF17 was shown to be important for learning and memory in mice, our data further suggest, for the first time, its involvement in DLB, and possibly in human dementia.

Isolated congenital diaphragm hernia associated with homozygous SLIT3 gene variant in dizygous twins

Congenital diaphragmatic hernia (CDH) is a serious life-threatening birth defect characterized by abnormal development in the muscular or tendinous portion of the diaphragm during embryogenesis. Despite its high incidence, the etiology of CDH hasn't been fully understood. Genetic factors are important in pathogenesis; however, few single genes have been definitively implicated in human CDH. SLIT1, SLIT2, and SLIT3 (slit guidance ligand) are three human homologs of the drosophila Slit gene. They interact with roundabout (Robo) homolog receptors to affect cell migration, adhesion, cell motility, and angiogenesis and play important roles in cell signaling pathways including the guidance of axons. In this report, we presented dizygous twin babies with CDH related to the SLIT3 gene variant. Previous studies showed that Slit3 null mice had congenital diaphragmatic hernias on or near the ventral midline portion of the central tendon. This is the first report of homozygous SLIT3 variant associated with CDH in humans.

Computational analysis of Cyclin D1 gene SNPs and association with breast cancer

CCND1 encodes for Cyclin D1 protein and single-nucleotide polymorphisms (SNPs) can modulate its activity. In the present study, the impact of CCND1 SNPs on structure and/or function of Cyclin D1 protein using in silico tools was investigated. Our analysis revealed only one splice site SNP (c.1988+5G<A) can effect CCND1 function. Subsequently, 78 out of 169 missense variants were predicted as pathogenic by Polyphen2, SIFT, PROVEAN, SNPs&GO, and PANTHER, and 4/78 missense SNPs were further evaluated because these four SNPs were found to be reside in highly conserved region of Cyclin D1. However, they did not show any major impact on tertiary structure and domain of Cyclin D1 but overall R15S and A190S has displayed a significant diseased phenotype and an altered molecular mechanism predicted by MutPred, FATHMM, SNPeffect, SNAP2, and PredictSNP. Consistently, A190S, R179L, and R15S may also cause a decrease in stability of Cyclin D1 anticipated by I-Mutant, HOPE and SNP effect. Furthermore, the Kaplan–Meier plotter has explained that high expression of CCND1 is associated with less survival rate of breast cancer patients. Altogether our study suggests that c.1988+5G<A, R15S, R179L, and A190S SNPs could directly or indirectly destabilize Cyclin D1.

Weakening of interaction networks with aging in tip-link protein induces hearing loss

Age-related hearing loss (ARHL) is a common condition in humans marking the gradual decrease in hearing with age. Perturbations in the tip-link protein cadherin-23 that absorbs the mechanical tension from sound and maintains the integrity of hearing is associated with ARHL. Here, in search of molecular origins for ARHL, we dissect the conformational behavior of cadherin-23 along with the mutant S47P that progresses the hearing loss drastically. Using an array of experimental and computational approaches, we highlight a lower thermodynamic stability, significant weakening in the hydrogen-bond network and inter-residue correlations among β-strands, due to the S47P mutation. The loss in correlated motions translates to not only a remarkable two orders of magnitude slower folding in the mutant but also to a proportionately complex unfolding mechanism. We thus propose that loss in correlated motions within cadherin-23 with aging may trigger ARHL, a molecular feature that likely holds true for other disease-mutations in β-strand-rich proteins.

Inferring the molecular and phenotypic impact of amino acid variants with MutPred2

Identifying pathogenic variants and underlying functional alterations is challenging. To this end, we introduce MutPred2, a tool that improves the prioritization of pathogenic amino acid substitutions over existing methods, generates molecular mechanisms potentially causative of disease, and returns interpretable pathogenicity score distributions on individual genomes. Whilst its prioritization performance is state-of-the-art, a distinguishing feature of MutPred2 is the probabilistic modeling of variant impact on specific aspects of protein structure and function that can serve to guide experimental studies of phenotype-altering variants. We demonstrate the utility of MutPred2 in the identification of the structural and functional mutational signatures relevant to Mendelian disorders and the prioritization of de novo mutations associated with complex neurodevelopmental disorders. We then experimentally validate the functional impact of several variants identified in patients with such disorders. We argue that mechanism-driven studies of human inherited disease have the potential to significantly accelerate the discovery of clinically actionable variants.

A genome-wide scan study identifies a single nucleotide substitution in MC1R gene associated with white coat colour in fallow deer (Dama dama)

Background

The coat colour of fallow deer is highly variable and even white animals can regularly be observed in game farming and in the wild. Affected animals do not show complete albinism but rather some residual pigmentation resembling a very pale beige dilution of coat colour. The eyes and claws of the animals are pigmented. To facilitate the conservation and management of such animals, it would be helpful to know the responsible gene and causative variant. We collected 102 samples from 22 white animals and from 80 animals with wildtype coat colour. The samples came from 12 different wild flocks or game conservations located in different regions of Germany, at the border to Luxembourg and in Poland. The genomes of one white hind and her brown calf were sequenced.

Results

Based on a list of colour genes of the International Federation of Pigment Cell Societies (http://www.ifpcs.org/albinism/), a variant in the MC1R gene (NM_174108.2:c.143 T > C) resulting in an amino acid exchange from leucine to proline at position 48 of the MC1R receptor protein (NP_776533.1:p.L48P) was identified as a likely cause of coat colour dilution. A gene test revealed that all animals of the white phenotype were of genotype CC whereas all pigmented animals were of genotype TT or TC. The study showed that 14% of the pigmented (brown or dark pigmented) animals carried the white allele.

Conclusions

A genome-wide scan study led to a molecular test to determine the coat colour of fallow deer. Identification of the MC1R gene provides a deeper insight into the mechanism of dilution. The gene marker is now available for the conservation of white fallow deer in wild and farmed animals.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-020-00950-3.

Pyrazinamide resistance of novel mutations in pncA and their dynamic behavior

Pyrazinamide (PZA) is one of the essential anti-mycobacterium drugs, active against non-replicating Mycobacterium tuberculosis (MTB) isolates. PZA is converted into its active state, called pyrazinoic acid (POA), by action of pncA encoding pyrazinamidase (PZase). In the majority of PZA-resistance isolates, pncA harbored mutations in the coding region. In our recent report, we detected a number of novel variants in PZA-resistance (PZA^R) MTB isolates, whose resistance mechanisms were yet to be determined. Here we performed several analyses to unveil the PZA^R mechanism of R123P, T76P, G150A, and H71R mutants (MTs) through molecular dynamics (MD) simulations. In brief, culture positive MTB isolates were subjected to PZA susceptibility tests using the WHO recommended concentration of PZA (100 μg ml⁻¹). The PZA^R samples were screened for mutations in pncA along sensitive isolates through polymerase chain reactions and sequencing. A large number of variants (GeneBank accession no. MH461111), including R123P, T76P, G150A, and H71R, have been spotted in more than 70% of isolates. However, the mechanism of PZA^R for mutants (MTs) R123P, T76P, G150A, and H71R was unknown. For the MTs and native PZase structures (WT), thermodynamic properties were compared using molecular dynamics simulations for 100 ns. The MTs structural activity was compared to the WT. Folding effect and pocket volume variations have been detected when comparing between WT and MTs. Geometric matching further confirmed the effect of R123P, T76P, G150A, and H71R mutations on PZase dynamics, making them vulnerable for activating the pro-drug into POA. This study offers a better understanding for management of PZA^R TB. The results may be used as alternative diagnostic tools to infer PZA resistance at a structural dynamics level.

We performed several analyses to unveil the pyrazinamide-resistance mechanism of R123P, T76P, G150A, and H71R mutants through molecular dynamics simulations.

Molecular characterization of pyridoxine 5'-phosphate oxidase and its pathogenic forms associated with neonatal epileptic encephalopathy

Defects of vitamin B₆ metabolism are responsible for severe neurological disorders, such as pyridoxamine 5'-phosphate oxidase deficiency (PNPOD; OMIM: 610090), an autosomal recessive inborn error of metabolism that usually manifests with neonatal-onset severe seizures and subsequent encephalopathy. At present, 27 pathogenic mutations of the gene encoding human PNPO are known, 13 of which are homozygous missense mutations; however, only 3 of them have been characterised with respect to the molecular and functional properties of the variant enzyme forms. Moreover, studies on wild type and variant human PNPOs have so far largely ignored the regulation properties of this enzyme. Here, we present a detailed characterisation of the inhibition mechanism of PNPO by pyridoxal 5'-phosphate (PLP), the reaction product of the enzyme. Our study reveals that human PNPO has an allosteric PLP binding site that plays a crucial role in the enzyme regulation and therefore in the regulation of vitamin B₆ metabolism in humans. Furthermore, we have produced, recombinantly expressed and characterised several PNPO pathogenic variants responsible for PNPOD (G118R, R141C, R225H, R116Q/R225H, and X262Q). Such replacements mainly affect the catalytic activity of PNPO and binding of the enzyme substrate and FMN cofactor, leaving the allosteric properties unaltered.

Molecular modelling and dynamic simulations of sequestosome 1 (SQSTM1) missense mutations linked to Paget disease of bone

The Paget disease (PDB; OMIM is 167250) is a chronic bone disease caused by pathogenic mutations in Sequestome1/p62 (SQSTM1) gene. This study has aimed to interpret the relationship of PDB linked SQSTM1 mutations with protein structure and its molecular dynamic features. The disease causative missense mutations were initially collected, and then analyzed for their, exonic and domain distribution, impact on secondary and tertiary structures, and their ability on protein-ligand interactions, using a combination of systems biology approaches. Our results show that most PDB linked SQSTM1 missense mutations affect amino acid residues clustered within or near the UBA domain (aa 389-434), which participates in the ubiquitination of substrates. We also report that the majority mutations occurred in α-helices over β-strands but their effects on the secondary structure were mostly neutral. Global tertiary structure deviations were minimal; however, at amino acid residue level minor structural changes were evident. The molecular dynamics simulation analysis showed that both PB1 and UBA domains were under constant structural fluctuations resulting in closed form conformation of SQSMT1 protein structure, when it is bound to PRKCI ligand. We also found salt bridge conformation changes in the UBA domain of SQSTM1 mutants when they bound to the PRKCI interactor protein. This finding suggests the possibility that mutations in SQSTM1 could impair its ability to ubiquitinate the substrates, eventually affecting autophagy and apoptosis, especially in mature osteoclasts. This study presents the additional insight into structure and function relationship between SQSTM1 mutations and PDB pathogenesis. Communicated by Ramaswamy H. Sarma.

A genome-wide scan study identifies a single nucleotide substitution in the tyrosinase gene associated with white coat colour in a red deer (Cervus elaphus) population

Background

Red deer with very pale coat colour are observed sporadically. In the red deer (Cervus elaphus) population of Reinhardswald in Germany, about 5% of animals have a white coat colour that is not associated with albinism. In order to facilitate the conservation of the animals, it should be determined whether and to what extent brown animals carry the white gene. For this purpose, samples of one white hind and her brown calf were available for whole genome sequencing to identify the single nucleotide polymorphism(s) responsible for the white phenotype. Subsequently, samples from 194 brown and 11 white animals were genotyped.

Results

Based on a list of colour genes of the International Federation of Pigment Cell Societies, a non-synonymous mutation with exchange of a glycine residue at position 291 of the tyrosinase protein by arginine was identified as the cause of dilution of the coat colour. A gene test led to exactly matching genotypes in all examined animals. The study showed that 14% of the brown animals carry the white gene. This provides a simple and reliable way of conservation for the white animals. However, results could not be transferred to another, unrelated red deer population with white animals. Although no brown animals with a white tyrosinase genotype were detected, the cause for the white colouring in this population was different.

Conclusions

A gene test for the conservation of white red deer is available for the population of the Reinhardswald. While mutations in the tyrosinase are commonly associated with oculocutaneous albinism type 1, the amino acid exchange at position 291 was found to be associated with coat colour dilution in Cervus elaphus.

Genotyping of the EIF1AY Gene in Iranian Patients with Non-Obstructive Azoospermia

Background

EIF1AY is one of the genes essential for normal spermatogenesis and is located in azoospermic factors region.

Objective

The present study was designed to investigate the EIF1AY gene nucleotide variations, and correlate it with spermatogenic maturation arrest and azoospermia in Iranian population.

Methods

A total number of 30 Iranian idiopathic non-obstructive azoospermic patients were selected as case group and 30 fertile men served as a control group who had at least 1 child. Nucleotide variation was analyzed in exon 3 and exon 5 in EIF1AY gene of both groups. DNA extraction from peripheral blood samples of selected individuals was done followed by amplification by PCR and sequencing with Sangar method.

Results

Totally 3 single nucleotide variations were identified: one in the intronic region of exon 3, next one in non-coding transcript exon variant (rs13447352) and the third one in the exonic region of exon 5, all were registered in NCBI-Gene database.

Conclusion

There was no statistically significant difference in the incidence of nucleotide variation between 2 study populations (p > 0.05). Further studies are required to specify the effects of Y: T20588295G variation on modification of protein structure, as well as the expression pattern of the gene and its association with azoospermia.

Partial sequence analyses of exon 7 of the ABO locus of cynomolgus (Macaca fascicularis) and rhesus (M. mulatta) macaques: Indeterminate phenotypes show the presence of the O blood group

Compatibility tests to identify A, B, and O alleles are critical for establishing suitable donor-recipient matches among experimental animals. Using a qPCR-based SNP probe assay, we have identified A, B, AB, and indeterminate blood group phenotypes in cynomolgus and rhesus macaques. We have hypothesized, albeit without molecular confirmation, that the indeterminate phenotype represents homozygosity for the null O allele at the macaque ABO locus. The indeterminate phenotype represents the unsuccessful detection of either A or B alleles using primers targeting the A-specific and B-specific single nucleotide polymorphisms (SNPs) in a variable region of exon 7 of the ABO locus. These SNPs are associated with two functional sites, detected using two allele-specific probes in the qPCR assay where the codons leucine and methionine (at codon 266) and glycine and alanine (at codon 268) are required for the synthesis of the A and B transferases, respectively. While reference sequences for the A and B alleles exhibited no novel mutations in the functional exon, plasmid Sanger sequence analyses showed unique mutations within the diagnostic target sites in 10 macaques exhibiting the indeterminate phenotype. Eight of these indeterminate individuals exhibited SNPs at codon 268 that should prevent the syntheses of an A or B transferase. While the two other indeterminate samples had functional codons that were consistent with A or B alleles, mutations in either their probe- or primer-binding sites that altered their peptide sequences probably impeded their detection by our assay.

High Functioning Autism with Missense Mutations in Synaptotagmin-Like Protein 4 (SYTL4) and Transmembrane Protein 187 (TMEM187) Genes: SYTL4- Protein Modeling, Protein-Protein Interaction, Expression Profiling and MicroRNA Studies

We describe a 7-year-old male with high functioning autism spectrum disorder (ASD) and maternally-inherited rare missense variant of Synaptotagmin-like protein 4 (SYTL4) gene (Xq22.1; c.835C>T; p.Arg279Cys) and an unknown missense variant of Transmembrane protein 187 (TMEM187) gene (Xq28; c.708G>T; p. Gln236His). Multiple in-silico predictions described in our study indicate a potentially damaging status for both X-linked genes. Analysis of predicted atomic threading models of the mutant and the native SYTL4 proteins suggest a potential structural change induced by the R279C variant which eliminates the stabilizing Arg279-Asp60 salt bridge in the N-terminal half of the SYTL4, affecting the functionality of the protein's critical RAB-Binding Domain. In the European (Non-Finnish) population, the allele frequency for this variant is 0.00042. The SYTL4 gene is known to directly interact with several members of the RAB family of genes, such as, RAB27A, RAB27B, RAB8A, and RAB3A which are known autism spectrum disorder genes. The SYTL4 gene also directly interacts with three known autism genes: STX1A, SNAP25 and STXBP1. Through a literature-based analytical approach, we identified three of five (60%) autism-associated serum microRNAs (miRs) with high predictive power among the total of 298 mouse Sytl4 associated/predicted microRNA interactions. Five of 13 (38%) miRs were differentially expressed in serum from ASD individuals which were predicted to interact with the mouse equivalent Sytl4 gene. TMEM187 gene, like SYTL4, is a protein-coding gene that belongs to a group of genes which host microRNA genes in their introns or exons. The novel Q236H amino acid variant in the TMEM187 in our patient is near the terminal end region of the protein which is represented by multiple sequence alignments and hidden Markov models, preventing comparative structural analysis of the variant harboring region. Like SYTL4, the TMEM187 gene is expressed in the brain and interacts with four known ASD genes, namely, HCFC1; TMLHE; MECP2; and GPHN. TMM187 is in linkage with MECP2, which is a well-known determinant of brain structure and size and is a well-known autism gene. Other members of the TMEM gene family, TMEM132E and TMEM132D genes are associated with bipolar and panic disorders, respectively, while TMEM231 is a known syndromic autism gene. Together, TMEM187 and SYTL4 genes directly interact with recognized important ASD genes, and their mRNAs are found in extracellular vesicles in the nervous system and stimulate target cells to translate into active protein. Our evidence shows that both these genes should be considered as candidate genes for autism. Additional biological testing is warranted to further determine the pathogenicity of these gene variants in the causation of autism.

A novel AVPR2 missense mutation in an Asian family with inherited nephrogenic diabetes insipidus: A case report

RATIONALE

X-linked nephrogenic diabetes insipidus (NDI) is a rare inherited disease, and is characterized by renal resistance to arginine vasopressin (AVP). Its diagnosis can be clinically challenging. The application of molecular genetic analysis can provide a rapid and definitive diagnosis.

PATIENT CONCERNS

A 75-year-old woman presented with recurrent nausea and vomiting was admitted to the Department of Gastroenterology. The patient had a strong family history of polydipsia and polyuria. Sequencing analysis of the antidiuretic hormone arginine vasopressin receptor 2 (AVPR2) revealed the novel missense mutation p. Trp164Cys (c.492G>G/C) in exon 2. There was a heterozygous mutation in the patient's sister and niece, while there was a mutation in her sons, brother and nephews. The locus is located on the X chromosome Xq28, and its mutation can lead to X linked recessive NDI. The p. Trp164Cys mutation of AVPR2 gene has not been reported in literature before. The mutation was predicted to be probably damaging by several prediction methods, including SIFT and PolyPhen-2. There was no significant abnormal variation in other detection regions of the gene. And there was also no abnormal variation in AVP and AQP2 genes in this family.

DIAGNOSIS

X-linked NDI was diagnosed according to the patient's family history and DNA sequencing analysis.

INTERVENTIONS AND OUTCOMES

After treated with desmopressin, antiemetic drugs and massive infusion glucose transfusion, the patient's urine volume decreased and electrolyte disturbance was corrected, and the symptoms of nausea and vomiting gradually disappeared.

LESSONS

The patients with suspected congenital NDI should undergo genetic sequencing analysis of AVPR2, AVP and AQP2 genes. A definitive diagnosis can benefit patient and avoid unnecessary investigations.

Computational Molecular Phenotypic Analysis of PTPN22 (W620R), IL6R (D358A), and TYK2 (P1104A) Gene Mutations of Rheumatoid Arthritis

Rheumatoid arthritis (RA) is a chronic autoimmune disorder of bone joints caused by the complex interplay between several factors like body physiology, the environment with genetic background. The recent meta-analysis of GWAS has expanded the total number of RA-associated loci to more than 100, out of which approximately ∼97% (98 variants) loci are located in non-coding regions, and the other ∼3% (3 variants) are in three different non-HLA genes, i.e., TYK2 (Prp1104Ala), IL6R (Asp358Ala), and PTPN22 (Trp620Arg). However, whether these variants prompt changes in the protein phenotype with regards to its stability, structure, and interaction with other molecules, remains unknown. Thus, we selected the three clinically pathogenic variants described above, as positive controls and applied diverse computational methods to scrutinize if those mutations cause changes in the protein phenotype. Both wild type and mutant protein structures of PTPN22 (W620R), IL6R (D358A), and TYK2 (P1104A) were modeled and studied for structural deviations. Furthermore, we have also studied the secondary structure characteristics, solvent accessibility and stability, and the molecular interaction deformities caused by the amino acid substitutions. We observed that simple nucleotide predictions of SIFT, PolyPhen, CADD and FATHMM yields mixed findings in screening the RA-missense variants which showed a ≥P-value threshold of 5 × 10^-8 in genome wide association studies. However, structure-based analysis confirms that mutant structures shows subtle but significant changes at their core regions, but their functional domains seems to lose wild type like functional interaction. Our findings suggest that the multidirectional computational analysis of clinically potential RA-mutations could act as a primary screening step before undertaking functional biology assays.

Sequence variant analysis of RNA sequences in severe equine asthma

Background

Severe equine asthma is a chronic inflammatory disease of the lung in horses similar to low-Th2 late-onset asthma in humans. This study aimed to determine the utility of RNA-Seq to call gene sequence variants, and to identify sequence variants of potential relevance to the pathogenesis of asthma.

Methods

RNA-Seq data were generated from endobronchial biopsies collected from six asthmatic and seven non-asthmatic horses before and after challenge (26 samples total). Sequences were aligned to the equine genome with Spliced Transcripts Alignment to Reference software. Read preparation for sequence variant calling was performed with Picard tools and Genome Analysis Toolkit (GATK). Sequence variants were called and filtered using GATK and Ensembl Variant Effect Predictor (VEP) tools, and two RNA-Seq predicted sequence variants were investigated with both PCR and Sanger sequencing. Supplementary analysis of novel sequence variant selection with VEP was based on a score of <0.01 predicted with Sorting Intolerant from Tolerant software, missense nature, location within the protein coding sequence and presence in all asthmatic individuals. For select variants, effect on protein function was assessed with Polymorphism Phenotyping 2 and screening for non-acceptable polymorphism 2 software. Sequences were aligned and 3D protein structures predicted with Geneious software. Difference in allele frequency between the groups was assessed using a Pearson’s Chi-squared test with Yates’ continuity correction, and difference in genotype frequency was calculated using the Fisher’s exact test for count data.

Results

RNA-Seq variant calling and filtering correctly identified substitution variants in PACRG and RTTN. Sanger sequencing confirmed that the PACRG substitution was appropriately identified in all 26 samples while the RTTN substitution was identified correctly in 24 of 26 samples. These variants of uncertain significance had substitutions that were predicted to result in loss of function and to be non-neutral. Amino acid substitutions projected no change of hydrophobicity and isoelectric point in PACRG, and a change in both for RTTN. For PACRG, no difference in allele frequency between the two groups was detected but a higher proportion of asthmatic horses had the altered RTTN allele compared to non-asthmatic animals.

Discussion

RNA-Seq was sensitive and specific for calling gene sequence variants in this disease model. Even moderate coverage (<10–20 counts per million) yielded correct identification in 92% of samples, suggesting RNA-Seq may be suitable to detect sequence variants in low coverage samples. The impact of amino acid alterations in PACRG and RTTN proteins, and possible association of the sequence variants with asthma, is of uncertain significance, but their role in ciliary function may be of future interest.

NDDVD: an integrated and manually curated Neurodegenerative Diseases Variation Database

Neurodegenerative diseases (NDDs) are associated with genetic variations including point substitutions, copy number alterations, insertions and deletions. At present, a few genetic variation repositories for some individual NDDs have been created, however, these databases are needed to be integrated and expanded to all the NDDs for systems biological investigation. We here build a relational database termed as NDDVD to integrate all the variations of NDDs using Leiden Open Variation Database (LOVD) platform. The items in the NDDVD are collected manually from PubMed or extracted from the existed variation databases. The cross-disease database includes over 6374 genetic variations of 289 genes associated with 37 different NDDs. The patterns, conservations and biological functions for variations in different NDDs are statistically compared and a user-friendly interface is provided for NDDVD at: http://bioinf.suda.edu.cn/NDDvarbase/LOVDv.3.0.

URL: http://bioinf.suda.edu.cn/NDDvarbase/LOVDv.3.0

Uncovering the molecular mechanisms behind disease-associated leptin variants

The pleiotropic hormone leptin has a pivotal role in regulating energy balance by inhibiting hunger and increasing energy expenditure. Homozygous mutations found in the leptin gene are associated with extreme obesity, marked hyperphagia, and impaired immune function. Although these mutations have been characterized in vivo, a detailed understanding of how they affect leptin structure and function remains elusive. In the current work, we used NMR, differential scanning calorimetry, molecular dynamics simulations, and bioinformatics calculations to characterize the effects of these mutations on leptin structure and function and binding to its cognate receptor. We found that mutations identified in patients with congenital leptin deficiency not only cause leptin misfolding or aggregation, but also cause changes in the dynamics of leptin residues on the receptor-binding interface. Therefore, we infer that mutation-induced leptin deficiency may arise from several distinct mechanisms including (i) blockade of leptin receptor interface II, (ii) decreased affinity in the second step of leptin's interaction with its receptor, (iii) leptin destabilization, and (iv) unsuccessful threading through the covalent loop, leading to leptin misfolding/aggregation. We propose that this expanded framework for understanding the mechanisms underlying leptin deficiency arising from genetic mutations may be useful in designing therapeutics for leptin-associated disorders.

Sixty-five years of the long march in protein secondary structure prediction: the final stretch?

Protein secondary structure prediction began in 1951 when Pauling and Corey predicted helical and sheet conformations for protein polypeptide backbone even before the first protein structure was determined. Sixty-five years later, powerful new methods breathe new life into this field. The highest three-state accuracy without relying on structure templates is now at 82-84%, a number unthinkable just a few years ago. These improvements came from increasingly larger databases of protein sequences and structures for training, the use of template secondary structure information and more powerful deep learning techniques. As we are approaching to the theoretical limit of three-state prediction (88-90%), alternative to secondary structure prediction (prediction of backbone torsion angles and Cα-atom-based angles and torsion angles) not only has more room for further improvement but also allows direct prediction of three-dimensional fragment structures with constantly improved accuracy. About 20% of all 40-residue fragments in a database of 1199 non-redundant proteins have <6 Å root-mean-squared distance from the native conformations by SPIDER2. More powerful deep learning methods with improved capability of capturing long-range interactions begin to emerge as the next generation of techniques for secondary structure prediction. The time has come to finish off the final stretch of the long march towards protein secondary structure prediction.

Protein secondary structure analysis of dried blood serum using infrared spectroscopy to identify markers for colitis screening

There remains a great need for diagnosis of inflammatory bowel disease, for which the current technique, colonoscopy, is costly and also has risks for complications. Attenuated total reflectance Fourier transform infrared spectroscopy is a new screening technique to evaluate colitis. Using second derivative spectral deconvolution of the absorbance spectra, a full set of spectral markers were identified based on statistical analysis. Using this method, Amide I group frequencies, (specifically, α-helix to β-sheet ratio of the protein secondary structure) were identified in addition to the previously reported glucose and mannose signatures in sera of chronic and acute mice models of colitis. We also used the same technique to demonstrate that these spectral markers (α-helix/β-sheet ratio, glucose and mannose) are recovering to basal levels upon anti-TNFα therapy. Hence, this technique will be able to identify changes in the sera due to diseases.

Applications of Reactive Cysteine Profiling

Cysteine thiols are involved in a diverse set of biological transformations, including nucleophilic and redox catalysis, metal coordination and formation of both dynamic and structural disulfides. Often posttranslationally modified, cysteines are also frequently alkylated by electrophilic compounds, including electrophilic metabolites, drugs, and natural products, and are attractive sites for covalent probe and drug development. Quantitative proteomics combined with activity-based protein profiling has been applied to annotate cysteine reactivity, susceptibility to posttranslational modifications, and accessibility to chemical probes, uncovering thousands of functional and small-molecule targetable cysteines across a diverse set of proteins, proteome-wide in an unbiased manner. Reactive cysteines have been targeted by high-throughput screening and fragment-based ligand discovery efforts. New cysteine-reactive electrophiles and compound libraries have been synthesized to enable inhibitor discovery broadly and to minimize nonspecific toxicity and off-target activity of compounds. With the recent blockbuster success of several covalent inhibitors, and the development of new chemical proteomic strategies to broadly identify reactive, ligandable and posttranslationally modified cysteines, cysteine profiling is poised to enable the development of new potent and selective chemical probes and even, in some cases, new drugs.

Compound Heterozygous Inheritance of Mutations in Coenzyme Q8A Results in Autosomal Recessive Cerebellar Ataxia and Coenzyme Q10 Deficiency in a Female Sib-Pair

Autosomal recessive ataxias are characterised by a fundamental loss in coordination of gait with associated atrophy of the cerebellum. There is significant clinical and genetic heterogeneity amongst inherited ataxias; however, an early molecular diagnosis is essential with low-risk treatments available for some of these conditions. We describe two female siblings who presented early in life with unsteady gait and cerebellar atrophy. Whole exome sequencing revealed compound heterozygous inheritance of two pathogenic mutations (p.Leu277Pro, c.1506+1G>A) in the coenzyme Q8A gene (COQ8A), a gene central to biosynthesis of coenzyme Q (CoQ). The paternally derived p.Leu277Pro mutation is predicted to disrupt a conserved motif in the substrate-binding pocket of the protein, resulting in inhibition of CoQ₁₀ production. The maternal c.1506+1G>A mutation destroys a canonical splice donor site in exon 12 affecting transcript processing and subsequent protein translation. Mutations in this gene can result in primary coenzyme Q₁₀ deficiency type 4, which is characterized by childhood onset of cerebellar ataxia and exercise intolerance, both of which were observed in this sib-pair. Muscle biopsies revealed unequivocally low levels of CoQ_10, and the siblings were subsequently established on a therapeutic dose of CoQ₁₀ with distinct clinical evidence of improvement after 1 year of treatment. This case emphasises the importance of an early and accurate molecular diagnosis for suspected inherited ataxias, particularly given the availability of approved treatments for some subtypes.

Privileged Electrophile Sensors: A Resource for Covalent Drug Development

This Perspective delineates how redox signaling affects the activity of specific enzyme isoforms and how this property may be harnessed for rational drug design. Covalent drugs have resurged in recent years and several reports have extolled the general virtues of developing irreversible inhibitors. Indeed, many modern pharmaceuticals contain electrophilic appendages. Several invoke a warhead that hijacks active-site nucleophiles whereas others take advantage of spectator nucleophilic side chains that do not participate in enzymatic chemistry, but are poised to bind/react with electrophiles. The latest data suggest that innate electrophile sensing-which enables rapid reaction with an endogenous signaling electrophile-is a quintessential resource for the development of covalent drugs. For instance, based on recent work documenting isoform-specific electrophile sensing, isozyme non-specific drugs may be converted to isozyme-specific analogs by hijacking privileged first-responder electrophile-sensing cysteines. Because this approach targets functionally relevant cysteines, we can simultaneously harness previously untapped moonlighting roles of enzymes linked to redox sensing.

Subcellular Redox Targeting: Bridging in Vitro and in Vivo Chemical Biology

Networks of redox sensor proteins within discrete microdomains regulate the flow of redox signaling. Yet, the inherent reactivity of redox signals complicates the study of specific redox events and pathways by traditional methods. Herein, we review designer chemistries capable of measuring flux and/or mimicking subcellular redox signaling at the cellular and organismal level. Such efforts have begun to decipher the logic underlying organelle-, site-, and target-specific redox signaling in vitro and in vivo. These data highlight chemical biology as a perfect gateway to interrogate how nature choreographs subcellular redox chemistry to drive precision redox biology.

Mutational Spectrum Analysis of Neurodegenerative Diseases and Its Pathogenic Implication

One of the most conspicuous features of neurodegenerative diseases (NDs) is the occurrence of dramatic conformation change of individual proteins. We performed a mutational spectrum analysis of disease-causing missense mutations in seven types of NDs at nucleotide and amino acid levels, and compared the results with those of non-NDs. The main findings included: (i) The higher mutation ratio of G:C→T:A transversion to G:C→A:T transition was observed in NDs than in non-NDs, interpreting the excessive guanine-specific oxidative DNA damage in NDs; (ii) glycine and proline had highest mutability in NDs than in non-NDs, which favor the protein conformation change in NDs; (iii) surprisingly low mutation frequency of arginine was observed in NDs. These findings help to understand how mutations may cause NDs.