Effect of ethanol on insulin dimer dissociation

Martini 3 OliGo̅mers: A Scalable Approach for Multimers and Fibrils in GROMACS

Martini 3 is a widely used coarse-grained simulation method for large-scale biomolecular simulations. It can be combined with a Go̅ model to realistically describe higher-order protein structures while allowing the folding and unfolding events. However, as of today, this method has largely been used only for individual monomers. In this article, we describe how the Go̅ model can be implemented within the framework of Martini 3 for a multimer system, taking into account both intramolecular and intermolecular interactions in an oligomeric protein system. We demonstrate the method by showing how it can be applied to both structural stability maintenance and assembly/disassembly of protein oligomers, using aquaporin tetramer, insulin dimer, and amyloid-β fibril as examples. We find that addition of intermolecular Go̅ potentials stabilizes the quaternary structure of proteins. The strength of the Go̅ potentials can be tuned so that the internal fluctuations of proteins match the behavior of atomistic simulation models, however, the results also show that the use of too strong intermolecular Go̅ potentials weakens the chemical specificity of oligomerization. The Martini-Go̅ model presented here enables the use of Go̅ potentials in oligomeric molecular systems in a computationally efficient and parallelizable manner, especially in the case of homopolymers, where the number of identical protein monomers is high. This paves the way for coarse-grained simulations of large protein complexes, such as viral protein capsids and prion fibrils, in complex biological environments.

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.

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

New Insights into the Cooperativity and Dynamics of Dimeric Enzymes

A survey of protein databases indicates that the majority of enzymes exist in oligomeric forms, with about half of those found in the UniProt database being homodimeric. Understanding why many enzymes are in their dimeric form is imperative. Recent developments in experimental and computational techniques have allowed for a deeper comprehension of the cooperative interactions between the subunits of dimeric enzymes. This review aims to succinctly summarize these recent advancements by providing an overview of experimental and theoretical methods, as well as an understanding of cooperativity in substrate binding and the molecular mechanisms of cooperative catalysis within homodimeric enzymes. Focus is set upon the beneficial effects of dimerization and cooperative catalysis. These advancements not only provide essential case studies and theoretical support for comprehending dimeric enzyme catalysis but also serve as a foundation for designing highly efficient catalysts, such as dimeric organic catalysts. Moreover, these developments have significant implications for drug design, as exemplified by Paxlovid, which was designed for the homodimeric main protease of SARS-CoV-2.

Molecular Mechanisms of Inhibition of Protein Amyloid Fibril Formation: Evidence and Perspectives Based on Kinetic Models

Inhibition of fibril formation is considered a possible treatment strategy for amyloid-related diseases. Understanding the molecular nature of inhibitor action is crucial for the design of drug candidates. In the present review, we describe the common kinetic models of fibril formation and classify known inhibitors by the mechanism of their interactions with the aggregating protein and its oligomers. This mechanism determines the step or steps of the aggregation process that become inhibited and the observed changes in kinetics and equilibrium of fibril formation. The results of numerous studies indicate that possible approaches to antiamyloid inhibitor discovery include the search for the strong binders of protein monomers, cappers blocking the ends of the growing fibril, or the species absorbing on the surface of oligomers preventing nucleation. Strongly binding inhibitors stabilizing the native state can be promising for the structured proteins while designing the drug candidates targeting disordered proteins is challenging.

Formulation Excipients and Their Role in Insulin Stability and Association State in Formulation

While  excipients are often overlooked as the "inactive" ingredients in pharmaceutical formulations, they often play a critical role in protein stability and absorption kinetics. Recent work has identified an ultrafast absorbing insulin formulation that is the result of excipient modifications. Specifically, the insulin monomer can be isolated by replacing zinc and the phenolic preservative metacresol with phenoxyethanol as an antimicrobial agent and an amphiphilic acrylamide copolymer excipient for stability. A greater understanding is needed of the interplay between excipients, insulin association state, and stability in order to optimize this formulation. Here, we formulated insulin with different preservatives and stabilizing excipient concentrations using both insulin lispro and regular human insulin and assessed the insulin association states using analytical ultracentrifugation as well as formulation stability. We determined that phenoxyethanol is required to eliminate hexamers and promote a high monomer content even in a zinc-free lispro formulation. There is also a concentration dependent relationship between the concentration of polyacrylamide-based copolymer excipient and insulin stability, where a concentration greater than 0.1 g/mL copolymer is required for a mostly monomeric zinc-free lispro formulation to achieve stability exceeding that of Humalog in a stressed aging assay. Further, we determined that under the formulation conditions tested zinc-free regular human insulin remains primarily hexameric and is not at this time a promising candidate for rapid-acting formulations.

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.

Coacervation of poly-electrolytes in the presence of lipid bilayers: mutual alteration of structure and morphology

Many intrinsically disordered peptides have been shown to undergo liquid-liquid phase separation and form complex coacervates, which play various regulatory roles in the cell. Recent experimental studies found that such phase separation processes may also occur at the lipid membrane surface and help organize biomolecules during signaling events; in some cases, phase separation of proteins at the membrane surface was also observed to lead to significant remodeling of the membrane morphology. The molecular mechanisms that govern the interactions between complex coacervates and lipid membranes and the impacts of such interactions on their structure and morphology, however, remain unclear. Here we study the coacervation of poly-glutamate (E30) and poly-lysine (K30) in the presence of lipid bilayers of different compositions. We carry out explicit-solvent coarse-grained molecular dynamics simulations by using the MARTINI (v3.0) force-field. We find that more than 20% anionic lipids are required for the coacervate to form stable contact with the bilayer. Upon wetting, the coacervate induces negative curvature to the bilayer and facilitates local lipid demixing, without any peptide insertion. The magnitude of negative curvature, extent of lipid demixing, and asphericity of the coacervate increase with the concentration of anionic lipids. Overall, we observe a decrease in the number of contacts among the polyelectrolytes as the droplet spreads over the bilayer. Therefore, unlike previous suggestions, interactions among polyelectrolytes do not constitute a driving force for the membrane bending upon wetting by the coacervate. Rather, analysis of interaction energy components suggests that bending of the membrane is favored by enhanced interactions between polyelectrolytes with lipids as well as with counterions. Kinetic studies reveal that, at the studied polyelectrolyte concentrations, the coacervate formation precedes bilayer wetting.

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.

Quantitative molecular simulations

All-atom simulations can provide molecular-level insights into the dynamics of gas-phase, condensed-phase and surface processes. One important requirement is a sufficiently realistic and detailed description of the underlying intermolecular interactions. The present perspective provides an overview of the present status of quantitative atomistic simulations from colleagues' and our own efforts for gas- and solution-phase processes and for the dynamics on surfaces. Particular attention is paid to direct comparison with experiment. An outlook discusses present challenges and future extensions to bring such dynamics simulations even closer to reality.

Atomistic Simulations for Reactions and Vibrational Spectroscopy in the Era of Machine Learning─Quo Vadis?

Atomistic simulations using accurate energy functions can provide molecular-level insight into functional motions of molecules in the gas and in the condensed phase. This Perspective delineates the present status of the field from the efforts of others and some of our own work and discusses open questions and future prospects. The combination of physics-based long-range representations using multipolar charge distributions and kernel representations for the bonded interactions is shown to provide realistic models for the exploration of the infrared spectroscopy of molecules in solution. For reactions, empirical models connecting dedicated energy functions for the reactant and product states allow statistically meaningful sampling of conformational space whereas machine-learned energy functions are superior in accuracy. The future combination of physics-based models with machine-learning techniques and integration into all-purpose molecular simulation software provides a unique opportunity to bring such dynamics simulations closer to reality.

Influence of Aqueous Arginine Solution on Regulating Conformational Stability and Hydration Properties of the Secondary Structural Segments of a Protein at Elevated Temperatures: A Molecular Dynamics Study

The effects of aqueous arginine solution on the conformational stability of the secondary structural segments of a globular protein, ubiquitin, and the structure and dynamics of the surrounding water and arginine were examined by performing atomistic molecular dynamics (MD) simulations. Attempts have been made to identify the osmolytic efficacy of arginine solution, and its influence in guiding the hydration properties of the protein at an elevated temperature of 450 K. The similar properties of the protein in pure water at elevated temperatures were computed and compared. Replica exchange MD simulation was performed to explore the arginine solution's sensitivity in stabilizing the protein conformations for a wide range of temperatures (300-450 K). It was observed that although all the helices and strands of the protein undergo unfolding at elevated temperature in pure water, they exhibited native-like conformational dynamics in the presence of arginine at both ambient and elevated temperatures. We find that the higher free energy barrier between the folded native and unfolded states of the protein primarily arises from the structural transformation of α-helix, relative to the strands. Our study revealed that the water structure around the secondary segments depends on the nature of amino acid compositions of the helices and strands. The reorientation of water dipoles around the helices and strands was found hindered due to the presence of arginine in the solution; such hindrance reduces the possibility of exchange of hydrogen bonds that formed between the secondary segments of protein and water (PW), and as a result, PW hydrogen bonds take longer time to relax than in pure water. On the other hand, the origin of slow relaxation of protein-arginine (PA) hydrogen bonds was identified to be due to the presence of different types of protein-bound arginine molecules, where arginine interacts with the secondary structural segments of the protein through multiple/bifurcated hydrogen bonds. These protein-bound arginine formed different kinds of bridged PA hydrogen bonds between amino acid residues of the same secondary segments or among multiple bonds and helped protein to conserve its native folded form firmly.

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.

Thermodynamic Origin of Differential Excipient-Lysozyme Interactions

Understanding the intricate interplay of interactions between proteins, excipients, ions and water is important to achieve the effective purification and stable formulation of protein therapeutics. The free energy of lysozyme interacting with two kinds of polyanionic excipients, citrate and tripolyphosphate, together with sodium chloride and TRIS-buffer, are analysed in multiple-walker metadynamics simulations to understand why tripolyphosphate causes lysozyme to precipitate but citrate does not. The resulting multiscale decomposition of energy and entropy components for water, sodium chloride, excipients and lysozyme reveals that lysozyme is more stabilised by the interaction of tripolyphosphate with basic residues. This is accompanied by more sodium ions being released into solution from tripolyphosphate than for citrate, whilst the latter instead has more water molecules released into solution. Even though lysozyme aggregation is not directly probed in this study, these different mechanisms are suspected to drive the cross-linking between lysozyme molecules with vacant basic residues, ultimately leading to precipitation.

Unfolding bovine α-lactalbumin with T-jump: Characterizing disordered intermediates via time-resolved x-ray solution scattering and molecular dynamics simulations

The protein folding process often proceeds through partially folded transient states. Therefore, a structural understanding of these disordered states is crucial for developing mechanistic models of the folding process. Characterization of unfolded states remains challenging due to their disordered nature, and incorporating multiple methods is necessary. Combining the time-resolved x-ray solution scattering (TRXSS) signal with molecular dynamics (MD), we are able to characterize transient partially folded states of bovine α-lactalbumin, a model system widely used for investigation of molten globule states, during its unfolding triggered by a temperature jump. We track the unfolding process between 20 µs and 70 ms and demonstrate that it passes through three distinct kinetic states. The scattering signals associated with these transient species are then analyzed with TRXSS constrained MD simulations to produce protein structures that are compatible with the input signals. Without utilizing any experimentally extracted kinetic information, the constrained MD simulation successfully drove the protein to an intermediate molten globule state; signals for two later disordered states are refined to terminal unfolded states. From our examination of the structural characteristics of these disordered states, we discuss the implications disordered states have on the folding process, especially on the folding pathway. Finally, we discuss the potential applications and limitations of this method.

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.

Revealing the Dynamical Role of Co-solvents in the Coupled Folding and Dimerization of Insulin

Solvent-protein interactions are important for protein biological functions, especially for a coupled folding and binding system such as insulin. By monitoring the change in the conformation of insulin dimers during dissociation with temperature-jump infrared spectroscopy, we show that co-solvents can significantly destabilize the dimers by perturbing their hydrophobic center. The transition from the native to intermediate dimer state is observed as the buried residues are exposed to solvents in the presence of 10% dimethyl sulfoxide and with α-helices unfolding when ethanol is present, which reduces the dissociation time dramatically to 50% and 20% of the value in a D₂O solution, respectively. We propose a self-consistent analysis using complementary methods to resolve this coupled folding and binding process and obtain a much higher rate of monomer association than of intermediate folding. Our results demonstrate that the conformational changes are critical in dimer formation and strongly affected by co-solvents.

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.

Stable Monomeric Insulin Formulations Enabled by Supramolecular PEGylation of Insulin Analogues

Current "fast-acting" insulin analogues contain amino acid modifications meant to inhibit dimer formation and shift the equilibrium of association states toward the monomeric state. However, the insulin monomer is highly unstable and current formulation techniques require insulin to primarily exist as hexamers to prevent aggregation into inactive and immunogenic amyloids. Insulin formulation excipients have thus been traditionally selected to promote insulin association into the hexameric form to enhance formulation stability. This study exploits a novel excipient for the supramolecular PEGylation of insulin analogues, including aspart and lispro, to enhance the stability and maximize the prevalence of insulin monomers in formulation. Using multiple techniques, it is demonstrated that judicious choice of formulation excipients (tonicity agents and parenteral preservatives) enables insulin analogue formulations with 70-80% monomer and supramolecular PEGylation imbued stability under stressed aging for over 100 h without altering the insulin association state. Comparatively, commercial "fast-acting" formulations contain less than 1% monomer and remain stable for only 10 h under the same stressed aging conditions. This simple and effective formulation approach shows promise for next-generation ultrafast insulin formulations with a short duration of action that can reduce the risk of post-prandial hypoglycemia in the treatment of diabetes.