Parmbsc1: Parameterization and Validation of a new State-of-the-art Force Field for DNA Simulations

Classical force fields are the core of classical simulations, particularly of molecular dynamics (MD), a technique that is changing our view on the structure, flexibility and function of biological macromolecules. Originated from the pioneering work of Lifson’s group in the sixties, force fields hav...

Descripción completa

Detalles Bibliográficos
Autor: Ivani, Iván
Tipo de recurso: tesis doctoral
Estado:Versión publicada
Fecha de publicación:2016
País:España
Institución:CBUC, CESCA
Repositorio:TDR. Tesis Doctorales en Red
OAI Identifier:oai:www.tdx.cat:10803/402737
Acceso en línea:http://hdl.handle.net/10803/402737
Access Level:acceso abierto
Palabra clave:Simulació per ordinador
Simulación por ordenador
Computer simulation
ADN
DNA
Dinàmica molecular
Dinámica molecular
Molecular dynamics
RNA
ARN
Ciències Experimentals i Matemàtiques
577
id ES_65dbe25692647ebe6d7154d6ebfdbb8e
oai_identifier_str oai:www.tdx.cat:10803/402737
network_acronym_str ES
network_name_str España
repository_id_str
dc.title.none.fl_str_mv Parmbsc1: Parameterization and Validation of a new State-of-the-art Force Field for DNA Simulations
title Parmbsc1: Parameterization and Validation of a new State-of-the-art Force Field for DNA Simulations
spellingShingle Parmbsc1: Parameterization and Validation of a new State-of-the-art Force Field for DNA Simulations
Ivani, Iván
Simulació per ordinador
Simulación por ordenador
Computer simulation
ADN
DNA
Dinàmica molecular
Dinámica molecular
Molecular dynamics
RNA
ARN
Ciències Experimentals i Matemàtiques
577
title_short Parmbsc1: Parameterization and Validation of a new State-of-the-art Force Field for DNA Simulations
title_full Parmbsc1: Parameterization and Validation of a new State-of-the-art Force Field for DNA Simulations
title_fullStr Parmbsc1: Parameterization and Validation of a new State-of-the-art Force Field for DNA Simulations
title_full_unstemmed Parmbsc1: Parameterization and Validation of a new State-of-the-art Force Field for DNA Simulations
title_sort Parmbsc1: Parameterization and Validation of a new State-of-the-art Force Field for DNA Simulations
dc.creator.none.fl_str_mv Ivani, Iván
author Ivani, Iván
author_facet Ivani, Iván
author_role author
dc.contributor.none.fl_str_mv Orozco López, Modesto
Orozco López, Modesto
Universitat de Barcelona. Facultat de Biologia
dc.subject.none.fl_str_mv Simulació per ordinador
Simulación por ordenador
Computer simulation
ADN
DNA
Dinàmica molecular
Dinámica molecular
Molecular dynamics
RNA
ARN
Ciències Experimentals i Matemàtiques
577
topic Simulació per ordinador
Simulación por ordenador
Computer simulation
ADN
DNA
Dinàmica molecular
Dinámica molecular
Molecular dynamics
RNA
ARN
Ciències Experimentals i Matemàtiques
577
description Classical force fields are the core of classical simulations, particularly of molecular dynamics (MD), a technique that is changing our view on the structure, flexibility and function of biological macromolecules. Originated from the pioneering work of Lifson’s group in the sixties, force fields have been in continuous evolution, improving in each generation the accuracy in the representation of proteins and nucleic acid. Force field development is tightly connected to the refinement of simulation procedures and to the extension of simulation time scales. Thus, as simulation time passed the microsecond barrier, MD simulations have revealed the existence of some errors in the default force field for DNA simulations, parmbsc0 (developed in the group). The goal of this thesis is to address these problems by a reparameterization of AMBER force field that aims to represent a wide range of DNA structures under physiological and non-physiological conditions. Keeping α/γ parmbsc0 corrections and parm99 non-bonded parameters, we systematically reparameterized sugar puckering, ε, ζ and χ torsions using high level QM calculations both in gas phase and solution. The refined force field has been tested for more than 3 years to an unprecedented level of detail, considering a large variety of DNAs, and analyzing structural, mechanical and dynamical properties of the DNAs resulting from the corresponding MD simulations. The refined force field parameters have been also subjected for more than 1 year of β-testing by different groups, finding to our knowledge no major drawbacks. In the world of RNA simulations, despite the recent efforts to improve the description of RNA in MD simulations, RNA force fields are still far in accuracy from those of DNA. A probable cause could be the incomplete understanding of the mechanism of 2’-OH orientation, which in big extent determines the RNA conformation and most probably serves as the molecular switch. THESIS ORGANIZATION This thesis is compiled of five publications (or in the process of publication) works; first three consider DNA force field development and following validation and benchmark while the last two are focused on RNA efforts. For better understanding of this work Chapter 1 introduces the central concepts related to nucleic acids, their structures and ways to study them. Chapter 2 goes into more details of the methodology employed here, briefly explaining basic QM formalism and MD simulations with an emphasis on force fields. Chapter 3 is a small handbook of methods employed in the analysis in this work. All together first three chapters should provide a solid ground to better understand the details and the relevance of the five publications in the following two chapters. Chapter 4 is based on the development of new force field, called parmbsc1, its further testing on the Drew-Dickerson sequence and benchmarking. Chapter 5 focuses on efforts to understand the mechanism of complexity of RNA structures studying 2’-OH rotation, and computational design of a new RNA dumbbell structure. A summary of the major results and a general discussion that reflects on the five projects and future work are presented in Chapter 6, with the main conclusions at the end of this work.
publishDate 2016
dc.date.none.fl_str_mv 2016
2017
2017
dc.type.none.fl_str_mv info:eu-repo/semantics/doctoralThesis
info:eu-repo/semantics/publishedVersion
format doctoralThesis
status_str publishedVersion
dc.identifier.none.fl_str_mv http://hdl.handle.net/10803/402737
url http://hdl.handle.net/10803/402737
dc.language.none.fl_str_mv Inglés
language_invalid_str_mv Inglés
dc.rights.none.fl_str_mv http://creativecommons.org/licenses/by-nc-nd/4.0/
info:eu-repo/semantics/openAccess
rights_invalid_str_mv http://creativecommons.org/licenses/by-nc-nd/4.0/
eu_rights_str_mv openAccess
dc.format.none.fl_str_mv 341 p.
application/pdf
application/pdf
dc.publisher.none.fl_str_mv Universitat de Barcelona
publisher.none.fl_str_mv Universitat de Barcelona
dc.source.none.fl_str_mv TDX (Tesis Doctorals en Xarxa)
reponame:TDR. Tesis Doctorales en Red
instname:CBUC, CESCA
instname_str CBUC, CESCA
reponame_str TDR. Tesis Doctorales en Red
collection TDR. Tesis Doctorales en Red
repository.name.fl_str_mv
repository.mail.fl_str_mv
_version_ 1869409774110507008
spelling Parmbsc1: Parameterization and Validation of a new State-of-the-art Force Field for DNA SimulationsIvani, IvánSimulació per ordinadorSimulación por ordenadorComputer simulationADNDNADinàmica molecularDinámica molecularMolecular dynamicsRNAARNCiències Experimentals i Matemàtiques577Classical force fields are the core of classical simulations, particularly of molecular dynamics (MD), a technique that is changing our view on the structure, flexibility and function of biological macromolecules. Originated from the pioneering work of Lifson’s group in the sixties, force fields have been in continuous evolution, improving in each generation the accuracy in the representation of proteins and nucleic acid. Force field development is tightly connected to the refinement of simulation procedures and to the extension of simulation time scales. Thus, as simulation time passed the microsecond barrier, MD simulations have revealed the existence of some errors in the default force field for DNA simulations, parmbsc0 (developed in the group). The goal of this thesis is to address these problems by a reparameterization of AMBER force field that aims to represent a wide range of DNA structures under physiological and non-physiological conditions. Keeping α/γ parmbsc0 corrections and parm99 non-bonded parameters, we systematically reparameterized sugar puckering, ε, ζ and χ torsions using high level QM calculations both in gas phase and solution. The refined force field has been tested for more than 3 years to an unprecedented level of detail, considering a large variety of DNAs, and analyzing structural, mechanical and dynamical properties of the DNAs resulting from the corresponding MD simulations. The refined force field parameters have been also subjected for more than 1 year of β-testing by different groups, finding to our knowledge no major drawbacks. In the world of RNA simulations, despite the recent efforts to improve the description of RNA in MD simulations, RNA force fields are still far in accuracy from those of DNA. A probable cause could be the incomplete understanding of the mechanism of 2’-OH orientation, which in big extent determines the RNA conformation and most probably serves as the molecular switch. THESIS ORGANIZATION This thesis is compiled of five publications (or in the process of publication) works; first three consider DNA force field development and following validation and benchmark while the last two are focused on RNA efforts. For better understanding of this work Chapter 1 introduces the central concepts related to nucleic acids, their structures and ways to study them. Chapter 2 goes into more details of the methodology employed here, briefly explaining basic QM formalism and MD simulations with an emphasis on force fields. Chapter 3 is a small handbook of methods employed in the analysis in this work. All together first three chapters should provide a solid ground to better understand the details and the relevance of the five publications in the following two chapters. Chapter 4 is based on the development of new force field, called parmbsc1, its further testing on the Drew-Dickerson sequence and benchmarking. Chapter 5 focuses on efforts to understand the mechanism of complexity of RNA structures studying 2’-OH rotation, and computational design of a new RNA dumbbell structure. A summary of the major results and a general discussion that reflects on the five projects and future work are presented in Chapter 6, with the main conclusions at the end of this work.Universitat de BarcelonaOrozco López, ModestoOrozco López, ModestoUniversitat de Barcelona. Facultat de Biologia201720172016info:eu-repo/semantics/doctoralThesisinfo:eu-repo/semantics/publishedVersion341 p.application/pdfapplication/pdfhttp://hdl.handle.net/10803/402737TDX (Tesis Doctorals en Xarxa)reponame:TDR. Tesis Doctorales en Redinstname:CBUC, CESCAInglésL'accés als continguts d'aquesta tesi queda condicionat a l'acceptació de les condicions d'ús establertes per la següent llicència Creative Commons: http://creativecommons.org/licenses/by-nc-nd/4.0/http://creativecommons.org/licenses/by-nc-nd/4.0/info:eu-repo/semantics/openAccessoai:www.tdx.cat:10803/4027372026-06-14T12:46:07Z
score 15,300719