A Single Synonymous Variant (c.354G>A [p.P118P]) in ADAMTS13 Confers Enhanced Specific Activity

How synonymous mutations alter enzyme structure and function over long timescales

The specific activity of enzymes can be altered over long timescales in cells by synonymous mutations that alter a messenger RNA molecule's sequence but not the encoded protein's primary structure. How this happens at the molecular level is unknown. Here, we use multiscale modelling of three Escherichia coli enzymes (type III chloramphenicol acetyltransferase, D-alanine-D-alanine ligase B and dihydrofolate reductase) to understand experimentally measured changes in specific activity due to synonymous mutations. The modelling involves coarse-grained simulations of protein synthesis and post-translational behaviour, all-atom simulations to test robustness and quantum mechanics/molecular mechanics calculations to characterize enzymatic function. We show that changes in codon translation rates induced by synonymous mutations cause shifts in co-translational and post-translational folding pathways that kinetically partition molecules into subpopulations that very slowly interconvert to the native, functional state. Structurally, these states resemble the native state, with localized misfolding near the active sites of the enzymes. These long-lived states exhibit reduced catalytic activity, as shown by their increased activation energies for the reactions they catalyse.

frDSM: An Ensemble Predictor With Effective Feature Representation for Deleterious Synonymous Mutation in Human Genome

With the discovery of causality between synonymous mutations and diseases, it has become increasingly important to identify deleterious synonymous mutations for better understanding of their functional mechanisms. Although several machine learning methods have been proposed to solve the task, an effective feature representation method that can make use of the inner difference and relevance between deleterious and benign synonymous mutations is still challenging considering the vast number of synonymous mutations in human genome. In this work, we developed a robust and accurate predictor called frDSM for deleterious synonymous mutation prediction using logistic regression. More specifically, we introduced an effective feature representation learning method which exploits multiple feature descriptors from different perspectives including functional scores obtained from previously computational methods, evolutionary conservation, splicing and sequence feature descriptors, and these features descriptors were input into the 76 XGBoost classifiers to obtain the predictive probabilities values. These probabilities were concatenated to generate the 76-dimension new feature vector, and feature selection method was used to remove redundant and irrelevant features. Experimental results show that frDSM enables robust and accurate prediction than the competing prediction methods with 31 optimal features, which demonstrated the effectiveness of the feature representation learning method. frDSM is freely available at http://frdsm.xialab.info.

The High Mutational Sensitivity of ccdA Antitoxin Is Linked to Codon Optimality

Deep mutational scanning studies suggest that synonymous mutations are typically silent and that most exposed, nonactive-site residues are tolerant to mutations. Here, we show that the ccdA antitoxin component of the Escherichia coli ccdAB toxin-antitoxin system is unusually sensitive to mutations when studied in the operonic context. A large fraction (∼80%) of single-codon mutations, including many synonymous mutations in the ccdA gene shows inactive phenotype, but they retain native-like binding affinity towards cognate toxin, CcdB. Therefore, the observed phenotypic effects are largely not due to alterations in protein structure/stability, consistent with a large region of CcdA being intrinsically disordered. E. coli codon preference and strength of ribosome-binding associated with translation of downstream ccdB gene are found to be major contributors of the observed ccdA mutant phenotypes. In select cases, proteomics studies reveal altered ratios of CcdA:CcdB protein levels in vivo, suggesting that the ccdA mutations likely alter relative translation efficiencies of the two genes in the operon. We extend these results by studying single-site synonymous mutations that lead to loss of function phenotypes in the relBE operon upon introduction of rarer codons. Thus, in their operonic context, genes are likely to be more sensitive to both synonymous and nonsynonymous point mutations than inferred previously.

When a Synonymous Variant Is Nonsynonymous

Term synonymous variation is widely used, but frequently in a wrong or misleading meaning and context. Twenty three point eight % of possible nucleotide substitution types in the universal genetic code are for synonymous amino acid changes, but when these variants have a phenotype and functional effect, they are very seldom synonymous. Such variants may manifest changes at DNA, RNA and/or protein levels. Large numbers of variations are erroneously annotated as synonymous, which causes problems e.g., in clinical genetics and diagnosis of diseases. To facilitate precise communication, novel systematics and nomenclature are introduced for variants that when looking only at the genetic code seem like synonymous, but which have phenotypes. A new term, unsense variant is defined as a substitution in the mRNA coding region that affects gene expression and protein production without introducing a stop codon in the variation site. Such variants are common and need to be correctly annotated. Proper naming and annotation are important also to increase awareness of these variants and their consequences.

Structural, functional, and immunogenicity implications of F9 gene recoding

Hemophilia B is a blood clotting disorder caused by deficient activity of coagulation factor IX (FIX). Multiple recombinant FIX proteins are currently approved to treat hemophilia B, and several gene therapy products are currently being developed. Codon optimization is a frequently used technique in the pharmaceutical industry to improve recombinant protein expression by recoding a coding sequence using multiple synonymous codon substitutions. The underlying assumption of this gene recoding is that synonymous substitutions do not alter protein characteristics because the primary sequence of the protein remains unchanged. However, a critical body of evidence shows that synonymous variants can affect cotranslational folding and protein function. Gene recoding could potentially alter the structure, function, and in vivo immunogenicity of recoded therapeutic proteins. Here, we evaluated multiple recoded variants of F9 designed to further explore the effects of codon usage bias on protein properties. The detailed evaluation of these constructs showed altered conformations, and assessment of translation kinetics by ribosome profiling revealed differences in local translation kinetics. Assessment of wild-type and recoded constructs using a major histocompatibility complex (MHC)-associated peptide proteomics assay showed distinct presentation of FIX-derived peptides bound to MHC class II molecules, suggesting that despite identical amino acid sequence, recoded proteins could exhibit different immunogenicity risks. Posttranslational modification analysis indicated that overexpression from gene recoding results in suboptimal posttranslational processing. Overall, our results highlight potential functional and immunogenicity concerns associated with gene-recoded F9 products. These findings have general applicability and implications for other gene-recoded recombinant proteins.

Research for type 2 diabetes mellitus in endemic arsenism areas in central China: role of low level of arsenic exposure and KEAP1 rs11545829 polymorphism

Type 2 diabetes mellitus (T2DM) is one of the major public health problems worldwide; both genetic and environmental factors are its risk factors. Arsenic, an environmental pollutant, might be a risk factor for T2DM, but the association of low-to-moderate level arsenic exposure with the risk of T2DM is still inconsistent. Single nucleotide polymorphisms (SNPs) can affect the development of T2DM, but the study on KEAP1 rs11545829 (G>A) SNP is few. In this paper, we explored the effect of KEAP1 rs11545829 (G>A) SNP and low-to-moderate level arsenic exposure on risk of T2DM in a cross-sectional case-control study conducted in Shanxi, China. Total of 938 participants, including 318 T2DM cases and 618 controls, were enrolled. Blood glycosylated haemoglobin (HbA1c) was detected by Automatic Biochemical Analyzer, and participants with HbA1c≧6.5% were diagnosed as T2DM. Urinary total arsenic (tAs, mg/L), as the indicator of arsenic exposure, was detected by liquid chromatography-atomic fluorescence spectrometry (LC-AFS). Genomic DNA was extracted and the genotypes of KEAP1 rs11545829 SNP were examined by multiplex polymerase chain reaction (PCR). The urinary tAs concentration in recruited participants was 0.075 (0.03-0.15) mg/L, and was associated with an increased risk of T2DM (OR = 8.45, 95% CI 2.63-27.17); rs11545829 mutation homozygote AA genotype had a protective effect on risk of T2DM (OR = 0.42, 95 % CI 0.25-0.73). Although this protective effect of AA genotype was found in participants with higher urinary tAs level (>0.032 mg/L) (OR = 0.48, 95% CI 0.26-0.86), there was no interaction effect for arsenic exposure and rs11545829 SNP on risk of T2DM. In addition, BMI modified the association between rs11545829 SNP and the risk of T2DM (RERI = -1.11, 95% CI -2.18-0.04). The present study suggest that low-to-moderate level arsenic exposure may be a risk factor, while KEAP1 rs11545829 SNP mutation homozygote AA genotype may be a protective factor for risk of T2DM, especially for T2DM patients with urinary tAs level>0.032 mg/L.

The Single Nucleotide Polymorphisms of AP1S1 are Associated with Risk of Esophageal Squamous Cell Carcinoma in Chinese Population

Background

The σ1A subunit of the adaptor protein 1 (AP1S1) participates in various intracellular transport pathways, especially the maintenance of copper homeostasis, which is pivotal in carcinogenesis. It is therefore rational to presume that AP1S1 might also be involved in carcinogenesis. In this hospital-based case-control study, we investigated the genetic susceptibility to ESCC in relation to SNPs of AP1S1 among Chinese population.

Methods

A database containing a total of 1303 controls and 1043 ESCC patients were retrospectively studied. The AP1S1 SNPs were analyzed based on ligation detection reaction (LDR) method. Then, the relationship between ESCC and SNPs of AP1S1 was determined with a significant crude P<0.05. Then the logistic regression analysis was used for the calculation for adjusted P in the demographic stratification comparison if a significant difference was observed in the previous step.

Results

AP1S1 rs77387752 C>T genotype TT was an independent risk factor for ESCC, while rs4729666 C>T genotype TC and rs35208462 C>T genotype TC were associated with a lower risk for ESCC, especially in co-dominant model and allelic test for younger, male subjects who are not alcohol-drinkers nor cigarette smokers.

Conclusion

AP1S1 rs77387752, rs4729666 and rs35208462 polymorphisms are associated with susceptibility to ESCC in Chinese individuals. AP1S1 SNPs may exert an important role in esophageal carcinogenesis and could serve as potential diagnostic biomarkers.

An Optimized Purification Design for Extracting Active ADAMTS13 from Conditioned Media

ADAMTS13 is a hemostatic enzyme that breaks down pro-thrombotic ultra-large multimers of von Willebrand factor (VWF). The deficiency of ADAMTS13 increases VWF-mediated thrombogenic potential and may lead to thrombotic thrombocytopenic purpura (TTP). Recently, clinical studies have shown the development of acquired TTP after COVID-19 infection and a correlation between low ADAMTS13 plasma levels and increased mortality. As a result, investigating ADAMTS13 as a potential recombinant therapeutic is of broad interest in the field of hematology. ADAMTS13 is considered challenging to purify in its biologically active state. Current purification methods utilize immobilized metal ions, which can interfere with ADAMTS13 metalloprotease activity. For this reason, we optimized an alternative strategy to isolate milligram quantities of highly active recombinant ADAMTS13 (rADAMTS13) from conditioned media after exogenous expression in human cell line, HEK293. HEK293 cells stably expressing C-terminal V5-His-tagged ADAMTS13 were grown in two parallel systems, culture bottles and flasks, for identifying an optimal cultivation strategy. Subsequently, we employed anion exchange followed by anti-V5-tag affinity chromatography to purify rADAMTS13, and extracted rADAMTS13 of high specific activity while preserving its native post-translational modifications. In addition, this process has been optimized and scaled up to produce active rADAMTS13 at levels sufficient for laboratory-scale structural, enzymatic, and biochemical studies.

A Code Within a Code: How Codons Fine-Tune Protein Folding in the Cell

The genetic code sets the correspondence between the sequence of a given nucleotide triplet in an mRNA molecule, called a codon, and the amino acid that is added to the growing polypeptide chain during protein synthesis. With four bases (A, G, U, and C), there are 64 possible triplet codons: 61 sense codons (encoding amino acids) and 3 nonsense codons (so-called, stop codons that define termination of translation). In most organisms, there are 20 common/standard amino acids used in protein synthesis; thus, the genetic code is redundant with most amino acids (with the exception of Met and Trp) are being encoded by more than one (synonymous) codon. Synonymous codons were initially presumed to have entirely equivalent functions, however, the finding that synonymous codons are not present at equal frequencies in mRNA suggested that the specific codon choice might have functional implications beyond coding for amino acid. Observation of nonequivalent use of codons in mRNAs implied a possibility of the existence of auxiliary information in the genetic code. Indeed, it has been found that genetic code contains several layers of such additional information and that synonymous codons are strategically placed within mRNAs to ensure a particular translation kinetics facilitating and fine-tuning co-translational protein folding in the cell via step-wise/sequential structuring of distinct regions of the polypeptide chain emerging from the ribosome at different points in time. This review summarizes key findings in the field that have identified the role of synonymous codons and their usage in protein folding in the cell.

Synonymous but Not Silent: The Codon Usage Code for Gene Expression and Protein Folding

Codon usage bias, the preference for certain synonymous codons, is found in all genomes. Although synonymous mutations were previously thought to be silent, a large body of evidence has demonstrated that codon usage can play major roles in determining gene expression levels and protein structures. Codon usage influences translation elongation speed and regulates translation efficiency and accuracy. Adaptation of codon usage to tRNA expression determines the proteome landscape. In addition, codon usage biases result in nonuniform ribosome decoding rates on mRNAs, which in turn influence the cotranslational protein folding process that is critical for protein function in diverse biological processes. Conserved genome-wide correlations have also been found between codon usage and protein structures. Furthermore, codon usage is a major determinant of mRNA levels through translation-dependent effects on mRNA decay and translation-independent effects on transcriptional and posttranscriptional processes. Here, we discuss the multifaceted roles and mechanisms of codon usage in different gene regulatory processes.

Nonoptimal Codon Usage Is Critical for Protein Structure and Function of the Master General Amino Acid Control Regulator CPC-1

Under amino acid starvation conditions, eukaryotic organisms activate a general amino acid control response. In Neurospora crassa, Cross Pathway Control Protein 1 (CPC-1), the ortholog of the Saccharomyces cerevisiae bZIP transcription factor GCN4, functions as the master regulator of the general amino acid control response. Codon usage biases are a universal feature of eukaryotic genomes and are critical for regulation of gene expression. Although codon usage has also been implicated in the regulation of protein structure and function, genetic evidence supporting this conclusion is very limited. Here, we show that Neurospora cpc-1 has a nonoptimal NNU-rich codon usage profile that contrasts with the strong NNC codon preference in the genome. Although substitution of the cpc-1 NNU codons with synonymous NNC codons elevated CPC-1 expression in Neurospora, it altered the CPC-1 degradation rate and abolished its amino acid starvation-induced protein stabilization. The codon-manipulated CPC-1 protein also exhibited different sensitivity to limited protease digestion. Furthermore, CPC-1 functions in rescuing the cell growth of the cpc-1 deletion mutant and activation of the expression of its target genes were impaired by the synonymous codon changes. Together, these results reveal the critical role of codon usage in regulation of CPC-1 expression and function and establish a genetic example of the importance of codon usage in protein folding.IMPORTANCE The general amino acid control response is critical for adaptation of organisms to amino acid starvation conditions. The preference to use certain synonymous codons is a universal feature of all genomes. Synonymous codon changes were previously thought to be silent mutations. In this study, we showed that the Neurospora cpc-1 gene has an unusual codon usage profile compared to other genes in the genome. We found that codon optimization of the cpc-1 gene without changing its amino acid sequence resulted in elevated CPC-1 expression, an altered protein degradation rate, and impaired protein functions due to changes in protein structure. Together, these results reveal the critical role of synonymous codon usage in regulation of CPC-1 expression and function and establish a genetic example of the importance of codon usage in protein structure.

A code within the genetic code: codon usage regulates co-translational protein folding

The genetic code is degenerate, and most amino acids are encoded by two to six synonymous codons. Codon usage bias, the preference for certain synonymous codons, is a universal feature of all genomes examined. Synonymous codon mutations were previously thought to be silent; however, a growing body evidence now shows that codon usage regulates protein structure and gene expression through effects on co-translational protein folding, translation efficiency and accuracy, mRNA stability, and transcription. Codon usage regulates the speed of translation elongation, resulting in non-uniform ribosome decoding rates on mRNAs during translation that is adapted to co-translational protein folding process. Biochemical and genetic evidence demonstrate that codon usage plays an important role in regulating protein folding and function in both prokaryotic and eukaryotic organisms. Certain protein structural types are more sensitive than others to the effects of codon usage on protein folding, and predicted intrinsically disordered domains are more prone to misfolding caused by codon usage changes than other domain types. Bioinformatic analyses revealed that gene codon usage correlates with different protein structures in diverse organisms, indicating the existence of a codon usage code for co-translational protein folding. This review focuses on recent literature on the role and mechanism of codon usage in regulating translation kinetics and co-translational protein folding.

In silico features of ADAMTS13 contributing to plasmatic ADAMTS13 levels in neonates with congenital heart disease

INTRODUCTION

Risk factors contributing to heightened thrombosis in pediatric congenital heart disease (CHD) patients are not fully understood. Among the neonatal CHD population, those presenting with single ventricular physiology are at the highest risk for perioperative thrombosis. The von Willebrand factor and ADAMTS13 interactions have emerged as causative risk factors for pediatric stroke and could contribute to heightened thrombosis in CHD neonates.

METHODS

This study investigates a cohort of children with single ventricle physiology and undergoing cardiac surgery, during which some patients developed thrombosis. In this cohort, we analyzed the relationship of several molecular features of ADAMTS13 with the plasma and activity levels in patients at risk of thrombosis. Additionally, in light of the natural antithrombotic activity of ADAMTS13, we have sequenced the ADAMTS13 gene for each patient and evaluated the role of genetic variants in determining the plasma ADAMTS13 levels using a series of in silico tools including Hidden Markov Models, EVmutation, and Rosetta.

RESULTS

Lower ADAMTS13 levels were found in patients that developed thrombosis. A novel in silico analysis to assess haplotype effect of co-occurring variants identified alterations in relative surface area and solvation energy as important contributors. Our analysis suggested that beneficial or deleterious effect of a variant can be reasonably predicted by comprehensive analysis of in silico assessment and in vitro and/or in vivo data.

CONCLUSION

Findings from this study add to our understanding the role of genetic features of ADAMTS13 in patients at high risk of thrombosis related to an imbalanced relation between VWF and ADAMTS13.

Ribosome profiling of HEK293T cells overexpressing codon optimized coagulation factor IX

Ribosome profiling provides the opportunity to evaluate translation kinetics at codon level resolution. Here, we describe ribosome profiling data, generated from two HEK293T cell lines. The ribosome profiling data are composed of Ribo-seq (mRNA sequencing data from ribosome protected fragments) and RNA-seq data (total RNA sequencing). The two HEK293T cell lines each express a version of the F9 gene, both of which are translated into identical proteins in terms of their amino acid sequences. However, these F9 genes vary drastically in their codon usage and predicted mRNA structure. We also provide the pipeline that we used to analyze the data. Further analyzing this dataset holds great potential as it can be used i) to unveil insights into the composition and regulation of the transcriptome, ii) for comparison with other ribosome profiling datasets, iii) to measure the rate of protein synthesis across the proteome and identify differences in elongation rates, iv) to discover previously unidentified translation of peptides, v) to explore the effects of codon usage or codon context in translational kinetics and vi) to investigate cotranslational folding. Importantly, a unique feature of this dataset, compared to other available ribosome profiling data, is the presence of the F9 gene in two very distinct coding sequences.