Deep phenotyping, including quantitative ciliary beating parameters, and extensive genotyping in primary ciliary dyskinesia.

BACKGROUND: Primary ciliary dyskinesia (PCD) is a rare genetic disorder resulting in abnormal ciliary motility/structure, extremely heterogeneous at genetic and ultrastructural levels. We aimed, in light of extensive genotyping, to identify specific and quantitative ciliary beating anomalies, according to the ultrastructural phenotype. METHODS: We prospectively included 75 patients with PCD exhibiting the main five ultrastructural phenotypes (n=15/group), screened all corresponding PCD genes and measured quantitative beating parameters by high-speed video-microscopy (HSV). RESULTS: Sixty-eight (91%) patients carried biallelic mutations. Combined outer/inner dynein arms (ODA/IDA) defect induces total ciliary immotility, regardless of the gene involved. ODA defect induces a residual beating with dramatically low ciliary beat frequency (CBF) related to increased recovery stroke and pause durations, especially in case of DNAI1 mutations. IDA defect with microtubular disorganisation induces a low percentage of beating cilia with decreased beating angle and, in case of CCDC39 mutations, a relatively conserved mean CBF with a high maximal CBF. Central complex defect induces nearly normal beating parameters, regardless of the gene involved, and a gyrating motion in a minority of ciliated edges, especially in case of RSPH1 mutations. PCD with normal ultrastructure exhibits heterogeneous HSV values, but mostly an increased CBF with an extremely high maximal CBF. CONCLUSION: Quantitative HSV analysis in PCD objectives beating anomalies associated with specific ciliary ultrastructures and genotypes. It represents a promising approach to guide the molecular analyses towards the best candidate gene(s) to be analysed or to assess the pathogenicity of the numerous sequence variants identified by next-generation-sequencing.

Pathogenesis of chronic rhinosinusitis with nasal polyps: role of IL-6 in airway epithelial cell dysfunction.

BACKGROUND: Chronic rhinosinusitis with nasal polyps (CRSwNP) is characterized by an alteration in airway epithelial cell functions including barrier function, wound repair mechanisms, mucociliary clearance. The mechanisms leading to epithelial cell dysfunction in nasal polyps (NPs) remain poorly understood. Our hypothesis was that among the inflammatory cytokines involved in NPs, IL-6 could alter epithelial repair mechanisms and mucociliary clearance. The aim of this study was to evaluate the in vitro effects of IL-6 on epithelial repair mechanisms in a wound repair model and on ciliary beating in primary cultures of Human Nasal Epithelial Cells (HNEC). METHODS: Primary cultures of HNEC taken from 38 patients during surgical procedures for CRSwNP were used in an in vitro model of wound healing. Effects of increasing concentrations of IL-6 (1 ng/mL, 10 ng/mL, and 100 ng/mL) and other ILs (IL-5, IL-9, IL-10) on wound closure kinetics were compared to cultures without IL-modulation. After wound closure, the differentiation process was characterized under basal conditions and after IL supplementation using cytokeratin-14, MUC5AC, and ßIV tubulin as immunomarkers of basal, mucus, and ciliated cells, respectively. The ciliated edges of primary cultures were analyzed on IL-6 modulation by digital high-speed video-microscopy to measure: ciliary beating frequency (CBF), ciliary length, relative ciliary density, metachronal wavelength and the ciliary beating efficiency index. RESULTS: Our results showed that: (i) IL-6 accelerated airway wound repair in vitro, with a dose-response effect whereas no effect was observed after other ILs-stimulation. After 24 h, 79% of wounded wells with IL6-100 were fully repaired, vs 46% in the IL6-10 group, 28% in the IL6-1 group and 15% in the control group; (ii) specific migration analyses of closed wound at late repair stage (Day 12) showed IL-6 had the highest migration compared with other ILs (iii) The study of the IL-6 effect on ciliary function showed that CBF and metachronal wave increased but without significant modifications of ciliary density, length of cilia and efficiency index. CONCLUSION: The up-regulated epithelial cell proliferation observed in polyps could be induced by IL-6 in the case of prior epithelial damage. IL-6 could be a major cytokine in NP physiopathology.

Surfactant delivery in rat lungs: Comparing 3D geometrical simulation model with experimental instillation.

Surfactant Replacement Therapy (SRT), which involves instillation of a liquid-surfactant mixture directly into the lung airway tree, is a major therapeutic treatment in neonatal patients with respiratory distress syndrome (RDS). This procedure has proved to be remarkably effective in premature newborns, inducing a five-fold decrease of mortality in the past 35 years. Disappointingly, its use in adults for treating acute respiratory distress syndrome (ARDS) experienced initial success followed by failures. Our recently developed numerical model has demonstrated that transition from success to failure of SRT in adults could, in fact, have a fluid mechanical origin that is potentially reversible. Here, we present the first numerical simulations of surfactant delivery into a realistic asymmetric conducting airway tree of the rat lung and compare them with experimental results. The roles of dose volume (VD), flow rate, and multiple aliquot delivery are investigated. We find that our simulations of surfactant delivery in rat lungs are in good agreement with our experimental data. In particular, we show that the monopodial architecture of the rat airway tree plays a major role in surfactant delivery, contributing to the poor homogeneity of the end distribution of surfactant. In addition, we observe that increasing VD increases the amount of surfactant delivered to the acini after losing a portion to coating the involved airways, the coating cost volume, VCC. Finally, we quantitatively assess the improvement resulting from a multiple aliquot delivery, a method sometimes employed clinically, and find that a much larger fraction of surfactant reaches the alveolar regions in this case. This is the first direct qualitative and quantitative comparison of our numerical model with experimental studies, which enhances our previous predictions in adults and neonates while providing a tool for predicting, engineering, and optimizing patient-specific surfactant delivery in complex situations.

A new index for characterizing micro-bead motion in a flow induced by ciliary beating: Part II, modeling.

Mucociliary clearance is one of the major lines of defense of the human respiratory system. The mucus layer coating the airways is constantly moved along and out of the lung by the activity of motile cilia, expelling at the same time particles trapped in it. The efficiency of the cilia motion can experimentally be assessed by measuring the velocity of micro-beads traveling through the fluid surrounding the cilia. Here we present a mathematical model of the fluid flow and of the micro-beads motion. The coordinated movement of the ciliated edge is represented as a continuous envelope imposing a periodic moving velocity boundary condition on the surrounding fluid. Vanishing velocity and vanishing shear stress boundary conditions are applied to the fluid at a finite distance above the ciliated edge. The flow field is expanded in powers of the amplitude of the individual cilium movement. It is found that the continuous component of the horizontal velocity at the ciliated edge generates a 2D fluid velocity field with a parabolic profile in the vertical direction, in agreement with the experimental measurements. Conversely, we show than this model can be used to extract microscopic properties of the cilia motion by extrapolating the micro-bead velocity measurement at the ciliated edge. Finally, we derive from these measurements a scalar index providing a direct assessment of the cilia beating efficiency. This index can easily be measured in patients without any modification of the current clinical procedures.

A new index for characterizing micro-bead motion in a flow induced by ciliary beating: Part I, experimental analysis.

Mucociliary clearance is one of the major lines of defense of the respiratory system. The mucus layer coating the pulmonary airways is moved along and out of the lung by the activity of motile cilia, thus expelling the particles trapped in it. Here we compare ex vivo measurements of a Newtonian flow induced by cilia beating (using micro-beads as tracers) and a mathematical model of this fluid flow, presented in greater detail in a second companion article. Samples of nasal epithelial cells placed in water are recorded by high-speed video-microscopy and ciliary beat pattern is inferred. Automatic tracking of micro-beads, used as markers of the flow generated by cilia motion, enables us also to assess the velocity profile as a function of the distance above the cilia. This profile is shown to be essentially parabolic. The obtained experimental data are used to feed a 2D mathematical and numerical model of the coupling between cilia, fluid, and micro-bead motion. From the model and the experimental measurements, the shear stress exerted by the cilia is deduced. Finally, this shear stress, which can easily be measured in the clinical setting, is proposed as a new index for characterizing the efficiency of ciliary beating.

Characterisation of cellular adhesion reinforcement by multiple bond force spectroscopy in alveolar epithelial cells.

BACKGROUND INFORMATION: Integrin-mediated adhesion is a key process by which cells physically connect with their environment, and express sensitivity and adaptation through mechanotransduction. A critical step of cell adhesion is the formation of the first bonds which individually generate weak contacts (∼tens pN) but can sustain thousand times higher forces (∼tens nN) when associated. RESULTS: We propose an experimental validation by multiple bond force spectroscopy (MFS) of a stochastic model predicting adhesion reinforcement permitted by non-cooperative, multiple bonds on which force is homogeneously distributed (called parallel bond configuration). To do so, spherical probes (diameter: 6.6 µm), specifically coated by RGD-peptide to bind integrins, are used to statically indent and homogenously stretch the multiple bonds created for short contact times (2 s) between the bead and the surface of epithelial cells (A549). Using different separation speeds (v = 2, 5, 10 µm/s) and measuring cellular Young's modulus as well as the local stiffness preceding local rupture events, we obtain cell-by-cell the effective loading rates both at the global cell level and at the local level of individual constitutive bonds. Local rupture forces are in the range: f*=60-115 pN , whereas global rupture (detachment) forces reach F*=0.8-1.7 nN . Global and local rupture forces both exhibit linear dependencies with the effective loading rate, the slopes of these two linear relationships providing an estimate of the number of independent integrin bonds constituting the tested multiple bond structure (∼12). CONCLUSIONS: The MFS method enables to validate the reinforcement of integrin-mediated adhesion induced by the multiple bond configuration in which force is homogeneously distributed amongst parallel bonds. Local rupture events observed in the course of a spectroscopy manoeuver (MFS) lead to rupture force values considered in the literature as single-integrin bonds. SIGNIFICANCE: Adhesion reinforcement permitted by the parallel multiple bond association is particularly challenging to verify for two reasons: first, it is difficult to control precisely the direction of forces experimentally, and second, because both global and local bond rupture forces depend on the effective loading rate applied to the bond. Here, we propose an integrin-specific MFS method capable of detecting bond number and characterising bond configuration and its impact on adhesion strength.

Exposure to Bordetella pertussis adenylate cyclase toxin affects integrin-mediated adhesion and mechanics in alveolar epithelial cells.

BACKGROUND INFORMATION: The adenylate cyclase (CyaA) toxin is a major virulent factor of Bordetella pertussis, the causative agent of whooping cough. CyaA toxin is able to invade eukaryotic cells where it produces high levels of cyclic adenosine monophosphate (cAMP) affecting cellular physiology. Whether CyaA toxin can modulate cell matrix adhesion and mechanics of infected cells remains largely unknown. RESULTS: In this study, we use a recently proposed multiple bond force spectroscopy (MFS) with an atomic force microscope to assess the early phase of cell adhesion (maximal detachment and local rupture forces) and cell rigidity (Young's modulus) in alveolar epithelial cells (A549) for toxin exposure <1 h. At 30 min of exposure, CyaA toxin has a minimal effect on cell viability (>95%) at CyaA concentration of 0.5 nM, but a significant effect (≈81%) at 10 nM. MFS performed on A549 for three different concentrations (0.5, 5 and 10 nM) demonstrates that CyaA toxin significantly affects both cell adhesion (detachment forces are decreased) and cell mechanics (Young's modulus is increased). CyaA toxin (at 0.5 nM) assessed at three indentation/retraction speeds (2, 5 and 10 µm/s) significantly affects global detachment forces, local rupture events and Young modulus compared with control conditions, while an enzymatically inactive variant CyaAE5 has no effect. These results reveal the loading rate dependence of the multiple bonds newly formed between the cell and integrin-specific coated probe as well as the individual bond kinetics which are only slightly affected by the patho-physiological dose of CyaA toxin. Finally, theory of multiple bond force rupture enables us to deduce the bond number N which is reduced by a factor of 2 upon CyaA exposure (N ≈ 6 versus N ≈ 12 in control conditions). CONCLUSIONS: MFS measurements demonstrate that adhesion and mechanical properties of A549 are deeply affected by exposure to the CyaA toxin but not to an enzymatically inactive variant. This indicates that the alteration of cell mechanics triggered by CyaA is a consequence of the increase in intracellular cAMP in these target cells. SIGNIFICANCE: These results suggest that mechanical and adhesion properties of the cells appear as pertinent markers of cytotoxicity of CyaA toxin.

Three-dimensional model of surfactant replacement therapy.

Surfactant replacement therapy (SRT) involves instillation of a liquid-surfactant mixture directly into the lung airway tree. It is widely successful for treating surfactant deficiency in premature neonates who develop neonatal respiratory distress syndrome (NRDS). However, when applied to adults with acute respiratory distress syndrome (ARDS), early successes were followed by failures. This unexpected and puzzling situation is a vexing issue in the pulmonary community. A pressing question is whether the instilled surfactant mixture actually reaches the adult alveoli/acinus in therapeutic amounts. In this study, to our knowledge, we present the first mathematical model of SRT in a 3D lung structure to provide insight into answering this and other questions. The delivery is computed from fluid mechanical principals for 3D models of the lung airway tree for neonates and adults. A liquid plug propagates through the tree from forced inspiration. In two separate modeling steps, the plug deposits a coating film on the airway wall and then splits unevenly at the bifurcation due to gravity. The model generates 3D images of the resulting acinar distribution and calculates two global indexes, efficiency and homogeneity. Simulating published procedural methods, we show the neonatal lung is a well-mixed compartment, whereas the adult lung is not. The earlier, successful adult SRT studies show comparatively good index values implying adequate delivery. The later, failed studies used different protocols resulting in very low values of both indexes, consistent with inadequate acinar delivery. Reasons for these differences and the evolution of failure from success are outlined and potential remedies discussed.

Effective Confining Potential of Quantum States in Disordered Media.

The amplitude of localized quantum states in random or disordered media may exhibit long-range exponential decay. We present here a theory that unveils the existence of an effective potential which finely governs the confinement of these states. In this picture, the boundaries of the localization subregions for low energy eigenfunctions correspond to the barriers of this effective potential, and the long-range exponential decay characteristic of Anderson localization is explained as the consequence of multiple tunneling in the dense network of barriers created by this effective potential. Finally, we show that Weyl's formula based on this potential turns out to be a remarkable approximation of the density of states for a large variety of one-dimensional systems, periodic or random.

One Single Static Measurement Predicts Wave Localization in Complex Structures.

A recent theoretical breakthrough has brought a new tool, called the localization landscape, for predicting the localization regions of vibration modes in complex or disordered systems. Here, we report on the first experiment which measures the localization landscape and demonstrates its predictive power. Holographic measurement of the static deformation under uniform load of a thin plate with complex geometry provides direct access to the landscape function. When put in vibration, this system shows modes precisely confined within the subregions delineated by the landscape function. Also the maxima of this function match the measured eigenfrequencies, while the minima of the valley network gives the frequencies at which modes become extended. This approach fully characterizes the low frequency spectrum of a complex structure from a single static measurement. It paves the way for controlling and engineering eigenmodes in any vibratory system, especially where a structural or microscopic description is not accessible.

Studies of the B-Z transition of DNA: The temperature dependence of the free-energy difference, the composition of the counterion sheath in mixed salt, and the preparation of a sample of the 5'-d[T-(m(5) C-G)12 -T] duplex in pure B-DNA or Z-DNA form.

It is often envisioned that cations might coordinate at specific sites of nucleic acids and play an important structural role, for instance in the transition between B-DNA and Z-DNA. However, nucleic acid models explicitly devoid of specific sites may also exhibit features previously considered as evidence for specific binding. Such is the case of the "composite cylinder" (or CC) model which spreads out localized features of DNA structure and charge by cylindrical averaging, while sustaining the main difference between the B and Z structures, namely the better immersion of the B-DNA phosphodiester charges in the solution. Here, we analyze the non-electrostatic component of the free-energy difference between B-DNA and Z-DNA. We also compute the composition of the counterion sheath in a wide range of mixed-salt solutions and of temperatures: in contrast with the large difference of composition between the B-DNA and Z-DNA forms, the temperature dependence of sheath composition, previously unknown, is very weak. In order to validate the model, the mixed-salt predictions should be compared to experiment. We design a procedure for future measurements of the sheath composition based on Anomalous Small-Angle X-ray Scattering and complemented by (31) P NMR. With due consideration for the kinetics of the B-Z transition and for the capacity of generating at will the B or Z form in a single sample, the 5'-d[T-(m(5) C-G)12 -T] 26-mer emerges as a most suitable oligonucleotide for this study. Finally, the application of the finite element method to the resolution of the Poisson-Boltzmann equation is described in detail. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 369-384, 2016.

A Macroscopic Model for Simulating the Mucociliary Clearance in a Bronchial Bifurcation: The Role of Surface Tension.

The mucociliary clearance in the bronchial tree is the main mechanism by which the lungs clear themselves of deposited particulate matter. In this work, a macroscopic model of the clearance mechanism is proposed. Lubrication theory is applied for thin films with both surface tension effects and a moving wall boundary. The flow field is computed by the use of a finite-volume scheme on an unstructured grid that replicates a bronchial bifurcation. The carina in bronchial bifurcations is of special interest because it is a location of increased deposition of inhaled particles. In this study, the mucus flow is computed for different values of the surface tension. It is found that a minimal surface tension is necessary for efficiently removing the mucus while maintaining the mucus film thickness at physiological levels.

Universal mechanism for Anderson and weak localization.

Localization of stationary waves occurs in a large variety of vibrating systems, whether mechanical, acoustical, optical, or quantum. It is induced by the presence of an inhomogeneous medium, a complex geometry, or a quenched disorder. One of its most striking and famous manifestations is Anderson localization, responsible for instance for the metal-insulator transition in disordered alloys. Yet, despite an enormous body of related literature, a clear and unified picture of localization is still to be found, as well as the exact relationship between its many manifestations. In this paper, we demonstrate that both Anderson and weak localizations originate from the same universal mechanism, acting on any type of vibration, in any dimension, and for any domain shape. This mechanism partitions the system into weakly coupled subregions. The boundaries of these subregions correspond to the valleys of a hidden landscape that emerges from the interplay between the wave operator and the system geometry. The height of the landscape along its valleys determines the strength of the coupling between the subregions. The landscape and its impact on localization can be determined rigorously by solving one special boundary problem. This theory allows one to predict the localization properties, the confining regions, and to estimate the energy of the vibrational eigenmodes through the properties of one geometrical object. In particular, Anderson localization can be understood as a special case of weak localization in a very rough landscape.

Optimisations and evolution of the mammalian respiratory system : A suggestion of possible gene sharing in evolution.

The respiratory system of mammalians is made of two successive branched structures with different physiological functions. The upper structure, or bronchial tree, is a fluid transportation system made of approximately 15 generations of bifurcations leading to the order of about 2(15) = 30, 000 terminal bronchioles with a diameter of approximately 0.5mm in the human lung. The branching pattern continues up to generation 23 but the structure and function of each of the subsequent structures, called acini, is different. Each acinus consists in a branched system of ducts surrounded by alveoli and plays the role of a diffusion cell where oxygen and carbon dioxide are exchanged with blood across the alveolar membrane. We show here that the bronchial tree simultaneously presents several different optimal properties. It is first energy efficient, second, it is space filling and third it is also "rapid". This physically based multi-optimality suggests that, in the course of evolution, an organ selected against one criterion could have been used later for a totally different purpose. For example, once selected for its energetic efficiency for the transport of a viscous fluid like blood, the same genetic material could have been used for its optimized rapidity. This would have allowed the emergence of atmospheric respiration made of inspiration-expiration cycles. For this phenomenon to exist, rapidity is essential as fresh air has to reach the gas exchange organs, the pulmonary acini, before the beginning of expiration. We finally show that the pulmonary acinus is optimized in the sense that the acinus morphology is directly related to the notion of a "best possible" extraction of entropic energy by a diffusion exchanger that has to feed oxygen efficiently from air to blood across a membrane of finite permeability.

The Electronic Disorder Landscape of Mixed Halide Perovskites.

Band gap tunability of lead mixed halide perovskites makes them promising candidates for various applications in optoelectronics. Here we use the localization landscape theory to reveal that the static disorder due to iodide:bromide compositional alloying contributes at most 3 meV to the Urbach energy. Our modeling reveals that the reason for this small contribution is due to the small effective masses in perovskites, resulting in a natural length scale of around 20 nm for the "effective confining potential" for electrons and holes, with short-range potential fluctuations smoothed out. The increase in Urbach energy across the compositional range agrees well with our optical absorption measurements. We model systems of sizes up to 80 nm in three dimensions, allowing us to accurately reproduce the experimentally observed absorption spectra of perovskites with halide segregation. Our results suggest that we should look beyond static contribution and focus on the dynamic temperature dependent contribution to the Urbach energy.

Optimal branching asymmetry of hydrodynamic pulsatile trees.

Most of the studies on optimal transport are done for steady state regime conditions. Yet, there exists numerous examples in living systems where supply tree networks have to deliver products in a limited time due to the pulsatile character of the flow, as it is the case for mammalian respiration. We report here that introducing a systematic branching asymmetry allows the tree to reduce the average delivery time of the products. It simultaneously increases its robustness against the inevitable variability of sizes related to morphogenesis. We then apply this approach to the human tracheobronchial tree. We show that in this case all extremities are supplied with fresh air, provided that the asymmetry is smaller than a critical threshold which happens to match the asymmetry measured in the human lung. This could indicate that the structure is tuned at the maximum asymmetry level that allows the lung to feed all terminal units with fresh air.

Functional and structural consequences of epithelial cell invasion by Bordetella pertussis adenylate cyclase toxin.

Bordetella pertussis, the causative agent of whopping cough, produces an adenylate cyclase toxin (CyaA) that plays a key role in the host colonization by targeting innate immune cells which express CD11b/CD18, the cellular receptor of CyaA. CyaA is also able to invade non-phagocytic cells, via a unique entry pathway consisting in a direct translocation of its catalytic domain across the cytoplasmic membrane of the cells. Within the cells, CyaA is activated by calmodulin to produce high levels of cyclic adenosine monophosphate (cAMP) and alter cellular physiology. In this study, we explored the effects of CyaA toxin on the cellular and molecular structure remodeling of A549 alveolar epithelial cells. Using classical imaging techniques, biochemical and functional tests, as well as advanced cell mechanics method, we quantify the structural and functional consequences of the massive increase of intracellular cyclic AMP induced by the toxin: cell shape rounding associated to adhesion weakening process, actin structure remodeling for the cortical and dense components, increase in cytoskeleton stiffness, and inhibition of migration and repair. We also show that, at low concentrations (0.5 nM), CyaA could significantly impair the migration and wound healing capacities of the intoxicated alveolar epithelial cells. As such concentrations might be reached locally during B. pertussis infection, our results suggest that the CyaA, beyond its major role in disabling innate immune cells, might also contribute to the local alteration of the epithelial barrier of the respiratory tract, a hallmark of pertussis.

Strong localization induced by one clamped point in thin plate vibrations.

We discover a strong localization of flexural (bi-Laplacian) waves in rigid thin plates. We show that clamping just one point inside such a plate not only perturbs its spectral properties, but essentially divides the plate into two independently vibrating regions. This effect progressively appears when increasing the plate eccentricity. Such a localization is qualitatively and quantitatively different from the results known for the Laplacian waves in domains of irregular boundary. It would allow us to control the confinement of mechanical vibrations in rigid plates and of eddies in the slow Stokes flow.

Universality of fold-encoded localized vibrations in enzymes.

Enzymes speed up biochemical reactions at the core of life by as much as 15 orders of magnitude. Yet, despite considerable advances, the fine dynamical determinants at the microscopic level of their catalytic proficiency are still elusive. In this work, we use a powerful mathematical approach to show that rate-promoting vibrations in the picosecond range, specifically encoded in the 3D protein structure, are localized vibrations optimally coupled to the chemical reaction coordinates at the active site. Remarkably, our theory also exposes an hithertho unknown deep connection between the unique localization fingerprint and a distinct partition of the 3D fold into independent, foldspanning subdomains that govern long-range communication. The universality of these features is demonstrated on a pool of more than 900 enzyme structures, comprising a total of more than 10,000 experimentally annotated catalytic sites. Our theory provides a unified microscopic rationale for the subtle structure-dynamics-function link in proteins.