Adsorption of water in Na-LTA zeolites: an ab initio molecular dynamics investigation.

The very wide range of applications of LTA zeolites, including the storage of tritiated water, implies that a detailed and accurate atomic-scale description of the adsorption processes taking place in their structure is crucial. To unravel with an unprecedented accuracy the mechanisms behind the water filling in NaA, we have conducted a systematic ab initio molecular dynamics investigation. Two LTA structural models, the conventional Z4A and the reduced one ZK4, have been used for static and dynamic ab initio calculations, respectively. After assessing this reduced model with comparative static DFT calculations, we start the filling of the α and ß cages by water, molecule by molecule. This allowed us to thoroughly study the interaction of water molecules with the zeolite structure and between water molecules, progressively forming H-bond chains and ring patterns as the cage is being filled. The adsorption energies could then be calculated with an unprecedented accuracy, which showed that the interaction of the molecules with the zeolite weakens as their number increases. By these methods, we have been able to highlight the primary role of Na+ cations in the interaction of water with zeolite, and inversely, the role of water in the displacement of cations when it is sufficiently solvated, allowing the passage between the α and ß cages. This phenomenon is possible thanks to the inhomogeneous distribution of water molecules on the cationic sites, as shown by our AIMD simulations, which allows the formation of water clusters. These results are important because they help in understanding how the coverage of cationic sites by water will affect the adsorption of other molecules inside the Na-LTA zeolite.

Theoretical insights into the mechanism of redox switch in heat shock protein Hsp33.

Heat shock protein 33 (Hsp33) is activated in the presence of H2O2 by a very interesting redox switch based on a tetra-coordinated zinc-cysteine complex present in the fully reduced and inactive protein form. The oxidation of this zinc center by H2O2 induces formation of two S-S bridges and the zinc release followed by the protein unfolding. We report here a theoretical study of the step-by-step sequence of the overall process starting with the oxidation of the first cysteine residue and ending with the zinc release. Each reaction step is characterized by its Gibbs free energy barrier (∆G (‡)). It is predicted that the first reaction step consists in the oxidation of Cys263 by H2O2 which is by far the most reactive cysteine (∆G (‡) = 15.4 kcal mol(-1)). The next two reaction steps are the formation of the first S-S bridge between Cys263 and Cys266 (∆G (‡) = 13.6 kcal mol(-1)) and the oxidation of Cys231 by H2O2 (∆G (‡) = 20.4 kcal mol(-1)). It is then shown that the formation of the second S-S bridge (Cys231-Cys233) before the zinc release is most unlikely (∆G (‡) = 34.8 kcal mol(-1)). Instead, the release of zinc just after the oxidation of the third cysteine (Cys231) is shown to be thermodynamically (dissociation Gibbs free energy ∆G d = 6.0 kcal mol(-1)) and kinetically (reaction rate constant k d ≈ 10(6) s(-1)) favored. This result is in good agreement with the experimental data on the oxidation mechanism of Hsp33 zinc center available to date.

A colorimetric assay adapted to fragment screening revealing aurones and chalcones as new arginase inhibitors.

Arginase, a difficult-to-target metalloenzyme, is implicated in a wide range of diseases, including cancer, infectious, and cardiovascular diseases. Despite the medical need, existing inhibitors have limited structural diversity, consisting predominantly of amino acids and their derivatives. The search for innovative arginase inhibitors has now extended to screening approaches. Due to the small and narrow active site of arginase, screening must meet the criteria of fragment-based screening. However, the limited binding capacity of fragments requires working at high concentrations, which increases the risk of interference and false positives. In this study, we investigated three colorimetric assays and selected one based on interference for screening under these challenging conditions. The subsequent adaptation and application to the screening a library of metal chelator fragments resulted in the identification of four compounds with moderate activity. The synthesis and evaluation of a series of compounds from one of the hits led to compound 21a with an IC50 value of 91.1 µM close to the reference compound piceatannol. Finally, molecular modelling supports the potential binding of aurones and chalcones to the active site of arginase, suggesting them as new candidates for the development of novel arginase inhibitors.

Hemoglobin as a nitrite anhydrase: modeling methemoglobin-mediated N2O3 formation.

Nitrite has recently been recognized as a storage form of NO in blood and as playing a key role in hypoxic vasodilation. The nitrite ion is readily reduced to NO by hemoglobin in red blood cells, which, as it happens, also presents a conundrum. Given NO's enormous affinity for ferrous heme, a key question concerns how it escapes capture by hemoglobin as it diffuses out of the red cells and to the endothelium, where vasodilation takes place. Dinitrogen trioxide (N(2)O(3)) has been proposed as a vehicle that transports NO to the endothelium, where it dissociates to NO and NO(2). Although N(2)O(3) formation might be readily explained by the reaction Hb-Fe(3+)+NO(2)(-)+NO⇌Hb-Fe(2+)+N(2)O(3), the exact manner in which methemoglobin (Hb-Fe(3+)), nitrite and NO interact with one another is unclear. Both an "Hb-Fe(3+)-NO(2)(-)+NO" pathway and an "Hb-Fe(3+)-NO+NO(2)(-) " pathway have been proposed. Neither pathway has been established experimentally. Nor has there been any attempt until now to theoretically model N(2)O(3) formation, the so-called nitrite anhydrase reaction. Both pathways have been examined here in a detailed density functional theory (DFT, B3LYP/TZP) study and both have been found to be feasible based on energetics criteria. Modeling the "Hb-Fe(3+)-NO(2)(-)+NO" pathway proved complex. Not only are multiple linkage-isomeric (N- and O-coordinated) structures conceivable for methemoglobin-nitrite, multiple isomeric forms are also possible for N(2)O(3) (the lowest-energy state has an N-N-bonded nitronitrosyl structure, O(2)N-NO). We considered multiple spin states of methemoglobin-nitrite as well as ferromagnetic and antiferromagnetic coupling of the Fe(3+) and NO spins. Together, the isomerism and spin variables result in a diabolically complex combinatorial space of reaction pathways. Fortunately, transition states could be successfully calculated for the vast majority of these reaction channels, both M(S)=0 and M(S)=1. For a six-coordinate Fe(3+)-O-nitrito starting geometry, which is plausible for methemoglobin-nitrite, we found that N(2)O(3) formation entails barriers of about 17-20 kcal mol(-1) , which is reasonable for a physiologically relevant reaction. For the "Hb-Fe(3+) -NO+NO(2) (-) " pathway, which was also found to be energetically reasonable, our calculations indicate a two-step mechanism. The first step involves transfer of an electron from NO(2)(-) to the Fe(3+)-heme-NO center ({FeNO}(6)) , resulting in formation of nitrogen dioxide and an Fe(2+)-heme-NO center ({FeNO}(7)). Subsequent formation of N(2)O(3) entails a barrier of only 8.1 kcal mol(-1) . From an energetics point of view, the nitrite anhydrase reaction thus is a reasonable proposition. Although it is tempting to interpret our results as favoring the "{FeNO}(6)+NO(2)(-) " pathway over the "Fe(3+)-nitrite+NO" pathway, both pathways should be considered energetically reasonable for a biological reaction and it seems inadvisable to favor a unique reaction channel based solely on quantum chemical modeling.

Negative Impact of Citral on Susceptibility of Pseudomonas aeruginosa to Antibiotics.

Essential oils (EOs) or their components are widely used by inhalation or nebulization to fight mild respiratory bacterial infections. However, their interaction with antibiotics is poorly known. In this study we evaluated the effects of citral, the main component of lemongrass oil, on in vitro susceptibility of Pseudomonas aeruginosa to antibiotics. Exposure of strain PA14 to subinhibitory concentrations of citral increased expression of operons encoding the multidrug efflux systems MexEF-OprN and MexXY/OprM, and bacterial resistance to anti-pseudomonal antibiotics including imipenem (twofold), gentamicin (eightfold), tobramycin (eightfold), ciprofloxacin (twofold), and colistin (≥128-fold). Use of pump deletion mutants showed that in addition to efflux other mechanisms were involved in this citral-induced phenotype. Determination of Zeta potential suggested that citral impairs the cell surface binding of aminoglycosides and colistin used at low concentrations (≤10 µg/mL). Moreover, experiments based on Raman spectroscopy and high-resolution mass spectrometry demonstrated formation of a Schiff base between the aldehyde group of citral and amino-groups of tobramycin and colistin. Chemical synthesis of tobracitryl, the imine compound resulting from condensation of citral and tobramycin, confirmed the loss of antibiotic activity due to adduct formation. Altogether these data point to the potential risk concern of self-medication with EOs containing citral in patients suffering from P. aeruginosa chronic lung infections and being treated with aerosols of aminoglycoside or colistin.

Cysteine oxidation by the superoxide radical: a theoretical study.

The cysteine residue oxidation by the superoxide radical in the gas phase and in aqueous solution is studied using the integrated molecular orbital+molecular orbital (IMOMO) method combining the quadratic configuration interaction [QCISD(T)] and density functional (DFT) methods. The molecular environment effects are systematically investigated by considering two alternative directions of attack of the superoxide radical on the thiol and two different cysteine residue conformations. It is found that hydrogen bonding and the electrostatic interactions between the superoxide radical and cysteine side chain significantly affect the reaction energy barrier, as compared to that derived for the simple thiol model methanethiol. Among the two possible reaction channels, the one involving the sulfinyl radical formation is predicted to be the dominant channel in aqueous solution. In a highly hydrophobic environment the thiyl radical formation channel becomes the main cysteine oxidation channel.

Fate of cisplatin and its main hydrolysed forms in the presence of thiolates: a comprehensive computational and experimental study.

Interaction of platinum-based drugs with proteins containing sulphur amino acids is usually argued as one of the major reasons for the observed resistance to these drugs, mainly due to the deactivation of the native compounds by very efficient thiolation processes in the organism. In this work, we have investigated the detailed thermodynamics and kinetics of reaction between cisplatin cis-[PtCl2(NH3)2] and its major hydrolysed forms (monohydroxocisplatin cis-[PtCl(OH)(NH3)2] and monoaquacisplatin cis-[PtCl(H2O)(NH3)2]+) with various thiolates (methanethiolate, cysteine and glutathione) and methionine. We have used a demanding quantum chemistry approach at the MP2 and DFT levels of theory to determine the Gibbs free energies and the barrier of reactions of the most possible reaction paths. The substitution of the four ligands of the complexes studied here (Cl-, OH-, H2O and NH3) can either proceed by direct thiolations or bidentations. Our Raman spectroscopy measurements show that only two thiolations actually occur, although four are possible in principle. The reason could lie in the bidentation reactions eventually taking place after each thiolation, which is backed up by our computational results. The observed lability scale of the ligands under thiolate exposure was found to be in the following order H2O > Cl- ≈ NH3(trans) > NH3(cis) > OH-, the difference between ammine ligands being induced by a significant trans-labilization by thiolates. Finally, the S,N bidentation is shown to be preferred with respect to the S,O one.

Methionine oxidation by hydrogen peroxide in peptides and proteins: A theoretical and Raman spectroscopy study.

The oxidation of proteins results in their deterioration via the oxidation of reactive amino acids. Oxidation of the amino acid, methionine plays an important role during biological conditions of oxidative stress, and equally a role in protein stability. In this study the oxidation of the methionine residue using the tripeptide GlyMetGly with respect to hydrogen peroxide has been studied using both Raman spectroscopy and DFT calculations. Spectral modifications following the formation of methionine sulfoxide are shown with the appearance of the SO vibration whilst there is also the modification of the CS vibrations at approximately 700 cm-1. The changes in the intensity of the CS stretching band were used to calculate the kinetic rate constant as 7.9 ±â€¯0.6 × 10-3 dm3 mol-1 s-1. The energy barrier for the reaction. is determined both experimentally and using DFT calculations. The reaction of the dairy protein beta-lactoglobulin with hydrogen peroxide is equally studied using the same technique. The solvent accessible surface area of the methionine residues within the protein were also determined and a comparison of the reaction rate constant and the energy barriers of reaction for the oxidation of the tripeptide and for the protein respectively thus, provides information about the role of the protein environment in the oxidation process.

Empirical force field for cisplatin based on quantum dynamics data: case study of new parameterization scheme for coordination compounds.

Parameterization of molecular complexes containing a metallic compound, such as cisplatin, is challenging due to the unconventional coordination nature of the bonds which involve platinum atoms. In this work, we develop a new methodology of parameterization for such compounds based on quantum dynamics (QD) calculations. We show that the coordination bonds and angles are more flexible than in normal covalent compounds. The influence of explicit solvent is also shown to be crucial to determine the flexibility of cisplatin in quantum dynamics simulations. Two empirical topologies of cisplatin were produced by fitting its atomic fluctuations against QD in vacuum and QD with explicit first solvation shell of water molecules respectively. A third topology built in a standard way from the static optimized structure was used for comparison. The later one leads to an excessively rigid molecule and exhibits much smaller fluctuations of the bonds and angles than QD reveals. It is shown that accounting for the high flexibility of cisplatin molecule is needed for adequate description of its first hydration shell. MD simulations with flexible QD-based topology also reveal a significant decrease of the barrier of passive diffusion of cisplatin accross the model lipid bilayer. These results confirm that flexibility of organometallic compounds is an important feature to be considered in classical molecular dynamics topologies. Proposed methodology based on QD simulations provides a systematic way of building such topologies.

An experimental and theoretical study of the amino acid side chain Raman bands in proteins.

The Raman spectra of a series of tripeptides with the basic formula GlyAAGly where the central amino acid (AA) was tryptophan, tyrosine, phenylalanine, glycine, methionine, histidine, lysine and leucine were measured in H2O. The theoretical Raman spectra obtained using density functional theory (DFT) calculations at the B3LYP/6-311+G(2df,2pd) level of theory allows a precise attribution of the vibrational bands. The experimental results show that there is a blue shift in the frequencies of several bands of the amino acid side chains in tripeptides compared to free amino acids, especially in the case of AAs containing aromatic rings. On the other hand, a very good agreement was found between the Raman bands of AA residues in tripeptides and those measured on three model proteins: bovine serum albumin, ß-lactoglobulin and lysozyme. The present analysis contributes to an unambiguous interpretation of the protein Raman spectra that is useful in monitoring the biological reactions involving AA side chains alteration.

Selenocysteine versus cysteine reactivity: a theoretical study of their oxidation by hydrogen peroxide.

The cysteine and selenocysteine oxidation by H2O2 in vacuo and in aqueous solution was studied using the integrated molecular orbital + molecular orbital (IMOMO) method combining the quadratic configuration method QCISD(T) and the spin projection of second-order perturbation theory PMP2. It is shown that including in the model system of cysteine (selenocysteine) residue up to 20 atoms has significant consequences upon the calculated reaction energy barrier. On the other hand, it is demonstrated that free cysteine and selenocysteine have very similar reaction energy barriers, 77-79 kJ mol(-1) in aqueous solution. It is thus concluded that the high experimental reaction rate constant reported for the oxidation of the selenocysteine residue in the glutathione peroxidase (GPx) active center is due to an important interaction between selenocysteine and its molecular environment. The sensitivity of the calculated energy barrier to the dielectric constant of the molecular environment observed for both cysteine and selenocysteine as well as the catalytic effect of the NH group emphasized in the case of cysteine supports this hypothesis.

Mechanism of thiol oxidation by the superoxide radical.

In spite of the large quantity of experimental work that deals with the oxidation of thiols by superoxide, the mechanism of this reaction is still controversial. The ab initio molecular orbital calculations reported here predict that the main reaction pathway includes the formation of a three-electron-bonded adduct followed by the elimination of the hydroxide anion, giving the sulfinyl radical as the reaction product. The alternative reaction pathway consisting of hydrogen atom transfer from the thiol to the protonated superoxide radical involves a reaction energy barrier that is significantly higher. The difference between the two reaction energy barriers is clearly beyond the expected computational uncertainty. The systematic scanning of the potential energy surface reveals no other competitive reaction pathways. The present results provide a useful basis for the interpretation of the complex experimental data related to thiol oxidation by superoxide radical in a biological environment.

Mechanism of cysteine oxidation by a hydroxyl radical: a theoretical study.

Cysteine oxidation by HO(.) was studied at a high level of ab initio theory in both gas phase and aqueous solution. Potential energy surface scans in the gas phase performed for the model system methanethiol+HO(.) indicate that the reactants can form two intermediate states: a sulfur-oxygen adduct and a hydrogen bound reactant complex. However these states appear to play a minor role in the reaction mechanism as long as they are fast dissociating states. Thus the main reaction channel predicted at the QCISD(T)/6-311+G(2df,2pd) level of theory is the direct hydrogen atom abstraction. The reaction mechanism is not perturbed by solvation which was found to induce only small variations in the Gibbs free energy of different reactant configurations. The larger size reactant system cysteine+HO* was treated by the integrated molecular orbital+molecular orbital (IMOMO) hybrid method mixing the QCISD(T)/6-311+G(2df,2pd) and the UMP2/6-311+G(d,p) levels of theory. The calculated potential energy, enthalpy, and Gibbs free energy barriers are slightly different from those of methanethiol. The method gave a rate constant for cysteine oxidation in aqueous solution, k=2.4 x 10(9) mol(-1) dm(3) s(-1), which is in good agreement with the experimental rate constant. Further analysis showed that the reaction is not very sensitive to hydrogen bonding and electrical polarity of the molecular environment.

A computational study of thiolate and selenolate oxidation by hydrogen peroxide.

Ab initio molecular orbital calculations have been used to study the effects of the molecular environment on the oxidation of thiolate and selenolate by hydrogen peroxide. The reaction was first examined in vacuo at the QCISD(T)/6-311+G(2df,2pd)//MP2/6-311+G(d,p) level of theory. It was found for both thiolate and selenolate that a reactant aggregate is formed, which has a dissociation rate constant comparable to the activation rate constant (about 10(-3) s(-1) for thiolate and 10(-1) s(-1) for selenolate). Using the polarizable continuum model (PCM) it was then found that the dissociation barrier energy decreases dramatically in water giving a dissociation rate constant of the order of 10(9) s(-1). In this case, the predicted overall rate constant of the thiolate reaction was about 10.2 mol(-1) dm3 s(-1), which is in good agreement with the experimental rate constant of cysteine oxidation in aqueous solution. The calculated rate constant for the selenolate reaction was somewhat higher (about 35.4 mol(-1) dm3 s(-1)). However, this value is several orders of magnitude smaller than the experimental value reported for the oxidation of selenocysteine in glutathione peroxidase. By considering the effect of the PCM dielectric constant on the reaction rate constant it was concluded that the high reactivity of the selenocysteine in glutathione peroxidase, as compared with cysteine, could be mainly due to the molecular environment of the selenocysteine residue.