Search

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 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; 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 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; 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.

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; 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; 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; Pseudomonas aeruginosa is an important opportunistic human pathogen that can establish bacterial communication by synchronizing the behavior of individual cells in a molecular phenomenon known as "quorum sensing". Through an elusive mechanism involving gene products of the pqs operon, the PqsE enzyme is absolutely required for the synthesis of extracellular phenazines, including the toxic blue pigment pyocyanin, effectively allowing cells to achieve full-fledged virulence. Despite several functional and structural attempts at deciphering the role of this relevant enzymatic drug target, no molecular function has yet been ascribed to PqsE. In the present study, we report a series of alanine scanning experiments aimed at altering the biological function of PqsE, allowing us to uncover key amino acid positions involved in the molecular function of this enzyme. We use sequence analysis and structural overlays with members of homologous folds to pinpoint critical positions located in the vicinity of the ligand binding cleft and surrounding environment, revealing the importance of a unique C-terminal α-helical motif in the molecular function of PqsE. Our results suggest that the active site of the enzyme involves residues that extend further into the hydrophobic core of the protein, advocating for a lid-like movement of the two terminal helices. This information should help design virtual libraries of PqsE inhibitors, providing means to counter P. aeruginosa virulence acquisition and helping to reduce nosocomial infections.

International audience; Herein, we discuss how, why, and when cascade complexation reactions produce stable, mononuclear, luminescent ternary complexes, by considering the binding of hexafluoroacetylacetonate anions (hfac(-)) and neutral, semi-rigid, tridentate 2,6-bis(benzimidazol-2-yl)pyridine ligands (Lk) to trivalent lanthanide atoms (Ln(III)). The solid-state structures of [Ln(Lk)(hfac)(3)] (Ln=La, Eu, Lu) showed that [Ln(hfac)(3)] behaved as a neutral six-coordinate lanthanide carrier with remarkable properties: 1) the strong cohesion between the trivalent cation and the didentate hfac anions prevented salt dissociation; 2) the electron-withdrawing trifluoromethyl substituents limited charge-neutralization and favored cascade complexation with Lk; 3) nine-coordination was preserved for [Ln(Lk)(hfac)(3)] for the complete lanthanide series, whilst a counterintuitive trend showed that the complexes formed with the smaller lanthanide elements were destabilized. Thermodynamic and NMR spectroscopic studies in solution confirmed that these characteristics were retained for solvated molecules, but the operation of concerted anion/ligand transfers with the larger cations induced subtle structural variations. Combined with the strong red photoluminescence of [Eu(Lk)(hfac)(3)], the ternary system Ln(III)/hfac(-)/Lk is a promising candidate for the planned metal-loading of preformed multi-tridentate polymers.

International audience; In the objective of developing ligands that simultaneously satisfy the requirements for MRI contrast agents and near-infrared emitting optical probes that are suitable for imaging, three isoquinoline-based polyaminocarboxylate ligands, L1, L2 and L3, have been synthesized and the corresponding Gd(3+), Nd(3+) and Yb(3+) complexes investigated. The specific challenge of the present work was to create NIR emitting agents which (i) have excitation wavelengths compatible with biological applications and (ii) are able to emit a sufficient number of photons to ensure sensitive NIR detection for microscopic imaging. Here we report the first observation of a NIR signal arising from a Ln(3+) complex in aqueous solution in a microscopy setup. The lanthanide complexes have high thermodynamic stability (log K(LnL) =17.7-18.7) and good selectivity for lanthanide ions versus the endogenous cations Zn(2+), Cu(2+), and Ca(2+) thus preventing transmetalation. A variable temperature and pressure (17)O NMR study combined with nuclear magnetic relaxation dispersion measurements yielded the microscopic parameters characterizing water exchange and rotation. Bishydration of the lanthanide cation in the complexes, an important advantage to obtain high relaxivity for the Gd(3+) chelates, has been demonstrated by (17)O chemical shifts for the Gd(3+) complexes and by luminescence lifetime measurements for the Yb(3+) analogues. The water exchange on the three Gd(3+) complexes is considerably faster (k(ex)(298) = (13.9-15.4) × 10(6) s(-1)) than on commercial Gd(3+)-based contrast agents and proceeds via a dissociative mechanism, as evidenced by the large positive activation volumes for GdL1 and GdL2 (+10.3 ± 0.9 and +10.6 ± 0.9 cm(3) mol(-1), respectively). The relaxivity of GdL1 is doubled at 40 MHz and 298 K in fetal bovine serum (r(1) = 16.1 vs 8.5 mM(-1) s(-1) in HEPES buffer), due to hydrophobic interactions between the chelate and serum proteins. The isoquinoline core allows for the optimization of the optical properties of the luminescent lanthanide complexes in comparison to the pyridinic analogues and provides significant shifts of the excitation energies toward lower values which therefore become more adapted for biological applications. L2 and L3 bear two methoxy substituents on the aromatic core in ortho and para positions, respectively, that further modulate their electronic structure. The Nd(3+) and Yb(3+) complexes of the ligand L3, which incorporates the p-dimethoxyisoquinoline moiety, can be excited up to 420 nm. This wavelength is shifted over 100 nm toward lower energy in comparison to the pyridine-based analogue. The luminescence quantum yields of the Nd(3+) (0.013-0.016%) and Yb(3+) chelates (0.028-0.040%) are in the range of the best nonhydrated complexes, despite the presence of two inner sphere water molecules. More importantly, the 980 nm NIR emission band of YbL3 was detected with a good sensitivity in a proof of concept microscopy experiment at a concentration of 10 μM in fetal bovine serum. Our results demonstrate that even bishydrated NIR lanthanide complexes can emit a sufficient number of photons to ensure sensitive detection in practical applications. In particular, these ligands containing an aromatic core with coordinating pyridine nitrogen can be easily modified to tune the optical properties of the NIR luminescent lanthanide complexes while retaining good complex stability and MRI characteristics for the Gd(3+) analogues. They constitute a highly versatile platform for the development of bimodal MR and optical imaging probes based on a simple mixture of Gd(3+) and Yb(3+)/Nd(3+) complexes using an identical chelator. Given the presence of two inner sphere water molecules, important for MRI applications of the corresponding Gd(3+) analogues, this result is particularly exciting and opens wide perspectives not only for NIR imaging based on Ln(3+) ions but also for the design of combined NIR optical and MRI probes.

The protein ERp57 (also known as PDIA3) is a widely distributed protein, mainly localized in the endoplasmic reticulum, where it acts as disulfide isomerase, oxidoreductase and chaperone, in concert with the lectins calreticulin (CRT) and calnexin. The ERp57/CRT complex has been detected on the cell surface and previous studies have suggested its involvement in programmed cell death. Although the ERp57-CRT complex has been characterized, little is known about its role in different cellular compartments as well as inhibitors of this interaction. We focused on the kinetic, extent and stability of the ERp57-CRT complex, using the surface plasmon resonance spectroscopy, investigating the possible role as inhibitor of the antibiotic vancomycin. Equilibrium thermodynamic data suggested that vancomycin may hinder the interaction between the two proteins and could interfere with the ERp57 conformational changes that stabilize the complex. Furthermore, by means of confocal microscopy, we evaluated the effect of the in vivo administration of vancomycin on the ERp57/CRT complex on the surface of HeLa cells. The model presented here could be used for the search of other specific inhibitors/interactors of ERp57, which can be extremely helpful to understand the biological pathways where the protein is involved and to modulate its activity.

International audience; Current knowledge on the reaction whereby a protein acquires its native three-dimensional structure was obtained by and large through characterization of the folding mechanism of simple systems. Given the multiplicity of amino acid sequences and unique folds, it is not so easy, however, to draw general rules by comparing folding pathways of different proteins. In fact, quantitative comparison may be jeopardized not only because of the vast repertoire of sequences but also in view of a multiplicity of structures of the native and denatured states. We have tackled the problem of the relationships between the sequence information and the folding pathway of a protein, using a combination of kinetics, protein engineering and computational methods, applied to relatively simple systems. Our strategy has been to investigate the folding mechanism determinants using two complementary approaches, i.e. (i) the study of members of the same family characterized by a common fold, but substantial differences in amino acid sequence, or (ii) heteromorphic pairs characterized by largely identical sequences but with different folds. We discuss some recent data on protein-folding mechanisms by presenting experiments on different members of the PDZ domain family and their circularly permuted variants. Characterization of the energetics and structures of intermediates and TSs (transition states), obtained by Φ-value analysis and restrained MD (molecular dynamics) simulations, provides a glimpse of the malleability of the dynamic states and of the role of the topology of the native states and of the denatured states in dictating folding and misfolding pathways.

Protein-protein interactions mediated by modular protein domains are critical for cell scaffolding, differentiation, signaling, and ultimately, evolution. Given the vast number of ligands competing for binding to a limited number of domain families, it is often puzzling how specificity can be achieved. Selectivity may be modulated by intradomain allostery, whereby a remote residue is energetically connected to the functional binding site via side chain or backbone interactions. Whereas several energetic pathways, which could mediate intradomain allostery, have been predicted in modular protein domains, there is a paucity of experimental data to validate their existence and roles. Here, we have identified such functional energetic networks in one of the most common protein-protein interaction modules, the PDZ domain. We used double mutant cycles involving site-directed mutagenesis of both the PDZ domain and the peptide ligand, in conjunction with kinetics to capture the fine energetic details of the networks involved in peptide recognition. We performed the analysis on two homologous PDZ-ligand complexes and found that the energetically coupled residues differ for these two complexes. This result demonstrates that amino acid sequence rather than topology dictates the allosteric pathways. Furthermore, our data support a mechanism whereby the whole domain and not only the binding pocket is optimized for a specific ligand. Such cross-talk between binding sites and remote residues may be used to fine tune target selectivity.

International audience; Understanding the driving forces governing protein assembly requires the characterization of interactions at molecular level. We focus on two homologous oppositely charged proteins, lysozyme and α-lactalbumin, which can assemble into microspheres. The assembly early steps were characterized through the identification of interacting surfaces monitored at residue level by NMR chemical shift perturbations by titrating one 15N-labeled protein with its unlabeled partner. While α-lactalbumin has a narrow interacting site, lysozyme has interacting sites scattered on a broad surface. The further assembly of these rather unspecific heterodimers into tetramers leads to the establishment of well-defined interaction sites. Within the tetramers, most of the electrostatic charge patches on the protein surfaces are shielded. Then, hydrophobic interactions, which are possible because α-lactalbumin is in a partially folded state, become preponderant, leading to the formation of larger oligomers. This approach will be particularly useful for rationalizing the design of protein assemblies as nanoscale devices.View: ACS ActiveView PDF | PDF | PDF w/ Links | Full Text HTMLCiting Articles

International audience; Archaeal initiation factor 2 (aIF2) is a protein involved in the initiation of protein biosynthesis. In its GTP-bound, "ON" conformation, aIF2 binds an initiator tRNA and carries it to the ribosome. In its GDP-bound, "OFF" conformation, it dissociates from tRNA. To understand the specific binding of GTP and GDP and its dependence on the ON or OFF conformational state of aIF2, molecular dynamics free energy simulations (MDFE) are a tool of choice. However, the validity of the computed free energies depends on the simulation model, including the force field and the boundary conditions, and on the extent of conformational sampling in the simulations. aIF2 and other GTPases present specific difficulties; in particular, the nucleotide ligand coordinates a divalent Mg(2+) ion, which can polarize the electronic distribution of its environment. Thus, a force field with an explicit treatment of electronic polarizability could be necessary, rather than a simpler, fixed charge force field. Here, we begin by comparing a fixed charge force field to quantum chemical calculations and experiment for Mg(2+):phosphate binding in solution, with the force field giving large errors. Next, we consider GTP and GDP bound to aIF2 and we compare two fixed charge force fields to the recent, polarizable, AMOEBA force field, extended here in a simple, approximate manner to include GTP. We focus on a quantity that approximates the free energy to change GTP into GDP. Despite the errors seen for Mg(2+):phosphate binding in solution, we observe a substantial cancellation of errors when we compare the free energy change in the protein to that in solution, or when we compare the protein ON and OFF states. Finally, we have used the fixed charge force field to perform MDFE simulations and alchemically transform GTP into GDP in the protein and in solution. With a total of about 200 ns of molecular dynamics, we obtain good convergence and a reasonable statistical uncertainty, comparable to the force field uncertainty, and somewhat lower than the predicted GTP/GDP binding free energy differences. The sign and magnitudes of the differences can thus be interpreted at a semiquantitative level, and are found to be consistent with the experimental binding preferences of ON- and OFF-aIF2.

We report the results of molecular dynamics simulations of an off-lattice protein model featuring a physical force-field and amino-acid sequence. We show that localized modes of nonlinear origin (discrete breathers) emerge naturally as continuations of a subset of high-frequency normal modes residing at specific sites dictated by the native fold. In the case of the small $\beta$-barrel structure that we consider, localization occurs on the turns connecting the strands. At high energies, discrete breathers stabilize the structure by concentrating energy on few sites, while their collapse marks the onset of large-amplitude fluctuations of the protein. Furthermore, we show how breathers develop as energy-accumulating centres following perturbations even at distant locations, thus mediating efficient and irreversible energy transfers. Remarkably, due to the presence of angular potentials, the breather induces a local static distortion of the native fold. Altogether, the combination of this two nonlinear effects may provide a ready means for remotely controlling local conformational changes in proteins., Comment: Submitted to Physical Biology

Background Disrupting protein-protein interactions by small organic molecules is nowadays a promising strategy employed to block protein targets involved in different pathologies. However, structural changes occurring at the binding interfaces make difficult drug discovery processes using structure-based drug design/virtual screening approaches. Here we focused on two homologous calcium binding proteins, calmodulin and human centrin 2, involved in different cellular functions via protein-protein interactions, and known to undergo important conformational changes upon ligand binding. Results In order to find suitable protein conformations of calmodulin and centrin for further structure-based drug design/virtual screening, we performed in silico structural/energetic analysis and molecular docking of terphenyl (a mimicking alpha-helical molecule known to inhibit protein-protein interactions of calmodulin) into X-ray and NMR ensembles of calmodulin and centrin. We employed several scoring methods in order to find the best protein conformations. Our results show that docking on NMR structures of calmodulin and centrin can be very helpful to take into account conformational changes occurring at protein-protein interfaces. Conclusions NMR structures of protein-protein complexes nowadays available could efficiently be exploited for further structure-based drug design/virtual screening processes employed to design small molecule inhibitors of protein-protein interactions.

International audience; Actin filaments are semiflexible polymers that display large-scale conformational twisting and bending motions. Modulation of filament bending and twisting dynamics has been linked to regulatory actin-binding protein function, filament assembly and fragmentation, and overall cell motility. The relationship between actin filament bending and twisting dynamics has not been evaluated. The numerical and analytical experiments presented here reveal that actin filaments have a strong intrinsic twist-bend coupling that obligates the reciprocal interconversion of bending energy and twisting stress. We developed a mesoscopic model of actin filaments that captures key documented features, including the subunit dimensions, interaction energies, helicity, and geometrical constraints coming from the double-stranded structure. The filament bending and torsional rigidities predicted by the model are comparable to experimental values, demonstrating the capacity of the model to assess the mechanical properties of actin filaments, including the coupling between twisting and bending motions. The predicted actin filament twist-bend coupling is strong, with a persistence length of 0.15-0.4 μm depending on the actin-bound nucleotide. Twist-bend coupling is an emergent property that introduces local asymmetry to actin filaments and contributes to their overall elasticity. Up to 60% of the filament subunit elastic free energy originates from twist-bend coupling, with the largest contributions resulting under relatively small deformations. A comparison of filaments with different architectures indicates that twist-bend coupling in actin filaments originates from their double protofilament and helical structure.

The acid/base properties of proteins are essential in biochemistry, and proton binding is a valuable reporter on electrostatic interactions. We propose a constant-pH Monte Carlo strategy to compute protonation free energies and pKas. The solvent is described implicitly, through a generalized Born model. The electronic polarizability and backbone motions of the protein are included through the protein dielectric constant. Side chain motions are described explicitly, by the Monte Carlo scheme. An efficient computational algorithm is described, which allows us to treat the fluctuating shape of the protein/solvent boundary in a way that is numerically exact (within the GB framework) this contrasts with several previous constant-pH approaches. For a test set of six proteins and 78 titratable groups, the model performs well, with an rms error of 1.2 pH units. While this is slightly greater than a simple Null model (rms error of 1.1) and a fully empirical model (rms error of 0.9), it is obtained using physically meaningful model parameters, including a low protein dielectric of four. Importantly, similar performance is obtained for side chains with large and small pKa shifts (relative to a standard model compound). The titration curve slopes and the conformations sampled are reasonable. Several directions to improve the method further are discussed. © 2010 American Chemical Society. 114 32 10634 10648 Cited By :20

International audience; Imatinib is an important anticancer drug, which binds specifically to the Abl kinase and blocks its signalling activity. To model imatinib:protein interactions, we have developed a molecular mechanics force field for imatinib and four close analogues, which is consistent with the CHARMM force field for proteins and nucleic acids. 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 imatinib fragments at 32 different positions. We considered both a neutral and a protonated imatinib. The final RMS deviation between the ab initio and force field energies, averaged over both forms, was 0.2 kcal/mol. The model also reproduces the ab initio geometry and flexibility of imatinib. To apply the force field to imatinib:Abl simulations, it is also necessary to determine the most likely imatinib protonation state when it binds to Abl. This was done using molecular dynamics free energy simulations, where imatinib is reversibly protonated during a series of MD simulations, both in solution and in complex with Abl. The simulations indicate that imatinib binds to Abl in its protonated, positively-charged form. To help test the force field and the protonation prediction, we did MD free energy simulations that compare the Abl binding affinities of two imatinib analogs, obtaining good agreement with experiment. Finally, two new imatinib variants were considered, one of which is predicted to have improved Abl binding. This variant could be of interest as a potential drug.

International audience; Many enzymes possess, besides their native function, additional promiscuous activities. Proteins with several activities (multipurpose catalysts) may have a wide range of biotechnological and biomedical applications. Natural promiscuity, however, appears to be of limited scope in this context, because the latent (promiscuous) function is often related to the evolved one (sharing the active site and even the chemical mechanism) and its enhancement upon suitable mutations usually brings about a decrease in the native activity. Here we explore the use of computational protein design to overcome these limitations. The high-plasticity positions close to the original ("native") active-site are the most promising candidates for mutations that create a second active-site associated to a new function. To avoid compromising protein folding and native activity, we propose a minimal-perturbation approach based on the combinatorial optimization of, both the de novo catalytic activity and the folding free-energy: essentially, we construct the Pareto Set of optimal stability/promiscuous-function solutions. We validate our approach by introducing a promiscuous esterase activity in E. coli thioredoxin on the basis of mutations at positions close to the native-active-site disulfide-bridge. Native oxidoreductase activity is not compromised and it is, in fact, found to be 1.5-fold enhanced, as determined by an insulin-reduction assay. This work provides general guidelines as to how computational design can be used to expand the scope and applications of protein promiscuity. From a more general viewpoint, it illustrates the potential of multi-objective optimization as the computational analogue of multi-feature natural selection.

The proteome of the radiation- and desiccation-resistant bacterium D. radiodurans features a group of proteins that contain significant intrinsically disordered regions that are not present in non-extremophile homologues. Interestingly, this group includes a number of housekeeping and repair proteins such as DNA polymerase III, nudix hydrolase and rotamase. Here, we focus on a member of the nudix hydrolase family from D. radiodurans possessing low-complexity N- and C-terminal tails, which exhibit sequence signatures of intrinsic disorder and have unknown function. The enzyme catalyzes the hydrolysis of oxidatively damaged and mutagenic nucleotides, and it is thought to play an important role in D. radiodurans during the recovery phase after exposure to ionizing radiation or desiccation. We use molecular dynamics simulations to study the dynamics of the protein, and study its hydration free energy using the GB/SA formalism. We show that the presence of disordered tails significantly decreases the hydration free energy of the whole protein. We hypothesize that the tails increase the chances of the protein to be located in the remaining water patches in the desiccated cell, where it is protected from the desiccation effects and can function normally. We extrapolate this to other intrinsically disordered regions in proteins, and propose a novel function for them: intrinsically disordered regions increase the “surface-properties” of the folded domains they are attached to, making them on the whole more hydrophilic and potentially influencing, in this way, their localization and cellular activity., PLoS Computational Biology, 6 (7), ISSN:1553-734X, ISSN:1553-7358

International audience; Free energy simulations compare multiple ligand:receptor complexes by "alchemically" transforming one into another, yielding binding free energy differences. Since their introduction in the 1980s, many technical and theoretical obstacles were surmounted, and the method ("MDFE," since molecular dynamics are often used) has matured into a powerful tool. We describe its current status, its effectiveness, and the challenges it faces. MDFE has provided chemical accuracy for many systems but remains expensive, with significant human overhead costs. The bottlenecks have shifted, partly due to increased computer power. To study diverse sets of ligands, force field availability and accuracy can be a major difficulty. Another difficulty is the frequent need to consider multiple states, related to sidechain protonation or buried waters, for example. Sophisticated, automated methods to sample these states are maturing, such as constant pH simulations. Meanwhile, combinations of MDFE and simpler approaches, like continuum dielectric models, can be very effective. As illustrations, we show how, with careful force field parameterization, MDFE accurately predicts binding specificities between complex tetracycline ligands and their targets. We describe substrate binding to the aspartyl-tRNA synthetase enzyme, where many distinct electrostatic states play a role, and a histidine and a Mg(2+) ion act as coupled switches that help enforce a strict preference for the aspartate substrate, relative to several analogs. Overall, MDFE has achieved a predictive status, where novel ligands can be studied and molecular recognition elucidated in depth. It should play an increasing role in the analysis of complex cellular processes and biomolecular engineering.

International audience; A method for computational design of protein-ligand interactions is implemented and tested on the asparaginyl- and aspartyl-tRNA synthetase enzymes (AsnRS, AspRS). The substrate specificity of these enzymes is crucial for the accurate translation of the genetic code. The method relies on a molecular mechanics energy function and a simple, continuum electrostatic, implicit solvent model. As test calculations, we first compute AspRS-substrate binding free energy changes due to nine point mutations, for which experimental data are available; we also perform large-scale redesign of the entire active site of each enzyme (40 amino acids) and compare to experimental sequences. We then apply the method to engineer an increased binding of aspartyl-adenylate (AspAMP) into AsnRS. Mutants are obtained using several directed evolution protocols, where four or five amino acid positions in the active site are randomized. Promising mutants are subjected to molecular dynamics simulations; Poisson-Boltzmann calculations provide an estimate of the corresponding, AspAMP, binding free energy changes, relative to the native AsnRS. Several of the mutants are predicted to have an inverted binding specificity, preferring to bind AspAMP rather than the natural substrate, AsnAMP. The computed binding affinities are significantly weaker than the native, AsnRS:AsnAMP affinity, and in most cases, the active site structure is significantly changed, compared to the native complex. This almost certainly precludes catalytic activity. One of the designed sequences has a higher affinity and more native-like structure and may represent a valid candidate for Asp activity.

International audience; Phosphatidyl-myo-inositol mannosyltransferase A (PimA) is an essential glycosyltransferase (GT) involved in the biosynthesis of phosphatidyl-myo-inositol mannosides (PIMs), which are key components of the mycobacterial cell envelope. PimA is the paradigm of a large family of peripheral membrane-binding GTs for which the molecular mechanism of substrate/membrane recognition and catalysis is still unknown. Strong evidence is provided showing that PimA undergoes significant conformational changes upon substrate binding. Specifically, the binding of the donor GDP-Man triggered an important interdomain rearrangement that stabilized the enzyme and generated the binding site for the acceptor substrate, phosphatidyl-myo-inositol (PI). The interaction of PimA with the beta-phosphate of GDP-Man was essential for this conformational change to occur. In contrast, binding of PI had the opposite effect, inducing the formation of a more relaxed complex with PimA. Interestingly, GDP-Man stabilized and PI destabilized PimA by a similar enthalpic amount, suggesting that they formed or disrupted an equivalent number of interactions within the PimA complexes. Furthermore, molecular docking and site-directed mutagenesis experiments provided novel insights into the architecture of the myo-inositol 1-phosphate binding site and the involvement of an essential amphiphatic alpha-helix in membrane binding. Altogether, our experimental data support a model wherein the flexibility and conformational transitions confer the adaptability of PimA to the donor and acceptor substrates, which seems to be of importance during catalysis. The proposed mechanism has implications for the comprehension of the peripheral membrane-binding GTs at the molecular level.

International audience; A new continuum model is presented for computing the solvation free energies of cations in water. It combines in a single formalism based on statistical thermodynamics the Poisson model for electrostatics with the Langevin dipole model to account for nonuniform water dipole distribution around the ions. An excellent match between experimental and computed solvation free energies is obtained for 10 monovalent and divalent ions.

International audience; Protein design has many applications not only in biotechnology but also in basic science. It uses our current knowledge in structural biology to predict, by computer simulations, an amino acid sequence that would produce a protein with targeted properties. As in other examples of synthetic biology, this approach allows the testing of many hypotheses in biology. The recent development of automated computational methods to design proteins has enabled proteins to be designed that are very different from any known ones. Moreover, some of those methods mostly rely on a physical description of atomic interactions, which allows the designed sequences not to be biased towards known proteins. In this paper, we will describe the use of energy functions in computational protein design, the use of atomic models to evaluate the free energy in the unfolded and folded states, the exploration and optimization of amino acid sequences, the problem of negative design and the design of biomolecular function. We will also consider its use together with the experimental techniques such as directed evolution. We will end by discussing the challenges ahead in computational protein design and some of their future applications.

International audience; A wealth of information on RNA folding and ribonucleoprotein assembly has emerged from analyses of structures and from the use of innovative biophysical tools. Although integrating data obtained from static structures with dynamic measurements presents major challenges, such efforts are opening new vistas on the RNA folding landscape.

International audience; Tetracycline (Tc) is a broad-spectrum antibiotic that kills bacteria by interrupting protein biosynthesis. It is thought that the bacteriostatic action of Tc is associated with its binding to the acceptor site (or A site) in the bacterial ribosome, interfering with the attachment of aminoacyl-tRNA. Recently, however, the crystal structure of a complex between Tc and trypsin-modified elongation factor Tu (tm-EF-Tu) was determined, raising the question of whether Tc binding to EF-Tu has a role in its inhibition of protein synthesis. We address this question using computer simulations. As controls, we first compute relative ribosome binding free energies for seven Tc variants for which experimental data are available, obtaining good agreement. We then consider the binding of Tc to both the trypsin-modified and unmodified EF-Tu-GDP complexes. We show that the direct contribution of EF-Tu to the binding free energy is negligible; rather, the binding can be solely attributed to interactions of Tc with a bridging Mg(2+) ion and the GDP phosphate groups. The effects of trypsin modification are modest. Further, our calculations show that EF-Tu does not exhibit any binding preference for Tc over the nonantibiotic, 4-dedimethyl-Tc, and EF-Tu does not bind the Tc analogue tigecycline, which is a potent antibiotic. In contrast, both the ribosome and the Tet Repressor protein (involved in Tc resistance) do show a binding preference for Tc over 4-dedimethyl-Tc, and the ribosome prefers to bind tigecycline over Tc. Overall, our results provide insights into the binding properties of tetracyclines and support the idea that EF-Tu is not their primary target.

International audience; TonB is a cytoplasmic membrane protein required for active transport of various essential substrates such as heme and iron siderophores through the outer membrane receptors of Gram-negative bacteria. This protein spans the periplasm, contacts outer membrane transporters by its C-terminal domain, and transduces energy from the protonmotive force to the transporters. The TonB box, a relatively conserved sequence localized on the periplasmic side of the transporters, has been shown to directly contact TonB. While Serratia marcescens TonB functions with various transporters, HasB, a TonB-like protein, is dedicated to the HasR transporter. HasR acquires heme either freely or via an extracellular heme carrier, the hemophore HasA, that binds to HasR and delivers heme to the transporter. Here, we study the interaction of HasR with a HasB C-terminal domain and compare it with that obtained with a TonB C-terminal fragment. Analysis of the thermodynamic parameters reveals that the interaction mode of HasR with HasB differs from that with TonB, the difference explaining the functional specificity of HasB for HasR. We also demonstrate that the presence of the substrate on the extracellular face of the transporter modifies, via enthalpy-entropy compensation, the interaction with HasB on the periplasmic face. The transmitted signal depends on the nature of the substrate. While the presence of heme on the transporter modifies only slightly the nature of interactions involved between HasR and HasB, hemophore binding on the transporter dramatically changes the interactions and seems to locally stabilize some structural motifs. In both cases, the HasR TonB box is the target for those modifications.

International audience; Specific recognitions of GNRA tetraloops by small helical receptors are among the most widespread long-range packing interactions in large ribozymes. However, in contrast to GYRA and GAAA tetraloops, very few GNRA/receptor interactions have yet been identified to involve GGAA tetraloops in nature. A novel in vitro selection scheme based on a rigid self-assembling tectoRNA scaffold designed for isolation of intermolecular interactions with A-minor motifs has yielded new GGAA tetraloop-binding receptors with affinity in the nanomolar range. One of the selected receptors is a novel 12 nt RNA motif, (CCUGUG ... AUCUGG), that recognizes GGAA tetraloop hairpin with a remarkable specificity and affinity. Its physical and chemical characteristics are comparable to those of the well-studied '11nt' GAAA tetraloop receptor motif. A second less specific motif (CCCAGCCC ... GAUAGGG) binds GGRA tetraloops and appears to be related to group IC3 tetraloop receptors. Mutational, thermodynamic and comparative structural analysis suggests that natural and in vitro selected GNRA receptors can essentially be grouped in two major classes of GNRA binders. New insights about the evolution, recognition and structural modularity of GNRA and A-minor RNA-RNA interactions are proposed.

International audience; Tetracyclines (Tc) are important antibiotics that inhibit bacterial ribosomes. Two and six Tc binding sites, respectively, were seen in two X-ray structures of the Thermus thermophilus 30S ribosome subunit; the exact functional role of the various sites remains unclear. We study the two consensus sites, seen in both structures: a primary site, which is positioned to block tRNA A-site binding, and a secondary site, which has a weaker electron density. We combine molecular dynamics simulations and continuum electrostatic calculations to estimate the relative affinities of the two sites. The dielectric constant of the ribosome is set to 8, to reproduce the experimental binding free energy differences between Tc and its analogues minocycline and doxycycline, as well as more rigorous free energy simulations of Mg2+ binding. We find that both sites include a prebound Mg2+ ion, present before Tc binds. Using long simulations and comparing eight structural models for each site, we then show that primary site Tc binding is stronger by 1−4 kcal/mol; this range appears consistent with the crystallographically observed occupancies of the two sites. With this free energy range, TET5 is largely unoccupied under physiological conditions. Thus, we propose that the primary site is the inhibitory site and that allosteric effects may not be essential for tetracycline function. Copyright (2008) American Chemical Society.

International audience; Tetracyclines (Tc) are important antibiotics that inhibit bacterial ribosomes. Two and six Tc binding sites, respectively, were seen in two X-ray structures of the Thermus thermophilus 30S ribosome subunit; the exact functional role of the various sites remains unclear. We study the two consensus sites, seen in both structures: a primary site, which is positioned to block tRNA A-site binding, and a secondary site, which has a weaker electron density. We combine molecular dynamics simulations and continuum electrostatic calculations to estimate the relative affinities of the two sites. The dielectric constant of the ribosome is set to 8, to reproduce the experimental binding free energy differences between Tc and its analogues minocycline and doxycycline, as well as more rigorous free energy simulations of Mg2+ binding. We find that both sites include a prebound Mg2+ ion, present before Tc binds. Using long simulations and comparing eight structural models for each site, we then show that primary site Tc binding is stronger by 1−4 kcal/mol; this range appears consistent with the crystallographically observed occupancies of the two sites. With this free energy range, TET5 is largely unoccupied under physiological conditions. Thus, we propose that the primary site is the inhibitory site and that allosteric effects may not be essential for tetracycline function. Copyright (2008) American Chemical Society.

Epigénome er paléogénomie; International audience; Nucleosome positioning plays an essential role in cellular processes by modulating accessibility of DNA to proteins. Here, using only sequence-dependent DNA flexibility and intrinsic curvature, we predict the nucleosome occupancy along the genomes of Saccharomyces cerevisiae and Drosophila melanogaster and demonstrate the predictive power and universality of our model through its correlation with experimentally determined nucleosome occupancy data. In yeast promoter regions, the computed average nucleosome occupancy closely superimposes with experimental data, exhibiting a

International audience; Owing to a striking, and most likely fortuitous, structural and sequence similarity with the bacterial 16 S ribosomal A site, the RNA kissing-loop complex formed by the HIV-1 genomic RNA dimerization initiation site (DIS) specifically binds 4,5-disubstituted 2-deoxystreptamine (2-DOS) aminoglycoside antibiotics. We used chemical probing, molecular modeling, isothermal titration calorimetry (ITC) and UV melting to investigate aminoglycoside binding to the DIS loop-loop complex. We showed that apramycin, an aminoglycoside containing a bicyclic moiety, also binds the DIS, but in a different way than 4,5-disubstituted 2-DOS aminoglycosides. The determination of thermodynamic parameters for various aminoglycosides revealed the role of the different rings in the drug-RNA interaction. Surprisingly, we found that the affinity of lividomycin and neomycin for the DIS (K(d) approximately 30 nM) is significantly higher than that obtained in the same experimental conditions for their natural target, the bacterial A site (K(d) approximately 1.6 microM). In good agreement with their respective affinity, aminoglycoside increase the melting temperature of the loop-loop interaction and also block the conversion from kissing-loop complex to extended duplex. Taken together, our data might be useful for selecting new molecules with improved specificity and affinity toward the HIV-1 DIS RNA.

International audience; The non-linear problem of simulating the structural transition between two known forms of a macromolecule still remains a challenge in structural biology. The problem is usually addressed in an approximate way using 'morphing' techniques, which are linear interpolations of either the Cartesian or the internal coordinates between the initial and end states, followed by energy minimization. Here we describe a web tool that implements a new method to calculate the most probable trajectory that is exact for harmonic potentials; as an illustration of the method, the classical Calpha-based Elastic Network Model (ENM) is used both for the initial and the final states but other variants of the ENM are also possible. The Langevin equation under this potential is solved analytically using the Onsager and Machlup action minimization formalism on each side of the transition, thus replacing the original non-linear problem by a pair of linear differential equations joined by a non-linear boundary matching condition. The crossover between the two multidimensional energy curves around each state is found numerically using an iterative approach, producing the most probable trajectory and fully characterizing the transition state and its energy. Jobs calculating such trajectories can be submitted on-line at: http://lorentz.dynstr.pasteur.fr/joel/index.php.

International audience; Low density lipoproteins (LDL) from egg yolk have a classical structure of lipoprotein with a core of neutral lipids surrounded by a monolayer of apoproteins and phospholipids. This structure collapses during adsorption and all constituents spread at the interface. To understand better the nature of the interactions between apoproteins and lipids at the interface, we have deposited LDL at an air-water interface and analysed the isotherms during their compression on a Langmuir trough. Then, these LDL films were studied by atomic force microscopy (AFM) imaging. To identify the protein and lipid structures, we imaged films before and after lipid solubilisation by butanol. To study the interactions in the LDL films, we have varied the pH, ionic strength and used simplified model systems. We also studied the correlation between observed structures and interfacial rheology of the film. The isotherms of interfacial LDL films were similar for pH 3 and 7, but their structures observed in AFM were different. At surface pressures below the transition corresponding to the demixion of apoprotein-neutral lipid complexes, the LDL film structure was not governed by electrostatic interactions. However, above this surface pressure transition (45mN/m), there was an effect of charge on this structure. Around the transition zone, the rheological properties of LDL films at pH 3 were different as a function of pH (viscous at pH 3 and visco-elastic at pH 7). So, the rheological properties of LDL films could be linked to the structures formed by apoproteins and observed in AFM.

International audience; In this review, we summarize the main structural and functional data on the role of the phosphocreatine (PCr)--creatine kinase (CK) pathway for compartmentalized energy transfer in cardiac cells. Mitochondrial creatine kinase, MtCK, fixed by cardiolipin molecules in the vicinity of the adenine nucleotide translocator, is a key enzyme in this pathway. Direct transfer of ATP and ADP between these proteins has been revealed both in experimental studies on the kinetics of the regulation of mitochondrial respiration and by mathematical modelling as a main mechanism of functional coupling of PCr production to oxidative phosphorylation. In cells in vivo or in permeabilized cells in situ, this coupling is reinforced by limited permeability of the outer membrane of the mitochondria for adenine nucleotides due to the contacts with cytoskeletal proteins. Due to these mechanisms, at least 80% of total energy is exported from mitochondria by PCr molecules. Mathematical modelling of intracellular diffusion and energy transfer shows that the main function of the PCr-CK pathway is to connect different pools (compartments) of ATP and, by this way, to overcome the local restrictions and diffusion limitation of adenine nucleotides due to the high degree of structural organization of cardiac cells.

International audience; In bacteria, P1-type ATPases are responsible for resistance to di- and monovalent toxic heavy metals by taking them out of the cell. These ATPases have a cytoplasmic N terminus comprising metal binding domains defined by a betaalphabetabetaalphabeta fold and a CXXC metal binding motif. To check how the structural properties of the metal binding site in the N terminus can influence the metal specificity of the ATPase, the first structure of a Cd(II)-ATPase N terminus was determined by NMR and its coordination sphere was investigated by X-ray absorption spectroscopy. A novel metal binding environment was found, comprising the two conserved Cys residues of the metal binding motif and a Glu in loop 5. A bioinformatic search identifies an ensemble of highly homologous sequences presumably with the same function. Another group of highly homologous sequences is found which can be referred to as zinc-detoxifying P1-type ATPases with the metal binding pattern DCXXC in the N terminus. Because no carboxylate groups participate in Cu(I) or Ag(I) binding sites, we suggest that the acidic residue plays a key role in the coordination properties of divalent cations, hence conferring a function to the N terminus in the metal specificity of the ATPase.

International audience; Specific recognition of their cognate amino acid substrates by the aminoacyl-tRNA synthetase enzymes is essential for the correct translation of the genetic code. For aspartyl-tRNA synthetase (AspRS), electrostatic interactions are expected to play an important role, since its three substrates (aspartate, ATP, tRNA) are all electrically charged. We used molecular-dynamics free-energy simulations and experiments to compare the binding of the substrate Asp and its electrically neutral analogue Asn to AspRS. The preference for Asp is found to be very strong, with good agreement between simulations and experiment. The simulations reveal long-range interactions that electrostatically couple the amino acid ligand, ATP, and its associated Mg2+ cations, a histidine side chain (His448) next to the amino acid ligand and a flexible loop that closes over the active site in response to amino acid binding. Closing this loop brings a negatively charged glutamate into the active site; this causes His448 to recruit a labile proton, which interacts favorably with Asp and accounts for most of the Asp/Asn discrimination. Cobinding of the second substrate, ATP, increases specificity for Asp further and makes the system robust towards removal of His448, which is mutated to a neutral amino acid in many organisms. Thus, AspRS specificity is assisted by a labile proton and a cosubstrate, and ATP acts as a mobile discriminator for specific Asp binding to AspRS. In asparaginyl-tRNA synthetase, a close homologue of AspRS, a few binding-pocket differences modify the charge balance so that asparagine binding predominates.

International audience; Increasingly complex schemes for representing solvent effects in an implicit fashion are being used in computational analyses of biological macromolecules. These schemes speed up the calculations by orders of magnitude and are assumed to compromise little on essential features of the solvation phenomenon. In this work we examine this assumption. Five implicit solvation models, a surface area-based empirical model, two models that approximate the generalized Born treatment and a finite difference Poisson-Boltzmann method are challenged in situations differing from those where these models were calibrated. These situations are encountered in automatic protein design procedures, whose job is to select sequences, which stabilize a given protein 3D structure, from a large number of alternatives. To this end we evaluate the energetic cost of burying amino acids in thousands of environments with different solvent exposures belonging, respectively, to decoys built with random sequences and to native protein crystal structures. In addition we perform actual sequence design calculations. Except for the crudest surface area-based procedure, all the tested models tend to favor the burial of polar amino acids in the protein interior over nonpolar ones, a behavior that leads to poor performance in protein design calculations. We show, on the other hand, that three of the examined models are nonetheless capable of discriminating between the native fold and many nonnative alternatives, a test commonly used to validate force fields. It is concluded that protein design is a particularly challenging test for implicit solvation models because it requires accurate estimates of the solvation contribution of individual residues. This contrasts with native recognition, which depends less on solvation and more on other nonbonded contributions.

Background Lectins are carbohydrate-binding proteins which potentially bind to cell surface glycoconjugates. They are found in various organisms including fungi. A lectin from the mushroom Xerocomus chrysenteron (XCL) has been isolated recently. It shows insecticidal activity and has antiproliferative properties. Results As the monosaccharide binding specificity is an important determinant of lectin function, we determined the affinity of XCL for the galactose moiety. Isothermal titration calorimetry studies revealed a dissociation constant Kd of 5.2 μM for the XCL:N-acetylgalactosamine interaction at 27degreesC. Higher affinities were observed at lower temperatures and higher osmotic pressures. The dissociation constant was five hundred times higher for the disaccharide beta-D-Gal(1–3)-D-GalNAc, Thomsen-Friedenreich (TF) antigen (Kd of 0.94 μM). By using fetuin and asialofetuin in interaction with the XCL, we revealed its ability to recognize the Thomsen-Friedenreich motif on glycoproteins. Conclusion The XCL antiproliferative effect and the TF antigen specificity presented in this work suggest that XCL and ABL may have similar binding mechanisms. The recent structure determination of these two proteins lead us to analyse these interactions in the light of our thermodynamic data. The understanding of this type of interaction may be a useful tool for the regulation of cell proliferation.

Optical melting was used to determine the stabilities of three series of RNA oligomers containing tandem XU base pairs, GGCXUGCC (5'XU3'), GGCUXGCC (5'UX3') and GGCXXGGC/CCGUUCCG (5'XX3'), where X is either A, G or I (inosine). The helices containing tandem AU base pairs were the most stable in the first two series (5'XU3' and 5'UX3'), with an average melting temperature approximately 11 degrees C higher than the helices with tandem 5'GU3' base pairs and 25 degrees C higher than the helices with tandem 5'IU3' base pairs. For the third series (5'XX3'), the helix containing tandem GG is the most stable, with an average melting temperature approximately 2 degrees C higher than the helix with tandem AA base pairs and approximately 24 degrees C higher than the helix with tandem II base pairs. The thermodynamic stability of the oligomers with tandem IU base pairs was also investigated as a function of magnesium ion concentration. As with normal A-U or G-U tandem duplexes, the data could best be interpreted as non-specific binding of magnesium ions to the inosine-containing RNA oligonucleotides.

International audience; beta-Carotene is, when cleaved, an important source of flavor and aroma compounds in fruits and flowers. Among these aroma compounds, the main degradation products are beta-ionone, 5,6-epoxy-beta-ionone, and dihydroactinidiolide (DHA), which are associated by flavorists and perfumers with fruity, floral, and woody notes. These three species can be produced by degradation of beta-carotene through an attack by enzyme-generated free radicals and a cleavage at the C9-C10 bond. This study investigated the influence of cis/trans isomerism at the C9-C10 bond on the production of beta-carotene degradation compounds, first with a predictive approach and then experimentally with different isomer mixtures. beta-Carotene solutions containing high ratios of 9-cis-isomers produced more DHA, suggesting a different pathway than for the transformation of all-trans-beta-carotene to ionone and DHA. These results are important in the search for financially viable processes to produce natural carotene-derived aroma compounds.