Search NASA⌕ Search

SEARCH · Search NASA

Results for “Molecular biophysics”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Editorial: Advances in computational molecular biophysics

This special issue of BBA is dedicated to Dr. Jeremy C. Smith on the occasion of his 60th birthday and his more than three decades as a practitioner of computational molecular biophysics. The special issue showcases cutting-edge work in computational molecular biophysics. The senior authors of the manuscripts in this issue are Dr. Smith’s previous or present collaborators. Their articles reflect Dr. Smith’s distinguished career integrating experimental and molecular simulation approaches into a unified description of the function, structure, and dynamics of biomolecules. Each of the articles presented here illustrates one aspect of the broad nature of computational molecular biophysics.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

OpenAWSEM with Open3SPN2: A fast, flexible, and accessible framework for large-scale coarse-grained biomolecular simulations

We present OpenAWSEM and Open3SPN2, new cross-compatible implementations of coarse-grained models for protein (AWSEM) and DNA (3SPN2) molecular dynamics simulations within the OpenMM framework. These new implementations retain the chemical accuracy and intrinsic efficiency of the original models while adding GPU acceleration and the ease of forcefield modification provided by OpenMM’s Custom Forces software framework. By utilizing GPUs, we achieve around a 30-fold speedup in protein and protein-DNA simulations over the existing LAMMPS-based implementations running on a single CPU core. We showcase the benefits of OpenMM’s Custom Forces framework by devising and implementing two new potentials that allow us to address important aspects of protein folding and structure prediction and by testing the ability of the combined OpenAWSEM and Open3SPN2 to model protein-DNA binding. The first potential is used to describe the changes in effective interactions that occur as a protein becomes partially buried in a membrane. We also introduced an interaction to describe proteins with multiple disulfide bonds. Using simple pairwise disulfide bonding terms results in unphysical clustering of cysteine residues, posing a problem when simulating the folding of proteins with many cysteines. We now can computationally reproduce Anfinsen’s early Nobel prize winning experiments by using OpenMM’s Custom Forces framework to introduce a multi-body disulfide bonding term that prevents unphysical clustering. Our protein-DNA simulations show that the binding landscape is funneled towards structures that are quite similar to those found using experiments. In summary, this paper provides a simulation tool for the molecular biophysics community that is both easy to use and sufficiently efficient to simulate large proteins and large protein-DNA systems that are central to many cellular processes. These codes should facilitate the interplay between molecular simulations and cellular studies, which have been hampered by the large mismatch between the time and length scales accessible to molecular simulations and those relevant to cell biology.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Epidermal growth factor receptor cascade prioritizes the maximization of signal transduction

Many studies have been performed to quantify cell signaling. Cell signaling molecules are phosphorylated in response to extracellular stimuli, with the phosphorylation sequence forming a signal cascade. The information gain during a signal event is given by the logarithm of the phosphorylation molecule ratio. The average information gain can be regarded as the signal transduction quantity (ST), which is identical to the Kullback–Leibler divergence (KLD), a relative entropy. We previously reported that if the total ST value in a given signal cascade is maximized, the ST rate (STR) of each signaling molecule per signal duration (min) approaches a constant value. To experimentally verify this theoretical conclusion, we measured the STR of the epidermal growth factor (EGF)-related cascade in A431 skin cancer cells following stimulation with EGF using antibody microarrays against phosphorylated signal molecules. The results were consistent with those from the theoretical analysis. Thus, signaling transduction systems may adopt a strategy that prioritizes the maximization of ST. Furthermore, signal molecules with similar STRs may form a signal cascade. In conclusion, ST and STR are promising properties for quantitative analysis of signal transduction.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Structure of active methyl-CoM reductase, Earth’s main methane producer

Summary Our work reveals the structure of the active state of Methyl-Coenzyme M Reductase (MCR), the key and rate-limiting enzyme in biological methane formation. We find large differences between the active Ni(I) and inactive Ni(II) proteins and provide insight into how nature makes and breaks the C-H bond of methane. The Ni(II)-F430 center in inactive MCR contains four planar nitrogen ligands, a lower axial glutamine oxo, and an upper axial thiolate. The Ni(I)-enzyme replaces the axial ligands with a single water. The one-electron redox change results in movement of the Ni ion and upward swing of the β-lactam ring in the tetrapyrrole coupled to a domino-like protein quake through second sphere residues, inter-subunit interactions, a substrate tunnel, affecting even the dimensions of the unit cell. These structural changes lead Ni(I)-MCR to release a charge clamp that, in the Ni(II) state, locks down substrate Coenzyme B. Determining the Ni(I)-MCR structure required development of rigorous anaerobic crystallographic techniques. Validation of the MCR redox state was accomplished by in-line and parallel spectroscopic and unit cell analyses. This structure has large implications for developing technologies to limit methane emissions and efficiently produce biofuels. Methodology described here will enhance structural biology for other oxygen-sensitive enzymes.

09 BIOMASS FUELS↗
DOE OSTI2024

Ligand‐Mediated Quantum Yield Enhancement in 1‐D Silver Organothiolate Metal–Organic Chalcogenolates

X-ray free electron laser (XFEL) microcrystallography and synchrotron single-crystal crystallography are used to evaluate the role of organic substituent position on the optoelectronic properties of metal–organic chalcogenolates (MOChas). MOChas are crystalline 1D and 2D semiconducting hybrid materials that have varying optoelectronic properties depending on composition, topology, and structure. While MOChas have attracted much interest, small crystal sizes impede routine crystal structure determination. A series of constitutional isomers where the aryl thiol is functionalized by either methoxy or methyl ester are solved by small molecule serial femtosecond X-ray crystallography (smSFX) and single crystal rotational crystallography. While all the methoxy examples have a low quantum yield (0-1%), the methyl ester in the ortho position yields a high quantum yield of 22%. Here, the proximity of the oxygen atoms to the silver inorganic core correlates to a considerable enhancement of quantum yield. Four crystal structures are solved at a resolution range of 0.8–1.0 Å revealing a collapse of the 2D topology for functional groups in the 2- and 3- positions, resulting in needle-like crystals. Further analysis using density functional theory (DFT) and many-body perturbation theory (MBPT) enables the exploration of complex excitonic phenomena within easily prepared material systems.

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

Exploring the impact of nucleotide length on lipid nanoparticle structure and properties

Lipid nanoparticles (LNPs) are versatile carriers for nucleic acid (NA) therapeutics, including ASOs, siRNA, mRNA, and poly-IC. While lipid composition is known to influence LNP properties, the impact of NA length on morphology and internal structure is less understood, particularly during the stages of carrier–cargo assembly. Here, we examine NA length and lipid composition immediately after mixing using high-throughput SAXS, dynamic light scattering, and cryogenic electron microscopy. All LNPs form ordered NA/lipid compartments, with longer NAs promoting inverse hexagonal (H II ) phases and larger intercompartment distances. In contrast, short NAs, especially in formulations with SM102 ionizable lipid, favor lamellar phases. SAXS peak deconvolution quantifies ordered versus disordered phases via a Robustness of Ordered Phase factor, which correlates with particle size and encapsulation efficiency. Formulations with MC3 ionizable and DOPE helper lipids exhibit the most stable H II -phase packing, highlighting the role of helper-lipid curvature in compartment stabilization. Variations in NA compartmentalization indicate differences in payload capacity, offering a framework for rational LNP design across diverse nucleic acid cargos.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Solvent-Mediated Control of Nanocellulose Dispersion: An Integrated Computational and Experimental Investigation

Fibrillated cellulose derived from forestry feedstocks represents a renewable and high-strength materials platform for circular bioeconomies. However, its practical implementation is hindered by the irreversible aggregation of nanocellulose architectures, including cellulose nanofibers (CNFs). Solvent-based dispersion offers a simple and practical route to prevent CNF aggregation. Here, in this work, we integrate classical and enhanced sampling molecular dynamics (MD) simulations with experimental suspension rheology and atomic force microscopy (AFM) to elucidate how solvent environments tune CNF–CNF interactions and dispersion stability. CNF–CNF contact free energies computed from MD simulations reveal reduced aggregation in acetone/water, γ-valerolactone (GVL)/water, and tetrahydrofuran (THF)/water and pure acetone compared with pure water, reflecting stronger CNF-solvent relative to inter-CNF interactions. Correspondingly, CNF-solvent suspensions in these solvent systems exhibit stronger inter-fibril network structures and enhanced recovery compared to water, indicating improved CNF-solvent affinity. Liquid cell AFM imaging in acetone–water mixtures and in pure acetone further confirm the presence of well-dispersed CNFs. By combining multiscale computation with targeted experiments, this study establishes a rational framework for solvent design to achieve stable nanocellulose dispersions for high-strength biobased materials and efficient bioenergy conversion.

cellulose↗
DOE OSTI2025

Naturally ornate RNA-only complexes revealed by cryo-EM

The structures of natural RNAs remain poorly characterized and may hold numerous surprises. Here we report three-dimensional structures of three large ornate bacterial RNAs using cryo-electron microscopy (cryo-EM). GOLLD (Giant, Ornate, Lake- and Lactobacillales-Derived), ROOL (Rumen-Originating, Ornate, Large) and OLE (Ornate Large Extremophilic) RNAs form homo-oligomeric complexes whose stoichiometries are retained at lower concentrations than measured in cells. OLE RNA forms a dimeric complex with long co-axial pipes spanning two monomers. Both GOLLD and ROOL form distinct RNA-only multimeric nanocages with diameters larger than the ribosome, each empty except for a disordered loop. Extensive intramolecular and intermolecular A-minor interactions, kissing loops, an unusual A–A helix and other interactions stabilize the three complexes. Sequence covariation analysis of these large RNAs reveals evolutionary conservation of intermolecular interactions, supporting the biological importance of large, ornate RNA quaternary structures that can assemble without any involvement of proteins.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Detection of non‐native species formed during fibrillization of the myocilin olfactomedin domain

Abstract Glaucoma is a group of neurodegenerative diseases that together are the leading cause of irreversible blindness worldwide. Myocilin‐associated glaucoma is an inherited form of this disease, caused by intracellular aggregation of misfolded mutant myocilin. In vitro, the myocilin C‐terminal olfactomedin domain (OLF), the relevant domain for glaucoma pathogenesis, can be driven to form amyloid‐like fibrils under mild conditions. Here we characterize a species present during in vitro fibrillization. Purified OLF was subjected to fibrillization at concentrations required for downstream electron microscopy imaging and NMR spectroscopy. Additional biophysical techniques, including analytical ultracentrifugation and X‐ray crystallography, were employed to further characterize the multicomponent mixture. Negative stain transmission electron microscopy (TEM) shows a non‐native species reminiscent of known prefibrillar oligomers from other amyloid systems, NMR indicates a minor population of partially misfolded species is present in solution, and cryo‐EM imaging shows two‐dimensional protein arrays. The predominant soluble species remaining in solution after the fibril reaction is natively folded, as evidenced by X‐ray crystallography. In summary, after incubating OLF under fibrillization‐promoting conditions, there is a heterogeneous mixture consisting of soluble folded protein, mature amyloid‐like fibrils, and partially misfolded intermediate species that at present belie additional molecular detail. The characterization of OLF fibrillar species illustrates the challenges associated with developing a comprehensive understanding of the fibrillization process for large, non‐model amyloidogenic proteins.

Scelsi, Hailee F. [School of Chemistry and Biochem↗
DOE OSTI2020

The trajectory of intrahelical lesion recognition and extrusion by the human 8-oxoguanine DNA glycosylase

Efficient search for DNA damage embedded in vast expanses of the DNA genome presents one of the greatest challenges to DNA repair enzymes. We report here crystal structures of human 8-oxoguanine (oxoG) DNA glycosylase, hOGG1, that interact with the DNA containing the damaged base oxoG and the normal base G while they are nested in the DNA helical stack. The structures reveal that hOGG1 engages the DNA using different protein- DNA contacts from those observed in the previously determined lesion recognition complex and other hOGG1-DNA complexes. By applying molecular dynamics simulations, we have determined the pathways taken by the lesion and normal bases when extruded from the DNA helix and their associated free energy profiles. These results reveal how the human oxoG DNA glycosylase hOGG1 locates the lesions inside the DNA helix and facilitates their extrusion for repair.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Additive Manufacturing with Cellulose-Based Composites: Materials, Modeling, and Applications

Recent advances in large-scale additive manufacturing (AM) with polymer-based composites have enabled efficient production of high-performance materials. Cellulose nanomaterials (CNMs) have emerged as bio-based feedstocks due to their exceptional strength and sustainability. However, challenges such as hornification and poor dispersion in polymer matrices still limit large-scale CNM–polymer composite manufacturing, requiring novel strategies. Here, this review outlines an approach starting with atomic-level simulations to link molecular composition to key parameters like bulk density, viscosity, and modulus. These simulations provide data for finite element analysis (FEA), which informs large-scale experiments and reduces the need for extensive trials. The strategy explores how atomic interactions impact the morphology, adhesion, and mechanical properties of CNM-based composites in AM processes. The review also discusses current developments in AM, along with predictions of mechanical and thermal properties for structural applications, packaging, flexible electronics, and hydrogel scaffolds. By integrating experimental findings with molecular dynamics (MD) simulations and finite element modeling (FEM), valuable insights for material design, process optimization, and performance enhancement in CNM-based AM are provided to address ongoing challenges.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Phytosulfokine downregulates defense‐related WRKY transcription factors and attenuates pathogen‐associated molecular pattern‐triggered immunity

SUMMARY Phytosulfokine (PSK) is a plant growth‐promoting peptide hormone that is perceived by its cell surface receptors PSKR1 and PSKR2 in Arabidopsis. Plants lacking the PSK receptors show phenotypes consistent with PSK signaling repressing some plant defenses. To gain further insight into the PSK signaling mechanism, comprehensive transcriptional profiling of Arabidopsis treated with PSK was performed, and the effects of PSK treatment on plant defense readouts were monitored. Our study indicates that PSK's major effect is to downregulate defense‐related genes; it has a more modest effect on the induction of growth‐related genes. WRKY transcription factors (TFs) emerged as key regulators of PSK‐responsive genes, sharing commonality with a pathogen‐associated molecular pattern (PAMP) responses, flagellin 22 (flg22), but exhibiting opposite regulatory directions. These PSK‐induced transcriptional changes were accompanied by biochemical and physiological changes that reduced PAMP responses, notably mitogen‐activated protein kinase (MPK) phosphorylation (previously implicated in WRKY activation) and the cell wall modification of callose deposition. Comparison with previous studies using other growth stimuli (the sulfated plant peptide containing sulfated tyrosine [PSY] and Pseudomonas simiae strain WCS417) also reveals WRKY TFs' overrepresentations in these pathways, suggesting a possible shared mechanism involving WRKY TFs for plant growth–defense trade‐off.

Liu, Dian [Biochemistry and Molecular Biophysics T↗
DOE OSTI2021

Ligand modulation of the conformational dynamics of the A 2A adenosine receptor revealed by single-molecule fluorescence

G protein-coupled receptors (GPCRs) are the largest class of transmembrane proteins, making them an important target for therapeutics. Activation of these receptors is modulated by orthosteric ligands, which stabilize one or several states within a complex conformational ensemble. The intra- and inter-state dynamics, however, is not well documented. Here, we used single-molecule fluorescence to measure ligand-modulated conformational dynamics of the adenosine A 2A receptor (A 2A R) on nanosecond to millisecond timescales. Experiments were performed on detergent-purified A 2R in either the ligand-free (apo) state, or when bound to an inverse, partial or full agonist ligand. Single-molecule Förster resonance energy transfer (smFRET) was performed on detergent-solubilized A 2A R to resolve active and inactive states via the separation between transmembrane (TM) helices 4 and 6. The ligand-dependent changes of the smFRET distributions are consistent with conformational selection and with inter-state exchange lifetimes ≥ 3 ms. Local conformational dynamics around residue 229 6.31 on TM6 was measured using fluorescence correlation spectroscopy (FCS), which captures dynamic quenching due to photoinduced electron transfer (PET) between a covalently-attached dye and proximal aromatic residues. Global analysis of PET-FCS data revealed fast (150–350 ns), intermediate (50–60 μs) and slow (200–300 μs) conformational dynamics in A 2A R, with lifetimes and amplitudes modulated by ligands and a G-protein mimetic (mini-G s ). Most notably, the agonist binding and the coupling to mini-G s accelerates and increases the relative contribution of the sub-microsecond phase. Molecular dynamics simulations identified three tyrosine residues (Y112, Y288 7.53 , and Y290 7.55 ) as being responsible for the dynamic quenching observed by PET-FCS and revealed associated helical motions around residue 229 6.31 on TM6. As a result, this study provides a quantitative description of conformational dynamics in A 2A R and supports the idea that ligands bias not only GPCR conformations but also the dynamics within and between distinct conformational states of the receptor.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Low cost and massively parallel force spectroscopy with fluid loading on a chip

Current approaches for single molecule force spectroscopy are typically constrained by low throughput and high instrumentation cost. Herein, a low-cost, high throughput technique is demonstrated using microfluidics for multiplexed mechanical manipulation of up to ~4000 individual molecules via molecular fluid loading on-a-chip (FLO-Chip). The FLO-Chip consists of serially connected microchannels with varying width, allowing for simultaneous testing at multiple loading rates. Molecular force measurements are demonstrated by dissociating Biotin-Streptavidin and Digoxigenin-AntiDigoxigenin interactions along with unzipping of double stranded DNA of varying sequence under different dynamic loading rates and solution conditions. Rupture force results under varying loading rates and solution conditions are in good agreement with prior studies, verifying a versatile approach for single molecule biophysics and molecular mechanobiology. FLO-Chip enables straightforward, rapid, low-cost, and portable mechanical testing of single molecules that can be implemented on a wide range of microscopes to broaden access and may enable new applications of molecular force spectroscopy.

42 ENGINEERING↗
DOE OSTI2023

Structural basis for severe pain caused by mutations in the voltage sensors of sodium channel NaV1.7

Voltage-gated sodium channels in peripheral nerves conduct nociceptive signals from nerve endings to the spinal cord. Mutations in voltage-gated sodium channel NaV1.7 are responsible for a number of severe inherited pain syndromes, including inherited erythromelalgia (IEM). Here, we describe the negative shifts in the voltage dependence of activation in the bacterial sodium channel NaVAb as a result of the incorporation of four different IEM mutations in the voltage sensor, which recapitulate the gain-of-function effects observed with these mutations in human NaV1.7. Crystal structures of NaVAb with these IEM mutations revealed that a mutation in the S1 segment of the voltage sensor facilitated the outward movement of S4 gating charges by widening the pathway for gating charge translocation. In contrast, mutations in the S4 segments modified hydrophobic interactions with surrounding amino acid side chains or membrane phospholipids that would enhance the outward movement of the gating charges. These results provide key structural insights into the mechanisms by which these IEM mutations in the voltage sensors can facilitate outward movements of the gating charges in the S4 segment and cause hyperexcitability and severe pain in IEM. Our work gives new insights into IEM pathogenesis at the near-atomic level and provides a molecular model for mutation-specific therapy of this debilitating disease.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Exploration of Nirmatrelvir Derivatives as Optimized SARS‐CoV‐2 Antivirals

Nirmatrelvir (NMV) is a SARS‐CoV‐2 antiviral component of the approved COVID‐19 therapeutic Paxlovid. It is a reversible covalent inhibitor of SARS‐CoV‐2 main protease (M Pro ) that is effluxed from human cells by P‐glycoprotein (P‐gp). To identify NMV analogs with improved potency and reduced P‐gp efflux, a structure–activity relationship campaign was conducted. Warheads alternative to nitrile for engaging the active site cysteine were tested showing aldehyde and dichloroacetamide with better enzyme inhibition potency. Crystal structure of MPI‐136−M Pro shows its aldehyde warhead forming a thiohemiacetal with active Cys145 of M Pro . Several S4 binders were explored revealing that an O‐to‐S shift at the N ‐terminal amide leads to better enzyme inhibition. By exploring different combinations of S2, S3, and S4 binders, two inhibitors with better enzyme inhibition potency than NMV were found. Crystal structure of MPI‐148, with ( S )‐2‐azaspiro[4,5]decane‐3‐carboxylate as an alternative S2 binder, shows extensive hydrogen‐bond networks for locking the inhibitor in active site, explaining high affinity of NMV analogs. Further characterization of cellular M Pro engagement and antiviral potency against SARS‐CoV‐2 revealed four inhibitors with greater potency than NMV in P‐gp‐expressing cells. Studies with the P‐gp inhibitor CP‐100356 showed that these compounds were less sensitive to P‐gp inhibition than NMV, consistent with reduced P‐gp‐mediated efflux.

Alugubelli, Yugendar R. [Texas A&M Drug Discovery ↗
DOE OSTI2022

Structural features of Dnase1L3 responsible for serum antigen clearance

Autoimmunity develops when extracellular DNA released from dying cells is not cleared from serum. While serum DNA is primarily digested by Dnase1 and Dnase1L3, Dnase1 cannot rescue autoimmunity arising from Dnase1L3 deficiencies. Dnase1L3 uniquely degrades antigenic forms of cell-free DNA, including DNA complexed with lipids and proteins. The distinct activity of Dnase1L3 relies on its unique C-terminal Domain (CTD), but the mechanism is unknown. We used multiple biophysical techniques and functional assays to study the interplay between the core catalytic domain and the CTD. While the core domain resembles Dnase1, there are key structural differences between the two enzymes. First, Dnase1L3 is not inhibited by actin due to multiple differences in the actin recognition site. Second, the CTD augments the ability of the core to bind DNA, thereby facilitating the degradation of complexed DNA. Together, these structural insights will inform the development of Dnase1L3-based therapies for autoimmunity.

59 BASIC BIOLOGICAL SCIENCES↗