Manali Joshi

Small molecule modulators of eukaryotic initiation factor 2α kinases, the key regulators of protein synthesis.

Eukaryotic initiation factor 2 alpha kinases (eIF-2α kinases) are key mediators of stress response in cells. In mammalian cells, there are four eIF-2α kinases, namely HRI (Heme-Regulated Inhibitor), PKR (RNA-dependent Protein Kinase), PERK (PKR-like ER Kinase) and GCN2 (General Control Non-derepressible 2). These kinases get activated during diverse cytoplasmic stress conditions and phosphorylate the alpha-subunit of eIF2, leading to global protein synthesis inhibition. Therefore, eIF-2α kinases play a vital role in various cellular processes such as proliferation, differentiation, apoptosis and cell signaling. Deregulation of eIF-2α kinases and protein synthesis has been linked to numerous pathological conditions such as certain cancers, anemia and neurodegenerative disorders. Thus, modulation of these kinases by small molecules holds a great therapeutic promise. In this review we have compiled the available information on inhibitors and activators of these four eIF-2α kinases. The review concludes with a note on the selectivity issue of currently available modulators and future perspectives for the design of specific small molecule probes.

Molecular insights into the dynamics of pharmacogenetically important N-terminal variants of the human β2-adrenergic receptor.

The human β2-adrenergic receptor (β2AR), a member of the G-protein coupled receptor (GPCR) family, is expressed in bronchial smooth muscle cells. Upon activation by agonists, β2AR causes bronchodilation and relief in asthma patients. The N-terminal polymorphism of β2AR at the 16th position, Arg16Gly, has warranted a lot of attention since it is linked to variations in response to albuterol (agonist) treatment. Although the β2AR is one of the well-studied GPCRs, the N-terminus which harbors this mutation, is absent in all available experimental structures. The goal of this work was to study the molecular level differences between the N-terminal variants using structural modeling and atomistic molecular dynamics simulations. Our simulations reveal that the N-terminal region of the Arg variant shows greater dynamics than the Gly variant, leading to differential placement. Further, the position and dynamics of the N-terminal region, further, affects the ligand binding-site accessibility. Interestingly, long-range effects are also seen at the ligand binding site, which is marginally larger in the Gly as compared to the Arg variant resulting in the preferential docking of albuterol to the Gly variant. This study thus reveals key differences between the variants providing a molecular framework towards understanding the variable drug response in asthma patients.

Pharmacophore-Based Virtual Screening to Aid in the Identification of Unknown Protein Function

A pharmacophore‐based virtual screening method was developed and validated for use in predicting the function of a novel protein in terms of small metabolite binding. Five test cases were used for the validation study which spanned two different folds, four superfamilies, and three enzyme classes. Binding sites were predicted using a combination of two methods (CASTp and THEMATICS). The binding site was mapped with chemical probes representing hydrogen‐bond donor, acceptor, negative ionizable, positive ionizable, and hydrophobe. The interaction maps were converted to three or four feature pharmacophore models and used to search a database containing 80 018 tautomers/protomers/conformers of 10 535 metabolites. The pharmacophore‐based virtual screening eliminated >92% of the database as potential substrates and retrieved specific hits, which were ranked using a physics‐based scoring function. The known substrate or product was ranked within the top 0.7% and substrate‐like compounds within the top 1% of the metabolite database for all of the five test cases. The results suggest that using this pharmacophore‐based virtual screening is a time‐efficient strategy that can be applied to screen large databases to help predict the function of small metabolite binding proteins.

Ensemble-Based Virtual Screening and Experimental Validation of Inhibitors Targeting a Novel Site of Human DNMT1

Human DNA methyltransferase1 (hDNMT1) is responsible for preserving DNA methylation patterns that play important regulatory roles in differentiation and development. Misregulation of DNA methylation has thus been linked to many syndromes, life style diseases, and cancers. Developing specific inhibitors of hDNMT1 is an important challenge in the area since the currently targeted cofactor and substrate binding site share structural features with various proteins. In this work, we generated a structural model of the active form of hDNMT1 and identified that the 5‐methylcytosine (5‐mC) binding site of the hDNMT1 is structurally unique to the protein. This site has been previously demonstrated to be critical for methylation activity. We further performed multiple nanosecond time scale atomistic molecular dynamics simulations of the structural model followed by virtual screening of the Asinex database to identify inhibitors targeting the 5‐mC site. Two compounds were discovered that inhibited hDNMT1 in vitro, one of which also showed inhibition in vivo corroborating the screening procedure. This study thus identifies and attempts to validate for the first time a unique site of hDNMT1 that could be harnessed for rationally designing highly selective and potent hypomethylating agents.

Structural insights and functional implications of inter-individual variability in β2-adrenergic receptor.

The human β2-adrenergic receptor (β2AR) belongs to the G protein-coupled receptor (GPCR) family and due to its central role in bronchodilation, is an important drug target. The inter-individual variability in β2AR has been implicated in disease susceptibility and differential drug response. In this work, we identified nine potentially deleterious non-synonymous single nucleotide polymorphisms (nsSNPs) using a consensus approach. The deleterious nsSNPs were found to cluster near the ligand binding site and towards the G-protein binding site. To assess their molecular level effects, we built structural models of these receptors and performed atomistic molecular dynamics simulations. Most notably, in the Phe290Ser variant we observed the rotameric flip of Trp2866.48, a putative activation switch that has not been reported in β2AR thus far. In contrast, the variant Met82Lys was found to be the most detrimental to epinephrine binding. Additionally, a few of the nsSNPs were seen to cause perturbations to the lipid bilayer, while a few lead to differences at the G-protein coupling site. We are thus able to classify the variants as ranging from activating to damaging, prioritising them for experimental studies.

A high affinity phage-displayed peptide as a recognition probe for the detection of Salmonella Typhimurium

Salmonellosis is one of the most common and widely distributed foodborne diseases. A sensitive and robust detection method of Salmonella Typhimurium (S. Typhimurium) in food can critically prevent a disease outbreak. In this work, the use of phage displayed peptides was explored for the detection of S. Typhimurium. A phage-displayed random dodecapeptide library was subjected to biopanning against lipopolysaccharide (LPS) of S. Typhimurium. The peptide NFMESLPRLGMH (pep49) derived from biopanning displayed a high affinity (25.8nM) for the LPS of S. Typhimurium and low cross-reactivity with other strains of Salmonella and related Gram-negative bacteria. Molecular insights into the interaction of pep49 with the LPS of S. Typhimurium was gleaned using atomistic molecular dynamics simulations and docking. It was deduced that the specificity of pep49 with S. Typhimurium LPS originated from the interactions of pep49 with abequose that is found only in the O-antigen of S. Typhimurium. Further, pep49 was able to detect S. Typhimurium at a LOD of 10(3) CFU/mL using ELISA, and may be a potential cost efficient alternative to antibodies.

Unravelling the structural interactions between PKR kinase domain and its small molecule inhibitors using computational approaches

The RNA-dependent protein kinase (PKR), an eIF2α kinase plays an important role in anti-viral response, apoptosis and cell survival. It is also implicated to play a role in several cancers, metabolic and neurodegenerative disorders. A few ATP competitive inhibitors of the PKR have been reported in the literature with promising results in vitro and in vivo. The aim of this study was to unravel the structural interactions between these inhibitors and the PKR kinase domain using molecular simulations and docking. Our study reveals that the reported inhibitors bind in the adenine pocket and form hydrogen bonds with the hinge region and vdW interactions with non-polar residues in the binding site. The most potent inhibitor has several favorable interactions with the binding site and induces the P-loop to fold inward, creating a significant hydrophobic enclosure for itself. The computed binding free energies of these inhibitors are in accord with experimental data (IC50). Strategies to design potent and selective PKR inhibitors are discussed to overcome the reported promiscuity.

Computational insights into the interaction of small molecule inhibitors with HRI kinase domain

The Heme-Regulated Inhibitor (HRI) kinase regulates globin synthesis in a heme-dependent manner in reticulocytes and erythroid cells in bone marrow. Inhibitors of HRI have been proposed to lead to an increased amount of haemoglobin, benefitting anaemia patients. A series of indeno[1,2-c]pyrazoles were discovered to be the first known in vitro inhibitors of HRI. However, the structural mechanism of inhibition is yet to be understood. The aim of this study was to unravel the binding mechanism of these inhibitors using molecular dynamic simulations and docking. The docking scores were observed to correlate well with experimentally determined pIC50 values. The inhibitors were observed to bind in the ATP-binding site forming hydrogen bonds with the hinge region and van der Waals interactions with non-polar residues in the binding site. Further, quantitative structure–activity relationship (QSAR) studies were performed to correlate the structural features of the inhibitors with their biological activity. The developed QSAR models were found to be statistically significant in terms of internal and external predictabilities. The presence of chlorine atoms and the hydroxymethyl groups were found to correlate with higher activity. The identified binding modes and the descriptors can support future rational identification of more potent and selective small molecule inhibitors for this kinase which are of therapeutic importance in the context of various human pathological disorders.

Characterizing clinically relevant natural variants of GPCRs using computational approaches

G protein-coupled receptors (GPCRs) are an important class of drug targets owing to their physiological role. A large number of clinically relevant single nucleotide polymorphisms (SNPs) have been observed in GPCRs that are linked to disease susceptibility and adverse drug response. It is therefore important to characterize the variants in order to improve GPCR therapeutics. Here, we discuss computational methods coupling molecular dynamics simulations with docking and free energy calculations to characterize the functional differences in GPCR variants. The hallmark of this approach is the explicit incorporation of receptor and membrane dynamics that allows us to analyze short- and long-range effects in the variant receptors. We use the SNPs reported in β2-adrenergic receptor (β2AR) as a test case and highlight the recent successes in analyzing structural and dynamic differences in a series of population variants. The computational approach we discuss here has a twofold benefit: it helps to unravel the molecular mechanisms underlying hypo- or hyperfunctionality of variant receptors as well as prioritizing novel variants that must be experimentally tested.

Membrane-induced organization and dynamics of the N-terminal domain of chemokine receptor CXCR1: insights from atomistic simulations

The CXC chemokine receptor 1 (CXCR1) is an important member of the G protein-coupled receptor (GPCR) family in which the extracellular N-terminal domain has been implicated in ligand binding and selectivity. The structure of this domain has not yet been elucidated due to its inherent dynamics, but experimental evidence points toward membrane-dependent organization and dynamics. To gain molecular insight into the interaction of the N-terminal domain with the membrane bilayer, we performed a series of microsecond time scale atomistic simulations of the N-terminal domain of CXCR1 in the presence and absence of POPC bilayers. Our results show that the peptide displays a high propensity to adopt a β-sheet conformation in the presence of the membrane bilayer. The interaction of the peptide with the membrane bilayer was found to be transient in our simulations. Interestingly, a scrambled peptide, containing the same residues in a randomly varying sequence, did not exhibit membrane-modulated structural dynamics. These results suggest that sequence-dependent electrostatics, modulated by the membrane, could play an important role in folding of the N-terminal domain. We believe that our results reinforce the emerging paradigm that cellular membranes could be important modulators of function of G protein-coupled receptors such as CXCR1.