The ribosome structure controls and directs mRNA entry, translocation and exit dynamics

The energy landscape of the ribosome

The ribosome is a prototypical assembly that can be used to establish general principles and techniques for the study of biological molecular machines. Motivated by the fact that the dynamics of every biomolecule is governed by an underlying energy landscape, there has been great interest to understand and quantify ribosome energetics. In the present review, we will focus on theoretical and computational strategies for probing the interactions that shape the energy landscape of the ribosome, with an emphasis on more recent studies of the elongation cycle. These efforts include the application of quantum mechanical methods for describing chemical kinetics, as well as classical descriptions to characterize slower (microsecond to millisecond) large-scale (10-100 Å) rearrangements, where motion is described in terms of diffusion across an energy landscape. Together, these studies provide broad insights into the factors that control a diverse range of dynamical processes in this assembly.

Antiviral responses versus virus-induced cellular shutoff: a game of thrones between influenza A virus NS1 and SARS-CoV-2 Nsp1

Following virus recognition of host cell receptors and viral particle/genome internalization, viruses replicate in the host via hijacking essential host cell machinery components to evade the provoked antiviral innate immunity against the invading pathogen. Respiratory viral infections are usually acute with the ability to activate pattern recognition receptors (PRRs) in/on host cells, resulting in the production and release of interferons (IFNs), proinflammatory cytokines, chemokines, and IFN-stimulated genes (ISGs) to reduce virus fitness and mitigate infection. Nevertheless, the game between viruses and the host is a complicated and dynamic process, in which they restrict each other via specific factors to maintain their own advantages and win this game. The primary role of the non-structural protein 1 (NS1 and Nsp1) of influenza A viruses (IAV) and the pandemic severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), respectively, is to control antiviral host-induced innate immune responses. This review provides a comprehensive overview of the genesis, spatial structure, viral and cellular interactors, and the mechanisms underlying the unique biological functions of IAV NS1 and SARS-CoV-2 Nsp1 in infected host cells. We also highlight the role of both non-structural proteins in modulating viral replication and pathogenicity. Eventually, and because of their important role during viral infection, we also describe their promising potential as targets for antiviral therapy and the development of live attenuated vaccines (LAV). Conclusively, both IAV NS1 and SARS-CoV-2 Nsp1 play an important role in virus-host interactions, viral replication, and pathogenesis, and pave the way to develop novel prophylactic and/or therapeutic interventions for the treatment of these important human respiratory viral pathogens.

Ribosomal Protein uS5 and Friends: Protein-Protein Interactions Involved in Ribosome Assembly and Beyond

Ribosomal proteins are fundamental components of the ribosomes in all living cells. The ribosomal protein uS5 (Rps2) is a stable component of the small ribosomal subunit within all three domains of life. In addition to its interactions with proximal ribosomal proteins and rRNA inside the ribosome, uS5 has a surprisingly complex network of evolutionarily conserved non-ribosome-associated proteins. In this review, we focus on a set of four conserved uS5-associated proteins: the protein arginine methyltransferase 3 (PRMT3), the programmed cell death 2 (PDCD2) and its PDCD2-like (PDCD2L) paralog, and the zinc finger protein, ZNF277. We discuss recent work that presents PDCD2 and homologs as a dedicated uS5 chaperone and PDCD2L as a potential adaptor protein for the nuclear export of pre-40S subunits. Although the functional significance of the PRMT3-uS5 and ZNF277-uS5 interactions remain elusive, we reflect on the potential roles of uS5 arginine methylation by PRMT3 and on data indicating that ZNF277 and PRMT3 compete for uS5 binding. Together, these discussions highlight the complex and conserved regulatory network responsible for monitoring the availability and the folding of uS5 for the formation of 40S ribosomal subunits and/or the role of uS5 in potential extra-ribosomal functions.

Sarecycline inhibits protein translation in Cutibacterium acnes 70S ribosome using a two-site mechanism

Acne vulgaris is a chronic disfiguring skin disease affecting ∼1 billion people worldwide, often having persistent negative effects on physical and mental health. The Gram-positive anaerobe, Cutibacterium acnes is implicated in acne pathogenesis and is, therefore, a main target for antibiotic-based acne therapy. We determined a 2.8-Å resolution structure of the 70S ribosome of Cutibacterium acnes by cryogenic electron microscopy and discovered that sarecycline, a narrow-spectrum antibiotic against Cutibacterium acnes, may inhibit two active sites of this bacterium's ribosome in contrast to the one site detected previously on the model ribosome of Thermus thermophilus. Apart from the canonical binding site at the mRNA decoding center, the second binding site for sarecycline exists at the nascent peptide exit tunnel, reminiscent of the macrolides class of antibiotics. The structure also revealed Cutibacterium acnes-specific features of the ribosomal RNA and proteins. Unlike the ribosome of the Gram-negative bacterium Escherichia coli, Cutibacterium acnes ribosome has two additional proteins, bS22 and bL37, which are also present in the ribosomes of Mycobacterium smegmatis and Mycobacterium tuberculosis. We show that bS22 and bL37 have antimicrobial properties and may be involved in maintaining the healthy homeostasis of the human skin microbiome.

Unravelling How Single-Stranded DNA Binding Protein Coordinates DNA Metabolism Using Single-Molecule Approaches

Single-stranded DNA-binding proteins (SSBs) play vital roles in DNA metabolism. Proteins of the SSB family exclusively and transiently bind to ssDNA, preventing the DNA double helix from re-annealing and maintaining genome integrity. In the meantime, they interact and coordinate with various proteins vital for DNA replication, recombination, and repair. Although SSB is essential for DNA metabolism, proteins of the SSB family have been long described as accessory players, primarily due to their unclear dynamics and mechanistic interaction with DNA and its partners. Recently-developed single-molecule tools, together with biochemical ensemble techniques and structural methods, have enhanced our understanding of the different coordination roles that SSB plays during DNA metabolism. In this review, we discuss how single-molecule assays, such as optical tweezers, magnetic tweezers, Förster resonance energy transfer, and their combinations, have advanced our understanding of the binding dynamics of SSBs to ssDNA and their interaction with other proteins partners. We highlight the central coordination role that the SSB protein plays by directly modulating other proteins' activities, rather than as an accessory player. Many possible modes of SSB interaction with protein partners are discussed, which together provide a bigger picture of the interaction network shaped by SSB.

Modeling the ribosome as a bipartite graph

Developing mathematical representations of biological systems that can allow predictions is a challenging and important research goal. It is demonstrated here how the ribosome, the nano-machine responsible for synthesizing all proteins necessary for cellular life, can be represented as a bipartite network. Ten ribosomal structures from Bacteria and six from Eukarya are explored. Ribosomal networks are found to exhibit unique properties despite variations in the nodes and edges of the different graphs. The ribosome is shown to exhibit very large topological redundancies, demonstrating mathematical resiliency. These results can potentially explain how it can function consistently despite changes in composition and connectivity. Furthermore, this representation can be used to analyze ribosome function within the large machinery of network theory, where the degrees of freedom are the possible interactions, and can be used to provide new insights for translation regulation and therapeutics.

Identifying Strategies to Experimentally Probe Multidimensional Dynamics in the Ribosome

The ribosome is a complex biomolecular machine that utilizes large-scale conformational rearrangements to synthesize proteins. For example, during the elongation cycle, the "head" domain of the ribosomal small subunit (SSU) is known to undergo transient rotation events that allow for movement of tRNA molecules (i.e., translocation). While the head may exhibit rigid-body-like properties, the precise relationship between experimentally accessible probes and multidimensional rotations has yet to be established. To address this gap, we perform molecular dynamics simulations of the translocation step of the elongation cycle in the ribosome, where the SSU head spontaneously undergoes rotation and tilt-like motions. With this data set (1250 simulated events), we used statistical and information-theory-based measures to identify possible single-molecule probes that can isolate SSU head rotation and head tilting. This analysis provides a molecular interpretation for previous single-molecule measurements, while establishing a framework for the design of next-generation experiments that may precisely probe the mechanistic and kinetic aspects of the ribosome.

The energetics of subunit rotation in the ribosome

Protein synthesis in the cell is controlled by an elaborate sequence of conformational rearrangements in the ribosome. The composition of a ribosome varies by species, though they typically contain ∼ 50-100 RNA and protein molecules. While advances in structural techniques have revolutionized our understanding of long-lived conformational states, a vast range of transiently visited configurations can not be directly observed. In these cases, computational/simulation methods can be used to understand the mechanical properties of the ribosome. Insights from these approaches can then help guide next-generation experimental measurements. In this short review, we discuss theoretical strategies that have been deployed to quantitatively describe the energetics of collective rearrangements in the ribosome. We focus on efforts to probe large-scale subunit rotation events, which involve the coordinated displacement of large numbers of atoms (tens of thousands). These investigations are revealing how the molecular structure of the ribosome encodes the mechanical properties that control large-scale dynamics.

Comparative Analysis of Structural and Dynamical Features of Ribosome Upon Association With mRNA Reveals Potential Role of Ribosomal Proteins

Ribosomes play a critical role in maintaining cellular proteostasis. The binding of messenger RNA (mRNA) to the ribosome regulates kinetics of protein synthesis. To generate an understanding of the structural, mechanistic, and dynamical features of mRNA recognition in the ribosome, we have analysed mRNA-protein interactions through a structural comparison of the ribosomal complex in the presence and absence of mRNA. To do so, we compared the 3-Dimensional (3D) structures of components of the two assembly structures and analysed their structural differences because of mRNA binding, using elastic network models and structural network-based analysis. We observe that the head region of 30S ribosomal subunit undergoes structural displacement and subunit rearrangement to accommodate incoming mRNA. We find that these changes are observed in proteins that lie far from the mRNA-protein interface, implying allostery. Further, through perturbation response scanning, we show that the proteins S13, S19, and S20 act as universal sensors that are sensitive to changes in the inter protein network, upon binding of 30S complex with mRNA and other initiation factors. Our study highlights the significance of mRNA binding in the ribosome complex and identifies putative allosteric sites corresponding to alterations in structure and/or dynamics, in regions away from mRNA binding sites in the complex. Overall, our work provides fresh insights into mRNA association with the ribosome, highlighting changes in the interactions and dynamics of the ribosome assembly because of the binding.

Programmed -1 ribosomal frameshifting from the perspective of the conformational dynamics of mRNA and ribosomes

Programmed -1 ribosomal frameshifting (-1 PRF) is a translation mechanism that regulates the relative expression level of two proteins encoded on the same messenger RNA (mRNA). This regulation is commonly used by viruses such as coronaviruses and retroviruses but rarely by host human cells, and for this reason, it has long been considered as a therapeutic target for antiviral drug development. Understanding the molecular mechanism of -1 PRF is one step toward this goal. Minus-one PRF occurs with a certain efficiency when translating ribosomes encounter the specialized mRNA signal consisting of the frameshifting site and a downstream stimulatory structure, which impedes translocation of the ribosome. The impeded ribosome can still undergo profound conformational changes to proceed with translocation; however, some of these changes may be unique and essential to frameshifting. In addition, most stimulatory structures exhibit conformational dynamics and sufficient mechanical strength, which, when under the action of ribosomes, may in turn further promote -1 PRF efficiency. In this review, we discuss how the dynamic features of ribosomes and mRNA stimulatory structures may influence the occurrence of -1 PRF and propose a hypothetical frameshifting model that recapitulates the role of conformational dynamics.

Local and Global Motions Underlying Antibiotic Binding in Bacterial Ribosome

The bacterial ribosome is one of the most important targets in the treatment of infectious diseases. As antibiotic resistance in bacteria poses a growing threat, a significant amount of effort is concentrated on exploring new drug-binding sites where testable predictions are of significance. Here, we study the dynamics of a ribosomal complex and 67 small and large subunits of the ribosomal crystal structures (64 antibiotic-bound, 3 antibiotic-free) from Deinococcus radiodurans, Escherichia coli, Haloarcula marismortui, and Thermus thermophilus by the Gaussian network model. Interestingly, a network of nucleotides coupled in high-frequency fluctuations reveals known antibiotic-binding sites. These sites are seen to locate at the interface of dynamic domains that have an intrinsic dynamic capacity to interfere with functional globular motions. The nucleotides and the residues fluctuating in the fast and slow modes of motion thus have promise for plausible antibiotic-binding and allosteric sites that can alter antibiotic binding and resistance. Overall, the present analysis brings a new dynamic perspective to the long-discussed link between small-molecule binding and large conformational changes of the supramolecule.

Network theory of the bacterial ribosome

In the past two decades, research into the biochemical, biophysical and structural properties of the ribosome have revealed many different steps of protein translation. Nevertheless, a complete understanding of how they lead to a rapid and accurate protein synthesis still remains a challenge. Here we consider a coarse network analysis in the bacterial ribosome formed by the connectivity between ribosomal (r) proteins and RNAs at different stages in the elongation cycle. The ribosomal networks are found to be dis-assortative and small world, implying that the structure allows for an efficient exchange of information between distant locations. An analysis of centrality shows that the second and fifth domains of 23S rRNA are the most important elements in all of the networks. Ribosomal protein hubs connect to much fewer nodes but are shown to provide important connectivity within the network (high closeness centrality). A modularity analysis reveals some of the different functional communities, indicating some known and some new possible communication pathways Our mathematical results confirm important communication pathways that have been discussed in previous research, thus verifying the use of this technique for representing the ribosome, and also reveal new insights into the collective function of ribosomal elements.

Selective translation by alternative bacterial ribosomes

Many organisms, including bacteria, code for multiple paralogues of some ribosomal protein subunits. The relative contribution of these alternative subunits to ribosome function and protein synthesis is unknown and controversial. Many studies on alternative ribosomes have been confounded by isolation of alternative and canonical ribosomes from different strains or growth conditions, potentially confounding results. Here, we show that one form of alternative ribosome from Mycobacterium smegmatis has a distinct translational profile compared with canonical ribosomes purified from an identical cellular context. We also identify a role for alternative ribosomes in iron homeostasis. Given the prevalence of alternative ribosomal genes in diverse organisms, our study suggests that alternative ribosomes may represent a further layer of regulation of gene translation.

Alternative ribosome subunit proteins are prevalent in the genomes of diverse bacterial species, but their functional significance is controversial. Attempts to study microbial ribosomal heterogeneity have mostly relied on comparing wild-type strains with mutants in which subunits have been deleted, but this approach does not allow direct comparison of alternate ribosome isoforms isolated from identical cellular contexts. Here, by simultaneously purifying canonical and alternative RpsR ribosomes from Mycobacterium smegmatis, we show that alternative ribosomes have distinct translational features compared with their canonical counterparts. Both alternative and canonical ribosomes actively take part in protein synthesis, although they translate a subset of genes with differential efficiency as measured by ribosome profiling. We also show that alternative ribosomes have a relative defect in initiation complex formation. Furthermore, a strain of M. smegmatis in which the alternative ribosome protein operon is deleted grows poorly in iron-depleted medium, uncovering a role for alternative ribosomes in iron homeostasis. Our work confirms the distinct and nonredundant contribution of alternative bacterial ribosomes for adaptation to hostile environments.

Resolution-exchanged structural modeling and simulations jointly unravel that subunit rolling underlies the mechanism of programmed ribosomal frameshifting

MOTIVATION

Programmed ribosomal frameshifting (PRF) is widely used by viruses and bacteria to produce different proteins from a single mRNA template. How steric hindrance of a PRF-stimulatory mRNA structure transiently modifies the conformational dynamics of the ribosome, and thereby allows tRNA slippage, remains elusive.

RESULTS

Here, we leverage linear response theories and resolution-exchanged simulations to construct a structural/dynamics model that connects and rationalizes existing structural, single-molecule and mutagenesis data by resolution-exchanged structural modelling and simulations. Our combined theoretical techniques provide a temporal and spatial description of PRF with unprecedented mechanistic details. We discover that ribosomal unfolding of the PRF-stimulating pseudoknot exerts resistant forces on the mRNA entrance of the ribosome, and thereby drives 30S subunit rolling. Such motion distorts tRNAs, leads to tRNA slippage, and in turn serves as a delicate control of cis-element's unwinding forces over PRF.

AVAILABILITY AND IMPLEMENTATION

All the simulation scripts and computational implementations of our methods/analyses (including linear response theory) are included in the bioStructureM suite, provided through GitHub at https://github.com/Yuan-Yu/bioStructureM.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Unique structural features of the Mycobacterium ribosome

Protein synthesis in all the living cells is mediated by a large protein-RNA complex called the ribosome. These macromolecular complexes can range from 2.5 (prokaryotes) to 4.2 MDa. (eukaryotes) in size and undergo various conformational transitions during protein synthesis to translate the genetic code into the nascent polypeptide chains. Recent advances in cryo-electron microscopy (cryo-EM) and image processing methods have provided numerous detailed structures of ribosomes from diverse sources and in different conformational states resolved to near-atomic resolutions. These structures have not only helped in better understanding of the translational mechanism but also revealed species-specific variations or adaptations in the ribosome structures. Structural investigations of the ribosomes from Mycobacterium smegmatis (Msm) and its closely related pathogenic Mycobacterium tuberculosis (Mtb) lead to the identification of two additional ribosomal proteins named as bS22 and bL37 and several unique extensions in ribosomal-protein and ribosomal-RNA. Hibernation Promoting Factor (HPF) bound structure of Msm ribosome, termed as the hibernating ribosome, possibly indicates a new mechanism of ribosome protection during dormancy. These studies enabled the identification of the mycobacteria-specific ribosomal features and provides an opportunity to understand their function and target them for further drug-discovery purposes. Here we review the unique structural features identified in Msm ribosome and their possible implications in comparison to a well-studied Escherichia coli (Ec) ribosome.

An Efficient Timer and Sizer of Biomacromolecular Motions

Life ticks as fast as how proteins move. Computationally expensive molecular dynamics simulation has been the only theoretical tool to gauge the time and sizes of these motions, though barely to their slowest ends. Here, we convert a computationally cheap elastic network model (ENM) into a molecular timer and sizer to gauge the slowest functional motions of structured biomolecules. Quasi-harmonic analysis, fluctuation profile matching, and the Wiener-Khintchine theorem are used to define the "time periods," t, for anharmonic principal components (PCs), which are validated by nuclear magnetic resonance (NMR) order parameters. The PCs with their respective "time periods" are mapped to the eigenvalues (λ_ENM) of the corresponding ENM modes. Thus, the power laws t(ns) = 56.1λ_ENM^-1.6 and σ²(Ų) = 32.7λ_ENM^-3.0 can be established allowing the characterization of the timescales of NMR-resolved conformers, crystallographic anisotropic displacement parameters, and important ribosomal motions, as well as motional sizes of the latter.

Large-scale simulations of nucleoprotein complexes: ribosomes, nucleosomes, chromatin, chromosomes and CRISPR

Recent advances in biotechnology such as Hi-C, CRISPR/Cas9 and ribosome display have placed nucleoprotein complexes at center stage. Understanding the structural dynamics of these complexes aids in optimizing protocols and interpreting data for these new technologies. The integration of simulation and experiment has helped advance mechanistic understanding of these systems. Coarse-grained simulations, reduced-description models, and explicit solvent molecular dynamics simulations yield useful complementary perspectives on nucleoprotein complex structural dynamics. When combined with Hi-C, cryo-EM, and single molecule measurements, these simulations integrate disparate forms of experimental data into a coherent mechanism.

Dissecting the Energetics of Subunit Rotation in the Ribosome

The accurate expression of proteins requires the ribosome to efficiently undergo elaborate conformational rearrangements. The most dramatic of these motions is subunit rotation, which is necessary for tRNA molecules to transition between ribosomal binding sites. While rigid-body descriptions provide a qualitative picture of the process, obtaining quantitative mechanistic insights requires one to account for the relationship between molecular flexibility and collective dynamics. Using simulated rotation events, we assess the quality of experimentally accessible measures for describing the collective displacement of the ∼4000-residue small subunit. For this, we ask whether each coordinate is able to identify the underlying free-energy barrier and transition state ensemble (TSE). We find that intuitive structurally motivated coordinates (e.g., rotation angle, interprotein distances) can distinguish between the endpoints, though they are poor indicators of barrier-crossing events, and they underestimate the free-energy barrier. In contrast, coordinates based on intersubunit bridges can identify the TSE. We additionally verify that the committor probability for the putative TSE configurations is 0.5, a hallmark feature of any transition state. In terms of structural properties, these calculations implicate a transition state in which flexibility allows for asynchronous rearrangements of the bridges, as the ribosome adopts a partially rotated orientation. This provides a theoretical foundation, upon which experimental techniques may precisely quantify the energy landscape of the ribosome.

Modeling of Protein Structural Flexibility and Large-Scale Dynamics: Coarse-Grained Simulations and Elastic Network Models

Fluctuations of protein three-dimensional structures and large-scale conformational transitions are crucial for the biological function of proteins and their complexes. Experimental studies of such phenomena remain very challenging and therefore molecular modeling can be a good alternative or a valuable supporting tool for the investigation of large molecular systems and long-time events. In this minireview, we present two alternative approaches to the coarse-grained (CG) modeling of dynamic properties of protein systems. We discuss two CG representations of polypeptide chains used for Monte Carlo dynamics simulations of protein local dynamics and conformational transitions, and highly simplified structure-based elastic network models of protein flexibility. In contrast to classical all-atom molecular dynamics, the modeling strategies discussed here allow the quite accurate modeling of much larger systems and longer-time dynamic phenomena. We briefly describe the main features of these models and outline some of their applications, including modeling of near-native structure fluctuations, sampling of large regions of the protein conformational space, or possible support for the structure prediction of large proteins and their complexes.

Dual DNA rulers reveal an 'mRNA looping' intermediate state during ribosome translocation

The precise 3-nucleotide movement of mRNA is critical for translation fidelity. One mRNA translocation error propagates to all of the following codons, which is detrimental to the cell. However, none of the current methods can reveal the mRNA dynamics near the ribosome entry site, which limits the understanding of this important issue. We have developed an assay of dual DNA rulers that provides such capability. By uniquely probing both the 3'- and 5'-ends of mRNA, we observed an antibiotic-trapped intermediate state that is consistent with a ribosomal conformation containing mRNA asymmetric partial displacements at its entry and exit sites. Based on the available ribosome structures and computational simulations, we proposed a 'looped' mRNA conformation, which suggested a stepwise 'inchworm' mechanism for ribosomal translocation. The same 'looped' intermediate state identified with the dual rulers persists with a '-1' frameshifting motif, indicating that the branching point of normal and frameshifting translocations occurs at a later stage of translocation.

Structures of Mycobacterium smegmatis 70S ribosomes in complex with HPF, tmRNA, and P-tRNA

Ribosomes are the dynamic protein synthesis machineries of the cell. They may exist in different functional states in the cell. Therefore, it is essential to have structural information on these different functional states of ribosomes to understand their mechanism of action. Here, we present single particle cryo-EM reconstructions of the Mycobacterium smegmatis 70S ribosomes in the hibernating state (with HPF), trans-translating state (with tmRNA), and the P/P state (with P-tRNA) resolved to 4.1, 12.5, and 3.4 Å, respectively. A comparison of the P/P state with the hibernating state provides possible functional insights about the Mycobacteria-specific helix H54a rRNA segment. Interestingly, densities for all the four OB domains of bS1 protein is visible in the hibernating 70S ribosome displaying the molecular details of bS1-70S interactions. Our structural data shows a Mycobacteria-specific H54a-bS1 interaction which seems to prevent subunit dissociation and degradation during hibernation without the formation of 100S dimer. This indicates a new role of bS1 protein in 70S protection during hibernation in Mycobacteria in addition to its conserved function during translation initiation.

Unique localization of the plastid-specific ribosomal proteins in the chloroplast ribosome small subunit provides mechanistic insights into the chloroplastic translation

Chloroplastic translation is mediated by a bacterial-type 70S chloroplast ribosome. During the evolution, chloroplast ribosomes have acquired five plastid-specific ribosomal proteins or PSRPs (cS22, cS23, bTHXc, cL37 and cL38) which have been suggested to play important regulatory roles in translation. However, their exact locations on the chloroplast ribosome remain elusive due to lack of a high-resolution structure, hindering our progress to understand their possible roles. Here we present a cryo-EM structure of the 70S chloroplast ribosome from spinach resolved to 3.4 Å and focus our discussion mainly on the architecture of the 30S small subunit (SSU) which is resolved to 3.7 Å. cS22 localizes at the SSU foot where it seems to compensate for the deletions in 16S rRNA. The mRNA exit site is highly remodeled due to the presence of cS23 suggesting an alternative mode of translation initiation. bTHXc is positioned at the SSU head and appears to stabilize the intersubunit bridge B1b during thermal fluctuations. The translation factor plastid pY binds to the SSU on the intersubunit side and interacts with the conserved nucleotide bases involved in decoding. Most of the intersubunit bridges are conserved compared to the bacteria, except for a new bridge involving uL2c and bS6c.

Identification of potential allosteric communication pathways between functional sites of the bacterial ribosome by graph and elastic network models

BACKGROUND

Accumulated evidence indicates that bacterial ribosome employs allostery throughout its structure for protein synthesis. The nature of the allosteric communication between remote functional sites remains unclear, but the contact topology and dynamics of residues may play role in transmission of a perturbation to distant sites.

METHODS/RESULTS

We employ two computationally efficient approaches - graph and elastic network modeling to gain insights about the allosteric communication in ribosome. Using graph representation of the structure, we perform k-shortest pathways analysis between peptidyl transferase center-ribosomal tunnel, decoding center-peptidyl transferase center - previously reported functional sites having allosteric communication. Detailed analysis on intact structures points to common and alternative shortest pathways preferred by different states of translation. All shortest pathways capture drug target sites and allosterically important regions. Elastic network model further reveals that residues along all pathways have the ability of quickly establishing pair-wise communication and to help the propagation of a perturbation in long-ranges during functional motions of the complex.

CONCLUSIONS

Contact topology and inherent dynamics of ribosome configure potential communication pathways between functional sites in different translation states. Inter-subunit bridges B2a, B3 and P-tRNA come forward for their high potential in assisting allostery during translation. Especially B3 emerges as a potential druggable site.

GENERAL SIGNIFICANCE

This study indicates that the ribosome topology forms a basis for allosteric communication, which can be disrupted by novel drugs to kill drug-resistant bacteria. Our computationally efficient approach not only overlaps with experimental evidence on allosteric regulation in ribosome but also proposes new druggable sites.

Distributions of experimental protein structures on coarse-grained free energy landscapes

Predicting conformational changes of proteins is needed in order to fully comprehend functional mechanisms. With the large number of available structures in sets of related proteins, it is now possible to directly visualize the clusters of conformations and their conformational transitions through the use of principal component analysis. The most striking observation about the distributions of the structures along the principal components is their highly non-uniform distributions. In this work, we use principal component analysis of experimental structures of 50 diverse proteins to extract the most important directions of their motions, sample structures along these directions, and estimate their free energy landscapes by combining knowledge-based potentials and entropy computed from elastic network models. When these resulting motions are visualized upon their coarse-grained free energy landscapes, the basis for conformational pathways becomes readily apparent. Using three well-studied proteins, T4 lysozyme, serum albumin, and sarco-endoplasmic reticular Ca(2+) adenosine triphosphatase (SERCA), as examples, we show that such free energy landscapes of conformational changes provide meaningful insights into the functional dynamics and suggest transition pathways between different conformational states. As a further example, we also show that Monte Carlo simulations on the coarse-grained landscape of HIV-1 protease can directly yield pathways for force-driven conformational changes.

Ribosome Mechanics Informs about Mechanism

The essential aspects of the ribosome's mechanism can be extracted from coarse-grained simulations, including the ratchet motion, the movement together of critical bases at the decoding center, and movements of the peptide tunnel lining that assist in the expulsion of the synthesized peptide. Because of its large size, coarse graining helps to simplify and to aid in the understanding of its mechanism. Results presented here utilize coarse-grained elastic network modeling to extract the dynamics, and both RNAs and proteins are coarse grained. We review our previous results, showing the well-known ratchet motions and the motions in the peptide tunnel and in the mRNA tunnel. The motions of the lining of the peptide tunnel appear to assist in the expulsion of the growing peptide chain, and clamps at the ends of the mRNA tunnel with three proteins ensure that the mRNA is held tightly during decoding and essential for the helicase activity at the entrance. The entry clamp may also assist in base recognition to ensure proper selection of the incoming tRNA. The overall precision of the ribosome machine-like motions is remarkable.

Bridging between NMA and Elastic Network Models: Preserving All-Atom Accuracy in Coarse-Grained Models

Dynamics can provide deep insights into the functional mechanisms of proteins and protein complexes. For large protein complexes such as GroEL/GroES with more than 8,000 residues, obtaining a fine-grained all-atom description of its normal mode motions can be computationally prohibitive and is often unnecessary. For this reason, coarse-grained models have been used successfully. However, most existing coarse-grained models use extremely simple potentials to represent the interactions within the coarse-grained structures and as a result, the dynamics obtained for the coarse-grained structures may not always be fully realistic. There is a gap between the quality of the dynamics of the coarse-grained structures given by all-atom models and that by coarse-grained models. In this work, we resolve an important question in protein dynamics computations—how can we efficiently construct coarse-grained models whose description of the dynamics of the coarse-grained structures remains as accurate as that given by all-atom models? Our method takes advantage of the sparseness of the Hessian matrix and achieves a high efficiency with a novel iterative matrix projection approach. The result is highly significant since it can provide descriptions of normal mode motions at an all-atom level of accuracy even for the largest biomolecular complexes. The application of our method to GroEL/GroES offers new insights into the mechanism of this biologically important chaperonin, such as that the conformational transitions of this protein complex in its functional cycle are even more strongly connected to the first few lowest frequency modes than with other coarse-grained models.

Proteins and other biomolecules are not static but are constantly in motion. Moreover, they possess intrinsic collective motion patterns that are tightly linked to their functions. Thus, an accurate and detailed description of their motions can provide deep insights into their functional mechanisms. For large protein complexes with hundreds of thousands of atoms or more, an atomic level description of the motions can be computationally prohibitive, and so coarse-grained models with fewer structural details are often used instead. However, there can be a big gap between the quality of motions derived from atomic models and those from coarse-grained models. In this work, we solve an important problem in protein dynamics studies: how to preserve the atomic-level accuracy in describing molecular motions while using coarse-grained models? We accomplish this by developing a novel iterative matrix projection method that dramatically speeds up the computations. This method is significant since it promises accurate descriptions of protein motions approaching an all-atom level even for the largest biomolecular complexes. Results shown here for a large molecular chaperonin demonstrate how this can provide new insights into its functional process.

Functional Importance of Mobile Ribosomal Proteins

Although the dynamic motions and peptidyl transferase activity seem to be embedded in the rRNAs, the ribosome contains more than 50 ribosomal proteins (r-proteins), whose functions remain largely elusive. Also, the precise forms of some of these r-proteins, as being part of the ribosome, are not structurally solved due to their high flexibility, which hinders the efforts in their functional elucidation. Owing to recent advances in cryo-electron microscopy, single-molecule techniques, and theoretical modeling, much has been learned about the dynamics of these r-proteins. Surprisingly, allosteric regulations have been found in between spatially separated components as distant as those in the opposite sides of the ribosome. Here, we focus on the functional roles and intricate regulations of the mobile L1 and L12 stalks and L9 and S1 proteins. Conformational flexibility also enables versatile functions for r-proteins beyond translation. The arrangement of r-proteins may be under evolutionary pressure that fine-tunes mass distributions for optimal structural dynamics and catalytic activity of the ribosome.

The ribosome's energy landscape: Recent insights from computation

The ever-increasing capacity of computing resources has extended ribosome calculations from the study of small-scale fluctuations to large-scale barrier-crossing processes. As the field of computational/theoretical biophysics shifts focus to large-scale conformational transitions, there is a growing need for a systematic framework to interpret and analyze ribosome dynamics. To this end, energy landscape principles, largely developed for the study of biomolecular folding, have proven to be invaluable. These tools not only provide a foundation for describing simulations but can be used to reconcile experimental results, as well. In this review, I will discuss recent efforts to employ computational methods to reveal the characteristics of the ribosome's landscape and how these studies can help guide a new generation of experiments that more closely probe the underlying energetics. As a result of these investigations, general principles about ribosome function are beginning to emerge, including that: (1) small-scale fluctuations are the result of structure, rather than detailed energetics, (2) molecular flexibility leads to entropically favored rearrangements, and (3) tRNA dynamics may be accurately described as diffusive movement across an energy landscape.

Using Molecular Simulation to Model High-Resolution Cryo-EM Reconstructions

An explosion of new data from high-resolution cryo-electron microscopy (cryo-EM) studies has produced a large number of data sets for many species of ribosomes in various functional states over the past few years. While many methods exist to produce structural models for lower resolution cryo-EM reconstructions, high-resolution reconstructions are often modeled using crystallographic techniques and extensive manual intervention. Here, we present an automated fitting technique for high-resolution cryo-EM data sets that produces all-atom models highly consistent with the EM density. Using a molecular dynamics approach, atomic positions are optimized with a potential that includes the cross-correlation coefficient between the structural model and the cryo-EM electron density, as well as a biasing potential preserving the stereochemistry and secondary structure of the biomolecule. Specifically, we use a hybrid structure-based/ab initio molecular dynamics potential to extend molecular dynamics fitting. In addition, we find that simulated annealing integration, as opposed to straightforward molecular dynamics integration, significantly improves performance. We obtain atomistic models of the human ribosome consistent with high-resolution cryo-EM reconstructions of the human ribosome. Automated methods such as these have the potential to produce atomistic models for a large number of ribosome complexes simultaneously that can be subsequently refined manually.

RNA-protein distance patterns in ribosomes reveal the mechanism of translational attenuation

Elucidating protein translational regulation is crucial for understanding cellular function and drug development. A key molecule in protein translation is ribosome, which is a super-molecular complex extensively studied for more than a half century. The structure and dynamics of ribosome complexes were resolved recently thanks to the development of X-ray crystallography, Cryo-EM, and single molecule biophysics. Current studies of the ribosome have shown multiple functional states, each with a unique conformation. In this study, we analyzed the RNA-protein distances of ribosome (2.5 MDa) complexes and compared these changes among different ribosome complexes. We found that the RNA-protein distance is significantly correlated with the ribosomal functional state. Thus, the analysis of RNA-protein binding distances at important functional sites can distinguish ribosomal functional states and help understand ribosome functions. In particular, the mechanism of translational attenuation by nascent peptides and antibiotics was revealed by the conformational changes of local functional sites.

Connecting the kinetics and energy landscape of tRNA translocation on the ribosome

Functional rearrangements in biomolecular assemblies result from diffusion across an underlying energy landscape. While bulk kinetic measurements rely on discrete state-like approximations to the energy landscape, single-molecule methods can project the free energy onto specific coordinates. With measures of the diffusion, one may establish a quantitative bridge between state-like kinetic measurements and the continuous energy landscape. We used an all-atom molecular dynamics simulation of the 70S ribosome (2.1 million atoms; 1.3 microseconds) to provide this bridge for specific conformational events associated with the process of tRNA translocation. Starting from a pre-translocation configuration, we identified sets of residues that collectively undergo rotary rearrangements implicated in ribosome function. Estimates of the diffusion coefficients along these collective coordinates for translocation were then used to interconvert between experimental rates and measures of the energy landscape. This analysis, in conjunction with previously reported experimental rates of translocation, provides an upper-bound estimate of the free-energy barriers associated with translocation. While this analysis was performed for a particular kinetic scheme of translocation, the quantitative framework is general and may be applied to energetic and kinetic descriptions that include any number of intermediates and transition states.

Elastic network normal modes provide a basis for protein structure refinement

It is well recognized that thermal motions of atoms in the protein native state, the fluctuations about the minimum of the global free energy, are well reproduced by the simple elastic network models (ENMs) such as the anisotropic network model (ANM). Elastic network models represent protein dynamics as vibrations of a network of nodes (usually represented by positions of the heavy atoms or by the C(α) atoms only for coarse-grained representations) in which the spatially close nodes are connected by harmonic springs. These models provide a reliable representation of the fluctuational dynamics of proteins and RNA, and explain various conformational changes in protein structures including those important for ligand binding. In the present paper, we study the problem of protein structure refinement by analyzing thermal motions of proteins in non-native states. We represent the conformational space close to the native state by a set of decoys generated by the I-TASSER protein structure prediction server utilizing template-free modeling. The protein substates are selected by hierarchical structure clustering. The main finding is that thermal motions for some substates, overlap significantly with the deformations necessary to reach the native state. Additionally, more mobile residues yield higher overlaps with the required deformations than do the less mobile ones. These findings suggest that structural refinement of poorly resolved protein models can be significantly enhanced by reduction of the conformational space to the motions imposed by the dominant normal modes.

Dynamics of endoglucanase catalytic domains: implications towards thermostability

Thermostable endoglucanases play a crucial role in the production of biofuels to breakdown plant cellulose. Analyzing their structure-dynamics relationship can inform about the origins of their thermostability. Although tertiary structures of many endoglucanase proteins are available, the relationship between thermostability, structure, and dynamics is not explored fully. We have generated elastic network models for thermostable and mesostable endoglucanases with the (αβ)₈ fold in substrate bound and unbound states. The comparative analyses shed light on the relation between protein dynamics, thermostability, and substrate binding. We observed specific differences in the dynamic behavior of catalytic residues in slow modes: while both the nucleophile and the acid/base donor residues show positively correlated motions in the thermophile, their dynamics is uncoupled in the mesophile. Our proof-of-concept comparison study suggests that global dynamics can be harnessed to further our understanding of thermostability.

Mechanisms of SecM-mediated stalling in the ribosome

Nascent-peptide modulation of translation is a common regulatory mechanism of gene expression. In this mechanism, while the nascent peptide is still in the exit tunnel of the ribosome, it induces translational pausing, thereby controlling the expression of downstream genes. One example is SecM, which inhibits peptide-bond formation in the ribosome's peptidyl transferase center (PTC) during its own translation, upregulating the expression of the protein translocase SecA. Although biochemical experiments and cryo-electron microscopy data have led to the identification of some residues involved in SecM recognition, the full pathway of interacting residues that connect SecM to the PTC through the ribosome has not yet been conclusively established. Here, using the cryo-electron microscopy data, we derived the first (to our knowledge) atomic model of the SecM-stalled ribosome via molecular-dynamics flexible fitting, complete with P- and A-site tRNAs. Subsequently, we carried out simulations of native and mutated SecM-stalled ribosomes to investigate possible interaction pathways between a critical SecM residue, R163, and the PTC. In particular, the simulations reveal the role of SecM in altering the position of the tRNAs in the ribosome, and thus demonstrate how the presence of SecM in the exit tunnel induces stalling. Finally, steered molecular-dynamics simulations in which SecM was pulled toward the tunnel exit suggest how SecA interacting with SecM from outside the ribosome relieves stalling.

Biomolecular dynamics: order-disorder transitions and energy landscapes

While the energy landscape theory of protein folding is now a widely accepted view for understanding how relatively weak molecular interactions lead to rapid and cooperative protein folding, such a framework must be extended to describe the large-scale functional motions observed in molecular machines. In this review, we discuss (1) the development of the energy landscape theory of biomolecular folding, (2) recent advances toward establishing a consistent understanding of folding and function and (3) emerging themes in the functional motions of enzymes, biomolecular motors and other biomolecular machines. Recent theoretical, computational and experimental lines of investigation have provided a very dynamic picture of biomolecular motion. In contrast to earlier ideas, where molecular machines were thought to function similarly to macroscopic machines, with rigid components that move along a few degrees of freedom in a deterministic fashion, biomolecular complexes are only marginally stable. Since the stabilizing contribution of each atomic interaction is on the order of the thermal fluctuations in solution, the rigid body description of molecular function must be revisited. An emerging theme is that functional motions encompass order-disorder transitions and structural flexibility provides significant contributions to the free energy. In this review, we describe the biological importance of order-disorder transitions and discuss the statistical-mechanical foundation of theoretical approaches that can characterize such transitions.

A computational investigation on the connection between dynamics properties of ribosomal proteins and ribosome assembly

Assembly of the ribosome from its protein and RNA constituents has been studied extensively over the past 50 years, and experimental evidence suggests that prokaryotic ribosomal proteins undergo conformational changes during assembly. However, to date, no studies have attempted to elucidate these conformational changes. The present work utilizes computational methods to analyze protein dynamics and to investigate the linkage between dynamics and binding of these proteins during the assembly of the ribosome. Ribosomal proteins are known to be positively charged and we find the percentage of positive residues in r-proteins to be about twice that of the average protein: Lys+Arg is 18.7% for E. coli and 21.2% for T. thermophilus. Also, positive residues constitute a large proportion of RNA contacting residues: 39% for E. coli and 46% for T. thermophilus. This affirms the known importance of charge-charge interactions in the assembly of the ribosome. We studied the dynamics of three primary proteins from E. coli and T. thermophilus 30S subunits that bind early in the assembly (S15, S17, and S20) with atomic molecular dynamic simulations, followed by a study of all r-proteins using elastic network models. Molecular dynamics simulations show that solvent-exposed proteins (S15 and S17) tend to adopt more stable solution conformations than an RNA-embedded protein (S20). We also find protein residues that contact the 16S rRNA are generally more mobile in comparison with the other residues. This is because there is a larger proportion of contacting residues located in flexible loop regions. By the use of elastic network models, which are computationally more efficient, we show that this trend holds for most of the 30S r-proteins.

Ribosomes are complex cellular machines that synthesize new proteins in the cell. The accurate and efficient assembly of ribosomal proteins (r-proteins) and ribosomal RNA (rRNA) to form a functional ribosome is important for cell growth, metabolic reactions, and other cellular processes. Additionally, some antibacterial drugs are believed to target the bacterial ribosome during its construction. Hence, ribosomal assembly has been an active research topic for many years because understanding the assembly mechanisms can provide insight into protein/RNA recognitions important in many other cellular processes, as well as optimize the development of antibacterial therapeutics. Experimental studies thus far have provided still limited understanding about the assembly process. To further understand the assembly process, we have computationally studied the dynamic properties that r-proteins exhibit during assembly and the relationship between dynamics, physical properties, and binding propensity. We observe significant charged interactions between r-proteins and rRNA. We also detect a strong correlation between contact residues and their dynamic mobilities. Protein residues contacting with rRNA are observed to be more mobile in comparison with other residues. We also relate the location of the r-protein in the fully assembled ribosome to its susceptibility for large conformational changes prior to binding.

Revealing -1 programmed ribosomal frameshifting mechanisms by single-molecule techniques and computational methods

Programmed ribosomal frameshifting (PRF) serves as an intrinsic translational regulation mechanism employed by some viruses to control the ratio between structural and enzymatic proteins. Most viral mRNAs which use PRF adapt an H-type pseudoknot to stimulate −1 PRF. The relationship between the thermodynamic stability and the frameshifting efficiency of pseudoknots has not been fully understood. Recently, single-molecule force spectroscopy has revealed that the frequency of −1 PRF correlates with the unwinding forces required for disrupting pseudoknots, and that some of the unwinding work dissipates irreversibly due to the torsional restraint of pseudoknots. Complementary to single-molecule techniques, computational modeling provides insights into global motions of the ribosome, whose structural transitions during frameshifting have not yet been elucidated in atomic detail. Taken together, recent advances in biophysical tools may help to develop antiviral therapies that target the ubiquitous −1 PRF mechanism among viruses.

Computational studies of molecular machines: the ribosome

The past decade has produced an avalanche of experimental data on the structure and dynamics of the ribosome. Groundbreaking studies in structural biology and kinetics have placed important constraints on ribosome structural dynamics. However, a gulf remains between static structures and time dependent data. In particular, X-ray crystallography and cryo-EM studies produce static models of the ribosome in various states, but lack dynamic information. Single molecule studies produce information on the rates of transitions between these states but do not have high-resolution spatial information. Computational studies have aided in bridging this gap by providing atomic resolution simulations of structural fluctuations and transitions between configurations.

Key intermolecular interactions in the E. coli 70S ribosome revealed by coarse-grained analysis

The ribosome is a very large complex that consists of many RNA and protein molecules and plays a central role in protein biosynthesis in all organisms. Extensive interactions between different molecules are critical to ribosomal functional dynamics. In this work, intermolecular interactions in the Escherichia coli 70S ribosome are investigated by coarse-grained (CG) analysis. CG models are defined to preserve dynamic domains in RNAs and proteins and to capture functional motions in the ribosome, and then the CG sites are connected by harmonic springs, and spring constants are obtained by matching the computed fluctuations to those of an all-atom molecular dynamics (MD) simulation. Those spring constants indicate how strong the interactions are between the ribosomal components, and they are in good agreement with various experimental data. Nearly all the bridges between the small and large ribosomal subunits are indicated by CG interactions with large spring constants. The head of the small subunit is very mobile because it has minimal CG interactions with the rest of the subunit; however, a large number of small subunit proteins bind to maintain the internal structure of the head. The results show a clear connection between the intermolecular interactions and the structural and functional properties of the ribosome because of the reduced complexity in domain-based CG models. The present approach also provides a useful strategy to map interactions between molecules within large biomolecular complexes since it is not straightforward to investigate these by either atomistic MD simulations or residue-based elastic network models.

Brownian dynamics study of the association between the 70S ribosome and elongation factor G

Protein synthesis on the ribosome involves a number of external protein factors that bind at its functional sites. One key factor is the elongation factor G (EF-G) that facilitates the translocation of transfer RNAs between their binding sites, as well as advancement of the messenger RNA by one codon. The details of the EF-G/ribosome diffusional encounter and EF-G association pathway still remain unanswered. Here, we applied Brownian dynamics methodology to study bimolecular association in the bacterial EF-G/70S ribosome system. We estimated the EF-G association rate constants at 150 and 300 mM monovalent ionic strengths and obtained reasonable agreement with kinetic experiments. We have also elucidated the details of EF-G/ribosome association paths and found that positioning of the L11 protein of the large ribosomal subunit is likely crucial for EF-G entry to its binding site.

Phosphorylation of ribosomal proteins influences subunit association and translation of poly (U) in Streptomyces coelicolor

The occurrence of phosphorylated proteins in ribosomes of Streptomyces coelicolor was investigated. Little is known about which biological functions these posttranslational modifications might fulfil. A protein kinase associated with ribosomes phosphorylated six ribosomal proteins of the small subunit (S3, S4, S12, S13, S14 and S18) and seven ribosomal proteins of the large subunit (L2, L3, L7/L12, L16, L17, L23 and L27). The ribosomal proteins were phosphorylated mainly on the Ser/Thr residues. Phosphorylation of the ribosomal proteins influences ribosomal subunits association. Ribosomes with phosphorylated proteins were used to examine poly (U) translation activity. Phosphorylation induced about 50% decrease in polyphenylalanine synthesis. After preincubation of ribosomes with alkaline phosphatase the activity of ribosomes was greatly restored. Small differences were observed between phosphorylated and unphosphorylated ribosomes in the kinetic parameters of the binding of Phe-tRNA to the A-site of poly (U) programmed ribosomes, suggesting that the initial binding of Phe-tRNA is not significantly affected by phosphorylation. On contrary, the rate of peptidyl transferase was about two-fold lower than that in unphosphorylated ribosomes. The data presented demonstrate that phosphorylation of ribosomal proteins affects critical steps of protein synthesis.

Coarse-grained models to study dynamics of nanoscale biomolecules and their applications to the ribosome

Biopolymers are of dynamic nature and undergo functional motions spanning a large spectrum of timescales. To study the internal dynamics of nano-sized molecular complexes that exceed hundred thousands of atoms with atomic detail is computationally inefficient. Therefore, to achieve both the spatial and temporal scales of biological interest coarse-grained models of macromolecules are often used. By uniting groups of atoms into single interacting centers one decreases the resolution of the system and gets rid of the irrelevant degrees of freedom. This simplification, even though it requires parameterization, makes the studies of biomolecular dynamics computationally tractable and allows us to reach beyond the microsecond time frame. Here, I review the coarse-grained models of macromolecules composed of proteins and nucleic acids. I give examples of one-bead models that were developed to investigate the internal dynamics and focus on their applications to the ribosome--the nanoscale protein synthesis machine.

Models to Approximate the Motions of Protein Loops

We approximate the loop motions of various proteins by using a coarse-grained model and the theory of rubberlike elasticity of polymer chains. The loops are considered as chains where only the first and the last residues thereof are tethered by their connections to the main structure; while within the loop, the loop residues are connected only to their sequence neighbors. We applied these approximate models to five proteins. Our approximation shows that the loop motions can usually be computed locally which shows these motions are robust and not random. But most interestingly, the new method presented here can be used to compute the likely motions of loops that are missing in the structures.

Coarse-Grained Representations of Large Biomolecular Complexes from Low-Resolution Structural Data

High-resolution atomistic structures of many large biomolecular complexes have not yet been solved by experiments, such as X-ray crystallography or NMR. Often however low-resolution information is obtained by alternative techniques, such as cryo-electron microscopy or small-angle X-ray scattering. Coarse-grained (CG) models are an appropriate choice to computationally study these complexes given the limited resolution experimental data. One of the important questions therefore is how to define CG representations from these low-resolution density maps. This work provides a space-based essential dynamics coarse-graining (ED-CG) method to define a CG representation from a density map without detailed knowledge of its underlying atomistic structure and primary sequence information. This method is demonstrated on G-actin (both the atomic structure and its density map). It is then applied to the density maps of the Escherichia coli 70S ribosome and the microtubule. The results indicate that the method can define highly CG models that still preserve functionally important dynamics of large biomolecular complexes.

Mechanism of cohesin loading onto chromosomes: a conformational dynamics study

The structure-function relationship of cohesin, an essential chromosome maintenance protein, is investigated by analyzing its collective dynamics and conformational flexibility, enhancing our understanding of the sister chromatid cohesion process. A three-dimensional model of cohesin has been constructed by homology modeling using both crystallographic and electron microscopy image data. The harmonic dynamics of the cohesin structure are calculated with a coarse-grained elastic network model. The model shows that the bending motion of the cohesin ring is able to adopt a head-to-tail conformation, in agreement with experimental data. Low-frequency conformational changes are observed to deform the highly conserved glycine residues at the interface of the cohesin heterodimer. Normal mode analysis further reveals that, near large globular structures such as nucleosome and accessory proteins docked to cohesin, the mobility of the coiled-coil regions is notably affected. Moreover, fully solvated molecular dynamics calculations, performed specifically on the hinge region, indicate that hinge opening starts from one side of the dimerization interface, and is coordinated by highly conserved glycine residues.