The polymerization of actin: Thermodynamics near the polymerization line

Stability matters, too - the thermodynamics of amyloid fibril formation

Amyloid fibrils are supramolecular homopolymers of proteins that play important roles in biological functions and disease. These objects have received an exponential increase in attention during the last few decades, due to their role in the aetiology of a range of severe disorders, most notably some of a neurodegenerative nature. While an overwhelming number of experimental studies exist that investigate how, and how fast, amyloid fibrils form and how their formation can be inhibited, a much more limited body of experimental work attempts to answer the question as to why these types of structures form (i.e. the thermodynamic driving force) and how stable they actually are. In this review, I attempt to give an overview of the types of experiments that have been performed to-date to answer these questions, and to summarise our current understanding of amyloid thermodynamics.

Comparison of Huggins Coefficients and Osmotic Second Virial Coefficients of Buffered Solutions of Monoclonal Antibodies

The Huggins coefficient kH is a well-known metric for quantifying the increase in solution viscosity arising from intermolecular interactions in relatively dilute macromolecular solutions, and there has been much interest in this solution property in connection with developing improved antibody therapeutics. While numerous kH measurements have been reported for select monoclonal antibodies (mAbs) solutions, there has been limited study of kH in terms of the fundamental molecular interactions that determine this property. In this paper, we compare measurements of the osmotic second virial coefficient B22, a common metric of intermolecular and interparticle interaction strength, to measurements of kH for model antibody solutions. This comparison is motivated by the seminal work of Russel for hard sphere particles having a short-range "sticky" interparticle interaction, and we also compare our data with known results for uncharged flexible polymers having variable excluded volume interactions because proteins are polypeptide chains. Our observations indicate that neither the adhesive hard sphere model, a common colloidal model of globular proteins, nor the familiar uncharged flexible polymer model, an excellent model of intrinsically disordered proteins, describes the dependence of kH of these antibodies on B22. Clearly, an improved understanding of protein and ion solvation by water as well as dipole-dipole and charge-dipole effects is required to understand the significance of kH from the standpoint of fundamental protein-protein interactions. Despite shortcomings in our theoretical understanding of kH for antibody solutions, this quantity provides a useful practical measure of the strength of interprotein interactions at elevated protein concentrations that is of direct significance for the development of antibody formulations that minimize the solution viscosity.

Hydration properties of a protein at low and high pressures: Physics of pressure denaturation

Using experimentally determined structures of ubiquitin at 1 and 3000 bar, we generate sufficiently large ensembles of model structures in the native and pressure-induced (denatured) states by means of molecular dynamics simulations with explicit water. We calculate the values of a free-energy function (FEF), which comprises the hydration free energy (HFE) and the intramolecular (conformational) energy and entropy, for the two states at 1 and 3000 bar. The HFE and the conformational entropy, respectively, are calculated using our statistical-mechanical method, which has recently been shown to be accurate, and the Boltzmann-quasi-harmonic method. The HFE is decomposed into a variety of physically insightful components. We show that the FEF of the native state is lower than that of the denatured state at 1 bar, whereas the opposite is true at 3000 bar, thus being successful in reproducing the pressure denaturation. We argue that the following two quantities of hydration play essential roles in the denaturation: the WASA-dependent term in the water-entropy loss upon cavity creation for accommodating the protein (WASA is the water-accessible surface area of the cavity) and the protein-water Lennard-Jones interaction energy. At a high pressure, the mitigation of the serious water crowding in the system is the most important, and the WASA needs to be sufficiently enlarged with the increase in the excluded-volume being kept as small as possible. The denatured structure thus induced is characterized by the water penetration into the protein interior. The pressure denaturation is accompanied by a significantly large gain of water entropy.

Supramolecular Self-assembly of a Model Hydrogelator: Characterization of Fiber Formation and Morphology

Hydrogels are of intense recent interest in connection with biomedical applications ranging from 3-D cell cultures and stem cell differentiation to regenerative medicine, controlled drug delivery and tissue engineering. This prototypical form of soft matter has many emerging material science applications outside the medical field. The physical processes underlying this type of solidification are incompletely understood and this limits design efforts aimed at optimizing these materials for applications. We address this general problem by applying multiple techniques (e.g., NMR, dynamic light scattering, small angle neutron scattering, rheological measurements) to the case of a peptide derivative hydrogelator (molecule 1, NapFFKYp) over a broad range of concentration and temperature to characterize both the formation of individual nanofibers and the fiber network. We believe that a better understanding of the hierarchical self-assembly process and control over the final morphology of this kind of material should have broad significance for biological and medicinal applications utilizing hydrogels.

Highly effective photonic cue for repulsive axonal guidance

In vivo nerve repair requires not only the ability to regenerate damaged axons, but most importantly, the ability to guide developing or regenerating axons along paths that will result in functional connections. Furthermore, basic studies in neuroscience and neuro-electronic interface design require the ability to construct in vitro neural circuitry. Both these applications require the development of a noninvasive, highly effective tool for axonal growth-cone guidance. To date, a myriad of technologies have been introduced based on chemical, electrical, mechanical, and hybrid approaches (such as electro-chemical, optofluidic flow and photo-chemical methods). These methods are either lacking in desired spatial and temporal selectivity or require the introduction of invasive external factors. Within the last fifteen years however, several attractive guidance cues have been developed using purely light based cues to achieve axonal guidance. Here, we report a novel, purely optical repulsive guidance technique that uses low power, near infrared light, and demonstrates the guidance of primary goldfish retinal ganglion cell axons through turns of up to 120 degrees and over distances of ∼90 µm.

Phase field method for nonequilibrium dynamics of reversible self-assembly systems

Phase field methods are extended to describe the nonequilibrium dynamics of reversible self-assembly systems, an extension that is complicated by the mutual coupling of many non-conserved order parameters into a set of highly nonlinear partial differential equations. Further complications arise because the sum of all non-conserved order parameters equals a conserved order parameter. The theory is developed for the simplest model of reversible self-assembly in which no additional constraints are imposed on the self-assembly process since the extension to treat more complex self-assembly models is straightforward. Specific calculations focus on the time evolution of the cluster size distribution for a free association system that is rapidly dropped from one ordered state to a more ordered state within the one-phase region. The dynamics proceed as expected, thereby providing validation of the theory which is also capable of treating systems with spatial inhomogeneities.

Response of the cell membrane-cytoskeleton complex to osmotic and freeze/thaw stresses. Part 2: The link between the state of the membrane-cytoskeleton complex and the cellular damage

In an earlier paper [35], we examined the mutual interaction between the actin cytoskeleton and the cell membrane and explored the role this interaction plays during freeze/thaw. In this follow-up paper, we investigate the physical and chemical stresses induced by freeze/thaw and explore the different mechanisms of damage caused by these stresses. Our results showed that changes in cell volume during freeze/thaw and the unfrozen water content in the solution alter the cytoskeleton stiffness, and the available membrane material. Combined with unfavorable ice-membrane interactions and increasing membrane stiffness, increased de-structuring of the membrane (such as bleb and microvilli formation) synergistically act on the membrane-cytoskeleton system generating irreversible damage.

Equilibrium polymerization models of re-entrant self-assembly

As is well known, liquid-liquid phase separation can occur either upon heating or cooling, corresponding to lower and upper critical solution phase boundaries, respectively. Likewise, self-assembly transitions from a monomeric state to an organized polymeric state can proceed either upon increasing or decreasing temperature, and the concentration dependent ordering temperature is correspondingly called the "floor" or "ceiling" temperature. Motivated by the fact that some phase separating systems exhibit closed loop phase boundaries with two critical points, the present paper analyzes self-assembly analogs of re-entrant phase separation, i.e., re-entrant self-assembly. In particular, re-entrant self-assembly transitions are demonstrated to arise in thermally activated equilibrium self-assembling systems, when thermal activation is more favorable than chain propagation, and in equilibrium self-assembly near an adsorbing boundary where strong competition exists between adsorption and self-assembly. Apparently, the competition between interactions or equilibria generally underlies re-entrant behavior in both liquid-liquid phase separation and self-assembly transitions.

Pressure effects on structures formed by entropically driven self-assembly: illustration for denaturation of proteins

We propose a general framework of pressure effects on the structures formed by the self-assembly of solute molecules immersed in solvent. The integral equation theory combined with the morphometric approach is employed for a hard-body model system. Our picture is that protein folding and ordered association of proteins are driven by the solvent entropy: At low pressures, the structures almost minimizing the excluded volume (EV) generated for solvent particles are stabilized. Such structures appear to be even more stabilized at high pressures. However, it is experimentally known that the native structure of a protein is unfolded, and ordered aggregates such as amyloid fibrils and actin filaments are dissociated by applying high pressures. This initially puzzling result can also be elucidated in terms of the solvent entropy. A clue to the basic mechanism is in the phenomenon that, when a large hard-sphere solute is immersed in small hard spheres forming the solvent, the small hard spheres are enriched near the solute and this enrichment becomes greater as the pressure increases. We argue that "attraction" is entropically provided between the solute surface and solvent particles, and the attraction becomes higher with rising pressure. Due to this effect, at high pressures, the structures possessing the largest possible solvent-accessible surface area together with sufficiently small EV become more stable in terms of the solvent entropy. To illustrate this concept, we perform an analysis of pressure denaturation of three different proteins. It is shown that only the structures that have the characteristics described above exhibit interesting behavior. They first become more destabilized relative to the native structure as the pressure increases, but beyond a threshold pressure the relative instability begins to decrease and they eventually become more stable than the native structure.

Clathrin triskelia show evidence of molecular flexibility

The clathrin triskelion, which is a three-legged pinwheel-shaped heteropolymer, is a major component in the protein coats of certain post-Golgi and endocytic vesicles. At low pH, or at physiological pH in the presence of assembly proteins, triskelia will self-assemble to form a closed clathrin cage, or "basket". Recent static light scattering and dynamic light scattering studies of triskelia in solution showed that an individual triskelion has an intrinsic pucker similar to, but differing from, that inferred from a high resolution cryoEM structure of a triskelion in a clathrin basket. We extend the earlier solution studies by performing small-angle neutron scattering (SANS) experiments on isolated triskelia, allowing us to examine a higher q range than that probed by static light scattering. Results of the SANS measurements are consistent with the light scattering measurements, but show a shoulder in the scattering function at intermediate q values (0.016 A(-1)), just beyond the Guinier regime. This feature can be accounted for by Brownian dynamics simulations based on flexible bead-spring models of a triskelion, which generate time-averaged scattering functions. Calculated scattering profiles are in good agreement with the experimental SANS profiles when the persistence length of the assumed semiflexible triskelion is close to that previously estimated from the analysis of electron micrographs.

The polymerization of actin: structural changes from small-angle neutron scattering

We present a new analysis of small-angle neutron-scattering data from rabbit muscle actin in the course of the polymerization from G-actin to F-actin as a function of temperature. The data, from Ivkov et al. [J. Chem. Phys. 108, 5599 (1998)], were taken in D2O buffer with Ca2+ as the divalent cation on the G-actin in the presence of ATP and with KCl as the initiating salt. The new analysis of the data using modeling and the method of generalized indirect fourier transform (O. Glatter, GIFT, University of Graz, Austria, http://physchem.kfunigraz.ac.at/sm/) provide shapes and dimensions of the G-actin monomer and of the growing actin oligomer in solution as a function of temperature and salt concentration. This analysis indicates that the G-actin monomer, under the conditions given above, is a sphere 50-54 A in diameter as opposed to the oblate ellipsoid seen by x-ray crystallography. The F-actin dimensions are consistent with x-ray crystal structure determinations.

The polymerization of actin: extent of polymerization under pressure, volume change of polymerization, and relaxation after temperature jumps

The protein actin can polymerize from monomeric globular G-actin to polymeric filamentary F-actin, under the regulation of thermodynamic variables such as temperature, pressure, and compositions of G-actin and salts. We present here new measurements of the extent of polymerization (phi) of actin under pressure (P), for rabbit skeletal muscle actin in H2O buffer in the presence of adenosine triposphate and calcium ions and at low (5-15 mM) KCl concentrations. We measured phi using pyrene-labeled actin, as a function of time (t) and temperature (T), for samples of fixed concentrations of initial G-actin and KCl and at fixed pressure. The phi(T,P) measurements at equilibrium have the same form as reported previously at 1 atm: low levels of polymerization at low temperatures, representing dimerization of the actin; an increase in phi at the polymerization temperature (Tp); a maximum in phi(T) above Tp) with a decrease in phi(T) beyond the maximum, indicating a depolymerization at higher T. From phi(T,P) at temperatures below Tp, we estimate the change in volume for the dimerization of actin, DeltaVdim, to be -307+/-10 ml/mol at 279 K. The change of Tp with pressure dTp/dP=(0.3015+/-0.0009) K/MPa=(30.15+/-0.09) mK/atm. The phi(T,P) data at higher T indicate the change in volume on propagation, DeltaVprop, to be +401+/-48 ml/mol at 301 K. The phi(t) measurements yield initial relaxation times rp(T) that reflect the behavior of phi(T) and support the presence of a depolymerization temperature. We also measured the density of polymerizing actin with a vibrating tube density meter, the results of which confirm that the data from this instrument are affected by viscosity changes and can be erroneous.

Actin polymerization under pressure: A theoretical study

An extended Flory-Huggins-type equilibrium polymerization theory for compressible systems is used to describe experimental data for the unusual pressure and temperature dependence of the equilibrium polymerization of G-actin to F-actin. The calculations provide rich insights into the reaction mechanism and the thermodynamics of actin polymerization at the molecular level. Volume changes associated with individual steps of the mechanism are calculated to be DeltaVactiv=(s1*-s1)upsilon0=+1553 mlmol for the activation reaction, DeltaVdim=(s2-s1*)upsilon0=-3810 mlmol for dimerization, and DeltaVprop=(sP-s1)upsilon0=+361 mlmol for the propagation reaction, where s1upsilon0, s1*upsilon0, s2upsilon0, and sPupsilon0 are the monomer volumes in the G-actin monomer, the activated G-action, the dimer, and higher polymers, respectively. Comparison with experimental measurements is made, and discrepancies are discussed.

Minimal model of relaxation in an associating fluid: Viscoelastic and dielectric relaxations in equilibrium polymer solutions

Cluster formation and disintegration greatly complicate the description of relaxation processes in complex fluids. We systematically contrast the viscoelastic and dielectric properties for models of equilibrium polymers whose thermodynamic properties have previously been established. In particular, the monomer-mediated model allows chain growth to proceed only by monomer addition, while the scission-recombination model enables all particles to associate democratically, so that chain scission and fusion occur at the interior segments as well as at chain ends. The minimal models neglect hydrodynamic and entanglement interactions and are designed to explore systematically the competition between chemical reaction and internal chain relaxation and how this coupling modifies the dynamics from that of a polydisperse solution of Rouse chains with fixed lengths (i.e., "frozen" chains). As expected, the stress relaxation is nearly single exponential when the assembly-disassembly reaction is fast on the time scale of structural chain rearrangements, while multiexponential or nearly stretched exponential relaxation is obtained when this reaction rate is slow compared to the broad relaxation spectrum of almost unperturbed, nearly "dead" chains of intrinsically polydisperse equilibrium polymer solutions. More generally, a complicated intermediate behavior emerges from the interplay between the chemical kinetic events and internal chain motions.

Symmetry, equivalence, and molecular self-assembly

Molecular self-assembly at equilibrium is fundamental to the fields of biological self-organization, the development of novel environmentally responsive polymeric materials, and nanofabrication. Our approach to understanding the principles governing this process is inspired by existing models and measurements for the self-assembly of actin, tubulin, and the ubiquitous icosahedral shell structures of viral capsids. We introduce a family of simple potentials that give rise to the self-assembly of linear polymeric, random surface ("membrane"), tubular ("nanotube"), and hollow icosahedral structures that are similar in many respects to their biological counterparts. The potentials involve equivalent particles and an interplay between directional (dipolar, multipolar) and short-range (van der Waals) interactions. Specifically, we find that the dipolar potential, having a continuous rotational symmetry about the dipolar axis, gives rise to chain formation, while particles with multipolar potentials, having discrete rotational symmetries (square quadrupole or triangular ring of dipoles or "hexapole"), lead to the self-assembly of open sheet, nanotube, and hollow icosahedral geometries. These changes in the geometry of self-assembly are accompanied by significant changes in the kinetics of the organization.

Compressible models of equilibrium polymerization

Flory-Huggins-type models of equilibrium polymerization are extended to describe compressible systems and, hence, the pressure dependence of thermodynamic properties. The theory is developed for three different mechanisms of equilibrium polymerization (the free association, monomer-activated polymerization, and chemically initiated polymerization models). In contrast to previous approaches for describing the pressure dependence, the theory delineates the thermodynamic consequences of the size disparities between solvent molecules, unpolymerized monomers, and the monomers within polymers. Basic thermodynamic properties (the extent of polymerization, density, heat capacities C(P) and C(V), etc.) are calculated analytically as functions of pressure, temperature, and composition of the associating species. Illustrative calculations refer to systems that polymerize upon cooling and demonstrate general agreement with numerous experimental trends. Comparisons with results from other theories are also discussed.

Equilibrium polymerization in the Stockmayer fluid as a model of supermolecular self-organization

A diverse range of molecular self-organization processes arises from a competition between directional and isotropic van der Waals intermolecular interactions. We conduct Monte Carlo simulations of the Stockmayer fluid (SF) with a large dipolar interaction as a minimal self-organization model and focus on basic thermodynamic properties that are needed to characterize the polymerization transition that occurs in this fluid. In particular, we determine the polymerization transition lines from the maximum in the specific heat, C(v), and the inflection point in the extent of polymerization, Phi. We also characterize the geometry (radius of gyration R(g), chain length L, chain topology) of the clusters that form in this associating fluid as a function of temperature, T, and concentration, rho . The pressure, P, and the second virial coefficient, B2, were determined, since these properties contain essential information about the strength of the isotropic (van der Waals) interactions. Our simulations indicate that the locations of the polymerization lines are quantitatively consistent with a model of equilibrium polymerization with the enthalpy of polymerization ("sticking energy") fixed by the minimum in the intermolecular potential. The polymerization transition in the SF is accompanied by a topological transition from predominantly linear to ring polymers upon cooling that is driven by the minimization of the dipolar energy of the clusters. We also find that the basic interaction parameters describing polymerization and phase separation in the SF can be estimated based on the existing theory of equilibrium polymerization, but the theory must be refined to account for ring formation in order to accurately describe the configurational properties of this model self-organizing fluid.