1
|
Knappeová B, Mlýnský V, Pykal M, Šponer J, Banáš P, Otyepka M, Krepl M. Comprehensive Assessment of Force-Field Performance in Molecular Dynamics Simulations of DNA/RNA Hybrid Duplexes. J Chem Theory Comput 2024. [PMID: 39012172 DOI: 10.1021/acs.jctc.4c00601] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/17/2024]
Abstract
Mixed double helices formed by RNA and DNA strands, commonly referred to as hybrid duplexes or hybrids, are essential in biological processes like transcription and reverse transcription. They are also important for their applications in CRISPR gene editing and nanotechnology. Yet, despite their significance, the hybrid duplexes have been seldom modeled by atomistic molecular dynamics methodology, and there is no benchmark study systematically assessing the force-field performance. Here, we present an extensive benchmark study of polypurine tract (PPT) and Dickerson-Drew dodecamer hybrid duplexes using contemporary and commonly utilized pairwise additive and polarizable nucleic acid force fields. Our findings indicate that none of the available force-field choices accurately reproduces all the characteristic structural details of the hybrid duplexes. The AMBER force fields are unable to populate the C3'-endo (north) pucker of the DNA strand and underestimate inclination. The CHARMM force field accurately describes the C3'-endo pucker and inclination but shows base pair instability. The polarizable force fields struggle with accurately reproducing the helical parameters. Some force-field combinations even demonstrate a discernible conflict between the RNA and DNA parameters. In this work, we offer a candid assessment of the force-field performance for mixed DNA/RNA duplexes. We provide guidance on selecting utilizable force-field combinations and also highlight potential pitfalls and best practices for obtaining optimal performance.
Collapse
Affiliation(s)
- Barbora Knappeová
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, Brno 612 00, Czech Republic
| | - Vojtěch Mlýnský
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, Brno 612 00, Czech Republic
| | - Martin Pykal
- Czech Advanced Technology and Research Institute, CATRIN, Palacký University, Křížkovského 511/8, Olomouc 779 00, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, Brno 612 00, Czech Republic
| | - Pavel Banáš
- Czech Advanced Technology and Research Institute, CATRIN, Palacký University, Křížkovského 511/8, Olomouc 779 00, Czech Republic
| | - Michal Otyepka
- Czech Advanced Technology and Research Institute, CATRIN, Palacký University, Křížkovského 511/8, Olomouc 779 00, Czech Republic
- IT4Innovations, VSB-Technical University of Ostrava, 17. listopadu 2172/15, Ostrava-Poruba 708 00, Czech Republic
| | - Miroslav Krepl
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, Brno 612 00, Czech Republic
| |
Collapse
|
2
|
Pirojsirikul T, Lee VS, Nimmanpipug P. Unraveling Bacterial Single-Stranded Sequence Specificities: Insights from Molecular Dynamics and MMPBSA Analysis of Oligonucleotide Probes. Mol Biotechnol 2024; 66:582-591. [PMID: 38374320 DOI: 10.1007/s12033-024-01082-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2022] [Accepted: 01/10/2024] [Indexed: 02/21/2024]
Abstract
We utilized molecular dynamics (MD) simulations and Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) free energy calculations to investigate the specificity of two oligonucleotide probes, namely probe B and probe D, in detecting single-stranded DNA (ssDNA) within three bacteria families: Enterobacteriaceae, Pasteurellaceae, and Vibrionaceae. Due to the limited understanding of molecular mechanisms in the previous research, we have extended the discussion to focus specifically on investigating the binding process of bacteria-probe DNA duplexes, with an emphasis on analyzing the binding free energy. The role of electrostatic contributions in the specificity between the oligonucleotide probes and the bacterial ssDNAs was investigated and found to be crucial. Our calculations yielded results that were highly consistent with the experimental data. Through our study, we have successfully exhibited the benefits of utilizing in-silico approaches as a powerful virtual-screening tool, particularly in research areas that demand a thorough comprehension of molecular interactions.
Collapse
Affiliation(s)
- Teerapong Pirojsirikul
- Division of Physical Science, Faculty of Science, Prince of Songkla University, Songkhla, 90110, Thailand.
| | - Vannajan Sanghiran Lee
- Department of Chemistry, Center of Theoretical and Computational Physics, Faculty of Science, University of Malaya, 50603, Kuala Lumpur, Malaysia
| | - Piyarat Nimmanpipug
- Department of Chemistry, Faculty of Science, Chiang Mai University, Chiang Mai, 50200, Thailand
| |
Collapse
|
3
|
Yushin II, Golyshev VM, Pyshnyi DV, Lomzov AA. Application of the weighted histogram method for calculating the thermodynamic parameters of the formation of oligodeoxyribonucleotide duplexes. Vavilovskii Zhurnal Genet Selektsii 2023; 27:807-814. [PMID: 38213713 PMCID: PMC10777287 DOI: 10.18699/vjgb-23-93] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2023] [Revised: 09/20/2023] [Accepted: 09/21/2023] [Indexed: 01/13/2024] Open
Abstract
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 ΔG° 37). 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 ΔG° 37 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 ΔG° 37, 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.
Collapse
Affiliation(s)
- I I Yushin
- Institute of Chemical Biology and Fundamental Medicine of the Siberian Branch of the Russian Academy of Sciences, Novosibirsk, Russia Novosibirsk State University, Novosibirsk, Russia
| | - V M Golyshev
- Institute of Chemical Biology and Fundamental Medicine of the Siberian Branch of the Russian Academy of Sciences, Novosibirsk, Russia Novosibirsk State University, Novosibirsk, Russia
| | - D V Pyshnyi
- Institute of Chemical Biology and Fundamental Medicine of the Siberian Branch of the Russian Academy of Sciences, Novosibirsk, Russia
| | - A A Lomzov
- Institute of Chemical Biology and Fundamental Medicine of the Siberian Branch of the Russian Academy of Sciences, Novosibirsk, Russia Novosibirsk State University, Novosibirsk, Russia
| |
Collapse
|
4
|
Oxenfarth A, Kümmerer F, Bottaro S, Schnieders R, Pinter G, Jonker HRA, Fürtig B, Richter C, Blackledge M, Lindorff-Larsen K, Schwalbe H. Integrated NMR/Molecular Dynamics Determination of the Ensemble Conformation of a Thermodynamically Stable CUUG RNA Tetraloop. J Am Chem Soc 2023. [PMID: 37479220 PMCID: PMC10401711 DOI: 10.1021/jacs.3c03578] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/23/2023]
Abstract
Both experimental and theoretical structure determinations of RNAs have remained challenging due to the intrinsic dynamics of RNAs. We report here an integrated nuclear magnetic resonance/molecular dynamics (NMR/MD) structure determination approach to describe the dynamic structure of the CUUG tetraloop. We show that the tetraloop undergoes substantial dynamics, leading to averaging of the experimental data. These dynamics are particularly linked to the temperature-dependent presence of a hydrogen bond within the tetraloop. Interpreting the NMR data by a single structure represents the low-temperature structure well but fails to capture all conformational states occurring at a higher temperature. We integrate MD simulations, starting from structures of CUUG tetraloops within the Protein Data Bank, with an extensive set of NMR data, and provide a structural ensemble that describes the dynamic nature of the tetraloop and the experimental NMR data well. We thus show that one of the most stable and frequently found RNA tetraloops displays substantial dynamics, warranting such an integrated structural approach.
Collapse
Affiliation(s)
- Andreas Oxenfarth
- Institute for Organic Chemistry and Chemical Biology, Center for Biomolecular Magnetic Resonance (BMRZ), Goethe University Frankfurt am Main, Max-von-Laue-Str. 7, 60438 Frankfurt/Main, Hessen, Germany
| | - Felix Kümmerer
- Structural Biology and NMR Laboratory, Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200 Copenhagen N, Denmark
| | - Sandro Bottaro
- Structural Biology and NMR Laboratory, Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200 Copenhagen N, Denmark
- IRCCS Humanitas Research Hospital, Department of Biomedical Sciences, Humanitas University, Milan 20089, Italy
| | - Robbin Schnieders
- Institute for Organic Chemistry and Chemical Biology, Center for Biomolecular Magnetic Resonance (BMRZ), Goethe University Frankfurt am Main, Max-von-Laue-Str. 7, 60438 Frankfurt/Main, Hessen, Germany
| | - György Pinter
- Institute for Organic Chemistry and Chemical Biology, Center for Biomolecular Magnetic Resonance (BMRZ), Goethe University Frankfurt am Main, Max-von-Laue-Str. 7, 60438 Frankfurt/Main, Hessen, Germany
| | - Hendrik R A Jonker
- Institute for Organic Chemistry and Chemical Biology, Center for Biomolecular Magnetic Resonance (BMRZ), Goethe University Frankfurt am Main, Max-von-Laue-Str. 7, 60438 Frankfurt/Main, Hessen, Germany
| | - Boris Fürtig
- Institute for Organic Chemistry and Chemical Biology, Center for Biomolecular Magnetic Resonance (BMRZ), Goethe University Frankfurt am Main, Max-von-Laue-Str. 7, 60438 Frankfurt/Main, Hessen, Germany
| | - Christian Richter
- Institute for Organic Chemistry and Chemical Biology, Center for Biomolecular Magnetic Resonance (BMRZ), Goethe University Frankfurt am Main, Max-von-Laue-Str. 7, 60438 Frankfurt/Main, Hessen, Germany
| | - Martin Blackledge
- Institut de Biologie Structurale (IBS), CEA, CNRS, University Grenoble Alpes, Grenoble 38000, France
| | - Kresten Lindorff-Larsen
- Structural Biology and NMR Laboratory, Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200 Copenhagen N, Denmark
| | - Harald Schwalbe
- Institute for Organic Chemistry and Chemical Biology, Center for Biomolecular Magnetic Resonance (BMRZ), Goethe University Frankfurt am Main, Max-von-Laue-Str. 7, 60438 Frankfurt/Main, Hessen, Germany
| |
Collapse
|
5
|
Dowerah D, V. N. Uppuladinne M, Sarma PJ, Biswakarma N, Sonavane UB, Joshi RR, Ray SK, Namsa ND, Deka RC. Design of LNA Analogues Using a Combined Density Functional Theory and Molecular Dynamics Approach for RNA Therapeutics. ACS OMEGA 2023; 8:22382-22405. [PMID: 37396274 PMCID: PMC10308574 DOI: 10.1021/acsomega.2c07860] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/09/2022] [Accepted: 02/17/2023] [Indexed: 07/04/2023]
Abstract
Antisense therapeutics treat a wide spectrum of diseases, many of which cannot be addressed with the current drug technologies. In the quest to design better antisense oligonucleotide drugs, we propose five novel LNA analogues (A1-A5) for modifying antisense oligonucleotides and establishing each with the five standard nucleic acids: adenine (A), guanine (G), cytosine (C), thymine (T), and uracil (U). Monomer nucleotides of these modifications were considered for a detailed Density Functional Theory (DFT)-based quantum chemical analysis to determine their molecular-level structural and electronic properties. A detailed MD simulation study was done on a 14-mer ASO (5'-CTTAGCACTGGCCT-3') containing these modifications targeting PTEN mRNA. Results from both molecular- and oligomer-level analysis clearly depicted LNA-level stability of the modifications, the ASO/RNA duplexes maintaining stable Watson-Crick base pairing preferring RNA-mimicking A-form duplexes. Notably, monomer MO isosurfaces for both purines and pyrimidines were majorly distributed on the nucleobase region in modifications A1 and A2 and in the bridging unit in modifications A3, A4, and A5, suggesting that A3/RNA, A4/RNA, and A5/RNA duplexes interact more with the RNase H and solvent environment. Accordingly, solvation of A3/RNA, A4/RNA, and A5/RNA duplexes was higher compared to that of LNA/RNA, A1/RNA, and A2/RNA duplexes. This study has resulted in a successful archetype for creating advantageous nucleic acid modifications tailored for particular needs, fulfilling a useful purpose of designing novel antisense modifications, which may overcome the drawbacks and improve the pharmacokinetics of existing LNA antisense modifications.
Collapse
Affiliation(s)
- Dikshita Dowerah
- CMML—Catalysis
and Molecular Modelling Lab, Department of Chemical Sciences, Tezpur University, Napaam, Sonitpur, Assam 784 028, India
| | - Mallikarjunachari V. N. Uppuladinne
- HPC—Medical
& Bioinformatics Applications Group, Centre for Development of Advanced Computing (C-DAC), Panchavati, Pashan, Pune 411008, India
| | - Plaban J. Sarma
- CMML—Catalysis
and Molecular Modelling Lab, Department of Chemical Sciences, Tezpur University, Napaam, Sonitpur, Assam 784 028, India
- Department
of Chemistry, Gargaon College, Sivasagar, Assam 785685, India
| | - Nishant Biswakarma
- CMML—Catalysis
and Molecular Modelling Lab, Department of Chemical Sciences, Tezpur University, Napaam, Sonitpur, Assam 784 028, India
| | - Uddhavesh B. Sonavane
- HPC—Medical
& Bioinformatics Applications Group, Centre for Development of Advanced Computing (C-DAC), Panchavati, Pashan, Pune 411008, India
| | - Rajendra R. Joshi
- HPC—Medical
& Bioinformatics Applications Group, Centre for Development of Advanced Computing (C-DAC), Panchavati, Pashan, Pune 411008, India
| | - Suvendra K. Ray
- Department
of Molecular Biology and Biotechnology, Tezpur University, Napaam, Sonitpur, Assam 784028, India
- Center
for Multidisciplinary Research, Tezpur University, Napaam, Sonitpur, Assam 784028, India
| | - Nima D. Namsa
- Department
of Molecular Biology and Biotechnology, Tezpur University, Napaam, Sonitpur, Assam 784028, India
- Center
for Multidisciplinary Research, Tezpur University, Napaam, Sonitpur, Assam 784028, India
| | - Ramesh Ch. Deka
- CMML—Catalysis
and Molecular Modelling Lab, Department of Chemical Sciences, Tezpur University, Napaam, Sonitpur, Assam 784 028, India
- Center
for Multidisciplinary Research, Tezpur University, Napaam, Sonitpur, Assam 784028, India
| |
Collapse
|
6
|
Goodman J, Attwood D, Kiely J, Coladas Mato P, Luxton R. Modeling Peptide Nucleic Acid Binding Enthalpies Using MM-GBSA. J Phys Chem B 2022; 126:9528-9538. [PMID: 36375178 PMCID: PMC9706560 DOI: 10.1021/acs.jpcb.2c05547] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
Abstract
The binding enthalpies of peptide nucleic acid (PNA) homoduplexes were predicted using a molecular mechanics generalized Born surface area approach. Using the nucleic acid nearest-neighbor model, these were decomposed into sequence parameters which could replicate the enthalpies from thermal melting experiments with a mean error of 8.7%. These results present the first systematic computational investigation into the relationship between sequence and binding energy for PNA homoduplexes and identified a stabilizing helix initiation enthalpy not observed for nucleic acids with phosphoribose backbones.
Collapse
Affiliation(s)
- Jack Goodman
- University
of the West of England, BristolBS16 1QY, U.K.,
| | - David Attwood
- University
of the West of England, BristolBS16 1QY, U.K.
| | - Janice Kiely
- University
of the West of England, BristolBS16 1QY, U.K.
| | | | - Richard Luxton
- University
of the West of England, BristolBS16 1QY, U.K.
| |
Collapse
|
7
|
Baranova SV, Zhdanova PV, Lomzov AA, Koval VV, Chernonosov AA. Structure- and Content-Dependent Efficiency of Cas9-Assisted DNA Cleavage in Genome-Editing Systems. Int J Mol Sci 2022; 23:ijms232213889. [PMID: 36430368 PMCID: PMC9693425 DOI: 10.3390/ijms232213889] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2022] [Revised: 11/02/2022] [Accepted: 11/09/2022] [Indexed: 11/16/2022] Open
Abstract
Genome-editing systems, being some of the key tools of molecular biologists, represent a reasonable hope for progress in the field of personalized medicine. A major problem with such systems is their nonideal accuracy and insufficient selectivity. The selectivity of CRISPR-Cas9 systems can be improved in several ways. One efficient way is the proper selection of the consensus sequence of the DNA to be cleaved. In the present work, we attempted to evaluate the effect of formed non-Watson-Crick pairs in a DNA duplex on the efficiency of DNA cleavage in terms of the influence of the structure of the formed partially complementary pairs. We also studied the effect of the location of such pairs in DNA relative to the PAM (protospacer-adjacent motif) on the cleavage efficiency. We believe that the stabilization of the Cas9-sgRNA complex with a DNA substrate containing noncomplementary pairs is due to loop reorganization in the RuvC domain of the enzyme. In addition, PAM-proximal mismatches in the DNA substrate lower enzyme efficiency because the "seed" region is involved in binding and cleavage, whereas PAM-distal mismatches have no significant impact on target DNA cleavage. Our data suggest that in the case of short duplexes with mismatches, the stages of recognition and binding of dsDNA substrates by the enzyme determine the reaction rate and time rather than the thermodynamic parameters affected by the "unwinding" of DNA. The results will provide a theoretical basis for predicting the efficiency and accuracy of CRISPR-Cas9 systems at cleaving target DNA.
Collapse
Affiliation(s)
- Svetlana V. Baranova
- Institute of Chemical Biology and Fundamental Medicine, Siberian Branch of Russian Academy of Sciences (ICBFM SB RAS), 630090 Novosibirsk, Russia
| | - Polina V. Zhdanova
- Institute of Chemical Biology and Fundamental Medicine, Siberian Branch of Russian Academy of Sciences (ICBFM SB RAS), 630090 Novosibirsk, Russia
| | - Alexander A. Lomzov
- Institute of Chemical Biology and Fundamental Medicine, Siberian Branch of Russian Academy of Sciences (ICBFM SB RAS), 630090 Novosibirsk, Russia
- Department of Natural Sciences, Novosibirsk State University, 630090 Novosibirsk, Russia
| | - Vladimir V. Koval
- Institute of Chemical Biology and Fundamental Medicine, Siberian Branch of Russian Academy of Sciences (ICBFM SB RAS), 630090 Novosibirsk, Russia
- Department of Natural Sciences, Novosibirsk State University, 630090 Novosibirsk, Russia
- Correspondence:
| | - Alexander A. Chernonosov
- Institute of Chemical Biology and Fundamental Medicine, Siberian Branch of Russian Academy of Sciences (ICBFM SB RAS), 630090 Novosibirsk, Russia
| |
Collapse
|
8
|
Zhdanova PV, Lomzov AA, Prokhorova DV, Stepanov GA, Chernonosov AA, Koval VV. Thermodynamic Swings: How Ideal Complex of Cas9-RNA/DNA Forms. Int J Mol Sci 2022; 23:8891. [PMID: 36012157 PMCID: PMC9408429 DOI: 10.3390/ijms23168891] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2022] [Revised: 08/02/2022] [Accepted: 08/09/2022] [Indexed: 12/26/2022] Open
Abstract
Most processes of the recognition and formation of specific complexes in living systems begin with collisions in solutions or quasi-solutions. Then, the thermodynamic regulation of complex formation and fine tuning of complexes come into play. Precise regulation is very important in all cellular processes, including genome editing using the CRISPR-Cas9 tool. The Cas9 endonuclease is an essential component of the CRISPR-Cas-based genome editing systems. The attainment of high-specificity and -efficiency Cas9 during targeted DNA cleavage is the main problem that limits the practical application of the CRISPR-Cas9 system. In this study, we analyzed the thermodynamics of interaction of a complex's components of Cas9-RNA/DNA through experimental and computer simulation methods. We found that there is a small energetic preference during Cas9-RNA/DNA formation from the Cas9-RNA and DNA/DNA duplex. The small difference in binding energy is relevant for biological interactions and could be part of the sequence-specific recognition of double-stranded DNA by the CRISPR-Cas9 system.
Collapse
Affiliation(s)
- Polina V. Zhdanova
- Institute of Chemical Biology and Fundamental Medicine, Siberian Branch of the Russian Academy of Sciences (SB RAS), 630090 Novosibirsk, Russia
- Department of Natural Sciences, Novosibirsk State University, 630090 Novosibirsk, Russia
| | - Alexander A. Lomzov
- Institute of Chemical Biology and Fundamental Medicine, Siberian Branch of the Russian Academy of Sciences (SB RAS), 630090 Novosibirsk, Russia
- Department of Natural Sciences, Novosibirsk State University, 630090 Novosibirsk, Russia
| | - Daria V. Prokhorova
- Institute of Chemical Biology and Fundamental Medicine, Siberian Branch of the Russian Academy of Sciences (SB RAS), 630090 Novosibirsk, Russia
| | - Grigory A. Stepanov
- Institute of Chemical Biology and Fundamental Medicine, Siberian Branch of the Russian Academy of Sciences (SB RAS), 630090 Novosibirsk, Russia
| | - Alexander A. Chernonosov
- Institute of Chemical Biology and Fundamental Medicine, Siberian Branch of the Russian Academy of Sciences (SB RAS), 630090 Novosibirsk, Russia
| | - Vladimir V. Koval
- Institute of Chemical Biology and Fundamental Medicine, Siberian Branch of the Russian Academy of Sciences (SB RAS), 630090 Novosibirsk, Russia
- Department of Natural Sciences, Novosibirsk State University, 630090 Novosibirsk, Russia
| |
Collapse
|