An insight to the dynamics of conserved water molecular triad in IMPDH II (human): recognition of cofactor and substrate to catalytic Arg 322

Dynamics of nucleoplasm in human leukemia cells: A thrust towards designing anti-leukemic agents

The human inosine monophosphate dehydrogenase (hIMPDH) is a metabolic enzyme that possesses a unique ability to self-assemble into higher-order structures, forming cytoophidia. The hIMPDH II isoform is more active in chronic myeloid leukemia (CML) cancer cells, making it a promising target for anti-leukemic therapy. However, the structural details and molecular mechanisms of the dynamics of hIMPDHcytoophidia assembly in vitro need to be better understood, and it is crucial to reconstitute the computational nucleoplasm model with cytophilic-like polymers in vitro to characterize their structure and function. Finally, a computational model and its dynamics of the nucleoplasm for CML cells have been proposed in this short review. This research on nucleoplasm aims to aid the scientific community's understanding of how metabolic enzymes like hIMPDH function in cancer and normal cells. However, validating and justifying the computational results from modeling and simulation with experimental data is essential. The new insights gained from this research could explain the structure/topology, geometrical, and electronic consequences of hIMPDH inhibitors on leukemic and normal cells. They could lead to further advancements in the knowledge of nucleoplasmic chemical reaction dynamics.

Structural and Dynamical Impact of Water Molecules at Substrate- or Product-Binding Sites in Human GMPR Enzyme: A Study by Molecular Dynamics Simulations

Human guanosine monophosphate reductase (hGMPR) enzyme maintains the intracellular balance between adenine and guanine nucleotide pools, and it is an excellent target for the design of isoform-specific antileukemic agents. In the present study, we have investigated solvation properties of substrate GMP or product inosine-5'-monophosphate (IMP)-binding pocket of hGMPR by employing molecular dynamics simulations on conformations A (substrate GMP), B [substrate GMP with cofactor nicotinamide adenine dinucleotide phosphate (NDP)], C (product IMP with cofactor NDP), and D (product IMP). Nineteen water sites are identified precisely; they are responsible for the catalytic activity of this site, control structural and dynamical integrity, and electronic consequences of GMP or IMP in the binding site of hGMPR. The water sites of category-1 (W1, W4, W5, W6, W13, and W15) in normal protein and category-2 (W2, W3, W7, W8, W10, W17, and W18) in cancerous protein are unique and stabilize the guanosine or inosine group of GMP or IMP for participation in the enzymatic reaction, whereas the remaining water centers either stabilize pentose sugar ribose or the phosphate group of GMP or IMP. Furthermore, water sites of category-4 (W11, W14, and W16) appear to be conserved in all conformations during the entire simulation. The GMP-binding site in cancerous protein 2C6Q is significantly expanded, and its dynamics are very different from normal protein 2BLE. Furthermore, unique interactions of GMP(N1)···W2···Asp129/Asn158, IMP(N1)···W3···Glu289, and IMP(O6)···W10···Ser270 might be used in a water mimic drug design for hGMPR-II. In this context, water finding probability, relative interaction energy (J) associated with water site W, entropy, and topologies of these three water sites are thermodynamically acceptable for the water displacement method by the modified ligand. Hence, their positions in the catalytic pocket may also facilitate future drug discovery for chronic myelogenous leukemia by the design of appropriately oriented chemical groups that may displace these water molecules to mimic their structural, electronic, and thermodynamic properties.

New biochemical insight of conserved water molecules at catalytic and structural Zn2+ ions in human matrix metalloproteinase-I: a study by MD-simulation

Human matrix metalloproteinase (MMP)-1 or collagenase-1 plays a significant role in embryonic development, tissue remodeling, and is also involved in several diseases like arthritis, metastasis, etc. Molecular dynamics simulation studies on hMMP-1 X-ray structures (PDB Id. 1CGE, 1CGF, 1CGL, 1HFC, and 2TCL) suggest that the three conserved water molecules (W_H/1, W_I, W_S) are coordinated with catalytic zinc (Zn_C), and one water molecule (W) is associated at structural zinc ion (Zn_S). Transition of the coordination geometry around Zn_C from tetrahedral to octahedral and tetrahedral to trigonal bipyramidal at Zn_S are also observed during the dynamics. Recognition of two zinc ions through water mediated bridges (Zn_C - W_H (W₁)…W₂….H₁₈₃ - Zn_S) and stabilization of secondary coordination zone around the metal ions indicates the possibility of Zn_C…Zn_S coupled catalytic mechanism in hMMP-I. This study not only reveals a functionally important role of conserved water molecules in hMMP-I but also highlights the involvement of other non catalytic residues, such as S172 and D170 in the catalytic mechanism. The results obtained in this study could be relevant for importance of conserved water mediated recognition site of the sequence residue id. 202(RWTNNFREY)210, interaction of W(tryptophan)203 to zinc bound histidine, their influence on the water molecules that are involved in bridging between Zn_C and Zn_S, and structure-based design of specific hMMP inhibitors. Graphical abstract Water mediated recognition of structural and catalytic zinc ions of hMMP-1 structure (MD simulatated conformation).

Elevated expression of IMPDH2 is associated with progression of kidney and bladder cancer

Novel molecular markers for cancer progression are valuable for the diagnosis and evaluation of treatment efficacies of the diseases. Expression of inosine 5'-monophosphate dehydrogenase type II (IMPDH2), a rate-limiting enzyme in the de novo guanine nucleotide biosynthesis, is up-regulated in various neoplasms, including prostate cancer and patient serum. However, whether IMPDH2 can serve as a biomarker for other urologic cancers is unknown. Paired patient tissue macroarrays were analyzed by immunohistochemistry, the IMPDH2 protein expression in these tissues was quantitated and expressed as immunoreactivity scores. Compared with non-cancerous tissues, IMPDH2 protein expression levels were significantly upregulated in kidney and bladder cancer, but no difference in testis cancer. In addition, expression of IMPDH2 was not associated with the disease clinical stages and pathological features. The findings suggest that overexpressed IMPDH2 can be used as a biomarker for kidney and bladder cancer diagnosis and is a potential therapeutic target for the diseases.

Trend of histone deacetylase inhibitors in cancer therapy: isoform selectivity or multitargeted strategy

Pharmacological inhibition of histone deacetylases (HDACs) has been successfully applied in the treatment of a wide range of disorders, including Parkinson's disease, infection, cardiac diseases, inflammation, and especially cancer. HDAC inhibitors (HDACIs) have been proved to be effective antitumor agents by various stages of investigation. At present, there are two opposite focuses of HDACI design in the cancer therapy, highly selective inhibitor strategy and dual- or multitargeted inhibitors. The former method, which is supposed to elucidate the function of individual HDAC and provide candidate inhibitors with fewer side effects, has been widely accepted by the inhibitor developer. The latter approach, though less practiced, has promising potential for the antitumor therapy based on HDACIs. Effective HDACIs, some of which are in clinic anticancer research, have been developed by both methods. In order to gain insight into HDACI design, the strategies and achievements of the two diverse methods are reviewed.

Conserved water-mediated recognition and dynamics of NAD+ (carboxamide group) to hIMPDH enzyme: water mimic approach toward the design of isoform-selective inhibitor

Inosine monophosphate dehydrogenase (IMPDH) enzyme involves in GMP biosynthesis pathway. Type I hIMPDH is expressed at lower levels in all cells, whereas type II is especially observed in acute myelogenous leukemia, chronic myelogenous leukemia cancer cells, and 10 ns simulation of the IMP-NAD(+) complex structures (PDB ID. 1B3O and 1JCN) have revealed the presence of a few conserved hydrophilic centers near carboxamide group of NAD(+). Three conserved water molecules (W1, W, and W1') in di-nucleotide binding pocket of enzyme have played a significant role in the recognition of carboxamide group (of NAD(+)) to D274 and H93 residues. Based on H-bonding interaction of conserved hydrophilic (water molecular) centers within IMP-NAD(+)-enzyme complexes and their recognition to NAD(+), some covalent modification at carboxamide group of di-nucleotide (NAD(+)) has been made by substituting the -CONH2group by -CONHNH2 (carboxyl hydrazide group) using water mimic inhibitor design protocol. The modeled structure of modified ligand may, though, be useful for the development of antileukemic agent or it could be act as better inhibitor for hIMPDH-II.

Role of salt bridge dynamics in inter domain recognition of human IMPDH isoforms: an insight to inhibitor topology for isoform-II

Inosine monophosphate dehydrogenase (IMPDH) enzyme involves in the biosynthesis pathway of guanosine nucleotide. Type II isoform of the enzyme is selectively upregulated in neoplastic fast replicating lymphocytes and CML cancer cells. The hIMPDH-II is an excellent target for antileukemic agent. The detailed investigation during MD-Simulation (15 ns) of three different unliganded structures (1B3O, 1JCN and 1JR1) have clearly explored the salt bridge mediated stabilization of inter or intra domain (catalytic domains I(N), I(C) with res. Id. 28-111 and 233-504, whereas two CBS domains C₁, C₂ are 112-171 and 172-232) in IMPDH enzyme which are mostly inaccessible in their X-rays structures. The salt bridge interaction in I(N)---C₁ inter-domain of hIMPDH-I, I(N)---C₂ of IMPDH-II and C₁---I(C) of nhIMPDH-II are discriminative features among the isoforms. The I(N)---C₂ recognition in hIMPDH-II (1B3O) is missing in type-I isoform (1JCN). The salt bridge interaction D232---K238 at the surface of protein and the involvement of three conserved water molecules or the hydrophilic centers (WA²³²(OD1), WB ²³²(OD2) and W²³⁸(NZ)) to those acidic and basic residues seem to be unique in hIMPDH-II. The hydrophilic susceptibility, geometrical and electronic consequences of this salt bridge interaction could be useful to design the topology of specific inhibitor for hIMPDH-II which may not be effective for hIMPDH-I. Possibly, the aliphatic ligand containing carboxyl, amide or hydrophilic groups with flexible structure may be implicated for hIMPDH-II inhibitor design using the conserved water mimic drug design protocol.

Biological evaluation and interaction of a newly designed anti-cancer Pd(II) complex and human serum albumin

The pharmacokinetics and pharmacodynamics of any drug will depend, largely, on the interaction that has with human serum albumin (HSA), the most abundant plasma protein. The interaction between newly synthesized Pd(II) complexe, 2,2'-bipyridin Butylglycinato Pd(II) nitrate, an anti-tumor component, with HSA was studied at different temperatures by fluorescence, far UV circular dichroism (CD), UV-visible spectrophotometry and theoretical approaches. The Pd(II) complex has a strong ability to quench the intrinsic fluorescence of HSA through a dynamic quenching procedure. The binding parameters and thermodynamic parameters, including δH°, δS° and δG° were calculated by fluorescence quenching method, indicated that hydrophobic forces play a major role in the interaction of Pd(II) complex with HSA. Based on Autodock, FRET (fluorescence resonance energy transfer) and fluorescence quenching data, it may be concluded that one of the binding sites in the complex of HSA is near the only one Trp of HSA (Trp214) in sub domain IIA of the protein. Far-UV-CD results indicated that Pd(II)-complex induced increase in the α-helical content of the protein. The anti-tumor property of the synthesized Pd(II) complex was studied by testing it on human tumor cell line K562. The 50% cytotoxic concentration (Cc₅₀) of complex was determined using MTT (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl tetrazolium bromide) assay. Also, fluorescence staining with DAPI (4,6-diamidino-2-phenylindole) revealed some typical nuclear changes that are characteristic of apoptosis which is induced at Cc₅₀ concentration of Pd(II) complex in K562 cell line after 24 h incubation. Our results suggest that Pd(II) complex is a promising anti-proliferative agent and should execute its biological effects by inducing apoptosis.

Molecular dynamics and docking studies on cardiac troponin C

Cardiac troponin C (cTnC) is the Ca²⁺ dependent switch for contraction in heart muscle making it a potential target for drug research in the therapy of heart failure. Calcium binding on Troponin C (TnC) triggers a series of conformational changes exposing a hydrophobic pocket in the N-domain of TnC (cNTnC), which leads to force generation. Mutations and acidic pH have been related to altering the sensitivity of TnC affecting the efficiency of the heart. Bepridil, identified as a calcium sensitizer to TnC, has been experimentally found to bind to the N-domain pocket of TnC but with negative cooperativity. Screening and de novo design were carried out using LUDI and AUTOLUDI programs in this work to identify and design potential ligands that can bind to the hydrophobic pocket of TnC. Two docking centers and multiple searching radii including 5 Å, 5.5 Å, 6 Å, 6.5 Å, 7.0 Å and 7.5 Å were used in LUDI to screen the ZINC database. Based on the LUDI docking results, 8 molecules were identified from the database with good potential to bind into the binding pocket and they were used as template molecules to generate a series of new molecules by AUTOLUDI design. Out of all the newly-designed molecules, 14 new ligands were recognized to be potential ligands that can bind and fit well into the binding pocket. These molecules can be used as starting molecules to develop TnC ligands. The binding stability and binding affinity of these molecules to the protein was further analyzed by molecular dynamics simulations. The results show that the binding energies, interactions and complex stabilities of 6 ligands are comparable to or better than bepridil.

Biphasic association of T7 RNA polymerase and a nucleotide analogue, cibacron blue as a model to understand the role of initiating nucleotide in the mechanism of enzyme action

T7 RNA polymerase (T7 RNAP) is an enzyme that utilizes ribonucleotides to synthesize the nascent RNA chain in a template-dependent manner. Here we have studied the interaction of T7 RNAP with cibacron blue, an anthraquinone monochlorotriazine dye, its effect on the function of the enzyme and the probable mode of binding of the dye. We have used difference absorption spectroscopy and isothermal titration calorimetry to show that the dye binds T7 RNAP in a biphasic manner. The first phase of the binding is characterized by inactivation of the enzyme. The second binding site overlaps with the common substrate-binding site of the enzyme. We have carried out docking experiment to map the binding site of the dye in the promoter bound protein. Competitive displacement of the dye from the high affinity site by labeled GTP and isothermal titration calorimetry of high affinity GTP bound enzyme with the dye suggests a strong correlation between the high affinity dye binding and the high affinity GTP binding in T7 RNAP reported earlier from our laboratory.

Using computer-aided drug design and medicinal chemistry strategies in the fight against diabetes

The aim of this work is to present a simple, practical and efficient protocol for drug design, in particular Diabetes, which includes selection of the illness, good choice of a target as well as a bioactive ligand and then usage of various computer aided drug design and medicinal chemistry tools to design novel potential drug candidates in different diseases. We have selected the validated target dipeptidyl peptidase IV (DPP-IV), whose inhibition contributes to reduce glucose levels in type 2 diabetes patients. The most active inhibitor with complex X-ray structure reported was initially extracted from the BindingDB database. By using molecular modification strategies widely used in medicinal chemistry, besides current state-of-the-art tools in drug design (including flexible docking, virtual screening, molecular interaction fields, molecular dynamics, ADME and toxicity predictions), we have proposed 4 novel potential DPP-IV inhibitors with drug properties for Diabetes control, which have been supported and validated by all the computational tools used herewith.

Thermal stability and unfolding pathways of Sso7d and its mutant F31A: insight from molecular dynamics simulation

The thermo-stability and unfolding behaviors of a small hyperthermophilic protein Sso7d as well as its single-point mutation F31A are studied by molecular dynamics simulation at temperatures of 300 K, 371 K and 500 K. Simulations at 300 K show that the F31A mutant displays a much larger flexibility than the wild type, which implies that the mutation obviously decreases the protein's stability. In the simulations at 371 K, although larger fluctuations were observed, both of these two maintain their stable conformations. High temperature simulations at 500 K suggest that the unfolding of these two proteins evolves along different pathways. For the wild-type protein, the C-terminal alpha-helix is melted at the early unfolding stage, whereas it is destroyed much later in the unfolding process of the F31A mutant. The results also show that the mutant unfolds much faster than its parent protein. The deeply buried aromatic cluster in the F31A mutant dissociates quickly relative to the wild-type protein at high temperature. Besides, it is found that the triple-stranded antiparallel β-sheet in the wild-type protein plays an important role in maintaining the stability of the entire structure.

Microsecond scale replica exchange molecular dynamic simulation of villin headpiece: an insight into the folding landscape

Reaching the experimental time scale of millisecond is a grand challenge for protein folding simulations. The development of advanced Molecular Dynamics techniques like Replica Exchange Molecular Dynamics (REMD) makes it possible to reach these experimental timescales. In this study, an attempt has been made to reach the multi microsecond simulation time scale by carrying out folding simulations on a three helix bundle protein, Villin, by combining REMD and Amber United Atom model. Twenty replicas having different temperatures ranging from 295 K to 390 K were simulated for 1.5 µs each. The lowest Root Mean Square Deviation (RMSD) structure of 2.5 Å was obtained with respect to native structure (PDB code 1VII), with all the helices formed. The folding population landscapes were built using segment-wise RMSD and Principal Components as reaction coordinates. These analyses suggest the two-stage folding for Villin. The combination of REMD and Amber United Atom model may be useful to understand the folding mechanism of various fast folding proteins.

The key residues of active sites on the catalytic fragment for paclitaxel interacting with poly (ADP-ribose) polymerase

Poly(ADP-ribose) polymerase (PARP) is regarded as a target protein for paclitaxel (PTX) to bind. An important issue is to identify the key residues as active sites for PTX interacting with PARP, which will help to understand the potential drug activity of PTX against cancer cells. Using docking method and MD simulation, we have constructed a refined structure of PTX docked on the catalytic function domain of PARP (PDB code: 1A26). The residues Glu327(988), Tyr246(907), Lys242(903), His165(826), Asp105(766), Gln102(763) and Gln98(759) in PARP are identified as potential sites involved in interaction with PTX according to binding energy (E(b)) between PTX and single residue calculated with B3LYP/6-31G(d,p). These residues form an active binding pocket located on the surface of the catalytic fragment, possibly interacting with the required groups of PTX leading to its activity against cancer cells. It is noted that most of the active sites make conatct with the "southern hemisphere" of PTX except for one residue, Tyr246(907), which interacts with the "northern hemisphere" of PTX. The conformation of PTX in complex with the catalytic fragment is observed as being T-shaped, similar to that complexed with β-tubulin. The total Eb of -269.9 kJ/mol represents the potent interaction between PTX and the catalytic fragment, implying that PTX can readily bind to the active pocket. The tight association of PTX with the catalytic fragment would inhibit PARP activation, suggesting a potential application of PTX as an effective antineoplastic agent.

Potent inhibitor design against H1N1 swine influenza: structure-based and molecular dynamics analysis for M2 inhibitors from traditional Chinese medicine database

The rapid spread of influenza virus subtype H1N1 poses a great threat to million lives worldwide. To search for new anti-influenza compounds, we performed molecular docking and molecular dynamics simulation to identify potential traditional Chinese medicine (TCM) constituents that could block influenza M2 channel activity. Quinic acid, genipin, syringic acid, cucurbitine, fagarine, and methyl isoferulate all have extremely well docking results as compared to control amantadine. Further de novo drug design suggests that derivatives of genipin and methyl isoferulate could have enhanced binding affinity towards M2 channel. Selected molecular dynamics simulations of M2-derivative complexes show stable hydrogen bond interactions between the derivatives and M2 residues, Ser10 and Ala9. To our best knowledge, this is the first study on the anti-viral activity of the above listed TCM compounds.

Conformational analysis of genotoxic benzo[a]pyrene-7,8-dione-duplex DNA adducts using a molecular dynamics method (II)

The conformations of the benzo[a]pyrene-7,8-quinone (BPQ) modified oligonucleotide were investigated using molecular dynamic simulation. In the initial structures, the central guanine base was modified with BPQ resulting in the formation of four structurally distinguishable 10-(N2-deoxyguanosyl)-9,10-dihydro-9-hydroxy benzo[a]pyrene-7,8-dione adducts (BPQ-G3,4). Each of the oligonucleotide adduct consisted of two conformers, namely syn and anti conformations, depending on the rotation around the glycosidic bond between BPQ and the guanine base. The results revealed that the BPQ moiety was located in the major groove for all four syn conformers. The relative energies of these conformers were high, and the backbone largely deviated from the B-form. On the other hand, BPQ was located in the minor groove with relatively low energies, and backbone was retained in all of the anti conformer cases. The most conceivable BPQ-modified double stranded oligonucleotide structure was proposed from the energy calculation and the structural analysis.

RNA stability under different combinations of amber force fields and solvation models

The proper matching of force field and solvent is critical to obtain correct result in molecular dynamics simulation of bio-molecules. This problem has been intensively investigated for protein but not for RNA yet. In this paper, we use standard molecular dynamics and replica exchange molecular dynamics to take a series of tests on the RNA stability under different combinations of Amber force field parameters (ff98, ff99 and ff99bsc0) and the general Born implicit solvent models (igb1, igb2 and igb5). It is found that only ff98 and ff99bsc0 with igb1 can keep the native conformations of RNA hairpin and duplex. Our results suggest that ff98 plus igb1 may be reasonable choice for molecular dynamics simulation of RNA dynamics.

A conformationally constrained peptidomimetic binds to the extracellular region of HER2 protein

Human epidermal growth factor receptor 2 (HER2) is a member of the human epidermal growth factor receptor kinases (other members include EGFR or HER1, HER3, and HER4) that are involved in signaling cascades for cell growth and differentiation. It is well established that HER2-mediated heterodimerization has important implications in cancer. Deregulation of signaling pathways and overexpression of HER2 is known to occur in cancer cells, indicating a role of HER2 in tumorigenesis. Therefore, blocking HER2-mediated signaling has potential therapeutic value. We have designed several peptidomimetics to inhibit HER2-mediated signaling for cell growth. One of the compounds (HERP5, Arg-beta Naph-Phe) exhibited antiproliferative activity with IC(50) values in the micromolar-to-nanomolar range in breast cancer cell lines. Binding of fluorescently labeled HERP5 to HER2 protein was evaluated by fluorescence assay, microscopy, and circular dichroism spectroscopy. Results indicated that HERP5 binds to the extracellular region of the HER2 protein. Structure of the peptidomimetic HERP5 was studied by NMR and molecular dynamics simulations. Based on these results a model was proposed for HER2-EGFR dimerization and possible blocking by HERP5 peptidomimetic using a protein-protein docking method.

Influence of the acetylcholinesterase active site protonation on omega loop and active site dynamics

Existence of alternative entrances in acetylcholinesterase (AChE) could explain the contrast between the very high AChE catalytic efficiency and the narrow and long access path to the active site revealed by X-ray crystallography. Alternative entrances could facilitate diffusion of the reaction products or at least water and ions from the active site. Previous molecular dynamics simulations identified side door and back door as the most probable alternative entrances. The simulations of non-inhibited AChE suggested that the back door opening events occur only rarely (0.8% of the time in the 10ns trajectory). Here we present a molecular dynamics simulation of non-inhibited AChE, where the back door opening appears much more often (14% of the time in the 12ns trajectory) and where the side door opening was observed quite frequently (78% of trajectory time). We also present molecular dynamics, where the back door does not open at all, or where large conformational changes of the AChE omega loop occur together with alternative passage opening events. All these differences in AChE dynamical behavior are caused by different protonation states of two glutamate residues located on bottom of the active site gorge (Glu202 and G450 in Mus musculus AChE). Our results confirm the results of previous molecular dynamics simulations, expand the view and suggest the probable reasons for the overall conformational behavior of AChE omega loop.

Anti-infectious bronchitis virus (IBV) activity of 1,8-cineole: effect on nucleocapsid (N) protein

In the present study, anti-IBV (infectious bronchitis virus) activity of 1,8-cineole was studied by MTT assay, as well as docking and molecular dynamic (MD) simulations. The CC50 of 1,8-cineole was above 10 mM. And the maximum noncytotoxic concentration (TD0) of 1,8-cineole was determined to be 3.90 ± 0.22 mM, which was much higher than that of ribavirin (0.78 ± 0.15 mM). 1,8-cineole could inhibit IBV with an IC(50) of 0.61 mM. MTT assay showed that the inhibition of IBV by 1, 8-cineole appears to occur moderately before entering the cell but much strongly after penetration of the virus into the cell. In silico simulations indicated that the binding site of 1,8-cineole was located at the N terminus of phosphorylated nucleocapsid (N) protein, with interaction energy equaling -40.33 kcal mol(-1). The residues TyrA92, ProA134, PheA137, AspA138 and TyrA140 had important roles during the binding process and are fully or partially conserved in various IBV strains. Based on spatial and energetic criteria, 1,8-cineole interfered with the binding between RNA and IBV N-protein. Results presented here may suggest that 1,8-cineole possesses anti-IBV properties, and therefore is a potential source of anti-IBV ingredients for the pharmaceutical industry.

A novel strategy for designing the selective PPAR agonist by the "sum of activity" model

Peroxisome proliferator-activated receptors alpha, delta and gamma are a collection of ligand-activated transcription factors crucial in lipid and glucose homeostasis. The involvement of these receptors in lipid metabolism makes them perfect therapeutic target for treating obesity and stroke. In this study, 'sum of activity' model was employed to design multi-target agonists. We used a new strategy to design agonists that fit both alpha and delta but not gamma, to avoid side effect. The CoMFA and CoMSIA models were used to explore the pharmacophore features by constructing three individual models: (a) alpha-model, (b) delta-model and (c) gamma-model, and two sum models: (d) alpha, delta- model, and (e) alpha, delta and gamma-model. The CoMFA model yielded a significant cross validation value, q(2), of 0.729 and non-cross validation value, r(2), of 0.933 in the alpha, delta-model. The CoMSIA studies yielded the best predictive models with q(2) of 0.622 in A+S and with r(2) of 0.911 in the alpha, delta-model. Finally, we proposed that distinct features shown in models (a), (b), (d) but not (c) and (e) should be accounted in designing weight-controlling drugs.

A theoretical elucidation of glucose interaction with HSA's domains

The interaction of different domains belonging to Human Serum Albumin (HSA) with open form of glucose have been investigated using molecular dynamics simulation methods. Applying docking, primary structures involving interaction of some residues with glucose have been obtained. Subsequently, equilibrium geometries at 300 K and minimum geometries have been determined for each of aforementioned structures by employing MD simulation and simulated annealing. The stability of species has been evaluated using a SAWSA v2.0 model. Ultimately, NBO analysis has been carried out to specify possible hydrogen bonding regarding the HSA interaction with glucose. Results obtained show that glucose can interact with Lys195, Lys199, and Glu153. In these interactions, each lysine forms an H-bonding with glucose. The H-bonding is obtained by stretching of N-H bond belonging to NH(3)(+) group of lysine along an oxygen atom of glucose. In addition, the above mentioned lysines are protonated, and there is an electrostatic interaction between glucose with Lys195 or Lys199. In addition, an H-bonding is formed between O atom of -COO group belonging to Glu153 and H atom of OH group belonging to glucose. Because, the N-H group of Lys195 interacts with the O atom of latter OH group, reaction of Lys195 is more desirable than that of Lys199. In fact, glucose is placed in the vicinity of Lys195 along with electrostatic interaction and H-bonding to Lys195 and Lys199 as well as H-bonding with Glu153, which subsequently reacts with Lys195. Thus, Lys195 is the primary site in reaction of glucose with HSA.

Insight derived from molecular dynamics simulation into substrate-induced changes in protein motions of proteinase K

Because of the significant industrial, agricultural and biotechnological importance of serine protease proteinase K, it has been extensively investigated using experimental approaches such as X-ray crystallography, site-directed mutagenesis and kinetic measurement. However, detailed aspects of enzymatic mechanism such as substrate binding, release and relevant regulation remain unstudied. Molecular dynamics (MD) simulations of the proteinase K alone and in complex with the peptide substrate AAPA were performed to investigate the effect of substrate binding on the dynamics/molecular motions of proteinase K. The results indicate that during simulations the substrate-complexed proteinase K adopt a more compact and stable conformation than the substrate-free form. Further essential dynamics (ED) analysis reveals that the major internal motions are confined within a subspace of very small dimension. Upon substrate binding, the overall flexibility of the protease is reduced; and the noticeable displacements are observed not only in substrate-binding regions but also in regions opposite the substrate-binding groove/pockets. The dynamic pockets caused by the large concerted motions are proposed to be linked to the substrate recognition, binding, orientation and product release; and the significant displacements in regions opposite the binding groove/pockets are considered to play a role in modulating the dynamics of enzyme-substrate interaction. Our simulation results complement the biochemical and structural studies, highlighting the dynamic mechanism of the functional properties of proteinase K.