Search NASA⌕ Search

SEARCH · Search NASA

Results for “Protein modeling”

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

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

At least 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Light-induced electron spin qubit coherences in the purple bacteria reaction center protein

Photosynthetic reaction center proteins (RCs) provide ideal model systems for studying quantum entanglement between multiple spins, a quantum mechanical phenomenon wherein the properties of the entangled particles become inherently correlated. Following light-generated sequential electron transfer, RCs generate spin-correlated radical pairs (SCRPs), also referred to as entangled spin qubit (radical) pairs (SQPs). Understanding and controlling coherence mechanisms in SCRP/SQPs is important for realizing practical uses of electron spin qubits in quantum sensing applications. The bacterial RC (bRC) provides an experimental system for exploring quantum effects in the SCRP P 865 + Q A − , where P 865 , a special pair of bacteriochlorophylls, is the primary donor, and Q A is the primary quinone acceptor. In this study, we focus on understanding how local molecular environments and isotopic substitution, particularly deuteration, influence spin coherence times (T M ). Using high-frequency electron paramagnetic resonance (EPR) spectroscopy, we observed that the local environment surrounding P 865 and Q A plays a significant role in determining T M . Our findings show that while deuteration led to a modest increase in T M , particularly at low temperatures, but the effect was substantially smaller than predicted by classical nuclear spin diffusion alone. This result is in contrast to our previous study of the photosystem I (PSI) RC, where no increase in T M was observed upon deuteration. Theoretical modeling identified several methyl groups at key distances from the spin centers of both bRC and PSI, and methyl group tunneling at low temperatures has been previously suggested as a mechanism for enhanced spin decoherence. Additionally, our study revealed a strong dependence of spin coherence on the orientation of the external magnetic field, highlighting the influence of the protein microenvironment on spin dynamics. In conclusion, these results offer new insights for optimizing coherence times in quantum system design for quantum information science and sensing applications.

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

Subcutaneous administration of insulin-like growth factor (IGF)-II/IGF binding protein-2 complex stimulates bone formation and prevents loss of bone mineral density in a rat model of disuse osteoporosis

Elevated serum levels of insulin-like growth factor binding protein-2 (IGFBP-2) and a precursor form of IGF-II are associated with marked increases in bone formation and skeletal mass in patients with hepatitis C-associated osteosclerosis. In vitro studies indicate that IGF-II in complex with IGFBP-2 has high affinity for bone matrix and is able to stimulate osteoblast proliferation. The purpose of this study was to determine the ability of the IGF-II/IGFBP-2 complex to increase bone mass in vivo. Osteopenia of the femur was induced by unilateral sciatic neurectomy in rats. At the time of surgery, 14-day osmotic minipumps containing vehicle or 2 microg IGF-II+9 microg IGFBP-2/100g body weight/day were implanted subcutaneously in the neck. Bone mineral density (BMD) measurements were taken the day of surgery and 14 days later using a PIXImus small animal densitometer. Neurectomy of the right hindlimb resulted in a 9% decrease in right femur BMD (P<0.05 vs. baseline). This loss in BMD was completely prevented by treatment with IGF-II/IGFBP-2. On the control limb, there was no loss of BMD over the 14 days and IGF-II/IGFBP-2 treatment resulted in a 9% increase in left femur BMD (P<0.05). Bone histomorphometry indicated increases in endocortical and cancellous bone formation rates and in trabecular thickness. These results demonstrate that short-term administration of the IGF-II/IGFBP-2 complex can prevent loss of BMD associated with disuse osteoporosis and stimulate bone formation in adult rats. Furthermore, they provide proof of concept for a novel anabolic approach to increasing bone mass in humans with osteoporosis.

NASA Discipline Regulatory Physiology↗
NASA NTRS1992

Pathways of proton release in the bacteriorhodopsin photocycle

The pH dependencies of the rate constants in the photocycles of recombinant D96N and D115N/D96N bacteriorhodopsins were determined from time-resolved difference spectra between 70 ns and 420 ms after photoexcitation. The results were consistent with the model suggested earlier for proteins containing D96N substitution: BR hv----K----L----M1----M2----BR. Only the M2----M1 back-reaction was pH-dependent: its rate increased with increasing [H+] between pH 5 and 8. We conclude from quantitative analysis of this pH dependency that its reverse, the M1----M2 reaction, is linked to the release of a proton from a group with a pKa = 5.8. This suggests a model for wild-type bacteriorhodopsin in which at pH greater than 5.8 the transported proton is released on the extracellular side from this as yet unknown group and on the 100-microseconds time scale, but at pH less than 5.8, the proton release occurs from another residue and later in the photocycle most likely directly from D85 during the O----BR reaction. We postulate, on the other hand, that proton uptake on the cytoplasmic side will be by D96 and during the N----O reaction regardless of pH. The proton kinetics as measured with indicator dyes confirmed the unique prediction of this model: at pH greater than 6, proton release preceded proton uptake, but at pH less than 6, the release was delayed until after the uptake. The results indicated further that the overall M1----M2 reaction includes a second kinetic step in addition to proton release; this is probably the earlier postulated extracellular-to-cytoplasmic reorientation switch in the proton pump.

NASA Discipline Exobiology↗
DOE OSTI2025

A compendium of human gene functions derived from evolutionary modelling

A comprehensive, computable representation of the functional repertoire of all macromolecules encoded within the human genome is a foundational resource for biology and biomedical research. The Gene Ontology Consortium has been working towards this goal by generating a structured body of information about gene functions, which now includes experimental findings reported in more than 175,000 publications for human genes and genes in experimentally tractable model organisms 1,2 . Here, we describe the results of a large, international effort to integrate all of these findings to create a representation of human gene functions that is as complete and accurate as possible. Specifically, we apply an expert-curated, explicit evolutionary modelling approach to all human protein-coding genes. This approach integrates available experimental information across families of related genes into models that reconstruct the gain and loss of functional characteristics over evolutionary time. The models and the resulting set of 68,667 integrated gene functions cover approximately 82% of human protein-coding genes. The functional repertoire reveals a marked preponderance of molecular regulatory functions, and the models provide insights into the evolutionary origins of human gene functions. We show that our set of descriptions of functions can improve the widely used genomic technique of Gene Ontology enrichment analysis. The experimental evidence for each functional characteristic is recorded, thereby enabling the scientific community to help review and improve the resource, which we have made publicly available.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2004

Radiation Risk Assessment of the Individual Astronaut: A Complement to Radiation Interests at the NIH

Predicting human risks following exposure to space radiation is uncertain in part because of unpredictable distribution of high-LET and low-dose-derived damage amongst cells in tissues, unknown synergistic effects of microgravity upon gene- and protein-expression, and inadequately modeled processing of radiation-induced damage within cells to produce rare and late-appearing malignant cancers. Furthermore, estimation of risks of radiogenic outcome within small numbers of astronauts is not possible using classic epidemiologic study. It therefore seems useful to develop strategies of risk-assessment based upon large datasets acquired from correlated biological models useful for resolving radiogenic risk-assessment for irradiated individuals. In this regard, it is suggested that sensitive cellular biodosimeters that simultaneously report 1) the quantity of absorbed dose after exposure to ionizing radiation, 2) the quality of radiation delivering that dose, and 3) the biomolecular risk of malignant transformation be developed in order to resolve these NASA-specific challenges. Multiparametric cellular biodosimeters could be developed using analyses of gene-expression and protein-expression whereby large datasets of cellular response to radiation-induced damage are analyzed for markers predictive for acute response as well as cancer-risk. A new paradigm is accordingly addressed wherein genomic and proteomic datasets are registered and interrogated in order to provide statistically significant dose-dependent risk estimation in individual astronauts. This evaluation of the individual for assessment of radiogenic outcomes connects to NIH program in that such a paradigm also supports assignment of a given patient to a specific therapy, the diagnosis of response of that patient to therapy, and the prediction of risks accumulated by that patient during therapy - such as risks incurred by scatter and neutrons produced during high-energy Intensity-Modulated Radiation Therapy. Value of assessment of radiogenic outcome for individuals exposed to radiation is suggested to be common to both NASA and NIH.

Richmond, Robert C.↗
DOE OSTI2025

Plant sulfate transporter protein sequences for phylogenetic analysis

Sulfur is an essential macronutrient that supports plant growth, development, and responses to environmental stress. Sulfate is the predominant inorganic form of sulfur in soils, and its uptake by roots and translocation to shoots are facilitated by the sulfate transporter (SULTR) family of proteins. Although the first plant SULTR gene was identified nearly three decades ago, several subfamily members, particularly those in the expansive and angiosperm-specific SULTR3 group, remain poorly characterized. To support comprehensive phylogenetic and sequence-based analyses, we compiled a curated dataset of 262 SULTR protein sequences from 22 plant species spanning the evolutionary breadth of land plants. This collection includes representatives from two basal lineages, two early-divergent angiosperms, six monocots, and ten dicots. All sequences were extracted from genome assemblies available in Phytozome v13 (Joint Genome Institute) and manually curated, with cross-referencing to additional databases such as NCBI when needed. This dataset provides a valuable resource for reconstructing the evolutionary history of the SULTR family, with particular emphasis on the diversification of SULTR3 transporters in flowering plants. This resource may also support functional annotation, comparative genomics, and structural modeling of sulfate transport proteins.

CBI↗
DOE OSTI2025

Machine learning approaches for integrating multi-omics data to expand microbiome annotation (Final Technical Report)

We fulfilled all original three aims of the proposal. Following the earlier release (during the first phase of the project at Montana) of software that identifies and fills gaps in the annotation of metabolic proteins within bacterial genomes, we have nearly completed a second gap-filling tool that improves accuracy and explainability. We completed software for alignment-based annotation of protein coding DNA, allowing for coding frameshifts caused by sequencing error. Finally, we completed a neural embedding model for identifying similarities between protein sequences based on amino-wise latent vectors.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1991

Macromolecular recognition: Structural aspects of the origin of the genetic system

Theoretical simulation of prebiotic chemical processes is an invaluable tool for probing the phenomenon of the evolution of life. Using computational and modeling techniques and guided by analogies from present day systems, we seek to understand the emergence of the genetic apparatus, enzymatic catalysis and protein synthesis under prebiotic conditions. Modeling of the ancestral aminoacyl-tRNA-synthetases (aRS) may provide important clues to the emergence of the genetic code and the protein synthetic machinery. The minimal structural requirements for the catalysis of tRNA aminoacylation are being explored. A formation of an aminoacyl adenylate was studied in the framework of ab initio molecular orbital theory. The role of individual residues in the vicinity of the TyrRS active site was examined, and the effect of all possible amino acids substitutions near the active site was examined. A formation of aminoacyl tRNA was studied by the molecular modeling system SYBYL with the high resolution crystallographic structures of the present day tRNA, aRS's complexes. The ultimate goal is to propose a simple RNA segment that is small enough to be build in the primordial chemical environment but maintains the specificity and catalytic activity of the contemporary RNA enzyme. To understand the mechanism of ribozyme catalyzed reactions, ab initio and semi-empirical (ZINDO) programs were used to investigate the reaction path of transphosphorylation. A special emphasis was placed on the possible catalytic and structural roles played by the coordinated magnesium cation. Both the inline and adjacent mechanisms of transphosphorylation were studied. The structural characteristics of the target helices, particularly a possible role for the G-T pair, is also studied by a molecular dynamics (MD) simulation technique.

Rein, Robert↗
NASA NTRS1994

Protein structure, electron transfer and evolution of prokaryotic photosynthetic reaction centers

Photosynthetic reaction centers from a variety of organisms have been isolated and characterized. The groups of prokaryotic photosynthetic organisms include the purple bacteria, the filamentous green bacteria, the green sulfur bacteria and the heliobacteria as anoxygenic representatives as well as the cyanobacteria and prochlorophytes as oxygenic representatives. This review focuses on structural and functional comparisons of the various groups of photosynthetic reaction centers and considers possible evolutionary scenarios to explain the diversity of existing photosynthetic organisms.

NASA Discipline Exobiology↗
NASA NTRS1993

Deriving the intermediate spectra and photocycle kinetics from time-resolved difference spectra of bacteriorhodopsin. The simpler case of the recombinant D96N protein

The bacteriorhodopsin photocycle contains more than five spectrally distinct intermediates, and the complexity of their interconversions has precluded a rigorous solution of the kinetics. A representation of the photocycle of mutated D96N bacteriorhodopsin near neutral pH was given earlier (Varo, G., and J. K. Lanyi. 1991. Biochemistry. 30:5008-5015) as BRhv-->K<==>L<==>M1-->M2--> BR. Here we have reduced a set of time-resolved difference spectra for this simpler system to three base spectra, each assumed to consist of an unknown mixture of the pure K, L, and M difference spectra represented by a 3 x 3 matrix of concentration values between 0 and 1. After generating all allowed sets of spectra for K, L, and M (i.e., M1 + M2) at a 1:50 resolution of the matrix elements, invalid solutions were eliminated progressively in a search based on what is expected, empirically and from the theory of polyene excited states, for rhodopsin spectra. Significantly, the average matrix values changed little after the first and simplest of the search criteria that disallowed negative absorptions and more than one maximum for the M intermediate. We conclude from the statistics that during the search the solutions strongly converged into a narrow region of the multidimensional space of the concentration matrix. The data at three temperatures between 5 and 25 degrees C yielded a single set of spectra for K, L, and M; their fits are consistent with the earlier derived photocycle model for the D96N protein.

NASA Program Exobiology↗
DOE OSTI2025

Consequences of the failure of equipartition for the p – V behavior of liquid water and the hydration free energy components of a small protein

Earlier we showed that in the molecular dynamics simulation of a rigid model of water it is necessary to use an integration time-step δ t ≤ 0.5 fs to ensure equipartition between translational and rotational modes. Here we extend that study in the NVT ensemble to NpT conditions and to an aqueous protein. We study neat liquid water with the rigid, SPC/E model and the protein BBA (PDB ID: 1FME) solvated in the rigid, TIP3P model. We examine integration time-steps ranging from 0.5 fs to 4.0 fs for various thermostat plus barostat combinations. We find that a small δ t is necessary to ensure consistent prediction of the simulation volume. Hydrogen mass repartitioning alleviates the problem somewhat, but is ineffective for the typical time-step used with this approach. The compressibility, a measure of volume fluctuations, and the dielectric constant, a measure of dipole moment fluctuations, are also seen to be sensitive to δ t . Using the mean volume estimated from the NpT simulation, we examine the electrostatic and van der Waals contribution to the hydration free energy of the protein in the NVT ensemble. These contributions are also sensitive to δ t . In going from δ t = 2 fs to δ t = 0.5 fs, the change in the net electrostatic plus van der Waals contribution to the hydration of BBA is already in excess of the folding free energy reported for this protein.

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

Structural models of the MscL gating mechanism

Three-dimensional structural models of the mechanosensitive channel of large conductance, MscL, from the bacteria Mycobacterium tuberculosis and Escherichia coli were developed for closed, intermediate, and open conformations. The modeling began with the crystal structure of M. tuberculosis MscL, a homopentamer with two transmembrane alpha-helices, M1 and M2, per subunit. The first 12 N-terminal residues, not resolved in the crystal structure, were modeled as an amphipathic alpha-helix, called S1. A bundle of five parallel S1 helices are postulated to form a cytoplasmic gate. As membrane tension induces expansion, the tilts of M1 and M2 are postulated to increase as they move away from the axis of the pore. Substantial expansion is postulated to occur before the increased stress in the S1 to M1 linkers pulls the S1 bundle apart. During the opening transition, the S1 helices and C-terminus amphipathic alpha-helices, S3, are postulated to dock parallel to the membrane surface on the perimeter of the complex. The proposed gating mechanism reveals critical spatial relationships between the expandable transmembrane barrel formed by M1 and M2, the gate formed by S1 helices, and "strings" that link S1s to M1s. These models are consistent with numerous experimental results and modeling criteria.

NASA Discipline Cell Biology↗
DOE OSTI2025

Machine learning prediction of enzyme optimum pH

The relationship between pH and enzyme catalytic activity, especially the optimal pH (pH opt ) at which enzymes function, is critical for biotechnological applications. Hence, computational methods to predict pH opt will enhance enzyme discovery and design by facilitating accurate identification of enzymes that function optimally at specific pH levels, and by elucidating sequence-function relationships. Here, in this study, we proposed and evaluated various machine learning methods for predicting pH opt , conducting extensive hyperparameter optimization and training over 11,000 model instances. Our results demonstrate that models utilizing language model embeddings markedly outperform other methods in predicting pHopt. We present EpHod, the best-performing model, to predict pHopt, making it publicly available to researchers. From sequence data, EpHod directly learns structural and biophysical features that relate to pH opt , including proximity of residues to the catalytic centre and the accessibility of solvent molecules. Overall, EpHod presents a promising advancement in pH opt prediction and will potentially speed up the development of enzyme technologies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Dataset for manuscript "Consequences of the failure of equipartition for the p-V behavior of liquid water and the hydration free energy components of a small protein"

Previously, we showed that in the molecular dynamics simulation of a rigid model of water it is necessary to use an integration time-step dt that is less than or equal to 0.5 fs to ensure equipartition between translational and rotational modes. We extended that study in the NVT ensemble to NpT conditions and to an aqueous protein. We study neat liquid water with the rigid, SPC/E model and the protein BBA (PDB ID: 1FME) solvated in the rigid, TIP3P model. We examined integration time-steps ranging from 0.5 fs to 4.0 fs for various thermostat plus barostat combinations. We find that a small time-step, dt, is necessary to ensure consistent prediction of the simulation volume. Hydrogen mass repartitioning alleviates the problem somewhat, but is ineffective for the typical time-step used with this approach. The compressibility, a measure of volume fluctuations, is seen to be sensitive to dt. Using the mean volume estimated from the NpT simulation, we examined the electrostatic and van der Waals contribution to the hydration free energy of the protein in the NVT ensemble. These contributions are also sensitive to dt. In going from a time-step of 2 fs to a time-step of 0.5 fs, the change in the net electrostatic plus van der Waals contribution to the hydration of BBA is already in excess of the folding free energy reported for this protein. The data-set contains the simulation metadata and log files that support the claims noted above.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Engineering PHL7 for Improved Poly(Ethylene Terephthalate) Depolymerization via Rational Design and Directed Evolution

Enzymatic depolymerization of poly(ethylene terephthalate) (PET) has emerged as a promising approach for polyester recycling, and, to date, many natural and engineered PET hydrolase enzymes have been reported. For industrial use, PET hydrolases must achieve high depolymerization extent and exhibit excellent thermostability. Here, we engineered a natural PET hydrolase, Polyester Hydrolase Leipzig #7 (PHL7), through rational design and directed evolution using a high-throughput screening platform. Four new enzymes were engineered with enhanced properties compared with the parent enzyme, wild-type PHL7 (PHL7-WT), and other benchmark PET hydrolases, under the tested conditions. In bioreactors, the exemplary engineered enzyme, PHL7-Jemez, exhibited improved ability to depolymerize amorphous PET film compared with PHL7-WT at 2.9% and 20% substrate loadings, with 37% and 270% higher hydrolysis, respectively, after 48 h. This study develops several state-of-the-art PET hydrolases and demonstrates a directed evolution platform to engineer high-performance enzymes, which can accelerate enzyme discovery toward improved biocatalytic recycling.

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

Structural Insights into the Mechanism of a Polyketide Synthase Thiocysteine Lyase Domain

Polyketide synthases (PKSs) are renowned for the structural diversity of the polyketide natural products they produce, but sulfur-containing functionalities are rarely installed by PKSs. We previously characterized thiocysteine lyase (SH) domains involved in the biosynthesis of the leinamycin (LNM) family of natural products, exemplified by LnmJ-SH and guangnanmycin (GnmT-SH). Here we report a detailed investigation into the PLP-dependent reaction catalyzed by the SH domains, guided by a 1.8 Å resolution crystal structure of GnmT-SH. A series of elaborate substrate mimics were synthesized to answer specific questions garnered from the crystal structure and from the biosynthetic logic of the LNM family of natural products. Here, through a combination of bioinformatics, molecular modeling, in vitro assays, and mutagenesis, we have developed a detailed model of acyl carrier protein (ACP)-tethered substrate-SH, and interdomain interactions, that contribute to the observed substrate specificity. Comparison of the GnmT-SH structure with archetypical PLP-dependent enzyme structures revealed how Nature, via evolution, has modified a common protein structural motif to accommodate an ACP-tethered substrate, which is significantly larger than any of those previously characterized. Overall, this study demonstrates how PLP-dependent chemistry can be incorporated into the context of PKS assembly lines and sets the stage for engineering PKSs to produce sulfur-containing polyketides.

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

Electrochemical cofactor recycling of bacterial microcompartments

Bacterial microcompartments (BMCs) are prokaryotic organelles that consist of a protein shell which sequesters metabolic reactions in its interior. While most of the substrates and products are relatively small and can permeate the shell, many of the encapsulated enzymes require cofactors that must be regenerated inside. We have analyzed the occurrence of an enzyme previously assigned as a cobalamin (vitamin B 12 ) reductase and, curiously, found it in many unrelated BMC types that do not employ B 12 cofactors. We propose Nicotinamide adenine dinucleotide (NAD+) regeneration as the function of this enzyme and name it Metabolosome Nicotinamide Adenine Dinucleotide Hydrogen (NADH) dehydrogenase (MNdh). Its partner shell protein BMC-T SE (tandem domain BMC shell protein of the single layer type for electron transfer) assists in passing the generated electrons to the outside. We support this hypothesis with bioinformatic analysis, functional assays, Electron Paramagnetic Resonance spectroscopy, protein voltammetry, and structural modeling verified with X-ray footprinting. This finding represents a paradigm for the BMC field, identifying a new, widely occurring route for cofactor recycling and a new function for the shell as separating redox environments.

bacterial microcompartment↗
DOE OSTI2024

Heterogeneous flexibility can contribute to chromatin segregation in the cell nucleus

The highly and slightly condensed forms of chromatin, heterochromatin and euchromatin, respectively, segregate in the cell nucleus. Heterochromatin is more abundant in the nucleus periphery. Here, in this work, we study the mechanism of heterochromatin segregation by modeling interphase chromosomes as diblock ring copolymers confined in a rigid spherical shell using molecular dynamics simulations. In our model, heterochromatin and euchromatin are distinguished by their bending stiffnesses only, while an interaction potential between the spherical shell and chromatin is used to model lamin-associated proteins. Our simulations indicate that in the absence of attractive interactions between the nuclear shell and the chromatin, most heterochromatin segregates towards the nuclear interior due to the depletion of less flexible heterochromatin segments from the nuclear periphery. This inverted chromatin distribution, which is opposite to the conventional case with heterochromatin dominating at the periphery, is in accord with experimental observations in rod cells. This “inversion” is also found to be independent of the heterochromatin concentration and chromosome number. The chromatin distribution at the periphery found in vivo can be recovered by further increasing the bending stiffness of heterochromatin segments or by turning on attractive interactions between the nuclear shell and heterochromatin. Our results indicate that the bending stiffness of chromatin could be a contributor to chromosome organization along with differential effects of HP ⁢1⁢𝛼-driven phase segregation and of loop extruders and interactions with the nuclear envelope and topological constraints.

biomolecular & subcellular processes↗