Drive Protein Folding Research Articles

A Numerical Integrator for Kinetostatic Folding of Protein Molecules Modeled as Robots with Hyper Degrees of Freedom

The kinetostatic compliance method (KCM) models protein molecules as nanomechanisms consisting of numerous rigid peptide plane linkages. These linkages articulate with respect to each other through changes in the molecule dihedral angles, resulting in a kinematic mechanism with hyper degrees of freedom. Within the KCM framework, nonlinear interatomic forces drive protein folding by guiding the molecule’s dihedral angle vector towards its lowest energy state in a kinetostatic manner. This paper proposes a numerical integrator that is well suited to KCM-based protein folding and overcomes the limitations of traditional explicit Euler methods with fixed step size. Our proposed integration scheme is based on pseudo-transient continuation with an adaptive step size updating rule that can efficiently compute protein folding pathways, namely, the transient three-dimensional configurations of protein molecules during folding. Numerical simulations utilizing the KCM approach on protein backbones confirm the effectiveness of the proposed integrator.

Macromolecular condensation buffers intracellular water potential

Optimum protein function and biochemical activity critically depends on water availability because solvent thermodynamics drive protein folding and macromolecular interactions1. Reciprocally, macromolecules restrict the movement of ‘structured’ water molecules within their hydration layers, reducing the available ‘free’ bulk solvent and therefore the total thermodynamic potential energy of water, or water potential. Here, within concentrated macromolecular solutions such as the cytosol, we found that modest changes in temperature greatly affect the water potential, and are counteracted by opposing changes in osmotic strength. This duality of temperature and osmotic strength enables simple manipulations of solvent thermodynamics to prevent cell death after extreme cold or heat shock. Physiologically, cells must sustain their activity against fluctuating temperature, pressure and osmotic strength, which impact water availability within seconds. Yet, established mechanisms of water homeostasis act over much slower timescales2,3; we therefore postulated the existence of a rapid compensatory response. We find that this function is performed by water potential-driven changes in macromolecular assembly, particularly biomolecular condensation of intrinsically disordered proteins. The formation and dissolution of biomolecular condensates liberates and captures free water, respectively, quickly counteracting thermal or osmotic perturbations of water potential, which is consequently robustly buffered in the cytoplasm. Our results indicate that biomolecular condensation constitutes an intrinsic biophysical feedback response that rapidly compensates for intracellular osmotic and thermal fluctuations. We suggest that preserving water availability within the concentrated cytosol is an overlooked evolutionary driver of protein (dis)order and function.

Chaperonin: Co-chaperonin Interactions.

Co-chaperonins function together with chaperonins to mediate ATP-dependent protein folding in a variety of cellular compartments. Chaperonins are evolutionarily conserved and form two distinct classes, namely, group I and group II chaperonins. GroEL and its co-chaperonin GroES form part of group I and are the archetypal members of this family of protein folding machines. The unique mechanism used by GroEL and GroES to drive protein folding is embedded in the complex architecture of double-ringed complexes, forming two central chambers that undergo conformational rearrangements that enable protein folding to occur. GroES forms a lid over the chamber and in doing so dislodges bound substrate into the chamber, thereby allowing non-native proteins to fold in isolation. GroES also modulates allosteric transitions of GroEL. Group II chaperonins are functionally similar to group I chaperonins but differ in structure and do not require a co-chaperonin. A significant number of bacteria and eukaryotes house multiple chaperonin and co-chaperonin proteins, many of which have acquired additional intracellular and extracellular biological functions. In some instances, co-chaperonins display contrasting functions to those of chaperonins. Human HSP60 (HSPD) continues to play a key role in the pathogenesis of many human diseases, in particular autoimmune diseases and cancer. A greater understanding of the fascinating roles of both intracellular and extracellular Hsp10 on cellular processes will accelerate the development of techniques to treat diseases associated with the chaperonin family.

Collective Variable for Metadynamics Derived From AlphaFold Output.

AlphaFold is a neural network–based tool for the prediction of 3D structures of proteins. In CASP14, a blind structure prediction challenge, it performed significantly better than other competitors, making it the best available structure prediction tool. One of the outputs of AlphaFold is the probability profile of residue–residue distances. This makes it possible to score any conformation of the studied protein to express its compliance with the AlphaFold model. Here, we show how this score can be used to drive protein folding simulation by metadynamics and parallel tempering metadynamics. Using parallel tempering metadynamics, we simulated the folding of a mini-protein Trp-cage and β hairpin and predicted their folding equilibria. We observe the potential of the AlphaFold-based collective variable in applications beyond structure prediction, such as in structure refinement or prediction of the outcome of a mutation.

Crocin attenuates inflammation of lung tissue in ovalbumin-sensitized mice by altering the expression of endoplasmic reticulum stress markers.

Endoplasmic reticulum (ER) stress plays a pivotal role in the pathogenesis of asthma. The present study aimed to investigate the reducing or suppressing effects of crocin in ovalbumin (OVA)-sensitized mice on ER stress markers. Mice were divided into six groups (n=5 per group) including control, OVA-sensitized (OVA), OVA-treated crocin (OVA-Cr25, OVA-Cr50, and OVA-Cr100 mg/kg), and OVA-treated dexamethasone (1mg/kg), (OVA-Dexa) groups. Animals 5 later groups were sensitized to OVA and the treatment groups received intraperitoneally crocin/dexamethasone in the last 5 days of the model. At the end of the study, lung tissue was evaluated for airway inflammation, caspase 12 and CHOP protein levels, and expression of ER stress markers using real-time-PCR. Sensitization with OVA significantly caused airway inflammation and induction of ER stress in mice compared to the control group based on the elevated inflammatory cells and ER stress markers in the lung tissue. Treatment with crocin and dexamethasone reduced airway inflammation and suppressed ER stress markers. Interestingly, in the OVA-Cr100 group, the suppressive effects on ER stress apoptotic markers were comparable to the OVA-Dexa group. The results suggest that crocin mediates maladaptive ER stress conditions possibly by creating adaptive ER stress status and driving protein folding correctly.

Functional Interplay between P5 and PDI/ERp72 to Drive Protein Folding.

Simple SummaryThe physiological functions of proteins are destined by their unique three-dimensional structures. Almost all biological kingdoms share conserved disulfide-catalysts and chaperone networks that assist in correct protein folding and prevent aggregation. Disruption of these networks is implicated in pathogenesis, including neurodegenerative disease. In the mammalian endoplasmic reticulum (ER), more than 20 members of the protein disulfide isomerase family (PDIs) are believed to cooperate in the client folding pathway, but it remains unclear whether complex formation among PDIs via non-covalent interaction is involved in regulating their enzymatic and chaperone functions. Herein, we report novel functional hetero complexes between PDIs that promote oxidative folding and inhibit aggregation along client folding. The findings provide insight into the physiological significance of disulfide-catalyst and chaperone networks and clues for understanding pathogenesis associated with disruption of the networks.P5 is one of protein disulfide isomerase family proteins (PDIs) involved in endoplasmic reticulum (ER) protein quality control that assists oxidative folding, inhibits protein aggregation, and regulates the unfolded protein response. P5 reportedly interacts with other PDIs via intermolecular disulfide bonds in cultured cells, but it remains unclear whether complex formation between P5 and other PDIs is involved in regulating enzymatic and chaperone functions. Herein, we established the far-western blot method to detect non-covalent interactions between P5 and other PDIs and found that PDI and ERp72 are partner proteins of P5. The enzymatic activity of P5-mediated oxidative folding is up-regulated by PDI, while the chaperone activity of P5 is stimulated by ERp72. These findings shed light on the mechanism by which the complex formations among PDIs drive to synergistically accelerate protein folding and prevents aggregation. This knowledge has implications for understanding misfolding-related pathology.

Sibe: a computation tool to apply protein sequence statistics to predict folding and design in silico

BackgroundEvolutionary information contained in the amino acid sequences of proteins specifies the biological function and fold, but exactly what information contained in the protein sequence drives both of these processes? Considerable progress has been made to answer this fundamental question, but it remains challenging to explore the potential space of cooperative interactions between amino acids. Statistical analysis plays a significant role in studying such interactions and its use has expanded in recent years to studies ranging from coevolution-guided rational protein design to protein folding in silico.ResultsHere we describe a computational tool named Sibe for use in studies of protein sequence, folding, and design using evolutionary coupling between amino acids as a driving factor. In this study, Sibe is used to identify positionally conserved couplings between pairwise amino acids and aid rational protein design. In this process, pairwise couplings are filtered according to the relative entropy computed from the positional conservations and grouped into several ’blocks’, which could contribute to driving protein folding and design. A human β2-adrenergic receptor (β2AR) was used to demonstrate that those ’blocks’ contribute the rational design for specifying functional residues. Sibe also provides folding modules based on both the positionally conserved couplings and well-established statistical potentials for simulating protein folding in silico and predicting tertiary structure. Our results show that statistically inferences of basic evolutionary principles, such as conservations and coupled-mutations, can be used to rapidly design a diverse set of proteins and study protein folding.ConclusionsThe developed software Sibe provides a computational tool for systematical analysis from protein primary to its tertiary structure using the evolutionary couplings as a driving factor. Sibe, written in C++, accounts for compatibility with the ’big data’ era in biological science, and it primarily focuses on protein sequence analysis, but it is also applicable to extend to other modeling and predictions of experimental measurements

Structural properties and cellular expression of AfrLEA6, a group 6 late embryogenesis abundant protein from embryos of Artemia franciscana.

Late embryogenesis abundant (LEA) proteins are intrinsically disordered proteins (IDPs) commonly found in anhydrobiotic organisms and are frequently correlated with desiccation tolerance. Herein we report new findings on AfrLEA6, a novel group 6 LEA protein from embryos of Artemia franciscana. Assessment of secondary structure in aqueous and dried states with circular dichroism (CD) reveals 89% random coil in the aqueous state, thus supporting classification of AfrLEA6 as an IDP. Removal of water from the protein by drying or exposure to trifluoroethanol (a chemical de-solvating agent) promotes a large gain in secondary structure of AfrLEA6, predominated by α-helix and exhibiting minimal β-sheet structure. We evaluated the impact of physiological concentrations (up to 400mM) of the disaccharide trehalose on the folding of LEA proteins in solution. CD spectra for AfrLEA2, AfrLEA3m, and AfrLEA6 are unaffected by this organic solute noted for its ability to drive protein folding. AfrLEA6 exhibits its highest concentration in vivo during embryonic diapause, drops acutely at diapause termination, and then declines during development to undetectable values at the larval stage. Maximum cellular titer of AfrLEA6 was 10-fold lower than for AfrLEA2 or AfrLEA3, both group 3 LEA proteins. Acute termination of diapause with H2O2 (a far more effective terminator than desiccation in this Great Salt Lake, UT, population) fostered a rapid 38% decrease in AfrLEA6 content of embryos. While the ultimate mechanism of diapause termination is unknown, disruption of key macromolecules could initiate physiological signaling events necessary for resumption of development and metabolism.

How Hydrophobic Hydration Destabilizes Surfactant Micelles at Low Temperature: A Coarse-Grained Simulation Study.

Micelles are self-assembled aggregates of amphiphilic surfactant molecules that are important in a variety of applications, including drug delivery, detergency, and catalysis. It is known that the micellization process is driven by the same physiochemical forces that drive protein folding, aggregation, and biological membrane self-assembly. Nevertheless, the molecular details of how micelle stability changes in water at low temperature are not fully clear. We develop and use a coarse-grained model to investigate how the interplay between nonionic surfactants and the surrounding water at the nanoscale affects the stability of micelles at high and low temperatures. Simulations of preformed C12E5 micelles in explicit water at a range of temperatures reveal the existence of two distinct surfactant conformations within the micelle, a bent structure and an extended structure, the latter being more prevalent at low temperature. Favorable interactions of the surfactant with more ordered solvation water stabilizes the extended configuration, allowing nanoscale wetting of the dry, hydrophobic core of the micelle, leading to the micelle breaking. Taken together, our coarse-grained simulations unravel how energetic and structural changes of the surfactant and the surrounding water destabilize micelles at low temperature, which is a direct consequence of the weakened hydrophobicity. Our approach thus provides an effective mean for extracting the molecular-level changes during hydrophobicity-driven destabilization of molecular self-assembly, which is important in a wide range of fields, including biology, polymer science, and nanotechnology.

Hydrophobic Collapse in N-Methylacetamide-Water Mixtures.

Aqueous N-methylacetamide solutions were investigated by polarization-resolved pump–probe and 2D infrared spectroscopy (2D IR), using the amide I mode as a reporter. The 2D IR results are compared with molecular dynamics simulations and spectral calculations to gain insight into the molecular structures in the mixture. N-Methylacetamide and water molecules tend to form clusters with “frozen” amide I dynamics. This is driven by a hydrophobic collapse as the methyl groups of the N-methylacetamide molecules cluster in the presence of water. Since the studied system can be considered as a simplified model for the backbone of proteins, the present study forms a convenient basis for understanding the structural and vibrational dynamics in proteins. It is particularly interesting to find out that a hydrophobic collapse as the one driving protein folding is observed in such a simple system.

Mechanistic insights on the reduction of glutathione disulfide by protein disulfide isomerase

We explore the enzymatic mechanism of the reduction of glutathione disulfide (GSSG) by the reduced a domain of human protein disulfide isomerase (hPDI) with atomistic resolution. We use classical molecular dynamics and hybrid quantum mechanics/molecular mechanics calculations at the mPW1N/6-311+G(2d,2p):FF99SB//mPW1N/6-31G(d):FF99SB level. The reaction proceeds in two stages: (i) a thiol-disulfide exchange through nucleophilic attack of the Cys53-thiolate to the GSSG-disulfide followed by the deprotonation of Cys56-thiol by Glu47-carboxylate and (ii) a second thiol-disulfide exchange between the Cys56-thiolate and the mixed disulfide intermediate formed in the first step. The Gibbs activation energy for the first stage was 18.7 kcal·mol-1, and for the second stage, it was 7.2 kcal·mol-1, in excellent agreement with the experimental barrier (17.6 kcal·mol-1). Our results also suggest that the catalysis by protein disulfide isomerase (PDI) and thiol-disulfide exchange is mostly enthalpy-driven (entropy changes below 2 kcal·mol-1 at all stages of the reaction). Hydrogen bonds formed between the backbone of His55 and Cys56 and the Cys56-thiol result in an increase in the Gibbs energy barrier of the first thiol-disulfide exchange. The solvent plays a key role in stabilizing the leaving glutathione thiolate formed. This role is not exclusively electrostatic, because an explicit inclusion of several water molecules at the density-functional theory level is a requisite to form the mixed disulfide intermediate. In the intramolecular oxidation of PDI, a transition state is only observed if hydrogen bond donors are nearby the mixed disulfide intermediate, which emphasizes that the thermochemistry of thiol-disulfide exchange in PDI is influenced by the presence of hydrogen bond donors.

Molecular dynamics analysis of the aggregation propensity of polyglutamine segments.

Protein misfolding and aggregation is a pathogenic feature shared among at least ten polyglutamine (polyQ) neurodegenerative diseases. While solvent-solution interaction is a key factor driving protein folding and aggregation, the solvation properties of expanded polyQ tracts are not well understood. By using GPU-enabled all-atom molecular dynamics simulations of polyQ monomers in an explicit solvent environment, this study shows that solvent-polyQ interaction propensity decreases as the lengths of polyQ tract increases. This study finds a predominance in long-distance interactions between residues far apart in polyQ sequences with longer polyQ segments, that leads to significant conformational differences. This study also indicates that large loops, comprised of parallel β-structures, appear in long polyQ tracts and present new aggregation building blocks with aggregation driven by long-distance intra-polyQ interactions. Finally, consistent with previous observations using coarse-grain simulations, this study demonstrates that there is a gain in the aggregation propensity with increased polyQ length, and that this gain is correlated with decreasing ability of solvent-polyQ interaction. These results suggest the modulation of solvent-polyQ interactions as a possible therapeutic strategy for treating polyQ diseases.

> π* Interaction in Designed Peptides

An n‐‐>π* interaction is a non‐covalent interaction that orders two consecutive residues in a peptide/protein by delocalizing the lone pair of electrons from the oxygen of a donor carbonyl to the antibonding orbital of an acceptor carbonyl carbon. This interaction gives insight into the inherent biases that drive protein folding, as it stabilizes secondary structures such as α‐helices and polyproline helices.We synthesized a series of hydroxyproline‐based peptides with variable N‐terminal substituents, modulating the electron‐donating capabilities on the donor carbonyl. By changing the electronics of each donor carbonyl we can tune the strength of this n‐‐>π* interaction and explore its structural implications. The n‐‐>π* interaction stabilizes the trans conformation of a peptide bond. Because proline can readily populate both trans and cis conformations, the ratio of the two species (Ktrans/cis) indicates the strength of the n‐‐>π* interaction.The series of derivatives synthesized were selected to compare the electron‐donating capabilities of their donor carbonyl to its n‐‐>π* interaction strength. The derivatives synthesized for analysis were, X‐Hyp(4‐NO2Bz)‐OMe where X = Boc, pivaloyl, formyl, acetyl, propionyl, monochloroacetyl, monofluoroacetyl, dichloroacetyl, difluoroacetyl, trichloroacetyl and trifluoroacetyl. These derivatives provide a range of electronic properties to be analyzed for the strength of the nàπ* interaction in each compound.Combined analysis via 1H and 13C NMR and X‐ray crystallography reveals the characteristics of the n‐‐>π* interaction. Cis and trans species can be differentiated by 1H NMR, providing the Ktrans/cis for each derivative. Electron donation involved in the n‐‐>π* interaction is seen via the carbonyl chemical shifts using 13C NMR. The chemical shift of the acceptor carbonyl moves downfield with stronger electron‐donating substituents, while carbonyls uninvolved in the interaction remain constant. X‐ray crystallography reveals the distances and angles of the n‐‐>π* interaction, including torsion angles seen in protein secondary structures. Specifically, the propionyl donor carbonyl results in the crystallographic observation of an α‐helical conformation, inducing compact torsion angles (φ = −58.6° and ψ = −37.8°), without the presence of a hydrogen bond. An electron rich donor results in the strongest interaction. Our comparison of this series of donor‐acceptor carbonyl pairs reveals the energetics and geometry of the n‐‐>π* interaction, granting new and unique insights into the problem of protein folding.Support or Funding InformationI'd like to extend my gratitude towards David Plastino and the NSF for providing funding for this research.

Detection of Cysteine Redox States in Mitochondrial Proteins in Intact Mammalian Cells.

Import, folding, and activity regulation of mitochondrial proteins are important for mitochondrial function. Cysteine residues play crucial roles in these processes as their thiol groups can undergo (reversible) oxidation reactions. For example, during import of many intermembrane space (IMS) proteins, cysteine oxidation drives protein folding and translocation over the outer membrane. Mature mitochondrial proteins can undergo changes in the redox state of specific cysteine residues, for example, as part of their enzymatic reaction cycle or as adaptations to changes of the local redox environment which might influence their activity. Here we describe methods to study changes in cysteine residue redox states in intact cells. These approaches allow to monitor oxidation-driven protein import as well as changes of cysteine redox states in mature proteins during oxidative stress or during the reaction cycle of thiol-dependent enzymes like oxidoreductases.

Basis of Protein Stabilization by K Glutamate: Unfavorable Interactions with Carbon, Oxygen Groups

Potassium glutamate (KGlu) is the primary Escherichia coli cytoplasmic salt. After sudden osmotic upshift, cytoplasmic KGlu concentration increases, initially because of water efflux and subsequently by K+ transport and Glu− synthesis, allowing water uptake and resumption of growth at high osmolality. In vitro, KGlu ranks with Hofmeister salts KF and K2SO4 in driving protein folding and assembly. Replacement of KCl by KGlu stabilizes protein-nucleic acid complexes. To interpret and predict KGlu effects on protein processes, preferential interactions of KGlu with 15 model compounds displaying six protein functional groups—sp3 (aliphatic) C; sp2 (aromatic, amide, carboxylate) C; amide and anionic (carboxylate) O; and amide and cationic N—were determined by osmometry or solubility assays. Analysis of these data yields interaction potentials (α-values) quantifying non-Coulombic chemical interactions of KGlu with unit area of these six groups. Interactions of KGlu with the 15 model compounds predicted from these six α-values agree well with experimental data. KGlu interactions with all carbon groups and with anionic (carboxylate) and amide oxygen are unfavorable, while KGlu interactions with cationic and amide nitrogen are favorable. These α-values, together with surface area information, provide quantitative predictions of why KGlu is an effective E. coli cytoplasmic osmolyte (because of the dominant effect of unfavorable interactions of KGlu with anionic and amide oxygens and hydrocarbon groups on the water-accessible surface of cytoplasmic biopolymers) and why KGlu is a strong stabilizer of folded proteins (because of the dominant effect of unfavorable interactions of KGlu with hydrocarbon groups and amide oxygens exposed in unfolding).

How the hydrophobic factor drives protein folding

How hydrophobicity (HY) drives protein folding is studied. The 1971 Nozaki-Tanford method of measuring HY is modified to use gases as solutes, not crystals, and this makes the method easy to use. Alkanes are found to be much more hydrophobic than rare gases, and the two different kinds of HY are termed intrinsic (rare gases) and extrinsic (alkanes). The HY values of rare gases are proportional to solvent-accessible surface area (ASA), whereas the HY values of alkanes depend on special hydration shells. Earlier work showed that hydration shells produce the hydration energetics of alkanes. Evidence is given here that the transfer energetics of alkanes to cyclohexane [Wolfenden R, Lewis CA, Jr, Yuan Y, Carter CW, Jr (2015) Proc Natl Acad Sci USA 112(24):7484-7488] measure the release of these shells. Alkane shells are stabilized importantly by van der Waals interactions between alkane carbon and water oxygen atoms. Thus, rare gases cannot form this type of shell. The very short (approximately picoseconds) lifetime of the van der Waals interaction probably explains why NMR efforts to detect alkane hydration shells have failed. The close similarity between the sizes of the opposing energetics for forming or releasing alkane shells confirms the presence of these shells on alkanes and supports Kauzmann's 1959 mechanism of protein folding. A space-filling model is given for the hydration shells on linear alkanes. The model reproduces the n values of Jorgensen et al. [Jorgensen WL, Gao J, Ravimohan C (1985) J Phys Chem 89:3470-3473] for the number of waters in alkane hydration shells.

Chaperonin-co-chaperonin interactions.

Co-chaperonins function together with chaperonins to mediate ATP-dependant protein folding in a variety of cellular compartments. GroEL and its co-chaperonin GroES are the only essential chaperones in Escherichia coli and are the archetypal members of this family of protein folding machines. The unique mechanism used by GroEL and GroES to drive protein folding is embedded in the complex architecture of double-ringed complexes, forming two central chambers that undergo structural rearrangements as part of the folding mechanism. GroES forms a lid over the chamber, and in doing so dislodges bound substrate into the chamber, thereby allowing non-native proteins to fold in isolation. GroES also modulates allosteric transitions of GroEL. A significant number of bacteria and eukaryotes house multiple chaperonin and co-chaperonin proteins, many of which have acquired additional intracellular and extracellular biological functions. In some instances co-chaperonins display contrasting functions to those of chaperonins. Human Hsp60 continues to play a key role in the pathogenesis of many human diseases, in particular autoimmune diseases and cancer. A greater understanding of the fascinating roles of both intracellular and extracellular Hsp10, in addition to its role as a co-chaperonin, on cellular processes will accelerate the development of techniques to treat diseases associated with the chaperonin family.

Dynamic hydration shell restores Kauzmann's 1959 explanation of how the hydrophobic factor drives protein folding.

Kauzmann's explanation of how the hydrophobic factor drives protein folding is reexamined. His explanation said that hydrocarbon hydration shells are formed, possibly of clathrate water, and they explain why hydrocarbons have uniquely low solubilities in water. His explanation was not universally accepted because of skepticism about the clathrate hydration shell. A revised version is given here in which a dynamic hydration shell is formed by van der Waals (vdw) attraction, as proposed in 1985 by Jorgensen et al. [Jorgensen WL, Gao J, Ravimohan C (1985) J Phys Chem 89:3470-3473]. The vdw hydration shell is implicit in theories of hydrophobicity that contain the vdw interaction between hydrocarbon C and water O atoms. To test the vdw shell model against the known hydration energetics of alkanes, the energetics should be based on the Ben-Naim standard state (solute transfer between fixed positions in the gas and liquid phases). Then the energetics are proportional to n, the number of water molecules correlated with an alkane by vdw attraction, given by the simulations of Jorgensen et al. The energetics show that the decrease in entropy upon hydration is the root cause of hydrophobicity; it probably results from extensive ordering of water molecules in the vdw shell. The puzzle of how hydrophobic free energy can be proportional to nonpolar surface area when the free energy is unfavorable and the only known interaction (the vdw attraction) is favorable, is resolved by finding that the unfavorable free energy is produced by the vdw shell.

Water promotes the sealing of nanoscale packing defects in folding proteins

A net dipole moment is shown to arise from a non-Debye component of water polarization created by nanoscale packing defects on the protein surface. Accordingly, the protein electrostatic field exerts a torque on the induced dipole, locally impeding the nucleation of ice at the protein–water interface. We evaluate the solvent orientation steering (SOS) as the reversible work needed to align the induced dipoles with the Debye electrostatic field and computed the SOS for the variable interface of a folding protein. The minimization of the SOS is shown to drive protein folding as evidenced by the entrainment of the total free energy by the SOS energy along trajectories that approach a Debye limit state where no torque arises. This result suggests that the minimization of anomalous water polarization at the interface promotes the sealing of packing defects, thereby maintaining structural integrity and committing the protein chain to fold.

Distance-dependent hydrophobic-hydrophobic contacts in protein folding simulations.

Successful prediction of protein folding from an amino acid sequence is a challenge in computational biology. In order to reveal the geometric constraints that drive protein folding, highlight those constraints kept or missed by distinct lattices and for establishing which class of intra- and inter-secondary structure element interactions is the most relevant for the correct folding of proteins, we have calculated inter-alpha carbon distances in a set of 42 crystal structures consisting of mainly helix, sheet or mixed conformations. The inter-alpha carbon distances were also calculated in several lattice "hydrophobic-polar" models built from the same protein set. We found that helix structures are more prone to form "hydrophobic-hydrophobic" contacts than beta-sheet structures. At a distance lower than or equal to 3.8 Å (very short-range interactions), "hydrophobic-hydrophobic" contacts are almost absent in the native structures, while they are frequent in all the analyzed lattice models. At distances in-between 3.8 and 9.5 Å (short-/medium-range interactions), the best performing lattice for reproducing mainly helix structures is the body-centered-cubic lattice. If protein structures contain sheet portions, lattice performances get worse, with few exceptions observed for double-tetrahedral and body-centered-cubic lattices. Finally, we can observe that ab initio protein folding algorithms, i.e. those based on the employment of lattices and Monte Carlo simulated annealings, can be improved simply and effectively by preventing the generation of "hydrophobic-hydrophobic" contacts shorter than 3.8 Å, by monitoring the "hydrophobic-hydrophobic/polar-polar" contact ratio in short-/medium distance ranges and by using preferentially a body-centered-cubic lattice.