Hydration structure and stability of Met-enkephalin studied by a three-dimensional reference interaction site model with a repulsive bridge correction and a thermodynamic perturbation method

Effect of caffeine on the aggregation of amyloid-β-A 3D RISM study

Alzheimer's disease is a detrimental neurological disorder caused by the formation of amyloid fibrils due to the aggregation of amyloid-β peptide. The primary therapeutic approaches for treating Alzheimer's disease are targeted to prevent this amyloid fibril formation using potential inhibitor molecules. The discovery of such inhibitor molecules poses a formidable challenge to the design of anti-amyloid drugs. This study investigates the effect of caffeine on dimer formation of the full-length amyloid-β using a combined approach of all-atom, explicit water molecular dynamics simulations and the three-dimensional reference interaction site model theory. The change in the hydration free energy of amyloid-β dimer, with and without the inhibitor molecules, is calculated with respect to the monomeric amyloid-β, where the hydration free energy is decomposed into energetic and entropic components, respectively. Dimerization is accompanied by a positive change in the partial molar volume. Dimer formation is spontaneous, which implies a decrease in the hydration free energy. However, a reverse trend is observed for the dimer with inhibitor molecules. It is observed that the negatively charged residues primarily contribute for the formation of the amyloid-β dimer. A residue-wise decomposition reveals that hydration/dehydration of the side-chain atoms of the charged amino acid residues primarily contribute to dimerization.

Aggregation propensities of proteins with varying degrees of disorder

The hydration thermodynamics of a globular protein (AcP), three intrinsically disordered protein regions (1CD3, 1MVF, 1F0R) and a fully disordered protein (α-synuclein) is studied by an approach that combines an all-atom explicit water molecular dynamics simulations and three-dimensional reference interaction site model (3D-RISM) theory. The variation in hydration free energy with percentage disorder of the selected proteins is investigated through its nonelectrostatic and electrostatic components. A decrease in hydration free energy is observed with an increase in percentage disorder, indicating favorable interactions of the disordered proteins with the solvent. This confirms the role of percentage disorder in determining the aggregation propensity of proteins which is measured in terms of the hydration free energy in addition to their respective mean net charge and mean hydrophobicity. The hydration free energy is decoupled into energetic and entropic terms. A residue-wise decomposition analysis of the hydration free energy for the selected proteins is evaluated. The decomposition shows that the disordered regions contribute more than the ordered ones for the intrinsically disordered protein regions. The dominant role of electrostatic interactions is confirmed from the residue-wise decomposition of the hydration free energy. The results depict that the negatively charged residues contribute more to the total hydration free energy for the proteins with negative mean net charge, while the positively charged residues contribute more for proteins with positive mean net charge.

Analysis of molecular dynamics simulations of 10-residue peptide, chignolin, using statistical mechanics: Relaxation mode analysis and three-dimensional reference interaction site model theory

Molecular dynamics simulation is a fruitful tool for investigating the structural stability, dynamics, and functions of biopolymers at an atomic level. In recent years, simulations can be performed on time scales of the order of milliseconds using special purpose systems. Since the most stable structure, as well as meta-stable structures and intermediate structures, is included in trajectories in long simulations, it is necessary to develop analysis methods for extracting them from trajectories of simulations. For these structures, methods for evaluating the stabilities, including the solvent effect, are also needed. We have developed relaxation mode analysis to investigate dynamics and kinetics of simulations based on statistical mechanics. We have also applied the three-dimensional reference interaction site model theory to investigate stabilities with solvent effects. In this paper, we review the results for designing amino-acid substitution of the 10-residue peptide, chignolin, to stabilize the misfolded structure using these developed analysis methods.

Theoretical approaches for dynamical ordering of biomolecular systems

BACKGROUND

Living systems are characterized by the dynamic assembly and disassembly of biomolecules. The dynamical ordering mechanism of these biomolecules has been investigated both experimentally and theoretically. The main theoretical approaches include quantum mechanical (QM) calculation, all-atom (AA) modeling, and coarse-grained (CG) modeling. The selected approach depends on the size of the target system (which differs among electrons, atoms, molecules, and molecular assemblies). These hierarchal approaches can be combined with molecular dynamics (MD) simulation and/or integral equation theories for liquids, which cover all size hierarchies.

SCOPE OF REVIEW

We review the framework of quantum mechanical/molecular mechanical (QM/MM) calculations, AA MD simulations, CG modeling, and integral equation theories. Applications of these methods to the dynamical ordering of biomolecular systems are also exemplified.

MAJOR CONCLUSIONS

The QM/MM calculation enables the study of chemical reactions. The AA MD simulation, which omits the QM calculation, can follow longer time-scale phenomena. By reducing the number of degrees of freedom and the computational cost, CG modeling can follow much longer time-scale phenomena than AA modeling. Integral equation theories for liquids elucidate the liquid structure, for example, whether the liquid follows a radial distribution function.

GENERAL SIGNIFICANCE

These theoretical approaches can analyze the dynamic behaviors of biomolecular systems. They also provide useful tools for exploring the dynamic ordering systems of biomolecules, such as self-assembly. This article is part of a Special Issue entitled "Biophysical Exploration of Dynamical Ordering of Biomolecular Systems" edited by Dr. Koichi Kato.

Integral equations in the study of polar and ionic interaction site fluids

In this review article we consider some of the current integral equation approaches and application to model polar liquid mixtures. We consider the use of multidimensional integral equations and in particular progress on the theory and applications of three dimensional integral equations. The IEs we consider may be derived from equilibrium statistical mechanical expressions incorporating a classical Hamiltonian description of the system. We give example including salt solutions, inhomogeneous solutions and systems including proteins and nucleic acids.

Protein solvation from theory and simulation: Exact treatment of Coulomb interactions in three-dimensional theories

Solvation forces dominate protein structure and dynamics. Integral equation theories allow a rapid and accurate evaluation of the effect of solvent around a complex solute, without the sampling issues associated with simulations of explicit solvent molecules. Advances in integral equation theories make it possible to calculate the angle dependent average solvent structure around an irregular molecular solution. We consider two methodological problems here: the treatment of long-ranged forces without the use of artificial periodicity or truncations and the effect of closures. We derive a method for calculating the long-ranged Coulomb interaction contributions to three-dimensional distribution functions involving only a rewriting of the system of integral equations and introducing no new formal approximations. We show the comparison of the exact forms with those implied by the supercell method. The supercell method is shown to be a good approximation for neutral solutes whereas the new method does not exhibit the known problems of the supercell method for charged solutes. Our method appears more numerically stable with respect to thermodynamic starting state. We also compare closures including the Kovalenko-Hirata closure, the hypernetted-chain (HNC) and an approximate three-dimensional bridge function combined with the HNC closure. Comparisons to molecular dynamics results are made for water as well as for the protein solute bovine pancreatic trypsin inhibitor. The proposed equations have less severe approximations and often provide results which compare favorably to molecular dynamics simulation where other methods fail.

Morphometric approach to the solvation free energy of complex molecules

We show that the solvation free energy of a complex molecule such as a protein can be calculated using only four geometrical measures of the molecular structure and corresponding thermodynamical coefficients. We compare results from this morphometric approach to those obtained by an elaborate statistical-mechanical theory in liquid state physics for a large variety of different structures of protein G and find excellent agreement. Since the computational time is drastically reduced, the new approach provides a practical and efficient way for calculating the solvation free energy which can be employed when this quantity has to be calculated for a large number of structures, as in a simulation study of protein folding.

Crucial importance of translational entropy of water in pressure denaturation of proteins

We present statistical thermodynamics of pressure denaturation of proteins, in which the three-dimensional integral equation theory is employed. It is applied to a simple model system focusing on the translational entropy of the solvent. The partial molar volume governing the pressure dependence of the structural stability of a protein is expressed for each structure in terms of the excluded volume for the solvent molecules, the solvent-accessible surface area (ASA), and a parameter related to the solvent-density profile formed near the protein surface. It is argued that the entropic effect originating from the translational movement of water molecules plays critical roles in the pressure-induced denaturation. We also show that the exceptionally small size of water molecules among dense liquids in nature is crucial for pressure denaturation. An unfolded structure, which is only moderately less compact than the native structure but has much larger ASA, is shown to turn more stable than the native one at an elevated pressure. The water entropy for the native structure is higher than that for the unfolded structure in the low-pressure region, whereas the opposite is true in the high-pressure region. Such a structure is characterized by the cleft and/or swelling and the water penetration into the interior. In another solvent whose molecular size is 1.5 times larger than that of water, however, the inversion of the stability does not occur any longer. The random coil becomes relatively more destabilized with rising pressure, irrespective of the molecular size of the solvent. These theoretical predictions are in qualitatively good agreement with the experimental observations.

Translational-entropy gain of solvent upon protein folding

We show that even in the complete absence of potential energies among the atoms in a protein-aqueous solution system, there is a physical factor that favors the folded state of the protein. It is a gain in the translational entropy (TE) of water originating from the translational movement of water molecules. An elaborate statistical-mechanical theory is employed to analyze the TE of water in which a protein or peptide with a prescribed conformation is immersed. It is shown that if the number of residues is sufficiently large, the TE gain is powerful enough to compete with the conformational-entropy loss upon folding. For protein G we have tested over 100 compact conformations generated by a computer simulation with the all-atom potentials as well as the native structure. A significant finding is that the largest TE is attained in the native structure. The translational movement of water molecules is quite effective in achieving the tight packing in the interior of a natural protein. These results are true only when the solvent is water whose molecular size is the smallest among the ordinary liquids in nature.