The self assembly of proteins; probing patchy protein interactions

Insights into the solubility of γ $$ \gamma $$ D-crystallin from multiscale atomistic simulations

The molecular basis underlying the rich phase behavior of globular proteins remains poorly understood. We use atomistic multiscale molecular simulations to model the solution-state conformational dynamics and interprotein interactions of γ $$ \gamma $$ D-crystallin and its P23T-R36S mutant, which drastically limits the protein solubility, at both infinite dilution and at a concentration where the mutant fluid phase and crystalline phase coexist. We find that while the mutant conserves the protein fold, changes to the surface exposure of residues in the neighborhood of residue-36 enhance protein-protein interactions and develop specific protein-protein contacts found in the protein crystal lattice.

Polymorphic protein phase transitions driven by surface anisotropy

Phase transitions of proteins are strongly influenced by surface chemical modifications or mutations. Human γD-crystallin (HGD) single-mutants have been extensively studied because they are associated with the onset of juvenile cataract. However, they have also provided a rich library of molecules to examine how specific inter-protein interactions direct protein assembly, providing new insights and valuable experimental data for coarse-grained patchy-particle models. Here, we demonstrate that the addition of new inter-protein interactions by mutagenesis is additive and increases the number and variety of condensed phases formed by proteins. When double mutations incorporating two specific single point mutations are made, the properties of both single mutations are retained in addition to the formation of a new condensed phase. We find that the HGD double-mutant P23VC110M self-assembles into spherical particles with retrograde solubility, orthorhombic crystals, and needle/plate shape crystals, while retaining the ability to undergo liquid-liquid phase separation. This rich polymorphism is only partially predicted by the experimental data on the constituent single mutants. We also report a previously un-characterized amorphous protein particle, with unique properties that differ from those of protein spherulites, protein particulates previously described. The particles we observe are amorphous, reversible with temperature, tens of microns in size, and perfectly spherical. When they are grown on pristine surfaces, they appear to form by homogeneous nucleation, making them unique, and we believe a new form of protein condensate. This work highlights the challenges in predicting protein behavior, which has frustrated rational assembly and crystallization but also provides rich data to develop new coarse-grained models to explain the observed polymorphism.

Surface Exposed Free Cysteine Suppresses Crystallization of Human γD-Crystallin

Human γD-crystallin (HGD) has remarkable stability against condensation in the human lens, sometimes over a whole lifetime. The native protein has a surface exposed free cysteine that forms dimers (Benedek, 1997; Ramkumar et al., 1864)^1,2 without specific biological function and leads to further protein association and/or aggregation, which creates a paradox for understanding its stability. Previous work has demonstrated that chemical modification of the protein at the free cysteine (C110), increases the temperature at which liquid-liquid phase separation occurs (LLPS), lowers protein solubility and suggests an important role for this amino acid in maintaining its long-term resistance to condensation. Here we demonstrate that mutation of the cysteine does not alter the structure or solubility (liquidus) line for the protein, but dramatically increases the protein crystal nucleation rate following LLPS, suggesting that the free cysteine has a vital role in suppressing crystallization in the human lens.

Protein-polymer mixtures in the colloid limit: Aggregation, sedimentation, and crystallization

While proteins have been treated as particles with a spherically symmetric interaction, of course in reality, the situation is rather more complex. A simple step toward higher complexity is to treat the proteins as non-spherical particles and that is the approach we pursue here. We investigate the phase behavior of the enhanced green fluorescent protein (eGFP) under the addition of a non-adsorbing polymer, polyethylene glycol. From small angle x-ray scattering, we infer that the eGFP undergoes dimerization and we treat the dimers as spherocylinders with aspect ratio L/D - 1 = 1.05. Despite the complex nature of the proteins, we find that the phase behavior is similar to that of hard spherocylinders with an ideal polymer depletant, exhibiting aggregation and, in a small region of the phase diagram, crystallization. By comparing our measurements of the onset of aggregation with predictions for hard colloids and ideal polymers [S. V. Savenko and M. Dijkstra, J. Chem. Phys. 124, 234902 (2006) and Lo Verso et al., Phys. Rev. E 73, 061407 (2006)], we find good agreement, which suggests that the behavior of the eGFP is consistent with that of hard spherocylinders and ideal polymers.

Self-assembled peptide and protein nanostructures for anti-cancer therapy: Targeted delivery, stimuli-responsive devices and immunotherapy

Self-assembled peptides and proteins possess tremendous potential as targeted drug delivery systems and key applications of these well-defined nanostructures reside in anti-cancer therapy. Peptides and proteins can self-assemble into nanostructures of diverse sizes and shapes in response to changing environmental conditions such as pH, temperature, ionic strength, as well as host and guest molecular interactions; their countless benefits include good biocompatibility and high loading capacity for hydrophobic and hydrophilic drugs. These self-assembled nanomaterials can be adorned with functional moieties to specifically target tumor cells. Stimuli-responsive features can also be incorporated with respect to the tumor microenvironment. This review sheds light on the growing interest in self-assembled peptides and proteins and their burgeoning applications in cancer treatment and immunotherapy.

Aggregation, liquid-liquid phase separation, and percolation behaviour of a model antibody fluid constrained by hard-sphere obstacles

This study is concerned with the behaviour of proteins within confinement created by hard-sphere obstacles. An individual antibody molecule is depicted as an assembly of seven hard spheres, organized to resemble a Y-shaped (on average) antibody (7-bead model) protein. For comparison with other studies we, in one case, model the protein as a hard sphere decorated by three short-range attractive sites. The antibody has two Fab and one Fc domains located in the corners of the letter Y. In this calculation, only the Fab-Fab and Fab-Fc attractive pair interactions are possible. The confinement is formed by the randomly distributed hard-sphere obstacles fixed in space. Aside from size exclusion, the obstacles do not interact with antibodies, but they affect the protein-protein correlation. We used a combination of the scaled-particle theory, Wertheim's thermodynamic perturbation theory and the Flory-Stockmayer theory to calculate: (i) the second virial coefficient of the protein fluid, (ii) the percolation threshold, (iii) cluster size distributions, and (iv) the liquid-liquid phase separation as a function of the strength of the various pair interactions of the protein and the model parameters, such as protein concentration and the packing fraction of obstacles. The conclusion is that hard-sphere obstacles strongly decrease the critical density and also, but to a much lesser extent, the critical temperature. Also, the confinement enhances clustering, making the percolating region broader. The effect depends on the model parameters, such as the packing fraction of obstacles η0, the inter-site interaction strength εIJ, and the ratio between the size of the obstacle σ0 and the size of one bead of the model antibody σhs; the value of this ratio is varied here from 2 to 5. Interestingly, at low to moderate packing fractions of obstacles, the second virial coefficient first slightly decreases (destabilization), and the slope depends on the observation temperature, but then at higher values of η0 it increases. The calculated values of the second virial coefficient also depend on the size of the obstacles.

Melting and solid-solid transitions of two-dimensional crystals composed of Janus spheres

Colloidal model systems have been extensively used in the studies of various phase transitions, but melting and solid-solid transitions have rarely been explored in monolayer colloidal crystals with anisotropic attractions. Patchy colloidal particles have served as important model systems of atoms and molecules with anisotropic interactions. In this work, we study the melting and solid-solid transitions of two-dimensional crystals composed of Janus colloidal spheres using Langevin dynamics simulation. We discovered a first-order solid-solid transition from a single crystal with uniform stripes to a novel crystal with polycrystalline domains of stripes. The centers of masses of the particles maintain the morphology of a single crystal with long-range translational and bond-orientational orders, but particle orientations form polycrystalline domains of stripes. The stripe domains form by a strain-induced nucleation process via the collective rotation of particles. In addition to this solid-solid transition, the melting transition at a higher temperature follows a two-step Kosterlitz-Thouless-Halperin-Nelson-Young (KTHNY) scenario, similar to most isotropic particle systems.

Using Schematic Models to Understand the Microscopic Basis for Inverted Solubility in γD-Crystallin

Inverted solubility-melting a crystal by cooling-is observed in a handful of proteins, such as carbomonoxy hemoglobin C and γD-crystallin. In human γD-crystallin, the phenomenon is associated with the mutation of the 23rd residue, a proline, to a threonine, serine, or valine. One proposed microscopic mechanism entails an increase in surface hydrophobicity upon mutagenesis. Recent crystal structures of a double mutant that includes the P23T mutation allow for a more careful investigation of this proposal. Here, we first measure the surface hydrophobicity of various mutant structures of γD-crystallin and discern no notable increase in hydrophobicity upon mutating the 23rd residue. We then investigate the solubility inversion regime with a schematic patchy particle model that includes one of three variants of temperature-dependent patch energies: two of the hydrophobic effect, and one of a more generic nature. We conclude that, while solubility inversion due to the hydrophobic effect may be possible, microscopic evidence to support it in γD-crystallin is weak. More generally, we find that solubility inversion requires a fine balance between patch strengths and their temperature-dependent component, which may explain why inverted solubility is not commonly observed in proteins. We also find that the temperature-dependent interaction has only a negligible impact on liquid-liquid phase boundaries of γD-crystallin, in line with previous experimental observations.

Self-assembly in amphiphilic macromolecules with solvent exposed hydrophobic moieties

Self-assembly by amphiphilic molecules with solvent exposed hydrophobic groups are relevant in biomolecular systems as well as in technological applications. Here we study such self-assembly in these systems using a model system of spherical particles having charge at core but solvent repelling surface, using Monte-Carlo simulations and mean field treatment. We find that solvophobicity mediated attraction leads aggregation, while electrostatic repulsions control stability of finite clusters. The aggregation threshold relates the parameters of two interactions through an algebraic dependence. The study also qualitatively explains experimental observations on aggregation of misfolded proteins and can be useful guide to tune stability of nm sized self-assembly in systems with exposed hydrophobic groups.

Temperature-Dependent Interactions Explain Normal and Inverted Solubility in a γD-Crystallin Mutant

Protein crystal production is a major bottleneck in the structural characterization of proteins. To advance beyond large-scale screening, rational strategies for protein crystallization are crucial. Understanding how chemical anisotropy (or patchiness) of the protein surface, due to the variety of amino-acid side chains in contact with solvent, contributes to protein-protein contact formation in the crystal lattice is a major obstacle to predicting and optimizing crystallization. The relative scarcity of sophisticated theoretical models that include sufficient detail to link collective behavior, captured in protein phase diagrams, and molecular-level details, determined from high-resolution structural information, is a further barrier. Here, we present two crystal structures for the P23T + R36S mutant of γD-crystallin, each with opposite solubility behavior: one melts when heated, the other when cooled. When combined with the protein phase diagram and a tailored patchy particle model, we show that a single temperature-dependent interaction is sufficient to stabilize the inverted solubility crystal. This contact, at the P23T substitution site, relates to a genetic cataract and reveals at a molecular level the origin of the lowered and retrograde solubility of the protein. Our results show that the approach employed here may present a productive strategy for the rationalization of protein crystallization.

Chemical Modification Alters Protein-Protein Interactions and Can Lead to Lower Protein Solubility

The chemical modification of proteins is at the frontier of developments in biological imaging and biopharmaceutics. With the advent of more sensitive and higher resolution imaging techniques,  researchers increasingly rely on the functionalization of proteins to enable these techniques to capture cellular processes. For biopharmaceutical therapies, chemically modified proteins, for example, antibody-drug conjugates (ADCs) offer the possibility of more tailored treatments for the disease with lower toxicities than traditional small molecule therapies. However, relatively little consideration is paid to how chemical modifications impact protein-protein interactions and solution stability. Using human γD-crystallin as a model, we demonstrate that chemical modification of the protein surface alters protein-protein interactions, which can result in lower solubility depending on the chemical nature of the modifier and the position on the protein where the modification is made. Understanding these effects is essential to ensure that modifying proteins effectively occurs with minimum self-association and that studies carried out using labeled proteins accurately reflect those of unmodified proteins.

Calixarene-mediated assembly of a small antifungal protein

Synthetic macrocycles such as calixarenes and cucurbiturils are increasingly applied as mediators of protein assembly and crystallization. The macrocycle can facilitate assembly by providing a surface on which two or more proteins bind simultaneously. This work explores the capacity of the sulfonato-calix[n]arene (sclx n ) series to effect crystallization of PAF, a small, cationic antifungal protein. Co-crystallization with sclx4, sclx6 or sclx8 led to high-resolution crystal structures. In the absence of sclx n , diffraction-quality crystals of PAF were not obtained. Interestingly, all three sclx n were bound to a similar patch on PAF. The largest and most flexible variant, sclx8, yielded a dimer of PAF. Complex formation was evident in solution via NMR and ITC experiments, showing more pronounced effects with increasing macrocycle size. In agreement with the crystal structure, the ITC data suggested that sclx8 acts as a bidentate ligand. The contributions of calixarene size/conformation to protein recognition and assembly are discussed. Finally, it is suggested that the conserved binding site for anionic calixarenes implicates this region of PAF in membrane binding, which is a prerequisite for antifungal activity.

Theory for the Liquid-Liquid Phase Separation in Aqueous Antibody Solutions

This study presents the theory for liquid-liquid phase separation for systems of molecules modeling monoclonal antibodies. Individual molecule is depicted as an assembly of seven hard spheres, organized to mimic the Y-shaped antibody. We consider the antibody-antibody interactions either through Fab, Fab' (two Fab fragments may be different), or Fc domain. Interaction between these three domains of the molecule (hereafter denoted as A, B, and C, respectively) is modeled by a short-range square-well attraction. To obtain numerical results for the model under study, we adapt Wertheim's thermodynamic perturbation theory. We use this model to calculate the liquid-liquid phase separation curve and the second virial coefficient B2. Various interaction scenarios are examined to see how the strength of the site-site interactions and their range shape the coexistence curve. In the asymmetric case, where an attraction between two sites is favored and the interaction energies for the other sites kept constant, critical temperature first increases and than strongly decreases. Some more microscopic information, for example, the probability for the particular two sites to be connected, has been calculated. Analysis of the experimental liquid-liquid phase diagrams, obtained from literature, is presented. In addition, we calculate the second virial coefficient under conditions leading to the liquid-liquid phase separation and present this quantity on the graph B2 versus protein concentration.

Nonclassical nucleation pathways in protein crystallization

Classical nucleation theory (CNT), which was established about 90 years ago, has been very successful in many research fields, and continues to be the most commonly used theory in describing the nucleation process. For a fluid-to-solid phase transition, CNT states that the solute molecules in a supersaturated solution reversibly form small clusters. Once the cluster size reaches a critical value, it becomes thermodynamically stable and favored for further growth. One of the most important assumptions of CNT is that the nucleation process is described by one reaction coordinate and all order parameters proceed simultaneously. Recent studies in experiments, computer simulations and theory have revealed nonclassical features in the early stage of nucleation. In particular, the decoupling of order parameters involved during a fluid-to-solid transition leads to the so-called two-step nucleation mechanism, in which a metastable intermediate phase (MIP) exists between the initial supersaturated solution and the final crystals. Depending on the exact free energy landscapes, the MIPs can be a high density liquid phase, mesoscopic clusters, or a pre-ordered state. In this review, we focus on the studies of nonclassical pathways in protein crystallization and discuss the applications of the various scenarios of two-step nucleation theory. In particular, we focus on protein solutions in the presence of multivalent salts, which serve as a model protein system to study the nucleation pathways. We wish to point out the unique features of proteins as model systems for further studies.

Protein self-assembly following in situ expression in artificial and mammalian cells

The importance of in vitro measurements in explaining the mechanisms underlying protein self-assembly in physiologically relevant conditions has been demonstrated in solution and in artificial and mammalian cells.

Modeling phase transitions in mixtures of β-γ lens crystallins

We analyze the experimentally determined phase diagram of a γD-βB1 crystallin mixture. Proteins are described as dumbbells decorated with attractive sites to allow inter-particle interaction. We use thermodynamic perturbation theory to calculate the free energy of such mixtures and, by applying equilibrium conditions, also the compositions and concentrations of the co-existing phases. Initially we fit the Tcloudversus packing fraction η measurements for a pure (x2 = 0) γD solution in 0.1 M phosphate buffer at pH = 7.0. Another piece of experimental data, used to fix the model parameters, is the isotherm x2vs. η at T = 268.5 K, at the same pH and salt content. We use the conventional Lorentz-Berthelot mixing rules to describe cross interactions. This enables us to determine: (i) model parameters for pure βB1 crystallin protein and to calculate; (ii) complete equilibrium surface (Tcloud-x2-η) for the crystallin mixtures. (iii) We present the results for several isotherms, including the tie-lines, as also the temperature-packing fraction curves. Good agreement with the available experimental data is obtained. An interesting result of these calculations is evidence of the coexistence of three phases. This domain appears for the region of temperatures just out of the experimental range studied so far. The input parameters, leading good description of experimental data, revealed a large difference between the numbers of the attractive sites for γD and βB1 proteins. This interesting result may be related to the fact that γD has a more than nine times smaller quadrupole moment than its partner in the mixture.

Relaxation dynamics of functionalized colloids on attractive substrates

Particle-based simulations are performed to study the post-relaxation dynamics of functionalized (patchy) colloids adsorbed on an attractive substrate. Kinetically arrested structures that depend on the number of adsorbed particles and the strength of the particle-particle and particle-substrate interactions are identified. The radial distribution function is characterized by a sequence of peaks, with relative intensities that depend on the number of adsorbed particles. The first-layer coverage is a non-monotonic function of the number of particles, with an optimal value around one layer of adsorbed particles. The initial relaxation towards these structures is characterized by a fast (exponential) and a slow (power-law) dynamics. The fast relaxation timescale is a linearly increasing function of the number of adsorbed particles in the submonolayer regime, but it saturates for more than one adsorbed layer. The slow dynamics exhibits two characteristic exponents, depending on the surface coverage.

How fluorescent labelling alters the solution behaviour of proteins

The impact of adding a fluorescent dye to a protein is examined using protein phase diagrams and numerical simulations. The addition of the dye is the equivalent of adding a large attractive patch to the surface of the protein, which results in significant changes to the protein solution behaviour.