Conformationally Constrained Sialyl Analogues as New Potential Binders of h-CD22.

Here, two conformationally constrained sialyl analogues were synthesized and characterized in their interaction with the inhibitory Siglec, human CD22 (h-CD22). An orthogonal approach, including biophysical assays (SPR and fluorescence), ligand-based NMR techniques, and molecular modelling, was employed to disentangle the interaction mechanisms at a molecular level. The results showed that the Sialyl-TnThr antigen analogue represents a promising scaffold for the design of novel h-CD22 inhibitors. Our findings also suggest that the introduction of a biphenyl moiety at position 9 of the sialic acid hampers canonical accommodation of the ligand in the protein binding pocket, even though the affinity with respect to the natural ligand is increased. Our results address the search for novel modifications of the Neu5Ac-α(2-6)-Gal epitope, outline new insights for the design and synthesis of high-affinity h-CD22 ligands, and offer novel prospects for therapeutic intervention to prevent autoimmune diseases and B-cell malignancies.

SMARTCyp 3.0: enhanced cytochrome P450 site-of-metabolism prediction server.

MOTIVATION: Cytochromes P450 are the most important class of drug metabolizing enzymes. Prediction of drug metabolism is important in development of new drugs, to understand and reduce adverse drug reactions and to reduce animal testing. RESULTS: SMARTCyp 3.0 is an updated version of our previous web server for prediction of site-of-metabolism for Cytochrome P450-mediated metabolism, now in Python 3 with increased structural coverage and new features. The SMARTCyp program is a first principle-based method using density functional theory determined activation energies for more than 250 molecules to identify the most likely site-of-metabolism. New features include a similarity measure between the query molecule and the model fragment, a new graphical interface and additional parameters expanding the structural coverage of the SMARTCyp program. AVAILABILITY AND IMPLEMENTATION: The SMARTCyp server is freely available for use on the web at smartcyp.sund.ku.dk. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Impact of Mutations on NPAC Structural Dynamics: Mechanistic Insights from MD Simulations.

NPAC is a cytokine-like nuclear factor involved in chromatin modification and regulation of gene expression. In humans, the C-terminal domain of NPAC has the conserved structure of the ß-hydroxyacid dehydrogenases (ß-HAD) protein superfamily, which forms a stable tetrameric core scaffold for demethylase enzymes and organizes multiple sites for chromatin interactions. In spite of the close structural resemblance to other ß-HAD family members, the human NPAC dehydrogenase domain lacks a highly conserved catalytic lysine, substituted by a methionine. The reintroduction of the catalytic lysine by M437 K mutation results in a significant decrease of stability of the tetramer. Here, we have computationally investigated the molecular determinants of the functional differences between methionine and lysine-containing NPAC proteins. We find that the single mutation can determine strong consequences in terms of dynamics, stability, and ultimately ability to assemble in supramolecular complexes: the higher stability and lower flexibility of the methionine variant structurally preorganizes the monomer for tetramerization, whereas lysine increases flexibility and favors conformations that, while catalytically active, are not optimal for tetrameric assembly. We combine structure-dynamics analysis to an evolutionary study of NPAC sequences, showing that the methionine mutation occurs in a specifically flexible region of the lysine-containing protein, flanked by two domains that concentrate most of the stabilizing interactions. In our model, such separation of stability nuclei and flexible regions appears to favor the functional innovability of the protein.

A new calculation method for the free electron fraction of an ionization chamber in the ultra-high-dose-per-pulse regimen.

The free electron fraction is the fraction of electrons, produced inside the cavity of an ionization chamber after irradiation, which does not bind to gas molecules and thereby reaches the electrode as free electrons. It is a fundamental quantity to describe the recombination processes of an ionization chamber, as it generates a gap of positive charges compared to negative ones, which certainly will not undergo recombination. The free electron fraction depends on the specific chamber geometry, the polarizing applied voltage and the gas thermodynamic properties. Therefore, it is necessary to evaluate such fraction in an accurate and easy way for any measurement condition. In this paper, a simple and direct method for evaluating the free electron fraction of ionization chambers is proposed. We first model the capture process of the electrons produced inside an ionization chamber after the beam pulse; then we present a method to evaluate the free electron fraction based on simple measurements of collected charge, by varying the applied voltage. Finally, the results obtained using an Advanced Markus chamber irradiated with a Flash Radiotherapy dedicated research Linac (ElectronFlash) to estimate the free electron fraction are presented. The proposed method allows the use of a conventional ionization chamber for measurements in ultra-high-dose-per-pulse (UHDP) conditions, up to values of dose-per-pulse at which the perturbation of the electric field due to the generated charge can be considered negligible.

Behavior of glycolylated sialoglycans in the binding pockets of murine and human CD22.

Siglecs (sialic acid binding immunoglobulin (Ig)-like lectins) constitute a group of 15 human and 9 murine cell-surface transmembrane receptors belonging to the I-type lectin family, mostly expressed on innate immune cells and characterized by broadly similar structural features. Here, the prominent inhibitory CD22 (Siglec-2), well known in maintaining tolerance and preventing autoimmune responses on B cells, is studied in its human and murine forms in complex with sialoglycans. In detail, the role of the N-glycolyl neuraminic acid (Neu5Gc) moiety in the interaction with both orthologues was explored. The analysis of the binding mode was carried out by the combination of NMR spectroscopy, computational approaches, and CORCEMA-ST calculations. Our findings provide a first model of Neu5Gc recognition by h-CD22 and show a comparable molecular recognition profile by h- and m-CD22. These data open the way to innovative diagnostic and/or therapeutic methodologies to be used in the modulation of the immune responses.

Fast Methods for Prediction of Aldehyde Oxidase-Mediated Site-of-Metabolism.

Aldehyde Oxidase (AO) is an enzyme involved in the metabolism of aldehydes and N-containing heterocyclic compounds. Many drug compounds contain heterocyclic moieties, and AO metabolism has lead to failure of several late-stage drug candidates. Therefore, it is important to take AO-mediated metabolism into account early in the drug discovery process, and thus, to have fast and reliable models to predict the site of metabolism (SOM). We have collected a dataset of 78 substrates of human AO with a total of 89 SOMs and 347 non-SOMs and determined atomic descriptors for each compound. The descriptors comprise NMR shielding and ESP charges from density functional theory (DFT), NMR chemical shift from ChemBioDraw, and Gasteiger charges from RDKit. Additionally, atomic accessibility was considered using 2D-SASA and relative span descriptors from SMARTCyp. Finally, stability of the product, the metabolite, was determined with DFT and also used as a descriptor. All descriptors have AUC larger than 0.75. In particular, descriptors related to the chemical shielding and chemical shift (AUC = 0.96) and ESP charges (AUC = 0.96) proved to be good descriptors. We recommend two simple methods to identify the SOM for a given molecule: 1) use ChemBioDraw to calculate the chemical shift or 2) calculate ESP charges or chemical shift using DFT. The first approach is fast but somewhat difficult to automate, while the second is more time-consuming, but can easily be automated. The two methods predict correctly 93% and 91%, respectively, of the 89 experimentally observed SOMs.

Author Correction: The oxidative inactivation of FeFe hydrogenase reveals the flexibility of the H-cluster.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Potent Inhibitors against Newcastle Disease Virus Hemagglutinin-Neuraminidase.

Neuraminidase activity is essential for the infection and propagation of paramyxoviruses, including human parainfluenza viruses (hPIVs) and the Newcastle disease virus (NDV). Thus, many inhibitors have been developed based on the 2-deoxy-2,3-didehydro-d-N-acetylneuraminic acid inhibitor (DANA) backbone. Along this line, herein we report a series of neuraminidase inhibitors, having C4 (p-toluenesulfonamido and azido substituents) and C5 (N-perfluorinated chains) modifications to the DANA backbone, resulting in compounds with 5- to 15-fold greater potency than the currently most active compound, the N-trifluoroacetyl derivative of DANA (FANA), toward the NDV hemagglutinin-neuraminidase (NDV-HN). Remarkably, these inhibitors were found to be essentially inactive against the human sialidase NEU3, which is present on the outer layer of the cell membrane and is highly affected by the current NDV inhibitor FANA.

Aldehyde Oxidase: Reaction Mechanism and Prediction of Site of Metabolism.

Aldehyde oxidase (AO) is a molybdenum-containing enzyme involved in the clearance of drug compounds containing aldehydes and N-containing heterocyclic fragments. AO has gained considerable interest in recent years because of examples of too fast clearance of drug compounds in development. Thus, it is important to be able to predict AO-mediated drug metabolism. Therefore, we have characterized the structural and energetic aspects of different mechanisms with density functional theory using the molybdenum cofactor as a model for the reactive part of the enzyme. For a series of 6-substituted 4-quinazolinones, the trend in activation energies is the same for three tested reaction mechanisms. Using the concerted mechanism as a model for the enzymatic reaction, the transition states (TSs) for the formation of all possible metabolites for a series of known AO substrates were determined. The lowest activation energies correspond in all cases to the experimentally observed sites of metabolism (SOMs). Various molecular properties were calculated and investigated as more easily determinable markers for reactivity. The stabilities of both intermediates and products correlate to some extent with the TS energies and may be used to predict the SOM. The electrostatic-potential-derived charges are also good markers for the prediction of the experimental SOM for this set of compounds and may pave the way for the development of fast methods for the prediction of SOM for AO substrates.

The CYP79A1 catalyzed conversion of tyrosine to (E)-p-hydroxyphenylacetaldoxime unravelled using an improved method for homology modeling.

The vast diversity and membrane-bound nature of plant P450s makes it challenging to study the structural characteristics of this class of enzymes especially with respect to accurate intermolecular enzyme-substrate interactions. To address this problem we here apply a modified hybrid structure strategy for homology modeling of plant P450s. This allows for structural elucidation based on conserved motifs in the protein sequence and secondary structure predictions. We modeled the well-studied Sorghum bicolor cytochrome P450 CYP79A1 catalyzing the first step in the biosynthesis of the cyanogenic glucoside dhurrin. Docking experiments identified key regions of the active site involved in binding of the substrate and facilitating catalysis. Arginine 152 and threonine 534 were identified as key residues interacting with the substrate. The model was validated experimentally using site-directed mutagenesis. The new CYP79A1 model provides detailed insights into the mechanism of the initial steps in cyanogenic glycoside biosynthesis. The approach could guide functional characterization of other membrane-bound P450s and provide structural guidelines for elucidation of key structure-function relationships of other plant P450s.

The oxidative inactivation of FeFe hydrogenase reveals the flexibility of the H-cluster.

Nature is a valuable source of inspiration in the design of catalysts, and various approaches are used to elucidate the mechanism of hydrogenases, the enzymes that oxidize or produce H2. In FeFe hydrogenases, H2 oxidation occurs at the H-cluster, and catalysis involves H2 binding on the vacant coordination site of an iron centre. Here, we show that the reversible oxidative inactivation of this enzyme results from the binding of H2 to coordination positions that are normally blocked by intrinsic CO ligands. This flexibility of the coordination sphere around the reactive iron centre confers on the enzyme the ability to avoid harmful reactions under oxidizing conditions, including exposure to O2. The versatile chemistry of the diiron cluster in the natural system might inspire the design of novel synthetic catalysts for H2 oxidation.