Stability of Iowa mutant and wild type Aβ-peptide aggregates

Impact of A2T and D23N mutations on C99 homodimer conformations

The proteolytic cleavage of C99 by γ-secretase is the last step in the production of amyloid-β (Aβ) peptides. Previous studies have shown that membrane lipid composition, cholesterol concentration, and mutation in the transmembrane helix modified the structures and fluctuations of C99. In this study, we performed atomistic molecular dynamics simulations of the homodimer of the 55-residue congener of the C-terminal domain of the amyloid protein precursor, C99(1-55), in a POPC-cholesterol lipid bilayer, and we compared the conformational ensemble of WT sequence to those of the A2T and D23N variants. These mutations are particularly interesting as the protective Alzheimer's disease (AD) A2T mutation is known to decrease Aβ production, whereas the early onset AD D23N mutation does not affect Aβ production. We found noticeable differences in the structural ensembles of the three sequences. In particular, A2T varies from both WT and D23N by having long-range effects on the population of the extracellular justamembrane helix, the interface between the G29xxx-G33xxx-G37 motifs and the fluctuations of the transmembrane helical topologies.

Pulsed Hydrogen-Deuterium Exchange Reveals Altered Structures and Mechanisms in the Aggregation of Familial Alzheimer's Disease Mutants

Mutations of the Amyloid Precursor Protein, from which the amyloid β peptide Aβ42 is cleaved, are associated with familial Alzheimer's disease. The disease-relevant familial mutations include the Arctic (E22G), Iowa (D23N), Italian (E22K), Dutch (E22Q), Japanese (D7N), English (D6R), and Flemish (A21G) variants. A detailed mechanistic understanding of the aggregation behavior of the mutant peptides at the residue level is, however, still lacking. We report here a study of the aggregation kinetics of these mutants in vitro by pulsed hydrogen-deuterium exchange mass spectrometry (HDX-MS) to obtain a temporally and sequence resolved picture of their self-assembly. For all variants, HDX occurs to give a bimodal distribution representing two soluble classes of aggregates, one protected and one solvent-exposed. There is no evidence of other classes of structural intermediates within the detection limits of the HDX approach. The fractional changes in the bimodal exchange profiles for several regions of Aβ42 reveal that the central and C-terminal peptides gain protection upon fibril formation, whereas the N-terminal regions remain largely solvent-accessible. For these mutants, all peptide fragments follow the same kinetics, acquiring solvent protection at the same time, further supporting that there are no significant populations of intermediate species under our experimental conditions. The results demonstrate the potential of pulsed HDX-MS for resolving the region-specific aggregation behavior of Aβ42 isoforms in solution where X-ray crystallography and solid-state NMR (ssNMR) are challenged.

Structural analyses and force fields comparison for NACore (68-78) and SubNACore (69-77) fibril segments of Parkinson's disease

The α-synuclein fibrils are a pathological hallmark of Parkinson's disease (PD) and are abundant in the brains of PD patients. These amyloid fibrils can aggregate into distinct polymorphism under different physical conditions. Therefore, these different fibril polymorph formations should be considered in drug design studies targeting amyloid fibrils. Recently, the atomic structures of two small fibril segments of α-synuclein, named NACore (68-78) and SubNACore (69-77), have been crystallized. These segments are critical for cytotoxicity and fibril formation. Therefore, elucidation of interface interactions between pair sheets of the NACore and SubNACore is significant for the clarification of the mechanism of fibril formation in PD. In this context, molecular dynamics (MD) simulation technique is a convenient tool to investigate interface interactions of these segments at the atomic level. However, the accuracy of these simulations depends on the utilized force fields. Therefore, we have tested the dependence of interface interactions and stabilities of these small amyloid fibrils on various force fields. From the results of triple long (100 ns) MD simulations, we inferred for the stability investigations of the NACore and SubNACore that CHARMM27 and GROMOS53A6 are the most convenient force fields whereas AMBER99SB-ILDN is the most unfavorable one. Consequently, it is expected that our findings will guide the selection of the appropriate force field for simulations between these segments and possible inhibitors of this disease.

Fundamentals of cross-seeding of amyloid proteins: an introduction

Misfolded protein aggregates formed by the same (homologous) or different (heterologous/cross) sequences are the pathological hallmarks of many protein misfolding diseases (PMDs) including Alzheimer's disease (AD) and type 2 diabetes (T2D). Different from homologous-amyloid aggregation that is solely associated with a specific PMD, cross-amyloid aggregation (i.e. cross-seeding) of different amyloid proteins is more fundamentally and biologically important for understanding and untangling not only the pathological process of each PMD, but also a potential molecular cross-talk between different PMDs. However, the cross-amyloid aggregation is still a subject poorly explored and little is known about its sequence/structure-dependent aggregation mechanisms, as compared to the widely studied homo-amyloid aggregation. Here, we review the most recent and important findings of amyloid cross-seeding behaviors from in vitro, in vivo, and in silico studies. Some typical cross-seeding phenomena between Aβ/hIAPP, Aβ/tau, Aβ/α-synuclein, and tau/α-synuclein are selected and presented, and the underlying specific or general cross-seeding mechanisms are also discussed to better reveal their sequence-structure-property relationships. The potential use of the cross-seeding concept to design amyloid inhibitors is also proposed. Finally, we offer some personal perspectives on current major challenges and future research directions in this less-studied yet important field, and hopefully this work will stimulate more research to explore all possible fundamental and practical aspects of amyloid cross-seeding.

Oncogenic G12D mutation alters local conformations and dynamics of K-Ras

K-Ras is the most frequently mutated oncoprotein in human cancers, and G12D is its most prevalent mutation. To understand how G12D mutation impacts K-Ras function, we need to understand how it alters the regulation of its dynamics. Here, we present local changes in K-Ras structure, conformation and dynamics upon G12D mutation, from long-timescale Molecular Dynamics simulations of active (GTP-bound) and inactive (GDP-bound) forms of wild-type and mutant K-Ras, with an integrated investigation of atomistic-level changes, local conformational shifts and correlated residue motions. Our results reveal that the local changes in K-Ras are specific to bound nucleotide (GTP or GDP), and we provide a structural basis for this. Specifically, we show that G12D mutation causes a shift in the population of local conformational states of K-Ras, especially in Switch-II (SII) and α3-helix regions, in favor of a conformation that is associated with a catalytically impaired state through structural changes; it also causes SII motions to anti-correlate with other regions. This detailed picture of G12D mutation effects on the local dynamic characteristics of both active and inactive protein helps enhance our understanding of local K-Ras dynamics, and can inform studies on the development of direct inhibitors towards the treatment of K-Ras^G12D-driven cancers.

Review: Precision medicine and driver mutations: Computational methods, functional assays and conformational principles for interpreting cancer drivers

At the root of the so-called precision medicine or precision oncology, which is our focus here, is the hypothesis that cancer treatment would be considerably better if therapies were guided by a tumor’s genomic alterations. This hypothesis has sparked major initiatives focusing on whole-genome and/or exome sequencing, creation of large databases, and developing tools for their statistical analyses—all aspiring to identify actionable alterations, and thus molecular targets, in a patient. At the center of the massive amount of collected sequence data is their interpretations that largely rest on statistical analysis and phenotypic observations. Statistics is vital, because it guides identification of cancer-driving alterations. However, statistics of mutations do not identify a change in protein conformation; therefore, it may not define sufficiently accurate actionable mutations, neglecting those that are rare. Among the many thematic overviews of precision oncology, this review innovates by further comprehensively including precision pharmacology, and within this framework, articulating its protein structural landscape and consequences to cellular signaling pathways. It provides the underlying physicochemical basis, thereby also opening the door to a broader community.

Molecular mechanisms for the destabilization of model membranes by islet amyloid polypeptide

Misfolding of human islet amyloid polypeptide (hIAPP) into insoluble aggregates is associated with Type 2 diabetes. It has been suggested that hIAPP toxicity may be due to its accumulation in pancreatic islets, causing membrane disruption and cell permeabilization, however the molecular basis underlying its lipid association are still unclear. Here, we combine solid-state NMR, fluorescence and bright field microscopy to investigate hIAPP - lipid membrane interactions. Real-time microscopy highlights a time-dependent penetration of hIAPP oligomers toward the most buried layers of the lipid vesicles until the membrane disrupts. Deuterium NMR was conducted on liposomes at different hIAPP concentration to probe lipid internal order and thermotropism. The gel-to-fluid phase transition of the lipids is decreased by the presence of hIAPP, and site-specific analysis of the order parameter showed a significant increase of lipid order for the first eight positions of the acyl chain, suggesting a partial insertion of the peptide inside the bilayer. These results offer experimental insight into the membrane destabilization of hIAPP on model membrane vesicles.

Conformational Dynamics and Stability of U-Shaped and S-Shaped Amyloid β Assemblies

Alzheimer’s disease is the most fatal neurodegenerative disorder characterized by the aggregation and deposition of Amyloid β (Aβ) oligomers in the brain of patients. Two principal variants of Aβ exist in humans: Aβ_1–40 and Aβ_1–42. The former is the most abundant in the plaques, while the latter is the most toxic species and forms fibrils more rapidly. Interestingly, fibrils of Aβ_1–40 peptides can only assume U-shaped conformations while Aβ_1–42 can also arrange as S-shaped three-stranded chains, as recently discovered. As alterations in protein conformational arrangement correlate with cell toxicity and speed of disease progression, it is important to characterize, at molecular level, the conformational dynamics of amyloid fibrils. In this work, Replica Exchange Molecular Dynamics simulations were carried out to compare the conformational dynamics of U-shaped and S-shaped Aβ_17–42 small fibrils. Our computational results provide support for the stability of the recently proposed S-shaped model due to the maximized interactions involving the C-terminal residues. On the other hand, the U-shaped motif is characterized by significant distortions resulting in a more disordered assembly. Outcomes of our work suggest that the molecular architecture of the protein aggregates might play a pivotal role in formation and conformational stability of the resulting fibrils.

Molecular Dynamics Simulations and Dynamic Network Analysis Reveal the Allosteric Unbinding of Monobody to H-Ras Triggered by R135K Mutation

Ras proteins, as small GTPases, mediate cell proliferation, survival and differentiation. Ras mutations have been associated with a broad spectrum of human cancers and thus targeting Ras represents a potential way forward for cancer therapy. A recently reported monobody NS1 allosterically disrupts the Ras-mediated signaling pathway, but its efficacy is reduced by R135K mutation in H-Ras. However, the detailed mechanism is unresolved. Here, using molecular dynamics (MD) simulations and dynamic network analysis, we explored the molecular mechanism for the unbinding of NS1 to H-Ras and shed light on the underlying allosteric network in H-Ras. MD simulations revealed that the overall structures of the two complexes did not change significantly, but the H-Ras–NS1 interface underwent significant conformational alteration in the mutant Binding free energy analysis showed that NS1 binding was unfavored after R135K mutation, which resulted in the unfavorable binding of NS1. Furthermore, the critical residues on H-Ras responsible for the loss of binding of NS1 were identified. Importantly, the allosteric networks for these important residues were revealed, which yielded a novel insight into the allosteric regulatory mechanism of H-Ras.

Stability differences in the NMR ensembles of amyloid β fibrils

Alzheimer’s Disease is characterized by the formation of amyloid beta (A[Formula: see text] fibril plaques in the brain. These fibrils can be probed by solid state NMR (ssNMR), which leads to an ensemble of configurations that are compatible with the NMR signals. Typically, only the lowest energy conformer is considered in computer simulations that probe the stability of fibrils and their binding with drug candidates. This restriction could produce data that are not physiologically relevant if the NMR entries differ significantly in stability or binding affinities. In order to study this effect, we have investigated the variance in stability between members of NMR ensembles. Our test cases are a patient-derived A[Formula: see text]-fibril model and two in vitro A[Formula: see text]-fibril models from a previous study we performed on comparative stability. The latter two models allow us also to compare different staggering patterns. We observe significant variations in molecular flexibility, compactness and secondary structure, suggesting that the full NMR ensemble must be considered for a physiologically relevant description of A[Formula: see text] fibrils.

Binding of ACE-inhibitors to in vitro and patient-derived amyloid-β fibril models

Currently, no drugs exist that can prevent or reverse Alzheimer's disease, a neurodegenerative disease associated with the presence, in the brain, of plaques that are composed of β-amyloid (Aβ) peptides. Recent studies suggest that angiotensin-converting enzyme (ACE) inhibitors, a set of drugs used to treat hypertension, may inhibit amyloid formation in vitro. In the present study, we investigate through computer simulations the binding of ACE inhibitors to patient-derived Aβ fibrils and contrast it with that of ACE inhibitors binding to in vitro generated fibrils. The binding affinities of the ACE inhibitors are compared with that of Congo red, a dye that is used to identify amyloid structures and that is known to be a weak inhibitor of Aβ aggregation. We find that ACE inhibitors have a lower binding affinity to the patient-derived fibrils than to in vitro generated ones. For patient-derived fibrils, their binding affinities are even lower than that of Congo red. Our observations raise doubts on the hypothesis that these drugs inhibit fibril formation in Alzheimer patients by interacting directly with the amyloids.

The Structural Basis of Oncogenic Mutations G12, G13 and Q61 in Small GTPase K-Ras4B

Ras mediates cell proliferation, survival and differentiation. Mutations in K-Ras4B are predominant at residues G12, G13 and Q61. Even though all impair GAP-assisted GTP → GDP hydrolysis, the mutation frequencies of K-Ras4B in human cancers vary. Here we aim to figure out their mechanisms and differential oncogenicity. In total, we performed 6.4 μs molecular dynamics simulations on the wild-type K-Ras4B (K-Ras4B(WT)-GTP/GDP) catalytic domain, the K-Ras4B(WT)-GTP-GAP complex, and the mutants (K-Ras4B(G12C/G12D/G12V)-GTP/GDP, K-Ras4B(G13D)-GTP/GDP, K-Ras4B(Q61H)-GTP/GDP) and their complexes with GAP. In addition, we simulated 'exchanged' nucleotide states. These comprehensive simulations reveal that in solution K-Ras4B(WT)-GTP exists in two, active and inactive, conformations. Oncogenic mutations differentially elicit an inactive-to-active conformational transition in K-Ras4B-GTP; in K-Ras4B(G12C/G12D)-GDP they expose the bound nucleotide which facilitates the GDP-to-GTP exchange. These mechanisms may help elucidate the differential mutational statistics in K-Ras4B-driven cancers. Exchanged nucleotide simulations reveal that the conformational transition is more accessible in the GTP-to-GDP than in the GDP-to-GTP exchange. Importantly, GAP not only donates its R789 arginine finger, but stabilizes the catalytically-competent conformation and pre-organizes catalytic residue Q61; mutations disturb the R789/Q61 organization, impairing GAP-mediated GTP hydrolysis. Together, our simulations help provide a mechanistic explanation of key mutational events in one of the most oncogenic proteins in cancer.

Effect of single point mutations in a form of systemic amyloidosis

Amyloid deposits of light-chain proteins are associated with the most common form of systemic amyloidosis. We have studied the effects of single point mutations on amyloid formation of these proteins using explicit solvent model molecular dynamics simulations. For this purpose, we compare the stability of the wild-type immunoglobulin light-chain protein REI in its native and amyloid forms with that of four mutants: R61N, G68D, D82I, and A84T. We argue that the experimentally observed differences in the propensity for amyloid formation result from two effects. First, the mutant dimers have a lower stability than the wild-type dimer due to increase exposure of certain hydrophobic residues. The second effect is a shift in equilibrium between monomers with amyloid-like structure and such with native structures. Hence, when developing drugs against light-chain associated systemic amyloidosis, one should look for components that either stabilize the dimer by binding to the dimer interface or reduce for the monomers the probability of the amyloid form.

On the lack of polymorphism in Aβ-peptide aggregates derived from patient brains

The amyloid beta (Aβ) oligomers and fibrils that are found in neural tissues of patients suffering from Alzheimer's disease may either cause or contribute to the pathology of the disease. In vitro, these Aβ-aggregates are characterized by structural polymorphism. However, recent solid state NMR data of fibrils acquired post mortem from the brains of two Alzheimer's patients indicate presence of only a single, patient-specific structure. Using enhanced molecular dynamic simulations we investigate the factors that modulate the stability of Aβ-fibrils. We find characteristic differences in molecular flexibility, dynamics of interactions, and structural behavior between the brain-derived Aβ-fibril structure and in vitro models. These differences may help to explain the lack of polymorphism in fibrils collected from patient brains, and have to be taken into account when designing aggregation inhibitors and imaging agents for Alzheimer's disease.