Destabilization of Insulin Hexamer in Water–Ethanol Binary Mixture

Counteractive Effects of Choline Geranate (CAGE) ILs and Ethanol on Insulin's Stability-A Leap Forward towards Oral Insulin Formulation

Choline geranate (CAGE) ionic liquids (ILs) stabilize insulin, thereby aiding its oral delivery, whereas ethanol (EtOH) affects its stability by disrupting the hydrophobic interactions. In this study, cognizance of the stabilization mechanism of insulin dimer in the presence of both CAGE ILs and EtOH mixtures is achieved through biased and unbiased molecular dynamics (MD) simulations. Here, two order parameters are employed to study the insulin dimer dissociation using well-tempered metadynamics (WT-MetaD). The stability of insulin is found to be strongly maintained until a 0.20 mole fraction of EtOH. Besides, higher concentrations of EtOH marginally affect the insulin stability. Moreover, geranate anions form a higher number of H-bonding interactions with water molecules, which aids insulin stabilization. Conversely, the addition of EtOH minimizes the water-mediated H-bonding interactions of geranate. Additionally, geranate traps the EtOH molecules, thereby preventing the interactions between insulin and EtOH. Furthermore, the free energy landscape (FEL) reveals the absence of dimer dissociation along with noticeable deviations in the distances R and the number of contacts Q. The dimerization free energy of insulin was calculated to be −16.1 kcal/mol at a 0.20 mole fraction of EtOH. Moreover, increments in mole fractions of EtOH effectuate a decrease in the insulin stability. Thus, the present study represents CAGE ILs as efficient insulin dimer stabilizes at low concentrations of EtOH.

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.

Alcohol-perturbed self-assembly of the tobacco mosaic virus coat protein

The self-assembly of the tobacco mosaic virus coat protein is significantly altered in alcohol-water mixtures. Alcohol cosolvents stabilize the disk aggregate and prevent the formation of helical rods at low pH. A high alcohol content favours stacked disk assemblies and large rafts, while a low alcohol concentration favours individual disks and short stacks. These effects appear to be caused by the hydrophobicity of the alcohol additive, with isopropyl alcohol having the strongest effect and methanol the weakest. We discuss several effects that may contribute to preventing the protein-protein interactions between disks that are necessary to form helical rods.

Molecular Picture of the Effect of Cosolvent Crowding on Ligand Binding and Dispersed Solvation Dynamics in G-Quadruplex DNA

Understanding molecular interactions and dynamics of proteins and DNA in a cell-like crowded environment is crucial for predicting their functions within the cell. Noncanonical G-quadruplex DNA (GqDNA) structures adopt various topologies that were shown to be strongly affected by molecular crowding. However, it is unknown how such crowding affects the solvation dynamics in GqDNA. Here, we study the effect of cosolvent (acetonitrile) crowding on ligand (DAPI) solvation dynamics within human telomeric antiparallel GqDNA through direct comparison of time-resolved fluorescence Stokes shift (TRFSS) experiments and molecular dynamics (MD) simulations results. We show that ligand binding affinity to GqDNA is drastically affected by acetonitrile (ACN). Solvation dynamics probed by DAPI in GqDNA groove show dispersed dynamics from ∼100 fs to 10 ns in the absence and presence of 20% and 40% (v/v) ACN. The nature of dynamics remain similar in buffer and 20% ACN, although in 40% ACN, distinct dynamics is observed in <100 ps. MD simulations performed on GqDNA/DAPI complex reveal preferential solvation of ligand by ACN, particularly in 40% ACN. Simulated solvation time-correlation functions calculated from MD trajectories compare very well to the overall solvation dynamics of DAPI in GqDNA, observed in experiments. Linear response decomposition of simulated solvation correlation functions unfolds the origin of dispersed dynamics, showing that the slower dynamics is dominated by DNA-motion in the presence of ACN (and also by the ACN dynamics at higher concentration). However, water-DNA coupled motion controls the slow dynamics in the absence of ACN. Our data, thus, unravel a detailed molecular picture showing that though ACN crowding affect ligand binding affinity to GqDNA significantly, the overall dispersed solvation dynamics in GqDNA remain similar in the absence and the presence of 20% ACN, albeit with a small effect on the dynamics in the presence of 40% ACN due to preferential solvation of ligand by ACN.

Kinetics of Phenol Escape from the Insulin R6 Hexamer

Therapeutic preparations of insulin often contain phenolic molecules, which can impact both pharmacokinetics and shelf life. Thus, understanding the interactions of insulin and phenolic molecules can aid in designing improved therapeutics. In this study, we use molecular dynamics to investigate phenol release from the insulin hexamer. Leveraging recent advances in methods for analyzing molecular dynamics data, we expand on existing simulation studies to identify and quantitatively characterize six phenol binding/unbinding pathways for wild-type and A10 Ile → Val and B13 Glu → Gln mutant insulins. A number of these pathways involve large-scale opening of the primary escape channel, suggesting that the hexamer is much more dynamic than previously appreciated. We show that phenol unbinding is a multipathway process, with no single pathway representing more than 50% of the reactive current and all pathways representing at least 10%. We use the mutant simulations to show how the contributions of specific pathways can be rationally manipulated. Predicting the net effects of mutations is more challenging because the kinetics depend on all of the pathways, demanding quantitatively accurate simulations and experiments.

Rate of Insulin Dimer Dissociation: Interplay between Memory Effects and Higher Dimensionality

We calculate the rate of dissociation of an insulin dimer into two monomers in water. The rate of this complex reaction is determined by multiple factors that are elucidated. By employing advanced sampling techniques, we first obtain the reaction free energy surface for the dimer dissociation as a function of two order parameters, namely, the distance between the center-of-mass of two monomers (R) and the number of cross-contacts (Q) among the backbone C_α atoms of two monomers. We then construct an orthogonal 2D reaction energy surface by introducing the reaction coordinate X to denote the minimum energy pathway and a conjugate coordinate Y that spans the orthogonal direction. The free energy landscape is rugged with multiple maxima and minima. We calculate the rate by employing not only the non-Markovian multidimensional rate theory but also several other theoretical approaches. The necessary reaction frequencies and the frictions are calculated from the time correlation function formalism. Our best estimate of the rate is 0.4 μs^-1. Our study reveals interesting opposite influences of dimensionality and memory in determining the rate constant of the reaction. We gain interesting insights into the dimer dissociation process by looking directly at the trajectories obtained from molecular dynamics simulation.

Unfolding of Dynamical Events in the Early Stage of Insulin Dimer Dissociation

The dissociation of an insulin dimer is an important biochemical event that could also serve as a prototype of dissociations in similar biomolecular assemblies. We use a recently developed multidimensional free energy landscape for insulin dimer dissociation to unearth the microscopic and mechanistic aspects of the initial stages of the process that could hold the key to understanding the stability and the rate. The following sequence of events occurs in the initial stages: (i) The backbone hydrogen bonds break partially at the antiparallel β-sheet junction, (ii) the two α-helices (chain B) move away from each other while several residues (chain A) move closer, and (iii) a flow of adjacent water molecules occurs into the junction region. Interestingly, the intermonomeric center-to-center distance does not increase, but the number of native contacts exhibits a sharp decrease. Subsequent steps involve further disengagement of hydrophobic groups. This process is slow because of an entropic bottleneck created by the existence of the large configuration space available in the native state (NS), which is inhabited by low-frequency conformational fluctuations. We carry out a density-of-states analyses in the dimer NS to unearth distinctive features not present in the monomers. These low-frequency modes are also responsible for a large entropic stabilization of the NS. Hydrophobic disengagement in the early stage leads to the formation of a twisted intermediate state which itself is a metastable minimum (IS-1). The subsequent progress leads to another dimeric complex (IS-2), which is on the dissociative pathway and characterized by a further decrease in the native contacts. The dissociation process provides insights into the workings of a biomolecular assembly.

Water in biodegradable glucose–water–urea deep eutectic solvent: modifications of structure and dynamics in a crowded environment

Molecular dynamics simulations have been performed on a highly viscous (η ∼ 255 cP) naturally abundant deep eutectic solvent (NADES) composed of glucose, urea and water in a weight ratio of 6 : 4 : 1 at 328 K.

Modulation of Membrane Fluidity to Control Interfacial Water Structure and Dynamics in Saturated and Unsaturated Phospholipid Vesicles

The structure and dynamics of interfacial water in biological systems regulate the biochemical reactions. But, it is still enigmatic how the behavior of the interfacial water molecule is controlled. Here, we have investigated the effect of membrane fluidity on the structure and dynamics of interfacial water molecules in biologically relevant phopholipid vesicles. This study delineates that modulation of membrane fluidity through interlipid separation and unsaturation not only mitigate membrane rigidity but also disrupt the strong hydrogen bond (H-bond) network around the lipid bilayer interface. As a result, a disorder in H-bonding between water molecules arises several layers beyond the first hydration shell of the polar headgroup, which essentially modifies the interfacial water structure and dynamics. Furthermore, we have also provided evidence of increasing transportation through these modulated membranes, which enhance the membrane mediated isomerization reaction rate.