Robustness of hen lysozyme monitored by random mutations

Chemoselective Caging of Carboxyl Groups for On-Demand Protein Activation with Small Molecules

Tools for on-demand protein activation enable impactful gain-of-function studies in biological settings. Thus far, however, proteins have been chemically caged at primarily Lys, Tyr, and Sec, typically through the genetic encoding of unnatural amino acids. Herein, we report that the preferential reactivity of diazo compounds with protonated acids can be used to expand this toolbox to solvent-accessible carboxyl groups with an elevated pKa value. As a model protein, we employed lysozyme (Lyz), which has an active-site Glu35 residue with a pKa value of 6.2. A diazo compound with a bioorthogonal self-immolative handle esterified Glu35 selectively, inactivating Lyz. The hydrolytic activity of the caged Lyz on bacterial cell walls was restored with two small-molecule triggers. The decaging was more efficient by small molecules than by esterases. This simple chemical strategy was also applied to a hemeprotein and an aspartyl protease, setting the stage for broad applicability.

Fragment-Directed Random Mutagenesis by the Reverse Kunkel Method

Two fundamentally different approaches are routinely used for protein engineering: user-defined mutagenesis and random mutagenesis, each with its own strengths and weaknesses. Here, we invent a unique mutagenesis protocol, which combines the advantages of user-defined mutagenesis and random mutagenesis. The new method, termed the reverse Kunkel method, allows the user to create random mutations at multiple specified regions in a one-pot reaction. We demonstrated the reverse Kunkel method by mimicking the somatic hypermutation in antibodies that introduces random mutations concentrated in complementarity-determining regions. Coupling with the phage display and yeast display selections, we successfully generated dramatically improved antibodies against a model protein and a neurotransmitter peptide in terms of affinity and immunostaining performance. The reverse Kunkel method is especially suitable for engineering proteins whose activities are determined by multiple variable regions, such as antibodies and adeno-associated virus capsids, or whose functional domains are composed of several discontinuous sequences, such as Cas9 and Cas12a.

Global regulator engineering enhances bioelectricity generation in Pseudomonas aeruginosa-inoculated MFCs

The electricigens with high-electroactivity is essential for resolving the low electricity power output (EPT) of microbial fuel cells (MFCs). However, the manipulation by single functional genes shows limitation because electroactivity is a complex phenotype controlled by multiple genes. Herein, global regulator engineering (GRE) was developed to optimize the electroactivity of an isolated strain (Pseudomonas aeruginosa P3-A-11) using an exogenous global regulator IrrE (ionizing radiation resistance E linkage group) as an object. The GRE was implemented through in vitro random mutagenesis by error-prone PCR and in vivo high-through screening comprised of cultures color assay, PYO measurement and MFCs operation. Four mutants with higher electroactivity were obtained, among which, the mutant 11/M2-59 not only displayed the maximal power density, but also exhibited stronger salt tolerance, consequently showing good performance of MFCs in the presence of salt. Apart from the reduced internal resistance, the increase in phenazines amounts primarily contributed to EPT improvement, which was realized by enhancing the core biosynthesis pathway and affecting other pathways (such as central metabolism pathway, quorum sensing system, regulatory network). Notably, IrrE exerted its positive effect on electroactivity even without native regulators (such as PmpR and RpoS). In addition, the significant fluctuations in expression levels of stress-responsive genes mediated by GRE were closely associated with the enhanced salt tolerance. This work demonstrated that GRE was an effective approach for simultaneously optimizing multiple phenotypes (such as electroactivity and stress tolerance), and thus would provide more opportunities to create high-efficiency electricigens and further promoted the practical application of MFCs.

Requirements for efficient correction of ΔF508 CFTR revealed by analyses of evolved sequences

Misfolding of ΔF508 cystic fibrosis (CF) transmembrane conductance regulator (CFTR) underlies pathology in most CF patients. F508 resides in the first nucleotide-binding domain (NBD1) of CFTR near a predicted interface with the fourth intracellular loop (ICL4). Efforts to identify small molecules that restore function by correcting the folding defect have revealed an apparent efficacy ceiling. To understand the mechanistic basis of this obstacle, positions statistically coupled to 508, in evolved sequences, were identified and assessed for their impact on both NBD1 and CFTR folding. The results indicate that both NBD1 folding and interaction with ICL4 are altered by the ΔF508 mutation and that correction of either individual process is only partially effective. By contrast, combination of mutations that counteract both defects restores ΔF508 maturation and function to wild-type levels. These results provide a mechanistic rationale for the limited efficacy of extant corrector compounds and suggest approaches for identifying compounds that correct both defective steps.

The stability effects of protein mutations appear to be universally distributed

How the thermodynamic stability effects of protein mutations (DeltaDeltaG) are distributed is a fundamental property related to the architecture, tolerance to mutations (mutational robustness), and evolutionary history of proteins. The stability effects of mutations also dictate the rate and dynamics of protein evolution, with deleterious mutations being the main inhibitory factor. Using the FoldX algorithm that attempts to computationally predict DeltaDeltaG effects of mutations, we deduced the overall distributions of stability effects for all possible mutations in 21 different globular, single domain proteins. We found that these distributions are strikingly similar despite a range of sizes and folds, and largely follow a bi-Gaussian function: The surface residues exhibit a narrow distribution with a mildly destabilizing mean DeltaDeltaG ( approximately 0.6 kcal/mol), whereas the core residues exhibit a wider distribution with a stronger destabilizing mean ( approximately 1.4 kcal/mol). Since smaller proteins have a higher fraction of surface residues, the relative weight of these single distributions correlates with size. We also found that proteins evolved in the laboratory follow an essentially identical distribution, whereas de novo designed folds show markedly less destabilizing distributions (i.e. they seem more robust to the effects of mutations). This bi-Gaussian model provides an analytical description of the predicted distributions of mutational stability effects. It comprises a novel tool for analyzing proteins and protein models, for simulating the effect of mutations under evolutionary processes, and a quantitative description of mutational robustness.

Why High-error-rate Random Mutagenesis Libraries are Enriched in Functional and Improved Proteins

The fraction of proteins that retain wild-type function after mutation has long been observed to decline exponentially as the average number of mutations per gene increases. Recently, several groups have used error-prone polymerase chain reactions (PCR) to generate libraries with 15 to 30 mutations per gene, on average, and have reported that orders of magnitude more proteins retain function than would be expected from the low-mutation-rate trend. Proteins with improved or novel function were isolated disproportionately from these high-error-rate libraries, leading to claims that high mutation rates unlock regions of sequence space that are enriched in positively coupled mutations. Here, we show experimentally that error-prone PCR produces a broader non-Poisson distribution of mutations consistent with a detailed model of PCR. As error rates increase, this distribution leads directly to the observed excesses in functional clones. We then show that while very low mutation rates result in many functional sequences, only a small number are unique. By contrast, very high mutation rates produce mostly unique sequences, but few retain function. Thus an optimal mutation rate exists that balances uniqueness and retention of function. Overall, high-error-rate mutagenesis libraries are enriched in improved sequences because they contain more unique, functional clones. Our findings demonstrate how optimal error-prone PCR mutation rates may be calculated, and indicate that "optimal" rates depend on both the protein and the mutagenesis protocol.

On the conservative nature of intragenic recombination

Intragenic recombination rapidly creates protein sequence diversity compared with random mutation, but little is known about the relative effects of recombination and mutation on protein function. Here, we compare recombination of the distantly related beta-lactamases PSE-4 and TEM-1 to mutation of PSE-4. We show that, among beta-lactamase variants containing the same number of amino acid substitutions, variants created by recombination retain function with a significantly higher probability than those generated by random mutagenesis. We present a simple model that accurately captures the differing effects of mutation and recombination in real and simulated proteins with only four parameters: (i) the amino acid sequence distance between parents, (ii) the number of substitutions, (iii) the average probability that random substitutions will preserve function, and (iv) the average probability that substitutions generated by recombination will preserve function. Our results expose a fundamental functional enrichment in regions of protein sequence space accessible by recombination and provide a framework for evaluating whether the relative rates of mutation and recombination observed in nature reflect the underlying imbalance in their effects on protein function.

PCR-Induced Sequence Alterations Hamper the Typing of Prehistoric Bone Samples for Diagnostic Achondroplasia Mutations

Achondroplasia (ACH) is a skeletal disorder (MIM100800) with an autosomal dominant Mendelian inheritance and complete penetrance. Here we report the screening of ancient bone samples for diagnostic ACH mutations. The diagnostic G-->A transition in the FGFR3 gene at cDNA position 1138 was detected in cloned polymerase chain reaction (PCR) products obtained from the dry mummy of the Semerchet tomb, Egypt (first dynasty, approximately 4,890-5,050 BP [before present]), and from an individual from Kirchheim, Germany (Merovingian period, approximately 1,300-1,500 BP), both of which had short stature. However, these mutations were also reproducibly observed in four ancient control samples from phenotypically healthy individuals (false-positives), rendering the reliable molecular typing of ancient bones for ACH impossible. The treatment of a false-positive DNA extract with uracil N-glycosylase (UNG) to minimize type 2 transitions (G-->A/C-->T) did not reduce the frequency of the false-positive diagnostic ACH mutations. Recently, it was suggested that ancient DNA extracts may induce mutations under PCR. Contemporary human template DNA from a phenotypically healthy individual was therefore spiked with an ancient DNA extract from a cave bear. Again, sequences with the diagnostic G-->A transition in the FGFR3 gene were observed, and it is likely that the false-positive G-->A transitions result from errors introduced during the PCR reaction. Amplifications in the presence of MnCl(2) indicate that position 1138 of the FGFR3 gene is particularly sensitive for mutations. Our data are in line with previously published results on the occurrence of nonrandom mutations in PCR products of contemporary human mitochondrial HVRI template DNA spiked with ancient DNA extracts.