Pablo Campomanes

Ligand substitutions between ruthenium–cymene compounds can control protein versus DNA targeting and anticancer activity

Ruthenium compounds have become promising alternatives to platinum drugs by displaying specific activities against different cancers and favourable toxicity and clearance properties. Nonetheless, their molecular targeting and mechanism of action are poorly understood. Here we study two prototypical ruthenium-arene agents—the cytotoxic antiprimary tumour compound [(η6-p-cymene)Ru(ethylene-diamine)Cl]PF6 and the relatively non-cytotoxic antimetastasis compound [(η6-p-cymene)Ru(1,3,5-triaza-7-phosphaadamantane)Cl2]—and discover that the former targets the DNA of chromatin, while the latter preferentially forms adducts on the histone proteins. Using a novel ‘atom-to-cell’ approach, we establish the basis for the surprisingly site-selective adduct formation behaviour and distinct cellular impact of these two chemically similar anticancer agents, which suggests that the cytotoxic effects arise largely from DNA lesions, whereas the protein adducts may be linked to the other therapeutic activities. Our study shows promise for developing new ruthenium drugs, via ligand-based modulation of DNA versus protein binding and thus cytotoxic potential, to target distinguishing epigenetic features of cancer cells.

Studies of glutathione transferase P1-1 bound to a platinum(IV)-based anticancer compound reveal the molecular basis of its activation.

Platinum-based cancer drugs, such as cisplatin, are highly effective chemotherapeutic agents used extensively for the treatment of solid tumors. However, their effectiveness is limited by drug resistance, which, in some cancers, has been associated with an overexpression of pi class glutathione S-transferase (GST P1-1), an important enzyme in the mercapturic acid detoxification pathway. Ethacraplatin (EA-CPT), a trans-Pt(IV) carboxylate complex containing ethacrynate ligands, was designed as a platinum cancer metallodrug that could also target cytosolic GST enzymes. We previously reported that EA-CPT was an excellent inhibitor of GST activity in live mammalian cells compared to either cisplatin or ethacrynic acid. In order to understand the nature of the drug-protein interactions between EA-CPT and GST P1-1, and to obtain mechanistic insights at a molecular level, structural and biochemical investigations were carried out, supported by molecular modeling analysis using quantum mechanical/molecular mechanical methods. The results suggest that EA-CPT preferentially docks at the dimer interface at GST P1-1 and subsequent interaction with the enzyme resulted in docking of the ethacrynate ligands at both active sites (in the H-sites), with the Pt moiety remaining bound at the dimer interface. The activation of the inhibitor by its target enzyme and covalent binding accounts for the strong and irreversible inhibition of enzymatic activity by the platinum complex.

Rhodopsin absorption from first principles: Bypassing common pitfalls

Bovine rhodopsin is the most extensively studied retinal protein and is considered the prototype of this important class of photosensitive biosystems involved in the process of vision. Many theoretical investigations have attempted to elucidate the role of the protein matrix in modulating the absorption of retinal chromophore in rhodopsin, but, while generally agreeing in predicting the correct location of the absorption maximum, they often reached contradicting conclusions on how the environment tunes the spectrum. To address this controversial issue, we combine here a thorough structural and dynamical characterization of rhodopsin with a careful validation of its excited-state properties via the use of a wide range of state-of-the-art quantum chemical approaches including various flavors of time-dependent density functional theory (TDDFT), different multireference perturbative schemes (CASPT2 and NEVPT2), and quantum Monte Carlo (QMC) methods. Through extensive quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations, we obtain a comprehensive structural description of the chromophore-protein system and sample a wide range of thermally accessible configurations. We show that, in order to obtain reliable excitation properties, it is crucial to employ a sufficient number of representative configurations of the system. In fact, the common use of a single, ad hoc structure can easily lead to an incorrect model and an agreement with experimental absorption spectra due to cancelation of errors. Finally, we show that, to properly account for polarization effects on the chromophore and to quench the large blue-shift induced by the counterion on the excitation energies, it is necessary to adopt an enhanced description of the protein environment as given by a large quantum region including as many as 250 atoms.

Ion Binding and Internal Hydration in the Multidrug Resistance Secondary Active Transporter NorM Investigated by Molecular Dynamics Simulations

Recently, a 3.65 Å resolution structure of the transporter NorM from the multidrug and toxic compound extrusion family has been determined in the outward-facing conformation. This antiporter uses electrochemical gradients to drive substrate export of a large class of antibiotic and toxic compounds in exchange for small monovalent cations (H(+) and Na(+)), but the molecular details of this mechanism are still largely unknown. Here we report all-atom molecular dynamics simulations of NorM, with and without the bound Na(+) cation and at different ion concentrations. Spontaneous binding of Na(+) is observed in several independent simulations with transient ion binding to D36 being necessary to reach the final binding site for which two competitive binding modes occur. Finally, the simulations indicate that the extracellular vestibule of the transporter invariably loses its characteristic V shape indicated by the crystallographic data, and it is reduced to a narrow permeation pathway lined by polar residues that can act as a specific pore for the transport of small cations. This event, together with the available structures of evolutionarily related transporters of the major facilitator superfamily (MFS), suggests that differences in the hydrophobic content of the extracellular vestibule may be characteristic of multidrug resistance transporters in contrast to substrate-selective members of the MFS.

Origin of the spectral shifts among the early intermediates of the rhodopsin photocycle.

A combined strategy based on the computation of absorption energies, using the ZINDO/S semiempirical method, for a statistically relevant number of thermally sampled configurations extracted from QM/MM trajectories is used to establish a one-to-one correspondence between the structures of the different early intermediates (dark, batho, BSI, lumi) involved in the initial steps of the rhodopsin photoactivation mechanism and their optical spectra. A systematic analysis of the results based on a correlation-based feature selection algorithm shows that the origin of the color shifts among these intermediates can be mainly ascribed to alterations in intrinsic properties of the chromophore structure, which are tuned by several residues located in the protein binding pocket. In addition to the expected electrostatic and dipolar effects caused by the charged residues (Glu113, Glu181) and to strong hydrogen bonding with Glu113, other interactions such as π-stacking with Ala117 and Thr118 backbone atoms, van der Waals contacts with Gly114 and Ala292, and CH/π weak interactions with Tyr268, Ala117, Thr118, and Ser186 side chains are found to make non-negligible contributions to the modulation of the color tuning among the different rhodopsin photointermediates.

Anandamide Hydrolysis in FAAH Reveals a Dual Strategy for Efficient Enzyme-Assisted Amide Bond Cleavage via Nitrogen Inversion

Herein, we combined classical molecular dynamics (MD) and quantum mechanical/molecular mechanics (QM/MM) simulations to unravel the whole catalytic cycle of fatty acid amide hydrolase (FAAH) in complex with anandamide, the main neurotransmitters involved in the control of pain. While microsecond MD simulations of FAAH in a realistic membrane/water environment provided a solid model for the reactant state of the enzymatic complex (Palermo et al. J. Chem. Theory Comput. 2013, 9, 1202-1213.), QM/MM simulations depict now a highly concerted two-step catalytic mechanism characterized by (1) acyl-enzyme formation after hydrolysis of the substrate amide bond and (2) deacylation reaction with restoration of the catalytic machinery. We found that a crucial event for anandamide hydrolysis is the inversion of the reactive nitrogen of the scissile amide bond, which occurs during the acylation rate-limiting step. We show that FAAH uses an exquisite catalytic strategy to induce amide bond distortion, reactive nitrogen inversion, and amide bond hydrolysis, promoting catalysis to completion. This new strategy is likely to be of general applicability to other amidases/peptidases that show similar catalytic site architectures, providing crucial insights for de novo enzyme design or drug discovery efforts.

Generalized QM/MM Force Matching Approach Applied to the 11-cis Protonated Schiff Base Chromophore of Rhodopsin

We extended a previously developed force matching approach to systems with covalent QM/MM boundaries and describe its user-friendly implementation in the publicly available software package CPMD. We applied this approach to the challenging case of the retinal protonated Schiff base in dark state bovine rhodopsin. We were able to develop a highly accurate force field that is able to capture subtle structural changes within the chromophore that have a pronounced influence on the optical properties. The optical absorption spectrum calculated from configurations extracted from a MD trajectory using the new force field is in excellent agreement with QM/MM and experimental references.

Keys to Lipid Selection in Fatty Acid Amide Hydrolase Catalysis: Structural Flexibility, Gating Residues and Multiple Binding Pockets

The fatty acid amide hydrolase (FAAH) regulates the endocannabinoid system cleaving primarily the lipid messenger anandamide. FAAH has been well characterized over the years and, importantly, it represents a promising drug target to treat several diseases, including inflammatory-related diseases and cancer. But its enzymatic mechanism for lipid selection to specifically hydrolyze anandamide, rather than similar bioactive lipids, remains elusive. Here, we clarify this mechanism in FAAH, examining the role of the dynamic paddle, which is formed by the gating residues Phe432 and Trp531 at the boundary between two cavities that form the FAAH catalytic site (the “membrane-access” and the “acyl chain-binding” pockets). We integrate microsecond-long MD simulations of wild type and double mutant model systems (Phe432Ala and Trp531Ala) of FAAH, embedded in a realistic membrane/water environment, with mutagenesis and kinetic experiments. We comparatively analyze three fatty acid substrates with different hydrolysis rates (anandamide > oleamide > palmitoylethanolamide). Our findings identify FAAH’s mechanism to selectively accommodate anandamide into a multi-pocket binding site, and to properly orient the substrate in pre-reactive conformations for efficient hydrolysis that is interceded by the dynamic paddle. Our findings therefore endorse a structural framework for a lipid selection mechanism mediated by structural flexibility and gating residues between multiple binding cavities, as found in FAAH. Based on the available structural data, this exquisite catalytic strategy for substrate specificity seems to be shared by other lipid-degrading enzymes with similar enzymatic architecture. The mechanistic insights for lipid selection might assist de-novo enzyme design or drug discovery efforts.

Assigning the EPR Fine Structure Parameters of the Mn(II) Centers in Bacillus subtilis Oxalate Decarboxylase by Site-Directed Mutagenesis and DFT/MM Calculations

Oxalate decarboxylase (OxDC) catalyzes the Mn-dependent conversion of the oxalate monoanion into CO2 and formate. EPR-based strategies for investigating the catalytic mechanism of decarboxylation are complicated by the difficulty of assigning the signals associated with the two Mn(II) centers located in the N- and C-terminal cupin domains of the enzyme. We now report a mutational strategy that has established the assignment of EPR fine structure parameters to each of these Mn(II) centers at pH 8.5. These experimental findings are also used to assess the performance of a multistep strategy for calculating the zero-field splitting parameters of protein-bound Mn(II) ions. Despite the known sensitivity of calculated D and E values to the computational approach, we demonstrate that good estimates of these parameters can be obtained using cluster models taken from carefully optimized DFT/MM structures. Overall, our results provide new insights into the strengths and limitations of theoretical methods for understanding electronic properties of protein-bound Mn(II) ions, thereby setting the stage for future EPR studies on the electronic properties of the Mn(II) centers in OxDC and site-specific variants.

Understanding regioselective cleavage in peptide hydrolysis by a palladium(II) aqua complex: a theoretical point of view.

Hydrolytic cleavage of the oligopeptides Ace-Ala-Lys-Tyr-Gly approximately Gly-Met-Ala-Ala-Arg-Ala and Ace-Lys-Gly-Gly-Ala-Gly approximately Pro-Met-Ala-Ala-Arg-Gly by [Pd(H(2)O)(4)](2+) was theoretically investigated by using molecular dynamics simulations and quantum mechanical calculations. The Pd anchorage to the peptide sequence is crucial to provoke the cleavage of the second bond upstream from the anchored methionine. For both cases, the most favorable reaction mechanism is a three-step route. The first step coincides with the experimental suggestion found for the Gly approximately Pro-Met sequence on a cleavage caused by an external attack of a water molecule to a complex in trans conformation of the scissile Gly approximately Gly and Gly approximately Pro peptide bonds. However, our results uncover the important role played by the presence of a Pd-coordinated water molecule, which simultaneously interacts with the carbonyl oxygen atom of the Gly amino acid in the Gly approximately Gly and Gly approximately Pro bonds. In accordance with experimental facts, the rise of the hydrolysis reaction rate when the Pro amino acid is located in the scissile peptide bond was also corroborated. The findings obtained at a molecular level from the present computations not only are relevant to rationalize the previously reported experiments but also could be of importance in designing new Pd(II) complexes for the regioselective cleavage of peptides and proteins.