Beyond Poisson-Boltzmann: Numerical sampling of charge density fluctuations

We present a method aimed at sampling charge density fluctuations in Coulomb systems. The derivation follows from a functional integral representation of the partition function in terms of charge density fluctuations. Starting from the mean-field solution given by the Poisson-Boltzmann equation, an original approach is proposed to numerically sample fluctuations around it, through the propagation of a Langevin like stochastic partial differential equation (SPDE). The diffusion tensor of the SPDE can be chosen so as to avoid the numerical complexity linked to long-range Coulomb interactions, effectively rendering the theory completely local. A finite-volume implementation of the SPDE is described, and the approach is illustrated with preliminary results on the study of a system made of two like-charge ions immersed in a bath of counter-ions

Dealing with structural variability in molecular replacement and crystallographic refinement through normal-mode analysis

The possibility of taking into account large-amplitude collective movements of a given model by using a subset of low-frequency normal modes is evaluated for molecular replacement and refinement using X-ray data or cryo-EM maps

NOMAD-Ref: visualization, deformation and refinement of macromolecular structures based on all-atom normal mode analysis

Normal mode analysis (NMA) is an efficient way to study collective motions in biomolecules that bypasses the computational costs and many limitations associated with full dynamics simulations. The NOMAD-Ref web server presented here provides tools for online calculation of the normal modes of large molecules (up to 100 000 atoms) maintaining a full all-atom representation of their structures, as well as access to a number of programs that utilize these collective motions for deformation and refinement of biomolecular structures. Applications include the generation of sets of decoys with correct stereochemistry but arbitrary large amplitude movements, the quantification of the overlap between alternative conformations of a molecule, refinement of structures against experimental data, such as X-ray diffraction structure factors or Cryo-EM maps and optimization of docked complexes by modeling receptor/ligand flexibility through normal mode motions. The server can be accessed at the URL

Insights into the origin of carbonaceous chondrite organics from their triple oxygen isotope composition

International audienceDust grains of organic matter were the main reservoir of C and N in the forming Solar System and are thus considered to be an essential ingredient for the emergence of life. However, the physical environment and the chemical mechanisms at the origin of these organic grains are still highly debated. In this study, we report high-precision triple oxygen isotope composition for insoluble organic matter isolated from three emblematic carbonaceous chondrites, Orgueil, Murchison, and Cold Bokkeveld. These results suggest that the O isotope composition of carbonaceous chondrite insoluble organic matter falls on a slope 1 correlation line in the triple oxygen isotope diagram. The lack of detectable mass-dependent O isotopic fractionation, indicated by the slope 1 line, suggests that the bulk of carbonaceous chondrite organics did not form on asteroi-dal parent bodies during low-temperature hydrothermal events. On the other hand, these O isotope data, together with the H and N isotope characteristics of insoluble organic matter, may indicate that parent bodies of different carbonaceous chondrite types largely accreted organics formed locally in the protosolar nebula, possibly by photochemical dissociation of C-rich precursors

Structural Basis for Potentiation by Alcohols and Anaesthetics in a Ligand-gated Ion Channel

Ethanol alters nerve signalling by interacting with proteins in the central nervous system, particularly pentameric ligand-gated ion channels. A recent series of mutagenesis experiments on Gloeobacter violaceus ligand-gated ion channel, a prokaryotic member of this family, identified a single-site variant that is potentiated by pharmacologically relevant concentrations of ethanol. Here we determine crystal structures of the ethanol-sensitized variant in the absence and presence of ethanol and related modulators, which bind in a transmembrane cavity between channel subunits and may stabilize the open form of the channel. Structural and mutagenesis studies defined overlapping mechanisms of potentiation by alcohols and anaesthetics via the inter-subunit cavity. Furthermore, homology modelling show this cavity to be conserved in human ethanol-sensitive glycine and GABA(A) receptors, and to involve residues previously shown to influence alcohol and anaesthetic action on these proteins. These results suggest a common structural basis for ethanol potentiation of an important class of targets for neurological actions of ethanol

MinActionPath: maximum likelihood trajectory for large-scale structural transitions in a coarse-grained locally harmonic energy landscape

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

AquaSAXS: a web server for computation and fitting of SAXS profiles with non-uniformally hydrated atomic models

Small Angle X-ray Scattering (SAXS) techniques are becoming more and more useful for structural biologists and biochemists, thanks to better access to dedicated synchrotron beamlines, better detectors and the relative easiness of sample preparation. The ability to compute the theoretical SAXS profile of a given structural model, and to compare this profile with the measured scattering intensity, yields crucial structural informations about the macromolecule under study and/or its complexes in solution. An important contribution to the profile, besides the macromolecule itself and its solvent-excluded volume, is the excess density due to the hydration layer. AquaSAXS takes advantage of recently developed methods, such as AquaSol, that give the equilibrium solvent density map around macromolecules, to compute an accurate SAXS/WAXS profile of a given structure and to compare it to the experimental one. Here, we describe the interface architecture and capabilities of the AquaSAXS web server (http://lorentz.dynstr.pasteur.fr/aquasaxs.php)

Enhanced Amino Acid Selection in Fully Evolved Tryptophanyl-tRNA Synthetase, Relative to Its Urzyme, Requires Domain Motion Sensed by the D1 Switch, a Remote Dynamic Packing Motif

We previously showed (Li, L., and Carter, C. W., Jr. (2013) J. Biol. Chem. 288, 34736–34745) that increased specificity for tryptophan versus tyrosine by contemporary Bacillus stearothermophilus tryptophanyl-tRNA synthetase (TrpRS) over that of TrpRS Urzyme results entirely from coupling between the anticodon-binding domain and an insertion into the Rossmann-fold known as Connecting Peptide 1. We show that this effect is closely related to a long range catalytic effect, in which side chain repacking in a region called the D1 Switch, accounts fully for the entire catalytic contribution of the catalytic Mg2+ ion. We report intrinsic and higher order interaction effects on the specificity ratio, (kcat/Km)Trp/(kcat/Km)Tyr, of 15 combinatorial mutants from a previous study (Weinreb, V., Li, L., and Carter, C. W., Jr. (2012) Structure 20, 128–138) of the catalytic role of the D1 Switch. Unexpectedly, the same four-way interaction both activates catalytic assist by Mg2+ ion and contributes −4.4 kcal/mol to the free energy of the specificity ratio. A minimum action path computed for the induced-fit and catalytic conformation changes shows that repacking of the four residues precedes a decrease in the volume of the tryptophan-binding pocket. We suggest that previous efforts to alter amino acid specificities of TrpRS and glutaminyl-tRNA synthetase (GlnRS) by mutagenesis without extensive, modular substitution failed because mutations were incompatible with interdomain motions required for catalysis

Structural Studies of HNA Substrate Specificity in Mutants of an Archaeal DNA Polymerase Obtained by Directed Evolution.

Archaeal DNA polymerases from the B-family (polB) have found essential applications in biotechnology. In addition, some of their variants can accept a wide range of modified nucleotides or xenobiotic nucleotides, such as 1,5-anhydrohexitol nucleic acid (HNA), which has the unique ability to selectively cross-pair with DNA and RNA. This capacity is essential to allow the transmission of information between different chemistries of nucleic acid molecules. Variants of the archaeal polymerase from Thermococcus gorgonarius, TgoT, that can either generate HNA from DNA (TgoT_6G12) or DNA from HNA (TgoT_RT521) have been previously identified. To understand how DNA and HNA are recognized and selected by these two laboratory-evolved polymerases, we report six X-ray structures of these variants, as well as an in silico model of a ternary complex with HNA. Structural comparisons of the apo form of TgoT_6G12 together with its binary and ternary complexes with a DNA duplex highlight an ensemble of interactions and conformational changes required to promote DNA or HNA synthesis. MD simulations of the ternary complex suggest that the HNA-DNA hybrid duplex remains stable in the A-DNA helical form and help explain the presence of mutations in regions that would normally not be in contact with the DNA if it were not in the A-helical form. One complex with two incorporated HNA nucleotides is surprisingly found in a one nucleotide-backtracked form, which is new for a DNA polymerase. This information can be used for engineering a new generation of more efficient HNA polymerase variants

Conferring a template-dependent polymerase activity to terminal deoxynucleotidyltransferase by mutations in the Loop1 region

Terminal deoxynucleotidyltransferase (Tdt) and DNA polymerase μ (pol μ) are two eukaryotic highly similar proteins involved in DNA processing and repair. Despite their high sequence identity, they differ widely in their activity: pol μ has a templated polymerase activity, whereas Tdt has a non-templated one. Loop1, first described when the Tdt structure was solved, has been invoked as the major structural determinant of this difference. Here we describe attempts to transform Tdt into pol μ with the minimal number of mutations in and around Loop1. First we describe the effect of mutations on six different positions chosen to destabilize Tdt Loop1 structure, either by alanine substitution or by deletion; they result at most in a reduction of Tdt activity, but adding Co++ restores most of this Tdt activity. However, a deletion of the entire Loop1 as in pol λ does confer a limited template-dependent polymerase behavior to Tdt while a chimera bearing an extended pol μ Loop1 reproduces pol μ behavior. Finally, 16 additional substitutions are reported, targeted at the two so-called ‘sequence determinant’ regions located just after Loop1 or underneath. Among them, the single-point mutant F401A displays a sequence-specific replicative polymerase phenotype that is stable upon Co++ addition. These results are discussed in light of the available crystal structures