Regulatory networks,genes and glycerophospholipid biosynthesis pathway in schistosomiasis: A systems biology view for pharmacological intervention

Understanding network topology through embracing the global dynamical regulation of genes in an active state space rather than traditional one-gene–one trait approach facilitates the rational drug development process. Schistosomiasis, a neglected tropical disease, has glycerophospholipids as abundant molecules present on its surface. Lack of effective clinical solutions to treat pathogens encourages us to carry out systems-level studies that could contribute to the development of an effective therapy. Development of a strategy for identifying drug targets by combined genome-scale metabolic network and essentiality analyses through in silico approaches provides tantalizing opportunity to investigate the role of protein/substrate metabolism. A genome-scale metabolic network model reconstruction represents choline–phosphate cytidyltransferase as the rate limiting enzyme and regulates the rate of phosphatidylcholine (PC) biosynthesis. The uptake of choline was regulated by choline concentration, promoting the regulation of phosphocholine synthesis. In Schistosoma, the change in developmental stage could result from the availability of choline, hampering its developmental cycle. There are no structural reports for this protein. In order to inhibit the activity of choline–phosphate cytidyltransferase (CCT), it was modeled by homology modeling using 1COZ as the template from Bacillus subtilis. The transition-state stabilization and catalytic residues were mapped as ‘HXGH’ and ‘RTEGISTT’ motif. CCT catalyzes the formation of CDP-choline from phosphocholine in which nucleotidyltransferase adds CTP to phosphocholine. The presence of phosphocholine permits the parasite to survive in an immunologically hostile environment. This feature endeavors development of an inhibitor specific for cytidyltransferase in Schistosoma. Flavonolignans were used to inhibit this activity in which hydnowightin showed the highest affinity as compared to miltefosine.     

Impact of the green tea ingredient epigallocatechin gallate and a short pentapeptide (Ile-Ile-Ala-Glu-Lys) on the structural organization of mixed micelles and the related uptake of cholesterol

Background

High levels of blood cholesterol are conventionally linked to an increased risk of developing cardiovascular disease (Grundy, 1986). Here we examine the molecular mode of action of natural products with known cholesterol-lowering activity, such as for example the green tea ingredient epigallocatechin gallate and a short pentapeptide, Ile-Ile-Ala-Glu-Lys.

A selective autophagy pathway takes an unconventional route

Selective autophagy transports specific cytoplasmic materials into lysosomes/vacuoles. In the case of macroautophagy the selectivity is mediated by receptors, which usually link the cargos to the machinery that sequesters them into the forming autophagosome. In our recent work, we found that fission yeast Nbr1, a homolog of the mammalian macroautophagy receptor NBR1, acts together with an unconventional autophagy-associated cargo sequestration apparatus, the endosomal sorting complexes required for transport (ESCRTs), to deliver 2 hydrolytic enzymes from the cytosol to the vacuole lumen. In this pathway, which we term the Nbr1-mediated vacuolar targeting (NVT) pathway, soluble cargos transit through the multi-vesicular body (MVB), rather than the autophagosome, on their way to the vacuole. Our findings reveal a novel mode of action of macroautophagy receptors and broaden our understanding of ESCRT-mediated autophagy.     

Immunophysical properties and prediction of activities for vaccinia virus complement control protein and smallpox inhibitor of complement enzymes using molecular dynamics and electrostatics

We present immunophysical modeling for VCP, SPICE, and three mutants using MD simulations and Poisson-Boltzmann-type electrostatic calculations. VCP and SPICE are homologous viral proteins that control the complement system by imitating, structurally and functionally, natural regulators of complement activation. VCP and SPICE consist of four CCP modules connected with short flexible loops. MD simulations demonstrate that the rather complex modules of VCP/SPICE and their mutants exhibit a high degree of intermodular spatial mobility, which is affected by surface mutations. Electrostatic calculations using snapshots from the MD trajectories demonstrate variable spatial distribution of the electrostatic potentials, which suggests dynamic binding properties. We use covariance analysis to identify correlated modular oscillations. We also use electrostatic similarity indices to cluster proteins with common electrostatic properties. Our results are compared with experimental data to form correlations between the overall positive electrostatic potential of VCP/SPICE with binding and activity. We show how these correlations can be used to predict binding and activity properties. This work is expected to be useful for understanding the function of native CCP-containing regulators of complement activation and receptors and for the design of antiviral therapeutics and complement inhibitors.     

Unraveling microalgal molecular interactions using evolutionary and structural bioinformatics

Microalgae are unicellular microorganisms indispensible for environmental stability and life on earth, because they produce approximately half of the atmospheric oxygen, with simultaneously feeding on the harmful greenhouse gas carbon dioxide. Using gene fusion analysis, a series of five fusion/fission events was identified, that provided the basis for critical insights to their evolutionary history. Moreover, the three-dimensional structures of both the fused and the component proteins were predicted, allowing us to envisage putative protein–protein interactions that are invaluable for the efficient usage, handling and exploitation of microalgae. Collectively, our proposed approach on the five fusion/fission alga protein events contributes towards the expansion of the microalgae knowledgebase, bridging protein evolution of the ancient microalgal species and the rapidly evolving, modern, bioinformatics field.     

Stabilization of the catalytic thiolate in a mammalian glutaredoxin: structure, dynamics and electrostatics of reduced pig glutaredoxin and its mutants

The variety of functions performed by proteins of the thioredoxin superfamily, including glutaredoxins, involves the wide range of redox potential associated with the -Cys-X-X-Cys- motif found in their active sites. The determinants of these differences in redox potential are still obscure. A better understanding requires a detailed characterization of the reduced state of these enzymes, especially because the lowered pK(a) of the reduced N-terminal active-site cysteine is a key feature of these enzymes' chemistry, including their redox potential. Analysis of the factors controlling this pK(a) is complicated by the apparent structural heterogeneity of the reduced active sites across glutaredoxins. In this family, pig glutaredoxin (pGrx) was one of the first to be functionally characterized, including some intriguing mutagenesis data, but a structure of its reduced state has been lacking. We used long molecular dynamics simulations and electrostatic calculations to analyze the structure, dynamics and electrostatics of reduced pGrx and some of its mutants. Comparison with experimental data is drawn whenever possible. It is shown that a dynamic model is essential to capture the structural properties of the cationic side-chains around the -Cys22-Pro23-Phe24-Cys25- sequence in the pGrx active site. Examples include Arg26, which can swing to stack on this sequence, and Lys19 which can contact the thiolate. However, contrary to a commonly held hypothesis, these cationic side-chains provide little stabilization for the thiolate, implying that they affect the enzymatic activity via other mechanisms. The pK(a) value of nucleophilic cysteine 22 (pK(a)(22)) is dominated by local hydrogen-bonds, formed only in a well-defined active-site conformation, supported by a comparison between the calculated and experimental values of pK(a)(22). The edge of the aromatic ring of Phe24 is polar enough to contribute to stabilize the thiolate, consistent with the conserved aromatic side-chain at this position in the glutaredoxin motif. The locality and directionality of the hydrogen bonds in the active site suffice to explain the vast difference between the pK(a) values of its two cysteine residues. A control of the cysteine pK(a) values by local hydrogen bonds implies that the peripheral ionized side-chains can evolve independently of the maintenance of these pK(a) values, maybe guided instead by substrate recognition. Comparison with other glutaredoxins indicates that the calculated pK(a) values of the N-terminal cysteine are better conserved than those of the C-terminal cysteine. Overall, a methodological strategy to systematically compare all reduced enzymes of this family emerges.     

Evaluation of blood-brain barrier penetration and examination of binding to human serum albumin of 7-O-arylpiperazinylcoumarins as potential antipsychotic agents

The delivery of drugs to the brain is complicated by the multiple factors including low blood–brain barrier (BBB) passive permeability, active BBB efflux systems, and plasma protein binding. Thus, a detailed understanding of the transport of the new potent substances through the membranes is vitally important and their physico-chemical characteristics should be analyzed at first. This work presents an evaluation of drug likeness of eight 7-O-arylpiperazinylcoumarin derivatives with high affinity towards serotoninergic receptors 5-HT_1A and 5-HT_2A with particular analysis of the requirements for the CNS chemotherapeutics. The binding constants to human serum albumin (HSA) were determined at physiological pH using fluorescence spectroscopy, and then their mode of action was explained by analysis of theoretical HSA complexes. Dynamic simulation of systems allowed for reliable evaluation of the interaction strength. The analyzed coumarins were able to pass BBB, and they present good drug likeness properties. They showed high affinities to HSA (log K_Q = 5.3–6.0 which corresponds to −8.12 to −7.15 kcalmol⁻¹ of Gibbs free energy). The changes of the emission intensity upon binding to HSA were scrutinized showing the different mode of action for 4-phenylpiperazinylcoumarins. The values of computed Gibbs free energy and determined on the basis of experimentally obtained binding constants log K_Q coincide suggesting a good quality of the theoretical model. Overall the 8-acetyl-7-O-arylpiperazinyl-4-methylcoumarin derivatives represent valuable lead compounds to be further tested in various preclinical assays as a possible chemotherapeutics against CNS diseases. Studied coumarins can be metabolized by cytochrome P450 to aldehydes and hydroxy derivatives. The existence of other binding sites inside HSA than Sudlow’s site 1 was postulated. The longer aliphatic linker between coumarin and piperazine moieties favored binding to HSA in other than Sudlow site 1 pocket.