Molecular Origin of DNA Kinking by Transcription Factors

C-B-A Test of DNA Force Fields

The DNA duplex may be locally strongly bent in complexes with proteins, for example, with polymerases or in a nucleosome. At such bends, the DNA helix is locally in the noncanonical forms A (with a narrow major groove and a large amount of north sugars) or C (with a narrow minor groove and a large share of BII phosphates). To model the formation of such complexes by molecular dynamics methods, the force field is required to reproduce these conformational transitions for a naked DNA. We analyzed the available experimental data on the B-C and B-A transitions under the conditions easily implemented in modeling: in an aqueous NaCl solution. We selected six DNA duplexes which conformations at different salt concentrations are known reliably enough. At low salt concentrations, poly(GC) and poly(A) are in the B-form, classical and slightly shifted to the A-form, respectively. The duplexes ATAT and GGTATACC have a strong and salt concentration dependent bias toward the A-form. The polymers poly(AC) and poly(G) take the C- and A-forms, respectively, at high salt concentrations. The reproduction of the behavior of these oligomers can serve as a test for the balance of interactions between the base stacking and the conformational flexibility of the sugar-phosphate backbone in a DNA force field. We tested the AMBER bsc1 and CHARMM36 force fields and their hybrids, and we failed to reproduce the experiment. In all the force fields, the salt concentration dependence is very weak. The known B-philicity of the AMBER force field proved to result from the B-philicity of its excessively strong base stacking. In the CHARMM force field, the B-form is a result of a fragile balance between the A-philic base stacking (especially for G:C pairs) and the C-philic backbone. Finally, we analyzed some recent simulations of the LacI-, SOX-4-, and Sac7d-DNA complex formation in the framework of the AMBER force field.

A Simulation Study of the Effect of Naturally Occurring Point Mutations on the SRY-DNA Complex

Molecular dynamics (MD) simulations were conducted in order to investigate the effect of the naturally occurring point mutations of the transcription factor (TF) sex-determining region Y (SRY) on the structure and dynamics of the SRY-DNA complex. The normal SRY, along with the two mutants I13T and G40R, comprising point mutations on the SRY chain, which have been clinically identified in patients with sex developmental disorders, were modeled as DNA complexes. Our modeling work aims at elucidating atomic-level structural determinants of the aberrant SRY-DNA complexation by means of μs-long MD. The results suggest that the observed disorders brought about by the G40R-DNA and I13T-DNA may arise predominantly from the destabilization of the complex being in accord with in vitro assays found elsewhere and from modifications of the DNA bending as revealed in this study. Comparative potential of mean force computations, over a sequence of short separation distances for the three complexes, verified a higher stability of the normal SRY-DNA. Examining the way the SRY mutations modulate the SRY-DNA complex dynamics at the microscopic level is important also toward elucidating molecular determinants of function for proteins capable of binding to DNA.

Doxorubicin impacts chromatin binding of HMGB1, Histone H1 and retinoic acid receptor

Doxorubicin (Dox), a widely used anticancer DNA-binding drug, affects chromatin in multiple ways, and these effects contribute to both its efficacy and its dose-limiting side effects, especially cardiotoxicity. Here, we studied the effects of Dox on the chromatin binding of the architectural proteins high mobility group B1 (HMGB1) and the linker histone H1, and the transcription factor retinoic acid receptor (RARα) by fluorescence recovery after photobleaching (FRAP) and fluorescence correlation spectroscopy (FCS) in live cells. At lower doses, Dox increased the binding of HMGB1 to DNA while decreasing the binding of the linker histone H1. At higher doses that correspond to the peak plasma concentrations achieved during chemotherapy, Dox reduced the binding of HMGB1 as well. This biphasic effect is interpreted in terms of a hierarchy of competition between the ligands involved and Dox-induced local conformational changes of nucleosome-free DNA. Combined, FRAP and FCS mobility data suggest that Dox decreases the overall binding of RARα to DNA, an effect that was only partially overcome by agonist binding. The intertwined interactions described are likely to contribute to both the effects and side effects of Dox.

Atomistic De-novo Inhibitor Generation-Guided Drug Repurposing for SARS-CoV-2 Spike Protein with Free-Energy Validation by Well-Tempered Metadynamics

Computational drug design is increasingly becoming important with new and unforeseen diseases like COVID-19. In this study, we present a new computational de novo drug design and repurposing method and applied it to find plausible drug candidates for the receptor binding domain (RBD) of SARS-CoV-2 (COVID-19). Our study comprises three steps: atom-by-atom generation of new molecules around a receptor, structural similarity mapping to existing approved and investigational drugs, and validation of their binding strengths to the viral spike proteins based on rigorous all-atom, explicit-water well-tempered metadynamics free energy calculations. By choosing the receptor binding domain of the viral spike protein, we showed that some of our new molecules and some of the repurposable drugs have stronger binding to RBD than hACE2. To validate our approach, we also calculated the free energy of hACE2 and RBD, and found it to be in an excellent agreement with experiments. These pool of drugs will allow strategic repurposing against COVID-19 for a particular prevailing conditions.

Molecular Mechanism of the Intercalation of the SOX-4 Protein into DNA Inducing Bends and Kinks

DNA-protein interactions regulate several biophysical functions, yet the mechanism of only a few is investigated in molecular detail. An important example is the intercalation of transcription factor proteins into DNA that produce bent and kinked DNA. Here, we have studied the molecular mechanism of the intercalation of a transcription factor SOX4 into DNA with a goal to understand the sequence of molecular events that precede the bending and kinking of the DNA. Our long well-tempered metadynamics and molecular dynamics (MD) simulations show that the protein primarily binds to the backbone of DNA and rotates around it to form an intercalative native state. We show that although there are multiple pathways for intercalation, the deintercalation pathway matches with the most probable intercalation pathway. In both cases, bending and kinking happen simultaneously, driven by the onset of the intercalation of the amino acid.

Mapping of scaffold/matrix attachment regions in human genome: a data mining exercise

Scaffold/matrix attachment regions (S/MARs) are DNA elements that serve to compartmentalize the chromatin into structural and functional domains. These elements are involved in control of gene expression which governs the phenotype and also plays role in disease biology. Therefore, genome-wide understanding of these elements holds great therapeutic promise. Several attempts have been made toward identification of S/MARs in genomes of various organisms including human. However, a comprehensive genome-wide map of human S/MARs is yet not available. Toward this objective, ChIP-Seq data of 14 S/MAR binding proteins were analyzed and the binding site coordinates of these proteins were used to prepare a non-redundant S/MAR dataset of human genome. Along with co-ordinate (location) details of S/MARs, the dataset also revealed details of S/MAR features, namely, length, inter-SMAR length (the chromatin loop size), nucleotide repeats, motif abundance, chromosomal distribution and genomic context. S/MARs identified in present study and their subsequent analysis also suggests that these elements act as hotspots for integration of retroviruses. Therefore, these data will help toward better understanding of genome functioning and designing effective anti-viral therapeutics. In order to facilitate user friendly browsing and retrieval of the data obtained in present study, a web interface, MARome (http://bioinfo.net.in/MARome), has been developed.

Intercalation of small molecules into DNA in chromatin is primarily controlled by superhelical constraint

The restricted access of regulatory factors to their binding sites on DNA wrapped around the nucleosomes is generally interpreted in terms of molecular shielding exerted by nucleosomal structure and internucleosomal interactions. Binding of proteins to DNA often includes intercalation of hydrophobic amino acids into the DNA. To assess the role of constrained superhelicity in limiting these interactions, we studied the binding of small molecule intercalators to chromatin in close to native conditions by laser scanning cytometry. We demonstrate that the nucleosome-constrained superhelical configuration of DNA is the main barrier to intercalation. As a result, intercalating compounds are virtually excluded from the nucleosome-occupied regions of the chromatin. Binding of intercalators to extranucleosomal regions is limited to a smaller degree, in line with the existence of net supercoiling in the regions comprising linker and nucleosome free DNA. Its relaxation by inducing as few as a single nick per ~50 kb increases intercalation in the entire chromatin loop, demonstrating the possibility for long-distance effects of regulatory potential.