As of today, influenza viruses remain a relevant target for the development of antiviral compounds due to their rapid evolution and acquisition of the resistance to existing drugs. Fullerene derivatives have already shown the ability to successfully interact with viruses, and polyhydrated fullerenes (or fullerenols) are particularly attractive due to their compatibility with biological fluids and low toxicity. Therefore, the goal of this work was to study the effect of two batches of a mixture of polyhydrated fullerenes with a mass ratio of 78.1% C60/C70 and 21.9% C76/C78/C84 on the influenza A (H1N1) virus. It was determined that the mixture of fullerenols, along with the low toxicity, showed high antiviral activity with a decrease in the viral infectious titer up to 4 orders of magnitude. In addition, studied fullerenols did not affect the hemagglutination process and did not show any significant prophylactic activity. With the help of molecular docking and molecular dynamics simulation, the likely target of fullerenols' action was determined—the binding site of the RNA primer of the viral RNA-dependent RNA polymerase. Therefore, we assume that the high antiviral effect of polyhydrated fullerenes on influenza A virus is related to their interaction with the viral RNA polymerase.

Cancer as an acquired genetic disease is based on changes both in the genome itself and in transcription processes. Accordingly, it is at the DNA level that it makes sense to search for and design agents capable of effective and selective anticancer action. In this study, we used an iterative approach based on a molecular dynamics simulation to design a highly selective DNA-intercalating agent called HASDI. To confirm its selective affinity to DNA, we conducted two simulation experiments: HASDI in a complex with a DNA fragment of the EBNA1 gene (it targets 16 nucleotide pairs of this gene) and HASDI in a complex with a random DNA fragment of the KCNH2 gene. The molecular dynamics simulation was carried out in the GROMACS 2019 package. The binding energy was calculated by gmx_MMPBSA 1.5.2. The further analysis was performed using the built-in utilities of GROMACS, gmx_MMPBSA and also XMGRACE and Pymol 1.8. As a result, we determined that the EBNA1-50nt/HASDI complex was stable throughout the whole simulation trajectory. HASDI, due to the presence of a linker modified depending on a specific pair of nitrogenous bases, formed an average of 32 hydrogen bonds with a sequence of 16 nucleotide pairs. Phenazine rings were stably intercalated every 2 base pairs. The root-mean-square deviation of HASDI in such a complex fluctuated around the value of 6.5 Å and had no tendency to increase. The calculated value of the binding free energy was − 235.3 ± 7.77 kcal/mol. The KCNH2-50nt/HASDI complex, as an example of the intercalation of the designed structure into a random part of the human genome, maintained the stability of its position at a level comparable to the EBNA1-50nt/HASDI complex. The phenazine rings were constantly intercalated in their original positions, and the root-mean-square deviation fluctuated around one value, although it had a tendency to chaotic changes. At the same time, this complex was characterized by 17–19 hydrogen bonds, on average, and the binding free energy was − 193.47 ± 14.09 kcal/mol. Moreover, the DNA duplex had local single-nucleotide melting in the region of the 4th linker. According to a significant decrease in the number of hydrogen bonds, a decrease in energy gain, as well as a decrease in the stability of the DNA duplex characteristic of the KCNH2-50nt/HASDI complex compared to the target EBNA1-50nt/HASDI complex, the molecule we designed can be considered a potentially selective DNA polyintercalating agent capable of relatively accurate recognition of 16 base pairs.

The possibility of precise external control of transcription opens up many prospects. One of them is dealing with disorders based on enhanced transcription of a certain part of the genome or its mutant fragment. Therefore, the aim of our study is to prove the fundamental possibility of developing DNA-intercalators capable of highly specific binding to certain pre-selected areas of DNA by using a number of in silico methods. A positively charged bis-intercalator-XR5944 (PDB ID: 2MG8) was chosen as the basis. The target 50-nucleotide DNA duplex (a part of the viral protein EBNA1 gene) was generated along with preliminary intercalation of XR5944 in Avogadro program. Then by manual modification and data analysis, structure of the DNA-intercalator was created. Thus, we developed a structural model of the DNA-intercalator which during a molecular dynamic experiment, along with constant intercalation, selectively recognized the nucleotide sequence of the DNA duplex. This process is based on the specific placement of the hydrogen bond donors and acceptors directed to each pair of nitrogenous bases in the direction of the large DNA furrow. So, we have proved the possibility of developing substances-DNA-intercalators-that can bind highly specifically to certain pre-selected areas of DNA.

No abstract

No abstract