A spatially resolved elemental nanodomain organization within acidocalcisomes in Trypanosoma cruzi.

How are ions distributed in the three-dimensional (3D) volume confined in a nanoscale compartment? Regulation of ionic flow in the intracellular milieu has been explained by different theoretical models and experimentally demonstrated for several compartments with microscale dimensions. Most of these models predict a homogeneous distribution of ions seconds or milliseconds after an initial diffusion step formed at the ion translocation site, leaving open questions when it comes to ion/element distribution in spaces/compartments with nanoscale dimensions. Due to the influence of compartment size on the regulation of ionic flow, theoretical variations of classical models have been proposed, suggesting heterogeneous distributions of ions/elements within nanoscale compartments. Nonetheless, such assumptions have not been fully proven for the 3D volume of an organelle. In this work, we used a combination of cutting-edge electron microscopy techniques to map the 3D distribution of diffusible elements within the whole volume of acidocalcisomes in trypanosomes. Cryofixed cells were analyzed by scanning transmission electron microscopy tomography combined with elemental mapping using a high-performance setup of X-ray detectors. Results showed the existence of elemental nanodomains within the acidocalcisomes, where cationic elements display a self-excluding pattern. These were validated by Pearson correlation analysis and in silico molecular dynamic simulations. Formation of element domains within the 3D space of an organelle is demonstrated. Distribution patterns that support the electrodiffusion theory proposed for nanophysiology models have been found. The experimental pipeline shown here can be applied to a variety of models where ion mobilization plays a crucial role in physiological processes.

Insights into nucleoside hydrolase from Leishmania donovani inhibition: A new bioaffinity chromatography-based screening assay and docking studies.

Leishmaniasis, a group of neglected infectious diseases, encompasses a serious health concern, particularly with visceral leishmaniasis exhibiting potentially fatal outcomes. Nucleoside hydrolase (NH) has a fundamental role in the purine salvage pathway, crucial for Leishmania donovani survival, and presents a promising target for developing new drugs for visceral leishmaniasis treatment. In this study, LdNH was immobilized into fused silica capillaries, resulting in immobilized enzyme reactors (IMERs). The LdNH-IMER activity was monitored on-flow in a multidimensional liquid chromatography system, with the IMER in the first dimension. A C18 analytical column in the second dimension furnished the rapid separation of the substrate (inosine) and product (hypoxanthine), enabling direct enzyme activity monitoring through product quantification. LdNH-IMER exhibited high stability and was characterized by determining the Michaelis-Menten constant. A known inhibitor (1-(ß-d-Ribofuranosyl)-4-quinolone derivative) was used as a model to validate the established method in inhibitor recognition. Screening of three additional derivatives of 1-(ß-d-Ribofuranosyl)-4-quinolone led to the discovery of novel inhibitors, with compound 2a exhibiting superior inhibitory activity (Ki = 23.37 ± 3.64 µmol/L) compared to the employed model inhibitor. Docking and Molecular Dynamics studies provided crucial insights into inhibitor interactions at the enzyme active site, offering valuable information for developing new LdNH inhibitors. Therefore, this study presents a novel screening assay and contributes to the development of potent LdNH inhibitors.

Enzymatic and structural properties of human glutamine:fructose-6-phosphate amidotransferase 2 (hGFAT2).

Glycoconjugates play a central role in several cellular processes, and alteration in their composition is associated with numerous human pathologies. Substrates for cellular glycosylation are synthesized in the hexosamine biosynthetic pathway, which is controlled by the glutamine:fructose-6-phosphate amidotransfera-se (GFAT). Human isoform 2 GFAT (hGFAT2) has been implicated in diabetes and cancer; however, there is no information about structural and enzymatic properties of this enzyme. Here, we report a successful expression and purification of a catalytically active recombinant hGFAT2 (rhGFAT2) in Escherichia coli cells fused or not to a HisTag at the C-terminal end. Our enzyme kinetics data suggest that hGFAT2 does not follow the expected ordered bi-bi mechanism, and performs the glucosamine-6-phosphate synthesis much more slowly than previously reported for other GFATs. In addition, hGFAT2 is able to isomerize fructose-6-phosphate into glucose-6-phosphate even in the presence of equimolar amounts of glutamine, which results in unproductive glutamine hydrolysis. Structural analysis of a three-dimensional model of rhGFAT2, corroborated by circular dichroism data, indicated the presence of a partially structured loop in the glutaminase domain, whose sequence is present in eukaryotic enzymes but absent in the E. coli homolog. Molecular dynamics simulations suggest that this loop is the most flexible portion of the protein and plays a key role on conformational states of hGFAT2. Thus, our study provides the first comprehensive set of data on the structure, kinetics, and mechanics of hGFAT2, which will certainly contribute to further studies on the (patho)physiology of hGFAT2.

A novel receptor for platelet-activating factor and lysophosphatidylcholine in Trypanosoma cruzi.

The lipid mediators, platelet-activating factor (PAF) and lysophosphatidylcholine (LPC), play relevant pathophysiological roles in Trypanosoma cruzi infection. Several species of LPC, including C18:1 LPC, which mimics the effects of PAF, are synthesized by T. cruzi. The present study identified a receptor in T. cruzi, which was predicted to bind to PAF, and found it to be homologous to members of the progestin and adiponectin family of receptors (PAQRs). We constructed a three-dimensional model of the T. cruzi PAQR (TcPAQR) and performed molecular docking to predict the interactions of the TcPAQR model with C16:0 PAF and C18:1 LPC. We knocked out T. cruzi PAQR (TcPAQR) gene and confirmed the identity of the expressed protein through immunoblotting and immunofluorescence assays using an anti-human PAQR antibody. Wild-type and knockout (KO) parasites were also used to investigate the in vitro cell differentiation and interactions with peritoneal mouse macrophages; TcPAQR KO parasites were unable to react to C16:0 PAF or C18:1 LPC. Our data are highly suggestive that PAF and LPC act through TcPAQR in T. cruzi, triggering its cellular differentiation and ability to infect macrophages.

Staphylococcus aureus Exfoliative Toxin E, Oligomeric State and Flip of P186: Implications for Its Action Mechanism.

Staphylococcal exfoliative toxins (ETs) are glutamyl endopeptidases that specifically cleave the Glu381-Gly382 bond in the ectodomains of desmoglein 1 (Dsg1) via complex action mechanisms. To date, four ETs have been identified in different Staphylococcus aureus strains and ETE is the most recently characterized. The unusual properties of ETs have been attributed to a unique structural feature, i.e., the 180° flip of the carbonyl oxygen (O) of the nonconserved residue 192/186 (ETA/ETE numbering), not conducive to the oxyanion hole formation. We report the crystal structure of ETE determined at 1.61 Å resolution, in which P186(O) adopts two conformations displaying a 180° rotation. This finding, together with free energy calculations, supports the existence of a dynamic transition between the conformations under the tested conditions. Moreover, enzymatic assays showed no significant differences in the esterolytic efficiency of ETE and ETE/P186G, a mutant predicted to possess a functional oxyanion hole, thus downplaying the influence of the flip on the activity. Finally, we observed the formation of ETE homodimers in solution and the predicted homodimeric structure revealed the participation of a characteristic nonconserved loop in the interface and the partial occlusion of the protein active site, suggesting that monomerization is required for enzymatic activity.

In silico identification of noncompetitive inhibitors targeting an uncharacterized allosteric site of falcipain-2.

Falcipain-2 (FP-2) is a Plasmodium falciparum hemoglobinase widely targeted in the search for antimalarials. FP-2 can be allosterically modulated by various noncompetitive inhibitors that have been serendipitously identified. Moreover, the crystal structures of two inhibitors bound to an allosteric site, termed site 6, of the homolog enzyme human cathepsin K (hCatK) suggest that the equivalent region in FP-2 might play a similar role. Here, we conduct the rational identification of FP-2 inhibitors through virtual screenings (VS) of compounds into several pocket-like conformations of site 6, sampled during molecular dynamics (MD) simulations of the free enzyme. Two noncompetitive inhibitors, ZINC03225317 and ZINC72290660, were confirmed using in vitro enzymatic assays and their poses into site 6 led to calculated binding free energies matching the experimental ones. Our results provide strong evidence about the allosteric inhibition of FP-2 through binding of small molecules to site 6, thus opening the way toward the discovery of new inhibitors against this enzyme.

Improving Blind Docking in DOCK6 through an Automated Preliminary Fragment Probing Strategy.

Probing protein surfaces to accurately predict the binding site and conformation of a small molecule is a challenge currently addressed through mainly two different approaches: blind docking and cavity detection-guided docking. Although cavity detection-guided blind docking has yielded high success rates, it is less practical when a large number of molecules must be screened against many detected binding sites. On the other hand, blind docking allows for simultaneous search of the whole protein surface, which however entails the loss of accuracy and speed. To bridge this gap, in this study, we developed and tested BLinDPyPr, an automated pipeline which uses FTMap and DOCK6 to perform a hybrid blind docking strategy. Through our algorithm, FTMap docked probe clusters are converted into DOCK6 spheres for determining binding regions. Because these spheres are solely derived from FTMap probes, their locations are contained in and specific to multiple potential binding pockets, which become the regions that are simultaneously probed and chosen by the search algorithm based on the properties of each candidate ligand. This method yields pose prediction results (45.2-54.3% success rates) comparable to those of site-specific docking with the classic DOCK6 workflow (49.7-54.3%) and is half as time-consuming as the conventional blind docking method with DOCK6.

Water Bridges Play a Key Role in Affinity and Selectivity for Malarial Protease Falcipain-2.

Falcipain-2 (FP-2) is hemoglobinase considered an attractive drug target of Plasmodium falciparum. Recently, it has been shown that peptidomimetic nitriles containing a 3-pyridyl (3Pyr) moiety at P2 display high affinity and selectivity for FP-2 with respect to human cysteine cathepsins (hCats), outperforming other P2-Pyr isomers and analogs. Further characterization demonstrated that certain P3 variants of these compounds possess micromolar inhibition of parasite growth in vitro and no cytotoxicity against human cell lines. However, the structural determinants underlying the selectivity of the 3Pyr-containing nitriles for FP-2 remain unknown. In this work, we conduct a thorough computational study combining MD simulations and free energy calculations to decipher the bases of the selectivity of the aforementioned nitriles. Our results reveal that water bridges involving the nitrogen and one carboxyl oxygen of I85 and D234 of FP-2, respectively, and the nitrogen of the neutral 3Pyr moiety, which are either less prevalent or nonexistent in the other complexes, explain the experimental activity profiles. The presence of crystallographic waters close to the bridging water positions in the studied proteases strongly supports the occurrence of such interactions. Overall, our findings suggest that selective FP-2 inhibitors can be designed by promoting water bridge formation at the bottom of the S2 subsite and/or by introducing complementary groups that displace the bridging water.

Pisum sativum Defensin 1 Eradicates Mouse Metastatic Lung Nodules from B16F10 Melanoma Cells.

Psd1 is a pea plant defensin which can be actively expressed in Pichia pastoris and shows broad antifungal activity. This activity is dependent on fungal membrane glucosylceramide (GlcCer), which is also important for its internalization, nuclear localization, and endoreduplication. Certain cancer cells present a lipid metabolism imbalance resulting in the overexpression of GlcCer in their membrane. In this work, in vitroassays using B16F10 cells showed that labeled fluorescein isothiocyanate FITC-Psd1 internalized into live cultured cells and targeted the nucleus, which underwent fragmentation, exhibiting approximately 60% of cells in the sub-G0/G1 stage. This phenomenon was dependent on GlcCer, and the participation of cyclin-F was suggested. In a murine lung metastatic melanoma model, intravenous injection of Psd1 together with B16F10 cells drastically reduced the number of nodules at concentrations above 0.5 mg/kg. Additionally, the administration of 1 mg/kg Psd1 decreased the number of lung inflammatory cells to near zero without weight loss, unlike animals that received melanoma cells only. It is worth noting that 1 mg/kg Psd1 alone did not provoke inflammation in lung tissue or weight or vital signal losses over 21 days, inferring no whole animal cytotoxicity. These results suggest that Psd1 could be a promising prototype for human lung anti-metastatic melanoma therapy.

Investigation of the binding mode of a novel cruzain inhibitor by docking, molecular dynamics, ab initio and MM/PBSA calculations.

Chagas disease remains a major health problem in South America, and throughout the world. The two drugs clinically available for its treatment have limited efficacy and cause serious adverse effects. Cruzain is an established therapeutic target of Trypanosoma cruzi, the protozoan that causes Chagas disease. Our group recently identified a competitive cruzain inhibitor (compound 1) with an IC50 = 15 µM that is also more synthetically accessible than the previously reported lead, compound 2. Prior studies, however, did not propose a binding mode for compound 1, hindering understanding of the structure-activity relationship and optimization. Here, the cruzain binding mode of compound 1 was investigated using docking, molecular dynamics (MD) simulations with ab initio derived parameters, ab initio calculations, and MM/PBSA. Two ligand protonation states and four binding poses were evaluated. A careful ligand parameterization method was employed to derive more physically meaningful parameters than those obtained by automated tools. The poses of unprotonated 1 were unstable in MD, showing large conformational changes and diffusing away from the binding site, whereas the protonated form showed higher stability and interaction with negatively charged residues Asp161 and Cys25. MM/PBSA also suggested that these two residues contribute favorably to binding of compound 1. By combining results from MD, ab initio calculations, and MM/PBSA, a binding mode of 1 is proposed. The results also provide insights for further optimization of 1, an interesting lead compound for the development of new cruzain inhibitors.

Predicting binding modes of reversible peptide-based inhibitors of falcipain-2 consistent with structure-activity relationships.

Falcipain-2 (FP-2) is a major hemoglobinase of Plasmodium falciparum, considered an important drug target for the development of antimalarials. A previous study reported a novel series of 20 reversible peptide-based inhibitors of FP-2. However, the lack of tridimensional structures of the complexes hinders further optimization strategies to enhance the inhibitory activity of the compounds. Here we report the prediction of the binding modes of the aforementioned inhibitors to FP-2. A computational approach combining previous knowledge on the determinants of binding to the enzyme, docking, and postdocking refinement steps, is employed. The latter steps comprise molecular dynamics simulations and free energy calculations. Remarkably, this approach leads to the identification of near-native ligand conformations when applied to a validation set of protein-ligand structures. Overall, we proposed substrate-like binding modes of the studied compounds fulfilling the structural requirements for FP-2 binding and yielding free energy values that correlated well with the experimental data. Proteins 2017; 85:1666-1683. © 2017 Wiley Periodicals, Inc.

Insights into cytochrome bc1 complex binding mode of antimalarial 2-hydroxy-1,4-naphthoquinones through molecular modelling.

BACKGROUND: Malaria persists as a major public health problem. Atovaquone is a drug that inhibits the respiratory chain of Plasmodium falciparum, but with serious limitations like known resistance, low bioavailability and high plasma protein binding. OBJECTIVES: The aim of this work was to perform molecular modelling studies of 2-hydroxy-1,4-naphthoquinones analogues of atovaquone on the Qo site of P. falciparum cytochrome bc1 complex (Pfbc1) to suggest structural modifications that could improve their antimalarial activity. METHODS: We have built the homology model of the cytochrome b (CYB) and Rieske iron-sulfur protein (ISP) subunits from Pfbc1 and performed the molecular docking of 41 2-hydroxy-1,4-naphthoquinones with known in vitro antimalarial activity and predicted to act on this target. FINDINGS: Results suggest that large hydrophobic R2 substituents may be important for filling the deep hydrophobic Qo site pocket. Moreover, our analysis indicates that the H-donor 2-hydroxyl group may not be crucial for efficient binding and inhibition of Pfbc1 by these atovaquone analogues. The C1 carbonyl group (H-acceptor) is more frequently involved in the important hydrogen bonding interaction with His152 of the Rieske ISP subunit. MAIN CONCLUSIONS: Additional interactions involving residues such as Ile258 and residues required for efficient catalysis (e.g., Glu261) could be explored in drug design to avoid development of drug resistance by the parasite.

Unraveling HIV protease flaps dynamics by Constant pH Molecular Dynamics simulations.

The active site of HIV protease (HIV-PR) is covered by two flaps. These flaps are known to be essential for the catalytic activity of the HIV-PR, but their exact conformations at the different stages of the enzymatic pathway remain subject to debate. Understanding the correct functional dynamics of the flaps might aid the development of new HIV-PR inhibitors. It is known that, the HIV-PR catalytic efficiency is pH-dependent, likely due to the influence of processes such as charge transfer and protonation/deprotonation of ionizable residues. Several Molecular Dynamics (MD) simulations have reported information about the HIV-PR flaps. However, in MD simulations the protonation of a residue is fixed and thus it is not possible to study the correlation between conformation and protonation state. To address this shortcoming, this work attempts to capture, through Constant pH Molecular Dynamics (CpHMD), the conformations of the apo, substrate-bound and inhibitor-bound HIV-PR, which differ drastically in their flap arrangements. The results show that the HIV-PR flaps conformations are defined by the protonation of the catalytic residues Asp25/Asp25' and that these residues are sensitive to pH changes. This study suggests that the catalytic aspartates can modulate the opening of the active site and substrate binding.

Allosteric regulation of the Plasmodium falciparum cysteine protease falcipain-2 by heme.

During the erythrocytic cycle of Plasmodium falciparum malaria parasites break down host hemoglobin, resulting in the release of free heme (ferriprotoporphyrin IX). Heme is a generator of free radicals that cause oxidative stress, but it is detoxified by crystallization into hemozoin inside the food vacuole. We evaluated the interaction of heme and heme analogues with falcipain-2, a P. falciparum food vacuole cysteine protease that plays a key role in hemoglobin digestion. Heme bound to falcipain-2 with a 1:1 stoichiometry, and heme inhibited falcipain-2 activity against both human hemoglobin and chromogenic peptide substrates through a noncompetitive-like mechanism. A series of porphyrin analogues was screened for inhibition of falcipain-2, demonstrating a minor contribution of iron to heme-falcipain-2 interaction, and revealing dependence on both propionic and vinyl groups for inhibition of falcipain-2 by heme. Docking and molecular dynamics simulation unveiled a novel, inducible heme-binding moiety in falcipain-2 adjacent to the catalytic site. Kinetic data suggested that the noncompetitive-like inhibition was substrate inhibition induced by heme. Collectively these data suggest that binding of heme to falcipain-2 may limit the accumulation of free heme in the parasite food vacuole, providing a means of heme detoxification in addition to crystallization into hemozoin.

Conformational changes in human Hsp70 induced by high hydrostatic pressure produce oligomers with ATPase activity but without chaperone activity.

We investigated the folding of the 70 kDa human cytosolic inducible protein (Hsp70) in vitro using high hydrostatic pressure as a denaturing agent. We followed the structural changes in Hsp70 induced by high hydrostatic pressure using tryptophan fluorescence, molecular dynamics, circular dichroism, high-performance liquid chromatography gel filtration, dynamic light scattering, ATPase activity, and chaperone activity. Although monomeric, Hsp70 is very sensitive to hydrostatic pressure; after pressure had been removed, the protein did not return to its native sate but instead formed oligomeric species that lost chaperone activity but retained ATPase activity.

Impact of M36I polymorphism on the interaction of HIV-1 protease with its substrates: insights from molecular dynamics.

BACKGROUND: Over the last decades, a vast structural knowledge has been gathered on the HIV-1 protease (PR). Noticeably, most of the studies focused the B-subtype, which has the highest prevalence in developed countries. Accordingly, currently available anti-HIV drugs target this subtype, with considerable benefits for the corresponding patients. RESULTS: Herein, we used molecular dynamics simulations to investigate the role of this polymorphism on the interaction of PR with six of its natural cleavage-sites substrates. CONCLUSIONS: With multiple approaches and analyses we identified structural and dynamical determinants associated with the changes found in the binding affinity of the M36I variant. This mutation influences the flexibility of both PR and its complexed substrate. The observed impact of M36I, suggest that combination with other non-B subtype polymorphisms, could lead to major effects on the interaction with the 12 known cleavage sites, which should impact the virion maturation.

Structure-based discovery of novel cruzain inhibitors with distinct trypanocidal activity profiles.

Over 110 years after the first formal description of Chagas disease, the trypanocidal drugs thus far available have limited efficacy and several side effects. This encourages the search for novel treatments that inhibit T. cruzi targets. One of the most studied anti-T. cruzi targets is the cysteine protease cruzain; it is associated with metacyclogenesis, replication, and invasion of the host cells. We used computational techniques to identify novel molecular scaffolds that act as cruzain inhibitors. First, with a docking-based virtual screening, we identified compound 8, a competitive cruzain inhibitor with a Ki of 4.6 µM. Then, aided by molecular dynamics simulations, cheminformatics, and docking, we identified the analog compound 22 with a Ki of 27 µM. Surprisingly, despite sharing the same isoquinoline scaffold, compound 8 presented higher trypanocidal activity against the epimastigote forms, while compound 22, against the trypomastigotes and amastigotes. Taken together, compounds 8 and 22 represent a promising scaffold for further development of trypanocidal compounds as drug candidates for treating Chagas disease.

Three decades targeting falcipains to develop antiplasmodial agents: what have we learned and what can be done next?

Malaria is a devastating infectious disease that affects large swathes of human populations across the planet's tropical regions. It is caused by parasites of the genus Plasmodium, with Plasmodium falciparum being responsible for the most lethal form of the disease. During the intraerythrocytic stage in the human hosts, malaria parasites multiply and degrade hemoglobin (Hb) using a battery of proteases, which include two cysteine proteases, falcipains 2 and 3 (FP-2 and FP-3). Due to their role as major hemoglobinases, FP-2 and FP-3 have been targeted in studies aiming to discover new antimalarials and numerous inhibitors with activity against these enzymes, and parasites in culture have been identified. Nonetheless, cross-inhibition of human cysteine cathepsins remains a serious hurdle to overcome for these compounds to be used clinically. In this article, we have reviewed key functional and structural properties of FP-2/3 and described different compound series reported as inhibitors of these proteases during decades of active research in the field. Special attention is also paid to the wide range of computer-aided drug design (CADD) techniques successfully applied to discover new active compounds. Finally, we provide guidelines that, in our understanding, will help advance the rational discovery of new FP-2/3 inhibitors.

GSAFold: a new application of GSA to protein structure prediction.

The folding process defines three-dimensional protein structures from their amino acid chains. A protein's structure determines its activity and properties; thus knowing such conformation on an atomic level is essential for both basic and applied studies of protein function and dynamics. However, the acquisition of such structures by experimental methods is slow and expensive, and current computational methods mostly depend on previously known structures to determine new ones. Here we present a new software called GSAFold that applies the generalized simulated annealing (GSA) algorithm on ab initio protein structure prediction. The GSA is a stochastic search algorithm employed in energy minimization and used in global optimization problems, especially those that depend on long-range interactions, such as gravity models and conformation optimization of small molecules. This new implementation applies, for the first time in ab initio protein structure prediction, an analytical inverse for the Visitation function of GSA. It also employs the broadly used NAMD Molecular Dynamics package to carry out energy calculations, allowing the user to select different force fields and parameterizations. Moreover, the software also allows the execution of several simulations simultaneously. Applications that depend on protein structures include rational drug design and structure-based protein function prediction. Applying GSAFold in a test peptide, it was possible to predict the structure of mastoparan-X to a root mean square deviation of 3.00 Å.

Genetic diversity of the melanocortin-1 receptor in an admixed population of Rio de Janeiro: Structural and functional impacts of Cys35Tyr variant.

The melanocortin-1 receptor (MC1R) is one of the key proteins involved in the regulation of melanin production and several polymorphisms have been associated with different phenotypes of skin and hair color in human and nonhuman species. Most of the knowledge is centered on more homogeneous populations and studies involving an admixed group of people should be encouraged due to the great importance of understanding the human color variation. This work evaluates the MC1R diversity and the possible impacts of MC1R variants in an admixed sample population of Rio de Janeiro, Brazil, which is a product of Native American, African, and European miscegenation. Sequencing of complete coding region and part of the 3´UTR of MC1R gene identified 31 variants including one insertion and three novel synonymous substitutions in sample population grouped according to skin, hair and eye pigmentation levels. In nonmetric multidimensional scaling analysis (NMDS), three main clusters were identified, in which the Brazilian dark skin group remained in the African cluster whereas the intermediate and the light skin color phenotype in the European one. None gathered with Asians since their immigration to Brazil was a recent event. In silico analyses demonstrated that Cys35Tyr, Ile155Thr and Pro256Ser, found in our population, have a negative effect on receptor function probably due to changes on the receptor structure. Notably, Cys35Tyr mutation could potentially impair agonist binding. Altogether, this work contributes to the understanding of the genetic background of color variation on an admixed population and gives insights into the damaging effects of MC1R variants.