Noninvasive Microcomputerized X-ray Tomography Visualization of

Jan 24, 2018 - Using PhotoShop CS5 (Adobe Systems, Inc., USA), we selected mineral deposits in each xy-plane frame using the color range selection too...
0 downloads 6 Views 32MB Size
Subscriber access provided by UNIV OF NEW ENGLAND ARMIDALE

Article

Non-invasive µCT visualization of mineralization directed by sea urchin- and nacre-specific proteins Martin Pendola, and John Spencer Evans Cryst. Growth Des., Just Accepted Manuscript • DOI: 10.1021/acs.cgd.7b01668 • Publication Date (Web): 24 Jan 2018 Downloaded from http://pubs.acs.org on January 30, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Crystal Growth & Design is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

Non-invasive µCT visualization of mineralization directed by sea urchin- and nacre-specific proteins. Martin Pendola and John Spencer Evans*

Laboratory for Chemical Physics, Center for Skeletal and Craniofacial Biology, New York University, 345 E. 24th Street, NY, NY, 10010 USA. *To whom correspondence should be addressed: John Spencer Evans, Laboratory for Chemical Physics, Division of Basic Sciences and Center for Skeletal and Craniofacial Medicine, New York University College of Dentistry, 345 E. 24th Street, New York, NY, 10010. Tel.: (212) 998-9605; Fax: (212) 995-4087. Email: [email protected].

ABBREVIATIONS rSpSM50 = recombinant version of the Strongylocentrotus purpuratus purple sea urchin spicule matrix protein 50; rPif97 = recombinant version of the Pinctada fucata Japanese pearl oyster shell framework nacre protein Pif97; µCT = microcomputerized X-ray tomography

1 ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ABSTRACT The biomineralization process offers novel principles for crystal engineering and solid-state chemistry, but to achieve this we must first understand how organisms such as the mollusk and sea urchin craft skeletal elements such as the shell and embryonic spicule, respectively. In vitro studies of mollusk- and sea urchin-associated proteins reveals that these proteins form hydrogel particles that control the nucleation process, assemble mineral nanoparticles, and modify the surfaces and interiors of existing crystals. However, visualization of these processes is hampered by destructive and invasive methods such as those used for conventional SEM/TEM. In this report we detail the novel use of µCT imaging to non-destructively investigate the in vitro calcium carbonate mineralization process in the presence of a recombinant sea urchin spicule matrix protein, rSpSM50, and a recombinant mollusk shell nacre protein, rPif97. Relative to the protein-free control scenario, both proteins generate calcite crystals that are consistent the results obtained from previous studies, but with µCT we discover new features; each protein generates a different number of mineral deposits, exhibits unique domains or motifs, and creates a stratification of mineral phases into layers. These results coincide with the known function of these proteins in vitro and in situ, and provide new information regarding biomineralization protein hydrogels and how they influence nucleation and crystal growth. We foresee that µCT imaging could, for appropriately sized systems, offer a bridge between in vitro experiments and in situ investigations of time-resolved nucleation and crystal growth phenomena in organisms, tissues, and synthetic materials.

2 ACS Paragon Plus Environment

Page 2 of 24

Page 3 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

INTRODUCTION The formation of calcium carbonate-based skeletal elements by oceanic organisms is a fascinating process that offers new and important insights into composite material formation. Two model organisms that have been extensively studied are the sea urchin1-9 and the mollusk.1020

In both cases, the mineralized skeletal elements that form follow a developmental plan that

involves the creation of a protein-based matrix into which nanoscale mineral deposits are formed and assembled into mesoscale single crystals of fracture-resistant calcite (sea urchin)1-8 or aragonite (mollusk).10-14 This process involves the formation of an amorphous precursor, amorphous calcium carbonate (ACC),2-5,16 which transforms into the appropriate crystalline polymorph (spicule = calcite; shell nacre = aragonite) during skeletal formation.

In both

organisms, there are protein families or proteomes6-9;17-20 that are expressed during this mineralization process that guide various aspects of single crystal formation, such as nanoparticle formation, ACC stabilization/transformation, and nanoparticle assembly. Thus, to better understand the formation of biological composites, we need to understand the role(s) that sea urchin and pearl oyster protein families play in mineralization. Previously, studies have shown that sea urchin spicule matrix21-23 and mollusk shell nacre layer24-26 proteins form ion-responsive “smart” mesoscale hydrogel particles that control the early events in in vitro non-classical nucleation processes27-29 and organize and assemble nanoscale mineral deposits within their matrices. This is an important event, in that subsequent assembly of mineral-impregnated hydrogel particles would position mineral nanoparticles to eventually form a crystal.30,31 Two protein hydrogelators which have been studied in this regard are the purple sea urchin (Strongylocentrotus purpuratus) embryonic spicule matrix protein, SpSM50,22 and the Japanese pearl oyster (Pinctada fucata) framework shell nacre layer protein, Pif97. Both 3 ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

proteins are strong hydrogelators; they form mesoscale protein hydrogel particles that contain mineral deposits.22,24 In addition, both protein hydrogel particles modify the surfaces of growing calcite crystals and form intracrystalline protein inclusions, thereby creating nanoporous crystals.22,24 In previous studies of both proteins, the organization and assembly of hydrogelmineral particles was analyzed using electron microscopic imaging (SEM, TEM),22,24 but because of sample dehydration requirements hydrogel collapse and distortion of mineral particle positioning occurred. Moreover, selective area imaging performed in these studies did not provide an overall picture of what occurs during the in vitro nucleation process. Thus, additional studies using appropriate non-destructive imaging methods and more inclusive sampling are necessary to advance our understanding of protein – directed mineral particle assembly and organization. In this present report we detail the use of microcomputerized X-ray tomography (microCT or µCT) to analyze and quantitate the in vitro calcium carbonate nucleation and crystal growth processes controlled by recombinant variants of SpSM50 (rSpSM50) and Pif97 (rPif97) proteins. Typically, µCT imaging has been used to examine a wide range of biological and material samples,32-38 but to our knowledge this is the first time that µCT has been applied to study and quantitate mineralization in aqueous solution. Using microvials, this non-invasive approach allows surveying of the entire mineralization volume and avoids dehydration distortions associated with past EM studies,22,24 thus preserving the positioning and organization of mineral deposits. We find that, compared to protein-free control scenarios, both rSpSM50 and rPif97 induce the formation of calcite mineral particles that are spatially distributed and exhibit multi-layer stratification, with rPif97 generating fewer crystals and rSpSM50 generating the largest number of mineral particles in this study. Further, motif or domain distributions for these 4 ACS Paragon Plus Environment

Page 4 of 24

Page 5 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

mineral particles differ from the control scenario, with rSpSM50 generating more motifs with smaller dimensions and rPif97 promoting mineral particle clustering leading to fewer motifs that are larger in size. In conclusion, we establish that both rSpSM50 and rPif97 hydrogel particles exhibit signature properties of mineral formation, distribution, and particle sizes.

These

differences may be important for the role that each protein plays in the formation of fractureresistant biological composites.

EXPERIMENTAL Protein preparation and in vitro mineralization assays. Recombinant S. purpuratus SpSM50 (rSpSM50) and P. fucata Pif97 (rPif97) proteins were subcloned and expressed in E. coli bacteria and purified as described in previous studies.22,24 Purified proteins were dissolved in 30 nm ultrapure Molecular Biology grade water (Fisher Scientific, USA) to create stock solutions for assay use. In vitro biomineralization assays (10 mM CaCl2 / NaHCO3 / Na2CO3, ~ pH 8.0, 60 min duration, 25 oC) were performed as described previously,22,24 using a final assay volume of 200 µL and hinged-lid sealable plastic polypropylene lab microvials (outside dimensions: 12.7mm diameter, 13.03mm height; inside dimensions: 7.37mm diameter, 5.72mm height)(LAContainer Store, Yorba Linda, CA, USA)(Figure 1). Prior to performing experiments with the microvials we obtained phantom images of the vials containing 200 µL of ultrapure 30nm filtered Molecular Biology grade water (Fisher Scientific, USA), which show typical water interference spots that occur in the upper z-axis region of the vials and not at the bottom of the vials (Figure S1, Supporting Information). Subsequently, three sample scenarios were created: a control scenario (no protein was added), and two protein scenarios using aliquots of rSpSM50 and rPif97 stock solutions (final protein concentrations = 1.5 µM). In parallel assays a 5mm x 5mm silicon wafer (Ted Pella, USA) was placed at the bottom of each vial to collect mineral and 5 ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 24

hydrogel deposits that formed by the 60 min period for subsequent microRaman and SEM analyses.

Figure 1. Snap-cap polypropylene vials utilized in the MicroCT study. Note the round region at the center of the vial bottom exterior (arrow); this is a manufacturing “stub” on the outside of the vial that remains from the sprue injection manufacturing process. This stub protrudes approximately 0.5mm into the interior of the vial as a hemispherical region. Outside dimensions: 12.7mm diameter, 13.03mm height; Inside dimensions: 7.37mm diameter, 5.72mm height. The vials hold approximately 200 µL of liquid.

µCT imaging and data processing.

At the end of 60 min, each mineralization assay vial was

mounted to the stage and scanned using a Bruker MicroCT SkyScan High Resolution Model 1172G (Bruker Scientific, Kontich, Belgium) with a 59kV source voltage and 167 µA current. The scanning process involved a round scanning trajectory and step-and-shoot protocol, with each image requiring 0.3 degrees/step and 780 msec exposure time. Approximately 1300 images with a final resolution of 4.99 µm per pixel were acquired for each sample over an 80 min period; thus, the total mineralization assay elapsed time was 140 minutes. For each sample the sequence of images was exported as a series of bmp image files. For every sample we identified the first xy-plane or transverse image showing mineral deposits at the bottom of the vial, then progressing in the +z-axis direction we selected the next 189 transverse images to complete a sequence of 190 images from the bottom of the vial up to approximately 950 µm in height. Collectively, these images encompassed the mineral deposits in toto generated in each assay. Using PhotoShop CS5 (Adobe Systems, Inc, USA) we selected mineral deposits in each xy-plane frame using the color range selection tool, with a low fuzziness 6 ACS Paragon Plus Environment

Page 7 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

factor to prevent selection of false positive objects.

Once the objects were selected, we

subtracted all the other elements in the image using the inverse selection tool to retain only mineral deposits. This process was repeated for all 190 images. The 190 images that were produced were then used for 3D reconstruction. These image stacks were imported into ImageJ 1.46r 64bit (NIH, USA) and rendered into 3D using the volume viewer tool with tricubic (sharp) interpolation. From this 3D reconstruction we selected 2 transverse images with a 0.25% (low magnification) and 50% (high magnification) scales. For further quantitative analyses of mineral deposits, we used the volume viewer reconstruction and analyze commands in ImageJ, and counted the number of mineral particles and measured their dimensions in each sample. In addition, the volume viewer images produced with ImageJ were analyzed using MountainsMap 7 (Digital Surf, France). Using the tool Binary segmentation, the events present in each image were grouped based on shapes respect the background, and represented as a single motif, with different colors assigned for each motif. MicroRaman and SEM analyses of mineral deposits. Using the deposits which collected onto Si wafers at the bottom of each assay vial, all samples were washed 3x with calcium carbonate – saturated methanol and dried overnight.22,24 Subsequently, these Si wafers were analyzed using micro-Raman spectroscopy (Thermo Scientific DXR Raman microscope in dark field mode, 532nm laser wavelength, 6mW power, 1800/900 grating, 25µM pinhole aperture, 16 exposures per sample at MP50x magnification. The spectra were processed using the OMNIC Atlus software designed for the instrument.22,24 Other parallel wafers were imaged using a Merlin (Carl Zeiss) field emission SEM (FESEM) with an Everhart-Thornley type secondary electron detector (SE2) at an accelerating voltage of 2kV and a probe current of 100 pA.22,24 Prior to analysis, samples were coated with iridium (3nm layer) using a Cressington 208HR sputter coater with

7 ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

thickness controller.22,24

RESULTS The purpose of this µCT imaging study is to non-invasively detect and analyze x-ray attenuating calcium carbonate mineral particles that form in the presence and absence of the two hydrogelator proteins within a polypropylene microvial system. The mineralization assay used in this study is tuned to stable calcite production.22-26 Our goal is to distinguish the quantity, size, distribution, and organization of mineral particles relative to each other and to a protein-free

scenario. Figure 2. (A) µCT images (along z-axis) of vials containing mineral deposits formed in control [(-) protein] and rSpSM50, rPif97 mineralization assays. Arrows denote presence of mineral particle clustering. Scalebars = 1mm. (B) 3D reconstruction of the same samples. The raised hemispherical region that represents the production spar can be clearly seen internally at the bottom of each vial. Orientation xyz axes shown.

The hydrogels themselves possess scarce x-ray attenuation features and thus due to low imaging contrast will not be directly detected by µCT. However, the influence of the hydrogels on mineral formation (dimensions, location, distribution) will be indirectly detected, as shown below. The use of sealable polypropylene plastic vials minimizes evaporation and carbon dioxide 8 ACS Paragon Plus Environment

Page 8 of 24

Page 9 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

outgassing from carbonate/bicarbonate species39 during the assay and subsequent imaging period and minimizes the backscattering x-ray beam artifacts (Figure S1, Supporting Information) associated with glass containers.

Moreover, given the dimensions of the vials, we can

completely image and quantitate the entire mineralization assay volume, something not possible in previous studies.22-26 Spicule and nacre-associated hydrogel particles differ in mineral nucleation. In Figure 2, we present the z-axis projection image (Figure 2A) and 3D reconstruction (Figure 2B) images of the liquid volumes contained within each sample vial with the surrounding plastic vial surfaces omitted. Two features are noted in this series of images. 1) Both images reveal the total population of mineral particles generated by each sample. We tabulated the total number of particles detected in each sample (Figure 3A) and the average area occupied by all the particles in each sample (Figure 3B). From these tabulations we find that the number of detected crystals and occupied area follows the relationship: rSpSM50 > control > rPif97. Both proteins are structurally similar and are nearly equal in aggregation propensity22,24 and thus these factors probably do not contribute significantly to these protein-specific results. However, it is known that Pif97 inhibits calcite crystal growth,18,24 and thus the lower number of mineral particles formed in the presence of rPif97 is logical. Likewise, rSpSM50, which forms a transient vateritic ACC phase in vitro that quickly converts to a more stable calcitic ACC phase, is expected to promote calcite formation.1,4,5,22 Secondly, we detected evidence of mineral particle clustering, i.e., the presence of particles which appear to be either in contact or closely spaced with one another (Figure 2A, arrows). We observed that the occurrence of clustering follows the relationship: rPif97 > control > rSpSM50, i.e., the inverse of the particle number relationship. In the case of the control, the presence of clustering may be attributable to epitaxial growth.

9 ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

However, for rPif97, given the millimeter dimensions of the formed clusters, it is unlikely that epitaxial growth contributed to cluster formation. More likely, the long-range clustering that we observe in the presence of rPif97 reflects protein hydrogel adhesion to mineral particles that leads to clustering during mineralization.24 From the 3D reconstruction and z-axis transverse images we conclude that relative to the control scenario rSpSM50 generates a larger number of mineral particles, whereas rPif97 generates fewer mineral particles (i.e., acts as a calcite inhibitor) over the same time interval but promotes particle clustering. Sagittal imaging plane reveals mineralized stratification. The 3D reconstruction images in Figure 2B suggest that in addition to mineral precipitation at the bottom of each vial, there exist mineral particles which are suspended in the assay liquid above the vial bottoms. To assess this, we

Figure 3. (A) Total mineral particle number and (B) average mineral particle area analyses for the three mineralization scenarios. For (B), one-way ANOVA (analysis of variance) shows that the differences between samples is significant (p > 0.0001).

examined cross-sectional xz-plane sagittal images taken at the midline of each vial (Figure 4; Figure S2, Supporting Information). In all samples we note the presence of mineral particles that are suspended in the liquid above the vial bottoms. Presumably, these are mineral particles whose masses or densities are lower than those which reside at the bottom of the vials. Thus, these particles were unable to settle out during the assay and imaging time periods. Of the three 10 ACS Paragon Plus Environment

Page 10 of 24

Page 11 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

samples, it appears that rPif97 possesses a larger number of suspended mineralized particles compared to either the control or rSpSM50 samples, and features evidence of cluster formation (white arrow, Figure 4). This is consistent with rPif97 hydrogel particles creating large micronsized protein inclusions within calcite during these assays,24 which would make these crystals more nanoporous and potentially more buoyant in aqueous solution. But perhaps the most intriguing feature noted in these sagittal images is that both the rPif97 and rSpSM50 samples form stratified mineralized phases that settle out over time (denoted as Layers 1, 2) and are distinguishable by their degree of x-ray attenuation. In comparison, the control sample possesses only a single layer, and thus the

11 ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 4. µCT XZ-sagittal plane volumetric sections (4mm x 4mm x 1.9mm) of microvial mineralization samples. In all samples, note the presence of suspended mineral particles above the bottom of the vials. In the case of the protein samples, note the presence of two layers of mineral particles (labeled as 1,2). Note also the presence of a large, suspended mineral cluster in the rPif97 image (white arrow). Scalebar = 1mm; height of each image = 1.9mm.

stratification effect is protein-driven. In the case of rPif97, Layer 1 appears to possess less x-ray attenuation compared to Layer 2, and interestingly, the reverse is true for the rSpSM50 sample. Note that the dimensions of the rSpSM50 layers are greater than those generated by either rPif97 or the control scenario, and that the rSpSM50 Layer 2 is significantly thicker and less dense than Layer 1. We believe that this stratification phenomenon originates from one or more of the following phenomena: 1) Protein differences in mineral particle organization, crystallinity (i.e., calcite crystals versus stabilized ACC particles), size, and assembly22,24 induced by the hydrogel

12 ACS Paragon Plus Environment

Page 12 of 24

Page 13 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

phases that lead to stratification differences; 2) Differences in protein content associated with the mineral particles, either through external surface associations or internal inclusions (i.e., intracrystalline protein-containing nanoporosities),2,3,14,15,22,24 which would alter the buoyant densities, sedimentation, and spatial positioning of the mineral phases, leading to the proteinspecific stratification effects. In conclusion, µCT imaging reveals a new in vitro protein-driven biomineralization phenomenon:

nacre and sea urchin matrix proteins generate a mineral

stratification effect that is protein-specific. We believe that this layering phenomenon reflects the influence of protein hydrogel particles on the calcium carbonate nucleation, transformation, stabilization, assembly, and crystal growth processes.22,24 Analyses of mineral particles and hydrogel phases generated by rSpSM50 and rPif97.

Using

image processing we visualized the mineral particles that were generated in the control, rSpSM50, and rPif97 assays (Figure 5A). Given that µCT resolution is approximately 5 µm, and that most mineral particles that form in the presence of both proteins are > 5 µm,22,24 we were able to detect both individual and clustered mineral particles formed in all assays. Obviously, aliasing becomes a problem at higher magnification, and thus we also performed SEM imaging on mineral deposits collected on Si wafers in parallel assays to determine crystal morphologies more precisely (Figure 5B). A comparison of the µCT and SEM images for the control and rSpSM50 scenarios reveal overall agreement:

we observe the typical rhombohedral

morphologies associated with calcite in these samples (Figure 5A,B). Moreover, as shown by the arrows in the SEM images, we also confirm the formation of protein phases or hydrogels depositing onto the Si wafers or onto the crystals themselves (Figure 5B). Using microRaman spectroscopy, the calcite polymorph was also confirmed for the modified mineral crystals

13 ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 24

generated by rPif97 samples (Figure S3, Table S1, Supporting Information). Note that we cannot confirm nor rule out the presence of ACC deposits

Figure 5. (A) Enlargement of representative µCT images of mineral particles that deposited in (-) protein, rSpSM50, and rPif97 mineralization assays. Scalebars = 4 microns. (B) Corresponding representative SEM images of mineral deposits collected from parallel mineralization assays. Note the close correspondence in morphologies between MicroCT and SEM images of crystals in each sample, and the presence of protein aggregates or hydrogel phases that have deposited onto the Si wafers or crystals (arrows). In SEM images scalebars correspond to either 2 or 1 microns.

in these assays. Hypothetically, any stabilized ACC that formed in the assay may be suspended in solution and probably associated with the protein hydrogel phase (Figure 4).

We also note

that for each protein hydrogel, signature crystal morphologies were observed that correlate with previous assay results (Figure 5A,B): The rounded, highly modified calcite crystals of rPif9724 are easily distinguishable from the smaller, rhombohedral crystals of rSpSM5022 using µCT. However, as one might expect, the nanotextured surface modifications created by rSpSM50 that are readily detected by SEM22 cannot be detected by µCT.

Thus, the microvial assays

qualitatively replicate the results obtained in past in vitro mineralization studies for each protein hydrogel22,24 and µCT imaging affirms that each protein differentially regulates calcite crystal

14 ACS Paragon Plus Environment

Page 15 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

size and morphology, but cannot provide information regarding protein-induced surface or subsurface modifications of crystals. Protein hydrogel-mediated organization of mineral particles. We now examine the arrangement or organization of mineral particles which form during the mineralization assay period. This can be done using binary segmentation analysis (MountainsMap software package), where mineral particles are selected by shape and then groups, “domains”, or “motifs” are defined and limited by the highest percentage of mineral particles that cover the smallest possible surface area.

Domain Figure 6: Domain mapping of microvials along the z-axis. Red is considered as empty or low density areas where mineral deposits have no relevant proximity to each other, other colors represent distinct domain types (colors were arbitrarily assigned). Each axis represents 1.7 cm in length. The number of mineral particle-containing domains are provided below each image.

boundaries are denoted by different colors, with red assigned to areas of lowest particle density or where no crystals exist (Figure 6). Using the control scenario as a reference (2969 domains) we can observe that rSpSM50-associated mineral particles occupy smaller but more numerous domains (4966, nearly 2x the number of the control scenario). Note that smaller mineral particle sizes have been reported to be efficient in space-filling arrangements within sea urchin skeletal elements,40 so the smaller, numerous domains generated by rSpSM50 are qualitatively consistent

15 ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

with these in situ findings. Conversely, the mineral particles generated in the presence of rPif97 exist in fewer domains (2264) than either the control or rSpSM50 samples but clearly possess larger domain regions and feature significant open regions where little or no mineral particles occur (Figure 6). We note that the domain phenomena (Figure 6) correlates with clustering phenomena, i.e., rPif97 > control > rSpSM50 (Figure 2). We presume that the protein hydrogel particles influence the formation of mineral particle domains in terms of number and size and this may have relevance for the formation of the crystalline phase from amorphous mineral nanoparticle precursors in both the nacre layer and the spicule matrix.

DISCUSSION µCT has been applied to the study of many different biological and material samples over the last 10 years, and has established itself as a valuable visualization tool for solid and liquid samples.32-38 In our present study, we adapted µCT for a unique purpose: a non-invasive study of real-time protein-driven nucleation of calcium carbonate mineral particles in an aqueous environment. As we have discussed above, this approach allowed us to completely image the entire mineralization assay volume. In doing so, the method yielded some very interesting information regarding the in vitro nucleation and crystal growth functionality of a recombinant mollusk shell nacre protein, rPif97,24 and a recombinant embryonic sea urchin spicule matrix protein, rSpSM50.22 Both proteins are strong hydrogelators associated with mesoscale crystal formation from the assembly of mineral nanoparticle precursors,1-5;10-14,16 and thus make excellent model systems for this type of study. We find that each protein has signature mineral particle formation and organizational characteristics. Specifically, rSpSM50, which oversees the transformation of efficient space-filling ACC nanoparticles40 to single crystal calcite in the developing spicule,1,4,5 generates the largest number of crystals with the smallest dimensions that 16 ACS Paragon Plus Environment

Page 16 of 24

Page 17 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

cover the largest area of the vials (Figures 2,3), produces the largest number of discrete crystal domains (Figure 6), and creates stratified layers of mineralized phases in the assay (Figure 4). In contrast, rPif97, a known calcite inhibitor,18,24 creates fewer mineral particles (Figures 2,3), but fosters larger mineral domains (Figure 6) with evidence of crystal clustering (Figures 2,4). This protein also creates stratified mineralized phase layers (Figure 4), but, compared to rSpSM50, this stratification effect is reduced in dimension and x-ray attenuation is not as pronounced compared to rSpSM50 (Figure 4). Thus, µCT imaging has detected new phenomena regarding rPif97- and rSpSM50-directed in vitro calcium carbonate nucleation and crystal growth that were previously undetected using invasive or destructive imaging methods.22,24 We believe that the in vitro mineralization behavior of rPif97 and rSpSM50 protein hydrogels provides new insights vis a vis the involvement of native Pif97 and SpSM50 in the mineralization processes that occur within the mollusk shell nacre layer and the sea urchin spicule, respectively. It is known that calcium carbonates such as calcite form from ACC nanoparticle assembly in vitro,30,31,41 and given that both skeletal elements arise from nanoparticle formation and assembly it is interesting to note that both proteins generate mineral particles of a specific quantity, size, and area (Figures 2,3). Given that SpSM50 is the most abundant matrix protein in the embryonic spicule,6-9 it may be the responsibility of this protein hydrogel to promote nucleation and limit the size of the mineral deposits (Figures 2-4) to facilitate space-filling arrangements40 and promote the assembly of single crystal calcite. Conversely, it is known that Pif97 acts as a calcite inhibitor and exists as a complex with Pif80 and n16.3 in the nacre layer of the oyster shell,17,18 and thus the formation of fewer calcite crystals, induction of crystal clustering, and creation of larger mineral domains may reflect Pif97 participation within this ternary protein complex. The other interesting feature is the stratification

17 ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

of mineral particles in each protein assay (Figure 4). This layering phenomenon suggests that the protein-mediated formation of mineral phases over time is not a homogeneous process, but a heterogeneous one that produces different particle dimensions and densities that reflect protein hydrogel content and/or association with the mineral phase. At this time we do not understand the significance of this in vitro phenomenon and its relevance to in situ mineralization events, and thus additional research into this process will be required in order to establish its importance. Finally, the clear strength of the µCT approach is that the systems under study remain nearly undisturbed and the instrumentation permits the non-invasive study and quantitation of nucleation and crystal growth. Potentially µCT could allow for the comparisons of in situ nucleation processes alongside parallel in vitro processes in a wide range of synthetic systems, materials, or organisms and tissues whose dimensions are compatible with µCT. This approach can also allow time-resolved monitoring or kinetics of nucleation processes, albeit at lower resolution compared to EM methods. This is a particularly attractive and powerful “bridging” scenario for linking crystal growth events that occur in vitro and in situ as a function of time, whether directed by proteins or other additives such as polymers. Finally, using appropriate contrasting agents,34,37 it may be possible to visualize the organic or protein-based hydrogel within in vitro and in situ environments, thereby gaining additional insight into hydrogel-directed mineral particle formation, assembly, and organization. We intend to pursue these avenues of research in future applications. Author Contributions The manuscript was written through contributions of all authors. All authors have given approval to the final version of the manuscript.

18 ACS Paragon Plus Environment

Page 18 of 24

Page 19 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

Funding Sources This research was supported by the Life Sciences Division, U.S. Army Research Office, under award W911NF-16-1-0262. ACKNOWLEDGMENT This report represents contribution number 90 from the Laboratory for Chemical Physics, New York University. ASSOCIATED CONTENT Supporting Information. 3D and XY-transverse images of microvials containing 200 µL of 30nm filtered Molecular Biology grade water, 25 oC (water phantom)(Figure S1); Sagittal XZ plane images obtained from control, rPif97, and rSpSM50 mineralization assays (Figure S2); MicroRaman spectra of rPif97 mineralization assay mineral deposits (Figure S3) and Raman band assignments for calcium carbonates (Table S1). This material is available free of charge via the Internet at http://pubs.acs.org.

References 1 Seto, J., Ma, Y., Davis, S.A., Meldrum, F., Schilde, U., Gourrier, A., Jaeger, C., Cölfen, H., Structure-property relationships of a biological mesocrystal in the adult sea urchin spine. Proc. Natl. Acad. Sci USA 2012, 109, 3699-3704. 2 Berman, A., Addadi, L., Kvick, A., Leiserowitz, L., Nelson, M., Weiner, S. Intercalation of sea urchin proteins in calcite: Study of a crystalline composite material. Science 1990, 250, 664667. 3 Aizenberg, J., Hanson, J., Koetzle, T.F., Weiner, S., Addadi, L. Control of macromolecular distribution within synthetic and biogenic single calcite crystals. J. Am. Chem. Soc. 1997, 119, 881-886. 4 Tester, C.C., Wu, C.-H., Krejci, M. R., Mueller, L., Park, A., Lai, B., Chen, S., Sun, C., Joester, D. Time-resolved evolution of short- and long-range order during the transformation of ACC to calcite in the sea urchin embryo. Adv. Fun. Mat. 2013, 23, 4185-4194. 5 Politi, Y, Metzler, R.A., Abrecht, M., Gilbert, B., Wilt F.H., Sagi, I, Addadi, L., Weiner, S., Gilbert, P. Transformation mechanism of amorphous calcium carbonate into calcite in the sea urchin larval spicule. Proc Natl Acad Sci U S A. 2008,105, 17362-17366. 6 Sea urchin genome sequencing consortium et. al., The genome of the sea urchin Strongylocentrotus purpuratus. Science 2006, 14, 941-952. 7 Mann, K., Poustka, A.J., Mann, M., The sea urchin (Strongylocentrotus purpuratus) test and spine proteomes. Proteome Science 2008, 6, 1-13. 8 Mann, K., Wilt, F.H., Poustka, A.J. Proteomic analysis of sea urchin (Strongylocentrotus purpuratus) spicule matrix. Proteome Science 2010, 8, 1-14.

19 ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

9 Livingston, B.T., Killian, C.E., Wilt, F., Cameron, A., Landrum, M.J., Ermolaeva, O., Sapojnikov, V., Maglott, D.R., Buchanan, A.M., Ettensohn, C.A. Genome-wide analysis of biomineralization-related proteins in the sea urchin Strongylocentrotus purpuratus. Devel. Biol. 2006, 300, 335-348. 10 Zhang, G., Li, X. Uncovering aragonite nanoparticle self-assembly in nacre – A natural armor. Cryst. Growth Design 2012, 12, 4306-4310. 11 Li, X., Chang, W.C., Chao, Y.J., Wang, R., Chang, M. Nanoscale structural and mechanical characterization of a natural nanocomposite material: The shell of red abalone. NanoLetters 2004, 4, 613-617. 12 Sun, J., Bhushan, B. Hierarchical structure and mechanical properties of nacre: A review. RSC Adv. 2012, 2, 7617-7632. 13 Li, X., Huang, Z. Unveiling the formation mechanism of pseudo-single-crystal aragonite platelets in nacre. Phys. Rev. Lett. 2009, 102, 075502-075506. 14 Zheng, G., Xu, J. From colloidal nanoparticles to a single crystal: New insights into the formation of nacre's aragonite tablets. J. Struct. Biol. 2013, 182, 36-43. 15 Bezares, J., Asaro, R.J., Hawley, M. Macromolecular structure of the organic framework of nacre in Haliotis rufescens: Implications for mechanical response. J. Struct. Biol. 2010, 170, 484-500. 16 DeVol, R.T., Sun, C.Y., Marcus, M.A., Coppersmith, S.N., Myneni, S.C.B., Gilbert, P.U.P.A. Nanoscale transforming mineral phases in fresh nacre. J. Am. Chem. Soc. 2015, 137, 1332513333. 17 Suzuki, M., Iwashima, A., Kimura, M., Kogure, T., Nagasawa, H. (2013) The molecular evolution of the Pif family proteins in various species of mollusks. Mar Biotech. 15, 145-158. 18 Suzuki, M., Saruwatari, K., Kogure, T., Yamamoto, Y., Nishimura, T., Kato, T., Nagasawa, H. (2009) An acidic matrix protein, Pif, is a key macromolecule for nacre formation. Science 325, 1388-1390. 19 Immel, F., Gaspard, D., Marie, A., Guichard, N., Cusack, M., Marin, F. (2015) Shell proteome of rhynchonelliform brachiopods. J. Struct. Biol. 190, 360-366. 20 Zhang, G., et al., (2012) The oyster genome reveals stress adaptation and complexity of shell formation. Nature 490, 49-54. 21 Jain, G., Pendola, M., Rao, A., Cölfen, H., Evans, J.S., A model sea urchin spicule matrix protein self-associates to form mineral-modifying protein hydrogels. Biochemistry 2016, 55, 4410-4421. 22 Jain, G., Pendola, M., Huang, Y.C., Gebauer, D., Evans, J.S. A model sea urchin spicule matrix protein, rSpSM50, is a hydrogelator that modifies and organizes the mineralization process. Biochemistry 2017, 56, 2663-2675. 20 ACS Paragon Plus Environment

Page 20 of 24

Page 21 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

23 Pendola, M., Davidyants, A., Jung, Y.S., Evans, J.S. (2017) Sea urchin spicule matrix proteins form mesoscale hydrogels that exhibit selective ion interactions. ACS Omega, 2, 61516158. 24 Chang, E.P., Evans, J.S. (2015) Pif97, a von Willebrand and Peritrophin biomineralization protein, organizes mineral nanoparticles and creates intracrystalline nanochambers. Biochemistry 54, 5348-5355. 25 Jain, G., Pendola, M., Huang, Y.C., Colas, J.J., Gebauer, D., Johnson, S., Evans, J.S. (2017) Functional prioritization and hydrogel regulation phenomena created by a combinatorial pearlassociated 2-protein biomineralization model system. Biochemistry 56, 3607-3618. 26 Perovic, I., Davidyants, A., Evans, J.S. (2016) Aragonite-associated mollusk shell protein aggregates to form mesoscale “smart” hydrogels. ACS Omega 1: 886-896. 27 Kellermeier, M., Raiteri, P., Berg, J.K., Kempter, A., Gale, J.D., Gebauer, D. Entropy drives calcium carbonate ion association. ChemPhysChem 2016, 17, 3535-3541. 28 44 Gebauer, D., Volkel, A., Cölfen, H., Stable prenucleation of calcium carbonate clusters. Science 2008, 322, 1819-1822. 29 Demichelis, R., Raiteri, P., Gale, J.D., Quigley, D., Gebauer, D., Stable prenucleation mineral clusters are liquid-like ionic polymers. Nature Commun. 2011, 2, 1-8. 30 De Yoreo, J.J., Gilbert, P.U.P.A., Sommerdijk, N.A.J.M., Penn, R.L., Whitelam, S., Joester, D. Zhang, H., Rimer, J.D., Navrotsky, A., Banfield, J.F. Wallace, A.F., Michel, F.M. Meldrum, F.C., Cölfen, H., Dove, P.M., Crystallization by particle attachment in synthetic, biogenic, and geologic environments. Science 2015, 349, 498-505. 31 Wallace, A.F.,Hedges, L.O., Fernandez-Martinez, A. Raiteri, P., Gale, J.D., Waychunas, G.A., Whitelam, S., Banfield, J.F., De Yoreo, J.J. Liquid-liquid separation in supersaturated calcium carbonate solutions. Science 2013, 341, 885-889. 32 Avanzi, F., Petrucci, G., Matzl, M., Scheebeli, M., De Michele C. Early formation of preferential flow in a homogeneous snowpack observed by micro-CT. AGU Water Resources Research 2017, 3, 3713-3729. 33 Islam, S.F., Sundara, R.V., Whitehouse, S., Althaus, T.O., Palzer, S., Hounslow, M.J., Salman, A.D. Movement of a secondary immiscible liquid in a suspension using a non-invasive technique. Chem. Eng. Res. Design 2016, 110, 160-170. 34 Crica, L.E., Wengenroth, J., Tiainen, H., Ionita, M., Haugen, H.J. Enhanced X-ray absorption for micro-CT analysis of low density polymers. J. Biomat. Sci. Polymer Ed. 2016, 27, 805-823. 35 Chiesura, G., Luykx, G., Voet, E., Lammens, N., Van Paepegem, W., Degrieck, J., Dierick, M., Van Hoorebeke, L., Vanderniepen, P., Sulejmani, S., Sonnenfield, C., Geernaert, T., Berghmans, F. A micro-computed tomography technique to study the quality of fibre optics embedded in composite materials. Sensors 2015, 15, 10852-10871

21 ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

36 Du Plessis, A., le Roux, S.G., Els, J., Booysen, G., Blaine, D.C. Application of microCT to the non-destructive testing of an additive manufactured titanium component. Case Studies in Non-destructive Testing and Evaluation 2015, 4, 1-17. 37 De S. e. Silva, J.M., Zanette, I., Noel, P.B., Cardoso, M.B., Kimm, M.A., Pfeffer, F. Threedimensional non-destructive soft-tissue visualization with X-ray staining micro-tomography. Nature Sci Rep. 2015, 5, 14088-14092. 38 Swart, P., Wicklein, M., Sykesm D., Ahmed, F., Krapp, H.G. A quantitative comparison of micro-CT preparations in Dipteran flies. Nature Sci. Reports 2016, 6, 39180-39186. 39 Mutch, G.A, Morandi, S., Walker, R., Anderson, J.A., Vega-Maza, D., Operti, L., Cerrato, G. Cation-dependent carbonate speciation and the effect of water. J. Phys. Chem. 2016, 120, 17570-17578. 40 Yang, L., Killian, C.E., Kunz, M., Tamura, N., Gilbert, P.U.P.A. Biomineral nanoparticles are space-filling. Nanoscale 2011, 3, 603-609. 41 Navarro, C.R., Cara, A.B., Elert, K., Putnis, C.V., Ruiz-Agudo, E. Direct nanoscale imaging reveals the growth of calcite crystals via amorphous nanoparticles. Crystal Growth and Design 2016, 16, 1850-1860.

22 ACS Paragon Plus Environment

Page 22 of 24

Page 23 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

For Table of Contents Use Only Non-invasive µCT visualization of mineralization directed by sea urchin- and nacre-specific proteins

Martin Pendola and John Spencer Evans Synopsis: Non-invasive µCT visualization of mineralization directed by sea urchin- and nacrespecific proteins reveals protein-specific mineralization phenomena (particle size, density, distribution) in vitro.

23 ACS Paragon Plus Environment

A

1 2 3 4 5 6 7 8 9 10

(-) Protein

Crystal Growth & Design

y x 1mm

B ACS Paragon Plus Environment z

y x

rSpSM50

rPif97

Page 24 of 24