José A. Gascón

Quantum Mechanics/Molecular Mechanics Study of the Catalytic Cycle of Water Splitting in Photosystem II

This paper investigates the mechanism of water splitting in photosystem II (PSII) as described by chemically sensible models of the oxygen-evolving complex (OEC) in the S0-S4 states. The reaction is the paradigm for engineering direct solar fuel production systems since it is driven by solar light and the catalyst involves inexpensive and abundant metals (calcium and manganese). Molecular models of the OEC Mn3CaO4Mn catalytic cluster are constructed by explicitly considering the perturbational influence of the surrounding protein environment according to state-of-the-art quantum mechanics/molecular mechanics (QM/MM) hybrid methods, in conjunction with the X-ray diffraction (XRD) structure of PSII from the cyanobacterium Thermosynechococcus elongatus. The resulting models are validated through direct comparisons with high-resolution extended X-ray absorption fine structure spectroscopic data. Structures of the S3, S4, and S0 states include an additional mu-oxo bridge between Mn(3) and Mn(4), not present in XRD structures, found to be essential for the deprotonation of substrate water molecules. The structures of reaction intermediates suggest a detailed mechanism of dioxygen evolution based on changes in oxidization and protonation states and structural rearrangements of the oxomanganese cluster and surrounding water molecules. The catalytic reaction is consistent with substrate water molecules coordinated as terminal ligands to Mn(4) and calcium and requires the formation of an oxyl radical by deprotonation of the substrate water molecule ligated to Mn(4) and the accumulation of four oxidizing equivalents. The oxyl radical is susceptible to nucleophilic attack by a substrate water molecule initially coordinated to calcium and activated by two basic species, including CP43-R357 and the mu-oxo bridge between Mn(3) and Mn(4). The reaction is concerted with water ligand exchange, swapping the activated water by a water molecule in the second coordination shell of calcium.

Interplay of Charge State, Lability, and Magnetism in the Molecule-like Au25(SR)18 Cluster

Au25(SR)18 (R = -CH2-CH2-Ph) is a molecule-like nanocluster displaying distinct electrochemical and optical features. Although it is often taken as an example of a particularly well-understood cluster, very recent literature has provided a quite unclear or even a controversial description of its properties. We prepared monodisperse Au25(SR)18(0) and studied by cyclic voltammetry, under particularly controlled conditions, the kinetics of its reduction or oxidation to a series of charge states, -2, -1, +1, +2, and +3. For each electrode process, we determined the standard heterogeneous electron-transfer (ET) rate constants and the reorganization energies. The latter points to a relatively large inner reorganization. Reduction to form Au25(SR)18(2-) and oxidation to form Au25(SR)18(2+) and Au25(SR)18(3+) are chemically irreversible. The corresponding decay rate constants and lifetimes are incompatible with interpretations of very recent literature reports. The problem of how ET affects the Au25 magnetism was addressed by comparing the continuous-wave electron paramagnetic resonance (cw-EPR) behaviors of radical Au25(SR)18(0) and its oxidation product, Au25(SR)18(+). As opposed to recent experimental and computational results, our study provides compelling evidence that the latter is a diamagnetic species. The DFT-computed optical absorption spectra and density of states of the -1, 0, and +1 charge states nicely reproduced the experimentally estimated dependence of the HOMO-LUMO energy gap on the actual charge carried by the cluster. The conclusions about the magnetism of the 0 and +1 charge states were also reproduced, stressing that the three HOMOs are not virtually degenerate as routinely assumed: In particular, the splitting of the HOMO manifold in the cation species is severe, suggesting that the usefulness of the superatom interpretation is limited. The electrochemical, EPR, and computational results thus provide a self-consistent picture of the properties of Au25(SR)18 as a function of its charge state and may furnish a methodology blueprint for understanding the redox and magnetic behaviors of similar molecule-like gold nanoclusters.

QM/MM Models of the O2-Evolving Complex of Photosystem II.

This paper introduces structural models of the oxygen-evolving complex of photosystem II (PSII) in the dark-stable S1 state, as well as in the reduced S0 and oxidized S2 states, with complete ligation of the metal-oxo cluster by amino acid residues, water, hydroxide, and chloride. The models are developed according to state-of-the-art quantum mechanics/molecular mechanics (QM/MM) hybrid methods, applied in conjunction with the X-ray crystal structure of PSII from the cyanobacterium Thermosynechococcus elongatus, recently reported at 3.5 Å resolution. Manganese and calcium ions are ligated consistently with standard coordination chemistry assumptions, supported by biochemical and spectroscopic data. Furthermore, the calcium-bound chloride ligand is found to be bound in a position consistent with pulsed electron paramagnetic resonance data obtained from acetate-substituted PSII. The ligation of protein ligands includes monodentate coordination of D1-D342, CP43-E354, and D1-D170 to Mn(1), Mn(3), and Mn(4), respectively; η(2) coordination of D1-E333 to both Mn(3) and Mn(2); and ligation of D1-E189 and D1-H332 to Mn(2). The resulting QM/MM structural models are consistent with available mechanistic data and also are compatible with X-ray diffraction models and extended X-ray absorption fine structure measurements of PSII. It is, therefore, conjectured that the proposed QM/MM models are particularly relevant to the development and validation of catalytic water-oxidation intermediates.

Au25(SEt)18, a Nearly Naked Thiolate-Protected Au25 Cluster: Structural Analysis by Single Crystal X-ray Crystallography and Electron Nuclear Double Resonance

X-ray crystallography has been fundamental in discovering fine structural features of ultrasmall gold clusters capped by thiolated ligands. For still unknown structures, however, new tools capable of providing relevant structural information are sought. We prepared a 25-gold atom nanocluster protected by the smallest ligand ever used, ethanethiol. This cluster displays the electrochemistry, mass spectrometry, and UV-vis absorption spectroscopy features of similar Au25 clusters protected by 18 thiolated ligands. The anionic and the neutral form of Au25(SEt)18 were fully characterized by (1)H and (13)C NMR spectroscopy, which confirmed the monolayers properties and the paramagnetism of neutral Au25(SEt)18(0). X-ray crystallography analysis of the latter provided the first known structure of a gold cluster protected by a simple, linear alkanethiolate. Here, we also report the direct observation by electron nuclear double resonance (ENDOR) of hyperfine interactions between a surface-delocalized unpaired electron and the gold atoms of a nanocluster. The advantages of knowing the exact molecular structure and having used such a small ligand allowed us to compare the experimental values of hyperfine couplings with DFT calculations unaffected by structures approximations or omissions.

Effect of the Charge State (z = −1, 0, +1) on the Nuclear Magnetic Resonance of Monodisperse Au25[S(CH2)2Ph]18z Clusters

Monodisperse Au(25)L(18)(0) (L = S(CH(2))(2)Ph) and [n-Oct(4)N(+)][Au(25)L(18)(-)] clusters were synthesized in tetrahydrofuran. An original strategy was then devised to oxidize them: in the presence of bis(pentafluorobenzoyl) peroxide, the neutral or the negatively charged clusters react as efficient electron donors in a dissociative electron-transfer (ET) process, in the former case yielding [Au(25)L(18)(+)][C(6)F(5)CO(2)(-)]. As opposed to other reported redox methods, this dissociative ET approach is irreversible, easily controllable, and clean, particularly for NMR purposes, as no hydrogen atoms are introduced. By using this approach, the -1, 0, and +1 charge states of Au(25)L(18) could be fully characterized by (1)H and (13)C NMR spectroscopy, using one- and two-dimensional techniques, in various solvents, and as a function of temperature. For all charge states, the NMR results and analysis nicely match recent structural findings about the presence of two different ligand populations in the capping monolayer, each resonance of the two ligand families displaying distinct NMR patterns. The radical nature of Au(25)L(18)(0) is particularly evident in the (1)H and (13)C NMR patterns of the inner ligands. The NMR behavior of radical Au(25)L(18)(0) was also simulated by DFT calculations, and the interplay between theory and experiments revealed a fundamental paramagnetic contribution coming from Fermi contact shifts. Interestingly, the NMR patterns of Au(25)L(18)(-) and Au(25)L(18)(+) were found to be quite similar, pointing to the latter cluster form as a diamagnetic species.

QM/MM approaches in medicinal chemistry research.

One of the goals of medicinal chemistry concerns the ability to compute protein-ligand interactions based on the structural knowledge of the receptor. To this end, the majority of current approaches incorporate classical force field potentials to describe receptor-ligand interactions. One of the most critical problems of standard molecular mechanics (MM) force fields is their fixed-charge treatment of electrostatic interactions. Two problems are derived from this approximation, polarization and charge transfer. As an immediate step in computational complexity, it seems natural to incorporate Quantum Mechanics (QM) within a hybrid QM/MM approach, which has shown to be a useful tool to describe structural and mechanistic aspects of chromophores and prosthetic residues in proteins. In this review, we describe specifically the role of QM/MM methods and their various applications to computational drug design and medicinal chemistry research in general.

QM/MM computational studies of substrate water binding to the oxygen-evolving centre of photosystem II

This paper reports computational studies of substrate water binding to the oxygen-evolving centre (OEC) of photosystem II (PSII), completely ligated by amino acid residues, water, hydroxide and chloride. The calculations are based on quantum mechanics/molecular mechanics hybrid models of the OEC of PSII, recently developed in conjunction with the X-ray crystal structure of PSII from the cyanobacterium Thermosynechococcus elongatus. The model OEC involves a cuboidal Mn3CaO4Mn metal cluster with three closely associated manganese ions linked to a single μ4-oxo-ligated Mn ion, often called the ‘dangling manganese’. Two water molecules bound to calcium and the dangling manganese are postulated to be substrate molecules, responsible for dioxygen formation. It is found that the energy barriers for the Mn(4)-bound water agree nicely with those of model complexes. However, the barriers for Ca-bound waters are substantially larger. Water binding is not simply correlated to the formal oxidation states of the metal centres but rather to their corresponding electrostatic potential atomic charges as modulated by charge-transfer interactions. The calculations of structural rearrangements during water exchange provide support for the experimental finding that the exchange rates with bulk 18O-labelled water should be smaller for water molecules coordinated to calcium than for water molecules attached to the dangling manganese. The models also predict that the S1→S2 transition should produce opposite effects on the two water-exchange rates.

A Self-Consistent Space-Domain Decomposition Method for QM/MM Computations of Protein Electrostatic Potentials

This paper introduces a self-consistent computational protocol for modeling protein electrostatic potentials according to static point-charge model distributions. The protocol involves a simple space-domain decomposition scheme where individual molecular domains are modeled as Quantum-Mechanical (QM) layers embedded in the otherwise classical Molecular-Mechanics (MM) protein environment. ElectroStatic-Potential (ESP) atomic charges of the constituent molecular domains are computed, to account for mutual polarization effects, and iterated until obtaining a self-consistent point-charge model of the protein electrostatic potential. The novel protocol achieves quantitative agreement with full QM calculations in the description of electrostatic potentials of small polypeptides where polarization effects are significant, showing a remarkable improvement relative to the corresponding electrostatic potentials obtained with popular MM force fields. The capabilities of the method are demonstrated in several applications, including calculations of the electrostatic potential in the potassium channel protein and the description of protein-protein electrostatic interactions.

QM/MM Study of the NMR Spectroscopy of the Retinyl Chromophore in Visual Rhodopsin

The (1)H and (13)C nuclear magnetic resonance (NMR) spectra of the retinyl chromophore in rhodopsin are investigated by using quantum mechanics/molecular mechanics (QM/MM) hybrid methods at the density functional theory (DFT) B3LYP/6-31G*:Amber level of theory, in conjunction with the gauge independent atomic orbital (GIAO) method for the ab initio self-consistent-field (SCF) calculation of NMR chemical shifts. The study provides a first-principle interpretation of solid-state NMR experiments based on recently developed QM/MM computational models of rhodopsin and bathorhodopsin [Gascón, J. A.; Batista, V. S. Biophys. J. 2004, 87, 2931-2941]. The reported results are particularly relevant to the development and validation of atomistic models of prototypical G-protein-coupled receptors which regulate signal transduction across plasma membranes.

Cardiac Glycoside Activities Link Na+/K+ ATPase Ion-Transport to Breast Cancer Cell Migration via Correlative SAR

The cardiac glycosides ouabain and digitoxin, established Na+/K+ ATPase inhibitors, were found to inhibit MDA-MB-231 breast cancer cell migration through an unbiased chemical genetics screen for cell motility. The Na+/K+ ATPase acts both as an ion-transporter and as a receptor for cardiac glycosides. To delineate which function is related to breast cancer cell migration, structure–activity relationship (SAR) profiles of cardiac glycosides were established at the cellular (cell migration inhibition), molecular (Na+/K+ ATPase inhibition), and atomic (computational docking) levels. The SAR of cardiac glycosides and their analogs revealed a similar profile, a decrease in potency when the parent cardiac glycoside structure was modified, for each activity investigated. Since assays were done at the cellular, molecular, and atomic levels, correlation of SAR profiles across these multiple assays established links between cellular activity and specific protein–small molecule interactions. The observed antimigratory effects in breast cancer cells are directly related to the inhibition of Na+/K+ transport. Specifically, the orientation of cardiac glycosides at the putative cation permeation path formed by transmembrane helices αM1–M6 correlates with the Na+ pump activity and cell migration. Other Na+/K+ ATPase inhibitors that are structurally distinct from cardiac glycosides also exhibit antimigratory activity, corroborating the conclusion that the antiport function of Na+/K+ ATPase and not the receptor function is important for supporting the motility of MDA-MB-231 breast cancer cells. Correlative SAR can establish new relationships between specific biochemical functions and higher-level cellular processes, particularly for proteins with multiple functions and small molecules with unknown or various modes of action.