Characterising fitness landscapes of protein evolution by next-generation sequencing

A protein’s amino acid sequence determines its structural, chemical and physical properties, yet how sequence variation influences protein function is still incompletely understood. Protein fitness landscapes powerfully describe the sequence-function relationship by dividing sequence space into functional hills and valleys. This representation is often invoked yet lacks experimental evidence; the immense vastness of possible sequence space makes comprehensive high-quality datasets difficult to obtain. Laboratory directed evolution has focused on optimal utilisation of substitution libraries, however examination of functional innovation in Nature shows that short insertions and deletions (InDels) also play a key role. Beyond rare targeted studies of specific InDels, high-throughput data on fitness landscape for mutations other than substitutions are lacking entirely. In my PhD, I worked towards experimentally describing the fitness landscapes of InDels and substitutions in three systems: GFP, phosphotriesterase (PTE) and the kinase MKK1 docking domain. Towards this goal, I established two experimental assays (GFP, PTE) for deep mutational scanning and a new software toolkit, InDelScanner, for interpreting resulting data that contain InDels. With GFP, I sorted the deletions and substitution libraries into three activity fractions using FACS, then deep sequenced them with Illumina MiSeq to obtain a pilot dataset. The comparison of deletion effects between different lengths of deletions (-3, -6 and -9 bp) indicates that deletions are partially tolerated in eGFP, with tolerance improved for short deletions and in the stabilised starting point GFP8. Further interpretation of data was complicated by limited resolution in the sequencing dataset stemming from poor FACS separation, so I optimised the conditions for better sorting resolution using the mKate2 fluorescent protein as an expression reporter. In the second iteration of the activity sorting I additionally included UMIs in the plasmid design to improve the utilisation of NGS capacity. In the case of PTE, I performed proof-of-concept experiments for microfluidic droplet sorting in an integrated device with an in-line incubation line and a fluorescent sorting design. In parallel, testing of solubility and activity of random InDel variants showed that functional InDels do not necessarily suffer from a stability handicap, making InDel mutagenesis a viable strategy for gene randomisation in directed evolution. One challenge of InDel library data analysis is that InDels are not compatible with existing, substitution-focused software. Using the GFP deletions dataset, I developed the InDelScanner scripts which accurately detect, aggregate and filter insertions, deletions and substitutions. Using the scripts for composition analysis of TRIAD libraries in PTE showed these libraries are well balanced and highly diverse. Finally, I used the InDelScanner scripts to interpret a deep mutational scanning dataset that recorded the sequence preferences in the MKK1 docking domain, acting to activate ERK2. This experiment showed that the fitness landscape in this kinase pair is shaped by the activating effect of hydrophobic residues in the docking groove, as well as widespread positive epistasis. Together, the projects in this thesis demonstrate that deep mutational scanning experiments are a powerful method for exploring the sequence-function relationship in proteins, which can extend into comparison of different types of mutations as well as probing their (epistatic) interactions.BBSRC BB/M011194/

Accessing unexplored regions of sequence space in directed enzyme evolution via insertion/deletion mutagenesis

Abstract: Insertions and deletions (InDels) are frequently observed in natural protein evolution, yet their potential remains untapped in laboratory evolution. Here we introduce a transposon-based mutagenesis approach (TRIAD) to generate libraries of random variants with short in-frame InDels, and screen TRIAD libraries to evolve a promiscuous arylesterase activity in a phosphotriesterase. The evolution exhibits features that differ from previous point mutagenesis campaigns: while the average activity of TRIAD variants is more compromised, a larger proportion has successfully adapted for the activity. Different functional profiles emerge: (i) both strong and weak trade-off between activities are observed; (ii) trade-off is more severe (20- to 35-fold increased kcat/KM in arylesterase with 60-400-fold decreases in phosphotriesterase activity) and (iii) improvements are present in kcat rather than just in KM, suggesting adaptive solutions. These distinct features make TRIAD an alternative to widely used point mutagenesis, accessing functional innovations and traversing unexplored fitness landscape regions

Alcohol-derived DNA crosslinks are repaired by two distinct mechanisms

Acetaldehyde is a highly reactive, DNA-damaging metabolite that is produced upon alcohol consumption1. Impaired detoxification of acetaldehyde is common in the Asian population, and is associated with alcohol-related cancers1,2. Cells are protected against acetaldehyde-induced damage by DNA crosslink repair, which when impaired causes Fanconi anaemia (FA), a disease resulting in failure to produce blood cells and a predisposition to cancer3,4. The combined inactivation of acetaldehyde detoxification and the FA pathway induces mutation, accelerates malignancies and causes the rapid attrition of blood stem cells5,6,7. However, the nature of the DNA damage induced by acetaldehyde and how this is repaired remains a key question. Here we generate acetaldehyde-induced DNA interstrand crosslinks and determine their repair mechanism in Xenopus egg extracts. We find that two replication-coupled pathways repair these lesions. The first is the FA pathway, which operates using excision—analogous to the mechanism used to repair the interstrand crosslinks caused by the chemotherapeutic agent cisplatin. However, the repair of acetaldehyde-induced crosslinks results in increased mutation frequency and an altered mutational spectrum compared with the repair of cisplatin-induced crosslinks. The second repair mechanism requires replication fork convergence, but does not involve DNA incisions—instead the acetaldehyde crosslink itself is broken. The Y-family DNA polymerase REV1 completes repair of the crosslink, culminating in a distinct mutational spectrum. These results define the repair pathways of DNA interstrand crosslinks caused by an endogenous and alcohol-derived metabolite, and identify an excision-independent mechanism

Emerging Links between Microbiome Composition and Skin Immunology in Diaper Dermatitis: A Narrative Review

Diaper dermatitis is a common type of irritant contact dermatitis occurring in infants and toddlers. Its occurrence is triggered by an unfavorable environment under the diaper, damage to skin integrity by fecal enzyme degradation, overhydration and disruption of the lipid bilayer structure facilitating the entry of irritants and microorganisms. In diaper dermatitis development, the central proinflammatory cytokines are IL-1α, IL-8 and TNF-α. The initial release of IL-1α and TNF-α starts a further cascade of pro-inflammatory chemo- and cytokines, resulting in inflammation and erythema of the skin. A recently recognized factor in diaper dermatitis is the composition of the skin microbiome; common pathogenic strains Candida albicans and Staphylococcus aureus are associated with skin irritation. The resulting impaired microbiome composition produces a local inflammatory response and may thus worsen the initial dermatitis clinical presentation and subsequent healing. Introduction of probiotics is an attractive treatment for microbiome modulation, which has shown success in other skin conditions in adults and children. Probiotics are thought to work as a protective shield against irritants, maintain low skin pH, secrete beneficial metabolites, and block pathogen invasion. There is preliminary evidence that certain probiotics given orally or topically could be used as a gentle intervention in diaper dermatitis

Accessing unexplored regions of sequence space in directed enzyme evolution via insertion/deletion mutagenesis

Abstract: Insertions and deletions (InDels) are frequently observed in natural protein evolution, yet their potential remains untapped in laboratory evolution. Here we introduce a transposon-based mutagenesis approach (TRIAD) to generate libraries of random variants with short in-frame InDels, and screen TRIAD libraries to evolve a promiscuous arylesterase activity in a phosphotriesterase. The evolution exhibits features that differ from previous point mutagenesis campaigns: while the average activity of TRIAD variants is more compromised, a larger proportion has successfully adapted for the activity. Different functional profiles emerge: (i) both strong and weak trade-off between activities are observed; (ii) trade-off is more severe (20- to 35-fold increased kcat/KM in arylesterase with 60-400-fold decreases in phosphotriesterase activity) and (iii) improvements are present in kcat rather than just in KM, suggesting adaptive solutions. These distinct features make TRIAD an alternative to widely used point mutagenesis, accessing functional innovations and traversing unexplored fitness landscape regions

Droplet-based screening of phosphate transfer catalysis reveals how epistasis shapes MAP kinase interactions with substrates.

The combination of ultrahigh-throughput screening and sequencing informs on function and intragenic epistasis within combinatorial protein mutant libraries. Establishing a droplet-based, in vitro compartmentalised approach for robust expression and screening of protein kinase cascades (>107 variants/day) allowed us to dissect the intrinsic molecular features of the MKK-ERK signalling pathway, without interference from endogenous cellular components. In a six-residue combinatorial library of the MKK1 docking domain, we identified 29,563 sequence permutations that allow MKK1 to efficiently phosphorylate and activate its downstream target kinase ERK2. A flexibly placed hydrophobic sequence motif emerges which is defined by higher order epistatic interactions between six residues, suggesting synergy that enables high connectivity in the sequence landscape. Through positive epistasis, MKK1 maintains function during mutagenesis, establishing the importance of co-dependent residues in mammalian protein kinase-substrate interactions, and creating a scenario for the evolution of diverse human signalling networks.RCUK | Biotechnology and Biological Sciences Research Council (BBSRC) - BB/M011194/1 EC | Horizon 2020 Framework Programme (EU Framework Programme for Research and Innovation H2020) - 721613 EC | Horizon 2020 Framework Programme (EU Framework Programme for Research and Innovation H2020) - 659029 EC | EU Framework Programme for Research and Innovation H2020 | H2020 Priority Excellent Science | H2020 European Research Council (H2020 Excellent Science - European Research Council) - 69566

Accessing unexplored regions of sequence space in directed enzyme evolution via insertion/deletion mutagenesis

Insertions/Deletions (InDels) remain an untapped source of protein diversity in laboratory evolution. The method TRIAD generates libraries of random variants with short in-frame InDels using transposons, allowing a comparison of their evolutionary potential with widely-used point mutant libraries

Synthesis, molecular docking, and neuroprotective effect of 2-methylcinnamic acid amide in 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP) - an induced Parkinson’s disease model

Parkinson’s disease (PD) has emerged as the second most common form of human neurodegenerative disorders. However, due to the severe side effects of the current antiparkinsonian drugs, the design of novel and safe compounds is a hot topic amongst the medicinal chemistry community. Herein, a convenient peptide method, TBTU (O-(benzotriazole-1-yl)-N,N,N′,N′-tetramethyluronium tetrafluoroborate), was used for the synthesis of the amide (E)-N-(2-methylcinnamoyl)-amantadine (CA(2-Me)-Am3)) derived from amantadine and 2-methylcinnamic acid. The obtained hybrid was studied for its antiparkinsonian activity in an experimental model of PD induced by MPTP. Mice (C57BL/6,male, 8 weeks old) were divided into four groups as follows: (1) the control, treated with normal saline (i.p.) for 12 consecutive days(2) MPTP (30 mg/kg/day, i.p.), applied daily for 5 consecutive days(3) MPTP + CA(2-Me)-Am, applied for 12 consecutive days, 5 days simultaneously with MPTP and 7 days after MPTP(4) CA(2-Me)-Am +oleanoic acid (OA), applied daily for 12 consecutive days. Neurobehavioral parameters in all experimental groups of mice were evaluated by rotarod test and passive avoidance test. Our experimental data showed that CA(2-Me)-Am in parkinsonian mice significantly restored memory performance, while neuromuscular coordination approached the control level, indicating the ameliorating effects of the new compound. In conclusion, the newly synthesized hybrid might be a promising agent for treating motor disturbances and cognitive impairment in experimental PD

Synthesis, Molecular Docking, and Neuroprotective Effect of 2-Methylcinnamic Acid Amide in 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP)—An Induced Parkinson’s Disease Model

Parkinson’s disease (PD) has emerged as the second most common form of human neurodegenerative disorders. However, due to the severe side effects of the current antiparkinsonian drugs, the design of novel and safe compounds is a hot topic amongst the medicinal chemistry community. Herein, a convenient peptide method, TBTU (O-(benzotriazole-1-yl)-N,N,N′,N′-tetramethyluronium tetrafluoroborate), was used for the synthesis of the amide (E)-N-(2-methylcinnamoyl)-amantadine (CA(2-Me)-Am; 3)) derived from amantadine and 2-methylcinnamic acid. The obtained hybrid was studied for its antiparkinsonian activity in an experimental model of PD induced by MPTP. Mice (C57BL/6,male, 8 weeks old) were divided into four groups as follows: (1) the control, treated with normal saline (i.p.) for 12 consecutive days; (2) MPTP (30 mg/kg/day, i.p.), applied daily for 5 consecutive days; (3) MPTP + CA(2-Me)-Am, applied for 12 consecutive days, 5 days simultaneously with MPTP and 7 days after MPTP; (4) CA(2-Me)-Am +oleanoic acid (OA), applied daily for 12 consecutive days. Neurobehavioral parameters in all experimental groups of mice were evaluated by rotarod test and passive avoidance test. Our experimental data showed that CA(2-Me)-Am in parkinsonian mice significantly restored memory performance, while neuromuscular coordination approached the control level, indicating the ameliorating effects of the new compound. In conclusion, the newly synthesized hybrid might be a promising agent for treating motor disturbances and cognitive impairment in experimental PD