Atomic charge is crucial in drug design for analyzing reactive sites and interactions between ligands and targets. While quantum mechanical methods offer high accuracy, they are generally computationally costly. Conversely, empirical approaches, while computationally efficient, frequently suffer from lack of precision and generalizability. Recent a number of machine learning-based models have been developed for atomic charge predictions, but they struggle with accurately representing molecular structures and capturing the chemical environments affecting atomic charges, thus limiting their generalization and accuracy. To overcome these limitations, we propose LumiCharge, a novel atomic charge prediction framework that incorporates high-order spherical harmonics convolutions and explicitly models multibody interactions. In constructing this model, we employ a strategy that integrates both high- and low-order information, enhancing its geometric spatial perception capability, which is currently underexplored in the field. Benchmark evaluations demonstrate that LumiCharge outperforms state-of-the-art (SOTA) models by 30%-60% across diverse data sets. Additionally, in cross-scale experiments, LumiCharge demonstrates exceptional extrapolation capability and robustness across molecules of varying sizes, effectively overcoming the limitations imposed by molecular sizes. On an external halogen-containing test set, LumiCharge achieves an RMSE of 0.055e, meeting practical application requirements. Finally, a case study of virtual screening for the androgen receptor (AR) target further validates its outstanding accuracy compared to the OPLS3e force field and other deep learning (DL)-based baseline models, highlighting its exceptional generalization capacity and practical utility in real-world scenarios.
Qun Su, Hui Zhang, Qiaolin Gou et al.· Journal of Physical Chemistr...· 2 citations
Atomic charge is a fundamental quantum chemical property essential for advancing drug design and discovery. Although quantum mechanics (QM) methods offer the highest level of accuracy, their computational demands scale quadratically with the number of atoms, limiting their practicality for large-scale applications. In light of this, empirical and semiempirical methods have been introduced to improve computational efficiency, albeit often at the expense of accuracy. The advent of artificial intelligence has witnessed a growing application of machine learning (ML) techniques to accelerate atomic charge predictions. However, existing ML models often suffer from low accuracy and limited generalization capabilities. To address these challenges, we introduce an advanced equivariant graph attention neural network specifically engineered to model long-range atomic electrostatic interactions with high precision. This model introduces a sophisticated global graph attention mechanism, enabling it to capture charge contributions across multiple scales. By utilizing a combination of structural symmetry-preserving transformations and multiscale attention, our approach not only preserves the inherent symmetries of molecular structures but also substantially improves the model's accuracy, generalization, and robustness in complex scenarios. Our empirical analyses demonstrate that, compared to leading baseline models, the proposed model improves charge prediction accuracy by over 40% on average across various charge-calculation schemes. Remarkably, the model achieves superior performance on the external RESP (restrained electrostatic potential) test data sets, with a 54.6% improvement over the baseline. Additionally, we evaluated our charge model under the setting of virtual screening, where it outperforms both the OPLS3 charges and baseline deep learning models across all evaluation metrics, highlighting its extensive potential for scientific discovery.
Qiaolin Gou, Qun Su, Jike Wang et al.· Journal of Chemical Informat...· 1 citation
To fulfill functions for differentially regulating the downstream signaling pathways, functional ligands (i.e., agonists or antagonists) targeting nuclear receptors (NRs) are designed to stabilize different conformations (active or inactive) of the proteins. However, in practical applications, it is usually difficult to determine the molecular category of an NR ligand because these molecules all bind in the same location of an NR protein, namely, the ligand-binding pocket (LBP). Considering that ligands with different properties (agonists or antagonists) prefer to bind with differential conformations of NRs, it is possible to identify the molecular type of a given ligand through the differential binding environment (active or inactive conformations) of the protein-ligand interaction. Therefore, in this study, we established a unified model (NRIGN) based on the deep graphic architecture to discriminate agonists and antagonists targeting 26 successful or in-clinical-trial NR targets. Our result shows that NRIGN achieves an excellent prediction accuracy (ACC >0.95) and is robust enough to be applied in various real-world scenarios, such as predicting the molecular type of ligands in crystallized NR structures, ligands with multiple NR activities, and ligands with their types altered by target mutations. The proposed model is expected to promote rational design of drugs targeting NR proteins.
Kaimo Yang, Dejun Jiang, Qirui Deng et al.· Journal of Chemical Informat...· 2 citations
Colorectal cancer (CRC) exhibits substantial molecular heterogeneity, necessitating the inference of subtype-specific driver genes and their interactions for drug-target discovery and precision oncology. Prior studies often fail to capture subtle, latent nonlinear regulatory mechanisms (dark causal relationships) driving tumor progression in specific subtypes. Here, we develop an explainable intelligence computational framework, Symbolic Trajectory-Embedded Dark Causal Interaction Inference (STE-DC2I), which combines symbolic trajectory embedding with historical prediction mechanisms to model nonmonotonic oscillatory dependencies between genes. Integrating single-cell transcriptomic and multiomics profiles from malignant epithelial subpopulations, STE-DC2I classifies CRC subtypes, reconstructs developmental trajectories, and uncovers interpretable subtype-specific driver genes with functional relevance. Unlike correlation-based and explicit causal approaches, STE-DC2I captures weak yet biologically critical regulatory signals, outperforming state-of-the-art methods in predicting subtype-specific CRC driver genes. Functional assays in CRC cell lines (in vitro) validated nine predicted driver genes, highlighting their therapeutic potential.This work systematically explores dark causal interactions between genes in CRC subtypes. STE-DC2I offers interpretable insights and a generalizable strategy for CRC drug-target discovery.
Meng Huang, Huijin Hu, Ming Li et al.· Journal of Chemical Informat...· 0 citations
It is evidenced that many elaborately designed molecules that can interact well with the binding pocket of their target fail to exhibit activity in wet-lab experiments. This may associate with the interacting process of drug-target recognition. To efficiently characterize the drug-target interacting process, various enhanced sampling technologies have been proposed; yet, very few studies have systemically investigated whether the settings of these simulations are favorable to characterize the purposed tasks. Here, by comparing two popular enhanced sampling technologies, namely, the well-temped metadynamics and random acceleration molecular dynamics (RAMD), we systemically investigate the strategies to efficiently characterize the dissociating process of protein-ligand interactions. Two target families are employed for the analysis, including the kinase family (represented by TRK1) that represents the interaction-pathway obvious systems and the nuclear receptor family (represented by THRβ) that represents the interaction-pathway unobvious systems. Our results suggest that (1) in terms of maintaining stability of the protein structure, MetaD at various simulation conditions and RAMD with a large random force are good choice; (2) drug residence time derived from both MetaD and RAMD based on various parameters shows reasonable correlation to the experimental binding strength of the ligands, but RAMD usually runs with much less simulation time; and (3) both enhanced sampling methods result in reasonably consistent pathway preference for the two target families. Taken together, it will be much time-saving to utilize RAMD with high random force for interaction pathway exploration for both the pathway obvious and unobvious systems if the protein keeps stable in the simulation; otherwise, MetaD with a high bias factor is proposed to balance the computational accuracy and efficiency for the exploration.
Zhiliang Jiang, Mingyun Shen, Zhe Wang et al.· Journal of Chemical Physics· 1 citation
Molecular glues, including protein degraders and protein-protein interaction (PPI) stabilizers, have emerged as a new paradigm of drug design for regulating interactions between biomacromolecules; yet it is still a challenge for rational design of molecular glues. KRAS, as a prevalent oncogenic driver, is notoriously difficult to target by traditional small molecular drugs due to its challenging binding surface and frequent mutations. Although the small molecular drug RMC7977 has been designed as a PPI stabilizer for stabilizing the inherently weak RAS-CYPA interaction, the precise molecular mechanism underlying its stabilization effect and selectivity difference requires a deeper understanding. To this end, we leverage an integrated computational strategy combining molecular dynamics (MD) simulation, end-point binding free-energy calculation, and enhanced sampling technologies to elucidate the dynamic characteristics of RAS-ligand-CYPA interactions. Our result exhibits a high correlation between the predicted binding affinities and the experimental observations, demonstrating that RMC7977, acting as a strong PPI stabilizer, significantly enhances the stability of the KRAS-CYPA interaction, where, by delicately remodeling the protein-protein interface, the drug optimizes various interactions. Moreover, the results also uncover the dynamic process of stabilizer-mediated KRAS-CYPA stabilization and the mechanistic origin of the binding selectivity. This study provides essential molecular-level insights into RMC7977's function and offers a valuable computational framework for evaluating the stabilization effect of ligands targeting the KRAS-CYPA and other challenging PPI systems.
Kexin Xu, Mingyun Shen, Zhe Wang et al.· Journal of Chemical Informat...· 0 citations
Proteolysis-targeting chimeras (PROTACs) present a potentially effective strategy against various diseases via selective proteolysis. How to increase the efficacy of PROTACs remains challenging. Here, we explore the necessity of the linker, which has been deemed as an integral part of heterobifunctional PROTACs. Adopting single amino acid-based degradation signals, we find that the linker is not a required feature of the PROTACs. Notably, the linker-free PROTAC, Pro-BA, exhibits superior efficacy over its linker-bearing counterparts in degrading EML4-ALK and inhibiting lung cancer cell growth, as Pro-BA induces a stronger interaction between the target and the E3 ubiquitin ligase. Pro-BA is a water-soluble, orally administered degrader that significantly inhibits the tumor growth in a xenograft mouse model. The broad applicability of this linker-free PROTAC strategy is further validated through the development of BCR-ABL degrader. Our study introduces a design paradigm for PROTACs, potentially facilitating the advancement of more efficient therapeutic degraders. Linkers are traditionally seen as important for PROTAC activity. Here, the authors demonstrate that linker-free PROTACs can outperform traditional designs, marking a paradigm shift in PROTAC development for targeted protein degradation.
The discovery of CAP-Gly domain-containing linker protein 1(CLIP1)-Leukocyte tyrosine kinase (LTK) as an oncogenic fusion reveals a unique dependency not only on LTK kinase activity but also on CLIP1-mediated multimerization, a noncatalytic function that drives oncogenic signaling. While this fusion is currently targeted with anaplastic lymphoma kinase inhibitors, their exclusive focus on kinase inhibition leaves the scaffolding function intact, necessitating a complete protein clearance strategy. Here, we report the AI-guided development of a first-in-class proteolysis-targeting chimera (PROTAC) designed to selectively degrade the CLIP1-LTK fusion protein. By integrating deep learning models for ternary complex prediction with structure-based molecular optimization, we designed DCL05, an orally bioavailable degrader of CLIP1-LTK fusion protein, achieving picomolar degradation potency (DC50 = 40 pM) and robust antitumor activity. DCL05 consistently outperformed existing kinase inhibitors across a broad spectrum of LTK resistance-associated mutations, both in vitro and in vivo. Collectively, our study explores resistance-associated contexts of LTK and establishes a structure-guided PROTAC development pipeline, providing a promising therapeutic strategy for overcoming acquired resistance in kinase-driven cancers.
Shicheng Chen, Haiting Duan, S. Zhong et al.· Proceedings of the National...· 0 citations
Multi-target drugs hold great promise for treating complex diseases, yet existing methodologies predominantly rely on ligand-based approaches, which lack sufficient biological context and are often confined to specific target pairs, resulting in limited generalizability. Here, we introduce LaMGen, a general-purpose multi-target drug design framework powered by large language models (LLMs). Built on MTD2025, a dataset comprising over 600,000 quantum-accurate molecular conformations and 700,000 multi-target associations, LaMGen directly yields energy-favorable conformations with quantum-level accuracy. The framework integrates ESM-C protein embeddings, rotation-aware ligand tokens, and a TriCoupleAttention module to capture multi-level target–ligand interactions. Across independent benchmarks, LaMGen outperforms diffusion-based model across multiple properties, generating molecules in an average of 0.44 s, while preserving high conformational plausibility. Retrospective analyses demonstrate that LaMGen not only can reproduce molecules identical to known actives, but also consistently produces structurally novel candidates with conserved core scaffolds and superior binding affinities. Designing effective multi-target therapeutics remains a major challenge, as existing ligand- or protein-centric methods struggle to generate biologically contextualized, spatially valid 3D molecules, particularly for triple-target systems. This study introduces LaMGen, an LLM-powered framework that leverages large-scale protein-ligand data and rotation-aware molecular encoding to rapidly produce chemically plausible multi-target candidates, achieving strong zero-shot generalization, superior molecular quality, and robust performance across dual- and triple-target design tasks.
Qun Su, Qiaolin Gou, Hui Zhang et al.· Nature Communications· 1 citation