Progenitor identification and SARS-CoV-2 infection in human distal lung organoids.

The distal lung contains terminal bronchioles and alveoli that facilitate gas exchange. Three-dimensional in vitro human distal lung culture systems would strongly facilitate the investigation of pathologies such as interstitial lung disease, cancer and coronavirus disease 2019 (COVID-19) pneumonia caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Here we describe the development of a long-term feeder-free, chemically defined culture system for distal lung progenitors as organoids derived from single adult human alveolar epithelial type II (AT2) or KRT5+ basal cells. AT2 organoids were able to differentiate into AT1 cells, and basal cell organoids developed lumens lined with differentiated club and ciliated cells. Single-cell analysis of KRT5+ cells in basal organoids revealed a distinct population of ITGA6+ITGB4+ mitotic cells, whose offspring further segregated into a TNFRSF12Ahi subfraction that comprised about ten per cent of KRT5+ basal cells. This subpopulation formed clusters within terminal bronchioles and exhibited enriched clonogenic organoid growth activity. We created distal lung organoids with apical-out polarity to present ACE2 on the exposed external surface, facilitating infection of AT2 and basal cultures with SARS-CoV-2 and identifying club cells as a target population. This long-term, feeder-free culture of human distal lung organoids, coupled with single-cell analysis, identifies functional heterogeneity among basal cells and establishes a facile in vitro organoid model of human distal lung infections, including COVID-19-associated pneumonia.

Progenitor identification and SARS-CoV-2 infection in long-term human distal lung organoid cultures.

The distal lung contains terminal bronchioles and alveoli that facilitate gas exchange and is affected by disorders including interstitial lung disease, cancer, and SARS-CoV-2-associated COVID-19 pneumonia. Investigations of these localized pathologies have been hindered by a lack of 3D in vitro human distal lung culture systems. Further, human distal lung stem cell identification has been impaired by quiescence, anatomic divergence from mouse and lack of lineage tracing and clonogenic culture. Here, we developed robust feeder-free, chemically-defined culture of distal human lung progenitors as organoids derived clonally from single adult human alveolar epithelial type II (AT2) or KRT5 + basal cells. AT2 organoids exhibited AT1 transdifferentiation potential, while basal cell organoids progressively developed lumens lined by differentiated club and ciliated cells. Organoids consisting solely of club cells were not observed. Upon single cell RNA-sequencing (scRNA-seq), alveolar organoids were composed of proliferative AT2 cells; however, basal organoid KRT5 + cells contained a distinct ITGA6 + ITGB4 + mitotic population whose proliferation segregated to a TNFRSF12A hi subfraction. Clonogenic organoid growth was markedly enriched within the TNFRSF12A hi subset of FACS-purified ITGA6 + ITGB4 + basal cells from human lung or derivative organoids. In vivo, TNFRSF12A + cells comprised ~10% of KRT5 + basal cells and resided in clusters within terminal bronchioles. To model COVID-19 distal lung disease, we everted the polarity of basal and alveolar organoids to rapidly relocate differentiated club and ciliated cells from the organoid lumen to the exterior surface, thus displaying the SARS-CoV-2 receptor ACE2 on the outwardly-facing apical aspect. Accordingly, basal and AT2 apical-out organoids were infected by SARS-CoV-2, identifying club cells as a novel target population. This long-term, feeder-free organoid culture of human distal lung alveolar and basal stem cells, coupled with single cell analysis, identifies unsuspected basal cell functional heterogeneity and exemplifies progenitor identification within a slowly proliferating human tissue. Further, our studies establish a facile in vitro organoid model for human distal lung infectious diseases including COVID-19-associated pneumonia.

Advances in Proximity Ligation in situ Hybridization (PLISH).

Understanding tissues in the context of development, maintenance and disease requires determining the molecular profiles of individual cells within their native in vivo spatial context. We developed a Proximity Ligation in situ Hybridization technology (PLISH) that enables quantitative measurement of single cell gene expression in intact tissues, which we have now updated. By recording spatial information for every profiled cell, PLISH enables retrospective mapping of distinct cell classes and inference of their in vivo interactions. PLISH has high sensitivity, specificity and signal to noise ratio. It is also rapid, scalable, and does not require expertise in molecular biology so it can be easily adopted by basic and clinical researchers.

Recording and Analyzing Nucleic Acid Distance Distributions with X-Ray Scattering Interferometry (XSI).

Most structural techniques provide averaged information or information about a single predominant conformational state. However, biological macromolecules typically function through series of conformations. Therefore, a complete understanding of macromolecular structures requires knowledge of the ensembles that represent probabilities on a conformational free energy landscape. Here we describe an emerging approach, X-ray scattering interferometry (XSI), a method that provides instantaneous distance distributions for molecules in solution. XSI uses gold nanocrystal labels site-specifically attached to a macromolecule and measures the scattering interference from pairs of heavy metal labels. The recorded signal can directly be transformed into a distance distribution between the two probes. We describe the underlying concepts, present a detailed protocol for preparing samples and recording XSI data, and provide a custom-written graphical user interface to facilitate XSI data analysis. © 2018 by John Wiley & Sons, Inc.

Gold nanocrystal labels provide a sequence-to-3D structure map in SAXS reconstructions.

Small-angle x-ray scattering (SAXS) is a powerful technique to probe the structure of biological macromolecules and their complexes under virtually arbitrary solution conditions, without the need for crystallization. While it is possible to reconstruct molecular shapes from SAXS data ab initio, the resulting electron density maps have a resolution of ~1 nm and are often insufficient to reliably assign secondary structure elements or domains. We show that SAXS data of gold-labeled samples significantly enhance the information content of SAXS measurements, allowing the unambiguous assignment of macromolecular sequence motifs to specific locations within a SAXS structure. We first demonstrate our approach for site-specifically internally and end-labeled DNA and an RNA motif. In addition, we present a protocol for highly uniform and site-specific labeling of proteins with small (~1.4 nm diameter) gold particles and apply our method to the signaling protein calmodulin. In all cases, the position of the small gold probes can be reliably identified in low-resolution electron density maps. Enhancing low-resolution measurements by site-selective gold labeling provides an attractive approach to aid modeling of a large range of macromolecular systems.

Single-cell Wnt signaling niches maintain stemness of alveolar type 2 cells.

Alveoli, the lung's respiratory units, are tiny sacs where oxygen enters the bloodstream. They are lined by flat alveolar type 1 (AT1) cells, which mediate gas exchange, and AT2 cells, which secrete surfactant. Rare AT2s also function as alveolar stem cells. We show that AT2 lung stem cells display active Wnt signaling, and many of them are near single, Wnt-expressing fibroblasts. Blocking Wnt secretion depletes these stem cells. Daughter cells leaving the Wnt niche transdifferentiate into AT1s: Maintaining Wnt signaling prevents transdifferentiation, whereas abrogating Wnt signaling promotes it. Injury induces AT2 autocrine Wnts, recruiting "bulk" AT2s as progenitors. Thus, individual AT2 stem cells reside in single-cell fibroblast niches providing juxtacrine Wnts that maintain them, whereas injury induces autocrine Wnts that transiently expand the progenitor pool. This simple niche maintains the gas exchange surface and is coopted in cancer.

Automated cell-type classification in intact tissues by single-cell molecular profiling.

A major challenge in biology is identifying distinct cell classes and mapping their interactions in vivo. Tissue-dissociative technologies enable deep single cell molecular profiling but do not provide spatial information. We developed a proximity ligation in situ hybridization technology (PLISH) with exceptional signal strength, specificity, and sensitivity in tissue. Multiplexed data sets can be acquired using barcoded probes and rapid label-image-erase cycles, with automated calculation of single cell profiles, enabling clustering and anatomical re-mapping of cells. We apply PLISH to expression profile ~2900 cells in intact mouse lung, which identifies and localizes known cell types, including rare ones. Unsupervised classification of the cells indicates differential expression of 'housekeeping' genes between cell types, and re-mapping of two sub-classes of Club cells highlights their segregated spatial domains in terminal airways. By enabling single cell profiling of various RNA species in situ, PLISH can impact many areas of basic and medical research.

Ultrasensitive optical imaging with lanthanide lumiphores.

In principle, the millisecond emission lifetimes of lanthanide chelates should enable their ultrasensitive detection in biological systems by time-resolved optical microscopy. In practice, however, lanthanide imaging techniques have provided no better sensitivity than conventional fluorescence microscopy. Here, we identified three fundamental problems that have impeded lanthanide microscopy: low photon flux, inefficient excitation, and optics-derived background luminescence. We overcame these limitations with a new lanthanide imaging modality, transreflected illumination with luminescence resonance energy transfer (trLRET), which increases the time-integrated signal intensities of lanthanide lumiphores by 170-fold and the signal-to-background ratios by 75-fold. We demonstrate that trLRET provides at least an order-of-magnitude increase in detection sensitivity over that of conventional epifluorescence microscopy when used to visualize endogenous protein expression in zebrafish embryos. We also show that trLRET can be used to optically detect molecular interactions in vivo. trLRET promises to unlock the full potential of lanthanide lumiphores for ultrasensitive, autofluorescence-free biological imaging.

Determination of the conformational ensemble of the TAR RNA by X-ray scattering interferometry.

The conformational ensembles of structured RNA's are crucial for biological function, but they remain difficult to elucidate experimentally. We demonstrate with HIV-1 TAR RNA that X-ray scattering interferometry (XSI) can be used to determine RNA conformational ensembles. X-ray scattering interferometry (XSI) is based on site-specifically labeling RNA with pairs of heavy atom probes, and precisely measuring the distribution of inter-probe distances that arise from a heterogeneous mixture of RNA solution structures. We show that the XSI-based model of the TAR RNA ensemble closely resembles an independent model derived from NMR-RDC data. Further, we show how the TAR RNA ensemble changes shape at different salt concentrations. Finally, we demonstrate that a single hybrid model of the TAR RNA ensemble simultaneously fits both the XSI and NMR-RDC data set and show that XSI can be combined with NMR-RDC to further improve the quality of the determined ensemble. The results suggest that XSI-RNA will be a powerful approach for characterizing the solution conformational ensembles of RNAs and RNA-protein complexes under diverse solution conditions.

Directed Chemical Evolution with an Outsized Genetic Code.

The first demonstration that macromolecules could be evolved in a test tube was reported twenty-five years ago. That breakthrough meant that billions of years of chance discovery and refinement could be compressed into a few weeks, and provided a powerful tool that now dominates all aspects of protein engineering. A challenge has been to extend this scientific advance into synthetic chemical space: to enable the directed evolution of abiotic molecules. The problem has been tackled in many ways. These include expanding the natural genetic code to include unnatural amino acids, engineering polyketide and polypeptide synthases to produce novel products, and tagging combinatorial chemistry libraries with DNA. Importantly, there is still no small-molecule analog of directed protein evolution, i.e. a substantiated approach for optimizing complex (≥ 10^9 diversity) populations of synthetic small molecules over successive generations. We present a key advance towards this goal: a tool for genetically-programmed synthesis of small-molecule libraries from large chemical alphabets. The approach accommodates alphabets that are one to two orders of magnitude larger than any in Nature, and facilitates evolution within the chemical spaces they create. This is critical for small molecules, which are built up from numerous and highly varied chemical fragments. We report a proof-of-concept chemical evolution experiment utilizing an outsized genetic code, and demonstrate that fitness traits can be passed from an initial small-molecule population through to the great-grandchildren of that population. The results establish the practical feasibility of engineering synthetic small molecules through accelerated evolution.

Absolute Intramolecular Distance Measurements with Angstrom-Resolution Using Anomalous Small-Angle X-ray Scattering.

Accurate determination of molecular distances is fundamental to understanding the structure, dynamics, and conformational ensembles of biological macromolecules. Here we present a method to determine the full distance distribution between small (∼7 Å radius) gold labels attached to macromolecules with very high-precision (≤1 Å) and on an absolute distance scale. Our method uses anomalous small-angle X-ray scattering close to a gold absorption edge to separate the gold-gold interference pattern from other scattering contributions. Results for 10-30 bp DNA constructs achieve excellent signal-to-noise and are in good agreement with previous results obtained by single-energy SAXS measurements without requiring the preparation and measurement of single labeled and unlabeled samples. The use of small gold labels in combination with ASAXS read out provides an attractive approach to determining molecular distance distributions that will be applicable to a broad range of macromolecular systems.

The solution structural ensembles of RNA kink-turn motifs and their protein complexes.

With the growing number of crystal structures of RNA and RNA-protein complexes, a critical next step is understanding the dynamic solution behavior of these entities in terms of conformational ensembles and energy landscapes. To this end, we have used X-ray scattering interferometry (XSI) to probe the ubiquitous RNA kink-turn motif and its complexes with the canonical kink-turn binding protein L7Ae. XSI revealed that the folded kink-turn is best described as a restricted conformational ensemble. The ions present in solution alter the nature of this ensemble, and protein binding can perturb the kink-turn ensemble without collapsing it to a unique state. This study demonstrates how XSI can reveal structural and ensemble properties of RNAs and RNA-protein complexes and uncovers the behavior of an important RNA-protein motif. This type of information will be necessary to understand, predict and engineer the behavior and function of RNAs and their protein complexes.

Quantifying Nucleic Acid Ensembles with X-ray Scattering Interferometry.

The conformational ensemble of a macromolecule is the complete description of the macromolecule's solution structures and can reveal important aspects of macromolecular folding, recognition, and function. However, most experimental approaches determine an average or predominant structure, or follow transitions between states that each can only be described by an average structure. Ensembles have been extremely difficult to experimentally characterize. We present the unique advantages and capabilities of a new biophysical technique, X-ray scattering interferometry (XSI), for probing and quantifying structural ensembles. XSI measures the interference of scattered waves from two heavy metal probes attached site specifically to a macromolecule. A Fourier transform of the interference pattern gives the fractional abundance of different probe separations directly representing the multiple conformation states populated by the macromolecule. These probe-probe distance distributions can then be used to define the structural ensemble of the macromolecule. XSI provides accurate, calibrated distance in a model-independent fashion with angstrom scale sensitivity in distances. XSI data can be compared in a straightforward manner to atomic coordinates determined experimentally or predicted by molecular dynamics simulations. We describe the conceptual framework for XSI and provide a detailed protocol for carrying out an XSI experiment.

From a structural average to the conformational ensemble of a DNA bulge.

Direct experimental measurements of conformational ensembles are critical for understanding macromolecular function, but traditional biophysical methods do not directly report the solution ensemble of a macromolecule. Small-angle X-ray scattering interferometry has the potential to overcome this limitation by providing the instantaneous distance distribution between pairs of gold-nanocrystal probes conjugated to a macromolecule in solution. Our X-ray interferometry experiments reveal an increasing bend angle of DNA duplexes with bulges of one, three, and five adenosine residues, consistent with previous FRET measurements, and further reveal an increasingly broad conformational ensemble with increasing bulge length. The distance distributions for the AAA bulge duplex (3A-DNA) with six different Au-Au pairs provide strong evidence against a simple elastic model in which fluctuations occur about a single conformational state. Instead, the measured distance distributions suggest a 3A-DNA ensemble with multiple conformational states predominantly across a region of conformational space with bend angles between 24 and 85 degrees and characteristic bend directions and helical twists and displacements. Additional X-ray interferometry experiments revealed perturbations to the ensemble from changes in ionic conditions and the bulge sequence, effects that can be understood in terms of electrostatic and stacking contributions to the ensemble and that demonstrate the sensitivity of X-ray interferometry. Combining X-ray interferometry ensemble data with molecular dynamics simulations gave atomic-level models of representative conformational states and of the molecular interactions that may shape the ensemble, and fluorescence measurements with 2-aminopurine-substituted 3A-DNA provided initial tests of these atomistic models. More generally, X-ray interferometry will provide powerful benchmarks for testing and developing computational methods.

Structural ensemble and microscopic elasticity of freely diffusing DNA by direct measurement of fluctuations.

Precisely measuring the ensemble of conformers that a macromolecule populates in solution is highly challenging. Thus, it has been difficult to confirm or falsify the predictions of nanometer-scale dynamical modeling. Here, we apply an X-ray interferometry technique to probe the solution structure and fluctuations of B-form DNA on a length scale comparable to a protein-binding site. We determine an extensive set of intrahelix distance distributions between pairs of probes placed at distinct points on the surface of the DNA duplex. The distributions of measured distances reveal the nature and extent of the thermally driven mechanical deformations of the helix. We describe these deformations in terms of elastic constants, as is common for DNA and other polymers. The average solution structure and microscopic elasticity measured by X-ray interferometry are in striking agreement with values derived from DNA-protein crystal structures and measured by force spectroscopy, with one exception. The observed microscopic torsional rigidity of DNA is much lower than is measured by single-molecule twisting experiments, suggesting that torsional rigidity increases when DNA is stretched. Looking forward, molecular-level interferometry can provide a general tool for characterizing solution-phase structural ensembles.

Highly parallel translation of DNA sequences into small molecules.

A large body of in vitro evolution work establishes the utility of biopolymer libraries comprising 10(10) to 10(15) distinct molecules for the discovery of nanomolar-affinity ligands to proteins. Small-molecule libraries of comparable complexity will likely provide nanomolar-affinity small-molecule ligands. Unlike biopolymers, small molecules can offer the advantages of cell permeability, low immunogenicity, metabolic stability, rapid diffusion and inexpensive mass production. It is thought that such desirable in vivo behavior is correlated with the physical properties of small molecules, specifically a limited number of hydrogen bond donors and acceptors, a defined range of hydrophobicity, and most importantly, molecular weights less than 500 Daltons. Creating a collection of 10(10) to 10(15) small molecules that meet these criteria requires the use of hundreds to thousands of diversity elements per step in a combinatorial synthesis of three to five steps. With this goal in mind, we have reported a set of mesofluidic devices that enable DNA-programmed combinatorial chemistry in a highly parallel 384-well plate format. Here, we demonstrate that these devices can translate DNA genes encoding 384 diversity elements per coding position into corresponding small-molecule gene products. This robust and efficient procedure yields small molecule-DNA conjugates suitable for in vitro evolution experiments.

Mesofluidic devices for DNA-programmed combinatorial chemistry.

Hybrid combinatorial chemistry strategies that use DNA as an information-carrying medium are proving to be powerful tools for molecular discovery. In order to extend these efforts, we present a highly parallel format for DNA-programmed chemical library synthesis. The new format uses a standard microwell plate footprint and is compatible with commercially available automation technology. It can accommodate a wide variety of combinatorial synthetic schemes with up to 384 different building blocks per chemical step. We demonstrate that fluidic routing of DNA populations in the highly parallel format occurs with excellent specificity, and that chemistry on DNA arrayed into 384 well plates proceeds robustly, two requirements for the high-fidelity translation and efficient in vitro evolution of small molecules.

Structural and kinetic mapping of side-chain exposure onto the protein energy landscape.

Identification and characterization of structural fluctuations that occur under native conditions is crucial for understanding protein folding and function, but such fluctuations are often rare and transient, making them difficult to study. Native-state hydrogen exchange (NSHX) has been a powerful tool for identifying such rarely populated conformations, but it generally reveals no information about the placement of these species along the folding reaction coordinate or the barriers separating them from the folded state and provides little insight into side-chain packing. To complement such studies, we have performed native-state alkyl-proton exchange, a method analogous to NSHX that monitors cysteine modification rather than backbone amide exchange, to examine the folding landscape of Escherichia coli ribonuclease H, a protein well characterized by hydrogen exchange. We have chosen experimental conditions such that the rate-limiting barrier acts as a kinetic partition: residues that become exposed only upon crossing the unfolding barrier are modified in the EX1 regime (alkylation rates report on the rate of unfolding), while those exposed on the native side of the barrier are modified predominantly in the EX2 regime (alkylation rates report on equilibrium populations). This kinetic partitioning allows for identification and placement of partially unfolded forms along the reaction coordinate. Using this approach we detect previously unidentified, rarely populated conformations residing on the native side of the barrier and identify side chains that are modified only upon crossing the unfolding barrier. Thus, in a single experiment under native conditions, both sides of the rate-limiting barrier are investigated.

Expedient synthesis of a modular phosphate affinity reagent.

Isolation and identification of phosphorylated macromolecules is essential for the deconvolution of most biological regulatory networks. Koike and co-workers recently reported the application of a dinuclear zinc-(pyridylmethyl)amine complex to phosphate-specific affinity purifications and gave it the shorthand name "phos-tag". This complex is valuable for studying phosphorylation because it binds selectively to phosphate dianion in the presence of acidic functional groups at physiological pH, and because the binding is largely independent of molecular context. These properties of phos-tag recommend it for applications in phosphoproteomics, metabolomics, and nucleic acid biology. The catch has been that the molecule is difficult to make and prohibitively expensive to buy. Here, we describe an efficient and inexpensive synthesis of a phos-tag derivative with a versatile alkyne handle. The alkyne handle allows for attachment of phos-tag to alkyl azides via the copper(I)-catalyzed azide-alkyne cycloaddition reaction ("click chemistry"). We characterize the phosphate binding behavior of the new phos-tag derivative in a variety of experimental assays, including its conjugation to a fluorescent reporter, to acrylamide gels, and to sepharose chromatography resin. The synthesis we report should enable a broader use of phos-tag for phosphate-related biochemistry, as both an analytical and a preparative reagent.

Remeasuring the double helix.

DNA is thought to behave as a stiff elastic rod with respect to the ubiquitous mechanical deformations inherent to its biology. To test this model at short DNA lengths, we measured the mean and variance of end-to-end length for a series of DNA double helices in solution, using small-angle x-ray scattering interference between gold nanocrystal labels. In the absence of applied tension, DNA is at least one order of magnitude softer than measured by single-molecule stretching experiments. Further, the data rule out the conventional elastic rod model. The variance in end-to-end length follows a quadratic dependence on the number of base pairs rather than the expected linear dependence, indicating that DNA stretching is cooperative over more than two turns of the DNA double helix. Our observations support the idea of long-range allosteric communication through DNA structure.