DeMaSk: a deep mutational scanning substitution matrix and its use for variant impact prediction

Enhanced Biosynthesis of 6'-Sialyllactose in Escherichia coli via Systematic Metabolic Engineering

6'-Sialyllactose (6'-SL), a sialylated oligosaccharide belonging to human milk oligosaccharides (HMOs), finds extensive applications in food, pharmaceuticals, cosmetics, and particularly in infant nutrition and dietary supplements. In this study, Escherichia coli K12 MG1655 was subjected to metabolic engineering to enhance 6'-SL biosynthesis. Initial pathway optimization involved the deletion of competing metabolic pathway genes (nanATEK, nagAB, and lacZ). Following overexpression of NmCSS and Pd2,6ST, an initial titer of 0.98 g/L was achieved. Subsequent combinatorial expression of heterologous enzymes from diverse microbial sources increased the yield to 1.26 g/L. Through systematic optimization of promoter strength and ribosome binding site (RBS) elements, the 6'-SL titer was further elevated to 2.41 g/L. Implementation of structure-guided mutagenesis on both key enzymes resulted in a significant yield improvement to 4.27 g/L. Final engineering strategies included overexpression of rate-limiting enzymes in the CTP biosynthesis pathway. The engineered strain demonstrated a production capacity of 4.92 g/L in shake-flask cultures. We transferred the strain to a 5-L fermenter and achieved a final 6'-SL production of 12.82 g/L after 76 h of cultivation, demonstrating a 13-fold enhancement compared to the initial shake-flask culture.

Variant effect predictor correlation with functional assays is reflective of clinical classification performance

BACKGROUND

Understanding the relationship between protein sequence and function is crucial for accurate classification of missense variants. Variant effect predictors (VEPs) play a vital role in deciphering this complex relationship, yet evaluating their performance remains challenging for several reasons, including data circularity, where the same or related data is used for training and assessment. High-throughput experimental strategies like deep mutational scanning (DMS) offer a promising solution.

RESULTS

In this study, we extend upon our previous benchmarking approach, assessing the performance of 97 VEPs using missense DMS measurements from 36 different human proteins. In addition, a new pairwise, VEP-centric approach mitigates the impact of missing predictions on overall performance comparison. We observe a strong correspondence between VEP performance in DMS-based benchmarks and clinical variant classification, especially for predictors that have not been directly trained on human clinical variants.

CONCLUSIONS

Our results suggest that comparing VEP performance against diverse functional assays represents a reliable strategy for assessing their relative performance in clinical variant classification. However, major challenges in clinical interpretation of VEP scores persist, highlighting the need for further research to fully leverage computational predictors for genetic diagnosis. We also address practical considerations for end users in terms of choice of methodology.

Point mutations of the mitochondrial chaperone TRAP1 affect its functions and pro-neoplastic activity

The mitochondrial chaperone TRAP1 is a key regulator of cellular homeostasis and its activity has important implications in neurodegeneration, ischemia and cancer. Recent evidence has indicated that TRAP1 mutations are involved in several disorders, even though the structural basis for the impact of point mutations on TRAP1 functions has never been studied. By exploiting a modular structure-based framework and molecular dynamics simulations, we investigated the effect of five TRAP1 mutations on its structure and stability. Each mutation differentially impacts long-range interactions, intra and inter-protomer dynamics and ATPase activity. Changes in these parameters influence TRAP1 functions, as revealed by their effects on the activity of the TRAP1 interactor succinate dehydrogenase (SDH). In keeping with this, TRAP1 point mutations affect the growth and migration of aggressive sarcoma cells, and alter sensitivity to a selective TRAP1 inhibitor. Our work provides new insights on the structure-activity relationship of TRAP1, identifying crucial amino acid residues that regulate TRAP1 proteostatic functions and pro-neoplastic activity.

Engineering inhibitory repeat domains of Pin-II type proteinase inhibitors indicate their high structural-functional tolerance to mutagenesis

Plant proteinase inhibitors (PIs) are critical in defending against biotic stress. Most PIs contain an inhibitory repeat domain (IRD), which serves as the functional component, displaying a high degree of sequence and structural conservation. In this study, we examined the structural and functional resilience of IRDs using a combination of computational modeling and experimental validation. We have taken an evolution-based approach to enhance the PIs effectiveness of two previously identified Capsicum annuum IRDs, IRD4 and IRD10. Through in silico site-saturation mutagenesis of IRD4 and IRD10, we identified key sites associated with enhanced PI activity for targeted mutagenesis. Binding energy predictions for a mutant IRD library, tested against target proteases, suggested that positions R11 and N32 in IRD4 and N32 and H33 in IRD10 were promising candidates for further modification to improve inhibitory potential. Subsequent experimental validation revealed that the mutant proteins IRD4_R11K and IRD4_N32S exhibited stronger chymotrypsin inhibition than the wild-type (WT) IRD4. Similarly, the mutants IRD10_N32S and IRD10_H33 N demonstrated improved trypsin inhibition relative to the WT IRD10. These findings indicate that engineered IRD variants can tolerate structural changes while maintaining or enhancing their inhibitory activity against target proteases. Overall, this study demonstrates the potential of engineering PIs to increase their structural and functional resilience, offering new opportunities for biotechnological applications.

Screening and modification of (+)-germacrene A synthase for the production of the anti-tumor drug (-)-β-elemene in engineered Saccharomyces cerevisiae

(-)-β-Elemene is a primary bioactive compound derived from Curcuma wenyujin and has been widely utilized as an anti-tumor agent for various types of cancer. Due to the inefficiency of plant extraction methods for β-elemene, significant efforts have been directed toward the heterogeneous biosynthesis of β-elemene using microbial cell factories. However, there has been less emphasis on the stereochemical configuration of germacrene A and its rearranged product, β-elemene. In this study, we constructed a yeast cell factory to produce (-)-β-elemene by optimizing the mevalonate pathway and screening for germacrene A synthases (GASs) from both plant and microbial sources. Notably, we discovered that the rearranged products of GASs exhibited different conformations, and only (+)-germacrene A produced by plant-derived GASs could rearrange to form (-)-β-elemene. Building on this discovery, we further investigated the catalytic mechanisms of GASs and developed an efficient catalytic gene module for generating (+)-germacrene A. Ultimately, the engineered yeast produced 1152 mg/L of (-)-β-elemene, marking the highest titer reported in yeast to date. Overall, this work highlights the differences in the stereoconformations of catalytic products between plant- and microbial-derived germacrene A synthases and establishes a foundation for the green and efficient production of β-elemene with a specific stereochemical configuration.

ASM variants in the spotlight: A structure-based atlas for unraveling pathogenic mechanisms in lysosomal acid sphingomyelinase

Lysosomal acid sphingomyelinase (ASM), a critical enzyme in lipid metabolism encoded by the SMPD1 gene, plays a crucial role in sphingomyelin hydrolysis in lysosomes. ASM deficiency leads to acid sphingomyelinase deficiency, a rare genetic disorder with diverse clinical manifestations, and the protein can be found mutated in other diseases. We employed a structure-based framework to comprehensively understand the functional implications of ASM variants, integrating pathogenicity predictions with molecular insights derived from a molecular dynamics simulation in a lysosomal membrane environment. Our analysis, encompassing over 400 variants, establishes a structural atlas of missense variants of lysosomal ASM, associating mechanistic indicators with pathogenic potential. Our study highlights variants that influence structural stability or exert local and long-range effects at functional sites. To validate our predictions, we compared them to available experimental data on residual catalytic activity in 135 ASM variants. Notably, our findings also suggest applications of the resulting data for identifying cases suited for enzyme replacement therapy. This comprehensive approach enhances the understanding of ASM variants and provides valuable insights for potential therapeutic interventions.

AlzDiscovery: A computational tool to identify Alzheimer's disease-causing missense mutations using protein structure information

Alzheimer's disease (AD) is one of the most common forms of dementia and neurodegenerative diseases, characterized by the formation of neuritic plaques and neurofibrillary tangles. Many different proteins participate in this complicated pathogenic mechanism, and missense mutations can alter the folding and functions of these proteins, significantly increasing the risk of AD. However, many methods to identify AD-causing variants did not consider the effect of mutations from the perspective of a protein three-dimensional environment. Here, we present a machine learning-based analysis to classify the AD-causing mutations from their benign counterparts in 21 AD-related proteins leveraging both sequence- and structure-based features. Using computational tools to estimate the effect of mutations on protein stability, we first observed a bias of the pathogenic mutations with significant destabilizing effects on family AD-related proteins. Combining this insight, we built a generic predictive model, and improved the performance by tuning the sample weights in the training process. Our final model achieved the performance on area under the receiver operating characteristic curve up to 0.95 in the blind test and 0.70 in an independent clinical validation, outperforming all the state-of-the-art methods. Feature interpretation indicated that the hydrophobic environment and polar interaction contacts were crucial to the decision on pathogenic phenotypes of missense mutations. Finally, we presented a user-friendly web server, AlzDiscovery, for researchers to browse the predicted phenotypes of all possible missense mutations on these 21 AD-related proteins. Our study will be a valuable resource for AD screening and the development of personalized treatment.

Variant Impact Predictor database (VIPdb), version 2: trends from three decades of genetic variant impact predictors

BACKGROUND

Variant interpretation is essential for identifying patients' disease-causing genetic variants amongst the millions detected in their genomes. Hundreds of Variant Impact Predictors (VIPs), also known as Variant Effect Predictors (VEPs), have been developed for this purpose, with a variety of methodologies and goals. To facilitate the exploration of available VIP options, we have created the Variant Impact Predictor database (VIPdb).

RESULTS

The Variant Impact Predictor database (VIPdb) version 2 presents a collection of VIPs developed over the past three decades, summarizing their characteristics, ClinGen calibrated scores, CAGI assessment results, publication details, access information, and citation patterns. We previously summarized 217 VIPs and their features in VIPdb in 2019. Building upon this foundation, we identified and categorized an additional 190 VIPs, resulting in a total of 407 VIPs in VIPdb version 2. The majority of the VIPs have the capacity to predict the impacts of single nucleotide variants and nonsynonymous variants. More VIPs tailored to predict the impacts of insertions and deletions have been developed since the 2010s. In contrast, relatively few VIPs are dedicated to the prediction of splicing, structural, synonymous, and regulatory variants. The increasing rate of citations to VIPs reflects the ongoing growth in their use, and the evolving trends in citations reveal development in the field and individual methods.

CONCLUSIONS

VIPdb version 2 summarizes 407 VIPs and their features, potentially facilitating VIP exploration for various variant interpretation applications. VIPdb is available at  https://genomeinterpretation.org/vipdb.

Flattening the curve-How to get better results with small deep-mutational-scanning datasets

Proteins are used in various biotechnological applications, often requiring the optimization of protein properties by introducing specific amino-acid exchanges. Deep mutational scanning (DMS) is an effective high-throughput method for evaluating the effects of these exchanges on protein function. DMS data can then inform the training of a neural network to predict the impact of mutations. Most approaches use some representation of the protein sequence for training and prediction. As proteins are characterized by complex structures and intricate residue interaction networks, directly providing structural information as input reduces the need to learn these features from the data. We introduce a method for encoding protein structures as stacked 2D contact maps, which capture residue interactions, their evolutionary conservation, and mutation-induced interaction changes. Furthermore, we explored techniques to augment neural network training performance on smaller DMS datasets. To validate our approach, we trained three neural network architectures originally used for image analysis on three DMS datasets, and we compared their performances with networks trained solely on protein sequences. The results confirm the effectiveness of the protein structure encoding in machine learning efforts on DMS data. Using structural representations as direct input to the networks, along with data augmentation and pretraining, significantly reduced demands on training data size and improved prediction performance, especially on smaller datasets, while performance on large datasets was on par with state-of-the-art sequence convolutional neural networks. The methods presented here have the potential to provide the same workflow as DMS without the experimental and financial burden of testing thousands of mutants. Additionally, we present an open-source, user-friendly software tool to make these data analysis techniques accessible, particularly to biotechnology and protein engineering researchers who wish to apply them to their mutagenesis data.

Variant Impact Predictor database (VIPdb), version 2: Trends from 25 years of genetic variant impact predictors

Background

Variant interpretation is essential for identifying patients' disease-causing genetic variants amongst the millions detected in their genomes. Hundreds of Variant Impact Predictors (VIPs), also known as Variant Effect Predictors (VEPs), have been developed for this purpose, with a variety of methodologies and goals. To facilitate the exploration of available VIP options, we have created the Variant Impact Predictor database (VIPdb).

Results

The Variant Impact Predictor database (VIPdb) version 2 presents a collection of VIPs developed over the past 25 years, summarizing their characteristics, ClinGen calibrated scores, CAGI assessment results, publication details, access information, and citation patterns. We previously summarized 217 VIPs and their features in VIPdb in 2019. Building upon this foundation, we identified and categorized an additional 186 VIPs, resulting in a total of 403 VIPs in VIPdb version 2. The majority of the VIPs have the capacity to predict the impacts of single nucleotide variants and nonsynonymous variants. More VIPs tailored to predict the impacts of insertions and deletions have been developed since the 2010s. In contrast, relatively few VIPs are dedicated to the prediction of splicing, structural, synonymous, and regulatory variants. The increasing rate of citations to VIPs reflects the ongoing growth in their use, and the evolving trends in citations reveal development in the field and individual methods.

Conclusions

VIPdb version 2 summarizes 403 VIPs and their features, potentially facilitating VIP exploration for various variant interpretation applications.

Availability

VIPdb version 2 is available at https://genomeinterpretation.org/vipdb.

Integrated multiplexed assays of variant effect reveal determinants of catechol-O-methyltransferase gene expression

Multiplexed assays of variant effect are powerful methods to profile the consequences of rare variants on gene expression and organismal fitness. Yet, few studies have integrated several multiplexed assays to map variant effects on gene expression in coding sequences. Here, we pioneered a multiplexed assay based on polysome profiling to measure variant effects on translation at scale, uncovering single-nucleotide variants that increase or decrease ribosome load. By combining high-throughput ribosome load data with multiplexed mRNA and protein abundance readouts, we mapped the cis-regulatory landscape of thousands of catechol-O-methyltransferase (COMT) variants from RNA to protein and found numerous coding variants that alter COMT expression. Finally, we trained machine learning models to map signatures of variant effects on COMT gene expression and uncovered both directional and divergent impacts across expression layers. Our analyses reveal expression phenotypes for thousands of variants in COMT and highlight variant effects on both single and multiple layers of expression. Our findings prompt future studies that integrate several multiplexed assays for the readout of gene expression.

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.

Analysis of a Novel Peptide That Is Capable of Inhibiting the Enzymatic Activity of the Protein Kinase A Catalytic Subunit-Like Protein from Trypanosoma equiperdum

A 26-residue peptide possessing the αN-helix motif of the protein kinase A (PKA) regulatory subunit-like proteins from the Trypanozoom subgenera (VAP26, sequence = VAPYFEKSEDETALILKLLTYNVLFS), was shown to inhibit the enzymatic activity of the Trypanosoma equiperdum PKA catalytic subunit-like protein, in a similar manner that the mammalian heat-stable soluble PKA inhibitor known as PKI. However, VAP26 does not contain the PKI inhibitory sequence. Bioinformatics analyzes of the αN-helix motif from various Trypanozoon PKA regulatory subunit-like proteins suggested that the sequence could form favorable peptide-protein interactions of hydrophobic nature with the PKA catalytic subunit-like protein, which possibly may represent an alternative PKA inhibitory mechanism. The sequence of the αN-helix motif of the Trypanozoon proteins was shown to be highly homologous but significantly divergent from the corresponding αN-helix motifs of their Leishmania and mammalian counterparts. This sequence divergence contrasted with the proposed secondary structure of the αN-helix motif, which appeared conserved in every analyzed regulatory subunit-like protein. In silico mutation experiments at positions I234, L238 and F244 of the αN-helix motif from the Trypanozoon proteins destabilized both the specific motif and the protein. On the contrary, mutations at positions T239 and Y240 stabilized the motif and the protein. These results suggested that the αN-helix motif from the Trypanozoon proteins probably possessed a different evolutionary path than their Leishmania and mammalian counterparts. Moreover, finding stabilizing mutations indicated that new inhibitory peptides may be designed based on the αN-helix motif from the Trypanozoon PKA regulatory subunit-like proteins.

Structural and Catalytic Insight into the Unique Pentacyclic Triterpene Synthase TwOSC

The oxidosqualene cyclase (OSC) catalyzed cyclization of the linear substrate (3S)-2,3-oxidosqualene to form diverse pentacyclic triterpenoid (PT) skeletons is one of the most complex reactions in nature. Friedelin has a unique PT skeleton involving a fascinating nine-step cation shuttle run (CSR) cascade rearrangement reaction, in which the carbocation formed at C2 moves to the other side of the skeleton, runs back to C3 to yield a friedelin cation, which is finally deprotonated. However, as crystal structure data of plant OSCs are lacking, it remains unknown why the CSR cascade reactions occur in friedelin biosynthesis, as does the exact catalytic mechanism of the CSR. In this study, we determined the first cryogenic electron microscopy structure of a plant OSC, friedelin synthase, from Tripterygium wilfordii Hook. f (TwOSC). We also performed quantum mechanics/molecular mechanics simulations to reveal the energy profile for the CSR cascade reaction and identify key residues crucial for PT skeleton formation. Furthermore, we semirationally designed two TwOSC mutants, which significantly improved the yields of friedelin and β-amyrin, respectively.

Characterization on the oncogenic effect of the missense mutations of p53 via machine learning

Dysfunctions caused by missense mutations in the tumour suppressor p53 have been extensively shown to be a leading driver of many cancers. Unfortunately, it is time-consuming and labour-intensive to experimentally elucidate the effects of all possible missense variants. Recent works presented a comprehensive dataset and machine learning model to predict the functional outcome of mutations in p53. Despite the well-established dataset and precise predictions, this tool was trained on a complicated model with limited predictions on p53 mutations. In this work, we first used computational biophysical tools to investigate the functional consequences of missense mutations in p53, informing a bias of deleterious mutations with destabilizing effects. Combining these insights with experimental assays, we present two interpretable machine learning models leveraging both experimental assays and in silico biophysical measurements to accurately predict the functional consequences on p53 and validate their robustness on clinical data. Our final model based on nine features obtained comparable predictive performance with the state-of-the-art p53 specific method and outperformed other generalized, widely used predictors. Interpreting our models revealed that information on residue p53 activity, polar atom distances and changes in p53 stability were instrumental in the decisions, consistent with a bias of the properties of deleterious mutations. Our predictions have been computed for all possible missense mutations in p53, offering clinical diagnostic utility, which is crucial for patient monitoring and the development of personalized cancer treatment.

Integrated multiplexed assays of variant effect reveal cis-regulatory determinants of catechol-O-methyltransferase gene expression

Multiplexed assays of variant effect are powerful methods to profile the consequences of rare variants on gene expression and organismal fitness. Yet, few studies have integrated several multiplexed assays to map variant effects on gene expression in coding sequences. Here, we pioneered a multiplexed assay based on polysome profiling to measure variant effects on translation at scale, uncovering single-nucleotide variants that increase and decrease ribosome load. By combining high-throughput ribosome load data with multiplexed mRNA and protein abundance readouts, we mapped the cis-regulatory landscape of thousands of catechol-O-methyltransferase (COMT) variants from RNA to protein and found numerous coding variants that alter COMT expression. Finally, we trained machine learning models to map signatures of variant effects on COMT gene expression and uncovered both directional and divergent impacts across expression layers. Our analyses reveal expression phenotypes for thousands of variants in COMT and highlight variant effects on both single and multiple layers of expression. Our findings prompt future studies that integrate several multiplexed assays for the readout of gene expression.

Systematically testing human HMBS missense variants to reveal mechanism and pathogenic variation

Defects in hydroxymethylbilane synthase (HMBS) can cause acute intermittent porphyria (AIP), an acute neurological disease. Although sequencing-based diagnosis can be definitive, ∼⅓ of clinical HMBS variants are missense variants, and most clinically reported HMBS missense variants are designated as "variants of uncertain significance" (VUSs). Using saturation mutagenesis, en masse selection, and sequencing, we applied a multiplexed validated assay to both the erythroid-specific and ubiquitous isoforms of HMBS, obtaining confident functional impact scores for >84% of all possible amino acid substitutions. The resulting variant effect maps generally agreed with biochemical expectations and provide further evidence that HMBS can function as a monomer. Additionally, the maps implicated specific residues as having roles in active site dynamics, which was further supported by molecular dynamics simulations. Most importantly, these maps can help discriminate pathogenic from benign HMBS variants, proactively providing evidence even for yet-to-be-observed clinical missense variants.

Accurate proteome-wide missense variant effect prediction with AlphaMissense

The vast majority of missense variants observed in the human genome are of unknown clinical significance. We present AlphaMissense, an adaptation of AlphaFold fine-tuned on human and primate variant population frequency databases to predict missense variant pathogenicity. By combining structural context and evolutionary conservation, our model achieves state-of-the-art results across a wide range of genetic and experimental benchmarks, all without explicitly training on such data. The average pathogenicity score of genes is also predictive for their cell essentiality, capable of identifying short essential genes that existing statistical approaches are underpowered to detect. As a resource to the community, we provide a database of predictions for all possible human single amino acid substitutions and classify 89% of missense variants as either likely benign or likely pathogenic.

Identifying the Molecular Drivers of Pathogenic Aldehyde Dehydrogenase Missense Mutations in Cancer and Non-Cancer Diseases

Human aldehyde dehydrogenases (ALDHs) comprising 19 isoenzymes play a vital role on both endogenous and exogenous aldehyde metabolism. This NAD(P)-dependent catalytic process relies on the intact structural and functional activity of the cofactor binding, substrate interaction, and the oligomerization of ALDHs. Disruptions on the activity of ALDHs, however, could result in the accumulation of cytotoxic aldehydes, which have been linked with a wide range of diseases, including both cancers as well as neurological and developmental disorders. In our previous works, we have successfully characterised the structure-function relationships of the missense variants of other proteins. We, therefore, applied a similar analysis pipeline to identify potential molecular drivers of pathogenic ALDH missense mutations. Variants data were first carefully curated and labelled as cancer-risk, non-cancer diseases, and benign. We then leveraged various computational biophysical methods to describe the changes caused by missense mutations, informing a bias of detrimental mutations with destabilising effects. Cooperating with these insights, several machine learning approaches were further utilised to investigate the combination of features, revealing the necessity of the conservation of ALDHs. Our work aims to provide important biological perspectives on pathogenic consequences of missense mutations of ALDHs, which could be invaluable resources in the development of cancer treatment.

Integrating deep mutational scanning and low-throughput mutagenesis data to predict the impact of amino acid variants

BACKGROUND

Evaluating the impact of amino acid variants has been a critical challenge for studying protein function and interpreting genomic data. High-throughput experimental methods like deep mutational scanning (DMS) can measure the effect of large numbers of variants in a target protein, but because DMS studies have not been performed on all proteins, researchers also model DMS data computationally to estimate variant impacts by predictors.

RESULTS

In this study, we extended a linear regression-based predictor to explore whether incorporating data from alanine scanning (AS), a widely used low-throughput mutagenesis method, would improve prediction results. To evaluate our model, we collected 146 AS datasets, mapping to 54 DMS datasets across 22 distinct proteins.

CONCLUSIONS

We show that improved model performance depends on the compatibility of the DMS and AS assays, and the scale of improvement is closely related to the correlation between DMS and AS results.

Genome interpretation using in silico predictors of variant impact

Estimating the effects of variants found in disease driver genes opens the door to personalized therapeutic opportunities. Clinical associations and laboratory experiments can only characterize a tiny fraction of all the available variants, leaving the majority as variants of unknown significance (VUS). In silico methods bridge this gap by providing instant estimates on a large scale, most often based on the numerous genetic differences between species. Despite concerns that these methods may lack reliability in individual subjects, their numerous practical applications over cohorts suggest they are already helpful and have a role to play in genome interpretation when used at the proper scale and context. In this review, we aim to gain insights into the training and validation of these variant effect predicting methods and illustrate representative types of experimental and clinical applications. Objective performance assessments using various datasets that are not yet published indicate the strengths and limitations of each method. These show that cautious use of in silico variant impact predictors is essential for addressing genome interpretation challenges.

How Functional Genomics Can Keep Pace With VUS Identification

Over the last two decades, an exponentially expanding number of genetic variants have been identified associated with inherited cardiac conditions. These tremendous gains also present challenges in deciphering the clinical relevance of unclassified variants or variants of uncertain significance (VUS). This review provides an overview of the advancements (and challenges) in functional and computational approaches to characterize variants and help keep pace with VUS identification related to inherited heart diseases.

A guide to machine learning for biologists

The expanding scale and inherent complexity of biological data have encouraged a growing use of machine learning in biology to build informative and predictive models of the underlying biological processes. All machine learning techniques fit models to data; however, the specific methods are quite varied and can at first glance seem bewildering. In this Review, we aim to provide readers with a gentle introduction to a few key machine learning techniques, including the most recently developed and widely used techniques involving deep neural networks. We describe how different techniques may be suited to specific types of biological data, and also discuss some best practices and points to consider when one is embarking on experiments involving machine learning. Some emerging directions in machine learning methodology are also discussed.

Intelligent host engineering for metabolic flux optimisation in biotechnology

Optimising the function of a protein of length N amino acids by directed evolution involves navigating a 'search space' of possible sequences of some 20N. Optimising the expression levels of P proteins that materially affect host performance, each of which might also take 20 (logarithmically spaced) values, implies a similar search space of 20P. In this combinatorial sense, then, the problems of directed protein evolution and of host engineering are broadly equivalent. In practice, however, they have different means for avoiding the inevitable difficulties of implementation. The spare capacity exhibited in metabolic networks implies that host engineering may admit substantial increases in flux to targets of interest. Thus, we rehearse the relevant issues for those wishing to understand and exploit those modern genome-wide host engineering tools and thinking that have been designed and developed to optimise fluxes towards desirable products in biotechnological processes, with a focus on microbial systems. The aim throughput is 'making such biology predictable'. Strategies have been aimed at both transcription and translation, especially for regulatory processes that can affect multiple targets. However, because there is a limit on how much protein a cell can produce, increasing kcat in selected targets may be a better strategy than increasing protein expression levels for optimal host engineering.

Harnessing the Genetic Plasticity of Porcine Circovirus Type 2 to Target Suicidal Replication

Porcine circovirus type 2 (PCV2), the causative agent of a wasting disease in weanling piglets, has periodically evolved into several new subtypes since its discovery, indicating that the efficacy of current vaccines can be improved. Although a DNA virus, the mutation rates of PCV2 resemble RNA viruses. The hypothesis that recoding of selected serine and leucine codons in the PCV2b capsid gene could result in stop codons due to mutations occurring during viral replication and thus result in rapid attenuation was tested. Vaccination of weanling pigs with the suicidal vaccine constructs elicited strong virus-neutralizing antibody responses. Vaccination prevented lesions, body-weight loss, and viral replication on challenge with a heterologous PCV2d strain. The suicidal PCV2 vaccine construct was not detectable in the sera of vaccinated pigs at 14 days post-vaccination, indicating that the attenuated vaccine was very safe. Exposure of the modified virus to immune selection pressure with sub-neutralizing levels of antibodies resulted in 5 of the 22 target codons mutating to a stop signal. Thus, the described approach for the rapid attenuation of PCV2 was both effective and safe. It can be readily adapted to newly emerging viruses with high mutation rates to meet the current need for improved platforms for rapid-response vaccines.