Efficient and accurate binding free energy calculation of Aß9-40 protofilament propagation.

Self-assembled aggregation of peptides and proteins into regular amyloid fibrils is associated with several neurodegenerative diseases. In case of Alzheimer's disease proteolytic cleavage products of the amyloid precursor protein form pathological amyloid-beta fibrils in a nucleation and propagation phase. The molecular details and thermodynamic driving forces of amyloid formation are not well understood, but are of high relevance for potential pharmacological interference. We used atomistic binding free energy simulations to calculate the free energy of protofilament propagation by an additional Aß9-40 peptide binding to the protofilament tip. It requires sampling of relevant conformational transitions which is challenging since the monomeric Aß9-40 peptide is intrinsically disordered. However, the convergence of umbrella simulations can be enhanced by applying additional restraining potentials on the axial, orientational and conformational degrees of freedom. The improved convergence leads to a much closer agreement with experimental binding free energy data compared to unrestrained umbrella sampling. Moreover, the restraining approach results in a separation of contributions to the total binding free energy. The calculated contributions indicate that the free energy change associated with the restriction of conformational freedom upon propagation makes a large opposing contribution of higher magnitude than the total binding free energy. Finally, optimization of the approach leads to further significant reduction of the computational demand which is crucial for systematic studies on mutations, denaturants and inhibitors in the fibril propagation step.

Targeting Telomeres: Molecular Dynamics and Free Energy Simulation of Gold-Carbene Binding to DNA.

DNA sequences in regulatory regions and in telomers at the ends of chromosomes frequently contain tandem repeats of guanine nucleotides that can form stacked structures stabilized by Hoogsten pairing and centrally bound monovalent cations. The replication and elongation of telomeres requires the disruption of these G-quadruplex structures. Hence, drug molecules such as gold (Au)-carbene that stabilize G-quadruplexes may also interfere with the elongation of telomeres and, in turn, could be used to control cell replication and growth. To better understand the molecular mechanism of Au-carbene binding to G-quadruplexes, we employed molecular dynamics simulations and free energy simulations. Whereas very restricted mobility of two Au-carbene ligands was found upon binding as a doublet to one side of the G-quadruplex, much larger translational and orientational mobility was observed for a single Au-carbene binding at the second G-quadruplex surface. Comparative simulations on duplex DNA in the presence of Au-carbene ligands indicates a preference for the minor groove and weaker unspecific and more salt-dependent binding than to the G-quadruplex surface. Analysis of energetic contributions reveals a dominance of nonpolar and van der Waals interactions to drive binding. The simulations can also be helpful for proposing possible modifications that could improve Au-carbene affinity and specificity for G-quadruplex binding.

From monomer to fibril: Abeta-amyloid binding to Aducanumab antibody studied by molecular dynamics simulation.

Alzheimer's disease is one of the most common causes of dementia. It is believed that the aggregation of short Aß-peptides to form oligomeric and protofibrillar amyloid assemblies plays a central role for disease-relevant neurotoxicity. In recent years, passive immunotherapy has been introduced as a potential treatment strategy with anti-amyloid antibodies binding to Aß-amyloids and inducing their subsequent degradation by the immune system. Although so far mostly unsuccessful in clinical studies, the high-dosed application of the monoclonal antibody Aducanumab has shown therapeutic potential that might be attributed to its much greater affinity to Aß-aggregates vs monomeric Aß-peptides. In order to better understand how Aducanumab interacts with aggregated Aß-forms compared to monomers, we have generated structural model complexes based on the known structure of Aducanumab in complex with an Aß2 - 7 -eptitope. Structural models of Aducanumab bound to full-sequence Aß1 - 40 -monomers, oligomers, protofilaments and mature fibrils were generated and investigated using extensive molecular dynamics simulations to characterize the flexibility and possible additional interactions. Indeed, an aggregate-specific N-terminal binding motif was found in case of Aducanumab binding to oligomers, protofilaments and fibrils that is located next to but not overlapping with the epitope binding site found in the crystal structure with Aß2 - 7 . Analysis of binding energetics indicates that this motif binds weaker than the epitope but likely contributes to Aducanumab's preference for aggregated Aß-species. The predicted aggregate-specific binding motif could potentially serve as a basis to reengineer Aducanumab for further enhanced preference to bind Aß-aggregates vs monomers.

Folding and Unfolding of the Short Light-Triggered ß-Hairpin Peptide AzoChignolin Occurs within 100 ns.

To map the underlying molecular mechanisms of folding dynamics in proteins, light-operated peptides have emerged as promising tools. In this study, we reveal the complete sequence of light-induced structural changes of AzoChignolin, a short ß-hairpin peptide containing an azobenzene photoswitch in its loop region. Light-triggered structural changes were monitored by time-resolved IR spectroscopy. Formation and destruction of the hairpin structure is very fast and occurs within 100 ns for AzoChignolin in methanol. Atomistic molecular dynamics simulations using two explicit solvents, methanol and water, revealed the underlying molecular processes and allowed us to gain further insight into the reaction mechanism. Despite its rapid reaction time, hairpin formation in these solvents is not force-driven by the molecular switch but proceeded via formation of interstrand hydrogen bonds and contacts between aromatic residues. Moreover, the combined experimental and theoretical study demonstrates that the solvent (methanol vs water) does not dictate the velocity of ß-hairpin formation in the AzoChignolin peptide comprising only a few hydrophobic residues in the strands.

Designed Macrocyclic Peptides as Nanomolar Amyloid Inhibitors Based on Minimal Recognition Elements.

Amyloid self-assembly is linked to the pathogenesis of Alzheimer's disease (AD) and type 2 diabetes (T2D), but so far, no anti-amyloid compound has reached the clinic. Macrocyclic peptides belong to the most attractive drug candidates. Herein we present macrocyclic peptides (MCIPs) designed using minimal IAPP-derived recognition elements as a novel class of nanomolar amyloid inhibitors of both Aß40(42) and IAPP or Aß40(42) alone and show that chirality controls inhibitor selectivity. Sequence optimization led to the discovery of an Aß40(42)-selective MCIP exhibiting high proteolytic stability in human plasma and human blood-brain barrier (BBB) crossing ability in a cell model, two highly desirable properties for anti-amyloid AD drugs. Owing to their favorable properties, MCIPs should serve as leads for macrocyclic peptide-based anti-amyloid drugs and scaffolds for the design of small-molecule peptidomimetics for targeting amyloidogenesis in AD or in both AD and T2D.

Key aromatic/hydrophobic amino acids controlling a cross-amyloid peptide interaction versus amyloid self-assembly.

The interaction of the intrinsically disordered polypeptide islet amyloid polypeptide (IAPP), which is associated with type 2 diabetes (T2D), with the Alzheimer's disease amyloid-ß (Aß) peptide modulates their self-assembly into amyloid fibrils and may link the pathogeneses of these two cell-degenerative diseases. However, the molecular determinants of this interaction remain elusive. Using a systematic alanine scan approach, fluorescence spectroscopy, and other biophysical methods, including heterocomplex pulldown assays, far-UV CD spectroscopy, the thioflavin T binding assay, transmission EM, and molecular dynamics simulations, here we identified single aromatic/hydrophobic residues within the amyloid core IAPP region as hot spots or key residues of its cross-interaction with Aß40(42) peptide. Importantly, we also find that none of these residues in isolation plays a key role in IAPP self-assembly, whereas simultaneous substitution of four aromatic/hydrophobic residues with Ala dramatically impairs both IAPP self-assembly and hetero-assembly with Aß40(42). Furthermore, our experiments yielded several novel IAPP analogs, whose sequences are highly similar to that of IAPP but have distinct amyloid self- or cross-interaction potentials. The identified similarities and major differences controlling IAPP cross-peptide interaction with Aß40(42) versus its amyloid self-assembly offer a molecular basis for understanding the underlying mechanisms. We propose that these insights will aid in designing intervention strategies and novel IAPP analogs for the management of type 2 diabetes, Alzheimer's disease, or other diseases related to IAPP dysfunction or cross-amyloid interactions.

Optical control of a receptor-linked guanylyl cyclase using a photoswitchable peptidic hormone.

The optical control over biological function with small photoswitchable molecules has gathered significant attention in the last decade. Herein, we describe the design and synthesis of a small library of photoswitchable peptidomimetics based upon human atrial natriuretic peptide (ANP), in which the photochromic amino acid [3-(3-aminomethyl)phenylazo]phenylacetic acid (AMPP) is incorporated into the peptide backbone. The endogeneous hormone ANP signals via the natriuretic peptide receptor A (NPR-A) through raising intracellular cGMP concentrations, and is involved in blood pressure regulation and sodium homeostasis, as well as lipid metabolism and pancreatic function. The cis- and trans-isomers of one of our peptidomimetics, termed TOP271, exhibit a four-fold difference in NPR-A mediated cGMP synthesis in vitro. Despite this seemingly small difference, TOP271 enables large, optically-induced conformational changes ex vivo and transforms the NPR-A into an endogenous photoswitch. Thus, application of TOP271 allows the reversible generation of cGMP using light and remote control can be afforded over vasoactivity in explanted murine aortic rings, as well as pancreatic beta cell function in islets of Langerhans. This study demonstrates the broad applicability of TOP271 to enzyme-dependent signalling processes, extends the toolbox of photoswitchable molecules to all classes of transmembrane receptors and utilizes photopharmacology to deduce receptor activation on a molecular level.

From Aß Filament to Fibril: Molecular Mechanism of Surface-Activated Secondary Nucleation from All-Atom MD Simulations.

Secondary nucleation pathways in which existing amyloid fibrils catalyze the formation of new aggregates and neurotoxic oligomers are of immediate importance for the onset and progression of Alzheimer's disease. Here, we apply extensive all-atom molecular dynamics simulations in explicit water to study surface-activated secondary nucleation pathways at the extended lateral ß-sheet surface of a preformed Aß9-40 filament. Calculation of free-energy profiles allows us to determine binding free energies and conformational intermediates for nucleation complexes consisting of 1-4 Aß peptides. In addition, we combine the free-energy profiles with position-dependent diffusion profiles to extract complementary kinetic information and macroscopic growth rates. Single monomers bind to the ß-sheet surface in a disordered, hydrophobically collapsed conformation, whereas dimers and larger oligomers can retain a cross-ß conformation resembling a more ordered fibril structure. The association processes during secondary nucleation follow a dock/lock mechanism consisting of a fast initial encounter phase (docking) and a slow structural rearrangement phase (locking). The major driving forces for surface-activated secondary nucleation are the release of a large number of hydration water molecules and the formation of hydrophobic interface contacts, the latter being in contrast to the elongation process at filament tips, which is dominated by the formation of stable and highly specific interface hydrogen bonds. The calculated binding free energies and the association rates for the attachment of Aß monomers and oligomers to the extended lateral ß-sheet surface of the filament seed are higher compared to those for elongation at the filament tips, indicating that secondary nucleation pathways can become important once a critical concentration of filaments has formed.

Dynamics of Seeded Aß40-Fibril Growth from Atomistic Molecular Dynamics Simulations: Kinetic Trapping and Reduced Water Mobility in the Locking Step.

Filamentous ß-amyloid aggregates are crucial for the pathology of Alzheimer's disease. Despite the tremendous biomedical importance, the molecular pathway of growth propagation is not completely understood and remains challenging to investigate by simulations due to the long time scales involved. Here, we apply extensive all-atom molecular dynamics simulations in explicit water to obtain free energy profiles and kinetic information from position-dependent diffusion profiles for three different Aß9-40-growth processes: fibril elongation by single monomers at the structurally unequal filament tips and association of larger filament fragments. Our approach provides insight into the molecular steps of the kinetic pathway and allows close agreement with experimental binding free energies and macroscopic growth rates. Water plays a decisive role, and solvent entropy is identified as the main driving force for assembly. Fibril growth is disfavored energetically due to cancellation of direct peptide-peptide interactions and solvation effects. The kinetics of growth is consistent with the characteristic dock/lock mechanism, and docking is at least 2 orders of magnitude faster. During initial docking, interactions are mediated by transient non-native hydrogen bonds, which efficiently catch the incoming monomer or fragment already at separations of about 3 nm. In subsequent locking, the dynamics is much slower due to formation of kinetically trapped conformations caused by long-lived non-native hydrogen bonds. Fibril growth additionally requires collective motion of water molecules to create a dry binding interface. Fibril growth is further retarded due to reduced mobility of the involved hydration water, evident from a 2-fold reduction of the diffusion coefficient.