Search

The potential of hydroxypyridinones for in vivo iron sequestration, in both biological and medical contexts, has been extensively discussed in the literature. Different chelators can be designed, with distinct lipophilicities that should alter their cell permeability, distribution, and rates of metabolism. However, for effective iron scavenging in biological systems, the redox potential and binding affinity of iron must fall within a proper range. Our objective was to assess the impact of different hydroxypyridinone chelators in 3:1 iron(III) complexes through comparison of these thermodynamic properties. For that purpose, we employed a cluster-continuum approach using density functional theory, on a dataset of 25 iron complexes. Whenever possible, our results were compared with experimental stability constants (log β) and with electrode potentials. We observed a good qualitative agreement between computed free energies of binding and log β values. In addition, we described which substitutions to the 3-hydroxypyridin-4-one ring should not markedly affect the redox properties and metal ion affinity considering iron. Graphical abstract Iron complexes of hydroxypyridinones.

Understanding the effect of glycation on the function of transferrin, the systemic iron transporter, is fundamental to fully grasp the mechanisms leading to the loss of iron homeostasis observed in diabetes mellitus (DM). The spontaneous reaction with protein amino groups is one of the main causes of glucose toxicity, but the site specificity of this reaction is still poorly understood. Here in, an in vitro approach was used to study human holo-transferrin glycation in detail. Lysine residues 103, 312 and 380 proved to be the most reactive sites, and overall glycation specificity was found to be remarkably different from that described for apo-transferrin. A computational biochemistry approach was subsequently applied to rationalize lysine reactivity. Even though pK a values, solvent accessible surface area, hydrogen bonds or the presence of nearby charged/polar residues could be related to lysine reactivity, these parameters do not suffice to describe glycation site specificity in holo-transferrin. Furthermore, analysis of the most reactive residues suggests that the correct lysine side chain orientation may play a fundamental role in reactivity. Nevertheless, in holo-transferrin, glycation occurs away from the iron-binding sites and, despite the observed iron release, the modification of apo-transferrin should play a more relevant role for the loss of iron-binding capacity observed in the blood serum of DM patients., (Copyright © 2018 Elsevier Inc. All rights reserved.)

Tyrosine ammonia lyase (TAL) catalyzes the deamination of tyrosine to p-coumaric acid in purple phototropic bacteria and Actinomycetales. The enzyme is used in bioengineering and has the potential to be used industrially. It belongs to a family of enzymes that uses a 4-methylidene-imidazole-5-one (MIO) cofactor to catalyze the deamination amino acids. In the present work, we used a QM/MM and a QM cluster models of TAL to explore two putative reaction paths for its catalytic mechanism. Part of the N-MIO mechanism was previously studied by computational methods. We improved on previous studies by using a larger, more complete model of the enzyme, and by describing the complete reaction path. The activation energy for this mechanism, in agreement with the previous study, is 28.5 kcal/mol. We also found another reaction path that has overall better kinetics and reaches the products in a single reaction step. The barrier for this Single-Step mechanism is 16.6 kcal/mol, which agrees very well with the experimental kcat of 16.0 kcal/mol. The geometrical parameters obtained for the cluster and QM/MM models are very similar, despite differences in the relative energies. This means that both approaches are capable of describing the correct catalytic path of TAL., (Copyright © 2015 Elsevier Inc. All rights reserved.)

The enzyme 3-hydroxy-3-methyl-glutaryl-CoA reductase (HMG-CoA-R) is the fundamental target for the treatment of hypercholesterolemia nowadays. The HMG-CoA-R clinical active site inhibitors (statins) are among the most widespread and profitable drugs ever sold but their side effects (myopathies, sometimes severe) still limit their use, which makes the finding of alternatives to statins a field of intense research. In this line, we address here a new strategy for inhibiting the homotetrameric HMG-CoA-R. The enzyme consists of a "dimer of dimers", each dimer having two active sites. We pursue here the inhibition of enzyme oligomerization, through drug binding to the dimer interface. We have computationally mutated 232 interfacial residues by alanine and calculated the loss in binding free energy among the monomers that build up each dimer of the homotetramer. This led to the identification of the (ten) key residues for the formation of the active dimer (Glu528, Ile531, Met534, Tyr644, Glu665, Asn686, Lys692, Lys735, Met742, and Val863). The results show that these residues are located in two specific spots of the protein with a cleft shape, whose shape and size is favorable for small drug binding. It is expectable that small molecules specifically bound to these druggable pockets will have a major effect on the oligomerization of the protein or/and in active site formation. This paves the way for the discovery of new families of inhibitors of HMG-CoA-R.

Considering our interest in the design of innovative radiometal-based complexes for in vivo imaging of nitric oxide synthase (NOS), we have recently introduced a set of M(CO)3-complexes (M=(99m)Tc, Re) containing a pendant N(ω)-NO2-L-arginine moiety, a known inhibitor of the enzyme. Enzymatic assays with purified inducible NOS have shown that the non-radioactive surrogates with 3-(Re1; Ki=84 μM) or 6-carbon linkers (Re2; Ki=6 μM) are stronger inhibitors than the respective metal-free conjugates L1 (Ki=178 μM) and L2 (Ki=36 μM), with Re2 displaying the highest inhibitory potency. Aiming to rationalize the experimental results we have performed a molecular docking study combined with molecular dynamics (MD) simulations and free energy perturbation (FEP) calculations. The higher inhibitory potency of Re2 arises from the stronger electrostatic interactions observed between the "Re(CO)3" core and the residues Arg260 and Arg382. This interaction is only possible due to the higher flexibility of its C6-carbon spacer, which links the N(ω)-NO2-L-arginine moiety and the "Re(CO)3" organometallic core. Furthermore, FEP calculations were carried out and the resultant relative binding energies (ΔΔGbind(calc)=0.690±0.028 kcal/mol,Re1/L1 and 1.825±0.318 kcal/mol, Re2/L2) are in accordance with the experimental results (ΔΔGbind(exp)=0.461±0.009 kcal/mol,Re1/L1 and 1.129±0.210 kcal/mol, Re2/L2); there is an energetic penalty for the transformation of the Re complexes into the ligands and this penalization is higher for the pair Re2/L2., (Copyright © 2013 Elsevier Inc. All rights reserved.)

In the present investigation, a computational analysis was performed on a data set comprised of human ether-a-go-go-related gene (hERG) blockers (triethanolamine, 1,3-thiazol-2-yl and tetrasubstituted imidazoline derivatives) in order to investigate the structural features required to reduce the hERG-induced cardiotoxicity problems in an early stage of drug discovery. The results derived from the quantitative structure-activity relationship (QSAR) analysis showed that the volume, surface area and shape descriptors (vsurf_) contributed significantly in all the models. This reveals that the hydrogen-bonding and hydrophilicity properties (vsurf_HB1, vsurf_CW4 and a_acc) on the van der Waals (vdW) surface of the molecule is negatively contributed for the hERG blocking activity and the hydrophobic property (vsurf_D6) and the total polar volume (vsurf_Wp2) on the vdW surface of the molecule are favourable for the activity. Further, the pharmacophore analysis also shows that the Aro/Hyd/Acc contour is one of the important biophore sites for the hERG blocking activity. This suggests that the presence of aromatic, hydrophobic and hydrogen-bonding groups in the molecules is favourable for interaction. In comparison with our earlier works (explaining the role of topological and hydrophobicity properties for the hERG blocking activity), these studies provided additional information on the importance of vdW surface area properties for the hERG blocking activity. These results can be used with other molecular modelling studies for the design of novel molecules that are free of cardiotoxicity.

Structural feature analysis of chlorogenic acid derivatives made up of varying lengths of alkyl groups as α-glucosidases inhibitors were performed by QSAR techniques. The statistically significant models derived from the study were validated by leave one out, Y-randomization and test set methods. The predictive capacity of the models was assessed by its validation parameters such as crossvalidated correlation coefficients (Q(2)), predictive residual analysis and other correlation parameters. The results obtained from the study show that the models were constructed with vsurf like properties (vsurf_ID4, vsurf_ID7 and vsurf_CW8), partial charge (Q_VSA_FNEG) and conformation dependent charged (dipoleX) descriptors. The integy moments of hydrophobicity descriptors (ID4 and ID7) are contributed for the inhibitory activity of the α-glucosidases enzymes of both the species. The vsurf_ID7 descriptor has contributed significantly (negatively) for the inhibitory activity prediction of α-glucosidases enzymes of S. cerevisiae. The partial negative charge on the surface of the molecules is detrimental for the activity, which reveals that the active site of the enzymes may have negatively charged groups. The pharmacophore analysis results also confirm the presence of hydrophilic properties on the vdW surface of the molecules. These results explain that the active sites of α-glucosidase enzymes of both the species have the same environment for the interaction. The alkyl side chain on the molecules is important for the pharmacokinetic behavior of the molecules and reduces the interaction energy of the molecules with the water. Hence, these results will be useful for designing novel molecules with multiple activities., (Copyright © 2012 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Dealing with receptor flexibility in docking methodology is still a problem. The main reason behind this difficulty is the large number of degrees of freedom that have to be considered in this kind of calculations. In this paper, we present an automated procedure, called MADAMM, that allows flexibilization of both the receptor and the ligand during a multistaged docking with an automated molecular modeling protocol. We show that the orientation of particular residues at the interface between the protein and the ligand have a crucial influence on the way they interact during the docking process, and the standard docking methodologies failed to predict their correct mode of binding. We present some examples that demonstrate the capabilities of this approach when compared with traditional docking methodologies., ((c) 2008 Wiley-Liss, Inc.)

Alanine scanning of protein-protein interfaces has shown that there are some residues in the protein-protein interfaces, responsible for most of the binding free energy, which are called hot spots. Hot spots tend to exist in densely packed central clusters, and a hypothesis has been proposed that considers that inaccessibility to the solvent must be a necessary condition to define a residue as a binding hot spot. This O-ring hypothesis is mainly based on the analysis of the accessible surface area (ASA) of 23 static, crystallographic structures of protein complexes. It is known, however, that protein flexibility allows for temporary exposures of buried interfacial groups, and even though the ASA provides a general trend of the propensity for hydration, protein/solvent-specific interactions or hydrogen bonding cannot be considered here. Therefore, a microscopic level, atomistic picture of hot spot solvation is needed to support the O-ring hypothesis. In this study, we began by applying a computational alanine-scanning mutagenesis technique, which reproduces the experimental results and allows for decomposing the binding free energy difference in its different energetic factors. Subsequently, we calculated the radial distribution function and residence times of the water molecules near the hot/warm spots to study the importance of the water environment around those energetically important amino acid residues. This study shows that within a flexible, dynamic protein framework, the warm/hot spot residues are, indeed, kept sheltered from the bulk solvent during the whole simulation, which allows a better interacting microenvironment.

In spite of the enormous interest that has been devoted to its study, the mechanism of the enzyme farnesyltransferase (FTase) remains the subject of several crucial doubts. In this article, we shed a new light in one of the most fundamental dilemmas that characterize the mechanism of this puzzling enzyme commonly referred to as the "distances paradox", which arises from the existence of a large 8-A distance between the two reactive atoms in the reaction catalyzed by this enzyme: a Zn-bound cysteine sulphur atom from a peptidic substrate and the farnesyldiphosphate (FPP) carbon 1. This distance must be overcome for the reaction to occur. In this study, the two possible alternatives were evaluated by combining molecular mechanics (AMBER) and quantum chemical calculations (B3LYP). Basically, our results have shown that an activation of the Zn-bound cysteine thiolate with subsequent displacement from the zinc coordination sphere towards the FPP carbon 1 is not a realistic hypothesis of overcoming the large distance reported in the crystallographic structures of the ternary complexes between the two reactive atoms, but that a rotation involving the FPP molecule can bring the two atoms closer with moderate energetic cost, coherent with previous experimental data. This conclusion opens the door to an understanding of the chemical step in the farnesylation reaction., ((c) 2006 Wiley-Liss, Inc.)

HIV-1 RT is one of the most important antiviral targets in the treatment of acquired immunodeficiency syndrome (AIDS). Several crystallographic structures are available for this enzyme, mostly with bound inhibitors. Despite their importance for structure based drug design towards new anti-HIV retrovirals, the X-ray structures of the unliganded enzyme could only be obtained incomplete, with a low resolution and until recently even the conformation of the p66 thumb was controversial. In this work we have aligned different X-ray RT structures, and built up a computational model of RT using homology modeling, which was afterwards refined and validated through MD simulations with explicit solvent. The model enzyme was structurally stable through the whole MD simulation, showing a RMSD of 2 Angstrom from the starting geometry. The Ramanchandram plot has improved along the simulation. Both intra-domain and interdomain movements were observed. The thumb kept its closed conformation through the whole simulation. A contact map, hydration sites study and a detailed analysis of the solvation of the nucleotide binding site are also presented.

This article focuses on the second step of the catalytic mechanism for the reduction of ribonucleotides catalyzed by the enzyme Ribonucleotide Reductase (RNR). This step corresponds to the protonation/elimination of the substrate's C-2' hydroxyl group. Protonation is accomplished by the neighbor Cys-225, leading to the formation of one water molecule. This is a very relevant step since most of the known inhibitors of this enzyme, which are already used in the fight against certain forms of cancer, are 2'-substituted substrate analogs. Even though some theoretical studies have been performed in the past, they have modeled the enzyme with minimal gas-phase models, basically represented by a part of the side chain of the relevant amino acids, disconnected from the protein backbone. This procedure resulted in a limited accuracy in the position and/or orientation of the participating residues, which can result in erroneous energetics and even mistakes in the choice of the correct mechanism for this step. To overcome these limitations we have used a very large model, including a whole R1 model with 733 residues plus the substrate and 10 A thick shell of water molecules, instead of the minimal gas-phase models used in previous works. The ONIOM method was employed to deal with such a large system. This model can efficiently account for the restrained mobility of the reactive residues, as well as the long-range enzyme-substrate interactions. The results gave additional information about this step, which previous small models could not provide, allowing a much clearer evaluation of the role of the enzyme. The interaction energy between the enzyme and the substrate along the reaction coordinate and the substrate steric strain energy have been obtained. The conclusion was that the barrier obtained with the present model was very similar to the one previously determined with minimal gas-phase models. Therefore, the role of the enzyme in this step was concluded to be mainly entropic, rather than energetic, by placing the substrate and the two reactive residues in a position that allows for the highly favorable concerted trimolecular reaction, and to protect the enzyme radical from the solvent.