CafeMol: A Coarse-Grained Biomolecular Simulator for Simulating Proteins at Work

Allosteric Regulation of Enzymatic Catalysis through Mechanical Force

Mechanical force has been increasingly recognized to play crucial roles in regulating various cellular processes, which has inspired wide interest in elucidating the biophysical mechanism underlying these mechanobiological processes. In this work, we investigate the mechanical regulation of enzyme catalysis by developing a residue-resolved computational model capable of describing the full catalytic cycle of enzymes under mechanical force. Intriguingly, for a model enzyme, adenylate kinase, we showed that applying tensile forces with biologically relevant strength can increase the enzymatic activity. Further analysis showed that mechanical tensile force allosterically modifies the global free energy landscape and conformational dynamics of the protein, which then promotes the rate-limiting product release step of the enzymatic cycle. The effect of mechanical allostery on enzyme catalysis depends on the intrinsic conformational propensity of the enzymes. The crucial role of mechanical allostery in enzymatic catalysis elucidated in this work sheds important insights into the biophysical principle of enzymatic regulation and suggests a possible strategy for fine-tuning the functioning dynamics of biological enzymes.

Solution AFM Imaging and Coarse-grained Molecular Modeling of Yeast Condensin Structural Variation Coupled to the ATP Hydrolysis Cycle

Condensin is a protein complex that regulates chromatin structural changes during mitosis. It varies the molecular conformation through the ATP hydrolysis cycle and extrudes DNA loops into its ring-like structure as a molecular motor. Condensin contains Smc2 and Smc4, in which a coiled-coil arm tethers the hinge and head domains and dimerizes at the hinge. ATPs bind between the heads, induce their engagement, and are hydrolyzed to promote their disengagement. Previously, we performed solution atomic force microscopy (AFM) imaging of yeast condensin holo-complex with ATP and conducted flexible molecular fitting, obtaining the hinge structure with open conformation. However, it has yet to be clarified how the opening/closing of the hinge is coupled to the ATP hydrolysis cycle. In this study, we performed solution AFM imaging in the presence and absence of varying nucleotides, including AMP-PNP, ATPγS, and ADP. Furthermore, we conducted molecular dynamics simulations of an Smc2/4 heterodimer and selected the structure best representing each AFM image. Our results suggested that head engagement upon ATP binding is coupled to hinge opening and that the N-terminal region of Brn1, one of the accessory subunits, re-associates to the Smc2 head after ADP release.

Molecular dynamics simulations of human cohesin subunits identify DNA binding sites and their potential roles in DNA loop extrusion

The SMC complex cohesin mediates interphase chromatin structural formation in eukaryotic cells through DNA loop extrusion. Here, we sought to investigate its mechanism using molecular dynamics simulations. To achieve this, we first constructed the amino-acid-residue-resolution structural models of the cohesin subunits, SMC1, SMC3, STAG1, and NIPBL. By simulating these subunits with double-stranded DNA molecules, we predicted DNA binding patches on each subunit and quantified the affinities of these patches to DNA using their dissociation rate constants as a proxy. Then, we constructed the structural model of the whole cohesin complex and mapped the predicted high-affinity DNA binding patches on the structure. From the spatial relations of the predicted patches, we identified that multiple patches on the SMC1, SMC3, STAG1, and NIPBL subunits form a DNA clamping patch group. The simulations of the whole complex with double-stranded DNA molecules suggest that this patch group facilitates DNA bending and helps capture a DNA segment in the cohesin ring formed by the SMC1 and SMC3 subunits. In previous studies, these have been identified as critical steps in DNA loop extrusion. Therefore, this study provides experimentally testable predictions of DNA binding sites implicated in previously proposed DNA loop extrusion mechanisms and highlights the essential roles of the accessory subunits STAG1 and NIPBL in the mechanism.

Challenges in simulating whole virus particles and how to fix them with the SIRAH force field

Current developments in specialized software and computer power make the simulation of large molecular assemblies a technical possibility despite their computational cost. Coarse-grained (CG) approaches simplify molecular complexity and reduce computational costs while preserving intermolecular physical/chemical interactions. These methods enable virus simulations, making them more accessible to research groups with limited supercomputing resources. However, setting up and running molecular dynamics simulations of multimillion systems requires specialized molecular modeling, editing, and visualization skills. Moreover, many issues related to the computational setup, the choice of simulation engines, and the force fields that rule the intermolecular interactions require particular attention and are key to attaining a realistic description of viral systems at the fully atomistic or CG levels. Here, we provide an overview of the current challenges in simulating entire virus particles and the potential of the SIRAH force field to address these challenges through its implementations for CG and multiscale simulations.

The ATPase asymmetry: Novel computational insight into coupling diverse FO motors with tripartite F1

ATP synthase, a crucial enzyme for cellular bioenergetics, operates via the coordinated coupling of an FO motor, which presents variable symmetry, and a tripartite F1 motor. Despite extensive research, the understanding of their coupling dynamics, especially with non-10-fold symmetrical FO motors, remains incomplete. This study investigates the coupling patterns between eightfold and ninefold FO motors and the constant threefold F1 motor using coarse-grained molecular dynamics simulations. We unveil that in the case of a ninefold FO motor, a 3-3-3 motion is most likely to occur, whereas a 3-3-2 motion predominates with an eightfold FO motor. Furthermore, our findings propose a revised model for the coupling method, elucidating that the pathways' energy usage is primarily influenced by F1 rotation and conformational changes hindered by the b-subunits. Our results present a crucial step toward comprehending the energy landscape and mechanisms governing ATP synthase operation.

Chromosomal domain formation by archaeal SMC, a roadblock protein, and DNA structure

In eukaryotes, structural maintenance of chromosomes (SMC) complexes form topologically associating domains (TADs) by extruding DNA loops and being stalled by roadblock proteins. It remains unclear whether a similar mechanism of domain formation exists in prokaryotes. Using high-resolution chromosome conformation capture sequencing, we show that an archaeal homolog of the bacterial Smc-ScpAB complex organizes the genome of Thermococcus kodakarensis into TAD-like domains. We find that TrmBL2, a nucleoid-associated protein that forms a stiff nucleoprotein filament, stalls the T. kodakarensis SMC complex and establishes a boundary at the site-specific recombination site dif. TrmBL2 stalls the SMC complex at tens of additional non-boundary loci with lower efficiency. Intriguingly, the stalling efficiency is correlated with structural properties of underlying DNA sequences. Our study illuminates a eukaryotic-like mechanism of domain formation in archaea and a role of intrinsic DNA structure in large-scale genome organization.

FlowBack: A Generalized Flow-Matching Approach for Biomolecular Backmapping

Coarse-grained models have become ubiquitous in biomolecular modeling tasks aimed at studying slow dynamical processes such as protein folding and DNA hybridization. These models can considerably accelerate sampling but it remains challenging to accurately and efficiently restore all-atom detail to the coarse-grained trajectory, which can be vital for detailed understanding of molecular mechanisms and calculation of observables contingent on all-atom coordinates. In this work, we introduce FlowBack as a deep generative model employing a flow-matching objective to map samples from a coarse-grained prior distribution to an all-atom data distribution. We construct our prior distribution to be agnostic to the coarse-grained map and molecular type. A protein-specific model trained on ∼65k structures from the Protein Data Bank achieves state-of-the-art performance on structural metrics compared to previous generative and rules-based approaches in applications to static PDB structures, all-atom simulations of fast-folding proteins, and coarse-grained trajectories generated by a machine-learned force field. A DNA-protein model trained on ∼1.5k DNA-protein complexes achieves excellent reconstruction and generative capabilities on static DNA-protein complexes from the Protein Data Bank as well as on out-of-distribution coarse-grained dynamical simulations of DNA-protein complexation. FlowBack offers an accurate, efficient, and easy-to-use tool to recover all-atom structures from coarse-grained molecular simulations with higher robustness and fewer steric clashes than previous approaches. We make FlowBack freely available to the community as an open source Python package.

The conformational space of RNase P RNA in solution

RNA conformational diversity has fundamental biological roles1-5, but direct visualization of its full conformational space in solution has not been possible using traditional biophysical techniques. Using solution atomic force microscopy, a deep neural network and statistical analyses, we show that the ribonuclease P (RNase P) RNA adopts heterogeneous conformations consisting of a conformationally invariant core and highly flexible peripheral structural elements that sample a broad conformational space, with amplitudes as large as 20-60 Å in a multitude of directions, with very low net energy cost. Increasing Mg2+ drives compaction and enhances enzymatic activity, probably by narrowing the conformational space. Moreover, analyses of the correlations and anticorrelations between spatial flexibility and sequence conservation suggest that the functional roles of both the structure and dynamics of key regions are embedded in the primary sequence. These findings reveal the structure-dynamics basis for the embodiment of both enzymatic precision and substrate promiscuity in the RNA component of the RNase P. Mapping the conformational space of the RNase P RNA demonstrates a new general approach to studying RNA structure and dynamics.

Determining structures of RNA conformers using AFM and deep neural networks

Much of the human genome is transcribed into RNAs1, many of which contain structural elements that are important for their function. Such RNA molecules-including those that are structured and well-folded2-are conformationally heterogeneous and flexible, which is a prerequisite for function3,4, but this limits the applicability of methods such as NMR, crystallography and cryo-electron microscopy for structure elucidation. Moreover, owing to the lack of a large RNA structure database, and no clear correlation between sequence and structure, approaches such as AlphaFold5 for protein structure prediction do not apply to RNA. Therefore, determining the structures of heterogeneous RNAs remains an unmet challenge. Here we report holistic RNA structure determination method using atomic force microscopy, unsupervised machine learning and deep neural networks (HORNET), a novel method for determining three-dimensional topological structures of RNA using atomic force microscopy images of individual molecules in solution. Owing to the high signal-to-noise ratio of atomic force microscopy, this method is ideal for capturing structures of large RNA molecules in distinct conformations. In addition to six benchmark cases, we demonstrate the utility of HORNET by determining multiple heterogeneous structures of RNase P RNA and the HIV-1 Rev response element (RRE) RNA. Thus, our method addresses one of the major challenges in determining heterogeneous structures of large and flexible RNA molecules, and contributes to the fundamental understanding of RNA structural biology.

Effect of Pregenomic RNA on the Mechanical Stability of HBV Capsid by Coarse-Grained Molecular Simulations

Hepatitis B virus (HBV) is a double-stranded DNA virus, but its life cycle involves an intermediate stage, during which pregenomic RNA (pgRNA) is encapsulated in the capsid and then reverse-transcribed into the minus DNA strand. These immature HBV virions are the key target for antiviral drug discovery. In this study, we investigate the flexibility and mechanical stability of the HBV capsid containing pgRNA by employing residue-resolved coarse-grained molecular dynamics simulations. The results showed that the presence of pgRNA tends to decrease the overall flexibility of the capsid. In addition, the symmetrically arranged subunits of the capsid show asymmetry in the dominant modes of the conformational fluctuations with or without the presence of pgRNA. Furthermore, the simulations revealed that the presence of pgRNA enhances the overall mechanical stability of the virion particle. Electrostatic interactions between the disordered CTD of capsid and pgRNA were found to play a crucial role in modulating viral mechanical stability. Decreasing the electrostatic interactions by CTD phosphorylation or high salt concentration significantly reduces the mechanical stability of the HBV capsid containing pgRNA. Finally, the 2-fold symmetric sites have been proposed to be the most vulnerable to rupture during the initial stages of capsid disassembly. These findings could enhance our understanding of the physical basis of viral invasion and provide valuable insights into the development of antiviral drugs.

Molecular mechanism of parental H3/H4 recycling at a replication fork

In chromatin replication, faithful recycling of histones from parental DNA to replicated strands is essential for maintaining epigenetic information across generations. A previous experiment has revealed that disrupting interactions between the N-terminal tail of Mcm2, a subunit in DNA replication machinery, and a histone H3/H4 tetramer perturb the recycling. However, the molecular pathways and the factors that regulate the ratio recycled to each strand and the destination location are yet to be revealed. Here, we performed molecular dynamics simulations of yeast DNA replication machinery, an H3/H4 tetramer, and replicated DNA strands. The simulations demonstrated that histones are recycled via Cdc45-mediated and unmediated pathways without histone chaperones, as our in vitro biochemical assays supported. Also, RPA binding regulated the ratio recycled to each strand, whereas DNA bending by Pol ε modulated the destination location. Together, the simulations provided testable hypotheses, which are vital for elucidating the molecular mechanisms of histone recycling.

Sodium Ion-Induced Structural Transition on the Surface of a DNA-Interacting Protein

Protein surfaces have pivotal roles in interactions between proteins and other biological molecules. However, the structural dynamics of protein surfaces have rarely been explored and are poorly understood. Here, the surface of a single-stranded DNA (ssDNA) binding protein (SSB) with four DNA binding domains that bind ssDNA in binding site sizes of 35, 56, and 65 nucleotides per tetramer is investigated. Using oligonucleotides as probes to sense the charged surface, NaCl induces a two-state structural transition on the SSB surface even at moderate concentrations. Chelation of sodium ions with charged amino acids alters the network of hydrogen bonds and/or salt bridges on the surface. Such changes are associated with changes in the electrostatic potential landscape and interaction mode. These findings advance the understanding of the molecular mechanism underlying the enigmatic salt-induced transitions between different DNA binding site sizes of SSBs. This work demonstrates that monovalent salt is a key regulator of biomolecular interactions that not only play roles in non-specific electrostatic screening effects as usually assumed but also may configure the surface of proteins to contribute to the effective regulation of biomolecular recognition and other downstream events.

Asymmetric fluctuation of overlapping dinucleosome studied by cryoelectron microscopy and small-angle X-ray scattering

Nucleosome remodelers modify the local structure of chromatin to release the region from nucleosome-mediated transcriptional suppression. Overlapping dinucleosomes (OLDNs) are nucleoprotein complexes formed around transcription start sites as a result of remodeling, and they consist of two nucleosome moieties: a histone octamer wrapped by DNA (octasome) and a histone hexamer wrapped by DNA (hexasome). While OLDN formation alters chromatin accessibility to proteins, the structural mechanism behind this process is poorly understood. Thus, this study investigated the characteristics of structural fluctuations in OLDNs. First, multiple structures of the OLDN were visualized through cryoelectron microscopy (cryoEM), providing an overview of the tilting motion of the hexasome relative to the octasome at the near-atomistic resolution. Second, small-angle X-ray scattering (SAXS) revealed the presence of OLDN conformations with a larger radius of gyration than cryoEM structures. A more complete description of OLDN fluctuation was proposed by SAXS-based ensemble modeling, which included possible transient structures. The ensemble model supported the tilting motion of the OLDN outlined by the cryoEM models, further suggesting the presence of more diverse conformations. The amplitude of the relative tilting motion of the hexasome was larger, and the nanoscale fluctuation in distance between the octasome and hexasome was also proposed. The cryoEM models were found to be mapped in the energetically stable region of the conformational distribution of the ensemble. Exhaustive complex modeling using all conformations that appeared in the structural ensemble suggested that conformational and motional asymmetries of the OLDN result in asymmetries in the accessibility of OLDN-binding proteins.

Bidirectional Allostery Mechanism in Catch-Bond Formation of CD44 Mediated Cell Adhesion

Catch-bonds, whereby noncovalent ligand-receptor interactions are counterintuitively reinforced by tensile forces, play a major role in cell adhesion under mechanical stress. A basic prerequisite for catch-bond formation, as implicated in classic catch-bond models, is that force-induced remodeling of the ligand binding interface occurs prior to bond rupture. However, what strategy receptor proteins utilize to meet such specific kinetic control remains elusive. Here we report a bidirectional allostery mechanism of catch-bond formation based on theoretical and molecular dynamics simulation studies. Binding of ligand allosterically reduces the threshold force for unlocking of otherwise stably folded force-sensing element (i.e., forward allostery), so that a much smaller tensile force can trigger the conformational switching of receptor protein to high binding-strength state via backward allosteric coupling before bond rupture. Such bidirectional allostery fulfills the specific kinetic control required by catch-bond formation and is likely to be commonly utilized in cell adhesion. The essential thermodynamic and kinetic features of receptor proteins essential for catch-bond formation were identified.

Theoretical insights into rotary mechanism of MotAB in the bacterial flagellar motor

Many bacteria enable locomotion by rotating their flagellum. It has been suggested that this rotation is realized by the rotary motion of the stator unit, MotAB, which is driven by proton transfer across the membrane. Recent cryo-electron microscopy studies have revealed a 5:2 MotAB configuration, in which a MotB dimer is encircled by a ring-shaped MotA pentamer. Although the structure implicates the rotary motion of the MotA wheel around the MotB axle, the molecular mechanisms of rotary motion and how they are coupled with proton transfer across the membrane remain elusive. In this study, we built a structure-based computational model for Campylobacter jejuni MotAB, conducted comprehensive protonation-state-dependent molecular dynamics simulations, and revealed a plausible proton-transfer-coupled rotation pathway. The model assumes rotation-dependent proton transfer, in which proton uptake from the periplasmic side to the conserved aspartic acid in MotB is followed by proton hopping to the MotA proton-carrying site, followed by proton export to the CP. We suggest that, by maintaining two of the proton-carrying sites of MotA in the deprotonated state, the MotA pentamer robustly rotates by ∼36° per proton transfer across the membrane. Our results provide a structure-based mechanistic model of the rotary motion of MotAB in bacterial flagellar motors and provide insights into various ion-driven rotary molecular motors.

Nucleosomal DNA unwinding pathway through canonical and non-canonical histone disassembly

The nucleosome including H2A.B, a mammalian-specific H2A variant, plays pivotal roles in spermatogenesis, embryogenesis, and oncogenesis, indicating unique involvement in transcriptional regulation distinct from canonical H2A nucleosomes. Despite its significance, the exact regulatory mechanism remains elusive. This study utilized solid-state nanopores to investigate DNA unwinding dynamics, applying local force between DNA and histones. Comparative analysis of canonical H2A and H2A.B nucleosomes demonstrated that the H2A.B variant required a lower voltage for complete DNA unwinding. Furthermore, synchronization analysis and molecular dynamics simulations indicate that the H2A.B variant rapidly unwinds DNA, causing the H2A-H2B dimer to dissociate from DNA immediately upon disassembly of the histone octamer. In contrast, canonical H2A nucleosomes unwind DNA at a slower rate, suggesting that the H2A-H2B dimer undergoes a state of stacking at the pore. These findings suggest that nucleosomal DNA in the H2A.B nucleosomes undergoes a DNA unwinding process involving histone octamer disassembly distinct from that of canonical H2A nucleosomes, enabling smoother unwinding. The integrated approach of MD simulations and nanopore measurements is expected to evolve into a versatile tool for studying molecular interactions, not only within nucleosomes but also through the forced dissociation of DNA-protein complexes.

The molecular basis for hydrodynamic properties of PEGylated human serum albumin

Polyethylene glycol (PEG) conjugation provides a protective modification that enhances the pharmacokinetics and solubility of proteins for therapeutic use. A knowledge of the structural ensemble of these PEGylated proteins is necessary to understand the molecular details that contribute to their hydrodynamic and colligative properties. Because of the large size and dynamic flexibility of pharmaceutically important PEGylated proteins, the determination of structure is challenging. In addition, the hydration of these conjugates that contain large polymers is difficult to determine with traditional methods that identify only first shell hydration water, which does not account for the complete hydrodynamic volume of a macromolecule. Here, we demonstrate that structural ensembles, generated by coarse-grained simulations, can be analyzed with HullRad and used to predict sedimentation coefficients and concentration-dependent hydrodynamic and diffusion nonideality coefficients of PEGylated proteins. A knowledge of these concentration-dependent properties enhances the ability to design and analyze new modified protein therapeutics. HullRad accomplishes this analysis by effectively accounting for the complete hydration of a macromolecule, including that of flexible polymers.

GENESIS 2.1: High-Performance Molecular Dynamics Software for Enhanced Sampling and Free-Energy Calculations for Atomistic, Coarse-Grained, and Quantum Mechanics/Molecular Mechanics Models

GENeralized-Ensemble SImulation System (GENESIS) is a molecular dynamics (MD) software developed to simulate the conformational dynamics of a single biomolecule, as well as molecular interactions in large biomolecular assemblies and between multiple biomolecules in cellular environments. To achieve the latter purpose, the earlier versions of GENESIS emphasized high performance in atomistic MD simulations on massively parallel supercomputers, with or without graphics processing units (GPUs). Here, we implemented multiscale MD simulations that include atomistic, coarse-grained, and hybrid quantum mechanics/molecular mechanics (QM/MM) calculations. They demonstrate high performance and are integrated with enhanced conformational sampling algorithms and free-energy calculations without using external programs except for the QM programs. In this article, we review new functions, molecular models, and other essential features in GENESIS version 2.1 and discuss ongoing developments for future releases.

GENESIS CGDYN: large-scale coarse-grained MD simulation with dynamic load balancing for heterogeneous biomolecular systems

Residue-level coarse-grained (CG) molecular dynamics (MD) simulation is widely used to investigate slow biological processes that involve multiple proteins, nucleic acids, and their complexes. Biomolecules in a large simulation system are distributed non-uniformly, limiting computational efficiency with conventional methods. Here, we develop a hierarchical domain decomposition scheme with dynamic load balancing for heterogeneous biomolecular systems to keep computational efficiency even after drastic changes in particle distribution. These schemes are applied to the dynamics of intrinsically disordered protein (IDP) droplets. During the fusion of two droplets, we find that the changes in droplet shape correlate with the mixing of IDP chains. Additionally, we simulate large systems with multiple IDP droplets, achieving simulation sizes comparable to those observed in microscopy. In our MD simulations, we directly observe Ostwald ripening, a phenomenon where small droplets dissolve and their molecules redeposit into larger droplets. These methods have been implemented in CGDYN of the GENESIS software, offering a tool for investigating mesoscopic biological processes using the residue-level CG models.

Effect of α-tubulin acetylation on the doublet microtubule structure

Acetylation of α-tubulin at the lysine 40 residue (αK40) by αTAT1/MEC-17 acetyltransferase modulates microtubule properties and occurs in most eukaryotic cells. Previous literatures suggest that acetylated microtubules are more stable and damage resistant. αK40 acetylation is the only known microtubule luminal post-translational modification site. The luminal location suggests that the modification tunes the lateral interaction of protofilaments inside the microtubule. In this study, we examined the effect of tubulin acetylation on the doublet microtubule (DMT) in the cilia of Tetrahymena thermophila using a combination of cryo-electron microscopy, molecular dynamics, and mass spectrometry. We found that αK40 acetylation exerts a small-scale effect on the DMT structure and stability by influencing the lateral rotational angle. In addition, comparative mass spectrometry revealed a link between αK40 acetylation and phosphorylation in cilia.

Computational Approaches to Predict Protein-Protein Interactions in Crowded Cellular Environments

Investigating protein-protein interactions is crucial for understanding cellular biological processes because proteins often function within molecular complexes rather than in isolation. While experimental and computational methods have provided valuable insights into these interactions, they often overlook a critical factor: the crowded cellular environment. This environment significantly impacts protein behavior, including structural stability, diffusion, and ultimately the nature of binding. In this review, we discuss theoretical and computational approaches that allow the modeling of biological systems to guide and complement experiments and can thus significantly advance the investigation, and possibly the predictions, of protein-protein interactions in the crowded environment of cell cytoplasm. We explore topics such as statistical mechanics for lattice simulations, hydrodynamic interactions, diffusion processes in high-viscosity environments, and several methods based on molecular dynamics simulations. By synergistically leveraging methods from biophysics and computational biology, we review the state of the art of computational methods to study the impact of molecular crowding on protein-protein interactions and discuss its potential revolutionizing effects on the characterization of the human interactome.

Observing Dynamic Conformational Changes within the Coiled-Coil Domain of Different Laminin Isoforms Using High-Speed Atomic Force Microscopy

Laminins are trimeric glycoproteins with important roles in cell-matrix adhesion and tissue organization. The laminin α, ß, and γ-chains have short N-terminal arms, while their C-termini are connected via a triple coiled-coil domain, giving the laminin molecule a well-characterized cross-shaped morphology as a result. The C-terminus of laminin alpha chains contains additional globular laminin G-like (LG) domains with important roles in mediating cell adhesion. Dynamic conformational changes of different laminin domains have been implicated in regulating laminin function, but so far have not been analyzed at the single-molecule level. High-speed atomic force microscopy (HS-AFM) is a unique tool for visualizing such dynamic conformational changes under physiological conditions at sub-second temporal resolution. After optimizing surface immobilization and imaging conditions, we characterized the ultrastructure of laminin-111 and laminin-332 using HS-AFM timelapse imaging. While laminin-111 features a stable S-shaped coiled-coil domain displaying little conformational rearrangement, laminin-332 coiled-coil domains undergo rapid switching between straight and bent conformations around a defined central molecular hinge. Complementing the experimental AFM data with AlphaFold-based coiled-coil structure prediction enabled us to pinpoint the position of the hinge region, as well as to identify potential molecular rearrangement processes permitting hinge flexibility. Coarse-grained molecular dynamics simulations provide further support for a spatially defined kinking mechanism in the laminin-332 coiled-coil domain. Finally, we observed the dynamic rearrangement of the C-terminal LG domains of laminin-111 and laminin-332, switching them between compact and open conformations. Thus, HS-AFM can directly visualize molecular rearrangement processes within different laminin isoforms and provide dynamic structural insight not available from other microscopy techniques.

Enhanced Coarse-Grained Molecular Dynamics Simulation with a Smoothed Hybrid Potential Using a Neural Network Model

In all-atom (AA) molecular dynamics (MD) simulations, the rugged energy profile of the force field makes it challenging to reproduce spontaneous structural changes in biomolecules within a reasonable calculation time. Existing coarse-grained (CG) models, in which the energy profile is set to a global minimum around the initial structure, are unsuitable to explore the structural dynamics between metastable states far away from the initial structure without any bias. In this study, we developed a new hybrid potential composed of an artificial intelligence (AI) potential and minimal CG potential related to the statistical bond length and excluded volume interactions to accelerate the transition dynamics while maintaining the protein character. The AI potential is trained by energy matching using a diverse structural ensemble sampled via multicanonical (Mc) MD simulation and the corresponding AA force field energy, profile of which is smoothed by energy minimization. By applying the new methodology to chignolin and TrpCage, we showed that the AI potential can predict the AA energy with significantly high accuracy, as indicated by a correlation coefficient (R-value) between the true and predicted energies exceeding 0.89. In addition, we successfully demonstrated that CGMD simulation based on the smoothed hybrid potential can significantly enhance the transition dynamics between various metastable states while preserving protein properties compared to those obtained with conventional CGMD and AAMD.

Genome-wide mapping and cryo-EM structural analyses of the overlapping tri-nucleosome composed of hexasome-hexasome-octasome moieties

The nucleosome is a fundamental unit of chromatin in which about 150 base pairs of DNA are wrapped around a histone octamer. The overlapping di-nucleosome has been proposed as a product of chromatin remodeling around the transcription start site, and previously found as a chromatin unit, in which about 250 base pairs of DNA continuously bind to the histone core composed of a hexamer and an octamer. In the present study, our genome-wide analysis of human cells suggests another higher nucleosome stacking structure, the overlapping tri-nucleosome, which wraps about 300-350 base-pairs of DNA in the region downstream of certain transcription start sites of actively transcribed genes. We determine the cryo-electron microscopy (cryo-EM) structure of the overlapping tri-nucleosome, in which three subnucleosome moieties, hexasome, hexasome, and octasome, are associated by short connecting DNA segments. Small angle X-ray scattering and coarse-grained molecular dynamics simulation analyses reveal that the cryo-EM structure of the overlapping tri-nucleosome may reflect its structure in solution. Our findings suggest that nucleosome stacking structures composed of hexasome and octasome moieties may be formed by nucleosome remodeling factors around transcription start sites for gene regulation.

A new twist on PIFE: photoisomerisation-related fluorescence enhancement

PIFE was first used as an acronym for protein-induced fluorescence enhancement, which refers to the increase in fluorescence observed upon the interaction of a fluorophore, such as a cyanine, with a protein. This fluorescence enhancement is due to changes in the rate ofcis/transphotoisomerisation. It is clear now that this mechanism is generally applicable to interactions with any biomolecule. In this review, we propose that PIFE is thereby renamed according to its fundamental working principle as photoisomerisation-related fluorescence enhancement, keeping the PIFE acronym intact. We discuss the photochemistry of cyanine fluorophores, the mechanism of PIFE, its advantages and limitations, and recent approaches to turning PIFE into a quantitative assay. We provide an overview of its current applications to different biomolecules and discuss potential future uses, including the study of protein-protein interactions, protein-ligand interactions and conformational changes in biomolecules.

Extracellular calcium functions as a molecular glue for transmembrane helices to activate the scramblase Xkr4

The "eat me" signal, phosphatidylserine is exposed on the surface of dying cells by phospholipid scrambling. Previously, we showed that the Xkr family protein Xkr4 is activated by caspase-mediated cleavage and binding of the XRCC4 fragment. Here, we show that extracellular calcium is an additional factor needed to activate Xkr4. The constitutively active mutant of Xkr4 is found to induce phospholipid scrambling in an extracellular, but not intracellular, calcium-dependent manner. Importantly, other Xkr family members also require extracellular calcium for activation. Alanine scanning shows that D123 and D127 of TM1 and E310 of TM3 coordinate calcium binding. Moreover, lysine scanning demonstrates that the E310K mutation-mediated salt bridge between TM1 and TM3 bypasses the requirement of calcium. Cysteine scanning proves that disulfide bond formation between TM1 and TM3 also activates phospholipid scrambling without calcium. Collectively, this study shows that extracellular calcium functions as a molecular glue for TM1 and TM3 of Xkr proteins for activation, thus demonstrating a regulatory mechanism for multi-transmembrane region-containing proteins.

The SIRAH force field: A suite for simulations of complex biological systems at the coarse-grained and multiscale levels

The different combinations of molecular dynamics simulations with coarse-grained representations have acquired considerable popularity among the scientific community. Especially in biocomputing, the significant speedup granted by simplified molecular models opened the possibility of increasing the diversity and complexity of macromolecular systems, providing realistic insights on large assemblies for more extended time windows. However, a holistic view of biological ensembles' structural and dynamic features requires a self-consistent force field, namely, a set of equations and parameters that describe the intra and intermolecular interactions among moieties of diverse chemical nature (i.e., nucleic and amino acids, lipids, solvent, ions, etc.). Nevertheless, examples of such force fields are scarce in the literature at the fully atomistic and coarse-grained levels. Moreover, the number of force fields capable of handling simultaneously different scales is restricted to a handful. Among those, the SIRAH force field, developed in our group, furnishes a set of topologies and tools that facilitate the setting up and running of molecular dynamics simulations at the coarse-grained and multiscale levels. SIRAH uses the same classical pairwise Hamiltonian function implemented in the most popular molecular dynamics software. In particular, it runs natively in AMBER and Gromacs engines, and porting it to other simulation packages is straightforward. This review describes the underlying philosophy behind the development of SIRAH over the years and across families of biological molecules, discussing current limitations and future implementations.

Molecular Design of FRET Probes Based on Domain Rearrangement of Protein Disulfide Isomerase for Monitoring Intracellular Redox Status

Multidomain proteins can exhibit sophisticated functions based on cooperative interactions and allosteric regulation through spatial rearrangements of the multiple domains. This study explored the potential of using multidomain proteins as a basis for Förster resonance energy transfer (FRET) biosensors, focusing on protein disulfide isomerase (PDI) as a representative example. PDI, a well-studied multidomain protein, undergoes redox-dependent conformational changes, enabling the exposure of a hydrophobic surface extending across the b' and a' domains that serves as the primary binding site for substrates. Taking advantage of the dynamic domain rearrangements of PDI, we developed FRET-based biosensors by fusing the b' and a' domains of thermophilic fungal PDI with fluorescent proteins as the FRET acceptor and donor, respectively. Both experimental and computational approaches were used to characterize FRET efficiency in different redox states. In vitro and in vivo evaluations demonstrated higher FRET efficiency of this biosensor in the oxidized form, reflecting the domain rearrangement and its responsiveness to intracellular redox environments. This novel approach of exploiting redox-dependent domain dynamics in multidomain proteins offers promising opportunities for designing innovative FRET-based biosensors with potential applications in studying cellular redox regulation and beyond.

FO-F1 coupling and symmetry mismatch in ATP synthase resolved in every FO rotation step

FOF1 ATP synthase, a ubiquitous enzyme that synthesizes most ATP in living cells, is composed of two rotary motors: a membrane-embedded proton-driven FO motor and a catalytic F1 motor. These motors share both central and peripheral stalks. Although both FO and F1 have pseudo-symmetric structures, their symmetries do not match. How symmetry mismatch is solved remains elusive because of the missing intermediate structures of the rotational steps. Here, for the case of Bacillus PS3 ATP synthases with three- and 10-fold symmetries in F1 and FO, respectively, we uncovered the mechanical couplings between FO and F1 at every 36° rotation step via molecular dynamics simulations and comparative studies of cryoelectron microscopy (cryo-EM) structures from three species. We found that the mismatch could be solved using several elements: 1) the F1 head partially rotates relative to the FO a subunit via elastic distortion of the b subunits, 2) the rotor is twisted, and 3) comparisons of cryo-EM structures further suggest that the c ring rotary angles can deviate from the symmetric ones. In addition, the F1 motor may have non-canonical structures, relieving stronger frustration. Thus, we provide new insights for solving the symmetry mismatch problem.

Modeling the effects of phosphorylation on phase separation of the FUS low-complexity domain

Aggregation of the RNA-binding protein fused in sarcoma (FUS) is a hallmark of neurodegenerative diseases. Phosphorylation of Ser/Thr in the FUS low-complexity domain (FUS-LC) may regulate phase separation of FUS and prevent pathological aggregation in cells. However, many details of this process remain elusive to date. In this work, we systematically investigated the phosphorylation of FUS-LC and the underlying molecular mechanism by molecular dynamics (MD) simulations and free energy calculations. The results clearly show that phosphorylation can destroy the fibril core structure of FUS-LC by breaking interchain interactions, particularly contacts involving residues like Tyr, Ser, and Gln. Among the six phosphorylation sites, Ser61 and Ser84 may have more important effects on the stability of the fibril core. Our study reveals structural and dynamic details of FUS-LC phase separation modulated by phosphorylation.

Crystal structure of Escherichia coli thiamine pyrophosphate-sensing riboswitch in the apo state

The thiamine pyrophosphate (TPP)-sensing riboswitch is one of the earliest discovered and most widespread riboswitches. Numerous structural studies have been reported for this riboswitch bound with various ligands. However, the ligand-free (apo) structure remains unknown. Here, we report a 3.1 Å resolution crystal structure of Escherichia coli TPP riboswitch in the apo state, which exhibits an extended, Y-shaped conformation further supported by small-angle X-ray scattering data and driven molecular dynamics simulations. The loss of ligand interactions results in helical uncoiling of P5 and disruption of the key tertiary interaction between the sensory domains. Opening of the aptamer propagates to the gene-regulatory P1 helix and generates the key conformational flexibility needed for the switching behavior. Much of the ligand-binding site at the three-way junction is unaltered, thereby maintaining a partially preformed pocket. Together, these results paint a dynamic picture of the ligand-induced conformational changes in TPP riboswitches that confer conditional gene regulation.

Generation of unfolded outer membrane protein ensembles defined by hydrodynamic properties

Outer membrane proteins (OMPs) must exist as an unfolded ensemble while interacting with a chaperone network in the periplasm of Gram-negative bacteria. Here, we developed a method to model unfolded OMP (uOMP) conformational ensembles using the experimental properties of two well-studied OMPs. The overall sizes and shapes of the unfolded ensembles in the absence of a denaturant were experimentally defined by measuring the sedimentation coefficient as a function of urea concentration. We used these data to model a full range of unfolded conformations by parameterizing a targeted coarse-grained simulation protocol. The ensemble members were further refined by short molecular dynamics simulations to reflect proper torsion angles. The final conformational ensembles have polymer properties different from unfolded soluble and intrinsically disordered proteins and reveal inherent differences in the unfolded states that necessitate further investigation. Building these uOMP ensembles advances the understanding of OMP biogenesis and provides essential information for interpreting structures of uOMP-chaperone complexes.

Micelle-like clusters in phase-separated Nanog condensates: A molecular simulation study

The phase separation model for transcription suggests that transcription factors (TFs), coactivators, and RNA polymerases form biomolecular condensates around active gene loci and regulate transcription. However, the structural details of condensates remain elusive. In this study, for Nanog, a master TF in mammalian embryonic stem cells known to form protein condensates in vitro, we examined protein structures in the condensates using residue-level coarse-grained molecular simulations. Human Nanog formed micelle-like clusters in the condensate. In the micelle-like cluster, the C-terminal disordered domains, including the tryptophan repeat (WR) regions, interacted with each other near the cluster center primarily via hydrophobic interaction. In contrast, hydrophilic disordered N-terminal and DNA-binding domains were exposed on the surface of the clusters. Electrostatic attractions of these surface residues were responsible for bridging multiple micelle-like structures in the condensate. The micelle-like structure and condensate were dynamic and liquid-like. Mutation of tryptophan residues in the WR region which was implicated to be important for a Nanog function resulted in dissolution of the Nanog condensate. Finally, to examine the impact of Nanog cluster to DNA, we added DNA fragments to the Nanog condensate. Nanog DNA-binding domains exposed to the surface of the micelle-like cluster could recruit more than one DNA fragments, making DNA-DNA distance shorter.

Histone chaperone Nap1 dismantles an H2A/H2B dimer from a partially unwrapped nucleosome

DNA translocases, such as RNA polymerases, inevitably collide with nucleosomes on eukaryotic chromatin. Upon these collisions, histone chaperones are suggested to facilitate nucleosome disassembly and re-assembly. In this study, by performing in vitro transcription assays and molecular simulations, we found that partial unwrapping of a nucleosome by an RNA polymerase dramatically facilitates an H2A/H2B dimer dismantling from the nucleosome by Nucleosome Assembly Protein 1 (Nap1). Furthermore, the results uncovered molecular mechanisms of Nap1 functions in which the highly acidic C-terminal flexible tails of Nap1 contribute to the H2A/H2B binding by associating with the binding interface buried and not accessible to Nap1 globular domains, supporting the penetrating fuzzy binding mechanism seemingly shared across various histone chaperones. These findings have broad implications for the mechanisms by which histone chaperones process nucleosomes upon collisions with translocases in transcription, histone recycling and nucleosomal DNA repair.

Particle Smoother to Assimilate Asynchronous Movie Data of High-Speed AFM with MD Simulations

High-speed (HS) atomic force microscopy (AFM) can be used to observe structural dynamics of biomolecules under near-physiological conditions. In the AFM measurement, the probe tip scans an area of interest and acquires height data pixel by pixel so that the obtained AFM image contains a measurement time difference. In this study, to integrate molecular dynamics simulations with asynchronous HS-AFM movie data, we developed a particle smoother (PS) method for Bayesian data assimilation, one of the machine learning approaches, by extending the previous particle filter method. With a twin experiment with an asynchronous pseudo HS-AFM movie of a nucleosome, we found that the PS method with the pixel-by-pixel data acquisition reproduced the dynamic behavior of a nucleosome better than the previous particle filter method that ignored the data asynchronicity. We examined several frequencies of particle resampling in the PS method, and found that resampling once per one frame was optimal for reproducing the dynamic behavior. Thus, we found that the PS method with an appropriate resampling frequency is a powerful method for estimating the dynamic behavior of a target molecule from HS-AFM data with low spatiotemporal resolution.

Protein dynamics by the combination of high-speed AFM and computational modeling

High-speed atomic force microscopy (HS-AFM) allows direct observation of biological molecules in dynamic action. However, HS-AFM has no atomic resolution. This article reviews recent progress of computational methods to infer high-resolution information, including the construction of 3D atomistic structures, from experimentally acquired resolution-limited HS-AFM images.

Native doublet microtubules from Tetrahymena thermophila reveal the importance of outer junction proteins

Cilia are ubiquitous eukaryotic organelles responsible for cellular motility and sensory functions. The ciliary axoneme is a microtubule-based cytoskeleton consisting of two central singlets and nine outer doublet microtubules. Cryo-electron microscopy-based studies have revealed a complex network inside the lumen of both tubules composed of microtubule-inner proteins (MIPs). However, the functions of most MIPs remain unknown. Here, we present single-particle cryo-EM-based analyses of the Tetrahymena thermophila native doublet microtubule and identify 42 MIPs. These data shed light on the evolutionarily conserved and diversified roles of MIPs. In addition, we identified MIPs potentially responsible for the assembly and stability of the doublet outer junction. Knockout of the evolutionarily conserved outer junction component CFAP77 moderately diminishes Tetrahymena swimming speed and beat frequency, indicating the important role of CFAP77 and outer junction stability in cilia beating generation and/or regulation.

BioAFMviewer software for simulation atomic force microscopy of molecular structures and conformational dynamics

Atomic force microscopy (AFM) and high-speed scanning have significantly advanced real time observation of biomolecular dynamics, with applications ranging from single molecules to the cellular level. To facilitate the interpretation of resolution-limited imaging, post-experimental computational analysis plays an increasingly important role to understand AFM measurements. Data-driven simulation of AFM, computationally emulating experimental scanning, and automatized fitting has recently elevated the understanding of measured AFM topographies by inferring the underlying full 3D atomistic structures. Providing an interactive user-friendly interface for simulation AFM, the BioAFMviewer software has become an established tool within the Bio-AFM community, with a plethora of applications demonstrating how the obtained full atomistic information advances molecular understanding beyond topographic imaging. This graphical review illustrates the BioAFMviewer capacities and further emphasizes the importance of simulation AFM to complement experimental observations.

Regulatory mechanisms of the dynein-2 motility by post-translational modification revealed by MD simulation

Intraflagellar transport for ciliary assembly and maintenance is driven by dynein and kinesins specific to the cilia. It has been shown that anterograde and retrograde transports run on different regions of the doublet microtubule, i.e., separate train tracks. However, little is known about the regulatory mechanism of this selective process. Since the doublet microtubule is known to display specific post-translational modifications of tubulins, i.e., "tubulin code", for molecular motor regulations, we investigated the motility of ciliary specific dynein-2 under different post-translational modification by coarse-grained molecular dynamics. Our setup allows us to simulate the landing behaviors of dynein-2 on un-modified, detyrosinated, poly-glutamylated and poly-glycylated microtubules in silico. Our study revealed that poly-glutamylation can play an inhibitory effect on dynein-2 motility. Our result indicates that poly-glutamylation of the B-tubule of the doublet microtubule can be used as an efficient means to target retrograde intraflagellar transport onto the A-tubule.

End-to-end differentiable blind tip reconstruction for noisy atomic force microscopy images

Observing the structural dynamics of biomolecules is vital to deepening our understanding of biomolecular functions. High-speed (HS) atomic force microscopy (AFM) is a powerful method to measure biomolecular behavior at near physiological conditions. In the AFM, measured image profiles on a molecular surface are distorted by the tip shape through the interactions between the tip and molecule. Once the tip shape is known, AFM images can be approximately deconvolved to reconstruct the surface geometry of the sample molecule. Thus, knowing the correct tip shape is an important issue in the AFM image analysis. The blind tip reconstruction (BTR) method developed by Villarrubia (J Res Natl Inst Stand Technol 102:425, 1997) is an algorithm that estimates tip shape only from AFM images using mathematical morphology operators. While the BTR works perfectly for noise-free AFM images, the algorithm is susceptible to noise. To overcome this issue, we here propose an alternative BTR method, called end-to-end differentiable BTR, based on a modern machine learning approach. In the method, we introduce a loss function including a regularization term to prevent overfitting to noise, and the tip shape is optimized with automatic differentiation and backpropagations developed in deep learning frameworks. Using noisy pseudo-AFM images of myosin V motor domain as test cases, we show that our end-to-end differentiable BTR is robust against noise in AFM images. The method can also detect a double-tip shape and deconvolve doubled molecular images. Finally, application to real HS-AFM data of myosin V walking on an actin filament shows that the method can reconstruct the accurate surface geometry of actomyosin consistent with the structural model. Our method serves as a general post-processing for reconstructing hidden molecular surfaces from any AFM images. Codes are available at https://github.com/matsunagalab/differentiable_BTR .

Efficient Multiple Domain Ligation for Proteins Using Asparaginyl Endopeptidase by Selection of Appropriate Ligation Sites Based on Steric Hindrance

Three domain fragments of a multi-domain protein, ER-60, were ligated in two short linker regions using asparaginyl endopeptidase not involving denaturation. To identify appropriate ligation sites, by selecting several potential ligation sites with fewer mutations around two short linker regions, their ligation efficiencies and the functions of the ligated ER-60s were examined experimentally. To evaluate the dependence of ligation efficiencies on the ligation sites computationally, steric hinderances around the sites for the ligation were calculated through molecular dynamics simulations. Utilizing the steric hindrance, a site-dependent ligation potential index was introduced as reproducing the experimental ligation efficiency. Referring to this index, the reconstruction of ER-60 was succeeded by the ligation of the three domains for the first time. In addition, the new ligation potential index well-worked for application to other domain ligations. Therefore, the index may serve as a more time-effective tool for multi-site ligations.

Combining Biophysical Methods for Structure-Function Analyses of RNA in Solution

RNA-level regulation by riboswitches relies on the specific binding of small metabolites to the aptamer domain to trigger substantial conformational changes that affect transcription or translation. Although several biophysical methods have been employed to study such RNAs, the utility of any one single method is limited. Hybrid approaches, therefore, are essential to better characterize these intrinsically dynamic molecules and elucidate their regulatory mechanisms driven by ligand-induced conformational changes. This chapter outlines procedures for biochemical and biophysical characterization of RNA that employs a combination of solution-based methods: isothermal titration calorimetry (ITC), small-angle X-ray scattering (SAXS), and atomic force microscopy (AFM). Collectively, these tools provide a semi-quantitative assessment of the thermodynamics associated with ligand binding and subsequent conformational changes.

Molecular dynamics simulations for the study of chromatin biology

The organization of Eukaryotic DNA into chromatin has profound implications for the processing of genetic information. In the past years, molecular dynamics (MD) simulations proved to be a powerful tool to investigate the mechanistic basis of chromatin biology. We review recent all-atom and coarse-grained MD studies revealing how the structure and dynamics of chromatin underlie its biological functions. We describe the latest method developments; the structural fluctuations of nucleosomes and the various factors affecting them; the organization of chromatin fibers, with particular emphasis on its liquid-like character; the interactions and dynamics of transcription factors on chromatin; and how chromatin organization is modulated by molecular motors acting on DNA.

Reconstruction and validation of entire virus model with complete genome from mixed resolution cryo-EM density

It is very difficult to reconstruct computationally a large biomolecular complex in its biological entirety from experimental data. The resulting atomistic model should not contain gaps structurally and it should yield stable dynamics. We, for the first time, reconstruct from the published incomplete cryo-EM density a complete MS2 virus at atomistic resolution, that is, the capsid with the genome, and validate the result by all-atom molecular dynamics with explicit water. The available experimental data includes a high resolution protein capsid and an inhomogeneously resolved genome map. For the genomic RNA, apart from 16 hairpins with atomistic resolution, the strands near the capsid's inner surface were resolved up to the nucleic backbone level, and the innermost density was completely unresolved. As a result, only 242 nucleotides (out of 3569) were positioned, while only a fragmented backbone was outlined for the rest of the genome, making a detailed model reconstruction necessary. For model reconstruction, in addition to the available atomistic structure information, we extensively used the predicted secondary structure of the genome (base pairing). The technique was based on semi-automatic building of relatively large strands of RNA with subsequent manual positioning over the traced backbone. The entire virus structure (capsid + genome) was validated by a molecular dynamics run in physiological solution with ions at standard conditions confirming the stability of the model.

Role of Repeated Conformational Transitions in Substrate Binding of Adenylate Kinase

The catalytic cycle of the enzyme adenylate kinase involves large conformational motions between open and closed states. A previous single-molecule experiment showed that substrate binding tends to accelerate both the opening and the closing rates and that a single turnover event often involves multiple rounds of conformational switching. In this work, we showed that the repeated conformational transitions of adenylate kinase are essential for the relaxation of incorrectly bound substrates into the catalytically competent conformation by combining all-atom and coarse-grained molecular simulations. In addition, free energy calculations based on all-atom and coarse-grained models demonstrated that the enzyme with incorrectly bound substrates has much a lower free energy barrier for domain opening compared to that with the correct substrate conformation, which may explain the the acceleration of the domain opening rate by substrate binding. The results of this work provide mechanistic understanding to previous experimental observations and shed light onto the interplay between conformational dynamics and enzyme catalysis.

Self-Effected Allosteric Coupling and Cooperativity in Hypoxic Response Regulation with Disordered Proteins

Hypersensitive regulation of cellular hypoxic response relies on cooperative displacement of one disordered protein (HIF-1α) by another disordered protein (CITED2) from the target in a negative feedback loop. Considering the weak intramolecule coupling in disordered proteins, the molecular mechanism of high cooperativity in the molecular displacement event remains elusive. Herein, we show that disordered proteins utilize a "self-effected allostery" mechanism to achieve high binding cooperativity. Different from the conventional allostery mechanisms shown by many structured or disordered proteins, this mechanism utilizes one part of the disordered protein as the effector to trigger the allosteric coupling and enhance the binding of the remaining part of the same disordered protein, contributing to high cooperativity of the displacement event. The conserved charge motif of CITED2 is the key determinant of the molecular displacement event by serving as the effector of allosteric coupling. Such self-effected allostery provides an efficient strategy to achieve high cooperativity in the molecular events involving disordered proteins.

Frustration and the Kinetic Repartitioning Mechanism of Substrate Inhibition in Enzyme Catalysis

Substrate inhibition, whereby enzymatic activity decreases with excess substrate after reaching a maximum turnover rate, is among the most elusive phenomena in enzymatic catalysis. Here, based on a dynamic energy landscape model, we investigate the underlying mechanism by performing molecular simulations and frustration analysis for a model enzyme adenylate kinase (AdK), which catalyzes the phosphoryl transfer reaction ATP + AMP ⇋ ADP + ADP. Intriguingly, these reveal a kinetic repartitioning mechanism of substrate inhibition, whereby excess substrate AMP suppresses the population of an energetically frustrated, but kinetically activated, catalytic pathway going through a substrate (ATP)-product (ADP) cobound complex with steric incompatibility. Such a frustrated pathway plays a crucial role in facilitating the bottleneck product ADP release, and its suppression by excess substrate AMP leads to a slow down of product release and overall turnover. The simulation results directly demonstrate that substrate inhibition arises from the rate-limiting product-release step, instead of the steps for populating the catalytically competent complex as often suggested in previous works. Furthermore, there is a tight interplay between the enzyme conformational equilibrium and the extent of substrate inhibition. Mutations biasing to more closed conformations tend to enhance substrate inhibition. We also characterized the key features of single-molecule enzyme kinetics with substrate inhibition effect. We propose that the above molecular mechanism of substrate inhibition may be relevant to other multisubstrate enzymes in which product release is the bottleneck step.

Structural dynamics of SARS-CoV-2 nucleocapsid protein induced by RNA binding

The nucleocapsid (N) protein of the SARS-CoV-2 virus, the causal agent of COVID-19, is a multifunction phosphoprotein that plays critical roles in the virus life cycle, including transcription and packaging of the viral RNA. To play such diverse roles, the N protein has two globular RNA-binding modules, the N- (NTD) and C-terminal (CTD) domains, which are connected by an intrinsically disordered region. Despite the wealth of structural data available for the isolated NTD and CTD, how these domains are arranged in the full-length protein and how the oligomerization of N influences its RNA-binding activity remains largely unclear. Herein, using experimental data from electron microscopy and biochemical/biophysical techniques combined with molecular modeling and molecular dynamics simulations, we show that, in the absence of RNA, the N protein formed structurally dynamic dimers, with the NTD and CTD arranged in extended conformations. However, in the presence of RNA, the N protein assumed a more compact conformation where the NTD and CTD are packed together. We also provided an octameric model for the full-length N bound to RNA that is consistent with electron microscopy images of the N protein in the presence of RNA. Together, our results shed new light on the dynamics and higher-order oligomeric structure of this versatile protein.

The nucleocapsid (N) protein of the SARS-CoV-2 virus plays an essential role in virus particle assembly as it specifically binds to and wraps the virus genomic RNA into a well-organized structure known as the ribonucleoprotein. Understanding how the N protein wraps around the virus RNA is critical for the development of strategies to inhibit virus assembly within host cells. One of the limitations regarding the molecular structure of the ribonucleoprotein, however, is that the N protein has several unstructured and mobile regions that preclude the resolution of its full atomic structure. Moreover, the N protein can form higher-order oligomers, both in the presence and absence of RNA. Here we employed computational methods, supported by experimental data, to simulate the N protein structural dynamics in the absence and presence of RNA. Our data suggest that the N protein forms structurally dynamic dimers in the absence of RNA, with its structured N- and C-terminal domains oriented in extended conformations. In the presence of RNA, however, the N protein assumes a more compact conformation. Our model for the oligomeric structure of the N protein bound to RNA helps to understand how N protein dimers interact to each other to form the ribonucleoprotein.

Inferring Conformational State of Myosin Motor in an Atomic Force Microscopy Image via Flexible Fitting Molecular Simulations

High-speed atomic force microscopy (HS-AFM) is a powerful technique to image the structural dynamics of biomolecules. We can obtain atomic-resolution structural information from the measured AFM image by superimposing a structural model on the image. We previously developed a flexible fitting molecular dynamics (MD) simulation method that allows for modest conformational changes when superimposed on an AFM image. In this study, for a molecular motor, myosin V (which changes its chemical state), we examined whether the conformationally distinct state in each HS-AFM image could be inferred via flexible fitting MD simulation. We first built models of myosin V bound to the actin filament in two conformational states, the “down-up” and “down-down” states. Then, for the previously obtained HS-AFM image of myosin bound to the actin filament, we performed flexible-fitting MD simulations using the two states. By comparing the fitting results, we inferred the conformational and chemical states from the AFM image.

Implementation of residue-level coarse-grained models in GENESIS for large-scale molecular dynamics simulations

Residue-level coarse-grained (CG) models have become one of the most popular tools in biomolecular simulations in the trade-off between modeling accuracy and computational efficiency. To investigate large-scale biological phenomena in molecular dynamics (MD) simulations with CG models, unified treatments of proteins and nucleic acids, as well as efficient parallel computations, are indispensable. In the GENESIS MD software, we implement several residue-level CG models, covering structure-based and context-based potentials for both well-folded biomolecules and intrinsically disordered regions. An amino acid residue in protein is represented as a single CG particle centered at the Cα atom position, while a nucleotide in RNA or DNA is modeled with three beads. Then, a single CG particle represents around ten heavy atoms in both proteins and nucleic acids. The input data in CG MD simulations are treated as GROMACS-style input files generated from a newly developed toolbox, GENESIS-CG-tool. To optimize the performance in CG MD simulations, we utilize multiple neighbor lists, each of which is attached to a different nonbonded interaction potential in the cell-linked list method. We found that random number generations for Gaussian distributions in the Langevin thermostat are one of the bottlenecks in CG MD simulations. Therefore, we parallelize the computations with message-passing-interface (MPI) to improve the performance on PC clusters or supercomputers. We simulate Herpes simplex virus (HSV) type 2 B-capsid and chromatin models containing more than 1,000 nucleosomes in GENESIS as examples of large-scale biomolecular simulations with residue-level CG models. This framework extends accessible spatial and temporal scales by multi-scale simulations to study biologically relevant phenomena, such as genome-scale chromatin folding or phase-separated membrane-less condensations.