Exploring the pH- and Ligand-Dependent Flap Dynamics of Malarial Plasmepsin II

Design, synthesis and modelling of photoreactive chemical probes for investigating target engagement of plasmepsin IX and X in Plasmodium falciparum

The emergence of Plasmodium parasite resistance to current front-line antimalarial treatments poses a serious threat to global malaria control and highlights the necessity for the development of therapeutics with novel targets and mechanisms of action. Plasmepsins IX and X (PMIX/PMX) have been recognised as highly promising targets in Plasmodium due to their contribution to parasite's pathogenicity. Recent research has demonstrated that dual PMIX/PMX inhibition results in the impairment of multiple parasite's life cycle stages, which is an important feature in drug resistance prevention. Herein we report novel hydroxyethylamine photoaffinity labelling (PAL) probes, designed for PMIX/PMX target engagement and proteomics experiments in Plasmodium parasites. The prepared probes have both a photoreactive group (diazirine or benzophenone) for covalent attachment to target proteins, and a terminal alkyne handle allowing their use in bioorthogonal ligation. One of the synthesised benzophenone probes was shown to be highly promising as demonstrated by its outstanding antimalarial potency (IC50 = 15 nM versus D10 P. falciparum) and its inhibitory effect against PfPMX in an enzymatic assay. Molecular docking and molecular dynamics studies show that the inclusion of the benzophenone and alkyne handle does not alter the binding mode compared to the parent compound. The photoaffinity probe can be used in future chemical proteomics studies to allow hydroxyethylamine drug scaffold target identification and validation in Plasmodium. We expect our findings to act as a tool for future investigations on PMIX/PMX inhibition in antimalarial drug discovery.

PKAD-2: New entries and expansion of functionalities of the database of experimentally measured pKa's of proteins

Almost all biological reactions are pH dependent and understanding the origin of pH dependence requires knowledge of the pKa's of ionizable groups. Here we report a new edition of PKAD, the PKAD-2, which is a database of experimentally measured pKa's of proteins, both wild type and mutant proteins. The new additions include 117 wild type and 54 mutant pKa values, resulting in total 1742 experimentally measured pKa's. The new edition of PKAD-2 includes 8 new wild type and 12 new mutant proteins, resulting in total of 220 proteins. This new edition incorporates a visual 3D image of the highlighted residue of interest within the corresponding protein or protein complex. Hydrogen bonds were identified, counted, and implemented as a search feature. Other new search features include the number of neighboring residues <4A from the heaviest atom of the side chain of a given amino acid. Here, we present PKAD-2 with the intention to continuously incorporate novel features and current data with the goal to be used as benchmark for computational methods.

Elucidating the Protonation State of the γ-Secretase Catalytic Dyad

γ-Secretase (GS) is an intramembrane aspartyl protease that participates in the sequential cleavage of C99 to generate different isoforms of the amyloid-β (Aβ) peptides that are associated with the development of Alzheimer's disease. Due to its importance in the proteolytic processing of C99 by GS, we performed pH replica exchange molecular dynamics (pH-REMD) simulations of GS in its apo and substrate-bound forms to sample the protonation states of the catalytic dyad. We found that the catalytic dyad is deprotonated at physiological pH in our apo form, but the presence of the substrate at the active site displaces its monoprotonated state toward physiological pH. Our results show that Asp257 acts as the general base and Asp385 as the general acid during the cleavage mechanism. We identified different amino acids such as Lys265, Arg269, and the PAL motif interacting with the catalytic dyad and promoting changes in its acid-base behavior. Finally, we also found a significant pK_a shift of Glu280 related to the internalization of TM6-CT in the GS-apo form. Our study provides critical mechanistic insight into the GS mechanism and the basis for future research on the genesis of Aβ peptides and the development of Alzheimer's disease.

Constant pH molecular dynamics simulations: Current status and recent applications

Many important protein functions are carried out through proton-coupled conformational dynamics. Thus, the ability to accurately model protonation states dynamically has wide-ranging implications. Over the past two decades, two main types of constant pH methods (discrete and continuous) have been developed to enable proton-coupled molecular dynamics (MD) simulations. In this short review, we discuss the current status of the development and highlight recent applications that have advanced our understanding of protein structure-function relationships. We conclude the review by outlining the remaining challenges in the method development and projecting important areas for future applications.

A Guide to the Continuous Constant pH Molecular Dynamics Methods in Amber and CHARMM [Article v1.0]

Like temperature and pressure, solution pH is an important environmental variable in biomolecular simulations. Virtually all proteins depend on pH to maintain their structure and function. In conventional molecular dynamics (MD) simulations of proteins, pH is implicitly accounted for by assigning and fixing protonation states of titratable sidechains. This is a significant limitation, as the assigned protonation states may be wrong and they may change during dynamics. In this tutorial, we guide the reader in learning and using the various continuous constant pH MD methods in Amber and CHARMM packages, which have been applied to predict pK a values and elucidate proton-coupled conformational dynamics of a variety of proteins including enzymes and membrane transporters.

The truncated AaActin1 promoter is a candidate tool for metabolic engineering of artemisinin biosynthesis in Artemisia annua L

Malaria is a devastating parasitic disease with high levels of morbidity and mortality worldwide. Artemisinin, the active substance against malaria, is a sesquiterpenoid produced by Artemisia annua. To improve artemisinin content in the native A. annua plants, considerable efforts have been attempted, with genetic transformation serving as an effective strategy. Although, the most frequently-used cauliflower mosaic virus (CaMV) 35S (CaMV35S) promoter has proved to be efficient in A. annua transgenic studies, it appears to show weak activity in peltate glandular secretory trichomes (GSTs) of A. annua plants. Here, we characterized the 1727 bp fragment upstream from the translation start codon (ATG) of AaActin1, however, found it was inactive in tobacco. After removal of the 5' intron, the truncated AaActin1 promoter (tpACT) showed 69% and 50% activity of CaMV35S promoter in transiently transformed tobacco and stably transformed A. annua, respectively. β-glucuronidase (GUS) staining analysis showed that the tpACT promoter was capable of directing the constant expression of a foreign gene in peltate GSTs of transgenic A. annua, representing higher activity than CaMV35S promoter. Collectively, our study provided a novel promoter available for metabolic engineering of artemisinin biosynthesis in A. annua.

Exploring Aspartic Protease Inhibitor Binding to Design-Selective Antimalarials

Selectivity is a major issue in the development of drugs targeting pathogen aspartic proteases. Here, we explore the selectivity-determining factors by studying specifically designed malaria aspartic protease (plasmepsin) open-flap inhibitors. Metadynamics simulations are used to uncover the complex binding/unbinding pathways of these inhibitors and describe the critical transition states in atomistic resolution. The simulation results are compared with experimentally determined enzymatic activities. Our findings demonstrate that plasmepsin inhibitor selectivity can be achieved by targeting the flap loop with hydrophobic substituents that enable ligand binding under the flap loop, as such a behavior is not observed for several other aspartic proteases. The ability to estimate the selectivity of compounds before they are synthesized is of considerable importance in drug design; therefore, we expect that our approach will be useful in selective inhibitor designs against not only aspartic proteases but also other enzyme classes.