ADMET and Solubility Analysis of New 5-Nitroisatine-Based Inhibitors of CDK2 Enzymes.

The development of new substances with the ability to interact with a biological target is only the first stage in the process of the creation of new drugs. The 5-nitroisatin derivatives considered in this study are new inhibitors of cyclin-dependent kinase 2 (CDK2) intended for anticancer therapy. The research, carried out based on the ADMET (absorption, distribution, metabolism, excretion, toxicity) methods, allowed a basic assessment of the physicochemical parameters of the tested drugs to be made. The collected data clearly showed the good oral absorption, membrane permeability, and bioavailability of the tested substances. The analysis of the metabolite activity and toxicity of the tested drugs did not show any critical hazards in terms of the toxicity of the tested substances. The substances' low solubility in water meant that extended studies tested compounds were required, which helped to select solvents with a high dissolving capacity of the examined substances, such as DMSO or NMP. The use of aqueous binary mixtures based on these two solvents allowed a relatively high solubility with significantly reduced toxicity and environmental index compared to pure solvents to be maintained, which is important in the context of the search for green solvents for pharmaceutical use.

Will the Interactions of Some Platinum (II)-Based Drugs with B-Vitamins Reduce Their Therapeutic Effect in Cancer Patients? Comparison of Chemotherapeutic Agents such as Cisplatin, Carboplatin and Oxaliplatin-A Review.

Pt (II) derivatives show anti-cancer activity by interacting with nucleobases of DNA, thus causing some spontaneous and non-spontaneous reactions. As a result, mono- and diaqua products are formed which further undergo complexation with guanine or adenine. Consequently, many processes are triggered, which lead to the death of the cancer cell. The theoretical and experimental studies confirm that such types of interactions can also occur with other chemical compounds. The vitamins from B group have a similar structure to the nucleobases of DNA and have aromatic rings with single-pair orbitals. Theoretical and experimental studies were performed to describe the interactions of B vitamins with Pt (II) derivatives such as cisplatin, oxaliplatin and carboplatin. The obtained results were compared with the values for guanine. Two levels of simulations were implemented at the theoretical level, namely, B3LYP/6-31G(d,p) with LANL2DZ bases set for platinum atoms and MN15/def2-TZVP. The polarizable continuum model (IEF-PCM preparation) and water as a solvent were used. UV-Vis spectroscopy was used to describe the drug-nucleobase and drug-B vitamin interactions. Values of the free energy (ΔGr) show spontaneous reactions with mono- and diaqua derivatives of cisplatin and oxaliplatin; however, interactions with diaqua derivatives are more preferable. The strength of these interactions was also compared. Carboplatin products have the weakest interaction with the studied structures. The presence of non-covalent interactions was demonstrated in the tested complexes. A good agreement between theory and experiment was also demonstrated.

Affinities to Oxaliplatin: Vitamins from B Group vs. Nucleobases.

Oxaliplatin, similar to Cisplatin, exhibits anticancer activity by interacting with DNA and inducing programmed cell death. It is biotransformed through a number of spontaneous and non-enzymatic processes. In this way, several transient reactive species are formed, including dichloro-, monochloro-, and diaqua-DACH platin, which can complex with DNA and other macromolecules. The molecular level suggests that such interactions can also take place with vitamins containing aromatic rings with lone pair orbitals. Theoretical and experimental studies were performed to investigate interactions of vitamins from the B group with Oxaliplatin, and the results were compared with values characterizing native purines. Quantum-chemical simulations were carried out at the B3LYP/6-31G(d,p) level, with the LANL2DZ basis set representing atomic orbitals of platinum atom, and at the MN15/def2-TZVP levels of theory with the use of Polarizable Continuum Model (IEF-PCM formulation) and water as a solvent. Additionally, time-dependent density functional theory (TD-DFT) was employed to study molecular properties in the electronic excited state. Interactions of vitamins and Oxaliplatin were investigated using UV-Vis spectroscopy. Values of the free energy (ΔGr) indicate spontaneous reactions with monoaqua [PtH2OClDACH]+ and diaqua [Pt(H2O)2DACH]2+ derivatives of Oxaliplatin. However, diaqua derivatives were found to be preferable. The free energy (ΔGr) values obtained for vitamins from the B group indicate lower affinity of Oxaliplatin compared with values characterizing complexes formed by guanine, adenine, and cytosine. The exception is the monoaqua form of vitamin B1 (thiamine) at the MN15/def2-TZVP levels of calculations. An application of atoms in molecules (AIM) theory revealed non-covalent interactions present in the complexes studied. The comparison of computed and experimental spectroscopic properties showed a good agreement.

Designing and Synthesis of New Isatin Derivatives as Potential CDK2 Inhibitors.

Tumors are still one of the main causes of death; therefore, the search for new therapeutic agents that will enable the implementation of effective treatment is a significant challenge for modern pharmacy. One of the important factors contributing to the development of neoplastic diseases is the overexpression of enzymes responsible for the regulation of cell division processes such as cyclin-dependent kinases. Numerous studies and examples of already-developed drugs confirm that isatin is a convenient basis for the development of new groups of inhibitors for this class of enzyme. Therefore, in this work, a new group of potential inhibitors of the CDK2 enzyme, utilizing isatin derivatives and substituted benzoylhydrazines, has been designed based on the application of computational chemistry methods, such as docking and molecular dynamics, and their inhibiting ability was assessed. In the cases of the selected compounds, a synthesis method was developed, and the selected physicochemical properties of the newly synthesized derivatives were estimated. As part of the completed project, new compounds are developed which are potential inhibitors of the CDK2 enzyme.

The assessment of physicochemical properties of Cisplatin complexes with purines and vitamins B group.

The positively charged products of Cisplatin hydrolysis can form bonds with Guanine and Adenine, showing the ability to crosslink with nucleobases within the double helical DNA, and leading to apoptosis of the neoplastic cell. It has been proved that the presence of chemicals other than nucleobases, compound of aromatic rings with a nitrogen lone pair on the ring, such as B vitamins, may have a competitive character in relation to a chemotherapeutic drug. A theoretical study confirms the stability of bonds formed not only between Cisplatin and Guanine/Adenine but also Cisplatin and B vitamins, namely Thiamine (vit.B1), Niacin (vit.B3), Riboflavin (vit.B2) and Pyridoxal phosphate (vit.B6). That is why it seems to be justified to conduct the research on the physicochemical, thermochemical and optical properties of mono and diaqua complexes of Cisplatin with nucleobases and B vitamins. Equally important is carrying out such research as spectroscopic measurements, bond order analysis and vibrational analysis of the studied complexes. The complexation reaction is spontaneous and thermodynamically favored with the high electronegativity value, a shift of the electron density from the metal zone towards organic compounds, a reduction the value of gap parameter and a shift of the maximum absorbance ΔλABS towards longer wavelengths. Moreover, the performed density variation upon photoexcitation showed the contributions from HOMO→LUMO transition, where the maximum of the absorption band shifts towards shorter wavelengths compared to the cisPt derivatives, thereby lowering the photoexcitation energy. The formation of complexes causes the reduction of energy gap (ΔEgap) values and show higher kinetic stability with high values of gap. The bonding energies between Cisplatin and the target molecules were performed. Cisplatin forms the strongest bonds with Guanine and Pyridoxal phosphate (vit.B6). The vibrations of Pt-N7 bond in complexes occur in the low range of frequencies with low intensities. Only in case of nucleobases appeared vibrations of Pt-N7 in high frequencies. The highest intensity is shown by symmetry and asymmetry stretching vibrations of H-O(H2O) bonds in the high range of frequencies. The vibrations related with ammonium groups (NH3) of great importance appear in the medium and high range of frequencies. The complexes with Guanine show the highest intensities. A theoretical IR spectra were studied. The obtained calculated IR spectra for native nucleobases and vitamins from B group were compared with their experimental FT-IR. So the conducted research provides theoretical knowledge in practical terms regarding the physicochemical and spectral properties of the formed Cisplatin complexes with both nucleobases and B vitamins. Theoretical and partially experimental studies of the molecular and electronic structure or prediction of spectroscopic characteristics are quite useful for better understanding of the reactivity of this drug with physiological target molecules.

Substituent and Solvent Polarity on the Spectroscopic Properties in Azo Derivatives of 2-Hydroxynaphthalene and Their Difluoroboranes Complexes.

Novel fluorescent dyes such as difluoroborane complexes of 1-phenylazonaphthalen-2-ol derivatives were successfully synthesized and characterized with a focus on the influence of a substituent and a solvent on the basic photophysical properties. 1H, 11B, 13C, 15N, and 19F nuclear magnetic resonance (NMR) spectra of substituted 1-phenylazonaphthalen-2-ol difluoroboranes and their parent azo dyes were recorded and discussed. The absorption and emission properties of synthesized compounds were investigated in solvents of varying polarity. They were found to be fluorescent despite the presence of the azo group. The azo group rotation was blocked by complexing with -BF2 to get a red shift in absorption. Solvent-dependent spectral properties of compounds were investigated using Lipper-Mataga and Bakhshiev plot. The calculated DFT energies and Frontier Molecular Orbitals calculations of the studied compounds were proved to be consistent with the experimental observations.

The Affinity of Carboplatin to B-Vitamins and Nucleobases.

Platinum compounds have found wide application in the treatment of various types of cancer and carboplatin is one of the main platinum-based drugs used as antitumor agents. The anticancer activity of carboplatin arises from interacting with DNA and inducing programmed cell death. However, such interactions may occur with other chemical compounds, such as vitamins containing aromatic rings with lone-pair orbitals, which reduces the anti-cancer effect of carboplatin. The most important aspect of the conducted research was related to the evaluation of carboplatin affinity to vitamins from the B group and the potential impact of such interactions on the reduction of therapeutic capabilities of carboplatin in anticancer therapy. Realized computations, including estimation of Gibbs Free Energies, allowed for the identification of the most reactive molecule, namely vitamin B6 (pyridoxal phosphate). In this case, the computational estimations indicating carboplatin reactivity were confirmed by spectrophotometric measurements.

The Oxindole Derivatives, New Promising GSK-3ß Inhibitors as One of the Potential Treatments for Alzheimer's Disease-A Molecular Dynamics Approach.

The glycogen synthase kinase 3ß (GSK-3ß) is a protein kinase involved in regulating numerous physiological processes such as embryonic development, transcription, insulin action, cell division cycle and multiple neuronal functions. The overexpression of this enzyme is related to many diseases such as schizophrenia, Alzheimer's disease, diabetes and cancer. One of the basic methods of treatment in these cases is the usage of ATP-competitive inhibitors. A significant group of such compounds are indirubin and its analogs, e.g., oxindole derivatives. The compounds considered in this work are 112 newly designed oxindole derivatives. In the first stage, such molecular properties of considered compounds as toxicity and LogP were estimated. The preliminary analysis of the binding capabilities of considered compounds towards the GSK-3ß active site was conducted with the use of the docking procedure. Based on obtained molecular properties and docking simulations, a selected group of complexes that were analyzed in the molecular dynamics stage was nominated. The proposed procedure allowed for the identification of compounds such as Oxind_4_9 and Oxind_13_10, which create stable complexes with GSK-3ß enzyme and are characterized by the highest values of binding affinity. The key interactions responsible for stabilization of considered ligand-protein complexes were identified, and their dynamic stability was also determined. Comparative analysis including analyzed compounds and reference molecule 3a, which is also an oxindole derivative with a confirmed inhibitory potential towards GSK3B protein, clearly indicates that the proposed compounds exhibit an analogous binding mechanism, and the obtained binding enthalpy values indicate a slightly higher binding potential than the reference molecule.

Substituted 2-Phenacylbenzoxazole Difluoroboranes: Synthesis, Structure and Properties.

Novel fluorescent dyes such as benzoxazole-boron complexes, bearing ß-ketoiminate ligands, have been synthesized and characterized with a focus on the influence of a substituent on the basic photophysical properties. 1H, 11B, 13C, 15N, and 19F nuclear magnetic resonance (NMR) spectra of substituted 2-phenacylbenzoxazole difluoroboranes have been recorded and discussed. It is worth mentioning that a high correlation coefficient was found between 15N-NMR parameters and substituent constants. The photophysical properties of these new dyes have been investigated by fluorescence and ultraviolet-visible (UV-Vis) absorption spectroscopy. The geometry optimization, vibrational spectra, and the HOMO and LUMO energies were calculated based on density functional theory with the use of the B3LYP functional and 6-311++G(d,p) basis set.

Does the Affinity of Cisplatin to B-Vitamins Impair the Therapeutic Effect in the Case of Patients with Lung Cancer-consuming Carrot or Beet Juice?

BACKGROUND: Cisplatin (CisPt) has a well-recognized anticancer activity by interacting with DNA and inducing programmed cell death. However, theoretical studies performed on the molecular level suggest that such nonspecific interactions can also take place with many competitive compounds, such as vitamins containing aromatic rings with lone-pair orbitals. OBJECTIVE: This work is a theoretical study on the initial Pt-N7(N1) bond formation with vitamins from B group and their comparison with values characterizing native purines. METHODS: Geometries of studied structures were optimized with an aid of Gaussian 09 using the B3LYP functional with the 6-31G** basis set. Atomic orbitals of platinum were represented by the lanl2dz basis. Solvation free energies were evaluated by a self-consistent reaction field (SCRF) approach. A dielectric constant of 78 for water was used in the PCM continuum model computations along with radii Bondii. RESULTS: The affinities of mono-aqua cis~[Pt(NH3)2Cl(H2O)]+ and di-aqua cis~ [Pt(NH3)2(H2O)2]2+ derivatives of Cisplatinum toward compounds belonging to the group of eight B vitamins were studied and compared to interactions with canonical purines. All the values of ΔGr unambiguously indicate that reactions with cisPt-diaqua are more preferable, but the comparison of ΔGr values obtained for compounds from vitamin B group and the ones characterizing complexes created by Guanine molecules indicates higher affinity of cisPt monomers toward purines. CONCLUSION: Based on the observations, the regular intake of vitamin-rich beetroot or carrot juices is strongly discouraged during anticancer therapy using CisPt drug. To confirm the results of the performed computational study, detailed clinical trials should be performed.

Biomolecular interactions of lysosomotropic surfactants with cytochrome c and its effect on the protein conformation: A biophysical approach.

The molecular interactions between two single-chain lysosomotropic surfactants DMM-11 (2-Dodecanoyloxyethyl)trimethylammonium bromide) and DMPM-11 (2-Dodecanoyloxypropyl)trimethylammonium bromide) with a small heme-protein (cytochrome c (cyt-c)) in Hepes buffer (pH = 7.4) were extensively investigated by surface tension, dynamic light scattering (DLS), circular dichroism (CD) and fluorescence spectroscopy in combination with molecular dynamic simulation techniques. The results demonstrated that surfactants can destroy the hydrophobic cavity of cyt-c, make the α-helical become loose and convert it into the ß-sheet structure. The interactions between surfactants and cyt-c are mainly hydrophobic. Molecular modelling approaches were also used to gather a deeper insight on the binding of lysosomotropic surfactants with cyt-c and the in silico results were found to be in good agreement with the experimental ones. This study provides a molecular basis for the applications of protein-surfactant complexes in biological, food, pharmaceutical, industrial and cosmetic systems.

Reactive group effects on the photophysical and biological properties of 2-phenyl-1H-phenanthro[9,10-d]imidazole derivatives as fluorescent markers.

The presented research focuses on the theoretical design and procedures for preparing protein conjugates with markers. For this purpose a series of phenanthroimidazole (PhI) analogous compounds was designed and investigated by means of first principle methods. Through the judicious choice of cross-linking reagents and the selection of reactive groups, five target fluorescent probes were selected, one of which was previously described using in vitro tests. For the best cognitive purpose and understanding of the nature of the protein conjugation, the studies describe the impact of the reactive group on the solvatochromism, the polarity of the charge transfer of the excited states, the Stokes' shift, ECD spectra and two-photon cross sections. The research is also extended to an analysis of PhI-Concanavalin A biocomplexes and changes in photophysical properties after conjugation. In order to identify valuable alternatives to commercial probes designed for cellular labelling in biological and biomedical imaging, biological properties were described such as ecotoxicity, log P and log BCF, and dye-protein binding was quantified by means of AutoDock and molecular dynamics simulations. The study showed that for phenanthroimidazole derivatives the factor which limits the possibility of their use in medical imaging is the presence of a pyridyl disulfide group, while the introduction of an N-hydroxysuccinimide ester may be used to create stable and valuable fluorescent probes with a wide spectrum for applications in biomedical imaging.

Inhibition mechanism of CDK-2 and GSK-3ß by a sulfamoylphenyl derivative of indoline-a molecular dynamics study.

A good understanding of the inhibition mechanism of enzymes exhibiting high levels of similarity is the first step to the discovery of new drugs with selective potential. Examples of such proteins include glycogen synthase kinase-3 (GSK-3ß) and cyclin-dependent kinase 2 (CDK-2). This article reports the mechanism of such enzyme inhibition as analyzed by an indoline sulfamylophenyl derivative (CHEMBL410072). Previous work has shown that such compounds exhibit selective properties towards their biological targets. This study used a combined procedure involving docking and molecular dynamics simulations, which allowed identification of interactions responsible for stabilization of complexes, and analysis of the dynamic stability of the systems obtained. The initial data obtained during the molecular docking stage show that the ligand molecule exhibits a similar affinity towards both active sites, which was confirmed by quantification of identified interactions and energy values. However, the data do not cover dynamic aspects of the considered systems. Molecular dynamics simulations realized for both complexes indicate significant dissimilarities in dynamics properties of both side chains of the considered ligands, especially in the case of the part containing the sulfamide group. Such increased mobility of the analyzed systems disrupts the stability of binding in the stabilized complex with GSK-3ß protein, which finally affects also the binding affinity of the ligand molecule towards this enzyme.

Potential inhibitory effect of indolizine derivatives on the two enzymes: nicotinamide phosphoribosyltransferase and beta lactamase, a molecular dynamics study.

Nicotinamide phosphoribosyl-transferases (NAMPT) are enzymes that play a role in targeting cancer metabolism, while beta lactamases are involved in bacterial resistance to beta-lactam antibiotics. Many protein inhibitors exhibit such property which is often correlated with their cellular potency. In order to understand such a phenomenon, the present article conducts an analysis of the dynamic behavior of complexes formed by the inhibitors, that is indolizine derivatives, with the studied enzymes. Both docking and molecular dynamics led to identification of their interactions and showed the mechanism of inhibition of the two studied enzymes. The differences in the behavior of ligand at the active sites of beta lactamases and nicotinamide phosphoribosyl-transferases are indicated by structural and enthalpy values.

Understanding the Action of Indolizines as Biologically Active Moieties: A Molecular Dynamics Study.

BACKGROUND: Indolizines represent a class of heteroaromatic compounds, of pharmacological importance, containing two condensed 5- and 6-memebered rings bridged by a nitrogen atom. Despite indolizine is an important medicinal moiety, a detailed view on the mechanism of action of biologically active indolizines is unavailable. OBJECTIVE: The study of ligand-enzyme affinity is of high interest; description of characteristics (energetic and geometric ones) of ligand binding to the active sites of an enzyme could be useful in understanding the action mechanism of a given ligand on the concerned enzyme. METHOD: After conducting a QSAR study, to predict IC50 (on 15-LO protein from soybeans) of indolizine derivatives and a docking study of indolizines on Beta lactamase and Nicotinamide phosphoribosyltransferase proteins [1], a molecular dynamics analysis was performed on one of the indolizine derivatives, complexed to the above proteins. RESULTS: The performed molecular dynamics study led to the identification of interactions responsible for the stabilization of complexes of the chosen ligand (i.e., indolizine derivative) with the considered enzymes and the specificity of the ligand interaction as well. The structural data and enthalpy values clearly indicate the differences in the behavior of ligand at the active sites of the three investigated enzymes. Among the studied proteins, the hydrophobicity of the active site of Nicotinamide phosphoribosyltransferase seems to be the main factor in promoting the interaction enzyme-ligand, much more manifested in this case, in comparison to the other two proteins Beta lactamase and Nicotinamide phosphoribosyltransferase. CONCLUSION: The present paper discusses a possible mechanism of interaction of an indolizine derivative with three enzyme proteins, providing information for future work in this topic.

Physical nature of intermolecular interactions inside Sir2 homolog active site: molecular dynamics and ab initio study.

In the present study, we analyze the interactions of NAD+-dependent deacetylase (Sir2 homolog yeast Hst2) with carba-nicotinamide-adenine-dinucleotide (ADP-HPD). For the Sir2 homolog, a yeast Hst2 docking procedure was applied. The structure of the protein-ADP-HPD complex obtained during the docking procedure was used as a starting point for molecular dynamics simulation. The intermolecular interaction energy partitioning was performed for protein-ADP-HPD complex resulting from molecular dynamics simulation. The analysis was performed for ADP-HPD and 15 amino acids forming a deacetylase binding pocket. Although the results indicate that the first-order electrostatic interaction energy is substantial, the presence of multiple hydrogen bonds in investigated complexes can lead to significant value of induction component.

Molecular dynamics study of the inhibitory effects of ChEMBL474807 on the enzymes GSK-3ß and CDK-2.

Indirubin derivatives and analogs comprise a significant group of ATP-competitive inhibitors. The inhibitory effects of ChEMBL474807 (1-(4-amino-1,2,5-oxadiazol-3-yl)-5-(piperidin-1-ylmethyl)-N'-(pyridin-4-ylmethylene)-1H-1,2,3-triazole-4-carbohydrazide) on two enzymes, namely glycogen synthase kinase-3ß (GSK-3ß) and cyclin-dependent kinase-2 (CDK-2), were analyzed. The close resemblance of the amino acid sequences of these two enzymes (with 25% identity and 41% similarity) explains why indirubin derivatives are inhibitors of both of the enzymes studied. The docking and molecular dynamics investigation performed here led to the identification of the interactions responsible for stabilizing the ligand ChEMBL474807 at the active sites of the enzymes considered. The structural and energetic data collected during our investigations clearly indicate that there are important differences in the behavior of the ligand at the two active sites investigated here.

Structural and energetic properties of canonical and oxidized telomeric complexes studied by molecular dynamics simulations.

The structural and energetic properties of native and oxidized telomeric complexes were defined by means of molecular dynamic (MD) simulations. As a starting point, the experimental conformation of B-DNA d(GpTpTpApGpGpGpTpTpApGpGpG) oligomer bound to human protein telomeric repeat binding factor 1 (TRF1) was used. The influence on the stability of the telomeric complex of the presence of 8-oxoguanine (8oxoG) in the central telomeric triad (CTT) was estimated based on trajectories collected during 130 ns MD runs. The data obtained indicate that the system analyzed is highly sensitive to the presence of oxidative damage in the CTT of the B-DNA telomeric sequence. The most important changes were observed in the immediate vicinity of the 8-oxoguanine molecule. The significantly higher mobility of arginine 425 interacting directly with the oxidized guanine molecule has a large influence on the structural, dynamic and energetic properties of neighboring amino acids. Local changes observed for individual hydrogen bonded interactions localized in the major groove of B-DNA also have significant impact on the properties of hydrophobic clusters, which are the second type of force responsible for stability of the studied bio-system. All the changes reported in detail here unambiguously indicate a significant decrease in telomer binding affinity after oxidation.

Synthesis and structural characterization of substituted 2-phenacylbenzoxazoles.

1H and 13C NMR spectra of eleven 2-phenacylbenzoxazoles (ketimine form) show that their CDCl(3)-solutions contains also (Z)-2-(benzo[d]oxazol-2-yl)-1-phenylethenols (enolimine form). Intramolecular hydrogen bonding in the latter tautomer was found to be significantly weaker than that one in respective (Z)-2-(2-hydroxy-2-phenylvinyl)pyridines. Integrals of the 1H NMR signals were used to evaluate the molar ratio of the tautomers. Strong electron-donating substituents were found to stabilize the ketimine tautomer. pK(T) (negative logarithm of the equilibrium constant, K(T) = [ketimine]/[enolimine]) was found to be linearly dependent on the Hammett substituent constant σ. The results of the MP2 ab initio calculations reveal enolimine including an intramolecular OH…N hydrogen bond to be the most stable form both with electron-donor and electron-acceptor substituents. The stability of ketimines is an intermediate of those found for enolimines and enaminones i.e., (E)-2-(benzo[d]oxazol-2(3H)-ylidene)-1-phenylethanones. 13C CPMAS NMR spectral data reveal that in the crystalline state the ketimine tautomer is predominant in p-NMe(2) substituted congener. On the other hand, enolimine forms were detected there when the substituent has less electron-donating character or when it is an electron-acceptor by character.

Structural and energetic consequences of oxidation of d(ApGpGpGpTpT) telomere repeat unit in complex with TRF1 protein.

The configuration hyperspace of canonical and oxidized 14-mers of B-DNA comprising telomere repeat units d(ApGpGpGpTpT) was sampled over 40 ns via molecular dynamic (MD) simulations. The energetic and structural consequences of TRF1 binding to telomere B-DNA were compared with non-complexed systems. Energetic properties of analyzed pairs, di- and tri-nucleotide steps occurring in central telomere repeat unit were estimated by means of advanced quantum chemistry computations including not only BSSE corrections, electron correlation contributions but also non-negligible many-body terms. These data along with bases pair and base step parameters distributions allow for quantization of consequences of oxidation and/or TRF1 binding to telomere repeat units. Occurrence of 8-oxoguanine in central telomeric triad (CTT) is the source of high stiffness if compared to non-modified oligomer. The origin of this property comes from significantly alteration of intermolecular interactions introduced by 8-oxoguanine. The increased stability observed for base-base interactions are accumulated and characterizes also di- and tri-nucleotides. The observed changes in the intermolecular interactions originate from structural alterations imposed by TRF1 binding to canonical and oxidized telomere B-DNA. First and most direct consequence of TRF1 binding to oxidized telomere repeat unit is alteration of shift-slide correlations if compared to canonical system. This in turn leads to large differences in purine-purine overlapping in oxidized structures. Thus, oxidized telomere B-DNA double strands are sensitive to interactions with protein ligands and numerous structural and energetic changes are imposed on base pairs forming CTT.