RESUMEN

The aerobic oxidation of carbon monoxide to carbon dioxide is catalysed by the Mo/Cu-containing CO-dehydrogenase enzyme in the soil bacterium Oligotropha carboxidovorans, enabling the organism to grow on the small gas molecule as carbon and energy source. It was shown experimentally that silver can be substituted for copper in the active site of Mo/Cu CODH to yield a functional enzyme. In this study, we employed QM/MM calculations to investigate whether the reaction mechanism of the silver-substituted enzyme is similar to that of the native enzyme. Our results suggest that the Ag-substituted enzyme can oxidize CO and release CO2 following the same reaction steps as the native enzyme, with a computed rate-limiting step of 10.4 kcal/mol, consistent with experimental findings. Surprisingly, lower activation energies for C-O bond formation have been found in the presence of silver. Furthermore, comparison of rate constants for reduction of copper- and silver-containing enzymes suggests a discrepancy in the transition state stabilization upon silver substitution. We also evaluated the effects that differences in the water-active site interaction may exert on the overall energy profile of catalysis. Finally, the formation of a thiocarbonate intermediate along the catalytic pathway was found to be energetically unfavorable for the Ag-substituted enzyme. This finding aligns with the hypothesis proposed for the wild-type form, suggesting that the creation of such species may not be necessary for the enzymatic catalysis of CO oxidation.

Asunto(s)

RESUMEN

Affitins constitute a class of small proteins belonging to Sul7d family, which, in microorganisms such as Sulfolobus acidocaldarius, bind DNA preventing its denaturation. Thanks to their stability and small size (60-66 residues in length) they have been considered as ideal candidates for engineering and have been used for more than 10 years now, for different applications. The individuation of a mutant able to recognize a specific target does not imply the knowledge of the binding geometry between the two proteins. However, its identification is of undoubted importance but not always experimentally accessible. For this reason, computational approaches such as protein-protein docking can be helpful for an initial structural characterization of the complex. This method, which produces tens of putative binding geometries ordered according to a binding score, needs to be followed by a further reranking procedure for finding the most plausible one. In the present paper, we use the server ClusPro for generating docking models of affitins with different protein partners whose experimental structures are available in the Protein Data Bank. Then, we apply two protocols for reranking the docking models. The first one investigates their stability by means of Molecular Dynamics simulations; the second one, instead, compares the docking models with the interacting residues predicted by the Matrix of Local Coupling Energies method. Results show that the more efficient way to deal with the reranking problem is to consider the information given by the two protocols together, i.e. employing a consensus approach.

RESUMEN

The Suzuki-Miyaura reaction between the aryl halide (1) and the phenyl boronic acid (2), in the presence of the palladium(0) complex (3) as catalyst, gives the cross-coupling product (4) in quantitative yield when performed in basic aqueous solution of the nonionic surfactant Kolliphor-EL (K-EL). The partition between the aqueous and micellar environments of the species of this reaction has been investigated by means of Molecular Dynamics (MD) simulations. Starting from the K-EL molecules dispersed in water, a micelle model has been generated by MD simulations, adopting the 2016H66 force field. Reagent and product species have been described with the same force field, once the reliability of this force field has been tested comparing the n-octanol/water partition free energies calculated from the MD and Free Energy Perturbation (FEP) method with those obtained from the quantum-mechanical SMD method. The potential of mean force for the transfer process between water and the micellar phase of the different species has been calculated by the MD simulations and the Umbrella Sampling (US) method. The overall picture that emerges from these results confirms that the molecular species involved in this reaction prefers the micellar environment and concentrates in different but close zones of the micelle. This supports the experimental evidence that the use of suitable surfactant agents promotes reactivity, allowing micelles to behave as nanoreactors in which reactive species are solubilized and enhance their local concentration.

Asunto(s)

RESUMEN

Tailoring the activity of an organometallic catalyst usually requires a targeted ligand design. Tuning the ligand bulkiness and tuning the electronic properties are popular approaches, which are somehow interdependent because substituents of different sizes within ligands can determine inter alia the occurrence of different degrees of inductive effects. Ligand basicity, in particular, turned out to be a key property for the modulation of protonation reactions occurring in vacuo at the metals in complexes bearing organophosphorus ligands; however, when the same reactions take place in a polar organic solvent, their energetics becomes dependent on the trade-off between ligand basicity and bulkiness, with the polarity of the solvent playing a key role in this regard [Bancroft et al., Inorg. Chem., 1986, 25, 3675; Rovaletti et al., J. Phys. Org. Chem., 2018, 31, e3748]. In the present contribution, we carried out molecular dynamics and density functional theory calculations on water-soluble Mo-based catalysts for proton reduction, in order to study the energetics of protonation reactions in complexes where the incipient proton binds a catalytically active ligand (i.e., an oxide or a disulphide). We considered complexes either soaked in water or in a vacuum, and featuring N-based ancillary ligands of different bulkiness (i.e. cages constituted either by pyridine or isoquinoline moieties). Our results show that the energetics of protonation events can be affected by ancillary ligand bulkiness even when the metal center does not play the role of the H+ acceptor. In vacuo, protonation at the O or S atom in the α position relative to the metal in complexes featuring the bulky isoquinoline-based ligand is more favored by around 10 kcal mol-1 when compared to the case of the pyridine-based counterparts, a difference that is almost zero when the same reactions occur in water. Such an outcome is rationalized in light of the different electrostatic properties of complexes bearing ancillary ligands of different sizes. The overall picture from theory indicates that such effects of ligand bulkiness can be relevant for the design of green chemistry catalysts that undergo protonation steps in water solutions.

RESUMEN

The photon upconversion based on triplet-triplet annihilation (TTA) is a mechanism that converts the absorbed low-energy electromagnetic radiation into higher energy photons also at extremely low excitation intensities, but its use in actual technologies is still hindered by the limited availability of efficient annihilator moieties. We present here the results obtained by the synthesis and application of two new fluorinated chromophores based on phenazine and acridine structures, respectively. Both compounds show upconverted emission demonstrating their ability as TTA annihilator. More interesting, the acridine-based chromophore shows an excellent TTA yield that overcomes the one of some of best model systems. By correlating the experimental data and the quantum mechanical modeling of the investigated compound, we propose an alternative efficient pathway for the generation of the upconverted emissive states involving the peculiar high-energy triplet levels of the dye, thus suggesting a new development strategy for TTA annihilators based on the fine tuning of their high-energy excited states properties.

Asunto(s)

RESUMEN

The aerobic CO dehydrogenase from Oligotropha carboxidovorans is an environmentally crucial bacterial enzyme for maintenance of subtoxic concentration of CO in the lower atmosphere, as it allows for the oxidation of CO to CO2 which takes place at its Mo-Cu heterobimetallic active site. Despite extensive experimental and theoretical efforts, significant uncertainties still concern the reaction mechanism for the CO oxidation. In this work, we used the hybrid quantum mechanical/molecular mechanical approach to evaluate whether a water molecule present in the active site might act as a nucleophile upon formation of the new C-O bond, a hypothesis recently suggested in the literature. Our study shows that activation of H2 O can be favoured by the presence of the Mo=Oeq group. However, overall our results suggest that mechanisms other than the nucleophilic attack by Mo=Oeq to the activated carbon of the CO substrate are not likely to constitute reactive channels for the oxidation of CO by the enzyme.

Asunto(s)

RESUMEN

Intrinsically Disordered Peptides and Proteins (IDPs) in solution can span a broad range of conformations that often are hard to characterize by both experimental and computational methods. However, obtaining a significant representation of the conformational space is important to understand mechanisms underlying protein functions such as partner recognition. In this work, we investigated the behavior of the Sic1 Kinase-Inhibitor Domain (KID) in solution by Molecular Dynamics (MD) simulations. Our results point out that application of common descriptors of molecular shape such as Solvent Accessible Surface (SAS) area can lead to misleading outcomes. Instead, more appropriate molecular descriptors can be used to define 3D structures. In particular, we exploited Weighted Holistic Invariant Molecular (WHIM) descriptors to get a coarse-grained but accurate definition of the variegated Sic1 KID conformational ensemble. We found that Sic1 is able to form a variable amount of folded structures even in absence of partners. Among them, there were some conformations very close to the structure that Sic1 is supposed to assume in the binding with its physiological complexes. Therefore, our results support the hypothesis that this protein relies on the conformational selection mechanism to recognize the correct molecular partners.

RESUMEN

The chemical groups present at the surface of graphite have been thought for a long time to be mainly responsible for its catalytic activity in the oxygen reduction reaction. Recently, it was proposed that the surface defects of graphite also significantly contribute to promote this reaction. Although the behaviour of surface defects has been reported, only few comments have been dedicated to their involvement in the mechanism and the possible intermediate species in the oxygen reduction reaction. Herein, we aim to present a more detailed explanation of the catalytic activity of graphite particles based on the structure of their defects and their size. Structural, spectroscopic and magnetic investigation (X-ray diffraction, Raman and electron spin resonance) and electrochemical measurements were performed to describe the nature of the defects and their aptitude to transfer electrons. Computational description supplied precise details of the energy of the different defects and their ability to promote the reduction, also suggesting the structure of the intermediate adduct in the oxygen reduction. The results indicated that molecular oxygen preferentially interacts with graphite defects, which involve the π-electron system and accumulation of the spin density on the edges of the grains, in particular, on the zig-zag edges present on ball-milled graphite. This promotes the reactivity of this nanomaterial. Furthermore, the activation increases by decreasing the particle size.

RESUMEN

Carbon monoxide (CO) is a highly toxic gas to many living organisms. However, some microorganisms are able to use this molecule as the sole source of carbon and energy. Soil bacteria such as the aerobic Oligotropha carboxidovorans are responsible for the annual removal of about 2x108 tons of CO from the atmosphere. Detoxification through oxidation of CO to CO2 is enabled by the MoCu-dependent CO-dehydrogenase enzyme (MoCu-CODH) which-differently from other enzyme classes with similar function-retains its catalytic activity in the presence of atmospheric O2. In the last few years, targeted advancements have been described in the field of bioengineering and biomimetics, which is functional for future technological exploitation of the catalytic properties of MoCu-CODH and for the reproduction of its reactivity in synthetic complexes. Notably, a growing interest for the quantum chemical investigation of this enzyme has recently also emerged. This mini-review compiles the current knowledge of the MoCu-CODH catalytic cycle, with a specific focus on the outcomes of theoretical studies on this enzyme class. Rather controversial aspects from different theoretical studies will be highlighted, thus illustrating the challenges posed by this system as far as the application of density functional theory and hybrid quantum-classical methods are concerned.

RESUMEN

The Mo/Cu-dependent CO dehydrogenase from Oligotropha carboxidovorans is an enzyme that is able to catalyze CO oxidation to CO2; moreover, it can also oxidize H2, thus eliciting a characteristic EPR signal. Interestingly, the Ag-substituted enzyme form proved unable to catalyze H2 oxidation. In the present contribution, we characterized the reactivity of the enzyme with H2 by quantum-chemical calculations. It was found that dihydrogen binding to the wild-type enzyme requires significant structural rearrangements of the active site Theoretical EPR spectra for plausible H2-bound models of the partially reduced, paramagnetic active site are also presented and compared with the experimental counterpart. Finally, density functional theory modeling shows that Ag substitution impairs H2 binding at the active site.

Asunto(s)

RESUMEN

Ad hoc quantum chemical modeling of the acetyl-CoA synthase local structure and folding allowed us to identify an unprecedented coordination mode of histidine sidechain to protein-embedded metal ions.

Asunto(s)

RESUMEN

It is well known that transition metal ions are often bound to proteins, conveying very specific functional properties. In fact, metalloproteins play crucial biological roles in the transport and activation of small molecules such as H2, O2, and N2, as well as in several other biochemical processes. However, even if the presence of transition metals in the active site of proteins allows a very rich biochemistry, the experimental disclosure of structure-activity relationships in metalloproteins is generally difficult exactly because of the presence of transition metals, which are intrinsically characterized by a very versatile and often elusive chemistry. For this reason, computational methods are becoming very popular tools in the characterization of metalloproteins. In particular, since computing power is becoming less and less expensive, due to the continuous technological development of CPUs, the computational tools suited to investigate metalloproteins are becoming more accessible and therefore more commonly used also in molecular biology and biochemistry laboratories. Here, we present the main procedures and computational methods based on quantum mechanics, which are commonly used to study the structural, electronic, and reactivity properties of metalloproteins and related bioinspired compounds, with a specific focus on the practical and technical aspects that must be generally tackled to properly study such biomolecular systems.

Asunto(s)

RESUMEN

The adiabatic transition energies from the lowest triplet states of four Europium tris ß-diketonate/phenantroline complexes have been determined in vacuo and in dicholomethane solution by the ΔSCF approach at the density functional theory level, using the PBE1PBE and the CAM-B3LYP hybrid functionals. The calculated adiabatic transition energies have been compared with the experimental 0-0 transitions of each complex determined from phosphorescence spectra of the corresponding Gd(3+) complexes and followed by direct comparison between simulated and experimental spectra line shapes. For compound 1, the Eu(TTA)3Phen system, triplet states other than the lowest one and conformational isomers other than the one present in the crystallographic structure have been considered. In the crystallographic structure, this compound presents three quasi-degenerate low energy triplet states, differing for the TTA ligand where the two unpaired electrons are localized and showing close adiabatic transition energies. For compound 1, the lowest triplet states of the four investigated conformational isomers show similar characteristics and close adiabatic transition energies. On the basis of these results, an investigation of compounds 2-4 (Eu(Br-TTA)3Phen, Eu(DTDK)3Phen, and Eu(MeT-TTA)3) has been performed by considering only the isomer present in the crystallographic structure and only the lowest triplet state of each compound. For compounds 1-3, the energies of the lowest triplet states calculated by both functionals in solution including zero-point energy corrections well reproduce the experimental trends as well as the values of the adiabatic transition energies: CAM-B3LYP, the best performing functional, provides energies of the lowest triplet state with deviations from experiments lower than 1200 cm(-1). Also, the calculated vibrationally resolved phosphorescence spectra and UV-vis absorptions well reproduce the main features of their experimental counterparts. Significant differences between calculated and experimental results are observed for compound 4, for which difficulties in the experimental determination of the triplet state energy were encountered: our results show that the negligible photoluminescence quantum yield of this compound is due to the fact that the energy of the most stable triplet state is significantly lower than that of the resonance level of the Europium ion, and thus the energy transfer process is prevented. These results confirm the reliability of the adopted computational approach in calculating the energy of the lowest triplet state energy of these systems, a key parameter in the design of new ligands for lanthanide complexes presenting large photoluminescence quantum yields.

RESUMEN

[FeFe] hydrogenases are H2-evolving enzymes that feature a diiron cluster in their active site (the [2Fe]H cluster). One of the iron atoms has a vacant coordination site that directly interacts with H2, thus favoring its splitting in cooperation with the secondary amine group of a neighboring, flexible azadithiolate ligand. The vacant site is also the primary target of the inhibitor O2. The [2Fe]H cluster can span various redox states. The active-ready form (Hox) attains the Fe(II)Fe(I) state. States more oxidized than Hox were shown to be inactive and/or resistant to O2. In this work, we used density functional theory to evaluate whether azadithiolate-to-iron coordination is involved in oxidative inhibition and protection against O2, a hypothesis supported by recent results on biomimetic compounds. Our study shows that Fe-N(azadithiolate) bond formation is favored for an Fe(II)Fe(II) active-site model which disregards explicit treatment of the surrounding protein matrix, in line with the case of the corresponding Fe(II)Fe(II) synthetic system. However, the study of density functional theory models with explicit inclusion of the amino acid environment around the [2Fe]H cluster indicates that the protein matrix prevents the formation of such a bond. Our results suggest that mechanisms other than the binding of the azadithiolate nitrogen protect the active site from oxygen in the so-called H ox (inact) state.

Asunto(s)

RESUMEN

Magnetic Resonance Imaging (MRI) using paramagnetic systems as contrast agents is receiving increased attention as diagnostic tool in the clinic. At the same time, NMR of paramagnetic systems can also be applied in biochemical fields; for example, the use of Paramagnetic Relaxation Enhancement (PRE) allows structure refinement and the analysis of transient dynamic processes involved in macromolecular complex formation. Herein we report the synthesis and computational characterization of a new DOTA-like sialic acid conjugate, which can be used both in MRI and PRE applications when coordinated to a suitable paramagnetic metal.

Asunto(s)

RESUMEN

A new oxyiminopyrazole-based ytterbium chelate enables NIR emission upon UV excitation in colorless single layer luminescent solar concentrators for building integrated photovoltaics.

Asunto(s)

RESUMEN

The conformational behaviour in aqueous solution of the EgadMe complex, a conditional gadolinium-based contrast agent sensitive to beta-galactosidase enzymatic activity, is investigated by means of ab initio calculations and classical molecular dynamics simulations. Furthermore, force field parameterization of gadolinium-ligand interactions is performed, and its reliability is tested on the bench mark [Gd(DOTA)](-) system by MD simulations. Both computational methods highlight the presence in EgadMe of two main conformational isomers. The lowest energy conformation is a "close" form, corresponding to a state of low-relaxivity (MRI "inactive"), in which the ninth coordination site of the gadolinium ion is occupied by one oxygen atom of the galactopyranose residue. The second isomer, which is 2.9 (at ab initio level) and 4.2 (at MD level) kcal mol(-1) above the global minimum, presents an "open" form, corresponding to a state of high-relaxivity (MRI "active") in which one water molecule coordinates the ion. These results are consistent with experimental findings reported for EgadMe, and show that competition at the ninth coordination site of gadolinium ion, between the intra (the galactopyranose residue) and inter (water molecules) molecular interactions, affects the relaxivity of this system.

Asunto(s)

RESUMEN

The presence of amyloid is a hallmark of Gerstmann-Sträussler-Scheinker (GSS) disease, which is a prion disease caused by germ line mutations in the PRNP gene. The major component of amyloid is a fragment spanning residues from 81-82 to 144-153, part of the minimal sequence thought to play a crucial role in the conversion reaction and to sustain prion replication. We present here a molecular dynamics study on the 82-146 peptide from the human prion protein. The aim is to identify its aggregation-prone folds. The 82-146 prion sequence corresponds to a naturally occurring prion peptide able to form fibrils rich in parallel beta-sheets. A spontaneous right-handed beta-helical arrangement with 13 residues per turn can be observed in the 103-135 segment of the 82-146 peptide. The observed fold is in accordance with the evidence of a parallel beta-sheet organization in amyloid and with experiments on 82-146 discussed in the literature. To elucidate the conformational properties that trigger this peptide's aggregation propensity, the conformational behavior of peptides of different length (106-126 and 113-120 prion segments) was also investigated. Simulation analysis has led to some interesting considerations on sequence specific flexibility and the effects of growth. Comparing peptides of different length allows the localization of the origin of the beta-helix conformational propensity in the 106-126 segment, though longer sequences appear necessary for a clear beta-helical arrangement. Structural features of the observed 82-146 beta-helical fold are compatible with the "dock and lock" mechanism proposed to interpret peptide aggregation kinetics.

Asunto(s)

RESUMEN

There is evidence that Tetracyclines are potentially useful drugs to treat prion disease, the fatal neurodegenerative disease in which cellular prion proteins change in conformation to become a disease-specific species (PrP(Sc)). Based on an in vitro anti-fibrillogenesis test, and using the peptide PrP106-126 in the presence of tetracycline and 14 derivatives, we carried out a three-dimensional quantitative structure-activity relationship (3D-QSAR) study to investigate the stereoelectronic features required for anti-fibrillogenic activity. A preliminary variable reduction technique was used to search for grid points where statistical indexes of interaction potential distributions present local maximum (or minimum) values. Variable selection genetic algorithms were then used to search for the best 3D-QSAR models. A 6-variable model showed the best predictability of the anti-fibrillogenic activity that highlighted the best tetracycline substitution patterns: hydroxyl group presence in positions 5 and 6, electrodonor substituents on the aromatic D-ring, alkylamine substituent at the amidic group in position 2 and non-epi configuration of the NMe2 group.

Asunto(s)

RESUMEN

AIMS: The treatment of schizophrenic disorders is the most important challenge for community care. The analysis focuses on packages of care provided to 23.602 patients with a ICD-10 diagnosis of schizophrenic disorder and treated in 2001 by the Departments of Mental Health in Lombardy, Italy. METHODS: Packages of care refer to a mix of treatments provided to each patient during the year by different settings. Direct costs of the packages were calculated. Linear Discriminant Analysis has been used to link socio-demographic and diagnostic sub-groups of the patients to packages of care. RESULTS: People with schizophrenic disorders received relatively few care packages: only four packages involved more than 5%. Two thirds of the patients received only care provided by Community Mental Health Centres. In the other two packages with a percentage over 5%, the activity was provided by CMHCs, jointly with General Hospitals or Day Care Facilities. Complex care packages were rare (only 6%). As well as the intensity, also the variety of care provided by CMHCs increased with the complexity of care packages. In Lombardy more than half of the resources were spent for schizophrenia. The range of the costs per package was very wide. LDA failed to link characteristics of the patients to packages of care. CONCLUSIONS: Care packages are useful tools to understand better how mental health system works, how resources have been spent and to point out problems in the quality of care.

Asunto(s)