A three-layer ONIOM model for the outside binding of cationic porphyrins and nucleotide pair DNA

In this work we investigated the outside binding mode between a cationic porphyrin and a nucleotide pair of DNA, adenine-thymine and guanine-cytosine, in a supramolecular assembly. We used two structural models (semi-extended, extended) that differ in the size of porphyrin, two kinds of theoretical methods: a three layer ONIOM (B3LYP/6-31G(d)/PM3/UFF), and DFT B3LYP/6-31G(d,p), and three cationic porphyrins. ONIOM method was first tested on the semi-extended model that was calculated in four environments: gas phase, solution phase using an explicit solvent model (H₂O), in the presence of a sodium cation (Na⁺) and in both (H₂O + Na⁺). From interaction energy results, we found that the affinity of the cationic substituent by the adenine nucleotide is favored upon the thymine nucleotide. The extended model that considers the whole porphyrin was applied in the gas phase to the four nucleotides. All the cationic porphyrins showed affinity by the nucleotides in the order adenine > guanine > thymine > cytosine. The interaction energy values for outside binding showed a strong porphyrin-nucleotide interaction (≈-90 kcal?mol^-1), that slightly varies between the nucleotides suggesting that this kind of cationic porphyrin has a little selectivity for some of them. We also found that the effect of the nature of the cationic substituent (chain length) in the porphyrin on the outside binding is small (≈2–13 kcal?mol^-1). Coherence between the results showed that ONIOM is a useful tool to get a reasonable molecular geometry to be used as a starting point in calculations of density functional theory.

Insight on the interaction of polychlorobiphenyl with nucleic acid–base

The interaction between one polychlorobiphenyl (3,3′,4,4′,-tetrachlorobiphenyl, coded PCB77) and the four DNA nucleic acid–base is studied by means of quantum mechanics calculations in stacked conformations. It is shown that even if the intermolecular dispersion energy is the largest component of the total interaction energy, some other contributions play a non negligible role. In particular the electrostatic dipole-dipole interaction and the charge transfer from the nucleobase to the PCB are responsible for the relative orientation of the monomers in the complexes. In addition, the charge transfer tends to flatten the PCB, which could therefore intercalate more easily between DNA base pairs. From these seminal results, we predict that PCB could intercalate completely between two base pairs, preferably between Guanine:Cytosine pairs.

Sequence selectivity of azinomycin B in DNA alkylation and cross-linking: a QM/MM study

Azinomycin B—a well-known antitumor drug—forms cross-links with DNA through alkylation of purine bases and blocks tumor cell growth. This reaction has been modeled using the ONIOM (B3LYP/6-31?+?g(d):UFF) method to understand the mechanism and sequence selectivity. ONIOM results have been checked for reliability by comparing them with full quantum mechanics calculations for selected paths. Calculations reveal that, among the purine bases, guanine is more reactive and is alkylated by aziridine ring through the C10 position, followed by alkylation of the epoxide ring through the C21 position of Azinomycin B. While the mono alkylation is controlled kinetically, bis-alkylation is controlled thermodynamically. Solvent effects were included using polarized-continuum-model calculations and no significant change from gas phase results was observed.

Synthesis and biological evaluation of cationic fullerene quinazolinone conjugates and their binding mode with modeled Mycobacterium tuberculosis hypoxanthine-guanine phosphoribosyltransferase enzyme

The present work reports a series of novel cationic fullerene derivatives bearing a substituted-quinazolinone moiety as a side arm. Fullerene-quinazolinone conjugates synthesized using the 1,3-dipolar cycloaddition reaction of C₆₀ with azomethine ylides generated from the corresponding Schiff bases of substituted quinazolinone were characterized by elemental analysis, FT-IR, ¹H NMR, ¹³C NMR and ESI-MS and screened for their antibacterial activity against Mycobacterium tuberculosis (H ₃₇ Rv strain). All the compounds exhibited significant activity with the most effective having MIC in the range of 1.562–3.125 μg/mL. Compound 9f exhibited good biological activity compared to standard drugs. We developed a computational strategy based on the modeled M. tuberculosis hypoxanthine-guanine phosphoribosyltransferase (HGPRT) using homology modeling techniques and studied its binding pattern with synthesized fullerene derivatives. We then explored the surface geometry of the protein to place the cage adjacent to the active site while optimizing its quinazolinone side arm to establish H bonding with active site residues.

A comparative DFT study on aquation and nucleobase binding of ruthenium (II) and osmium (II) arene complexes

The potential energy surfaces of the reactions of organometallic arene complexes of the type [(η ⁶-arene)M^II(pic)Cl] (where pic = 2-picolinic acid, M = Ru or Os) were examined by a DFT computational study. Among the seven density functional methods, hybrid exchange functional B3LYP outperforms the others to explain the aquation of the complexes. The reactions and binding energies of Ru^II and Os^II arene complexes with both 9EtG and 9EtA were studied to gain insight into the reactivity of these types of organometallic complexes with DNA. The obtained data rationalize experimental observation, contributing to partly understanding the potential biological and medical applications of organometallic complexes.

Penta- and heteropentadienyl ligands coordinated to beryllium

In this work we have performed a systematic study of new organometallic complexes containing penta- and heteropentadienyl (CH₂CHCHCHX, X?=?CH₂, O, NH, S) ligands coordinated to beryllium. Calculated complexes were studied using the density functional theory (PBE) in combination with the 6-311++G(3d,2p) basis set. The coordination number on the beryllium atom varies according to the type of ligand. Pentadienyl ligand shows hapticities η¹ and η⁵, while heteropentadienyl ligands display η¹ and η² hapticities. A Wiberg bond indices study was performed in order to get information about their bond orders.

Structural and energetic insights into sequence-specific interaction in DNA–drug recognition: development of affinity predictor and analysis of binding selectivity

Although the molecular mechanism and thermodynamic profile of a wide variety of chemical agents have been examined intensively in the past decades in terms of specific recognition of their protein receptors, to date the physicochemical nature of DNA–drug recognition and association still remains largely unexplored. The present study focused on understanding the structural basis, energetic landscape, and biological implications underlying the binding of small-molecule ligands to their cognate or non-cognate DNA receptors. First, a new method to capture the structural features of DNA–drug complex architecture was proposed and then used to correlate the extracted features with binding affinity of the complexes. By employing this method, a statistical regression-based predictor was developed to quantitatively evaluate the interaction potency of drug compounds with DNA in a fast and reliable manner. Subsequently, we use the predictor to examine the binding behavior of a number of structure-available, affinity-known DNA–drug complexes as well as a large pool of randomly generated DNA decoys in complex with the same drugs. It was found that (1) as compared with protein–DNA recognition, small-molecule agents have relatively low specificity in selecting their cognate DNA targets from the background of numerous random decoys; (2) the abundance of A–T base pairs in the DNA core motif exhibits a significant positive correlation with the affinity of drug ligand binding to the DNA receptor; and (3) high affinity seems not to be closely related to high selectivity for a DNA-targeting drug, although high-affinity drug entities have a greater possibility of being ranked computationally as top binders. We hope that this work will provide a preliminary insight into the molecular origin of sequence-specific interactions in DNA–drug recognition.

DFT studies of the conversion of four mesylate esters during reaction with ammonia

The energetics of the Menshutkin-like reaction between four mesylate derivatives and ammonia have been computed using B3LYP functional with the 6-31+G** basis set. Additionally, MPW1K/6-31+G** level calculations were carried out to estimate activation barrier heights in the gas phase. Solvent effect corrections were computed using PCM/B3LYP/6-31+G** level. The conversion of the reactant complexes into ion pairs is accompanied by a strong energy decrease in the gas phase and in all solvents. The ion pairs are stabilized with two strong hydrogen bonds in the gas phase. The bifurcation at C2 causes a significant activation barrier increase. Also, bifurcation at C5 leads to noticeable barrier height differentiation. Both B3LYP/6-31+G** and MPW1K/6-31+G** activation barriers suggest the reaction 2 (2a?+?NH₃) to be the fastest in the gas phase. The reaction 4 is the slowest one in all environments.

First-principles study of ammonium ions and their hydration in montmorillonites

Density functional theory calculations were performed to investigate the adsorption and hydration of an ammonium ion (NH₄ ⁺) confined in the interlayer space of montmorillonites (MMT). NH₄ ⁺ is trapped in the six-oxygen-ring on the internal surface and forms a strong binding with the surface O atoms. The hydration of NH₄ ⁺ is affected significantly by the surface. Water molecules prefer the surface sites, and do not bind with the NH₄ ⁺ unless enough water molecules are supplied. Moreover, the water molecules involved in NH₄ ⁺ hydration tend to bind with the surface simultaneously. The hydration energy increases with the intercalated water molecules, in contrast to that in gas phase. In addition, the hydration leads to the extension of MMT basal spacing.

Non-covalent interactions – QTAIM and NBO analysis

MP2(full)/6-311++G(3df,3pd) calculations were carried out on complexes linked through various non-covalent Lewis acid – Lewis base interactions. These are: hydrogen bond, dihydrogen bond, hydride bond and halogen bond. The quantum theory of ´atoms in molecules´ (QTAIM) as well as the natural bond orbitals (NBO) method were applied to analyze properties of these interactions. It was found that for the A-H…B hydrogen bond as well as for the A-X…B halogen bond (X designates halogen) the complex formation leads to the increase of s-character in the A-atom hybrid orbital aimed toward the H or X atom. In opposite, for the A…H-B hydride bond, where the H-atom possesses negative charge, the decrease of s-character in the B-atom orbital is observed. All these changes connected with the redistribution of the electron charge being the effect of the complex formation are in line with Bent´s rule. The numerous correlations between energetic, geometrical, NBO and QTAIM parameters were also found.

Modulating weak intramolecular interactions through the formation of beryllium bonds: complexes between squaric acid and BeH2

The electronic structure of the two most stable isomers of squaric acid and their complexes with BeH₂ were investigated at the B3LYP/6-311?+?G(3df,2p)// B3LYP/6-31?+?G(d,p) level of theory. Squaric acid forms rather strong beryllium bonds with BeH₂, with binding energies of the order of 60 kJ?mol^?1. The preferential sites for BeH₂ attachment are the carbonyl oxygen atoms, but the global minima of the potential energy surfaces of both EZ and ZZ isomers are extra-stabilized through the formation of a BeH···HO dihydrogen bond. More importantly, analysis of the electron density of these complexes shows the existence of significant cooperative effects between the beryllium bond and the dihydrogen bond, with both becoming significantly reinforced. The charge transfer involved in the formation of the beryllium bond induces a significant electron density redistribution within the squaric acid subunit, affecting not only the carbonyl group interacting with the BeH₂ moiety but significantly increasing the electron delocalization within the four membered ring. Accordingly the intrinsic properties of squaric acid become perturbed, as reflected in its ability to self-associate.

A new interaction mechanism of LiNH2 with MgH2: magnesium bond

Quantum chemical calculations were performed for LiNH₂–HMgX (X?=?H, F, Cl, Br, CH₃, OH, and NH₂) complexes to propose a new interaction mechanism between them. This theoretical survey showed that the complexes are stabilized through the combinative interaction of magnesium and lithium bonds. The binding energies are in the range of 63.2–66.5 kcal mol^?1, i.e., much larger than that of the lithium bond. Upon complexation, both Mg–H and Li–N bonds are lengthened. Substituents increase Mg-H bond elongation and at the same time decrease Li-N bond elongation. These cyclic complexes were characterized with the presence of a ring critical point and natural population analysis charges.

Prodrugs of fumarate esters for the treatment of psoriasis and multiple sclerosis—a computational approach

Density functional theory (DFT) calculations at B3LYP/6-31 G (d,p) and B3LYP/6-311?+?G(d,p) levels for the substituted pyridine-catalyzed isomerization of monomethyl maleate revealed that isomerization proceeds via four steps, with the rate-limiting step being proton transfer from the substituted pyridinium ion to the C=C double bond in INT1. In addition, it was found that the isomerization rate (maleate to fumarate) is solvent dependent. Polar solvents, such as water, tend to accelerate the isomerization rate, whereas apolar solvents, such as chloroform, act to slow down the reaction. A linear correlation was obtained between the isomerization activation energy and the dielectric constant of the solvent. Furthermore, linearity was achieved when the activation energy was plotted against the pK _a value of the catalyst. Substituted-pyridine derivatives with high pK _a values were able to catalyze isomerization more efficiently than those with low pK _a values. The calculated relative rates for prodrugs 1–6 were: 1 (406.7), 2 (7.6?×?10⁶), 3 (1.0), 4 (20.7), 5 (13.5) and 6 (2.2?×?10³). This result indicates that isomerizations of prodrugs 1 and 3–5 are expected to be slow and that of prodrugs 2 and 6 are expected to be relatively fast. Hence, prodrugs 2 and 3–5 have the potential to be utilized as prodrugs for the slow release of monomethylfumarate in the treatment of psoriasis and multiple sclerosis.

Insights into the structural determinants for selective inhibition of nitric oxide synthase isoforms

Selective inhibition of the nitric oxide synthase isoforms (NOS) is a promising approach for the treatment of various disorders. However, given the high active site conservation among all NOS isoforms, the design of selective inhibitors is a challenging task. Analysis of the X-ray crystal structures of the NOS isoforms complexed with known inhibitors most often gives no clues about the structural determinants behind the selective inhibition since the inhibitors share the same binding conformation. Aimed at a better understanding of the structural factors responsible for selective inhibition of NOS isoforms we have performed MD simulations for iNOS, nNOS and eNOS complexed with N^ω-NO₂-L-Arg (1), and with the aminopyridine derivatives 2 and 3. The slightly better selectivity of 1 for nNOS may be assigned to the presence of extra charge–charge interactions due to its “extended” conformation. While the high affinity of 2 for iNOS can be explained by the formation of an iNOS-specific subpocket upon binding, the lack of affinity for eNOS is associated to a conformational change in Glu363. The strong van der Waals and electrostatic interactions between 3 and the active site of nNOS are most likely responsible for its higher affinity for this isoform. Owing to the elongated and narrow binding pocket of iNOS, the correct positioning of 3 over the heme group is difficult, which may account for its lower affinity toward this isoform. Brought together, our results might help to rationalize the design of selective NOS inhibitors.

DFT and docking studies of rhodostreptomycins A and B and their interactions with solvated/nonsolvated Mg2+ and Ca2+ ions

The interactions of L-aminoglucosidic stereoisomers such as rhodostreptomycins A (Rho A) and B (Rho B) with cations (Mg²⁺, Ca²⁺, and H⁺) were studied by a quantum mechanical method that utilized DFT with B3LYP/6-311G**. Docking studies were also carried out in order to explore the surface recognition properties of L-aminoglucoside with respect to Mg²⁺ and Ca²⁺ ions under solvated and nonsolvated conditions. Although both of the stereoisomers possess similar physicochemical/antibiotic properties against Helicobacter pylori, the thermochemical values for these complexes showed that its high affinity for Mg²⁺ cations caused the hydration of Rho B. According to the results of the calculations, for Rho A–Ca²⁺(H₂O)₆, ΔH = ?72.21 kcal?mol^?1; for Rho B–Ca²⁺(H₂O)₆, ΔH = ?72.53 kcal?mol^?1; for Rho A–Mg²⁺(H₂O)₆, ΔH = ?72.99  kcal?mol^?1 and for Rho B–Mg²⁺(H₂O)₆, ΔH = ?95.00  kcal?mol^?1, confirming that Rho B binds most strongly with hydrated Mg²⁺, considering the energy associated with this binding process. This result suggests that Rho B forms a more stable complex than its isomer does with magnesium ion. Docking results show that both of these rhodostreptomycin molecules bind to solvated Ca²⁺ or Mg²⁺ through hydrogen bonding. Finally, Rho B is more stable than Rho A when protonation occurs.