Search

International audience; In all organisms, translational initiation takes place on the small ribosomal subunit and two classes of methionine tRNA are present. The initiator is used exclusively for initiation of protein synthesis while the elongator is used for inserting methionine internally in the nascent polypeptide chain. The crystal structure of Escherichia coli initiator tRNA(f)(Met) has been solved at 3.1 A resolution. The anticodon region is well-defined and reveals a unique structure, which has not been described in any other tRNA. It encompasses a Cm32*A38 base pair with a peculiar geometry extending the anticodon helix, a base triple between A37 and the G29-C41 pair in the major groove of the anticodon stem and a modified stacking organization of the anticodon loop. This conformation is associated with the three GC basepairs in the anticodon stem, characteristic of initiator tRNAs and suggests a mechanism by which the translation initiation machinery could discriminate the initiator tRNA from all other tRNAs.

Recently, we observed the first example of a left‐handed G‐quadruplex structure formed by natural DNA, named Z‐G4 . We analysed the Z‐G4 structure and inspected its primary 28‐nt sequence in order to identify motifs that convey the unique left‐handed twist. Using circular dichroism spectroscopy, NMR spectroscopy, and X‐ray crystallography, we revealed a minimal sequence motif of 12 nt (GTGGTGGTGGTG) for formation of the left‐handed DNA G‐quadruplex, which is found to be highly abundant in the human genome. A systematic analysis of thymine loop mutations revealed a moderate sequence tolerance, which would further broaden the space of sequences prone to left‐handed G‐quadruplex formation. NRF (Natl Research Foundation, S’pore) MOE (Min. of Education, S’pore) Accepted version

International audience; Computational protein design (CPD) aims to predict amino acid sequences that fold to specific structures and perform desired functions. CPD depends on a rotamer library, an energy function, and an algorithm to search the sequence/conformation space. Variable neighborhood search (VNS) with cost function networks is a powerful framework that can provide tight upper bounds on the global minimum energy. We propose a new CPD heuristic based on VNS in which a subset of the solution space (a “neighborhood”) is explored, whose size is gradually increased with a dedicated probabilistic heuristic. The algorithm was tested on 99 protein designs with fixed backbones involving nine proteins from the SH2, SH3, and PDZ families. The number of mutating positions was 20, 30, or all of the amino acids, while the rest of the protein explored side-chain rotamers. VNS was more successful than Monte Carlo (MC), replica-exchange MC, and a heuristic steepest-descent energy minimization, providing solutions with equal or lower best energies in most cases. For complete protein redesign, it gave solutions that were 2.5 to 11.2 kcal/mol lower in energy than those obtained with the other approaches. VNS is implemented in the toulbar2 software. It could be very helpful for large and/or complex design problems.

Une correction de l'article est publiée sur le site de l'éditeur : https://pubs.acs.org/doi/abs/10.1021/acs.jpca.8b12337; International audience; Phosphate groups are essential components of nucleic acids and proteins, whose interactions with solvent, metal ions, and ionic side chains help control folding and binding. Methyl phosphate (MP) represents a simple analog of phosphate moieties that are post-translation modifications in proteins and present at the termini of nucleic acids, among other environments. In the present study, we optimized parameters for use in polarizable molecular dynamics simulations of MP in its mono- and dianionic forms, MP- ≡ CH3HPO4- and MP2- ≡ CH3PO42-, along with P i2- ≡ HPO42-, in the context of the classical Drude oscillator model. Parameter optimization was done in a manner consistent with the remainder of the Drude molecular mechanics force field, choosing atomic charges and polarizabilities to reproduce molecular properties from quantum mechanics as well as experimental hydration free energies. Optimized parameters were similar to existing dimethyl phosphate parameters, with a few significant differences. The developed parameters were then used to compute magnesium binding affinities in aqueous solution, using alchemical molecular dynamics free energy simulations. Good agreement with experiment was obtained, and outer sphere binding was shown to be predominant for MP- and MP2-.

International audience; One proposed scenario for the emergence of biochemical oscillations is that they may have provided the basic mechanism behind cellular self-replication by growth and division. However, alternative scenarios not requiring any chemical oscillation have also been proposed. Each of the various protocell models proposed to support one or another scenario comes with its own set of specific assumptions, which makes it difficult to ascertain whether chemical oscillations are required or not for cellular self-replication. This article compares these two cases within a single whole-cell model framework. This model relies upon a membrane embedding a chemical reaction network (CRN) synthesizing all the cellular constituents, including the membrane, by feeding from an external nutrient. Assuming the osmolarity is kept constant, the system dynamics are governed by a set of nonlinear differential equations coupling the chemical concentrations and the surface-area-to-volume ratio. The resulting asymptotic trajectories are used to determine the cellular shape by minimizing the membrane bending energy (within an approximate predefined family of shapes). While the stationary case can be handled quite generally, the oscillatory one is investigated using a simple oscillating CRN example, which is used to identify features that are expected to hold for any network. It is found that cellular self-replication can be reached with or without chemical oscillations, and that a requirement common to both stationary and oscillatory cases is that a minimum spontaneous curvature of the membrane is required for the cell to divide once its area and volume are both doubled. The oscillatory case can result in a greater variety of cellular shape trajectories but raises additional constraints for cellular division and self-replication: (i) the ratio of doubling time to oscillation period should be an integer, and (ii) if the oscillation amplitude is sufficiently high, then the spontaneous curvature must be below a maximum value to avoid early division before the end of the cycle. Because of these additional stringent constraints, it is likely that early protocells did not rely upon chemical oscillations. Biochemical oscillations typical of modern evolved cells may have emerged later through evolution for other reasons (e.g., metabolic advantage) and must have required additional feedback mechanisms for such a self-replicating system to be robust against even slight environmental variations (e.g., temperature fluctuations).

International audience; Condensed-phase simulations commonly use periodic boundary conditions (PBCs) to represent the thermodynamic limit. For the vapor to liquid transfer of an ion, the gas/liquid boundary and its associated potential change are then missing. Furthermore, the electric potential and field at a given point are given by conditionally convergent infinite series, for which different summation schemes give different results. Nevertheless, standard simulation protocols can be used to compute experimental quantities unambiguously. In particular, using an auxiliary test particle and a multistep solvation path, we show that particle-based, Ewald, and common molecule-based summation schemes for the potential and field are all essentially equivalent. However, all methods require prior knowledge of the gas/liquid boundary potential to compute ionic solvation free energies using PBC protocols for both force-field and quantum-mechanical models.

Calcium is a second messenger crucial to a myriad of cellular processes ranging from regulation of metabolism and cell survival to vesicle release and motility. Current strategies to directly manipulate endogenous calcium signals lack cellular and subcellular specificity. We developed SpiCee, a versatile and genetically-encoded chelator combining low and high affinity sites for calcium. This scavenger enables altering endogenous calcium signaling and functions in single cells in vitro and in vivo with biochemically-controlled subcellular resolution. SpiCee paves the way to investigate local calcium signaling in vivo and directly manipulate this second messenger for therapeutic use.

International audience; More than 100 naturally occurring modified nucleotides have been found in RNA molecules, in particular in tRNAs. We have determined molecular mechanics force field parameters compatible with the CHARMM36 all-atom additive force field for all these modifications using the CHARMM force field parametrization strategy. Emphasis was placed on fine tuning of the partial atomic charges and torsion angle parameters. Quantum mechanics calculations on model compounds provided the initial set of target data, and extensive molecular dynamics simulations of nucleotides and oligonucleotides in aqueous solutions were used for further refinement against experimental data. The presented parameters will allow for computational studies of a wide range of RNAs containing modified nucleotides, including the ribosome and transfer RNAs. (C) 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

International audience; For the high throughput design of protein:peptide binding, one must explore a vast space of amino acid sequences in search of low binding free energies. This complex problem is usually addressed with either simple heuristic scoring or expensive sequence enumeration schemes. Far more efficient than enumeration is a recent Monte Carlo approach that adaptively flattens the energy landscape in sequence space of the unbound peptide and provides formally exact binding free energy differences. The method allows the binding free energy to be used directly as the design criterion. We propose several improvements that allow still more efficient sampling and can address larger design problems. They include the use of Replica Exchange Monte Carlo and landscape flattening for both the unbound and bound peptides. We used the method to design peptides that bind to the PDZ domain of the Tiam1 signaling protein and could serve as inhibitors of its activity. Four peptide positions were allowed to mutate freely. Almost 75 000 peptide variants were processed in two simulations of 109 steps each that used 1 CPU hour on a desktop machine. 96% of the theoretical sequence space was sampled. The relative binding free energies agreed qualitatively with values from experiment. The sampled sequences agreed qualitatively with an experimental library of Tiam1-binding peptides. The main assumption limiting accuracy is the fixed backbone approximation, which could be alleviated in future work by using increased computational resources and multi-backbone designs.

International audience; Cell migration is essential to many biological processes such as embryonic development, immune surveillance and wound healing. Random cell migration refers to the intrinsic ability of cells to migrate, often called cell motility. This basal condition contrasts with directed cell migration, where cells migrate toward a chemical or physical cue. Unlike Brownian particles, however, randomly migrating cells exhibit a directional persistence, i.e., they are more likely to sustain the movement in the direction they previously took than to change, even if this direction is randomly chosen in an isotropic environment. Here we describe how to set up time-lapse recording of mammalian cells freely moving on a two-dimensional surface coated with extracellular matrix proteins, how to acquire single cell trajectories from movies and how to extract key parameters that characterize cell motility, such as cell speed, directionality, mean square displacement, and directional persistence.

The potentially structured core domain of the intrinsically disordered protein Knr4 fromSaccharomyces cerevisiae, comprising residues 80–340, was expressed inEscherichia coliand crystallized using the hanging-drop vapour-diffusion method. Selenomethionine-containing (SeMet) protein was also purified and crystallized. Crystals of both proteins belonged to space groupP6522, with unit-cell parametersa=b= 112.44,c= 265.21 Å for the native protein anda = b = 112.49,c= 262.21 Å for the SeMet protein, and diffracted to 3.50 and 3.60 Å resolution, respectively. There are two molecules in the asymmetric unit related by a twofold axis. The anomalous signal of selenium was recorded and yielded an electron-density map of sufficient quality to allow the identification of secondary-structure elements.

International audience; GTPases typically switch between an inactive, OFF conformation and an active, ON conformation when a GDP ligand is replaced by GTP. Their ON/OFF populations and activity thus depend on the stabilities of four protein complexes, two apo-protein forms, and GTP/GDP in solution. A complete characterization is usually not possible experimentally and poses major challenges for simulations. We review the most important methodological challenges and we review thermodynamic data for two GTPases involved in translation of the genetic code: archaeal Initiation Factors 2 and 5B (aIF2, aIF5B). One main challenge is the multiplicity of states and conformations, including those of GTP/GDP in solution. Another is force field accuracy, especially for interactions of GTP/GDP with co-bound divalent Mg2 + ions. The calculation of electrostatic free energies also poses specific challenges, and requires careful protocols. For aIF2, experiments and earlier simulations showed that it is a “classic” GTPase, with distinct ON/OFF conformations that prefer to bind GTP and GDP, respectively. For aIF5B, we recently proposed a non-classic mechanism, where the ON/OFF states differ only in the protonation state of Glu81 in the nucleotide binding pocket. This model is characterized here using free energy simulations. The methodological analysis should help future studies, while the aIF2, aIF5B examples illustrate the diversity of ATPase/GTPase mechanisms. This article is part of a Special Issue entitled Recent developments of molecular dynamics.

International audience; Eukaryotic and archaeal translation initiation processes involve a heterotrimeric GTPase e/aIF2 crucial for accuracy of start codon selection. In eu-karyotes, the GTPase activity of eIF2 is assisted by a GTPase-activating protein (GAP), eIF5. In ar-chaea, orthologs of eIF5 are not found and aIF2 GT-Pase activity is thought to be non-assisted. However , no in vitro GTPase activity of the archaeal factor has been reported to date. Here, we show that aIF2 significantly hydrolyses GTP in vitro. Within aIF2␥, H97, corresponding to the catalytic histidine found in other translational GTPases, and D19, from the GKT loop, both participate in this activity. Several high-resolution crystal structures were determined to get insight into GTP hydrolysis by aIF2␥. In particular, a crystal structure of the H97A mutant was obtained in the presence of non-hydrolyzed GTP. This structure reveals the presence of a second magnesium ion bound to GTP and D19. Quantum chemical/molecular mechanical simulations support the idea that the second magnesium ion may assist GTP hydrolysis by helping to neutralize the developing negative charge in the transition state. These results are discussed in light of the absence of an identified GAP in archaea to assist GTP hydrolysis on aIF2.

International audience; Computational protein design aims to create proteins with novel properties. A key element is the energy or scoring function used to select the sequences and conformations. We study the performance of an "MMGBSA" energy function, which combines molecular mechanics terms, a generalized Born and surface area (GBSA) solvent model, with approximations that make the model pairwise additive. Our approach is implemented in the Proteus software. The use of a physics-based energy function ensures a certain model transferability and explanatory power. As a first test, we redesign the sequence of nine proteins, one position at a time, with the rest of the protein having its native sequence and crystallographic conformation. As a second test, all positions are designed together. The contributions of individual energy terms are evaluated, and various parametrizations are compared. We find that the GB term significantly improves the results compared to simple Coulomb electrostatics but is affected by pairwise decomposition errors when all positions are designed together. The SA term, with distinct energy coefficients for nonpolar and polar atoms, makes a decisive contribution to obtain realistic protein sequences and can partially compensate for the absence of a GB term. With the best GBSA protocol, we obtain nativelike protein cores and Superfamily recognition of almost all of our sequences.

Selenium is an essential trace element due to its incorporation into selenoproteins with important biological functions. However, at high doses it is toxic. Selenium toxicity is generally attributed to the induction of oxidative stress. However, it has become apparent that the mode of action of seleno-compounds varies, depending on its chemical form and speciation. Recent studies in various eukaryotic systems, in particular the model organism Saccharomyces cerevisiae, provide new insights on the cytotoxic mechanisms of selenomethionine and selenocysteine. This review first summarizes current knowledge on reactive oxygen species (ROS)-induced genotoxicity of inorganic selenium species. Then, we discuss recent advances on our understanding of the molecular mechanisms of selenocysteine and selenomethionine cytotoxicity. We present evidences indicating that both oxidative stress and ROS-independent mechanisms contribute to selenoamino acids cytotoxicity. These latter mechanisms include disruption of protein homeostasis by selenocysteine misincorporation in proteins and/or reaction of selenols with protein thiols.

© 2014 Elsevier Ltd. All rights reserved. Upf1, Upf2, and Upf3 are the principal regulators of nonsense-mediated mRNA decay (NMD), a cytoplasmic surveillance pathway that accelerates the degradation of mRNAs undergoing premature translation termination. These three proteins interact with each other, the ribosome, the translation termination machinery, and multiple mRNA decay factors, but the precise mechanism allowing the selective detection and degradation of nonsense-containing transcripts remains elusive. Here, we have determined the crystal structure of the N-terminal mIF4G domain from Saccharomyces cerevisiae Upf2 and identified a highly conserved region in this domain that is essential for NMD and independent of Upf2's binding sites for Upf1 and Upf3. Mutations within this conserved region not only inactivate NMD but also disrupt Upf2 binding to specific proteins, including Dbp6, a DEAD-box helicase. Although current models indicate that Upf2 functions principally as an activator of Upf1 and a bridge between Upf1 and Upf3, our data suggest that it may also serve as a platform for the association of additional factors that play roles in premature translation termination and NMD., This work was supported by the Centre National pour la Recherche Scientifique ATIP-AVENIR program (to M.G.), the Agence Nationale pour la Recherche (grant ANR-06-BLAN-0075-02 to H.v.T., S.K., and M.G.), the Association Française Contre les Myopathies (to H.v.T., S.K., and M.G.), the Fondation pour la Recherche Médicale (to Z.F.), the European Union Sixth Framework program “3D-Repertoire” (LSHG-CT-2005-512028, to H.v.T. and M.G.), the National Institutes of Health (R37 GM27757 to A.J.), the Ministerio de Economía y Competitividad (BFU2012-35367 to I.U.), the Centro para el Desarrollo Tecnológico Industrial (IDC-2010-1173 to I.U.), and the Human Frontiers Science Program (grant RGP0018, to H.v.T., S.K., A.J., and M.G.)

Analysis of infrared spectra of palladium nanoparticles (NPs) immersed in the tri-tert-butyl-R-phosphonium-based ionic liquids (ILs) demonstrates that both cations and anions of the ILs interact with the NPs. According to quantum-chemical simulations of these interactions, the binding energy of anions to the Pd-6 cluster, taken as a minimal-size model of the NPs, increases from similar to 6 to similar to 27 kcal mol(-1) in the order [PF6](-) approximate to [BF4](-) < [Tf2N](-) < [OTf](-) < [Br](-) << [TFA](-). In contrast, the binding energy for all types of the [(Bu3PR)-P-t](+) cations slightly varies at about similar to 22 kcal mol(-1) only moderately depending on the choice of the R moiety (n-pentyl, 2-hydroxyethyl, 2-methoxyethyl, 2-ethoxy-2-oxoethyl). As a result, the energies of interaction between a Pd-6 cluster and various ion pairs, formed by the abovementioned counter-ions, follow the order found for the anions and vary from similar to 24 to similar to 47 kcal mol(-1). These values are smaller than the energy of addition of a Pd atom to a Pd-n cluster (similar to 58 kcal mol(-1)), which suggests kinetic stabilization of the NPs in phosphonium-based ILs rather than thermodynamic stabilization. The results are qualitatively similar to the trends found earlier for interactions between palladium clusters and components of imidazolium-based ILs, in spite of much larger contributions of the London dispersion forces to the binding of the [(Bu3PR)-P-t](+) cations to the cluster (up to 80%) relative to the case of 1-R-3-methylimidazolium cations (up to 40%).

International audience; The control of gene expression is a multi-layered process occurring at the level of DNA, RNA, and proteins. With the emergence of highly sensitive techniques, new aspects of RNA regulation have been uncovered leading to the emerging field of epitranscriptomics dealing with RNA modifications. Among those post-transcriptional modifications, N6-methyladenosine (m 6 A) is the most prevalent in messenger RNAs (mRNAs). This mark can either prevent or stimulate the formation of RNA-protein complexes, thereby influencing mRNA-related mechanisms and cellular processes. This review focuses on proteins containing a YTH domain (for YT521-B Homology), a small building block, that selectively detects the m 6 A nucleotide embedded within a consensus motif. Thereby, it contributes to the recruitment of various effectors involved in the control of mRNA fates through adjacent regions present in the different YTH-containing proteins.

International audience; Titratable residues determine the acid/base behavior of proteins, strongly influencing their function; in addition, proton binding is a valuable reporter on electrostatic interactions. We describe a method for pK(a) calculations, using constant-pH Monte Carlo (MC) simulations to explore the space of sidechain conformations and protonation states, with an efficient and accurate generalized Born model (GB) for the solvent effects. To overcome the many-body dependency of the GB model, we use a "Native Environment" approximation, whose accuracy is shown to be good. It allows the precalculation and storage of interactions between all sidechain pairs, a strategy borrowed from computational protein design, which makes the MC simulations themselves very fast. The method is tested for 12 proteins and 167 titratable sidechains. It gives an rms error of 1.1 pH units, similar to the trivial "Null" model. The only adjustable parameter is the protein dielectric constant. The best accuracy is achieved for values between 4 and 8, a range that is physically plausible for a protein interior. For sidechains with large pKa shifts, ≥2, the rms error is 1.6, compared to 2.5 with the Null model and 1.5 with the empirical PROPKA method.

International audience; Characterizing the variability within an ensemble of protein structures is a common requirement in structural biology and bioinformatics. With the increasing number of protein structures becoming available, there is a need for new tools capable of automating the structural comparison of large ensemble of structures. We present Protein Structural Statistics (PSS), a command-line program written in Perl for Unix-like environments, dedicated to the calculation of structural statistics for a set of proteins. PSS can perform multiple sequence alignments, structure superpositions, calculate Cartesian and dihedral coordinate statistics, and execute cluster analyses. An HTML report that contains a convenient summary of results with figures, tables, and hyperlinks can also be produced. PSS is a new tool providing an automated way to compare multiple structures. It integrates various types of structural analyses through an user-friendly and flexible interface, facilitating the access to powerful but more specialized programs. PSS is easy to modify and extend and is distributed under a free and open source license. The relevance of PSS is illustrated by examples of application to pertinent biological problems.

International audience; Free energy simulations are a powerful tool to study molecular recognition. The most rigorous variants can provide in depth understanding for a particular system, but are not suited for high throughput application to large libraries of compounds. Related, but less expensive methods are increasingly popular, including continuum electrostatic methods like PBSA (''Poisson-Boltzmann Surface Area'') and Linear Response or Linear Interaction Energy methods (LRA, LIE). Here, we review the theoretical background of these methods and provide a unified framework. We focus on the electrostatic contributions to the binding free energy, analyzing nonpolar contributions more briefly. The methods reviewed introduce a multi-step pathway for ligand unbinding, with distinct steps that uncharge the bound ligand, then recharge the unbound ligand. They assume that the system responds to the charging/uncharging in a linear way. With this approximation, the free energy can be described by its one or two first derivatives with respect to a progress variable. The methods can then be classified according to which states of the system are actually simulated and the number of free energy derivatives (one or two) that are employed. The analysis should help clarify the relations between several important free energy methods and the approximations they make. It can suggest new ways to test them, and provide routes for their improvement.

International audience; Di- and tri-phosphate nucleotides are essential cofactors for many proteins, usually in an Mg(2+) -bound form. Proteins like GTPases often detect the difference between NDP and NTP and respond by changing conformations. To study such complexes, simple, fixed charge force fields have been used, which allow long simulations and precise free energy calculations. The preference for NTP or NDP binding depends on many factors, including ligand structure and Mg(2+) coordination and the changes they undergo upon binding. Here, we use a simple force field to examine two Mg(2+) coordination modes for the unbound GDP and GTP: direct, or "Inner Sphere" (IS) coordination by one or more phosphate oxygens and indirect, "Outer Sphere" (OS) coordination involving one or more bridging waters. We compare GTP: and GDP:Mg binding with OS and IS coordination; combining the results with experimental data then indicates that GTP prefers the latter. We also examine different kinds of IS coordination and their sensitivity to a key force field parameter: the optimal Mg:oxygen van der Waals distance Rmin . Increasing Rmin improves the Mg:oxygen distances, the GTP: and GDP:Mg binding affinities, and the fraction of GTP:Mg with β + γ phosphate coordination, but does not improve or change the GTP/GDP affinity difference, which remains much larger than experiment. It has no effect on the free energy of GDP binding to a GTPase. © 2012 Wiley Periodicals, Inc. Copyright © 2012 Wiley Periodicals, Inc.

International audience; Multistate protein design explores side chain mutations, with the backbone allowed to sample a small, predetermined library of conformations. To achieve Boltzmann sampling of sequences and conformations, we use a hybrid Monte Carlo (MC) scheme: a trial hop between backbone models is followed by a short MC segment where side chain rotamers adjust to the new backbone, before applying a Metropolis-like acceptance test. The theoretical form and a practical approximation for the acceptance test are derived. We then compute backbone conformational free energies for two SH2 and SH3 proteins using different routes and protocols, and verify that for simple test problems, the free energy behaves like a state function, a hallmark of Boltzmann sampling.

International audience; Leucyl-tRNA synthetase catalyzes attachment of leucine amino acid to its cognate tRNA. During the second, aminoacetylation, step of the reaction, the leucyl moiety is transferred from leucyl-adenylate to the terminal A76 adenosine of tRNA. In this work, we have investigated the aminoacetylation step catalyzed by leucyl-tRNA synthase, using ab initio quantum chemical/molecular mechanical hybrid potentials in conjunction with reaction-path-location algorithms and molecular dynamics free energy simulations. We have modeled reaction mechanisms arising from both crystallographic studies and computational work. We invoke various groups as potential proton acceptors namely, the phosphate and leucyl amino groups of leucyl-adenylate, the A76 base of tRNA, and the Asp80 and Glu532 residues of the protein and consider both metal-assisted and metal-free reactions. Free energy calculations indicate that both the phosphate group of leucyl adenylate and Glu532 are not strong bases. This agrees with the results of the quantum chemical/molecular mechanical reaction path calculations which give high free energy barriers for the studied pathways involving these groups. A self-assisted mechanism with the leucyl amino group and Asp80 as proton acceptors is the most likely. Furthermore, in this mechanism the presence of a metal ion coordinated by the phosphate group and Glu532 strongly activates the reaction.

International audience; With the increasing number of protein structures available, there is a need for tools capable of automating the comparison of ensembles of structures, a common requirement in structural biology and bioinformatics. PSSweb is a web server for protein structural statistics. It takes as input an ensemble of PDB files of protein structures, performs a multiple sequence alignment and computes structural statistics for each position of the alignment. Different optional functionalities are proposed: structure superposition, Cartesian coordinate statistics, dihedral angle calculation and statistics, and a cluster analysis based on dihedral angles. An interactive report is generated, containing a summary of the results, tables, figures and 3D visualization of superposed structures. The server is available at http://pssweb.org.

International audience; The GTPase aIF5B is a universally conserved initiation factor that assists ribosome assembly. Crystal structures of its nucleotide complexes, X-ray(GTP) and X-ray(GDP), are similar in the nucleotide vicinity, but differ in the orientation of a distant domain IV. This has led to two, contradictory, mechanistic models. One postulates that X-ray(GTP) and X-ray(GDP) are, respectively, the active, "ON" and the inactive, "OFF" states; the other postulates that both structures are OFF, whereas the ON state is still uncharacterized. We study GTP/GDP binding using molecular dynamics and a continuum electrostatic free energy method. We predict that X-ray(GTP) has a ≈ 3 kcal/mol preference to bind GDP, apparently contradicting its assignment as ON. However, the preference arises mainly from a single, nearby residue from the switch 2 motif: Glu81, which becomes protonated upon GTP binding, with a free energy cost of about 4 kcal/mol. We then propose a different model, where Glu81 protonation/deprotonation defines the ON/OFF states. With this model, the X-ray(GTP):GTP complex, with its protonated Glu81, is ON, whereas X-ray(GTP):GDP is OFF. The model postulates that distant conformational changes such as domain IV rotation are "uncoupled" from GTP/GDP exchange and do not affect the relative GTP/GDP binding affinities. We analyze the model using a general thermodynamic framework for GTPases. It yields rather precise predictions for the nucleotide specificities of each state, and the state specificities of each nucleotide, which are roughly comparable to the homologues IF2 and aIF2, despite the lack of any conformational switching in the model. © 2012 Wiley Periodicals, Inc.

International audience; Proteins exist as a dynamic ensemble of interconverting substates, which defines their conformational energy landscapes. Recent work has indicated that mutations that shift the balance between conformational substates (CSs) are one of the main mechanisms by which proteins evolve new functions. In the present study, we probe this assertion by examining phenotypic protein adaptation to extreme conditions, using the allosteric tetrameric lactate dehydrogenase (LDH) from the hyperthermophilic bacterium Thermus thermophilus (Tt) as a model enzyme. In the presence of fructose 1, 6 bis-phosphate (FBP), allosteric LDHs catalyze the conversion of pyruvate to lactate with concomitant oxidation of nicotinamide adenine dinucleotide, reduced form (NADH). The catalysis involves a structural transition between a low-affinity inactive "T-state" and a high-affinity active "R-state" with bound FBP. During this structural transition, two important residues undergo changes in their side chain conformations. These are R171 and H188, which are involved in substrate and FBP binding, respectively. We designed two mutants of Tt-LDH with one ("1-Mut") and five ("5-Mut") mutations distant from the active site and characterized their catalytic, dynamical, and structural properties. In 1-Mut Tt-LDH, without FBP, the K(m)(Pyr) is reduced compared with that of the wild type, which is consistent with a complete shifting of the CS equilibrium of H188 to that observed in the R-state. By contrast, the CS populations of R171, k(cat) and protein stability are little changed. In 5-Mut Tt-LDH, without FBP, K(m)(Pyr) approaches the values it has with FBP and becomes almost temperature independent, k(cat) increases substantially, and the CS populations of R171 shift toward those of the R-state. These changes are accompanied by a decrease in protein stability at higher temperature, which is consistent with an increased flexibility at lower temperature. Together, these results show that the thermal properties of an enzyme can be strongly modified by only a few or even a single mutation, which serve to alter the equilibrium and, hence, the relative populations of functionally important native-state CSs, without changing the nature of the CSs themselves. They also provide insights into the types of mutational pathways by which protein adaptation to temperature is achieved.

International audience; Archaeal initiation factor 2 (aIF2) is a GTPase involved in protein biosynthesis. In its GTP-bound, "ON" conformation, it binds an initiator tRNA and carries it to the ribosome. In its GDP-bound, "OFF" conformation, it dissociates from tRNA. To improve our understanding of the role of each conformational state in the aIF2 "life cycle", we start from the state immediately after GTP hydrolysis, ON:GDP:P(i) (where P(i) is inorganic phosphate), and consider the possible next steps on the pathway to the OFF:GDP product. The first possibility is P(i) dissociation, leading to ON:GDP, which could then relax into OFF:GDP. We use molecular dynamics simulations to compute the P(i) dissociation free energy and show that dissociation is highly favorable. The second possibility is conformational relaxation into the OFF state before P(i) dissociation, to form OFF:GDP:P(i). We estimate the corresponding free energy approximately, 2 ± 3.5 kcal/mol, so that this is an uphill or weakly downhill process. A third possibility is relaxation into another conformation, neither ON nor OFF. Indeed, a third, "MIXED" conformation was seen recently in a crystal structure of the aIF2:GDP:P(i) complex. For this conformational state, P(i) dissociation is weakly unfavorable, in contrast to the ON and OFF states. From this, we will deduce that if the MIXED:GDP complex is not too unstable, the ON:GDP:P(i) → MIXED:GDP:P(i) transformation is a downhill process, which can occur spontaneously. This suggests that the MIXED state could be a functional intermediate.

Computational Protein Design (CPD) is a promising method for high throughput protein and ligand mutagenesis. Recently, we developed a CPD method that used a polar-hydrogen energy function for protein interactions and a Coulomb/Accessible Surface Area (CASA) model for solvent effects. We applied this method to engineer aspartyl-adenylate (AspAMP) specificity into Asparaginyl-tRNA synthetase (AsnRS), whose substrate is asparaginyl-adenylate (AsnAMP). Here, we implement a more accurate function, with an all-atom energy for protein interactions and a residue-pairwise generalized Born model for solvent effects. As a first test, we compute aminoacid affinities for several point mutants of Aspartyl-tRNA synthetase (AspRS) and Tyrosyl-tRNA synthetase and stability changes for three helical peptides and compare with experiment. As a second test, we readdress the problem of AsnRS aminoacid engineering. We compare three design criteria, which optimize the folding free-energy, the absolute AspAMP affinity, and the relative (AspAMP-AsnAMP) affinity. The sequences and conformations are improved with respect to our previous, polar-hydrogen/CASA study: For several designed complexes, the AspAMP carboxylate forms three interactions with a conserved arginine and a designed lysine, as in the active site of the AspRS:AspAMP complex. The conformations and interactions are well maintained in molecular dynamics simulations and the sequences have an inverted specificity, favoring AspAMP over AsnAMP. The method is not fully successful, since experimental measurements with the seven most promising sequences show that they do not catalyze at a detectable level the adenylation of Asp (or Asn) with ATP. This may be due to weak AspAMP binding and/or disruption of transition-state stabilization. © 2011 Wiley-Liss, Inc. 79 12 3448 3468 Cited By :12

Computational methods are needed to help characterize the structure and function of protein–protein complexes. To develop and improve such methods, standard test problems are essential. One important test is to identify experimental structures from among large sets of decoys. Here, a flexible docking procedure was used to produce such a large ensemble of decoy complexes. In addition to their use for structure prediction, they can serve as a proxy for the nonspecific, protein–protein complexes that occur transiently in the cell, which are hard to characterize experimentally, yet biochemically important. For 202 homodimers and 41 heterodimers with known X-ray structures, we produced an average of 1217 decoys each. The structures were characterized in detail. The decoys have rather large protein–protein interfaces, with at least 45 residue–residue contacts for every 100 contacts found in the experimental complex. They have limited intramonomer deformation and limited intermonomer steric conflicts. The decoys thoroughly sample each monomer's surface, with all the surface amino acids being part of at least one decoy interface. The decoys with the lowest intramonomer deformation were analyzed separately, as proxies for nonspecific protein–protein complexes. Their interfaces are less hydrophobic than the experimental ones, with an amino acid composition similar to the overall surface composition. They have a poorer shape complementarity and a weaker association energy, but are no more fragmented than the experimental interfaces, with 2.1 distinct patches of interacting residues on average, compared to 2.6 for the experimental interfaces. The decoys should be useful for testing and parameterizing docking methods and scoring functions; they are freely available as PDB files at http://biology.polytechnique.fr/decoys. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010

The RNA silencing pathway constitutes a defense mechanism highly conserved in eukaryotes, especially in plants, where the underlying working principle relies on the repressive action triggered by the intracellular presence of double-stranded RNAs. This immune system performs a post-transcriptional suppression of aberrant mRNAs or viral RNAs by small interfering RNAs that are directed towards their target in a sequencespecific manner. However, viruses have evolved strategies to escape from silencing surveillance while promoting their own replication. Several viruses encode suppressor proteins that interact with different elements of the RNA silencing pathway and block it. The different suppressors are not phylogenetically nor structurally related and also differ in their mechanism of action. Here, we adopt a model-driven forward-engineering approach to understand the evolution of suppressor proteins and, in particular, why viral suppressors preferentially target some components of the silencing pathway. We analyzed three strategies characterized by different design principles: replication in the absence of a suppressor, suppressors targeting the first protein component of the pathway and suppressors targeting the small interfering RNAs. Our results seed light into the question of whether a virus must opt for devoting more time into transcription or into translation and on which would be the optimal step of the silencing pathway to be targeted by suppressors. In addition, we discussed the evolutionary implications of such designing principles., This work was supported by the Spanish Ministerio de Ciencia e Innovación grants BFU2009- 06993 to S.F.E. and TIN-2006-12860 to A.J. And by FP6-NEST-043340 (BioModularH2), FP7- ICT-043338 (Bactocom), FP7-KBBE-212894 (Tarpol), the Structural Funds of the European Regional Development Fund, the ATIGE-Genopole and the Foundation pour la Recherche Medicale grants (all to A.J.). J.C, G.R. and A.J. also acknowledge the HPC-Europa program (RII3- CT-2003-506079). G.R. was supported by a graduate fellowship from the Generalitat Valenciana and an EMBO Short-term fellowship.

International audience; We review the current status of expression of heterologous systems for bioenergy and bioproduction in bacteria using a model-based approach. As an aim for synthetic biology, it requires mathematical models of genetic modules that could be characterized independently of their context. This fastens the design of metabolic circuits using a combinatorial design approach, where given pathways could be optimized for maximal bioproduction, while being nontoxic for the chassis. We show how recent characterization of genetic parts, such as promoters, RBS or sRNAs could be used to fine-tune the expression of individual genes to achieve that goal. We also present lists of enzymes that are used for bioproduction, enlarging such set of biological parts.

International audience; The correct amino acid sequence of E. coli isoleucyl-tRNA synthetase (IleRS) was established by means of peptide mapping by MALDI mass spectrometry, using a set of four endoproteases (trypsin, LysC, AspN and GluC). Thereafter, the active site of IleRS was mapped by affinity labeling with reactive analogs of the substrates. For the ATP binding site, the affinity labeling reagent was pyridoxal 5'-diphospho-5'-adenosine (ADP-PL), whereas periodate-oxidized tRNA(Ile), the 2',3'-dialdehyde derivative of tRNA(Ile) was used to label the binding site for the 3'-end of tRNA on the synthetase. Incubation of either reagent with IleRS resulted in a rapid loss of both the tRNA(Ile) aminoacylation and isoleucinedependent isotopic ATP-PPi exchange activities. The stoichiometries of IleRS labeling by ADP-PL or tRNA(Ile)ox corresponded to 1 mol of reagent incorporated per mol of enzyme. Altogether, the oxidized 3'-end of tRNA(Ile) and the pyridoxal moiety of the ATP analog ADP-PL react with the lysyl residues 601 and 604 of the consensus sequence (601)KMSKS(605). Identification of the binding site for L-isoleucine or for non cognate amino acids on E. coli IleRS was achieved by qualitative comparative labeling of the synthetase with bromomethyl ketone derivatives of L-isoleucine (IBMK) or of the non-cognate amino acids valine (VBMK), phenylalanine (FBMK) and norleucine (NleBMK). Labeling of the enzyme with IBMK resulted in a complete loss of isoleucine-dependent isotopic [(32)P]PPi-ATP exchange activity. VBMK, NleBMK and FBMK were also capable of abolishing the activity of IleRS, FBMK being the less efficient in inactivating the synthetase. Analysis by MALDI mass spectrometry designated cysteines-462 and -718 as the target residues of the substrate analog IBMK on E. coli IleRS, whereas VBMK, NleBMK and FBMK labeled in common His-394, His-478 and Cys-718. In addition, VBMK and NleBMK, which are chemically similar to IBMK, were found covalently bound to Cys-462, and VBMK was specifically attached to His-332 (or His-337) of the synthetase. The amino acid residues labeled by the substrate analogs are mainly distributed between three regions in the primary structure of E. coli IleRS: these are segments [325-394], [451-479] and [591-604]. In the 3-D structures of IleRS from T. thermophilus and S. aureus, the [325-394] stretch is part of the editing domain, while fragments [451-479] and [591-604] representing the isoleucine binding domain and the dinucleotide (or Rossmann) fold domain, respectively, are located in the catalytic core. His-332 of E. coli IleRS, that is strictly conserved among all the available IleRS sequences is located in the editing active site of the synthetase. It is proposed that His-332 of E. coli IleRS participates directly in hydrolysis, or helps to deprotonate the hydroxyl group of threonine at the hydrolytic site.

International audience; Methionyl-adenylate, the mixed carboxylic-phosphoric acid anhydride synthesized by methionyl-tRNA synthetase (MetRS) is capable of reacting with this synthetase or other proteins, by forming an isopeptide bond with the epsilon-NH2 group of lysyl residues. It is proposed that the mechanism for the in vitro methionylation of MetRS might be accounted for by the in situ covalent reaction of methionyl-adenylate with lysine side chains surrounding the active center of the enzyme, as well as by exchange of the label between donor and acceptor proteins. Following the incorporation of 7.0 +/- 0.5 mol of methionine per mol of a monomeric truncated methionyl-tRNA synthetase species, the enzymic activities of [P-32]PP1-ATP isotopic exchange and tRNA(Met) aminoacylation were lowered by 75% and more than 90%, respectively. The addition of tRNA(Met) protected the enzyme against inactivation and methionine incorporation. Matrix-assisted laser desorption-ionization mass spectrometry designated lysines-114, -132, -142 (or -147), -270, -282. -335, -362, -402, -439, -465, and -547 of truncated methionyl-tRNA synthetase as the target residues for covalent binding of methionine. These lysyl residues are distributed at the surface of the enzyme between three regions [114-150], [270-362], and [402-465], all of which were previously shown to be involved in catalysis or to be located in the binding sites of the three substrates, methionine, ATP, and tRNA.

International audience; The Tet repressor (TetR) mediates the most important mechanism of bacterial resistance against tetracycline (Tc) antibiotics. In the absence of Tc, TetR is tightly bound to its operator DNA; upon binding of Tc with an associated Mg(2+) ion, it dissociates from the DNA, allowing expression of the repressed genes. Its tight control by Tc makes TetR broadly useful in genetic engineering. The Tc binding site is over 20 A from the DNA, so the binding signal must propagate a long distance. We use molecular dynamics simulations and continuum electrostatic calculations to test two models of the allosteric mechanism. We simulate the TetR:DNA complex, the Tc-bound, "induced" TetR, and the transition pathway between them. The simulations support the model inferred previously from the crystal structures and reveal new details. When [Tc:Mg](+) binds, the Mg(2+) ion makes direct and water-mediated interactions with helix 8 of one TetR monomer and helix 6 of the other monomer, and helix 6 is pulled in towards the central core of the structure. Hydrophobic interactions with helix 6 then pull helix 4 in a pendulum motion, with a maximal displacement at its N-terminus: the DNA interface. The crystal structure of an additional TetR reported here corroborates this motion. The N-terminal residue of helix 4, Lys48, is highly conserved in DNA-binding regulatory proteins of the TetR class and makes the largest contribution of any amino acid to the TetR:DNA binding free energy. Thus, the conformational changes lead to a drastic reduction in the TetR:DNA binding affinity, allowing TetR to detach itself from the DNA. Tc plays the role of a specific Mg(2+) carrier, whereas the Mg(2+) ion itself makes key interactions that trigger the allosteric transition in the TetR:Tc complex.

International audience; Computational protein design will continue to improve as new implementations and parameterizations are explored. An automated protein design procedure is implemented and applied to the full redesign of 16 globular proteins. We combine established but simple ingredients: a molecular mechanics description of the protein where nonpolar hydrogens are implicit, a simple solvent model, a folded state where the backbone is fixed, and a tripeptide model of the unfolded state. Sequences are selected to optimize the folding free energy, using a simple heuristic algorithm to explore sequence and conformational space. We show that a balanced parametrization, obtained here and in our previous work, makes this procedure effective, despite the simplicity of the ingredients. Calculations were done using our Proteins @ Home distributed computing platform, with the help of several thousand volunteers. We describe the software implementation, the optimization of selected terms in the energy function, and the performance of the method. We allowed all amino acids to mutate except glycines, prolines, and cysteines. For 15 of the 16 test proteins, the scores of the computed sequences were comparable to those of natural homologues. Using the low energy computed sequences in a BLAST search of the SWISSPROT database, we could retrieve natural sequences for all protein families considered, with no high-ranking false-positives. The good stability of the designed sequences was supported by molecular dynamics simulations of selected sequences, which gave structures close to the experimental native structure.

International audience; Yor1p, a Saccharomyces cerevisiae plasma membrane ABC-transporter, is associated to oligomycin resistance and to rhodamine B transport. Here, by using the overexpressing strain Superyor [A. Decottignies, A.M. Grant, J.W. Nichols, H. de Wet, D.B. McIntosh, A. Goffeau, ATPase and multidrug transport activities of the overexpressed yeast ABC protein Yor1p, J. Biol. Chem. 273 (1998) 12612-12622], we show that Yor1p also confers resistance to rhodamine 6G and to doxorubicin. In addition, Yor1p protects cells, although weakly, against tetracycline, verapamil, eosin Y and ethidium bromide. The basal ATPase activity of the overexpressed form of Yor1p was studied in membrane preparations. This activity is quenched upon addition of micromolar amounts of vanadate. Vmax and Km values of approximately 0.8 s(-1) and 50+/-8 microM are measured. Mutations of essential residues in the nucleotide binding domain 2 reduces the activity to that measured with a Deltayor1 strain. ATP hydrolysis is strongly inhibited by the addition of potential substrates of the transporter. Covalent reaction of 8-azido-[alpha-(32)P]ATP with Yor1p is not sensitive to the presence of excess oligomycin. Thus, competition of the drug with ATP binding is unlikely. Finally, we inspect possible hypotheses accounting for substrate inhibition, rather than stimulation, of ATP hydrolysis by the membrane preparation.

Faithful genetic code translation requires that each aminoacyl-tRNA synthetase recognise its cognate amino acid ligand specifically. Aspartyl-tRNA synthetase (AspRS) distinguishes between its negatively-charged Asp substrate and two competitors, neutral Asn and di-negative succinate, using a complex network of electrostatic interactions. Here, we used molecular dynamics simulations and site-directed mutagenesis experiments to probe these interactions further. We attempt to decrease the Asp/Asn binding free energy difference via single, double and triple mutations that reduce the net positive charge in the active site of Escherichia coli AspRS. Earlier, Glutamine 199 was changed to a negatively-charged glutamate, giving a computed reduction in Asp affinity in good agreement with experiment. Here, Lysine 198 was changed to a neutral leucine; then, Lys198 and Gln199 were mutated simultaneously. Both mutants are predicted to have reduced Asp binding and improved Asn binding, but the changes are insufficient to overcome the initial, high specificity of the native enzyme, which retains a preference for Asp. Probing the aminoacyl-adenylation reaction through pyrophosphate exchange experiments, we found no detectable activity for the mutant enzymes, indicating weaker Asp binding and/or poorer transition state stabilization. The simulations show that the mutations' effect is partly offset by proton uptake by a nearby histidine. Therefore, we performed additional simulations where the nearby Histidines 448 and 449 were mutated to neutral or negative residues: (Lys198Leu, His448Gln, His449Gln), and (Lys198Leu, His448Glu, His449Gln). This led to unexpected conformational changes and loss of active site preorganization, suggesting that the AspRS active site has a limited structural tolerance for electrostatic modifications. The data give insights into the complex electrostatic network in the AspRS active site and illustrate the difficulty in engineering charged-to-neutral changes of the preferred ligand. Proteins 2008. © 2007 Wiley-Liss, Inc.

Initiation of translation in eukaryotes and in archaea involves eukaryotic/archaeal initiation factor (e/aIF)1 and the heterotrimeric initiation factor e/aIF2. In its GTP-bound form, e/aIF2 provides the initiation complex with Met–tRNA i Met . After recognition of the start codon by initiator tRNA, e/aIF1 leaves the complex. Finally, e/aIF2, now in a GDP-bound form, loses affinity for Met–tRNA i Met and dissociates from the ribosome. Here, we report a 3D structure of an aIF2 heterotrimer from the archeon Sulfolobus solfataricus obtained in the presence of GDP. Our report highlights how the two-switch regions involved in formation of the tRNA-binding site on subunit γ exchange conformational information with α and β. The zinc-binding domain of β lies close to the guanine nucleotide and directly contacts the switch 1 region. As a result, switch 1 adopts a not yet described conformation. Moreover, unexpectedly for a GDP-bound state, switch 2 has the “ON” conformation. The stability of these conformations is accounted for by a ligand, most probably a phosphate ion, bound near the nucleotide binding site. The structure suggests that this GDP–inorganic phosphate (Pi) bound state of aIF2 may be proficient for tRNA binding. Recently, it has been proposed that dissociation of eIF2 from the initiation complex is closely coupled to that of Pi from eIF2γ upon start codon recognition. The nucleotide state of aIF2 shown here is indicative of a similar mechanism in archaea. Finally, we consider the possibility that release of Pi takes place after e/aIF2γ has been informed of e/aIF1 dissociation by e/aIF2β.

Our project for iGEM 2006 consisted of designing a cellular vanillin biosensor. We used an EnvZ -E. coli strain as a chassis, and constructed two different devices: a sensor and an actuator, assembled using OmpR-P as a standardised mediator. The sensor device contained a computation- ally designed vanillin receptor and a synthetic two-component signal transduction protein (Trz). The receptor protein was based on a ribose-binding protein as scaffold. The Trz was built by fusion of the periplasmic and transmembrane domains of a Trg protein with an EnvZ kinase domain. When the receptor complex binds Trg, an allosteric motion is propagated to the cyto- plasmic EnvZ kinase domain, resulting in autophosphorylation and subsequent phosphate transfer to the OmpR transcription factor, which finally induces transcription of the ompC promoter. As actuator, we used a synthetic transcriptional circuit, which implements an OmpR-P band detector having GFP and RFP as an output. We designed this circuit using a synthetic promoter working as an AND gate, which is synergistically activated by cI and CRP. Our constructed Trg-EnvZ fusion and AND promoter will be very useful to future synthetic biology projects. 1 Aims of the project The main goal of this project was to introduce students to synthetic biology through the iGEM competition (1), while trying to do an interesting project that could generate parts to be reused by many future projects in synthetic biology. We chose to design a cellular biosensor having vanillin as an input and the expression of two reporters as output. To perform this, we constructed two devices: a sensor (based on a receptor and a two-component signal transductor) and an actuator (which received the intermedi- ate signal and expressed the reporters). We used a phos- phorylation mechanism from EnvZ as transmembranal pathway using OmpR-P as a standardised mediator (Fig. 1). To avoid interference with the natural pathways and have a better control of the kinase/phosphatase domain, we chose an EnvZ - E. coli strain as a chassis. The receptor domain of the sensor was inspired on the work by the group of Dr Hellinga, in which a sensor of

International audience; GEK1, an Arabidopsis thaliana gene product, was recently identified through its involvement in ethanol tolerance. Later, this protein was shown to display 26% strict identity with archaeal d-Tyr-tRNA(Tyr) deacylases. To determine whether it actually possessed deacylase activity, the product of the GEK1 open reading frame was expressed in Escherichia coli from a multi-copy plasmid. Purified GEK1 protein contains two zinc ions and proves to be a broad-specific, markedly active d-aminoacyl-tRNA deacylase in vitro. Moreover, GEK1 expression is capable of functionally compensating in E. coli for the absence of endogeneous d-Tyr- tRNA(Tyr) deacylase. Possible connections between exposure of plants to ethanol/acetaldehyde and misaminoacylation of tRNA by d-amino acids are considered.

International audience; Free energy simulations for electrostatic and charging processes in complex molecular systems encounter specific difficulties owing to the long-range, 1/r Coulomb interaction. To calculate the solvation free energy of a simple ion, it is essential to take into account the polarization of nearby solvent but also the electrostatic potential drop across the liquid-gas boundary, however distant. The latter does not exist in a simulation model based on periodic boundary conditions because there is no physical boundary to the system. An important consequence is that the reference value of the electrostatic potential is not an ion in a vacuum. Also, in an infinite system, the electrostatic potential felt by a perturbing charge is conditionally convergent and dependent on the choice of computational conventions. Furthermore, with Ewald lattice summation and tinfoil conducting boundary conditions, the charges experience a spurious shift in the potential that depends on the details of the simulation system such as the volume fraction occupied by the solvent. All these issues can be handled with established computational protocols, as reviewed here and illustrated for several small ions and three solvated proteins.

International audience; In this study, experimental far infrared measurements of L-serine, L-threonine, L-cysteine, and L-methionine are presented showing the spectra for the 1.0-13.0 pH range. In parallel, solid state DFT calculations were performed on the amino acid zwitterions in the crystalline form. We focused on the lowest frequency far infrared normal modes, which required the most precision and convergence of the calculations. Analysis of the computational results, which included the potential energy distribution of the vibrational modes, permitted a detailed and almost complete assignment of the experimental spectrum. In addition to characteristic signals of the two main acid-base couples, CO2H/CO2- and NH3+/NH2, specific side chain contributions for these amino acids, including CCO and CCS vibrational modes were analyzed. This study is in line with the growing application of FIR measurements to biomolecules. (c) 2015 Elsevier B.V. All rights reserved.

International audience; Sorting of cargoes in endosomes occurs through their concentration into sorting platforms, called microdomains, from which transport intermediates are formed. The WASH complex localizes to such endosomal microdomains and triggers localized branched actin nucleation by activating the Arp2/3 complex. These branched actin networks are required for both the lateral compartmentalization of endosome membranes into distinct microdomains and for the fission of transport intermediates from these sorting platforms. In this chapter, we provide experimental protocols to study these two aspects of WASH physiology. We first describe how to image the dynamic membrane tubules resulting from the defects of WASH-mediated fission. We then describe how to study quantitatively the microdomain localization of WASH in live and fixed cells. Since microdomains are below the resolution limit of conventional light-microscopy techniques, this required the development of specific image analysis pipelines, which are detailed. The guidelines presented in this chapter can apply to other endomembrane microdomains beyond WASH in order to increase our understanding of trafficking in molecular and quantitative terms.

International audience; Molecular recognition between the aminoacyl-tRNA synthetase enzymes and their cognate amino acid ligands is essential for the faithful translation of the genetic code. In aspartyl-tRNA synthetase (AspRS), the co-substrate ATP binds preferentially with three associated Mg2+ cations in an unusual, bent geometry. The Mg2+ cations play a structural role and are thought to also participate catalytically in the enzyme reaction. Co-binding of the ATP x Mg3(2+) complex was shown recently to increase the Asp/Asn binding free energy difference, indicating that amino acid discrimination is substrate-assisted. Here, we used molecular dynamics free energy simulations and continuum electrostatic calculations to resolve two related questions. First, we showed that if one of the Mg2+ cations is removed, the Asp/Asn binding specificity is strongly reduced. Second, we computed the relative stabilities of the three-cation complex and the 2-cation complexes. We found that the 3-cation complex is overwhelmingly favored at ordinary magnesium concentrations, so that the protein is protected against the 2-cation state. In the homologous LysRS, the 3-cation complex was also strongly favored, but the third cation did not affect Lys binding. In tRNA-bound AspRS, the single remaining Mg2+ cation strongly favored the Asp-adenylate substrate relative to Asn-adenylate. Thus, in addition to their structural and catalytic roles, the Mg2+ cations contribute to specificity in AspRS through long range electrostatic interactions with the Asp side chain in both the pre- and post-adenylation states.

International audience; Tetracycline (Tc) is an important antibiotic, which binds specifically to the ribosome and several proteins, in the form of a Tc-:Mg2+ complex. To model Tc:protein and Tc:RNA interactions, we have developed a molecular mechanics force field model of Tc, which is consistent with the CHARMM force field for proteins and nucleic acids. We used structures from the Cambridge Crystallographic Data Base to identify the main Tc conformations that are likely to be present in solution and in biomolecular complexes. A conformational search was also done, using the MM3 force field to perform simulated annealing of Tc. Several resulting, low-energy structures were optimized with an ab initio model and used in developing the new Tc force field. Atomic charges and Lennard-Jones parameters were derived from a supermolecule ab initio approach. We considered the ab initio energies and geometries of a probe water molecule interacting with Tc at 36 different positions. We considered both a neutral and a zwitterionic Tc form, with and without bound Mg2+. The final rms deviation between the ab initio and force field energies, averaged over all forms, was just 0.35 kcal/mol. The model also reproduces the ab initio geometry and flexibility of Tc. As further tests, we did simulations of a Tc crystal, of Tc:Mg2+ and Tc:Ca2+ complexes in aqueous solution, and of a solvated complex between Tc:Mg2+ and the Tet repressor protein (TetR). With slight, ad hoc adjustments, the model can reproduce the experimental, relative, Tc binding affinities of Mg2+ and Ca2+. It performs well for the structure and fluctuations of the Tc:Mg2+:TetR complex. The model should therefore be suitable to investigate the interactions of Tc with proteins and RNA. It provides a starting point to parameterize other compounds in the large Tc family.

Introduction The production of antimicrobial peptides by airway epithelial cells is an important component of the innate immune response to pulmonary infection and inflammation. Hepcidin is a β-defensin-like antimicrobial peptide and acts as a principal iron regulatory hormone. Hepcidin is mostly produced by hepatocytes, but is also expressed by other cells, such as airway epithelial cells. However, nothing is known about its function in lung infections and inflammatory diseases. We therefore sought to investigate the role of airway epithelial cell-derived hepcidin in sepsis-induced acute lung injury. Methods Acute lung injury was induced by polymicrobial sepsis via cecal ligation and puncture (CLP) surgery. Adenovirus-mediated short hairpin RNA specific for the mouse hepcidin gene hepc1 and control adenovirus were intratracheally injected into mice. The adenovirus-mediated knockdown of hepcidin in airway epithelial cells was evaluated in vivo. Lung injury and the seven-day survival rate were assessed. The levels of hepcidin-related iron export protein ferroportin were measured, and the iron content and function of alveolar macrophages were evaluated. Results The hepcidin level in airway epithelial cells was upregulated during polymicrobial sepsis. The knockdown of airway epithelial cell-derived hepcidin aggravated the polymicrobial sepsis-induced lung injury and pulmonary bacterial infection and increased mortality (53.33% in Ad-shHepc1-treated mice versus 12.5% in Ad-shNeg-treated mice, P

We describe an efficient generalized Born (GB) approximation for proteins, in which the interaction energy between two amino acids depends on the whole protein structure, but can be accurately computed from residue-pairwise information. Two results make the scheme pairwise. First, an accurate expression exists for the interaction energy between two residues R and R′ that depends on the product B = BRBR′ of their residue Born solvation radii. Second, this expression is accurately fitted by a parabolic function of B the (three) fitting coefficients depend only on the pair RR′, not on its environment. In effect, the quantity B captures all the information that is relevant about the pair's dielectric environment. The method is tested with calculations on several hundred structures of the proteins trpcage, BPTI, ubiqutin, and thoredoxin. It yields solvation energies in better agreement with Poisson calculations than a traditional GB formulation. We also compute the effect of the protein/solvent environment on the interactions between pairs of charged residues in the active site of the enzyme aspartyl-tRNA synthetase. Our method captures this effect as accurately as traditional GB. Because it is residue-pairwise, the method can be incorporated into efficient protocols for rotamer placement and computational protein design. © 2005 American Chemical Society. 109 47 22667 22673 Cited By :26

Proton binding plays a critical role in protein structure and function. We report pKa calculations for three aspartates in two proteins, using a linear response approach, as well as a "standard" Poisson-Boltzmann approach. Averaging over conformations from the two endpoints of the proton-binding reaction, the protein's atomic degrees of freedom are explicitly modeled. Treating macroscopically the protein's electronic polarizability and the solvent, a meaningful model is obtained, without adjustable parameters. It reproduces qualitatively the electrostatic potentials, proton-binding free energies, Marcus reorganization free energies, and pKa shifts from explicit solvent molecular dynamics simulations, and the pKa shifts from experiment. For thioredoxin Asp-26, which has a large pKa upshift, we correctly capture the balance between unfavorable carboxylate desolvation and favorable interactions with a nearby lysine similarly for RNase A Asp-14, which has a large pKa downshift. For the unshifted thioredoxin Asp-20, desolvation by the protein cavity is overestimated by 2.9 pKa units several effects could explain this. "Standard" Poisson-Boltzmann methods sidestep this problem by using a large, ad hoc protein dielectric but protein charge-charge interactions are then incorrectly downscaled, giving an unbalanced description of the reaction and a large error for the shifted pKa values of Asp-26 and Asp-14. © 2005 by the Biophysical Society. 88 6 3888 3904 Cited By :56

International audience; Eukaryotic and archaeal initiation factor 2 (e- and aIF2, respectively) are heterotrimeric proteins (alphabetagamma) supplying the small subunit of the ribosome with methionylated initiator tRNA. The gamma subunit forms the core of the heterotrimer. It resembles elongation factor EF1-A and ensures interaction with Met-tRNA(i)(Met). In the presence of the alpha subunit, which is composed of three domains, the gamma subunit expresses full tRNA binding capacity. This study reports the crystallographic structure of the intact aIF2alpha subunit from the archaeon Pyrococcus abyssi and that of a derived C-terminal fragment containing domains 2 and 3. The obtained structures are compared with those of N-terminal domains 1 and 2 of yeast and human eIF2alpha and with the recently determined NMR structure of human eIF2alpha. We show that the three-domain organization in the alpha subunit is conserved in archaea and eukarya. Domains 1 and 2 form a rigid body linked to a mobile third domain. Sequence comparisons establish that the most conserved regions in the aIF2alpha polypeptide lie at opposite sides of the protein, within domain 1 and domain 3, respectively. These two domains are known to exhibit RNA binding capacities. We propose that domain 3, which is known to glue the alpha subunit onto the gamma subunit, participates in Met-tRNA(i)(Met) binding while domain 1 recognizes either rRNA or mRNA on the ribosome. Thereby, the observed structural mobility within the e- and aIF2alpha molecules would be an integral part of the biological function of this subunit in the heterotrimeric e- and aIF2alphabetagamma factors.