Computational approaches to study the effects of small genomic variations

Untapped Potential of Poly(ADP-Ribose) Polymerase Inhibitors: Lessons Learned From the Real-World Clinical Homologous Recombination Repair Mutation Testing

Background

Testing for homologous recombination deficiency (HRD) mutations is pivotal to assess individual risk, to proact preventive measures in healthy carriers and to tailor treatments for cancer patients. Increasing prominence of poly(ADP-ribose) polymerase (PARP) inhibitors with remarkable impact on molecular-selected patient survival across diverse nosologies, ingrains testing for BRCA genes and beyond in clinical practice. Nevertheless, testing strategies remain a question of debate. While several pathogenic BRCA1/2 gene variants have been described as founder pathogenic mutations frequently found in patients from Russia, other homologous recombination repair (HRR) genes have not been sufficiently explored. In this study, we present real-world data of routine HRR gene testing in Russia.

Methods

We evaluated clinical and sequencing data from cancer patients who had germline/somatic next-generation sequencing (NGS) HRR gene testing in Russia (BRCA1/2/ATM/CHEK2, or 15 HRR genes). The primary objectives of this study were to evaluate the frequency of BRCA1/2 and non-BRCA gene mutations in real-world unselected patients from Russia, and to determine whether testing beyond BRCA1/2 is feasible.

Results

Data of 2,032 patients were collected from February 2021 to February 2023. Most had breast (n = 715, 35.2%), ovarian (n = 259, 12.7%), pancreatic (n = 85, 4.2%), or prostate cancer (n = 58, 2.9%). We observed 586 variants of uncertain significance (VUS) and 372 deleterious variants (DVs) across 487 patients, with 17.6% HRR-mutation positivity. HRR testing identified 120 (11.8%) BRCA1/2-positive, and 172 (16.9%) HRR-positive patients. With 51 DVs identified in 242 formalin-fixed paraffin-embedded (FFPE), testing for variant origin clarification was required in one case (0.4%). Most BRCA1/2 germline variants were DV (121 DVs, 26 VUS); in non-BRCA1/2 genes, VUS were ubiquitous (53 DVs, 132 VUS). In silico prediction identified additional 4.9% HRR and 1.2% BRCA1/2/ATM/CHEK2 mutation patients.

Conclusions

Our study represents one of the first reports about the incidence of DV and VUS in HRR genes, including genes beyond BRCA1/2, identified in cancer patients from Russia, assessed by NGS. In silico predictions of the observed HRR gene variants suggest that non-BRCA gene testing is likely to result in higher frequency of patients who are candidates for PARP inhibitor therapy. Continuing sequencing efforts should clarify interpretation of frequently observed non-BRCA VUS.

OptiMo-LDLr: An Integrated In Silico Model with Enhanced Predictive Power for LDL Receptor Variants, Unraveling Hot Spot Pathogenic Residues

Familial hypercholesterolemia (FH) is an inherited metabolic disease affecting cholesterol metabolism, with 90% of cases caused by mutations in the LDL receptor gene (LDLR), primarily missense mutations. This study aims to integrate six commonly used predictive software to create a new model for predicting LDLR mutation pathogenicity and mapping hot spot residues. Six predictive-software are selected: Polyphen-2, SIFT, MutationTaster, REVEL, VARITY, and MLb-LDLr. Software accuracy is tested with the characterized variants annotated in ClinVar and, by bioinformatic and machine learning techniques all models are integrated into a more accurate one. The resulting optimized model presents a specificity of 96.71% and a sensitivity of 98.36%. Hot spot residues with high potential of pathogenicity appear across all domains except for the signal peptide and the O-linked domain. In addition, translating this information into 3D structure of the LDLr highlights potentially pathogenic clusters within the different domains, which may be related to specific biological function. The results of this work provide a powerful tool to classify LDLR pathogenic variants. Moreover, an open-access guide user interface (OptiMo-LDLr) is provided to the scientific community. This study shows that combination of several predictive software results in a more accurate prediction to help clinicians in FH diagnosis.

Exploiting public databases of genomic variation to quantify evolutionary constraint on the branch point sequence in 30 plant and animal species

The branch point sequence is a degenerate intronic heptamer required for the assembly of the spliceosome during pre-mRNA splicing. Disruption of this motif may promote alternative splicing and eventually cause phenotype variation. Despite its functional relevance, the branch point sequence is not included in most genome annotations. Here, we predict branch point sequences in 30 plant and animal species and attempt to quantify their evolutionary constraints using public variant databases. We find an implausible variant distribution in the databases from 16 of 30 examined species. Comparative analysis of variants from whole-genome sequencing shows that variants submitted from exome sequencing or false positive variants are widespread in public databases and cause these irregularities. We then investigate evolutionary constraint with largely unbiased public variant databases in 14 species and find that the fourth and sixth position of the branch point sequence are more constrained than coding nucleotides. Our findings show that public variant databases should be scrutinized for possible biases before they qualify to analyze evolutionary constraint.

In-silico assessment of high-risk non-synonymous SNPs in ADAMTS3 gene associated with Hennekam syndrome and their impact on protein stability and function

Hennekam Lymphangiectasia-Lymphedema Syndrome 3 (HKLLS3) is a rare genetical disorder caused by mutations in a few genes including ADAMTS3. It is characterized by lymphatic dysplasia, intestinal lymphangiectasia, severe lymphedema and distinctive facial appearance. Up till now, no extensive studies have been conducted to elucidate the mechanism of the disease caused by various mutations. As a preliminary investigation of HKLLS3, we sorted out the most deleterious nonsynonymous single nucleotide polymorphisms (nsSNPs) that might affect the structure and function of ADAMTS3 protein by using a variety of in silico tools. A total of 919 nsSNPs in the ADAMTS3 gene were identified. 50 nsSNPs were predicted to be deleterious by multiple computational tools. 5 nsSNPs (G298R, C567Y, A370T, C567R and G374S) were found to be the most dangerous and can be associated with the disease as predicted by different bioinformatics tools. Modelling of the protein shows it can be divided into segments 1, 2 and 3, which are connected by short loops. Segment 3 mainly consists of loops without substantial secondary structures. With prediction tools and molecular dynamics simulation, some SNPs were found to significantly destabilize the protein structure and disrupt the secondary structures, especially in segment 2. The deleterious effects of mutations in segment 1 are possibly not from destabilization but from other factors such as the change in phosphorylation as suggested by post-translational modification (PTM) studies. This is the first-ever study of ADAMTS3 gene polymorphism, and the predicted nsSNPs in ADAMST3, some of which have not been reported yet in patients, will serve for diagnostic purposes and further therapeutic implications in Hennekam syndrome, contributing to better diagnosis and treatment.

In silico analyses of Wnt1 nsSNPs reveal structurally destabilizing variants, altered interactions with Frizzled receptors and its deregulation in tumorigenesis

Wnt1 is the first mammalian Wnt gene, which is discovered as proto-oncogene and in human the gene is located on the chromosome 12q13. Mutations in Wnt1 are reported to be associated with various cancers and other human diseases. The structural and functional consequences of most of the non-synonymous SNPs (nsSNPs), present in the human Wnt1 gene, are not known. In the present work, extensive bioinformatics analyses are used to screen 292 nsSNPs of Wnt1 for predicting pathogenic and harmless polymorphisms. We have identified 10 highly deleterious nsSNPs among which 7 are located within the highly conserved areas. These 10 nsSNPs are also predicted to affect the post-translational modifications of Wnt1. Further, structure based stability analyses of these 10 highly deleterious nsSNPs revealed 8 variants as highly destabilizing. These 8 highly destabilizing variants were shown to have high BC score and high RMSIP score from normal mode analyses. Based on the deformation energies, obtained from the normal mode analyses, variants like G169A, G169S, G331R and G331S were found to be unstable. Molecular Dynamics (MD) simulations revealed structural stability and fluctuation of WT Wnt1 and its prioritized variants. RMSD remained fluctuating mostly between 4 and 5 Å and occasionally between 3.5 and 5.5 Å ranges. RMSF in the CTD region (residues 330–360) of the binding pocket were lower compared to that of WT. Studying the impacts of nsSNPs on the binding interface of Wnt1 and seven Frizzled receptors have predicted substitutions which can stabilize or destabilize the binding interface. We have found that Wnt1 and FZD8-CRD is the best docked complex in our study. MD simulation based analyses of wild type Wnt1-FZD8-CRD complex and the 8 prioritized variants revealed that RMSF was higher in the unstructured regions and RMSD remained fluctuating in the region of 5 Å ± 1 Å. We have also observed differential Wnt1 gene expression pattern in normal, tumor and metastatic conditions across different tissues. Wnt1 gene expression was significantly higher in metastatic tissues of lungs, colon and skin; and was significantly lower in metastatic tissues of breast, esophagus and kidney. We have also found that Wnt1 deregulation is associated with survival outcome in patients with gastric and breast cancer. Furthermore, these computationally screened highly deleterious nsSNPs of Wnt1 can be analyzed in population based genetic studies and may help understand the Wnt1 associated diseases.

A correction into "Theoretical prediction of the trigger linkage, cage strain and explosive sensitivity of CL-20 in the external electric fields"

In order to obtain the more reliable impact sensitivities of CL-20 in the external electric fields, the calculation scheme for the sensitivities shown in the paper published in Journal of Molecular Modeling (entitled "Theoretical prediction of the trigger linkage, cage strain and explosive sensitivity of CL-20 in the external electric fields") was re-evaluated. We found that the model with the averages of the surface electrostatic potentials (ESPs) ([Formula: see text] and [Formula: see text]) may be more suitable for predicting the impact sensitivity of the cage-shaped CL-20 than those containing the variabilities of the surface ESPs ([Formula: see text] and [Formula: see text]) or the balance of charges (ν).

Predicting the Most Deleterious Missense Nonsynonymous Single-Nucleotide Polymorphisms of Hennekam Syndrome-Causing CCBE1 Gene, In Silico Analysis

Hennekam lymphangiectasia-lymphedema syndrome has been linked to single-nucleotide polymorphisms in the CCBE1 (collagen and calcium-binding EGF domains 1) gene. Several bioinformatics methods were used to find the most dangerous nsSNPs that could affect CCBE1 structure and function. Using state-of-the-art in silico tools, this study examined the most pathogenic nonsynonymous single-nucleotide polymorphisms (nsSNPs) that disrupt the CCBE1 protein and extracellular matrix remodeling and migration. Our results indicate that seven nsSNPs, rs115982879, rs149792489, rs374941368, rs121908254, rs149531418, rs121908251, and rs372499913, are deleterious in the CCBE1 gene, four (G330E, C102S, C174R, and G107D) of which are the highly deleterious, two of them (G330E and G107D) have never been seen reported in the context of Hennekam syndrome. Twelve missense SNPs, rs199902030, rs267605221, rs37517418, rs80008675, rs116596858, rs116675104, rs121908252, rs147974432, rs147681552, rs192224843, rs139059968, and rs148498685, are found to revert into stop codons. Structural homology-based methods and sequence homology-based tools revealed that 8.8% of the nsSNPs are pathogenic. SIFT, PolyPhen2, M-CAP, CADD, FATHMM-MKL, DANN, PANTHER, Mutation Taster, LRT, and SNAP2 had a significant score for identifying deleterious nsSNPs. The importance of rs374941368 and rs200149541 in the prediction of post-translation changes was highlighted because it impacts a possible phosphorylation site. Gene-gene interactions revealed CCBE1's association with other genes, showing its role in a number of pathways and coexpressions. The top 16 deleterious nsSNPs found in this research should be investigated further in the future while researching diseases caused CCBE1 gene specifically HS. The FT web server predicted amino acid residues involved in the ligand-binding site of the CCBE1 protein, and two of the substitutions (R167W and T153N) were found to be involved. These highly deleterious nsSNPs can be used as marker pathogenic variants in the mutational diagnosis of the HS syndrome, and this research also offers potential insights that will aid in the development of precision medicines. CCBE1 proteins from Hennekam syndrome patients should be tested in animal models for this purpose.

Molecular dynamics approach to identification of new OGG1 cancer-associated somatic variants with impaired activity

DNA of living cells is always exposed to damaging factors. To counteract the consequences of DNA lesions, cells have evolved several DNA repair systems, among which base excision repair is one of the most important systems. Many currently used antitumor drugs act by damaging DNA, and DNA repair often interferes with chemotherapy and radiotherapy in cancer cells. Tumors are usually extremely genetically heterogeneous, often bearing mutations in DNA repair genes. Thus, knowledge of the functionality of cancer-related variants of proteins involved in DNA damage response and repair is of great interest for personalization of cancer therapy. Although computational methods to predict the variant functionality have attracted much attention, at present, they are mostly based on sequence conservation and make little use of modern capabilities in computational analysis of 3D protein structures. We have used molecular dynamics (MD) to model the structures of 20 clinically observed variants of a DNA repair enzyme, 8-oxoguanine DNA glycosylase. In parallel, we have experimentally characterized the activity, thermostability, and DNA binding in a subset of these mutant proteins. Among the analyzed variants of 8-oxoguanine DNA glycosylase, three (I145M, G202C, and V267M) were significantly functionally impaired and were successfully predicted by MD. Alone or in combination with sequence-based methods, MD may be an important functional prediction tool for cancer-related protein variants of unknown significance.

Electrochemical Detection of Solution Phase Hybridization Related to Single Nucleotide Mutation by Carbon Nanofibers Enriched Electrodes

In the present study, a sensitive and selective impedimetric detection of solution-phase nucleic acid hybridization related to Factor V Leiden (FV Leiden) mutation was performed by carbon nanofibers (CNF) modified screen printed electrodes (SPE). The microscopic and electrochemical characterization of CNF-SPEs was explored in comparison to the unmodified electrodes. Since the FV Leiden mutation is a widespread inherited risk factor predisposing to venous thromboembolism, this study herein aimed to perform the impedimetric detection of FV Leiden mutation by a zip nucleic acid (ZNA) probe-based assay in combination with CNF-SPEs. The selectivity of the assay was then examined against the mutation-free DNA sequences as well as the synthetic PCR samples.

Understanding molecular consequences of putative drug resistant mutations in Mycobacterium tuberculosis

Genomic studies of Mycobacterium tuberculosis bacteria have revealed loci associated with resistance to anti-tuberculosis drugs. However, the molecular consequences of polymorphism within these candidate loci remain poorly understood. To address this, we have used computational tools to quantify the effects of point mutations conferring resistance to three major anti-tuberculosis drugs, isoniazid (n = 189), rifampicin (n = 201) and D-cycloserine (n = 48), within their primary targets, katG, rpoB, and alr. Notably, mild biophysical effects brought about by high incidence mutations were considered more tolerable, while different structural effects brought about by haplotype combinations reflected differences in their functional importance. Additionally, highly destabilising mutations such as alr Y388, highlighted a functional importance of the wildtype residue. Our qualitative analysis enabled us to relate resistance mutations onto a theoretical landscape linking enthalpic changes with phenotype. Such insights will aid the development of new resistance-resistant drugs and, via an integration into predictive tools, in pathogen surveillance.

Targeted sequencing reveals complex, phenotype-correlated genotypes in cystic fibrosis

Background

Cystic fibrosis (CF) is one of the most common life-threatening genetic disorders. Around 2000 variants in the CFTR gene have been identified, with some proportion known to be pathogenic and 300 disease-causing mutations have been characterized in detail by CFTR2 database, which complicates its analysis with conventional methods.

Methods

We conducted next-generation sequencing (NGS) in a cohort of 89 adult patients negative for p.Phe508del homozygosity. Complete clinical and demographic information were available for 84 patients.

Results

By combining MLPA with NGS, we identified disease-causing alleles in all the CF patients. Importantly, in 10% of cases, standard bioinformatics pipelines were inefficient in identifying causative mutations. Class IV-V mutations were observed in 38 (45%) cases, predominantly ones with pancreatic sufficient CF disease; rest of the patients had Class I-III mutations. Diabetes was seen only in patients homozygous for class I-III mutations. We found that 12% of the patients were heterozygous for more than two pathogenic CFTR mutations. Two patients were observed with p.[Arg1070Gln, Ser466*] complex allele which was associated with milder pulmonary obstructions (FVC 107 and 109% versus 67%, CI 95%: 63-72%; FEV 90 and 111% versus 47%, CI 95%: 37-48%). For the first time p.[Phe508del, Leu467Phe] complex allele was reported, observed in four patients (5%).

Conclusion

NGS can be a more information-gaining technology compared to standard methods. Combined with its equivalent diagnostic performance, it can therefore be implemented in the clinical practice, although careful validation is still required.

Electronic supplementary material

The online version of this article (10.1186/s12920-018-0328-z) contains supplementary material, which is available to authorized users.

Vermont: a multi-perspective visual interactive platform for mutational analysis

BACKGROUND

A huge amount of data about genomes and sequence variation is available and continues to grow on a large scale, which makes experimentally characterizing these mutations infeasible regarding disease association and effects on protein structure and function. Therefore, reliable computational approaches are needed to support the understanding of mutations and their impacts. Here, we present VERMONT 2.0, a visual interactive platform that combines sequence and structural parameters with interactive visualizations to make the impact of protein point mutations more understandable.

RESULTS

We aimed to contribute a novel visual analytics oriented method to analyze and gain insight on the impact of protein point mutations. To assess the ability of VERMONT to do this, we visually examined a set of mutations that were experimentally characterized to determine if VERMONT could identify damaging mutations and why they can be considered so.

CONCLUSIONS

VERMONT allowed us to understand mutations by interpreting position-specific structural and physicochemical properties. Additionally, we note some specific positions we believe have an impact on protein function/structure in the case of mutation.

Theoretical insights into the stabilities, detonation performance, and electrostatic potentials of cocrystals containing α- or β-HMX and TATB, FOX-7, NTO, or DMF in various molar ratios

A molecular dynamics method was employed to study the binding energies associated with the cocrystallization (at selected crystal planes) of either 1,3,5-triamino-2,4,6-trinitro-benzene (TATB), 1,1-diamino-2,2-dinitroethylene, 3-nitro-1,2,4-triazol-5-one (TATB, FOX-7, and NTO, respectively, all of which are explosives), or N,N-dimethylformamide (DMF, a nonenergetic solvent) in various molar ratios with 1,3,5,7-tetranitro-1,3,5,7-tetrazacyclooctane in its α and β conformations (α-HMX and β-HMX, respectively). The results showed that the cocrystals with low molar ratios (2:1, 1:1, 1:2, and 1:3) were the most stable. The binding energies of HMX/NTO and HMX/DMF were larger than those of HMX/TATB and HMX/FOX-7. According to the calculated stabilities, HMX prefers to adopt its α form in HMX/TATB and its β form in HMX/NTO, whereas the two forms coexist in HMX/FOX-7. For HMX/TATB, HMX/NTO, and α-HMX/FOX-7, increasing the proportion of the cocrystal component with the highest detonation heat (HMX in the first two cases, FOX-7 in the latter) increases the detonation heat, velocity, and pressure of the cocrystal. However, increasing the proportion of the component with the highest detonation heat in β-HMX/FOX-7 and γ-CL-20/FOX-7 increases the detonation heat of the cocrystal but decreases its detonation velocity. An investigation of the surface electrostatic potential revealed how the sensitivity changes upon cocrystal formation. Graphical Abstract Surface electrostatic potential of HMX/TATB.

The impact of structural genomics: the first quindecennial

The period 2000-2015 brought the advent of high-throughput approaches to protein structure determination. With the overall funding on the order of $2 billion (in 2010 dollars), the structural genomics (SG) consortia established worldwide have developed pipelines for target selection, protein production, sample preparation, crystallization, and structure determination by X-ray crystallography and NMR. These efforts resulted in the determination of over 13,500 protein structures, mostly from unique protein families, and increased the structural coverage of the expanding protein universe. SG programs contributed over 4400 publications to the scientific literature. The NIH-funded Protein Structure Initiatives alone have produced over 2000 scientific publications, which to date have attracted more than 93,000 citations. Software and database developments that were necessary to handle high-throughput structure determination workflows have led to structures of better quality and improved integrity of the associated data. Organized and accessible data have a positive impact on the reproducibility of scientific experiments. Most of the experimental data generated by the SG centers are freely available to the community and has been utilized by scientists in various fields of research. SG projects have created, improved, streamlined, and validated many protocols for protein production and crystallization, data collection, and functional analysis, significantly benefiting biological and biomedical research.

In-Silico Computing of the Most Deleterious nsSNPs in HBA1 Gene

Background

α-Thalassemia (α-thal) is a genetic disorder caused by the substitution of single amino acid or large deletions in the HBA1 and/or HBA2 genes.

Method

Using modern bioinformatics tools as a systematic in-silico approach to predict the deleterious SNPs in the HBA1 gene and its significant pathogenic impact on the functions and structure of HBA1 protein was predicted.

Results and Discussion

A total of 389 SNPs in HBA1 were retrieved from dbSNP database, which includes: 201 non-coding synonymous (nsSNPs), 43 human active SNPs, 16 intronic SNPs, 11 mRNA 3′ UTR SNPs, 9 coding synonymous SNPs, 9 5′ UTR SNPs and other types. Structural homology-based method (PolyPhen) and sequence homology-based tool (SIFT), SNPs&Go, PROVEAN and PANTHER revealed that 2.4% of the nsSNPs are pathogenic.

Conclusions

A total of 5 nsSNPs (G60V, K17M, K17T, L92F and W15R) were predicted to be responsible for the structural and functional modifications of HBA1 protein. It is evident from the deep comprehensive in-silico analysis that, two nsSNPs such as G60Vand W15R in HBA1 are highly deleterious. These “2 pathogenic nsSNPs” can be considered for wet-lab confirmatory analysis.