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.

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.

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.

Role of bacterial RNA polymerase gate opening dynamics in DNA loading and antibiotics inhibition elucidated by quasi-Markov State Model.

To initiate transcription, the holoenzyme (RNA polymerase [RNAP] in complex with σ factor) loads the promoter DNA via the flexible loading gate created by the clamp and ß-lobe, yet their roles in DNA loading have not been characterized. We used a quasi-Markov State Model (qMSM) built from extensive molecular dynamics simulations to elucidate the dynamics of Thermus aquaticus holoenzyme's gate opening. We showed that during gate opening, ß-lobe oscillates four orders of magnitude faster than the clamp, whose opening depends on the Switch 2's structure. Myxopyronin, an antibiotic that binds to Switch 2, was shown to undergo a conformational selection mechanism to inhibit clamp opening. Importantly, we reveal a critical but undiscovered role of ß-lobe, whose opening is sufficient for DNA loading even when the clamp is partially closed. These findings open the opportunity for the development of antibiotics targeting ß-lobe of RNAP. Finally, we have shown that our qMSMs, which encode non-Markovian dynamics based on the generalized master equation formalism, hold great potential to be widely applied to study biomolecular dynamics.

Postsynaptic protein assembly in three and two dimensions studied by mesoscopic simulations.

Recently, cellular biomolecular condensates formed via phase separation have received considerable attention. While they can be formed either in cytosol (denoted as 3D) or beneath the membrane (2D), the underlying difference between the two has not been well clarified. To compare the phase behaviors in 3D and 2D, postsynaptic density (PSD) serves as a model system. PSD is a protein condensate located under the postsynaptic membrane that influences the localization of glutamate receptors and thus contributes to synaptic plasticity. Recent in vitro studies have revealed the formation of droplets of various soluble PSD proteins via liquid-liquid phase separation. However, it is unclear how these protein condensates are formed beneath the membrane and how they specifically affect the localization of glutamate receptors in the membrane. In this study, focusing on the mixture of a glutamate receptor complex, AMPAR-TARP, and a ubiquitous scaffolding protein, PSD-95, we constructed a mesoscopic model of protein-domain interactions in PSD and performed comparative molecular simulations. The results showed a sharp contrast in the phase behaviors of protein assemblies in 3D and those under the membrane (2D). A mixture of a soluble variant of the AMPAR-TARP complex and PSD-95 in the 3D system resulted in a phase-separated condensate, which was consistent with the experimental results. However, with identical domain interactions, AMPAR-TARP embedded in the membrane formed clusters with PSD-95, but did not form a stable separated phase. Thus, the cluster formation behaviors of PSD proteins in the 3D and 2D systems were distinct. The current study suggests that, more generally, stable phase separation can be more difficult to achieve in and beneath the membrane than in 3D systems.

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.

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.

Extension of the iSoLF implicit-solvent coarse-grained model for multicomponent lipid bilayers.

iSoLF is a coarse-grained (CG) model for lipid molecules with the implicit-solvent approximation used in molecular dynamics (MD) simulations of biological membranes. Using the original iSoLF (iSoLFv1), MD simulations of lipid bilayers consisting of either POPC or DPPC and these bilayers, including membrane proteins, can be performed. Here, we improve the original model, explicitly treating the electrostatic interactions between different lipid molecules and adding CG particle types. As a result, the available lipid types increase to 30. To parameterize the potential functions of the new model, we performed all-atom MD simulations of each lipid at three different temperatures using the CHARMM36 force field and the modified TIP3P model. Then, we parameterized both the bonded and non-bonded interactions to fit the area per lipid and the membrane thickness of each lipid bilayer by using the multistate Boltzmann Inversion method. The final model reproduces the area per lipid and the membrane thickness of each lipid bilayer at the three temperatures. We also examined the applicability of the new model, iSoLFv2, to simulate the phase behaviors of mixtures of DOPC and DPPC at different concentrations. The simulation results with iSoLFv2 are consistent with those using Dry Martini and Martini 3, although iSoLFv2 requires much fewer computations. iSoLFv2 has been implemented in the GENESIS MD software and is publicly available.

The lane-switch mechanism for nucleosome repositioning by DNA translocase.

Translocases such as DNA/RNA polymerases, replicative helicases, and exonucleases are involved in eukaryotic DNA transcription, replication, and repair. Since eukaryotic genomic DNA wraps around histone octamers and forms nucleosomes, translocases inevitably encounter nucleosomes. A previous study has shown that a nucleosome repositions downstream when a translocase collides with the nucleosome. However, the molecular mechanism of the downstream repositioning remains unclear. In this study, we identified the lane-switch mechanism for downstream repositioning with molecular dynamics simulations and validated it with restriction enzyme digestion assays and deep sequencing assays. In this mechanism, after a translocase unwraps nucleosomal DNA up to the site proximal to the dyad, the remaining wrapped DNA switches its binding lane to that vacated by the unwrapping, and the downstream DNA rewraps, completing downstream repositioning. This mechanism may have broad implications for transcription through nucleosomes, histone recycling, and nucleosome remodeling.

Testing mechanisms of DNA sliding by architectural DNA-binding proteins: dynamics of single wild-type and mutant protein molecules in vitro and in vivo.

Architectural DNA-binding proteins (ADBPs) are abundant constituents of eukaryotic or bacterial chromosomes that bind DNA promiscuously and function in diverse DNA reactions. They generate large conformational changes in DNA upon binding yet can slide along DNA when searching for functional binding sites. Here we investigate the mechanism by which ADBPs diffuse on DNA by single-molecule analyses of mutant proteins rationally chosen to distinguish between rotation-coupled diffusion and DNA surface sliding after transient unbinding from the groove(s). The properties of yeast Nhp6A mutant proteins, combined with molecular dynamics simulations, suggest Nhp6A switches between two binding modes: a static state, in which the HMGB domain is bound within the minor groove with the DNA highly bent, and a mobile state, where the protein is traveling along the DNA surface by means of its flexible N-terminal basic arm. The behaviors of Fis mutants, a bacterial nucleoid-associated helix-turn-helix dimer, are best explained by mobile proteins unbinding from the major groove and diffusing along the DNA surface. Nhp6A, Fis, and bacterial HU are all near exclusively associated with the chromosome, as packaged within the bacterial nucleoid, and can be modeled by three diffusion modes where HU exhibits the fastest and Fis the slowest diffusion.

Nucleosome allostery in pioneer transcription factor binding.

While recent experiments revealed that some pioneer transcription factors (TFs) can bind to their target DNA sequences inside a nucleosome, the binding dynamics of their target recognitions are poorly understood. Here we used the latest coarse-grained models and molecular dynamics simulations to study the nucleosome-binding procedure of the two pioneer TFs, Sox2 and Oct4. In the simulations for a strongly positioning nucleosome, Sox2 selected its target DNA sequence only when the target was exposed. Otherwise, Sox2 entropically bound to the dyad region nonspecifically. In contrast, Oct4 plastically bound on the nucleosome mainly in two ways. First, the two POU domains of Oct4 separately bound to the two parallel gyres of the nucleosomal DNA, supporting the previous experimental results of the partial motif recognition. Second, the POUS domain of Oct4 favored binding on the acidic patch of histones. Then, simulating the TFs binding to a genomic nucleosome, the LIN28B nucleosome, we found that the recognition of a pseudo motif by Sox2 induced the local DNA bending and shifted the population of the rotational position of the nucleosomal DNA. The redistributed DNA phase, in turn, changed the accessibility of a distant TF binding site, which consequently affected the binding probability of a second Sox2 or Oct4. These results revealed a nucleosomal DNA-mediated allosteric mechanism, through which one TF binding event can change the global conformation, and effectively regulate the binding of another TF at distant sites. Our simulations provide insights into the binding mechanism of single and multiple TFs on the nucleosome.

The stoichiometric interaction model for mesoscopic MD simulations of liquid-liquid phase separation.

Liquid-liquid phase separation (LLPS) has received considerable attention in recent years for explaining the formation of cellular biomolecular condensates. The fluidity and the complexity of their components make molecular simulation approaches indispensable for gaining structural insights. Domain-resolution mesoscopic model simulations have been explored for cases in which condensates are formed by multivalent proteins with tandem domains. One problem with this approach is that interdomain pairwise interactions cannot regulate the valency of the binding domains. To overcome this problem, we propose a new potential, the stoichiometric interaction (SI) potential. First, we verified that the SI potential maintained the valency of the interacting domains for the test systems. We then examined a well-studied LLPS model system containing tandem repeats of SH3 domains and proline-rich motifs. We found that the SI potential alone cannot reproduce the phase diagram of LLPS quantitatively. We had to combine the SI and a pairwise interaction; the former and the latter represent the specific and nonspecific interactions, respectively. Biomolecular condensates with the mixed SI and pairwise interaction exhibited fluidity, whereas those with the pairwise interaction alone showed no detectable diffusion. We also compared the phase diagrams of the systems containing different numbers of tandem domains with those obtained from the experiments and found quantitative agreement in all but one case.

Rigid-body fitting to atomic force microscopy images for inferring probe shape and biomolecular structure.

Atomic force microscopy (AFM) can visualize functional biomolecules near the physiological condition, but the observed data are limited to the surface height of specimens. Since the AFM images highly depend on the probe tip shape, for successful inference of molecular structures from the measurement, the knowledge of the probe shape is required, but is often missing. Here, we developed a method of the rigid-body fitting to AFM images, which simultaneously finds the shape of the probe tip and the placement of the molecular structure via an exhaustive search. First, we examined four similarity scores via twin-experiments for four test proteins, finding that the cosine similarity score generally worked best, whereas the pixel-RMSD and the correlation coefficient were also useful. We then applied the method to two experimental high-speed-AFM images inferring the probe shape and the molecular placement. The results suggest that the appropriate similarity score can differ between target systems. For an actin filament image, the cosine similarity apparently worked best. For an image of the flagellar protein FlhAC, we found the correlation coefficient gave better results. This difference may partly be attributed to the flexibility in the target molecule, ignored in the rigid-body fitting. The inferred tip shape and placement results can be further refined by other methods, such as the flexible fitting molecular dynamics simulations. The developed software is publicly available.

The kinetic landscape of nucleosome assembly: A coarse-grained molecular dynamics study.

The organization of nucleosomes along the Eukaryotic genome is maintained over time despite disruptive events such as replication. During this complex process, histones and DNA can form a variety of non-canonical nucleosome conformations, but their precise molecular details and roles during nucleosome assembly remain unclear. In this study, employing coarse-grained molecular dynamics simulations and Markov state modeling, we characterized the complete kinetics of nucleosome assembly. On the nucleosome-positioning 601 DNA sequence, we observe a rich transition network among various canonical and non-canonical tetrasome, hexasome, and nucleosome conformations. A low salt environment makes nucleosomes stable, but the kinetic landscape becomes more rugged, so that the system is more likely to be trapped in off-pathway partially assembled intermediates. Finally, we find that the co-operativity between DNA bending and histone association enables positioning sequence motifs to direct the assembly process, with potential implications for the dynamic organization of nucleosomes on real genomic sequences.

Modeling of DNA binding to the condensin hinge domain using molecular dynamics simulations guided by atomic force microscopy.

The condensin protein complex compacts chromatin during mitosis using its DNA-loop extrusion activity. Previous studies proposed scrunching and loop-capture models as molecular mechanisms for the loop extrusion process, both of which assume the binding of double-strand (ds) DNA to the hinge domain formed at the interface of the condensin subunits Smc2 and Smc4. However, how the hinge domain contacts dsDNA has remained unknown. Here, we conducted atomic force microscopy imaging of the budding yeast condensin holo-complex and used this data as basis for coarse-grained molecular dynamics simulations to model the hinge structure in a transient open conformation. We then simulated the dsDNA binding to open and closed hinge conformations, predicting that dsDNA binds to the outside surface when closed and to the outside and inside surfaces when open. Our simulations also suggested that the hinge can close around dsDNA bound to the inside surface. Based on these simulation results, we speculate that the conformational change of the hinge domain might be essential for the dsDNA binding regulation and play roles in condensin-mediated DNA-loop extrusion.

The HMGB chromatin protein Nhp6A can bypass obstacles when traveling on DNA.

DNA binding proteins rapidly locate their specific DNA targets through a combination of 3D and 1D diffusion mechanisms, with the 1D search involving bidirectional sliding along DNA. However, even in nucleosome-free regions, chromosomes are highly decorated with associated proteins that may block sliding. Here we investigate the ability of the abundant chromatin-associated HMGB protein Nhp6A from Saccharomyces cerevisiae to travel along DNA in the presence of other architectural DNA binding proteins using single-molecule fluorescence microscopy. We observed that 1D diffusion by Nhp6A molecules is retarded by increasing densities of the bacterial proteins Fis and HU and by Nhp6A, indicating these structurally diverse proteins impede Nhp6A mobility on DNA. However, the average travel distances were larger than the average distances between neighboring proteins, implying Nhp6A is able to bypass each of these obstacles. Together with molecular dynamics simulations, our analyses suggest two binding modes: mobile molecules that can bypass barriers as they seek out DNA targets, and near stationary molecules that are associated with neighboring proteins or preferred DNA structures. The ability of mobile Nhp6A molecules to bypass different obstacles on DNA suggests they do not block 1D searches by other DNA binding proteins.

Modeling lipid-protein interactions for coarse-grained lipid and Cα protein models.

Biological membranes that play major roles in diverse functions are composed of numerous lipids and proteins, making them an important target for coarse-grained (CG) molecular dynamics (MD) simulations. Recently, we have developed the CG implicit solvent lipid force field (iSoLF) that has a resolution compatible with the widely used Cα protein representation [D. Ugarte La Torre and S. Takada, J. Chem. Phys. 153, 205101 (2020)]. In this study, we extended it and developed a lipid-protein interaction model that allows the combination of the iSoLF and the Cα protein force field, AICG2+. The hydrophobic-hydrophilic interaction is modeled as a modified Lennard-Jones potential in which parameters were tuned partly to reproduce the experimental transfer free energy and partly based on the free energy profile normal to the membrane surface from previous all-atom MD simulations. Then, the obtained lipid-protein interaction is tested for the configuration and placement of transmembrane proteins, water-soluble proteins, and peripheral proteins, showing good agreement with prior knowledge. The interaction is generally applicable and is implemented in the publicly available software, CafeMol.

How Cytoplasmic Dynein Couples ATP Hydrolysis Cycle to Diverse Stepping Motions: Kinetic Modeling.

Cytoplasmic dynein is a two-headed molecular motor that moves to the minus end of a microtubule by ATP hydrolysis free energy. By employing its two heads (motor domains), cytoplasmic dynein exhibits various bipedal stepping motions: inchworm and hand-over-hand motions, as well as nonalternating steps of one head. However, the molecular basis to achieve such diverse stepping manners remains unclear because of the lack of an experimental method to observe stepping and the ATPase reaction of dynein simultaneously. Here, we propose a kinetic model for bipedal motions of cytoplasmic dynein and perform Gillespie Monte Carlo simulations that qualitatively reproduce most experimental data obtained to date. The model represents the status of each motor domain as five states according to conformation and nucleotide- and microtubule-binding conditions of the domain. In addition, the relative positions of the two domains were approximated by three discrete states. Accompanied by ATP hydrolysis cycles, the model dynein stochastically and processively moved forward in multiple steps via diverse pathways, including inchworm and hand-over-hand motions, similarly to experimental data. The model reproduced key experimental motility-related properties, including velocity and run length, as functions of the ATP concentration and external force, therefore providing a plausible explanation of how dynein achieves various stepping manners with explicit characterization of nucleotide states. Our model highlights the uniqueness of dynein in the coupling of ATPase with its movement during both inchworm and hand-over-hand stepping.

Coarse-grained implicit solvent lipid force field with a compatible resolution to the Cα protein representation.

Biological membranes have been prominent targets for coarse-grained (CG) molecular dynamics simulations. While minimal CG lipid models with three beads per lipid and quantitative CG lipid models with >10 beads per lipid have been well studied, in between them, CG lipid models with a compatible resolution to residue-level CG protein models are much less developed. Here, we extended a previously developed three-bead lipid model into a five-bead model and parameterized it for two phospholipids, POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) and DPPC (1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine). The developed model, iSoLF, reproduced the area per lipid, hydrophobic thickness, and phase behaviors of the target phospholipid bilayer membranes at the physiological temperature. The model POPC and DPPC membranes were in liquid and gel phases, respectively, in accordance with experiments. We further examined the spontaneous formation of a membrane bilayer, the temperature dependence of physical properties, the vesicle dynamics, and the POPC/DPPC two-component membrane dynamics of the CG lipid model, showing some promise. Once combined with standard Cα protein models, the iSoLF model will be a powerful tool to simulate large biological membrane systems made of lipids and proteins.

DNA sliding in nucleosomes via twist defect propagation revealed by molecular simulations.

While nucleosomes are highly stable structures as fundamental units of chromatin, they also slide along the DNA, either spontaneously or by active remodelers. Here, we investigate the microscopic mechanisms of nucleosome sliding by multiscale molecular simulations, characterizing how the screw-like motion of DNA proceeds via the formation and propagation of twist defects. Firstly, coarse-grained molecular simulations reveal that the sliding dynamics is highly dependent on DNA sequence. Depending on the sequence and the nucleosome super-helical location, we find two distinct types of twist defects: a locally under-twisted DNA region, previously observed in crystal structures, and a locally over-twisted DNA, an unprecedented feature. The stability of the over-twist defect was confirmed via all-atom simulations. Analysis of our trajectories via Markov state modeling highlights how the sequence-dependence of the sliding dynamics is due to the different twist defect energy costs, and in particular how nucleosome regions where defects cannot easily form introduce the kinetic bottlenecks slowing down repositioning. Twist defects can also mediate sliding of nucleosomes made with strong positioning sequences, albeit at a much lower diffusion coefficient, due to a high-energy intermediate state. Finally, we discuss how chromatin remodelers may exploit these spontaneous fluctuations to induce unidirectional sliding of nucleosomes.