Effect of in vivo post-translational modifications of the HMGB1 protein upon binding to platinated DNA: a molecular simulation study

High performance-oriented computer aided drug design approaches in the exascale era

INTRODUCTION

In 2023, the first exascale supercomputer was opened to the public in the US. With a demonstrated 1.1 exaflops of performance, Frontier represents an unprecedented breakthrough in high-performance computing (HPC). Currently, more (and more powerful) machines are being installed worldwide. Computer-aided drug design (CADD) is one of the fields of computational science that can greatly benefit from exascale computing for the benefit of the whole society. However, scaling CADD approaches to exploit exascale machines require new algorithmic and software solutions.

AREAS COVERED

Here, the authors consider physics-based and machine learning (ML)-aided techniques for the design of small molecule binders capable of leveraging modern parallel computer architectures. Specifically, the authors focus on HPC-oriented large-scale applications from the past 3 years that were enabled by (pre)exascale supercomputers by running on up tothousands of accelerated nodes.

EXPERT OPINION

In the area of ML, exascale computers can enable the training of generative models with unprecedented predictive power to design novel ligands, provided large amounts of high-quality data are available. Exascale computers could also unlock the potential of accurate ML-aided physics-based methods to boost the success rate of structure-based drug design campaigns. Currently, however, methodological developments are still required to allow routine large-scale applications of such rigorous approaches.

Click Chemistry-Based Force Spectroscopy Revealed Enhanced Binding Dynamics of Phosphorylated HMGB1 to Cisplatin-DNA

Cisplatin, a cornerstone in cancer chemotherapy, is known for its DNA-binding capacity and forms lesions that lead to cancer cell death. However, the repair of these lesions compromises cisplatin's effectiveness. This study investigates how phosphorylation of HMGB1, a nuclear protein, modifies its binding to cisplatin-modified DNA (CP-DNA) and thus protects it from repair. Despite numerous methods for detecting protein-DNA interactions, quantitative approaches for understanding their molecular mechanism remain limited. Here, we applied click chemistry-based single-molecule force spectroscopy, achieving high-precision quantification of the interaction between phosphorylated HMGB1 and CP-DNA. This method utilizes a synergy of click chemistry and enzymatic ligation for precise DNA-protein immobilization and interaction in the system. Our results revealed that HMGB1 binds to CP-DNA with a significantly high rupture force of ∼130 pN, stronger than most natural DNA-protein interactions and varying across different DNA sequences. Moreover, Ser14 is identified as the key phosphorylation site, enhancing the interaction's kinetic stability by 35-fold. This increase in stability is attributed to additional hydrogen bonding suggested by molecular dynamics (MD) simulations. Our findings not only reveal the important role of phosphorylated HMGB1 in potentially improving cisplatin's therapeutic efficacy but also provide a precise method for quantifying protein-DNA interactions.

Cisplatin Dependent Secretion of Immunomodulatory High Mobility Group Box 1 (HMGB1) Protein from Lung Cancer Cells

High mobility group box 1 (HMGB1) is secreted from activated immune cells, necrotic cells, and certain cancers. Previous studies have reported that different patterns of post-translational modification, particularly acetylation and oxidation, mediate HMGB1 release and confer distinct extracellular HMGB1 signaling activity. Here we report that cisplatin but not carboplatin induces secretion of HMGB1 from human A549 non-small cell lung cancer (NSCLC) cells. Cisplatin-mediated HMGB1 secretion was dose-dependent and was regulated by nuclear exportin 1 (XPO1) also known as chromosomal maintenance 1 (CRM1) rather than adenosine diphosphate (ADP)-ribosylation, acetylation, or oxidation. HMGB1, as well as lysine acetylation and cysteine disulfide oxidation of secreted HMGB1, were monitored by sensitive and specific assays using immunoprecipitation, stable isotope dilution, differential alkylation, and nano liquid chromatography parallel reaction monitoring/high-resolution mass spectrometry (nano-LC-PRM/HRMS). A major fraction of the HMGB1 secreted by low-dose cisplatin treatment of A549 NSCLC cells was found to be in the fully reduced form. In contrast, mainly oxidized forms of HMGB1 were secreted by dimethyl sulfoxide (DMSO)-mediated apoptosis. These findings suggest that inhibition of XPO1 could potentiate the anti-tumor activity of cisplatin by increasing the nuclear accumulation of HMGB1 protein, an inhibitor of cisplatin DNA-adduct repair. Furthermore, low-dose cisplatin therapy could modulate the immune response in NSCLC through the established chemokine activity of extracellular reduced HMGB1. This could potentially enhance the efficacy of subsequent immunotherapy treatment.

Thirty years of molecular dynamics simulations on posttranslational modifications of proteins

Posttranslational modifications (PTMs) are an integral component to how cells respond to perturbation. While experimental advances have enabled improved PTM identification capabilities, the same throughput for characterizing how structural changes caused by PTMs equate to altered physiological function has not been maintained. In this Perspective, we cover the history of computational modeling and molecular dynamics simulations which have characterized the structural implications of PTMs. We distinguish results from different molecular dynamics studies based upon the timescales simulated and analysis approaches used for PTM characterization. Lastly, we offer insights into how opportunities for modern research efforts on in silico PTM characterization may proceed given current state-of-the-art computing capabilities and methodological advancements.

Conformational Changes of Glutamine 5'-Phosphoribosylpyrophosphate Amidotransferase for Two Substrates Analogue Binding: Insight from Conventional Molecular Dynamics and Accelerated Molecular Dynamics Simulations

Glutamine 5′-phosphoribosylpyrophosphate amidotransferase (GPATase) catalyzes the synthesis of phosphoribosylamine, pyrophosphate, and glutamate from phosphoribosylpyrophosphate, as well as glutamine at two sites (i.e., glutaminase and phosphoribosylpyrophosphate sites), through a 20 Å NH₃ channel. In this study, conventional molecular dynamics (cMD) simulations and enhanced sampling accelerated molecular dynamics (aMD) simulations were integrated to characterize the mechanism for coordination catalysis at two separate active sites in the enzyme. Results of cMD simulations illustrated the mechanism by which two substrate analogues, namely, DON and cPRPP, affect the structural stability of GPATase from the perspective of dynamic behavior. aMD simulations obtained several key findings. First, a comparison of protein conformational changes in the complexes of GPATase–DON and GPATase–DON–cPRPP showed that binding cPRPP to the PRTase flexible loop (K326 to L350) substantially effected the formation of the R73-DON salt bridge. Moreover, only the PRTase flexible loop in the GPATase–DON–cPRPP complex could remain closed and had sufficient space for cPRPP binding, indicating that binding of DON to the glutamine loop had an impact on the PRTase flexible loop. Finally, both DON and cPRPP tightly bonded to the two domains, thereby inducing the glutamine loop and the PRTase flexible loop to move close to each other. This movement facilitated the transfer of NH3 via the NH3 channel. These theoretical results are useful to the ongoing research on efficient inhibitors related to GPATase.

The Potential Therapeutic Role of the HMGB1-TLR Pathway in Epilepsy

Epilepsy is one of the most common serious neurological disorders, affecting over 70 million people worldwide. For the treatment of epilepsy, antiepileptic drugs (AEDs) and surgeries are widely used. However, drug resistance and adverse effects indicate the need to develop targeted AEDs based on further exploration of the epileptogenic mechanism. Currently, many efforts have been made to elucidate the neuroinflammation theory in epileptogenesis, which may show potential in the treatment of epilepsy. In this respect, an important target protein, high mobility group box 1 (HMGB1), has received increased attention and has been developed rapidly. HMGB1 is expressed in various eukaryotic cells and localized in the cell nucleus. When HMGB1 is released by injuries or diseases, it participates in inflammation. Recent studies suggest that HMGB1 via Toll-like receptor (TLR) pathways can trigger inflammatory responses and play an important role in epilepsy. In addition, studies of HMGB1 have shown its potential in the treatment of epilepsy. Herein, the authors analyzed the experimental and clinical evidence of the HMGB1-TLR pathway in epilepsy to summarize the theory of epileptogenesis and provide insights into antiepileptic therapy in this novel field.

All-atom simulations to studying metallodrugs/target interactions

Metallodrugs are extensively used to treat and diagnose distinct disease types. The unique physical-chemical properties of metal ions offer tantalizing opportunities to tailor effective scaffolds for selectively targeting specific biomolecules. Modern experimental techniques have collected a large body of structural data concerning the interactions of metallodrugs with their biomolecular targets, although being unable to exhaustively assess the molecular basis of their mechanism of action. In this scenario, the complementary use of accurate computational methods allows uncovering the minutiae of metallodrugs/targets interactions and their underlying mechanism of action at an atomic-level of detail. This knowledge is increasingly perceived as an invaluable requirement to rationally devise novel and selective metallodrugs. Building on literature studies, selected largely from the last 2 years, this compendium encompasses a cross-section of the current role, advances, and challenges met by computer simulations to decipher the mechanistic intricacies of prototypical metallodrugs.

Use of Molecular Dynamics Simulations in Structure-Based Drug Discovery

BACKGROUND

Traditional drug discovery is a lengthy process which involves a huge amount of resources. Modern-day drug discovers various multidisciplinary approaches amongst which, computational ligand and structure-based drug designing methods contribute significantly. Structure-based drug designing techniques require the knowledge of structural information of drug target and drug-target complexes. Proper understanding of drug-target binding requires the flexibility of both ligand and receptor to be incorporated. Molecular docking refers to the static picture of the drug-target complex(es). Molecular dynamics, on the other hand, introduces flexibility to understand the drug binding process.

OBJECTIVE

The aim of the present study is to provide a systematic review on the usage of molecular dynamics simulations to aid the process of structure-based drug design.

METHOD

This review discussed findings from various research articles and review papers on the use of molecular dynamics in drug discovery. All efforts highlight the practical grounds for which molecular dynamics simulations are used in drug designing program. In summary, various aspects of the use of molecular dynamics simulations that underline the basis of studying drug-target complexes were thoroughly explained.

RESULTS

This review is the result of reviewing more than a hundred papers. It summarizes various problems that use molecular dynamics simulations.

CONCLUSION

The findings of this review highlight how molecular dynamics simulations have been successfully implemented to study the structure-function details of specific drug-target complexes. It also identifies the key areas such as stability of drug-target complexes, ligand binding kinetics and identification of allosteric sites which have been elucidated using molecular dynamics simulations.

Cisplatin reacts with histone H1 and the adduct forms a ternary complex with DNA

Cisplatin is an anticancer drug widely used in clinics; it induces the apoptosis of cancer cells by targeting DNA. However, its interaction with proteins has been found to be crucial in modulating the pre and post-target activity. Nuclear DNA is tightly assembled with histone proteins to form nucleosomes in chromatin; this can impede the drug to access DNA. On the other hand, the linker histone H1 is considered 'the gate to nucleosomal DNA' due to its exposed location and dynamic conformation; therefore, this protein can influence the platination of DNA. In this study, we performed a reaction of cisplatin with histone H1 and investigated the interaction of the H1/cisplatin adduct with DNA. The reactions were conducted on the N-terminal domains of H1.4 (sequence 1-90, H1N90) and H1.0 (sequence 1-7, H1N7). The results show that H1 readily reacts with cisplatin and generates bidentate and tridentate adducts, with methionine and glutamate residues as the preferential binding sites. Chromatographic and NMR analyses show that the platination rate of H1 is slightly higher than that of DNA and the platinated H1 can form H1-cisplatin-DNA ternary complexes. Interestingly, cisplatin is more prone to form H1-Pt-DNA ternary complexes than trans-oriented platinum agents. The formation of H1-cisplatin-DNA ternary complexes and their preference for cis- over trans-oriented platinum agents suggest an important role of histone H1 in the mechanism of action of cisplatin.

Oxidation State Dependent Conformational Changes of HMGB1 Regulate the Formation of the CXCL12/HMGB1 Heterocomplex

High-mobility Group Box 1 (HMGB1) is an abundant protein present in all mammalian cells and involved in several processes. During inflammation or tissue damage, HMGB1 is released in the extracellular space and, depending on its redox state, can form a heterocomplex with CXCL12. The heterocomplex acts exclusively via the chemokine receptor CXCR4 enhancing leukocyte recruitment.

Here, we used multi-microsecond molecular dynamics (MD) simulations to elucidate the effect of the disulfide bond on the structure and dynamics of HMGB1.

The results of the MD simulations show that the presence or lack of the disulfide bond between Cys23 and Cys45 modulates the conformational space explored by HMGB1, making the reduced protein more suitable to form a complex with CXCL12.

Clinical significance of high-mobility group box-1 (HMGB1) in subjects with type 2 diabetes mellitus (T2DM) combined with chronic obstructive pulmonary disease (COPD)

Background

Simple method to predict type 2 diabetes mellitus (T2DM) combined with chronic obstructive pulmonary disease (COPD) is in great need clinically. This study aims to assess the clinical significance of high‐mobility group box‐1 (HMGB1) in predicting T2DM combined with COPD in Chinese patients with T2DM or COPD.

Methods

Serum concentrations of glycated hemoglobin (HbA1C), fasting plasma glucose (FPG), fasting insulin (FINS), total cholesterol (TC), triglyceride (TG), low‐density lipoprotein cholesterol (LDL‐c), high‐density lipoprotein cholesterol (HDL‐c), C‐reactive protein (CRP), fibrinogen (FIB), HMGB1, white blood cell count (WBC), neutrophil% (NEU%), and lung function text such as forced expiratory volume 1/forced vital capacity (FEV1/FVC) and forced expiratory volume 1% predicted value (FEV1%pred) were measured in 126 T2DM patients, 118 COPD patients, 112 T2DM combined with COPD patients, and 120 healthy controls. Logistic regression was used to estimate the risk factors for T2DM combined with COPD.

Results

High‐mobility group box‐1 elevated in patients with T2DM combined with COPD, significantly higher than other subjects (P < 0.05), and differences in HMGB1 also existed between patients with T2DM or COPD and healthy individuals (P < 0.01). HMGB1 was positively correlated with HOMA‐IR, FBG, and HbA1c (P < 0.01) and negatively correlated with FEV1/FVC and FEV1%pred (P < 0.01). Logistic regression showed that HMGB1 was identified to be independent risk factor for T2DM combined with COPD.

Conclusion

High‐mobility group box‐1 was independent risk factor for T2DM combined with COPD and can be served to predict the occurrence of T2DM combined with COPD.