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

Exploring Structural Insights of Aβ42 and α-Synuclein Monomers and Heterodimer: A Comparative Study Using Implicit and Explicit Solvent Simulations

Protein misfolding, aggregation, and fibril formation play a central role in the development of severe neurological disorders, including Alzheimer's and Parkinson's diseases. The structural stability of mature fibrils in these diseases is of great importance, as organisms struggle to effectively eliminate amyloid plaques. To address this issue, it is crucial to investigate the early stages of fibril formation when monomers aggregate into small, toxic, and soluble oligomers. However, these structures are inherently disordered, making them challenging to study through experimental approaches. Recently, it has been shown experimentally that amyloid-β 42 (Aβ42) and α-synuclein (α-Syn) can coassemble. This has motivated us to investigate the interaction between their monomers as a first step toward exploring the possibility of forming heterodimeric complexes. In particular, our study involves the utilization of various Amber and CHARMM force-fields, employing both implicit and explicit solvent models in replica exchange and conventional simulation modes. This comprehensive approach allowed us to assess the strengths and weaknesses of these solvent models and force fields in comparison to experimental and theoretical findings, ensuring the highest level of robustness. Our investigations revealed that Aβ42 and α-Syn monomers can indeed form stable heterodimers, and the resulting heterodimeric model exhibits stronger interactions compared to the Aβ42 dimer. The binding of α-Syn to Aβ42 reduces the propensity of Aβ42 to adopt fibril-prone conformations and induces significant changes in its conformational properties. Notably, in AMBER-FB15 and CHARMM36m force fields with the use of explicit solvent, the presence of Aβ42 significantly increases the β-content of α-Syn, consistent with the experiments showing that Aβ42 triggers α-Syn aggregation. Our analysis clearly shows that although the use of implicit solvent resulted in too large compactness of monomeric α-Syn, structural properties of monomeric Aβ42 and the heterodimer were preserved in explicit-solvent simulations. We anticipate that our study sheds light on the interaction between α-Syn and Aβ42 proteins, thus providing the atom-level model required to assess the initial stage of aggregation mechanisms related to Alzheimer's and Parkinson's diseases.

Modeling structural interconversion in Alzheimers' amyloid beta peptide with classical and intrinsically disordered protein force fields

A comprehensive understanding of the aggregation mechanism in amyloid beta 42 (Aβ42) peptide is imperative for developing therapeutic drugs to prevent or treat Alzheimer's disease. Because of the high flexibility and lack of native tertiary structures of Aβ42, molecular dynamics (MD) simulations may help elucidate the peptide's dynamics with atomic details and collectively improve ensembles not seen in experiments. We applied microsecond-timescale MD simulations to investigate the dynamics and conformational changes of Aβ42 by using a newly developed Amber force field (ff14IDPSFF). We compared the ff14IDPSFF and the regular ff14SB force field by examining the conformational changes of two distinct Aβ42 monomers in explicit solvent. Conformational ensembles obtained by simulations depend on the force field and initial structure, Aβ42^α-helix or Aβ42^β-strand. The ff14IDPSFF sampled a high ratio of disordered structures and diverse β-strand secondary structures; in contrast, ff14SB favored helicity during the Aβ42^α-helix simulations. The conformations obtained from Aβ42^β-strand simulations maintained a balanced content in the disordered and helical structures when simulated by ff14SB, but the conformers clearly favored disordered and β-sheet structures simulated by ff14IDPSFF. The results obtained with ff14IDPSFF qualitatively reproduced the NMR chemical shifts well. In-depth peptide and cluster analysis revealed some characteristic features that may be linked to early onset of the fibril-like structure. The C-terminal region (mainly M35-V40) featured in-registered anti-parallel β-strand (β-hairpin) conformations with tested systems. Our work should expand the knowledge of force field and structure dependency in MD simulations and reveals the underlying structural mechanism-function relationship in Aβ42 peptides. Communicated by Ramaswamy H. Sarma.

Monitoring early-stage β-amyloid dimer aggregation by histidine site-specific two-dimensional infrared spectroscopy in a simulation study

Monitoring early-stage β-amyloid (Aβ) dimerization is a formidable challenge for understanding neurological diseases. We compared β-sheet formation and histidine site-specific two-dimensional infrared (2D IR) spectroscopic signatures of Aβ dimers with different histidine states (δ; N^δ1-H, ε; N^ε2-H, or π; both protonated). Molecular dynamics (MD) simulations revealed that β-sheet formation is favored for the δδδ:δδδ and πππ:πππ tautomeric isomers showing strong couplings and frequent contacts between the central hydrophobic core and C-terminus compared with the εεε:εεε isomer. Characteristic blue-shifts in the 2D IR central bands were observed upon monomer-dimer transformation. The εεε:εεε dimer exhibited larger frequency shifts than δδδ:δδδ and πππ:πππ implying that the red-shift may have a correlation with N^δ1-H(δ) protonation. Our results support the tautomerization/protonation hypothesis that attributes Aβ misfolding to histidine tautomers as a possible primary initiator for Aβ aggregation and facilitates the application of histidine site-specific 2D IR spectroscopy for studying early-stage Aβ self-assembly.

Structural and Molecular Dynamics of Mycobacterium tuberculosis Malic Enzyme, a Potential Anti-TB Drug Target

Tuberculosis (TB) is the most lethal bacterial infectious disease worldwide. It is notoriously difficult to treat, requiring a cocktail of antibiotics administered over many months. The dense, waxy outer membrane of the TB-causing agent, Mycobacterium tuberculosis (Mtb), acts as a formidable barrier against uptake of antibiotics. Subsequently, enzymes involved in maintaining the integrity of the Mtb cell wall are promising drug targets. Recently, we demonstrated that Mtb lacking malic enzyme (MEZ) has altered cell wall lipid composition and attenuated uptake by macrophages. These results suggest that MEZ contributes to lipid biosynthesis by providing reductants in the form of NAD(P)H. Here, we present the X-ray crystal structure of MEZ to 3.6 Å. We use biochemical assays to demonstrate MEZ is dimeric in solution and to evaluate the effects of pH and allosteric regulators on its kinetics and thermal stability. To assess the interactions between MEZ and its substrate malate and cofactors, Mn²⁺ and NAD(P)⁺, we ran a series of molecular dynamics (MD) simulations. First, the MD analysis corroborates our empirical observations that MEZ is unusually flexible, which persists even with the addition of substrate and cofactors. Second, the MD simulations reveal that dimeric MEZ subunits alternate between open and closed states, and that MEZ can stably bind its NAD(P)⁺ cofactor in multiple conformations, including an inactive, compact NAD⁺ form. Together the structure of MEZ and insights from its dynamics can be harnessed to inform the design of MEZ inhibitors that target Mtb and not human malic enzyme homologues.

An α-helix mimetic oligopyridylamide, ADH-31, modulates Aβ42 monomer aggregation and destabilizes protofibril structures: insights from molecular dynamics simulations

Alzheimer's disease (AD), an epidemic growing worldwide due to no effective medical aid available in the market, is a neurological disorder. AD is known to be directly associated with the toxicity of amyloid-β (Aβ) aggregates. In search of potent inhibitors of Aβ aggregation, Hamilton and co-workers reported an α-helix mimetic, ADH-31, which acts as a powerful antagonist of Aβ42 aggregation. To identify the key interactions between protein-ligand complexes and to gain insights into the inhibitory mechanism of ADH-31 against Aβ42 aggregation, molecular dynamics (MD) simulations were performed in the present study. The MD simulations highlighted that ADH-31 showed distinct binding capabilities with residues spanning from the N-terminal to the central hydrophobic core (CHC) region of Aβ42 and restricted the conformational transition of the helix-rich structure of Aβ42 into another form of secondary structures (coil/turn/β-sheet). Hydrophobic contacts, hydrogen bonding and π-π interaction contribute to the strong binding between ADH-31 and Aβ42 monomer. The Dictionary of Secondary Structure of Proteins (DSSP) analysis highlighted that the probability of helical content increases from 38.5% to 50.2% and the turn content reduces from 14.7% to 6.2% with almost complete loss of the β-sheet structure (4.5% to 0%) in the Aβ42 monomer + ADH-31 complex. The per-residue binding free energy analysis demonstrated that Arg5, Tyr10, His14, Gln15, Lys16, Val18, Phe19 and Lys28 residues of Aβ42 are responsible for the favourable binding free energy in Aβ42 monomer + ADH-31 complex, which is consistent with the 2D HSQC NMR of the Aβ42 monomer that depicted a change in the chemical shift of residues spanning from Glu11 to Phe20 in the presence of ADH-31. The MD simulations highlighted the prevention of sampling of amyloidogenic β-strand conformations in Aβ42 trimer in the presence of ADH-31 as well as the ability of ADH-31 to destabilize Aβ42 trimer and protofibril structures. The lower binding affinity between Aβ42 trimer chains in the presence of ADH-31 highlights the destabilization of the Aβ42 trimer structure. Overall, MD results highlighted that ADH-31 inhibited Aβ42 aggregation by constraining Aβ peptides into helical conformation and destabilized Aβ42 trimer as well as protofibril structures. The present study provides a theoretical insight into the atomic level details of the inhibitory mechanism of ADH-31 against Aβ42 aggregation as well as protofibril destabilization and could be implemented in the structure-based drug design of potent therapeutic agents for AD.

Spontaneous self-assembly of amyloid β (1-40) into dimers

The self-assembly and fibrillation of amyloid β (Aβ) proteins is the neuropathological hallmark of Alzheimer's disease. However, the molecular mechanism of how disordered monomers assemble into aggregates remains largely unknown. In this work, we characterize the assembly of Aβ (1-40) monomers into dimers using long-time molecular dynamics simulations. Upon interaction, the monomers undergo conformational transitions, accompanied by change of the structure, leading to the formation of a stable dimer. The dimers are stabilized by interactions in the N-terminal region (residues 5-12), in the central hydrophobic region (residues 16-23), and in the C-terminal region (residues 30-40); with inter-peptide interactions focused around the N- and C-termini. The dimers do not contain long β-strands that are usually found in fibrils.

Replica exchange molecular dynamics simulation of the coordination of Pt(ii)-Phenanthroline to amyloid-β

We report replica exchange molecular dynamics (REMD) simulations of the complex formed between amyloid-β peptides and platinum bound to a phenanthroline ligand, Pt(phen). After construction of an AMBER-style forcefield for the Pt complex, REMD simulation employing temperatures between 270 and 615 K was used to provide thorough sampling of the conformational freedom available to the peptide. We find that the full length peptide Aβ42, in particular, frequently adopts a compact conformation with a large proportion of α- and 3,10-helix content, with smaller amounts of β-strand in the C-terminal region of the peptide. Helical structures are more prevalent than in the metal-free peptide, while turn and strand conformations are markedly less common. Non-covalent interactions, including salt-bridges, hydrogen bonds, and π-stacking between aromatic residues and the phenanthroline ligand, are common, and markedly different from those seen in the amyloid-β peptides alone.

Replica exchange molecular dynamics are used to explore the conformational freedom of amyloid-βbound to Pt(phenanthroline), highlighting important differences in secondary and tertiary structure from the metal-free peptide.

Perspective: Identification of collective variables and metastable states of protein dynamics

The statistical analysis of molecular dynamics simulations requires dimensionality reduction techniques, which yield a low-dimensional set of collective variables (CVs) {x _i } = x that in some sense describe the essential dynamics of the system. Considering the distribution P( x ) of the CVs, the primal goal of a statistical analysis is to detect the characteristic features of P( x ), in particular, its maxima and their connection paths. This is because these features characterize the low-energy regions and the energy barriers of the corresponding free energy landscape ΔG( x ) = -k _B T ln P( x ), and therefore amount to the metastable states and transition regions of the system. In this perspective, we outline a systematic strategy to identify CVs and metastable states, which subsequently can be employed to construct a Langevin or a Markov state model of the dynamics. In particular, we account for the still limited sampling typically achieved by molecular dynamics simulations, which in practice seriously limits the applicability of theories (e.g., assuming ergodicity) and black-box software tools (e.g., using redundant input coordinates). We show that it is essential to use internal (rather than Cartesian) input coordinates, employ dimensionality reduction methods that avoid rescaling errors (such as principal component analysis), and perform density based (rather than k-means-type) clustering. Finally, we briefly discuss a machine learning approach to dimensionality reduction, which highlights the essential internal coordinates of a system and may reveal hidden reaction mechanisms.

Principal component analysis on a torus: Theory and application to protein dynamics

A dimensionality reduction method for high-dimensional circular data is developed, which is based on a principal component analysis (PCA) of data points on a torus. Adopting a geometrical view of PCA, various distance measures on a torus are introduced and the associated problem of projecting data onto the principal subspaces is discussed. The main idea is that the (periodicity-induced) projection error can be minimized by transforming the data such that the maximal gap of the sampling is shifted to the periodic boundary. In a second step, the covariance matrix and its eigendecomposition can be computed in a standard manner. Adopting molecular dynamics simulations of two well-established biomolecular systems (Aib₉ and villin headpiece), the potential of the method to analyze the dynamics of backbone dihedral angles is demonstrated. The new approach allows for a robust and well-defined construction of metastable states and provides low-dimensional reaction coordinates that accurately describe the free energy landscape. Moreover, it offers a direct interpretation of covariances and principal components in terms of the angular variables. Apart from its application to PCA, the method of maximal gap shifting is general and can be applied to any other dimensionality reduction method for circular data.

Effects of force fields on the conformational and dynamic properties of amyloid β(1-40) dimer explored by replica exchange molecular dynamics simulations

The conformational space and structural ensembles of amyloid beta (Aβ) peptides and their oligomers in solution are inherently disordered and proven to be challenging to study. Optimum force field selection for molecular dynamics (MD) simulations and the biophysical relevance of results are still unknown. We compared the conformational space of the Aβ(1-40) dimers by 300 ns replica exchange MD simulations at physiological temperature (310 K) using: the AMBER-ff99sb-ILDN, AMBER-ff99sb*-ILDN, AMBER-ff99sb-NMR, and CHARMM22* force fields. Statistical comparisons of simulation results to experimental data and previously published simulations utilizing the CHARMM22* and CHARMM36 force fields were performed. All force fields yield sampled ensembles of conformations with collision cross sectional areas for the dimer that are statistically significantly larger than experimental results. All force fields, with the exception of AMBER-ff99sb-ILDN (8.8 ± 6.4%) and CHARMM36 (2.7 ± 4.2%), tend to overestimate the α-helical content compared to experimental CD (5.3 ± 5.2%). Using the AMBER-ff99sb-NMR force field resulted in the greatest degree of variance (41.3 ± 12.9%). Except for the AMBER-ff99sb-NMR force field, the others tended to under estimate the expected amount of β-sheet and over estimate the amount of turn/bend/random coil conformations. All force fields, with the exception AMBER-ff99sb-NMR, reproduce a theoretically expected β-sheet-turn-β-sheet conformational motif, however, only the CHARMM22* and CHARMM36 force fields yield results compatible with collapse of the central and C-terminal hydrophobic cores from residues 17-21 and 30-36. Although analyses of essential subspace sampling showed only minor variations between force fields, secondary structures of lowest energy conformers are different.

Bending energy penalty enhances the adhesive strength of functional amyloid curli to surfaces

The functional amyloid curli fiber, a major proteinaceous component of biofilm extracellular matrices, plays an important role in biofilm formation and enterobacteriaceae adhesion. Curli nanofibers exhibit exceptional underwater adhesion to various surfaces, have high rigidity and strong tensile mechanical properties, and thus hold great promise in biomaterials. The mechanisms of how curli fibers strongly attach to surfaces and detach under force remain elusive. To investigate curli fiber adhesion to surfaces, we developed a coarse-grained curli fiber model, in which the protein subunit CsgA (curli specific gene A) self-assembles into the fiber. The coarse-grained model yields physiologically relevant and tunable bending rigidity and persistence length. The force-induced desorption of a single curli fiber is examined using coarse-grained modeling and theoretical analysis. We find that the bending energy penalty arising from high persistence length enhances the resistance of the curli fiber against desorption and thus strengthens the adhesion of the curli fiber to surfaces. The CsgA-surface adhesion energy and the curli fiber bending rigidity both play crucial roles in the resistance of curli fiber against desorption from surfaces. To enable the desorption process, the applied peeling force must overcome both the interfacial adhesion energy and the energy barrier for bending the curli fiber at the peeling front. We show that the energy barrier to desorption increases with the interfacial adhesion energy, however, the bending induced failure of a single curli fiber limits the work of adhesion if the proportion of the CsgA-surface adhesion energy to the CsgA-CsgA cohesive energy becomes large. These results illustrate that the optimal adhesion performance of nanofibers is dictated by the interplay between bending, surface energy and cohesive energy. Our model provides timely insight into enterobacteriaceae adhesion mechanisms as well as future designs of engineered curli fiber based adhesives.

Major Reaction Coordinates Linking Transient Amyloid-β Oligomers to Fibrils Measured at Atomic Level

The structural underpinnings for the higher toxicity of the oligomeric intermediates of amyloidogenic peptides, compared to the mature fibrils, remain unknown at present. The transient nature and heterogeneity of the oligomers make it difficult to follow their structure. Here, using vibrational and solid-state nuclear magnetic resonance spectroscopy, and molecular dynamics simulations, we show that freely aggregating Aβ₄₀ oligomers in physiological solutions have an intramolecular antiparallel configuration that is distinct from the intermolecular parallel β-sheet structure observed in mature fibrils. The intramolecular hydrogen-bonding network flips nearly 90°, and the two β-strands of each monomeric unit move apart, to give rise to the well-known intermolecular in-register parallel β-sheet structure in the mature fibrils. Solid-state nuclear magnetic resonance distance measurements capture the interstrand separation within monomer units during the transition from the oligomer to the fibril form. We further find that the D23-K28 salt-bridge, a major feature of the Aβ₄₀ fibrils and a focal point of mutations linked to early onset Alzheimer's disease, is not detectable in the small oligomers. Molecular dynamics simulations capture the correlation between changes in the D23-K28 distance and the flipping of the monomer secondary structure between antiparallel and parallel β-sheet architectures. Overall, we propose interstrand separation and salt-bridge formation as key reaction coordinates describing the structural transition of the small Aβ₄₀ oligomers to fibrils.