Carlos Navarro-Retamal

Flavonoids as CDK1 Inhibitors: Insights in Their Binding Orientations and Structure-Activity Relationship

In the last years, the interactions of flavonoids with protein kinases (PKs) have been described by using crystallographic experiments. Interestingly, different orientations have been found for one flavonoid inside different PKs and different chemical substitutions lead to different orientations of the flavonoid scaffold inside one PK. Accordingly, orientation predictions of novel analogues could help to the design of flavonoids with high PK inhibitory activities. With this in mind, we studied the binding modes of 37 flavonoids (flavones and chalcones) inside the cyclin-dependent PK CDK1 using docking experiments. We found that the compounds under study adopted two different orientations into the active site of CDK1 (orientations I and II in the manuscript). In addition, quantitative structure–activity relationship (QSAR) models using CoMFA and CoMSIA methodologies were constructed to explain the trend of the CDK1 inhibitory activities for the studied flavonoids. Template-based and docking-based alignments were used. Models developed starting from docking-based alignment were applied for describing the whole dataset and compounds with orientation I. Adequate R2 and Q2 values were obtained by each method; interestingly, only hydrophobic and hydrogen bond donor fields describe the differential potency of the flavonoids as CDK1 inhibitors for both defined alignments and subsets. Our current application of docking and QSAR together reveals important elements to be drawn for the design of novel flavonoids with increased PK inhibitory activities.

Structural and Affinity Determinants in the Interaction between Alcohol Acyltransferase from F. x ananassa and Several Alcohol Substrates: A Computational Study

Aroma and flavor are important factors of fruit quality and consumer preference. The specific pattern of aroma is generated during ripening by the accumulation of volatiles compounds, which are mainly esters. Alcohol acyltransferase (AAT) (EC 2.3.1.84) catalyzes the esterification reaction of aliphatic and aromatic alcohols and acyl-CoA into esters in fruits and flowers. In Fragaria x ananassa, there are different volatiles compounds that are obtained from different alcohol precursors, where octanol and hexanol are the most abundant during fruit ripening. At present, there is not structural evidence about the mechanism used by the AAT to synthesize esters. Experimental data attribute the kinetic role of this enzyme to 2 amino acidic residues in a highly conserved motif (HXXXD) that is located in the middle of the protein. With the aim to understand the molecular and energetic aspects of volatiles compound production from F. x ananassa, we first studied the binding modes of a series of alcohols, and also different acyl-CoA substrates, in a molecular model of alcohol acyltransferase from Fragaria x ananassa (SAAT) using molecular docking. Afterwards, the dynamical behavior of both substrates, docked within the SAAT binding site, was studied using routine molecular dynamics (MD) simulations. In addition, in order to correlate the experimental and theoretical data obtained in our laboratories, binding free energy calculations were performed; which previous results suggested that octanol, followed by hexanol, presented the best affinity for SAAT. Finally, and concerning the SAAT molecular reaction mechanism, it is suggested from molecular dynamics simulations that the reaction mechanism may proceed through the formation of a ternary complex, in where the Histidine residue at the HXXXD motif deprotonates the alcohol substrates. Then, a nucleophilic attack occurs from alcohol charged oxygen atom to the carbon atom at carbonyl group of the acyl CoA. This mechanism is in agreement with previous results, obtained in our group, in alcohol acyltransferase from Vasconcellea pubescens (VpAAT1).

Direct and Auger Electron-Induced, Single- and Double-Strand Breaks on Plasmid DNA Caused by 99mTc-Labeled Pyrene Derivatives and the Effect of Bonding Distance.

It is evident that 99mTc causes radical-mediated DNA damage due to Auger electrons, which were emitted simultaneously with the known γ-emission of 99mTc. We have synthesized a series of new 99mTc-labeled pyrene derivatives with varied distances between the pyrene moiety and the radionuclide. The pyrene motif is a common DNA intercalator and allowed us to test the influence of the radionuclide distance on damages of the DNA helix. In general, pUC 19 plasmid DNA enables the investigation of the unprotected interactions between the radiotracers and DNA that results in single-strand breaks (SSB) or double-strand breaks (DSB). The resulting DNA fragments were separated by gel electrophoresis and quantified by fluorescent staining. Direct DNA damage and radical-induced indirect DNA damage by radiolysis products of water were evaluated in the presence or absence of the radical scavenger DMSO. We demonstrated that Auger electrons directly induced both SSB and DSB in high efficiency when 99mTc was tightly bound to the plasmid DNA and this damage could not be completely prevented by DMSO, a free radical scavenger. For the first time, we were able to minimize this effect by increasing the carbon chain lengths between the pyrene moiety and the 99mTc nuclide. However, a critical distance between the 99mTc atom and the DNA helix could not be determined due to the significantly lowered DSB generation resulting from the interaction which is dependent on the type of the 99mTc binding motif. The effect of variable DNA damage caused by the different chain length between the pyrene residue and the Tc-core as well as the possible conformations of the applied Tc-complexes was supplemented with molecular dynamics (MD) calculations. The effectiveness of the DNA-binding 99mTc-labeled pyrene derivatives was demonstrated by comparison to non-DNA-binding 99mTcO4–, since nearly all DNA damage caused by 99mTcO4– was prevented by incubating with DMSO.

Energetic differences between non-domain-swapped and domain-swapped chain connectivities in the K2P potassium channel TRAAK

Two-pore domain (K2P) channels are twofold symmetric K+ channels which control cell excitability by enabling the leak of potassium ions from cells in response to physicochemical stimuli. Crystallization of K2P channels revealed the presence of several structural features, which include an external cap. In the available crystallographic structures, the cap is present as non-domain-swapped (NDS) and domain-swapped (DS) chain conformations, where DS chain conformation exchanges two opposing outer helices 180° around the channel. In this work, energy differences between the residues located at the highest point of the cap in NDS and DS conformations were evaluated for TRAAK, a K2P channel that was crystallized in both conformations. Results indicated a preference for DS conformation, but this result is not extensible to TASK K2P channels.

Mycobacterium tuberculosis serine/threonine protein kinases: structural information for the design of their specific ATP-competitive inhibitors

In the last decades, human protein kinases (PKs) have been relevant as targets in the development of novel therapies against many diseases, but the study of Mycobacterium tuberculosis PKs (MTPKs) involved in tuberculosis pathogenesis began much later and has not yet reached an advanced stage of development. To increase knowledge of these enzymes, in this work we studied the structural features of MTPKs, with focus on their ATP-binding sites and their interactions with inhibitors. PknA, PknB, and PknG are the most studied MTPKs, which were previously crystallized; ATP-competitive inhibitors have been designed against them in the last decade. In the current work, reported PknA, PknB, and PknG inhibitors were extracted from literature and their orientations inside the ATP-binding site were proposed by using docking method. With this information, interaction fingerprints were elaborated, which reveal the more relevant residues for establishing chemical interactions with inhibitors. The non-crystallized MTPKs PknD, PknF, PknH, PknJ, PknK, and PknL were also studied; their three-dimensional structural models were developed by using homology modeling. The main characteristics of MTPK ATP-binding sites (the non-crystallized and crystallized MTPKs, including PknE and PknI) were accounted; schemes of the main polar and nonpolar groups inside their ATP-binding sites were constructed, which are suitable for a major understanding of these proteins as antituberculotic targets. These schemes could be used for establishing comparisons between MTPKs and human PKs in order to increase selectivity of MTPK inhibitors. As a key tool for guiding medicinal chemists interested in the design of novel MTPK inhibitors, our work provides a map of the structural elements relevant for the design of more selective ATP-competitive MTPK inhibitors.

Insights into the Structural Requirements of 2(S)-Amino-6-Boronohexanoic Acid Derivatives as Arginase I Inhibitors: 3D-QSAR, Docking, and Interaction Fingerprint Studies

Human arginase I (hARGI) is an important enzyme involved in the urea cycle; its overexpression has been associated to cardiovascular and cerebrovascular diseases. In the last years, several congeneric sets of hARGI inhibitors have been reported with possible beneficial roles for the cardiovascular system. At the same time, crystallographic data have been reported including hARGI–inhibitor complexes, which can be considered for the design of novel inhibitors. In this work, the structure–activity relationship (SAR) of Cα substituted 2(S)-amino-6-boronohexanoic acid (ABH) derivatives as hARGI inhibitors was studied by using a three-dimensional quantitative structure–activity relationships (3D-QSAR) method. The predictivity of the obtained 3D-QSAR model was demonstrated by using internal and external validation experiments. The best model revealed that the differential hARGI inhibitory activities of the ABH derivatives can be described by using steric and electrostatic fields; the local effects of these fields in the activity are presented. In addition, binding modes of the above-mentioned compounds inside the hARGI binding site were obtained by using molecular docking. It was found that ABH derivatives adopted the same orientation reported for ABH within the hARGI active site, with the substituents at Cα exposed to the solvent with interactions with residues at the entrance of the binding site. The hARGI residues involved in chemical interactions with inhibitors were identified by using an interaction fingerprints (IFPs) analysis.

Folding and Lipid Composition Determine Membrane Interaction of the Disordered Protein COR15A

Plants from temperate climates, such as the model plant Arabidopsis thaliana, are challenged with seasonal low temperatures that lead to increased freezing tolerance in fall in a process termed cold acclimation. Among other adaptations, this involves the accumulation of cold-regulated (COR) proteins, such as the intrinsically disordered chloroplast-localized protein COR15A. Together with its close homolog COR15B, it stabilizes chloroplast membranes during freezing. COR15A folds into amphipathic α-helices in the presence of high concentrations of low-molecular-mass crowders or upon dehydration. Under these conditions, the (partially) folded protein binds peripherally to membranes. In our study, we have used coarse-grained molecular dynamics simulations to elucidate the details of COR15A-membrane binding and its effects on membrane structure and dynamics. Simulation results indicate that at least partial folding of COR15A and the presence of highly unsaturated galactolipids in the membranes are necessary for efficient membrane binding. The bound protein is stabilized on the membrane by interactions of charged and polar amino acids with galactolipid headgroups and by interactions of hydrophobic amino acids with the upper part of the fatty acyl chains. Experimentally, the presence of liposomes made from a mixture of lipids mimicking chloroplast membranes induces additional folding in COR15A under conditions of partial dehydration, in agreement with the simulation results.