Structural studies of a hyperthermophilic thymidylate kinase enzyme reveal conformational substates along the reaction coordinate

Adenylate-driven equilibration of both ribo- and deoxyribonucleotides is under magnesium control: Quantification of the Mg2+-signal

Nucleoside mono-, di- and triphosphates (NMP, NDP, and NTP) and their deoxy-counterparts (dNMP, dNDP, dNTP) are involved in energy metabolism and are the building blocks of RNA and DNA, respectively. The production of NTP and dNTP is carried out by several NMP kinases (NMPK) and NDP kinases (NDPK). All NMPKs are fully reversible and use defined Mg-free and Mg-complexed nucleotides in both directions of their reactions, with Mg2+ controlling the ratios of Mg-free and Mg-complexed reactants. Their activities are driven by adenylates produced by adenylate kinase which controls the direction of NMPK and NDPK reactions, depending on the energy status of a cell. This enzymatic machinery is localized in the cytosol, mitochondria, and plastids, i.e. compartments with high energy budgets and where (except for cytosol) RNA and DNA synthesis occur. Apparent equilibrium constants of NMPKs, based on total nucleotide contents, are [Mg2+]-dependent. This allows for an indirect estimation of internal [Mg2+], which constitutes a signal of the energetic status of a given tissue/cell/compartment. Adenylates contribute the most to this Mg2+-signal, followed by uridylates, guanylates, and cytidylates, with deoxynucleotides' contribution deemed negligible. A method to quantify the Mg2+-signal, using nucleotide datasets, is discussed.

The archaeal family 3 polyphosphate kinase reveals a function of polyphosphate as energy buffer under low energy charge

Inorganic polyphosphate, a linear polymer of orthophosphate residues linked by phosphoanhydride bonds, occurs in all three domains of life and plays a diverse and prominent role in metabolism and cellular regulation. While the polyphosphate metabolism and its physiological significance have been well studied in bacteria and eukaryotes including human, there are only few studies in archaea available so far. In Crenarchaeota including members of Sulfolobaceae , the presence of polyphosphate and degradation via exopolyphosphatase has been reported and there is some evidence for a functional role in metal ion chelation, biofilm formation, adhesion and motility, however, the nature of the crenarchaeal polyphosphate kinase is still unknown. Here we used the crenarchaeal model organism Sulfolobus acidocaldarius to study the enzymes involved in polyphosphate synthesis. The two genes annotated as thymidylate kinase ( saci_2019 and saci_2020 ), localized downstream of the exopolyphosphatase, were identified as the missing polyphosphate kinase in S. acidocaldarius ( Sa PPK3). Thymidylate kinase activity was confirmed for Saci_0893. Notably Saci_2020 showed no polyphosphate kinase activity on its own but served as regulatory subunit (rPPK3) and was able to enhance polyphosphate kinase activity of the catalytically active subunit Saci_2019 (cPPK3). Heteromeric polyphosphate kinase activity is reversible and shows a clear preference for polyP-dependent nucleotide kinase activity, i.e. polyP-dependent formation of ATP from ADP (12.4 U/mg) and to a lower extent of GDP to GTP whereas AMP does not serve as substrate. PPK activity in the direction of ATP-dependent polyP synthesis is rather low (0.25 U/mg); GTP was not used as phosphoryl donor. A combined experimental modelling approach using quantitative 31 P NMR allowed to follow the reversible enzyme reaction for both ATP and polyP synthesis. PolyP synthesis was only observed when the ATP/ADP ratio was kept high, using an ATP recycling system. In absence of such a recycling system, all incubations with polyP and PPK would reach an equilibrium state with an ATP/ADP ratio between 3 and 4, independent of the initial conditions. Structural and sequence comparisons as well as phylogenetic analysis reveal that the S. acidocaldarius PPK is a member of a new PPK family, named PPK3, within the thymidylate kinase family of the P-loop kinase superfamily, clearly separated from PPK2. Our studies show that polyP, in addition to its function as phosphate storage, has a special importance for the energy homeostasis of S. acidocaldarius and due to its reversibility serves as energy buffer under low energy charge enabling a quick response to changes in cellular demand.

Identification of mycobacterial Thymidylate kinase inhibitors: a comprehensive pharmacophore, machine learning, molecular docking, and molecular dynamics simulation studies

Thymidylate kinase (TMK) is a pivotal enzyme in Mycobacterium tuberculosis (Mtb), crucial for phosphorylating thymidine monophosphate (dTMP) to thymidine diphosphate (dTDP), thereby playing a critical role in DNA biosynthesis. Dysregulation or inhibition of TMK activity disrupts DNA replication and cell division, making it an attractive target for anti-tuberculosis drug development. In this study, the statistically validated pharmacophore mode was developed from a set of known TMK inhibitors. Further, the robust pharmacophore was considered for screening the Enamine database. The chemical space was reduced through multiple molecular docking approaches, pharmacokinetics, and absolute binding energy estimation. Two different molecular docking algorithms favor the strong binding affinity of the proposed molecules towards TMK. Machine learning-based absolute binding energy also showed the potentiality of the proposed molecules. The binding interactions analysis exposed the strong binding affinity between the proposed molecules and active site amino residues of TMK. Several statistical parameters from all atoms MD simulation explained the stability between proposed molecules and TMK in the dynamic states. The MM-GBSA approach also found a strong binding affinity for each proposed molecule. Therefore, the proposed molecules might be crucial TMK inhibitors for managing Mtb inhibition subjected to in vitro/in vivo validations.

Multifunctional in vitro, in silico and DFT analyses on antimicrobial BagremycinA biosynthesized by Micromonospora chokoriensis CR3 from Hieracium canadense

Among the actinomycetes in the rare genera, Micromonospora is of great interest since it has been shown to produce novel therapeutic compounds. Particular emphasis is now on its isolation from plants since its population from soil has been extensively explored. The strain CR3 was isolated as an endophyte from the roots of Hieracium canadense, and it was identified as Micromonospora chokoriensis through 16S gene sequencing and phylogenetic analysis. The in-vitro analysis of its extract revealed it to be active against the clinical isolates of methicillin-resistant Staphylococcus aureus (MRSA) and Candida tropicalis (15 mm). No bioactivity was observed against Gram-negative bacteria, Escherichia coli ATCC 25922, and Klebsiella pneumoniae ATCC 706003. The Micromonospora chokoriensis CR3 extract was also analyzed through the HPLC-DAD-UV-VIS resident database, and it gave a maximum match factor of 997.334 with the specialized metabolite BagremycinA (BagA). The in-silico analysis indicated that BagA strongly interacted with the active site residues of the sterol 14-α demethylase and thymidylate kinase enzymes, with the lowest binding energies of - 9.7 and - 8.3 kcal/mol, respectively. Furthermore, the normal mode analysis indicated that the interaction between these proteins and BagA was stable. The DFT quantum chemical properties depicted BagA to be reasonably reactive with a HOMO-LUMO gap of (ΔE) of 4.390 eV. BagA also passed the drug-likeness test with a synthetic accessibility score of 2.06, whereas Protox-II classified it as a class V toxicity compound with high LD50 of 2644 mg/kg. The current study reports an endophytic actinomycete, M. chokoriensis, associated with H. canadense producing the bioactive metabolite BagA with promising antimicrobial activity, which can be further modified and developed into a safe antimicrobial drug.

Combined Plasma and Urinary Metabolomics Uncover Metabolic Perturbations Associated with Severe Respiratory Syncytial Viral Infection and Future Development of Asthma in Infant Patients

A large percentage of infants develop viral bronchiolitis needing medical intervention and often develop further airway disease such as asthma. To characterize metabolic perturbations in acute respiratory syncytial viral (RSV) bronchiolitis, we compared metabolomic profiles of moderate and severe RSV patients versus sedation controls. RSV patients were classified as moderate or severe based on the need for invasive mechanical ventilation. Whole blood and urine samples were collected at two time points (baseline and 72 h). Plasma and urinary metabolites were extracted in cold methanol and analyzed by liquid chromatography coupled with tandem mass spectrometry (LC-MS/MS), and data from the two biofluids were combined for multivariate data analysis. Metabolite profiles were clustered according to severity, characterized by unique metabolic changes in both plasma and urine. Plasma metabolites that correlated with severity included intermediates in the sialic acid biosynthesis, while urinary metabolites included citrate as well as multiple nucleotides. Furthermore, metabolomic profiles were predictive of future development of asthma, with urinary metabolites exhibiting higher predictive power than plasma. These metabolites may offer unique insights into the pathology of RSV bronchiolitis and may be useful in identifying patients at risk for developing asthma.

Modeling a unit cell: crystallographic refinement procedure using the biomolecular MD simulation platform Amber

A procedure has been developed for the refinement of crystallographic protein structures based on the biomolecular simulation program Amber. The procedure constructs a model representing a crystal unit cell, which generally contains multiple protein molecules and is fully hydrated with TIP3P water. Periodic boundary conditions are applied to the cell in order to emulate the crystal lattice. The refinement is conducted in the form of a specially designed short molecular-dynamics run controlled by the Amber ff14SB force field and the maximum-likelihood potential that encodes the structure-factor-based restraints. The new Amber-based refinement procedure has been tested on a set of 84 protein structures. In most cases, the new procedure led to appreciably lower R free values compared with those reported in the original PDB depositions or obtained by means of the industry-standard phenix.refine program. In particular, the new method has the edge in refining low-accuracy scrambled models. It has also been successful in refining a number of molecular-replacement models, including one with an r.m.s.d. of 2.15 Å. In addition, Amber-refined structures consistently show superior MolProbity scores. The new approach offers a highly realistic representation of protein-protein interactions in the crystal, as well as of protein-water interactions. It also offers a realistic representation of protein crystal dynamics (akin to ensemble-refinement schemes). Importantly, the method fully utilizes the information from the available diffraction data, while relying on state-of-the-art molecular-dynamics modeling to assist with those elements of the structure that do not diffract well (for example mobile loops or side chains). Finally, it should be noted that the protocol employs no tunable parameters, and the calculations can be conducted in a matter of several hours on desktop computers equipped with graphical processing units or using a designated web service.

Insights into product release dynamics through structural analyses of thymidylate kinase

Several studies on enzyme catalysis have pointed out that the product release event could be a rate limiting step. In this study, we have compared the release event of two products, Adenosine di-phosphate (ADP) and Thymidine di-phosphate (TDP) from the active-site of human and Thermus thermophilus thymidine mono-phosphate kinase (TMPK), referred to as hTMPK and ttTMPK, respectively. TMPK catalyses the conversion of Thymidine mono-phosphate (TMP) to TDP using ATP as phosphoryl donor in the presence of Mg²⁺ ion. Most of the earlier studies on this enzyme have focused on understanding substrate binding and catalysis, but the critical product release event remains elusive. Competitive binding experiments of the substrates and the products using ttTMPK apo crystals have indicated that the substrate (TMP) can replace the bound product (TDP), even in the presence of an ADP molecule. Further, the existing random accelerated molecular dynamics (RAMD) simulation program was modified to study the release of both the products simultaneously from the active site. The RAMD simulations on product-bound structures of both ttTMPK and hTMPK, revealed that while several exit patterns of the products are permissible, the sequential exit mode is the most preferred pattern for both ttTMPK and hTMPK enzymes. Additionally, the product release from the hTMPK was found to be faster and more directional as compared to ttTMPK. Structural investigation revealed that the critical changes in the residue composition in the LID-region of ttTMPK and hTMPK have an effect on the product release and can be attributed to the observed differences during product release event. Understanding of these dissimilarities is of considerable utility in designing potent inhibitors or prodrugs that can distinguish between eukaryotic and prokaryotic homologues of thymidylate kinase.

The Binding Mechanism Between Inositol Phosphate (InsP) and the Jasmonate Receptor Complex: A Computational Study

Jasmonates are critical plant hormones, mediating stress response in plants and regulating plant growth and development. The jasmonate receptor is a multi-component complex, composed of Arabidopsis SKP-LIKE PROTEIN1 (ASK1), CORONATINE INSENSITIVE 1 (COI1), inositol phosphate (InsP), and jasmonate ZIM-domain protein (JAZ). COI1 acts as multi-component signaling hub that binds with each component. InsP is suggested to play important roles in the hormone perception. How InsP binds with COI1 and the structural changes in COI1 upon binding with InsP, JA-Ile, and JAZ are not well understood. In this study, we integrated multiple computational methods, such as molecular docking, molecular dynamics simulations, residue interaction network analysis and binding free energy calculation, to explore the effect of InsP on the dynamic behavior of COI1 and the recognition mechanism of each component of the jasmonate receptor complex. We found that upon binding with InsP, JA-Ile, and JAZ1, the structure of COI1 becomes more compact. The binding of InsP with COI1 stabilizes the conformation of COI1 and promotes the binding between JA-Ile or JAZ1 and COI1. Analysis of the network parameters led to the identification of some hub nodes in this network, including Met88, His118, Arg120, Arg121, Arg346, Tyr382, Arg409, Trp467, and Lys492. The structural and dynamic details will be helpful for understanding the recognition mechanism of each component and the discovery and design of novel jasmonate signaling pathway modulators.

Structural and functional roles of dynamically correlated residues in thymidylate kinase

Thymidylate kinase is an important enzyme in DNA synthesis. It catalyzes the conversion of thymidine monophosphate to thymidine diphosphate, with ATP as the preferred phosphoryl donor, in the presence of Mg²⁺. In this study, the dynamics of the active site and the communication paths between the substrates, ATP and TMP, are reported for thymidylate kinase from Thermus thermophilus. Conformational changes upon ligand binding and the path for communication between the substrates and the protein are important in understanding the catalytic mechanism of the enzyme. High-resolution X-ray crystal structures of thymidylate kinase in apo and ligand-bound states were solved. This is the first report of structures of binary and ternary complexes of thymidylate kinase with its natural substrates ATP and ATP-TMP, respectively. Distinct conformations of the active-site residues, the P-loop and the LID region observed in the apo and ligand-bound structures revealed that their concerted motion is required for the binding and proper positioning of the substrate TMP. Structural analyses provide an insight into the mode of substrate binding at the active site. The residues involved in communication between the substrates were identified through network analysis using molecular-dynamics simulations. The residues identified showed high sequence conservation across species. Biochemical analyses show that mutations of these residues either resulted in a loss of activity or affected the thermal stability of the protein. Further, molecular-dynamics analyses of mutants suggest that the proper positioning of TMP is important for catalysis. These data also provide an insight into the phosphoryl-transfer mechanism.