Circular permutation of betaB2-crystallin changes the hierarchy of domain assembly

Cataract-causing Y204X mutation of crystallin protein CRYβB1 promotes its C-terminal degradation and higher-order oligomerization

Crystallin proteins are a class of main structural proteins of the vertebrate eye lens, and their solubility and stability directly determine transparency and refractive power of the lens. Mutation in genes that encode these crystallin proteins is the most common cause for congenital cataracts. Despite extensive studies, the pathogenic and molecular mechanisms that effect congenital cataracts remain unclear. In this study, we identified a novel mutation in CRYBB1 from a congenital cataract family, and demonstrated that this mutation led to an early termination of mRNA translation, resulting in a 49-residue C-terminally truncated CRYβB1 protein. We show this mutant is susceptible to proteolysis, which allowed us to determine a 1.2-Å resolution crystal structure of CRYβB1 without the entire C-terminal domain. In this crystal lattice, we observed that two N-terminal domain monomers form a dimer that structurally resembles the WT monomer, but with different surface characteristics. Biochemical analyses and cell-based data also suggested that this mutant is significantly more liable to aggregate and degrade compared to WT CRYβB1. Taken together, our results provide an insight into the mechanism regarding how a mutant crystalin contributes to the development of congenital cataract possibly through alteration of inter-protein interactions that result in protein aggregation.

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.

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.

Exploring the folding process of human βB2-crystallin using multiscale molecular dynamics and the Markov state model

Adequate knowledge of protein conformations is crucial for understanding their function and their association properties with other proteins. The cataract disease is correlated with conformational changes in key proteins called crystallins. These changes are due to mutations or post-translational modifications that may lead to protein unfolding, and thus the formation of aggregate states. Human βB2-crystallin (HβB2C) is found in high proportion in the eye lens, and its mutations are related to some cataracts. HβB2C also associates into dimers, tetramers, and other higher-order supramolecular complexes. However, it is the only protein of the βγ-crystallin family that has been found in an extended conformation. Therefore, we hypothesize that the extended conformation is not energetically favourable and that HβB2C may adopt a closed (completely folded) conformation, similar to the other members of the βγ-crystallin family. To corroborate this hypothesis, we performed extensive molecular dynamics simulations of HβB2C in its monomeric and dimeric conformations, using all-atom and coarse-grained scales. We employed Markov state model (MSM) analysis to characterize the conformational and kinetically relevant states in the folding process of monomeric HβB2C. The MSM analysis clearly shows that HβB2C adopts a completely folded structure, and this conformation is the most kinetically and energetically favourable one. In contrast, the extended conformations are kinetically unstable and energetically unfavourable. Our MSM analysis also reveals a key metastable state, which is particularly interesting because it is from this state that the folded state is reached. The folded state is stabilized by the formation of two salt bridges between the residue-pairs E74-R187 and R97-E166 and the two hydrophobic residue-pairs V59-L164 and V72-V151. Furthermore, free energy surface (FES) analysis revealed that the HβB2C dimer with both monomers in a closed conformation (face-en-face dimer) is energetically more stable than the domain-swapped dimer (crystallographic structure). The results presented in this report shed light on the molecular details of the folding mechanism of HβB2C in an aqueous environment and may contribute to interpreting different experimental findings. Finally, a detailed knowledge of HβB2C folding may be key to the rational design of potential molecules to treat cataract disease.

Primary sequence contribution to the optical function of the eye lens

The crystallins have relatively high refractive increments compared to other proteins. The Greek key motif in βγ-crystallins was compared with that in other proteins, using predictive analysis from a protein database, to see whether this may be related to the refractive increment. Crystallins with Greek keys motifs have significantly higher refractive increments and more salt bridges than other proteins with Greek key domains. Specific amino acid substitutions: lysine and glutamic acid residues are replaced by arginine and aspartic acid, respectively as refractive increment increases. These trends are also seen in S-crystallins suggesting that the primary sequence of crystallins may be specifically enriched with amino acids with appropriate values of refractive increment to meet optical requirements. Comparison of crystallins from five species: two aquatic and three terrestrial shows that the lysine/arginine correlation with refractive increment occurs in all species investigated. This may be linked with formation and maintenance of salt bridges.

Microbial βγ-crystallins

βγ-Crystallins have emerged as a superfamily of structurally homologous proteins with representatives across the domains of life. A major portion of this superfamily is constituted by members from microorganisms. This superfamily has also been recognized as a novel group of Ca(2+)-binding proteins with huge diversity. The βγ domain shows variable properties in Ca(2+) binding, stability and association with other domains. The various members present a series of evolutionary adaptations culminating in great diversity in properties and functions. Most of the predicted βγ-crystallins are yet to be characterized experimentally. In this review, we outline the distinctive features of microbial βγ-crystallins and their position in the βγ-crystallin superfamily.

Evolution of crystallins for a role in the vertebrate eye lens

The camera eye lens of vertebrates is a classic example of the re-engineering of existing protein components to fashion a new device. The bulk of the lens is formed from proteins belonging to two superfamilies, the α-crystallins and the βγ-crystallins. Tracing their ancestry may throw light on the origin of the optics of the lens. The α-crystallins belong to the ubiquitous small heat shock proteins family that plays a protective role in cellular homeostasis. They form enormous polydisperse oligomers that challenge modern biophysical methods to uncover the molecular basis of their assembly structure and chaperone-like protein binding function. It is argued that a molecular phenotype of a dynamic assembly suits a chaperone function as well as a structural role in the eye lens where the constraint of preventing protein condensation is paramount. The main cellular partners of α-crystallins, the β- and γ-crystallins, have largely been lost from the animal kingdom but the superfamily is hugely expanded in the vertebrate eye lens. Their structures show how a simple Greek key motif can evolve rapidly to form a complex array of monomers and oligomers. Apart from remaining transparent, a major role of the partnership of α-crystallins with β- and γ-crystallins in the lens is to form a refractive index gradient. Here, we show some of the structural and genetic features of these two protein superfamilies that enable the rapid creation of different assembly states, to match the rapidly changing optical needs among the various vertebrates.

Three-dimensional domain swapping in nitrollin, a single-domain betagamma-crystallin from Nitrosospira multiformis, controls protein conformation and stability but not dimerization

The betagamma-crystallin superfamily has a well-characterized protein fold, with several members found in both prokaryotic and eukaryotic worlds. A majority of them contain two betagamma-crystallin domains. A few examples, such as ciona crystallin and spherulin 3a exist that represent the eukaryotic single-domain proteins of this superfamily. This study reports the high-resolution crystal structure of a single-domain betagamma-crystallin protein, nitrollin, from the ammonium-oxidizing soil bacterium Nitrosospira multiformis. The structure retains the characteristic betagamma-crystallin fold despite a very low sequence identity. The protein exhibits a unique case of homodimerization in betagamma-crystallins by employing its N-terminal extension to undergo three-dimensional (3D) domain swapping with its partner. Removal of the swapped strand results in partial loss of structure and stability but not dimerization per se as determined using gel filtration and equilibrium unfolding studies. Overall, nitrollin represents a distinct single-domain prokaryotic member that has evolved a specialized mode of dimerization hitherto unknown in the realm of betagamma-crystallins.

A comprehensive analysis of non-sequential alignments between all protein structures

Background

The majority of relations between proteins can be represented as a conventional sequential alignment. Nevertheless, unusual non-sequential alignments with different connectivity of the aligned fragments in compared proteins have been reported by many researchers. It is interesting to understand those non-sequential alignments; are they unique, sporadic cases or they occur frequently; do they belong to a few specific folds or spread among many different folds, as a common feature of protein structure. We present here a comprehensive large-scale study of non-sequential alignments between available protein structures in Protein Data Bank.

Results

The study has been conducted on a non-redundant set of 8,865 protein structures aligned with the aid of the TOPOFIT method. It has been estimated that between 17.4% and 35.2% of all alignments are non-sequential depending on variations in the parameters. Analysis of the data revealed that non-sequential relations between proteins do occur systematically and in large quantities. Various sizes and numbers of non-sequential fragments have been observed with all possible complexities of fragment rearrangements found for alignments consisting of up to 12 fragments. It has been found that non-sequential alignments are not limited to proteins of any particular fold and are present in more than two hundred of them. Moreover, many of them are found between proteins with different fold assignments. It has been shown that protein structure symmetry does not explain non-sequential alignments. Therefore, compelling evidences have been provided that non-sequential alignments between proteins are systematic and widespread across the protein universe.

Conclusion

The phenomenon of the widespread occurrence of non-sequential alignments between proteins might represent a missing rule of protein structure organization. More detailed study of this phenomenon will enhance our understanding of protein stability, folding, and evolution.

Protein reconstitution and three-dimensional domain swapping: benefits and constraints of covalency

The phenomena of protein reconstitution and three-dimensional domain swapping reveal that highly similar structures can be obtained whether a protein is comprised of one or more polypeptide chains. In this review, we use protein reconstitution as a lens through which to examine the range of protein tolerance to chain interruptions and the roles of the primary structure in related features of protein structure and folding, including circular permutation, natively unfolded proteins, allostery, and amyloid fibril formation. The results imply that noncovalent interactions in a protein are sufficient to specify its structure under the constraints imposed by the covalent backbone.

Ageing and vision: structure, stability and function of lens crystallins

The alpha-, beta- and gamma-crystallins are the major protein components of the vertebrate eye lens, alpha-crystallin as a molecular chaperone as well as a structural protein, beta- and gamma-crystallins as structural proteins. For the lens to be able to retain life-long transparency in the absence of protein turnover, the crystallins must meet not only the requirement of solubility associated with high cellular concentration but that of longevity as well. For proteins, longevity is commonly assumed to be correlated with long-term retention of native structure, which in turn can be due to inherent thermodynamic stability, efficient capture and refolding of non-native protein by chaperones, or a combination of both. Understanding how the specific interactions that confer intrinsic stability of the protein fold are combined with the stabilizing effect of protein assembly, and how the non-specific interactions and associations of the assemblies enable the generation of highly concentrated solutions, is thus of importance to understand the loss of transparency of the lens with age. Post-translational modification can have a major effect on protein stability but an emerging theme of the few studies of the effect of post-translational modification of the crystallins is one of solubility and assembly. Here we review the structure, assembly, interactions, stability and post-translational modifications of the crystallins, not only in isolation but also as part of a multi-component system. The available data are discussed in the context of the establishment, the maintenance and finally, with age, the loss of transparency of the lens. Understanding the structural basis of protein stability and interactions in the healthy eye lens is the route to solve the enormous medical and economical problem of cataract.

Crystal structure of truncated human betaB1-crystallin

Crystallins are long-lived proteins packed inside eye lens fiber cells that are essential in maintaining the transparency and refractive power of the eye lens. Members of the two-domain betagamma-crystallin family assemble into an array of oligomer sizes, forming intricate higher-order networks in the lens cell. Here we describe the 1.4 angstroms resolution crystal structure of a truncated version of human betaB1 that resembles an in vivo age-related truncation. The structure shows that unlike its close homolog, betaB2-crystallin, the homodimer is not domain swapped, but its domains are paired intramolecularly, as in more distantly related monomeric gamma-crystallins. However, the four-domain dimer resembles one half of the crystallographic bovine betaB2 tetramer and is similar to the engineered circular permuted rat betaB2. The crystal structure shows that the truncated betaB1 dimer is extremely well suited to form higher-order lattice interactions using its hydrophobic surface patches, linker regions, and sequence extensions.

Laser light-scattering evidence for an altered association of beta B1-crystallin deamidated in the connecting peptide

Deamidation is a prevalent modification of crystallin proteins in the vertebrate lens. The effect of specific sites of deamidation on crystallin stability in vivo is not known. Using mass spectrometry, a previously unreported deamidation in beta B1-crystallin was identified at Gln146. Another deamidation was investigated at Asn157. It was determined that whole soluble beta B1 contained 13%-17% deamidation at Gln146 and Asn157. Static and quasi-elastic laser light scattering, circular dichroism, and heat aggregation studies were used to explore the structure and associative properties of recombinantly expressed wild-type (wt) beta B1 and the deamidated beta B1 mutants, Q146E and N157D. Dimer formation occurred for wt beta B1, Q146E, and N157D in a concentration-dependent manner, but only Q146E showed formation of higher ordered oligomers at the concentrations studied. Deamidation at Gln146, but not Asn157, led to an increased tendency of beta B1 to aggregate upon heating. We conclude that deamidation creates unique effects depending upon where the deamidation is introduced in the crystallin structure.

Rapid purification of recombinant betaB2-crystallin using hydrophobic interaction chromatography

BetaB2-crystallin, the major subunit of beta-crystallins, is difficult to purify either from lens homogenate or from betaH-or betaL-crystallins. It has been prepared by heterologous expression in Escherichia coli. Most often, the methods used for purifying a recombinant globular protein employ the combination of ion-exchange with gel filtration chromatography. In the case of betaB2-crystallin too, different approaches have been used to obtain the purified protein, majority of which use a combination of ion-exchange and gel filtration chromatography. We present a new approach to purify betaB2-crystallin using hydrophobic interaction chromatography. In this method, the protein is bound to the hydrophobic matrix in the presence of high concentration of a non-chaotropic salt and eluted by decreasing the salt concentration. The method that we have used for the purification of this globular protein has definite advantages over the earlier methods in its simplicity and efficiency. The most noted advantage of this procedure is the rapid purification with a relatively purified product and a comparatively high yield (>20 mg/L of culture). Over all, the present protocol provides a rapid, efficient and simplified procedure for the preparation of betaB2-crystallin in large yield, sufficient for structural and functional studies.

The evolution of monomeric and oligomeric betagamma-type crystallins. Facts and hypotheses

The case of homologous monomeric gamma-type and oligomeric beta-type crystallins has been described and analyzed in evolutionary terms. Data and hypotheses from molecular genetics and structural investigations converge and suggest a novel three-phase model for the evolutionary history of crystallin-type proteins. In the divergent cascades of monomeric and oligomeric crystallins, a pivotal role was played by alterations in the gene segments encoding the C-terminal extensions and the intermotif or interdomain linker peptides. These were genomic hot spots where evolution experimented to produce the modern variety of betagamma-crystallin-type quaternary structures.

In vivo assembly of aspartate transcarbamoylase from fragmented and circularly permuted catalytic polypeptide chains

Previous studies on Escherichia coli aspartate transcarbamoylase (ATCase) demonstrated that active, stable enzyme was formed in vivo from complementing polypeptides of the catalytic (c) chain encoded by gene fragments derived from the pyrBI operon. However, the enzyme lacked the allosteric properties characteristic of wild-type ATCase. In order to determine whether the loss of homotropic and heterotropic properties was attributable to the location of the interruption in the polypeptide chain rather than to the lack of continuity, we constructed a series of fragmented genes so that the breaks in the polypeptide chains would be dispersed in different domains and diverse regions of the structure. Also, analogous molecules containing circularly permuted c chains with altered termini were constructed for comparison with the ATCase molecules containing fragmented c chains. Studies were performed on four sets of ATCase molecules containing cleaved c chains at positions between residues 98 and 99, 121 and 122, 180 and 181, and 221 and 222; the corresponding circularly permuted chains had N termini at positions 99, 122, 181, and 222. All of the ATCase molecules containing fragmented or circularly permuted c chains exhibited the homotropic and heterotropic properties characteristic of the wild-type enzyme. Hill coefficients (n(H:)) and changes in them upon the addition of ATP and CTP were similar to those observed with wild-type ATCase. In addition, the conformational changes revealed by the decrease in sedimentation coefficient upon the addition of a bisubstrate analog were virtually identical to that for the wild-type enzyme. Differential scanning calorimetry showed that neither the breakage of the polypeptide chains nor the newly formed covalent bond between the termini in the wild-type enzyme had a significant impact on the thermal stability of the assembled dodecamers. The studies demonstrate that continuity of the polypeptide chain within structural domains is not essential for the assembly, activity, and allosteric properties of ATCase.

Structure of the crystallins

The lens is formed from two protein superfamilies, the alpha- and beta gamma-crystallins. Representative three-dimensional structures show they both have a basic 2-beta-sheet domain fold, with the beta gamma-domain being made from two intercalating Greek keys. X-ray structures of monomeric gamma-crystallins and simple oligomeric beta-crystallins show how multiple gene duplications can give rise to highly symmetrical assemblies based on paired domains. These protein folds have been engineered by directed mutagenesis to investigate the roles of the critical region in domain pairing and assembly. Inherited human cataracts have been described that are associated with representatives of each of the crystallin protein families. Mutations to certain beta- and gamma-crystallin genes cause expression of truncated polypeptides that would not be expected to fold properly; instead they would randomly aggregate causing light scattering. As crystallin proteins are not renewed, age-related cataract is a gradual accumulation of small changes to pre-existing normal proteins. The precise sites of post-translational modifications are now being mapped to the various crystallins.

Eye lens betaB2-crystallin: circular permutation does not influence the oligomerization state but enhances the conformational stability

The related vertebrate eye lens polypeptides, betaB2- and gammaB-crystallin, each fold into two similar beta-sheet domains. The main difference is the state of oligomerization resulting from intermolecular domain interactions in the oligomeric beta-crystallins and intramolecular contacts in the monomeric gamma-crystallins. The question arises whether it is possible to create a monomeric gammaB-like betaB2-molecule by protein engineering methods. We wanted to produce such a molecule by circularly permuting the domains of betaB2-crystallin. The new termini were created from the original connecting peptide, and the new linker from stumps of the original extensions, while the rest of the flexible extensions were deleted. As judged by circular dichroism and fluorescence, the permutation causes little change in the structure of the protein. The circularly permuted protein forms dimers as wild-type betaB2-crystallin. On the other hand, cpbetaB2 shows a slightly enhanced stability against urea with a midpoint of transition of 2.1 M urea versus 1.9 M for the wild-type protein lacking N and C-terminal arms, thus indicating stronger domain interactions. To our knowledge this is the first circularly permuted protein which exhibits a higher stability than the corresponding wild-type protein.