Orkid Coskuner

Amyloid-β peptide structure in aqueous solution varies with fragment size

Various fragment sizes of the amyloid-β (Aβ) peptide have been utilized to mimic the properties of the full-length Aβ peptide in solution. Among these smaller fragments, Aβ16 and Aβ28 have been investigated extensively. In this work, we report the structural and thermodynamic properties of the Aβ16, Aβ28, and Aβ42 peptides in an aqueous solution environment. We performed replica exchange molecular dynamics simulations along with thermodynamic calculations for investigating the conformational free energies, secondary and tertiary structures of the Aβ16, Aβ28, and Aβ42 peptides. The results show that the thermodynamic properties vary from each other for these peptides. Furthermore, the secondary structures in the Asp1-Lys16 and Asp1-Lys28 regions of Aβ42 cannot be completely captured by the Aβ16 and Aβ28 fragments. For example, the β-sheet structures in the N-terminal region of Aβ16 and Aβ28 are either not present or the abundance is significantly decreased in Aβ42. The α-helix and β-sheet abundances in Aβ28 and Aβ42 show trends--to some extent--with the potential of mean forces but no such trend could be obtained for Aβ16. Interestingly, Arg5 forms salt bridges with large abundances in all three peptides. The formation of a salt bridge between Asp23-Lys28 is more preferred over the Glu22-Lys28 salt bridge in Aβ28 but this trend is vice versa for Aβ42. This study shows that the Asp1-Lys16 and Asp1-Lys28 regions of the full length Aβ42 peptide cannot be completely mimicked by studying the Aβ16 and Aβ28 peptides.

Structures of the E46K Mutant-Type α-Synuclein Protein and Impact of E46K Mutation on the Structures of the Wild-Type α-Synuclein Protein

The E46K genetic missense mutation of the wild-type α-synuclein protein was recently identified in a family of Spanish origin with hereditary Parkinsons disease. Detailed understanding of the structures of the monomeric E46K mutant-type α-synuclein protein as well as the impact of the E46K missense mutation on the conformations and free energy landscapes of the wild-type α-synuclein are required for gaining insights into the pathogenic mechanism of Parkinsons disease. In this study, we use extensive parallel tempering molecular dynamics simulations along with thermodynamic calculations to assess the secondary and tertiary structural properties as well as the conformational preferences of the monomeric wild-type and E46K mutant-type α-synuclein proteins in an aqueous solution environment. We also present the residual secondary structure component conversion stabilities with dynamics using a theoretical strategy, which we most recently developed. To the best of our knowledge, this study presents the first detailed comparison of the structural and thermodynamic properties of the wild-type and E46K mutant-type α-synuclein proteins in an aqueous solution environment at the atomic level with dynamics. We find that the E46K mutation results not only in local but also in long-range changes in the structural properties of the wild-type α-synuclein protein. The mutation site shows a significant decrease in helical content as well as a large increase in β-sheet structure formation upon E46K mutation. In addition, the β-sheet content of the C-terminal region increases significantly in the E46K mutant-type αS in comparison to the wild-type αS. Our theoretical strategy developed to assess the thermodynamic preference of secondary structure transitions indicates that this shift in secondary structure is the result of a decrease in the thermodynamic preference of turn to helix conversions while the coil to β-sheet preference increases for these residues. Long-range intramolecular protein interactions of the C-terminal with the N-terminal and NAC regions increase upon E46K mutation, resulting in more compact structures for the E46K mutant-type rather than wild-type αS. However, the E46K mutant-type αS structures are less stable than the wild-type αS. Overall, our results show that the E46K mutant-type αS has a higher propensity to aggregate than the wild-type αS and that the N-terminal and C-terminal regions are reactive toward fibrillization and aggregation upon E46K mutation and we explain the associated reasons based on the structural properties herein. Small molecules or drugs that can block the specific residues forming abundant β-sheet structure, which we report here, might help to reduce the reactivity of these intrinsically disordered fibrillogenic proteins toward aggregation and their toxicity.

Structures and Free Energy Landscapes of the A53T Mutant-Type α-Synuclein Protein and Impact of A53T Mutation on the Structures of the Wild-Type α-Synuclein Protein with Dynamics

The A53T genetic missense mutation of the wild-type α-synuclein (αS) protein was initially identified in Greek and Italian families with familial Parkinsons disease. Detailed understanding of the structures and the changes induced in the wild-type αS structure by the A53T mutation, as well as establishing the direct relationships between the rapid conformational changes and free energy landscapes of these intrinsically disordered fibrillogenic proteins, helps to enhance our fundamental knowledge and to gain insights into the pathogenic mechanism of Parkinsons disease. We employed extensive parallel tempering molecular dynamics simulations along with thermodynamic calculations to determine the secondary and tertiary structural properties as well as the conformational free energy surfaces of the wild-type and A53T mutant-type αS proteins in an aqueous solution medium using both implicit and explicit water models. The confined aqueous volume effect in the simulations of disordered proteins using an explicit model for water is addressed for a model disordered protein. We also assessed the stabilities of the residual secondary structure component interconversions in αS based on free energy calculations at the atomic level with dynamics using our recently developed theoretical strategy. To the best of our knowledge, this study presents the first detailed comparison of the structural properties linked directly to the conformational free energy landscapes of the monomeric wild-type and A53T mutant-type α-synuclein proteins in an aqueous solution environment. Results demonstrate that the β-sheet structure is significantly more altered than the helical structure upon A53T mutation of the monomeric wild-type αS protein in aqueous solution. The β-sheet content close to the mutation site in the N-terminal region is more abundant while the non-amyloid-β component (NAC) and C-terminal regions show a decrease in β-sheet abundance upon A53T mutation. Obtained results utilizing our new theoretical strategy show that the residual secondary structure conversion stabilities resulting in α-helix formation are not significantly affected by the mutation. Interestingly, the residual secondary structure conversion stabilities show that secondary structure conversions resulting in β-sheet formation are influenced by the A53T mutation and the most stable residual transition yielding β-sheet occurs directly from the coil structure. Long-range interactions detected between the NAC region and the N- or C-terminal regions of the wild-type αS disappear upon A53T mutation. The A53T mutant-type αS structures are thermodynamically more stable than those of the wild-type αS protein structures in aqueous solution. Overall, the higher propensity of the A53T mutant-type αS protein to aggregate in comparison to the wild-type αS protein is related to the increased β-sheet formation and lack of strong intramolecular long-range interactions in the N-terminal region in comparison to its wild-type form. The specific residual secondary structure component stabilities reported herein provide information helpful for designing and synthesizing small organic molecules that can block the β-sheet forming residues, which are reactive toward aggregation.

The structures of the E22Δ mutant-type amyloid-β alloforms and the impact of E22Δ mutation on the structures of the wild-type amyloid-β alloforms.

Structural differences between the intrinsically disordered fibrillogenic wild-type Aβ40 and Aβ42 peptides are linked to Alzheimers disease. Recently, the E22Δ genetic missense mutation was detected in patients exhibiting Alzheimers-disease type dementia. However, detailed knowledge about the E22Δ mutant-type Aβ40 and Aβ42 alloform structures as well as the differences from the wild-type Aβ40 and Aβ42 alloform structures is currently lacking. In this study, we present the structures of the E22Δ mutant-type Aβ40 and Aβ42 alloforms as well as the impact of E22Δ mutation on the wild-type Aβ40 and Aβ42 alloform structures. For this purpose, we performed extensive microsecond-time scale parallel tempering molecular dynamics simulations coupled with thermodynamic calculations. For studying the residual secondary structure component transition stabilities, we developed and applied a new theoretical strategy in our studies. We find that the E22Δ mutant-type Aβ40 might have a higher tendency toward aggregation due to more abundant β-sheet formation in the C-terminal region in comparison to the E22Δ mutant-type Aβ42 peptide. More abundant α-helix is formed in the mid-domain regions of the E22Δ mutant-type Aβ alloforms rather than in their wild-type forms. The turn structure at Ala21-Ala30 of the wild-type Aβ, which has been linked to the aggregation process, is less abundant upon E22Δ mutation of both Aβ alloforms. Intramolecular interactions between the N-terminal and central hydrophobic core (CHC), N- and C-terminal, and CHC and C-terminal regions are less abundant or disappear in the E22Δ mutant-type Aβ alloform structures. The thermodynamic trends indicate that the wild-type Aβ42 tends to aggregate more than the wild-type Aβ40 peptide, which is in agreement with experiments. However, this trend is vice versa for the E22Δ mutant-type alloforms. The structural properties of the E22Δ mutant-type Aβ40 and Aβ42 peptides reported herein may prove useful for the development of new drugs to block the formation of toxic E22Δ mutant-type oligomers by either stabilizing helical or destabilizing β-sheet structure in the C-terminal region of these two mutant alloforms.

Coordination Studies of Al-EDTA in Aqueous Solution

The degree of aluminum toxicity is based on its complexation with organic ligands. One of these complexes is AlEDTA- (Al = aluminum, EDTA = ethylenediaminetetraacetate), the structure of which in aqueous solution has been debated on the basis of X-ray absorption and NMR measurements with different interpretations proposing different coordination. In addition, there is a lack of consensus regarding the relationship of crystalline AlEDTA- and its geometry in solution. This debate must be resolved, not merely for scientific interest, but because the use of an incorrect coordination might lead to the wrong interpretation of bioactivity and kinetics data. In this work, we predict the coordination of Al in aqueous AlEDTA- by employing ab initio calculations and Car-Parrinello molecular dynamics simulations. Our results indicate that AlEDTA- favors Al in octahedral coordination in aqueous solution. Furthermore, the predicted crystalline and solution-phase structures of AlEDTA- are similar and agree well with recent X-ray measurements, supporting the strong chelating nature of this metal-organic complex in aqueous solution.

Structures and free energy landscapes of the wild-type and A30P mutant-type α-synuclein proteins with dynamics.

The genetic missense A30P mutation of the wild-type α-synuclein protein results in the replacement of the 30th amino acid residue from alanine (Ala) to proline (Pro) and was initially found in the members of a German family who developed Parkinsons disease. Even though the structures of these proteins have been measured before, detailed understanding about the structures and their relationships with free energy landscapes is lacking, which is of interest to provide insights into the pathogenic mechanism of Parkinsons disease. We report the secondary and tertiary structures and conformational free energy landscapes of the wild-type and A30P mutant-type α-synuclein proteins in an aqueous solution environment via extensive parallel tempering molecular dynamics simulations along with thermodynamic calculations. In addition, we present the residual secondary structure component transition stabilities at the atomic level with dynamics in terms of free energy change calculations using a new strategy that we reported most recently. Our studies yield new interesting results; for instance, we find that the A30P mutation has local as well as long-range effects on the structural properties of the wild-type α-synuclein protein. The helical content at Ala18-Gly31 is less prominent in comparison to the wild-type α-synuclein protein. The β-sheet structure abundance decreases in the N-terminal region upon A30P mutation of the wild-type α-synuclein, whereas the NAC and C-terminal regions possess larger tendencies for β-sheet structure formation. Long-range intramolecular protein interactions are less abundant upon A30P mutation, especially between the NAC and C-terminal regions, which is linked to the less compact and less stable structures of the A30P mutant-type rather than the wild-type α-synuclein protein. Results including the usage of our new strategy for secondary structure transition stabilities show that the A30P mutant-type α-synuclein tendency toward aggregation is higher than the wild-type α-synuclein but we also find that the C-terminal and NAC regions of the A30P mutant-type α-synuclein are reactive toward fibrillzation and aggregation based on atomic level studies with dynamics in an aqueous solution environment. Therefore, we propose that small molecules or drugs blocking the specific residues, which we report herein, located in the NAC- and C-terminal regions of the A30P mutant-type α-synuclein protein might help to reduce the toxicity of the A30P mutant-type α-synuclein protein.

Structures and free energy landscapes of aqueous zinc(II)-bound amyloid-β(1–40) and zinc(II)-bound amyloid-β(1–42) with dynamics

Binding of divalent metal ions with intrinsically disordered fibrillogenic proteins, such as amyloid-β (Aβ), influences the aggregation process and the severity of neurodegenerative diseases. The Aβ monomers and oligomers are the building blocks of the aggregates. In this work, we report the structures and free energy landscapes of the monomeric zinc(II)-bound Aβ40 (Zn:Aβ40) and zinc(II)-bound Aβ42 (Zn:Aβ42) intrinsically disordered fibrillogenic metallopeptides in an aqueous solution by utilizing an approach that employs first principles calculations and parallel tempering molecular dynamics simulations. The structural and thermodynamic properties, including the secondary and tertiary structures and conformational Gibbs free energies of these intrinsically disordered metallopeptide alloforms, are presented. The results show distinct differing characteristics for these metallopeptides. For example, prominent β-sheet formation in the N-terminal region (Asp1, Arg5, and Tyr10) of Zn:Aβ40 is significantly decreased or lacking in Zn:Aβ42. Our findings indicate that blocking multiple reactive residues forming abundant β-sheet structure located in the central hydrophobic core and C-terminal regions of Zn:Aβ42 via antibodies or small organic molecules might help to reduce the aggregation of Zn(II)-bound Aβ42. Furthermore, we find that helix formation increases but β-sheet formation decreases in the C-terminal region upon Zn(II) binding to Aβ. This depressed β-sheet formation in the C-terminal region (Gly33–Gly38) in monomeric Zn:Aβ42 might be linked to the formation of amorphous instead of fibrillar aggregates of Zn:Aβ42.Graphical abstract

Hydrophobic Interactions by Monte Carlo Simulations

Summary The structural and thermodynamic properties of liquid water and of the dilute solutions of methane and ethane in water were calculated by Monte Carlo simulations in the temperature range 298 K to 318 K and 298 K to 333 K, respectively. The nonpolar molecules were modeled as one- and two-center Lennard–Jones particles; for the interaction potential of water a modified TIP5P model was used. The results indicate that the nonpolar solutes tend to aggregate with increasing temperature. Methane molecules preferably form water-separated pairs, even at higher temperatures, whereas for ethane contact pairs are more likely. For the thermodynamic conditions studied here, the residual chemical potential of water is a linear function of temperature.

Hydrophobic Interactions of Xenon by Monte Carlo Simulations

Hydrophobic interactions of xenon atoms dissolved in liquid water were studied by NpT Monte Carlo simulations in the temperature range 298.15 to 333K and at ambient pressure. Structural properties of dilute xenon solutions were calculated and compared to those of bulk water in order to show the influence of the hydrophobic solute. It was found that the xenon atoms tend to aggregate with increasing temperature. At low temperatures the aggregates are predominantly solvent-separated pairs; at higher temperatures the quota of contact pairs increases. Furthermore, the residual chemical potentials of xenon and water were calculated with different methods; it was found that the Widom insertion methods works best for this system. For the thermodynamic conditions of this work, the residual chemical potential of water in the presence of xenon was found to be a linear function of temperature.

Preferred conformation of the glycosidic linkage of methyl-β-mannose

The conformational preference of the glycosidic linkage of methyl-β-mannose was studied in the gas phase and in aqueous solution by ab initio calculations, and by molecular dynamics (MD) and Car-Parrinello molecular dynamics (CPMD) simulations. MD simulations were performed with various water potential functions to study the impact of the chosen water potential on the predicted conformational preference of the glycosidic linkage of the carbohydrate in solution. This study shows that the trans (t) orientation of the glycosidic linkage of methyl-β-mannose is preferred over its gauche clockwise (g+) orientation in solution. CPMD simulations clearly indicate that this preference is due to intermolecular hydrogen bonding with surrounding water molecules, whereas no such information could be demonstrated by MD simulations. This study demonstrates the importance of ab initio molecular dynamics simulations in studying the structural properties of carbohydrate-water interactions.