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,1113. 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 factor1416. In vitro and in silico studies have consistently revealed that mutations, including R47H, destabilize TREM2’s CDRs8,1719, 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,2325, 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 progression3537. 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 progression3840.

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-GUI4852, 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 implementation6466. 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.

sTREM2s 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 emerged7981, 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

Supplemental Information

Additional information

Funding

National Science Foundation Graduate Research Fellowship Program (2024361899)