Next Article in Journal
Molecular Mechanisms of Treatment Resistance in Glioblastoma
Next Article in Special Issue
Structure and Dynamics of Meprin β in Complex with a Hydroxamate-Based Inhibitor
Previous Article in Journal
Enhanced Suppression of Immune Cells In Vitro by MSC Overexpressing FasL
Previous Article in Special Issue
Solution Structures and Dynamic Assembly of the 24-Meric Plasmodial Pdx1–Pdx2 Complex
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 22
Issue 1
10.3390/ijms22010349
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pick’s Tau Fibril Shows Multiple Distinct PET Probe Binding Sites: Insights from Computational Modelling

by
Sushil K. Mishra
1,
Yoshiki Yamaguchi
2,*,
Makoto Higuchi
3 and
Naruhiko Sahara
3,*
1
Advance Glycoscience Research Cluster, National University of Ireland Galway, H91 W2TY Galway, Ireland
2
Faculty of Pharmaceutical Sciences, Tohoku Medical and Pharmaceutical University, Sendai, Miyagi 981-8558, Japan
3
National Institute of Radiological Sciences, National Institutes for Quantum and Radiological Science and Technology, Chiba 263-8555, Japan
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2021, 22(1), 349; https://doi.org/10.3390/ijms22010349
Submission received: 8 December 2020 / Revised: 27 December 2020 / Accepted: 27 December 2020 / Published: 31 December 2020
(This article belongs to the Special Issue Proteins and Protein-Ligand Interactions)
Browse Figures
Review Reports Versions Notes

Abstract

:
In recent years, it has been realized that the tau protein is a key player in multiple neurodegenerative diseases. Positron emission tomography (PET) radiotracers that bind to tau filaments in Alzheimer’s disease (AD) are in common use, but PET tracers binding to tau filaments of rarer, age-related dementias, such as Pick’s disease, have not been widely explored. To design disease-specific and tau-selective PET tracers, it is important to determine where and how PET tracers bind to tau filaments. In this paper, we present the first molecular modelling study on PET probe binding to the structured core of tau filaments from a patient with Pick’s disease (TauPiD). We have used docking, molecular dynamics simulations, binding-affinity and tunnel calculations to explore TauPiD binding sites, binding modes, and binding energies of PET probes (AV-1451, MK-6240, PBB3, PM-PBB3, THK-5351 and PiB) with TauPiD. The probes bind to TauPiD at multiple surface binding sites as well as in a cavity binding site. The probes show unique surface binding patterns, and, out of them all, PM-PBB3 proves to bind the strongest. The findings suggest that our computational workflow of structural and dynamic details of the tau filaments has potential for the rational design of TauPiD specific PET tracers.

1. Introduction

Filamentous tau aggregates are the histopathological hallmark of a group of neurodegenerative diseases including Alzheimer’s disease (AD), progressive supranuclear palsy (PSP), corticobasal degeneration (CBD), Pick’s disease (PiD), and familial frontotemporal lobar degeneration with underlying tau pathology (FTLD-Tau) collectively referred as tauopathies [1]. The conformational diversity of tau assemblies derived from distinct isoform compositions, posttranslational modifications and interactions with other molecules results in variations in pathological tau depositions in different subcellular compartments, cell types, and brain regions [2]. The composition of tau fibrils in each of the tauopathies is well characterized biochemically and by pathological examination. Paired helical filaments (PHFs) in AD are composed of all six tau isoforms, whereas straight filaments (SFs) and twisted ribbons are made up of isoforms with four repeat (4R) domains in PSP and CBD, and with three repeat (3R) domains in PiD [3]. Electron cryo-microscopy (cryo-EM) has revealed high-resolution structures of the C-terminal core region in tau filaments from AD, PiD, CBD and chronic traumatic encephalopathy (CTE) brains [4,5,6]. Analysis of samples, from multiple cases of sporadic and inherited AD, has revealed a common structure of the filament cores [7]. The TauAD filaments are composed of two C-shaped protofilaments arranged either symmetrically base-to-base (PHF) or back-to-base (SF). The cores in both types of filament consist of 73 residues (V306–F378) with eight β-strands (β1–β8). On the other hand, distinct tau filament structures are found in PiD, CBD or CTE subjects [5,6,8]. TauPiD core is formed of 94 residues (K254 to F378), and is, therefore, longer than TauAD. A total of nine β-strands (β1–β9) is arranged in four cross-β packing stacks in TauPiD. These cryo-EM studies raise the interesting question as to whether different tau filament structures can be distinguished by tau tracers, which have recently been developed for in vivo positron emission tomography (PET) imaging.
Tau tracers are mostly designed according to β-sheet binding properties and theoretically label all filamentous tau aggregates (e.g., neurofibrillary tangles (NFTs), neuropil threads, tufted astrocytes, astrocytic plaques, and coiled bodies). Several PET tracers for detecting filamentous tau inclusions, such as [11C]PBB3, [18F]PM-PBB3, [18F]AV1451, [18F]THK5351, [18F]MK-6240, [18F]R06958948, [18F]GTP-1, and [18F]PI-2620, were developed and have been applied to human subjects [9,10,11,12,13,14,15,16]. The distribution of the bound tracers recapitulated Braak NFT staging [17] in AD brains, suggesting in vivo imaging-based staging of tau pathologies. On the other hand, [18F]AV1451 and [18F]THK5351 exhibit off-target effects such as binding to monoamine oxidase (MAO)-A and MAO-B, respectively [18,19]. In particular, the binding of [18F]THK5351 to MAO-B is strongly suggested in the detection of MAO-B-positive astrogliosis [20,21,22]. Second generation tracers ([18F]PM-PBB3, [18F]MK-6240, [18F]R06958948, [18F]GTP-1, and [18F]PI-2620) seem to have less off-target binding [12,23,24]. More recently, the binding properties of these PET tracers with PHFs and SFs were examined by computational modeling using structural information from cryo-EM studies [25]. TauPiD folds in a J-shaped structure, instead of the C shape seen in patients with AD [4,8]. The modeling identified three to four potential high-affinity binding sites for these tracers with some preference for certain sites over others for each tracer [25,26]. Importantly, validation for the off-target interaction of tau tracers with amyloid β fibrils was also undertaken [25]. However, modeling of PET tracer binding sites in non-AD tauopathies has not been explored.
Of the various neurodegenerative diseases, Pick’s disease is a rare type of age-related dementia that affects the frontal lobes of the brain, resulting in speech problems like aphasia, behavior difficulties and eventually death. Pick’s tau folds quite differently from Alzheimer’s tau filaments [4,8], and its surface residues where the tracers potentially bind are very different too. Since the binding specificity of current tau PET probes have not been completely validated against tau pathologies, approaches employing the structure-guided design of tracers for different filaments can be extremely useful. In order to elucidate the binding mechanism of PET probes in TauPiD, we performed computational modeling to characterize the binding of various PET probes to TauPiD, and obtained insights into their unique binding properties.

2. Results and Discussion

2.1. Assessing Binding Regions and Affinity of PET Tracers to TauPiD by Docking

Current tau PET tracers such as [18F]AV-1451 ([18F]flortaucipir), [18F]THK5351, [11C]PBB3, [18F]PM-PBB3 and [18F]MK-6240 (Figure 1A) have been successful in the in vivo evaluation of tau pathology of AD brains [9,11,16,27,28]. [18F]PM-PBB3 PET imaging also was demonstrated to be useful for the differential diagnosis of non-AD dementia patients (e.g., PSP, CBD, and Pick disease) [16]. However, binding of other tracers to tau filaments derived from non-AD brains is still under investigation. For evaluation of each of the tracers, we first tested our docking strategy on TauPiD using the structural information on TauPiD from the Protein Data Bank (PDB ID: 6GX5) [8]. PET tracers are designed to bind to cross-β structure [29]; therefore, we assumed that PET tracers bind to the C-terminal core region of tau filaments and not to the other flexible region [4,8]. The docking results show that a number of surface binding sites (S1–S11 in Figure 1B) in TauPiD can potentially bind PET probes. AutoDock Vina docking scores (hereafter called docking scores) and X-Score rescoring energies (hereafter called rescoring energies) of the highest-scoring pose in each binding site are listed in Table 1. The inner cavity formed between P312-V313-D314 and P332-G333-G334 (C1 in Figure 1B) showed the strongest affinity for all PET probes.
Below, we primarily focus on the solvent-accessible surface binding of PET probes and then examine the cavity binding using a bigger tau assembly. In the case of the surface binding sites, docking scores and rescoring energies of the same binding pose show a similar trend for all the PET tracers.
The predicted docking scores of the PET probes at surface sites are in the micromolar range (binding free energy of −8.2 kcal·mol−1 corresponds to Kd of 10 μM), and are in agreement with other docking studies [25]. Docking scores for MK-6240 (−7.4 kcal·mol−1 in S6), PM-PBB3 (−6.1 kcal·mol−1 in S8), and PiB (−7.1 kcal·mol−1 in S7) at surface binding sites are weaker compared to AV-1451 (−9.4 kcal·mol−1 in S7), PBB3 (−8.4 kcal·mol−1 in S7) and THK-5351 (−8.3 kcal·mol−1 in S7) (Table 1). Although PiB is used for imaging β-amyloid (Aβ) plaques, this compound has a docking score similar to that of the Tau PET probes (MK-6240 and PM-PBB3).
Rescoring energies predict PBB3 to be the most potent binder at a surface binding site (S7; −8.7 kcal·mol−1). Otherwise, rescoring energies for the surface binding sites trend like docking scores, in that AV-1451 (−8.4 kcal·mol−1 in S7) and THK-5351 (−8.3 kcal·mol−1 in S7) bind with similar strength and somewhat more strongly than PM-PBB3 (−7.5 kcal·mol−1 in S8), MK-6240 (−7.7 kcal·mol−1 in S6) and PiB (−7.8 kcal·mol−1 in S7) (Table 1). Among all surface sites, PM-PBB3, THK-5351 and PiB have the weakest rescoring energies, which is almost in line with the docking scores (Table 1). Four PET probes, not MK-6240 and PM-PBB3, prefer surface binding site S7. PM-PBB3 and MK-6240 bind neighboring S8 and S6, respectively, with the highest affinity. Only PBB3 and PiB bind to S1 and S3, but with weak affinity. S11 uniquely binds only PM-PBB3 and PiB.
Docking scores and rescoring energies are in good agreement when the scores of individual PET probes are compared in multiple binding sites. However, there are some inconsistencies when comparing the scores of different PET probes with each other. Since docking is performed to a rigid protein, it fails to take into account protein flexibility and solvent effects, and the resulting docking poses can be artificial [30]. Therefore, we further validated the predicted docking poses by MD simulations and binding free-energy calculations.

2.2. Molecular Dynamics and MM/GBSA Calculations

We performed MD simulations of each PET probe bound to TauPiD simultaneously in a single simulation over a range of high to low affinity binding sites as predicted by docking. Since in each individual MD simulation, the same PET probe is present in all predicted binding sites, neighboring PET probes can influence binding, as expected in reality. In the TauPiD/PET probe complexes, compounds present at surface binding sites sometimes drifted out from their binding site to bulk solvent and then started interacting with another probe that is already bound to tau (Supplementary Figures S1 and S2). Therefore, for the Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) binding energy calculation in each site, PET probes in all the remaining binding sites were stripped off by pre-processing the MD trajectory, otherwise the energy contribution arising from probe-probe interaction could be misleading. Thus, there is no influence of neighboring probes in the MM/GBSA binding-energy calculation and the presented values solely reflect intermolecular interactions between PET probe and tau at that particular binding site.
TauPiD was quite stable during the MD simulations, except that the fibrils at both edges start to unfold during at extended times. Unfolding could be due to interaction of the fibril edges with the solvent. In a proper tau assembly, filaments would not be exposed to solvent; instead, they are stacked between other fibrils. Since TauPiD fibrils exist in different polymorphic forms, fibrils are expected to be highly flexible. The MD-based MM/GBSA approach used here for binding free-energy calculations takes into account tracer and receptor flexibility issues and, therefore, is expected to provide more realistic binding energies than docking or rescoring functions. We refer to binding enthalpies as binding energies here, because we assume that entropic contributions do not play a major role in predicting relative binding order for the same system. The computed binding energies of the PET probes in each binding site are shown in Figure 2 and listed in Table 2.
MM/GBSA binding calculations on 1000 frames suggest that Pick’s tau core has multiple surface binding sites (see Figure 1B) where S7 is the most potent for AV-1451 (−25.5 kcal·mol−1), THK-5351 (−22.8 kcal·mol−1) and PiB (−25.9 kcal·mol−1), and S8 for PBB3 (−27.9 kcal·mol−1) and PM-PBB3 (−34.2 kcal·mol−1) (Table 2). PBB3, which has been used to visualize tau deposition in AD and PiD, shows binding to multiple TauPiD surface sites with almost equal strength (−26.0 kcal·mol−1 in S1; −26.2 kcal·mol−1 in S5; −24.3 kcal·mol−1 in S6; −27.6 kcal·mol−1 in S7 and −27.9 kcal·mol−1 in S8) and could be a potent TauPiD tracer (Figure 2). PM-PBB3 also binds in a site (S2n) neighboring site S2 with binding energy −30.2 ± 3.2 kcal·mol−1 (Table 2). PM-PBB3 also docked to site S3n (neighbor of S3) but binding energies could not be extracted due to van der Waals clashes in some of the snapshots. Site S8 is unique for both PBB3 and PM-PBB3 as only these show binding here (all other PET tracers either failed to bind here or they drifted away from the binding site during the MD simulation) (Table 2, Figure 2).
MK-6240 showed very weak affinity for surface sites, which is due to the fact that the probe drifted away from all the predicted surface sites during the simulation (Figure 2). Unlike other PET probes, the structure of MK-6240 is not planar in 3D and cannot bind to narrow hydrophobic patches formed between two amino acid side chains. MK-6240 seems an unlikely candidate as a TauPiD tracer. PBB3 and PM-PBB3 bind most strongly to TauPiD, either at S7 or neighboring sites. AV-1451, THK-5351 and PiB effectively bind to S7 only; other binding sites are too weak to be significant. It has been reported that PBB3 strongly binds to Pick bodies in brain slices while AV1451 only weakly binds [31]. Thus, S8 might be the high-affinity site for PBB3 and PM-PBB3 under physiological condition, which needs to be verified in future structural study.
MM/GBSA calculations have the capability to provide per-residue decomposition of the energy contribution to the binding free-energy, and this can provide insight into the binding of the PET tracers to tau, namely that it is primarily derived from CH/π stacking interactions (Supplementary Table S1). For example, in sites S7 and S8, PET probes bind to hydrophobic patches formed by sidechains of K347 and R349, respectively, via stacking interactions between hydrogen atoms of the amino-acid side chain and π electrons in the PET tracers (Figure 3). The CH/π stacking interactions are a strong attractive molecular force occurring between a soft acid and a soft base and are frequently seen in protein–glycan complexes [32,33,34].
In general, the binding pattern shown by the MM/GBSA energies is quite similar to the docking scores in most cases, but there are significant differences in some. Surface binding of PET probes at S7 is predicted to be the strongest by both docking and MM/GBSA calculations. On the other hand, PBB3 binding to S1 is predicted to be very weak by docking but quite strong by MM/GBSA energies. Two factors likely cause these discrepancies. Firstly, due to the narrow range of binding scores of AutoDock Vina, computational differentiation between weak and strong binding sites is constrained. Thus, a small error in docking scores can affect the rankings. Secondly, the differences in MM/GBSA energies are within the range of standard deviation of the energies over the 1000 frames used in the calculations. In this context, MM/GBSA energies are a more useful parameter as they also provide a meaningful statistic that can be taken into account to rank the binding strength of PET tracers in multiple binding sites.

2.3. Comparison of Probe Binding Affinities between TauAD and TauPiD

Figure 4A shows PET probe binding sites in TauAD first reported by Goedert et al., [26] and soon after by Murugan et al. [25], as well as the sites on TauPiD from this study. From a structural stance, it can be seen that the binding sites of the majority of the PET probes in TauAD and TauPiD are formed in cavities surrounded by a basic residue (Lys, Arg) either from one or from both sides (Figure 4A,B). Docking and MM/GBSA energies of PBB3, THK5351 and MK-6240 with TauAD were published by Murugan et al. [25] and are here compared with TauPiD. Docking scores and MM/GBSA approach-based binding energies of PBB3, THK5351 and MK-6240 at their strongest surface binding site within TauAD and TauPiD are shown in Figure 4C–F. The figures suggest that these PET probes have a greater binding preference for the TauAD fibril than that of TauPiD (Figure 4E,F). The study by Murugan et al. [25] lacks the name of the docking software and details of the Generalized Born (GB) model used in the MM/GBSA calculation. Since different models in AMBER18 lead to a significantly different range of absolute binding energies, the deviations between binding energies for two different types of tau filaments could be due to a difference in the length of the MD simulation (100 ns) or to a different GB model chosen for the MM/GBSA calculations. We chose longer MD simulations (100 ns) to allow sufficient time for dissociation. Despite expected differences in absolute energies, the relative binding order of three common PET tracers is the same in both studies (PBB3 > THK5351 > MK-6240) (Figure 4C,D). Binding energies from the MM/GBSA calculation also show the same binding preference (PBB3 > THK5351 > MK-6240) in TauAD and TauPiD (Figure 4E,F), suggesting that a discussion of relative binding order is still valid despite a possible different choice of MM/GBSA parameters.

2.4. PM-PBB3 as a Potent PET Tracer for Detecting TauPiD

Most recent findings with both clinical and pre-clinical studies of [18F]PM-PBB3 imaging revealed the feasibility of using PM-PBB3 to recognize tau pathologies in both AD and non-AD brains including Pick’s disease [16]. To confirm PM-PBB3 binding to TauAD and TauPiD and to determine preference through binding affinity, we calculated binding affinities of PM-PBB3 for both TauAD and TauPiD by docking and MM/GBSA calculations. Absolute MM/GBSA energies of PM-PBB3 are similar for both TauAD and TauPiD (−34.2 kcal·mol−1 for TauAD and −34.3 kcal·mol−1 for TauPiD) (Figure 5A,B), confirming that PM-PBB3 binds to both kinds of tau filaments. This is in the agreement with recent experimental study that indicates binding of PM-PBB3 to TauPiD and TauAD is nearly similar [16]. On the other hand, the relative binding order of PM-PBB3 illustrates unique features of the TauPiD protofibril when compared to that of TauAD. Among PET tracers used in this study, the binding energy of PM-PBB3 relative to PBB3, THK-5351, MK-6240 and other PETs is relatively strong for TauAD compared to TauPiD. Since docking cannot predict experimental binding trends of PM-PBB3 so accurately, MD based free-energy approaches are important to consider in the rational design of probes for disease specific tau filaments. In addition, this approach can be applied to the interaction of PET tracers with non-tauopathy off-targets, such as MAO-A/-B and other dementia/neurodegenerative proteins (amyloid-β or α-synuclein). Since previous studies confirmed relatively low affinities of tau tracers, [11C]PBB3 and [18F]PM-PBB3, for MAO-A/B and α-synuclein by in vitro radioligand binding assay [16,35,36], it will be interesting to perform side-by-side validation with both an in vitro assay and an MM/GBSA approach.
Our simulation studies demonstrate that PM-PBB3 has the highest MM/GBSA binding energy of the seven compounds tested, and revealed that the surface binding sites in TauPiD are rather unique. S8 site is highly selective for PM-PBB3 over PBB3 (Figure 2). In addition to these predicted features, PM-PBB3 has several practical advantages, e.g., greater metabolic stability compared with PBB3 [16]. Furthermore, a 18F-labeled tracer has the advantage of broader availability and higher throughput than 11C-labeled tracers. Taking all these points together, PM-PBB3 stands out as the most promising tracer for diagnostic imaging of diverse neurodegenerative tauopathies.

2.5. Does Pick’s Tau have a Potential Cavity Binding Site?

Unlike TauAD, TauPiD shows a potential cavity binding-site (C1) where PET probes bind within a cavity formed between P312-V313-D314 and P332-G333-G334 (Figure 1B). The PET tracers AV-1451 (−11.0 kcal·mol−1), MK-6240 (−11.0 kcal·mol−1), PBB3 (−10.1 kcal·mol−1), PM-PBB3 (−10.0 kcal·mol−1), and THK-5351 (−10.3 kcal·mol−1) have comparable docking scores in the C1 site, but PiB (−8.9 kcal·mol−1) binds comparatively weakly. Of all the docking scores, those of PiB in C1 and at the surface sites are weakest. AV-1451, PBB3, and THK-5351 bind with similar strength to C1, but MK-6240 and PiB show weak affinity. Binding energies from MD simulations followed by MM/GBSA calculations on 1000 frames suggest that binding site C1, which is formed due to folding of the fibril, is a favored binding site. All seven compounds show stable binding to C1, and the binding is stronger than it is to surface binding sites S1 to S11 (Figure 6). This observation is similar to PET probe binding to β -amyloid, where cavity sites also bind to PET probers with higher affinity than surface sites [25]. Of the tracers, PBB3 binds strongest to C1—an interesting observation as the binding is stronger than PiB. The C1 site primarily consists of hydrophobic interactions mediated by the hydrogen atoms of P312 and P332.
To further explore the cavity site C1 unique to TauPiD, we performed 1-μs MD simulation of a bigger tau assembly composed of a 49 polypeptide (Tau49PiD) without any PET probes in the binding sites. We further performed a Tunnel analysis around C1, in a 1-μs MD trajectory to see if the tunnel is realistic and maintained during the MD simulation. The tunnel profile in the cryo-EM structure of TauPiD shows two potential tunnels with some slightly narrow entry paths between filaments to C1 (Figure 7A and Supplementary Figures S3–S5). We further evaluated protein tunnels on 500 conformations of a bigger TauPiD assembly (Tau49PiD) taken at every 2 ns from 1-μs long MD simulation (Figure 7A–F). Tunnel lengths and radii in the MD-simulated structure vary during the simulation, but are very narrow compared to the cryo-EM structure. A cluster analysis of all the tunnels in those 500 snapshots shows clusters mainly in two distinct regions in the plot (Figure 7G). The majority of tunnels cluster in a region varying in length from 150 to 300 Å and number 60 to 150. The distance between P213 residues at both ends of the Tau49PiD is ~230 Å, indicating that some of the tunnels are >230 Å long. There may be interlinking tunnels (Supplementary Figure S4). The other cluster shows smaller tunnels of less than 50 Å in length but higher in number. Therefore, the majority of conformers possess a tunnel of about 150 to 250 Å in length. The radius profile suggests that tunnels shrink during MD simulation (Supplementary Figure S4). However, a number of snapshots still show continuous tunnels connecting both ends of the assembly, often merging to one or exiting to the surface.
It is essential to mention here that tunnel finding algorithms have limitations and it is expected that a program will identify an exit towards the surface where the sphere radius for moving forward in the tunnel is lower than the radius to exit towards the surface. Some of the MD frames show tunnels starting at one end of the assembly and finishing at the other end (Figure 7E). However, very often a tunnel either ends inside the assembly or finds an alternative path through a neighboring tunnel (Figure 7F) and then connects back to the original tunnel in the C1 region.
Our overall analysis of tunnels in the MD simulation snapshots suggests that there may be a potential cavity binding site in TauPiD. However, very narrow circumferences in several snapshots and reductions in tunnel radii in the large tau assembly also suggest that the broad cavernous region around C1 could be a result of limitations in the cryo-EM structure. The cryo-EM structure of TauPiD was solved with a resolution of 3.2 Å and it has a distance of 8.5 Å between the P312 and P332 sidechains (Figure 7H) where C1 is formed. Ideally, the cavity binding site needs to be verified experimentally. Recent cryo-EM studies of CBD and CTE tau filaments indicate the presence of additional non-proteinous densities inside the cavity [5,6]. Our computational approach will also contribute to the understanding of such internal binding of non-proteinous molecules once the structure is identified.
In summary, the computational workflow presented in this study can be a basis for the rational design of tracers specific to a particular tau fiber. Currently, no atomic data are available to experimentally show the binding mode between PET tracer and tau fibril. Docking studies are complementary to 3D structural analyses, such as cryo-EM, which usually provide averaged electron densities at low-to-medium resolutions. The strong point of the MM/GBSA approach is the ability to estimate individual binding energies (binding constants) of multiple binding sites such as those that exist in the tau fiber. Experimental approaches cannot provide binding constants for each of such multiple binding sites. Although a gap still exists between theoretical predictions and experimental results, a computational approach will ensure advances in elucidating tracer-protein interactions of multiple binding sites.

3. Materials and Methods

3.1. Structure Preparation

The structures of TauPiD (PDB ID: 6GX5) and TauAD (PDB ID: 5O3L) were taken from the Protein Data Bank [8]. TauPiD and TauAD structures comprise 94 and 73 amino acids, respectively, from the total 441 residues in the 2N4R tau sequence. The reported TauPiD assembly has three tau fibrils (K254–F378), so we prepared a five-fibril assembly by superimposing structures over each other to have a structure comparable to the one used in the TauAD studies [25]. Five model compounds representing PET tracers, AV-1451, MK-6240, PBB3, PM-PBB3, and THK-5351, and the Pittsburgh compound B (PiB), were prepared using Maestro [37]. PET structures were optimized geometrically initially using the quantum mechanics (QM) approach (HF/6-31+G(d,p) level of theory and then further at B3LYP/6-31G*) in Gaussian09 [38]. The Restrained Electrostatic Potential (RESP) atomic partial charges were calculated using program antechamber of Amber18 [39] with ESP potentials from Gaussian.

3.2. Docking

The tau assembly was centered in a conformational search space of 30 Å × 30 Å × 70 Å. We chose this size to ensure that it was large enough to capture all the surface sites in the fibrils. A total of 100 docking conformations of each PET probe was obtained by software AutoDock Vina-carb using chi_coeff = 0 and chi_cutoff = 12. The used values of chi_coeff of and chi_cutoff switch Vina-carb to AutoDock Vina [40] are denoted as Vina throughout the manuscript. An exhaustiveness of 20 and energy range cutoff of 8 kcal·mol−1 were used. All the binding poses were refined by a rescoring function of X-Score [41]. The docking population of PET probes was calculated for each site, and all top scoring docking poses from each cluster were further analyzed by extensive MD simulations and binding free energy calculations. Possible binding sites were labelled as cavity site 1 (C1), and surface sites 1 to 11 (S1 to S11).

3.3. MD Simulation

A TauPiD assembly with bound PET probe in all potential binding sites was prepared. For each PET probe, we performed MD simulation of fully solvated assemblies where the PET probe is simultaneously bound to TauPiD in all the binding sites predicted by docking. For each binding site, the highest-scoring binding pose from each cluster was selected for MD simulation. The protein was treated with Amber ff14SB force field [42], whereas PET probes were modelled using General Amber Force Field v.2.0 (GAFF2) [43] with RESP charges. All the complexes were solvated by an octahedral TIP3P water box extending 12 Å from each edge of the protein. A total of 45 Cl− ions was added to neutralize the whole system. A multi-step protocol published elsewhere [44] was used to equilibrate the complexes. Finally, a 100 nanosecond MD of each complex was performed at NPT, using MD settings: temperature at 300 K, temperature scaling by Langevin dynamics (collision frequency = 2), pressure relaxation every 1.2 picosecond, SHAKE constraints, nonbonded interaction cutoff of 10 Å, and 2 fs integration time step. Similarly, MD of TauAD complexed with PM-PBB3 in five binding sites was performed.
The TauPiD cryo-EM structure possesses a core cavity at a site between β-sheets containing residues P312 and G332, which can accommodate PET tracers. To check the validity of the existence of the core cavity (whether it is real or an artefact resulting from the cryo-electron microscopy approach), a 230 Å long 49 fibril TauPiD assembly (Tau49PiD) was prepared for MD simulation. The MD of Tau49PiD in apo form was performed on the microsecond timescale. In contrast to the MD of five fibril TauPiD, this Tau49PiD assembly was solvated with a rectangular water box, extending 10 Å from the protein surface, and comprised a total of ~80,000 water molecules and 441 Cl− ions. Prepared systems were quite large (320,000 and 400,000 atoms for the apo and PET probe-bound TauPiD structures, respectively). MD of Tau49PiD fibrils was performed for 1 µs in explicit solvent using cuda version of pmemd in Amber18 [39]. The non-bonded interaction cut-off was set to 9 Å for these systems. The other MD settings were as described before. We analyzed the MD trajectories using cpptraj [45].

3.4. MM/GBSA Calculations

To calculate binding free energy of the PET probe in each site, solvent molecules and tracers present in remaining binding sites were removed from the trajectory. A total of 1000 snapshots were extracted after each 10 ps of a 100-ns MD simulation and used for the PET probe binding energy calculation with the MM/GBSA approach [46]. A similar procedure was repeated for the PET probe in other binding sites one by one. We equate binding enthalpies as binding energies assuming that entropic contributions do not play a key role in predicting relative order of binding for the same system [25,47]. The GBOBC generalized-born model [48] was used with set mbondi2 radii. Salt concentration was 0.15 M. The surface tension and non-polar solvation free energy correction terms were set to 0.005 kcal·mol−1 and 0.0, respectively, for the solvent accessible surface area (SASA) calculation. An offset value of 0.0 was used to correct non-polar solvation free energy contributions. The interior and exterior dielectric constants were set to 1 and 80, respectively. Other parameters were their default values in Amber18.

3.5. Tunnel Calculations

Tunnels around the core-site in TauPiD were calculated using the CAVER Analyst 2.0 program [49] with a probe radius of 0.9 Å and shell radius of 3.0 Å, shell depth of 4.0 Å and a clustering threshold of 3.5 Å. Tunnels were calculated on 500 frames taken after every 2 ns from the 1-µs simulation.

Supplementary Materials

The following are available online at https://www.mdpi.com/1422-0067/22/1/349/s1.

Author Contributions

Conceptualization, Y.Y. and N.S.; methodology, S.K.M.; data analysis, S.K.M.; writing—original draft preparation, S.K.M.; writing—review and editing, S.K.M., Y.Y., M.H. and N.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by Grants-in-And for Scientific Research on Innovative Areas “Brain Protein Aging and Dementia Control” (KAKENHI grant number JP26117008 to N.S. and Y.Y.), “Singularity Biology” (19H05437 to N.S.), for Scientific Research (C) (19K06896 to N.S.) and for Brain Mapping by Integrated Neurotechnologies for Disease Studies from the Japan Agency for Medical Research and Development (19dm0207072 and JP20dm0207072 to M.H.).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available upon request from the corresponding authors.

Acknowledgments

S.K.M. acknowledges the Irish Centre for High-End Computing (ICHEC) for the provision of computational facilities for part of this work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lee, V.M.; Goedert, M.; Trojanowski, J.Q. Neurodegenerative Tauopathies. Annu. Rev. Neurosci. 2001, 24, 1121–1159. [Google Scholar] [CrossRef] [PubMed]
  2. Jadhav, S.; Avila, J.; Schöll, M.; Kovacs, G.G.; Kövari, E.; Skrabana, R.; Evans, L.D.; Kontsekova, E.; Malawska, B.; de Silva, R.; et al. A Walk through Tau Therapeutic Strategies. Acta Neuropathol. Commun. 2019, 7, 22. [Google Scholar] [CrossRef] [PubMed]
  3. Goedert, M.; Eisenberg, D.S.; Crowther, R.A. Propagation of Tau Aggregates and Neurodegeneration. Annu. Rev. Neurosci. 2017, 40, 189–210. [Google Scholar] [CrossRef] [PubMed]
  4. Fitzpatrick, A.W.P.; Falcon, B.; He, S.; Murzin, A.G.; Murshudov, G.; Garringer, H.J.; Crowther, R.A.; Ghetti, B.; Goedert, M.; Scheres, S.H.W. Cryo-EM Structures of Tau Filaments from Alzheimer’s Disease. Nature 2017, 547, 185–190. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Falcon, B.; Zivanov, J.; Zhang, W.; Murzin, A.G.; Garringer, H.J.; Vidal, R.; Crowther, R.A.; Newell, K.L.; Ghetti, B.; Goedert, M.; et al. Novel Tau Filament Fold in Chronic Traumatic Encephalopathy Encloses Hydrophobic Molecules. Nature 2019, 568, 420–423. [Google Scholar] [CrossRef] [PubMed]
  6. Zhang, W.; Tarutani, A.; Newell, K.L.; Murzin, A.G.; Matsubara, T.; Falcon, B.; Vidal, R.; Garringer, H.J.; Shi, Y.; Ikeuchi, T.; et al. Novel Tau Filament Fold in Corticobasal Degeneration. Nature 2020, 580, 283–287. [Google Scholar] [CrossRef]
  7. Falcon, B.; Zhang, W.; Schweighauser, M.; Murzin, A.G.; Vidal, R.; Garringer, H.J.; Ghetti, B.; Scheres, S.H.W.; Goedert, M. Tau Filaments from Multiple Cases of Sporadic and Inherited Alzheimer’s Disease Adopt a Common Fold. Acta Neuropathol. 2018, 136, 699–708. [Google Scholar] [CrossRef] [Green Version]
  8. Falcon, B.; Zhang, W.; Murzin, A.G.; Murshudov, G.; Garringer, H.J.; Vidal, R.; Crowther, R.A.; Ghetti, B.; Scheres, S.H.W.; Goedert, M. Structures of Filaments from Pick’s Disease Reveal a Novel Tau Protein Fold. Nature 2018, 561, 137–140. [Google Scholar] [CrossRef]
  9. Maruyama, M.; Shimada, H.; Suhara, T.; Shinotoh, H.; Ji, B.; Maeda, J.; Zhang, M.-R.; Trojanowski, J.Q.; Lee, V.M.-Y.; Ono, M.; et al. Imaging of Tau Pathology in a Tauopathy Mouse Model and in Alzheimer Patients Compared to Normal Controls. Neuron 2013, 79, 1094–1108. [Google Scholar] [CrossRef] [Green Version]
  10. Xia, C.-F.; Arteaga, J.; Chen, G.; Gangadharmath, U.; Gomez, L.F.; Kasi, D.; Lam, C.; Liang, Q.; Liu, C.; Mocharla, V.P.; et al. [18F]T807, a Novel Tau Positron Emission Tomography Imaging Agent for Alzheimer’s Disease. Alzheimers Dement. 2013, 9, 666–676. [Google Scholar] [CrossRef]
  11. Hostetler, E.D.; Walji, A.M.; Zeng, Z.; Miller, P.; Bennacef, I.; Salinas, C.; Connolly, B.; Gantert, L.; Haley, H.; Holahan, M.; et al. Preclinical Characterization of 18F-MK-6240, a Promising PET Tracer for In Vivo Quantification of Human Neurofibrillary Tangles. J. Nucl. Med. 2016, 57, 1599–1606. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Walji, A.M.; Hostetler, E.D.; Selnick, H.; Zeng, Z.; Miller, P.; Bennacef, I.; Salinas, C.; Connolly, B.; Gantert, L.; Holahan, M.; et al. Discovery of 6-(Fluoro-18F)-3-(1H-Pyrrolo[2,3-c]Pyridin-1-Yl)Isoquinolin-5-Amine ([18F]-MK-6240): A Positron Emission Tomography (PET) Imaging Agent for Quantification of Neurofibrillary Tangles (NFTs). J. Med. Chem. 2016, 59, 4778–4789. [Google Scholar] [CrossRef] [PubMed]
  13. Hall, B.; Mak, E.; Cervenka, S.; Aigbirhio, F.I.; Rowe, J.B.; O’Brien, J.T. In Vivo Tau PET Imaging in Dementia: Pathophysiology, Radiotracer Quantification, and a Systematic Review of Clinical Findings. Ageing Res. Rev. 2017, 36, 50–63. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Mueller, A.; Kroth, H.; Schieferstein, H.; Berndt, M.; Oden, F.; Capotosti, F.; Molette, J.; Juergens, T.; Darmency, V.; Schmitt-Willich, H.; et al. Preclinical Characterization of PI-2620, A Novel Tau PET Tracer for Detection of Tau in AD and Other Tauopathies. Alzheimer’s Dement. J. Alzheimer’s Assoc. 2017, 13, P141–P142. [Google Scholar] [CrossRef]
  15. Saint-Aubert, L.; Lemoine, L.; Chiotis, K.; Leuzy, A.; Rodriguez-Vieitez, E.; Nordberg, A. Tau PET Imaging: Present and Future Directions. Mol. Neurodegener. 2017, 12, 19. [Google Scholar] [CrossRef] [Green Version]
  16. Tagai, K.; Ono, M.; Kubota, M.; Kitamura, S.; Takahata, K.; Seki, C.; Takado, Y.; Shinotoh, H.; Sano, Y.; Yamamoto, Y.; et al. High-Contrast In Vivo Imaging of Tau Pathologies in Alzheimer’s and Non-Alzheimer’s Disease Tauopathies. Neuron 2020. [Google Scholar] [CrossRef]
  17. Braak, H.; Braak, E. Diagnostic Criteria for Neuropathologic Assessment of Alzheimer’s Disease. Neurobiol. Aging 1997, 18, S85–S88. [Google Scholar] [CrossRef]
  18. Lowe, V.J.; Curran, G.; Fang, P.; Liesinger, A.M.; Josephs, K.A.; Parisi, J.E.; Kantarci, K.; Boeve, B.F.; Pandey, M.K.; Bruinsma, T.; et al. An Autoradiographic Evaluation of AV-1451 Tau PET in Dementia. Acta Neuropathol. Commun. 2016, 4, 58. [Google Scholar] [CrossRef] [Green Version]
  19. Ng, K.P.; Pascoal, T.A.; Mathotaarachchi, S.; Therriault, J.; Kang, M.S.; Shin, M.; Guiot, M.-C.; Guo, Q.; Harada, R.; Comley, R.A.; et al. Monoamine Oxidase B Inhibitor, Selegiline, Reduces 18F-THK5351 Uptake in the Human Brain. Alzheimers Res. 2017, 9, 25. [Google Scholar] [CrossRef] [Green Version]
  20. Harada, R.; Ishiki, A.; Kai, H.; Sato, N.; Furukawa, K.; Furumoto, S.; Tago, T.; Tomita, N.; Watanuki, S.; Hiraoka, K.; et al. Correlations of 18F-THK5351 PET with Postmortem Burden of Tau and Astrogliosis in Alzheimer Disease. J. Nucl. Med. 2018, 59, 671–674. [Google Scholar] [CrossRef] [Green Version]
  21. Ishiki, A.; Harada, R.; Kai, H.; Sato, N.; Totsune, T.; Tomita, N.; Watanuki, S.; Hiraoka, K.; Ishikawa, Y.; Funaki, Y.; et al. Neuroimaging-Pathological Correlations of [18F]THK5351 PET in Progressive Supranuclear Palsy. Acta Neuropathol. Commun. 2018, 6, 53. [Google Scholar] [CrossRef]
  22. Ishibashi, K.; Kameyama, M.; Tago, T.; Toyohara, J.; Ishii, K. Potential Use of 18F-THK5351 PET to Identify Wallerian Degeneration of the Pyramidal Tract Caused by Cerebral Infarction. Clin. Nucl. Med. 2017, 42, e523–e524. [Google Scholar] [CrossRef] [PubMed]
  23. Gobbi, L.C.; Knust, H.; Körner, M.; Honer, M.; Czech, C.; Belli, S.; Muri, D.; Edelmann, M.R.; Hartung, T.; Erbsmehl, I.; et al. Identification of Three Novel Radiotracers for Imaging Aggregated Tau in Alzheimer’s Disease with Positron Emission Tomography. J. Med. Chem. 2017, 60, 7350–7370. [Google Scholar] [CrossRef] [PubMed]
  24. Kroth, H.; Oden, F.; Molette, J.; Schieferstein, H.; Capotosti, F.; Mueller, A.; Berndt, M.; Schmitt-Willich, H.; Darmency, V.; Gabellieri, E.; et al. Discovery and Preclinical Characterization of [18F]PI-2620, a next-Generation Tau PET Tracer for the Assessment of Tau Pathology in Alzheimer’s Disease and Other Tauopathies. Eur. J. Nucl. Med. Mol. Imaging 2019, 46, 2178–2189. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Murugan, N.A.; Nordberg, A.; Ågren, H. Different Positron Emission Tomography Tau Tracers Bind to Multiple Binding Sites on the Tau Fibril: Insight from Computational Modeling. ACS Chem. Neurosci. 2018, 9, 1757–1767. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Goedert, M.; Yamaguchi, Y.; Mishra, S.K.; Higuchi, M.; Sahara, N. Tau Filaments and the Development of Positron Emission Tomography Tracers. Front. Neurol. 2018, 9, 70. [Google Scholar] [CrossRef] [Green Version]
  27. Chien, D.T.; Bahri, S.; Szardenings, A.K.; Walsh, J.C.; Mu, F.; Su, M.-Y.; Shankle, W.R.; Elizarov, A.; Kolb, H.C. Early Clinical PET Imaging Results with the Novel PHF-Tau Radioligand [F-18]-T807. J. Alzheimers Dis. 2013, 34, 457–468. [Google Scholar] [CrossRef] [Green Version]
  28. Betthauser, T.J.; Hillmer, A.T.; Lao, P.J.; Ehlerding, E.; Mukherjee, J.; Stone, C.K.; Christian, B.T. Human Biodistribution and Dosimetry of [18F]Nifene, an A4β2* Nicotinic Acetylcholine Receptor PET Tracer. Nucl. Med. Biol. 2017, 55, 7–11. [Google Scholar] [CrossRef]
  29. Harada, R.; Okamura, N.; Furumoto, S.; Yanai, K. Imaging Protein Misfolding in the Brain Using β-Sheet Ligands. Front. Neurosci. 2018, 12. [Google Scholar] [CrossRef] [Green Version]
  30. Chen, Y.-C. Beware of Docking! Trends Pharm. Sci. 2015, 36, 78–95. [Google Scholar] [CrossRef]
  31. Ono, M.; Sahara, N.; Kumata, K.; Ji, B.; Ni, R.; Koga, S.; Dickson, D.W.; Trojanowski, J.Q.; Lee, V.M.-Y.; Yoshida, M.; et al. Distinct Binding of PET Ligands PBB3 and AV-1451 to Tau Fibril Strains in Neurodegenerative Tauopathies. Brain 2017, 140, 764–780. [Google Scholar] [CrossRef] [Green Version]
  32. Hudson, K.L.; Bartlett, G.J.; Diehl, R.C.; Agirre, J.; Gallagher, T.; Kiessling, L.L.; Woolfson, D.N. Carbohydrate-Aromatic Interactions in Proteins. J. Am. Chem. Soc. 2015, 137, 15152–15160. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Wimmerova, M.; Kozmon, S.; Nečasová, I.; Mishra, S.K.; Komárek, J.; Koca, J. Stacking Interactions between Carbohydrate and Protein Quantified by Combination of Theoretical and Experimental Methods. PLoS ONE 2012, 7, e46032. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Houser, J.; Kozmon, S.; Mishra, D.; Mishra, S.K.; Romano, P.R.; Wimmerová, M.; Koča, J. Influence of Trp Flipping on Carbohydrate Binding in Lectins. An Example on Aleuria Aurantia Lectin AAL. PLoS ONE 2017, 12, e0189375. [Google Scholar] [CrossRef] [Green Version]
  35. Ni, R.; Ji, B.; Ono, M.; Sahara, N.; Zhang, M.-R.; Aoki, I.; Nordberg, A.; Suhara, T.; Higuchi, M. Comparative In Vitro and In Vivo Quantifications of Pathologic Tau Deposits and Their Association with Neurodegeneration in Tauopathy Mouse Models. J. Nucl. Med. 2018, 59, 960–966. [Google Scholar] [CrossRef] [Green Version]
  36. Ono, M.; Takahashi, M.; Shimozawa, A.; Fujinaga, M.; Mori, W.; Nagai, Y.; Mimura, K.; Minamihisamatsu, T.; Uchida, S.; Shimojo, M.; et al. In Vivo Visualization of Propagating α-Synuclein Pathologies in Mouse and Marmoset Models by a Bimodal Imaging Probe, C05-05. bioRxiv 2020. [Google Scholar] [CrossRef]
  37. Maestro; Schrödinger, LLC: New York, NY, USA, 2018.
  38. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian˜16 Revision B.01; Gaussian, Inc.: Wallingford, CT, USA, 2016. [Google Scholar]
  39. Case, D.A.; Belfon, K.; Ben-Shalom, I.Y.; Brozell, S.R.; Cerutti, D.S.; Cheatham, T.E., III; Cruzeiro, V.W.D.; Darden, T.A.; Duke, R.E.; Giambasu, G.; et al. AMBER 18; University of California: San Fransisco, CA, USA, 2018. [Google Scholar]
  40. Trott, O.; Olson, A.J. AutoDock Vina: Improving the Speed and Accuracy of Docking with a New Scoring Function, Efficient Optimization, and Multithreading. J. Comput. Chem. 2010, 31, 455–461. [Google Scholar] [CrossRef] [Green Version]
  41. Wang, R.; Lai, L.; Wang, S. Further Development and Validation of Empirical Scoring Functions for Structure-Based Binding Affinity Prediction. J. Comput. Aided Mol. Des. 2002, 16, 11–26. [Google Scholar] [CrossRef]
  42. Maier, J.A.; Martinez, C.; Kasavajhala, K.; Wickstrom, L.; Hauser, K.E.; Simmerling, C. Ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from Ff99SB. J. Chem. Theory Comput. 2015, 11, 3696–3713. [Google Scholar] [CrossRef] [Green Version]
  43. Wang, J.; Wolf, R.M.; Caldwell, J.W.; Kollman, P.A.; Case, D.A. Development and Testing of a General Amber Force Field. J. Comput. Chem. 2004, 25, 1157–1174. [Google Scholar] [CrossRef]
  44. Nagae, M.; Mishra, S.K.; Hanashima, S.; Tateno, H.; Yamaguchi, Y. Distinct Roles for Each N-Glycan Branch Interacting with Mannose-Binding Type Jacalin-Related Lectins Orysata and Calsepa. Glycobiology 2017, 27, 1120–1133. [Google Scholar] [CrossRef] [PubMed]
  45. Roe, D.R.; Cheatham, T.E. PTRAJ and CPPTRAJ: Software for Processing and Analysis of Molecular Dynamics Trajectory Data. J. Chem. Theory Comput. 2013, 9, 3084–3095. [Google Scholar] [CrossRef] [PubMed]
  46. Mishra, S.K.; Koča, J. Assessing the Performance of MM/PBSA, MM/GBSA, and QM–MM/GBSA Approaches on Protein/Carbohydrate Complexes: Effect of Implicit Solvent Models, QM Methods, and Entropic Contributions. J. Phys. Chem. B 2018, 122, 8113–8121. [Google Scholar] [CrossRef] [PubMed]
  47. Genheden, S.; Ryde, U. The MM/PBSA and MM/GBSA Methods to Estimate Ligand-Binding Affinities. Expert Opin. Drug Discov. 2015, 10, 449–461. [Google Scholar] [CrossRef] [PubMed]
  48. Okamoto, Y. Generalized-Ensemble Algorithms: Enhanced Sampling Techniques for Monte Carlo and Molecular Dynamics Simulations. J. Mol. Graph. Model. 2004, 22, 425–439. [Google Scholar] [CrossRef] [Green Version]
  49. Jurcik, A.; Bednar, D.; Byska, J.; Marques, S.M.; Furmanova, K.; Daniel, L.; Kokkonen, P.; Brezovsky, J.; Strnad, O.; Stourac, J.; et al. CAVER Analyst 2.0: Analysis and Visualization of Channels and Tunnels in Protein Structures and Molecular Dynamics Trajectories. Bioinformatics 2018, 34, 3586–3588. [Google Scholar] [CrossRef] [Green Version]
Ijms 22 00349 g001 550
Figure 1. Positron emission tomography (PET) probes and their binding sites in tau predicted by molecular docking. (A) Chemical structures of unlabeled PET probes studied in this paper. (B) Secondary structure of tau filaments (K254 to F378) from Pick’s patients showing surface binding sites S1 to S11 (green), cavity site C1 (blue) and neighboring amino acid residues. Each β strand is labeled from β1 to β9.
Figure 1. Positron emission tomography (PET) probes and their binding sites in tau predicted by molecular docking. (A) Chemical structures of unlabeled PET probes studied in this paper. (B) Secondary structure of tau filaments (K254 to F378) from Pick’s patients showing surface binding sites S1 to S11 (green), cavity site C1 (blue) and neighboring amino acid residues. Each β strand is labeled from β1 to β9.
Ijms 22 00349 g001
Ijms 22 00349 g002 550
Figure 2. Graph of MM/GBSA binding energies of PET probes in surface binding site (S1 to S11).
Figure 2. Graph of MM/GBSA binding energies of PET probes in surface binding site (S1 to S11).
Ijms 22 00349 g002
Ijms 22 00349 g003 550
Figure 3. Molecular interactions of PET probes to TauPiD taken after 100-ns MD simulation. (A) AV-1451 and (B) PBB3 binding to R349 and K347 sidechain hydrogens in site S7 and (C) PM-PBB3 binding to R349 and Q351 sidechain hydrogens in site S8.
Figure 3. Molecular interactions of PET probes to TauPiD taken after 100-ns MD simulation. (A) AV-1451 and (B) PBB3 binding to R349 and K347 sidechain hydrogens in site S7 and (C) PM-PBB3 binding to R349 and Q351 sidechain hydrogens in site S8.
Ijms 22 00349 g003
Ijms 22 00349 g004 550
Figure 4. Comparison of the surface binding sites in TauAD reported previously (A) and TauPiD studied in this paper (B). Surface representation for tau structures is colored as residue type (acidic-red, basic-blue, nonpolar-grey, polar-green, cysteine-yellow). Docking scores of the strongest docking pose for PBB3, THK-5351 and MK-6240 in the best scoring site in TauAD reported in Murugan et al. [25] (C) and TauPiD (D). Similarity of MM/GBSA binding energy (kcal·mol−1) of the strongest binding modes of PET probes in TauAD (E) and TauPiD (F).
Figure 4. Comparison of the surface binding sites in TauAD reported previously (A) and TauPiD studied in this paper (B). Surface representation for tau structures is colored as residue type (acidic-red, basic-blue, nonpolar-grey, polar-green, cysteine-yellow). Docking scores of the strongest docking pose for PBB3, THK-5351 and MK-6240 in the best scoring site in TauAD reported in Murugan et al. [25] (C) and TauPiD (D). Similarity of MM/GBSA binding energy (kcal·mol−1) of the strongest binding modes of PET probes in TauAD (E) and TauPiD (F).
Ijms 22 00349 g004
Ijms 22 00349 g005 550
Figure 5. Docking scores and binding affinities of PM-PBB3 in TauAD (A) and TauPiD (B). Standard deviation in MM/GBSA binding energy over all snapshots is shown as error bars.
Figure 5. Docking scores and binding affinities of PM-PBB3 in TauAD (A) and TauPiD (B). Standard deviation in MM/GBSA binding energy over all snapshots is shown as error bars.
Ijms 22 00349 g005
Ijms 22 00349 g006 550
Figure 6. MM/GBSA binding energies of PET probes in site C1 and surface binding sites that show the most potent binding (S1–S11).
Figure 6. MM/GBSA binding energies of PET probes in site C1 and surface binding sites that show the most potent binding (S1–S11).
Ijms 22 00349 g006
Ijms 22 00349 g007 550
Figure 7. Tunnels predicted by the CAVER program in (A) cryo-EM structure and (B–F) structure taken from MD simulation. The snapshots taken after (B) 200 ns, (C) 400 ns, (D) 600 ns, (E) 800 ns, (F) 1 microsecond simulation showing tunnels around the cavity site. (G) Plot of the tunnel length and number of tunnels found by CAVER on 500 frames (extracted each at 2 ns from 1 microsecond MD simulation) of the unbound tau filament. (H) Surface of the TauPiD structure showing cavities in the structure and distance between neighboring Pro residues.
Figure 7. Tunnels predicted by the CAVER program in (A) cryo-EM structure and (B–F) structure taken from MD simulation. The snapshots taken after (B) 200 ns, (C) 400 ns, (D) 600 ns, (E) 800 ns, (F) 1 microsecond simulation showing tunnels around the cavity site. (G) Plot of the tunnel length and number of tunnels found by CAVER on 500 frames (extracted each at 2 ns from 1 microsecond MD simulation) of the unbound tau filament. (H) Surface of the TauPiD structure showing cavities in the structure and distance between neighboring Pro residues.
Ijms 22 00349 g007
Table
Table 1. Calculated docking scores (kcal·mol−1) from molecular docking of various PET probes in the different binding sites of TauPiD protofibril. Values in brackets are the binding energies from the rescoring function X-Score for the same binding pose. Site C1 corresponds to a cavity site while the remainder (S1 to S11) correspond to surface binding sites. Corresponding binding sites are shown in Figure 1B.
Table 1. Calculated docking scores (kcal·mol−1) from molecular docking of various PET probes in the different binding sites of TauPiD protofibril. Values in brackets are the binding energies from the rescoring function X-Score for the same binding pose. Site C1 corresponds to a cavity site while the remainder (S1 to S11) correspond to surface binding sites. Corresponding binding sites are shown in Figure 1B.
Binding
Site		AV-1451MK-6240PBB3PM-PBB3THK5351PiB
C1−11.0 (−8.9)−11.0 (−9.0)−10.1 (−9.4)−10.0 (−9.4)−10.3 (−9.0)−8.9 (−8.5)
S1NBNB−4.9 (−7.0)NBNB−4.4 (−6.6)
S2−6.0 (−7.3)NBNB−5.0 (−7.2)−5.7 (−7.0)−4.4 (−6.9)
S3NBNB−5.7 (−7.6)−5.4 (−7.4)NB−4.6 (−7.2)
S4NB−5.1 (−6.7)−5.0 (−6.9)NBNB−4.5 (−6.5)
S5−6.5 (−7.4)−5.5 (−6.7)−6.0 (−7.7)NB−5.8 (−7.1)−5.2 (−7.2)
S6−8.3 (−7.9)−7.4 (−7.7)−7.3 (−8.0)NB−7.3 (−7.8)−6.3 (−7.5)
S7−9.4 (−8.4)NB−8.4 (−8.7)NB−8.3 (−8.3)−7.1 (−7.8)
S8−7.1 (−7.5)NB−6.2 (−7.6)−6.1 (−7.5)NBNB
S9−5.9 (−6.7)−5.7 (−6.9)NB−4.5 (−6.5)−5.6 (−6.7)−4.6 (−6.6)
S10−6.0 (−7.1)NBNB−4.6 (−6.6)−5.2 (−6.8)−4.6 (−6.7)
S11NBNBNB−4.5 (−6.9)NB−4.5 (−6.9)
PM-PBB3 also binds to a site in the left neighbor of site S2 (S2n) with a docking score of −4.6 (−6.7) and a site in the right neighbor of site S3 (S3n) with a docking score of −3.7 (−6.8). NB; did not bind. Lowest scores among all the surface binding sites are in bold.
Table
Table 2. The calculated MM/GBSA binding energies (kcal·mol−1) of PET probes in TauPiD binding sites. The label “x” represents sites where PET probe was not docked and a label “-” shows cases where MM/GBSA binding energies could not be calculated to due to tracer drift or steric clashes in certain snapshots. All the energy values are reported in kcal·mol−1.
Table 2. The calculated MM/GBSA binding energies (kcal·mol−1) of PET probes in TauPiD binding sites. The label “x” represents sites where PET probe was not docked and a label “-” shows cases where MM/GBSA binding energies could not be calculated to due to tracer drift or steric clashes in certain snapshots. All the energy values are reported in kcal·mol−1.
Binding
Site		AV-1451MK-6240PBB3PM-PBB3THK-5351PiB
C1−25.8 ± 2.8−32.1 ± 4.1−40.9 ± 2.7−50.0 ± 3.1−33.2 ± 5.1−27.0 ± 3.0
S1xx−26.0 ± 5.5-x−21.4 ± 2.9
S2−10.7 ± 1.6xx−12.7 ± 4.6−12.6 ± 1.8−13.1 ± 2.8
S3xx−13.6 ± 1.6 −11.2 ± 5.1x−10.7 ± 1.7
S4x−3.9 ± 3.7−16.9 ± 2.6−18.9 ± 2.4x−5.5 ± 5.7
S5−19.9 ± 1.9−12.5 ± 3.3−26.2 ± 3.3x−11.3 ± 19.3−13.8 ± 3.4
S6−13.8 ± 6.9-−24.3 ± 2.4x−19.3 ± 2.9−13.7 ± 4.8
S7−25.5 ± 2.5x−27.6 ± 2.6x−22.8 ± 3.2−25.9 ± 2.6
S8−3.2 ± 5.0x−27.9 ± 2.9−34.2 ± 3.3xx
S9-−7.1 ± 5.7x -−10.7 ± 6.5
S10-xx-−17.4 ± 8.7−5.6 ± 6.9
S11xxx-x−16.8 ± 25.8
MM/GBSA binding energy of PM-PBB3 in S2n site (binding site in the left neighbor of S2) (S2n) is −30.2 ± 3.2 kcal·mol−1. PM-PBB3 was docked to a site S3n (a site in the right neighbor of site S3 but MM/GBSA energies could not be calculated due to steric clashes. Lowest scores among all the surface binding sites are in bold.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mishra, S.K.; Yamaguchi, Y.; Higuchi, M.; Sahara, N. Pick’s Tau Fibril Shows Multiple Distinct PET Probe Binding Sites: Insights from Computational Modelling. Int. J. Mol. Sci. 2021, 22, 349. https://doi.org/10.3390/ijms22010349

AMA Style

Mishra SK, Yamaguchi Y, Higuchi M, Sahara N. Pick’s Tau Fibril Shows Multiple Distinct PET Probe Binding Sites: Insights from Computational Modelling. International Journal of Molecular Sciences. 2021; 22(1):349. https://doi.org/10.3390/ijms22010349

Chicago/Turabian Style

Mishra, Sushil K., Yoshiki Yamaguchi, Makoto Higuchi, and Naruhiko Sahara. 2021. "Pick’s Tau Fibril Shows Multiple Distinct PET Probe Binding Sites: Insights from Computational Modelling" International Journal of Molecular Sciences 22, no. 1: 349. https://doi.org/10.3390/ijms22010349

APA Style

Mishra, S. K., Yamaguchi, Y., Higuchi, M., & Sahara, N. (2021). Pick’s Tau Fibril Shows Multiple Distinct PET Probe Binding Sites: Insights from Computational Modelling. International Journal of Molecular Sciences, 22(1), 349. https://doi.org/10.3390/ijms22010349

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop