Abstract
The microglial surface protein Triggering Receptor Expressed on Myeloid Cells 2 (TREM2) plays a critical role in mediating brain homeostasis and inflammatory responses in Alzheimer’s disease (AD). The soluble form of TREM2 (sTREM2) exhibits neuroprotective effects in AD, though the underlying mechanisms remain elusive. Moreover, differences in ligand binding between TREM2 and sTREM2, which have major implications for their roles in AD pathology, remain unexplained. To address these knowledge gaps, we conducted the most computationally intensive molecular dynamics simulations to date of (s)TREM2, exploring their interactions with key damage- and lipoprotein-associated phospholipids and the impact of the AD-risk mutation R47H. Our results demonstrate that the flexible stalk domain of sTREM2 serves as the molecular basis for differential ligand binding between sTREM2 and TREM2, facilitated by its role in modulating the dynamics of the Ig-like domain and altering the accessibility of canonical ligand binding sites. We identified a novel ligand binding site on sTREM2, termed the ‘Expanded Surface 2’, which emerges due to competitive binding of the stalk with the Ig-like domain. Additionally, we observed that the stalk domain itself functions as a site for ligand binding, with increased binding frequency in the presence of R47H. This suggests that sTREM2’s neuroprotective role in AD may, at least in part, arise from the stalk domain’s ability to rescue dysfunctional ligand binding caused by AD-risk mutations. Lastly, our findings indicate that R47H-induced dysfunction in TREM2 may result from both diminished ligand binding due to restricted complementarity-determining region 2 loop motions and an impaired ability to differentiate between ligands, proposing a novel mechanism for loss-of-function. In summary, these results provide valuable insights into the role of sTREM2 in AD pathology, laying the groundwork for the design of new therapeutic approaches targeting (s)TREM2 in AD.
Introduction
Alzheimer’s disease (AD) is a fatal neurodegenerative condition marked by progressive memory loss, cognitive decline, and impaired daily functioning1. Despite extensive research, the root causes remain unknown and curative treatments remain elusive. AD is hallmarked by the presence of extracellular amyloid-β (Aβ) plaques, intracellular neurofibrillary tau tangles, lipid droplet accumulation, and neuroinflammation2,3. Microglia, the brain’s macrophages, play a pivotal yet complex role in AD pathology. While their phagocytic activity toward Aβ is considered neuroprotective, they also release cytokines that increase neuroinflammation, ultimately harming neighboring neurons and glia4,5. Two newly approved monoclonal antibodies targeting Aβ, while transformative, have shown limited efficacy, only modestly slowing disease progression6. This underscores the need for a deeper understanding of the intricate interactions between the neuroimmune system and AD-relevant proteins.
At the forefront of this investigation is Triggering Receptor Expressed on Myeloid Cells 2 (TREM2), a transmembrane receptor expressed on microglial surfaces. TREM2 propagates a downstream signal through interactions with its co-signaling partner, DNAX-activating protein 12 (DAP12). Recent research highlights TREM2’s crucial role in modulating microglial responses and maintaining brain homeostasis7. TREM2 contains an extracellular immunoglobulin (Ig)-like domain (19-130 amino acid (aa)), a short extracellular stalk domain (131-174 aa), a helical transmembrane domain (175-195 aa), and an intracellular cytosolic tail domain (196-230 aa)8. N-linked glycosylation of TREM2 can occur at residues 20 and 79, however its effects are debated9,10. The Ig-like domain contains several ligand binding sites. These include a hydrophobic tip, which is characterized by a highly positive electrostatic potential and contains several aromatic residues and three complementarity-determining regions (CDRs), the latter of which (CDR2) also span part of a positively-charged basic patch known as the putative ligand interacting region, or ‘Surface 1’ (Fig. 1)8,11–13. Mutations in TREM2 correlate with altered risks of developing AD, with the R47H mutation, found on Surface 1, standing out as a significant genetic risk factor14–16. In vitro and in silico studies have consistently revealed that mutations, including R47H, destabilize TREM2’s CDRs8,17–19, exposing once-buried negatively charged residues8,18 and disrupting homeostatic TREM2-ligand binding behavior11,12,18,20.

Overview of structural domains in sTREM2 and full-length TREM2
(A) Pre-MD structure of sTREM2WT. (B) Full sequence of TREM2, indicating significant structural domains.
TREM2 binds diverse anionic and lipidic ligands, including Aβ species12,21,22, lipoproteins12,23–25, nucleic acids26, carbohydrates27, and phospholipids (PLs)11,20,28—the focus of our study. PLs play a crucial role in lipid metabolism and maintaining brain homeostasis29. TREM2 clears excess PLs during demyelination and interacts with PLs when they are bound to lipoproteins with a core of neutral lipids surrounded by a monolayer of PLs, free cholesterol, and apolipoproteins30,31. Many PLs bind to TREM2, including phosphatidyl-choline (PC), -serine (PS), -inositol (PI), - glycerol (PG), -ethanolamine (PE), phosphatidic acid (PA), sphingomyelin (SM), cholesterol, and sulfatide11,18,20,31. PLs primarily bind to TREM2’s hydrophobic tip and Surface 1, with varying affinities observed in different contexts11,18,20. Direct binding assays show stronger TREM2 binding to anionic moieties (PS, PE, PA) and weaker binding to PC and SM. Conversely, TREM2-expressing reporter cells reveal high TREM2 stimulation from PC and SM, especially with the TREM2R47H variant20. Collectively, however, these results emphasize TREM2’s broad binding capabilities for PLs. Yet, one review suggested effective TREM2 stimulation by PLs may require co-presentation with other molecules, potentially reflecting the nature of lipoprotein endocytosis32. Another study observed minimal changes in TREM2-PL interactions despite TREM2 mutations (R47H, R62H, T96K)11. Ultimately, there remains a major gap in understanding the mechanism by which PLs differentially interact with TREM2 and how these interactions are altered in the presence of disease-associated TREM2 variants.
The activation of TREM2, mediated by the binding of ligands such as PLs, shapes key microglial functions, including proliferation, phagocytosis, and lipid metabolism. Notably, in pathological conditions, ligand-induced TREM2 activation triggers microglial phenotype switching to Disease-Associated Microglia (DAM), characterized by the activation of inflammatory, phagocytic, and lipid metabolic pathways33. Additionally, TREM2’s extracellular domain can undergo cleavage from ADAM 10/17 sheddase at residue H157, yielding soluble TREM2 (sTREM2). The role and relevance of sTREM2 in disease pathology has been heavily debated34. In the cerebrospinal fluid of individuals with early-stage AD, elevated sTREM2 levels have been detected and linked to slower AD progression35–37. Further, there is a strong correlation between sTREM2 levels in cerebrospinal fluid and those of Tau, although correlation with Aβ is inconclusive. These findings have established sTREM2 as a long-time biomarker for AD diagnosis and progression38–40.
Many studies have indicated a neuroprotective role for sTREM2 in disease pathology41. It has been suggested, for instance, that sTREM2 may function as a “dummy receptor” in AD states, preventing disease-associated ligands from binding TREM234. Moreover, in vivo AD mouse models evaluating the therapeutic potential of recombinant sTREM2 have observed the suppression of microglial apoptosis, reduced Aβ plaque load, and improved learning and memory abilities42,43. More recent studies have indicated that sTREM2 not only serves as an activator for microglial uptake of Aβ but also directly inhibits Aβ aggregation12,13,44. Specifically, the binding of Aβ to TREM2 has been shown to increase shedding of sTREM244, which can then bind to Aβ oligomers and fibrils to inhibit their secondary nucleation13,44. Interestingly, the effect of R47H on Aβ aggregation is unclear, highlighting the need to study mechanistic aspects of ligand binding13,44.
Some anionic ligands, including Aβ, predominantly bind to Surface 1 on TREM2. Intriguingly, recent observations revealed that Aβ binds to an alternative binding region, termed ‘Surface 2’, on sTREM2, situated opposite Surface 1 (Fig. 1). Surface 2 features a group of positively charged residues surrounded by acidic residues, creating a variegated electrostatic potential13. Herein, we aimed to unravel the molecular basis behind this functionally significant distinction in ligand binding between soluble and membrane-bound TREM2, utilizing molecular dynamics (MD) simulations. We focused on (s)TREM2-PL interactions, establishing a controlled framework to assess the impacts of various PL chemistries on ligand binding, specifically comparing the binding behavior of anionic PS and neutral PC. We hypothesized that the oft-overlooked flexible stalk domain of sTREM2, minimally explored in previous in silico studies, may play a pivotal role in mediating the observed variations in binding. Furthermore, we sought to understand the impact of the AD-risk mutation R47H on ligand binding, thereby unraveling fundamental roles of (s)TREM2 in AD pathology. To our knowledge, this study represents the second-ever application of MD to investigate sTREM2, totaling an unprecedented 31.2 µs of simulation time. Ultimately, this research may unveil new insights into the mechanistic and therapeutic roles of sTREM2 in AD.
Methods
Preparation of simulated structures
We employed the AlphaFold45,46 model (AF-Q9NZC2-F1) as the initial structure for wildtype (WT) sTREM2 in our simulations, chosen for its inclusion of the unstructured flexible stalk domain. The partial stalk domain spans residues 130 through 157, while the Ig-like domain consists of residues 19 through 130. For WT simulations, we constructed two protein systems: one with just the Ig-like domain (“IgWT”) and another containing both the partial stalk and Ig-like domain (“sTREM2WT”). Similarly, two protein systems were constructed for the variant simulations: one with just the Ig-like domain containing the R47H mutation (“IgR47H”) and another containing both the partial stalk and mutant Ig-like domain (“sTREM2R47H”). In contrast to the use of the AlphaFold model for WT, we utilized a crystal structure of TREM2R47H (Protein Data Bank (PDB) code 5UD818), to which the unstructured stalk domain from the AlphaFold model was added using alignment tools in Visual Molecular Dynamics (VMD). Missing residues were incorporated into the TREM2R47H Ig domain using MODELLER47. The initial molecular structures for the PLs, stearoyl-oleoyl-PC (SOPC) and stearoyl-oleoyl-PS (SOPS), were obtained from CHARMM-GUI48–52, with each considered as a singular PL. All eight protein-ligand systems (IgWT, sTREM2WT, IgR47H, and sTREM2R47H, each with SOPC or SOPS) and six pure-component systems were solvated with explicit water and with counterions added as needed to neutralize the charge.
Molecular dynamics simulations
All proteins, PLs, and counterions were parameterized using the CHARMM36 force field53,54, while the TIP3P model was used to describe water55. Prior to subsequent docking studies, each pure-component system underwent steepest-descent energy minimization, initially in vacuum and then in a solvated state. This was followed by a multi-step equilibration protocol, which included a 1 ns NVT equilibration simulation at 310K using the Bussi-Donadio-Parrinello thermostat55, followed by a 1 ns NPT equilibration simulation at 310K and 1 bar using the same thermostat and Berendsen barostat56. Finally, production simulations were carried out at 310 K and 1 bar, utilizing the same thermostat and Parrinello-Rahman barostat57. The duration of the production simulations was 150 ns for each pure-component PL system and 1 µs for each pure-component protein system. Each pure protein system was run with six replicates (see Results).
All simulations were conducted using the GROMACS MD engine58. The LINCS algorithm59 was used to constrain bonds involving hydrogen atoms, and particle-mesh Ewald (PME) summations60 were employed for calculating long-range electrostatics with a cutoff of 1.2 nm. Lennard-Jones interactions were evaluated up to 1.2 nm and shifted to eliminate energy discontinuities. Neighbor lists were reconstructed in 10-step intervals with a 1.4 nm cutoff. A timestep of 2 fs was implemented in all simulations, and periodic boundary conditions were applied in the x, y, and z directions. Configurations from these production simulations were used as inputs in ensuing docking calculations (see next section). After the docking calculations, we conducted additional 150 ns production simulations on the combined post-docking models, employing the same parameters as described above for the pre-docking production simulations.
Molecular docking calculations
For each of the four pre-docking, pure-component protein systems, we clustered our initial 1 µs production simulation trajectory using the gromos method implemented in GROMACS. Representative structures from the top two clusters in each case were selected and prepared for subsequent docking calculations using AutoDock Tools61. Docking calculations were carried out with AutoDock Vina62,63, treating the proteins as rigid receptors. Given that AutoDock Vina employs a flexible ligand docking procedure, the final PL conformation from each pure-component simulation served as the ligand in the docking calculations. Grids with dimensions of 30 Å x 30 Å x 30 Å were constructed, redundantly covering the complete surface of each receptor protein. The exhaustiveness parameter for docking was set to eight. Docked complexes were analyzed based on the AutoDock score, the number of highly similar complexes, and biological relevance. Unique structures across grids and clusters for each ligand-receptor system were identified for post-docking MD simulations. From molecular docking, we obtained 7 unique SOPS/IgWT models, 5 SOPS/sTREM2WT models, 5 SOPC/IgWT models, 6 SOPC/sTREM2WT models, 4 SOPS/IgR47H models, 8 SOPS/sTREM2R47H models, 6 SOPC/IgR47H models, and 7 SOPC/sTREM2R47H models.
Trajectory analysis protocols
The Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) approach was used to calculate PL-protein binding free energies, utilizing the gmx_MMPBSA implementation64–66. The analysis focused on temporal regions of the simulation trajectories where PL-protein complexes demonstrated stability, as indicated by relatively constant root-mean-square deviation (RMSD) values for at least 50 ns. Over the same temporal regions, we calculated Interaction Entropy, as implemented in gmx_MMPBSA, to estimate entropic trends in PL binding65,67. Significance testing for these values was performed using two-tailed heteroscedastic T-tests. Fractional occupancy values, characterizing the frequency and location of PL binding on each protein surface, were determined by calculating the fraction of simulation trajectory frames in which a given protein residue was within 4 Å of the PL across all simulations. We calculated conformational changes of the CDR2 loop using MDAnalysis by measuring the distance between residues 45 and 70 at the tops of the CDR1 and CDR2 loops, respectively.
Results and Discussion
The partial stalk domain of sTREM2 differentially modulates the CDR2 loop and broader Ig-like domain dynamics via ‘Surface 1’ binding in WT and R47H models
Prior TREM2 research identified an open CDR2 loop in R47H models, which disrupts ligand binding to Surface 18,17,18,68. To validate our in silico approach, we first sought to recapitulate these findings by performing a 1-μs simulation each of IgR47H and sTREM2R47H. Consistent with prior studies, these “initial” simulations showed persistently open CDR2 loops in both constructs (Fig. S1A, slate blue and teal lines, respectively; see Methods for loop characterization details). Surprisingly, however—and in contrast to prior experimental studies17,18—a 1-μs simulation each of IgWT and sTREM2WT (Fig. S1A, pink and red lines, respectively) revealed that the former transitioned from a closed to an open CDR2 loop midway through the simulation.
To understand this unexpected result, we examined whether variation in the protonation state of histidine residue H43 adjacent to the CDR2 loop of IgWT and IgR47H could account for the observed differences in CDR2 loop dynamics. Initial protonation states (HSE in IgWT, HDE in IgR47H) were assigned by the pdb2gmx algorithm in GROMACS, based on optimal hydrogen-bonding conformations. Two additional 1-μs simulations of IgWT were performed with the original protonation state (HSE) and three with the alternate state (HDE), totaling six replicates of IgWT. The open CDR2 loop conformation was sampled in two of the six simulations (Fig. S2A-i, purple and brown lines): the original IgWT simulation with HSE, and one replicate with HDE. We note that the CDR2 loop remained open in a 1-μs simulation of IgR47H with the alternate protonation state (HSE), as well as in four additional replicates performed with the original protonation state (HDE), resulting in six replicates of IgR47H that sampled an open CDR2 loop (Fig. S2B-i). These results indicate that the protonation state of H43 does not explain the differences in loop dynamics. Rather, the R47H mutation appears to stabilize the open CDR2 loop conformation, whereas loop opening in the WT construct may be more stochastic or transient. These conclusions are further supported by Fig. 2A, which shows the average CDR2 loop distance and standard error of the mean (SEM) across all six replicates of both IgWT and IgR47H. This analysis reveals that the average loop distance increased in the second half of the IgWT simulations, albeit to a lesser extent than in the R47H models.

The stalk of sTREM2WT modulates the Ig-like domain via binding to Surface 1
(A) Distance between CDR loop residues 45 and 70 over time for WT and R47H models of sTREM2 and TREM2 (Ig), averaged across six replicates. Error bars represent the SEM. (B) Temporally averaged Cα RMSF of the isolated Ig-like domain of WT and R47H sTREM2 compared to the Ig-like domain of TREM2, averaged across six replicates. Error bars represent the SEM. (C-D) Cα RMSD of (C) WT and (D) R47H models of sTREM2 and TREM2 over time during the initial simulations, with corresponding snapshots shown for the starting and equilibrated structures in each case. The stalk domain of sTREM2 is shown in dark purple, and the CDR2 loop of TREM2 is shown in light purple. (E-F) Normalized fractional occupancy of residues in the Ig-like domain of (E) sTREM2WT and (F) sTREM2R47H by the stalk during the initial simulations.
Notably, sTREM2WT maintained a closed-loop conformation throughout the initial 1-μs simulation (Fig. S1A). As with the Ig constructs, we again performed five additional 1-μs simulations of sTREM2WT (and of sTREM2R47H) to further probe this behavior and more robustly sample protein dynamics. Similar to IgWT, an open CDR2 loop conformation was sampled in two of the six sTREM2WT simulations (Fig. S2C-i, yellow and brown lines). However, averaging the loop distance across replicates suggests that CDR2 loop opening occurs with a delay of approximately 200 ns in sTREM2WT compared to IgWT (Fig. 2A). This observation hints at a stronger thermodynamic preference for the closed-loop conformation in sTREM2WT, potentially due to the presence of the partial stalk domain. In sTREM2WT, the stalk may help stabilize the closed state through intramolecular interactions that are sterically hindered in membrane-bound TREM2 (represented by IgWT), where it is tethered to the transmembrane domain and less conformationally flexible—a hypothesis supported by additional structural analyses, described below.
Given the potential stochastic nature of CDR2 loop opening in the absence of the R47H mutation, we sought to examine its impact on protein stability more rigorously using root-mean-square fluctuation (RMSF) analysis. These calculations included the six IgWT simulations (Fig. S2A-ii) and the isolated Ig-like domain from the six sTREM2WT simulations (sTREM2WT-Ig; Fig. S2C-ii), which were subsequently averaged across trajectories for each system (Fig. 2B). RMSF was calculated relative to the post-minimization structure, which adopted a closed CDR2 loop in WT structures and an open CDR2 loop in R47H structures. On average, sTREM2WT-Ig exhibited slightly reduced, but statistically overlapping, residue-level fluctuations across most of the Ig-like domain compared to IgWT. This suggests that subtle differences in sTREM2WT dynamics may stem from stalk-mediated effects on the sampling of specific conformational states. In light of these findings, we conducted a similar RMSF analysis to compare the dynamics of the six IgR47H simulations (Fig. S2B-ii) with those of the isolated Ig-like domain from the six sTREM2R47H simulations (sTREM2R47H-Ig; Fig. S2D-ii). In contrast to the WT model, the presence of the stalk in sTREM2R47H-Ig led to substantially increased residue-level fluctuations across much of the Ig-like domain, particularly in residues 70-80, which comprise most of the CDR2 loop (Fig. 2B). Collectively, these findings suggest a dual role for the partial stalk in modulating Ig-like domain dynamics: reducing the likelihood of stochastic CDR2 loop opening in WT constructs, while destabilizing the open loop conformation in the R47H mutant. This stalk-mediated destabilization may contribute to sTREM2’s neuroprotective function by preventing prolonged stabilization of the open CDR2 loop state, which is promoted by R47H and may underlie its pathogenic effects.
We next aimed to elucidate the molecular basis of interactions between the Ig-like domain and partial stalk domain in sTREM2. To this end, we calculated the RMSD of Cα atoms across the six replicate simulations for each construct (Figs. S2A-D-iii), tracking temporal deviations relative to the post-equilibration configuration in each case. To highlight specific CDR2 loop dynamics of interest, we compared the RMSD of the initial 1-μs simulations of IgWT and sTREM2WT, and then for IgR47H and sTREM2R47H, as shown in Figs. 2C and 2D, respectively. Stalk convergence analyses across all sTREM2 replicates confirmed that these initial simulations captured a dominant high-stability conformation for each construct (Fig. S1B; WT1 (left) and R47H1 (right)). The RMSD of each structure converged relatively quickly during the simulations. IgWT and IgR47H displayed behavior consistent with the CDR2 loop distances depicted in Fig. 2A. IgR47H exhibited a consistently low RMSD profile throughout the simulation, indicative of a persistently open CDR2 loop (Fig. 2D, inset iii). In contrast, the RMSD profile for IgWT remained low until approximately 400 ns, at which point it began to increase, signaling a transition from a closed to an open CDR2 loop (Fig. 2C, insets iii-iv). Although the loop remained open for the rest of the simulation, IgWT showed markedly greater RMSD fluctuations than IgR47H, suggesting that the R47H mutation likely imparts some stability to the open-loop conformation.
Sharp initial increases in RMSD are observed with both sTREM2WT and sTREM2R47H. In both cases, the starting structures (inset panels i of Figs. 2C-D) exhibited a generally linear configuration of the stalk, potentially resembling its membrane-bound form. After the initial RMSD increases, the partial stalk domains were observed to consistently interact with the Ig-like domains of sTREM2WT and sTREM2R47H (inset panels ii of Figs. 2C-D). To identify the specific regions of the protein surface with which the partial stalk interacted during the representative simulations, we generated per-residue fractional occupancy maps characterizing the contact frequency (atom-atom distance within 4 Å) between the stalk and each residue in sTREM2WT (Fig. 2E) and sTREM2R47H (Fig. 2F). We observed a notable increase in the interaction frequency between residues in Surface 1 and the partial stalk domain of sTREM2R47H compared to sTREM2WT. This difference is likely due to the more persistently open CDR2 loop of sTREM2R47H, providing greater accessibility of Surface 1 residues for interaction with the stalk. We generated similar occupancy maps, averaged across replicate simulations and observed similar, albeit more diffuse, patterns (Figs. S1C-D).
Given the highly negative overall charge of the stalk domain (−8), its recurrent interactions with the positively charged residues in Surface 1 are perhaps unsurprising. In addition to these contacts, however, the stalk domain in sTREM2WT also frequently interacted with residues 20-35, though typically with lower occupancy. While this region has not been directly noted in previous experimental studies, it lies directly adjacent to residues comprising Surface 2 as defined by Belsare et al13. We therefore refer to it here as ‘Expanded Surface 2’. Taken together, the diversity of surface regions contacted by the stalk domain—and the overlap in many cases with known TREM2 ligand-binding sites—suggests a potential mechanism for differences in ligand binding between sTREM2 and full-length TREM2, which we explore in detail in the following section.
sTREM2’s partial stalk domain modulates PL binding by promoting interactions with a new site on the Ig domain and creating an additional binding site within the stalk itself
In vitro studies examining molecular interactions between (s)TREM2 and ligands remain limited. One study reported that sTREM2 exhibited slightly lower, but statistically insignificant, affinity for monomeric Aβ1-42 compared to TREM212. Based on the similar binding affinities, the authors proposed that the Ig-like domain constitutes the principal binding surface of sTREM2. However, another recent report from the same group found that the TREM2 Ig-like domain does not engage ligands at ‘Surface 2’69, contrasting earlier cross-linked mass spectrometry data showing Aβ fibrils binding to this region in sTREM270. These conflicting results suggest that differences in ligand binding between sTREM2 and TREM2 likely stem from more complex factors than simply affinity and ligand type. Another study presented a crystal structure of a TREM2 trimer bound by three PS molecules, revealing extensive PS interactions with residues in the CDR2 loop and hydrophobic tip18. While the study did not examine PL/sTREM2 interactions, we hypothesized that similar differential patterns of (s)TREM2 binding may occur with PLs—as previously seen with Aβ—due to stalk-induced alterations in ligand accessibility, which may in turn modulate binding affinity.
To test our hypothesis, we conducted 150-ns simulations of both IgWT and sTREM2WT bound by the PLs stearoyl-oleoyl-PC (SOPC) and stearoyl-oleoyl-PS (SOPS) in various favorable starting docking configurations (see Methods). As discussed previously, the initial 1-μs simulation of each protein construct was deemed representative of the replicate-averaged behavior and was thus used for these docking studies. We then analyzed PL binding by calculating the fractional occupancy of each protein residue, temporally averaged across all simulations for each protein/PL system. Across all four systems, the highest occupancy by both PLs was observed at residues within the CDR2 loop and hydrophobic tip (Figs. 3A-D, primarily green regions), consistent with earlier in vitro findings and supporting the validity of our approach. Nevertheless, distinctions between the SOPS and SOPC models were also evident, underscoring the influence of ligand charge on binding interactions. Specifically, SOPS, being negatively charged, showed higher relative occupancy of the positively charged residues on Surface 1 compared to SOPC, which is neutrally charged (Figs. 3A-B vs. 3C-D, respectively, blue regions). Conversely, SOPC showed higher relative occupancy of residues within the hydrophobic tip than SOPS (Figs. 3C-D vs. 3A-B, green regions).

The stalk of sTREM2WT modulates phospholipid binding
(A-D) Normalized fractional occupancy of residues in (A) IgWT by SOPS, (B) sTREM2WT by SOPS, (C) IgWT by SOPC, and (D) sTREM2WT by SOPC, averaged across seven, five, five, and six 150-ns simulations, respectively. (E-H) VMD snapshots of representative structures from the PL/(s)TREM2WT simulations, with PLs and corresponding complex binding free energies shown in matching colors.
Upon comparing the fractional occupancy plots of the PL/sTREM2WT simulations (Figs. 3B,D) with those of the PL/IgWT simulations (Figs. 3A,C), we observed a noticeable decrease in the relative occupancy of residues within Surface 1, the hydrophobic tip, and particularly the CDR2 loop in the sTREM2WT systems. Instead, both PL/sTREM2WT simulations showed ligand binding at the newly defined Expanded Surface 2 (Figs. 3B,D, light orange regions), which was previously shown to interact with the sTREM2WT stalk, albeit with lower occupancy than Surface 1 residues (Fig. 2E). In contrast, PL binding to Expanded Surface 2 was almost entirely absent in the IgWT simulations. Collectively, these results suggest that Expanded Surface 2 functions as a secondary ligand-binding surface in sTREM2, becoming preferentially engaged when Surface 1 is occupied by the flexible stalk (and when Expanded Surface 2 is not itself engaged by the stalk). Finally, we observed frequent interactions between both PLs and residues within the stalk domain of sTREM2WT (Figs. 3B,D, pink regions), suggesting that the stalk itself may independently serve as an additional binding site for diverse ligands.
We gained deeper insights into PL/(s)TREM2 interactions through visual analysis of our simulation trajectories using VMD and binding free energy calculations performed with the MM-PBSA approach (see Methods). Representative snapshots from converged portions of the simulations for each protein/PL system (defined by stable RMSD values for at least 50 ns) are shown in Figs. 3E-H, alongside their corresponding binding free energies. Among the WT systems, a SOPS/IgWT complex exhibited the most favorable (lowest) binding free energy, consistent with TREM2’s known affinity for anionic ligands. Notably, in the lowest-energy binding models, both SOPS and SOPC directly engaged the CDR2 loop in IgWT (Fig. 3E, gray model; Fig. 3G, orange model, respectively). In contrast, both SOPS and SOPC bound most favorably to the newly defined Expanded Surface 2 on sTREM2WT (Fig. 3F, orange model; Fig. 3H, purple model). RMSF calculations revealed heightened dynamics in the stalk domain of sTREM2WT for these particular models (Figs. S3B-iii and S3D-iii), suggesting that PL binding perturbs homeostatic stalk/Ig-like domain interactions. This is further supported by the observation that, in the absence of a ligand, the stalk most frequently interacts with Surface 1 (Fig. 2E)—the same region where PLs, especially SOPS, bind in the absence of the stalk (Figs. 3A,C), with reduced occupancy when the stalk is present (Figs. 3B,D). Together, these results indicate that the PLs compete with the stalk domain for access to the CDR2 loop and Surface 1, likely due to their shared amphipathic and negatively charged character.
To gain deeper insights into the thermodynamics of PL/(s)TREM2 interactions, we performed Interaction Entropy (IE) calculations (see Methods). While many free energy estimates neglect entropy, recent studies have highlighted its essential role in characterizing the full Gibbs free energy landscape in biomolecular systems71,72. IE calculations capture entropic changes in ligand, protein, and solvent interactions, and efficiently estimate these contributions directly from MD simulations, making them particularly useful for comparing relative entropies67. Our results broadly indicate that the entropic contributions (-TΔS) were positive across all systems, with higher values for SOPS models than SOPC models, although the differences were not statistically significant (Fig. S4I). Entropic loss upon ligand binding reflects reduced conformational freedom in the ligand, protein, and/or surrounding solvent67,71. Thus, the larger entropic penalties observed for SOPS models suggest greater ‘snugness’ of binding73, implying that binding in these systems is more enthalpically driven. However, interpretation of these results is limited by large standard deviations, likely due to the highly dynamic nature of PL/(s)TREM2 interactions74.
Finally, some PL/TREM2 complexes were observed to undergo rapid conformational changes in the CDR2 loop during the simulations, indicated by marked increases in the RMSD of the protein (Figs. S3A-i,C-i) combined with high RMSF values in CDR2 loop residues (Figs. S3A-iii,C-iii). These changes—likely facilitated by elevated, direct PL binding to the CDR2 loop of IgWT compared to sTREM2WT, as discussed earlier—occurred on much shorter timescales than in the ligand-free simulations and included dynamic opening and closing of the loop along with shifts in its α -helical character. This suggests a mechanism by which the CDR2 loop may function as a dynamically responsive element during ligand binding. However, further experimental studies are needed to clarify the functional significance of CDR2 loop remodeling in response to ligand binding and to define the conditions under which such conformational transitions occur.
The AD-risk mutation R47H diminishes ligand discrimination by and binding to the Ig-like domain of (s)TREM2, while sTREM2’s stalk domain partially restores overall ligand binding
To investigate the roles of (s)TREM2 in disease pathology, we examined PL binding in the presence of the AD-risk mutation R47H. Previous studies have suggested that the decreased ligand-binding capabilities of TREM2R47H may underlie its observed loss-of-function11,12,18,20. Notably, one study showed reduced reporter cell activity of TREM2R47H compared to TREM2WT when binding anionic lipids such as PS, while no significant difference was observed in their binding to PC20. Separately, sTREM2R47H and sTREM2WT were shown to bind Aβ and inhibit fibrillization to a similar extent13. Taken together, these studies suggest that the effects of R47H extend beyond reduced binding affinity and may also be influenced by ligand type, binding site, membrane context, and the presence of the stalk.
Upon examining the fractional occupancy plots of the PL/(s)TREM2R47H simulations, we observed broadly similar patterns of binding to the Ig-like domain, regardless of stalk presence or PL charge (Figs. 4A-D). Across all four simulations, PLs predominantly occupied the CDR2 loop and Surface 1 of the Ig-like domain. This uniformity in binding contrasts sharply with the previously discussed behavior in the PL/(s)TREM2WT simulations (Figs. 3A-D), where distinct differences were observed between SOPS and SOPC binding to both IgWT and sTREM2WT, as well as between the two proteins for a given PL. These results suggest that the presence of R47H diminishes the ability of IgR47H and sTREM2R47H to distinguish between crucial brain-derived ligands, proposing a novel mechanism of loss-of-function for this key AD-risk mutation.

The R47H mutation decreases ligand discrimination but increases binding to the flexible stalk domain
(A-D) Normalized fractional occupancy of residues in (A) IgR47H by SOPS, (B) sTREM2R47H by SOPS, (C) IgR47H by SOPC, and (D) sTREM2R47H by SOPC, averaged across four, seven, six, and eight 150-ns simulations, respectively. (E-H) VMD snapshots of representative structures from the PL/(s)TREM2R47H simulations, with PLs and corresponding complex binding free energies shown in matching colors. (I) Comparison of the binding free energies averaged across all WT and R47H models for each PL/protein system. Errors bars represent the standard error of the mean.
Compared to the PL/(s)TREM2WT simulations, all four PL/(s)TREM2R47H simulations demonstrated markedly reduced PL binding to the hydrophobic tip. These reductions likely stem from conformational effects induced by R47H, which promotes a more open and rigid CDR2 loop adjacent to the hydrophobic tip, thereby altering its accessibility and ligand-binding properties. In addition, changes in stalk-Ig-like domain interactions observed in the presence of R47H (Figs. 2E-F) may further influence PL binding by reshaping the local structural landscape. Consistent with these effects, we also observed diminished PL binding to Expanded Surface 2 in the sTREM2R47H simulations, suggesting that R47H not only impairs canonical ligand recognition surfaces but may also disrupt secondary binding modes unique to sTREM2.
Compensating for these reductions in binding, increased PL binding to Surface 1 was observed in all four R47H simulations, along with enhanced PL binding to the flexible stalk of sTREM2R47H. These results provide further support for the notion that sTREM2’s stalk domain functions as an independent ligand-binding site. Moreover, they suggest the stalk domain may possess an inherent capacity to ‘rescue’ sTREM2 deficiencies in ligand binding to the Ig-like domain caused by disease-associated mutations, offering a mechanistic explanation for sTREM2’s neuroprotective role in AD. Indeed, this may help explain why, as noted earlier, R47H does not significantly impair sTREM2R47H’s ability to bind Aβ and inhibit fibrillization to the same extent as sTREM2WT.
Following the methodologies used to evaluate the PL/(s)TREM2WT simulations, we conducted visual analysis of our trajectories paired with binding free energy calculations. Representative structures with corresponding binding free energies are shown in Figs. 4E-H. Broadly, the trends align closely with those observed for the WT models (Figs. 3E-H), with a SOPS/IgR47H complex again exhibiting the lowest (most favorable) binding free energy (Fig. 4E). Unlike the WT models, however, these results do not indicate competitive binding between the PLs and the stalk, as evidenced by high PL interaction frequencies with the CDR2 loop and Surface 1 that persist even when the stalk is present (Figs. 4B,D vs. 4A,C), despite the stalk’s increased interaction frequency with Surface 1 in the context of R47H (Fig. 2F vs. 2E). This is likely due to the more persistently open CDR2 loop in the R47H complexes, which alters surface accessibility and stalk dynamics.
Our binding free energy results show that across all models, sTREM2 and TREM2 bind PS more favorably than PC (Fig. 4I), affirming previous findings that TREM2 favors anionic ligands32. However, these differences were not statistically significant. We also observed no significant difference or trend in binding free energies between PS/(s)TREM2WT and PS/(s)TREM2R47H complexes, contrasting with previous experimental findings that R47H reduces TREM2’s affinity for endogenous ligands12,18,20. Similar to WT models, SOPS binding poses in the R47H models exhibited higher average entropic contributions than their SOPC counterparts, with statistically significant differences for sTREM2R47H but not for IgR47H (Fig. S4I). Together, these observations suggest that reduced ligand affinity alone does not fully account for the signaling deficits associated with the R47H variant. More broadly, our results indicate that differences in binding across a range of ligands, (s)TREM2 variants, and between sTREM2 and TREM2 are likely driven by differences in accessibility of protein surface residues and ligand binding patterns, factors that extend beyond mere differences in affinity.
Lastly, we observed a prominent difference in the dynamic behavior of the CDR2 loop between IgR47H and IgWT. In most simulations, the CDR2 loop in IgR47H remained stably open regardless of PL binding (Fig. S5), whereas in IgWT, the loop underwent dynamic opening and closing in response to ligand engagement. These findings suggest that, beyond impairing ligand affinity and selectivity, the R47H mutation may also hinder the CDR2 loop’s ability to respond adaptively to ligand binding, potentially disrupting physiological TREM2 signaling. Future investigation of the TREM2-DAP12 complex will be essential to clarify how ligand binding to the Ig-like domain propagates downstream signaling, and how this process may be disrupted by the R47H variant.
Discussion and Conclusions
Main Conclusions
We utilized long-timescale MD simulations to investigate relevant structural domains of sTREM2 and TREM2, along with their interactions with key PLs in the brain. Through the analysis of RMSF, RMSD, and PL/protein residue occupancy calculations, we established the flexible stalk domain of sTREM2 as the molecular basis for differential ligand binding between sTREM2 and TREM2. This difference in binding arises from the stalk’s role in modulating dynamics of the Ig-like domain and altering the ligand accessibility of its surface residues. By integrating free energy, interaction entropy, and ligand occupancy calculations, we quantified the energetics of these interactions, confirming the presence of an alternate ligand binding site on sTREM2WT, which we termed the ‘Expanded Surface 2’. Binding of PLs to this site arises from competitive binding of the flexible stalk to Surface 1 on the Ig-like domain. These stalk-Ig-like domain interactions were disrupted in the presence of the AD-risk mutation R47H and entirely absent in the Ig (TREM2) models, indicating occupancy of Expanded Surface 2 occurs solely with sTREM2WT. These observations underscore our conclusion that, rather than (or in addition to) sTREM2’s previously conceived role as a dummy receptor for TREM2, the flexible stalk confers sTREM2 with unique ligand binding preferences and patterns that facilitate distinct endogenous functions compared to TREM2.
Furthermore, we found that the stalk domain itself serves as an independent site for ligand binding, with heightened PL occupancy observed in the presence of R47H compared to the wildtype model. This suggests that the stalk domain may have the capacity to partially ‘rescue’ dysfunctional ligand binding to the Ig-like domain of sTREM2 caused by disease-associated mutations like R47H. Moreover, our observations for both sTREM2 and TREM2 indicate that R47H-induced dysfunction may result not only from diminished ligand binding but also an impaired ability to discriminate between different ligands in the brain, proposing a novel mechanism for loss-of-function. In summary, the findings of this study reveal the endogenous structural and dynamical mechanisms of (s)TREM2, a critical component in AD pathology. These insights offer new fundamental knowledge that can serve as guiding principles for the design of future therapeutics, paving the way for potential advancements in AD treatment strategies.
Ideas and Speculation
Our findings suggest that loss-of-function in sTREM2 and TREM2 occurs not only through reductions in ligand binding affinity, but also through a change in ligand binding patterns and a loss of ligand discrimination capacity. Altered ligand binding may hinder the ability of TREM2 and its co-signaling partner DAP12 to transmit signals across the cell membrane. This deficiency implicates an impaired microglial response to ligand binding as a key mechanism for dysfunction in AD. Pathologically, this would lead to reduced lipid uptake75, diminished Aβ plaque clearance21, decreased microglial activation75, and inhibited intracellular lipid metabolism76. Previous studies speculate that reduced microglial lipid metabolism can trigger neurotoxic activation states or loss of neuroprotective functions77. Furthermore, the inability to discriminate among PLs may result in an invariable response from microglia when presented with diverse ligands. While PC is the most abundant PL in cell membranes, PS expression on the outer membrane leaflet increases in apoptotic cells, acting as a damage-associated signal78. Failure to differentiate between these PLs may lead to chronic over-activation of microglia, or at the very least, loss of endogenous protective functions. Additionally, there are two known alternatively spliced isoforms of sTREM2 which constitute the minority of CSF sTREM2 levels. These isoforms occur due to alternate splice sites and skipping within exon 470. While our study suggests that the flexible stalk domain of sTREM2 is highly relevant, the specific conclusions do not extend to these isoforms, as they have entirely different stalk compositions and lengths. A separate study of these isoforms is pertinent.
Recently, monoclonal antibodies targeting TREM2 have emerged79–81, designed to bind to the stalk region above its cleavage site80,81. Consequently, the primary binding epitope of at least one of these antibodies also resides on sTREM2. Notably, this antibody also reduces the shedding of sTREM281. It is conceivable that treatment with this and similar antibodies may compromise the endogenous function of sTREM2, given that the stalk, a domain our study has identified as functionally significant, may be sequestered by the antibody. This may ultimately result in off-target effects for these therapeutics and prompts a consideration of how to separately target sTREM2 and TREM2.
Limitations
As with all models, particularly computational ones, it is crucial to recognize their limitations. In our study, the initial positions of PLs bound to (s)TREM2 were determined exclusively by AutoDock Vina. To mitigate potential bias from the docking process, we conducted MD simulations on all unique structures observed in the top 20 docked models. Rigorous, yet computationally efficient free energy calculations remain a challenge for computational protein-ligand studies. Despite employing the latest MM-PBSA free energy calculation algorithms and novel interaction entropy calculations, achieving statistical significance for many of the trends we observed that align with experiments remained challenging. This difficulty was further compounded by our approach of averaging statistics across multiple models of PLs docked to different initial binding locations on (s)TREM2. While this approach was intended to holistically capture PL/(s)TREM2 interactions, it effectively diluted our sample size against any specific epitope, potentially impacting the statistical significance of our findings. Future studies could address these challenges by utilizing larger sample sizes of independent MD simulations and more intensive energy calculations as computational power and resources continue to advance.
Moreover, our study, like previous investigations of TREM2, partially assumed that the Ig-like domain is representative of membrane-bound TREM2. This assumption overlooks potential contributions from the membrane, membrane-bound flexible stalk domain, DAP12, or TREM2 glycosylation. Ultimately, more biologically representative simulations of membrane-bound TREM2 are warranted to test this assumption. Additionally, it is unlikely that PLs are presented as individual soluble ligands in the brain in an in vivo setting. Instead, they are more likely co-presented with other components in lipoproteins or in contexts such as demyelination, cell debris, or Aβ. This co-presentation would influence the charge and steric profiles of lipid presentation. While our study provides valuable insights and fills existing knowledge gaps, further research is necessary to fully understand TREM2 binding in more biologically relevant ligand-binding contexts.
Acknowledgements
This work utilized the Alpine high performance computing resource and Blanca condo computing resource at the University of Colorado Boulder. Alpine is jointly funded by the University of Colorado Boulder, the University of Colorado Anschutz, and Colorado State University. Blanca is jointly funded by computing users and the University of Colorado Boulder. This material is based upon work supported by the National Science Foundation Graduate Research Fellowship Program under Fellow ID (2024361899). Any opinions, findings, and conclusions or recommendations expressed in this material are those of the author(s) and do not necessarily reflect the views of the National Science Foundation.
Additional files
Additional information
Funding
National Science Foundation Graduate Research Fellowship Program (2024361899)
References
- 1.Alzheimer’s Association 2024 Alzheimer’s Disease Facts and Figureshttps://www.alz.org/media/Documents/alzheimers-facts-and-figures.pdf
- 2.Inflammation as a central mechanism in Alzheimer’s diseaseAlzheimer’s and Dementia: Translational Research and Clinical Interventions 4https://doi.org/10.1016/j.trci.2018.06.014Google Scholar
- 3.Lipid Droplets in Neurodegenerative DisordersFrontiers in Neuroscience 14https://doi.org/10.3389/fnins.2020.00742Google Scholar
- 4.Microglia in Alzheimer’s disease: pathogenesis, mechanisms, and therapeutic potentialsFrontiers in Aging Neuroscience 15https://doi.org/10.3389/fnagi.2023.1201982Google Scholar
- 5.Astrocytes and microglia in neurodegenerative diseases: Lessons from human in vitro modelsProgress in Neurobiology 200https://doi.org/10.1016/j.pneurobio.2020.101973Google Scholar
- 6.Emerging Alzheimer’s disease therapeutics: promising insights from lipid metabolism and microglia-focused interventionsFrontiers in Aging Neuroscience 15https://doi.org/10.3389/fnagi.2023.1259012Google Scholar
- 7.The Physiology, Pathology, and Potential Therapeutic Applications of the TREM2 Signaling PathwayCell 181https://doi.org/10.1016/j.cell.2020.05.003Google Scholar
- 8.Neurodegenerative Disease–Associated Variants in TREM2 Destabilize the Apical Ligand-Binding Region of the Immunoglobulin DomainFront Neurol 10Google Scholar
- 9.The Alzheimer’s Disease-Associated R47H Variant of TREM2 Has an Altered Glycosylation Pattern and Protein StabilityFront Neurosci 10Google Scholar
- 10.The role of TREM2 N-glycans in trafficking to the cell surface and signal transduction of TREM2The Journal of Biochemistry 172:347–353Google Scholar
- 11.Neurodegenerative disease mutations in TREM2 reveal a functional surface and distinct loss-of-function mechanismseLife 5Google Scholar
- 12.Functional insights from biophysical study of TREM2 interactions with apoE and Aβ1-42Alzheimer’s and Dementia 17:475–488Google Scholar
- 13.Soluble TREM2 inhibits secondary nucleation of Aβ fibrillization and enhances cellular uptake of fibrillar AβProc Natl Acad Sci U S A 119Google Scholar
- 14.Variant of TREM2 Associated with the Risk of Alzheimer’s DiseaseNew England Journal of Medicine 368Google Scholar
- 15.TREM2 Variants in Alzheimer’s DiseaseNew England Journal of Medicine 368Google Scholar
- 16.Evidence of Trem2 variant associated with triple risk of alzheimer’s diseasePLoS One 9Google Scholar
- 17.Molecular Dynamics simulations of Alzheimer’s variants, R47H and R62H, in TREM2 provide evidence for structural alterations behind functional changesbioRxiv https://doi.org/10.1101/536540Google Scholar
- 18.Molecular basis for the loss-of-function effects of the Alzheimer’s disease–associated R47H variant of the immune receptor TREM2Journal of Biological Chemistry 293Google Scholar
- 19.Dynamic insights into the effects of nonsynonymous polymorphisms (nsSNPs) on loss of TREM2 functionSci Rep 12Google Scholar
- 20.TREM2 lipid sensing sustains the microglial response in an Alzheimer’s disease modelCell 160:1061–1071Google Scholar
- 21.TREM2 Is a Receptor for β-Amyloid that Mediates Microglial FunctionNeuron 97Google Scholar
- 22.Amyloid-beta modulates microglial responses by binding to the triggering receptor expressed on myeloid cells 2 (TREM2)Mol Neurodegener 13Google Scholar
- 23.The triggering receptor expressed on myeloid cells 2 binds apolipoprotein EJournal of Biological Chemistry 290:26033–26042Google Scholar
- 24.TREM2 Binds to Apolipoproteins, Including APOE and CLU/APOJ, and Thereby Facilitates Uptake of Amyloid-Beta by MicrogliaNeuron 91Google Scholar
- 25.Apolipoprotein E is a ligand for triggering receptor expressed on myeloid cells 2 (TREM2)Journal of Biological Chemistry 290Google Scholar
- 26.Triggering receptor expressed on myeloid cells 2 (TREM2) deficiency attenuates phagocytic activities of microglia and exacerbates ischemic damage in experimental strokeJournal of Neuroscience 35Google Scholar
- 27.Pattern Recognition by TREM-2: Binding of Anionic LigandsThe Journal of Immunology 171Google Scholar
- 28.Specific lipid recognition is a general feature of CD300 and TREM moleculesImmunogenetics 64Google Scholar
- 29.Lipid processing in the brain: A key regulator of systemic metabolismFrontiers in Endocrinology 8https://doi.org/10.3389/fendo.2017.00060Google Scholar
- 30.Lipid and Lipoprotein Metabolism in MicrogliaFrontiers in Physiology 11https://doi.org/10.3389/fphys.2020.00393Google Scholar
- 31.TREM2 in the pathogenesis of AD: a lipid metabolism regulator and potential metabolic therapeutic targetMolecular Neurodegeneration 17https://doi.org/10.1186/s13024-022-00542-yGoogle Scholar
- 32.TREM2-Ligand Interactions in Health and DiseaseJ Mol Biol 429:1607–1629Google Scholar
- 33.A Unique Microglia Type Associated with Restricting Development of Alzheimer’s DiseaseCell 169Google Scholar
- 34.Soluble TREM2: Innocent bystander or active player in neurological diseases?Neurobiology of Disease 165https://doi.org/10.1016/j.nbd.2022.105630Google Scholar
- 35.A high cerebrospinal fluid soluble TREM2 level is associated with slow clinical progression of Alzheimer’s diseaseAlzheimer’s and Dementia: Diagnosis, Assessment and Disease Monitoring 12Google Scholar
- 36.Higher CSF sTREM2 attenuates ApoE4-related risk for cognitive decline and neurodegenerationMol Neurodegener 15Google Scholar
- 37.Increased soluble TREM2 in cerebrospinal fluid is associated with reduced cognitive and clinical decline in Alzheimer’s diseaseSci Transl Med 11Google Scholar
- 38.sTREM 2 cerebrospinal fluid levels are a potential biomarker for microglia activity in early stage Alzheimer’s disease and associate with neuronal injury markersEMBO Mol Med 8Google Scholar
- 39.Early increase of CSF sTREM2 in Alzheimer’s disease is associated with tau related-neurodegeneration but not with amyloid-β pathologyMol Neurodegener 14Google Scholar
- 40.The relationship of soluble TREM2 to other biomarkers of sporadic Alzheimer’s diseaseSci Rep 11Google Scholar
- 41.Does Soluble TREM2 Protect Against Alzheimer’s Disease?Frontiers in Aging Neuroscience 13https://doi.org/10.3389/fnagi.2021.834697Google Scholar
- 42.Soluble TREM2 induces inflammatory responses and enhances microglial survivalJournal of Experimental Medicine 214Google Scholar
- 43.Soluble TREM2 ameliorates pathological phenotypes by modulating microglial functions in an Alzheimer’s disease modelNat Commun 10Google Scholar
- 44.Wild-type sTREM2 blocks Aβ aggregation and neurotoxicity, but the Alzheimer’s R47H mutant increases Aβ aggregationJournal of Biological Chemistry 296Google Scholar
- 45.AlphaFold Protein Structure Database: Massively expanding the structural coverage of protein-sequence space with high-accuracy modelsNucleic Acids Res 50Google Scholar
- 46.Highly accurate protein structure prediction with AlphaFoldNature 596Google Scholar
- 47.Comparative protein modelling by satisfaction of spatial restraintsJ Mol Biol 234Google Scholar
- 48.CHARMM-GUI membrane builder toward realistic biological membrane simulationsJournal of Computational Chemistry 35https://doi.org/10.1002/jcc.23702Google Scholar
- 49.CHARMM-GUI membrane builder for mixed bilayers and its application to yeast membranesBiophys J 97Google Scholar
- 50.Automated builder and database of protein/membrane complexes for molecular dynamics simulationsPLoS One 2Google Scholar
- 51.CHARMM-GUI Membrane Builder for Complex Biological Membrane Simulations with Glycolipids and LipoglycansJ Chem Theory Comput 15Google Scholar
- 52.CHARMM-GUI Membrane Builder for Lipid Nanoparticles with Ionizable Cationic Lipids and PEGylated LipidsJ Chem Inf Model 61Google Scholar
- 53.CHARMM: The biomolecular simulation programJ Comput Chem 30Google Scholar
- 54.CHARMM-GUI Input Generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM Simulations Using the CHARMM36 Additive Force FieldJ Chem Theory Comput 12Google Scholar
- 55.Canonical sampling through velocity rescalingJournal of Chemical Physics 126Google Scholar
- 56.Molecular dynamics with coupling to an external bathJ Chem Phys 81Google Scholar
- 57.Polymorphic transitions in single crystals: A new molecular dynamics methodJ Appl Phys 52Google Scholar
- 58.Gromacs: High performance molecular simulations through multi-level parallelism from laptops to supercomputersSoftwareX 1–2Google Scholar
- 59.LINCS: A Linear Constraint Solver for molecular simulationsJ Comput Chem 18Google Scholar
- 60.Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems TheStatistical Mechanics of Fluid Mixtures The Journal of Chemical Physics 115Google Scholar
- 61.Software news and updates AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibilityJ Comput Chem 30Google Scholar
- 62.AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreadingJ Comput Chem 31Google Scholar
- 63.AutoDock Vina 1.2.0: New Docking Methods, Expanded Force Field, and Python BindingsJ Chem Inf Model 61Google Scholar
- 64.Electrostatics of nanosystems: Application to microtubules and the ribosomeProc Natl Acad Sci U S A 98Google Scholar
- 65.Gmx_MMPBSA: A New Tool to Perform End-State Free Energy Calculations with GROMACSJ Chem Theory Comput 17:6281–6291Google Scholar
- 66.MMPBSA.py: An efficient program for end-state free energy calculationsJ Chem Theory Comput 8:3314–3321Google Scholar
- 67.Interaction entropy: A new paradigm for highly efficient and reliable computation of protein-ligand binding free energyJ Am Chem Soc 138:5722–5728Google Scholar
- 68.Synergistic reduction in interfacial flexibility of TREM2R47H and ApoE4 may underlie AD pathologyAlzheimer’s and Dementia 21Google Scholar
- 69.Biophysical mapping of TREM2-ligand interactions reveals shared surfaces for engagement of multiple Alzheimer’s disease ligandsMol Neurodegener 20Google Scholar
- 70.TREM2 splicing emerges as crucial aspect to understand TREM2 biologyJournal of Leukocyte Biology 110:827–828https://doi.org/10.1002/JLB.2CE0521-264RGoogle Scholar
- 71.Ligand entropy is hard but should not be ignoredJournal of Chemical Information and Modeling 60https://doi.org/10.1021/acs.jcim.0c01146Google Scholar
- 72.Estimation of conformational entropy in protein-ligand interactions: A computational perspectiveMethods in Molecular Biology 819Google Scholar
- 73.Ligand configurational entropy and protein bindingProc Natl Acad Sci U S A 104Google Scholar
- 74.On the Use of Interaction Entropy and Related Methods to Estimate Binding EntropiesJ Chem Theory Comput 17:5379–5391Google Scholar
- 75.Gene expression and functional deficits underlie TREM2-knockout microglia responses in human models of Alzheimer’s diseaseNat Commun 11Google Scholar
- 76.TREM2 regulates microglial lipid droplet formation and represses post-ischemic brain injuryBiomedicine and Pharmacotherapy 170Google Scholar
- 77.TREM2: Modulator of Lipid Metabolism in MicrogliaNeuron 105https://doi.org/10.1016/j.neuron.2020.02.008Google Scholar
- 78.Exposure of phosphatidylserine on the cell surfaceCell Death and Differentiation 23https://doi.org/10.1038/cdd.2016.7Google Scholar
- 79.TREM2 Agonism with a Monoclonal Antibody Attenuates Tau Pathology and NeurodegenerationCells 12Google Scholar
- 80.Enhancing protective microglial activities with a dual function TREM 2 antibody to the stalk regionEMBO Mol Med 12Google Scholar
- 81.A TREM2-activating antibody with a blood–brain barrier transport vehicle enhances microglial metabolism in Alzheimer’s disease modelsNat Neurosci 26Google Scholar
Article and author information
Author information
Version history
- Sent for peer review:
- Preprint posted:
- Reviewed Preprint version 1:
- Reviewed Preprint version 2:
Cite all versions
You can cite all versions using the DOI https://doi.org/10.7554/eLife.102269. This DOI represents all versions, and will always resolve to the latest one.
Copyright
© 2024, Saeb et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
- views
- 1,706
- downloads
- 37
- citations
- 0
Views, downloads and citations are aggregated across all versions of this paper published by eLife.