Search

MolecularWebXR is a new web-based platform for education, science communication and scientific peer discussion in chemistry and biology, based on modern web-based Virtual Reality (VR) and Augmented Reality (AR). With no installs as it is all web-served, MolecularWebXR enables multiple users to simultaneously explore, communicate and discuss concepts about chemistry and biology in immersive 3D environments, by manipulating and passing around objects with their bare hands and pointing at different elements with natural hand gestures. Users may either be present in the same physical space or distributed around the world, in the latter case talking naturally with each other thanks to built-in audio. While MolecularWebXR offers the most immersive experience on high-end AR/VR headsets, its WebXR core also supports participation on consumer devices such as smartphones (with optional cardboard goggles for enhanced immersion), computers, and tablets. MolecularWebXR includes preset VR rooms covering topics in general, inorganic, and organic chemistry, as well as biophysics, structural biology, and general biology. Users can also add new content via the PDB2AR tool. We demonstrate MolecularWebXR's versatility and ease of use across a wide age range (12-80) in fully virtual and mixed real-virtual sessions at science outreach events, undergraduate and graduate courses, scientific collaborations, and conference presentations. MolecularWebXR is available for free use without registration at https://molecularwebxr.org. A blog post version of this preprint with embedded videos is available at https://go.epfl.ch/molecularwebxr-blog-post., Competing Interests: Declaration of competing interest The authors declare no conflict of interest., (Copyright © 2024. Published by Elsevier Inc.)

Protein design and engineering are evolving at an unprecedented pace leveraging the advances in deep learning. Current models nonetheless cannot natively consider non-protein entities within the design process. Here, we introduce a deep learning approach based solely on a geometric transformer of atomic coordinates and element names that predicts protein sequences from backbone scaffolds aware of the restraints imposed by diverse molecular environments. To validate the method, we show that it can produce highly thermostable, catalytically active enzymes with high success rates. This concept is anticipated to improve the versatility of protein design pipelines for crafting desired functions., (© 2024. The Author(s).)

Cardiolipin (CL) is a mitochondria-specific phospholipid that forms heterotypic interactions with membrane-shaping proteins and regulates the dynamic remodeling and function of mitochondria. However, the precise mechanisms through which CL influences mitochondrial morphology are not well understood. In this study, employing molecular dynamics (MD) simulations, we observed CL localize near the membrane-binding sites of the mitochondrial fusion protein Optic Atrophy 1 (OPA1). To validate these findings experimentally, we developed a bromine-labeled CL probe to enhance cryoEM contrast and characterize the structure of OPA1 assemblies bound to the CL-brominated lipid bilayers. Our images provide direct evidence of interactions between CL and two conserved motifs within the paddle domain (PD) of OPA1, which control membrane-shaping mechanisms. We further observed a decrease in membrane remodeling activity for OPA1 in lipid compositions with increasing concentrations of monolyso-cardiolipin (MLCL). Suggesting that the partial replacement of CL by MLCL accumulation, as observed in Barth syndrome-associated mutations of the tafazzin phospholipid transacylase, compromises the stability of protein-membrane interactions. Our analyses provide insights into how biological membranes regulate the mechanisms governing mitochondrial homeostasis., Competing Interests: Competing Interest Statement A.F. and F.R.M. are shareholders and employees of Altos Labs.

Distinct morphologies of the mitochondrial network support divergent metabolic and regulatory processes that determine cell function and fate 1-3 . The mechanochemical GTPase optic atrophy 1 (OPA1) influences the architecture of cristae and catalyses the fusion of the mitochondrial inner membrane 4,5 . Despite its fundamental importance, the molecular mechanisms by which OPA1 modulates mitochondrial morphology are unclear. Here, using a combination of cellular and structural analyses, we illuminate the molecular mechanisms that are key to OPA1-dependent membrane remodelling and fusion. Human OPA1 embeds itself into cardiolipin-containing membranes through a lipid-binding paddle domain. A conserved loop within the paddle domain inserts deeply into the bilayer, further stabilizing the interactions with cardiolipin-enriched membranes. OPA1 dimerization through the paddle domain promotes the helical assembly of a flexible OPA1 lattice on the membrane, which drives mitochondrial fusion in cells. Moreover, the membrane-bending OPA1 oligomer undergoes conformational changes that pull the membrane-inserting loop out of the outer leaflet and contribute to the mechanics of membrane remodelling. Our findings provide a structural framework for understanding how human OPA1 shapes mitochondrial morphology and show us how human disease mutations compromise OPA1 functions., (© 2023. The Author(s), under exclusive licence to Springer Nature Limited.)

Neural differentiation, synaptic transmission, and action potential propagation depend on membrane sphingolipids, whose metabolism is tightly regulated. Mutations in the ceramide transporter CERT (CERT1), which is involved in sphingolipid biosynthesis, are associated with intellectual disability, but the pathogenic mechanism remains obscure. Here, we characterize 31 individuals with de novo missense variants in CERT1. Several variants fall into a previously uncharacterized dimeric helical domain that enables CERT homeostatic inactivation, without which sphingolipid production goes unchecked. The clinical severity reflects the degree to which CERT autoregulation is disrupted, and inhibiting CERT pharmacologically corrects morphological and motor abnormalities in a Drosophila model of the disease, which we call ceramide transporter (CerTra) syndrome. These findings uncover a central role for CERT autoregulation in the control of sphingolipid biosynthetic flux, provide unexpected insight into the structural organization of CERT, and suggest a possible therapeutic approach for patients with CerTra syndrome.

Proteins are essential molecular building blocks of life, responsible for most biological functions as a result of their specific molecular interactions. However, predicting their  binding  interfaces remains a challenge. In this study, we present a geometric transformer that acts directly on atomic coordinates labeled only with element names. The resulting model-the Protein Structure Transformer, PeSTo-surpasses the current state of the art in predicting protein-protein interfaces and can also predict and differentiate between interfaces involving nucleic acids, lipids, ions, and small molecules with high confidence. Its low computational cost enables processing high volumes of structural data, such as molecular dynamics ensembles allowing for the discovery of interfaces that remain otherwise inconspicuous in static experimentally solved structures. Moreover, the growing foldome provided by de novo structural predictions can be easily analyzed, providing new opportunities to uncover unexplored biology., (© 2023. The Author(s).)

The complex architecture of the endoplasmic reticulum (ER) comprises distinct dynamic features, many at the nanoscale, that enable the coexistence of the nuclear envelope, regions of dense sheets and a branched tubular network that spans the cytoplasm. A key player in the formation of ER sheets is cytoskeleton-linking membrane protein 63 (CLIMP-63). The mechanisms by which CLIMP-63 coordinates ER structure remain elusive. Here, we address the impact of S-acylation, a reversible post-translational lipid modification, on CLIMP-63 cellular distribution and function. Combining native mass-spectrometry, with kinetic analysis of acylation and deacylation, and data-driven mathematical modelling, we obtain in-depth understanding of the CLIMP-63 life cycle. In the ER, it assembles into trimeric units. These occasionally exit the ER to reach the plasma membrane. However, the majority undergoes S-acylation by ZDHHC6 in the ER where they further assemble into highly stable super-complexes. Using super-resolution microscopy and focused ion beam electron microscopy, we show that CLIMP-63 acylation-deacylation controls the abundance and fenestration of ER sheets. Overall, this study uncovers a dynamic lipid post-translational regulation of ER architecture., (© 2023. The Author(s).)

Coenzyme Q (CoQ) is a redox-active lipid essential for core metabolic pathways and antioxidant defense. CoQ is synthesized upon the mitochondrial inner membrane by an ill-defined "complex Q" metabolon. Here, we present structure-function analyses of a lipid-, substrate-, and NADH-bound complex comprising two complex Q subunits: the hydroxylase COQ7 and the lipid-binding protein COQ9. We reveal that COQ7 adopts a ferritin-like fold with a hydrophobic channel whose substrate-binding capacity is enhanced by COQ9. Using molecular dynamics, we further show that two COQ7:COQ9 heterodimers form a curved tetramer that deforms the membrane, potentially opening a pathway for the CoQ intermediates to translocate from the bilayer to the proteins' lipid-binding sites. Two such tetramers assemble into a soluble octamer with a pseudo-bilayer of lipids captured within. Together, these observations indicate that COQ7 and COQ9 cooperate to access hydrophobic precursors within the membrane and coordinate subsequent synthesis steps toward producing CoQ., Competing Interests: Declaration of interests J.J.C. is a consultant for Thermo Scientific. A.F. is a shareholder and employee of Altos Labs and a shareholder and consultant for Relay Therapeutics, LLC., (Copyright © 2022 Elsevier Inc. All rights reserved.)

Several groups developed in the last years augmented and virtual reality (AR/VR) software to visualize 3D molecules, most rather static, limited in content, and requiring software installs, some even requiring expensive hardware. We launched in 2020 moleculARweb (https://molecularweb.epfl.ch), a website that offers interactive content for chemistry and structural biology education through commodity web-based AR that works on consumer devices like smartphones, tablets and laptops. Among thousands of users, teachers increasingly request more biological macromolecules to be available, a demand that we cannot address individually. Therefore, to allow users to build their own material, we built a web interface where they can create online AR experiences in few steps starting from Protein Data Bank, AlphaFold or custom uploaded structures, or from virtual objects/scenes exported from the Visual Molecular Dynamics program, without any programming knowledge. The web tool also returns WebXR sessions for viewing and manipulating the models in WebXR-compatible devices including smartphones, tablets, and also immersive VR headsets with WebXR-capable browsers, where models can be manipulated even with bare hands when supported by the device. The tool is accessible for free at https://molecularweb.epfl.ch/pages/pdb2ar.html., (Copyright © 2022 The Authors. Published by Elsevier Inc. All rights reserved.)

moleculARweb (https://molecularweb.epfl.ch) began as a website for education and outreach in chemistry and structural biology through augmented reality (AR) content that runs in the web browsers of regular devices like smartphones, tablets, and computers. Here we present two evolutions of moleculARweb's Virtual Modeling Kits (VMK), tools where users can build and view molecules, and explore their mechanics, in 3D AR by handling the molecules in full 3D with custom-printed cube markers (VMK 2.0) or by moving around a simulated scene with mouse or touch gestures (VMK 3.0). Upon simulation the molecules experience visually realistic torsions, clashes, and hydrogen-bonding interactions that the user can manually switch on and off to explore their effects. Moreover, by manually tuning a fictitious temperature the users can accelerate conformational transitions or 'freeze' specific conformations for careful inspection in 3D. Even some phase transitions and separations can be simulated. We here showcase these and other features of the new VMKs connecting them to possible specific applications to teaching and self-learning of concepts from general, organic, biological and physical chemistry; and in assisting with small tasks in molecular modelling for research. Last, in a short discussion section we overview what future developments are needed for the 'dream tool' for the future of chemistry education and work., (Copyright 2022 Fabio Cortes Rodriguez, Lucien Krapp, Matteo Dal Peraro, Luciano Abriata. License: This work is licensed under a Creative Commons Attribution 4.0 International License.)

SARS-CoV-2 virions are surrounded by a lipid bilayer that contains membrane proteins such as spike, responsible for target-cell binding and virus fusion. We found that during SARS-CoV-2 infection, spike becomes lipid modified, through the sequential action of the S-acyltransferases ZDHHC20 and 9. Particularly striking is the rapid acylation of spike on 10 cytosolic cysteines within the ER and Golgi. Using a combination of computational, lipidomics, and biochemical approaches, we show that this massive lipidation controls spike biogenesis and degradation, and drives the formation of localized ordered cholesterol and sphingolipid-rich lipid nanodomains in the early Golgi, where viral budding occurs. Finally, S-acylation of spike allows the formation of viruses with enhanced fusion capacity. Our study points toward S-acylating enzymes and lipid biosynthesis enzymes as novel therapeutic anti-viral targets., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2021 The Authors. Published by Elsevier Inc. All rights reserved.)

Post-translational modifications (PTMs) within the first 17 amino acids (Nt17) of the Huntingtin protein (Htt) have been shown to inhibit the aggregation and attenuate the toxicity of mutant Htt proteins in vitro and in various models of Huntington's disease. Here, we expand on these studies by investigating the effect of methionine eight oxidation (oxM8) and its crosstalk with lysine 6 acetylation (AcK6) or threonine 3 phosphorylation (pT3) on the aggregation of mutant Httex1 (mHttex1). We show that M8 oxidation delays but does not inhibit the aggregation and has no effect on the final morphologies of mHttex1aggregates. The presence of both oxM8 and AcK6 resulted in dramatic inhibition of Httex1 fibrillization. Circular dichroism spectroscopy and molecular dynamics simulation studies show that PTMs that lower the mHttex1 aggregation rate (oxM8, AcK6/oxM8, pT3, pT3/oxM8, and pS13) result in increased population of a short N-terminal helix (first eight residues) in Nt17 or decreased abundance of other helical forms, including long helix and short C-terminal helix. PTMs that did not alter the aggregation rate (AcK6) of mHttex1 exhibit a similar distribution of helical conformation as the unmodified peptides. These results show that the relative abundance of N- vs. C-terminal helical conformations and long helices, rather than the overall helicity of Nt17, better explains the effect of different Nt17 PTMs on mHttex1; thus, explaining the lack of correlation between the effect of PTMs on the overall helicity of Nt17 and mHttex1 aggregation in vitro . Taken together, our results provide novel structural insight into the differential effects of single PTMs and crosstalk between different PTMs in regulating mHttex1 aggregation., Competing Interests: HL was the Founder and Chief Scientific Office of ND BioSciences SA. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Chiki, Zhang, Rajasekhar, Abriata, Rostami, Krapp, Boudeffa, Dal Peraro and Lashuel.)

Glutamate dehydrogenases (GDHs) are widespread metabolic enzymes that play key roles in nitrogen homeostasis. Large glutamate dehydrogenases composed of 180 kDa subunits (L-GDHs 180 ) contain long N- and C-terminal segments flanking the catalytic core. Despite the relevance of L-GDHs 180 in bacterial physiology, the lack of structural data for these enzymes has limited the progress of functional studies. Here we show that the mycobacterial L-GDH 180 (mL-GDH 180 ) adopts a quaternary structure that is radically different from that of related low molecular weight enzymes. Intersubunit contacts in mL-GDH 180 involve a C-terminal domain that we propose as a new fold and a flexible N-terminal segment comprising ACT-like and PAS-type domains that could act as metabolic sensors for allosteric regulation. These findings uncover unique aspects of the structure-function relationship in the subfamily of L-GDHs.

Residue coevolution estimations coupled to machine learning methods are revolutionizing the ability of protein structure prediction approaches to model proteins that lack clear homologous templates in the Protein Data Bank (PDB). This has been patent in the last round of the Critical Assessment of Structure Prediction (CASP), which presented several very good models for the hardest targets. Unfortunately, literature reporting on these advances often lacks digests tailored to lay end users; moreover, some of the top-ranking predictors do not provide webservers that can be used by nonexperts. How can then end users benefit from these advances and correctly interpret the predicted models? Here we review the web resources that biologists can use today to take advantage of these state-of-the-art methods in their research, including not only the best de novo modeling servers but also datasets of models precomputed by experts for structurally uncharacterized protein families. We highlight their features, advantages and pitfalls for predicting structures of proteins without clear templates. We present a broad number of applications that span from driving forward biochemical investigations that lack experimental structures to actually assisting experimental structure determination in X-ray diffraction, cryo-EM and other forms of integrative modeling. We also discuss issues that must be considered by users yet still require further developments, such as global and residue-wise model quality estimates and sources of residue coevolution other than monomeric tertiary structure., (© The Author(s) 2020. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

Developing molecules that emulate the properties of naturally occurring ice-binding proteins (IBPs) is a daunting challenge. Rather than relying on the (limited) existing structure-property relationships that have been established for IBPs, here we report the use of phage display for the identification of short peptide mimics of IBPs. To this end, an ice-affinity selection protocol is developed, which enables the selection of a cyclic ice-binding peptide containing just 14 amino acids. Mutational analysis identifies three residues, Asp8, Thr10 and Thr14, which are found to be essential for ice binding. Molecular dynamics simulations reveal that the side chain of Thr10 hydrophobically binds to ice revealing a potential mechanism. To demonstrate the biotechnological potential of this peptide, it is expressed as a fusion ('Ice-Tag') with mCherry and used to purify proteins directly from cell lysate.

Continuous assessment of transferable forcefields for molecular simulations is essential to identify their weaknesses and direct improvement efforts. The latest efforts focused on better describing disordered proteins while retaining proper description of folded domains, important because forcefields of the previous generations produce overly compact disordered states. Such improvements should additionally alleviate the related problem of over-stabilized protein-protein interactions, which has been largely overlooked. Here we evaluated three state-of-the-art forcefields, current flagships of their respective developers, optimized for ordered and disordered proteins: CHARMM36m with its recommended corrected TIP3P* water, ff19SB with the recommended OPC water, and the 2019 a99SBdisp forcefield by D. E. Shaw Research with its modified TIP4P water; plus ff14SB with TIP3P as an example of the former generation of forcefields. Our evaluation entailed simulations of (i) multiple copies of a protein that is highly soluble yet undergoes weak dimerization, (ii) a disordered peptide with low, well-characterized alpha helical propensity, and (iii) a peptide known to form insoluble β-aggregates. Our results recapitulate ff14SB-TIP3P over-stabilizing aggregates and secondary structures and place a99SBdisp-TIP4PD at the other end i.e. predicting overly weak intermolecular interactions despite reasonably predicting secondary structure propensities. In-between, CHARMM36m-TIP3P* still over-stabilizes aggregates but predicts residue-wise alpha helical propensities in solution slightly better than ff19SB-OPC, while ff19SB-OPC poses the best prediction of weak dimerization of the soluble protein still predicting aggregation of the β-peptides. This independent assessment shows that the claimed forcefield improvements are real, but also that a right balance between noncovalent attraction and repulsion has not yet been reached. We thus propose developers to consider systems like those tested here in their forcefield tuning protocols. Last, the good performance of CHARMM36m-TIP3P* further shows that tuning 3-point water models might still be an alternative to the more costly 4-point models like OPC and TIP4PD., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2021 The Author(s).)

Phosphorus is an essential nutrient taken up by organisms in the form of inorganic phosphate (Pi). Eukaryotes have evolved sophisticated Pi sensing and signaling cascades, enabling them to stably maintain cellular Pi concentrations. Pi homeostasis is regulated by inositol pyrophosphate signaling molecules (PP-InsPs), which are sensed by SPX domain-containing proteins. In plants, PP-InsP-bound SPX receptors inactivate Myb coiled-coil (MYB-CC) Pi starvation response transcription factors (PHRs) by an unknown mechanism. Here we report that a InsP 8 -SPX complex targets the plant-unique CC domain of PHRs. Crystal structures of the CC domain reveal an unusual four-stranded anti-parallel arrangement. Interface mutations in the CC domain yield monomeric PHR1, which is no longer able to bind DNA with high affinity. Mutation of conserved basic residues located at the surface of the CC domain disrupt interaction with the SPX receptor in vitro and in planta, resulting in constitutive Pi starvation responses. Together, our findings suggest that InsP 8 regulates plant Pi homeostasis by controlling the oligomeric state and hence the promoter binding capability of PHRs via their SPX receptors.

Capillary morphogenesis gene 2 (CMG2/ANTXR2) is a cell surface receptor for both collagen VI and anthrax toxin. Biallelic loss-of-function mutations in CMG2 lead to a severe condition, hyaline fibromatosis syndrome (HFS). We have here dissected a network of dynamic interactions between CMG2 and various actin interactors and regulators, describing a different behavior from other extracellular matrix receptors. CMG2 binds talin, and thereby the actin cytoskeleton, only in its ligand-free state. Extracellular ligand binding leads to src-dependent talin release and recruitment of the actin cytoskeleton regulator RhoA and its effectors. These sequential interactions of CMG2 are necessary for the control of oriented cell division during fish development. Finally, we demonstrate that effective switching between talin and RhoA binding is required for the intracellular degradation of collagen VI in human fibroblasts, which explains why HFS mutations in the cytoskeleton-binding domain lead to dysregulation of extracellular matrix homeostasis., Competing Interests: Declaration of Interests The authors declare no competing interests., (Copyright © 2020 Elsevier Inc. All rights reserved.)

For years, immersive interfaces using virtual and augmented reality (AR) for molecular visualization and modeling have promised a revolution in the way how we teach, learn, communicate and work in chemistry, structural biology and related areas. However, most tools available today for immersive modeling require specialized hardware and software, and are costly and cumbersome to set up. These limitations prevent wide use of immersive technologies in education and research centers in a standardized form, which in turn prevents large-scale testing of the actual effects of such technologies on learning and thinking processes. Here, I discuss building blocks for creating marker-based AR applications that run as web pages on regular computers, and explore how they can be exploited to develop web content for handling virtual molecular systems in commodity AR with no more than a webcam- and internet-enabled computer. Examples span from displaying molecules, electron microscopy maps and molecular orbitals with minimal amounts of HTML code, to incorporation of molecular mechanics, real-time estimation of experimental observables and other interactive resources using JavaScript. These web apps provide virtual alternatives to physical, plastic-made molecular modeling kits, where the computer augments the experience with information about spatial interactions, reactivity, energetics, etc. The ideas and prototypes introduced here should serve as starting points for building active content that everybody can utilize online at minimal cost, providing novel interactive pedagogic material in such an open way that it could enable mass-testing of the effect of immersive technologies on chemistry education., Competing Interests: The author declares that he has no competing interests., (© 2020 Abriata.)

The treatment of hospital- and community-associated infections by methicillin-resistant Staphylococcus aureus (MRSA) is a perpetual challenge. This Gram-positive bacterium is resistant specifically to β-lactam antibiotics, and generally to many other antibacterial agents. Its resistance mechanisms to β-lactam antibiotics are activated only when the bacterium encounters a β-lactam. This activation is regulated by the transmembrane sensor/signal transducer proteins BlaR1 and MecR1. Neither the transmembrane/metalloprotease domain, nor the complete MecR1 and BlaR1 proteins, are isolatable for mechanistic study. Here we propose a model for full-length MecR1 based on homology modeling, residue coevolution data, a new extensive experimental mapping of transmembrane topology, partial structures, molecular simulations, and available NMR data. Our model defines the metalloprotease domain as a hydrophilic transmembrane chamber effectively sealed by the apo-sensor domain. It proposes that the amphipathic helices inserted into the gluzincin domain constitute the route for transmission of the β-lactam-binding event in the extracellular sensor domain, to the intracellular and membrane-embedded zinc-containing active site. From here, we discuss possible routes for subsequent activation of proteolytic action. This study provides the first coherent model of the structure of MecR1, opening routes for future functional investigations on how β-lactam binding culminates in the proteolytic degradation of MecI.

Environmental awareness is an essential attribute of all organisms. The homeoviscous adaptation system of Bacillus subtilis provides a powerful experimental model for the investigation of stimulus detection and signaling mechanisms at the molecular level. These bacteria sense the order of membrane lipids with the transmembrane (TM) protein DesK, which has an N-terminal sensor domain and an intracellular catalytic effector domain. DesK exhibits autokinase activity as well as phosphotransferase and phosphatase activities toward a cognate response regulator, DesR, that controls the expression of an enzyme that remodels membrane fluidity when the temperature drops below ∼30°C. Membrane fluidity signals are transmitted from the DesK sensor domain to the effector domain via rotational movements of a connecting 2-helix coiled coil (2-HCC). Previous molecular dynamic simulations suggested important roles for TM prolines in transducing the initial signals of membrane fluidity status to the 2-HCC. Here, we report that individual replacement of prolines in DesKs TM1 and TM5 helices by alanine (DesKPA) locked DesK in a phosphatase-ON state, abrogating membrane fluidity responses. An unbiased mutagenic screen identified the L174P replacement in the internal side of the repeated heptad of the 2-HCC structure that alleviated the signaling defects of every transmembrane DesKPA substitution. Moreover, substitutions by proline in other internal positions of the 2-HCC reestablished the kinase-ON state of the DesKPA mutants. These results imply that TM prolines are essential for finely tuned signal generation by the N-terminal sensor helices, facilitating a conformational control by the metastable 2-HCC domain of the DesK signaling state. IMPORTANCE Signal sensing and transduction is an essential biological process for cell adaptation and survival. Histidine kinases (HK) are the sensory proteins of two-component systems that control many bacterial responses to different stimuli, like environmental changes. Here, we focused on the HK DesK from Bacillus subtilis , a paradigmatic example of a transmembrane thermosensor suited to remodel membrane fluidity when the temperature drops below 30°C. DesK provides a tractable system for investigating the mechanism of transmembrane signaling, one of the majors interrogates in biology to date. Our studies demonstrate that transmembrane proline residues modulate the conformational switch of a 2-helix coiled-coil (2-HCC) structural motif that controls input-output in a variety of HK. Our results highlight the relevance of proline residues within sensor domains and could inspire investigations of their role in different signaling proteins., (Copyright © 2019 Fernández et al.)

A recent analysis of evolutionary rates in >500 globular soluble enzymes revealed pervasive conservation gradients toward catalytic residues. By looking at amino acid preference profiles rather than evolutionary rates in the same data set, we quantified the effects of active sites on site-specific constraints for physicochemical traits. We found that conservation gradients respond to constraints for polarity, hydrophobicity, flexibility, rigidity and structure in ways consistent with fold polarity principles; while sites far from active sites seem to experience no physicochemical constraint, rather being highly variable and favoring amino acids of low metabolic cost. Globally, our results highlight that amino acid variation contains finer information about protein structure than usually regarded in evolutionary models, and that this information is retrievable automatically with simple fits. We propose that analyses of the kind presented here incorporated into models of protein evolution should allow for better description of the physical chemistry that underlies molecular evolution., (© The Author(s) 2019. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Stenotrophomonas maltophilia is a Gram-negative, nonfermenting, environmental bacillus that is an important cause of nosocomial infections, primarily associated with the respiratory tract in the immunocompromised population. Aiming to understand the population structure, microbiological characteristics and impact of allelic variation on β-lactamase structure and function, we collected 130 clinical isolates from across the United States. Identification of 90 different sequence types (STs), of which 63 are new allelic combinations, demonstrates the high diversity of this species. The majority of the isolates (45%) belong to genomic group 6. We also report excellent activity of the ceftazidime-avibactam and aztreonam combination, especially against strains recovered from blood and respiratory infections for which the susceptibility is higher than the susceptibility to trimethoprim-sulfamethoxazole, considered the "first-line" antibiotic to treat S. maltophilia Analysis of 73 bla L1 and 116 bla L2 genes identified 35 and 43 novel variants of L1 and L2 β-lactamases, respectively. Investigation of the derived amino acid sequences showed that substitutions are mostly conservative and scattered throughout the protein, preferentially affecting positions that do not compromise enzyme function but that may have an impact on substrate and inhibitor binding. Interestingly, we detected a probable association between a specific type of L1 and L2 and genomic group 6. Taken together, our results provide an overview of the molecular epidemiology of S. maltophilia clinical strains from the United States. In particular, the discovery of new L1 and L2 variants warrants further study to fully understand the relationship between them and the β-lactam resistance phenotype in this pathogen. IMPORTANCE Multiple antibiotic resistance mechanisms, including two β-lactamases, L1, a metallo-β-lactamase, and L2, a class A cephalosporinase, make S. maltophilia naturally multidrug resistant. Thus, infections caused by S. maltophilia pose a big therapeutic challenge. Our study aims to understand the microbiological and molecular characteristics of S. maltophilia isolates recovered from human sources. A highlight of the resistance profile of this collection is the excellent activity of the ceftazidime-avibactam and aztreonam combination. We hope this result prompts controlled and observational studies to add clinical data on the utility and safety of this therapy. We also identify 35 and 43 novel variants of L1 and L2, respectively, some of which harbor novel substitutions that could potentially affect substrate and/or inhibitor binding. We believe our results provide valuable knowledge to understand the epidemiology of this species and to advance mechanism-based inhibitor design to add to the limited arsenal of antibiotics active against this pathogen.

Acidithiobacillus ferrooxidans , a chemolithoautotrophic Gram-negative bacterium, has a remarkable ability to obtain energy from ferrous iron oxidation at pH 2. Several metalloproteins have been described as being involved in this respiratory chain coupling iron oxidation with oxygen reduction. However, their properties and physiological functions remain largely unknown, preventing a clear understanding of the global mechanism. In this work, we focus on two metalloproteins of this respiratory pathway, a diheme cytochrome c 4 (Cyt c 4 ) and a green copper protein (AcoP) of unknown function. We first demonstrate the formation of a complex between these two purified proteins, which allows homogeneous intermolecular electron-transfer in solution. We then mimic the physiological interaction between the two partners by replacing one at a time with electrodes displaying different chemical functionalities. From the electrochemical behavior of individual proteins, we show that, while electron transfer on AcoP requires weak electrostatic interaction, electron transfer on Cyt c 4 tolerates different charge and hydrophobicity conditions, suggesting a pivotal role of this protein in the metabolic chain. The electrochemical study of the proteins incubated together demonstrates an intermolecular electron transfer involving the protein complex, in which AcoP is reduced through the high potential heme of Cyt c 4 . Modelling of the electrochemical signals at different scan rates allows us to estimate the rate constant of this intermolecular electron transfer in the range of a few s -1 . Possible routes for electron transfer in the acidophilic bacterium are deduced.

Unlike in metazoans, plant microRNAs (miRNAs) undergo stepwise nuclear maturation before engaging cytosolic, sequence-complementary transcripts in association with the silencing effector protein ARGONAUTE1 (AGO1). Since their discovery, how and under which form plant miRNAs translocate to the cytosol has remained unclear, as has their sub-cellular AGO1 loading site(s). Here, we show that the N termini of all plant AGO1s contain a nuclear-localization (NLS) and nuclear-export signal (NES) that, in Arabidopsis thaliana (At), enables AtAGO1 nucleo-cytosolic shuttling in a Leptomycin-B-inhibited manner, diagnostic of CRM1(EXPO1)/NES-dependent nuclear export. Nuclear-only AtAGO1 contains the same 2'O-methylated miRNA cohorts as its nucleo-cytosolic counterpart, but it preferentially interacts with the miRNA loading chaperone HSP90. Furthermore, mature miRNA translocation and miRNA-mediated silencing both require AtAGO1 nucleo-cytosolic shuttling. These findings lead us to propose a substantially revised view of the plant miRNA pathway in which miRNAs are matured, methylated, loaded into AGO1 in the nucleus, and exported to the cytosol as AGO1:miRNA complexes in a CRM1(EXPO1)/NES-dependent manner., (Copyright © 2018 Elsevier Inc. All rights reserved.)

Progressive myoclonus epilepsies (PMEs) are inherited disorders characterized by myoclonus, generalized tonic-clonic seizures, and ataxia. One of the genes that is associated with PME is the ER-to-Golgi Q b -SNARE GOSR2, which forms a SNARE complex with syntaxin-5, Bet1 and Sec22b. Most PME patients are homo-zygous for a p.Gly144Trp mutation and develop similar clinical presentations. Recently, a patient who was compound heterozygous for p.Gly144Trp and a previously unseen p.Lys164del mutation was identified. Because this patient presented with a milder disease phenotype, we hypothesized that the p.Lys164del mutation may be less severe compared to p.Gly144Trp. To characterize the effect of the p.Gly144Trp and p.Lys164del mutations, both of which are present in the SNARE motif of GOSR2, we examined the corresponding mutations in the yeast ortholog Bos1. Yeasts expressing the orthologous mutants in Bos1 showed impaired growth, suggesting a partial loss of function, which was more severe for the Bos1 p.Gly176Trp mutation. Using anisotropy and gel filtration, we report that Bos1 p.Gly176Trp and p.Arg196del are capable of complex formation, but with partly reduced activity. Molecular dynamics (MD) simulations showed that the hydrophobic core, which triggers SNARE complex formation, is compromised due to the glycine-to-tryptophan substitution in both GOSR2 and Bos1. In contrast, the deletion of residue p.Lys164 (or p.Arg196del in Bos1) interferes with the formation of hydrogen bonds between GOSR2 and syntaxin-5. Despite these perturbations, all SNARE complexes stayed intact during longer simulations. Thus, our data suggest that the milder course of disease in compound heterozygous PME is due to less severe impairment of the SNARE function., Competing Interests: Competing interestsThe authors declare no competing or financial interests., (© 2017. Published by The Company of Biologists Ltd.)

Aerolysin is the founding member of a major class of β-pore-forming toxins (β-PFTs) found throughout all kingdoms of life. PFTs are cytotoxic proteins produced as soluble monomers, which oligomerize at the membrane of target host cells forming pores that may lead to osmotic lysis and cell death. Besides their role in microbial infection, they have become interesting for their potential as biotechnological sensors and delivery systems. Using an approach that integrates bioinformatics with molecular modeling and simulation, we looked for conserved features across this large toxin family. The cell surface-binding domains present high variability within the family to provide membrane receptor specificity. On the contrary, the novel concentric double β-barrel structure found in aerolysin is highly conserved in terms of sequence, structure and conformational dynamics, which likely contribute to preserve a common transition mechanism from the prepore to the mature pore within the family.Our results point to the key role of several amino acids in the conformational changes needed for oligomerization and further pore formation, such as Y221, W227, P248, Q263 and L277, which we propose are involved in the release of the stem loop and the two adjacent β-strands to form the transmembrane β-barrel.

This Briefing reviews the widely used, currently active, up-to-date databases derived from the worldwide Protein Data Bank (PDB) to facilitate browsing, finding and exploring its entries. These databases contain visualization and analysis tools tailored to specific kinds of molecules and interactions, often including also complex metrics precomputed by experts or external programs, and connections to sequence and functional annotation databases. Importantly, updates of most of these databases involves steps of curation and error checks based on specific expertise about the subject molecules or interactions, and removal of sequence redundancy, both leading to better data sets for mining studies compared with the full list of raw PDB entries. The article presents the databases in groups such as those aimed to facilitate browsing through PDB entries, their molecules and their general information, those built to link protein structure with sequence and dynamics, those specific for transmembrane proteins, nucleic acids, interactions of biomacromolecules with each other and with small molecules or metal ions, and those concerning specific structural features or specific protein families. A few webservers directly connected to active databases, and a few databases that have been discontinued but would be important to have back, are also briefly commented on. Along the Briefing, sample cases where these databases have been used to aid structural studies or advance our knowledge about biological macromolecules are referenced. A few specific examples are also given where using these databases is easier and more informative than using raw PDB data., (© The Author 2016. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

Adhesion to host tissue is one of the key steps of the bacterial pathogenic process. Xanthomonas citri ssp. citri possesses a non-fimbrial adhesin protein, XacFhaB, required for bacterial attachment, which we have previously demonstrated to be an important virulence factor for the development of citrus canker. XacFhaB is a 4753-residue-long protein with a predicted β-helical fold structure, involved in bacterial aggregation, biofilm formation and adhesion to the host. In this work, to further characterize this protein and considering its large size, XacFhaB was dissected into three regions based on bioinformatic and structural analyses for functional studies. First, the capacity of these protein regions to aggregate bacterial cells was analysed. Two of these regions were able to form bacterial aggregates, with the most amino-terminal region being dispensable for this activity. Moreover, XacFhaB shows features resembling pathogen-associated molecular patterns (PAMPs), which are recognized by plants. As PAMPs activate plant basal immune responses, the role of the three XacFhaB regions as elicitors of these responses was investigated. All adhesin regions were able to induce basal immune responses in host and non-host plants, with a stronger activation by the carboxyl-terminal region. Furthermore, pre-infiltration of citrus leaves with XacFhaB regions impaired X. citri ssp. citri growth, confirming the induction of defence responses and restraint of citrus canker. This work reveals that adhesins from plant pathogens trigger plant defence responses, opening up new pathways for the development of protective strategies for disease control., (© 2016 BSPP and John Wiley & Sons Ltd.)

Most studies of protein structure and function are performed in dilute conditions, but proteins typically experience high solute concentrations in their physiological scenarios and biotechnological applications. High solute concentrations have well-known effects on coarse protein traits like stability, diffusion, and shape, but likely also perturb other traits through finer effects pertinent at the residue and atomic levels. Here, NMR and molecular dynamics investigations on ubiquitin disclose variable interactions with concentrated solutes that lead to localized perturbations of the protein's surface, hydration, electrostatics, and dynamics, all dependent on solute size and chemical properties. Most strikingly, small polar uncharged molecules are sticky on the protein surface, whereas charged small molecules are not, but the latter still perturb the internal protein electrostatics as they diffuse nearby. Meanwhile, interactions with macromolecular crowders are favored mainly through hydrophobic, but not through polar, surface patches. All the tested small solutes strongly slow down water exchange at the protein surface, whereas macromolecular crowders do not exert such strong perturbation. Finally, molecular dynamics simulations predict that unspecific interactions slow down microsecond- to millisecond-timescale protein dynamics despite having only mild effects on pico- to nanosecond fluctuations as corroborated by NMR. We discuss our results in the light of recent advances in understanding proteins inside living cells, focusing on the physical chemistry of quinary structure and cellular organization, and we reinforce the idea that proteins should be studied in native-like media to achieve a faithful description of their function., (Copyright © 2016 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Parainfluenza virus 5 (PIV5) is an enveloped, single-stranded, negative-sense RNA virus of the Paramyxoviridae family. PIV5 fusion and entry are mediated by the coordinated action of the receptor-binding protein, hemagglutinin-neuraminidase (HN), and the fusion protein (F). Upon triggering by HN, F undergoes an irreversible ATP- and pH-independent conformational change, going down an energy gradient from a metastable prefusion state to a highly stable postfusion state. Previous studies have highlighted key conformational changes in the F-protein refolding pathway, but a detailed understanding of prefusion F-protein metastability remains elusive. Here, using two previously described F-protein mutations (S443D or P22L), we examine the capacity to modulate PIV5 F stability and the mechanisms by which these point mutants act. The S443D mutation destabilizes prefusion F proteins by disrupting a hydrogen bond network at the base of the F-protein globular head. The introduction of a P22L mutation robustly rescues destabilized F proteins through a local hydrophobic interaction between the N-terminal helix and a hydrophobic pocket. Prefusion stabilization conferred by a P22L-homologous mutation is demonstrated in the F protein of Newcastle disease virus, a paramyxovirus of a different genus, suggesting a conserved stabilizing structural element within the paramyxovirus family. Taken together, the available data suggest that movement of the N-terminal helix is a necessary early step for paramyxovirus F-protein refolding and presents a novel target for structure-based drug design.

The understanding of protein evolution depends on the ability to relate the impact of mutations on molecular traits to organismal fitness. Biological activity and robustness have been regarded as important features in shaping protein evolutionary landscapes. Conformational dynamics, which is essential for protein function, has received little attention in the context of evolutionary analyses. Here we employ NMR spectroscopy, the chief experimental tool to describe protein dynamics at atomic level in solution at room temperature, to study the intrinsic dynamic features of a metallo- Β: -lactamase enzyme and three variants identified during a directed evolution experiment that led to an expanded substrate profile. We show that conformational dynamics in the catalytically relevant microsecond to millisecond timescale is optimized along the favored evolutionary trajectory. In addition, we observe that the effects of mutations on dynamics are epistatic. Mutation Gly262Ser introduces slow dynamics on several residues that surround the active site when introduced in the wild-type enzyme. Mutation Asn70Ser removes the slow dynamics observed for few residues of the wild-type enzyme, but increases the number of residues that undergo slow dynamics when introduced in the Gly262Ser mutant. These effects on dynamics correlate with the epistatic interaction between these two mutations on the bacterial phenotype. These findings indicate that conformational dynamics is an evolvable trait, and that proteins endowed with more dynamic active sites also display a larger potential for promoting evolution., (© The Author 2016. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Background: Protein variability can now be studied by measuring high-resolution tolerance-to-substitution maps and fitness landscapes in saturated mutational libraries. But these rich and expensive datasets are typically interpreted coarsely, restricting detailed analyses to positions of extremely high or low variability or dubbed important beforehand based on existing knowledge about active sites, interaction surfaces, (de)stabilizing mutations, etc., Results: Our new webserver PsychoProt (freely available without registration at http://psychoprot.epfl.ch or at http://lucianoabriata.altervista.org/psychoprot/index.html ) helps to detect, quantify, and sequence/structure map the biophysical and biochemical traits that shape amino acid preferences throughout a protein as determined by deep-sequencing of saturated mutational libraries or from large alignments of naturally occurring variants., Discussion: We exemplify how PsychoProt helps to (i) unveil protein structure-function relationships from experiments and from alignments that are consistent with structures according to coevolution analysis, (ii) recall global information about structural and functional features and identify hitherto unknown constraints to variation in alignments, and (iii) point at different sources of variation among related experimental datasets or between experimental and alignment-based data. Remarkably, metabolic costs of the amino acids pose strong constraints to variability at protein surfaces in nature but not in the laboratory. This and other differences call for caution when extrapolating results from in vitro experiments to natural scenarios in, for example, studies of protein evolution., Conclusion: We show through examples how PsychoProt can be a useful tool for the broad communities of structural biology and molecular evolution, particularly for studies about protein modeling, evolution and design.

Unlabelled: Mycobacterium tuberculosis possesses a thick and highly hydrophobic cell wall principally composed of a mycolyl-arabinogalactan-peptidoglycan complex, which is critical for survival and virulence. DprE1 is a well-characterized component of decaprenyl-phospho-ribose epimerase, which produces decaprenyl-phospho-arabinose (DPA) for the biosynthesis of mycobacterial arabinans. Upstream of dprE1 lies rv3789, which encodes a short transmembrane protein of the GtrA family, whose members are often involved in the synthesis of cell surface polysaccharides. We demonstrate that rv3789 and dprE1 are cotranscribed from a common transcription start site situated 64 bp upstream of rv3789. Topology mapping revealed four transmembrane domains in Rv3789 and a cytoplasmic C terminus consistent with structural models built using analysis of sequence coevolution. To investigate its role, we generated an unmarked rv3789 deletion mutant in M. tuberculosis. The mutant was characterized by impaired growth and abnormal cell morphology, since the cells were shorter and more swollen than wild-type cells. This phenotype likely stems from the decreased incorporation of arabinan into arabinogalactan and was accompanied by an accumulation of DPA. A role for Rv3789 in arabinan biosynthesis was further supported by its interaction with the priming arabinosyltransferase AftA, as demonstrated by a two-hybrid approach. Taken together, the data suggest that Rv3789 does not act as a DPA flippase but, rather, recruits AftA for arabinogalactan biosynthesis., Importance: Upstream of the essential dprE1 gene, encoding a key enzyme of the decaprenyl phospho-arabinose (DPA) pathway, lies rv3789, coding for a short transmembrane protein of unknown function. In this study, we demonstrated that rv3789 and dprE1 are cotranscribed from a common transcription start site located 64 bp upstream of rv3789 in M. tuberculosis. Furthermore, the deletion of rv3789 led to a reduction in arabinan content and to an accumulation of DPA, confirming that Rv3789 plays a role in arabinan biosynthesis. Topology mapping, structural modeling, and protein interaction studies suggest that Rv3789 acts as an anchor protein recruiting AftA, the first arabinosyl transferase. This investigation provides deeper insight into the mechanism of arabinan biosynthesis in mycobacteria., (Copyright © 2015, American Society for Microbiology. All Rights Reserved.)

The thermosensor histidine kinase DesK from Bacillus subtilis senses changes in membrane fluidity initiating an adaptive response. Structural changes in DesK have been implicated in transmembrane signaling, but direct evidence is still lacking. On the basis of structure-guided mutagenesis, we now propose a mechanism of DesK-mediated signal sensing and transduction. The data indicate that stabilization/destabilization of a 2-helix coiled coil, which connects the transmembrane sensory domain of DesK to its cytosolic catalytic region, is crucial to control its signaling state. Computational modeling and simulations reveal couplings between protein, water and membrane mechanics. We propose that membrane thickening is the main driving force for signal sensing and that it acts by inducing helix stretching and rotation prompting an asymmetric kinase-competent state. Overall, the known structural changes of the sensor kinase, as well as further dynamic rearrangements that we now predict, consistently link structure determinants to activity modulation., (© 2015 John Wiley & Sons Ltd.)

Maturation of cytochrome oxidases is a complex process requiring assembly of several subunits and adequate uptake of the metal cofactors. Two orthologous Sco proteins (Sco1 and Sco2) are essential for the correct assembly of the dicopper CuA site in the human oxidase, but their function is not fully understood. Here, we report an in vitro biochemical study that shows that Sco1 is a metallochaperone that selectively transfers Cu(I) ions based on loop recognition, whereas Sco2 is a copper-dependent thiol reductase of the cysteine ligands in the oxidase. Copper binding to Sco2 is essential to elicit its redox function and as a guardian of the reduced state of its own cysteine residues in the oxidizing environment of the mitochondrial intermembrane space (IMS). These results provide a detailed molecular mechanism for CuA assembly, suggesting that copper and redox homeostasis are intimately linked in the mitochondrion.

Protein-protein recognition and binding are governed by diffusion, noncovalent forces and conformational flexibility, entangled in a way that only molecular dynamics simulations can dissect at high resolution. Here we exploited ubiquitin's noncovalent dimerization equilibrium to assess the potential of atomistic simulations to reproduce reversible protein-protein binding, by running submicrosecond simulations of systems with multiple copies of the protein at millimolar concentrations. The simulations essentially fail because they lead to aggregates, yet they reproduce some specificity in the binding interfaces as observed in known covalent and noncovalent ubiquitin dimers. Following similar observations in literature we hint at electrostatics and water descriptions as the main liable force field elements, and propose that their optimization should consider observables relevant to multi-protein systems and unfolded proteins. Within limitations, analysis of binding events suggests salient features of protein-protein recognition and binding, to be retested with improved force fields. Among them, that specific configurations of relative direction and orientation seem to trigger fast binding of two molecules, even over 50 Å distances; that conformational selection can take place within surface-to-surface distances of 10 to 40 Å i.e. well before actual intermolecular contact; and that establishment of contacts between molecules further locks their conformations and relative orientations.

The need for interfacing structural biology and biophysics to molecular evolution is being increasingly recognized. One part of the big problem is to understand how physics and chemistry shape the sequence space available to functional proteins, while satisfying the needs of biology. Here we present a quantitative, structure-based analysis of a high-resolution map describing the tolerance to all substitutions in all positions of a functional enzyme, namely a TEM lactamase previously studied through deep sequencing of mutants growing in competition experiments with selection against ampicillin. Substitutions are rarely observed within 7 Å of the active site, a stringency that is relaxed slowly and extends up to 15-20 Å, with buried residues being especially sensitive. Substitution patterns in over one third of the residues can be quantitatively modeled by monotonic dependencies on amino acid descriptors and predictions of changes in folding stability. Amino acid volume and steric hindrance shape constraints on the protein core; hydrophobicity and solubility shape constraints on hydrophobic clusters underneath the surface, and on salt bridges and polar networks at the protein surface together with charge and hydrogen bonding capacity. Amino acid solubility, flexibility and conformational descriptors also provide additional constraints at many locations. These findings provide fundamental insights into the chemistry underlying protein evolution and design, by quantitating links between sequence and different protein traits, illuminating subtle and unexpected sequence-trait relationships and pinpointing what traits are sacrificed upon gain-of-function mutation.

We present herein a thorough description of the effects of high glucose concentrations on the diffusion, hydration and internal dynamics of ubiquitin, as predicted from extensive molecular dynamics simulations on several systems described at fully atomistic level. We observe that the protein acts as a seed that speeds up the natural propensity of glucose to cluster at high concentration; the sugar molecules thus aggregate around the protein trapping it inside a dynamic cage. This process extensively dehydrates the protein surface, restricts the motions of the remaining water molecules, and drags the large-scale, collective motions of protein atoms slowing down the rate of exploration of the conformational space despite only a slight dampening of fast, local dynamics. We discuss how these effects could be relevant to the function of sugars as preservation agents in biological materials, and how crowding by small sticky molecules could modulate proteins across different reaction coordinates inside the cellular cytosol.

The influence of postharvest fruit ripening in the composition of metabolites, transcripts and enzymes in tomato (Solanum lycopersicum L.) is poorly understood. The goal of this work was to study the changes in the metabolite composition of the tomato fruit ripened off-the-vine using the cultivar Micro-Tom as model system. Proton nuclear magnetic resonance (1H NMR) was used for analysis of the metabolic profile of tomato fruits ripened on- and off-the-vine. Significant differences under both ripening conditions were observed principally in the contents of fructose, glucose, aspartate and glutamate. Transcript levels and enzyme activities of -amino butyrate transaminase (EC 2.6.1.19) and glutamate decarboxylase (EC 4.1.1.15) showed differences in fruits ripened under these two conditions. These data indicate that the contents of metabolites involved in primary metabolism, and conferring the palatable properties of fruits, are altered when fruits are ripened off-the-vine.

Protein X-ray structures with non-corrin cobalt(II)-containing sites, either natural or substituting another native ion, were downloaded from the Protein Data Bank and explored to (i) describe which amino acids are involved in their first ligand shells and (ii) analyze cobalt(II)-donor bond lengths in comparison with previously reported target distances, CSD data and EXAFS data. The set of amino acids involved in Co(II) binding is similar to that observed for catalytic Zn(II) sites, i.e. with a large fraction of carboxylate O atoms from aspartate and glutamate and aromatic N atoms from histidine. The computed Co(II)-donor bond lengths were found to depend strongly on structure resolution, an artifact previously detected for other metal-donor distances. Small corrections are suggested for the target bond lengths to the aromatic N atoms of histidines and the O atoms of water and hydroxide. The available target distance for cysteine (Scys) is confirmed; those for backbone O and other donors remain uncertain and should be handled with caution in refinement and modeling protocols. Finally, a relationship between both Co(II)-O bond lengths in bidentate carboxylates is quantified.

Electron transfer is the simplest chemical reaction and constitutes the basis of a large variety of biological processes, such as photosynthesis and cellular respiration. Nature has evolved specific proteins and cofactors for these functions. The mechanisms optimizing biological electron transfer have been matter of intense debate, such as the role of the protein milieu between donor and acceptor sites. Here we propose a mechanism regulating long-range electron transfer in proteins. Specifically, we report a spectroscopic, electrochemical, and theoretical study on WT and single-mutant Cu(A) redox centers from Thermus thermophilus, which shows that thermal fluctuations may populate two alternative ground-state electronic wave functions optimized for electron entry and exit, respectively, through two different and nearly perpendicular pathways. These findings suggest a unique role for alternative or "invisible" electronic ground states in directional electron transfer. Moreover, it is shown that this energy gap and, therefore, the equilibrium between ground states can be fine-tuned by minor perturbations, suggesting alternative ways through which protein-protein interactions and membrane potential may optimize and regulate electron-proton energy transduction.

A dataset of TEM lactamase variants with different substrate and inhibition profiles was compiled and analyzed. Trends show that loops are the main evolvable regions in these enzymes, gradually accumulating mutations to generate increasingly complex functions. Notably, many mutations present in evolved enzymes are also found in simpler variants, probably originating functional promiscuity. Following a function-stability tradeoff, the increase in functional complexity driven by accumulation of mutations fosters the incorporation of other stability-restoring substitutions, although our analysis suggests they might not be as "global" as generally accepted and seem instead specific to different networks of protein sites. Finally, we show how this dataset can be used to model functional changes in TEMs based on the physicochemical properties of the amino acids., (Copyright © 2012 Federation of European Biochemical Societies. Published by Elsevier B.V. All rights reserved.)

Protein-mediated electron transfer is an essential event in many biochemical processes. Efficient electron transfer requires the reorganization energy of the redox event to be minimized, which is ensured by the presence of rigid donor and acceptor sites. Electron transfer copper sites are present in the ubiquitous cupredoxin fold, able to bind one or two copper ions. The low reorganization energy in these metal centers has been accounted for by assuming that the protein scaffold creates an entatic/rack-induced state, which gives rise to a rigid environment by means of a preformed metal chelating site. However, this notion is incompatible with the need for an exposed metal-binding site and protein-protein interactions enabling metallochaperone-mediated assembly of the copper site. Here we report an NMR study that reveals a high degree of structural heterogeneity in the metal-binding region of the nonmetallated Cu(A)-binding cupredoxin domain, arising from microsecond to second dynamics that are quenched upon metal binding. We also report similar dynamic features in apo-azurin, a paradigmatic blue copper protein, suggesting a general behavior. These findings reveal that the entatic/rack-induced state, governing the features of the metal center in the copper-loaded protein, does not require a preformed metal-binding site. Instead, metal binding is a major contributor to the rigidity of electron transfer copper centers. These results reconcile the seemingly contradictory requirements of a rigid, occluded center for electron transfer, and an accessible, dynamic site required for in vivo copper uptake.

Cytochrome c, a mitochondrial electron transfer protein containing a hexacoordinated heme, is involved in other physiologically relevant events, such as the triggering of apoptosis, and the activation of a peroxidatic activity. The latter occurs secondary to interactions with cardiolipin and/or post-translational modifications, including tyrosine nitration by peroxynitrite and other nitric oxide-derived oxidants. The gain of peroxidatic activity in nitrated cytochrome c has been related to a heme site transition in the physiological pH region, which normally occurs at alkaline pH in the native protein. Herein, we report a spectroscopic characterization of two nitrated variants of horse heart cytochrome c by using optical spectroscopy studies and NMR. Highly pure nitrated cytochrome c species modified at solvent-exposed Tyr-74 or Tyr-97 were generated after treatment with a flux of peroxynitrite, separated, purified by preparative high pressure liquid chromatography, and characterized by mass spectrometry-based peptide mapping. It is shown that nitration of Tyr-74 elicits an early alkaline transition with a pKa = 7.2, resulting in the displacement of the sixth and axial iron ligand Met-80 and replacement by a weaker Lys ligand to yield an alternative low spin conformation. Based on the study of site-specific Tyr to Phe mutants in the four conserved Tyr residues, we also show that this transition is not due to deprotonation of nitro-Tyr-74, but instead we propose a destabilizing steric effect of the nitro group in the mobile Omega-loop of cytochrome c, which is transmitted to the iron center via the nearby Tyr-67. The key role of Tyr-67 in promoting the transition through interactions with Met-80 was further substantiated in the Y67F mutant. These results therefore provide new insights into how a remote post-translational modification in cytochrome c such as tyrosine nitration triggers profound structural changes in the heme ligation and microenvironment and impacts in protein function.