Peroxidase-mediated mucin cross-linking drives pathologic mucus gel formation in IL-13-stimulated airway epithelial cells.

Mucus plugs occlude airways to obstruct airflow in asthma. Studies in patients and in mouse models show that mucus plugs occur in the context of type 2 inflammation, and studies in human airway epithelial cells (HAECs) show that interleukin 13 (IL-13) activated cells generate pathologic mucus independently of immune cells. To determine how HAECs autonomously generate pathologic mucus, we used a magnetic microwire rheometer to characterize the viscoelastic properties of mucus secreted under varying conditions. We found that normal HAEC mucus exhibits viscoelastic liquid behavior and that mucus secreted by IL-13 activated HAECs exhibits solid-like behavior caused by mucin cross-linking. In addition, IL-13 activated HAECs show increased peroxidase activity in apical secretions, and an overlaid thiolated polymer (thiomer) solution shows an increase in solid behavior that is prevented by peroxidase inhibition. Furthermore, gene expression for thyroid peroxidase (TPO), but not lactoperoxidase (LPO), is increased in IL-13 activated HAECs and both TPO and LPO catalyze the formation of oxidant acids that cross-link thiomer solutions. Finally, gene expression for TPO in airway epithelial brushings is increased in asthma patients with high airway mucus plug scores. Together, our results show that IL-13 activated HAECs autonomously generate pathologic mucus via peroxidase-mediated cross-linking of mucin polymers.

Air-liquid intestinal cell culture allows in situ rheological characterization of intestinal mucus.

Intestinal health heavily depends on establishing a mucus layer within the gut with physical properties that strike a balance between being sufficiently elastic to keep out harmful pathogens yet viscous enough to flow and turnover the contents being digested. Studies investigating dysfunction of the mucus layer in the intestines are largely confined to animal models, which require invasive procedures to collect the mucus fluid. In this work, we develop a nondestructive method to study intestinal mucus. We use an air-liquid interface culture of primary human intestinal epithelial cells that exposes their apical surface to allow in situ analysis of the mucus layer. Mucus collection is not only invasive but also disrupts the mucus microstructure, which plays a crucial role in the interaction between mucus and the gut microbiome. Therefore, we leverage a noninvasive rheology technique that probes the mechanical properties of the mucus without removal from the culture. Finally, to demonstrate biomedical uses for this cell culture system, we characterize the biochemical and biophysical properties of intestinal mucus due to addition of the cytokine IL-13 to recapitulate the gut environment of Nippostrongylus brasiliensis infection.

Response of lymphatic endothelial cells to combined spatial and temporal variations in fluid flow.

One-way valves within lymphatic vessels are required for the efficient drainage of lymphatic fluids. Fluid flow is proposed to be a key cue in regulating both the formation and maintenance of lymphatic valves. However, to our knowledge, no previous study has systematically examined the response of LECs to the complex combination of spatially and temporally varying fluid flows that occur at lymphatic valves in vivo. We built an in vitro microfluidic device that reproduces key aspects of the flow environment found at lymphatic valves. Using this device, we found that a combination of spatially and temporally varying wall shear stresses (WSSs) led to upregulated transcription of PROX1 and FOXC2. In addition, we observed that combined spatial and temporal variations in WSS-modulated Ca2+ signaling and led to increased cellular levels of NFATc1. These observations suggest that the physical cues generated by the flow environment present within lymphatic valves may act to activate key regulatory pathways that contribute to valve maintenance.

The Role of Membrane-Tethered Mucins in Axial Epithelial Adhesion in Controlled Normal Stress Environments.

The collective adhesive behavior of epithelial cell layers mediated by complex macromolecular fluid environments plays a vital role in many biological processes. Mucins, a family of highly glycosylated proteins, are known to lubricate cell-on-cell contacts in the shear direction. However, the role of mucins mediating axial epithelial adhesion in the direction perpendicular to the plane of the cell sheet has received less attention. This article subjects cell-on-cell layers of live ocular epithelia that express mucins on their apical surfaces to compression/decompression cycles and tensile loading using a customized instrument. In addition to providing compressive moduli of native cell-on-cell layers, it is found that the mucin layer between the epithelia acts as a soft cushion between the epithelial cell layers. Decompression experiments reveal mucin layers act as soft, nonlinear springs in the axial direction. The cell-on-cell layers withstand decompression before fracturing by a cohesive failure within the mucin layer. When mucin deficiency is induced via a protease treatment, it is found that the axial adhesion between the cell layers is increased. The findings which correlate changes in biological factors with changes in mechanical properties might be of interest to challenges in ophthalmology, vision care, and mucus research.

Coalescence of surface bubbles: The crucial role of motion-induced dynamic adsorption layer.

The formation of motion-induced dynamic adsorption layers of surfactants at the surface of rising bubbles is a widely accepted phenomenon. Although their existence and formation kinetics have been theoretically postulated and confirmed in many experimental reports, the investigations primarily remain qualitative in nature. In this paper we present results that, to the best of our knowledge, provide a first quantitative proof of the influence of the dynamic adsorption layer on drainage dynamics of a single foam film formed under dynamic conditions. This is achieved by measuring the drainage dynamics of single foam films, formed by air bubbles of millimetric size colliding against the interface between n-octanol solutions and air. This was repeated for a total of five different surfactant concentrations and two different liquid column heights. All three steps preceding foam film rupture, namely the rising, bouncing and drainage steps, were sequentially examined. In particular, the morphology of the single film formed during the drainage step was analyzed considering the rising and bouncing history of the bubble. It was found that, depending on the motion-induced state of adsorption layer at the bubble surface during the rising and the bouncing steps, single foam film drainage dynamics can be spectacularly different. Using Direct Numerical Simulations (DNS), it was revealed that surfactant redistribution can occur at the bubble surface as a result of the bouncing dynamics (approach-bounce cycles), strongly affecting the interfacial mobility, and leading to slower rates of foam film drainage. Since the bouncing amplitude directly depends on the rising velocity, which correlates in turn with the adsorption layer of surfactants at the bubble surface during the rising step, it is demonstrated that the lifetime of surface bubbles should intimately be related to the history of their formation.

Stable High-Concentration Monoclonal Antibody Formulations Enabled by an Amphiphilic Copolymer Excipient.

Monoclonal antibodies are a staple in modern pharmacotherapy. Unfortunately, these biopharmaceuticals are limited by their tendency to aggregate in formulation, resulting in poor stability and often requiring low concentration drug formulations. Moreover, existing excipients designed to stabilize these formulations are often limited by their toxicity and tendency to form particles such as micelles. Here, we demonstrate the ability of a simple "drop-in", amphiphilic copolymer excipient to enhance the stability of high concentration formulations of clinically-relevant monoclonal antibodies without altering their pharmacokinetics or injectability. Through interfacial rheology and surface tension measurements, we demonstrate that the copolymer excipient competitively adsorbs to formulation interfaces. Further, through determination of monomeric composition and retained bioactivity through stressed aging, we show that this excipient confers a significant stability benefit to high concentration antibody formulations. Finally, we demonstrate that the excipient behaves as an inactive ingredient, having no significant impact on the pharmacokinetic profile of a clinically relevant antibody in mice. This amphiphilic copolymer excipient demonstrates promise as a simple formulation additive to create stable, high concentration antibody formulations, thereby enabling improved treatment options such as a route-of-administration switch from low concentration intravenous (IV) to high concentration subcutaneous (SC) delivery while reducing dependence on the cold chain.

Effect of Recombinant Human Lubricin on Model Tear Film Stability.

Purpose: To investigate and quantify the effect of recombinant human lubricin (rh-lubricin) on model tear film stability. Methods: A custom-built, interferometry-based instrument called the Interfacial Dewetting and Drainage Optical Platform was used to create and record the spatiotemporal evolution of model acellular tear films. Image segmentation and analysis was performed in MATLAB to extract the most essential features from the wet area fraction versus time curve, namely the evaporative break-up time and the final wet area fraction (A10). These two parameters indicate the tear film stability in the presence of rh-lubricin in its unstressed and stressed forms. Results: Our parameters successfully captured the trend of increasing tear film stability with increasing rh-lubricin concentration, and captured differences in rh-lubricin efficacy after various industrially relevant stresses. Specifically, aggregation and fragmentation caused by a 4-week, high temperature stress condition negatively impacted rh-lubricin's ability to maintain model tear film stability. Adsorbed rh-lubricin alone was not sufficient to resist break-up and maintain full area coverage of the model tear film surface. Conclusions: Our results demonstrate that fragmentation and aggregation can negatively impact rh-lubricin's ability to maintain a stable tear film. In addition, the ability of rh-lubricin to maintain wetted area coverage is due to both freely dispersed and adsorbed rh-lubricin. Translational Relevance: Our platform and analysis method provide a facile, intuitive, and clinically relevant means to quantify the effect of ophthalmic drugs and formulations intended for improving tear film stability, as well as capture differences between variants related to drug stability and efficacy.

Tear Film Stability as a Function of Tunable Mucin Concentration Attached to Supported Lipid Bilayers.

In this work, we describe the development of a tunable, acellular in vitro model of the mucin layer of the human tear film. First, supported lipid bilayers (SLBs) comprised of the phospholipid DOPC (1,2-dioleoyl-sn-glycero-3-phosphocholine) and biotinyl cap PE (1,2-dioleoyl-sn-glycero-3-phosphoethanolamine-N-(cap biotinyl)) are created on the surface of a glass dome with radius of curvature comparable to the human eye. Next, biotinylated bovine submaxillary mucins (BSM) are tethered onto the SLB using streptavidin protein. The mucin presentation can be tuned by altering the concentration of biotinylated BSM, which we confirm using fluorescence microscopy. Due to the optically smooth surface that results, this model is compatible with interferometry for monitoring film thickness. Below a certain level of mucin coverage, we observe short model tear film breakup times, mimicking a deficiency in membrane-associated mucins. In contrast, the breakup time is significantly delayed for SLBs with high mucin coverage. Because no differences in mobility or wettability were observed, we hypothesize that higher mucin coverage provides a thicker hydrated layer that can protect against external disturbances to thin film stability. This advance paves the way for a more physiological, interferometry-based in vitro model for investigating tear film breakup.

Systematic characterization of effect of flow rates and buffer compositions on double emulsion droplet volumes and stability.

Double emulsion droplets (DEs) are water/oil/water droplets that can be sorted via fluorescence-activated cell sorting (FACS), allowing for new opportunities in high-throughput cellular analysis, enzymatic screening, and synthetic biology. These applications require stable, uniform droplets with predictable microreactor volumes. However, predicting DE droplet size, shell thickness, and stability as a function of flow rate has remained challenging for monodisperse single core droplets and those containing biologically-relevant buffers, which influence bulk and interfacial properties. As a result, developing novel DE-based bioassays has typically required extensive initial optimization of flow rates to find conditions that produce stable droplets of the desired size and shell thickness. To address this challenge, we conducted systematic size parameterization quantifying how differences in flow rates and buffer properties (viscosity and interfacial tension at water/oil interfaces) alter droplet size and stability, across 6 inner aqueous buffers used across applications such as cellular lysis, microbial growth, and drug delivery, quantifying the size and shell thickness of >22 000 droplets overall. We restricted our study to stable single core droplets generated in a 2-step dripping-dripping formation regime in a straightforward PDMS device. Using data from 138 unique conditions (flow rates and buffer composition), we also demonstrated that a recent physically-derived size law of Wang et al. can accurately predict double emulsion shell thickness for >95% of observations. Finally, we validated the utility of this size law by using it to accurately predict droplet sizes for a novel bioassay that requires encapsulating growth media for bacteria in droplets. This work has the potential to enable new screening-based biological applications by simplifying novel DE bioassay development.

A Mucin-Deficient Ocular Surface Mimetic Platform for Interrogating Drug Effects on Biolubrication, Antiadhesion Properties, and Barrier Functionality.

Dry eye disease (DED) affects more than 100 million people worldwide, causing significant patient discomfort and imposing a multi-billion-dollar burden on global health care systems. In DED patients, the natural biolubrication process that facilitates pain-free blinking goes awry due to an imbalance of lipids, aqueous medium, and mucins in the tear film, resulting in ocular surface damage. Identifying strategies to reduce adhesion and shear stresses between the ocular surface and the conjunctival cells lining the inside of the eyelid during blink cycles is a promising approach to improve the signs and symptoms of DED. However, current preclinical models for screening ocular lubricants rely on scarce, heterogeneous tissue samples or model substrates that do not capture the complex biochemical and biophysical cues present at the ocular surface. To recapitulate the hierarchical architecture and phenotype of the ocular interface for preclinical drug screening, we developed an in vitro mucin-deficient DED model platform that mimics the complexity of the ocular interface and investigated its utility in biolubrication, antiadhesion, and barrier protection studies using recombinant human lubricin, a promising investigational therapy for DED. The biomimetic platform recapitulated the pathological changes in biolubrication, adhesion, and barrier functionality often observed in mucin-deficient DED patients and demonstrated that recombinant human lubricin can reverse the damage induced by mucin loss in a dose- and conformation-dependent manner. Taken together, these results highlight the potential of the platform─and recombinant human lubricin─in advancing the standard of care for mucin-deficient DED patients.

Engineering Insulin Cold Chain Resilience to Improve Global Access.

There are 150 million people with diabetes worldwide who require insulin replacement therapy, and the prevalence of diabetes is rising the fastest in middle- and low-income countries. The current formulations require costly refrigerated transport and storage to prevent loss of insulin integrity. This study shows the development of simple "drop-in" amphiphilic copolymer excipients to maintain formulation integrity, bioactivity, pharmacokinetics, and pharmacodynamics for over 6 months when subjected to severe stressed aging conditions that cause current commercial formulation to fail in under 2 weeks. Further, when these copolymers are added to Humulin R (Eli Lilly) in original commercial packaging, they prevent insulin aggregation for up to 4 days at 50 °C compared to less than 1 day for Humulin R alone. These copolymers demonstrate promise as simple formulation additives to increase the cold chain resilience of commercial insulin formulations, thereby expanding global access to these critical drugs for treatment of diabetes.

Flowering in bursting bubbles with viscoelastic interfaces.

The lifetime of bubbles, from formation to rupture, attracts attention because bubbles are often present in natural and industrial processes, and their geometry, drainage, coarsening, and rupture strongly affect those operations. Bubble rupture happens rapidly, and it may generate a cascade of small droplets or bubbles. Once a hole is nucleated within a bubble, it opens up with a variety of shapes and velocities depending on the liquid properties. A range of bubble rupture modes are reported in literature in which the reduction of a surface energy drives the rupture against inertial and viscous forces. The role of surface viscoelasticity of the liquid film in this colorful scenario is, however, still unknown. We found that the presence of interfacial viscoelasticity has a profound effect in the bubble bursting dynamics. Indeed, we observed different bubble bursting mechanisms upon the transition from viscous-controlled to surface viscoelasticity-controlled rupture. When this transition occurs, a bursting bubble resembling the blooming of a flower is observed. A simple modeling argument is proposed, leading to the prediction of the characteristic length scales and the number and shape of the bubble flower petals, thus paving the way for the control of liquid formulations with surface viscoelasticity as a key ingredient. These findings can have important implications in the study of bubble dynamics, with consequences for the numerous processes involving bubble rupture. Bubble flowering can indeed impact phenomena such as the spreading of nutrients in nature or the life of cells in bioreactors.

Mucin-Like Glycoproteins Modulate Interfacial Properties of a Mimetic Ocular Epithelial Surface.

Dry eye disease (DED) has high personal and societal costs, but its pathology remains elusive due to intertwined biophysical and biochemical processes at the ocular surface. Specifically, mucin deficiency is reported in a subset of DED patients, but its effects on ocular interfacial properties remain unclear. Herein a novel in vitro mucin-deficient mimetic ocular surface (Mu-DeMOS) with a controllable amount of membrane-tethered mucin molecules is developed to represent the diseased ocular surfaces. Contact angle goniometry on mimetic ocular surfaces reveals that high surface roughness, but not the presence of hydrophilic mucin molecules, delivers constant hydration over native ocular surface epithelia. Live-cell rheometry confirms that the presence of mucin-like glycoproteins on ocular epithelial cells reduces shear adhesive strength at cellular interfaces. Together, optimal surface roughness and surface chemistry facilitate sustainable lubrication for healthy ocular surfaces, while an imbalance between them contributes to lubrication-related dysfunction at diseased ocular epithelial surfaces. Furthermore, the restoration of low adhesive strength at Mu-DeMOS interfaces through a mucin-like glycoprotein, recombinant human lubricin, suggests that increased frictional damage at mucin-deficient cellular surfaces may be reversible. More broadly, these results demonstrate that Mu-DeMOS is a promising platform for drug screening assays and fundamental studies on ocular physiology.

Adsorption and Aggregation of Monoclonal Antibodies at Silicone Oil-Water Interfaces.

Monoclonal antibody (mAb) therapies are rapidly growing for the treatment of various diseases like cancer and autoimmune disorders. Many mAb drug products are sold as prefilled syringes and vials with liquid formulations. Typically, the walls of prefilled syringes are coated with silicone oil to lubricate the surfaces during use. MAbs are surface-active and adsorb to these silicone oil-solution interfaces, which is a potential source of aggregation. We studied formulations containing two different antibodies, mAb1 and mAb2, where mAb1 aggregated more when agitated in the presence of an oil-water interface. This directly correlated with differences in surface activity of the mAbs, studied with interfacial tension, surface mass adsorption, and interfacial rheology. The difference in interfacial properties between the mAbs was further reinforced in the coalescence behavior of oil droplets laden with mAbs. We also looked at the efficacy of surfactants, typically added to stabilize mAb formulations, in lowering adsorption and aggregation of mAbs at oil-water interfaces. We showed the differences between poloxamer-188 and polysorbate-20 in competing with mAbs for adsorption to interfaces and in lowering particulate and overall aggregation. Our results establish a direct correspondence between the adsorption of mAbs at oil-water interfaces and aggregation and the effect of surfactants in lowering aggregation by competitively adsorbing to these interfaces.

In-Use Interfacial Stability of Monoclonal Antibody Formulations Diluted in Saline i.v. Bags.

The use of monoclonal antibodies (mAbs) for the treatment of a variety of diseases is rapidly growing each year. Many mAbs are administered intravenously using i.v. bags containing 0.9% NaCl (normal saline). We studied the aggregation propensity of these antibody solutions in saline and compared it with a low ionic strength formulation buffer. The mAb studied in this work is prone to aggregate, and is known to form a viscoelastic network at the air-solution interface. We observed that this interfacial elasticity increased when formulated in saline. In the bulk, the mAbs exhibited a tendency to self-associate that was higher in saline. We also studied the aggregation of the mAbs in the presence of polysorbate-20, typically added to formulations to mitigate interfacial aggregation. We observed that with surfactants, the presence of salt in the buffer led to a greater mAb adsorption at the interface and resulted in the formation of more particulate aggregates. Our results show that the addition of salt to the buffer led to differences in the interfacial aggregation in mAb formulations, showing that stress studies used to screen for mAb aggregation intended for i.v. administration should be performed in conditions representative of their intended route of administration.

Tuning corneal epithelial cell adhesive strength with varying crosslinker content in silicone hydrogel materials.

Purpose: To quantify the effect of silicone hydrogel crosslink density on the adhesion at corneal epithelial cells/silicone hydrogel contact lens interface. Methods: A custom-built rheometer, referred to as the live cell monolayer rheometer, was used to measure the adhesive strengths between corneal epithelial cell monolayers and silicone hydrogel lens surfaces. The resulting stress relaxations of senofilcon A-derived silicone hydrogel materials with varying crosslinking densities and delefilcon A were tested. Senofilcon A-like materials labeled L1, L2, L3, L4, and L5 contained crosslinker concentrations of 1.2, 1.35, 1.5, 1.65, and 1.8 wt%, respectively. The residual modulus measured from the live cell monolayer rheometer provided a direct indication of adhesive attachment. Results: Within the senofilcon-derived series, the adhesive strength shows a surprising minimum with respect to crosslink density. Specifically, L1 (1.20%) has the highest adhesive strength of 39.5 ± 11.2 Pa. The adhesive strength diminishes to a minimum of 11.2 ± 2.1 Pa for L3, whereafter it increases to 14.5 ± 2.5 Pa and 18.1 ± 5.1 Pa for L4 and L5, respectively. The delefilcon A lens exhibits a comparable adhesive strength of 27.8 ± 6.3 Pa to L1. Conclusions: These results demonstrated that increasing the crosslink density has a nonmonotonic influence on the adherence of lenses to mucin-expressing corneal epithelial cells, which suggests a competition mechanism at the cell/lens interface. Translational Relevance: Because the adhesiveness of contact lenses to ocular tissues may impact the comfort level for lens wearers and affect ease of removal, this study suggests that lens adhesion can be optimized through the control of crosslink density.

Understanding the adsorption and potential tear film stability properties of recombinant human lubricin and bovine submaxillary mucins in an in vitro tear film model.

The wetting and adsorption properties for two glycoproteins, recombinant human lubricin and bovine submaxillary mucins (BSM) were evaluated on hydrophilic and hydrophobic glass dome surfaces in a simplified in vitro tear film model. We show that both recombinant human lubricin (rh-lubricin) and BSM solutions render surfaces hydrophilic and when the fluid films reach 500 nm or less, the fluids resist evaporation-driven breakup through a volumetric flux across the surface, which we believe is due to evaporation-driven solutocapillary flows. rh-Lubricin was able to maintain a wet film without spontaneous breakup for longer periods of time than BSM at lower concentrations, which we attribute to differences in adsorption properties, measured by QCM-D, that result from surface charge and structural differences (confirmed by zeta potential, DLS, and SAXS measurements).

Evaporation-induced Rayleigh-Taylor instabilities in polymer solutions.

Understanding the mechanics of detrimental convective instabilities in drying polymer solutions is crucial in many applications such as the production of film coatings. It is well known that solvent evaporation in polymer solutions can lead to Rayleigh-Bénard or Marangoni-type instabilities. Here, we reveal another mechanism, namely that evaporation can cause the interface to display Rayleigh-Taylor instabilities due to the build-up of a dense layer at the air-liquid interface. We study experimentally the onset time (tp) of the instability as a function of the macroscopic properties of aqueous polymer solutions, which we tune by varying the polymer concentration (c0), molecular weight and polymer type. In dilute solutions, tp shows two limiting behaviours depending on the polymer diffusivity. For high diffusivity polymers (low molecular weight), the pluming time scales as [Formula: see text]. This result agrees with previous studies on gravitational instabilities in miscible systems where diffusion stabilizes the system. On the other hand, in low diffusivity polymers the pluming time scales as [Formula: see text]. The stabilizing effect of an effective interfacial tension, similar to those in immiscible systems, explains this strong concentration dependence. Above a critical concentration, [Formula: see text], viscosity delays the growth of the instability, allowing time for diffusion to act as the dominant stabilizing mechanism. This results in tp scaling as (ν/c0)2/3. This article is part of the theme issue 'Stokes at 200 (Part 1)'.

Polymeric-nanofluids stabilized emulsions: Interfacial versus bulk rheology.

HYPOTHESIS: The properties of oil-in-water emulsions are influenced by the rheology of the aqueous phase (continuous phase) and the rheology of the oil-water interfaces. The bulk and interfacial rheological parameters can be tuned by incorporating nanoparticles (NPs) featuring different surface chemistries and polymers with different chemical or physical structures. Therefore, NPs and polymers can be used to formulate emulsions with different properties. EXPERIMENTS: The viscoelasticity at the oil-(aqueous phase) interface and the bulk viscoelasticity of aqueous phase were investigated in the presence of different fumed silica NPs (i.e., hydrophilic, hydrophobic, and slightly hydrophobic) and polymers with two different molecular weights. Bulk and interfacial viscoelastic properties were investigated, employing oscillatory rheological techniques. Furthermore, morphology and stability of the oil-in-(aqueous nanofluid) emulsions were explored utilizing bulk emulsification and single drop coalescence experiments. FINDINGS: Introducing polymers into the aqueous nanofluids had opposite effects on bulk and interfacial viscoelasticity. Despite the significant increase in bulk viscoelasticity upon addition of polymers into the aqueous nanofluids, the interfacial viscoelasticity and emulsion stability considerably decreased. The slightly hydrophobic NP nanofluids without polymers showed no bulk viscoelasticity, but displayed the highest interfacial viscoelasticity and emulsion stability. This provided us a unique opportunity to unravel the importance of bulk and interfacial viscoelasticity on oil-in-water emulsification and proved the dominant role of interfacial viscoelasticity over bulk viscoelasticity on emulsion stability.

Correction: Viscoelastic interfaces comprising of cellulose nanocrystals and lauroyl ethyl arginate for enhanced foam stability.

Correction for 'Viscoelastic interfaces comprising of cellulose nanocrystals and lauroyl ethyl arginate for enhanced foam stability' by Agnieszka Czakaj et al., Soft Matter, 2020, 16, 3981-3990, DOI: .