Excited States in Single-Stranded and i-Motif DNA with Silver Ions.

We have studied the excited states and structural properties for the complexes of cytosine (dC)10 chains with silver ions (Ag+) in a wide range of the Ag+ to DNA ratio (r) and pH conditions using circular dichroism, steady-state absorption, and fluorescence spectroscopy along with the ultrafast fluorescence upconversion technique. We also calculated vertical electronic transition energies and determined the nature of the corresponding excited states in some models of the cytosine-Ag+ complexes. We show that (dC)10 chains in the presence of silver ions form a duplex stabilized by C-Ag+-C bonds. It is also shown that the i-motif structure formed by (dC)10 chains is destabilized in the presence of Ag+ ions. The excited-state properties in the studied complexes depend on the amount of binding ions and the binding sites, which is supported by the calculations. In particular, new low-lying excited states appear when the second Ag+ ion interacts with the O atom of cytosine in the C-Ag+-C pairs. A similar picture is observed in the case when one Ag+ ion interacts with one cytosine via the N7 atom.

Interactions of deprotonated phenylalanine with gold Clusters: Theoretical study with prospects for amino acid detection.

Nanomaterials are widely used nowadays in industry and medicine. The specific properties of gold nanoclusters (Au NCs) are chemical stability, low cytotoxicity, low photobleaching, high sensitivity to the molecular environment. This set of properties allows to use Au NCs as nanosensors in bioimaging and diagnostics. We have investigated gold cluster complexes with proteinogenic amino acid phenylalanine (Phe). Detection of phenylalanine is essential for diagnostics of phenylketonuria, vitiligo, sclerosis, cancer, tuberculosis, etc. We have studied the complexes of Phe with Aunq clusters with atomic number equal 1-6, 8, 20 and a charge equal 0-2. We have established that the clusters Au40, Au21+ and Au32+ form the most stable complexes with Phe among NCs with charge 0, +1 and + 2, respectively. Intracomplex interactions have been studied using Atoms-In-Molecules (AIM) theory and Natural Bond Orbital (NBO) analysis. It has been shown that metal-ligand intracomplex interactions are partially covalent and partially electrostatic. Also, we have simulated the UV-vis absorption and Raman spectra of the Phe-Au NCs. We have established that the clusters possess prospective features if being used for colorimetric and Raman detection of Phe. Au20 cluster is remarkable for its six-times enhancement of the Raman signal. Moreover, our study provides insights into metal-ligand interactions for clusters synthesized inside a polypeptide globula. Hence, to the best of our knowledge this is a first attempt to perform a detailed analysis of Phe interactions with gold using quantum chemical calculations.

Photosensitized isomerization of resveratrol: Evaluation of energy and electron transfer pathways.

Resveratrol (3,5,4'-trihydroxystilbene, RSV) is a natural stilbene synthetized as trans-isomer in plants exposed to oxidative stress. In order to understand the mechanism involved during photosensitized degradation of trans-resveratrol, steady-state and time-resolved experiments were performed and compared with quantum-chemical calculations using density functional theory (DFT). Pterin (Ptr), a well-known photosensitizer, under UV-A radiation induces the oxidation of several biomolecules mainly through electron-transfer mechanisms. On the one hand, it was observed that trans-RSV participates in an energy-transfer pathway with Ptr triplet excited state (3Ptr*) forming 3trans-RSV*, which dissipates the energy by isomerization to cis-RSV. On the other hand, RSV neutral radical (trans-RSV(-H)•) was detected in laser flash photolysis experiments, evidencing an electron-transfer mechanism. The electron-transfer from 3Ptr* to trans-RSV is a barely feasible reaction, however, more favorable is the formation of trans-RSV(-H)• in a reaction between trans-RSV and Ptr radical cation (Ptr•+), which is produced during irradiation. The combination of experimental and theoretical approaches evidences the capability of trans-RSV to undergo energy-transfer (feasible by DFT calculations) and/or one-electron transfer pathways with 3Ptr*. These findings reveal the mechanisms involved in the interaction of trans-RSV and pterin excited states and provide information on the antioxidant action of resveratrol during photosensitized oxidation of biomolecules.

Tetrahydrobiopterin as a Trigger for Vitiligo: Phototransformation during UV Irradiation.

Vitiligo is a type of hypomelanosis. Tetrahydrobiopterin (H4Bip), the coenzyme of the initial stage of melanogenesis, appears to be a trigger for vitiligo. H4Bip is present in vitiligo in 3-5-fold excess and causes oxidative stress by triggering an autocatalytic cycle of excess hydrogen peroxide synthesis. Using quantum-chemical calculations, we have evaluated the possibility of H4Bip reactions occurring in the dark and under ultraviolet (UV) irradiation, including the formation of dihydropterin dimers. In order to simulate the oxidative stress, oxidative modification of human serum albumin (HSA) has been carried out in the presence of excessive H4Bip using the fluorescence method. The fraction of oxidized protein (FOP) has been calculated. It has been established that there is a strong oxidative modification of amino acids chromophores (tryptophan and tyrosine) in the protein (FOP 0.64). Under UV irradiation of the system (HSA + H4Bip), FOP is reduced to 0.39. Apparently, a part of H4Bip transforms into dihydropterin dimers and does not participate in the oxidative modification of the protein. The data on oxidative modification of HSA are consistent with dynamic light scattering: H4Bip promotes HSA aggregation with the formation of particles with a hydrodynamic radius Rh ≥ 2000 nm, which can become immunogenic.

Quantitative determination of albumin and immunoglobulin in human serum using gold nanoclusters.

In this experimental study, we developed a simple and selective approach to determine the concentrations of human serum albumin (HSA) and total amount of immunoglobulins (Ig) in real human serum (HS) sample using luminescent gold nanoclusters (Au NCs). In doing so, Au NCs were grown directly on the HS proteins without any sample pretreatment. We synthesized Au NCs on HSA and Ig and studied their photophysical properties. Using combined fluorescent and colorimetric assay we were able to obtain protein concentrations with a high degree of accuracy relative to techniques currently used in clinical diagnostics. We used method of standard additions to determine both HSA and Ig concentrations in HS by the Au NCs absorbance and fluorescence signals. A simple and cost-effective method developed in this work represents an excellent alternative to the techniques currently used in clinical diagnostics.

Rapid and selective colorimetric determination of L-DOPA in human serum with silver nanoparticles.

L-DOPA, or l-3,4-dihydroxyphenylalanine is an aromatic amino acid, which plays a significant role in human metabolism as a precursor of important neurotransmitters. We develop a fast and simple colorimetric method for the detection of L-DOPA in biological fluids. The method is based on the reduction of silver ions with L-DOPA and the subsequent formation of L-DOPA stabilized silver nanoparticles (Ag NPs). In this novel approach, L-DOPA works as both reducing and stabilizing agent, which provides selectivity and simplifies the procedure. HR-TEM images show very narrow Ag NPs distribution with an average size of 24 nm. Such sensor design is suggested for the first time. We also calculate vertical ionization potential, vertical electron affinity, and Gibbs free energy change of different ionic forms of L-DOPA and amino acids at the M06-2X/def2-TZVP level for the gas phase in comparison with that of silver. A model of silver ions reduction by aromatic amino acids is proposed: the ionic forms with charge -1 are suggested to reduce silver ions. High selectivity against aromatic amino acids, dopamine and serotonin is achieved by tuning pH and involving two L-DOPA forms with charged both hydroxyphenolate and carboxylate groups in the stabilization of uniform-sized Ag NPs. The method is applicable for the determination of L-DOPA in human serum with the 50 nM limit of detection and the linear range up to 5 µM. Ag NPs formation and coloring the solution proceeds in a few minutes. The suggested colorimetric method has potential application in clinical trials.

Insights into Molecular Structure of Pterins Suitable for Biomedical Applications.

Pterins are an inseparable part of living organisms. Pterins participate in metabolic reactions mostly as tetrahydropterins. Dihydropterins are usually intermediates of these reactions, whereas oxidized pterins can be biomarkers of diseases. In this review, we analyze the available data on the quantum chemistry of unconjugated pterins as well as their photonics. This gives a comprehensive overview about the electronic structure of pterins and offers some benefits for biomedicine applications: (1) one can affect the enzymatic reactions of aromatic amino acid hydroxylases, NO synthases, and alkylglycerol monooxygenase through UV irradiation of H4pterins since UV provokes electron donor reactions of H4pterins; (2) the emission properties of H2pterins and oxidized pterins can be used in fluorescence diagnostics; (3) two-photon absorption (TPA) should be used in such pterin-related infrared therapy because single-photon absorption in the UV range is inefficient and scatters in vivo; (4) one can affect pathogen organisms through TPA excitation of H4pterin cofactors, such as the molybdenum cofactor, leading to its detachment from proteins and subsequent oxidation; (5) metal nanostructures can be used for the UV-vis, fluorescence, and Raman spectroscopy detection of pterin biomarkers. Therefore, we investigated both the biochemistry and physical chemistry of pterins and suggested some potential prospects for pterin-related biomedicine.

Silver cluster interactions with Pterin: Complex structure, binding energies and spectroscopy.

Metal nanoclusters (NCs) are widely present today in biosensing, bioimaging, and diagnostics due to their small size, great biocompatibility, and sensitivity to the biomolecular environment. Silver (Ag) NCs often possess intense fluorescence, photostability, and low photobleaching, which is in high demand during the detection of organic molecules. Pterins are small compounds, which are used in medicine as biomarkers of oxidative stress, cardiovascular diseases, neurotransmitter synthesis, inflammation and immune system activation. It is experimentally possible to detect pterin (Ptr) through the adsorption on Ag colloid. We optimized geometries and evaluated the binding energy in Ptr-Agnq complexes (n = 1-6; q = 0, +1, +2) using quantum chemistry methods. Different Ptr atoms were preferential for silver attachment depending on NC charge and size. The highest Eb was obtained for the complexes between the Ptr0 and Ag32+ (-50.8 kcal mol-1), between Ptr-1 and Ag32+ (-64.8 kcal mol-1), which means that these complexes should be formed preferably in aqueous solutions in acidic and alkaline media, respectively. The colorimetric detection of pterin with silver clusters does not seem to be promising. However, intense S0→S1 transitions of Ag5+ complexes look promising for luminescent Ptr detection. SERS detection of pterin is better to be done at pH > 8 since deprotonated pterin Raman undergo more dramatic changes upon addition of Ag than the neutral pterin. The characteristics of absorption and vibrational spectra of silver-pterin should be exploited during biosensor development.

Silver Cluster Interactions with Tyrosine: Towards Amino Acid Detection.

Tyrosine (Tyr) is involved in the synthesis of neurotransmitters, catecholamines, thyroid hormones, etc. Multiple pathologies are associated with impaired Tyr metabolism. Silver nanoclusters (Ag NCs) can be applied for colorimetric, fluorescent, and surface-enhanced Raman spectroscopy (SERS) detection of Tyr. However, one should understand the theoretical basics of interactions between Tyr and Ag NCs. Thereby, we calculated the binding energy (Eb) between Tyr and Agnq (n = 1-8; q = 0-2) NCs using the density functional theory (DFT) to find the most stable complexes. Since Ag NCs are synthesized on Tyr in an aqueous solution at pH 12.5, we studied Tyr-1, semiquinone (SemiQ-1), and Tyr-2. Ag32+ and Ag5+ had the highest Eb. The absorption spectrum of Tyr-2 significantly red-shifts with the attachment of Ag32+, which is prospective for colorimetric Tyr detection. Ag32+ interacts with all functional groups of SemiQ-1 (phenolate, amino group, and carboxylate), which makes detection of Tyr possible due to band emergence at 1324 cm-1 in the vibrational spectrum. The ground state charge transfer between Ag and carboxylate determines the band emergence at 1661 cm-1 in the Raman spectrum of the SemiQ-1-Ag32+ complex. Thus, the prospects of Tyr detection using silver nanoclusters were demonstrated.

UV Radiation in DNA Damage and Repair Involving DNA-Photolyases and Cryptochromes.

Prolonged exposure to ultraviolet radiation on human skin can lead to mutations in DNA, photoaging, suppression of the immune system, and other damage up to skin cancer (melanoma, basal cell, and squamous cell carcinoma). We reviewed the state of knowledge of the damaging action of UVB and UVA on DNA, and also the mechanisms of DNA repair with the participation of the DNA-photolyase enzyme or of the nucleotide excision repair (NER) system. In the course of evolution, most mammals lost the possibility of DNA photoreparation due to the disappearance of DNA photolyase genes, but they retained closely related cryptochromes that regulate the transcription of the NER system enzymes. We analyze the published relationships between DNA photolyases/cryptochromes and carcinogenesis, as well as their possible role in the prevention and treatment of diseases caused by UV radiation.

Quantitative Structure-Property Relationship Modelling for the Prediction of Singlet Oxygen Generation by Heavy-Atom-Free BODIPY Photosensitizers*.

Heavy-atom-free sensitizers forming long-living triplet excited states via the spin-orbit charge transfer intersystem crossing (SOCT-ISC) process have recently attracted attention due to their potential to replace costly transition metal complexes in photonic applications. The efficiency of SOCT-ISC in BODIPY donor-acceptor dyads, so far the most thoroughly investigated class of such sensitizers, can be finely tuned by structural modification. However, predicting the triplet state yields and reactive oxygen species (ROS) generation quantum yields for such compounds in a particular solvent is still very challenging due to a lack of established quantitative structure-property relationship (QSPR) models. In this work, the available data on singlet oxygen generation quantum yields (ΦΔ ) for a dataset containing >70 heavy-atom-free BODIPY in three different solvents (toluene, acetonitrile, and tetrahydrofuran) were analyzed. In order to build reliable QSPR model, a series of new BODIPYs were synthesized that bear different electron donating aryl groups in the meso position, their optical and structural properties were studied along with the solvent dependence of singlet oxygen generation, which confirmed the formation of triplet states via the SOCT-ISC mechanism. For the combined dataset of BODIPY structures, a total of more than 5000 quantum-chemical descriptors was calculated including quantum-chemical descriptors using density functional theory (DFT), namely M06-2X functional. QSPR models predicting ΦΔ values were developed using multiple linear regression (MLR), which perform significantly better than other machine learning methods and show sufficient statistical parameters (R=0.88-0.91 and q2 =0.62-0.69) for all three solvents. A small root mean squared error of 8.2 % was obtained for ΦΔ values predicted using MLR model in toluene. As a result, we proved that QSPR and machine learning techniques can be useful for predicting ΦΔ values in different media and virtual screening of new heavy-atom-free BODIPYs with improved photosensitizing ability.

Autoxidation and photooxidation of tetrahydrobiopterin: a theoretical study.

Pterins are naturally occurring pigments and enzyme cofactors widespread in living organisms. Tetrahydrobiopterin (H4Bip) is a coenzyme of aromatic amino acid hydroxylases, NO-synthases, and alkylglycerol monooxygenases. This coenzyme is prone to oxidation in the presence of molecular oxygen, a so-called autoxidation. The reactions participating in H4Bip autoxidation are well known. However, our study is an attempt to evaluate theoretically the feasibility of reactions participating in autoxidation. To do so, we have calculated the Gibbs free energy of elementary reactions between H4Bip, its derivatives, molecular oxygen, and reactive oxygen species (ROS). In the last few years, we have established the photosensitized oxidation of H4Bip experimentally. Thus, we have also evaluated the feasibility of H4Bip photooxidation reactions, which may occur according to both type-I and type-II photosensitized oxidation. We calculated Fukui indices for H4Bip and found particular atoms in the molecule that interact with nucleophiles (for example, singlet oxygen 1O2) and radicals (in particular, molecular oxygen 3O2). Therefore, we evaluated the probability of H4Bip autoxidation reactions, photooxidation reactions, and the reactivity of particular atoms in H4Bip molecule using the theoretical methods of quantum chemistry.

Methods and Applications of In Silico Aptamer Design and Modeling.

Aptamers are nucleic acid analogues of antibodies with high affinity to different targets, such as cells, viruses, proteins, inorganic materials, and coenzymes. Empirical approaches allow the design of in vitro aptamers that bind particularly to a target molecule with high affinity and selectivity. Theoretical methods allow significant expansion of the possibilities of aptamer design. In this study, we review theoretical and joint theoretical-experimental studies dedicated to aptamer design and modeling. We consider aptamers with different targets, such as proteins, antibiotics, organophosphates, nucleobases, amino acids, and drugs. During nucleic acid modeling and in silico design, a full set of in silico methods can be applied, such as docking, molecular dynamics (MD), and statistical analysis. The typical modeling workflow starts with structure prediction. Then, docking of target and aptamer is performed. Next, MD simulations are performed, which allows for an evaluation of the stability of aptamer/ligand complexes and determination of the binding energies with higher accuracy. Then, aptamer/ligand interactions are analyzed, and mutations of studied aptamers made. Subsequently, the whole procedure of molecular modeling can be reiterated. Thus, the interactions between aptamers and their ligands are complex and difficult to understand using only experimental approaches. Docking and MD are irreplaceable when aptamers are studied in silico.

Comparative study of gold and silver interactions with amino acids and nucleobases.

Metal nanoclusters (NCs) have gained much attention in the last decade. In solution, metal nanoclusters can be stabilized by proteins, and, thus, exhibit many advantages in biocatalysis, biosensing, and bioimaging. In spite of much progress in the synthesis of polypeptide-stabilized gold (Au) clusters, their structure, as well as amino acid-cluster and amino acid-Au+ interactions, remain poorly understood. It is not entirely clear which amino acid (AA) residues and sites in the protein are preferred for binding. The understanding of NC-protein interactions and how they evolve in the polypeptide templates is the key to designing Au NCs. In this work, binding of gold ion Au+ and diatomic neutral gold nanocluster Au2 with a full set of α-proteinogenic amino acids is studied using Density Functional Theory (DFT) and the ab initio RI-MP2 method in order to find the preferred sites of gold interaction in proteins. We demonstrated that the interaction of gold cations and clusters with protonated and deprotonated amino acid residues do not differ greatly. The binding affinity of AAs to the Au2 cluster increases in the following order: Cys(-H+) > Asp(-H+) > Tyr(-H+) > Glu(-H+) > Arg > Gln, His, Met ≫ Asn, Pro, Trp > Lys, Tyr, Phe > His(+H+) > Asp > Lys(+H+) > Glu, Leu > Arg(+H+) > Ile, Val, Ala > Thr, Ser > Gly, Cys, which agrees with the available experimental data that gold cluster synthesis occurs in a wide range of pH - amino acid residues with different protonation states are involved in this process. The significant difference in the binding energy of metal atoms with nucleobases and amino acids apparently means that unlike on DNA templates, neutral metal atoms are strongly bound to amino acid residues and can't freely diffuse in a polypeptide globula. This fact allows one to conclude that formation of metal NCs in proteins occurs through the nucleation of reduced Au atoms bound to the neighboring amino acid residues, and the flexibility of the amino acid residue side-chains and protein chain as a whole plays a significant role in this process.

Nano-(Q)SAR for Cytotoxicity Prediction of Engineered Nanomaterials.

Although nanotechnology is a new and rapidly growing area of science, the impact of nanomaterials on living organisms is unknown in many aspects. In this regard, it is extremely important to perform toxicological tests, but complete characterization of all varying preparations is extremely laborious. The computational technique called quantitative structure-activity relationship, or QSAR, allows reducing the cost of time- and resource-consuming nanotoxicity tests. In this review, (Q)SAR cytotoxicity studies of the past decade are systematically considered. We regard here five classes of engineered nanomaterials (ENMs): Metal oxides, metal-containing nanoparticles, multi-walled carbon nanotubes, fullerenes, and silica nanoparticles. Some studies reveal that QSAR models are better than classification SAR models, while other reports conclude that SAR is more precise than QSAR. The quasi-QSAR method appears to be the most promising tool, as it allows accurately taking experimental conditions into account. However, experimental artifacts are a major concern in this case.

Exciton Absorption and Luminescence in i-Motif DNA.

We have studied the excited-state dynamics for the i-motif form of cytosine chains (dC)10, using the ultrafast fluorescence up-conversion technique. We have also calculated vertical electronic transition energies and determined the nature of the corresponding excited states in a model tetramer i-motif structure. Quantum chemical calculations of the excitation spectrum of a tetramer i-motif structure predict a significant (0.3 eV) red shift of the lowest-energy transition in the i-motif form relative to its absorption maximum, which agrees with the experimental absorption spectrum. The lowest excitonic state in i-(dC)10 is responsible for a 2 ps red-shifted emission at 370 nm observed in the decay-associated spectra obtained on the femtosecond time-scale. This delocalized (excitonic) excited state is likely a precursor to a long-lived excimer state observed in previous studies. Another fast 310 fs component at 330 nm is assigned to a monomer-like locally excited state. Both emissive states form within less than the available time resolution of the instrument (100 fs). This work contributes to the understanding of excited-state dynamics of DNA within the first few picoseconds, which is the most interesting time range with respect to unraveling the photodamage mechanism, including the formation of the most dangerous DNA lesions such as cyclobutane pyrimidine dimers.

QSAR analysis of immune recognition for triazine herbicides based on immunoassay data for polyclonal and monoclonal antibodies.

A common task in the immunodetection of structurally close compounds is to analyze the selectivity of immune recognition; it is required to understand the regularities of immune recognition and to elucidate the basic structural elements which provide it. Triazines are compounds of particular interest for such research due to their high variability and the necessity of their monitoring to provide safety for agricultural products and foodstuffs. We evaluated the binding of 20 triazines with polyclonal (pAb) and monoclonal (mAb) antibodies obtained using atrazine as the immunogenic hapten. A total of over 3000 descriptors were used in the quantitative structure-activity relationship (QSAR) analysis of binding activities (pIC50). A comparison of the two enzyme immunoassay systems showed that the system with pAb is much easier to describe using 2D QSAR methodology, while the system with mAb can be described using the 3D QSAR CoMFA. Thus, for the 3D QSAR model of the polyclonal antibodies, the main statistical parameter q2 ('leave-many-out') is equal to 0.498, and for monoclonal antibodies, q2 is equal to 0.566. Obviously, in the case of pAb, we deal with several targets, while in the case of mAb the target is one, and therefore it is easier to describe it using specific fields of molecular interactions distributed in space.

Silver cluster-amino acid interactions: a quantum-chemical study.

Binding of silver ion (Ag+) and two atomic neutral silver cluster (Ag2) with a set of amino acids has been studied using Density Functional Theory (DFT) and ab initio MP2 method. We show that binding energy with Ag2 is higher for deprotonated anionic amino acids. Cysteine, aspartic acid, and tyrosine with deprotonated side chain exhibit the highest binding energy (Gbind) values among all the amino acids: - 30.1 kcal mol-1, - 30.7 kcal mol-1, and - 30.9 kcal mol-1, respectively. Binding energies of deprotonated cysteine, glutamic acid, tyrosine, and aspartic acid with silver ion Ag+ are reported here for the first time. Natural bond orbital (NBO) analysis has also been performed to calculate charge transfer, natural populations, occupancies, and Wiberg bond indices of the amino acid-Ag2 complexes. Atoms-in-molecules (AIM) theory was used to establish the nature of interactions. It was shown that, in most cases, the bonds between amino acid and Ag2 cluster are partially electrostatic and partially covalent.

Ciprofloxacin and Clinafloxacin Antibodies for an Immunoassay of Quinolones: Quantitative Structure⁻Activity Analysis of Cross-Reactivities.

A common problem in the immunodetection of structurally close compounds is understanding the regularities of immune recognition, and elucidating the basic structural elements that provide it. Correct identification of these elements would allow for select immunogens to obtain antibodies with either wide specificity to different representatives of a given chemical class (for class-specific immunoassays), or narrow specificity to a unique compound (mono-specific immunoassays). Fluoroquinolones (FQs; antibiotic contaminants of animal-derived foods) are of particular interest for such research. We studied the structural basis of immune recognition of FQs by antibodies against ciprofloxacin (CIP) and clinafloxacin (CLI) as the immunizing hapten. CIP and CLI possess the same cyclopropyl substituents at the N1 position, while their substituents at C7 and C8 are different. Anti-CIP antibodies were specific to 22 of 24 FQs, while anti-CLI antibodies were specific to 11 of 26 FQs. The molecular size was critical for the binding between the FQs and the anti-CIP antibody. The presence of the cyclopropyl ring at the N1 position was important for the recognition between fluoroquinolones and the anti-CLI antibody. The anti-CIP quantitative structure⁻activity relationship (QSAR) model was well-equipped to predict the test set (pred_R² = 0.944). The statistical parameters of the anti-CLI model were also high (R² = 0.885, q² = 0.864). Thus, the obtained QSAR models yielded sufficient correlation coefficients, internal stability, and predictive ability. This work broadens our knowledge of the molecular mechanisms of FQs' interaction with antibodies, and it will contribute to the further development of antibiotic immunoassays.

A quantitative structure-property relationship (QSPR) study of singlet oxygen generation by pteridines.

The QSPR method is used in photochemistry for the prediction of the absorption wavelength, fluorescence intensity, photolysis quantum yield, etc. However, to our knowledge, no attempts have been made to use the quantum yield of singlet oxygen ((1)O2) generation (ΦΔ) as an analyzed parameter in a QSPR study. We performed QSPR analysis of 29 pteridine compounds (including pterin and flavin sensitizers) for their ability to produce singlet oxygen in aqueous (D2O) solutions. Pteridines are ubiquitously present in living systems (mostly as coenzymes), possess high photochemical activity and have multiple applications as photosensitizers. Our goal was to develop a QSPR model for the fast virtual screening and prediction of the (1)O2 generation quantum yield of pteridines. Quantum-chemical descriptors were calculated using the AM1 semi-empirical method. The ability of pteridines to generate singlet oxygen was found to be significantly correlated with the HOMO orbital energy (R(2) = 0.806) and electronegativity (R(2) = 0.840). The best QSPR model obtained using electronegativity, dipole density and electrostatic charge of the N3 atom of the pteridine system allows us to predict ΦΔ of pterin and flavin photosensitizers. The model possesses high internal stability (q(2) = 0.881), as well as high predicting ability for the external dataset (pred_R(2) = 0.873). More QSPR analysis is needed for the prediction of ΦΔ of pteridines and other groups of sensitizers in aqueous as well as in non-polar solutions.