How do water-mediated interactions and osmotic second virial coefficients vary with particle size?

We examine quantitatively the solute-size dependences of the effective interactions between nonpolar solutes in water and in a simple liquid. The potential w(r) of mean force and the osmotic second virial coefficients B are calculated with high accuracy from molecular dynamics simulations. As the solute diameter increases from methane's to C60's with the solute-solute and solute-solvent attractive interaction parameters fixed to those for the methane-methane and methane-water interactions, the first minimum of w(r) lowers from -1.1 to -4.7 in units of the thermal energy kT. Correspondingly, the magnitude of B (<0) increases proportional to σα with some power close to 6 or 7, which reinforces the solute-size dependence of B found earlier for a smaller range of σ [H. Naito, R. Okamoto, T. Sumi and K. Koga, J. Chem. Phys., 2022, 156, 221104]. We also demonstrate that the strength of the attractive interactions between solute and solvent molecules can qualitatively change the characteristics of the effective pair interaction between solute particles, both in water and in a simple liquid. If the solute-solvent attractive force is set to be weaker (stronger) than a threshold, the effective interaction becomes increasingly attractive (repulsive) with increasing solute size.

Solvation free energies of alcohols in water: temperature and pressure dependences.

Solvation free energies µ* of amphiphilic species, methanol and 1,2-hexanediol, are obtained as a function of temperature or pressure based on molecular dynamics simulations combined with efficient free-energy calculation methods. In general, µ* of an amphiphile can be divided into and , the nonpolar and electrostatic contributions, and the former is further divided into and which are the work of cavity formation process and the free energy change due to weak, attractive interactions between the solute molecule and surrounding solvent molecules. We demonstrate that µ* of the two amphiphilic solutes can be obtained accurately using a perturbation combining method, which relies on the exact expressions for and and requires no simulations of intermediate systems between the solute with strong, repulsive interactions and the solute with the van der Waals and electrostatic interactions. The decomposition of µ* gives us several physical insights including that µ* is an increasing function of T due to , that the contributions of hydrophilic groups to the temperature dependence of µ* are additive, and that the contribution of the van der Waals attraction to the solvation volume is greater than that of the electrostatic interactions.

Molecular mechanism of the common and opposing cosolvent effects of fluorinated alcohol and urea on a coiled coil protein.

Alcohols and urea are widely used as effective protein denaturants. Among monohydric alcohols, 2,2,2-trifluoroethanol (TFE) has large cosolvent effects as a helix stabilizer in proteins. In contrast, urea efficiently denatures ordered native structures, including helices, into coils. These opposing cosolvent effects of TFE and urea are well known, even though both preferentially bind to proteins; however, the underlying molecular mechanism remains controversial. Cosolvent-dependent relative stability between native and denatured states is rigorously related to the difference in preferential binding parameters (PBPs) between these states. In this study, GCN4-p1 with two-stranded coiled coil helices was employed as a model protein, and molecular dynamics simulations for the helix dimer and isolated coil were conducted in aqueous solutions with 2 M TFE and urea. As 2 M cosolvent aqueous solutions did not exhibit clustering of cosolvent molecules, we were able to directly investigate the molecular origin of the excess PBP without considering the enhancement effect of PBPs arising from the concentration fluctuations. The calculated excess PBPs of TFE for the helices and those of urea for the coils were consistent with experimentally observed stabilization of helix by TFE and that of coil by urea. The former was caused by electrostatic interactions between TFE and side chains of the helices, while the latter was attributed to both electrostatic and dispersion interactions between urea and the main chains. Unexpectedly, reverse-micelle-like orientations of TFE molecules strengthened the electrostatic interactions between TFE and the side chains, resulting in strengthening of TFE solvation.

Muscarinic acetylcholine receptor-dependent and NMDA receptor-dependent LTP and LTD share the common AMPAR trafficking pathway.

The forebrain cholinergic system promotes higher brain function in part by signaling through the M1 muscarinic acetylcholine receptor (mAChR). Long-term potentiation (LTP) and long-term depression (LTD) of excitatory synaptic transmission in the hippocampus are also induced by mAChR. An AMPA receptor (AMPAR) trafficking model for hippocampal neurons has been proposed to simulate N-methyl-D-aspartate receptor (NMDAR)-dependent synaptic plasticity in the early phase. In this study, we demonstrated the validity of the hypothesis that the mAChR-dependent LTP/LTD shares a common AMPAR trafficking pathway associated with NMDAR-dependent LTP/LTD. However, unlike NMDAR, Ca2+ influx into the spine cytosol occurs owing to the Ca2+ stored inside the ER and is induced via the activation of inositol 1,4,5-trisphosphate (IP3) receptors during M1 mAChR activation. Moreover, the AMPAR trafficking model implies that alterations in LTP and LTD observed in Alzheimer's disease could be attributed to age-dependent reductions in AMPAR expression levels.

Immune response to SARS-CoV-2 in severe disease and long COVID-19.

COVID-19 is mild to moderate in otherwise healthy individuals but may nonetheless cause life-threatening disease and/or a wide range of persistent symptoms. The general determinant of disease severity is age mainly because the immune response declines in aging patients. Here, we developed a mathematical model of the immune response to SARS-CoV-2 and revealed that typical age-related risk factors such as only a several 10% decrease in innate immune cell activity and inhibition of type-I interferon signaling by autoantibodies drastically increased the viral load. It was reported that the numbers of certain dendritic cell subsets remained less than half those in healthy donors even seven months after infection. Hence, the inflammatory response was ongoing. Our model predicted the persistent DC reduction and showed that certain patients with severe and even mild symptoms could not effectively eliminate the virus and could potentially develop long COVID.

Improvement of Protein Solubility in Macromolecular Crowding during Myoglobin Evolution.

The inside of living cells is crowded by extremely high concentrations of biomolecules, and thus globular proteins should have been developed to increase their solubility under such crowding conditions during organic evolution. The O2-storage protein myoglobin (Mb) is known to be expressed in myocytes of diving mammals in much larger quantities than those of land mammals. We have previously resurrected ancient whale and pinniped Mbs and experimentally demonstrated that the diving animal Mbs have evolved to maintain high solubility under the crowding conditions or to increase their tolerance against macromolecular precipitants, rather than solubility in a dilute buffer solution. However, the detail of chemical mechanisms of the precipitant tolerance remains unclear. Here, we investigated pH dependence of the precipitant tolerance (ß, slope of the solubility against precipitant concentration) of extant Mbs and plotted the ß values, as well as those of ancestral Mbs, against their surface net charges (ZMb). The results demonstrated that the precipitant tolerance was approximated by the square of ZMb, that is, ß = aZMb2 + b, in which a and b are constants. This effect of ZMb against the precipitation is not predicted by a classical excluded volume theory that gives constant ß for Mbs but can be explained by electrostatic repulsion between Mb molecules. The present study elucidates how Mb molecules have evolved to increase their in vivo solubility and shows the physiological significance of either neutral or basic isoelectric points (pI) of the natural Mbs, rather than acidic pI.

Osmotic second virial coefficients for hydrophobic interactions as a function of solute size.

To gain quantitative insight into how the overall strength of the hydrophobic interaction varies with the molecular size, we calculate osmotic second virial coefficients B for hydrophobic spherical molecules of different diameters σ in water based on molecular simulation with corrections to the finite-size and finite-concentration effects. It is shown that B (<0) changes by two orders of magnitude greater as σ increases twofold and its solute-size dependence is best fit by a power law B ∝ σα with the exponent α ≃ 6, which contrasts with the cubic power law that the second virial coefficients of gases obey. It is also found that values of B for the solutes in a nonpolar solvent are positive but they obey the same power law as in water. A thermodynamic identity for B derived earlier [K. Koga, V. Holten, and B. Widom, J. Phys. Chem. B 119, 13391 (2015)] indicates that if B is asymptotically proportional to a power of σ, the exponent α must be equal to or greater than 6.

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.

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.

Ion Size Dependences of the Salting-Out Effect: Reversed Order of Sodium and Lithium Ions.

A general trend of the salting-out effect on hydrophobic solutes in aqueous solution is that the smaller the size of a dissolved ion, the larger the effect of reducing the solubility of a hydrophobe. An exception is that Li+, the smallest in alkali metal ions, has a notably weaker effect than Na+. To understand the reversed order in the cation series, we performed molecular dynamics simulations of aqueous solutions of salt ions and calculated the Setschenow coefficient of methane with the ionic radius of either a cation or an anion varied in a wide range. It is confirmed that the Setschenow coefficient is correlated with the packing fraction of salt solution, as observed in earlier studies, and also correlated with the partial molar volume of an ion. Analyses of correlation function integrals, packing fractions of solvation spheres, and orientations of water molecules surrounding an ion reveal the key differences in microscopic properties between the cation and anion series, which give rise to the reversed order in the cation series of the partial molar volumes of ions and ultimately that of the Setschenow coefficients.

Unveiling the Interaction Potential Surface between Drug-Entrapped Polymeric Micelles Clarifying the High Drug Nanocarrier Efficiency.

Polymeric micelles are invaluable media as drug nanocarriers. Although knowledge of an interaction between the micelles is a key to understanding the mechanisms and developing the superior functions, the interaction potential surface between drug-incorporated polymeric micelles has not yet been quantitatively evaluated due to the extremely complex structure. Here, the interaction potential surface between drug-entrapped polymeric micelles was unveiled by combining a small-angle scattering experiment and a model-potential-free liquid-state theory. Triblock copolymer composed of poly(ethylene oxide) and poly(propylene oxide) was investigated over a wide concentration range (0.5-10.0 wt %). Effects of the entrapment of a water-insoluble hydrophobic drug, cyclosporin A, on the interaction were explored by comparing the interactions with and without the drug. The results directly clarified the high drug carrier efficiency in terms of the interaction between the micelles. In addition, an investigation based on density functional theory provided a deeper insight into the monomer contribution to the extremely stable dispersion of the nanocarrier.

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.

Kinetics of the ancestral carbon metabolism pathways in deep-branching bacteria and archaea.

The origin of life is believed to be chemoautotrophic, deriving all biomass components from carbon dioxide, and all energy from inorganic redox couples in the environment. The reductive tricarboxylic acid cycle (rTCA) and the Wood-Ljungdahl pathway (WL) have been recognized as the most ancient carbon fixation pathways. The rTCA of the chemolithotrophic Thermosulfidibacter takaii, which was recently demonstrated to take place via an unexpected reverse reaction of citrate synthase, was reproduced using a kinetic network model, and a competition between reductive and oxidative fluxes on rTCA due to an acetyl coenzyme A (ACOA) influx upon acetate uptake was revealed. Avoiding ACOA direct influx into rTCA from WL is, therefore, raised as a kinetically necessary condition to maintain a complete rTCA. This hypothesis was confirmed for deep-branching bacteria and archaea, and explains the kinetic factors governing elementary processes in carbon metabolism evolution from the last universal common ancestor.

Mechanism underlying hippocampal long-term potentiation and depression based on competition between endocytosis and exocytosis of AMPA receptors.

N-methyl-D-aspartate (NMDA) receptor-dependent long-term potentiation (LTP) and long-term depression (LTD) of signal transmission form neural circuits and thus are thought to underlie learning and memory. These mechanisms are mediated by AMPA receptor (AMPAR) trafficking in postsynaptic neurons. However, the regulatory mechanism of bidirectional plasticity at excitatory synapses remains unclear. We present a network model of AMPAR trafficking for adult hippocampal pyramidal neurons, which reproduces both LTP and LTD. We show that the induction of both LTP and LTD is regulated by the competition between exocytosis and endocytosis of AMPARs, which are mediated by the calcium-sensors synaptotagmin 1/7 (Syt1/7) and protein interacting with C-kinase 1 (PICK1), respectively. Our result indicates that recycling endosomes containing AMPAR are always ready for Syt1/7-dependent exocytosis of AMPAR at peri-synaptic/synaptic membranes. This is because molecular motor myosin Vb constitutively transports the recycling endosome toward the membrane in a Ca2+-independent manner.

Interaction potential surface between Raman scattering enhancing nanoparticles conjugated with a functional copolymer.

A novel Raman scattering enhancement was discovered using colloid nanoparticles conjugated with an amine-based copolymer. The interaction potential surface between Raman scattering enhancing nanoparticles was clarified by combining a small-angle scattering method and a model-potential-free liquid-state theory as an in situ observation in the solution state. The potential surface indicates that the most stable position is located around 0.9 nm from the particle surface, suggesting the existence of a nanogap structure between the nanocomposites. The change in Raman scattering enhancement was also acquired during the dispersion process of the aggregated nanocomposites through a glutathione-triggered nanosensing reaction.

Is F1-ATPase a Rotary Motor with Nearly 100% Efficiency? Quantitative Analysis of Chemomechanical Coupling and Mechanical Slip.

We present a chemomechanical network model of the rotary molecular motor F1-ATPase which quantitatively describes not only the rotary motor dynamics driven by ATP hydrolysis but also the ATP synthesis caused by forced reverse rotations. We observe a high reversibility of F1-ATPase, that is, the main cycle of ATP synthesis corresponds to the reversal of the main cycle in the hydrolysis-driven motor rotation. However, our quantitative analysis indicates that torque-induced mechanical slip without chemomechanical coupling occurs under high external torque and reduces the maximal efficiency of the free energy transduction to 40-80% below the optimal efficiency. Heat irreversibly dissipates not only through the viscous friction of the probe but also directly from the motor due to torque-induced mechanical slip. Such irreversible heat dissipation is a crucial limitation for achieving a 100% free-energy transduction efficiency with biological nanomachines because biomolecules are easily deformed by external torque.

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.

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

Recently, we proposed a reference-modified density functional theory (RMDFT) to calculate solvation free energy (SFE), in which a hard-sphere fluid was introduced as the reference system instead of an ideal molecular gas. Through the RMDFT, using an optimal diameter for the hard-sphere reference system, the values of the SFE calculated at room temperature and normal pressure were in good agreement with those for more than 500 small organic molecules in water as determined by experiments. In this study, we present an application of the RMDFT for calculating the temperature and pressure dependences of the SFE for solute molecules in water. We demonstrate that the RMDFT has high predictive ability for the temperature and pressure dependences of the SFE for small solute molecules in water when the optimal reference hard-sphere diameter determined for each thermodynamic condition is used. We also apply the RMDFT to investigate the temperature and pressure dependences of the thermodynamic stability of an artificial small protein, chignolin, and discuss the mechanism of high-temperature and high-pressure unfolding of the protein. © 2017 Wiley Periodicals, Inc.

Myosin V: Chemomechanical-coupling ratchet with load-induced mechanical slip.

A chemomechanical-network model for myosin V is presented on the basis of both the nucleotide-dependent binding affinity of the head to an actin filament (AF) and asymmetries and similarity relations among the chemical transitions due to an intramolecular strain of the leading and trailing heads. The model allows for branched chemomechanical cycles and takes into account not only two different force-generating mechanical transitions between states wherein the leading head is strongly bound and the trailing head is weakly bound to the AF but also load-induced mechanical-slip transitions between states in which both heads are strongly bound. The latter is supported by the fact that ATP-independent high-speed backward stepping has been observed for myosin V, although such motility has never been for kinesin. The network model appears as follows: (1) the high chemomechanical-coupling ratio between forward step and ATP hydrolysis is achieved even at low ATP concentrations by the dual mechanical transitions; (2) the forward stepping at high ATP concentrations is explained by the front head-gating mechanism wherein the power stroke is triggered by the inorganic-phosphate (Pi) release from the leading head; (3) the ATP-binding or hydrolyzed ADP.Pi-binding leading head produces a stable binding to the AF, especially against backward loading.