Computational Insights Into the Inhibition Mechanism of Proanthocyanidin B2 on Tau Hexapeptide (PHF6) Oligomer

Deciphering the Inhibitory Mechanism of Naphthoquinone-Dopamine on the Aggregation of Tau Core Fragments PHF6* and PHF6

The formation of neurofibrillary tangles by abnormal aggregation of tau protein is considered to be an important pathological characteristic of tauopathies, including Alzheimer's disease and chronic traumatic encephalopathy. Two hexapeptides 275VQIINK280 and 306VQIVYK311 in the microtubule binding region, named PHF6* and PHF6, are known to be aggregation-prone and responsible for tau fibrillization. Previous experiments reported that naphthoquinone-dopamine (NQDA) could effectively inhibit the aggregation of PHF6* and PHF6 and disrupt the fibrillar aggregates into nontoxic species, displaying a dual effect on the amyloid aggregation. However, the underlying molecular mechanism remains mostly elusive. Herein, we performed all-atom molecular dynamics (MD) simulations for 114 μs in total to systematically investigate the impacts of NQDA on the oligomerization of PHF6* and PHF6. The conformational ensembles of PHF6* and PHF6 peptides generated by replica exchange MD simulations show that NQDA could effectively prevent the hydrogen bond formation, reduce the ability of peptides to self-assemble into long β-strand and large β-sheets, and induce peptides to form a loosely packed and coil-rich oligomer. The interaction analysis shows that the binding of NQDA to PHF6* is mainly through hydrophobic interactions with residue I277 and hydrogen bonding interactions with Q276; for the PHF6 peptides, NQDA displays a strong π-π stacking interaction with residue Y310, thus impeding the Y310-Y310 π-π stacking and I308-Y310 CH-π interactions. The DA group of NQDA displays a stronger cation-π interaction than the NQ group, while the NQ group exhibits a stronger π-π stacking interaction. MD simulations demonstrate that NQDA prevents the conformational conversion to β-sheet-rich aggregates and displays an inhibitory effect on the oligomerization dynamics of PHF6* and PHF6. Our results provide a complete picture of inhibitory mechanisms of NQDA on PHF6* and PHF6 oligomerization, which may pave the way for designing drug candidates for the treatment of tauopathies.

On the Tracks of the Aggregation Mechanism of the PHF6 Peptide from Tau Protein: Molecular Dynamics, Energy, and Interaction Network Investigations

The formation of neurofibrillary tangles (NFTs), composed of tau protein aggregates, is a hallmark of some neurodegenerative diseases called tauopathies. NFTs are composed of paired helical filaments (PHFs) of tau protein with a dominant β-sheet secondary structuration. The NFT formation mechanism is not known yet. This study focuses on PHF6, a crucial hexapeptide responsible for tau aggregation. A 2 μs molecular dynamics simulation was launched to determine the keys of the PHF6 aggregation mechanism. Hydrogen bonding, van der Waals, and other non-covalent interactions as π-stacking were investigated. Parallel aggregation was slightly preferred due to its adaptability, but antiparallel aggregation remained widely present during the PHF6 aggregation. The analysis highlighted the leading role of hydrogen bonds identified at the atomic level for each aggregation process. The aggregation study emphasized the importance of Tyr310 during the β-sheets' complexation through π-stacking.

Molecular insights into the allosteric coupling mechanism between an agonist and two different transducers for μ-opioid receptors

G protein-coupled receptors (GPCRs) as the most important class of pharmacological targets regulate G-protein and β-arrestin-mediated signaling through allosteric interplay, which are responsible for different biochemical and physiological actions like therapeutic efficacy and side effects. However, the allosteric mechanism underlying preferentially recruiting one transducer versus the other has been poorly understood, limiting drug design. Motivated by this issue, we utilize accelerated molecular dynamics simulation coupled with potential of mean force (PMF), molecular mechanics Poisson Boltzmann surface area (MM/PBSA) and protein structure network (PSN) to study two ternary complex systems of a representative class A GPCR (μ-opioid receptor (μOR)) bound by an agonist and one specific transducer (G-protein and β-arrestin). The results show that no significant difference exists in the whole structure of μOR between two transducer couplings, but displays transducer-dependent changes in the intracellular binding region of μOR, where the β-arrestin coupling results in a narrower crevice with TM7 inward movement compared with the G-protein. In addition, both the G-protein and β-arrestin coupling can increase the binding affinity of the agonist to the receptor. However, the interactions between the agonist and μOR also exhibit transducer-specific changes, in particular for the interaction with ECL2 that plays an important role in recruiting β-arrestin. The allosteric network analysis further indicates that Y148^3.33, F152^3.37, F156^3.41, N191^4.49, T160^3.45, Y106^2.42, W293^6.48, F289^6.44, I248^5.54 and Y252^5.58 play important roles in equally activating G-protein and β-arrestin. In contrast, M161^3.46 and R165^3.50 devote important contributions to preferentially recruit G-protein while D164^3.49 and R179^ICL2 are revealed to be important for selectively activating β-arrestin. The observations provide useful information for understanding the biased activation mechanism.