The elementary reactions for incorporation into crystals.

The growth rates of crystals are largely dictated by the chemical reaction between solute and kinks, in which a solute molecule severs its bonds with the solvent and establishes new bonds with the kink. Details on this sequence of bond breaking and rebuilding remain poorly understood. To elucidate the reaction at the kinks we employ four solvents with distinct functionalities as reporters on the microscopic structures and their dynamics along the pathway into a kink. We combine time-resolved in situ atomic force microscopy and x-ray and optical methods with molecular dynamics simulations. We demonstrate that in all four solvents the solute, etioporphyrin I, molecules reach the steps directly from the solution; this finding identifies the measured rate constant for step growth as the rate constant of the reaction between a solute molecule and a kink. We show that the binding of a solute molecule to a kink divides into two elementary reactions. First, the incoming solute molecule sheds a fraction of its solvent shell and attaches to molecules from the kink by bonds distinct from those in its fully incorporated state. In the second step, the solute breaks these initial bonds and relocates to the kink. The strength of the preliminary bonds with the kink determines the free energy barrier for incorporation into a kink. The presence of an intermediate state, whose stability is controlled by solvents and additives, may illuminate how minor solution components guide the construction of elaborate crystal architectures in nature and the search for solution compositions that suppress undesirable or accelerate favored crystallization in industry.

Antagonistic cooperativity between crystal growth modifiers.

Ubiquitous processes in nature and the industry exploit crystallization from multicomponent environments1-5; however, laboratory efforts have focused on the crystallization of pure solutes6,7 and the effects of single growth modifiers8,9. Here we examine the molecular mechanisms employed by pairs of inhibitors in blocking the crystallization of haematin, which is a model organic compound with relevance to the physiology of malaria parasites10,11. We use a combination of scanning probe microscopy and molecular modelling to demonstrate that inhibitor pairs, whose constituents adopt distinct mechanisms of haematin growth inhibition, kink blocking and step pinning12,13, exhibit both synergistic and antagonistic cooperativity depending on the inhibitor combination and applied concentrations. Synergism between two crystal growth modifiers is expected, but the antagonistic cooperativity of haematin inhibitors is not reflected in current crystal growth models. We demonstrate that kink blockers reduce the line tension of step edges, which facilitates both the nucleation of crystal layers and step propagation through the gates created by step pinners. The molecular viewpoint on cooperativity between crystallization modifiers provides guidance on the pairing of modifiers in the synthesis of crystalline materials. The proposed mechanisms indicate strategies to understand and control crystallization in both natural and engineered systems, which occurs in complex multicomponent media1-3,8,9. In a broader context, our results highlight the complexity of crystal-modifier interactions mediated by the structure and dynamics of the crystal interface.

A minimal colloid model of solution crystallization nucleates crystals classically.

A fundamental assumption of the classical theories of crystal nucleation is that the individual molecules from the "old" phase associate to an emerging nucleus individually and sequentially. Numerous recent studies of crystal nucleation in solution have revealed nonclassical pathways, whereby crystal nuclei are hosted and fed by amorphous clusters pre-formed in the solution. A sizable knowledge gap has persisted, however, in the definition of the molecular-level parameters that direct a solute towards classical or nonclassical nucleation. Here we construct a suspension of colloid particles of hydrodynamic diameter 1.1 µm and monitor their individual motions towards a quasi-two-dimensional crystal by scanning confocal microscopy. We combine electrostatic repulsion and polymer-induced attraction to obtain a simple isotropic pair interaction potential with a single attractive minimum of tunable depth between 1.2kBT and 2.7kBT. We find that even the smallest aggregates that form in this system structure as hexagonal two-dimensional crystals and grow and maturate by the association and exchange of single particles from the solution, signature behaviors during classical nucleation. The particles in the suspension equilibrate with those in the clusters and the volume fractions of suspensions at equilibrium correspond to straightforward thermodynamic predictions based on depth of the interparticle attraction. These results demonstrate that classical nucleation is selected by particles interacting with a minimal potential and present a benchmark for future modifications of the molecular interactions that may induce nonclassical nucleation.

Frustrated peptide chains at the fibril tip control the kinetics of growth of amyloid-ß fibrils.

Amyloid fibrillization is an exceedingly complex process in which incoming peptide chains bind to the fibril while concertedly folding. The coupling between folding and binding is not fully understood. We explore the molecular pathways of association of Aß40 monomers to fibril tips by combining time-resolved in situ scanning probe microscopy with molecular modeling. The comparison between experimental and simulation results shows that a complex supported by nonnative contacts is present in the equilibrium structure of the fibril tip and impedes fibril growth in a supersaturated solution. The unraveling of this frustrated state determines the rate of fibril growth. The kinetics of growth of freshly cut fibrils, in which the bulk fibril structure persists at the tip, complemented by molecular simulations, indicate that this frustrated complex comprises three or four monomers in nonnative conformations and likely is contained on the top of a single stack of peptide chains in the fibril structure. This pathway of fibril growth strongly deviates from the common view that the conformational transformation of each captured peptide chain is templated by the previously arrived peptide. The insights into the ensemble structure of the frustrated complex may guide the search for suppressors of Aß fibrillization. The uncovered dynamics of coupled structuring and assembly during fibril growth are more complex than during the folding of most globular proteins, as they involve the collective motions of several peptide chains that are not guided by a funneled energy landscape.

Mesoscopic protein-rich clusters host the nucleation of mutant p53 amyloid fibrils.

The protein p53 is a crucial tumor suppressor, often called "the guardian of the genome"; however, mutations transform p53 into a powerful cancer promoter. The oncogenic capacity of mutant p53 has been ascribed to enhanced propensity to fibrillize and recruit other cancer fighting proteins in the fibrils, yet the pathways of fibril nucleation and growth remain obscure. Here, we combine immunofluorescence three-dimensional confocal microscopy of human breast cancer cells with light scattering and transmission electron microscopy of solutions of the purified protein and molecular simulations to illuminate the mechanisms of phase transformations across multiple length scales, from cellular to molecular. We report that the p53 mutant R248Q (R, arginine; Q, glutamine) forms, both in cancer cells and in solutions, a condensate with unique properties, mesoscopic protein-rich clusters. The clusters dramatically diverge from other protein condensates. The cluster sizes are decoupled from the total cluster population volume and independent of the p53 concentration and the solution concentration at equilibrium with the clusters varies. We demonstrate that the clusters carry out a crucial biological function: they host and facilitate the nucleation of amyloid fibrils. We demonstrate that the p53 clusters are driven by structural destabilization of the core domain and not by interactions of its extensive unstructured region, in contradistinction to the dense liquids typical of disordered and partially disordered proteins. Two-step nucleation of mutant p53 amyloids suggests means to control fibrillization and the associated pathologies through modifying the cluster characteristics. Our findings exemplify interactions between distinct protein phases that activate complex physicochemical mechanisms operating in biological systems.

Suppression of amyloid-ß fibril growth by drug-engineered polymorph transformation.

Fibrillization of the protein amyloid ß is assumed to trigger Alzheimer's pathology. Approaches that target amyloid plaques, however, have garnered limited clinical success, and their failures may relate to the scarce understanding of the impact of potential drugs on the intertwined stages of fibrillization. Here, we demonstrate that bexarotene, a T-cell lymphoma medication with known antiamyloid activity both in vitro and in vivo, suppresses amyloid fibrillization by promoting an alternative fibril structure. We employ time-resolved in situ atomic force microscopy to quantify the kinetics of growth of individual fibrils and supplement it with structure characterization by cryo-EM. We show that fibrils with structure engineered by the drug nucleate and grow substantially slower than "normal" fibrils; remarkably, growth remains stunted even in drug-free solutions. We find that the suppression of fibril growth by bexarotene is not because of the drug binding to the fibril tips or to the peptides in the solution. Kinetic analyses attribute the slow growth of drug-enforced fibril polymorph to the distinctive dynamics of peptide chain association to their tips. As an additional benefit, the bexarotene fibrils kill primary rat hippocampal neurons less efficiently than normal fibrils. In conclusion, the suggested drug-driven polymorph transformation presents a mode of action to irreversibly suppress toxic aggregates not only in Alzheimer's but also potentially in myriad diverse pathologies that originate with protein condensation.

A second mechanism employed by artemisinins to suppress Plasmodium falciparum hinges on inhibition of hematin crystallization.

Malaria is a pervasive disease that affects millions of lives each year in equatorial regions of the world. During the erythrocytic phase of the parasite life cycle, Plasmodium falciparum invades red blood cells, where it catabolizes hemoglobin and sequesters the released toxic heme as innocuous hemozoin crystals. Artemisinin (ART)-class drugs are activated in vivo by newly released heme, which creates a carbon-centered radical that markedly reduces parasite density. Radical damage to parasite lipids and proteins is perceived to be ARTs' dominant mechanism of action. By contrast, quinoline-class antimalarials inhibit the formation of hemozoin and in this way suppress heme detoxification. Here, we combine malaria parasite assays and scanning probe microscopy of growing ß-hematin crystals to elucidate an unexpected mechanism employed by two widely administered antimalarials, ART, and artesunate to subdue the erythrocytic phase of the parasite life cycle. We demonstrate that heme-drug adducts, produced after the radical activation of ARTs and largely believed to be benign bystanders, potently kills P. falciparum at low exogenous concentrations. We show that these adducts inhibit ß-hematin crystallization and heme detoxification, a pathway which complements the deleterious effect of radicals generated via parent drug activation. Our findings reveal an irreversible mechanism of heme-ART adduct inhibition of heme crystallization, unique among antimalarials and common crystal growth inhibitors, that opens new avenues for evaluating drug dosing regimens and understanding growing resistance of P. falciparum to ART.

Precrystallization solute assemblies and crystal symmetry.

Solution crystallization is a part of the synthesis of materials ranging from geological and biological minerals to pharmaceuticals, fine chemicals, and advanced electronic components. Attempts to predict the structure, growth rates and properties of emerging crystals have been frustrated, in part, by the poor understanding of the correlations between the oligomeric state of the solute, the growth unit, and the crystal symmetry. To explore how a solute monomer or oligomer is selected as the unit that incorporates into kinks and how crystal symmetry impacts this selection, we combine scanning probe microscopy, optical spectroscopy, and all-atom molecular simulations using as examples two organic materials, olanzapine (OZPN) and etioporphyrin I (EtpI). The dominance of dimeric structures in OZPN crystals has spurred speculation that the dimers preform in the solution, where they capture the majority of the solute, and then assemble into crystals. By contrast, EtpI in crystals aligns in parallel stacks of flat EtpI monomers unrelated by point symmetry. Raman and absorption spectroscopies show that solute monomers are the majority solute species in solutions of both compounds. Surprisingly, the kinetics of incorporation of OZPN into kinks is bimolecular, indicating that the growth unit is a solute dimer, a minority solution component. The disconnection between the dominant solute species, the growth unit, and the crystal symmetry is even stronger with EtpI, for which the (010) face grows by incorporating monomers, whereas the growth unit of the (001) face is a dimer. Collectively, the crystallization kinetics results with OZPN and EtpI establish that the structures of the dominant solute species and of the incorporating solute complex do not correlate with the symmetry of the crystal lattice. In a broader context, these findings illuminate the immense complexity of crystallization scenarios that need to be explored on the road to the understanding and control of crystallization.

Antimalarials inhibit hematin crystallization by unique drug-surface site interactions.

In malaria pathophysiology, divergent hypotheses on the inhibition of hematin crystallization posit that drugs act either by the sequestration of soluble hematin or their interaction with crystal surfaces. We use physiologically relevant, time-resolved in situ surface observations and show that quinoline antimalarials inhibit ß-hematin crystal surfaces by three distinct modes of action: step pinning, kink blocking, and step bunch induction. Detailed experimental evidence of kink blocking validates classical theory and demonstrates that this mechanism is not the most effective inhibition pathway. Quinolines also form various complexes with soluble hematin, but complexation is insufficient to suppress heme detoxification and is a poor indicator of drug specificity. Collectively, our findings reveal the significance of drug-crystal interactions and open avenues for rationally designing antimalarial compounds.

Two types of amorphous protein particles facilitate crystal nucleation.

Nucleation, the primary step in crystallization, dictates the number of crystals, the distribution of their sizes, the polymorph selection, and other crucial properties of the crystal population. We used time-resolved liquid-cell transmission electron microscopy (TEM) to perform an in situ examination of the nucleation of lysozyme crystals. Our TEM images revealed that mesoscopic clusters, which are similar to those previously assumed to consist of a dense liquid and serve as nucleation precursors, are actually amorphous solid particles (ASPs) and act only as heterogeneous nucleation sites. Crystalline phases never form inside them. We demonstrate that a crystal appears within a noncrystalline particle assembling lysozyme on an ASP or a container wall, highlighting the role of heterogeneous nucleation. These findings represent a significant departure from the existing formulation of the two-step nucleation mechanism while reaffirming the role of noncrystalline particles. The insights gained may have significant implications in areas that rely on the production of protein crystals, such as structural biology, pharmacy, and biophysics, and for the fundamental understanding of crystallization mechanisms.

Mechanisms of hematin crystallization and inhibition by the antimalarial drug chloroquine.

Hematin crystallization is the primary mechanism of heme detoxification in malaria parasites and the target of the quinoline class of antimalarials. Despite numerous studies of malaria pathophysiology, fundamental questions regarding hematin growth and inhibition remain. Among them are the identity of the crystallization medium in vivo, aqueous or organic; the mechanism of crystallization, classical or nonclassical; and whether quinoline antimalarials inhibit crystallization by sequestering hematin in the solution, or by blocking surface sites crucial for growth. Here we use time-resolved in situ atomic force microscopy (AFM) and show that the lipid subphase in the parasite may be a preferred growth medium. We provide, to our knowledge, the first evidence of the molecular mechanisms of hematin crystallization and inhibition by chloroquine, a common quinoline antimalarial drug. AFM observations demonstrate that crystallization strictly follows a classical mechanism wherein new crystal layers are generated by 2D nucleation and grow by the attachment of solute molecules. We identify four classes of surface sites available for binding of potential drugs and propose respective mechanisms of drug action. Further studies reveal that chloroquine inhibits hematin crystallization by binding to molecularly flat {100} surfaces. A 2-µM concentration of chloroquine fully arrests layer generation and step advancement, which is ∼10(4)× less than hematin's physiological concentration. Our results suggest that adsorption at specific growth sites may be a general mode of hemozoin growth inhibition for the quinoline antimalarials. Because the atomic structures of the identified sites are known, this insight could advance the future design and/or optimization of new antimalarials.

Early Onset of Kinetic Roughening due to a Finite Step Width in Hematin Crystallization.

The structure of the interface of a growing crystal with its nutrient phase largely determines the growth dynamics. We demonstrate that hematin crystals, crucial for the survival of malaria parasites, transition from faceted to rough growth interfaces at increasing thermodynamic supersaturation Δµ. Contrary to theoretical predictions and previous observations, this transition occurs at moderate values of Δµ. Moreover, surface roughness varies nonmonotonically with Δµ, and the rate constant for rough growth is slower than that resulting from nucleation and spreading of layers. We attribute these unexpected behaviors to the dynamics of step growth dominated by surface diffusion and the loss of identity of nuclei separated by less than the step width w. We put forth a general criterion for the onset of kinetic roughening using w as a critical length scale.

Deconstructing Quinoline-Class Antimalarials to Identify Fundamental Physicochemical Properties of Beta-Hematin Crystal Growth Inhibitors.

A versatile approach to control crystallization involves the use of modifiers, which are additives that interact with crystal surfaces and alter their growth rates. Elucidating a modifier's binding specificity to anisotropic crystal surfaces is a ubiquitous challenge that is critical to their design. In this study, we select hematin, a byproduct of malaria parasites, as a model system to examine the complementarity of modifiers (i.e., antimalarial drugs) to ß-hematin crystal surfaces. We divide two antimalarials, chloroquine and amodiaquine, into segments consisting of a quinoline base, common to both drugs, and side chains that differentiate their modes of action. Using a combination of scanning probe microscopy, bulk crystallization, and analytical techniques, we show that the base and side chain work synergistically to reduce the rate of hematin crystallization. In contrast to general observations that modifiers retain their function upon segmentation, we show that the constituents do not act as modifiers. A systematic study of quinoline isomers and analogues shows how subtle rearrangement and removal of functional moieties can create effective constituents from previously ineffective modifiers, along with tuning their inhibitory modes of action. These findings highlight the importance of specific functional moieties in drug compounds, leading to an improved understanding of modifier-crystal interactions that could prove to be applicable to the design of new antimalarials.

Probing the Twisted Structure of Sickle Hemoglobin Fibers via Particle Simulations.

Polymerization of sickle hemoglobin (HbS) is the primary pathogenic event of sickle cell disease. For insight into the nature of the HbS polymer fiber formation, we develop a particle model-resembling a coarse-grained molecular model-constructed to match the intermolecular contacts between HbS molecules. We demonstrate that the particle model predicts the formation of HbS polymer fibers by attachment of monomers to rough fiber ends and the growth rate increases linearly with HbS concentration. We show that the characteristic 14-molecule fiber cross section is preserved during growth. We also correlate the asymmetry of the contact sites on the HbS molecular surface with the structure of the polymer fiber composed of seven helically twisted double strands. Finally, we show that the same asymmetry mediates the mechanical and structural properties of the HbS polymer fiber.

Characterization of the diffusive dynamics of particles with time-dependent asymmetric microscopy intensity profiles.

We put forth an algorithm to track isolated micron-size solid and liquid particles that produce time-dependent asymmetric intensity patterns. This method quantifies the displacement of a particle in the image plane from the peak of a spatial cross-correlation function with a reference image. The peak sharpness results in subpixel resolution. We demonstrate the utility of the method for tracking liquid droplets with changing shapes and micron-size particles producing images with exaggerated asymmetry. We compare the accuracy of diffusivity determination with particles of known size by this method to that by common tracking techniques and demonstrate that our algorithm is superior. We address several open questions on the characterization of diffusive behaviors. We show that for particles, diffusing with a root-mean-square displacement of 0.6 pixel widths in the time between two successive recorded frames, more accurate diffusivity determinations result from mean squared displacement (MSD) for lag times up to 5 time intervals and that MSDs determined from non-overlapping displacements do not yield more accurate diffusivities. We discuss the optimal length of image sequences and demonstrate that lower frame rates do not affect the accuracy of the estimated diffusivity.

Lack of Dependence of the Sizes of the Mesoscopic Protein Clusters on Electrostatics.

Protein-rich clusters of steady submicron size and narrow size distribution exist in protein solutions in apparent violation of the classical laws of phase equilibrium. Even though they contain a minor fraction of the total protein, evidence suggests that they may serve as essential precursors for the nucleation of ordered solids such as crystals, sickle-cell hemoglobin polymers, and amyloid fibrils. The cluster formation mechanism remains elusive. We use the highly basic protein lysozyme at nearly neutral and lower pH as a model and explore the response of the cluster population to the electrostatic forces, which govern numerous biophysical phenomena, including crystallization and fibrillization. We tune the strength of intermolecular electrostatic forces by varying the solution ionic strength I and pH and find that despite the weaker repulsion at higher I and pH, the cluster size remains constant. Cluster responses to the presence of urea and ethanol demonstrate that cluster formation is controlled by hydrophobic interactions between the peptide backbones, exposed to the solvent after partial protein unfolding that may lead to transient protein oligomers. These findings reveal that the mechanism of the mesoscopic clusters is fundamentally different from those underlying the two main classes of ordered protein solid phases, crystals and amyloid fibrils, and partial unfolding of the protein chain may play a significant role.

The free heme concentration in healthy human erythrocytes.

Heme, the prosthetic group of hemoglobin, may be released from its host due to an intrinsic instability of hemoglobin and accumulate in the erythrocytes. Free heme is in the form of hematin (Fe(3+) protoporphyrin IX OH) and follows several pathways of biochemical toxicity to tissues, cells, and organelles since it catalyzes the production of reactive oxygen species. To determine concentration of soluble free heme in human erythrocytes, we develop a new method. We lyse the red blood cells and isolate free heme from hemoglobin by dialysis. We use the heme to reconstitute horseradish peroxidase (HRP) from an excess of the apoenzyme and determine the HRP reaction rate from the evolution of the emitted luminescence. We find that in a population of five healthy adults the average free heme concentration in the erythrocytes is 21±2µM, ca. 100× higher than previously determined. Tests suggest that the lower previous value was due to the use of elevated concentrations of NaCl, which drive hematin precipitation and re-association with apoglobin. We show that the found hematin concentration is significantly higher than estimates based on equilibrium release and the known hematin dimerization. The factors that lead to enhanced heme release remain an open question.

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(-3) 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.