1
|
Farré-Gil D, Arcon JP, Laughton CA, Orozco M. CGeNArate: a sequence-dependent coarse-grained model of DNA for accurate atomistic MD simulations of kb-long duplexes. Nucleic Acids Res 2024; 52:6791-6801. [PMID: 38813824 PMCID: PMC11229373 DOI: 10.1093/nar/gkae444] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2023] [Revised: 05/01/2024] [Accepted: 05/14/2024] [Indexed: 05/31/2024] Open
Abstract
We present CGeNArate, a new model for molecular dynamics simulations of very long segments of B-DNA in the context of biotechnological or chromatin studies. The developed method uses a coarse-grained Hamiltonian with trajectories that are back-mapped to the atomistic resolution level with extreme accuracy by means of Machine Learning Approaches. The method is sequence-dependent and reproduces very well not only local, but also global physical properties of DNA. The efficiency of the method allows us to recover with a reduced computational effort high-quality atomic-resolution ensembles of segments containing many kilobases of DNA, entering into the gene range or even the entire DNA of certain cellular organelles.
Collapse
Affiliation(s)
- David Farré-Gil
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, Baldiri Reixac 10-12, E-08028 Barcelona, Spain
| | - Juan Pablo Arcon
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, Baldiri Reixac 10-12, E-08028 Barcelona, Spain
| | - Charles A Laughton
- School of Pharmacy and Biodiscovery Institute, University of Nottingham, University Park, Nottingham NG7 2RD, UK
| | - Modesto Orozco
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, Baldiri Reixac 10-12, E-08028 Barcelona, Spain
- Department of Biochemistry and Biomedicine, University of Barcelona, E-08028 Barcelona, Spain
| |
Collapse
|
2
|
Roldán-Piñero C, Luengo-Márquez J, Assenza S, Pérez R. Systematic Comparison of Atomistic Force Fields for the Mechanical Properties of Double-Stranded DNA. J Chem Theory Comput 2024; 20:2261-2272. [PMID: 38411091 PMCID: PMC10938644 DOI: 10.1021/acs.jctc.3c01089] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/02/2023] [Revised: 02/14/2024] [Accepted: 02/14/2024] [Indexed: 02/28/2024]
Abstract
The response of double-stranded DNA to external mechanical stress plays a central role in its interactions with the protein machinery in the cell. Modern atomistic force fields have been shown to provide highly accurate predictions for the fine structural features of the duplex. In contrast, and despite their pivotal function, less attention has been devoted to the accuracy of the prediction of the elastic parameters. Several reports have addressed the flexibility of double-stranded DNA via all-atom molecular dynamics, yet the collected information is insufficient to have a clear understanding of the relative performance of the various force fields. In this work, we fill this gap by performing a systematic study in which several systems, characterized by different sequence contexts, are simulated with the most popular force fields within the AMBER family, bcs1 and OL15, as well as with CHARMM36. Analysis of our results, together with their comparison with previous work focused on bsc0, allows us to unveil the differences in the predicted rigidity between the newest force fields and suggests a roadmap to test their performance against experiments. In the case of the stretch modulus, we reconcile these differences, showing that a single mapping between sequence-dependent conformation and elasticity via the crookedness parameter captures simultaneously the results of all force fields, supporting the key role of crookedness in the mechanical response of double-stranded DNA.
Collapse
Affiliation(s)
- Carlos Roldán-Piñero
- Departamento
de Física Teórica de la Materia Condensada, Universidad Autónoma de Madrid, E-28049 Madrid, Spain
| | - Juan Luengo-Márquez
- Departamento
de Física Teórica de la Materia Condensada, Universidad Autónoma de Madrid, E-28049 Madrid, Spain
- Instituto
Nicolás Cabrera, Universidad Autónoma
de Madrid, E-28049 Madrid, Spain
| | - Salvatore Assenza
- Departamento
de Física Teórica de la Materia Condensada, Universidad Autónoma de Madrid, E-28049 Madrid, Spain
- Instituto
Nicolás Cabrera, Universidad Autónoma
de Madrid, E-28049 Madrid, Spain
- Condensed
Matter Physics Center (IFIMAC), Universidad
Autónoma de Madrid, E-28049 Madrid, Spain
| | - Rubén Pérez
- Departamento
de Física Teórica de la Materia Condensada, Universidad Autónoma de Madrid, E-28049 Madrid, Spain
- Condensed
Matter Physics Center (IFIMAC), Universidad
Autónoma de Madrid, E-28049 Madrid, Spain
| |
Collapse
|
3
|
Coshic K, Maffeo C, Winogradoff D, Aksimentiev A. The structure and physical properties of a packaged bacteriophage particle. Nature 2024; 627:905-914. [PMID: 38448589 PMCID: PMC11196859 DOI: 10.1038/s41586-024-07150-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2023] [Accepted: 02/01/2024] [Indexed: 03/08/2024]
Abstract
A string of nucleotides confined within a protein capsid contains all the instructions necessary to make a functional virus particle, a virion. Although the structure of the protein capsid is known for many virus species1,2, the three-dimensional organization of viral genomes has mostly eluded experimental probes3,4. Here we report all-atom structural models of an HK97 virion5, including its entire 39,732 base pair genome, obtained through multiresolution simulations. Mimicking the action of a packaging motor6, the genome was gradually loaded into the capsid. The structure of the packaged capsid was then refined through simulations of increasing resolution, which produced a 26 million atom model of the complete virion, including water and ions confined within the capsid. DNA packaging occurs through a loop extrusion mechanism7 that produces globally different configurations of the packaged genome and gives each viral particle individual traits. Multiple microsecond-long all-atom simulations characterized the effect of the packaged genome on capsid structure, internal pressure, electrostatics and diffusion of water, ions and DNA, and revealed the structural imprints of the capsid onto the genome. Our approach can be generalized to obtain complete all-atom structural models of other virus species, thereby potentially revealing new drug targets at the genome-capsid interface.
Collapse
Affiliation(s)
- Kush Coshic
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL, USA
| | - Christopher Maffeo
- Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, USA
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, IL, USA
| | - David Winogradoff
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, IL, USA
| | - Aleksei Aksimentiev
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL, USA.
- Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, USA.
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, IL, USA.
| |
Collapse
|
4
|
Chen YT, Yang H, Chu JW. Mechanical codes of chemical-scale specificity in DNA motifs. Chem Sci 2023; 14:10155-10166. [PMID: 37772098 PMCID: PMC10529945 DOI: 10.1039/d3sc01671d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2023] [Accepted: 08/28/2023] [Indexed: 09/30/2023] Open
Abstract
In gene transcription, certain sequences of double-stranded (ds)DNA play a vital role in nucleosome positioning and expression initiation. That dsDNA is deformed to various extents in these processes leads us to ask: Could the genomic DNA also have sequence specificity in its chemical-scale mechanical properties? We approach this question using statistical machine learning to determine the rigidity between DNA chemical moieties. What emerges for the polyA, polyG, TpA, and CpG sequences studied here is a unique trigram that contains the quantitative mechanical strengths between bases and along the backbone. In a way, such a sequence-dependent trigram could be viewed as a DNA mechanical code. Interestingly, we discover a compensatory competition between the axial base-stacking interaction and the transverse base-pairing interaction, and such a reciprocal relationship constitutes the most discriminating feature of the mechanical code. Our results also provide chemical-scale understanding for experimental observables. For example, the long polyA persistence length is shown to have strong base stacking while its complement (polyAc) exhibits high backbone rigidity. The mechanical code concept enables a direct reading of the physical interactions encoded in the sequence which, with further development, is expected to shed new light on DNA allostery and DNA-binding drugs.
Collapse
Affiliation(s)
- Yi-Tsao Chen
- Institute of Bioinformatics and Systems Biology, National Yang Ming Chiao Tung University Hsinchu 30010 Taiwan Republic of China
| | - Haw Yang
- Department of Chemistry, Princeton University Princeton NJ 08544 USA
| | - Jhih-Wei Chu
- Institute of Bioinformatics and Systems Biology, National Yang Ming Chiao Tung University Hsinchu 30010 Taiwan Republic of China
- Department of Biological Science and Technology, Institute of Molecular Medicine and Bioengineering, Center for Intelligent Drug Systems and Smart Bio-devices (IDS2B), National Yang Ming Chiao Tung University Hsinchu 30010 Taiwan Republic of China
| |
Collapse
|
5
|
Hörberg J, Reymer A. Decoding the dual recognition mechanism of the glucocorticoid receptor for DNA and RNA: sequence versus shape. Sci Rep 2023; 13:16125. [PMID: 37752333 PMCID: PMC10522765 DOI: 10.1038/s41598-023-43244-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2023] [Accepted: 09/21/2023] [Indexed: 09/28/2023] Open
Abstract
Transcription factors (TFs) regulate eukaryotic transcription through selective DNA-binding, can also specifically interact with RNA, which may present another layer of transcriptional control. The mechanisms of the TFs-DNA recognition are often well-characterised, while the details of TFs-RNA complexation are less understood. Here we investigate the dual recognition mechanism of the glucocorticoid receptor (GR), which interacts with similar affinities with consensus DNA and diverse RNA hairpin motifs but discriminates against uniform dsRNA. Using atomic molecular dynamics simulations, we demonstrate that the GR binding to nucleic acids requires a wide and shallow groove pocket. The protein effectively moulds its binding site within DNA major groove, which enables base-specific interactions. Contrary, the GR binding has little effect on the grooves geometry of RNA systems, most notably in uniform dsRNA. Instead, a hairpin motif in RNA yields a wide and shallow major groove pocket, allowing the protein to anchor itself through nonspecific electrostatic contacts with RNA backbone. Addition of a bulge increases RNA hairpin flexibility, which leads to a greater number of GR-RNA contacts and, thus, higher affinity. Thus, the combination of structural motifs defines the GR-RNA selective binding: a recognition mechanism, which may be shared by other zinc finger TFs.
Collapse
Affiliation(s)
- Johanna Hörberg
- Department of Chemistry and Molecular Biology, University of Gothenburg, 405 30, Göteborg, Sweden
| | - Anna Reymer
- Department of Chemistry and Molecular Biology, University of Gothenburg, 405 30, Göteborg, Sweden.
| |
Collapse
|
6
|
Battistini F, Sala A, Hospital A, Orozco M. Sequence-Dependent Properties of the RNA Duplex. J Chem Inf Model 2023; 63:5259-5271. [PMID: 37577978 DOI: 10.1021/acs.jcim.3c00741] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/15/2023]
Abstract
Sequence-dependent properties of the DNA duplex have been accurately described using extensive molecular dynamics simulations. The RNA duplex meanwhile─which is typically represented as a sequence-averaged rigid rod─does not benefit from having equivalent molecular dynamics simulations. In this paper, we present a massive simulation effort using a set of ABC-optimized duplexes from which we derived tetramer-resolution properties of the RNA duplex and a simple mesoscopic model that can represent elastic properties of long RNA duplexes. Despite the extreme chemical similarity between DNA and RNA, the local and global elastic properties of the duplexes are very different. DNA duplexes show a complex and nonelastic pattern of flexibility, for instance, while RNA duplexes behave as an elastic system whose deformations can be represented by simple harmonic potentials. In RNA duplexes (RNA2), not only are intra- and interbase pair parameters (equilibrium and mechanical) different from those in the equivalent DNA duplex sequences (DNA2) but the correlations between movements also differ. Simple statements on the relative flexibility or stability of both polymers are meaningless and should be substituted by a more detailed description depending on the sequence and the type of deformation considered.
Collapse
Affiliation(s)
- Federica Battistini
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology (BIST), Baldiri Reixac 10, Barcelona 08028, Spain
- Departament de Bioquímica i Biomedicina. Facultat de Biologia, Universitat de Barcelona, Avgda Diagonal 647, Barcelona 08028, Spain
| | - Alba Sala
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology (BIST), Baldiri Reixac 10, Barcelona 08028, Spain
| | - Adam Hospital
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology (BIST), Baldiri Reixac 10, Barcelona 08028, Spain
| | - Modesto Orozco
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology (BIST), Baldiri Reixac 10, Barcelona 08028, Spain
- Departament de Bioquímica i Biomedicina. Facultat de Biologia, Universitat de Barcelona, Avgda Diagonal 647, Barcelona 08028, Spain
| |
Collapse
|
7
|
Ashwood B, Jones MS, Radakovic A, Khanna S, Lee Y, Sachleben JR, Szostak JW, Ferguson AL, Tokmakoff A. Thermodynamics and kinetics of DNA and RNA dinucleotide hybridization to gaps and overhangs. Biophys J 2023; 122:3323-3339. [PMID: 37469144 PMCID: PMC10465710 DOI: 10.1016/j.bpj.2023.07.009] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2023] [Revised: 06/27/2023] [Accepted: 07/17/2023] [Indexed: 07/21/2023] Open
Abstract
Hybridization of short nucleic acid segments (<4 nt) to single-strand templates occurs as a critical intermediate in processes such as nonenzymatic nucleic acid replication and toehold-mediated strand displacement. These templates often contain adjacent duplex segments that stabilize base pairing with single-strand gaps or overhangs, but the thermodynamics and kinetics of hybridization in such contexts are poorly understood because of the experimental challenges of probing weak binding and rapid structural dynamics. Here we develop an approach to directly measure the thermodynamics and kinetics of DNA and RNA dinucleotide dehybridization using steady-state and temperature-jump infrared spectroscopy. Our results suggest that dinucleotide binding is stabilized through coaxial stacking interactions with the adjacent duplex segments as well as from potential noncanonical base-pairing configurations and structural dynamics of gap and overhang templates revealed using molecular dynamics simulations. We measure timescales for dissociation ranging from 0.2-40 μs depending on the template and temperature. Dinucleotide hybridization and dehybridization involve a significant free energy barrier with characteristics resembling that of canonical oligonucleotides. Together, our work provides an initial step for predicting the stability and kinetics of hybridization between short nucleic acid segments and various templates.
Collapse
Affiliation(s)
- Brennan Ashwood
- Department of Chemistry, The University of Chicago, Chicago, Illinois; The James Franck Institute and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois
| | - Michael S Jones
- Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois
| | | | - Smayan Khanna
- Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois
| | - Yumin Lee
- Department of Chemistry, The University of Chicago, Chicago, Illinois; The James Franck Institute and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois
| | - Joseph R Sachleben
- Biomolecular NMR Core Facility, Biological Sciences Division, The University of Chicago, Chicago, Illinois
| | - Jack W Szostak
- Department of Chemistry, The University of Chicago, Chicago, Illinois
| | - Andrew L Ferguson
- Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois
| | - Andrei Tokmakoff
- Department of Chemistry, The University of Chicago, Chicago, Illinois; The James Franck Institute and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois.
| |
Collapse
|
8
|
Ashwood B, Jones MS, Radakovic A, Khanna S, Lee Y, Sachleben JR, Szostak JW, Ferguson AL, Tokmakoff A. Direct monitoring of the thermodynamics and kinetics of DNA and RNA dinucleotide dehybridization from gaps and overhangs. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.04.10.536266. [PMID: 37090657 PMCID: PMC10120721 DOI: 10.1101/2023.04.10.536266] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/25/2023]
Abstract
Hybridization of short nucleic acid segments (<4 nucleotides) to single-strand templates occurs as a critical intermediate in processes such as non-enzymatic nucleic acid replication and toehold-mediated strand displacement. These templates often contain adjacent duplex segments that stabilize base pairing with single-strand gaps or overhangs, but the thermodynamics and kinetics of hybridization in such contexts are poorly understood due to experimental challenges of probing weak binding and rapid structural dynamics. Here we develop an approach to directly measure the thermodynamics and kinetics of DNA and RNA dinucleotide dehybridization using steady-state and temperature-jump infrared spectroscopy. Our results suggest that dinucleotide binding is stabilized through coaxial stacking interactions with the adjacent duplex segments as well as from potential non-canonical base pairing configurations and structural dynamics of gap and overhang templates revealed using molecular dynamics simulations. We measure timescales for dissociation ranging from 0.2 to 40 µs depending on the template and temperature. Dinucleotide hybridization and dehybridization involves a significant free energy barrier with characteristics resembling that of canonical oligonucleotides. Together, our work provides an initial step for predicting the stability and kinetics of hybridization between short nucleic acid segments and various templates.
Collapse
Affiliation(s)
- Brennan Ashwood
- Department of Chemistry, The University of Chicago, 5735 S. Ellis Avenue, Chicago, IL 60637
- The James Franck Institute and Institute for Biophysical Dynamics, The University of Chicago, 929 East 57 Street, Chicago, Illinois 60637, United States
| | - Michael S Jones
- Pritzker School of Molecular Engineering, The University of Chicago, 5640 South Ellis Avenue, Chicago, Illinois 60637, United States
| | | | - Smayan Khanna
- Pritzker School of Molecular Engineering, The University of Chicago, 5640 South Ellis Avenue, Chicago, Illinois 60637, United States
| | - Yumin Lee
- Department of Chemistry, The University of Chicago, 5735 S. Ellis Avenue, Chicago, IL 60637
| | - Joseph R Sachleben
- Biomolecular NMR Core Facility, Biological Sciences Division, The University of Chicago, Chicago, IL 60637, United States
| | - Jack W Szostak
- Department of Chemistry, The University of Chicago, 5735 S. Ellis Avenue, Chicago, IL 60637
| | - Andrew L Ferguson
- Pritzker School of Molecular Engineering, The University of Chicago, 5640 South Ellis Avenue, Chicago, Illinois 60637, United States
| | - Andrei Tokmakoff
- Department of Chemistry, The University of Chicago, 5735 S. Ellis Avenue, Chicago, IL 60637
- The James Franck Institute and Institute for Biophysical Dynamics, The University of Chicago, 929 East 57 Street, Chicago, Illinois 60637, United States
| |
Collapse
|
9
|
Sharma D, Sharma K, Mishra A, Siwach P, Mittal A, Jayaram B. Molecular dynamics simulation-based trinucleotide and tetranucleotide level structural and energy characterization of the functional units of genomic DNA. Phys Chem Chem Phys 2023; 25:7323-7337. [PMID: 36825435 DOI: 10.1039/d2cp04820e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/12/2023]
Abstract
Genomes of most organisms on earth are written in a universal language of life, made up of four units - adenine (A), thymine (T), guanine (G), and cytosine (C), and understanding the way they are put together has been a great challenge to date. Multiple efforts have been made to annotate this wonderfully engineered string of DNA using different methods but they lack a universal character. In this article, we have investigated the structural and energetic profiles of both prokaryotes and eukaryotes by considering two essential genomic sites, viz., the transcription start sites (TSS) and exon-intron boundaries. We have characterized these sites by mapping the structural and energy features of DNA obtained from molecular dynamics simulations, which considers all possible trinucleotide and tetranucleotide steps. For DNA, these physicochemical properties show distinct signatures at the TSS and intron-exon boundaries. Our results firmly convey the idea that DNA uses the same dialect for prokaryotes and eukaryotes and that it is worth going beyond sequence-level analyses to physicochemical space to determine the functional destiny of DNA sequences.
Collapse
Affiliation(s)
- Dinesh Sharma
- Supercomputing Facility for Bioinformatics & Computational Biology, Kusuma School of Biological Sciences, Indian Institute of Technology, Delhi, India
| | - Kopal Sharma
- Supercomputing Facility for Bioinformatics & Computational Biology, Kusuma School of Biological Sciences, Indian Institute of Technology, Delhi, India
| | - Akhilesh Mishra
- Supercomputing Facility for Bioinformatics & Computational Biology, Kusuma School of Biological Sciences, Indian Institute of Technology, Delhi, India
| | - Priyanka Siwach
- Department of Biotechnology, Chaudhary Devi Lal University, Sirsa, Haryana, India
| | - Aditya Mittal
- Supercomputing Facility for Bioinformatics & Computational Biology, Kusuma School of Biological Sciences, Indian Institute of Technology, Delhi, India
| | - B Jayaram
- Supercomputing Facility for Bioinformatics & Computational Biology, Kusuma School of Biological Sciences, Indian Institute of Technology, Delhi, India.,Department of Chemistry, Indian Institute of Technology, Delhi, India.
| |
Collapse
|
10
|
Basu A, Bobrovnikov DG, Cieza B, Arcon JP, Qureshi Z, Orozco M, Ha T. Deciphering the mechanical code of the genome and epigenome. Nat Struct Mol Biol 2022; 29:1178-1187. [PMID: 36471057 PMCID: PMC10142808 DOI: 10.1038/s41594-022-00877-6] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2021] [Accepted: 10/18/2022] [Indexed: 12/12/2022]
Abstract
Diverse DNA-deforming processes are impacted by the local mechanical and structural properties of DNA, which in turn depend on local sequence and epigenetic modifications. Deciphering this mechanical code (that is, this dependence) has been challenging due to the lack of high-throughput experimental methods. Here we present a comprehensive characterization of the mechanical code. Utilizing high-throughput measurements of DNA bendability via loop-seq, we quantitatively established how the occurrence and spatial distribution of dinucleotides, tetranucleotides and methylated CpG impact DNA bendability. We used our measurements to develop a physical model for the sequence and methylation dependence of DNA bendability. We validated the model by performing loop-seq on mouse genomic sequences around transcription start sites and CTCF-binding sites. We applied our model to test the predictions of all-atom molecular dynamics simulations and to demonstrate that sequence and epigenetic modifications can mechanically encode regulatory information in diverse contexts.
Collapse
Affiliation(s)
- Aakash Basu
- Department of Biosciences, Durham University, Durham, UK. .,Department of Biophysics and Biophysical Chemistry, Johns Hopkins University School of Medicine, Baltimore, MD, USA.
| | - Dmitriy G Bobrovnikov
- Department of Biophysics and Biophysical Chemistry, Johns Hopkins University School of Medicine, Baltimore, MD, USA
| | - Basilio Cieza
- Department of Biophysics, Johns Hopkins University, Baltimore, MD, USA
| | - Juan Pablo Arcon
- Institute for Research in Biomedicine (IRB Barcelona), Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Zan Qureshi
- Department of Biophysics, Johns Hopkins University, Baltimore, MD, USA
| | - Modesto Orozco
- Institute for Research in Biomedicine (IRB Barcelona), Barcelona Institute of Science and Technology, Barcelona, Spain.,Department of Biochemistry and Biomedicine, Universitat de Barcelona, Barcelona, Spain
| | - Taekjip Ha
- Department of Biophysics and Biophysical Chemistry, Johns Hopkins University School of Medicine, Baltimore, MD, USA. .,Department of Biophysics, Johns Hopkins University, Baltimore, MD, USA. .,Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, USA. .,Howard Hughes Medical Institute, Baltimore, MD, USA.
| |
Collapse
|
11
|
Neguembor MV, Arcon JP, Buitrago D, Lema R, Walther J, Garate X, Martin L, Romero P, AlHaj Abed J, Gut M, Blanc J, Lakadamyali M, Wu CT, Brun Heath I, Orozco M, Dans PD, Cosma MP. MiOS, an integrated imaging and computational strategy to model gene folding with nucleosome resolution. Nat Struct Mol Biol 2022; 29:1011-1023. [PMID: 36220894 PMCID: PMC9627188 DOI: 10.1038/s41594-022-00839-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2022] [Accepted: 08/23/2022] [Indexed: 11/06/2022]
Abstract
The linear sequence of DNA provides invaluable information about genes and their regulatory elements along chromosomes. However, to fully understand gene function and regulation, we need to dissect how genes physically fold in the three-dimensional nuclear space. Here we describe immuno-OligoSTORM, an imaging strategy that reveals the distribution of nucleosomes within specific genes in super-resolution, through the simultaneous visualization of DNA and histones. We combine immuno-OligoSTORM with restraint-based and coarse-grained modeling approaches to integrate super-resolution imaging data with Hi-C contact frequencies and deconvoluted micrococcal nuclease-sequencing information. The resulting method, called Modeling immuno-OligoSTORM, allows quantitative modeling of genes with nucleosome resolution and provides information about chromatin accessibility for regulatory factors, such as RNA polymerase II. With Modeling immuno-OligoSTORM, we explore intercellular variability, transcriptional-dependent gene conformation, and folding of housekeeping and pluripotency-related genes in human pluripotent and differentiated cells, thereby obtaining the highest degree of data integration achieved so far to our knowledge.
Collapse
Affiliation(s)
- Maria Victoria Neguembor
- Center for Genomic Regulation (CRG), Barcelona Institute of Science and Technology, Barcelona, Spain.
| | - Juan Pablo Arcon
- Institute for Research in Biomedicine (IRB Barcelona), Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Diana Buitrago
- Institute for Research in Biomedicine (IRB Barcelona), Barcelona Institute of Science and Technology, Barcelona, Spain
- Departamento de Física y Matemáticas, Universidad Autónoma de Manizales, Manizales, Colombia
| | - Rafael Lema
- Institute for Research in Biomedicine (IRB Barcelona), Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Jürgen Walther
- Institute for Research in Biomedicine (IRB Barcelona), Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Ximena Garate
- Center for Genomic Regulation (CRG), Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Laura Martin
- Center for Genomic Regulation (CRG), Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Pablo Romero
- Institute for Research in Biomedicine (IRB Barcelona), Barcelona Institute of Science and Technology, Barcelona, Spain
| | | | - Marta Gut
- CNAG-CRG, Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology, Barcelona, Spain
- Universitat Pompeu Fabra (UPF), Barcelona, Spain
| | - Julie Blanc
- CNAG-CRG, Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Melike Lakadamyali
- Department of Physiology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - Chao-Ting Wu
- Department of Genetics, Harvard Medical School, Boston, MA, USA
| | - Isabelle Brun Heath
- Institute for Research in Biomedicine (IRB Barcelona), Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Modesto Orozco
- Institute for Research in Biomedicine (IRB Barcelona), Barcelona Institute of Science and Technology, Barcelona, Spain.
- Faculty of Biology, University of Barcelona, Barcelona, Spain.
- ICREA, Barcelona, Spain.
| | - Pablo D Dans
- Institute for Research in Biomedicine (IRB Barcelona), Barcelona Institute of Science and Technology, Barcelona, Spain.
- Department of Biological Sciences, CENUR Litoral Norte, Universidad de la República (UdelaR), Salto, Uruguay.
- Bioinformatics Unit, Institut Pasteur de Montevideo, Montevideo, Uruguay.
| | - Maria Pia Cosma
- Center for Genomic Regulation (CRG), Barcelona Institute of Science and Technology, Barcelona, Spain.
- Universitat Pompeu Fabra (UPF), Barcelona, Spain.
- ICREA, Barcelona, Spain.
- CAS Key Laboratory of Regenerative Biology, Guangdong Provincial Key Laboratory of Stem Cell and Regenerative Medicine, Guangzhou Institutes of Biomedicine and Health, Chinese Academy of Sciences, Guangzhou, China.
| |
Collapse
|
12
|
Barissi S, Sala A, Wieczór M, Battistini F, Orozco M. DNAffinity: a machine-learning approach to predict DNA binding affinities of transcription factors. Nucleic Acids Res 2022; 50:9105-9114. [PMID: 36018808 PMCID: PMC9458447 DOI: 10.1093/nar/gkac708] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2022] [Revised: 07/21/2022] [Accepted: 08/08/2022] [Indexed: 12/24/2022] Open
Abstract
We present a physics-based machine learning approach to predict in vitro transcription factor binding affinities from structural and mechanical DNA properties directly derived from atomistic molecular dynamics simulations. The method is able to predict affinities obtained with techniques as different as uPBM, gcPBM and HT-SELEX with an excellent performance, much better than existing algorithms. Due to its nature, the method can be extended to epigenetic variants, mismatches, mutations, or any non-coding nucleobases. When complemented with chromatin structure information, our in vitro trained method provides also good estimates of in vivo binding sites in yeast.
Collapse
Affiliation(s)
| | | | - Miłosz Wieczór
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and Technology. Baldiri Reixac 10–12, 08028 Barcelona, Spain,Department of Physical Chemistry. Gdansk University of Technology, 80-233 Gdańsk, Poland
| | | | - Modesto Orozco
- Correspondence may also be addressed to Modesto Orozco. Tel: +34 934 037 156;
| |
Collapse
|
13
|
Wieczór M, Genna V, Aranda J, Badia RM, Gelpí JL, Gapsys V, de Groot BL, Lindahl E, Municoy M, Hospital A, Orozco M. Pre-exascale HPC approaches for molecular dynamics simulations. Covid-19 research: A use case. WILEY INTERDISCIPLINARY REVIEWS. COMPUTATIONAL MOLECULAR SCIENCE 2022; 13:e1622. [PMID: 35935573 PMCID: PMC9347456 DOI: 10.1002/wcms.1622] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 02/02/2022] [Revised: 04/25/2022] [Accepted: 04/28/2022] [Indexed: 06/15/2023]
Abstract
Exascale computing has been a dream for ages and is close to becoming a reality that will impact how molecular simulations are being performed, as well as the quantity and quality of the information derived for them. We review how the biomolecular simulations field is anticipating these new architectures, making emphasis on recent work from groups in the BioExcel Center of Excellence for High Performance Computing. We exemplified the power of these simulation strategies with the work done by the HPC simulation community to fight Covid-19 pandemics. This article is categorized under:Data Science > Computer Algorithms and ProgrammingData Science > Databases and Expert SystemsMolecular and Statistical Mechanics > Molecular Dynamics and Monte-Carlo Methods.
Collapse
Affiliation(s)
- Miłosz Wieczór
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and TechnologyBarcelonaSpain
- Department of Physical ChemistryGdansk University of TechnologyGdańskPoland
| | - Vito Genna
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and TechnologyBarcelonaSpain
| | - Juan Aranda
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and TechnologyBarcelonaSpain
| | | | - Josep Lluís Gelpí
- Barcelona Supercomputing CenterBarcelonaSpain
- Department of Biochemistry and BiomedicineUniversity of BarcelonaBarcelonaSpain
| | - Vytautas Gapsys
- Max Planck Institute for Multidisciplinary SciencesComputational Biomolecular Dynamics GroupGoettingenGermany
| | - Bert L. de Groot
- Max Planck Institute for Multidisciplinary SciencesComputational Biomolecular Dynamics GroupGoettingenGermany
| | - Erik Lindahl
- Department of Applied PhysicsSwedish e‐Science Research Center, KTH Royal Institute of TechnologyStockholmSweden
- Department of Biochemistry and Biophysics, Science for Life LaboratoryStockholm UniversityStockholmSweden
| | | | - Adam Hospital
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and TechnologyBarcelonaSpain
| | - Modesto Orozco
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and TechnologyBarcelonaSpain
- Department of Biochemistry and BiomedicineUniversity of BarcelonaBarcelonaSpain
| |
Collapse
|
14
|
Bayarri G, Andrio P, Hospital A, Orozco M, Gelpí JL. BioExcel Building Blocks Workflows (BioBB-Wfs), an integrated web-based platform for biomolecular simulations. Nucleic Acids Res 2022; 50:W99-W107. [PMID: 35639735 PMCID: PMC9252775 DOI: 10.1093/nar/gkac380] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2022] [Revised: 04/19/2022] [Accepted: 05/02/2022] [Indexed: 11/15/2022] Open
Abstract
We present BioExcel Building Blocks Workflows, a web-based graphical user interface (GUI) offering access to a collection of transversal pre-configured biomolecular simulation workflows assembled with the BioExcel Building Blocks library. Available workflows include Molecular Dynamics setup, protein-ligand docking, trajectory analyses and small molecule parameterization. Workflows can be launched in the platform or downloaded to be run in the users’ own premises. Remote launching of long executions to user's available High-Performance computers is possible, only requiring configuration of the appropriate access credentials. The web-based graphical user interface offers a high level of interactivity, with integration with the NGL viewer to visualize and check 3D structures, MDsrv to visualize trajectories, and Plotly to explore 2D plots. The server requires no login but is recommended to store the users’ projects and manage sensitive information such as remote credentials. Private projects can be made public and shared with colleagues with a simple URL. The tool will help biomolecular simulation users with the most common and repetitive processes by means of a very intuitive and interactive graphical user interface. The server is accessible at https://mmb.irbbarcelona.org/biobb-wfs.
Collapse
Affiliation(s)
- Genís Bayarri
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and Technology. Baldiri Reixac 10-12, 08028 Barcelona, Spain
| | - Pau Andrio
- Barcelona Supercomputing Center (BSC), Jordi Girona 29, 08034, Barcelona, Spain
| | - Adam Hospital
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and Technology. Baldiri Reixac 10-12, 08028 Barcelona, Spain
| | - Modesto Orozco
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and Technology. Baldiri Reixac 10-12, 08028 Barcelona, Spain.,Department of Biochemistry and Molecular Biology, University of Barcelona, 08028 Barcelona, Spain
| | - Josep Lluís Gelpí
- Barcelona Supercomputing Center (BSC), Jordi Girona 29, 08034, Barcelona, Spain.,Department of Biochemistry and Molecular Biology, University of Barcelona, 08028 Barcelona, Spain
| |
Collapse
|
15
|
Esmaeeli R, Andal B, Perez A. Searching for Low Probability Opening Events in a DNA Sliding Clamp. Life (Basel) 2022; 12:life12020261. [PMID: 35207548 PMCID: PMC8876151 DOI: 10.3390/life12020261] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/30/2021] [Revised: 02/01/2022] [Accepted: 02/03/2022] [Indexed: 11/27/2022] Open
Abstract
The β subunit of E. coli DNA polymererase III is a DNA sliding clamp associated with increasing the processivity of DNA synthesis. In its free form, it is a circular homodimer structure that can accomodate double-stranded DNA in a nonspecific manner. An open state of the clamp must be accessible before loading the DNA. The opening mechanism is still a matter of debate, as is the effect of bound DNA on opening/closing kinetics. We use a combination of atomistic, coarse-grained, and enhanced sampling strategies in both explicit and implicit solvents to identify opening events in the sliding clamp. Such simulations of large nucleic acid and their complexes are becoming available and are being driven by improvements in force fields and the creation of faster computers. Different models support alternative opening mechanisms, either through an in-plane or out-of-plane opening event. We further note some of the current limitations, despite advances, in modeling these highly charged systems with implicit solvent.
Collapse
|
16
|
Whitford PC. Overview of the Biomolecular Association and Dynamics session at the 20th IUPAB congress, 45th Brazilian congress of SBBF, and the 50th annual meeting of SBBq. Biophys Rev 2022; 13:863-865. [PMID: 35059009 DOI: 10.1007/s12551-021-00884-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2021] [Accepted: 11/01/2021] [Indexed: 05/31/2023] Open
Abstract
In this session, experts in molecular biophysics described the dynamics of biopolymers across a wide range of length and time scales. This discussion highlighted numerous techniques that span from highly detailed simulations, to coarse-grained theoretical models, as well as high-resolution structural analysis. The topics were equally diverse, where there was discussion of biological processes at small (individual atoms), intermediate (assemblies) and very large scales (phase separation).
Collapse
Affiliation(s)
- Paul Charles Whitford
- Department of Physics and Center for Theoretical Biological Physics, Northeastern University, Boston, MA USA
| |
Collapse
|
17
|
Martin B, Dans PD, Wieczór M, Villegas N, Brun-Heath I, Battistini F, Terrazas M, Orozco M. Molecular basis of Arginine and Lysine DNA sequence-dependent thermo-stability modulation. PLoS Comput Biol 2022; 18:e1009749. [PMID: 35007284 PMCID: PMC8782489 DOI: 10.1371/journal.pcbi.1009749] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2021] [Revised: 01/21/2022] [Accepted: 12/13/2021] [Indexed: 11/19/2022] Open
Abstract
We have used a variety of theoretical and experimental techniques to study the role of four basic amino acids-Arginine, Lysine, Ornithine and L-2,4-Diaminobutyric acid-on the structure, flexibility and sequence-dependent stability of DNA. We found that the presence of organic ions stabilizes the duplexes and significantly reduces the difference in stability between AT- and GC-rich duplexes with respect to the control conditions. This suggests that these amino acids, ingredients of the primordial soup during abiogenesis, could have helped to equalize the stability of AT- and GC-rich DNA oligomers, facilitating a general non-catalysed self-replication of DNA. Experiments and simulations demonstrate that organic ions have an effect that goes beyond the general electrostatic screening, involving specific interactions along the grooves of the double helix. We conclude that organic ions, largely ignored in the DNA world, should be reconsidered as crucial structural elements far from mimics of small inorganic cations.
Collapse
Affiliation(s)
- Benjamin Martin
- Institute for Research in Biomedicine (IRB Barcelona), the Barcelona Institute of Science and Technology, Barcelona, Spain
- Institute of Bioengineering, School of Life Sciences, Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
| | - Pablo D. Dans
- Institute for Research in Biomedicine (IRB Barcelona), the Barcelona Institute of Science and Technology, Barcelona, Spain
- Department of Biological Sciences, CENUR Litoral Norte, Universidad de la República (UdelaR), Salto, Uruguay
- Functional Genomics Unit, Institut Pasteur de Montevideo, Montevideo, Uruguay
| | - Milosz Wieczór
- Institute for Research in Biomedicine (IRB Barcelona), the Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Nuria Villegas
- Institute for Research in Biomedicine (IRB Barcelona), the Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Isabelle Brun-Heath
- Institute for Research in Biomedicine (IRB Barcelona), the Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Federica Battistini
- Institute for Research in Biomedicine (IRB Barcelona), the Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Montserrat Terrazas
- Institute for Research in Biomedicine (IRB Barcelona), the Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Modesto Orozco
- Institute for Research in Biomedicine (IRB Barcelona), the Barcelona Institute of Science and Technology, Barcelona, Spain
- Department of Biochemistry and Biomedicine, Faculty of Biology, University of Barcelona, Barcelona, Spain
| |
Collapse
|
18
|
Homologous basic helix–loop–helix transcription factors induce distinct deformations of torsionally-stressed DNA: a potential transcription regulation mechanism. QRB DISCOVERY 2022. [PMID: 37529292 PMCID: PMC10392670 DOI: 10.1017/qrd.2022.5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
Abstract
Abstract
Changing torsional restraints on DNA is essential for the regulation of transcription. Torsional stress, introduced by RNA polymerase, can propagate along chromatin facilitating topological transitions and modulating the specific binding of transcription factors (TFs) to DNA. Despite the importance, the mechanistic details on how torsional stress impacts the TFs-DNA complexation remain scarce. Herein, we address the impact of torsional stress on DNA complexation with homologous human basic helix–loop–helix (BHLH) hetero- and homodimers: MycMax, MadMax and MaxMax. The three TF dimers exhibit specificity towards the same DNA consensus sequence, the E-box response element, while regulating different transcriptional pathways. Using microseconds-long atomistic molecular dynamics simulations together with the torsional restraint that controls DNA total helical twist, we gradually over- and underwind naked and complexed DNA to a maximum of ± 5°/bp step. We observe that the binding of the BHLH dimers results in a similar increase in DNA torsional rigidity. However, under torsional stress the BHLH dimers induce distinct DNA deformations, characterised by changes in DNA grooves geometry and a significant asymmetric DNA bending. Supported by bioinformatics analyses, our data suggest that torsional stress may contribute to the execution of differential transcriptional programs of the homologous TFs by modulating their collaborative interactions.
Collapse
|
19
|
Yoo J, Park S, Maffeo C, Ha T, Aksimentiev A. DNA sequence and methylation prescribe the inside-out conformational dynamics and bending energetics of DNA minicircles. Nucleic Acids Res 2021; 49:11459-11475. [PMID: 34718725 PMCID: PMC8599915 DOI: 10.1093/nar/gkab967] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2020] [Revised: 09/27/2021] [Accepted: 10/11/2021] [Indexed: 11/13/2022] Open
Abstract
Eukaryotic genome and methylome encode DNA fragments' propensity to form nucleosome particles. Although the mechanical properties of DNA possibly orchestrate such encoding, the definite link between 'omics' and DNA energetics has remained elusive. Here, we bridge the divide by examining the sequence-dependent energetics of highly bent DNA. Molecular dynamics simulations of 42 intact DNA minicircles reveal that each DNA minicircle undergoes inside-out conformational transitions with the most likely configuration uniquely prescribed by the nucleotide sequence and methylation of DNA. The minicircles' local geometry consists of straight segments connected by sharp bends compressing the DNA's inward-facing major groove. Such an uneven distribution of the bending stress favors minimum free energy configurations that avoid stiff base pair sequences at inward-facing major grooves. Analysis of the minicircles' inside-out free energy landscapes yields a discrete worm-like chain model of bent DNA energetics that accurately account for its nucleotide sequence and methylation. Experimentally measuring the dependence of the DNA looping time on the DNA sequence validates the model. When applied to a nucleosome-like DNA configuration, the model quantitatively reproduces yeast and human genomes' nucleosome occupancy. Further analyses of the genome-wide chromatin structure data suggest that DNA bending energetics is a fundamental determinant of genome architecture.
Collapse
Affiliation(s)
- Jejoong Yoo
- Department of Physics, Sungkyunkwan University, Suwon 16419, Republic of Korea
| | - Sangwoo Park
- Department of Biophysics and Biophysical Chemistry, Johns Hopkins University, Baltimore, MD 21205, USA
| | - Christopher Maffeo
- Department of Physics and the Center for the Physics of Living Cells, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Taekjip Ha
- Department of Biophysics and Biophysical Chemistry, Johns Hopkins University, Baltimore, MD 21205, USA
- Howard Hughes Medical Institute, Baltimore, MD 21218, USA
| | - Aleksei Aksimentiev
- Department of Physics and the Center for the Physics of Living Cells, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
- Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| |
Collapse
|
20
|
Sequence-dependent structural properties of B-DNA: what have we learned in 40 years? Biophys Rev 2021; 13:995-1005. [DOI: 10.1007/s12551-021-00893-8] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/14/2021] [Accepted: 11/01/2021] [Indexed: 11/27/2022] Open
|
21
|
Battistini F, Dans PD, Terrazas M, Castellazzi CL, Portella G, Labrador M, Villegas N, Brun-Heath I, González C, Orozco M. The Impact of the HydroxyMethylCytosine epigenetic signature on DNA structure and function. PLoS Comput Biol 2021; 17:e1009547. [PMID: 34748533 PMCID: PMC8601608 DOI: 10.1371/journal.pcbi.1009547] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2021] [Revised: 11/18/2021] [Accepted: 10/10/2021] [Indexed: 12/30/2022] Open
Abstract
We present a comprehensive, experimental and theoretical study of the impact of 5-hydroxymethylation of DNA cytosine. Using molecular dynamics, biophysical experiments and NMR spectroscopy, we found that Ten-Eleven translocation (TET) dioxygenases generate an epigenetic variant with structural and physical properties similar to those of 5-methylcytosine. Experiments and simulations demonstrate that 5-methylcytosine (mC) and 5-hydroxymethylcytosine (hmC) generally lead to stiffer DNA than normal cytosine, with poorer circularization efficiencies and lower ability to form nucleosomes. In particular, we can rule out the hypothesis that hydroxymethylation reverts to unmodified cytosine physical properties, as hmC is even more rigid than mC. Thus, we do not expect dramatic changes in the chromatin structure induced by differences in physical properties between d(mCpG) and d(hmCpG). Conversely, our simulations suggest that methylated-DNA binding domains (MBDs), associated with repression activities, are sensitive to the substitution d(mCpG) ➔ d(hmCpG), while MBD3 which has a dual activation/repression activity is not sensitive to the d(mCpG) d(hmCpG) change. Overall, while gene activity changes due to cytosine methylation are the result of the combination of stiffness-related chromatin reorganization and MBD binding, those associated to 5-hydroxylation of methylcytosine could be explained by a change in the balance of repression/activation pathways related to differential MBD binding. In Eukaryotic cells, DNA epigenetic modifications play an important role in gene expression and regulation, and protein recognition. In this work we investigate the physical implications of cytosine 5-hydroxymethylation on DNA, its structural and flexibility differences with methylated and unmodified cytosine using molecular dynamics, biophysical experiments and NMR spectroscopy. In particular the effect of hydroxyl group on free energy of nucleosome and Methyl binding Protein (MBD) binding, comparing in silico and experimental data to shed light on the effect of the reduced flexibility and the direct protein-DNA recognition.
Collapse
Affiliation(s)
- Federica Battistini
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Pablo D. Dans
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and Technology, Barcelona, Spain
- Department of Biological Sciences, CENUR Litoral Norte, Universidad de la República (UdelaR), Salto, Uruguay
- Functional Genomics Lab., Institut Pasteur of Montevideo, Montevideo, Uruguay
| | - Montserrat Terrazas
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Chiara L. Castellazzi
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Guillem Portella
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and Technology, Barcelona, Spain
- Chemistry Department, University of Cambridge, Cambridge, United Kingdom
| | - Mireia Labrador
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Núria Villegas
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Isabelle Brun-Heath
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Carlos González
- Instituto Química Física Rocasolano, Consejo Superior de Investigaciones Científicas (CSIC), Madrid, Spain
| | - Modesto Orozco
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and Technology, Barcelona, Spain
- Department of Biochemistry and Molecular Biology, University of Barcelona, Barcelona, Spain
- * E-mail:
| |
Collapse
|
22
|
Dohnalová H, Lankaš F. Deciphering the mechanical properties of
B‐DNA
duplex. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2021. [DOI: 10.1002/wcms.1575] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Hana Dohnalová
- Department of Informatics and Chemistry University of Chemistry and Technology Prague Praha 6 Czech Republic
| | - Filip Lankaš
- Department of Informatics and Chemistry University of Chemistry and Technology Prague Praha 6 Czech Republic
| |
Collapse
|
23
|
Hörberg J, Moreau K, Tamás MJ, Reymer A. Sequence-specific dynamics of DNA response elements and their flanking sites regulate the recognition by AP-1 transcription factors. Nucleic Acids Res 2021; 49:9280-9293. [PMID: 34387667 PMCID: PMC8450079 DOI: 10.1093/nar/gkab691] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2020] [Revised: 07/21/2021] [Accepted: 07/27/2021] [Indexed: 11/28/2022] Open
Abstract
Activator proteins 1 (AP-1) comprise one of the largest families of eukaryotic basic leucine zipper transcription factors. Despite advances in the characterization of AP-1 DNA-binding sites, our ability to predict new binding sites and explain how the proteins achieve different gene expression levels remains limited. Here we address the role of sequence-specific DNA flexibility for stability and specific binding of AP-1 factors, using microsecond-long molecular dynamics simulations. As a model system, we employ yeast AP-1 factor Yap1 binding to three different response elements from two genetic environments. Our data show that Yap1 actively exploits the sequence-specific flexibility of DNA within the response element to form stable protein–DNA complexes. The stability also depends on the four to six flanking nucleotides, adjacent to the response elements. The flanking sequences modulate the conformational adaptability of the response element, making it more shape-efficient to form specific contacts with the protein. Bioinformatics analysis of differential expression of the studied genes supports our conclusions: the stability of Yap1–DNA complexes, modulated by the flanking environment, influences the gene expression levels. Our results provide new insights into mechanisms of protein–DNA recognition and the biological regulation of gene expression levels in eukaryotes.
Collapse
Affiliation(s)
- Johanna Hörberg
- Department of Chemistry and Molecular Biology, University of Gothenburg, Gothenburg 40530, Sweden
| | - Kevin Moreau
- Department of Chemistry and Molecular Biology, University of Gothenburg, Gothenburg 40530, Sweden
| | - Markus J Tamás
- Department of Chemistry and Molecular Biology, University of Gothenburg, Gothenburg 40530, Sweden
| | - Anna Reymer
- Department of Chemistry and Molecular Biology, University of Gothenburg, Gothenburg 40530, Sweden
| |
Collapse
|
24
|
Beyond the double helix: DNA structural diversity and the PDB. J Biol Chem 2021; 296:100553. [PMID: 33744292 PMCID: PMC8063756 DOI: 10.1016/j.jbc.2021.100553] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2020] [Revised: 01/15/2021] [Accepted: 03/16/2021] [Indexed: 12/11/2022] Open
Abstract
The determination of the double helical structure of DNA in 1953 remains the landmark event in the development of modern biological and biomedical science. This structure has also been the starting point for the determination of some 2000 DNA crystal structures in the subsequent 68 years. Their structural diversity has extended to the demonstration of sequence-dependent local structure in duplex DNA, to DNA bending in short and long sequences and in the DNA wound round the nucleosome, and to left-handed duplex DNAs. Beyond the double helix itself, in circumstances where DNA sequences are or can be induced to unwind from being duplex, a wide variety of topologies and forms can exist. Quadruplex structures, based on four-stranded cores of stacked G-quartets, are prevalent though not randomly distributed in the human and other genomes and can play roles in transcription, translation, and replication. Yet more complex folds can result in DNAs with extended tertiary structures and enzymatic/catalytic activity. The Protein Data Bank is the depository of all these structures, and the resource where structures can be critically examined and validated, as well as compared one with another to facilitate analysis of conformational and base morphology features. This review will briefly survey the major structural classes of DNAs and illustrate their significance, together with some examples of how the use of the Protein Data Bank by for example, data mining, has illuminated DNA structural concepts.
Collapse
|
25
|
Mondal M, Yang L, Cai Z, Patra P, Gao YQ. A perspective on the molecular simulation of DNA from structural and functional aspects. Chem Sci 2021; 12:5390-5409. [PMID: 34168783 PMCID: PMC8179617 DOI: 10.1039/d0sc05329e] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Abstract
As genetic material, DNA not only carries genetic information by sequence, but also affects biological functions ranging from base modification to replication, transcription and gene regulation through its structural and dynamic properties and variations. The motion and structural properties of DNA involved in related biological processes are also multi-scale, ranging from single base flipping to local DNA deformation, TF binding, G-quadruplex and i-motif formation, TAD establishment, compartmentalization and even chromosome territory formation, just to name a few. The sequence-dependent physical properties of DNA play vital role in all these events, and thus it is interesting to examine how simple sequence information affects DNA and the formation of the chromatin structure in these different hierarchical orders. Accordingly, molecular simulations can provide atomistic details of interactions and conformational dynamics involved in different biological processes of DNA, including those inaccessible by current experimental methods. In this perspective, which is mainly based on our recent studies, we provide a brief overview of the atomistic simulations on how the hierarchical structure and dynamics of DNA can be influenced by its sequences, base modifications, environmental factors and protein binding in the context of the protein-DNA interactions, gene regulation and structural organization of chromatin. We try to connect the DNA sequence, the hierarchical structures of DNA and gene regulation.
Collapse
Affiliation(s)
- Manas Mondal
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory 518055 Shenzhen China
| | - Lijiang Yang
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University 100871 Beijing China
| | - Zhicheng Cai
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University 100871 Beijing China.,Biomedical Pioneering Innovation Center, Peking University 100871 Beijing China
| | - Piya Patra
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory 518055 Shenzhen China .,Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University 100871 Beijing China
| | - Yi Qin Gao
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory 518055 Shenzhen China .,Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University 100871 Beijing China.,Biomedical Pioneering Innovation Center, Peking University 100871 Beijing China.,Beijing Advanced Innovation Center for Genomics, Peking University 100871 Beijing China
| |
Collapse
|
26
|
Zerze GH, Stillinger FH, Debenedetti PG. Thermodynamics of DNA Hybridization from Atomistic Simulations. J Phys Chem B 2021; 125:771-779. [PMID: 33434025 DOI: 10.1021/acs.jpcb.0c09237] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Studying DNA hybridization equilibrium at atomistic length scales, either via molecular dynamics (MD) or through commonly used advanced sampling approaches, is notoriously difficult. In this work, we describe an order-parameter-based advanced sampling technique to calculate the free energy of hybridization, and estimate the melting temperature of DNA oligomers at atomistic resolution. The free energy landscapes are reported as a function of a native-topology-based order parameter for the Drew-Dickerson dodecamer and for a range of DNA decamer sequences of different GC content. Our estimated melting temperatures match the experimental numbers within ±15 °C. As a test of the numerical reliability of the procedures employed, it was verified that the predicted free energy surfaces and melting temperatures of the d- and l-enantiomers of the Drew-Dickerson dodecamer were indistinguishable within numerical accuracy.
Collapse
Affiliation(s)
- Gül H Zerze
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Frank H Stillinger
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
| | - Pablo G Debenedetti
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| |
Collapse
|
27
|
Ben Imeddourene A, Zargarian L, Buckle M, Hartmann B, Mauffret O. Slow motions in A·T rich DNA sequence. Sci Rep 2020; 10:19005. [PMID: 33149183 PMCID: PMC7642443 DOI: 10.1038/s41598-020-75645-x] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2020] [Accepted: 10/12/2020] [Indexed: 01/09/2023] Open
Abstract
In free B-DNA, slow (microsecond-to-millisecond) motions that involve equilibrium between Watson-Crick (WC) and Hoogsteen (HG) base-pairing expand the DNA dynamic repertoire that could mediate DNA-protein assemblies. R1ρ relaxation dispersion NMR methods are powerful tools to capture such slow conformational exchanges in solution using 13C/15 N labelled DNA. Here, these approaches were applied to a dodecamer containing a TTAAA element that was assumed to facilitate nucleosome formation. NMR data and inferred exchange parameters assign HG base pairs as the minor, transient conformers specifically observed in three successive A·T base pairs forming the TAA·TTA segment. The abundance of these HG A·T base pairs can be up to 1.2% which is high compared to what has previously been observed. Data analyses support a scenario in which the three adenines undergo non-simultaneous motions despite their spatial proximity, thus optimising the probability of having one HG base pair in the TAA·TTA segment. Finally, revisiting previous NMR data on H2 resonance linewidths on the basis of our results promotes the idea of there being a special propensity of A·T base pairs in TAA·TTA tracts to adopt HG pairing. In summary, this study provides an example of a DNA functional element submitted to slow conformational exchange. More generally, it strengthens the importance of the role of the DNA sequence in modulating its dynamics, over a nano- to milli-second time scale.
Collapse
Affiliation(s)
- A Ben Imeddourene
- LBPA, ENS de Paris-Saclay, UMR 8113 CNRS, Institut D'Alembert, Université Paris-Saclay, 4, avenue des Sciences, 91190, Gif-sur-Yvette, France
| | - L Zargarian
- LBPA, ENS de Paris-Saclay, UMR 8113 CNRS, Institut D'Alembert, Université Paris-Saclay, 4, avenue des Sciences, 91190, Gif-sur-Yvette, France
| | - M Buckle
- LBPA, ENS de Paris-Saclay, UMR 8113 CNRS, Institut D'Alembert, Université Paris-Saclay, 4, avenue des Sciences, 91190, Gif-sur-Yvette, France
| | - B Hartmann
- LBPA, ENS de Paris-Saclay, UMR 8113 CNRS, Institut D'Alembert, Université Paris-Saclay, 4, avenue des Sciences, 91190, Gif-sur-Yvette, France
| | - O Mauffret
- LBPA, ENS de Paris-Saclay, UMR 8113 CNRS, Institut D'Alembert, Université Paris-Saclay, 4, avenue des Sciences, 91190, Gif-sur-Yvette, France.
| |
Collapse
|
28
|
Retureau R, Foloppe N, Elbahnsi A, Oguey C, Hartmann B. A dynamic view of DNA structure within the nucleosome: Biological implications. J Struct Biol 2020; 211:107511. [PMID: 32311461 DOI: 10.1016/j.jsb.2020.107511] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2020] [Revised: 04/10/2020] [Accepted: 04/13/2020] [Indexed: 01/21/2023]
Abstract
Most of eukaryotic cellular DNA is packed in nucleosome core particles (NCPs), in which the DNA (DNANCP) is wrapped around histones. The influence of this organization on the intrinsic local dynamics of DNA is largely unknown, in particular because capturing such information from experiments remains notoriously challenging. Given the importance of dynamical properties in DNA functions, we addressed this issue using CHARMM36 MD simulations of a nucleosome containing the NCP positioning 601 sequence and four related free dodecamers. Comparison between DNANCP and free DNA reveals a limited impact of the dense DNA-histone interface on correlated motions of dinucleotide constituents and on fluctuations of inter base pair parameters. A characteristic feature intimately associated with the DNANCP super-helical path is a set of structural periodicities that includes a marked alternation of regions enriched in backbone BI and BII conformers. This observation led to uncover a convincing correspondence between the sequence effect on BI/BII propensities in both DNANCP and free DNA, strengthening the idea that the histone preference for particular DNA sequences relies on those intrinsic structural properties. These results offer for the first time a detailed view of the DNA dynamical behavior within NCP. They show in particular that the DNANCP dynamics is substantial enough to preserve the ability to structurally adjust to external proteins, for instance remodelers. Also, fresh structural arguments highlight the relevance of relationships between DNA sequence and structural properties for NCP formation. Overall, our work offers a more rational framework to approach the functional, biological roles of NCP.
Collapse
Affiliation(s)
- Romain Retureau
- Université Paris-Saclay, ENS Paris-Saclay, CNRS, Laboratoire de biologie et pharmacologie appliquée, 61 avenue du Président Wilson, 94235 Cachan cedex, France
| | | | - Ahmad Elbahnsi
- Université Paris-Saclay, ENS Paris-Saclay, CNRS, Laboratoire de biologie et pharmacologie appliquée, 61 avenue du Président Wilson, 94235 Cachan cedex, France; LPTM, UMR8089, CNRS, CY Cergy Paris Université, 2 avenue Adolphe Chauvin, 95302 Cergy-Pontoise, France
| | - Christophe Oguey
- LPTM, UMR8089, CNRS, CY Cergy Paris Université, 2 avenue Adolphe Chauvin, 95302 Cergy-Pontoise, France
| | - Brigitte Hartmann
- Université Paris-Saclay, ENS Paris-Saclay, CNRS, Laboratoire de biologie et pharmacologie appliquée, 61 avenue du Président Wilson, 94235 Cachan cedex, France.
| |
Collapse
|
29
|
Walther J, Dans PD, Balaceanu A, Hospital A, Bayarri G, Orozco M. A multi-modal coarse grained model of DNA flexibility mappable to the atomistic level. Nucleic Acids Res 2020; 48:e29. [PMID: 31956910 PMCID: PMC7049737 DOI: 10.1093/nar/gkaa015] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2019] [Revised: 12/10/2019] [Accepted: 01/04/2020] [Indexed: 12/13/2022] Open
Abstract
We present a new coarse grained method for the simulation of duplex DNA. The algorithm uses a generalized multi-harmonic model that can represent any multi-normal distribution of helical parameters, thus avoiding caveats of current mesoscopic models for DNA simulation and representing a breakthrough in the field. The method has been parameterized from accurate parmbsc1 atomistic molecular dynamics simulations of all unique tetranucleotide sequences of DNA embedded in long duplexes and takes advantage of the correlation between helical states and backbone configurations to derive atomistic representations of DNA. The algorithm, which is implemented in a simple web interface and in a standalone package reproduces with high computational efficiency the structural landscape of long segments of DNA untreatable by atomistic molecular dynamics simulations.
Collapse
Affiliation(s)
- Jürgen Walther
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, 08028 Barcelona, Spain
| | - Pablo D Dans
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, 08028 Barcelona, Spain
- Department of Biological Sciences, CENUR North Coast, University of the Republic, 50000 Salto, Uruguay
| | - Alexandra Balaceanu
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, 08028 Barcelona, Spain
| | - Adam Hospital
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, 08028 Barcelona, Spain
| | - Genís Bayarri
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, 08028 Barcelona, Spain
| | - Modesto Orozco
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, 08028 Barcelona, Spain
- Department of Biochemistry and Biomedicine, The University of Barcelona, 08028 Barcelona, Spain
| |
Collapse
|
30
|
Velasco-Berrelleza V, Burman M, Shepherd JW, Leake MC, Golestanian R, Noy A. SerraNA: a program to determine nucleic acids elasticity from simulation data. Phys Chem Chem Phys 2020; 22:19254-19266. [DOI: 10.1039/d0cp02713h] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Abstract
AT-rich motifs can generate extreme mechanical properties, which are critical for creating strong global bends when phased properly.
Collapse
Affiliation(s)
| | | | | | - Mark C. Leake
- Department of Physics
- University of York
- York
- UK
- Department of Biology
| | - Ramin Golestanian
- Max Planck Institute for Dynamics and Self-Organization (MPIDS)
- Göttingen
- Germany
- Rudolf Peierls Center for Theoretical Physics
- University of Oxford
| | - Agnes Noy
- Department of Physics
- University of York
- York
- UK
| |
Collapse
|
31
|
Hospital A, Battistini F, Soliva R, Gelpí JL, Orozco M. Surviving the deluge of biosimulation data. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2019. [DOI: 10.1002/wcms.1449] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- Adam Hospital
- Institut de Recerca Biomèdica, IRB Barcelona, The Barcelona Institute of Science and Technology Joint IRB‐BSC Program in Computational Biology Barcelona Spain
| | - Federica Battistini
- Institut de Recerca Biomèdica, IRB Barcelona, The Barcelona Institute of Science and Technology Joint IRB‐BSC Program in Computational Biology Barcelona Spain
| | | | - Josep Lluis Gelpí
- Barcelona Supercomputing Center Join IRB‐BSC Program in Computational Biology Barcelona Spain
- Departament de Bioquímica i Biomedicina Universitat de Barcelona Barcelona Spain
| | - Modesto Orozco
- Institut de Recerca Biomèdica, IRB Barcelona, The Barcelona Institute of Science and Technology Joint IRB‐BSC Program in Computational Biology Barcelona Spain
- Departament de Bioquímica i Biomedicina Universitat de Barcelona Barcelona Spain
| |
Collapse
|