The experiment used here shows a method of molecular docking combined with cellular thermal shift assay to predict and validate the interaction between small molecules and protein targets.
Proteins are fundamental to human physiology, with their targets being crucial in research and drug development. The identification and validation of crucial protein targets have become integral to drug development. Molecular docking is a computational tool widely utilized to investigate protein-ligand binding, especially in the context of drug and protein target interactions. For the experimental verification of the binding and to access the binding of the drug and its target directly, the cellular thermal shift assay (CETSA) method is used. This study aimed to integrate molecular docking with CETSA to predict and validate interactions between drugs and vital protein targets. Specifically, we predicted the interaction between xanthatin and Keap1 protein as well as its binding mode through molecular docking analysis, followed by verification of the interaction using the CETSA assay. Our results demonstrated that xanthatin could establish hydrogen bonds with specific amino acid residues of Keap1 protein and reduce the thermostability of Keap1 protein, indicating that xanthatin could directly interact with Keap1 protein.
Proteins are highly important macromolecules in living organisms and possess a diverse range of unique functions within cells, such as membrane composition, cytoskeleton formation, enzyme activity, transportation, cell signaling, and involvement in both intracellular and extracellular mechanisms1,2,3. Proteins manifest their biological functions primarily through specific interactions with a variety of molecules, including other proteins, nucleic acids, small molecule ligands, and metal ions1,4. Ligands are small molecular compounds that specifically bind to proteins in an organism. The interaction between proteins and ligands occurs at specific sites on the protein, called the binding sites, also known as the binding pockets5. In medicinal chemistry research, the focus lies in identifying key proteins that are clearly associated with diseases, which serve as targets for drugs6. Therefore, gaining a deep understanding of the binding sites between proteins and ligands is of utmost importance in promoting drug discovery, design, and research7,8.
Molecular docking is a widely used computational tool for studying protein-ligand binding, which employs the three-dimensional structures of proteins and ligands to explore their primary binding modes and affinities when forming stable complexes9,10,11. The application of molecular docking technology originated in the 1970s. Based on the lock and key pairing principle and utilizing the algorithms of molecular docking software, one can determine the interaction between compounds and molecular targets by analyzing docking results. This approach enables the prediction of active binding sites for both the compound and the target molecule. Consequently, it facilitates the identification of an optimal binding conformation (here called the binding model) for ligand-receptor interactions, which is crucial for understanding the mechanics of these molecular engagements12,13,14,15. While molecular docking provides valuable computer-based predictions of ligand-receptor interactions, it is important to note that these are preliminary findings. Consequently, further experimental verification is essential to confirm these interactions.
The cellular thermal shift assay (CETSA), initially proposed by Pär Nordlund's research team in 2013, serves as a method for validating drug-target protein interactions. This technique specifically tests the thermal stability of target proteins induced by drug binding, providing a practical approach for confirming molecular interactions16,17,18. This approach is based on the fundamental principle that ligand binding initiates a thermal shift within target proteins and is applicable to a wide array of biological samples, including cell lysates, intact living cells, and tissues19,20. CETSA supports direct target engagement of small molecules in intact cells by detecting thermodynamic stabilization of proteins due to ligand binding and linking the observed phenotypic response to the target compound21,22. Among the various methodologies derived from CETSA, Western Blot-CETSA (WB-CETSA) is considered a classical approach. Following sample preparation using the CETSA method, western blot analysis is utilized to detect alterations in the thermal stability of the target protein. This allows for the precise determination of drug-protein interactions within cellular systems17,23.
Xanthatin is a bioactive compound isolated from the plant Xanthium L. with properties such as anti-inflammatory, which has been used in traditional Chinese medicine to treat diseases like nasal sinusitis and arthritis24,25. The kelch-like ECH-associated protein 1 (Keap1) is a component of the Cullin3-based Cullin-RING E3 ubiquitin ligase multi-subunit protein complex and an important regulator of intracellular redox homeostasis, which influences the intensity and duration of the inflammatory response by modulating the intracellular redox state26. In this study, we first utilized molecular docking to investigate the interaction between xanthatin (small molecule) and Keap1 protein, aiming to predict their binding mode. Subsequently, we employed the CETSA method to validate this interaction by assessing the impact of xanthatin on the thermal stability of the Keap1 protein.
1. Downloading the structures of xanthatin and Keap1
2. Molecular docking
3. Cell culture and CETSA sample preparation
4. Western blot
Molecular docking analysis predicted the interaction between xanthatin and Keap1 protein. Figure 2 demonstrates the formation of hydrogen bonds between xanthatin and amino acid residues Gly-367 and Val-606 of Keap1 protein, with a hydrogen bond length of 2.17 Å for Gly-367 and 2.13 Å for Val-606. In addition, the calculated docking score of -5.69 kcal/mol signifies a good binding affinity between xanthatin and Keap1 protein.
The CETSA method showed that xanthatin binding shifted the thermal stability of Keap1 protein, confirming the interaction between xanthatin and Keap1 protein21. Figure 3 indicates that xanthatin remarkably shifts the thermal stability of Keap1 protein within the temperature range of 48 °C to 57 °C, as compared to the DMSO group. To ensure equal sample loading, firstly, the density of cell inoculation in which we performed the experiments was consistent; secondly, the cell samples from the DMSO group and the xanthatin group were mixed well in the cell culture dishes and then were equally distributed into two microcentrifuge tubes, which ensured the consistent amount of cell samples, and, the volume of supernatant added to each PCR tube was 60 µL, and lastly, the volume of each sample uploaded in the Western blot was 20 µL, which ensured equal loading for heated samples. Protein degradation in samples from around 51 °C is mainly due to changes in protein structure and accelerated chemical reactions caused by high temperatures. The above results demonstrate that xanthatin could directly bind to Keap1 protein.
Figure 1: Loading order of western blot. The temperature in order from left to right were 45 °C, 48 °C, 51 °C, 54 °C, 57 °C, 60 °C, and 63 °C. The first lane was added with a marker; two lanes were used for each temperature; the first lane was the DMSO group, and the second lane was the xanthatin group. Please click here to view a larger version of this figure.
Figure 2: Interaction between xanthatin and Keap1 (PDB: 2FLU) protein. (A) The 2D representation of xanthatin binding to Keap1. (B) The 3D visualization of xanthatin and Keap1 protein. Please click here to view a larger version of this figure.
Figure 3: Effect of xanthatin on the thermal stability of Keap1 protein. (A) The effects of xanthatin on the thermal stability of Keap1 were detected by CETSA. (B) The optical densities of Keap1 were normalized to those obtained at 45 °C. Data are represented as mean ± SD. Statistical analysis was performed by one-way analysis of variance (one-way ANOVA). n = 3. *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001. Please click here to view a larger version of this figure.
The identification of disease targets and the discovery and development of drugs are closely interconnected27. By precisely targeting specific targets, drug candidates can be developed to treat particular diseases more effectively while concurrently minimizing the side effects associated with the drugs28,29. The most commonly used targets are protein targets30. However, the identification of special protein targets poses an immense challenge due to the extensive diversity of proteins within cells30,31. In this paper, we adopted an approach that combines molecular docking with CETSA to identify and validate protein targets. The results revealed a direct interaction between xanthatin and Keap1 protein. The integration of computer technology and biological experiments can also be employed in the investigation of other drugs and protein targets to ascertain the presence of a direct interaction between the drug and the target. This approach, which combines prediction with verification, effectively reduces both time consumption and economic costs. Furthermore, the experimental instruments and materials are easy to obtain, and the operation process is not complicated. Finally, The CETSA samples are easy to prepare and test by the western blot method, which ensures the reproducibility of the samples.
In our experience, there are several critical aspects that require special attention during experimental procedures. Firstly, it is essential to select the crystal structure of the protein resolved by X-ray analysis when choosing the protein structure before molecular docking. Secondly, cell density and freezing/thawing time are crucial factors in CETSA experiments. Cells in the logarithmic growth phase should be selected for CETSA experiments as their growth state is stable and has minimal impact on experimental results. During repeated freezing and thawing of cells, timing should be well-controlled with liquid nitrogen used for freezing until a white solid form, followed by immediate thawing at room temperature before subsequent freezing cycles. Lastly, western blot experiments must be completed promptly after preparing CETSA samples to prevent protein degradation.
Molecular docking and CETSA are both well-established techniques that play a crucial role in the field of drug development. Molecular docking is instrumental in predicting molecular-level interactions between compounds and their targets, offering valuable insights into potential binding modes10,11. CETSA, in contrast, focuses on assessing the impact of drugs on protein thermal stability and serves as a tool for validating drug-protein interactions17. While CETSA is a relatively straightforward method for verification, it is also a low throughput way by using western blot for the study of thermal shift of drug-protein complex, and it should not be the sole technique relied upon. In addition, the steps of CETSA are complicated, and any step errors during the operation can have a great impact on the results. Furthermore, the CETSA method cannot determine the specific binding site and binding mode of the compound to the protein target and cannot indicate whether the binding is covalent or non-covalent17,18,23. Beyond thermal shift assay, other kinds of energetic assay, such as pH-dependent protein precipitation (pHDPP), can also be used for the study of ligand-protein interactions32. Moreover, thermal proteome profiling is an energy-based method for revealing ligand-protein interactions, allowing high-throughput, large-scale analysis of protein-drug interactions by combining proteomics with CETSA33.
In addition, localized surface plasmon resonance (LSPR) can provide further qualitative and quantitative analysis34,35. It not only determines the existence of interactions between a drug and its target but also measures the affinity parameters and kinetic parameters of these intermolecular interactions, offering a more comprehensive understanding.
Beyond thermal shift assay, pH-dependent protein precipitation (pHDPP) can also be used for the study of ligand-protein interactions32. Molecular docking is a physico-mechanics-based simulation method for predicting interactions between small molecules and biomolecules. This method aims to achieve binding by optimizing the conformation and orientation of small molecules to produce the best complementarity between them and biomolecules. Target prediction utilizes artificial intelligence techniques to predict the binding sites between small molecules and biomolecules. The two approaches have key differences in method application and prediction accuracy. Method application: molecular docking is mainly used for drug design and optimization, providing a theoretical basis for new drug discovery. Target prediction, with the help of more advanced artificial intelligence, is more often applied to high-throughput screening and rapid discovery of lead compounds to improve the efficiency of drug discovery. Prediction accuracy: Since molecular docking is based on physical simulations, its prediction accuracy is higher, but the computational cost is also relatively high. While target prediction has a lower computational cost, the prediction accuracy depends on the abundance of training data and model selection.
The authors have nothing to disclose.
This work was supported by National Natural Science Foundation of China (82004031) and Sichuan Science and Technology Program (2022NSFSC1303). We express our great appreciation to Jiayi Sun at Innovative Institute of Chinese Medicine and Pharmacy, Chengdu University of Traditional Chinese Medicine, for the assistance with western blot.
0.45 μm Polyvinylidene fluoride membrane | Millipore | PR05509 | |
Anhydrous ethanol | Chron chemicals | 64-17-5 | |
Bovine serum albumin | BioFroxx | 4240GR100 | |
Broad-spectrum protease inhibitor mixtures | Boster Biological Technology Co., Ltd | AR1193 | |
DMSO | Boster Biological Technology Co., Ltd | PYG0040 | |
Enhanced chemiluminescence reagent | Beyotime Biotechnology Co., Ltd | P0018S | |
GAPDH antibody | ProteinTech Group Co., Ltd | 10494-1-AP | |
Gel Imaging Instrument | E-BLOT | Touch Imager Pro | |
Gradient PCR instrument | Biometra TADVANCED | Biometra Tadvanced 96SG | |
High-speed freezing centrifuge | Beckman Coulter | Allegra X-30R | |
Horseradish peroxidase-conjugated affiniPure goat antibody | ProteinTech Group Co., Ltd | SA00001-2 | |
Isopropyl alcohol | Chron chemicals | 67-63-0 | |
Keap1 antibody | Zen BioScience Co., Ltd | R26935 | |
Metal bath | Analytik Jena | TSC | |
Methanol | Chron chemicals | 67-56-1 | |
Ncmblot rapid transfer buffer (20×) | NCM Biotech Co., Ltd | WB4600 | |
Omni-Easy OneStep PAGE gel fast preparation kie | Epizyme Biotech Co., Ltd | PG212 | |
Phosphate buffer saline | Boster Biological Technology Co., Ltd | PYG0021 | |
Prestained Color Protein Marker | Biosharp | BL741A | |
Protein Blotting Electrophoresis System | Bio-Rad | MiniPROTEANÒTetra Cell | |
RAW264.7 cell | Beyotime Biotechnology Co., Ltd | C7505 | |
RAW264.7 cell-specific medium | Procell Life Science&Technology Co., Ltd | CM-0597 | |
SDS-PAGE protein loading buffer | Boster Biological Technology Co., Ltd | AR1112-10 | |
SDS-PAGE running buffer powder | Servicebio | G2018 | |
Tris buffered saline powder | Servicebio | G0001 | |
Tween 20 | BioFroxx | 1247ML100 | |
Water bath | Memmert | WNE10 | |
Water purifier | Millipore | Milli- IQ 7005 | |
Xanthatin | ChemConst Biotechnology Co., Ltd | CONST210706 |