Investigation of glucose binding sites on insulin

Understanding the Reversible Binding of a Multichain Protein-Protein Complex through Free-Energy Calculations

We demonstrate that the binding affinity of a multichain protein-protein complex, insulin dimer, can be accurately predicted using a streamlined route of standard binding free-energy calculations. We find that chains A and C, which do not interact directly during binding, stabilize the insulin monomer structures and reduce the binding affinity of the two monomers, therefore enabling their reversible association. Notably, we confirm that although classical methods can estimate the binding affinity of the insulin dimer, conventional molecular dynamics, enhanced sampling algorithms, and classical geometrical routes of binding free-energy calculations may not fully capture certain aspects of the role played by the noninteracting chains in the binding dynamics. Therefore, this study not only elucidates the role of noninteracting chains in the reversible binding of the insulin dimer but also offers a methodological guide for investigating the reversible binding of multichain protein-protein complexes utilizing streamlined free-energy calculations.

Insulin facilitates corneal wound healing in the diabetic environment through the RTK-PI3K/Akt/mTOR axis in vitro

Diabetic patients can develop degenerative corneal changes, termed diabetic keratopathy, during the course of their disease. Topical insulin has been shown to reduce corneal wound area and restore sensitivity in diabetic rats, and both the insulin receptor (IR) and insulin-like growth factor 1 receptor (IGF-1R) stimulate cell signaling of the PI3K-Akt pathway. The purpose of this study was to assess a mechanism by which improved wound healing occurs by characterizing expression within the PI3K-Akt pathway in corneal epithelial and stromal cells. In vitro scratch tests were used to evaluate wound healing outcomes under variable glucose conditions in the presence or absence of insulin. Protein expression of intracellular kinases in the PI3K pathway, stromal cell markers, and GLUT-1 was evaluated by immunoblotting.TGF-β1 expression was evaluated by ELISA. Insulin promoted in vitro wound healing in all cell types. In human corneal epithelial cells, insulin did not induce PI3K-Akt signaling; however, in all other cell types evaluated, insulin increased expression of PI3K-Akt signaling proteins compared to vehicle control. Fibroblasts variably expressed α-SMA under all treatment conditions, with significant increases in α-SMA and TGF-β1 occurring in a dose-dependent manner with glucose concentration. These results indicate that insulin can promote corneal cellular migration and proliferation by inducing Akt signaling. Exogenous insulin therapy may serve as a novel target of therapeutic intervention for diabetic keratopathy.

Progress in Simulation Studies of Insulin Structure and Function

Insulin is a peptide hormone known for chiefly regulating glucose level in blood among several other metabolic processes. Insulin remains the most effective drug for treating diabetes mellitus. Insulin is synthesized in the pancreatic β-cells where it exists in a compact hexameric architecture although its biologically active form is monomeric. Insulin exhibits a sequence of conformational variations during the transition from the hexamer state to its biologically-active monomer state. The structural transitions and the mechanism of action of insulin have been investigated using several experimental and computational methods. This review primarily highlights the contributions of molecular dynamics (MD) simulations in elucidating the atomic-level details of conformational dynamics in insulin, where the structure of the hormone has been probed as a monomer, dimer, and hexamer. The effect of solvent, pH, temperature, and pressure have been probed at the microscopic scale. Given the focus of this review on the structure of the hormone, simulation studies involving interactions between the hormone and its receptor are only briefly highlighted, and studies on other related peptides (e.g., insulin-like growth factors) are not discussed. However, the review highlights conformational dynamics underlying the activities of reported insulin analogs and mimetics. The future prospects for computational methods in developing promising synthetic insulin analogs are also briefly highlighted.

Discovery of Flazin, an Alkaloid Isolated from Cherry Tomato Juice, As a Novel Non-Enzymatic Protein Glycation Inhibitor via in Vitro and in Silico Studies

Both overproduced reactive oxygen species/reactive nitrogen species and hyperglycemic conditions accompany a significant increase in protein glycation and nitration that contribute to the initiation and progression of diabetic complications and neuronal disorders. In this study, 19 compounds, including steroidal saponins, alkaloids, cerebroside, phenolic compounds, sterols, and nucleosides, were isolated from cherry tomato (Solanum lycopersicum var. cerasiforme) juice, of which flazin showed good inhibition on monosaccharide-induced non-enzymatic bovine pancreas insulin and bovine serum albumin (BSA) glycation. Molecular dynamics simulations revealed that flazin continuously interacts with Phe1, Val2, Tyr26, and Lys29 insulin residues, which play a key role in insulin glycation/dimerization. In addition, depending upon the flazin dose, this blocked the tyrosine nitration of BSA via scavenging peroxynitrite anions. Taken together, our novel findings suggest that flazin could be a lead compound for the treatment of diabetes and neuronal disorders via the inhibition of non-enzymatic protein glycation and the elimination of peroxynitrite.

Identifying Phlorofucofuroeckol-A as a Dual Inhibitor of Amyloid-β25-35 Self-Aggregation and Insulin Glycation: Elucidation of the Molecular Mechanism of Action

Both amyloid-β (Aβ) and insulin are amyloidogenic peptides, and they play a critical role in Alzheimer's disease (AD) and type-2 diabetes (T2D). Misfolded or aggregated Aβ and glycated insulin are commonly found in AD and T2D patients, respectively, and exhibit neurotoxicity and oxidative stress. The present study examined the anti-Aβ25-35 aggregation and anti-insulin glycation activities of five phlorotannins isolated from Ecklonia stolonifera. Thioflavin-T assay results suggest that eckol, dioxinodehydroeckol, dieckol, and phlorofucofuroeckol-A (PFFA) significantly inhibit Aβ25-35 self-assembly. Molecular docking and dynamic simulation analyses confirmed that these phlorotannins have a strong potential to interact with Aβ25-35 peptides and interrupt their self-assembly and conformational transformation, thereby inhibiting Aβ25-35 aggregation. In addition, PFFA dose-dependently inhibited d-ribose and d-glucose induced non-enzymatic insulin glycation. To understand the molecular mechanism for insulin glycation and its inhibition, we predicted the binding site of PFFA in insulin via computational analysis. Interestingly, PFFA strongly interacted with the Phe1 in insulin chain-B, and this interaction could block d-glucose access to the glycation site of insulin. Taken together, our novel findings suggest that phlorofucofuroeckol-A could be a new scaffold for AD treatment by inhibiting the formation of β-sheet rich structures in Aβ25-35 and advanced glycation end-products (AGEs) in insulin.

Tracking Insulin Glycation in Real Time by Time-Resolved Emission Spectroscopy

The application of time-resolved fluorescence sensing to the study of heterogenic biomolecular systems remains challenging because of the complexity of the resulting photophysics. Measuring the time-resolved emission spectroscopy (TRES) spectra can provide a more informative alternative to the modeling of the fluorescence decay that is currently employed. Here, we demonstrate this approach by monitoring real-time changes in intrinsic insulin fluorescence by TRES as a straightforward probe to directly measure kinetics of insulin aggregation and glycation. Our findings hold promise for monitoring the storage of insulin and its application in the control of diabetes and may support the development of more effective therapeutics against amyloidosis.

Targeting the Pentose Phosphate Pathway: Characterization of a New 6PGL Inhibitor

Human African trypanosomiasis, or sleeping sickness, is a lethal disease caused by the protozoan parasite Trypanosoma brucei. However, although many efforts have been made to understand the biochemistry of this parasite, drug development has led to treatments that are of limited efficiency and of great toxicity. To develop new drugs, new targets must be identified, and among the several metabolic processes of trypanosomes that have been proposed as drug targets, carbohydrate metabolism (glycolysis and the pentose phosphate pathway (PPP)) appears as a promising one. As far as the PPP is concerned, a limited number of studies are related to the glucose-6-phosphate dehydrogenase. In this work, we have focused on the activity of the second PPP enzyme (6-phospho-gluconolactonase (6PGL)) that transforms 6-phosphogluconolactone into 6-phosphogluconic acid. A lactam analog of the natural substrate has been synthesized, and binding of the ligand to 6PGL has been investigated by NMR titration. The ability of this ligand to inhibit 6PGL has also been demonstrated using ultraviolet experiments, and protein-inhibitor interactions have been investigated through docking calculations and molecular dynamics simulations. In addition, a marginal inhibition of the third enzyme of the PPP (6-phosphogluconate dehydrogenase) was also demonstrated. Our results thus open new prospects for targeting T. brucei.

Rapid Non-Enzymatic Glycation of the Insulin Receptor under Hyperglycemic Conditions Inhibits Insulin Binding In Vitro: Implications for Insulin Resistance

The causes of insulin resistance are not well-understood in either type 1 or type 2 diabetes. Insulin (INS) is known to undergo rapid non-enzymatic covalent conjugation to glucose or other sugars (glycation). Because the insulin receptor (IR) has INS-like regions associated with both glucose and INS binding, we hypothesize that hyperglycemic conditions may rapidly glycate the IR, chronically interfering with INS binding. IR peptides were synthesized spanning IR- associated INS-binding regions. Glycation rates of peptides under hyperglycemic conditions were followed over six days using matrix assisted laser desorption/ionization-time of flight (MALDI-TOF) mass spectrometry. INS conjugated to horse-radish peroxidase was used to determine INS binding to IR peptides in glycated and non-glycated forms. Several IR peptides were glycated up to 14% within days of exposure to 20-60 mM glucose. Rates of IR-peptide glycation were comparable to those of insulin. Glycation of four IR peptides significantly inhibits INS binding to them. Glycation of intact IR also decreases INS binding by about a third, although it was not possible to confirm the glycation sites on the intact IR. Glycation of the IR may therefore provide a mechanism by which INS resistance develops in diabetes. Demonstration of glycation of intact IR in vivo is needed.

An Insulin-Like Modular Basis for the Evolution of Glucose Transporters (GLUT) with Implications for Diabetes

Glucose transporters (GLUT) are twelve-transmembrane spanning proteins that contain two pores capable of transporting glucose and dehydroascorbate in and out of cells. The mechanism by which transport is effected is unknown. An evolutionarily-based hypothesis for the mechanism of glucose transport is presented here based on reports that insulin has multiple binding sites for glucose. It is proposed that insulin-like peptides were incorporated as modular elements into transmembrane proteins during evolution, resulting in glucose transporting capacity. Homology searching reveals that all GLUT contain multiple copies of insulin-like regions. These regions map onto a model of GLUT in positions that define the glucose transport cores. This observation provides a mechanism for glucose transport involving the diffusion of glucose from one insulin-like glucose-binding region to another. It also suggests a mechanism by which glucose disregulation may occur in both type 1 and type 2 diabetes: insulin rapidly self-glycates under hyperglycemic conditions. Insulin-like regions of GLUT may also self-glycate rapidly, thereby interfering with transport of glucose into cells and disabling GLUT sensing of blood glucose levels. All aspects of the hypothesis are experimentally testable.

Molecular recognition between insulin and dextran encapsulated gold nanoparticles

Insulin is a peptide hormone that can regulate the metabolism of carbohydrates and lipids. This hormone is closely related to glucose-uptake in cells and can control blood glucose levels. Dextran is a polysaccharide composed of glucose units. In this study, we discovered that dextran-encapsulated gold nanoparticles (AuNPs@Dextran) and nanoclusters (AuNCs@Dextran) can be used to recognize insulin. The dissociation constant of insulin toward AuNPs@Dextran was estimated to be ∼5.3 × 10^-6  M. The binding site on insulin toward the dextran on the nanoprobes was explored as well. It was found that the sequence of numbers 1-22 on the insulin B chain can interact with the dextran encapsulated nanoprobes. Additionally, we also demonstrated that the dextran-encapsulated nanoprobes could be used as concentration probes to selectively enrich trace amounts of insulin (∼1 pM) from serum samples. Copyright © 2016 John Wiley & Sons, Ltd.

A modular hierarchy-based theory of the chemical origins of life based on molecular complementarity

Albert Szent-Gyorgyi once defined discovery as seeing what everyone else sees and thinking what no one else thinks. I often find that phenomena that are obvious to other people are not obvious to me. Molecular complementarity is one of these phenomena: while rare among any random set of compounds, it is ubiquitous in living systems. Because every molecule in a living system binds more or less specifically to several others, we now speak of "interactomes". What explains the ubiquity of molecular complementarity in living systems? What might such an explanation reveal about the chemical origins of life and the principles that have governed its evolution? Beyond this, what might complementarity tell us about the optimization of integrated systems in general? My research combines theoretical and experimental approaches to molecular complementarity relating to evolution from prebiotic chemical systems to superorganismal interactions. Experimentally, I have characterized complementarity involving specific binding between small molecules and explored how these small-molecule modules have been incorporated into macromolecular systems such as receptors and transporters. Several general principles have emerged from this research. Molecules that bind to each other almost always alter each other's physiological effects; and conversely, molecules that have antagonistic or synergistic physiological effects almost always bind to each other. This principle suggests a chemical link between biological structure and function. Secondly, modern biological systems contain an embedded molecular paleontology based on complementarity that can reveal their chemical origins. This molecular paleontology is often manifested through modules involving small, molecularly complementary subunits that are built into modern macromolecular structures such as receptors and transporters. A third principle is that complementary modules are conserved and repurposed at every stage of evolution. Molecular complementarity plays critical roles in the evolution of chemical systems and resolves a significant number of outstanding problems in the emergence of complex systems. All physical and mathematical models of organization within complex systems rely upon nonrandom linkage between components. Molecular complementarity provides a naturally occurring nonrandom linker. More importantly, the formation of hierarchically organized stable modules vastly improves the probability of achieving self-organization, and molecular complementarity provides a mechanism by which hierarchically organized stable modules can form. Finally, modularity based on molecular complementarity produces a means for storing and replicating information. Linear replicating molecules such as DNA or RNA are not required to transmit information from one generation of compounds to the next: compositional replication is as ubiquitous in living systems as genetic replication and is equally important to its functions. Chemical systems composed of complementary modules mediate this compositional replication and gave rise to linear replication schemes. In sum, I propose that molecular complementarity is ubiquitous in living systems because it provides the physicochemical basis for modular, hierarchical ordering and replication necessary for the evolution of the chemical systems upon which life is based. I conjecture that complementarity more generally is an essential agent that mediates evolution at every level of organization.

Carbohydrate force fields

Carbohydrates present a special set of challenges to the generation of force fields. First, the tertiary structures of monosaccharides are complex merely by virtue of their exceptionally high number of chiral centers. In addition, their electronic characteristics lead to molecular geometries and electrostatic landscapes that can be challenging to predict and model. The monosaccharide units can also interconnect in many ways, resulting in a large number of possible oligosaccharides and polysaccharides, both linear and branched. These larger structures contain a number of rotatable bonds, meaning they potentially sample an enormous conformational space. This article briefly reviews the history of carbohydrate force fields, examining and comparing their challenges, forms, philosophies, and development strategies. Then it presents a survey of recent uses of these force fields, noting trends, strengths, deficiencies, and possible directions for future expansion.

Role of computational methods in pharmaceutical sciences

Over the past two decades computational methods have eased up the financial and experimental burden of early drug discovery process. The in silico methods have provided support in terms of databases, data mining of large genomes, network analysis, systems biology on the bioinformatics front and structure-activity relationship, similarity analysis, docking, and pharmacophore methods for lead design and optimization. This review highlights some of the applications of bioinformatics and chemoinformatics methods that have enriched the field of drug discovery. In addition, the review also provided insights into the use of free energy perturbation methods for efficiently computing binding energy. These in silico methods are complementary and can be easily integrated into the traditional in vitro and in vivo methods to test pharmacological hypothesis.

Software and resources for computational medicinal chemistry

Computer-aided drug design plays a vital role in drug discovery and development and has become an indispensable tool in the pharmaceutical industry. Computational medicinal chemists can take advantage of all kinds of software and resources in the computer-aided drug design field for the purposes of discovering and optimizing biologically active compounds. This article reviews software and other resources related to computer-aided drug design approaches, putting particular emphasis on structure-based drug design, ligand-based drug design, chemical databases and chemoinformatics tools.

Controlling the aggregation and rate of release in order to improve insulin formulation: molecular dynamics study of full-length insulin amyloid oligomer models

Insulin is a hormone that regulates the physiological glucose level in human blood. Insulin injections are used to treat diabetic patients. The amyloid aggregation of insulin may cause problems during the production, storage, and delivery of insulin formulations. Several modifications to the C-terminus of the B chain have been suggested in order to improve the insulin formulation. The central fragments of the A and B chains (LYQLENY and LVEALYL) have recently been identified as β-sheet-forming regions, and their microcrystalline structures have been used to build a high-resolution amyloid fibril model of insulin. Here we report on a molecular dynamics (MD) study of single-layer oligomers of the full-length insulin which aimed to identify the structural elements that are important for amyloid stability, and to suggest single glycine mutants in the β-sheet region that may improve the formulation. Structural stability, aggregation behavior and the thermodynamics of association were studied for the wild-type and mutant aggregates. A comparison of the oligomers of different sizes revealed that adding strands enhances the internal stability of the wild-type aggregates. We call this "dynamic cooperativity". The secondary structure content and clustering analysis of the MD trajectories show that the largest aggregates retain the fibril conformation, while the monomers and dimers lose their conformations. The degree of structural similarity between the oligomers in the simulation and the fibril conformation is proposed as a possible explanation for the experimentally observed shortening of the nucleation lag phase of insulin with oligomer seeding. Decomposing the free energy into electrostatic, van der Waals and solvation components demonstrated that electrostatic interactions contribute unfavorably to the binding, while the van der Waals and especially solvation effects are favorable for it. A per-atom decomposition allowed us to identify the residues that contribute most to the binding free energy. Residues in the β-sheet regions of chains A and B were found to be the key residues as they provided the largest favorable contributions to single-layer association. The positive ∆∆G (mut) values of 37.3 to 1.4 kcal mol(-1) of the mutants in the β-sheet region indicate that they have a lower tendency to aggregate than the wild type. The information obtained by identifying the parts of insulin molecules that are crucial to aggregate formation and stability can be used to design new analogs that can better control the blood glucose level. The results of our simulation may help in the rational design of new insulin analogs with a decreased propensity for self-association, thus avoiding injection amyloidosis. They may also be used to design new fast-acting and delayed-release insulin formulations.

Glucose binds to the insulin receptor affecting the mutual affinity of insulin and its receptor

Insulin activity is sensitive to glucose concentration but the mechanisms are still unclear. An unexamined possibility is that the insulin receptor (IR) is sensitive to glucose concentration. We demonstrate here that insulin-like peptides derived from the IR bind glucose at low millimolar, and cytochalasin B at low micromolar, concentrations; several insulin-like IR peptides bind insulin at nanomolar Kd; and this binding is antagonized by increasing glucose concentrations. In addition, glucose and cytochalasin B bind to IR isolated from rat liver and increasing glucose decreases insulin binding to this IR preparation. The presence of GLUT 1 in our IR preparation suggests the possibility of additional glucose-mediated allosteric control. We propose a model in which glucose binds to insulin, the IR, and GLUT; insulin binds to the IR; and the IR binds to GLUT. This set of interactions produces an integrated system of insulin-dependent interactions that is highly sensitive to glucose concentration.

Collective motions in glucosamine-6-phosphate synthase: influence of ligand binding and role in ammonia channelling and opening of the fructose-6-phosphate binding site

The large protein motions of the bacterial enzyme glucosamine-6-phosphate synthase have been addressed using full atom normal modes analysis for the empty, the glucose-6-phosphate and the glucose-6-phosphate+glutamate bound proteins. The approach that was used involving energy minimizations along the normal modes coordinates identified functional motions of the protein, some of which were characterized earlier by X-ray diffraction studies. This method made it possible for the first time to highlight significant energy differences according to whether none, only one or both of the active sites of the protein were occupied. Our data favoured a specific motion of the glutamine binding domain following the fixation of fructose-6-phosphate and suggested a rigidified structure with both sites occupied. Here, we show that most of the collective large amplitude motions of glucosamine-6-phosphate synthase that are modulated by ligand binding are crucial for the enzyme catalytic cycle, as they strongly modify the geometry of both the ammonia channel and the C-tail, demonstrating their role in ammonia transfer and ligand binding.

Insulin: a model system for nanomedicine?

One of the predominant aims of insulin therapy for diabetes is to appropriately mimic physiological insulin secretion levels and their correlation with glucose concentration in healthy individuals. This report outlines current methods and their limitations in glycemic control and their possible relationship to insufficient knowledge about the structure and dynamics of the insulin hormone itself. Based on recent experimental and computational work, a possible approach to less-invasive insulin administration is sketched.

Ammonia channeling in bacterial glucosamine-6-phosphate synthase (Glms): molecular dynamics simulations and kinetic studies of protein mutants

Ammonia transfer from the glutamine site to the fructose-6P site of bacterial glucosamine-6-phosphate synthase was studied by molecular dynamics simulations. The studies suggest a key role for Trp74, in the sealing of the hydrophobic channel connecting the two binding sites, as well as for the two Ala602 and Val605 residues, which form a narrow passage whose opening/closing constitutes an essential event in ammonia transfer. Kinetic analyses of the corresponding protein mutants confirmed our predictions. The efficiency of ammonia transfer which was close to zero in the W74A mutant was partially restored by increasing the size of the corresponding side-chain; the simulations performed on the W74A mutant suggested the formation of a hole in the channel. In the case of A602L and V605L mutants, the efficiency of ammonia transfer decreased to approximately 50% of the value of the native protein. None of the mutants were, however, able to use exogenous ammonia as a substrate.

Compositional complementarity and prebiotic ecology in the origin of life

We hypothesize that life began not with the first self-reproducing molecule or metabolic network, but as a prebiotic ecology of co-evolving populations of macromolecular aggregates (composomes). Each composome species had a particular molecular composition resulting from molecular complementarity among environmentally available prebiotic compounds. Natural selection acted on composomal species that varied in properties and functions such as stability, catalysis, fission, fusion and selective accumulation of molecules from solution. Fission permitted molecular replication based on composition rather than linear structure, while fusion created composomal variability. Catalytic functions provided additional chemical novelty resulting eventually in autocatalytic and mutually catalytic networks within composomal species. Composomal autocatalysis and interdependence allowed the Darwinian co-evolution of content and control (metabolism). The existence of chemical interfaces within complex composomes created linear templates upon which self-reproducing molecules (such as RNA) could be synthesized, permitting the evolution of informational replication by molecular templating. Mathematical and experimental tests are proposed.

Combining docking and molecular dynamic simulations in drug design

A rational approach is needed to maximize the chances of finding new drugs, and to exploit the opportunities of potential new drug targets emerging from genomic and proteomic initiatives, and from the large libraries of small compounds now readily available through combinatorial chemistry. Despite a shaky early history, computer-aided drug design techniques can now be effective in reducing costs and speeding up drug discovery. This happy outcome results from development of more accurate and reliable algorithms, use of more thoughtfully planned strategies to apply them, and greatly increased computer power to allow studies with the necessary reliability to be performed. Our review focuses on applications and protocols, with the main emphasis on critical analysis of recent studies where docking calculations and molecular dynamics (MD) simulations were combined to dock small molecules into protein receptors. We highlight successes to demonstrate what is possible now, but also point out drawbacks and future directions. The review is structured to lead the reader from the simpler to more compute-intensive methods. Thus, while inexpensive and fast docking algorithms can be used to scan large compound libraries and reduce their size, more accurate but expensive MD simulations can be applied when a few selected ligand candidates remain. MD simulations can be used: during the preparation of the protein receptor before docking, to optimize its structure and account for protein flexibility; for the refinement of docked complexes, to include solvent effects and account for induced fit; to calculate binding free energies, to provide an accurate ranking of the potential ligands; and in the latest developments, during the docking process itself to find the binding site and correctly dock the ligand a priori.

Multiple ligand-binding modes in bacterial R67 dihydrofolate reductase

R67 dihydrofolate reductase (DHFR), a bacterial plasmid-encoded enzyme associated with resistance to the drug trimethoprim, shows neither sequence nor structural homology with the chromosomal DHFR. It presents a highly symmetrical toroidal structure, where four identical monomers contribute to the unique central active-site pore. Two reactants (dihydrofolate, DHF), two cofactors (NADPH) or one of each (R67*DHF*NADPH) can be found simultaneously within the active site, the last one being the reactive ternary complex. As the positioning of the ligands has proven elusive to empirical determination, we addressed the problem from a theoretical perspective. Several potential structures of the ternary complex were generated using the docking programs AutoDock and FlexX. The variability among the final poses, many of which conformed to experimental data, prompted us to perform a comparative scoring analysis and molecular dynamics simulations to assess the stability of the complexes. Analysis of ligand-ligand and ligand-protein interactions along the 4 ns trajectories of eight different structures allowed us to identify important inter-ligand contacts and key protein residues. Our results, combined with published empirical data, clearly suggest that multipe binding modes of the ligands are possible within R67 DHFR. While the pterin ring of DHF and the nicotinamide ring of NADPH assume a stacked endo-conformation at the centre of the pore, probably assisted by V66, Q67 and I68, the tails of the molecules extend towards opposite ends of the cavity, adopting multiple configurations in a solvent rich-environment where hydrogen-bond interactions with K32 and Y69 may play important roles.

Investigation of glucose binding sites on insulin

Possible insulin binding sites for D-glucose have been investigated theoretically by docking and molecular dynamics (MD) simulations. Two different docking programs for small molecules were used; Multiple Copy Simultaneous Search (MCSS) and Solvation Energy for Exhaustive Docking (SEED) programs. The configurations resulting from the MCSS search were evaluated with a scoring function developed to estimate the binding free energy. SEED calculations were performed using various values for the dielectric constant of the solute. It is found that scores emphasizing non-polar interactions gave a preferential binding site in agreement with that inferred from recent fluorescence and NMR NOESY experiments. The calculated binding affinity of -1.4 to -3.5 kcal/mol is within the measured range of -2.0 +/- 0.5 kcal/mol. The validity of the binding site is suggested by the dynamical stability of the bound glucose when examined with MD simulations with explicit solvent. Alternative binding sites were found in the simulations and their relative stabilities were estimated. The motions of the bound glucose during molecular dynamics simulations are correlated with the motions of the insulin side chains that are in contact with it and with larger scale insulin motions. These results raise the question of whether glucose binding to insulin could play a role in its activity. The results establish the complementarity of molecular dynamics simulations and normal mode analyses with the search for binding sites proposed with small molecule docking programs.