Metabolic priming by multiple enzyme systems supports glycolysis, HIF1α stabilisation, and human cancer cell survival in early hypoxia.

Adaptation to chronic hypoxia occurs through changes in protein expression, which are controlled by hypoxia-inducible factor 1α (HIF1α) and are necessary for cancer cell survival. However, the mechanisms that enable cancer cells to adapt in early hypoxia, before the HIF1α-mediated transcription programme is fully established, remain poorly understood. Here we show in human breast cancer cells, that within 3 h of hypoxia exposure, glycolytic flux increases in a HIF1α-independent manner but is limited by NAD+ availability. Glycolytic ATP maintenance and cell survival in early hypoxia rely on reserve lactate dehydrogenase A capacity as well as the activity of glutamate-oxoglutarate transaminase 1 (GOT1), an enzyme that fuels malate dehydrogenase 1 (MDH1)-derived NAD+. In addition, GOT1 maintains low α-ketoglutarate levels, thereby limiting prolyl hydroxylase activity to promote HIF1α stabilisation in early hypoxia and enable robust HIF1α target gene expression in later hypoxia. Our findings reveal that, in normoxia, multiple enzyme systems maintain cells in a primed state ready to support increased glycolysis and HIF1α stabilisation upon oxygen limitation, until other adaptive processes that require more time are fully established.

Erratum: Genome editing reveals a role for OCT4 in human embryogenesis.

This corrects the article DOI: 10.1038/nature24033.

Genome editing reveals a role for OCT4 in human embryogenesis.

Despite their fundamental biological and clinical importance, the molecular mechanisms that regulate the first cell fate decisions in the human embryo are not well understood. Here we use CRISPR-Cas9-mediated genome editing to investigate the function of the pluripotency transcription factor OCT4 during human embryogenesis. We identified an efficient OCT4-targeting guide RNA using an inducible human embryonic stem cell-based system and microinjection of mouse zygotes. Using these refined methods, we efficiently and specifically targeted the gene encoding OCT4 (POU5F1) in diploid human zygotes and found that blastocyst development was compromised. Transcriptomics analysis revealed that, in POU5F1-null cells, gene expression was downregulated not only for extra-embryonic trophectoderm genes, such as CDX2, but also for regulators of the pluripotent epiblast, including NANOG. By contrast, Pou5f1-null mouse embryos maintained the expression of orthologous genes, and blastocyst development was established, but maintenance was compromised. We conclude that CRISPR-Cas9-mediated genome editing is a powerful method for investigating gene function in the context of human development.

In vitro generation of neuromesodermal progenitors reveals distinct roles for wnt signalling in the specification of spinal cord and paraxial mesoderm identity.

Cells of the spinal cord and somites arise from shared, dual-fated precursors, located towards the posterior of the elongating embryo. Here we show that these neuromesodermal progenitors (NMPs) can readily be generated in vitro from mouse and human pluripotent stem cells by activating Wnt and Fgf signalling, timed to emulate in vivo development. Similar to NMPs in vivo, these cells co-express the neural factor Sox2 and the mesodermal factor Brachyury and differentiate into neural and paraxial mesoderm in vitro and in vivo. The neural cells produced by NMPs have spinal cord but not anterior neural identity and can differentiate into spinal cord motor neurons. This is consistent with the shared origin of spinal cord and somites and the distinct ontogeny of the anterior and posterior nervous system. Systematic analysis of the transcriptome during differentiation identifies the molecular correlates of each of the cell identities and the routes by which they are obtained. Moreover, we take advantage of the system to provide evidence that Brachyury represses neural differentiation and that signals from mesoderm are not necessary to induce the posterior identity of spinal cord cells. This indicates that the mesoderm inducing and posteriorising functions of Wnt signalling represent two molecularly separate activities. Together the data illustrate how reverse engineering normal developmental mechanisms allows the differentiation of specific cell types in vitro and the analysis of previous difficult to access aspects of embryo development.

The Phylogenetic Signature Underlying ATP Synthase c-Ring Compliance.

The proton-driven ATP synthase (FOF1) is comprised of two rotary, stepping motors (FO and F1) coupled by an elastic power transmission. The elastic compliance resides in the rotor module that includes the membrane-embedded FO c-ring. Proton transport by FO is firmly coupled to the rotation of the c-ring relative to other FO subunits (ab2). It drives ATP synthesis. We used a computational method to investigate the contribution of the c-ring to the total elastic compliance. We performed principal component analysis of conformational ensembles built using distance constraints from the bovine mitochondrial c-ring x-ray structure. Angular rotary twist, the dominant ring motion, was estimated to show that the c-ring accounted in part for the measured compliance. Ring rotation was entrained to rotation of the external helix within each hairpin-shaped c-subunit in the ring. Ensembles of monomer and dimers extracted from complete c-rings showed that the coupling between collective ring and the individual subunit motions was independent of the size of the c-ring, which varies between organisms. Molecular determinants were identified by covariance analysis of residue coevolution and structural-alphabet-based local dynamics correlations. The residue coevolution gave a readout of subunit architecture. The dynamic couplings revealed that the hinge for both ring and subunit helix rotations was constructed from the proton-binding site and the adjacent glycine motif (IB-GGGG) in the midmembrane plane. IB-GGGG motifs were linked by long-range couplings across the ring, while intrasubunit couplings connected the motif to the conserved cytoplasmic loop and adjacent segments. The correlation with principal collective motions shows that the couplings underlie both ring rotary and bending motions. Noncontact couplings between IB-GGGG motifs matched the coevolution signal as well as contact couplings. The residue coevolution reflects the physiological importance of the dynamics that may link proton transfer to ring compliance.

GSATools: analysis of allosteric communication and functional local motions using a structural alphabet.

MOTIVATION: GSATools is a free software package to analyze conformational ensembles and to detect functional motions in proteins by means of a structural alphabet. The software integrates with the widely used GROMACS simulation package and can generate a range of graphical outputs. Three applications can be supported: (i) investigation of the conformational variability of local structures; (ii) detection of allosteric communication; and (iii) identification of local regions that are critical for global functional motions. These analyses provide insights into the dynamics of proteins and allow for targeted design of functional mutants in theoretical and experimental studies. AVAILABILITY: The C source code of the GSATools, along with a set of pre-compiled binaries, is freely available under GNU General Public License from http://mathbio.nimr.mrc.ac.uk/wiki/GSATools.

Detection of allosteric signal transmission by information-theoretic analysis of protein dynamics.

Allostery offers a highly specific way to modulate protein function. Therefore, understanding this mechanism is of increasing interest for protein science and drug discovery. However, allosteric signal transmission is difficult to detect experimentally and to model because it is often mediated by local structural changes propagating along multiple pathways. To address this, we developed a method to identify communication pathways by an information-theoretical analysis of molecular dynamics simulations. Signal propagation was described as information exchange through a network of correlated local motions, modeled as transitions between canonical states of protein fragments. The method was used to describe allostery in two-component regulatory systems. In particular, the transmission from the allosteric site to the signaling surface of the receiver domain NtrC was shown to be mediated by a layer of hub residues. The location of hubs preferentially connected to the allosteric site was found in close agreement with key residues experimentally identified as involved in the signal transmission. The comparison with the networks of the homologues CheY and FixJ highlighted similarities in their dynamics. In particular, we showed that a preorganized network of fragment connections between the allosteric and functional sites exists already in the inactive state of all three proteins.

A branching model of lineage differentiation underpinning the neurogenic potential of enteric glia.

Glial cells have been proposed as a source of neural progenitors, but the mechanisms underpinning the neurogenic potential of adult glia are not known. Using single cell transcriptomic profiling, we show that enteric glial cells represent a cell state attained by autonomic neural crest cells as they transition along a linear differentiation trajectory that allows them to retain neurogenic potential while acquiring mature glial functions. Key neurogenic loci in early enteric nervous system progenitors remain in open chromatin configuration in mature enteric glia, thus facilitating neuronal differentiation under appropriate conditions. Molecular profiling and gene targeting of enteric glial cells in a cell culture model of enteric neurogenesis and a gut injury model demonstrate that neuronal differentiation of glia is driven by transcriptional programs employed in vivo by early progenitors. Our work provides mechanistic insight into the regulatory landscape underpinning the development of intestinal neural circuits and generates a platform for advancing glial cells as therapeutic agents for the treatment of neural deficits.

STRIKE: evaluation of protein MSAs using a single 3D structure.

MOTIVATION: Evaluating alternative multiple protein sequence alignments is an important unsolved problem in Biology. The most accurate way of doing this is to use structural information. Unfortunately, most methods require at least two structures to be embedded in the alignment, a condition rarely met when dealing with standard datasets. RESULT: We developed STRIKE, a method that determines the relative accuracy of two alternative alignments of the same sequences using a single structure. We validated our methodology on three commonly used reference datasets (BAliBASE, Homestrad and Prefab). Given two alignments, STRIKE manages to identify the most accurate one in 70% of the cases on average. This figure increases to 79% when considering very challenging datasets like the RV11 category of BAliBASE. This discrimination capacity is significantly higher than that reported for other metrics such as Contact Accepted mutation or Blosum. We show that this increased performance results both from a refined definition of the contacts and from the use of an improved contact substitution score. CONTACT: cedric.notredame@crg.eu AVAILABILITY: STRIKE is an open source freeware available from www.tcoffee.org SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Malaria parasite actin polymerization and filament structure.

A novel form of acto-myosin regulation has been proposed in which polymerization of new actin filaments regulates motility of parasites of the apicomplexan class of protozoa. In vivo and in vitro parasite F-actin is very short and unstable, but the structural basis and details of filament dynamics remain unknown. Here, we show that long actin filaments can be obtained by polymerizing unlabeled rabbit skeletal actin (RS-actin) onto both ends of the short rhodamine-phalloidin-stabilized Plasmodium falciparum actin I (Pf-actin) filaments. Following annealing, hybrid filaments of micron length and "zebra-striped" appearance are observed by fluorescence microscopy that are stable enough to move over myosin class II motors in a gliding filament assay. Using negative stain electron microscopy we find that pure Pf-actin stabilized by jasplakinolide (JAS) also forms long filaments, indistinguishable in length from RS-actin filaments, and long enough to be characterized structurally. To compare structures in near physiological conditions in aqueous solution we imaged Pf-actin and RS-actin filaments by atomic force microscopy (AFM). We found the monomer stacking to be distinctly different for Pf-actin compared with RS-actin, such that the pitch of the double helix of Pf-actin filaments was 10% larger. Our results can be explained by a rotational angle between subunits that is larger in the parasite compared with RS-actin. Modeling of the AFM data using high-resolution actin filament models supports our interpretation of the data. The structural differences reported here may be a consequence of weaker inter- and intra-strand contacts, and may be critical for differences in filament dynamics and for regulation of parasite motility.

Structural alphabets derived from attractors in conformational space.

BACKGROUND: The hierarchical and partially redundant nature of protein structures justifies the definition of frequently occurring conformations of short fragments as 'states'. Collections of selected representatives for these states define Structural Alphabets, describing the most typical local conformations within protein structures. These alphabets form a bridge between the string-oriented methods of sequence analysis and the coordinate-oriented methods of protein structure analysis. RESULTS: A Structural Alphabet has been derived by clustering all four-residue fragments of a high-resolution subset of the protein data bank and extracting the high-density states as representative conformational states. Each fragment is uniquely defined by a set of three independent angles corresponding to its degrees of freedom, capturing in simple and intuitive terms the properties of the conformational space. The fragments of the Structural Alphabet are equivalent to the conformational attractors and therefore yield a most informative encoding of proteins. Proteins can be reconstructed within the experimental uncertainty in structure determination and ensembles of structures can be encoded with accuracy and robustness. CONCLUSIONS: The density-based Structural Alphabet provides a novel tool to describe local conformations and it is specifically suitable for application in studies of protein dynamics.

Functional cross-talk between allosteric effects of activating and inhibiting ligands underlies PKM2 regulation.

Several enzymes can simultaneously interact with multiple intracellular metabolites, however, how the allosteric effects of distinct ligands are integrated to coordinately control enzymatic activity remains poorly understood. We addressed this question using, as a model system, the glycolytic enzyme pyruvate kinase M2 (PKM2). We show that the PKM2 activator fructose 1,6-bisphosphate (FBP) alone promotes tetramerisation and increases PKM2 activity, but addition of the inhibitor L-phenylalanine (Phe) prevents maximal activation of FBP-bound PKM2 tetramers. We developed a method, AlloHubMat, that uses eigenvalue decomposition of mutual information derived from molecular dynamics trajectories to identify residues that mediate FBP-induced allostery. Experimental mutagenesis of these residues identified PKM2 variants in which activation by FBP remains intact but cannot be attenuated by Phe. Our findings reveal residues involved in FBP-induced allostery that enable the integration of allosteric input from Phe and provide a paradigm for the coordinate regulation of enzymatic activity by simultaneous allosteric inputs.

MinSet: a general approach to derive maximally representative database subsets by using fragment dictionaries and its application to the SCOP database.

MOTIVATION: The size of current protein databases is a challenge for many Bioinformatics applications, both in terms of processing speed and information redundancy. It may be therefore desirable to efficiently reduce the database of interest to a maximally representative subset. RESULTS: The MinSet method employs a combination of a Suffix Tree and a Genetic Algorithm for the generation, selection and assessment of database subsets. The approach is generally applicable to any type of string-encoded data, allowing for a drastic reduction of the database size whilst retaining most of the information contained in the original set. We demonstrate the performance of the method on a database of protein domain structures encoded as strings. We used the SCOP40 domain database by translating protein structures into character strings by means of a structural alphabet and by extracting optimized subsets according to an entropy score that is based on a constant-length fragment dictionary. Therefore, optimized subsets are maximally representative for the distribution and range of local structures. Subsets containing only 10% of the SCOP structure classes show a coverage of >90% for fragments of length 1-4. AVAILABILITY: http://mathbio.nimr.mrc.ac.uk/~jkleinj/MinSet. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The TAB1-p38α complex aggravates myocardial injury and can be targeted by small molecules.

Inhibiting MAPK14 (p38α) diminishes cardiac damage in myocardial ischemia. During myocardial ischemia, p38α interacts with TAB1, a scaffold protein, which promotes p38α autoactivation; active p38α (pp38α) then transphosphorylates TAB1. Previously, we solved the X-ray structure of the p38α-TAB1 (residues 384-412) complex. Here, we further characterize the interaction by solving the structure of the pp38α-TAB1 (residues 1-438) complex in the active state. Based on this information, we created a global knock-in (KI) mouse with substitution of 4 residues on TAB1 that we show are required for docking onto p38α. Whereas ablating p38α or TAB1 resulted in early embryonal lethality, the TAB1-KI mice were viable and had no appreciable alteration in their lymphocyte repertoire or myocardial transcriptional profile; nonetheless, following in vivo regional myocardial ischemia, infarction volume was significantly reduced and the transphosphorylation of TAB1 was disabled. Unexpectedly, the activation of myocardial p38α during ischemia was only mildly attenuated in TAB1-KI hearts. We also identified a group of fragments able to disrupt the interaction between p38α and TAB1. We conclude that the interaction between the 2 proteins can be targeted with small molecules. The data reveal that it is possible to selectively inhibit signaling downstream of p38α to attenuate ischemic injury.

Interaction of malaria parasite-inhibitory antibodies with the merozoite surface protein MSP1(19) by computational docking.

Merozoite surface protein 1 (MSP1) of the malaria parasite Plasmodium falciparum is an important vaccine candidate antigen. Antibodies specific for the C-terminal maturation product, MSP1(19), have been shown to inhibit erythrocyte invasion and parasite growth. Specific monoclonal antibodies react with conformational epitopes contained within the two EGF-like domains that constitute the antigen MSP1(19). To gain greater insight into the inhibitory process, the authors selected two strongly inhibitory antibodies (designated 12.8 and 12.10) and modeled their structures by homology. Computational docking was used to generate antigen-antibody complexes and a selection filter based on NMR data was applied to obtain plausible models. Molecular Dynamics simulations of the selected complexes were performed to evaluate the role of specific side chains in the binding. Favorable complexes were obtained that complement the NMR data in defining specific binding sites. These models can provide valuable guidelines for future experimental work that is devoted to the understanding of the action mechanism of invasion-inhibitory antibodies.

POPSCOMP: an automated interaction analysis of biomolecular complexes.

Large-scale analysis of biomolecular complexes reveals the functional network within the cell. Computational methods are required to extract the essential information from the available data. The POPSCOMP server is designed to calculate the interaction surface between all components of a given complex structure consisting of proteins, DNA or RNA molecules. The server returns matrices and graphs of surface area burial that can be used to automatically annotate components and residues that are involved in complex formation, to pinpoint conformational changes and to estimate molecular interaction energies. The analysis can be performed on a per-atom level or alternatively on a per-residue level for low-resolution structures. Here, we present an analysis of ribosomal structures in complex with various antibiotics to exemplify the potential and limitations of automated complex analysis. The POPSCOMP server is accessible at http://ibivu.cs.vu.nl/programs/popscompwww/.

A Gene Regulatory Network Balances Neural and Mesoderm Specification during Vertebrate Trunk Development.

Transcriptional networks, regulated by extracellular signals, control cell fate decisions and determine the size and composition of developing tissues. One example is the network controlling bipotent neuromesodermal progenitors (NMPs) that fuel embryo elongation by generating spinal cord and trunk mesoderm tissue. Here, we use single-cell transcriptomics to identify the molecular signature of NMPs and reverse engineer the mechanism that regulates their differentiation. Together with genetic perturbations, this reveals a transcriptional network that integrates opposing retinoic acid (RA) and Wnt signals to determine the rate at which cells enter and exit the NMP state. RA, produced by newly generated mesodermal cells, provides feedback that initiates NMP generation and induces neural differentiation, thereby coordinating the production of neural and mesodermal tissue. Together, the data define a regulatory network architecture that balances the generation of different cell types from bipotential progenitors in order to facilitate orderly axis elongation.

Lineage-dependent spatial and functional organization of the mammalian enteric nervous system.

The enteric nervous system (ENS) is essential for digestive function and gut homeostasis. Here we show that the amorphous neuroglia networks of the mouse ENS are composed of overlapping clonal units founded by postmigratory neural crest-derived progenitors. The spatial configuration of ENS clones depends on proliferation-driven local interactions of ENS progenitors with lineally unrelated neuroectodermal cells, the ordered colonization of the serosa-mucosa axis by clonal descendants, and gut expansion. Single-cell transcriptomics and mutagenesis analysis delineated dynamic molecular states of ENS progenitors and identified RET as a regulator of neurogenic commitment. Clonally related enteric neurons exhibit synchronous activity in response to network stimulation. Thus, lineage relationships underpin the organization of the peripheral nervous system.

An engineered photoswitchable mammalian pyruvate kinase.

Changes in allosteric regulation of glycolytic enzymes have been linked to metabolic reprogramming involved in cancer. Remarkably, allosteric mechanisms control enzyme function at significantly shorter time-scales compared to the long-term effects of metabolic reprogramming on cell proliferation. It remains unclear if and how the speed and reversibility afforded by rapid allosteric control of metabolic enzymes is important for cell proliferation. Tools that allow specific, dynamic modulation of enzymatic activities in mammalian cells would help address this question. Towards this goal, we have used molecular dynamics simulations to guide the design of mPKM2 internal light/oxygen/voltage-sensitive domain 2 (LOV2) fusion at position D24 (PiL[D24]), an engineered pyruvate kinase M2 (PKM2) variant that harbours an insertion of the light-sensing LOV2 domain from Avena Sativa within a region implicated in allosteric regulation by fructose 1,6-bisphosphate (FBP). The LOV2 photoreaction is preserved in the PiL[D24] chimera and causes secondary structure changes that are associated with a 30% decrease in the Km of the enzyme for phosphoenolpyruvate resulting in increased pyruvate kinase activity after light exposure. Importantly, this change in activity is reversible upon light withdrawal. Expression of PiL[D24] in cells leads to light-induced increase in labelling of pyruvate from glucose. PiL[D24] therefore could provide a means to modulate cellular glucose metabolism in a remote manner and paves the way for studying the importance of rapid allosteric phenomena in the regulation of metabolism and enzyme control.

POPS: A fast algorithm for solvent accessible surface areas at atomic and residue level.

POPS (Parameter OPtimsed Surfaces) is a new method to calculate solvent accessible surface areas, which is based on an empirically parameterisable analytical formula and fast to compute. Atomic and residue areas (the latter represented by a single sphere centered on the C(alpha) atom of amino acids and at the P atom of nucleotides) have been optimised versus accurate all-atom methods. The parameterisation has been derived from a selected dataset of proteins and nucleic acids of different sizes and topologies. The residue based approach POPS-R, has been devised as a useful tool for the analysis of large macromolecular assemblies like the ribosome and it is specially suited for the refinement of low resolution structures. POPS-R also allows for estimates of the loss of free energy of solvation upon complex formation, which should be particularly useful for the design of new protein-protein and protein-nucleic acid complexes. The program POPS is available at http://mathbio.nimr.mrc.ac.uk/~ffranca/POPS and at the mirror site http://www.cs.vu.nl/~ibivu/programs/popswww.