Local fractal dimension of collagen detects increased spatial complexity in fibrosis.

Increase of collagen content and reorganization characterizes fibrosis but quantifying the latter remains challenging. Spatially complex structures are often analyzed via the fractal dimension; however, established methods for calculating this quantity either provide a single dimension for an entire object or a spatially distributed dimension that only considers binary images. These neglect valuable information related to collagen density in images of fibrotic tissue. We sought to develop a fractal analysis that can be applied to 3-dimensional (3D) images of fibrotic tissue. A fractal dimension map for each image was calculated by determining a single fractal dimension for a small area surrounding each image pixel, using fiber thickness as the third dimension. We found that this local fractal dimension increased with age and with progression of fibrosis regardless of collagen content. Our new method of distributed 3D fractal analysis can thus distinguish between changes in collagen content and organization induced by fibrosis.

Visualizing intratumoral injections in lung tumors by endobronchial ultrasound.

Real-time endobronchial ultrasound images are crucial for the accurate placement of the needle in peribronchial lung tumors and lymph nodes for diagnostic sampling. Beyond its role as a diagnostic tool, ultrasound-guided bronchoscopy can also aid the delivery of anti-cancer agents intratumorally, enabling diagnosis, staging, and treatment to occur within the same anesthesia, reducing the patient's burden. However, determining drug retention and distribution in situ remains challenging, albeit pivotal in assessing the success or failure of the therapeutic intervention. We hypothesized that ultrasound images acquired by the bronchoscope during the injection can provide qualitative and quantitative real-time information about drug transport. As a proof-of-concept, we retrospectively analyzed 13 videos of intratumoral cisplatin injections in advanced non-small cell lung cancers. We identified the injection and performed quantitative analysis through image processing and segmentation algorithms and mathematical models in 5 of them. We were able to infer the unlikeliness of a laminar flow through interstitial pores in favor of the emergence of tissue fractures. These data imply that the structural integrity of the tumor is a critical determinant of the ultimate distribution of an intratumorally delivered agent.

A computational modeling approach for dosing endoscopic intratumoral chemotherapy for advanced non-small cell lung cancer.

We recently developed a computational model of cisplatin pharmacodynamics in an endobronchial lung tumor following ultrasound-guided transbronchial needle injection (EBUS-TBNI). The model suggests that it is more efficacious to apportion the cisplatin dose between injections at different sites rather than giving it all in a single central injection, but the model was calibrated only on blood cisplatin data from a single patient. Accordingly, we applied a modified version of our original model in a set of 32 patients undergoing EBUS-TBNI for non-small cell lung cancer (NSCLC). We used the model to predict clinical responses and compared them retrospectively to actual patient outcomes. The model correctly predicted the clinical response in 72% of cases, with 80% accuracy for adenocarcinomas and 62.5% accuracy for squamous-cell lung cancer. We also found a power-law relationship between tumor volume and the minimal dose needed to induce a response, with the power-law exponent depending on the number of injections administered. Our results suggest that current injection strategies may be significantly over- or under-dosing the agent depending on tumor size, and that computational modeling can be a useful planning tool for EBUS-TBNI of cisplatin in lung cancer.

Signal processing to remove spurious contributions to the assessment of respiratory mechanics.

Signal disruptions in small animals during the realization of the Forced Oscillation Technique are a well-known cause of data loss as it leads to non-reliable estimations of the respiratory impedance. In this work, we assessed the effects of removing the disrupted epoch when a 3-seconds input signal composed of one and a half 2-seconds full cycle is used.We tested our hypothesis in 25 SAMR1 mice under different levels of bronchoconstriction due to methacholine administration by iv bolus injections in different doses (15 animals) and by iv continuous infusion in different infusion rates (10 animals). Signal disruptions were computationally simulated as sharp drops in the pressure signal within a short timescale, and signal processing was performed using own developed algorithms.We found that the model goodness of fit worsens when averaging techniques to estimate the input respiratory impedance are not used. However, no statistically significant differences were observed in the comparison between Constant Phase Model parameters of the full 3-s signal and the 2-s non disrupted epoch in all doses or infusion rates for both methacholine delivery strategies.The proposed technique presents reliable outcomes that can reduce animal use in Forced Oscillation Technique realization.

A simple assessment of lung nodule location for reduction in unnecessary invasive procedures.

BACKGROUND: CT screening for lung cancer results in a significant mortality reduction but is complicated by invasive procedures performed for evaluation of the many detected benign nodules. The purpose of this study was to evaluate measures of nodule location within the lung as predictors of malignancy. METHODS: We analyzed images and data from 3,483 participants in the National Lung Screening Trial (NLST). All nodules (4-20 mm) were characterized by 3D geospatial location using a Cartesian coordinate system and evaluated in logistic regression analysis. Model development and probability cutpoint selection was performed in the NLST testing set. The Geospatial test was then validated in the NLST testing set, and subsequently replicated in a new cohort of 147 participants from The Detection of Early Lung Cancer Among Military Personnel (DECAMP) Consortium. RESULTS: The Geospatial Test, consisting of the superior-inferior distance (Z distance), nodule diameter, and radial distance (carina to nodule) performed well in both the NLST validation set (AUC 0.85) and the DECAMP replication cohort (AUC 0.75). A negative Geospatial Test resulted in a less than 2% risk of cancer across all nodule diameters. The Geospatial Test correctly reclassified 19.7% of indeterminate nodules with a diameter over 6mm as benign, while only incorrectly classifying 1% of cancerous nodules as benign. In contrast, the parsimonious Brock Model applied to the same group of nodules correctly reclassified 64.5% of indeterminate nodules as benign but resulted in misclassification of a cancer as benign in 18.2% of the cases. Applying the Geospatial test would result in reducing invasive procedures performed for benign lesions by 11.3% with a low rate of misclassification (1.3%). In contrast, the Brock model applied to the same group of patients results in decreasing invasive procedures for benign lesion by 39.0% but misclassifying 21.1% of cancers as benign. CONCLUSIONS: Utilizing information about geospatial location within the lung improves risk assessment for indeterminate lung nodules and may reduce unnecessary procedures. TRIAL REGISTRATION: NCT00047385, NCT01785342.

Percolation of collagen stress in a random network model of the alveolar wall.

Fibrotic diseases are characterized by progressive and often irreversible scarring of connective tissue in various organs, leading to substantial changes in tissue mechanics largely as a result of alterations in collagen structure. This is particularly important in the lung because its bulk modulus is so critical to the volume changes that take place during breathing. Nevertheless, it remains unclear how fibrotic abnormalities in the mechanical properties of pulmonary connective tissue can be linked to the stiffening of its individual collagen fibers. To address this question, we developed a network model of randomly oriented collagen and elastin fibers to represent pulmonary alveolar wall tissue. We show that the stress-strain behavior of this model arises via the interactions of collagen and elastin fiber networks and is critically dependent on the relative fiber stiffnesses of the individual collagen and elastin fibers themselves. We also show that the progression from linear to nonlinear stress-strain behavior of the model is associated with the percolation of stress across the collagen fiber network, but that the location of the percolation threshold is influenced by the waviness of collagen fibers.

Central airway collapse is related to obesity independent of asthma phenotype.

BACKGROUND AND OBJECTIVE: Late-onset non-allergic asthma in obesity is characterized by an abnormally compliant, collapsible lung periphery; it is not known whether this abnormality exists in proximal airways. We sought to compare collapsibility of central airways between lean and obese individuals with and without asthma. METHODS: A cross-sectional study comparing luminal area and shape (circularity) of the trachea, left mainstem bronchus, right bronchus intermedius and right inferior lobar bronchus at RV and TLC by CT was conducted. RESULTS: In 11 lean controls (BMI: 22.4 (21.5, 23.8) kg/m2 ), 10 lean individuals with asthma (23.6 (22.0, 24.8) kg/m2 ), 10 obese controls (45.5 (40.3, 48.5) kg/m2 ) and 21 obese individuals with asthma (39.2 (35.8, 42.9) kg/m2 ), lumen area and circularity increased significantly with an increase in lung volume from RV to TLC for all four airways (P < 0.05 for all). Changes in area and circularity with lung volume were similar in obese individuals with and without asthma, and both obese groups had severe airway collapse at RV. In multivariate analysis, change in lumen area was related to BMI and change in circularity to waist circumference, but neither was related to asthma diagnosis. CONCLUSION: Excessive collapse of the central airways is related to obesity, and occurs in both obese controls and obese asthma. Increased airway collapse could contribute to ventilation abnormalities in obese individuals particularly at lower lung volumes, and complicate asthma in obese individuals.

Modeling Lung Derecruitment in VILI Due to Fluid-Occlusion: The Role of Emergent Behavior.

Ventilator-induced lung injury (VILI) is driven by the processes of volutrauma and atelectrauma, which can act synergistically to compromise the blood-gas barrier. We have postulated that this synergy arises through a rich-get-richer mechanism whereby atelectrauma causes holes to form in the blood-gas barrier while concomitant volutrauma causes susceptible holes to progressively enlarge as VILI worsens. We previously developed an analytical model based on this idea that accurately predicts the progressive increases in lung elastance seen immediately following a recruitment maneuver as VILI progresses over the course of hours. In the present study we extend this model to account for the rate of change of elastance, due to closure of lung units, in the minutes following a recruitment maneuver. We found that the distribution of unit closing velocities throughout the lung can be described by a power law with an exponent of -2 that matches previously published power laws associated with the dynamics of lung recruitment. Our model thus reveals lung collapse as an example of emergent complex behavior and links the dynamics of altered function in the injured lung to structural damage in a way that explains the mechanisms of injury progression arising from the ongoing stresses and strains applied by mechanical ventilation.

Modeling lung perfusion abnormalities to explain early COVID-19 hypoxemia.

Early stages of the novel coronavirus disease (COVID-19) are associated with silent hypoxia and poor oxygenation despite relatively minor parenchymal involvement. Although speculated that such paradoxical findings may be explained by impaired hypoxic pulmonary vasoconstriction in infected lung regions, no studies have determined whether such extreme degrees of perfusion redistribution are physiologically plausible, and increasing attention is directed towards thrombotic microembolism as the underlying cause of hypoxemia. Herein, a mathematical model demonstrates that the large amount of pulmonary venous admixture observed in patients with early COVID-19 can be reasonably explained by a combination of pulmonary embolism, ventilation-perfusion mismatching in the noninjured lung, and normal perfusion of the relatively small fraction of injured lung. Although underlying perfusion heterogeneity exacerbates existing shunt and ventilation-perfusion mismatch in the model, the reported hypoxemia severity in early COVID-19 patients is not replicated without either extensive perfusion defects, severe ventilation-perfusion mismatch, or hyperperfusion of nonoxygenated regions.

Can Hyperperfusion of Nonaerated Lung Explain COVID-19 Hypoxia?

Early stages of the novel coronavirus disease (COVID-19) have been associated with 'silent hypoxia' and poor oxygenation despite relatively small fractions of afflicted lung. Although it has been speculated that such paradoxical findings may be explained by impairment of hypoxic pulmonary vasoconstriction in infected lungs regions, no studies have confirmed this hypothesis nor determined whether such extreme degrees of perfusion redistribution are physiologically plausible. Here, we present a mathematical model which provides evidence that the extreme amount of pulmonary shunt observed in patients with early COVID-19 is not plausible without hyperperfusion of the relatively small fraction of injured lung, with three-fold increases in regional perfusion to afflicted regions. Although underlying perfusion heterogeneity (e.g., due to gravity or pulmonary emboli) exacerbated existing shunt in the model, the reported severity of hypoxia in early COVID-19 patients could not be replicated without considerable reduction of vascular resistance in nonoxygenated regions.

Evidence for Expanding Invasive Mediastinal Staging for Peripheral T1 Lung Tumors.

BACKGROUND: Guidelines recommend invasive mediastinal staging for patients with non-small cell lung cancer and a "central" tumor. However, there is no consensus definition for central location. As such, the decision to perform invasive staging largely remains on an empirical foundation. RESEARCH QUESTION: Should patients with peripheral T1 lung tumors undergo invasive mediastinal staging? STUDY DESIGN AND METHODS: All participants with a screen-detected cancer with a solid component between 8 and 30 mm were identified from the National Lung Screening Trial. After translation of CT data, cancer location was identified and the X, Y, Z coordinates were determined as well as distance from the main carina. A multivariable logistic regression model was constructed to evaluate for predictors associated with lymph node metastasis. RESULTS: Three hundred thirty-two participants were identified, of which 69 had lymph node involvement (20.8%). Of those with lymph node metastasis, 39.1% were N2. There was no difference in rate of lymph node metastasis based on tumor size (OR, 1.03; P = .248). There was also no statistical difference in rate of lymph node metastasis based on location, either by distance from the carina (OR, 0.99; P = .156) or tumor coordinates (X: P = .180; Y: P = .311; Z: P = .292). When adjusted for age, sex, histology, and smoking history, there was no change in the magnitude of the risk, and tests of significance were not altered. INTERPRETATION: Our data indicate a high rate of N2 metastasis among T1 tumors and no significant relationship between tumor diameter or location. This suggests that patients with small, peripheral lung cancers may benefit from invasive mediastinal staging.

Methacholine dose response curve and acceptability criteria of respiratory mechanics modeling.

Aim: This study aimed to analyze the Constant Phase Model (CPM) Coefficient of Determination (COD) and an index of harmonic distortion ([Formula: see text]) behavior in intravenous methacholine dose response curve. We studied the COD and [Formula: see text] behavior of Control and Lung Inflammation (OVA) groups of mice and we proposed an alternative for moments when the CPM should not be applied. Methods: 9-week female BALB/c mice were studied, 8 of the control group (23.11 ± 1.27 g) and 11 of the lung inflammation group (OVA) (21.45 ± 2.16 g). The COD values were obtained during the respiratory mechanics assessment via Forced Oscillation Technique (FOT) and the [Formula: see text] was estimated a posteriori. Both control and OVA groups were submitted to 4 doses of Methacholine (MCh) protocol. Results: A strong correlation between COD and [Formula: see text] was present at the last two doses (0.3 mg/kg: r = -0.75, p = 0.0013 and 1 mg/kg: r = -0.91; p < 0.0001) in the OVA group. Differences were found in doses of 0.3 mg/kg between control and OVA for the maximum values of Rn (Newtonian Resistance) and G (tissue viscous); and between groups at PBS and doses of 0.03, 0.1 and 0.3 mg/kg for H (Elastance). A similar behavior was observed for the analysis of Area Under the Curve with the exclusion of the 3 first measurements of each dose. However, in this scenario, the comparison with the maximum value presented a higher discriminatory capacity of the parameters associated with the parenchyma. Conclusions: During severe bronchoconstriction there is a strong negative correlation between model goodness of fit and nonlinearities levels, reinforcing that COD is a robust acceptance criterion, whether still simple and easily obtained from the ventilator. We also pointed out the area under the CPM parameters dose response curve is a useful and can be used as a complementary analysis to peak comparison following bolus injections of methacholine.

Cisplatin Pharmacodynamics Following Endobronchial Ultrasound-Guided Transbronchial Needle Injection into Lung Tumors.

Intratumoral delivery of cisplatin by endobronchial ultrasound-guided transbronchial needle injection (EBUS-TBNI) has recently emerged as a therapy for treating peribronchial lung cancers. It remains unclear, however, where best to inject drug into a tumor, and at how many sites, so current cisplatin delivery strategies remain empirical. Motivated by the need to put EBUS-TBNI treatment of lung cancer on a more objective footing, we developed a computational model of cisplatin pharmacodynamics following EBUS-TBNI. The model accounts for diffusion of cisplatin within and between the intracellular and extracellular spaces of a tumor, as well as clearance of cisplatin from the tumor via the vasculature and clearance from the body via the kidneys. We matched the tumor model geometry to that determined from a thoracic CT scan of a patient with lung cancer. The model was calibrated by fitting its predictions of cisplatin blood concentration versus time to measurements made up to 2 hrs following EBUS-TBNI of cisplatin into the patient's lung tumor. This gave a value for the systemic volume of distribution for cisplatin of 12.2 L and a rate constant of clearance from the tumor into the systemic compartment of 1.46 × 10-4 s-1. Our model indicates that the minimal dose required to kill all cancerous cells in a lung tumor can be reduced by roughly 3 orders of magnitude if the cisplatin is apportioned between 5 optimally spaced locations throughout the tumor rather than given as a single bolus to the tumor center. Our findings suggest that optimizing the number and location of EBUS-TBNI sites has a dramatic effect on the dose of cisplatin required for efficacious treatment of lung cancer.

Linking Physiological Biomarkers of Ventilator-Induced Lung Injury to a Rich-Get-Richer Mechanism of Injury Progression.

Mechanical ventilation is a crucial tool in the management of acute respiratory distress syndrome, yet it may itself also further damage the lung in a phenomenon known as ventilator-induced lung injury (VILI). We have previously shown in mice that volutrauma and atelectrauma act synergistically to cause VILI. We have also postulated that this synergy arises because of a rich-get-richer mechanism in which repetitive lung recruitment generates initial small holes in the blood-gas barrier which are then expanded by over-distension in a manner that favors large holes over small ones. In order to understand the causal link between this process and the derangements in lung mechanics associated with VILI, we developed a mathematical model that incorporates both atelectrauma and volutrauma to predict how the propensity of the lung to derecruit depends on the accumulation of plasma-derived fluid and proteins in the airspaces. We found that the model accurately predicts derecruitment in mice with experimentally induced VILI.