PECULIARITY OF ELECTRON DENSITY CALCULATION DURING INTERACTION OF ULTRASHORT LASER PULSE WITH NITROGENOUS BASE OF DNA MOLECULE ADENINE Текст научной статьи по специальности «Физика»

i i St. Petersburg Polytechnic University Journal: Physics and Mathematics. 2022 Vol. 15, No. 3.2 Научно-технические ведомости СПбГПУ. Физико-математические науки. 15 (3.2) 2022

BIOPHYSICS AND MEDICAL PHYSICS

Conference materials UDC 539.192

DOI: https://doi.org/10.18721/JPM.153.253

Peculiarity of electron density calculation during interaction of ultrashort laser pulse with nitrogenous base of DNA molecule adenine

A. A. Kharlamova lH, D. N. Makarov 1

1 Northern (Arctic) Federal University named after M.V. Lomonosov, Arkhangelsk, Russia H kharlamova.anastasya2015@yandex.ru

Abstract. Today, ultrashort laser pulses are used to study and analyze biomolecules. Their advantages over X-rays lie in the speed of measurement and accuracy, which allows studying the structure of a molecule on a time scale commensurate with the time of various biochemical processes in it. However, the practical possibility of studying chemical and biological processes in facilities using existing and promising X-ray free-electron lasers requires the development of a theory of the interaction of laser pulses with matter. In this article we will consider the interaction of an ultrashort laser pulse (USPs) with the base of the DNA (Deoxyribonucleic acid) molecule - adenine. In more detail we will consider the calculation of the photon birth probability, after the interaction of the pulse with the molecule. We will build a model of the base of the molecule and propose a method for calculating the electron density for the case when the substance has a different elemental composition.

Keywords: Deoxyribonucleic acid, ultrashort laser pulse, photon, electron density, elemental composition

Funding: The research was supported by State Task of the Russian Federation FSRU-2021-0008.

Citation: Kharlamova A. A., Makarov D. N., Peculiarity of electron density calculation during interaction of ultrashort laser pulse with nitrogenous base of DNA molecule adenine, St. Petersburg State Polytechnical University Journal. Physics and Mathematics. 15 (3.2) (2022) 291-295. DOI: https://doi.org/10.18721/JPM.153.253

This is an open access article under the CC BY-NC 4.0 license (https://creativecommons. org/licenses/by-nc/4.0/)

Материалы конференции УДК 539.192

DOI: https://doi.org/10.18721/JPM.153.253

Особенность расчета электронной плотности при взаимодействии ультракороткого лазерного импульса с азотистым основанием молекулы ДНК аденином

А. А. Харламова 1Н, Д. Н. Макаров 1

1 Северный (Арктический) федеральный университет им. М.В. Ломоносова, г. Архангельск, Россия

н kharlamova.anastasya2015@yandex.ru

Аннотация. Сегодня применяются для изучения и анализа биомолекул ультракороткие лазерные импульсы. Их преимущества по сравнению с рентгеновским излучением заключается в быстроте измерения и точности. Однако, практическая возможность изучения химических и биологических процессов требует разработки теории взаимодействия лазерных импульсов с веществом. В данной статье мы рассмотрим взаимодействие ультракороткого лазерного импульса (УКИ) с основанием молекулы ДНК - аденином. Подробнее рассмотрим расчет вероятности рождения фотона, предложим способ расчета электронной плотности для случая, когда вещество имеет разный элементный состав.

© Kharlamova A. A., Makarov D. N., 2022. Published by Peter the Great St. Petersburg Polytechnic University.

^ St. Petersburg Polytechnic University Journal. Physics and Mathematics. 2022 Vol. 15, No. 3.2

Ключевые слова: Дезоксирибонуклеиновая кислота, ультракороткий лазерный импульс, фотон, электронная плотность, элементный состав

Финансирование: Исследование выполнено при поддержке государственного задания Российской Федерации № FSRU-2020-0008.

Ссылка при цитировании: Харламова А. А., Макаров Д. Н. Особенность расчета электронной плотности при взаимодействии ультракороткого лазерного импульса с азотистым основанием молекулы ДНК аденином // Научно-технические ведомости СПбГПУ. Физико-математические науки. Т. 15. № 3.2. С. 291-295. DOI: https://doi. org/10.18721/JPM.153.253

Статья открытого доступа, распространяемая по лицензии CC BY-NC 4.0 (https:// creativecommons.org/licenses/by-nc/4.0/)

Introduction

Imaging techniques from diffraction data obtained using ultrashort laser pulses are being actively developed for applications in nanotechnology and structural biology. Methods of X-ray diffraction analysis of matter (XRD) as well as X-ray crystallography (XRC) are well studied and underlie many modern discoveries in physics, chemistry, biology, medicine [1,2]. The structure of a molecule is determined from the time-domain infinite plane wave scattering spectrum [1,2]. But recently scientists have started to apply ultrashort laser pulses for investigation. In this case the scattering spectrum moves into the femto and atto-second time range, which is poorly understood up to now.

The use of (USPs) to analyze the structure of molecules is a very promising area of research. As noted, ultrashort laser pulses (USPs) operate on a small time scale. This is very convenient for the description of organic structures, as such materials can be both stationary and non-stationary. Non-stationary objects, such as DNA, for example, can only be studied using (USPs), because the characteristic time in such systems is smaller or comparable to the femtosecond time scale. The non-stationarity in them may be caused by the formation and rearrangement of chemical bonds, intra-atomic and molecular processes, etc. [3].

Experimental studies of molecules using (USPs) on European X-Ray Free Electron Laser (XFEL) are already in progress [3]. The technique of introducing the molecules into the X-ray beam (e.g. XFEL) is being actively researched and improved. In order to analyse the experimentally obtained data, theoretical descriptions of the interaction processes between (USPs) and biostructures are urgently needed.

So far, a lot of theoretical studies on the specificity of scattered (USPs) on various multi-atomic systems have been carried out. However, despite this fact the theoretical base is still weak and does not take into account all specific features of diffuse (USPs).

In this paper we will consider the DNA molecule as a promising and poorly studied biomolecule at the atomic level. The molecule stores biological information in the form of a genetic code consisting of a sequence of nucleotides. Consequently, a defect on a nucleotide can indicate the presence of a defect in the molecule. In [4] we considered the interaction of an ultrashort laser pulse with a DNA molecule in a simplified form, considering each nucleotide as a separate unit of matter, within the whole molecule. We now focus in detail on a single nucleotide with the nitrogenous base adenine. The difficulty of the new calculation is that now it is necessary to take into account the elemental composition of the molecule. It is possible to take it into account through the electron density of the molecule, which we will consider in this work.

In the following, the atomic system of units will be used, where h = 1, |e| = 1, me = 1.

Materials and Methods

A nucleotide is a compound of a phosphate group, a carbohydrate, and a nitrogenous base — adenine. Let us imagine it as an arbitrary polyatomic system, on which the (USPs) falls in the direction n0. The duration of such a pulse will be considered much less than the characteristic atomic time та = 1, i.e., т/та = 1. This will allow us to use the sudden perturbation approximation.

© Харламова А. А., Макаров Д. Н., 2022. Издатель: Санкт-Петербургский политехнический университет Петра Великого.

In this approximation, the proper Hamiltonian of the system can be neglected, since the electron in the atom does not have time to evolve under the action of the atom's field due to the too fast interaction of the momentum with the electron in the atom [5]. It should be added that the condition x/xa = 1 for using our approximation is not strict. To do this, it is sufficient to assume, as was shown in [5], that, ra0xa = 1, where ra0 is the carrier frequency of the incident (USPs). Obviously, the condition ra0xa = 1 holds for X-ray frequencies. The intensity of the electromagnetic field of the (USPs) E(r,i) is considered in general terms, i.e., we will consider it to depend on the coordinates and time. We are interested in electromagnetic fields that are not strong enough to account for relativistic effects. As was obtained in [6], taking into account these approximations, the probability of the birth of a photon W of a given frequency ra per unit solid angle Qk with the simultaneous transition of the polyatomic system from the ground state to all possible final states (hereinafter, the scattering spectrum of the USPs)

d 2W

1 1

dradQ (2n) c ra

(NaNbG(&, n, no) + 8n (p)(Nb -1)F(ra, n, n0),

(1)

where Na is the total number of atoms in the system, Nb is the number of electrons in the atom, n is the speed of light, G(ra,n,n0) is the direction of escape of the born photon, /(ra,n,n0) and the mean values expressed in terms of electron density p(r) as

If I i2 G(ra,n,n0) = — lp(r) f (r) dr,

N J

F (ra, n, n o) = N- |Jp(r )f (r )e-'kr dr| \

(2) (3)

where k = nra / c, f(ra) = (ra) x n

E (ra) = j

É(ra,e , t) - 1 Va

É(ra,e, t)

eimdt,

where V = .

a,e /dra,e

Furthermore, 5N(p) in (2) is a factor which completely determines the geometrical arrangement of the atoms in the target, which is calculated in general form

SN (p) = Ze

ÍP(Ra -Rb ) _

I

ipRa

(4)

2

Fig. 1. 3D model of the nucleotide: C corresponds to carbohydrate, N to nitrogen, O to oxygen, P to phosphorus, H to hydrogen

[source: Avogadro software]

St. Petersburg Polytechnic University Journal. Physics and Mathematics. 2022 Vol. 15, No. 3.2

Fig. 2. Arrangement of atoms in a nucleotide: C corresponds to carbohydrate, N to nitrogen, O to oxygen, P to phosphorus, H to hydrogen

[source: Avogadro software]

The summation in (3) is carried out for all atoms of the system in question. The electron density can be taken in different models. We will use the Dirac—Hartree—Fock—Slater model from [7]:

a; r

p(r) = £14 a,V

where A,, a. are the table coefficients giving the electron density in the atom [7]. In our case we cannot use this form of electron density since the base of the molecule consists of different elements. Let us present a 3D model of the nucleotide (Fig. 1.). Using the given distances of 1.372 A between the elements and their coordinates in space we will find the electron density.

In general, phosphate consists of 5 atoms bonded together by covalent bonds. The covalent radius of oxygen is 73 pm, the radius of the atomic nucleus is 3.27 fm, for phosphorus respectively 106 pm and 3.27 fm. The radius of the atomic nucleus is calculated as , where A is the number of nucleons. We can assume that the nuclei of atoms are at a distance many times larger than their own radius.

Results and Discussion

Let us set in Avogadro the coordinates of all elements. Then the electron density of the phosphate group will take the form.

P(rph) = Pi 0i) +P2 (r2) +P3(r3) + P4(r4) + P(r), (5)

where is the electron of the relevant element

Pi(ri) = a P V<a>T- (6)

We will use the following data for further calculations (Table 1).

Then the expression G\(ro,n,n0) for each atom, taking formula (4) into account, will take the form:

if ii2

Gt (q n, no) = — |P,. (r) f (r) dr, (7)

Ne

For the whole phosphate group, let us do the summation:

5

Gph (Q n n o) = Z G ^ n o X (8)

i=i

Using this model it is possible to find expressions G.(ra,n,n0) for and go on to further calculate the interaction of USPs with the molecule.

Table 1

Distribution of atoms at the base of adenines by coordinates

Number Element X(A) Y(A) Z(A)

Atom 1 H -1.547 9.377 -1.216

Atom 2 0 -1.44 8.896 -0.401

Atom 3 P 0.038 8.907 0.116

Atom 4 0 0.220 10.175 0.734

Atom 5 0 0.790 8.733 -1.239

Atom 34 H 2.147 0.165 2.718

Atom 35 H 1.815 6.312 5.385

Conclusion

The paper considers the probability of photon emission of a given frequency during the interaction of USPs with the nitrogenous base of DNA molecule, adenine. The peculiarities of the calculation of complex molecule with different elemental chemical composition were pointed out. A method of quantitative finding of electron density to refine the probability calculation is proposed.

Acknowledgments

The research was supported by the state task of the Russian Federation No. FSRU-2021-0008.

REFERENCES

1. Suryanarayana C., Norton M., X-Ray Diffraction: A Practical Approach. Springer Science & Business Media. ISBN 9781489901484.

2. Poletto F.C., Optical Technologies for Extreme-Ultraviolet and Soft X-ray Coherent Sources, Springer-Verlag Berlin Heidelberg, 2015.

3. Zewail A., Voyage through Time: Walks of Life to the Nobel Prize. — World Scientific, 2002.

4. Makarov D.N., Kharlamova A.A., Scattering of ultrashort electromagnetic field pulses by biomolecules, AIP Conference Proceedings 2313, 030070, 2020.

5. Goshev A.A., Eseev M.K., Makarov D.N., Scattering of Attosecond Electromagnetic Field Pulses by Molecular Anions Including the Magnetic Field Component, J. Exp. Theor. Phys. 130 (1), 28—34, 2020.

6. Eseev M.K., Goshev A.A., Makarov D.N., Scattering of Ultrashort X-ray Pulses by Various Nanosystems, Nanomaterials, 10(7), 1355, 2020.

7. Salvat F., Martinez J. D., Mayol R., Parellada J., Analytical Dirac-Hartree-Fock-Slater screening function for atoms (Z = 1-92), Phys. Rev. A, 36 (2), 467-474, 1987.

THE AUTHORS

KHARLAMOVA Anastasya MAKAROV Dmitry

Kharlamova.anastasya2015@yandex.ru d.makarov@narfu.ru

ORCID: 0000-0002-8192-615X ORCID: 0000-0003-1999-1433

Received 11.07.2022. Approved after reviewing 24.07.2022. Accepted 25.07.2022.

© Peter the Great St. Petersburg Polytechnic University, 2022