Building and deploying a cyberinfrastructure for the data-driven design of chemical systems and the exploration of chemical space

Abstract

The use of modern data science has recently emerged as a promising new path to tackling the complex challenges involved in the creation of next-generation chemistry and materials. However, despite the appeal of this potentially transformative development, the chemistry community has yet to incorporate it as a central tool in every-day work. Our research program is designed to enable and advance this emerging research approach. It is centred around the creation of a software ecosystem that brings together physics-based modelling, high-throughput in silico screening and data analytics (i.e. the use of machine learning and informatics for the validation, mining and modelling of chemical data). This cyberinfrastructure is devised to offer a comprehensive set of data science techniques and tools as well as a general-purpose scope to make it as versatile and widely applicable as possible. It also emphasises user-friendliness to make it accessible to the community at large. It thus provides the means for the large-scale exploration of chemical space and for a better understanding of the hidden mechanisms that determine the properties of complex chemical systems. Such insights can dramatically accelerate, streamline and ultimately transform the way chemical research is conducted. Aside from serving as a production-level tool, our cyberinfrastructure is also designed to facilitate and assess methodological innovation. Both the software and method development work are driven by concrete molecular design problems, which also allow us to assess the efficacy of the overall cyberinfrastructure.     

Induced fit docking,free energy calculation and molecular dynamics studies on Mycobacterium tuberculosis alanine racemase inhibitor

Alar, a Pyridoxal 5′-phosphate (PLP)-dependent bacterial enzyme is responsible for the racemisation of L-alanine into D-alanine which is essential for the peptidoglycan biosynthesis in both Gram-positive and Gram-negative bacteria. In the present study, we performed induced fit docking, binding free energy calculation and molecular dynamics simulation to elucidate the Alar inhibition potential of 1,2,4-thiadiazolidine-3,5-dione-based inhibitor 1. The inhibitor binds to the hydrophobic groove of Alar and the binding was found to be stable throughout 20-ns MD simulation. Induced fit docking result showed that Lys42, Tyr46, Tyr175 and Tyr364 residues are primarily responsible for the stabilisation of inhibitor–protein complex. Further, high negative van der Waals binding free energy value of –38.88 kcal/mol, indicated it as the main driving force for the inhibitor binding. Based on the information obtained from this study, we designed few molecules as potent Alar inhibitor. In order to gain structural insight and to validate the stability of complex, we performed 20-ns MD simulation of the designed molecule D1. Results obtained from this study can be used for the design of M. tuberculosis Alar potent inhibitors lacking affinity for the co-factor PLP.     

Potential roles of YCF54 and ferredoxin‐NADPH reductase for magnesium protoporphyrin monomethylester cyclase

Chlorophyll is synthesized from activated glutamate in the tetrapyrrole biosynthesis pathway through at least 20 different enzymatic reactions. Among these, the MgProto monomethylester (MgProtoME) cyclase catalyzes the formation of a fifth isocyclic ring to tetrapyrroles to form protochlorophyllide. The enzyme consists of two proteins. The CHL27 protein is proposed to be the catalytic component, while LCAA/YCF54 likely acts as a scaffolding factor. In comparison to other reactions of chlorophyll biosynthesis, this enzymatic step lacks clear elucidation and it is hardly understood, how electrons are delivered for the NADPH‐dependent cyclization reaction. The present study intends to elucidate more precisely the role of LCAA/YCF54. Transgenic Arabidopsis lines with inactivated and overexpressed YCF54 reveal the mutual stability of YCF54 and CHL27. Among the YCF54‐interacting proteins, the plastidal ferredoxin‐NADPH reductase (FNR) was identified. We showed in N. tabacum and A. thaliana that a deficit of FNR1 or YCF54 caused MgProtoME accumulation, the substrate of the cyclase, and destabilization of the cyclase subunits. It is proposed that FNR serves as a potential donor for electrons required in the cyclase reaction and connects chlorophyll synthesis with photosynthetic activity.     

<Emphasis Type="Italic">CBS</Emphasis> mutations and <Emphasis Type="Italic">MTFHR</Emphasis> SNPs causative of hyperhomocysteinemia in Pakistani children

Three index patients with hyperhomocysteinemia and ocular anomalies were screened for cystathionine beta synthase (CBS) and methylenetetrahydrofolate reductase (MTHFR) polymorphisms. Genotyping of hyperhomocysteinemia associated MTHFR polymorphisms C677T (rs1801133) and A1298C (rs1801131) was done by PCR-restriction fragment length polymorphism. Sanger sequencing was performed for CBS exonic sequences along with consensus splice sites. In the case of MTHFR polymorphisms, all the patients were heterozygous CT for the single nucleotide polymorphism (SNP) C677T and were therefore carriers of the risk allele (T), while the patients were homozygous CC for the risk genotype of the SNP A1298C. CBS sequencing resulted in the identification of two novel mutations, a missense change (c.467T>C; p.Leu156Pro) in exon 7 and an in-frame deletion (c.808_810del; p.Glu270del) in exon 10. In addition, a recurrent missense mutation (c.770C>T; p.Thr257Met) in exon 10 of the gene was also identified. The mutations were present homozygously in the patients and were inherited from the carrier parents. This is the first report from Pakistan where novel as well as recurrent CBS mutations causing hyperhomocysteinemia and lens dislocation in three patients from different families are being reported with the predicted effect of the risk allele of the MTHFR SNP in causing hyperhomocysteinemia.     

Regulation of aging and oxidative stress pathways in aged pancreatic islets using alpha-lipoic acid

Molecular and Cellular Biochemistry - Oxidative stress has been involved in the aging process and the pathogenesis of type-2 diabetes, which is a serious health problem worldwide. This study...     

Discovery of potent polyphosphate kinase 1 (PPK1) inhibitors using structure‐based exploration of PPK1Pharmacophoric space coupled with docking analyses

Inorganic polyphosphate (polyP) is present in all living forms of life. Studied mainly in prokaryotes, polyP and its associated enzymes are vital in diverse metabolic activities, in some structural functions, and most importantly in stress responses. Bacterial species, including many pathogens, encode a homolog of a major polyP synthesis enzyme, Poly Phosphate Kinase (PPK) with 2 different genes coding for PPK1 and PPK2. Genetic deletion of the ppk1 gene leads to reduced polyP levels and the consequent loss of virulence and stress adaptation responses. This far, no PPK1 homolog has been identified in higher‐order eukaryotes, and, therefore, PPK1 represents a novel target for chemotherapy. The aim of the current study is to investigate PPK1 from Escherichia coli with comprehensive understanding of the enzyme's structure and binding sites, which were used to design pharmacophores and screen a library of compounds for potential discovery of selective PPK1 inhibitors. Verification of the resultant inhibitors activities was conducted using a combination of mutagenic and chemical biological approaches. The metabolic phenotypic maps of the wild type E. coli (WT) and ppk1 knockout mutant were generated and compared with the metabolic map of the chemically inhibited WT. In addition, biofilm formation ability was measured in WT, ppk1 knockout mutant, and the chemically inhibited WT. The results demonstrated that chemical inhibition of PPK1, with the designed inhibitors, was equivalent to gene deletion in altering specific metabolic pathways, changing the metabolic fingerprint, and suppressing the ability of E. coli to form a biofilm.