Active Site Breathing of Human Alkbh5 Revealed by Solution NMR and Accelerated Molecular Dynamics

The role of m6A demethylases in lung cancer: diagnostic and therapeutic implications

m6A is the most prevalent internal modification of eukaryotic mRNA, and plays a crucial role in tumorigenesis and various other biological processes. Lung cancer is a common primary malignant tumor of the lungs, which involves multiple factors in its occurrence and progression. Currently, only the demethylases FTO and ALKBH5 have been identified as associated with m6A modification. These demethylases play a crucial role in regulating the growth and invasion of lung cancer cells by removing methyl groups, thereby influencing stability and translation efficiency of mRNA. Furthermore, they participate in essential biological signaling pathways, making them potential targets for intervention in lung cancer treatment. Here we provides an overview of the involvement of m6A demethylase in lung cancer, as well as their potential application in the diagnosis, prognosis and treatment of the disease.

Solution Structure Ensembles of the Open and Closed Forms of the ∼130 kDa Enzyme I via AlphaFold Modeling, Coarse Grained Simulations, and NMR

Large-scale interdomain rearrangements are essential to protein function, governing the activity of large enzymes and molecular machineries. Yet, obtaining an atomic-resolution understanding of how the relative domain positioning is affected by external stimuli is a hard task in modern structural biology. Here, we show that combining structural modeling by AlphaFold2 with coarse-grained molecular dynamics simulations and NMR residual dipolar coupling data is sufficient to characterize the spatial domain organization of bacterial enzyme I (EI), a ∼130 kDa multidomain oligomeric protein that undergoes large-scale conformational changes during its catalytic cycle. In particular, we solve conformational ensembles for EI at two different experimental temperatures and demonstrate that a lower temperature favors sampling of the catalytically competent closed state of the enzyme. These results suggest a role for conformational entropy in the activation of EI and demonstrate the ability of our protocol to detect and characterize the effect of external stimuli (such as mutations, ligand binding, and post-translational modifications) on the interdomain organization of multidomain proteins. We expect the ensemble refinement protocol described here to be easily transferrable to the investigation of the structure and dynamics of other uncharted multidomain systems and have assembled a Google Colab page (https://potoyangroup.github.io/Seq2Ensemble/) to facilitate implementation of the presented methodology elsewhere.

The regulation of N6-methyladenosine modification in PD-L1-induced anti-tumor immunity

There is growing evidence that programmed death ligand-1 (PD-L1) has exciting therapeutic efficacy in hematological malignancy and partial solid tumors. However, many patients still face failure with the treatment of immune checkpoint blockade because of PD-L1 expression regulation during transcription and post-transcription processes, including N6-methyladenosine (m6A). Similar to the epigenetic regulation in DNA and histones, recent research has revealed the essential regulation of m6A modification in RNA nuclear export, metabolism and translation. Recent studies have shown that m6A-induced PD-L1 expression emerges as one of the main reasons for the immunological alteration in this process and contributes to the failure of T cell-induced anti-tumor immunity. The results of preclinical studies demonstrate the potential of m6A-targeted therapy in combination with immune checkpoint blockade. The comprehensive expression of m6A-related genes also provided the possibility to indicate the prognosis and to optimize the treatment for patients of various cancer types. In this review, we focus on the m6A modification in PD-L1 mRNA as well as the regulation of PD-L1 expression in cancer cells and summarize its clinical value in anti-PD-L1 cancer immune therapy.

Spectroscopic studies reveal details of substrate-induced conformational changes distant from the active site in isopenicillin N synthase

Isopenicillin N synthase (IPNS) catalyzes formation of the β-lactam and thiazolidine rings of isopenicillin N from its linear tripeptide l-δ-(α-aminoadipoyl)-l-cysteinyl-d-valine (ACV) substrate in an iron- and dioxygen (O₂)-dependent four-electron oxidation without precedent in current synthetic chemistry. Recent X-ray free-electron laser studies including time-resolved serial femtosecond crystallography show that binding of O₂ to the IPNS–Fe(II)–ACV complex induces unexpected conformational changes in α-helices on the surface of IPNS, in particular in α3 and α10. However, how substrate binding leads to conformational changes away from the active site is unknown. Here, using detailed ¹⁹F NMR and electron paramagnetic resonance experiments with labeled IPNS variants, we investigated motions in α3 and α10 induced by binding of ferrous iron, ACV, and the O₂ analog nitric oxide, using the less mobile α6 for comparison. ¹⁹F NMR studies were carried out on singly and doubly labeled α3, α6, and α10 variants at different temperatures. In addition, double electron–electron resonance electron paramagnetic resonance analysis was carried out on doubly spin-labeled variants. The combined spectroscopic and crystallographic results reveal that substantial conformational changes in regions of IPNS including α3 and α10 are induced by binding of ACV and nitric oxide. Since IPNS is a member of the structural superfamily of 2-oxoglutarate-dependent oxygenases and related enzymes, related conformational changes may be of general importance in nonheme oxygenase catalysis.

Pan-Cancer Analysis Shows That ALKBH5 Is a Potential Prognostic and Immunotherapeutic Biomarker for Multiple Cancer Types Including Gliomas

Background

AlkB homolog 5 (ALKBH5) is a N⁶-methyladenosine (m⁶A) demethylase associated with the development, growth, and progression of multiple cancer types. However, the biological role of ALKBH5 has not been investigated in pan-cancer datasets. Therefore, in this study, comprehensive bioinformatics analysis of pan-cancer datasets was performed to determine the mechanisms through which ALKBH5 regulates tumorigenesis.

Methods

Online websites and databases such as NCBI, UCSC, CCLE, HPA, TIMER2, GEPIA2, cBioPortal, UALCAN, STRING, SangerBox, ImmuCellAl, xCell, and GenePattern were used to extract data of ALKBH5 in multiple cancers. The pan-cancer patient datasets were analyzed to determine the relationship between ALKBH5 expression, genetic alterations, methylation status, and tumor immunity. Targetscan, miRWalk, miRDB, miRabel, LncBase databases and Cytoscape tool were used to identify microRNAs (miRNAs) and long non-coding RNAs (lncRNAs) that regulate expression of ALKBH5 and construct the lncRNA-miRNA-ALKBH5 network. In vitro CCK-8, wound healing, Transwell and M2 macrophage infiltration assays as well as in vivo xenograft animal experiments were performed to determine the biological functions of ALKBH5 in glioma cells.

Results

The pan-cancer analysis showed that ALKBH5 was upregulated in several solid tumors. ALKBH5 expression significantly correlated with the prognosis of cancer patients. Genetic alterations including duplications and deep mutations of the ALKBH5 gene were identified in several cancer types. Alterations in the ALKBH5 gene correlated with tumor prognosis. GO and KEGG enrichment analyses showed that ALKBH5-related genes were enriched in the inflammatory, metabolic, and immune signaling pathways in glioma. ALKBH5 expression correlated with the expression of immune checkpoint (ICP) genes, and influenced sensitivity to immunotherapy. We constructed a lncRNA-miRNA network that regulates ALKBH5 expression in tumor development and progression. In vitro and in vivo experiments showed that ALKBH5 promoted proliferation, migration, and invasion of glioma cells and recruited the M2 macrophage to glioma cells.

Conclusions

ALKBH5 was overexpressed in multiple cancer types and promoted the development and progression of cancers through several mechanisms including regulation of the tumor-infiltration of immune cells. Our study shows that ALKBH5 is a promising prognostic and immunotherapeutic biomarker in some malignant tumors.

Gaussian accelerated molecular dynamics (GaMD): principles and applications

Gaussian accelerated molecular dynamics (GaMD) is a robust computational method for simultaneous unconstrained enhanced sampling and free energy calculations of biomolecules. It works by adding a harmonic boost potential to smooth biomolecular potential energy surface and reduce energy barriers. GaMD greatly accelerates biomolecular simulations by orders of magnitude. Without the need to set predefined reaction coordinates or collective variables, GaMD provides unconstrained enhanced sampling and is advantageous for simulating complex biological processes. The GaMD boost potential exhibits a Gaussian distribution, thereby allowing for energetic reweighting via cumulant expansion to the second order (i.e., "Gaussian approximation"). This leads to accurate reconstruction of free energy landscapes of biomolecules. Hybrid schemes with other enhanced sampling methods, such as the replica exchange GaMD (rex-GaMD) and replica exchange umbrella sampling GaMD (GaREUS), have also been introduced, further improving sampling and free energy calculations. Recently, new "selective GaMD" algorithms including the ligand GaMD (LiGaMD) and peptide GaMD (Pep-GaMD) enabled microsecond simulations to capture repetitive dissociation and binding of small-molecule ligands and highly flexible peptides. The simulations then allowed highly efficient quantitative characterization of the ligand/peptide binding thermodynamics and kinetics. Taken together, GaMD and its innovative variants are applicable to simulate a wide variety of biomolecular dynamics, including protein folding, conformational changes and allostery, ligand binding, peptide binding, protein-protein/nucleic acid/carbohydrate interactions, and carbohydrate/nucleic acid interactions. In this review, we present principles of the GaMD algorithms and recent applications in biomolecular simulations and drug design.

N 6-methyladenosine binding induces a metal-centered rearrangement that activates the human RNA demethylase Alkbh5

Alkbh5 catalyzes demethylation of the N ⁶-methyladenosine (m⁶A), an epigenetic mark that controls several physiological processes including carcinogenesis and stem cell differentiation. The activity of Alkbh5 comprises two coupled reactions. The first reaction involves decarboxylation of α-ketoglutarate (αKG) and formation of a Fe⁴⁺═O species. This oxyferryl intermediate oxidizes the m⁶A to reestablish the canonical base. Despite coupling between the two reactions being required for the correct Alkbh5 functioning, the mechanisms linking dioxygen activation to m⁶A binding are not fully understood. Here, we use solution NMR to investigate the structure and dynamics of apo and holo Alkbh5. We show that binding of m⁶A to Alkbh5 induces a metal-centered rearrangement of αKG that increases the exposed area of the metal, making it available for binding O₂ Our study reveals the molecular mechanisms underlying activation of Alkbh5, therefore opening new perspectives for the design of novel strategies to control gene expression and cancer progression.

A Single Point Mutation Controls the Rate of Interconversion Between the g + and g - Rotamers of the Histidine 189 χ2 Angle That Activates Bacterial Enzyme I for Catalysis

Enzyme I (EI) of the bacterial phosphotransferase system (PTS) is a master regulator of bacterial metabolism and a promising target for development of a new class of broad-spectrum antibiotics. The catalytic activity of EI is mediated by several intradomain, interdomain, and intersubunit conformational equilibria. Therefore, in addition to its relevance as a drug target, EI is also a good model for investigating the dynamics/function relationship in multidomain, oligomeric proteins. Here, we use solution NMR and protein design to investigate how the conformational dynamics occurring within the N-terminal domain (EIN) affect the activity of EI. We show that the rotameric g⁺-to-g⁻ transition of the active site residue His¹⁸⁹ χ2 angle is decoupled from the state A-to-state B transition that describes a ∼90° rigid-body rearrangement of the EIN subdomains upon transition of the full-length enzyme to its catalytically competent closed form. In addition, we engineered EIN constructs with modulated conformational dynamics by hybridizing EIN from mesophilic and thermophilic species, and used these chimeras to assess the effect of increased or decreased active site flexibility on the enzymatic activity of EI. Our results indicate that the rate of the autophosphorylation reaction catalyzed by EI is independent from the kinetics of the g⁺-to-g⁻ rotameric transition that exposes the phosphorylation site on EIN to the incoming phosphoryl group. In addition, our work provides an example of how engineering of hybrid mesophilic/thermophilic chimeras can assist investigations of the dynamics/function relationship in proteins, therefore opening new possibilities in biophysics.

Structure elucidation of the elusive Enzyme I monomer reveals the molecular mechanisms linking oligomerization and enzymatic activity

Enzyme I (EI) is a phosphotransferase enzyme responsible for converting phosphoenolpyruvate (PEP) into pyruvate. This reaction initiates a five-step phosphorylation cascade in the bacterial phosphotransferase (PTS) transduction pathway. Under physiological conditions, EI exists in an equilibrium between a functional dimer and an inactive monomer. The monomer-dimer equilibrium is a crucial factor regulating EI activity and the phosphorylation state of the overall PTS. Experimental studies of EI's monomeric state have yet been hampered by the dimer's high thermodynamic stability, which prevents its characterization by standard structural techniques. In this study, we modified the dimerization domain of EI (EIC) by mutating three amino acids involved in the formation of intersubunit salt bridges. The engineered variant forms an active dimer in solution that can bind and hydrolyze PEP. Using hydrostatic pressure as an additional perturbation, we were then able to study the complete dissociation of the variant from 1 bar to 2.5 kbar in the absence and the presence of EI natural ligands. Backbone residual dipolar couplings collected under high-pressure conditions allowed us to determine the conformational ensemble of the isolated EIC monomeric state in solution. Our calculations reveal that three catalytic loops near the dimerization interface become unstructured upon monomerization, preventing the monomeric enzyme from binding its natural substrate. This study provides an atomic-level characterization of EI's monomeric state and highlights the role of the catalytic loops as allosteric connectors controlling both the activity and oligomerization of the enzyme.

15N CPMG Relaxation Dispersion for the Investigation of Protein Conformational Dynamics on the µs-ms Timescale

Protein conformational dynamics play fundamental roles in regulation of enzymatic catalysis, ligand binding, allostery, and signaling, which are important biological processes. Understanding how the balance between structure and dynamics governs biological function is a new frontier in modern structural biology and has ignited several technical and methodological developments. Among these, CPMG relaxation dispersion solution NMR methods provide unique, atomic-resolution information on the structure, kinetics, and thermodynamics of protein conformational equilibria occurring on the µs-ms timescale. Here, the study presents detailed protocols for acquisition and analysis of a ¹⁵N relaxation dispersion experiment. As an example, the pipeline for the analysis of the µs-ms dynamics in the C-terminal domain of bacteria Enzyme I is shown.

Multi-substrate selectivity based on key loops and non-homologous domains: new insight into ALKBH family

AlkB homologs (ALKBH) are a family of specific demethylases that depend on Fe2+ and α-ketoglutarate to catalyze demethylation on different substrates, including ssDNA, dsDNA, mRNA, tRNA, and proteins. Previous studies have made great progress in determining the sequence, structure, and molecular mechanism of the ALKBH family. Here, we first review the multi-substrate selectivity of the ALKBH demethylase family from the perspective of sequence and structural evolution. The construction of the phylogenetic tree and the comparison of key loops and non-homologous domains indicate that the paralogs with close evolutionary relationship have similar domain compositions. The structures show that the lack and variations of four key loops change the shape of clefts to cause the differences in substrate affinity, and non-homologous domains may be related to the compatibility of multiple substrates. We anticipate that the new insights into selectivity determinants of the ALKBH family are useful for understanding the demethylation mechanisms.

Structural determinants of nucleobase modification recognition in the AlkB family of dioxygenases

Iron-dependent dioxygenases of the AlkB protein family found in most organisms throughout the tree of life play a major role in oxidative dealkylation processes. Many of these enzymes have attracted the attention of researchers across different fields and have been subjected to thorough biochemical characterization because of their link to human health and disease. For example, several mammalian AlkB homologues are involved in the direct reversal of alkylation damage in DNA, while others have been shown to play a regulatory role in epigenetic or epitranscriptomic nucleic acid methylation or in post-translational modifications such as acetylation of actin filaments. These studies show that that divergence in amino acid sequence and structure leads to different characteristics and substrate specificities. In this review, we aim to summarize current insights in the structural features involved in the substrate selection of AlkB homologues, with focus on nucleic acid interactions.

The biological function of m6A demethylase ALKBH5 and its role in human disease

Human AlkB homolog H5 (ALKBH5) is a primary m6A demethylase, which is dysregulated and acts as a biological and pharmacological role in human cancers or non-cancers. ALKBH5 plays a dual role in various cancers through regulating kinds of biological processes, such as proliferation, migration, invasion, metastasis and tumor growth. In addition, it takes a great part in human non-cancer, including reproductive system diseases. The underlying regulatory mechanisms of ALKBH5 that relys on m6A-dependent modification are implicated with long non-coding RNA, cancer stem cell, autophagy and hypoxia. ALKBH5 is also an independent prognostic indicator in various cancers. In this review, we summarized the current evidence on ALKBH5 in diverse human cancers or non-cancers and its potential as a prognostic target.

Molecular Dynamics Simulations with NAMD2

X-ray diffraction crystallography is the primary technique to determine the three-dimensional structures of biomolecules. Although a robust method, X-ray crystallography is not able to access the dynamical behavior of macromolecules. To do so, we have to carry out molecular dynamics simulations taking as an initial system the three-dimensional structure obtained from experimental techniques or generated using homology modeling. In this chapter, we describe in detail a tutorial to carry out molecular dynamics simulations using the program NAMD2. We chose as a molecular system to simulate the structure of human cyclin-dependent kinase 2.

Hybrid Thermophilic/Mesophilic Enzymes Reveal a Role for Conformational Disorder in Regulation of Bacterial Enzyme I

Conformational disorder is emerging as an important feature of biopolymers, regulating a vast array of cellular functions, including signaling, phase separation, and enzyme catalysis. Here we combine NMR, crystallography, computer simulations, protein engineering, and functional assays to investigate the role played by conformational heterogeneity in determining the activity of the C-terminal domain of bacterial Enzyme I (EIC). In particular, we design chimeric proteins by hybridizing EIC from thermophilic and mesophilic organisms, and we characterize the resulting constructs for structure, dynamics, and biological function. We show that EIC exists as a mixture of active and inactive conformations and that functional regulation is achieved by tuning the thermodynamic balance between active and inactive states. Interestingly, we also present a hybrid thermophilic/mesophilic enzyme that is thermostable and more active than the wild-type thermophilic enzyme, suggesting that hybridizing thermophilic and mesophilic proteins is a valid strategy to engineer thermostable enzymes with significant low-temperature activity.

Role of Structural Dynamics in Selectivity and Mechanism of Non-heme Fe(II) and 2-Oxoglutarate-Dependent Oxygenases Involved in DNA Repair

AlkB and its human homologue AlkBH2 are Fe(II)- and 2-oxoglutarate (2OG)-dependent oxygenases that repair alkylated DNA bases occurring as a consequence of reactions with mutagenic agents. We used molecular dynamics (MD) and combined quantum mechanics/molecular mechanics (QM/MM) methods to investigate how structural dynamics influences the selectivity and mechanisms of the AlkB- and AlkBH2-catalyzed demethylation of 3-methylcytosine (m3C) in single (ssDNA) and double (dsDNA) stranded DNA. Dynamics studies reveal the importance of the flexibility in both the protein and DNA components in determining the preferences of AlkB for ssDNA and of AlkBH2 for dsDNA. Correlated motions, including of a hydrophobic β-hairpin, are involved in substrate binding in AlkBH2-dsDNA. The calculations reveal that 2OG rearrangement prior to binding of dioxygen to the active site Fe is preferred over a ferryl rearrangement to form a catalytically productive Fe(IV)=O intermediate. Hydrogen atom transfer proceeds via a σ-channel in AlkBH2-dsDNA and AlkB-dsDNA; in AlkB-ssDNA, there is a competition between σ- and π-channels, implying that the nature of the complexed DNA has potential to alter molecular orbital interactions during the substrate oxidation. Our results reveal the importance of the overall protein-DNA complex in determining selectivity and how the nature of the substrate impacts the mechanism.