1
|
Rezaei M, Sakong S, Groß A. Sodium Triflate Water-in-Salt Electrolytes in Advanced Battery Applications: A First-Principles-Based Molecular Dynamics Study. ACS APPLIED MATERIALS & INTERFACES 2024; 16:32169-32188. [PMID: 38862108 PMCID: PMC11212028 DOI: 10.1021/acsami.4c01449] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/25/2024] [Revised: 04/03/2024] [Accepted: 05/29/2024] [Indexed: 06/13/2024]
Abstract
Offering a compelling combination of safety and cost-effectiveness, water-in-salt (WiS) electrolytes have emerged as promising frontiers in energy storage technology. Still, there is a strong demand for research and development efforts to make these electrolytes ripe for commercialization. Here, we present a first-principles-based molecular dynamics (MD) study addressing in detail the properties of a sodium triflate WiS electrolyte for Na-ion batteries. We have developed a workflow based on a machine learning (ML) potential derived from ab initio MD simulations. As ML potentials are typically restricted to the interpolation of the data points of the training set and have hardly any predictive properties, we subsequently optimize a classical force field based on physics principles to ensure broad applicability and high performance. Performing and analyzing detailed MD simulations, we identify several very promising properties of the sodium triflate as a WiS electrolyte but also indicate some potential stability challenges associated with its use as a battery electrolyte.
Collapse
Affiliation(s)
- Majid Rezaei
- Institute
of Theoretical Chemistry, Ulm University, Oberberghof 7, 89081 Ulm, Germany
| | - Sung Sakong
- Institute
of Theoretical Chemistry, Ulm University, Oberberghof 7, 89081 Ulm, Germany
| | - Axel Groß
- Institute
of Theoretical Chemistry, Ulm University, Oberberghof 7, 89081 Ulm, Germany
- Helmholtz
Institute Ulm (HIU) for Electrochemical Energy Storage, Helmholtzstraße 11, 89069 Ulm, Germany
| |
Collapse
|
2
|
Khavani M, Mehranfar A, Mofrad MRK. On the potentials of sialic acid derivatives as inhibitors for the mumps virus: A molecular dynamics and quantum chemistry investigation. Virus Res 2023; 326:199050. [PMID: 36682462 DOI: 10.1016/j.virusres.2023.199050] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2022] [Revised: 12/21/2022] [Accepted: 01/18/2023] [Indexed: 01/22/2023]
Abstract
Mumps virus is an infectious pathogen causing major health problems for humans such as encephalitis, orchitis, and parotitis. Therefore, designing an inhibitor for this virus is of great medical and public health importance. With this goal in mind, we investigate the affinity of different sialic acid-based compounds (ligands) against the hemagglutinin-neuraminidase (HN) protein of the mumps virus, using a combination of molecular dynamics (MD) simulations and quantum chemistry calculations. Our MD simulation results indicate that the ligands form stable complexes with the HN protein through a combination of electrostatic, van der Waals (vdW), and hydrogen bond (H-bond) interactions, which the electrostatic interactions play a more important role in the complexation process. Based on the obtained results from the structural analysis Arg381, Arg291, and Arg49 play a key role in the binding site interactions with the different ligands, in comparison with other residues. There are some candidates such as Neu5Acα2-6Galβ1-4GlcNAcβ, Neu5Acα2-3Galβ1-3GlcNacβ1-3Galβ1-4Glc, and Neu5Acα2-6Galβ1-4GlcNAcβ1-3Galβ1-4Glc that form more stable complexes with the HN than the α2-3-Sialyllactose confirmed by the calculated Gibbs binding energies (-39.65, -46.93, and -36.49 kcal.mol-1, respectively). To investigate the relationship between the molecular properties of the selected compounds and their affinity to the HN receptor, density functional theory dispersion corrected (DFT-D3) calculations were employed. According to our DFT-D3 results, neutral sialic acid-based compounds have lower reactivity to the mumps virus than the negativity charge structures. Moreover, by increasing the electronic chemical potential (μ) the vdW and H-bond interactions between drugs and the HN protein increase. In other words, by elevating the electron tendency of the selected ligands their affinity to the mumps virus increases. Our quantum chemistry calculations reveal that in addition to the structural features the molecular properties of the drugs can play important roles in their affinity and reactivity against the virus. The results of this study can provide useful details to design new compounds or improve their properties against the mumps virus.
Collapse
Affiliation(s)
- Mohammad Khavani
- Departments of Bioengineering and Mechanical Engineering, Molecular Cell Biomechanics Laboratory, University of California Berkeley, Berkeley, CA 94720, USA
| | - Aliyeh Mehranfar
- Departments of Bioengineering and Mechanical Engineering, Molecular Cell Biomechanics Laboratory, University of California Berkeley, Berkeley, CA 94720, USA
| | - Mohammad R K Mofrad
- Departments of Bioengineering and Mechanical Engineering, Molecular Cell Biomechanics Laboratory, University of California Berkeley, Berkeley, CA 94720, USA.
| |
Collapse
|
3
|
Khavani M, Mehranfar A, Mofrad MRK. On the Sensitivity and Affinity of Gold, Silver, and Platinum Surfaces against the SARS-CoV-2 Virus: A Comparative Computational Study. J Chem Inf Model 2023; 63:1276-1292. [PMID: 36735895 PMCID: PMC9924083 DOI: 10.1021/acs.jcim.2c01378] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2022] [Indexed: 02/05/2023]
Abstract
The novel coronavirus disease and its complications have motivated the design of new sensors with the highest sensitivity, and affinity for the detection of the SARS-CoV-2 virus is considered in many research studies. In this research article, we employ full atomistic molecular dynamics (MD) models to study the interactions between the receptor binding domain (RBD) and spike protein of the coronavirus and different metals such as gold (Au), platinum (Pt), and silver (Ag) to analyze their sensitivity against this virus. The comparison between the RBD interactions with ACE2 (angiotensin-converting enzyme 2) and different metals indicates that metals have remarkable effects on the structural features and dynamical properties of the RBD. The binding site of the RBD has more affinity to the surfaces of gold, platinum, and silver than to the other parts of the protein. Moreover, the initial configuration of the RBD relative to the metal surface plays an important role in the stability of metal complexes with the RBD. The binding face of the protein to the metal surface has been changed in the presence of different metals. In other words, the residues of the RBD that participate in RBD interactions with the metals are different irrespective of the initial configurations in which the [Asn, Thr, Tyr], [Ser, Thr, Tyr], and [Asn, Asp, Tyr] residues of the protein have a greater affinity to Ag, Au, and Pt, respectively. The corresponding metals have a considerable affinity to the RBD, which due to strong interactions with the protein can change the secondary structure and structural features. Based on the obtained results during the complexation process between the protein and metals, the helical structure of the protein changes to the bend and antiparallel β-sheets. The calculated binding energies for the RBD complexes with silver, gold, and platinum are -95.03, -138.03, and -133.96 kcal·mol-1, respectively. The adsorption process of the spike protein on the surfaces of different metals represents similar results and indicates that the entire spike protein of the coronavirus forms a more stable complex with the gold surface compared with other metals. Moreover, the RBD of the spike protein has more interactions with the surfaces than with the other parts of the protein. Therefore, it is possible to predict the properties of the coronavirus on the metal surface based on the dynamical behavior of the RBD. Overall, our computational results confirm that the gold surface can be considered as an outstanding substrate for developing new sensors with the highest sensitivity against SARS-CoV-2.
Collapse
Affiliation(s)
- Mohammad Khavani
- Molecular Cell Biomechanics Laboratory, Departments of Bioengineering and
Mechanical Engineering, University of California Berkeley,
Berkeley, California94720, United States
| | - Aliyeh Mehranfar
- Molecular Cell Biomechanics Laboratory, Departments of Bioengineering and
Mechanical Engineering, University of California Berkeley,
Berkeley, California94720, United States
| | - Mohammad R. K. Mofrad
- Molecular Cell Biomechanics Laboratory, Departments of Bioengineering and
Mechanical Engineering, University of California Berkeley,
Berkeley, California94720, United States
| |
Collapse
|
4
|
Wilson KA, Jeong YER, Wetmore SD. Multiscale computational investigations of the translesion synthesis bypass of tobacco-derived DNA adducts: critical insights that complement experimental biochemical studies. Phys Chem Chem Phys 2022; 24:10667-10683. [PMID: 35502640 DOI: 10.1039/d2cp00481j] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2023]
Abstract
Among the numerous agents that damage DNA, tobacco products remain one of the most lethal and result in the most diverse set of DNA lesions. This perspective aims to provide an overview of computational work conducted to complement experimental biochemical studies on the mutagenicity of adducts derived from the most potent tobacco carcinogen, namely 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (nicotine-derived nitrosaminoketone or NNK). Lesions ranging from the smallest methylated thymine derivatives to the larger, flexible pyridyloxobutyl (POB) guanine adducts are considered. Insights are obtained from density functional theory (DFT) calculations and molecular dynamics (MD) simulations into the damaged nucleobase and nucleoside structures, the accommodation of the lesions in the active site of key human polymerases, the intrinsic base pairing potentials of the adducts, and dNTP incorporation opposite the lesions. Overall, the computational data provide atomic level information that can rationalize the differential mutagenic properties of tobacco-derived lesions and uncover important insights into the impact of adduct size, nucleobase, position, and chemical composition of the bulky moiety.
Collapse
Affiliation(s)
- Katie A Wilson
- Department of Chemistry and Biochemistry, Alberta RNA Research and Training Institute (ARRTI) and Southern Alberta Genome Sciences Center (SAGSC), University of Lethbridge, 4401 University Drive West, Lethbridge, Alberta, T1K 3M4, Canada.
| | - Ye Eun Rebecca Jeong
- Department of Chemistry and Biochemistry, Alberta RNA Research and Training Institute (ARRTI) and Southern Alberta Genome Sciences Center (SAGSC), University of Lethbridge, 4401 University Drive West, Lethbridge, Alberta, T1K 3M4, Canada.
| | - Stacey D Wetmore
- Department of Chemistry and Biochemistry, Alberta RNA Research and Training Institute (ARRTI) and Southern Alberta Genome Sciences Center (SAGSC), University of Lethbridge, 4401 University Drive West, Lethbridge, Alberta, T1K 3M4, Canada.
| |
Collapse
|
5
|
Mehranfar A, Khavani M, Izadyar M. A molecular dynamic study on the ability of phosphorene for designing new sensor for SARS-CoV-2 detection. J Mol Liq 2022; 345:117852. [PMID: 34690390 PMCID: PMC8520178 DOI: 10.1016/j.molliq.2021.117852] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2021] [Revised: 10/11/2021] [Accepted: 10/13/2021] [Indexed: 12/04/2022]
Abstract
Due to the dramatic increase in the number of patients with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), designing new selective and sensitive sensors for the detection of this virus is of importance. In this research, by employing full atomistic molecular dynamics (MD) simulations, the interactions of the receptor-binding domain (RBD) of the SARS-CoV-2 with phosphorene and graphene nanosheets were analyzed to investigate their sensing ability against this protein. Based on the obtained results, the RBD interactions with the surface of graphene and phosphorene nanosheets do not have important effects on the folding properties of the RBD but this protein has unique dynamical behavior against each nanostructure. In the presence of graphene and phosphorene, the RBD has lower stability because due to the strong interactions between RBD and these nanostructures. This protein spreads on the surface and has lower structural compaction, but in comparison with graphene, RBD shows greater stability on the surface of the phosphorene nanosheet. Moreover, RBD forms a more stable complex with phosphorene nanosheet in comparison with graphene due to greater electrostatic and van der Waals interactions. The calculated Gibbs binding energy for the RBD complexation process with phosphorene and graphene are −200.37 and −83.65 kcal mol−1, respectively confirming that phosphorene has higher affinity and sensitivity against this protein than graphene. Overall, the obtained results confirm that phosphorene can be a good candidate for designing new nanomaterials for selective detection of SARS-CoV-2.
Collapse
Affiliation(s)
- Aliyeh Mehranfar
- Research Center for Modeling and Computational Sciences, Department of Chemistry, Faculty of Science, Ferdowsi University of Mashhad, Mashhad, Iran
| | - Mohammad Khavani
- Department of Chemistry and Materials Science, School of Chemical Engineering, Aalto University, P.O. Box 16100, FI-00076 Aalto, Finland
| | - Mohammad Izadyar
- Research Center for Modeling and Computational Sciences, Department of Chemistry, Faculty of Science, Ferdowsi University of Mashhad, Mashhad, Iran
| |
Collapse
|
6
|
Salsbury AM, Dean TJ, Lemkul JA. Polarizable Molecular Dynamics Simulations of Two c-kit Oncogene Promoter G-Quadruplexes: Effect of Primary and Secondary Structure on Loop and Ion Sampling. J Chem Theory Comput 2020; 16:3430-3444. [PMID: 32307997 DOI: 10.1021/acs.jctc.0c00191] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
G-quadruplexes (GQs) are highly ordered nucleic acid structures that play fundamental roles in regulating gene expression and maintaining genomic stability. GQs are topologically diverse and enriched in promoter sequences of growth regulatory genes and proto-oncogenes, suggesting that they may serve as attractive targets for drug design at the level of transcription rather than inhibiting the activity of the protein products of these genes. The c-kit promoter contains three adjacent GQ-forming sequences that have proposed antagonistic effects on gene expression and thus are promising drug targets for diseases such as gastrointestinal stromal tumors, mast cell disease, and leukemia. Because GQ stability is influenced by primary structure, secondary structure, and ion interactions, a greater understanding of GQ structure, dynamics, and ion binding properties is needed to develop novel, GQ-targeting therapeutics. Here, we performed molecular dynamics simulations to systematically study the c-kit2 and c-kit* GQs, evaluating nonpolarizable and polarizable force fields (FFs) and examining the effects of base substitutions and cation type (K+, Na+, and Li+) on the dynamics of their isolated and linked structures. We found that the Drude polarizable FF outperformed the additive CHARMM36 FF in two- and three-tetrad GQs and solutions of KCl, NaCl, and LiCl. Drude simulations with different cations agreed with the known GQ stabilization preference (K+ > Na+ > Li+) and illustrated that tetrad core-ion coordination differs as a function of cation type. Finally, we showed that differences in primary and secondary structure influence loop sampling, ion binding, and core-ion energetics of GQs.
Collapse
Affiliation(s)
- Alexa M Salsbury
- Department of Biochemistry and Center for Drug Discovery, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Tanner J Dean
- Department of Biochemistry and Center for Drug Discovery, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Justin A Lemkul
- Department of Biochemistry and Center for Drug Discovery, Virginia Tech, Blacksburg, Virginia 24061, United States
| |
Collapse
|
7
|
Okamoto D, Watanabe Y, Yoshida N, Nakano H. Implementation of state-averaged MCSCF method to RISM- and 3D-RISM-SCF schemes. Chem Phys Lett 2019. [DOI: 10.1016/j.cplett.2019.05.051] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
8
|
Galindo-Murillo R, Cheatham TE. Lessons learned in atomistic simulation of double-stranded DNA: Solvation and salt concerns [Article v1.0]. ACTA ACUST UNITED AC 2019; 1. [PMID: 33073182 DOI: 10.33011/livecoms.1.2.9974] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Nucleic acids are highly charged macromolecules sensitive to their surroundings of water, salt, and other biomolecules. Molecular dynamics simulations with accurate biomolecular force fields provide a detailed atomistic view into DNA and RNA that has been useful to study the structure and dynamics of these molecules and their biological relevance. In this work we study the Drew-Dickerson dodecamer duplex with the sequence d(GCGCAATTGCGC)2 in three different salt concentrations and using different monvalent salt types to detect possible structural influence. Overall, the DNA shows no major structural changes regardless of amount or type of monovalent ions used. Our results show that only at very high salt conditions (5M) is a small structural effect observed in the DNA duplex, which mainly consist of narrowing of the grooves due to increased residence of ions. We also present the importance of sampling time to achieve a converged ensemble, which is of major relevance in any simulation to avoid biased or non-meaningful results.
Collapse
Affiliation(s)
- Rodrigo Galindo-Murillo
- Department of Medicinal Chemistry, L. S. Skaggs Pharmacy Institute, University of Utah, Salt Lake City, UT 84112
| | - Thomas E Cheatham
- Department of Medicinal Chemistry, L. S. Skaggs Pharmacy Institute, University of Utah, Salt Lake City, UT 84112
| |
Collapse
|
9
|
Thermodynamics of helix formation in small peptides of varying length in vacuo, in implicit solvent, and in explicit solvent. J Mol Model 2018; 25:3. [DOI: 10.1007/s00894-018-3886-2] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2018] [Accepted: 11/28/2018] [Indexed: 10/27/2022]
|
10
|
Jasiński M, Kulik M, Wojciechowska M, Stolarski R, Trylska J. Interactions of 2'-O-methyl oligoribonucleotides with the RNA models of the 30S subunit A-site. PLoS One 2018; 13:e0191138. [PMID: 29351348 PMCID: PMC5774723 DOI: 10.1371/journal.pone.0191138] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2017] [Accepted: 12/28/2017] [Indexed: 12/15/2022] Open
Abstract
Synthetic oligonucleotides targeting functional regions of the prokaryotic rRNA could be promising antimicrobial agents. Indeed, such oligonucleotides were proven to inhibit bacterial growth. 2’-O-methylated (2’-O-Me) oligoribonucleotides with a sequence complementary to the decoding site in 16S rRNA were reported as inhibitors of bacterial translation. However, the binding mode and structures of the formed complexes, as well as the level of selectivity of the oligonucleotides between the prokaryotic and eukaryotic target, were not determined. We have analyzed three 2’-O-Me oligoribonucleotides designed to hybridize with the models of the prokaryotic rRNA containing two neighboring aminoglycoside binding pockets. One pocket is the paromomycin/kanamycin binding site corresponding to the decoding site in the small ribosomal subunit and the other one is the close-by hygromycin B binding site whose dynamics has not been previously reported. Molecular dynamics (MD) simulations, as well as isothermal titration calorimetry, gel electrophoresis and spectroscopic studies have shown that the eukaryotic rRNA model is less conformationally stable (in terms of hydrogen bonds and stacking interactions) than the corresponding prokaryotic one. In MD simulations of the eukaryotic construct, the nucleotide U1498, which plays an important role in correct positioning of mRNA during translation, is flexible and spontaneously flips out into the solvent. In solution studies, the 2’-O-Me oligoribonucleotides did not interact with the double stranded rRNA models but all formed stable complexes with the single-stranded prokaryotic target. 2’-O-Me oligoribonucleotides with one and two mismatches bound less tightly to the eukaryotic target. This shows that at least three mismatches between the 2’-O-Me oligoribonucleotide and eukaryotic rRNA are required to ensure target selectivity. The results also suggest that, in the ribosome environment, the strand invasion is the preferred binding mode of 2’-O-Me oligoribonucleotides targeting the aminoglycoside binding sites in 16S rRNA.
Collapse
Affiliation(s)
- Maciej Jasiński
- Centre of New Technologies, University of Warsaw, Warsaw, Poland
- College of Inter-Faculty Individual Studies in Mathematics and Natural Sciences, University of Warsaw, Warsaw, Poland
| | - Marta Kulik
- Centre of New Technologies, University of Warsaw, Warsaw, Poland
| | | | - Ryszard Stolarski
- Department of Biophysics, Institute of Experimental Physics, Faculty of Physics, University of Warsaw, Warsaw, Poland
| | - Joanna Trylska
- Centre of New Technologies, University of Warsaw, Warsaw, Poland
- * E-mail:
| |
Collapse
|
11
|
Ivanović MT, Bruetzel LK, Shevchuk R, Lipfert J, Hub JS. Quantifying the influence of the ion cloud on SAXS profiles of charged proteins. Phys Chem Chem Phys 2018; 20:26351-26361. [DOI: 10.1039/c8cp03080d] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
MD simulations and Poisson–Boltzmann calculations predict ion cloud effects on SAXS experiments.
Collapse
Affiliation(s)
- Miloš T. Ivanović
- Georg-August-Universität Göttingen, Institute for Microbiology and Genetics
- 37077 Göttingen
- Germany
| | - Linda K. Bruetzel
- Ludwig-Maximilian-Universität München, Department of Physics
- 80799 München
- Germany
| | - Roman Shevchuk
- Georg-August-Universität Göttingen, Institute for Microbiology and Genetics
- 37077 Göttingen
- Germany
| | - Jan Lipfert
- Ludwig-Maximilian-Universität München, Department of Physics
- 80799 München
- Germany
| | - Jochen S. Hub
- Georg-August-Universität Göttingen, Institute for Microbiology and Genetics
- 37077 Göttingen
- Germany
| |
Collapse
|
12
|
Choi JH, Choi HR, Jeon J, Cho M. Ion aggregation in high salt solutions. VII. The effect of cations on the structures of ion aggregates and water hydrogen-bonding network. J Chem Phys 2017; 147:154107. [DOI: 10.1063/1.4993479] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Jun-Ho Choi
- Center for Molecular Spectroscopy and Dynamics, Institute for Basic Science, Korea University, Seoul 02841, South Korea
| | - Hyung Ran Choi
- Center for Molecular Spectroscopy and Dynamics, Institute for Basic Science, Korea University, Seoul 02841, South Korea
- Department of Chemistry, Korea University, Seoul 02841, South Korea
| | - Jonggu Jeon
- Center for Molecular Spectroscopy and Dynamics, Institute for Basic Science, Korea University, Seoul 02841, South Korea
| | - Minhaeng Cho
- Center for Molecular Spectroscopy and Dynamics, Institute for Basic Science, Korea University, Seoul 02841, South Korea
- Department of Chemistry, Korea University, Seoul 02841, South Korea
| |
Collapse
|
13
|
Havrila M, Stadlbauer P, Islam B, Otyepka M, Šponer J. Effect of Monovalent Ion Parameters on Molecular Dynamics Simulations of G-Quadruplexes. J Chem Theory Comput 2017; 13:3911-3926. [PMID: 28657760 DOI: 10.1021/acs.jctc.7b00257] [Citation(s) in RCA: 42] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
G-quadruplexes (GQs) are key noncanonical DNA and RNA architectures stabilized by desolvated monovalent cations present in their central channels. We analyze extended atomistic molecular dynamics simulations (∼580 μs in total) of GQs with 11 monovalent cation parametrizations, assessing GQ overall structural stability, dynamics of internal cations, and distortions of the G-tetrad geometries. Majority of simulations were executed with the SPC/E water model; however, test simulations with TIP3P and OPC water models are also reported. The identity and parametrization of ions strongly affect behavior of a tetramolecular d[GGG]4 GQ, which is unstable with several ion parametrizations. The remaining studied RNA and DNA GQs are structurally stable, though the G-tetrad geometries are always deformed by bifurcated H-bonding in a parametrization-specific manner. Thus, basic 10-μs-scale simulations of fully folded GQs can be safely done with a number of cation parametrizations. However, there are parametrization-specific differences and basic force-field errors affecting the quantitative description of ion-tetrad interactions, which may significantly affect studies of the ion-binding processes and description of the GQ folding landscape. Our d[GGG]4 simulations indirectly suggest that such studies will also be sensitive to the water models. During exchanges with bulk water, the Na+ ions move inside the GQs in a concerted manner, while larger relocations of the K+ ions are typically separated. We suggest that the Joung-Cheatham SPC/E K+ parameters represent a safe choice in simulation studies of GQs, though variation of ion parameters can be used for specific simulation goals.
Collapse
Affiliation(s)
- Marek Havrila
- Institute of Biophysics, Academy of Sciences of the Czech Republic , Královopolská 135, 612 65 Brno, Czech Republic.,CEITEC - Central European Institute of Technology, Masaryk University , Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Petr Stadlbauer
- Institute of Biophysics, Academy of Sciences of the Czech Republic , Královopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , 17. listopadu 12, 771 46 Olomouc, Czech Republic
| | - Barira Islam
- Institute of Biophysics, Academy of Sciences of the Czech Republic , Královopolská 135, 612 65 Brno, Czech Republic
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , 17. listopadu 12, 771 46 Olomouc, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic , Královopolská 135, 612 65 Brno, Czech Republic.,CEITEC - Central European Institute of Technology, Masaryk University , Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| |
Collapse
|
14
|
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
|
15
|
Atzori A, Liggi S, Laaksonen A, Porcu M, Lyubartsev AP, Saba G, Mocci F. Base sequence specificity of counterion binding to DNA: what can MD simulations tell us? CAN J CHEM 2016. [DOI: 10.1139/cjc-2016-0296] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Nucleic acids are highly charged biopolymers whose secondary structure is strongly dependent on electrostatic interactions. Solvent molecules and ions are also believed to play an important role in mediating and directing both sequence recognition and interactions with other molecules, such as proteins and a variety of ligands. Therefore, to fully understand the biological functions of DNA, it is necessary to understand the interactions with the surrounding counterions. It is well known that monovalent counterions can bind to the minor groove of DNA with consecutive sequences of four, or more, adenine and thymine (A-tracts) with relatively long residence times. However, much less is known about their binding to the backbone and to the major groove. In this work, we used molecular dynamics simulations to both investigate the interactions between the backbone and major groove of DNA and one of its physiological counterions (Na+) and evaluate the relationship between these interactions and the nucleotide sequence. Three dodecamers, namely CGAAAATTTTCG, CGCTCTAGAGCG, and CGCGAATTCGCG, were simulated using the Toukan–Rahman flexible SPC water model and Smith and Dang parameters for Na+, revealing a significant sequence dependence on the ion binding to both backbone and major groove. In the absence of experimental data on the atomistic details of the studied interactions, the reliability of the results was evaluated performing the simulations with additional sets of potential parameters for ions and solvent, namely the Ȧqvist or the Joung and Cheatham ion parameters and the TIP3P water model. This allowed us to evaluate the results by verifying which features are preserved independently from the parameters adopted.
Collapse
Affiliation(s)
- Alessio Atzori
- Dipartimento di Scienze Chimiche e Geologiche, Università di Cagliari, Cittadella Universitaria di Monserrato, I-09042 Monserrato (CA), Italy
| | - Sonia Liggi
- Dipartimento di Scienze Chimiche e Geologiche, Università di Cagliari, Cittadella Universitaria di Monserrato, I-09042 Monserrato (CA), Italy
| | - Aatto Laaksonen
- Dipartimento di Scienze Chimiche e Geologiche, Università di Cagliari, Cittadella Universitaria di Monserrato, I-09042 Monserrato (CA), Italy
- Division of Physical Chemistry, Arrhenius Laboratory, Stockholm University, S-10691 Stockholm, Sweden
- Science for Life Laboratory, 171 21 Solna, Sweden
| | - Massimiliano Porcu
- Dipartimento di Scienze Chimiche e Geologiche, Università di Cagliari, Cittadella Universitaria di Monserrato, I-09042 Monserrato (CA), Italy
- Division of Physical Chemistry, Arrhenius Laboratory, Stockholm University, S-10691 Stockholm, Sweden
| | - Alexander P. Lyubartsev
- Division of Physical Chemistry, Arrhenius Laboratory, Stockholm University, S-10691 Stockholm, Sweden
| | - Giuseppe Saba
- Dipartimento di Scienze Chimiche e Geologiche, Università di Cagliari, Cittadella Universitaria di Monserrato, I-09042 Monserrato (CA), Italy
| | - Francesca Mocci
- Dipartimento di Scienze Chimiche e Geologiche, Università di Cagliari, Cittadella Universitaria di Monserrato, I-09042 Monserrato (CA), Italy
- Division of Physical Chemistry, Arrhenius Laboratory, Stockholm University, S-10691 Stockholm, Sweden
- Science for Life Laboratory, 171 21 Solna, Sweden
| |
Collapse
|
16
|
Choi JH, Cho M. Ion aggregation in high salt solutions. VI. Spectral graph analysis of chaotropic ion aggregates. J Chem Phys 2016; 145:174501. [DOI: 10.1063/1.4966246] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jun-Ho Choi
- Center for Molecular Spectroscopy and Dynamics, Institute for Basic Science, Korea University, Seoul 02841, South Korea
- Department of Chemistry, Korea University, Seoul 02841, South Korea
| | - Minhaeng Cho
- Center for Molecular Spectroscopy and Dynamics, Institute for Basic Science, Korea University, Seoul 02841, South Korea
- Department of Chemistry, Korea University, Seoul 02841, South Korea
| |
Collapse
|
17
|
Öztürk MA, Pachov GV, Wade RC, Cojocaru V. Conformational selection and dynamic adaptation upon linker histone binding to the nucleosome. Nucleic Acids Res 2016; 44:6599-613. [PMID: 27270081 PMCID: PMC5001602 DOI: 10.1093/nar/gkw514] [Citation(s) in RCA: 35] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2016] [Revised: 05/06/2016] [Accepted: 05/30/2016] [Indexed: 01/25/2023] Open
Abstract
Linker histones are essential for DNA compaction in chromatin. They bind to nucleosomes in a 1:1 ratio forming chromatosomes. Alternative configurations have been proposed in which the globular domain of the linker histone H5 (gH5) is positioned either on- or off-dyad between the nucleosomal and linker DNAs. However, the dynamic pathways of chromatosome assembly remain elusive. Here, we studied the conformational plasticity of gH5 in unbound and off-dyad nucleosome-bound forms with classical and accelerated molecular dynamics simulations. We find that the unbound gH5 converts between open and closed conformations, preferring the closed form. However, the open gH5 contributes to a more rigid chromatosome and restricts the motion of the nearby linker DNA through hydrophobic interactions with thymidines. Moreover, the closed gH5 opens and reorients in accelerated simulations of the chromatosome. Brownian dynamics simulations of chromatosome assembly, accounting for a range of amplitudes of nucleosome opening and different nucleosome DNA sequences, support the existence of both on- and off-dyad binding modes of gH5 and reveal alternative, sequence and conformation-dependent chromatosome configurations. Taken together, these findings suggest that the conformational dynamics of linker histones and nucleosomes facilitate alternative chromatosome configurations through an interplay between induced fit and conformational selection.
Collapse
Affiliation(s)
- Mehmet Ali Öztürk
- Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies (HITS), Heidelberg 69118, Germany The Hartmut Hoffmann-Berling International Graduate School of Molecular and Cellular Biology (HBIGS), Heidelberg University, Heidelberg 69120, Germany
| | - Georgi V Pachov
- Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies (HITS), Heidelberg 69118, Germany
| | - Rebecca C Wade
- Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies (HITS), Heidelberg 69118, Germany Center for Molecular Biology (ZMBH), DKFZ-ZMBH Alliance, Heidelberg University, Heidelberg 69120, Germany Interdisciplinary Center for Scientific Computing (IWR), Heidelberg 69120, Germany
| | - Vlad Cojocaru
- Computational Structural Biology Laboratory, Department of Cellular and Developmental Biology, Max Planck Institute for Molecular Biomedicine, Münster 48149, Germany Center for Multiscale Theory and Computation, Westfälische Wilhelms University, Münster 48149, Germany
| |
Collapse
|
18
|
Cruz FJ, Mota JP. Conformational Thermodynamics of DNA Strands in Hydrophilic Nanopores. J Phys Chem B 2016. [DOI: 10.1021/acs.jpcb.6b06234] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
19
|
Rebič M, Laaksonen A, Šponer J, Uličný J, Mocci F. Molecular Dynamics Simulation Study of Parallel Telomeric DNA Quadruplexes at Different Ionic Strengths: Evaluation of Water and Ion Models. J Phys Chem B 2016; 120:7380-91. [PMID: 27379924 DOI: 10.1021/acs.jpcb.6b06485] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Most molecular dynamics (MD) simulations of DNA quadruplexes have been performed under minimal salt conditions using the Åqvist potential parameters for the cation with the TIP3P water model. Recently, this combination of parameters has been reported to be problematic for the stability of quadruplex DNA, especially caused by the ion interactions inside or near the quadruplex channel. Here, we verify how the choice of ion parameters and water model can affect the quadruplex structural stability and the interactions with the ions outside the channel. We have performed a series of MD simulations of the human full-parallel telomeric quadruplex by neutralizing its negative charge with K(+) ions. Three combinations of different cation potential parameters and water models have been used: (a) Åqvist ion parameters, TIP3P water model; (b) Joung and Cheatham ion parameters, TIP3P water model; and (c) Joung and Cheatham ion parameters, TIP4Pew water model. For the combinations (b) and (c), the effect of the ionic strength has been evaluated by adding increasing amounts of KCl salt (50, 100, and 200 mM). Two independent simulations using the Åqvist parameters with the TIP3P model show that this combination is clearly less suited for the studied quadruplex with K(+) as counterions. In both simulations, one ion escapes from the channel, followed by significant deformation of the structure, leading to deviating conformation compared to that in the reference crystallographic data. For the other combinations of ion and water potentials, no tendency is observed for the channel ions to escape from the quadruplex channel. In addition, the internal mobility of the three loops, torsion angles, and counterion affinity have been investigated at varied salt concentrations. In summary, the selection of ion and water models is crucial as it can affect both the structure and dynamics as well as the interactions of the quadruplex with its counterions. The results obtained with the TIP4Pew model are found to be closest to the experimental data at all of the studied ion concentrations.
Collapse
Affiliation(s)
- Matúš Rebič
- Department of Materials and Environmental Chemistry, Division of Physical Chemistry, Arrhenius Laboratory, Stockholm University , 10691 Stockholm, Sweden.,Science for Life Laboratory (SciLifelab) , 17121 Solna, Sweden
| | - Aatto Laaksonen
- Department of Materials and Environmental Chemistry, Division of Physical Chemistry, Arrhenius Laboratory, Stockholm University , 10691 Stockholm, Sweden.,Science for Life Laboratory (SciLifelab) , 17121 Solna, Sweden.,Stellenbosch Institute of Advanced Study (STIAS), Wallenberg Research Centre at Stellenbosch University , 7600 Stellenbosch, South Africa.,Department of Chemical and Geological Sciences, University of Cagliari , I-09042 Monserrato, Italy
| | - Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic (AVČR) , Kralovopolska 135, 612 65 Brno, Czech Republic.,Central European Institute of Technology (CEITEC) , Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | | | - Francesca Mocci
- Department of Materials and Environmental Chemistry, Division of Physical Chemistry, Arrhenius Laboratory, Stockholm University , 10691 Stockholm, Sweden.,Science for Life Laboratory (SciLifelab) , 17121 Solna, Sweden.,Department of Chemical and Geological Sciences, University of Cagliari , I-09042 Monserrato, Italy
| |
Collapse
|
20
|
Choi JH, Cho M. Ion aggregation in high salt solutions. V. Graph entropy analyses of ion aggregate structure and water hydrogen bonding network. J Chem Phys 2016; 144:204126. [DOI: 10.1063/1.4952648] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Affiliation(s)
- Jun-Ho Choi
- Center for Molecular Spectroscopy and Dynamics, Institute for Basic Science, Korea University, Seoul 02841, South Korea and Department of Chemistry, Korea University, Seoul 02841, South Korea
| | - Minhaeng Cho
- Center for Molecular Spectroscopy and Dynamics, Institute for Basic Science, Korea University, Seoul 02841, South Korea and Department of Chemistry, Korea University, Seoul 02841, South Korea
| |
Collapse
|
21
|
Dans PD, Danilāne L, Ivani I, Dršata T, Lankaš F, Hospital A, Walther J, Pujagut RI, Battistini F, Gelpí JL, Lavery R, Orozco M. Long-timescale dynamics of the Drew-Dickerson dodecamer. Nucleic Acids Res 2016; 44:4052-66. [PMID: 27084952 PMCID: PMC4872116 DOI: 10.1093/nar/gkw264] [Citation(s) in RCA: 69] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2015] [Accepted: 03/31/2016] [Indexed: 12/24/2022] Open
Abstract
We present a systematic study of the long-timescale dynamics of the Drew–Dickerson dodecamer (DDD: d(CGCGAATTGCGC)2) a prototypical B-DNA duplex. Using our newly parameterized PARMBSC1 force field, we describe the conformational landscape of DDD in a variety of ionic environments from minimal salt to 2 M Na+Cl− or K+Cl−. The sensitivity of the simulations to the use of different solvent and ion models is analyzed in detail using multi-microsecond simulations. Finally, an extended (10 μs) simulation is used to characterize slow and infrequent conformational changes in DDD, leading to the identification of previously uncharacterized conformational states of this duplex which can explain biologically relevant conformational transitions. With a total of more than 43 μs of unrestrained molecular dynamics simulation, this study is the most extensive investigation of the dynamics of the most prototypical DNA duplex.
Collapse
Affiliation(s)
- Pablo D Dans
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, Baldiri Reixac 10-12, 08028 Barcelona, Spain Joint BSC-IRB Research Program in Computational Biology, Baldiri Reixac 10-12, 08028 Barcelona, Spain
| | - Linda Danilāne
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, Baldiri Reixac 10-12, 08028 Barcelona, Spain Joint BSC-IRB Research Program in Computational Biology, Baldiri Reixac 10-12, 08028 Barcelona, Spain School of Chemistry, University of East Anglia (UEA), Norwich Research Park, Norwich NR4 7TJ, UK
| | - Ivan Ivani
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, Baldiri Reixac 10-12, 08028 Barcelona, Spain Joint BSC-IRB Research Program in Computational Biology, Baldiri Reixac 10-12, 08028 Barcelona, Spain
| | - Tomáš Dršata
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Flemingovo nám 2, 166 10 Prague, Czech Republic
| | - Filip Lankaš
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Flemingovo nám 2, 166 10 Prague, Czech Republic Laboratory of Informatics and Chemistry, University of Chemistry and Technology Prague, Technická 5, 166 28 Prague, Czech Republic
| | - Adam Hospital
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, Baldiri Reixac 10-12, 08028 Barcelona, Spain Joint BSC-IRB Research Program in Computational Biology, Baldiri Reixac 10-12, 08028 Barcelona, Spain
| | - Jürgen Walther
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, Baldiri Reixac 10-12, 08028 Barcelona, Spain Joint BSC-IRB Research Program in Computational Biology, Baldiri Reixac 10-12, 08028 Barcelona, Spain
| | - Ricard Illa Pujagut
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, Baldiri Reixac 10-12, 08028 Barcelona, Spain Joint BSC-IRB Research Program in Computational Biology, Baldiri Reixac 10-12, 08028 Barcelona, Spain
| | - Federica Battistini
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, Baldiri Reixac 10-12, 08028 Barcelona, Spain Joint BSC-IRB Research Program in Computational Biology, Baldiri Reixac 10-12, 08028 Barcelona, Spain
| | - Josep Lluis Gelpí
- Department of Biochemistry and Molecular Biology, University of Barcelona, 08028 Barcelona, Spain
| | - Richard Lavery
- Bases Moléculaires et Structurales des Systèmes Infectieux, Université Lyon I/CNRS UMR 5086, IBCP, 7 Passage du Vercors, Lyon 69367, France
| | - Modesto Orozco
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, Baldiri Reixac 10-12, 08028 Barcelona, Spain Joint BSC-IRB Research Program in Computational Biology, Baldiri Reixac 10-12, 08028 Barcelona, Spain Department of Biochemistry and Molecular Biology, University of Barcelona, 08028 Barcelona, Spain
| |
Collapse
|
22
|
Abstract
Metal ions are essential cofactors for the structure and functions of nucleic acids. Yet, the early discovery in the 70s of the crucial role of Mg(2+) in stabilizing tRNA structures has occulted for a long time the importance of monovalent cations. Renewed interest in these ions was brought in the late 90s by the discovery of specific potassium metal ions in the core of a group I intron. Their importance in nucleic acid folding and catalytic activity is now well established. However, detection of K(+) and Na(+) ions is notoriously problematic and the question about their specificity is recurrent. Here we review the different methods that can be used to detect K(+) and Na(+) ions in nucleic acid structures such as X-ray crystallography, nuclear magnetic resonance or molecular dynamics simulations. We also discuss specific versus non-specific binding to different structures through various examples.
Collapse
Affiliation(s)
- Pascal Auffinger
- Architecture et Réactivité de l'ARN, Université de Strasbourg, IBMC, CNRS, 15 rue René Descartes, F-67084, Strasbourg, France.
| | - Luigi D'Ascenzo
- Architecture et Réactivité de l'ARN, Université de Strasbourg, IBMC, CNRS, 15 rue René Descartes, F-67084, Strasbourg, France.
| | - Eric Ennifar
- Architecture et Réactivité de l'ARN, Université de Strasbourg, IBMC, CNRS, 15 rue René Descartes, F-67084, Strasbourg, France.
| |
Collapse
|
23
|
Cruz FJ, de Pablo JJ, Mota JP. Nanoscopic Characterization of DNA within Hydrophobic Pores: Thermodynamics and Kinetics. Biochem Eng J 2015. [DOI: 10.1016/j.bej.2015.04.027] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
24
|
Havrila M, Zgarbová M, Jurečka P, Banáš P, Krepl M, Otyepka M, Šponer J. Microsecond-Scale MD Simulations of HIV-1 DIS Kissing-Loop Complexes Predict Bulged-In Conformation of the Bulged Bases and Reveal Interesting Differences between Available Variants of the AMBER RNA Force Fields. J Phys Chem B 2015; 119:15176-90. [DOI: 10.1021/acs.jpcb.5b08876] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Marek Havrila
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Královopolská
135, 612 65 Brno, Czech Republic
| | - Marie Zgarbová
- Regional
Centre of Advanced Technologies and Materials, Department of Physical
Chemistry, Faculty of Science, Palacký University, tř.
17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Petr Jurečka
- Regional
Centre of Advanced Technologies and Materials, Department of Physical
Chemistry, Faculty of Science, Palacký University, tř.
17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Pavel Banáš
- Regional
Centre of Advanced Technologies and Materials, Department of Physical
Chemistry, Faculty of Science, Palacký University, tř.
17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Miroslav Krepl
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Královopolská
135, 612 65 Brno, Czech Republic
| | - Michal Otyepka
- Regional
Centre of Advanced Technologies and Materials, Department of Physical
Chemistry, Faculty of Science, Palacký University, tř.
17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Jiří Šponer
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Královopolská
135, 612 65 Brno, Czech Republic
- CEITEC
- Central European Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| |
Collapse
|
25
|
Stadlbauer P, Kührová P, Banáš P, Koča J, Bussi G, Trantírek L, Otyepka M, Šponer J. Hairpins participating in folding of human telomeric sequence quadruplexes studied by standard and T-REMD simulations. Nucleic Acids Res 2015; 43:9626-44. [PMID: 26433223 PMCID: PMC4787745 DOI: 10.1093/nar/gkv994] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2015] [Accepted: 09/22/2015] [Indexed: 11/13/2022] Open
Abstract
DNA G-hairpins are potential key structures participating in folding of human telomeric guanine quadruplexes (GQ). We examined their properties by standard MD simulations starting from the folded state and long T-REMD starting from the unfolded state, accumulating ∼130 μs of atomistic simulations. Antiparallel G-hairpins should spontaneously form in all stages of the folding to support lateral and diagonal loops, with sub-μs scale rearrangements between them. We found no clear predisposition for direct folding into specific GQ topologies with specific syn/anti patterns. Our key prediction stemming from the T-REMD is that an ideal unfolded ensemble of the full GQ sequence populates all 4096 syn/anti combinations of its four G-stretches. The simulations can propose idealized folding pathways but we explain that such few-state pathways may be misleading. In the context of the available experimental data, the simulations strongly suggest that the GQ folding could be best understood by the kinetic partitioning mechanism with a set of deep competing minima on the folding landscape, with only a small fraction of molecules directly folding to the native fold. The landscape should further include non-specific collapse processes where the molecules move via diffusion and consecutive random rare transitions, which could, e.g. structure the propeller loops.
Collapse
Affiliation(s)
- Petr Stadlbauer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic
| | - Petra Kührová
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Jaroslav Koča
- CEITEC - Central European Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic National Center for Biomolecular Research, Faculty of Science, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati, Via Bonomea 265, 34136 Trieste, Italy
| | - Lukáš Trantírek
- CEITEC - Central European Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic CEITEC - Central European Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| |
Collapse
|
26
|
Antila HS, Härkönen M, Sammalkorpi M. Chemistry specificity of DNA-polycation complex salt response: a simulation study of DNA, polylysine and polyethyleneimine. Phys Chem Chem Phys 2015; 17:5279-89. [PMID: 25607687 DOI: 10.1039/c4cp04967e] [Citation(s) in RCA: 32] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
In this work, the chemistry specific stability determining factors of DNA-polycation complexes are examined by performing all-atom molecular dynamics simulations. To this end, we conduct a systematic variation of polycation line charge through polyethyleneimine (PEI) protonation and polycation chemistry via comparison with poly-l-lysine (PLL). Our simulations show that increasing line charge of the polycation alone does not lead to more salt tolerant complexes. Instead, the effective charge compensation by the polycation correlates with the increased stability of the complex against additional salt. The salt stability of PEI-DNA complexes also links to the proton sponge property of weak polycations, commonly assumed to be behind the effectivity of PEI as a gene delivery vector. Examination of the complexes reveals the mechanism behind this behaviour; more Cl(-) ions are attracted by the protonated complexes but, in contrast to the common depiction of the proton sponge behaviour, the ion influx does not cause swelling of the complex structure itself. However, PEI protonation leads to release of PEI while DNA remains tightly bound to the complex. Jointly, these findings shed light on the stability determining factors of DNA-polycation complexes, raise charge distribution as an important stability determining contributor, and indicate that the effectivity of PEI in gene delivery is likely to result from the freed PEI facilitating gene transfection.
Collapse
Affiliation(s)
- Hanne S Antila
- Aalto University School of Chemical Technology, Department of Chemistry, P.O. Box 16100, FI-00076, Aalto, Finland.
| | | | | |
Collapse
|
27
|
Choi JH, Cho M. Ion aggregation in high salt solutions. IV. Graph-theoretical analyses of ion aggregate structure and water hydrogen bonding network. J Chem Phys 2015; 143:104110. [DOI: 10.1063/1.4930608] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jun-Ho Choi
- Department of Chemistry, Korea University, Seoul 136-713, South Korea
| | - Minhaeng Cho
- Department of Chemistry, Korea University, Seoul 136-713, South Korea
- Center for Molecular Spectroscopy and Dynamics, Institute for Basic Science, Korea University, Seoul 136-713, South Korea
| |
Collapse
|
28
|
Islam B, Stadlbauer P, Krepl M, Koca J, Neidle S, Haider S, Sponer J. Extended molecular dynamics of a c-kit promoter quadruplex. Nucleic Acids Res 2015; 43:8673-93. [PMID: 26245347 PMCID: PMC4605300 DOI: 10.1093/nar/gkv785] [Citation(s) in RCA: 42] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2015] [Accepted: 07/21/2015] [Indexed: 01/29/2023] Open
Abstract
The 22-mer c-kit promoter sequence folds into a parallel-stranded quadruplex with a unique structure, which has been elucidated by crystallographic and NMR methods and shows a high degree of structural conservation. We have carried out a series of extended (up to 10 μs long, ∼50 μs in total) molecular dynamics simulations to explore conformational stability and loop dynamics of this quadruplex. Unfolding no-salt simulations are consistent with a multi-pathway model of quadruplex folding and identify the single-nucleotide propeller loops as the most fragile part of the quadruplex. Thus, formation of propeller loops represents a peculiar atomistic aspect of quadruplex folding. Unbiased simulations reveal μs-scale transitions in the loops, which emphasizes the need for extended simulations in studies of quadruplex loops. We identify ion binding in the loops which may contribute to quadruplex stability. The long lateral-propeller loop is internally very stable but extensively fluctuates as a rigid entity. It creates a size-adaptable cleft between the loop and the stem, which can facilitate ligand binding. The stability gain by forming the internal network of GA base pairs and stacks of this loop may be dictating which of the many possible quadruplex topologies is observed in the ground state by this promoter quadruplex.
Collapse
Affiliation(s)
- Barira Islam
- Central European Institute of Technology (CEITEC), Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Petr Stadlbauer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic
| | - Miroslav Krepl
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic
| | - Jaroslav Koca
- Central European Institute of Technology (CEITEC), Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic National Center for Biomolecular Research, Faculty of Science, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Stephen Neidle
- UCL School of Pharmacy, University College London, 29-39 Brunswick Square, London WC1N 1AX, UK
| | - Shozeb Haider
- UCL School of Pharmacy, University College London, 29-39 Brunswick Square, London WC1N 1AX, UK
| | - Jiri Sponer
- Central European Institute of Technology (CEITEC), Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic
| |
Collapse
|
29
|
Choi JH, Kim H, Kim S, Lim S, Chon B, Cho M. Ion aggregation in high salt solutions. III. Computational vibrational spectroscopy of HDO in aqueous salt solutions. J Chem Phys 2015; 142:204102. [DOI: 10.1063/1.4920972] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jun-Ho Choi
- Center for Molecular Spectroscopy and Dynamics, Institute for Basic Science, Seoul 136-713, South Korea
- Department of Chemistry, Korea University, Seoul 136-713, South Korea
| | - Heejae Kim
- Department of Chemistry, Korea University, Seoul 136-713, South Korea
| | - Seongheun Kim
- Department of Chemistry, Korea University, Seoul 136-713, South Korea
| | - Sohee Lim
- Center for Molecular Spectroscopy and Dynamics, Institute for Basic Science, Seoul 136-713, South Korea
- Department of Chemistry, Korea University, Seoul 136-713, South Korea
| | - Bonghwan Chon
- Center for Molecular Spectroscopy and Dynamics, Institute for Basic Science, Seoul 136-713, South Korea
- Department of Chemistry, Korea University, Seoul 136-713, South Korea
| | - Minhaeng Cho
- Center for Molecular Spectroscopy and Dynamics, Institute for Basic Science, Seoul 136-713, South Korea
- Department of Chemistry, Korea University, Seoul 136-713, South Korea
| |
Collapse
|
30
|
Cruz FJAL, de Pablo JJ, Mota JPB. Endohedral confinement of a DNA dodecamer onto pristine carbon nanotubes and the stability of the canonical B form. J Chem Phys 2015; 140:225103. [PMID: 24929415 DOI: 10.1063/1.4881422] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
Although carbon nanotubes are potential candidates for DNA encapsulation and subsequent delivery of biological payloads to living cells, the thermodynamical spontaneity of DNA encapsulation under physiological conditions is still a matter of debate. Using enhanced sampling techniques, we show for the first time that, given a sufficiently large carbon nanotube, the confinement of a double-stranded DNA segment, 5'-D(*CP*GP*CP*GP*AP*AP*TP*TP*CP*GP*CP*G)-3', is thermodynamically favourable under physiological environments (134 mM, 310 K, 1 bar), leading to DNA-nanotube hybrids with lower free energy than the unconfined biomolecule. A diameter threshold of 3 nm is established below which encapsulation is inhibited. The confined DNA segment maintains its translational mobility and exhibits the main geometrical features of the canonical B form. To accommodate itself within the nanopore, the DNA's end-to-end length increases from 3.85 nm up to approximately 4.1 nm, due to a ~0.3 nm elastic expansion of the strand termini. The canonical Watson-Crick H-bond network is essentially conserved throughout encapsulation, showing that the contact between the DNA segment and the hydrophobic carbon walls results in minor rearrangements of the nucleotides H-bonding. The results obtained here are paramount to the usage of carbon nanotubes as encapsulation media for next generation drug delivery technologies.
Collapse
Affiliation(s)
- Fernando J A L Cruz
- Requimte/CQFB, Departamento de Química, Faculdade de Ciências e Tecnologia, Universidade Nova de Lisboa, Caparica 2829-516, Portugal
| | - Juan J de Pablo
- Department of Chemical and Biological Engineering, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| | - José P B Mota
- Requimte/CQFB, Departamento de Química, Faculdade de Ciências e Tecnologia, Universidade Nova de Lisboa, Caparica 2829-516, Portugal
| |
Collapse
|
31
|
Kathuria P, Sharma P, Abendong MN, Wetmore SD. Conformational Preferences of DNA following Damage by Aristolochic Acids: Structural and Energetic Insights into the Different Mutagenic Potential of the ALI and ALII-N6-dA Adducts. Biochemistry 2015; 54:2414-28. [PMID: 25761009 DOI: 10.1021/bi501484m] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Preetleen Kathuria
- Department of Chemistry and
Biochemistry, University of Lethbridge, 4401 University Drive West, Lethbridge, Alberta, Canada T1K 3M4
| | - Purshotam Sharma
- Department of Chemistry and
Biochemistry, University of Lethbridge, 4401 University Drive West, Lethbridge, Alberta, Canada T1K 3M4
| | - Minette N. Abendong
- Department of Chemistry and
Biochemistry, University of Lethbridge, 4401 University Drive West, Lethbridge, Alberta, Canada T1K 3M4
| | - Stacey D. Wetmore
- Department of Chemistry and
Biochemistry, University of Lethbridge, 4401 University Drive West, Lethbridge, Alberta, Canada T1K 3M4
| |
Collapse
|
32
|
Panecka J, Šponer J, Trylska J. Conformational dynamics of bacterial and human cytoplasmic models of the ribosomal A-site. Biochimie 2015; 112:96-110. [PMID: 25748164 DOI: 10.1016/j.biochi.2015.02.021] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2014] [Accepted: 02/23/2015] [Indexed: 01/12/2023]
Abstract
The aminoacyl-tRNA binding site (A-site) is located in helix 44 of small ribosomal subunit. The mobile adenines 1492 and 1493 (Escherichia coli numbering), forming the A-site bulge, act as a functional switch that ensures mRNA decoding accuracy. Structural data on the oligonucleotide models mimicking the ribosomal A-site with sequences corresponding to bacterial and human cytoplasmic sites confirm that this RNA motif forms also without the ribosome context. We performed all-atom molecular dynamics simulations of these crystallographic A-site models to compare their conformational properties. We found that the human A-site bulge is more internally flexible than the bacterial one and has different base pairing preferences, which result in the overall different shapes of these bulges and cation density distributions. Also, in the human A-site model we observed repetitive destacking of A1492, while A1493 was more stably paired than in the bacterial variant. Based on the dynamics of the A-sites we suggest why aminoglycoside antibiotics, which target the bacterial A-site, have lower binding affinities and anti-translational activities toward the human variant.
Collapse
Affiliation(s)
- Joanna Panecka
- Division of Biophysics, Institute of Experimental Physics, University of Warsaw, Żwirki i Wigury 93, 02-089 Warsaw, Poland; Interdisciplinary Centre for Mathematical and Computational Modelling, University of Warsaw, Pawińskiego 5a, 02-106 Warsaw, Poland
| | - Jiří Šponer
- CEITEC - Central European Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic; Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic.
| | - Joanna Trylska
- Centre of New Technologies, University of Warsaw, Banacha 2c, 02-097 Warsaw, Poland.
| |
Collapse
|
33
|
Choi JH, Cho M. Ion aggregation in high salt solutions. II. Spectral graph analysis of water hydrogen-bonding network and ion aggregate structures. J Chem Phys 2014; 141:154502. [DOI: 10.1063/1.4897638] [Citation(s) in RCA: 42] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023] Open
Affiliation(s)
- Jun-Ho Choi
- Department of Chemistry, Korea University, Seoul 136-713, South Korea
| | - Minhaeng Cho
- Department of Chemistry, Korea University, Seoul 136-713, South Korea
| |
Collapse
|
34
|
Kim S, Kim H, Choi JH, Cho M. Ion aggregation in high salt solutions: Ion network versus ion cluster. J Chem Phys 2014; 141:124510. [DOI: 10.1063/1.4896227] [Citation(s) in RCA: 46] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Seongheun Kim
- Department of Chemistry, Korea University, Seoul 136-713, South Korea
| | - Heejae Kim
- Department of Chemistry, Korea University, Seoul 136-713, South Korea
| | - Jun-Ho Choi
- Department of Chemistry, Korea University, Seoul 136-713, South Korea
| | - Minhaeng Cho
- Department of Chemistry, Korea University, Seoul 136-713, South Korea
| |
Collapse
|
35
|
Giambaşu GM, Luchko T, Herschlag D, York DM, Case DA. Ion counting from explicit-solvent simulations and 3D-RISM. Biophys J 2014; 106:883-94. [PMID: 24559991 PMCID: PMC3944826 DOI: 10.1016/j.bpj.2014.01.021] [Citation(s) in RCA: 83] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2013] [Revised: 12/19/2013] [Accepted: 01/10/2014] [Indexed: 12/30/2022] Open
Abstract
The ionic atmosphere around nucleic acids remains only partially understood at atomic-level detail. Ion counting (IC) experiments provide a quantitative measure of the ionic atmosphere around nucleic acids and, as such, are a natural route for testing quantitative theoretical approaches. In this article, we replicate IC experiments involving duplex DNA in NaCl(aq) using molecular dynamics (MD) simulation, the three-dimensional reference interaction site model (3D-RISM), and nonlinear Poisson-Boltzmann (NLPB) calculations and test against recent buffer-equilibration atomic emission spectroscopy measurements. Further, we outline the statistical mechanical basis for interpreting IC experiments and clarify the use of specific concentration scales. Near physiological concentrations, MD simulation and 3D-RISM estimates are close to experimental results, but at higher concentrations (>0.7 M), both methods underestimate the number of condensed cations and overestimate the number of excluded anions. The effect of DNA charge on ion and water atmosphere extends 20-25 Å from its surface, yielding layered density profiles. Overall, ion distributions from 3D-RISMs are relatively close to those from corresponding MD simulations, but with less Na(+) binding in grooves and tighter binding to phosphates. NLPB calculations, on the other hand, systematically underestimate the number of condensed cations at almost all concentrations and yield nearly structureless ion distributions that are qualitatively distinct from those generated by both MD simulation and 3D-RISM. These results suggest that MD simulation and 3D-RISM may be further developed to provide quantitative insight into the characterization of the ion atmosphere around nucleic acids and their effect on structure and stability.
Collapse
Affiliation(s)
- George M Giambaşu
- Department of Chemistry and Chemical Biology and BioMaPS Institute, Rutgers University, Piscataway, New Jersey
| | - Tyler Luchko
- Department of Chemistry and Chemical Biology and BioMaPS Institute, Rutgers University, Piscataway, New Jersey
| | - Daniel Herschlag
- Department of Biochemistry, Stanford University, Stanford, California
| | - Darrin M York
- Department of Chemistry and Chemical Biology and BioMaPS Institute, Rutgers University, Piscataway, New Jersey.
| | - David A Case
- Department of Chemistry and Chemical Biology and BioMaPS Institute, Rutgers University, Piscataway, New Jersey.
| |
Collapse
|
36
|
Kesselheim S, Müller W, Holm C. Origin of current blockades in nanopore translocation experiments. PHYSICAL REVIEW LETTERS 2014; 112:018101. [PMID: 24483933 DOI: 10.1103/physrevlett.112.018101] [Citation(s) in RCA: 45] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/07/2013] [Indexed: 06/03/2023]
Abstract
We present a detailed investigation of the ionic current in a cylindrical model nanopore in the absence and the presence of a double stranded DNA homopolymer. Our atomistic simulations are capable of reproducing almost exactly the experimental data obtained by Smeets et al., including notably the crossover salt concentration that yields equal current measurements in both situations. We can rule out that the observed current blockade is due to the steric exclusion of charge carriers from the DNA, since for all investigated salt concentrations the charge carrier density is higher when the DNA is present. Calculations using a mean-field electrokinetic model proposed by van Dorp et al. fail quantitatively in predicting this effect. We can relate the shortcomings of the mean-field model to a surface related molecular drag that the ions feel in the presence of the DNA. This drag is independent of the salt concentration and originates from electrostatic, hydrodynamic, and excluded volume interactions.
Collapse
Affiliation(s)
- Stefan Kesselheim
- Institut für Computerphysik, Universität Stuttgart, Allmandring 3, D-70569 Stuttgart, Germany
| | - Wojciech Müller
- Institut für Computerphysik, Universität Stuttgart, Allmandring 3, D-70569 Stuttgart, Germany
| | - Christian Holm
- Institut für Computerphysik, Universität Stuttgart, Allmandring 3, D-70569 Stuttgart, Germany
| |
Collapse
|
37
|
Cruz FJAL, de Pablo JJ, Mota JPB. Free energy landscapes of the encapsulation mechanism of DNA nucleobases onto carbon nanotubes. RSC Adv 2014. [DOI: 10.1039/c3ra45033c] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023] Open
|
38
|
Noy A, Golestanian R. Noy and Golestanian reply:. PHYSICAL REVIEW LETTERS 2013; 111:179802. [PMID: 24206522 DOI: 10.1103/physrevlett.111.179802] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/03/2013] [Indexed: 06/02/2023]
Affiliation(s)
- Agnes Noy
- Rudolf Peierls Centre for Theoretical Physics, University of Oxford, 1 Keble Road, Oxford OX1 3NP, United Kingdom
| | | |
Collapse
|
39
|
Do TN, Carloni P, Varani G, Bussi G. RNA/Peptide Binding Driven by Electrostatics-Insight from Bidirectional Pulling Simulations. J Chem Theory Comput 2013; 9:1720-30. [PMID: 26587630 DOI: 10.1021/ct3009914] [Citation(s) in RCA: 41] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
RNA/protein interactions play crucial roles in controlling gene expression. They are becoming important targets for pharmaceutical applications. Due to RNA flexibility and to the strength of electrostatic interactions, standard docking methods are insufficient. We here present a computational method which allows studying the binding of RNA molecules and charged peptides with atomistic, explicit-solvent molecular dynamics. In our method, a suitable estimate of the electrostatic interaction is used as an order parameter (collective variable) which is then accelerated using bidirectional pulling simulations. Since the electrostatic interaction is only used to enhance the sampling, the approximations used to compute it do not affect the final accuracy. The method is employed to characterize the binding of TAR RNA from HIV-1 and a small cyclic peptide. Our simulation protocol allows blindly predicting the binding pocket and pose as well as the binding affinity. The method is general and could be applied to study other electrostatics-driven binding events.
Collapse
Affiliation(s)
- Trang N Do
- SISSA/ISAS - International School for Advanced Studies, Trieste 34136, Italy
| | - Paolo Carloni
- Computational Biophysics, German Research School for Simulation Sciences, D-52425 Jülich, Germany and Institute for Advanced Simulation IAS-5, Computational Biomedicine, Forschungszentrum Jülich, D-52425 Jülich, Germany
| | - Gabriele Varani
- Department of Chemistry and Department of Biochemistry, University of Washington, Seattle, Washington 98195-1700, United States
| | - Giovanni Bussi
- SISSA/ISAS - International School for Advanced Studies, Trieste 34136, Italy
| |
Collapse
|
40
|
Foloppe N, Guéroult M, Hartmann B. Simulating DNA by molecular dynamics: aims, methods, and validation. Methods Mol Biol 2013; 924:445-468. [PMID: 23034759 DOI: 10.1007/978-1-62703-017-5_17] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/01/2023]
Abstract
The structure and dynamics of the B-DNA double helix involves subtle sequence-dependent effects which are decisive for its function, but difficult to characterize. These structural and dynamic effects can be addressed by simulations of DNA sequences in explicit solvent. Here, we present and discuss the state-of-art of B-DNA molecular dynamics simulations with the major force fields in use today. We explain why a critical analysis of the MD trajectories is required to assess their reliability, and estimate the value and limitations of these models. Overall, simulations of DNA bear great promise towards deciphering the structural and physical subtleties of this biopolymer, where much remains to be understood.
Collapse
|
41
|
Robbins TJ, Wang Y. Effect of initial ion positions on the interactions of monovalent and divalent ions with a DNA duplex as revealed with atomistic molecular dynamics simulations. J Biomol Struct Dyn 2012; 31:1311-23. [PMID: 23153112 DOI: 10.1080/07391102.2012.732344] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
Monovalent (Na(+)) and divalent (Mg(2+)) ion distributions around the Dickerson-Drew dodecamer were studied by atomistic molecular dynamics (MD) simulations with AMBER molecular modeling software. Different initial placements of ions were tried and the resulting effects on the ion distributions around DNA were investigated. For monovalent ions, results were found to be nearly independent of initial cation coordinates. However, Mg(2+) ions demonstrated a strong initial coordinate dependent behavior. While some divalent ions initially placed near the DNA formed essentially permanent direct coordination complexes with electronegative DNA atoms, Mg(2+) ions initially placed further away from the duplex formed a full, nonexchanging, octahedral first solvation shell. These fully solvated cations were still capable of binding with DNA with events lasting up to 20 ns, and in comparison were bound much longer than Na(+) ions. Force field parameters were also investigated with modest and little differences arising from ion (ions94 and ions08) and nucleic acid description (ff99, ff99bsc0, and ff10), respectively. Based on known Mg(2+) ion solvation structure, we conclude that in most cases Mg(2+) ions retain their first solvation shell, making only solvent-mediated contacts with DNA duplex. The proper way to simulate Mg(2+) ions around DNA duplex, therefore, should begin with ions placed in the bulk water.
Collapse
Affiliation(s)
- Timothy J Robbins
- a Department of Chemistry , University of Memphis , Memphis , TN , 38152 , USA
| | | |
Collapse
|
42
|
Dršata T, Pérez A, Orozco M, Morozov AV, Sponer J, Lankaš F. Structure, Stiffness and Substates of the Dickerson-Drew Dodecamer. J Chem Theory Comput 2012; 9:707-721. [PMID: 23976886 DOI: 10.1021/ct300671y] [Citation(s) in RCA: 68] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
The Dickerson-Drew dodecamer (DD) d-[CGCGAATTCGCG]2 is a prototypic B-DNA molecule whose sequence-specific structure and dynamics have been investigated by many experimental and computational studies. Here, we present an analysis of DD properties based on extensive atomistic molecular dynamics (MD) simulations using different ionic conditions and water models. The 0.6-2.4-µs-long MD trajectories are compared to modern crystallographic and NMR data. In the simulations, the duplex ends can adopt an alternative base-pairing, which influences the oligomer structure. A clear relationship between the BI/BII backbone substates and the basepair step conformation has been identified, extending previous findings and exposing an interesting structural polymorphism in the helix. For a given end pairing, distributions of the basepair step coordinates can be decomposed into Gaussian-like components associated with the BI/BII backbone states. The nonlocal stiffness matrices for a rigid-base mechanical model of DD are reported for the first time, suggesting salient stiffness features of the central A-tract. The Riemann distance and Kullback-Leibler divergence are used for stiffness matrix comparison. The basic structural parameters converge very well within 300 ns, convergence of the BI/BII populations and stiffness matrices is less sharp. Our work presents new findings about the DD structural dynamics, mechanical properties, and the coupling between basepair and backbone configurations, including their statistical reliability. The results may also be useful for optimizing future force fields for DNA.
Collapse
Affiliation(s)
- Tomáš Dršata
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Flemingovo nám. 2, 166 10 Prague, Czech Republic ; Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, Brno, Czech Republic
| | | | | | | | | | | |
Collapse
|
43
|
Beššeová I, Banáš P, Kührová P, Košinová P, Otyepka M, Šponer J. Simulations of A-RNA Duplexes. The Effect of Sequence, Solute Force Field, Water Model, and Salt Concentration. J Phys Chem B 2012; 116:9899-916. [DOI: 10.1021/jp3014817] [Citation(s) in RCA: 59] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Ivana Beššeová
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska
135, 612 65 Brno, Czech Republic
| | - Pavel Banáš
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska
135, 612 65 Brno, Czech Republic
- Regional Centre of Advanced
Technologies and Materials, Department of Physical Chemistry, Faculty
of Science, Palacky University, tr. 17
listopadu 12, 771 46, Olomouc, Czech Republic
| | - Petra Kührová
- Regional Centre of Advanced
Technologies and Materials, Department of Physical Chemistry, Faculty
of Science, Palacky University, tr. 17
listopadu 12, 771 46, Olomouc, Czech Republic
| | - Pavlína Košinová
- Regional Centre of Advanced
Technologies and Materials, Department of Physical Chemistry, Faculty
of Science, Palacky University, tr. 17
listopadu 12, 771 46, Olomouc, Czech Republic
| | - Michal Otyepka
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska
135, 612 65 Brno, Czech Republic
- Regional Centre of Advanced
Technologies and Materials, Department of Physical Chemistry, Faculty
of Science, Palacky University, tr. 17
listopadu 12, 771 46, Olomouc, Czech Republic
| | - Jiří Šponer
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska
135, 612 65 Brno, Czech Republic
- CEITEC - Central European Institute of Technology, Campus Bohunice, Kamenice
5, 625 00 Brno, Czech Republic
| |
Collapse
|
44
|
Spasic A, Serafini J, Mathews DH. The Amber ff99 Force Field Predicts Relative Free Energy Changes for RNA Helix Formation. J Chem Theory Comput 2012; 8:2497-2505. [PMID: 23112748 PMCID: PMC3482406 DOI: 10.1021/ct300240k] [Citation(s) in RCA: 33] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
The ability of the Amber ff99 force field to predict relative free energies of RNA helix formation was investigated. The test systems were three hexaloop RNA hairpins with identical loops and varying stems. The potential of mean force of stretching the hairpins from the native state to an extended conformation was calculated with umbrella sampling. Because the hairpins have identical loop sequence, the differences in free energy changes are only from the stem composition. The Amber ff99 force field was able to correctly predict the order of stabilities of the hairpins, although the magnitude of the free energy change is larger than that determined by optical melting experiments. The two measurements cannot be compared directly because the unfolded state in the optical melting experiments is a random coil, while the end state in the umbrella sampling simulations was an elongated chain. The calculations can be compared to reference data by using a thermodynamic cycle. By applying the thermodynamic cycle to the transitions between the hairpins using simulations and nearest neighbor data, agreement was found to be within the sampling error of simulations, thus demonstrating that ff99 force field is able to accurately predict relative free energies of RNA helix formation.
Collapse
Affiliation(s)
- Aleksandar Spasic
- Department of Biochemistry & Biophysics and Center for RNA Biology, University of Rochester Medical Center, Rochester, New York
| | - John Serafini
- Department of Biochemistry & Biophysics and Center for RNA Biology, University of Rochester Medical Center, Rochester, New York
| | - David H. Mathews
- Department of Biochemistry & Biophysics and Center for RNA Biology, University of Rochester Medical Center, Rochester, New York
- Department of Biostatistics & Computational Biology, University of Rochester Medical Center, Rochester, New York
| |
Collapse
|
45
|
Arcella A, Portella G, Ruiz ML, Eritja R, Vilaseca M, Gabelica V, Orozco M. Structure of Triplex DNA in the Gas Phase. J Am Chem Soc 2012; 134:6596-606. [DOI: 10.1021/ja209786t] [Citation(s) in RCA: 47] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Affiliation(s)
- Annalisa Arcella
- Joint IRB BSC Research Program
in Computational Biology, Institute for Research in Biomedicine, Baldiri Reixach 10, Barcelona 08028, Spain
| | - Guillem Portella
- Joint IRB BSC Research Program
in Computational Biology, Institute for Research in Biomedicine, Baldiri Reixach 10, Barcelona 08028, Spain
| | - Maria Luz Ruiz
- Chemistry and Molecular Pharmacology
Program, Institute for Research in Biomedicine, IQAC-CSIC, CIBER-BBN, Barcelona 08028, Spain
| | - Ramon Eritja
- Chemistry and Molecular Pharmacology
Program, Institute for Research in Biomedicine, IQAC-CSIC, CIBER-BBN, Barcelona 08028, Spain
| | - Marta Vilaseca
- Mass Spectrometry Core Facility, Institute for Research in Biomedicine, Barcelona 08028,
Spain
| | - Valérie Gabelica
- Department of Chemistry, University of Liège, Allée de la Chimie,
Building B6c, B-4000 Liège, Belgium
| | - Modesto Orozco
- Joint IRB BSC Research Program
in Computational Biology, Institute for Research in Biomedicine, Baldiri Reixach 10, Barcelona 08028, Spain
- Departament de Bioquímica
i Biología Molecular, Facultat de Biología, Universitat de Barcelona, Avgda Diagonal 645, Barcelona
08028, Spain
| |
Collapse
|
46
|
Cerón-Carrasco JP, Requena A, Jacquemin D. Impact of DFT functionals on the predicted magnesium–DNA interaction: an ONIOM study. Theor Chem Acc 2012. [DOI: 10.1007/s00214-012-1188-9] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
|
47
|
Abstract
It has been known for decades that DNA is extremely flexible and polymorphic, but our knowledge of its accessible conformational space remains limited. Structural data, primarily from X-ray diffraction studies, is sparse in comparison to the manifold configurations possible, and direct experimental examinations of DNA's flexibility still suffer from many limitations. In the face of these shortcomings, molecular dynamics (MD) is now an essential tool in the study of DNA. It affords detailed structural and dynamical insights, which explains its recent transition from a small number of highly specialized laboratories to a large variety of groups dealing with challenging biological problems. MD is now making an irreversible journey to the mainstream of research in biology, with the attendant opportunities and challenges. But given the speed with which MD studies of DNA have spread, the roots remain somewhat shallow: in many cases, there is a lack of deep knowledge about the foundations, strengths, and limits of the technique. In this Account, we discuss how MD has become the most important source of structural and flexibility data on DNA, focusing on advances since 2007 of atomistic MD in the description of DNA under near-physiological conditions and highlighting the possibilities and shortcomings of the technique. The evolution in the field over the past four years is a prelude to the ongoing revolution. The technique has gained in robustness and predictive power, which when coupled with the spectacular improvements in software and hardware has enabled the tackling of systems of increasing complexity. Simulation times of microseconds have now been achieved, with even longer times when specialized hardware is used. As a result, we have seen the first real-time simulation of large conformational transitions, including folding and unfolding of short DNA duplexes. Noteworthy advances have also been made in the study of DNA-ligand interactions, and we predict that a global thermodynamic and kinetic picture of the binding landscape of DNA will become available in a few years. MD will become a crucial tool in areas such as biomolecular engineering and synthetic biology. MD has also been shown to be an excellent source of parameters for mesoscopic models of DNA flexibility. Such models can be refined through atomistic MD simulations on small duplexes and then applied to the study of entire chromosomes. Recent evidence suggests that MD-derived elastic models can successfully predict the position of regulatory regions in DNA and can help advance our understanding of nucleosome positioning and chromatin plasticity. If these results are confirmed, MD simulations can become the ultimate tool to decipher a physical code that can contribute to gene regulation. We are entering the golden age of MD simulations of DNA. Undoubtedly, the expectations are high, but the challenges are also enormous. These include the need for more accurate potential energy functionals and for longer and more complex simulations in more realistic systems. The joint research effort of several groups will be crucial for adapting the technique to the requirements of the coming decade.
Collapse
Affiliation(s)
- Alberto Pérez
- Joint IRB-BSC Program in Computational Biology, Institute of Research in Biomedicine Barcelona, Baldiri i Reixac 10, Barcelona 08028, Spain
| | - F. Javier Luque
- Department de Fisicoquímica and Institut de Biomedicina (IBUB), Facultat de Farmàcia, Universitat de Barcelona, Avgda Diagonal 643, Barcelona 08028, Spain
| | - Modesto Orozco
- Joint IRB-BSC Program in Computational Biology, Institute of Research in Biomedicine Barcelona, Baldiri i Reixac 10, Barcelona 08028, Spain
- Departament de Bioquímica, Universitat de Barcelona, Avgda Diagonal 647, Barcelona 08028, Spain, and Instituto Nacional de Bioinformàtica, Parc Científic de Barcelona, Baldiri i Reixac 10, Barcelona 08028, Spain
| |
Collapse
|
48
|
|
49
|
Kučerka N, Holland BW, Gray CG, Tomberli B, Katsaras J. Scattering density profile model of POPG bilayers as determined by molecular dynamics simulations and small-angle neutron and X-ray scattering experiments. J Phys Chem B 2011; 116:232-9. [PMID: 22107350 DOI: 10.1021/jp208920h] [Citation(s) in RCA: 84] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
We combine molecular dynamics (MD) simulations and experiment, both small-angle neutron (SANS) and small-angle X-ray scattering (SAXS), to determine the precise structure of bilayers composed of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylglycerol (POPG), a lipid commonly encountered in bacterial membranes. Experiment and simulation are used to develop a one-dimensional scattering density profile (SDP) model suitable for the analysis of experimental data. The joint refinement of such data (i.e., SANS and SAXS) results in the area per lipid that is then used in the fixed-area simulations. In the final step, the direct comparison of simulated-to-experimental data gives rise to the detailed structure of POPG bilayers. From these studies we conclude that POPG's molecular area is 66.0 ± 1.3 Å(2), its overall bilayer thickness is 36.7 ± 0.7 Å, and its hydrocarbon region thickness is 27.9 ± 0.6 Å, assuming a simulated value of 1203 Å(3) for the total lipid volume.
Collapse
Affiliation(s)
- Norbert Kučerka
- Canadian Neutron Beam Centre, National Research Council, Chalk River, Ontario, Canada.
| | | | | | | | | |
Collapse
|
50
|
Réblová K, Šponer JE, Špačková N, Beššeová I, Šponer J. A-minor tertiary interactions in RNA kink-turns. Molecular dynamics and quantum chemical analysis. J Phys Chem B 2011; 115:13897-910. [PMID: 21999672 DOI: 10.1021/jp2065584] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
The RNA kink-turn is an important recurrent RNA motif, an internal loop with characteristic consensus sequence forming highly conserved three-dimensional structure. Functional arrangement of RNA kink-turns shows a sharp bend in the phosphodiester backbone. Among other signature interactions, kink-turns form A-minor interaction between their two stems. Most kink-turns possess extended A-minor I (A-I) interaction where adenine of the second A•G base pair of the NC-stem interacts with the first canonical pair of the C-stem (i.e., the receptor pair) via trans-sugar-edge/sugar-edge (tSS) and cis-sugar-edge/sugar-edge (cSS) interactions. The remaining kink-turns have less compact A-minor 0 (A-0) interaction with just one tSS contact. We show that kink-turns with A-I in ribosomal X-ray structures keep G═C receptor base pair during evolution while the inverted pair (C═G) is not realized. In contrast, kink-turns with A-0 in the observed structures alternate G═C and C═G base pairs in sequences. We carried out an extended set (~5 μs) of explicit-solvent molecular dynamics simulations of kink-turns to rationalize this structural/evolutionary pattern. The simulations were done using a net-neutral Na(+) cation atmosphere (with ~0.25 M cation concentration) supplemented by simulations with either excess salt KCl atmosphere or inclusion of Mg(2+). The results do not seem to depend on the treatment of ions. The simulations started with X-ray structures of several kink-turns while we tested the response of the simulated system to base substitutions, modest structural perturbations and constraints. The trends seen in the simulations reveal that the A-I/G═C arrangement is preferred over all three other structures. The A-I/C═G triple appears structurally entirely unstable, consistent with the covariation patterns seen during the evolution. The A-0 arrangements tend to shift toward the A-I pattern in simulations, which suggests that formation of the A-0 interaction is likely supported by the surrounding protein and RNA molecules. A-0 may also be stabilized by additional kink-turn nucleotides not belonging to the kink-turn consensus, as shown for the kink-turn from ribosomal Helix 15. Quantum-chemical calculations on all four A-minor triples suggest that there is a different balance of electrostatic and dispersion stabilization in the A-I/G═C and A-I/C═G triples, which may explain different behavior of these otherwise isosteric triples in the context of kink-turns.
Collapse
Affiliation(s)
- Kamila Réblová
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 61265 Brno, Czech Republic.
| | | | | | | | | |
Collapse
|