Search

International audience; Intrinsically disordered proteins are ubiquitous throughout all known proteomes, playing essential roles in all aspects of cellular and extracellular biochemistry. To understand their function, it is necessary to determine their structural and dynamic behavior and to describe the physical chemistry of their interaction trajectories. Nuclear magnetic resonance is perfectly adapted to this task, providing ensemble averaged structural and dynamic parameters that report on each assigned resonance in the molecule, unveiling otherwise inaccessible insight into the reaction kinetics and thermodynamics that are essential for function. In this review, we describe recent applications of NMR-based approaches to understanding the conformational energy landscape, the nature and time scales of local and long-range dynamics and how they depend on the environment, even in the cell. Finally, we illustrate the ability of NMR to uncover the mechanistic basis of functional disordered molecular assemblies that are important for human health.

International audience; Intercalation into DNA is the interaction mode of some anthracycline antibiotics. Recently, the molecular mechanism of this process was explored using the static free energy landscape. Here we explore the dynamical effects in the intercalation of proflavine into DNA by calculating the transmission coefficient κ-providing a measure of the departure from transition state theory for the reaction rate constant-by examination of the recrossing events at the transition state. For that purpose, we first found the accurate transition state of this complex system-as judged by a committor analysis-using a set of all-atom simulations of total length 6.3 ms. In a subsequent calculation of the transmission coefficient κ in another extensive set of simulations the small value κ = 0.1 was found, indicating a significant departure from TST. Comparison of this result with Grote-Hynes and Kramers theories shows that neither theory is able to capture this complex system's recrossing events; the source of this striking failure is discussed, as are related aspects of the mechanism. This study suggests that, for biomolecular processes similar to this, dynamical effects essential for the process are complex in nature and require novel approaches for their elucidation.

13 págs., 6 figs., 3 tabs., Nucleoside 2′-deoxyribosyltransferases (NDTs) catalyze the cleavage of glycosidic bonds of 2′-deoxynucleosides and the following transfer of the 2′-deoxyribose moiety to acceptor nucleobases. Here, we report the crystal structures and biochemical properties of the first tetrameric NDTs: the type I NDT from the mesophilic bacterium Enterococcus faecalis V583 (EfPDT) and the type II NDT from the bacterium Desulfotalea psychrophila (DpNDT), the first psychrophilic NDT. This novel structural and biochemical data permitted an exhaustive comparative analysis aimed to shed light into the basis of the high global stability of the psychrophilic DpNDT, which has a higher melting temperature than EfPDT (58.5 °C versus 54.4 °C) or other mesophilic NDTs. DpNDT possesses a combination of unusual structural motifs not present neither in EfPDT nor any other NDT that most probably contribute to its global stability, in particular, a large aliphatic isoleucine-leucine-valine (ILV) bundle accompanied by a vicinal disulfide bridge and also an intersubunit disulfide bridge, the first described for an NDT. The functional and structural features of DpNDT do not fit the standard features of psychrophilic enzymes, which lead us to consider the implication of (sub)cellular levels together with the protein level in the adaptation of enzymatic activity to low temperatures., This work was supported by grants CTQ2009-11543 (IDLM), PID2020-117025RB-I00 (JF-L) and AGL2017-84614-C2-2-R (JMM) from the Spanish Ministry of Science and Innovation, S2009/PPQ-1752 (IDLM) from Comunidad Autonoma de Madrid and XSAN192006 (IDLM) from the Santander Foundation. This research was funded in part by National Institutes of Health grant (GM115586 to AJ). The use of SBC beamlines at the Advanced Photon Source is supported by the U.S. Department of Energy (DOE) Office of Science and operated for the DOE Office of Science by Argonne National Laboratory under Contract No. DE-AC02-06CH11357 (AJ).

International audience; Microbial transition state theory (MTS) offers a theoretically explicit mathematical model for substrate limited microbial growth. By considering a first order approximation of the MTS equation one recovers the well-known Monod’s expression for growth, which was regarded as a purely empirical function. The harvest volume of a cell as defined in MTS theory can then be related to the affinity concept, giving a new physical interpretation to it, and a new way to determine its value. Consequences of such a relationship are discussed.

6 p.-5 fig., Protein (un)folding rates depend on the free-energy barrier separating the native and unfolded states and a prefactor term, which sets the timescale for crossing such barrier or folding speed limit. Because extricating these two factors is usually unfeasible, it has been common to assume a constant prefactor and assign all rate variability to the barrier. However, theory and simulations postulate a protein-specific prefactor that contains key mechanistic information. Here, we exploit the special properties of fast-folding proteins to experimentally resolve the folding rate prefactor and investigate how much it varies among structural homologs. We measure the ultrafast (un)folding kinetics of five natural WW domains using nanosecond laser-induced temperature jumps. All five WW domains fold in microseconds, but with a 10-fold difference between fastest and slowest. Interestingly, they all produce biphasic kinetics in which the slower phase corresponds to reequilibration over the small barrier (, This research has been funded by Grant ERC-2012-ADG-323059 from the European Research Council and Grant NSF-MCB-1616759 from the National Science Foundation. V.M. also acknowledges support from the W. M. Keck Foundation and the Center for Cellular and Biomolecular Machines at University of California, Merced (Grant NSF-CREST-1547848). I.L. acknowledges support from the Spanish Ministry of Science and Education through Grants BIO2012-39922-CO2-01 and BIO2016-787746-C2-1-R, and the European Union through Fonds Européen de Développement Économique et Régional.

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; Cryo-electron microscopy is rapidly emerging as a powerful technique to determine the structures of complex macromolecular systems elusive to other techniques. Because many of these systems are highly dynamical, characterizing their movements is also a crucial step to unravel their biological functions. To achieve this goal, we report an integrative modeling approach to simultaneously determine structure and dynamics of macromolecular systems from cryo-electron microscopy density maps. By quantifying the level of noise in the data and dealing with their ensemble-averaged nature, this approach enables the integration of multiple sources of information to model ensembles of structures and infer their populations. We illustrate the method by characterizing structure and dynamics of the integral membrane receptor STRA6, thus providing insights into the mechanisms by which it interacts with retinol binding protein and translocates retinol across the membrane.

International audience; Hearing relies on the transduction of sound-evoked vibrations into electric signals, occurring in the stereocilia bundle of hair cells. The bundle is organized in a staircase pattern formed by rows of packed stereocilia. This architecture is pivotal to transduction and involves a network of scaffolding proteins with hitherto uncharacterized features. Key interactions in this network are mediated by PDZ domains. Here, we describe the architecture of the first two PDZ domains of whirlin, a protein involved in these assemblies and associated with congenital deaf-blindness. C-terminal hairpin extensions of the PDZ domains mediate the transient supramodular assembly, which improves the binding capacity of the first domain. We determined a detailed structural model of the closed conformation of the PDZ tandem and characterized its equilibrium with an ensemble of open conformations. The structural and dynamic behavior of this PDZ tandem provides key insights into the regulatory mechanisms involved in the hearing machinery.

PorM is a membrane protein that is involved in the assembly of the type IX secretion system (T9SS) inPorphyromonas gingivalis, a major bacterial pathogen that is responsible for periodontal disease in humans. In the context of structural studies of PorM to better understand T9SS assembly, four camelid nanobodies were selected, produced and purified, and their specific interaction with the N-terminal or C-terminal part of the periplasmic domain of PorM was investigated. Diffracting crystals were also obtained, and the structures of the four nanobodies were solved by molecular replacement. Furthermore, two nanobodies were used as crystallization chaperones and turned out to be valuable tools in the structure-determination process of the periplasmic domain of PorM.

International audience; When CnrX, the periplasmic sensor protein in the CnrYXH transmembrane signal transduction complex of Cupriavidus metallidurans CH34, binds the cognate metal ions Ni(II) or Co(II), the ECF-type sigma factor CnrH is made available in the cytoplasm for the RNA-polymerase to initiate transcription at the cnrYp and cnrCp promoters. Ni(II) or Co(II) are sensed by a metal-binding site with a N3O2S coordination sphere with octahedral geometry, where S stands for the thioether sulfur of the only methionine (Met123) residue of CnrX. The M123A-CnrX derivative has dramatically reduced signal propagation in response to metal sensing while the X-ray structure of Ni-bound M123A-CnrXs showed that the metal-binding site was not affected by the mutation. Ni(II) remained six-coordinate in M123A-CnrXs, with a water molecule replacing the sulfur as the sixth ligand. H32A-CnrXs, the soluble model of the wild-type membrane-anchored CnrX, was compared to the double mutants H32A-M123A-CnrXs and H32A-M123C-CnrXs to spectroscopically evaluate the role of this unique ligand in the binding site of Ni or Co. The Co- and Ni-bound forms of the protein display unusually blue-shifted visible spectra. TD-DFT calculations using structure-based models allowed identification and assignment of the electronic transitions of Co-bound form of the protein and its M123A derivative. Among them, the signature of the S-Co transition is distinguishable in the shoulder at 530 nm. In vitro affinity measurements point out the crucial role of Met123 in the selectivity for Ni or Co, and in vivo data support the conclusion that Met123 is a trigger of the signal transduction.

International audience; This work rationalizes the glucuronidation process (one of the reactions of the phase II metabolism) for drugs having a carboxylic acid moiety. At this stage, acylglucuronides (AG) metabolites are produced, that have largely been reported in the literature for various drugs (e.g., mycophenolic acid (MPA), diclofenac, ibuprofen, phenylacetic acids). The competition between migration and hydrolysis is rationalized by adequate quantum calculations, combing MP2 and density functional theory (DFT) methods. At the molecular scale, the former process is a real rotation of the drug around the glucuconic acid. This chemical-engine provides four different metabolites with various toxicities. Migration definitely appears feasible under alkaline conditions, making proton release from the OH groups. The latter reaction (hydrolysis) releases the free drug, so the competition is of crucial importance to tackle drug action and elimination. From the theoretical data, both migration and hydrolysis appear kinetically and thermodynamically favored, respectively.

[Departement_IRSTEA]Ecotechnologies [TR1_IRSTEA]TED [TR2_IRSTEA]INSPIRE; International audience; Although mathematical modelling has reached a degree of maturity in the last decades, microbial ecology is still developing, albeit at a rapid pace thanks to new insights provided by modern molecular tools. However, whilst microbiologists have long enjoyed the perspectives that particular mathematical frameworks can provide, there remains a reluctance to fully embrace the potential of models, which appear too complex, esoteric or distant from the “real-world”. Nevertheless there is a strong case for pursuing the development of mathematical models to describe microbial behaviour and interactions, dynamically, spatially and across scales. Here we put forward perspectives on the current state of mathematical modelling in microbial ecology, looking back at the developments that have defined the synergies between the disciplines, and outline some of the existing challenges that motivate us to provide practical models in the hope that greater engagement with empiricists and practitioners in the microbiological domain may be achieved. We also indicate recent advances in modelling that have had impact in both the fundamental understanding of microbial ecology and its practical application in engineered biological systems. In this way, it is anticipated that interest can be garnered from across the microbiological spectrum resulting in a broader uptake of mathematical concepts in lecture theatres, laboratories and industrial systems.

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; The hematopoietic colony stimulating factor-1 receptor (CSF-1R or FMS) is essential for the cellular repertoire of the mammalian immune system. Here, we report a structural and mechanistic consensus for the assembly of human and mouse CSF-1:CSF-1R complexes. The EM structure of the complete extracellular assembly of the human CSF-1:CSF-1R complex reveals how receptor dimerization by CSF-1 invokes a ternary complex featuring extensive homotypic receptor contacts and striking structural plasticity at the extremities of the complex. Studies by small-angle X-ray scattering of unliganded hCSF-1R point to large domain rearrangements upon CSF-1 binding, and provide structural evidence for the relevance of receptor predimerization at the cell surface. Comparative structural and binding studies aiming to dissect the assembly principles of human and mouse CSF-1R complexes, including a quantification of the CSF-1/CSF-1R species cross-reactivity, show that bivalent cytokine binding to receptor coupled to ensuing receptor-receptor interactions are common denominators in extracellular complex formation.

International audience; This paper reports on the properties of bilayers composed of dioctadecyldimethylammonium bromide (DODAB) and oleic acid (OA) at various molar ratios. The mole fraction of OA, X(OA), was varied in the range of 0-1 and the total lipid content was constant and equal to 10 mM. The DODAB/OA dispersions were extruded at a temperature higher than that of the gel-liquid transition of DODAB. The morphology of bilayer structures formed in the dispersions was inspected using a cryogenic transmission electron microscopy (cryo-TEM) and a differential interference contrast microscopy (DIC). The observations revealed that the incorporation of OA into DODAB bilayer results in a decrease of the membrane curvature. Anisotropy measurements using 1,6-diphenylhexatriene (DPH) as a rotator probe demonstrated that the DODAB/OA membrane microviscosity decreased considerably for X(OA)>0.4. The thermal behavior of DODAB/OA membranes has been studied by differential scanning calorimetry (DSC). In the case of the systems in which X(OA)

International audience; The multimeric DNA sliding clamps confer high processivity to replicative DNA polymerases and are also binding platforms for various enzymes involved in DNA metabolism. These enzymes interact with the clamp through a small peptide that binds into a hydrophobic pocket which is a potential target for the development of new antibacterial compounds. Starting from a generic heptapeptide, we used a structure-based strategy to improve the design of new peptide ligands. Chemical modifications at specific residues result in a dramatic increase of the interaction as measured by SPR and ITC. The affinity of our best hits was improved by 2 orders of magnitude as compared to the natural ligand, reaching 10(-8) M range. The molecular basis of the interactions was analyzed by solving the co-crystal structures of the most relevant peptides bound to the clamp and reveals how chemical modifications establish new contacts and contributes to an increased affinity of the ligand.

International audience; The fast reaction of (*)NO and O(2)(*-) to give ONOO(-) has been extensively studied at irreversible conditions, but the reasons for the wide variations in observed forward rate constants (3.8 ≤ k(f) ≤ 20 × 10(9) M(-1) s(-1)) remain unexplained. We characterized the diffusion-dependent aqueous (pH > 12) chemical equilibrium of the form (*)NO + O(2)(*-) = ONOO(-) with respect to its dependence on temperature, viscosity, and [ONOO(-)](eq) by determining [ONOO(-)](eq) and [(*)NO](eq). The equilibrium forward reaction rate constant (k(f)(eq)) has negative activation energy, in contrast to that found under irreversible conditions. In contradiction to the law of mass action, we demonstrate that the equilibrium constant depends on ONOO(-) concentration. Therefore, a wide range of k(f)(eq) values could be derived (7.5-21 × 10(9) M(-1) s(-1)). Of general interest, the variations in k(f) can thus be explained by its dependence on the distance between ONOO(-) particles (sites of generation of (*)NO and O(2)(*-)).

International audience; The theoretical and experimental description of fluid phase endocytosis (FPE) requires an asymmetry in phospholipid number between the two leaflets of the cell membrane, which provides the biomechanical torque needed to generate membrane budding. Although the motor force behind FPE is defined, its kinetic has yet to be determined. Based on a body of evidences suggesting that the mean surface tension is unlikely to be involved in endocytosis we decided to determine whether the cytosolic hydrostatic pressure could be involved, by considering a constant energy exchanged between the cytosol and the cell membrane. The theory is compared to existing experimental data obtained from FPE kinetic studies in living cells where altered phospholipid asymmetry or changes in the extracellular osmotic pressure have been investigated. The model demonstrates that FPE is dependent on the influx and efflux of vesicular volumes (i.e. vesicular volumes recycling) rather than the membrane tension of cells. We conclude that: (i) a relationship exists between membrane lipid number asymmetry and resting cytosolic pressure and (ii) the validity of Laplace's law is limited to cells incubated in a definite hypotonic regime. Finally, we discuss how the model could help clarifying elusive observations obtained from different fields and including: (a) the non-canonical shuttling of aquaporin in cells, (b) the relationship between high blood pressure and inflammation and (c) the mechanosensitivity of the sodium/proton exchanger.

International audience; The assembly of cytochrome c oxidase (CcO) in yeast mitochondria is dependent on a new assembly factor designated Coa2. Coa2 was identified from its ability to suppress the respiratory deficiency of coa1Delta and shy1Delta cells. Coa1 and Shy1 function at an early step in maturation of the Cox1 subunit of CcO. Coa2 functions downstream of the Mss51-Coa1 step in Cox1 maturation and likely concurrent with the Shy1-related heme a(3) insertion into Cox1. Coa2 interacts with Shy1. Cells lacking Coa2 show a rapid degradation of newly synthesized Cox1. Rapid Cox1 proteolysis also occurs in shy1Delta cells, suggesting that in the absence of Coa2 or Shy1, Cox1 forms an unstable conformer. Overexpression of Cox10 or Cox5a and Cox6 or attenuation of the proteolytic activity of the m-AAA protease partially restores respiration in coa2Delta cells. The matrix-localized Coa2 protein may aid in stabilizing an early Cox1 intermediate containing the nuclear subunits Cox5a and Cox6.

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; TEM-1 is a class A beta-lactamase that contributes to the primary defensive measure used by bacteria to hydrolyze the clinically-relevant beta-lactam antibiotics. Several crystal structures of this enzyme complexed with inhibitors display the active-site residue Tyr105 in an alternate orientation relative to that assigned in the free or in the substrate-bound forms. Thus, the alternate conformation may not be favored in the free enzyme and may be adopted only in the presence of inhibitor. As the residue at position 105 is a determinant of substrate specificity, we sought a better understanding of the relation between its conformation and its function in ligand binding. Here, we perform a molecular dynamics simulated annealing protocol to identify stable orientations adopted by Tyr105 in free TEM-1. Our results demonstrate that, in the absence of substrate, structurally validated conformers of Tyr105 predominantly adopt either of the two rotameric orientations observed in the crystal structures. This suggests that adoption of either conformation in the free enzyme is energetically favored and is not strictly promoted by ligand binding. We propose that free TEM-1 alternates between these two conformations of Tyr105 and that a dynamically heterogeneous population of both rotamers exists in solution. The conformational change significantly reshapes the active-site cavity and modifies the potential for forming specific ligand contacts. Our results add to the body of evidence suggesting that Tyr105 displays a dynamical behavior resulting in alternate ligand binding modes and are consistent with the lower affinity of TEM-1 for cephalosporins relative to penicillins.

The molecular mechanics Poisson–Boltzmann surface area (MM/PBSA) method was applied to the study of the protein–protein complex between a camelid single chain variable domain (cAb-Lys3) and hen egg white lysozyme (HEL), and between cAb-Lys3 and turkey egg white lysozyme (TEL). The electrostatic energy was estimated by solving the linear Poisson–Boltzmann equation. A free energy decomposition scheme was developed to determine binding energy hot spots of each complex. The calculations identified amino acids of the antibody that make important contributions to the interaction with lysozyme. They further showed the influence of small structural variations on the energetics of binding and they showed that the antibody amino acids that make up the hot spots are organized in such a way as to mimic the lysozyme substrate. Through further analysis of the results, we define the concept of “efficient amino acids,” which can provide an assessment of the binding potential of a particular hot spot interaction. This information, in turn, can be useful in the rational design of small molecules that mimic the antibody. The implications of using free energy decomposition to identify regions of a protein–protein complex that could be targeted by small molecules inhibitors are discussed. Proteins 2007. © 2007 Wiley-Liss, Inc.

International audience; The activity of a neural network is defined by patterns of spiking and silence from the individual neurons. Because spikes are (relatively) sparse, patterns of activity with increasing numbers of spikes are less probable, but, with more spikes, the number of possible patterns increases. This tradeoff between probability and numerosity is mathematically equivalent to the relationship between entropy and energy in statistical physics. We construct this relationship for populations of up to N = 160 neurons in a small patch of the vertebrate retina, using a combination of direct and model-based analyses of experiments on the response of this network to naturalistic movies. We see signs of a thermodynamic limit, where the entropy per neuron approaches a smooth function of the energy per neuron as N increases. The form of this function corresponds to the distribution of activity being poised near an unusual kind of critical point. We suggest further tests of criticality, and give a brief discussion of its functional significance.

International audience; Recent experimental results based on multielectrode and imaging techniques have reinvigorated the idea that large neural networks operate near a critical point, between order and disorder. However, evidence for criticality has relied on the definition of arbitrary order parameters, or on models that do not address the dynamical nature of network activity. Here we introduce a novel approach to assess criticality that overcomes these limitations, while encompassing and generalizing previous criteria. We find a simple model to describe the global activity of large populations of ganglion cells in the rat retina, and show that their statistics are poised near a critical point. Taking into account the temporal dynamics of the activity greatly enhances the evidence for criticality, revealing it where previous methods would not. The approach is general and could be used in other biological networks.

International audience; The effect of increasing concentrations of 2,2,2-trifluoroethanol (TFE) on the conformational stability of the Shiga toxin B-subunit (STxB), a bacterial homopentameric protein involved in cell-surface binding and intracellular transport, has been studied by far-, near-UV circular dichroism (CD), intrinsic fluorescence, analytical ultracentrifugation, and differential scanning calorimetry (DSC) under equilibrium conditions. Our data show that the native structure of STxB is highly perturbed by the presence of TFE. In fact, at concentrations of TFE above 20% (v/v), the native pentameric conformation of the protein is cooperatively transformed into a helix-rich monomeric and partially folded conformational state with no significant tertiary structure. Additionally, no cooperative transition was detected upon a further increase in the TFE concentration (above 40% (v/v)). The thermal stability of STxB was investigated at several different TFE concentrations using DSC and CD spectroscopy. Thermal transitions at TFE concentrations of up to 20% (v/v) were successfully fitted to the two-state folding/unfolding coupled to oligomerization model consistent with the transition between a pentameric folded conformation to a monomeric state of the protein, which the presence of TFE stabilizes as a partially folded conformation.

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.

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

The HIV-1 type dimerization initiation signal (DIS) loop was used as a starting point for the analysis of the stability of Watson-Crick (WC) base pairs in a tertiary structure context. We used ultraviolet melting to determine thermodynamic parameters for loop-loop tertiary interactions and compared them with regular secondary structure RNA helices of the same sequences. In 1 M Na+ the loop-loop interaction of a HIV-1 DIS type pairing is 4 kcal/mol more stable than its sequence in an equivalent regular and isolated RNA helix. This difference is constant and sequence independent, suggesting that the rules governing the stability of WC base pairs in the secondary structure context are also valid for WC base pairs in the tertiary structure context. Moreover, the effect of ion concentration on the stability of loop-loop tertiary interactions differs considerably from that of regular RNA helices. The stabilization by Na+ and Mg2+ is significantly greater if the base pairing occurs within the context of a loop-loop interaction. The dependence of the structural stability on salt concentration was defined via the slope of a T(m)/log [ion] plot. The short base-paired helices are stabilized by 8 degrees C/log [Mg2+] or 11 degrees C/log [Na+], whereas base-paired helices forming tertiary loop-loop interactions are stabilized by 16 degrees C/log [Mg2+] and 26 degrees C/log [Na+]. The different dependence on ionic strength that is observed might reflect the contribution of specific divalent ion binding to the preformation of the hairpin loops poised for the tertiary kissing loop-loop contacts.

International audience; The RecQ helicase family is highly conserved from bacteria to men and plays a conserved role in the preservation of genome integrity. Its deficiency in human cells leads to a marked genomic instability that is associated with premature aging and cancer. To determine the thermodynamic parameters for the interaction of Escherichia coli RecQ helicase with DNA, equilibrium binding studies have been performed using the thermodynamic rigorous fluorescence titration technique. Steady-state fluorescence anisotropy measurements of fluorescein-labeled oligonucleotides revealed that RecQ helicase bound to DNA with an apparent binding stoichiometry of 1 protein monomer/10 nucleotides. This stoichiometry was not altered in the presence of AMPPNP (adenosine 5'-(beta,gamma-imido) triphosphate) or ADP. Analyses of RecQ helicase interactions with oligonucleotides of different lengths over a wide range of pH, NaCl, and nucleic acid concentrations indicate that the RecQ helicase has a single strong DNA binding site with an association constant at 25 degrees C of K=6.7 +/- 0.95 x 10(6) M(-1) and a cooperativity parameter of omega=25.5 +/- 1.2. Both single-stranded DNA and double-stranded DNA bind competitively to the same site. The intrinsic affinities are salt-dependent, and the formation of DNA-helicase complex is accompanied by a net release of 3-4 ions. Allosteric effects of nucleotide cofactors on RecQ binding to DNA were observed only for single-stranded DNA in the presence of 1.5 mM AMPPNP, whereas both AMPPNP and ADP had no detectable effect on double-stranded DNA binding over a large range of nucleotide cofactor concentrations.

International audience; The structure and thermal stability of a hetero chiral decaoligodeoxyribonucleotide duplex d(C1m8 G2C3G4C5LG6LC7G8C9G10)d(C11m8G12C13G14C15LG16LC17G18C19G20) (O1) with two contiguous pairs of enantiomeric 2'-deoxy-L-ribonucleotides (C5LG6L/C15LG16L) at its centre and an 8-methylguanine at position 2/12 was analysed by circular dichroism, NMR and molecular modelling. O1 resolves in a left-handed helical structure already at low salt concentration (0.1 M NaCl). The central L2-sugar portion assumes a B* left-handed conformation (mirror-image of right-handed B-DNA) while its flanking D4-sugar portions adopt the known Z left-handed conformation. The resulting Z4-B2*-Z4 structure (left-handed helix) is the reverse of that of B4-Z2*-B4 (right-handed helix) displayed by the nearly related decaoligodeoxyribonucleotide d(mC1G2mC3G4C5L G6LmC7G8mC9G10)2, at the same low salt concentration (0.1 M NaCl). In the same experimental conditions, d(C1m8G2C3G4C5G6C7G8C9G10)2 (O2), the stereoregular version of O1, resolves into a right-handed B-DNA helix. Thus, both the 8-methylguanine and the enantiomeric step CLpGL at the centre of the molecule are needed to induce left-handed helicity. Remarkably, in the various heterochiral decaoligodeoxyribonucleotides so far analysed by us, when the central CLpGL adopts the B* (respectively Z*) conformation, then the adjacent steps automatically resolves in the Z (respectively B) conformation. This allows a good optimisation of the base-base stackings and base-sugar van der Waals interactions at the ZB*/B*Z (respectively BZ*/Z*B) junctions so that the Z4-B2*-Z4 (respectively B4-Z2*-B4) helix displays a Tm (approximately 65 degrees C) that is only 5 degrees C lower than the one of its homochiral counterpart. Here we anticipate that a large variety of DNA helices can be generated at low salt concentration by manipulating internal factors such as sugar configuration, duplex length, nucleotide composition and base methylation. These helices can constitute powerful tools for structural and biological investigations, especially as they can be used in physiological conditions.

International audience; Although the immunoglobulin-like beta-sandwich fold has no specifically conserved function, some common structural features have been observed, in particular a structural motif, the tyrosine corner. Such a motif was described in neocarzinostatin (NCS), a bacterial protein the structure of which is very similar to that of the immunoglobulin domain. Compared with the other beta-sheet proteins, the NCS 'tyrosine corner' presents non-standard structural features. To investigate the role of this motif in the NCS structure and stability, we studied the properties of a mutant where the H bond interaction had been eliminated by replacing the tyrosine with a phenylalanine. This mutation costs 4.0 kcal/mol showing that the NCS 'tyrosine corner' is involved in protein stability as in the other Greek key proteins. This destabilization is accompanied by remote structural effects, including modification of the binding properties, suggesting an increase in the internal flexibility of the protein. With a view to using this protein for drug targeting, these results along with those obtained previously allow us to define clearly the limitations of the modifications that can be performed on this scaffold.

International audience; Mutagenized dockerin domains of endoglucanase CelD (type I) and of the cellulosome-integrating protein CipA (type II) were constructed by swapping residues 10 and 11 of the first or the second duplicated segment between the two polypeptides. These residues have been proposed to determine the specificity of cohesin-dockerin interactions. The dockerin domain of CelD still bound to the seventh cohesin domain of CipA (CohCip7), provided that mutagenesis occurred in one segment only. Binding was no longer detected by nondenaturing gel electrophoresis when both segments were mutagenized. The dockerin domain of CipA bound to the cohesin domain of SdbA as long as the second segment was intact. None of the mutated dockerins displayed detectable binding to the noncognate cohesin domain. Isothermal titration calorimetry showed that binding of the CelD dockerin to CohCip7 occurred with a high affinity [K(a) = (2.6 +/- 0.5) x 10(9) M(-1)] and a 1:1 stoichiometry. The reaction was weakly exothermic (DeltaHdegrees = -2.22 +/- 0.2 kcal x mol(-1)) and largely entropy driven (TDeltaSdegrees = 10.70 +/- 0.5 kcal x mol(-1)). The heat capacity change on complexation was negative (DeltaC(p) = -305 +/- 15 cal x mol(-1) x K(-1)). These values show that cohesin-dockerin binding is mainly hydrophobic. Mutations in the first or the second dockerin segment reduced or enhanced, respectively, the hydrophobic character of the interaction. Due to partial enthalpy-entropy compensation, these mutations induced only small changes in binding affinity. However, the binding affinity was strongly decreased when both segments were mutated, indicating strong negative cooperativity between the two mutated sites.

International audience; The only experimental strategy to address the structure of folding transition states, the so-called Φ value analysis, relies on the synergy between site directed mutagenesis and the measurement of reaction kinetics. Despite its importance, the Φ value analysis has been often criticized and its power to pinpoint structural information has been questioned. In this hypothesis, we demonstrate that comparing the Φ values between proteins not only allows highlighting the robustness of folding pathways but also provides per se a strong validation of the method.

International audience; Our ability to model the growth of microbes only relies on empirical laws, fundamentally restricting our understanding and predictive capacity in many environmental systems. In particular, the link between energy balances and growth dynamics is still not understood. Here we demonstrate a microbial growth equation relying on an explicit theoretical ground sustained by Boltzmann statistics, thus establishing a relationship between microbial growth rate and available energy. The validity of our equation was then questioned by analyzing the microbial isotopic fractionation phenomenon, which can be viewed as a kinetic consequence of the differences in energy contents of isotopic isomers used for growth. We illustrate how the associated theoretical predictions are actually consistent with recent experimental evidences. Our work links microbial population dynamics to the thermodynamic driving forces of the ecosystem, which opens the door to many biotechnological and ecological developments.

Oligostilbenoids (e.g., ampelopsin F, viniferin, pallidol) result from homogeneous or heterogeneous coupling of monomeric stilbenoid units, leading to various chemical structures. Oligostilbenoid synthesis is regio- and stereocontrolled. To tackle this regio- and stereocontrol, a supramolecular chemistry approach is required that can be achieved by quantum chemistry. The stability of noncovalent π-stacks, formed between two stilbenoid units prior to oxidation, is accurately evaluated with density functional theory (DFT) including dispersive effects (within the DFT-D formalism). These noncovalent arrangements drive the regiocontrol. The rest of the chemical pathway is a succession of dearomatization and rearomatization stages. The thermodynamics and kinetics of the processes are calculated with classical hybrid functionals. This study allows discrimination between the two main possible chemical pathways, namely, radical–neutral and radical–radical reactions. The former appears more likely, thermodynamics and kinetics being in perfect agreement with the experimental 1:2 ratio obtained for ampelopsin F:pallidol analogues, respectively. The authors thank the “Conseil Régional du Limousin” for financial support and CALI (CAlcul en LImousin) for computing facilities. Research in Limoges is also supported by COST actions (FA1003 “East-West Collaboration for Grapevine Diversity Exploration and Mobilization of Adaptive Traits for Breeding” and 0804 “Chemical Biology with Natural Compounds”). Financial support by the French Embassy in Malaysia and the MICINN of Spain (project CTQ2011-27253) is gratefully acknowledged. The authors gratefully acknowledge the support by the Operational Program Research and Development for Innovations−European Regional Development Fund (project CZ.1.05/2.1.00/03.0058 of the Ministry of Education, Youth and Sports of the Czech Republic).

The role of the denatured state in protein folding represents a key issue for the proper evaluation of folding kinetics and mechanisms. The yeast ortholog of the human frataxin, a mitochondrial protein essential for iron homeostasis and responsible for Friedreich's ataxia, has been shown to undergo cold denaturation above 0 °C, in the absence of chemical denaturants. This interesting property provides the unique opportunity to explore experimentally the molecular mechanism of both the hot and cold denaturation. In this work, we present the characterization of the temperature and urea dependence of the folding kinetics of yeast frataxin, and show that while at neutral pH and in the absence of a denaturant a simple two-state model may satisfactorily describe the temperature dependence of the folding and unfolding rate constants, the results obtained in urea over a wide range of pH reveal an intriguing complexity, suggesting that folding of frataxin involves a broad smooth free energy barrier. © 2014 the Partner Organisations.

International audience; A cationic amphiphilic peptide made of 10 leucine and 10 lysine residues, and four of its fluorescent derivatives in which leucines were substituted by Trp residues at different locations on the primary sequence have been synthesized. The interactions of these five peptides with neutral anionic or cationic vesicles were investigated using circular dichroism, steady state and time-resolved fluorescence with a combination of Trp quenching by brominated lipid probes, monolayers, modeling with minimization and simulated annealing procedures. We show that all the five peptides interact with neutral and anionic DMPC, DMPG, DOPC or egg yolk PC vesicles. The binding takes place whatever the peptide conformation in solution is. In the case of DMPC bilayers the binding free energy DeltaG is estimated at -8 kcal mole-1 and the number of phospholipid molecules involved is about 20-25 per peptide molecule. Peptides are bound as single-stranded alpha helices orientated parallel to the bilayer surface. In the anchoring of phospholipid head groups around the peptides, the lipid molecules are not smeared out in a plane parallel to the membrane surface but are organized around the hydrophilic face of the alpha helices like 'wheat grains around an ear' and protrude outside the bilayer towards the solvent. We suggest that such a lipid arrangement generates transient structural defects responsible for the membrane permeability enhancement. When an electrical potential is applied, the axis of the peptide helices remains parallel to the membrane surface and does not reorient to give rise to a bundle of helix monomers that forms transmembrane channels via a 'barrel stave' mechanism. The penetration depth of alpha helices in relation to the position of phosphorus atoms in the unperturbed lipid leaflet is estimated at 3.2 A.

International audience; To study conformational parameters of ligands before and after complex formation with the galactoside-binding agglutinin of Viscum album L. (VAA) in solution, combined computer-assisted random walk molecular mechanics (RAMM) calculations extended by conformational clustering analysis (CCA), molecular dynamics (MD) simulations as well as two-dimensional rotating-frame nuclear Overhauser effect (ROE) and two-dimensional nuclear Overhauser effect (NOE) spectroscopy NMR experiments were employed. Derivatives of the naturally occurring disaccharides Galbeta1-3GlcNAcbeta1-R and Galbeta1-3GalNAcbeta1-R as well as of a synthetic high-affinity binding partner, i.e. the disaccharide Galbeta1-2Galbeta1-R', were chosen as ligands in this study. The disaccharides displayed inherent flexibility in the valley of the global minimum between phi/psi combinations of (40 degrees/60 degrees) and (40 degrees/-60 degrees). Calculations of the de-N-acetylated sugars revealed that presence of this group did not markedly influence the distribution of low-energy conformers in the phi, psi, epsilon plot. Occupation of side minima at phi/psi (180 degrees/0 degrees) or (0 degrees/180 degrees) is either unlikely or low according to the results of MD simulations and RAMM calculations extended by CCA. Notably, these side minima define conformations which are not stable during a MD simulation. Transitions to other minima occur already a few picoseconds after the start of the simulation. NMR experiments of the free-state ligand confirmed the validity of the data sets obtained by the calculations. Following the description of the conformational space in the free-state NMR experiments were performed for these disaccharides complexed with VAA. They yielded two interresidual contacts for Galbeta1-3GlcNAcbeta1-R and Galbeta1-2Galbeta1-R'. The ligand conformations in the complex did not deviate markedly from those of a minimum conformation in the free state. One- and two-dimensional transferred nuclear Overhauser enhancement (TRNOE) experiments at different mixing times excluded the influence of spin-diffusion effects. When the NOE build-up curves in the three studied cases were compared, the residual mobility of the penultimate carbohydrate unit of Galbeta1-3GalNAcbeta1-R was observed to be higher than that of the respective hexopyranose unit of the other two bound ligands. Due to the availability of the conformational parameters of Galbeta1-2Galbeta1-R' in association with a galectin, namely the beta-galactoside-binding protein from chicken liver, it is remarkable to note that this ligand displays different conformations in the binding sites of either the plant or the animal lectin. They correspond to local energy-minimum conformations in the phi,psi, epsilon plot and substantiate differential conformer selection by these two lectins with identical nominal monosaccharide specificity.

International audience; The gene dnph1 (or rcl) encodes a hydrolase that cleaves the 2'-deoxyribonucleoside 5'-monophosphate (dNMP) N-glycosidic bond to yield a free nucleobase and 2-deoxyribose 5-phosphate. Recently, the crystal structure of rat DNPH1, a potential target for anti-cancer therapies, suggested that various analogs of AMP may inhibit this enzyme. From this result, we asked whether N (6)-substituted AMPs, and among them, cytotoxic cytokinin riboside 5'-monophosphates, may inhibit DNPH1. Here, we characterized the structural and thermodynamic aspects of the interactions of these various analogs with DNPH1. Our results indicate that DNPH1 is inhibited by cytotoxic cytokinins at concentrations that inhibit cell growth.

International audience

International audience; Various models on membrane structure and organization of proteins and complexes in natural membranes emerged during the last years. However, the lack of systematic dynamical studies to complement structural investigations hindered the establishment of a more complete picture of these systems. Elastic incoherent neutron scattering gives access to the dynamics on a molecular level and was applied to natural membranes extracted from the hyperthermophile Aquifex aeolicus and the mesophile Wolinella succinogenes bacteria. The results permitted to extract a hierarchy of dynamic flexibility and atomic resilience within the samples, which correlated with the organization of proteins in bioenergetics complexes and the functionality of the membranes.

International audience; Burkholderia ambifaria is a bacterium member of the Burkholderia cepacia complex (BCC), a closely related group of Gram-negative bacteria responsible for "cepacia syndrome" in immunocompromised patients. B. ambifaria produces BambL, a fucose-binding lectin that displays fine specificity to human fucosylated epitopes. Here, we report the first example of a synthetic ligand able to selectively bind, in the micromolar range, the pathogen-lectin BambL. The synthetic routes for the preparation of the α conformationally constrained fucoside are described, focusing on a totally diastereoselective inverse electron demand [4 + 2] Diels-Alder reaction. Isothermal titration calorimetry (ITC) demonstrated that this compound binds to the pathogen-associated lectin BambL with an affinity comparable to that of natural fucose-containing oligosaccharides. No binding was observed by LecB, a fucose-binding lectin from Pseudomonas aeruginosa, and the differences in affinity between the two lectins could be rationalized by modeling. Furthermore, SPR analyses showed that this fucomimetic does not bind to the human fucose-binding lectin DC-SIGN, thus supporting the selective binding profile towards B. ambifaria lectin.

International audience; Proteins are dynamic molecular machines whose structure and function are modulated by environmental perturbations and natural selection. Allosteric regulation, discovered in 1963 as a novel molecular mechanism of enzymatic adaptation [Monod, Changeux & Jacob (1963). J. Mol. Biol.6, 306-329], seems to be the leit motiv of enzymes and metabolic pathways, enabling fine and quick responses toward external perturbations. Hemoglobin (Hb), the oxygen transporter of all vertebrates, has been for decades the paradigmatic system to test the validity of the conformational selection mechanism, the conceptual innovation introduced by Monod, Wyman and Changeux. We present hereby the results of a comparative analysis of structure, function and thermodynamics of two extensively investigated hemoglobins: human HbA and trout HbI. They represent a unique and challenging comparison to test the general validity of the stereochemical model proposed by Perutz. Indeed both proteins are ideal for the purpose being very similar yet very different. In fact, T-HbI is a low-ligand-affinity cooperative tetrameric Hb, insensitive to all allosteric effectors. This remarkable feature, besides being physiologically sound, supports the stereochemical model, given that the six residues identified in HbA as responsible for the Bohr and the 2,3-di-phosphoglycerate effects are all mutated. Comparison of the three-dimensional structures of HbA and T-HbI allows unveiling the molecular mechanism whereby the latter has a lower O2 affinity. Moreover, the energetic balance sheet shows that the salt bridges breaking upon allosteric quaternary transition are important yet insufficient to account for the free energy of heme-heme interactions in both hemoglobins.

International audience; Neutron scattering measurements on protein powders reveal two deviations from harmonic dynamics at low temperature, whose molecular origin, physical nature and biological relevance are still matter of discussion. In this study we present a new experimental and theoretical approach to evidence the resolution dependence of anharmonic onsets: the use of strategically selected homomeric polypeptides allows revealing the exact resolution dependence; a two-site energy landscape model, where resolution effects are explicitly taken into account, is able to interpret quantitatively the experimental data in terms of energy landscape parameters. The energetic description provided by this analysis, together with recent experimental evidences obtained on chemically and structurally different peptide systems, allows us to connect the protein/water energy landscape structure with the two-wells water interaction potential proposed to explain the low-density→high-density liquid-liquid transition observed in supercooled water.

International audience; Incoherent neutron scattering is one of the most powerful tools for studying dynamics in biological matter. Using the cold neutron backscattering spectrometer IN16 at the Institut Laue Langevin (ILL, Grenoble, France), temperature dependence of cholinesterases' dynamics (human butyrylcholinesterase from plasma: hBChE; recombinant human acetylcholinesterase: hAChE and recombinant mouse acetylcholinesterase: mAChE) was examined using elastic incoherent neutron scattering (EINS). The dynamics was characterized by the averaged atomic mean square displacement (MSD), associated with the sample flexibility at a given temperature. We found MSD values of hAChE above the dynamical transition temperature (around 200K) larger than for mAChE and hBChE, implying that hAChE is more flexible than the other ChEs. Activation energies for thermodynamical transition were extracted through the frequency window model (FWM) (Becker et al. 2004) [1] and turned out to increase from hBChE to mAChE and finally to hAChE, inversely to the MSDs relations. Between 280 and 316K, catalytic studies of these enzymes were carried out using thiocholine esters: at the same temperature, the hAChE activity was systematically higher than the mAChE or hBChE ones. Our results thus suggest a strong correlation between dynamics and activity within the ChE family. We also studied and compared the ChEs thermal inactivation kinetics. Here, no direct correlation with the dynamics was observed, thus suggesting that relations between enzyme dynamics and catalytic stability are more complex. Finally, the possible relation between flexibility and protein ability to grow in crystals is discussed.

The precision with which X-ray crystallography and nuclear magnetic resonance (NMR) have provided structural models of biologically active and inactive conformations of countless proteins belies an easily overlooked dilemma. Proteins are inherently dynamic, exhibiting conformational freedom on timescales from picoseconds to seconds, implicating structural rearrangements that are essential for their biological function. Classical structural biology determines static models, that afford little insight into the underlying conformational equilibrium. The role that structural dynamics play in biological processes can only be understood by characterizing all thermally accessible protein conformations and their populations. NMR spectroscopy is uniquely sensitive to the presence of conformational dynamics in solution. Residual dipolar couplings (RDCs) measured in weakly aligned proteins, scalar couplings, and chemical shifts, probe motions occurring on timescales faster than 100 s of microseconds. These parameters therefore offer general tools to characterize protein motion on physiologically important timescales. A common approach to the dynamic interpretation of RDCs is to combine experimental restraint terms with a classical potential-energy force field to develop a conformational ensemble in agreement with experimental data. RDCs have also been exploited to characterize the conformational space sampled by the protein backbone either by fitting experimental data to determine angular excursions of internuclear bond vectors, or in comparison with different levels of accelerated molecular dynamics (AMD) to describe the most appropriate ensemble. Comparison of motions modeled using the Gaussian axial fluctuation (GAF) model, with ensembles derived from restraint-free AMD, demonstrated that such methods can provide a convergent description of protein motion.

International audience; In the mammalian mitochondrial translation apparatus, the proteins and their partner RNAs are coded by two genomes. The proteins are nuclear-encoded and resemble their homologs, whereas the RNAs coming from the rapidly evolving mitochondrial genome have lost critical structural information. This raises the question of molecular adaptation of these proteins to their peculiar partner RNAs. The crystal structure of the homodimeric bacterial-type human mitochondrial aspartyl-tRNA synthetase (DRS) confirmed a 3D architecture close to that of Escherichia coli DRS. However, the mitochondrial enzyme distinguishes by an enlarged catalytic groove, a more electropositive surface potential and an alternate interaction network at the subunits interface. It also presented a thermal stability reduced by as much as 12°C. Isothermal titration calorimetry analyses revealed that the affinity of the mitochondrial enzyme for cognate and non-cognate tRNAs is one order of magnitude higher, but with different enthalpy and entropy contributions. They further indicated that both enzymes bind an adenylate analog by a cooperative allosteric mechanism with different thermodynamic contributions. The larger flexibility of the mitochondrial synthetase with respect to the bacterial enzyme, in combination with a preserved architecture, may represent an evolutionary process, allowing nuclear-encoded proteins to cooperate with degenerated organelle RNAs.

International audience; Hybridization of nucleic acids on solid surfaces is a key process involved in high-throughput technologies such as microarrays and, in some cases, next-generation sequencing (NGS). A physical understanding of the hybridization process helps to determine the accuracy of these technologies. The goal of a widespread research program is to develop reliable transformations between the raw signals reported by the technologies and individual molecular concentrations from an ensemble of nucleic acids. This research has inputs from many areas, from bioinformatics and biostatistics, to theoretical and experimental biochemistry and biophysics, to computer simulations. A group of leading researchers met in Ploen Germany in 2011 to discuss present knowledge and limitations of our physico-chemical understanding of high-throughput nucleic acid technologies. This meeting inspired us to write this summary, which provides an overview of the state-of-the-art approaches based on physico-chemical foundation to modeling of the nucleic acids hybridization process on solid surfaces. In addition, practical application of current knowledge is emphasized.

International audience; The Flavivirus genus includes widespread and severe human pathogens like the four serotypes of dengue virus (DENV1 to DENV4), yellow fever virus, Japanese encephalitis virus and West Nile virus. Domain III (ED3) of the viral envelope protein interacts with cell receptors and contains epitopes recognized by virus neutralizing antibodies. Its structural, antigenic and immunogenic properties have been thoroughly studied contrary to its physico-chemical properties. Here, the ED3 domains of the above pathogenic flaviviruses were produced in the periplasm of Escherichia coli. Their thermodynamic stabilities were measured and compared in experiments of unfolding equilibriums, induced with chemicals or heat and monitored through protein fluorescence. A designed ED3 domain, with the consensus sequence of DENV strains from all serotypes, was highly stable. The low stability of the ED3 domain from DENV3 was increased by three changes of residues in the protein core without affecting its reactivity towards DENV-infected human serums. Additional changes showed that the stability of ED3 varied with the DENV3 genotype. The T(m) of ED3 was higher than 69°C for all the tested viruses and reached 86°C for the consensus ED3. The latter, deprived of its disulfide bond by mutations, was predominantly unfolded at 20°C. These results will help better understand and design the properties of ED3 for its use as diagnostic, vaccine or therapeutic tools.