1
|
Junier I, Ghobadpour E, Espeli O, Everaers R. DNA supercoiling in bacteria: state of play and challenges from a viewpoint of physics based modeling. Front Microbiol 2023; 14:1192831. [PMID: 37965550 PMCID: PMC10642903 DOI: 10.3389/fmicb.2023.1192831] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2023] [Accepted: 09/25/2023] [Indexed: 11/16/2023] Open
Abstract
DNA supercoiling is central to many fundamental processes of living organisms. Its average level along the chromosome and over time reflects the dynamic equilibrium of opposite activities of topoisomerases, which are required to relax mechanical stresses that are inevitably produced during DNA replication and gene transcription. Supercoiling affects all scales of the spatio-temporal organization of bacterial DNA, from the base pair to the large scale chromosome conformation. Highlighted in vitro and in vivo in the 1960s and 1970s, respectively, the first physical models were proposed concomitantly in order to predict the deformation properties of the double helix. About fifteen years later, polymer physics models demonstrated on larger scales the plectonemic nature and the tree-like organization of supercoiled DNA. Since then, many works have tried to establish a better understanding of the multiple structuring and physiological properties of bacterial DNA in thermodynamic equilibrium and far from equilibrium. The purpose of this essay is to address upcoming challenges by thoroughly exploring the relevance, predictive capacity, and limitations of current physical models, with a specific focus on structural properties beyond the scale of the double helix. We discuss more particularly the problem of DNA conformations, the interplay between DNA supercoiling with gene transcription and DNA replication, its role on nucleoid formation and, finally, the problem of scaling up models. Our primary objective is to foster increased collaboration between physicists and biologists. To achieve this, we have reduced the respective jargon to a minimum and we provide some explanatory background material for the two communities.
Collapse
Affiliation(s)
- Ivan Junier
- CNRS, UMR 5525, VetAgro Sup, Grenoble INP, TIMC, Université Grenoble Alpes, Grenoble, France
| | - Elham Ghobadpour
- CNRS, UMR 5525, VetAgro Sup, Grenoble INP, TIMC, Université Grenoble Alpes, Grenoble, France
- École Normale Supérieure (ENS) de Lyon, CNRS, Laboratoire de Physique and Centre Blaise Pascal de l'ENS de Lyon, Lyon, France
| | - Olivier Espeli
- Center for Interdisciplinary Research in Biology (CIRB), Collège de France, CNRS, INSERM, Université PSL, Paris, France
| | - Ralf Everaers
- École Normale Supérieure (ENS) de Lyon, CNRS, Laboratoire de Physique and Centre Blaise Pascal de l'ENS de Lyon, Lyon, France
| |
Collapse
|
2
|
Gutiérrez Fosado YA, Landuzzi F, Sakaue T. Coarse Graining DNA: Symmetry, Nonlocal Elasticity, and Persistence Length. PHYSICAL REVIEW LETTERS 2023; 130:058402. [PMID: 36800451 DOI: 10.1103/physrevlett.130.058402] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/03/2022] [Accepted: 01/06/2023] [Indexed: 06/18/2023]
Abstract
While the behavior of double-stranded DNA at mesoscopic scales is fairly well understood, less is known about its relation to the rich mechanical properties in the base-pair scale, which is crucial, for instance, to understand DNA-protein interactions and the nucleosome diffusion mechanism. Here, by employing the rigid base-pair model, we connect its microscopic parameters to the persistence length. Combined with all-atom molecular dynamic simulations, our scheme identifies relevant couplings between different degrees of freedom at each coarse-graining step. This allows us to clarify how the scale dependence of the elastic moduli is determined in a systematic way encompassing the role of previously unnoticed off-site couplings between deformations with different parity.
Collapse
Affiliation(s)
- Yair Augusto Gutiérrez Fosado
- School of Physics and Astronomy, University of Edinburgh, Peter Guthrie Tait Road, Edinburgh, EH9 3FD, United Kingdom
| | - Fabio Landuzzi
- Centro CMP3VdA, Istituto Italiano di Tecnologia, via Lavoratori Vittime del Col du Mont 28, 11100, Aosta, Italy
| | - Takahiro Sakaue
- Department of Physics and Mathematics, Aoyama Gakuin University, 5-10-1 Fuchinobe, Chuo-ku, Sagamihara-shi, Kanagawa 252-5258, Japan
| |
Collapse
|
3
|
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
|
4
|
Li Z, Sun R, Bishop TC. Genome Dashboards: Framework and Examples. Biophys J 2020; 118:2077-2085. [PMID: 32171420 PMCID: PMC7203004 DOI: 10.1016/j.bpj.2020.02.018] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2019] [Revised: 01/28/2020] [Accepted: 02/18/2020] [Indexed: 12/11/2022] Open
Abstract
Genomics is a sequence-based informatics science and a three-dimensional-structure-based material science. However, in practice, most genomics researchers utilize sequence-based informatics approaches or three-dimensional-structure-based material science techniques, not both. This division is, at least in part, the result of historical developments rather than a fundamental necessity. The underlying computational tools, experimental techniques, and theoretical models were developed independently. The primary result presented here is a framework for the unification of informatics- and physics-based data associated with DNA, nucleosomes, and chromatin. The framework is based on the mathematical representation of geometrically exact rods and the generalization of DNA basepair step parameters. Data unification enables researchers to integrate computational, experimental, and theoretical approaches for the study of chromatin biology. The framework can be implemented using model-view-controller design principles, existing genome browsers, and existing molecular visualization tools. We developed a minimal, web-based genome dashboard, G-Dash-min, and applied it to two simple examples to demonstrate the usefulness of data unification and proof of concept. Genome dashboards developed using the framework and design principles presented here are extensible and customizable and are therefore more broadly applicable than the examples presented. We expect a number of purpose-specific genome dashboards to emerge as a novel means of investigating structure-function relationships for genomes that range from basepairs to entire chromosomes and for generating, validating, and testing mechanistic hypotheses.
Collapse
Affiliation(s)
- Zilong Li
- Chemistry and Physics, Louisiana Tech University, Ruston, Louisiana
| | - Ran Sun
- Chemistry and Physics, Louisiana Tech University, Ruston, Louisiana
| | - Thomas C Bishop
- Chemistry and Physics, Louisiana Tech University, Ruston, Louisiana.
| |
Collapse
|
5
|
Dohnalová H, Dršata T, Šponer J, Zacharias M, Lipfert J, Lankaš F. Compensatory Mechanisms in Temperature Dependence of DNA Double Helical Structure: Bending and Elongation. J Chem Theory Comput 2020; 16:2857-2863. [DOI: 10.1021/acs.jctc.0c00037] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Hana Dohnalová
- Department of Informatics and Chemistry, University of Chemistry and Technology Prague, Technická 5, 166 28 Prague, Czech Republic
| | - Tomáš Dršata
- Department of Informatics and Chemistry, University of Chemistry and Technology Prague, Technická 5, 166 28 Prague, Czech Republic
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Martin Zacharias
- Physics-Department T38, Technical University of Munich, James-Franck-Strasse 1, 85748 Garching, Germany
| | - Jan Lipfert
- Department of Physics and Center for Nanoscience, LMU Munich, Amalienstrasse 54, 80799 Munich, Germany
| | - Filip Lankaš
- Department of Informatics and Chemistry, University of Chemistry and Technology Prague, Technická 5, 166 28 Prague, Czech Republic
| |
Collapse
|
6
|
Chakraborty M, Xu C, White AD. Encoding and selecting coarse-grain mapping operators with hierarchical graphs. J Chem Phys 2018; 149:134106. [PMID: 30292213 DOI: 10.1063/1.5040114] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
Coarse-grained (CG) molecular dynamics (MD) can simulate systems inaccessible to fine-grained (FG) MD simulations. A CG simulation decreases the degrees of freedom by mapping atoms from an FG representation into agglomerate CG particles. The FG to CG mapping is not unique. Research into systematic selection of these mappings is challenging due to their combinatorial growth with respect to the number of atoms in a molecule. Here we present a method of reducing the total count of mappings by imposing molecular topology and symmetry constraints. The count reduction is illustrated by considering all mappings for nearly 50 000 molecules. The resulting number of mapping operators is still large, so we introduce a novel hierarchical graphical approach which encodes multiple CG mapping operators. The encoding method is demonstrated for methanol and a 14-mer peptide. With the test cases, we show how the encoding can be used for automated selection of reasonable CG mapping operators.
Collapse
Affiliation(s)
- Maghesree Chakraborty
- Department of Chemical Engineering, University of Rochester, Rochester, New York 14627, USA
| | - Chenliang Xu
- Department of Computer Science, University of Rochester, Rochester, New York 14627, USA
| | - Andrew D White
- Department of Chemical Engineering, University of Rochester, Rochester, New York 14627, USA
| |
Collapse
|
7
|
Silanteva I, Komolkin AV. Representation of DNA environment: Spiral staircase distribution function. J Comput Chem 2018; 39:2300-2306. [PMID: 30299550 DOI: 10.1002/jcc.25549] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2018] [Revised: 07/09/2018] [Accepted: 07/11/2018] [Indexed: 12/17/2022]
Abstract
In the present study, we investigated the local structure of DNA and its environment using a new visualization technique. The spiral staircase distribution function (SSDF) is determined as two-dimensional density distribution of atoms of water and ligands in local reference frames linked with each base pair of poly-DNA molecule, either GC or AT. This property of SSDF provides opportunity to study sequence-specific binding of ions, peptides, and other agents derived from a molecular dynamics computer simulation. The spatial structure of double-stranded DNA environment in water solution containing either Mg2+ or Na+ ions was investigated using of SSDF. The distributions of ions around GC and AT base pairs are shown separately. It is observed that Mg2+ ions interact with DNA atoms by means of the layer of water molecules and penetrate into the major groove only. Na+ ions have a direct contact with DNA atoms and penetrate both into the major and minor grooves of DNA. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Irina Silanteva
- Faculty of Physics, Saint Petersburg State University, Universitetsaya emb., 7-9, Saint Petersburg, 199034, Russian Federation
| | - Andrei V Komolkin
- Faculty of Physics, Saint Petersburg State University, Universitetsaya emb., 7-9, Saint Petersburg, 199034, Russian Federation
| |
Collapse
|
8
|
Kriegel F, Matek C, Dršata T, Kulenkampff K, Tschirpke S, Zacharias M, Lankaš F, Lipfert J. The temperature dependence of the helical twist of DNA. Nucleic Acids Res 2018; 46:7998-8009. [PMID: 30053087 PMCID: PMC6125625 DOI: 10.1093/nar/gky599] [Citation(s) in RCA: 38] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2017] [Revised: 06/08/2018] [Accepted: 07/20/2018] [Indexed: 01/11/2023] Open
Abstract
DNA is the carrier of all cellular genetic information and increasingly used in nanotechnology. Quantitative understanding and optimization of its functions requires precise experimental characterization and accurate modeling of DNA properties. A defining feature of DNA is its helicity. DNA unwinds with increasing temperature, even for temperatures well below the melting temperature. However, accurate quantitation of DNA unwinding under external forces and a microscopic understanding of the corresponding structural changes are currently lacking. Here we combine single-molecule magnetic tweezers measurements with atomistic molecular dynamics and coarse-grained simulations to obtain a comprehensive view of the temperature dependence of DNA twist. Experimentally, we find that DNA twist changes by ΔTw(T) = (-11.0 ± 1.2)°/(°C·kbp), independent of applied force, in the range of forces where torque-induced melting is negligible. Our atomistic simulations predict ΔTw(T) = (-11.1 ± 0.3)°/(°C·kbp), in quantitative agreement with experiments, and suggest that the untwisting of DNA with temperature is predominantly due to changes in DNA structure for defined backbone substates, while the effects of changes in substate populations are minor. Coarse-grained simulations using the oxDNA framework yield a value of ΔTw(T) = (-6.4 ± 0.2)°/(°C·kbp) in semi-quantitative agreement with experiments.
Collapse
Affiliation(s)
- Franziska Kriegel
- Department of Physics, Nanosystems Initiative Munich, and Center for Nanoscience, LMU Munich, Amalienstr. 54, 80799 Munich, Germany
| | - Christian Matek
- Technical University of Munich and Institute of Computational Biology, Helmholtz Zentrum München-German Research Center for Environmental Health, Ingolstädter Landstr. 1, 85764 Neuherberg, Germany
| | - Tomáš Dršata
- Department of Informatics and Chemistry, University of Chemistry and Technology Prague, Technická 5, 166 28 Prague, Czech Republic
| | - Klara Kulenkampff
- Department of Physics, Nanosystems Initiative Munich, and Center for Nanoscience, LMU Munich, Amalienstr. 54, 80799 Munich, Germany
| | - Sophie Tschirpke
- Department of Physics, Nanosystems Initiative Munich, and Center for Nanoscience, LMU Munich, Amalienstr. 54, 80799 Munich, Germany
| | - Martin Zacharias
- Physics-Department T38, Technical University of Munich, James-Franck-Strasse 1, 85748 Garching, Germany
| | - Filip Lankaš
- Department of Informatics and Chemistry, University of Chemistry and Technology Prague, Technická 5, 166 28 Prague, Czech Republic
| | - Jan Lipfert
- Department of Physics, Nanosystems Initiative Munich, and Center for Nanoscience, LMU Munich, Amalienstr. 54, 80799 Munich, Germany
| |
Collapse
|
9
|
Langenberg M, Jackson NE, de Pablo JJ, Müller M. Role of translational entropy in spatially inhomogeneous, coarse-grained models. J Chem Phys 2018. [DOI: 10.1063/1.5018178] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Marcel Langenberg
- Institut für Theoretische Physik, Universität Göttingen, 37077 Göttingen, Germany
| | - Nicholas E. Jackson
- Institute for Molecular Engineering, University of Chicago, Chicago, Illinois 60615, USA
- Argonne National Laboratory, Lemont, Illinois 06349, USA
| | - Juan J. de Pablo
- Institute for Molecular Engineering, University of Chicago, Chicago, Illinois 60615, USA
- Argonne National Laboratory, Lemont, Illinois 06349, USA
| | - Marcus Müller
- Institut für Theoretische Physik, Universität Göttingen, 37077 Göttingen, Germany
| |
Collapse
|
10
|
Pan K, Bricker WP, Ratanalert S, Bathe M. Structure and conformational dynamics of scaffolded DNA origami nanoparticles. Nucleic Acids Res 2017; 45:6284-6298. [PMID: 28482032 PMCID: PMC5499760 DOI: 10.1093/nar/gkx378] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2016] [Accepted: 04/25/2017] [Indexed: 12/22/2022] Open
Abstract
Synthetic DNA is a highly programmable nanoscale material that can be designed to self-assemble into 3D structures that are fully determined by underlying Watson–Crick base pairing. The double crossover (DX) design motif has demonstrated versatility in synthesizing arbitrary DNA nanoparticles on the 5–100 nm scale for diverse applications in biotechnology. Prior computational investigations of these assemblies include all-atom and coarse-grained modeling, but modeling their conformational dynamics remains challenging due to their long relaxation times and associated computational cost. We apply all-atom molecular dynamics and coarse-grained finite element modeling to DX-based nanoparticles to elucidate their fine-scale and global conformational structure and dynamics. We use our coarse-grained model with a set of secondary structural motifs to predict the equilibrium solution structures of 45 DX-based DNA origami nanoparticles including a tetrahedron, octahedron, icosahedron, cuboctahedron and reinforced cube. Coarse-grained models are compared with 3D cryo-electron microscopy density maps for these five DNA nanoparticles and with all-atom molecular dynamics simulations for the tetrahedron and octahedron. Our results elucidate non-intuitive atomic-level structural details of DX-based DNA nanoparticles, and offer a general framework for efficient computational prediction of global and local structural and mechanical properties of DX-based assemblies that are inaccessible to all-atom based models alone.
Collapse
Affiliation(s)
- Keyao Pan
- Department of Biological Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | - William P Bricker
- Department of Biological Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | - Sakul Ratanalert
- Department of Biological Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA.,Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | - Mark Bathe
- Department of Biological Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| |
Collapse
|
11
|
Mitchell JS, Glowacki J, Grandchamp AE, Manning RS, Maddocks JH. Sequence-Dependent Persistence Lengths of DNA. J Chem Theory Comput 2017; 13:1539-1555. [DOI: 10.1021/acs.jctc.6b00904] [Citation(s) in RCA: 64] [Impact Index Per Article: 9.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
| | - Jaroslaw Glowacki
- Ecole Polytechnique Fédérale de Lausanne, Lausanne CH 1273, Switzerland
| | | | | | - John H. Maddocks
- Ecole Polytechnique Fédérale de Lausanne, Lausanne CH 1273, Switzerland
| |
Collapse
|
12
|
Abstract
More than 20 coarse-grained (CG) DNA models have been developed for simulating the behavior of this molecule under various conditions, including those required for nanotechnology. However, none of these models reproduces the DNA polymorphism associated with conformational changes in the ribose rings of the DNA backbone. These changes make an essential contribution to the DNA local deformability and provide the possibility of the transition of the DNA double helix from the B-form to the A-form during interactions with biological molecules. We propose a CG representation of the ribose conformational flexibility. We substantiate the choice of the CG sites (six per nucleotide) needed for the "sugar" GC DNA model, and obtain the potentials of the CG interactions between the sites by the "bottom-up" approach using the all-atom AMBER force field. We show that the representation of the ribose flexibility requires one non-harmonic and one three-particle potential, the forms of both the potentials being different from the ones generally used. The model also includes (i) explicit representation of ions (in an implicit solvent) and (ii) sequence dependence. With these features, the sugar CG DNA model reproduces (with the same parameters) both the B- and A- stable forms under corresponding conditions and demonstrates both the A to B and the B to A phase transitions. Graphical Abstract The proposed coarse-grained DNA model allows to reproduce both the B- and A- DNA forms and the transitions between them under corresponding conditions.
Collapse
|
13
|
Salari H, Eslami-Mossallam B, Ranjbar HF, Ejtehadi MR. Stiffer double-stranded DNA in two-dimensional confinement due to bending anisotropy. Phys Rev E 2017; 94:062407. [PMID: 28085439 DOI: 10.1103/physreve.94.062407] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2016] [Indexed: 11/07/2022]
Abstract
Using analytical approach and Monte Carlo (MC) simulations, we study the elastic behavior of the intrinsically twisted elastic ribbons with bending anisotropy, such as double-stranded DNA (dsDNA), in two-dimensional (2D) confinement. We show that, due to the bending anisotropy, the persistence length of dsDNA in 2D conformations is always greater than three-dimensional (3D) conformations. This result is in consistence with the measured values for DNA persistence length in 2D and 3D in equal biological conditions. We also show that in two dimensions, an anisotropic, intrinsically twisted polymer exhibits an implicit twist-bend coupling, which leads to the transient curvature increasing with a half helical turn periodicity along the bent polymer.
Collapse
Affiliation(s)
- H Salari
- Department of Physics, Sharif University of Technology, P.O. Box 11155-9161, Tehran, Iran
| | - B Eslami-Mossallam
- Department of Bionanoscience, Kavli Institute of Nanoscience, Delft University of Technology, Lorentzweg 1, 2628 CJ Delft, the Netherlands
| | - H F Ranjbar
- Institute of Complex Systems (ICS-2), Forschungszentrum Jülich, Wilhelm-Johnen-Straße, 52425 Jülich, Germany
| | - M R Ejtehadi
- Department of Physics, Sharif University of Technology, P.O. Box 11155-9161, Tehran, Iran and School of Nano Science, Institute for Research in Fundamental Sciences (IPM), Tehran 19395-5531, Iran
| |
Collapse
|
14
|
Mondal M, Halder S, Chakrabarti J, Bhattacharyya D. Hybrid simulation approach incorporating microscopic interaction along with rigid body degrees of freedom for stacking between base pairs. Biopolymers 2015; 105:212-26. [PMID: 26600167 DOI: 10.1002/bip.22787] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2015] [Revised: 10/19/2015] [Accepted: 11/17/2015] [Indexed: 11/07/2022]
Abstract
Stacking interaction between the aromatic heterocyclic bases plays an important role in the double helical structures of nucleic acids. Considering the base as rigid body, there are total of 18 degrees of freedom of a dinucleotide step. Some of these parameters show sequence preferences, indicating that the detailed atomic interactions are important in the stacking. Large variants of non-canonical base pairs have been seen in the crystallographic structures of RNA. However, their stacking preferences are not thoroughly deciphered yet from experimental results. The current theoretical approaches use either the rigid body degrees of freedom where the atomic information are lost or computationally expensive all atom simulations. We have used a hybrid simulation approach incorporating Monte-Carlo Metropolis sampling in the hyperspace of 18 stacking parameters where the interaction energies using AMBER-parm99bsc0 and CHARMM-36 force-fields were calculated from atomic positions. We have also performed stacking energy calculations for structures from Monte-Carlo ensemble by Dispersion corrected density functional theory. The available experimental data with Watson-Crick base pairs are compared to establish the validity of the method. Stacking interaction involving A:U and G:C base pairs with non-canonical G:U base pairs also were calculated and showed that these structures were also sequence dependent. This approach could be useful to generate multiscale modeling of nucleic acids in terms of coarse-grained parameters where the atomic interactions are preserved. This method would also be useful to predict structure and dynamics of different base pair steps containing non Watson-Crick base pairs, as found often in the non-coding RNA structures. © 2015 Wiley Periodicals, Inc. Biopolymers 105: 212-226, 2016.
Collapse
Affiliation(s)
- Manas Mondal
- Computational Science Division, Saha Institute of Nuclear Physics, 1/AF Bidhannagar, Kolkata, 700 064, India
| | - Sukanya Halder
- Computational Science Division, Saha Institute of Nuclear Physics, 1/AF Bidhannagar, Kolkata, 700 064, India
| | - Jaydeb Chakrabarti
- Department of Chemical, Biological and Macro-Molecular Sciences, S.N. Bose National Center for Basic Sciences, Sector III, Salt Lake, Kolkata, 700 098, India
| | - Dhananjay Bhattacharyya
- Computational Science Division, Saha Institute of Nuclear Physics, 1/AF Bidhannagar, Kolkata, 700 064, India
| |
Collapse
|
15
|
Liebl K, Drsata T, Lankas F, Lipfert J, Zacharias M. Explaining the striking difference in twist-stretch coupling between DNA and RNA: A comparative molecular dynamics analysis. Nucleic Acids Res 2015; 43:10143-56. [PMID: 26464435 PMCID: PMC4666353 DOI: 10.1093/nar/gkv1028] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2015] [Accepted: 09/28/2015] [Indexed: 01/07/2023] Open
Abstract
Double stranded helical DNA and RNA are flexible molecules that can undergo global conformational fluctuations. Their bending, twisting and stretching deformabilities are of similar magnitude. However, recent single-molecule experiments revealed a striking qualitative difference indicating an opposite sign for the twist-stretch couplings of dsDNA and dsRNA [Lipfert et al. 2014. Proc. Natl. Acad. Sci. U.S.A. 111, 15408] that is not explained by existing models. Employing unconstrained Molecular Dynamics (MD) simulations we are able to reproduce the qualitatively different twist-stretch coupling for dsDNA and dsRNA in semi-quantitative agreement with experiment. Similar results are also found in simulations that include an external torque to induce over- or unwinding of DNA and RNA. Detailed analysis of the helical deformations coupled to twist indicate that the interplay of helical rise, base pair inclination and displacement from the helix axis upon twist changes are responsible for the different twist-stretch correlations. Overwinding of RNA results in more compact conformations with a narrower major groove and consequently reduced helical extension. Overwinding of DNA decreases the size of the minor groove and the resulting positive base pair inclination leads to a slender and more extended helical structure.
Collapse
Affiliation(s)
- Korbinian Liebl
- Physik-Department T38, Technische Universität München, James-Franck-Strasse, D-85748 Garching, Germany
| | - Tomas Drsata
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Flemingovo namesti 2, 166 10 Prague, Czech Republic Department of Physical and Macromolecular Chemistry, Faculty of Science, Charles University Prague, Albertov 6, 128 43 Prague, Czech Republic
| | - Filip Lankas
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Flemingovo namesti 2, 166 10 Prague, Czech Republic Laboratory of Informatics and Chemistry, University of Chemistry and Technology Prague, Technická 5, 166 28 Prague, Czech Republic
| | - Jan Lipfert
- Department of Physics, Center for Nanoscience (CeNS), and Nanosystems Initiative Munich (NIM), Ludwig-Maximilian-University Munich, 80799 Munich, Germany
| | - Martin Zacharias
- Physik-Department T38, Technische Universität München, James-Franck-Strasse, D-85748 Garching, Germany
| |
Collapse
|
16
|
Endeward B, Marko A, Denysenkov VP, Sigurdsson ST, Prisner TF. Advanced EPR Methods for Studying Conformational Dynamics of Nucleic Acids. Methods Enzymol 2015; 564:403-25. [PMID: 26477259 DOI: 10.1016/bs.mie.2015.06.007] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
Pulsed electron paramagnetic resonance (EPR) spectroscopy has become an important tool for structural characterization of biomolecules allowing measurement of the distances between two paramagnetic spin labels attached to a biomolecule in the 2-8 nm range. In this chapter, we will focus on applications of this approach to investigate tertiary structure elements as well as conformational dynamics of nucleic acid molecules. Both aspects take advantage of using specific spin labels that are rigidly attached to the nucleobases, as they allow obtaining not only the distance but also the relative orientation between both nitroxide moieties with high accuracy. Thus, not only the distance but additionally the three Euler angles between both the nitroxide axis systems and the two polar angles of the interconnecting vector with respect to the nitroxide axis systems can be extracted from a single pair of spin labels. To extract all these parameters independently and unambiguously, a set of multifrequency/multifield pulsed EPR experiments have to be performed. We will describe the experimental procedure as well as newly developed spin labels, which are helpful to disentangle all these parameters, and tools which we have developed to analyze such data sets. The procedures and analyses will be illustrated by examples from our laboratory.
Collapse
Affiliation(s)
- B Endeward
- Institute of Physical and Theoretical Chemistry and Center of Biomolecular Magnetic Resonance, Goethe University Frankfurt am Main, Frankfurt am Main, Germany
| | - A Marko
- Institute of Physical and Theoretical Chemistry and Center of Biomolecular Magnetic Resonance, Goethe University Frankfurt am Main, Frankfurt am Main, Germany
| | - V P Denysenkov
- Institute of Physical and Theoretical Chemistry and Center of Biomolecular Magnetic Resonance, Goethe University Frankfurt am Main, Frankfurt am Main, Germany
| | - S Th Sigurdsson
- Department of Chemistry, Science Institute, University of Iceland, Reykjavık, Iceland
| | - T F Prisner
- Institute of Physical and Theoretical Chemistry and Center of Biomolecular Magnetic Resonance, Goethe University Frankfurt am Main, Frankfurt am Main, Germany.
| |
Collapse
|
17
|
Stachiewicz A, Molski A. A coarse-grained MARTINI-like force field for DNA unzipping in nanopores. J Comput Chem 2015; 36:947-56. [DOI: 10.1002/jcc.23874] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2014] [Revised: 01/14/2015] [Accepted: 01/17/2015] [Indexed: 11/08/2022]
Affiliation(s)
- Anna Stachiewicz
- Department of Chemistry; Adam Mickiewicz University; Poznan Poland
| | - Andrzej Molski
- Department of Chemistry; Adam Mickiewicz University; Poznan Poland
| |
Collapse
|
18
|
Fathizadeh A, Schiessel H, Ejtehadi MR. Molecular Dynamics Simulation of Supercoiled DNA Rings. Macromolecules 2014. [DOI: 10.1021/ma501660w] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Affiliation(s)
- Arman Fathizadeh
- School
of Physics, Institute for Research in Fundamental Sciences (IPM), Tehran, Iran
- Institute
for Nanoscience and Nanotechnology, Sharif University of Technology, Tehran, Iran
| | - Helmut Schiessel
- Instituut-Lorentz
for Theoretical Physics, P.O. Box 9506, 2300 RA Leiden, The Netherlands
| | - Mohammad Reza Ejtehadi
- Department
of Physics, Sharif University of Technology, P.O. Box 11155-8639, Tehran, Iran
| |
Collapse
|
19
|
|
20
|
Shen L, Hu H. Resolution-Adapted All-Atomic and Coarse-Grained Model for Biomolecular Simulations. J Chem Theory Comput 2014; 10:2528-36. [DOI: 10.1021/ct401029k] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- Lin Shen
- Department
of Chemistry, University of Hong Kong, Pokfulam Road, Hong Kong, China
| | - Hao Hu
- Department
of Chemistry, University of Hong Kong, Pokfulam Road, Hong Kong, China
- The University of Hong Kong-Shenzhen Institute of Research and Innovation, Kejizhong Second Road, Shenzhen, China
| |
Collapse
|
21
|
Meyer S, Jost D, Theodorakopoulos N, Peyrard M, Lavery R, Everaers R. Temperature dependence of the DNA double helix at the nanoscale: structure, elasticity, and fluctuations. Biophys J 2014; 105:1904-14. [PMID: 24138866 DOI: 10.1016/j.bpj.2013.09.004] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2013] [Revised: 09/01/2013] [Accepted: 09/09/2013] [Indexed: 10/26/2022] Open
Abstract
Biological organisms exist over a broad temperature range of -15°C to +120°C, where many molecular processes involving DNA depend on the nanoscale properties of the double helix. Here, we present results of extensive molecular dynamics simulations of DNA oligomers at different temperatures. We show that internal basepair conformations are strongly temperature-dependent, particularly in the stretch and opening degrees of freedom whose harmonic fluctuations can be considered the initial steps of the DNA melting pathway. The basepair step elasticity contains a weaker, but detectable, entropic contribution in the roll, tilt, and rise degrees of freedom. To extend the validity of our results to the temperature interval beyond the standard melting transition relevant to extremophiles, we estimate the effects of superhelical stress on the stability of the basepair steps, as computed from the Benham model. We predict that although the average twist decreases with temperature in vitro, the stabilizing external torque in vivo results in an increase of ∼1°/bp (or a superhelical density of Δσ ≃ +0.03) in the interval 0-100°C. In the final step, we show that the experimentally observed apparent bending persistence length of torsionally unconstrained DNA can be calculated from a hybrid model that accounts for the softening of the double helix and the presence of transient denaturation bubbles. Although the latter dominate the behavior close to the melting transition, the inclusion of helix softening is important around standard physiological temperatures.
Collapse
Affiliation(s)
- Sam Meyer
- Université de Lyon, Laboratoire de Physique, Ecole Normale Supérieure de Lyon, Lyon, France; Centre Blaise Pascal, Ecole Normale Supérieure de Lyon, Lyon, France
| | | | | | | | | | | |
Collapse
|
22
|
|
23
|
Doye JPK, Ouldridge TE, Louis AA, Romano F, Šulc P, Matek C, Snodin BEK, Rovigatti L, Schreck JS, Harrison RM, Smith WPJ. Coarse-graining DNA for simulations of DNA nanotechnology. Phys Chem Chem Phys 2013; 15:20395-414. [PMID: 24121860 DOI: 10.1039/c3cp53545b] [Citation(s) in RCA: 139] [Impact Index Per Article: 12.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
To simulate long time and length scale processes involving DNA it is necessary to use a coarse-grained description. Here we provide an overview of different approaches to such coarse-graining, focussing on those at the nucleotide level that allow the self-assembly processes associated with DNA nanotechnology to be studied. OxDNA, our recently-developed coarse-grained DNA model, is particularly suited to this task, and has opened up this field to systematic study by simulations. We illustrate some of the range of DNA nanotechnology systems to which the model is being applied, as well as the insights it can provide into fundamental biophysical properties of DNA.
Collapse
Affiliation(s)
- Jonathan P K Doye
- Physical & Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford, OX1 3QZ, UK
| | | | | | | | | | | | | | | | | | | | | |
Collapse
|
24
|
Dršata T, Kara M, Zacharias M, Lankaš F. Effect of 8-oxoguanine on DNA structure and deformability. J Phys Chem B 2013; 117:11617-22. [PMID: 24028561 DOI: 10.1021/jp407562t] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
8-Oxoguanine (oxoG) is an abundant product of oxidative DNA damage. It is removed by repair glycosylases, but exactly how the enzymes recognize oxoG in the large surplus of undamaged bases is not fully understood. The lesion may induce changes in the properties of naked DNA that facilitate the recognition. In this work, we assess the effect of oxoG on DNA structure and mechanical deformability. We performed extensive unrestrained, atomic resolution molecular dynamics simulations to parametrize a nonlocal, rigid base mechanical model of DNA. Our data indicate that oxoG induces unwinding of the base pair step at the 5'-side of the lesion. This brings the damaged DNA closer to its conformation in the initial complex with bacterial glycosylase MutM. The untwisting is partially caused by different BII substate populations and is further enhanced by the base-sugar repulsion within oxoG. On the other hand, our analysis shows that damaged and undamaged DNA have very similar harmonic stiffness. These results suggest an indirect readout component of the MutM-DNA initial complex formation. They also help one to understand the effect of oxoG on the formation of nucleosomes and looped gene regulatory complexes.
Collapse
Affiliation(s)
- Tomáš Dršata
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic , Flemingovo náměstı́ 2, 166 10, Praha 6, Czech Republic
| | | | | | | |
Collapse
|
25
|
Fathizadeh A, Heidari M, Eslami-Mossallam B, Ejtehadi MR. Confinement dynamics of a semiflexible chain inside nano-spheres. J Chem Phys 2013; 139:044912. [PMID: 23902024 DOI: 10.1063/1.4816626] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
We study the conformations of a semiflexible chain, confined in nano-scaled spherical cavities, under two distinct processes of confinement. Radial contraction and packaging are employed as two confining procedures. The former method is performed by gradually decreasing the diameter of a spherical shell which envelopes a confined chain. The latter procedure is carried out by injecting the chain inside a spherical shell through a hole on the shell surface. The chain is modeled with a rigid body molecular dynamics simulation and its parameters are adjusted to DNA base-pair elasticity. Directional order parameter is employed to analyze and compare the confined chain and the conformations of the chain for two different sizes of the spheres are studied in both procedures. It is shown that for the confined chains in the sphere sizes of our study, they appear in spiral or tennis-ball structures, and the tennis-ball structure is more likely to be observed in more compact confinements. Our results also show that the dynamical procedure of confinement and the rate of the confinement are influential parameters of the structure of the chain inside spherical cavities.
Collapse
Affiliation(s)
- A Fathizadeh
- Sharif University of Technology, Institute for Nanoscience and Nanotechnology, P.O. Box 14588-89694, Tehran, Iran
| | | | | | | |
Collapse
|
26
|
Dršata T, Lankaš F. Theoretical models of DNA flexibility. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2013. [DOI: 10.1002/wcms.1144] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
|
27
|
Gonzalez O, Petkevičiūtė D, Maddocks JH. A sequence-dependent rigid-base model of DNA. J Chem Phys 2013; 138:055102. [DOI: 10.1063/1.4789411] [Citation(s) in RCA: 44] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
|
28
|
Hsu CW, Fyta M, Lakatos G, Melchionna S, Kaxiras E. Ab initio determination of coarse-grained interactions in double-stranded DNA. J Chem Phys 2013; 137:105102. [PMID: 22979896 DOI: 10.1063/1.4748105] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
We derive the coarse-grained interactions between DNA nucleotides from ab initio total-energy calculations based on density functional theory (DFT). The interactions take into account base and sequence specificity, and are decomposed into physically distinct contributions that include hydrogen bonding, stacking interactions, backbone, and backbone-base interactions. The interaction energies of each contribution are calculated from DFT for a wide range of configurations and are fitted by simple analytical expressions for use in the coarse-grained model, which reduces each nucleotide into two sites. This model is not derived from experimental data, yet it successfully reproduces the stable B-DNA structure and gives good predictions for the persistence length. It may be used to realistically probe dynamics of DNA strands in various environments at the μs time scale and the μm length scale.
Collapse
Affiliation(s)
- Chia Wei Hsu
- Department of Physics, Harvard University, Cambridge, Massachusetts 02138, USA
| | | | | | | | | |
Collapse
|
29
|
Chung M, Koo B, Boxer SG. Formation and analysis of topographical domains between lipid membranes tethered by DNA hybrids of different lengths. Faraday Discuss 2013; 161:333-45; discussion 419-59. [PMID: 23805748 PMCID: PMC3703934 DOI: 10.1039/c2fd20108a] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We recently described a strategy to prepare DNA-tethered lipid membranes either to fixed DNA on a surface or to DNA displayed on a supported bilayer [Boxer et al., J. Struct. Biol., 2009, 168, 190; Boxer et al., Langmuir, 2011, 27, 5492]. With the latter system, the DNA hybrids are laterally mobile; when orthogonal sense-antisense pairs of different lengths are used, the DNA hybrids segregate by height and the tethered membrane deforms to accommodate the height difference. This architecture is particularly useful for modelling interactions between membranes mediated by molecular recognition and resembles cell-to-cell junctions. The length, affinity and population of the DNA hybrids between the membranes are completely controllable. Interesting patterns of height segregation are observed by fluorescence interference contrast microscopy. Diverse behavior is observed in the segregation and pattern forming process and possible mechanisms are discussed. This model system captures some of the essential physics of synapse formation and is a step towards understanding lipid membrane behaviour in cell-to-cell junctions.
Collapse
Affiliation(s)
- Minsub Chung
- Department of Chemistry, Stanford University, Stanford, CA 94305-5012 (USA); fax +1 650-723-4817; Tel. +1 650 723 4482;
| | - Bonjun Koo
- Department of Chemistry, Stanford University, Stanford, CA 94305-5012 (USA); fax +1 650-723-4817; Tel. +1 650 723 4482;
| | - Steven G. Boxer
- Department of Chemistry, Stanford University, Stanford, CA 94305-5012 (USA); fax +1 650-723-4817; Tel. +1 650 723 4482;
| |
Collapse
|
30
|
Šulc P, Romano F, Ouldridge TE, Rovigatti L, Doye JPK, Louis AA. Sequence-dependent thermodynamics of a coarse-grained DNA model. J Chem Phys 2012; 137:135101. [DOI: 10.1063/1.4754132] [Citation(s) in RCA: 215] [Impact Index Per Article: 17.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
31
|
Dans PD, Pérez A, Faustino I, Lavery R, Orozco M. Exploring polymorphisms in B-DNA helical conformations. Nucleic Acids Res 2012; 40:10668-78. [PMID: 23012264 PMCID: PMC3510489 DOI: 10.1093/nar/gks884] [Citation(s) in RCA: 73] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
The traditional mesoscopic paradigm represents DNA as a series of base-pair steps whose energy response to equilibrium perturbations is elastic, with harmonic oscillations (defining local stiffness) around a single equilibrium conformation. In addition, base sequence effects are often analysed as a succession of independent XpY base-pair steps (i.e. a nearest-neighbour (NN) model with only 10 unique cases). Unfortunately, recent massive simulations carried out by the ABC consortium suggest that the real picture of DNA flexibility may be much more complex. The paradigm of DNA flexibility therefore needs to be revisited. In this article, we explore in detail one of the most obvious violations of the elastic NN model of flexibility: the bimodal distributions of some helical parameters. We perform here an in-depth statistical analysis of a very large set of MD trajectories and also of experimental structures, which lead to very solid evidence of bimodality. We then suggest ways to improve mesoscopic models to account for this deviation from the elastic regime.
Collapse
Affiliation(s)
- Pablo D Dans
- Joint IRB-BSC Program on Computational Biology, Institute for Research in Biomedicine, Parc Cientific de Barcelona, Josep Samitier 1-5, Barcelona 08028, Spain
| | | | | | | | | |
Collapse
|
32
|
Wolfe KC, Hastings WA, Dutta S, Long A, Shapiro BA, Woolf TB, Guthold M, Chirikjian GS. Multiscale modeling of double-helical DNA and RNA: a unification through Lie groups. J Phys Chem B 2012; 116:8556-72. [PMID: 22676719 PMCID: PMC4833121 DOI: 10.1021/jp2126015] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
Several different mechanical models of double-helical nucleic-acid structures that have been presented in the literature are reviewed here together with a new analysis method that provides a reconciliation between these disparate models. In all cases, terminology and basic results from the theory of Lie groups are used to describe rigid-body motions in a coordinate-free way, and when necessary, coordinates are introduced in a way in which simple equations result. We consider double-helical DNAs and RNAs which, in their unstressed referential state, have backbones that are either straight, slightly precurved, or bent by the action of a protein or other bound molecule. At the coarsest level, we consider worm-like chains with anisotropic bending stiffness. Then, we show how bi-rod models converge to this for sufficiently long filament lengths. At a finer level, we examine elastic networks of rigid bases and show how these relate to the coarser models. Finally, we show how results from molecular dynamics simulation at full atomic resolution (which is the finest scale considered here) and AFM experimental measurements (which is at the coarsest scale) relate to these models.
Collapse
Affiliation(s)
- Kevin C. Wolfe
- Department of Mechanical Engineering, Johns Hopkins University, Baltimore, Maryland, United States
| | | | - Samrat Dutta
- Department of Physics, Wake Forest University, Winston-Salem, North Carolina, United States
| | - Andrew Long
- Department of Mechanical Engineering, Northwestern University, Evanston, Illinois, United States
| | - Bruce A. Shapiro
- Center for Cancer Research Nanobiology Program, Frederick National Laboratory for Cancer Research, National Cancer Institute, Frederick, Maryland, United States
| | - Thomas B. Woolf
- Department of Physiology, Johns Hopkins University School of Medicine, Baltimore, Maryland, United States
| | - Martin Guthold
- Department of Physics, Wake Forest University, Winston-Salem, North Carolina, United States
| | - Gregory S. Chirikjian
- Department of Mechanical Engineering, Johns Hopkins University, Baltimore, Maryland, United States
| |
Collapse
|
33
|
Réblová K, Šponer J, Lankaš F. Structure and mechanical properties of the ribosomal L1 stalk three-way junction. Nucleic Acids Res 2012; 40:6290-303. [PMID: 22451682 PMCID: PMC3401443 DOI: 10.1093/nar/gks258] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2012] [Revised: 03/07/2012] [Accepted: 03/07/2012] [Indexed: 01/06/2023] Open
Abstract
The L1 stalk is a key mobile element of the large ribosomal subunit which interacts with tRNA during translocation. Here, we investigate the structure and mechanical properties of the rRNA H76/H75/H79 three-way junction at the base of the L1 stalk from four different prokaryotic organisms. We propose a coarse-grained elastic model and parameterize it using large-scale atomistic molecular dynamics simulations. Global properties of the junction are well described by a model in which the H76 helix is represented by a straight, isotropically flexible elastic rod, while the junction core is represented by an isotropically flexible spherical hinge. Both the core and the helix contribute substantially to the overall H76 bending fluctuations. The presence of wobble pairs in H76 does not induce any increased flexibility or anisotropy to the helix. The half-closed conformation of the L1 stalk seems to be accessible by thermal fluctuations of the junction itself, without any long-range allosteric effects. Bending fluctuations of H76 with a bulge introduced in it suggest a rationale for the precise position of the bulge in eukaryotes. Our elastic model can be generalized to other RNA junctions found in biological systems or in nanotechnology.
Collapse
Affiliation(s)
- Kamila Réblová
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, CEITEC—Central European Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno and Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Flemingovo nám. 2, 166 10 Praha 6, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, CEITEC—Central European Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno and Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Flemingovo nám. 2, 166 10 Praha 6, Czech Republic
| | - Filip Lankaš
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, CEITEC—Central European Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno and Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Flemingovo nám. 2, 166 10 Praha 6, Czech Republic
| |
Collapse
|
34
|
Lankaš F. Modelling Nucleic Acid Structure and Flexibility: From Atomic to Mesoscopic Scale. INNOVATIONS IN BIOMOLECULAR MODELING AND SIMULATIONS 2012. [DOI: 10.1039/9781849735056-00001] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Abstract
This chapter surveys some of the recent developments in coarse-grained modelling of nucleic acids. We first discuss models based on pseudoatoms, effective spherical particles representing groups of atoms. A major part of the chapter is devoted to models in which bases or base pairs are represented as independent, interacting rigid bodies. Two popular definitions of internal coordinates, as implemented in the programs 3DNA and Curves+, are outlined from a common perspective. Recently developed rigid base and basepair models with nonlocal quadratic interactions are presented. A statistical mechanical description of the models on their full phase space yields exact relations between model parameters and expected values of some state functions. We estimated shape and stiffness parameters for nonlocal rigid base and basepair models of a DNA oligomer containing A-tract. The parameterization is based on atomic-resolution molecular dynamics simulation data. We found that the rigid base model is consistent with a local interaction pattern, while interactions in the rigid basepair model are visibly non-local, in agreement with earlier findings. Differences in shape and stiffness parameters obtained using Curves+ and 3DNA coordinates are found to be small for structures within the B-DNA family. Anharmonic effects, coarser models, and other approaches to describe nucleic acid structure and flexibility are discussed.
Collapse
Affiliation(s)
- Filip Lankaš
- Centre for Complex Molecular Systems and Biomolecules Institute of Organic Chemistry and Biochemistry Academy of Sciences of the Czech Republic, Flemingovo nám. 2, 166 10 Praha 6 Czech Republic
| |
Collapse
|
35
|
Meyer S, Becker NB, Syed SH, Goutte-Gattat D, Shukla MS, Hayes JJ, Angelov D, Bednar J, Dimitrov S, Everaers R. From crystal and NMR structures, footprints and cryo-electron-micrographs to large and soft structures: nanoscale modeling of the nucleosomal stem. Nucleic Acids Res 2011; 39:9139-54. [PMID: 21835779 PMCID: PMC3241633 DOI: 10.1093/nar/gkr573] [Citation(s) in RCA: 37] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2010] [Revised: 06/17/2011] [Accepted: 06/25/2011] [Indexed: 01/18/2023] Open
Abstract
The interaction of histone H1 with linker DNA results in the formation of the nucleosomal stem structure, with considerable influence on chromatin organization. In a recent paper [Syed,S.H., Goutte-Gattat,D., Becker,N., Meyer,S., Shukla,M.S., Hayes,J.J., Everaers,R., Angelov,D., Bednar,J. and Dimitrov,S. (2010) Single-base resolution mapping of H1-nucleosome interactions and 3D organization of the nucleosome. Proc. Natl Acad. Sci. USA, 107, 9620-9625], we published results of biochemical footprinting and cryo-electron-micrographs of reconstituted mono-, di- and tri-nucleosomes, for H1 variants with different lengths of the cationic C-terminus. Here, we present a detailed account of the analysis of the experimental data and we include thermal fluctuations into our nano-scale model of the stem structure. By combining (i) crystal and NMR structures of the nucleosome core particle and H1, (ii) the known nano-scale structure and elasticity of DNA, (iii) footprinting information on the location of protected sites on the DNA backbone and (iv) cryo-electron micrographs of reconstituted tri-nucleosomes, we arrive at a description of a polymorphic, hierarchically organized stem with a typical length of 20 ± 2 base pairs. A comparison to linker conformations inferred for poly-601 fibers with different linker lengths suggests, that intra-stem interactions stabilize and facilitate the formation of dense chromatin fibers.
Collapse
Affiliation(s)
- Sam Meyer
- Laboratoire de Physique and Centre Blaise Pascal, École Normale Supérieure de Lyon, Université de Lyon, CNRS UMR 5672, 46 allée d'Italie 69364 Lyon cedex 07, France.
| | | | | | | | | | | | | | | | | | | |
Collapse
|
36
|
Marko A, Denysenkov V, Margraf D, Cekan P, Schiemann O, Sigurdsson ST, Prisner TF. Conformational Flexibility of DNA. J Am Chem Soc 2011; 133:13375-9. [DOI: 10.1021/ja201244u] [Citation(s) in RCA: 96] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Affiliation(s)
- Andriy Marko
- Institute of Physical and Theoretical Chemistry and Center of Biomolecular Magnetic Resonance, Goethe University Frankfurt, Max-von-Laue-Str. 7, 60438 Frankfurt, Germany
| | - Vasyl Denysenkov
- Institute of Physical and Theoretical Chemistry and Center of Biomolecular Magnetic Resonance, Goethe University Frankfurt, Max-von-Laue-Str. 7, 60438 Frankfurt, Germany
| | - Dominik Margraf
- Institute of Physical and Theoretical Chemistry and Center of Biomolecular Magnetic Resonance, Goethe University Frankfurt, Max-von-Laue-Str. 7, 60438 Frankfurt, Germany
| | - Pavol Cekan
- Science Institute, University of Iceland, Dunhaga 3, 107 Reykjavik, Iceland
| | - Olav Schiemann
- Centre for Biomolecular Sciences, Centre of Magnetic Resonance, University of St Andrews, North Haugh, KY16 9ST St Andrews, U.K
| | | | - Thomas F. Prisner
- Institute of Physical and Theoretical Chemistry and Center of Biomolecular Magnetic Resonance, Goethe University Frankfurt, Max-von-Laue-Str. 7, 60438 Frankfurt, Germany
| |
Collapse
|
37
|
Savin AV, Mazo MA, Kikot IP, Manevitch LI, Onufriev AV. Heat conductivity of DNA double helix. PHYSICAL REVIEW. B, CONDENSED MATTER AND MATERIALS PHYSICS 2011; 83:245406. [PMID: 26207085 PMCID: PMC4508875 DOI: 10.1103/physrevb.83.245406] [Citation(s) in RCA: 38] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
Thermal conductivity of isolated single molecule DNA fragments is of importance for nanotechnology, but has not yet been measured experimentally. Theoretical estimates based on simplified (1D) models predict anomalously high thermal conductivity. To investigate thermal properties of single molecule DNA we have developed a 3D coarse-grained (CG) model that retains the realism of the full all-atom description, but is significantly more efficient. Within the proposed model each nucleotide is represented by 6 particles or grains; the grains interact via effective potentials inferred from classical molecular dynamics (MD) trajectories based on a well-established all-atom potential function. Comparisons of 10 ns long MD trajectories between the CG and the corresponding all-atom model show similar root-mean-square deviations from the canonical B-form DNA, and similar structural fluctuations. At the same time, the CG model is 10 to 100 times faster depending on the length of the DNA fragment in the simulation. Analysis of dispersion curves derived from the CG model yields longitudinal sound velocity and torsional stiffness in close agreement with existing experiments. The computational efficiency of the CG model makes it possible to calculate thermal conductivity of a single DNA molecule not yet available experimentally. For a uniform (polyG-polyC) DNA, the estimated conductivity coefficient is 0.3 W/mK which is half the value of thermal conductivity for water. This result is in stark contrast with estimates of thermal conductivity for simplified, effectively 1D chains ("beads on a spring") that predict anomalous (infinite) thermal conductivity. Thus, full 3D character of DNA double-helix retained in the proposed model appears to be essential for describing its thermal properties at a single molecule level.
Collapse
Affiliation(s)
- Alexander V Savin
- Semenov Institute of Chemical Physics, Russian Academy of Sciences, Moscow 119991, Russia
| | - Mikhail A Mazo
- Semenov Institute of Chemical Physics, Russian Academy of Sciences, Moscow 119991, Russia
| | - Irina P Kikot
- Semenov Institute of Chemical Physics, Russian Academy of Sciences, Moscow 119991, Russia
| | - Leonid I Manevitch
- Semenov Institute of Chemical Physics, Russian Academy of Sciences, Moscow 119991, Russia
| | - Alexey V Onufriev
- Departments of Computer Science and Physics, 2160C Torgersen Hall, Virginia Tech, Blacksburg, VA 24061, USA
| |
Collapse
|
38
|
Affiliation(s)
- Juan J. de Pablo
- Department of Chemical and Biological Engineering, University of Wisconsin–Madison, Madison, Wisconsin 53706;
| |
Collapse
|
39
|
Dans PD, Zeida A, Machado MR, Pantano S. A Coarse Grained Model for Atomic-Detailed DNA Simulations with Explicit Electrostatics. J Chem Theory Comput 2010; 6:1711-25. [PMID: 26615701 DOI: 10.1021/ct900653p] [Citation(s) in RCA: 111] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Coarse-grain (CG) techniques allow considerable extension of the accessible size and time scales in simulations of biological systems. Although many CG representations are available for the most common biomacromolecules, very few have been reported for nucleic acids. Here, we present a CG model for molecular dynamics simulations of DNA on the multi-microsecond time scale. Our model maps the complexity of each nucleotide onto six effective superatoms keeping the "chemical sense" of specific Watson-Crick recognition. Molecular interactions are evaluated using a classical Hamiltonian with explicit electrostatics calculated under the framework of the generalized Born approach. This CG representation is able to accurately reproduce experimental structures, breathing dynamics, and conformational transitions from the A to the B form in double helical fragments. The model achieves a good qualitative reproduction of temperature-driven melting and its dependence on size, ionic strength, and sequence specificity. Reconstruction of atomistic models from CG trajectories give remarkable agreement with structural, dynamic, and energetic features obtained from fully atomistic simulation, opening the possibility to acquire nearly atomic detail data from CG trajectories.
Collapse
Affiliation(s)
- Pablo D Dans
- Institut Pasteur de Montevideo, Mataojo 2020, CP 11400 Montevideo, Uruguay
| | - Ari Zeida
- Institut Pasteur de Montevideo, Mataojo 2020, CP 11400 Montevideo, Uruguay
| | - Matías R Machado
- Institut Pasteur de Montevideo, Mataojo 2020, CP 11400 Montevideo, Uruguay
| | - Sergio Pantano
- Institut Pasteur de Montevideo, Mataojo 2020, CP 11400 Montevideo, Uruguay
| |
Collapse
|
40
|
Abstract
Mathew-Fenn et al. (Reports, 17 October 2008, p. 446) reported unexpected distance fluctuations in short end-labeled DNA constructs and interpreted them as evidence for cooperative DNA stretching modes. We show that when accounting for a subtle linker leverage effect, their data can be understood within standard noncooperative DNA elasticity.
Collapse
Affiliation(s)
- Nils B Becker
- Centre Blaise Pascal et Laboratoire de Physique, CNRS UMR 5672, Ecole Normale Supérieure, Université de Lyon, 46 Allée d'Italie, 69007 Lyon, France
| | | |
Collapse
|
41
|
Friedrich BM, Jülicher F. Steering chiral swimmers along noisy helical paths. PHYSICAL REVIEW LETTERS 2009; 103:068102. [PMID: 19792615 DOI: 10.1103/physrevlett.103.068102] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/03/2009] [Indexed: 05/28/2023]
Abstract
Chemotaxis along helical paths towards a target releasing a chemoattractant is found in sperm cells and many microorganisms. We discuss the stochastic differential geometry of the noisy helical swimming path of a chiral swimmer. A chiral swimmer equipped with a simple feedback system can navigate in a concentration gradient of chemoattractant. We derive an effective equation for the alignment of helical paths with a concentration gradient which is related to the alignment of a dipole in an external field and discuss the chemotaxis index.
Collapse
|
42
|
Eslami-Mossallam B, Ejtehadi MR. Asymmetric elastic rod model for DNA. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2009; 80:011919. [PMID: 19658741 DOI: 10.1103/physreve.80.011919] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/26/2008] [Revised: 06/02/2009] [Indexed: 05/28/2023]
Abstract
In this paper we consider the anharmonic corrections to the anisotropic elastic rod model for DNA. Our model accounts for the difference between the bending energies of positive and negative rolls, which comes from the asymmetric structure of the DNA molecule. We will show that the model can explain the high flexibility of DNA at small length scales, as well as kink formation at high deformation limit.
Collapse
Affiliation(s)
- B Eslami-Mossallam
- Department of Physics, Sharif University of Technology, P.O. Box 11365-8639, 14588-89694 Tehran, Iran
| | | |
Collapse
|
43
|
Becker NB, Everaers R. DNA nanomechanics in the nucleosome. Structure 2009; 17:579-89. [PMID: 19368891 DOI: 10.1016/j.str.2009.01.013] [Citation(s) in RCA: 28] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2008] [Revised: 12/16/2008] [Accepted: 01/28/2009] [Indexed: 10/20/2022]
Abstract
The sequence-dependent mechanics of DNA-histone binding are essential for nucleosome positioning and mobility. Here we reanalyze nucleosome crystal structures in terms of the well-characterized base-pair scale DNA elasticity, extracting the forces and torques acting on all bound DNA base pairs. We find that the strongest forces follow a characteristic repeating pattern that recovers the 12 known DNA backbone-histone contact sites. DNA twist defects and histone point mutations modify this pattern in interpretable ways. Additional, irregular forces between contact sites reveal histone tail-DNA interactions, whereas requiring the absence of external forces leads to a structural refinement of linker DNA. Based on these observations, we propose a simple, structure-based mechanical model of the nucleosome that is able to explain the placement of DNA twist defects in 146 base-pair nucleosomes and allows an estimate of the elastic energy spectrum of nucleosome twist defect states.
Collapse
Affiliation(s)
- Nils B Becker
- Laboratoire de Physique de l'Ecole Normale Supérieure, Université de Lyon, France.
| | | |
Collapse
|
44
|
Becker NB, Everaers R. DNA nanomechanics: How proteins deform the double helix. J Chem Phys 2009; 130:135102. [DOI: 10.1063/1.3082157] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
45
|
Mazur AK. Kinetic and Thermodynamic DNA Elasticity at Micro- and Mesoscopic Scales. J Phys Chem B 2009; 113:2077-89. [DOI: 10.1021/jp8098945] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Affiliation(s)
- Alexey K. Mazur
- CNRS UPR9080, Institut de Biologie Physico-Chimique 13, rue Pierre et Marie Curie, Paris 75005, France
| |
Collapse
|
46
|
Cascella M, Neri MA, Carloni P, Dal Peraro M. Topologically Based Multipolar Reconstruction of Electrostatic Interactions in Multiscale Simulations of Proteins. J Chem Theory Comput 2008; 4:1378-85. [DOI: 10.1021/ct800122x] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Michele Cascella
- Laboratory of Computational Chemistry and Biochemistry and Laboratory for Biomolecular Modeling, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland, International School for Advanced Studies (SISSA/ISAS) and CNR-INFM-DEMOCRITOS, I-34014 Trieste, Italy, and Italian Institute of Technology, Italy
| | - Marilisa A. Neri
- Laboratory of Computational Chemistry and Biochemistry and Laboratory for Biomolecular Modeling, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland, International School for Advanced Studies (SISSA/ISAS) and CNR-INFM-DEMOCRITOS, I-34014 Trieste, Italy, and Italian Institute of Technology, Italy
| | - Paolo Carloni
- Laboratory of Computational Chemistry and Biochemistry and Laboratory for Biomolecular Modeling, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland, International School for Advanced Studies (SISSA/ISAS) and CNR-INFM-DEMOCRITOS, I-34014 Trieste, Italy, and Italian Institute of Technology, Italy
| | - Matteo Dal Peraro
- Laboratory of Computational Chemistry and Biochemistry and Laboratory for Biomolecular Modeling, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland, International School for Advanced Studies (SISSA/ISAS) and CNR-INFM-DEMOCRITOS, I-34014 Trieste, Italy, and Italian Institute of Technology, Italy
| |
Collapse
|
47
|
Abstract
Optimal search strategies and their implementations in biological systems are a subject of active research. Here we study a search problem which is motivated by the hunt of sperm cells for the egg. We ask for the probability for an active swimmer to find a target under the condition that the swimmer starts at a certain distance from the target. We find that success probability is maximal for a certain level of fluctuations characterized by the persistence length of the swimming path of the swimmer. We derive a scaling law for the optimal persistence length as a function of the initial target distance and search time by mapping the search on a polymer physics problem.
Collapse
|
48
|
Pérez A, Lankas F, Luque FJ, Orozco M. Towards a molecular dynamics consensus view of B-DNA flexibility. Nucleic Acids Res 2008; 36:2379-94. [PMID: 18299282 PMCID: PMC2367714 DOI: 10.1093/nar/gkn082] [Citation(s) in RCA: 130] [Impact Index Per Article: 8.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2007] [Revised: 02/07/2008] [Accepted: 02/08/2008] [Indexed: 01/05/2023] Open
Abstract
We present a systematic study of B-DNA flexibility in aqueous solution using long-scale molecular dynamics simulations with the two more recent versions of nucleic acids force fields (CHARMM27 and parmbsc0) using four long duplexes designed to contain several copies of each individual base pair step. Our study highlights some differences between pambsc0 and CHARMM27 families of simulations, but also extensive agreement in the representation of DNA flexibility. We also performed additional simulations with the older AMBER force fields parm94 and parm99, corrected for non-canonical backbone flips. Taken together, the results allow us to draw for the first time a consensus molecular dynamics picture of B-DNA flexibility.
Collapse
Affiliation(s)
- Alberto Pérez
- Joint IRB-BSC Program on Computational Biology, Institute of Research in Biomedicine, Parc Científic de Barcelona, Josep Samitier 1-5, Barcelona 08028, Barcelona Supercomputing Centre, Jordi Girona 31, Edifici Torre Girona. Barcelona 08034, Departament de Fisicoquímica, Facultat de Farmàcia, Avgda Diagonal sn, Barcelona 08028, Spain, Laboratory for Computation and Visualization in Mathematics and Mechanics, Swiss Federal Institute of Technology (EPFL), CH-1015 Lausanne, Switzerland, Centre for Complex Molecular Systems and Biomolecues, Institute of Organic Chemistry and Biochemistry Flemingovo nam. 2, 166 10 Praha 6, Czech Republic, National Institute of Bioinformatics, Parc Científic de Barcelona, Josep Samitier 1-5 and Departament de Bioquímica, Facultat de Biología, Avgda Diagonal 647, Barcelona 08028, Spain
| | - Filip Lankas
- Joint IRB-BSC Program on Computational Biology, Institute of Research in Biomedicine, Parc Científic de Barcelona, Josep Samitier 1-5, Barcelona 08028, Barcelona Supercomputing Centre, Jordi Girona 31, Edifici Torre Girona. Barcelona 08034, Departament de Fisicoquímica, Facultat de Farmàcia, Avgda Diagonal sn, Barcelona 08028, Spain, Laboratory for Computation and Visualization in Mathematics and Mechanics, Swiss Federal Institute of Technology (EPFL), CH-1015 Lausanne, Switzerland, Centre for Complex Molecular Systems and Biomolecues, Institute of Organic Chemistry and Biochemistry Flemingovo nam. 2, 166 10 Praha 6, Czech Republic, National Institute of Bioinformatics, Parc Científic de Barcelona, Josep Samitier 1-5 and Departament de Bioquímica, Facultat de Biología, Avgda Diagonal 647, Barcelona 08028, Spain
| | - F. Javier Luque
- Joint IRB-BSC Program on Computational Biology, Institute of Research in Biomedicine, Parc Científic de Barcelona, Josep Samitier 1-5, Barcelona 08028, Barcelona Supercomputing Centre, Jordi Girona 31, Edifici Torre Girona. Barcelona 08034, Departament de Fisicoquímica, Facultat de Farmàcia, Avgda Diagonal sn, Barcelona 08028, Spain, Laboratory for Computation and Visualization in Mathematics and Mechanics, Swiss Federal Institute of Technology (EPFL), CH-1015 Lausanne, Switzerland, Centre for Complex Molecular Systems and Biomolecues, Institute of Organic Chemistry and Biochemistry Flemingovo nam. 2, 166 10 Praha 6, Czech Republic, National Institute of Bioinformatics, Parc Científic de Barcelona, Josep Samitier 1-5 and Departament de Bioquímica, Facultat de Biología, Avgda Diagonal 647, Barcelona 08028, Spain
| | - Modesto Orozco
- Joint IRB-BSC Program on Computational Biology, Institute of Research in Biomedicine, Parc Científic de Barcelona, Josep Samitier 1-5, Barcelona 08028, Barcelona Supercomputing Centre, Jordi Girona 31, Edifici Torre Girona. Barcelona 08034, Departament de Fisicoquímica, Facultat de Farmàcia, Avgda Diagonal sn, Barcelona 08028, Spain, Laboratory for Computation and Visualization in Mathematics and Mechanics, Swiss Federal Institute of Technology (EPFL), CH-1015 Lausanne, Switzerland, Centre for Complex Molecular Systems and Biomolecues, Institute of Organic Chemistry and Biochemistry Flemingovo nam. 2, 166 10 Praha 6, Czech Republic, National Institute of Bioinformatics, Parc Científic de Barcelona, Josep Samitier 1-5 and Departament de Bioquímica, Facultat de Biología, Avgda Diagonal 647, Barcelona 08028, Spain
| |
Collapse
|
49
|
|