Multiple Proton-Transfer Reactions in DNA Base Pairs by Coordination of Pt Complex

Binding energies and hydrogen bonds effects on DNA-cisplatin interactions: a DFT-xTB study

CONTEXT

A systematic study of hydrogen bonds in base pairs and the interaction of cisplatin with DNA fragments was carried out. Structure, binding energies, and electron density were analyzed. xTB has proven to be an accurate method for obtaining structures and binding energies in DNA structures. Our xTB values for DNA base binding energy were in the same order and in some cases better than CAM-B3LYP values compared to experimental values. Double-stranded DNA-cisplatin structures have been calculated and the hydrogen bonds of water molecules are a decisive factor contributing to the preference for the cisplatin-Guanine interaction. Higher values of the water hydrogen bonding energies were obtained in cisplatin-Guanine structures. Furthermore, the electrostatic potential was used to investigate and improve the analysis of DNA-cisplatin structures.

METHODS

We applied the xTB method and the CAM-B3LYP functional combined with def2-SVP basis set to perform and analyze of the bonding energies of the cisplatin interaction and the effects of the hydrogen bonds. Results were calculated employing the xTB and the ORCA software.

Systematic analysis of electronic barrier heights and widths for concerted proton transfer in cyclic hydrogen bonded clusters: (HF)n, (HCl)n and (H2O)n where n = 3, 4, 5

The MP2 and CCSD(T) methods are paired with correlation consistent basis sets as large as aug-cc-pVQZ to optimize the structures of the cyclic minima for (HF)n, (HCl)n and (H2O)n where n = 3-5, as well as the corresponding transition states (TSs) for concerted proton transfer (CPT). MP2 and CCSD(T) harmonic vibrational frequencies confirm the nature of each minimum and TS. Both conventional and explicitly correlated CCSD(T) computations are employed to assess the electronic dissociation energies and barrier heights for CPT near the complete basis (CBS) limit for all 9 clusters. Results for (HF)n are consistent with prior studies identifying Cnh and Dnh point group symmetry for the minima and TSs, respectively. Our computations also confirm that CPT proceeds through Cs TS structures for the C1 minima of (H2O)3 and (H2O)5, whereas the process goes through a TS with D2d symmetry for the S4 global minimum of (H2O)4. This work corroborates earlier findings that the minima for (HCl)3, (HCl)4 and (HCl)5 have C3h, S4 and C1 point group symmetry, respectively, and that the Cnh structures are not minima for n = 4 and 5. Moreover, our computations show the TSs for CPT in (HCl)3, (HCl)4 and (HCl)5 have D3h, D2d, and C2 point group symmetry, respectively. At the CCSD(T) CBS limit, (HF)4 and (HF)5 have the smallest electronic barrier heights for CPT (≈15 kcal mol-1 for both), followed by the HF trimer (≈21 kcal mol-1). The barriers are appreciably higher for the other clusters (around 27 kcal mol-1 for (H2O)4 and (HCl)3; roughly 30 kcal mol-1 for (H2O)3, (H2O)5 and (HCl)4; up to 38 kcal mol-1 for (HCl)5). At the CBS limit, MP2 significantly underestimates the CCSD(T) barrier heights (e.g., by ca. 2, 4 and 7 kcal mol-1 for the pentamers of HF, H2O and HCl, respectively), whereas CCSD overestimates these barriers by roughly the same magnitude. Scaling the barrier heights and dissociation energies by the number of fragments in the cluster reveals strong linear relationships between the two quantities and with the magnitudes of the imaginary vibrational frequency for the TSs.

Theoretical study of copper binding to GHK peptide

We report ligand field molecular mechanics, density functional theory and semi-empirical studies on the binding of Cu(II) to GlyHisLys (GHK) peptide. Following exhaustive conformational searching using molecular mechanics, we show that relative energy and geometry of conformations are in good agreement between GFN2-xTB semi-empirical and B3LYP-D DFT levels. Conventional molecular dynamics simulation of Cu-GHK shows the stability of the copper-peptide binding over 100 ps trajectory. Four equatorial bonds in 3N1O coordination remain stable throughout simulation, while a fifth in apical position from C-terminal carboxylate is more fluxional. We also show that the automated conformer and rotamer search algorithm CREST is able to correctly predict the metal binding position from a starting point consisting of separated peptide, copper and water.

Applications of hydrophobins: current state and perspectives

Hydrophobins are proteins exclusively produced by filamentous fungi. They self-assemble at hydrophilic-hydrophobic interfaces into an amphipathic film. This protein film renders hydrophobic surfaces of gas bubbles, liquids, or solid materials wettable, while hydrophilic surfaces can be turned hydrophobic. These properties, among others, make hydrophobins of interest for medical and technical applications. For instance, hydrophobins can be used to disperse hydrophobic materials; to stabilize foam in food products; and to immobilize enzymes, peptides, antibodies, cells, and anorganic molecules on surfaces. At the same time, they may be used to prevent binding of molecules. Furthermore, hydrophobins have therapeutic value as immunomodulators and can been used to produce recombinant proteins.

Assessing the importance of proton transfer reactions in DNA

Although engineered by millions of years of evolution, the cellular machinery is not flawless, and errors regularly appear during DNA replication. The subsequent alteration of the stored genetic message results in a mutation and might be the starting point of important health disorders. The question therefore is what causes DNA mutations? All living organisms are constantly exposed to a number of external agents such as free radicals and to radiation, which may lead to induced mutations. There are also mutations happening without invoking the action of any exogenous element, the so-called spontaneous mutations. The former can be partially controlled by avoiding exposure to high-risk environments, while the latter are more intriguing because their origin is unclear and difficult to determine. As noted by Watson and Crick when they first discovered the DNA structure, the correct replication of DNA rests on the assumption that the base pairs remain in their most stable, canonical form. However, protons along the interbase hydrogen-bond network are not static entities. They can in fact interchange their positions in DNA bases through proton transfer (PT) reactions before strands unwind, giving rise to noncanonical structures defined as rare tautomers. The importance of these rare tautomers was also cleverly anticipated by Watson and Crick and some years later claimed by Löwdin to be a source of spontaneous mutations. In Watson and Crick's words: "It would be of interest to know the precise difference in free energy between the various tautomeric forms under physiological conditions." Unfortunately, rare tautomeric forms are very difficult to detect, so no direct and accurate free energy measure has been discerned. In contrast, theoretical chemistry is making good progress toward the quantification of PT reactions in DNA and their biological consequences. This Account touches upon the theoretical studies devoted to appraising the importance of rare tautomers as promoters of spontaneous mutations. We focus in particular on the crucial role played by the biological environment on DNA stability. It has now been demonstrated that valuable macroscopic predictions require not only highly accurate theories but also refined chemical models. Hybrid quantum mechanics/molecular mechanics (QM/MM) simulations performed on short but complete DNA sequence fragments emerge in this context as the most adequate tools. In addition, these methods can be used to quantify the effect of different external agents on the PT tautomeric equilibria and, eventually, to conveniently handle them. This is the case for the possible alteration of the naturally observed mutation rate by exposure to intense electric fields. Theoretical predictions envision in this respect promising applications of ultrashort electric pulses in medicine to selectively modify the mutated/canonical ratio in DNA.

Theoretical studies on sulfur and metal cation (Cu(II), Ni(II), Pd(II), and Pt(II))-containing artificial DNA

We evaluated the stability, UV-vis spectra, and possibility of stacking of [S-M(II)-S] (M=Ni, Pd, Pt, S: hydroxypyridinethione) using density functional theory (DFT). We calculated the formation energies of modified bases with possible combinations of chalcogen atoms and metal cations. The results confirmed that [H-Ni(II)-H] (H: hydroxypyridone), [H-Cu(II)-H], and [S-Cu(II)-S] would form stable metal-base pairing; on the other hand, [H-Zn(II)-H], [S-Zn(II)-S], and [H-Pt(II)-H] would not. We predicted UV-vis excitations at 400-410 nm, mainly dominated by a d-pi* transition accompanied by a pi-pi* transition in [S-M(II)-S], where a metal-to-ligand charge transfer shifts the peak of S (without metal cations) by performing time-dependent density functional theory (TDDFT) calculations. By adding the Andersson-Langreth-Lundqvist van der Waals correction to the hybrid DFT results, we estimated the interaction energy between base pairs of [S-M(II)-S]. According to the results, the interaction between [S-M(II)-S] monomers for all cases is attractive, so that the selectivity of formation of the metal-containing DNA is governed mainly by the formation of the base pair in the aqueous phase and only slightly by the pi-pi stacking interaction between the modified bases.

Density functional theory for transition metals and transition metal chemistry

We introduce density functional theory and review recent progress in its application to transition metal chemistry. Topics covered include local, meta, hybrid, hybrid meta, and range-separated functionals, band theory, software, validation tests, and applications to spin states, magnetic exchange coupling, spectra, structure, reactivity, and catalysis, including molecules, clusters, nanoparticles, surfaces, and solids.

An exploration of mechanisms for the transformation of 8-oxoguanine to guanidinohydantoin and spiroiminodihydantoin by density functional theory

The potential energy surface for formation of 2-amino-5-hydroxy-7,9-dihydropurine-6,8-dione (5-OH-OG), guanidinohydantoin (Gh) and spiroiminodihydantoin (Sp) from 8-oxoguanine (8-oxoG) has been mapped out using B3LYP density functional theory, the aug-cc-pVTZ and 6-31+G(d,p) basis sets and the IEF-polarizable continuum model (PCM) solvation model. Three pathways for formation of 5-OH-OG from 8-oxoG were evaluated: (A) stepwise loss of two electrons and two protons to form the quinonoid intermediate 2-amino-7,9-dihydro-purine-6,8-dione (8-oxoG(ox)) followed by hydration; (B) stepwise loss of two electrons and one proton and net addition of hydroxide, in which the key step is nucleophilic addition to the 8-oxoG radical cation; and (C) stepwise loss of one electron and one proton and addition of hydroxyl radical to the 8-oxoG radical cation. The data suggest that all three pathways are energetically feasible mechanisms for the formation of 5-OH-OG, however, Pathway A may be kinetically favored over Pathway B. Although lower in energy, Pathway C may be of limited biological significance since it depends on the local concentration of hydroxyl radical. Pathways for hydrolysis and decarboxylation of 5-OH-OG to form Gh via either a carboxylic acid or substituted carbamic acid intermediate have been evaluated with the result that cleavage of the N1-C6 bond is clearly favored over that of the C5-C6 bond. Formation of Sp from 5-OH-OG via stepwise proton transfer and acyl migration or ring opening followed by proton transfer and ring closure have also been explored and suggest that deprotonation of the hydroxyl group facilitates a 1,2 acyl shift. Results of the calculations are consistent with experimental studies showing dependence of the Gh/Sp product ratio on pH. Under neutral and basic conditions, the data predict that formation of Sp is kinetically favored over the pathways for formation of Gh. Under acidic conditions, Gh is predicted to be the kinetically favored product.