Evolution of protein function by domain swapping

Engineering of a Baeyer-Villiger monooxygenase to Improve Substrate Scope, Stereoselectivity and Regioselectivity

Baeyer-Villiger monooxygenases belong to a family of flavin-binding proteins that catalyze the Baeyer-Villiger (BV) oxidation of ketones to produce lactones or esters, which are important intermediates in pharmaceuticals or sustainable materials. Phenylacetone monooxygenase (PAMO) from Thermobifida fusca with moderate thermostability catalyzes the oxidation of aryl ketone substrates, but is limited by high specificity and narrow substrate scope. In the present study, we applied loop optimization by loop swapping followed by focused saturation mutagenesis in order to evolve PAMO mutants capable of catalyzing the regioselective BV oxidation of cyclohexanone and cyclobutanone derivatives with formation of either normal or abnormal esters or lactones. We further modulated PAMO to increase enantioselectivity. Crystal structure studies indicate that rotation occurs in the NADP-binding domain and that the high B-factor region is predominantly distributed in the catalytic pocket residues. Computational analyses further revealed dynamic character in the catalytic pocket and reshaped hydrogen bond interaction networks, which is more favorable for substrate binding. Our study provides useful insights for studying enzyme-substrate adaptations.

Conifers Concentrate Large Numbers of NLR Immune Receptor Genes on One Chromosome

Nucleotide-binding domain and leucine-rich repeat (NLR) immune receptor genes form a major line of defense in plants, acting in both pathogen recognition and resistance machinery activation. NLRs are reported to form large gene clusters in limber pine (Pinus flexilis), but it is unknown how widespread this genomic architecture may be among the extant species of conifers (Pinophyta). We used comparative genomic analyses to assess patterns in the abundance, diversity, and genomic distribution of NLR genes. Chromosome-level whole genome assemblies and high-density linkage maps in the Pinaceae, Cupressaceae, Taxaceae, and other gymnosperms were scanned for NLR genes using existing and customized pipelines. The discovered genes were mapped across chromosomes and linkage groups and analyzed phylogenetically for evolutionary history. Conifer genomes are characterized by dense clusters of NLR genes, highly localized on one chromosome. These clusters are rich in TNL-encoding genes, which seem to have formed through multiple tandem duplication events. In contrast to angiosperms and nonconiferous gymnosperms, genomic clustering of NLR genes is ubiquitous in conifers. NLR-dense genomic regions are likely to influence a large part of the plant's resistance, informing our understanding of adaptation to biotic stress and the development of genetic resources through breeding.

Synthaser: a CD-Search enabled Python toolkit for analysing domain architecture of fungal secondary metabolite megasynth(et)ases

BACKGROUND

Fungi are prolific producers of secondary metabolites (SMs), which are bioactive small molecules with important applications in medicine, agriculture and other industries. The backbones of a large proportion of fungal SMs are generated through the action of large, multi-domain megasynth(et)ases such as polyketide synthases (PKSs) and nonribosomal peptide synthetases (NRPSs). The structure of these backbones is determined by the domain architecture of the corresponding megasynth(et)ase, and thus accurate annotation and classification of these architectures is an important step in linking SMs to their biosynthetic origins in the genome.

RESULTS

Here we report synthaser, a Python package leveraging the NCBI's conserved domain search tool for remote prediction and classification of fungal megasynth(et)ase domain architectures. Synthaser is capable of batch sequence analysis, and produces rich textual output and interactive visualisations which allow for quick assessment of the megasynth(et)ase diversity of a fungal genome. Synthaser uses a hierarchical rule-based classification system, which can be extensively customised by the user through a web application ( http://gamcil.github.io/synthaser ). We show that synthaser provides more accurate domain architecture predictions than comparable tools which rely on curated profile hidden Markov model (pHMM)-based approaches; the utilisation of the NCBI conserved domain database also allows for significantly greater flexibility compared to pHMM approaches. In addition, we demonstrate how synthaser can be applied to large scale genome mining pipelines through the construction of an Aspergillus PKS similarity network.

CONCLUSIONS

Synthaser is an easy to use tool that represents a significant upgrade to previous domain architecture analysis tools. It is freely available under a MIT license from PyPI ( https://pypi.org/project/synthaser ) and GitHub ( https://github.com/gamcil/synthaser ).

CirPred, the first structure modeling and linker design system for circularly permuted proteins

Background

This work aims to help develop new protein engineering techniques based on a structural rearrangement phenomenon called circular permutation (CP), equivalent to connecting the native termini of a protein followed by creating new termini at another site. Although CP has been applied in many fields, its implementation is still costly because of inevitable trials and errors.

Results

Here we present CirPred, a structure modeling and termini linker design method for circularly permuted proteins. Compared with state-of-the-art protein structure modeling methods, CirPred is the only one fully capable of both circularly-permuted modeling and traditional co-linear modeling. CirPred performs well when the permutant shares low sequence identity with the native protein and even when the permutant adopts a different conformation from the native protein because of three-dimensional (3D) domain swapping. Linker redesign experiments demonstrated that the linker design algorithm of CirPred achieved subangstrom accuracy.

Conclusions

The CirPred system is capable of (1) predicting the structure of circular permutants, (2) designing termini linkers, (3) performing traditional co-linear protein structure modeling, and (4) identifying the CP-induced occurrence of 3D domain swapping. This method is supposed helpful for broadening the application of CP, and its web server is available at http://10.life.nctu.edu.tw/CirPred/ and http://lo.life.nctu.edu.tw/CirPred/.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04403-1.

frDriver: A Functional Region Driver Identification for Protein Sequence

Identifying cancer drivers is a crucial challenge to explain the underlying mechanisms of cancer development. There are many methods to identify cancer drivers based on the single mutation site or the entire gene. But they ignore a large number of functional elements with medium in size. It is hypothesized that mutations occurring in different regions of the protein sequence have different effects on the progression of cancer. Here, we develop a novel functional region driver(frDriver) identification method based on Bayesian probability and multiple linear regression models to identify protein regions that can regulate gene expression levels and have high functional impact potential. Combining gene expression data and somatic mutation data, with functional impact scores(SIFT, PROVEAN) as a priori knowledge, we identified cancer driver regions that are most accurate in predicting gene expression levels. We evaluated the performance of frDriver on the BRCA and GBM datasets from TCGA. The results showed that frDriver identified known cancer drivers and outperformed the other three state-of-the-art methods(eDriver, ActiveDriver and OncodriveCLUST). In addition, we performed KEGG pathway and GO term enrichment analysis, and the results indicated that the cancer drivers predicted by frDriver were related to processes such as cancer formation and gene regulation.

Structurally distributed surface sites tune allosteric regulation

Our ability to rationally optimize allosteric regulation is limited by incomplete knowledge of the mutations that tune allostery. Are these mutations few or abundant, structurally localized or distributed? To examine this, we conducted saturation mutagenesis of a synthetic allosteric switch in which Dihydrofolate reductase (DHFR) is regulated by a blue-light sensitive LOV2 domain. Using a high-throughput assay wherein DHFR catalytic activity is coupled to E. coli growth, we assessed the impact of 1548 viable DHFR single mutations on allostery. Despite most mutations being deleterious to activity, fewer than 5% of mutations had a statistically significant influence on allostery. Most allostery disrupting mutations were proximal to the LOV2 insertion site. In contrast, allostery enhancing mutations were structurally distributed and enriched on the protein surface. Combining several allostery enhancing mutations yielded near-additive improvements to dynamic range. Our results indicate a path toward optimizing allosteric function through variation at surface sites.

Many proteins exhibit a property called ‘allostery’. In allostery, an input signal at a specific site of a protein – such as a molecule binding, or the protein absorbing a photon of light – leads to a change in output at another site far away. For example, the protein might catalyze a chemical reaction faster or bind to another molecule more tightly in the presence of the input signal. This protein ‘remote control’ allows cells to sense and respond to changes in their environment. An ability to rapidly engineer new allosteric mechanisms into proteins is much sought after because this would provide an approach for building biosensors and other useful tools. One common approach to engineering new allosteric regulation is to combine a ‘sensor’ or input region from one protein with an ‘output’ region or domain from another.

When researchers engineer allostery using this approach of combining input and output domains from different proteins, the difference in the output when the input is ‘on’ versus ‘off’ is often small, a situation called ‘modest allostery’. McCormick et al. wanted to know how to optimize this domain combination approach to increase the difference in output between the ‘on’ and ‘off’ states.

More specifically, McCormick et al. wanted to find out whether swapping out or mutating specific amino acids (each of the individual building blocks that make up a protein) enhances or disrupts allostery. They also wanted to know if there are many possible mutations that change the effectiveness of allostery, or if this property is controlled by just a few amino acids. Finally, McCormick et al. questioned where in a protein most of these allostery-tuning mutations were located.

To answer these questions, McCormick et al. engineered a new allosteric protein by inserting a light-sensing domain (input) into a protein involved in metabolism (a metabolic enzyme that produces a biomolecule called a tetrahydrofolate) to yield a light-controlled enzyme. Next, they introduced mutations into both the ‘input’ and ‘output’ domains to see where they had a greater effect on allostery.

After filtering out mutations that destroyed the function of the output domain, McCormick et al. found that only about 5% of mutations to the ‘output’ domain altered the allosteric response of their engineered enzyme. In fact, most mutations that disrupted allostery were found near the site where the ‘input’ domain was inserted, while mutations that enhanced allostery were sprinkled throughout the enzyme, often on its protein surface. This was surprising in light of the commonly-held assumption that mutations on protein surfaces have little impact on the activity of the ‘output’ domain. Overall, the effect of individual mutations on allostery was small, but McCormick et al. found that these mutations can sometimes be combined to yield larger effects.

McCormick et al.’s results suggest a new approach for optimizing engineered allosteric proteins: by introducing mutations on the protein surface. It also opens up new questions: mechanically, how do surface sites affect allostery? In the future, it will be important to characterize how combinations of mutations can optimize allosteric regulation, and to determine what evolutionary trajectories to high performance allosteric ‘switches’ look like.

Blueprints for Biosensors: Design, Limitations, and Applications

Biosensors are enabling major advances in the field of analytics that are both facilitating and being facilitated by advances in synthetic biology. The ability of biosensors to rapidly and specifically detect a wide range of molecules makes them highly relevant to a range of industrial, medical, ecological, and scientific applications. Approaches to biosensor design are as diverse as their applications, with major biosensor classes including nucleic acids, proteins, and transcription factors. Each of these biosensor types has advantages and limitations based on the intended application, and the parameters that are required for optimal performance. Specifically, the choice of biosensor design must consider factors such as the ligand specificity, sensitivity, dynamic range, functional range, mode of output, time of activation, ease of use, and ease of engineering. This review discusses the rationale for designing the major classes of biosensor in the context of their limitations and assesses their suitability to different areas of biotechnological application.

Structural entanglements in protein complexes

We consider multi-chain protein native structures and propose a criterion that determines whether two chains in the system are entangled or not. The criterion is based on the behavior observed by pulling at both termini of each chain simultaneously in the two chains. We have identified about 900 entangled systems in the Protein Data Bank and provided a more detailed analysis for several of them. We argue that entanglement enhances the thermodynamic stability of the system but it may have other functions: burying the hydrophobic residues at the interface and increasing the DNA or RNA binding area. We also study the folding and stretching properties of the knotted dimeric proteins MJ0366, YibK, and bacteriophytochrome. These proteins have been studied theoretically in their monomeric versions so far. The dimers are seen to separate on stretching through the tensile mechanism and the characteristic unraveling force depends on the pulling direction.

Dual binding in cohesin-dockerin complexes: the energy landscape and the role of short, terminal segments of the dockerin module

The assembly of the polysaccharide degradating cellulosome machinery is mediated by tight binding between cohesin and dockerin domains. We have used an empirical model known as FoldX as well as molecular mechanics methods to determine the free energy of binding between a cohesin and a dockerin from Clostridium thermocellum in two possible modes that differ by an approximately 180° rotation. Our studies suggest that the full-length wild-type complex exhibits dual binding at room temperature, i.e., the two modes of binding have comparable probabilities at equilibrium. The ability to bind in the two modes persists at elevated temperatures. However, single-point mutations or truncations of terminal segments in the dockerin result in shifting the equilibrium towards one of the binding modes. Our molecular dynamics simulations of mechanical stretching of the full-length wild-type cohesin-dockerin complex indicate that each mode of binding leads to two kinds of stretching pathways, which may be mistakenly taken as evidence of dual binding.

Exploring modular allostery via interchangeable regulatory domains

Most proteins comprise two or more domains from a limited suite of protein families. These domains are often rearranged in various combinations through gene fusion events to evolve new protein functions, including the acquisition of protein allostery through the incorporation of regulatory domains. The enzyme 3-deoxy-d-arabino-heptulosonate 7-phosphate synthase (DAH7PS) is the first enzyme of aromatic amino acid biosynthesis and displays a diverse range of allosteric mechanisms. DAH7PSs adopt a common architecture with a shared (β/α)₈ catalytic domain which can be attached to an ACT-like or a chorismate mutase regulatory domain that operates via distinct mechanisms. These respective domains confer allosteric regulation by controlling DAH7PS function in response to ligand Tyr or prephenate. Starting with contemporary DAH7PS proteins, two protein chimeras were created, with interchanged regulatory domains. Both engineered proteins were catalytically active and delivered new functional allostery with switched ligand specificity and allosteric mechanisms delivered by their nonhomologous regulatory domains. This interchangeability of protein domains represents an efficient method not only to engineer allostery in multidomain proteins but to create a new bifunctional enzyme.

A combination of mutational and computational scanning guides the design of an artificial ligand-binding controlled lipase

Allostery, i.e. the control of enzyme activity by a small molecule at a location distant from the enzyme’s active site, represents a mechanism essential for sustaining life. The rational design of allostery is a non-trivial task but can be achieved by fusion of a sensory domain, which responds to environmental stimuli with a change in its structure. Hereby, the site of domain fusion is difficult to predict. We here explore the possibility to rationally engineer allostery into the naturally not allosterically regulated Bacillus subtilis lipase A, by fusion of the citrate-binding sensor-domain of the CitA sensory-kinase of Klebsiella pneumoniae. The site of domain fusion was rationally determined based on whole-protein site-saturation mutagenesis data, complemented by computational evolutionary-coupling analyses. Functional assays, combined with biochemical and biophysical studies suggest a mechanism for control, similar but distinct to the one of the parent CitA protein, with citrate acting as an indirect modulator of Triton-X100 inhibition of the fusion protein. Our study demonstrates that the introduction of ligand-dependent regulatory control by domain fusion is surprisingly facile, suggesting that the catalytic mechanism of some enzymes may be evolutionary optimized in a way that it can easily be perturbed by small conformational changes.

Linker Flexibility Facilitates Module Exchange in Fungal Hybrid PKS-NRPS Engineering

Polyketide synthases (PKSs) and nonribosomal peptide synthetases (NRPSs) each give rise to a vast array of complex bioactive molecules with further complexity added by the existence of natural PKS-NRPS fusions. Rational genetic engineering for the production of natural product derivatives is desirable for the purpose of incorporating new functionalities into pre-existing molecules, or for optimization of known bioactivities. We sought to expand the range of natural product diversity by combining modules of PKS-NRPS hybrids from different hosts, hereby producing novel synthetic natural products. We succeeded in the construction of a functional cross-species chimeric PKS-NRPS expressed in Aspergillus nidulans. Module swapping of the two PKS-NRPS natural hybrids CcsA from Aspergillus clavatus involved in the biosynthesis of cytochalasin E and related Syn2 from rice plant pathogen Magnaporthe oryzae lead to production of novel hybrid products, demonstrating that the rational re-design of these fungal natural product enzymes is feasible. We also report the structure of four novel pseudo pre-cytochalasin intermediates, niduclavin and niduporthin along with the chimeric compounds niduchimaeralin A and B, all indicating that PKS-NRPS activity alone is insufficient for proper assembly of the cytochalasin core structure. Future success in the field of biocombinatorial synthesis of hybrid polyketide-nonribosomal peptides relies on the understanding of the fundamental mechanisms of inter-modular polyketide chain transfer. Therefore, we expressed several PKS-NRPS linker-modified variants. Intriguingly, the linker anatomy is less complex than expected, as these variants displayed great tolerance with regards to content and length, showing a hitherto unreported flexibility in PKS-NRPS hybrids, with great potential for synthetic biology-driven biocombinatorial chemistry.

BluB/CobT2 fusion enzyme activity reveals mechanisms responsible for production of active form of vitamin B₁₂ by Propionibacterium freudenreichii

Background

Propionibacterium freudenreichii is a food grade bacterium that has gained attention as a producer of appreciable amounts of cobalamin, a cobamide with activity of vitamin B₁₂. Production of active form of vitamin is a prerequisite for attempts to naturally fortify foods with B₁₂ by microbial fermentation. Active vitamin B₁₂ is distinguished from the pseudovitamin by the presence of 5,6-dimethylbenzimidazole (DMBI) as the lower ligand. Genomic data indicate that P. freudenreichii possesses a fusion gene, bluB/cobT2, coding for a predicted phosphoribosyltransferase/nitroreductase, which is presumably involved in production of vitamin B₁₂. Understanding the mechanisms affecting the synthesis of different vitamin forms is useful for rational strain selection and essential for engineering of strains with improved B₁₂ production properties.

Results

Here, we investigated the activity of heterologously expressed and purified fusion enzyme BluB/CobT2. Our results show that BluB/CoBT2 is responsible for the biosynthesis of the DMBI base and its activation into α-ribazole phosphate, preparing it for attachment as the lower ligand of cobalamin. The fusion enzyme was found to be efficient in metabolite channeling and the enzymes’ inability to react with adenine, a lower ligand present in the pseudovitamin, revealed a mechanism favoring the production of the active form of the vitamin. P. freudenreichii did not produce cobalamin under strictly anaerobic conditions, confirming the requirement of oxygen for DMBI synthesis. In vivo experiments also revealed a clear preference for incorporating DMBI over adenine into cobamide under both microaerobic and anaerobic conditions.

Conclusions

The herein described BluB/CobT2 is responsible for the production and activation of DMBI. Fusing those two activities results in high pressure towards production of the true vitamin B₁₂ by efficiently activating DMBI formed within the same enzymatic complex. This indicates that BluB/CobT2 is the crucial enzyme in the B₁₂ biosynthetic pathway of P. freudenreichii. The GRAS organism status and the preference for synthesizing active vitamin form make P. freudenreichii a unique candidate for the in situ production of vitamin B₁₂ within food products.

Electronic supplementary material

The online version of this article (doi:10.1186/s12934-015-0363-9) contains supplementary material, which is available to authorized users.

Design of hyperthermophilic lipase chimeras by key motif-directed recombination

Recombination of diverse natural evolved domains within a superfamily offers greater opportunity for enzyme function leaps. How to recombine protein modules from distant parents with less disruption in cross-interfaces is a challenging issue. Here, we identified the existence of a key motif, the sequence VVSVN(D)YR, within a structural motif ψ loop in the α/β-hydrolase fold superfamily, by using a MEME server and the PROMOTIF program. To obtain thermostable lipase-like enzymes, two chimeras were engineered at the key motif regions through recombination of domains from a mesophilic lipase and a hyperthermophilic esterase/peptidase with amino acid identity less than 21 %. The chimeras retained the desirable substrate preference of their mesophilic parent and exhibited more than 100-fold increased thermostability at 50 °C. Through site-directed mutation, we further improved activity of the chimera by 4.6-fold. The recombination strategy presented here enables the creation of novel catalysts.

C-terminal domain swapping of SSB changes the size of the ssDNA binding site

Single-stranded DNA-binding protein (SSB) plays an important role in DNA metabolism, including DNA replication, repair, and recombination, and is therefore essential for cell survival. Bacterial SSB consists of an N-terminal ssDNA-binding/oligomerization domain and a flexible C-terminal protein-protein interaction domain. We characterized the ssDNA-binding properties of Klebsiella pneumoniae SSB (KpSSB), Salmonella enterica Serovar Typhimurium LT2 SSB (StSSB), Pseudomonas aeruginosa PAO1 SSB (PaSSB), and two chimeric KpSSB proteins, namely, KpSSBnStSSBc and KpSSBnPaSSBc. The C-terminal domain of StSSB or PaSSB was exchanged with that of KpSSB through protein chimeragenesis. By using the electrophoretic mobility shift assay, we characterized the stoichiometry of KpSSB, StSSB, PaSSB, KpSSBnStSSBc, and KpSSBnPaSSBc, complexed with a series of ssDNA homopolymers. The binding site sizes were determined to be 26 ± 2, 21 ± 2, 29 ± 2, 21 ± 2, and 29 ± 2 nucleotides (nt), respectively. Comparison of the binding site sizes of KpSSB, KpSSBnStSSBc, and KpSSBnPaSSBc showed that the C-terminal domain swapping of SSB changes the size of the binding site. Our observations suggest that not only the conserved N-terminal domain but also the C-terminal domain of SSB is an important determinant for ssDNA binding.

Uncoupling of allosteric and oligomeric regulation in a functional hybrid enzyme constructed from Escherichia coli and human ribonucleotide reductase

An N-terminal-domain (NTD) and adjacent catalytic body (CB) make up subunit-α of ribonucleotide reductase (RNR), the rate-limiting enzyme for de novo dNTP biosynthesis. A strong linkage exists between ligand binding at the NTD and oligomerization-coupled RNR inhibition, inducible by both dATP and nucleotide chemotherapeutics. These observations have distinguished the NTD as an oligomeric regulation domain dictating the assembly of inactive RNR oligomers. Inactive states of RNR differ between eukaryotes and prokaryotes (α6 in human versus α4β4 in Escherichia coli , wherein β is RNR's other subunit); however, the NTD structurally interconnects individual α2 or α2 and β2 dimeric motifs within the respective α6 or α4β4 complexes. To elucidate the influence of NTD ligand binding on RNR allosteric and oligomeric regulation, we engineered a human- E. coli hybrid enzyme (HE) where human-NTD is fused to E. coli -CB. Both the NTD and the CB of the HE bind dATP. The HE specifically partners with E. coli -β to form an active holocomplex. However, although the NTD is the sole physical tether to support α2 and/or β2 associations in the dATP-bound α6 or α4β4 fully inhibited RNR complexes, the binding of dATP to the HE NTD only partially suppresses HE activity and fully precludes formation of higher-order HE oligomers. We postulate that oligomeric regulation is the ultimate mechanism for potent RNR inhibition, requiring species-specific NTD-CB interactions. Such interdomain cooperativity in RNR oligomerization is unexpected from structural studies alone or biochemical studies of point mutants.

Dramatic structural changes resulting from the loss of a crucial hydrogen bond in the hinge region involved in C-terminal helix swapping in SurE: a survival protein from Salmonella typhimurium

Domain swapping is an interesting feature of some oligomeric proteins in which each protomer of the oligomer provides an identical surface for exclusive interaction with a segment or domain belonging to another protomer. Here we report results of mutagenesis experiments on the structure of C-terminal helix swapped dimer of a stationary phase survival protein from Salmonella typhimurium (StSurE). Wild type StSurE is a dimer in which a large helical segment at the C-terminus and a tetramerization loop comprising two β strands are swapped between the protomers. Key residues in StSurE that might promote C-terminal helix swapping were identified by sequence and structural comparisons. Three mutants in which the helix swapping is likely to be avoided were constructed and expressed in E. coli. Three-dimensional X-ray crystal structures of the mutants H234A and D230A/H234A could be determined at 2.1 Å and 2.35 Å resolutions, respectively. Contrary to expectations, helix swapping was mostly retained in both the mutants. The loss of the crucial D230 OD2– H234 NE2 hydrogen bond (2.89 Å in the wild type structure) in the hinge region was compensated by new inter and intra-chain interactions. However, the two fold molecular symmetry was lost and there were large conformational changes throughout the polypeptide. In spite of these changes, the dimeric structure and an approximate tetrameric organization were retained, probably due to the interactions involving the tetramerization loop. Mutants were mostly functionally inactive, highlighting the importance of precise inter-subunit interactions for the symmetry and function of StSurE.

Engineering allosteric control to an unregulated enzyme by transfer of a regulatory domain

Allosteric regulation of protein function is a critical component of metabolic control. Its importance is underpinned by the diversity of mechanisms and its presence in all three domains of life. The first enzyme of the aromatic amino acid biosynthesis, 3-deoxy-D-arabino-heptulosonate 7-phosphate synthase, shows remarkable variation in allosteric response and machinery, and both contemporary regulated and unregulated orthologs have been described. To examine the molecular events by which allostery can evolve, we have generated a chimeric protein by joining the catalytic domain of an unregulated 3-deoxy-D-arabino-heptulosonate 7-phosphate synthase with the regulatory domain of a regulated enzyme. We demonstrate that this simple gene fusion event on its own is sufficient to confer functional allostery to the unregulated enzyme. The fusion protein shares structural similarities with its regulated parent protein and undergoes an analogous major conformational change in response to the binding of allosteric effector tyrosine to the regulatory domain. These findings help delineate a remarkably facile mechanism for the evolution of modular allostery by domain recruitment.

Structure-based optimization of designed Armadillo-repeat proteins

The armadillo domain is a right-handed super-helix of repeating units composed of three α-helices each. Armadillo repeat proteins (ArmRPs) are frequently involved in protein-protein interactions, and because of their modular recognition of extended peptide regions they can serve as templates for the design of artificial peptide binding scaffolds. On the basis of sequential and structural analyses, different consensus-designed ArmRPs were synthesized and show high thermodynamic stabilities, compared to naturally occurring ArmRPs. We determined the crystal structures of four full-consensus ArmRPs with three or four identical internal repeats and two different designs for the N- and C-caps. The crystal structures were refined at resolutions ranging from 1.80 to 2.50 Å for the above mentioned designs. A redesign of our initial caps was required to obtain well diffracting crystals. However, the structures with the redesigned caps caused domain swapping events between the N-caps. To prevent this domain swap, 9 and 6 point mutations were introduced in the N- and C-caps, respectively. Structural and biophysical analysis showed that this subsequent redesign of the N-cap prevented domain swapping and improved the thermodynamic stability of the proteins. We systematically investigated the best cap combinations. We conclude that designed ArmRPs with optimized caps are intrinsically stable and well-expressed monomeric proteins and that the high-resolution structures provide excellent structural templates for the continuation of the design of sequence-specific modular peptide recognition units based on armadillo repeats.

Phosphorothioate-based DNA recombination: an enzyme-free method for the combinatorial assembly of multiple DNA fragments

Rational guided generation of protein chimeras has developed into an attractive approach in protein engineering for tailoring catalytic and biophysical properties of enzymes. Combinatorial recombination of structural elements or whole protein domains is still technically challenging due to sequence dependent biases diminishing the overall quality of resulting chimeric libraries. Since methods for generating such libraries are often limited by a low frequency of crossover points and suffer from challenges in handling, we developed the phosphorothioate-based DNA recombination method (PTRec). PTRec is an enzyme-free method and only requires a short stretch of four amino acids that is identical among the proteins to be recombined in order to define a single crossover point. In a PTRec-generated chimeric library that shuffled five domains of phytase using genes from three different species, 88% of 42 randomly picked and sequenced genes were efficiently recombined. Furthermore, PTRec is a technically simple, fast, and reliable method that can be used for domain-and exon-shuffling or recombination of DNA fragments in general.

Functional repertoire, molecular pathways and diseases associated with 3D domain swapping in the human proteome

Background

3D domain swapping is a novel structural phenomenon observed in diverse set of protein structures in oligomeric conformations. A distinct structural feature, where structural segments in a protein dimer or higher oligomer were shared between two or more chains of a protein structure, characterizes 3D domain swapping. 3D domain swapping was observed as a key mediator of numerous functional mechanisms and play pathogenic role in various diseases including conformational diseases like amyloidosis, Alzheimer's disease, Parkinson's disease and prion diseases. We report the first study with a focus on identifying functional classes, pathways and diseases mediated by 3D domain swapping in the human proteome.

Methods

We used a panel of four enrichment tools with two different ontologies and two annotations database to derive biological and clinical relevant information associated with 3D domain swapping. Protein domain enrichment analysis followed by Gene Ontology (GO) term enrichment analysis revealed the functional repertoire of proteins involved in swapping. Pathway analysis using KEGG annotations revealed diverse pathway associations of human proteins involved in 3D domain swapping. Disease Ontology was used to find statistically significant associations with proteins in swapped conformation and various disease categories (P-value < 0.05).

Results

We report meta-analysis results of a literature-curated dataset of human gene products involved in 3D domain swapping and discuss new insights about the functional repertoire, pathway associations and disease implications of proteins involved in 3D domain swapping.

Conclusions

Our integrated bioinformatics pipeline comprising of four different enrichment tools, two ontologies and two annotations revealed new insights into the functional and disease correlations with 3D domain swapping. GO term enrichment were used to infer terms associated with three different GO categories. Protein domain enrichment was used to identify conserved domains enriched in swapped proteins. Pathway enrichment analysis using KEGG annotations revealed that proteins with swapped conformations are present in all six classes of KEGG BRITE hierarchy and significantly enriched KEGG pathways were observed in five classes. Five major classes of disease were found to be associated with 3D domain swapping using functional disease ontology based enrichment analysis. Five classes of human diseases: cancer, diseases of the respiratory or pulmonary system, degenerative diseases of the central nervous system, vascular disease and encephalitis were found to be significant. In conclusion, our study shows that bioinformatics based analytical approaches using curated data can enhance the understanding of functional and disease implications of 3D domain swapping.

Analysis of osmotin, a PR protein as metabolic modulator in plants

Osmotin is an abundant cationic multifunctional protein discovered in cells of tobacco (Nicotiana tabacum L. var Wisconsin 38) adapted to an environment of low osmotic potential. Beside its role as osmoregulator, it provides plants protection from pathogens, hence also placed in the PRP family of proteins. The osmotin induced proline accumulation has been reported to confer tolerance against both biotic and abiotic stresses in plants including transgenic tomato and strawberry overexpressing osmotin gene. The exact mechanism of induction of proline by osmotin is however, not known till date. These observations have led us to hypothesize that osmotin could be regulating these plant responses through its involvement either as transcription factor, cell signal pathway modulator or both in plants. We have therefore, undertaken the present investigation to analyze the osmotin protein as transcription factor using bioinformatics tools. The results of available online DNA binding motif search programs revealed that osmotin does not contain DNAbinding motifs. The alignment results of osmotin protein with the protein sequence from DATF showed the homology in the range of 0-20%, suggesting that it might not contain a DNA binding motif. Further to find unique DNA-binding domain, the superimposition of osmotin 3D structure on modeled Arabidopsis transcription factors using Chimera also suggested absence of the same. However, evidence implicating osmotin in cell signaling were found during the study. With these results, we therefore, concluded that osmotin is not a transcription factor, but regulating plant responses to biotic and abiotic stresses through cell signaling.

Exploring and exploiting allostery: Models, evolution, and drug targeting

The concept of allostery was elaborated almost 50years ago by Monod and coworkers to provide a framework for interpreting experimental studies on the regulation of protein function. In essence, binding of a ligand at an allosteric site affects the function at a distant site exploiting protein flexibility and reshaping protein energy landscape. Both monomeric and oligomeric proteins can be allosteric. In the past decades, the behavior of allosteric systems has been analyzed in many investigations while general theoretical models and variations thereof have been steadily proposed to interpret the experimental data. Allostery has been established as a fundamental mechanism of regulation in all organisms, governing a variety of processes that range from metabolic control to receptor function and from ligand transport to cell motility. A number of studies have shed light on how evolutionary pressures have favored and molded the development of allosteric features in specific macromolecular systems. The widespread occurrence of allostery has been recently exploited for the development and design of allosteric drugs that bind to either physiological or non-physiological allosteric sites leading to gain of function or loss of function. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

USER friendly DNA recombination (USERec): a simple and flexible near homology-independent method for gene library construction

USER friendly DNA recombination (USERec) is introduced as a near homology-independent method that allows the simultaneous recombination of an unprecedented number of 10 DNA fragments (approximately 40-400 bp) within a day. The large number of fragments and their ease of preparation enables the creation of libraries of much larger genetic diversity (potentially approximately 10(10)-10(11) sequences) than current alternative methods based on DNA truncation (ITCHY, SCRATCHY and SHIPREC) or type IIb restriction enzymes (SISDC). At the same time, the frequency of frameshifts in the recombined library is low (90% of the recombined sequences are in frame). Compared to overlap extension PCR, USERec also requires much reduced crossover sequence constraints (only a 5'-AN(4-8)T-3' motif) and fewer experimental steps. Based on its simplicity and flexibility, and the accessibility of large and high quality recombined DNA libraries, USERec is established as a convenient alternative for the combinatorial assembly of gene fragments (e.g. exon or domain shuffling) and for a number of applications in gene library construction, such as loop grafting and multi-site-directed or random mutagenesis.

Molecular mechanisms of recombination restriction in the envelope gene of the human immunodeficiency virus

The ability of pathogens to escape the host's immune response is crucial for the establishment of persistent infections and can influence virulence. Recombination has been observed to contribute to this process by generating novel genetic variants. Although distinctive recombination patterns have been described in many viral pathogens, little is known about the influence of biases in the recombination process itself relative to selective forces acting on newly formed recombinants. Understanding these influences is important for determining how recombination contributes to pathogen genome and proteome evolution. Most previous research on recombination-driven protein evolution has focused on relatively simple proteins, usually in the context of directed evolution experiments. Here, we study recombination in the envelope gene of HIV-1 between primary isolates belonging to subtypes that recombine naturally in the HIV/AIDS pandemic. By characterizing the early steps in the generation of recombinants, we provide novel insights into the evolutionary forces that shape recombination patterns within viral populations. Specifically, we show that the combined effects of mechanistic processes that determine the locations of recombination breakpoints across the HIV-1 envelope gene, and purifying selection acting against dysfunctional recombinants, can explain almost the entire distribution of breakpoints found within this gene in nature. These constraints account for the surprising paucity of recombination breakpoints found in infected individuals within this highly variable gene. Thus, the apparent randomness of HIV evolution via recombination may in fact be relatively more predictable than anticipated. In addition, the dominance of purifying selection in localized areas of the HIV genome defines regions where functional constraints on recombinants appear particularly strong, pointing to vulnerable aspects of HIV biology.

Recombination allows mixing portions of genomes of different origins, generating chimeric genes and genomes. With respect to the random generation of new mutations, it can lead to the simultaneous insertion of several substitutions, introducing more drastic changes in the genome. Furthermore, recombination is expected to yield a higher proportion of functional products since it combines variants that already exist in the population and that are therefore compatible with the survival of the organism. However, when recombination involves genetically distant strains, it can be constrained by the necessity to retain the functionality of the resulting products. In pathogens, which are subjected to strong selective pressures, recombination is particularly important, and several viruses, such as the human immunodeficiency virus (HIV), readily recombine. Here, we demonstrate the existence of preferential regions for recombination in the HIV-1 envelope gene when crossing sequences representative of strains observed to recombine in vivo. Furthermore, some recombinants give a decreased proportion of functional products. When considering these factors, one can retrace the history of most natural HIV recombinants. Recombination in HIV appears not so unpredictable, therefore, and the existence of recombinants that frequently generate nonfunctional products highlights previously unappreciated limits of the genetic flexibility of HIV.

Structural and functional modularity of proteins in the de novo purine biosynthetic pathway

It is generally accepted that naturally existing functional domains can serve as building blocks for complex protein structures, and that novel functions can arise from assembly of different combinations of these functional domains. To inform our understanding of protein evolution and explore the modular nature of protein structure, two model enzymes were chosen for study, purT-encoded glycinamide ribonucleotide formyltransferase (PurT) and purK-encoded N(5)-carboxylaminoimidazole ribonucleotide synthetase (PurK). Both enzymes are found in the de novo purine biosynthetic pathway of Escherichia coli. In spite of their low sequence identity, PurT and PurK share significant similarity in terms of tertiary structure, active site organization, and reaction mechanism. Their characteristic three domain structures categorize both PurT and PurK as members of the ATP-grasp protein superfamily. In this study, we investigate the exchangeability of individual protein domains between these two enzymes and the in vivo and in vitro functional properties of the resulting hybrids. Six domain-swapped hybrids were unable to catalyze full wild-type reactions, but each hybrid protein could catalyze partial reactions. Notably, an additional loop replacement in one of the domain-swapped hybrid proteins was able to restore near wild-type PurK activity. Therefore, in this model system, domain-swapped proteins retained the ability to catalyze partial reactions, but further modifications were required to efficiently couple the reaction intermediates and achieve catalysis of the full reaction. Implications for understanding the role of domain swapping in protein evolution are discussed.

Ligand binding and allostery can emerge simultaneously

A heterotropic allosteric effect involves an effector molecule that is distinct from the substrate or ligand of the protein. How heterotropic allostery originates is an unanswered question. We have previously created several heterotropic allosteric enzymes by recombining the genes for TEM1 beta-lactamase (BLA) and maltose binding protein (MBP) to create BLAs that are positively or negatively regulated by maltose. We show here that one of these engineered enzymes has approximately 10(6) M(-1) affinity for Zn(2+), a property that neither of the parental proteins possesses. Furthermore, Zn(2+) is a negative effector that noncompetitively switches off beta-lactam hydrolysis activity. Mutagenesis experiments indicate that the Zn(2+)-binding site does not involve a histidine or a cysteine, which is atypical of natural Zn(2+)-binding sites. These studies also implicate helices 1 and 12 of the BLA domain in allosteric signal propagation. These results support a model for the evolution of heterotropic allostery in which effector affinity and allosteric signaling emerge simultaneously.

Directed Evolution and Identification of Control Regions of ColE1 Plasmid Replication Origins Using Only Nucleotide Deletions

Genes can be mutated by altering DNA content (base changes) or DNA length (insertions or deletions). Most in vitro directed evolution processes utilize nucleotide content changes to produce DNA libraries. We tested whether gain of function mutations could be identified using a mutagenic process that produced only nucleotide deletions. Short nucleotide stretches were deleted in a plasmid encoding lacZ, and screened for increased beta-galactosidase activity. Several mutations were found in the origin of replication that quantitatively and qualitatively altered plasmid behavior in vivo. Some mutations allowed co-residence of ColE1 plasmids in Escherichia coli, and implicate hairpin structures II and III of the ColE1 RNA primer as determinants of plasmid compatibility. Thus, useful and unexpected mutations can be found from libraries containing only deletions.

A Molecular Switch Created by In Vitro Recombination of Nonhomologous Genes

We have created a molecular switch by the in vitro recombination of nonhomologous genes and subjecting the recombined genes to evolutionary pressure. The gene encoding TEM1 beta-lactamase was circularly permuted in a random fashion and subsequently randomly inserted into the gene encoding Escherichia coli maltose binding protein. From this library, a switch (RG13) was identified in which its beta-lactam hydrolysis activity was compromised in the absence of maltose but increased 25-fold in the presence of maltose. Upon removal of maltose, RG13's catalytic activity returned to its premaltose level, illustrating that the switching is reversible. The modularity of RG13 was demonstrated by increasing maltose affinity while preserving switching activity. RG13 gave rise to a novel cellular phenotype, illustrating the potential of molecular switches to rewire the cellular circuitry.

Probing the S100 protein family through genomic and functional analysis

The EF-hand superfamily of calcium binding proteins includes the S100, calcium binding protein, and troponin subfamilies. This study represents a genome, structure, and expression analysis of the S100 protein family, in mouse, human, and rat. We confirm the high level of conservation between mammalian sequences but show that four members, including S100A12, are present only in the human genome. We describe three new members of the S100 family in the three species and their locations within the S100 genomic clusters and propose a revised nomenclature and phylogenetic relationship between members of the EF-hand superfamily. Two of the three new genes were induced in bone-marrow-derived macrophages activated with bacterial lipopolysaccharide, suggesting a role in inflammation. Normal human and murine tissue distribution profiles indicate that some members of the family are expressed in a specific manner, whereas others are more ubiquitous. Structure-function analysis of the chemotactic properties of murine S100A8 and human S100A12, particularly within the active hinge domain, suggests that the human protein is the functional homolog of the murine protein. Strong similarities between the promoter regions of human S100A12 and murine S100A8 support this possibility. This study provides insights into the possible processes of evolution of the EF-hand protein superfamily. Evolution of the S100 proteins appears to have occurred in a modular fashion, also seen in other protein families such as the C2H2-type zinc-finger family.

Structural evidence for a hydride transfer mechanism of catalysis in phosphoglucose isomerase from Pyrococcus furiosus

In the Euryarchaeota species Pyrococcus furiosus and Thermococcus litoralis, phosphoglucose isomerase (PGI) activity is catalyzed by an enzyme unrelated to the well known family of PGI enzymes found in prokaryotes, eukaryotes, and some archaea. We have determined the crystal structure of PGI from Pyrococcus furiosus in native form and in complex with two active site ligands, 5-phosphoarabinonate and gluconate 6-phosphate. In these structures, the metal ion, which in vivo is presumed to be Fe2+, is located in the core of the cupin fold and is immediately adjacent to the C1-C2 region of the ligands, suggesting that Fe2+ is involved in catalysis rather than serving a structural role. The active site contains a glutamate residue that contacts the substrate, but, because it is also coordinated to the metal ion, it is highly unlikely to mediate proton transfer in a cis-enediol mechanism. Consequently, we propose a hydride shift mechanism of catalysis. In this mechanism, Fe2+ is responsible for proton transfer between O1 and O2, and the hydride shift between C1 and C2 is favored by a markedly hydrophobic environment in the active site. The absence of any obvious enzymatic machinery for catalyzing ring opening of the sugar substrates suggests that pyrococcal PGI has a preference for straight chain substrates and that metabolism in extreme thermophiles may use sugars in both ring and straight chain forms.

Library analysis of SCHEMA-guided protein recombination

The computational algorithm SCHEMA was developed to estimate the disruption caused when amino acid residues that interact in the three-dimensional structure of a protein are inherited from different parents upon recombination. To evaluate how well SCHEMA predicts disruption, we have shuffled the distantly-related beta-lactamases PSE-4 and TEM-1 at 13 sites to create a library of 2(14) (16,384) chimeras and examined which ones retain lactamase function. Sequencing the genes from ampicillin-selected clones revealed that the percentage of functional clones decreased exponentially with increasing calculated disruption (E = the number of residue-residue contacts that are broken upon recombination). We also found that chimeras with low E have a higher probability of maintaining lactamase function than chimeras with the same effective level of mutation but chosen at random from the library. Thus, the simple distance metric used by SCHEMA to identify interactions and compute E allows one to predict which chimera sequences are most likely to retain their function. This approach can be used to evaluate crossover sites for recombination and to create highly mosaic, folded chimeras.

Exploring the structure and function of the mycobacterial KatG protein using trans-dominant mutants

In order to probe the structure and function of the mycobacterial catalase-peroxidase enzyme (KatG), we employed a genetic approach using dominant-negative analysis of katG merodiploids. Transformation of Mycobacterium bovis BCG with various katG point mutants (expressed from low-copy-number plasmids) resulted in reductions in peroxidase and catalase activities as measured in cell extracts. These reductions in enzymatic activity usually correlated with increased resistance to the antituberculosis drug isoniazid (INH). However, for the N138S trans-dominant mutant, the catalase-peroxidase activity was significantly decreased while the sensitivity to INH was retained. trans-dominance required katG expression from multicopy plasmids and could not be demonstrated with katG mutants integrated elsewhere on the wild-type M. bovis BCG chromosome. Reversal of the mutant phenotype through plasmid exchange suggested the catalase-peroxidase deficiency occurred at the protein level and that INH resistance was not due to a second site mutation(s). Electrophoretic analysis of KatG proteins from the trans-dominant mutants showed a reduction in KatG dimers compared to WT and formation of heterodimers with reduced activity. The mutants responsible for these defects cluster around proposed active site residues: N138S, T275P, S315T, and D381G. In an attempt to identify mutants that might delimit the region(s) of KatG involved in subunit interactions, C-terminal truncations were constructed (with and without the D381G dominant-negative mutation). None of the C-terminal deletions were able to complement a DeltakatG strain, nor could they cause a dominant-negative effect on the WT. Taken together, these results suggest an intricate association between the amino- and carboxy-terminal regions of KatG and may be consistent with a domain-swapping mechanism for KatG dimer formation.

An analysis of protein domain linkers: their classification and role in protein folding

Recent advances in protein engineering have come from creating multi-functional chimeric proteins containing modules from various proteins. These modules are typically joined via an oligopeptide linker, the correct design of which is crucial for the desired function of the chimeric protein. Here we analyse the properties of naturally occurring inter-domain linkers with the aim to design linkers for domain fusion. Two main types of linker were identified; helical and non-helical. Helical linkers are thought to act as rigid spacers separating two domains. Non-helical linkers are rich in prolines, which also leads to structural rigidity and isolation of the linker from the attached domains. This means that both linker types are likely to act as a scaffold to prevent unfavourable interactions between folding domains. Based on these results we have constructed a linker database intended for the rational design of linkers for domain fusion, which can be accessed via the Internet at http://mathbio.nimr.mrc.ac.uk.

Control of domain swapping in bovine odorant-binding protein

As revealed by the X-ray structure, bovine odorant-binding protein (OBPb) is a domain swapped dimer [Tegoni, Ramoni, Bignetti, Spinelli and Cambillau (1996) Nat. Struct. Biol. 3, 863-867; Bianchet, Bains, Petosi, Pevsner, Snyder, Monaco and Amzel (1996) Nat. Struct. Biol. 3, 934-939]. This contrasts with all known mammalian OBPs, which are monomers, and in particular with porcine OBP (OBPp), sharing 42.3% identity with OBPb. By the mechanism of domain swapping, monomers are proposed to evolve into dimers and oligomers, as observed in human prion. Comparison of bovine and porcine OBP sequences pointed at OBPp glycine 121, in the hinge linking the beta-barrel to the alpha-helix. The absence of this residue in OBPb might explain why the normal lipocalin beta-turn is not formed. In order to decipher the domain swapping determinants we have produced a mutant of OBPb in which a glycine residue was inserted after position 121, and a mutant of OBPp in which glycine 121 was deleted. The latter mutation did not result in dimerization, while OBPb-121Gly+ became monomeric, suggesting that domain swapping was reversed. Careful structural analysis revealed that besides the presence of a glycine in the hinge, the dimer interface formed by the C-termini and by the presence of the lipocalins conserved disulphide bridge may also control domain swapping.

Directed enzyme evolution

Laboratory evolutionists continue to generate better enzymes for industrial and research applications. Exciting developments include new biocatalysts for enantioselective carbon-carbon bond formation and fatty acid production in plants. Creative contributions to the repertoire of evolutionary methods will ensure further growth in applications and expand the scope and complexity of biological design problems that can be addressed. Researchers are also starting to elucidate mechanisms of enzyme adaptation and natural evolution by testing evolutionary scenarios in the laboratory.

Libraries of hybrid proteins from distantly related sequences

We introduce a method for sequence homology-independent protein recombination (SHIPREC) that can create libraries of single-crossover hybrids of unrelated or distantly related proteins. The method maintains the proper sequence alignment between the parents and introduces crossovers mainly at structurally related sites distributed over the aligned sequences. We used SHIPREC to create a library of interspecies hybrids of a membrane-associated human cytochrome P450 (1A2) and the heme domain of a soluble bacterial P450 (BM3). By fusing the hybrid gene library to the gene for chloramphenicol acetyl transferase (CAT), we were able to select for soluble and properly folded protein variants. Screening for 1A2 activity (deethylation of 7-ethoxyresorufin) identified two functional P450 hybrids that were more soluble in the bacterial cytoplasm than the wild-type 1A2 enzyme.

Combinatorial and computational challenges for biocatalyst design

Nature provides a fantastic array of catalysts extremely well suited to supporting life, but usually not so well suited for technology. Whether biocatalysis will have a significant technological impact depends on our finding robust routes for tailoring nature's catalysts or redesigning them anew. Laboratory evolution methods are now used widely to fine-tune the selectivity and activity of enzymes. The current rapid development of these combinatorial methods promises solutions to more complex problems, including the creation of new biosynthetic pathways. Computational methods are also developing quickly. The marriage of these approaches will allow us to generate the efficient, effective catalysts needed by the pharmaceutical, food and chemicals industries and should open up new opportunities for producing energy and chemicals from renewable resources.