Machine Learning Framework for Modeling Exciton-Polaritons in Molecular Materials
Explore the source record for details and available documents.
Material Science↗
SEARCH · Search NASA
Results for “computational tools”
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.
Structure-Based Discovery and Development of Highly Potent Dihydroorotate Dehydrogenase Inhibitors for Malaria Chemoprevention
Malaria remains a serious global health challenge, yet treatment and control programs are threatened by drug resistance. Dihydroorotate dehydrogenase (DHODH) was clinically validated as a target for treatment and prevention of malaria through human studies with DSM265, but currently no drugs against this target are in clinical use. We used structure-based computational tools including free energy perturbation (FEP+) to discover highly ligand efficient, potent, and selective pyrazolebased Plasmodium DHODH inhibitors through a scaffold hop from a pyrrole-based series. Optimized pyrazole-based compounds were identified with low nM-to-pM Plasmodium falciparum cell potency and oral activity in a humanized SCID mouse malaria infection model. The lead compound DSM1465 is more potent and has improved absorption, distribution, metabolism and excretion/pharmacokinetic (ADME/PK) properties compared to DSM265 that support the potential for once-monthly chemoprevention at a low dose. This compound meets the objective of identifying compounds with potential to be used for monthly chemoprevention in Africa to support malaria elimination efforts.
60 APPLIED LIFE SCIENCES↗
Exploring Allosteric Pathways of a V-Type Enzyme with Dynamical Perturbation Networks
Elucidation of the allosteric pathways in proteins is a computational challenge that strongly benefits from combination of atomistic molecular dynamics (MD) simulations and coarse-grained analysis of the complex dynamical network of chemical interactions based on graph theory. Here, we introduce and assess the performances of the dynamical perturbation network analysis of allosteric pathways in a prototypical V-type allosteric enzyme. Dynamical atomic contacts obtained from MD simulations are used to weight the allosteric protein graph, which involves an extended network of contacts perturbed by the effector binding in the allosteric site. The outcome showed good agreement with previously reported theoretical and experimental extended studies and it provided recognition of new potential allosteric spots that can be exploited in future mutagenesis experiments. Lastly, overall, the dynamical perturbation network analysis proved to be a powerful computational tool, complementary to other network-based approaches that can assist the full exploitation of allosteric phenomena for advances in protein engineering and rational drug design.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Effects of Selenium Doping in Zinc Telluride from First Principles
We report point defects and compensation in zinc telluride have been subjects of study for several decades due to the compound’s scientific/technological applications. Despite this, there remains some disagreement in the literature as to the fine details of its defect chemistry. Recently, it has also come under renewed interest due to the emerging Cd( 1–x )Zn x Te( 1–y )Se y material family, used in nuclear radiation detector applications. In this article, we use state-of-the-art computational tools to re-examine the native defect chemistry of ZnTe and explore how it changes upon addition of selenium using density functional theory with screened hybrid exchange–correlation functionals, thermochemical models of the equilibration environment, and statistical mechanics models of the point defect ensemble.
II-VI semiconductors↗
Tunability and Long-Range Enhancement of Resonance Energy Transfer Facilitated by Plasmonic Nanorods
Resonance energy transfer (RET) between molecules or quantum dots is an important process in many energy-related applications. Different environmental structures have been studied and demonstrated to be able to enhance the RET rates between a nearby donor–acceptor pair. In particular, cylindrical silver nanorods and nanowires have shown an extraordinary ability to transfer energy along their longitudinal axes over large distances. However, the detailed mechanism of such transfer and the effects on the molecular RET process are as yet elusive. In this study, we use the recently developed computational tool based on the plasmon-coupled RET method to systematically study the effects of nanorods with different dimensions on the RET rates. We find that highly frequency-dependent coupling factor (CF) spectra, whose amplitudes determine RET rates, can be obtained due to the localized surface plasmon polariton modes of the rods with nanoscale dimensions. Simple phenomenological models can be derived for the wavelengths of CF peaks in relation to the length and width of the nanorods, providing easy tunability for enhancing the RET rate in specific wavelength ranges. When coupled to longer rods with mesoscale lengths, exponential decay of the CF over long donor–acceptor distances with a small decay constant is observed, leading to the possibility of long-range RET processes. Furthermore, drumhead resonance modes emerge on the flat ends of the rod when the rod’s diameter reaches 300 nm, resulting in extra enhancement to RET rate compared to certain thinner rods. Furthermore, these findings shed new light on the mechanism of plasmonic enhancement with silver nanorods and establish design principles for how to optimally utilize these structures to manipulate RET processes for various applications.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
The Challenge of Characterizing High-Concentration Electrolytes at the Molecular Level: A Perspective
High-concentration electrolytes (HCEs) are promising materials composed of highly concentrated salt solutions in organic solvents. HCEs have many desirable properties and are particularly important in the field of batteries. However, the number of ways in which these materials can be tuned is very large, which is crucial for tailored electrolyte design. Moreover, the molecular characterization of HCEs is challenging both experimentally and computationally, but it is necessary for their rational design. Therefore, currently the structure–property–performance relationship of these electrolytes has not been directly derived from their characterization. Here, in this Perspective, we present a brief overview of the HCEs and discuss the state-of-the-art characterization methods used to study them at the molecular level. We also address the challenges associated with these methods, including both experimental techniques and computational tools currently available. Emphasis is placed on methods aimed at understanding the physical phenomena that govern the molecular structure and dynamics occurring on the subnanosecond and nanometer time and length scales. Finally, we discuss new strategies for obtaining a comprehensive characterization of HCEs at the molecular level.
electrolytes↗
Computational Redesign of Acyl-ACP Thioesterase with Improved Selectivity toward Medium-Chain-Length Fatty Acids
Enzyme and metabolic engineering offer the potential to develop biocatalysts for converting natural resources to a wide range of chemicals. To broaden the scope of potential products beyond natural metabolites, methods of engineering enzymes to accept alternative substrates and/or perform novel chemistries must be developed. DNA synthesis can create large libraries of enzyme-coding sequences, but most biochemistries lack a simple assay to screen for promising enzyme variants. Our solution to this challenge is structure-guided mutagenesis, in which optimization algorithms select the best sequences from libraries based on specified criteria (i.e., binding selectivity). We demonstrate this approach by identifying medium-chain (C8–C12) acyl-ACP thioesterases through structure-guided mutagenesis. Medium-chain fatty acids, which are products of thioesterase-catalyzed hydrolysis, are limited in natural abundance, compared to long-chain fatty acids; the limited supply leads to high costs of C6–C10 oleochemicals such as fatty alcohols, amines, and esters. Here, we applied computational tools to tune substrate binding of the highly active ‘TesA thioesterase in Escherichia coli. We used the IPRO algorithm to design thioesterase variants with enhanced C12 or C8 specificity, while maintaining high activity. After four rounds of structure-guided mutagenesis, we identified 3 variants with enhanced production of dodecanoic acid (C12) and 27 variants with enhanced production of octanoic acid (C8). The top variants reached up to 49% C12 and 50% C8 while exceeding native levels of total free fatty acids. A comparably sized library created by random mutagenesis failed to identify promising mutants. The chain length-preference of ‘TesA and the best mutant were confirmed in vitro using acyl-CoA substrates. Molecular dynamics simulations, confirmed by resolved crystal structures, of ‘TesA variants suggest that hydrophobic forces govern ‘TesA substrate specificity. Finally, we expect the design rules that we uncovered and the thioesterase variants that we identified will be useful to metabolic engineering projects aimed at sustainable production of medium-chain-length oleochemicals.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Beyond Local Solvation Structure: Nanometric Aggregates in Battery Electrolytes and their Effect on Electrolyte Properties
Electrolytes are an essential component of all electrochemical storage and conversion devices, such as batteries. In the history of battery development, the complex nature of electrolytes has often been a bottleneck. Fundamental knowledge of electrolyte systems encompasses elucidation of structure-property relationships of the solution species. Recently, nanometric aggregates have been observed in several classes of electrolytes, including super-concentrated, redox-flow, multivalent, polymer, and ionic liquid-based electrolytes. Compared with the well-studied local solvation structures such as contact ion pairs and solvent-separated ions, these aggregates impose unique effects on the ion distribution and transport both within bulk electrolytes and at electrode/electrolyte interfaces. This Perspective highlights the discovery of the aggregates in various battery electrolytes and their impact on electrolyte properties. We also present an outlook for future studies of this emerging field of nanometric aggregates and the need for the development of new experimental and computational tools to study their properties.
Yu, Zhou↗
Dual Enhancement of Thermostability and Activity of Xylanase through Computer-Aided Rational Design
In the realm of enzyme engineering, the dual enhancement of thermostability and activity remains a challenge. Herein, we employed a computer-aided approach integrating folding free energy calculations and evolutionary analysis to engineer Paecilomyces thermophila xylanase into a hyperthermophilic enzyme for application in the paper and pulp industry. Through the computational rational design, XynM9 with superior thermostability and enhanced activity was designed. Its optimal reaction temperature increases by 10 °C to 85 °C, its T m increases by 10 °C to 93 °C, and its half-life increases 11-fold to 5.8 h. Additionally, its catalytic efficiency improves by 57% to 3926 s –1 mM –1 . Molecular dynamics simulations revealed that XynM9 is stabilized by more hydrogen bonds and salt bridges than wild-type xylanase. The mutant’s narrower catalytic cleft enhances the substrate-binding affinity, thus improving the catalytic efficiency. In harsh conditions at 80 °C and pH 10, using XynM9 significantly reduced both hemicellulose and lignin, which makes it a good candidate for use in the paper and pulp process. Finally, our study presents an accurate and efficient strategy for the dual enhancement of enzyme properties, guiding further improvement of computational tools for protein stabilization.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Coarse-Grained Modeling Elucidates Differential Metabolism of Saccharomyces cerevisiae under Varied Nutrient Limitations
Microorganisms such as Saccharomyces cerevisiae have a native ability to adapt their metabolism to varying nutrient conditions. Understanding their responses to nutrient limitations is critical for decoding cellular physiology and designing strategies for metabolic engineering. While the influence of carbon availability on yeast metabolism has been extensively studied, the role of nitrogen availability remains relatively underexplored. In this study, we utilized a coarse-grained kinetic model to systematically analyze and compare the effects of carbon and nitrogen limitations on yeast metabolism. Our model successfully revealed the differential metabolic characteristics of S. cerevisiae under carbon- and nitrogen-limited chemostat conditions. It also highlighted the significance of protein activity regulation at varying carbon-to-nitrogen ratios, and elucidated distinct strategies employed to maintain ATP homeostasis. Furthermore, this study provides a computational tool for investigating yeast physiology under nutrient limitations and offers quantitative and mechanistic insights into yeast metabolism.
59 BASIC BIOLOGICAL SCIENCES↗
Identification of Ultramodified Proteins Using Top-Down Mass Spectra
Post-translational modifications (PTMs) play an important role in various biological processes through changing protein structure and function. Some ultramodified proteins (like histones) have multiple PTMs forming PTM patterns that define the functionality of a protein. While bottom-up mass spectrometry (MS) has been successful in identifying individual PTMs within short peptides, it is unable to identify PTM patterns spread along entire proteins in a coordinated fashion. In contrast, top-down MS analyzes intact proteins and reveals PTM patterns along the entire proteins. However, while recent advances in instrumentation have made top-down MS accessible to many laboratories, most computational tools for top-down MS focus on proteins with few PTMs and are unable to identify complex PTM patterns. We propose a new algorithm, MS-Align-E, that identifies both expected and unexpected PTMs in ultramodified proteins. We demonstrate that MS-Align-E identifies many protein forms of histone H4 and benchmark it against the currently accepted software tools.
Environmental Molecular Sciences Laboratory↗
Baseflow Identification via Explainable AI With Kolmogorov‐Arnold Networks
Abstract Hydrological models often involve constitutive laws that may not be optimal in every application. We propose to replace such laws with the Kolmogorov‐Arnold networks (KANs), a class of neural networks designed to identify symbolic expressions. We demonstrate KAN's potential on the problem of baseflow identification, a notoriously challenging task plagued by significant uncertainty. KAN‐derived functional dependencies of the baseflow components on the aridity index outperform their original counterparts; they demonstrate that water availability, rather than potential evapotranspiration, drives baseflow by constraining actual evapotranspiration under arid conditions. On a test set, they increase the Nash‐Sutcliffe efficiency (NSE) by 65%, decrease the root mean squared error by 29%, and increase the Kling‐Gupta efficiency by 34%. This superior performance is achieved while reducing the number of fitting parameters from three to two. Next, we use data from 378 catchments across the continental United States to refine the water‐balance equation at the mean‐annual scale. The KAN‐derived equations based on the refined water balance outperform both the current aridity index model, with up to a 105% increase in NSE, and the KAN‐derived equations based on the original water balance. While the performance of our model and tree‐based machine learning methods is similar, KANs offer the advantage of simplicity and transparency and require no specific software or computational tools. This case study focuses on the aridity index formulation, but the approach is flexible and transferable to other hydrological processes. Plain Language Summary Equations used in hydrologic model are often suboptimal, resulting in reduced prediction accuracy and efficiency. We implemented Kolmogorov‐Arnold networks (KAN), a machine learning algorithm for deriving symbolic formulations, to estimate groundwater recharge and showed that it outperforms an existing state‐of‐the‐art semi‐empirical formulation. In hydrology, Nash‐Sutcliffe efficiency (NSE), root mean squared error (RMSE), and Kling‐Gupta efficiency (KGE) are commonly used to evaluate model performance. Higher NSE and KGE values indicate better performance, while lower RMSE values are preferable. Our results show that NSE increased by 71%, RMSE decreased by 32%, and KGE improved by 25%. In addition, KAN identifies an optimal functional form and can be used to derive new analytical formulas using the prior knowledge. The KAN‐inspired equation outperformed the original formulation and reduced the fitting parameters. Furthermore, we refined the water‐balance equation at the mean‐annual scale and showed that, based on the new water‐balance equation, KAN can derive new formulations that are superior to the original aridity index formulations (up to 105% increase in NSE) and KAN‐derived equations based on the original water balance. These findings highlight the significant potential of KAN to advance the scientific understanding of a wide range of hydrologic processes. Key Points Kolmogorov‐Arnold networks (KANs) enhance interpretability of machine‐learned hydrological models KAN‐derived symbolic formulations outperform state‐of‐the‐art semi‐empirical aridity indices KAN‐identified functional form yields an analytical index with fewer fitting parameters and improved performance
baseflow↗
Establishing the precise evolutionary history of a gene improves prediction of disease-causing missense mutations
Here, predicting the phenotypic effects of mutations has become an important application in clinical genetic diagnostics. Computational tools evaluate the behavior of the variant over evolutionary time and assume that variations seen during the course of evolution are probably benign in humans. However, current tools do not take into account orthologous/paralogous relationships. Paralogs have dramatically different roles in Mendelian diseases. For example, whereas inactivating mutations in the NPC1 gene cause the neurodegenerative disorder Niemann-Pick C, inactivating mutations in its paralog NPC1L1 are not disease-causing and, moreover, are implicated in protection from coronary heart disease. Methods: We identified major events in NPC1 evolution and revealed and compared orthologs and paralogs of the human NPC1 gene through phylogenetic and protein sequence analyses. We predicted whether an amino acid substitution affects protein function by reducing the organism s fitness. As a result, removing the paralogs and distant homologs improved the overall performance of categorizing disease-causing and benign amino acid substitutions. In conclusion, the results show that a thorough evolutionary analysis followed by identification of orthologs improves the accuracy in predicting disease-causing missense mutations. We anticipate that this approach will be used as a reference in the interpretation of variants in other genetic diseases as well.
59 BASIC BIOLOGICAL SCIENCES↗
Single-cell genomics-based analysis of virus–host interactions in marine surface bacterioplankton
Viral infections dynamically alter the composition and metabolic potential of marine microbial communities and the evolutionary trajectories of host populations with resulting feedback on biogeochemical cycles. It is quite possible that all microbial populations in the ocean are impacted by viral infections. Our knowledge of virus–host relationships, however, has been limited to a minute fraction of cultivated host groups. Here, we utilized single-cell sequencing to obtain genomic blueprints of viruses inside or attached to individual bacterial and archaeal cells captured in their native environment, circumventing the need for host and virus cultivation. Furthermore, a combination of comparative genomics, metagenomic fragment recruitment, sequence anomalies and irregularities in sequence coverage depth and genome recovery were utilized to detect viruses and to decipher modes of virus–host interactions. Members of all three tailed phage families were identified in 20 out of 58 phylogenetically and geographically diverse single amplified genomes (SAGs) of marine bacteria and archaea. At least four phage–host interactions had the characteristics of late lytic infections, all of which were found in metabolically active cells. One virus had genetic potential for lysogeny. Our findings include first known viruses of Thaumarchaeota, Marinimicrobia, Verrucomicrobia and Gammaproteobacteria clusters SAR86 and SAR92. Viruses were also found in SAGs of Alphaproteobacteria and Bacteroidetes. A high fragment recruitment of viral metagenomic reads confirmed that most of the SAG-associated viruses are abundant in the ocean. This study demonstrates that single-cell genomics, in conjunction with sequence-based computational tools, enable in situ, cultivation-independent insights into host–virus interactions in complex microbial communities.
58 GEOSCIENCES↗
Microbiome-wide association studies link dynamic microbial consortia to disease
Rapid advances in DNA sequencing, metabolomics, proteomics and computational tools are dramatically increasing access to the microbiome and identification of its links with disease. In particular, time-series studies and multiple molecular perspectives are facilitating microbiome-wide association studies, which are analogous to genome-wide association studies. Early findings point to actionable outcomes of microbiome-wide association studies, although their clinical application has yet to be approved. An appreciation of the complexity of interactions among the microbiome and the host's diet, chemistry and health, as well as determining the frequency of observations that are needed to capture and integrate this dynamic interface, is paramount for developing precision diagnostics and therapies that are based on the microbiome.
Gilbert, Jack A.↗
Tracking HIV-1 recombination to resolve its contribution to HIV-1 evolution in natural infection
Recombination in HIV-1 is well documented, but its importance in the low-diversity setting of within-host diversification is less understood. Here we develop a novel computational tool (RAPR (Recombination Analysis PRogram)) to enable a detailed view of in vivo viral recombination during early infection, and we apply it to near-full-length HIV-1 genome sequences from longitudinal samples. Recombinant genomes rapidly replace transmitted/founder (T/F) lineages, with a median half-time of 27 days, increasing the genetic complexity of the viral population. We identify recombination hot and cold spots that differ from those observed in inter-subtype recombinants. Furthermore, RAPR analysis of longitudinal samples from an individual with well-characterized neutralizing antibody responses shows that recombination helps carry forward resistance-conferring mutations in the diversifying quasispecies. These findings provide insight into molecular mechanisms by which viral recombination contributes to HIV-1 persistence and immunopathogenesis and have implications for studies of HIV transmission and evolution in vivo.
59 BASIC BIOLOGICAL SCIENCES↗
Multiplex base editing to convert $\mathrm{TAG}$ into $\mathrm{TAA}$ codons in the human genome
Whole-genome recoding has been shown to enable nonstandard amino acids, biocontainment and viral resistance in bacteria. Here we take the first steps to extend this to human cells demonstrating exceptional base editing to convert TAG to TAA for 33 essential genes via a single transfection, and examine base-editing genome-wide (observing ~40 C-to-T off-target events in essential gene exons). We also introduce GRIT, a computational tool for recoding. This demonstrates the feasibility of recoding, and highly multiplex editing in mammalian cells.
59 BASIC BIOLOGICAL SCIENCES↗
An updated nuclear-physics and multi-messenger astrophysics framework for binary neutron star mergers
The multi-messenger detection of the gravitational-wave signal GW170817, the corresponding kilonova AT2017gfo and the short gamma-ray burst GRB170817A, as well as the observed afterglow has delivered a scientific breakthrough. For an accurate interpretation of all these different messengers, one requires robust theoretical models that describe the emitted gravitational-wave, the electromagnetic emission, and dense matter reliably. In addition, one needs efficient and accurate computational tools to ensure a correct cross-correlation between the models and the observational data. For this purpose, we have developed the Nuclear-physics and Multi-Messenger Astrophysics framework NMMA. The code allows incorporation of nuclear-physics constraints at low densities as well as X-ray and radio observations of isolated neutron stars. In previous works, the NMMA code has allowed us to constrain the equation of state of supranuclear dense matter, to measure the Hubble constant, and to compare dense-matter physics probed in neutron-star mergers and in heavy-ion collisions, and to classify electromagnetic observations and perform model selection. Here, we show an extension of the NMMA code as a first attempt of analyzing the gravitational-wave signal, the kilonova, and the gamma-ray burst afterglow simultaneously. Incorporating all available information, we estimate the radius of a 1.4M ⊙ neutron star to be $R$ $=$ $11.98^{+0.35}_{-0.04}$ km.
79 ASTRONOMY AND ASTROPHYSICS↗