Atomic decomposition of the protein solvation free energy and its application to amyloid-beta protein in water

Synergistic Modeling of Liquid Properties: Integrating Neural Network-Derived Molecular Features with Modified Kernel Models

A significant challenge in applying machine learning to computational chemistry, particularly considering the growing complexity of contemporary machine learning models, is the scarcity of available experimental data. To address this issue, we introduce an approach that derives molecular features from an intricate neural network-based model and applies them to a simpler conventional machine learning model that is robust to overfitting. This method can be applied to predict various properties of a liquid system, including viscosity or surface tension, based on molecular features drawn from the ab initio calculated free energy of solvation. Furthermore, we propose a modified kernel model that includes Arrhenius temperature dependence to incorporate theoretical principles and diminish extreme nonlinearity in the model. The modified kernel model demonstrated significant improvements in certain scenarios and possible extensions to various theoretical concepts of molecular systems.

Dynamic protein hydration water mediates the aggregation kinetics of amyloid β peptides at interfaces

Protein hydration water is essential for protein misfolding and amyloid formation, but how it directs the course of amyloid formation has yet to be elucidated. Here, we experimentally demonstrated that femtosecond sum frequency generation vibrational spectroscopy (SFG-VS) and the femtosecond IR pump-SFG probe technique can serve as powerful tools for addressing this issue. Using amyloid β(1-42) peptide as a model, we determined the transient misfolding intermediates by probing the amide band spectral features and the local hydration water changes by measuring the ultrafast vibrational dynamics of the amide I band. For the first time, we established a correlation between the dynamic change in protein hydration water and aggregation propensity. The aggregation propensity depends on the dynamic change in the hydration water, rather than the static hydration water content of the initial protein state. Water expulsion enhances the aggregation propensity and promotes amyloid formation, while protein hydration attenuates the aggregation propensity and inhibits amyloid formation. The suppression of water expulsion and protein hydration can prevent protein aggregation and stabilize proteins. These findings contribute to a better understanding of the underlying effect of hydration water on amyloid formation and protein structural stability and provide a strategy for maintaining long-term stabilization of biomolecules.

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.

Effect of Main and Side Chains on the Folding Mechanism of the Trp-Cage Miniprotein

Proteins that do not fold into their functional native state have been linked to diseases. In this study, the influence of the main and side chains of individual amino acids on the folding of the tryptophan cage (Trp-cage), a designed 20-residue miniprotein, was analyzed. For this purpose, we calculated the solvation free energy (SFE) contributions of individual atoms by using the 3D-reference interaction site model with the atomic decomposition method. The mechanism by which the Trp-cage is stabilized during the folding process was examined by calculating the total energy, which is the sum of the conformational energy and SFE. The folding process of the Trp-cage resulted in a stable native state, with a total energy that was 62.4 kcal/mol lower than that of the unfolded state. The solvation entropy, which is considered to be responsible for the hydrophobic effect, contributed 31.3 kcal/mol to structural stabilization. In other words, the contribution of the solvation entropy accounted for approximately half of the total contribution to Trp-cage folding. The hydrophobic core centered on Trp6 contributed 15.6 kcal/mol to the total energy, whereas the solvation entropy contribution was 6.3 kcal/mol. The salt bridge formed by the hydrophilic side chains of Asp9 and Arg16 contributed 10.9 and 5.0 kcal/mol, respectively. This indicates that not only the hydrophobic core but also the salt bridge of the hydrophilic side chains gain solvation entropy and contribute to stabilizing the native structure of the Trp-cage.

Atomic-level thermodynamics analysis of the binding free energy of SARS-CoV-2 neutralizing antibodies

Understanding how protein-protein binding affinity is determined from molecular interactions at the interface is essential in developing protein therapeutics such as antibodies, but this has not yet been fully achieved. Among the major difficulties are the facts that it is generally difficult to decompose thermodynamic quantities into contributions from individual molecular interactions and that the solvent effect-dehydration penalty-must also be taken into consideration for every contact formation at the binding interface. Here, we present an atomic-level thermodynamics analysis that overcomes these difficulties and illustrate its utility through application to SARS-CoV-2 neutralizing antibodies. Our analysis is based on the direct interaction energy computed from simulated antibody-protein complex structures and on the decomposition of solvation free energy change upon complex formation. We find that the formation of a single contact such as a hydrogen bond at the interface barely contributes to binding free energy due to the dehydration penalty. On the other hand, the simultaneous formation of multiple contacts between two interface residues favorably contributes to binding affinity. This is because the dehydration penalty is significantly alleviated: the total penalty for multiple contacts is smaller than a sum of what would be expected for individual dehydrations of those contacts. Our results thus provide a new perspective for designing protein therapeutics of improved binding affinity.

Mutation effects on FAS1 domain 4 based on structure and solubility

Mutations in the fasciclin 1 domain 4 (FAS1-4) of transforming growth factor β-induced protein (TGFBIp) are associated with insoluble extracellular deposits and corneal dystrophies (CDs). The decrease in solubility upon mutation has been implicated in CD; however, the exact molecular mechanisms are not well understood. Here, we performed molecular dynamics simulations followed by solvation thermodynamic analyses of the FAS1-4 domain and its three mutants-R555W, R555Q, and A546T-linked to granular corneal dystrophy type 1, Thiel-Behnke corneal dystrophy and lattice corneal dystrophy, respectively. We found that both R555W and R555Q mutants have less affinity toward solvent water relative to the wild-type protein. In the R555W mutant, a remarkable increase in solvation free energy was observed because of the structural changes near the mutation site. The mutation site W555 is buried in other hydrophobic residues, and R557 simultaneously forms salt bridges with E554 and D561. In the R555Q mutant, the increase in solvation free energy is caused by structural rearrangements far from the mutation site. R558 separately forms salt bridges with D575, E576, and E598. Thus, we thus identified the relationship between the decrease in solubility and conformational changes caused by mutations, which may be useful in designing potential therapeutics and in blocking FAS1 aggregation related to CD.

MLSolvA: solvation free energy prediction from pairwise atomistic interactions by machine learning

Recent advances in machine learning technologies and their applications have led to the development of diverse structure-property relationship models for crucial chemical properties. The solvation free energy is one of them. Here, we introduce a novel ML-based solvation model, which calculates the solvation energy from pairwise atomistic interactions. The novelty of the proposed model consists of a simple architecture: two encoding functions extract atomic feature vectors from the given chemical structure, while the inner product between the two atomistic feature vectors calculates their interactions. The results of 6239 experimental measurements achieve outstanding performance and transferability for enlarging training data owing to its solvent-non-specific nature. An analysis of the interaction map shows that our model has significant potential for producing group contributions on the solvation energy, which indicates that the model provides not only predictions of target properties but also more detailed physicochemical insights.

Site-Specific Backbone and Side-Chain Contributions to Thermodynamic Stabilizing Forces of the WW Domain

The native structure of a protein is stabilized by a number of interactions such as main-chain hydrogen bonds and side-chain hydrophobic contacts. However, it has been challenging to determine how these interactions contribute to protein stability at single amino acid resolution. Here, we quantified site-specific thermodynamic stability at the molecular level to extend our understanding of the stabilizing forces in protein folding. We derived the free energy components of individual amino acid residues separately for the folding of the human Pin WW domain based on simulated structures. A further decomposition of the thermodynamic properties into contributions from backbone and side-chain groups enabled us to identify the critical residues in the secondary structure and hydrophobic core formation, without introducing physical modifications to the system as in site-directed mutagenesis methods. By relating the structural and thermodynamic changes upon folding for each residue, we find that the simultaneous formation of the backbone hydrogen bonds and side-chain contacts cooperatively stabilizes the folded structure. The identification of stabilizing interactions in a folding protein at atomic resolution will provide molecular insights into understanding the origin of the protein structure and into engineering a more stable protein.

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.

Local environment effects on charged mutations for developing aggregation-resistant monoclonal antibodies

Protein aggregation is a major concern in biotherapeutic applications of monoclonal antibodies. Introducing charged mutations is among the promising strategies to improve aggregation resistance. However, the impact of such mutations on solubilizing activity depends largely on the inserting location, whose mechanism is still not well understood. Here, we address this issue from a solvation viewpoint, and this is done by analyzing how the change in solvation free energy upon charged mutation is composed of individual contributions from constituent residues. To this end, we perform molecular dynamics simulations for a number of antibody mutants and carry out the residue-wise decomposition of the solvation free energy. We find that, in addition to the previously identified “global” principle emphasizing the key role played by the protein total net charge, a local net charge within \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim$$\end{document}∼15 Å from the mutation site exerts significant effects. For example, when the net charge of an antibody is positive, the global principle states that introducing a positively charged mutation will lead to more favorable solvation. Our finding further adds that an even more optimal mutation can be done at the site around which more positively charged residues and fewer negatively charged residues are present. Such a “local” design principle accounts for the location dependence of charged mutations, and will be useful in producing aggregation-resistant antibodies.

Effect of linker on the binding free energy of stapled p53/HDM2 complex

Inactivation of the tumor suppressor p53 resulting from the binding with a negative regulator HDM2 is among the predominant defects in human cancers. p53-mimicking peptides whose conformational and proteolytic stability is enhanced by an all-hydrocarbon staple are being recognized as promising anticancer agents for disrupting the p53–HDM2 binding and reactivating p53. Herein, we conduct a computational modeling and thermodynamic characterization of stapled p53/HDM2 complex via molecular docking, simulations, and binding free energy analysis. The binding thermodynamics analysis is done based on the end-point calculation of the effective binding energy—a sum of the direct peptide–protein interaction energy and the dehydration penalty—and on its decomposition into contributions from specific groups constituting the complex. This allows us to investigate how individual amino acids in the stapled p53 and HDM2 contribute to the binding affinity. We find that not only the epitope residues (F19, W23 and L26), but also the hydrocarbon linker of the stapled p53 impart significant contributions. Our computational approach will be useful in designing new stapled peptides in which the staple location is also optimized to improve the binding affinity.

Gibbs Free Energy Calculation of Mutation in PncA and RpsA Associated With Pyrazinamide Resistance

A central approach for better understanding the forces involved in maintaining protein structures is to investigate the protein folding and thermodynamic properties. The effect of the folding process is often disturbed in mutated states. To explore the dynamic properties behind mutations, molecular dynamic (MD) simulations have been widely performed, especially in unveiling the mechanism of drug failure behind mutation. When comparing wild type (WT) and mutants (MTs), the structural changes along with solvation free energy (SFE), and Gibbs free energy (GFE) are calculated after the MD simulation, to measure the effect of mutations on protein structure. Pyrazinamide (PZA) is one of the first-line drugs, effective against latent Mycobacterium tuberculosis isolates, affecting the global TB control program 2030. Resistance to this drug emerges due to mutations in pncA and rpsA genes, encoding pyrazinamidase (PZase) and ribosomal protein S1 (RpsA) respectively. The question of how the GFE may be a measure of PZase and RpsA stabilities, has been addressed in the current review. The GFE and SFE of MTs have been compared with WT, which were already found to be PZA-resistant. WT structures attained a more stable state in comparison with MTs. The physiological effect of a mutation in PZase and RpsA may be due to the difference in energies. This difference between WT and MTs, depicted through GFE plots, might be useful in predicting the stability and PZA-resistance behind mutation. This study provides useful information for better management of drug resistance, to control the global TB problem.

Mutation-induced change in chignolin stability from π-turn to α-turn

A mutation from threonine to proline at the eighth residue in chignolin changes π-turn to α-turn.

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.

Folding Free Energy Landscape of Ordered and Intrinsically Disordered Proteins

Folding funnel is the essential concept of the free energy landscape for ordered proteins. How does this concept apply to intrinsically disordered proteins (IDPs)? Here, we address this fundamental question through the explicit characterization of the free energy landscapes of the representative α-helical (HP-35) and β-sheet (WW domain) proteins and of an IDP (pKID) that folds upon binding to its partner (KIX). We demonstrate that HP-35 and WW domain indeed exhibit the steep folding funnel: the landscape slope for these proteins is ca. −50 kcal/mol, meaning that the free energy decreases by ~5 kcal/mol upon the formation of 10% native contacts. On the other hand, the landscape of pKID is funneled but considerably shallower (slope of −24 kcal/mol), which explains why pKID is disordered in free environments. Upon binding to KIX, the landscape of pKID now becomes significantly steep (slope of −54 kcal/mol), which enables otherwise disordered pKID to fold. We also show that it is the pKID–KIX intermolecular interactions originating from hydrophobic residues that mainly confer the steep folding funnel. The present work not only provides the quantitative characterization of the protein folding free energy landscape, but also establishes the usefulness of the folding funnel concept to IDPs.

Thermodynamic Stability of Human γD-Crystallin Mutants Using Alchemical Free-Energy Calculations

γD-Crystallin (HγDC) is a key structural protein in the human lens, whose aggregation has been associated with the development of cataracts. Single-point mutations and post-translational modifications destabilize HγDC interactions, forming partially folded intermediates, where hydrophobic residues are exposed and thus triggering its aggregation. In this work, we used alchemical free-energy calculations to predict changes in thermodynamic stability (ΔΔG) of 10 alanine-scanning variants and 12 HγDC mutations associated with the development of congenital cataract. Our results show that W42R is the most destabilizing mutation in HγDC. This has been corroborated through experimental determination of ΔΔG employing differential scanning calorimetry. Calculations of hydration free energies from the HγDC WT and the W42R mutant suggested that the mutant has a higher aggregation propensity. Our combined theoretical and experimental results contribute to understand HγDC destabilization and aggregation mechanisms in age-onset cataracts.

Analysis of Structural Stability of Chignolin

We discuss the stability of an entire protein and the influence of main chains and side chains of individual amino acids to investigate the protein-folding mechanism. For this purpose, we calculated the solvation free-energy contribution of individual atoms using the three-dimensional reference interaction site model with the atomic decomposition method. We generated structures of chignolin miniprotein by a molecular dynamics simulation and classified them into six types: native 1, native 2, misfolded 1, misfolded 2, intermediate, and unfolded states. The total energies of the native (-171.1 kcal/mol) and misfolded (-171.2 kcal/mol) states were almost the same and lower than those of the intermediate (-158.5 kcal/mol) and unfolded (-148.1 kcal/mol) states; however, their components were different. In the native state, the side-chain interaction between Thr6 and Thr8 is important for the formation of π-turn. On the other hand, the hydrogen bonds between the atoms of the main chains in the misfolded state become stronger than those in the intermediate state.

Dynamics of Hydration Water Plays a Key Role in Determining the Binding Thermodynamics of Protein Complexes

Interfacial waters are considered to play a crucial role in protein–protein interactions, but in what sense and why are they important? Here, using molecular dynamics simulations and statistical thermodynamic analyses, we demonstrate distinctive dynamic characteristics of the interfacial water and investigate their implications for the binding thermodynamics. We identify the presence of extraordinarily slow (~1,000 times slower than in bulk water) hydrogen-bond rearrangements in interfacial water. We rationalize the slow rearrangements by introducing the “trapping” free energies, characterizing how strongly individual hydration waters are captured by the biomolecular surface, whose magnitude is then traced back to the number of water–protein hydrogen bonds and the strong electrostatic field produced at the binding interface. We also discuss the impact of the slow interfacial waters on the binding thermodynamics. We find that, as expected from their slow dynamics, the conventional approach to the water-mediated interaction, which assumes rapid equilibration of the waters’ degrees of freedom, is inadequate. We show instead that an explicit treatment of the extremely slow interfacial waters is critical. Our results shed new light on the role of water in protein–protein interactions, highlighting the need to consider its dynamics to improve our understanding of biomolecular bindings.

Structural properties of amyloid β(1-40) dimer explored by replica exchange molecular dynamics simulations

Replica exchange molecular dynamics simulations (300 ns) were used to study the dimerization of amyloid β(1-40) (Aβ(1-40)) polypeptide. Configurational entropy calculations revealed that at physiological temperature (310 K, 37°C) dynamic dimers are formed by randomly docked monomers. Free energy of binding of the two chains to each other was -93.56 ± 6.341 kJ mol^-1 . Prevalence of random coil conformations was found for both chains with the exceptions of increased β-sheet content from residues 16-21 and 29-32 of chain A and residues 15-21 and 30-33 of chain B with β-turn/β-bend conformations in both chains from residues 1-16, 21-29 of chain A, 1-16, and 21-29 of chain B. There is a mixed β-turn/β-sheet region from residues 33-38 of both chains. Analysis of intra- and interchain residue distances shows that, although the individual chains are highly flexible, the dimer system stays in a loosely packed antiparallel β-sheet configuration with contacts between residues 17-21 of chain A with residues 17-21 and 31-36 of chain B as well as residues 31-36 of chain A with residues 17-21 and 31-36 of chain B. Based on dihedral principal component analysis, the antiparallel β-sheet-loop-β-sheet conformational motif is favored for many low energy sampled conformations. Our results show that Aβ(1-40) can form dynamic dimers in aqueous solution that have significant conformational flexibility and are stabilized by collapse of the central and C-terminal hydrophobic cores with the expected β-sheet-loop-β-sheet conformational motif. Proteins 2017; 85:1024-1045. © 2017 Wiley Periodicals, Inc.

Chemical Fluorescent Probe for Detection of Aβ Oligomers

Aggregation of amyloid β-peptide (Aβ) is implicated in the pathology of Alzheimer's disease (AD), with the soluble, Aβ oligomeric species thought to be the critical pathological species. Identification and characterization of intermediate species formed during the aggregation process is crucial to the understanding of the mechanisms by which oligomeric species mediate neuronal toxicity and following disease progression. Probing these species proved to be extremely challenging, as evident by the lack of reliable sensors, due to their heterogeneous and transient nature. We describe here an oligomer-specific fluorescent chemical probe, BoDipy-Oligomer (BD-Oligo), developed through the use of the diversity-oriented fluorescent library approach (DOFLA) and high-content, imaging-based screening. This probe enables dynamic oligomer monitoring during fibrillogenesis in vitro and shows in vivo Aβ oligomers staining possibility in the AD mice model.

Structural and Thermodynamic Characteristics of Amyloidogenic Intermediates of β-2-Microglobulin

β-2-microglobulin (β2m) self-aggregates to form amyloid fibril in renal patients taking long-term dialysis treatment. Despite the extensive structural and mutation studies carried out so far, the molecular details on the factors that dictate amyloidogenic potential of β2m remain elusive. Here we report molecular dynamics simulations followed by the solvation thermodynamic analyses on the wild-type β2m and D76N, D59P, and W60C mutants at the native (N) and so-called aggregation-prone intermediate (I_T) states, which are distinguished by the native cis- and non-native trans-Pro32 backbone conformations. Three major structural and thermodynamic characteristics of the I_T-state relative to the N-state in β2m protein are detected that contribute to the increased amyloidogenic potential: (i) the disruption of the edge D-strand, (ii) the increased solvent-exposed hydrophobic interface, and (iii) the increased solvation free energy (less affinity toward solvent water). Mutation effects on these three factors are shown to exhibit a good correlation with the experimentally observed distinct amyloidogenic propensity of the D76N (+), D59P (+), and W60C (−) mutants (+/− for enhanced/decreased). Our analyses thus identify the structural and thermodynamic characteristics of the amyloidogenic intermediates, which will serve to uncover molecular mechanisms and driving forces in β2m amyloid fibril formation.

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.

Interaction with the surrounding water plays a key role in determining the aggregation propensity of proteins

Understanding the molecular determinants of the relative propensities of proteins to aggregate in a cellular environment is a central issue for treating protein-aggregation diseases and developing peptide-based therapeutics. Despite the expectation that protein aggregation can largely be attributed to direct protein-protein interactions, a crucial role the surrounding water in determining the aggregation propensity of proteins both in vitro and in vivo was identified. The overall protein hydrophobicity, defined solely by the hydration free energy of a protein in its monomeric state sampling its equilibrium structures, was shown to be the predominant determinant of protein aggregation propensity in aqueous solution. Striking discrimination of positively and negatively charged residues by the surrounding water was also found. This effect depends on the protein net charge and plays a crucial role in regulating the solubility of the protein. These results pave the way for the design of aggregation-resistant proteins as biotherapeutics.

Impact of chemical heterogeneity on protein self-assembly in water

Hydrophobicity is thought to underlie self-assembly in biological systems. However, the protein surface comprises hydrophobic and hydrophilic patches, and understanding the impact of such a chemical heterogeneity on protein self-assembly in water is of fundamental interest. Here, we report structural and thermodynamic investigations on the dimer formation of full-length amyloid-β proteins in water associated with Alzheimer's disease. Spontaneous dimerization process--from the individual diffusive regime at large separations, through the approach stage in which two proteins come close to each other, to the structural adjustment stage toward compact dimer formation--was captured in full atomic detail via unguided, explicit-water molecular dynamics simulations. The integral-equation theory of liquids was then applied to simulated protein structures to analyze hydration thermodynamic properties and the water-mediated interaction between proteins. We demonstrate that hydrophilic residues play a key role in initiating the dimerization process. A long-range hydration force of enthalpic origin acting on the hydrophilic residues provides the major thermodynamic force that drives two proteins to approach from a large separation to a contact distance. After two proteins make atomic contacts, the nature of the water-mediated interaction switches from a long-range enthalpic attraction to a short-range entropic one. The latter acts both on the hydrophobic and hydrophilic residues. Along with the direct protein-protein interactions that lead to the formation of intermonomer hydrogen bonds and van der Waals contacts, the water-mediated attraction of entropic origin brings about structural adjustment of constituent monomer proteins toward the formation of a compact dimer structure.

Molecular interactions of Alzheimer's Aβ protofilaments with lipid membranes

Amyloid fibrils and peptide oligomers play central roles in the pathology of Alzheimer's disease, type 2 diabetes, Parkinson's disease, Huntington's disease, and prion-related disease. Here, we investigate the molecular interactions between preformed amyloid β (Aβ) molecular protofilaments and lipid bilayer membranes, in the presence of explicit water molecules, using computational models and all-atom molecular dynamics. These interactions play an important role in the stability and function of both Aβ fibrils and the adjacent cellular membrane. Taking advantage of the symmetry-related and directional properties of the protofilaments, we build models that cover several relative protofilament-membrane orientations. Our molecular dynamics simulations reveal the relative contributions of different structural elements to the dynamics and stability of Aβ protofilament segments near membranes, and the first steps in the mechanism of fibril-membrane interactions. During this process, we observe a significant alteration of the side-chain contact pattern in protofilaments, although a fraction of the characteristic β-sheet content is preserved. As a major driving force, we identify the electrostatic interactions between Aβ charged side chains, including E22, D23, and K28, and lipid headgroups. Together with hydrogen bonding with atoms from lipid headgroups, these interactions can facilitate the penetration of hydrophobic C-terminal amino acids through the lipid headgroup region, which can finally lead both to further loss of the initial fibril structure and to local membrane-thinning effects. Our results may guide new experiments that could test the extent to which the structural features of water-formed amyloid fibrils are preserved, lost, or reshaped by membrane-mediated interactions.