Structural basis for a novel mechanism of DNA bridging and alignment in eukaryotic DSB DNA repair

Eukaryotic DNA polymerase mu of the PolX family can promote the association of the two 3'-protruding ends of a DNA double-strand break (DSB) being repaired (DNA synapsis) even in the absence of the core non-homologous end-joining (NHEJ) machinery. Here, we show that terminal deoxynucleotidyltransferase (TdT), a closely related PolX involved in V(D)J recombination, has the same property. We solved its crystal structure with an annealed DNA synapsis containing one micro-homology (MH) base pair and one nascent base pair. This structure reveals how the N-terminal domain and Loop 1 of Tdt cooperate for bridging the two DNA ends, providing a templating base in trans and limiting the MH search region to only two base pairs. A network of ordered water molecules is proposed to assist the incorporation of any nucleotide independently of the in trans templating base. These data are consistent with a recent model that explains the statistics of sequences synthesized in vivo by Tdt based solely on this dinucleotide step. Site-directed mutagenesis and functional tests suggest that this structural model is also valid for Pol mu during NHEJ.

Genuine open form of the pentameric ligand-gated ion channel GLIC

Pentameric ligand-gated ion channels (pLGICs) mediate fast chemical neurotransmission of nerve signalling in the central and peripheral nervous systems. GLIC is a bacterial homologue of eukaryotic pLGIC, the X-ray structure of which has been determined in three different conformations. GLIC is thus widely used as a model to study the activation and the allosteric transition of this family of receptors. The recently solved high-resolution structure of GLIC (2.4 Å resolution) in the active state revealed two bound acetate molecules in the extracellular domain (ECD). Here, it is shown that these two acetates exactly overlap with known sites of pharmacological importance in pLGICs, and their potential influence on the structure of the open state is studied in detail. Firstly, experimental evidence is presented for the correct assignment of these acetate molecules by using the anomalous dispersion signal of bromoacetate. Secondly, the crystal structure of GLIC in the absence of acetate was solved and it is shown that acetate binding induces local conformational changes that occur in strategic sites of the ECD. It is expected that this acetate-free structure will be useful in future computational studies of the gating transition in GLIC and other pLGICs.

Modified Poisson–Boltzmann equations for characterizing biomolecular solvation

Methods for computing electrostatic interactions often account implicitly for the solvent, due to the much smaller number of degrees of freedom involved. In the Poisson–Boltzmann (PB) approach the electrostatic potential is obtained by solving the Poisson–Boltzmann equation (PBE), where the solvent region is modeled as a homogeneous medium with a high dielectric constant. PB however is not exempt of problems. It does not take into account for example the sizes of the ions in the atmosphere surrounding the solute, nor does it take into account the inhomogeneous dielectric response of water due to the presence of a highly charged surface. In this paper we review two major modifications of PB that circumvent these problems, namely the size-modified PB (SMPB) equation and the Dipolar Poisson–Boltzmann Langevin (DPBL) model. In SMPB, steric effects between ions are accounted for with a lattice gas model. In DPBL, the solvent region is no longer modeled as a homogeneous dielectric media but rather as an assembly of self-orienting interacting dipoles of variable density. This model results in a dielectric profile that transits smoothly from the solute to the solvent region as well as in a variable solvent density that depends on the charges of the solute. We show successful applications of the DPBL formalism to computing the solvation free energies of isolated ions in water. Further developments of more accurately modified PB models are discussed.

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 Mg(2+) 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 Mg(2+) 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.

New nucleotide-competitive non-nucleoside inhibitors of terminal deoxynucleotidyl transferase: discovery, characterization, and crystal structure in complex with the target

Terminal deoxynucletidyl transferase (TdT) is overexpressed in some cancer types, where it might compete with pol μ during the mutagenic repair of double strand breaks (DSBs) through the nonhomologous end joining (NHEJ) pathway. Here we report the discovery and characterization of pyrrolyl and indolyl diketo acids that specifically target TdT and behave as nucleotide-competitive inhibitors. These compounds show a selective toxicity toward MOLT-4 compared to HeLa cells that correlate well with in vitro selectivity for TdT. The binding site of two of these inhibitors was determined by cocrystallization with TdT, explaining why these compounds are competitive inhibitors of the deoxynucleotide triphosphate (dNTP). In addition, because of the observed dual localization of the phenyl substituent, these studies open the possibility of rationally designing more potent compounds.

Computational assembly of polymorphic amyloid fibrils reveals stable aggregates

Amyloid proteins aggregate into polymorphic fibrils that damage tissues of the brain, nerves, and heart. Experimental and computational studies have examined the structural basis and the nucleation of short fibrils, but the ability to predict and precisely quantify the stability of larger aggregates has remained elusive. We established a complete classification of fibril shapes and developed a tool called CreateFibril to build such complex, polymorphic, modular structures automatically. We applied stability landscapes, a technique we developed to reveal reliable fibril structural parameters, to assess fibril stability. CreateFibril constructed HET-s, Aβ, and amylin fibrils up to 17 nm in length, and utilized a novel dipolar solvent model that captured the effect of dipole-dipole interactions between water and very large molecular systems to assess their aqueous stability. Our results validate experimental data for HET-s and Aβ, and suggest novel (to our knowledge) findings for amylin. In particular, we predicted the correct structural parameters (rotation angles, packing distances, hydrogen bond lengths, and helical pitches) for the one and three predominant HET-s protofilaments. We reveal and structurally characterize all known Aβ polymorphic fibrils, including structures recently classified as wrapped fibrils. Finally, we elucidate the predominant amylin fibrils and assert that native amylin is more stable than its amyloid form. CreateFibril and a database of all stable polymorphic fibril models we tested, along with their structural energy landscapes, are available at http://amyloid.cs.mcgill.ca.

Structures of intermediates along the catalytic cycle of terminal deoxynucleotidyltransferase: dynamical aspects of the two-metal ion mechanism

Terminal deoxynucleotidyltransferase (Tdt) is a non-templated eukaryotic DNA polymerase of the polX family that is responsible for the random addition of nucleotides at the V(D)J junctions of immunoglobulins and T-cell receptors. Here we describe a series of high-resolution X-ray structures that mimic the pre-catalytic state, the post-catalytic state and a competent state that can be transformed into the two other ones in crystallo via the addition of dAMPcPP and Zn(2+), respectively. We examined the effect of Mn(2+), Co(2+) and Zn(2+) because they all have a marked influence on the kinetics of the reaction. We demonstrate a dynamic role of divalent transition metal ions bound to site A: (i) Zn(2+) (or Co(2+)) in Metal A site changes coordination from octahedral to tetrahedral after the chemical step, which explains the known higher affinity of Tdt for the primer strand when these ions are present, and (ii) metal A has to leave to allow the translocation of the primer strand and to clear the active site, a typical feature for a ratchet-like mechanism. Except for Zn(2+), the sugar puckering of the primer strand 3' terminus changes from C2'-endo to C3'-endo during catalysis. In addition, our data are compatible with a scheme where metal A is the last component that binds to the active site to complete its productive assembly, as already inferred in human pol beta. The new structures have potential implications for modeling pol mu, a closely related polX implicated in the repair of DNA double-strand breaks, in a complex with a DNA synapsis.

Structure and pharmacology of pentameric receptor channels: from bacteria to brain

Orthologs of the pentameric receptor channels that mediate fast synaptic transmission in the central and peripheral nervous systems have been found in several bacterial species and in a single archaea genus. Recent X-ray structures of bacterial and invertebrate pentameric receptors point to a striking conservation of the structural features within the whole family, even between distant prokaryotic and eukaryotic members. These structural data reveal general principles of molecular organization that allow allosteric membrane proteins to mediate chemoelectric transduction. Notably, several conformations have been solved, including open and closed channels with distinct global tertiary and quaternary structure. The data reveal features of the ion channel architecture and of diverse categories of binding sites, such as those that bind orthosteric ligands, including neurotransmitters, and those that bind allosteric modulators, such as general anesthetics, ivermectin, or lipids. In this review, we summarize the most recent data, discuss insights into the mechanism of action in these systems, and elaborate on newly opened avenues for drug design.

2'-Deoxyribonucleoside phosphoramidate triphosphate analogues as alternative substrates for E. coli polymerase III

Thermostable bacterial polymerases like Taq, Therminator and Vent exo(-) are able to perform DNA synthesis by using modified DNA precursors, a property that is exploited in several therapeutic and biotechnological applications. Viral polymerases are also known to accept modified substrates, and this has proven crucial in the development of antiviral therapies. However, non-thermostable polymerases of bacterial origin, or engineered variants, that have similar substrate tolerance and could be used for synthetic biology purposes remain to be identified. We have identified the α subunit of Escherichia coli polymerase III (Pol III α) as a bacterial polymerase that is able to recognise and process as substrates several pyrophosphate-modified dATP analogues in place of its natural substrate dATP for template-directed DNA synthesis. A number of dATP analogues featuring a modified pyrophosphate group were able to serve as substrates during enzymatic DNA synthesis by Pol III α. Features such as the presence of potentially chelating chemical groups and the size and spatial flexibility of the chemical structure seem to be of major importance for the modified leaving group to play its role during the enzymatic reaction. In addition, we could establish that if the pyrophosphate group is altered, deoxynucleotide incorporation proceeds with an efficiency varying with the nature of the nucleobase. Our results represent a great step towards the achievement of a system of artificial DNA synthesis hosted by E. coli and involving the use of altered nucleotide precursors for nucleic acid synthesis.

Molecular recognition of canonical and deaminated bases by P. abyssi family B DNA polymerase

Euryarchaeal polymerase B can recognize deaminated bases on the template strand, effectively stalling the replication fork 4nt downstream the modified base. Using Pyrococcus abyssi DNA B family polymerase (PabPolB), we investigated the discrimination between deaminated and natural nucleotide(s) by primer extension assays, electrophoretic mobility shift assays, and X-ray crystallography. Structures of complexes between the protein and DNA duplexes with either a dU or a dH in position +4 were solved at 2.3Å and 2.9Å resolution, respectively. The PabPolB is found in the editing mode. A new metal binding site has been uncovered below the base-checking cavity where the +4 base is flipped out; it is fully hydrated in an octahedral fashion and helps guide the strongly kinked template strand. Four other crystal structures with each of the canonical bases were also solved in the editing mode, and the presence of three nucleotides in the exonuclease site caused a shift in the coordination state of its metal A from octahedral to tetrahedral. Surprisingly, we find that all canonical bases also enter the base-checking pocket with very small differences in the binding geometry and in the calculated binding free energy compared to deaminated ones. To explain how this can lead to stalling of the replication fork, the full catalytic pathway and its branches must be taken into account, during which the base is checked several times. Our results strongly suggest a switch from elongation to editing modes right after nucleotide insertion when the modified base is at position +5.

A locally closed conformation of a bacterial pentameric proton-gated ion channel

Pentameric ligand-gated ion channels mediate signal transduction through conformational transitions between closed-pore and open-pore states. To stabilize a closed conformation of GLIC, a bacterial proton-gated homolog from Gloeobacter violaceus whose open structure is known, we separately generated either four cross-links or two single mutations. We found all six mutants to be in the same 'locally closed' conformation using X-ray crystallography, sharing most of the features of the open form but showing a locally closed pore as a result of a concerted bending of all of its M2 helices. The mutants adopt several variant conformations of the M2-M3 loop, and in all cases an interacting lipid that is observed in the open form disappears. A single cross-linked mutant is functional, according to electrophysiology, and the locally closed structure of this mutant indicates that it has an increased flexibility. Further cross-linking, accessibility and molecular dynamics data suggest that the locally closed form is a functionally relevant conformation that occurs during allosteric gating transitions.

Adapting Poisson-Boltzmann to the self-consistent mean field theory: application to protein side-chain modeling

We present an extension of the self-consistent mean field theory for protein side-chain modeling in which solvation effects are included based on the Poisson-Boltzmann (PB) theory. In this approach, the protein is represented with multiple copies of its side chains. Each copy is assigned a weight that is refined iteratively based on the mean field energy generated by the rest of the protein, until self-consistency is reached. At each cycle, the variational free energy of the multi-copy system is computed; this free energy includes the internal energy of the protein that accounts for vdW and electrostatics interactions and a solvation free energy term that is computed using the PB equation. The method converges in only a few cycles and takes only minutes of central processing unit time on a commodity personal computer. The predicted conformation of each residue is then set to be its copy with the highest weight after convergence. We have tested this method on a database of hundred highly refined NMR structures to circumvent the problems of crystal packing inherent to x-ray structures. The use of the PB-derived solvation free energy significantly improves prediction accuracy for surface side chains. For example, the prediction accuracies for χ(1) for surface cysteine, serine, and threonine residues improve from 68%, 35%, and 43% to 80%, 53%, and 57%, respectively. A comparison with other side-chain prediction algorithms demonstrates that our approach is consistently better in predicting the conformations of exposed side chains.

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

Structural insights into the quinolone resistance mechanism of Mycobacterium tuberculosis DNA gyrase

Mycobacterium tuberculosis DNA gyrase, an indispensable nanomachine involved in the regulation of DNA topology, is the only type II topoisomerase present in this organism and is hence the sole target for quinolone action, a crucial drug active against multidrug-resistant tuberculosis. To understand at an atomic level the quinolone resistance mechanism, which emerges in extensively drug resistant tuberculosis, we performed combined functional, biophysical and structural studies of the two individual domains constituting the catalytic DNA gyrase reaction core, namely the Toprim and the breakage-reunion domains. This allowed us to produce a model of the catalytic reaction core in complex with DNA and a quinolone molecule, identifying original mechanistic properties of quinolone binding and clarifying the relationships between amino acid mutations and resistance phenotype of M. tuberculosis DNA gyrase. These results are compatible with our previous studies on quinolone resistance. Interestingly, the structure of the entire breakage-reunion domain revealed a new interaction, in which the Quinolone-Binding Pocket (QBP) is blocked by the N-terminal helix of a symmetry-related molecule. This interaction provides useful starting points for designing peptide based inhibitors that target DNA gyrase to prevent its binding to DNA.

Atomic structure and dynamics of pentameric ligand-gated ion channels: new insight from bacterial homologues

Pentameric ligand-gated ion channels (pLGICs) are widely expressed in the animal kingdom and are key players of neurotransmission by acetylcholine (ACh), gamma-amminobutyric acid (GABA), glycine and serotonin. It is now established that this family has a prokaryotic origin, since more than 20 homologues have been discovered in bacteria. In particular, the GLIC homologue displays a ligand-gated ion channel function and is activated by protons. The prokaryotic origin of these membrane proteins facilitated the X-ray structural resolution of the first members of this family. ELIC was solved at 3.3 A in a closed-pore conformation, and GLIC at up to 2.9 A in an apparently open-pore conformation. These data reveal many structural features, notably the architecture of the pore, including its gate and its selectivity filter, and the interactions between the protein and lipids. In addition, comparison of the structures of GLIC and ELIC hints at a mechanism of channel opening, which consists of both a quaternary twist and a tertiary deformation. This mechanism couples opening-closing motions of the channel with a global reorganization of the protein, including the subunit interface that holds the neurotransmitter binding sites in eukaryotic pLGICs.

AQUASOL: An efficient solver for the dipolar Poisson-Boltzmann-Langevin equation

The Poisson-Boltzmann (PB) formalism is among the most popular approaches to modeling the solvation of molecules. It assumes a continuum model for water, leading to a dielectric permittivity that only depends on position in space. In contrast, the dipolar Poisson-Boltzmann-Langevin (DPBL) formalism represents the solvent as a collection of orientable dipoles with nonuniform concentration; this leads to a nonlinear permittivity function that depends both on the position and on the local electric field at that position. The differences in the assumptions underlying these two models lead to significant differences in the equations they generate. The PB equation is a second order, elliptic, nonlinear partial differential equation (PDE). Its response coefficients correspond to the dielectric permittivity and are therefore constant within each subdomain of the system considered (i.e., inside and outside of the molecules considered). While the DPBL equation is also a second order, elliptic, nonlinear PDE, its response coefficients are nonlinear functions of the electrostatic potential. Many solvers have been developed for the PB equation; to our knowledge, none of these can be directly applied to the DPBL equation. The methods they use may adapt to the difference; their implementations however are PBE specific. We adapted the PBE solver originally developed by Holst and Saied [J. Comput. Chem. 16, 337 (1995)] to the problem of solving the DPBL equation. This solver uses a truncated Newton method with a multigrid preconditioner. Numerical evidences suggest that it converges for the DPBL equation and that the convergence is superlinear. It is found however to be slow and greedy in memory requirement for problems commonly encountered in computational biology and computational chemistry. To circumvent these problems, we propose two variants, a quasi-Newton solver based on a simplified, inexact Jacobian and an iterative self-consistent solver that is based directly on the PBE solver. While both methods are not guaranteed to converge, numerical evidences suggest that they do and that their convergence is also superlinear. Both variants are significantly faster than the solver based on the exact Jacobian, with a much smaller memory footprint. All three methods have been implemented in a new code named AQUASOL, which is freely available.

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

Terminal deoxynucleotidyltransferase (Tdt) and DNA polymerase mu (pol mu) are two eukaryotic highly similar proteins involved in DNA processing and repair. Despite their high sequence identity, they differ widely in their activity: pol mu 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 mu 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 lambda does confer a limited template-dependent polymerase behavior to Tdt while a chimera bearing an extended pol mu Loop1 reproduces pol mu 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.

Computing ion solvation free energies using the dipolar Poisson model

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.

Structure of the archaeal pab87 peptidase reveals a novel self-compartmentalizing protease family

Self-compartmentalizing proteases orchestrate protein turnover through an original architecture characterized by a central catalytic chamber. Here we report the first structure of an archaeal member of a new self-compartmentalizing protease family forming a cubic-shaped octamer with D₄ symmetry and referred to as CubicO. We solved the structure of the Pyrococcus abyssi Pab87 protein at 2.2 Å resolution using the anomalous signal of the high-phasing-power lanthanide derivative Lu-HPDO3A. A 20 Å wide channel runs through this supramolecular assembly of 0.4 MDa, giving access to a 60 Å wide central chamber holding the eight active sites. Surprisingly, activity assays revealed that Pab87 degrades specifically d-amino acid containing peptides, which have never been observed in archaea. Genomic context of the Pab87 gene showed that it is surrounded by genes involved in the amino acid/peptide transport or metabolism. We propose that CubicO proteases are involved in the processing of d-peptides from environmental origins.

Beyond the Poisson-Boltzmann model: modeling biomolecule-water and water-water interactions

We present an extension to the Poisson-Boltzmann model in which the solvent is modeled as an assembly of self-orienting dipoles of variable densities. Interactions between these dipoles are included implicitly using a Yukawa potential field. This model leads to a set of equations whose solutions give the dipole densities; we use the latter to study the organization of water around biomolecules. The computed water density profiles resemble those derived from molecular dynamics simulations. We also derive an excess free energy that discriminates correct from incorrect conformations of proteins.

Incorporating dipolar solvents with variable density in Poisson-Boltzmann electrostatics

We describe a new way to calculate the electrostatic properties of macromolecules that goes beyond the classical Poisson-Boltzmann treatment with only a small extra CPU cost. The solvent region is no longer modeled as a homogeneous dielectric media but rather as an assembly of self-orienting interacting dipoles of variable density. The method effectively unifies both the Poisson-centric view and the Langevin Dipole model. The model results in a variable dielectric constant epsilon(r) in the solvent region and also in a variable solvent density rho(r) that depends on the nature of the closest exposed solute atoms. The model was calibrated using small molecules and ions solvation data with only two adjustable parameters, namely the size and dipolar moment of the solvent. Hydrophobicity scales derived from the solvent density profiles agree very well with independently derived hydrophobicity scales, both at the atomic or residue level. Dimerization interfaces in homodimeric proteins or lipid-binding regions in membrane proteins clearly appear as poorly solvated patches on the solute accessible surface. Comparison of the thermally averaged solvent density of this model with the one derived from molecular dynamics simulations shows qualitative agreement on a coarse-grained level. Because this calculation is much more rapid than that from molecular dynamics, applications of a density-profile-based solvation energy to the identification of the true structure among a set of decoys become computationally feasible. Various possible improvements of the model are discussed, as well as extensions of the formalism to treat mixtures of dipolar solvents of different sizes.

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

Normal-mode analysis (NMA) can be used to generate multiple structural variants of a given template model, thereby increasing the chance of finding the molecular-replacement solution. Here, it is shown that it is also possible to directly refine the amplitudes of the normal modes against experimental data (X-ray or cryo-EM), generalizing rigid-body refinement methods by adding just a few additional degrees of freedom that sample collective and large-amplitude movements. It is also argued that the situation where several (conformations of) models are present simultaneously in the crystal can be studied with adjustable occupancies using techniques derived from statistical thermodynamics and already used in molecular modelling.

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.

An asymmetric underlying rule in the assignment of codons: possible clue to a quick early evolution of the genetic code via successive binary choices

Aminoacyl-tRNA synthetases (aaRSs) are responsible for creating the pool of correctly charged aminoacyl-tRNAs that are necessary for the translation of genetic information (mRNA) by the ribosome. Each aaRS belongs to either one of only two classes with two different mechanisms of aminoacylation, making use of either the 2'OH (Class I) or the 3'OH (Class II) of the terminal A76 of the tRNA and approaching the tRNA either from the minor groove (2'OH) or the major groove (3'OH). Here, an asymmetric pattern typical of differentiation is uncovered in the partition of the codon repertoire, as defined by the mechanism of aminoacylation of each corresponding tRNA. This pattern can be reproduced in a unique cascade of successive binary decisions that progressively reduces codon ambiguity. The deduced order of differentiation is manifestly driven by the reduction of translation errors. A simple rule can be defined, decoding each codon sequence in its binary class, thereby providing both the code and the key to decode it. Assuming that the partition into two mechanisms of tRNA aminoacylation is a relic that dates back to the invention of the genetic code in the RNA World, a model for the assignment of amino acids in the codon table can be derived. The model implies that the stop codon was always there, as the codon whose tRNA cannot be charged with any amino acid, and makes the prediction of an ultimate differentiation step, which is found to correspond to the codon assignment of the 22nd amino acid pyrrolysine in archaebacteria.

Three dimensional structure and implications for the catalytic mechanism of 6-phosphogluconolactonase from Trypanosoma brucei

Enzymes from the pentose phosphate pathway (PPP) are potential drug targets for the development of new drugs against Trypanosoma brucei, the causative agent of African sleeping disease: for instance, the 6-phosphogluconate dehydrogenase is currently studied actively for such purposes. Structural and functional studies are necessary to better characterize the associated enzymes and compare them to their human homologues, in order to undertake structure-based drug design studies on such targets. In this context, the crystal structure of 6-phosphogluconolactonase (6PGL) from T. brucei, the second enzyme from PPP, was determined at 2.1 Angstroms resolution. Comparison of its sequence and structure to other related proteins in the 6PGL family with a known structure (Thermotoga maritima Tm6GPL 1PBT and Vibrio cholerae Vc6PGL (1Y89), which have not been discussed in print), or in the glucosamine-6-phosphate-deaminase family (hexameric Escherichia coli 1DEA and monomeric Bacillus subtilis 2BKV), allowed the identification of the 6PGL active site. In addition to the analysis of the crystal structure, 3D NMR interaction studies and docking experiments are reported here. Key residues involved in substrate binding or in catalysis were identified.

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 .

PDB_Hydro: incorporating dipolar solvents with variable density in the Poisson-Boltzmann treatment of macromolecule electrostatics

We describe a new way to calculate the electrostatic properties of macromolecules which eliminates the assumption of a constant dielectric value in the solvent region, resulting in a Generalized Poisson–Boltzmann–Langevin equation (GPBLE). We have implemented a web server () that both numerically solves this equation and uses the resulting water density profiles to place water molecules at preferred sites of hydration. Surface atoms with high or low hydration preference can be easily displayed using a simple PyMol script, allowing for the tentative prediction of the dimerization interface in homodimeric proteins, or lipid binding regions in membrane proteins. The web site includes options that permit mutations in the sequence as well as reconstruction of missing side chain and/or main chain atoms. These tools are accessible independently from the electrostatics calculation, and can be used for other modeling purposes. We expect this web server to be useful to structural biologists, as the knowledge of solvent density should prove useful to get better fits at low resolution for X-ray diffraction data and to computational biologists, for whom these profiles could improve the calculation of interaction energies in water between ligands and receptors in docking simulations.

Determination of dihedral Ψ angles in large proteins by combining NHN/CαHα dipole/dipole cross-correlation and chemical shifts

We propose a strategy based on the combination of experimental NH(N)/C(alpha)H(alpha) dipole/dipole cross-correlated relaxation rates and chemical shift analysis for the determination of Psi torsion angles in proteins. The method allows the determination of a dihedral angle that is not easily accessible by nuclear magnetic resonance (NMR). The measurement of dihedral angle restraints can be used for structure calculation, which is known to improve the quality of NMR structures. The method is of particular interest in the case of large proteins, for which spectral assignment of the nuclear Overhauser effect spectra, and therefore straightforward structural determination, is out of reach. One advantage of the method is that it is reasonably simple to implement, and could be used in association with other methods aiming at obtaining structural information on complex systems, such as residual dipolar coupling measurements. An illustrative example is analyzed in the case of the 30-kDa protein 6-phosphogluconolactonase.

Refinement of docked protein-ligand and protein-DNA structures using low frequency normal mode amplitude optimization

Prediction of structural changes resulting from complex formation, both in ligands and receptors, is an important and unsolved problem in structural biology. In this work, we use all-atom normal modes calculated with the Elastic Network Model as a basis set to model structural flexibility during formation of macromolecular complexes and refine the non-bonded intermolecular energy between the two partners (protein–ligand or protein–DNA) along 5–10 of the lowest frequency normal mode directions. The method handles motions unrelated to the docking transparently by first applying the modes that improve non-bonded energy most and optionally restraining amplitudes; in addition, the method can correct small errors in the ligand position when the first six rigid-body modes are switched on. For a test set of six protein receptors that show an open-to-close transition when binding small ligands, our refinement scheme reduces the protein coordinate cRMS by 0.3–3.2 Å. For two test cases of DNA structures interacting with proteins, the program correctly refines the docked B-DNA starting form into the expected bent DNA, reducing the DNA cRMS from 8.4 to 4.8 Å and from 8.7 to 5.4 Å, respectively. A public web server implementation of the refinement method is available at .

Normal mode analysis suggests a quaternary twist model for the nicotinic receptor gating mechanism

We present a three-dimensional model of the homopentameric alpha7 nicotinic acetylcholine receptor (nAChR), that includes the extracellular and membrane domains, developed by comparative modeling on the basis of: 1), the x-ray crystal structure of the snail acetylcholine binding protein, an homolog of the extracellular domain of nAChRs; and 2), cryo-electron microscopy data of the membrane domain collected on Torpedo marmorata nAChRs. We performed normal mode analysis on the complete three-dimensional model to explore protein flexibility. Among the first 10 lowest frequency modes, only the first mode produces a structural reorganization compatible with channel gating: a wide opening of the channel pore caused by a concerted symmetrical quaternary twist motion of the protein with opposing rotations of the upper (extracellular) and lower (transmembrane) domains. Still, significant reorganizations are observed within each subunit, that involve their bending at the domain interface, an increase of angle between the two beta-sheets composing the extracellular domain, the internal beta-sheet being significantly correlated to the movement of the M2 alpha-helical segment. This global symmetrical twist motion of the pentameric protein complex, which resembles the opening transition of other multimeric ion channels, reasonably accounts for the available experimental data and thus likely describes the nAChR gating process.

On the use of low-frequency normal modes to enforce collective movements in refining macromolecular structural models

As more and more structures of macromolecular complexes get solved in different conditions, it has become apparent that flexibility is an inherent part of their biological function. Normal mode analysis using simplified models of proteins such as the elastic network model has proved very effective in showing that many of the structural transitions derived from a survey of the Protein Data Bank can be explained by just a few of the lowest-frequency normal modes. In this work, normal modes are used to carry out medium- or low-resolution structural refinement, enforcing collective and large-amplitude movements that are beyond the reach of existing methods. Refinement is carried out in reciprocal space with respect to the normal mode amplitudes, by using standard conjugate-gradient minimization. Several tests on synthetic diffraction data whose mode concentration follows the one of real movements observed in the Protein Data Bank have shown that the radius of convergence is larger than the one of rigid-body refinement. Tests with experimental diffraction data for the same protein in different environments also led to refined structural models showing drastic reduction of the rms deviation with the target model. Because the structural transition is described by very few parameters, over-fitting of real experimental data is easily detected by using a cross-validation test. The method has also been applied to the refinement of atomic models into molecular envelopes and could readily be used to fit large macromolecular complex rearrangements into cryo-electron microscopy-reconstructed images as well as small-angle x-ray scattering-derived envelopes.

Mycobacterium tuberculosis thymidylate kinase: structural studies of intermediates along the reaction pathway

Mycobacterium tuberculosis TMP kinase (TMPK(Mtub)) represents a promising target for developing drugs against tuberculosis because the configuration of its active site is unique in the TMPK family. To help elucidate the phosphorylation mechanism employed by this enzyme, structural changes occurring upon binding of substrates and subsequent catalysis were investigated by protein crystallography. Six new structures of TMPK(Mtub) were solved at a resolution better than 2.3A, including the first structure of an apo-TMPK, obtained by triggering catalysis in a crystal of a TMPK(Mtub)-TMP complex, which resulted in the release of the TDP product. A series of snapshots along the reaction pathway is obtained, revealing the closure of the active site in going from an empty to a fully occupied state, suggestive of an induced-fit mechanism typical of NMPKs. However, in TMPK(Mtub) the LID closure couples to the binding with an unusual location for a magnesium ion coordinating TMP in the active site. Our data suggest strongly that this ion is required for catalysis, acting as a clamp, possibly in concert with Arg95, to neutralise electrostatic repulsion between the anionic substrates, optimise their proper alignment and activate them through direct and water-mediated interactions. The 3'-hydroxyl moiety of TMP, critical to metal stabilisation, appears to be a target of choice for the design of potent inhibitors. On the other hand, the usual NTP-bound magnesium is not seen in our structures and Arg14, a P-loop residue unique to TMPK(Mtub), may take over its role. Therefore, TMPK(Mtub) seems to have swapped the use of a metal ion as compared with e.g. human TMPK. Finally, TTP was observed in crystals of TMPK(Mtub), locked by Arg14, thus providing a structural explanation for the observed inhibitory effect of TTP putatively involved in a mechanism of feedback regulation of the enzymatic activity.

Enzymatic and structural analysis of inhibitors designed against Mycobacterium tuberculosis thymidylate kinase. New insights into the phosphoryl transfer mechanism

The chemical synthesis of new compounds designed as inhibitors of Mycobacterium tuberculosis TMP kinase (TMPK) is reported. The synthesis concerns TMP analogues modified at the 5-position of the thymine ring as well as a novel compound with a six-membered sugar ring. The binding properties of the analogues are compared with the known inhibitor azido-TMP, which is postulated here to work by excluding the TMP-bound Mg(2+) ion. The crystallographic structure of the complex of one of the compounds, 5-CH(2)OH-dUMP, with TMPK has been determined at 2.0 A. It reveals a major conformation for the hydroxyl group in contact with a water molecule and a minor conformation pointing toward Ser(99). Looking for a role for Ser(99), we have identified an unusual catalytic triad, or a proton wire, made of strictly conserved residues (including Glu(6), Ser(99), Arg(95), and Asp(9)) that probably serves to protonate the transferred PO(3) group. The crystallographic structure of the commercially available bisubstrate analogue P(1)-(adenosine-5')-P(5)-(thymidine-5')-pentaphosphate bound to TMPK is also reported at 2.45 A and reveals an alternative binding pocket for the adenine moiety of the molecule compared with what is observed either in the Escherichia coli or in the yeast enzyme structures. This alternative binding pocket opens a way for the design of a new family of specific inhibitors.

Simplified normal mode analysis of conformational transitions in DNA-dependent polymerases: the elastic network model

The Elastic Network Model is used to investigate the open/closed transition in all DNA-dependent polymerases whose structure is known in both forms. For each structure the model accounts well for experimental crystallographic B-factors. It is found in all cases that the transition can be well described with just a handful of the normal modes. Usually, only the lowest and/or the second lowest frequency normal modes deduced from the open form give rise to calculated displacement vectors that have a correlation coefficient larger than 0.50 with the observed difference vectors between the two forms. This is true for every structural class of DNA-dependent polymerases where a direct comparison with experimental structural data is available. In cases where only one form has been observed by X-ray crystallography, it is possible to make predictions concerning the possible existence of another form in solution by carefully examining the vector displacements predicted for the lowest frequency normal modes. This simple model, which has the advantage to be computationally inexpensive, could be used to design novel kind of drugs directed against polymerases, namely drugs preventing the open/closed transition from occurring in bacterial or viral DNA-dependent polymerases.

Database searching for thymidine and thymidylate kinase inhibitors using three-dimensional structure-based methods

Structure-based drug design methods were used to search for novel inhibitors of herpes simplex virus type 1 (HSV-1) thymidine kinase and Mycobacterium tuberculosis thymidylate kinase. The method involved the use of crystal structure complexes to guide database searching for potential inhibitors. A number of weak inhibitors of HSV-2 were identified, one of which was found to inhibit HSV-1 TK and HSV-1 TK-deficient viral strains. Each compound tested against M. tuberculosis thymidylate kinase was found to have some activity. The best of these compounds was only 4.6-fold less potent than 3'-azido-3'-deoxythymidine-5'-monophosphate (AZTMP). This study demonstrates the utility of structure-based drug design methods in the search for novel enzyme inhibitors.

Cryophotolysis of caged compounds: a technique for trapping intermediate states in protein crystals

Caged compounds in combination with protein crystallography represent a valuable tool in studies of enzyme reaction intermediates. To date, photochemical triggering of reactions has been performed close to room temperature. Synchronous reaction initiation has only been achieved with enzymes of relatively slow turnover (<0.1 s(-1)) and caged compounds of high quantum yield. Here X-ray crystallography and microspectrophotometry were used to provide evidence that (nitrophenyl)ethyl (NPE) ester bonds can be photolyzed by UV light at cryotemperatures. NPE-caged ATP in flash-cooled crystals of Mycobacterium tuberculosis thymidylate kinase was photolyzed successfully at 100-150 K as assessed by the structural observation of ATP-dependent enzymatic conversion of TMP to TDP after temporarily warming the crystals to room temperature. A new method is proposed in which cryo-photolysis combined with temperature-controlled protein crystallography can be used to trap reaction intermediates even in some of the fastest enzymes and/or when only compounds of low quantum yield are available. Raising the temperature after cryophotolysis may allow a transition barrier to be passed and an intermediate to accumulate in the crystal. A comparable method has only been used so far with proteins displaying endogenous photosensitivity. The approach described here opens the way to studying the reaction mechanisms of a much larger number of crystalline enzymes. Furthermore, it is shown that X-ray-induced radiolysis of caged compounds occurs if high-intensity synchrotron beamlines are used. This caveat should be taken into account when deriving data-collection protocols. It could also be used potentially as a way to trigger reactions.

Crystal structures of a template-independent DNA polymerase: murine terminal deoxynucleotidyltransferase

The crystal structure of the catalytic core of murine terminal deoxynucleotidyltransferase (TdT) at 2.35 A resolution reveals a typical DNA polymerase beta-like fold locked in a closed form. In addition, the structures of two different binary complexes, one with an oligonucleotide primer and the other with an incoming ddATP-Co(2+) complex, show that the substrates and the two divalent ions in the catalytic site are positioned in TdT in a manner similar to that described for the human DNA polymerase beta ternary complex, suggesting a common two metal ions mechanism of nucleotidyl transfer in these two proteins. The inability of TdT to accommodate a template strand can be explained by steric hindrance at the catalytic site caused by a long lariat-like loop, which is absent in DNA polymerase beta. However, displacement of this discriminating loop would be sufficient to unmask a number of evolutionarily conserved residues, which could then interact with a template DNA strand. The present structure can be used to model the recently discovered human polymerase mu, with which it shares 43% sequence identity.

The growth of tomato (Lycopersicon esculentum Mill.) hypocotyls in the light and in darkness differentially involves auxin

Light and auxin antagonistically regulate hypocotyl elongation. We have investigated the physiological interactions of light and auxin in the control of tomato (Lycopersicon esculentum Mill.) hypocotyl elongation by studying the auxin-insensitive mutant diageotropica (dgt). The length of the hypocotyls of the dgt mutant is significantly reduced when compared to the wild type line Ailsa Craig (AC) in the dark and under red light, but not under the other light conditions tested, indicating that auxin sensitivity is involved in the elongation of hypocotyls only in these conditions. Similarly, the auxin transport inhibitor naphthylphthalamic [correction of naphtylphtalamic] acid (NPA) differentially affects elongation of dark- or light-grown hypocotyls of the MoneyMaker (MM) tomato wild type. Using different photomorphogenic mutants, we demonstrate that at least phytochrome A, phytochrome B1 and, to a much lesser extent [correction of extend], cryptochrome 1, are necessary for a switch from an auxin transport-dependent elongation of hypocotyls in the dark to an auxin transport-independent elongation in the light. Interestingly, the dgt mutant and NPA-treated seedlings exhibit a looped phenotype only under red light, indicating that the negative gravitropism of hypocotyls also differentially involves auxin in the various light conditions.