Mapping the Distribution of Conformational Information Throughout a Protein Sequence

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/.

Discovering the Ultimate Limits of Protein Secondary Structure Prediction

Secondary structure prediction (SSP) of proteins is an important structural biology technique with many applications. There have been ~300 algorithms published in the past seven decades with fierce competition in accuracy. In the first 60 years, the accuracy of three-state SSP rose from ~56% to 81%; after that, it has long stayed at 81-86%. In the 1990s, the theoretical limit of three-state SSP accuracy had been estimated to be 88%. Thus, SSP is now generally considered not challenging or too challenging to improve. However, we found that the limit of three-state SSP might be underestimated. Besides, there is still much room for improving segment-based and eight-state SSPs, but the limits of these emerging topics have not been determined. This work performs large-scale sequence and structural analyses to estimate SSP accuracy limits and assess state-of-the-art SSP methods. The limit of three-state SSP is re-estimated to be ~92%, 4-5% higher than previously expected, indicating that SSP is still challenging. The estimated limit of eight-state SSP is 84-87%. Several proposals for improving future SSP algorithms are made based on our results. We hope that these findings will help move forward the development of SSP and all its applications.

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.

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.

A secondary structure-based position-specific scoring matrix applied to the improvement in protein secondary structure prediction

Protein secondary structure prediction (SSP) has a variety of applications; however, there has been relatively limited improvement in accuracy for years. With a vision of moving forward all related fields, we aimed to make a fundamental advance in SSP. There have been many admirable efforts made to improve the machine learning algorithm for SSP. This work thus took a step back by manipulating the input features. A secondary structure element-based position-specific scoring matrix (SSE-PSSM) is proposed, based on which a new set of machine learning features can be established. The feasibility of this new PSSM was evaluated by rigid independent tests with training and testing datasets sharing <25% sequence identities. In all experiments, the proposed PSSM outperformed the traditional amino acid PSSM. This new PSSM can be easily combined with the amino acid PSSM, and the improvement in accuracy was remarkable. Preliminary tests made by combining the SSE-PSSM and well-known SSP methods showed 2.0% and 5.2% average improvements in three- and eight-state SSP accuracies, respectively. If this PSSM can be integrated into state-of-the-art SSP methods, the overall accuracy of SSP may break the current restriction and eventually bring benefit to all research and applications where secondary structure prediction plays a vital role during development. To facilitate the application and integration of the SSE-PSSM with modern SSP methods, we have established a web server and standalone programs for generating SSE-PSSM available at http://10.life.nctu.edu.tw/SSE-PSSM.

The structure of unliganded sterol carrier protein 2 from Yarrowia lipolytica unveils a mechanism for binding site occlusion

Isolated or as a part of multidomain proteins, Sterol Carrier Protein 2 (SCP2) exhibits high affinity and broad specificity for different lipidic and hydrophobic compounds. A wealth of structural information on SCP2 domains in all forms of life is currently available; however, many aspects of its ligand binding activity are poorly understood. ylSCP2 is a well-characterized single domain SCP2 from the yeast Yarrowia lipolytica. Herein, we report the X-ray structure of unliganded ylSCP2 refined to 2.0 Å resolution. Comparison with the previously solved liganded ylSCP2 structure unveiled a novel mechanism for binding site occlusion. The liganded ylSCP2 binding site is a large cavity with a volume of more than 800 Å³. In unliganded ylSCP2 the binding site is reduced to about 140 Å³. The obliteration is caused by a swing movement of the C-terminal α helix 5 and a subtle compaction of helices 2-4. Previous pairwise comparisons were between homologous SCP2 domains with a uncertain binding status. The reported unliganded ylSCP2 structure allows for the first time a fully controlled comparative analysis of the conformational effects of ligand occupation dispelling several doubts regarding the architecture of SCP2 binding site.

Equilibrium partially folded states of B. licheniformis[Formula: see text]-lactamase

[Formula: see text]-Lactamases (penicillinases) facilitate bacterial resistance to antibiotics and are excellent theoretical and experimental models in protein structure, dynamics and evolution. Bacillus licheniformis exo-small penicillinase (ESP) is a Class A [Formula: see text]-lactamase with three tryptophan residues located one in each of its two domains and one in the interface between domains. The conformational landscape of three well-characterized ESP Trp[Formula: see text]Phe mutants was characterized in equilibrium unfolding experiments by measuring tryptophan fluorescence, far-UV CD, activity, hydrodynamic radius, and limited proteolysis. The Trp[Formula: see text]Phe substitutions had little impact on the native conformation, but changed the properties of the partially folded states populated at equilibrium. The results were interpreted in the framework of modern theories of protein folding.

Engineering β-lactamase zymogens for use in protease activity assays

In this report, we engineered β-lactamase zymogens and developed a sensitive protease assay method based on the precursor enzymes.

Quantifying protein modularity and evolvability: a comparison of different techniques

Modularity increases evolvability by reducing constraints on adaptation and by allowing preexisting parts to function in new contexts for novel uses. Protein evolution provides an excellent context to study the causes and consequences of biological modularity. In order to address such questions, however, an index for protein modularity is necessary. This paper proposes a simple index for protein modularity-"module density"-which is the number of evolutionarily independent modules that compose a protein divided by the number of amino acids in the protein. The decomposition of proteins into constituent modules can be accomplished by either of two classes of methods. The first class of methods relies on "suppositional" criteria to assign amino acids to modules, whereas the second class of methods relies on "coevolutionary" criteria for this task. One simple and practical method from the first class consists of approximating the number of modules in a protein as the number of regular secondary structure elements (i.e., helices and sheets). Methods based on coevolutionary criteria require more elaborate data, but they have the advantage of being able to specify modules without prior assumptions about why they exist. Given the increasing availability of datasets sampling protein mutational spectra (e.g., from comparative genomics, experimental evolution, and computational prediction), methods based on coevolutionary criteria will likely become more promising in the near future. The ability to meaningfully quantify protein modularity via simple indices has the potential to aid future efforts to understand protein evolutionary rate determinants, improve molecular evolution models and engineer novel proteins.

CPred: a web server for predicting viable circular permutations in proteins

Circular permutation (CP) is a protein structural rearrangement phenomenon, through which nature allows structural homologs to have different locations of termini and thus varied activities, stabilities and functional properties. It can be applied in many fields of protein research and bioengineering. The limitation of applying CP lies in its technical complexity, high cost and uncertainty of the viability of the resulting protein variants. Not every position in a protein can be used to create a viable circular permutant, but there is still a lack of practical computational tools for evaluating the positional feasibility of CP before costly experiments are carried out. We have previously designed a comprehensive method for predicting viable CP cleavage sites in proteins. In this work, we implement that method into an efficient and user-friendly web server named CPred (CP site predictor), which is supposed to be helpful to promote fundamental researches and biotechnological applications of CP. The CPred is accessible at http://sarst.life.nthu.edu.tw/CPred.

X-ray evidence of a native state with increased compactness populated by tryptophan-less B. licheniformis β-lactamase

β-lactamases confer antibiotic resistance, one of the most serious world-wide health problems, and are an excellent theoretical and experimental model in the study of protein structure, dynamics and evolution. Bacillus licheniformis exo-small penicillinase (ESP) is a Class-A β-lactamase with three tryptophan residues located in the protein core. Here, we report the 1.7-Å resolution X-ray structure, catalytic parameters, and thermodynamic stability of ESP(ΔW), an engineered mutant of ESP in which phenylalanine replaces the wild-type tryptophan residues. The structure revealed no qualitative conformational changes compared with thirteen previously reported structures of B. licheniformis β-lactamases (RMSD = 0.4-1.2 Å). However, a closer scrutiny showed that the mutations result in an overall more compact structure, with most atoms shifted toward the geometric center of the molecule. Thus, ESP(ΔW) has a significantly smaller radius of gyration (R(g)) than the other B. licheniformis β-lactamases characterized so far. Indeed, ESP(ΔW) has the smallest R(g) among 126 Class-A β-lactamases in the Protein Data Bank (PDB). Other measures of compactness, like the number of atoms in fixed volumes and the number and average of noncovalent distances, confirmed the effect. ESP(ΔW) proves that the compactness of the native state can be enhanced by protein engineering and establishes a new lower limit to the compactness of the Class-A β-lactamase fold. As the condensation achieved by the native state is a paramount notion in protein folding, this result may contribute to a better understanding of how the sequence determines the conformational variability and thermodynamic stability of a given fold.

Deciphering the preference and predicting the viability of circular permutations in proteins

Circular permutation (CP) refers to situations in which the termini of a protein are relocated to other positions in the structure. CP occurs naturally and has been artificially created to study protein function, stability and folding. Recently CP is increasingly applied to engineer enzyme structure and function, and to create bifunctional fusion proteins unachievable by tandem fusion. CP is a complicated and expensive technique. An intrinsic difficulty in its application lies in the fact that not every position in a protein is amenable for creating a viable permutant. To examine the preferences of CP and develop CP viability prediction methods, we carried out comprehensive analyses of the sequence, structural, and dynamical properties of known CP sites using a variety of statistics and simulation methods, such as the bootstrap aggregating, permutation test and molecular dynamics simulations. CP particularly favors Gly, Pro, Asp and Asn. Positions preferred by CP lie within coils, loops, turns, and at residues that are exposed to solvent, weakly hydrogen-bonded, environmentally unpacked, or flexible. Disfavored positions include Cys, bulky hydrophobic residues, and residues located within helices or near the protein's core. These results fostered the development of an effective viable CP site prediction system, which combined four machine learning methods, e.g., artificial neural networks, the support vector machine, a random forest, and a hierarchical feature integration procedure developed in this work. As assessed by using the hydrofolate reductase dataset as the independent evaluation dataset, this prediction system achieved an AUC of 0.9. Large-scale predictions have been performed for nine thousand representative protein structures; several new potential applications of CP were thus identified. Many unreported preferences of CP are revealed in this study. The developed system is the best CP viability prediction method currently available. This work will facilitate the application of CP in research and biotechnology.

Optical studies of single-tryptophan B. licheniformis beta-lactamase variants

beta-lactamases (penicillinases) are important complicating factors in bacterial infections and excellent theoretical and experimental models in protein structure, dynamics and evolution. Bacillus licheniformis exo-small penicillinase (ESP) is a Class A beta-lactamase with three tryptophan residues, one located in each of the two protein domains and one located in the interface between domains. To determine the tryptophan contribution to the ESP UV-absorption, circular dichroism, and steady-state and time-resolved fluorescence, four Trp-->Phe mutants were prepared and characterized. The residue substitutions had little impact on the native conformation. UV-absorption and CD features were identified and ascribed to specific aromatic residues. Time-resolved fluorescence showed that most of the fluorescence decay of ESP tryptophans is due to a discrete exponential component with a lifetime of 5-6ns. Fluorescence polarization measurements indicated that fluorescence of Trp 210 is nearly independent of the fluorescence of Trp 229 and Trp 251, whereas a substantial energy homotransfer between the latter pair takes place. The spectroscopic information was rationalized on the basis of structural considerations and should help in the interpretation and monitoring of the changes at the sub domain level during the conformational transitions and fluctuations of ESP and other Class A beta-lactamases.

Analysis of p300/CBP histone acetyltransferase regulation using circular permutation and semisynthesis

The histone acetyltransferase (HAT) p300/CBP has been shown to undergo autoacetylation on lysines in an apparent regulatory loop that stimulates HAT activity. Here we have developed a strategy to introduce acetyl-Lys at up to six known modification sites in p300/CBP HAT using a combination of circular permutation and expressed protein ligation. We show that these semisynthetic, circularly permuted acetylated proteins retain high affinity for an acetyl-CoA substrate analogue and that HAT activity correlates positively with degree of acetylation. This study provides novel evidence for control of p300/CBP HAT activity by site-specific autoacetylation and outlines a potentially general strategy for using expressed protein ligation and circular permutation to chemically interrogate internal regions of proteins.

CPDB: a database of circular permutation in proteins

Circular permutation (CP) in a protein can be considered as if its sequence were circularized followed by a creation of termini at a new location. Since the first observation of CP in 1979, a substantial number of studies have concluded that circular permutants (CPs) usually retain native structures and functions, sometimes with increased stability or functional diversity. Although this interesting property has made CP useful in many protein engineering and folding researches, large-scale collections of CP-related information were not available until this study. Here we describe CPDB, the first CP DataBase. The organizational principle of CPDB is a hierarchical categorization in which pairs of circular permutants are grouped into CP clusters, which are further grouped into folds and in turn classes. Additions to CPDB include a useful set of tools and resources for the identification, characterization, comparison and visualization of CP. Besides, several viable CP site prediction methods are implemented and assessed in CPDB. This database can be useful in protein folding and evolution studies, the discovery of novel protein structural and functional relationships, and facilitating the production of new CPs with unique biotechnical or industrial interests. The CPDB database can be accessed at http://sarst.life.nthu.edu.tw/cpdb

CPSARST: an efficient circular permutation search tool applied to the detection of novel protein structural relationships

CPSARST (Circular Permutation Search Aided by Ramachandran Sequential Transformation) is an efficient database search tool that provides a new way for rapidly detecting novel relationships among proteins.

Circular permutation of a protein can be visualized as if the original amino- and carboxyl termini were linked and new ones created elsewhere. It has been well-documented that circular permutants usually retain native structures and biological functions. Here we report CPSARST (Circular Permutation Search Aided by Ramachandran Sequential Transformation) to be an efficient database search tool. In this post-genomics era, when the amount of protein structural data is increasing exponentially, it provides a new way to rapidly detect novel relationships among proteins.

Identification and thermodynamic characterization of molten globule states of periplasmic binding proteins

Molten globule-like intermediates have been shown to occur during protein folding and are thought to be involved in protein translocation and membrane insertion. However, the determinants of molten globule stability and the extent of specific packing in molten globules is currently unclear. Using far- and near-UV CD and intrinsic and ANS fluorescence, we show that four periplasmic binding proteins (LBP, LIVBP, MBP, and RBP) form molten globules at acidic pH values ranging from 3.0 to 3.4. Only two of these (LBP and LIVBP) have similar sequences, but all four proteins adopt similar three-dimensional structures. We found that each of the four molten globules binds to its corresponding ligand without conversion to the native state. Ligand binding affinity measured by isothermal titration calorimetry for the molten globule state of LIVBP was found to be comparable to that of the corresponding native state, whereas for LBP, MBP, and RBP, the molten globules bound ligand with approximately 5-30-fold lower affinity than the corresponding native states. All four molten globule states exhibited cooperative thermal unfolding assayed by DSC. Estimated values of DeltaCp of unfolding show that these molten globule states contain 28-67% of buried surface area relative to the native states. The data suggest that molten globules of these periplasmic binding proteins retain a considerable degree of long range order. The ability of these sequentially unrelated proteins to form highly ordered molten globules may be related to their large size as well as an intrinsic property of periplasmic binding protein folds.

Effects of serine-to-cysteine mutations on beta-lactamase folding

B. licheniformis exo-small beta-lactamase (ESBL) has two nonsequential domains and a complex architecture. We replaced ESBL serine residues 126 and 265 with cysteine to probe the conformation of buried regions in each domain. Spectroscopic, hydrodynamic, and chemical methods revealed that the mutations do not alter the native fold but distinctly change stability (S-126C > wild-type > S-126/265C > S-265C ESBL) and the features of partially folded states. The observed wild-type ESBL equilibrium intermediate has decreased fluorescence but full secondary structure. S-126C ESBL intermediate has the fluorescence of the unfolded state, no thiol reactivity, and partial secondary structure. S-265C and S-126/265C ESBL populate intermediate states unfolded by fluorescence and thiol reactivity but with full secondary structure. Mass analysis of S-126/265C ESBL in the partially folded state proved that both thiol groups become exposed simultaneously. None of the intermediates is compatible with sequential domain unfolding. Molecular dynamics simulation suggests that the stabilizing effect of the S-126C substitution is due to optimization of van der Waals interactions and packing. On the other hand, destabilization induced by the S-265C mutation results from alteration of the hydrogen-bond network. The results illustrate the large impact that seemingly conservative serine-to-cysteine changes can have on the energy landscape of proteins.