Dimerization of Full-length Aß -42 Peptide: A Comparison of Different Force Fields and Water Models.

Among the two isoforms of amyloid- i.e., Aß-40 and Aß-42, Aß-42 is more toxic due to its increased aggregation propensity. The oligomerization pathways of amyloid-ß may be investigated by studying its dimerization process at an atomic level. Intrinsically disordered proteins (IDPs) lack well-defined structures and are associated with numerous neurodegenerative disorders. Molecular dynamics simulations of these proteins are often limited by the choice of parameters due to inconsistencies in the empirically developed protein force fields and water models. To evaluate the accuracy of recently developed force fields for IDPs, we study the dimerization of full-length Aß-42 in aqueous solution with three different combinations of AMBER force field parameters and water models such as ff14SB/TIP3P, ff19SB/OPC, and ff19SB/TIP3P using classical MD and Umbrella Sampling method. This work may be used as a benchmark to compare the performance of different force fields for the simulations of IDPs.

Molecular Dynamics Simulation Study of the Self-Assembly of Tau-Derived PHF6 and Its Inhibition by Oleuropein Aglycone from Extra Virgin Olive Oil.

Alzheimer's disease (AD) and other taupathies are neurodegenerative disorders associated with the amyloid deposition of the Tau protein in the brain. This amyloid formation may be inhibited by small molecules, which is recognized as one of the best therapeutic strategies to stop the progression of the disease. This work focuses on the small nucleating segment, hexapeptide-paired helical filament 6 (PHF6), responsible for Tau aggregation. Using computational modeling and classical molecular dynamics simulations, we show that PHF6 monomers collapse in water to form ß-sheet rich structures, and the main olive oil polyphenol oleuropein aglycone (OleA) prevents peptide aggregation significantly. We gradually increase the ratio of the PHF6-OleA from 1:1 to 1:3 and find that for the 1:1 ratio, the peptide monomers are prone to form aggregated structures, while for the 1:2 ratio, the formation of the extended ß-sheet structure is significantly less. For a 1:3 ratio of protein/OleA, the peptide residues are sufficiently crowded by OleA molecules through hydrogen bonding, hydrophobic interactions, and π-π stacking; hence, the peptide chains prefer to exist in a monomeric random coil conformation.

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.

Intrinsic viscosity and dielectric relaxation of ring polymers in dilute solutions.

The absence of chain ends makes ring polymers distinctly different from their linear analogues. The intrinsic viscosity, complex viscosity and the dielectric relaxation of ring polymers are investigated within the tenets of the optimized Rouse-Zimm theory. The distance dependent excluded volume interactions (EVIs) are obtained from Flory's mean field theory. The hydrodynamic interactions (HIs) between the pairs of monomers are estimated using the preaveraged Oseen tensor. The intrinsic viscosity of linear and ring polymers both with and without EVI are compared as a function of ring size. A monotonically increasing trend of the intrinsic viscosity is observed in both cases. The intrinsic viscosity of both linear and ring polymers both with and without EVI show a very good agreement with the experimental results of polystyrene over a wide range of molecular weights in both good and theta solvents, respectively. The fractal dimensions of the ring polymers with EVI lie between that of a random walk and a self-avoiding walk model of linear polymers in three dimensions. The ring size increases with EVI and the effect of EVI is stronger on larger rings than that on smaller rings. The dielectric relaxation follow a connectivity independent universal scaling behavior at low and high frequency regions. The imaginary part of the complex dielectric susceptibility displays a local maxima in the intermediate frequency region, which reveals a structure dependent behavior of the rings. The theoretically calculated dielectric loss of ring polymers with HI matches well with those obtained from experiments.

Curvature induced structural changes of the chicken villin headpiece subdomain by single walled carbon nanotubes.

Carbon nanotubes (CNTs) are identified as potential candidates for drug and biomolecular loading and delivery. CNTs of different chiralities have different diameters, which may significantly affect their abilities to interact with different types of biomolecules. Herein, we employ classical molecular dynamics simulation to provide insight into the curvature-dependent interactions between a model protein, chicken villin headpiece subdomain (HP36), with CNTs having chiralities (8,8), (12,12), and (20,20). It is revealed that, with increasing radii, the protein encounters more aromatic carbon atoms on the surface of the CNT, leading to its increasing strength of adsorption. However, the extent of adsorption has a limiting magnitude, after which an increase in the radius of the nanotube has practically no effect on the extent of adsorption. Spontaneous encapsulation of the protein was demonstrated using a (28,28) CNT, where the protein is found to undergo insignificant structural perturbation. Finally, steered molecular dynamics simulations have been performed to mimic the force-induced release of the protein from within the nanotube cavity. It has been identified that a minimum force of ∼300 pN and a minimum velocity of 4 Šns-1 are required to release the protein from the CNT at 300 K. Any external force below the critical magnitude and inducing velocity less than 4 Šns-1 allows the translocation of the protein through the inner surface of the CNT; however, before being released, the protein undergoes unfolding, thereby losing the secondary structure and biological activity.

Exploring the aggregation of amyloid-ß 42 through Monte Carlo simulations.

Coarse-grained Monte Carlo simulations are performed for a disordered protein, amyloid-ß 42 to identify the interactions and understand the mechanism of its aggregation. A statistical potential is developed from a selected dataset of intrinsically disordered proteins, which accounts for the respective contributions of the bonded and non-bonded potentials. While, the bonded potential comprises the bond, bend, and dihedral constraints, the nonbonded interactions include van der Waals interactions, hydrogen bonds, and the two-body potential. The two-body potential captures the features of both hydrophobic and electrostatic interactions that brings the chains at a contact distance, while the repulsive van der Waals interactions prevent them from a collapse. Increased two-body hydrophobic interactions facilitate the formation of amorphous aggregates rather than the fibrillar ones. The formation of aggregates is validated from the interchain distances, and the total energy of the system. The aggregate is structurally characterized by the root-mean-square deviation, root-mean-square fluctuation and the radius of gyration. The aggregates are characterized by a decrease in SASA, an increase in the non-local interactions and a distinct free energy minimum relative to that of the monomeric state of amyloid-ß 42. The hydrophobic residues help in nucleation, while the charged residues help in oligomerization and aggregation.

Relaxation dynamics measure the aggregation propensity of amyloid-ß and its mutants.

Atomistic molecular dynamics simulations are employed to investigate the global and segmental relaxation dynamics of the amyloid-ß protein and its causative and protective mutants. Amyloid-ß exhibits significant global/local dynamics that span a broad range of length and time scales due to its intrinsically disordered nature. The relaxation dynamics of the amyloid-ß protein and its mutants is quantitatively correlated with its experimentally measured aggregation propensity. The protective mutant has slower relaxation dynamics, whereas the causative mutants exhibit faster global dynamics compared with that of the wild-type amyloid-ß. The local dynamics of the amyloid-ß protein or its mutants is governed by a complex interplay of the charge, hydrophobicity, and change in the molecular mass of the mutated residue.

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 the Passage Times for Unfolding/Folding of the Adenine Riboswitch Aptamer.

The conformational transitions of the adenosine deaminase A-riboswitch aptamer both with and without ligand binding are investigated within the tenets of the generalized Langevin equation in a complex viscoelastic cellular environment. Steered molecular dynamics (SMD) simulations are performed to evaluate and compare the results of the first passage times (FPTs) with those obtained from the theory for the unfold and fold transitions of the aptamer. The results of the distribution of Kramers's FPT reveal that the unfold-fold transitions are faster and hence more probable as compared to the fold-unfold transitions of the riboswitch aptamer for both ligand-bound and -unbound states. The transition path time is lower than Kramers's FPT for the riboswitch aptamer as the transition path times for the unfold-fold transition of both without and with ligand binding are insensitive to the details of the exact mechanism of the transition events. However, Kramers's FPTs show varied distributions which correspond to different transition pathways, unlike the transition path times. The mean FPT increases with an increase in the complexity of the cellular environment. The results of Kramers's FPT, transition path time distribution, and mean FPT obtained from our calculations qualitatively match with those obtained from the SMD simulations. Analytically derived values of the mean transition path time show good quantitative agreement with those estimated from the single-molecule force spectroscopy experiments for higher barrier heights.

Effect of ligand binding on riboswitch folding: Theory and simulations.

The effect of ligand binding on the conformational transitions of the add A-riboswitch in cellular environments is investigated theoretically within the framework of the generalized Langevin equation combined with steered molecular dynamics simulations. Results for the transition path time distribution provide an estimate of the transit times, which are difficult to determine experimentally. The time for the conformational transitions of the riboswitch aptamer is longer for the ligand bound state as compared to that of the unbound one. The transition path time of the riboswitch follows a counterintuitive trend as it decreases with an increase in the barrier height. The mean transition path time of either transitions of the riboswitch in the ligand bound/unbound state increases with an increase in the complexity of the surrounding environment due to the caging effect. The results of the probability density function, transition path time distribution, and mean transition path time obtained from the theory qualitatively agree with those obtained from the simulations and with earlier experimental and theoretical studies.

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.

Conformational properties of complexes of poly(propylene imine) dendrimers with linear polyelectrolytes in dilute solutions.

This study investigates the conformational properties of complexes of poly(propylene imine) dendrimers with a linear polyelectrolyte (LPE) at neutral pH in an aqueous solution via molecular dynamics simulations. Various conformational properties, such as the atomic density profile, counterion density distribution, charge distribution, cavity volume, and the static structure factor are studied as a function of the charge and chain length of the LPE. The lower generation dendrimer complexes encapsulate the shorter linear PE chains, while the longer PE chains are adsorbed on the dendrimer surface that screen the surface charge and prevent the penetration of the counterions and water molecules. However, the overall charge of the higher generation dendrimers is not neutralized by the charge of the PE chains, which results in chloride counterion penetration within the dendrimers. The adsorption of the PE chains on the dendrimers is also verified from the charge distribution of the dendrimer-PE complexes. The charge on the lower generation dendrimer complexes is overcompensated by the longer PE chains resulting in an overall negative charge on the complexes, while the PE chains do not completely neutralize the charge of the higher generation dendrimers and produce positively charged complexes. The results of the structure factor indicate a conformational transition of the dendrimer-PE complexes from a dense compact structure to an open one with an increase in the PE chain length. This transition is characterized by an increase in the cavity volume in dendrimers with an increase in the PE chain length.

Conformations of poly(propylene imine) dendrimers in an ionic liquid at different pH.

The conformational behavior of poly(propylene imine) (PPI) dendrimers at three different solution pH is studied in an ionic liquid (IL) solvent, 1-butyl-3-methylimidazolium chloride ([BMIM]Cl), through molecular dynamics (MD) simulations. The size, shape, density distribution, structure factor, and the scattering intensity are evaluated to probe the conformational transition of dendrimers as a function of pH. The results of the radial atomic and terminal amine group density distribution at low pH indicate a shift in the density towards the periphery of the dendrimers due to the electrostatic repulsion between the charged tertiary amine groups within the dendrimers. The [BMIM] cations are not encapsulated within dendrimers and predominantly reside near the periphery. The extensive back-folding of the outer branches due to the electrostatic repulsion between the solvent cations and the peripheral charged amine groups at neutral and low pH results in a dense compact structure in [BMIM]Cl as compared to that in water, as evident from the results of the structure factor and scattering intensity. The structural analysis in terms of the fractal dimension reveals that the lower generation dendrimers exhibit conformational transition as a function of pH, while the higher generations exhibit a highly compact structure at all solution pH. However, PPI dendrimers at low pH exhibit more free volume as compared to that at high pH, which may be utilized to accommodate specific guest molecules.

Interaction Volume Is a Measure of the Aggregation Propensity of Amyloid-ß.

This study highlights the significance of the partial molar volume of amino acids in predicting the aggregation propensity of an intrinsically disordered protein, amyloid-ß (Aß), and its mutants in aqueous solution. The change in the interaction volume of the protein or mutant is quantitatively correlated with its calculated experimental aggregation propensity. This method also reveals how the interaction volume may be tuned by changing the charge and hydrophobicity of Aß. While a positive change in the interaction volume and a higher aggregation propensity are observed for mutants with a decrease in the overall charge and/or an increase in hydrophobicity, a reverse trend is observed for the mutants with a decrease in the hydrophobicity and/or an increase in its charge. Hence, the interaction volume may be considered as a key parameter for monitoring protein aggregation that bridges the gap between the experimental aggregation kinetics and solvation thermodynamics.

Orientational Relaxation of Poly(propylene imine) Dendrimers at Different pH.

The dilute solution dynamics of poly(propylene imine) (PPI) dendrimers is investigated at three different solution pH through molecular dynamics (MD) simulations. The dynamics of PPI dendrimers is characterized by both global and local relaxations that occur at different time and length scales. While the global dynamics may be described in terms of rotational diffusion, the local motion may be characterized through orientational relaxation dynamics measured in terms of the time autocorrelation function (ACF), second-order orientational ACF, and the spin-lattice relaxation rate. The global motion of dendrimers decreases with an increase in the size from high pH to low pH with increasing generations of growth. The results reveal that the segments at low pH relax faster than those at high pH, and the local mobility of the segments near the periphery is higher than the core segments. This observation is also evident from the spectral density and spin-lattice relaxation rate. High values of the spectral density at higher frequencies imply higher segmental mobility of the dendrimer at low pH relative to that at high pH. A shift in the maximum of the spin-lattice relaxation rate toward lower frequencies with decreasing generations indicates the dependence of local mobility on the topological distance of the segment from the periphery at all pH conditions.

Effect of Alzheimer's Disease Causative and Protective Mutations on the Hydration Environment of Amyloid-ß.

This study investigates the local structure and dynamics of hydration water around the intrinsically disordered protein amyloid-ß (Aß) and its Alzheimer's disease causative N-terminus mutants, i.e., A2V, Taiwan (D7H), Tottori (D7N), and English (H6R), and the protective A2T mutant via atomistic MD simulations. The effect of mutations on the hydration environment around different domains of this protein is evaluated through the surface distribution function, tetrahedral order parameter, and the survival probability of the water molecules within the hydration shell. The water density around the hydrophobic hp1 (17-21) domain is found to be higher for the A2T mutant as compared to the wild-type Aß and its selected causative mutants. The average tetrahedral order parameter of water molecules around the hydrophobic hp1 (17-21) domain shows that water molecules are less ordered around the A2V, Taiwan, Tottori, and English mutants and more ordered around the A2T mutant than those of the wild-type protein. The survival probability decays rapidly for the A2V, Taiwan, Tottori, and English mutants, while it is comparatively slower for the protective A2T mutant. These results exhibit different hydration environments for the causative and protective mutants highlighting the differential roles of charge and hydrophobicity.

Orientational relaxation of ring polymers in dilute solutions.

The segmental relaxation dynamics of ring polymers in dilute solutions is investigated via optimized Rouse-Zimm theory. To the best of our knowledge, this is the first study that characterizes the orientational relaxation dynamics of ring polymers in dilute solutions. The orientational time autocorrelation functions are governed by two major processes that span a broad range of timescales: (i) local segmental motion at short times, independent of the ring size, and (ii) overall motion of the ring at long times that depends on the limiting ring size. Smaller rings relax faster than larger rings and their respective linear analogues. The hydrodynamic interactions decrease the higher relaxation rates corresponding to the local relaxation modes and increase the smaller relaxation rates which correspond to the collective relaxation modes. The spectral density is independent of frequency in the low frequency regime while it decreases with increasing frequency. Regardless of the ring size, the spin-lattice relaxation rate exhibits a single characteristic maximum as a function of frequency that shifts to a lower value with increasing strength of hydrodynamic interactions.

Effect of Correlated Pair Mutations in Protein Misfolding.

A Monte Carlo simulation based sequence design method is proposed to explore the effect of correlated pair mutations in proteins. In the designed sequences, the most correlated residue pairs are identified and mutated with all possible amino acid pairs except those already present. The cumulative correlated pair mutations generated an array of mutated sequences. Results show a significant increase in the probability of misfolding for correlated pair mutations as compared to that of the random pair mutations. The pair mutations of correlated residues that are in contact record a higher probability of misfolding as compared to the correlated residues that are not in contact. The probability of misfolding increases on pair mutation of nonlocally correlated residue pairs as compared to that of the locally correlated residue pairs. The choice of a compact or expanded conformation does not depend on the type of correlated pair mutations. Pair mutation of the most correlated residue pairs at the surface with hydrophobic amino acids results in higher misfolding probability as compared to that in the core. An exactly opposite behavior is observed on pair mutation with hydrophilic and charged amino acid pairs. The neutral amino acid pairs do not differentiate between core and surface sites. This study may be used for targeted mutation experiments to predict complex mutation patterns, reengineer the existing proteins, and design new proteins with reduced misfolding propensity.