Hydrophobic collapse in (in silico) protein folding

Heat-induced agglomeration of water-soluble cod proteins toward gelled structures

Cod proteins (CPs) have potential applications in designing desirable gel-based products, and this study aimed to unravel their heat-induced aggregation pattern and further probe the roles in protein gels. SDS-PAGE analysis indicated that high-precipitation-coefficient aggregates (HPCAs) of CPs aggregates were composed of considerable polymers of myosin heavy chains and actin, and their low-precipitation-coefficient aggregates (LPCAs) contained myosin light chains and tropomyosin. Studies from correlation analysis between the structure and aggregation kinetics revealed that the generation of β-sheet and SS bonds were responsible for their spontaneous thermal aggregation induced by heating temperature and protein concentration, respectively. Additionally, as protein denaturation ratio increased, more and larger HPCAs were formed, which was evidenced driving the network formation of protein gels and resulting in higher storage modulus (G') values. These novel findings may be applicable to other animal proteins for better tailoring the manufacturing of muscle gel-based products.

Diffusive dynamics of a model protein chain in solution

A Markov state model is a powerful tool that can be used to track the evolution of populations of configurations in an atomistic representation of a protein. For a coarse-grained linear chain model with discontinuous interactions, the transition rates among states that appear in the Markov model when the monomer dynamics is diffusive can be determined by computing the relative entropy of states and their mean first passage times, quantities that are unchanged by the specification of the energies of the relevant states. In this paper, we verify the folding dynamics described by a diffusive linear chain model of the crambin protein in three distinct solvent systems, each differing in complexity: a hard-sphere solvent, a solvent undergoing multi-particle collision dynamics, and an implicit solvent model. The predicted transition rates among configurations agree quantitatively with those observed in explicit molecular dynamics simulations for all three solvent models. These results suggest that the local monomer-monomer interactions provide sufficient friction for the monomer dynamics to be diffusive on timescales relevant to changes in conformation. Factors such as structural ordering and dynamic hydrodynamic effects appear to have minimal influence on transition rates within the studied solvent densities.

Constructing a strongly interacting Pea-Cod binary protein system by introducing metal cations toward enhanced gelling properties

Developing prospective plant-animal binary protein systems with desirable nutritional and rheological properties stands as a significant and challenging pursuit within the food industry. Our understanding of the effect of adding salt on the aggregation behavior of food proteins is currently based on single model protein systems, however, this knowledge is rather limited following binary protein systems. Herein, various ionic strength settings are used to mitigate the repulsive forces between pea-cod mixed proteins during the thermal process, which further benefits the construction of a strengthened gel network. Transmission electron microscopy (TEM) and dynamic light scattering (DLS) collectively demonstrated that larger heat-induced protein aggregates were formed, which increased in size with higher ionic strength. In the presence of 2.5 mM CaCl2 and 50 mM NaCl, the disulfide bonds significantly increased from 19.3 to 27.53 and 30.5 μM/g, respectively. Notably, similar aggregation behavior could be found when introducing 2.5 mM CaCl2 or 25 mM NaCl, due to the enhanced aggregation tendency by specific binding of Ca2+ to proteins. With relevance to the strengthened cross-links between protein molecules, salt endowed composite gels with preferable gelling properties, evidenced by increased storage modulus. Additionally, the gelling temperature of mixed proteins decreased below 50 °C at elevated ionic strength. Simultaneously, the proportion of network proteins in composite gels increased remarkably from 82.05 % to 93.61 % and 92.31 % upon adding 5.0 mM CaCl2 and 100 mM NaCl, respectively. The findings provide a valuable foundation for designing economically viable and health-oriented plant-animal binary protein systems.

Microwave treatment on structure and digestibility characteristics of Spirulina platensis protein

As a novel protein resource, the low digestibility of Spirulina platensis protein (SPP) limits its large-scale application. From the perspective of food processing methods, different heating treatments were explored to improve the structure and digestibility of SPP. In this study, SPP was heated by water bath and microwave at the same heating rate and heating temperature. Microwave accelerated protein denaturation and structure unfolded as the heating intensity increases, causing more exposed hydrophobic residues and enhancing surface hydrophobicity. The data of free sulfhydryl group, particle size, and gel electrophoresis, showed that microwave treatment promoted the formation of protein aggregates. The structural changes can potentially improve the accessibility of digestive enzymes, promote the in vitro digestibility rate, and further accelerate the production of small molecular peptides and the release of free amino acids. This study provided an innovative approach to improve the digestibility and therefore the utilization efficiency of SPP.

Mechanisms of heat-mediated aggregation behavior of water-soluble cod protein

Cod proteins (CPs) are considered potential functional ingredients for developing gel-based foods, but present studies on the aggregation behavior of CPs upon heating remain limited. With this respect, the heat-induced aggregation kinetics of CPs at a subunit level was investigated. Based on different centrifugal forces, CPs aggregates were divided into three fractions: large-sized, intermediary-sized, and small-sized aggregates. SDS-PAGE and diagonal SDS-PAGE indicated that myosin heavy chains exhibited a higher affinity with actin to form intermediary-sized and large-sized aggregates; tropomyosin and myosin light chains were hardly engaged in the thermal aggregation and formed small-sized aggregates. The highly-polymerized aggregates adopted considerable transitions of helix-to-sheet in protein structures, whereas the structure of small-sized aggregates featured substantial helix-coil transitions. Furthermore, molecular interactions at different heating stages were revealed. These novel insights might advance our knowledge on the heat-induced aggregation behavior of CPs and provide fundamental information for the application of CPs in gel-based foods.

Evaluation of the roles of hydrophobic residues in the N-terminal region of archaeal trehalase in its folding

Thermoplasma trehalase Tvn1315 is predicted to be composed of a β-sandwich domain (BD) and a catalytic domain (CD) based on the structure of the bacterial GH15 family glucoamylase (GA). Tvn1315 as well as Tvn1315 (Δ5), in which the 5 N-terminal amino acids are deleted, could be expressed in Escherichia coli as active enzymes, but deletion of 10 residues (Δ10) led to inclusion body formation. To further investigate the role of the N-terminal region of BD, we constructed five mutants of Δ5, in which each of the 5th to 10th residues of the N-terminus of Tvn1315 was mutated to Ala. Every mutant protein could be recovered in soluble form, but only a small fraction of the Y9A mutant was recovered in the soluble fraction. The Y9A mutant recovered in soluble form had similar specific activity to the other proteins. Subsequent mutation analysis at the 9th position of Tvn1315 in Δ5 revealed that aromatic as well as bulky hydrophobic residues could function properly, but residues with hydroxy groups impaired the solubility. Similar results were obtained with mutants based on untruncated Tvn1315. When the predicted BD, Δ5BD, Δ10BD, and BD mutants were expressed, the Δ10BD protein formed inclusion bodies, and the BD mutants behaved similarly to the Δ5 and full-length enzyme mutants. These results suggest that the hydrophobic region is involved in the solubilization of BD during the folding process. Taken together, these results indicate that the solubility of CD depends on BD folding. KEY POINTS: • N-terminal hydrophobic region of the BD is involved in the protein folding. • The N-terminal hydrophobic region of the BD is also involved in the BD folding. • BD is able to weakly interact with the insoluble β-glucan.

A benchmark of optimally folded protein structures using integer programming and the 3D-HP-SC model

The Protein Structure Prediction (PSP) problem comprises, among other issues, forecasting the three-dimensional native structure of proteins using only their primary structure information. Most computational studies in this area use synthetic data instead of real biological data. However, the closer to the real-world, the more the impact of results and their applicability. This work presents 17 real protein sequences extracted from the Protein Data Bank for a benchmark to the PSP problem using the tri-dimensional Hydrophobic-Polar with Side-Chains model (3D-HP-SC). The native structure of these proteins was found by maximizing the number of hydrophobic contacts between the side-chains of amino acids. The problem was treated as an optimization problem and solved by means of an Integer Programming approach. Although the method optimally solves the problem, the processing time has an exponential trend. Therefore, due to computational limitations, the method is a proof-of-concept and it is not applicable to large sequences. For unknown sequences, an upper bound of the number of hydrophobic contacts (using this model) can be found, due to a linear relationship with the number of hydrophobic residues. The comparison between the predicted and the biological structures showed that the highest similarity between them was found with distance thresholds around 5.2-8.2 Å. Both the dataset and the programs developed will be freely available to foster further research in the area.

eMatchSite: sequence order-independent structure alignments of ligand binding pockets in protein models

Detecting similarities between ligand binding sites in the absence of global homology between target proteins has been recognized as one of the critical components of modern drug discovery. Local binding site alignments can be constructed using sequence order-independent techniques, however, to achieve a high accuracy, many current algorithms for binding site comparison require high-quality experimental protein structures, preferably in the bound conformational state. This, in turn, complicates proteome scale applications, where only various quality structure models are available for the majority of gene products. To improve the state-of-the-art, we developed eMatchSite, a new method for constructing sequence order-independent alignments of ligand binding sites in protein models. Large-scale benchmarking calculations using adenine-binding pockets in crystal structures demonstrate that eMatchSite generates accurate alignments for almost three times more protein pairs than SOIPPA. More importantly, eMatchSite offers a high tolerance to structural distortions in ligand binding regions in protein models. For example, the percentage of correctly aligned pairs of adenine-binding sites in weakly homologous protein models is only 4–9% lower than those aligned using crystal structures. This represents a significant improvement over other algorithms, e.g. the performance of eMatchSite in recognizing similar binding sites is 6% and 13% higher than that of SiteEngine using high- and moderate-quality protein models, respectively. Constructing biologically correct alignments using predicted ligand binding sites in protein models opens up the possibility to investigate drug-protein interaction networks for complete proteomes with prospective systems-level applications in polypharmacology and rational drug repositioning. eMatchSite is freely available to the academic community as a web-server and a stand-alone software distribution at http://www.brylinski.org/ematchsite.

Independent of their localization in protein the hydrophobic amino acid residues have no effect on the molten globule state of apomyoglobin and the disulfide bond on the surface of apomyoglobin stabilizes this intermediate state

At present it is unclear which interactions in proteins reveal the presence of intermediate states, their stability and formation rate. In this study, we have investigated the effect of substitutions of hydrophobic amino acid residues in the hydrophobic core of protein and on its surface on a molten globule type intermediate state of apomyoglobin. It has been found that independent of their localization in protein, substitutions of hydrophobic amino acid residues do not affect the stability of the molten globule state of apomyoglobin. It has been shown also that introduction of a disulfide bond on the protein surface can stabilize the molten globule state. However in the case of apomyoglobin, stabilization of the intermediate state leads to relative destabilization of the native state of apomyoglobin. The result obtained allows us not only to conclude which mutations can have an effect on the intermediate state of the molten globule type, but also explains why the introduction of a disulfide bond (which seems to “strengthen” the protein) can result in destabilization of the protein native state of apomyoglobin.

Signature of hydrophobic hydration in a single polymer

Hydrophobicity underpins self-assembly in many natural and synthetic molecular and nanoscale systems. A signature of hydrophobicity is its temperature dependence. The first experimental evaluation of the temperature and size dependence of hydration free energy in a single hydrophobic polymer is reported, which tests key assumptions in models of hydrophobic interactions in protein folding. Herein, the hydration free energy required to extend three hydrophobic polymers with differently sized aromatic side chains was directly measured by single molecule force spectroscopy. The results are threefold. First, the hydration free energy per monomer is found to be strongly dependent on temperature and does not follow interfacial thermodynamics. Second, the temperature dependence profiles are distinct among the three hydrophobic polymers as a result of a hydrophobic size effect at the subnanometer scale. Third, the hydration free energy of a monomer on a macromolecule is different from a free monomer; corrections for the reduced hydration free energy due to hydrophobic interaction from neighboring units are required.

Microsecond scale replica exchange molecular dynamic simulation of villin headpiece: an insight into the folding landscape

Reaching the experimental time scale of millisecond is a grand challenge for protein folding simulations. The development of advanced Molecular Dynamics techniques like Replica Exchange Molecular Dynamics (REMD) makes it possible to reach these experimental timescales. In this study, an attempt has been made to reach the multi microsecond simulation time scale by carrying out folding simulations on a three helix bundle protein, Villin, by combining REMD and Amber United Atom model. Twenty replicas having different temperatures ranging from 295 K to 390 K were simulated for 1.5 µs each. The lowest Root Mean Square Deviation (RMSD) structure of 2.5 Å was obtained with respect to native structure (PDB code 1VII), with all the helices formed. The folding population landscapes were built using segment-wise RMSD and Principal Components as reaction coordinates. These analyses suggest the two-stage folding for Villin. The combination of REMD and Amber United Atom model may be useful to understand the folding mechanism of various fast folding proteins.

Intermediates in the protein folding process: a computational model

The paper presents a model for simulating the protein folding process in silico. The two-step model (which consists of the early stage—ES and the late stage—LS) is verified using two proteins, one of which is treated (according to experimental observations) as the early stage and the second as an example of the LS step. The early stage is based solely on backbone structural preferences, while the LS model takes into account the water environment, treated as an external hydrophobic force field and represented by a 3D Gauss function. The characteristics of 1ZTR (the ES intermediate, as compared with 1ENH, which is the LS intermediate) confirm the link between the gradual disappearance of ES characteristics in LS structural forms and the simultaneous emergence of LS properties in the 1ENH protein. Positive verification of ES and LS characteristics in these two proteins (1ZTR and 1ENH respectively) suggest potential applicability of the presented model to in silico protein folding simulations.

Free energy of contact formation in proteins: efficient computation in the elastic network approximation

Biomolecular simulations have become a major tool in understanding biomolecules and their complexes. However, one can typically only investigate a few mutants or scenarios due to the severe computational demands of such simulations, leading to a great interest in method development to overcome this restriction. One way to achieve this is to reduce the complexity of the systems by an approximation of the forces acting upon the constituents of the molecule. The harmonic approximation used in elastic network models simplifies the physical complexity to the most reduced dynamics of these molecular systems. The reduced polymer modeled this way is typically comprised of mass points representing coarse-grained versions of, e.g., amino acids. In this work, we show how the computation of free energy contributions of contacts between two residues within the molecule can be reduced to a simple lookup operation in a precomputable matrix. Being able to compute such contributions is of great importance: protein design or molecular evolution changes introduce perturbations to these pair interactions, so we need to understand their impact. Perturbation to the interactions occurs due to randomized and fixated changes (in molecular evolution) or designed modifications of the protein structures (in bioengineering). These perturbations are modifications in the topology and the strength of the interactions modeled by the elastic network models. We apply the new algorithm to (1) the bovine trypsin inhibitor, a well-known enzyme in biomedicine, and show the connection to folding properties and the hydrophobic collapse hypothesis and (2) the serine proteinase inhibitor CI-2 and show the correlation to Φ values to characterize folding importance. Furthermore, we discuss the computational complexity and show empirical results for the average case, sampled over a library of 77 structurally diverse proteins. We found a relative speedup of up to 10 000-fold for large proteins with respect to repeated application of the initial model.

Two-intermediate model to characterize the structure of fast-folding proteins

This paper introduces a new model that enables researchers to conduct protein folding simulations. A two-step in silico process is used in the course of structural analysis of a set of fast-folding proteins. The model assumes an early stage (ES) that depends solely on the backbone conformation, as described by its geometrical properties--specifically, by the V-angle between two sequential peptide bond planes (which determines the radius of curvature, also called R-radius, according to a second-degree polynomial form). The agreement between the structure under consideration and the assumed model is measured in terms of the magnitude of dispersion of both parameters with respect to idealized values. The second step, called late-stage folding (LS), is based on the "fuzzy oil drop" model, which involves an external hydrophobic force field described by a three-dimensional Gauss function. The degree of conformance between the structure under consideration and its idealized model is expressed quantitatively by means of the Kullback-Leibler entropy, which is a measure of disparity between the observed and expected hydrophobicity distributions. A set of proteins, representative of the fast-folding group - specifically, cold shock proteins - is shown to agree with the proposed model.

"Fuzzy oil drop" model applied to individual small proteins built of 70 amino acids

The proteins composed of short polypeptides (about 70 amino acid residues) representing the following functional groups (according to PDB notation): growth hormones, serine protease inhibitors, antifreeze proteins, chaperones and proteins of unknown function, were selected for structural and functional analysis. Classification based on the distribution of hydrophobicity in terms of deficiency/excess as the measure of structural and functional specificity is presented. The experimentally observed distribution of hydrophobicity in the protein body is compared to the idealized one expressed by a three-dimensional Gauss function. The differences between these two distributions reveal the specificity of structural/functional characteristics of the protein. The residues of hydrophobicity deficiency versus the idealized distribution are assumed to indicate cavities with the potential to bind ligands, while the residues of hydrophobicity excess are interpreted as potentially participating in protein-protein complexation. The distribution of hydrophobicity irregularity seems to be specific for particular structures and functions of proteins. A comparative analysis of such profiles is carried out to identify the potential biological activity of proteins of unknown function.

Application of the fuzzy-oil-drop model to membrane protein simulation

The analysis of structural properties and biological activity of membrane proteins requires long lasting simulation of molecular dynamics. The large number of atoms present in protein molecule, membrane (phospholipids), and water environment makes the simulation of large scale. The implementation of simplified model representing the natural environment for membrane proteins is presented and compared with the vacuum simulation and simulation in the presence of water molecules and membrane phospholipids presented explicite. The comparative structural analysis and computational times for these three models makes the simplified model promising.

On transversal hydrophobicity of some proteins and their modules

Hydrophobicity of proteins encoded in the genomes of diverse organisms was quantified using two novel concepts: (A) amino acid (AA) bulkiness-dependent hydrophobicity profiles and (B) spatial context of hydrophobicity distribution in AA triads. Both concepts were introduced into an algorithm that was used for extracting protein clusters from diverse genomic databases whose sequence attributes were similar to those in the multiple sequence alignment (MSA) of a given family of proteins. The sequences of the G protein-coupled receptors (GPCRs) encoded in different genomes were used as templates for testing the above concepts. The following sequence attributes were used for protein clustering: (A) sequence similarity scores (IDs); (B) amino acid composition (AAC); (C) hydrophobicity; (D) AA-bulkiness; and (E) alpha-helical propensity potentials. Diverse GPCRs display variable distributions of AA bulkiness-dependent buildups and declines in the hydrophobicity profiles that may be related to their function-dependent way of packing and allostery in the membrane. It is shown that intramolecular transversal nonbonded interactions between the TM segments in diverse GPCRs involve about 50% of hydrophobic atoms. Similar interaction networks exist between alpha-helices of tetratricopeptide (TPR) motifs-containing immunophilins and other proteins containing alpha-helical bundles.

Chaperonin structure: the large multi-subunit protein complex

The multi sub-unit protein structure representing the chaperonins group is analyzed with respect to its hydrophobicity distribution. The proteins of this group assist protein folding supported by ATP. The specific axial symmetry GroEL structure (two rings of seven units stacked back to back - 524 aa each) and the GroES (single ring of seven units - 97 aa each) polypeptide chains are analyzed using the hydrophobicity distribution expressed as excess/deficiency all over the molecule to search for structure-to-function relationships. The empirically observed distribution of hydrophobic residues is confronted with the theoretical one representing the idealized hydrophobic core with hydrophilic residues exposure on the surface. The observed discrepancy between these two distributions seems to be aim-oriented, determining the structure-to-function relation. The hydrophobic force field structure generated by the chaperonin capsule is presented. Its possible influence on substrate folding is suggested.

Never born proteins as a test case for ab initio protein structures prediction

The number of natural proteins although large is significantly smaller than the theoretical number of proteins that can be obtained combining the 20 natural amino acids, the so-called “never born proteins” (NBPs). The study of the structure and properties of these proteins allows to investigate the sources of the natural proteins being of unique characteristics or special properties. However the structural study of NPBs can also been intended as an ideal test for evaluating the efficiency of software packages for the ab initio protein structure prediction. In this research, 10.000 three-dimensional structures of proteins of completely random sequence generated according to ROSETTA and FOD model were compared. The results show the limits of these software packages, but at the same time indicate that in many cases there is a significant agreement between the prediction obtained.

Revised charge equilibration potential for liquid alkanes

We present a revised liquid alkane force field based on the charge equilibration formalism for incorporating electrostatic nonadditive effects arising from local polarization. The model is a revision of earlier work by Patel and Brooks, specifically addressing deficiencies in the dihedral potential, electrostatic, and Lennard-Jones (van der Waals) parameters of the force field. We discuss refinement of the alkane backbone torsion potential to match high-level ab initio relative conformational energetics for pentane, hexane, and heptane. We further discuss refinement of the electrostatic and Lennard-Jones (van der Waals) parameters to reproduce the experimental polarizability, liquid density, and vaporization enthalpy of hexane. Finally, we calculate bulk liquid properties including densities, vaporization enthalpies, self-diffusion constants, isothermal compressibilities, constant pressure heat capacities, and NMR T 1 relaxation times for a series of linear alkanes ranging from hexane to pentadecane based on the current revised model. We also compute free energies of hydration for pentane, hexane, and heptane. The revised force field offers a significantly improved overall description of these properties relative to the original parametrization. The current alkane force field represents a platform for ongoing development of a CHARMM (Chemistry at Harvard Molecular Mechanics) polarizable force field for lipids and integral membrane proteins.

Prediction of functional sites based on the fuzzy oil drop model

A description of many biological processes requires knowledge of the 3-D structure of proteins and, in particular, the defined active site responsible for biological function. Many proteins, the genes of which have been identified as the result of human genome sequencing, and which were synthesized experimentally, await identification of their biological activity. Currently used methods do not always yield satisfactory results, and new algorithms need to be developed to recognize the localization of active sites in proteins. This paper describes a computational model that can be used to identify potential areas that are able to interact with other molecules (ligands, substrates, inhibitors, etc.). The model for active site recognition is based on the analysis of hydrophobicity distribution in protein molecules. It is shown, based on the analyses of proteins with known biological activity and of proteins of unknown function, that the region of significantly irregular hydrophobicity distribution in proteins appears to be function related.

Localization of ligand binding site in proteins identified in silico

Knowledge-based models for protein folding assume that the early-stage structural form of a polypeptide is determined by the backbone conformation, followed by hydrophobic collapse. Side chain-side chain interactions, mostly of hydrophobic character, lead to the formation of the hydrophobic core, which seems to stabilize the structure of the protein in its natural environment. The fuzzy-oil-drop model is employed to represent the idealized hydrophobicity distribution in the protein molecule. Comparing it with the one empirically observed in the protein molecule reveals that they are not in agreement. It is shown in this study that the irregularity of hydrophobic distributions is aim-oriented. The character and strength of these irregularities in the organization of the hydrophobic core point to the specificity of a particular protein's structure/function. When the location of these irregularities is determined versus the idealized fuzzy-oil-drop, function-related areas in the protein molecule can be identified. The presented model can also be used to identify ways in which protein-protein complexes can possibly be created. Active sites can be predicted for any protein structure according to the presented model with the free prediction server at http://www.bioinformatics.cm-uj.krakow.pl/activesite. The implication based on the model presented in this work suggests the necessity of active presence of ligand during the protein folding process simulation.