Stereochemistry-Dependent Labeling of Organelles with a Near-Infrared-Emissive Phosphorus-Bridged Rhodamine Dye in Live-Cell Imaging.

The development of near-infrared (NIR) fluorophores that have both excellent chemical stability and photostability, as well as efficient cell permeability, is highly demanded. In this study, we present phospha-rhodamine (POR) dyes which display significantly improved performance for protein labeling. This is achieved by incorporating a 2-carboxy-3-benzothiophenyl group at the 9-position of the xanthene scaffold. The resulting cis and trans isomers were successfully isolated and structurally characterized using X-ray diffraction. The HaloTag ligand conjugates of the two isomers exhibited different staining abilities in live cells. While the cis isomer showed non-specific accumulation on the organelle membranes, the trans isomer selectively labeled the HaloTag-fused proteins, enabling the long-term imaging of cell division and the 5-color imaging of cell organelles. Molecular dynamics simulations of the HaloTag ligand conjugates within the lipid membrane suggested that the cis isomer is more prone to forming oligomers in the membrane. In contrast, the oligomerization of the trans isomer is effectively suppressed by its interaction with the lipid molecules. By taking advantage of the superior labeling performance of the trans isomer and its NIR-emissive properties, multi-color time-lapse super-resolution 3D imaging, namely super-resolution 5D-imaging, of the interconnected network between the endoplasmic reticulum and microtubules was achieved in living cells.

CRY2 isoform selectivity of a circadian clock modulator with antiglioblastoma efficacy.

The mammalian cryptochrome isoforms, CRY1 and CRY2, are core circadian clock regulators that work redundantly. Recent studies revealed distinct roles of these closely related homologs in clock output pathways. Isoform-selective control of CRY1 and CRY2 is critical for further understanding their redundant and distinct roles. KL001 was the first identified small-molecule CRY modulator that activates both CRY1 and CRY2. SHP656 is an orally available KL001 derivative and has shown efficacy in blood glucose control and inhibition of glioblastoma stem cell (GSC) growth in animal models. However, CRY isoform selectivity of SHP656 was uncharacterized, limiting understanding of the roles of CRY1 and CRY2. Here, we report the elucidation of CRY2 selectivity of SHP656. SHP656 lengthened cellular circadian period in a CRY2-dependent manner and selectively interacted with CRY2. By determining the X-ray crystal structure of CRY2 in complex with SHP656 and performing molecular dynamics simulations, we elucidated compound interaction mechanisms. SHP656 binding was compatible with the intrinsic CRY2 gatekeeper W417 "in" orientation and also a close "further in" conformation. Perturbation of W417 interaction with the lid loop resulted in a reduced effect of SHP656 on CRY2, supporting an important role of gatekeeper orientation in isoform selectivity. We also identified the R form of SHP656 (called SHP1703) as the active isomer. Treatment with SHP1703 effectively reduced GSC viability. Our results suggest a direct role of CRY2 in glioblastoma antitumorigenesis and provide a rationale for the selective modulation of CRY isoforms in the therapeutic treatment of glioblastoma and other circadian clock-related diseases.

Molecular insights into substrate binding mechanism of undecaprenyl pyrophosphate with membrane integrated phosphatidyl glycerophosphate phosphatase B (PgpB) using molecular dynamics simulation approach.

Undecaprenyl phosphate (C55-P) acts as carrier lipid in the synthesis of peptidoglycan, which is de novo synthesized from dephosphorylation of undecaprenyl pyrophosphate (C55-PP). The phosphatidylglycerol phosphate phosphatase B (PgpB) catalyzes the dephosphorylation of C55-PP and forms C55-P. As no structural study has been made regarding the binding of C55-PP to PgpB, in the current study, in silico molecular docking, followed by 150 ns molecular dynamics simulation of the putative binding complex in membrane/solvent environment has been performed to understand conformational dynamics. Results are compared with simulated apo form and PE inhibitor-bound form. Analysis of correlated residual fluctuation network in apo form, C55-PP bound and PE inhibitor-bound form suggests that difference in dynamic coupling between TM domain and α2 and α3 helix of periplasmic domain provides ligand binding to facilitate catalysis or to show inhibitory activity. Distance distribution in catalytic residual pair, H207-R104; H207-R201 and H207-D211 which stabilizes phosphate-enzyme intermediate shows a narrow peak in 2.4-3.6 Å in substrate-bound compared to apo form. Binding interactions and binding free energy analyses complement the partial inhibition of PE where PE has less binding free energy compared to the C55-PP substrate as well as the difference in binding interaction with catalytic pocket. Thus, the present study provides how substrate binding couples the movement in TM domain and periplasmic domain which might help in the understanding of active site communication in PgpB. C55-PP phosphatase interactions with a catalytic pocket of PgpB provide new insight for designing drugs against bacterial infection.

Conformational Dynamics of thiM Riboswitch To Understand the Gene Regulation Mechanism Using Markov State Modeling and the Residual Fluctuation Network Approach.

Thiamine pyrophosphate (TPP) riboswitch is a cis-regulatory element in the noncoding region of mRNA. The aptamer domain of TPP riboswitch detects the high abundance of coenzyme thiamine pyrophosphate (TPP) and modulates the gene expression for thiamine synthetic gene. The mechanistic understanding in recognition of TPP in aptamer domain and ligand-induced compactness for folding of expression platform are most important to designing novel modulators. To understand the dynamic behavior of TPP riboswitch upon TPP binding, molecular dynamics simulations were performed for 400 ns in both apo and TPP bound forms of thiM riboswitch from E. coli and analyzed in terms of eRMSD-based Markov state modeling and residual fluctuation network. Markov state models show good correlations in transition probability among metastable states from simulated trajectory and generated models. Structural compactness in TPP bound form is observed which is correlated with SAXS experiment. The importance of junction of P4 and P5 is evident during dynamics, which correlates with FRET analysis. The dynamic nature of two sensor forearms is due to the flexible P1 helix, which is its intrinsic property. The transient state in TPP-bound form was observed in the Markov state model, along with stable states. We believe that this transient state is responsible to assist the influx and outflux of ligand molecule by creating a solvent channel around the junction region of P4 and P5 and such a structure was anticipated in FRET analysis. The dynamic nature of riboswitch is dependent on the interaction between residues on distal loops L3 and L5/P3 and junction P4 and P5, J3/2 which stabilize the J2/4. It helps in the transfer of allosteric information between J2/4 and P3/L5 tertiary docking region through the active site residues. Understanding such information flow will benefit in highlighting crucial residues in highly dynamic and kinetic systems. Here, we report the residues and segments in riboswitch that play vital roles in providing stability and this can be exploited in designing inhibitors to regulate the functioning of riboswitches.

Drug-Target Interactions: Prediction Methods and Applications.

Identifying the interactions between drugs and target proteins is a key step in drug discovery. This not only aids to understand the disease mechanism, but also helps to identify unexpected therapeutic activity or adverse side effects of drugs. Hence, drug-target interaction prediction becomes an essential tool in the field of drug repurposing. The availability of heterogeneous biological data on known drug-target interactions enabled many researchers to develop various computational methods to decipher unknown drug-target interactions. This review provides an overview on these computational methods for predicting drug-target interactions along with available webservers and databases for drug-target interactions. Further, the applicability of drug-target interactions in various diseases for identifying lead compounds has been outlined.

Identification of novel Nicotinamide Phosphoribosyltransferase (NAMPT) inhibitors using computational approaches.

Nicotinamide Phosphoribosyltransferase (NAMPT) is a rate-limiting enzyme in the biosynthesis of NAD. Cancer cells have elevated poly [ADP-Ribose] polymerase 1 (PARP) activity as well as the immense necessity of ATP: thereby consuming NAD at a higher rate than normal tissues. The perturbation of these intracellular processes is more sensitive and highly dependent on NAMPT to maintain the required NAD levels. Functional inhibition of NAMPT is, therefore, a promising drug target in therapeutic oncology. In this study, the importance of intermolecular contacts was realized based on contact occupancy and favorable energetic from molecular dynamic simulation to discern non-critical contacts of four different classes of potential NAMPT inhibitor bound complexes. Further, pharmacophore modeling, molecular docking, a quantum mechanical properties and MD simulation, as well as active site residual network communication were employed to identify potential leads. Present studies identified two leads, 2 and 3 which have better binding free energy compared to known inhibitors and showed stable hydrogen bonding and hydrophobic contacts with ß barrel cavity lining residues in the active site of the dimer interface (A'B). Lead 2 containing fluorene as central core and lead 3 having phenyl-benzamide as a core showed stable moiety which was observed from electronic property analysis. Active site residual communication in identified leads bound complex also showed similarity to known inhibitor complexes. Compounds containing these moieties were not reported until now against NAMPT inhibition and can be considered as novel cores for future development of drugs to inhibit NAMPT function.

Molecular docking and molecular dynamics simulations of fumarate hydratase and its mutant H235N complexed with pyromellitic acid and citrate.

Fumarase catalyzes the reversible, stereospecific hydration/dehydration of fumarate to L-malate during the Kreb's cycle. In the crystal structure of the tetrameric fumarase, it was found that some of the active site residues S145, T147, N188 G364 and H235 had water-mediated hydrogen bonding interactions with pyromellitic acid and citrate which help to the protonation state for the conversion of fumarate to malate. When His 235 is mutated with Asn (H235N), water-mediated interactions were lost due to the shifting of active site water molecule by 0.7 Å away. Molecular dynamics (MD) simulations were also carried out by NAMD and analyzed using Assisted Model Building with Energy Refinement (AMBER) program to better understand the conformational stability and other aspects during the binding of pyromellitic acid and citrate with native and mutant FH. The role of hydrogen bonds and hydrophobic interactions was also analyzed. The present study confirms that the H235N mutation has a major effect on the catalytic activity of fumarase which is evident from the biochemical studies.

Identification of new BACE1 inhibitors using Pharmacophore and Molecular dynamics simulations approach.

Inhibition of ß-Secretase (BACE1) is crucial for the treatment of Alzheimer's disease (AD). Availability of BACE1 crystal structures in both apo and complexed forms enables to find structure-based BACE1 inhibitors for controlling AD. There are two catalytic aspartates (ASP32 and ASP228) presents in the active domain of BACE1. In order to understand the binding mechanism and structure-activity relationship of amidine-containing BACE1 inhibitors, molecular docking, and pharmacophore and 3D-QSAR studies have been carried out with 34 amidine derivatives to develop a pharmacophore model. Pharmacophore-based virtual screening (PBVS) has been performed against BACE1 (PDB ID: 2FDP), using three chemical databases (CoCoCo, Enamine, Zinc), which yielded 6000 hit compounds. These compounds were further analyzed using structure-based docking in hierarchical filtering approaches of Glide such as HTVS, SP, and XP precision modes. The docking results show that binding orientations of the inhibitors at Asp dyad active site amino acid residues of ß-Secretase. Results from glide XP docking and induced fit docking showed that four leads (Lead1, Lead3, Lead4 and Lead5) have good interactions with the target protein in comparison with cocrystal (amino-ethylene BACE1 inhibitor). Further, molecular dynamics (MD) simulation for these leads bound with BACE1 shows conformational stability and difference in dynamical flap behaviors of the active site with cocrystal inhibitor. Binding free energetic using MM-GB/SA approaches suggest lead 1 and lead 3 has comparably favorable binding to cocrystal inhibitor. Thus, the present study emphasizes these leads for an effective drug to treat Alzheimer disease.

Comparison of proteomic profiles of the venoms of two of the 'Big Four' snakes of India, the Indian cobra (Naja naja) and the common krait (Bungarus caeruleus), and analyses of their toxins.

Snake venoms are mixtures of biologically-active proteins and peptides, and several studies have described the characteristics of some of these toxins. However, complete proteomic profiling of the venoms of many snake species has not yet been done. The Indian cobra (Naja naja) and common krait (Bungarus caeruleus) are elapid snake species that are among the 'Big Four' responsible for the majority of human snake envenomation cases in India. As understanding the composition and complexity of venoms is necessary for successful treatment of envenomation in humans, we utilized three different proteomic profiling approaches to characterize these venoms: i) one-dimensional SDS-PAGE coupled with in-gel tryptic digestion and electrospray tandem mass spectrometry (ESI-LC-MS/MS) of individual protein bands; ii) in-solution tryptic digestion of crude venoms coupled with ESI-LC-MS/MS; and iii) separation by gel-filtration chromatography coupled with tryptic digestion and ESI-LC-MS/MS of separated fractions. From the generated data, 81 and 46 different proteins were identified from N. naja and B. caeruleus venoms, respectively, belonging to fifteen different protein families. Venoms from both species were found to contain a variety of phospholipases A2 and three-finger toxins, whereas relatively higher numbers of snake venom metalloproteinases were found in N. naja compared to B. caeruleus venom. The analyses also identified less represented venom proteins including L-amino acid oxidases, cysteine-rich secretory proteins, 5'-nucleotidases and venom nerve growth factors. Further, Kunitz-type serine protease inhibitors, cobra venom factors, phosphodiesterases, vespryns and aminopeptidases were identified in the N. naja venom, while acetylcholinesterases and hyaluronidases were found in the B. caeruleus venom. We further analyzed protein coverage (Lys/Arg rich and poor regions as well as potential glycosylation sites) using in-house software. These studies expand our understanding of the proteomes of the venoms of these two medically-important species.

Identification of novel natural inhibitor for NorM - a multidrug and toxic compound extrusion transporter - an insilico molecular modeling and simulation studies.

The emergence of bacterial multidrug resistance is an increasing problem in treatment of infectious diseases. An important cause for the multidrug resistance of bacteria is the expression of multidrug efflux transporters. The multidrug and toxic compound extrusion (MATE) transporters are most recently recognized as unique efflux system for extrusion of antimicrobials and therapeutic drugs due to energy stored in either Na+ or H+ electrochemical gradient. In the present study, high throughput virtual screening of natural compound collections against NorM - a MATE transporter from Neisseria gonorrhea (NorM-NG) has been carried out followed by flexible docking. The molecular simulation in membrane environment has been performed for understanding the stability and binding energetic of top lead compounds. Results identified a compound from the Indian medicinal plant "Terminalia chebula" which has good binding free energy compared to substrates (rhodamine 6 g, ethidium) and more favorable interactions with the central cavity forming active site residues. The compound has restricted movement in TM7, TM8, and TM1, thus blocking the disruption of Na+ - coordination along with equilibrium state bias towards occlude state of NorM transporter. Thus, this compound blocks the effluxing pathway of antimicrobial drugs and provides as a natural bioactive lead inhibitor against NorM transporter in drug-resistant gonorrhea.

Determination on the binding of thiadiazole derivative to human serum albumin: a spectroscopy and computational approach.

4-[3-acetyl-5-(acetylamino)-2,3-dihydro-1,3,4-thiadiazole-2-yl]phenyl benzoate from the family of thiadiazole derivative has been newly synthesized. It has good anticancer activity as well as antibacterial and less toxic in nature, its binding characteristics are therefore of huge interest for understanding pharmacokinetic mechanism of the drug. The binding of thiadiazole derivative to human serum albumin (HSA) has been investigated by studying its quenching mechanism, binding kinetics and the molecular distance, r between the donor (HSA) and acceptor (thiadiazole derivative) was estimated according to Forster's theory of non-radiative energy transfer. The Gibbs free energy (ΔG), enthalpy (ΔH) and entropy (ΔS) changes of temperature-dependent Kb was calculated, which explains that the reaction is spontaneous and exothermic. The microenvironment of HSA have also been studied using synchronous fluorescence spectroscopy, and the feature of thiadiazole derivative-induced structural changes of HSA have been carried using Fourier transform infrared spectroscopy and the Molecular modelling simulations explore the hydrophobic and hydrogen bonding interactions.

Fasciola gigantica thioredoxin glutathione reductase: Biochemical properties and structural modeling.

Platyhelminth thioredoxin glutathione reductase (TGR) is a multifunctional enzyme that crosstalk between the conventional thioredoxin (Trx) and glutathione (GSH) system. It has been validated as a potential drug target in blood flukes. In the present study, we have performed a biochemical study on Fasciola gigantica TGR with substrates DTNB and GSSG. The Michaelis constant (Km) with DTNB was found to be 4.34±0.12µM while it was 61.15±1.50µM with GSSG. The kinetic results were compared with the TGR activities of other helminths. FgTGR showed typical hysteretic behavior with GSSG as other TGRs. We also described a homology-based structure of FgTGR. The cofactors (NADPH and FAD) and substrates (GSSG and DTNB) were docked, and two possible binding sites for substrates were identified in a single chain. The substrates were found to bind more favorably in the second site of TrxR domains. We also presented the first report on binding interaction of DTNB with a TGR. DTNB forms H-bond with His204 and Arg450 of chain A, Sec597, and Gly598 from chain B, salt-bridge with Lys124, and numerous other hydrophobic interactions. Helminth TGR represents an important enzyme in the redox and antioxidant system; hence, its inhibition can be used as an effective strategy against liver flukes.

Insights into the binding of thiosemicarbazone derivatives with human serum albumin: spectroscopy and molecular modelling studies.

4-[(1Z)-1-(2-carbamothioylhydrazinylidene)ethyl]phenyl acetate [Ace semi],4-[(1Z)-1-(2-carbamothioylhydrazinylidene)ethyl]phenyl propanoate [Pro semi] from the family of thiosemicarbazones derivative has been newly synthesized. It has good anticancer activity as well as antibacterial and it is also less toxic in nature, its binding characteristics are therefore of huge interest for understanding pharmacokinetic mechanism of the drug. The binding of thiosemicarbazone derivative to human serum albumin (HSA) has been investigated by studying its quenching mechanism, binding kinetics and the molecular distance (r) between donor (HSA) and acceptor (thiosemicarbazone derivative) was estimated according to Forster's theory of non-radiative energy transfer using fluorescence spectroscopy. The binding dynamics has been elaborated using synchronous fluorescence spectroscopy, and the feature of thiosemicarbazone derivative induced structural changes of HSA has been studied by circular dichorism, Fourier transform infrared spectroscopy. Molecular modelling simulations explore the hydrophobic interaction and hydrogen bonding which stabilizes the interaction.

Synthesis and molecular modelling studies of novel sulphonamide derivatives as dengue virus 2 protease inhibitors.

Development of antivirals for dengue is now based on rational approach targeting the enzymes involved in its life cycle. Among the targets available for inhibition of dengue virus, non-structural protein NS2B-NS3 protease is considered as a promising target for the development of anti-dengue agents. In the current study we have synthesized a series of 4-(1,3-dioxo-2,3-dihydro-1H-isoindol-2-yl)benzene-1-sulphonamide derivatives and screened for DENV2 protease activity. Compounds 16 and 19 showed IC50 of DENV2 Protease activity with 48.2 and 121.9µM respectively. Molecular docking and molecular dynamic simulation studies were carried out to know the binding mode responsible for the activity. MD simulations revealed that, NS2B/NS3 protease was more stable when it binds with the active compound. Structure optimization of the lead compounds 16 and 19 and their co-crystallization studies are underway.