The Ethidium–UA/AU Intercalation Site: Effect of Model Fragmentation and Backbone Charge State

Intercalation Ability of Novel Monofunctional Platinum Anticancer Drugs: A Key Step in Their Biological Action

Phenanthriplatin (PtPPH) is a monovalent platinum(II)-based complex with a large cytotoxicity against cancer cells. Although the aqua-activated drug has been assumed to be the precursor for DNA damage, it is still under debate whether the way in which that metallodrug attacks to DNA is dominated by a direct binding to a guanine base or rather by an intercalated intermediate product. Aiming to capture the mechanism of action of PtPPH, the present contribution used theoretical tools to systematically assess the sequence of all possible mechanisms on drug activation and reactivity, for example, hydrolysis, intercalation, and covalent damage to DNA. Ab initio quantum mechanical (QM) methods, hybrid QM/QM′ schemes, and independent gradient model approaches are implemented in an unbiased protocol. The performed simulations show that the cascade of reactions is articulated in three well-defined stages: (i) an early and fast intercalation of the complex between the DNA bases, (ii) a subsequent hydrolysis reaction that leads to the aqua-activated form, and (iii) a final formation of the covalent bond between PtPPH and DNA at a guanine site. The permanent damage to DNA is consequently driven by that latter bond to DNA but with a simultaneous π–π intercalation of the phenanthridine into nucleobases. The impact of the DNA sequence and the lateral backbone was also discussed to provide a more complete picture of the forces that anchor the drug into the double helix.

Importance of model size in quantum mechanical studies of DNA intercalation

The convergence of DFT-computed interaction energies with increasing binding site model size was assessed. The data show that while accurate intercalator interaction energies can be derived from binding site models featuring only the flanking nucleotides for uncharged intercalators that bind parallel to the DNA base pairs, errors remain significant even when including distant nucleotides for intercalators that are charged, exhibit groove-binding tails that engage in noncovalent interactions with distant nucleotides, or that bind perpendicular to the DNA base pairs. Consequently, binding site models that include at least three adjacent nucleotides are required to consistently predict converged binding energies. The computationally inexpensive HF-3c method is shown to provide reliable interaction energies and can be routinely applied to such large models.

Optimization of the linear-scaling local natural orbital CCSD(T) method: Redundancy-free triples correction using Laplace transform

An improved algorithm is presented for the evaluation of the (T) correction as a part of our local natural orbital (LNO) coupled-cluster singles and doubles with perturbative triples [LNO-CCSD(T)] scheme [Z. Rolik et al., J. Chem. Phys. 139, 094105 (2013)]. The new algorithm is an order of magnitude faster than our previous one and removes the bottleneck related to the calculation of the (T) contribution. First, a numerical Laplace transformed expression for the (T) fragment energy is introduced, which requires on average 3 to 4 times fewer floating point operations with negligible compromise in accuracy eliminating the redundancy among the evaluated triples amplitudes. Second, an additional speedup factor of 3 is achieved by the optimization of our canonical (T) algorithm, which is also executed in the local case. These developments can also be integrated into canonical as well as alternative fragmentation-based local CCSD(T) approaches with minor modifications. As it is demonstrated by our benchmark calculations, the evaluation of the new Laplace transformed (T) correction can always be performed if the preceding CCSD iterations are feasible, and the new scheme enables the computation of LNO-CCSD(T) correlation energies with at least triple-zeta quality basis sets for realistic three-dimensional molecules with more than 600 atoms and 12 000 basis functions in a matter of days on a single processor.

Physical Nature of Fatty Acid Amide Hydrolase Interactions with Its Inhibitors: Testing a Simple Nonempirical Scoring Model

Fatty acid amide hydrolase (FAAH) is an enzyme responsible for the deactivating hydrolysis of fatty acid ethanolamide neuromodulators. FAAH inhibitors have gained considerable interest due to their possible application in the treatment of anxiety, inflammation, and pain. In the context of inhibitor design, the availability of reliable computational tools for predicting binding affinity is still a challenging task, and it is now well understood that empirical scoring functions have several limitations that in principle could be overcome by quantum mechanics. Herein, systematic ab initio analyses of FAAH interactions with a series of inhibitors belonging to the class of the N-alkylcarbamic acid aryl esters have been performed. In contrast to our earlier studies of other classes of enzyme-inhibitor complexes, reasonable correlation with experimental results required us to consider correlation effects along with electrostatic term. Therefore, the simplest comprehensive nonempirical model allowing for qualitative predictions of binding affinities for FAAH ligands consists of electrostatic multipole and second-order dispersion terms. Such a model has been validated against the relative stabilities of the benchmark S66 set of biomolecular complexes. As it does not involve parameters fitted to experimentally derived data, this model offers a unique opportunity for generally applicable inhibitor design and virtual screening.

Atomistic account of structural and dynamical changes induced by small binders in the double helix of a short DNA

Nucleic acids are flexible molecules and their dynamical properties play a key role in molecular recognition events. Small binders interacting with DNA fragments induce both structural and dynamical changes in the double helix. We study the dynamics of a DNA dodecamer and of its complexes with Hoechst 33258, which is a minor groove binder, and with the ethidium cation, which is an intercalator, by molecular dynamics simulation. The thermodynamics of DNA-drug interaction is evaluated in connection with the structure and the dynamics of the resulting complexes. We identify and characterize the relevant changes in the configurational distribution of the DNA helix and relate them to the corresponding entropic contributions to the binding free energy. The binder Hoechst locks the breathing motion of the minor groove inducing a reduction of the configurational entropy of the helix, which amounts to 20 kcal mol(-1). In contrast, intercalations with the ethidium cation enhance the flexibility of the double helix. We show that the balance between the energy required to deform the helix for the intercalation and the gain in configurational entropy is the origin of cooperativity in the binding of a second ethidium and of anti-cooperativity in the binding of a third one. The results of our study provide an understanding of the relation between structure, dynamics and energetics in the interaction between DNA fragments and small binders, highlighting the role of dynamical changes and consequent variation of the configurational entropy of the DNA double helix for both types of binders.

A computational proposal for the experimentally observed discriminatory behavior of hypoxanthine, a weak universal nucleobase

A computational model composed of six nucleobases was used to investigate why hypoxanthine does not yield duplexes of equal stability when paired opposite each of the natural DNA nucleobases. The magnitudes of all nearest-neighbor interactions in a DNA helix were calculated, including hydrogen-bonding, intra- and interstrand stacking interactions, as well as 1-3 intrastrand stacking interactions. Although the stacking interactions in DNA relevant arrangements are significant and account for at least one third of the total stabilization energy in our nucleobase complexes, the trends in the magnitude of the stacking interactions cannot explain the relative experimental melting temperatures previously reported in the literature. Furthermore, although the total hydrogen-bonding interactions explain why hypoxanthine preferentially pairs with cytosine, the experimental trend for the remaining nucleobases (A, T, G) is not explained. In fact, the calculated pairing preference of hypoxanthine matches that determined experimentally only when the sum of all types of nearest-neighbor interactions is considered. This finding highlights a strong correlation between the relative magnitude of the total nucleobase-nucleobase interactions and measured melting temperatures for DNA strands containing hypoxanthine despite the potential role of other factors (including hydration, temperature, sugar-phosphate backbone). By considering a large range of sequence combinations, we reveal that the binding preference of hypoxanthine is strongly dependent on the nucleobase sequence, which may explain the varied ability of hypoxanthine to universally bind to the natural bases. As a result, we propose that future work should closely examine the interplay between the dominant nucleobase-nucleobase interactions and the overall strand stability to fully understand how sequence context affects the universal binding properties of modified bases and to aid the design of new molecules with ambiguous pairing properties.