Search

Melanoma is the third most common malignant skin tumor and has increased in morbidity and mortality over the previous decade due to its rapid spread into the bloodstream or lymphatic system. This study used integrated bioinformatics and network-based methodologies to reliably identify molecular targets and small molecular medicines that may be more successful for Melanoma diagnosis, prognosis and treatment. The statistical LIMMA approach utilized for bioinformatics analysis in this study found 246 common differentially expressed genes (cDEGs) between case and control samples from two microarray gene-expression datasets (GSE130244 and GSE15605). Protein-protein interaction network study revealed 15 cDEGs (PTK2, STAT1, PNO1, CXCR4, WASL, FN1, RUNX2, SOCS3, ITGA4, GNG2, CDK6, BRAF, AGO2, GTF2H1 and AR) to be critical in the development of melanoma (KGs). According to regulatory network analysis, the most important transcriptional and post-transcriptional regulators of DEGs and hub-DEGs are ten transcription factors and three miRNAs. We discovered the pathogenetic mechanisms of MC by studying DEGs' biological processes, molecular function, cellular components and KEGG pathways. We used molecular docking and dynamics modeling to select the four most expressed genes responsible for melanoma malignancy to identify therapeutic candidates. Then, utilizing the Connectivity Map (CMap) database, we analyzed the top 4-hub-DEGs-guided repurposable drugs. We validated four melanoma cancer drugs (Fisetin, Epicatechin Gallate, 1237586-97-8 and PF 431396) using molecular dynamics simulation with their target proteins. As a result, the results of this study may provide resources to researchers and medical professionals for the wet-lab validation of MC diagnosis, prognosis and treatments.Communicated by Ramaswamy H. Sarma.

In the current era of biological big data, which are rapidly populating the biological chemical space, in silico polypharmacology drug design approaches help to decode structure-multiple activity relationships (SMARts). Current computational methods can predict or categorize multiple properties simultaneously, which aids the generation, identification, curation, prioritization, optimization, and repurposing of molecules. Computational methods have generated opportunities and challenges in medicinal chemistry, pharmacology, food chemistry, toxicology, bioinformatics, and chemoinformatics. It is anticipated that computer-guided SMARts could contribute to the full automatization of drug design and drug repurposing campaigns, facilitating the prediction of new biological targets, side and off-target effects, and drug-drug interactions., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

SwissDrugDesign is an important initiative led by the Molecular Modeling Group of the SIB Swiss Institute of Bioinformatics. This project provides a collection of freely available online tools for computer-aided drug design. Some of these web-based methods, i.e., SwissSimilarity and SwissTargetPrediction, were especially developed to perform virtual screening, while others such as SwissADME, SwissDock, SwissParam and SwissBioisostere can find applications in related activities. The present review aims at providing a short description of these methods together with examples of their application in virtual screening, where SwissDrugDesign tools successfully supported the discovery of bioactive small molecules.

Matrix metalloproteinases (MMPs) are a family of zinc-containing enzymes required for homeostasis. These enzymes are an important class of drug targets as their over expression is associated with many disease states. Most of the inhibitors reported against this class of proteins have failed in clinical trials due to lack of specificity. In order to assist in drug design endeavors for MMP targets, a computationally tractable pathway is presented, comprising, (1) docking of small molecule inhibitors against the target MMPs, (2) derivation of quantum mechanical charges on the zinc ion in the active site and the amino acids coordinating with zinc including the inhibitor molecule, (3) molecular dynamics simulations on the docked ligand-MMP complexes, and (4) evaluation of binding affinities of the ligand-MMP complexes via an accurate scoring function for zinc containing metalloprotein-ligand complexes. The above pathway was applied to study the interaction of the inhibitor Batimastat with MMPs, which resulted in a high correlation between the predicted and experimental binding free energies, suggesting the potential applicability of the pathway.

Charcot-Marie-Tooth (CMT) disease is an inherited peripheral neuromuscular disorder characterized by length-dependent and progressive degeneration of peripheral nerves, leading to muscular weakness. Research has shown that mutated HSPB8 may be responsible for depression, neurodegenerative disorders, and improper functioning of peripheral nerves, resulting in neuromuscular disorders like CMT. In the current work, a hybrid approach of virtual screening and molecular docking studies was followed by homology modeling and pharmacophore identification. Detailed screening analyses were carried out by 2-D similarity search against prescribed antidepressant drugs with physicochemical properties. LigandScout was employed to ascertain novel molecules and pharmacophore properties. In this study, we report three novel compounds that showed maximum binding affinity with HSPB8. Docking analysis elucidated that Met37, Ser57, Ser58, Trp60, Thr63, Thr114, Lys115, Asp116, Gly117, Val152, Val154, Leu186, Asp189, Ser190, Gln191, and Glu192 are critical residues for ligand-receptor interactions. Our analyses suggested paroxetine as a potent compound for targeting HSPB8. Selected compounds have more effective energy scores than the selected drug analogs. Additionally, site-directed mutagenesis could be significant for further analysis of the binding pocket. The novel findings based on an in silico approach may be momentous for potent drug design against depression and CMT.

Molecular recognition by protein mostly occurs in a local region on the protein surface. Thus, an efficient computational method for accurate characterization of protein local structural conservation is necessary to better understand biology and drug design. We present a novel local structure alignment tool, G-LoSA. G-LoSA aligns protein local structures in a sequence order independent way and provides a GA-score, a chemical feature-based and size-independent structure similarity score. Our benchmark validation shows the robust performance of G-LoSA to the local structures of diverse sizes and characteristics, demonstrating its universal applicability to local structure-centric comparative biology studies. In particular, G-LoSA is highly effective in detecting conserved local regions on the entire surface of a given protein. In addition, the applications of G-LoSA to identifying template ligands and predicting ligand and protein binding sites illustrate its strong potential for computer-aided drug design. We hope that G-LoSA can be a useful computational method for exploring interesting biological problems through large-scale comparison of protein local structures and facilitating drug discovery research and development. G-LoSA is freely available to academic users at http://im.compbio.ku.edu/GLoSA/., (© 2016 The Protein Society.)

The androgen receptor (AR) is an important regulator of genes responsible for the development and recurrence of prostate cancer. Current therapies for this disease rely on small-molecule inhibitors that block the transcriptional activity of the AR. Recently, major advances in the development of novel AR inhibitors resulted from X-ray crystallographic information on the receptor and utilization of in silico drug design synergized with rigorous experimental testing.Herein, we describe a drug-discovery pipeline for in silico screening for small molecules that target an allosteric region on the AR termed the binding-function 3 (BF3) site. Following the identification of potential candidates, the compounds are tested in cell culture and biochemical assays for their ability to interact with and inhibit the AR. The described pipeline is readily accessible and could be applied in drug design efforts toward any surface-exposed region on the AR or other related steroid nuclear receptor.

BackgroundProtein-protein interactions (PPIs) are intriguing targets for designing novel small-molecule inhibitors. The role of PPIs in various infectious and neurodegenerative disorders makes them potential therapeutic targets . Despite being portrayed as undruggable targets, due to their flat surfaces, disorderedness, and lack of grooves. Recent progresses in computational biology have led researchers to reconsider PPIs in drug discovery.Areas coveredIn this review, we introduce in-silico methods used to identify PPI interfaces and present an in-depth overview of various computational methodologies that are successfully applied to annotate the PPIs. We also discuss several successful case studies that use computational tools to understand PPIs modulation and their key roles in various physiological processes.Expert opinionComputational methods face challenges due to the inherent flexibility of proteins, which makes them expensive, and result in the use of rigid models. This problem becomes more significant in PPIs due to their flexible and flat interfaces. Computational methods like molecular dynamics (MD) simulation and machine learning can integrate the chemical structure data into biochemical and can be used for target identification and modulation. These computational methodologies have been crucial in understanding the structure of PPIs, designing PPI modulators, discovering new drug targets, and predicting treatment outcomes.

The Rat Sarcoma (RAS) family (NRAS, HRAS, and KRAS) is endowed with GTPase activity to regulate various signaling pathways in ubiquitous animal cells. As proto-oncogenes, RAS mutations can maintain activation, leading to the growth and proliferation of abnormal cells and the development of a variety of human cancers. For the fight against tumors, the discovery of RAS-targeted drugs is of high significance. On the one hand, the structural properties of the RAS protein make it difficult to find inhibitors specifically targeted to it. On the other hand, targeting other molecules in the RAS signaling pathway often leads to severe tissue toxicities due to the lack of disease specificity. However, computer-aided drug design (CADD) can help solve the above problems. As an interdisciplinary approach that combines computational biology with medicinal chemistry, CADD has brought a variety of advances and numerous benefits to drug design, such as the rapid identification of new targets and discovery of new drugs. Based on an overview of RAS features and the history of inhibitor discovery, this review provides insight into the application of mainstream CADD methods to RAS drug design. [ABSTRACT FROM AUTHOR]

Background: The investigation and knowledge of calcium handling mechanisms in the plasmodium has been considered as a potential biological target against malaria. Objective: This study deals with the evaluation of inhibitory activity of secondary metabolites of ethylacetate partitioned-fraction of Adansonia digitata stem bark extract on malaria-associated protein using in silico docking studies. Materials and methods: Molecular docking and virtual screening was performed to understand the mechanism of ligand binding and to identify potent calcium transporter inhibitors. The stem bark extracts of A. digitata contains rich sources of phytochemicals. The secondary metabolites were determined by HPLC-DAD and HRGC-MS analysis. The major chemical constituent present in the ethylacetate partitioned-fraction of A. digitata stem bark extract were examined for their antiplasmodial activity and were also involved in docking study. Results: The secondary metabolites, quercetin and apigenin inhibited the formation of β-hematin. The results showed that all the selected compounds in the A. digitata showed binding energy ranging between -6.5 kcal/mol and -7.1 kcal/mol. Among the two chemical constituents, apigenin has the highest docking score along with the highest number of hydrogen bonds formed when compared to quercetin. Analysis of the results suggests that apigenin and quercetin could act as an anti-malaria agent. Conclusion: Molecular docking analysis could lead to further development of potent calcium transporter inhibitors for the prevention and treatment of malaria and related conditions.

There are over 420 human solute carrier (SLC) transporters from 65 families that are expressed ubiquitously in the body. The SLCs mediate the movement of ions, drugs, and metabolites across membranes and their dysfunction has been associated with a variety of diseases, such as diabetes, cancer, and central nervous system (CNS) disorders. Thus, SLCs are emerging as important targets for therapeutic intervention. Recent technological advances in experimental and computational biology allow better characterization of SLC pharmacology. Here we describe recent approaches to modulate SLC transporter function, with an emphasis on the use of computational approaches and computer-aided drug design (CADD) to study nutrient transporters. Finally, we discuss future perspectives in the rational design of SLC drugs. Solute carrier transporters (SLCs) are biomedically important druggable targets but only a few SLCs are targeted by clinically approved drugs. Computer-aided drug design (CADD) has shown promise in the development of tool compounds for SLC transporters. Challenges in CADD for SLCs include the limited number of atomic resolution structures and small-molecule ligands for the development of accurate predictive models. Advances in computational power and machine learning will be powerful in the characterization of SLCs and the design of modulators. [ABSTRACT FROM AUTHOR]

Essential tremor (ET) is one of the most common movement disorders. The etiology of ET remains largely unexplained. Whole genome sequencing (WGS) is likely to be of value in understanding a large proportion of ET with Mendelian and complex disease inheritance patterns. In ET families with Mendelian inheritance patterns, WGS may lead to gene identification where WES analysis failed to identify the causative single nucleotide variant (SNV) or indel due to incomplete coverage of the entire coding region of the genome, in addition to accurate detection of larger structural variants (SVs) and copy number variants (CNVs). Alternatively, in ET families with complex disease inheritance patterns with gene x gene and gene x environment interactions enrichment of functional rare coding and non-coding variants may explain the heritability of ET. We performed WGS in eight ET families (n = 40 individuals) enrolled in the Family Study of Essential Tremor. The analysis included filtering WGS data based on allele frequency in population databases, rare SNV and indel classification and association testing using the Mixed-Model Kernel Based Adaptive Cluster (MM-KBAC) test. A separate analysis of rare SV and CNVs segregating within ET families was also performed. Prioritization of candidate genes identified within families was performed using phenolyzer. WGS analysis identified candidate genes for ET in 5/8 (62.5%) of the families analyzed. WES analysis in a subset of these families in our previously published study failed to identify candidate genes. In one family, we identified a deleterious and damaging variant (c.1367G>A, p.(Arg456Gln)) in the candidate gene, CACNA1G, which encodes the pore forming subunit of T-type Ca(2+) channels, CaV3.1, and is expressed in various motor pathways and has been previously implicated in neuronal autorhythmicity and ET. Other candidate genes identified include SLIT3 which encodes an axon guidance molecule and in three families, phenolyzer prioritized genes that are associated with hereditary neuropathies (family A, KARS, family B, KIF5A and family F, NTRK1). Functional studies of CACNA1G and SLIT3 suggest a role for these genes in ET disease pathogenesis. [ABSTRACT FROM AUTHOR]

Despite ongoing large-scale population-based whole-genome sequencing (WGS) projects such as the NIH NHLBI TOPMed program and the NHGRI Genome Sequencing Program, WGS-based association analysis of complex traits remains a tremendous challenge due to the large number of rare variants, many of which are non-trait-associated neutral variants. External biological knowledge, such as functional annotations based on the ENCODE, Epigenomics Roadmap and GTEx projects, may be helpful in distinguishing causal rare variants from neutral ones; however, each functional annotation can only provide certain aspects of the biological functions. Our knowledge for selecting informative annotations a priori is limited, and incorporating non-informative annotations will introduce noise and lose power. We propose FunSPU, a versatile and adaptive test that incorporates multiple biological annotations and is adaptive at both the annotation and variant levels and thus maintains high power even in the presence of noninformative annotations. In addition to extensive simulations, we illustrate our proposed test using the TWINSUK cohort (n = 1,752) of UK10K WGS data based on six functional annotations: CADD, RegulomeDB, FunSeq, Funseq2, GERP++, and GenoSkyline. We identified genome-wide significant genetic loci on chromosome 19 near gene TOMM40 and APOC4-APOC2 associated with low-density lipoprotein (LDL), which are replicated in the UK10K ALSPAC cohort (n = 1,497). These replicated LDL-associated loci were missed by existing rare variant association tests that either ignore external biological information or rely on a single source of biological knowledge. We have implemented the proposed test in an R package “FunSPU”. [ABSTRACT FROM AUTHOR]

Computational tools are widely used for interpreting variants detected in sequencing projects. The choice of these tools is critical for reliable variant impact interpretation for precision medicine and should be based on systematic performance assessment. The performance of the methods varies widely in different performance assessments, for example due to the contents and sizes of test datasets. To address this issue, we obtained 63,160 common amino acid substitutions (allele frequency ≥1% and <25%) from the Exome Aggregation Consortium (ExAC) database, which contains variants from 60,706 genomes or exomes. We evaluated the specificity, the capability to detect benign variants, for 10 variant interpretation tools. In addition to overall specificity of the tools, we tested their performance for variants in six geographical populations. PON-P2 had the best performance (95.5%) followed by FATHMM (86.4%) and VEST (83.5%). While these tools had excellent performance, the poorest method predicted more than one third of the benign variants to be disease-causing. The results allow choosing reliable methods for benign variant interpretation, for both research and clinical purposes, as well as provide a benchmark for method developers. [ABSTRACT FROM AUTHOR]

There is a growing need to develop variant prediction tools capable of assessing a wide spectrum of evidence. We present a Bayesian framework that involves aggregating pathogenicity data across multiple in silico scores on a gene-by-gene basis and multiple evidence statistics in both quantitative and qualitative forms, and performs 5-tiered variant classification based on the resulting probability credible interval. When evaluated in 1,161 missense variants, our gene-specific in silico model-based meta-predictor yielded an area under the curve (AUC) of 96.0% and outperformed all other in silico predictors. Multifactorial model analysis incorporating all available evidence yielded 99.7% AUC, with 22.8% predicted as variants of uncertain significance (VUS). Use of only 3 auto-computed evidence statistics yielded 98.6% AUC with 56.0% predicted as VUS, which represented sufficient accuracy to rapidly assign a significant portion of VUS to clinically meaningful classifications. Collectively, our findings support the use of this framework to conduct large-scale variant prioritization using in silico predictors followed by variant prediction and classification with a high degree of predictive accuracy. [ABSTRACT FROM AUTHOR]

Many genes have been suggested as candidate genes for keratoconus based on their function, their proximity to associated polymorphisms or due to the identification of putative causative variants within the gene. However, very few of these genes have been assessed for rare variation in keratoconus more broadly. In contrast, VSX1 and SOD1 have been widely assessed, however, the vast majority of studies have been small and the findings conflicting. In a cohort of Australians of European descent, consisting of 385 keratoconus cases and 396 controls, we screened 21 keratoconus candidate genes: BANP, CAST, COL4A3, COL4A4, COL5A1, FOXO1, FNDC3B, HGF, IL1A, IL1B, ILRN, IMMP2L, MPDZ, NFIB, RAB3GAP1, RAD51, RXRA, SLC4A11, SOD1, TF and VSX1. The candidate genes were sequenced in these individuals by either whole exome sequencing or targeted gene sequencing. Variants were filtered to identify rare (minor allele frequency <1%), potentially pathogenic variants. A total of 164 such variants were identified across the two groups with no variants fulfilling these criteria in cases in IL1RN, BANP, IL1B, RAD51 or SOD1. The frequency of variants was compared between cases and controls using chi-square or Fishers’ Exact tests for each gene with at least one rare potentially pathogenic variant identified in the case cohort. The number of rare potentially pathogenic variants per gene ranged from three (RXRA) to 102 (MPDZ), however for all genes, there was no difference in the frequency between the cases and controls. We conclude that rare potentially pathogenic variation in the 21 candidate genes assessed do not play a major role in keratoconus susceptibility and pathogenesis. [ABSTRACT FROM AUTHOR]

Many automatic classifiers were introduced to aid inference of phenotypical effects of uncategorised nsSNVs (nonsynonymous Single Nucleotide Variations) in theoretical and medical applications. Lately, several meta-estimators have been proposed that combine different predictors, such as PolyPhen and SIFT, to integrate more information in a single score. Although many advances have been made in feature design and machine learning algorithms used, the shortage of high-quality reference data along with the bias towards intensively studied in vitro models call for improved generalisation ability in order to further increase classification accuracy and handle records with insufficient data. Since a meta-estimator basically combines different scoring systems with highly complicated nonlinear relationships, we investigated how deep learning (supervised and unsupervised), which is particularly efficient at discovering hierarchies of features, can improve classification performance. While it is believed that one should only use deep learning for high-dimensional input spaces and other models (logistic regression, support vector machines, Bayesian classifiers, etc) for simpler inputs, we still believe that the ability of neural networks to discover intricate structure in highly heterogenous datasets can aid a meta-estimator. We compare the performance with various popular predictors, many of which are recommended by the American College of Medical Genetics and Genomics (ACMG), as well as available deep learning-based predictors. Thanks to hardware acceleration we were able to use a computationally expensive genetic algorithm to stochastically optimise hyper-parameters over many generations. Overfitting was hindered by noise injection and dropout, limiting coadaptation of hidden units. Although we stress that this work was not conceived as a tool comparison, but rather an exploration of the possibilities of deep learning application in ensemble scores, our results show that even relatively simple modern neural networks can significantly improve both prediction accuracy and coverage. We provide open-access to our finest model via the web-site: . [ABSTRACT FROM AUTHOR]

The deleterious genetic variants contributing to certain diseases may differ in terms of number and allele frequency from population to population depending on their evolutionary background. Here, we prioritize the deleterious variants from Pakistani population in manually curated gene list already reported to be associated with common, Mendelian, and congenital cardiovascular diseases (CVDs) using the genome/exome sequencing data of Pakistani individuals publically available in 1000 Genomes Project (PJL), and Exome Aggregation Consortium (ExAC) South Asia. By applying a set of tools such as Combined Annotation Dependent Depletion (CADD), ANNOVAR, and Variant Effect Predictor (VEP), we highlighted 561 potentially detrimental variants from PJL data, and 7374 variants from ExAC South Asian data. Likewise, filtration from ClinVar for CVDs revealed 03 pathogenic and 02 likely pathogenic variants from PJL and 112 pathogenic and 42 likely pathogenic variants from ExAC South Asians. The comparison of derived allele frequencies (DAF) revealed many of these prioritized variants having two fold and higher DAF in Pakistani individuals than in other populations. The highest number of deleterious variants contributing to common CVDs in descending order includes hypertension, atherosclerosis, heart failure, aneurysm, and coronary heart disease, and for Mendelian and congenital CVDs cardiomyopathies, cardiac arrhythmias, and atrioventricular septal defects. [ABSTRACT FROM AUTHOR]

Background: High throughput sequencing technologies have revolutionized the identification of mutations responsible for genetic diseases such as hypertrophic cardiomyopathy (HCM). However, approximately 50% of individuals with a clinical diagnosis of HCM have no causal mutation identified. This may be due to the presence of pathogenic mutations located deep within the introns, which are not detected by conventional sequencing analysis restricted to exons and exon-intron boundaries. Objective: The aim of this study was to develop a whole-gene sequencing strategy to prioritize deep intronic variants that may play a role in HCM pathogenesis. Methods and results: The full genomic DNA sequence of 26 genes previously associated with HCM was analysed in 16 unrelated patients. We identified likely pathogenic deep intronic variants in VCL, PRKAG2 and TTN genes. These variants, which are predicted to act through disruption of either splicing or transcription factor binding sites, are 3-fold more frequent in our cohort of probands than in normal European populations. Moreover, we found a patient that is compound heterozygous for a splice site mutation in MYBPC3 and the deep intronic VCL variant. Analysis of family members revealed that carriers of the MYBPC3 mutation alone do not manifest the disease, while family members that are compound heterozygous are clinically affected. Conclusion: This study provides a framework for scrutinizing variation along the complete intronic sequence of HCM-associated genes and prioritizing candidates for mechanistic and functional analysis. Our data suggest that deep intronic variation contributes to HCM phenotype. [ABSTRACT FROM AUTHOR]

Discriminating the causative disease variant(s) for individuals with inherited or de novo mutations presents one of the main challenges faced by the clinical genetics community today. Computational approaches for variant prioritization include machine learning methods utilizing a large number of features, including molecular information, interaction networks, or phenotypes. Here, we demonstrate the PhenomeNET Variant Predictor (PVP) system that exploits semantic technologies and automated reasoning over genotype-phenotype relations to filter and prioritize variants in whole exome and whole genome sequencing datasets. We demonstrate the performance of PVP in identifying causative variants on a large number of synthetic whole exome and whole genome sequences, covering a wide range of diseases and syndromes. In a retrospective study, we further illustrate the application of PVP for the interpretation of whole exome sequencing data in patients suffering from congenital hypothyroidism. We find that PVP accurately identifies causative variants in whole exome and whole genome sequencing datasets and provides a powerful resource for the discovery of causal variants. [ABSTRACT FROM AUTHOR]

Aims: Dilated cardiomyopathy (DCM) is an important cause of heart failure with a strong familial component. We performed an exome-wide array-based association study (EWAS) to assess the contribution of missense variants to sporadic DCM. Methods and results: 116,855 single nucleotide variants (SNVs) were analyzed in 2796 DCM patients and 6877 control subjects from 6 populations of European ancestry. We confirmed two previously identified associations with SNVs in BAG3 and ZBTB17 and discovered six novel DCM-associated loci (Q-value<0.01). The lead-SNVs at novel loci are common and located in TTN, SLC39A8, MLIP, FLNC, ALPK3 and FHOD3. In silico fine mapping identified HSPB7 as the most likely candidate at the ZBTB17 locus. Rare variant analysis (MAF<0.01) demonstrated significant association for TTN variants only (P = 0.0085). All candidate genes but one (SLC39A8) exhibit preferential expression in striated muscle tissues and mutations in TTN, BAG3, FLNC and FHOD3 are known to cause familial cardiomyopathy. We also investigated a panel of 48 known cardiomyopathy genes. Collectively, rare (n = 228, P = 0.0033) or common (n = 36, P = 0.019) variants with elevated in silico severity scores were associated with DCM, indicating that the spectrum of genes contributing to sporadic DCM extends beyond those identified here. Conclusion: We identified eight loci independently associated with sporadic DCM. The functions of the best candidate genes at these loci suggest that proteostasis regulation might play a role in DCM pathophysiology. [ABSTRACT FROM AUTHOR]

Background: Deep vein thrombosis (DVT) genetic predisposition is partially known. Objectives: This study aimed at assessing the functional impact of nine ADAMTS13 single nucleotide variants (SNVs) previously reported to be associated as a group with DVT in a burden test and the individual association of selected variants with DVT risk in two replication studies. Methods: Wild-type and mutant recombinant ADAMTS13 were transiently expressed in HEK293 cells. Antigen and activity of recombinant ADAMTS13 were measured by ELISA and FRETS-VWF73 assays, respectively. The replication studies were performed in an Italian case-control study (Milan study; 298/298 patients/controls) using a next-generation sequencing approach and in a Dutch case-control study (MEGA study; 4306/4887 patients/controls) by TaqMan assays. Results: In vitro results showed reduced ADAMTS13 activity for three SNVs (p.Val154Ile [15%; 95% confidence interval [CI] 14–16], p.Asp187His [19%; 95%[CI] 17–21], p.Arg421Cys [24%; 95%[CI] 22–26]) similar to reduced plasma ADAMTS13 levels of patients carriers for these SNVs. Therefore these three SNVs were interrogated for risk association. The first replication study identified 3 heterozygous carriers (2 cases, 1 control) of p.Arg421Cys (odds ratio [OR] 2, 95%[CI] 0.18–22.25). The second replication study identified 2 heterozygous carriers (1 case, 1 control) of p.Asp187His ([OR] 1.14, 95%[CI] 0.07–18.15) and 10 heterozygous carriers (4 cases, 6 controls) of p.Arg421Cys ([OR] 0.76, 95%[CI] 0.21–2.68). Conclusions: Three SNVs (p.Val154Ile, p.Asp187His and p.Arg421Cys) showed reduced ex vivo and in vitro ADAMTS13 levels. However, the low frequency of these variants makes it difficult to confirm their association with DVT. [ABSTRACT FROM AUTHOR]

Left-ventricular outflow tract obstructions (LVOTO) encompass a wide spectrum of phenotypically heterogeneous heart malformations which frequently cluster in families. We performed family based whole-exome and targeted re-sequencing on 182 individuals from 51 families with multiple affected members. Central to our approach is the family unit which serves as a reference to identify causal genotype-phenotype correlations. Screening a multitude of 10 overlapping phenotypes revealed disease associated and co-segregating variants in 12 families. These rare or novel protein altering mutations cluster predominantly in genes (NOTCH1, ARHGAP31, MAML1, SMARCA4, JARID2, JAG1) along the Notch signaling cascade. This is in line with a significant enrichment (Wilcoxon, p< 0.05) of variants with a higher pathogenicity in the Notch signaling pathway in patients compared to controls. The significant enrichment of novel protein truncating and missense mutations in NOTCH1 highlights the allelic and phenotypic heterogeneity in our pediatric cohort. We identified novel co-segregating pathogenic mutations in NOTCH1 associated with left and right-sided cardiac malformations in three independent families with a total of 15 affected individuals. In summary, our results suggest that a small but highly pathogenic fraction of family specific mutations along the Notch cascade are a common cause of LVOTO. [ABSTRACT FROM AUTHOR]

An important message taken from human genome sequencing projects is that the human population exhibits approximately 99.9% genetic similarity. Variations in the remaining parts of the genome determine our identity, trace our history and reveal our heritage. The precise delineation of phenotypically causal variants plays a key role in providing accurate personalized diagnosis, prognosis, and treatment of inherited diseases. Several computational methods for achieving such delineation have been reported recently. However, their ability to pinpoint potentially deleterious variants is limited by the fact that their mechanisms of prediction do not account for the existence of different categories of variants. Consequently, their output is biased towards the variant categories that are most strongly represented in the variant databases. Moreover, most such methods provide numeric scores but not binary predictions of the deleteriousness of variants or confidence scores that would be more easily understood by users. We have constructed three datasets covering different types of disease-related variants, which were divided across five categories: (i) regulatory, (ii) splicing, (iii) missense, (iv) synonymous, and (v) nonsense variants. These datasets were used to develop category-optimal decision thresholds and to evaluate six tools for variant prioritization: CADD, DANN, FATHMM, FitCons, FunSeq2 and GWAVA. This evaluation revealed some important advantages of the category-based approach. The results obtained with the five best-performing tools were then combined into a consensus score. Additional comparative analyses showed that in the case of missense variations, protein-based predictors perform better than DNA sequence-based predictors. A user-friendly web interface was developed that provides easy access to the five tools’ predictions, and their consensus scores, in a user-understandable format tailored to the specific features of different categories of variations. To enable comprehensive evaluation of variants, the predictions are complemented with annotations from eight databases. The web server is freely available to the community at . [ABSTRACT FROM AUTHOR]

Background: It has become common practice to analyse large scale sequencing data with statistical approaches based around the aggregation of rare variants within the same gene. We applied a novel approach to rare variant analysis by collapsing variants together using protein domain and family coordinates, regarded to be a more discrete definition of a biologically functional unit. Methods: Using Pfam definitions, we collapsed rare variants (Minor Allele Frequency ≤ 1%) together in three different ways 1) variants within single genomic regions which map to individual protein domains 2) variants within two individual protein domain regions which are predicted to be responsible for a protein-protein interaction 3) all variants within combined regions from multiple genes responsible for coding the same protein domain (i.e. protein families). A conventional collapsing analysis using gene coordinates was also undertaken for comparison. We used UK10K sequence data and investigated associations between regions of variants and lipid traits using the sequence kernel association test (SKAT). Results: We observed no strong evidence of association between regions of variants based on Pfam domain definitions and lipid traits. Quantile-Quantile plots illustrated that the overall distributions of p-values from the protein domain analyses were comparable to that of a conventional gene-based approach. Deviations from this distribution suggested that collapsing by either protein domain or gene definitions may be favourable depending on the trait analysed. Conclusion: We have collapsed rare variants together using protein domain and family coordinates to present an alternative approach over collapsing across conventionally used gene-based regions. Although no strong evidence of association was detected in these analyses, future studies may still find value in adopting these approaches to detect previously unidentified association signals. [ABSTRACT FROM AUTHOR]

Background: The success of collapsing methods which investigate the combined effect of rare variants on complex traits has so far been limited. The manner in which variants within a gene are selected prior to analysis has a crucial impact on this success, which has resulted in analyses conventionally filtering variants according to their consequence. This study investigates whether an alternative approach to filtering, using annotations from recently developed bioinformatics tools, can aid these types of analyses in comparison to conventional approaches. Methods & Results: We conducted a candidate gene analysis using the UK10K sequence and lipids data, filtering according to functional annotations using the resource CADD (Combined Annotation-Dependent Depletion) and contrasting results with ‘nonsynonymous’ and ‘loss of function’ consequence analyses. Using CADD allowed the inclusion of potentially deleterious intronic variants, which was not possible when filtering by consequence. Overall, different filtering approaches provided similar evidence of association, although filtering according to CADD identified evidence of association between ANGPTL4 and High Density Lipoproteins (P = 0.02, N = 3,210) which was not observed in the other analyses. We also undertook genome-wide analyses to determine how filtering in this manner compared to conventional approaches for gene regions. Results suggested that filtering by annotations according to CADD, as well as other tools known as FATHMM-MKL and DANN, identified association signals not detected when filtering by variant consequence and vice versa. Conclusion: Incorporating variant annotations from non-coding bioinformatics tools should prove to be a valuable asset for rare variant analyses in the future. Filtering by variant consequence is only possible in coding regions of the genome, whereas utilising non-coding bioinformatics annotations provides an opportunity to discover unknown causal variants in non-coding regions as well. This should allow studies to uncover a greater number of causal variants for complex traits and help elucidate their functional role in disease. [ABSTRACT FROM AUTHOR]

Diabetes mellitus is a debilitating metabolic disorder and remains a significant threat to public health. Herbal medicines have been proven to be effective anti-diabetic agents compared to synthetic drugs in terms of side effects. However, the complexity in their chemical constituents and mechanism of action, hinder the effort to discover novel anti-diabetic drugs. Hence, understanding the biological and chemical basis of pharmacological action of phytochemicals is essential for the discovery of potential anti-diabetic drugs. Identifying important active compounds, their protein targets and the pathways involved in diabetes would serve this purpose. In this context, the present study was aimed at exploring the mechanism of action of anti-diabetic plants phytochemicals through network and chemical-based approaches. This study also involves a focused and constructive strategy for preparing new effective anti-diabetic formulations. Further, a protocol for target enrichment was proposed, to identify novel protein targets for important active compounds. Therefore, the successive use of network analysis combined with target enrichment studies would accelerate the discovery of potential anti-diabetic phytochemicals. [ABSTRACT FROM AUTHOR]

Late-stage or post-market identification of adverse drug reactions (ADRs) is a significant public health issue and a source of major economic liability for drug development. Thus, reliable in silico screening of drug candidates for possible ADRs would be advantageous. In this work, we introduce a computational approach that predicts ADRs by combining the results of molecular docking and leverages known ADR information from DrugBank and SIDER. We employed a recently parallelized version of AutoDock Vina (VinaLC) to dock 906 small molecule drugs to a virtual panel of 409 DrugBank protein targets. L1-regularized logistic regression models were trained on the resulting docking scores of a 560 compound subset from the initial 906 compounds to predict 85 side effects, grouped into 10 ADR phenotype groups. Only 21% (87 out of 409) of the drug-protein binding features involve known targets of the drug subset, providing a significant probe of off-target effects. As a control, associations of this drug subset with the 555 annotated targets of these compounds, as reported in DrugBank, were used as features to train a separate group of models. The Vina off-target models and the DrugBank on-target models yielded comparable median area-under-the-receiver-operating-characteristic-curves (AUCs) during 10-fold cross-validation (0.60–0.69 and 0.61–0.74, respectively). Evidence was found in the PubMed literature to support several putative ADR-protein associations identified by our analysis. Among them, several associations between neoplasm-related ADRs and known tumor suppressor and tumor invasiveness marker proteins were found. A dual role for interstitial collagenase in both neoplasms and aneurysm formation was also identified. These associations all involve off-target proteins and could not have been found using available drug/on-target interaction data. This study illustrates a path forward to comprehensive ADR virtual screening that can potentially scale with increasing number of CPUs to tens of thousands of protein targets and millions of potential drug candidates. [ABSTRACT FROM AUTHOR]

Quantum mechanical (QM) calculations have been used to predict the binding affinity of a set of ligands towards HIV-1 RT associated RNase H (RNH). The QM based chelation calculations show improved binding affinity prediction for the inhibitors compared to using an empirical scoring function. Furthermore, full protein fragment molecular orbital (FMO) calculations were conducted and subsequently analysed for individual residue stabilization/destabilization energy contributions to the overall binding affinity in order to better understand the true and false predictions. After a successful assessment of the methods based on the use of a training set of molecules, QM based chelation calculations were used as filter in virtual screening of compounds in the ZINC database. By this, we find, compared to regular docking, QM based chelation calculations to significantly reduce the large number of false positives. Thus, the computational models tested in this study could be useful as high throughput filters for searching HIV-1 RNase H active-site molecules in the virtual screening process. [ABSTRACT FROM AUTHOR]