Molecular recognition of morphine and fentanyl by the human µ-opioid receptor.

Morphine and fentanyl are among the most used opioid drugs that confer analgesia and unwanted side effects through both G protein and arrestin signaling pathways of µ-opioid receptor (µOR). Here, we report structures of the human µOR-G protein complexes bound to morphine and fentanyl, which uncover key differences in how they bind the receptor. We also report structures of µOR bound to TRV130, PZM21, and SR17018, which reveal preferential interactions of these agonists with TM3 side of the ligand-binding pocket rather than TM6/7 side. In contrast, morphine and fentanyl form dual interactions with both TM3 and TM6/7 regions. Mutations at the TM6/7 interface abolish arrestin recruitment of µOR promoted by morphine and fentanyl. Ligands designed to reduce TM6/7 interactions display preferential G protein signaling. Our results provide crucial insights into fentanyl recognition and signaling of µOR, which may facilitate rational design of next-generation analgesics.

Dimerization and antidepressant recognition at noradrenaline transporter.

The noradrenaline transporter has a pivotal role in regulating neurotransmitter balance and is crucial for normal physiology and neurobiology1. Dysfunction of noradrenaline transporter has been implicated in numerous neuropsychiatric diseases, including depression and attention deficit hyperactivity disorder2. Here we report cryo-electron microscopy structures of noradrenaline transporter in apo and substrate-bound forms, and as complexes with six antidepressants. The structures reveal a noradrenaline transporter dimer interface that is mediated predominantly by cholesterol and lipid molecules. The substrate noradrenaline binds deep in the central binding pocket, and its amine group interacts with a conserved aspartate residue. Our structures also provide insight into antidepressant recognition and monoamine transporter selectivity. Together, these findings advance our understanding of noradrenaline transporter regulation and inhibition, and provide templates for designing improved antidepressants to treat neuropsychiatric disorders.

GPCR activation and GRK2 assembly by a biased intracellular agonist.

Phosphorylation of G-protein-coupled receptors (GPCRs) by GPCR kinases (GRKs) desensitizes G-protein signalling and promotes arrestin signalling, which is also modulated by biased ligands1-6. The molecular assembly of GRKs on GPCRs and the basis of GRK-mediated biased signalling remain largely unknown owing to the weak GPCR-GRK interactions. Here we report the complex structure of neurotensin receptor 1 (NTSR1) bound to GRK2, Gαq and the arrestin-biased ligand SBI-5537. The density map reveals the arrangement of the intact GRK2 with the receptor, with the N-terminal helix of GRK2 docking into the open cytoplasmic pocket formed by the outward movement of the receptor transmembrane helix 6, analogous to the binding of the G protein to the receptor. SBI-553 binds at the interface between GRK2 and NTSR1 to enhance GRK2 binding. The binding mode of SBI-553 is compatible with arrestin binding but clashes with the binding of Gαq protein, thus providing a mechanism for its arrestin-biased signalling capability. In sum, our structure provides a rational model for understanding the details of GPCR-GRK interactions and GRK2-mediated biased signalling.

Molecular Basis for Hormone Recognition and Activation of Corticotropin-Releasing Factor Receptors.

Corticotropin-releasing factor (CRF) and the three related peptides urocortins 1-3 (UCN1-UCN3) are endocrine hormones that control the stress responses by activating CRF1R and CRF2R, two members of class B G-protein-coupled receptors (GPCRs). Here, we present two cryoelectron microscopy (cryo-EM) structures of UCN1-bound CRF1R and CRF2R with the stimulatory G protein. In both structures, UCN1 adopts a single straight helix with its N terminus dipped into the receptor transmembrane bundle. Although the peptide-binding residues in CRF1R and CRF2R are different from other members of class B GPCRs, the residues involved in receptor activation and G protein coupling are conserved. In addition, both structures reveal bound cholesterol molecules to the receptor transmembrane helices. Our structures define the basis of ligand-binding specificity in the CRF receptor-hormone system, establish a common mechanism of class B GPCR activation and G protein coupling, and provide a paradigm for studying membrane protein-lipid interactions for class B GPCRs.

Structures of Gi-bound metabotropic glutamate receptors mGlu2 and mGlu4.

The metabotropic glutamate receptors (mGlus) have key roles in modulating cell excitability and synaptic transmission in response to glutamate (the main excitatory neurotransmitter in the central nervous system)1. It has previously been suggested that only one receptor subunit within an mGlu homodimer is responsible for coupling to G protein during receptor activation2. However, the molecular mechanism that underlies the asymmetric signalling of mGlus remains unknown. Here we report two cryo-electron microscopy structures of human mGlu2 and mGlu4 bound to heterotrimeric Gi protein. The structures reveal a G-protein-binding site formed by three intracellular loops and helices III and IV that is distinct from the corresponding binding site in all of the other G-protein-coupled receptor (GPCR) structures. Furthermore, we observed an asymmetric dimer interface of the transmembrane domain of the receptor in the two mGlu-Gi structures. We confirmed that the asymmetric dimerization is crucial for receptor activation, which was supported by functional data; this dimerization may provide a molecular basis for the asymmetric signal transduction of mGlus. These findings offer insights into receptor signalling of class C GPCRs.

Identification of CAF-related lncRNAs at the pan-cancer level represents a potential carcinogenic risk.

Cancer-associated fibroblasts (CAFs) are increasingly recognized as playing a crucial role in regulating cancer progression and metastasis. These cells can be activated by long non-coding RNAs (lncRNAs), promoting the malignant biological processes of tumor cells. Therefore, it is essential to understand the regulatory relationship between CAFs and lncRNAs in cancers. Here, we identified CAF-related lncRNAs at the pan-cancer level to systematically predict their potential regulatory functions. The identified lncRNAs were also validated using various external data at both tissue and cellular levels. This study has revealed that these CAF-related lncRNAs exhibit expression perturbations in cancers and are highly correlated with the infiltration of stromal cells, particularly fibroblasts and endothelial cells. By prioritizing a list of CAF-related lncRNAs, we can further distinguish patient subtypes that show survival and molecular differences. In addition, we have developed a web server, CAFLnc (https://46906u5t63.zicp.fun/CAFLnc/), to visualize our results. In conclusion, CAF-related lncRNAs hold great potential as a valuable resource for comprehending lncRNA functions and advancing the identification of biomarkers for cancer progression and therapeutic targets in cancer treatment.

Structure of GPR101-Gs enables identification of ligands with rejuvenating potential.

GPR101 is an orphan G protein-coupled receptor actively participating in energy homeostasis. Here we report the cryo-electron microscopy structure of GPR101 constitutively coupled to Gs heterotrimer, which reveals unique features of GPR101, including the interaction of extracellular loop 2 within the 7TM bundle, a hydrophobic chain packing-mediated activation mechanism and the structural basis of disease-related mutants. Importantly, a side pocket is identified in GPR101 that facilitates in silico screening to identify four small-molecule agonists, including AA-14. The structure of AA-14-GPR101-Gs provides direct evidence of the AA-14 binding at the side pocket. Functionally, AA-14 partially restores the functions of GH/IGF-1 axis and exhibits several rejuvenating effects in wild-type mice, which are abrogated in Gpr101-deficient mice. In summary, we provide a structural basis for the constitutive activity of GPR101. The structure-facilitated identification of GPR101 agonists and functional analysis suggest that targeting this orphan receptor has rejuvenating potential.

LATS2 condensates organize signalosomes for Hippo pathway signal transduction.

Biomolecular condensates have been proposed to mediate cellular signaling transduction. However, the mechanism and functional consequences of signal condensates are not well understood. Here we report that LATS2, the core kinase of the Hippo pathway, responds to F-actin cytoskeleton reduction and forms condensates. The proline-rich motif (PRM) of LATS2 mediates its condensation. LATS2 partitions with the main components of the Hippo pathway to assemble a signalosome for LATS2 activation and for its stability by physically compartmentalizing from E3 ligase FBXL16 complex-dependent degradation, which in turn mediates yes-associated protein (YAP)-transcriptional coactivator with PDZ-binding motif (TAZ) recruitment and inactivation. This oncogenic FBXL16 complex blocks LATS2 condensation by binding to the PRM region to promote its degradation. Disruption of LATS2 condensation leads to tumor progression. Thus, our study uncovers that the signalosomes assembled by LATS2 condensation provide a compartmentalized and reversible platform for Hippo signaling transduction and protein stability, which have potential implications in cancer diagnosis and therapeutics.

Activation of the GLP-1 receptor by a non-peptidic agonist.

Class B G-protein-coupled receptors are major targets for the treatment of chronic diseases, including diabetes and obesity1. Structures of active receptors reveal peptide agonists engage deep within the receptor core, leading to an outward movement of extracellular loop 3 and the tops of transmembrane helices 6 and 7, an inward movement of transmembrane helix 1, reorganization of extracellular loop 2 and outward movement of the intracellular side of transmembrane helix 6, resulting in G-protein interaction and activation2-6. Here we solved the structure of a non-peptide agonist, TT-OAD2, bound to the glucagon-like peptide-1 (GLP-1) receptor. Our structure identified an unpredicted non-peptide agonist-binding pocket in which reorganization of extracellular loop 3 and transmembrane helices 6 and 7 manifests independently of direct ligand interaction within the deep transmembrane domain pocket. TT-OAD2 exhibits biased agonism, and kinetics of G-protein activation and signalling that are distinct from peptide agonists. Within the structure, TT-OAD2 protrudes beyond the receptor core to interact with the lipid or detergent, providing an explanation for the distinct activation kinetics that may contribute to the clinical efficacy of this compound series. This work alters our understanding of the events that drive the activation of class B receptors.

Structural analysis of the dual agonism at GLP-1R and GCGR.

Glucagon-like peptide-1 receptor (GLP-1R) and glucagon receptor (GCGR), two members of class B1 G protein-coupled receptors, play important roles in glucose homeostasis and energy metabolism. They share a high degree of sequence homology but have different functionalities. Unimolecular dual agonists of both receptors developed recently displayed better clinical efficacies than that of monotherapy. To study the underlying molecular mechanisms, we determined high-resolution cryo-electron microscopy structures of GLP-1R or GCGR in complex with heterotrimeric Gs protein and three GLP-1R/GCGR dual agonists including peptide 15, MEDI0382 (cotadutide) and SAR425899 with variable activating profiles at GLP-1R versus GCGR. Compared with related structures reported previously and supported by our published pharmacological data, key residues responsible for ligand recognition and dual agonism were identified. Analyses of peptide conformational features revealed a difference in side chain orientations within the first three residues, indicating that distinct engagements in the deep binding pocket are required to achieve receptor selectivity. The middle region recognizes extracellular loop 1 (ECL1), ECL2, and the top of transmembrane helix 1 (TM1) resulting in specific conformational changes of both ligand and receptor, especially the dual agonists reshaped ECL1 conformation of GLP-1R relative to that of GCGR, suggesting an important role of ECL1 interaction in executing dual agonism. Structural investigation of lipid modification showed a better interaction between lipid moiety of MEDI0382 and TM1-TM2 cleft, in line with its increased potency at GCGR than SAR425899. Together, the results provide insightful information for the design and development of improved therapeutics targeting these two receptors simultaneously.

Structural insights into the interaction of three Y-shaped ligands with PI3Kα.

Class IA phosphoinositide 3-kinase alpha (PI3Kα) is an important drug target because it is one of the most frequently mutated proteins in human cancers. However, small molecule inhibitors currently on the market or under development have safety concerns due to a lack of selectivity. Therefore, other chemical scaffolds or unique mechanisms of catalytic kinase inhibition are needed. Here, we report the cryo-electron microscopy structures of wild-type PI3Kα, the dimer of p110α and p85α, in complex with three Y-shaped ligands [cpd16 (compound 16), cpd17 (compound 17), and cpd18 (compound 18)] of different affinities and no inhibitory effect on the kinase activity. Unlike ATP-competitive inhibitors, cpd17 adopts a Y-shaped conformation with one arm inserted into a binding pocket formed by R770 and W780 and the other arm lodged in the ATP-binding pocket at an angle that is different from that of the ATP phosphate tail. Such a special interaction induces a conformation of PI3Kα resembling that of the unliganded protein. These observations were confirmed with two isomers (cpd16 and cpd18). Further analysis of these Y-shaped ligands revealed the structural basis of differential binding affinities caused by stereo- or regiochemical modifications. Our results may offer a different direction toward the design of therapeutic agents against PI3Kα.

Molecular basis of signal transduction mediated by the human GIPR splice variants.

Glucose-dependent insulinotropic polypeptide receptor (GIPR) is a potential drug target for metabolic disorders. It works with glucagon-like peptide-1 receptor and glucagon receptor in humans to maintain glucose homeostasis. Unlike the other two receptors, GIPR has at least 13 reported splice variants (SVs), more than half of which have sequence variations at either C or N terminus. To explore their roles in endogenous peptide-mediated GIPR signaling, we determined the cryoelectron microscopy (cryo-EM) structures of the two N terminus-altered SVs (referred as GIPR-202 and GIPR-209 in the Ensembl database, SV1 and SV2 here, respectively) and investigated the outcome of coexpressing each of them in question with GIPR in HEK293T cells with respect to ligand binding, receptor expression, cAMP (adenosine 3,5-cyclic monophosphate) accumulation, ß-arrestin recruitment, and cell surface localization. It was found that while both N terminus-altered SVs of GIPR neither bound to the hormone nor elicited signal transduction per se, they suppressed ligand binding and cAMP accumulation of GIPR. Meanwhile, SV1 reduced GIPR-mediated ß-arrestin 2 responses. The cryo-EM structures of SV1 and SV2 showed that they reorganized the extracellular halves of transmembrane helices 1, 6, and 7 and extracellular loops 2 and 3 to adopt a ligand-binding pocket-occupied conformation, thereby losing binding ability to the peptide. The results suggest a form of signal bias that is constitutive and ligand-independent, thus expanding our knowledge of biased signaling beyond pharmacological manipulation (i.e., ligand specific) as well as constitutive and ligand-independent (e.g., SV1 of the growth hormone-releasing hormone receptor).

Single hormone or synthetic agonist induces Gs/Gi coupling selectivity of EP receptors via distinct binding modes and propagating paths.

To accomplish concerted physiological reactions, nature has diversified functions of a single hormone at at least two primary levels: 1) Different receptors recognize the same hormone, and 2) different cellular effectors couple to the same hormone-receptor pair [R.P. Xiao, Sci STKE 2001, re15 (2001); L. Hein, J. D. Altman, B.K. Kobilka, Nature 402, 181-184 (1999); Y. Daaka, L. M. Luttrell, R. J. Lefkowitz, Nature 390, 88-91 (1997)]. Not only these questions lie in the heart of hormone actions and receptor signaling but also dissecting mechanisms underlying these questions could offer therapeutic routes for refractory diseases, such as kidney injury (KI) or X-linked nephrogenic diabetes insipidus (NDI). Here, we identified that Gs-biased signaling, but not Gi activation downstream of EP4, showed beneficial effects for both KI and NDI treatments. Notably, by solving Cryo-electron microscope (cryo-EM) structures of EP3-Gi, EP4-Gs, and EP4-Gi in complex with endogenous prostaglandin E2 (PGE2)or two synthetic agonists and comparing with PGE2-EP2-Gs structures, we found that unique primary sequences of prostaglandin E2 receptor (EP) receptors and distinct conformational states of the EP4 ligand pocket govern the Gs/Gi transducer coupling selectivity through different structural propagation paths, especially via TM6 and TM7, to generate selective cytoplasmic structural features. In particular, the orientation of the PGE2 ω-chain and two distinct pockets encompassing agonist L902688 of EP4 were differentiated by their Gs/Gi coupling ability. Further, we identified common and distinct features of cytoplasmic side of EP receptors for Gs/Gi coupling and provide a structural basis for selective and biased agonist design of EP4 with therapeutic potential.

Variation of TaMyb10 and their function on grain color and pre-harvest sprouting resistance of wheat.

Pre-harvest sprouting (PHS) is a significant threat to global food security due to its association with losses in both yield and quality. Among the genes involved in PHS resistance in wheat, PHS-3D (TaMyb10-D) plays a crucial role. Here, we characterized the sequence variations of TaMyb10 genes in 416 bread wheat and 302 Aegilops tauschii accessions. Within TaMyb10-A sequences, we identified a deletion ranging from 214 to 305 bp in the signal and amino acid coding region, present in 61.3% of the accessions. Similarly, 79.3% of the TaMyb10-B sequences within the third exon region exhibited a 19 bp deletion. Additionally, 40.8% of the accessions lacked the 2.4 Mb fragment (in/del mutations) on Chr3D, where TaMyb10-D/PHS-3D was located. Interestingly, the geographical distribution of accessions showed little correlation with the divergence of TaMyb10. TaMyb10-A-IIIDele, TaMyb10-B-IVDele, and TaMyb10-D-VDele genotypes were prevalent in wheat populations across continents. Despite their structural variations, the five distinct protein types exhibited comparable ability to bind the promoters of downstream genes in the flavonoid and ABA pathways, such as CHS, DFR, and NCED. Furthermore, the combination of TaMyb10 homologs was significantly associated with grain color and germination percentages. Accessions exclusively harboring TaMyb10-D displayed red seed color and reduced germination percentages, indicating the predominant role of TaMyb10-D compared to TaMyb10-A and TaMyb10-B. This comprehensive investigation enhances our understanding of the structural variations and functional divergence of TaMyb10, providing valuable insights and resources for improving PHS resistance in wheat.

Cryo-EM structures of cancer-specific helical and kinase domain mutations of PI3Kα.

Phosphoinositide 3-kinases (PI3Ks) are a family of lipid kinases that perform multiple and important cellular functions. The protein investigated here belongs to class IA of the PI3Ks; it is a dimer consisting of a catalytic subunit, p110α, and a regulatory subunit, p85α, and is referred to as PI3Kα. The catalytic subunit p110α is frequently mutated in cancer. The mutations induce a gain of function and constitute a driving force in cancer development. About 80% of these mutations lead to single-amino-acid substitutions in one of three sites of p110α: two in the helical domain of the protein (E542K and E545K) and one at the C-terminus of the kinase domain (H1047R). Here, we report the cryo-electron microscopy structures of these mutants in complex with the p110α-specific inhibitor BYL-719. The H1047R mutant rotates its sidechain to a new position and weakens the kα11 activation loop interaction, thereby reducing the inhibitory effect of p85α on p110α. E542K and E545K completely abolish the tight interaction between the helical domain of p110α and the N-terminal SH2 domain of p85α and lead to the disruption of all p85α binding and a dramatic increase in flexibility of the adaptor-binding domain (ABD) in p110α. Yet, the dimerization of PI3Kα is preserved through the ABD-p85α interaction. The local and global structural features induced by these mutations provide molecular insights into the activation of PI3Kα, deepen our understanding of the oncogenic mechanism of this important signaling molecule, and may facilitate the development of mutant-specific inhibitors.

Structural basis of peptidomimetic agonism revealed by small- molecule GLP-1R agonists Boc5 and WB4-24.

Glucagon-like peptide-1 receptor (GLP-1R) agonists are effective in treating type 2 diabetes and obesity with proven cardiovascular benefits. However, most of these agonists are peptides and require subcutaneous injection except for orally available semaglutide. Boc5 was identified as the first orthosteric nonpeptidic agonist of GLP-1R that mimics a broad spectrum of bioactivities of GLP-1 in vitro and in vivo. Here, we report the cryoelectron microscopy structures of Boc5 and its analog WB4-24 in complex with the human GLP-1R and Gs protein. Bound to the extracellular domain, extracellular loop 2, and transmembrane (TM) helices 1, 2, 3, and 7, one arm of both compounds was inserted deeply into the bottom of the orthosteric binding pocket that is usually accessible by peptidic agonists, thereby partially overlapping with the residues A8 to D15 in GLP-1. The other three arms, meanwhile, extended to the TM1-TM7, TM1-TM2, and TM2-TM3 clefts, showing an interaction feature substantially similar to the previously known small-molecule agonist LY3502970. Such a unique binding mode creates a distinct conformation that confers both peptidomimetic agonism and biased signaling induced by nonpeptidic modulators at GLP-1R. Further, the conformational difference between Boc5 and WB4-24, two closed related compounds, provides a structural framework for fine-tuning of pharmacological efficacy in the development of future small-molecule therapeutics targeting GLP-1R.

Nanobodies and chemical cross-links advance the structural and functional analysis of PI3Kα.

Nanobodies and chemical cross-linking were used to gain information on the identity and positions of flexible domains of PI3Kα. The application of chemical cross-linking mass spectrometry (CXMS) facilitated the identification of the p85 domains BH, cSH2, and SH3 as well as their docking positions on the PI3Kα catalytic core. Binding of individual nanobodies to PI3Kα induced activation or inhibition of enzyme activity and caused conformational changes that could be correlated with enzyme function. Binding of nanobody Nb3-126 to the BH domain of p85α substantially improved resolution for parts of the PI3Kα complex, and binding of nanobody Nb3-159 induced a conformation of PI3Kα that is distinct from known PI3Kα structures. The analysis of CXMS data also provided mechanistic insights into the molecular underpinning of the flexibility of PI3Kα.

Abnormal global alternative RNA splicing in COVID-19 patients.

Viral infections can alter host transcriptomes by manipulating host splicing machinery. Despite intensive transcriptomic studies on SARS-CoV-2, a systematic analysis of alternative splicing (AS) in severe COVID-19 patients remains largely elusive. Here we integrated proteomic and transcriptomic sequencing data to study AS changes in COVID-19 patients. We discovered that RNA splicing is among the major down-regulated proteomic signatures in COVID-19 patients. The transcriptome analysis showed that SARS-CoV-2 infection induces widespread dysregulation of transcript usage and expression, affecting blood coagulation, neutrophil activation, and cytokine production. Notably, CD74 and LRRFIP1 had increased skipping of an exon in COVID-19 patients that disrupts a functional domain, which correlated with reduced antiviral immunity. Furthermore, the dysregulation of transcripts was strongly correlated with clinical severity of COVID-19, and splice-variants may contribute to unexpected therapeutic activity. In summary, our data highlight that a better understanding of the AS landscape may aid in COVID-19 diagnosis and therapy.

Highly Luminescent Chiral Double π-Helical Nanoribbons.

Unveiling the mechanism behind chirality propagation and dissymmetry amplification at the molecular level is of significance for the development of chiral systems with comprehensively outstanding chiroptical performances. Herein, we have presented a straightforward Cu-mediated Ullmann homocoupling approach to synthesize perylene diimide-entwined double π-helical nanoribbons encompassing dimer, trimer, and tetramer while producing homochiral or heterochiral linking of chiral centers. A significant dissymmetry amplification was achieved, with absorption dissymmetry factors (|gabs|) increasing from 0.009 to 0.017 and further to 0.019, and luminescence dissymmetry factors (|glum|) rising from 0.007 to 0.013 and eventually to 0.015 for homochiral double π-helical oligomers. The disparity of magnetic transition dipole moment (m) densities in homochiral and heterochiral tetramers by time-dependent density functional theory calculations confirmed that homochiral oligomerization can maximize the total m, which is favorable for achieving ever-increasing g factors. Notably, these double π-helices exhibited exceptional photoluminescence quantum yields (ΦPL) ranging from 83 to 95%. The circularly polarized luminescence brightness (BCPL) eventually reached a remarkable 575 M-1 cm-1 for the homochiral tetramer, which is among the highest values reported for chiral small molecules. This kind of linearly extended double π-helices offers a platform for a comprehensive understanding of the mechanism behind chirality propagation and dissymmetry amplification.

Regulated stepwise ESDPT mechanism associated with chalcogen substitutions in BDIBD derivatives.

Inspired by the brilliant photochemical and photophysical properties of organic molecules containing chalcogenide substitutions that could be potentially applied across various disciplines, in this work, the effects of the atomic electronegativity of chalcogens (O, S, and Se) on hydrogen bond interactions and excited state proton transfer (ESPT) are mainly focused. We present characteristic oxygen-hydroxybenzazole-substituted 2,5-bis(4,5-diphenyl-1H-imidazol-2-yl)benzene-1,4-diol (BDIBD) derivatives that contain intramolecular double hydrogen bonds. The main objective of this study was to explore in detail the influence of the change of chalcogen atomic electronegativity on dual hydrogen bond interaction and ESPT behavior. By comparing the structural changes and infrared (IR) vibrational spectra of BDIBD derivative (BDIBD-O, BDIBD-S and BDIBD-Se) fluorophores in S0 and S1 states, combined with the preliminary detection of hydrogen bond interaction via the core-valence bifurcation (CVB) index and predicted hydrogen bonding energy (EHB), we conclude that dual hydrogen bonds should be strengthened in the S1 state, which is favorable for the occurrence of ESPT reactions. The charge recombination behavior of hydrogen bonds, induced by photoexcitation, further illustrates this point. By constructing potential energy surfaces (PESs) based on restrictive optimization and by searching the transition state (TS) structure, we finally elucidate stepwise excited-state double proton transfer (ESDPT). Specifically, we confirm that a change in atomic electronegativity has a regulatory effect on the ESDPT behavior in BDIBD derivatives, that is, lower atomic electronegativity is more conducive to stepwise ESDPT.