Reprint of PSII manganese cluster: protonation of W2, O5, O4 and His337 in the S1 state explored by combined quantum chemical and electrostatic energy computations.

Photosystem II (PSII) is a membrane-bound protein complex that oxidizes water to produce energized protons, which are used to built up a proton gradient across the thylakoidal membrane in the leafs of plants. This light-driven reaction is catalyzed by withdrawing electrons from the Mn4CaO5-cluster (Mn-cluster) in four discrete oxidation steps [S1-(S4/S0)] characterized in the Kok-cycle. In order to understand in detail the proton release events and the subsequent translocation of such energized protons, the protonation pattern of the Mn-cluster need to be elucidated. The new high-resolution PSII crystal structure from Umena, Kawakami, Shen, and Kamiya is an excellent basis to make progress in solving this problem. Following our previous work on oxidation and protonation states of the Mn-cluster, in this work, quantum chemical/electrostatic calculations were performed in order to estimate the pKa of different protons of relevant groups and atoms of the Mn-cluster such as W2, O4, O5 and His337. In broad agreement with previous experimental and theoretical work, our data suggest that W2 and His337 are likely to be in hydroxyl and neutral form, respectively, O5 and O4 to be unprotonated. This article is part of a special issue entitled: photosynthesis research for sustainability: keys to produce clean energy.

PSII manganese cluster: protonation of W2, O5, O4 and His337 in the S1 state explored by combined quantum chemical and electrostatic energy computations.

Photosystem II (PSII) is a membrane-bound protein complex that oxidizes water to produce energized protons, which are used to built up a proton gradient across the thylakoidal membrane in the leafs of plants. This light-driven reaction is catalyzed by withdrawing electrons from the Mn4CaO5-cluster (Mn-cluster) in four discrete oxidation steps [S1-(S4/S0)] characterized in the Kok-cycle. In order to understand in detail the proton release events and the subsequent translocation of such energized protons, the protonation pattern of the Mn-cluster need to be elucidated. The new high-resolution PSII crystal structure from Umena, Kawakami, Shen, and Kamiya is an excellent basis to make progress in solving this problem. Following our previous work on oxidation and protonation states of the Mn-cluster, in this work, quantum chemical/electrostatic calculations were performed in order to estimate the pKa of different protons of relevant groups and atoms of the Mn-cluster such as W2, O4, O5 and His337. In broad agreement with previous experimental and theoretical work, our data suggest that W2 and His337 are likely to be in hydroxyl and neutral form, respectively, O5 and O4 to be unprotonated. This article is part of a Special Issue entitled: Photosynthesis Research for Sustainability: Keys to Produce Clean Energy.

Anion⋠⋠⋠Si interactions in an inverse sandwich complex: a computational study.

Combination of an electron-rich molecule (e.g. chloride anion or nitrile group) with a chlorinated cyclohexasilane ring produces a supramolecular inverse sandwich complex formed by two guests (Cl(-) or RC≡N) strongly bonded to both faces of a planar host (Si6 ring). In-depth theoretical studies were carried out to investigate the nature of the bonding interactions that generate such a stable complex. Second-order Møller-Plesset perturbation theory (MP2) calculations confirmed that the presence of the Cl substituents is fundamental to the stability of the supramolecular assemblies. The density functional theory (DFT) functional wB97XD gave an estimation of the contribution of dispersion interactions to the binding energy. These interactions become more important as the Cl atoms of the rings are systematically replaced by methyl groups or hydrogen atoms. Analysis of the topology of the electron density and the reduced density gradient gave insight into the binding of the studied supramolecular assemblies.

Computational analysis of the nature and strength of the supramolecular contacts involved in the binding of chloride anions by imidazolium-based cyclic receptors.

The supramolecular bonding contacts driving the recognition of chloride ions by macrocyclic imidazolium-based receptors have been investigated by density functional theory calculations, both in vacuo and in solution (DMSO). This computational study reveals that the most stable host-guest complexes in vacuo and solution are different. While the anion interacts by means of two C-H···Cl(-) hydrogen bonds with the host molecule in vacuo (in a similar manner to that observed in the published single-crystal X-ray structure of the Cl-host complex), four C-H···Cl(-) contacts are clearly present in solution, as observed experimentally by earlier studies. In addition, the computed optimal Cl-host complex in solution confirms that the cavity of the host macrocycle, formed by four aromatic rings, does not includes the anions, which are located outside the cyclic receptor with which they interact through two CPh-H···Cl(-) hydrogen bonds and two unconventional (CIm-H)(+)···X(-) interactions.

Oxygen-evolving Mn cluster in photosystem II: the protonation pattern and oxidation state in the high-resolution crystal structure.

Extensive quantum chemical DFT calculations were performed on the high-resolution (1.9 Å) crystal structure of photosystem II in order to determine the protonation pattern and the oxidation states of the oxygen-evolving Mn cluster. First, our data suggest that the experimental structure is not in the S(1)-state. Second, a rather complete set of possible protonation patterns is studied, resulting in very few alternative protonation patterns whose relevance is discussed. Finally, we show that the experimental structure is a mixture of states containing highly reduced forms, with the largest contribution (almost 60%) from the S(-3)-state, Mn(II,II,III,III).

Heme proteins: the role of solvent in the dynamics of gates and portals.

Water plays a pivotal role in the correct functioning of proteins. Hydration is fundamental to their stability and flexibility, to folding process and specific functions, and to protein-protein interactions. In this work, the effects of solvation on proteins dynamics have been investigated by employing molecular dynamics simulations and using myoglobin as a model system. The investigation has been focused on solvent waters residing around/inside the protein, with average times of up to tens of nanoseconds, revealing that these slow waters may have significant effects on biological functioning of the protein. Our study pointed out that water is able to interact with proteins in diverse ways, leading to different kinds of perturbations in their intrinsic dynamic behavior. In particular, for myoglobin it was found that a water molecule can (i) "block" entry/escape of ligands to/from a particular docking site, (ii) act as a "wedge" modulating the dynamics of internal cavities, or (iii) join a "flow" of waters taking a ligand into (or "washing" a ligand away from) the protein interior. The information gathered in this work allowed us to provide a fingerprint of protein solvation state, the hydration sites map, which may represent a novel tool for comparing different forms/species of globular proteins.

Manganese(III)-mediated cyclodimerization of a hydrazinyl derivative generating an unprecedented 1,2,3,5,6-substituted leuco-verdazyl ring.

The one-pot reaction of Mn(OAc)(3) with two equivalents of 6-(pyridine-2-ylhydrazonomethyl)phenol (Hphp) in methanol leads to a unique leuco-verdazyl-type ring. The remarkable cyclodimerization reaction is proposed to involve a number of one-electron oxidative steps mediated by Mn(3+) ions, as well as a very uncommon 1,2-radical aryl migration and a nucleophilic attack of the solvent, i.e. methanol.

Copper-1,10-phenanthroline complexes binding to DNA: structural predictions from molecular simulations.

Copper-1,10-phenanthroline (phen) complexes Cu(phen)(2)(), Cu(2-Clip-phen), and Cu(3-Clip-phen) (Clip = a serinol bridge between the phen parts) are typically employed as DNA-cleaving agents and are now becoming increasingly important for building multifunctional drugs with improved cytotoxic properties. For instance, Cu(3-Clip-phen) has been combined with distamycin-like minor-groove binders and cisplatin-derivatives, leading to promising results. Density Functional Theory (DFT) and docking calculations as well as molecular dynamics (MD) simulations were performed to describe the mode of binding to DNA of these complexes. Our data suggest the minor-groove binding to be more probable than (partial) intercalation and major-groove binding. In addition, it was found that a combination of factors including planarity, van der Waals interactions with DNA, and structural complementarities may be the key for the cleavage efficiency of these copper complexes.

Robust recognition of malonate and 2-amino-4-picolinium in conjunction with M(II) as a triad (M = Ni/Co/Mn): role of this highly stable hydrogen-bonded motif in driving supramolecular self-assembly.

The crystal lattice of the three isostructural compounds , (C(6)H(8)N(2)H)(2)[M(C(3)H(2)O(4))(2)(H(2)O)(2)].4H(2)O (C(6)H(8)N(2)H = protonated 2-amino-4-picoline, M = Ni/Co/Mn, C(3)H(4)O(4) = malonate dianion; hereafter, malonate) is formed by supramolecular 2D layers. Hydrogen-bonding, pi...pi and lone pair...pi interactions play crucial role in organizing monomeric [M(II)(mal)(2)(H(2)O)(2)] units into 2D sheets along the ab plane, through the self-association between two different supramolecular building blocks, namely a tetrameric water cluster including metal-coordinated water molecules, and R(2)(2)(8) and R(2)(2)(7) hydrogen-bonded recognition synthons between 2-amino-4-picolinium and malonate. DFT calculations clearly show that the robust 2-amino-4-picolinium/malonate hydrogen-bonded motif drives the self-assembly of the supramolecular network observed.

Breathing motions of a respiratory protein revealed by molecular dynamics simulations.

Internal cavities, which are central to the biological functions of myoglobin, are exploited by gaseous ligands (e.g., O(2), NO, CO, etc.) to migrate inside the protein matrix. At present, it is not clear whether the ligand makes its own way inside the protein or instead the internal cavities are an intrinsic feature of myoglobin. To address this issue, standard molecular dynamics simulations were performed on horse-heart met-myoglobin with no ligand migrating inside the protein matrix. To reveal intrinsic internal pathways, the use of a statistical approach was applied to the cavity calculation, with special emphasis on the major pathway from the distal pocket to Xe1. Our study points out the remarkable dynamical behavior of Xe4, whose "breathing motions" may facilitate migration of ligands through the distal region. Additionally, our results highlight a two-way path for a ligand to diffuse through the proximal region, possibly allowing an alternative route in case Xe1 is occupied. Finally, our approach has led us to the identification of key residues, such as leucines, that may work as switches between cavities.

Water-chloride and water-bromide hydrogen-bonded networks: influence of the nature of the halide ions on the stability of the supramolecular assemblies.

Two compounds, namely, [TTPH(2)](Cl)(2) x 4 H(2)O (1) and [TTPH(2)](Br)(2) x 4 H(2)O (2), (TTP = 4'-p-tolyl-2,2':6',2''-terpyridine) were synthesized from purely aqueous media and characterized by physical techniques. In the solid-state structures of these compounds, interesting supramolecular assemblies are observed. In 1, an unusual staircase-like architecture of the tape of edge-shared planar water hexamer is of importance, where the chloride ions are at the two edges of the tape. In 2, the polymeric nature of the water-bromide assembly is of interest, where discrete open-cube water octamers are doubly bridged by bromide ions. Semiempirical and DFT calculations confirm that the nature of the anion indeed affects the topology of the water-halide assemblies. We conclude that the protonated [TTPH(2)](2+) species can act as appropriate receptors for halide ions, which in turn act as a matrix for the formation of polymeric 1D water-halide assemblies.

Magnetic coupling between copper(II) ions mediated by hydrogen-bonded (neutral) water molecules.

A new hydrogen-bonded dinuclear copper(II) coordination compound has been synthesized from the Schiff-base ligand 6-(pyridine-2-ylhydrazonomethyl)phenol (Hphp). The molecular structure of [Cu(2)(php)(2)(H(2)O)(2)(ClO(4))](ClO(4))(H(2)O) (1), determined by single-crystal X-ray diffraction, reveals the presence of two copper(II) centers held together by means of two strong hydrogen bonds, with O...O contacts of only 2.60-2.68 A. Temperature-dependent magnetic susceptibility measurements down to 3 K show that the two metal ions are antiferromagnetically coupled (J = -19.8(2) cm(-1)). This exchange is most likely through two hydrogen-bonding pathways, where a coordinated water on the first Cu, donates a H bond to the O atoms of the coordinated php at the other Cu. This strong O...H(water) bonding interaction has been clearly evidenced by theoretical calculations. In the relatively few related cases from the literature, this exchange path, mediated by a (neutral) coordinated water molecule, was not recognized.

Supramolecular lone pair-pi/pi-pi/pi-anion assembly in a Mg(II)-malonate-2-aminopyridine-nitrate ternary system.

The solid-state structure of an alkaline-earth metal complex reveals the formation of a remarkable supramolecular framework based on concurrent lone pair-pi, pi-pi, and pi-anion interactions whose stability has been investigated by density functional theory.

Concurrent anion...pi interactions between a perchlorate ion and two pi-acidic aromatic rings, namely pentafluorophenol and 1,3,5-triazine.

The rational design of a ligand containing two electron-poor pi-rings, i.e. a triazine and a pentafluorophenoxy groups, has allowed the preparation of a copper complex where both the anticipated anion...pi interactions are present.

The hydrolysis mechanism of the anticancer ruthenium drugs NAMI-A and ICR investigated by DFT-PCM calculations.

(ImH)[trans-RuCl(4)(DMSO-S)(Im)], (Im = imidazole, DMSO-S = S-bonded dimethylsulfoxide), NAMI-A, is the first anticancer ruthenium compound that successfully completed Phase I clinical trials. NAMI-A shows a remarkable activity against lung metastases of solid tumors, but is not effective in the reduction of primary cancer. The structurally similar (ImH)[trans-RuCl(4)(Im)(2)], ICR (or KP418), and its indazole analog (KP1019) are promising candidate drugs in the treatment of colorectal cancers, but have no antimetastatic activity. Despite the pharmacological relevance of these compounds, no rationale has been furnished to explain their markedly different activity. While the nature of the chemical species responsible for their antimetastatic/anticancer activity has not been determined, it has been suggested that the difference between reduction potentials of NAMI-A and ICR may be the key to the different biological responses they induce. In this work, Density Functional Theory calculations were performed to investigate the hydrolysis of NAMI-A and ICR in both Ru(III) and Ru(II) oxidation states, up to the third aquation. In line with experimental findings, our calculations provide a picture of the hydrolysis of NAMI-A and ICR mainly as a stepwise loss of chloride ligands. While dissociation of Im is unlikely under neutral conditions, that of DMSO becomes competitive with the loss of chloride ions as the hydrolysis proceeds. Redox properties of NAMI-A and ICR and of their most relevant hydrolytic intermediates were also studied in order to monitor the effects of biological reductants on the mechanism of action. Our findings may contribute to the identification of the active compounds that interact with biological targets, and to explain the different biological activity of NAMI-A and ICR.

Density functional theory studies on copper phenanthroline complexes.

Density functional theory calculations have been employed to investigate the role of structural properties of copper phenanthroline complexes for DNA-cleavage activity. Structural changes imposed on the coordination geometries of Cu(phen)(2)(+,2+) (phen = 1,10-phenanthroline) linked by a serinol bridge (abbreviated as Clip) were studied, as well as their energetic profiles. Our calculations show that structures of these copper complexes (in this work named as clipped complexes) strongly depend on the position of the link, rather than on the copper oxidation state. Ionization energies slightly differ among the three selected complexes, while inner-sphere reorganization energies more markedly depend on the serinol link. However, the relative rates of the redox reaction of Cu(phen)(2), Cu(2-Clip-phen), and Cu(3-Clip-phen) were found not to correlate with their relative DNA-cleavage activity experimentally observed. Thus, the serinol link mainly affects the structural properties of copper phenanthroline complexes rather than their electronic properties. Docking simulations of clipped and nonclipped Cu(I) phenanthroline complexes on a DNA 16mer, d[CGCTCAACTGTGATAC](2), were finally performed to assess how different structural properties could affect the formation of DNA adducts. This analysis revealed that the most stable adducts of Cu(phen)(2+) and Cu(3-Clip-phen)(+) with DNA bind in the minor groove, whereas Cu(2-Clip-phen)(+) binds preferentially into the major groove.

Nitrogen Fixation by a Molybdenum Catalyst Mimicking the Function of the Nitrogenase Enzyme: A Critical Evaluation of DFT and Solvent Effects.

Compounds mimicking the enzyme nitrogenase represent promising alternative routes to the current Haber-Bosch industrial synthesis of ammonia from molecular hydrogen and nitrogen. In this work, we investigated the full catalytic cycle of one of such compounds, Mo(HIPTN3N) (with HIPT = hexaisopropylterphenyl), by means of DFT calculations. Our results suggest these large ligands to exert mainly a steric influence on the structural properties of the catalyst. In addition, we provided a structural and electronic characterization of the putative reaction intermediates along with a picture of the electronic mechanism of molecular nitrogen N-N bond breaking. A large discrepancy was observed between calculated and experimental reaction free energies, suggesting that in the present case the predictability of DFT reaction energies is limited. Investigation of explicit solvation of specific catalytic intermediates as well as of the protonation and reducing agents reveal the crucial role played by the solvent molecules (benzene and heptane) particularly for protonation steps. Furthermore, the analysis of several DFT functionals indicates that these have to be carefully chosen in order to reproduce the energetic profile of reduction steps. This study shows how DFT calculations may be a powerful tool to describe structural and electronic properties of the intermediates of the catalytic cycle, yet, due to the complexity of the system, reaction energies cannot be easily reproduced without a careful choice of the solvation model and the exchange-correlation functional.

A QM/MM study of cisplatin-DNA oligonucleotides: from simple models to realistic systems.

QM/MM calculations were employed to investigate the role of hydrogen bonding and pi stacking in several single- and double-stranded cisplatin-DNA structures. Computed geometrical parameters reproduce experimental structures of cisplatin and its complex with guanine-phosphate-guanine. Following QM/MM optimisation, single-point DFT calculations allowed estimation of intermolecular forces through atoms in molecules (AIM) analysis. Binding energies of platinated single-strand DNA qualitatively agree with myriad experimental and theoretical studies showing that complexes of guanine are stronger than those of adenine. The topology of all studied complexes confirms that platination strongly affects the stability of both single- and double-stranded DNAs: Pt-N-H...X (X = N or O) interactions are ubiquitous in these complexes and account for over 70 % of all H-bonding interactions. The pi stacking is greatly reduced by both mono- and bifunctional complexation: the former causes a loss of about 3-4 kcal mol(-1), whereas the latter leads to more drastic disruption. The effect of platination on Watson-Crick GC is similar to that found in previous studies: major redistribution of energy occurs, but the overall stability is barely affected. The BH&H/AMBER/AIM approach was also used to study platination of a double-stranded DNA octamer d(CCTG*G*TCC)d(GGACCAGG), for which an experimental structure is available. Comparison between theory and experiment is satisfactory, and also reproduces previous DFT-based studies of analogous structures. The effect of platination is similar to that seen in model systems, although the effect on GC pairing was more pronounced. These calculations also reveal weaker, secondary interactions of the form Pt...O and Pt...N, detected in several single- and double-stranded DNA.

Gas-phase DNA oligonucleotide structures. A QM/MM and atoms in molecules study.

QM/MM calculations have been employed to investigate the role of hydrogen bonding and pi-stacking in single- and double-stranded DNA oligonucleotides. DFT calculations and Atoms in Molecules analysis on QM/MM-optimized structures allow characterization and estimation of the energies of pi-stacking and hydrogen-bond interactions. This shows that pi-stacking interactions depend on the number and the nature of the DNA bases for single-stranded nucleotides; for instance, guanines are found to be involved in strong hydrogen bonds, whereas adenines interact mainly via stacking interactions. The role of interbase hydrogen bonding was explored: the -NH2 groups of guanine, adenine, and cytosine participate in N-H...O and N-H...N interactions. These are much stronger in single-strand oligonucleotides, where the -NH2 groups are highly nonplanar. In double-stranded DNA, the strong base-pairing hydrogen bonds of complementary bases lead to more planar -NH2 groups, which tend to be involved in pi-stacking interactions rather than H-bonds. The use of AIM also allows us to evaluate the interplay of pi-stacking and H-bonding, suggesting that cooperativity does occur, but is generally limited to about 1-2 kcal/mol.

Hybrid density functional theory for pi-stacking interactions: application to benzenes, pyridines, and DNA bases.

The suitability of a hybrid density functional to qualitatively reproduce geometric and energetic details of parallel pi-stacked aromatic complexes is presented. The hybrid functional includes an ad hoc mixture of half the exact (HF) exchange with half of the uniform electron gas exchange, plus Lee, Yang, and Parr's expression for correlation energy. This functional, in combination with polarized, diffuse basis sets, gives a binding energy for the parallel-displaced benzene dimer in good agreement with the best available high-level calculations reported in the literature, and qualitatively reproduces the local MP2 potential energy surface of the parallel-displaced benzene dimer. This method was further critically compared to high-level calculations recently reported in the literature for a range of pi-stacked complexes, including monosubstituted benzene-benzene dimers, along with DNA and RNA bases, and generally agrees with MP2 and/or CCSD(T) results to within +/-2 kJ mol(-1). We also show that the resulting BH&H binding energy is closely related to the electron density in the intermolecular region. The net result is that the BH&H functional, presumably due to fortuitous cancellation of errors, provides a pragmatic, computationally efficient quantum mechanical tool for the study of large pi-stacked systems such as DNA.