Noncanonical amino acids (ncAAs) for photo- and chemical crosslinking are powerful biochemical tools for studying and manipulating interactions between proteins both in vitro and in intact cells. Since the first crosslinking ncAAs were genetically encoded about 20 years ago, the technology has now ripened beyond the proof-of-principle demonstrations and is contributing to the study of relevant biological questions in the frame of modern integrative approaches. Here, we provide an overview of available photo-activatable ncAAs for photo-crosslinking and electrophilic ncAAs for genetically encoded chemical crosslinking (GECX), with a major focus on the most recent entries such as ncAAs for SuFEx click chemistry and photo-activatable ncAAs for chemical crosslinking. We present recent examples of the application of genetically encoded crosslinkers to capture protein-protein interactions and identify interaction partners in live cells, to investigate molecular mechanisms of protein function, to stabilize protein complexes for structural studies, to derive structural information about protein complexes from the physiological cell environment, up to perspective applications of GECX-ncAAs for the development of covalent drugs.

Post-translational modifications (PTMs) of proteins not only exponentially increase the diversity of proteoforms, but also contribute to dynamically modulating the localization, stability, activity, and interaction of proteins. Understanding the biological consequences and functions of specific PTMs has been challenging for many reasons, including the dynamic nature of many PTMs and the technical limitations to access homogenously modified proteins. The genetic code expansion technology has emerged to provide unique approaches for studying PTMs. Through site-specific incorporation of unnatural amino acids (UAAs) bearing PTMs or their mimics into proteins, genetic code expansion allows the generation of homogenous proteins with site-specific modifications and atomic resolution both in vitro and in vivo. With this technology, various PTMs and mimics have been precisely introduced into proteins. In this review, we summarize the UAAs and approaches that have been recently developed to site-specifically install PTMs and their mimics into proteins for functional studies of PTMs.

BAZ2A promotes migration and invasion in prostate cancer. Two chemical probes, the specific BAZ2-ICR, and the BAZ2/BRD9 cross-reactive GSK2801, interfere with the recognition of acetylated lysines in histones by the bromodomains of BAZ2A and of its BAZ2B paralog. The two chemical probes were tested in prostate cancer cell lines with opposite androgen susceptibility. BAZ2-ICR and GSK2801 showed different cellular efficacies in accordance with their unequal selectivity profiles. Concurrent inhibition of BAZ2 and BRD9 did not reproduce the effects observed with GSK2801, indicating possible off-targets for this chemical probe. On the other hand, the single BAZ2 inhibition by BAZ2-ICR did not phenocopy genetic ablation, demonstrating that bromodomain interference is not sufficient to strongly affect BAZ2A functionality and suggesting a PROTAC-based chemical ablation as an alternative optimization strategy and a possible therapeutic approach. In this context, we also present the crystallographic structures of BAZ2A in complex with the above chemical probes. Binding poses of TP-238 and GSK4027, chemical probes for the bromodomain subfamily I, and two ligands of the CBP/EP300 bromodomains identify additional headgroups for the development of BAZ2A ligands.

Recombinant human growth hormone (rhGH) and GH receptor antagonists (GHAs) are used clinically to treat a range of disorders associated with GH deficiency or hypersecretion, respectively. However, these biotherapeutics can be difficult and expensive to manufacture with multiple challenges from recombinant protein generation through to the development of long-acting formulations required to improve the circulating half-life of the drug. In this review, we summarize methodologies and approaches used for making and purifying recombinant GH and GHA proteins, and strategies to improve pharmacokinetic and pharmacodynamic properties, including PEGylation and fusion proteins. Therapeutics that are in clinical use or are currently under development are also discussed.

To date, more than 37 amyloidogenic proteins have been found to form toxic aggregates that are implicated in the progression of numerous debilitating protein misfolding diseases including Alzheimer's disease (AD). Extensive literature highlights the role of β-amyloid (Aβ) aggregates in causing excessive neuronal cell loss in the brains of AD patients. In fact, major advances in our understanding of Aβ aggregation process, including kinetics, toxicity, and structures of fibrillar aggregates have been revealed by examining in vitro preparations of synthetic Aβ peptides. However, ongoing research shows that brain-derived Aβ aggregates have specific characteristics that distinguish them from in vitro prepared species. Notably, the molecular structures of amyloid fibrils grown in the human brain were found to be markedly different than synthetic Aβ fibrils. In addition, recent findings report the existence of heterogeneous Aβ proteoforms in AD brain tissue in contrast to synthetically produced full-length aggregates. Despite their high relevance to AD progression, brain-derived Aβ species are less well-characterized compared with synthetic aggregates. The aim of this review is to provide an overview of the literature on brain-derived Aβ aggregates with particular focus on recent studies that report their structures as well as pathological roles in AD progression. The main motivation of this review is to highlight the importance of utilizing brain-derived amyloids for characterizing the structural and toxic effects of amyloid species. With this knowledge, brain-derived aggregates can be adopted to identify more relevant drug targets and validate potent aggregation inhibitors toward designing highly effective therapeutic strategies against AD.

Because steric crowding is most effective when the crowding agent is similar in size to the molecule that it acts upon and the average macromolecule inside cells is much larger than a small protein or peptide, steric crowding is not predicted to affect their folding inside cells. On the other hand, chemical interactions should perturb in-cell structure and stability because they arise from interactions between the surface of the small protein or peptide and its environment. Indeed, previous in vitro measurements of the λ-repressor fragment, λ6-85 , in crowding matrices comprised of Ficoll or protein crowders support these predictions. Here, we directly quantify the in-cell stability of λ6-85 and distinguish the contribution of steric crowding and chemical interactions to its stability. Using a FRET-labeled λ6-85 construct, we find that the fragment is stabilized by 5°C in-cells compared to in vitro. We demonstrate that this stabilization cannot be explained by steric crowding because, as anticipated, Ficoll has no effect on λ6-85 stability. We find that the in-cell stabilization arises from chemical interactions, mimicked in vitro by mammalian protein extraction reagent (M-PER™). Comparison between FRET values in-cell and in Ficoll confirms that U-2 OS cytosolic crowding is reproduced at macromolecule concentrations of 15% w/v. Our measurements validate the cytomimetic of 15% Ficoll and 20% M-PER™ that we previously developed for protein and RNA folding studies. However, because the in-cell stability of λ6-85 is reproduced by 20% v/v M-PER™ alone, we predict that this simplified mixture could be a useful tool to predict the in-cell behaviors of other small proteins and peptides.

The cosolvent 2,2,2-trifluoroethanol (TFE) is often used to mimic protein desiccation. We assessed the effects of TFE on cytosolic abundant heat soluble protein D (CAHS D) from tardigrades. CAHS D is a member of a unique protein class that is necessary and sufficient for tardigrades to survive desiccation. We find that the response of CAHS D to TFE depends on the concentration of both species. Dilute CAHS D remains soluble and, like most proteins exposed to TFE, gains α-helix. More concentrated solutions of CAHS D in TFE accumulate β-sheet, driving both gel formation and aggregation. At even higher TFE and CAHS D concentrations, samples phase separate without aggregation or increases in helix. Our observations show the importance of considering protein concentration when using TFE.

Protein aggregation results in an array of different size soluble oligomers and larger insoluble fibrils. Insoluble fibrils were originally thought to cause neuronal cell deaths in neurodegenerative diseases due to their prevalence in tissue samples and disease models. Despite recent studies demonstrating the toxicity associated with soluble oligomers, many therapeutic strategies still focus on fibrils or consider all types of aggregates as one group. Oligomers and fibrils require different modeling and therapeutic strategies, targeting the toxic species is crucial for successful study and therapeutic development. Here, we review the role of different-size aggregates in disease, and how factors contributing to aggregation (mutations, metals, post-translational modifications, and lipid interactions) may promote oligomers opposed to fibrils. We review two different computational modeling strategies (molecular dynamics and kinetic modeling) and how they are used to model both oligomers and fibrils. Finally, we outline the current therapeutic strategies targeting aggregating proteins and their strengths and weaknesses for targeting oligomers versus fibrils. Altogether, we aim to highlight the importance of distinguishing the difference between oligomers and fibrils and determining which species is toxic when modeling and creating therapeutics for protein aggregation in disease.

Control of eukaryotic cellular function is heavily reliant on the phosphorylation of proteins at specific amino acid residues, such as serine, threonine, tyrosine, and histidine. Protein kinases that are responsible for this process comprise one of the largest families of evolutionarily related proteins. Dysregulation of protein kinase signaling pathways is a frequent cause of a large variety of human diseases including cancer, autoimmune, neurodegenerative, and cardiovascular disorders. In this study, we mapped all pathogenic mutations in 497 human protein kinase domains from the ClinVar database to the reference structure of Aurora kinase A (AURKA) and grouped them by the relevance to the disease type. Our study revealed that the majority of mutation hotspots associated with cancer are situated within the catalytic and activation loops of the kinase domain, whereas non-cancer-related hotspots tend to be located outside of these regions. Additionally, we identified a hotspot at residue R371 of the AURKA structure that has the highest number of exclusively non-cancer-related pathogenic mutations (21) and has not been previously discussed.

De novo protein design enhances our understanding of the principles that govern protein folding and interactions, and has the potential to revolutionize biotechnology through the engineering of novel protein functionalities. Despite recent progress in computational design strategies, de novo design of protein structures remains challenging, given the vast size of the sequence-structure space. AlphaFold2 (AF2), a state-of-the-art neural network architecture, achieved remarkable accuracy in predicting protein structures from amino acid sequences. This raises the question whether AF2 has learned the principles of protein folding sufficiently for de novo design. Here, we sought to answer this question by inverting the AF2 network, using the prediction weight set and a loss function to bias the generated sequences to adopt a target fold. Initial design trials resulted in de novo designs with an overrepresentation of hydrophobic residues on the protein surface compared to their natural protein family, requiring additional surface optimization. In silico validation of the designs showed protein structures with the correct fold, a hydrophilic surface and a densely packed hydrophobic core. In vitro validation showed that 7 out of 39 designs were folded and stable in solution with high melting temperatures. In summary, our design workflow solely based on AF2 does not seem to fully capture basic principles of de novo protein design, as observed in the protein surface's hydrophobic vs. hydrophilic patterning. However, with minimal post-design intervention, these pipelines generated viable sequences as assessed experimental characterization. Thus, such pipelines show the potential to contribute to solving outstanding challenges in de novo protein design.

The BRICHOS protein superfamily is a diverse group of proteins associated with a wide variety of human diseases, including respiratory distress, COVID-19, dementia, and cancer. A key characteristic of these proteins-besides their BRICHOS domain present in the ER lumen/extracellular part-is that they harbor an aggregation-prone region, which the BRICHOS domain is proposed to chaperone during biosynthesis. All so far studied BRICHOS domains modulate the aggregation pathway of various amyloid-forming substrates, but not all of them can keep denaturing proteins in a folding-competent state, in a similar manner as small heat shock proteins. Current evidence suggests that the ability to interfere with the aggregation pathways of substrates with entirely different end-point structures is dictated by BRICHOS quaternary structure as well as specific surface motifs. This review aims to provide an overview of the BRICHOS protein family and a perspective of the diverse molecular chaperone-like functions of various BRICHOS domains in relation to their structure and conformational plasticity. Furthermore, we speculate about the physiological implication of the diverse molecular chaperone functions and discuss the possibility to use the BRICHOS domain as a blood-brain barrier permeable molecular chaperone treatment of protein aggregation disorders.

The fact that proteins can have their chain formed in a knot is known for almost 30 years. However, as they are not common, only a fraction of such proteins is available in the Protein Data Bank. It was not possible to assess their importance and versatility up until now because we did not have access to the whole proteome of an organism, let alone a human one. The arrival of efficient machine learning methods for protein structure prediction, such as AlphaFold and RoseTTaFold, changed that. We analyzed all proteins from the human proteome (over 20,000) determined with AlphaFold in search for knots and found them in less than 2% of the structures. Using a variety of methods, including homolog search, clustering, quality assessment, and visual inspection, we determined the nature of each of the knotted structures and classified it as either knotted, potentially knotted, or an artifact, and deposited all of them in a database available at: https://knotprot.cent.uw.edu.pl/alphafold. Overall, we found 51 credible knotted proteins (0.2% of human proteome). The set of potentially knotted structures includes a new complex type of a knot not reported in proteins yet. That knot type, denoted 63 in mathematical notation, would necessitate a more complex folding path than any knotted protein characterized to date.

Small Ubiquitin-like Modifier 1 (SUMO1) is an essential protein for many cellular functions, including regulation, signaling, etc., achieved by a process known as SUMOylation, which involves covalent attachment of SUMO1 to target proteins. SUMO1 also regulates the function of several proteins via non-covalent interactions involving the hydrophobic patch in the target protein identified as SUMO Binding or Interacting Motif (SBM/SIM). Here, we demonstrate a crucial functional potential of SUMO1 mediated by its non-covalent interactions with α-Synuclein, a protein responsible for many neurodegenerative diseases called α-Synucleinopathies. SUMO1 hinders the fibrillation of α-Synuclein, an intrinsically disordered protein (IDP) that undergoes a transition to β-structures during the fibrillation process. Using a plethora of biophysical techniques, we show that SUMO1 transiently binds to the N-terminus region of α-Synuclein non-covalently and causes structural compaction, which hinders the self-association process and thereby delays the fibrillation process. On the one hand, this study demonstrates an essential functional role of SUMO1 protein concerning neurodegeneration; it also illustrates the commonly stated mechanism that IDPs carry out multiple functions by structural adaptation to suit specific target proteins, on the other. Residue-level details about the SUMO1-α-Synuclein interaction obtained here also serve as a reliable approach for investigating the detailed mechanisms of IDP functions.

Neutral mutational drift is an important source of biological diversity that remains underexploited in fundamental studies of protein biophysics. This study uses a synthetic transcriptional circuit to study neutral drift in protein tyrosine phosphatase 1B (PTP1B), a mammalian signaling enzyme for which conformational changes are rate limiting. Kinetic assays of purified mutants indicate that catalytic activity, rather than thermodynamic stability, guides enrichment under neutral drift, where neutral or mildly activating mutations can mitigate the effects of deleterious ones. In general, mutants show a moderate activity-stability tradeoff, an indication that minor improvements in the activity of PTP1B do not require concomitant losses in its stability. Multiplexed sequencing of large mutant pools suggests that substitutions at allosterically influential sites are purged under biological selection, which enriches for mutations located outside of the active site. Findings indicate that the positional dependence of neutral mutations within drifting populations can reveal the presence of allosteric networks and illustrate an approach for using synthetic transcriptional systems to explore these mutations in regulatory enzymes.

A single gene yields many forms of proteins via combinations of posttranscriptional/posttranslational modifications. Proteins also fold into higher-order structures and interact with other molecules. The combined molecular diversity leads to the heterogeneity of proteins that manifests as distinct phenotypes. Structural biology has generated vast amounts of data, effectively enabling accurate structural prediction by computational methods. However, structures are often obtained heterologously under homogeneous states in vitro. The lack of native heterogeneity under cellular context creates challenges in precisely connecting the structural data to phenotypes. Mass spectrometry (MS) based proteomics methods can profile proteome composition of complex biological samples. Most MS methods follow the "bottom-up" approach, which denatures and digests proteins into short peptide fragments for ease of detection. Coupled with chemical biology approaches, higher-order structures can be probed via incorporation of covalent labels on native proteins that are maintained at the peptide level. Alternatively, native MS follows the "top-down" approach and directly analyzes intact proteins under nondenaturing conditions. Various tandem MS activation methods can dissect the intact proteins for in-depth structural elucidation. Herein, we review recent native MS applications for characterizing heterogeneous samples, including proteins binding to mixtures of ligands, homo/hetero-complexes with varying stoichiometry, intrinsically disordered proteins with dynamic conformations, glycoprotein complexes with mixed modification states, and active membrane protein complexes in near-native membrane environments. We summarize the benefits, challenges, and ongoing developments in native MS, with the hope to demonstrate an emerging technology that complements other tools by filling the knowledge gaps in understanding the molecular heterogeneity of proteins.

Present understanding of protein structure dynamics trails behind that of static structures. A torsion-angle-based approach, called the representation of protein entities, derives an interpretable conformational space that correlates with data collection temperature, resolution, and reaction coordinate. For more complex systems, atomic coordinates fail to separate functional conformational states, which are still preserved by torsion angle-derived space. This indicates that torsion angles are often a more sensitive and biologically relevant descriptor for protein conformational dynamics than atomic coordinates.

Within the cell, the trace element molybdenum (Mo) is only biologically active when complexed either within the nitrogenase-specific FeMo cofactor or within the molybdenum cofactor (Moco). Moco consists of an organic part, called molybdopterin (MPT) and an inorganic part, that is, the Mo-center. The enzyme which catalyzes the Mo-center formation is the molybdenum insertase (Mo-insertase). Mo-insertases consist of two functional domains called G- and E-domain. The G-domain catalyzes the formation of adenylated MPT (MPT-AMP), which is the substrate for the E-domain, that catalyzes the actual molybdate insertion reaction. Though the functions of E- and G-domain have been elucidated to great structural and mechanistic detail, their combined function is poorly characterized. In this work, we describe a structural model of the eukaryotic Mo-insertase Cnx1 complex that was generated based on cross-linking mass spectrometry combined with computational modeling. We revealed Cnx1 to form an asymmetric hexameric complex which allows the E- and G-domain active sites to align in a catalytic productive orientation toward each other.

We propose a high-throughput method for quantitatively measuring hundreds of protein-peptide binding affinities in parallel. In this assay, a solution of protein is dialyzed into a buffer containing a pool of potential binding peptides, such that upon equilibration the relative abundance of a peptide species is mathematically related to that peptide's dissociation constant, Kd . We use isobaric multiplexed quantitative proteomics to simultaneously determine the relative abundance, and hence the Kd and its associated error, for an entire peptide library. We apply this technique, which we call PEDAL (parallel equilibrium dialysis for affinity learning), to determine accurate Kd 's between a PDZ domain and hundreds of peptides, spanning an affinity range of multiple orders of magnitude in a single experiment. PEDAL is a convenient, fast, and low-cost method for measuring large numbers of protein-peptide affinities in parallel, providing a rare combination of true in-solution binding equilibria with the ability to multiplex.

In the decades following the discovery that genes encode proteins, scientists have tried to exhaustively and comprehensively characterize the human genome. Recent advances in computational methods along with transcriptomic and proteomic techniques have now shown that historically non-coding genomic regions may contain non-canonical open reading frames (ncORFs), which may encode functional miniproteins or otherwise exert regulatory activity through coding-independent functions. Increasingly, it is clear that these ncORFs may play critical roles in major human diseases such as cancer. In this review, we summarize the history and current progress of ncORF research and explore the known functions of ncORFs and the miniproteins they may encode. We particularly highlight the emerging body of evidence supporting a role for ncORFs and miniproteins contributions in cancer. Finally, we provide a blueprint for high-priority areas of future research for ncORFs in cancer, focusing on ncORF detection, functional characterization, and therapeutic intervention.

Comment in Protein Sci. 2023 Apr;32(4):e4601.

Cyanophycin is a natural polymer composed of a poly-aspartate backbone with arginine attached to each of the aspartate sidechains. Produced by a wide range of bacteria, which mainly use it as a store of fixed nitrogen, it has many promising industrial applications. Cyanophycin can be synthesized from the amino acids Asp and Arg by the widespread cyanophycin synthetase 1 (CphA1), or from the dipeptide β-Asp-Arg by the cyanobacterial enzyme cyanophycin synthetase 2 (CphA2). CphA2 enzymes display a range of oligomeric states, from dimers to dodecamers. Recently, the crystal structure of a CphA2 dimer was solved but could not be obtained in complex with substrate. Here, we report cryo-EM structures of the hexameric CphA2 from Stanieria sp. at ~2.8 Å resolution, both with and without ATP analog and cyanophycin. The structures show a two-fold symmetrical, trimer-of-dimers hexameric architecture, and substrate-binding interactions that are similar to those of CphA1. Mutagenesis experiments demonstrate the importance of several conserved substrate-binding residues. We also find that a Q416A/R528G double mutation prevents hexamer formation and use this double mutant to show that hexamerization augments the rate of cyanophycin synthesis. Together, these results increase our mechanistic understanding of how an interesting green polymer is biosynthesized.

Numerous age-linked diseases are rooted in protein misfolding; this has motivated the development of small molecules and therapeutic antibodies that target the aggregation of disease-linked proteins. Here we explore another approach: molecular chaperones with engineerable protein scaffolds such as the ankyrin repeat domain (ARD). We tested the ability of cpSRP43, a small, robust, ATP- and cofactor-independent plant chaperone built from an ARD, to antagonize disease-linked protein aggregation. cpSRP43 delays the aggregation of multiple proteins including the amyloid beta peptide (Aβ) associated with Alzheimer's disease and α-synuclein associated with Parkinson's disease. Kinetic modeling and biochemical analyses show that cpSRP43 targets early oligomers during Aβ aggregation, preventing their transition to a self-propagating nucleus on the fibril surface. Accordingly, cpSRP43 rescued neuronal cells from the toxicity of extracellular Aβ42 aggregates. The substrate-binding domain of cpSRP43, composed primarily of the ARD, is necessary and sufficient to prevent Aβ42 aggregation and protect cells against Aβ42 toxicity. This work provides an example in which an ARD chaperone non-native to mammalian cells harbors anti-amyloidal activity, which may be exploited for bioengineering.

Protein regions which are intrinsically disordered, exist as an ensemble of rapidly interconverting structures. Cooling proteins to cryogenic temperatures for dynamic nuclear polarization (DNP) magic angle spinning (MAS) NMR studies suspends most of the motions, resulting in peaks that are broad but not featureless. To demonstrate that detailed conformational restraints can be retrieved from the peak shapes of frozen proteins alone, we developed and used a simulation framework to assign peak features to conformers in the ensemble. We validated our simulations by comparing them to spectra of α-synuclein acquired under different experimental conditions. Our assignments of peaks to discrete dihedral angle populations suggest that structural constraints are attainable under cryogenic conditions. The ability to infer ensemble populations from peak shapes has important implications for DNP MAS NMR studies of proteins with regions of disorder in living cells because chemical shifts are the most accessible measured parameter.

Molecular chaperones are key components of protein quality control system, which plays an essential role in controlling protein homeostasis. Aha1 has been identified as a co-chaperone of Hsp90 known to strongly accelerate Hsp90's ATPase activity. Meanwhile, it is reported that Aha1 could also act as an autonomous chaperone and protect stressed or disordered proteins from aggregation. Here, in this article, a series of in vitro experiments were conducted to verify whether Aha1 has a non-Hsp90-dependent holdase activity and to elucidate the associated molecular mechanism for substrate recognition. According to the results of the refolding assay, the highly conserved N-terminal extension spanning M1 to R16 in Aha1 from higher eukaryotes is responsible for the holdase activity of the protein. As revealed by the NMR data, Aha1's N-terminal extension mainly adopts a disordered conformation in solution and shows no tight contacts with the core structure of Aha1's N-terminal domain. Based on the intrinsically disordered structure feature and the primary sequence of Aha1's N-terminal extension, the fuzzy-type protein-protein interactions involving this specific region and the unfolded substrate proteins are expected. The following mutation analysis data demonstrated that the Van der Waals contacts potentially involving two tryptophans including W4 and W11 do not play a dominant role in the interaction between Aha1 and unfolded maltose binding protein (MBP). Meanwhile, since the high concentration of NaCl could abolish the holdase activity of Aha1, the electrostatic interactions mediated by those charged residues in Aha1's N-terminal extension are thus indicated to play a crucial role in the substrate recognition.

Multisite λ-dynamics (MSλD) is a novel method for the calculation of relative free energies of binding for ligands to their targeted receptors. It can be readily used to examine a large number of molecules with multiple functional groups at multiple sites around a common core. This makes MSλD a powerful tool in structure-based drug design. In the present study, MSλD is applied to calculate the relative binding free energies of 1296 inhibitors to the testis specific serine kinase 1B (TSSK1B), a validated target for male contraception. For this system, MSλD requires significantly fewer computational resources compared to traditional free energy methods like free energy perturbation or thermodynamic integration. From MSλD simulations, we examined whether modifications of a ligand at two different sites are coupled or not. Based on our calculations, we established a quantitative structure-activity relationship (QSAR) for this set of molecules and identified a site in the ligand where further modification, such as adding more polar groups, may lead to increased binding affinity.

There is ample computational, but only sparse experimental data suggesting that pico-ns motions with 1 Å amplitude are pervasive in proteins in solution. Such motions, if present in reality, must deeply affect protein function and protein entropy. Several NMR relaxation experiments have provided insights into motions of proteins in solution, but they primarily report on azimuthal angle variations of vectors of covalently-linked atoms. As such, these measurements are not sensitive to distance fluctuations, and cannot but under-represent the dynamical properties of proteins. Here we analyze a novel NMR relaxation experiment to measure amide proton transverse relaxation rates in uniformly 15 N labeled proteins, and present results for protein domain GB1 at 283 and 303 K. These relaxation rates depend on fluctuations of dipolar interactions between 1 HN and many nearby protons on both the backbone and sidechains. Importantly, they also report on fluctuations in the distances between these protons. We obtained a large mismatch between rates computed from the crystal structure of GB1 and the experimental rates. But when the relaxation rates were calculated from a 200 ns molecular dynamics trajectory using a novel program suite, we obtained a substantial improvement in the correspondence of experimental and theoretical rates. As such, this work provides novel experimental evidence of widespread motions in proteins. Since the improvements are substantial, but not sufficient, this approach may also present a new benchmark to help improve the theoretical forcefields underlying the molecular dynamics calculations.

Methylenetetrahydrofolate reductase (MTHFR) is a key metabolic enzyme in colonization and virulence of Neisseria meningitidis, a causative agent of meningococcal diseases. Here, the biochemical and structural properties of MTHFR from a virulent strain of N. meningitidis serogroup B (NmMTHFR) were characterized. Unlike other orthologs, NmMTHFR functions as a unique homohexamer, composed of three homo-dimerization partners, as shown in our 2.7 Å resolution crystal structure. Six active sites were formed solely within monomers and located away from the oligomerization interfaces. Flavin adenine dinucleotide cofactor formed hydrogen bonds with conserved sidechains, positioning its isoalloxazine ring adjacent to the overlapping binding sites of nicotinamide adenine dinucleotide (NADH) coenzyme and CH2 -H4 folate substrate. NmMTHFR utilized NADH (Km  = 44 μM) as an electron donor in the NAD(P)H-CH2 -H4 folate oxidoreductase assay, but not nicotinamide adenine dinucleotide phosphate (NADPH) which is the donor required in human MTHFR. In silico analysis and mutagenesis studies highlighted the significant difference in orientation of helix α7A (Phe215-Thr225) with that in the human enzyme. The extended sidechain of Met221 on helix α7A plays a role in stabilizing the folded structure of NADH in the hydrophobic box. This supports the NADH specificity by restricting the phosphate group of NADPH that causes steric clashes with Glu26. The movement of Met221 sidechain allows the CH2 -H4 folate substrate to bind. The unique topology of its NADH and CH2 -H4 folate binding pockets makes NmMTHFR a promising drug target for the development of new antimicrobial agents that may possess reduced off-target side effects.

We recently converted the GAF domain of NpR3784 cyanobacteriochrome into near-infrared (NIR) fluorescent proteins (FPs). Unlike cyanobacterichrome, which incorporates phycocyanobilin tetrapyrrole, engineered NIR FPs bind biliverdin abundant in mammalian cells, thus being the smallest scaffold for it. Here, we determined the crystal structure of the brightest blue-shifted protein of the series, miRFP670nano3, at 1.8 Å resolution, characterized its chromophore environment and explained the molecular basis of its spectral properties. Using the determined structure, we have rationally designed a red-shifted NIR FP, termed miRFP704nano, with excitation at 680 nm and emission at 704 nm. miRFP704nano exhibits a small size of 17 kDa, enhanced molecular brightness, photostability and pH-stability. miRFP704nano performs well in various protein fusions in live mammalian cells and should become a versatile genetically-encoded NIR probe for multiplexed imaging across spatial scales in different modalities.

Parkinson's disease can manifest either as a sporadic form, which is common, or as an inherited autosomal dominant trait resulting from missense mutations. Recently, the novel α-synuclein variant V15A was identified in two Caucasian and two Japanese families with Parkinson's disease. Using a combination of NMR spectroscopy, membrane binding assays and aggregation assays we show that the V15A mutation does not strongly perturb the conformational ensemble of monomeric α-synuclein in solution, but weakens its affinity for membranes. Attenuated membrane binding raises the concentration of the aggregation-prone disordered α-synuclein in solution, allowing only the V15A variant but not wild-type α-synuclein to form amyloid fibrils in the presence of liposomes. These findings, together with earlier research on other missense mutations of α-synuclein, suggest that maintaining a balance between membrane-bound and free aggregation-competent α-synuclein is critical in α-synucleinopathies.

S-adenosyl methionine (SAM)-dependent methyl transferases (MTases) are a ubiquitous class of enzymes catalyzing dozens of essential life processes. Despite targeting a large space of substrates with diverse intrinsic reactivity, SAM MTases have similar catalytic efficiency. While understanding of MTase mechanism has grown tremendously through the integration of structural characterization, kinetic assays, and multiscale simulations, it remains elusive how these enzymes have evolved to fit the diverse chemical needs of their respective substrates. In this work, we performed a high-throughput molecular modeling analysis of 91 SAM MTases to better understand how their properties (i.e., electric field [EF] strength and active site volumes) help achieve similar catalytic efficiency toward substrates of different reactivity. We found that EF strengths have largely adjusted to make the target atom a better methyl acceptor. For MTases that target RNA/DNA and histone proteins, our results suggest that EF strength accommodates formal hybridization state and variation in cavity volume trends with diversity of substrate classes. Metal ions in SAM MTases contribute negatively to EF strength for methyl donation and enzyme scaffolds tend to offset these contributions.

The Immune Epitope Database (IEDB) catalogs T cell, B cell, and major histocompatibility complex ligand information in the context of infectious disease, allergy, autoimmunity, and transplantation. An important component of this information is three-dimensional structural data on T cell receptors, antibodies, and pairwise residue interactions between immune receptors and antigens, which we refer to as IEDB-3D. Such data is highly valuable for mechanically understanding receptor:ligand interactions. Here, we present IEDB-3D 2.0, which comprises a complete overhaul of how we obtain and present 3D structural data. A new 3D viewer experience that utilizes iCn3D has been implemented to replace outdated java-based technology. In addition, we have designed a new epitope mapping system that matches each epitope available in the IEDB with its antigen structural data. Finally, immunogenicity data retrieved from the IEDB's ImmunomeBrowser can now be used to highlight immunogenic regions of an antigen directly in iCn3D. Overall, the IEDB-3D 2.0 provides an updated tool platform to visualize epitope data cataloged in the IEDB.

Palladin is an actin binding protein that is specifically upregulated in metastatic cancer cells but also colocalizes with actin stress fibers in normal cells and is critical for embryonic development as well as wound healing. Of nine isoforms present in humans, only the 90 kDa isoform of palladin, comprising three immunoglobulin (Ig) domains and one proline-rich region, is ubiquitously expressed. Previous work has established that the Ig3 domain of palladin is the minimal binding site for F-actin. In this work, we compare functions of the 90 kDa isoform of palladin to the isolated actin binding domain. To understand the mechanism of action for how palladin can influence actin assembly, we monitored F-actin binding and bundling as well as actin polymerization, depolymerization, and copolymerization. Together, these results demonstrate that there are key differences between the Ig3 domain and full-length palladin in actin binding stoichiometry, polymerization, and interactions with G-actin. Understanding the role of palladin in regulating the actin cytoskeleton may help us develop means to prevent cancer cells from reaching the metastatic stage of cancer progression.

PADI4 is one of the human isoforms of a family of enzymes implicated in the conversion of arginine to citrulline. MDM2 is an E3 ubiquitin ligase which is crucial for down-regulation of degradation of the tumor suppressor gene p53. Given the relationship between both PADI4 and MDM2 with p53-signaling pathways, we hypothesized they may interact directly, and this interaction could be relevant in the context of cancer. Here, we showed their association in the nucleus and cytosol in several cancer cell lines. Furthermore, binding was hampered in the presence of GSK484, an enzymatic PADI4 inhibitor, suggesting that MDM2 could bind to the active site of PADI4, as confirmed by in silico experiments. In vitro and in silico studies showed that the isolated N-terminal region of MDM2, N-MDM2, interacted with PADI4, and residues Thr26, Val28, Phe91 and Lys98 were more affected by the presence of the enzyme. Moreover, the dissociation constant between N-MDM2 and PADI4 was comparable to the IC50 of GSK484 from in cellulo experiments. The interaction between MDM2 and PADI4 might imply MDM2 citrullination, with potential therapeutic relevance for improving cancer treatment, due to the generation of new antigens.

Slow conformational changes are often directly linked to protein function. It is however less clear how such processes may perturb the overall folding stability of a protein. We previously found that the stabilizing double mutant L49I/I57V in the small protein chymotrypsin inhibitor 2 from barley led to distributed increased nanosecond and faster dynamics. Here we asked what effects the L49I and I57V substitutions, either individually or together, have on the slow conformational dynamics of CI2. We used 15 N CPMG spin relaxation dispersion experiments to measure the kinetics, thermodynamics, and structural changes associated with slow conformational change in CI2. These changes result in an excited state that is populated to 4.3% at 1°C. As the temperature is increased the population of the excited state decreases. Structural changes in the excited state are associated with residues that interact with water molecules that have well defined positions and are found at these positions in all crystal structures of CI2. The substitutions in CI2 have only little effect on the structure of the excited state whereas the stability of the excited state to some extent follows the stability of the main state. The minor state is thus most populated for the most stable CI2 variant and least populated for the least stable variant. We hypothesize that the interactions between the substituted residues and the well-ordered water molecules links subtle structural changes around the substituted residues to the region in the protein that experience slow conformational changes.

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) nucleocapsid protein is the most abundantly expressed viral protein during infection where it targets both RNA and host proteins. However, identifying how a single viral protein interacts with so many different targets remains a challenge, providing the impetus here for identifying the interaction sites through multiple methods. Through a combination of nuclear magnetic resonance (NMR), electron microscopy, and biochemical methods, we have characterized nucleocapsid interactions with RNA and with three host proteins, which include human cyclophilin-A, Pin1, and 14-3-3τ. Regarding RNA interactions, the nucleocapsid protein N-terminal folded domain preferentially interacts with smaller RNA fragments relative to the C-terminal region, suggesting an initial RNA engagement is largely dictated by this N-terminal region followed by weaker interactions to the C-terminal region. The nucleocapsid protein forms 10 nm ribonuclear complexes with larger RNA fragments that include 200 and 354 nucleic acids, revealing its potential diversity in sequestering different viral genomic regions during viral packaging. Regarding host protein interactions, while the nucleocapsid targets all three host proteins through its serine-arginine-rich region, unstructured termini of the nucleocapsid protein also engage host cyclophilin-A and host 14-3-3τ. Considering these host proteins play roles in innate immunity, the SARS-CoV-2 nucleocapsid protein may block the host response by competing interactions. Finally, phosphorylation of the nucleocapsid protein quenches an inherent dynamic exchange process within its serine-arginine-rich region. Our studies identify many of the diverse interactions that may be important for SARS-CoV-2 pathology during infection.

Small molecules that regulate protein-protein interactions can be valuable drugs; however, the development of such small molecules is challenging as the molecule must interfere with an interaction that often involves a large surface area. Herein, we propose that modulating the conformational ensemble of the proteins participating in a given interaction, rather than blocking the interaction by directly binding to the interface, is a relevant strategy for interfering with a protein-protein interaction. In this study, we applied this concept to P-cadherin, a cell surface protein forming homodimers that are essential for cell-cell adhesion in various biological contexts. We first determined the crystal structure of P-cadherin with a small molecule inhibitor whose inhibitory mechanism was unknown. Molecular dynamics simulations suggest that the inhibition of cell adhesion by this small molecule results from modulation of the conformational ensemble of P-cadherin. Our study demonstrates the potential of small molecules altering the conformation ensemble of a protein as inhibitors of biological relevant protein-protein interactions.

DJ-1, a protein encoded by PARK7 plays a protective role against neurodegeneration. Since its glyoxalase III activity catalyzing methylglyoxal (MG) to lactate was discovered, DJ-1 has been re-established as a deglycase decomposing the MG-intermediates with amino acids and nucleotides (hemithioacetal and hemiaminal) rather than MG itself, but it is still debatable. Here, we have clarified that human DJ-1 directly recognizes MG, and not MG-intermediates, by monitoring the detailed catalytic processes and enantiomeric lactate products. The hemithioacetal intermediate between C106 of 15 N-labeled DJ-1 (15N DJ-1) and MG was also monitored by NMR. TRIS molecule formed stable diastereotopic complexes with MG (Kd , 1.57 ± 0.27 mM) by utilizing its three OH groups, which likely disturbed the assay of deglycase activity. The low kcat of DJ-1 for MG and its MG-induced structural perturbation may suggest that DJ-1 has a regulatory function as an in vivo sensor of reactive carbonyl stress.

DisProt is the primary repository of Intrinsically Disordered Proteins (IDPs). This database is manually curated and the annotations there have strong experimental support. Currently, DisProt contains a relatively small number of proteins highlighting the importance of transferring annotations regarding verified disorder state and corresponding functions to homologous proteins in other species. In such a way, providing them with highly valuable information to better understand their biological roles. While the principles and practicalities of homology transfer are well-established for globular proteins, these are largely lacking for disordered proteins. We used DisProt to evaluate the transferability of the annotation terms to orthologous proteins. For each protein, we looked for their orthologs, with the assumption that they will have a similar function. Then, for each protein and their orthologs, we made multiple sequence alignments (MSAs). Disordered sequences are fast evolving and can be hard to align, therefore, we implemented alignment quality control steps ensuring robust alignments before mapping the annotations. We have designed a pipeline to obtain good-quality MSAs and to transfer annotations from any protein to their orthologs. Applying the pipeline to DisProt proteins, from the 1731 entries with 5623 annotations, we can reach 97,555 orthologs and transfer a total of 301,190 terms by homology. We also provide a web server for consulting the results of DisProt proteins and execute the pipeline for any other protein. The server Homology Transfer IDP (HoTIDP) is accessible at http://hotidp.leloir.org.ar.

Nicotinic acetylcholine receptors (N-AChRs) mediate fast synaptic signaling and are members of the pentameric ligand-gated ion channel (pLGIC) family. They rely on a network of accessory proteins in vivo for correct formation and transport to the cell surface. Resistance to cholinesterase 3 (RIC-3) is an endoplasmic reticulum protein that physically interacts with nascent pLGIC subunits and promotes their oligomerization. It is not known why some N-AChRs require RIC-3 in heterologous expression systems, whereas others do not. Previously we reported that the ACR-16 N-AChR from the parasitic nematode Dracunculus medinensis does not require RIC-3 in Xenopus laevis oocytes. This is unusual because all other nematode ACR-16, like the closely related Ascaris suum ACR-16, require RIC-3. Their high sequence similarity limits the number of amino acids that may be responsible, and the goal of this study was to identify them. A series of chimeras and point mutations between A. suum and D. medinensis ACR-16, followed by functional characterization with electrophysiology, identified two residues that account for a majority of the receptor requirement for RIC-3. ACR-16 with R/K159 in the cys-loop and I504 in the C-terminal tail did not require RIC-3 for functional expression. Mutating either of these to R/K159E or I504T, residues found in other nematode ACR-16, conferred a RIC-3 requirement. Our results agree with previous studies showing that these regions interact and are involved in receptor synthesis. Although it is currently unclear what precise mechanism they regulate, these residues may be critical during specific subunit folding and/or assembly cascades that RIC-3 may promote.

Cellobiose dehydrogenase (CDH) is a bioelectrocatalyst that enables direct electron transfer (DET) in biosensors and biofuel cells. The application of this bidomain hemoflavoenzyme for physiological glucose measurements is limited by its acidic pH optimum and slow interdomain electron transfer (IET) at pH 7.5. The reason for this rate-limiting electron transfer step is electrostatic repulsion at the interface between the catalytic dehydrogenase domain and the electron mediating cytochrome domain (CYT). We applied rational interface engineering to accelerate the IET for the pH prevailing in blood or interstitial fluid. Phylogenetic and structural analyses guided the design of 17 variants in which acidic amino acids were mutated at the CYT domain. Five mutations (G71K, D160K, Q174K, D177K, M180K) increased the pH optimum and IET rate. Structure-based analysis of the variants suggested two mechanisms explaining the improvements: electrostatic steering and stabilization of the closed state by hydrogen bonding. Combining the mutations into six combinatorial variants with up to five mutations shifted the pH optimum from 4.5 to 7.0 and increased the IET at pH 7.5 over 12-fold from 0.1 to 1.24 s-1 . While the mutants sustained a high enzymatic activity and even surpassed the IET of the wild-type enzyme, the accumulated positive charges on the CYT domain decreased DET, highlighting the importance of CYT for IET and DET. This study shows that interface engineering is an effective strategy to shift the pH optimum and improve the IET of CDH, but future work needs to maintain the DET of the CYT domain for bioelectronic applications.

In an effort to investigate the molecular determinants of ligand recognition of the C-terminal SH2 domain of the SHP2 protein, we conducted extensive site-directed mutagenesis and kinetic binding experiments with a peptide mimicking a specific portion of a physiological ligand (the scaffold protein Gab2). Obtained data provided an in-depth characterization of the binding reaction, allowing us to pinpoint residues topologically far from the binding pocket of the SH2 domain to have a role in the recognition and binding of the peptide. The presence of a sparse energetic network regulating the interaction with Gab2 was identified and characterized through double mutant cycle analysis, performed by challenging all the designed site-directed variants of C-SH2 with a Gab2 peptide mutated at +3 position relative to its phosphorylated tyrosine, a key residue for C-SH2 binding specificity. Results highlighted non-optimized residues involved in the energetic network regulating the binding with Gab2, which may be at the basis of the ability of this SH2 domain to interact with different partners in the intracellular environment. Moreover, a detailed analysis of kinetic and thermodynamic parameters revealed the role of the residue at +3 position on Gab2 in the early and late events of the binding reaction with the C-SH2 domain.

β-lactam antibiotics target DD-transpeptidases, enzymes that perform the last step of bacterial cell-wall synthesis. To block the antimicrobial activity of these antibiotics, bacteria have evolved lactamases that render them inert. Among these, TEM-1, a class A lactamase, has been extensively studied. In 2004, Horn et al. described a novel allosteric TEM-1 inhibitor, FTA, that binds distant from the TEM-1 orthosteric (penicillin-binding) pocket. TEM-1 has subsequently become a model for the study of allostery. In the present work, we perform molecular dynamics simulations of FTA-bound and FTA-absent TEM-1, totaling ~3 μS, that provide new insight into TEM-1 inhibition. In one of the simulations, bound FTA assumed a conformation different than that observed crystallographically. We provide evidence that the alternate pose is physiologically plausible and describe how it impacts our understanding of TEM-1 allostery.

Antiviral therapeutics to treat SARS-CoV-2 are needed to diminish the morbidity of the ongoing COVID-19 pandemic. A well-precedented drug target is the main viral protease (MPro ), which is targeted by an approved drug and by several investigational drugs. Emerging viral resistance has made new inhibitor chemotypes more pressing. Adopting a structure-based approach, we docked 1.2 billion non-covalent lead-like molecules and a new library of 6.5 million electrophiles against the enzyme structure. From these, 29 non-covalent and 11 covalent inhibitors were identified in 37 series, the most potent having an IC50 of 29 and 20 μM, respectively. Several series were optimized, resulting in low micromolar inhibitors. Subsequent crystallography confirmed the docking predicted binding modes and may template further optimization. While the new chemotypes may aid further optimization of MPro inhibitors for SARS-CoV-2, the modest success rate also reveals weaknesses in our approach for challenging targets like MPro versus other targets where it has been more successful, and versus other structure-based techniques against MPro itself.

The availability of an expanded genetic code opens exciting new opportunities in enzyme design and engineering. In this regard histidine analogues have proven particularly versatile, serving as ligands to augment metalloenzyme function and as catalytic nucleophiles in designed enzymes. The ability to genetically encode multiple functional residues could greatly expand the range of chemistry accessible within enzyme active sites. Here, we develop mutually orthogonal translation components to selectively encode two structurally similar histidine analogues. Transplanting known mutations from a promiscuous Methanosarcina mazei pyrrolysyl-tRNA synthetase (MmPylRSIFGFF ) into a single domain PylRS from Methanomethylophilus alvus (MaPylRSIFGFF ) provided a variant with improved efficiency and specificity for 3-methyl-L-histidine (MeHis) incorporation. The MaPylRSIFGFF clone was further characterized using in vitro biochemical assays and x-ray crystallography. We subsequently engineered the orthogonal MmPylRS for activity and selectivity for 3-(3-pyridyl)-L-alanine (3-Pyr), which was used in combination with MaPylRSIFGFF to produce proteins containing both 3-Pyr and MeHis. Given the versatile roles played by histidine in enzyme mechanisms, we anticipate that the tools developed within this study will underpin the development of enzymes with new and enhanced functions.

The characterization of protein stability is essential for understanding the functions of proteins. Hydroxysteroid dehydrogenase is involved in the biosynthesis of steroid hormones and the detoxification of xenobiotic carbonyl compounds. However, the stability of hydroxysteroid dehydrogenases has not yet been characterized in detail. Here, we determined the changes in Gibbs free energy, enthalpy, entropy, and heat capacity of unfolding for 3α-hydroxysteroid dehydrogenase/carbonyl reductase (3α-HSD/CR) by varying the pH and urea concentration through differential scanning fluorimetry and presented pH-dependent protein stability as a function of temperature. 3α-HSD/CR shows the maximum stability of 30.79 kJ mol-1 at 26.4°C, pH 7.6 and decreases to 7.74 kJ mol-1 at 25.7°C, pH 4.5. The change of heat capacity of 30.25 ± 1.38 kJ mol-1  K-1 is obtained from the enthalpy of denaturation as a function of melting temperature at varied pH. Two proton uptakes are linked to protein unfolding from residues with differential pKa of 4.0 and 6.5 in the native and denatured states, respectively. The large positive heat capacity change indicated that hydrophobic interactions played an important role in the folding of 3α-HSD/CR. These studies reveal the mechanism of protein unfolding in HSD and provide a convenient method to extract thermodynamic parameters for characterizing protein stability using differential scanning fluorimetry.

BiP (immunoglobulin heavy-chain binding protein) is a Hsp70 monomeric ATPase motor that plays broad and crucial roles in maintaining proteostasis inside the cell. Structurally, BiP is formed by two domains, a nucleotide-binding domain (NBD) with ATPase activity connected by a flexible hydrophobic linker to the substrate-binding domain. While the ATPase and substrate binding activities of BiP are allosterically coupled, the latter is also dependent on nucleotide binding. Recent structural studies have provided new insights into BiP's allostery; however, the influence of temperature on the coupling between substrate and nucleotide binding to BiP remains unexplored. Here, we study BiP's binding to its substrate at the single molecule level using thermo-regulated optical tweezers which allows us to mechanically unfold the client protein and explore the effect of temperature and different nucleotides on BiP binding. Our results confirm that the affinity of BiP for its protein substrate relies on nucleotide binding, by mainly regulating the binding kinetics between BiP and its substrate. Interestingly, our findings also showed that the apparent affinity of BiP for its protein substrate in the presence of nucleotides remains invariable over a wide range of temperatures, suggesting that BiP may interact with its client proteins with similar affinities even when the temperature is not optimal. Thus, BiP could play a role as a "thermal buffer" in proteostasis.

Many methodologically diverse computational methods have been applied to the growing challenge of predicting and interpreting the effects of protein variants. As many pathogenic mutations have a perturbing effect on protein stability or intermolecular interactions, one highly interpretable approach is to use protein structural information to model the physical impacts of variants and predict their likely effects on protein stability and interactions. Previous efforts have assessed the accuracy of stability predictors in reproducing thermodynamically accurate values and evaluated their ability to distinguish between known pathogenic and benign mutations. Here, we take an alternate approach, and explore how well stability predictor scores correlate with functional impacts derived from deep mutational scanning (DMS) experiments. In this work, we compare the predictions of 9 protein stability-based tools against mutant protein fitness values from 49 independent DMS datasets, covering 170,940 unique single amino acid variants. We find that FoldX and Rosetta show the strongest correlations with DMS-based functional scores, similar to their previous top performance in distinguishing between pathogenic and benign variants. For both methods, performance is considerably improved when considering intermolecular interactions from protein complex structures, when available. Furthermore, using these two predictors, we derive a "Foldetta" consensus score, which improves upon the performance of both, and manages to match dedicated variant effect predictors in reflecting variant functional impacts. Finally, we also highlight that predicted stability effects show consistently higher correlations with certain DMS experimental phenotypes, particularly those based upon protein abundance, and, in certain cases, can significantly outcompete sequence-based variant effect prediction methodologies for predicting functional scores from DMS experiments.

Huntington's disease (HD) is a neurodegenerative disease resulting from an expansion of the polyglutamine (polyQ) domain within the huntingtin protein (htt). PolyQ expansion triggers toxic aggregation and alters htt/lipid interactions. The first 17 amino acids at the N-terminus of htt (Nt17) have a propensity to form an amphipathic α-helix crucial to aggregation and membrane binding. Htt interacts closely with a variety of membrane systems including those of the endoplasmic reticulum, mitochondria, nuclear envelope, and plasma membrane. Membrane composition heavily influences both htt aggregation and lipid interactions, and cholesterol is a crucial membrane component that modulates properties such as fluidity, permeability, and organization. In HD, cholesterol homeostasis is disrupted, and likely plays a role in toxicity. The objective of these studies was to identify the impact of cholesterol on htt aggregation and lipid interactions in various lipid systems. Lipid systems of POPC, DOPC, and POPG with varied levels of exogenously added cholesterol were exposed to htt, and the influences on aggregation, lipid binding, and htt/lipid complexation were evaluated using thioflavin-T aggregation assays, atomic force microscopy, colorimetric lipid binding assays, and mass spectrometry. The addition of cholesterol to DOPC vesicles enhanced htt aggregation. In the presence of vesicles of either POPC or POPG, the addition of cholesterol reduced htt aggregation. Htt/lipid binding decreased for POPC and increased for both DOPC and POPG with increasing cholesterol content, with observed differences in htt/lipid complexation. Altered cholesterol content influences htt aggregation, lipid binding, and complexation differently depending on overall lipid composition.

We have purified ledodin, a cytotoxic 22-kDa protein from shiitake mushroom (Lentinula edodes) consisting of a 197 amino acid chain. Ledodin possessed N-glycosylase activity on the sarcin-ricin loop of mammalian 28S rRNA and inhibited protein synthesis. However, it was not active against insect, fungal, and bacterial ribosomes. In vitro and in silico studies suggested that ledodin exhibits a catalytic mechanism like that of DNA glycosylases and plant ribosome-inactivating proteins. Moreover, the sequence and structure of ledodin was not related to any protein of known function, although ledodin-homologous sequences were found in the genome of several species of fungi, some edible, belonging to different orders of the class Agaricomycetes. Therefore, ledodin could be the first of a new family of enzymes widely distributed among this class of basidiomycetes. The interest of these proteins lies both, in the fact that they can be a toxic agent of some edible mushrooms and in their application in medicine and biotechnology.

We have developed an algorithm, ParSe, which accurately identifies from the primary sequence those protein regions likely to exhibit physiological phase separation behavior. Originally, ParSe was designed to test the hypothesis that, for flexible proteins, phase separation potential is correlated to hydrodynamic size. While our results were consistent with that idea, we also found that many different descriptors could successfully differentiate between three classes of protein regions: folded, intrinsically disordered, and phase-separating intrinsically disordered. Consequently, numerous combinations of amino acid property scales can be used to make robust predictions of protein phase separation. Built from that finding, ParSe 2.0 uses an optimal set of property scales to predict domain-level organization and compute a sequence-based prediction of phase separation potential. The algorithm is fast enough to scan the whole of the human proteome in minutes on a single computer and is equally or more accurate than other published predictors in identifying proteins and regions within proteins that drive phase separation. Here, we describe a web application for ParSe 2.0 that may be accessed through a browser by visiting https://stevewhitten.github.io/Parse_v2_FASTA to quickly identify phase-separating proteins within large sequence sets, or by visiting https://stevewhitten.github.io/Parse_v2_web to evaluate individual protein sequences.

Unlike laccases sensu stricto, which are usually monomeric enzymes, laccase-like enzymes recently re-classified as Novel Laccases (NLACs) are characterized by the formation of heterodimers with small proteins (subunits) of unknown function. Here the NLAC from Pleurotus eryngii (PeNL) and a small protein selected from the fungal genome, that is homologous to reported POXA3 from Pleurotus ostreatus, were produced in Aspergillus oryzae separately or together. The two proteins interacted regardless of whether the small subunit was co-expressed or exogenously added to the enzyme. The stability and catalytic activity of PeNL was significantly enhanced in the presence of the small subunit. Size exclusion chromatography-multi angle light scattering (SEC-MALS) analysis confirmed that the complex PeNL-ss is a heterodimer of 77.4 kDa. The crystallographic structure of the small protein expressed in Escherichia coli was solved at 1.6 Å resolution. This is the first structure elucidated of a small subunit of a NLAC. The helix bundle structure of the small subunit accommodates well with the enzyme model structure, including interactions with specific regions of NLACs and some amino acid residues of the substrate-binding loops.

ICA512/PTPRN is a receptor tyrosine-like phosphatase implicated in the biogenesis and turnover of the insulin secretory granules (SGs) in pancreatic islet beta cells. Previously we found biophysical evidence that its luminal RESP18 homology domain (RESP18HD) forms a biomolecular condensate and interacts with insulin in vitro at close-to-neutral pH, that is, in conditions resembling those present in the early secretory pathway. Here we provide further evidence for the relevance of these findings by showing that at pH 6.8 RESP18HD interacts also with proinsulin-the physiological insulin precursor found in the early secretory pathway and the major luminal cargo of β-cell nascent SGs. Our light scattering analyses indicate that RESP18HD and proinsulin, but also insulin, populate nanocondensates ranging in size from 15 to 300 nm and 10e2 to 10e6 molecules. Co-condensation of RESP18HD with proinsulin/insulin transforms the initial nanocondensates into microcondensates (size >1 μm). The intrinsic tendency of proinsulin to self-condensate implies that, in the ER, a chaperoning mechanism must arrest its spontaneous intermolecular condensation to allow for proper intramolecular folding. These data further suggest that proinsulin is an early driver of insulin SG biogenesis, in a process in which its co-condensation with RESP18HD participates in their phase separation from other secretory proteins in transit through the same compartments but destined to other routes. Through the cytosolic tail of ICA512, proinsulin co-condensation with RESP18HD may further orchestrate the recruitment of cytosolic factors involved in membrane budding and fission of transport vesicles and nascent SGs.

Conserved residues in protein homolog sequence alignments are structurally or functionally important. For intrinsically disordered proteins or proteins with intrinsically disordered regions (IDRs), however, alignment often fails because they lack a steric structure to constrain evolution. Although sequences vary, the physicochemical features of IDRs may be preserved in maintaining function. Therefore, a method to retrieve common IDR features may help identify functionally important residues. We applied unsupervised contrastive learning to train a model with self-attention neuronal networks on human IDR orthologs. Parameters in the model were trained to match sequences in ortholog pairs but not in other IDRs. The trained model successfully identifies previously reported critical residues from experimental studies, especially those with an overall pattern (e.g., multiple aromatic residues or charged blocks) rather than short motifs. This predictive model can be used to identify potentially important residues in other proteins, improving our understanding of their functions. The trained model can be run directly from the Jupyter Notebook in the GitHub repository using Binder (mybinder.org). The only required input is the primary sequence. The training scripts are available on GitHub (https://github.com/allmwh/IFF). The training datasets have been deposited in an Open Science Framework repository (https://osf.io/jk29b).

The Escherichia coli cytochrome bo3 ubiquinol oxidase is a four-subunit heme-copper oxidase that serves as a proton pump in the E. coli aerobic respiratory chain. Despite many mechanistic studies, it is unclear whether this ubiquinol oxidase functions as a monomer, or as a dimer in a manner similar to its eukaryotic counterparts-the mitochondrial electron transport complexes. In this study, we determined the monomeric and dimeric structures of the E. coli cytochrome bo3 ubiquinol oxidase reconstituted in amphipol by cryogenic electron microscopy single particle reconstruction (cryo-EM SPR) to a resolution of 3.15 and 3.46 Å, respectively. We have discovered that the protein can form a dimer with C2 symmetry, with the dimerization interface maintained by interactions between the subunit II of one monomer and the subunit IV of the other monomer. Moreover, the dimerization does not induce significant structural changes in the monomers, except the movement of a loop in subunit IV (residues 67-74).

Many bacterial genomes encode a transmembrane protein kinase belonging to the PASTA kinase family, which controls numerous processes in diverse bacterial pathogens, including antibiotic resistance, cell division, stress resistance, toxin production, and virulence. PASTA kinases share a conserved three-part domain architecture, consisting of an extracellular PASTA domain, proposed to sense the peptidoglycan layer status, a single transmembrane helix, and an intracellular Ser/Thr kinase domain. The crystal structures of the kinase domain from two homologous PASTA kinases reveal a characteristic two-lobed structure typical of eukaryotic protein kinases with a centrally located, but unresolved, activation loop that becomes phosphorylated and regulates downstream signaling pathways. We previously identified three sites of phosphorylation on the activation loop (T163, T166, and T168) of IreK, a PASTA kinase from the pathogen Enterococcus faecalis, as well as a distal phosphorylation site (T218) that each influence IreK activity in vivo. Still, the mechanism by which loop phosphorylation regulates PASTA kinase function is yet unknown. Therefore, we utilized site-directed spin labeling (SDSL) and continuous wave (CW) electron paramagnetic resonance (EPR) spectroscopy to assess the E. faecalis IreK kinase activation loop dynamics, including the effects of phosphorylation on activation loop motion, and the IreK-IreB interaction. Our results reveal that the IreK activation loop occupies a more immobile state when dephosphorylated, and that loop autophosphorylation shifts the loop to a more mobile state that can then enable interaction with IreB, a known substrate.

Several hydrolases have been described to degrade polyethylene terephthalate (PET) at moderate temperatures ranging from 25°C to 40°C. These mesophilic PET hydrolases (PETases) are less efficient in degrading this plastic polymer than their thermophilic homologs and have, therefore, been the subject of many protein engineering campaigns. However, enhancing their enzymatic activity through rational design or directed evolution poses a formidable challenge due to the need for exploring a large number of mutations. Additionally, evaluating the improvements in both activity and stability requires screening numerous variants, either individually or using high-throughput screening methods. Here, we utilize instead the design of chimeras as a protein engineering strategy to increase the activity and stability of Mors1, an Antarctic PETase active at 25°C. First, we obtained the crystal structure of Mors1 at 1.6 Å resolution, which we used as a scaffold for structure- and sequence-based chimeric design. Then, we designed a Mors1 chimera via loop exchange of a highly divergent active site loop from the thermophilic leaf-branch compost cutinase (LCC) into the equivalent region in Mors1. After restitution of an active site disulfide bond into this chimera, the enzyme exhibited a shift in optimal temperature for activity to 45°C and an increase in fivefold in PET hydrolysis when compared with wild-type Mors1 at 25°C. Our results serve as a proof of concept of the utility of chimeric design to further improve the activity and stability of PETases active at moderate temperatures.

Phosphodiesterase-5 (PDE5) is responsible for regulating the concentration of the second messenger molecule cGMP by hydrolyzing it into 5'-GMP. PDE5 is implicated in erectile dysfunction and cardiovascular diseases. The substrate binding site in the catalytic domain of PDE5 is surrounded by several dynamic structural motifs (including the α14 helix, M-loop, and H-loop) that are known to switch between inactive and active conformational states via currently unresolved structural intermediates. We evaluated the conformational dynamics of these structural motifs in the apo state and upon binding of an allosteric inhibitor (evodiamine) or avanafil, a competitive inhibitor. We employed enhanced sampling-based replica exchange solute scaling (REST2) method, principal component analysis (PCA), time-lagged independent component analysis (tICA), molecular dynamics (MD) simulations, and well-tempered metadynamics simulations to probe the conformational changes in these structural motifs. Our results support a regulatory mechanism for PDE5, where the α14 helix alternates between an inward (lower activity) conformation and an outward (higher activity) conformation that is accompanied by the folding/unfolding of the α8' and α8″ helices of the H-loop. When the allosteric inhibitor evodiamine is bound to PDE5, the inward (inactive) state of the α14 helix is preferred, thus preventing substrate access to the catalytic site. In contrast, competitive inhibitors of PDE5 block catalysis by occupying the active site accompanied by stabilization of the outward conformation of the α14 helix. Defining the conformational dynamics underlying regulation of PDE5 activation will be helpful in rational design of next-generation small molecules modulators of PDE5 activity.

5-Aminolevulinic acid synthase (ALAS) is a pyridoxal 5'-phosphate (PLP)-dependent enzyme that catalyzes the first and rate-limiting step of heme biosynthesis in α-proteobacteria and several non-plant eukaryotes. All ALAS homologs contain a highly conserved catalytic core, but eukaryotes also have a unique C-terminal extension that plays a role in enzyme regulation. Several mutations in this region are implicated in multiple blood disorders in humans. In Saccharomyces cerevisiae ALAS (Hem1), the C-terminal extension wraps around the homodimer core to contact conserved ALAS motifs proximal to the opposite active site. To determine the importance of these Hem1 C-terminal interactions, we determined the crystal structure of S. cerevisiae Hem1 lacking the terminal 14 amino acids (Hem1 ΔCT). With truncation of the C-terminal extension, we show structurally and biochemically that multiple catalytic motifs become flexible, including an antiparallel β-sheet important to Fold-Type I PLP-dependent enzymes. The changes in protein conformation result in an altered cofactor microenvironment, decreased enzyme activity and catalytic efficiency, and ablation of subunit cooperativity. These findings suggest that the eukaryotic ALAS C-terminus has a homolog-specific role in mediating heme biosynthesis, indicating a mechanism for autoregulation that can be exploited to allosterically modulate heme biosynthesis in different organisms.

Proteases are involved in many physiologic processes, and dysregulated proteolysis is basis of a variety of diseases. Specific inhibition of pathogenetic proteases via monoclonal antibodies therefore holds significant therapeutic promise. Inspired by the competitive mechanism utilized by many naturally occurring and man-made protease inhibitors, we hypothesized that substrate-like peptide sequences can act as protease subsite blocking motifs if they occupy only one side of the reaction center. To test this hypothesis, a degenerate codon library representing MMP-14 substrate profiles at P1-P5' positions was constructed in the context of an anti-MMP-14 Fab by replacing its inhibitory motif in CDR-H3 with MMP-14 substrate repertoires. After selection for MMP-14 active-site binders by phage panning, results indicated that diverse substrate-like sequences conferring antibodies inhibitory potencies were enriched in the isolated clones. Optimal residues at each of P1-P5' positions were then identified, and the corresponding mutation combinations showed improved characteristics as effective inhibitors of MMP-14. Insights on efficient library designs for inhibitory peptide motifs were further discussed. Overall, this study proved the concept that substrate-derived sequences were able to behave as the inhibitory motifs in protease-specific antibodies. With accumulating data available on protease substrate profiles, we expect the approach described here can be broadly applied to facilitate the generation of antibody inhibitors targeting biomedically important proteases.

The efficiency of drug biosynthesis depends on different transcriptional regulatory pathways in Streptomyces, and the protein degradation system adds another layer of complexity to the regulatory processes. AtrA, a transcriptional regulator in the A-factor regulatory cascade, stimulates the production of daptomycin by binding to the dptE promoter in Streptomyces roseosporus. Using pull-down assays, bacterial two-hybrid system and knockout verification, we demonstrated that AtrA is a substrate for ClpP protease. Furthermore, we showed that ClpX is necessary for AtrA recognition and subsequent degradation. Bioinformatics analysis, truncating mutation, and overexpression proved that the AAA motifs of AtrA were essential for initial recognition in the degradation process. Finally, overexpression of mutated atrA (AAA-QQQ) in S. roseosporus increased the yield of daptomycin by 225% in shake flask and by 164% in the 15 L bioreactor. Thus, improving the stability of key regulators is an effective method to promote the ability of antibiotic synthesis.

Proteins gain optimal fitness such as foldability and function through evolutionary selection. However, classical studies have found that evolutionarily designed protein sequences alone cannot guarantee foldability, or at least not without considering local contacts associated with the initial folding steps. We previously showed that foldability and function can be restored by removing frustration in the folding energy landscape of a model WW domain protein, CC16, which was designed based on Statistical Coupling Analysis (SCA). Substitutions ensuring the formation of five local contacts identified as "on-path" were selected using the closest homolog native folded sequence, N21. Surprisingly, the resulting sequence, CC16-N21, bound to Group I peptides, while N21 did not. Here, we identified single-point mutations that enable N21 to bind a Group I peptide ligand through structure and dynamic-based computational design. Comparison of the docked position of the CC16-N21/ligand complex with the N21 structure showed that residues at positions 9 and 19 are important for peptide binding, whereas the dynamic profiles identified position 10 as allosterically coupled to the binding site and exhibiting different dynamics between N21 and CC16-N21. We found that swapping these positions in N21 with matched residues from CC16-N21 recovers nature-like binding affinity to N21. This study validates the use of dynamic profiles as guiding principles for affecting the binding affinity of small proteins.

(2R)-4-thiaproline (Thp) is an analog of proline, replacing Cγ in the pyrrolidine ring with sulfur. Its thiazolidine ring easily interconverts between endo and exo puckers due to a small energy barrier, which leads to destabilize polyproline helices. Collagen, composed of three polyproline II helices, mainly consists of X-Y-Gly triplets, where X is often proline and Y is frequently (2S,4R)-hydroxyproline. In this study, we incorporated Thp into either position-X or position-Y to investigate the consequences of such a replacement on the triple helix. Circular dichroism and differential scanning calorimetry analyses showed that the Thp-containing collagen-mimetic peptides (CMPs) can fold into stable triple helices, in which the substitution at position-Y exhibits a larger destabilization effect. Additionally, we also prepared the derivative peptides by oxidizing Thp in the peptide to N-formyl-cysteine or S,S-dioxide Thp. The results showed that the oxidized derivatives at position-X only slightly affect collagen stability, but those at position-Y induce a large destabilization effect. The consequences of incorporating Thp and its oxidized derivatives into CMPs are position dependent. Computational results suggested that the ease of interconversion between exo and endo puckers for Thp and the twist conformation of S,S-dioxide Thp may cause the destabilization effect at position-Y. We have revealed new insights into the impacts of Thp and its oxidized derivatives on collagen and demonstrated that Thp can be used to design collagen-related biomaterials.

Electronegative clusters (ENCs) made up of acidic residues and/or phosphorylation sites are the most abundant repetitive sequences in RNA-binding proteins. Previous studies have indicated that ENCs inhibit RNA binding for structured RNA-binding domains (RBDs). However, this is not the case for the unstructured RBD in histone pre-mRNA stem-loop binding protein (SLBP). The SLBP RBD contains 70 amino acids and is followed by a phosphorylatable ENC. ENC phosphorylation increases RNA-binding affinity of SLBP to the sub-picomolar range. In this study, we use NMR and molecular dynamics simulations to elucidate the mechanism for this tight binding. Our NMR data demonstrate that the ENC transiently folds apo SLBP into an RNA-bound resembling state. We find that in the RNA-bound state, the phosphorylated ENC interacts with the loop region opposite to the RNA-binding site. This allosteric interaction stabilizes the complex and therefore enhances RNA binding. To evaluate the generality of our findings, we graft an ENC onto endoribonuclease homolog 1's first double-stranded RNA-binding motif (DRBM1), an unstructured RBD that shares no homology with SLBP. We find that the engineered ENC increases the folded species of DRBM1 and inhibits RNA binding. On the contrary, introducing basic residues to DRBM1 makes the domain more unfolded, enhances RNA binding, and mitigates the inhibitory effect of the engineered ENC. In summary, our study suggests that ENCs promote folding of unstructured RNA-binding domains, and their effects on RNA binding depend on the electropositive charges on the RBD surface.

Amadori rearrangement products are stable sugar-amino acid conjugates that are formed nonenzymatically during preparation, dehydration, and storage of foods. Because Amadori compounds such as fructose-lysine (F-Lys), an abundant constituent in processed foods, shape the animal gut microbiome, it is important to understand bacterial utilization of these fructosamines. In bacteria, F-Lys is first phosphorylated, either during or after uptake to the cytoplasm, to form 6-phosphofructose-lysine (6-P-F-Lys). FrlB, a deglycase, then converts 6-P-F-Lys to L-lysine and glucose-6-phosphate. Here, to elucidate the catalytic mechanism of this deglycase, we first obtained a 1.8-Å crystal structure of Salmonella FrlB (without substrate) and then used computational approaches to dock 6-P-F-Lys on this structure. We also took advantage of the structural similarity between FrlB and the sugar isomerase domain of Escherichia coli glucosamine-6-phosphate synthase (GlmS), a related enzyme for which a structure with substrate has been determined. An overlay of FrlB-6-P-F-Lys on GlmS-fructose-6-phosphate structures revealed parallels in their active-site arrangement and guided our selection of seven putative active-site residues in FrlB for site-directed mutagenesis. Activity assays with eight recombinant single-substitution mutants identified residues postulated to serve as the general acid and general base in the FrlB active site and indicated unexpectedly significant contributions from their proximal residues. By exploiting native mass spectrometry (MS) coupled to surface-induced dissociation, we distinguished mutations that impaired substrate binding versus cleavage. As demonstrated with FrlB, an integrated approach involving x-ray crystallography, in silico approaches, biochemical assays, and native MS can synergistically aid structure-function and mechanistic studies of enzymes.

Angiotensin-converting enzyme 2 (ACE2) has been investigated for its ability to beneficially modulate the angiotensin receptor (ATR) therapeutic axis to treat multiple human diseases. Its broad substrate scope and diverse physiological roles, however, limit its potential as a therapeutic agent. In this work, we address this limitation by establishing a yeast display-based liquid chromatography screen that enabled use of directed evolution to discover ACE2 variants that possess both wild-type or greater Ang-II hydrolytic activity and improved specificity toward Ang-II relative to the off-target peptide substrate Apelin-13. To obtain these results, we screened ACE2 active site libraries to reveal three substitution-tolerant positions (M360, T371, and Y510) that can be mutated to enhance ACE2's activity profile and followed up on these hits with focused double mutant libraries to further improve the enzyme. Relative to wild-type ACE2, our top variant (T371L/Y510Ile) displayed a sevenfold increase in Ang-II turnover number (kcat ), a sixfold diminished catalytic efficiency (kcat /Km ) on Apelin-13, and an overall decreased activity on other ACE2 substrates that were not directly assayed in the directed evolution screen. At physiologically relevant substrate concentrations, T371L/Y510Ile hydrolyzes as much or more Ang-II than wild-type ACE2 with concomitant Ang-II:Apelin-13 specificity improvements reaching 30-fold. Our efforts have delivered ATR axis-acting therapeutic candidates with relevance to both established and unexplored ACE2 therapeutic applications and provide a foundation for further ACE2 engineering efforts.

Insects have evolved a chemical communication system using terpenoids, a structurally diverse class of specialized metabolites, previously thought to be exclusively produced by plants and microbes. Gene discovery, bioinformatics, and biochemical characterization of multiple insect terpene synthases (TPSs) revealed that isopentenyl diphosphate synthases (IDS), enzymes from primary isoprenoid metabolism, are their likely evolutionary progenitors. However, the mutations underlying the emergence of the TPS function remain a mystery. To address this gap, we present the first structural and mechanistic model for the evolutionary emergence of TPS function in insects. Through identifying key mechanistic differences between IDS and TPS enzymes, we hypothesize that the loss of isopentenyl diphosphate (IPP) binding motifs strongly correlates with the gain of the TPS function. Based on this premise, we have elaborated the first explicit structural definition of isopentenyl diphosphate-binding motifs (IBMs) and used the IBM definitions to examine previously characterized insect IDSs and TPSs and to predict the functions of as yet uncharacterized insect IDSs. Consistent with our hypothesis, we observed a clear pattern of disruptive substitutions to IBMs in characterized insect TPSs. In contrast, insect IDSs maintain essential consensus residues for binding IPP. Extending our analysis, we constructed the most comprehensive phylogeny of insect IDS sequences (430 full length sequences from eight insect orders) and used IBMs to predict the function of TPSs. Based on our analysis, we infer multiple, independent TPS emergence events across the class of insects, paving the way for future gene discovery efforts.

Amino acids (AAs) with a noncanonical backbone would be a valuable tool for protein engineering, enabling new structural motifs and building blocks. To incorporate them into an expanded genetic code, the first, key step is to obtain an appropriate aminoacyl-tRNA synthetase. Currently, directed evolution is not available to optimize AAs with noncanonical backbones, since an appropriate selective pressure has not been discovered. Computational protein design (CPD) is an alternative. We used a new CPD method to redesign MetRS and increase its activity towards β-Met, which has an extra backbone methylene. The new method considered a few active site positions for design and used a Monte Carlo exploration of the corresponding sequence space. During the exploration, a bias energy was adaptively learned, such that the free energy landscape of the apo enzyme was flattened. Enzyme variants could then be sampled, in the presence of the ligand and the bias energy, according to their β-Met binding affinities. Eighteen predicted variants were chosen for experimental testing; 10 exhibited detectable activity for β-Met adenylation. Top predicted hits were characterized experimentally in detail. Dissociation constants, catalytic rates, and Michaelis constants for both α-Met and β-Met were measured. The best mutant retained a preference for α-Met over β-Met; however, the preference was reduced, compared to the wildtype, by a factor of 29. For this mutant, high resolution crystal structures were obtained in complex with both α-Met and β-Met, indicating that the predicted, active conformation of β-Met in the active site was retained.

Inosine 5'-monophosphate (IMP) dehydrogenase (IMPDH) is an ubiquitous enzyme that catalyzes the NAD+ -dependent oxidation of inosine 5'-monophosphate into xanthosine 5'-monophosphate. This enzyme is formed of two distinct domains, a core domain where the catalytic reaction occurs, and a less-conserved Bateman domain. Our previous studies gave rise to the classification of bacterial IMPDHs into two classes, according to their oligomeric and kinetic properties. MgATP is a common effector but cause to different effects when it binds within the Bateman domain: it is either an allosteric activator for Class I IMPDHs or a modulator of the oligomeric state for Class II IMPDHs. To get insight into the role of the Bateman domain in the dissimilar properties of the two classes, deleted variants of the Bateman domain and chimeras issued from the interchange of the Bateman domain between the three selected IMPDHs have been generated and characterized using an integrative structural biology approach. Biochemical, biophysical, structural, and physiological studies of these variants unveil the Bateman domain as being the carrier of the molecular behaviors of both classes.

Monotopic phosphoglycosyl transferases (monoPGTs) are an expansive superfamily of enzymes that catalyze the first membrane-committed step in the biosynthesis of bacterial glycoconjugates. MonoPGTs show a strong preference for their cognate nucleotide diphospho-sugar (NDP-sugar) substrates. However, despite extensive characterization of the monoPGT superfamily through previous development of a sequence similarity network comprising >38,000 nonredundant sequences, the connection between monoPGT sequence and NDP-sugar substrate specificity has remained elusive. In this work, we structurally characterize the C-terminus of a prototypic monoPGT for the first time and show that 19 C-terminal residues play a significant structural role in a subset of monoPGTs. This new structural information facilitated the identification of co-conserved sequence "fingerprints" that predict NDP-sugar substrate specificity for this subset of monoPGTs. A Hidden Markov model was generated that correctly assigned the substrate of previously unannotated monoPGTs. Together, these structural, sequence, and biochemical analyses have delivered new insight into the determinants guiding substrate specificity of monoPGTs and have provided a strategy for assigning the NDP-sugar substrate of a subset of enzymes in the superfamily that use UDP-di-N-acetyl bacillosamine. Moving forward, this approach may be applied to identify additional sequence motifs that serve as fingerprints for monoPGTs of differing UDP-sugar substrate specificity.

Polyketide metabolism-associated proteins in Mycobacterium tuberculosis play an essential role in the survival of the bacterium, which makes them potential drug targets for the treatment of tuberculosis (TB). The novel ribonuclease protein Rv1546 is predicted to be a member of the steroidogenic acute regulatory protein-related lipid-transfer (START) domain superfamily, which comprises bacterial polyketide aromatase/cyclases (ARO/CYCs). Here, we determined the crystal structure of Rv1546 in a V-shaped dimer. The Rv1546 monomer consists of four α-helices and seven antiparallel β-strands. Interestingly, in the dimeric state, Rv1546 forms a helix-grip fold, which is present in START domain proteins, via three-dimensional domain swapping. Structural analysis revealed that the conformational change of the C-terminal α-helix of Rv1546 might contribute to the unique dimer structure. Site-directed mutagenesis followed by in vitro ribonuclease activity assays was performed to identify catalytic sites of the protein. This experiment suggested that surface residues R63, K84, K88, and R113 are important in the ribonuclease function of Rv1546. In summary, this study presents the structural and functional characterization of Rv1546 and supplies new perspectives for exploiting Rv1546 as a novel drug target for TB treatment.

Orthogonal translation is an efficient tool that provides many valuable spectral probes capable of covering different parts of the electromagnetic spectrum and thus enabling parameterization of various structural and dynamic phenomena in proteins. In this context, nitrile-containing tryptophan analogs are very useful probes to study local electrostatics and hydrogen bonding in both rigid and dynamic environments. Here, we report a semi-rational approach to engineer a tyrosyl-tRNA synthetase (TyrRS) variant of Methanocaldococcus jannaschii capable of incorporating 5-cyanotryptophan (5CNW) via orthogonal translation. We combined one round of the well-established positive selection system with saturation mutagenesis at preselected TyrRS positions, resulting in a novel 5CNW-specific enzyme that also exhibits high substrate tolerance to other aromatic noncanonical amino acids. We demonstrated the utility of our orthogonal pair by inserting 5CNW into the cyanobacteriochrome Slr1393g3, a bilin-binding photosensor of the phytochrome superfamily. The nitrile (CN) group of the inserted 5CNW provides non-invasive labeling in the local structural context while yielding information on local electrostatics and hydrogen bonding by IR spectroscopy. 5CNW is a versatile probe that can be used for both static and dynamic measurements.

Hsp90 is a molecular chaperone involved in the refolding and activation of numerous protein substrates referred to as clients. While the molecular determinants of Hsp90 client specificity are poorly understood and limited to a handful of client proteins, strong clients are thought to be destabilized and conformationally extended. Here, we measured the phosphotransferase activity of 3929 variants of the tyrosine kinase Src in both the presence and absence of an Hsp90 inhibitor. We identified 84 previously unknown functionally dependent client variants. Unexpectedly, many destabilized or extended variants were not functionally dependent on Hsp90. Instead, functionally dependent client variants were clustered in the αF pocket and β1-β2 strand regions of Src, which have yet to be described in driving Hsp90 dependence. Hsp90 dependence was also strongly correlated with kinase activity. We found that a combination of activation, global extension, and general conformational flexibility, primarily induced by variants at the αF pocket and β1-β2 strands, was necessary to render Src functionally dependent on Hsp90. Moreover, the degree of activation and flexibility required to transform Src into a functionally dependent client varied with variant location, suggesting that a combination of regulatory domain disengagement and catalytic domain flexibility are required for chaperone dependence. Thus, by studying the chaperone dependence of a massive number of variants, we highlight factors driving Hsp90 client specificity and propose a model of chaperone-kinase interactions.

The introduction of an engineered aminoacyl-tRNA synthetase/tRNA pair enables site-specific incorporation of unnatural amino acids (uAAs) with functionalized side chains into proteins of interest. Genetic Code Expansion (GCE) via amber codon suppression confers functionalities to proteins but can also be used to temporally control the incorporation of genetically encoded elements into proteins. Here, we report an optimized GCE system (GCEXpress) for efficient and fast uAA incorporation. We demonstrate that GCEXpress can be used to efficiently alter the subcellular localization of proteins within living cells. We show that click labeling can resolve co-labeling problems of intercellular adhesive protein complexes. We apply this strategy to study the adhesion G protein-coupled receptor (aGPCR) ADGRE5/CD97 and its ligand CD55/DAF that play central roles in immune functions and oncological processes. Furthermore, we use GCEXpress to analyze the time course of ADGRE5-CD55 ligation and replenishment of mature receptor-ligand complexes. Supported by fluorescence recovery after photobleaching (FRAP) experiments our results show that ADGRE5 and CD55 form stable intercellular contacts that may support transmission of mechanical forces onto ADGRE5 in a ligand-dependent manner. We conclude that GCE in combination with biophysical measurements can be a useful approach to analyze the adhesive, mechanical and signaling properties of aGPCRs and their ligand interactions.

The mTORC1 signaling pathway regulates cell growth and metabolism in a variety of organisms from yeast to human, and inhibition of the mTORC1 pathway has the prospect to treat cancer or achieve longevity. The tuberous sclerosis protein complex (TSCC) is a master negative regulator of the mTORC1 signaling pathway through hydrolyzing the GTP loaded on the small GTPase Rheb, which is a key activator of mTOR. However, the large size (~700 kDa) and complex structural organization of TSCC render it vulnerable to degradation and inactivation, thus limiting its potential application. In this work, based on thorough analysis and understanding of the structural mechanism of how the stabilization domain of TSC2 secures the association of TSC2-GAP with Rheb and thus enhances its GAP activity, we designed two proteins, namely SSG-MTM (short stabilization domain and GAP domain-membrane targeting motif) and SSG-TSC1N, which were able to function like TSCC to negatively regulate Rheb and mTORC1, but with much-reduced sizes (~1/15 and ~ 1/9 of the size of TSCC, respectively). Biochemical and cell biological assays demonstrated that these designed proteins indeed could promote the GTPase activity of Rheb to hydrolyze GTP, inhibit the kinase activity of mTORC1, and prevent mTORC1 from down-regulating catabolism and autophagy.

The type III secretion system (T3SS) is a large, transmembrane protein machinery used by various pathogenic gram-negative bacteria to transport virulence factors into the host cell during infection. Understanding the structure of T3SSs is crucial for future developments of therapeutics that could target this system. However, much of the knowledge about the structure of T3SS is available only for Salmonella, and it is unclear how this large assembly is conserved across species. Here, we combined cryo-electron microscopy, cross-linking mass spectrometry, and integrative modeling to determine the structure of the T3SS needle complex from Shigella flexneri. We show that the Shigella T3SS exhibits unique features distinguishing it from other structurally characterized T3SSs. The secretin pore complex adopts a new fold of its C-terminal S domain and the pilotin MxiM[SctG] locates around the outer surface of the pore. The export apparatus structure exhibits a conserved pseudohelical arrangement but includes the N-terminal domain of the SpaS[SctU] subunit, which was not present in any of the previously published virulence-related T3SS structures. Similar to other T3SSs, however, the apparatus is anchored within the needle complex by a network of flexible linkers that either adjust conformation to connect to equivalent patches on the secretin oligomer or bind distinct surface patches at the same height of the export apparatus. The conserved and unique features delineated by our analysis highlight the necessity to analyze T3SS in a species-specific manner, in order to fully understand the underlying molecular mechanisms of these systems. The structure of the type III secretion system from Shigella flexneri delineates conserved and unique features, which could be used for the development of broad-range therapeutics.

Isoforms of heat shock protein 90 (HSP90) fold oncoproteins that facilitate all 10 hallmarks of cancer. However, its promise as a therapeutic target remains unfulfilled as there is still no FDA-approved drug targeting HSP90 in disease. Among the reasons hindering progress are side effects caused by pan-HSP90 inhibition. Selective targeting of the four isoforms is challenging due to high sequence and structural similarity. Surprisingly, while decades of drug discovery efforts have produced almost 400 human HSP90 structures, no single ligand has been structurally characterized across all four human isoforms to date, which could reveal structural differences to achieve selectivity. To better understand the HSP90 landscape relevant for ligand binding and design we take a three-pronged approach. First, we solved the first complete set of structures of a single ligand bound to all four human isoforms. This enabled a systematic comparison of how side-chains and water networks respond to ligand binding across isoforms. Second, we expanded our analysis to publicly available, incomplete isoform-ligand series with distinct ligand chemistry. This highlighted general trends of protein and water mobility that differ among isoforms and impact ligand binding. Third, we further probed the Hsp90α conformational landscape for accommodating a congeneric series containing the purine scaffold common to HSP90 inhibitors. This revealed how minor ligand modifications flip ligand poses and perturb water and protein conformations. Taken together, this work illustrates how a systematic approach can shed new light on an "old" target and reveal hidden isoform-specific accommodations of congeneric ligands that may be exploited in ligand discovery and design.

Recombinant spider silk proteins can be prepared in scalable fermentation processes and have been proven as sources of biomaterials for biomedical and technical applications. Nanofibrils, formed through the self-assembly of these proteins, possess unique structural and mechanical properties, serving as fundamental building blocks for the fabrication of micro- and nanostructured scaffolds. Despite significant progress in utilizing nanofibrils-based morphologies of recombinant spider silk proteins, a comprehensive understanding of the molecular mechanisms of nanofibrils self-assembly remains a challenge. Here, a detailed kinetic study of nanofibril formation from a recombinant spider silk protein eADF4(C16) in dependence on the protein concentration, seeding, and temperature is provided. For the global fitting of kinetic data obtained during the fibril formation, we utilized the online platform AmyloFit. Evaluation of the data revealed that the self-assembly mechanism of recombinant spider silk is dominated by secondary nucleation. Thermodynamic analyses show that both primary and secondary nucleations, as well as the elongation step of the eADF4(C16), are endothermic processes.

Human pre-mRNA processing protein 40 homolog A (hPrp40A) is a splicing factor that interacts with the Huntington's disease protein huntingtin (Htt). Evidence has accumulated that both Htt and hPrp40A are modulated by the intracellular Ca2+ sensor calmodulin (CaM). Here we report characterization of the interaction of human CM with the third FF domain (FF3 ) of hPrp40A using calorimetric, fluorescence and structural approaches. Homology modeling, differential scanning calorimetry and small angle X-ray scattering (SAXS) data show FF3 forms a folded globular domain. CaM was found to bind FF3 in a Ca2+ -dependent manner with a 1:1 stoichiometry and a dissociation constant (Kd ) of 25 ± 3 μM at 25°C. NMR studies showed that both domains of CaM are engaged in binding and SAXS analysis of the FF3 -CaM complex revealed CaM occupies an extended configuration. Analysis of the FF3 sequence showed that the anchors for CaM binding must be buried in its hydrophobic core, suggesting that binding to CaM requires unfolding of FF3 . Trp anchors were proposed based on sequence analysis and confirmed by intrinsic Trp fluorescence of FF3 upon binding of CaM and substantial reductions in affinity for Trp-Ala FF3 mutants. The consensus model of the complex showed that binding to CaM binding occurs to an extended, non-globular state of the FF3 , consistent with coupling to transient unfolding of the domain. The implications of these results are discussed in the context of the complex interplay of Ca2+ signaling and Ca2+ sensor proteins in modulating Prp40A-Htt function.

Transthyretin (TTR)-related amyloidosis (ATTR) is a syndrome of diseases characterized by the extracellular deposition of fibrillar materials containing TTR variants. Ala97Ser (A97S) is the major mutation reported in Taiwanese ATTR patients. Here, we combine atomic resolution structural information together with the biochemical data to demonstrate that substitution of polar Ser for a small hydrophobic side chain of Ala at residue 97 of TTR largely influences the local packing density of the FG-loop, thus leading to the conformational instability of native tetramer, the increased monomeric species, and thus the enhanced amyloidogenicity of apo-A97S. Based on calorimetric studies, the tetramer destabilization of A97S can be substantially altered by interacting with native stabilizers via similarly energetic patterns compared to that of wild-type (WT) TTR; however, stabilizer binding partially rearranges the networks of hydrogen bonding in TTR variants while FG-loops of tetrameric A97S still remain relatively flexible. Moreover, TTR in complexed with holo-retinol binding protein 4 is slightly influenced by the structural and dynamic changes of FG-loop caused by A97S substitution with an approximately five-fold difference in binding affinity. Collectively, our findings suggest that the amyloidogenic A97S mutation destabilizes TTR by increasing the flexibility of the FG-loop in the monomer, thus modulating the rate of amyloid fibrillization.

As a consequence of the fast resistance spreading, a limited number of drugs are available to treat fungal infections. Therefore, there is an urgent need to develop new antifungal treatment strategies. The features of a disulfide bond-stabilized antifungal protein, NFAP2 secreted by the mold Neosartorya (Aspergillus) fischeri render it to be a promising template for future protein-based antifungal drug design, which requires knowledge about the native disulfide linkage pattern as it is one of the prerequisites for biological activity. However, in the lack of tryptic and chymotryptic proteolytic sites in the ACNCPNNCK sequence, the determination of the disulfide linkage pattern of NFAP2 is not easy with traditional mass spectrometry-based methods. According to in silico predictions working with a preliminary nuclear magnetic resonance (NMR) solution structure, two disulfide isomers of NFAP2 (abbacc and abbcac) were possible. Both were chemically synthesized; and comparative reversed-phase high-performance liquid chromatography, electronic circular dichroism and NMR spectroscopy analyses, and antifungal susceptibility and efficacy tests indicated that the abbcac is the native pattern. This knowledge allowed rational modification of NAFP2 to improve the antifungal efficacy and spectrum through the modulation of the evolutionarily conserved γ-core region, which is responsible for the activity of several antimicrobial peptides. Disruption of the steric structure of NFAP2 upon γ-core modification led to the conclusions that this motif may affect the formation of the biologically active three-dimensional structure, and that the γ-core modulation is not an efficient tool to improve the antifungal efficacy or to change the antifungal spectrum of NFAP2.

Glucosamine-6-phosphate (GlcN6P) deaminases from Escherichia coli (EcNagBI) and Shewanella denitrificans (SdNagBII) are special examples of what constitute nonhomologous isofunctional enzymes due to their convergence, not only in catalysis, but also in cooperativity and allosteric properties. Additionally, we found that the sigmoidal kinetics of SdNagBII cannot be explained by the existing models of homotropic activation. This study describes the regulatory mechanism of SdNagBII using enzyme kinetics, isothermal titration calorimetry (ITC), and X-ray crystallography. ITC experiments revealed two different binding sites with distinctive thermodynamic signatures: a single binding site per monomer for the allosteric activator N-acetylglucosamine 6-phosphate (GlcNAc6P) and two binding sites per monomer for the transition-state analog 2-amino-2-deoxy-D-glucitol 6-phosphate (GlcNol6P). Crystallographic data demonstrated the existence of an unusual allosteric site that can bind both GlcNAc6P and GlcNol6P, implying that the homotropic activation of this enzyme arises from the occupation of the allosteric site by the substrate. In this work we describe the presence of this novel allosteric site in the SIS-fold deaminases, which is responsible for the homotropic and heterotropic activation of SdNagBII by GlcN6P and GlcNAc6P, respectively. This study unveils an original mechanism to generate a high degree of homotropic activation in SdNagBII, mimicking the allosteric and cooperative properties of hexameric EcNagBI but with a reduced number of subunits.

Ym1 (chitinase-like 3, Chil3) expressed in mice is a nonenzymatic chitinase-like protein, which shows 67% identity with mouse acidic chitinase (Chia). Similar to Chia, Ym1 is overexpressed in asthma and parasitic infections in mouse lungs. Due to the lack of chitin-degrading activity, the biomedical role of Ym1 under these pathophysiological conditions remains to be determined. In this study, we investigated what region and amino acid changes in Ym1 resulted in the loss of enzymatic activity. Replacing two amino acids at the catalytic motif to obtain a Chia-like sequence (N136D and Q140E; MT-Ym1) did not activate the protein. We conducted a comparative study of Ym1 and Chia. We found that three protein segments-(i) the catalytic motif residues, (ii) exons 6 and 7, and (iii) exon 10-are responsible for chitinase activity loss in Ym1. We show that replacing each of these three segments in Chia that are also involved in substrate recognition and binding by the Ym1 sequence can fully abolish the enzymatic activity. In addition, we show that there have been extensive gene duplication events at the Ym1 locus specific to the rodent lineages. Consistent with this result, Ym1 orthologs from the rodent genome were under positive selection when analyzed through the CODEML program. These data suggest that numerous amino acid substitutions in the regions involved in the chitin recognition, binding, and degradation ability of the ancestor Ym1 molecule lead to the irreversible inactivation of the protein.

Cystathionine γ-lyase (CGL) is a PLP-dependent enzyme that catalyzes the last step of the reverse transsulfuration route for endogenous cysteine biosynthesis. The canonical CGL-catalyzed process consists of an α,γ-elimination reaction that breaks down cystathionine into cysteine, α-ketobutyrate, and ammonia. In some species, the enzyme can alternatively use cysteine as a substrate, resulting in the production of hydrogen sulfide (H2 S). Importantly, inhibition of the enzyme and consequently of its H2 S production activity, makes multiresistant bacteria considerably more susceptible to antibiotics. Other organisms, such as Toxoplasma gondii, the causative agent of toxoplasmosis, encode a CGL enzyme (TgCGL) that almost exclusively catalyzes the canonical process, with only minor reactivity to cysteine. Interestingly, the substitution of N360 by a serine (the equivalent amino acid residue in the human enzyme) at the active site changes the specificity of TgCGL for the catalysis of cystathionine, resulting in an enzyme that can cleave both the CγS and the CβS bond of cystathionine. Based on these findings and to deepen the molecular basis underlying the enzyme-substrate specificity, we have elucidated the crystal structures of native TgCGL and the variant TgCGL-N360S from crystals grown in the presence of cystathionine, cysteine, and the inhibitor d,l-propargylglycine (PPG). Our structures reveal the binding mode of each molecule within the catalytic cavity and help explain the inhibitory behavior of cysteine and PPG. A specific inhibitory mechanism of TgCGL by PPG is proposed.

We investigated the relationship between mutations and dynamics in Escherichia coli dihydrofolate reductase (DHFR) using computational methods. Our study focused on the M20 and FG loops, which are known to be functionally important and affected by mutations distal to the loops. We used molecular dynamics simulations and developed position-specific metrics, including the dynamic flexibility index (DFI) and dynamic coupling index (DCI), to analyze the dynamics of wild-type DHFR and compared our results with existing deep mutational scanning data. Our analysis showed a statistically significant association between DFI and mutational tolerance of the DHFR positions, indicating that DFI can predict functionally beneficial or detrimental substitutions. We also applied an asymmetric version of our DCI metric (DCIasym ) to DHFR and found that certain distal residues control the dynamics of the M20 and FG loops, whereas others are controlled by them. Residues that are suggested to control the M20 and FG loops by our DCIasym metric are evolutionarily nonconserved; mutations at these sites can enhance enzyme activity. On the other hand, residues controlled by the loops are mostly deleterious to function when mutated and are also evolutionary conserved. Our results suggest that dynamics-based metrics can identify residues that explain the relationship between mutation and protein function or can be targeted to rationally engineer enzymes with enhanced activity.

Förster resonance energy transfer (FRET) is a valuable method for monitoring protein conformation and biomolecular interactions. Intrinsically fluorescent amino acids that can be genetically encoded, such as acridonylalanine (Acd), are particularly useful for FRET studies. However, quantitative interpretation of FRET data to derive distance information requires careful use of controls and consideration of photophysical effects. Here we present two case studies illustrating how Acd can be used in FRET experiments to study small molecule induced conformational changes and multicomponent biomolecular complexes.

The Escherichia coli enzyme EcAIII catalyzes the hydrolysis of L-Asn to L-Asp and ammonia. Using a nature-inspired mutagenesis approach, we designed and produced five new EcAIII variants (M200I, M200L, M200K, M200T, M200W). The modified proteins were characterized by spectroscopic and crystallographic methods. All new variants were enzymatically active, confirming that the applied mutagenesis procedure has been successful. The determined crystal structures revealed new conformational states of the EcAIII molecule carrying the M200W mutation and allowed a high-resolution observation of an acyl-enzyme intermediate with the M200L mutant. In addition, we performed structure prediction, substrate docking, and molecular dynamics simulations for 25 selected bacterial orthologs of EcAIII, to gain insights into how mutations at the M200 residue affect the active site and substrate binding mode. This comprehensive strategy, including both experimental and computational methods, can be used to guide further enzyme engineering and can be applied to the study of other proteins of medicinal or biotechnological importance.

[4Fe-4S]2+ cluster assembly in human cytosol requires both a [2Fe-2S] cluster chaperone being able to donate two [2Fe-2S]2+ clusters and an electron donor providing two electrons to reductively couple the two [2Fe-2S]2+ clusters into a [4Fe-4S]2+ cluster. The mechanism through which the cytosolic [4Fe-4S]2+ cluster assembly works is still not defined. Here, we show that a hetero-tetrameric complex formed by two molecules of cluster-reduced [2Fe-2S]+ 2 -anamorsin and one molecule of dimeric cluster-oxidized [2Fe-2S]2+ 2 -GLRX32 orchestrates the assembly of a [4Fe-4S]2+ cluster on the N-terminal cluster binding site of the cytosolic protein NUBP1. We demonstrate that the hetero-tetrameric complex is able to synergically provide two [2Fe-2S]2+ clusters from GLRX3 and two electrons from anamorsin for the assembly of the [4Fe-4S]2+ cluster on the N-terminal cluster binding site of NUBP1. We also showed that only one of the two [2Fe-2S] clusters bound to anamorsin, that is, that bound to the CX8 CX2 CXC motif, provides the electrons required to form the [4Fe-4S]2+ cluster. Our study contributes to the molecular understanding of the mechanism of [4Fe-4S] protein biogenesis in the cytosol.

Type Three Secretion System (T3SS) is a sophisticated nano-scale weapon utilized by several gram negative bacteria under stringent spatio-temporal regulation to manipulate and evade host immune systems in order to cause infection. To the best of our knowledge, this present study is the first report where we embark upon characterizing inherent features of native type three secretion effector protein PemB through biophysical techniques. Herein, first, we demonstrate binding affinity of PemB for phosphoinositides through isothermal calorimetric titrations. Second, we shed light on its strong homo-oligomerization propensity in aqueous solution through multiple biophysical methods. Third, we also employ several spectroscopic techniques to delineate its disordered and helical conformation. Lastly, we perform a phylogenetic analysis of this new effector to elucidate evolutionary relationship with other organisms. Taken together, our results shall surely contribute to our existing knowledge of Pseudomonas aeruginosa secretome.

The glycerophosphodiester phosphodiesterase (GDPD)-like SMaseD/PLD domain family, which includes phospholipase D (PLD) toxins in recluse spiders and actinobacteria, evolved anciently in bacteria from the GDPD. The PLD enzymes retained the core (β/α)8 barrel fold of GDPD, while gaining a signature C-terminal expansion motif and losing a small insertion domain. Using sequence alignments and phylogenetic analysis, we infer that the C-terminal motif derives from a segment of an ancient bacterial PLAT domain. Formally, part of a protein containing a PLAT domain repeat underwent fusion to the C terminus of a GDPD barrel, leading to attachment of a segment of a PLAT domain, followed by a second complete PLAT domain. The complete domain was retained only in some basal homologs, but the PLAT segment was conserved and repurposed as the expansion motif. The PLAT segment corresponds to strands β7-β8 of a β-sandwich, while the expansion motif as represented in spider PLD toxins has been remodeled as an α-helix, a β-strand, and an ordered loop. The GDPD-PLAT fusion led to two acquisitions in founding the GDPD-like SMaseD/PLD family: (1) a PLAT domain that presumably supported early lipase activity by mediating membrane association, and (2) an expansion motif that putatively stabilized the catalytic domain, possibly compensating for, or permitting, loss of the insertion domain. Of wider significance, messy domain shuffling events can leave behind scraps of domains that can be salvaged, remodeled, and repurposed.

Virtual screening (VS) is a routine method to evaluate chemical libraries for lead identification. Therefore, the selection of appropriate protein structures for VS is an essential prerequisite to identify true actives during docking. But the presence of several crystal structures of the same protein makes it difficult to select one or few structures rationally for screening. Therefore, a computational prioritization protocol has been developed for shortlisting crystal structures that identify true active molecules with better efficiency. As identification of small-molecule inhibitors is an important clinical requirement for the T790M/L858R (TMLR) EGFR mutant, it has been selected as a case study. The approach involves cross-docking of 21 co-crystal ligands with all the structures of the same protein to select structures that dock non-native ligands with lower RMSD. The cross docking performance was then correlated with ligand similarity and binding-site conformational similarity. Eventually, structures were shortlisted by integrating cross-docking performance, and ligand and binding-site similarity. Thereafter, binding pose metadynamics was employed to identify structures having stable co-crystal ligands in their respective binding pockets. Finally, different enrichment metrics like BEDROC, RIE, AUAC, and EF1% were evaluated leading to the identification of five TMLR structures (5HCX, 5CAN, 5CAP, 5CAS, and 5CAO). These structures docked a number of non-native ligands with low RMSD, contain structurally dissimilar ligands, have conformationally dissimilar binding sites, harbor stable co-crystal ligands, and also identify true actives early. The present approach can be implemented for shortlisting protein targets of any other important therapeutic kinases.

The HspB8-BAG3 complex plays an important role in the protein quality control acting alone or within multi-components complexes. To clarify the mechanism underlying its activity, in this work we used biochemical and biophysical approaches to study the tendency of both proteins to auto-assemble and to form the complex. Solubility and Thioflavin T assays, Fourier transform infrared spectroscopy and atomic force microscopy analyses clearly showed the tendency of HspB8 to self-assemble at high concentration and to form oligomers in a "native-like" conformation; otherwise, BAG3 aggregates poorly. Noteworthy, also HspB8 and BAG3 associate in a "native-like" conformation, forming a stable complex. Furthermore, the high difference between dissociation constant values of HspB8-HspB8 interaction with respect to the binding to BAG3 obtained by surface plasmon resonance confirms that HspB8 is an obligated partner of BAG3 in vivo. Lastly, both proteins alone or in the complex are able to bind and affect the aggregation of the Josephin domain, the structured domain that triggers the ataxin-3 fibrillation. In particular, the complex displayed higher activity than HspB8 alone. All this considered, we can assert that the two proteins form a stable assembly with chaperone-like activity that could contribute to the physiological role of the complex in vivo.

The irreversible aggregation of proteins or peptides greatly limits their bioavailability; therefore, effective inhibition using small molecules or biocompatible materials is very difficult. Human calcitonin (hCT), a hormone polypeptide with 32 residues, is secreted by the C-cells of the thyroid gland. The biological function of this hormone is to regulate calcium and phosphate concentrations in the blood via several different pathways. One of these is to inhibit the activity of osteoclasts; thus, calcitonin could be used to treat osteoporosis and Paget's disease of the bone. However, hCT is prone to aggregation in aqueous solution and forms amyloid fibrils. Salmon and eel calcitonin are currently used as clinical substitutes for hCT. In a previous study, we found that the replacement of two residues at positions 12 and 17 of hCT with amino acids that appear in the salmon sequence can greatly suppress peptide aggregation. The double mutations of hCT (DM hCT) also act as good inhibitors by disrupting wild-type hCT fibrillization, although the inhibition mechanism is not clear. More importantly, we demonstrated that DM hCT is biologically active in interacting with the calcitonin receptor. To further understand the inhibitory effect of DM hCT on hCT fibrillization, we created four relevant peptide fragments based on the DM hCT sequence. Our examination revealed that the formation of a helix of DM hCT was possibly a key component contributing to its inhibitory effect. This finding could help in the development of peptide-based inhibitors and in understanding the aggregation mechanism of hCT.

A sensor of blue-light using flavin adenine dinucleotide (BLUF) is a typical blue light photoreceptor domain that is found in many photosensor proteins in bacteria and some eukaryotic algae. SyPixD in Synechocystis is one of the well-studied BLUF proteins. In the dark state, it forms a decamer and, upon photoexcitation, a dissociation reaction takes place to yield dimers. Such change in the intermolecular interactions of the protomers is important for the biological function. The effect of the N- and C-terminal sequences on the stability of SyPixD oligomeric states and photoreactions of SyPixD were studied to understand how the oligomeric form is maintained with weak interaction. It was found that a few residues that frequently persist at the N-terminus after removing a tag for purification are sensitive to the stability of the decamer structure. Even two or three residues at the N-terminus considerably reduces decamer stability, whereas four or more residues completely prevent decamer formation. Unexpectedly, truncating C-terminal sequences, which locate far from any protomer interface and of which structure is undetermined in crystal structure, also destabilizes the decamer structure. This destabilization is also apparent from the dissociation reaction dynamics detected by the transient grating and transient absorption measurements. The dissociation reaction is faster and the yield increases when the C-terminus does not contain seven amino acid residues. Photoexcitation induces a conformational change in the C-terminus of the decamer but not the dimer.

Cell signaling by receptor protein tyrosine kinases (RTKs) is tightly controlled by the counterbalancing actions of receptor protein tyrosine phosphatases (RPTPs). Due to their role in attenuating the signal-initiating potency of RTKs, RPTPs have long been viewed as therapeutic targets. However, the development of activators of RPTPs has remained limited. We previously reported that the homodimerization of a representative member of the RPTP family (protein tyrosine phosphatase receptor J or PTPRJ) is regulated by specific transmembrane (TM) residues. Disrupting this interaction by single point mutations promotes PTPRJ access to its RTK substrates (e.g., EGFR and FLT3), reduces RTK's phosphorylation and downstream signaling, and ultimately antagonizes RTK-driven cell phenotypes. Here, we designed and tested a series of first-in-class pH-responsive TM peptide agonists of PTPRJ that are soluble in aqueous solution but insert as a helical TM domain in lipid membranes when the pH is lowered to match that of the acidic microenvironment of tumors. The most promising peptide reduced EGFR's phosphorylation and inhibited cancer cell EGFR-driven migration and proliferation, similar to the PTPRJ's TM point mutations. Developing tumor-selective and TM-targeting peptide binders of critical RPTPs could afford a potentially transformative approach to studying RPTP's selectivity mechanism without requiring less specific inhibitors and represent a novel class of therapeutics against RTK-driven cancers.

Intrinsically disordered proteins (IDPs) are often multifunctional and frequently posttranslationally modified. Deleted in split hand/split foot 1 (Dss1-Sem1 in budding yeast) is a highly multifunctional IDP associated with a range of protein complexes. However, it remains unknown if the different functions relate to different modified states. In this work, we show that Schizosaccharomyces pombe Dss1 is a substrate for casein kinase 2 in vitro, and we identify three phosphorylated threonines in its linker region separating two known disordered ubiquitin-binding motifs. Phosphorylations of the threonines had no effect on ubiquitin-binding but caused a slight destabilization of the C-terminal α-helix and mediated a direct interaction with the forkhead-associated (FHA) domain of the RING-FHA E3-ubiquitin ligase defective in mitosis 1 (Dma1). The phosphorylation sites are not conserved and are absent in human Dss1. Sequence analyses revealed that the Txx(E/D) motif, which is important for phosphorylation and Dma1 binding, is not linked to certain branches of the evolutionary tree. Instead, we find that the motif appears randomly, supporting the mechanism of ex nihilo evolution of novel motifs. In support of this, other threonine-based motifs, although frequent, are nonconserved in the linker, pointing to additional functions connected to this region. We suggest that Dss1 acts as an adaptor protein that docks to Dma1 via the phosphorylated FHA-binding motifs, while the C-terminal α-helix is free to bind mitotic septins, thereby stabilizing the complex. The presence of Txx(D/E) motifs in the disordered regions of certain septin subunits may be of further relevance to the formation and stabilization of these complexes.

Antibodies are used for many therapeutic and biotechnological purposes. Because the affinity of an antibody to the antigen is critical for clinical efficacy of pharmaceuticals, many affinity maturation strategies have been developed. Although we previously reported an affinity maturation strategy in which the association rate of the antibody toward its antigen is improved by introducing a cluster of arginine residues into the framework region of the antibody, the detailed molecular mechanism responsible for this improvement has been unknown. In this study, we introduced five arginine residues into an anti-hen egg white lysozyme antibody (HyHEL10) Fab fragment to create the R5-mutant and comprehensively characterized the interaction between antibody and antigen using thermodynamic analysis, X-ray crystallography, and molecular dynamics (MD) simulations. Our results indicate that introduction of charged residues strongly enhanced the association rate, as previously reported, and the antibody-antigen complex structure was almost the same for the R5-mutant and wild-type Fabs. The MD simulations indicate that the mutation increased conformational diversity in complementarity-determining region loops and thereby enhanced the association rate. These observations provide the molecular basis of affinity maturation by R5 mutation.

CTP synthases (CTPS) catalyze the de novo production of CTP using UTP, ATP, and l-glutamine with the anticancer drug metabolite gemcitabine-5'-triphosphate (dF-dCTP) being one of its most potent nucleotide inhibitors. To delineate the structural origins of this inhibition, we solved the structures of Escherichia coli CTPS (ecCTPS) in complex with CTP (2.0 Å), 2'-ribo-F-dCTP (2.0 Å), 2'-arabino-F-CTP (2.4 Å), dF-dCTP (2.3 Å), dF-dCTP and ADP (2.1 Å), and dF-dCTP and ATP (2.1 Å). These structures revealed that the increased binding affinities observed for inhibitors bearing the 2'-F-arabino group (dF-dCTP and F-araCTP), relative to CTP and F-dCTP, arise from interactions between the inhibitor's fluorine atom exploiting a conserved hydrophobic pocket formed by F227 and an interdigitating loop from an adjacent subunit (Q114-V115-I116). Intriguingly, crystal structures of ecCTPS•dF-dCTP complexes in the presence of select monovalent and divalent cations demonstrated that the in crystallo tetrameric assembly of wild-type ecCTPS was induced into a conformation similar to inhibitory ecCTPS filaments solely through the binding of Na+ -, Mg2+ -, or Mn2+ •dF-dCTP. However, in the presence of potassium, the dF-dCTP-bound structure is demetalated and in the low-affinity, non-filamentous conformation, like the conformation seen when bound to CTP and the other nucleotide analogues. Additionally, CTP can also induce the filament-competent conformation linked to high-affinity dF-dCTP binding in the presence of high concentrations of Mg2+ . This metal-dependent, compacted CTP pocket conformation therefore furnishes the binding environment responsible for the tight binding of dF-dCTP and provides insights for further inhibitor design.

l-Malate is a key flavor enhancer and acidulant in the food and beverage industry, particularly winemaking. Enzyme-based amperometric biosensors offer convenience for monitoring its concentration. However, only a small number of off-the-shelf malate-oxidizing enzymes have been used in previous devices. These typically have linear ranges poorly suited for the l-malate concentrations found in fruit processing and winemaking, making it necessary to use precisely diluted samples. Here, we describe a pipeline of database-mining, gene synthesis, recombinant expression, and spectrophotometric assays to characterize previously untested enzymes for their suitability in biosensors. The pipeline yielded a bespoke biocatalyst-the Ascaris suum malic enzyme carrying mutation R181Q [AsME(R181Q)]. Our first prototype with AsME(R181Q) had an ultra-wide linear range of 50-200 mM l-malate, corresponding to concentrations found in undiluted fruit juices (including grape). Changing the dication from Mg2+ to Mn2+ increased sensitivity five-fold and adding citrate (100 mM) increased it another six-fold, albeit decreasing the linear range to 1-10 mM. To our knowledge, this is the first time an l-malate biosensor with a tuneable combination of sensitivity and linear range has been described. The sensor response was also tested in the presence of various molecules abundant in juices and wines, with ascorbate shown to be a potent interferent. Interference was mitigated by the addition of ascorbate oxidase, allowing for differential measurements on an undiluted, untreated wine sample that corresponded well with commercial l-malate testing kits. Overall, this work demonstrates the power of an enzyme-centric approach for designing electrochemical biosensors with improved operational parameters and novel functionality.

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.

Organisms from all kingdoms of life synthesize L-serine (L-Ser) from 3-phosphoglycerate through the phosphorylated pathway, a three-step diversion of glycolysis. Phosphoserine aminotransferase (PSAT) catalyzes the intermediate step, the pyridoxal 5'-phosphate-dependent transamination of 3-phosphohydroxypyruvate and L-glutamate to O-phosphoserine (OPS) and α-ketoglutarate. PSAT is particularly relevant in the central nervous system of mammals because L-Ser is the metabolic precursor of D-serine, cysteine, phospholipids, and nucleotides. Several mutations in the human psat gene have been linked to serine deficiency disorders, characterized by severe neurological symptoms. Furthermore, PSAT is overexpressed in many tumors and this overexpression has been associated with poor clinical outcomes. Here, we report the detailed functional and structural characterization of the recombinant human PSAT. The reaction catalyzed by PSAT is reversible, with an equilibrium constant of about 10, and the enzyme is very efficient, with a kcat /Km of 5.9 × 106  M-1  s-1 , thus contributing in driving the pathway towards the products despite the extremely unfavorable first step catalyzed by 3-phosphoglycerate dehydrogenase. The 3D X-ray crystal structure of PSAT was solved in the substrate-free as well as in the OPS-bound forms. Both structures contain eight protein molecules in the asymmetric unit, arranged in four dimers, with a bound cofactor in each subunit. In the substrate-free form, the active site of PSAT contains a sulfate ion that, in the substrate-bound form, is replaced by the phosphate group of OPS. Interestingly, fast crystal soaking used to produce the substrate-bound form allowed the trapping of different intermediates along the catalytic cycle.

Peptidylglycine monooxygenase (PHM) is essential for the biosynthesis of many neuroendocrine peptides via a copper-dependent hydroxylation of a glycine-extended pro-peptide. The "canonical" mechanism requires the transfer of two electrons from one mononuclear copper (CuH, H-site) to a second mononuclear copper (CuM, M-site) which is the site of oxygen binding and catalysis. In most crystal structures the copper centers are separated by 11 Å of disordered solvent, but recent work has established that a PHM variant H108A forms a closed conformer in the presence of citrate with a reduced Cu-Cu site separation of ~4 Å. Here we report three new PHM structures where the H and M sites are separated by a longer distance of ~14 Å. Variation in Cu-Cu distance is the result of a rotation of the M subdomain about a hinge point centered on the pro199 -leu200 -ile201 triad which forms the linker between subdomains. The energetic cost of domain dynamics is likely small enough to allow free rotation of the subdomains relative to each other, adding credence to recent suggestions that an open-to-closed transition to form a binuclear oxygen binding intermediate is an essential element of catalysis. This inference would explain many experimental observations that are inconsistent with the current canonical mechanism including substrate-induced oxygen activation and isotope scrambling during the peroxide shunt.

Enhanced Disease Susceptibility 1 (EDS1), a key component of microbe-triggered immunity and effector-triggered immunity in most higher plants, forms functional heterodimeric complexes with its homologs Phytoalexin Deficient 4 (PAD4) or Senescence-associated Gene 101 (SAG101). Here, the crystal structure of VvEDS1Nterm , the N-terminal domain of EDS1 from Vitis vinifera, is reported, representing the first structure of an EDS1 entity beyond the model plant Arabidopsis thaliana. VvEDS1Nterm has an α/β-hydrolase fold, is similar to the N-terminal domain of A. thaliana EDS1 and forms stable homodimers in solution as well as in crystals. These VvEDS1Nterm homodimers are spatially incompatible with heterodimers with PAD4 or SAG101, they explain why VvEDS1Nterm does not interact with V. vinifera PAD4 according to gel filtration, and they serve as a guide to develop a plausible, albeit experimentally not verified model of full-length EDS1. VvEDS1Nterm is a splicing variant comprising two of three exons of the VvEDS1 gene. It originates from a naturally occurring mRNA, in which the first of two introns was removed while the second one containing a stop codon close to the exon/intron border was retained. This is a potential case of intron retention and the first report of this phenomenon in the context of EDS1. Its biological significance has not yet been clarified, nor has the question if a VvEDS1Nterm protein with a specific function can occur under physiological conditions.

Programmed death-1 (PD-1), an immune checkpoint receptor, is expressed on activated lymphocytes, macrophages, and some types of tumor cells. While PD-1+ cells have been implicated in outcomes of cancer immunity, autoimmunity, and chronic infections, the exact roles of these cells in various physiological and pathological processes remain elusive. Molecules that target and deplete PD-1+ cells would be instrumental in defining the roles unambiguously. Previously, an immunotoxin has been generated for the depletion of PD-1+ cells though its usage is impeded by its low production yield. Thus, a more practical molecular tool is desired to deplete PD-1+ cells and to examine functions of these cells. We designed and generated a novel anti-PD1 diphtheria immunotoxin, termed PD-1 DIT, targeting PD-1+ cells. PD-1 DIT is comprised of two single chain variable fragments (scFv) derived from an anti-PD-1 antibody, coupled with the catalytic and translocation domains of the diphtheria toxin. PD-1 DIT was produced using a yeast expression system that has been engineered to efficiently produce protein toxins. The yield of PD-1 DIT reached 1-2 mg/L culture, which is 10 times higher than the previously reported immunotoxin. Flow cytometry and confocal microscopy analyses confirmed that PD-1 DIT specifically binds to and enters PD-1+ cells. The binding avidities between PD-1 DIT and two PD-1+ cell lines are approximately 25 nM. Moreover, PD-1 DIT demonstrated potent cytotoxicity toward PD-1+ cells, with a half maximal effective concentration (EC50 ) value of 1 nM. In vivo experiments further showed that PD-1 DIT effectively depleted PD-1+ cells and enabled mice inoculated with PD-1+ tumor cells to survive throughout the study. Our findings using PD-1 DIT revealed the critical role of pancreatic PD-1+ T cells in the development of type-1 diabetes (T1D). Additionally, we observed that PD-1 DIT treatment ameliorated relapsing-remitting experimental autoimmune encephalomyelitis (RR-EAE), a mouse model of relapsing-remitting multiple sclerosis (RR-MS). Lastly, we did not observe significant hepatotoxicity in mice treated with PD-1 DIT, which had been reported for other immunotoxins derived from the diphtheria toxin. With its remarkable selective and potent cytotoxicity toward PD-1+ cells, coupled with its high production yield, PD-1 DIT emerges as a powerful biotechnological tool for elucidating the physiological roles of PD-1+ cells. Furthermore, the potential of PD-1 DIT to be developed into a novel therapeutic agent becomes evident.

The interaction between the receptor-binding domain (RBD) of the spike glycoprotein of SARS-CoV-2 and the peptidase domain of the human angiotensin-converting enzyme 2 (ACE2) allows the first specific contact at the virus-cell interface making it the main target of neutralizing antibodies. Here, we show a unique and cost-effective protocol using Drosophila S2 cells to produce both RBD and soluble human ACE2 peptidase domain (shACE2) as thermostable proteins, purified via Strep-tag with yields >40 mg L-1 in a laboratory scale. Furthermore, we demonstrate its binding with KD values in the lower nanomolar range (independently of Strep-tag removal) and its capability to be blocked by serum antibodies in a competition ELISA with Strep-Tactin-HRP as a proof-of-concept. In addition, we assess the capacity of RBD to bind native dimeric ACE2 overexpressed in human cells and its antigen properties with specific serum antibodies. Finally, for completeness, we analyzed RBD microheterogeneity associated with glycosylation and negative charges, with negligible effect on binding either with antibodies or shACE2. Our system represents an accessible and reliable tool for designing in-house surrogate virus neutralization tests (sVNTs), enabling the rapid characterization of neutralizing humoral responses elicited against vaccines or infection, especially in the absence of facilities to conduct virus neutralization tests. Moreover, our biophysical and biochemical characterization of RBD and shACE2 produced in S2 cells lays the groundwork for adapting to different variants of concern (VOCs) to study humoral responses elicited against different VOCs and vaccine formulations.

Methicillin-resistant Staphylococcus aureus (MRSA) causes life-threatening human infections. Bacteriophage-encoded endolysins degrade the cell walls of Gram-positive bacteria by selectively hydrolyzing the peptidoglycan layer and thus are promising candidates to combat bacterial infections. PlyGRCS, the S. aureus-specific bacteriophage endolysin, contains a catalytic CHAP domain and a cell-wall binding SH3_5 domain connected by a linker. Here, we show the crystal structure of full-length PlyGRCS refined to 2.1 Å resolution. In addition, a serendipitous finding revealed that PlyGRCS binds to cold-shock protein C (CspC) by interacting with its CHAP and SH3_5 domains. CspC is an RNA chaperone that plays regulatory roles by conferring bacterial adaptability to various stress conditions. PlyGRCS has substantial lytic activity against S. aureus and showed only minimal change in its lytic activity in the presence of CspC. Whereas the PlyGRCS-CspC complex greatly reduced CspC-nucleic acid binding, the aforesaid complex may downregulate the CspC function during bacterial infection. Overall, the crystal structure and biochemical results of PlyGRCS provide a molecular basis for the bacteriolytic activity of PlyGRCS against S. aureus.

The interaction between the transcription factor p53 and the ubiquitin ligase MDM2 results in the degradation of p53 and is well-studied in cancer biology and drug development. Available sequence data suggest that both p53 and MDM2-family proteins are present across the animal kingdom. However, the interacting regions are missing in some animal groups, and it is not clear whether MDM2 interacts with, and regulates p53 in all species. We used phylogenetic analyses and biophysical measurements to examine the evolution of affinity between the interacting protein regions: a conserved 12-residue intrinsically disordered binding motif in the p53 transactivation domain (TAD) and the folded SWIB domain of MDM2. The affinity varied significantly across the animal kingdom. The p53TAD/MDM2 interaction among jawed vertebrates displayed high affinity, in particular for chicken and human proteins (KD around 0.1 μM). The affinity of the bay mussel p53TAD/MDM2 complex was lower (KD  = 15 μM) and those from a placozoan, an arthropod, and a jawless vertebrate were very low or non-detectable (KD  > 100 μM). Binding experiments with reconstructed ancestral p53TAD/MDM2 variants suggested that a micromolar affinity interaction was present in the ancestral bilaterian animal and was later enhanced in tetrapods while lost in other linages. The different evolutionary trajectories of p53TAD/MDM2 affinity during speciation demonstrate high plasticity of motif-mediated interactions and the potential for rapid adaptation of p53 regulation during times of change. Neutral drift in unconstrained disordered regions may underlie the plasticity and explain the observed low sequence conservation in TADs such as p53TAD.

The outer leaflet of the outer membrane (OM) of bacteria such as Escherichia coli, Pseudomonas aeruginosa, and other important pathogens is largely composed of lipopolysaccharide (LPS), which is essential to nearly all Gram-negative bacteria. LPS is transported to the outer leaflet of the OM through a yet unknown mechanism by seven proteins that comprise the LPS transport system. LptA, the only entirely periplasmic Lpt protein, bridges the periplasmic space between the IM LptB2 FGC and the OM LptDE complexes. LptA is postulated to protect the hydrophobic acyl chains of LPS as it crosses the hydrophilic periplasm, is essential to cell viability, and contains many conserved residues distributed across the protein. To identify which side chains are required for function of E. coli LptA in vivo, we performed a systematic, unbiased, high-throughput screen of the effect of 172 single alanine substitutions on cell viability utilizing an engineered BL21 derivative with a chromosomal knockout of the lptA gene. Remarkably, LptA is highly tolerant to amino acid substitution with alanine. Only four alanine mutants could not complement the chromosomal knockout; CD spectroscopy showed that these substitutions resulted in proteins with significantly altered secondary structure. In addition, 29 partial loss-of-function mutants were identified that led to OM permeability defects; interestingly, these sites were solely located within β-strands of the central core of the protein and each resulted in misfolding of the protein. Therefore, no single residue within LptA is responsible for LPS binding, supporting previous EPR spectroscopy data indicating that sites across the entire protein work in concert to bind and transport LPS.

Membrane-embedded transporters impart essential functions to cells as they mediate sensing and the uptake and extrusion of nutrients, waste products, and effector molecules. Promiscuous multidrug exporters are implicated in resistance to drugs and antibiotics and are highly relevant for microbial engineers who seek to enhance the tolerance of cell factory strains to hydrophobic bioproducts. Here, we report on the identification of small multidrug resistance (SMR) transporters in early-branching anaerobic fungi (Neocallimastigomycetes). The SMR class of transporters is commonly found in bacteria but has not previously been reported in eukaryotes. In this study, we show that SMR transporters from anaerobic fungi can be produced heterologously in the model yeast Saccharomyces cerevisiae, demonstrating the potential of these proteins as targets for further characterization. The discovery of these novel anaerobic fungal SMR transporters offers a promising path forward to enhance bioproduction from engineered microbial strains.

UHRF1 is an essential chromatin protein required for DNA methylation maintenance, mammalian development, and gene regulation. We investigated the Tandem-Tudor domain (TTD) of human UHRF1 that is known to bind H3K9me2/3 histones and is a major driver of UHRF1 localization in cells. We verified binding to H3K9me2/3 but unexpectedly discovered stronger binding to H3 peptides and mononucleosomes containing K9me2/3 with additional K4me1. We investigated the combined binding of TTD to H3K4me1-K9me2/3 versus H3K9me2/3 alone, engineered mutants with specific and differential changes of binding, and discovered a novel read-out mechanism for H3K4me1 in an H3K9me2/3 context that is based on the interaction of R207 with the H3K4me1 methyl group and on counting the H-bond capacity of H3K4. Individual TTD mutants showed up to a 10,000-fold preference for the double-modified peptides, suggesting that after a conformational change, WT TTD could exhibit similar effects. The frequent appearance of H3K4me1-K9me2 regions in human chromatin demonstrated in our TTD chromatin pull-down and ChIP-western blot data suggests that it has specific biological roles. Chromatin pull-down of TTD from HepG2 cells and full-length murine UHRF1 ChIP-seq data correlate with H3K4me1 profiles indicating that the H3K4me1-K9me2/3 interaction of TTD influences chromatin binding of full-length UHRF1. We demonstrate the H3K4me1-K9me2/3 specific binding of UHRF1-TTD to enhancers and promoters of cell-type-specific genes at the flanks of cell-type-specific transcription factor binding sites, and provided evidence supporting an H3K4me1-K9me2/3 dependent and TTD mediated downregulation of these genes by UHRF1. All these findings illustrate the important physiological function of UHRF1-TTD binding to H3K4me1-K9me2/3 double marks in a cellular context.

Bacterial lytic transglycosylases (LTs) contribute to peptidoglycan cell wall metabolism and are potential drug targets to potentiate β-lactam antibiotics to overcome antibiotic resistance. Since LT inhibitor development is underexplored, we probed 15 N-acetyl-containing heterocycles in a structure-guided fashion for their ability to inhibit and bind to the Campylobacter jejuni LT Cj0843c. Ten GlcNAc analogs were synthesized with substitutions at the C1 position, with two having an additional modification at the C4 or C6 position. Most of the compounds showed weak inhibition of Cj0843c activity. Compounds with alterations at the C4 position, replacing the -OH with a -NH2 , and C6 position, the addition of a -CH3 , yielded improved inhibitory efficacy. All 10 GlcNAc analogs were crystallographically analyzed via soaking experiments using Cj0843c crystals and found to bind to the +1 +2 saccharide subsites with one of them additionally binding to the -2 -1 subsite region. We also probed other N-acetyl-containing heterocycles and found that sialidase inhibitors N-acetyl-2,3-dehydro-2-deoxyneuraminic acid and siastatin B inhibited Cj0843c weakly and crystallographically bound to the -2 -1 subsites. Analogs of the former also showed inhibition and crystallographic binding and included zanamivir amine. This latter set of heterocycles positioned their N-acetyl group in the -2 subsite with additional moieties interacting in the -1 subsite. Overall, these results could provide novel opportunities for LT inhibition via exploring different subsites and novel scaffolds. The results also increased our mechanistic understanding of Cj0843c regarding peptidoglycan GlcNAc subsite binding preferences and ligand-dependent modulation of the protonation state of the catalytic E390.

Abrupt aggregation of misfolded proteins is the underlying molecular cause of numerous severe pathologies including Alzheimer's and Parkinson's diseases. Protein aggregation yields small oligomers that can later propagate into amyloid fibrils, β-sheet-rich structures with a variety of topologies. A growing body of evidence suggests that lipids play an important role in abrupt aggregation of misfolded proteins. In this study, we investigate the roles of length and saturation of fatty acids (FAs) in phosphatidylserine (PS), an anionic lipid that is responsible for the recognition of apoptotic cells by macrophages, in lysozyme aggregation. We found that both the length and saturation of FAs in PS contribute to the aggregation rate of insulin. PS with 14-carbon-long FAs (14:0) enabled a much stronger acceleration of protein aggregation compared to PS with 18-carbon-long FAs (18:0). Our results demonstrate that the presence of double bonds in FAs accelerated the rate of insulin aggregation relative to PS with fully saturated FAs. Biophysical methods revealed morphological and structural differences in lysozyme aggregates grown in the presence of PS with varying lengths and FA saturation. We also found that such aggregates exerted diverse cell toxicities. These results demonstrate that the length and saturation of FAs in PS can uniquely alter the stability of misfolded proteins on lipid membranes.

Intein enzymes catalyze the splicing of their flanking polypeptide chains and have found tremendous biotechnological applications. Their terminal residues form the catalytic core and participate in the splicing reaction. Hence, the neighboring N- and C-terminal extein residues influence the catalytic rate. As these extein residues vary depending on the substrate identity, we tested the influence of 20 amino acids at these sites in the Spl DnaX intein and observed significant variation of spliced product as well as N- and C-terminus cleavage product formation. We investigated the dependence of these reactions on the extein residues by molecular dynamics (MD) simulations on eight extein variants, and found that the conformational sampling of the active-site residues of the intein enzyme differed among these extein variants. We found that the extein variants that sample higher population of near-attack conformers (NACs) of the active-site residues undergo higher product formation in our activity assays. Ground state conformers that closely resemble the transition state are referred to as NACs. Very good correlation was observed between the NAC populations from the MD simulations of eight extein variants and the corresponding product formation from our activity assays. Furthermore, this molecular detail enabled us to elucidate the mechanistic roles of several conserved active-site residues in the splicing reaction. Overall, this study shows that the catalytic power of Spl DnaX intein enzyme, and most likely other inteins, depends on the efficiency of formation of NACs in the ground state, which is further modulated by the extein residues.

Pulse EPR measurements provide information on distances and distance distributions in proteins but require the incorporation of pairs of spin labels that are usually attached to engineered cysteine residues. In previous work, we demonstrated that efficient in vivo labeling of the Escherichia coli outer membrane vitamin B12 transporter, BtuB, could only be achieved using strains defective in the periplasmic disulfide bond formation (Dsb) system. Here, we extend these in vivo measurements to FecA, the E. coli ferric citrate transporter. As seen for BtuB, pairs of cysteines cannot be labeled when the protein is present in a standard expression strain. However, incorporating plasmids that permit an arabinose induced expression of FecA into a strain defective in the thiol disulfide oxidoreductase, DsbA, enables efficient spin-labeling and pulse EPR of FecA in cells. A comparison of the measurements made on FecA in cells with measurements made in reconstituted phospholipid bilayers suggests that the cellular environment alters the behavior of the extracellular loops of FecA. In addition to these in situ EPR measurements, the use of a DsbA minus strain for the expression of BtuB improves the EPR signals and pulse EPR data obtained in vitro from BtuB that is labeled, purified, and reconstituted into phospholipid bilayers. The in vitro data also indicate the presence of intermolecular BtuB-BtuB interactions, which had not previously been observed in a reconstituted bilayer system. This result suggests that in vitro EPR measurements on other outer membrane proteins would benefit from protein expression in a DsbA minus strain.

We report the discovery of the androgen receptor missense mutation V770D, that was found in two sisters suffering from complete androgen insensitivity. Experimental validation of AR V770 variants demonstrated that AR V770D was transcriptionally inactive due to the inability to dimerize and a reduced ligand binding affinity. The more conservative AR V770A variant showed a dimerization defect at low levels of DHT with a partial recovery of the transcriptional activity and of the receptor's ability to dimerize when increasing the DHT levels. With V770 located outside of the proposed LBD dimerization interface of the AR LBD homodimer crystal structure, the effects of the V770A mutation on AR dimerization were unexpected. We therefore explored whether the AR LBD dimerization interface would be better described by an alternative dimerization mode based on available human homodimeric LBD crystal structures of other nuclear receptors. Superposition of the monomeric AR LBD in the homodimeric crystal structures of GR, PR, ER, CAR, TRβ, and HNF-4α showed that the GR-like LBD dimer model was energetically the most stable. Moreover, V770 was a key energy residue in the GR-like LBD dimer while it was not involved in the stabilization of the AR LBD homodimer according to the crystal structure. Additionally, the observation that 4 AIS mutations impacted the stability of the AR LBD dimer while 16 mutations affected the GR-like LBD dimer, suggested that the AR LBD dimer crystal is a snapshot of a breathing AR LBD homodimer that can transition into the GR-like LBD dimer model.

ADP-glucose pyrophosphorylase is a key regulatory enzyme involved in starch and glycogen synthesis in plants and bacteria, respectively. It has been hypothesized that inter-subunit communications are important for the allosteric effect in this enzyme. However, no specific interactions have been identified as part of the regulatory signal. The enzyme from Agrobacterium tumefaciens is a homotetramer allosterically regulated by fructose 6-phosphate and pyruvate. Three pairs of distinct subunit-subunit interfaces are present. Here we focus on an interface that features two symmetrical interactions between Arg11 and Asp141 from one subunit with residues Asp141 and Arg11 of the neighbor subunit, respectively. Previously, scanning mutagenesis showed that a mutation at the Arg11 position disrupted the activation of the enzyme. Considering the distance of these residues from the allosteric and catalytic sites, we hypothesized that the interaction between Arg11 and Asp141 is critical for allosteric signaling rather than effector binding. To prove our hypothesis, we mutated those two sites (D141A, D141E, D141N, D141R, R11D, and R11K) and performed kinetic and binding analysis. Mutations that altered the charge affected the regulation the most. To prove that the interaction per se (rather than the presence of specific residues) is critical, we partially rescued the R11D protein by introducing a second mutation (R11D/D141R). This could not restore the activator effect on kcat , but it did rescue the effect on substrate affinity. Our results indicate the critical functional role of Arg11 and Asp141 to relay the allosteric signal in this subunit interface.

Abrupt aggregation of misfolded proteins is the underlying molecular cause of Alzheimer disease (AD) and Parkinson disease (PD). Both AD and PD are severe pathologies that affect millions of people around the world. A small 42 amino acid long peptide, known as amyloid β (Aβ), aggregates in the frontal cortex of AD patients forming oligomers and fibrils, highly toxic protein aggregates that cause progressive neuron death. Similar aggregates of α-synuclein (α-Syn), a small protein that facilitates neurotransmitter release, are observed in the midbrain, hypothalamus, and thalamus of people with PD. In this study, we utilized the innovative nano-Infrared imaging technique to investigate the structural organization of individual Aβ and α-syn fibrils postmortem extracted from brains of AD and PD patients, respectively. We observed two morphologically different Aβ and α-Syn fibril polymorphs in each patient's brain. One had twisted topology, whereas another exhibited flat tape-like morphology. We found that both polymorphs shared the same parallel β-sheet-dominated secondary structure. These findings suggested that both fibril polymorphs were built from structurally similar if not identical filaments that coiled forming twisted fibrils or associated side-by-side in the case of straight Aβ and α-Syn fibrils. Nano-Infrared analysis of individual protein aggregates also revealed the presence of lipids in the structure of both twisted and tape-like α-Syn fibrils that were not observed in any of the Aβ fibril polymorphs. These findings demonstrate that lipid membranes can play a critically important role in the onset and progression of PD.

Human aromatic amino acid decarboxylase (AADC) is a pyridoxal 5'-phosphate-dependent enzyme responsible for the biosynthesis of dopamine and serotonin, essential neurotransmitters involved in motor and cognitive abilities. Mutations in its gene lead to AADC deficiency, a monogenic rare neurometabolic childhood parkinsonism characterized by severe motor and neurodevelopmental symptoms. Here, for the first time, we solved the crystal structure of human holoAADC in the internal aldimine (1.9 Å) and in the external aldimine (2.4 Å) of the substrate analog L-Dopa methylester. In this intermediate, the highly flexible AADC catalytic loop (CL) is captured in a closed state contacting all protein domains. In addition, each active site, composed by residues of both subunits, is connected to the other through weak interactions and a central cavity. By combining crystallographic analyses with all-atom and coarse-grained molecular dynamics simulations, SAXS investigations and limited proteolysis experiments, we realized that the functionally obligate homodimeric AADC enzyme in solution is an elongated, asymmetric molecule, where the fluctuations of the CL are coupled to flexibility at the edge between the N-terminal and C-terminal domains. The structural integrity of this peripheral protein region is essential to catalysis, as assessed by both artificial and 37 AADC deficiency pathogenic variants leading to the interpretation that structural dynamics in protein regions far from the active site is essential for CL flexibility and the acquirement of a correct catalytically competent structure. This could represent the molecular basis for pathogenicity prediction in AADC deficiency.