CHARMM-GUI PDB Reader and Manipulator: Covalent Ligand Modeling and Simulation

 Lingyang Kong, Sang-Jun Park Wonpil Im

Journal of Molecular Biology, 2024

https://doi.org/10.1016/j.jmb.2024.168554

Molecular modeling and simulation serve an important role in exploring biological functions of proteins at the molecular level, which is complementary to experiments. CHARMM-GUI (https://www.charmm-gui.org) is a web-based graphical user interface that generates molecular simulation systems and input files, and we have been continuously developing and extending its functionalities to facilitate various complex modeling and make molecular dynamics simulations more accessible to the scientific community. Currently, covalent drug discovery emerges as a popular and important field. Covalent drug forms a chemical bond with specific residues on the target protein, and it has advantages in potency for its prolonged inhibition effects. Even though there are higher demands in modeling PDB protein structures with various covalent ligand types, proper modeling of covalent ligands remains challenging. This work presents a new functionality in CHARMM-GUI PDB Reader & Manipulator that can handle a diversity of ligand-amino acid linkage types, which is validated by a careful benchmark study using over 1,000 covalent ligand structures in RCSB PDB. We hope that this new functionality can boost the modeling and simulation study of covalent ligands.

Restricted Rotational Flexibility of the C5α-Methyl-Substituted Carbapenem NA-1-157 Leads to Potent Inhibition of the GES-5 Carbapenemase

Nichole K. Stewart, Marta Toth, Pojun Quan, Michael Beer, John D. Buynak, Clyde A. Smith, and Sergei B. Vakulenko

ACS Infectious Diseases 2024

DOI: 10.1021/acsinfecdis.3c00683

Carbapenem antibiotics are used as a last-resort treatment for infections caused by multidrug-resistant bacteria. The wide spread of carbapenemases in Gram-negative bacteria has severely compromised the utility of these drugs and represents a serious public health threat. To combat carbapenemase-mediated resistance, new antimicrobials and inhibitors of these enzymes are urgently needed. Here, we describe the interaction of the atypically C5α-methyl-substituted carbapenem, NA-1-157, with the GES-5 carbapenemase. MICs of this compound against Escherichia coli, Klebsiella pneumoniae, and Acinetobacter baumannii producing the enzyme were reduced 4–16-fold when compared to MICs of the commercial carbapenems, reaching clinically sensitive breakpoints. When NA-1-157 was combined with meropenem, a strong synergistic effect was observed. Kinetic and ESI-LC/MS studies demonstrated that NA-1-157 is a potent inhibitor of GES-5, with a high inactivation efficiency of (2.9 ± 0.9) × 105 M–1 s–1. Acylation of GES-5 by NA-1-157 was biphasic, with the fast phase completing within seconds, and the slow phase taking several hours and likely proceeding through a reversible tetrahedral intermediate. Deacylation was extremely slow (k3 = (2.4 ± 0.3) × 10–7 s–1), resulting in a residence time of 48 ± 6 days. MD simulation of the GES-5-meropenem and GES-5-NA-1-157 acyl-enzyme complexes revealed that the C5α-methyl group in NA-1-157 sterically restricts rotation of the 6α-hydroxyethyl group preventing ingress of the deacylating water into the vicinity of the scissile bond of the acyl-enzyme intermediate. These data demonstrate that NA-1-157 is a potent irreversible inhibitor of the GES-5 carbapenemase.

Discovery of a Covalent Inhibitor Selectively Targeting the Autophosphorylation Site of c-Src Kinase

Huimin Zhang, Dounan Xu, Hongchan Huang, Hao Jiang, Linghao Hu, Liping Liu, Ge Sun, Jing Gao, Yuanqing Li, Cuicui Xia, Shijie Chen, Hu Zhou, Xiangqian Kong, Mingliang Wang, and Cheng Luo

ACS Chemical Biology 2024

DOI: 10.1021/acschembio.4c00048

Nonreceptor tyrosine kinase c-Src plays a crucial role in cell signaling and contributes to tumor progression. However, the development of selective c-Src inhibitors turns out to be challenging. In our previous study, we performed posttranslational modification-inspired drug design (PTMI-DD) to provide a plausible way for designing selective kinase inhibitors. In this study, after identifying a unique pocket comprising a less conserved cysteine and an autophosphorylation site in c-Src as well as a promiscuous covalent inhibitor, chemical optimization was performed to obtain (R)-LW-Srci-8 with nearly 75-fold improved potency (IC50 = 35.83 ± 7.21 nM). Crystallographic studies revealed the critical C–F···C═O interactions that may contribute to tight binding. The kinact and Ki values validated the improved binding affinity and decreased warhead reactivity of (R)-LW-Srci-8 for c-Src. Notably, in vitro tyrosine kinase profiling and cellular activity-based protein profiling (ABPP) cooperatively indicated a specific inhibition of c-Src by (R)-LW-Srci-8. Intriguingly, (R)-LW-Srci-8 preferentially binds to inactive c-Src with unphosphorylated Y419 both in vitro and in cells, subsequently disrupting the autophosphorylation. Collectively, our study demonstrated the feasibility of developing selective kinase inhibitors by cotargeting a nucleophilic residue and a posttranslational modification site and providing a chemical probe for c-Src functional studies.

Redirecting the pioneering function of FOXA1 with covalent small molecules

Sang Joon Won, Yuxiang Zhang, Christopher J. Reinhardt, Nicole S. MacRae, Kristen E. DeMeester, Evert Njomen, Lauren M. Hargis, Jarrett R. Remsberg, Bruno Melillo, Benjamin F. Cravatt, Michael A. Erb

bioRxiv 2024.03.21.586158; doi: https://doi.org/10.1101/2024.03.21.586158

Pioneer transcription factors (TFs) exhibit a specialized ability to bind to and open closed chromatin, facilitating engagement by other regulatory factors involved in gene activation or repression. Chemical probes are lacking for pioneer TFs, which has hindered their mechanistic investigation in cells. Here, we report the chemical proteomic discovery of electrophilic small molecules that stereoselectively and site-specifically bind the pioneer TF, FOXA1, at a cysteine (C258) within the forkhead DNA-binding domain. We show that these covalent ligands react with FOXA1 in a DNA-dependent manner and rapidly remodel its pioneer activity in prostate cancer cells reflected in redistribution of FOXA1 binding across the genome and directionally correlated changes in chromatin accessibility. Motif analysis supports a mechanism where the covalent ligands relax the canonical DNA binding preference of FOXA1 by strengthening interactions with suboptimal ancillary sequences in predicted proximity to C258. Our findings reveal a striking plasticity underpinning the pioneering function of FOXA1 that can be controlled by small molecules.

Advancing protein therapeutics through proximity-induced chemistry

Linqi Cheng Yixian Wang, Yiming Guo, Sophie S. Zhang Han Xiao
Cell Chemical Biology, 2024 Volume 31,  3, 428 - 445

https://doi.org/10.1016/j.chembiol.2023.09.004

Recent years have seen a remarkable growth in the field of protein-based medical treatments. Nevertheless, concerns have arisen regarding the cytotoxicity limitations, low affinity, potential immunogenicity, low stability, and challenges to modify these proteins. To overcome these obstacles, proximity-induced chemistry has emerged as a next-generation strategy for advancing protein therapeutics. This method allows site-specific modification of proteins with therapeutic agents, improving their effectiveness without extensive engineering. In addition, this innovative approach enables spatial control of the reaction based on proximity, facilitating the formation of irreversible covalent bonds between therapeutic proteins and their targets. This capability becomes particularly valuable in addressing challenges such as the low affinity frequently encountered between therapeutic proteins and their targets, as well as the limited availability of small molecules for specific protein targets. As a result, proximity-induced chemistry is reshaping the field of protein drug preparation and propelling the revolution in novel protein therapeutics.

Molecular Insights into the Impact of Mutations on the Binding Affinity of Targeted Covalent Inhibitors of BTK

 Ernest Awoonor-Williams and Abd Al-Aziz A. Abu-Saleh

The Journal of Physical Chemistry B 2024

DOI: 10.1021/acs.jpcb.4c00310

Targeted covalent inhibitors (TCIs) have witnessed a significant resurgence in recent years, particularly in the kinase drug discovery field for treating diverse clinical indications. The inhibition of Bruton’s tyrosine kinase (BTK) for treating B-cell cancers is a classic example where TCIs such as ibrutinib have had breakthroughs in targeted therapy. However, selectivity remains challenging, and the emergence of resistance mutations is a critical concern for clinical efficacy. Computational methods that can accurately predict the impact of mutations on inhibitor binding affinity could prove helpful in informing targeted approaches─providing insights into drug resistance mechanisms. In addition, such systems could help guide the systematic evaluation and impact of mutations in disease models for optimal experimental design. Here, we have employed in silico physics-based methods to understand the effects of mutations on the binding affinity and conformational dynamics of select TCIs of BTK. The TCIs studied include ibrutinib, acalabrutinib, and zanubrutinib─all of which are FDA-approved drugs for treating multiple forms of leukemia and lymphoma. Our results offer useful molecular insights into the structural determinants, thermodynamics, and conformational energies that impact ligand binding for this biological target of clinical relevance.

Covalent Fragment Screening and Optimization Identifies the Chloroacetohydrazide Scaffold as Inhibitors for Ubiquitin C-terminal Hydrolase L1

Ryan D. Imhoff, Rishi Patel, Muhammad Hassan Safdar, Hannah B. L. Jones, Adan Pinto-Fernandez, Iolanda Vendrell, Hao Chen, Christine S. Muli, Aaron D. Krabill, Benedikt M. Kessler, Michael K. Wendt, Chittaranjan Das, and Daniel P. Flaherty

Journal of Medicinal Chemistry 2024

DOI: 10.1021/acs.jmedchem.3c01661

Dysregulation of the ubiquitin-proteasome systems is a hallmark of various disease states including neurodegenerative diseases and cancer. Ubiquitin C-terminal hydrolase L1 (UCHL1), a deubiquitinating enzyme, is expressed primarily in the central nervous system under normal physiological conditions, however, is considered an oncogene in various cancers, including melanoma, lung, breast, and lymphoma. Thus, UCHL1 inhibitors could serve as a viable treatment strategy against these aggressive cancers. Herein, we describe a covalent fragment screen that identified the chloroacetohydrazide scaffold as a covalent UCHL1 inhibitor. Subsequent optimization provided an improved fragment with single-digit micromolar potency against UCHL1 and selectivity over the closely related UCHL3. The molecule demonstrated efficacy in cellular assays of metastasis. Additionally, we report a ligand-bound crystal structure of the most potent molecule in complex with UCHL1, providing insight into the binding mode and information for future optimization.

An orally bioavailable SARS-CoV-2 main protease inhibitor exhibits improved affinity and reduced sensitivity to mutations

Michael Westberg et al.

Sci. Transl. Med.16,eadi0979(2024).

DOI:10.1126/scitranslmed.adi0979

Inhibitors of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main protease (Mpro) such as nirmatrelvir (NTV) and ensitrelvir (ETV) have proven effective in reducing the severity of COVID-19, but the presence of resistance-conferring mutations in sequenced viral genomes raises concerns about future drug resistance. Second-generation oral drugs that retain function against these mutants are thus urgently needed. We hypothesized that the covalent hepatitis C virus protease inhibitor boceprevir (BPV) could serve as the basis for orally bioavailable drugs that inhibit SARS-CoV-2 Mpro more efficiently than existing drugs. Performing structure-guided modifications of BPV, we developed a picomolar-affinity inhibitor, ML2006a4, with antiviral activity, oral pharmacokinetics, and therapeutic efficacy similar or superior to those of NTV. A crucial feature of ML2006a4 is a derivatization of the ketoamide reactive group that improves cell permeability and oral bioavailability. Last, ML2006a4 was found to be less sensitive to several mutations that cause resistance to NTV or ETV and occur in the natural SARS-CoV-2 population. Thus, anticipatory design can preemptively address potential resistance mechanisms to expand future treatment options against coronavirus variants.

Strain-release alkylation of Asp12 enables mutant selective targeting of K-Ras-G12D

Qinheng Zheng, Ziyang Zhang, Keelan Z. Guiley & Kevan M. Shokat

Nat Chem Biol 2024

 https://doi.org/10.1038/s41589-024-01565-w

K-Ras is the most commonly mutated oncogene in human cancer. The recently approved non-small cell lung cancer drugs sotorasib and adagrasib covalently capture an acquired cysteine in K-Ras-G12C mutation and lock it in a signaling-incompetent state. However, covalent inhibition of G12D, the most frequent K-Ras mutation particularly prevalent in pancreatic ductal adenocarcinoma, has remained elusive due to the lack of aspartate-targeting chemistry. Here we present a set of malolactone-based electrophiles that exploit ring strain to crosslink K-Ras-G12D at the mutant aspartate to form stable covalent complexes. Structural insights from X-ray crystallography and exploitation of the stereoelectronic requirements for attack of the electrophile allowed development of a substituted malolactone that resisted attack by aqueous buffer but rapidly crosslinked with the aspartate-12 of K-Ras in both GDP and GTP state. The GTP-state targeting allowed effective suppression of downstream signaling, and selective inhibition of K-Ras-G12D-driven cancer cell proliferation in vitro and xenograft growth in mice.

Chemical tools to expand the ligandable proteome: diversity-oriented synthesis-based photoreactive stereoprobes

Daisuke Ogasawara, David Konrad, Zher Yin Tan, Kimberly Carey, Jessica Luo, Sang Joon Won, Haoxin Li, Trever Carter, Kristen DeMeester, Evert Njomen, Stuart Schreiber, Ramnik Xavier, Bruno Melillo, Benjamin Cravatt

bioRxiv 2024.02.27.582206; 

doi: https://doi.org/10.1101/2024.02.27.582206

Chemical proteomics enables the global assessment of small molecule-protein interactions in native biological systems and has emerged as a versatile approach for ligand discovery. The range of small molecules explored by chemical proteomics has, however, been limited. Here, we describe a diversity-oriented synthesis (DOS)-inspired library of stereochemically-defined compounds bearing diazirine and alkyne units for UV light-induced covalent modification and click chemistry enrichment of interacting proteins, respectively. We find that these photo-stereoprobes interact in a stereoselective manner with hundreds of proteins from various structural and functional classes in human cells and demonstrate that these interactions can form the basis for high-throughput screening-compatible nanoBRET assays. Integrated phenotypic analysis and chemical proteomics identified photo-stereoprobes that modulate autophagy by engaging the mitochondrial serine protease CLPP. Our findings show the utility of photo-stereoprobes for expanding the ligandable proteome, furnishing target engagement assays, and discovering and characterizing bioactive small molecules by cell-based screening.

Formaldehyde regulates S-adenosylmethionine biosynthesis and one-carbon metabolism

VANHA N. PHAM, KEVIN J. BRUEMMER, JOEL D. W. TOH, EVA J. GE, LOGAN TENNEY, CARL C. WARD, FELIX A. DINGLER, CHRISTOPHER L. MILLINGTON, CARLOS A. GARCIA-PRIETO  MIA C. PULOS-HOLMES  NICHOLAS T. INGOLIA LUCAS B. PONTEL MANEL ESTELLER  KETAN J. PATEL  DANIEL K. NOMURA AND CHRISTOPHER J. CHANG 

Science 2023 382, 6670

DOI: 10.1126/science.abp9201

INTRODUCTION

One-carbon metabolism manages cellular carbon pools by detoxifying highly reactive carbon species, such as aldehydes, and diverting their carbon toward the biosynthesis of useful products, including amino acids and nucleotides. Formaldehyde (FA) is a major one-carbon unit derived from exogenous environmental exposure and endogenous sources and is quickly scavenged in the cell through enzymatic oxidation to formate and carbon dioxide and/or metabolized through the folate cycle. S-adenosylmethionine (SAM) serves as the primary cellular methyl donor and harbors one-carbon units in a stable and accessible form. The ability to decipher the biochemical interplay between toxic reactive carbon species and stable physiological carbon units is essential for understanding fundamentals of one-carbon metabolism across all kingdoms of life. Especially important is understanding how aberrant carbon imbalances are connected to human diseases such as cancer, liver diseases, and asthma. Although the chronic exposure of toxic aldehydes is correlated to disease states, biological mechanisms of aldehyde signaling and their relation to carbon metabolism remain underexplored.

RATIONALE

Owing to its highly electrophilic nature, we hypothesized that FA could act as a one-carbon signal sensed by privileged cysteine sites across the proteome. FA reacts with cysteines on synthetic peptides, and we designed an unbiased, proteome-wide profiling study to systematically identify FA-sensitive cysteine residues. This work builds a biochemical framework for understanding global FA reactivity as a selective posttranslational modification of target proteins and downstream regulatory effects of such modifications.

RESULTS

Activity-based protein profiling identified FA modification of privileged cysteine sites across the proteome, including several enzymes responsible for FA metabolism, one-carbon metabolism, and amino acid biosynthesis. We focused on biochemical characterization of a key Cys120 residue on the SAM-generating enzyme S-adenosylmethionine synthase isoform type-1 (MAT1A) that is proximal to the MAT1A active site. FA exposure resulted in inhibition of MAT1A activity in an isoform-specific manner, which led to decreased SAM production. Cellular models containing only the MAT1A isoform displayed a reciprocal decrease in SAM levels with increasing doses of FA exposure. Moreover, an Adh5–/– mouse model of chronic FA elevation also showed SAM deficiency accompanied by lower levels of methylation on select histone methyl sinks. The chronic FA model also resulted in a decrease in methylation of the Mat1a promoter region, resulting in increased MAT1A expression as a compensatory mechanism to maintain available carbon units. We deciphered a compensatory biochemical feedback cycle where FA-dependent SAM deficiency led to an increase in MAT1A expression through genetic and epigenetic mechanisms regulated by FA-dependent transcription factors and DNA promoter hypomethylation, respectively.

CONCLUSION

In contrast to the traditional view of FA as an indiscriminate electrophile and toxic metabolite, we show that FA is sensed by specific cysteine sites in the proteome to regulate one-carbon metabolism feedback cycles through SAM biosynthesis. FA reacts with a key cysteine residue on MAT1A to inhibit its activity, resulting in SAM depletion and downstream changes in histone and DNA methylation. Under normal homeostatic conditions, FA is quickly sequestered into the folate cycle for conservation of one-carbon units to maintain balanced SAM biosynthesis. In response to FA overload, reciprocal SAM depletion through isoform-specific MAT1A inhibition results in changes to cellular methylation potential, epigenetic dysregulation, and perturbations in one-carbon metabolism, which in turns leads to compensatory up-regulation of MAT1A expression. This work provides a starting point for further exploration of aldehydes as signaling agents and the nexus between one-carbon metabolism and one-carbon signaling.