ATP-Magnesium Coordination: Protein Structure-Based Force Field Evaluation and Corrections

Investigating how HIV-1 antiretrovirals differentially behave as substrates and inhibitors of P-glycoprotein via molecular dynamics simulations

HIV-1 can rapidly infect the brain upon initial infection, establishing latent reservoirs that induce neuronal damage and/or death, resulting in HIV-Associated Neurocognitive Disorder. Though anti-HIV-1 antiretrovirals (ARVs) suppress viral load, the blood-brain barrier limits drug access to the brain, largely because of highly expressed efflux proteins like P-glycoprotein (P-gp). While no FDA-approved P-gp inhibitor currently exists, HIV-1 protease inhibitors show promise as partial P-gp inhibitors, potentially enhancing drug delivery to the brain. Herein, we employed docking and molecular dynamics simulations to elucidate key differences in P-gp's interactions with several antiretrovirals, including protease inhibitors, with known inhibitory or substrate-like behaviors towards P-gp. Our results led us to hypothesize new mechanistic details of small-molecule efflux by and inhibition of P-gp, where the "Lower Pocket" in P-gp's transmembrane domain serves as the primary initial site for small-molecule binding. Subsequently, this pocket merges with the more traditionally studied drug binding site-the "Upper Pocket"-thus funneling small-molecule drugs, such as ARVs, towards the Upper Pocket for efflux. Furthermore, our results reinforce the understanding that both binding energetics and changes in protein dynamics are crucial in discerning small molecules as non-substrates, substrates, or inhibitors of P-gp. Our findings indicate that interactions between P-gp and inhibitory ARVs induce bridging of transmembrane domain helices, impeding P-gp conformational changes and contributing to the inhibitory behavior of these ARVs. Overall, insights gained in this study could serve to guide the design of future P-gp-targeting therapeutics for a wide range of pathological conditions and diseases, including HIV-1.

Magnesium enhances aurintricarboxylic acid's inhibitory action on the plasma membrane Ca2+-ATPase

Our research aimed to elucidate the mechanism by which aurintricarboxylic acid (ATA) inhibits plasma membrane Ca2+-ATPase (PMCA), a crucial enzyme responsible for calcium transport. Given the pivotal role of PMCA in cellular calcium homeostasis, understanding how it is inhibited by ATA holds significant implications for potentially regulating physiopathological cellular processes in which this pump is involved. Our experimental findings revealed that ATA employs multiple modes of action to inhibit PMCA activity, which are influenced by ATP but also by the presence of calcium and magnesium ions. Specifically, magnesium appears to enhance this inhibitory effect. Our experimental and in-silico results suggest that, unlike those reported in other proteins, ATA complexed with magnesium (ATA·Mg) is the molecule that inhibits PMCA. In summary, our study presents a novel perspective and establishes a solid foundation for future research efforts aimed at the development of new pharmacological molecules both for PMCA and other proteins.

How Does Mg2+(aq) Interact with ATP(aq)? Biomolecular Structure through the Lens of Liquid-Jet Photoemission Spectroscopy

Liquid-jet photoemission spectroscopy (LJ-PES) allows for a direct probing of electronic structure in aqueous solutions. We show the applicability of the approach to biomolecules in a complex environment, exploring site-specific information on the interaction of adenosine triphosphate in the aqueous phase (ATP(aq)) with magnesium (Mg2+(aq)), highlighting the synergy brought about by the simultaneous analysis of different regions in the photoelectron spectrum. In particular, we demonstrate intermolecular Coulombic decay (ICD) spectroscopy as a new and powerful addition to the arsenal of techniques for biomolecular structure investigation. We apply LJ-PES assisted by electronic-structure calculations to study ATP(aq) solutions with and without dissolved Mg2+. Valence photoelectron data reveal spectral changes in the phosphate and adenine features of ATP(aq) due to interactions with the divalent cation. Chemical shifts in Mg 2p, Mg 2s, P 2p, and P 2s core-level spectra as a function of the Mg2+/ATP concentration ratio are correlated to the formation of [Mg(ATP) 2]6-(aq), [MgATP]2-(aq), and [Mg2ATP](aq) complexes, demonstrating the element sensitivity of the technique to Mg2+-phosphate interactions. The most direct probe of the intermolecular interactions between ATP(aq) and Mg2+(aq) is delivered by the emerging ICD electrons following ionization of Mg 1s electrons. ICD spectra are shown to sensitively probe ligand exchange in the Mg2+-ATP(aq) coordination environment. In addition, we report and compare P 2s data from ATP(aq) and adenosine mono- and diphosphate (AMP(aq) and ADP(aq), respectively) solutions, probing the electronic structure of the phosphate chain and the local environment of individual phosphate units in ATP(aq). Our results provide a comprehensive view of the electronic structure of ATP(aq) and Mg2+-ATP(aq) complexes relevant to phosphorylation and dephosphorylation reactions that are central to bioenergetics in living organisms.

Systematic QM/MM Study for Predicting 31P NMR Chemical Shifts of Adenosine Nucleotides in Solution and Stages of ATP Hydrolysis in a Protein Environment

NMR (nuclear magnetic resonance) spectroscopy allows for important atomistic insights into the structure and dynamics of biological macromolecules; however, reliable assignments of experimental spectra are often difficult. Herein, quantum mechanical/molecular mechanical (QM/MM) calculations can provide crucial support. A major problem for the simulations is that experimental NMR signals are time-averaged over much longer time scales, and since computed chemical shifts are highly sensitive to local changes in the electronic and structural environment, sufficiently large averages over representative structural ensembles are essential. This entails high computational demands for reliable simulations. For NMR measurements in biological systems, a nucleus of major interest is 31P since it is both highly present (e.g., in nucleic acids) and easily observable. The focus of our present study is to develop a robust and computationally cost-efficient framework for simulating 31P NMR chemical shifts of nucleotides. We apply this scheme to study the different stages of the ATP hydrolysis reaction catalyzed by p97. Our methodology is based on MM molecular dynamics (MM-MD) sampling, followed by QM/MM structure optimizations and NMR calculations. Overall, our study is one of the most comprehensive QM-based 31P studies in a protein environment and the first to provide computed NMR chemical shifts for multiple nucleotide states in a protein environment. This study sheds light on a process that is challenging to probe experimentally and aims to bridge the gap between measured and calculated NMR spectroscopic properties.

Biochemical and structural characterization of an inositol pyrophosphate kinase from a giant virus

Kinases that synthesize inositol phosphates (IPs) and pyrophosphates (PP-IPs) control numerous biological processes in eukaryotic cells. Herein, we extend this cellular signaling repertoire to viruses. We have biochemically and structurally characterized a minimalist inositol phosphate kinase (i.e., TvIPK) encoded by Terrestrivirus, a nucleocytoplasmic large ("giant") DNA virus (NCLDV). We show that TvIPK can synthesize inositol pyrophosphates from a range of scyllo- and myo-IPs, both in vitro and when expressed in yeast cells. We present multiple crystal structures of enzyme/substrate/nucleotide complexes with individual resolutions from 1.95 to 2.6 Å. We find a heart-shaped ligand binding pocket comprising an array of positively charged and flexible side chains, underlying the observed substrate diversity. A crucial arginine residue in a conserved "G-loop" orients the γ-phosphate of ATP to allow substrate pyrophosphorylation. We highlight additional conserved catalytic and architectural features in TvIPK, and support their importance through site-directed mutagenesis. We propose that NCLDV inositol phosphate kinases may have assisted evolution of inositol pyrophosphate signaling, and we discuss the potential biogeochemical significance of TvIPK in soil niches.

Polarizable AMOEBA Model for Simulating Mg2+·Protein·Nucleotide Complexes

Molecular mechanics (MM) simulations have the potential to provide detailed insights into the mechanisms of enzymes that utilize nucleotides as cofactors. In most cases, the activities of these enzymes also require the binding of divalent cations to catalytic sites. However, modeling divalent cations in MM simulations has been challenging. The inclusion of explicit polarization was considered promising, but despite improvements over nonpolarizable force fields and despite the inclusion of "Nonbonded-fix (NB-fix)" corrections, errors in interaction energies of divalent cations with proteins remain large. Importantly, the application of these models fails to reproduce the experimental structural data on Mg2+·Protein·ATP complexes. Focusing on these complexes, here we provide a systematic assessment of the polarizable AMOEBA model and recommend critical changes that substantially improve its predictive performance. Our key results are as follows. We first show that our recent revision of the AMOEBA protein model (AMOEBABIO18-HFC), which contains high field corrections (HFCs) to induced dipoles, dramatically improves Mg2+-protein interaction energies, reducing the mean absolute error (MAE) from 17 to 10 kcal/mol. This further supports the general applicability of AMOEBABIO18-HFC. The inclusion of many-body NB-fix corrections further reduces MAE to 6 kcal/mol, which amounts to less than 2% error. The errors are estimated with respect to vdW-inclusive density functional theory that we benchmark against CCSD(T) calculations and experiments. We also present a new model of ATP with revised polarization parameters to better capture its high field response, as well as new vdW and dihedral parameters. The ATP model accurately predicts experimental Mg2+-ATP binding free energy in the aqueous phase and provides new insights into how Mg2+ associates with ATP. Finally, we show that molecular dynamics (MD) simulations of Mg2+·Kinase·ATP complexes carried out with these improvements lead to a better agreement in global and local catalytic site structures between MD and X-ray crystallography.

MESPEUS: a database of metal coordination groups in proteins

MESPEUS is a freely accessible database which uses carefully selected metal coordination groups found in metalloprotein structures taken from the Protein Data Bank. The database contains geometrical information of metal sites within proteins, including 40 metal types. In order to completely determine the metal coordination, the symmetry-related units of a given protein structure are taken into account and are generated using the appropriate space group symmetry operations. This permits a more complete description of the metal coordination geometry by including all coordinating atoms. The user-friendly web interface allows users to directly search for a metal site of interest using several useful options, including searching for metal elements, metal-donor distances, coordination number, donor residue group, and structural resolution. These searches can be carried out singly or in combination. The details of a metal site and the metal site(s) in the whole structure can be graphically displayed using the interactive web interface. MESPEUS is automatically updated monthly by synchronizing with the PDB database. An investigation for the Mg-ATP interaction is given to demonstrate how MESPEUS can be used to extract information about metal sites by selecting structure and coordination features. MESPEUS is available at http://mespeus.nchu.edu.tw/.

Eukaryotic ribosome display for antibody discovery: A review

Monoclonal antibody biologics have significantly transformed the therapeutic landscape within the biopharmaceutical industry, partly due to the utilisation of discovery technologies such as the hybridoma method and phage display. While these established platforms have streamlined the development process to date, their reliance on cell transformation for antibody identification faces limitations related to library diversification and the constraints of host cell physiology. Cell-free systems like ribosome display offer a complementary approach, enabling antibody selection in a completely in vitro setting while harnessing enriched cellular molecular machinery. This review aims to provide an overview of the fundamental principles underlying the ribosome display method and its potential for advancing antibody discovery and development.

Comparative Atomistic Insights on Apo and ATP-I1171N/S/T in Nonsmall-Cell Lung Cancer

Anaplastic lymphoma kinase (ALK) rearrangements occur in about 5% of nonsmall cell lung cancer (NSCLC) patients. Despite being first recognized as EML4-ALK, fusions with several additional genes have been identified, all of which cause constitutive activation of the ALK kinase and subsequently lead to tumor development. ALK inhibitors first-line crizotinib, second-line ceritinib, and alectinib are effective against NSCLC patients with these rearrangements. Patients progressing on crizotinib had various mutations in the ALK kinase domain. ALK fusion proteins are activated by oligomerization through the fusion partner, which leads to the autophosphorylation of the kinase's domain and consequent downstream activation. The proposed computational study focuses on understanding the activation mechanism of ALK and ATP binding of wild-type (WT) and I1171N/S/T mutations. We analyzed the conformational change of ALK I1171N/S/T mutations and ATP binding using molecular docking and molecular dynamics simulation approaches. According to principal component analysis and free energy landscape, it is clear that I1171N/S/T mutations in Apo and ATP showed different energy minima/unstable structures compared to WT-Apo. The results revealed that I1171N/S/T mutations and ATP binding significantly supported a change toward an active-state conformation, whereas WT-Apo remained inactive. We demonstrated that I1171N/S/T mutations are persistent in an active state and independent of ATP. The I1171S/T mutations showed greater intermolecular H-bonds with ATP than WT-ATP. The molecular mechanics Poisson-Boltzmann surface area analysis revealed that the I1171N/S/T mutation binding energy was similar to that of WT-ATP. This study shows that I1171N/S/T can form stable bonds with ATP and may contribute to a constitutively active kinase. Based on the Y1278-C1097 H-bond and E1167-K1150 salt bridge interaction, I1171N strongly promotes the constitutively active kinase independent of ATP. This structural mechanism study will aid in understanding the oncogenic activity of ALK and the basis for improving the ALK inhibitors.

Purinergic system in cancer stem cells

Accumulating evidence supports the idea that cancer stem cells (CSCs) are those with the capacity to initiate tumors, generate phenotypical diversity, sustain growth, confer drug resistance, and orchestrate the spread of tumor cells. It is still controversial whether CSCs originate from normal stem cells residing in the tissue or cancer cells from the tumor bulk that have dedifferentiated to acquire stem-like characteristics. Although CSCs have been pointed out as key drivers in cancer, knowledge regarding their physiology is still blurry; thus, research focusing on CSCs is essential to designing novel and more effective therapeutics. The purinergic system has emerged as an important autocrine-paracrine messenger system with a prominent role at multiple levels of the tumor microenvironment, where it regulates cellular aspects of the tumors themselves and the stromal and immune systems. Recent findings have shown that purinergic signaling also participates in regulating the CSC phenotype. Here, we discuss updated information regarding CSCs in the purinergic system and present evidence supporting the idea that elements of the purinergic system expressed by this subpopulation of the tumor represent attractive pharmacological targets for proposing innovative anti-cancer therapies.

Architecture of the MKK6-p38α complex defines the basis of MAPK specificity and activation

The mitogen-activated protein kinase (MAPK) p38α is a central component of signaling in inflammation and the immune response and is, therefore, an important drug target. Little is known about the molecular mechanism of its activation by double phosphorylation from MAPK kinases (MAP2Ks), because of the challenge of trapping a transient and dynamic heterokinase complex. We applied a multidisciplinary approach to generate a structural model of p38α in complex with its MAP2K, MKK6, and to understand the activation mechanism. Integrating cryo-electron microscopy with molecular dynamics simulations, hydrogen-deuterium exchange mass spectrometry, and experiments in cells, we demonstrate a dynamic, multistep phosphorylation mechanism, identify catalytically relevant interactions, and show that MAP2K-disordered amino termini determine pathway specificity. Our work captures a fundamental step of cell signaling: a kinase phosphorylating its downstream target kinase.

Magnesium Deficiency and Cardiometabolic Disease

Magnesium (Mg²⁺) has many physiological functions within the body. These include important roles in maintaining cardiovascular functioning, where it contributes to the regulation of cardiac excitation-contraction coupling, endothelial functioning and haemostasis. The haemostatic roles of Mg²⁺ impact upon both the protein and cellular arms of coagulation. In this review, we examine how Mg²⁺ homeostasis is maintained within the body and highlight the various molecular roles attributed to Mg²⁺ in the cardiovascular system. In addition, we describe how nutritional and/or disease-associated magnesium deficiency, seen in some metabolic conditions, has the potential to influence cardiac and vascular outcomes. Finally, we also examine the potential for magnesium supplements to be employed in the prevention and treatment of cardiovascular disorders and in the management of cardiometabolic health.

Modeling of minimal systems based on ATP-Zn coordination for chemically fueled self-assembly

Following nature's example, there is currently strong interest in using adenosine 5'-triphosphate (ATP) as a fuel for the self-assembly of functional materials with transient/non-equilibrium behaviours. These hold great promise for applications, e.g. in catalysis and drug delivery. In a recent seminal work [Maiti et al., Nat. Chem., 2016, 8, 725], binding of ATP to the metallosurfactant zinc hexadecyl-1,4,7-triazacyclononane ([ZnC₁₆ TACN]²⁺) was exploited to produce ATP-fueled transient vesicles. Crucial to the complex formation is the ability of ATP to bind to the metal ion. As a first step to unveil the key elements underlying this process, we investigate the interaction of ATP with Zn²⁺ and with methyl-1,4,7-triazacyclononane ([ZnCH₃ TACN]²⁺), using all-atom molecular dynamics simulations. The free energy landscape of the complex formation is sampled using well-tempered metadynamics with three collective variables, corresponding to the coordination numbers of Zn²⁺ with the oxygen atoms of the three phosphate groups. We find that the structure of the ternary complex is controlled by direct triphosphate coordination to zinc, with a minor role played by the interactions between ATP and CH₃ TACN which, however, may be important for the build-up of supramolecular assemblies.

Computational Prediction of Resistance Induced Alanine-Mutation in ATP Site of Epidermal Growth Factor Receptor

Epidermal growth factor receptor (EGFR) resistance to tyrosine kinase inhibitors can cause low survival rates in mutation-positive non-small cell lung cancer patients. It is necessary to predict new mutations in the development of more potent EGFR inhibitors since classical and rare mutations observed were known to affect the effectiveness of the therapy. Therefore, this research aimed to perform alanine mutagenesis scanning on ATP binding site residues without COSMIC data, followed by molecular dynamic simulations to determine their molecular interactions with ATP and erlotinib compared to wild-type complexes. Based on the result, eight mutations were found to cause changes in the binding energy of the ATP analogue to become more negative. These included G779A, Q791A, L792A, R841A, N842A, V843A, I853A, and D855A, which were predicted to enhance the affinity of ATP and reduce the binding ability of inhibitors with the same interaction site. Erlotinib showed more positive energy among G779A, Q791A, I853A, and D855A, due to their weaker binding energy than ATP. These four mutations could be anticipated in the development of the next inhibitor to overcome the incidence of resistance in lung cancer patients.

Magnesium(II)‐ATP Complexes in 1‐Ethyl‐3‐Methylimidazolium Acetate Solutions Characterized by 31 Mg β‐Radiation‐Detected NMR Spectroscopy

The complexation of Mg^II with adenosine 5′‐triphosphate (ATP) is omnipresent in biochemical energy conversion, but is difficult to interrogate directly. Here we use the spin‐1/2 β‐emitter ³¹Mg to study Mg^II‐ATP complexation in 1‐ethyl‐3‐methylimidazolium acetate (EMIM‐Ac) solutions using β‐radiation‐detected nuclear magnetic resonance (β‐NMR). We demonstrate that (nuclear) spin‐polarized ³¹Mg, following ion‐implantation from an accelerator beamline into EMIM‐Ac, binds to ATP within its radioactive lifetime before depolarizing. The evolution of the spectra with solute concentration indicates that the implanted ³¹Mg initially bind to the solvent acetate anions, whereafter they undergo dynamic exchange and form either a mono‐ (³¹Mg‐ATP) or di‐nuclear (³¹MgMg‐ATP) complex. The chemical shift of ³¹Mg‐ATP is observed up‐field of ³¹MgMg‐ATP, in accord with quantum chemical calculations. These observations constitute a crucial advance towards using β‐NMR to probe chemistry and biochemistry in solution.

Using Atomistic Simulations to Explore the Role of Methylation and ATP in Chemotaxis Signal Transduction

A bacterial chemotaxis mechanism is activated when nutrients bind to surface receptors. The sequence of intra- and interprotein events in this signal cascade from the receptors to the eventual molecular motors has been clearly identified. However, the atomistic details remain elusive, as in general may be expected of intraprotein signal transduction pathways, especially when fibrillar proteins are involved. We performed atomistic calculations of the methyl accepting chemoprotein (MCP)-CheA-CheW multidomain complex from Escherichia coli, simulating the methylated and unmethylated conditions in the chemoreceptors and the ATP-bound and apo conditions of the CheA. Our results indicate that these atomistic simulations, especially with one of the two force fields we tried, capture several relevant features of the downstream effects, such as the methylation favoring an intermediate structure that is more toward a dipped state and increases the chance of ATP hydrolysis. The results thus suggest the sensitivity of the model to reflect the nutrient signal response, a nontrivial validation considering the complexity of the system, encouraging even more detailed studies on the thermodynamic quantification of the effects and the identification of the signaling networks.

Catalytic activity in vitro of the human protein kinase ASK1 mutants: Experimental and molecular simulation study

Kinases have become an important class of targets for drug discovery since the milestone approval of imatinib in 2001. Although a great success has been achieved for targeting kinases with over 70 inhibitors approved by the FDA, it is inevitable that drug resistance would emerge during treatment. Thus, assessment of the kinase mutations is an essential issue for the development of the next generation inhibitors. Apoptosis signal-regulating kinase 1 (ASK1) is a crucial regulator of classical mitogen-activated protein kinase cascade that is being explored under several clinical trials as a promising target. Herein, we investigate the catalytic activity in vitro of ASK1 by constructing two mutants: M754T and H729L, from gatekeeper and αC-helix, respectively. Compared to wild type, the mutation of M754T and H729L results in a roughly 3-fold and 2-fold decrease in binding affinity experimentally. In addition, their binding modes with substrate are theoretically predicted and compared by molecular dynamics. Trajectory analyses of simulations indicate that the decrease of binding affinity should be attributed to the loss of H-bond interaction with gatekeeper methionine. Unexpectedly, the conformation of αC-helix in H729L mutant did not alter significantly during the simulations, although the putatively important H-bond with H729 is lost. These simulations showed the regulatory role of H729 in αC-helix is maintained by leucine residue through the interaction with non-polar residues around H729 site.

Free energy and kinetics of cAMP permeation through connexin26 via applied voltage and milestoning

The connexin family is a diverse group of highly regulated wide-pore channels permeable to biological signaling molecules. Despite the critical roles of connexins in mediating selective molecular signaling in health and disease, the basis of molecular permeation through these pores remains unclear. Here, we report the thermodynamics and kinetics of binding and transport of a second messenger, adenosine-3',5'-cyclophosphate (cAMP), through a connexin26 hemichannel (Cx26). First, inward and outward fluxes of cAMP molecules solvated in KCl solution were obtained from 4 μs of ± 200 mV simulations. These fluxes data yielded a single-channel permeability of cAMP and cAMP/K⁺ permeability ratio consistent with experimentally measured values. The results from voltage simulations were then compared with the potential of mean force (PMF) and the mean first passage times (MFPTs) of a single cAMP without voltage, obtained from a total of 16.5 μs of Voronoi-tessellated Markovian milestoning simulations. Both the voltage simulations and the milestoning simulations revealed two cAMP-binding sites, for which the binding constants K_D and dissociation rates k_off were computed from PMF and MFPTs. The protein dipole inside the pore produces an asymmetric PMF, reflected in unequal cAMP MFPTs in each direction once within the pore. The free energy profiles under opposite voltages were derived from the milestoning PMF and revealed the interplay between voltage and channel polarity on the total free energy. In addition, we show how these factors influence the cAMP dipole vector during permeation, and how cAMP affects the local and nonlocal pore diameter in a position-dependent manner.