1
|
Chakrabarti R, Verma L, Hadjiev VG, Palmer JC, Vekilov PG. The elementary reactions for incorporation into crystals. Proc Natl Acad Sci U S A 2024; 121:e2320201121. [PMID: 38315836 PMCID: PMC10873555 DOI: 10.1073/pnas.2320201121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2023] [Accepted: 12/26/2023] [Indexed: 02/07/2024] Open
Abstract
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.
Collapse
Affiliation(s)
- Rajshree Chakrabarti
- William A. Brookshire Department of Chemical and Biomolecular Engineering, University of Houston, Houston, TX77204-4004
| | - Lakshmanji Verma
- William A. Brookshire Department of Chemical and Biomolecular Engineering, University of Houston, Houston, TX77204-4004
| | - Viktor G. Hadjiev
- Texas Center for Superconductivity, University of Houston, Houston, TX77004-50024
| | - Jeremy C. Palmer
- William A. Brookshire Department of Chemical and Biomolecular Engineering, University of Houston, Houston, TX77204-4004
| | - Peter G. Vekilov
- William A. Brookshire Department of Chemical and Biomolecular Engineering, University of Houston, Houston, TX77204-4004
- Department of Chemistry, University of Houston, Houston, TX77204-5003
| |
Collapse
|
2
|
Warzecha M, Verma L, Chakrabarti R, Hadjiev VG, Florence AJ, Palmer JC, Vekilov PG. Precrystallization solute assemblies and crystal symmetry. Faraday Discuss 2022; 235:307-321. [PMID: 35393981 DOI: 10.1039/d1fd00080b] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
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.
Collapse
Affiliation(s)
- Monika Warzecha
- EPSRC, CMAC, Future Manufacturing Research Hub, c/o Strathclyde Institute of Pharmacy and Biomedical Sciences, Technology and Innovation Centre, 99 George Street, Glasgow, G1 1RD, UK
| | - Lakshmanji Verma
- Department of Chemical and Biomolecular Engineering, University of Houston, 4726 Calhoun Rd., Houston, TX 77204-4004, USA.
| | - Rajshree Chakrabarti
- Department of Chemical and Biomolecular Engineering, University of Houston, 4726 Calhoun Rd., Houston, TX 77204-4004, USA.
| | - Viktor G Hadjiev
- Texas Center for Superconductivity, University of Houston, 3369 Cullen Blvd., Suite 202, Houston, Texas 77004-50024, USA
| | - Alastair J Florence
- EPSRC, CMAC, Future Manufacturing Research Hub, c/o Strathclyde Institute of Pharmacy and Biomedical Sciences, Technology and Innovation Centre, 99 George Street, Glasgow, G1 1RD, UK
| | - Jeremy C Palmer
- Department of Chemical and Biomolecular Engineering, University of Houston, 4726 Calhoun Rd., Houston, TX 77204-4004, USA.
| | - Peter G Vekilov
- Department of Chemical and Biomolecular Engineering, University of Houston, 4726 Calhoun Rd., Houston, TX 77204-4004, USA. .,Department of Chemistry, University of Houston, 3585 Cullen Blvd., Houston, TX 77204-5003, USA
| |
Collapse
|
3
|
Verma L, Warzecha M, Chakrabarti R, Hadjiev VG, Palmer JC, Vekilov PG. How to Identify the Crystal Growth Unit. Isr J Chem 2021. [DOI: 10.1002/ijch.202100081] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Affiliation(s)
- Lakshmanji Verma
- Department of Chemical and Biomolecular Engineering University of Houston 4726 Calhoun Rd. Houston, TX 77204-4004 USA
| | - Monika Warzecha
- EPSRC CMAC Future Manufacturing Research Hub c/o Strathclyde Institute of Pharmacy and Biomedical Sciences Technology and Innovation Centre 99 George Street Glasgow G1 1RD U.K
| | - Rajshree Chakrabarti
- Department of Chemical and Biomolecular Engineering University of Houston 4726 Calhoun Rd. Houston, TX 77204-4004 USA
| | - Viktor G. Hadjiev
- Texas Center for Superconductivity University of Houston 3369 Cullen Blvd., Suite 202 Houston Texas 77004-50024 USA
| | - Jeremy C. Palmer
- Department of Chemical and Biomolecular Engineering University of Houston 4726 Calhoun Rd. Houston, TX 77204-4004 USA
| | - Peter G. Vekilov
- Department of Chemical and Biomolecular Engineering University of Houston 4726 Calhoun Rd. Houston, TX 77204-4004 USA
- Department of Chemistry University of Houston 3585 Cullen Blvd. Houston, TX 77204-5003 USA
| |
Collapse
|
4
|
The Ambiguous Functions of the Precursors That Enable Nonclassical Modes of Olanzapine Nucleation and Growth. CRYSTALS 2021. [DOI: 10.3390/cryst11070738] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
One of the most consequential assumptions of the classical theories of crystal nucleation and growth is the Szilard postulate, which states that molecules from a supersaturated phase join a nucleus or a growing crystal individually. In the last 20 years, observations in complex biological, geological, and engineered environments have brought to light violations of the Szilard rule, whereby molecules assemble into ordered or disordered precursors that then host and promote nucleation or contribute to fast crystal growth. Nonclassical crystallization has risen to a default mode presumed to operate in the majority of the inspected crystallizing systems. In some cases, the existence of precursors in the growth media is admitted as proof for their role in nucleation and growth. With the example of olanzapine, a marketed drug for schizophrenia and bipolar disorder, we demonstrate that molecular assemblies in the solution selectively participate in crystal nucleation and growth. In aqueous and organic solutions, olanzapine assembles into both mesoscopic solute-rich clusters and dimers. The clusters facilitate nucleation of crystals and crystal form transformations. During growth, however, the clusters land on the crystal surface and transform into defects, but do not support step growth. The dimers are present at low concentrations in the supersaturated solution, yet the crystals grow by the association of dimers, and not of the majority monomers. The observations with olanzapine emphasize that detailed studies of the crystal and solution structures and the dynamics of molecular association may empower classical and nonclassical models that advance the understanding of natural crystallization, and support the design and manufacture of promising functional materials.
Collapse
|
5
|
Seckfort D, Montgomery Pettitt B. Price of disorder in the lac repressor hinge helix. Biopolymers 2018; 110:e23239. [PMID: 30485404 DOI: 10.1002/bip.23239] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2018] [Revised: 09/12/2018] [Accepted: 10/04/2018] [Indexed: 12/26/2022]
Abstract
The Lac system of genes has been pivotal in understanding gene regulation. When the lac repressor protein binds to the correct DNA sequence, the hinge region of the protein goes through a disorder to order transition. The structure of this region of the protein is well understood when it is in this bound conformation, but less so when it is not. Structural studies show that this region is flexible. Our simulations show this region is extremely flexible in solution; however, a high concentration of salt can help kinetically trap the hinge helix. Thermodynamically, disorder is more favorable without the DNA present.
Collapse
Affiliation(s)
- Danielle Seckfort
- Structural and Computational Biology and Molecular Biophysics, Baylor College of Medicine, Houston, Texas
| | - B Montgomery Pettitt
- Structural and Computational Biology and Molecular Biophysics, Baylor College of Medicine, Houston, Texas.,Department of Biochemistry and Molecular Biology, Sealy Center for Structural Biology, University of Texas Medical Branch, Galveston, Texas
| |
Collapse
|
6
|
Pereverzev YV, Prezhdo OV. Universal laws in the force-induced unraveling of biological bonds. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2007; 75:011905. [PMID: 17358182 DOI: 10.1103/physreve.75.011905] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/11/2005] [Revised: 08/09/2006] [Indexed: 05/14/2023]
Abstract
Universal laws in the force-induced unbinding of receptor-ligand complexes are established for a general functional dependence of the dissociation rate constant on the applied force and are detailed with the two-pathway model that describes the recently discovered biological catch bond. The relationships link the data obtained with constant and time-dependent forces in different regimes, provide common representation for the previously unrelated data sets, and, thereby, greatly facilitate analysis and interpretation of experiments. The universal laws are demonstrated with the monomeric and dimeric catch-slip bonds between P-selectins and P-selectin glycoprotein ligands-1, and the slip bond between E-selectin and sialyl Lewis;{x} antigen.
Collapse
Affiliation(s)
- Yuriy V Pereverzev
- Department of Chemistry, University of Washington, Seattle, Washington 98195-1700, USA
| | | |
Collapse
|
7
|
Gasperino D, Yeckel A, Olmsted BK, Ward MD, Derby JJ. Mass transfer limitations at crystallizing interfaces in an atomic force microscopy fluid cell: a finite element analysis. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2006; 22:6578-86. [PMID: 16831000 DOI: 10.1021/la060592k] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/10/2023]
Abstract
Although atomic force microscopy (AFM) has emerged as the preeminent experimental tool for real-time in situ measurements of crystal growth processes in solution, relatively little is known about the mass transfer limitations that may impact these measurements. We present a continuum analysis of flow and mass transfer in an atomic force microscope fluid cell during crystal growth, using data acquired from calcium oxalate monohydrate (COM) crystal growth measurements as a comparison. Steady-state flows and solute concentration fields are computed using a three-dimensional, finite element method implemented on a parallel supercomputer. Steady-state flow results are compared with flow visualization experiments to validate the model. Computations of the flow field demonstrate how nonlinear momentum transport alters the spatial structure of the flow with increasing flow volume, altering mass transport conditions near the AFM cantilever and tip. The simulations demonstrate that the combination of solute depletion from crystal growth and mass transfer resistance lowers the solute concentration in the region between the tip and the crystal compared with the solute concentration at the inlet of the AFM cell. For example, using experimentally measured growth rates for COM, the solute concentration in this region is 3.1% lower than the inlet value because the solute consumed by crystal growth beneath the AFM tip cannot be replenished fully due to mass transport limitations. The simulations also reveal that increasing the flow rate through the cell does not affect this difference significantly because of the inherent shielding by the AFM tip in proximity with the crystal surface. Models such as the one presented here, used in conjunction with AFM measurements, promise more precise interpretations of measurement data.
Collapse
Affiliation(s)
- David Gasperino
- Department of Chemical Engineering and Materials Science, University of Minnesota, Minneapolis, MN 55455-0132, USA
| | | | | | | | | |
Collapse
|
8
|
Georgiou DK, Vekilov PG. A fast response mechanism for insulin storage in crystals may involve kink generation by association of 2D clusters. Proc Natl Acad Sci U S A 2006; 103:1681-6. [PMID: 16446456 PMCID: PMC1413625 DOI: 10.1073/pnas.0506526103] [Citation(s) in RCA: 35] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2005] [Indexed: 11/18/2022] Open
Abstract
Crystals that are likely rhombohedral of Zn-insulin hexamers form in the islets of Langerhans in the pancreases of many mammals. The suggested functions of crystal formation is to protect the insulin from proteases and increase the degree of conversion of soluble proinsulin. To accomplish these ends, crystal growth should be fast and adaptable to rate fluctuations in the conversion reaction. Zn-insulin crystals grow layer by layer. Each layer spreads by the attachment of molecules to kinks located at the layers' edges, also called steps. The kinks are thought to be generated either by thermal fluctuations, as postulated by Gibbs, or by 1D nucleation of new crystalline rows. The kink density determines the rate at which steps advance, and these two kink-generation mechanisms lead to weak near-linear responses of the growth rate to concentration variations. We demonstrate for the crystallization of Zn-insulin a mechanism of kink generation whereby 2D clusters of several insulin molecules preformed on the terraces between steps associate to the steps. This mechanism results in several-fold-higher kink density, a faster rate of crystallization, and a high sensitivity of the kinetics to small increases of the solute concentration. If the found mechanism operates during insulin crystallization in vivo, it could be a part of the biological regulation of insulin production and function. For other crystallizing materials in biological and nonbiological systems, this mechanism provides an understanding of the often seen nonlinear acceleration of the kinetics.
Collapse
Affiliation(s)
- Dimitra K. Georgiou
- Department of Chemical Engineering, University of Houston, Houston, TX 77204-4004
| | - Peter G. Vekilov
- Department of Chemical Engineering, University of Houston, Houston, TX 77204-4004
| |
Collapse
|
9
|
Reviakine I, Georgiou DK, Vekilov PG. Capillarity effects on crystallization kinetics: insulin. J Am Chem Soc 2003; 125:11684-93. [PMID: 13129373 DOI: 10.1021/ja030194t] [Citation(s) in RCA: 51] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
During layerwise growth of crystals, capillarity governs the generation of new crystal layers. Theory predicts that the line tension of the layer edge determines, via the characteristic two-dimensional capillary length L(c), the rates of generation and initial growth of the new layers. To test the correlation between L(c) and the rate of layer generation, we used in situ Tapping Mode Atomic Force Microscopy (TM-AFM) to study the generation and spreading of layers during crystallization of rhombohedral, R3, porcine insulin. We show that crystallization of this insulin form is uniquely suitable for such an investigation due to the linear kinetics of step growth it exhibits. This linear kinetics reflects the abundance of the incorporation sites along the rough steps, the lack of long-range step-step interactions, and the transport control of the growth kinetics. The kinetic coefficients are 7 x 10(-)(3) and 4 x 10(-)(2) cm s(-)(1), respectively, in the absence and presence of the cosolvent acetone-somewhat high for proteins and comparable to values for inorganic systems. We show that (i). the relevant capillary length, the size of a critical quadrangular 2D nucleus L(c), is the main scaling factor for the density of growth steps, while (ii). all steps longer than L(c) grow with a rate determined only by the supersaturation and independent of their length. We explain the divergence of (ii). from theoretical predictions with the high supersaturations typical of the growth of this protein system.
Collapse
Affiliation(s)
- Ilya Reviakine
- Department of Chemical Engineering, University of Houston, Houston, Texas 77204-4004, USA
| | | | | |
Collapse
|