Role of methyl groups in dynamics and evolution of biomolecules

Recent studies have discovered strong differences between the dynamics of nucleic acids (RNA and DNA) and proteins, especially at low hydration and low temperatures. This difference is caused primarily by dynamics of methyl groups that are abundant in proteins, but are absent or very rare in RNA and DNA. In this paper, we present a hypothesis regarding the role of methyl groups as intrinsic plasticizers in proteins and their evolutionary selection to facilitate protein dynamics and activity. We demonstrate the profound effect methyl groups have on protein dynamics relative to nucleic acid dynamics, and note the apparent correlation of methyl group content in protein classes and their need for molecular flexibility. Moreover, we note the fastest methyl groups of some enzymes appear around dynamical centers such as hinges or active sites. Methyl groups are also of tremendous importance from a hydrophobicity/folding/entropy perspective. These significant roles, however, complement our hypothesis rather than preclude the recognition of methyl groups in the dynamics and evolution of biomolecules.     

Biological synthesis of circular polypeptides

Here, we review the use of different biochemical approaches for biological synthesis of circular or backbone-cyclized proteins and peptides. These methods allow the production of circular polypeptides either in vitro or in vivo using standard recombinant DNA expression techniques. Protein circularization can significantly impact protein engineering and research in protein folding. Basic polymer theory predicts that circularization should lead to a net thermodynamic stabilization of a folded protein by reducing the entropy associated with the unfolded state. Protein cyclization also provides a valuable tool for exploring the effects of topology on protein folding kinetics. Furthermore, the biological production of cyclic polypeptides makes possible the production of cyclic polypeptide libraries. The generation of such libraries, which was previously restricted to the domain of synthetic chemists, now offers biologists access to highly diverse and stable molecular libraries for probing protein structure and function.     

Two different late embryogenesis abundant proteins from Arabidopsis thaliana contain specific domains that inhibit Escherichia coli growth

Late embryogenesis abundant (LEA) proteins constitute a set of proteins widespread in the plant kingdom that show common physicochemical properties such as high hydrophilicity and high content of small amino acid residues such as glycine, alanine, and serine. Typically, these proteins accumulate in response to water deficit conditions imposed by the environment or during plant normal development. In this work, we show that the over-expression in Escherichia coli of proteins of the LEA 2 and the LEA 4 families from Arabidopsis thaliana leads to inhibition of bacterial growth and that this effect is dependent on discrete regions of the proteins. Our data indicate that their antimicrobial effect is achieved through their interaction with intracellular targets. The relevance of the cationic nature and the predicted structural organization of particular protein domains in this detrimental effect on the bacteria growth process is discussed.     

Medium throughput cage state stability screen of conditions for the generation of gold nanoparticles encapsulated within a mini-ferritin

Methods for the generation of nanoparticles encapsulated within cage proteins, such as ferritins, provide particles with low polydispersities due to size constraint by the cage. The proteins can provide enhanced water solubility to enable biological applications and affinity and identification tags to facilitate delivery or the assembly of advanced materials. Many effective methods have been developed, however, they are often impeded by cage protein instability in the presence of reagents or conditions for formation of the nanoparticles. Although the stability of ferritin cage quaternary structure can be enhanced, application of ferritins to materials science remains limited by unpredictable behaviour. Recently, we reported a medium throughput technique to directly detect the ferritin cage state. Herein, we expand this strategy to screen conditions commonly used for the formation of gold nanoparticles. Not only do we report nanoparticle formation conditions that permit ferritin stability, we establish a general screening strategy based on protein cage stability that could be applied to other protein cages or for the generation of other types of particles.     

Enhanced crystal packing due to solvent reorganization through reductive methylation of lysine residues in oxidoreductase from <Emphasis Type="Italic">Streptococcus pneumoniae</Emphasis>

Protein crystallization is in part driven by the changes in the entropy of the system, but opinions differ as to whether the solute (protein) or solvent (water) molecules make more of a contribution to the overall entropic change. Methylation of lysine residues in proteins has been used to enhance protein crystallization. We investigated using molecular dynamics simulations with explicit solvent molecules, the behavior of several native proteins and their methylated counterparts chosen from an earlier large-scale study. Methylated lysines are capable of making a variety of interactions including H-bonds with protein residues and solvent. We demonstrate that methylation on the lysine slightly increases its side chain conformational entropy by about 3.5 J mol⁻¹ K⁻¹. Analysis of the radial and spatial distributions of the water molecules around the methylated lysine surface in oxidoreductase from Streptococcus pneumoniae revealed a larger sphere of water molecules with low entropy, as compared with solvent associated with unmethylated lysine. If methylated lysine were to make interactions at the protein–protein interface, the low-entropy water molecules associated with methylated lysines would be released, resulting in a gain of entropy. We show that this gain more than compensates for the loss of protein entropy. Therefore, we propose that lysine methylation favors the formation of crystals through solvent entropic gain.     

Heat capacity of protein folding.

We construct a Hamiltonian for a single domain protein where the contact enthalpy and the chain entropy decrease linearly with the number of native contacts. The hydration effect upon protein unfolding is included by modeling water as ideal dipoles that are ordered around the unfolded surfaces, where the influence of these surfaces, covered with an "ice-like" shell of water, is represented by an effective field that directs the water dipoles. An intermolecular pair interaction between water molecules is also introduced. The heat capacity of the model exhibits, the common feature of small globular proteins, two peaks corresponding to cold and warm unfolding, respectively. By introducing ad hoc vibrational modes, we obtain quantitatively good accordance with experiments on myoglobin.     

Quantifying the Entropy of Binding for Water Molecules in Protein Cavities by Computing Correlations

Protein structural analysis demonstrates that water molecules are commonly found in the internal cavities of proteins. Analysis of experimental data on the entropies of inorganic crystals suggests that the entropic cost of transferring such a water molecule to a protein cavity will not typically be greater than 7.0 cal/mol/K per water molecule, corresponding to a contribution of approximately +2.0 kcal/mol to the free energy. In this study, we employ the statistical mechanical method of inhomogeneous fluid solvation theory to quantify the enthalpic and entropic contributions of individual water molecules in 19 protein cavities across five different proteins. We utilize information theory to develop a rigorous estimate of the total two-particle entropy, yielding a complete framework to calculate hydration free energies. We show that predictions from inhomogeneous fluid solvation theory are in excellent agreement with predictions from free energy perturbation (FEP) and that these predictions are consistent with experimental estimates. However, the results suggest that water molecules in protein cavities containing charged residues may be subject to entropy changes that contribute more than +2.0 kcal/mol to the free energy. In all cases, these unfavorable entropy changes are predicted to be dominated by highly favorable enthalpy changes. These findings are relevant to the study of bridging water molecules at protein-protein interfaces as well as in complexes with cognate ligands and small-molecule inhibitors.     

Enthalpy–entropy compensation phenomena in water solutions of proteins and small molecules: A ubiquitous properly of water

This article presents evidence for the existence of a specific linear relationship between the entropy change and the enthalpy change in a variety of processes of small solutes in water solution. The processes include solvation of ions and nonelectrolytes, hydrolysis, oxidation–reduction, ionization of weak electrolytes, and quenching of indole fluorescence among others. The values of the proportionality constant, called the compensation temperature, lie in a relatively narrow range, from about 250 to 315 °K, for all these processes. Such behavior can be a consequence of experimental errors but for a number of the processes the precision of the data is sufficient to show that the enthalpy–entropy compensation pattern is real. It is tentatively concluded that the pattern is real, very common and a consequence of the properties of liquid water as a solvent regardless of the solutes and the solute processes studied. As such the phenomenon requires that theoretical treatments of solute processes in water be expanded by inclusion of a specific treatment of the characteristic of water responsible for compensation behavior. The possible bases of the effect are proposed to be temperature-independent heat-capacity changes and/or shifts in concentrations of the two phenomenologically significant species of water. The relationship of these alternatives to the two-state process of water suggested by spectroscopic and relaxation studies is examined. The existence of a similar and probably identical relationship between enthalpy and entropy change in a variety of protein reactions suggests that liquid water plays a direct role in many protein processes and may be a common participant in the physiological function of proteins. It is proposed that the linear enthalpy–entropy relationship be used as a diagnostic test for the participation of water in protein processes. On this basis the catalytic processes of chymotrypsin and acetylcholinesterase are dominated by the properties of bulk water. The binding of oxygen by hemoglobin may fall in the same category. Similarities and differences in the behavior of small-solute and protein processes are examined to show how they may be related. No positive conclusions are established, but it is possible that protein processes are coupled to water via expansions and contractions of the protein and that in general the special pattern of enthalpy–entropy compensation is a consequent of the properties of water which require that expansions and contractions of solutes effect changes in the free volume of the nearby liquid water. It is shown that proteins can be expected to respond to changes in nearby water and interfacial free energy by expansions and contractions. Such responses may explain a variety of currently unexplained characteristics of protein solutions. More generally, the enthalpy–entropy compensation pattern appears to be the thermodynamic manifestation of “structure making” and “structure breaking,” operationally defined terms much used in discussions of water solutions. If so, the compensation pattern is ubiquitous and requires re-examination of a large body of molecular interpretations derived from quantitative studies of processes in water. Theories of processes in water may have to be expanded to accommodate this aspect of water behavior.     

Allostery without conformational change

A general model is presented whereby lignand-induced changes in protein dynamics could produce allosteric communication between distinct binding sites, even in the absence of a macromolecular conformational change. Theoretical analysis, based on the statistical thermodynamics of ligand binding, shows that cooperative interaction free energies amounting to several kJ · mol^-1 may be generated by this means. The effect arises out of the possible changes in frequencies and amplitudes of macromolecular thermal fluctuations in response to ligand attachment, and can involve all forms of dynamic behaviour, ranging from highly correlated, low-frequency normal mode vibrations to random local anharmonic motions of individual atoms or groups. Dynamic allostery of this form is primarily an entropy effect, and we derive approximate expressions which might allow the magnitude of the interaction in real systems to be calculated directly from experimental observations such as changes in normal mode frequencies and mean-square atomic displacements. Long-range influence of kinetic processes at different sites might also be mediated by a similar mechanism. We suggest that proteins and other biological macromolecules may have evolved to take functional advantage not only of mean conformational states but also of the inevitable thermal fluctuations about the mean.     

Multiple Post-translational Modifications Affect Heterologous Protein Synthesis

Post-translational modifications (PTMs) are required for proper folding of many proteins. The low capacity for PTMs hinders the production of heterologous proteins in the widely used prokaryotic systems of protein synthesis. Until now, a systematic and comprehensive study concerning the specific effects of individual PTMs on heterologous protein synthesis has not been presented. To address this issue, we expressed 1488 human proteins and their domains in a bacterial cell-free system, and we examined the correlation of the expression yields with the presence of multiple PTM sites bioinformatically predicted in these proteins. This approach revealed a number of previously unknown statistically significant correlations. Prediction of some PTMs, such as myristoylation, glycosylation, palmitoylation, and disulfide bond formation, was found to significantly worsen protein amenability to soluble expression. The presence of other PTMs, such as aspartyl hydroxylation, C-terminal amidation, and Tyr sulfation, did not correlate with the yield of heterologous protein expression. Surprisingly, the predicted presence of several PTMs, such as phosphorylation, ubiquitination, SUMOylation, and prenylation, was associated with the increased production of properly folded soluble proteins. The plausible rationales for the existence of the observed correlations are presented. Our findings suggest that identification of potential PTMs in polypeptide sequences can be of practical use for predicting expression success and optimizing heterologous protein synthesis. In sum, this study provides the most compelling evidence so far for the role of multiple PTMs in the stability and solubility of heterologously expressed recombinant proteins.     

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. Proteins 2015; 83:2137–2146. © 2015 Wiley Periodicals, Inc.     

Identification of proteins interacting with the catalytic subunit of PP2A by proteomics

The protein phosphatase 2A (PP2A) is a serine/threonine phosphatase involved in the regulation of multiple signaling pathways including the Wnt/beta-catenin and the ERK pathways. To understand the complex signaling networking associated with PP2A, we searched proteins interacting with the catalytic subunit of protein phosphatase 2A (PP2Ac) by a pull-down analysis followed by 2-D gel electrophoresis and proteomic analyses. The probability of identification of the proteins interacting with PP2Ac was increased by searching proteins differently interacting with PP2Ac according to stimulation of Wnt3a, which regulates both the Wnt/beta-catenin and the ERK pathways. Around 100 proteins, pulled-down by His-tagged PP2Ac, were identified in 2-D gels stained with CBB. By MALDI-TOF-MS analyses of 45 protein spots, we identified several proteins that were previously known to interact with PP2A, such as Axin and CaMK IV. In addition, we also identified many proteins that potentially interact with PP2Ac. The interactions of several candidate proteins, such as tuberous sclerosis complex 2, RhoB, R-Ras, and Nm23H2, with PP2Ac, were confirmed by in vitro binding analyses and/or coimmunoprecipitation experiments.     

The effect of context on the folding of β-hairpins

Small β-hairpin peptides have been widely used as models for the folding of β-sheets. But how applicable is the folding of such models to β-structure in larger proteins with conventional hydrophobic cores? Here we present multiple unfolding simulations of three such proteins that contain the WW domain double hairpin β-sheet motif: cold shock protein A (CspA), cold shock protein B (CspB) and glucose permease IIA domain. We compare the behavior of the free motif in solution and in the context of proteins of different size and architecture. Both Csp proteins lost contacts between the double-hairpin motif and the protein core as the first step of unfolding and proceeded to unfold with loss of the third β-strand, similar to the isolated WW domain. The glucose permease IIA domain is a larger protein and the contacts between the motif and the core were not lost as quickly. Instead the unfolding pathway of glucose permease IIA followed a different pathway with β1 pulling away from the sheet first. Interestingly, when the double hairpin motif was excised from the glucose permease IIA domain and simulated in isolation in water it unfolded by the same pathway as the WW domain, indicating that it is tertiary interactions with the protein that alter the motif’s unfolding not a sequence dependent effect on its intrinsic unfolding behavior. With respect to the unfolding of the hairpins, there was no consistent order to the loss of hydrogen bonds between the β-strands in the hairpins in any of the systems. Our results show that while the folding behavior of the isolated WW domain is generally consistent with the double hairpin motif’s behavior in the cold shock proteins, it is not the case for the glucose permease IIA domain. So, one must be cautious in extrapolating findings from model systems to larger more complicated proteins where tertiary interactions can overwhelm intrinsic behavior.     

Tracing Primordial Protein Evolution through Structurally Guided Stepwise Segment Elongation

The understanding of how primordial proteins emerged has been a fundamental and longstanding issue in biology and biochemistry. For a better understanding of primordial protein evolution, we synthesized an artificial protein on the basis of an evolutionary hypothesis, segment-based elongation starting from an autonomously foldable short peptide. A 10-residue protein, chignolin, the smallest foldable polypeptide ever reported, was used as a structural support to facilitate higher structural organization and gain-of-function in the development of an artificial protein. Repetitive cycles of segment elongation and subsequent phage display selection successfully produced a 25-residue protein, termed AF.2A1, with nanomolar affinity against the Fc region of immunoglobulin G. AF.2A1 shows exquisite molecular recognition ability such that it can distinguish conformational differences of the same molecule. The structure determined by NMR measurements demonstrated that AF.2A1 forms a globular protein-like conformation with the chignolin-derived β-hairpin and a tryptophan-mediated hydrophobic core. Using sequence analysis and a mutation study, we discovered that the structural organization and gain-of-function emerged from the vicinity of the chignolin segment, revealing that the structural support served as the core in both structural and functional development. Here, we propose an evolutionary model for primordial proteins in which a foldable segment serves as the evolving core to facilitate structural and functional evolution. This study provides insights into primordial protein evolution and also presents a novel methodology for designing small sized proteins useful for industrial and pharmaceutical applications.     

Water's potential role: Insights from studies of the p53 core domain

Soluble proteins with amyloidogenic propensity such as the tumor suppressor protein p53 have high proportion of incompletely desolvated backbone H bonds (HB). Such bonds are vulnerable to water attack, thus potentially leading to the misfolding of these proteins. However, it is still not clear how the surrounding solvent influences the protein native states. To address this, systematic surveys by molecular dynamics simulations and entropy analysis were performed on the p53 core domain in this work. We examined seven wild/mutant X-ray structures and observed two types of water-network hydration in three "hot hydration centers" (DNA- or small molecule- binding surfaces of the p53 core domain). The "tight" water, resulting from the local collective hydrogen-bond interactions, is probably fundamental to the protein structural stability. The second type of water is highly "dynamical" and exchanges very fast within the bulk solution, which is unambiguously assisted by the local protein motions. An entropy mapping of the solvent around the protein and a temperature perturbation analysis further present the main features of the p53 hydration network. The particular environment created by different water molecules around the p53 core domain also partly explains the structural vulnerabilities of this protein.     

A Group 6 Late Embryogenesis Abundant Protein from Common Bean Is a Disordered Protein with Extended Helical Structure and Oligomer-forming Properties

Late embryogenesis-abundant proteins accumulate to high levels in dry seeds. Some of them also accumulate in response to water deficit in vegetative tissues, which leads to a remarkable association between their presence and low water availability conditions. A major sub-group of these proteins, also known as typical LEA proteins, shows high hydrophilicity and a high percentage of glycine and other small amino acid residues, distinctive physicochemical properties that predict a high content of structural disorder. Although all typical LEA proteins share these characteristics, seven groups can be distinguished by sequence similarity, indicating structural and functional diversity among them. Some of these groups have been extensively studied; however, others require a more detailed analysis to advance in their functional understanding. In this work, we report the structural characterization of a group 6 LEA protein from a common bean (Phaseolus vulgaris L.) (PvLEA6) by circular dichroism and nuclear magnetic resonance showing that it is a disordered protein in aqueous solution. Using the same techniques, we show that despite its unstructured nature, the addition of trifluoroethanol exhibited an intrinsic potential in this protein to gain helicity. This property was also promoted by high osmotic potentials or molecular crowding. Furthermore, we demonstrate that PvLEA6 protein is able to form soluble homo-oligomeric complexes that also show high levels of structural disorder. The association between PvLEA6 monomers to form dimers was shown to occur in plant cells by bimolecular fluorescence complementation, pointing to the in vivo functional relevance of this association.     

Small molecules that target phosphorylation dependent protein–protein interaction

Protein–protein interaction is one of the key events in the signal transduction pathway. The interaction changes the conformations, activities, localization and stabilities of the proteins, and transduces the signal to the next step. Frequently, this interaction occurs upon the protein phosphorylation. When upstream signals are stimulated, protein kinase(s) is/are activated and phosphorylate(s) their substrates, and induce the phosphorylation dependent protein–protein interaction. For this interaction, several domains in proteins are known to specifically recognize the phosphorylated residues of target proteins. These specific domains for interaction are important in the progression of the diseases caused by disordered signal transduction such as cancer. Thus small molecules that modulate this interaction are attractive lead compounds for the treatment of such diseases. In this review, we focused on three examples of phosphorylation dependent protein–protein interaction modules (14-3-3, polo box domain of Plk1 and F-box proteins in SCF ubiquitin ligases) and summarize small molecules that modulate their interaction. We also introduce our original screening system to identify such small molecules.     

The crowd you're in with: Effects of different types of crowding agents on protein aggregation

The intracellular environment contains high concentrations of macromolecules occupying up to 30% of the total cellular volume. Presence of these macromolecules decreases the effective volume available for the proteins in the cell and thus increases the effective protein concentrations and stabilizes the compact protein conformations. Macromolecular crowding created by various macromolecules such as proteins, nucleic acids, and carbohydrates has been shown to have a significant effect on a variety of cellular processes including protein aggregation. Most studies of macromolecular crowding have used neutral, flexible polysaccharides that function primarily via excluded volume effect as model crowding agents. Here we have examined the effects of more rigid polysaccharides on protein structure and aggregation. Our results indicate that rigid and flexible polysaccharides influence protein aggregation via different mechanisms and suggest that, in addition to excluded volume effect, changes in solution viscosity and non-specific protein–polymer interactions influence the structure and dynamics of proteins in crowded environments.