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.

Insights into the Structural Basis of Amyloid Resistance Provided by Cryo-EM Structures of AApoAII Amyloid Fibrils.

Amyloid resistance is the inability or the reduced susceptibility of an organism to develop amyloidosis. In this study we have analysed the molecular basis of the resistance to systemic AApoAII amyloidosis, which arises from the formation of amyloid fibrils from apolipoprotein A-II (ApoA-II). The disease affects humans and animals, including SAMR1C mice that express the C allele of ApoA-II protein, whereas other mouse strains are resistant to development of amyloidosis due to the expression of other ApoA-II alleles, such as ApoA-IIF. Using cryo-electron microscopy, molecular dynamics simulations and other methods, we have determined the structures of pathogenic AApoAII amyloid fibrils from SAMR1C mice and analysed the structural effects of ApoA-IIF-specific mutational changes. Our data show that these changes render ApoA-IIF incompatible with the specific fibril morphologies, with which ApoA-II protein can become pathogenic in vivo.

RNA G-quadruplex folding is a multi-pathway process driven by conformational entropy.

The kinetics of folding is crucial for the function of many regulatory RNAs including RNA G-quadruplexes (rG4s). Here, we characterize the folding pathways of a G-quadruplex from the telomeric repeat-containing RNA by combining all-atom molecular dynamics and coarse-grained simulations with circular dichroism experiments. The quadruplex fold is stabilized by cations and thus, the ion atmosphere forming a double layer surrounding the highly charged quadruplex guides the folding process. To capture the ionic double layer in implicit solvent coarse-grained simulations correctly, we develop a matching procedure based on all-atom simulations in explicit water. The procedure yields quantitative agreement between simulations and experiments as judged by the populations of folded and unfolded states at different salt concentrations and temperatures. Subsequently, we show that coarse-grained simulations with a resolution of three interaction sites per nucleotide are well suited to resolve the folding pathways and their intermediate states. The results reveal that the folding progresses from unpaired chain via hairpin, triplex and double-hairpin constellations to the final folded structure. The two- and three-strand intermediates are stabilized by transient Hoogsteen interactions. Each pathway passes through two on-pathway intermediates. We hypothesize that conformational entropy is a hallmark of rG4 folding. Conformational entropy leads to the observed branched multi-pathway folding process for TERRA25. We corroborate this hypothesis by presenting the free energy landscapes and folding pathways of four rG4 systems with varying loop length.

Cryo-EM Analysis of the Effect of Seeding with Brain-derived Aß Amyloid Fibrils.

Aß amyloid fibrils from Alzheimer's brain tissue are polymorphic and structurally different from typical in vitro formed Aß fibrils. Here, we show that brain-derived (ex vivo) fibril structures can be proliferated by seeding in vitro. The proliferation reaction is only efficient for one of the three abundant ex vivo Aß fibril morphologies, which consists of two peptide stacks, while the inefficiently proliferated fibril morphologies contain four or six peptide stacks. In addition to the seeded fibril structures, we find that de novo nucleated fibril structures can emerge in seeded samples if the seeding reaction is continued over multiple generations. These data imply a competition between de novo nucleation and seed extension and suggest further that seeding favours the outgrowth of fibril morphologies that contain fewer peptide stacks.

pH-dependent structural transitions in cationic ionizable lipid mesophases are critical for lipid nanoparticle function.

Lipid nanoparticles (LNPs) are advanced core-shell particles for messenger RNA (mRNA) based therapies that are made of polyethylene glycol (PEG) lipid, distearoylphosphatidylcholine (DSPC), cationic ionizable lipid (CIL), cholesterol (chol), and mRNA. Yet the mechanism of pH-dependent response that is believed to cause endosomal release of LNPs is not well understood. Here, we show that eGFP (enhanced green fluorescent protein) protein expression in the mouse liver mediated by the ionizable lipids DLin-MC3-DMA (MC3), DLin-KC2-DMA (KC2), and DLinDMA (DD) ranks MC3 ≥ KC2 > DD despite similar delivery of mRNA per cell in all cell fractions isolated. We hypothesize that the three CIL-LNPs react differently to pH changes and hence study the structure of CIL/chol bulk phases in water. Using synchrotron X-ray scattering a sequence of ordered CIL/chol mesophases with lowering pH values are observed. These phases show isotropic inverse micellar, cubic Fd3m inverse micellar, inverse hexagonal [Formula: see text] and bicontinuous cubic Pn3m symmetry. If polyadenylic acid, as mRNA surrogate, is added to CIL/chol, excess lipid coexists with a condensed nucleic acid lipid [Formula: see text] phase. The next-neighbor distance in the excess phase shows a discontinuity at the Fd3m inverse micellar to inverse hexagonal [Formula: see text] transition occurring at pH 6 with distinctly larger spacing and hydration for DD vs. MC3 and KC2. In mRNA LNPs, DD showed larger internal spacing, as well as retarded onset and reduced level of DD-LNP-mediated eGFP expression in vitro compared to MC3 and KC2. Our data suggest that the pH-driven Fd3m-[Formula: see text] transition in bulk phases is a hallmark of CIL-specific differences in mRNA LNP efficacy.

Adsorbing DNA to Mica by Cations: Influence of Valency and Ion Type.

Ion-mediated attraction between DNA and mica plays a crucial role in biotechnological applications and molecular imaging. Here, we combine molecular dynamics simulations and single-molecule atomic force microscopy experiments to characterize the detachment forces of single-stranded DNA at mica surfaces mediated by the metal cations Li+, Na+, K+, Cs+, Mg2+, and Ca2+. Ion-specific adsorption at the mica/water interface compensates (Li+ and Na+) or overcompensates (K+, Cs+, Mg2+, and Ca2+) the bare negative surface charge of mica. In addition, direct and water-mediated contacts are formed between the ions, the phosphate oxygens of DNA, and mica. The different contact types give rise to low- and high-force pathways and a broad distribution of detachment forces. Weakly hydrated ions, such as Cs+ and water-mediated contacts, lead to low detachment forces and high mobility of the DNA on the surface. Direct ion-DNA or ion-surface contacts lead to significantly higher forces. The comprehensive view gained from our combined approach allows us to highlight the most promising cations for imaging in physiological conditions: K+, which overcompensates the negative mica charge and induces long-ranged attractions. Mg2+ and Ca2+, which form a few specific and long-lived contacts to bind DNA with high affinity.

Combining molecular dynamics simulations and x-ray scattering techniques for the accurate treatment of protonation degree and packing of ionizable lipids in monolayers.

The pH-dependent change in protonation of ionizable lipids is crucial for the success of lipid-based nanoparticles as mRNA delivery systems. Despite their widespread application in vaccines, the structural changes upon acidification are not well understood. Molecular dynamics simulations support structure prediction but require an a priori knowledge of the lipid packing and protonation degree. The presetting of the protonation degree is a challenging task in the case of ionizable lipids since it depends on pH and on the local lipid environment and often lacks experimental validation. Here, we introduce a methodology of combining all-atom molecular dynamics simulations with experimental total-reflection x-ray fluorescence and scattering measurements for the ionizable lipid Dlin-MC3-DMA (MC3) in POPC monolayers. This joint approach allows us to simultaneously determine the lipid packing and the protonation degree of MC3. The consistent parameterization is expected to be useful for further predictive modeling of the action of MC3-based lipid nanoparticles.

Cryo-EM Structure of the Full-length hnRNPA1 Amyloid Fibril.

Heterogeneous nuclear ribonucleoprotein A1 (hnRNPA1) is a multifunctional RNA-binding protein that is associated with neurodegenerative diseases, such as amyotrophic lateral sclerosis and multisystem proteinopathy. In this study, we have used cryo-electron microscopy to investigate the three-dimensional structure of amyloid fibrils from full-length hnRNPA1 protein. We find that the fibril core is formed by a 45-residue segment of the prion-like low-complexity domain of the protein, whereas the remaining parts of the protein (275 residues) form a fuzzy coat around the fibril core. The fibril consists of two fibril protein stacks that are arranged into a pseudo-21 screw symmetry. The ordered core harbors several of the positions that are known to be affected by disease-associated mutations, but does not encompass the most aggregation-prone segments of the protein. These data indicate that the structures of amyloid fibrils from full-length proteins may be more complex than anticipated by current theories on protein misfolding.

Structural insights on ionizable Dlin-MC3-DMA lipids in DOPC layers by combining accurate atomistic force fields, molecular dynamics simulations and neutron reflectivity.

Ionizable lipids such as the promising Dlin-MC3-DMA (MC3) are essential for the successful design of lipid nanoparticles (LNPs) as drug delivery agents. Combining molecular dynamics simulations with experimental data, such as neutron reflectivity experiments and other scattering techniques, is essential to provide insights into the internal structure of LNPs, which is not fully understood to date. However, the accuracy of the simulations relies on the choice of force field parameters and high-quality experimental data is indispensable to verify the parametrization. For MC3, different parameterizations in combination with the CHARMM and the Slipids force fields have recently emerged. Here, we complement the existing efforts by providing parameters for cationic and neutral MC3 compatible with the AMBER Lipid17 force field. Subsequently, we carefully assess the accuracy of the different force fields by providing a direct comparison to neutron reflectivity experiments of mixed lipid bilayers consisting of MC3 and DOPC at different pHs. At low pH (cationic MC3) and at high pH (neutral MC3) the newly developed MC3 parameters in combination with AMBER Lipid17 for DOPC give good agreement with the experiments. Overall, the agreement is similar compared to the Park-Im parameters for MC3 in combination with the CHARMM36 force field for DOPC. The Ermilova-Swenson MC3 parameters in combination with the Slipids force field underestimate the bilayer thickness. While the distribution of cationic MC3 is very similar, the different force fields for neutral MC3 reveal distinct differences ranging from strong accumulation in the membrane center (current MC3/AMBER Lipid17 DOPC), over mild accumulation (Park-Im MC3/CHARMM36 DOPC) to surface accumulation (Ermilova-Swenson MC3/Slipids DOPC). These pronounced differences highlight the importance of accurate force field parameters and their experimental validation.

Remembering the Work of Phillip L. Geissler: A Coda to His Scientific Trajectory.

Phillip L. Geissler made important contributions to the statistical mechanics of biological polymers, heterogeneous materials, and chemical dynamics in aqueous environments. He devised analytical and computational methods that revealed the underlying organization of complex systems at the frontiers of biology, chemistry, and materials science. In this retrospective we celebrate his work at these frontiers.

RNA Captures More Cations than DNA: Insights from Molecular Dynamics Simulations.

The distribution of cations around nucleic acids is essential for a broad variety of processes ranging from DNA condensation and RNA folding to the detection of biomolecules in biosensors. Predicting the exact distribution of ions remains challenging since the distribution and, hence, a broad variety of nucleic acid properties depend on the salt concentration, the valency of the ions, and the ion type. Despite the importance, a general theory to quantify ion-specific effects for highly charged biomolecules is still lacking. Moreover, recent experiments reveal that despite their similar building blocks, DNA and RNA duplexes can react differently to the same ionic conditions. The aim of our current work is to provide a comprehensive set of molecular dynamics simulations using more than 180 µs of simulation time. For the mono- and divalent cations Li+, Na+, K+, Cs+, Ca2+, Sr2+, and Ba2+, the simulations allow us to reveal the ion-specific distributions and binding patterns for DNA and RNA duplexes. The microscopic insights from the simulations display the origin of ion-specificity and shed light on the question of why DNA and RNA show opposing behavior in the same ionic conditions. Finally, the detailed binding patterns from the simulations reveal why RNA can capture more cations than DNA.

SAA fibrils involved in AA amyloidosis are similar in bulk and by single particle reconstitution: A MAS solid-state NMR study.

AA amyloidosis is one of the most prevalent forms of systemic amyloidosis and affects both humans and other vertebrates. In this study, we compare MAS solid-state NMR data with a recent cryo-EM study of fibrils involving full-length murine SAA1.1. We address the question whether the specific requirements for the reconstitution of an amyloid fibril structure by cryo-EM can potentially yield a bias towards a particular fibril polymorph. We employ fibril seeds extracted from in to vivo material to imprint the fibril structure onto the biochemically produced protein. Sequential assignments yield the secondary structure elements in the fibril state. Long-range DARR and PAR experiments confirm largely the topology observed in the ex-vivo cryo-EM study. We find that the ß-sheets identified in the NMR experiments are similar to the ß-sheets found in the cryo-EM study, with the exception of amino acids 33-42. These residues cannot be assigned by solid-state NMR, while they adopt a stable ß-sheet in the cryo-EM structure. We suggest that the differences between MAS solid-state NMR and cryo-EM data are a consequence of a second conformer involving residues 33-42. Moreover, we were able to characterize the dynamic C-terminal tail of SAA in the fibril state. The C-terminus is flexible, remains detached from the fibrils, and does not affect the SAA fibril structure as confirmed further by molecular dynamics simulations. As the C-terminus can potentially interact with other cellular components, binding to cellular targets can affect its accessibility for protease digestion.

Twisting DNA by salt.

The structure and properties of DNA depend on the environment, in particular the ion atmosphere. Here, we investigate how DNA twist -one of the central properties of DNA- changes with concentration and identity of the surrounding ions. To resolve how cations influence the twist, we combine single-molecule magnetic tweezer experiments and extensive all-atom molecular dynamics simulations. Two interconnected trends are observed for monovalent alkali and divalent alkaline earth cations. First, DNA twist increases monotonously with increasing concentration for all ions investigated. Second, for a given salt concentration, DNA twist strongly depends on cation identity. At 100 mM concentration, DNA twist increases as Na+ < K+ < Rb+ < Ba2+ < Li+ ≈ Cs+ < Sr2+ < Mg2+ < Ca2+. Our molecular dynamics simulations reveal that preferential binding of the cations to the DNA backbone or the nucleobases has opposing effects on DNA twist and provides the microscopic explanation of the observed ion specificity. However, the simulations also reveal shortcomings of existing force field parameters for Cs+ and Sr2+. The comprehensive view gained from our combined approach provides a foundation for understanding and predicting cation-induced structural changes both in nature and in DNA nanotechnology.

Combining Coarse-Grained Simulations and Single Molecule Analysis Reveals a Three-State Folding Model of the Guanidine-II Riboswitch.

Riboswitch RNAs regulate gene expression by conformational changes induced by environmental conditions and specific ligand binding. The guanidine-II riboswitch is proposed to bind the small molecule guanidinium and to subsequently form a kissing loop interaction between the P1 and P2 hairpins. While an interaction was shown for isolated hairpins in crystallization and electron paramagnetic resonance experiments, an intrastrand kissing loop formation has not been demonstrated. Here, we report the first evidence of this interaction in cis in a ligand and Mg2+ dependent manner. Using single-molecule FRET spectroscopy and detailed structural information from coarse-grained simulations, we observe and characterize three interconvertible states representing an open and kissing loop conformation as well as a novel Mg2+ dependent state for the guanidine-II riboswitch from E. coli. The results further substantiate the proposed switching mechanism and provide detailed insight into the regulation mechanism for the guanidine-II riboswitch class. Combining single molecule experiments and coarse-grained simulations therefore provides a promising perspective in resolving the conformational changes induced by environmental conditions and to yield molecular insights into RNA regulation.

Magnesium force fields for OPC water with accurate solvation, ion-binding, and water-exchange properties: Successful transfer from SPC/E.

Magnesium plays a vital role in a large variety of biological processes. To model such processes by molecular dynamics simulations, researchers rely on accurate force field parameters for Mg2+ and water. OPC is one of the most promising water models yielding an improved description of biomolecules in water. The aim of this work is to provide force field parameters for Mg2+ that lead to accurate simulation results in combination with OPC water. Using 12 different Mg2+ parameter sets that were previously optimized with different water models, we systematically assess the transferability to OPC based on a large variety of experimental properties. The results show that the Mg2+ parameters for SPC/E are transferable to OPC and closely reproduce the experimental solvation free energy, radius of the first hydration shell, coordination number, activity derivative, and binding affinity toward the phosphate oxygens on RNA. Two optimal parameter sets are presented: MicroMg yields water exchange in OPC on the microsecond timescale in agreement with experiments. NanoMg yields accelerated exchange on the nanosecond timescale and facilitates the direct observation of ion binding events for enhanced sampling purposes.

Cryo-EM demonstrates the in vitro proliferation of an ex vivo amyloid fibril morphology by seeding.

Several studies showed that seeding of solutions of monomeric fibril proteins with ex vivo amyloid fibrils accelerated the kinetics of fibril formation in vitro but did not necessarily replicate the seed structure. In this research we use cryo-electron microscopy and other methods to analyze the ability of serum amyloid A (SAA)1.1-derived amyloid fibrils, purified from systemic AA amyloidosis tissue, to seed solutions of recombinant SAA1.1 protein. We show that 98% of the seeded fibrils remodel the full fibril structure of the main ex vivo fibril morphology, which we used for seeding, while they are notably different from unseeded in vitro fibrils. The seeded fibrils show a similar proteinase K resistance as ex vivo fibrils and are substantially more stable to proteolytic digestion than unseeded in vitro fibrils. Our data support the view that the fibril morphology contributes to determining proteolytic stability and that pathogenic amyloid fibrils arise from proteolytic selection.

Artificial Intelligence Resolves Kinetic Pathways of Magnesium Binding to RNA.

Magnesium is an indispensable cofactor in countless vital processes. In order to understand its functional role, the characterization of the binding pathways to biomolecules such as RNA is crucial. Despite the importance, a molecular description is still lacking since the transition from the water-mediated outer-sphere to the direct inner-sphere coordination is on the millisecond time scale and therefore out of reach for conventional simulation techniques. To fill this gap, we use transition path sampling to resolve the binding pathways and to elucidate the role of the solvent in the binding process. The results reveal that the molecular void provoked by the leaving phosphate oxygen of the RNA is immediately filled by an entering water molecule. In addition, water molecules from the first and second hydration shell couple to the concerted exchange. To capture the intimate solute-solvent coupling, we perform a committor analysis as the basis for a machine learning algorithm that derives the optimal deep learning model from thousands of scanned architectures using hyperparameter tuning. The results reveal that the properly optimized deep network architecture recognizes the important solvent structures, extracts the relevant information, and predicts the commitment probability with high accuracy. Our results provide detailed insights into the solute-solvent coupling which is ubiquitous for kosmotropic ions and governs a large variety of biochemical reactions in aqueous solutions.

Optimized Magnesium Force Field Parameters for Biomolecular Simulations with Accurate Solvation, Ion-Binding, and Water-Exchange Properties in SPC/E, TIP3P-fb, TIP4P/2005, TIP4P-Ew, and TIP4P-D.

Magnesium is essential in many vital processes. To correctly describe Mg2+ in physiological processes by molecular dynamics simulations, accurate force fields are fundamental. Despite the importance, force fields based on the commonly used 12-6 Lennard-Jones potential showed significant shortcomings. Recently progress was made by an optimization procedure that implicitly accounts for polarizability. The resulting microMg and nanoMg force fields (J. Chem. Theory Comput. 2021, 17, 2530-2540) accurately reproduce a broad range of experimental solution properties and the binding affinity to nucleic acids in TIP3P water. Since countless simulation studies rely on available water models and ion force fields, we here extend the optimization and provide Mg2+ parameters in combination with the SPC/E, TIP3P-fb, TIP4P/2005, TIP4P-Ew, and TIP4P-D water models. For each water model, the Mg2+ force fields reproduce the solvation free energy, the distance to oxygens in the first hydration shell, the hydration number, the activity coefficient derivative in MgCl2 solutions, and the binding affinity and distance to the phosphate oxygens on nucleic acids. We present two parameter sets: MicroMg yields water exchange on the microsecond time scale and matches the experimental exchange rate. Depending on the water model, nanoMg yields accelerated water exchange in the range of 106 to 108 exchanges per second. The nanoMg parameters can be used to enhance the sampling of binding events, to obtain converged distributions of Mg2+, or to predict ion binding sites in biomolecular simulations. The parameter files are freely available at https://github.com/bio-phys/optimizedMgFFs.

Kinetic pathways of water exchange in the first hydration shell of magnesium: Influence of water model and ionic force field.

Water exchange between the first and second hydration shell is essential for the role of Mg2+ in biochemical processes. In order to provide microscopic insights into the exchange mechanism, we resolve the exchange pathways by all-atom molecular dynamics simulations and transition path sampling. Since the exchange kinetics relies on the choice of the water model and the ionic force field, we systematically investigate the influence of seven different polarizable and non-polarizable water and three different Mg2+ models. In all cases, water exchange can occur either via an indirect or direct mechanism (exchanging molecules occupy different/same position on the water octahedron). In addition, the results reveal a crossover from an interchange dissociative (Id) to an associative (Ia) reaction mechanism dependent on the range of the Mg2+-water interaction potential of the respective force field. Standard non-polarizable force fields follow the Id mechanism in agreement with experimental results. By contrast, polarizable and long-ranged non-polarizable force fields follow the Ia mechanism. Our results provide a comprehensive view on the influence of the water model and the ionic force field on the exchange dynamics and the foundation to assess the choice of the force field in biomolecular simulations.

Extended magnesium and calcium force field parameters for accurate ion-nucleic acid interactions in biomolecular simulations.

Magnesium and calcium play an essential role in the folding and function of nucleic acids. To correctly describe their interactions with DNA and RNA in biomolecular simulations, an accurate parameterization is crucial. In most cases, the ion parameters are optimized based on a set of experimental solution properties such as solvation free energies, radial distribution functions, water exchange rates, and activity coefficient derivatives. However, the transferability of such bulk-optimized ion parameters to quantitatively describe biomolecular systems is limited. Here, we extend the applicability of our previous bulk-optimized parameters by including experimental binding affinities toward the phosphate oxygen on nucleic acids. In particular, we systematically adjust the combination rules that are an integral part of the pairwise interaction potentials of classical force fields. This allows us to quantitatively describe specific ion binding to nucleic acids without changing the solution properties in the most simple and efficient way. We show the advancement of the optimized Lorentz combination rule for two representative nucleic acid systems. For double-stranded DNA, the optimized combination rule for Ca2+ significantly improves the agreement with experiments, while the standard combination rule leads to unrealistically distorted DNA structures. For the add A-riboswitch, the optimized combination rule for Mg2+ improves the structure of two specifically bound Mg2+ ions as judged by the experimental distance to the binding site. Including experimental binding affinities toward specific ion binding sites on biomolecules, therefore, provides a promising perspective to develop a more accurate description of metal cations for biomolecular simulations.