Circular permutation in proteins

Protein engineering using circular permutation - structure, function, stability, and applications

Protein engineering is important for creating novel variants from natural proteins, enabling a wide range of applications. Approaches such as rational design and directed evolution are routinely used to make new protein variants. Computational tools like de novo design can introduce new protein folds. Expanding the amino acid repertoire to include unnatural amino acids with non-canonical side chains in vitro by native chemical ligation and in vivo via codon expansion methods broadens sequence and structural possibilities. Circular permutation (CP) is an invaluable approach to redesigning a protein by rearranging the amino acid sequence, where the connectivity of the secondary structural elements is altered without changing the overall structure of the protein. Artificial CP proteins (CPs) are employed in various applications such as biocatalysis, sensing of small molecules by fluorescence, genome editing, ligand-binding protein switches, and optogenetic engineering. Many studies have shown that CP can lead to either reduced or enhanced stability or catalytic efficiency. The effects of CP on a protein's energy landscape cannot be predicted a priori. Thus, it is important to understand how CP can affect the thermodynamic and kinetic stability of a protein. In this review, we discuss the discovery and advancement of techniques to create protein CP, and existing reviews on CP. We delve into the plethora of biological applications for designed CP proteins. We subsequently discuss the experimental and computational reports on the effects of CP on the thermodynamic and kinetic stabilities of proteins of various topologies. An understanding of the various aspects of CP will allow the reader to design robust CP proteins for their specific purposes.

State-of-the-art in engineering small molecule biosensors and their applications in metabolic engineering

Genetically encoded biosensors are crucial for enhancing our understanding of how molecules regulate biological systems. Small molecule biosensors, in particular, help us understand the interaction between chemicals and biological processes. They also accelerate metabolic engineering by increasing screening throughput and eliminating the need for sample preparation through traditional chemical analysis. Additionally, they offer significantly higher spatial and temporal resolution in cellular analyte measurements. In this review, we discuss recent progress in in vivo biosensors and control systems-biosensor-based controllers-for metabolic engineering. We also specifically explore protein-based biosensors that utilize less commonly exploited signaling mechanisms, such as protein stability and induced degradation, compared to more prevalent transcription factor and allosteric regulation mechanism. We propose that these lesser-used mechanisms will be significant for engineering eukaryotic systems and slower-growing prokaryotic systems where protein turnover may facilitate more rapid and reliable measurement and regulation of the current cellular state. Lastly, we emphasize the utilization of cutting-edge and state-of-the-art techniques in the development of protein-based biosensors, achieved through rational design, directed evolution, and collaborative approaches.

The Biochemical Impact of Extracting an Embedded Adenylate Kinase Domain Using Circular Permutation

Adenylate kinases (AKs) have evolved AMP-binding and lid domains that are encoded as continuous polypeptides embedded at different locations within the discontinuous polypeptide encoding the core domain. A prior study showed that AK homologues of different stabilities consistently retain cellular activity following circular permutation that splits a region with high energetic frustration within the AMP-binding domain into discontinuous fragments. Herein, we show that mesophilic and thermophilic AKs having this topological restructuring retain activity and substrate-binding characteristics of the parental AK. While permutation decreased the activity of both AK homologues at physiological temperatures, the catalytic activity of the thermophilic AK increased upon permutation when assayed >30 °C below the melting temperature of the native AK. The thermostabilities of the permuted AKs were uniformly lower than those of native AKs, and they exhibited multiphasic unfolding transitions, unlike the native AKs, which presented cooperative thermal unfolding. In addition, proteolytic digestion revealed that permutation destabilized each AK in differing manners, and mass spectrometry suggested that the new termini within the AMP-binding domain were responsible for the increased proteolysis sensitivity. These findings illustrate how changes in contact order can be used to tune enzyme activity and alter folding dynamics in multidomain enzymes.

Optimization strategies and advances in the research and development of AAV-based gene therapy to deliver large transgenes

Adeno-associated virus (AAV)-based therapies are recognized as one of the most potent next-generation treatments for inherited and genetic diseases. However, several biological and technological aspects of AAV vectors remain a critical issue for their widespread clinical application. Among them, the limited capacity of the AAV genome significantly hinders the development of AAV-based gene therapy. In this context, genetically modified transgenes compatible with AAV are opening up new opportunities for unlimited gene therapies for many genetic disorders. Recent advances in de novo protein design and remodelling are paving the way for new, more efficient and targeted gene therapeutics. Using computational and genetic tools, AAV expression cassette and transgenic DNA can be split, miniaturized, shuffled or created from scratch to mediate efficient gene transfer into targeted cells. In this review, we highlight recent advances in AAV-based gene therapy with a focus on its use in translational research. We summarize recent research and development in gene therapy, with an emphasis on large transgenes (>4.8 kb) and optimizing strategies applied by biomedical companies in the research pipeline. We critically discuss the prospects for AAV-based treatment and some emerging challenges. We anticipate that the continued development of novel computational tools will lead to rapid advances in basic gene therapy research and translational studies.

Circular permutation at azurin's active site slows down its folding

Circular permutation (CP) is a technique by which the primary sequence of a protein is rearranged to create new termini. The connectivity of the protein is altered but the overall protein structure generally remains unperturbed. Understanding the effect of CP can help design robust proteins for numerous applications such as in genetic engineering, optoelectronics, and improving catalytic activity. Studies on different protein topologies showed that CP usually affects protein stability as well as unfolding rates. Though a significant number of proteins contain metals or other cofactors, reports of metalloprotein CPs are rare. Thus, we chose a bacterial metalloprotein, azurin, and its CP within the metal-binding site (cpF114). We studied the stabilities, folding, and unfolding rates of apo- and Zn2+-bound CP azurin using fluorescence and circular dichroism. The introduced CP had destabilizing effects on the protein. Also, the folding of the Zn2+-CP protein was much slower than that of the Zn2+-WT or apo-protein. We compared this study to our previously reported azurin-cpN42, where we had observed an equilibrium and kinetic intermediate. cpF114 exhibits an apparent two-state equilibrium unfolding but has an off-pathway kinetic intermediate. Our study hinted at CP as a method to modify the energy landscape of proteins to alter their folding pathways. WT azurin, being a faster folder, may have evolved to optimize the folding rate of metal-bound protein compared to its CPs, albeit all of them have the same structure and function. Our study underscores that protein sequence and protein termini positions are crucial for metalloproteins. TOC Figure. (Top) Zn2+-azurin WT structure (PDB code: 1E67) and 2-D topology diagram of Zn2+-cpF114 azurin. (Bottom) Cartoon diagram representing folding (red arrows) and unfolding (blue arrows) of apo- and Zn2+- WT and cpF114 azurins. The width of the arrows represents the rate of the corresponding processes.

A single-domain protein catenane of dihydrofolate reductase

A single-domain protein catenane refers to two mechanically interlocked polypeptide rings that fold synergistically into a compact and integrated structure, which is extremely rare in nature. Here, we report a single-domain protein catenane of dihydrofolate reductase (cat-DHFR). This design was achieved by rewiring the connectivity between secondary motifs to introduce artificial entanglement and synthesis was readily accomplished through a series of programmed and streamlined post-translational processing events in cells without any additional in vitro reactions. The target molecule contained few exogenous motifs and was thoroughly characterized using a combination of ultra-performance liquid chromatography-mass spectrometry, sodium dodecyl sulfate-polyacrylamide gel electrophoresis, protease cleavage experiments and ion mobility spectrometry-mass spectrometry. Compared with the linear control, cat-DHFR retained its catalytic capability and exhibited enhanced stability against thermal or chemical denaturation due to conformational restriction. These results suggest that linear proteins may be converted into their concatenated single-domain counterparts with almost identical chemical compositions, well-preserved functions and elevated stabilities, representing an entirely new horizon in protein science.

Elucidation of the folding pathway of a circular permutant of topologically knotted YbeA by tryptophan substitutions

CP74 is an engineered circular permutant of a deep trefoil knotted SpoU-TrmD (SPOUT) RNA methyl transferase protein YbeA from E. coli. We have previously established that the circular permutation unties the knotted topology of YbeA and CP74 forms a domain-swapped dimer with a large dimeric interface of ca. 4600 Å². To understand the impact of domain-swapping and the newly formed hinge region joining the two folded domains on the folding and stability of CP74, the five equally spaced tryptophan residues were individually substituted into phenylalanine to monitor their conformational and stability changes by a battery of biophysical tools. Far-UV circular dichroism, intrinsic fluorescence, and small-angle X-ray scattering dictated minimal global conformational perturbations to the native structures in the tryptophan variants. The structures of the tryptophan variants also showed the conservation of the domain-swapped ternary structure with the exception that the W72F exhibited significant asymmetry in the α-helix 5. Comparative global thermal and chemical stability analyses indicated the pivotal role of W100 in the folding of CP74 followed by W19 and W72. Solution-state NMR spectroscopy and hydrogen-deuterium exchange mass spectrometry further revealed the accumulation of a native-like intermediate state in which the hinge region made important contributions to maintain the domain-swapped ternary structure of CP74.

Ribonuclease T2 represents a distinct circularly permutated version of the BECR RNases

Detection of homologous relationships among proteins and understanding their mechanisms of diversification are major topics in the fields of protein science, bioinformatics, and phylogenetics. Recent developments in sequence/profile-based and structural similarity-based methods have greatly facilitated the unification and classification of many protein families into superfamilies or folds, yet many proteins remain unclassified in current protein databases. As one of the three earliest identified RNases in biology, ribonuclease T2, also known as RNase I in Escherichia coli, RNase Rh in fungi, or S-RNase in plant, is thought to be an ancient RNase family due to its widespread distribution and distinct structure. In this study, we present evidence that RNase T2 represents a circularly permutated version of the BECR (Barnase-EndoU-Colicin E5/D-RelE) fold RNases. This subtle relationship cannot be detected by traditional methods such as sequence/profile-based comparisons, structure-similarity searches, and circular permutation detections. However, we were able to identify the structural similarity using rational reconstruction of a theoretical RNase T2 ancestor via a reverse circular permutation process, followed by structural modeling using AlphaFold2, and structural comparisons. This relationship is further supported by the fact that RNase T2 and other typical BECR RNases, namely Colicin D, RNase A, and BrnT, share similar catalytic site configurations, all involving an analogous set of conserved residues on the α0 helix and the β4 strand of the BECR fold. This study revealed a hidden root of RNase T2 in bacterial toxin systems and demonstrated that reconstruction and modeling of ancestral topology is an effective strategy to identify remote relationship between proteins.

SeqCP: A sequence-based algorithm for searching circularly permuted proteins

Circular permutation (CP) is a protein sequence rearrangement in which the amino- and carboxyl-termini of a protein can be created in different positions along the imaginary circularized sequence. Circularly permutated proteins usually exhibit conserved three-dimensional structures and functions. By comparing the structures of circular permutants (CPMs), protein research and bioengineering applications can be approached in ways that are difficult to achieve by traditional mutagenesis. Most current CP detection algorithms depend on structural information. Because there is a vast number of proteins with unknown structures, many CP pairs may remain unidentified. An efficient sequence-based CP detector will help identify more CP pairs and advance many protein studies. For instance, some hypothetical proteins may have CPMs with known functions and structures that are informative for functional annotation, but existing structure-based CP search methods cannot be applied when those hypothetical proteins lack structural information. Despite the considerable potential for applications, sequence-based CP search methods have not been well developed. We present a sequence-based method, SeqCP, which analyzes normal and duplicated sequence alignments to identify CPMs and determine candidate CP sites for proteins. SeqCP was trained by data obtained from the Circular Permutation Database and tested with nonredundant datasets from the Protein Data Bank. It shows high reliability in CP identification and achieves an AUC of 0.9. SeqCP has been implemented into a web server available at: http://pcnas.life.nthu.edu.tw/SeqCP/.

Comparative phylogeny and evolutionary analysis of Dicer-like protein family in two plant monophyletic lineages

BACKGROUND

Small RNAs (sRNAs) that do not get untranslated into proteins exhibit a pivotal role in the expression regulation of their cognate gene(s) in almost all eukaryotic lineages, including plants. Hitherto, numerous protein families such as Dicer, a unique class of Ribonuclease III, have been reported to be involved in sRNAs processing pathways and silencing. In this study, we aimed to investigate the phylogenetic relationship and evolutionary history of the DCL protein family.

RESULTS

Our results illustrated the DCL family of proteins grouped into four main subfamilies (DCLs 1-4) presented in either Eudicotyledons or Liliopsids. The accurate observation of the phylogenetic trees supports the independent expansion of DCL proteins among the Eudicotyledons and Liliopsids species. They share the common origin, and the main duplication events for the formation of the DCL subfamilies occurred before the Eudicotyledons/Liliopsids split from their ancestral DCL. In addition, shreds of evidence revealed that the divergence happened when multicellularization started and since the need for complex gene regulation considered being a necessity by organisms. At that time, they have evolved independently among the monophyletic lineages. The other finding was that the combination of DCL protein subfamilies bears several highly conserved functional domains in plant species that originated from their ancestor architecture. The conservation of these domains happens to be both lineage-specific and inter lineage-specific.

CONCLUSIONS

DCL subfamilies (i.e., DCL1-DCL4) distribute in their single clades after diverging from their common ancestor and before emerging into higher plants. Therefore, it seems that the main duplication events for the formation of the DCL subfamilies occurred before the Eudicotyledons/Liliopsida split and before the appearance of moss, and after the single-cell green algae. We also observed the same trends among the main DCL subfamilies from functional unit composition and architecture. Despite the long evolutionary course from the divergence of Liliopsida lineage from the Eudicotyledons, a significant diversifying force to domain composition and orientation was absent. The results of this study provide a deeper insight into DCL protein evolutionary history and possible sequence and structural relationships between DCL protein subfamilies in the main higher plant monophyletic lineages; i.e., Eudicotyledons and Liliopsida.

Obtaining protein foldability information from computational models of AlphaFold2 and RoseTTAFold

The recent breakthrough from AlphaFold2 and RoseTTAFold set a profound milestone for solving the protein folding problem, but they were not explicitly trained to predict protein foldability, i.e., if a protein can really fold into the predicted 3D structure. We wondered if the computational models from AlphaFold2 and RoseTTAFold might carry protein foldability information. Therefore, we predicted the structural models of 159 circular permutants and 158 alanine insertion mutants of the 159-residue dihydrofolate reductase. Our data showed that although AlphaFold2 and RoseTTAFold cannot directly identify unfoldable proteins, the RMSD values of computational models are correlated with protein foldability, with higher RMSD values indicating lower protein foldability. Furthermore, this correlation is independent of secondary structures, and the RMSD values of computational models are quantitatively correlated with protein foldability but not protein functions. Additionally, using a dataset of 129 de novo designed proteins, we showed that inter-model RMSD values between AlphaFold2 models and RoseTTAFold models are a good indicator of protein foldability. At last, we showed that inter-model RMSD values are also useful for evaluating protein solubility by modeling 1664 natural proteins. Our work could be of great value to the design of novel proteins and the prediction of protein foldability.

Application of Restriction Free (RF) Cloning in Circular Permutation

The restriction free (RF) cloning has emerged as one of the highly efficient techniques in the area of genetic engineering. RF cloning has wide range of applications in plasmid DNA manipulation including cloning of a single gene, simultaneous assembly of multiple DNA fragments, and mutagenesis from single to multiple simultaneous alterations of a target DNA. Recently, we have developed a new technique of circular permutation using RF cloning. Circular permutation is widely used to investigate the mechanisms of protein folding and function. Previously, restriction enzyme based cloning was used to introduce circular permutation. Our RF cloning method made the protocol faster and more cost-effective. In this chapter, we describe a step-by-step protocol for generating circular permutants using RF methodology.

Cryogenian Origin and Subsequent Diversification of the Plant Cell-Wall Enzyme XTH Family

All land plants encode large multigene families of xyloglucan endotransglucosylase/hydrolases (XTHs), plant-specific enzymes that cleave and reconnect plant cell-wall polysaccharides. Despite the ubiquity of these enzymes, considerable uncertainty remains regarding the evolutionary history of the XTH family. Phylogenomic and comparative analyses in this study traced the non-plant origins of the XTH family to Alphaproteobacteria ExoKs, bacterial enzymes involved in loosening biofilms, rather than Firmicutes licheninases, plant biomass digesting enzymes, as previously supposed. The relevant horizontal gene transfer (HGT) event was mapped to the divergence of non-swimming charophycean algae in the Cryogenian geological period. This HGT event was the likely origin of charophycean EG16-2s, which are putative intermediates between ExoKs and XTHs. Another HGT event in the Cryogenian may have led from EG16-2s or ExoKs to fungal Congo Red Hypersensitive proteins (CRHs) to fungal CRHs, enzymes that cleave and reconnect chitin and glucans in fungal cell walls. This successive transfer of enzyme-encoding genes may have supported the adaptation of plants and fungi to the ancient icy environment by facilitating their sessile lifestyles. Furthermore, several protein evolutionary steps, including coevolution of substrate-interacting residues and putative intra-family gene fusion, occurred in the land plant lineage and drove diversification of the XTH family. At least some of those events correlated with the evolutionary gain of broader substrate specificities, which may have underpinned the expansion of the XTH family by enhancing duplicated gene survival. Together, this study highlights the Precambrian evolution of life and the mode of multigene family expansion in the evolutionary history of the XTH family.

Purification and characterization of a highly thermostable GlcNAc-binding lectin from Collaea speciosa seeds

Lectins from plants of the Diocleinae subtribe often exhibit specificity towards mannose/glucose and derived sugars, with some plants also displaying a second lectin specific to lactose/GalNAc. Here, we present a novel lectin from Collaea speciosa, named CsL, that displays specificity for GlcNAc/glucose. The lectin was extracted from Collaea speciosa seeds and purified by a single chromatographic step on a Sephadex G-50 matrix. In solution, the lectin appears as a dimeric protein composed of 25 kDa monomers. The protein is stable at pH 7-8 and dependent on divalent cations. CsL maintained its agglutination activity after heating to 90 °C for 1 h. Glycan array studies revealed that CsL binds to N-glycans with terminal GlcNAc residues, chitobiose and chitotriose moieties. The partial amino acid sequence of the lectin is similar to that of some lactose-specific lectins from the same subtribe. In contrast to other ConA-like lectins, CsL is not toxic to Artemia. Because of its remarkably different properties and specificity, this lectin could be the first member of a new group inside the Diocleinae lectins.

Structural and biochemical analyses of concanavalin A circular permutation by jack bean asparaginyl endopeptidase

Over 30 years ago, an intriguing posttranslational modification was found responsible for creating concanavalin A (conA), a carbohydrate-binding protein from jack bean (Canavalia ensiformis) seeds and a common carbohydrate chromatography reagent. ConA biosynthesis involves what was then an unprecedented rearrangement in amino-acid sequence, whereby the N-terminal half of the gene-encoded conA precursor (pro-conA) is swapped to become the C-terminal half of conA. Asparaginyl endopeptidase (AEP) was shown to be involved, but its mechanism was not fully elucidated. To understand the structural basis and consequences of circular permutation, we generated recombinant jack bean pro-conA plus jack bean AEP (CeAEP1) and solved crystal structures for each to 2.1 and 2.7 Å, respectively. By reconstituting conA biosynthesis in vitro, we prove CeAEP1 alone can perform both cleavage and cleavage-coupled transpeptidation to form conA. CeAEP1 structural analysis reveals how it is capable of carrying out both reactions. Biophysical assays illustrated that pro-conA is less stable than conA. This observation was explained by fewer intermolecular interactions between subunits in the pro-conA crystal structure and consistent with a difference in the prevalence for tetramerization in solution. These findings elucidate the consequences of circular permutation in the only posttranslation example known to occur in nature.

Roles of Plant-Specific Inserts in Plant Defense

Ubiquitously expressed in plants, the plant-specific insert (PSI) of typical plant aspartic proteases (tpAPs) has been associated with plant development, stress response, and defense processes against invading pathogens. Despite sharing high sequence identity, structural studies revealed possible different mechanisms of action among species. The PSI induces signaling pathways of defense hormones in vivo and demonstrates broad-spectrum activity against phytopathogens in vitro. Recent characterization of the PSI-tpAP relationship uncovered novel, nonconventional intracellular protein transport pathways and improved tpAP production yields for industrial applications. In spite of research to date, relatively little is known about the structure-function relationships of PSIs. A comprehensive understanding of their biological roles may benefit plant protection strategies against virulent phytopathogens.

Versatile format of minichaperone-based protein fusion system

Hydrophobic recombinant proteins often tend to aggregate upon expression into inclusion bodies and are difficult to refold. Producing them in soluble forms constitutes a common bottleneck problem. A fusion system for production of insoluble hydrophobic proteins in soluble stable forms with thermophilic minichaperone, GroEL apical domain (GrAD) as a carrier, has recently been developed. To provide the utmost flexibility of the system for interactions between the carrier and various target protein moieties a strategy of making permutated protein variants by gene engineering has been applied: the original N- and C-termini of the minichaperone were linked together by a polypeptide linker and new N- and C-termini were made at desired parts of the protein surface. Two permutated GrAD forms were created and analyzed. Constructs of GrAD and both of its permutated forms fused with the initially insoluble N-terminal fragment of hepatitis C virus' E2 protein were tested. Expressed fusions formed inclusion bodies. After denaturation, all fusions were completely renatured in stable soluble forms. A variety of permutated GrAD variants can be created. The versatile format of the system provides opportunities for choosing an optimal pair between particular target protein moiety and the best-suited original or specific permutated carrier.

Thermostable Lichenase from Clostridium thermocellum as a Host Protein in the Domain Insertion Approach

Clostridium thermocellum lichenase (endo-β-1,3;1,4-glucan-D-glycosyl hydrolase, EC 3.2.1.73 (P29716)) has been tested for the insertion of two model fluorescent proteins (EGFP and TagRFP) into two regions of this enzyme. Functional folding of the resulting proteins was confirmed by retention of lichenase activity and EGFP and TagRFP fluorescence. These results convincingly demonstrate that (i) the two experimentally selected lichenase loop regions may serve as the areas for domain insertion without disturbing enzyme folding in vivo; (ii) lichenase permits not only single but also tandem insertions of large protein domains. High specific activity, outstanding thermostability, and efficient in vitro refolding of thermostable lichenase make it an attractive new host protein for the insertional fusion of domains in the engineering of multifunctional proteins.

Tandem domain swapping: determinants of multidomain protein misfolding

•

Domain swapping refers to the exchange of structural elements between protein domains.

•

Experiments show that tandem homologous domains are prone to domain swapping.

•

Recent studies establish a framework to understand the formation of tandem domain swaps.

•

Prediction of tandem domain swaps is possible but hindered by the amount of available data.

Tandem homologous domains in proteins are susceptible to misfolding through the formation of domain swaps, non-native conformations involving the exchange of equivalent structural elements between adjacent domains. Cutting-edge biophysical experiments have recently allowed the observation of tandem domain swapping events at the single molecule level. In addition, computer simulations have shed light into the molecular mechanisms of domain swap formation and serve as the basis for methods to systematically predict them. At present, the number of studies on tandem domain swaps is still small and limited to a few domain folds, but they offer important insights into the folding and evolution of multidomain proteins with applications in the field of protein design.

Facilitating circular permutation using Restriction Free (RF) cloning

Circular permutation is a powerful tool to test the role of topology in protein folding and function. Previous methods for generating circular permutants were based on rearranging gene elements using restriction enzymes-based cloning. Here, we present a Restriction Free (RF) approach to achieve circular permutation which is faster and more cost-effective.

Circular Permutation Obscures Universality of a Ribosomal Protein

Functions, origins, and evolution of the translation system are best understood in the context of unambiguous and phylogenetically based taxonomy and nomenclature. Here, we map ribosomal proteins onto the tree of life and provide a nomenclature for ribosomal proteins that is consistent with phylogenetic relationships. We have increased the accuracy of homology relationships among ribosomal proteins, providing a more informative picture of their lineages. We demonstrate that bL33 (bacteria) and eL42 (archaea/eukarya) are homologs with common ancestry and acute similarities in sequence and structure. Their similarities were previously obscured by circular permutation. The most likely mechanism of permutation between bL33 and eL42 is duplication followed by fusion and deletion of both the first and last β-hairpins. bL33 and eL42 are composed of zinc ribbon protein folds, one of the most common zinc finger fold-groups of, and most frequently observed in translation-related domains. Bacterial-specific ribosomal protein bL33 and archaeal/eukaryotic-specific ribosomal protein eL42 are now both assigned the name of uL33, indicating a universal ribosomal protein. We provide a phylogenetic naming scheme for all ribosomal proteins that is based on phylogenetic relationships to be used as a tool for studying the systemics, evolution, and origins of the ribosome.

Protein topology determines substrate-binding mechanism in homologous enzymes

During evolution, some homologs proteins appear with different connectivity between secondary structures (different topology) but conserving the tridimensional arrangement of them (same architecture). These events can produce two types of arrangements; circular permutation or non-cyclic permutations. The first one results in the N and C terminus transferring to a different position on a protein sequence while the second refers to a more complex arrangement of the structural elements. In ribokinase superfamily, two different topologies can be identified, which are related to each other as a non-cyclic permutation occurred during the evolution. Interestingly, this change in topology is correlated with the nucleotide specificity of its members. Thereby, the connectivity of the secondary elements allows us to distinguish an ATP-dependent and an ADP-dependent topology. Here we address the impact of introducing the topology of a homologous ATP-dependent kinase in an ADP-dependent kinase (Thermococcus litoralis glucokinase) in the structure, nucleotide specificity, and substrate binding order of the engineered enzyme. Structural evidence demonstrates that rewiring the topology of TlGK leads to an active and soluble enzyme without modifications on its three-dimensional architecture. The permuted enzyme (PerGK) retains the nucleotide preference of the parent TlGK enzyme but shows a change in the substrate binding order. Our results illustrate how the rearrangement of the protein folding topology during the evolution of the ribokinase superfamily enzymes may have dictated the substrate-binding order in homologous enzymes of this superfamily.

Circularly permuted variants of two CG-specific prokaryotic DNA methyltransferases

The highly similar prokaryotic DNA (cytosine-5) methyltransferases (C5-MTases) M.MpeI and M.SssI share the specificity of eukaryotic C5-MTases (5'-CG), and can be useful research tools in the study of eukaryotic DNA methylation and epigenetic regulation. In an effort to improve the stability and solubility of complementing fragments of the two MTases, genes encoding circularly permuted (CP) variants of M.MpeI and M.SssI were created, and cloned in a plasmid vector downstream of an arabinose-inducible promoter. MTase activity of the CP variants was tested by digestion of the plasmids with methylation-sensitive restriction enzymes. Eleven of the fourteen M.MpeI permutants and six of the seven M.SssI permutants had detectable MTase activity as indicated by the full or partial protection of the plasmid carrying the cpMTase gene. Permutants cp62M.MpeI and cp58M.SssI, in which the new N-termini are located between conserved motifs II and III, had by far the highest activity. The activity of cp62M.MpeI was comparable to the activity of wild-type M.MpeI. Based on the location of the split sites, the permutants possessing MTase activity can be classified in ten types. Although most permutation sites were designed to fall outside of conserved motifs, and the MTase activity of the permutants measured in cell extracts was in most cases substantially lower than that of the wild-type enzyme, the high proportion of circular permutation topologies compatible with MTase activity is remarkable, and is a new evidence for the structural plasticity of C5-MTases. A computer search of the REBASE database identified putative C5-MTases with CP arrangement. Interestingly, all natural circularly permuted C5-MTases appear to represent only one of the ten types of permutation topology created in this work.

Transcriptomic analysis of human norovirus NS1-2 protein highlights a multifunctional role in murine monocytes

Background

The GII.4 Sydney 2012 strain of human norovirus (HuNoV) is a pandemic strain that is responsible for the majority of norovirus outbreaks in healthcare settings. The function of the non-structural (NS)1-2 protein from HuNoV is unknown.

Results

In silico analysis of human norovirus NS1-2 protein showed that it shares features with the murine NS1-2 protein, including a disordered region, a transmembrane domain and H-box and NC sequence motifs. The proteins also contain caspase cleavage and phosphorylation sites, indicating that processing and phosphorylation may be a conserved feature of norovirus NS1-2 proteins. In this study, RNA transcripts of human and murine norovirus full-length and the disordered region of NS1-2 were transfected into monocytes, and next generation sequencing was used to analyse the transcriptomic profile of cells expressing virus proteins. The profiles were then compared to the transcriptomic profile of MNV-infected cells.

Conclusions

RNAseq analysis showed that NS1-2 proteins from human and murine noroviruses affect multiple immune systems (chemokine, cytokine, and Toll-like receptor signaling) and intracellular pathways (NFκB, MAPK, PI3K-Akt signaling) in murine monocytes. Comparison to the transcriptomic profile of MNV-infected cells indicated the pathways that NS1-2 may affect during norovirus infection.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3417-4) contains supplementary material, which is available to authorized users.

Structure of AP205 Coat Protein Reveals Circular Permutation in ssRNA Bacteriophages

AP205 is a single-stranded RNA bacteriophage that has a coat protein sequence not similar to any other known single-stranded RNA phage. Here, we report an atomic-resolution model of the AP205 virus-like particle based on a crystal structure of an unassembled coat protein dimer and a cryo-electron microscopy reconstruction of the assembled particle, together with secondary structure information from site-specific solid-state NMR data. The AP205 coat protein dimer adopts the conserved Leviviridae coat protein fold except for the N-terminal region, which forms a beta-hairpin in the other known single-stranded RNA phages. AP205 has a similar structure at the same location formed by N- and C-terminal beta-strands, making it a circular permutant compared to the other coat proteins. The permutation moves the coat protein termini to the most surface-exposed part of the assembled particle, which explains its increased tolerance to long N- and C-terminal fusions.

Structural Determinants of Misfolding in Multidomain Proteins

Recent single molecule experiments, using either atomic force microscopy (AFM) or Förster resonance energy transfer (FRET) have shown that multidomain proteins containing tandem repeats may form stable misfolded structures. Topology-based simulation models have been used successfully to generate models for these structures with domain-swapped features, fully consistent with the available data. However, it is also known that some multidomain protein folds exhibit no evidence for misfolding, even when adjacent domains have identical sequences. Here we pose the question: what factors influence the propensity of a given fold to undergo domain-swapped misfolding? Using a coarse-grained simulation model, we can reproduce the known propensities of multidomain proteins to form domain-swapped misfolds, where data is available. Contrary to what might be naively expected based on the previously described misfolding mechanism, we find that the extent of misfolding is not determined by the relative folding rates or barrier heights for forming the domains present in the initial intermediates leading to folded or misfolded structures. Instead, it appears that the propensity is more closely related to the relative stability of the domains present in folded and misfolded intermediates. We show that these findings can be rationalized if the folded and misfolded domains are part of the same folding funnel, with commitment to one structure or the other occurring only at a relatively late stage of folding. Nonetheless, the results are still fully consistent with the kinetic models previously proposed to explain misfolding, with a specific interpretation of the observed rate coefficients. Finally, we investigate the relation between interdomain linker length and misfolding, and propose a simple alchemical model to predict the propensity for domain-swapped misfolding of multidomain proteins.

Detection of circular permutations within protein structures using CE-CP

MOTIVATION

Circular permutation is an important type of protein rearrangement. Natural circular permutations have implications for protein function, stability and evolution. Artificial circular permutations have also been used for protein studies. However, such relationships are difficult to detect for many sequence and structure comparison algorithms and require special consideration.

RESULTS

We developed a new algorithm, called Combinatorial Extension for Circular Permutations (CE-CP), which allows the structural comparison of circularly permuted proteins. CE-CP was designed to be user friendly and is integrated into the RCSB Protein Data Bank. It was tested on two collections of circularly permuted proteins. Pairwise alignments can be visualized both in a desktop application or on the web using Jmol and exported to other programs in a variety of formats.

AVAILABILITY AND IMPLEMENTATION

The CE-CP algorithm can be accessed through the RCSB website at http://www.rcsb.org/pdb/workbench/workbench.do. Source code is available under the LGPL 2.1 as part of BioJava 3 (http://biojava.org; http://github.com/biojava/biojava).

CONTACT

sbliven@ucsd.edu or info@rcsb.org.

Cloud4Psi: cloud computing for 3D protein structure similarity searching

Summary: Popular methods for 3D protein structure similarity searching, especially those that generate high-quality alignments such as Combinatorial Extension (CE) and Flexible structure Alignment by Chaining Aligned fragment pairs allowing Twists (FATCAT) are still time consuming. As a consequence, performing similarity searching against large repositories of structural data requires increased computational resources that are not always available. Cloud computing provides huge amounts of computational power that can be provisioned on a pay-as-you-go basis. We have developed the cloud-based system that allows scaling of the similarity searching process vertically and horizontally. Cloud4Psi (Cloud for Protein Similarity) was tested in the Microsoft Azure cloud environment and provided good, almost linearly proportional acceleration when scaled out onto many computational units.

Availability and implementation: Cloud4Psi is available as Software as a Service for testing purposes at: http://cloud4psi.cloudapp.net/. For source code and software availability, please visit the Cloud4Psi project home page at http://zti.polsl.pl/dmrozek/science/cloud4psi.htm.

Contact:dariusz.mrozek@polsl.pl

Systematic detection of internal symmetry in proteins using CE-Symm

Symmetry is an important feature of protein tertiary and quaternary structures that has been associated with protein folding, function, evolution, and stability. Its emergence and ensuing prevalence has been attributed to gene duplications, fusion events, and subsequent evolutionary drift in sequence. This process maintains structural similarity and is further supported by this study. To further investigate the question of how internal symmetry evolved, how symmetry and function are related, and the overall frequency of internal symmetry, we developed an algorithm, CE-Symm, to detect pseudo-symmetry within the tertiary structure of protein chains. Using a large manually curated benchmark of 1007 protein domains, we show that CE-Symm performs significantly better than previous approaches. We use CE-Symm to build a census of symmetry among domain superfamilies in SCOP and note that 18% of all superfamilies are pseudo-symmetric. Our results indicate that more domains are pseudo-symmetric than previously estimated. We establish a number of recurring types of symmetry-function relationships and describe several characteristic cases in detail. With the use of the Enzyme Commission classification, symmetry was found to be enriched in some enzyme classes but depleted in others. CE-Symm thus provides a methodology for a more complete and detailed study of the role of symmetry in tertiary protein structure [availability: CE-Symm can be run from the Web at http://source.rcsb.org/jfatcatserver/symmetry.jsp. Source code and software binaries are also available under the GNU Lesser General Public License (version 2.1) at https://github.com/rcsb/symmetry. An interactive census of domains identified as symmetric by CE-Symm is available from http://source.rcsb.org/jfatcatserver/scopResults.jsp].

Circular permutation of E. coli EPSP synthase: increased inhibitor resistance, improved catalytic activity, and an indicator for protein fragment complementation

We performed the first circular permutation analysis for E. coli 5-enolpyruvylshikimate-3-phosphate synthase, and identified one circular permutant with notably increased resistance to its specific inhibitor and several others with moderately improved catalytic activity. Valid circular permutation sites can be used as effective split sites of protein fragment complementation.

An empirical test of convergent evolution in rhodopsins

Rhodopsins are photochemically reactive membrane proteins that covalently bind retinal chromophores. Type I rhodopsins are found in both prokaryotes and eukaryotic microbes, whereas type II rhodopsins function as photoactivated G-protein coupled receptors (GPCRs) in animal vision. Both rhodopsin families share the seven transmembrane α-helix GPCR fold and a Schiff base linkage from a conserved lysine to retinal in helix G. Nevertheless, rhodopsins are widely cited as a striking example of evolutionary convergence, largely because the two families lack detectable sequence similarity and differ in many structural and mechanistic details. Convergence entails that the shared rhodopsin fold is so especially suited to photosensitive function that proteins from separate origins were selected for this architecture twice. Here we show, however, that the rhodopsin fold is not required for photosensitive activity. We engineered functional bacteriorhodopsin variants with novel folds, including radical noncircular permutations of the α-helices, circular permutations of an eight-helix construct, and retinal linkages relocated to other helices. These results contradict a key prediction of convergence and thereby provide an experimental attack on one of the most intractable problems in molecular evolution: how to establish structural homology for proteins devoid of discernible sequence similarity.

Evolutionary optimization of protein folding

Nature has shaped the make up of proteins since their appearance, 3.8 billion years ago. However, the fundamental drivers of structural change responsible for the extraordinary diversity of proteins have yet to be elucidated. Here we explore if protein evolution affects folding speed. We estimated folding times for the present-day catalog of protein domains directly from their size-modified contact order. These values were mapped onto an evolutionary timeline of domain appearance derived from a phylogenomic analysis of protein domains in 989 fully-sequenced genomes. Our results show a clear overall increase of folding speed during evolution, with known ultra-fast downhill folders appearing rather late in the timeline. Remarkably, folding optimization depends on secondary structure. While alpha-folds showed a tendency to fold faster throughout evolution, beta-folds exhibited a trend of folding time increase during the last 1.5 billion years that began during the “big bang” of domain combinations. As a consequence, these domain structures are on average slow folders today. Our results suggest that fast and efficient folding of domains shaped the universe of protein structure. This finding supports the hypothesis that optimization of the kinetic and thermodynamic accessibility of the native fold reduces protein aggregation propensities that hamper cellular functions.

Nature has come up with an enormous variety of protein three-dimensional structures, each of which is thought to be optimized for its specific function. A fundamental biological endeavor is to uncover the driving evolutionary forces for discovering and optimizing new folds. A long-standing hypothesis is that fold evolution obeys constraints to properly fold into native structure. We here test this hypothesis by analyzing trends of proteins to fold fast during evolution. Using phylogenomic and structural analyses, we observe an overall decrease in folding times between 3.8 and 1.5 billion years ago, which can be interpreted as an evolutionary optimization for rapid folding. This trend towards fast folding probably resulted in manifold advantages, including high protein accessibility for the cell and a reduction of protein aggregation during misfolding.

A primer in macromolecular linguistics

Polymeric macromolecules, when viewed abstractly as strings of symbols, can be treated in terms of formal language theory, providing a mathematical foundation for characterizing such strings both as collections and in terms of their individual structures. In addition this approach offers a framework for analysis of macromolecules by tools and conventions widely used in computational linguistics. This article introduces the ways that linguistics can be and has been applied to molecular biology, covering the relevant formal language theory at a relatively nontechnical level. Analogies between macromolecules and human natural language are used to provide intuitive insights into the relevance of grammars, parsing, and analysis of language complexity to biology.

BioJava: an open-source framework for bioinformatics in 2012

Motivation: BioJava is an open-source project for processing of biological data in the Java programming language. We have recently released a new version (3.0.5), which is a major update to the code base that greatly extends its functionality.

Results: BioJava now consists of several independent modules that provide state-of-the-art tools for protein structure comparison, pairwise and multiple sequence alignments, working with DNA and protein sequences, analysis of amino acid properties, detection of protein modifications and prediction of disordered regions in proteins as well as parsers for common file formats using a biologically meaningful data model.

Availability: BioJava is an open-source project distributed under the Lesser GPL (LGPL). BioJava can be downloaded from the BioJava website (http://www.biojava.org). BioJava requires Java 1.6 or higher. All inquiries should be directed to the BioJava mailing lists. Details are available at http://biojava.org/wiki/BioJava:MailingLists

Contact: andreas.prlic@gmail.com