SEARCH · Search NASA
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.
There exists a broad class of sequencing problems in soft materials such as proteins and polymers that can be formulated as a heuristic search that involves decision making akin to a computer game.
SUMMARY Populus species play a foundational role in diverse ecosystems and are important renewable feedstocks for bioenergy and bioproducts. Hybrid aspen Populus tremula × P. alba INRA 717‐1B4 is a widely used transformation model in tree functional genomics and biotechnology research. As an outcrossing interspecific hybrid, its genome is riddled with sequence polymorphisms which present a challenge for sequence‐sensitive analyses. Here we report a telomere‐to‐telomere genome for this hybrid aspen with two chromosome‐scale, haplotype‐resolved assemblies. We performed a comprehensive analysis of the repetitive landscape and identified both tandem repeat array‐based and array‐less centromeres. Unexpectedly, the most abundant satellite repeats in both haplotypes lie outside of the centromeres, consist of a 147 bp monomer PtaM147, frequently span >1 megabases, and form heterochromatic knobs. PtaM147 repeats are detected exclusively in aspens (section Populus ) but PtaM147‐like sequences occur in LTR‐retrotransposons of closely related species, suggesting their origin from the retrotransposons. The genomic resource generated for this transformation model genotype has greatly improved the design and analysis of genome editing experiments that are highly sensitive to sequence polymorphisms. The work should motivate future hypothesis‐driven research to probe into the function of the abundant and aspen‐specific PtaM147 satellite DNA.
Neurodegenerative diseases are characterized by the pathologic accumulation of aggregated proteins. Known as amyloid, these fibrillar aggregates include proteins such as tau and amyloid-β (Aβ) in Alzheimer’s disease (AD) and alpha-synuclein (αSyn) in Parkinson’s disease (PD). The development and spread of amyloid fibrils within the brain correlates with disease onset and progression, and inhibiting amyloid formation is a possible route toward therapeutic development. Recent advances have enabled the determination of amyloid fibril structures to atomic-level resolution, improving the possibility of structure-based inhibitor design. In this work, we use these amyloid structures to design inhibitors that bind to the ends of fibrils, “capping” them so as to prevent further growth. Using de novo protein design, we develop a library of miniprotein inhibitors of 35 to 48 residues that target the amyloid structures of tau, Aβ, and αSyn. Biophysical characterization of top in silico designed inhibitors shows they form stable folds, have no sequence similarity to naturally occurring proteins, and specifically prevent the aggregation of their targeted amyloid-prone proteins in vitro. The inhibitors also prevent the seeded aggregation and toxicity of fibrils in cells. In vivo evaluation reveals their ability to reduce aggregation and rescue motor deficits in Caenorhabditis elegans models of PD and AD.
Peptide materials offer a broad platform to design biomimetic soft matter, and filamentous networks that emulate those in extracellular matrices and the cytoskeleton are among the important targets. Given the vast sequence space, a combination of computational approaches and readily accessible experimental techniques is required to design peptide materials efficiently. Here, we report here on a strategy that utilizes this combination to predict supramolecular cohesion within filaments of peptide amphiphiles, a property recently linked to supramolecular dynamics and consequently bioactivity. Using established coarse-grained simulations on 10,000 randomly generated peptide sequences, we identified 3500 likely to self-assemble in water into nanoscale filaments. Atomistic simulations of small clusters were used to further analyze this subset of sequences and identify mathematical descriptors that are predictive of intermolecular cohesion, which was the main purpose of this work. We arbitrarily selected a small cohort of these sequences for chemical synthesis and verified their fiber morphology. With further characterization, we were able to link the latent heat associated with fiber to micelle transitions, an indicator of cohesion and potential supramolecular dynamicity within the filaments, to calculated hydrogen bond densities in the simulation clusters. Based on validation from in situ synchrotron X-ray scattering and differential scanning calorimetry, we conclude that the phase transitions can be easily observed by very simple polarized light microscopy experiments. We are encouraged by the methodology explored here as a relatively low-cost and fast way to design potential functions of peptide materials.
Transmitted/founder (T/F) HIV-1 envelope proteins (Envs) from infected individuals that developed neutralization breadth are likely to possess inherent features desirable for vaccine immunogen design. To explore this premise, we conducted an immunization study in rhesus macaques (RM) using T/F Env sequences from two human subjects, one of whom developed potent and broad neutralizing antibodies (Z1800M) while the other developed little to no neutralizing antibody responses (R66M) during HIV-1 infection. Using a DNA/MVA/protein immunization protocol, 10 RM were immunized with each T/F Env. Within each T/F Env group, the protein boosts were administered as either monomeric gp120 or stabilized trimeric gp140 protein. All vaccination regimens elicited high titers of antigen-specific IgG, and two animals that received monomeric Z1800M Env gp120 developed autologous neutralizing activity. Using early Env escape variants isolated from subject Z1800M as guides, the serum neutralizing activity of the two immunized RM was found to be dependent on the gp120 V5 region. Interestingly, the exact same residues of V5 were also targeted by a neutralizing monoclonal antibody (nmAb) isolated from the subject Z1800M early in infection. Glycan profiling and computational modeling of the Z1800M Env gp120 immunogen provided further evidence that the V5 loop is exposed in this T/F Env and was a dominant feature that drove neutralizing antibody targeting during infection and immunization. An expanded B cell clonotype was isolated from one of the neutralization-positive RM and nmAbs corresponding to this group demonstrated V5-dependent neutralization similar to both the RM serum and the human Z1800M nmAb. The results demonstrate that neutralizing antibody responses elicited by the Z1800M T/F Env in RM converged with those in the HIV-1 infected human subject, illustrating the potential of using immunogens based on this or other T/F Envs with well-defined immunogenicity as a starting point to drive breadth.
With the advances in next generation sequencing technologies, the number of sequenced genomes is growing exponentially, resulting in a technology bottleneck for the translation of sequence information into usable hypotheses about the function of each gene. We have proposed leveraging our leadership high-performance computing (HPC) resources to help break this annotation bottleneck. Here we design an HPC-based framework to infer gene function from gene sequence by incorporating information about protein structure and interactions predicted by deep learning approaches. Accurate functional prediction and gene annotation using computational methods will facilitate breakthroughs in the genomic sciences essential to understanding and harnessing life processes in bacteria, fungi and plants. The development and applications of the state-of-the-art deep neural networks to protein structural modeling, interaction prediction, sequence comparison, and quality assessment of protein structural models will be made possible by leadership computational resources. These HPC-enabled bioinformatics and molecular modeling tools will provide powerful insights into molecular functions of genes.
Background: HIV-1 envelope (Env) variable loops 1 and 2 (V1V2) directed non-neutralizing antibodies were a correlate of decreased transmission risk in the RV144 vaccine trial. Thus, the elicitation and breadth of antibody responses against the V1V2 of HIV-1 Env are important considerations for HIV-1 vaccine candidates. The V1V2 region’s highly variable nature and the extensive diversity of subtype C HIV-1 Envelopes (Envs) make the V1V2 response breadth a high priority for HIV-1 vaccine regimens aiming for V1V2-mediated protection in Southern Africa. Here, we determined whether the breadth of the anti-V1V2 vaccine response can be broadened by including HIV-1 Env strains computationally designed to enhance the coverage of subtype C V1V2 sequence diversity. Methods: Three subtype C Env strains were selected to maximize antibody binding coverage while complementing subtype C vaccine gp120s that were given in human clinical trials in South Africa, as well as to improve epitope accessibility. Humoral immunogenicity of a novel trivalent gp120 vaccine immunogen, a bivalent gp120 boost already in clinical trials (1086C and TV1), and a pentavalent (all five gp120s combined) were evaluated in a preclinical immunization study in guinea pigs. The pentavalent combination was further evaluated with alum versus glucopyranosyl lipid adjuvants formulated in squalene-in-water emulsion (GLA-SE) adjuvants in non-human primates. The breadth of the anti-V1V2 response was assessed using an array of cross-subtype variable loops 1&2 (V1V2) scaffold proteins and linear V2 peptides. Results: The breadth of the IgG response against V1V2 antigens of the trivalent and pentavalent groups was comparable, and both were greater than the breadth of the bivalent group. Linear epitope mapping showed that two linear epitopes in V2 were targeted by the vaccinated animals: the V2 hotspot focused at 169K that potentially correlated with decreased HIV-1 risk in RV144 and the V2.2 site (179LDV/I181) that is part of the integrin α4β7 binding site. The bivalent vaccine elicited a significantly higher magnitude of binding to the V2 hotspot compared to the trivalent vaccine whereas the trivalent vaccine elicited significantly higher binding to the V2.2 epitope compared to the bivalent vaccine, while the pentavalent recognized both regions. Conclusions: These results demonstrate that the three new computationally selected subtype C Envs successfully complemented 1086C and TV1 for broader V1V2 antibody responses, and, in concert with adjuvants that stimulate V1V2 responses, can be considered as part of a rationale immunogen design to improve V1V2 IgG coverage in future vaccine trials in South Africa.
Work in the project was spread across 4 main themes:1) Improved Understanding of Why Cracks Form and Cause Degradation in Solar Cells: We developed a finite element model which correlated well with experimental determined deflection vs load profiles and which predicted stress profiles vs load level and temperature. These profiles helped to guide brace designs and showed the high stress levels seen in the silicon near the interconnect wire locations as the temperature is cooled below -30C, likely forming "invisible" microcracks. Experiments showed how for polymer backsheet modules, a brief exposure to such low temperatures made modules extremely sensitive to propagation of these microcracks at relatively low load levels, and how thicker and softer encapsulant helped to reduce this crack sensitivity. We showed for some modules that most cracks that are seen under loading had their origin at short, hard to see, V-cracks from soldering damage, and how the shunting from cells with high densities of cracks reduces energy delivery more than Pmax. 2) Improved Durability Testing Methods With Respect to Cracked Cells: We lobbied the PV community throughout the project that IEC 61215 was insufficient to probe the sensitivity of module Pmax degradation due to cracks in that in did not have a test leg that included a crack creation step (static loading) followed by a crack opening step (cyclic loading and/or thermal cycling). We published extensively on how pre- existing cracks opened up with cyclic loading. The new IEC TS 63209-1 ED1 test sequence now include our recommendations, and this should lead to improved module designs and reduced risk for investors. We also published recommendations for a new "Cold Climate Test" sequence to probe the sensitivity of modules to degradation from short microcracks formed when modules see temperatures below -30C followed by front side loading. We participated in several USTAG meetings related to cracked cells, EL testing, and mechanical load testing and provided input to the group to improve new standards under consideration. We also showed how the act of performing an EL test can affect the test results due to heating from the injected current. 3) New Module Designs and Related Manufacturing Methods & Equipment: We demonstrated reduced crack sensitivity on encapsulated cell coupons the benefits of building protective compressive stress into the cells by laminating them in a bent state and then allowing the coupon to flatten out after lamination. Our attempts to scale this to the module level with a more manufacturable sequence did not demonstrate the same positive effects, but we hope that others will follow up with other variations on this bent lamination or bent cooling concept to make more durable modules. 4) New Module Mounting Designs: We explored a range of bracing concepts which pressed on the backsheet to slightly bow the glass outward and place cells in a state of protective compressive stress and to reduce the deflection under front side loads. These successfully reduced cell cracking and crack opening, but high static loads permanently deformed cost-effective braces, after which they provided little benefit. We settled on an elegant & effective solution – the RailPad, where lightweight bracing elements are placed between the mounting rail(s) and the backsheet, and the existing strength of the rails is used to prevent module deflection under load. New construction and retrofit designs were made, and test sites were implemented in FL and MA. We applied for RailPad patent protection.
Bacterial mRNA sequencing is inefficient due to the abundance of ribosomal RNA that is challenging to deplete. While commercial kits target rRNA from common bacterial species, they are frequently inefficient when applied to divergent species, including those from environmental isolates. Similarly, other methods typically employ large probe sets that tile the entire length of rRNAs; however, such approaches are infeasible when applied to many species. Therefore, we present EMBR-seq+, which requires fewer than 10 oligonucleotides per rRNA by combining rRNA blocking primers with RNase H-mediated depletion to achieve rRNA removal efficiencies of up to 99% in diverse bacterial species. Furthermore, in more complex microbial co-cultures between Fibrobacter succinogenes strain UWB7 and anaerobic fungi, EMBR-seq+ depleted both bacterial and fungal rRNA, with a fourfold improvement in bacterial rRNA depletion compared with a commercial kit, thereby demonstrating that the method can be applied to non-model microbial mixtures. Notably, for microbes with unknown rRNA sequences, EMBR-seq+ enables rapid iterations in probe design without requiring to start experiments from total RNA. Finally, efficient depletion of rRNA enabled systematic quantification of the reprogramming of the bacterial transcriptome when cultured in the presence of the anaerobic fungi Anaeromyces robustus or Caecomyces churrovis. We observed that F. succinogenes strain UWB7 downregulated several lignocellulose-degrading carbohydrate-active enzymes in the presence of anaerobic gut fungi, suggesting close interactions between two cellulolytic species that specialize in different aspects of biomass breakdown. Thus, EMBR-seq+ enables efficient, cost-effective, and rapid quantification of the transcriptome to gain insights into non-model microbial systems.
Caffeine (CFF) is a common antagonist to the four subtypes of adenosine G-protein-coupled receptors (GPCRs), which are critical drug targets for treating heart failure, cancer, and neurological diseases. However, the pathways and mechanism of CFF binding to the target receptors remain unclear. In this study, we have performed all-atom-enhanced sampling simulations using a robust Gaussian-accelerated molecular dynamics (GaMD) method to elucidate the binding mechanism of CFF to human adenosine A 2A receptor (A 2A AR). Multiple 500–1,000 ns GaMD simulations captured both binding and dissociation of CFF in the A 2A AR. The GaMD-predicted binding poses of CFF were highly consistent with the x-ray crystal conformations with a characteristic hydrogen bond formed between CFF and residue N6.55 in the receptor. In addition, a low-energy intermediate binding conformation was revealed for CFF at the receptor extracellular mouth between ECL2 and TM1. While the ligand-binding pathways of the A 2A AR were found similar to those of other class A GPCRs identified from previous studies, the ECL2 with high sequence divergence serves as an attractive target site for designing allosteric modulators as selective drugs of the A 2A AR.
Molecular modeling of a supramolecular catalytic system is conducted resulting from the assembling between a small peptide and the surface of cationic self-assembled monolayers on gold nanoparticles, through a multiscale iterative approach including atomistic force field development, flexible docking with Brownian Dynamics and µs-long Molecular Dynamics simulations. Self-assembly is a prerequisite for the catalysis, since the catalytic peptides do not display any activity in the absence of the gold nanocluster. Atomistic simulations reveal details of the association dynamics as regulated by defined conformational changes of the peptide due to peptide length and sequence. Our results show the importance of a rational design of the peptide to enhance the catalytic activity of peptide–nanoparticle conjugates and present a viable computational approach toward the design of enzyme mimics having a complex structure–function relationship, for technological and nanomedical applications.
Understanding how a macromolecule’s primary sequence governs its conformational landscape is crucial for elucidating its function, yet these design principles are still emerging for macromolecules with intrinsic disorder. Herein, we introduce a high-throughput workflow that implements a practical colorimetric conformational assay, introduces a semi-automated sequencing protocol using matrix-assisted laser desorption/ionization and tandem mass spectrometry (MALDI-MS/MS), and develops a generalizable sequence-structure algorithm. Using a model system of 20mer peptidomimetics containing polar glycine and hydrophobic N-butylglycine residues, we identified nine classifications of conformational disorder and isolated 122 unique sequences across varied compositions and conformations. Conformational distributions of three compositionally identical library sequences were corroborated through atomistic simulations and ion mobility spectrometry coupled with liquid chromatography. A data-driven strategy was developed using existing sequence variables and data-derived “motifs” to inform a machine-learning algorithm toward conformation prediction. Here, this multifaceted approach enhances our understanding of sequence-conformation relationships and offers a powerful tool for accelerating the discovery of materials with conformational control.
This paper presents an optimization framework for sequential reconfiguration using an assortment of switching devices and repair process in distribution system restoration. Compared to existing studies, this paper considers types, capabilities and operational limits of different switching devices, making it applicable in practice. We develop a novel multi-phase method to find the optimal sequential operation of various switching devices and repair faulted areas. We consider circuit breakers, reclosers, sectionalizers, load breaker switches, and fuses. The switching operation problem is decomposed into two mixed-integer linear programming (MILP) subproblems. The first subproblem determines the optimal network topology and estimates the number of steps to reach that topology, while the second subproblem generates a sequence of switching operations to coordinate the switches. For repairing the faults, we design an MILP model that dispatches repair crews to clear faults and replace melted fuses. After clearing a fault, we update the topology of the network by generating a new sequence of switching operations, and the process continues until all faults are cleared. To improve the computational efficiency, a network reduction algorithm is developed to group line sections, such that only switchable sections are present in the reduced network. The proposed method is validated on the IEEE 123-bus and 8500-bus systems.
Underground hydrogen storage (UHS) in geological formations is a key element of the clean energy transition as it enables the decarbonization of the transportation and industrial sectors by decoupling hydrogen production and storage. UHS has many benefits, including low cost, much wider availability, large storage capacity, well-established infrastructure, and increased safety because of geological sealing capabilities. However, the impact of hydrogen (H2) biogeochemical interactions in the presence of subsurface microorganisms is largely neglected from UHS perspectives. These interactions might affect the effectiveness of storage and can even cause H2 to leak into the shallow aquifers. Leakage of H2 into groundwater can change the geochemistry and induce several microbial-driven processes. Microorganisms, such as sulfate-reducers, are naturally abundant in groundwater and consume H2 to produce hydrogen sulfide (H2S), which can contaminate the freshwater drinking groundwater and cause damage to infrastructure. Hydrogen leakage can also trigger microbial reactions responsible for metal mobility, which can impact the water quality. However, the kinetics of these reactions and the temporal impact of hydrogen leakage in groundwater are still unknown. Therefore, a time series hydrogen-groundwater interaction experiment was conducted, and the changes in fluid chemistry and headspace gas composition will be analyzed along with DNA sequencing results to understand the extent and kinetics of biogeochemical reactions that occur if hydrogen leaks into groundwater. In the experiments, Ultra High Purity (UHP) hydrogen gas will be injected into glass vials with groundwater samples, for a designated time period. For each glass-sealed vial, 16S rRNA gene sequencing, IC, ICP-MS, and GC-TCD will be performed. The experiments provide insights into plausible impacts of hydrogen leakage into shallow drinking water aquifers.
Safety analysis simulations have been carried out using the SAS4A/SASSYS-1 safety analysis code system to update the baseline set of results for the Versatile Test Reactor conceptual design. Results from a previous revision of this work provided input for the Conceptual Safety Design Report summarizing the overall safety basis for the VTR. This work will also support future analyses that will be required for Chapter 13, “Accident Analysis,” of a preliminary safety analysis report. A number of transients have been evaluated based on categories of postulated accident sequences suggested for consideration in NUREG-1537. Ultimately, the full spectrum of event initiators and subsequent accident sequences will be defined through a probabilistic risk assessment of the plant design following the implementation of the Licensing Modernization Project guidelines as described in the VTR Safety Design Strategy. Since that process is ongoing, the present analysis addresses several transients that are expected to be bounding accident scenarios that represent the three key ways to perturb a reactor: through changes to the core inlet temperature, mass flow rate, or reactivity. These correspond to a loss of heat sink (LOHS), loss of flow (LOF) or station blackout (SBO), and transient overpower (TOP), respectively. Both Protected and Unprotected versions of these transients have been evaluated, along with additional transients that include transient overpower with one stuck rod, pump coast-down failures, seismic events, and overcooling events. At the current stage of design, transient simulation results for the Versatile Test Reactor indicate that large safety margins exist for many event initiators. However, several “enabling” assumptions have been made, in terms of both design features and design limits, in order to perform the analyses; and these assumptions will need to be revised as the design matures. In addition, a larger spectrum of events needs to be evaluated as part of the on-going probabilistic risk assessment process.
Free-standing two-dimensional (2D) organic nanomaterials are highly attractive for biological applications because of their unique structural properties and high biocompatibility. Herein, we designed and synthesized a new class of highly bright and photostable membrane-mimetic 2D nanosheets from sequence-defined peptoids. These nanosheets exhibited the high quantum yield and photostability as a result of the precise placement and ordering of dansyl dye molecules within crystalline nanosheets. We further showed the use of these nanosheets as biocompatible and programmable probes for live cell imaging and cell labeling. Furthermore, by programing these nanosheets with different surface charges, we achieved the enhanced lysosome escape of these nanosheets, showing their great potential as nanocarriers for the efficient intracellular delivery of macromolecular drugs.
Apoptosis is important for development and tissue homeostasis, and its dysregulation can lead to diseases, including cancer. As an apoptotic effector, BAK undergoes conformational changes that promote mitochondrial outer membrane disruption, leading to cell death. This is termed “activation” and can be induced by peptides from the human proteins BID, BIM, and PUMA. To identify additional peptides that can regulate BAK, we used computational protein design, yeast surface display screening, and structure-based energy scoring to identify 10 diverse new binders. We discovered peptides from the human proteins BNIP5 and PXT1 and three non-native peptides that activate BAK in liposome assays and induce cytochrome c release from mitochondria. Crystal structures and binding studies reveal a high degree of similarity among peptide activators and inhibitors, ruling out a simple function-determining property. Our results shed light on the vast peptide sequence space that can regulate BAK function and will guide the design of BAK-modulating tools and therapeutics.
Microbiomes are extremely important for their host organisms, providing many vital functions and extending their hosts’ phenotypes. Natural studies of host-associated microbiomes can be difficult to interpret due to the high complexity of microbial communities, which hinders our ability to track and identify individual members along with the many factors that structure or perturb those communities. For this reason, researchers have turned to synthetic or constructed communities in which the identities of all members are known. However, due to the lack of tracking methods and the difficulty of creating a more diverse and identifiable community that can be distinguished through next-generation sequencing, most such in vivo studies have used only a few strains. To address this issue, we developed DISCo-microbe, a program for the design of an identifiable synthetic community of microbes for use in in vivo experimentation. The program is composed of two modules; (1) create, which allows the user to generate a highly diverse community list from an input DNA sequence alignment using a custom nucleotide distance algorithm, and (2) subsample, which subsamples the community list to either represent a number of grouping variables, including taxonomic proportions, or to reach a user-specified maximum number of community members. As an example, we demonstrate the generation of a synthetic microbial community that can be distinguished through amplicon sequencing. The synthetic microbial community in this example consisted of 2,122 members from a starting DNA sequence alignment of 10,000 16S rRNA sequences from the Ribosomal Database Project. We generated simulated Illumina sequencing data from the constructed community and demonstrate that DISCo-microbe is capable of designing diverse communities with members distinguishable by amplicon sequencing. Using the simulated data we were able to recover sequences from between 97–100% of community members using two different post-processing workflows. Furthermore, 97–99% of sequences were assigned to a community member with zero sequences being misidentified. We then subsampled the community list using taxonomic proportions to mimic a natural plant host–associated microbiome, ultimately yielding a diverse community of 784 members. DISCo-microbe can create a highly diverse community list of microbes that can be distinguished through 16S rRNA gene sequencing, and has the ability to subsample (i.e., design) the community for the desired number of members and taxonomic proportions. Although developed for bacteria, the program allows for any alignment input from any taxonomic group, making it broadly applicable. The software and data are freely available from GitHub and Python Package Index (PYPI).