Accurate protein stability predictions from homology models

Prospective Evaluation of Structure-Based Simulations Reveal Their Ability to Predict the Impact of Kinase Mutations on Inhibitor Binding

Small molecule kinase inhibitors are critical in the modern treatment of cancers, evidenced by the existence of over 80 FDA-approved small-molecule kinase inhibitors. Unfortunately, intrinsic or acquired resistance, often causing therapy discontinuation, is frequently caused by mutations in the kinase therapeutic target. The advent of clinical tumor sequencing has opened additional opportunities for precision oncology to improve patient outcomes by pairing optimal therapies with tumor mutation profiles. However, modern precision oncology efforts are hindered by lack of sufficient biochemical or clinical evidence to classify each mutation as resistant or sensitive to existing inhibitors. Structure-based methods show promising accuracy in retrospective benchmarks at predicting whether a kinase mutation will perturb inhibitor binding, but comparisons are made by pooling disparate experimental measurements across different conditions. We present the first prospective benchmark of structure-based approaches on a blinded dataset of in-cell kinase inhibitor affinities to Abl kinase mutants using a NanoBRET reporter assay. We compare NanoBRET results to structure-based methods and their ability to estimate the impact of mutations on inhibitor binding (measured as ΔΔG). Comparing physics-based simulations, Rosetta, and previous machine learning models, we find that structure-based methods accurately classify kinase mutations as inhibitor-resistant or inhibitor-sensitizing, and each approach has a similar degree of accuracy. We show that physics-based simulations are best suited to estimate ΔΔG of mutations that are distal to the kinase active site. To probe modes of failure, we retrospectively investigate two clinically significant mutations poorly predicted by our methods, T315A and L298F, and find that starting configurations and protonation states significantly alter the accuracy of our predictions. Our experimental and computational measurements provide a benchmark for estimating the impact of mutations on inhibitor binding affinity for future methods and structure-based models. These structure-based methods have potential utility in identifying optimal therapies for tumor-specific mutations, predicting resistance mutations in the absence of clinical data, and identifying potential sensitizing mutations to established inhibitors.

Prospective evaluation of structure-based simulations reveal their ability to predict the impact of kinase mutations on inhibitor binding

Small molecule kinase inhibitors are critical in the modern treatment of cancers, evidenced by the existence of over 80 FDA-approved small-molecule kinase inhibitors. Unfortunately, intrinsic or acquired resistance, often causing therapy discontinuation, is frequently caused by mutations in the kinase therapeutic target. The advent of clinical tumor sequencing has opened additional opportunities for precision oncology to improve patient outcomes by pairing optimal therapies with tumor mutation profiles. However, modern precision oncology efforts are hindered by lack of sufficient biochemical or clinical evidence to classify each mutation as resistant or sensitive to existing inhibitors. Structure-based methods show promising accuracy in retrospective benchmarks at predicting whether a kinase mutation will perturb inhibitor binding, but comparisons are made by pooling disparate experimental measurements across different conditions. We present the first prospective benchmark of structure-based approaches on a blinded dataset of in-cell kinase inhibitor affinities to Abl kinase mutants using a NanoBRET reporter assay. We compare NanoBRET results to structure-based methods and their ability to estimate the impact of mutations on inhibitor binding (measured as ΔΔG). Comparing physics-based simulations, Rosetta, and previous machine learning models, we find that structure-based methods accurately classify kinase mutations as inhibitor-resistant or inhibitor-sensitizing, and each approach has a similar degree of accuracy. We show that physics-based simulations are best suited to estimate ΔΔG of mutations that are distal to the kinase active site. To probe modes of failure, we retrospectively investigate two clinically significant mutations poorly predicted by our methods, T315A and L298F, and find that starting configurations and protonation states significantly alter the accuracy of our predictions. Our experimental and computational measurements provide a benchmark for estimating the impact of mutations on inhibitor binding affinity for future methods and structure-based models. These structure-based methods have potential utility in identifying optimal therapies for tumor-specific mutations, predicting resistance mutations in the absence of clinical data, and identifying potential sensitizing mutations to established inhibitors.

An allelic atlas of immunoglobulin heavy chain variable regions reveals antibody binding epitope preference resilient to SARS-CoV-2 mutation escape

Background

Although immunoglobulin (Ig) alleles play a pivotal role in the antibody response to pathogens, research to understand their role in the humoral immune response is still limited.

Methods

We retrieved the germline sequences for the IGHV from the IMGT database to illustrate the amino acid polymorphism present within germline sequences of IGHV genes. We aassembled the sequences of IgM and IgD repertoire from 130 people to investigate the genetic variations in the population. A dataset comprising 10,643 SARS-CoV-2 spike-specific antibodies, obtained from COV-AbDab, was compiled to assess the impact of SARS-CoV-2 infection on allelic gene utilization. Binding affinity and neutralizing activity were determined using bio-layer interferometry and pseudovirus neutralization assays. Primary docking was performed using ZDOCK (3.0.2) to generate the initial conformation of the antigen-antibody complex, followed by simulations of the complete conformations using Rosetta SnugDock software. The original and simulated structural conformations were visualized and presented using ChimeraX (v1.5).

Results

We present an allelic atlas of immunoglobulin heavy chain (IgH) variable regions, illustrating the diversity of allelic variants across 33 IGHV family germline sequences by sequencing the IgH repertoire of in the population. Our comprehensive analysis of SARS-CoV-2 spike-specific antibodies revealed the preferential use of specific Ig alleles among these antibodies. We observed an association between Ig alleles and antibody binding epitopes. Different allelic genotypes binding to the same RBD epitope on the spike show different neutralizing potency and breadth. We found that antibodies carrying the IGHV1-69*02 allele tended to bind to the RBD E2.2 epitope. The antibodies carrying G50 and L55 amino acid residues exhibit potential enhancements in binding affinity and neutralizing potency to SARS-CoV-2 variants containing the L452R mutation on RBD, whereas R50 and F55 amino acid residues tend to have reduced binding affinity and neutralizing potency. IGHV2-5*02 antibodies using the D56 allele bind to the RBD D2 epitope with greater binding and neutralizing potency due to the interaction between D56 on HCDR2 and K444 on RBD of most Omicron subvariants. In contrast, IGHV2-5*01 antibodies using the N56 allele show increased binding resistance to the K444T mutation on RBD.

Discussion

This study provides valuable insights into humoral immune responses from the perspective of Ig alleles and population genetics. These findings underscore the importance of Ig alleles in vaccine design and therapeutic antibody development.

Exploration of Fingerprints and Data Mining-based Prediction of Some Bioactive Compounds from Allium sativum as Histone Deacetylase 9 (HDAC9) Inhibitors

BACKGROUND

Histone deacetylase 9 (HDAC9) is an important member of the class IIa family of histone deacetylases. It is well established that over-expression of HDAC9 causes various types of cancers including gastric cancer, breast cancer, ovarian cancer, liver cancer, lung cancer, lymphoblastic leukaemia, etc. The important role of HDAC9 is also recognized in the development of bone, cardiac muscles, and innate immunity. Thus, it will be beneficial to find out the important structural attributes of HDAC9 inhibitors for developing selective HDAC9 inhibitors with higher potency.

METHODS

The classification QSAR-based methods namely Bayesian classification and recursive partitioning method were applied to a dataset consisting of HADC9 inhibitors. The structural features strongly suggested that sulphur-containing compounds can be a good choice for HDAC9 inhibition. For this reason, these models were applied further to screen some natural compounds from Allium sativum. The screened compounds were further accessed for the ADME properties and docked in the homology-modelled structure of HDAC9 in order to find important amino acids for the interaction. The best-docked compound was considered for molecular dynamics (MD) simulation study.

RESULTS

The classification models have identified good and bad fingerprints for HDAC9 inhibition. The screened compounds like ajoene, 1,2 vinyl dithiine, diallyl disulphide and diallyl trisulphide had been identified as compounds having potent HDAC9 inhibitory activity. The results from ADME and molecular docking study of these compounds show the binding interaction inside the active site of the HDAC9. The best-docked compound ajoene shows satisfactory results in terms of different validation parameters of MD simulation study.

CONCLUSION

This in-silico modelling study has identified the natural potential lead (s) from Allium sativum. Specifically, the ajoene with the best in-silico features can be considered for further in-vitro and in-vivo investigation to establish as potential HDAC9 inhibitors.

An overview on olfaction in the biological, analytical, computational, and machine learning fields

Recently, the comprehension of odor perception has advanced, unveiling the mysteries of the molecular receptors within the nasal passages and the intricate mechanisms governing signal transmission between these receptors, the olfactory bulb, and the brain. This review provides a comprehensive panorama of odors, encompassing various topics ranging from the structural and molecular underpinnings of odorous substances to the physiological intricacies of olfactory perception. It extends to elucidate the analytical methods used for their identification and explores the frontiers of computational methodologies.

Vaccine development using artificial intelligence and machine learning: A review

The COVID-19 pandemic has underscored the critical importance of effective vaccines, yet their development is a challenging and demanding process. It requires identifying antigens that elicit protective immunity, selecting adjuvants that enhance immunogenicity, and designing delivery systems that ensure optimal efficacy. Artificial intelligence (AI) can facilitate this process by using machine learning methods to analyze large and diverse datasets, suggest novel vaccine candidates, and refine their design and predict their performance. This review explores how AI can be applied to various aspects of vaccine development, such as predicting immune response from protein sequences, discovering adjuvants, optimizing vaccine doses, modeling vaccine supply chains, and predicting protein structures. We also address the challenges and ethical issues that emerge from the use of AI in vaccine development, such as data privacy, algorithmic bias, and health data sensitivity. We contend that AI has immense potential to accelerate vaccine development and respond to future pandemics, but it also requires careful attention to the quality and validity of the data and methods used.

Using AI-predicted protein structures as a reference to predict loss-of-function activity in tumor suppressor breast cancer genes

Background

The loss-of-function (LOF) classification of most missense variants in tumor suppressor breast cancer genes BRCA1, BRCA2, PALB2, and RAD51C remains unclassified and confounds clinical actionability. Classifying these variants is challenging due to their rarity, leading clinicians to rely on in silico predictive methods. Protein stability changes are associated with function, making stability predictors valuable. Stability predictions upon missense variant perturbations require high-resolution protein structures. However, the availability of these high-resolution structures is lacking. This study explores using generative AI to predict high-resolution protein structures, which can then be analyzed with in silico protein stability prediction methods to assess LOF activity in ordered regions of the protein. This study also determines the appropriate in silico protein stability and dedicated in silico missense prediction methods in dbNSFP v4.7 database to predict LOF activity in ordered regions of these four genes. Functional classifications from homology recombination DNA repair (HDR) assays and variant classifications from the ClinVar database provide a reliable dataset for evaluating the performance of these in silico prediction methods.

Results

Complex AlphaFold2 structures of the BRCA1-C terminal (BRCT) domain and the DNA-binding (DB) domain of BRCA2, analyzed using protein stability tool FoldX predicts LOF activity from missense variants significantly better than experimentally-derived structures in ordered regions. The BRCT domain achieved an Area Under the Curve (AUC)= 0.861 (95 % CI:0.858-0.863) and AUC= 0.842 (95 % CI:0.840-0.845), while the DB domain achieved an AUC= 0.836 (95 % CI:0.8322-0.841), compared to AUC= 0.847 (95 % CI:0.844-0.850) and AUC= 0.835 (95 % CI:0.832-0.837) from the BRCT domain, and AUC= 0.830 (95 % CI:0.821-0.8320) from the DB domain from experimentally-derived structures. Protein stability does not predict LOF activity from missense variants better than dedicated in silico missense predictors. Overall, we find that AlphaMissense ranks highly, with an average AUC= 0.890 (95 % CI 0.886-0.895) from ordered regions across these four cancer genes, compared to all other in silico missense predictors present in the dbNSFP database.

Conclusions

The study reveals that generative AI protein predicted structures can outperform experimentally-derived structures in evaluating LOF activity from predicted protein stability in ordered regions of genes BRCA1, BRCA2, PALB2 and RAD51C. The study also highlights the predictive performance of AlphaMissense as the premier in silico missense prediction method to predict LOF activity from missense variants in these four tumor suppressor breast cancer genes. The code for this study can be downloaded for free on GitHub (https://github.com/rohandavidg/CarePred).

Integrative analysis of pathogenic variants in glucose-6-phosphatase based on an AlphaFold2 model

Mediating the terminal reaction of gluconeogenesis and glycogenolysis, the integral membrane protein glucose-6-phosphate catalytic subunit 1 (G6PC1) regulates hepatic glucose production by catalyzing hydrolysis of glucose-6-phosphate (G6P) within the lumen of the endoplasmic reticulum. Consistent with its vital contribution to glucose homeostasis, inactivating mutations in G6PC1 causes glycogen storage disease (GSD) type 1a characterized by hepatomegaly and severe hypoglycemia. Despite its physiological importance, the structural basis of G6P binding to G6PC1 and the molecular disruptions induced by missense mutations within the active site that give rise to GSD type 1a are unknown. In this study, we determine the atomic interactions governing G6P binding as well as explore the perturbations imposed by disease-linked missense variants by subjecting an AlphaFold2 G6PC1 structural model to molecular dynamics simulations and in silico predictions of thermodynamic stability validated with robust in vitro and in situ biochemical assays. We identify a collection of side chains, including conserved residues from the signature phosphatidic acid phosphatase motif, that contribute to a hydrogen bonding and van der Waals network stabilizing G6P in the active site. The introduction of GSD type 1a mutations modified the thermodynamic landscape, altered side chain packing and substrate-binding interactions, and induced trapping of catalytic intermediates. Our results, which corroborate the high quality of the AF2 model as a guide for experimental design and to interpret outcomes, not only confirm the active-site structural organization but also identify previously unobserved mechanistic contributions of catalytic and noncatalytic side chains.

Prediction of cancer driver genes and mutations: the potential of integrative computational frameworks

The vast amount of available sequencing data allows the scientific community to explore different genetic alterations that may drive cancer or favor cancer progression. Software developers have proposed a myriad of predictive tools, allowing researchers and clinicians to compare and prioritize driver genes and mutations and their relative pathogenicity. However, there is little consensus on the computational approach or a golden standard for comparison. Hence, benchmarking the different tools depends highly on the input data, indicating that overfitting is still a massive problem. One of the solutions is to limit the scope and usage of specific tools. However, such limitations force researchers to walk on a tightrope between creating and using high-quality tools for a specific purpose and describing the complex alterations driving cancer. While the knowledge of cancer development increases daily, many bioinformatic pipelines rely on single nucleotide variants or alterations in a vacuum without accounting for cellular compartments, mutational burden or disease progression. Even within bioinformatics and computational cancer biology, the research fields work in silos, risking overlooking potential synergies or breakthroughs. Here, we provide an overview of databases and datasets for building or testing predictive cancer driver tools. Furthermore, we introduce predictive tools for driver genes, driver mutations, and the impact of these based on structural analysis. Additionally, we suggest and recommend directions in the field to avoid silo-research, moving towards integrative frameworks.

Systematic identification of conditionally folded intrinsically disordered regions by AlphaFold2

The AlphaFold Protein Structure Database contains predicted structures for millions of proteins. For the majority of human proteins that contain intrinsically disordered regions (IDRs), which do not adopt a stable structure, it is generally assumed that these regions have low AlphaFold2 confidence scores that reflect low-confidence structural predictions. Here, we show that AlphaFold2 assigns confident structures to nearly 15% of human IDRs. By comparison to experimental NMR data for a subset of IDRs that are known to conditionally fold (i.e., upon binding or under other specific conditions), we find that AlphaFold2 often predicts the structure of the conditionally folded state. Based on databases of IDRs that are known to conditionally fold, we estimate that AlphaFold2 can identify conditionally folding IDRs at a precision as high as 88% at a 10% false positive rate, which is remarkable considering that conditionally folded IDR structures were minimally represented in its training data. We find that human disease mutations are nearly fivefold enriched in conditionally folded IDRs over IDRs in general and that up to 80% of IDRs in prokaryotes are predicted to conditionally fold, compared to less than 20% of eukaryotic IDRs. These results indicate that a large majority of IDRs in the proteomes of human and other eukaryotes function in the absence of conditional folding, but the regions that do acquire folds are more sensitive to mutations. We emphasize that the AlphaFold2 predictions do not reveal functionally relevant structural plasticity within IDRs and cannot offer realistic ensemble representations of conditionally folded IDRs.

Molecular mechanisms of catalytic inhibition for active site mutations in glucose-6-phosphatase catalytic subunit 1 linked to glycogen storage disease

Mediating the terminal reaction of gluconeogenesis and glycogenolysis, the integral membrane protein G6PC1 regulates hepatic glucose production by catalyzing hydrolysis of glucose-6-phosphate (G6P) within the lumen of the endoplasmic reticulum. Consistent with its vital contribution to glucose homeostasis, inactivating mutations in G6PC1 cause glycogen storage disease (GSD) type 1a characterized by hepatomegaly and severe hypoglycemia. Despite its physiological importance, the structural basis of G6P binding to G6PC1 and the molecular disruptions induced by missense mutations within the active site that give rise to GSD type 1a are unknown. Exploiting a computational model of G6PC1 derived from the groundbreaking structure prediction algorithm AlphaFold2 (AF2), we combine molecular dynamics (MD) simulations and computational predictions of thermodynamic stability with a robust in vitro screening platform to define the atomic interactions governing G6P binding as well as explore the energetic perturbations imposed by disease-linked variants. We identify a collection of side chains, including conserved residues from the signature phosphatidic acid phosphatase motif, that contribute to a hydrogen bonding and van der Waals network stabilizing G6P in the active site. Introduction of GSD type 1a mutations into the G6PC1 sequence elicits changes in G6P binding energy, thermostability and structural properties, suggesting multiple pathways of catalytic impairment. Our results, which corroborate the high quality of the AF2 model as a guide for experimental design and to interpret outcomes, not only confirm active site structural organization but also suggest novel mechanistic contributions of catalytic and non-catalytic side chains.

Rapid protein stability prediction using deep learning representations

Predicting the thermodynamic stability of proteins is a common and widely used step in protein engineering, and when elucidating the molecular mechanisms behind evolution and disease. Here, we present RaSP, a method for making rapid and accurate predictions of changes in protein stability by leveraging deep learning representations. RaSP performs on-par with biophysics-based methods and enables saturation mutagenesis stability predictions in less than a second per residue. We use RaSP to calculate ∼ 230 million stability changes for nearly all single amino acid changes in the human proteome, and examine variants observed in the human population. We find that variants that are common in the population are substantially depleted for severe destabilization, and that there are substantial differences between benign and pathogenic variants, highlighting the role of protein stability in genetic diseases. RaSP is freely available-including via a Web interface-and enables large-scale analyses of stability in experimental and predicted protein structures.