Modeling allosteric mechanisms of eukaryotic type II topoisomerases.

Type II topoisomerases (TopoIIs) are ubiquitous enzymes that are involved in crucial nuclear processes such as genome organization, chromosome segregation, and other DNA metabolic processes. These enzymes function as large, homodimeric complexes that undergo a complex cycle of binding and hydrolysis of two ATP molecules in their ATPase domains, which regulates the capture and passage of one DNA double-helix through a second, cleaved DNA molecule. This process requires the transmission of information about the state of the bound nucleotide over vast ranges in the TopoII complex. How this information is transmitted at the molecular level to regulate TopoII functions and how protein substitutions disrupt these mechanisms remains largely unknown. Here, we employed extensive microsecond-scale molecular dynamics simulations of the yeast TopoII enzyme in multiple nucleotide-bound states and with amino acid substitutions near both the N and C termini of the complex. Simulation results indicate that the ATPase domains are remarkably flexible on the sub-microsecond timescale and that these dynamics are modulated by the identity of the bound nucleotides and both local and distant amino acid substitutions. Network analyses point toward specific allosteric networks that transmit information about the hydrolysis cycle throughout the complex, which include residues in both the protein and the bound DNA molecule. Amino acid substitutions weaken many of these pathways. Together, our results provide molecular level details on how the TopoII catalytic cycle is controlled through nucleotide binding and hydrolysis and how mutations may disrupt this process.

The effects of RNA.DNA-DNA triple helices on nucleosome structures and dynamics.

Noncoding RNAs (ncRNAs) are an emerging epigenetic factor and have been recognized as playing a key role in many gene expression pathways. Structurally, binding of ncRNAs to isolated DNA is strongly dependent on sequence complementary and results in the formation of an RNA.DNA-DNA (RDD) triple helix. However, in vivo DNA is not isolated but is rather packed in chromatin fibers, the fundamental unit of which is the nucleosome. Biochemical experiments have shown that ncRNA binding to nucleosomal DNA is elevated at DNA entry and exit sites and is dependent on the presence of the H3 N-terminal tails. However, the structural and dynamical bases for these mechanisms remain unknown. Here, we have examined the mechanisms and effects of RDD formation in the context of the nucleosome using a series of all-atom molecular dynamics simulations. Results highlight the importance of DNA sequence on complex stability, elucidate the effects of the H3 tails on RDD structures, show how RDD formation impacts the structure and dynamics of the H3 tails, and show how RNA alters the local and global DNA double-helical structure. Together, our results suggest ncRNAs can modify nucleosome, and potentially higher-order chromatin, structures and dynamics as a means of exerting epigenetic control.

Conformational and Thermodynamic Differences Underlying Wild-Type and Mutant Eleven-Nineteen-Leukemia YEATS Domain Specificity for Epigenetic Marks.

Histone post-translational modifications (PTMs) are interpreted by multiple reader domains and proteins to regulate gene expression. The eleven-nineteen-leukemia (ENL) YEATS domain is a prototypical PTM reader that recognizes multiple lysine acetylation marks on the histone H3 tails as a way of recruiting chromatin remodellers. Two ENL YEATS mutations have been identified which have been linked with leukemia, Wilms tumor, and other forms of cancer and result in either an insertion or deletion of residues in the loop connecting beta sheets distant from the protein active site. In vitro experiments have shown that these mutations modulate the selectivities of YEATS domains for various lysine acetylation marks, although different experiments have provided contrasting views on the abilities of the insertion and deletion mutants to discern specific PTMs. Here, we have performed multiple molecular dynamics simulations of wild-type and insertion and deletion mutant YEATS domains free from and in complex with two PTM peptides: one that is acetylated at K9 of H3 and the other that is acetylated at residue K27 of H3. Results show that these two peptides have distinct flexibilities and binding energetics when bound to YEATS domains and that these properties are affected by interactions with residues within and outside of the peptide consensus motif. Furthermore, these properties are modulated by the YEATS insertion and deletion mutants, which results in disparate binding effects in these systems. Together, these results suggest that only the partial exposure of histone tails is sufficient in the context of nucleosomes for YEATS-mediated recognition of acetylation marks on histone tails. They also caution against the overinterpretation of results obtained from experiments on reader domain-histone peptide binding in isolation and not in the full-length nucleosome context.

Nucleosome composition regulates the histone H3 tail conformational ensemble and accessibility.

Hexasomes and tetrasomes are intermediates in nucleosome assembly and disassembly. Their formation is promoted by histone chaperones, ATP-dependent remodelers, and RNA polymerase II. In addition, hexasomes are maintained in transcribed genes and could be an important regulatory factor. While nucleosome composition has been shown to affect the structure and accessibility of DNA, its influence on histone tails is largely unknown. Here, we investigate the conformational dynamics of the H3 tail in the hexasome and tetrasome. Using a combination of NMR spectroscopy, MD simulations, and trypsin proteolysis, we find that the conformational ensemble of the H3 tail is regulated by nucleosome composition. As has been found for the nucleosome, the H3 tails bind robustly to DNA within the hexasome and tetrasome, but upon loss of the H2A/H2B dimer, we determined that the adjacent H3 tail has an altered conformational ensemble, increase in dynamics, and increase in accessibility. Similar to observations of DNA dynamics, this is seen to be asymmetric in the hexasome. Our results indicate that nucleosome composition has the potential to regulate chromatin signaling and ultimately help shape the chromatin landscape.

Elucidating the influence of linker histone variants on chromatosome dynamics and energetics.

Linker histones are epigenetic regulators that bind to nucleosomes and alter chromatin structures and dynamics. Biophysical studies have revealed two binding modes in the linker histone/nucleosome complex, the chromatosome, where the linker histone is either centered on or askew from the dyad axis. Each has been posited to have distinct effects on chromatin, however the molecular and thermodynamic mechanisms that drive them and their dependence on linker histone compositions remain poorly understood. We present molecular dynamics simulations of chromatosomes with the globular domain of two linker histone variants, generic H1 (genGH1) and H1.0 (GH1.0), to determine how their differences influence chromatosome structures, energetics and dynamics. Results show that both unbound linker histones adopt a single compact conformation. Upon binding, DNA flexibility is reduced, resulting in increased chromatosome compaction. While both variants enthalpically favor on-dyad binding, energetic benefits are significantly higher for GH1.0, suggesting that GH1.0 is more capable than genGH1 of overcoming the large entropic reduction required for on-dyad binding which helps rationalize experiments that have consistently demonstrated GH1.0 in on-dyad states but that show genGH1 in both locations. These simulations highlight the thermodynamic basis for different linker histone binding motifs, and details their physical and chemical effects on chromatosomes.

Effects of H2A.B incorporation on nucleosome structures and dynamics.

The H2A.B histone variant is an epigenetic regulator involved in transcriptional upregulation, DNA synthesis, and splicing that functions by replacing the canonical H2A histone in the nucleosome core particle. Introduction of H2A.B results in less compact nucleosome states with increased DNA unwinding and accessibility at the nucleosomal entry and exit sites. Despite being well characterized experimentally, the molecular mechanisms by which H2A.B incorporation alters nucleosome stability and dynamics remain poorly understood. To study the molecular mechanisms of H2A.B, we have performed a series of conventional and enhanced sampling molecular dynamics simulation of H2A.B- and canonical H2A-containing nucleosomes. Results of conventional simulations show that H2A.B weakens protein-protein and protein-DNA interactions at specific locations throughout the nucleosome. These weakened interactions result in significantly more DNA opening from both the entry and exit sites in enhanced sampling simulations. Furthermore, free energy profiles show that H2A.B-containing nucleosomes have significantly broader free wells and that H2A.B allows for sampling of states with increased DNA breathing, which are shown to be stable on the hundreds of nanoseconds timescale with further conventional simulations. Together, our results show the molecular mechanisms by which H2A.B creates less compacted nucleosome states as a means of increasing genetic accessibility and gene transcription.

BEES: Bayesian Ensemble Estimation from SAS.

Many biomolecular complexes exist in a flexible ensemble of states in solution that is necessary to perform their biological function. Small-angle scattering (SAS) measurements are a popular method for characterizing these flexible molecules because of their relative ease of use and their ability to simultaneously probe the full ensemble of states. However, SAS data is typically low dimensional and difficult to interpret without the assistance of additional structural models. In theory, experimental SAS curves can be reconstituted from a linear combination of theoretical models, although this procedure carries a significant risk of overfitting the inherently low-dimensional SAS data. Previously, we developed a Bayesian-based method for fitting ensembles of model structures to experimental SAS data that rigorously avoids overfitting. However, we have found that these methods can be difficult to incorporate into typical SAS modeling workflows, especially for users that are not experts in computational modeling. To this end, we present the Bayesian Ensemble Estimation from SAS (BEES) program. Two forks of BEES are available, the primary one existing as a module for the SASSIE web server and a developmental version that is a stand-alone Python program. BEES allows users to exhaustively sample ensemble models constructed from a library of theoretical states and to interactively analyze and compare each model's performance. The fitting routine also allows for secondary data sets to be supplied, thereby simultaneously fitting models to both SAS data as well as orthogonal information. The flexible ensemble of K63-linked ubiquitin trimers is presented as an example of BEES' capabilities.

Empirical and Computational Comparison of Alternative Therapeutic Exon Skip Repairs for Duchenne Muscular Dystrophy.

Duchenne muscular dystrophy (DMD) is a common and devastating genetic disease primarily caused by exon deletions that create a genetic frameshift in dystrophin. Exon skipping therapy seeks to correct this by masking an exon during the mRNA maturation process, restoring dystrophin expression, but creating an edited protein missing both the original defect and the therapeutically skipped region. Crucially, it is possible to correct many defects in alternative ways, by skipping an exon either before or after the patient's defect. This results in alternatively edited, hybrid proteins that might have different properties and therapeutic consequences. We examined three such dystrophin exon-skipped edits, Δe45-53, Δe46-54, and Δe47-55, comprising two pairs of alternative repairs of Δe46-53 and Δe47-54 DMD defects. We found that in both cases, Δe46-54 was the more stable repair as determined by a variety of thermodynamic and biochemical measurements. We also examined the origin of these differences with molecular dynamics simulations, which showed that these stability differences were the result of different types of structural perturbations. For example, in one edit there was partial unfolding at the edit site that caused domain-localized perturbations while in another there was unfolding at the protein domain junctions distal to the edit site that increased molecular flexibility. These results demonstrate that alternative exon skip repairs of the same underlying defect can have very different consequences at the level of protein structure and stability and furthermore that these can arise by different mechanisms, either locally or by more subtle long-range perturbations.

Energetics underlying hemin extraction from human hemoglobin by Staphylococcus aureus.

Staphylococcus aureus is a leading cause of life-threatening infections in the United States. It actively acquires the essential nutrient iron from human hemoglobin (Hb) using the iron-regulated surface-determinant (Isd) system. This process is initiated when the closely related bacterial IsdB and IsdH receptors bind to Hb and extract its hemin through a conserved tri-domain unit that contains two NEAr iron Transporter (NEAT) domains that are connected by a helical linker domain. Previously, we demonstrated that the tri-domain unit within IsdH (IsdHN2N3) triggers hemin release by distorting Hb's F-helix. Here, we report that IsdHN2N3 promotes hemin release from both the α- and ß-subunits. Using a receptor mutant that only binds to the α-subunit of Hb and a stopped-flow transfer assay, we determined the energetics and micro-rate constants of hemin extraction from tetrameric Hb. We found that at 37 °C, the receptor accelerates hemin release from Hb up to 13,400-fold, with an activation enthalpy of 19.5 ± 1.1 kcal/mol. We propose that hemin removal requires the rate-limiting hydrolytic cleavage of the axial HisF8 Nϵ-Fe3+ bond, which, based on molecular dynamics simulations, may be facilitated by receptor-induced bond hydration. Isothermal titration calorimetry experiments revealed that two distinct IsdHN2N3·Hb protein·protein interfaces promote hemin release. A high-affinity receptor·Hb(A-helix) interface contributed ∼95% of the total binding standard free energy, enabling much weaker receptor interactions with Hb's F-helix that distort its hemin pocket and cause unfavorable changes in the binding enthalpy. We present a model indicating that receptor-introduced structural distortions and increased solvation underlie the IsdH-mediated hemin extraction mechanism.

Atomistic Scale Effects of Lipopolysaccharide Modifications on Bacterial Outer Membrane Defenses.

Lipopolysaccharides (LPS) are a main constituent of the outer membrane of Gram-negative bacteria. Salmonella enterica, like many other bacterial species, are able to chemically modify the structure of their LPS molecules through the PhoPQ pathway as a defense mechanism against the host immune response. These modifications make the outer membrane more resistant to antimicrobial peptides (AMPs), large lipophilic drugs, and cation depletion, and are crucial for survival within a host organism. It is believed that these LPS modifications prevent the penetration of large molecules and AMPs through a strengthening of lateral interactions between neighboring LPS molecules. Here, we performed a series of long-timescale molecular dynamics simulations to study how each of three key S. enterica lipid A modifications affect bilayer properties, with a focus on membrane structural characteristics, lateral interactions, and the divalent cation bridging network. Our results discern the unique impact each modification has on strengthening the bacterial outer membrane through effects such as increased hydrogen bonding and tighter lipid packing. Additionally, one of the modifications studied shifts Ca2+ from the lipid A region, replacing it as a major cross-linking agent between adjacent lipids and potentially making bacteria less susceptible to AMPs that competitively displace cations from the membrane surface. These results further improve our understanding of outer membrane chemical properties and help elucidate how outer membrane modification systems, such as PhoPQ in S. enterica, are able to alter bacterial virulence.

Effects of MacroH2A and H2A.Z on Nucleosome Dynamics as Elucidated by Molecular Dynamics Simulations.

Eukaryotes tune the transcriptional activity of their genome by altering the nucleosome core particle through multiple chemical processes. In particular, replacement of the canonical H2A histone with the variants macroH2A and H2A.Z has been shown to affect DNA accessibility and nucleosome stability; however, the processes by which this occurs remain poorly understood. In this study, we elucidate the molecular mechanisms of these variants with an extensive molecular dynamics study of the canonical nucleosome along with three variant-containing structures: H2A.Z, macroH2A, and an H2A mutant with macroH2A-like L1 loops. Simulation results show that variant L1 loops play a pivotal role in stabilizing DNA binding to the octamer through direct interactions, core structural rearrangements, and altered allosteric networks in the nucleosome. All variants influence dynamics; however, macroH2A-like systems have the largest effect on energetics. In addition, we provide a comprehensive analysis of allosteric networks in the nucleosome and demonstrate that variants take advantage of stronger interactions between L1 loops to propagate dynamics throughout the complex. Furthermore, we show that posttranslational modifications are enriched at key locations in these networks. Taken together, these results provide, to our knowledge, new insights into the relationship between the structure, dynamics, and function of the nucleosome core particle and chromatin fibers, and how they are influenced by chromatin remodeling factors.

Free-energy landscape and characteristic forces for the initiation of DNA unzipping.

DNA unzipping, the separation of its double helix into single strands, is crucial in modulating a host of genetic processes. Although the large-scale separation of double-stranded DNA has been studied with a variety of theoretical and experimental techniques, the minute details of the very first steps of unzipping are still unclear. Here, we use atomistic molecular-dynamics simulations, coarse-grained simulations, and a statistical-mechanical model to study the initiation of DNA unzipping by an external force. Calculation of the potential of mean force profiles for the initial separation of the first few terminal basepairs in a DNA oligomer revealed that forces ranging between 130 and 230 pN are needed to disrupt the first basepair, and these values are an order of magnitude larger than those needed to disrupt basepairs in partially unzipped DNA. The force peak has an echo of ∼50 pN at the distance that unzips the second basepair. We show that the high peak needed to initiate unzipping derives from a free-energy basin that is distinct from the basins of subsequent basepairs because of entropic contributions, and we highlight the microscopic origin of the peak. To our knowledge, our results suggest a new window of exploration for single-molecule experiments.

Structural and computational studies of the Staphylococcus aureus sortase B-substrate complex reveal a substrate-stabilized oxyanion hole.

Sortase cysteine transpeptidases covalently attach proteins to the bacterial cell wall or assemble fiber-like pili that promote bacterial adhesion. Members of this enzyme superfamily are widely distributed in Gram-positive bacteria that frequently utilize multiple sortases to elaborate their peptidoglycan. Sortases catalyze transpeptidation using a conserved active site His-Cys-Arg triad that joins a sorting signal located at the C terminus of their protein substrate to an amino nucleophile located on the cell surface. However, despite extensive study, the catalytic mechanism and molecular basis of substrate recognition remains poorly understood. Here we report the crystal structure of the Staphylococcus aureus sortase B enzyme in a covalent complex with an analog of its NPQTN sorting signal substrate, revealing the structural basis through which it displays the IsdC protein involved in heme-iron scavenging from human hemoglobin. The results of computational modeling, molecular dynamics simulations, and targeted amino acid mutagenesis indicate that the backbone amide of Glu(224) and the side chain of Arg(233) form an oxyanion hole in sortase B that stabilizes high energy tetrahedral catalytic intermediates. Surprisingly, a highly conserved threonine residue within the bound sorting signal substrate facilitates construction of the oxyanion hole by stabilizing the position of the active site arginine residue via hydrogen bonding. Molecular dynamics simulations and primary sequence conservation suggest that the sorting signal-stabilized oxyanion hole is a universal feature of enzymes within the sortase superfamily.

Nucleotide-dependent mechanism of Get3 as elucidated from free energy calculations.

The unique topology of tail-anchored (TA) proteins precludes them from utilizing the well-studied cotranslational translocation mechanism of most transmembrane proteins, forcing them into a distinct, posttranslational pathway. In yeast, this process is the guided entry of TA-proteins (GET) pathway, which utilizes a combination of cytosolic and transmembrane proteins to identify a TA protein, transfer it, and insert it into the endoplasmic reticulum membrane. At the center of this mechanism is the Get3 homodimer, which transfers a TA protein between the two GET phases by leveraging energy gained in ATP binding and hydrolysis to undergo significant structural changes from "open" to "closed" conformations. We present all-atom molecular dynamics simulations of Get3 in multiple nucleotide states, and through rigorous potential of mean force calculations, compute the free energy landscape of the Get3 opening/closing pathway. Results agree well with experiments on the nucleotide bias of Get3 open and closed structures in the crystallographically observed no-nucleotide, two ATP, and two ADP states, and also reveal their populations in the asymmetric one ATP and one ADP cases. Structures also compare well with the recently observed "semiopen" conformation and suggest that Get3 may sample this state free in solution and not just when bound to Get1, as observed in experiments. Finally, we present evidence for a unique, "wide-open" conformation of Get3. These calculations describe the nucleotide-dependent thermodynamics of Get3 in solution, and improve our understanding of its mechanism in each phase of the GET cycle.

Statistical mechanics and molecular dynamics in evaluating thermodynamic properties of biomolecular recognition.

Molecular recognition plays a central role in biochemical processes. Although well studied, understanding the mechanisms of recognition is inherently difficult due to the range of potential interactions, the molecular rearrangement associated with binding, and the time and length scales involved. Computational methods have the potential for not only complementing experiments that have been performed, but also in guiding future ones through their predictive abilities. In this review, we discuss how molecular dynamics (MD) simulations may be used in advancing our understanding of the thermodynamics that drive biomolecular recognition. We begin with a brief review of the statistical mechanics that form a basis for these methods. This is followed by a description of some of the most commonly used methods: thermodynamic pathways employing alchemical transformations and potential of mean force calculations, along with end-point calculations for free energy differences, and harmonic and quasi-harmonic analysis for entropic calculations. Finally, a few of the fundamental findings that have resulted from these methods are discussed, such as the role of configurational entropy and solvent in intermolecular interactions, along with selected results of the model system T4 lysozyme to illustrate potential and current limitations of these methods.

The Influence of Ionic Environment on Nucleosome-Mica Interactions Revealed via Molecular Dynamics Simulations.

Mica serves as a crucial substrate in Atomic Force Microscopy (AFM) studies for visualizing and characterizing nucleosomes. Nucleosomes interact with the negatively charged mica surface via adsorbed cations. However, the specific influences of monovalent and divalent cations on nucleosome adsorption to the mica surface remain unclear. In this study, we investigated the binding of nucleosomes to the mica surface in the presence of monovalent potassium ions and divalent magnesium ions using molecular dynamics simulations. We also explored the impact of pre-treated mica surfaces on nucleosome binding and structure. Our findings reveal that nucleosome-mica interactions vary depending on the cations present, resulting in distinct effects on nucleosome structure. Notably, nucleosomes bind effectively to a mica surface in the presence of potassium ions with minimal structural perturbations.

Modeling Allosteric Mechanisms of Eukaryotic Type II Topoisomerases.

Type II topoisomerases (TopoIIs) are ubiquitous enzymes that are involved in crucial nuclear processes such as genome organization, chromosome segregation, and other DNA metabolic processes. These enzymes function as large, homodimeric complexes that undergo a complex cycle of binding and hydrolysis of two ATP molecules in their ATPase domains, which regulates the capture and passage of one DNA double-helix through a second, cleaved DNA molecule. This process requires the transmission of information about the state of the bound nucleotide over vast ranges in the TopoII complex. How this information is transmitted at the molecular level to regulate TopoII functions and how protein substitutions disrupt these mechanisms remains largely unknown. Here, we employed extensive microsecond scale molecular dynamics simulations of the yeast TopoII enzyme in multiple nucleotide-bound states and with amino acid substitutions near both the N- and C-terminals of the complex. Simulation results indicate that the ATPase domains are remarkably flexible on the sub-microsecond timescale and that these dynamics are modulated by the identity of the bound nucleotides and both local and distant amino acid substitutions. Network analyses point towards specific allosteric networks that transmit information about the hydrolysis cycle throughout the complex, which include residues in both the protein and the bound DNA molecule. Amino acid substitutions weaken many of these pathways. Together, our results provide molecular-level details on how the TopoII catalytic cycle is controlled through nucleotide binding and hydrolysis and how mutations may disrupt this process.

Simulations of biased agonists in the ß(2) adrenergic receptor with accelerated molecular dynamics.

The biased agonism of the G protein-coupled receptors (GPCRs), where in addition to a traditional G protein-signaling pathway a GPCR promotes intracellular signals though ß-arrestin, is a novel paradigm in pharmacology. Biochemical and biophysical studies have suggested that a GPCR forms a distinct ensemble of conformations signaling through the G protein and ß-arrestin. Here we report on the dynamics of the ß2 adrenergic receptor bound to the ß-arrestin and G protein-biased agonists and the empty receptor to further characterize the receptor conformational changes caused by biased agonists. We use conventional and accelerated molecular dynamics (aMD) simulations to explore the conformational transitions of the GPCR from the active state to the inactive state. We found that aMD simulations enable monitoring of the transition within the nanosecond time scale while capturing the known microscopic characteristics of the inactive states, such as the ionic lock, the inward position of F6.44, and water clusters. Distinct conformational states are shown to be stabilized by each biased agonist. In particular, in simulations of the receptor with the ß-arrestin-biased agonist N-cyclopentylbutanepherine, we observe a different pattern of motions in helix 7 when compared to simulations with the G protein-biased agonist salbutamol that involves perturbations of the network of interactions within the NPxxY motif. Understanding the network of interactions induced by biased ligands and the subsequent receptor conformational shifts will lead to development of more efficient drugs.

Thermodynamic integration to predict host-guest binding affinities.

An alchemical free energy method with explicit solvent molecular dynamics simulations was applied as part of the blind prediction contest SAMPL3 to calculate binding free energies for seven guests to an acyclic cucurbit-[n]uril host. The predictions included determination of protonation states for both host and guests, docking pose generation, and binding free energy calculations using thermodynamic integration. We found a root mean square error (RMSE) of 3.6 kcal mol(-1) from the reference experimental results, with an R(2) correlation of 0.51. The agreement with experiment for the largest contributor to this error, guest 6, is improved by 1.7 kcal mol(-1) when a periodicity-induced free energy correction is applied. The corrections for the other ligands were significantly smaller, and altogether the RMSE was reduced by 0.4 kcal mol(-1). We link properties of the host-guest systems during simulation to errors in the computed free energies. Overall, we show that charged host-guest systems studied here, initialized in unconfirmed docking poses, present a challenge to accurate alchemical simulation methods.