C. David Andersson

Determinants for Simultaneous Binding of Copper and Platinum to Human Chaperone Atox1: Hitchhiking not Hijacking

Cisplatin (CisPt) is an anticancer agent that has been used for decades to treat a variety of cancers. CisPt treatment causes many side effects due to interactions with proteins that detoxify the drug before reaching the DNA. One key player in CisPt resistance is the cellular copper-transport system involving the uptake protein Ctr1, the cytoplasmic chaperone Atox1 and the secretory path ATP7A/B proteins. CisPt has been shown to bind to ATP7B, resulting in vesicle sequestering of the drug. In addition, we and others showed that the apo-form of Atox1 could interact with CisPt in vitro and in vivo. Since the function of Atox1 is to transport copper (Cu) ions, it is important to assess how CisPt binding depends on Cu-loading of Atox1. Surprisingly, we recently found that CisPt interacted with Cu-loaded Atox1 in vitro at a position near the Cu site such that unique spectroscopic features appeared. Here, we identify the binding site for CisPt in the Cu-loaded form of Atox1 using strategic variants and a combination of spectroscopic and chromatographic methods. We directly prove that both metals can bind simultaneously and that the unique spectroscopic signals originate from an Atox1 monomer species. Both Cys in the Cu-site (Cys12, Cys15) are needed to form the di-metal complex, but not Cys41. Removing Met10 in the conserved metal-binding motif makes the loop more floppy and, despite metal binding, there are no metal-metal electronic transitions. In silico geometry minimizations provide an energetically favorable model of a tentative ternary Cu-Pt-Atox1 complex. Finally, we demonstrate that Atox1 can deliver CisPt to the fourth metal binding domain 4 of ATP7B (WD4), indicative of a possible drug detoxification mechanism.

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.

Discovery of Ligands for ADP-Ribosyltransferases via Docking-Based Virtual Screening.

The diphtheria toxin-like ADP-ribosyltransferases (ARTDs) are an enzyme family that catalyzes the transfer of ADP-ribose units onto substrate proteins by using nicotinamide adenine dinucleotide (NAD(+)) as a cosubstrate. They have a documented role in chromatin remodelling and DNA repair, and inhibitors of ARTD1 and 2 (PARP1 and 2) are currently in clinical trials for the treatment of cancer. The detailed function of most other ARTDs is still unknown. By using virtual screening, we identified small ligands of ARTD7 (PARP15/BAL3) and ARTD8 (PARP14/BAL2). Thermal-shift assays confirmed that 16 compounds, belonging to eight structural classes, bound to ARTD7/ARTD8. Affinity measurements with isothermal titration calorimetry for two isomers of the most promising hit compound confirmed binding in the low micromolar range to ARTD8. Crystal structures showed anchoring of the hits in the nicotinamide pocket. These results form a starting point in the development of chemical tools for the study of the role and function of ARTD7 and ARTD8.

Targeting Acetylcholinesterase: Identification of Chemical Leads by High Throughput Screening, Structure Determination and Molecular Modeling

Acetylcholinesterase (AChE) is an essential enzyme that terminates cholinergic transmission by rapid hydrolysis of the neurotransmitter acetylcholine. Compounds inhibiting this enzyme can be used (inter alia) to treat cholinergic deficiencies (e.g. in Alzheimers disease), but may also act as dangerous toxins (e.g. nerve agents such as sarin). Treatment of nerve agent poisoning involves use of antidotes, small molecules capable of reactivating AChE. We have screened a collection of organic molecules to assess their ability to inhibit the enzymatic activity of AChE, aiming to find lead compounds for further optimization leading to drugs with increased efficacy and/or decreased side effects. 124 inhibitors were discovered, with considerable chemical diversity regarding size, polarity, flexibility and charge distribution. An extensive structure determination campaign resulted in a set of crystal structures of protein-ligand complexes. Overall, the ligands have substantial interactions with the peripheral anionic site of AChE, and the majority form additional interactions with the catalytic site (CAS). Reproduction of the bioactive conformation of six of the ligands using molecular docking simulations required modification of the default parameter settings of the docking software. The results show that docking-assisted structure-based design of AChE inhibitors is challenging and requires crystallographic support to obtain reliable results, at least with currently available software. The complex formed between C5685 and Mus musculus AChE (C5685•mAChE) is a representative structure for the general binding mode of the determined structures. The CAS binding part of C5685 could not be structurally determined due to a disordered electron density map and the developed docking protocol was used to predict the binding modes of this part of the molecule. We believe that chemical modifications of our discovered inhibitors, biochemical and biophysical characterization, crystallography and computational chemistry provide a route to novel AChE inhibitors and reactivators.

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. Proteins 2010.

Towards small molecule inhibitors of mono-ADP-ribosyltransferases

Protein ADP-ribosylation is a post-translational modification involved in DNA repair, protein degradation, transcription regulation, and epigenetic events. Intracellular ADP-ribosylation is catalyzed predominantly by ADP-ribosyltransferases with diphtheria toxin homology (ARTDs). The most prominent member of the ARTD family, poly(ADP-ribose) polymerase-1 (ARTD1/PARP1) has been a target for cancer drug development for decades. Current PARP inhibitors are generally non-selective, and inhibit the mono-ADP-ribosyltransferases with low potency. Here we describe the synthesis of acylated amino benzamides and screening against the mono-ADP-ribosyltransferases ARTD7/PARP15, ARTD8/PARP14, ARTD10/PARP10, and the poly-ADP-ribosyltransferase ARTD1/PARP1. The most potent compound inhibits ARTD10 with sub-micromolar IC50.

Design of Glycopeptides Used to Investigate Class II MHC Binding and T-Cell Responses Associated with Autoimmune Arthritis

The glycopeptide fragment CII259–273 from type II collagen (CII) binds to the murine Aq and human DR4 class II Major Histocompatibility Complex (MHC II) proteins, which are associated with development of murine collagen-induced arthritis (CIA) and rheumatoid arthritis (RA), respectively. It has been shown that CII259–273 can be used in therapeutic vaccination of CIA. This glycopeptide also elicits responses from T-cells obtained from RA patients, which indicates that it has an important role in RA as well. We now present a methodology for studies of (glyco)peptide-receptor interactions based on a combination of structure-based virtual screening, ligand-based statistical molecular design and biological evaluations. This methodology included the design of a CII259–273 glycopeptide library in which two anchor positions crucial for binding in pockets of Aq and DR4 were varied. Synthesis and biological evaluation of the designed glycopeptides provided novel structure-activity relationship (SAR) understanding of binding to Aq and DR4. Glycopeptides that retained high affinities for these MHC II proteins and induced strong responses in panels of T-cell hybridomas were also identified. An analysis of all the responses revealed groups of glycopeptides with different response patterns that are of high interest for vaccination studies in CIA. Moreover, the SAR understanding obtained in this study provides a platform for the design of second-generation glycopeptides with tuned MHC affinities and T-cell responses.

Similar but Different: Thermodynamic and Structural Characterization of a Pair of Enantiomers Binding to Acetylcholinesterase†

Take a closer look: Unexpectedly, a pair of enantiomeric ligands proved to have similar binding affinities for acetylcholinesterase. Further studies indicated that the enantiomers exhibit different ...

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.

Structural Basis for Lack of ADP-ribosyltransferase Activity in Poly(ADP-ribose) Polymerase-13/Zinc Finger Antiviral Protein

Background: PARP13 contains a divergent PARP homology ADP-ribosyltransferase domain of unknown function. Results: The consensus NAD+ pocket of PARP13 is occluded by interacting protein side chains. Conclusion: PARP13 lacks the structural requirements for NAD+ binding. Significance: Evolutionary conservation of enzymatic inactivity suggests a need for a rigid domain structure. 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.