The dipeptidyl peptidase-4 inhibitory effect of flavonoids is hindered in protein rich environments

Dipeptidyl peptidase-4 (DPP-4) inhibitors present a unique approach for the management of type 2 diabetes (T2D). In the present study, the inhibition of DPP-4 was evaluated for a large panel of flavonoids, important components of the human diet, using in vitro and ex vivo models. The activity of the isolated enzyme was assayed in vitro. Subsequently, the most active flavonoids were tested ex vivo in human whole blood and plasma. In this study, contrary to the in vitro fluorometric tests, flavonoids did not show inhibitory activity against DPP-4. Due to the discrepancy in the results between the in vitro and ex vivo approaches, plasma protein binding values were determined, presenting values from 43.9 to 100.0%. This work provides a new insight into the inhibitory activity for DPP-4, based on the flavonoid scaffold. Additionally, the obtained results showed that the inhibitory effect of flavonoids against DPP-4 was hindered in protein rich environments, like that occurring in blood, and indicated the need for experimental refinement in drug discovery for blood targets.

Prevalence and incidence of surgical site infections in the European Union/European Economic Area: how do these measures relate?

BACKGROUND

In 2011-2012, the European Centre for Disease Prevention and Control (ECDC) initiated the first European point prevalence survey (PPS) of healthcare-associated infections (HCAIs) in addition to targeted surveillance of the incidence of specific types of HCAI such as surgical site infections (SSIs).

AIM

To investigate whether national and multi-country SSI incidence can be estimated from ECDC PPS data.

METHODS

In all, 159 hospitals were included from 15 countries that participated in both ECDC surveillance modules, aligning surgical procedures in the incidence surveillance to corresponding specialties from the PPS. National daily prevalence of SSIs was simulated from the incidence surveillance data, the Rhame and Sudderth (R&S) formula was used to estimate national and multi-country SSI incidence from the PPS data, and national incidence per specialty was predicted using a linear model including data from the PPS.

FINDINGS

The simulation of daily SSI prevalence from incidence surveillance of SSIs showed that prevalence fluctuated randomly depending on the day of measurement. The correlation between the national aggregated incidence estimated with R&S formula and observed SSI incidence was low (correlation coefficient = 0.24), but specialty-specific incidence results were more reliable, especially when the number of included patients was large (correlation coefficients ranging from 0.40 to 1.00). The linear prediction model including PPS data had low proportion of explained variance (0.40).

CONCLUSION

Due to a lack of accuracy, use of PPS data to estimate SSI incidence is recommended only in situations where incidence surveillance of SSIs is not performed, and where sufficiently large samples of PPS data are available.

A computational study on the redox properties and binding affinities of iron complexes of hydroxypyridinones

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.

Evaluation of a flavonoids library for inhibition of pancreatic α-amylase towards a structure-activity relationship

α-Amylase has been considered an important therapeutic target for the management of type 2 diabetes mellitus (T2DM), decreasing postprandial hyperglycaemia (PPHG). In the present work, a panel of 40 structurally related flavonoids was tested, concerning their ability to inhibit α-amylase activity, using a microanalysis screening system, an inhibitory kinetic analysis and molecular docking calculations. From the obtained results, it was possible to observe that the flavone with a -Cl ion at 3-position of C-ring, an –OH group at 3′- and 4′- positions of B-ring and at 5- and 7- positions of A-ring and the C2 = C3 double bond, was the most active tested flavonoid, through competitive inhibition. In conclusion, some of the tested flavonoids have shown promising inhibition of α-amylase and may be considered as possible alternatives to the modulation of T2DM.

Structural and mechanistic aspects of S-S bonds in the thioredoxin-like family of proteins

Disulfide bonds play a critical role in a variety of structural and mechanistic processes associated with proteins inside the cells and in the extracellular environment. The thioredoxin family of proteins like thioredoxin (Trx), glutaredoxin (Grx) and protein disulfide isomerase, are involved in the formation, transfer or isomerization of disulfide bonds through a characteristic thiol-disulfide exchange reaction. Here, we review the structural and mechanistic determinants behind the thiol-disulfide exchange reactions for the different enzyme types within this family, rationalizing the known experimental data in light of the results from computational studies. The analysis sheds new atomic-level insight into the structural and mechanistic variations that characterize the different enzymes in the family, helping to explain the associated functional diversity. Furthermore, we review here a pattern of stabilization/destabilization of the conserved active-site cysteine residues presented beforehand, which is fully consistent with the observed roles played by the thioredoxin family of enzymes.

β-Adrenoceptors Trigger Melatonin Synthesis in Phagocytes

Melatonin (5-methoxy-N-acetylserotonin), the pineal hormone, is also synthesized by immune-competent cells. The pineal hormone signals darkness, while melatonin synthesized on demand by activated macrophages at any hour of the day acts locally, favoring regulatory/tolerant phenotypes. Activation of β-adrenoceptors in pinealocytes is the main route for triggering melatonin synthesis. However, despite the well-known role of β-adrenoceptors in the resolution macrophage phenotype (M2), and the relevance of macrophage synthesized melatonin in facilitating phagocytic activity, there is no information regarding whether activation of β-adrenoceptors would induce melatonin synthesis by monocytes. Here we show that catecholamines stimulate melatonin synthesis in bone marrow-derived dendritic cells and RAW 264.7 macrophages. Activation of β-adrenoceptors promotes the synthesis of melatonin by stimulating cyclic AMP/protein kinase A (PKA) pathway and by activating the nuclear translocation of NF-κB. Considering the great number of macrophages around sympathetic nerve terminals, and the relevance of this system for maintaining macrophages in stages compatible to low-grade inflammation, our data open the possibility that extra-pineal melatonin acts as an autocrine/paracrine signal in macrophages under resolution or tolerant phenotypes.

Determining the glycation site specificity of human holo-transferrin

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.

A QM/MM study of the reaction mechanism of human β-ketoacyl reductase

Human fatty acid synthase (hFAS) is a multifunctional enzyme involved in a wide diversity of biological functions. For instance, it is a precursor of phospholipids and other complex processes such as the de novo synthesis of long chain fatty acid. Human FAS is also a component of biological membranes and it is implicated in the overexpression of several types of cancers. In this work, we describe the catalytic mechanism of β-ketoreductase (KR), which is a catalytic domain of the hFAS enzyme that catalyzes the reduction of β-ketoacyl to β-hydroxyacyl with the concomitant oxidation of the NADPH cofactor. The catalysis by KR is an intermediate step in the cycle of reactions that elongate the substrate's carbon chain until the final product is obtained. We study and propose the catalytic mechanism of the KR domain determined using the hybrid QM/MM methodology, at the ONIOM(B3LYP/6-311+G(2d,2p):AMBER) level of theory. The results indicate that the reaction mechanism occurs in two stages: (i) nucleophilic attack by a NADPH hydride to the β-carbon of the substrate, together with an asynchronous deprotonation of the Tyr2034 by the oxygen of the β-alkoxide to hold the final alcohol product; and (ii) an asynchronous deprotonation of the hydroxyl in the NADP⁺'s ribose by Tyr2034, and of the Lys1995 by the resulting alkoxide in the former ribose to restore the protonation state of Tyr2034. The reduction step occurs with a Gibbs energy barrier of 11.7 kcal mol^-1 and a Gibbs reaction energy of -10.6 kcal mol^-1. These results have provided an understanding of the catalytic mechanism of the KR hFAS domain, a piece of the heavy hFAS biosynthetic machinery.

Mechanistic Insights on Human Phosphoglucomutase Revealed by Transition Path Sampling and Molecular Dynamics Calculations

Human α-phosphoglucomutase 1 (α-PGM) catalyzes the isomerization of glucose-1-phosphate into glucose-6-phosphate (G6P) through two sequential phosphoryl transfer steps with a glucose-1,6-bisphosphate (G16P) intermediate. Given that the release of G6P in the gluconeogenesis raises the glucose output levels, α-PGM represents a tempting pharmacological target for type 2 diabetes. Here, we provide the first theoretical study of the catalytic mechanism of human α-PGM. We performed transition-path sampling simulations to unveil the atomic details of the two catalytic chemical steps, which could be key for developing transition state (TS) analogue molecules with inhibitory properties. Our calculations revealed that both steps proceed through a concerted S_N 2-like mechanism, with a loose metaphosphate-like TS. Even though experimental data suggests that the two steps are identical, we observed noticeable differences: 1) the transition state ensemble has a well-defined TS region and a late TS for the second step, and 2) larger coordinated protein motions are required to reach the TS of the second step. We have identified key residues (Arg23, Ser117, His118, Lys389), and the Mg²⁺ ion that contribute in different ways to the reaction coordinate. Accelerated molecular dynamics simulations suggest that the G16P intermediate may reorient without leaving the enzymatic binding pocket, through significant conformational rearrangements of the G16P and of specific loop regions of the human α-PGM.

Properties that rank protein:protein docking poses with high accuracy

The development of docking algorithms to predict near-native structures of protein:protein complexes from the structure of the isolated monomers is of paramount importance for molecular biology and drug discovery.

A QM/MM approach on the structural and stereoelectronic factors governing glycosylation by GTF-SI fromStreptococcus mutans

This manuscript contains novel insights into the reaction mechanism catalyzed by GTF-SI. Structural and electronic features of the system are revealed, such as the strong hydrogen bond depicted above.

Membrane partition of bis-(3-hydroxy-4-pyridinonato) zinc(ii) complexes revealed by molecular dynamics simulations

The class of 3-hydroxy-4-pyridinone ligands is widely known and valuable for biomedical and pharmaceutical purposes. Their chelating properties towards biologically-relevant transition metal ions highlight their potential biomedical utility. A set of 3-hydroxy-4-pyridinone Zn(ii) complexes at different concentrations was studied for their ability to interact with lipid phases. We employed umbrella sampling simulations to attain the potential-of-mean force for a set of ligands and one Zn(ii) complex, as these permeated a 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) hydrated bilayer system. In addition, we used conventional molecular dynamics simulations to study the behavior of various Zn(ii) complexes in hydrated bilayer systems. This work discusses: (i) the partition of 3-hydroxy-4-pyridinone ligands to bilayer phases; (ii) self-aggregation in crowded environments of Zn(ii) complexes; and (iii) possible mechanisms for the membrane translocation of Zn(ii) complexes. We observed distinct interactions for the studied complexes, and distinct membrane partition coefficients (K_mem) depending on the considered ligand. The more hydrophobic ligand, 1-hexyl-3-hydroxy-2-methyl-4(1H)-pyridinone, partitioned more favorably to lipid phases (at least two orders of magnitude higher K_mem when compared to the other ligands), and the corresponding Zn(ii) complex was also prone to self-aggregation when an increased concentration of the complex was employed. We also observed that the inclusion of a coordinated water molecule in the parameterization of the Zn(ii) coordination sphere, as proposed in the available crystallographic structure of the complex, decreased the partition coefficient and membrane permeability for the tested complex.

The membrane partition of hydroxypyridinones and of zinc complexes explored by molecular dynamics.

Immune-pineal axis - acute inflammatory responses coordinate melatonin synthesis by pinealocytes and phagocytes

Melatonin is well known for its circadian production by the pineal gland, and there is a growing body of data showing that it is also produced by many other cells and organs, including immune cells. The chronobiotic role of pineal melatonin, as well as its protective effects in vitro and in vivo, have been extensively explored. However, the interaction between the chronobiotic and defence functions of endogenous melatonin has been little investigated. This review details the current knowledge regarding the coordinated shift in melatonin synthesis from the pineal gland (circadian and monitoring roles) to the regulation of acute immune responses via immune cell production and autocrine effects, producing systemic interactions termed the immune-pineal axis. An acute inflammatory response drives the transcription factor, NFκB, to switch melatonin synthesis from pinealocytes to macrophages/microglia and, upon acute inflammatory resolution, back to pinealocytes. The potential pathophysiological relevance of immune-pineal axis dysregulation is highlighted, with both research and clinical implications, across several medical conditions, including host/parasite interaction, neurodegenerative diseases and cancer. LINKED ARTICLES: This article is part of a themed section on Recent Developments in Research of Melatonin and its Potential Therapeutic Applications. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v175.16/issuetoc.

Improving the Catalytic Power of the DszD Enzyme for the Biodesulfurization of Crude Oil and Derivatives

The enhancement of the catalytic power of enzymes is a subject of enormous interest both for science and for industry. The latter, in particular, due to the vast applications enzymes can have in industrial processes, for instance in the desulfurization of crude oil, which is mandatory by law in many developed countries and is currently performed using costly chemical processes. In this work we sought to enhance the turnover rate of DszD from Rhodococcus erythropolis, a NADH-FMN oxidoreductase responsible for supplying FMNH₂ to DszA and DszC in the biodesulfurization process of crude oil, the 4S pathway. For this purpose, we replaced the wild type spectator residue of the rate limiting step of the reduction of FMN to FMNH₂ , a process catalysed by DszD and known to play an important role in the reaction energy profile. As replacements, we used all the naturally occurring amino acids, one at a time, using computational methodologies, and repeated the above-mentioned reaction with each mutant. To calculate the different free energy profiles, one for each mutated model, we applied quantum mechanics/molecular mechanics (QM/MM) methods within an ONIOM scheme. The free energy barriers obtained varied between 15.1 and 29.9 kcal mol^-1 . Multiple factors contributed to the different ΔG values. The most relevant were electrostatic interactions and the induction of a favourable alignment between substrate and cofactor. These results confirm the great potential that chirurgic mutations have for increasing the catalytic power of DszD in relation to the wild type (wt) enzyme.

Inhibition of protein tyrosine phosphatase 1B by flavonoids: A structure - activity relationship study

The classical non-transmembrane protein tyrosine phosphatase 1B (PTP1B) has emerged as a key negative regulator of insulin signaling pathways that leads to insulin resistance, turning this enzyme a promising therapeutic target in the management of type 2 diabetes mellitus (T2DM). In the present work, the in vitro inhibitory activity of a panel of structurally related flavonoids, for recombinant human PTP1B was studied and the type of inhibition of the most active compounds further evaluated. The majority of the studied flavonoids was tested in this work for the first time, including flavonoid C13, which was the most potent inhibitor. It was observed that the ability to inhibit PTP1B depends on the nature, position and number of substituents in the flavonoid structure, as the presence of both 7- and 8-OBn groups in the A ring, together with the presence of both 3' and 4'-OMe groups in the B ring and the 3-OH group in the C ring; these substituents increase the flavonoids' ability to inhibit PTP1B. In conclusion, some of the tested flavonoids seem to be promising PTP1B inhibitors and potential effective agents in the management of T2DM, by increasing insulin sensitivity.

α-Glucosidase inhibition by flavonoids: an in vitro and in silico structure-activity relationship study

α-Glucosidase inhibitors are described as the most effective in reducing post-prandial hyperglycaemia (PPHG) from all available anti-diabetic drugs used in the management of type 2 diabetes mellitus. As flavonoids are promising modulators of this enzyme’s activity, a panel of 44 flavonoids, organised in five groups, was screened for their inhibitory activity of α-glucosidase, based on in vitro structure–activity relationship studies. Inhibitory kinetic analysis and molecular docking calculations were also applied for selected compounds. A flavonoid with two catechol groups in A- and B-rings, together with a 3-OH group at C-ring, was the most active, presenting an IC₅₀ much lower than the one found for the most widely prescribed α-glucosidase inhibitor, acarbose. The present work suggests that several of the studied flavonoids have the potential to be used as alternatives for the regulation of PPHG.

The mechanism of the Ser-(cis)Ser-Lys catalytic triad of peptide amidases

In this paper, we report a theoretical investigation of the catalytic mechanism of peptide amidases that involve a Ser-(cis)Ser-Lys catalytic triad. Previous suggestions propose that these enzymes should follow a distinct catalytic mechanism from the one that is present in the classic Ser-His-Asp catalytic triad. The theoretical and computational results obtained in this work indicate the opposite idea, showing that both mechanisms are very similar and only few differences are observed between both reactions. The results reveal that the different alignment of the Ser-(cis)Ser-Lys catalytic triad in relation to the classical Ser-His-Asp triad may provide a better stabilisation of the reaction intermediates, and therefore make these enzymes catalytically more efficient. The catalytic mechanism has been determined at the M06-2X/6-311++G**//B3LYP/6-31G* level of theory and requires five sequential steps instead of the two that are generally proposed: (i) nucleophilic attack of serine on the carbonyl group of the substrate, forming the first tetrahedral intermediate, (ii) formation of an acyl-enzyme complex, (ii) release of an ammonia product, (iv) nucleophilic attack of a water molecule forming the second tetrahedral intermediate, and (iv) the release of the product of the reaction, the carboxylic acid. The computational results suggest that the rate-limiting step is the first one that requires an activation free energy of 15.93 kcal mol^-1. This result agrees very well with the available experimental data that predict a reaction rate of 2200 s^-1, which corresponds to a free energy barrier of 14 kcal mol^-1.

Revisiting Partition in Hydrated Bilayer Systems

We discuss potential-of-mean force convergence issues involving the methods for generating starting structures of umbrella sampling simulations for the study of the water-bilayer permeation of drug-like compounds. We provide a proof of concept of the benefits of potential-energy modification methods, in particular the flooding technique, to enhance the sampling and improve the convergence of the potential-of-mean force. In addition, we also discuss how these potential modifiers could introduce the necessary bias to induce the correct description of other relevant internal degrees of freedom for the compounds, thus providing a more accurate description of the permeation process. Effective methodological approximations are taken into consideration, more specifically the inclusion of atomic polarization, and experimental data are used in the validation of computed values whenever possible. The parametrization of the atomic point charges is always a relevant issue in biomolecular simulations. We discuss the influence of a molecule's conformation-dependent atomic point charges and their impact on the predicted potential-of-mean force.

Insights into the reaction mechanism of 3-O-sulfotransferase through QM/MM calculations

3-O-Sulfotransferase (3-OST) is one of the enzymes involved in heparan sulfate (HS) biosynthesis. HSs are polysaccharides with variable patterns of sulfation and acetylation that serve as entry receptors for herpes simplex virus type 1 (HSV-1). 3-OST is responsible for the transfer of a sulfate group from 3'-phosphoadenosine-5'-phosphosulfate (PAPS) to glucosamine units of HS. In this work, the catalytic mechanism of 3-OST was studied with atomic detail, using computational methods. We investigated the protonation state of key residues using the H++ web-based pKa prediction tool and molecular dynamics (MD) simulations and estimated the most relevant protonation state of the catalytic residues during catalysis. Catalytic histidine (His186) is predominantly protonated, while catalytic aspartate and glutamate (Asp189 and Glu184) are predominantly deprotonated. Subsequently, to study the catalytic mechanism, we applied a QM/MM method at the ONIOM(B3LYP/6-31G(d):ff94) level, starting from three geometries extracted from the 3, 6 and 8 ns point on the MD simulation. The results show that the reaction mechanism of 3-OST occurs by a single elementary step, consisting of an associative SN2 transfer of the sulfate group from PAPS to the HS glucosamine units, with the transfer of a proton from glucosamine to the catalytic Glu184. The activation free energies for this reaction were determined at the ONIOM(M06-2X-D3/6-311++G(2d,2p):ff94//B3LYP/6-31G(d):ff94) level of theory. Despite the free energy differences among the three conformations (10.2, 20.9 and 16.1 kcal mol(-1)), our results are consistent with the upper limit determined experimentally for the full cycle (20.4 kcal mol(-1)). The data obtained in this study will be useful for further studies on the inhibition of this enzyme, which is a useful target for drugs that block HSV-1 viral infections.

Unique Triphenylphosphonium Derivatives for Enhanced Mitochondrial Uptake and Photodynamic Therapy

In this study, unique methyl-functionalized derivatives (T*PP⁺) of the drug carrier triphenylphosphonium (TPP⁺) that exhibit significant enhancement of the accumulation of both the cation and its conjugated cargo in cell mitochondria are designed. We show that the presence of methyl group(s) at key positions within the phenyl ring results in an increase in the hydrophobicity and solvent accessible surface area of T*PP⁺. In particular, when the para position of the phenyl ring in T*PP⁺ is functionalized with a methyl group, the cation is most exposed to the surrounding environment, leading to a large decrease in water entropy and an increase in the level of van der Waals interaction with and partition into a nonpolar solvent. Therefore, stronger binding between the hydrophobic T*PP⁺ and mitochondrial membrane occurs. This is exemplified in a (hexachloro-fluorescein)-TPP⁺ conjugate system, where an ∼12 times increase in the rate of mitochondrial uptake and a 2 times increase in photodynamic therapy (PDT) efficacy against HeLa and FU97 cancer cells are achieved when TPP⁺ is replaced with T*PP⁺. Importantly, nearly all the FU97 cells treated with the (hexachloro-fluorescein)-T*PP⁺ conjugate are killed as compared to only half the population of cells in the case of the (hexachloro-fluorescein)-TPP⁺ conjugate at a similar PDT light dosage. This study thus forms a platform for the healthcare community to explore alternative TPP⁺ derivatives that can act as optimal drug transporters for enhanced mitochondrially targeted therapies.

Dual Effect of Catecholamines and Corticosterone Crosstalk on Pineal Gland Melatonin Synthesis

BACKGROUND/AIM

The nocturnal production of melatonin by the pineal gland is triggered by sympathetic activation of adrenoceptors and may be modulated by immunological signals. The effect of glucocorticoids on nocturnal melatonin synthesis is controversial; both stimulatory and inhibitory effects have been reported. During pathophysiological processes, an increased sympathetic tonus could result in different patterns of adrenoceptor activation in the pineal gland. Therefore, in this investigation, we evaluated whether the pattern of adrenergic stimulation of the pineal gland drives the direction of the glucocorticoid effect on melatonin production.

METHODS

The corticosterone effect on the pineal hormonal production induced by β-adrenoceptor or β+α1-adrenoceptor activation was evaluated in cultured glands. We also investigated whether the in vivo lipopolysaccharide (LPS)-induced inhibition of melatonin is dependent on the interaction of glucocorticoids and the α1-adrenoceptor in adrenalectomized animals and on the in vivo blockade of glucocorticoid receptors (GRs) or the α1-adrenoceptor.

RESULTS

Corticosterone potentiated β-adrenoceptor-induced pineal melatonin synthesis, whilst corticosterone-dependent inhibition was observed when melatonin production was induced by β+α1-adrenoceptors agonists. The inhibitory effect of corticosterone is mediated by GR, as it was abolished in the presence of a GR antagonist. Moreover, LPS-induced reduction in melatonin nocturnal plasma content was reversed by adrenalectomy and by antagonizing GR or α1-adrenoceptors.

CONCLUSIONS

The dual effect of corticosterone on pineal melatonin synthesis is determined by the activation pattern of adrenoceptors (β or β+α1) in the gland during GR activation, suggesting that increased activation of the sympathetic system and the hypothalamic-pituitary-adrenal axis are necessary for the control of melatonin production during defense responses.

Cholesterol Biosynthesis: A Mechanistic Overview

Cholesterol is an essential component of cell membranes and the precursor for the synthesis of steroid hormones and bile acids. The synthesis of this molecule occurs partially in a membranous world (especially the last steps), where the enzymes, substrates, and products involved tend to be extremely hydrophobic. The importance of cholesterol has increased in the past half-century because of its association with cardiovascular diseases, which are considered one of the leading causes of death worldwide. In light of the current need for new drugs capable of controlling the levels of cholesterol in the bloodstream, it is important to understand how cholesterol is synthesized in the organism and identify the main enzymes involved in this process. Taking this into account, this review presents a detailed description of several enzymes involved in the biosynthesis of cholesterol. In this regard, the structure and catalytic mechanism of the enzymes involved in cholesterol biosynthesis, from the initial two-carbon acetyl-CoA building block, will be reviewed and their current pharmacological importance discussed. We believe that this review may contribute to a deeper level of understanding of cholesterol metabolism and that it will serve as a useful resource for future studies of the cholesterol biosynthesis pathway.

Glutamine Synthetase Drugability beyond Its Active Site: Exploring Oligomerization Interfaces and Pockets

BACKGROUND

Glutamine synthetase (GS) is a crucial enzyme to the nitrogen cycle with great commercial and pharmaceutical value. Current inhibitors target the active site, affecting GS activity indiscriminately in all organisms. As the active site is located at the interface between two monomers, the protein-protein interface (PPI) of GSs gains a new role, by providing new targets for enzyme inhibition. Exploring GSs PPI could allow for the development of inhibitors selective for specific organisms. Here we map the PPI of three GSs-human (hsGS), maize (zmGS) and Mycobacterium tuberculosis (mtGS)-and unravel new drugable pockets.

METHODS

The PPI binding free energy coming from key residues on three GSs from different organisms were mapped by computational alanine scan mutagenesis, applying a multiple dielectric constant MM-PBSA methodology. The most relevant residues for binding are referred as hot-spots. Drugable pockets on GS were detected with the Fpocket software.

RESULTS AND CONCLUSIONS

A total of 23, 19 and 30 hot-spots were identified on hsGS, zmGS and mtGS PPI. Even possessing differences in the hot-spots, hsGS and zmGS PPI are overall very similar. On the other hand, mtGS PPI differs greatly from hsGS and zmGS PPI. A novel drugable pocket was detected on the mtGS PPI. It seems particularly promising for the development of selective anti-tuberculosis drugs given its location on a PPI region that is highly populated with hot-spots and is completely different from the hsGS and zmGS PPIs. Drugs targeting this pockets should be inactive on eukaryotic GS II enzymes.

The Catalytic Mechanism of the Marine-Derived Macrocyclase PatGmac

Cyclic peptides are a class of compounds with high therapeutic potential, possessing bioactivities including antitumor and antiviral (including anti-HIV). Despite their desirability, efficient design and production of these compounds has not been achieved to date. The catalytic mechanism of patellamide macrocyclization by the PatG macrocyclase domain has been computationally investigated by using quantum mechanics/molecular mechanics methodology, specifically ONIOM(M06/6-311++G(2d,2p):ff94//B3LYP/6-31G(d):ff94). The mechanism proposed herein begins with a proton transfer from Ser783 to His 618 and from the latter to Asp548. Nucleophilic attack of Ser783 on the substrate leads to the formation of an acyl-enzyme covalent complex. The leaving group Ala-Tyr-Asp-Gly (AYDG) of the substrate is protonated by the substrate's N terminus, leading to the breakage of the P1-P1' bond. Finally, the substrate's N terminus attacks the P1 residue, decomposing the acyl-enzyme complex forming the macrocycle. The formation and decomposition of the acyl-enzyme complex have the highest activation free energies (21.1 kcal mol(-1) and 19.8 kcal mol(-1) respectively), typical of serine proteases. Understanding the mechanism behind the macrocyclization of patellamides will be important to the application of the enzymes in the pharmaceutical and biotechnological industries.

Establishing the catalytic mechanism of human pancreatic α-amylase with QM/MM methods

In this work, we studied the catalytic mechanism of human pancreatic α-amylase (HPA). Our goal was to determine the catalytic mechanism of HPA with atomic detail using computational methods. We demonstrated that the HPA catalytic mechanism consists of two steps, the first of which (glycosylation step) involves breaking the glycosidic bond to culminate in the formation of a covalent intermediate. The second (deglycosylation step) consists of the addition of a water molecule to release the enzyme/substrate covalent intermediate, completing the hydrolysis of the sugar. The active site was very open to the solvent. Our mechanism basically differs from the previously proposed mechanism by having two water molecules instead of only one near the active site that participate in the mechanism. We also demonstrate the relevant role of the three catalytic amino acids, two aspartate residues and a glutamate (D197, E233, and D300), during catalysis. It was also shown that the rate limiting step was glycosylation, and its activation energy was in agreement with experimental values obtained for HPA. The experimental activation energy was 14.4 kcal mol(-1), and the activation energy obtained computationally was 15.1 kcal mol(-1).

Computational alanine scanning mutagenesis--an improved methodological approach

Alanine scanning mutagenesis of protein-protein interfacial residues can be applied to a wide variety of protein complexes to understand the structural and energetic characteristics of the hot-spots. Binding free energies have been estimated with reasonable accuracy with empirical methods, such as Molecular Mechanics/Poisson-Boltzmann surface area (MM-PBSA), and with more rigorous computational approaches like Free Energy Perturbation (FEP) and Thermodynamic Integration (TI). The main objective of this work is the development of an improved methodological approach, with less computational cost, that predicts accurately differences in binding free energies between the wild-type and alanine mutated complexes (DeltaDeltaG(binding)). The method was applied to three complexes, and a mean unsigned error of 0.80 kcal/mol was obtained in a set of 46 mutations. The computational method presented here achieved an overall success rate of 80% and an 82% success rate in residues for which alanine mutation causes an increase in the binding free energy > 2.0 kcal/mol (warm- and hot-spots). This fully atomistic computational methodological approach consists in a computational Molecular Dynamics simulation protocol performed in a continuum medium using the Generalized Born model. A set of three different internal dielectric constants, to mimic the different degree of relaxation of the interface when different types of amino acids are mutated for alanine, have to be used for the proteins, depending on the type of amino acid that is mutated. This method permits a systematic scanning mutagenesis of protein-protein interfaces and it is capable of anticipating the experimental results of mutagenesis, thus guiding new experimental investigations.

Including sheath effects in the interpretation of planar retarding potential analyzer's low-energy ion data

The interpretation of planar retarding potential analyzers (RPA) during ionospheric sounding rocket missions requires modeling the thick 3D plasma sheath. This paper overviews the theory of RPAs with an emphasis placed on the impact of the sheath on current-voltage (I-V) curves. It then describes the Petite Ion Probe (PIP) which has been designed to function in this difficult regime. The data analysis procedure for this instrument is discussed in detail. Data analysis begins by modeling the sheath with the Spacecraft Plasma Interaction System (SPIS), a particle-in-cell code. Test particles are traced through the sheath and detector to determine the detector's response. A training set is constructed from these simulated curves for a support vector regression analysis which relates the properties of the I-V curve to the properties of the plasma. The first in situ use of the PIPs occurred during the MICA sounding rocket mission which launched from Poker Flat, Alaska in February of 2012. These data are presented as a case study, providing valuable cross-instrument comparisons. A heritage top-hat thermal ion electrostatic analyzer, called the HT, and a multi-needle Langmuir probe have been used to validate both the PIPs and the data analysis method. Compared to the HT, the PIP ion temperature measurements agree with a root-mean-square error of 0.023 eV. These two instruments agree on the parallel-to-B plasma flow velocity with a root-mean-square error of 130 m/s. The PIP with its field of view aligned perpendicular-to-B provided a density measurement with an 11% error compared to the multi-needle Langmuir Probe. Higher error in the other PIP's density measurement is likely due to simplifications in the SPIS model geometry.

Structure of a truncated form of leucine zipper II of JIP3 reveals an unexpected antiparallel coiled-coil arrangement

JIP3 and JIP4, two highly related scaffolding proteins for MAP kinases, are binding partners for two molecular motors as well as for the small G protein ARF6. The leucine zipper II (LZII) region of JIP3/4 is the binding site for these three partners. Previously, the crystal structure of ARF6 bound to JIP4 revealed LZII in a parallel coiled-coil arrangement. Here, the crystal structure of an N-terminally truncated form of LZII of JIP3 alone shows an unexpected antiparallel arrangement. Using molecular dynamics and modelling, the stability of this antiparallel LZII arrangement, as well as its specificity for ARF6, were investigated. This study highlights that N-terminal truncation of LZII can change its coiled-coil orientation without affecting its overall stability. Further, a conserved buried asparagine residue was pinpointed as a possible structural determinant for this dramatic structural rearrangement. Thus, LZII of JIP3/4 is a versatile structural motif, modifications of which can impact partner recognition and thus biological function.