IKKα-deficient lung adenocarcinomas generate an immunosuppressive microenvironment by overproducing Treg-inducing cytokines.

The tumor microenvironment (TME) provides potential targets for cancer therapy. However, how signals originating in cancer cells affect tumor-directed immunity is largely unknown. Deletions in the CHUK locus, coding for IκB kinase α (IKKα), correlate with reduced lung adenocarcinoma (ADC) patient survival and promote KrasG12D-initiated ADC development in mice, but it is unknown how reduced IKKα expression affects the TME. Here, we report that low IKKα expression in human and mouse lung ADC cells correlates with increased monocyte-derived macrophage and regulatory T cell (Treg) scores and elevated transcription of genes coding for macrophage-recruiting and Treg-inducing cytokines (CSF1, CCL22, TNF, and IL-23A). By stimulating recruitment of monocyte-derived macrophages from the bone marrow and enforcing a TNF/TNFR2/c-Rel signaling cascade that stimulates Treg generation, these cytokines promote lung ADC progression. Depletion of TNFR2, c-Rel, or TNF in CD4+ T cells or monocyte-derived macrophages dampens Treg generation and lung tumorigenesis. Treg depletion also attenuates carcinogenesis. In conclusion, reduced cancer cell IKKα activity enhances formation of a protumorigenic TME through a pathway whose constituents may serve as therapeutic targets for KRAS-initiated lung ADC.

Accelerating antibody discovery and design with artificial intelligence: Recent advances and prospects.

Therapeutic antibodies are the largest class of biotherapeutics and have been successful in treating human diseases. However, the design and discovery of antibody drugs remains challenging and time-consuming. Recently, artificial intelligence technology has had an incredible impact on antibody design and discovery, resulting in significant advances in antibody discovery, optimization, and developability. This review summarizes major machine learning (ML) methods and their applications for computational predictors of antibody structure and antigen interface/interaction, as well as the evaluation of antibody developability. Additionally, this review addresses the current status of ML-based therapeutic antibodies under preclinical and clinical phases. While many challenges remain, ML may offer a new therapeutic option for the future direction of fully computational antibody design.

Computational investigation of the inhibitory interaction of IRF3 and SARS-CoV-2 accessory protein ORF3b.

ORF3b is one of the SARS-CoV-2 accessory proteins. Previous experimental study suggested that ORF3b prevents IRF3 translocating to nucleus. However, the biophysical mechanism of ORF3b-IRF3 interaction is elusive. Here, we explored the conformation ensemble of ORF3b using all-atom replica exchange molecular dynamics simulation. Disordered ORF3b has mixed α-helix, ß-turn and loop conformers. The potential ORF3b-IRF3 binding modes were searched by docking representative ORF3b conformers with IRF3, and 50 ORF3b-IRF3 complex poses were screened using molecular dynamics simulations ranging from 500 to 1000 ns. We found that ORF3b binds IRF3 predominantly on its CBP binding and phosphorylated pLxIS motifs, with CBP binding site has the highest binding affinity. The ORF3b-IRF3 binding residues are highly conserved in SARS-CoV-2. Our results provided biophysics insights into ORF3b-IRF3 interaction and explained its interferon antagonism mechanism.

Re-Engineering Therapeutic Anti-Aß Monoclonal Antibody to Target Amyloid Light Chain.

Amyloidosis involves the deposition of misfolded proteins. Even though it is caused by different pathogenic mechanisms, in aggregate, it shares similar features. Here, we tested and confirmed a hypothesis that an amyloid antibody can be engineered by a few mutations to target a different species. Amyloid light chain (AL) and ß-amyloid peptide (Aß) are two therapeutic targets that are implicated in amyloid light chain amyloidosis and Alzheimer's disease, respectively. Though crenezumab, an anti-Aß antibody, is currently unsuccessful, we chose it as a model to computationally design and prepare crenezumab variants, aiming to discover a novel antibody with high affinity to AL fibrils and to establish a technology platform for repurposing amyloid monoclonal antibodies. We successfully re-engineered crenezumab to bind both Aß42 oligomers and AL fibrils with high binding affinities. It is capable of reversing Aß42-oligomers-induced cytotoxicity, decreasing the formation of AL fibrils, and alleviating AL-fibrils-induced cytotoxicity in vitro. Our research demonstrated that an amyloid antibody could be engineered by a few mutations to bind new amyloid sequences, providing an efficient way to reposition a therapeutic antibody to target different amyloid diseases.

Precision Therapy of Recurrent Breast Cancer through Targeting Different Malignant Tumor Cells with a HER2/CD44-Targeted Hydrogel Nanobot.

Heterogeneity and drug resistance of tumor cells are the leading causes of incurability and poor survival for patients with recurrent breast cancer. In order to accurately deliver the biological anticancer drugs to different subtypes of malignant tumor cells for omnidirectional targeted treatment of recurrent breast cancer, a distinct design is demonstrated by embedding liposome-based nanocomplexes containing pro-apoptotic peptide and survivin siRNA drugs (LPR) into Herceptin/hyaluronic acid cross-linked nanohydrogels (Herceptin-HA) to fabricate a HER2/CD44-targeted hydrogel nanobot (named as ALPR). ALPR delivered cargoes to the cells overexpressing CD44 and HER2, followed by Herceptin-HA biodegradation, subsequently, the exposed lipid component containing DOPE fused with the endosomal membrane and released peptide and siRNA into the cytoplasm. These experiments indicated that ALPR can specifically deliver Herceptin, peptide, and siRNA drugs to HER2-positive SKBR-3, triple-negative MDA-MB-231, and HER2-negative drug-resistant MCF-7 human breast cancer cells. ALPR completely inhibited the growth of heterogeneous breast tumors via multichannel synergistic effects: disrupting mitochondria, downregulating the survivin gene, and blocking HER2 receptors on the surface of HER2-positive cells. The present design overcomes the chemical drug resistance and opens a feasible route for the combinative treatment of recurrent breast cancer, even other solid tumors, utilizing different kinds of biological drugs.

Amyloid Oligomers: A Joint Experimental/Computational Perspective on Alzheimer's Disease, Parkinson's Disease, Type II Diabetes, and Amyotrophic Lateral Sclerosis.

Protein misfolding and aggregation is observed in many amyloidogenic diseases affecting either the central nervous system or a variety of peripheral tissues. Structural and dynamic characterization of all species along the pathways from monomers to fibrils is challenging by experimental and computational means because they involve intrinsically disordered proteins in most diseases. Yet understanding how amyloid species become toxic is the challenge in developing a treatment for these diseases. Here we review what computer, in vitro, in vivo, and pharmacological experiments tell us about the accumulation and deposition of the oligomers of the (Aß, tau), α-synuclein, IAPP, and superoxide dismutase 1 proteins, which have been the mainstream concept underlying Alzheimer's disease (AD), Parkinson's disease (PD), type II diabetes (T2D), and amyotrophic lateral sclerosis (ALS) research, respectively, for many years.

Data Mining Suggests That CXCL14 Gene Silencing in Colon Cancer Is Due to Promoter Methylation.

CXCL14 is one of the most evolutionarily conserved members of the chemokine family and is constitutionally expressed in multiple organs, suggesting that it is involved in the homeostasis maintenance of the system. CXCL14 is highly expressed in colon epithelial cells and shows obvious gene silencing in clinical colon cancer samples, suggesting that its silencing is related to the immune escape of cancer cells. In this paper, we analyzed the expression profiles of multiple human clinical colon cancer datasets and mouse colon cancer models to reveal the variation trend of CXCL14 expression during colitis, colon polyps, primary colon cancer, and liver metastases. The relationship between CXCL14 gene silencing and promoter hypermethylation was revealed through the colorectal carcinoma methylation database. The results suggest that CXCL14 is a tumor suppressor gene in colorectal carcinoma which is activated first and then silenced during the process of tumor occurrence and deterioration. Promoter hypermethylation is the main cause of CXCL14 silencing. The methylation level of CXCL14 is correlated with the anatomic site of tumor occurrence, positively correlated with patient age, and associated with prognosis. Reversing the hypermethylation of CXCL14 may be an epigenetic therapy for colon cancer.

IKKα inactivation promotes Kras-initiated lung adenocarcinoma development through disrupting major redox regulatory pathways.

Lung adenocarcinoma (ADC) and squamous cell carcinoma (SCC) are two distinct and predominant types of human lung cancer. IκB kinase α (IKKα) has been shown to suppress lung SCC development, but its role in ADC is unknown. We found inactivating mutations and homologous or hemizygous deletions in the CHUK locus, which encodes IKKα, in human lung ADCs. The CHUK deletions significantly reduced the survival time of patients with lung ADCs harboring KRAS mutations. In mice, lung-specific Ikkα ablation (IkkαΔLu ) induces spontaneous ADCs and promotes KrasG12D-initiated ADC development, accompanied by increased cell proliferation, decreased cell senescence, and reactive oxygen species (ROS) accumulation. IKKα deletion up-regulates NOX2 and down-regulates NRF2, leading to ROS accumulation and blockade of cell senescence induction, which together accelerate ADC development. Pharmacologic inhibition of NADPH oxidase or ROS impairs KrasG12D-mediated ADC development in IkkαΔLu mice. Therefore, IKKα modulates lung ADC development by controlling redox regulatory pathways. This study demonstrates that IKKα functions as a suppressor of lung ADC in human and mice through a unique mechanism that regulates tumor cell-associated ROS metabolism.

Conformational stability and dynamics of the cancer-associated isoform Δ133p53ß are modulated by p53 peptides and p53-specific DNA.

p53 is a tumor suppressor protein that maintains genome stability, but its Δ133p53ß and Δ160p53ß isoforms promote breast cancer cell invasion. The sequence truncations in the p53 core domain raise key questions related to their physicochemical properties, including structural stabilities, interaction mechanisms, and DNA-binding abilities. Herein, we investigated the conformational dynamics of Δ133p53ß and Δ160p53ß with and without binding to p53-specific DNA by using molecular dynamics simulations. We observed that the core domains of the 2 truncated isoforms are much less stable than wild-type (wt) p53ß, and the increased solvent exposure of their aggregation-triggering segment indicates their higher aggregation propensities than wt p53. We also found that Δ133p53ß stability is modulable by peptide or DNA interactions. Adding a p53 peptide (derived from truncated p53 sequence 107-129) may help stabilize Δ133p53. Most importantly, our simulations of p53 isomer-DNA complexes indicate that Δ133p53ß dimer, but not Δ160p53ß dimer, could form a stable complex with p53-specific DNA, which is consistent with recent experiments. This study provides physicochemical insight into Δ133p53ß, Δ133p53ß-DNA complexes, Δ133p53ß's pathologic mechanism, and peptide-based inhibitor design against p53-related cancers.-Lei, J., Qi, R., Tang, Y., Wang, W., Wei, G., Nussinov, R., Ma, B. Conformational stability and dynamics of the cancer-associated isoform Δ133p53ß are modulated by p53 peptides and p53-specific DNA.

Structural disorder in four-repeat Tau fibrils reveals a new mechanism for barriers to cross-seeding of Tau isoforms.

The intracellular deposition of fibrils composed of the microtubule-associated protein Tau is a characteristic feature of Alzheimer's disease (AD) and other fatal neurodegenerative disorders collectively known as tauopathies. Short Tau fibrils spread intracerebrally through transfer between interconnected neurons. Once taken up by a recipient cell, Tau fibrils recruit Tau monomers onto their ends. Based on the number of microtubule-binding repeats, there are two distinct groups of Tau isoforms: three-repeat (3R) Tau and four-repeat (4R) Tau. In AD, all Tau isoforms are deposited, whereas in other tauopathies, only 3R or 4R Tau isoforms are deposited. The molecular basis for these isoform-specific depositions is poorly understood, although conformation-based cross-seeding barriers are key. Here, we used sedimentation assays, EPR spectroscopy, and other structural readouts to better understand the cross-seeding barriers of 4R Tau fibrils. We observed that fibrils formed from truncated Tau (K18), but not full-length Tau (htau40), exhibit a barrier that inhibits 3R Tau recruitment. Investigating an array of differently sized fragments, we found that the Tau C terminus modulates the cross-seeding barrier and that the N terminus plays a synergistic role. Two disease-associated Tau variants, P301S and P301L, also established strong cross-seeding barriers. EPR analysis indicated that fibrils seeded with truncated and mutated Tau, but not htau40, are structurally disordered in the second half of repeat four and onward. These findings suggest that the disorder in this region diminishes the ability of 4R Tau fibrils to recruit 3R Tau monomers, revealing a new mechanism for Tau cross-seeding barriers.

Peptide-MHC (pMHC) binding to a human antiviral T cell receptor induces long-range allosteric communication between pMHC- and CD3-binding sites.

T cells generate adaptive immune responses mediated by the T cell receptor (TCR)-CD3 complex comprising an αß TCR heterodimer noncovalently associated with three CD3 dimers. In early T cell activation, αß TCR engagement by peptide-major histocompatibility complex (pMHC) is first communicated to the CD3 signaling apparatus of the TCR-CD3 complex, but the underlying mechanism is incompletely understood. It is possible that pMHC binding induces allosteric changes in TCR conformation or dynamics that are then relayed to CD3. Here, we carried out NMR analysis and molecular dynamics (MD) simulations of both the α and ß chains of a human antiviral TCR (A6) that recognizes the Tax antigen from human T cell lymphotropic virus-1 bound to the MHC class I molecule HLA-A2. We observed pMHC-induced NMR signal perturbations in the TCR variable (V) domains that propagated to three distinct sites in the constant (C) domains: 1) the Cß FG loop projecting from the Vß/Cß interface; 2) a cluster of Cß residues near the Cß αA helix, a region involved in interactions with CD3; and 3) the Cα AB loop at the membrane-proximal base of the TCR. A biological role for each of these allosteric sites is supported by previous mutational and functional studies of TCR signaling. Moreover, the pattern of long-range, ligand-induced changes in TCR A6 revealed by NMR was broadly similar to that predicted by the MD simulations. We propose that the unique structure of the TCR ß chain enables allosteric communication between the TCR-binding sites for pMHC and CD3.

Molecular dynamics based improvement of the solubilizing self-cleavable tag Zbasic-ΔI-CM application in the preparation of recombinant proteins in Escherichia coli.

Zbasic-ΔI-CM is a novel intein-based self-cleavable tag we developed to accelerate the soluble expression of recombinant proteins in Escherichia coli (E. coli). Previously we found that intein activity could be interfered by its flanking exteins, and thus reducing the production efficiency and final yield. In this work, we used CXC-chemokine 9 (CXCL9) as a model C-extein, which fusion with Zbasic-ΔI-CM showed high intein activity. When the fusion protein got soluble expression, CXCL9 was released immediately and purified directly from cell lysis supernatant. The results demonstrated that Zbasic-ΔI-CM tag had successfully mediated the efficient production of high-quality CXCL9 with reduced time and resources consumption in comparison with inclusion bodies expression. Molecular dynamics simulations suggested that the improved cleavage activity of Zbasic-ΔI-CM upon fusion with CXCL9 may be due to the higher dynamics of the first half loop and stabilization of the second half loop of intein. Our results proved that the self-cleavable Zbasic-ΔI-CM mediated soluble expression could be a feasible process for cytokines like CXCL9, thus of attractive potentials for production of therapeutic proteins using E. coli expression system.

Mechanisms of recognition of amyloid-ß (Aß) monomer, oligomer, and fibril by homologous antibodies.

Alzheimer's disease is one of the most devastating neurodegenerative diseases without effective therapies. Immunotherapy is a promising approach, but amyloid antibody structural information is limited. Here we simulate the recognition of monomeric, oligomeric, and fibril amyloid-ß (Aß) by three homologous antibodies (solanezumab, crenezumab, and their chimera, CreneFab). Solanezumab only binds the monomer, whereas crenezumab and CreneFab can recognize different oligomerization states; however, the structural basis for this observation is not understood. We successfully identified stable complexes of crenezumab with Aß pentamer (oligomer model) and 16-mer (fibril model). It is noteworthy that solanezumab targets Aß residues 16-26 preferentially in the monomeric state; conversely, crenezumab consistently targets residues 13-16 in different oligomeric states. Unlike the buried monomeric peptide in solanezumab's complementarity-determining region, crenezumab binds the oligomer's lateral and edge residues. Surprisingly, crenezumab's complementarity-determining region loops can effectively bind the Aß fibril lateral surface around the same 13-16 region. The constant domain influences antigen recognition through entropy redistribution. Different constant domain residues in solanezumab/crenezumab/chimera influence the binding of Aß aggregates. Collectively, we provide molecular insight into the recognition mechanisms facilitating antibody design.

Structure and energetic basis of overrepresented λ light chain in systemic light chain amyloidosis patients.

Amyloid formation and deposition of immunoglobulin light-chain proteins in systemic amyloidosis (AL) cause major organ failures. While the κ light-chain is dominant (λ/κ=1:2) in healthy individuals, λ is highly overrepresented (λ/κ=3:1) in AL patients. The structural basis of the amyloid formation and the sequence preference are unknown. We examined the correlation between sequence and structural stability of dimeric variable domains of immunoglobulin light chains using molecular dynamics simulations of 24 representative dimer interfaces, followed by energy evaluation of conformational ensembles for 20 AL patients' light chain sequences. We identified a stable interface with displaced N-terminal residues, provides the structural basis for AL protein fibrils formation. Proline isomerization may cause the N-terminus to adopt amyloid-prone conformations. We found that λ light-chains prefer misfolded dimer conformation, while κ chain structures are stabilized by a natively folded dimer. Our study may facilitate structure-based small molecule and antibody design to inhibit AL. This article is part of a Special Issue entitled: Accelerating Precision Medicine through Genetic and Genomic Big Data Analysis edited by Yudong Cai & Tao Huang.

Local and global anatomy of antibody-protein antigen recognition.

Deciphering antibody-protein antigen recognition is of fundamental and practical significance. We constructed an antibody structural dataset, partitioned it into human and murine subgroups, and compared it with nonantibody protein-protein complexes. We investigated the physicochemical properties of regions on and away from the antibody-antigen interfaces, including net charge, overall antibody charge distributions, and their potential role in antigen interaction. We observed that amino acid preference in antibody-protein antigen recognition is entropy driven, with residues having low side-chain entropy appearing to compensate for the high backbone entropy in interaction with protein antigens. Antibodies prefer charged and polar antigen residues and bridging water molecules. They also prefer positive net charge, presumably to promote interaction with negatively charged protein antigens, which are common in proteomes. Antibody-antigen interfaces have large percentages of Tyr, Ser, and Asp, but little Lys. Electrostatic and hydrophobic interactions in the Ag binding sites might be coupled with Fab domains through organized charge and residue distributions away from the binding interfaces. Here we describe some features of antibody-antigen interfaces and of Fab domains as compared with nonantibody protein-protein interactions. The distributions of interface residues in human and murine antibodies do not differ significantly. Overall, our results provide not only a local but also a global anatomy of antibody structures.

Protein Ensembles: How Does Nature Harness Thermodynamic Fluctuations for Life? The Diverse Functional Roles of Conformational Ensembles in the Cell.

All soluble proteins populate conformational ensembles that together constitute the native state. Their fluctuations in water are intrinsic thermodynamic phenomena, and the distributions of the states on the energy landscape are determined by statistical thermodynamics; however, they are optimized to perform their biological functions. In this review we briefly describe advances in free energy landscape studies of protein conformational ensembles. Experimental (nuclear magnetic resonance, small-angle X-ray scattering, single-molecule spectroscopy, and cryo-electron microscopy) and computational (replica-exchange molecular dynamics, metadynamics, and Markov state models) approaches have made great progress in recent years. These address the challenging characterization of the highly flexible and heterogeneous protein ensembles. We focus on structural aspects of protein conformational distributions, from collective motions of single- and multi-domain proteins, intrinsically disordered proteins, to multiprotein complexes. Importantly, we highlight recent studies that illustrate functional adjustment of protein conformational ensembles in the crowded cellular environment. We center on the role of the ensemble in recognition of small- and macro-molecules (protein and RNA/DNA) and emphasize emerging concepts of protein dynamics in enzyme catalysis. Overall, protein ensembles link fundamental physicochemical principles and protein behavior and the cellular network and its regulation.

Conformational selection in amyloid-based immunotherapy: Survey of crystal structures of antibody-amyloid complexes.

BACKGROUND: The dominant feature in neurodegenerative diseases is protein aggregations that lead to neuronal loss. Immunotherapies using antibodies or antibody fragments to target the aggregations are a highly perused approach. The molecular mechanisms underlying the amyloid-based immunotherapy are complex. Deciphering the properties of amyloidogenic proteins responsible for these diseases is essential to obtain insights into antibody recognition of the amyloid antigens. SCOPE OF REVIEW: We systematically explore all available crystal structures of antibody-amyloid complexes related to neurodegenerative diseases, including antibodies that recognize the Aß peptide, tau protein, prion protein, alpha-synuclein, huntingtin protein (mHTT), and polyglutamine. MAJOR CONCLUSIONS: We found that antibodies mostly use the conformational selection mechanism to recognize the highly flexible amyloid antigens. In particular, solanezumab bound to Aß12-28 tripeptide motif conformation (F19F20A21), which is shared with the Aß42 fibril. This motif, which is trapped by the antibody, may provide the missing link in amyloid formation. Water molecules often bridge between the antibody and amyloid, contributing to the recognition. GENERAL SIGNIFICANCE: This paper provides the structural basis for antibody recognition of amyloidogenic proteins. The analysis and discussion of known structures are expected to help in the design and optimization of antibodies in neurodegenerative diseases. This article is part of a Special Issue entitled "System Genetics" Guest Editor: Dr. Yudong Cai and Dr. Tao Huang.

Dimerization of the SP1 Region of HIV-1 Gag Induces a Helical Conformation and Association into Helical Bundles: Implications for Particle Assembly.

UNLABELLED: HIV-1 immature particle (virus-like particle [VLP]) assembly is mediated largely by interactions between the capsid (CA) domains of Gag molecules but is facilitated by binding of the nucleocapsid (NC) domain to nucleic acid. We previously investigated the role of SP1, a "spacer" between CA and NC, in VLP assembly. We found that small changes in SP1 drastically disrupt assembly and that a peptide representing the sequence around the CA-SP1 junction is helical at high but not low concentrations. We suggested that by virtue of such a concentration-dependent change, this region could act as a molecular switch to activate HIV-1 Gag for VLP assembly. A leucine zipper domain can replace NC in Gag and still lead to the efficient assembly of VLPs. We find that SP1 mutants also disrupt assembly by these Gag-Zip proteins and have now studied a small fragment of this Gag-Zip protein, i.e., the CA-SP1 junction region fused to a leucine zipper. Dimerization of the zipper places SP1 at a high local concentration, even at low total concentrations. In this context, the CA-SP1 junction region spontaneously adopts a helical conformation, and the proteins associate into tetramers. Tetramerization requires residues from both CA and SP1. The data suggest that once this region becomes helical, its propensity to self-associate could contribute to Gag-Gag interactions and thus to particle assembly. There is complete congruence between CA/SP1 sequences that promote tetramerization when fused to zippers and those that permit the proper assembly of full-length Gag; thus, equivalent interactions apparently participate in VLP assembly and in SP1-Zip tetramerization. IMPORTANCE: Assembly of HIV-1 Gag into virus-like particles (VLPs) appears to require an interaction with nucleic acid, but replacement of its principal nucleic acid-binding domain with a dimerizing leucine zipper domain leads to the assembly of RNA-free VLPs. It has not been clear how dimerization triggers assembly. Results here show that the SP1 region spontaneously switches to a helical state when fused to a leucine zipper and that these helical molecules further associate into tetramers, mediated by interactions between hydrophobic faces of the helices. Thus, the correct juxtaposition of the SP1 region makes it "association competent." Residues from both capsid and SP1 contribute to tetramerization, while mutations disrupting proper assembly in Gag also prevent tetramerization. Thus, this region is part of an associating interface within Gag, and its intermolecular interactions evidently help stabilize the immature Gag lattice. These interactions are disrupted by proteolysis of the CA-SP1 junction during virus maturation.

Principles and Overview of Sampling Methods for Modeling Macromolecular Structure and Dynamics.

Investigation of macromolecular structure and dynamics is fundamental to understanding how macromolecules carry out their functions in the cell. Significant advances have been made toward this end in silico, with a growing number of computational methods proposed yearly to study and simulate various aspects of macromolecular structure and dynamics. This review aims to provide an overview of recent advances, focusing primarily on methods proposed for exploring the structure space of macromolecules in isolation and in assemblies for the purpose of characterizing equilibrium structure and dynamics. In addition to surveying recent applications that showcase current capabilities of computational methods, this review highlights state-of-the-art algorithmic techniques proposed to overcome challenges posed in silico by the disparate spatial and time scales accessed by dynamic macromolecules. This review is not meant to be exhaustive, as such an endeavor is impossible, but rather aims to balance breadth and depth of strategies for modeling macromolecular structure and dynamics for a broad audience of novices and experts.

Binding of protofibrillar Aß trimers to lipid bilayer surface enhances Aß structural stability and causes membrane thinning.

Alzheimer's disease, a common neurodegenerative disease, is characterized by the aggregation of amyloid-ß (Aß) peptides. The interactions of Aß with membranes cause changes in membrane morphology and ion permeation, which are responsible for its neurotoxicity and can accelerate fibril growth. However, the Aß-lipid interactions and how these induce membrane perturbation and disruption at the atomic level and the consequences for the Aß organization are not entirely understood. Here, we perform multiple atomistic molecular dynamics simulations on three protofibrillar Aß9-40 trimers. Our simulations show that, regardless of the morphologies and the initial orientations of the three different protofibrillar Aß9-40 trimers, the N-terminal ß-sheet of all trimers preferentially binds to the membrane surface. The POPG lipid bilayers enhance the structural stability of protofibrillar Aß trimers by stabilizing inter-peptide ß-sheets and D23-K28 salt-bridges. The interaction causes local membrane thinning. We found that the trimer structure related to Alzheimer's disease brain tissue () is the most stable both in water solution and at membrane surface, and displays slightly stronger membrane perturbation capability. These results provide mechanistic insights into the membrane-enhanced structural stability of protofibrillar Aß oligomers and the first step of Aß-induced membrane disruption at the atomic level.