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

Embracing exascale computing in nucleic acid simulations

This mini-review reports the recent advances in biomolecular simulations, particularly for nucleic acids, and provides the potential effects of the emerging exascale computing on nucleic acid simulations, emphasizing the need for advanced computational strategies to fully exploit this technological frontier. Specifically, we introduce recent breakthroughs in computer architectures for large-scale biomolecular simulations and review the simulation protocols for nucleic acids regarding force fields, enhanced sampling methods, coarse-grained models, and interactions with ligands. We also explore the integration of machine learning methods into simulations, which promises to significantly enhance the predictive modeling of biomolecules and the analysis of complex data generated by the exascale simulations. Finally, we discuss the challenges and perspectives for biomolecular simulations as we enter the dawning exascale computing era.

Structural dynamics in chromatin unraveling by pioneer transcription factors

Pioneer transcription factors are proteins with a dual function. First, they regulate transcription by binding to nucleosome-free DNA regulatory elements. Second, they bind to DNA while wrapped around histone proteins in the chromatin and mediate chromatin opening. The molecular mechanisms that connect the two functions are yet to be discovered. In recent years, pioneer factors received increased attention mainly because of their crucial role in promoting cell fate transitions that could be used for regenerative therapies. For example, the three factors required to induce pluripotency in somatic cells, Oct4, Sox2, and Klf4 were classified as pioneer factors and studied extensively. With this increased attention, several structures of complexes between pioneer factors and chromatin structural units (nucleosomes) have been resolved experimentally. Furthermore, experimental and computational approaches have been designed to study two unresolved, key scientific questions: First, do pioneer factors induce directly local opening of nucleosomes and chromatin fibers upon binding? And second, how do the unstructured tails of the histones impact the structural dynamics involved in such conformational transitions? Here we review the current knowledge about transcription factor-induced nucleosome dynamics and the role of the histone tails in this process. We discuss what is needed to bridge the gap between the static views obtained from the experimental structures and the key structural dynamic events in chromatin opening. Finally, we propose that integrating nuclear magnetic resonance spectroscopy with molecular dynamics simulations is a powerful approach to studying pioneer factor-mediated dynamics of nucleosomes and perhaps small chromatin fibers using native DNA sequences.

Conformational changes of ribosomes during translation elongation resolved by molecular dynamics simulations

Molecular dynamics simulations have emerged as a powerful set of tools to unravel the intricate dynamics of ribosomes during protein synthesis. Recent advancements in this field have enabled simulations to delve deep into the conformational rearrangements of ribosomes and associated factors, providing invaluable insights into the intricacies of translation. Emphasis on simulations has recently been on translation elongation, such as tRNA selection, translocation, and ribosomal head-swivel motions. These studies have offered crucial structural interpretations of how genetic information is faithfully translated into proteins. This review outlines recent discoveries concerning ribosome conformational changes occurring during translation elongation, as elucidated through molecular dynamics simulations.

Structure and interaction of therapeutic proteins in solution: a combined simulation and experimental study

The aggregation of therapeutic proteins in solution has attracted significant interest, driving efforts to understand the relationship between microscopic structural changes and protein-protein interactions determining aggregation processes in solution. Additionally, there is substantial interest in being able to predict aggregation based on protein structure as part of molecular developability assessments. Molecular Dynamics provides theoretical tools to complement experimental studies and to interrogate and identify the microscopic mechanisms determining aggregation. Here we perform all-atom MD simulations to study the structure and inter-protein interaction of the Fab and Fc fragments of the monoclonal antibody (mAb) COE3. We unravel the role of ion-protein interactions in building the ionic double layer and determining effective inter-protein interaction. Further, we demonstrate, using various state-of-the-art force fields (charmm, gromos, amber, opls/aa), that the protein solvation, ionic structure and protein-protein interaction depend significantly on the force field parameters. We perform SANS and Static Light Scattering experiments to assess the accuracy of the different forcefields. Comparison of the simulated and experimental results reveal significant differences in the forcefields' performance, particularly in their ability to predict the protein size in solution and inter-protein interactions quantified through the second virial coefficients. In addition, the performance of the forcefields is correlated with the protein hydration structure.

Dynamic exchange controls the assembly structure of nucleic-acid-peptide chimeras

Recent attempts to develop the next generation of functional biomaterials focus on systems chemistry approaches exploiting dynamic networks of hybrid molecules. This task is often found challenging, but we herein present ways for profiting from the multiple interaction interfaces forming Nucleic-acid-Peptide assemblies and tuning their formation. We demonstrate that the formation of well-defined structures by double-stranded DNA-peptide conjugates (dsCon) is restricted to a specific range of environmental conditions and that precise DNA hybridization, satisfying the interaction interfaces, is a crucial factor in this process. We further reveal the impact of external stimuli, such as competing free DNA elements or salt additives, which initiate dynamic interconversions, resulting in hybrid structures exhibiting spherical and fibrillar domains or a mixture of spherical and fibrillar particles. This extensive analysis of the co-assembly systems chemistry offers new insights into prebiotic hybrid assemblies that may now facilitate the design of new functional materials. We discuss the implications of these findings for the emergence of function in synthetic materials and during early chemical evolution.

Decoding and Recoding of mRNA Sequences by the Ribosome

Faithful translation of messenger RNA (mRNA) into protein is essential to maintain protein homeostasis in the cell. Spontaneous translation errors are very rare due to stringent selection of cognate aminoacyl transfer RNAs (tRNAs) and the tight control of the mRNA reading frame by the ribosome. Recoding events, such as stop codon readthrough, frameshifting, and translational bypassing, reprogram the ribosome to make intentional mistakes and produce alternative proteins from the same mRNA. The hallmark of recoding is the change of ribosome dynamics. The signals for recoding are built into the mRNA, but their reading depends on the genetic makeup of the cell, resulting in cell-specific changes in expression programs. In this review, I discuss the mechanisms of canonical decoding and tRNA-mRNA translocation; describe alternative pathways leading to recoding; and identify the links among mRNA signals, ribosome dynamics, and recoding.

Engineering and Design of Programmable Genome Editors

Programmable genome editors are enzymes that can be targeted to a specific location in the genome for making site-specific alterations or deletions. The engineering, design, and development of sequence-specific editors has resulted in a dramatic increase in the precision of editing for nucleotide sequences. These editors can target specific locations in a genome, in vivo. The genome editors are being deployed for the development of genetically modified organisms for agriculture and industry, and for gene therapy of inherited human genetic disorders, cancer, and immunotherapy. Experimental and computational studies of structure, binding, activity, dynamics, and folding, reviewed here, have provided valuable insights that have the potential for increasing the functional efficiency of these gene/genome editors. Biochemical and biophysical studies of the specificities of natural and engineered genome editors reveal that increased binding affinity can be detrimental because of the increase of off-target effects and that the engineering and design of genome editors with higher specificity may require modulation and control of the conformational dynamics.

Hyper-swivel head domain motions are required for complete mRNA-tRNA translocation and ribosome resetting

Translocation of messenger RNA (mRNA) and transfer RNA (tRNA) substrates through the ribosome during protein synthesis, an exemplar of directional molecular movement in biology, entails a complex interplay of conformational, compositional, and chemical changes. The molecular determinants of early translocation steps have been investigated rigorously. However, the elements enabling the ribosome to complete translocation and reset for subsequent protein synthesis reactions remain poorly understood. Here, we have combined molecular simulations with single-molecule fluorescence resonance energy transfer imaging to gain insights into the rate-limiting events of the translocation mechanism. We find that diffusive motions of the ribosomal small subunit head domain to hyper-swivelled positions, governed by universally conserved rRNA, can maneuver the mRNA and tRNAs to their fully translocated positions. Subsequent engagement of peptidyl-tRNA and disengagement of deacyl-tRNA from mRNA, within their respective small subunit binding sites, facilitate the ribosome resetting mechanism after translocation has occurred to enable protein synthesis to resume.

Collective mechanical responses of cadherin-based adhesive junctions as predicted by simulations

Cadherin-based adherens junctions and desmosomes help stabilize cell-cell contacts with additional function in mechano-signaling, while clustered protocadherin junctions are responsible for directing neuronal circuits assembly. Structural models for adherens junctions formed by epithelial cadherin (CDH1) proteins indicate that their long, curved ectodomains arrange to form a periodic, two-dimensional lattice stabilized by tip-to-tip trans interactions (across junction) and lateral cis contacts. Less is known about the exact architecture of desmosomes, but desmoglein (DSG) and desmocollin (DSC) cadherin proteins are also thought to form ordered junctions. In contrast, clustered protocadherin (PCDH)-based cell-cell contacts in neuronal tissues are thought to be responsible for self-recognition and avoidance, and structural models for clustered PCDH junctions show a linear arrangement in which their long and straight ectodomains form antiparallel overlapped trans complexes. Here, we report all-atom molecular dynamics simulations testing the mechanics of minimalistic adhesive junctions formed by CDH1, DSG2 coupled to DSC1, and PCDHγB4, with systems encompassing up to 3.7 million atoms. Simulations generally predict a favored shearing pathway for the adherens junction model and a two-phased elastic response to tensile forces for the adhesive adherens junction and the desmosome models. Complexes within these junctions first unbend at low tensile force and then become stiff to unbind without unfolding. However, cis interactions in both the CDH1 and DSG2-DSC1 systems dictate varied mechanical responses of individual dimers within the junctions. Conversely, the clustered protocadherin PCDHγB4 junction lacks a distinct two-phased elastic response. Instead, applied tensile force strains trans interactions directly, as there is little unbending of monomers within the junction. Transient intermediates, influenced by new cis interactions, are observed after the main rupture event. We suggest that these collective, complex mechanical responses mediated by cis contacts facilitate distinct functions in robust cell-cell adhesion for classical cadherins and in self-avoidance signaling for clustered PCDHs.

Searching for Low Probability Opening Events in a DNA Sliding Clamp

The β subunit of E. coli DNA polymererase III is a DNA sliding clamp associated with increasing the processivity of DNA synthesis. In its free form, it is a circular homodimer structure that can accomodate double-stranded DNA in a nonspecific manner. An open state of the clamp must be accessible before loading the DNA. The opening mechanism is still a matter of debate, as is the effect of bound DNA on opening/closing kinetics. We use a combination of atomistic, coarse-grained, and enhanced sampling strategies in both explicit and implicit solvents to identify opening events in the sliding clamp. Such simulations of large nucleic acid and their complexes are becoming available and are being driven by improvements in force fields and the creation of faster computers. Different models support alternative opening mechanisms, either through an in-plane or out-of-plane opening event. We further note some of the current limitations, despite advances, in modeling these highly charged systems with implicit solvent.

Mechanism for the Unfolding of the TOP7 Protein in Steered Molecular Dynamics Simulations as Revealed by Mutual Information Analysis

We employed mutual information (MI) analysis to detect motions affecting the mechanical resistance of the human-engineered protein Top7. The results are based on the MI analysis of pair contact correlations measured in steered molecular dynamics (SMD) trajectories and their statistical dependence on global unfolding. This study is the first application of the MI analysis to SMD forced unfolding, and we furnish specific SMD recommendations for the utility of parameters and options in the TimeScapes package. The MI analysis provided a global overview of the effect of perturbation on the stability of the protein. We also employed a more conventional trajectory analysis for a detailed description of the mechanical resistance of Top7. Specifically, we investigated 1) the hydropathy of the interactions of structural segments, 2) the H₂O concentration near residues relevant for unfolding, and 3) the changing hydrogen bonding patterns and main chain dihedral angles. The results show that the application of MI in the study of protein mechanical resistance can be useful for the engineering of more resistant mutants when combined with conventional analysis. We propose a novel mutation design based on the hydropathy of residues that would stabilize the unfolding region by mimicking its more stable symmetry mate. The proposed design process does not involve the introduction of covalent crosslinks, so it has the potential to preserve the conformational space and unfolding pathway of the protein.

Getting to the bottom of lncRNA mechanism: structure-function relationships

While long non-coding RNAs are known to play key roles in disease and development, relatively few structural studies have been performed for this important class of RNAs. Here, we review functional studies of long non-coding RNAs and expose the need for high-resolution 3-D structural studies, discussing the roles of long non-coding RNAs in the cell and how structure–function relationships might be used to elucidate further understanding. We then describe structural studies of other classes of RNAs using chemical probing, nuclear magnetic resonance, small-angle X-ray scattering, X-ray crystallography, and cryogenic electron microscopy (cryo-EM). Next, we review early structural studies of long non-coding RNAs to date and describe the way forward for the structural biology of long non-coding RNAs in terms of cryo-EM.

Primitive selection of the fittest emerging through functional synergy in nucleopeptide networks

Many fundamental cellular and viral functions, including replication and translation, involve complex ensembles hosting synergistic activity between nucleic acids and proteins/peptides. There is ample evidence indicating that the chemical precursors of both nucleic acids and peptides could be efficiently formed in the prebiotic environment. Yet, studies on nonenzymatic replication, a central mechanism driving early chemical evolution, have focused largely on the activity of each class of these molecules separately. We show here that short nucleopeptide chimeras can replicate through autocatalytic and cross-catalytic processes, governed synergistically by the hybridization of the nucleobase motifs and the assembly propensity of the peptide segments. Unequal assembly-dependent replication induces clear selectivity toward the formation of a certain species within small networks of complementary nucleopeptides. The selectivity pattern may be influenced and indeed maximized to the point of almost extinction of the weakest replicator when the system is studied far from equilibrium and manipulated through changes in the physical (flow) and chemical (template and inhibition) conditions. We postulate that similar processes may have led to the emergence of the first functional nucleic-acid-peptide assemblies prior to the origin of life. Furthermore, spontaneous formation of related replicating complexes could potentially mark the initiation point for information transfer and rapid progression in complexity within primitive environments, which would have facilitated the development of a variety of functions found in extant biological assemblies.

Is the H4 histone tail intrinsically disordered or intrinsically multifunctional?

The structural versatility of histone tails is one of the key elements in the organisation of chromatin, which allows for the compact storage of genomic information. However, this structural diversity also complicates experimental and computational studies. Here, the potential and free energy landscape for the isolated and bound H4 histone tail are explored. The landscapes exhibit a set of distinct structural ensembles separated by high energy barriers, with little difference between isolated and bound tails. This consistency is a desirable feature that facilitates the formation of transient interactions, which are required for the liquid-like chromatin organisation. The existence of multiple, distinct structures on a multifunnel energy landscape is likely to be associated with multifunctionality, i.e. a set of evolved, distinct functions. Contrasting it with previously reported results for other disordered peptides, this type of landscape may be associated with a conformational selection based binding mechanism. Given the similarity to other systems exhibiting similar multifunnel energy landscapes, the disorder in histone tails might be better described in context of multifunctionality.

Structural Heterogeneities of the Ribosome: New Frontiers and Opportunities for Cryo-EM

The extent of ribosomal heterogeneity has caught increasing interest over the past few years, as recent studies have highlighted the presence of structural variations of the ribosome. More precisely, the heterogeneity of the ribosome covers multiple scales, including the dynamical aspects of ribosomal motion at the single particle level, specialization at the cellular and subcellular scale, or evolutionary differences across species. Upon solving the ribosome atomic structure at medium to high resolution, cryogenic electron microscopy (cryo-EM) has enabled investigating all these forms of heterogeneity. In this review, we present some recent advances in quantifying ribosome heterogeneity, with a focus on the conformational and evolutionary variations of the ribosome and their functional implications. These efforts highlight the need for new computational methods and comparative tools, to comprehensively model the continuous conformational transition pathways of the ribosome, as well as its evolution. While developing these methods presents some important challenges, it also provides an opportunity to extend our interpretation and usage of cryo-EM data, which would more generally benefit the study of molecular dynamics and evolution of proteins and other complexes.