MinActionPath2: path generation between different conformations of large macromolecular assemblies by action minimization.

Recent progress in solving macromolecular structures and assemblies by cryogenic electron microscopy techniques enables sampling of their conformations in different states that are relevant to their biological function. Knowing the transition path between these conformations would provide new avenues for drug discovery. While the experimental study of transition paths is intrinsically difficult, in-silico methods can be used to generate an initial guess for those paths. The Elastic Network Model (ENM), along with a coarse-grained representation (CG) of the structures are among the most popular models to explore such possible paths. Here we propose an update to our software platform MinActionPath that generates non-linear transition paths based on ENM and CG models, using action minimization to solve the equations of motion. The new website enables the study of large structures such as ribosomes or entire virus envelopes. It provides direct visualization of the trajectories along with quantitative analyses of their behaviors at http://dynstr.pasteur.fr/servers/minactionpath/minactionpath2_submission.

Reclassification of family A DNA polymerases reveals novel functional subfamilies and distinctive structural features.

Family A DNA polymerases (PolAs) form an important and well-studied class of extant polymerases participating in DNA replication and repair. Nonetheless, despite the characterization of multiple subfamilies in independent, dedicated works, their comprehensive classification thus far is missing. We therefore re-examine all presently available PolA sequences, converting their pairwise similarities into positions in Euclidean space, separating them into 19 major clusters. While 11 of them correspond to known subfamilies, eight had not been characterized before. For every group, we compile their general characteristics, examine their phylogenetic relationships and perform conservation analysis in the essential sequence motifs. While most subfamilies are linked to a particular domain of life (including phages), one subfamily appears in Bacteria, Archaea and Eukaryota. We also show that two new bacterial subfamilies contain functional enzymes. We use AlphaFold2 to generate high-confidence prediction models for all clusters lacking an experimentally determined structure. We identify new, conserved features involving structural alterations, ordered insertions and an apparent structural incorporation of a uracil-DNA glycosylase (UDG) domain. Finally, genetic and structural analyses of a subset of T7-like phages indicate a splitting of the 3'-5' exo and pol domains into two separate genes, observed in PolAs for the first time.

Structural dynamics and determinants of 2-aminoadenine specificity in DNA polymerase DpoZ of vibriophage ϕVC8.

All genetic information in cellular life is stored in DNA copolymers composed of four basic building blocks (ATGC-DNA). In contrast, a group of bacteriophages belonging to families Siphoviridae and Podoviridae has abandoned the usage of one of them, adenine (A), replacing it with 2-aminoadenine (Z). The resulting ZTGC-DNA is more stable than its ATGC-DNA counterpart, owing to the additional hydrogen bond present in the 2-aminoadenine:thymine (Z:T) base pair, while the additional amino group also confers resistance to the host endonucleases. Recently, two classes of replicative proteins found in ZTGC-DNA-containing phages were characterized and one of them, DpoZ from DNA polymerase A (PolA) family, was shown to possess significant Z-vs-A specificity. Here, we present the crystallographic structure of the apo form of DpoZ of vibriophage ϕVC8, composed of the 3'-5' exonuclease and polymerase domains. We captured the enzyme in two conformations that involve the tip of the thumb subdomain and the exonuclease domain. We highlight insertions and mutations characteristic of ϕVC8 DpoZ and its close homologues. Through mutagenesis and functional assays we suggest that the preference of ϕVC8 DpoZ towards Z relies on a polymerase backtracking process, more efficient when the nascent base pair is A:T than when it is Z:T.

Structural basis for allosteric transitions of a multidomain pentameric ligand-gated ion channel.

Pentameric ligand-gated ion channels (pLGICs) are allosteric receptors that mediate rapid electrochemical signal transduction in the animal nervous system through the opening of an ion pore upon binding of neurotransmitters. Orthologs have been found and characterized in prokaryotes and they display highly similar structure-function relationships to eukaryotic pLGICs; however, they often encode greater architectural diversity involving additional amino-terminal domains (NTDs). Here we report structural, functional, and normal-mode analysis of two conformational states of a multidomain pLGIC, called DeCLIC, from a Desulfofustis deltaproteobacterium, including a periplasmic NTD fused to the conventional ligand-binding domain (LBD). X-ray structure determination revealed an NTD consisting of two jelly-roll domains interacting across each subunit interface. Binding of Ca2+ at the LBD subunit interface was associated with a closed transmembrane pore, with resolved monovalent cations intracellular to the hydrophobic gate. Accordingly, DeCLIC-injected oocytes conducted currents only upon depletion of extracellular Ca2+; these were insensitive to quaternary ammonium block. Furthermore, DeCLIC crystallized in the absence of Ca2+ with a wide-open pore and remodeled periplasmic domains, including increased contacts between the NTD and classic LBD agonist-binding sites. Functional, structural, and dynamical properties of DeCLIC paralleled those of sTeLIC, a pLGIC from another symbiotic prokaryote. Based on these DeCLIC structures, we would reclassify the previous structure of bacterial ELIC (the first high-resolution structure of a pLGIC) as a "locally closed" conformation. Taken together, structures of DeCLIC in multiple conformations illustrate dramatic conformational state transitions and diverse regulatory mechanisms available to ion channels in pLGICs, particularly involving Ca2+ modulation and periplasmic NTDs.

Structure of the DP1-DP2 PolD complex bound with DNA and its implications for the evolutionary history of DNA and RNA polymerases.

PolD is an archaeal replicative DNA polymerase (DNAP) made of a proofreading exonuclease subunit (DP1) and a larger polymerase catalytic subunit (DP2). Recently, we reported the individual crystal structures of the DP1 and DP2 catalytic cores, thereby revealing that PolD is an atypical DNAP that has all functional properties of a replicative DNAP but with the catalytic core of an RNA polymerase (RNAP). We now report the DNA-bound cryo-electron microscopy (cryo-EM) structure of the heterodimeric DP1-DP2 PolD complex from Pyrococcus abyssi, revealing a unique DNA-binding site. Comparison of PolD and RNAPs extends their structural similarities and brings to light the minimal catalytic core shared by all cellular transcriptases. Finally, elucidating the structure of the PolD DP1-DP2 interface, which is conserved in all eukaryotic replicative DNAPs, clarifies their evolutionary relationships with PolD and sheds light on the domain acquisition and exchange mechanism that occurred during the evolution of the eukaryotic replisome.

Parameterizing elastic network models to capture the dynamics of proteins.

Coarse-grained normal mode analyses of protein dynamics rely on the idea that the geometry of a protein structure contains enough information for computing its fluctuations around its equilibrium conformation. This geometry is captured in the form of an elastic network (EN), namely a network of edges between its residues. The normal modes of a protein are then identified with the normal modes of its EN. Different approaches have been proposed to construct ENs, focusing on the choice of the edges that they are comprised of, and on their parameterizations by the force constants associated with those edges. Here we propose new tools to guide choices on these two facets of EN. We study first different geometric models for ENs. We compare cutoff-based ENs, whose edges have lengths that are smaller than a cutoff distance, with Delaunay-based ENs and find that the latter provide better representations of the geometry of protein structures. We then derive an analytical method for the parameterization of the EN such that its dynamics leads to atomic fluctuations that agree with experimental B-factors. To limit overfitting, we attach a parameter referred to as flexibility constant to each atom instead of to each edge in the EN. The parameterization is expressed as a non-linear optimization problem whose parameters describe both rigid-body and internal motions. We show that this parameterization leads to improved ENs, whose dynamics mimic MD simulations better than ENs with uniform force constants, and reduces the number of normal modes needed to reproduce functional conformational changes.

Extracting Dynamical Correlations and Identifying Key Residues for Allosteric Communication in Proteins by correlationplus.

Extracting dynamical pairwise correlations and identifying key residues from large molecular dynamics trajectories or normal-mode analysis of coarse-grained models are important for explaining various processes like ligand binding, mutational effects, and long-distance interactions. Efficient and flexible tools to perform this task can provide new insights about residues involved in allosteric regulation and protein function. In addition, combining and comparing dynamical coupling information with sequence coevolution data can help to understand better protein function. To this aim, we developed a Python package called correlationplus to calculate, visualize, and analyze pairwise correlations. In this way, the package aids to identify key residues and interactions in proteins. The source code of correlationplus is available under LGPL version 3 at https://github.com/tekpinar/correlationplus. The current version of the package (0.2.0) can be installed with common installation methods like conda or pip in addition to source code installation. Moreover, docker images are also available for usage of the code without installation.

Crystal structures of a pentameric ion channel gated by alkaline pH show a widely open pore and identify a cavity for modulation.

Pentameric ligand-gated ion channels (pLGICs) constitute a widespread class of ion channels, present in archaea, bacteria, and eukaryotes. Upon binding of their agonists in the extracellular domain, the transmembrane pore opens, allowing ions to go through, via a gating mechanism that can be modulated by a number of drugs. Even though high-resolution structural information on pLGICs has increased in a spectacular way in recent years, both in bacterial and in eukaryotic systems, the structure of the open channel conformation of some intensively studied receptors whose structures are known in a nonactive (closed) form, such as Erwinia chrysanthemi pLGIC (ELIC), is still lacking. Here we describe a gammaproteobacterial pLGIC from an endo-symbiont of Tevnia jerichonana (sTeLIC), whose sequence is closely related to the pLGIC from ELIC with 28% identity. We provide an X-ray crystallographic structure at 2.3 Å in an active conformation, where the pore is found to be more open than any current conformation found for pLGICs. In addition, two charged restriction rings are present in the vestibule. Functional characterization shows sTeLIC to be a cationic channel activated at alkaline pH. It is inhibited by divalent cations, but not by quaternary ammonium ions, such as tetramethylammonium. Additionally, we found that sTeLIC is allosterically potentiated by aromatic amino acids Phe and Trp, as well as their derivatives, such as 4-bromo-cinnamate, whose cocrystal structure reveals a vestibular binding site equivalent to, but more deeply buried than, the one already described for benzodiazepines in ELIC.

Electrostatics, proton sensor, and networks governing the gating transition in GLIC, a proton-gated pentameric ion channel.

The pentameric ligand-gated ion channel (pLGIC) from Gloeobacter violaceus (GLIC) has provided insightful structure-function views on the permeation process and the allosteric regulation of the pLGICs family. However, GLIC is activated by pH instead of a neurotransmitter and a clear picture for the gating transition driven by protons is still lacking. We used an electrostatics-based (finite difference Poisson-Boltzmann/Debye-Hückel) method to predict the acidities of all aspartic and glutamic residues in GLIC, both in its active and closed-channel states. Those residues with a predicted pKa close to the experimental pH50 were individually replaced by alanine and the resulting variant receptors were titrated by ATR/FTIR spectroscopy. E35, located in front of loop F far away from the orthosteric site, appears as the key proton sensor with a measured individual pKa at 5.8. In the GLIC open conformation, E35 is connected through a water-mediated hydrogen-bond network first to the highly conserved electrostatic triad R192-D122-D32 and then to Y197-Y119-K248, both located at the extracellular domain-transmembrane domain interface. The second triad controls a cluster of hydrophobic side chains from the M2-M3 loop that is remodeled during the gating transition. We solved 12 crystal structures of GLIC mutants, 6 of them being trapped in an agonist-bound but nonconductive conformation. Combined with previous data, this reveals two branches of a continuous network originating from E35 that reach, independently, the middle transmembrane region of two adjacent subunits. We conclude that GLIC's gating proceeds by making use of loop F, already known as an allosteric site in other pLGICs, instead of the classic orthosteric site.

Structural evidence for an in trans base selection mechanism involving Loop1 in polymerase µ at an NHEJ double-strand break junction.

Eukaryotic DNA polymerase (Pol) X family members such as Pol µ and terminal deoxynucleotidyl transferase (TdT) are important components for the nonhomologous DNA end-joining (NHEJ) pathway. TdT participates in a specialized version of NHEJ, V(D)J recombination. It has primarily nontemplated polymerase activity but can take instructions across strands from the downstream dsDNA, and both activities are highly dependent on a structural element called Loop1. However, it is unclear whether Pol µ follows the same mechanism, because the structure of its Loop1 is disordered in available structures. Here, we used a chimeric TdT harboring Loop1 of Pol µ that recapitulated the functional properties of Pol µ in ligation experiments. We solved three crystal structures of this TdT chimera bound to several DNA substrates at 1.96-2.55 Å resolutions, including a full DNA double-strand break (DSB) synapsis. We then modeled the full Pol µ sequence in the context of one these complexes. The atomic structure of an NHEJ junction with a Pol X construct that mimics Pol µ in a reconstituted system explained the distinctive properties of Pol µ compared with TdT. The structure suggested a mechanism of base selection relying on Loop1 and taking instructions via the in trans templating base independently of the primer strand. We conclude that our atomic-level structural observations represent a paradigm shift for the mechanism of base selection in the Pol X family of DNA polymerases.

Full mutational mapping of titratable residues helps to identify proton-sensors involved in the control of channel gating in the Gloeobacter violaceus pentameric ligand-gated ion channel.

The Gloeobacter violaceus ligand-gated ion channel (GLIC) has been extensively studied by X-ray crystallography and other biophysical techniques. This provided key insights into the general gating mechanism of pentameric ligand-gated ion channel (pLGIC) signal transduction. However, the GLIC is activated by lowering the pH and the location of its putative proton activation site(s) still remain(s) unknown. To this end, every Asp, Glu, and His residue was mutated individually or in combination and investigated by electrophysiology. In addition to the mutational analysis, key mutations were structurally resolved to address whether particular residues contribute to proton sensing, or alternatively to GLIC-gating, independently of the side chain protonation. The data show that multiple residues located below the orthosteric site, notably E26, D32, E35, and D122 in the lower part of the extracellular domain (ECD), along with E222, H235, E243, and H277 in the transmembrane domain (TMD), alter GLIC activation. D122 and H235 were found to also alter GLIC expression. E35 is identified as a key proton-sensing residue, whereby neutralization of its side chain carboxylate stabilizes the active state. Thus, proton activation occurs allosterically to the orthosteric site, at the level of multiple loci with a key contribution of the coupling interface between the ECD and TMD.

Rapid enzymatic synthesis of long RNA polymers: A simple protocol to generate RNA libraries with random sequences.

RNA aptamers have several advantages over DNA aptamers due to their propensity to fold into three-dimensional structures. However, the synthesis of large RNA libraries remains a challenge as it requires more precautions to conserve their functional integrity, especially when such libraries are intended for aptamers or ribozymes selection. Here, we present an enzymatic method that enables the rapid synthesis of RNA polymers thanks to the efficient incorporation of ribonucleotides (NTPs) as well as chemically modified ribonucleotides by human DNA polymerase Theta (θ) mutants. These mutants have the ability to generate long single-stranded RNA polynucleotides of random sequences due to their improved template-free terminal nucleotidyltransferase activity. Here we describe the detailed protocols to produce large and diverse libraries of RNA, to make them ready to use in repeated cycles of Systematic Evolution of Ligands by Exponential enrichment (SELEX) and to synthesize C2'-modified nucleic acids with higher nuclease resistance.

Enzymatic synthesis of random sequences of RNA and RNA analogues by DNA polymerase theta mutants for the generation of aptamer libraries.

Nucleic acid aptamers, especially RNA, exhibit valuable advantages compared to protein therapeutics in terms of size, affinity and specificity. However, the synthesis of libraries of large random RNAs is still difficult and expensive. The engineering of polymerases able to directly generate these libraries has the potential to replace the chemical synthesis approach. Here, we start with a DNA polymerase that already displays a significant template-free nucleotidyltransferase activity, human DNA polymerase theta, and we mutate it based on the knowledge of its three-dimensional structure as well as previous mutational studies on members of the same polA family. One mutant exhibited a high tolerance towards ribonucleotides (NTPs) and displayed an efficient ribonucleotidyltransferase activity that resulted in the assembly of long RNA polymers. HPLC analysis and RNA sequencing of the products were used to quantify the incorporation of the four NTPs as a function of initial NTP concentrations and established the randomness of each generated nucleic acid sequence. The same mutant revealed a propensity to accept other modified nucleotides and to extend them in long fragments. Hence, this mutant can deliver random natural and modified RNA polymers libraries ready to use for SELEX, with custom lengths and balanced or unbalanced ratios.

String method solution of the gating pathways for a pentameric ligand-gated ion channel.

Pentameric ligand-gated ion channels control synaptic neurotransmission by converting chemical signals into electrical signals. Agonist binding leads to rapid signal transduction via an allosteric mechanism, where global protein conformational changes open a pore across the nerve cell membrane. We use all-atom molecular dynamics with a swarm-based string method to solve for the minimum free-energy gating pathways of the proton-activated bacterial GLIC channel. We describe stable wetted/open and dewetted/closed states, and uncover conformational changes in the agonist-binding extracellular domain, ion-conducting transmembrane domain, and gating interface that control communication between these domains. Transition analysis is used to compute free-energy surfaces that suggest allosteric pathways; stabilization with pH; and intermediates, including states that facilitate channel closing in the presence of an agonist. We describe a switching mechanism that senses proton binding by marked reorganization of subunit interface, altering the packing of ß-sheets to induce changes that lead to asynchronous pore-lining M2 helix movements. These results provide molecular details of GLIC gating and insight into the allosteric mechanisms for the superfamily of pentameric ligand-gated channels.

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.

An updated structural classification of replicative DNA polymerases.

Replicative DNA polymerases are nano-machines essential to life, which have evolved the ability to copy the genome with high fidelity and high processivity. In contrast with cellular transcriptases and ribosome machines, which evolved by accretion of complexity from a conserved catalytic core, no replicative DNA polymerase is universally conserved. Strikingly, four different families of DNA polymerases have evolved to perform DNA replication in the three domains of life. In Bacteria, the genome is replicated by DNA polymerases belonging to the A- and C-families. In Eukarya, genomic DNA is copied mainly by three distinct replicative DNA polymerases, Polα, Polδ, and Polε, which all belong to the B-family. Matters are more complicated in Archaea, which contain an unusual D-family DNA polymerase (PolD) in addition to PolB, a B-family replicative DNA polymerase that is homologous to the eukaryotic ones. PolD is a heterodimeric DNA polymerase present in all Archaea discovered so far, except Crenarchaea. While PolD is an essential replicative DNA polymerase, it is often underrepresented in the literature when the diversity of DNA polymerases is discussed. Recent structural studies have shown that the structures of both polymerase and proofreading active sites of PolD differ from other structurally characterized DNA polymerases, thereby extending the repertoire of folds known to perform DNA replication. This review aims to provide an updated structural classification of all replicative DNAPs and discuss their evolutionary relationships, both regarding the DNA polymerase and proofreading active sites.

Statistical Physics Approach to the Optimal Transport Problem.

Originally defined for the optimal allocation of resources, optimal transport (OT) has found many theoretical and practical applications in multiple domains of science and physics. In this Letter we develop a new method for solving the discrete version of this problem using techniques derived from statistical physics. We derive a strongly concave free energy function that captures the constraints of the OT problem at a finite temperature. Its maximum defines an optimal transport plan, or registration between the two discrete probability measures that are compared, as well as a pseudodistance between those measures that satisfies the triangular inequalities. The computation of this pseudodistance is fast and numerically stable. The temperature dependent OT pseudodistance is shown to decrease monotonically with respect to the inverse of the temperature and to converge to the standard OT distance at zero temperature, providing a robust framework for temperature annealing. We illustrate applications of this framework to the problem of image comparison.

Meet-U: Educating through research immersion.

We present a new educational initiative called Meet-U that aims to train students for collaborative work in computational biology and to bridge the gap between education and research. Meet-U mimics the setup of collaborative research projects and takes advantage of the most popular tools for collaborative work and of cloud computing. Students are grouped in teams of 4-5 people and have to realize a project from A to Z that answers a challenging question in biology. Meet-U promotes "coopetition," as the students collaborate within and across the teams and are also in competition with each other to develop the best final product. Meet-U fosters interactions between different actors of education and research through the organization of a meeting day, open to everyone, where the students present their work to a jury of researchers and jury members give research seminars. This very unique combination of education and research is strongly motivating for the students and provides a formidable opportunity for a scientific community to unite and increase its visibility. We report on our experience with Meet-U in two French universities with master's students in bioinformatics and modeling, with protein-protein docking as the subject of the course. Meet-U is easy to implement and can be straightforwardly transferred to other fields and/or universities. All the information and data are available at www.meet-u.org.

Barbiturates Bind in the GLIC Ion Channel Pore and Cause Inhibition by Stabilizing a Closed State.

Barbiturates induce anesthesia by modulating the activity of anionic and cationic pentameric ligand-gated ion channels (pLGICs). Despite more than a century of use in clinical practice, the prototypic binding site for this class of drugs within pLGICs is yet to be described. In this study, we present the first X-ray structures of barbiturates bound to GLIC, a cationic prokaryotic pLGIC with excellent structural homology to other relevant channels sensitive to general anesthetics and, as shown here, to barbiturates, at clinically relevant concentrations. Several derivatives of barbiturates containing anomalous scatterers were synthesized, and these derivatives helped us unambiguously identify a unique barbiturate binding site within the central ion channel pore in a closed conformation. In addition, docking calculations around the observed binding site for all three states of the receptor, including a model of the desensitized state, showed that barbiturates preferentially stabilize the closed state. The identification of this pore binding site sheds light on the mechanism of barbiturate inhibition of cationic pLGICs and allows the rationalization of several structural and functional features previously observed for barbiturates.

Allosteric and hyperekplexic mutant phenotypes investigated on an α1 glycine receptor transmembrane structure.

The glycine receptor (GlyR) is a pentameric ligand-gated ion channel (pLGIC) mediating inhibitory transmission in the nervous system. Its transmembrane domain (TMD) is the target of allosteric modulators such as general anesthetics and ethanol and is a major locus for hyperekplexic congenital mutations altering the allosteric transitions of activation or desensitization. We previously showed that the TMD of the human α1GlyR could be fused to the extracellular domain of GLIC, a bacterial pLGIC, to form a functional chimera called Lily. Here, we overexpress Lily in Schneider 2 insect cells and solve its structure by X-ray crystallography at 3.5 Å resolution. The TMD of the α1GlyR adopts a closed-channel conformation involving a single ring of hydrophobic residues at the center of the pore. Electrophysiological recordings show that the phenotypes of key allosteric mutations of the α1GlyR, scattered all along the pore, are qualitatively preserved in this chimera, including those that confer decreased sensitivity to agonists, constitutive activity, decreased activation kinetics, or increased desensitization kinetics. Combined structural and functional data indicate a pore-opening mechanism for the α1GlyR, suggesting a structural explanation for the effect of some key hyperekplexic allosteric mutations. The first X-ray structure of the TMD of the α1GlyR solved here using GLIC as a scaffold paves the way for mechanistic investigation and design of allosteric modulators of a human receptor.