Biochemistry Biochemistry最新文献

Transketolases (TKs) are key enzymes of the pentose phosphate pathway, regulating several other critical pathways in cells. Considering their metabolic importance, TKs are expected to be conserved throughout evolution. However, Tittmann et al. (J Biol Chem, 2010, 285(41): 31559–31570) demonstrated that Homo sapiens TK (hsTK) possesses several structural and kinetic differences compared to bacterial TKs. Here, we study 14 TKs from pathogenic bacteria, fungi, and parasites and compare them with hsTK using biochemical, bioinformatic, and structural approaches. For this purpose, six new TK structures are solved by X-ray crystallography, including the TK of Plasmodium falciparum. All of these TKs have the same general fold as bacterial TKs. This comparative study shows that hsTK greatly differs from TKs from pathogens in terms of enzymatic activity, spatial positions of the active site, and monomer–monomer interface residues. An ubiquitous structural pattern is identified in all TKs as a six-residue histidyl crown around the TK cofactor (thiamine pyrophosphate), except for hsTK containing only five residues in the crown. Residue mapping of the monomer–monomer interface and the active site reveals that hsTK contains more unique residues than other TKs. From an evolutionary standpoint, TKs from animals (including H. sapiens) and Schistosoma sp. belong to a distinct structural group from TKs of bacteria, plants, fungi, and parasites, mostly based on a different linker between domains, raising hypotheses regarding evolution and regulation.

Global substitution of leucine for analogues containing CH₂F instead of methyl groups delivers proteins with multiple sites for monitoring by ¹⁹F nuclear magnetic resonance (NMR) spectroscopy. The 19 kDa Escherichia coli peptidyl–prolyl cis–trans isomerase B (PpiB) was prepared with uniform high-level substitution of leucine by (2S,4S)-5-fluoroleucine, (2S,4R)-5-fluoroleucine, or 5,5′-difluoroleucine. The stability of the samples toward thermal denaturation was little altered compared to the wild-type protein. ¹⁹F nuclear magnetic resonance (NMR) spectra showed large chemical shift dispersions between 6 and 17 ppm. The ¹⁹F chemical shifts correlate with the three-bond ¹H–¹⁹F couplings (³J_HF), providing the first experimental verification of the γ-gauche effect predicted by [Feeney, J. J. Am. Chem. Soc. 1996, 118, 8700–8706] and establishing the effect as the predominant determinant of the ¹⁹F chemical shifts of CH₂F groups. Individual CH₂F groups can be confined to single rotameric states by the protein environment, but most CH₂F groups exchange between different rotamers at a rate that is fast on the NMR chemical shift scale. Interactions between fluorine atoms in 5,5′-difluoroleucine bias the CH₂F rotamers in agreement with results obtained previously for 1,3-difluoropropane. The sensitivity of the ¹⁹F chemical shift to the rotameric state of the CH₂F groups potentially renders them particularly sensitive for detecting allosteric effects.

Nitric oxide synthase (NOS) in mammals is a family of multidomain proteins in which interdomain electron transfer (IET) is controlled by domain–domain interactions. Calmodulin (CaM) binds to the canonical CaM-binding site in the linker region between the FMN and heme domains of NOS and allows tethered FMN domain motions, enabling an intersubunit FMN-heme IET in the output state for NO production. Our previous cross-linking mass spectrometric (XL MS) results demonstrated site-specific protein dynamics in the CaM-responsive regions of rat neuronal NOS (nNOS) reductase construct, a monomeric protein [Jiang et al., Biochemistry, 2023, 62, 2232–2237]. In this work, we have extended our combined approach of XL MS structural mapping and AlphaFold structural prediction to examine the homodimeric nNOS oxygenase/FMN (oxyFMN) construct, an established model of the NOS output state. We employed parallel reaction monitoring (PRM) based quantitative XL MS (qXL MS) to assess the CaM-induced changes in interdomain dynamics and interactions. Intersubunit cross-links were identified by mapping the cross-links onto top AlphaFold structural models, which was complemented by comparing their relative abundances in the cross-linked dimeric and monomeric bands. Furthermore, contrasting the CaM-free and CaM-bound nNOS samples shows that CaM enables the formation of the intersubunit FMN-heme docking complex and that CaM binding induces extensive, allosteric conformational changes across the NOS regions. Moreover, the observed cross-links sites specifically respond to changes in ionic strength. This indicates that interdomain salt bridges are responsible for stabilizing and orienting the output state for efficient FMN-heme IET. Taken together, our targeted qXL MS results have revealed that CaM and ionic strength modulate specific dynamic changes in the CaM/FMN/heme complexes, particularly in the context of intersubunit interdomain FMN-heme interactions.

Exiguobacterium sibiricum rhodopsin (ESR) functions as a light-driven proton pump utilizing Lys96 for proton uptake and maintaining its activity over a wide pH range. Using a combination of methodologies including the linear Poisson–Boltzmann equation and a quantum mechanical/molecular mechanical approach with a polarizable continuum model, we explore the microscopic mechanisms underlying its pumping activity. Lys96, the primary proton uptake site, remains deprotonated owing to the loss of solvation in the ESR protein environment. Asp85, serving as a proton acceptor group for Lys96, does not form a low-barrier H-bond with His57. Instead, deprotonated Asp85 forms a salt-bridge with protonated His57, and the proton is predominantly located at the His57 moiety. Glu214, the only acidic residue at the end of the H-bond network exhibits a pK_a value of ∼6, slightly elevated due to solvation loss. It seems likely that the H-bond network [Asp85···His57···H₂O···Glu214] serves as a proton-conducting pathway toward the protein bulk surface.

DNA gyrases catalyze negative supercoiling of DNA, are essential for bacterial DNA replication, transcription, and recombination, and are important antibacterial targets in multiple pathogens, including Mycobacterium tuberculosis, which in 2021 caused >1.5 million deaths worldwide. DNA gyrase is a tetrameric (A₂B₂) protein formed from two subunit types: gyrase A (GyrA) carries the breakage-reunion active site, whereas gyrase B (GyrB) catalyzes ATP hydrolysis required for energy transduction and DNA translocation. The GyrB ATPase domains dimerize in the presence of ATP to trap the translocated DNA (T-DNA) segment as a first step in strand passage, for which hydrolysis of one of the two ATPs and release of the resulting inorganic phosphate is rate-limiting. Here, dynamical-nonequilibrium molecular dynamics (D-NEMD) simulations of the dimeric 43 kDa N-terminal fragment of M. tuberculosis GyrB show how events at the ATPase site (dissociation/hydrolysis of bound nucleotides) are propagated through communication pathways to other functionally important regions of the GyrB ATPase domain. Specifically, our simulations identify two distinct pathways that respectively connect the GyrB ATPase site to the corynebacteria-specific C-loop, thought to interact with GyrA prior to DNA capture, and to the C-terminus of the GyrB transduction domain, which in turn contacts the C-terminal GyrB topoisomerase-primase (TOPRIM) domain responsible for interactions with GyrA and the centrally bound G-segment DNA. The connection between the ATPase site and the C-loop of dimeric GyrB is consistent with the unusual properties of M. tuberculosis DNA gyrase relative to those from other bacterial species.

Proteins produced with leucine analogues, where CH₂F groups substitute specific methyl groups, can readily be probed by ¹⁹F NMR spectroscopy. As CF and CH groups are similar in hydrophobicity and size, fluorinated leucines are expected to cause minimal structural perturbation, but the impact of fluorine on the rotational freedom of CH₂F groups is unclear. We present high-resolution crystal structures of Escherichia coli peptidyl-prolyl cis–trans isomerase B (PpiB) prepared with uniform high-level substitution of leucine by (2S,4S)-5-fluoroleucine, (2S,4R)-5-fluoroleucine, or 5,5′-difluoroleucine. Apart from the fluorinated leucine residues, the structures show complete structural conservation of the protein backbone and the amino acid side chains except for a single isoleucine side chain located next to a fluorine atom in the hydrophobic core of the protein. The carbon skeletons of the fluorinated leucine side chains are also mostly conserved. The CH₂F groups show a strong preference for staggered rotamers and often appear locked into single rotamers. Substitution of leucine CH₃ groups for CH₂F groups is thus readily tolerated in the three-dimensional (3D) structure of a protein, and the rotation of CH₂F groups can be halted at cryogenic temperatures.

Allostery is a fundamental mechanism driving biomolecular processes that holds significant therapeutic concern. Our study rigorously investigates how two distinct machine-learning algorithms uniquely classify two already close-to-active DFG-in states of TAK1, differing just by the presence or absence of its allosteric activator TAB1, from an ensemble mixture of conformations (obtained from 2.4 μs molecular dynamics (MD) simulations). The novelty, however, lies in understanding the deeper algorithmic potentials to systematically derive a diverse set of differential residue connectivity features that reconstruct the essential mechanistic architecture for TAK1-TAB1 allostery in such a close-to-active biochemical scenario. While the recursive, random forest-based workflow displays the potential of conducting discretized, hierarchical derivation of allosteric features, a multilayer perceptron-based approach gains considerable efficacy in revealing fluid connected patterns of features when hybridized with mutual information scoring. Interestingly, both pipelines benchmark similar directions of functional conformational changes for TAK1′s activation. The findings significantly advance the depth of mechanistic understanding by highlighting crucial activation signatures along a directed C-lobe → activation loop → ATP pocket channel of information flow, including (1) the αF-αE biterminal alignments and (2) the “catalytic” drift of the activation loop toward kinase active site. Besides, some novel allosteric hotspots (K253, Y206, N189, etc.) are further recognized as TAB1 sensors, transducers, and responders, including a benchmark E70 mutation site, precisely mapping the important structural segments for sequential allosteric execution. Hence, our work demonstrates how to navigate through greater structural depths and dimensions of dynamic allosteric machineries just by leveraging standard ML methods in suitable streamlined workflows adaptive to the specific system and objectives.

Lysine specific demethylase-1 (LSD1) serves as a regulator of transcription and represents a promising epigenetic target for anticancer treatment. LSD1 inhibitors are in clinical trials for the treatment of Ewing’s sarcoma (EWS), acute myeloid leukemia, and small cell lung cancer, and the development of robust inhibitors requires accurate methods for probing demethylation, potency, and selectivity. Here, the inhibition kinetics on the H3K4me2 peptide and nucleosome substrates was examined, comparing the rates of demethylation in the presence of reversible [CC-90011 (PD) and SP-2577 (SD)] and irreversible [ORY-1001 (ID) and tranylcypromine (TCP)] inhibitors. Inhibitors were also subject to viability studies in three human cell lines and Western blot assays to monitor H3K4me2 nucleosome levels in EWS (TC-32) cells, enabling a correlation of drug potency, inhibition in vitro, and cell-based studies. For example, SP-2577, a drug in clinical trials for EWS, inhibits activity on small peptide substrates (K_i = 60 ± 20 nM) using an indirect coupled assay but does not inhibit demethylation on H3K4me2 peptides or nucleosomes using direct Western blot approaches. In addition, the drug has no effect on H3K4me2 levels in TC-32 cells. These data show that SP-2577 is not an LSD1 enzyme inhibitor, although the drug may function independent of demethylation due to its cytotoxic selectivity in TC-32 cells. Taken together, this work highlights the pitfalls of using coupled assays to ascribe a drug’s mode of action, emphasizes the use of physiologically relevant substrates in epigenetic drug targeting strategies, and provides insight into the development of substrate-selective inhibitors of LSD1.

PGM1-linked congenital disorder of glycosylation (PGM1-CDG) is an autosomal recessive disease characterized by several phenotypes, some of which are life-threatening. Research focusing on the disease-related variants of the α-D-phosphoglucomutase 1 (PGM1) protein has shown that several are insoluble in vitro and expressed at low levels in patient fibroblasts. Due to these observations, we hypothesized that some disease-linked PGM1 protein variants are structurally destabilized and subject to protein quality control (PQC) and rapid intracellular degradation. Employing yeast-based assays, we show that a disease-associated human variant, PGM1 L516P, is insoluble, inactive, and highly susceptible to ubiquitylation and rapid degradation by the proteasome. In addition, we show that PGM1 L516P forms aggregates in S. cerevisiae and that both the aggregation pattern and the abundance of PGM1 L516P are chaperone-dependent. Finally, using computational methods, we perform saturation mutagenesis to assess the impact of all possible single residue substitutions in the PGM1 protein. These analyses identify numerous missense variants with predicted detrimental effects on protein function and stability. We suggest that many disease-linked PGM1 variants are subject to PQC-linked degradation and that our in silico site-saturated data set may assist in the mechanistic interpretation of PGM1 variants.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) frameshift stimulatory element (FSE) is necessary for programmed −1 ribosomal frameshifting (−1 PRF) and optimized viral efficacy. The FSE has an abundance of context-dependent alternate conformations, but two of the structures most crucial to −1 PRF are an attenuator hairpin and a three-stem H-type pseudoknot structure. A crystal structure of the pseudoknot alone features three RNA stems in a helically stacked linear structure, whereas a 6.9 Å cryo-EM structure including the upstream heptameric slippery site resulted in a bend between two stems. Our previous research alluded to an extended upstream multibranch loop that includes both the attenuator hairpin and the slippery site–a conformation not previously modeled. We aim to provide further context to the SARS-CoV-2 FSE via computational and medium resolution cryo-EM approaches, by presenting a 6.1 Å cryo-EM structure featuring a linear pseudoknot structure and a dynamic upstream multibranch loop.