Reducing mutational bias in random protein libraries

The success of protein optimization through directed molecular evolution depends to a large extent on the size and quality of the displayed library. Current low-fidelity DNA polymerases that are commonly used during random mutagenesis and recombination in vitro display strong mutational preferences, favoring the substitution of certain nucleotides over others. The result is a biased and reduced functional diversity in the library under selection. In an effort to reduce mutational bias, we combined two different low-fidelity DNA polymerases, Taq and Mutazyme, which have opposite mutational spectra. As a first step, random mutants of the Bacillus thuringiensis cry9Ca1 gene were generated by separate error-prone polymerase chain reactions (PCRs) with each of the two polymerases. Subsequent shuffling by staggered extension process (StEP) of the PCR products resulted in intermediate numbers of AT and GC substitutions, compared to the Taq or Mutazyme error-prone PCR libraries. This strategy should allow generating unbiased libraries or libraries with a specific degree of mutational bias by applying optimal mutagenesis frequencies during error-prone PCR and controlling the concentration of template in the shuffling reaction while taking into account the GC content of the target gene.     

用于蛋白质体外分子进化研究的DNA随机突变技术

蛋白质体外分子进化是模拟自然的进化过程,利用基因随机突变和定向筛选（选择）技术,以获得具有预期新功能的突变体分子。虽然体外进化近几年才产生,但已成为医药和工业领域中筛选具有特殊催化性质的酶的最重要的方法之一。ＤＮＡ随机突变技术是蛋白质体外分子进化研究的基础,本文将对几种最重要的突变方法：倾向错误的ＰＣＲ、ＤＮＡ重排、模板交错延伸反应和随机延伸突变的原理和应用等加以介绍。     

通过原位易错PCR一步构建基因突变文库

主要介绍一种通过原位易错PCR构建随机突变文库的新技术。本实验室最近发表的一项国际专利中,利用来源于海栖热孢菌的极耐热性DNA连接酶,在传统PCR循环中加入一个连接步骤,即变性—退火—延伸—连接的四步循环法PCR,从而实现环状质粒的PCR指数扩增(PPCP)。原位易错PCR中所用引物为一段线性双链DNA,它含有与模板质粒不同的筛选标记,产物转化宿主菌后,模板质粒在筛选平板上被直接剔除。筛选到的阳性突变子可用作模板直接进入突变文库的再次构建,通过筛选获得二级或多级累加的正突变。利用这种方法构建了一个木聚糖酶基因和一个纤维素酶基因的随机突变文库,并筛选出具有正向突变的蛋白,证明以PPCP为基础的原位易错PCR技术,为基因定向进化提供了一种快速有效的随机突变文库构建的新方法。     

Improved PCR method for the creation of saturation mutagenesis libraries in directed evolution: application to difficult-to-amplify templates

Saturation mutagenesis constitutes a powerful method in the directed evolution of enzymes. Traditional protocols of whole plasmid amplification such as Stratagene’s QuikChange™ sometimes fail when the templates are difficult to amplify. In order to overcome such restrictions, we have devised a simple two-primer, two-stage polymerase chain reaction (PCR) method which constitutes an improvement over existing protocols. In the first stage of the PCR, both the mutagenic primer and the antiprimer that are not complementary anneal to the template. In the second stage, the amplified sequence is used as a megaprimer. Sites composed of one or more residues can be randomized in a single PCR reaction, irrespective of their location in the gene sequence.The method has been applied to several enzymes successfully, including P450-BM3 from Bacillus megaterium, the lipases from Pseudomonas aeruginosa and Candida antarctica and the epoxide hydrolase from Aspergillus niger. Here, we show that megaprimer size as well as the direction and design of the antiprimer are determining factors in the amplification of the plasmid. Comparison of the results with the performances of previous protocols reveals the efficiency of the improved method. Joaquin Sanchis, Layla Fernández, and J. Daniel Carballeira contributed equally.     

一个新的L-半胱亚磺酸脱羧酶基因的分析和突变

【目的】完成一个来源于碱性污染土壤宏基因组文库的新的L-半胱亚磺酸脱羧酶基因undec1A的鉴定,研究其酶学性质并利用非理性设计技术对其进行分子改造。【方法】以p ETBlue-2为表达载体构建包含undec1A基因的重组表达质粒,转化至宿主细胞E.coli Tuner(DE3)p Lac?中构建重组表达克隆,采用镍亲和层析完成了酶蛋白的分离纯化,完成其生化特征研究,利用连续易错PCR技术对Undec1A蛋白进行分子改造。【结果】生物信息学分析结果揭示Undec1A蛋白与已知的L-半胱亚磺酸脱羧酶存在类似的辅酶结合位点和底物识别催化基序等。分子对接结果显示氨基酸残基Val237、Asp239、Asp266、Ile267、Ala268和Lys298等决定了与底物分子L-半胱亚磺酸的识别和结合催化。以L-半胱亚磺酸作为底物,重组Undec1A蛋白的最适作用p H为7.0,最适作用温度为35°C;分子动力学常数K_m为(1.557±0.015)mmol/L,V_(max)为(49.07±3.19)μmol/(L·min),k_(cat)为(45.80±1.32)/min。利用连续易错PCR技术完成了亲本酶的分子改造,分离筛选到了一个酶活力更高的突变酶Undec1A-1180。在优化条件下,Undec1A-1180的比活力较亲本酶提高了约5.62倍。【结论】本研究为构建牛磺酸的生物合成工艺提供了理论参考,因而具有重要的实践意义。     

Engineering of a bacterial tyrosinase for improved catalytic efficiency towards D‐tyrosine using random and site directed mutagenesis approaches

The tyrosinase gene from Ralstonia solanacearum (GenBank NP518458) was subjected to random mutagenesis resulting in tyrosinase variants (RVC10 and RV145) with up to 3.2‐fold improvement in k_cat, 5.2‐fold lower K_m and 16‐fold improvement in catalytic efficiency for D ‐tyrosine. Based on RVC10 and RV145 mutated sequences, single mutation variants were generated with all variants showing increased k_cat for D ‐tyrosine compared to the wild type (WT). All single mutation variants based on RV145 had a higher k_cat and K_m value compared to the RV145 and thus the combination of four mutations in RV145 was antagonistic for turnover, but synergistic for affinity of the enzyme for D ‐tyrosine. Single mutation variant 145_V153A exhibited the highest (6.9‐fold) improvement in k_cat and a 2.4‐fold increase in K_m compared to the WT. Two single mutation variants, C10_N322S and C10_T183I reduced the K_m up to 2.6‐fold for D ‐tyrosine but one variant 145_V153A increased the K_m 2.4‐fold compared to the WT. Homology based modeling of R. solanacearum tyrosinase showed that mutation V153A disrupts the van der Waals interactions with an α‐helix providing one of the conserved histidine residues of the active site. The k_cat and K_m values for L ‐tyrosine decreased for RV145 and RVC10 compared to the WT. RV145 exhibited a 2.1‐fold high catalytic efficiency compared to the WT which is a 7.6‐fold lower improvement compared to D ‐tyrosine. RV145 exhibited a threefold higher monophenolase:diphenolase activity ratio for D ‐tyrosine:D ‐DOPA and a 1.4‐fold higher L ‐tyrosine:L ‐DOPA activity ratio compared to the WT. Biotechnol. Bioeng. 2013; 110: 1849–1857. © 2013 Wiley Periodicals, Inc.     

Asymmetric overlap extension PCR method bypassing intermediate purification and the amplification of wild-type template in site-directed mutagenesis

By combining asymmetric PCR and overlap extension, we developed a novel asymmetric overlap extension PCR (AOE-PCR) method for site-directed mutagenesis which bypassed the need for intermediate purification and excluded the amplification of a wild-type template. This method was used to introduce single base mutations into a small GTPase gene from cotton and to simultaneously introduce two mutations just by repeating this method using the first round AOE-PCR products as template. Our results suggested that the AOE-PCR method represents a valuable improvement of the original overlap extension PCR for site-directed mutagenesis.     

A Deterministic Approximation to the Number of Mutants in Growing Clones

The use of an approximation to the median of the Poisson distribution to represent each occurrence of mutations in a growing clone permits the prediction of the number of mutants per clone without the limitations imposed by more heuristic expressions. Its application to the evaluation of mutation rates yields results comparable to those obtained by fluctuation analysis.     

Mutational fitness effects in RNA and single-stranded DNA viruses: common patterns revealed by site-directed mutagenesis studies

The fitness effects of mutations are central to evolution, yet have begun to be characterized in detail only recently. Site-directed mutagenesis is a powerful tool for achieving this goal, which is particularly suited for viruses because of their small genomes. Here, I discuss the evolutionary relevance of mutational fitness effects and critically review previous site-directed mutagenesis studies. The effects of single-nucleotide substitutions are standardized and compared for five RNA or single-stranded DNA viruses infecting bacteria, plants or animals. All viruses examined show very low tolerance to mutation when compared with cellular organisms. Moreover, for non-lethal mutations, the mean fitness reduction caused by single mutations is remarkably constant (0.10–0.13), whereas the fraction of lethals varies only modestly (0.20–0.41). Other summary statistics are provided. These generalizations about the distribution of mutational fitness effects can help us to better understand the evolution of RNA and single-stranded DNA viruses.     

Improvement in thermal stability and reactivity of pyranose oxidase from Coriolus versicolor by random mutagenesis

A mutant producing a pyranose oxidase, which has a higher thermal stability and lower Km values for d-glucose and 1,5-anhydro-d-glucitol than those of the wild type enzyme, was obtained. A single amino acid substitution, Lys for Glu at position 542, had occurred. This altered enzyme, E542K, was not only stable at 55°C, which was 5°C higher than the wild-type enzyme, but was stable in alkaline solution at pH 8.0–11.0. Km values of E542K for d-glucose and 1,5-anhydro-d-glucitol were 0.7 mM and 14.3 mM, respectively, in contrast with 1.4 mm and 35.3 mM for the wild-type enzyme. A little effect was observed in kcat values, and improvement in reactivity was mainly due to the decreases in Km values. This altered pyranose oxidase is useful for food analysis and diagnosis.     

High-throughput insertion mutagenesis and functional screening in the entomopathogenic fungus Beauveria bassiana

The entomopathogenic fungus Beauveria bassiana displays a broad insect host range and serves as a model for examining host-pathogen interactions. Rapid construction and screening of random-insertion mutants of B. bassiana provides a powerful tool to dissect the molecular mechanisms of fungal virulence. LiAc/DMSO treated B. bassiana blastospores were found to be highly competent to transformation using linear DNA and a polyethylene glycol-based method. Selection on cellophane-layered Czapek-Dox agar at a lowered pH (from 7.5 to 6.3) greatly decreased background growth of non-transformed cells and improved screening of transformants. Optimization of the protocol using integration of the bar phosphinothricin resistance gene resulted in high transformation rates (200-250 transformants/μg DNA/10⁸ cells). A collection of ∼4000 insertion mutants was examined via high-throughput screens for hydrocarbon utilization. One mutant was isolated that grew poorly on both n-hexadecane and tributyrin. The random insertion site was mapped to a gene that displayed homology to vitamin H (biotin)/tartrate transporters. Insect bioassays using Galleria mellonella as the target host revealed decreased virulence in the mutant. This system provides a simple and rapid method for the generation and screening of insertion mutants and should expand our ability to genetically analyze the B. bassiana lifestyle.     

拟南芥钙调素定点突变基因分离及其在钙不依赖钙调素结合蛋白检测中的应用

动植物系统研究表明,钙调素不仅在结合钙离子时调节多种靶酶或靶蛋白的活性,而且没有钙离子结合时,还可以通过结合钙不依赖的钙调素结合蛋白,发挥多种生物学作用．然而,目前却没有体内分析钙调素与钙不依赖钙调素结合蛋白相互作用的方法．首先,采用定点突变的方式,得到了拟南芥钙调素亚型2的多个突变基因mCaM2,随后,大肠杆菌重组表达突变蛋白的电泳迁移率及⁴⁵Ca²⁺覆盖分析表明,得到了编码失去钙结合能力的钙调素的突变基因mCaM2₁₂₃₄, mCaM2₁₂₃₄突变钙调素中所有4个钙结合EF-hand结构域中的关键氨基酸谷氨酸均突变为谷氨酰胺．在酵母双杂交体系中,作为诱饵蛋白的突变钙调素mCaM2₁₂₃₄与我们前期体外方法报道的钙不依赖性钙调素结合蛋白AtIQD26存在相互作用．这将为钙不依赖性钙调素结合蛋白提供有用的体内研究工具,有利于我们全面认识钙-钙调素-钙调素结合蛋白信号途径．     

Patterns of nucleotide substitution in pseudogenes and functional genes

Summary The pattern of point mutations is inferred from nucleotide substitutions in pseudogenes. The pattern obtained suggests that transition mutations occur somewhat more frequently than transversion mutations and that mutations result more often in A or T than in G or C. Our results are discussed with respect to the predictions from Topal and Fresco's model for the molecular basis of point (substitution) mutations (Nature 263:285–289, 1976). The pattern of nucleotide substitution at the first and second positions of codons in functional genes is quite similar to that in pseudogenes, but the relative frequency of the transition CT in the sense strand is drastically reduced and those of the transversions CG and GC are doubled. The differences between the two patterns can be explained by the observation that in the protein evolution amino acid substitutions occur mainly between amino acids with similar biochemical properties (Grantham, Science 185:862–864, 1974). Our results for the patterns of nucleotide substitutions in pseudogenes and in functional genes lead to the prediction that both the coding and non-coding regions of protein coding genes should have high frequencies of A and T. Available data show that the non-coding regions are indeed high in A and T but the coding regions are low in T, though high in A.     

Molecular cloning of a cDNA encoding mouse D-aspartate oxidase and functional characterization of its recombinant proteins by site-directed mutagenesis

Summary. The cDNA encoding D-aspartate oxidase (DASPO) was cloned from mouse kidney RNA by RT–PCR. Sequence analysis showed that it contained a 1023-bp open reading frame encoding a protein of 341 amino acid residues. The protein was expressed in Escherichia coli with or without an N-terminal His-tag and had functional DASPO activity that was highly specific for D-aspartate and N-methyl-D-aspartate. To investigate the roles of the Arg-216 and Arg-237 residues of the mouse DASPO (mDASPO), we generated clones with several single amino acid substitutions of these residues in an N-terminally His-tagged mDASPO. These substitutions significantly reduced the activity of the recombinant enzyme against acidic D-amino acids and did not confer any additional specificity to other amino acids. These results suggest that the Arg-216 and Arg-237 residues of mDASPO are catalytically important for full enzyme activity.     

The pathogenicity,structural and functional exploration of human HMGB1 single nucleotide polymorphisms using in silico study

Abstract

The human HMGB1 gene mutations have a major impact on several immune-related diseases and cancer. The detrimental effect of non-synonymous mutations of HMGB1 has not been investigated yet, hence the present study aims to examine single nucleotide polymorphisms and their implications on the structure-function of human HMGB1. The multifaceted HMGB1 protein acts as pleiotropic cytokine and regulates essential genes for coordinated cellular functions. The mutational effect on HMGB1 was analyzed by sequence-based homology methods, supervised learning methods, and structure-based methods. The study identified 58 non-synonymous mutations in human HMGB1, out of which only 2 mutations; R10T (rs61742222) and F103C (rs61733675) were classified as the SNPs with highest deleterious and disease-causing mutants. The effect of these mutations in structure of HMGB1 was scrutinized and the R10T mutant found to have a distinct structural behaviour in the B-box domain. In addition, R10T mutant predicted that it affects the MoRF function of HMGB1 and it could disrupt the DNA binding or/and protein partner interaction activity by HMGB1. F103C mutation takes place at the TLR binding and cytokine inducing region of HMGB1, hence it could affect the protein binding activity which involves in many cellular signaling. The study identified potent mutations R10T (a cancer-causing somatic mutation) and F103C (a novel mutation) and these mutations either directly or indirectly hinder DNA binding activity and TLR and cytokine binding of HMGB1. These findings will help in understanding the molecular basis of these promising mutations and functional role of human HMGB1 in cancer and immunological diseases.

Abbreviations AGER Advanced glycosylation end product-specific receptor

CXCL Chemokine (C-X-C motif) ligand

dbSNP The single nucleotide polymorphism database

HMGB1 High mobility group box 1

LINCS LINear Constraint Solver

MDS Molecular dynamics simulation

MoRF Molecular recognition features

NPT Number of particle, Pressure and Temperature

NVT Number of particle, Volume and Temperature

nsSNP Non-synonymous SNP

PBC Partial boundary condition

PCA Principal component analysis

PME Partial mesh Ewald

RMSD Root mean square deviation

RMSF Root mean square fluctuation

SNP Single nucleotide polymorphism

SPC Single-point charge

TLR Toll-like receptor

UTR Un-translated Region

Communicated by Ramaswamy H. Sarma     

Different efficiency of UmuDC and MucAB proteins in UV light induced mutagenesis in Escherichia coli

Summary Two multicopy plasmids carrying either the umuDC or the mucAB operon were used to compare the efficiency of UmuDC and MucAB proteins in UV mutagenesis of Escherichia coli K12. It was found that in recA ⁺ uvr ⁺bacteria, plasmid pIC80, mucAB ⁺mediated UV mutagenesis more efficiently than did plasmid pSE117, umuDC ⁺. A similar result was obtained in lexA51(Def) cells, excluding the possibility that this was due to a differential regulation by LexA of the umuDC and mucAB operons. We conclude that some structural characteristic of the UmuDC and MucAB proteins determines their different efficiency in UV mutagenesis. This characteristic could be also responsible for the observation that in the recA430 mutant, pIC80 but no pSE117 can mediate UV mutagenesis. In the recA142 mutant, pIC80 also promoted UV mutagenesis more efficiently than pSE117. In this mutant, the recombination proficiency, the protease activity toward LexA and the mutation frequency were increased by the presence of adenine in the medium. In recA ⁺ uvrB5 bacteria, plasmid pSE117,umuDC caused both an increase in UV sensitivity as well as a reduction in the mutation frequency. These nagative effects resulting from the overproduction of UmuDC proteins were higher in recA142 uvrB5 than in recA ⁺ uvrB5 cells. In contrast, overproduction of MucAB proteins in excision-deficient bacteria containing pIC80 led to a large increase in the mutation frequency. We suggest that the functional differences between UmuDC and MucAB proteins might be due to their different dependence on the direct role of RecA protease in UV mutagenesis.     

The ecology and genetics of fitness in Chlamydomonas. IX. The rate of accumulation of variation of fitness under selection

Populations of Chlamydomonas founded by single cells were cultured in chemostats for 50 days, representing about 125 generations. The mean and variance of division rate was measured daily by withdrawing cells from the effluent and culturing them for 24 h on filtered effluent medium solidified with agar. Mean fitness did not change during the period of culture, and the behavior of neutral markers indicated that no substitutions of novel beneficial mutations occurred. However, the variance of fitness increased markedly at about the same rate in two replicate populations. The standardized rate, or mutational heritability, was Vm/VE = 4-5 x 10(-3) per generation. This is substantially greater than most other estimates for characters closely correlated with fitness. Moreover, it seems difficult to reconcile with the absence of any change in mean fitness. We investigated the possibility that frequency-dependent selection was created by spatial heterogeneity within the culture vessel by testing cell populations with different phenotypes from the top, bottom, and surface of the chemostats. However, the differentiation of these populations seemed to be attributable to phenotypic plasticity, with no evidence that their characteristics were heritable. Finally, we report an experiment in which lines were selected for about 100 generations on solid or liquid medium. These lines became specifically adapted to the medium on which they were cultured, showing that liquid and solid media, even when chemically identical, provide different conditions of growth for Chlamydomonas. The genetic variance appearing in the cultures was therefore attributed to conditionally neutral mutations that were not expressed in the chemostat. This implies that rates of accumulation of mutational variance measured in the culture environment itself (where this can be done) may greatly underestimate the variation available for a response through selection to environmental change. Moreover, it suggests that chemostat populations may be more dynamic and more diverse than is usually thought.