Stability and folding of domain proteins

Targeting the protein folding transition state by mutation: Large scale (un)folding rate accelerations without altering native stability

Proteins are constantly undergoing folding and unfolding transitions, with rates that determine their homeostasis in vivo and modulate their biological function. The ability to optimize these rates without affecting overall native stability is hence highly desirable for protein engineering and design. The great challenge is, however, that mutations generally affect folding and unfolding rates with inversely complementary fractions of the net free energy change they inflict on the native state. Here we address this challenge by targeting the folding transition state (FTS) of chymotrypsin inhibitor 2 (CI2), a very slow and stable two-state folding protein with an FTS known to be refractory to change by mutation. We first discovered that the CI2's FTS is energetically taxed by the desolvation of several, highly conserved, charges that form a buried salt bridge network in the native structure. Based on these findings, we designed a CI2 variant that bears just four mutations and aims to selectively stabilize the FTS. This variant has >250-fold faster rates in both directions and hence identical native stability, demonstrating the success of our FTS-centric design strategy. With an optimized FTS, CI2 also becomes 250-fold more sensitive to proteolytic degradation by its natural substrate chymotrypsin, and completely loses its activity as inhibitor. These results indicate that CI2 has been selected through evolution to have a very unstable FTS in order to attain the kinetic stability needed to effectively function as protease inhibitor. Moreover, the CI2 case showcases that protein (un)folding rates can critically pivot around a few key residues-interactions, which can strongly modify the general effects of known structural factors such as domain size and fold topology. From a practical standpoint, our results suggest that future efforts should perhaps focus on identifying such critical residues-interactions in proteins as best strategy to significantly improve our ability to predict and engineer protein (un)folding rates.

Conflicting Interfacial Electrostatic Interactions as a Design Principle to Modulate Long-Range Interdomain Communication

The extent and molecular basis of interdomain communication in multidomain proteins, central to understanding allostery and function, is an open question. One simple evolutionary strategy could involve the selection of either conflicting or favorable electrostatic interactions across the interface of two closely spaced domains to tune the magnitude of interdomain connectivity. Here, we study a bilobed domain FF34 from the eukaryotic p190A RhoGAP protein to explore one such design principle that mediates interdomain communication. We find that while the individual structural units in wild-type FF34 are marginally coupled, they exhibit distinct intrinsic stabilities and low cooperativity, manifesting as slow folding. The FF3-FF4 interface harbors a frustrated network of highly conserved electrostatic interactions-a charge troika-that promotes the population of multiple, decoupled, and non-native structural modes on a rugged native landscape. Perturbing this network via a charge-reversal mutation not only enhances stability and cooperativity but also dampens the fluctuations globally and speeds up the folding rate by at least an order of magnitude. Our work highlights how a conserved but nonoptimal network of interfacial electrostatic interactions shapes the native ensemble of a bilobed protein, a feature that could be exploited in designing molecular systems with long-range connectivity and enhanced cooperativity.

Living in trinity of extremes: Genomic and proteomic signatures of halophilic, thermophilic, and pH adaptation

Since nucleic acids and proteins of unicellular prokaryotes are directly exposed to extreme environmental conditions, it is possible to explore the genomic-proteomic compositional determinants of molecular mechanisms of adaptation developed by them in response to harsh environmental conditions. Using a wealth of currently available complete genomes/proteomes we were able to explore signatures of adaptation to three environmental factors, pH, salinity, and temperature, observing major trends in compositions of their nucleic acids and proteins. We derived predictors of thermostability, halophilic, and pH adaptations and complemented them by the principal components analysis. We observed a clear difference between thermophilic and salinity/pH adaptations, whereas latter invoke seemingly overlapping mechanisms. The genome-proteome compositional trade-off reveals an intricate balance between the work of base paring and base stacking in stabilization of coding DNA and r/tRNAs, and, at the same time, universal requirements for the stability and foldability of proteins regardless of the nucleotide biases. Nevertheless, we still found hidden fingerprints of ancient evolutionary connections between the nucleotide and amino acid compositions indicating their emergence, mutual evolution, and adjustment. The evolutionary perspective on the adaptation mechanisms is further studied here by means of the comparative analysis of genomic/proteomic traits of archaeal and bacterial species. The overall picture of genomic/proteomic signals of adaptation obtained here provides a foundation for future engineering and design of functional biomolecules resistant to harsh environments.

Bioinformatic analysis of the effect of SNPs in the pig TERT gene on the structural and functional characteristics of the enzyme to develop new genetic markers of productivity traits

BACKGROUND

Telomerase reverse transcriptase (TERT) plays a crucial role in synthesizing telomeric repeats that safeguard chromosomes from damage and fusion, thereby maintaining genome stability. Mutations in the TERT gene can lead to a deviation in gene expression, impaired enzyme activity, and, as a result, abnormal telomere shortening. Genetic markers of productivity traits in livestock can be developed based on the TERT gene polymorphism for use in marker-associated selection (MAS). In this study, a bioinformatic-based approach is proposed to evaluate the effect of missense single-nucleotide polymorphisms (SNPs) in the pig TERT gene on enzyme function and structure, with the prospect of developing genetic markers.

RESULTS

A comparative analysis of the coding and amino acid sequences of the pig TERT was performed with corresponding sequences of other species. The distribution of polymorphisms in the pig TERT gene, with respect to the enzyme's structural-functional domains, was established. A three-dimensional model of the pig TERT structure was obtained through homological modeling. The potential impact of each of the 23 missense SNPs in the pig TERT gene on telomerase function and stability was assessed using predictive bioinformatic tools utilizing data on the amino acid sequence and structure of pig TERT.

CONCLUSIONS

According to bioinformatic analysis of 23 missense SNPs of the pig TERT gene, a predictive effect of rs789641834 (TEN domain), rs706045634 (TEN domain), rs325294961 (TRBD domain) and rs705602819 (RTD domain) on the structural and functional parameters of the enzyme was established. These SNPs hold the potential to serve as genetic markers of productivity traits. Therefore, the possibility of their application in MAS should be further evaluated in associative analysis studies.

Ordered-domain unfolding of thermophilic isolated β subunit ATP synthase

The flexibility of the ATP synthase's β subunit promotes its role in the ATP synthase rotational mechanism, but its domains stability remains unknown. A reversible thermal unfolding of the isolated β subunit (Tβ) of the ATP synthase from Bacillus thermophilus PS3, tracked through circular dichroism and molecular dynamics, indicated that Tβ shape transits from an ellipsoid to a molten globule through an ordered unfolding of its domains, preserving the β-sheet residual structure at high temperature. We determined that part of the stability origin of Tβ is due to a transversal hydrophobic array that crosses the β-barrel formed at the N-terminal domain and the Rossman fold of the nucleotide-binding domain (NBD), while the helix bundle of the C-terminal domain is the less stable due to the lack of hydrophobic residues, and thus the more flexible to trigger the rotational mechanism of the ATP synthase.

Nano-Formulation of Antioxidants as Effective Inhibitors of γD-Crystallin Aggregation

The aggregation of crystallin proteins is related to cataracts and age-related macular degeneration. Apart from surgical replacement of the cataract lens, no other alternative treatment is available till date for this ailment. In the current work, we carried out an in-depth investigation of the effect of polyphenol-loaded nano-formulations on the aggregation of γD-crystallin. At first, the protein was allowed to form amorphous aggregates under denaturing conditions. Several polyphenols were then tried to inhibit the aggregation of the protein. Among the polyphenols tested, resveratrol and quercetin were found to be the most effective. Since polyphenols are prone to degradation, they were encapsulated in chitosan nanoparticles in order to provide ambient conditions for them to function effectively. The loading efficiency and polyphenol release kinetics were subsequently tested. Finally, the efficacy of resveratrol/quercetin-loaded chitosan nano-particles as inhibitors of γD-crystallin aggregation was confirmed in a series of experiments demonstrating the potency of the system in the prospective therapeutic intervention of eye ailments concerning self-assembly of γD-crystallin proteins.

Engineering functional thermostable proteins using ancestral sequence reconstruction

Natural proteins are often only slightly more stable in the native state than the denatured state, and an increase in environmental temperature can easily shift the balance towards unfolding. Therefore, the engineering of proteins to improve protein stability is an area of intensive research. Thermostable proteins are required to withstand industrial process conditions, for increased shelf-life of protein therapeutics, for developing robust 'biobricks' for synthetic biology applications, and for research purposes (e.g. structure determination). In addition, thermostability buffers the often destabilizing effects of mutations introduced to improve other properties. Rational design approaches to engineering thermostability require structural information, but even with advanced computational methods, it is challenging to predict or parameterize all the relevant structural factors with sufficient precision to anticipate the results of a given mutation. Directed evolution is an alternative when structures are unavailable but requires extensive screening of mutant libraries. Recently however, bioinspired approaches based on phylogenetic analyses have shown great promise. Leveraging the rapid expansion in sequence data and bioinformatic tools, ancestral sequence reconstruction (ASR) can generate highly stable folds for novel applications in industrial chemistry, medicine, and synthetic biology. This review provides an overview of the factors important for successful inference of thermostable proteins by ASR and what it can reveal about the determinants of stability in proteins.

Co-Translational Folding of Multi-Domain Proteins

The majority of proteins in nature are composed of multiple domains connected in a single polypeptide. How these long sequences fold into functional structures without forming toxic misfolds or aggregates is poorly understood. Their folding is inextricably linked to protein synthesis and interactions with cellular machinery, making mechanistic studies challenging. Recent progress has revealed critical features of multi-domain protein folding in isolation and in the context of translation by the ribosome. In this review, we discuss challenges and progress in understanding multi-domain protein folding, and highlight how molecular interactions shape folding and misfolding pathways. With the development of new approaches and model systems, the stage is now set for mechanistically exploring the folding of large multi-domain proteins.

Insights to Human γD-Crystallin Unfolding by NMR Spectroscopy and Molecular Dynamics Simulations

Human γD-crystallin (HGDC) is an abundant lens protein residing in the nucleus of the human lens. Aggregation of this and other structural proteins within the lens leads to the development of cataract. Much has been explored on the stability and aggregation of HGDC and where detailed investigation at the atomic resolution was needed, the X-ray structure was used as an initial starting conformer for molecular modeling. In this study, we implemented NMR-solution HGDC structures as starting conformers for molecular dynamics simulations to provide the missing pieces of the puzzle on the very early stages of HGDC unfolding leading up to the domain swap theories proposed by past studies. The high-resolution details of the conformational dynamics also revealed additional insights to possible early intervention for cataractogenesis.

Resurrection and characterization of ancestral CYP11A1 enzymes

Mitochondrial cytochromes P450 presumably originated from a common microsomal P450 ancestor. However, it is still unknown how ancient mitochondrial P450s were able to retain their oxygenase function following relocation to the mitochondrial matrix and later emerged as enzymes specialized for steroid hormone biosynthesis in vertebrates. Here, we used the approach of ancestral sequence reconstruction (ASR) to resurrect ancient CYP11A1 enzymes and characterize their unique biochemical properties. Two ancestral CYP11A1 variants, CYP11A_Mammal_N101 and CYP11A_N1, as well as an extant bovine form were recombinantly expressed and purified to homogeneity. All enzymes showed characteristic P450 spectral properties and were able to convert cholesterol as well as other sterol substrates to pregnenolone, yet with different specificities. The vertebrate CYP11A_N1 ancestor preferred the cholesterol precursor, desmosterol, as substrate suggesting a convergent evolution of early cholesterol metabolism and CYP11A1 enzymes. Both ancestors were able to withstand increased levels of hydrogen peroxide but only the ancestor CYP11A_N1 showed increased thermostability (~ 25 °C increase in T₅₀ ) compared with the extant CYP11A1. The extraordinary robustness of ancient mitochondrial P450s, as demonstrated for CYP11A_N1, may have allowed them to stay active when presented with poorly compatible electron transfer proteins and resulting harmful ROS in the new environment of the mitochondrial matrix. To the best of our knowledge, this work represents the first study that describes the resurrection of ancient mitochondrial P450 enzymes. The results will help to understand and gain fundamental functional insights into the evolutionary origins of steroid hormone biosynthesis in animals.

The concept of protein folding/unfolding and its impacts on human health

Proteins have evolved in specific 3D structures and play different functions in cells and determine various reactions and pathways. The newly synthesized amino acid chains once depart ribosome must crumple into three-dimensional structures so can be biologically active. This process of protein that makes a functional molecule is called protein folding. The protein folding is both a biological and a physicochemical process that depends on the sequence of it. In fact, this process occurs more complicated and in some cases and in exposure to some molecules like glucose (glycation), mistaken folding leads to amyloid structures and fatal disorders called conformational diseases. Such conditions are detected by the quality control system of the cell and these abnormal proteins undergo renovation or degradation. This scenario takes place by the chaperones, chaperonins, and Ubiquitin-proteasome complex. Understanding of protein folding mechanisms from different views including experimental and computational approaches has revealed some intermediate ensembles such as molten globule and has been subjected to biophysical and molecular biology attempts to know more about prevalent conformational diseases.

Characterization and immobilization of protease secreted by the fungus Moorella speciosa

Protease was extracellularly produced in submerged fermentation by the fungus Moorella speciosa with maximum activity of 8.6 × 10³ U/mL. The optimal pH and temperature for enzyme activity were 6.78 and 60.88 °C, respectively. The enzyme was incubated in the presence of several ions at concentrations of 0.1 M and 0.01 M to address the effect on enzyme activity. Enzyme activity was increased by 56% and 130% in the presence of 0.1 M BaCl₂ and of 0.01 M Na₂SO₄, respectively. The V _max and K _m values were 0.01474 U/min/mg protein and 0.04190 mg/mL, respectively. The enzyme retained about 90% of enzymatic activity at 90 °C. Among the methods tested for enzyme immobilization, adsorption onto MAT540 carrier led to the most promising results, since after 15 reuse cycles up to 60% of the initial catalytic activity was retained. Entrapment in calcium alginate matrix allowed to retain up to 51% of the initial catalytic activity after 8 reuse cycles. This protease from M. speciosa, in either free or immobilized form, can be foreseen as a useful biocatalytic tool in process design by reducing operating costs, decreasing the use of chemical processing and, consequently, meeting the global demand for clean technologies.

Arachnid Hemocyanins

Hemocyanin (Hc), a copper-containing extracellular multimeric protein, is the major protein component of hemolymph in different arachnid groups. Hc possesses 7 or 8 very well-characterized types of monomers with molecular weights ranging from 70 to 85 kDa, organized in hexamers or multiple of hexamers. The present chapter compiles the existing data with relation to the function of this protein in the arachnids. Hc has as main function the reversible transport of O₂, but it shows many secondary though not less important functions. With reference to this, it has been described that Hc can transport hydrophobic molecules (lipid-derived hormones and lipids) to the different organs, having a key role in the lipid transport system. In arachnids, like in other arthropods and invertebrates, Hc has phenoloxidase function which is related to different metabolic processes such as melanin formation and defense against pathogens. In addition, Hc has additional defensive functions since it can serve as precursor for the production of antimicrobial peptides. In short, the evolution of this protein has led to the development of multiple functions essential for organisms possessing this protein.

Crystallin Fusion Proteins Improve the Thermal Properties of Hair

Styling hair with straightening irons is a popular daily hair routine that significantly damage the hair keratin fiber due to the high temperature applied. In this study, we investigate the effect of two fusion proteins based on the human eye γD-crystallin conjugated with a keratin-based peptide (KP-Cryst Wt and KP-Cryst Mut) on hair exposed to thermal damage. The mutant form was designed to improve protein stability and promote interaction with the hair. Through the study, it was demonstrated the protection of Asian and Caucasian virgin hair's structure by the pretreatments with the KP-Cryst fusion proteins. After hair thermal exposure, a higher water content was quantified by TGA on the hair fibers pretreated with the fusion proteins (about 38% for the KP-Cryst Wt and 44% for the KP-Cryst Mut). Also, negligible alterations in hair fibers' stiffness were observed after iron application, demonstrating the proteins capacity to effectively prevent the conversion of keratin α-helix structure into β-sheets. The results proved the capacity of the fusion proteins to bind to hair and protect it against high temperatures', supporting the development of new formulations based on the KP-Cryst proteins.

Vitellogenesis in spiders: first analysis of protein changes in different reproductive stages of Polybetes pythagoricus

Vitellogenesis represents one of the most vital processes of oviparous species during which various proteins, carbohydrates, and lipids are synthesized and stored inside the developing oocytes. Through analyzing protein changes in the midgut diverticula, hemolymph, and ovaries of females throughout the different vitellogenic stages of the spider Polybetes pythagoricus, we determined the origin of the different proteins involved in the formation of lipovitellins (LVs) along with the existence of a linkage between the hemocyanin and this vital process. An increase in the total protein content of the midgut diverticula, hemolymph, and ovary occurred throughout vitellogenesis followed by a decrease in those levels after laying. The presence of hemocyanin in egg and in LV2, as well as its accumulation in the ovary throughout the vitellogenesis process, was determined. Considering that all biologic processes depend on the correct structure and function of proteins, this study establishes, for the first time for the Order Araneae, the coexistence of three different origins of vitellogenesis-related proteins: one predominantly ovarian involving peptides of 120, 75, 46, and 30 kDa; another extraovarian one originated from the midgut diverticula and represented by a 170 kDa peptide, and a third hemolymphatic one, represented by the 67 kDa peptide.

The Proteome of Cataract Markers: Focus on Crystallins

Cataract is a major cause of blindness worldwide. It is characterized by lens opacification and is accompanied by extensive posttranslational modifications (PTMs) in various proteins. PTMs play an essential role in lens opacification. Several PTMs have been described in proteins isolated from relatively old human lenses, including phosphorylation, deamidation, racemization, truncation, acetylation, and methylation. An overwhelming majority of previous cataract proteomic studies have exclusively focused on crystallin proteins, which are the most abundant proteome components of the lens. To investigate the proteome of cataract markers, this chapter focuses on the proteomic research on the functional relevance of the major PTMs in crystallins of human cataractous lenses. Elucidating the role of these modifications in cataract formation has been a challenging task because they are among the most difficult PTMs to study analytically. The proteomic status of some amides presents similar properties in normal aged and cataractous lenses, whereas some may undergo greater PTMs in cataract. Therefore, it is of great importance to review the current proteomic research on crystallins, the major protein markers in different types of cataract, to elucidate the pathogenesis of this major human-blinding condition.

A Combined NMR and SAXS Analysis of the Partially Folded Cataract-Associated V75D γD-Crystallin

A cataract is a pathological condition characterized by the clouding of the normally clear eye lens brought about by deposition of crystallin proteins in the lens fiber cells. These protein aggregates reduce visual acuity by scattering or blocking incoming light. Chemical damage to proteins of the crystallin family, accumulated over a lifetime, leads to age-related cataract, whereas inherited mutations are associated with congenital or early-onset cataract. The V75D mutant of γD-crystallin is associated with congenital cataract in mice and was previously shown to un/fold via a partially folded intermediate. Here, we structurally characterized the stable equilibrium urea unfolding intermediate of V75D at the ensemble level using solution NMR and small-angle x-ray scattering. Our data show that, in the intermediate, the C-terminal domain retains a folded conformation that is similar to the native wild-type protein, whereas the N-terminal domain is unfolded and comprises an ensemble of random conformers, without any detectable residual structural propensities.

Accelerated Stability Modeling for Peptides: a Case Study with Bacitracin

The Accelerated Stability Assessment Program (ASAP) was applied for the first time to a peptide, the antibiotic active pharmaceutical ingredient bacitracin. Bacitracin and its complex with zinc were exposed to temperature and relative humidity conditions from 50 to 80°C and from 0 to 63% for up to 21 days. High-performance liquid chromatography was used to analyze the stressed samples for both degradant formation and loss of the active (bacitracin A) and two inactive isoforms, with identities confirmed by mass spectrometry. These data were then analyzed using a humidity-corrected Arrhenius equation and isoconversion approach to create a shelf-life predicting model for typical storage conditions. Model fitting was found to be good with low residuals in both temperature and relative humidity axes for all parameters examined. The generated model's predictions for both the native and zinc complex of bacitracin for both formation of the major degradation product (F) and loss of the active isoform (A) were consistent with longer-term measured values at 30°C/53%RH and 40°C/75%RH, validating this approach for accelerating the determination of long-term stability of a peptide.

Thyroglobulin From Molecular and Cellular Biology to Clinical Endocrinology

Thyroglobulin (Tg) is a vertebrate secretory protein synthesized in the thyrocyte endoplasmic reticulum (ER), where it acquires N-linked glycosylation and conformational maturation (including formation of many disulfide bonds), leading to homodimerization. Its primary functions include iodide storage and thyroid hormonogenesis. Tg consists largely of repeating domains, and many tyrosyl residues in these domains become iodinated to form monoiodo- and diiodotyrosine, whereas only a small portion of Tg structure is dedicated to hormone formation. Interestingly, evolutionary ancestors, dependent upon thyroid hormone for development, synthesize thyroid hormones without the complete Tg protein architecture. Nevertheless, in all vertebrates, Tg follows a strict pattern of region I, II-III, and the cholinesterase-like (ChEL) domain. In vertebrates, Tg first undergoes intracellular transport through the secretory pathway, which requires the assistance of thyrocyte ER chaperones and oxidoreductases, as well as coordination of distinct regions of Tg, to achieve a native conformation. Curiously, regions II-III and ChEL behave as fully independent folding units that could function as successful secretory proteins by themselves. However, the large Tg region I (bearing the primary T4-forming site) is incompetent by itself for intracellular transport, requiring the downstream regions II-III and ChEL to complete its folding. A combination of nonsense mutations, frameshift mutations, splice site mutations, and missense mutations in Tg occurs spontaneously to cause congenital hypothyroidism and thyroidal ER stress. These Tg mutants are unable to achieve a native conformation within the ER, interfering with the efficiency of Tg maturation and export to the thyroid follicle lumen for iodide storage and hormonogenesis.

Ca2+ and βγ-crystallins: An affair that did not last?

BACKGROUND

During the last three decades, lens β- and γ-crystallins have found a huge number of kin from numerous taxonomical sources. Most of these proteins from invertebrates and microbes have been demonstrated or predicted to bind Ca2+ involving a distinct double-clamp motif, which is largely degenerated in lens homologues.

SCOPE OF REVIEW

The various aspects of transformation of βγ-crystallins from a quintessential Ca2+-binding protein into a primarily structural molecule have been reviewed.

MAJOR CONCLUSIONS

In lens members of βγ-crystallins, the residues involved in Ca2+ binding have diverged considerably from the classical consensus with consequent reduction in their Ca2+-binding properties. This evolutionary change is congenial to their new role as robust constituents of lens. The exact functions of the residual affinity for Ca2+ are yet to be established.

GENERAL SIGNIFICANCE

This review highlights the significance of reduction in Ca2+-binding ability of the βγ-crystallins for lens physiology and why this residual affinity may be functionally important. This article is part of a Special Issue entitled Crystallin Biochemistry in Health and Disease.

Thermodynamic Protein Destabilization by GFP Tagging: A Case of Interdomain Allostery

The Engrailed Homeodomain (EnHD) transcription factor of Drosophila melanogaster was fused to the enhanced green fluorescent protein (eGFP) either at its C- or N-terminus via three- or ten-residue flexible linkers. Here, we show that EnHD undergoes destabilization upon fusing it to eGFP regardless of the linker length used and whether the tethering is to its N- or C-terminus. The destabilization is reflected in melting points that are lower by up to 9°C. Thermodynamic analysis and coarse-grained molecular dynamic simulations indicate that this destabilization is due to eGFP-promoted entropic stabilization of the denatured state ensemble of EnHD. Our results provide, therefore, an example for destabilizing interdomain allostery. They are also important given the widespread use of eGFP tagging in cell biology, as they indicate that such tagging can cause unintended protein destabilization and concomitant effects.

The fundamental tradeoff in genomes and proteomes of prokaryotes established by the genetic code, codon entropy, and physics of nucleic acids and proteins

Background

Mutations in nucleotide sequences provide a foundation for genetic variability, and selection is the driving force of the evolution and molecular adaptation. Despite considerable progress in the understanding of selective forces and their compositional determinants, the very nature of underlying mutational biases remains unclear.

Results

We explore here a fundamental tradeoff, which analytically describes mutual adjustment of the nucleotide and amino acid compositions and its possible effect on the mutational biases. The tradeoff is determined by the interplay between the genetic code, optimization of the codon entropy, and demands on the structure and stability of nucleic acids and proteins.

Conclusion

The tradeoff is the unifying property of all prokaryotes regardless of the differences in their phylogenies, life styles, and extreme environments. It underlies mutational biases characteristic for genomes with different nucleotide and amino acid compositions, providing foundation for evolution and adaptation.

Reviewers

This article was reviewed by Eugene Koonin, Michael Gromiha, and Alexander Schleiffer.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-014-0029-2) contains supplementary material, which is available to authorized users.

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.

Electrostatic origin of in vitro aggregation of human γ-crystallin

The proteins α-, β-, and γ-crystallins are the major components of the lens in the human eye. Using dynamic light scattering method, we have performed in vitro investigations of protein-protein interactions in dilute solutions of human γ-crystallin and α-crystallin. We find that γ-crystallin spontaneously aggregates into finite-sized clusters in phosphate buffer solutions. There are two distinct populations of unaggregated and aggregated γ-crystallins in these solutions. On the other hand, α-crystallin molecules are not aggregated into large clusters in solutions of α-crystallin alone. When α-crystallin and γ-crystallin are mixed in phosphate buffer solutions, we demonstrate that the clusters of γ-crystallin are prevented. By further investigating the roles of temperature, protein concentration, pH, salt concentration, and a reducing agent, we show that the aggregation of γ-crystallin under our in vitro conditions arises from non-covalent electrostatic interactions. In addition, we show that aggregation of γ-crystallin occurs under the dilute in vitro conditions even in the absence of oxidizing agents that can induce disulfide cross-links, long considered to be responsible for human cataracts. Aggregation of γ-crystallin when maintained under reducing conditions suggests that oxidation does not contribute to the aggregation in dilute solutions.

Molecular mechanisms of adaptation emerging from the physics and evolution of nucleic acids and proteins

DNA, RNA and proteins are major biological macromolecules that coevolve and adapt to environments as components of one highly interconnected system. We explore here sequence/structure determinants of mechanisms of adaptation of these molecules, links between them, and results of their mutual evolution. We complemented statistical analysis of genomic and proteomic sequences with folding simulations of RNA molecules, unraveling causal relations between compositional and sequence biases reflecting molecular adaptation on DNA, RNA and protein levels. We found many compositional peculiarities related to environmental adaptation and the life style. Specifically, thermal adaptation of protein-coding sequences in Archaea is characterized by a stronger codon bias than in Bacteria. Guanine and cytosine load in the third codon position is important for supporting the aerobic life style, and it is highly pronounced in Bacteria. The third codon position also provides a tradeoff between arginine and lysine, which are favorable for thermal adaptation and aerobicity, respectively. Dinucleotide composition provides stability of nucleic acids via strong base-stacking in ApG dinucleotides. In relation to coevolution of nucleic acids and proteins, thermostability-related demands on the amino acid composition affect the nucleotide content in the second codon position in Archaea.

Investigations on the role of π-π interactions and π-π networks in eNOS and nNOS proteins

π-π Interactions play an important role in the stability of protein structures. In the present study, we have analyzed the influence of π-π interactions in eNOS and nNOS proteins. The contribution of these π-π interacting residues in sequential separation, secondary structure involvement, solvent accessibility and stabilization centers has been evaluated. π-π interactions stabilize the core regions within eNOS and nNOS proteins. π-π interacting residues are evolutionary conserved. There is a significant number of π-π interactions in spite of the lesser natural occurrences of π-residues in eNOS and nNOS proteins. In addition to π-π interactions, π residues also form π-π networks in both eNOS and nNOS proteins which might play an important role in the structural stability of these protein structures.

Potential of human γD-crystallin for hair damage repair: insights into the mechanical properties and biocompatibility

OBJECTIVES

The objective of this work was to develop a new strategy to physically 'repair' chemically damaged hair. Hence the human eye γD-crystallin, a protein from the superfamily characterized structurally by the Greek key motif, was studied. The human γD-crystallin was chosen based on the ability of proteins belonging to this superfamily to be involved in the coating of specific structures. Two crystallins were used on the study, the wild type (Protein Data Bank ID: 1HK0) and the mutant protein. The mutant form was intended to induce a strong and quick protein polymerization as well to have new possible points of anchorage to hair.

METHODS

The ability of both crystallins to bind to damaged hair and even penetrate into its cortex was checked by fluorescence microscopy, confocal microscopy and scanning electron microscopy. Furthermore the reinforcement of hair mechanical resistance, the potential cytotoxic/inflammatory effect of crystallins were studied in order to have a fully comprehension about the protein based formulation.

RESULTS

Although the chemical over-bleaching treatment induced a decrease of 20% on the resistance of the hair, the crystallins which bind and penetrate the hair fibre were able to recover and even to improve its mechanical properties when compared to the virgin hair. Moreover none of the crystallins displayed a toxic effect in fibroblasts for all the range of tested concentrations upon 72 h of exposure. The active aggregation process of mutant crystallin induced an inflammatory response in fibroblasts in the first 24 h of contact, measured by the amount of released pro-inflammatory cytokine IL-6 to the medium. In contrast contact with wild type crystallin did not lead to significant inflammation.

CONCLUSION

Outcome from protein formulation characterization supports the hypothesis that the γD-crystallin it is able to recover and improve the mechanical properties of chemical damaged hair. Therefore it can be considered as a very promising strengthening agent for the development of new restorative hair care products.

Structure and dynamics of the fish eye lens protein, γM7-crystallin

The vertebrate eye lens contains high concentrations of crystallins. The dense lenses of fish are particularly abundant in a class called γM-crystallin whose members are characterized by an unusually high methionine content and partial loss of the four tryptophan residues conserved in all γ-crystallins from mammals which are proposed to contribute to protection from UV-damage. Here, we present the structure and dynamics of γM7-crystallin from zebrafish (Danio rerio). The solution structure shares the typical two-domain, four-Greek-key motif arrangement of other γ-crystallins, with the major difference noted in the final loop of the N-terminal domain, spanning residues 65-72. This is likely due to the absence of the conserved tryptophans. Many of the methionine residues are exposed on the surface but are mostly well-ordered and frequently have contacts with aromatic side chains. This may contribute to the specialized surface properties of these proteins that exist under high molecular crowding in the fish lens. NMR relaxation data show increased backbone conformational motions in the loop regions of γM7 compared to those of mouse γS-crystallin and show that fast internal motion of the interdomain linker in γ-crystallins correlates with linker length. Unfolding studies monitored by tryptophan fluorescence confirm results from mutant mouse γS-crystallin and show that unfolding of a βγ-crystallin domain likely starts from unfolding of the variable loop containing the more fluorescently quenched tryptophan residue, resulting in a native-like unfolding intermediate.

Dissimilarity in the folding of human cytosolic creatine kinase isoenzymes

Creatine kinase (CK, EC 2.7.3.2) plays a key role in the energy homeostasis of excitable cells. The cytosolic human CK isoenzymes exist as homodimers (HMCK and HBCK) or a heterodimer (MBCK) formed by the muscle CK subunit (M) and/or brain CK subunit (B) with highly conserved three-dimensional structures composed of a small N-terminal domain (NTD) and a large C-terminal domain (CTD). The isoforms of CK provide a novel system to investigate the sequence/structural determinants of multimeric/multidomain protein folding. In this research, the role of NTD and CTD as well as the domain interactions in CK folding was investigated by comparing the equilibrium and kinetic folding parameters of HMCK, HBCK, MBCK and two domain-swapped chimeric forms (BnMc and MnBc). Spectroscopic results indicated that the five proteins had distinct structural features depending on the domain organizations. MBCK BnMc had the smallest CD signals and the lowest stability against guanidine chloride-induced denaturation. During the biphasic kinetic refolding, three proteins (HMCK, BnMc and MnBc), which contained either the NTD or CTD of the M subunit and similar microenvironments of the Trp fluorophores, refolded about 10-fold faster than HBCK for both the fast and slow phase. The fast folding of these three proteins led to an accumulation of the aggregation-prone intermediate and slowed down the reactivation rate thereby during the kinetic refolding. Our results suggested that the intra- and inter-subunit domain interactions modified the behavior of kinetic refolding. The alternation of domain interactions based on isoenzymes also provides a valuable strategy to improve the properties of multidomain enzymes in biotechnology.

The molecular refractive function of lens γ-Crystallins

γ-Crystallins constitute the major protein component in the nucleus of the vertebrate eye lens. Present at very high concentrations, they exhibit extreme solubility and thermodynamic stability to prevent scattering of light and formation of cataracts. However, functions beyond this structural role have remained mostly unclear. Here, we calculate molecular refractive index increments of crystallins. We show that all lens γ-crystallins have evolved a significantly elevated molecular refractive index increment, which is far above those of most proteins, including nonlens members of the βγ-crystallin family from different species. The same trait has evolved in parallel in crystallins of different phyla, including S-crystallins of cephalopods. A high refractive index increment can lower the crystallin concentration required to achieve a suitable refractive power of the lens and thereby reduce their propensity to aggregate and form cataracts. To produce a significant increase in the refractive index increment, a substantial global shift in amino acid composition is required, which can naturally explain the highly unusual amino acid composition of γ-crystallins and their functional homologues. This function provides a new perspective for interpreting their molecular structure.

Differential scanning calorimetry: An invaluable tool for a detailed thermodynamic characterization of macromolecules and their interactions

Differential Scanning Calorimetry (DSC) is a highly sensitive technique to study the thermotropic properties of many different biological macromolecules and extracts. Since its early development, DSC has been applied to the pharmaceutical field with excipient studies and DNA drugs. In recent times, more attention has been applied to lipid-based drug delivery systems and drug interactions with biomimetic membranes. Highly reproducible phase transitions have been used to determine values, such as, the type of binding interaction, purity, stability, and release from a drug delivery mechanism. This review focuses on the use of DSC for biochemical and pharmaceutical applications.

Impact of inter-subunit interactions on the dimeric arginine kinase activity and structural stability

Arginine kinase (AK) is a key enzyme for cellular energy metabolism, catalyzing the reversible phosphoryl transfer from phosphoarginine to ADP in invertebrates. In this study, the inter-subunit hydrogen bonds between the Q53 and D200 and between D57 and D200 were disrupted to explore their roles in the activity and structural stability of Stichopus japonicus (S. japonicus) AK. Mutating Q53 and/or D57 to alanine (A) can cause pronounced loss of activity and substrate synergism, and cause distinct conformational changes. Spectroscopic experiments indicated that mutations destroying the inter-subunit hydrogen bonds impaired the structure of dimer AK, and resulted in a partially unfolded state. The inability to fold to the functional compact state made the mutants prone to be inactivated and aggregate under environmental stresses. Restoring hydrogen bonds in Q53E and D57E mutants could rescue the loss of activity and substrate synergism, and conformational changes. All those results suggested that the inter-subunit interactions played a key role in keeping the activity, substrate synergism and structural stability of dimer AK. The result herein may provide a clue in understanding the folding and self-assembly processes of oligomeric proteins.

The interaction of ammonia and xenon with the imidazole glycerol phosphate synthase from Thermotoga maritima as detected by NMR spectroscopy

The imidazole glycerol phosphate (ImGP) synthase from the hyperthermophilic bacterium Thermotoga maritima is a 1:1 complex of the glutaminase subunit HisH and the cyclase subunit HisF. It has been proposed that ammonia generated by HisH is transported through a channel to the active site of HisF, which generates intermediates of histidine (ImGP) and de novo biosynthesis of 5-aminoimidazole-4-carboxamideribotide. Solution NMR spectroscopy of ammonium chloride-titrated samples was used to study the interaction of NH(3) with amino acids inside this channel. Although numerous residues showed (15)N chemical shift changes, most of these changes were caused by nonspecific ionic strength effects. However, several interactions appeared to be specific. Remarkably, the amino acid residue Thr 78-which is located in the central channel-shows a large chemical shift change upon titration with ammonium chloride. This result and the reduced catalytic activity of the Thr78Met mutant indicate a special role of this residue in ammonia channeling. To detect and further characterize internal cavities in HisF, which might for example contribute to ammonia channeling, the interaction of HisF with the noble gas xenon was analyzed by solution NMR spectroscopy using (1)H-(15)N HSQC experiments. The results indicate that HisF contains three distinct internal cavities, which could be identified by xenon-induced chemical shift changes of the neighboring amino acid residues. Two of these cavities are located at the active site at opposite ends of the substrate N'-[(5'-phosphoribulosyl)formimino]-5-aminoimidazole-4-carboxamide-ribonucleotide (PRFAR) binding groove. The third cavity is located in the interior of the central β-barrel of HisF and overlaps with the putative ammonia transport channel.

Structural characterization and unfolding mechanism of human 4F2hc ectodomain

4F2hc (CD98hc) is a multifunctional type II membrane glycoprotein involved in several functions as amino acid transport, cell fusion, β1-integrin-signaling and transformation. 4F2hc ectodomain has been crystallized and its three-dimensional structure determined. We have carried out a spectroscopical/structural characterization of the recombinant ectodomain in order to obtain information on its dynamic structure in solution and on its ability to form homodimers by itself in the absence of the transmembrane helix and of the potential interactions with the plasma membrane. Analytical ultracentrifugation and crosslinking experiments showed that the ectodomain is monomeric in solution. The secondary structure determined by far-UV circular dichroism (CD) spectroscopy (around 30% α-helix and 20% β-sheets, 12% antiparallel and 8% parallel) reveals a compact and thermally stable structure with a high melting temperature (57-59°C). Tryptophan residues are mainly buried and immobilized in the hydrophobic core of the protein as suggested by near-UV CD spectrum, the position of the Trp maximum fluorescence emission (323nm) and from the acrylamide quenching constant (2.6M(-1)). Urea unfolding equilibrium has been studied by far-UV CD and fluorescence spectroscopy to gain information on the folding/unfolding process of the ectodomain. The analyses suggest the existence of two intermediate states as reported for other TIM barrel-containing proteins rather than an independent unfolding of each domain [A, (βα)(8) barrel; C, antiparallel β(8) sandwich]. Folding seems to be directed by the initial formation of hydrophobic clusters within the first strands of the β-barrel of domain A followed by additional hydrophobic interactions in domain C.

Diversity of function-related conformational changes in proteins: coordinate uncertainty, fragment rigidity, and stability

It was found that the variety of function-related conformational changes ("movements") in proteins is beyond the earlier simple classifications. Here we offer biochemists a more comprehensive, transparent, and easy-to-use approach allowing a detailed and accurate interpretation of such conformational changes. It makes possible a more multifaceted characterization of protein flexibility via identification of rigidly and nonrigidly repositioned fragments, stable and nonstable fragments, and domain and nondomain repositioning. "Coordinate uncertainty thresholds" derived from computed differences between independently determined coordinates of the same molecules are used as the criteria for conformational identity. "Identical" rigid substructures are localized in the distance difference matrices (DDMs). A sequence of simple transformations determines whether a structural change occurs by rigid-body movements of fragments or largely through non-rigid-body deformations. We estimate the stability of protein fragments and compare stable and rigidly moving fragments. The motions computed with the coarse-grained elastic networks are also compared to those of their DDM analogues. We study and suggest a classification for 17 structural pairs, differing in their functional states. For five of the 17 proteins, conformational change cannot be accomplished by rigid-body transformations and requires significant non-rigid-body deformations. Stable fragments rarely coincide with rigidly moving fragments and often disagree with the CATH identifications of domains. Almost all monomeric apo chains, containing stable fragments and/or domains, indicate instability of the entire molecule, suggesting the importance of fragments and domains motions prior to stabilization by substrate binding or crystallization. Notably, kinases exhibit the greatest extent of nonrigidity among the proteins investigated.

Buffer-dependent fragmentation of a humanized full-length monoclonal antibody

During storage stability studies of a monoclonal antibody (mAb) it was determined that the primary route of degradation involved fragmentation into lower molecular weight species. The fragmentation was characterized with size-exclusion high performance liquid chromatography (SE-HPLC), SDS-PAGE, and matrix-assisted laser desorption/ionization time of flight (MALDI-TOF) mass spectrometry. Fragmentation proceeded via hydrolysis, likely catalyzed by trace metal ions, of a peptide bond in the hinge region of the mAb's heavy chain, which produced two prominent low molecular weight species during storage: a single, free Fab fragment and a Fab + Fc fragment. The fragmentation is observed in phosphate-buffered solutions at two ionic strengths but not in histidine-buffered solutions at identical ionic strengths. Chaotrope-induced and thermally induced unfolding studies of the mAb indicated differences in the unfolding pathways between the two buffer solutions. The folding intermediate observed during chaotrope-induced unfolding was further characterized by intrinsic fluorescence quenching, which suggested that a small portion of the molecule is resistant to chaotrope-induced unfolding in histidine buffer systems. The thermally induced unfolding indicates a reduction in cooperativity of the unfolding process in the presence of histidine relative to phosphate. A relationship between the histidine-induced effects on unfolding pathway and the relative resistance to fragmentation is suggested.

Sequential events in the irreversible thermal denaturation of human brain-type creatine kinase by spectroscopic methods

The non-cooperative or sequential events which occur during protein thermal denaturation are closely correlated with protein folding, stability, and physiological functions. In this research, the sequential events of human brain-type creatine kinase (hBBCK) thermal denaturation were studied by differential scanning calorimetry (DSC), CD, and intrinsic fluorescence spectroscopy. DSC experiments revealed that the thermal denaturation of hBBCK was calorimetrically irreversible. The existence of several endothermic peaks suggested that the denaturation involved stepwise conformational changes, which were further verified by the discrepancy in the transition curves obtained from various spectroscopic probes. During heating, the disruption of the active site structure occurred prior to the secondary and tertiary structural changes. The thermal unfolding and aggregation of hBBCK was found to occur through sequential events. This is quite different from that of muscle-type CK (MMCK). The results herein suggest that BBCK and MMCK undergo quite dissimilar thermal unfolding pathways, although they are highly conserved in the primary and tertiary structures. A minor difference in structure might endow the isoenzymes dissimilar local stabilities in structure, which further contribute to isoenzyme-specific thermal stabilities.

A single destabilizing mutation (F9S) promotes concerted unfolding of an entire globular domain in gammaS-crystallin

Conformational change and aggregation of native proteins are associated with many serious age-related and neurological diseases. gammaS-Crystallin is a highly stable, abundant structural component of vertebrate eye lens. A single F9S mutation in the N-terminal domain of mouse gammaS-crystallin causes the severe Opj cataract, with disruption of cellular organization and appearance of fibrillar structures in the lens. Although the mutant protein has a near-native fold at room temperature, significant increases in hydrogen/deuterium exchange rates were observed by NMR for all the well-protected beta-sheet core residues throughout the entire N-terminal domain of the mutant protein, resulting in up to a 3.5-kcal/mol reduction in the free energy of the folding/unfolding equilibrium. No difference was detected for the C-terminal domain. At a higher temperature, this effect further increases to allow for a much more uniform exchange rate among the N-terminal core residues and those of the least well-structured surface loops. This suggests a concerted unfolding intermediate of the N-terminal domain, while the C-terminal domain stays intact. Increasing concentrations of guanidinium chloride produced two transitions for the Opj mutant, with an unfolding intermediate at approximately 1 M guanidinium chloride. The consequence of this partial unfolding, whether by elevated temperature or by denaturant, is the formation of thioflavin T staining aggregates, which demonstrated fibril-like morphology by atomic force microscopy. Seeding with the already unfolded protein enhanced the formation of fibrils. The Opj mutant protein provides a model for stress-related unfolding of an essentially normally folded protein and production of aggregates with some of the characteristics of amyloid fibrils.

Mapping the stability clusters in bovine pancreatic ribonuclease A

In the present work, we have thermodynamically characterized the thermally induced unfolding of 20 variants of bovine pancreatic ribonuclease A (RNase A) to experimentally describe the residues and the regions that are critical for the stability of the enzyme. The achieved results, complemented with previous studies by our group, allowed us to define the significance of the two hydrophobic nuclei present in the RNase A structure, as well as the contribution of the participating residues within each nucleus, to the global enzyme stability. We propose a structural model for the major and the minor hydrophobic nuclei of RNase A. The major nucleus is composite and located in the cavity delimited by alpha-helices 1 and 3, and the beta-sheet that is formed by strands 2, 3, 5, and 6. It consists of a central tight packed part constituted by residues Phe8, Met13, Val54, Val57, Ile106, Val108, and Phe120. This central part is surrounded by a layer formed by residues Val63, Tyr73, Met79, Ile107, Val116, and Val118. The minor nucleus, although less complex, is also constituted by a tight packing that involves the side chains of residues Tyr25, Met29, Met30, Leu35, Phe46, and Tyr97, which fill the cavity that originates the beta-sheet formed by beta-strands 1, 4, and 5 together with alpha-helix2.

Residue-specific analysis of frustration in the folding landscape of repeat beta/alpha protein apoflavodoxin

Flavodoxin adopts the common repeat beta/alpha topology and folds in a complex kinetic reaction with intermediates. To better understand this reaction, we analyzed a set of Desulfovibrio desulfuricans apoflavodoxin variants with point mutations in most secondary structure elements by in vitro and in silico methods. By equilibrium unfolding experiments, we first revealed how different secondary structure elements contribute to overall protein resistance to heat and urea. Next, using stopped-flow mixing coupled with far-UV circular dichroism, we probed how individual residues affect the amount of structure formed in the experimentally detected burst-phase intermediate. Together with in silico folding route analysis of the same point-mutated variants and computation of growth in nucleation size during early folding, computer simulations suggested the presence of two competing folding nuclei at opposite sides of the central beta-strand 3 (i.e., at beta-strands 1 and 4), which cause early topological frustration (i.e., misfolding) in the folding landscape. Particularly, the extent of heterogeneity in folding nuclei growth correlates with the in vitro burst-phase circular dichroism amplitude. In addition, phi-value analysis (in vitro and in silico) of the overall folding barrier to apoflavodoxin's native state revealed that native-like interactions in most of the beta-strands must form in transition state. Our study reveals that an imbalanced competition between the two sides of apoflavodoxin's central beta-sheet directs initial misfolding, while proper alignment on both sides of beta-strand 3 is necessary for productive folding.