Application of the weighted histogram method for calculating the thermodynamic parameters of the formation of oligodeoxyribonucleotide duplexes

To date, many derivatives and analogs of nucleic acids (NAs) have been developed. Some of them have found uses in scientific research and biomedical applications. Their effective use is based on the data about their properties. Some of the most important physicochemical properties of oligonucleotide...

Full description

Saved in:
Bibliographic Details
Main Authors: I. I. Yushin, V. M. Golyshev, D. V. Pyshnyi, A. A. Lomzov
Format: Article
Language:English
Published: Siberian Branch of the Russian Academy of Sciences, Federal Research Center Institute of Cytology and Genetics, The Vavilov Society of Geneticists and Breeders 2023-12-01
Series:Вавиловский журнал генетики и селекции
Subjects:
Online Access:https://vavilov.elpub.ru/jour/article/view/3980
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1832575021743079424
author I. I. Yushin
V. M. Golyshev
D. V. Pyshnyi
A. A. Lomzov
author_facet I. I. Yushin
V. M. Golyshev
D. V. Pyshnyi
A. A. Lomzov
author_sort I. I. Yushin
collection DOAJ
description To date, many derivatives and analogs of nucleic acids (NAs) have been developed. Some of them have found uses in scientific research and biomedical applications. Their effective use is based on the data about their properties. Some of the most important physicochemical properties of oligonucleotides are thermodynamic parameters of the formation of their duplexes with DNA and RNA. These parameters can be calculated only for a few NA derivatives: locked NAs, bridged oligonucleotides, and peptide NAs. Existing predictive approaches are based on an analysis of experimental data and the consequent construction of predictive models. The ongoing pilot studies aimed at devising methods for predicting the properties of NAs by computational modeling techniques are based only on knowledge about the structure of oligonucleotides. In this work, we studied the applicability of the weighted histogram analysis method (WHAM) in combination with umbrella sampling to the calculation of thermodynamic parameters of DNA duplex formation (changes in enthalpy ∆H°, entropy ∆S°, and Gibbs free energy          ∆G37° ). A procedure was designed involving WHAM for calculating the hybridization properties of oligodeoxyribonucleotides. Optimal parameters for modeling and calculation of thermodynamic parameters were determined. The feasibility of calculation of ∆H°, ∆S°, and          ∆G37° was demonstrated using a representative sample of 21 oligonucleotides 4–16 nucleotides long with a GC content of 14–100 %. Error of the calculation of the thermodynamic parameters was 11.4, 12.9, and 11.8 % for ∆H°, ∆S°, and          ∆G37° , respectively, and the melting temperature was predicted with an average error of 5.5 °C. Such high accuracy of computations is comparable with the accuracy of the experimental approach and of other methods for calculating the energy of NA duplex formation. In this paper, the use of WHAM for computation of the energy of DNA duplex formation was systematically investigated for the first time. Our results show that a reliable calculation of the hybridization parameters of new NA derivatives is possible, including derivatives not yet synthesized. This work opens up new horizons for a rational design of constructs based on NAs for solving problems in biomedicine and biotechnology.
format Article
id doaj-art-6d2ca25d2f9b458e9dc3d32bd7a58bca
institution Kabale University
issn 2500-3259
language English
publishDate 2023-12-01
publisher Siberian Branch of the Russian Academy of Sciences, Federal Research Center Institute of Cytology and Genetics, The Vavilov Society of Geneticists and Breeders
record_format Article
series Вавиловский журнал генетики и селекции
spelling doaj-art-6d2ca25d2f9b458e9dc3d32bd7a58bca2025-02-01T09:58:12ZengSiberian Branch of the Russian Academy of Sciences, Federal Research Center Institute of Cytology and Genetics, The Vavilov Society of Geneticists and BreedersВавиловский журнал генетики и селекции2500-32592023-12-0127710.18699/VJGB-23-931408Application of the weighted histogram method for calculating the thermodynamic parameters of the formation of oligodeoxyribonucleotide duplexesI. I. Yushin0V. M. Golyshev1D. V. Pyshnyi2A. A. Lomzov3Institute of Chemical Biology and Fundamental Medicine of the Siberian Branch of the Russian Academy of Sciences; Novosibirsk State UniversityInstitute of Chemical Biology and Fundamental Medicine of the Siberian Branch of the Russian Academy of Sciences; Novosibirsk State UniversityInstitute of Chemical Biology and Fundamental Medicine of the Siberian Branch of the Russian Academy of SciencesInstitute of Chemical Biology and Fundamental Medicine of the Siberian Branch of the Russian Academy of Sciences; Novosibirsk State UniversityTo date, many derivatives and analogs of nucleic acids (NAs) have been developed. Some of them have found uses in scientific research and biomedical applications. Their effective use is based on the data about their properties. Some of the most important physicochemical properties of oligonucleotides are thermodynamic parameters of the formation of their duplexes with DNA and RNA. These parameters can be calculated only for a few NA derivatives: locked NAs, bridged oligonucleotides, and peptide NAs. Existing predictive approaches are based on an analysis of experimental data and the consequent construction of predictive models. The ongoing pilot studies aimed at devising methods for predicting the properties of NAs by computational modeling techniques are based only on knowledge about the structure of oligonucleotides. In this work, we studied the applicability of the weighted histogram analysis method (WHAM) in combination with umbrella sampling to the calculation of thermodynamic parameters of DNA duplex formation (changes in enthalpy ∆H°, entropy ∆S°, and Gibbs free energy          ∆G37° ). A procedure was designed involving WHAM for calculating the hybridization properties of oligodeoxyribonucleotides. Optimal parameters for modeling and calculation of thermodynamic parameters were determined. The feasibility of calculation of ∆H°, ∆S°, and          ∆G37° was demonstrated using a representative sample of 21 oligonucleotides 4–16 nucleotides long with a GC content of 14–100 %. Error of the calculation of the thermodynamic parameters was 11.4, 12.9, and 11.8 % for ∆H°, ∆S°, and          ∆G37° , respectively, and the melting temperature was predicted with an average error of 5.5 °C. Such high accuracy of computations is comparable with the accuracy of the experimental approach and of other methods for calculating the energy of NA duplex formation. In this paper, the use of WHAM for computation of the energy of DNA duplex formation was systematically investigated for the first time. Our results show that a reliable calculation of the hybridization parameters of new NA derivatives is possible, including derivatives not yet synthesized. This work opens up new horizons for a rational design of constructs based on NAs for solving problems in biomedicine and biotechnology.https://vavilov.elpub.ru/jour/article/view/3980dnahybridizationthermodynamic parametersgibbs free energyweighted histogram analysis methodwhammolecular dynamics
spellingShingle I. I. Yushin
V. M. Golyshev
D. V. Pyshnyi
A. A. Lomzov
Application of the weighted histogram method for calculating the thermodynamic parameters of the formation of oligodeoxyribonucleotide duplexes
Вавиловский журнал генетики и селекции
dna
hybridization
thermodynamic parameters
gibbs free energy
weighted histogram analysis method
wham
molecular dynamics
title Application of the weighted histogram method for calculating the thermodynamic parameters of the formation of oligodeoxyribonucleotide duplexes
title_full Application of the weighted histogram method for calculating the thermodynamic parameters of the formation of oligodeoxyribonucleotide duplexes
title_fullStr Application of the weighted histogram method for calculating the thermodynamic parameters of the formation of oligodeoxyribonucleotide duplexes
title_full_unstemmed Application of the weighted histogram method for calculating the thermodynamic parameters of the formation of oligodeoxyribonucleotide duplexes
title_short Application of the weighted histogram method for calculating the thermodynamic parameters of the formation of oligodeoxyribonucleotide duplexes
title_sort application of the weighted histogram method for calculating the thermodynamic parameters of the formation of oligodeoxyribonucleotide duplexes
topic dna
hybridization
thermodynamic parameters
gibbs free energy
weighted histogram analysis method
wham
molecular dynamics
url https://vavilov.elpub.ru/jour/article/view/3980
work_keys_str_mv AT iiyushin applicationoftheweightedhistogrammethodforcalculatingthethermodynamicparametersoftheformationofoligodeoxyribonucleotideduplexes
AT vmgolyshev applicationoftheweightedhistogrammethodforcalculatingthethermodynamicparametersoftheformationofoligodeoxyribonucleotideduplexes
AT dvpyshnyi applicationoftheweightedhistogrammethodforcalculatingthethermodynamicparametersoftheformationofoligodeoxyribonucleotideduplexes
AT aalomzov applicationoftheweightedhistogrammethodforcalculatingthethermodynamicparametersoftheformationofoligodeoxyribonucleotideduplexes