Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Homology modelling”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Analysis of meiosis in Pristionchus pacificus reveals plasticity in homolog pairing and synapsis in the nematode lineage

Meiosis is conserved across eukaryotes yet varies in the details of its execution. Here we describe a new comparative model system for molecular analysis of meiosis, the nematode Pristionchus pacificus, a distant relative of the widely studied model organism Caenorhabditis elegans. P. pacificus shares many anatomical and other features that facilitate analysis of meiosis in C. elegans. However, while C. elegans has lost the meiosis-specific recombinase Dmc1 and evolved a recombination-independent mechanism to synapse its chromosomes, P. pacificus expresses both DMC-1 and RAD-51. We find that SPO-11 and DMC-1 are required for stable homolog pairing, synapsis, and crossover formation, while RAD-51 is dispensable for these key meiotic processes. RAD-51 and DMC-1 localize sequentially to chromosomes during meiotic prophase and show nonoverlapping functions. We also present a new genetic map for P. pacificus that reveals a crossover landscape very similar to that of C. elegans, despite marked divergence in the regulation of synapsis and crossing-over between these lineages.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Initiation of fatty acid biosynthesis in Pseudomonas putida $\mathrm{KT2440}$

Deciphering the mechanisms of bacterial fatty acid biosynthesis is crucial for both the engineering of bacterial hosts to produce fatty acid-derived molecules and the development of new antibiotics. However, gaps in our understanding of the initiation of fatty acid biosynthesis remain. Here, we demonstrate that the industrially relevant microbe Pseudomonas putida KT2440 contains three distinct pathways to initiate fatty acid biosynthesis. The first two routes employ conventional β-ketoacyl-ACP synthase III enzymes, FabH1 and FabH2, that accept short- and medium-chain-length acyl-CoAs, respectively. The third route utilizes a malonyl-ACP decarboxylase enzyme, MadB. A combination of exhaustive in vivo alanine-scanning mutagenesis, in vitro biochemical characterization, X-ray crystallography, and computational modeling elucidate the presumptive mechanism of malonyl-ACP decarboxylation via MadB. Given that functional homologs of MadB are widespread throughout domain Bacteria, this ubiquitous alternative fatty acid initiation pathway provides new opportunities to target a range of biotechnology and biomedical applications.

09 BIOMASS FUELS↗
DOE OSTI2025

ThinCurr: An open-source 3D thin-wall eddy current modeling code for the analysis of large-scale systems of conducting structures

In this paper we present a new thin-wall eddy current modeling code, ThinCurr, for studying inductively-coupled currents in 3D conducting structures -- with primary application focused on the interaction between currents flowing in coils, plasma, and conducting structures of magnetically-confined plasma devices. The code utilizes a boundary finite element method on an unstructured, triangular grid to accurately capture device structures. The new code, part of the broader Open FUSION Toolkit, is open-source and designed for ease of use without sacrificing capability and speed through a combination of Python, Fortran, and C/C++ components. Scalability to large models is enabled through use of hierarchical off-diagonal low-rank compression of the inductance matrix, which is otherwise dense. Ease of handling large models of complicated geometry is further supported by automatic determination of supplemental elements through a greedy homology approach. Here, a detailed description of the numerical methods of the code and verification of the implementation of those methods using cross-code comparisons against the VALEN code and Ansys commercial analysis software is shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Learning-based approach to plasticity in athermal sheared amorphous packings: Improving softness

The plasticity of amorphous solids undergoing shear is characterized by quasi-localized rearrangements of particles. While many models of plasticity exist, the precise relationship between the plastic dynamics and the structure of a particle’s local environment remains an open question. Previously, machine learning was used to identify a structural predictor of rearrangements called “softness.” Although softness has been shown to predict which particles will rearrange with high accuracy, the method can be difficult to implement in experiments where data are limited and the combinations of descriptors it identifies are often difficult to interpret physically. Here, we address both of these weaknesses, presenting two major improvements to the standard softness method. First, we present a natural representation of each particle’s observed mobility, allowing for the use of statistical models that are both simpler and provide greater accuracy in limited datasets. Second, we employ persistent homology as a systematic means of identifying simple, topologically informed, structural quantities that are easy to interpret and measure experimentally. We test our methods on two-dimensional athermal packings of soft spheres under quasi-static shear. We find that the same structural information that predicts small variations in the response is also predictive of where plastic events will localize. We also find that an excellent accuracy is achieved in athermal sheared packings using simply a particle’s species and the number of nearest neighbor contacts.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Structural basis of the amidase ClbL central to the biosynthesis of the genotoxin colibactin

Colibactin is a genotoxic natural product produced by select commensal bacteria in the human gut microbiota. The compound is a bis-electrophile that is predicted to form interstrand DNA cross-links in target cells, leading to double-strand DNA breaks. The biosynthesis of colibactin is carried out by a mixed NRPS–PKS assembly line with several noncanonical features. An amidase, ClbL, plays a key role in the pathway, catalyzing the final step in the formation of the pseudodimeric scaffold. ClbL couples α-aminoketone and β-ketothioester intermediates attached to separate carrier domains on the NRPS–PKS assembly. Here, the 1.9 Å resolution structure of ClbL is reported, providing a structural basis for this key step in the colibactin biosynthetic pathway. The structure reveals an open hydrophobic active site surrounded by flexible loops, and comparison with homologous amidases supports its unusual function and predicts macromolecular interactions with pathway carrier-protein substrates. Modeling protein–protein interactions supports a predicted molecular basis for enzyme–carrier domain interactions. Overall, the work provides structural insight into this unique enzyme that is central to the biosynthesis of colibactin.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Transport of metformin metabolites by guanidinium exporters of the small multidrug resistance family

Proteins from the small multidrug resistance (SMR) family are frequently associated with horizontally transferred multidrug resistance gene arrays found in bacteria from wastewater and the human-adjacent biosphere. Recent studies suggest that a subset of SMR transporters might participate in the metabolism of the common pharmaceutical metformin by bacterial consortia. Here, we show that both genomic and plasmid-associated transporters of the SMRGdx functional subtype export byproducts of microbial metformin metabolism, with particularly high export efficiency for guanylurea. We use solid-supported membrane electrophysiology to evaluate the transport kinetics for guanylurea and native substrate guanidinium by four representative SMRGdx homologs. Using an internal reference to normalize independent electrophysiology experiments, we show that transport rates are comparable for genomic and plasmid-associated SMRGdx homologs, and using a proteoliposome-based transport assay, we show that 2 proton:1 substrate transport stoichiometry is maintained. Additional characterization of guanidinium and guanylurea export properties focuses on the structurally characterized homolog, Gdx-Clo, for which we examined the pH dependence and thermodynamics of substrate binding and solved an x-ray crystal structure with guanylurea bound. Together, these experiments contribute in two main ways. By providing the first detailed kinetic examination of the structurally characterized SMRGdx homolog Gdx-Clo, they provide a functional framework that will inform future mechanistic studies of this model transport protein. Second, this study casts light on a potential role for SMRGdx transporters in microbial handling of metformin and its microbial metabolic byproducts, providing insight into how native transport physiologies are co-opted to contend with new selective pressures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Anomalies in the topology of the temperature fluctuations in the cosmic microwave background: An analysis of the NPIPE and FFP10 data releases

We present a topological analysis of the temperature fluctuation maps from the Planck 2020 Data Release 4 NPIPE dataset and the Planck 2018 Data Release 3 FFP10 dataset. We performed a multiscale analysis in terms of the homology characteristics of the maps, invoking relative homology to account for the analysis in the presence of masks. We performed our analysis for a range of smoothing scales spanning sub- and super-horizon scales corresponding to a full width at half maximum (FWHM) of 5',10',20',40',80',160',320', and 640', and employed simulations based on the standard model for comparison, which assumes the initial fluctuation field to be an isotropic and homogeneous Gaussian random field. Examining the behavior of topological components, represented by the 0D homology group, we find the observations to be approximately 2σ or less deviant from the simulations for all resolutions and scales for the NPIPE dataset. For the FFP10 dataset, we detect a 2.96σ deviation between the observations and simulations at N = 128, FWHM = 80'. For the topological loops, represented by the first homology group, the simulations and observations are consistent within 2σ for most resolutions and scales for both the datasets. However, for the NPIPE dataset, we observe a high deviation between the observation and simulations in the number of loops at FWHM = 320', but at a low dimensionless threshold ν = –2.5. Under a Gaussian assumption, this would amount to a deviation of ~4σ. However, the distribution in this bin is manifestly non-Gaussian and does not obey Poisson statistics either. In the absence of a true theoretical understanding, we simply note that the significance is higher than what may be resolved by 600 simulations, yielding an empirical p-value of at most 0.0016. Specifically in this case, our tests indicate that the numbers arise from a statistically stable regime, despite being based on small numbers. For the FFP10 dataset, the differences are not as strong as for the NPIPE dataset, indicating a 2.77σ deviation at this resolution and threshold. The Euler characteristic, which is the alternating sum of the ranks of relative homology groups, reflects the deviations in the components and loops. To assess the significance of combined levels for a given scale, we employed the empirical and theoretical versions of the χ 2 test as well as the nonparametric Tukey depth test. Although all statistics exhibit a stable distribution, we favor the empirical version of the χ 2 test in the final interpretation, as it indicates the most conservative differences. For the NPIPE dataset, we find that the components and loops differ at more than 95%, but agree within the 99% confidence level with respect to the base model at N = 32, FWHM = 320'. The Euler characteristic at this resolution displays a per mil deviation. In contrast, the FFP10 dataset shows that the observations are consistent with the base model within the 95% confidence level, at this and smaller scales. This is consistent with the observations of the Planck analysis pipeline via Minkowski functionals. For the largest smoothing scale, N = 16, FWHM = 640', both datasets exhibit an anomalous behavior of the loops, where FFP10 data exhibit a deviation that is larger by an order of magnitude than that of the NPIPE dataset. In contrast, the values for the topological components and the Euler characteristic agree between observations and model to within a confidence level of 99%. However, for the largest scales, the statistics are based on low numbers and may have to be regarded with caution. Even though both datasets exhibit mild to significant discrepancies, they also exhibit contrasting behaviors at various instances. Therefore, we do not find it feasible to convincingly accept or reject the null hypothesis. Disregarding the large-scale anomalies that persist at similar scales in WMAP and Planck, observations of the cosmic microwave background are largely consistent with the standard cosmological model within 2σ.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Functional protein mining with conformal guarantees

Molecular structure prediction and homology detection offer promising paths to discovering protein function and evolutionary relationships. However, current approaches lack statistical reliability assurances, limiting their practical utility for selecting proteins for further experimental and in-silico characterization. To address this challenge, we introduce a statistically principled approach to protein search leveraging principles from conformal prediction, offering a framework that ensures statistical guarantees with user-specified risk and provides calibrated probabilities (rather than raw ML scores) for any protein search model. Our method (1) lets users select many biologically-relevant loss metrics (i.e. false discovery rate) and assigns reliable functional probabilities for annotating genes of unknown function; (2) achieves state-of-the-art performance in enzyme classification without training new models; and (3) robustly and rapidly pre-filters proteins for computationally intensive structural alignment algorithms. Our framework enhances the reliability of protein homology detection and enables the discovery of uncharacterized proteins with likely desirable functional properties.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2020

Shedding Light on Microbial Dark Matter with A Universal Language of Life

The majority of microbial genomes have yet to be cultured, and most proteins predicted from microbial genomes or sequenced from the environment cannot be functionally annotated. As a result, current computational approaches to describe microbial systems rely on incomplete reference databases that cannot adequately capture the full functional diversity of the microbial tree of life, limiting our ability to model high-level features of biological sequences. The scientific community needs a means to capture the functionally and evolutionarily relevant features underlying biology, independent of our incomplete reference databases. Such a model can form the basis for transfer learning tasks, enabling downstream applications in environmental microbiology, medicine, and bioengineering. Here we present LookingGlass, a deep learning model capturing a “universal language of life”. LookingGlass encodes contextually-aware, functionally and evolutionarily relevant representations of short DNA reads, distinguishing reads of disparate function, homology, and environmental origin. We demonstrate the ability of LookingGlass to be fine-tuned to perform a range of diverse tasks: to identify novel oxidoreductases, to predict enzyme optimal temperature, and to recognize the reading frames of DNA sequence fragments. LookingGlass is the first contextually-aware, general purpose pre-trained “biological language” representation model for short-read DNA sequences. LookingGlass enables functionally relevant representations of otherwise unknown and unannotated sequences, shedding light on the microbial dark matter that dominates life on Earth.

A Hoarfrost↗
DOE OSTI2024

Electron-Withdrawing meso -Substituents Turn On Magneto-Optical Activity in Porphyrins

A series of square planar metalloporphyrins (M(TPP), TPP is 5,10,15,20-tetraphenylporphyrin and M(TPFPP), TPFPP is 5,10,15,20-tetrapentafluorophenylporphyrin; M is Zn 2+ , Ni 2+ , Pd 2+ , or Pt 2+ ) with distinct meso-substituents were prepared, and their magneto-optical activity (MOA) was characterized by magnetic circular dichroism (MCD) and magneto-optical rotary dispersion spectroscopy (MORD; also known as Faraday rotation spectroscopy). MOA is crucial in the development of next-generation magneto-optical devices and quantum computing. Here, the data show that the presence of meso-pentafluorophenyl substituents results in significant increase in MOA in comparison to the homologous phenyl group. Differences in the MOA of these metalloporphyrins are rationalized using the Gouterman four-orbital model and pave the way for rational design of improved and tailorable magneto-optical materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1976

H alpha observations of the 12 August 1975 type 3-RS bursts

H alpha filtergram observations of a number of the Type III-RS (reverse slope) bursts that occurred on August 12, 1975 are presented. Solar radio emission was peculiar on that date in that a large number, and proportion, of the usually rare reverse slope bursts were observed. The radio bursts are shown to coincide in time with a homologous set of H alpha flares located at the limbward edge of spot group Mt. Wilson 19598. A model is proposed in which the reverse slope bursts are the downward branches of U bursts, whose upward branches are hidden behind the coronal density enhancement over the spot group.

Labonte, B. J.↗
DOE OSTI2023

Structural and Biophysical Properties of a [4Fe-4S] Ferredoxin-Like Protein from Synechocystis sp. PCC 6803 with a Unique Two Domain Structure

Electron carrier proteins (ECPs), binding iron-sulfur clusters, are vital components within the intricate network of metabolic and photosynthetic reactions. They play a crucial role in the distribution of reducing equivalents. In Synechocystis sp. PCC 6803, the ECP network includes at least nine ferredoxins. Previous research, including global expression analyses and protein binding studies, has offered initial insights into the functional roles of individual ferredoxins within this network. This study primarily focuses on Ferredoxin 9 (slr2059). Through sequence analysis and computational modeling, Ferredoxin 9 emerges as a unique ECP with a distinctive two-domain architecture. It consists of a C-terminal iron-sulfur binding domain and an N-terminal domain with homology to Nil-domain proteins, connected by a structurally rigid 4-amino acid linker. Notably, in contrast to canonical [2Fe-2S] ferredoxins exemplified by PetF (ssl0020), which feature highly acidic surfaces facilitating electron transfer with photosystem I reaction centers, models of Ferredoxin 9 reveal a more neutral to basic protein surface. Using a combination of electron paramagnetic resonance spectroscopy and square-wave voltammetry on heterologously produced Ferredoxin 9, this study demonstrates that the protein coordinates 2x[4Fe-4S]2+/1+ redox-active and magnetically interacting clusters, with measured redox potentials of -420 +/- 9 mV and -516 +/- 10 mV vs SHE. A more in-depth analysis of Fdx9's unique structure and protein sequence suggests that this type of Nil-2[4Fe-4S] multi-domain ferredoxin is well conserved in cyanobacteria, bearing structural similarities to proteins involved in homocysteine synthesis in methanogens.

cyanobacteria↗
DOE OSTI2022

A separations and purification process for improving yields and meeting fuel contaminant specifications for high-octane gasoline produced from dimethyl-ether over a Cu/BEA catalyst

In this work, we have been developing a three-step conversion of biomass-derived syngas to methanol to dimethyl-ether (DME) to non-aromatic hydrocarbons for use as high-octane gasoline and sustainable aviation fuel. This process produces branched alkanes from DME using a Cu/BEA catalyst and is a promising alternative to other syngas conversion processes such as Fischer-Tropsch to linear alkanes and traditional ZSM-5 catalyzed methanol to aromatic gasoline. In this short article we describe some advances in our understanding related to separations and purification via the use of more detailed experimental speciation in an updated process model involving multiple phase equilibrium-based separation steps. Primary modeled reactor outlet constituents (and weight %) are: C3 and lighter hydrocarbon gases (11.1%), C4s (54.5%), H 2 (1.2%), CO 2 (2.9%), water (5.0%), unreacted DME (16.5%), methanol (2.3%), and C5+ hydrocarbons (6.4%). DME (the primary reactant) and H 2 recycle and reuse are important for the overall process efficiency, and the recycle of C4s is important to increase the C5+ yield via reactivation and homologation. Thus H 2 , C4s, and DME are targeted for recycle, while methanol and water need to be removed from the product to conform with fuel specifications. Model predictions from Aspen Plus using the NRTL-RK property method indicate a fuel composition with C5+ content of 97.1 wt%, with minor constituents: 2.4 wt% C4s, 0.3 wt% methanol, 0.1 wt% DME, 0.03 wt% water, and 0.01 wt% C3s. These ranges of minor components conform with fuel quality requirements, and the modeled product is amenable for unconstrained blending to boost gasoline octane ratings.

09 BIOMASS FUELS↗
NASA NTRS1994

Analysis of the early spectra and light curve of SN 1987A

Numerical modeling of supernova spectra, light curves, and hydrodynamics requires physical inputs, numerical techniques, approximations, and assumptions which must be thoroughly understood in order to study the details of supernova explosions. Here, we discuss some of these in the context of the early evolution of supernova 1987A. Gray radiation-hydrodynamics is used to calculate the bolometric light curve and the hydrodynamic evolution of the supernova. Synthetic spectra are then obtained for the resulting density and velocity structure. The spectrum calculations are performed using a special-relativistic treatment of the radiative transfer equation in the comoving frame, line blanketing by about 10(exp 5) spectral lines, and departures from local thermodynamic equilibrium (LTE) for H I, He I, Mg II, and Ca II. We find that we are able to simultaneously fit the early light curve and spectra reasonably well, using a progenitor model from Arnett (1991a), without fine-tuning the free parameters. Temperature structures and radiative equilibrium, non-LTE effects, homologous expansion, and mean opacities are discussed.

Hauschildt, Peter H.↗
NASA NTRS2002

Predicting functional divergence in protein evolution by site-specific rate shifts

Most modern tools that analyze protein evolution allow individual sites to mutate at constant rates over the history of the protein family. However, Walter Fitch observed in the 1970s that, if a protein changes its function, the mutability of individual sites might also change. This observation is captured in the "non-homogeneous gamma model", which extracts functional information from gene families by examining the different rates at which individual sites evolve. This model has recently been coupled with structural and molecular biology to identify sites that are likely to be involved in changing function within the gene family. Applying this to multiple gene families highlights the widespread divergence of functional behavior among proteins to generate paralogs and orthologs.

Review↗
DOE OSTI2025

Bottom-Up Simulation, Reconstruction, and Quantification of Macromolecule Sequences from Experimental Polymerizations

Motivated by the canonical sequence–structure–function paradigm, tools to characterize chemical patterning in natural biomacromolecules, from proteins to nucleic acids, have grown exponentially in recent years. However, analogous strategies for synthetic macromolecules remain in nascent stages, complicated by sequence polydispersity and analytical limitations. To address this, we have developed a comprehensive and open-source Python package, PRISM (polymer rate insights and sequence modeling), an end-to-end workflow that provides a path from experimental kinetics measurements to quantitative and qualitative metrics for describing chemical patterning in stochastic polymers. First, a numerical integration strategy was constructed to simulate and fit experimental data from reversible addition–fragmentation chain transfer (RAFT) polymerization kinetics, enabling the facile estimation of relevant reactivity ratios. These ratios were then used in a mechanism-specific stochastic kinetic simulation strategy to simulate sequence ensembles corresponding to model systems spanning experimental copolymers, classes of statistical polymers (e.g., alternating, block, and gradient), and multiblock copolymers. Lastly, inspired by sequence homology metrics from bioinformatics, we introduce visualization strategies and quantitative metrics to facilitate comparisons of different sequence ensembles. As the sequence–structure–function paradigm becomes increasingly central in de novo design of synthetic macromolecules, this toolkit provides a first step toward accurate and representative sequence description and featurization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1976

H-alpha observations of the August 12, 1975 Type III-RS bursts

We present H-alpha filtergram observations of a number of the Type III-RS (reverse slope) bursts that occurred on August 12, 1975. Solar radio emission was peculiar on that date in that a large number, and proportion, of the usually rare reverse slope bursts were observed (Tarnstrom and Zehntner, 1975). We show that the radio bursts coincide in time with a homologous set of H-alpha flares located at the limbward edge of spot group Mt. Wilson 19598. We propose a model in which the reverse slope bursts are the downward branches of U bursts, whose upward branches are hidden behind the coronal density enhancement over the spot group.

Labonte, B. J.↗