Insulin dimer dissociation in aqueous solution: A computational study of free energy landscape and evolving microscopic structure along the reaction pathway

Saponins Effect on Human Insulin Amyloid Aggregation

The misfolding and amyloid aggregation of proteins have been attracting scientific interest for a few decades, due to their link with several diseases, particularly neurodegenerative diseases. Proteins can assemble and result in insoluble aggregates that, together with intermediate oligomeric species, modify the extracellular environment. Many efforts have been and are devoted to the search for cosolvents and cosolutes able to interfere with amyloid aggregation. In this work, we intensively study the effect of saponins, bioactive compounds, on human insulin aggregation. To monitor the kinetic of amyloid aggregation following secondary structure changes, we perform fluorescence and UV-Visible absorption spectroscopies, using Thioflavin T and Congo Red as amyloid specific probes, and Circular Dichroism. To study the overall structural features and size of aggregates, we perform Synchrotron Small-Angle X-ray Scattering and Dynamic Light Scattering experiments. The morphology of the aggregates was assessed by Atomic Force Microscopy. To deepen the understanding of the saponins interaction with insulin, a Molecular Dynamics investigation is performed, too. The reported data demonstrate that saponins interfere with the amyloid aggregation by inducing a strong inhibition on the formation of insulin fibrils, likely through specific interactions with insulin monomers. A dose-dependent effect is evident, and amyloid inhibition is already clear when saponins are just 0.01% w/w in solution. We suggest that saponins, which are natural metabolites present in a wide range of foods ranging from grains, pulses, and green leaves to sea stars and cucumbers, can be promising metabolites to inhibit human insulin aggregation. This basic research work can pave the way to further investigations concerning insulin amyloidosis, suggesting the use of saponins as amyloid inhibitors and/or stabilizing agents in solution.

Analysis of the Dynamics of a Complex, Multipathway Reaction: Insulin Dimer Dissociation

The protein hormone insulin forms a homodimer that must dissociate to bind to its receptor. Understanding the kinetics and mechanism of dissociation is essential for the rational design of therapeutic analogs. In addition to its physiological importance, this dissociation process serves as a paradigm for coupled (un)folding and (un)binding. Based on previous free energy simulations, insulin dissociation is thought to involve multiple pathways with comparable free energy barriers. Here, we analyze the mechanism of insulin dimer dissociation using a recently developed computational framework for estimating kinetic statistics from short-trajectory data. These statistics indicate that the likelihood of dissociation (the committor) closely tracks the decrease in the number of (native and nonnative) intermonomer contacts and the increase in the number of water contacts at the dimer interface; the transition state with equal likelihood of association and dissociation corresponds to an encounter complex with relatively few native contacts and many nonnative contacts. We identify four pathways out of the dimer state and quantify their contributions to the rate as well as their exchange by computing reactive fluxes. We show that both the pathways and their extents of exchange can be understood in terms of rotations around three axes of the dimer structure. Our results provide insights into the kinetics of insulin analogs and, more generally, how to characterize complex, multipathway processes.

Analysis of the dynamics of a complex, multipathway reaction: Insulin dimer dissociation

The protein hormone insulin forms a homodimer that must dissociate to bind to its receptor. Understanding the kinetics and mechanism of dissociation is essential for rational design of therapeutic analogs. In addition to its physiological importance, this dissociation process serves as a paradigm for coupled (un)folding and (un)binding. Based on previous free energy simulations, insulin dissociation is thought to involve multiple pathways with comparable free energy barriers. Here, we analyze the mechanism of insulin dimer dissociation using a recently developed computational framework for estimating kinetic statistics from short-trajectory data. These statistics indicate that the likelihood of dissociation (the committor) closely tracks the decrease in the number of (native and nonnative) intermonomer contacts and the increase in the number of water contacts at the dimer interface; the transition state with equal likelihood of association and dissociation corresponds to an encounter complex with relatively few native contacts and many nonnative contacts. We identify four pathways out of the dimer state and quantify their contributions to the rate, as well as their exchange, by computing reactive fluxes. We show that both the pathways and their extents of exchange can be understood in terms of rotations around three axes of the dimer structure. Our results provide insights into the kinetics of insulin analogues and, more generally, how to characterize complex, multipathway processes.

Evolution of biosynthetic human insulin and its analogues for diabetes management

Hormones play a crucial role in maintaining the normal human physiology. By acting as chemical messengers that facilitate the communication between different organs, tissues and cells of the body hormones assist in responding appropriately to external and internal stimuli that trigger growth, development and metabolic activities of the body. Any abnormalities in the hormonal composition and balance can lead to devastating health consequences. Hormones have been important therapeutic agents since the early 20th century, when it was realized that their exogenous supply could serve as a functional substitution for those hormones which are not produced enough or are completely lacking, endogenously. Insulin, the pivotal anabolic hormone in the body, was used for the treatment of diabetes mellitus, a metabolic disorder due to the absence or intolerance towards insulin, since 1921 and is the trailblazer in hormone therapeutics. At present the largest market share for therapeutic hormones is held by insulin. Many other hormones were introduced into clinical practice following the success with insulin. However, for the six decades following the introduction the first therapeutic hormone, there was no reliable method for producing human hormones. The most common source for hormones were animals, although semisynthetic and synthetic hormones were also developed. However, none of these were optimal because of their allergenicity, immunogenicity, lack of consistency in purity and most importantly, scalability. The advent of recombinant DNA technology was a game changer for hormone therapeutics. This revolutionary molecular biology tool made it possible to synthesize human hormones in microbial cell factories. The approach allowed for the synthesis of highly pure hormones which were structurally and biochemically identical to the human hormones. Further, the fermentation techniques utilized to produce recombinant hormones were highly scalable. Moreover, by employing tools such as site directed mutagenesis along with recombinant DNA technology, it became possible to amend the molecular structure of the hormones to achieve better efficacy and mimic the exact physiology of the endogenous hormone. The first recombinant hormone to be deployed in clinical practice was insulin. It was called biosynthetic human insulin to reflect the biological route of production. Subsequently, the biochemistry of recombinant insulin was modified using the possibilities of recombinant DNA technology and genetic engineering to produce analogues that better mimic physiological insulin. These analogues were tailored to exhibit pharmacokinetic and pharmacodynamic properties of the prandial and basal human insulins to achieve better glycemic control. The present chapter explores the principles of genetic engineering applied to therapeutic hormones by reviewing the evolution of therapeutic insulin and its analogues. It also focuses on how recombinant analogues account for the better management of diabetes mellitus.

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.

Cooperative Membrane Binding of HIV-1 Matrix Proteins

The HIV-1 assembly process begins with a newly synthesized Gag polyprotein being targeted to the inner leaflet of the plasma membrane of the infected cells to form immature viral particles. Gag-membrane interactions are mediated through the myristoylated (Myr) N-terminal matrix (MA) domain of Gag, which eventually multimerize on the membrane to form trimers and higher order oligomers. The study of the structure and dynamics of peripheral membrane proteins like MA has been challenging for both experimental and computational studies due to the complex transient dynamics of protein-membrane interactions. Although the roles of anionic phospholipids (PIP2, PS) and the Myr group in the membrane targeting and stable membrane binding of MA are now well-established, the cooperative interactions between the MA monomers and MA-membrane remain elusive in the context of viral assembly and release. Our present study focuses on the membrane binding dynamics of a higher order oligomeric structure of MA protein (a dimer of trimers), which has not been explored before. Employing time-lagged independent component analysis (tICA) to our microsecond-long trajectories, we investigate conformational changes of the matrix protein induced by membrane binding. Interestingly, the Myr switch of an MA monomer correlates with the conformational switch of adjacent monomers in the same trimer. Together, our findings suggest complex protein dynamics during the formation of the immature HIV-1 lattice; while MA trimerization facilitates Myr insertion, MA trimer-trimer interactions in the immature lattice can hinder the same.

Molecular dynamics simulations of HIV-1 matrix-membrane interactions at different stages of viral maturation

Although the structural rearrangement of the membrane-bound matrix (MA) protein trimers upon HIV-1 maturation has been reported, the consequences of MA maturation on the MA-lipid interactions are not well understood. Long-timescale molecular dynamics simulations of the MA multimeric assemblies of immature and mature virus particles with our realistic asymmetric membrane model have explored MA-lipid interactions and lateral organization of lipids around MA complexes. The number of stable MA-phosphatidylserine and MA-phosphatidylinositol 4,5-bisphosphate (PIP2) interactions at the trimeric interface of the mature MA complex is observed to be greater compared to that of the immature MA complex. Our simulations identified an alternative PIP2-binding site in the immature MA complex where the multivalent headgroup of a PIP2 lipid with a greater negative charge binds to multiple basic amino acid residues such as ARG3 residues of both the MA monomers at the trimeric interface and highly basic region (HBR) residues (LYS29, LYS31) of one of the MA monomers. Our enhanced sampling simulations have explored the conformational space of phospholipids at different binding sites of the trimer-trimer interface of MA complexes that are not accessible by conventional unbiased molecular dynamics. Unlike the immature MA complex, the 2' acyl tail of two PIP2 lipids at the trimeric interface of the mature MA complex is observed to sample stable binding pockets of MA consisting of helix-4 residues. Together, our results provide molecular-level insights into the interactions of MA trimeric complexes with membrane and different lipid conformations at the specific binding sites of MA protein before and after viral maturation.

Conformational transitions of the HIV-1 Gag polyprotein upon multimerization and gRNA binding

During the HIV-1 assembly process, the Gag polyprotein multimerizes at the producer cell plasma membrane, resulting in the formation of spherical immature virus particles. Gag-genomic RNA (gRNA) interactions play a crucial role in the multimerization process, which is yet to be fully understood. We performed large-scale all-atom molecular dynamics simulations of membrane-bound full-length Gag dimer, hexamer, and 18-mer. The inter-domain dynamic correlation of Gag, quantified by the heterogeneous elastic network model applied to the simulated trajectories, is observed to be altered by implicit gRNA binding, as well as the multimerization state of the Gag. The lateral dynamics of our simulated membrane-bound Gag proteins, with and without gRNA binding, agree with prior experimental data and help to validate our simulation models and methods. The gRNA binding is observed to affect mainly the SP1 domain of the 18-mer and the matrix-capsid linker domain of the hexamer. In the absence of gRNA binding, the independent dynamical motion of the nucleocapsid domain results in a collapsed state of the dimeric Gag. Unlike stable SP1 helices in the six-helix bundle, without IP6 binding, the SP1 domain undergoes a spontaneous helix-to-coil transition in the dimeric Gag. Together, our findings reveal conformational switches of Gag at different stages of the multimerization process and predict that the gRNA binding reinforces an efficient binding surface of Gag for multimerization, and also regulates the dynamic organization of the local membrane region itself.

Cooperative Membrane Binding of HIV-1 Matrix Proteins

The HIV-1 assembly process begins with a newly synthesized Gag polyprotein being targeted to the inner leaflet of the plasma membrane of the infected cells to form immature viral particles. Gag-membrane interactions are mediated through the myristoylated(Myr) N-terminal matrix (MA) domain of Gag which eventually multimerize on the membrane to form trimers and higher-order oligomers. The study of the structure and dynamics of peripheral membrane proteins like MA has been challenging for both experimental and computational studies due to the complex dynamics of protein-membrane interactions. Although the roles of anionic phospholipids (PIP2, PS) and the Myr group in the membrane targeting and stable membrane binding of MA are now well-established, the cooperative interactions between MA monomers and MA-membrane still remain elusive. Our present study focuses on the membrane binding dynamics of a higher-order oligomeric structure of MA protein (a dimer of trimers), which has not been explored before. Employing time-lagged independent component analysis (tICA) to our microsecond-long trajectories, we investigate conformational changes of the matrix protein induced by membrane binding. Interestingly, the Myr switch of a MA monomer correlates with the conformational switch of adjacent monomers in the same trimer. Together, our findings suggest that MA trimerization facilitates Myr insertion, but MA trimer-trimer interactions in the lattice of immature HIV-1 particles can hinder the same. Additionally, local lipid density patterns of different lipid species provide a signature of the initial stage of lipid-domain formation upon membrane binding of the protein complex.

TOC

Atomic Mechanisms of Transthyretin Tetramer Dissociation Studied by Molecular Dynamics Simulations

The dissociation of the transthyretin (TTR) tetramer into a monomer is closely related to various TTR amyloidoses in humans. While the tetramer dissociation has been reported to be the rate-limiting step for TTR aggregation, few details are known about the mechanism. Here, molecular dynamics (MD) simulations were performed by combining conventional MD and biased metadynamics to investigate the mechanism for the wild-type (WT) and mutant (T119M) structures. Both were found to have a great deal in common. Conventional MD simulations reveal that interfacial hydrophobic interactions contribute significantly to stabilize the tetramer. Interfacial residues including L17, V20, L110, and V121 with close contacts form a hydrophobic channel. Metadynamics simulations indicate that the mouth opening of the hydrophobic channel is the first and the most difficult step for dissociation. Interactions of V20 between opposing dimers lock four monomers into the tetramer, and disruption of the interactions is found to be involved in the final step. During the dissociation, an increasing extent of solvation was observed by calculating the radial distribution functions of water around interfacial hydrophobic residues, suggesting that water plays a role in driving the tetramer dissociation. Moreover, compared to T119, residue M119 has a longer side chain that extends into the hydrophobic channel, making solvation more difficult, consistent with a higher energy barrier for dissociation of the T119M tetramer. This result provides a good explanation for the protective role of the T119M mutation. Overall, this study can provide atomic-level insights to better understand the pathogenesis of TTR amyloidosis and guide rational drug design in the future.

Metal-Organic Framework-Derived Carbon as a Photoacoustic Modulator of Alzheimer's β-Amyloid Aggregate Structure

Photoacoustic materials emit acoustic waves into the surrounding by absorbing photon energy. In an aqueous environment, light-induced acoustic waves form cavitation bubbles by altering the localized pressure to trigger the phase transition of liquid water into vapor. In this study, we report photoacoustic dissociation of beta-amyloid (Aβ) aggregates, a hallmark of Alzheimer's disease, by metal-organic framework-derived carbon (MOFC). MOFC exhibits a near-infrared (NIR) light-responsive photoacoustic characteristic that possesses defect-rich and entangled graphitic layers that generate intense cavitation bubbles by absorbing tissue-penetrable NIR light. According to our video analysis, the photoacoustic cavitation by MOFC occurs within milliseconds in the water, which was controllable by NIR light dose. The photoacoustic cavitation successfully transforms robust, β-sheet-dominant neurotoxic Aβ aggregates into nontoxic debris by changing the asymmetric distribution of water molecules around the Aβ's amino acid residues. This work unveils the therapeutic potential of NIR-triggered photoacoustic cavitation as a modulator of the Aβ aggregate structure.

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.

Tug-of-War between Internal and External Frictions and Viscosity Dependence of Rate in Biological Reactions

The role of water in biological processes is studied in three reactions, namely, the Fe-CO bond rupture in myoglobin, GB1 unfolding, and insulin dimer dissociation. We compute both internal and external components of friction on relevant reaction coordinates. In all of the three cases, the cross-correlation between forces from protein and water is found to be large and negative that serves to reduce the total friction significantly, increase the calculated reaction rate, and weaken solvent viscosity dependence. The computed force spectrum reveals bimodal 1/f noise, suggesting the use of a non-Markovian rate theory.

Protein-Protein Interaction-Gaussian Accelerated Molecular Dynamics (PPI-GaMD): Characterization of Protein Binding Thermodynamics and Kinetics

Protein-protein interactions (PPIs) play key roles in many fundamental biological processes such as cellular signaling and immune responses. However, it has proven challenging to simulate repetitive protein association and dissociation in order to calculate binding free energies and kinetics of PPIs due to long biological timescales and complex protein dynamics. To address this challenge, we have developed a new computational approach to all-atom simulations of PPIs based on a robust Gaussian accelerated molecular dynamics (GaMD) technique. The method, termed "PPI-GaMD", selectively boosts interaction potential energy between protein partners to facilitate their slow dissociation. Meanwhile, another boost potential is applied to the remaining potential energy of the entire system to effectively model the protein's flexibility and rebinding. PPI-GaMD has been demonstrated on a model system of the ribonuclease barnase interactions with its inhibitor barstar. Six independent 2 μs PPI-GaMD simulations have captured repetitive barstar dissociation and rebinding events, which enable calculations of the protein binding thermodynamics and kinetics simultaneously. The calculated binding free energies and kinetic rate constants agree well with the experimental data. Furthermore, PPI-GaMD simulations have provided mechanistic insights into barstar binding to barnase, which involves long-range electrostatic interactions and multiple binding pathways, being consistent with previous experimental and computational findings of this model system. In summary, PPI-GaMD provides a highly efficient and easy-to-use approach for binding free energy and kinetics calculations of PPIs.

Challenges and frontiers of computational modelling of biomolecular recognition

Biomolecular recognition including binding of small molecules, peptides and proteins to their target receptors plays a key role in cellular function and has been targeted for therapeutic drug design. However, the high flexibility of biomolecules and slow binding and dissociation processes have presented challenges for computational modelling. Here, we review the challenges and computational approaches developed to characterise biomolecular binding, including molecular docking, molecular dynamics simulations (especially enhanced sampling) and machine learning. Further improvements are still needed in order to accurately and efficiently characterise binding structures, mechanisms, thermodynamics and kinetics of biomolecules in the future.

The molecular mechanism behind the stabilization of insulin by choline and geranate (CAGE) ionic liquids - computational insights into oral insulin drug formulation

Insulin is a principal hormone that is involved in the regulation of glucose levels in the blood. Oral insulin formulation is a recent development in drug delivery systems. Biocompatible choline-based ionic liquids (ILs) show promising antibacterial activity and are useful for oral and transdermal drug delivery applications. Choline and geranate (CAGE) ILs enhance the stability and oral efficacy of insulin delivery. The molecular mechanism behind insulin formulation in the oral form is at issue. In the present work, the molecular-level understanding of CAGE ILs in insulin is scrutinized by employing atomistic molecular dynamics (MD) simulations. To identify the stability of insulin in an IL medium, we have studied a series of concentration (mole fraction 0.05-1.00) of CAGE ILs with an insulin dimer. It can be well evidenced from the experimental reports that in an aqueous medium, there is a refashioning of CAGE nanostructures at 0.50 mole fraction. It is found from our calculations that the first solvation shell of insulin is readily occupied by choline and geranate ions in the presence of water. Moreover, the geranate ions strongly interacted with the water molecules and thereby, eliminating the intermolecular hydrogen bonding (H-bonding) interactions towards the insulin at 0.30-0.50 mole fraction of CAGE ILs. The most desirable 0.30-0.50 mole fraction of CAGE invigorates water-mediated H-bonding interactions with geranate ions, which also enhances the electrostatic behavior around the vicinity of the insulin dimer. These important findings can help in the development of oral insulin drug delivery and related applications.

Structural Stability of Insulin Oligomers and Protein Association-Dissociation Processes: Free Energy Landscape and Universal Role of Water

Association and dissociation of proteins are important biochemical events. In this Feature Article, we analyze the available studies of these processes for insulin oligomers in aqueous solution. We focus on the solvation of the insulin monomer in water, stability and dissociation of its dimer, and structural integrity of the hexamer. The intricate role of water in solvation of the dimer- and hexamer-forming surfaces, in long-range interactions between the monomers and the stability of the oligomers, is discussed. Ten water molecules inside the central cavity stabilize the structure of the insulin hexamer. We discuss how different order parameters can be used to understand the dissociation of the insulin dimer. The calculation of the rate using a recently computed multidimensional free energy provides considerable insight into the interplay between protein and water dynamics.

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.

Computational IR Spectroscopy of Insulin Dimer Structure and Conformational Heterogeneity

We have investigated the structure and conformational dynamics of insulin dimer using a Markov state model (MSM) built from extensive unbiased atomistic molecular dynamics simulations and performed infrared spectral simulations of the insulin MSM to describe how structural variation within the dimer can be experimentally resolved. Our model reveals two significant conformations to the dimer: a dominant native state consistent with other experimental structures of the dimer and a twisted state with a structure that appears to reflect a ∼55° clockwise rotation of the native dimer interface. The twisted state primarily influences the contacts involving the C-terminus of insulin's B chain, shifting the registry of its intermolecular hydrogen bonds and reorganizing its side-chain packing. The MSM kinetics predict that these configurations exchange on a 14 μs time scale, largely passing through two Markov states with a solvated dimer interface. Computational amide I spectroscopy of site-specifically ¹³C¹⁸O labeled amides indicates that the native and twisted conformation can be distinguished through a series of single and dual labels involving the B24F, B25F, and B26Y residues. Additional structural heterogeneity and disorder is observed within the native and twisted states, and amide I spectroscopy can also be used to gain insight into this variation. This study will provide important interpretive tools for IR spectroscopic investigations of insulin structure and transient IR kinetics experiments studying the conformational dynamics of insulin dimer.

Evaluating the Efficiency of the Martini Force Field to Study Protein Dimerization in Aqueous and Membrane Environments

Protein-protein complex assembly is one of the major drivers of biological response. Understanding the mechanisms of protein oligomerization/dimerization would allow one to elucidate how these complexes participate in biological activities and could ultimately lead to new approaches in designing novel therapeutic agents. However, determining the exact association pathways and structures of such complexes remains a challenge. Here, we use parallel tempering metadynamics simulations in the well-tempered ensemble to evaluate the performance of Martini 2.2P and Martini open-beta 3 (Martini 3) force fields in reproducing the structure and energetics of the dimerization process of membrane proteins and proteins in an aqueous solution in reasonable accuracy and throughput. We find that Martini 2.2P systematically overestimates the free energy of association by estimating large barriers in distinct areas, which likely leads to overaggregation when multiple monomers are present. In comparison, the less viscous Martini 3 results in a systematic underestimation of the free energy of association for proteins in solution, while it performs well in describing the association of membrane proteins. In all cases, the near-native dimer complexes are identified as minima in the free energy surface albeit not always as the lowest minima. In the case of Martini 3, we find that the spurious supramolecular protein aggregation present in Martini 2.2P multimer simulations is alleviated and thus this force field may be more suitable for the study of protein oligomerization. We propose that the use of enhanced sampling simulations with a refined coarse-grained force field and appropriately defined collective variables is a robust approach for studying the protein dimerization process, although one should be cautious of the ranking of energy minima.

Dynamics and Infrared Spectrocopy of Monomeric and Dimeric Wild Type and Mutant Insulin

The infrared spectroscopy and dynamics of -CO labels in wild type and mutant insulin monomer and dimer are characterized from molecular dynamics simulations using validated force fields. It is found that the spectroscopy of monomeric and dimeric forms in the region of the amide-I vibration differs for residues B24-B26 and D24-D26, which are involved in dimerization of the hormone. Also, the spectroscopic signatures change for mutations at position B24 from phenylalanine, which is conserved in many organisms and is known to play a central role in insulin aggregation, to alanine or glycine. Using three different methods to determine the frequency trajectories (solving the nuclear Schrödinger equation on an effective 1-dimensional potential energy curve, using instantaneous normal modes, and using parametrized frequency maps) leads to the same overall conclusions. The spectroscopic response of monomeric WT and mutant insulin differs from that of their respective dimers, and the spectroscopy of the two monomers in the dimer is also not identical. For the WT and F24A and F24G monomers, spectroscopic shifts are found to be ∼20 cm^-1 for residues (B24-B26) located at the dimerization interface. Although the crystal structure of the dimer is that of a symmetric homodimer, dynamically the two monomers are not equivalent on the nanosecond time scale. Together with earlier work on the thermodynamic stability of the WT and the same mutants, it is concluded that combining computational and experimental infrared spectroscopy provides a potentially powerful way to characterize the aggregation state and dimerization energy of modified insulins.

MEPSAnd: minimum energy path surface analysis over n-dimensional surfaces

SUMMARY

n-dimensional energy surfaces are becoming computationally accessible, yet interpreting their information is not straightforward. We present minimum energy path surface analysis over n-dimensional surfaces (MEPSAnd), an open source GUI-based program that natively calculates minimum energy paths across energy surfaces of any number of dimensions. Among other features, MEPSAnd can compute the path through lowest barriers and automatically provide a set of alternative paths. MEPSAnd offers distinct plotting solutions as well as direct python scripting.

AVAILABILITY AND IMPLEMENTATION

MEPSAnd is freely available (under GPLv3 license) at: http://bioweb.cbm.uam.es/software/MEPSAnd/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Insulin Dissociates by Diverse Mechanisms of Coupled Unfolding and Unbinding

The protein hormone insulin exists in various oligomeric forms, and a key step in binding its cellular receptor is dissociation of the dimer. This dissociation process and its corresponding association process have come to serve as paradigms of coupled (un)folding and (un)binding more generally. Despite its fundamental and practical importance, the mechanism of insulin dimer dissociation remains poorly understood. Here, we use molecular dynamics simulations, leveraging recent developments in umbrella sampling, to characterize the energetic and structural features of dissociation in unprecedented detail. We find that the dissociation is inherently multipathway with limiting behaviors corresponding to conformational selection and induced fit, the two prototypical mechanisms of coupled folding and binding. Along one limiting path, the dissociation leads to detachment of the C-terminal segment of the insulin B chain from the protein core, a feature believed to be essential for receptor binding. We simulate IR spectroscopy experiments to aid in interpreting current experiments and identify sites where isotopic labeling can be most effective for distinguishing the contributions of the limiting mechanisms.

Dynamical control by water at a molecular level in protein dimer association and dissociation

Water, often termed as the "lubricant of life," is expected to play an active role in navigating protein dissociation-association reactions. In order to unearth the molecular details, we first compute the free-energy surface (FES) of insulin dimer dissociation employing metadynamics simulation, and then carry out analyses of insulin dimerization and dissociation using atomistic molecular-dynamics simulation in explicit water. We select two sets of initial configurations from 1) the dissociated state and 2) the transition state, and follow time evolution using several long trajectories (∼1-2 μs). During the process we not only monitor configuration of protein monomers, but also the properties of water. Although the equilibrium structural properties of water between the two monomers approach bulklike characteristics at a separation distance of ∼5 nm, the dynamics differ considerably. The complex association process is observed to be accompanied by several structural and dynamical changes of the system, such as large-scale correlated water density fluctuations, coupled conformational fluctuation of protein monomers, a dewettinglike transition with the change of intermonomeric distance R_MM from ∼4 to ∼2 nm, orientation of monomers and hydrophobic hydration in the monomers. A quasistable, solvent-shared, protein monomer pair (SSPMP) forms at around 2 nm during association process which is a local free-energy minimum having ∼50-60% of native contacts. Simulations starting with arrangements sampled from the transition state (TS) of the dimer dissociation reveal that the final outcome depends on relative orientation of the backbone in the "hotspot" region.

Calculating the absolute binding free energy of the insulin dimer in an explicit solvent

Insulin is a significant hormone in the regulation of glucose level in the blood. Its monomers bind to each other to form dimers or hexamers through a complex process. To study the binding of the insulin dimer, we first calculate its absolute binding free energy by the steered molecular dynamics method and the confinement method based on a fictitious thermodynamic cycle. After considering some special correction terms, the final calculated binding free energy at 298 K is −8.97 ± 1.41 kcal mol⁻¹, which is close to the experimental value of −7.2 ± 0.8 kcal mol⁻¹. Furthermore, we discuss the important residue–residue interactions between the insulin monomers, including hydrophobic interactions, π–π interactions and hydrogen bond interactions. The analysis reveals five key residues, Vla^B12, Tyr^B16, Phe^B24, Phe^B25, and Tyr^B26, for the dimerization of the insulin. We also perform MM-PBSA calculations for the wild-type dimer and some mutants and study the roles of the key residues by the change of the binding energy of the insulin dimer.

In this paper, we calculate the absolute binding free energy of an insulin dimer by steered MD method. The result of −8.97 kcal mol⁻¹ is close to the experimental value −7.2 kcal mol⁻¹. We also analyze the residue–residue interactions.

Effect of ethanol on insulin dimer dissociation

Insulin-dimer dissociation is an essential biochemical process required for the activity of the hormone. We investigate this dissociation process at the molecular level in water and at the same time, in 5% and 10% water-ethanol mixtures. We compute the free energy surface of the protein dissociation processes by employing biased molecular dynamics simulation. In the presence of ethanol (EtOH), we observe a marked lowering in the free energy barrier of activation of dimer dissociation from that in the neat water, by as much as ∼50%, even in the 5% water-ethanol solution. In addition, ethanol is found to induce significant changes in the dissociation pathway. We extract the most probable conformations of the intermediate states along the minimum energy pathway in the case of all the three concentrations (EtOH mole fractions 0, 5, and 10). We explore the change in microscopic structures that occur in the presence of ethanol. Interestingly, we discover a stable intermediate state in the water-ethanol binary mixture where the centers of the monomers are separated by about 3 nm and the contact order parameter is close to zero. This intermediate is stabilized by the wetting of the interface between the two monomers by the preferential distribution of ethanol and water molecules. This wetting serves to reduce the free energy barrier significantly and thus results in an increase in the rate of dimer dissociation. We also analyze the solvation of the two monomers during the dissociation and both the proteins' departure from the native state configuration to obtain valuable insights into the dimer dissociation processes.