2.1 Silk fibroin continuously self-assembles into nanoglobular coatings in the presence of potassium phosphate. Snapshots of coating formation over time were generated by immersing a TiO2 substrate, chosen for its common use in medical implants, in an aqueous solution containing dilute silk fibroin (0.5 mg/ml) and potassium phosphate (200mM, pH 5.0). Figure 1A shows monotonically increasing coating thickness, reaching ~ 20 nm (dry), over the 24 hr observation period. Dry height profile AFM images at discrete time points in the coating process (Fig. 1B) show that coatings are comprised of globular aggregates which continuously accumulate on the substrate surface over time. Individual globules at early timepoints (i.e. 15 min) exhibit slightly larger apparent diameters of 60 ± 7 nm compared to 50 ± 7 nm for multilayer coatings made for 24 hrs likely due to spreading of the globule structures on the substrate surface. Globular size evolution over 24 hr coating growth is provided in (Figure S1).
The relatively narrow size distribution of globular aggregates is surprising considering the polydispersity of the fibroin molecular weight due to the degumming process33 and the natural tendency of silk proteins to assemble into fibrillar structures through a nucleation and growth mechanism18,22,26,34–39. At early coating timepoints (15 min – 1 hr), portions of the underlying TiO2 surface are visible while silk fibroin nanoglobules appear to be randomly distributed across the surface. Analysis of AFM images yields surface coverage of 30 ± 5%, 47 ± 4%, and 63 ± 3% for coatings made for 15 minutes, 30 minutes, and 1 hour, respectively. After 3 hours, the underlying substrate becomes no longer observable as globule density increases past full surface coverage. As growth continues, the coating surface remains smooth, globular, and uniform while thickness increases. At long coating timepoints (i.e. 24 hours) the coatings exhibit extremely smooth surfaces, with calculated surface roughness of ~ 2.8 nm indicating that the continuous deposition of the silk fibroin globules occurs homogeneously over the substrate surface. This result suggests that under these solution conditions the interfacial self-assembly process is well-controlled, without autocatalytic formation of precipitates or large mounds. The high surface density of silk fibroin is atypical for protein adsorption, which are generally do not reach full coverage (e.g. random sequential adsorption model predicts saturation at ~ 55% surface coverage)10. The ability of silk fibroin to form strong supramolecular interactions in the coating buffer is likely responsible for its high packing density.
Our AFM images suggest that our coatings form though accumulation of globular nanoaggregates onto the surface, rat her than nanofibrils previously observed with recombinant spidroins, such as eADF4(C16)18 and 4Rep-CT25. Analysis of isolated nanoglobules at the 15 minute timepoint allows for an estimation of the number of silk fibroin chains comprising each globule. Analyzing 100 individual nanoglobules using the Asylum AFM software, we calculate their average volume to be 1.31x104 ± 5.35x103 nm3. Based on an average degummed silk fibroin molecular weight of 100 kDa and a protein density of 1.0–1.3 g/cm340, each nanoglobule is estimated to be comprised of 80–100 individual silk fibroin chains. The silk fibroin coatings swell when imaged in ultrapure water (Fig. 1C), suggesting that these coatings are hydrogel-like, which is consistent with previous reports of silk fibroin-based materials41. In the hydrated state, nanoglobules can be observed to be comprised of smaller sub-globules, which supports their hierarchically assembled nature. It is difficult to estimate the number of silk fibroin molecules present in these sub-globules due to the difficulty of isolating individual sub-globules from the larger globule complex. However, assuming the subglobules exhibit a half ellipsoid-like configuration on the TiO2 surface with an average diameter of ~ 30 nm and average height of ~ 12 nm via AFM analysis of 5 distinct subglobules, their estimated volume is ~ 5.65 x 103 nm3. Assuming an estimated hydrated protein density of 0.74–0.97 g/cm3 (assuming a swelling ratio of 1.34 by comparing the dry-state globule volume via AFM to the hydrated-state globule volume via DLS) this would suggest that the subglobules are comprised of 25–33 individual silk fibroin chains.
2.2 Coating growth relies on a balance of protein-protein and protein-surface interactions. The rate of coating formation is highly dependent on both the phosphate concentration and solution pH of the coating solution, as evidenced by coating thickness after 24 hrs on TiO2 using different coating solution compositions (Fig. 2A). Coatings could not be grown above pH 6.5 at any phosphate concentration, nor above a phosphate concentration of 300 mM at any solution pH. Even within a pH range of 5–6.5 and phosphate range of 100 mM – 300 mM, a regime where coating growth is facilitated, there are stark differences in coating thickness. A solution containing 200 mM phosphate at pH 5 appears to optimally facilitate coating growth. At this pH, silk fibroin is near its isoelectric point of 4.3942 and thus has decreased protein-protein charge repulsion, as confirmed by V-potential analysis (Figure S2B). Bare TiO2 has a slight negative charge, so protein-surface repulsion is also decreased at pH 5, thus further facilitating adsorption. However, better coating growth at intermediate phosphate concentrations is counterintuitive, as phosphate is a kosmotropic ion that promotes protein-protein interactions and chain folding by removing the water hydration layer around protein backbones. This result suggests that the interfacial assembly of nano-thin silk fibroin coatings relies on a careful balance of protein-protein and protein-surface interactions.
Solution-phase self-assembly occurs concurrently with interfacial coating growth in our system, thus we investigated the solution-phase aggregation behavior of silk fibroin under different phosphate concentrations at pH 5 via dynamic light scattering (DLS) (Fig. 2B & C). In the presence of phosphate, we observe two distinct populations of silk fibroin aggregates - a small nanoscale species with a narrow size distribution that participates in coating formation as visualized by AFM (herein referred to as the “active” species), and a large heterogeneous aggregate species that does not participate in coating formation. The “active” nanoaggregates are predominantly observed at lower phosphate concentrations (100–250 mM), reaching a maximum size around 30–35 nm in diameter at 200 mM phosphate, which coincides with the optimal concentration for coating growth. At these conditions, active nanoaggregates comprise roughly 90% of the species in solution by volume. Furthermore, Figure S3 shows these active nanoaggregates only form at pH 5, where protein surface charges are minimized. Upon further increasing the phosphate concentration (> 200 mM), we observe a reduction in both the diameter and the population of the active nanoaggregate, likely due to salting-out effects previously documented for silk fibroin in the presence of high phosphate concentrations37. The large inactive aggregates are favored at higher phosphate concentrations and are several hundred nanometers to micrometers in diameter. Both the diameter and population of this inactive species increase with phosphate concentration until all observed species in solution are inactive protein precipitates at ≥ 400 mM phosphate, where no coating growth is observed. Here, the kosmotropic phosphate ions likely shift the balance of interactions to favor solution-phase aggregation rather interfacial assembly, inhibiting coating formation at higher phosphate concentrations. It should be noted that our studies are performed at a silk fibroin concentration of 0.05 wt%, which is 10–100 times lower than concentrations typically explored in literature, and thus our self-assembly behavior may vary from previous reports where hydrogelation or nanofiber formation was observed43,44.
To investigate the interplay of protein-protein and protein-surface interactions during coating formation, AFM force spectroscopy measurements between a silk fibroin-coated silica particle and a bare or silk fibroin-coated TiO2 surface were conducted. Not all previous solution pH and phosphate concentration combinations (Fig. 2A) were tested, rather, solution compositions were chosen to compare the interaction strengths at the optimum coating conditions to the most extreme conditions (i.e. the highest and lowest combinations of phosphate and solution pH). Force maps generated between the silk fibroin-coated probe and bare TiO2 are displayed in (Fig. 3A) for the different conditions tested. The force maps show a clear dependence of the interaction strength between silk fibroin and TiO2 on both solution pH and phosphate concentration, with the strongest interactions (white) observed at pH 5 and 200 mM phosphate. Protein-protein interactions measured using the same buffer conditions between the silk fibroin-coated silica particle and an existing silk fibroin coating on TiO2 show similar trends, where the highest interactions are observed at pH 5 and 200 mM phosphate (Fig. 3B). Additional analysis of AFM force measurements plotted as histogram distributions with mean and STDEV of force measurements for each of the tested conditions can be found in (Figure S4). These results support previous insights gained by DLS and suggest that pH near the isoelectric point of silk fibroin facilitates protein-protein and protein-surface interactions by reducing electrostatic repulsion, while phosphate ions play a more complex role in protein assembly in solution and at an interface. Higher phosphate concentrations may promote intramolecular or intra-aggregate interactions (e.g. hydrogen bonding, hydrophobic interactions) at the cost of intermolecular and inter-aggregate interactions. Alternatively, lower phosphate concentrations may not have a strong enough kosmotropic effect to promote intermolecular or intramolecular interactions at all. Continuous coating growth appears to require an intermediate phosphate concentration that balances interactions in a manner that promotes interfacial accumulation without massive solution-phase aggregation.
2.3 β-sheet formation plays a critical role in coating growth. The self-assembly of silk fibroin is characterized by a transition in secondary structure from being predominantly random coil to β-sheet rich, which is often experimentally triggered by kosmotropic salts or dehydrating alcohols22,24,45,46. We investigated the change in secondary structure of silk fibroin over time by ATR-FTIR to elucidate the intermolecular forces at play during coating formation (Fig. 4). Multipeak deconvolution of the amide I region (Figure S5) between 1600 cm-1 − 1700 cm-1 show that early coating timepoints exhibit higher β-sheet content (40.0 ± 0.4% at 15 min) compared to coatings at later timepoints (which stabilize at 33.3 ± 1.2% at 24 hrs). Over the same period, β-turn content increases from 31.3 ± 1.1% to 37.0 ± 2.5% while the α-helix and random coil contributions remain relatively constant. This result is counterintuitive, as we may expect β-sheet content to increase over time as the coatings are continuously exposed to a kosmotropic ion. Secondary structure characterization of the solution-phase silk fibroin nanoaggregates formed in 200 mM phosphate for 24 hours yield β-sheet content of 30.2 ± 7.5%, which is similar to the content observed in coatings after 12 and 24 hours, though the β-turn content is lower than that of the coating (Figure S6, Fitted Results Figure S7). Since the β-sheet content of the coatings decreases asymptotically towards the β-sheet content of solution-phase nanoaggregates, we can surmise that the initial stage of coating formation is dominated by absorption of a relatively β-sheet-rich species while interfacial accumulation of the nanoaggregates observed in solution dominates the later stages.
Infrared photo-induced force microscopy (PiFM) measurements were conducted on the silk fibroin coatings to investigate the spatial distribution of the secondary structures at various timepoints in the coating process (Fig. 5A). PiFM detects infrared adsorption at nano-scale resolution using mechanical forces to measure near-field optical interactions between an AFM tip and sample47. After 15 minutes of coating time, surface-bound silk fibroin nanoglobules appear to be comprised mainly of α-helix/random coil structures (red channel, 1648 cm-1), although bulk measurement of the coating via ATR-FTIR shows a much higher β-sheet content. Surprisingly, the β-sheet signal (green, 1610 cm-1) appears to be homogeneously distributed across the surface rather than isolated within the silk fibroin nanoglobules. These β-sheet structures are much smaller than the silk fibroin nanoglobules and are not detectable by height or phase AFM, likely because they exist as a molecular monolayer. These small β-sheet structures may be formed by monomeric silk fibroin chains (potentially the shorter hydrolyzed fractions of the silk fibroin that results from the degumming process), which diffuse and adsorb to the surface first due to their smaller size. The rapid adsorption of a small β-sheet rich species followed by slower accumulation of relatively β-sheet poor nanoaggregates would be consistent with our ATR-FTIR results (Fig. 4). We hypothesize that the adsorbed β-sheet structures serve as anchoring points, acting as a primer layer for subsequent attachment of the nanoglobules. To test this hypothesis, we removed low molecular weight species from our silk fibroin coating solution using a 100 kDa centrifuge spin filter to compare coating kinetics with and without this primer (Figure S8). Removal of the low molecular weight species appears to reduce coating formation, suggesting that the β-sheet structures play a key role in facilitating coating growth.
Silk fibroin coatings formed in 206 mM NaCl and 500 mM phosphate were also investigated using PiFM to visualize the effect of solution composition on early-stage coating formation and secondary structure distribution (Figure S9). For both conditions, no β-sheet primer layer was detected, which is consistent with the lack of coating formation in these conditions. In the case of NaCl, there is likely insufficient self-assembly to promote the formation of β-sheet rich protein structures, as NaCl is not a kosmotropic salt. Virtually no surface-bound protein was observed for the 500 mM phosphate condition (data not shown). Here the solution-phase self-assembly is likely so strong that all available protein in solution forms large aggregates, as indicated by DLS (Fig. 2C), rather than adsorb to the surface. It is likely that intermediate phosphate concentrations, such as 200 mM, promotes sufficient self-assembly to generate small β-sheet structures at an interface without aggregating all silk fibroin chains into large species that cannot effectively attach to the surface.
To investigate spatial distribution of secondary structures within and between neighboring surface-bound nanoglobules, a 1.5-hour coating was analyzed as the substrate is nearly fully covered by protein at this timepoint (Fig. 5B). Similar to the 15-minute coating, the nanoglobules in the 1.5-hour coating contain a high degree of α-helix/random coil structures. Interestingly, an enrichment of β-sheet and/or β-turn structures (1695 cm-1) is observed at the interfaces between the silk fibroin nanoglobules, suggesting these secondary structures contribute to the cohesive interactions between neighboring globules. The role of hydrogen-bonded β-sheets in providing attractive interactions between silk fibroin chains is supported by a combination of DLS and Thioflavin T assay of our coating solution, which shows that the diameter of solution-phase assemblies increases proportionally with β-sheet content (Figure S10).
2.4 Coating growth occurs in distinct stages, with kinetics that strongly depend on phosphate concentration. Quartz crystal microbalance with dissipation monitoring (QCM-D) was used to investigate silk fibroin coating kinetics in real time. Silk fibroin adsorption kinetics and coating stiffness were determined by simultaneously measuring frequency loss of the oscillating crystal (adsorbed mass) and dissipation of the shear propagating wave through the adsorbed layer (stiffness of the adsorbed layer), respectively. Mass deposition measured via QCM-D within a range of phosphate concentrations at pH 5 are displayed in Fig. 6A. These results are compared with NaCl, which does not promote interfacial self-assembly. In this case, protein accumulation occurs rapidly then saturates, resulting in sub-monolayer coverage, as shown by dry-state AFM (Figure S11A). However, with intermediate phosphate concentrations, silk fibroin self-assembly drives interfacial accumulation past the saturation mass observed with NaCl. Here, mass accumulation initially occurs rapidly then transitions to a slower linear rate at later times. We can divide coating growth three stages: 1) an early-stage, where the substrate is accessible for adsorption of solution-phase species via protein-surface interactions; 2) a late-stage, where the surface is fully covered by protein and the coating process is driven by protein-protein interactions; and 3) a transition region between the early and late stages. Lower phosphate concentrations (< 175 mM) promote faster early-stage and slower late-stage coating growth. Intermediate phosphate concentrations (175–225 mM) promote slower early-stage kinetics but faster late-state growth. Higher phosphate concentrations (> 225 mM) show both slower early-stage and late-state coating growth. In agreement with our ellipsometry measurements, the optimum phosphate concentration for fast, continuous coating growth is 200 mM, where we observe slightly reduced early-stage growth but the fastest rate of late-stage accumulation.
To ensure that the continuous coating growth observed via QCM-D represents protein adsorption rather than coating swelling in the hydrated environment, QCM-D was performed alongside ellipsometry in a specialized flow chamber. This data allows us to decouple dry protein mass (obtained by ellipsometry) from total hydrated mass (obtained by QCM-D). Figure S12 displays the silk fibroin dry mass deposition over time in coating solutions containing different phosphate concentrations. For all phosphate concentrations tested, except for the no-salts condition, dry mass monotonically increases over time, confirming the QCM-D data shows protein mass deposition rather than simply solvent swelling. However, due to the significant variations in flow geometry and dead volume above the substrate, the adsorption kinetics measured in the QCM-D/ellipsometry module cannot be directly compared to that of the standard flow module.
2.5 A structure change marks the transition between early- and late-stage coating growth. To understand the relationship between early and late-stage coating growth, we define a “transition point” where the coating mechanism shifts from being predominantly dependent on protein-surface interactions to protein-protein interactions. To understand this transition, the dissipation-to-frequency (D/F) ratio over time was analyzed from the raw data, which measures the loss of energy through the coating (D) as a function of mass deposited on the surface (F). This ratio gives insight into structural properties of the silk fibroin coating as it forms. A low D/F ratio represents a stiff coating, whereas a high D/F ratio represents a softer coating. To this end, a D/F ratio of lower than 0.05–0.1 is defined as a rigid coating, whereas a value above 0.1 is considered soft. The magnitude of the D/F ratios of our silk fibroin coatings indicate a soft, hydrogel like coating (Fig. 6B), which is supported by liquid-mode AFM imaging (Fig. 1C). In the early stages of coating growth, the coating starts very soft (high D/F ratio), but then rapidly stiffens. The minimum in the D/F ratio can be interpreted as the timepoint in which the coating is at its stiffest, representing a more rigidly bound interfacial layer that is possibly due to protein spreading, a commonly observed phenomenon48–51. After this point, the coating slowly becomes softer as more nanoaggregates accumulate on the existing silk layer rather than the substrate surface during late-stage coating growth. Since the time at which this D/F minimum occurs depends on phosphate concentration, we consider this minimum in the D/F ratio to be the transition point between the early- and late-stage adsorption regions.
Figure 7B & C displays the timepoint and surface mass deposition corresponding to the minimal D/F ratio. As phosphate concentration increases, the time when the minimal D/F is observed also increases, indicating a slower protein-surface adsorption when more phosphate is present (Fig. 7B). However, the total mass absorbed to the surface at this transition point for all phosphate concentrations tested lies in a narrow range from 2750–3100 ng/cm2, despite occurring at different timepoints for each coating condition (Fig. 7C). Comparing these values to the NaCl sample, which undergoes protein-surface adsorption then rapid saturation without further growth, the mass accumulated at the transition point is nearly identical to the mass accumulated at the point of saturation in NaCl (~ 3000 ng/cm2, Fig. 6A, red). This finding suggests that the transition between the early and late-stage coating growth occurs at the apparent “jamming limit” of protein adsorption without the influence of self-assembly. At the “jamming limit”, all accessible surface binding sites are effectively occupied, and without attractive interspecies interactions, protein adsorption is expected to reach equilibrium. However, in our system, silk fibroin can surpass this common limitation of protein adsorption and continuously accumulate at the interface via protein-protein interactions at optimal phosphate concentrations.
2.6 A two-stage modified Langmuir kinetics adsorption model mathematically describes coating growth. Most existing models of protein adsorption, such as Langmuir, Random Sequential Adsorption (RSA), and Brunauer–Emmett–Teller, include a term for the maximum possible surface coverage at equilibrium (“jamming limit”), which changes with protein concentration in solution according to their isotherms 10,48,49,52. Our kinetics studies show that unlike conventional protein adsorption, our self-assembled coatings do not saturate in time and instead experience linear growth at later times. While some existing empirical adsorption models may be designed to accommodate multilayer formation, such as the Freundlich or BET models53, they do not provide parameters that capture our two-stage coating mechanism, nor can they effectively represent the transition from early stage to late-stage coating growth.
Here, we use the Langmuir adsorption equation as a starting point to derive a model mathematically describing our coating process. In conditions that do not promote self-assembly, silk fibroin follows Langmuir-like adsorption, where mass accumulation stops after reaching saturation at sub-monolayer surface coverage (Figure S13A). Although our coating assembly does not meeting all the criteria for true Langmuir adsorption (e.g. in our system, adsorbed species interact laterally and the adsorption process is not fully reversible), we find Langmuir better represents our data compared to RSA (Figure S13B & C) when these models are fit to silk fibroin mass flux to the surface over time in conditions that do not promote self-assembly. While erroneous conclusions can be made by interpreting protein adsorption events as Langmuir when the criteria are not satisfied49, we do not attempt to extract thermodynamic properties such as equilibrium constants or adsorption free energies from this model. Our modified kinetic adsorption model simply serves as a basis for separating early and late-stage coating growth rates.
We start with the original Langmuir kinetic adsorption equation
$$\frac{d\theta }{dt}={k}_{a}{c}^{*}\left(1-\frac{\theta }{{\theta }_{max}}\right)-{k}_{d}\theta$$
1
where \(\theta\) is surface coverage, \({\theta }_{max}\) is the maximum coverage level of the protein, C* is the effective concentration of the bulk protein solution above the surface, and ka and kd are the adsorption and desorption rate constants, respectively. Assuming mass deposition is proportional to surface coverage in the early stages of adsorption, we replace surface coverage with total specific mass deposition in the model, yielding:
$$\frac{dM}{dt}={k}_{a}{c}^{*}\left(1-\frac{M}{{M}_{max}}\right) -{k}_{d}M$$
2
The term dM/dt, the change in specific mass over time, represents protein flux to the surface. To capture our linear late-stage coating growth, a new parameter, kss, is introduced to represent silk fibroin adsorption to surface-bound silk assemblies (protein-protein interactions). To avoid overfitting the data, the desorption term was removed as desorption of silk fibroin after coating formation was experimentally shown to be negligible (Figure S13A).
$$\frac{dM}{dt}={k}_{a}{c}^{*}\left(1-\frac{M}{{M}_{max}}\right)+{k}_{ss}$$
3
The first part of the equation, which resembles the original Langmuir kinetic adsorption equation, represents protein-surface interactions. However, as the surface becomes completely covered with protein globules (Fig. 1), the effect of the surface on adsorbing proteins diminishes. Here, we introduce an exponential function with a fitted time constant (β) to extinguish the protein-surface interaction term as the coating approaches full surface coverage, yielding
$$\frac{dM}{dt}=\frac{{k}_{a}{c}^{*}\left(1-\frac{M}{{M}_{max}}\right)}{{e}^{\beta \times t}} +{k}_{ss}$$
4
As the protein-surface interaction term is extinguished, the model smoothly converges to the kss term, describing late-stage coating growth governed entirely by protein-protein interactions.
Taking the derivative of our coating kinetics curves produced by QCM-D yields protein flux (ng/cm2*s) to the surface over time, allowing us to fit our model to the data and extract early- and late-stage adsorption rate constants, ka (ng/cm2*M*s) and kss (ng/cm2*s), respectively. Since we continuously replenish our concentration boundary layer in a QCM-D flow module, we assume C* to be constant and consider the full term \({k}_{a}{c}^{*}\) to be our early-stage adsorption rate constant. This allows for direct comparison between early- and late-stage adsorption rate constants as they are both in units of (ng/cm2*s). Figure 7A shows the early- and late-stage adsorption rate constants derived from fitting our modified Eq. 4 to the QCM-D data (Fig. 6A) as a function of phosphate concentration (Model Fits and Residuals, see Figure S14). The early-stage adsorption rate constant, \({k}_{a}{c}^{*}\), is maximized at low phosphate concentration, then steadily decreases as phosphate concentration is increased, suggesting that protein-protein assembly negatively impacts the rate of protein-surface adsorption in the initial stages of coating formation. Conversely, the late-stage adsorption rate constant, kss, increases with phosphate concentration until the optimum coating concentration (200 mM) is reached. However, further increasing phosphate concentration decreases kss, suggesting that high phosphate concentrations negatively impact interfacial accumulation regardless of whether it is driven by protein-protein or protein-surface interactions. These fitted results support the DLS results, which show that low phosphate concentrations do not adequately promote self-assembly, thus leading to more conventional adsorption behavior with minimal late-stage growth, while high phosphate concentrations generate mostly large solution-phase that cannot attach to the surface.
We demonstrated the utility of our model by predicting when the transition between early-stage and late-stage coating growth occurs, which we hypothesize to be represented by the inverse time constant (1/β) in our model. We compared 1/β from our fitted data to our experimentally determined transition point, which is where D/F is minimized. Figure 7B shows that the predicted transition time aligns closely with experimentally determined values. Additionally, the total masses accumulated on the surface at these predicted times fall within the range of 2300–3300 ng/cm2, which is in close agreement with experimentally determined values of mass accumulated at the transition point (Fig. 7C). We also utilize our model to quantify variations in coating kinetics when alterations are made to the coating process. As discussed previously, Figure S8 shows a reduction in overall coating formation when the smallest population in the coating solution, the species likely responsible for the β-sheet primer layer, is removed by centrifuge spin filters. Here we find that the early-stage adsorption rate constant decreases from 13.1 to 5.5 ng/cm2s, a 58% decrease, while the late-stage adsorption rate constant decreases 25% from 0.55 to 0.41 ng/cm2s.
2.7 Self-assembled silk fibroin coating provide a facile “universal” method for enhancing surface functionality. Previous work in our lab has demonstrated that our silk fibroin coatings can be used functionalize high aspect ratio poly-L-lactic acid (PLLA) scaffolds for enhanced nerve tissue regeneration 16. To highlight the versatility of our self-assembled coatings, we explored coating properties and stability on a variety of substrates and in a variety of solvents. Figure 8A shows static water contact angle measurements on a variety of surfaces with and without silk coatings. Regardless of the underlying properties of the substrates (e.g., hydrophobic, hydrophilic, organic, inorganic), the coated substrates exhibit a narrow distribution of water contact angles that correspond to the contact angle of bulk silk fibroin, which suggests the ability of our coatings to completely cover a wide range of surfaces. Our coatings are tenaciously bound to surfaces, owing to a multiplexing of adhesive interactions by lateral attractive forces within the coating. Figure S15 shows that these coatings are stable (e.g. no delamination or dissolution) in a variety of harsh solution conditions such as ethanol, toluene, and dichloromethane, and physiologically relevant buffers such as MOPS and phosphate buffer. The high packing density of surface-bound nanoglobules (Fig. 1F) suggests the potential for these silk fibroin coatings to act as non-fouling surfaces. Using QCM-D, Fig. 8B shows that a silk fibroin coating formed on a TiO2 sensor can reduce bovine serum albumin (BSA) adsorption by ~ 85% when compared to a bare TiO2 surface under identical conditions (Fig. 8B, inset) at low BSA concentrations. An anti-fouling study was additionally conducted at a higher BSA concentration (10 mg/mL) on silica resin, where adsorption of FITC-labeled BSA was imaged on both silk fibroin-coated and bare resin (Fig. 8C & D). We observe less overall FITC-BSA adsorption onto the silk-coated resin, demonstrating the nanothin silk fibroin coating can prevent BSA fouling on non-planar surfaces in biologically relevant serum protein concentrations. The ability of silk fibroin to act as a non-fouling surface likely stems from the density of silk fibroin nanoglobules on the surface, leaving no available surface area for unwanted protein adsorption. Additionally, these hydrogel-like silk fibroin coatings are expected to retain bound water, contributing entropic and osmotic repulsive forces to approaching proteins such as BSA54.