spotter: a single-nucleotide resolution stochastic simulation model of supercoiling-mediated transcription and translation in prokaryotes.

Stochastic simulation models have played an important role in efforts to understand the mechanistic basis of prokaryotic transcription and translation. Despite the fundamental linkage of these processes in bacterial cells, however, most simulation models have been limited to representations of either transcription or translation. In addition, the available simulation models typically either attempt to recapitulate data from single-molecule experiments without considering cellular-scale high-throughput sequencing data or, conversely, seek to reproduce cellular-scale data without paying close attention to many of the mechanistic details. To address these limitations, we here present spotter (Simulation of Prokaryotic Operon Transcription & Translation Elongation Reactions), a flexible, user-friendly simulation model that offers highly-detailed combined representations of prokaryotic transcription, translation, and DNA supercoiling. In incorporating nascent transcript and ribosomal profiling sequencing data, spotter provides a critical bridge between data collected in single-molecule experiments and data collected at the cellular scale. Importantly, in addition to rapidly generating output that can be aggregated for comparison with next-generation sequencing and proteomics data, spotter produces residue-level positional information that can be used to visualize individual simulation trajectories in detail. We anticipate that spotter will be a useful tool in exploring the interplay of processes that are crucially linked in prokaryotes.

The native 3D organization of bacterial polysomes.

Recent advances have led to insights into the structure of the bacterial ribosome, but little is known about the 3D organization of ribosomes in the context of translating polysomes. We employed cryoelectron tomography and a template-matching approach to map 70S ribosomes in vitrified bacterial translation extracts and in lysates of active E. coli spheroplasts. In these preparations, polysomal arrangements were observed in which neighboring ribosomes are densely packed and exhibit preferred orientations. Analysis of characteristic examples of polysomes reveals a staggered or pseudohelical organization of ribosomes along the mRNA trace, with the transcript being sequestered on the inside, the tRNA entrance sites being accessible, and the polypeptide exit sites facing the cytosol. Modeling of elongating nascent polypeptide chains suggests that this arrangement maximizes the distance between nascent chains on adjacent ribosomes, thereby reducing the probability of intermolecular interactions that would give rise to aggregation and limit productive folding.

The C-terminal region of translesion synthesis DNA polymerase η is partially unstructured and has high conformational flexibility.

Eukaryotic DNA polymerase η catalyzes translesion synthesis of thymine dimers and 8-oxoguanines. It is comprised of a polymerase domain and a C-terminal region, both of which are required for its biological function. The C-terminal region mediates interactions with proliferating cell nuclear antigen (PCNA) and other translesion synthesis proteins such as Rev1. This region contains a ubiquitin-binding/zinc-binding (UBZ) motif and a PCNA-interacting protein (PIP) motif. Currently little structural information is available for this region of polymerase η. Using a combination of approaches-including genetic complementation assays, X-ray crystallography, Langevin dynamics simulations, and small-angle X-ray scattering-we show that the C-terminal region is partially unstructured and has high conformational flexibility. This implies that the C-terminal region acts as a flexible tether linking the polymerase domain to PCNA thereby increasing its local concentration. Such tethering would facilitate the sampling of translesion synthesis polymerases to ensure that the most appropriate one is selected to bypass the lesion.

Features of genomic organization in a nucleotide-resolution molecular model of the Escherichia coli chromosome.

We describe structural models of the Escherichia coli chromosome in which the positions of all 4.6 million nucleotides of each DNA strand are resolved. Models consistent with two basic chromosomal orientations, differing in their positioning of the origin of replication, have been constructed. In both types of model, the chromosome is partitioned into plectoneme-abundant and plectoneme-free regions, with plectoneme lengths and branching patterns matching experimental distributions, and with spatial distributions of highly-transcribed chromosomal regions matching recent experimental measurements of the distribution of RNA polymerases. Physical analysis of the models indicates that the effective persistence length of the DNA and relative contributions of twist and writhe to the chromosome's negative supercoiling are in good correspondence with experimental estimates. The models exhibit characteristics similar to those of 'fractal globules,' and even the most genomically-distant parts of the chromosome can be physically connected, through paths combining linear diffusion and inter-segmental transfer, by an average of only ∼10 000 bp. Finally, macrodomain structures and the spatial distributions of co-expressed genes are analyzed: the latter are shown to depend strongly on the overall orientation of the chromosome. We anticipate that the models will prove useful in exploring other static and dynamic features of the bacterial chromosome.

Dynamic binding of replication protein a is required for DNA repair.

Replication protein A (RPA), the major eukaryotic single-stranded DNA (ssDNA) binding protein, is essential for replication, repair and recombination. High-affinity ssDNA-binding by RPA depends on two DNA binding domains in the large subunit of RPA. Mutation of the evolutionarily conserved aromatic residues in these two domains results in a separation-of-function phenotype: aromatic residue mutants support DNA replication but are defective in DNA repair. We used biochemical and single-molecule analyses, and Brownian Dynamics simulations to determine the molecular basis of this phenotype. Our studies demonstrated that RPA binds to ssDNA in at least two modes characterized by different dissociation kinetics. We also showed that the aromatic residues contribute to the formation of the longer-lived state, are required for stable binding to short ssDNA regions and are needed for RPA melting of partially duplex DNA structures. We conclude that stable binding and/or the melting of secondary DNA structures by RPA is required for DNA repair, including RAD51 mediated DNA strand exchange, but is dispensable for DNA replication. It is likely that the binding modes are in equilibrium and reflect dynamics in the RPA-DNA complex. This suggests that dynamic binding of RPA to DNA is necessary for different cellular functions.

Human heterochromatin protein 1α promotes nucleosome associations that drive chromatin condensation.

HP1(Hsα)-containing heterochromatin is located near centric regions of chromosomes and regulates DNA-mediated processes such as DNA repair and transcription. The higher-order structure of heterochromatin contributes to this regulation, yet the structure of heterochromatin is not well understood. We took a multidisciplinary approach to determine how HP1(Hsα)-nucleosome interactions contribute to the structure of heterochromatin. We show that HP1(Hsα) preferentially binds histone H3K9Me3-containing nucleosomal arrays in favor of non-methylated nucleosomal arrays and that nonspecific DNA interactions and pre-existing chromatin compaction promote binding. The chromo and chromo shadow domains of HP1(Hsα) play an essential role in HP1(Hsα)-nucleosome interactions, whereas the hinge region appears to have a less significant role. Electron microscopy of HP1(Hsα)-associated nucleosomal arrays showed that HP1(Hsα) caused nucleosome associations within an array, facilitating chromatin condensation. Differential sedimentation of HP1(Hsα)-associated nucleosomal arrays showed that HP1(Hsα) promotes interactions between arrays. These strand-to-strand interactions are supported by in vivo studies where tethering the Drosophila homologue HP1a to specific sites promotes interactions with distant chromosomal sites. Our findings demonstrate that HP1(Hsα)-nucleosome interactions cause chromatin condensation, a process that regulates many chromosome events.

ATP and AMP mutually influence their interaction with the ATP-binding cassette (ABC) adenylate kinase cystic fibrosis transmembrane conductance regulator (CFTR) at separate binding sites.

Cystic fibrosis transmembrane conductance regulator (CFTR) is an anion channel in the ATP-binding cassette (ABC) transporter protein family. In the presence of ATP and physiologically relevant concentrations of AMP, CFTR exhibits adenylate kinase activity (ATP + AMP &lrarr2; 2 ADP). Previous studies suggested that the interaction of nucleotide triphosphate with CFTR at ATP-binding site 2 is required for this activity. Two other ABC proteins, Rad50 and a structural maintenance of chromosome protein, also have adenylate kinase activity. All three ABC adenylate kinases bind and hydrolyze ATP in the absence of other nucleotides. However, little is known about how an ABC adenylate kinase interacts with ATP and AMP when both are present. Based on data from non-ABC adenylate kinases, we hypothesized that ATP and AMP mutually influence their interaction with CFTR at separate binding sites. We further hypothesized that only one of the two CFTR ATP-binding sites is involved in the adenylate kinase reaction. We found that 8-azidoadenosine 5'-triphosphate (8-N3-ATP) and 8-azidoadenosine 5'-monophosphate (8-N3-AMP) photolabeled separate sites in CFTR. Labeling of the AMP-binding site with 8-N3-AMP required the presence of ATP. Conversely, AMP enhanced photolabeling with 8-N3-ATP at ATP-binding site 2. The adenylate kinase active center probe P(1),P(5)-di(adenosine-5') pentaphosphate interacted simultaneously with an AMP-binding site and ATP-binding site 2. These results show that ATP and AMP interact with separate binding sites but mutually influence their interaction with the ABC adenylate kinase CFTR. They further indicate that the active center of the adenylate kinase comprises ATP-binding site 2.

Improving Signal and Transit Peptide Predictions Using AlphaFold2-predicted Protein Structures.

Many proteins contain cleavable signal or transit peptides that direct them to their final subcellular locations. Such peptides are usually predicted from sequence alone using methods such as TargetP 2.0 and SignalP 6.0. While these methods are usually very accurate, we show here that an analysis of a protein's AlphaFold2-predicted structure can often be used to identify false positive predictions. We start by showing that when given a protein's full-length sequence, AlphaFold2 builds experimentally annotated signal and transit peptides in orientations that point away from the main body of the protein. This indicates that AlphaFold2 correctly identifies that a signal is not destined to be part of the mature protein's structure and suggests, as a corollary, that predicted signals that AlphaFold2 folds with high confidence into the main body of the protein are likely to be false positives. To explore this idea, we analyzed predicted signal peptides in 48 proteomes made available in DeepMind's AlphaFold2 database (https://alphafold.ebi.ac.uk). Applying TargetP 2.0 and SignalP 6.0 to the 561,562 proteins in the database results in 95,236 being predicted to contain a cleavable signal or transit peptide. In 95.1% of these cases, the AlphaFold2 structure of the full-length protein is fully consistent with the prediction of TargetP 2.0 or SignalP 6.0. In the remaining 4.9% of cases where the AlphaFold2 structure does not appear consistent with the prediction, the signal is often only predicted with low confidence. The potential false positives identified here may be useful for training even more accurate signal prediction methods.

AutoRNC: An automated modeling program for building atomic models of ribosome-nascent chain complexes.

The interpretation of experimental studies of co-translational protein folding often benefits from the use of computational methods that seek to model or simulate the nascent chain and its interactions with the ribosome. Building realistic 3D models of ribosome-nascent chain (RNC) constructs often requires expert knowledge, so to circumvent this issue, we describe here AutoRNC, an automated modeling program capable of constructing large numbers of plausible atomic models of RNCs within minutes. AutoRNC takes input from the user specifying any regions of the nascent chain that contain secondary or tertiary structure and attempts to build conformations compatible with those specifications-and with the constraints imposed by the ribosome-by sampling and progressively piecing together dipeptide conformations extracted from the Research Collaboratory for Structural Bioinformatics Protein Data Bank (RCSB). Despite using only modest computational resources, we show here that AutoRNC can build plausible conformations for a wide range of RNC constructs for which experimental data have already been reported.

Flexibility of the bacterial chaperone trigger factor in microsecond-timescale molecular dynamics simulations.

The bacterial chaperone trigger factor (TF) is the first chaperone to be encountered by a nascent protein chain as it emerges from the ribosome exit tunnel. Experimental results suggest that TF possesses considerable conformational flexibility, and in an attempt to provide an atomic-level view of this flexibility, we have performed independent 1.5-µs molecular dynamics simulations of TF in explicit solvent using two different simulation force fields (OPLS-AA/L and AMBER ff99SB-ILDN). Both simulations indicate that TF possesses tremendous flexibility, with huge excursions from the crystallographic conformation caused by reorientations of the protein's constituent domains; both simulations also predict the formation of extensive contacts between TF's PPIase domain and the Arm 1 domain that is involved in nascent-chain binding. In the OPLS simulation, however, TF rapidly settles into a very compact conformation that persists for at least 1 µs, whereas in the AMBER simulation, it remains highly dynamic; additional simulations in which the two force fields were swapped suggest that these differences are at least partly attributable to sampling issues. The simulation results provide potential rationalizations of a number of experimental observations regarding TF's conformational behavior and have implications for using simulations to model TF's function on translating ribosomes.

Easy Removal of Steric Clashes and Entanglements in Macromolecular Systems by Temporary Addition of a Fourth Spatial Dimension.

When models of complicated macromolecular systems are constructed, it is common to inadvertently include either gross steric clashes or entanglements of extended loop regions. Removing these problems with conventional energy minimization or dynamics algorithms can often be difficult. Here I show that one easy alternative is to temporarily add an extra spatial dimension and to displace atoms or molecules along this fourth dimension such that the distances between atoms, when measured in 4D, are no longer considered clashing. Adding in half-harmonic potential functions to mimic walls in this 4 th dimension, and then moving these walls toward each other, has the effect of decreasing the space available in the 4 th dimension and drives atoms to avoid each other in 3D. I illustrate the method with three examples: two showing how interlocked ring polymers can be easily disentangled from each other in both 2D and 3D, and one showing how ten identical coarse-grained protein models, all placed at the same point in 3D space, can be separated from each other, without distorting their structures, during the course of a single energy minimization. A sample program implementing the method is available that can be easily adapted to other situations.

AutoRNC: an automated modeling program for building atomic models of ribosome-nascent chain complexes.

The interpretation of experimental studies of co-translational protein folding often benefits from the use of computational methods that seek to model the nascent chain and its interactions with the ribosome. Ribosome-nascent chain (RNC) constructs studied experimentally can vary significantly in size and the extent to which they contain secondary and tertiary structure, and building realistic 3D models of them therefore often requires expert knowledge. To circumvent this issue, we describe here AutoRNC, an automated modeling program capable of constructing large numbers of plausible atomic models of RNCs within minutes. AutoRNC takes input from the user specifying any regions of the nascent chain that contain secondary or tertiary structure and attempts to build conformations compatible with those specifications - and with the constraints imposed by the ribosome - by sampling and progressively piecing together dipeptide conformations extracted from the RCSB. We first show that conformations of completely unfolded proteins built by AutoRNC in the absence of the ribosome have radii of gyration that match well with the corresponding experimental data. We then show that AutoRNC can build plausible conformations for a wide range of RNC constructs for which experimental data have already been reported. Since AutoRNC requires only modest computational resources, we anticipate that it will prove to be a useful hypothesis generator for experimental studies, for example, in providing indications of whether designed constructs are likely to be capable of folding, as well as providing useful starting points for downstream atomic or coarse-grained simulations of the conformational dynamics of RNCs.

Orientationally Averaged Version of the Rotne-Prager-Yamakawa Tensor Provides a Fast but Still Accurate Treatment of Hydrodynamic Interactions in Brownian Dynamics Simulations of Biological Macromolecules.

The Brownian dynamics (BD) simulation technique is widely used to model the diffusive and conformational dynamics of complex systems comprising biological macromolecules. For the diffusive properties of macromolecules to be described correctly by BD simulations, it is necessary to include hydrodynamic interactions (HIs). When modeled at the Rotne-Prager-Yamakawa (RPY) level of theory, for example, the translational and rotational diffusion coefficients of isolated macromolecules can be accurately reproduced; when HIs are neglected, however, diffusion coefficients can be underestimated by an order of magnitude or more. The principal drawback to the inclusion of HIs in BD simulations is their computational expense, and several previous studies have sought to accelerate their modeling by developing fast approximations for the calculation of the correlated random displacements. Here, we explore the use of an alternative way to accelerate the calculation of HIs, i.e., by replacing the full RPY tensor with an orientationally averaged (OA) version which retains the distance dependence of the HIs but averages out their orientational dependence. We seek here to determine whether such an approximation can be justified in application to the modeling of typical proteins and RNAs. We show that the use of an OA-RPY tensor allows translational diffusion of macromolecules to be modeled with very high accuracy at the cost of rotational diffusion being underestimated by ∼25%. We show that this finding is independent of the type of macromolecule simulated and the level of structural resolution employed in the models. We also show, however, that these results are critically dependent on the inclusion of a non-zero term that describes the divergence of the diffusion tensor: when this term is omitted from simulations that use the OA-RPY model, unfolded macromolecules undergo rapid collapse. Our results indicate that the orientationally averaged RPY tensor is likely to be a useful, fast, approximate way of including HIs in BD simulations of intermediate-scale systems.

An Orientationally Averaged Version of the Rotne-Prager-Yamakawa Tensor Provides A Fast But Still Accurate Treatment Of Hydrodynamic Interactions In Brownian Dynamics Simulations Of Biological Macromolecules.

The Brownian dynamics (BD) simulation technique is widely used to model the diffusive and conformational dynamics of complex systems comprising biological macromolecules. For the diffusive properties of macromolecules to be described correctly by BD simulations, it is necessary to include hydrodynamic interactions (HI). When modeled at the Rotne-Prager-Yamakawa (RPY) level of theory, for example, the translational and rotational diffusion coefficients of isolated macromolecules can be accurately reproduced; when HIs are neglected, however, diffusion coefficients can be underestimated by an order of magnitude or more. The principal drawback to the inclusion of HIs in BD simulations is their computational expense, and several previous studies have sought to accelerate their modeling by developing fast approximations for the calculation of the correlated random displacements. Here we explore the use of an alternative way to accelerate calculation of HIs, i.e., by replacing the full RPY tensor with an orientationally averaged (OA) version which retains the distance dependence of the HIs but averages out their orientational dependence. We seek here to determine whether such an approximation can be justified in application to the modeling of typical proteins and RNAs. We show that the use of an OA RPY tensor allows translational diffusion of macromolecules to be modeled with very high accuracy at the cost of rotational diffusion being underestimated by ∼25%. We show that this finding is independent of the type of macromolecule simulated and the level of structural resolution employed in the models. We also show, however, that these results are critically dependent on the inclusion of a non-zero term that describes the divergence of the diffusion tensor: when this term is omitted from simulations that use the OA RPY model, unfolded macromolecules undergo rapid collapse. Our results indicate that the orientationally averaged RPY tensor is likely to be a useful, fast approximate way of including HIs in BD simulations of intermediate-scale systems.

spotter : A single-nucleotide resolution stochastic simulation model of supercoiling-mediated transcription and translation in prokaryotes.

Stochastic simulation models have played an important role in efforts to understand the mechanistic basis of prokaryotic transcription and translation. Despite the fundamental linkage of these processes in bacterial cells, however, most simulation models have been limited to representations of either transcription or translation. In addition, the available simulation models typically either attempt to recapitulate data from single-molecule experiments without considering cellular-scale high-throughput sequencing data or, conversely, seek to reproduce cellular-scale data without paying close attention to many of the mechanistic details. To address these limitations, we here present spotter (Simulation of Prokaryotic Operon Transcription & Translation Elongation Reactions), a flexible, user-friendly simulation model that offers highly-detailed combined representations of prokaryotic transcription, translation, and DNA supercoiling. In incorporating nascent transcript and ribosomal profiling sequencing data, spotter provides a critical bridge between data collected in single-molecule experiments and data collected at the cellular scale. Importantly, in addition to rapidly generating output that can be aggregated for comparison with next-generation sequencing and proteomics data, spotter produces residue-level positional information that can be used to visualize individual simulation trajectories in detail. We anticipate that spotter will be a useful tool in exploring the interplay of processes that are crucially linked in prokaryotes.

Modeling the 3D structure and conformational dynamics of very large RNAs using coarse-grained molecular simulations.

We describe a computational approach to building and simulating realistic 3D models of very large RNA molecules (>1000 nucleotides) at a resolution of one "bead" per nucleotide. The method starts with a predicted secondary structure and uses several stages of energy minimization and Brownian dynamics (BD) simulation to build 3D models. A key step in the protocol is the temporary addition of a 4 th spatial dimension that allows all predicted helical elements to become disentangled from each other in an effectively automated way. We then use the resulting 3D models as input to Brownian dynamics simulations that include hydrodynamic interactions (HIs) that allow the diffusive properties of the RNA to be modelled as well as enabling its conformational dynamics to be simulated. To validate the dynamics part of the method, we first show that when applied to small RNAs with known 3D structures the BD-HI simulation models accurately reproduce their experimental hydrodynamic radii (Rh). We then apply the modelling and simulation protocol to a variety of RNAs for which experimental Rh values have been reported ranging in size from 85 to 3569 nucleotides. We show that the 3D models, when used in BD-HI simulations, produce hydrodynamic radii that are usually in good agreement with experimental estimates for RNAs that do not contain tertiary contacts that persist even under very low salt conditions. Finally, we show that sampling of the conformational dynamics of large RNAs on timescales of 100 µs is computationally feasible with BD-HI simulations.

Atomic Models of All Major Trans-Envelope Complexes Involved in Lipid Trafficking in Escherichia Coli Constructed Using a Combination of AlphaFold2, AF2Complex, and Membrane Morphing Simulations.

In Gram-negative bacteria, several trans-envelope complexes (TECs) have been identified that span the periplasmic space in order to facilitate lipid transport between the inner- and outer- membranes. While partial or near-complete structures of some of these TECs have been solved by conventional experimental techniques, most remain incomplete. Here we describe how a combination of computational approaches, constrained by experimental data, can be used to build complete atomic models for four TECs implicated in lipid transport in Escherichia coli . We use DeepMind's protein structure prediction algorithm, AlphaFold2, and a variant of it designed to predict protein complexes, AF2Complex, to predict the oligomeric states of key components of TECs and their likely interfaces with other components. After obtaining initial models of the complete TECs by superimposing predicted structures of subcomplexes, we use the membrane orientation prediction algorithm OPM to predict the likely orientations of the inner- and outer- membrane components in each TEC. Since, in all cases, the predicted membrane orientations in these initial models are tilted relative to each other, we devise a novel molecular mechanics-based strategy that we call "membrane morphing" that adjusts each TEC model until the two membranes are properly aligned with each other and separated by a distance consistent with estimates of the periplasmic width in E. coli . The study highlights the potential power of combining computational methods, operating within limits set by both experimental data and by cell physiology, for producing useable atomic structures of very large protein complexes.

A PEG-based oligomer as a backbone replacement for surface-exposed loops in a protein tertiary structure.

PEGged out: Poly(ethylene glycol), a simple biocompatible polymer, can replace natural loop segments in a 56-residue protein domain with a well-defined tertiary structure. Biophysical characterization of chimeras of the protein GB1 coupled with molecular dynamics simulations show that PEG enhances local backbone torsional freedom without compromising the overall protein fold or function.

Identification of Iron-Sulfur (Fe-S) Cluster and Zinc (Zn) Binding Sites Within Proteomes Predicted by DeepMind's AlphaFold2 Program Dramatically Expands the Metalloproteome.

DeepMind's AlphaFold2 software has ushered in a revolution in high quality, 3D protein structure prediction. In very recent work by the DeepMind team, structure predictions have been made for entire proteomes of twenty-one organisms, with >360,000 structures made available for download. Here we show that thousands of novel binding sites for iron-sulfur (Fe-S) clusters and zinc (Zn) ions can be identified within these predicted structures by exhaustive enumeration of all potential ligand-binding orientations. We demonstrate that AlphaFold2 routinely makes highly specific predictions of ligand binding sites: for example, binding sites that are comprised exclusively of four cysteine sidechains fall into three clusters, representing binding sites for 4Fe-4S clusters, 2Fe-2S clusters, or individual Zn ions. We show further: (a) that the majority of known Fe-S cluster and Zn binding sites documented in UniProt are recovered by the AlphaFold2 structures, (b) that there are occasional disputes between AlphaFold2 and UniProt with AlphaFold2 predicting highly plausible alternative binding sites, (c) that the Fe-S cluster binding sites that we identify in E. coli agree well with previous bioinformatics predictions, (d) that cysteines predicted here to be part of ligand binding sites show little overlap with those shown via chemoproteomics techniques to be highly reactive, and (e) that AlphaFold2 occasionally appears to build erroneous disulfide bonds between cysteines that should instead coordinate a ligand. These results suggest that AlphaFold2 could be an important tool for the functional annotation of proteomes, and the methodology presented here is likely to be useful for predicting other ligand-binding sites.

Design principles that protect the proteasome from self-destruction.

The proteasome is a powerful intracellular protease that can degrade effectively any protein, self or foreign, for regulation, quality control, or immune response. Proteins are targeted for degradation by localizing them to the proteasome, typically by ubiquitin tags. At the same time, the proteasome is built from ~33 subunits, and their assembly into the complex and activity are tuned by post-translational modifications on long disordered regions on the subunits. Molecular modeling and biochemical experiments show that some of the disordered regions of proteasomal subunits can access the substrate recognition sites. All disordered regions tested, independent of location, are constructed from amino acid sequences that escape recognition. Replacing a disordered region with a sequence that is recognized by the proteasome leads to self-degradation and, in the case of an essential subunit, cell death.