Entropic stabilization of proteins and its proteomic consequences

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.

Titin domains with reduced core hydrophobicity cause dilated cardiomyopathy

The underlying genetic defect in most cases of dilated cardiomyopathy (DCM), a common inherited heart disease, remains unknown. Intriguingly, many patients carry single missense variants of uncertain pathogenicity targeting the giant protein titin, a fundamental sarcomere component. To explore the deleterious potential of these variants, we first solved the wild-type and mutant crystal structures of I21, the titin domain targeted by pathogenic variant p.C3575S. Although both structures are remarkably similar, the reduced hydrophobicity of deeply buried position 3575 strongly destabilizes the mutant domain, a scenario supported by molecular dynamics simulations and by biochemical assays that show no disulfide involving C3575. Prompted by these observations, we have found that thousands of similar hydrophobicity-reducing variants associate specifically with DCM. Hence, our results imply that titin domain destabilization causes DCM, a conceptual framework that not only informs pathogenicity assessment of gene variants but also points to therapeutic strategies counterbalancing protein destabilization.

Evolution of Organic Solvent-Resistant DNA Polymerases

The introduction of thermostable polymerases revolutionized the polymerase chain reaction (PCR) and biotechnology. However, many GC-rich genes cannot be PCR-amplified with high efficiency in water, irrespective of temperature. Although polar organic cosolvents can enhance nucleic acid polymerization and amplification by destabilizing duplex DNA and secondary structures, nature has not selected for the evolution of solvent-tolerant polymerase enzymes. Here, we used ultrahigh-throughput droplet-based selection and deep sequencing along with computational free-energy and binding affinity calculations to evolve Taq polymerase to generate enzymes that are both stable and highly active in the presence of organic cosolvents, resulting in up to 10% solvent resistance and over 100-fold increase in stability at 97.5 °C in the presence of 1,4-butanediol, as well as tolerance to up to 10 times higher concentrations of the potent cosolvents sulfolane and 2-pyrrolidone. Using these polymerases, we successfully amplified a broad spectrum of GC-rich templates containing regions with over 90% GC content, including templates recalcitrant to amplification with existing polymerases, even in the presence of cosolvents. We also demonstrated dramatically reduced GC bias in the amplification of genes with widely varying GC content in quantitative polymerase chain reaction (qPCR). By expanding the scope of solvent systems compatible with nucleic acid polymerization, these organic solvent-resistant polymerases enable a dramatic reduction of sequence bias not achievable through thermal resistance alone, with significant implications for a wide range of applications including sequencing and synthetic biology in mixed aqueous-organic media.

A proteome-scale analysis of vertebrate protein amino acid occurrence: Thermoadaptation shows a correlation with protein solvation but less so with dynamics

Despite differences in behaviors and living conditions, vertebrate organisms share the great majority of proteins, often with subtle differences in amino acid sequence. Here, we present a simple way to analyze the difference in amino acid occurrence by comparing highly homologous proteins on a subproteome level between several vertebrate model organisms. Specifically, we use this method to identify a pattern of amino acid conservation as well as a shift in amino acid occurrence between homeotherms (warm-blooded species) and poikilotherms (cold-blooded species). Importantly, this general analysis and a specific example further establish a broad correlation, if not likely connection between the thermal adaptation of protein sequences and two of their physical features: on average a change in their protein dynamics and, even more strongly, in their solvation. For poikilotherms, such as frog and fish, the lower body temperature is expected to increase the protein-protein interaction due to a decrease in protein internal dynamics. In order to counteract the tendency for enhanced binding caused by low temperatures, poikilotherms enhance the solvation of their proteins by favoring polar amino acids. This feature appears to dominate over possible changes in dynamics for some proteins. The results suggest that a general trend for amino acid choice is part of the mechanism for thermoadaptation of vertebrate organisms at the molecular level.

Thermostabilization mechanisms in thermophilic versus mesophilic three-helix bundle proteins

The engineered three‐helix bundle, UVF, is thermostabilized entropically due to heightened, native‐state dynamics. However, it is unclear whether this thermostabilization strategy is observed in natural proteins from thermophiles. We performed all‐atom, explicit solvent molecular dynamics simulations of two three‐helix bundles from thermophilic H. butylicus (2lvsN and 2lvsC) and compared their dynamics to a mesophilic three‐helix bundle, the Engrailed homeodomain (EnHD). Like UVF, 2lvsC had heightened native dynamics, which it maintained without unfolding at 100°C. Shortening and rigidification of loops in 2lvsN and 2lvsC and increased surface hydrogen bonds in 2lvsN were observed, as is common in thermophilic proteins. A buried disulfide and salt bridge in 2lvsN and 2lvsC, respectively, provided some stabilization, and addition of a homologous disulfide bond in EnHD slowed unfolding. The transferability and commonality of stabilization strategies among members of the three‐helix bundle fold suggest that these strategies may be general and deployable in designing thermostable proteins.

Determination of the Amino Acid Recruitment Order in Early Life by Genome-Wide Analysis of Amino Acid Usage Bias

The mechanisms shaping the amino acids recruitment pattern into the proteins in the early life history presently remains a huge mystery. In this study, we conducted genome-wide analyses of amino acids usage and genetic codons structure in 7270 species across three domains of life. The carried-out analyses evidenced ubiquitous usage bias of amino acids that were likely independent from codon usage bias. Taking advantage of codon usage bias, we performed pseudotime analysis to re-determine the chronological order of the species emergence, which inspired a new species relationship by tracing the imprint of codon usage evolution. Furthermore, the multidimensional data integration showed that the amino acids A, D, E, G, L, P, R, S, T and V might be the first recruited into the last universal common ancestry (LUCA) proteins. The data analysis also indicated that the remaining amino acids most probably were gradually incorporated into proteogenesis process in the course of two long-timescale parallel evolutionary routes: I→F→Y→C→M→W and K→N→Q→H. This study provides new insight into the origin of life, particularly in terms of the basic protein composition of early life. Our work provides crucial information that will help in a further understanding of protein structure and function in relation to their evolutionary history.

Double Mutant of Chymotrypsin Inhibitor 2 Stabilized through Increased Conformational Entropy

The conformational heterogeneity of a folded protein can affect not only its function but also stability and folding. We recently discovered and characterized a stabilized double mutant (L49I/I57V) of the protein CI2 and showed that state-of-the-art prediction methods could not predict the increased stability relative to the wild-type protein. Here, we have examined whether changed native-state dynamics, and resulting entropy changes, can explain the stability changes in the double mutant protein, as well as the two single mutant forms. We have combined NMR relaxation measurements of the ps-ns dynamics of amide groups in the backbone and the methyl groups in the side chains with molecular dynamics simulations to quantify the native-state dynamics. The NMR experiments reveal that the mutations have different effects on the conformational flexibility of CI2: a reduction in conformational dynamics (and entropy estimated from this) of the native state of the L49I variant correlates with its decreased stability, while increased dynamics of the I57V and L49I/I57V variants correlates with their increased stability. These findings suggest that explicitly accounting for changes in native-state entropy might be needed to improve the predictions of the effect of mutations on protein stability.

Order-disorder transitions of cytoplasmic N-termini in the mechanisms of P-type ATPases

Membrane protein structure and function are modulated via interactions with their lipid environment. This is particularly true for integral membrane pumps, the P-type ATPases. These ATPases play vital roles in cell physiology, where they are associated with the transport of cations and lipids, thereby generating and maintaining crucial (electro-)chemical potential gradients across the membrane. Several pumps (Na⁺, K⁺-ATPase, H⁺, K⁺-ATPase and the plasma membrane Ca²⁺-ATPase) which are located in the asymmetric animal plasma membrane have been found to possess polybasic (lysine-rich) domains on their cytoplasmic surfaces, which are thought to act as phosphatidylserine (PS) binding domains. In contrast, the sarcoplasmic reticulum Ca²⁺-ATPase, located within an intracellular organelle membrane, does not possess such a domain. Here we focus on the lysine-rich N-termini of the plasma-membrane-bound Na⁺, K⁺- and H⁺, K⁺-ATPases. Synthetic peptides corresponding to the N-termini of these proteins were found, via quartz crystal microbalance and circular dichroism measurements, to interact via an electrostatic interaction with PS-containing membranes, thereby undergoing an increase in helical or other secondary structure content. As well as influencing ion pumping activity, it is proposed that this interaction could provide a mechanism for sensing the lipid asymmetry of the plasma membrane, which changes drastically when a cell undergoes apoptosis, i.e. programmed cell death. Thus, polybasic regions of plasma membrane-bound ion pumps could potentially perform the function of a "death sensor", signalling to a cell to reduce pumping activity and save energy.

The Molecular Basis for Life in Extreme Environments

Sampling and genomic efforts over the past decade have revealed an enormous quantity and diversity of life in Earth's extreme environments. This new knowledge of life on Earth poses the challenge of understandingits molecular basis in such inhospitable conditions, given that such conditions lead to loss of structure and of function in biomolecules from mesophiles. In this review, we discuss the physicochemical properties of extreme environments. We present the state of recent progress in extreme environmental genomics. We then present an overview of our current understanding of the biomolecular adaptation to extreme conditions. As our current and future understanding of biomolecular structure-function relationships in extremophiles requires methodologies adapted to extremes of pressure, temperature, and chemical composition, advances in instrumentation for probing biophysical properties under extreme conditions are presented. Finally, we briefly discuss possible future directions in extreme biophysics.

Entropy of Proteins Using Multiscale Cell Correlation

A new multiscale method is presented to calculate the entropy of proteins from molecular dynamics simulations. Termed Multiscale Cell Correlation (MCC), the method decomposes the protein into sets of rigid-body units based on their covalent-bond connectivity at three levels of hierarchy: molecule, residue, and united atom. It evaluates the vibrational and topographical entropy from forces, torques, and dihedrals at each level, taking into account correlations between sets of constituent units that together make up a larger unit at the coarser length scale. MCC gives entropies in close agreement with normal-mode analysis and smaller than those using quasiharmonic analysis as well as providing much faster convergence. Moreover, MCC provides an insightful decomposition of entropy at each length scale and for each type of amino acid according to their solvent exposure and whether they are terminal residues. While the residue entropy depends weakly on solvent exposure, there is greater variation in entropy components for larger, more polar amino acids, which have increased conformational entropy but reduced vibrational entropy with greater solvent exposure.

Nearest-neighbor amino acids of specificity-determining residues influence the activity of engineered Cre-type recombinases

The tyrosine-type site-specific DNA recombinase Cre recombines its target site, loxP, with high activity and specificity without cross-recombining the target sites of highly related recombinases. Understanding how Cre achieves this precision is key to be able to rationally engineer site-specific recombinases (SSRs) for genome editing applications. Previous work has revealed key residues for target site selectivity in the Cre/loxP and the related Dre/rox recombinase systems. However, enzymes in which these residues were changed to the respective counterpart only showed weak activity on the foreign target site. Here, we use molecular modeling and dynamics simulation techniques to comprehensively explore the mechanisms by which these residues determine target recognition in the context of their flanking regions in the protein-DNA interface, and we establish a structure-based rationale for the design of improved recombination activities. Our theoretical models reveal that nearest-neighbors to the specificity-determining residues are important players for enhancing SSR activity on the foreign target site. Based on the established rationale, we design new Cre variants with improved rox recombination activities, which we validate experimentally. Our work provides new insights into the target recognition mechanisms of Cre-like recombinases and represents an important step towards the rational design of SSRs for applied genome engineering.

Effects of Single Mutations on Protein Stability Are Gaussian Distributed

The distribution of protein stability effects is known to be well approximated by a Gaussian distribution from previous empirical fits. Starting from first-principles statistical mechanics, we more rigorously motivate this empirical observation by deriving per-residue-position protein stability effects to be Gaussian. Our derivation requires the number of amino acids to be large, which is satisfied by the standard set of 20 amino acids found in nature. No assumption is needed on the number of residues in close proximity in space, in contrast to previous applications of the central limit theorem to protein energetics. We support our derivation results with computational and experimental data on mutant protein stabilities across all types of protein residues.

Protein adaptation to high hydrostatic pressure: Computational analysis of the structural proteome

Hydrostatic pressure has a vital role in the biological adaptation of the piezophiles, organisms that live under high hydrostatic pressure. However, the mechanisms by which piezophiles are able to adapt their proteins to high hydrostatic pressure is not well understood. One proposed hypothesis is that the volume changes of unfolding (ΔVTot ) for proteins from piezophiles is distinct from those of nonpiezophilic organisms. Since ΔVTot defines pressure dependence of stability, we performed a comprehensive computational analysis of this property for proteins from piezophilic and nonpiezophilic organisms. In addition, we experimentally measured the ΔVTot of acylphosphatases and thioredoxins belonging to piezophilic and nonpiezophilic organisms. Based on this analysis we concluded that there is no difference in ΔVTot for proteins from piezophilic and nonpiezophilic organisms. Finally, we put forward the hypothesis that increased concentrations of osmolytes can provide a systemic increase in pressure stability of proteins from piezophilic organisms and provide experimental thermodynamic evidence in support of this hypothesis.

Protein thermal stability

Proteins, in general, fold to a well-organized three-dimensional structure in order to function. The stability of this functional shape can be perturbed by external environmental conditions, such as temperature. Understanding the molecular factors underlying the resistance of proteins to the thermal stress has important consequences. First of all, it can aid the design of thermostable enzymes able to perform efficient catalysis in the high-temperature regime. Second, it is an essential brick of knowledge required to decipher the evolutionary pathways of life adaptation on Earth. Thanks to the development of atomistic simulations and ad hoc enhanced sampling techniques, it is now possible to investigate this problem in silico, and therefore provide support to experiments. After having described the methodological aspects, the chapter proposes an extended discussion on two problems. First, we focus on thermophilic proteins, a perfect model to address the issue of thermal stability and molecular evolution. Second, we discuss the issue of how protein thermal stability is affected by crowded in vivo-like conditions.

Finding the generalized molecular principles of protein thermal stability

Are there any generalized molecular principles of thermal adaptation? Here, integrating the concepts of structural bioinformatics, sequence analysis, and classical knot theory, we develop a robust computational framework that seeks for mechanisms of thermal adaptation by comparing orthologous mesophilic-thermophilic and mesophilic-hyperthermophilic proteins of remarkable structural and topological similarities, and still leads us to context-independent results. A comprehensive analysis of 4741 high-resolution, non-redundant X-ray crystallographic structures collected from 11 hyperthermophilic, 32 thermophilic and 53 mesophilic prokaryotes unravels at least five "nearly universal" signatures of thermal adaptation, irrespective of the enormous sequence, structure, and functional diversity of the proteins compared. A careful investigation further extracts a set of amino acid changes that can potentially enhance protein thermal stability, and remarkably, these mutations are overrepresented in protein crystallization experiments, in disorder-to-order transitions and in engineered thermostable variants of existing mesophilic proteins. These results could be helpful to find a precise, global picture of thermal adaptation.

Voltage vs. Ligand II: Structural insights of the intrinsic flexibility in cyclic nucleotide-gated channels

In the preceding article, we present a flexibility analysis of the voltage-gated ion channel (VGIC) superfamily. In this study, we describe in detail the flexibility profile of the voltage-sensor domain (VSD) and the pore domain (PD) concerning the evolution of 6TM ion channels. In particular, we highlight the role of flexibility in the emergence of CNG channels and describe a significant level of sequence similarity between the archetypical VSD and the TolQ proteins. A highly flexible S4-like segment exhibiting Lys instead Arg for these membrane proteins is reported. Sequence analysis indicates that, in addition to this S4-like segment, TolQ proteins also show similarity with specific motifs in S2 and S3 from typical V-sensors. Notably, S3 flexibility profiles from typical VSDs and S3-like in TolQ proteins are also similar. Interestingly, TolQ from early divergent prokaryotes are comparatively more flexible than those in modern counterparts or true V-sensors. Regarding the PD, we also found that 2TM K+-channels in early prokaryotes are considerably more flexible than the ones in modern microbes, and such flexibility is comparable to the one present in CNG channels. Voltage dependence is mainly exhibited in prokaryotic CNG channels whose VSD is rigid whereas the eukaryotic CNG channels are considerably more flexible and poorly V-dependent. The implication of the flexibility present in CNG channels, their sensitivity to cyclic nucleotides and the cation selectivity are discussed. Finally, we generated a structural model of the putative cyclic nucleotide-modulated ion channel, which we coined here as AqK, from the thermophilic bacteria Aquifex aeolicus, one of the earliest diverging prokaryotes known. Overall, our analysis suggests that V-sensors in CNG-like channels were essentially rigid in early prokaryotes but raises the possibility that this module was probably part of a very flexible stator protein of the bacterial flagellum motor complex.

An NMR Confirmation for Increased Folded State Entropy Following Loop Truncation

Previous studies conducted on flexible loop regions in proteins revealed that the energetic consequences of changing loop length predominantly arise from the entropic cost of ordering a loop during folding. However, in an earlier study of human acylphosphatase (hmAcP) using experimental and computational approaches, we showed that thermodynamic stabilization upon loop truncation can be attributed mainly to the increased entropy of the folded state. Here, using 15N NMR spectroscopy, we studied the effect of loop truncation on hmAcP backbone dynamics on the picosecond-nanosecond timescale with the aim of confirming the effect of folded state entropy on protein stability. NMR-relaxation-derived N-H squared generalized order parameters reveal that loop truncation results in a significant increase in protein conformational flexibility. Comparison of these results with previously acquired all-atom molecular dynamics simulation, analyzed here in terms of squared generalized NMR order parameters, demonstrates general agreement between the two methods. The NMR study not only provides direct evidence for the enhanced conformational entropy of the folded state of hmAcP upon loop truncation but also gives a quantitative measure of the observed effects.

Functional Properties of Amino Acid Side Chains as Biomarkers of Extraterrestrial Life

The present study proposes to search our solar system (Mars, Enceladus, Europa) for patterns of organic molecules that are universally associated with biological functions and structures. The functions are primarily catalytic because life could only have originated within volume/space-constrained compartments containing chemical reactions catalyzed by certain polymers. The proposed molecular structures are specific groups in the side chains of amino acids with the highest catalytic propensities related to life on Earth, that is, those that most frequently participate as key catalytic groups in the active sites of enzymes such as imidazole, thiol, guanidinium, amide, and carboxyl. Alternatively, these or other catalytic groups can be searched for on non-amino-acid organic molecules, which can be tested for certain hydrolytic catalytic activities. The first scenario assumes that life may have originated in a similar manner as the terrestrial set of α-amino acids, while the second scenario does not set such a requirement. From the catalytic propensity perspective proposed in the first scenario, life must have invented amino acids with high catalytic propensity (His, Cys, Arg) in order to overcome, and be complemented by, the low catalytic propensity of the initially available abiogenic amino acids. The abiogenic and the metabolically invented amino acids with the lowest catalytic propensity can also serve as markers of extraterrestrial life when searching for patterns on the basis of the following functional propensities related to protein secondary/quaternary structure: (1) amino acids that are able to form α-helical intramembrane peptide domains, which can serve as primitive transporters in protocell membrane bilayers and catalysts of simple biochemical reactions; (2) amino acids that tend to accumulate in extremophile proteins of Earth and possibly extraterrestrial life. The catalytic/structural functional propensity approach offers a new perspective in the search for extraterrestrial life and could help unify previous amino acid-based approaches.

Thermophilic Proteins as Versatile Scaffolds for Protein Engineering

Literature from the past two decades has outlined the existence of a trade-off between protein stability and function. This trade-off creates a unique challenge for protein engineers who seek to introduce new functionality to proteins. These engineers must carefully balance the mutation-mediated creation and/or optimization of function with the destabilizing effect of those mutations. Subsequent research has shown that protein stability is positively correlated with "evolvability" or the ability to support mutations which bestow new functionality on the protein. Since the ultimate goal of protein engineering is to create and/or optimize a protein's function, highly stable proteins are preferred as potential scaffolds for protein engineering. This review focuses on the application potential for thermophilic proteins as scaffolds for protein engineering. The relatively high inherent thermostability of these proteins grants them a great deal of mutational robustness, making them promising scaffolds for various protein engineering applications. Comparative studies on the evolvability of thermophilic and mesophilic proteins have strongly supported the argument that thermophilic proteins are more evolvable than mesophilic proteins. These findings indicate that thermophilic proteins may represent the scaffold of choice for protein engineering in the future.

Entropic Stabilization of Cas4 Protein SSO0001 Predicted with Popcoen

Popcoen is a method for configurational entropy estimation of proteins based on machine-learning. Entropy is predicted with an artificial neural network which was trained on simulation trajectories of a large set of representative proteins. Popcoen is extremely fast compared to other approaches based on the sampling of a multitude of microstates. Consequently, Popcoen can be incorporated into a large class of protein software which currently neglects configurational entropy for performance reasons. Here, we apply Popcoen to various conformations of the Cas4 protein SSO0001 of Sulfolobus solfataricus, a protein that assembles to a decamer of known toroidal shape. We provide numerical evidence that the native state (NAT) of a SSO0001 monomer has a similar structure to the protomers of the oligomer, where NAT of the monomer is stabilized mainly entropically. Due to its large amount of configurational entropy, NAT has lower free energy than alternative conformations of very low enthalpy and solvation free-energy. Hence, SSO0001 serves as an example case where neglecting configurational entropy leads to incorrect conclusion. Our results imply that no refolding of the subunits is required during oligomerization which suggests that configurational entropy is employed by nature to largely enhance the rate of assembly.

Domain communication in Thermotoga maritima Arginine Binding Protein unraveled through protein dissection

Substrate binding proteins represent a large protein family that plays fundamental roles in selective transportation of metabolites across membrane. The function of these proteins relies on the relative motions of their two domains. Insights into domain communication in this class of proteins have been here collected using Thermotoga maritima Arginine Binding Protein (TmArgBP) as model system. TmArgBP was dissected into two domains (D1 and D2) that were exhaustively characterized using a repertoire of different experimental and computational techniques. Indeed, stability, crystalline structure, ability to recognize the arginine substrate, and dynamics of the two individual domains have been here studied. Present data demonstrate that, although in the parent protein both D1 and D2 cooperate for the arginine anchoring; only D1 is intrinsically able to bind the substrate. The implications of this finding on the mechanism of arginine binding and release by TmArgBP have been discussed. Interestingly, both D1 and D2 retain the remarkable thermal/chemical stability of the parent protein. The analysis of the structural and dynamic properties of TmArgBP and of the individual domains highlights possible routes of domain communication. Finally, this study generated two interesting molecular tools, the two stable isolated domains that could be used in future investigations.

Testing the mutual information expansion of entropy with multivariate Gaussian distributions

The mutual information expansion (MIE) represents an approximation of the configurational entropy in terms of low-dimensional integrals. It is frequently employed to compute entropies from simulation data of large systems, such as macromolecules, for which brute-force evaluation of the full configurational integral is intractable. Here, we test the validity of MIE for systems consisting of more than m = 100 degrees of freedom (dofs). The dofs are distributed according to multivariate Gaussian distributions which were generated from protein structures using a variant of the anisotropic network model. For the Gaussian distributions, we have semi-analytical access to the configurational entropy as well as to all contributions of MIE. This allows us to accurately assess the validity of MIE for different situations. We find that MIE diverges for systems containing long-range correlations which means that the error of consecutive MIE approximations grows with the truncation order n for all tractable n ≪ m. This fact implies severe limitations on the applicability of MIE, which are discussed in the article. For systems with correlations that decay exponentially with distance, MIE represents an asymptotic expansion of entropy, where the first successive MIE approximations approach the exact entropy, while MIE also diverges for larger orders. In this case, MIE serves as a useful entropy expansion when truncated up to a specific truncation order which depends on the correlation length of the system.

Prediction of Protein Configurational Entropy (Popcoen)

A knowledge-based method for configurational entropy prediction of proteins is presented; this methodology is extremely fast, compared to previous approaches, because it does not involve any type of configurational sampling. Instead, the configurational entropy of a query fold is estimated by evaluating an artificial neural network, which was trained on molecular-dynamics simulations of ∼1000 proteins. The predicted entropy can be incorporated into a large class of protein software based on cost-function minimization/evaluation, in which configurational entropy is currently neglected for performance reasons. Software of this type is used for all major protein tasks such as structure predictions, proteins design, NMR and X-ray refinement, docking, and mutation effect predictions. Integrating the predicted entropy can yield a significant accuracy increase as we show exemplarily for native-state identification with the prominent protein software FoldX. The method has been termed Popcoen for Prediction of Protein Configurational Entropy. An implementation is freely available at http://fmc.ub.edu/popcoen/ .

Perplexing cooperative folding and stability of a low-sequence complexity, polyproline 2 protein lacking a hydrophobic core

The burial of hydrophobic side chains in a protein core generally is thought to be the major ingredient for stable, cooperative folding. Here, we show that, for the snow flea antifreeze protein (sfAFP), stability and cooperativity can occur without a hydrophobic core, and without α-helices or β-sheets. sfAFP has low sequence complexity with 46% glycine and an interior filled only with backbone H-bonds between six polyproline 2 (PP2) helices. However, the protein folds in a kinetically two-state manner and is moderately stable at room temperature. We believe that a major part of the stability arises from the unusual match between residue-level PP2 dihedral angle bias in the unfolded state and PP2 helical structure in the native state. Additional stabilizing factors that compensate for the dearth of hydrophobic burial include shorter and stronger H-bonds, and increased entropy in the folded state. These results extend our understanding of the origins of cooperativity and stability in protein folding, including the balance between solvent and polypeptide chain entropies.

Shedding light on the extra thermal stability of thermophilic proteins

An entropic stabilization mechanism has recently gained attention and credibility as the physical ground for the extra thermal stability of globular proteins from thermophilic microorganisms. An empirical result, obtained from the analysis of thermodynamic data for a large set of proteins, strengthens the general reliability of the theoretical approach originally devised to rationalize the occurrence of cold denaturation [Graziano, PCCP 2014, 16, 21755-21767]. It is shown that this theoretical approach can readily account for the entropic stabilization mechanism. On decreasing the conformational entropy gain associated with denaturation, the thermal stability of a model globular protein increases markedly.

Basic units of protein structure, folding, and function

Study of the hierarchy of domain structure with alternative sets of domains and analysis of discontinuous domains, consisting of remote segments of the polypeptide chain, raised a question about the minimal structural unit of the protein domain. The hypothesis on the decisive role of the polypeptide backbone in determining the elementary units of globular proteins have led to the discovery of closed loops. It is reviewed here how closed loops form the loop-n-lock structure of proteins, providing the foundation for stability and designability of protein folds/domain and underlying their co-translational folding. Simplified protein sequences are considered here with the aim to explore the basic principles that presumably dominated the folding and stability of proteins in the early stages of structural evolution. Elementary functional loops (EFLs), closed loops with one or few catalytic residues, are, in turn, units of the protein function. They are apparent descendants of the prebiotic ring-like peptides, which gave rise to the first functional folds/domains being fused in the beginning of the evolution of protein structure. It is also shown how evolutionary relations between protein functional superfamilies and folds delineated with the help of EFLs can contribute to establishing the rules for design of desired enzymatic functions. Generalized descriptors of the elementary functions are proposed to be used as basic units in the future computational design.

Vibrational entropy of a protein: large differences between distinct conformations

In this article, it is investigated whether vibrational entropy (VE) is an important contribution to the free energy of globular proteins at ambient conditions. VE represents the major configurational-entropy contribution of these proteins. By definition, it is an average of the configurational entropies of the protein within single minima of the energy landscape, weighted by their occupation probabilities. Its large part originates from thermal motion of flexible torsion angles giving rise to the finite peak widths observed in torsion angle distributions. While VE may affect the equilibrium properties of proteins, it is usually neglected in numerical calculations as its consideration is difficult. Moreover, it is sometimes believed that all well-packed conformations of a globular protein have similar VE anyway. Here, we measure explicitly the VE for six different conformations from simulation data of a test protein. Estimates are obtained using the quasi-harmonic approximation for three coordinate sets, Cartesian, bond-angle-torsion (BAT), and a new set termed rotamer-degeneracy lifted BAT coordinates by us. The new set gives improved estimates as it overcomes a known shortcoming of the quasi-harmonic approximation caused by multiply populated rotamer states, and it may serve for VE estimation of macromolecules in a very general context. The obtained VE values depend considerably on the type of coordinates used. However, for all coordinate sets we find large entropy differences between the conformations, of the order of the overall stability of the protein. This result may have important implications on the choice of free energy expressions used in software for protein structure prediction, protein design, and NMR refinement.

Role of Proteome Physical Chemistry in Cell Behavior

We review how major cell behaviors, such as bacterial growth laws, are derived from the physical chemistry of the cell’s proteins. On one hand, cell actions depend on the individual biological functionalities of their many genes and proteins. On the other hand, the common physics among proteins can be as important as the unique biology that distinguishes them. For example, bacterial growth rates depend strongly on temperature. This dependence can be explained by the folding stabilities across a cell’s proteome. Such modeling explains how thermophilic and mesophilic organisms differ, and how oxidative damage of highly charged proteins can lead to unfolding and aggregation in aging cells. Cells have characteristic time scales. For example, E. coli can duplicate as fast as 2–3 times per hour. These time scales can be explained by protein dynamics (the rates of synthesis and degradation, folding, and diffusional transport). It rationalizes how bacterial growth is slowed down by added salt. In the same way that the behaviors of inanimate materials can be expressed in terms of the statistical distributions of atoms and molecules, some cell behaviors can be expressed in terms of distributions of protein properties, giving insights into the microscopic basis of growth laws in simple cells.

Shortening a loop can increase protein native state entropy

Protein loops are essential structural elements that influence not only function but also protein stability and folding rates. It was recently reported that shortening a loop in the AcP protein may increase its native state conformational entropy. This effect on the entropy of the folded state can be much larger than the lower entropic penalty of ordering a shorter loop upon folding, and can therefore result in a more pronounced stabilization than predicted by polymer model for loop closure entropy. In this study, which aims at generalizing the effect of loop length shortening on native state dynamics, we use all-atom molecular dynamics simulations to study how gradual shortening a very long or solvent-exposed loop region in four different proteins can affect their stability. For two proteins, AcP and Ubc7, we show an increase in native state entropy in addition to the known effect of the loop length on the unfolded state entropy. However, for two permutants of SH3 domain, shortening a loop results only with the expected change in the entropy of the unfolded state, which nicely reproduces the observed experimental stabilization. Here, we show that an increase in the native state entropy following loop shortening is not unique to the AcP protein, yet nor is it a general rule that applies to all proteins following the truncation of any loop. This modification of the loop length on the folded state and on the unfolded state may result with a greater effect on protein stability.

Rigidity versus flexibility: the dilemma of understanding protein thermal stability

The role of fluctuations in protein thermostability has recently received considerable attention. In the current literature a dualistic picture can be found: thermostability seems to be associated with enhanced rigidity of the protein scaffold in parallel with the reduction of flexible parts of the structure. In contradiction to such arguments it has been shown by experimental studies and computer simulation that thermal tolerance of a protein is not necessarily correlated with the suppression of internal fluctuations and mobility. Both concepts, rigidity and flexibility, are derived from mechanical engineering and represent temporally insensitive features describing static properties, neglecting that relative motion at certain time scales is possible in structurally stable regions of a protein. This suggests that a strict separation of rigid and flexible parts of a protein molecule does not describe the reality correctly. In this work the concepts of mobility/flexibility versus rigidity will be critically reconsidered by taking into account molecular dynamics calculations of heat capacity and conformational entropy, salt bridge networks, electrostatic interactions in folded and unfolded states, and the emerging picture of protein thermostability in view of recently developed network theories. Last, but not least, the influence of high temperature on the active site and activity of enzymes will be considered.

Organization of the multiaminoacyl-tRNA synthetase complex and the cotranslational protein folding

Aminoacyl-tRNA synthetases (ARSs) play an essential role in the protein synthesis by catalyzing an attachment of their cognate amino acids to tRNAs. Unlike their prokaryotic counterparts, ARSs in higher eukaryotes form a multiaminoacyl-tRNA synthetase complex (MARS), consisting of the subset of ARS polypeptides and three auxiliary proteins. The intriguing feature of MARS complex is the presence of only nine out of twenty ARSs, specific for Arg, Asp, Gln, Glu, Ile, Leu, Lys, Met, and Pro, regardless of the organism, cell, or tissue types. Although existence of MARSs complex in higher eukaryotes has been already known for more than four decades, its functional significance remains elusive. We found that seven of the nine corresponding amino acids (Arg, Gln, Glu, Ile, Leu, Lys, and Met) together with Ala form a predictor of the protein α-helicity. Remarkably, all amino acids (besides Ala) in the predictor have the highest possible number of side-chain rotamers. Therefore, compositional bias of a typical α-helix can contribute to the helix's stability by increasing the entropy of the folded state. It also appears that position-specific α-helical propensity, specifically periodic alternation of charged and hydrophobic residues in the helices, may well be provided by the structural organization of the complex. Considering characteristics of MARS complex from the perspective of the α-helicity, we hypothesize that specific composition and structure of the complex represents a functional mechanism for coordination of translation with the fast and correct folding of amphiphilic α-helices.

Nonspecific yet decisive: Ubiquitination can affect the native-state dynamics of the modified protein

Ubiquitination is one of the most common post-translational modifications of proteins, and mediates regulated protein degradation among other cellular processes. A fundamental question regarding the mechanism of protein ubiquitination is whether and how ubiquitin affects the biophysical nature of the modified protein. For some systems, it was shown that the position of ubiquitin within the attachment site is quite flexible and ubiquitin does not specifically interact with its substrate. Nevertheless, it was revealed that polyubiquitination can decrease the thermal stability of the modified protein in a site-specific manner because of alterations of the thermodynamic properties of the folded and unfolded states. In this study, we used detailed atomistic simulations to focus on the molecular effects of ubiquitination on the native structure of the modified protein. As a model, we used Ubc7, which is an E2 enzyme whose in vivo ubiquitination process is well characterized and known to lead to degradation. We found that, despite the lack of specific direct interactions between the ubiquitin moiety and Ubc7, ubiquitination decreases the conformational flexibility of certain regions of the substrate Ubc7 protein, which reduces its entropy and thus destabilizes it. The strongest destabilizing effect was observed for systems in which Lys48-linked tetra-ubiquitin was attached to sites used for in vivo degradation. These results reveal how changes in the configurational entropy of the folded state may modulate the stability of the protein's native state. Overall, our results imply that ubiquitination can modify the biophysical properties of the attached protein in the folded state and that, in some proteins, different ubiquitination sites will lead to different biophysical outcomes. We propose that this destabilizing effect of polyubiquitin on the substrate is linked to the functions carried out by the modification, and in particular, regulatory control of protein half-life through proteasomal degradation.

Mining the entire Protein DataBank for frequent spatially cohesive amino acid patterns

BACKGROUND

The three-dimensional structure of a protein is an essential aspect of its functionality. Despite the large diversity in protein structures and functionality, it is known that there are common patterns and preferences in the contacts between amino acid residues, or between residues and other biomolecules, such as DNA. The discovery and characterization of these patterns is an important research topic within structural biology as it can give fundamental insight into protein structures and can aid in the prediction of unknown structures.

RESULTS

Here we apply an efficient spatial pattern miner to search for sets of amino acids that occur frequently in close spatial proximity in the protein structures of the Protein DataBank. This allowed us to mine for a new class of amino acid patterns, that we term FreSCOs (Frequent Spatially Cohesive Component sets), which feature synergetic combinations. To demonstrate the relevance of these FreSCOs, they were compared in relation to the thermostability of the protein structure and the interaction preferences of DNA-protein complexes. In both cases, the results matched well with prior investigations using more complex methods on smaller data sets.

CONCLUSIONS

The currently characterized protein structures feature a diverse set of frequent amino acid patterns that can be related to the stability of the protein molecular structure and that are independent from protein function or specific conserved domains.

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.

A role for A-to-I RNA editing in temperature adaptation

A-to-I RNA editing can recode mRNAs, giving organisms the option to express diverse, functionally distinct protein isoforms. Here, we propose that RNA editing is inherently geared for temperature adaptation because it tends to recode to smaller, less stabilizing amino acids. Studies on how editing affects protein function support this idea.

Stability and rigidity/flexibility-two sides of the same coin?

Protein molecules require both flexibility and rigidity for functioning. The fast and accurate prediction of protein rigidity/flexibility is one of the important problems in protein science. We have determined flexible regions for four homologous pairs from thermophilic and mesophilic organisms by two methods: the fast FoldUnfold which uses amino acid sequence and the time consuming MDFirst which uses three-dimensional structures. We demonstrate that both methods allow determining flexible regions in protein structure. For three of the four thermophile-mesophile pairs of proteins, FoldUnfold predicts practically the same flexible regions which have been found by the MD/First method. As expected, molecular dynamics simulations show that thermophilic proteins are more rigid in comparison to their mesophilic homologues. Analysis of rigid clusters and their decomposition provides new insights into protein stability. It has been found that the local networks of salt bridges and hydrogen bonds in thermophiles render their structure more stable with respect to fluctuations of individual contacts. Such network includes salt bridge triads Agr-Glu-Lys and Arg-Glu-Arg, or salt bridges (such as Arg-Glu) connected with hydrogen bonds. This ionic network connects alpha helices and rigidifies the structure. Mesophiles can be characterized by stand alone salt bridges and hydrogen bonds or small ionic clusters. Such difference in the network of salt bridges results in different flexibility of homologous proteins. Combining both approaches allows characterizing structural features in atomic detail that determine the rigidity/flexibility of a protein structure. This article is a part of a Special Issue entitled: The emerging dynamic view of proteins: Protein plasticity in allostery, evolution and self-assembly.

The effect of protein composition on hydration dynamics

Water dynamics at the surface of two homologous proteins with different thermal resistances is found to be unaffected by the different underlying amino-acid compositions, and when proteins are folded it responds similarly to temperature variations. Upon unfolding the water dynamics slowdown with respect to bulk decreases by a factor of two. Our findings are explained by the dominant topological perturbation induced by the protein on the water hydrogen bond dynamics.

Computational and statistical analyses of amino acid usage and physico-chemical properties of the twelve late embryogenesis abundant protein classes

Late Embryogenesis Abundant Proteins (LEAPs) are ubiquitous proteins expected to play major roles in desiccation tolerance. Little is known about their structure - function relationships because of the scarcity of 3-D structures for LEAPs. The previous building of LEAPdb, a database dedicated to LEAPs from plants and other organisms, led to the classification of 710 LEAPs into 12 non-overlapping classes with distinct properties. Using this resource, numerous physico-chemical properties of LEAPs and amino acid usage by LEAPs have been computed and statistically analyzed, revealing distinctive features for each class. This unprecedented analysis allowed a rigorous characterization of the 12 LEAP classes, which differed also in multiple structural and physico-chemical features. Although most LEAPs can be predicted as intrinsically disordered proteins, the analysis indicates that LEAP class 7 (PF03168) and probably LEAP class 11 (PF04927) are natively folded proteins. This study thus provides a detailed description of the structural properties of this protein family opening the path toward further LEAP structure - function analysis. Finally, since each LEAP class can be clearly characterized by a unique set of physico-chemical properties, this will allow development of software to predict proteins as LEAPs.

Computational design of a PDZ domain peptide inhibitor that rescues CFTR activity

The cystic fibrosis transmembrane conductance regulator (CFTR) is an epithelial chloride channel mutated in patients with cystic fibrosis (CF). The most prevalent CFTR mutation, ΔF508, blocks folding in the endoplasmic reticulum. Recent work has shown that some ΔF508-CFTR channel activity can be recovered by pharmaceutical modulators (“potentiators” and “correctors”), but ΔF508-CFTR can still be rapidly degraded via a lysosomal pathway involving the CFTR-associated ligand (CAL), which binds CFTR via a PDZ interaction domain. We present a study that goes from theory, to new structure-based computational design algorithms, to computational predictions, to biochemical testing and ultimately to epithelial-cell validation of novel, effective CAL PDZ inhibitors (called “stabilizers”) that rescue ΔF508-CFTR activity. To design the “stabilizers”, we extended our structural ensemble-based computational protein redesign algorithm to encompass protein-protein and protein-peptide interactions. The computational predictions achieved high accuracy: all of the top-predicted peptide inhibitors bound well to CAL. Furthermore, when compared to state-of-the-art CAL inhibitors, our design methodology achieved higher affinity and increased binding efficiency. The designed inhibitor with the highest affinity for CAL (kCAL01) binds six-fold more tightly than the previous best hexamer (iCAL35), and 170-fold more tightly than the CFTR C-terminus. We show that kCAL01 has physiological activity and can rescue chloride efflux in CF patient-derived airway epithelial cells. Since stabilizers address a different cellular CF defect from potentiators and correctors, our inhibitors provide an additional therapeutic pathway that can be used in conjunction with current methods.

Cystic fibrosis (CF) is an inherited disease that causes the body to produce thick mucus that clogs the lungs and obstructs the breakdown and absorption of food. The cystic fibrosis transmembrane conductance regulator (CFTR) is mutated in CF patients, and the most common mutation causes three defects in CFTR: misfolding, decreased function, and rapid degradation. Drugs are currently being studied to correct the first two CFTR defects, but the problem of rapid degradation remains. Recently, key protein-protein interactions have been discovered that implicate the protein CAL in CFTR degradation. Here we have developed new computational protein design algorithms and used them to successfully predict peptide inhibitors of the CAL-CFTR interface. Our algorithm uses a structural ensemble-based evaluation of protein sequences and conformations to calculate accurate predictions of protein-peptide binding affinities. The algorithm is general and can be applied to a wide variety of protein-protein interface designs. All of our designed inhibitors bound CAL with high affinity. We tested our top binding peptide and observed that the inhibitor could successfully rescue CFTR function in CF patient-derived epithelial cells. Our designed inhibitors provide a novel therapeutic path which could be used in combination with existing CF therapeutics for additive benefit.

Folding without charges

Surface charges of proteins have in several cases been found to function as "structural gatekeepers," which avoid unwanted interactions by negative design, for example, in the control of protein aggregation and binding. The question is then if side-chain charges, due to their desolvation penalties, play a corresponding role in protein folding by avoiding competing, misfolded traps? To find out, we removed all 32 side-chain charges from the 101-residue protein S6 from Thermus thermophilus. The results show that the charge-depleted S6 variant not only retains its native structure and cooperative folding transition, but folds also faster than the wild-type protein. In addition, charge removal unleashes pronounced aggregation on longer timescales. S6 provides thus an example where the bias toward native contacts of a naturally evolved protein sequence is independent of charges, and point at a fundamental difference in the codes for folding and intermolecular interaction: specificity in folding is governed primarily by hydrophobic packing and hydrogen bonding, whereas solubility and binding relies critically on the interplay of side-chain charges.

Proteome adaptation to high temperatures in the ectothermic hydrothermal vent Pompeii worm

Taking advantage of the massive genome sequencing effort made on thermophilic prokaryotes, thermal adaptation has been extensively studied by analysing amino acid replacements and codon usage in these unicellular organisms. In most cases, adaptation to thermophily is associated with greater residue hydrophobicity and more charged residues. Both of these characteristics are positively correlated with the optimal growth temperature of prokaryotes. In contrast, little information has been collected on the molecular 'adaptive' strategy of thermophilic eukaryotes. The Pompeii worm A. pompejana, whose transcriptome has recently been sequenced, is currently considered as the most thermotolerant eukaryote on Earth, withstanding the greatest thermal and chemical ranges known. We investigated the amino-acid composition bias of ribosomal proteins in the Pompeii worm when compared to other lophotrochozoans and checked for putative adaptive changes during the course of evolution using codon-based Maximum likelihood analyses. We then provided a comparative analysis of codon usage and amino-acid replacements from a greater set of orthologous genes between the Pompeii worm and Paralvinella grasslei, one of its closest relatives living in a much cooler habitat. Analyses reveal that both species display the same high GC-biased codon usage and amino-acid patterns favoring both positively-charged residues and protein hydrophobicity. These patterns may be indicative of an ancestral adaptation to the deep sea and/or thermophily. In addition, the Pompeii worm displays a set of amino-acid change patterns that may explain its greater thermotolerance, with a significant increase in Tyr, Lys and Ala against Val, Met and Gly. Present results indicate that, together with a high content in charged residues, greater proportion of smaller aliphatic residues, and especially alanine, may be a different path for metazoans to face relatively 'high' temperatures and thus a novelty in thermophilic metazoans.

Physical limits of cells and proteomes

What are the physical limits to cell behavior? Often, the physical limitations can be dominated by the proteome, the cell's complement of proteins. We combine known protein sizes, stabilities, and rates of folding and diffusion, with the known protein-length distributions P(N) of proteomes (Escherichia coli, yeast, and worm), to formulate distributions and scaling relationships in order to address questions of cell physics. Why do mesophilic cells die around 50 °C? How can the maximal growth-rate temperature (around 37 °C) occur so close to the cell-death temperature? The model shows that the cell's death temperature coincides with a denaturation catastrophe of its proteome. The reason cells can function so well just a few degrees below their death temperature is because proteome denaturation is so cooperative. Why are cells so dense-packed with protein molecules (about 20% by volume)? Cells are packed at a density that maximizes biochemical reaction rates. At lower densities, proteins collide too rarely. At higher densities, proteins diffuse too slowly through the crowded cell. What limits cell sizes and growth rates? Cell growth is limited by rates of protein synthesis, by the folding rates of its slowest proteins, and--for large cells--by the rates of its protein diffusion. Useful insights into cell physics may be obtainable from scaling laws that encapsulate information from protein knowledge bases.