Structure-function study of a Ca2+-independent metacaspase involved in lateral root emergence.

Metacaspases are part of an evolutionarily broad family of multifunctional cysteine proteases, involved in disease and normal development. As the structure-function relationship of metacaspases remains poorly understood, we solved the X-ray crystal structure of an Arabidopsis thaliana type II metacaspase (AtMCA-IIf) belonging to a particular subgroup not requiring calcium ions for activation. To study metacaspase activity in plants, we developed an in vitro chemical screen to identify small molecule metacaspase inhibitors and found several hits with a minimal thioxodihydropyrimidine-dione structure, of which some are specific AtMCA-IIf inhibitors. We provide mechanistic insight into the basis of inhibition by the TDP-containing compounds through molecular docking onto the AtMCA-IIf crystal structure. Finally, a TDP-containing compound (TDP6) effectively hampered lateral root emergence in vivo, probably through inhibition of metacaspases specifically expressed in the endodermal cells overlying developing lateral root primordia. In the future, the small compound inhibitors and crystal structure of AtMCA-IIf can be used to study metacaspases in other species, such as important human pathogens, including those causing neglected diseases.

Broad-Spectrum Antidote Discovery by Untangling the Reactivation Mechanism of Nerve-Agent-Inhibited Acetylcholinesterase.

Reactivators are vital for the treatment of organophosphorus nerve agent (OPNA) intoxication but new alternatives are needed due to their limited clinical applicability. The toxicity of OPNAs stems from covalent inhibition of the essential enzyme acetylcholinesterase (AChE), which reactivators relieve via a chemical reaction with the inactivated enzyme. Here, we present new strategies and tools for developing reactivators. We discover suitable inhibitor scaffolds by using an activity-independent competition assay to study non-covalent interactions with OPNA-AChEs and transform these inhibitors into broad-spectrum reactivators. Moreover, we identify determinants of reactivation efficiency by analysing reactivation and pre-reactivation kinetics together with structural data. Our results show that new OPNA reactivators can be discovered rationally by exploiting detailed knowledge of the reactivation mechanism of OPNA-inhibited AChE.

In Situ Assembly of Choline Acetyltransferase Ligands by a Hydrothiolation Reaction Reveals Key Determinants for Inhibitor Design.

The potential drug target choline acetyltransferase (ChAT) catalyses the production of the neurotransmitter acetylcholine in cholinergic neurons, T-cells, and B-cells. Herein, we show that arylvinylpyridiniums (AVPs), the most widely studied class of ChAT inhibitors, act as substrate in an unusual coenzyme A-dependent hydrothiolation reaction. This in situ synthesis yields an adduct that is the actual enzyme inhibitor. The adduct is deeply buried in the active site tunnel of ChAT and interactions with a hydrophobic pocket near the choline binding site have major implications for the molecular recognition of inhibitors. Our findings clarify the inhibition mechanism of AVPs, establish a drug modality that exploits a target-catalysed reaction between exogenous and endogenous precursors, and provide new directions for the development of ChAT inhibitors with improved potency and bioactivity.

Nucleation of an Activating Conformational Change by a Cation-π Interaction.

As a key molecule in biology, adenosine triphosphate (ATP) has numerous crucial functions in, for instance, energetics, post-translational modifications, nucleotide biosynthesis, and cofactor metabolism. Here, we have discovered an intricate interplay between the enzyme adenylate kinase and its substrate ATP. The side chain of an arginine residue was found to be an efficient sensor of the aromatic moiety of ATP through the formation of a strong cation-π interaction. In addition to recognition, the interaction was found to have dual functionality. First, it nucleates the activating conformational transition of the ATP binding domain and also affects the specificity in the distant AMP binding domain. In light of the functional consequences resulting from the cation-π interaction, it is possible that the mode of ATP recognition may be a useful tool in enzyme design.

Structure of a prereaction complex between the nerve agent sarin, its biological target acetylcholinesterase, and the antidote HI-6.

Organophosphorus nerve agents interfere with cholinergic signaling by covalently binding to the active site of the enzyme acetylcholinesterase (AChE). This inhibition causes an accumulation of the neurotransmitter acetylcholine, potentially leading to overstimulation of the nervous system and death. Current treatments include the use of antidotes that promote the release of functional AChE by an unknown reactivation mechanism. We have used diffusion trap cryocrystallography and density functional theory (DFT) calculations to determine and analyze prereaction conformers of the nerve agent antidote HI-6 in complex with Mus musculus AChE covalently inhibited by the nerve agent sarin. These analyses reveal previously unknown conformations of the system and suggest that the cleavage of the covalent enzyme-sarin bond is preceded by a conformational change in the sarin adduct itself. Together with data from the reactivation kinetics, this alternate conformation suggests a key interaction between Glu202 and the O-isopropyl moiety of sarin. Moreover, solvent kinetic isotope effect experiments using deuterium oxide reveal that the reactivation mechanism features an isotope-sensitive step. These findings provide insights into the reactivation mechanism and provide a starting point for the development of improved antidotes. The work also illustrates how DFT calculations can guide the interpretation, analysis, and validation of crystallographic data for challenging reactive systems with complex conformational dynamics.

Retention-time prediction in comprehensive two-dimensional gas chromatography to aid identification of unknown contaminants.

Comprehensive two-dimensional (2D) gas chromatography (GC×GC) coupled to mass spectrometry (MS, GC×GC-MS), which enhances selectivity compared to GC-MS analysis, can be used for non-directed analysis (non-target screening) of environmental samples. Additional tools that aid in identifying unknown compounds are needed to handle the large amount of data generated. These tools include retention indices for characterizing relative retention of compounds and prediction of such. In this study, two quantitative structure-retention relationship (QSRR) approaches for prediction of retention times (1tR and 2tR) and indices (linear retention indices (LRIs) and a new polyethylene glycol-based retention index (PEG-2I)) in GC × GC were explored, and their predictive power compared. In the first method, molecular descriptors combined with partial least squares (PLS) analysis were used to predict times and indices. In the second method, the commercial software package ChromGenius (ACD/Labs), based on a "federation of local models," was employed. Overall, the PLS approach exhibited better accuracy than the ChromGenius approach. Although average errors for the LRI prediction via ChromGenius were slightly lower, PLS was superior in all other cases. The average deviations between the predicted and the experimental value were 5% and 3% for the 1tR and LRI, and 5% and 12% for the 2tR and PEG-2I, respectively. These results are comparable to or better than those reported in previous studies. Finally, the developed model was successfully applied to an independent dataset and led to the discovery of 12 wrongly assigned compounds. The results of the present work represent the first-ever prediction of the PEG-2I. Graphical abstract ᅟ.

An Unusual Dimeric Inhibitor of Acetylcholinesterase: Cooperative Binding of Crystal Violet.

Acetylcholinesterase (AChE) is an essential enzyme that terminates cholinergic transmission by a rapid hydrolysis of the neurotransmitter acetylcholine. AChE is an important target for treatment of various cholinergic deficiencies, including Alzheimer's disease and myasthenia gravis. In a previous high throughput screening campaign, we identified the dye crystal violet (CV) as an inhibitor of AChE. Herein, we show that CV displays a significant cooperativity for binding to AChE, and the molecular basis for this observation has been investigated by X-ray crystallography. Two monomers of CV bind to residues at the entrance of the active site gorge of the enzyme. Notably, the two CV molecules have extensive intermolecular contacts with each other and with AChE. Computational analyses show that the observed CV dimer is not stable in solution, suggesting the sequential binding of two monomers. Guided by the structural analysis, we designed a set of single site substitutions, and investigated their effect on the binding of CV. Only moderate effects on the binding and the cooperativity were observed, suggesting a robustness in the interaction between CV and AChE. Taken together, we propose that the dimeric cooperative binding is due to a rare combination of chemical and structural properties of both CV and the AChE molecule itself.

Structural basis for lack of ADP-ribosyltransferase activity in poly(ADP-ribose) polymerase-13/zinc finger antiviral protein.

The mammalian poly(ADP-ribose) polymerase (PARP) family includes ADP-ribosyltransferases with diphtheria toxin homology (ARTD). Most members have mono-ADP-ribosyltransferase activity. PARP13/ARTD13, also called zinc finger antiviral protein, has roles in viral immunity and microRNA-mediated stress responses. PARP13 features a divergent PARP homology domain missing a PARP consensus sequence motif; the domain has enigmatic functions and apparently lacks catalytic activity. We used x-ray crystallography, molecular dynamics simulations, and biochemical analyses to investigate the structural requirements for ADP-ribosyltransferase activity in human PARP13 and two of its functional partners in stress granules: PARP12/ARTD12, and PARP15/BAL3/ARTD7. The crystal structure of the PARP homology domain of PARP13 shows obstruction of the canonical active site, precluding NAD(+) binding. Molecular dynamics simulations indicate that this closed cleft conformation is maintained in solution. Introducing consensus side chains in PARP13 did not result in 3-aminobenzamide binding, but in further closure of the site. Three-dimensional alignment of the PARP homology domains of PARP13, PARP12, and PARP15 illustrates placement of PARP13 residues that deviate from the PARP family consensus. Introducing either one of two of these side chains into the corresponding positions in PARP15 abolished PARP15 ADP-ribosyltransferase activity. Taken together, our results show that PARP13 lacks the structural requirements for ADP-ribosyltransferase activity.

The Nature of Activated Non-classical Hydrogen Bonds: A Case Study on Acetylcholinesterase-Ligand Complexes.

Molecular recognition events in biological systems are driven by non-covalent interactions between interacting species. Here, we have studied hydrogen bonds of the CH⋅⋅⋅Y type involving electron-deficient CH donors using dispersion-corrected density functional theory (DFT) calculations applied to acetylcholinesterase-ligand complexes. The strengths of CH⋅⋅⋅Y interactions activated by a proximal cation were considerably strong; comparable to or greater than those of classical hydrogen bonds. Significant differences in the energetic components compared to classical hydrogen bonds and non-activated CH⋅⋅⋅Y interactions were observed. Comparison between DFT and molecular mechanics calculations showed that common force fields could not reproduce the interaction energy values of the studied hydrogen bonds. The presented results highlight the importance of considering CH⋅⋅⋅Y interactions when analysing protein-ligand complexes, call for a review of current force fields, and opens up possibilities for the development of improved design tools for drug discovery.

Hydroxyethylene isosteres introduced in type II collagen fragments substantially alter the structure and dynamics of class II MHC A(q)/glycopeptide complexes.

Class II major histocompatibility complex (MHC) proteins are involved in initiation of immune responses to foreign antigens via presentation of peptides to receptors of CD4(+) T-cells. An analogous presentation of self-peptides may lead to autoimmune diseases, such as rheumatoid arthritis (RA). The glycopeptide fragment CII259-273, derived from type II collagen, is presented by A(q) MHCII molecules in the mouse and has a key role in development of collagen induced arthritis (CIA), a validated model for RA. We have introduced hydroxyethylene amide bond isosteres at the Ala(261)-Gly(262) position of CII259-273. Biological evaluation showed that A(q) binding and T cell recognition were dramatically reduced for the modified glycopeptides, although static models predicted similar binding modes as the native type II collagen fragment. Molecular dynamics (MD) simulations demonstrated that introduction of the hydroxyethylene isosteres disturbed the entire hydrogen bond network between the glycopeptides and A(q). As a consequence the hydroxyethylene isosteric glycopeptides were prone to dissociation from A(q) and unfolding of the ß1-helix. Thus, the isostere induced adjustment of the hydrogen bond network altered the structure and dynamics of A(q)/glycopeptide complexes leading to the loss of A(q) affinity and subsequent T cell response.

Benefits of statistical molecular design, covariance analysis, and reference models in QSAR: a case study on acetylcholinesterase.

Scientific disciplines such as medicinal- and environmental chemistry, pharmacology, and toxicology deal with the questions related to the effects small organic compounds exhort on biological targets and the compounds' physicochemical properties responsible for these effects. A common strategy in this endeavor is to establish structure-activity relationships (SARs). The aim of this work was to illustrate benefits of performing a statistical molecular design (SMD) and proper statistical analysis of the molecules' properties before SAR and quantitative structure-activity relationship (QSAR) analysis. Our SMD followed by synthesis yielded a set of inhibitors of the enzyme acetylcholinesterase (AChE) that had very few inherent dependencies between the substructures in the molecules. If such dependencies exist, they cause severe errors in SAR interpretation and predictions by QSAR-models, and leave a set of molecules less suitable for future decision-making. In our study, SAR- and QSAR models could show which molecular sub-structures and physicochemical features that were advantageous for the AChE inhibition. Finally, the QSAR model was used for the prediction of the inhibition of AChE by an external prediction set of molecules. The accuracy of these predictions was asserted by statistical significance tests and by comparisons to simple but relevant reference models.

Structure-Activity Relationships Reveal Beneficial Selectivity Profiles of Inhibitors Targeting Acetylcholinesterase of Disease-Transmitting Mosquitoes.

Insecticide resistance jeopardizes the prevention of infectious diseases such as malaria and dengue fever by vector control of disease-transmitting mosquitoes. Effective new insecticidal compounds with minimal adverse effects on humans and the environment are therefore urgently needed. Here, we explore noncovalent inhibitors of the well-validated insecticidal target acetylcholinesterase (AChE) based on a 4-thiazolidinone scaffold. The 4-thiazolidinones inhibit AChE1 from the mosquitoes Anopheles gambiae and Aedes aegypti at low micromolar concentrations. Their selectivity depends primarily on the substitution pattern of the phenyl ring; halogen substituents have complex effects. The compounds also feature a pendant aliphatic amine that was important for activity; little variation of this group is tolerated. Molecular docking studies suggested that the tight selectivity profiles of these compounds are due to competition between two binding sites. Three 4-thiazolidinones tested for in vivo insecticidal activity had similar effects on disease-transmitting mosquitoes despite a 10-fold difference in their in vitro activity.

Bioisosteric replacement based on 1,2,4-oxadiazoles in the discovery of 1H-indazole-bearing neuroprotective MAO B inhibitors.

Following a hybridization strategy, a series of 5-substituted-1H-indazoles were designed and evaluated in vitro as inhibitors of human monoamine oxidase (hMAO) A and B. Among structural modifications, the bioisostere-based introduction of 1,2,4-oxadiazole ring returned the most potent and selective human MAO B inhibitor (compound 20, IC50 = 52 nM, SI > 192). The most promising inhibitors were studied in cell-based neuroprotection models of SH-SY5Y and astrocytes line against H2O2. Moreover, preliminary drug-like features (aqueous solubility at pH 7.4; hydrolytic stability at acidic and neutral pH) were assessed for selected 1,2,4-oxadiazoles and compared to amide analogues through RP-HPLC methods. Molecular docking simulations highlighted the crucial role of molecular flexibility in providing a better shape complementarity for compound 20 within MAO B enzymatic cleft than rigid analogue 18. Enzymatic kinetics analysis along with thermal stability curves (Tm shift = +2.9 °C) provided clues of a tight-binding mechanism for hMAO B inhibition by 20.

Mapping pilicide anti-virulence effect in Escherichia coli, a comprehensive structure-activity study.

Pilicides prevent pili formation and thereby the development of bacterial biofilms in Escherichia coli. We have performed a comprehensive structure activity relationship (SAR) study of the dihydrothiazolo ring-fused 2-pyridone pilicide central fragment by varying all open positions. Orthogonal projections to latent structures discriminant analysis (OPLS-DA) was used to distinguish active from inactive compounds in which polarity proved to be the most important factor for discrimination. A quantitative SAR (QSAR) partial least squares (PLS) model was calculated on the active compounds for prediction of biofilm inhibition activity. In this model, compounds with high inhibitory activity were generally larger, more lipophilic, more flexible and had a lower HOMO. Overall, this resulted in both highly valuable SAR information and potent inhibitors of type 1 pili dependent biofilm formation. The most potent biofilm inhibitor had an EC(50) of 400 nM.

Probing Fluorinated Motifs onto Dual AChE-MAO B Inhibitors: Rational Design, Synthesis, Biological Evaluation, and Early-ADME Studies.

Bioisosteric H/F or CH2OH/CF2H replacement was introduced in coumarin derivatives previously characterized as dual AChE-MAO B inhibitors to probe the effects on both inhibitory potency and drug-likeness. Along with in vitro screening, we investigated early-ADME parameters related to solubility and lipophilicity (Sol7.4, CHI7.4, log D7.4), oral bioavailability and central nervous system (CNS) penetration (PAMPA-HDM and PAMPA-blood-brain barrier (BBB) assays, Caco-2 bidirectional transport study), and metabolic liability (half-lives and clearance in microsomes, inhibition of CYP3A4). Both specific and nonspecific tissue toxicities were determined in SH-SY5Y and HepG2 lines, respectively. Compound 15 bearing a -CF2H motif emerged as a water-soluble, orally bioavailable CNS-permeant potent inhibitor of both human AChE (IC50 = 550 nM) and MAO B (IC50 = 8.2 nM, B/A selectivity > 1200). Moreover, 15 behaved as a safe and metabolically stable neuroprotective agent, devoid of cytochrome liability.

(E)-alkene and ethylene isosteres substantially alter the hydrogen-bonding network in class II MHC A(q)/glycopeptide complexes and affect T-cell recognition.

The structural basis for antigen presentation by class II major histocompatibility complex (MHC) proteins to CD4(+) T-cells is important for understanding and possibly treating autoimmune diseases. In the work described in this paper, (E)-alkene and ethylene amide-bond isosteres were used to investigate the effect of removing hydrogen-bonding possibilities from the CII259-270 glycopeptide, which is bound by the arthritis-associated murine A(q) class II MHC protein. The isostere-modified glycopeptides showed varying and unexpectedly large losses of A(q) binding that could be linked to the dynamics of the system. Molecular dynamics (MD) simulations revealed that the backbone of CII259-270 and the A(q) protein are able to form up to 11 hydrogen bonds, but fewer than this number are present at any one time. Most of the strong hydrogen-bond interactions were formed by the N-terminal part of the glycopeptide, i.e., in the region where the isosteric replacements were made. The structural dynamics also revealed that hydrogen bonds were strongly coupled to each other; the loss of one hydrogen-bond interaction had a profound effect on the entire hydrogen-bonding network. The A(q) binding data revealed that an ethylene isostere glycopeptide unexpectedly bound more strongly to A(q) than the corresponding (E)-alkene, which is in contrast to the trend observed for the other isosteres. Analysis of the MD trajectories revealed that the complex conformation of this ethylene isostere was structurally different and had an altered molecular interaction pattern compared to the other A(q)/glycopeptide complexes. The introduced amide-bond isosteres also affected the interactions of the glycopeptide/A(q) complexes with T-cell receptors. The dynamic variation of the patterns and strengths of the hydrogen-bond interactions in the class II MHC system is of critical importance for the class II MHC/peptide/TCR signaling system.

Postprocessing of docked protein-ligand complexes using implicit solvation models.

Molecular docking plays an important role in drug discovery as a tool for the structure-based design of small organic ligands for macromolecules. Possible applications of docking are identification of the bioactive conformation of a protein-ligand complex and the ranking of different ligands with respect to their strength of binding to a particular target. We have investigated the effect of implicit water on the postprocessing of binding poses generated by molecular docking using MM-PB/GB-SA (molecular mechanics Poisson-Boltzmann and generalized Born surface area) methodology. The investigation was divided into three parts: geometry optimization, pose selection, and estimation of the relative binding energies of docked protein-ligand complexes. Appropriate geometry optimization afforded more accurate binding poses for 20% of the complexes investigated. The time required for this step was greatly reduced by minimizing the energy of the binding site using GB solvation models rather than minimizing the entire complex using the PB model. By optimizing the geometries of docking poses using the GB(HCT+SA) model then calculating their free energies of binding using the PB implicit solvent model, binding poses similar to those observed in crystal structures were obtained. Rescoring of these poses according to their calculated binding energies resulted in improved correlations with experimental binding data. These correlations could be further improved by applying the postprocessing to several of the most highly ranked poses rather than focusing exclusively on the top-scored pose. The postprocessing protocol was successfully applied to the analysis of a set of Factor Xa inhibitors and a set of glycopeptide ligands for the class II major histocompatibility complex (MHC) A(q) protein. These results indicate that the protocol for the postprocessing of docked protein-ligand complexes developed in this paper may be generally useful for structure-based design in drug discovery.

Mapping of ligand-binding cavities in proteins.

The complex interactions between proteins and small organic molecules (ligands) are intensively studied because they play key roles in biological processes and drug activities. Here, we present a novel approach to characterize and map the ligand-binding cavities of proteins without direct geometric comparison of structures, based on Principal Component Analysis of cavity properties (related mainly to size, polarity, and charge). This approach can provide valuable information on the similarities and dissimilarities, of binding cavities due to mutations, between-species differences and flexibility upon ligand-binding. The presented results show that information on ligand-binding cavity variations can complement information on protein similarity obtained from sequence comparisons. The predictive aspect of the method is exemplified by successful predictions of serine proteases that were not included in the model construction. The presented strategy to compare ligand-binding cavities of related and unrelated proteins has many potential applications within protein and medicinal chemistry, for example in the characterization and mapping of "orphan structures", selection of protein structures for docking studies in structure-based design, and identification of proteins for selectivity screens in drug design programs.

Oxazole-modified glycopeptides that target arthritis-associated class II MHC A(q) and DR4 proteins.

The glycopeptide CII259-273, a fragment from type II collagen (CII), can induce tolerance in mice susceptible to collagen-induced arthritis (CIA), which is a validated disease model for rheumatoid arthritis (RA). Here, we describe the design and synthesis of a small series of modified CII259-273 glycopeptides with oxazole heterocycles replacing three potentially labile peptide bonds. These glycopeptidomimetics were evaluated for binding to murine CIA-associated A(q) and human RA-associated DR4 class II major histocompatibility complex (MHC) proteins. The oxazole modifications drastically reduced or completely abolished binding to A(q). Two of the glycopeptidomimetics were, however, well tolerated in binding to DR4 and they also induced strong responses by one or two DR4-restricted T-cell hybridomas. This work contributes to the development of an altered glycopeptide for inducing immunological tolerance in CIA, with the long-term goal of developing a therapeutic vaccine for treatment of RA.