Microcanonical versus canonical analysis of protein folding

Microcanonical insights into the physicochemical stability of the coformulation of insulin with amylin analogues

Injections of insulin are the main treatment for diabetes, but in the long run this therapy can induce serious drawbacks. This has inspired new drugs able to decrease insulin requirements. For instance, human amylin (hIAPP) is a small hormone cosecreted by pancreatic β-cells with insulin to which is a synergistic partner. However, the high amyloidogenicity of hIAPP precluded it as a therapeutics and led to the design of pramlintide (sIAPP), a chimeric analogue with substitutions (A25P, S28P, and S29P) inherited from the aggregation-resistant rat isoform (rIAPP). Despite sIAPP advantages, it still shares with hIAPP a poorly soluble profile at physiological pH that hampers its mixture with insulin. Recent improvements, as charge-enhanced mutants, have been proposed. For instance, sIAPP⁺ was screened in silico by purely microcanonical thermostatistical methods and adds to sIAPP an S20R mutation to uplift its solubility. This suggests that such physically inspired computational approach may also be auspicious on devising effective coformulations of insulin with amylin analogues. In this seminal attempt, we make comparative multicanonical simulations of regular acting human insulin coformulated with hIAPP, sIAPP, or sIAPP⁺. To assess the respective physicochemical stabilities against aggregation, we characterize the structural-phase transitions through the microcanonical thermodynamic formalism and evaluate their time lags using the classical nucleation theory. These results are then correlated with estimates of solvation free energies, modeled by the Poisson-Boltzmann equation, and structural propensities. Experimental essays are compared to our simulations and support our methodology.

Potential energy landscape and thermodynamic transitions of coarse-grained protein models revealed by the multicanonical generalized hybrid Monte Carlo method

In the present study, thermodynamic properties of coarse-grained protein models have been studied by an extended ensemble method. Two types of protein model were analyzed; one is categorized into a fast folder and the other into a slow folder. Both models exhibit the following thermodynamic transitions: the collapse transition between random coil states and spatially compact, but non-native states and the folding transition between the collapsed states and the folded native states. Caloric curve for the fast folder shows strong statistical ensemble dependence, while almost no ensemble dependence is found for the slow folder. Microcanonical caloric curve for the fast folder exhibits S-shaped temperature dependence on the internal energy around the collapse transition which is reminiscent of the van der Waals loop observed for the first order transition; at the transition temperature, the collapsed and random coil states coexist dynamically. The corresponding microcanonical heat capacity is found to have negative region around the transition. This kind of exotic behaviors could be utilized to distinguish fast folding proteins.

Diagrams of States of Single Flexible-Semiflexible Multi-Block Copolymer Chains: A Flat-Histogram Monte Carlo Study

The combination of flexibility and semiflexibility in a single molecule is a powerful design principle both in nature and in materials science. We present results on the conformational behavior of a single multiblock-copolymer chain, consisting of equal amounts of Flexible (F) and Semiflexible (S) blocks with different affinity to an implicit solvent. We consider a manifold of macrostates defined by two terms in the total energy: intermonomer interaction energy and stiffness energy. To obtain diagrams of states (pseudo-phase diagrams), we performed flat-histogram Monte Carlo simulations using the Stochastic Approximation Monte Carlo algorithm (SAMC). We have accumulated two-Dimensional Density of States (2D DoS) functions (defined on the 2D manifold of macrostates) for a SF-multiblock-copolymer chain of length N=64 with block lengths b = 4, 8, 16, and 32 in two different selective solvents. In an analysis of the canonical ensemble, we calculated the heat capacity and determined its maxima and the most probable morphologies in different regions of the state diagrams. These are rich in various, non-trivial morphologies, which are formed without any specific interactions, and depend on the block length and the type of solvent selectivity (preferring S or F blocks, respectively). We compared the diagrams with those for the non-selective solvent and reveal essential changes in some cases. Additionally, we implemented microcanonical analysis in the “conformational” microcanonical (NVU, where U is the potential energy) and the true microcanonical (NVE, where E is the total energy) ensembles with the aim to reveal and classify pseudo-phase transitions, occurring under the change of temperature.

Be positive: optimizing pramlintide from microcanonical analysis of amylin isoforms

Amylin, or human islet amyloid polypeptide (hIAPP), is a 37-residue hormone synergistic to insulin and co-secreted with it by β-cells in the pancreas. The deposition of its cytotoxic amyloid fibrils is strongly related to the progression of Type II diabetes (T2D) and islet graft failures. Notably, isoforms from some mammalian species, such as rats (rIAPP) and porcine (pIAPP), present a few key mutations preventing aggregation. This has lead to biotechnological development of drugs for adjunct therapies of T2D, such as pramlintide, a variant of hIAPP inspired by rIAPP whose proline substitutions have β-strand fibril-breaking properties. Ideally, such a drug should be formulated with insulin and co-administered, but this has been prevented by a poor solubility profile at the appropriate pH. Hopefully, this could be improved with appropriate point mutations, increasing the molecular net charge. Despite experimental progress, preliminary screening during rational drug design can greatly benefit from thermodynamic insight derived from molecular simulations. So we introduce microcanonical thermostatistics analysis of multicanonical (MUCA) simulations of wild-type amylin isoforms as a systematic assessment of protein thermostability. As a consequence of this comprehensive investigation, the most suitable single-point mutations able to optimize pramlintide are located among the wild-type amylin isoforms. In particular, we find that aggregation inhibition and increased solubility are inherited by pramlintide through further S20R substitution typical of pIAPP. Thus, we provide a consistent thermostatistical methodology to aid the design of improved adjunct therapies for T2D according to current clinical knowledge.

Microcanonical analysis of a finite-size nonequilibrium system

Microcanonical analysis is a powerful method that can be used to generalize the concept of phase transitions to finite-size systems. However, microcanonical analysis has only been applied to equilibrium systems. I show that it is possible to conduct the microcanonical analysis of a finite-size nonequilibrium system by generalizing the concept of microcanonical entropy. A one-dimensional asymmetric diffusion process is studied as an example for which such a generalized entropy can be explicitly found, and the microcanonical method is used to define a generalized phase transition for the finite-size nonequilibrium system.

Breakout character of islet amyloid polypeptide hydrophobic mutations at the onset of type-2 diabetes

Toxic fibrillar aggregates of islet amyloid polypeptide (IAPP) appear as the physical outcome of a peptidic phase transition signaling the onset of type-2 diabetes mellitus in different mammalian species. In particular, experimentally verified mutations on the amyloidogenic segment 20-29 in humans, cats, and rats are highly correlated with the molecular aggregation propensities. Through a microcanonical analysis of the aggregation of IAPP_{20-29} isoforms, we show that a minimalist one-bead hydrophobic-polar continuum model for protein interactions properly quantifies those propensities from free-energy barriers. Our results highlight the central role of sequence-dependent hydrophobic mutations on hot spots for stabilization, and thus for the engineering, of such biological peptides.

Thermodynamic and conformational insights into the phase transition of a single flexible homopolymer chain using replica exchange Monte Carlo method

The phase transition of a single flexible homopolymer chain in the limit condition of dilute solution is systematically investigated using a coarse-grained model. Replica exchange Monte Carlo method is used to enhance the performance of the conformation space exploration, and thus detailed investigation of phase behavior of the system can be provided. With the designed potentials, the coil-globule transition and the liquid-solid-like transition are identified, and the transition temperatures are measured with the conformational and thermodynamic analyses. Additionally, by extrapolating the coil-globule transition temperature, T Θ , and the liquid-solid-like transition temperature T(L → S) to the thermodynamic limit, N → ∞, we found no "tri-critical" point in the current model.

First-order coil-globule transition driven by vibrational entropy

By shifting the balance between conformational entropy and internal energy, polymers modify their shape under external stimuli, such as changes in temperature. Prominent among such transformations is the coil-globule transition, whereby a polymer can switch from an entropy-dominated coil conformation to a globular one, governed by energy. The nature of the coil-globule transition has remained elusive, with evidence for both continuous and discontinuous transitions, with the two-state behaviour of proteins as an instance of the latter. Theoretical models mostly predict second-order transitions. Here we introduce a model that takes into consideration hitherto neglected features common to any polymer. We show that a first-order phase transition smoothly appears as a function of the model parameters. Our results can relieve part of the conflicts between theory and experiments in the field of protein folding, in the wake of recent studies tracing back the remarkable properties of proteins to basic polymer physics.

Selecting fast-folding proteins by their rate of convergence

We propose a general method for predicting potentially good folders from a given number of amino acid sequences. Our approach is based on the calculation of the rate of convergence of each amino acid chain towards the native structure using only the very initial parts of the dynamical trajectories. It does not require any preliminary knowledge of the native state and can be applied to different kinds of models, including atomistic descriptions. We tested the method within both the lattice and off-lattice model frameworks and obtained several so far unknown good folders. The unbiased algorithm also allows one to determine the optimal folding temperature and takes at least 3-4 orders of magnitude fewer time steps than those needed to compute folding times.

Calculation of microcanonical entropy differences from configurational averages

A simple expression is derived, enabling the calculation of the entropy difference between two microcanonical equilibrium states at different energies in atomistic computer simulations. This expression only requires potential energy samples from molecular dynamics or Monte Carlo simulations at the relevant energies. This presents an alternative to switching methods such as thermodynamic integration or nonequilibrium work relations, as well as flat-histogram random walks, all of which involve sampling in between the relevant states. The method is especially suited for small (nanoscopic) systems such as clusters and proteins, and is applicable to first-principles data directly.

Microcanonical entropy inflection points: key to systematic understanding of transitions in finite systems

We introduce a systematic classification method for the analogs of phase transitions in finite systems. This completely general analysis, which is applicable to any physical system and extends toward the thermodynamic limit, is based on the microcanonical entropy and its energetic derivative, the inverse caloric temperature. Inflection points of this quantity signal cooperative activity and thus serve as distinct indicators of transitions. We demonstrate the power of this method through application to the long-standing problem of liquid-solid transitions in elastic, flexible homopolymers.

Structural basis of folding cooperativity in model proteins: insights from a microcanonical perspective

Two-state cooperativity is an important characteristic in protein folding. It is defined by a depletion of states that lie energetically between folded and unfolded conformations. There are different ways to test for two-state cooperativity; however, most of these approaches probe indirect proxies of this depletion. Generalized-ensemble computer simulations allow us to unambiguously identify this transition by a microcanonical analysis on the basis of the density of states. Here, we present a detailed characterization of several helical peptides obtained by coarse-grained simulations. The level of resolution of the coarse-grained model allowed to study realistic structures ranging from small α-helices to a de novo three-helix bundle without biasing the force field toward the native state of the protein. By linking thermodynamic and structural features, we are able to show that whereas short α-helices exhibit two-state cooperativity, the type of transition changes for longer chain lengths because the chain forms multiple helix nucleation sites, stabilizing a significant population of intermediate states. The helix bundle exhibits signs of two-state cooperativity owing to favorable helix-helix interactions, as predicted from theoretical models. A detailed analysis of secondary and tertiary structure formation fits well into the framework of several folding mechanisms and confirms features that up to now have been observed only in lattice models.

Generalized simulated tempering for exploring strong phase transitions

An extension of the simulation tempering algorithm is proposed. It is shown to be particularly suited to the exploration of first-order phase transition systems characterized by the backbending or S-loop in the statistical temperature or a microcanonical caloric curve. A guided Markov process in an auxiliary parameter space systematically combines a set of parametrized Tsallis-weight ensemble simulations, which are targeted to transform unstable or metastable energy states of canonical ensembles into stable ones and smoothly join ordered and disordered phases across phase transition regions via a succession of unimodal energy distributions. The inverse mapping between the sampling weight and the effective temperature enables an optimal selection of relevant Tsallis-weight parameters. A semianalytic expression for the biasing weight in parameter space is adaptively updated "on the fly" during the simulation to achieve rapid convergence. Accelerated tunneling transitions with a comprehensive sampling for phase-coexistent states are explicitly demonstrated in systems subject to strong hysteresis including Potts and Ising spin models and a 147 atom Lennard-Jones cluster.

Generalized replica exchange method

We present a powerful replica exchange method, particularly suited to first-order phase transitions associated with the backbending in the statistical temperature, by merging an optimally designed generalized ensemble sampling with replica exchanges. The key ingredients of our method are parametrized effective sampling weights, smoothly joining ordered and disordered phases with a succession of unimodal energy distributions by transforming unstable or metastable energy states of canonical ensembles into stable ones. The inverse mapping between the sampling weight and the effective temperature provides a systematic way to design the effective sampling weights and determine a dynamic range of relevant parameters. Illustrative simulations on Potts spins with varying system size and simulation conditions demonstrate a comprehensive sampling for phase-coexistent states with a dramatic acceleration of tunneling transitions. A significant improvement over the power-law slowing down of mean tunneling times with increasing system size is obtained, and the underlying mechanism for accelerated tunneling is discussed.

Phase transitions of a single polymer chain: A Wang-Landau simulation study

A single flexible homopolymer chain can assume a variety of conformations which can be broadly classified as expanded coil, collapsed globule, and compact crystallite. Here we study transitions between these conformational states for an interaction-site polymer chain comprised of N=128 square-well-sphere monomers with hard-sphere diameter sigma and square-well diameter lambdasigma. Wang-Landau sampling with bond-rebridging Monte Carlo moves is used to compute the density of states for this chain and both canonical and microcanonical analyses are used to identify and characterize phase transitions in this finite size system. The temperature-interaction range (i.e., T-lambda) phase diagram is constructed for lambda<or=1.30. Chains assume an expanded coil conformation at high temperatures and a crystallite structure at low temperatures. For lambda>1.06 these two states are separated by an intervening collapsed globule phase and thus, with decreasing temperature a chain undergoes a continuous coil-globule (collapse) transition followed by a discontinuous globule-crystal (freezing) transition. For well diameters lambda<1.06 the collapse transition is pre-empted by the freezing transition and thus there is a direct first-order coil-crystal phase transition. These results confirm the recent prediction, based on a lattice polymer model, that a collapsed globule state is unstable with respect to a solid phase for flexible polymers with sufficiently short-range monomer-monomer interactions.

All-or-none proteinlike folding transition of a flexible homopolymer chain

Here we report a first-order all-or-none transition from an expanded coil to a compact crystallite for a flexible polymer chain. Wang-Landau sampling is used to construct the complete density of states for square-well chains up to length 256. Analysis within both the microcanonical and canonical ensembles shows a direct freezing transition for finite length chains with sufficiently short-range interactions. This type of transition is a distinctive feature of "one-step" protein folding and our findings demonstrate that a simple homopolymer model can exhibit protein-folding thermodynamics.

Relationship between protein folding thermodynamics and the energy landscape

The origin of protein folding thermodynamics is examined in terms of the energy landscape, employing an off-lattice protein model with scaled non-native attractions, which is continuously tunable between a Go-like model and a highly frustrated system. Extensive statistical temperature molecular dynamics simulations, combined with inherent structure analysis, reveal the intimate connection between the global geometric properties of the energy landscape and the statistical temperature. The basin depth of the energy landscape is shown to play a key role in the first-order-like characteristics of the statistical temperature, which are easily identified by the squared modulus of the potential energy gradient in the microcanonical ensemble.