Abstract: Computational fluid dynamics (CFD) is a technique that is used increasingly in the biomedical field. Solving the flow equations numerically provides a convenient way to assess the efficiency of therapies and devices, ranging from cardiovascular stents and heart valves to hemodialysis workflows. Also in the respiratory field CFD has gained increasing interest, especially through the combination of three dimensional image reconstruction which results in highend patient-specific models. This paper provides an overview of clinical applications of CFD through image based modeling, resulting from recent studies performed in our center. We focused on two applications: assessment of the efficiency of inhalation medication and analysis of endobronchial valve placement. In the first application we assessed the mode of action of a novel bronchodilator in 10 treated patients and 4 controls. We assessed the local volume increase and resistance change based on the combination of imaging and CFD. We found a good correlation between the changes in volume and resistance coming from the CFD results and the clinical tests. In the second application we assessed the placement and effect of one way endobronchial valves on respiratory function in 6 patients. We found a strong patientspecific result of the therapy where in some patients the therapy resulted in complete atelectasis of the target lobe while in others the lobe remained inflated. We concluded from these applications that CFD can provide a better insight into clinically relevant therapies.

Abstract: Introduction: Noninvasive ventilation (NIV) is a well-established treatment for acute-on-chronic respiratory failure in hypercapnic COPD patients. Less is known about the effects of a long-term treatment with NIV in hypercapnic COPD patients and about the factors that may predict response in terms of improved oxygenation and lowered CO2 retention.Methods: In this study, we randomized 15 patients to a routine pharmacological treatment (n = 5, age 66 [standard deviation ± 6] years, FEV1 30.5 [±5.1] %pred, PaO2 65 [±6] mmHg, PaCO2 52.4 [±6.0] mmHg) or to a routine treatment and NIV (using the Synchrony BiPAP device [Respironics, Inc, Murrsville, PA]) (n = 10, age 65 [±7] years, FEV1 29.5 [±9.0] %pred, PaO2 59 [±13] mmHg, PaCO2 55.4 [±7.7] mmHg) for 6 months. We looked at arterial blood gasses, lung function parameters and performed a low-dose computed tomography of the thorax, which was later used for segmentation (providing lobe and airway volumes, iVlobe and iVaw) and post-processing with computer methods (providing airway resistance, iRaw) giving overall a functional image of the separate airways and lobes.Results: In both groups there was a nonsignificant change in FEV1 (NIV group 29.5 [9.0] to 38.5 [14.6] %pred, control group 30.5 [5.1] to 36.8 [8.7] mmHg). PaCO2 dropped significantly only in the NIV group (NIV: 55.4 [7.7] → 44.5 [4.70], P = 0.0076; control: 52.4 [6.0] → 47.6 [8.2], NS). Patients actively treated with NIV developed a more inhomogeneous redistribution of mass flow than control patients. Subsequent analysis indicated that in NIV-treated patients that improve their blood gases, mass flow was also redistributed towards areas with higher vessel density and less emphysema, indicating that flow was redistributed towards areas with better perfusion. There was a highly significant correlation between the % increase in mass flow towards lobes with a blood vessel density of >9% and the increase in PaO2. Improved ventilation–perfusion match and recruitment of previously occluded small airways can explain the improvement in blood gases.Conclusion: We can conclude that in hypercapnic COPD patients treated with long-term NIV over 6 months, a mass flow redistribution occurs, providing a better ventilation–perfusion match and hence better blood gases and lung function. Control patients improve homogeneously in iVaw and iRaw, without improvement in gas exchange since there is no improved ventilation/perfusion ratio or increased alveolar ventilation. These differences in response can be detected through functional imaging, which gives a more detailed report on regional lung volumes and resistances than classical lung function tests do. Possibly only patients with localized small airway disease are good candidates for long-term NIV treatment. To confirm this and to see if better arterial blood gases also lead to better health related quality of life and longer survival, we have to study a larger population.

Abstract: Background: Salbutamol and ipratropium bromide improve lung function in patients with chronic obstructive pulmonary disease (COPD). However, their bronchodilating effect has not yet been compared in the central and distal airways. Functional imaging using computational fluid dynamics offers the possibility of making such a comparison. The objective of this study was to assess the effects of salbutamol and ipratropium bromide on the geometry and computational fluid dynamics-based resistance of the central and distal airways.Methods: Five patients with Global Initiative for Chronic Obstructive Lung Disease Stage III COPD were randomized to a single dose of salbutamol or ipratropium bromide in a crossover manner with a 1-week interval between treatments. Patients underwent lung function testing and a multislice computed tomography scan of the thorax that was used for functional imaging. Two hours after dosing, the patients again underwent lung function tests and repeat computed tomography.Results: Lung function parameters, including forced expiratory volume in 1 second, vital capacity, overall airway resistance, and specific airway resistance, changed significantly after administration of each product. On functional imaging, the bronchodilating effect was greater in the distal airways, with a corresponding drop in airway resistance, compared with the central airways. Salbutamol and ipratropium bromide were equally effective at first glance when looking at lung function tests, but when viewed in more detail with functional imaging, hyporesponsiveness could be shown for salbutamol in one patient. Salbutamol was more effective in the other patients.Conclusion: This pilot study gives an innovative insight into the modes of action of salbutamol and ipratropium bromide in patients with COPD, using the new techniques of functional imaging and computational fluid dynamics.

Abstract: Quantitative atomic resolution annular dark field scanning transmission electron microscopy (ADF STEM) has become a powerful technique for nanoparticle atom-counting. However, a lot of nanoparticles provide a severe characterisation challenge because of their limited size and beam sensitivity. Therefore, quantitative ADF STEM may greatly benefit from statistical detection theory in order to optimise the instrumental microscope settings such that the incoming electron dose can be kept as low as possible whilst still retaining single-atom precision. The principles of detection theory are used to quantify the probability of error for atom-counting. This enables us to decide between different image performance measures and to optimise the experimental detector settings for atom-counting in ADF STEM in an objective manner. To demonstrate this, ADF STEM imaging of an industrial catalyst has been conducted using the near-optimal detector settings. For this experiment, we discussed the limits for atomcounting diagnosed by combining a thorough statistical method and detailed image simulations.

Abstract: Since its discovery 20 years ago, many studies have been performed to gain insight into the functional role of cytoglobin (Cygb). However, Cygb has been proven to be a promiscuous protein. Yet, there is a consensus that Cygb is a cytoprotective protein involved in redox homeostasis. CYGB is a ubiquitously expressed hexacoordinated globin that is highly expressed in melanocytes and is often found to be downregulated during melanocyte-to-melanoma transition. In Chapter III, we investigated the molecular mechanism through which CYGB could be involved in redox regulation. Here, we showed that CYGB contains two redox-sensitive cysteine residues and that the formation of an intramolecular disulfide bridge resulted in the heme group becoming more accessible to external ligands. This supports the hypothesis that Cys38 and Cys83 serve as sensitive redox sensors. In Chapter IV we showed that CYGB mRNA and protein levels were elevated upon exposure to hypoxia. Interestingly, this upregulation was most likely HIF-2α-dependent. We propose that in melanoma, HIF-2α, rather than HIF-1α, positively regulates CYGB under hypoxic conditions in a cell type specific way. In Chapter V, the cytotoxic effect of indirect NTP treatment in two melanoma cell lines with divergent endogenous CYGB expression levels was investigated. We confirmed that NTP endows cytotoxicity that induces cell death through apoptosis and that this was mediated through the production of ROS. Moreover, we showed that CYGB protects melanoma cells from ROS-induced apoptosis by the scavenging of ROS. Interestingly, CYGB expression influenced the expression of NRF2 and HO-1. We identified the lncRNA MEG3 as a possible mechanism through which NRF2 expression and its downstream target HO-1 can be regulated by CYGB. In chapter VI, increased basal ROS levels and higher degree of lipid peroxidation upon RSL3 treatment contributed to the increased sensitivity of CYGB knockdown G361 cells to ferroptosis. Furthermore, transcriptome analysis demonstrates the enrichment of multiple cancer malignancy pathways upon CYGB knockdown, supporting a tumor-suppressive role for CYGB. Remarkably, CYGB expression regulation was identified as a critical determinant of the ferroptosis–pyroptosis therapy response. This suggests that CYGB is involved in the regulation of multiple modes of programmed cell death. FInally, we sought to delineate the RONS that are responsible for plasma-induced ICD. Our results highlight the importance of the short-lived species. Furthermore, we are first to demonstrate that NTP-created vaccine is safely prepared and offers complete protection. Moreover, we provide conclusive evidence that direct application of NTP induces ICD in melanoma.

Abstract: In the present paper, a statistical model-based method to count the number of atoms of monotype crystalline nanostructures from high resolution high-angle annular dark-field (HAADF) scanning transmission electron microscopy (STEM) images is discussed in detail together with a thorough study on the possibilities and inherent limitations. In order to count the number of atoms, it is assumed that the total scattered intensity scales with the number of atoms per atom column. These intensities are quantitatively determined using model-based statistical parameter estimation theory. The distribution describing the probability that intensity values are generated by atomic columns containing a specific number of atoms is inferred on the basis of the experimental scattered intensities. Finally, the number of atoms per atom column is quantified using this estimated probability distribution. The number of atom columns available in the observed STEM image, the number of components in the estimated probability distribution, the width of the components of the probability distribution, and the typical shape of a criterion to assess the number of components in the probability distribution directly affect the accuracy and precision with which the number of atoms in a particular atom column can be estimated. It is shown that single atom sensitivity is feasible taking the latter aspects into consideration.

Abstract: Statistical parameter estimation theory is proposed as a quantitative method to measure unknown structure parameters from electron microscopy images. Images are then purely considered as data planes from which structure parameters have to be determined as accurately and precisely as possible using a parametric statistical model of the observations. For this purpose, an efficient algorithm is proposed for the estimation of atomic column positions and intensities from high angle annular dark field (HAADF) scanning transmission electron microscopy (STEM) images. Furthermore, the so-called CramérRao lower bound (CRLB) is reviewed to determine the limits to the precision with which continuous parameters such as atomic column positions and intensities can be estimated. Since this lower bound can only be derived for continuous parameters, alternative measures using the principles of detection theory are introduced for problems concerning the estimation of discrete parameters such as atomic numbers. An experimental case study is presented to show the practical use of these measures for the optimization of the experiment design if the purpose is to decide between the presence of specific atom types using STEM images.

Abstract: An efficient model-based estimation algorithm is introduced in order to quantify the atomic column positions and intensities from atomic resolution (scanning) transmission electron microscopy ((S)TEM) images. This algorithm uses the least squares estimator on image segments containing individual columns fully accounting for the overlap between neighbouring columns, enabling the analysis of a large field of view. For this algorithm, the accuracy and precision with which measurements for the atomic column positions and scattering cross-sections from annular dark field (ADF) STEM images can be estimated, is investigated. The highest attainable precision is reached even for low dose images. Furthermore, advantages of the model- based approach taking into account overlap between neighbouring columns are highlighted. To provide end-users this well-established quantification method, a user friendly program, StatSTEM, is developed which is freely available under a GNU public license.

Abstract: Many current anti-cancer therapies rely on increasing the intracellular reactive oxygen and nitrogen species (RONS) contents with the aim to induce irreparable damage, which subsequently results in tumor cell death. A novel tool in cancer therapy is the use of cold atmospheric plasma (CAP), which has been found to be very effective in the treatment of many different cancer cell types in vitro as well as in vivo, mainly through the vast generation of RONS. One of the key determinants of the cell's fate will be the interaction of RONS, generated by CAP, with important proteins, i.e. redox-regulatory proteins. One such protein is cytoglobin (CYGB), a recently discovered globin proposed to be involved in the protection of the cell against oxidative stress. In this study, the effect of plasma-produced RONS on CYGB was investigated through the treatment of CYGB with CAP for different treatment times. Spectroscopic analysis of CYGB showed that although chemical modifications occur, its secondary structure remains intact. Mass spectrometry experiments identified these modifications as oxidations of mainly sulfur-containing and aromatic amino acids. With longer treatment time, the treatment was also found to induce nitration of the heme. Furthermore, the two surface-exposed cysteine residues of CYGB were oxidized upon treatment, leading to the formation of intermolecular disulfide bridges, and potentially also intramolecular disulfide bridges. In addition, molecular dynamics and docking simulations confirmed, and further show, that the formation of an intramolecular disulfide bond, due to oxidative conditions, affects the CYGB 3D structure, thereby opening the access to the heme group, through gate functioning of His117. Altogether, the results obtained in this study (1) show that plasma-produced RONS can extensively oxidize proteins and (2) that the oxidation status of two redox-active cysteines lead to different conformations of CYGB.

Abstract: We present a quantitative method which allows us to reliably measure dynamic changes in the atomic structure of monatomic crystalline nanomaterials from a time series of atomic resolution annular dark field scanning transmission electron microscopy images. The approach is based on the so-called hidden Markov model and estimates the number of atoms in each atomic column of the nanomaterial in each frame of the time series. We discuss the origin of the improved performance for time series atom-counting as compared to the current state-of-the-art atom-counting procedures, and show that the so-called transition probabilities that describe the probability for an atomic column to lose or gain one or more atoms from frame to frame are particularly important. Using these transition probabilities, we show that the method can also be used to estimate the probability and cross section related to structural changes. Furthermore, we explore the possibilities for applying the method to time series recorded under variable environmental conditions. The method is shown to be promising for a reliable quantitative analysis of dynamic processes such as surface diffusion, adatom dynamics, beam effects, or in situ experiments.

Abstract: Bacterial and mammalian proteins, such as lysozyme, are gaining increasing interest as anticancer drugs. This study aims to modify the lysozyme structure using cold atmospheric plasma to boost its cancer cell killing effect. We investigated the structure at acidic and neutral pH using various experimental techniques (circular dichroism, fluorescence, and mass spectrometry) and molecular dynamics simulations. The controlled structural modification of lysozyme at neutral pH enhances its activity, while the activity was lost at acidic pH at the same treatment conditions. Indeed, a larger number of amino acids were oxidized at acidic pH after plasma treatment, which results in a greater distortion of the lysozyme structure, whereas only limited structural changes were observed in lysozyme after plasma treatment at neutral pH. We found that the plasma-treated lysozyme significantly induced apoptosis to the cancer cells. Our results reveal that plasma-treated lysozyme could have potential as a new cancer cell killing drug.

Abstract: By combining statistical parameter estimation and model-order selection using a Bayesian framework, the maximum a posteriori (MAP) probability rule is proposed in this chapter as an objective and quantitative method to detect atom columns from high-resolution scanning transmission electron microscopy (HRSTEM) images. The validity and usefulness of this approach is demonstrated to both simulated and experimental annular dark-field (ADF) STEM images, but also to simultaneously acquired annular bright-field (ABF) and ADF STEM image data.

Abstract: An efficient model-based estimation algorithm is introduced to quantify the atomic column positions and intensities from atomic-resolution (scanning) transmission electron microscopy ((S)TEM) images. This algorithm uses the least squares estimator on image segments containing individual columns fully accounting for overlap between neighboring columns, enabling the analysis of a large field of view. To provide end-users with this well-established quantification method, a user friendly program, StatSTEM, is developed which is freely available under a GNU public license. In this chapter, this efficient algorithm is applied to three different nanostructures for which the analysis of a large field of view is required.

Abstract: The maximum a posteriori (MAP) probability rule for atom column detection can also be used as a tool to evaluate the relation between scanning transmission electron microscopy (STEM) image quality and atom detectability. In this chapter, a new image-quality measure is proposed that correlates well with atom detectability, namely the integrated contrast-to-noise ratio (ICNR). Furthermore, the working principle of the MAP probability rule is described in detail showing a close relation to the principles of model-selection methods.