Search

Voltage-sensing phosphatases (VSPs) dephosphorylate phosphoinositide (PIP) signaling lipids in response to membrane depolarization. VSPs possess an S4-containing voltage sensor domain (VSD), resembling that of voltage-gated cation channels, and a lipid phosphatase domain (PD). The mechanism by which voltage turns on enzyme activity is unclear. Structural analysis and modeling suggest several sites of VSD-PD interaction that could couple voltage sensing to catalysis. Voltage clamp fluorometry reveals voltage-driven rearrangements in three sites implicated earlier in enzyme activation-the VSD-PD linker, gating loop and R loop-as well as the N-terminal domain, which has not yet been explored. N-terminus mutations perturb both rearrangements in the other segments and enzyme activity. Our results provide a model for a dynamic assembly by which S4 controls the catalytic site.

Heparan sulfate (HS) is degraded in lysosome by a series of glycosidases. Before the glycosidases can act, the terminal glucosamine of HS must be acetylated by the integral lysosomal membrane enzyme heparan-α-glucosaminide N-acetyltransferase (HGSNAT). Mutations of HGSNAT cause HS accumulation and consequently mucopolysaccharidosis IIIC, a devastating lysosomal storage disease characterized by progressive neurological deterioration and early death where no treatment is available. HGSNAT catalyzes a unique transmembrane acetylation reaction where the acetyl group of cytosolic acetyl-CoA is transported across the lysosomal membrane and attached to HS in one reaction. However, the reaction mechanism remains elusive. Here we report six cryo-EM structures of HGSNAT along the reaction pathway. These structures reveal a dimer arrangement and a unique structural fold, which enables the elucidation of the reaction mechanism. We find that a central pore within each monomer traverses the membrane and controls access of cytosolic acetyl-CoA to the active site at its luminal mouth where glucosamine binds. A histidine-aspartic acid catalytic dyad catalyzes the transfer reaction via a ternary complex mechanism. Furthermore, the structures allow the mapping of disease-causing variants and reveal their potential impact on the function, thus creating a framework to guide structure-based drug discovery efforts.

Crystal structures of human long-chain acyl-CoA dehydrogenase (LCAD) and the catalytically inactive Glu291Gln mutant, have been determined. These structures suggest that LCAD harbors functions beyond its historically defined role in mitochondrial β-oxidation of long and medium-chain fatty acids. LCAD is a homotetramer containing one FAD per 43 kDa subunit with Glu291 as the catalytic base. The substrate binding cavity of LCAD reveals key differences which makes it specific for longer and branched chain substrates. The presence of Pro132 near the start of the E helix leads to helix unwinding that, together with adjacent smaller residues, permits binding of bulky substrates such as 3α, 7α, l2α-trihydroxy-5β-cholestan-26-oyl-CoA. This structural element is also utilized by ACAD11, a eucaryotic ACAD of unknown function, as well as bacterial ACADs known to metabolize sterol substrates. Sequence comparison suggests that ACAD10, another ACAD of unknown function, may also share this substrate specificity. These results suggest that LCAD, ACAD10, ACAD11 constitute a distinct class of eucaryotic acyl CoA dehydrogenases.

Insertion sequence (IS) elements are the simplest autonomous transposable elements found in prokaryotic genomes1. We recently discovered that IS110 family elements encode a recombinase and a non-coding bridge RNA (bRNA) that confers modular specificity for target DNA and donor DNA through two programmable loops2. Here we report the cryo-electron microscopy structures of the IS110 recombinase in complex with its bRNA, target DNA and donor DNA in three different stages of the recombination reaction cycle. The IS110 synaptic complex comprises two recombinase dimers, one of which houses the target-binding loop of the bRNA and binds to target DNA, whereas the other coordinates the bRNA donor-binding loop and donor DNA. We uncovered the formation of a composite RuvC-Tnp active site that spans the two dimers, positioning the catalytic serine residues adjacent to the recombination sites in both target and donor DNA. A comparison of the three structures revealed that (1) the top strands of target and donor DNA are cleaved at the composite active sites to form covalent 5-phosphoserine intermediates, (2) the cleaved DNA strands are exchanged and religated to create a Holliday junction intermediate, and (3) this intermediate is subsequently resolved by cleavage of the bottom strands. Overall, this study reveals the mechanism by which a bispecific RNA confers target and donor DNA specificity to IS110 recombinases for programmable DNA recombination.

The retinal light response in animals originates from the photoisomerization of an opsin-coupled 11-cis-retinaldehyde chromophore. This visual chromophore is enzymatically produced through the action of carotenoid cleavage dioxygenases. Vertebrates require two carotenoid cleavage dioxygenases, β-carotene oxygenase 1 and retinal pigment epithelium 65 (RPE65), to form 11-cis-retinaldehyde from carotenoid substrates, whereas invertebrates such as insects use a single enzyme known as Neither Inactivation Nor Afterpotential B (NinaB). RPE65 and NinaB couple trans-cis isomerization with hydrolysis and oxygenation, respectively, but the mechanistic relationship of their isomerase activities remains unknown. Here we report the structure of NinaB, revealing details of its active site architecture and mode of membrane binding. Structure-guided mutagenesis studies identify a residue cluster deep within the NinaB substrate-binding cleft that controls its isomerization activity. Our data demonstrate that isomerization activity is mediated by distinct active site regions in NinaB and RPE65-an evolutionary convergence that deepens our understanding of visual system diversity.

Unlocking the potential of protein arginine deiminase 4 (PAD4) as a drug target for rheumatoid arthritis requires a deeper understanding of its regulation. In this study, we use unbiased antibody selections to identify functional antibodies capable of either activating or inhibiting PAD4 activity. Through cryogenic-electron microscopy, we characterized the structures of these antibodies in complex with PAD4 and revealed insights into their mechanisms of action. Rather than steric occlusion of the substrate-binding catalytic pocket, the antibodies modulate PAD4 activity through interactions with allosteric binding sites adjacent to the catalytic pocket. These binding events lead to either alteration of the active site conformation or the enzyme oligomeric state, resulting in modulation of PAD4 activity. Our study uses antibody engineering to reveal new mechanisms for enzyme regulation and highlights the potential of using PAD4 agonist and antagonist antibodies for studying PAD4-dependency in disease models and future therapeutic development.

CRISPR-Cas are adaptive immune systems in bacteria and archaea that utilize CRISPR RNA-guided surveillance complexes to target complementary RNA or DNA for destruction1-5. Target RNA cleavage at regular intervals is characteristic of type III effector complexes6-8. Here, we determine the structures of the Synechocystis type III-Dv complex, an apparent evolutionary intermediate from multi-protein to single-protein type III effectors9,10, in pre- and post-cleavage states. The structures show how multi-subunit fusion proteins in the effector are tethered together in an unusual arrangement to assemble into an active and programmable RNA endonuclease and how the effector utilizes a distinct mechanism for target RNA seeding from other type III effectors. Using structural, biochemical, and quantum/classical molecular dynamics simulation, we study the structure and dynamics of the three catalytic sites, where a 2-OH of the ribose on the target RNA acts as a nucleophile for in line self-cleavage of the upstream scissile phosphate. Strikingly, the arrangement at the catalytic residues of most type III complexes resembles the active site of ribozymes, including the hammerhead, pistol, and Varkud satellite ribozymes. Our work provides detailed molecular insight into the mechanisms of RNA targeting and cleavage by an important intermediate in the evolution of type III effector complexes.

The innate immune system relies on molecular sensors to detect distinctive molecular patterns, including viral double-stranded RNA (dsRNA), which triggers responses resulting in apoptosis and immune infiltration. Adenosine Deaminases Acting on RNA (ADARs) catalyze the deamination of adenosine (A) to inosine (I), serving as a mechanism to distinguish self from non-self RNA and prevent aberrant immune activation. Loss-of-function mutations in the ADAR1 gene are one cause of Aicardi Goutières Syndrome (AGS), a severe autoimmune disorder in children. Although seven out of the eight AGS-associated mutations in ADAR1 occur within the catalytic domain of the ADAR1 protein, their specific effects on the catalysis of adenosine deamination remain poorly understood. In this study, we carried out a biochemical investigation of four AGS-causing mutations (G1007R, R892H, K999N, and Y1112F) in ADAR1 p110 and truncated variants. These studies included adenosine deamination rate measurements with two different RNA substrates derived from human transcripts known to be edited by ADAR1 p110 (glioma-associated oncogene homologue 1 (hGli1), 5-hydroxytryptamine receptor 2C (5-HT2cR)). Our results indicate that AGS-associated mutations at two amino acid positions directly involved in stabilizing the base-flipped conformation of the ADAR-RNA complex (G1007R and R892H) had the most detrimental impact on catalysis. The K999N mutation, positioned near the RNA binding interface, altered catalysis contextually. Finally, the Y1112F mutation had small effects in each of the assays described here. These findings shed light on the differential effects of disease-associated mutations on adenosine deamination by ADAR1, thereby advancing our structural and functional understanding of ADAR1-mediated RNA editing.

The long interspersed element-1 (LINE-1, hereafter L1) retrotransposon has generated nearly one-third of the human genome and serves as an active source of genetic diversity and human disease1. L1 spreads through a mechanism termed target-primed reverse transcription, in which the encoded enzyme (ORF2p) nicks the target DNA to prime reverse transcription of its own or non-self RNAs2. Here we purified full-length L1 ORF2p and biochemically reconstituted robust target-primed reverse transcription with template RNA and target-site DNA. We report cryo-electron microscopy structures of the complete human L1 ORF2p bound to structured template RNAs and initiating cDNA synthesis. The template polyadenosine tract is recognized in a sequence-specific manner by five distinct domains. Among them, an RNA-binding domain bends the template backbone to allow engagement of an RNA hairpin stem with the L1 ORF2p C-terminal segment. Moreover, structure and biochemical reconstitutions demonstrate an unexpected target-site requirement: L1 ORF2p relies on upstream single-stranded DNA to position the adjacent duplex in the endonuclease active site for nicking of the longer DNA strand, with a single nick generating a staggered DNA break. Our research provides insights into the mechanism of ongoing transposition in the human genome and informs the engineering of retrotransposon proteins for gene therapy.

Complex II (CII) activity controls phenomena that require crosstalk between metabolism and signaling, including neurodegeneration, cancer metabolism, immune activation, and ischemia-reperfusion injury. CII activity can be regulated at the level of assembly, a process that leverages metastable assembly intermediates. The nature of these intermediates and how CII subunits transfer between metastable complexes remains unclear. In this work, we identify metastable species containing the SDHA subunit and its assembly factors, and we assign a preferred temporal sequence of appearance of these species during CII assembly. Structures of two species show that the assembly factors undergo disordered-to-ordered transitions without the appearance of significant secondary structure. The findings identify that intrinsically disordered regions are critical in regulating CII assembly, an observation that has implications for the control of assembly in other biomolecular complexes.

Protein functions are dynamically regulated by allostery, which enables conformational communication even between faraway residues, and expresses itself in many forms, akin to different languages: allosteric control pathways predominating in an unperturbed protein are often unintuitively reshaped whenever biochemical perturbations arise (e.g., mutations). To accurately model allostery, unbiased molecular dynamics (MD) simulations require integration with a reliable method able to, e.g., detect incipient allosteric changes or likely perturbation pathways; this is because allostery can operate at longer time scales than those accessible by plain MD. Such methods are typically applied singularly, but we here argue their joint application─as a multilingual approach─could work significantly better. We successfully prove this through unbiased MD simulations (∼100 μs) of the widely studied, allosterically active oncotarget K-Ras4B, solvated and embedded in a phospholipid membrane, from which we decrypt allostery using four showcase languages: Distance Fluctuation analysis and the Shortest Path Map capture allosteric hotspots at equilibrium; Anisotropic Thermal Diffusion and Dynamical Non-Equilibrium MD simulations assess perturbations upon, respectively, either superheating or hydrolyzing the GTP that oncogenically activates K-Ras4B. Chosen languages work synergistically, providing an articulate, mutually coherent, experimentally consistent picture of K-Ras4B allostery, whereby distinct traits emerge at equilibrium and upon GTP cleavage. At equilibrium, combined evidence confirms prominent allosteric communication from the membrane-embedded hypervariable region, through a hub comprising helix α5 and sheet β5, and up to the active site, encompassing allosteric switches I and II (marginally), and two proposed pockets. Upon GTP cleavage, allosteric perturbations mostly accumulate on the switches and documented interfaces.

Protein kinase A (PKA) neuronal function is controlled by the interaction of a regulatory (R) subunit dimer with two catalytic subunits. Recently, the L50R variant in the gene encoding the RIβ subunit was identified in individuals with a novel neurodegenerative disease. However, the mechanisms driving the disease phenotype remained unknown. In this study, we generated a mouse model carrying the RIβ-L50R mutation to replicate the human disease phenotype and study its progression with age. We examined post-mortem brains of affected individuals as well as live cell cultures. Employing biochemical assays, immunohistochemistry and behavioural assessments, we investigated the impact of the mutation on PKA complex assembly, protein aggregation and neuronal degeneration. We reveal that RIβ is an aggregation-prone protein that progressively accumulates in wildtype and Alzheimer's mouse models with age, while aggregation is accelerated in the RIβ-L50R mouse model. We define RIβ-L50R as a causal mutation driving an age-dependent behavioural and disease phenotype in human and mouse models. Mechanistically, this mutation disrupts RIβ dimerization, leading to aggregation of its monomers. Intriguingly, interaction with the catalytic subunit protects the RIβ-L50R from self-aggregating, in a dose-dependent manner. Furthermore, cAMP signaling induces RIβ-L50R aggregation. The pathophysiological mechanism elucidated here for a newly recognized neurodegenerative disease, in which protein aggregation is the result of disrupted homodimerization, sheds light on a remarkably under-appreciated but potentially common mechanism across several neurodegenerative diseases. [ABSTRACT FROM AUTHOR]

Esophageal cancer is one of the most common malignant tumors, and the 5-year overall survival rate is only 20%. Esophageal squamous cell carcinoma (ESCC) is the primary histological type of esophageal carcinoma in China. Protein phosphatase 1 regulatory subunit 18 (PPP1r18) is one of the actin-regulatory proteins and is able to bind to protein phosphatase 1 catalytic subunit alpha (PPP1CA). Yet, little is known about the role of PPP1r18 in ESCC. This study aimed to elucidate the biological functions of PPP1r18 in the ESCC progression. Clinical samples first confirmed that PPP1r18 expression was upregulated in ESCC, and PPP1r18 was correlated with tumor invasion depth, lymph node metastasis, distant metastasis and reduced overall survival. We then observed that PPP1r18 overexpression enhanced cell proliferation in vitro and in vivo. Mechanistically, PPP1r18 regulated tumor progression of ESCC through activating the calcineurin-mediated ERK pathway, rather than binding to PPP1CA. Collectively, our results suggest that PPP1r18 promotes ESCC progression by regulating the calcineurin-mediated ERK pathway. PPP1r18 might be a potential target for the diagnosis and treatment of ESCC. [ABSTRACT FROM AUTHOR]

Aerobic ribonucleotide reductases (RNRs) initiate synthesis of DNA building blocks by generating a free radical within the R2 subunit; the radical is subsequently shuttled to the catalytic R1 subunit through proton-coupled electron transfer (PCET). We present a high-resolution room temperature structure of the class Ie R2 protein radical captured by x-ray free electron laser serial femtosecond crystallography. The structure reveals conformational reorganization to shield the radical and connect it to the translocation path, with structural changes propagating to the surface where the protein interacts with the catalytic R1 subunit. Restructuring of the hydrogen bond network, including a notably short O-O interaction of 2.41 angstroms, likely tunes and gates the radical during PCET. These structural results help explain radical handling and mobilization in RNR and have general implications for radical transfer in proteins.

Human cytochrome P450 3A4 (CYP3A4) is a major drug-metabolizing enzyme that shows extreme substrate promiscuity. Moreover, its large and malleable active site can simultaneously accommodate several substrate molecules of the same or different nature, which may lead to cooperative binding and allosteric behavior. Due to difficulty of crystallization of CYP3A4-substrate complexes, it remains unknown how multiple substrates can arrange in the active site. We determined crystal structures of CYP3A4 bound to three and six molecules of caffeine, a psychoactive alkaloid serving as a substrate and modulator of CYP3A4. In the ternary complex, one caffeine binds to the active site suitably for C8-hydroxylation, most preferable for CYP3A4. In the senary complex, three caffeine molecules stack parallel to the heme with the proximal ligand poised for 3-N-demethylation. However, the caffeine stack forms extensive hydrophobic interactions that could preclude product dissociation and multiple turnovers. In both complexes, caffeine is also bound in the substrate channel and on the outer surface known as a peripheral site. At all sites, aromatic stacking with the caffeine ring(s) is likely a dominant interaction, while direct and water-mediated polar contacts provide additional stabilization for the substrate-bound complexes. Protein-ligand interactions via the active site R212, intrachannel T224, and peripheral F219 were experimentally confirmed, and the latter two residues were identified as important for caffeine association. Collectively, the structural, spectral, and mutagenesis data provide valuable insights on the ligand binding mechanism and help better understand how purine-based pharmaceuticals and other aromatic compounds could interact with CYP3A4 and mediate drug-drug interactions.

PIK3R1 (also known as p85α) is a regulatory subunit of phosphoinositide 3-kinases (PI3Ks). PI3K, a heterodimer of a regulatory subunit and a catalytic subunit, phosphorylates phosphatidylinositol into secondary signaling molecules involved in regulating metabolic homeostasis. PI3K converts phosphatidylinositol 4,5-bisphosphate (PIP2) to phosphatidylinositol 3,4,5-triphosphate (PIP3), which recruits protein kinase AKT to the inner leaflet of the cell membrane to be activated and to participate in various metabolic functions. PIK3R1 stabilizes and inhibits p110 catalytic activity and serves as an adaptor to interact with insulin receptor substrate (IRS) proteins and growth factor receptors. Thus, mutations in PIK3R1 or altered expression of PIK3R1 could modulate the activity of PI3K and result in significant metabolic outcomes. Interestingly, recent studies also found PI3K-independent functions of PIK3R1. Overall, in this article, we will provide an updated review of the metabolic functions of PIK3R1 that includes studies of PIK3R1 in various metabolic tissues using animal models, the mechanisms modulating PIK3R1 activity, and studies on the mutations of human PIK3R1 gene.

The Protein Phosphatase type 1 catalytic subunit (PP1c) (PF3D7_1414400) operates in combination with various regulatory proteins to specifically direct and control its phosphatase activity. However, there is little information about this phosphatase and its regulators in the human malaria parasite, Plasmodium falciparum. To address this knowledge gap, we conducted a comprehensive investigation into the structural and functional characteristics of a conserved Plasmodium-specific regulator called Gametocyte EXported Protein 15, GEXP15 (PF3D7_1031600). Through in silico analysis, we identified three significant regions of interest in GEXP15: an N-terminal region housing a PP1-interacting RVxF motif, a conserved domain whose function is unknown, and a GYF-like domain that potentially facilitates specific protein-protein interactions. To further elucidate the role of GEXP15, we conducted in vitro interaction studies that demonstrated a direct interaction between GEXP15 and PP1 via the RVxF-binding motif. This interaction was found to enhance the phosphatase activity of PP1. Additionally, utilizing a transgenic GEXP15-tagged line and live microscopy, we observed high expression of GEXP15 in late asexual stages of the parasite, with localization predominantly in the nucleus. Immunoprecipitation assays followed by mass spectrometry analyses revealed the interaction of GEXP15 with ribosomal- and RNA-binding proteins. Furthermore, through pull-down analyses of recombinant functional domains of His-tagged GEXP15, we confirmed its binding to the ribosomal complex via the GYF domain. Collectively, our study sheds light on the PfGEXP15-PP1-ribosome interaction, which plays a crucial role in protein translation. These findings suggest that PfGEXP15 could serve as a potential target for the development of malaria drugs.

An increased understanding of how the acceptor site in Gcn5-related N-acetyltransferase (GNAT) enzymes recognizes various substrates provides important clues for GNAT functional annotation and their use as chemical tools. In this study, we explored how the PA3944 enzyme from Pseudomonas aeruginosa recognizes three different acceptor substrates, including aspartame, NANMO, and polymyxin B, and identified acceptor residues that are critical for substrate specificity. To achieve this, we performed a series of molecular docking simulations and tested methods to identify acceptor substrate binding modes that are catalytically relevant. We found that traditional selection of best docking poses by lowest S scores did not reveal acceptor substrate binding modes that were generally close enough to the donor for productive acetylation. Instead, sorting poses based on distance between the acceptor amine nitrogen atom and donor carbonyl carbon atom placed these acceptor substrates near residues that contribute to substrate specificity and catalysis. To assess whether these residues are indeed contributors to substrate specificity, we mutated seven amino acid residues to alanine and determined their kinetic parameters. We identified several residues that improved the apparent affinity and catalytic efficiency of PA3944, especially for NANMO and/or polymyxin B. Additionally, one mutant (R106A) exhibited substrate inhibition toward NANMO, and we propose scenarios for the cause of this inhibition based on additional substrate docking studies with R106A. Ultimately, we propose that this residue is a key gatekeeper between the acceptor and donor sites by restricting and orienting the acceptor substrate within the acceptor site.

SARS-CoV-2 RNA dependent RNA polymerase (RdRp) serves as a highly promising antiviral drug target such as for a Remdesivir nucleotide analogue (RDV-TP or RTP). In this work, we mainly used alchemical all-atom simulations to characterize relative binding free energetics between the nucleotide analogue RTP and natural cognate substrate ATP upon initial binding and pre-catalytic insertion into the active site of SARS-CoV-2 RdRp. Natural non-cognate substrate dATP and mismatched GTP were also examined for computation control. We first identified significant differences in dynamical responses between nucleotide initial binding and subsequent insertion configurations to the open and closed active sites of the RdRp, respectively, though the RdRp protein conformational changes between the active site's open and closed states are subtle. Our alchemical simulations indicated that upon initial binding (active site open), RTP and ATP show similar binding free energies to the active sites while in the insertion state (active site closed), ATP is more stabilized (∼-2.4 kcal mol-1) than RTP in free energetics. Additional analyses show, however, that RTP is more stabilized in binding energetics than ATP, in both the insertion and initial binding states, with RTP more stabilized due to the electrostatic energy in the insertion state and due to vdW energy in the initial binding state. Hence, it appears that natural cognate ATP still excels at association stability with the RdRp active site due to that ATP maintains sufficient flexibilities e.g., in base pairing with the template, which exemplifies an entropic contribution to the cognate substrate stabilization. These findings highlight the importance of substrate flexibilities in addition to energetic stabilization in antiviral nucleotide analogue design.

The lipolytic processing of triglyceride-rich lipoproteins (TRLs) by lipoprotein lipase (LPL) is crucial for the delivery of dietary lipids to the heart, skeletal muscle, and adipose tissue. The processing of TRLs by LPL is regulated in a tissue-specific manner by a complex interplay between activators and inhibitors. Angiopoietin-like protein 4 (ANGPTL4) inhibits LPL by reducing its thermal stability and catalyzing the irreversible unfolding of LPLs α/β-hydrolase domain. We previously mapped the ANGPTL4 binding site on LPL and defined the downstream unfolding events resulting in LPL inactivation. The binding of LPL to glycosylphosphatidylinositol-anchored high-density lipoprotein-binding protein 1 protects against LPL unfolding. The binding site on LPL for an activating cofactor, apolipoprotein C2 (APOC2), and the mechanisms by which APOC2 activates LPL have been unclear and controversial. Using hydrogen-deuterium exchange/mass spectrometry, we now show that APOC2s C-terminal α-helix binds to regions of LPL surrounding the catalytic pocket. Remarkably, APOC2s binding site on LPL overlaps with that for ANGPTL4, but their effects on LPL conformation are distinct. In contrast to ANGPTL4, APOC2 increases the thermal stability of LPL and protects it from unfolding. Also, the regions of LPL that anchor the lid are stabilized by APOC2 but destabilized by ANGPTL4, providing a plausible explanation for why APOC2 is an activator of LPL, while ANGPTL4 is an inhibitor. Our studies provide fresh insights into the molecular mechanisms by which APOC2 binds and stabilizes LPL-and properties that we suspect are relevant to the conformational gating of LPLs active site.

Modular polyketide synthases (PKSs) are polymerases that employ α-carboxyacyl-CoAs as extender substrates. This enzyme family contains several catalytic modules, where each module is responsible for a single round of polyketide chain extension. Although PKS modules typically use malonyl-CoA or methylmalonyl-CoA for chain elongation, many other malonyl-CoA analogues are used to diversify polyketide structures in nature. Previously, we developed a method to alter an extension substrate of a given module by exchanging an acyltransferase (AT) domain while maintaining protein folding. Here, we report in vitro polyketide biosynthesis by 13 PKSs (the wild-type PKS and 12 AT-exchanged PKSs with unusual ATs) and 14 extender substrates. Our ∼200 in vitro reactions resulted in 13 structurally different polyketides, including several polyketides that have not been reported. In some cases, AT-exchanged PKSs produced target polyketides by >100-fold compared to the wild-type PKS. These data also indicate that most unusual AT domains do not incorporate malonyl-CoA and methylmalonyl-CoA but incorporate various rare extender substrates that are equal to in size or slightly larger than natural substrates. We developed a computational workflow to predict the approximate AT substrate range based on active site volumes to support the selection of ATs. These results greatly enhance our understanding of rare AT domains and demonstrate the benefit of using the proposed PKS engineering strategy to produce novel chemicals in vitro.

Transition metal chemistry is essential to life, where metal binding to DNA, RNA, and proteins underpins all facets of the central dogma of biology. In this context, metals in proteins are typically studied as static active site cofactors. However, the emergence of transition metal signaling, where mobile metal pools can transiently bind to biological targets beyond active sites, is expanding this conventional view of bioinorganic chemistry. This Minireview focuses on the concept of metalloallostery, using copper as a canonical example of how metals can regulate protein function by binding to remote allosteric sites (e.g., exosites). We summarize advances in and prospects for the field, including imaging dynamic transition metal signaling pools, allosteric inhibition or activation of protein targets by metal binding, and metal-dependent signaling pathways that underlie nutrient vulnerabilities in diseases spanning obesity, fatty liver disease, cancer, and neurodegeneration.

Macrophage migration inhibitory factor (MIF) and D-dopachrome tautomerase (D-DT) are two pleotropic cytokines, which are coexpressed in various cell types to activate the cell surface receptor CD74. Via the MIF/CD74 and D-DT/CD74 axes, the two proteins exhibit either beneficial or deleterious effect on human diseases. In this study, we report the identification of 2,5-pyridinedicarboxylic acid (a.k.a. 1) that effectively blocks the D-DT-induced activation of CD74 and demonstrates an impressive 79-fold selectivity for D-DT over MIF. Crystallographic characterization of D-DT-1 elucidates the binding features of 1 and reveals previously unrecognized differences between the MIF and D-DT active sites that explain the ligand's functional selectivity. The commercial availability, low cost, and high selectivity make 1 the ideal tool for studying the pathophysiological functionality of D-DT in disease models. At the same time, our comprehensive biochemical, computational, and crystallographic analyses serve as a guide for generating highly potent and selective D-DT inhibitors.

The transcription factor TFEB is a master regulator of lysosomal biogenesis and autophagy1. The phosphorylation of TFEB by the mechanistic target of rapamycin complex 1 (mTORC1)2-5 is unique in its mTORC1 substrate recruitment mechanism, which is strictly dependent on the amino acid-mediated activation of the RagC GTPase activating protein FLCN6,7. TFEB lacks the TOR signalling motif responsible for the recruitment of other mTORC1 substrates. We used cryogenic-electron microscopy to determine the structure of TFEB as presented to mTORC1 for phosphorylation, which we refer to as the 'megacomplex'. Two full Rag-Ragulator complexes present each molecule of TFEB to the mTOR active site. One Rag-Ragulator complex is bound to Raptor in the canonical mode seen previously in the absence of TFEB. A second Rag-Ragulator complex (non-canonical) docks onto the first through a RagC GDP-dependent contact with the second Ragulator complex. The non-canonical Rag dimer binds the first helix of TFEB with a RagCGDP-dependent aspartate clamp in the cleft between the Rag G domains. In cellulo mutation of the clamp drives TFEB constitutively into the nucleus while having no effect on mTORC1 localization. The remainder of the 108-amino acid TFEB docking domain winds around Raptor and then back to RagA. The double use of RagC GDP contacts in both Rag dimers explains the strong dependence of TFEB phosphorylation on FLCN and the RagC GDP state.

Heterotrimeric guanine nucleotide-binding proteins (G proteins) that function as molecular switches for cellular growth and metabolism are activated by GTP and inactivated by GTP hydrolysis. In uveal melanoma, a conserved glutamine residue critical for GTP hydrolysis in the G protein α subunit is often mutated in Gαq or Gα11 to either leucine or proline. In contrast, other glutamine mutations or mutations in other Gα subtypes are rare. To uncover the mechanism of the genetic selection and the functional role of this glutamine residue, we analyzed all possible substitutions of this residue in multiple Gα isoforms. Through cell-based measurements of activity, we showed that some mutants were further activated and inactivated by G protein-coupled receptors. Through biochemical, molecular dynamics, and nuclear magnetic resonance-based structural studies, we showed that the Gα mutants were functionally distinct and conformationally diverse, despite their shared inability to hydrolyze GTP. Thus, the catalytic glutamine residue contributes to functions beyond GTP hydrolysis, and these functions include subtype-specific, allosteric modulation of receptor-mediated subunit dissociation. We conclude that G proteins do not function as simple on-off switches. Rather, signaling emerges from an ensemble of active states, a subset of which are favored in disease and may be uniquely responsive to receptor-directed ligands.

Enzymes host active sites inside protein macromolecules, which have diverse, often incredibly complex, and atom-expensive structures. It is an outstanding question what the role of these expensive scaffolds might be in enzymatic catalysis. Answering this question is essential to both enzymology and the design of artificial enzymes with proficiencies that will match those of the best natural enzymes. Protein rigidifying the active site, contrasted with the dynamics and vibrational motion promoting the reaction, as well as long-range electrostatics (also known as electrostatic preorganization) were all proposed as central contributions of the scaffold to the catalysis. Here, we show that all these effects inevitably produce changes in the quantum mechanical electron density in the active site, which in turn defines the reactivity. The phenomena are therefore fundamentally inseparable. The geometry of the electron density-a scalar field characterized by a number of mathematical features such as critical points-is a rigorous and convenient descriptor of enzymatic catalysis and a reporter on the role of the protein. We show how this geometry can be analyzed, linked to the reaction barriers, and report in particular on intramolecular electric fields in enzymes. We illustrate these tools on the studies of electrostatic preorganization in several representative enzyme classes, both natural and artificial. We highlight the forward-looking aspects of the approach.

De novo enzyme design has sought to introduce active sites and substrate-binding pockets that are predicted to catalyse a reaction of interest into geometrically compatible native scaffolds1,2, but has been limited by a lack of suitable protein structures and the complexity of native protein sequence-structure relationships. Here we describe a deep-learning-based 'family-wide hallucination' approach that generates large numbers of idealized protein structures containing diverse pocket shapes and designed sequences that encode them. We use these scaffolds to design artificial luciferases that selectively catalyse the oxidative chemiluminescence of the synthetic luciferin substrates diphenylterazine3 and 2-deoxycoelenterazine. The designed active sites position an arginine guanidinium group adjacent to an anion that develops during the reaction in a binding pocket with high shape complementarity. For both luciferin substrates, we obtain designed luciferases with high selectivity; the most active of these is a small (13.9 kDa) and thermostable (with a melting temperature higher than 95 °C) enzyme that has a catalytic efficiency on diphenylterazine (kcat/Km = 106 M-1 s-1) comparable to that of native luciferases, but a much higher substrate specificity. The creation of highly active and specific biocatalysts from scratch with broad applications in biomedicine is a key milestone for computational enzyme design, and our approach should enable generation of a wide range of luciferases and other enzymes.

Objective To identify the TMPRSS6 gene variants in Mexican patients with iron treatment refractoriness, to describe hematological and iron profile parameters, and to use bioinformatic prediction and protein modeling tools to assess a possible biological impact for the detected missense variants. Methods Nineteen patients referred with iron treatment refractoriness were studied. Peripheral blood was collected to determine hematic cytometry, iron profile, hemoglobin electrophoresis, and quantification. Molecular screening was carried out for exons 15 through 18 of the TMPRSS6 gene by Sanger sequencing and for frequent thalassemia variants by amplification-refractory mutation system-polymerase chain reaction (PCR) and gap-PCR. The biological impact of the detected missense variants was assessed using bioinformatic prediction and protein modeling tools. Results We found 5 genetic variants in the matriptase-2 catalytic domain: 1 at intron-15/exon-16 junction (rs60484081) and 4 exonic, 3 missense (rs377054987, p.Gly626Asp; rs1384127820, p.Ser672Thr; rs855791, p.Val727Ala) and 1 synonymous (rs2235321, p.Tyr730=), with frequencies ranging from 0.18 to 0.53. No significant differences were observed in the hematological parameters or iron profile, considering type and number of variants. Bioinformatic predictions suggested a possible biological impact only for rs377054987. Conclusions The TMPRSS6 variants observed in Mexican patients with oral iron treatment refractoriness have high frequencies; nevertheless, their relationship with hematological and iron profile parameters needs further research. The possible biological impact for rs377054987 is due to size and amino acid hydrophobicity changes and hydrogen bond modifications. [ABSTRACT FROM AUTHOR]

Bond bundle analysis is used to investigate enzymatic catalysis in the ketosteroid isomerase (KSI) active site. We identify the unique bonding regions in five KSI systems, including those exposed to applied oriented electric fields and those with amino acid mutations, and calculate the precise redistribution of electron density and other regional properties that accompanies either enhancement or inhibition of KSI catalytic activity. We find that catalytic enhancement results from promoting both inter- and intra-molecular electron density redistribution, between bond bundles and bond wedges within the KSI-docked substrate molecule, in the forward direction of the catalyzed reaction. Though the redistribution applies to both types of perturbed systems and is thus suggestive of a general catalytic role, we observe that bond properties (e.g., volume vs energy vs electron count) can respond independently and disproportionately depending on the type of perturbation. We conclude that the resulting catalytic enhancement/inhibition proceeds via different mechanisms, where some bond properties are utilized more by one type of perturbation than the other. Additionally, we find that the correlations between bond wedge properties and catalyzed reaction barrier energies are additive to predict those of bond bundles and atomic basins, providing a rigorous grounding for connecting changes in local charge density to resulting shifts in reaction barrier energy.

Ancestral sequence reconstruction (ASR) infers predicted ancestral states for sites within sequences and can constrain the functions and properties of ancestors of extant protein families. Here, we compare the likely sequences of inferred nitrogenase ancestors to extant nitrogenase sequence diversity. We show that the most-likely combinations of ancestral states for key substrate channel residues are not represented in extant sequence space, and rarely found within a more broadly defined physiochemical space-supporting that the earliest ancestors of extant nitrogenases likely had alternative substrate channel composition. These differences may indicate differing environmental selection pressures acting on nitrogenase substrate specificity in ancient environments. These results highlight ASR's potential as an in silico tool for developing hypotheses about ancestral enzyme functions, as well as improving hypothesis testing through more targeted in vitro and in vivo experiments.

β-lactamases inactivate β-lactam antibiotics leading to drug resistance. Consequently, inhibitors of β-lactamases can combat this resistance, and the β-lactamase inhibitory protein (BLIP) is a naturally occurring inhibitor. The widespread CTX-M-14 and CTX-M-15 β-lactamases have an 83% sequence identity. In this study, we show that BLIP weakly inhibits CTX-M-14 but potently inhibits CTX-M-15. The structure of the BLIP/CTX-M-15 complex reveals that binding is associated with a conformational change of an active site loop of β-lactamase. Surprisingly, the loop structure in the complex is similar to that in a drug-resistant variant (N106S) of CTX-M-14. We hypothesized that the pre-established favorable loop conformation of the N106S mutant would facilitate binding. The N106S substitution results in a ~100- and 10-fold increase in BLIP inhibition potency for CTX-M-14 and CTX-M-15, respectively. Thus, this indicates that an active site loop in β-lactamase toggles between conformations that control antibiotic hydrolysis and inhibitor susceptibility. These findings highlight the role of accessible active site conformations in controlling enzyme activity and inhibitor susceptibility as well as the influence of mutations in selectively stabilizing discrete conformations.

Non-ribosomal peptides play a critical role in the clinic as therapeutic agents. To access more chemically diverse therapeutics, non-ribosomal peptide synthetases (NRPSs) have been targeted for engineering through combinatorial biosynthesis; however, this has been met with limited success in part due to the lack of proper protein-protein interactions between non-cognate proteins. Herein, we report our use of chemical biology to enable X-ray crystallography, molecular dynamics (MD) simulations, and biochemical studies to elucidate binding specificities between peptidyl carrier proteins (PCPs) and adenylation (A) domains. Specifically, we determined X-ray crystal structures of a type II PCP crosslinked to its cognate A domain, PigG and PigI, and of PigG crosslinked to a non-cognate PigI homologue, PltF. The crosslinked PCP-A domain structures possess large protein-protein interfaces that predominantly feature hydrophobic interactions, with specific electrostatic interactions that orient the substrate for active site delivery. MD simulations of the PCP-A domain complexes and unbound PCP structures provide a dynamical evaluation of the transient interactions formed at PCP-A domain interfaces, which confirm the previously hypothesized role of a PCP loop as a crucial recognition element. Finally, we demonstrate that the interfacial interactions at the PCP loop 1 region can be modified to control PCP binding specificity through gain-of-function mutations. This work suggests that loop conformational preferences and dynamism account for improved shape complementary in the PCP-A domain interactions. Ultimately, these studies show how crystallographic, biochemical, and computational methods can be used to rationally re-engineer NRPSs for non-cognate interactions.

Multiple psychiatric disorders have been associated with abnormalities in both the innate and adaptive immune systems. The role of these abnormalities in pathogenesis, and whether they are driven by psychiatric risk variants, remains unclear. We test for enrichment of GWAS variants associated with multiple psychiatric disorders (cross-disorder or trans-diagnostic risk), or 5 specific disorders (cis-diagnostic risk), in regulatory elements in immune cells. We use three independent epigenetic datasets representing multiple organ systems and immune cell subsets. Trans-diagnostic and cis-diagnostic risk variants (for schizophrenia and depression) are enriched at epigenetically active sites in brain tissues and in lymphoid cells, especially stimulated CD4+ T cells. There is no evidence for enrichment of either trans-risk or cis-risk variants for schizophrenia or depression in myeloid cells. This suggests a possible model where environmental stimuli activate T cells to unmask the effects of psychiatric risk variants, contributing to the pathogenesis of mental health disorders.

Cytochromes P450 (CYPs) are a superfamily of enzymes responsible for biosynthesis and drug metabolism. Monitoring the activity of CYP3A4, the major human drug-metabolizing enzyme, is vital for assessing the metabolism of pharmaceuticals and identifying harmful drug-drug interactions. Existing probes for CYP3A4 are irreversible turn-on substrates that monitor activity at specific time points in end-point assays. To provide a more dynamic approach, we designed, synthesized, and characterized emissive Ir(III) and Ru(II) complexes that allow monitoring of the CYP3A4 active-site occupancy in real time. In the bound state, probe emission is quenched by the active-site heme. Upon displacement from the active site by CYP3A4-specific inhibitors or substrates, these probes show high emission turn-on. Direct probe binding to the CYP3A4 active site was confirmed by X-ray crystallography. The lead Ir(III)-based probe has nanomolar Kd and high selectivity for CYP3A4, efficient cellular uptake, and low toxicity in CYP3A4-overexpressing HepG2 cells.

With the growing acceptance of the contribution of protein conformational ensembles to enzyme catalysis and regulation, research in the field of protein dynamics has shifted toward an understanding of the atomistic properties of protein dynamical networks and the mechanisms and time scales that control such behavior. A full description of an enzymatic reaction coordinate is expected to extend beyond the active site and include site-specific networks that communicate with the protein/water interface. Advances in experimental tools for the spatial resolution of thermal activation pathways are being complemented by biophysical methods for visualizing dynamics in real time. An emerging multidimensional model integrates the impacts of bound substrate/effector on the distribution of protein substates that are in rapid equilibration near room temperature with reaction-specific protein embedded heat transfer conduits.

Reversible phosphorylation is the basis for signal transduction in eukaryotic cells, and this is tightly controlled by the complex interplay of kinases and phosphatases. Many malignancies are characterized by dysregulation of the delicate protein phosphorylation balance. The targeting of protein phosphatases has been gaining attention as their role in cancer development and progression has been elucidated. The protein phosphatase-2A (PP2A) family of phosphatases are the primary cellular serine/threonine phosphatases. PP2A heterotrimers containing the B55α (PR55α) regulatory subunit have been associated with oncogenic signaling, and B55 subunits are found exclusively in forms of PP2A in which the carboxyl terminus of the catalytic subunit (PP2Ac) is methylated. Methylation of PP2Ac is primarily mediated by leucine carboxyl methyltransferase-1 (LCMT-1). Demethylation is controlled by an esterase, PP2A methylesterase (PME-1). We tested two potential LCMT-1 small molecule inhibitors and found that methyl 4-methyl-2-[(2-methylbenzoyl)amino]-5-[[(3-methylphenyl)amino]carbonyl]-3-thiophenecarboxylate (henceforth referred to as Compound 2) significantly inhibited PP2Ac methylation in vitro (p = 0.0024), and in the MDA-MB-231 breast carcinoma (p = 0.0431) and Rosi melanoma (p = 0.0335) cell lines. Compound 2 significantly reduced survival in HEK-293, HS-5, MDA-MB-231 and Rosi cells; and constrained clonogenic colony formation in MCF7, MDA-MB-231 and Rosi cells. The LCMT-1inhibitor induced G0/G1 cell cycle arrest in Rosi cells (p = 0.0193) and induced apoptosis in MDA-MB-231 cells (p < 0.0001). Increased phosphorylation of the receptor-interacting serine/threonine protein kinase-1 (RIPK1) was detected in MDA-MB-231 (p = 0.0273) and Rosi cells (p = 0.0179) in response to treatment with Compound 2. These data add to the body of evidence pointing to LCMT-1 as an oncogenic target.

Improvements in operando spectroscopy have enabled the catalysis community to investigate the dynamic nature of catalysts under operating conditions with increasing detail. Still, the highly dynamic nature of some catalysts, such as fluxional supported subnano clusters, presents a formidable challenge even for the most state-of-the-art techniques. The reason is that such fluxional catalytic interfaces contain a variety of thermally accessible states. Operando spectroscopies used in catalysis generally fall into two categories: ensemble-based techniques, which provide spectra containing the signals of the entire ensemble of states of the catalyst and are not necessarily dominated by the most active species, and localized techniques, which provide atomistic-level information about the dynamics of active sites in a very small area, which might not include the most active species. Combining many different kinds of techniques can provide detailed insight; however, we propose that effective utilization of specific computational techniques and approaches within the fluxionality paradigm can fill the gap and enable atomistic characterization of the most relevant catalytic sites.

The development of new catalyst materials for energy-efficient chemical synthesis is critical as over 80% of industrial processes rely on catalysts, with many of the most energy-intensive processes specifically using heterogeneous catalysis. Catalytic performance is a complex interplay of phenomena involving temperature, pressure, gas composition, surface composition, and structure over multiple length and time scales. In response to this complexity, the integrated approach to heterogeneous dilute alloy catalysis reviewed here brings together materials synthesis, mechanistic surface chemistry, reaction kinetics, in situ and operando characterization, and theoretical calculations in a coordinated effort to develop design principles to predict and improve catalytic selectivity. Dilute alloy catalysts─in which isolated atoms or small ensembles of the minority metal on the host metal lead to enhanced reactivity while retaining selectivity─are particularly promising as selective catalysts. Several dilute alloy materials using Au, Ag, and Cu as the majority host element, including more recently introduced support-free nanoporous metals and oxide-supported nanoparticle "raspberry colloid templated (RCT)" materials, are reviewed for selective oxidation and hydrogenation reactions. Progress in understanding how such dilute alloy catalysts can be used to enhance selectivity of key synthetic reactions is reviewed, including quantitative scaling from model studies to catalytic conditions. The dynamic evolution of catalyst structure and composition studied in surface science and catalytic conditions and their relationship to catalytic function are also discussed, followed by advanced characterization and theoretical modeling that have been developed to determine the distribution of minority metal atoms at or near the surface. The integrated approach demonstrates the success of bridging the divide between fundamental knowledge and design of catalytic processes in complex catalytic systems, which can accelerate the development of new and efficient catalytic processes.

Selective metal coordination is central to the functions of metalloproteins:1,2 each metalloprotein must pair with its cognate metallocofactor to fulfil its biological role3. However, achieving metal selectivity solely through a three-dimensional protein structure is a great challenge, because there is a limited set of metal-coordinating amino acid functionalities and proteins are inherently flexible, which impedes steric selection of metals3,4. Metal-binding affinities of natural proteins are primarily dictated by the electronic properties of metal ions and follow the Irving-Williams series5 (Mn2+  Zn2+) with few exceptions6,7. Accordingly, metalloproteins overwhelmingly bind Cu2+ and Zn2+ in isolation, regardless of the nature of their active sites and their cognate metal ions1,3,8. This led organisms to evolve complex homeostatic machinery and non-equilibrium strategies to achieve correct metal speciation1,3,8-10. Here we report an artificial dimeric protein, (AB)2, that thermodynamically overcomes the Irving-Williams restrictions in vitro and in cells, favouring the binding of lower-Irving-Williams transition metals over Cu2+, the most dominant ion in the Irving-Williams series. Counter to the convention in molecular design of achieving specificity through structural preorganization, (AB)2 was deliberately designed to be flexible. This flexibility enabled (AB)2 to adopt mutually exclusive, metal-dependent conformational states, which led to the discovery of structurally coupled coordination sites that disfavour Cu2+ ions by enforcing an unfavourable coordination geometry. Aside from highlighting flexibility as a valuable element in protein design, our results illustrate design principles for constructing selective metal sequestration agents.

High-fidelity replication of the large RNA genome of coronaviruses (CoVs) is mediated by a 3'-to-5' exoribonuclease (ExoN) in nonstructural protein 14 (nsp14), which excises nucleotides including antiviral drugs misincorporated by the low-fidelity viral RNA-dependent RNA polymerase (RdRp) and has also been implicated in viral RNA recombination and resistance to innate immunity. Here, we determined a 1.6-Å resolution crystal structure of severe acute respiratory syndrome CoV 2 (SARS-CoV-2) ExoN in complex with its essential cofactor, nsp10. The structure shows a highly basic and concave surface flanking the active site, comprising several Lys residues of nsp14 and the N-terminal amino group of nsp10. Modeling suggests that this basic patch binds to the template strand of double-stranded RNA substrates to position the 3' end of the nascent strand in the ExoN active site, which is corroborated by mutational and computational analyses. We also show that the ExoN activity can rescue a stalled RNA primer poisoned with sofosbuvir and allow RdRp to continue its extension in the presence of the chain-terminating drug, biochemically recapitulating proofreading in SARS-CoV-2 replication. Molecular dynamics simulations further show remarkable flexibility of multidomain nsp14 and suggest that nsp10 stabilizes ExoN for substrate RNA binding to support its exonuclease activity. Our high-resolution structure of the SARS-CoV-2 ExoN-nsp10 complex serves as a platform for future development of anticoronaviral drugs or strategies to attenuate the viral virulence.

The biosynthesis of the active site of the [FeFe]-hydrogenases (HydA1), the H-cluster, is of interest because these enzymes are highly efficient catalysts for the oxidation and production of H2. The biosynthesis of the [2Fe]H subcluster of the H-cluster proceeds from simple precursors, which are processed by three maturases: HydG, HydE, and HydF. Previous studies established that HydG produces an Fe(CO)2(CN) adduct of cysteine, which is the substrate for HydE. In this work, we show that by using the synthetic cluster [Fe2(μ-SH)2(CN)2(CO)4]2- active HydA1 can be biosynthesized without maturases HydG and HydE.

This Account summarizes recent findings centered on the role that redox partner binding, allostery, and conformational dynamics plays in cytochrome P450 proton coupled electron transfer. P450s are one of Natures largest enzyme families and it is not uncommon to find a P450 wherever substrate oxidation is required in the formation of essential molecules critical to the life of the organism or in xenobiotic detoxification. P450s can operate on a remarkably large range of substrates from the very small to the very large, yet the overall P450 three-dimensional structure is conserved. Given this conservation of structure, it is generally assumed that the basic catalytic mechanism is conserved. In nearly all P450s, the O2 O-O bond must be cleaved heterolytically enabling one oxygen atom, the distal oxygen, to depart as water and leave behind a heme iron-linked O atom as the powerful oxidant that is used to activate the nearby substrate. For this process to proceed efficiently, externally supplied electrons and protons are required. Two protons must be added to the departing O atom while an electron is transferred from a redox partner that typically contains either a Fe2S2 or FMN redox center. The paradigm P450 used to unravel the details of these mechanisms has been the bacterial CYP101A1 or P450cam. P450cam is specific for its own Fe2S2 redox partner, putidaredoxin or Pdx, and it has long been postulated that Pdx plays an effector/allosteric role by possibly switching P450cam to an active conformation. Crystal structures, spectroscopic data, and direct binding experiments of the P450cam-Pdx complex provide some answers. Pdx shifts the conformation of P450cam to a more open state, a transition that is postulated to trigger the proton relay network required for O2 activation. An essential part of this proton relay network is a highly conserved Asp (sometimes Glu) that is known to be critical for activity in a number of P450s. How this Asp and proton delivery networks are connected to redox partner binding is quite simple. In the closed state, this Asp is tied down by salt bridges, but these salt bridges are ruptured when Pdx binds, leaving the Asp free to serve its role in proton transfer. An alternative hypothesis suggests that a specific proton relay network is not really necessary. In this scenario, the Asp plays a structural role in the open/close transition and merely opening the active site access channel is sufficient to enable solvent protons in for O2 protonation. Experiments designed to test these various hypotheses have revealed some surprises in both P450cam and other bacterial P450s. Molecular dynamics and crystallography show that P450cam can undergo rather significant conformational gymnastics that result in a large restructuring of the active site requiring multiple cis/trans proline isomerizations. It also has been found that X-ray driven substrate hydroxylation is a useful tool for better understanding the role that the essential Asp and surrounding residues play in catalysis. Here we summarize these recent results which provide a much more dynamic picture of P450 catalysis.

m6A methylation is the most abundant and reversible chemical modification on mRNA with approximately one-fourth of eukaryotic mRNAs harboring at least one m6A-modified base. The recruitment of the mRNA m6A methyltransferase writer complex to phase-separated nuclear speckles is likely to be crucial in its regulation; however, control over the activity of the complex remains unclear. Supported by our observation that a core catalytic subunit of the methyltransferase complex, METTL3, is endogenously colocalized within nuclear speckles as well as in noncolocalized puncta, we tracked the components of the complex with a Cry2-METTL3 fusion construct to disentangle key domains and interactions necessary for the phase separation of METTL3. METTL3 is capable of self-interaction and likely provides the multivalency to drive condensation. Condensates in cells necessarily contain myriad components, each with partition coefficients that establish an entropic barrier that can regulate entry into the condensate. In this regard, we found that, in contrast to the constitutive binding of METTL14 to METTL3 in both the diffuse and the dense phase, WTAP only interacts with METTL3 in dense phase and thereby distinguishes METTL3/METTL14 single complexes in the dilute phase from METTL3/METTL14 multicomponent condensates. Finally, control over METTL3/METTL14 condensation is determined by its small molecule cofactor, S-adenosylmethionine (SAM), which regulates conformations of two gate loops, and some cancer-associated mutations near gate loops can impair METTL3 condensation. Therefore, the link between SAM binding and the control of writer complex phase state suggests that the regulation of its phase state is a potentially critical facet of its functional regulation.

NMR-assisted crystallography-the integrated application of solid-state NMR, X-ray crystallography, and first-principles computational chemistry-holds significant promise for mechanistic enzymology: by providing atomic-resolution characterization of stable intermediates in enzyme active sites, including hydrogen atom locations and tautomeric equilibria, NMR crystallography offers insight into both structure and chemical dynamics. Here, this integrated approach is used to characterize the tryptophan synthase α-aminoacrylate intermediate, a defining species for pyridoxal-5'-phosphate-dependent enzymes that catalyze β-elimination and replacement reactions. For this intermediate, NMR-assisted crystallography is able to identify the protonation states of the ionizable sites on the cofactor, substrate, and catalytic side chains as well as the location and orientation of crystallographic waters within the active site. Most notable is the water molecule immediately adjacent to the substrate β-carbon, which serves as a hydrogen bond donor to the ε-amino group of the acid-base catalytic residue βLys87. From this analysis, a detailed three-dimensional picture of structure and reactivity emerges, highlighting the fate of the L-serine hydroxyl leaving group and the reaction pathway back to the preceding transition state. Reaction of the α-aminoacrylate intermediate with benzimidazole, an isostere of the natural substrate indole, shows benzimidazole bound in the active site and poised for, but unable to initiate, the subsequent bond formation step. When modeled into the benzimidazole position, indole is positioned with C3 in contact with the α-aminoacrylate Cβ and aligned for nucleophilic attack. Here, the chemically detailed, three-dimensional structure from NMR-assisted crystallography is key to understanding why benzimidazole does not react, while indole does.

Lipoprotein-associated phospholipase A2 (Lp-PLA2) associates with low- and high-density lipoproteins in human plasma and specifically hydrolyzes circulating oxidized phospholipids involved in oxidative stress. The association of this enzyme with the lipoprotein's phospholipid monolayer to access its substrate is the most crucial first step in its catalytic cycle. The current study demonstrates unequivocally that a significant movement of a major helical peptide region occurs upon membrane binding, resulting in a large conformational change upon Lp-PLA2 binding to a phospholipid surface. This allosteric regulation of an enzyme's activity by a large membrane-like interface inducing a conformational change in the catalytic site defines a unique dimension of allosterism. The mechanism by which this enzyme associates with phospholipid interfaces to select and extract a single phospholipid substrate molecule and carry out catalysis is key to understanding its physiological functioning. A lipidomics platform was employed to determine the precise substrate specificity of human recombinant Lp-PLA2 and mutants. This study uniquely elucidates the association mechanism of this enzyme with membranes and its resulting conformational change as well as the extraction and binding of specific oxidized and short acyl-chain phospholipid substrates. Deuterium exchange mass spectrometry coupled with molecular dynamics simulations was used to define the precise specificity of the subsite for the oxidized fatty acid at the sn-2 position of the phospholipid backbone. Despite the existence of several crystal structures of this enzyme cocrystallized with inhibitors, little was understood about Lp-PLA2's specificity toward oxidized phospholipids.

The 2 major molecular switches in biology, kinases and GTPases, are both contained in the Parkinson disease-related leucine-rich repeat kinase 2 (LRRK2). Using hydrogen-deuterium exchange mass spectrometry (HDX-MS) and molecular dynamics (MD) simulations, we generated a comprehensive dynamic allosteric portrait of the C-terminal domains of LRRK2 (LRRK2RCKW). We identified 2 helices that shield the kinase domain and regulate LRRK2 conformation and function. One helix in COR-B (COR-B Helix) tethers the COR-B domain to the αC helix of the kinase domain and faces its activation loop, while the C-terminal helix (Ct-Helix) extends from the WD40 domain and interacts with both kinase lobes. The Ct-Helix and the N-terminus of the COR-B Helix create a "cap" that regulates the N-lobe of the kinase domain. Our analyses reveal allosteric sites for pharmacological intervention and confirm the kinase domain as the central hub for conformational control.

The botulinum neurotoxin serotype A (BoNT/A) cuts a single peptide bond in SNAP25, an activity used to treat a wide range of diseases. Reengineering the substrate specificity of BoNT/A's protease domain (LC/A) could expand its therapeutic applications; however, LC/A's extended substrate recognition (≈ 60 residues) challenges conventional approaches. We report a directed evolution method for retargeting LC/A and retaining its exquisite specificity. The resultant eight-mutation LC/A (omLC/A) has improved cleavage specificity and catalytic efficiency (1300- and 120-fold, respectively) for SNAP23 versus SNAP25 compared to a previously reported LC/A variant. Importantly, the BoNT/A holotoxin equipped with omLC/A retains its ability to form full-length holotoxin, infiltrate neurons, and cleave SNAP23. The identification of substrate control loops outside BoNT/A's active site could guide the design of improved BoNT proteases and inhibitors.

Assembly-line polyketide synthases, such as the 6-deoxyerythronolide B synthase (DEBS), are large enzyme factories prized for their ability to produce specific and complex polyketide products. By channeling protein-tethered substrates across multiple active sites in a defined linear sequence, these enzymes facilitate programmed small-molecule syntheses that could theoretically be harnessed to access countless polyketide product structures. Using cryogenic electron microscopy to study DEBS module 1, we present a structural model describing this substrate-channeling phenomenon. Our 3.2- to 4.3-angstrom-resolution structures of the intact module reveal key domain-domain interfaces and highlight an unexpected module asymmetry. We also present the structure of a product-bound module that shines light on a recently described “turnstile” mechanism for transient gating of active sites along the assembly line.

Molecular switch proteins whose cycling between states is controlled by opposing regulators1,2 are central to biological signal transduction. As switch proteins function within highly connected interaction networks3, the fundamental question arises of how functional specificity is achieved when different processes share common regulators. Here we show that functional specificity of the small GTPase switch protein Gsp1 in Saccharomyces cerevisiae (the homologue of the human protein RAN)4 is linked to differential sensitivity of biological processes to different kinetics of the Gsp1 (RAN) switch cycle. We make 55 targeted point mutations to individual protein interaction interfaces of Gsp1 (RAN) and show through quantitative genetic5 and physical interaction mapping that Gsp1 (RAN) interface perturbations have widespread cellular consequences. Contrary to expectation, the cellular effects of the interface mutations group by their biophysical effects on kinetic parameters of the GTPase switch cycle and not by the targeted interfaces. Instead, we show that interface mutations allosterically tune the GTPase cycle kinetics. These results suggest a model in which protein partner binding, or post-translational modifications at distal sites, could act as allosteric regulators of GTPase switching. Similar mechanisms may underlie regulation by other GTPases, and other biological switches. Furthermore, our integrative platform to determine the quantitative consequences of molecular perturbations may help to explain the effects of disease mutations that target central molecular switches.