Molecular basis for proofreading by the unique exonuclease domain of Family-D DNA polymerases

Replicative DNA polymerases duplicate entire genomes at high fidelity. This feature is shared among the three domains of life and is facilitated by their dual polymerase and exonuclease activities. Family D replicative DNA polymerases (PolD), found exclusively in Archaea, contain an unusual RNA polymerase-like catalytic core, and a unique Mre11-like proofreading active site. Here, we present cryo-EM structures of PolD trapped in a proofreading mode, revealing an unanticipated correction mechanism that extends the repertoire of protein domains known to be involved in DNA proofreading. Based on our experimental structures, mutants of PolD were designed and their contribution to mismatch bypass and exonuclease kinetics was determined. This study sheds light on the convergent evolution of structurally distinct families of DNA polymerases, and the domain acquisition and exchange mechanism that occurred during the evolution of the replisome in the three domains of life.

Analyzing the Geometry and Dynamics of Viral Structures: A Review of Computational Approaches Based on Alpha Shape Theory, Normal Mode Analysis, and Poisson-Boltzmann Theories

The current SARS-CoV-2 pandemic highlights our fragility when we are exposed to emergent viruses either directly or through zoonotic diseases. Fortunately, our knowledge of the biology of those viruses is improving. In particular, we have more and more structural information on virions, i.e., the infective form of a virus that includes its genomic material and surrounding protective capsid, and on their gene products. It is important to have methods that enable the analyses of structural information on such large macromolecular systems. We review some of those methods in this paper. We focus on understanding the geometry of virions and viral structural proteins, their dynamics, and their energetics, with the ambition that this understanding can help design antiviral agents. We discuss those methods in light of the specificities of those structures, mainly that they are huge. We focus on three of our own methods based on the alpha shape theory for computing geometry, normal mode analyses to study dynamics, and modified Poisson-Boltzmann theories to study the organization of ions and co-solvent and solvent molecules around biomacromolecules. The corresponding software has computing times that are compatible with the use of regular desktop computers. We show examples of their applications on some outer shells and structural proteins of the West Nile Virus.

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.

DNA-binding mechanism and evolution of replication protein A

Replication Protein A (RPA) is a heterotrimeric single stranded DNA-binding protein with essential roles in DNA replication, recombination and repair. Little is known about the structure of RPA in Archaea, the third domain of life. By using an integrative structural, biochemical and biophysical approach, we extensively characterize RPA from Pyrococcus abyssi in the presence and absence of DNA. The obtained X-ray and cryo-EM structures reveal that the trimerization core and interactions promoting RPA clustering on ssDNA are shared between archaea and eukaryotes. However, we also identified a helical domain named AROD (Acidic Rpa1 OB-binding Domain), and showed that, in Archaea, RPA forms an unanticipated tetrameric supercomplex in the absence of DNA. The four RPA molecules clustered within the tetramer could efficiently coat and protect stretches of ssDNA created by the advancing replisome. Finally, our results provide insights into the evolution of this primordial replication factor in eukaryotes.

Brain charts for the human lifespan

Over the past few decades, neuroimaging has become a ubiquitous tool in basic research and clinical studies of the human brain. However, no reference standards currently exist to quantify individual differences in neuroimaging metrics over time, in contrast to growth charts for anthropometric traits such as height and weight1. Here we assemble an interactive open resource to benchmark brain morphology derived from any current or future sample of MRI data ( http://www.brainchart.io/ ). With the goal of basing these reference charts on the largest and most inclusive dataset available, acknowledging limitations due to known biases of MRI studies relative to the diversity of the global population, we aggregated 123,984 MRI scans, across more than 100 primary studies, from 101,457 human participants between 115 days post-conception to 100 years of age. MRI metrics were quantified by centile scores, relative to non-linear trajectories2 of brain structural changes, and rates of change, over the lifespan. Brain charts identified previously unreported neurodevelopmental milestones3, showed high stability of individuals across longitudinal assessments, and demonstrated robustness to technical and methodological differences between primary studies. Centile scores showed increased heritability compared with non-centiled MRI phenotypes, and provided a standardized measure of atypical brain structure that revealed patterns of neuroanatomical variation across neurological and psychiatric disorders. In summary, brain charts are an essential step towards robust quantification of individual variation benchmarked to normative trajectories in multiple, commonly used neuroimaging phenotypes.

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.

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.

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.

Characterization of a triad of genes in cyanophage S-2L sufficient to replace adenine by 2-aminoadenine in bacterial DNA

Cyanophage S-2L is known to profoundly alter the biophysical properties of its DNA by replacing all adenines (A) with 2-aminoadenines (Z), which still pair with thymines but with a triple hydrogen bond. It was recently demonstrated that a homologue of adenylosuccinate synthetase (PurZ) and a dATP triphosphohydrolase (DatZ) are two important pieces of the metabolism of 2-aminoadenine, participating in the synthesis of ZTGC-DNA. Here, we determine that S-2L PurZ can use either dATP or ATP as a source of energy, thereby also depleting the pool of nucleotides in dATP. Furthermore, we identify a conserved gene (mazZ) located between purZ and datZ genes in S-2L and related phage genomes. We show that it encodes a (d)GTP-specific diphosphohydrolase, thereby providing the substrate of PurZ in the 2-aminoadenine synthesis pathway. High-resolution crystal structures of S-2L PurZ and MazZ with their respective substrates provide a rationale for their specificities. The Z-cluster made of these three genes - datZ, mazZ and purZ - was expressed in E. coli, resulting in a successful incorporation of 2-aminoadenine in the bacterial chromosomal and plasmidic DNA. This work opens the possibility to study synthetic organisms containing ZTGC-DNA.

PI3K functions as a hub in mechanotransduction

Mammalian cells integrate different types of stimuli that govern their fate. These stimuli encompass biochemical as well as biomechanical cues (shear, tensile, and compressive stresses) that are usually studied separately. The phosphatidylinositol 3-kinase (PI3K) enzymes, producing signaling phosphoinositides at plasma and intracellular membranes, are key in intracellular signaling and vesicular trafficking pathways. Recent evidence in cancer research demonstrates that these enzymes are essential in mechanotransduction. Despite this, the importance of the integration of biomechanical cues and PI3K-driven biochemical signals is underestimated. In this opinion article, we make the hypothesis that modeling of biomechanical cues is critical to understand PI3K oncogenicity. We also identify known/missing knowledge in terms of isoform specificity and molecular pathways of activation, knowledge that is needed for clinical applications.

Simultaneous Identification of Multiple Binding Sites in Proteins: A Statistical Mechanics Approach

We present an extension of the Poisson-Boltzmann model in which the solute of interest is immersed in an assembly of self-orienting Langevin water dipoles, anions, cations, and hydrophobic molecules, all of variable densities. Interactions between charges are controlled by electrostatics, while hydrophobic interactions are modeled with a Yukawa potential. We impose steric constraints by assuming that the system is represented on a cubic lattice. We also assume incompressibility; i.e., all sites of the lattice are occupied. This model, which we refer to as the Hydrophobic Dipolar Poisson-Boltzmann Langevin (HDPBL) model, leads to a system of two equations whose solutions give the water dipole, salt, and hydrophobic molecule densities, all of them in the presence of the others in a self-consistent way. We use those to study the organization of the ions, cosolvent, and solvent molecules around proteins. In particular, peaks of densities are expected to reveal, simultaneously, the presence of compatible binding sites of different kinds on a protein. We have tested and validated the ability of HDPBL to detect pockets in proteins that bind to hydrophobic ligands, polar ligands, and charged small probes as well as to characterize the binding sites of lipids for membrane proteins.

Fast and efficient purification of SARS-CoV-2 RNA dependent RNA polymerase complex expressed in Escherichia coli

To stop the COVID-19 pandemic due to the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), which caused more than 2.5 million deaths to date, new antiviral molecules are urgently needed. The replication of SARS-CoV-2 requires the RNA-dependent RNA polymerase (RdRp), making RdRp an excellent target for antiviral agents. RdRp is a multi-subunit complex composed of 3 viral proteins named nsp7, nsp8 and nsp12 that ensure the ~30 kb RNA genome’s transcription and replication. The main strategies employed so far for the overproduction of RdRp consist of expressing and purifying the three subunits separately before assembling the complex in vitro. However, nsp12 shows limited solubility in bacterial expression systems and is often produced in insect cells. Here, we describe an alternative strategy to co-express the full SARS-CoV-2 RdRp in E. coli, using a single plasmid. Characterization of the purified recombinant SARS-CoV-2 RdRp shows that it forms a complex with the expected (nsp7)(nsp8)₂(nsp12) stoichiometry. RNA polymerization activity was measured using primer-extension assays showing that the purified enzyme is functional. The purification protocol can be achieved in one single day, surpassing in speed all other published protocols. Our construct is ideally suited for screening RdRp and its variants against very large chemical compounds libraries and has been made available to the scientific community through the Addgene plasmid depository (Addgene ID: 165451).

Physics approach to the variable-mass optimal-transport problem

Optimal transport (OT) has become a discipline by itself that offers solutions to a wide range of theoretical problems in probability and mathematics with applications in several applied fields such as imaging sciences, machine learning, and in data sciences in general. The traditional OT problem suffers from a severe limitation: its balance condition imposes that the two distributions to be compared be normalized and have the same total mass. However, it is important for many applications to be able to relax this constraint and allow for mass creation and/or destruction. This is true, for example, in all problems requiring partial matching. In this paper, we propose an approach to solving a generalized version of the OT problem, which we refer to as the discrete variable-mass optimal-transport (VMOT) problem, using techniques adapted from statistical physics. Our first contribution is to fully describe this formalism, including all the proofs of its main claims. In particular, we derive a strongly concave effective free-energy function that captures the constraints of the VMOT problem at a finite temperature. From its maximum we derive a weak distance (i.e., a divergence) between possibly unbalanced distribution functions. The temperature-dependent OT distance decreases monotonically to the standard variable-mass OT distance, providing a robust framework for temperature annealing. Our second contribution is to show that the implementation of this formalism has the same properties as the regularized OT algorithms in time complexity, making it a competitive approach to solving the VMOT problem. We illustrate applications of the framework to the problem of partial two- and three-dimensional shape-matching problems.

Structural evidence for the binding of monocarboxylates and dicarboxylates at pharmacologically relevant extracellular sites of a pentameric ligand-gated ion channel

GLIC is a bacterial homologue of the pentameric ligand-gated ion channels (pLGICs) that mediate the fast chemical neurotransmission of nerve signalling in eukaryotes. Because the activation and allosteric modulation features are conserved among prokaryotic and eukaryotic pLGICs, GLIC is commonly used as a model to study the allosteric transition and structural pharmacology of pLGICs. It has previously been shown that GLIC is inhibited by some carboxylic acid derivatives. Here, experimental evidence for carboxylate binding to GLIC is provided by solving its X-ray structures with a series of monocarboxylate and dicarboxylate derivatives, and two carboxylate-binding sites are described: (i) the `intersubunit' site that partially overlaps the canonical pLGIC orthosteric site and (ii) the `intrasubunit' vestibular site, which is only occupied by a subset of the described derivatives. While the intersubunit site is widely conserved in all pLGICs, the intrasubunit site is only conserved in cationic eukaryotic pLGICs. This study sheds light on the importance of these two extracellular modulation sites as potential drug targets in pLGICs.

Structural basis for the increased processivity of D-family DNA polymerases in complex with PCNA

Replicative DNA polymerases (DNAPs) have evolved the ability to copy the genome with high processivity and fidelity. In Eukarya and Archaea, the processivity of replicative DNAPs is greatly enhanced by its binding to the proliferative cell nuclear antigen (PCNA) that encircles the DNA. We determined the cryo-EM structure of the DNA-bound PolD–PCNA complex from Pyrococcus abyssi at 3.77 Å. Using an integrative structural biology approach — combining cryo-EM, X-ray crystallography, protein–protein interaction measurements, and activity assays — we describe the molecular basis for the interaction and cooperativity between a replicative DNAP and PCNA. PolD recruits PCNA via a complex mechanism, which requires two different PIP-boxes. We infer that the second PIP-box, which is shared with the eukaryotic Polα replicative DNAP, plays a dual role in binding either PCNA or primase, and could be a master switch between an initiation and a processive phase during replication.

Replicative DNA polymerases (DNAPs) have evolved the ability to copy the genome with high processivity and fidelity. Here, the authors present a cryo-EM structure of the DNA-bound PolD–PCNA complex from Pyrococcus abyssi to reveal the molecular basis for the interaction and cooperativity between a replicative DNAP and PCNA.

Structural Basis for a Bimodal Allosteric Mechanism of General Anesthetic Modulation in Pentameric Ligand-Gated Ion Channels

Ion channel modulation by general anesthetics is a vital pharmacological process with implications for receptor biophysics and drug development. Functional studies have implicated conserved sites of both potentiation and inhibition in pentameric ligand-gated ion channels, but a detailed structural mechanism for these bimodal effects is lacking. The prokaryotic model protein GLIC recapitulates anesthetic modulation of human ion channels, and it is accessible to structure determination in both apparent open and closed states. Here, we report ten X-ray structures and electrophysiological characterization of GLIC variants in the presence and absence of general anesthetics, including the surgical agent propofol. We show that general anesthetics can allosterically favor closed channels by binding in the pore or favor open channels via various subsites in the transmembrane domain. Our results support an integrated, multi-site mechanism for allosteric modulation, and they provide atomic details of both potentiation and inhibition by one of the most common general anesthetics.

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.

Optimal transport at finite temperature

Optimal transport (OT) has become a discipline by itself that offers solutions to a wide range of theoretical problems in probability and mathematics. Despite its appealing theoretical properties, solving the OT problem involves the resolution of a linear program whose computational cost can quickly become prohibitive whenever the size of the problem exceeds a few hundred points. The recent introduction of entropy regularization, however, has led to the development of fast algorithms for solving an approximate OT problem. The successes of those algorithms have resulted in a popularization of the applications of OT in several applied fields such as imaging sciences and machine learning, and in data sciences in general. Problems remain, however, as to the numerical convergence of those regularized approximations towards the actual OT solution. In addition, the physical meaning of this regularization is unclear. In this paper, we propose an approach to solving the discrete OT problem using techniques adapted from statistical physics. Our first contribution is to fully describe this formalism, including all the proofs of its main claims. In particular we derive a strongly concave effective free energy function that captures the constraints of the optimal transport problem at a finite temperature. Its maximum defines a pseudo distance between the two set of weighted points that are compared, which satisfies the triangular inequalities. The temperature dependent OT pseudo distance decreases monotonically to the standard OT distance, providing a robust framework for temperature annealing. Our second contribution is to show that the implementation of this formalism has the same properties as the regularized OT algorithms in time complexity, making it a competitive approach to solving the OT problem. We illustrate applications of the framework to the problem of protein fold recognition based on sequence information only.

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.

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.

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.

Combined approaches from physics, statistics, and computer science for ab initio protein structure prediction: ex unitate vires (unity is strength)?

Connecting the dots among the amino acid sequence of a protein, its structure, and its function remains a central theme in molecular biology, as it would have many applications in the treatment of illnesses related to misfolding or protein instability. As a result of high-throughput sequencing methods, biologists currently live in a protein sequence-rich world. However, our knowledge of protein structure based on experimental data remains comparatively limited. As a consequence, protein structure prediction has established itself as a very active field of research to fill in this gap. This field, once thought to be reserved for theoretical biophysicists, is constantly reinventing itself, borrowing ideas informed by an ever-increasing assembly of scientific domains, from biology, chemistry, (statistical) physics, mathematics, computer science, statistics, bioinformatics, and more recently data sciences. We review the recent progress arising from this integration of knowledge, from the development of specific computer architecture to allow for longer timescales in physics-based simulations of protein folding to the recent advances in predicting contacts in proteins based on detection of coevolution using very large data sets of aligned protein sequences.

mTORC1 Controls Phase Separation and the Biophysical Properties of the Cytoplasm by Tuning Crowding

Macromolecular crowding has a profound impact on reaction rates and the physical properties of the cell interior, but the mechanisms that regulate crowding are poorly understood. We developed genetically encoded multimeric nanoparticles (GEMs) to dissect these mechanisms. GEMs are homomultimeric scaffolds fused to a fluorescent protein that self-assemble into bright, stable particles of defined size and shape. By combining tracking of GEMs with genetic and pharmacological approaches, we discovered that the mTORC1 pathway can modulate the effective diffusion coefficient of particles ≥20 nm in diameter more than 2-fold by tuning ribosome concentration, without any discernable effect on the motion of molecules ≤5 nm. This change in ribosome concentration affected phase separation both in vitro and in vivo. Together, these results establish a role for mTORC1 in controlling both the mesoscale biophysical properties of the cytoplasm and biomolecular condensation.

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.

mTORC1 Controls Phase Separation and the Biophysical Properties of the Cytoplasm by Tuning Crowding

Macromolecular crowding has a profound impact on reaction rates and the physical properties of the cell interior, but the mechanisms that regulate crowding are poorly understood. We developed genetically encoded multimeric nanoparticles (GEMs) to dissect these mechanisms. GEMs are homomultimeric scaffolds fused to a fluorescent protein that self-assemble into bright, stable particles of defined size and shape. By combining tracking of GEMs with genetic and pharmacological approaches, we discovered that the mTORC1 pathway can modulate the effective diffusion coefficient of particles ≥20 nm in diameter more than 2-fold by tuning ribosome concentration, without any discernable effect on the motion of molecules ≤5 nm. This change in ribosome concentration affected phase separation both in vitro and in vivo. Together, these results establish a role for mTORC1 in controlling both the mesoscale biophysical properties of the cytoplasm and biomolecular condensation.

Terminal deoxynucleotidyltransferase: the story of an untemplated DNA polymerase capable of DNA bridging and templated synthesis across strands

Terminal deoxynucleotidyltransferase (TdT) is a member of the polX family which is involved in DNA repair. It has been known for years as an untemplated DNA polymerase used during V(D)J recombination to generate diversity at the CDR3 region of immunoglobulins and T-cell receptors. Recently, however, TdT was crystallized in the presence of a complete DNA synapsis made of two double-stranded DNA (dsDNA), each with a 3' protruding end, and overlapping with only one micro-homology base-pair, thus giving structural insight for the first time into DNA synthesis across strands. It was subsequently shown that TdT indeed has an in trans template-dependent activity in the presence of an excess of the downstream DNA duplex. A possible biological role of this dual activity is discussed.

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.

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.

Pentameric ligand-gated ion channels are an important class of receptors that are involved in many neurological diseases. They have been extensively studied but a full understanding of their mechanism of action has yet to be achieved. In an effort to bypass obstacles in the research of human receptors, bacterial versions have been used to characterize the family’s structure-function relationship. One key bacterial receptor, known as GLIC, has lead the way in structural resolution of various mechanistic states along the gating pathway, yet its activation by protons is significantly less understood than its human counterparts. To define the site(s) involved in proton gating, we systematically mutated all titratable residues near the pH₅₀ of activation: Asp, Glu, and His. We determined that a previously established His residue in the transmembrane domain is structurally important but likely plays little or no role in proton gating. We instead found that proton activation is a complex multiple loci mechanism, with the key contribution stemming from the coupling interface between the extracellular and transmembrane domain, with E35 acting as a key proton-sensing residue.

Ab initio sampling of transition paths by conditioned Langevin dynamics

We propose a novel stochastic method to generate Brownian paths conditioned to start at an initial point and end at a given final point during a fixed time t_f under a given potential U(x). These paths are sampled with a probability given by the overdamped Langevin dynamics. We show that these paths can be exactly generated by a local stochastic partial differential equation. This equation cannot be solved in general but we present several approximations that are valid either in the low temperature regime or in the presence of barrier crossing. We show that this method warrants the generation of statistically independent transition paths. It is computationally very efficient. We illustrate the method first on two simple potentials, the two-dimensional Mueller potential and the Mexican hat potential, and then on the multi-dimensional problem of conformational transitions in proteins using the "Mixed Elastic Network Model" as a benchmark.

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.

Identification of a pre-active conformation of a pentameric channel receptor

Pentameric ligand-gated ion channels (pLGICs) mediate fast chemical signaling through global allosteric transitions. Despite the existence of several high-resolution structures of pLGICs, their dynamical properties remain elusive. Using the proton-gated channel GLIC, we engineered multiple fluorescent reporters, each incorporating a bimane and a tryptophan/tyrosine, whose close distance causes fluorescence quenching. We show that proton application causes a global compaction of the extracellular subunit interface, coupled to an outward motion of the M2-M3 loop near the channel gate. These movements are highly similar in lipid vesicles and detergent micelles. These reorganizations are essentially completed within 2 ms and occur without channel opening at low proton concentration, indicating that they report a pre-active intermediate state in the transition pathway toward activation. This provides a template to investigate the gating of eukaryotic neurotransmitter receptors, for which intermediate states also participate in activation.

DOI:http://dx.doi.org/10.7554/eLife.23955.001

In the nervous system, proteins of the pLGIC family are found in the membrane that surrounds each neuron. These proteins have channels that can allow ions to pass through the membrane and are responsible for transmitting electrical signals from one neuron to the next. Small molecules called neurotransmitters interact with the pLGICs to open or close the ion channel. If the ability of the pLGIC channels to open is altered, it can lead to behavioral changes like addiction, or diseases such as schizophrenia or epilepsy.

For a pLGIC channel to switch between the “open” and “closed” states, specific parts of the protein need to move in relation to each other. However, to study these transitions researchers have previously relied on comparing the three-dimensional structures of open and closed pLGICs extracted out of the cell membrane. Different techniques are needed to directly follow these movements within membranes.

Bacteria also have proteins belonging to the pLGIC family, and Menny et al. have now investigated one such bacterial protein to understand how pLGICs open. First, a small fluorescent molecule that glows differently if the environment around it changes was attached to various parts of the bacterial channel. These fluorescent markers revealed how several parts of the protein move and they also made it possible to measure how quickly these movements take place. Some of these movements happen before the channel opens, suggesting that the activation of this pLGIC protein happens in stages and involves the protein adopting a temporary intermediate state.

The next step will be to better understand the structure of the intermediate state, which could help us to understand how pLGICs work in the nervous systems of animals. In future this may aid the design of new drugs that can modify the activity of these channels in patients with neurological conditions or addictions.

DOI:http://dx.doi.org/10.7554/eLife.23955.002

The Renormalization Group and Its Applications to Generating Coarse-Grained Models of Large Biological Molecular Systems

Understanding the dynamics of biomolecules is the key to understanding their biological activities. Computational methods ranging from all-atom molecular dynamics simulations to coarse-grained normal-mode analyses based on simplified elastic networks provide a general framework to studying these dynamics. Despite recent successes in studying very large systems with up to a 100,000,000 atoms, those methods are currently limited to studying small- to medium-sized molecular systems due to computational limitations. One solution to circumvent these limitations is to reduce the size of the system under study. In this paper, we argue that coarse-graining, the standard approach to such size reduction, must define a hierarchy of models of decreasing sizes that are consistent with each other, i.e., that each model contains the information of the dynamics of its predecessor. We propose a new method, Decimate, for generating such a hierarchy within the context of elastic networks for normal-mode analysis. This method is based on the concept of the renormalization group developed in statistical physics. We highlight the details of its implementation, with a special focus on its scalability to large systems of up to millions of atoms. We illustrate its application on two large systems, the capsid of a virus and the ribosome translation complex. We show that highly decimated representations of those systems, containing down to 1% of their original number of atoms, still capture qualitatively and quantitatively their dynamics. Decimate is available as an OpenSource resource.

Cell-like pressure sensors reveal increase of mechanical stress towards the core of multicellular spheroids under compression

The surrounding microenvironment limits tumour expansion, imposing a compressive stress on the tumour, but little is known how pressure propagates inside the tumour. Here we present non-destructive cell-like microsensors to locally quantify mechanical stress distribution in three-dimensional tissue. Our sensors are polyacrylamide microbeads of well-defined elasticity, size and surface coating to enable internalization within the cellular environment. By isotropically compressing multicellular spheroids (MCS), which are spherical aggregates of cells mimicking a tumour, we show that the pressure is transmitted in a non-trivial manner inside the MCS, with a pressure rise towards the core. This observed pressure profile is explained by the anisotropic arrangement of cells and our results suggest that such anisotropy alone is sufficient to explain the pressure rise inside MCS composed of a single cell type. Furthermore, such pressure distribution suggests a direct link between increased mechanical stress and previously observed lack of proliferation within the spheroids core.

Comparative Normal Mode Analysis of the Dynamics of DENV and ZIKV Capsids

Key steps in the life cycle of a virus, such as the fusion event as the virus infects a host cell and its maturation process, relate to an intricate interplay between the structure and the dynamics of its constituent proteins, especially those that define its capsid, much akin to an envelope that protects its genomic material. We present a comprehensive, comparative analysis of such interplay for the capsids of two viruses from the flaviviridae family, Dengue (DENV) and Zika (ZIKV). We use for that purpose our own software suite, DD-NMA, which is based on normal mode analysis. We describe the elements of DD-NMA that are relevant to the analysis of large systems, such as virus capsids. In particular, we introduce our implementation of simplified elastic networks and justify their parametrization. Using DD-NMA, we illustrate the importance of packing interactions within the virus capsids on the dynamics of the E proteins of DENV and ZIKV. We identify differences between the computed atomic fluctuations of the E proteins in DENV and ZIKV and relate those differences to changes observed in their high resolution structures. We conclude with a discussion on additional analyses that are needed to fully characterize the dynamics of the two viruses.

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.

Plant programmed cell death from a chromatin point of view

Programmed cell death (PCD) is a ubiquitous genetically regulated process consisting of the activation of finely controlled signalling pathways that lead to cellular suicide. PCD can be part of a developmental programme (dPCD) or be triggered by environmental conditions (ePCD). In plant cells, as in animal cells, extensive chromatin condensation and degradation of the nuclear DNA are among the most conspicuous features of cells undergoing PCD. Changes in chromatin condensation could either reflect the structural changes required for internucleosomal fragmentation of nuclear DNA or relate to large-scale chromatin rearrangements associated with a major transcriptional switch occurring during cell death. The aim of this review is to give an update on plant PCD processes from a chromatin point of view. The first part will be dedicated to chromatin conformational changes associated with cell death observed in various developmental and physiological conditions, whereas the second part will be devoted to histone dynamics and DNA modifications associated with critical changes in genome expression during the cell death process.

Shared active site architecture between archaeal PolD and multi-subunit RNA polymerases revealed by X-ray crystallography

Archaeal replicative DNA polymerase D (PolD) constitute an atypical class of DNA polymerases made of a proofreading exonuclease subunit (DP1) and a larger polymerase catalytic subunit (DP2), both with unknown structures. We have determined the crystal structures of Pyrococcus abyssi DP1 and DP2 at 2.5 and 2.2 Å resolution, respectively, revealing a catalytic core strikingly different from all other known DNA polymerases (DNAPs). Rather, the PolD DP2 catalytic core has the same ‘double-psi β-barrel' architecture seen in the RNA polymerase (RNAP) superfamily, which includes multi-subunit transcriptases of all domains of life, homodimeric RNA-silencing pathway RNAPs and atypical viral RNAPs. This finding bridges together, in non-viral world, DNA transcription and DNA replication within the same protein superfamily. This study documents further the complex evolutionary history of the DNA replication apparatus in different domains of life and proposes a classification of all extant DNAPs.

The structures of many DNA polymerases is known, but PolD was a missing piece. Here, the authors report the crystal structure of this protein and use it to connect the DNA replication machinery with the transcription machinery in the same protein family.

Structural Basis for a New Templated Activity by Terminal Deoxynucleotidyl Transferase: Implications for V(D)J Recombination

Eukaryotic DNA polymerase of the polX family, such as pol μ and terminal deoxynucleotidyl transferase (TdT), are key components of the non-homologous end-joining or V(D)J recombination machinery, respectively. The established role of TdT is to add random nucleotides during V(D)J recombination. Here we show that TdT also has a templated-polymerase activity, similar to pol μ, in the presence of higher concentrations of a downstream DNA duplex, and performs a micro-homology single base-pair search to align the DNA synapsis. To understand the molecular basis of this alignment, we solve crystal structures of TdT with four DNA strands and study the influence of the 3' protruding end. Two mutations in TdT inspired by sequence alignments with pol μ further improve the templated activity. We propose that both templated and untemplated activities of TdT are needed to explain the distributions of lengths of N regions observed experimentally in T cell receptors and antibodies.

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

Structural Basis for Xenon Inhibition in a Cationic Pentameric Ligand-Gated Ion Channel

GLIC receptor is a bacterial pentameric ligand-gated ion channel whose action is inhibited by xenon. Xenon has been used in clinical practice as a potent gaseous anaesthetic for decades, but the molecular mechanism of interactions with its integral membrane receptor targets remains poorly understood. Here we characterize by X-ray crystallography the xenon-binding sites within both the open and "locally-closed" (inactive) conformations of GLIC. Major binding sites of xenon, which differ between the two conformations, were identified in three distinct regions that all belong to the trans-membrane domain of GLIC: 1) in an intra-subunit cavity, 2) at the interface between adjacent subunits, and 3) in the pore. The pore site is unique to the locally-closed form where the binding of xenon effectively seals the channel. A putative mechanism of the inhibition of GLIC by xenon is proposed, which might be extended to other pentameric cationic ligand-gated ion channels.