Studying the phosphoryl transfer mechanism of the E. coli phosphofructokinase-2: from X-ray structure to quantum mechanics/molecular mechanics simulations

Structural characterization of a novel cyclic 2,3-diphosphoglycerate synthetase involved in extremolyte production in the archaeon Methanothermus fervidus

The enzyme cyclic di-phosphoglycerate synthetase that is involved in the production of the osmolyte cyclic 2,3-diphosphoglycerate has been studied both biochemically and structurally. Cyclic 2,3-diphosphoglycerate is found exclusively in the hyperthermophilic archaeal methanogens, such as Methanothermus fervidus, Methanopyrus kandleri, and Methanothermobacter thermoautotrophicus. Its presence increases the thermostability of archaeal proteins and protects the DNA against oxidative damage caused by hydroxyl radicals. The cyclic 2,3-diphosphoglycerate synthetase enzyme has been crystallized and its structure solved to 1.7 Å resolution by experimental phasing. It has also been crystallized in complex with its substrate 2,3 diphosphoglycerate and the co-factor ADP and this structure has been solved to 2.2 Å resolution. The enzyme structure has two domains, the core domain shares some structural similarity with other NTP-dependent enzymes. A significant proportion of the structure, including a 127 amino acid N-terminal domain, has no structural similarity to other known enzyme structures. The structure of the complex shows a large conformational change that occurs in the enzyme during catalytic turnover. The reaction involves the transfer of the γ-phosphate group from ATP to the substrate 2,3 -diphosphoglycerate and the subsequent SN2 attack to form a phosphoanhydride. This results in the production of the unusual extremolyte cyclic 2,3 -diphosphoglycerate which has important industrial applications.

Activation of p53: How phosphorylated Ser15 triggers sequential phosphorylation of p53 at Thr18 by CK1δ

The N-terminal transactivation domain (TAD) of p53 is a disordered region with multiple phosphorylation sites. Phosphorylation at Thr18 is crucial for the release of p53 from its negative regulator, MDM2. In stressed cells, CK1δ is responsible for phosphorylating Thr18, but requires Ser15 to be phosphorylated. To understand the mechanistic underpinnings of this sequential phosphorylation, molecular modeling and molecular dynamics simulation studies of these phosphorylation events were carried out. Our models suggest that a positively charged region on CK1δ near the adenosine triphosphate (ATP) binding pocket, which is conserved across species, sequesters the negatively charged pSer15, thereby constraining the positioning of the rest of the peptide, such that the side chain of Thr18 is positioned close to the γ-phosphate of ATP. Furthermore, our studies show that the phosphorylated p53 TAD1 (p53pSer15) peptide binds more strongly to CK1δ than does p53. p53 adopts a helical structure when bound to CK1δ, which is lost upon phosphorylation at Ser15, thus gaining higher flexibility and ability to morph into the binding site. We propose that upon phosphorylation at Ser15 the p53 TAD1 peptide binds to CK1δ through an electrostatically driven induced fit mechanism resulting in a flanking fuzzy complex.

3D matrix promotes cell dedifferentiation into colorectal cancer stem cells via integrin/cytoskeleton/glycolysis signaling

The potential for tumor occurrence triggered by cancer stem cells (CSCs) has emerged as a significant challenge for human colorectal cancer therapy. However, the underlying mechanism of CSC development remains controversial. Our study provided evidence that the bulk of tumor cells could dedifferentiate to CSCs and reacquire CSCs-like phenotypes if cultured in the presence of extracellular matrix reagents, such as Matrigel and fibrin gels. In these 3D gels, CD133^- colorectal cancer cells can regain tumorigenic potential and stem-like phenotypes. Mechanistically, the 3D extracellular matrix could mediate cytoskeletal F-actin bundling through biomechanical force associated receptors integrin β1 (ITGB1), contributing to the release of E3 ligase Tripartite motif protein 11 (TRIM11) from cytoskeleton and degradation of the glycolytic rate-limiting enzyme phosphofructokinase (PFK). Consequently, PFK inhibition resulted in enhanced glycolysis and upregulation of hypoxia-inducible factor 1 (HIF1α), thereby promoting the reprogramming of stem cell transcription factors and facilitating tumor progression in patients. This study provided novel insights into the role of the extracellular matrix in the regulation of CSC dedifferentiation in a cytoskeleton/glycolysis-dependent manner.

Catalytic Mechanism of ATP Hydrolysis in the ATPase Domain of Human DNA Topoisomerase IIα

Human DNA topoisomerase IIα is a biological nanomachine that regulates the topological changes of the DNA molecule and is considered a prime target for anticancer drugs. Despite intensive research, many atomic details about its mechanism of action remain unknown. We investigated the ATPase domain, a segment of the human DNA topoisomerase IIα, using all-atom molecular simulations, multiscale quantum mechanics/molecular mechanics (QM/MM) calculations, and a point mutation study. The results suggested that the binding of ATP affects the overall dynamics of the ATPase dimer. Reaction modeling revealed that ATP hydrolysis favors the dissociative substrate-assisted reaction mechanism with the catalytic Glu87 serving to properly position and polarize the lytic water molecule. The point mutation study complemented our computational results, demonstrating that Lys378, part of the important QTK loop, acts as a stabilizing residue. The work aims to pave the way to a deeper understanding of these important molecular motors and to advance the development of new therapeutics.

Molecular Basis of Sulfosugar Selectivity in Sulfoglycolysis

The sulfosugar sulfoquinovose (SQ) is produced by essentially all photosynthetic organisms on Earth and is metabolized by bacteria through the process of sulfoglycolysis. The sulfoglycolytic Embden-Meyerhof-Parnas pathway metabolizes SQ to produce dihydroxyacetone phosphate and sulfolactaldehyde and is analogous to the classical Embden-Meyerhof-Parnas glycolysis pathway for the metabolism of glucose-6-phosphate, though the former only provides one C3 fragment to central metabolism, with excretion of the other C3 fragment as dihydroxypropanesulfonate. Here, we report a comprehensive structural and biochemical analysis of the three core steps of sulfoglycolysis catalyzed by SQ isomerase, sulfofructose (SF) kinase, and sulfofructose-1-phosphate (SFP) aldolase. Our data show that despite the superficial similarity of this pathway to glycolysis, the sulfoglycolytic enzymes are specific for SQ metabolites and are not catalytically active on related metabolites from glycolytic pathways. This observation is rationalized by three-dimensional structures of each enzyme, which reveal the presence of conserved sulfonate binding pockets. We show that SQ isomerase acts preferentially on the β-anomer of SQ and reversibly produces both SF and sulforhamnose (SR), a previously unknown sugar that acts as a derepressor for the transcriptional repressor CsqR that regulates SQ-utilization. We also demonstrate that SF kinase is a key regulatory enzyme for the pathway that experiences complex modulation by the metabolites SQ, SLA, AMP, ADP, ATP, F6P, FBP, PEP, DHAP, and citrate, and we show that SFP aldolase reversibly synthesizes SFP. This body of work provides fresh insights into the mechanism, specificity, and regulation of sulfoglycolysis and has important implications for understanding how this biochemistry interfaces with central metabolism in prokaryotes to process this major repository of biogeochemical sulfur.

Temperature effect on water dynamics in tetramer phosphofructokinase matrix and the super-arrhenius respiration rate

Advances in understanding the temperature effect on water dynamics in cellular respiration are important for the modeling of integrated energy processes and metabolic rates. For more than half a century, experimental studies have contributed to the understanding of the catalytic role of water in respiration combustion, yet the detailed water dynamics remains elusive. We combine a super-Arrhenius model that links the temperature-dependent exponential growth rate of a population of plant cells to respiration, and an experiment on isotope labeled ¹⁸O₂ uptake to H₂¹⁸O transport role and to a rate-limiting step of cellular respiration. We use Phosphofructokinase (PFK-1) as a prototype because this enzyme is known to be a pacemaker (a rate-limiting enzyme) in the glycolysis process of respiration. The characterization shows that PFK-1 water matrix dynamics are crucial for examining how respiration (PFK-1 tetramer complex breathing) rates respond to temperature change through a water and nano-channel network created by the enzyme folding surfaces, at both short and long (evolutionary) timescales. We not only reveal the nano-channel water network of PFK-1 tetramer hydration topography but also clarify how temperature drives the underlying respiration rates by mapping the channels of water diffusion with distinct dynamics in space and time. The results show that the PFK-1 assembly tetramer possesses a sustainable capacity in the regulation of the water network toward metabolic rates. The implications and limitations of the reciprocal-activation-reciprocal-temperature relationship for interpreting PFK-1 tetramer mechanisms are briefly discussed.

Molecular Insights into the Trapping Effect of Ca2+ in Protein Kinase A: A Molecular Dynamics Study

Protein kinase A has become a model system for the study of kinases, and therefore, a comprehensive understanding of the underlying molecular mechanisms in its catalytic cycle is of crucial importance. One of the aspects that has received recent attention is the role that metal cofactors play in the catalytic cycle. Although Mg²⁺ is the well-known physiological ion used by protein kinases, Ca²⁺ ions can also assist the phosphoryl transfer reaction but with lower catalytic activities. This inhibitory effect has been attributed to the ability of Ca²⁺ to trap the reaction products at the active site, and it has been proposed as a possible regulatory mechanism of the enzyme. Thus, in order to get a clearer understanding of these molecular events, computational simulations in the product state of PKA, in the presence of Mg²⁺ and Ca²⁺ ions, were performed through molecular dynamics (MD). Different protonation states of the active site were considered in order to model the different mechanistic pathways that have been proposed. Our results show that different protonation states of the phosphorylated serine residue at the peptide substrate (pSer21), as well as the protonation state of residue Asp166, can have a marked influence on the flexibility of regions surrounding the active site. This is the case of the glycine-rich loop, a structural motif that is directly involved in the release of the products from the PKA active site. MD simulations were capable to reproduce the crystallographic conformations but also showed other conformations not previously reported in the crystal structures that may be involved in enhancing the affinity of pSP20 to PKA in the presence of Ca²⁺. Hydrogen bonding interactions at the PKA-pSP20 interface were influenced whether by the protonation state of the active site or by the metal cofactor used by the enzyme. Altogether, our results provide molecular aspects into the inhibitory mechanism of Ca²⁺ in PKA and suggest which is the most probable protonation state of the phosphorylated product at the active site.

Design and SAR Analysis of Covalent Inhibitors Driven by Hybrid QM/MM Simulations

Quantum mechanics/molecular mechanics (QM/MM) hybrid technique is emerging as a reliable computational method to investigate and characterize chemical reactions occurring in enzymes. From a drug discovery perspective, a thorough understanding of enzyme catalysis appears pivotal to assist the design of inhibitors able to covalently bind one of the residues belonging to the enzyme catalytic machinery. Thanks to the current advances in computer power, and the availability of more efficient algorithms for QM-based simulations, the use of QM/MM methodology is becoming a viable option in the field of covalent inhibitor design. In the present review, we summarized our experience in the field of QM/MM simulations applied to drug design problems which involved the optimization of agents working on two well-known drug targets, namely fatty acid amide hydrolase (FAAH) and epidermal growth factor receptor (EGFR). In this context, QM/MM simulations gave valuable information in terms of geometry (i.e., of transition states and metastable intermediates) and reaction energetics that allowed to correctly predict inhibitor binding orientation and substituent effect on enzyme inhibition. What is more, enzyme reaction modelling with QM/MM provided insights that were translated into the synthesis of new covalent inhibitor featured by a unique combination of intrinsic reactivity, on-target activity, and selectivity.

A Transient and Flexible Cation-π Interaction Promotes Hydrolysis of Nucleic Acids in DNA and RNA Nucleases

Metal-dependent DNA and RNA nucleases are enzymes that cleave nucleic acids with great efficiency and precision. These enzyme-mediated hydrolytic reactions are fundamental for the replication, repair, and storage of genetic information within the cell. Here, extensive classical and quantum-based free-energy molecular simulations show that a cation-π interaction is transiently formed in situ at the metal core of Bacteriophage-λ Exonuclease (Exo-λ), during catalysis. This noncovalent interaction (Lys131-Tyr154) triggers nucleophile activation for nucleotide excision. Then, our simulations also show the oscillatory dynamics and swinging of the newly formed cation-π dyad, whose conformational change may favor proton release from the cationic Lys131 to the bulk solution, thus restoring the precatalytic protonation state in Exo-λ. Altogether, we report on the novel mechanistic character of cation-π interactions for catalysis. Structural and bioinformatic analyses support that flexible orientation and transient formation of mobile cation-π interactions may represent a common catalytic strategy to promote nucleic acid hydrolysis in DNA and RNA nucleases.