Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “homology modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

Improved Biofuel Production through Discovery and Engineering of Terpene Metabolism in Switchgrass

Project Objectives - Of the myriad specialized metabolites that plants deploy to adapt to environmental challenges, terpenes form the largest group. In many major crops, unique terpene blends serve as key stress defenses that directly impact plant fitness and yield. In addition, terpenes, such as bisabolene and pinene, are used for producing renewable biofuels. Essential to advancing a broader use of terpenes for biofuel feedstock engineering is a system-wide knowledge of the diverse biosynthetic machinery and defensive potential of often species-specific terpene blends. The proposed project would merge genome-wide enzyme discovery with comparative –omics, protein structural and plant microbiome studies to define the biosynthesis and stress-defensive functions of the switchgrass (Panicum virgatum) terpene network. These insights would be combined with developing and applying non-transgenic genome editing tools to design plants with desirable terpene blends for higher productivity and biofuel production on marginal lands. As a dedicated lignocellulosic feedstock for U.S. biofuel production with high net energy yield, stress tolerance, and available genome resources, switchgrass is well-suited for devising new avenues for biofuel production. Project Description – The diversity of plant terpene defenses is governed by species-specific families of terpene synthase (TPS) and cytochrome P450 monooxygenase (P450) enzymes. Mining of the switchgrass genome (genotype Alamo) identified ~100 TPS and P450 candidate genes, and combinatorial biochemical analysis of synthesized TPSs and P450s revealed more than a dozen enzymes with common and novel activities. In addition, several identified terpene metabolites and the corresponding transcripts were up-regulated in response to abiotic stressors. These findings demonstrate a unique switchgrass terpene network with probable importance to abiotic stress tolerance, thus providing a large chemical portfolio for optimizing crop resistance, yield, and biofuel composition. Leveraging these preliminary data, we propose to generate a genome-wide map of the switchgrass terpene metabolic network through multi-gene co-expression analyses that allow the efficient cross-validation of TPS and P450 functions. Key enzymes would further be applied to structure-function studies via X-ray protein crystallography, homology modeling and site-directed mutagenesis to gain mechanistic insight into the catalytic specificity of switchgrass terpene metabolism and provide gene and amino acid targets for genome editing. In tandem with terpene pathway discovery, system-wide metabolomics, transcriptomics and proteomics studies in switchgrass accessions of contrasting drought tolerance would define the role of switchgrass terpene metabolism in conferring abiotic stress resilience. Metabolic changes would be assessed in a combined approach of targeted (terpenes) and untargeted metabolite profiling using a high-resolution LC-MS/MS approach, differential gene expression analyses through multiplexed Illumina RNA sequencing, and quantitative analysis of high-priority pathway enzymes using multiple reaction monitoring (MRM). Drawing on these insights, knock-down/out mutants of stress-associated pathway nodes would be generated by optimizing transient virus-induced gene silencing (VIGS) and CRISPR/Cas9 systems under control of the Tobacco Rattle Virus (TRV). The resulting mutant lines would then be analyzed for stress susceptibility and the impact on the root microbiome to define gene functions in planta. Knowledge of terpene pathways, enzyme mechanisms and bioactivities would be applied to enhance switchgrass stress resilience and to tailor-make terpene blends for biofuel production. Here, TRV-enabled CRISPR/Cas9 genome editing, including allele-specific knock-out of redundant genes, engineering of enzyme specificity via structure-guided point mutations, and overexpression of terpene genes relevant to stress-protection or biofuel production, would be used to increase metabolic flux toward desired pathways. Broader Impacts - Integrating the system-wide discovery, mechanistic analysis and non-transgenic genome engineering of the switchgrass terpene network aligns the required steps to unlock the chemical potential of this important metabolite class to generate crops that are more resistant to stress and provide advanced biofuel production in light of rising climate pressures as foreseeable challenges for bioenergy crop cultivation. The proposed project would further offer interdisciplinary student training through active involvement in the project and integration of research concepts and outcomes into newly-developed graduate and undergraduate courses on Plant Biotechnology.

09 BIOMASS FUELS↗

Cytochrome Nanowire XRD

Microbial metabolism plays an important role in the global cycling of carbon, nutrients, and metals. Recent discovery of microorganisms sharing energy using direct electrical connections might represent a prominent strategy by which anaerobic microorganisms interact with the diversity of environments. The purpose of the proposed research is to apply these newly developed principles of electron exchange in microbial communities to the study of biogeochemical cycling in climate and subsurface systems. Previous studies revealed that soil bacteria Geobacter sulfurreducens produce conductive pili nanofilaments that facilitate long-range electron transport to extracellular electron acceptors and other cells. This discovery has already transformed our understanding of the function of microbial communities in diverse environments, making basic contributions to microbial ecology and helping to improve practical applications such as increasing current output of microbial fuel cells and enhancing the conversion of organic waste to methane. Microbial electron exchange via pili may also contribute to methane production in terrestrial environments that are major sources of atmospheric methane. However, the mechanism of electron transport through pili is not well understood. Initial studies have suggested that pili exhibit metallic-like conductivity similar to conducting polymer polyaniline. However, these studies were limited to electrical measurements. X-ray diffraction studies have suggested that aromatic amino acids tyrosine and phenylalanine are closely packed in pili. But these studies used homology modeling to predict the pilus assembly and information about the number of monomers in each pilus assembly unit and their organization are still lacking. The central hypothesis of this proposal is aromatics in pili enable intermolecular electron delocalization due to pi-stacking that can give rise to metallic-like conductivity. By analyzing pili structure using infrared scattering-scanning near field optical microscope (IR s-SNOM) and combining with solution NMR, Focused Ion Beam (FIB and computational capabilities, this work will provide comprehensive understanding of the conduction mechanism from a structure-function perspective. The specific aims are: 1) With infrared scattering scanning near field optical microscope (IR s-SNOM), image the organization and density of pilin monomers and pH-induced conformational changes in a pilus filament. 2) With solution NMR, resolve pH-induced structural changes in the pilin monomer. 3) Using Focused Ion Beam (FIB), draw electrodes on single pili to measure their conductivity. 4) With computational modeling using the NWChem, simulate the effect of environment on the electronic properties of pili as a function of pH. The proposed work will gain the nanoscale insight into the structure of the pili, as the basis for understanding interspecies electrical communications. The improved understanding will be harnessed to accelerate bioenergy production and bioremediation using pili as well as for predictive modeling of carbon cycling in diverse environments to reach the central goals of DOE-BER to develop sustainable energy sources, regulate the contaminants in the subsurface and understand the effects of greenhouse gas emissions on biosphere. The proposed studies will provide design principles for engineering of electrical interactions among these communities by manipulating the amino acid composition of pili. Moreover, these studies will identify the forces that maintain the function of these communities in the face of environmental perturbation. Therefore, these studies will help to address the need "to incorporate process understanding of biogeochemical cycling into a scalable hierarchy of predictive capabilities" specified in EMSL's Science Themes.

59 BASIC BIOLOGICAL SCIENCES↗

Mechanosensitive Ion Channels in Bacteria: Functional Domains and Mechanisms of Gating

The past funding period was productive for the group. The progress in the mechanosensitive channel field was critically affected in the end of 1998 by the solution of the crystal structure of the mycobacterial homolog of MscL by our colleagues from Caltech. Having the structure of TbMscL in the closed state, we developed a detailed homology model of EcoMscL, and related the structural model with the wealth of functional phenomenology available for the E. coli version of the channel (EcoMscL). The biophysical properties of the open MscL helped to model the open conformation and infer the pathway for the entire gating transition. The following experiments provided strong support to the atomic model of the gating process, and allowed to make further predictions. The work has advanced our understanding of tension-driven conformational transitions in membrane-embedded mechanosensory proteins, determine major energetic contributions and set the stage for further exploration of the whole family of mechanosensitive channels. The results have been published in seven experimental and theoretical papers, with three other papers currently in press or in preparation.

Sukharev, Sergei↗

Calcium-stimulated autophosphorylation site of plant chimeric calcium/calmodulin-dependent protein kinase

The existence of two molecular switches regulating plant chimeric Ca(2+)/calmodulin-dependent protein kinase (CCaMK), namely the C-terminal visinin-like domain acting as Ca(2+)-sensitive molecular switch and calmodulin binding domain acting as Ca(2+)-stimulated autophosphorylation-sensitive molecular switch, has been described (Sathyanarayanan, P. V., Cremo, C. R., and Poovaiah, B. W. (2000) J. Biol. Chem. 275, 30417-30422). Here we report the identification of Ca(2+)-stimulated autophosphorylation site of CCaMK by matrix-assisted laser desorption ionization time of flight-mass spectrometry. Thr(267) was confirmed as the Ca(2+)-stimulated autophosphorylation site by post-source decay experiments and by site-directed mutagenesis. The purified T267A mutant form of CCaMK did not show Ca(2+)-stimulated autophosphorylation, autophosphorylation-dependent variable calmodulin affinity, or Ca(2+)/calmodulin stimulation of kinase activity. Sequence comparison of CCaMK from monocotyledonous plant (lily) and dicotyledonous plant (tobacco) suggests that the autophosphorylation site is conserved. This is the first identification of a phosphorylation site specifically responding to activation by second messenger system (Ca(2+) messenger system) in plants. Homology modeling of the kinase and calmodulin binding domain of CCaMK with the crystal structure of calcium/calmodulin-dependent protein kinase 1 suggests that the Ca(2+)-stimulated autophosphorylation site is located on the surface of the kinase and far from the catalytic site. Analysis of Ca(2+)-stimulated autophosphorylation with increasing concentration of CCaMK indicates the possibility that the Ca(2+)-stimulated phosphorylation occurs by an intermolecular mechanism.

Non-NASA Center↗

Modeling Non-homologous End Joining

Non-homologous end joining (NHEJ) is the dominant DNA double strand break (DSB) repair pathway and involves several NHEJ proteins such as Ku, DNA-PKcs, XRCC4, Ligase IV and so on. Once DSBs are generated, Ku is first recruited to the DNA end, followed by other NHEJ proteins for DNA end processing and ligation. Because of the direct ligation of break ends without the need for a homologous template, NHEJ turns out to be an error-prone but efficient repair pathway. Some mechanisms have been proposed of how the efficiency of NHEJ repair is affected. The type of DNA damage is an important factor of NHEJ repair. For instance, the length of DNA fragment may determine the recruitment efficiency of NHEJ protein such as Ku [1], or the complexity of the DNA breaks [2] is accounted for the choice of NHEJ proteins and subpathway of NHEJ repair. On the other hand, the chromatin structure also plays a role of the accessibility of NHEJ protein to the DNA damage site. In this talk, some mathematical models of NHEJ, that consist of series of biochemical reactions complying with the laws of chemical reaction (e.g. mass action, etc.), will be introduced. By mathematical and numerical analysis and parameter estimation, the models are able to capture the qualitative biological features and show good agreement with experimental data. As conclusions, from the viewpoint of modeling, how the NHEJ proteins are recruited will be first discussed for connection between the classical sequential model [4] and recently proposed two-phase model [5]. Then how the NHEJ repair pathway is affected, by the length of DNA fragment [6], the complexity of DNA damage [7] and the chromatin structure [8], will be addressed

Li, Yongfeng↗

Three-Dimensional Modeling of Quasi-Homologous Solar Jets

Recent solar observations (e.g., obtained with Hinode and STEREO) have revealed that coronal jets are a more frequent phenomenon than previously believed. This higher frequency results, in part, from the fact that jets exhibit a homologous behavior: successive jets recur at the same location with similar morphological features. We present the results of three-dimensional (31)) numerical simulations of our model for coronal jets. This study demonstrates the ability of the model to generate recurrent 3D untwisting quasi-homologous jets when a stress is constantly applied at the photospheric boundary. The homology results from the property of the 3D null-point system to relax to a state topologically similar to its initial configuration. In addition, we find two distinct regimes of reconnection in the simulations: an impulsive 3D mode involving a helical rotating current sheet that generates the jet, and a quasi-steady mode that occurs in a 2D-like current sheet located along the fan between the sheared spines. We argue that these different regimes can explain the observed link between jets and plumes.

Pariat, E.↗

Cannabis monoterpene synthases: evaluating structure–function relationships

Terpene synthases catalyze the first committed step in the biosynthesis of terpenes, a structurally diverse class of natural products that also encompasses volatiles derived from precursors in the C10 to C15 range (termed monoterpenes and sesquiterpenes, respectively). In the review section of this article, we are providing information about all functionally characterized monoterpene synthases (MTSs) and sesquiterpene synthases (STSs) of Cannabis sativa L. We are also exploring the locations of MTSs and STSs in the chromosome-level assembly of the reference chemovar CBDRx. A follow-up computational structure–function analysis focuses on MTSs, as there is already a rich literature available on the topic. More specifically, by employing sequence comparisons and homology structural modeling, we infer which amino acid residues are likely to constrain the available space in the active site of cannabis MTSs. The emphasis of these studies was to investigate why some MTSs accept only a C10 diphosphate as substrate, while mixed MTS/STS enzymes also accommodate a C15 diphosphate. Here, by combining a literature review and computational analyses in a hybrid format, we are laying the foundation for future studies to better understand the determinants of substrate and product specificity in these fascinating enzymes.

59 BASIC BIOLOGICAL SCIENCES↗

Osteocalcin binds to a GPRC6A Venus fly trap allosteric site to positively modulate GPRC6A signaling

GPRC6A, a member of the Family C G-protein coupled receptors, regulates energy metabolism and sex hormone production and is activated by diverse ligands, including cations, L-amino acids, the osteocalcin (Ocn) peptide and the steroid hormone testosterone. We sought a structural framework for the ability of multiple distinct classes of ligands to active GPRC6A. We created a structural model of GPRC6A using Alphafold2. Using this model we explored a putative orthosteric ligand binding site in the bilobed Venus fly trap (VFT) domain of GPRC6A and two positive allosteric modulator (PAM) sites, one in the VFT and the other in the 7 transmembrane (7TM) domain. We provide evidence that Ocn peptides act as a PAM for GPRC6A by binding to a site in the VFT that is distinct from the orthosteric site for calcium and L-amino acids. In agreement with this prediction, alternatively spliced GPRC6A isoforms 2 and 3, which lack regions of the VFT, and mutations in the computationally predicted Ocn binding site, K352E and H355P, prevent Ocn activation of GPRC6A. These observations explain how dissimilar ligands activate GPRC6A and set the stage to develop novel molecules to activate and inhibit this previously poorly understood receptor.

59 BASIC BIOLOGICAL SCIENCES↗

Expanded Very Large Array Nova Project Observations of the Classical NovaV1723 Aquilae

We present radio light curves and spectra of the classical nova VI723 Aql obtained with the Expanded Very Large Array (EVLA). This is the first paper to showcase results from the EVLA Nova Project, which comprises a team of observers and theorists utilizing the greatly enhanced sensitivity and frequency coverage of EVLA radio observations, along with observations at other wavelengths, to reach a deeper understanding of the energetics, morphology, and temporal characteristics of nova explosions. Our observations of VI723 Aql span 1-37 GHz in frequency, and we report on data from 14 to 175 days following the time of the nova explosion. The broad frequency coverage and frequent monitoring show that the radio behavior of VI723 Aql does not follow the classic Hubble-flow model of homologous spherically expanding thermal ejecta. The spectra are always at least partially optically thin, and the flux rises on faster timescales than can be reproduced with linear expansion. Therefore, any description of the underlying physical processes must go beyond this simple picture. The unusual spectral properties and light curve evolution might be explained by multiple emitting regions or shocked material. Indeed, X-ray observations from Swift reveal that shocks are likely present.

Krauss, Miriam I.↗

A hybrid surrogate modeling framework for the Digital Twin of a Fluoride-salt-cooled High-temperature Reactor (FHR)

While nuclear energy is a non-greenhouse-gas emitting energy source, expensive operational costs due to the high-level of safety requirements decreases their competitiveness in the sustainable energy market. Advanced reactor concepts paired with Digital Twins aim to increase the commercialization gains of nuclear energy by reducing operational costs, increasing reactor reliability and enhancing power generation. To support Digital Twin tasks such as real-time autonomous control, proactive maintenance monitoring or optimizing power demand operations, a fast and accurate virtual representation of the Nuclear Power Plant (NPP) is required. The computational cost of high-fidelity, physics-based models are unsuitable for real-time analysis or scalability. Here, in this work, a hybrid surrogate modeling framework is developed fora Fluoride-salt-cooled High-temperature Reactor (FHR) that leverages physics-inspired models for key reactor components and uses data-driven methods for rapid system state space prediction. The Xenon reactivity feedback model is integrated to inform the surrogate model about the reactor core and the homologous pump theory model is the basis for representing pump degradation. Using a detailed, two dimensional thermal hydraulics model to generate data on the FHR, we train a network of Vectorized Autoregressive Moving-Average with eXogenous input (VARMAX) models to predict the remaining state values. The result is a surrogate model that provides a detailed reactor state representation of 41 system states and a pump degradation analysis. The framework is applied to Load Follows profiles, yielding high accuracy and a speedup that is more than 4000x faster compared to the higher- fidelity thermal hydraulics model, enabling real-time operational intelligence and applications in long horizon predictions. While the surrogate model framework is demonstrated for the particular case of FHR, the hybrid physical/data-driven modeling approach including the network of surrogates and the underlying modularity has the potential to be applied to other physical asset systems.

Digital Twins↗

A comparative modeling of supernova 1993J

The light curve of Supernova 1993J is calculated using two computational radiation transport approaches. The two approaches are represented by the computer codes STELLA and EDDINGTON. The emphasis is on the shock breakout and the photometry in the U, B and V bands during the first 120 days of the supernova. The STELLA model includes implicit hydrodynamics and is able to model early supernova evolution before the expansion is homologous. The STELLA model employs multi-group photonics and is able to follow the radiation as it decouples from the matter. The EDDINGTON code uses an algorithm for integrating the transport equation which assumes homologous expansion and uses a finer frequency resolution. The agreement between the two codes is considered to be satisfactory only in the case where compatible physical assumptions are made concerning the opacity. The assumptions are justified. The continuum spectrum for SN 1993J is predicted near the shock breakout to be superior to that predicted by standard single energy group hydrocodes. The uncertainties involved in current time dependent models of supernova light curves are discussed.

Blinnikov, Sergei↗

An FPGA-based hardware accelerator supporting sensitive sequence homology filtering with profile hidden Markov models

Abstract Background Sequence alignment lies at the heart of genome sequence annotation. While the BLAST suite of alignment tools has long held an important role in alignment-based sequence database search, greater sensitivity is achieved through the use of profile hidden Markov models (pHMMs). Here, we describe an FPGA hardware accelerator, called HAVAC, that targets a key bottleneck step (SSV) in the analysis pipeline of the popular pHMM alignment tool, HMMER. Results The HAVAC kernel calculates the SSV matrix at 1739 GCUPS on a $$\sim$$ ∼ $3000 Xilinx Alveo U50 FPGA accelerator card, $$\sim$$ ∼ 227× faster than the optimized SSV implementation in nhmmer . Accounting for PCI-e data transfer data processing, HAVAC is 65× faster than nhmmer’s SSV with one thread and 35× faster than nhmmer with four threads, and uses $$\sim$$ ∼ 31% the energy of a traditional high end Intel CPU. Conclusions HAVAC demonstrates the potential offered by FPGA hardware accelerators to produce dramatic speed gains in sequence annotation and related bioinformatics applications. Because these computations are performed on a co-processor, the host CPU remains free to simultaneously compute other aspects of the analysis pipeline.

59 BASIC BIOLOGICAL SCIENCES↗

Structural determination of a full-length plant cellulose synthase informed by experimental and in silico methods

Three-dimensional structure determination and prediction of proteins with intrinsically disordered regions, unstructured regions, conformational flexibility, and lacking homologous structures are challenging. We previously predicted and refined an in silico structure of a plant cellulose synthase from cotton (GhCESA1), and more recently, cryo-electron microscopy (cryo-EM) has resolved a majority of the lengths of two CESA structures from poplar (PttCESA8) and cotton (GhCESA7). However, 26–30% of these cryo-EM structures remain unresolved, including the N-terminal domain, half of the class-specific region, the gating loop region, and the C-terminal domain. Here, we describe the generation and evaluation of a full-length hybrid GhCESA1 model based on this cryo-EM PttCESA8 structure, with unresolved regions completed using this in silico refined GhCESA1 model. All-atom molecular dynamics simulations and subsequent energy minimizations were performed for the in silico and hybrid GhCESA1 models in a lipid bilayer-water-ion environment, and structural stability, dynamics, energetics, contacts, and quality were evaluated. The unresolved regions were found to be the most dynamic, in agreement with their poor electron density with cryo-EM. The hybrid model exhibited a higher total secondary structure content, more favorable intra-protein and protein-lipid interaction energies, and improved quality metrics. Moreover, hydrogen bonding was revealed to be a primary mechanism for intra-protein and protein-lipid contacts. These results demonstrate that in silico structure prediction and refinement may be useful to augment experimental structure determination, especially for disordered and unstructured regions. Furthermore, this hybrid model can serve as a steppingstone to derive full-length homology models of other CESAs found in more experimentally tractable organisms.

59 BASIC BIOLOGICAL SCIENCES↗

Empirical evidence that glucan-interacting amino acid side chains within the transmembrane channel collectively facilitate cellulose synthase function

The fundamental mechanism of cellulose synthesis is widely conserved across Kingdoms and depends on cellulose synthases, which are processive, dual-function, family 2 glycosyltransferases (GT-2). These enzymes polymerize glucose on the cytoplasmic side of the plasma membrane and export the glucan chain to the cell surface through an integral transmembrane (TM) channel. Structural studies of active plant cellulose synthases (CESAs) have revealed interactions between the nascent glucan chain and the side chains of polar, charged, and aromatic amino acid residues that line the TM channel. However, the functional consequences of modifying these side chains have not been tested in vivo in CESAs or other processive GT-2s. To test this, we used an established in vivo assay based on genetic complementation of CESA5 in the moss, Physcomitrium patens. For accurate prediction of glucan-interacting amino acid residues, we generated a complete homotrimeric molecular model of PpCESA5 using a combination of homology and de novo modeling. All-atom molecular dynamics-based analyses of contact metrics and interaction energy identified 23 amino acid residues with high propensity to interact with the nascent glucan chain within the TM channel or on the apoplastic surface of PpCESA5. Mutating any one of 18 of these amino acid residues to alanine, thereby removing their side chains, abolished or impaired CESA function, with the strongest effects observed upon the loss of charged amino acid side chains. This provides direct evidence to support the hypothesis that multiple amino acid residues collectively maintain a smooth energy landscape within the TM channel to facilitate glucan translocation.

59 BASIC BIOLOGICAL SCIENCES↗

Biochemical Kinetics Model of DSB Repair and GammaH2AX FOCI by Non-homologous End Joining

We developed a biochemical kinetics approach to describe the repair of double strand breaks (DSB) produced by low LET radiation by modeling molecular events associated with the mechanisms of non-homologous end-joining (NHEJ). A system of coupled non-linear ordinary differential equations describes the induction of DSB and activation pathways for major NHEJ components including Ku(sub 70/80), DNA-PK(sub cs), and the Ligase IV-XRCC4 hetero-dimer. The autophosphorylation of DNA-PK(sub cs and subsequent induction of gamma-H2AX foci observed after ionizing radiation exposure were modeled. A two-step model of DNA-PK(sub cs) regulation of repair was developed with the initial step allowing access of other NHEJ components to breaks, and a second step limiting access to Ligase IV-XRCC4. Our model assumes that the transition from the first to second-step depends on DSB complexity, with a much slower-rate for complex DSB. The model faithfully reproduced several experimental data sets, including DSB rejoining as measured by pulsed-field electrophoresis (PFGE), quantification of the induction of gamma-H2AX foci, and live cell imaging of the induction of Ku(sub 70/80). Predictions are made for the behaviors of NHEJ components at low doses and dose-rates, where a steady-state is found at dose-rates of 0.1 Gy/hr or lower.

Cucinotta, Francis, A.↗

No ordinary proteins: Adsorption and molecular orientation of monoclonal antibodies

The interaction of monoclonal antibodies (mAbs) with air/water interfaces plays a crucial role in their overall stability in solution. We aim to understand this behavior using pendant bubble measurements to track the dynamic tension reduction and x-ray reflectivity to obtain the electron density profiles (EDPs) at the surface. Native immunoglobulin G mAb is a rigid molecule with a flat, “Y” shape, and simulated EDPs are obtained by rotating a homology construct at the surface. Comparing simulations with experimental EDPs, we obtain surface orientation probability maps showing mAbs transition from flat-on Y-shape configurations to side-on or end-on configurations with increasing concentration. The modeling also shows the presence of β sheets at the surface. Overall, the experiments and the homology modeling elucidate the orientational phase space during different stages of adsorption of mAbs at the air/water interface. These finding will help define new strategies for the manufacture and storage of antibody-based therapeutics.

59 BASIC BIOLOGICAL SCIENCES↗

A Wox3 -patterning module organizes planar growth in grass leaves and ligules

Grass leaves develop from a ring of primordial initial cells within the periphery of the shoot apical meristem, a pool of organogenic stem cells that generates all of the organs of the plant shoot. At maturity, the grass leaf is a flattened, strap-like organ comprising a proximal supportive sheath surrounding the stem and a distal photosynthetic blade. The sheath and blade are partitioned by a hinge-like auricle and the ligule, a fringe of epidermally derived tissue that grows from the adaxial (top) leaf surface. Together, the ligule and auricle comprise morphological novelties that are specific to grass leaves. Understanding how the planar outgrowth of grass leaves and their adjoining ligules is genetically controlled can yield insight into their evolutionary origins. Here we use single-cell RNA-sequencing analyses to identify a ‘rim’ cell type present at the margins of maize leaf primordia. Cells in the leaf rim have a distinctive identity and share transcriptional signatures with proliferating ligule cells, suggesting that a shared developmental genetic programme patterns both leaves and ligules. Moreover, we show that rim function is regulated by genetically redundant Wuschel-like homeobox3 (WOX3) transcription factors. Higher-order mutations in maize Wox3 genes greatly reduce leaf width and disrupt ligule outgrowth and patterning. Together, these findings illustrate the generalizable use of a rim domain during planar growth of maize leaves and ligules, and suggest a parsimonious model for the homology of the grass ligule as a distal extension of the leaf sheath margin.

59 BASIC BIOLOGICAL SCIENCES↗