Dissecting Protein Function: An Efficient Protocol for Identifying Separation-of-Function Mutations That Encode Structurally Stable Proteins

Mutations that confer the loss of a single biochemical property (separation-of-function mutations) can often uncover a previously unknown role for a protein in a particular biological process. However, most mutations are identified based on loss-of-function phenotypes, which cannot differentiate between separation-of-function alleles vs. mutations that encode unstable/unfolded proteins. An alternative approach is to use overexpression dominant-negative (ODN) phenotypes to identify mutant proteins that disrupt function in an otherwise wild-type strain when overexpressed. This is based on the assumption that such mutant proteins retain an overall structure that is comparable to that of the wild-type protein and are able to compete with the endogenous protein (Herskowitz 1987). To test this, the in vivo phenotypes of mutations in the Est3 telomerase subunit from Saccharomyces cerevisiae were compared with the in vitro secondary structure of these mutant proteins as analyzed by circular-dichroism spectroscopy, which demonstrates that ODN is a more sensitive assessment of protein stability than the commonly used method of monitoring protein levels from extracts. Reverse mutagenesis of EST3, which targeted different categories of amino acids, also showed that mutating highly conserved charged residues to the oppositely charged amino acid had an increased likelihood of generating a severely defective est3⁻ mutation, which nevertheless encoded a structurally stable protein. These results suggest that charge-swap mutagenesis directed at a limited subset of highly conserved charged residues, combined with ODN screening to eliminate partially unfolded proteins, may provide a widely applicable and efficient strategy for generating separation-of-function mutations.     

Understanding the functional difference between growth arrest-specific protein 6 and protein S: an evolutionary approach

Although protein S (PROS1) and growth arrest-specific protein 6 (GAS6) proteins are homologous with a high degree of structural similarity, they are functionally different. The objectives of this study were to identify the evolutionary origins from which these functional differences arose. Bioinformatics methods were used to estimate the evolutionary divergence time and to detect the amino acid residues under functional divergence between GAS6 and PROS1. The properties of these residues were analysed in the light of their three-dimensional structures, such as their stability effects, the identification of electrostatic patches and the identification potential protein–protein interaction. The divergence between GAS6 and PROS1 probably occurred during the whole-genome duplications in vertebrates. A total of 78 amino acid sites were identified to be under functional divergence. One of these sites, Asn463, is involved in N-glycosylation in GAS6, but is mutated in PROS1, preventing this post-translational modification. Sites experiencing functional divergence tend to express a greater diversity of stabilizing/destabilizing effects than sites that do not experience such functional divergence. Three electrostatic patches in the LG1/LG2 domains were found to differ between GAS6 and PROS1. Finally, a surface responsible for protein–protein interactions was identified. These results may help researchers to analyse disease-causing mutations in the light of evolutionary and structural constraints, and link genetic pathology to clinical phenotypes.     

New insights into the melanophilin (MLPH) gene controlling coat color phenotypes in American mink

The mutation causing the Silverblue color type (pp) is one of the most used recessive mutations within American mink (Neovison vison) fur farming, since it is involved in some of the popular color types such as Violet and Saphire which originate from a combination of recessive mutations. In the present study, the genomic and mRNA sequences of the melanophilin (MLPH) gene were studied in Violet, Silverblue and wild-type (wt) mink animals. Although breeding schemes and previous literature indicates that the Violet (aammpp) phenotype is a triple recessive color type involving the same locus as the Silverblue (pp) color type, our findings indicate different genotypes at the MLPH locus. Upon comparison at genomic level, we identified two deletions of the entire intron 7 and of the 5′ end of intron 8 in the sequence of the Silverblue MLPH gene. When investigating the mRNA, the Silverblue animals completely lack exon 8, which encodes 65 residues, of which 47 define the Myosin Va (MYO5A) binding domain. This may cause the incorrect anchoring of the MLPH protein to MYO5A in Silverblue animals, resulting in an improper pigmentation as seen in diluted phenotypes. Additionally, in the MLPH mRNA of wt, Violet and Silverblue phenotypes, part of intron 8 is retained resulting in a truncated MLPH protein, which is 359 residues long in wt and Violet and 284 residues long in Silverblue. Subsequently, our findings point out that the missing actin-binding domain, in neither of the 3 analyzed phenotypes affects the transport of melanosomes or the consequent final pigmentation. Moreover, the loss of the major part of the MYO5A domain in the Silverblue MLPH protein seems to be the responsible for the dilute phenotype. Based on our genomic DNA data, genetic tests for selecting Silverblue and Violet carrier animals can be performed in American mink.     

Recessive Mutations in COL25A1 Are a Cause of Congenital Cranial Dysinnervation Disorder

Abnormal ocular motility is a common clinical feature in congenital cranial dysinnervation disorder (CCDD). To date, eight genes related to neuronal development have been associated with different CCDD phenotypes. By using linkage analysis, candidate gene screening, and exome sequencing, we identified three mutations in collagen, type XXV, alpha 1 (COL25A1) in individuals with autosomal-recessive inheritance of CCDD ophthalmic phenotypes. These mutations affected either stability or levels of the protein. We further detected altered levels of sAPP (neuronal protein involved in axon guidance and synaptogenesis) and TUBB3 (encoded by TUBB3, which is mutated in CFEOM3) as a result of null mutations in COL25A1. Our data suggest that lack of COL25A1 might interfere with molecular pathways involved in oculomotor neuron development, leading to CCDD phenotypes.     

Genomic Perspectives on the Emerging SARS-CoV-2 Omicron Variant

A new variant of concern for SARS-CoV-2,Omicron(B.1.1.529),was designated by the World Health Organization on November 26,2021.This study analyzed the viral genome sequencing data of 108 samples collected from patients infected with Omicron.First,we found that the enrichment efficiency of viral nucleic acids was reduced due to mutations in the region where the primers anneal to.Second,the Omicron variant possesses an excessive number of mutations compared to other variants circulating at the sam...     

Cysteine scanning of MscL transmembrane domains reveals residues critical for mechanosensitive channel gating

The mechanosensitive channel of large conductance (MscL), a bacterial channel, is perhaps the best characterized mechanosensitive protein. A structure of the Mycobacterium tuberculosis ortholog has been solved by x-ray crystallography, but details of how the channel gates remain obscure. Here, cysteine scanning was used to identify residues within the transmembrane domains of Escherichia coli MscL that are crucial for normal function. Utilizing genetic screens, we identified several mutations that induced gain-of-function or loss-of-function phenotypes in vivo. Mutants that exhibited the most severe phenotypes were further characterized using electrophysiological techniques and chemical modifications of the substituted cysteines. Our results verify the importance of residues in the putative primary gate in the first transmembrane domain, corroborate other residues previously noted as critical for normal function, and identify new ones. In addition, evaluation of disulfide bridging in native membranes suggests alterations of existing structural models for the “fully closed” state of the channel.     

Dissecting mutational allosteric effects in alkaline phosphatases associated with different Hypophosphatasia phenotypes: An integrative computational investigation

Hypophosphatasia (HPP) is a rare inherited disorder characterized by defective bone mineralization and is highly variable in its clinical phenotype. The disease occurs due to various loss-of-function mutations in ALPL, the gene encoding tissue-nonspecific alkaline phosphatase (TNSALP). In this work, a data-driven and biophysics-based approach is proposed for the large-scale analysis of ALPL mutations-from nonpathogenic to severe HPPs. By using a pipeline of synergistic approaches including sequence-structure analysis, network modeling, elastic network models and atomistic simulations, we characterized allosteric signatures and effects of the ALPL mutations on protein dynamics and function. Statistical analysis of molecular features computed for the ALPL mutations showed a significant difference between the control, mild and severe HPP phenotypes. Molecular dynamics simulations coupled with protein structure network analysis were employed to analyze the effect of single-residue variation on conformational dynamics of TNSALP dimers, and the developed machine learning model suggested that the topological network parameters could serve as a robust indicator of severe mutations. The results indicated that the severity of disease-associated mutations is often linked with mutation-induced modulation of allosteric communications in the protein. This study suggested that ALPL mutations associated with mild and more severe HPPs can exert markedly distinct effects on the protein stability and long-range network communications. By linking the disease phenotypes with dynamic and allosteric molecular signatures, the proposed integrative computational approach enabled to characterize and quantify the allosteric effects of ALPL mutations and role of allostery in the pathogenesis of HPPs.     

Predict prokaryotic proteins through detecting N-formylmethionine residues in protein sequences using support vector machine

Identifying prokaryotes in silico is commonly based on DNA sequences. In experiments where DNA sequences may not be immediately available, we need to have a different approach to detect prokaryotes based on RNA or protein sequences. N-formylmethionine (fMet) is known as a typical characteristic of prokaryotes. A web tool has been implemented here for predicting prokaryotes through detecting the N-formylmethionine residues in protein sequences. The predictor is constructed using support vector machine. An online predictor has been implemented using Python. The implemented predictor is able to achieve the total prediction accuracy 80% with the specificity 80% and the sensitivity 81%.     

Carbohydrate Kinase (RhaK)-Dependent ABC Transport of Rhamnose in Rhizobium leguminosarum Demonstrates Genetic Separation of Kinase and Transport Activities

In Rhizobium leguminosarum the ABC transporter responsible for rhamnose transport is dependent on RhaK, a sugar kinase that is necessary for the catabolism of rhamnose. This has led to a working hypothesis that RhaK has two biochemical functions: phosphorylation of its substrate and affecting the activity of the rhamnose ABC transporter. To address this hypothesis, a linker-scanning random mutagenesis of rhaK was carried out. Thirty-nine linker-scanning mutations were generated and mapped. Alleles were then systematically tested for their ability to physiologically complement kinase and transport activity in a strain carrying an rhaK mutation. The rhaK alleles generated could be divided into three classes: mutations that did not affect either kinase or transport activity, mutations that eliminated both transport and kinase activity, and mutations that affected transport activity but not kinase activity. Two genes of the last class (rhaK72 and rhaK73) were found to have similar biochemical phenotypes but manifested different physiological phenotypes. Whereas rhaK72 conferred a slow-growth phenotype when used to complement rhaK mutants, the rhaK73 allele did not complement the inability to use rhamnose as a sole carbon source. To provide insight to how these insertional variants might be affecting rhamnose transport and catabolism, structural models of RhaK were generated based on the crystal structure of related sugar kinases. Structural modeling suggests that both rhaK72 and rhaK73 affect surface-exposed residues in two distinct regions that are found on one face of the protein, suggesting that this protein''s face may play a role in protein-protein interaction that affects rhamnose transport.     

Effector-Triggered Immune Response in Arabidopsis thaliana Is a Quantitative Trait

We identified loci responsible for natural variation in Arabidopsis thaliana (Arabidopsis) responses to a bacterial pathogen virulence factor, HopAM1. HopAM1 is a type III effector protein secreted by the virulent Pseudomonas syringae strain Pto DC3000. Delivery of HopAM1 from disarmed Pseudomonas strains leads to local cell death, meristem chlorosis, or both, with varying intensities in different Arabidopsis accessions. These phenotypes are not associated with differences in bacterial growth restriction. We treated the two phenotypes as quantitative traits to identify host loci controlling responses to HopAM1. Genome-wide association (GWA) of 64 Arabidopsis accessions identified independent variants highly correlated with response to each phenotype. Quantitative trait locus (QTL) mapping in a recombinant inbred population between Bur-0 and Col-0 accessions revealed genetic linkage to regions distinct from the top GWA hits. Two major QTL associated with HopAM1-induced cell death were also associated with HopAM1-induced chlorosis. HopAM1-induced changes in Arabidopsis gene expression showed that rapid HopAM1-dependent cell death in Bur-0 is correlated with effector-triggered immune responses. Studies of the effect of mutations in known plant immune system genes showed, surprisingly, that both cell death and chlorosis phenotypes are enhanced by loss of EDS1, a regulatory hub in the plant immune-signaling network. Our results reveal complex genetic architecture for response to this particular type III virulence effector, in contrast to the typical monogenic control of cell death and disease resistance triggered by most type III effectors.     

Multiple gene mutations, not the type of mutation, are the modifier of left ventricle hypertrophy in patients with hypertrophic cardiomyopathy

Genotype-phenotype correlation of hypertrophic cardiomyopathy (HCM) has been challenging because of the genetic and clinical heterogeneity. To determine the mutation profile of Chinese patients with HCM and to correlate genotypes with phenotypes, we performed a systematic mutation screening of the eight most commonly mutated genes encoding sarcomere proteins in 200 unrelated Chinese adult patients using direct DNA sequencing. A total of 98 mutations were identified in 102 mutation carriers. The frequency of mutations in MYH7, MYBPC3, TNNT2 and TNNI3 was 26.0, 18.0, 4.0 and 3.5 % respectively. Among the 200 genotyped HCM patients, 83 harbored a single mutation, and 19 (9.5 %) harbored multiple mutations. The number of mutations was positively correlated with the maximum wall thickness. We found that neither particular gene nor specific mutation was correlated to clinical phenotype. In summary, the frequency of multiple mutations was greater in Chinese HCM patients than in the Caucasian population. Multiple mutations in sarcomere protein may be a risk factor for left ventricular wall thickness.     

Catalytic mechanism of ancestral L-lysine oxidase assigned by sequence data mining

A large number of protein sequences are registered in public databases such as PubMed. Functionally uncharacterized enzymes are included in these databases, some of which likely have potential for industrial applications. However, assignment of the enzymes remained difficult tasks for now. In this study, we assigned a total of 28 original sequences to uncharacterized enzymes in the FAD-dependent oxidase family expressed in some species of bacteria including Chryseobacterium, Flavobacterium, and Pedobactor. Progenitor sequence of the assigned 28 sequences was generated by ancestral sequence reconstruction, and the generated sequence exhibited L-lysine oxidase activity; thus, we named the enzyme AncLLysO. Crystal structures of ligand-free and ligand-bound forms of AncLLysO were determined, indicating that the enzyme recognizes L-Lys by hydrogen bond formation with R76 and E383. The binding of L-Lys to AncLLysO induced dynamic structural change at a plug loop formed by residues 251 to 254. Biochemical assays of AncLLysO variants revealed the functional importance of these substrate recognition residues and the plug loop. R76A and E383D variants were also observed to lose their activity, and the k_cat/K_m value of G251P and Y253A mutations were approximately 800- to 1800-fold lower than that of AncLLysO, despite the indirect interaction of the substrates with the mutated residues. Taken together, our data demonstrate that combinational approaches to sequence classification from database and ancestral sequence reconstruction may be effective not only to find new enzymes using databases of unknown sequences but also to elucidate their functions.     

Enhanced Mutability Associated with a Temperature-Sensitive Mutant of Vesicular Stomatitis Virus

Temperature-sensitive (ts) mutant tsD1 of vesicular stomatitis virus, New Jersey serotype, is the sole representative of complementation group D. Clones derived from this mutant exhibited three different phenotypes with respect to electrophoretic mobility of the G and N polypeptides of the virion in sodium dodecyl sulfate-polyacrylamide gel. Analysis of non-ts pseudorevertants showed that none of the three phenotypes was associated with the temperature sensitivity of mutant tsD1. Additional phenotypes, some also involving the NS polypeptide, appeared during sequential cloning, indicating that mutations were generated at high frequency during replication of tsD1. Furthermore, mutations altering the electrophoretic mobility of the G, N, NS, and M polypeptides were induced in heterologous viruses multiplying in the same cells as tsD1. These heterologous viruses included another complementing ts mutant of vesicular stomatitis virus New Jersey and ts mutants of vesicular stomatitis virus Indiana and Chandipura virus. Complete or incomplete virions of tsD1 appeared to be equally efficient inducers of mutations in heterologous viruses. Analysis of the progeny of a mixed infection of two complementing ts mutants of vesicular stomatitis virus New Jersey with electrophoretically distinguishable G, N, NS, and M proteins yielded no recombinants and excluded recombination as a factor in the generation of the electrophoretic mobility variants. In vitro translation of total cytoplasmic RNA from BHK cells indicated that post-translational processing was not responsible for the aberrant electrophoretic mobility of the N, NS, and M protein mutants. Aberrant glycosylation could account for three of four G protein mutants, however. Some clones of tsD1 had an N polypeptide which migrated faster in sodium dodecyl sulfate-polyacrylamide gel than did the wild type, suggesting that the polypeptide might be shorter by about 10 amino acids. Determination of the nucleotide sequence to about 200 residues from each terminus of the N gene of one of these clones, a revertant, and the wild-type parent revealed no changes compatible with synthesis of a shorter polypeptide by premature termination or late initiation of translation. The sequence data indicated, however, that the N-protein mutant and its revertant differed from the parental wild type in two of the 399 nucleotides determined. These sequencing results and the phenomenon of enhanced mutability associated with mutant tsD1 reveal that rapid and extensive evolution of the viral genome can occur during the course of normal cytolytic infection of cultured cells.     

PPIcons: identification of protein-protein interaction sites in selected organisms

The physico-chemical properties of interaction interfaces have a crucial role in characterization of protein–protein interactions (PPI). In silico prediction of participating amino acids helps to identify interface residues for further experimental verification using mutational analysis, or inhibition studies by screening library of ligands against given protein. Given the unbound structure of a protein and the fact that it forms a complex with another known protein, the objective of this work is to identify the residues that are involved in the interaction. We attempt to predict interaction sites in protein complexes using local composition of amino acids together with their physico-chemical characteristics. The local sequence segments (LSS) are dissected from the protein sequences using a sliding window of 21 amino acids. The list of LSSs is passed to the support vector machine (SVM) predictor, which identifies interacting residue pairs considering their inter-atom distances. We have analyzed three different model organisms of Escherichia coli, Saccharomyces Cerevisiae and Homo sapiens, where the numbers of considered hetero-complexes are equal to 40, 123 and 33 respectively. Moreover, the unified multi-organism PPI meta-predictor is also developed under the current work by combining the training databases of above organisms. The PPIcons interface residues prediction method is measured by the area under ROC curve (AUC) equal to 0.82, 0.75, 0.72 and 0.76 for the aforementioned organisms and the meta-predictor respectively.     

Caenorhabditis elegans SUR-5, a Novel but Conserved Protein, Negatively Regulates LET-60 Ras Activity during Vulval Induction

The let-60 ras gene acts in a signal transduction pathway to control vulval differentiation in Caenorhabditis elegans. By screening suppressors of a dominant negative let-60 ras allele, we isolated three loss-of-function mutations in the sur-5 gene which appear to act as negative regulators of let-60 ras during vulval induction. sur-5 mutations do not cause an obvious mutant phenotype of their own, and they appear to specifically suppress only one of the two groups of let-60 ras dominant negative mutations, suggesting that the gene may be involved in a specific aspect of Ras activation. Consistent with its negative function, overexpressing sur-5 from an extragenic array partially suppresses the Multivulva phenotype of an activated let-60 ras mutation and causes synergistic phenotypes with a lin-45 raf mutation. We have cloned sur-5 and shown that it encodes a novel protein. We have also identified a potential mammalian SUR-5 homolog that is about 35% identical to the worm protein. SUR-5 also has some sequence similarity to acetyl coenzyme A synthetases and is predicted to contain ATP/GTP and AMP binding sites. Our results suggest that sur-5 gene function may be conserved through evolution.     

Insightful directed evolution of Escherichia coli quorum sensing promoter region of the lsrACDBFG operon: a tool for synthetic biology systems and protein expression

Quorum sensing (QS) regulates many natural phenotypes (e.q. virulence, biofilm formation, antibiotic resistance), and its components, when incorporated into synthetic genetic circuits, enable user-directed phenotypes. We created a library of Escherichia coli lsr operon promoters using error-prone PCR (ePCR) and selected for promoters that provided E. coli with higher tetracycline resistance over the native promoter when placed upstream of the tet(C) gene. Among the fourteen clones identified, we found several mutations in the binding sites of QS repressor, LsrR. Using site-directed mutagenesis we restored all p-lsrR-box sites to the native sequence in order to maintain LsrR repression of the promoter, preserving the other mutations for analysis. Two promoter variants, EP01rec and EP14rec, were discovered exhibiting enhanced protein expression. In turn, these variants retained their ability to exhibit the LsrR-mediated QS switching activity. Their sequences suggest regulatory linkage between CytR (CRP repressor) and LsrR. These promoters improve upon the native system and exhibit advantages over synthetic QS promoters previously reported. Incorporation of these promoters will facilitate future applications of QS-regulation in synthetic biology and metabolic engineering.