Publicacions

Simulations based on molecular dynamics coupled to excitation energy calculations were used to generate simulated absorption spectra for a family of halide derivatives of azobenzene, a family of photoswitch molecules with a weak absorption band around 400–600 nm and potential uses in living tissue. This is a case where using the conventional approach in theoretical spectroscopy (estimation of absorption maxima based on the vertical transition from the potential energy minimum on the ground electronic state) does not provide valid results that explain how the observed band shape extends towards the low energy region of the spectrum. The method affords a reasonable description of the main features of the low-energy UV-Vis spectra of these compounds. A bathochromic trend was detected linked to the size of the halide atom. Analysis of the excitation reveals a correlation between the energy of the molecular orbital where excitation starts and the energy of the highest occupied atomic orbital of the free halide atom. This was put to the test with a new brominated compound with good results. The energy level of the highest occupied orbital on the free halide was identified as a key factor that strongly affects the energy gap in the photoswitch. This opens the way for the design of bathochromically shifted variants of the photoswitch with possible applications.

In order to theoretically design multi-state photoswitches with specific properties, an exhaustive computational study is first carried out for an azobenzene dimer that has been recently synthesized and experimentally studied. This study allows for a full comprehension of the factors that govern the photoactivated isomerization processes of these molecules so to provide a conceptual/computational protocol that can be applied to generic multi-state photoswitches. From this knowledge a new dimer with a similar chemical design is designed and also fully characterized. Our theoretical calculations predict that the new dimer proposed is one step further in the quest for a double photoswitch, where the four metastable isomers could be selectively interconverted through the use of different irradiation sequences.

The 3-phosphoinositide-dependent protein kinase 1 (PDK1) K465E mutant kinase can still activate protein kinase B (PKB) at the membrane in a phosphatidylinositol-3,4,5-trisphosphate (PIP₃, PtdIns(3,4,5)P₃) independent manner. To understand this new PDK1 regulatory mechanism, docking and molecular dynamics calculations were performed for the first time to simulate the wild-type kinase domain–pleckstrin homology (PH) domain complex with PH-in and PH-out conformations. These simulations were then compared to the PH-in model of the KD–PH(mutant K465E) PDK1 complex. Additionally, three KD–PH complexes were simulated, including a substrate analogue bound to a hydrophobic pocket (denominated the PIF-pocket) substrate-docking site. We find that only the PH-out conformation, with the PH domain well-oriented to interact with the cellular membrane, is active for wild-type PDK1. In contrast, the active conformation of the PDK1 K465E mutant is PH-in, being ATP-stable at the active site while the PIF-pocket is more accessible to the peptide substrate. We corroborate that both the docking-site binding and the catalytic activity are in fact enhanced in knock-in mouse samples expressing the PDK1 K465E protein, enabling the phosphorylation of PKB in the absence of PIP₃ binding.

LTA₄H is a bifunctional zinc metalloenzyme that converts leukotriene A₄ (LTA₄) into leukotriene B₄ (LTB₄), one of the most potent chemotactic agents involved in acute and chronic inflammatory diseases. In this reaction, LTA₄H acts as an epoxide hydrolase with a unique and fascinating mechanism, which includes the stereoselective attachment of one water molecule to the carbon backbone of LTA₄ several methylene units away from the epoxide moiety. By combining Molecular Dynamics simulations and Quantum Mechanics/Molecular Mechanics calculations, we obtained a very detailed molecular picture of the different consecutive steps of that mechanism. By means of a rather unusual 1,7-nucleophilic substitution through a clear S_N1 mechanism, the epoxide opens and the triene moiety of the substrate twists in such a way that the bond C₆-C₇ adopts its cis (Z) configuration, thus exposing the R face of C₁₂ to the addition of a water molecule hydrogen-bonded to ASP375. Thus, the two stereochemical features that are required for the bioactivity of LTB₄ appear to be closely related. The noncovalent π-π stacking interactions between the triene moiety and two tyrosines (TYR267 and, especially, TYR378) that wrap the triene system along the whole reaction explain the preference for the cis configuration inside LTA₄H.

Here, we describe the first systematic study on the mechanism of substrate-selective inhibition of mammalian ALOX15 orthologs. For this purpose, we prepared a series of N-substituted 5-(1H-indol-2-yl)anilines and found that (N-(5-(1H-indol-2-yl)-2-methoxyphenyl)sulfamoyl)carbamates and their monofluorinated analogues are potent and selective inhibitors of the linoleate oxygenase activity of rabbit and human ALOX15. Introduction of a 2-methoxyaniline moiety into the core pharmacophore plays a crucial role in substrate-selective inhibition of ALOX15-catalyzed oxygenation of linoleic acid at submicromolar concentrations without affecting arachidonic acid oxygenation. Steady-state kinetics, mutagenesis studies, and molecular dynamics (MD) simulations suggested an allosteric mechanism of action. Using a dimer model of ALOX15, our MD simulations suggest that the binding of the inhibitor at the active site of one monomer induces conformational alterations in the other monomer so that the formation of a productive enzyme–linoleic acid complex is energetically compromised.

Specialized pro-resolving lipid mediators (SPMs) are natural bioactive agents actively involved in inflammation resolution. SPMs act when uncontrolled inflammatory processes are developed, for instance, in patients of COVID-19 or other diseases. The so-called resolution pharmacology aims at developing new treatments based on the use of SPMs as agonists, which promote inflammation resolution without unwanted side effects. It has been shown that the biosynthesis of SPMs called eicosapentaenoic acid (EPA)-derived E-series resolvins is initiated by aspirin-acetylated COX-2 from EPA, leading to 18-hydroperoxy-eicosapentaenoic acid (18-HpEPE). However, there are many open questions concerning the intriguing role of aspirin in the molecular mechanism of resolvin formation. Our MD simulations, combined with QM/MM calculations, show that the potential energy barriers for the H₁₆-abstraction from EPA, required for forming 18-HpEPE, are higher than for the H₁₃-abstraction, thus explaining why 18-HpEPE is a marginal product of COX-2 catalysis. By contrast, in the aspirin-acetylated COX-2/EPA complex, the H_16proS-abstraction energy barriers are somewhat lower than the H_13proS energy barriers and much smaller than the H₁₆-transfer barriers in the wild type COX-2/EPA system. Those results agree with the experimental observation that aspirin favours the synthesis of several SPMs known as aspirin-triggered resolvins. In the following step of the catalytic mechanism, the calculated O₂ addition to C₁₈ is preferred versus the addition to C₁₄ which also agrees with 18R-HEPE and 18S-HEPE being the main products from EPA in aspirin-acetylated COX-2.