Structural and Thermodynamic Characteristics That Seed Aggregation of Amyloid-β Protein in Water

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.

Total free energy analysis of fully hydrated proteins

The total free energy of a hydrated biomolecule and its corresponding decomposition of energy and entropy provides detailed information about regions of thermodynamic stability or instability. The free energies of four hydrated globular proteins with different net charges are calculated from a molecular dynamics simulation, with the energy coming from the system Hamiltonian and entropy using multiscale cell correlation. Water is found to be most stable around anionic residues, intermediate around cationic and polar residues, and least stable near hydrophobic residues, especially when more buried, with stability displaying moderate entropy-enthalpy compensation. Conversely, anionic residues in the proteins are energetically destabilized relative to singly solvated amino acids, while trends for other residues are less clear-cut. Almost all residues lose intraresidue entropy when in the protein, enthalpy changes are negative on average but may be positive or negative, and the resulting overall stability is moderate for some proteins and negligible for others. The free energy of water around single amino acids is found to closely match existing hydrophobicity scales. Regarding the effect of secondary structure, water is slightly more stable around loops, of intermediate stability around β strands and turns, and least stable around helices. An interesting asymmetry observed is that cationic residues stabilize a residue when bonded to its N-terminal side but destabilize it when on the C-terminal side, with a weaker reversed trend for anionic residues.

Conformational Transitions of Amyloid-β: A Langevin and Generalized Langevin Dynamics Simulation Study

The dynamics of conformational transitions of the disordered protein, amyloid-β, is studied via Langevin and generalized Langevin dynamics simulations. The transmission coefficient for the unfold-misfold transition of amyloid-β is calculated from multiple independent trajectories that originate at the transition state with different initial velocities and are directly correlated to Kramers and Grote-Hynes theories. For lower values of the frictional coefficient, a well-defined rate constant is obtained, whereas, for higher values, the transmission coefficient decays with time, indicating a breakdown of the Kramers and Grote-Hynes theories and the emergence of a dynamic disorder, which demonstrates the presence of multiple local minima in the misfolding potential energy surface. The calculated free energy profile describes a two-state transition of amyloid-β in the energy landscape. The transition path time distribution computed from these simulations is compared with the related experimental and theoretical results for the unfold-misfold transition of amyloid-β. The high free energy barrier for this transition confirms the misfolding of amyloid-β. These findings offer an insight into the dynamics of the unfold-misfold transition of this protein.

Hydration Thermodynamics of Familial Parkinson's Disease-Linked Mutants of α-Synuclein

The hydration thermodynamics of different mutants of α-synuclein (α-syn) related to familial Parkinson's disease (PD) is explored using a computational approach that combines both molecular dynamics simulations in water and integral equation theory of molecular liquids. This analysis focuses on the change in conformational entropy, hydration free energy (HFE), and partial molar volume of α-syn upon mutation. The results show that A53T, A30P, E46K, and H50Q mutants aggregate more readily and display increased HFE and less negative interaction volume than the wild-type α-syn. In contrast, an opposite trend is observed for the G51D mutant with a lower experimental aggregation rate. The residuewise decomposition analysis of the HFE highlights that the dehydration/hydration of the hydrophilic residue-rich N- and C-termini of α-syn majorly contributes to the change upon mutation. The hydration shell contributions of different residues to the interaction volume are consistent with its increase/decrease upon mutation. This work shows that both HFE and interaction volume determine the aggregation kinetics of α-syn upon mutation and may serve as an appropriate benchmark for the treatment of PD.

Hydration Thermodynamics of the N-Terminal FAD Mutants of Amyloid-β

The hydration thermodynamics of amyloid-β (Aβ) and its pathogenic familial Alzheimer's disease (FAD) mutants such as A2V, Taiwan (D7H), Tottori (D7N), and English (H6R) and the protective A2T mutant is investigated by a combination of all-atom, explicit water molecular dynamics (MD) simulations and the three-dimensional reference interaction site model (3D-RISM) theory. The change in the hydration free energy on mutation is decomposed into the energetic and entropic components, which comprise electrostatic and nonelectrostatic contributions. An increase in the hydration free energy is observed for A2V, D7H, D7N, and H6R mutations that increase the aggregation propensity of Aβ and lead to an early onset of Alzheimer's disease, while a reverse trend is noted for the protective A2T mutation. An antiphase correlation is found between the change in the hydration energy and the internal energy of Aβ upon mutation. A residue-wise decomposition analysis shows that the change in the hydration free energy of Aβ on mutation is primarily due to the hydration/dehydration of the side-chain atoms of the negatively charged residues. The decrease in the hydration of the negatively charged residues on mutation may decrease the solubility of the mutant, which increases the observed aggregation propensity of the FAD mutants. Results obtained from the theory show an excellent match with the experimentally reported data.

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.

The Formation Mechanism of Segmented Ring-Shaped Aβ Oligomers and Protofibrils

A clear understanding of amyloid formation with diverse morphologies is critical to overcoming the fatal disease amyloidosis. Studies have revealed that monomer concentration is a crucial factor for determining amyloid morphologies, such as protofibrils, annular, or spherical oligomers. However, gaining a complete understanding of the mechanism of formation of the various amyloid morphologies has been limited by the lack of experimental devices and insufficient knowledge. In this study, we demonstrate that the monomer concentration is an essential factor in determining the morphology of beta-amyloid (Aβ) oligomers or protofibrils. By computational and experimental approaches, we investigated the strategies for structural stabilization of amyloid protein, the morphological changes, and amyloid aggregation. In particular, we found unprecedented conformations, e.g., single bent oligomers and segmented ring-shaped protofibrils, the formation of which was explained by the computational analysis. Our findings provide insight into the structural features of amyloid molecules formed at low concentrations of monomer, which will help determine the clinical targets (in therapy) to effectively inhibit amyloid formation in the early stages of the amyloid growth phase.

Structural Modulation of Human Amylin Protofilaments by Naturally Occurring Mutations

Human islet amyloid polypeptide (hIAPP), also known as amylin, is a 37-amino-acid peptide, co-secreted with insulin, and widely found in fibril form in type-2 diabetes patients. By using all-atom molecular dynamics simulations, we study hIAPP fibril segments (i.e., fibrillar oligomers) formed with sequences of naturally occurring variants from cat, rat, and pig, presenting different aggregation propensities. We characterize the effect of mutations on the structural dynamics of solution-formed hIAPP fibril models built from solid-state NMR data. Results from this study are in agreement with experimental observations regarding their respective relative aggregation propensities. We analyze in detail the specific structural characteristics and infer mechanisms that modulate the conformational stability of amylin fibrils. Results provide a platform for further studies and the design of new drugs that could interfere with amylin aggregation and its cytotoxicity. One particular mutation, N31K, has fibril-destabilizing properties, and could potentially improve the solubility of therapeutic amylin analogs.

Inhibition of GNNQQNY prion peptide aggregation by trehalose: a mechanistic view

Deposition of amyloid fibrils is the seminal event in the pathogenesis of numerous neurodegenerative diseases. The formation of this amyloid assembly is the manifestation of a cascade of structural transitions including toxic oligomer formation in the early stages of aggregation. Thus a viable therapeutic strategy involves the use of small molecular ligands to interfere with this assembly. In this perspective, we have explored the kinetics of aggregate formation of the fibril forming GNNQQNY peptide fragment from the yeast prion protein SUP35 using multiple all atom MD simulations with explicit solvent and provided mechanistic insights into the way trehalose, an experimentally known aggregation inhibitor, modulates the aggregation pathway. The results suggest that the assimilation process is impeded by different barriers at smaller and larger oligomeric sizes: the initial one being easily surpassed at higher temperatures and peptide concentrations. The kinetic profile demonstrates that trehalose delays the aggregation process by increasing both these activation barriers, specifically the latter one. It increases the sampling of small-sized aggregates that lack the beta sheet conformation. Analysis reveals that the barrier in the growth of larger stable oligomers causes the formation of multiple stable small oligomers which then fuse together bimolecularly. The PCA of 26 properties was carried out to deconvolute the events within the temporary lag phases, which suggested dynamism in lags involving an increase in interchain contacts and burial of SASA. The predominant growth route is monomer addition, which changes to condensation on account of a large number of depolymerisation events in the presence of trehalose. The favourable interaction of trehalose specifically with the sidechain of the peptide promotes crowding of trehalose molecules in its vicinity - the combination of both these factors imparts the observed behaviour. Furthermore, increasing trehalose concentration leads to faster expulsion of water molecules than interpeptide interactions. These expelled water molecules have larger translational movement, suggesting an entropy factor to favor the assembly process. Different conformations observed under this condition suggest the role of water molecules in guiding the morphology of the aggregates as well. A similar scenario exists on increasing peptide concentration.

The impact of water on the ambivalent behavior and paradoxical phenomenon of the amyloid-β fibril protein

The crucial role of water in amyloid-β(Aβ) fibril proteins is evaluated in several ways including the water's thermodynamic and kinetic solvation effects. As regards the water's character, its hindered-rotation barriers are also considered. The following protein molecules considered here are: the Aβ40 (PDB ID: 2LMN), Aβ42 (PDB ID: 5KK3 and 2NAO) and the double-layered Aβ17-42 fibril. We discuss: (i) extracellular Aβ40 and Aβ42 fibril monomers exhibit an ambivalent propensity to transform into a helical form toward the N-term region and a β-strand-like form near the C-terminal; (ii) interfacial water molecules play a crucial role in protein-protein interactions, as molecular dynamics simulations have shown a significant impact on the protein-protein binding; (iii) it is shown that the spontaneous dimerization process of the Aβ42 fibril protein in water occurs via a two-step nucleation-accommodation mechanism; (iv) MD simulations of the double-layered Aβ17-42 fibril model show that the C↔C interface appears more energetically favorable than the N↔N interface due to large hydrophobic contacts; (v) the water's role in the HET-s prion and in the Aβ fibrillar aggregates; (vi) it was found that the monomer-oligomer equilibrium spontaneously dissociates into stable monomeric species when they are incubated up to 3 μm for a longer time (>1 week) in a physiological buffer.

Statistical Thermodynamics for Actin-Myosin Binding: The Crucial Importance of Hydration Effects

Actomyosin is an important molecular motor, and the binding of actin and myosin is an essential research target in biophysics. Nevertheless, the physical factors driving or opposing the binding are still unclear. Here, we investigate the role of water in actin-myosin binding using the most reliable statistical-mechanical method currently available for assessing biomolecules immersed in water. This method is characterized as follows: water is treated not as a dielectric continuum but as an ensemble of molecules; the polyatomic structures of proteins are taken into consideration; and the binding free energy is decomposed into physically insightful entropic and energetic components by accounting for the hydration effect to its full extent. We find that the actin-myosin binding brings large gains of electrostatic and Lennard-Jones attractive interactions. However, these gains are accompanied by even larger losses of actin-water and myosin-water electrostatic and LJ attractive interactions. Although roughly half of the energy increase due to the losses is cancelled out by the energy decrease arising from structural reorganization of the water released upon binding, the remaining energy increase is still larger than the energy decrease brought by the gains mentioned above. Hence, the net change in system energy is positive, which opposes binding. Importantly, the binding is driven by a large gain of configurational entropy of water, which surpasses the positive change in system energy and the conformational entropy loss occurring for actin and myosin. The principal physical origin of the large water-entropy gain is as follows: the actin-myosin interface is closely packed with the achievement of high shape complementarity on the atomic level, leading to a large increase in the total volume available to the translational displacement of water molecules in the system and a resultant reduction of water crowding (i.e., entropic correlations among water molecules).

Dimerization Mechanism of Alzheimer Aβ40 Peptides: The High Content of Intrapeptide-Stabilized Conformations in A2V and A2T Heterozygous Dimers Retards Amyloid Fibril Formation

Amyloid beta (Aβ) oligomerization is associated with the origin and progression of Alzheimer's disease (AD). While the A2V mutation enhances aggregation kinetics and toxicity, mixtures of wild-type (WT) and A2V, and also WT and A2T, peptides retard fibril formation and protect against AD. In this study, we simulate the equilibrium ensemble of WT:A2T Aβ₄₀ dimer by means of extensive atomistic replica exchange molecular dynamics and compare our results with previous equivalent simulations of A2V:A2V, WT:WT, and WT:A2V Aβ₄₀ dimers for a total time scale of nearly 0.1 ms. Qualitative comparison of the resulting thermodynamic properties, such as the relative binding free energies, with the reported experimental kinetic and thermodynamic data affords us important insight into the conversion from slow-pathway to fast-pathway dimer conformations. The crucial reaction coordinate or driving force of such transformation turns out to be related to hydrophobic interpeptide interactions. Analysis of the equilibrium ensembles shows that the fast-pathway conformations contain interpeptide out-of-register antiparallel β-sheet structures at short interpeptide distances. In contrast, the slow-pathway conformations are formed by the association of peptides at large interpeptide distances and high intrapeptide compactness, such as conformations containing intramolecular three-stranded β-sheets which sharply distinguish fast (A2V:A2V and WT:WT) and slow (WT:A2T and WT:A2V) amyloid-forming sequences. Also, this analysis leads us to predict that a molecule stabilizing the intramolecular three-stranded β-sheet or inhibiting the formation of an interpeptide β-sheet spanning residues 17-20 and 31-37 would further reduce fibril formation and probably the cytotoxicity of Aβ species.

Peptide self-assembly: thermodynamics and kinetics

Self-assembling systems play a significant role in physiological functions and have therefore attracted tremendous attention due to their great potential for applications in energy, biomedicine and nanotechnology. Peptides, consisting of amino acids, are among the most popular building blocks and programmable molecular motifs. Nanostructures and materials assembled using peptides exhibit important potential for green-life new technology and biomedical applications mostly because of their bio-friendliness and reversibility. The formation of these ordered nanostructures pertains to the synergistic effect of various intermolecular non-covalent interactions, including hydrogen-bonding, π-π stacking, electrostatic, hydrophobic, and van der Waals interactions. Therefore, the self-assembly process is mainly driven by thermodynamics; however, kinetics is also a critical factor in structural modulation and function integration. In this review, we focus on the influence of thermodynamic and kinetic factors on structural assembly and regulation based on different types of peptide building blocks, including aromatic dipeptides, amphiphilic peptides, polypeptides, and amyloid-relevant peptides.

New Computational Approach for External Entropy in Protein-Protein Binding

Molecular recognition through the noncovalent association of biomolecules is of central importance in biology and pharmacology. Developing reliable computational methods for estimating binding thermodynamic parameters is therefore of great practical value. However, considerable uncertainty remains regarding the external entropy that is associated with the reduction in the external (positional and orientational) degrees of freedom upon complex formation. Here, we present a novel statistical mechanical method for computing the external entropy by extending the energetic approach we have developed for unimolecular processes to association processes. We find that, in contrary to what is postulated in most of the previous methods, intrinsic couplings between the internal and external degrees of freedom of bound complex cannot in general be neglected in the determination of the external entropy. Nevertheless, there exists the best choice of the external coordinates with which those couplings are minimized. With such a choice of the external coordinates, the lowest upper bound of the external entropy is obtained from a tractable expression, which serves as an estimate of the external entropy. Our method can be implemented in a straightforward manner with molecular dynamics simulations, and its applicability is demonstrated through the application to the barnase-barstar complex.

Statistical thermodynamics of aromatic–aromatic interactions in aqueous solution

To elucidate the interactions between aromatic rings, which are believed to play essential roles in a variety of biological processes, we analyze the water-mediated interactions between toluene molecules along face-to-face stacked (FF) and point-to-face T-shaped (TS) paths using a statistical-mechanical theory of liquids combined with a molecular model for water.

Mechanism of One-to-Many Molecular Recognition Accompanying Target-Dependent Structure Formation: For the Tumor Suppressor p53 Protein as an Example

The new type of molecular recognition, in which an intrinsically disordered region (IDR) of a protein binds to many different target proteins with target-dependent structure formation, is indispensable to the expression of life phenomena and also implicated in a number of diseases. According to the prevailing view, the physicochemical factors responsible for the binding are also target dependent. Here we consider an IDR of the tumor suppressor p53 protein, p53CTD, as an important example related to carcinogenesis and analyze its binding to four targets accompanying the formation of target-dependent structures (i.e., helix, sheet, and two different coils) using our statistical-mechanical method combined with molecular models for water. We find that all of the seemingly different binding processes are driven by a large gain of the translational, configurational entropy of water in the system. The gain originates from sufficiently high shape complementarity on the atomic level within the p53CTD-target interface. It is also required that the electrostatic complementarity be ensured as much as possible to compensate for the dehydration. Such complementarities are achieved in harmony with the portion of the target to which p53CTD binds, leading to a large diversity of structures of p53CTD formed upon binding: If they are not achievable, the binding does not occur. This finding is made possible only by calculating the changes in thermodynamic quantities upon binding and decomposing them into physically insightful components.

Distinct role of hydration water in protein misfolding and aggregation revealed by fluctuating thermodynamics analysis

Protein aggregation in aqueous cellular environments is linked to diverse human diseases. Protein aggregation proceeds through a multistep process initiated by conformational transitions, called protein misfolding, of monomer species toward aggregation-prone structures. Various forms of aggregate species are generated through the association of misfolded monomers including soluble oligomers and amyloid fibrils. Elucidating the molecular mechanisms and driving forces involved in the misfolding and subsequent association has been a central issue for understanding and preventing protein aggregation diseases such as Alzheimer's, Parkinson's, and type II diabetes. In this Account, we provide a thermodynamic perspective of the misfolding and aggregation of the amyloid-beta (Aβ) protein implicated in Alzheimer's disease through the application of fluctuating thermodynamics. This approach "dissects" the conventional thermodynamic characterization of the end states into the one of the fluctuating processes connecting them, and enables one to analyze variations in the thermodynamic functions that occur during the course of protein conformational changes. The central quantity in this approach is the solvent-averaged effective energy, f = Eu + Gsolv, comprising the protein potential energy (Eu) and the solvation free energy (Gsolv), whose time variation reflects the protein dynamics on the free energy landscape. Protein configurational entropy is quantified by the magnitude of fluctuations in f. We find that misfolding of the Aβ monomer when released from a membrane environment to an aqueous phase is driven by favorable changes in protein potential energy and configurational entropy, but it is also accompanied by an unfavorable increase in solvation free energy. The subsequent dimerization of the misfolded Aβ monomers occurs in two steps. The first step, where two widely separated monomers come into contact distance, is driven by water-mediated attraction, that is, by a decrease in solvation free energy, harnessing the monomer solvation free energy earned during the misfolding. The second step, where a compact dimer structure is formed, is driven by direct protein-protein interactions, but again it is accompanied by an increase in solvation free energy. The increased solvation free energy of the dimer will function as the driving force to recruit another Aβ protein in the approach stage of subsequent oligomerizations. The fluctuating thermodynamics analysis of the misfolding and dimerization of the Aβ protein indicates that the interaction of the protein with surrounding water plays a critical role in protein aggregation. Such a water-centric perspective is further corroborated by demonstrating that, for a large number of Aβ mutants and mutants of other protein systems, the change in the experimental aggregation propensity upon mutation has a significant correlation with the protein solvation free energy change. We also find striking discrimination between the positively and negatively charged residues on the protein surface by surrounding water molecules, which is shown to play a crucial role in determining the protein aggregation propensity. We argue that the protein total charge dictates such striking behavior of the surrounding water molecules. Our results provide new insights for understanding and predicting the protein aggregation propensity, thereby offering novel design principles for producing aggregation-resistant proteins for biotherapeutics.

Crucial role of nonspecific interactions in amyloid nucleation

Protein oligomers have been implicated as toxic agents in a wide range of amyloid-related diseases. However, it has remained unsolved whether the oligomers are a necessary step in the formation of amyloid fibrils or just a dangerous byproduct. Analogously, it has not been resolved if the amyloid nucleation process is a classical one-step nucleation process or a two-step process involving prenucleation clusters. We use coarse-grained computer simulations to study the effect of nonspecific attractions between peptides on the primary nucleation process underlying amyloid fibrillization. We find that, for peptides that do not attract, the classical one-step nucleation mechanism is possible but only at nonphysiologically high peptide concentrations. At low peptide concentrations, which mimic the physiologically relevant regime, attractive interpeptide interactions are essential for fibril formation. Nucleation then inevitably takes place through a two-step mechanism involving prefibrillar oligomers. We show that oligomers not only help peptides meet each other but also, create an environment that facilitates the conversion of monomers into the β-sheet-rich form characteristic of fibrils. Nucleation typically does not proceed through the most prevalent oligomers but through an oligomer size that is only observed in rare fluctuations, which is why such aggregates might be hard to capture experimentally. Finally, we find that the nucleation of amyloid fibrils cannot be described by classical nucleation theory: in the two-step mechanism, the critical nucleus size increases with increases in both concentration and interpeptide interactions, which is in direct contrast with predictions from classical nucleation theory.

ATP-induced conformational changes of nucleotide-binding domains in an ABC transporter. Importance of the water-mediated entropic force

ATP binding cassette (ABC) proteins belong to a superfamily of active transporters. Recent experimental and computational studies have shown that binding of ATP to the nucleotide binding domains (NBDs) of ABC proteins drives the dimerization of NBDs, which, in turn, causes large conformational changes within the transmembrane domains (TMDs). To elucidate the active substrate transport mechanism of ABC proteins, it is first necessary to understand how the NBD dimerization is driven by ATP binding. In this study, we selected MalKs (NBDs of a maltose transporter) as a representative NBD and calculated the free-energy change upon dimerization using molecular mechanics calculations combined with a statistical thermodynamic theory of liquids, as well as a method to calculate the translational, rotational, and vibrational entropy change. This combined method is applied to a large number of snapshot structures obtained from molecular dynamics simulations containing explicit water molecules. The results suggest that the NBD dimerization proceeds with a large gain of water entropy when ATP molecules bind to the NBDs. The energetic gain arising from direct NBD-NBD interactions is canceled by the dehydration penalty and the configurational-entropy loss. ATP hydrolysis induces a loss of the shape complementarity between the NBDs, which leads to the dissociation of the dimer, due to a decrease in the water-entropy gain and an increase in the configurational-entropy loss. This interpretation of the NBD dimerization mechanism in concert with ATP, especially focused on the water-mediated entropy force, is potentially applicable to a wide variety of the ABC transporters.

Binding of an RNA aptamer and a partial peptide of a prion protein: crucial importance of water entropy in molecular recognition

It is a central issue to elucidate the new type of molecular recognition accompanied by a global structural change of a molecule upon binding to its targets. Here we investigate the driving force for the binding of R12 (a ribonucleic acid aptamer) and P16 (a partial peptide of a prion protein) during which P16 exhibits the global structural change. We calculate changes in thermodynamic quantities upon the R12–P16 binding using a statistical-mechanical approach combined with molecular models for water which is currently best suited to studies on hydration of biomolecules. The binding is driven by a water-entropy gain originating primarily from an increase in the total volume available to the translational displacement of water molecules in the system. The energy decrease due to the gain of R12–P16 attractive (van der Waals and electrostatic) interactions is almost canceled out by the energy increase related to the loss of R12–water and P16–water attractive interactions. We can explain the general experimental result that stacking of flat moieties, hydrogen bonding and molecular-shape and electrostatic complementarities are frequently observed in the complexes. It is argued that the water-entropy gain is largely influenced by the geometric characteristics (overall shapes, sizes and detailed polyatomic structures) of the biomolecules.