Application of reference‐modified density functional theory: Temperature and pressure dependences of solvation free energy

Structural Stability Analysis of Proteins Using End-to-End Distance: A 3D-RISM Approach

The stability of a protein is determined from its properties and surrounding solvent. In our previous study, the total energy as a sum of the conformational and solvation free energies was demonstrated to be an appropriate energy function for evaluating the stability of a protein in a protein folding system. We plotted the various energies against the root mean square deviation, required as a reference structure. Herein, we replotted the various energies against the end-to-end distance between the N- and C-termini, which is not a required reference and is experimentally measurable. The solvation free energies for all proteins tend to be low as the end-to-end distance increases, whereas the conformational energies tend to be low as the end-to-end distance decreases. The end-to-end distance is one of interesting measures to study the behavior of proteins.

Common and unique strategies of myoglobin evolution for deep-sea adaptation of diving mammals

Myoglobin (Mb) is highly concentrated in the myocytes of diving mammals such as whales and seals, in comparison with land animals, and its molecular evolution has played a crucial role in their deep-sea adaptation. We previously resurrected ancestral whale Mbs and demonstrated the evolutional strategies for higher solubility under macromolecular crowding conditions. Pinnipeds, such as seals and sea lions, are also expert diving mammals with Mb-rich muscles. In the present study, we resurrected ancestral pinniped Mbs and investigated their biochemical and structural properties. Comparisons between pinniped and whale Mbs revealed the common and distinctive strategies for the deep-sea adaptation. The overall evolution processes, gaining precipitant tolerance and improving thermodynamic stability, were commonly observed. However, the strategies for improving the folding stability differed, and the pinniped Mbs exploited the shielding of hydrophobic surfaces more effectively than the whale Mbs.

Water-mediated interactions destabilize proteins

Proteins are folded to avoid exposure of the nonpolar groups to water because water‐mediated interactions between nonpolar groups are a promising factor in the thermodynamic stabilities of proteins—which is a well‐accepted view as one of the unique effects of hydrophobic interactions. This article poses a critical question for this classical view by conducting an accurate solvation free‐energy calculation for a thermodynamic cycle of a protein folding using a liquid‐state density functional theory. Here, the solvation‐free energy for a leucine zipper formation was examined in the coiled‐coil protein GCN4‐p1, a typical model for hydrophobic interactions, which demonstrated that water‐mediated interactions were unfavorable for the association of nonpolar groups in the native state, while the dispersion forces between them were, instead, responsible for the association. Furthermore, the present analysis well predicted the isolated helical state stabilized by pressure, which was previously observed in an experiment. We reviewed the problems in the classical concept and semiempirical presumption that the energetic cost of the hydration of nonpolar groups is a driving force of folding.

Two different regimes in alcohol-induced coil-helix transition: effects of 2,2,2-trifluoroethanol on proteins being either independent of or enhanced by solvent structural fluctuations

Inhomogeneous distribution of constituent molecules in a mixed solvent has been known to give remarkable effects on the solute, e.g., conformational changes of biomolecules in an alcohol-water mixture. We investigated the general effects of 2,2,2-trifluoroethanol (TFE) on proteins/peptides in a mixture of water and TFE using melittin as a model protein. Fluctuations and Kirkwood-Buff integrals (KBIs) in the TFE-H2O mixture, quantitative descriptions of inhomogeneity, were determined by small-angle X-ray scattering investigation and compared with those in the aqueous solutions of other alcohols. The concentration fluctuation for the mixtures ranks as methanol < ethanol ≪ TFE < tert-butanol < 1-propanol, indicating that the inhomogeneity of molecular distribution in the TFE-H2O mixture is unexpectedly comparable to those in the series of mono-ols. On the basis of the concentration dependence of KBIs between the TFE molecules, it was found that a strong attraction between the TFE molecules is not necessarily important to induce helix conformation, which is inconsistent with the previously proposed mechanism. To address this issue, by combining the KBIs and the helix contents reported by the experimental spectroscopic studies, we quantitatively evaluated the change in the preferential binding parameter of TFE to melittin attributed to the coil-helix transition. As a result, we found two different regimes on TFE-induced helix formation. In the dilute concentration region of TFE below ∼2 M, where the TFE molecules are not aggregated among themselves, the excess preferential binding of TFE to the helix occurs due to the direct interaction between them, namely independent of the solvent fluctuation. In the higher concentration region above ∼2 M, in addition to the former effect, the excess preferential binding is significantly enhanced by the solvent fluctuation. This scheme should be held as general cosolvent effects of TFE on proteins/peptides.

Theoretical analysis on thermodynamic stability of chignolin

Understanding the dominant factor in thermodynamic stability of proteins remains an open challenge. Kauzmann’s hydrophobic interaction hypothesis, which considers hydrophobic interactions between nonpolar groups as the dominant factor, has been widely accepted for about sixty years and attracted many scientists. The hypothesis, however, has not been verified or disproved because it is difficult, both theoretically and experimentally, to quantify the solvent effects on the free energy change in protein folding. Here, we developed a computational method for extracting the dominant factor behind thermodynamic stability of proteins and applied it to a small, designed protein, chignolin. The resulting free energy profile quantitatively agreed with the molecular dynamics simulations. Decomposition of the free energy profile indicated that intramolecular interactions predominantly stabilized collapsed conformations, whereas solvent-induced interactions, including hydrophobic ones, destabilized them. These results obtained for chignolin were consistent with the site-directed mutagenesis and calorimetry experiments for globular proteins with hydrophobic interior cores.

Tracing whale myoglobin evolution by resurrecting ancient proteins

Extant cetaceans, such as sperm whale, acquired the great ability to dive into the ocean depths during the evolution from their terrestrial ancestor that lived about 50 million years ago. Myoglobin (Mb) is highly concentrated in the myocytes of diving animals, in comparison with those of land animals, and is thought to play a crucial role in their adaptation as the molecular aqualung. Here, we resurrected ancestral whale Mbs, which are from the common ancestor between toothed and baleen whales (Basilosaurus), and from a further common quadrupedal ancestor between whale and hippopotamus (Pakicetus). The experimental and theoretical analyses demonstrated that whale Mb adopted two distinguished strategies to increase the protein concentration in vivo along the evolutionary history of deep sea adaptation; gaining precipitant tolerance in the early phase of the evolution, and increase of folding stability in the late phase.