Inhibition of protein crystallization by evolutionary negative design

Temperature-dependent intracellular crystallization of firefly luciferase in mammalian cells is suppressed by D-luciferin and stabilizing inhibitors

Firefly luciferase (Fluc) from Photinus pyralis is one of the most widely used reporter proteins in biomedical research. Despite its widespread use, Fluc's protein phase transition behaviors and phase separation characteristics have not received much attention. Current research uncovers Fluc's intrinsic property to phase separate in mammalian cells upon a simple cell culture temperature change. Specifically, Fluc spontaneously produced needle-shaped crystal-like inclusion bodies upon temperature shift to the hypothermic temperatures ranging from 25 °C to 31 °C. The crystal-like inclusion bodies were not associated with or surrounded by membranous organelles and were likely built from the cytosolic pool of Fluc. Furthermore, the crystal-like inclusion formation was suppressed when cells were cultured in the presence of D-luciferin and its synthetic analog, as well as the benzothiazole family of so-called stabilizing inhibitors. These two classes of compounds inhibited intracellular Fluc crystallization by different modes of action as they had contrasting effects on steady-state luciferase protein accumulation levels. This study suggests that, under substrate insufficient conditions, the excess Fluc phase separates into a crystal-like state that can modulate intracellular soluble enzyme availability and protein turnover rate.

Assembly Requirements for the Construction of Large-Scale Binary Protein Structures

The precise assembly of multiple biomacromolecules into well-defined structures and materials is of great importance for various biomedical and nanobiotechnological applications. In this study, we investigate the assembly requirements for two-component materials using charged protein nanocages as building blocks. To achieve this, we designed several variants of ferritin nanocages to determine the surface characteristics necessary for the formation of large-scale binary three-dimensional (3D) assemblies. These nanocage variants were employed in protein crystallization experiments and macromolecular crystallography analyses, complemented by computational methods. Through the screening of nanocage variant combinations at various ionic strengths, we identified three essential features for successful assembly: (1) the presence of a favored crystal contact region, (2) the presence of a charged patch not involved in crystal contacts, and (3) sufficient distinctiveness between the nanocages. Surprisingly, the absence of noncrystal contact mediating patches had a detrimental effect on the assemblies, highlighting their unexpected importance. Intriguingly, we observed the formation of not only binary structures but also both negatively and positively charged unitary structures under previously exclusively binary conditions. Overall, our findings will inform future design strategies by providing some design rules, showcasing the utility of supercharging symmetric building blocks in facilitating the assembly of biomacromolecules into large-scale binary 3D assemblies.

Agglomeration: when folded proteins clump together

Protein self-association is a widespread phenomenon that results in the formation of multimeric protein structures with critical roles in cellular processes. Protein self-association can lead to finite protein complexes or open-ended, and potentially, infinite structures. This review explores the concept of protein agglomeration, a process that results from the infinite self-assembly of folded proteins. We highlight its differences from other better-described processes with similar macroscopic features, such as aggregation and liquid-liquid phase separation. We review the sequence, structural, and biophysical factors influencing protein agglomeration. Lastly, we briefly discuss the implications of agglomeration in evolution, disease, and aging. Overall, this review highlights the need to study protein agglomeration for a better understanding of cellular processes.

Characterization of insulin crystalline form in isolated β-cell secretory granules

Insulin is stored in vivo inside the pancreatic β-cell insulin secretory granules. In vitro studies have led to an assumption that high insulin and Zn²⁺ concentrations inside the pancreatic β-cell insulin secretory granules should promote insulin crystalline state in the form of Zn²⁺-stabilized hexamers. Electron microscopic images of thin sections of the pancreatic β-cells often show a dense, regular pattern core, suggesting the presence of insulin crystals. However, the structural features of the storage forms of insulin in native preparations of secretory granules are unknown, because of their small size, fragile character and difficult handling. We isolated and investigated the secretory granules from MIN6 cells under near-native conditions, using cryo-electron microscopic (Cryo-EM) techniques. The analysis of these data from multiple intra-granular crystals revealed two different rhomboidal crystal lattices. The minor lattice has unit cell parameters (a ≃ b ≃ 84.0 Å, c ≃ 35.2 Å), similar to in vitro crystallized human 4Zn²⁺-insulin hexamer, whereas the largely prevalent unit cell has more than double c-axis (a ≃ b ≃ c ≃ 96.5 Å) that probably corresponds to two or three insulin hexamers in the asymmetric unit. Our experimental data show that insulin can be present in pancreatic MIN6 cell granules in a microcrystalline form, probably consisting of 4Zn²⁺-hexamers of this hormone.

Mutant libraries reveal negative design shielding proteins from supramolecular self-assembly and relocalization in cells

Genetic mutations fuel organismal evolution but can also cause disease. As proteins are the cell’s workhorses, the ways in which mutations can disrupt their structure, stability, function, and interactions have been studied extensively. However, proteins evolve and function in a cellular context, and our ability to relate changes in protein sequence to cell-level phenotypes remains limited. In particular, the molecular mechanism underlying most disease-associated mutations is unknown. Here, we show that mutations changing a protein’s surface chemistry can dramatically impact its supramolecular self-assembly and localization in the cell. These results highlight the complex nature of genotype–phenotype relationships with a simple system.

Understanding the molecular consequences of mutations in proteins is essential to map genotypes to phenotypes and interpret the increasing wealth of genomic data. While mutations are known to disrupt protein structure and function, their potential to create new structures and localization phenotypes has not yet been mapped to a sequence space. To map this relationship, we employed two homo-oligomeric protein complexes in which the internal symmetry exacerbates the impact of mutations. We mutagenized three surface residues of each complex and monitored the mutations’ effect on localization and assembly phenotypes in yeast cells. While surface mutations are classically viewed as benign, our analysis of several hundred mutants revealed they often trigger three main phenotypes in these proteins: nuclear localization, the formation of puncta, and fibers. Strikingly, more than 50% of random mutants induced one of these phenotypes in both complexes. Analyzing the mutant’s sequences showed that surface stickiness and net charge are two key physicochemical properties associated with these changes. In one complex, more than 60% of mutants self-assembled into fibers. Such a high frequency is explained by negative design: charged residues shield the complex from self-interacting with copies of itself, and the sole removal of the charges induces its supramolecular self-assembly. A subsequent analysis of several other complexes targeted with alanine mutations suggested that such negative design is common. These results highlight that minimal perturbations in protein surfaces’ physicochemical properties can frequently drive assembly and localization changes in a cellular context.

Neuroevolutionary Learning of Particles and Protocols for Self-Assembly

Within simulations of molecules deposited on a surface we show that neuroevolutionary learning can design particles and time-dependent protocols to promote self-assembly, without input from physical concepts such as thermal equilibrium or mechanical stability and without prior knowledge of candidate or competing structures. The learning algorithm is capable of both directed and exploratory design: it can assemble a material with a user-defined property, or search for novelty in the space of specified order parameters. In the latter mode it explores the space of what can be made, rather than the space of structures that are low in energy but not necessarily kinetically accessible.

Light chain subunit of a poorly soluble human IgG2λ crystallizes in physiological pH environment both in cellulo and in vitro

Prominent inclusion bodies can develop in the endoplasmic reticulum (ER) when overexpressed antibodies possess intrinsically high condensation propensities. These observations suggest that antibodies deemed to show notable solubility problems may reveal such characteristics preemptively in the form of ER-associated inclusion bodies during antibody overexpression. To define the relationships between solubility problems and inclusion body phenotypes, we investigated the biosynthesis of a model human IgG2λ that shows severe opalescence in an acidic formulation buffer yet retains high solubility at physiological pH. Consistent with the pH-dependent solubility characteristics, the model antibody did not induce notable inclusion body in the physiological pH environment of the ER lumen. However, when individual subunit chains of the antibody were expressed separately, the light chain (LC) spontaneously induced notable crystal-like inclusion bodies in the ER. The LC crystallization event was readily reproducible in vitro by simply concentrating the purified LC protein at physiological pH. Two independent structural determinants for the LC crystallization were identified through rational mutagenesis approach by monitoring the effect of amino acid substitutions on intracellular LC crystallogenesis. The effect of mutations on crystallization was also recapitulated in vitro using purified LC proteins. Importantly, when introduced directly into the model antibody, a mutation that prevents the LC crystallization remediated the antibody's solubility problem without compromising the secretory output or antigen binding. These results illustrate that the ER can serve as a "physiological test tube" that not only reports secretory cargo's high condensation propensity at physiological pH, but also provides an orthogonal method that guides antibody engineering strategy.

Controlling Protein Crystallization by Free Energy Guided Design of Interactions at Crystal Contacts

Protein crystallization can function as an effective method for protein purification or formulation. Such an application requires a comprehensive understanding of the intermolecular protein–protein interactions that drive and stabilize protein crystal formation to ensure a reproducible process. Using alcohol dehydrogenase from Lactobacillus brevis (LbADH) as a model system, we probed in our combined experimental and computational study the effect of residue substitutions at the protein crystal contacts on the crystallizability and the contact stability. Increased or decreased contact stability was calculated using molecular dynamics (MD) free energy simulations and showed excellent qualitative correlation with experimentally determined increased or decreased crystallizability. The MD simulations allowed us to trace back the changes to their physical origins at the atomic level. Engineered charge–charge interactions as well as engineered hydrophobic effects could be characterized and were found to improve crystallizability. For example, the simulations revealed a redesigning of a water mediated electrostatic interaction (“wet contact”) into a water depleted hydrophobic effect (“dry contact”) and the optimization of a weak hydrogen bonding contact towards a strong one. These findings explained the experimentally found improved crystallizability. Our study emphasizes that it is difficult to derive simple rules for engineering crystallizability but that free energy simulations could be a very useful tool for understanding the contribution of crystal contacts for stability and furthermore could help guide protein engineering strategies to enhance crystallization for technical purposes.

Cell-penetrating Alphabody protein scaffolds for intracellular drug targeting

The therapeutic scope of antibody and nonantibody protein scaffolds is still prohibitively limited against intracellular drug targets. Here, we demonstrate that the Alphabody scaffold can be engineered into a cell-penetrating protein antagonist against induced myeloid leukemia cell differentiation protein MCL-1, an intracellular target in cancer, by grafting the critical B-cell lymphoma 2 homology 3 helix of MCL-1 onto the Alphabody and tagging the scaffold's termini with designed cell-penetration polypeptides. Introduction of an albumin-binding moiety extended the serum half-life of the engineered Alphabody to therapeutically relevant levels, and administration thereof in mouse tumor xenografts based on myeloma cell lines reduced tumor burden. Crystal structures of such a designed Alphabody in complex with MCL-1 and serum albumin provided the structural blueprint of the applied design principles. Collectively, we provide proof of concept for the use of Alphabodies against intracellular disease mediators, which, to date, have remained in the realm of small-molecule therapeutics.

Learning to grow: Control of material self-assembly using evolutionary reinforcement learning

We show that neural networks trained by evolutionary reinforcement learning can enact efficient molecular self-assembly protocols. Presented with molecular simulation trajectories, networks learn to change temperature and chemical potential in order to promote the assembly of desired structures or choose between competing polymorphs. In the first case, networks reproduce in a qualitative sense the results of previously known protocols, but faster and with higher fidelity; in the second case they identify strategies previously unknown, from which we can extract physical insight. Networks that take as input the elapsed time of the simulation or microscopic information from the system are both effective, the latter more so. The evolutionary scheme we have used is simple to implement and can be applied to a broad range of examples of experimental self-assembly, whether or not one can monitor the experiment as it proceeds. Our results have been achieved with no human input beyond the specification of which order parameter to promote, pointing the way to the design of synthesis protocols by artificial intelligence.

Infinite Assembly of Folded Proteins in Evolution, Disease, and Engineering

Mutations and changes in a protein's environment are well known for their potential to induce misfolding and aggregation, including amyloid formation. Alternatively, such perturbations can trigger new interactions that lead to the polymerization of folded proteins. In contrast to aggregation, this process does not require misfolding and, to highlight this difference, we refer to it as agglomeration. This term encompasses the amorphous assembly of folded proteins as well as the polymerization in one, two, or three dimensions. We stress the remarkable potential of symmetric homo-oligomers to agglomerate even by single surface point mutations, and we review the double-edged nature of this potential: how aberrant assemblies resulting from agglomeration can lead to disease, but also how agglomeration can serve in cellular adaptation and be exploited for the rational design of novel biomaterials.

Protein crystallization in living cells

Protein crystallization in living cells has been observed surprisingly often as a native assembly process during the past decades, and emerging evidence indicates that this phenomenon is also accessible for recombinant proteins. But only recently the advent of high-brilliance synchrotron sources, X-ray free-electron lasers, and improved serial data collection strategies has allowed the use of these micrometer-sized crystals for structural biology. Thus, in cellulo crystallization could offer exciting new possibilities for proteins that do not crystallize applying conventional approaches. In this review, we comprehensively summarize the current knowledge of intracellular protein crystallization. This includes an overview of the cellular functions, the physical properties, and, if known, the mode of regulation of native in cellulo crystal formation, complemented with a discussion of the reported crystallization events of recombinant proteins and the current method developments to successfully collect X-ray diffraction data from in cellulo crystals. Although the intracellular protein self-assembly mechanisms are still poorly understood, regulatory differences between native in cellulo crystallization linked to a specific function and accidently crystallizing proteins, either disease associated or recombinantly introduced, become evident. These insights are important to systematically exploit living cells as protein crystallization chambers in the future.

Peculiarities of Protein Crystal Nucleation and Growth

This paper reviews investigations on protein crystallization. It aims to present a comprehensive rather than complete account of recent studies and efforts to elucidate the most intimate mechanisms of protein crystal nucleation. It is emphasized that both physical and biochemical factors are at play during this process. Recently-discovered molecular scale pathways for protein crystal nucleation are considered first. The bond selection during protein crystal lattice formation, which is a typical biochemically-conditioned peculiarity of the crystallization process, is revisited. Novel approaches allow us to quantitatively describe some protein crystallization cases. Additional light is shed on the protein crystal nucleation in pores and crevices by employing the so-called EBDE method (equilibration between crystal bond and destructive energies). Also, protein crystal nucleation in solution flow is considered.

IDeAS: automated design tool for hetero-chiral protein folds

Incorporating D amino acids in the protein design alphabet can in principle multiply the design space by many orders of magnitude. All native proteins are polymers composed of L chiral amino acids. Practically limitless in diversity over amino acid sequences, protein structure is limited in folds and thus shapes, principally due to the poly L stereochemistry of their backbone. To diversify shapes, we introduced both L- and D α-amino acids as design alphabets to explore the possibility of generating novel folds, varied in chemical as well as stereo-chemical sequence. Now, to have stereochemically-defined proteins tuned chemically, we present the Inverse Design and Automation Software, IDeAS. Retro-fitting side chains on a backbone with L and D stereochemistry, the software demonstrate functional fits over stereo-chemically diverse folds in a range of applications of interest in protein design.

Availability of NHS-biotin labeling to identify free protein lysine revealed by experiment and MD simulation

It is known that the crystallizability of protein molecules may be improved by replacing their surface lysine residues with other residue types. Here an experimental method to identify surface lysine residues by NHS-biotin chemical modification combined with MALDI-TOF MS was proposed and was evaluated using PH1033 protein from Pyrococcus horikoshii. Interestingly, the biotinylation experiment with a protein-reagent molar ratio of 1:1 revealed that only seven of twenty-two lysine residues in the protein comprising 144 residues were labeled. To investigate the result, we analyzed structures from a molecular-dynamics simulation mimicking the experiment. A logistic regression analysis revealed that the biotinylation was significantly correlated with four factors relevant to the local environment of lysine residues: the solvent accessibility, the electrostatic energy, the number of hydrogen bonds, and the estimated pKa value. This result is overall in agreement with that from the same analysis on the crystal structure. However, reflecting the flexibility of the protein molecule in solution state, the factors except for the electrostatic energy were highly variable in the MD structures depending upon the protonation state of Tyr87. The present procedure of biotin-labeling can avoid lysine residues with extensive intramolecular interactions that are incompatible with the rational design of protein crystals.

Protein Crystallizability

Obtaining diffracting quality crystals remains a major challenge in protein structure research. We summarize and compare methods for selecting the best protein targets for crystallization, construct optimization and crystallization condition design. Target selection methods are divided into algorithms predicting the chance of successful progression through all stages of structural determination (from cloning to solving the structure) and those focusing only on the crystallization step. We tried to highlight pros and cons of different approaches examining the following aspects: data size, redundancy and representativeness, overfitting during model construction, and results evaluation. In summary, although in recent years progress was made and several sequence properties were reported to be relevant for crystallization, the successful prediction of protein crystallization behavior and selection of corresponding crystallization conditions continue to challenge structural researchers.

Packing topology in crystals of proteins and small molecules: a comparison

We compared the topologies of protein and small molecule crystals, which have many common features - both are molecular crystals with intermolecular interactions much weaker than intramolecular interactions. They also have different features - a considerably large fraction of the volume of protein crystals is occupied by liquid water while no room is available to other molecules in small molecule crystals. We analyzed the overall and local topology and performed multilevel topological analyses (with the software package ToposPro) of carefully selected high quality sets of protein and small molecule crystal structures. Given the suboptimal packing of protein crystals, which is due the special shape and size of proteins, it would be reasonable to expect that the topology of protein crystals is different from the topology of small molecule crystals. Surprisingly, we discovered that these two types of crystalline compounds have strikingly similar topologies. This might suggest that molecular crystal formations share symmetry rules independent of molecular dimension.

Proteins evolve on the edge of supramolecular self-assembly

The self-association of proteins into symmetric complexes is ubiquitous in all kingdoms of life. Symmetric complexes possess unique geometric and functional properties, but their internal symmetry can pose a risk. In sickle-cell disease, the symmetry of haemoglobin exacerbates the effect of a mutation, triggering assembly into harmful fibrils. Here we examine the universality of this mechanism and its relation to protein structure geometry. We introduced point mutations solely designed to increase surface hydrophobicity among 12 distinct symmetric complexes from Escherichia coli. Notably, all responded by forming supramolecular assemblies in vitro, as well as in vivo upon heterologous expression in Saccharomyces cerevisiae. Remarkably, in four cases, micrometre-long fibrils formed in vivo in response to a single point mutation. Biophysical measurements and electron microscopy revealed that mutants self-assembled in their folded states and so were not amyloid-like. Structural examination of 73 mutants identified supramolecular assembly hot spots predictable by geometry. A subsequent structural analysis of 7,471 symmetric complexes showed that geometric hot spots were buffered chemically by hydrophilic residues, suggesting a mechanism preventing mis-assembly of these regions. Thus, point mutations can frequently trigger folded proteins to self-assemble into higher-order structures. This potential is counterbalanced by negative selection and can be exploited to design nanomaterials in living cells.

Guidelines for the successful generation of protein-ligand complex crystals

With continuous technical improvements at synchrotron facilities, data-collection rates have increased dramatically. This makes it possible to collect diffraction data for hundreds of protein-ligand complexes within a day, provided that a suitable crystal system is at hand. However, developing a suitable crystal system can prove challenging, exceeding the timescale of data collection by several orders of magnitude. Firstly, a useful crystallization construct of the protein of interest needs to be chosen and its expression and purification optimized, before screening for suitable crystallization and soaking conditions can start. This article reviews recent publications analysing large data sets of crystallization trials, with the aim of identifying factors that do or do not make a good crystallization construct, and gives guidance in the design of an expression construct. It provides an overview of common protein-expression systems, addresses how ligand binding can be both help and hindrance for protein purification, and describes ligand co-crystallization and soaking, with an emphasis on troubleshooting.

The "Sticky Patch" Model of Crystallization and Modification of Proteins for Enhanced Crystallizability

Crystallization of macromolecules has long been perceived as a stochastic process, which cannot be predicted or controlled. This is consistent with another popular notion that the interactions of molecules within the crystal, i.e., crystal contacts, are essentially random and devoid of specific physicochemical features. In contrast, functionally relevant surfaces, such as oligomerization interfaces and specific protein-protein interaction sites, are under evolutionary pressures so their amino acid composition, structure, and topology are distinct. However, current theoretical and experimental studies are significantly changing our understanding of the nature of crystallization. The increasingly popular "sticky patch" model, derived from soft matter physics, describes crystallization as a process driven by interactions between select, specific surface patches, with properties thermodynamically favorable for cohesive interactions. Independent support for this model comes from various sources including structural studies and bioinformatics. Proteins that are recalcitrant to crystallization can be modified for enhanced crystallizability through chemical or mutational modification of their surface to effectively engineer "sticky patches" which would drive crystallization. Here, we discuss the current state of knowledge of the relationship between the microscopic properties of the target macromolecule and its crystallizability, focusing on the "sticky patch" model. We discuss state-of-the-art in silico methods that evaluate the propensity of a given target protein to form crystals based on these relationships, with the objective to design variants with modified molecular surface properties and enhanced crystallization propensity. We illustrate this discussion with specific cases where these approaches allowed to generate crystals suitable for structural analysis.

Minimal Positive Design for Self-Assembly of the Archimedean Tilings

A challenge of molecular self-assembly is to understand how to design particles that self-assemble into a desired structure and not any of a potentially large number of undesired structures. Here we use simulation to show that a strategy of minimal positive design allows the self-assembly of networks equivalent to the 8 semiregular Archimedean tilings of the plane, structures not previously realized in simulation. This strategy consists of identifying the fewest distinct types of interparticle interaction that appear in the desired structure, and does not require enumeration of the many possible undesired structures. The resulting particles, which self-assemble into the desired networks, possess DNA-like selectivity of their interactions. Assembly of certain molecular networks may therefore require such selectivity.

Waters in room temperature and cryo protein crystal structures

Since it has been observed that low temperature protein crystal structures may differ from room temperature structures, it is necessary to compare systematically the protein hydration structure in low and room protein crystal structures. High quality data sets of protein structures were built in an extremely rigorous manner and crystal symmetry was included in the identification of four types of water molecules (buried in the protein core, deeply inserted into crevices at the protein surface, first and second hydration layers). More water molecules are observed at low temperature only if the resolution is better than 2.1–2.3 Å. At worse resolution, temperature does not play any role. The numerous water molecules that become detectable at low temperature and at higher resolution are more mobile, relative to the protein average flexibility. Despite that, the occupancy does not depend on temperature. It can be hypothesized that water structure and around proteins and hydrogen bond network do not depend on the temperature, at least in the temperature range examined here. At low temperature more water molecules are detected because the average flexibility of all the atoms decreases, so that also water molecules that are considerably more mobile than the average atoms become observable in the electron density maps.

The role of negative selection in protein evolution revealed through the energetics of the native state ensemble

Knowing the determinants of conformational specificity is essential for understanding protein structure, stability, and fold evolution. To address this issue, a novel statistical measure of energetic compatibility between sequence and structure was developed using an experimentally validated model of the energetics of the native state ensemble. This approach successfully matched sequences from a diverse subset of the human proteome to their respective folds. Unexpectedly, significant energetic compatibility between ostensibly unrelated sequences and structures was also observed. Interrogation of these matches revealed a general framework for understanding the origins of conformational specificity within a proteome: specificity is a complex function of both the ability of a sequence to adopt folds other than the native, and ability of a fold to accommodate sequences other than the native. The regional variation in energetic compatibility indicates that the compatibility is dominated by incompatibility of sequence for alternative fold segments, suggesting that evolution of protein sequences has involved substantial negative selection, with certain segments serving as "gatekeepers" that presumably prevent alternative structures. Beyond these global trends, a size dependence exists in the degree to which the energetic compatibility is determined from negative selection, with smaller proteins displaying more negative selection. This partially explains how short sequences can adopt unique folds, despite the higher probability in shorter proteins for small numbers of mutations to increase compatibility with other folds. In providing evolutionary ground rules for the thermodynamic relationship between sequence and fold, this framework imparts valuable insight for rational design of unique folds or fold switches.

Simulation of protein association: Kinetic pathways towards crystal contacts

We conducted molecular dynamics simulations combined with distance-based umbrella sampling and forward flux sampling to investigate the early stages of protein crystallization. Formation of contacts with long-range interactions and/or an exposed position on the protein surface was kinetically preferred over more stable hydrophobic contacts with a shorter attractive range, while the thermodynamic stability of the protein crystal was provided by hydrophobic interactions. Contacts with a large interaction area showed complex dissociation pathways that were not detected by distance-based umbrella sampling. Instead, forward flux sampling simulations of contact dissociation identified long-range attractive interactions.

Recent advances in the understanding of two-step nucleation of protein crystals

The two-step mechanism of nucleation of crystals in solutions posits that the formation of crystal nuclei occurs within structures of extended lifetimes, in which the nucleating solute is at high concentration. The validity of this mechanism has been demonstrated for proteins, small-molecule organic and inorganic materials, colloids, and polymers. Due to large molecule sizes, proteins are an ideal system to study the details of this nucleation pathway, in particular the formation mechanisms of the nucleation precursors and the associated physico-chemical rules. The precursors of protein crystal nuclei are protein-rich clusters of sizes ∼100 nm that contain 10 000–100 000 molecules and occupy less than 10⁻³of the total solution volume. Here we demonstrate, using oblique illumination microscopy, the liquid nature of the clusters of the protein lysozyme and reveal their inhomogeneous structure. We test a hypothesis put forth by theory that clusters primarily consist of transient protein oligomers. For this, we explore how varying the strength of the Coulomb interaction affects the cluster characteristics. We find that the cluster’s size is insensitive to variations of pH and ionic strength. In contrast, the addition of urea, a chaotropic agent that leads to protein unfolding, strongly decreases the cluster size. Shear stress, a known protein denaturant, induced by bubbling of the solutions with an inert gas, elicits a similar response. These observations support partial protein unfolding, followed by dimerization, as the mechanism of cluster formation. The amide hydrogen–deuterium exchange, monitored by nuclear magnetic resonance, highlights that lysozyme conformational flexibility is a condition for the formation of the protein-rich clusters and facilitates the nucleation of protein crystals.

In cellulo structure determination of a novel cypovirus polyhedrin

This work demonstrates that with the use of a microfocus synchrotron beam the structure of a novel viral polyhedrin could be successfully determined from microcrystals within cells, removing the preparatory step of sample isolation and maintaining a favourable biological environment. The data obtained are of high quality, comparable to that obtained from isolated crystals, and enabled a facile structure determination. A small but significant difference is observed between the unit-cell parameters and the mosaic spread of in cellulo and isolated crystals, suggesting that even these robust crystals are adversely affected by removal from the cell.

Characterizing protein crystal contacts and their role in crystallization: rubredoxin as a case study

The fields of structural biology and soft matter have independently sought out fundamental principles to rationalize protein crystallization. Yet the conceptual differences and the limited overlap between the two disciplines have thus far prevented a comprehensive understanding of the phenomenon to emerge. We conduct a computational study of proteins from the rubredoxin family that bridges the two fields. Using atomistic simulations, we characterize the protein crystal contacts, and accordingly parameterize patchy particle models. Comparing the phase diagrams of these schematic models with experimental results enables us to critically examine the assumptions behind the two approaches. The study also reveals features of protein–protein interactions that can be leveraged to crystallize proteins more generally.

Mixing normal and anomalous diffusion

In the densely filled biological cells often subdiffusion is observed, where the average squared displacement increases slower than linear with the length of the observation interval. One reason for such subdiffusive behavior is attractive interactions between the diffusing particles that lead to temporary complex formation. Here, we show that such transient binding is not an average state of the particles but that intervals of free diffusion alternate with slower displacement when bound to neighboring particles. The observed macroscopic behavior is then the weighted average of these two contributions. Interestingly, even at very high concentrations, the unbound fraction still exhibits essentially normal diffusion.

Cellular crowding imposes global constraints on the chemistry and evolution of proteomes

In living cells, functional protein-protein interactions compete with a much larger number of nonfunctional, or promiscuous, interactions. Several cellular properties contribute to avoiding unwanted protein interactions, including regulation of gene expression, cellular compartmentalization, and high specificity and affinity of functional interactions. Here we investigate whether other mechanisms exist that shape the sequence and structure of proteins to favor their correct assembly into functional protein complexes. To examine this question, we project evolutionary and cellular abundance information onto 397, 196, and 631 proteins of known 3D structure from Escherichia coli, Saccharomyces cerevisiae, and Homo sapiens, respectively. On the basis of amino acid frequencies in interface patches versus the solvent-accessible protein surface, we define a propensity or "stickiness" scale for each of the 20 amino acids. We find that the propensity to interact in a nonspecific manner is inversely correlated with abundance. In other words, high abundance proteins have less sticky surfaces. We also find that stickiness constrains protein evolution, whereby residues in sticky surface patches are more conserved than those found in nonsticky patches. Finally, we find that the constraint imposed by stickiness on protein divergence is proportional to protein abundance, which provides mechanistic insights into the correlation between protein conservation and protein abundance. Overall, the avoidance of nonfunctional interactions significantly influences the physico-chemical and evolutionary properties of proteins. Remarkably, the effects observed are consistently larger in E. coli and S. cerevisiae than in H. sapiens, suggesting that promiscuous protein-protein interactions may be freer to accumulate in the human lineage.

Folding without charges

Surface charges of proteins have in several cases been found to function as "structural gatekeepers," which avoid unwanted interactions by negative design, for example, in the control of protein aggregation and binding. The question is then if side-chain charges, due to their desolvation penalties, play a corresponding role in protein folding by avoiding competing, misfolded traps? To find out, we removed all 32 side-chain charges from the 101-residue protein S6 from Thermus thermophilus. The results show that the charge-depleted S6 variant not only retains its native structure and cooperative folding transition, but folds also faster than the wild-type protein. In addition, charge removal unleashes pronounced aggregation on longer timescales. S6 provides thus an example where the bias toward native contacts of a naturally evolved protein sequence is independent of charges, and point at a fundamental difference in the codes for folding and intermolecular interaction: specificity in folding is governed primarily by hydrophobic packing and hydrogen bonding, whereas solubility and binding relies critically on the interplay of side-chain charges.

Generic pathways to stability in concentrated protein mixtures

We present a series of experimental results that disclose the crucial role of ionic strength and partial volume fractions in the control of the phase behaviour of binary protein mixtures. Our findings can be understood as that the ionic strength determines the relative contribution of the entropy of the protein counter-ions to the overall thermodynamics of the system. Associative phase separation and crystallization observed at, respectively, low and high ionic strength are suppressed at intermediate salt concentrations, where the entropy gain upon releasing the counter-ions from the double layer of the proteins is negligible and the entropy loss upon confining the counter-ions within the protein crystal phase significant. Moreover, we find that the partial volume fraction of the protein prone to crystallize determines the crystallization boundary and that the presence of other proteins strongly delays crystallization, leading to temporarily stable mixtures. These findings suggest that stability in more complex protein mixtures, such as the cytosol, relates to the ionic strength and protein composition rather than to protein specific properties.

Genetic construct design and recombinant protein expression for structural biology

Obtaining diffraction quality crystals is frequently an iterative process which traditionally has involved screening large numbers of crystallization conditions. Due to advances in high-throughput gene engineering, recombinant expression, and purification, the protein of interest has now become one of the many variables routinely investigated during crystallization trials. As such, construct design is a critical step in the path toward successful crystallization. In this chapter will we address construct design strategies frequently employed to improve the solution and crystallization behavior of proteins. Topics covered include choosing a recombinant expression system and reducing disorder through truncations and surface mutagenesis. Also covered are strategies to reduce heterogeneity from posttranslational modifications, impurities, and aggregation.

Participation of protein sequence termini in crystal contacts

The analysis of the crystal packing interactions, in a nonredundant set of high resolution and monomeric globular protein crystal structures, shows that the residues located at the N- and C-termini of the sequence tend to participate in packing interaction more often than expected and that often they interact with each other. Since the sequence termini are, in general, conformationally very flexible and since they host electrical charges of opposite sign, it can be hypothesized that they play a crucial role in the early formation of the nonphysiological contacts that bring to protein crystallization. It is thus not surprising that modest lengthening/shortening of the sequence termini have often a dramatic effect on protein crystallogenesis.

Evolutionary dynamics in a simple model of self-assembly

We investigate the evolutionary dynamics of an idealized model for the robust self-assembly of two-dimensional structures called polyominoes. The model includes rules that encode interactions between sets of square tiles that drive the self-assembly process. The relationship between the model's rule set and its resulting self-assembled structure can be viewed as a genotype-phenotype map and incorporated into a genetic algorithm. The rule sets evolve under selection for specified target structures. The corresponding complex fitness landscape generates rich evolutionary dynamics as a function of parameters such as the population size, search space size, mutation rate, and method of recombination. Furthermore, these systems are simple enough that in some cases the associated model genome space can be completely characterized, shedding light on how the evolutionary dynamics depends on the detailed structure of the fitness landscape. Finally, we apply the model to study the emergence of the preference for dihedral over cyclic symmetry observed for homomeric protein tetramers.

Systematic comparison of crystal and NMR protein structures deposited in the protein data bank

Nearly all the macromolecular three-dimensional structures deposited in Protein Data Bank were determined by either crystallographic (X-ray) or Nuclear Magnetic Resonance (NMR) spectroscopic methods. This paper reports a systematic comparison of the crystallographic and NMR results deposited in the files of the Protein Data Bank, in order to find out to which extent these information can be aggregated in bioinformatics. A non-redundant data set containing 109 NMR – X-ray structure pairs of nearly identical proteins was derived from the Protein Data Bank. A series of comparisons were performed by focusing the attention towards both global features and local details. It was observed that: (1) the RMDS values between NMR and crystal structures range from about 1.5 Å to about 2.5 Å; (2) the correlation between conformational deviations and residue type reveals that hydrophobic amino acids are more similar in crystal and NMR structures than hydrophilic amino acids; (3) the correlation between solvent accessibility of the residues and their conformational variability in solid state and in solution is relatively modest (correlation coefficient = 0.462); (4) beta strands on average match better between NMR and crystal structures than helices and loops; (5) conformational differences between loops are independent of crystal packing interactions in the solid state; (6) very seldom, side chains buried in the protein interior are observed to adopt different orientations in the solid state and in solution.

Protein crystallizability

Obtaining well-diffracting crystals remains a major challenge in protein structure research. In this chapter, we review currently available computational methods to estimate the crystallization potential of a protein, to optimize amino acid sequences toward improved crystallization likelihood, and to design optimal crystal screen conditions.

Application of protein engineering to enhance crystallizability and improve crystal properties

Until recently, protein crystallization has mostly been regarded as a stochastic event over which the investigator has little or no control. With the dramatic technological advances in synchrotron-radiation sources and detectors and the equally impressive progress in crystallographic software, including automated model building and validation, crystallization has increasingly become the rate-limiting step in X-ray diffraction studies of macromolecules. However, with the advent of recombinant methods it has also become possible to engineer target proteins and their complexes for higher propensity to form crystals with desirable X-ray diffraction qualities. As most proteins that are under investigation today are obtained by heterologous overexpression, these techniques hold the promise of becoming routine tools with the potential to transform classical crystallization screening into a more rational high-success-rate approach. This article presents an overview of protein-engineering methods designed to enhance crystallizability and discusses a number of examples of their successful application.

Introduction of a leucine half-zipper engenders multiple high-quality crystals of a recalcitrant tRNA synthetase

Although Escherichia coli alanyl-tRNA synthetase was among the first tRNA synthetases to be sequenced and extensively studied by functional analysis, it has proved to be recalcitrant to crystallization. This challenge remained even for crystallization of the catalytic fragment. By mutationally introducing three stacked leucines onto the solvent-exposed side of an alpha-helix, an engineered catalytic fragment of the synthetase was obtained that yielded multiple high-quality crystals and cocrystals with different ligands. The engineered alpha-helix did not form a leucine zipper that interlocked with the same alpha-helix from another molecule. Instead, using the created hydrophobic spine, it interacted with other surfaces of the protein as a leucine half-zipper (LHZ) to enhance the crystal lattice interactions. The LHZ made crystal lattice contacts in all crystals of different space groups. These results illustrate the power of introducing an LHZ into helices to facilitate crystallization. The authors propose that the method can be unified with surface-entropy reduction and can be broadly used for protein-surface optimization in crystallization.

Physicochemical principles that regulate the competition between functional and dysfunctional association of proteins

To maintain protein homeostasis, a variety of quality control mechanisms, such as the unfolded protein response and the heat shock response, enable proteins to fold and to assemble into functional complexes while avoiding the formation of aberrant and potentially harmful aggregates. We show here that a complementary contribution to the regulation of the interactions between proteins is provided by the physicochemical properties of their amino acid sequences. The results of a systematic analysis of the protein-protein complexes in the Protein Data Bank (PDB) show that interface regions are more prone to aggregate than other surface regions, indicating that many of the interactions that promote the formation of functional complexes, including hydrophobic and electrostatic forces, can potentially also cause abnormal intermolecular association. We also show, however, that aggregation-prone interfaces are prevented from triggering uncontrolled assembly by being stabilized into their functional conformations by disulfide bonds and salt bridges. These results indicate that functional and dysfunctional association of proteins are promoted by similar forces but also that they are closely regulated by the presence of specific interactions that stabilize native states.

The self-assembly of DNA Holliday junctions studied with a minimal model

In this paper, we explore the feasibility of using coarse-grained models to simulate the self-assembly of DNA nanostructures. We introduce a simple model of DNA where each nucleotide is represented by two interaction sites corresponding to the sugar-phosphate backbone and the base. Using this model, we are able to simulate the self-assembly of both DNA duplexes and Holliday junctions from single-stranded DNA. We find that assembly is most successful in the temperature window below the melting temperatures of the target structure and above the melting temperature of misbonded aggregates. Furthermore, in the case of the Holliday junction, we show how a hierarchical assembly mechanism reduces the possibility of becoming trapped in misbonded configurations. The model is also able to reproduce the relative melting temperatures of different structures accurately and allows strand displacement to occur.

Automated technologies and novel techniques to accelerate protein crystallography for structural genomics

The sequence infrastructure that has arisen through large-scale genomic projects dedicated to protein analysis, has provided a wealth of information and brought together scientists and institutions from all over the world. As a consequence, the development of novel technologies and methodologies in proteomics research is helping to unravel the biochemical and physiological mechanisms of complex multivariate diseases at both a functional and molecular level. In the late sixties, when X-ray crystallography had just been established, the idea of determining protein structure on an almost universal basis was akin to an impossible dream or a miracle. Yet only forty years after, automated protein structure determination platforms have been established. The widespread use of robotics in protein crystallography has had a huge impact at every stage of the pipeline from protein cloning, over-expression, purification, crystallization, data collection, structure solution, refinement, validation and data management- all of which have become more or less automated with minimal human intervention necessary. Here, recent advances in protein crystal structure analysis in the context of structural genomics will be discussed. In addition, this review aims to give an overview of recent developments in high throughput instrumentation, and technologies and strategies to accelerate protein structure/function analysis.

Protein crystallization in drug design: towards a rational approach

X-ray crystallography is the method of choice for the detailed characterization of stereochemistry of interactions of drug leads and potential chemotherapeutics with their protein targets. The resulting atomic models allow for rational enhancement of the lead properties and consequently for the design of high-affinity inhibitors. However, a major bottleneck of the technique is the requirement for the protein and its complexes to yield high quality single crystals. Furthermore, it is highly desirable that such crystals diffract to high resolution, preferably ≥ 1.2 Å, revealing the structures in atomic detail. Unfortunately, only a small portion of proteins readily crystallize in that fashion. New approaches are being developed to circumvent this problem. One proposed option includes rational protein surface engineering to systematically improve the crystallizability of the protein. This is accomplished by creating surface patches readily mediating weak, but specific, intermolecular interactions that take on the role of crystal contacts during nucleation and crystal growth phase.

Why are proteins charged? Networks of charge-charge interactions in proteins measured by charge ladders and capillary electrophoresis

Almost all proteins contain charged amino acids. While the function in catalysis or binding of individual charges in the active site can often be identified, it is less clear how to assign function to charges beyond this region. Are they necessary for solubility? For reasons other than solubility? Can manipulating these charges change the properties of proteins? A combination of capillary electrophoresis (CE) and protein charge ladders makes it possible to study the roles of charged residues on the surface of proteins outside the active site. This method involves chemical modification of those residues to generate a large number of derivatives of the protein that differ in charge. CE separates those derivatives into groups with the same number of modified charged groups. By studying the influence of charge on the properties of proteins using charge ladders, it is possible to estimate the net charge and hydrodynamic radius and to infer the role of charged residues in ligand binding and protein folding.