Characterization of the Domain Orientations of E. coli 5'-Nucleotidase by Fitting an Ensemble of Conformers to DEER Distance Distributions

Molecular Dynamics Simulations of the Human Ecto-5'-Nucleotidase (h-ecto-5'-NT, CD73): Insights into Protein Flexibility and Binding Site Dynamics

Human ecto-5'-nucleotidase (h-ecto-5'-NT, CD73) is a homodimeric Zn2+-binding metallophosphoesterase that hydrolyzes adenosine 5'-monophosphate (5'-AMP) to adenosine and phosphate. h-Ecto-5'-NT is a key enzyme in purinergic signaling pathways and has been recognized as a promising biological target for several diseases, including cancer and inflammatory, infectious, and autoimmune diseases. Despite its importance as a biological target, little is known about h-ecto-5'-NT dynamics, which poses a considerable challenge to the design of inhibitors of this target enzyme. Here, to explore h-ecto-5'-NT flexibility, all-atom unbiased molecular dynamics (MD) simulations were performed. Remarkable differences in the dynamics of the open (catalytically inactive) and closed (catalytically active) conformations of the apo-h-ecto-5'-NT were observed during the simulations, and the nucleotide analogue inhibitor AMPCP was shown to stabilize the protein structure in the closed conformation. Our results suggest that the large and complex domain motion that enables the h-ecto-5'-NT open/closed conformational switch is slow, and therefore, it could not be completely captured within the time scale of our simulations. Nonetheless, we were able to explore the faster dynamics of the h-ecto-5'-NT substrate binding site, which is mainly located at the C-terminal domain and well conserved among the protein's open and closed conformations. Using the TRAPP ("Transient Pockets in Proteins") approach, we identified transient subpockets close to the substrate binding site. Finally, conformational states of the substrate binding site with higher druggability scores than the crystal structure were identified. In summary, our study provides valuable insights into h-ecto-5'-NT structural flexibility, which can guide the structure-based design of novel h-ecto-5'-NT inhibitors.

A novel approach to modeling side chain ensembles of the bifunctional spin label RX

We introduce a novel approach to modeling side chain ensembles of bifunctional spin labels. This approach utilizes rotamer libraries to generate side chain conformational ensembles. Because the bifunctional label is constrained by two attachment sites, the label is split into two monofunctional rotamers which are first attached to their respective sites, then rejoined by a local optimization in dihedral space. We validate this method against a set of previously published experimental data using the bifunctional spin label, RX. This method is relatively fast and can readily be used for both experimental analysis and protein modeling, providing significant advantages over modeling bifunctional labels with molecular dynamics simulations. Use of bifunctional labels for site directed spin labeling (SDSL) electron paramagnetic resonance (EPR) spectroscopy dramatically reduces label mobility, which can significantly improve resolution of small changes in protein backbone structure and dynamics. Coupling the use of bifunctional labels with side chain modeling methods allows for improved quantitative application of experimental SDSL EPR data to protein modeling.

Benchmark Test and Guidelines for DEER/PELDOR Experiments on Nitroxide-Labeled Biomolecules

Distance distribution information obtained by pulsed dipolar EPR spectroscopy provides an important contribution to many studies in structural biology. Increasingly, such information is used in integrative structural modeling, where it delivers unique restraints on the width of conformational ensembles. In order to ensure reliability of the structural models and of biological conclusions, we herein define quality standards for sample preparation and characterization, for measurements of distributed dipole-dipole couplings between paramagnetic labels, for conversion of the primary time-domain data into distance distributions, for interpreting these distributions, and for reporting results. These guidelines are substantiated by a multi-laboratory benchmark study and by analysis of data sets with known distance distribution ground truth. The study and the guidelines focus on proteins labeled with nitroxides and on double electron-electron resonance (DEER aka PELDOR) measurements and provide suggestions on how to proceed analogously in other cases.

Reconstruction of Coupled Intra- and Interdomain Protein Motion from Nuclear and Electron Magnetic Resonance

Proteins composed of multiple domains allow for structural heterogeneity and interdomain dynamics that may be vital for function. Intradomain structures and dynamics can influence interdomain conformations and vice versa. However, no established structure determination method is currently available that can probe the coupling of these motions. The protein Pin1 contains separate regulatory and catalytic domains that sample "extended" and "compact" states, and ligand binding changes this equilibrium. Ligand binding and interdomain distance have been shown to impact the activity of Pin1, suggesting interdomain allostery. In order to characterize the conformational equilibrium of Pin1, we describe a novel method to model the coupling between intra- and interdomain dynamics at atomic resolution using multistate ensembles. The method uses time-averaged nuclear magnetic resonance (NMR) restraints and double electron-electron resonance (DEER) data that resolve distance distributions. While the intradomain calculation is primarily driven by exact nuclear Overhauser enhancements (eNOEs), J couplings, and residual dipolar couplings (RDCs), the relative domain distribution is driven by paramagnetic relaxation enhancement (PREs), RDCs, interdomain NOEs, and DEER. Our data support a 70:30 population of the compact and extended states in apo Pin1. A multistate ensemble describes these conformations simultaneously, with distinct conformational differences located in the interdomain interface stabilizing the compact or extended states. We also describe correlated conformations between the catalytic site and interdomain interface that may explain allostery driven by interdomain contact.

Substrate binding modes of purine and pyrimidine nucleotides to human ecto-5'-nucleotidase (CD73) and inhibition by their bisphosphonic acid derivatives

Human ecto-5-nucleotidase (CD73) is involved in purinergic signalling, which influences a diverse range of biological processes. CD73 hydrolyses AMP and is the major control point for the levels of extracellular adenosine. Inhibitors of CD73 thus block the immunosuppressive action of adenosine, a promising approach for cancer immunotherapy. Interestingly, ADP and ATP are competitive inhibitors of CD73, with the most potent small-molecule inhibitors to date being non-hydrolysable ADP analogues. While AMP is the major substrate of the enzyme, CD73 has been reported to hydrolyse other 5′-nucleoside monophosphates. Based on a fragment screening campaign at the BESSY II synchrotron, we present the binding modes of various deoxyribo- and ribonucleoside monophosphates and of four additional fragments binding to the nucleoside binding site of the open form of the enzyme. Kinetic analysis of monophosphate hydrolysis shows that ribonucleotide substrates are favoured over their deoxyribose equivalents with AMP being the best substrate. We characterised the initial step of AMP hydrolysis, the binding mode of AMP to the open conformation of CD73 and compared that to other monophosphate substrates. In addition, the inhibitory activity of various bisphosphonic acid derivatives of nucleoside diphosphates was determined. Although AMPCP remains the most potent inhibitor, replacement of the adenine base with other purines or with pyrimidines increases the K_i value only between twofold and sixfold. On the other hand, these nucleobases offer new opportunities to attach substituents for improved pharmacological properties.

Methodology for rigorous modeling of protein conformational changes by Rosetta using DEER distance restraints

We describe an approach for integrating distance restraints from Double Electron-Electron Resonance (DEER) spectroscopy into Rosetta with the purpose of modeling alternative protein conformations from an initial experimental structure. Fundamental to this approach is a multilateration algorithm that harnesses sets of interconnected spin label pairs to identify optimal rotamer ensembles at each residue that fit the DEER decay in the time domain. Benchmarked relative to data analysis packages, the algorithm yields comparable distance distributions with the advantage that fitting the DEER decay and rotamer ensemble optimization are coupled. We demonstrate this approach by modeling the protonation-dependent transition of the multidrug transporter PfMATE to an inward facing conformation with a deviation to the experimental structure of less than 2Å Cα RMSD. By decreasing spin label rotamer entropy, this approach engenders more accurate Rosetta models that are also more closely clustered, thus setting the stage for more robust modeling of protein conformational changes.

Substrate Specificity of Chimeric Enzymes Formed by Interchange of the Catalytic and Specificity Domains of the 5'-Nucleotidase UshA and the 3'-Nucleotidase CpdB

The 5′-nucleotidase UshA and the 3′-nucleotidase CpdB from Escherichia coli are broad-specificity phosphohydrolases with similar two-domain structures. Their N-terminal domains (UshA_Ndom and CpdB_Ndom) contain the catalytic site, and their C-terminal domains (UshA_Cdom and CpdB_Cdom) contain a substrate-binding site responsible for specificity. Both enzymes show only partial overlap in their substrate specificities. So, it was decided to investigate the catalytic behavior of chimeras bearing the UshA catalytic domain and the CpdB specificity domain, or vice versa. UshA_Ndom–CpdB_Cdom and CpdB_Ndom–UshA_Cdom were constructed and tested on substrates specific to UshA (5′-AMP, CDP-choline, UDP-glucose) or to CpdB (3′-AMP), as well as on 2′,3′-cAMP and on the common phosphodiester substrate bis-4-NPP (bis-4-nitrophenylphosphate). The chimeras did show neither 5′-nucleotidase nor 3′-nucleotidase activity. When compared to UshA, UshA_Ndom–CpdB_Cdom conserved high activity on bis-4-NPP, some on CDP-choline and UDP-glucose, and displayed activity on 2′,3′-cAMP. When compared to CpdB, CpdB_Ndom–UshA_Cdom conserved phosphodiesterase activities on 2′,3′-cAMP and bis-4-NPP, and gained activity on the phosphoanhydride CDP-choline. Therefore, the non-nucleotidase activities of UshA and CpdB are not fully dependent on the interplay between domains. The specificity domains may confer the chimeras some of the phosphodiester or phosphoanhydride selectivity displayed when associated with their native partners. Contrarily, the nucleotidase activity of UshA and CpdB depends strictly on the interplay between their native catalytic and specificity domains.

Inference of Joint Conformational Distributions from Separately Acquired Experimental Measurements

Flexible proteins serve vital roles in a multitude of biological processes. However, determining their full conformational ensembles is extremely difficult because this requires detailed knowledge about the heterogeneity of the protein's degrees of freedom. Label-based experiments such as double electron-electron resonance (DEER) are very useful in studying flexible proteins, as they provide distributional data on heterogeneity. These experiments are typically performed separately, so information about correlation between distributions is lost. We have developed a method to recover correlation information using nonequilibrium work estimates in molecular dynamics refinement. We tested this method on a simple model of an alternating-access transporter for which the true joint distributions are known, and it successfully recovered the true joint distribution. We also applied our method to the protein syntaxin-1a, where it discarded physically implausible conformations. Our method thus provides a way to recover correlation structure in separate experimental measurements of conformational ensembles and refines the resulting structural ensemble.

Molecular Dissection of Escherichia coli CpdB: Roles of the N Domain in Catalysis and Phosphate Inhibition, and of the C Domain in Substrate Specificity and Adenosine Inhibition

CpdB is a 3′-nucleotidase/2′3′-cyclic nucleotide phosphodiesterase, active also with reasonable efficiency on cyclic dinucleotides like c-di-AMP (3′,5′-cyclic diadenosine monophosphate) and c-di-GMP (3′,5′-cyclic diadenosine monophosphate). These are regulators of bacterial physiology, but are also pathogen-associated molecular patterns recognized by STING to induce IFN-β response in infected hosts. The cpdB gene of Gram-negative and its homologs of gram-positive bacteria are virulence factors. Their protein products are extracytoplasmic enzymes (either periplasmic or cell–wall anchored) and can hydrolyze extracellular cyclic dinucleotides, thus reducing the innate immune responses of infected hosts. This makes CpdB(-like) enzymes potential targets for novel therapeutic strategies in infectious diseases, bringing about the necessity to gain insight into the molecular bases of their catalytic behavior. We have dissected the two-domain structure of Escherichia coli CpdB to study the role of its N-terminal and C-terminal domains (CpdB_Ndom and CpdB_Cdom). The specificity, kinetics and inhibitor sensitivity of point mutants of CpdB, and truncated proteins CpdB_Ndom and CpdB_Cdom were investigated. CpdB_Ndom contains the catalytic site, is inhibited by phosphate but not by adenosine, while CpdB_Cdom is inactive but contains a substrate-binding site that determines substrate specificity and adenosine inhibition of CpdB. Among CpdB substrates, 3′-AMP, cyclic dinucleotides and linear dinucleotides are strongly dependent on the CpdB_Cdom binding site for activity, as the isolated CpdB_Ndom showed much-diminished activity on them. In contrast, 2′,3′-cyclic mononucleotides and bis-4-nitrophenylphosphate were actively hydrolyzed by CpdB_Ndom, indicating that they are rather independent of the CpdB_Cdom binding site.

Significant influence of phosphorus resources on the growth and alkaline phosphatase activities of Microcystis aeruginosa

It is well-accepted that phosphorus, particularly orthophosphate, is a determinant factor in aquatic eutrophication. However, numerous kinds of phosphorus sources exist in real world scenario, and limited studies have characterized the pairwise relationships among abundant different phosphorus sources and the physiological behaviour of algae. The present study developed a high-throughput assay to investigate the effects of 59 different phosphorus sources (equal initial concentration of total phosphorus) on the growth and alkaline phosphatase (AKP) activities of Microcystis aeruginosa, a model cyanobacteria whose predominance holds sway in lake eutrophication. M. aeruginosa cultivated with nucleoside monophosphates (NMPs) had higher growth, relative AKP activities and residual orthophosphate, which were positively intercorrelated. Oppositely, non-NMPs cultivation of M. aeruginosa led to negative relationships between the relative AKP activities and their growth or residual orthophosphate. These results indicated distinct mechanisms for M. aeruginosa to utilize different phosphorus sources in real-world scenario, and both phosphorus source and content are determinant factors on the growth and physiological behaviour of M. aeruginosa. Given the complicated and vast phosphorus pool in the natural environment, phosphorus resources might significantly alter the abundance and physiological behaviour of M. aeruginosa and other bloom-forming algae, then influence the phytoplanktonic community structure and affect the possibility and intensity of algal bloom. Our work hints the underestimation of the restriction factors in lake eutrophication and provides a new tool to study the driven forces of phytoplanktonic community dynamics as phosphorus from both internal and external sources.

Rapid Simulation of Unprocessed DEER Decay Data for Protein Fold Prediction

Despite advances in sampling and scoring strategies, Monte Carlo modeling methods still struggle to accurately predict de novo the structures of large proteins, membrane proteins, or proteins of complex topologies. Previous approaches have addressed these shortcomings by leveraging sparse distance data gathered using site-directed spin labeling and electron paramagnetic resonance spectroscopy to improve protein structure prediction and refinement outcomes. However, existing computational implementations entail compromises between coarse-grained models of the spin label that lower the resolution and explicit models that lead to resource-intense simulations. These methods are further limited by their reliance on distance distributions, which are calculated from a primary refocused echo decay signal and contain uncertainties that may require manual refinement. Here, we addressed these challenges by developing RosettaDEER, a scoring method within the Rosetta software suite capable of simulating double electron-electron resonance spectroscopy decay traces and distance distributions between spin labels fast enough to fold proteins de novo. We demonstrate that the accuracy of resulting distance distributions match or exceed those generated by more computationally intensive methods. Moreover, decay traces generated from these distributions recapitulate intermolecular background coupling parameters even when the time window of data collection is truncated. As a result, RosettaDEER can discriminate between poorly folded and native-like models by using decay traces that cannot be accurately converted into distance distributions using regularized fitting approaches. Finally, using two challenging test cases, we demonstrate that RosettaDEER leverages these experimental data for protein fold prediction more effectively than previous methods. These benchmarking results confirm that RosettaDEER can effectively leverage sparse experimental data for a wide array of modeling applications built into the Rosetta software suite.

Roles of Phosphorus Sources in Microbial Community Assembly for the Removal of Organic Matters and Ammonia in Activated Sludge

Various phosphorus sources are utilized by microbes in WWTPs, eventually affecting microbial assembly and functions. This study identified the effects of phosphorus source on microbial communities and functions in the activated sludge. By cultivation with 59 phosphorus sources, including inorganic phosphates (IP), nucleoside-monophosphates (NMP), cyclic-nucleoside-monophosphates (cNMP), and other organophosphates (OP), we evaluated the change in removal efficiencies of total organic carbon (TOC) and ammonia, microbial biomass, alkaline phosphatase (AKP) activity, microbial community structure, and AKP-associated genes. TOC and ammonia removal efficiency was highest in IP (64.8%) and cNMP (52.3%) treatments. Microbial community structure changed significantly across phosphorus sources that IP and cNMP encouraged Enterobacter and Aeromonas, respectively. The abundance of phoA and phoU genes was higher in IP treatments, whereas phoD and phoX genes dominated OP treatments. Our findings suggested that the performance of WWTPs was dependent on phosphorus sources and provided new insights into effective WWTP management.

Signaling through non-membrane nuclear phosphoinositide binding proteins in human health and disease

Phosphoinositide membrane signaling is critical for normal physiology, playing well-known roles in diverse human pathologies. The basic mechanisms governing phosphoinositide signaling within the nucleus, however, have remained deeply enigmatic owing to their presence outside the nuclear membranes. Over 40% of nuclear phosphoinositides can exist in this non-membrane state, held soluble in the nucleoplasm by nuclear proteins that remain largely unidentified. Recently, two nuclear proteins responsible for solubilizing phosphoinositides were identified, steroidogenic factor-1 (SF-1; NR5A1) and liver receptor homolog-1 (LRH-1; NR5A2), along with two enzymes that directly remodel these phosphoinositide/protein complexes, phosphatase and tensin homolog (PTEN; MMAC) and inositol polyphosphate multikinase (IPMK; ipk2). These new footholds now permit the assignment of physiological functions for nuclear phosphoinositides in human diseases, such as endometriosis, nonalcoholic fatty liver disease/steatohepatitis, glioblastoma, and hepatocellular carcinoma. The unique nature of nuclear phosphoinositide signaling affords extraordinary clinical opportunities for new biomarkers, diagnostics, and therapeutics. Thus, phosphoinositide biology within the nucleus may represent the next generation of low-hanging fruit for new drugs, not unlike what has occurred for membrane phosphatidylinositol 3-kinase drug development. This review connects recent basic science discoveries in nuclear phosphoinositide signaling to clinical pathologies, with the hope of inspiring development of new therapies.

A probabilistic network model for structural transitions in biomolecules

Biological macromolecules often undergo large conformational rearrangements during a functional cycle. To simulate these structural transitions with full atomic detail typically demands extensive computational resources. Moreover, it is unclear how to incorporate, in a principled way, additional experimental information that could guide the structural transition. This article develops a probabilistic model for conformational transitions in biomolecules. The model can be viewed as a network of anharmonic springs that break, if the experimental data support the rupture of bonds. Hamiltonian Monte Carlo in internal coordinates is used to infer structural transitions from experimental data, thereby sampling large conformational transitions without distorting the structure. The model is benchmarked on a large set of conformational transitions. Moreover, we demonstrate the use of the probabilistic network model for integrative modeling of macromolecular complexes based on data from crosslinking followed by mass spectrometry.

The contribution of modern EPR to structural biology

Electron paramagnetic resonance (EPR) spectroscopy combined with site-directed spin labelling is applicable to biomolecules and their complexes irrespective of system size and in a broad range of environments. Neither short-range nor long-range order is required to obtain structural restraints on accessibility of sites to water or oxygen, on secondary structure, and on distances between sites. Many of the experiments characterize a static ensemble obtained by shock-freezing. Compared with characterizing the dynamic ensemble at ambient temperature, analysis is simplified and information loss due to overlapping timescales of measurement and system dynamics is avoided. The necessity for labelling leads to sparse restraint sets that require integration with data from other methodologies for building models. The double electron–electron resonance experiment provides distance distributions in the nanometre range that carry information not only on the mean conformation but also on the width of the native ensemble. The distribution widths are often inconsistent with Anfinsen's concept that a sequence encodes a single native conformation defined at atomic resolution under physiological conditions.

Natural Conformational Sampling of Human TNFα Visualized by Double Electron-Electron Resonance

Double electron-electron resonance in conjunction with site-directed spin labeling has been used to probe natural conformational sampling of the human tumor necrosis factor α trimer. We suggest a previously unreported, predeoligomerization conformation of the trimer that has been shown to be sampled at low frequency. A model of this trimeric state has been constructed based on crystal structures using the double-electron-electron-resonance distances. The model shows one of the protomers to be rotated and tilted outward at the tip end, leading to a breaking of the trimerous symmetry and distortion at a receptor-binding interface. The new structure offers opportunities to modulate the biological activity of tumor necrosis factor α through stabilization of the distorted trimer with small molecules.

Key features of an Hsp70 chaperone allosteric landscape revealed by ion-mobility native mass spectrometry and double electron-electron resonance

Proteins are dynamic entities that populate conformational ensembles, and most functions of proteins depend on their dynamic character. Allostery, in particular, relies on ligand-modulated shifts in these conformational ensembles. Hsp70s are allosteric molecular chaperones with conformational landscapes that involve large rearrangements of their two domains (viz. the nucleotide-binding domain and substrate-binding domain) in response to adenine nucleotides and substrates. However, it remains unclear how the Hsp70 conformational ensemble is populated at each point of the allosteric cycle and how ligands control these populations. We have mapped the conformational species present under different ligand-binding conditions throughout the allosteric cycle of the Escherichia coli Hsp70 DnaK by two complementary methods, ion-mobility mass spectrometry and double electron-electron resonance. Our results obtained under biologically relevant ligand-bound conditions confirm the current picture derived from NMR and crystallographic data of domain docking upon ATP binding and undocking in response to ADP and substrate. Additionally, we find that the helical lid of DnaK is a highly dynamic unit of the structure in all ligand-bound states. Importantly, we demonstrate that DnaK populates a partially docked state in the presence of ATP and substrate and that this state represents an energy minimum on the DnaK allosteric landscape. Because Hsp70s are emerging as potential drug targets for many diseases, fully mapping an allosteric landscape of a molecular chaperone like DnaK will facilitate the development of small molecules that modulate Hsp70 function via allosteric mechanisms.