Search NASA⌕ Search

SEARCH · Search NASA

Results for “Docking”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Design of Diverse Asymmetric Pockets in De Novo Homo-oligomeric Proteins

A challenge for design of protein–small-molecule recognition is that incorporation of cavities with size, shape, and composition suitable for specific recognition can considerably destabilize protein monomers. This challenge can be overcome through binding pockets formed at homo-oligomeric interfaces between folded monomers. Interfaces surrounding the central homo-oligomer symmetry axes necessarily have the same symmetry and so may not be well suited to binding asymmetric molecules. To enable general recognition of arbitrary asymmetric substrates and small molecules, we developed an approach to designing asymmetric interfaces at off-axis sites on homo-oligomers, analogous to those found in native homo-oligomeric proteins such as glutamine synthetase. We symmetrically dock curved helical repeat proteins such that they form pockets at the asymmetric interface of the oligomer with sizes ranging from several angstroms, appropriate for binding a single ion, to up to more than 20 Å across. Of the 133 proteins tested, 84 had soluble expression in E. coli, 47 had correct oligomeric states in solution, 35 had small-angle X-ray scattering (SAXS) data largely consistent with design models, and 8 had negative-stain electron microscopy (nsEM) 2D class averages showing the structures coming together as designed. Both an X-ray crystal structure and a cryogenic electron microscopy (cryoEM) structure are close to the computational design models. The nature of these proteins as homo-oligomers allows them to be readily built into higher-order structures such as nanocages, and the asymmetric pockets of these structures open rich possibilities for small-molecule binder design free from the constraints associated with monomer destabilization.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

The Conformational Change in Elongation Factor Tu Involves Separation of Its Domains

Elongation factor Tu (EF-Tu) is a highly conserved GTPase that is responsible for supplying the aminoacylated tRNA to the ribosome. Upon binding to the ribosome, EF-Tu undergoes GTP hydrolysis, which drives a major conformational change, triggering the release of aminoacylated tRNA to the ribosome. Using a combination of molecular simulation techniques, we studied the transition between the pre- and post-hydrolysis structures through two distinct pathways. Here, we show that the transition free energy is minimal along a non-intuitive pathway that involves “separation” of the GTP binding domain (domain 1) from the OB folds (domains 2 and 3), followed by domain 1 rotation, and, eventually, locking the EF-Tu conformation in the post-hydrolysis state. The domain separation also leads to a slight extension of the linker connecting domain 1 to domain 2. Using docking tools and correlation-based analysis, we identified and characterized the EF-Tu conformations that release the tRNA. These calculations suggest that EF-Tu can release the tRNA before the domains separate and after domain 1 rotates by 25°. Lastly, we also examined the EF-Tu conformations in the context of the ribosome. Given the high degrees of sequence similarity with other translational GTPases, we predict a similar separation mechanism is followed.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

Biochemical Characterization and Structure-Based Mutational Analysis Provide Insight into the Binding and Mechanism of Action of Novel Aspartate Aminotransferase Inhibitors

Pancreatic cancer cells are characterized by deregulated metabolic programs that facilitate growth and resistance to oxidative stress. Among these programs, pancreatic cancers preferentially utilize a metabolic pathway through the enzyme aspartate aminotransferase 1 [also known as glutamate oxaloacetate transaminase 1 (GOT1)] to support cellular redox homeostasis. As such, small molecule inhibitors that target GOT1 could serve as starting points for the development of new therapies for pancreatic cancer. We ran a high-throughput screen for inhibitors of GOT1 and identified a small molecule, iGOT1-01, with in vitro GOT1 inhibitor activity. Application in pancreatic cancer cells revealed metabolic and growth inhibitory activity reflecting a promiscuous inhibitory profile. We then performed an in silico docking analysis to study inhibitor–GOT1 interactions with iGOT1-01 analogues that possess improved solubility and potency properties. Here, these results suggested that the GOT1 inhibitor competed for binding to the pyridoxal 5-phosphate (PLP) cofactor site of GOT1. To analyze how the GOT1 inhibitor bound to GOT1, a series of GOT1 mutant enzymes that abolished PLP binding were generated. Application of the mutants in X-ray crystallography and thermal shift assays again suggested but were unable to formally conclude that the GOT1 inhibitor bound to the PLP site. Mutational studies revealed the relationship between PLP binding and the thermal stability of GOT1 while highlighting the essential nature of several residues for GOT1 catalytic activity. Insight into the mode of action of GOT1 inhibitors may provide leads to the development of drugs that target redox balance in pancreatic cancer.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Structural and Functional Characterization of YdjI, an Aldolase of Unknown Specificity in Escherichia coli K12

The ydj gene cluster is found in 80% of sequenced Escherichia coli genomes and other closely related species in the human microbiome. On the basis of the annotations of the enzymes located in this cluster, it is expected that together they catalyze the catabolism of an unknown carbohydrate. The focus of this investigation is on YdjI, which is in the ydj gene cluster of E. coli K-12. It is predicted to be a class II aldolase of unknown function. Here we describe a structural and functional characterization of this enzyme. YdjI catalyzes the hydrogen/deuterium exchange of the pro- S hydrogen at C3 of dihydroxyacetone phosphate (DHAP). In the presence of DHAP, YdjI catalyzes an aldol condensation with a variety of aldo sugars. YdjI shows a strong preference for higher-order (seven-, eight-, and nine-carbon) monosaccharides with specific hydroxyl stereochemistries and a negatively charged terminus (carboxylate or phosphate). The best substrate is l -arabinuronic acid with an apparent k cat of 3.0 s –1 . The product, l - glycero - l - galacto -octuluronate-1-phosphate, has a k cat / K m value of 2.1 × 10 3 M –1 s –1 in the retro-aldol reaction with YdjI. This is the first recorded synthesis of l - glycero - l - galacto -octuluronate-1-phosphate and six similar carbohydrates. The crystal structure of YdjI, determined to a nominal resolution of 1.75 Å (Protein Data Bank entry 6OFU ), reveals unusual positions for two arginine residues located near the active site. Furthermore, computational docking was utilized to distinguish preferable binding orientations for l - glycero - l - galacto -octuluronate-1-phosphate. These results indicate a possible alternative binding orientation for l - glycero - l - galacto -octuluronate-1-phosphate compared to that observed in other class II aldolases, which utilize shorter carbohydrate molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Structural, Chemical, and Dynamical Frustration: Origins of Superionic Conductivity in closo -Borate Solid Electrolytes

Li 2 B 12 H 12 , Na 2 B 12 H 12 , and their closo-borate relatives exhibit unusually high ionic conductivity, making them attractive as a new class of candidate electrolytes in solid-state Li- and Na-ion batteries. However, further optimization of these materials requires a deeper understanding of the fundamental mechanisms underlying ultrafast ion conduction. To this end, we use ab initio molecular dynamics simulations and density-functional calculations to explore the motivations for cation diffusion. We find that superionic behavior in Li 2 B 12 H 12 and Na 2 B 12 H 12 results from a combination of key structural, chemical, and dynamical factors that introduce intrinsic frustration and disorder. A statistical metric is used to show that the structures exhibit a high density of accessible interstitial sites and site types, which correlates with the flatness of the energy landscape and the observed cation mobility. Furthermore, cations are found to dock to specific anion sites, leading to a competition between the geometric symmetry of the anion and the symmetry of the lattice itself, which can facilitate cation hopping. Finally, facile anion reorientations and other low-frequency thermal vibrations lead to fluctuations in the local potential that enhance cation mobility by creating a local driving force for hopping. In conclusion, we discuss the relevance of each factor for developing new ionic conductivity descriptors that can be used for discovery and optimization of closo-borate solid electrolytes, as well as superionic conductors more generally.

25 ENERGY STORAGE↗
DOE OSTI2016

Structure–Function Studies of Naphthalene, Phenanthrene, Biphenyl, and Their Derivatives in Interaction with and Oxidation by Cytochromes P450 2A13 and 2A6

In this study, naphthalene, phenanthrene, biphenyl, and their derivatives having different ethynyl, propynyl, butynyl, and propargyl ether substitutions were examined for their interaction with and oxidation by cytochromes P450 (P450) 2A13 and 2A6. Spectral interaction studies suggested that most of these chemicals interacted with P450 2A13 to induce Type I binding spectra more readily than with P450 2A6. Among the various substituted derivatives examined, 2-ethynylnaphthalene, 2-naphthalene propargyl ether, 3-ethynylphenanthrene, and 4-biphenyl propargyl ether had larger ΔA max /K s values in inducing Type I binding spectra with P450 2A13 than their parent compounds. P450 2A13 was found to oxidize naphthalene, phenanthrene, and biphenyl to 1-naphthol, 9-hydroxyphenanthrene, and 2- and/or 4-hydroxybiphenyl, respectively, at much higher rates than P450 2A6. Other human P450 enzymes including P450s 1A1, 1A2, 1B1, 2C9, and 3A4 had lower rates of oxidation of naphthalene, phenanthrene, and biphenyl than P450s 2A13 and 2A6. Additionally, those alkynylated derivatives that strongly induced Type I binding spectra with P450s 2A13 and 2A6 were extensively oxidized by these enzymes upon analysis with HPLC. Molecular docking studies supported the hypothesis that ligand-interaction energies (U values) obtained with reported crystal structures of P450 2A13 and 2A6 bound to 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone, indole, pilocarpine, nicotine, and coumarin are of use in understanding the basis of possible molecular interactions of these xenobiotic chemicals with the active sites of P450 2A13 and 2A6 enzymes. In fact, the ligand-interaction energies with P450 2A13 4EJG bound to these chemicals were found to relate to their induction of Type I binding spectra.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

In-Use Passenger Vessel Emission Rates of Black Carbon and Nitrogen Oxides

This study quantified emission factors of black carbon (BC) and nitrogen oxides (NO x ) from 21 engines on in-use excursion vessels and ferries operating in California's San Francisco Bay, including EPA uncertified and Tier 1-4 engines and across engine operating modes. On average, ~60 fuel-based emission factors per engine were measured using a novel combination of exhaust plume capture combined with GPS location and speed data that can be more readily deployed than common portable emissions measurement systems. Here, BC and NO x emission factors (g kg -1 ) were lowest and least variable during fast cruising and highest during maneuvering and docked operation. Selective catalytic reduction (SCR) reduced NO x emissions by ~80% when functional. However, elevated NO x emissions that exceeded corresponding exhaust standards were measured on most Tier 3 and Tier 4 engines sampled, which can be attributed to inactive SCR during frequent low engine load operation. In contrast, BC emissions exceeded the PM emission standard for only one engine, and SCR systems employed as a NO x reduction technology also reduced emitted BC. Using these measured emission factors to compare commuting options, we show that the CO 2 -equivalent emissions per passenger-kilometer are comparable when commuting by car and ferry, but BC and NO x emissions can be several to more than ten times larger when commuting by ferry.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Paired Simulations and Experimental Investigations into the Calcium-Dependent Conformation of Albumin

Serum albumin is the most abundant protein in blood plasma, and it is involved in multiple biological processes. Serum albumin has recently been adapted for improving biomaterial integration with bone tissue, and studies have shown the importance of this protein in bone repair and regeneration. However, the mechanism of action is not yet clear. In stark contrast, other studies have demonstrated that albumin blocks cell adhesion to surfaces, which is seen as a limitation to its bone healing role. These apparent contradictions suggest that the conformation of albumin facilitates its bioactivity, leading to enhanced bone repair. Serum albumin is known to play a major role in maintaining the calcium ion concentration in blood plasma. Due to the prevalence of calcium at bone repair and regeneration sites, it has been hypothesized that calcium binding to serum albumin triggers a conformational change, leading to bioactivity. In the current study, molecular modeling approaches including molecular docking, atomic molecular dynamics (MD) simulation, and coarse-grained MD simulation were used to test this hypothesis by investigating the conformational changes induced in bovine serum albumin by interaction with calcium ions. The computational results were qualitatively validated with experimental Fourier-transform infrared spectroscopy analysis. We find that free calcium ions in solution transiently bind with the three major loops in albumin, triggering a conformational change where N-terminal and C-terminal domains separate from each other in a partial unfolding process. The separation distance between these domains was found to correlate with the calcium ion concentration. The experimental data support the simulation results showing that albumin has enhanced conformational heterogeneity upon exposure to intermediate levels of calcium, without any significant secondary structure changes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

CoarsenConf: Equivariant Coarsening with Aggregated Attention for Molecular Conformer Generation

Molecular conformer generation (MCG) is an important task in cheminformatics and drug discovery. The ability to efficiently generate low-energy 3D structures can avoid expensive quantum mechanical simulations, leading to accelerated virtual screenings and enhanced structural exploration. Several generative models have been developed for MCG, but many struggle to consistently produce high-quality conformers for meaningful downstream applications. To address these issues, we introduce CoarsenConf, which coarse-grains molecular graphs based on torsional angles and integrates them into an SE(3)-equivariant hierarchical variational autoencoder. Through equivariant coarse-graining, we aggregate the fine-grained atomic coordinates of subgraphs connected via rotatable bonds, creating a variable-length coarse-grained latent representation. Our model uses a novel aggregated attention mechanism to restore fine-grained coordinates from the coarse-grained latent representation, enabling efficient generation of accurate conformers. Furthermore, we evaluate the chemical and biochemical quality of our generated conformers on multiple downstream applications, including property prediction and large-scale oracle-based protein docking. Overall, CoarsenConf generates more accurate conformer ensembles compared to prior generative models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Comparative Analysis of TCR and TCR-pMHC Complex Structure Prediction Tools

The rapid development of computational approaches for predicting the structures of T cell receptors (TCRs) and TCR-peptide-major histocompatibility (TCR-pMHC) complexes, accelerated by AI breakthroughs such as AlphaFold, has made it feasible to calculate these structures with increasing accuracy. Although these tools show great potential, their relative accuracy and limitations remain unclear due to the lack of standardized benchmarks. Here, we systematically evaluate seven tools for predicting isolated TCR structures together with six tools for predicting TCR-pMHC complex structures. The methods include homology-based approaches, general prediction tools using AlphaFold, TCR-specific tools derived from AlphaFold2, and the newly developed tFold-TCR model. The evaluation uses a post-training data set comprising 40 αβ TCRs and 27 TCR-pMHC complexes (21 Class I and 6 Class II). Model accuracy is assessed at global, local, and interface levels using a variety of metrics. We find that each tool offers distinct advantages in various aspects of its predictions. AlphaFold2, AlphaFold3, and tFold-TCR excel in overall accuracy of TCR structure prediction, and TCRmodel2 and AlphaFold2 perform well in overall accuracy of TCR-pMHC structure prediction. However, TCR-specific tools derived from AlphaFold2 show lower accuracy in the framework region than both homology-based methods and general-purpose tools such as AlphaFold, and challenges remain for all in modeling CDR3 loops, docking orientations, TCR-peptide interfaces, and Class II MHC-peptide interfaces. Furthermore, these findings will guide researchers in selecting appropriate tools, emphasize the importance of using multiple evaluation metrics to assess model performance, and offer suggestions for improving TCR and TCR-pMHC structure prediction tools.

Chemical structure↗
DOE OSTI2025

Design, Synthesis, and Biological Activity of Novel Ornithine Decarboxylase (ODC) Inhibitors

We here describe the design, synthesis, and biological activity of novel ornithine decarboxylase (ODC) inhibitors that show significantly higher potency in vitro than α-difluoromethylornithine (DFMO), a U.S. Food and Drug Administration (FDA) approved drug. We report two X-ray structures of ODC complexed with new ODC inhibitors, computational docking, molecular dynamics, and binding free energy calculations to validate the experimental models. The X-ray structures reveal that covalent adducts with pyridoxal phosphate (PLP) are formed in the active site of the human ODC enzyme, as verified by their preparation and enzymatic testing. Finally, we verified that the cellular activity of endogenous ODC was inhibited, and polyamine levels were reduced. Given that ODC is a clinically validated target, combined with the fact that DFMO is currently the only ODC inhibitor in clinical use for several indications, the further development of more potent ODC inhibitors with superior activity and physical properties is warranted.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Discovery of Tarantula Venom-Derived Na V 1.7-Inhibitory JzTx-V Peptide 5-Br-Trp24 Analogue AM-6120 with Systemic Block of Histamine-Induced Pruritis

Inhibitors of the voltage-gated sodium channel Na V 1.7 are being investigated as pain therapeutics due to compelling human genetics. We previously identified Na V 1.7-inhibitory peptides GpTx-1 and JzTx-V from tarantula venom screens. Potency and selectivity were modulated through attribute-based positional scans of native residues via chemical synthesis. In this paper, we report JzTx-V lead optimization to identify a pharmacodynamically active peptide variant. Molecular docking of peptide ensembles from NMR into a homology model-derived Na V 1.7 structure supported prioritization of key residues clustered on a hydrophobic face of the disulfide-rich folded peptide for derivatization. Replacing Trp24 with 5-Br- Trp24 identified lead peptides with activity in electrophysiology assays in engineered and neuronal cells. 5-Br-Trp24 containing peptide AM-6120 was characterized in X-ray crystallography and pharmacokinetic studies and blocked histamine-induced pruritis in mice after subcutaneous administration, demonstrating systemic Na V 1.7-dependent pharmacodynamics. Our data suggests a need for high target coverage based on plasma exposure for impacting in vivo end points with selectivity-optimized peptidic Na V 1.7 inhibitors.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

3D-Scaffold: A Deep Learning Framework to Generate 3D Coordinates of Drug-like Molecules with Desired Scaffolds

The prerequisite of therapeutic drug design is to identify novel molecules with desired biophysical and biochemical properties. Deep generative models have demonstrated their ability to find such molecules by exploring a huge chemical space efficiently. An effective way to obtain molecules with desired target properties is the preservation of critical scaffolds in the generation process. To this end, we propose a domain aware generative framework called 3D-Scaffold that takes 3D coordinates of a desired scaffold as an input and generates 3D coordinates of novel therapeutic candidates as an output while always preserving the desired scaffolds in generated structures. We show that our framework generates predominantly valid, unique, novel, and experimentally synthesizable molecules that have drug-like properties similar to the molecules in the training set. Using domain specific datasets, we generate covalent and non-covalent antiviral inhibitors. Therefore, to measure the success of our framework in generating therapeutic candidates, generated structures were subjected to high throughput virtual screening via docking simulations, which shows favorable interaction against SARS-CoV-2 main protease and non-structural protein endoribonuclease (NSP15) targets. Most importantly, our model performs well with relatively small volumes of training data and generalizes to new scaffolds, making it applicable to other domain.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Deciphering the Cofilin Oligomers via Intermolecular Disulfide Bond Formation: A Coarse-Grained Molecular Dynamics Approach to Understanding Cofilin’s Regulation on Actin Filaments

Cofilin, a key actin-binding protein, orchestrates the dynamics of the actomyosin network through its actin-severing activity and by promoting the recycling of actin monomers. Recent experimental work suggests that cofilin also forms functionally distinct oligomers through thiol post-translational modification (PTM) that encourages actin nucleation and assembly. Despite these advances, the structural conformations of cofilin oligomers that modulate actin activity remain elusive because there are combinatorial ways to oxidize thiols in cysteines to form disulfide bonds rapidly. This study employs molecular dynamics simulations to investigate human cofilin 1 as a case study for exploring cofilin dimers via disulfide bond formation. Using the free energy profiling, our simulations unveil a range of probable cofilin dimer structures not represented in current Protein Data Bank entries. These candidate dimers are characterized by their distinct population distributions and relative free energies. Of particular note is a dimer featuring an interface between cysteines 139 and 147 residues, which demonstrates stable free energy characteristics and intriguingly symmetrical geometry. In contrast, the experimentally proposed dimer structure exhibits a less stable free energy profile. Here, we also evaluate frustration quantification based on the energy landscape theory in the protein-protein interactions at the dimer interfaces. Notably, the 39-39 dimer configuration emerges as a promising candidate for forming cofilin tetramers, as substantiated by frustration analysis. Additionally, docking simulations with actin filaments further evaluate the stability of these cofilin dimer-actin complexes. Our findings thus offer a computational framework for understanding the role of thiol post-translational modification cofilin proteins in regulating oligomerization, and the subsequent cofilin-mediated actin dynamics in the actomyosin network.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Revealing Surface Waters on an Antifreeze Protein by Fusion Protein Crystallography Combined with Molecular Dynamic Simulations

Antifreeze proteins (AFPs) adsorb to ice through an extensive, flat, relatively hydrophobic surface. It has been suggested that this ice-binding site (IBS) organizes surface waters into an ice-like clathrate arrangement that matches and fuses to the quasi-liquid layer on the ice surface. On cooling, these waters join the ice lattice and freeze the AFP to its ligand. Evidence for the generality of this binding mechanism is limited because AFPs tend to crystallize with their IBS as a preferred protein–protein contact surface, which displaces some bound waters. Type III AFP is a 7 kDa globular protein with an IBS made up two adjacent surfaces. In the crystal structure of the most active isoform (QAE1), the part of the IBS that docks to the primary prism plane of ice is partially exposed to solvent and has clathrate waters present that match this plane of ice. The adjacent IBS, which matches the pyramidal plane of ice, is involved in protein–protein crystal contacts with few surface waters. Here we have changed the protein–protein contacts in the ice-binding region by crystallizing a fusion of QAE1 to maltose-binding protein. In this 1.9 Å structure, the IBS that fits the pyramidal plane of ice is exposed to solvent. By combining crystallography data with MD simulations, the surface waters on both sides of the IBS were revealed and match well with the target ice planes. Here, the waters on the pyramidal plane IBS were loosely constrained, which might explain why other isoforms of type III AFP that lack the prism plane IBS are less active than QAE1. The AFP fusion crystallization method can potentially be used to force the exposure to solvent of the IBS on other AFPs to reveal the locations of key surface waters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Relative binding affinities of monolignols to horseradish peroxidase

Monolignol binding to the peroxidase active site is the first step in lignin polymerization in plant cell walls. Using molecular dynamics, docking, and free energy perturbation calculations, we investigate the binding of monolignols to horseradish peroxidase C. Our results suggest that p-coumaryl alcohol has the strongest binding affinity followed by sinapyl and coniferyl alcohol. Stacking interactions between the monolignol aromatic rings and nearby phenylalanine residues play an important role in determining the calculated relative binding affinities. p-Coumaryl and coniferyl alcohols bind in a pose productive for reaction in which a direct H-bond is formed between the phenolic –OH group and a water molecule (W2) that may facilitate proton transfer during oxidation. In contrast, in the case of sinapyl alcohol there is no such direct interaction, the phenolic –OH group instead interacting with Pro139. Furthermore, since proton and electron transfer is the rate-limiting step in monolignol oxidation by peroxidase, the binding pose (and thus the formation of near attack conformation) appears to play a more important role than the overall binding affinity in determining the oxidation rate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

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.

Schwierz, Nadine↗