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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A hybrid surrogate modeling framework for the Digital Twin of a Fluoride-salt-cooled High-temperature Reactor (FHR)

While nuclear energy is a non-greenhouse-gas emitting energy source, expensive operational costs due to the high-level of safety requirements decreases their competitiveness in the sustainable energy market. Advanced reactor concepts paired with Digital Twins aim to increase the commercialization gains of nuclear energy by reducing operational costs, increasing reactor reliability and enhancing power generation. To support Digital Twin tasks such as real-time autonomous control, proactive maintenance monitoring or optimizing power demand operations, a fast and accurate virtual representation of the Nuclear Power Plant (NPP) is required. The computational cost of high-fidelity, physics-based models are unsuitable for real-time analysis or scalability. Here, in this work, a hybrid surrogate modeling framework is developed fora Fluoride-salt-cooled High-temperature Reactor (FHR) that leverages physics-inspired models for key reactor components and uses data-driven methods for rapid system state space prediction. The Xenon reactivity feedback model is integrated to inform the surrogate model about the reactor core and the homologous pump theory model is the basis for representing pump degradation. Using a detailed, two dimensional thermal hydraulics model to generate data on the FHR, we train a network of Vectorized Autoregressive Moving-Average with eXogenous input (VARMAX) models to predict the remaining state values. The result is a surrogate model that provides a detailed reactor state representation of 41 system states and a pump degradation analysis. The framework is applied to Load Follows profiles, yielding high accuracy and a speedup that is more than 4000x faster compared to the higher- fidelity thermal hydraulics model, enabling real-time operational intelligence and applications in long horizon predictions. While the surrogate model framework is demonstrated for the particular case of FHR, the hybrid physical/data-driven modeling approach including the network of surrogates and the underlying modularity has the potential to be applied to other physical asset systems.

Digital Twins↗
NASA NTRS1997

A comparative modeling of supernova 1993J

The light curve of Supernova 1993J is calculated using two computational radiation transport approaches. The two approaches are represented by the computer codes STELLA and EDDINGTON. The emphasis is on the shock breakout and the photometry in the U, B and V bands during the first 120 days of the supernova. The STELLA model includes implicit hydrodynamics and is able to model early supernova evolution before the expansion is homologous. The STELLA model employs multi-group photonics and is able to follow the radiation as it decouples from the matter. The EDDINGTON code uses an algorithm for integrating the transport equation which assumes homologous expansion and uses a finer frequency resolution. The agreement between the two codes is considered to be satisfactory only in the case where compatible physical assumptions are made concerning the opacity. The assumptions are justified. The continuum spectrum for SN 1993J is predicted near the shock breakout to be superior to that predicted by standard single energy group hydrocodes. The uncertainties involved in current time dependent models of supernova light curves are discussed.

Blinnikov, Sergei↗
DOE OSTI2024

An FPGA-based hardware accelerator supporting sensitive sequence homology filtering with profile hidden Markov models

Abstract Background Sequence alignment lies at the heart of genome sequence annotation. While the BLAST suite of alignment tools has long held an important role in alignment-based sequence database search, greater sensitivity is achieved through the use of profile hidden Markov models (pHMMs). Here, we describe an FPGA hardware accelerator, called HAVAC, that targets a key bottleneck step (SSV) in the analysis pipeline of the popular pHMM alignment tool, HMMER. Results The HAVAC kernel calculates the SSV matrix at 1739 GCUPS on a $$\sim$$ ∼ $3000 Xilinx Alveo U50 FPGA accelerator card, $$\sim$$ ∼ 227× faster than the optimized SSV implementation in nhmmer . Accounting for PCI-e data transfer data processing, HAVAC is 65× faster than nhmmer’s SSV with one thread and 35× faster than nhmmer with four threads, and uses $$\sim$$ ∼ 31% the energy of a traditional high end Intel CPU. Conclusions HAVAC demonstrates the potential offered by FPGA hardware accelerators to produce dramatic speed gains in sequence annotation and related bioinformatics applications. Because these computations are performed on a co-processor, the host CPU remains free to simultaneously compute other aspects of the analysis pipeline.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Structural determination of a full-length plant cellulose synthase informed by experimental and in silico methods

Three-dimensional structure determination and prediction of proteins with intrinsically disordered regions, unstructured regions, conformational flexibility, and lacking homologous structures are challenging. We previously predicted and refined an in silico structure of a plant cellulose synthase from cotton (GhCESA1), and more recently, cryo-electron microscopy (cryo-EM) has resolved a majority of the lengths of two CESA structures from poplar (PttCESA8) and cotton (GhCESA7). However, 26–30% of these cryo-EM structures remain unresolved, including the N-terminal domain, half of the class-specific region, the gating loop region, and the C-terminal domain. Here, we describe the generation and evaluation of a full-length hybrid GhCESA1 model based on this cryo-EM PttCESA8 structure, with unresolved regions completed using this in silico refined GhCESA1 model. All-atom molecular dynamics simulations and subsequent energy minimizations were performed for the in silico and hybrid GhCESA1 models in a lipid bilayer-water-ion environment, and structural stability, dynamics, energetics, contacts, and quality were evaluated. The unresolved regions were found to be the most dynamic, in agreement with their poor electron density with cryo-EM. The hybrid model exhibited a higher total secondary structure content, more favorable intra-protein and protein-lipid interaction energies, and improved quality metrics. Moreover, hydrogen bonding was revealed to be a primary mechanism for intra-protein and protein-lipid contacts. These results demonstrate that in silico structure prediction and refinement may be useful to augment experimental structure determination, especially for disordered and unstructured regions. Furthermore, this hybrid model can serve as a steppingstone to derive full-length homology models of other CESAs found in more experimentally tractable organisms.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Empirical evidence that glucan-interacting amino acid side chains within the transmembrane channel collectively facilitate cellulose synthase function

The fundamental mechanism of cellulose synthesis is widely conserved across Kingdoms and depends on cellulose synthases, which are processive, dual-function, family 2 glycosyltransferases (GT-2). These enzymes polymerize glucose on the cytoplasmic side of the plasma membrane and export the glucan chain to the cell surface through an integral transmembrane (TM) channel. Structural studies of active plant cellulose synthases (CESAs) have revealed interactions between the nascent glucan chain and the side chains of polar, charged, and aromatic amino acid residues that line the TM channel. However, the functional consequences of modifying these side chains have not been tested in vivo in CESAs or other processive GT-2s. To test this, we used an established in vivo assay based on genetic complementation of CESA5 in the moss, Physcomitrium patens. For accurate prediction of glucan-interacting amino acid residues, we generated a complete homotrimeric molecular model of PpCESA5 using a combination of homology and de novo modeling. All-atom molecular dynamics-based analyses of contact metrics and interaction energy identified 23 amino acid residues with high propensity to interact with the nascent glucan chain within the TM channel or on the apoplastic surface of PpCESA5. Mutating any one of 18 of these amino acid residues to alanine, thereby removing their side chains, abolished or impaired CESA function, with the strongest effects observed upon the loss of charged amino acid side chains. This provides direct evidence to support the hypothesis that multiple amino acid residues collectively maintain a smooth energy landscape within the TM channel to facilitate glucan translocation.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2007

Biochemical Kinetics Model of DSB Repair and GammaH2AX FOCI by Non-homologous End Joining

We developed a biochemical kinetics approach to describe the repair of double strand breaks (DSB) produced by low LET radiation by modeling molecular events associated with the mechanisms of non-homologous end-joining (NHEJ). A system of coupled non-linear ordinary differential equations describes the induction of DSB and activation pathways for major NHEJ components including Ku(sub 70/80), DNA-PK(sub cs), and the Ligase IV-XRCC4 hetero-dimer. The autophosphorylation of DNA-PK(sub cs and subsequent induction of gamma-H2AX foci observed after ionizing radiation exposure were modeled. A two-step model of DNA-PK(sub cs) regulation of repair was developed with the initial step allowing access of other NHEJ components to breaks, and a second step limiting access to Ligase IV-XRCC4. Our model assumes that the transition from the first to second-step depends on DSB complexity, with a much slower-rate for complex DSB. The model faithfully reproduced several experimental data sets, including DSB rejoining as measured by pulsed-field electrophoresis (PFGE), quantification of the induction of gamma-H2AX foci, and live cell imaging of the induction of Ku(sub 70/80). Predictions are made for the behaviors of NHEJ components at low doses and dose-rates, where a steady-state is found at dose-rates of 0.1 Gy/hr or lower.

Cucinotta, Francis, A.↗
DOE OSTI2023

A Wox3 -patterning module organizes planar growth in grass leaves and ligules

Grass leaves develop from a ring of primordial initial cells within the periphery of the shoot apical meristem, a pool of organogenic stem cells that generates all of the organs of the plant shoot. At maturity, the grass leaf is a flattened, strap-like organ comprising a proximal supportive sheath surrounding the stem and a distal photosynthetic blade. The sheath and blade are partitioned by a hinge-like auricle and the ligule, a fringe of epidermally derived tissue that grows from the adaxial (top) leaf surface. Together, the ligule and auricle comprise morphological novelties that are specific to grass leaves. Understanding how the planar outgrowth of grass leaves and their adjoining ligules is genetically controlled can yield insight into their evolutionary origins. Here we use single-cell RNA-sequencing analyses to identify a ‘rim’ cell type present at the margins of maize leaf primordia. Cells in the leaf rim have a distinctive identity and share transcriptional signatures with proliferating ligule cells, suggesting that a shared developmental genetic programme patterns both leaves and ligules. Moreover, we show that rim function is regulated by genetically redundant Wuschel-like homeobox3 (WOX3) transcription factors. Higher-order mutations in maize Wox3 genes greatly reduce leaf width and disrupt ligule outgrowth and patterning. Together, these findings illustrate the generalizable use of a rim domain during planar growth of maize leaves and ligules, and suggest a parsimonious model for the homology of the grass ligule as a distal extension of the leaf sheath margin.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

A multi-ancestry GWAS of Fuchs corneal dystrophy highlights the contributions of laminins, collagen, and endothelial cell regulation

Fuchs endothelial corneal dystrophy (FECD) is a leading indication for corneal transplantation, but its molecular etiology remains poorly understood. We performed genome-wide association studies (GWAS) of FECD in the Million Veteran Program followed by multi-ancestry meta-analysis with the previous largest FECD GWAS, for a total of 3970 cases and 333,794 controls. We confirm the previous four loci, and identify eight novel loci: SSBP3, THSD7A, LAMB1, PIDD1, RORA, HS3ST3B1, LAMA5, and COL18A1. We further confirm the TCF4 locus in GWAS for admixed African and Hispanic/Latino ancestries and show an enrichment of European-ancestry haplotypes at TCF4 in FECD cases. Among the novel associations are low frequency missense variants in laminin genes LAMA5 and LAMB1 which, together with previously reported LAMC1, form laminin-511 (LM511). AlphaFold 2 protein modeling, validated through homology, suggests that mutations at LAMA5 and LAMB1 may destabilize LM511 by altering inter-domain interactions or extracellular matrix binding. Finally, phenome-wide association scans and colocalization analyses suggest that the TCF4 CTG18.1 trinucleotide repeat expansion leads to dysregulation of ion transport in the corneal endothelium and has pleiotropic effects on renal function.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Topological Signatures of Adversaries in Multimodal Alignments

Topological Data Analysis for Adversarial Detection (LANL O4937) - Detects adversarial examples in vision-language models using persistent homology and two-sample testing. Combines TDA features from CLIP embeddings with statistical methods (ME, SCF, SAMMD, C2ST) for robust detection across ImageNet, CIFAR-10/100.

Bhattarai, Manish↗
DOE OSTI2024

Bacterial hemophilin homologs and their specific type eleven secretor proteins have conserved roles in heme capture and are diversifying as a family

Cellular life relies on enzymes that require metals, which must be acquired from extracellular sources. Bacteria utilize surface and secreted proteins to acquire such valuable nutrients from their environment. These include the cargo proteins of the type eleven secretion system (T11SS), which have been connected to host specificity, metal homeostasis, and nutritional immunity evasion. This Sec-dependent, Gram-negative secretion system is encoded by organisms throughout the phylum Proteobacteria, including human pathogens Neisseria meningitidis, Proteus mirabilis, Acinetobacter baumannii, and Haemophilus influenzae. Experimentally verified T11SS-dependent cargo include transferrin-binding protein B (TbpB), the hemophilin homologs heme receptor protein C (HrpC), hemophilin A (HphA), the immune evasion protein factor-H binding protein (fHbp), and the host symbiosis factor nematode intestinal localization protein C (NilC). Here, we examined the specificity of T11SS systems for their cognate cargo proteins using taxonomically distributed homolog pairs of T11SS and hemophilin cargo and explored the ligand binding ability of those hemophilin cargo homologs. In vivo expression in Escherichia coli of hemophilin homologs revealed that each is secreted in a specific manner by its cognate T11SS protein. Sequence analysis and structural modeling suggest that all hemophilin homologs share an N-terminal ligand-binding domain with the same topology as the ligand-binding domains of the Haemophilus haemolyticus heme binding protein (Hpl) and HphA. We term this signature feature of this group of proteins the hemophilin ligand-binding domain. Network analysis of hemophilin homologs revealed five subclusters and representatives from four of these showed variable heme-binding activities, which, combined with sequence-structure variation, suggests that hemophilins are diversifying in function.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2004

Intricate Crystal Structure of Dihydrolipoamide Dehydrogenase (E3) with its Binding Protein: Multiple Copies, Dynamic and Static Disorders

Human E3 and binding protein E3BP are two components of the pyruvate dehydrogenase complex. Crystallization of E3 with 221-amino acid fragment of E3BP (E3BPdd) led to crystals that diffracted to a resolution of 2.6 Angstroms. Structure determination involved molecular replacement using a dimer of E3 homolog as a search model and de novo building of the E3BPdd peptide. Solution was achieved by inclusion of one E3 dimer at a time, followed by refinement until five E3 dimers were located. This complete content of E3 provided electron density maps suitable for tracing nine peptide chains of E3BPdd, eight of them being identified with partial occupancies. Final content of the asymmetric unit consists of five E3 dimers, each binding one E3BPdd molecule. In four of these molecular complexes, E3BPdd is in static disorder resulting in E3BPdd binding to either one or the other monomer of the E3 dimer. However, E3BPdd of the fifth E3 dimer forms specific contacts that lock it at one monomer. In addition to this static disorder, E3BPdd reveals high mobility in the limited space of the crystal lattice. Support from NIH and NASA.

Makal, A.↗
DOE OSTI2026

Quantitative SANS and multi-model analysis of spacer-dependent micellization of urea-based gemini surfactants

The micellization behavior of urea-based cationic gemini surfactants was investigated using small-angle neutron scattering (SANS) with multi-model form factor analysis. A homologous series of surfactants with urea group included in the hydrophobic tail and polymethylene spacers consisting of two to ten methylene units was analyzed using three form factor models: a core–shell ellipsoid and two variants of homogeneous ellipsoids. The results from all models show a consistent trend of the micelle structures, confirming that the spacer length critically influences micellar geometry, aggregation number, and hydration. The surfactant with four CH 2 groups in the spacer formed the largest micelles with the highest aggregation number, while longer spacers led to progressively smaller, more compact aggregates. The shell hydration—quantified as the volume fraction of heavy water within the hydrophilic region—decreased systematically with increasing spacer length due to enhanced hydrophobicity of the headgroup-spacer region. Intermicellar interactions, modeled as screened Coulomb interaction using the rescaled mean spherical approximation (RMSA), revealed the strongest electrostatic repulsion for the case of four methylene groups in the spacer, corresponding to the highest micellar charge and largest interparticle spacing. The observed spacer-dependent trends were robust across all modeling approaches, demonstrating that the spacer length serves as a key structural determinant of self-assembly in this type of urea-based gemini systems. These findings provide insight into the design of gemini surfactants with tailored aggregation behavior for applications in drug delivery, nanostructure templating, and solubilization technologies.

Core–shell ellipsoid model↗
NASA NTRS1992

The structure of horizontal-branch models. I - The zero-age horizontal branch

A detailed study of the structure of zero-age horizontal-branch (ZAHB) models is studied in order to show how the hydrostatic structure of these models changes with the input parameters and determines the H-R diagram location of a given model. The properties of composite polytropes on the homology-invariant (U,V)-plane are demonstrated. A variety of test models and sequences were constructed to elucidate the underlying factors that give rise to the wide variation in HB model properties with composition. The roles of the CNO elements as nuclear catalysts and of the envelope sources, as well as the envelope helium abundance are reexamined. It is found that, for stars of a fixed range of mass arriving on the HB, the stellar distribution is determined mainly by CNO for low metallicities (Fe/H of less than about -1), but mainly by opacity sources for high metallicities. The value of Fe/H where CNO ceases to dominate depends significantly on the adopted opacity and will decrease if and when opacity estimates are revised upward.

Dorman, Ben↗
DOE OSTI2024

Harnessing Machine Learning and Data Fusion for Accurate Undocumented Well Identification in Satellite Images

This study utilizes satellite data to detect undocumented oil and gas wells, which pose significant environmental concerns, including greenhouse gas emissions. Three key findings emerge from the study. Firstly, the problem of imbalanced data is addressed by recommending oversampling techniques like Rotation–GaussianBlur–Solarization data augmentation (RGS), the Synthetic Minority Over-Sampling Technique (SMOTE), or ADASYN (an extension of SMOTE) over undersampling techniques. The performance of borderline SMOTE is less effective than that of the rest of the oversampling techniques, as its performance relies heavily on the quality and distribution of data near the decision boundary. Secondly, incorporating pre-trained models trained on large-scale datasets enhances the models’ generalization ability, with models trained on one county’s dataset demonstrating high overall accuracy, recall, and F1 scores that can be extended to other areas. This transferability of models allows for wider application. Lastly, including persistent homology (PH) as an additional input improves performance for in-distribution testing but may affect the model’s generalization for out-of-distribution testing. A careful consideration of PH’s impact on overall performance and generalizability is recommended. Overall, this study provides a robust approach to identifying undocumented oil and gas wells, contributing to the acceleration of a net-zero economy and supporting environmental sustainability efforts.

SMOTE↗
DOE OSTI2025

Modeling Protein–Protein and Protein–Ligand Interactions by the ClusPro Team in CASP16

ABSTRACT In the CASP16 experiment, our team employed hybrid computational strategies to predict both protein–protein and protein–ligand complex structures. For protein–protein docking, we combined physics‐based sampling—using ClusPro FFT docking and molecular dynamics—with AlphaFold (AF)‐based sampling, followed by AF‐based refinement. Our method produced numerous high‐accuracy complex models, including cases where AF alone failed, underscoring the critical role of physics‐based sampling alongside deep learning‐based refinement. For protein–ligand docking, we integrated the ClusPro LigTBM template‐based approach with a machine learning‐based confidence model for rescoring. The method preserves conserved interaction fragments derived from homologous complexes, followed by local resampling using physics‐based sampling and a diffusion model. Our template‐based strategy achieved a mean lDDT‐PLI of 0.69 across 233 targets, which was highly competitive. These results demonstrate that combining physics‐based modeling with AI‐driven refinement can significantly enhance the accuracy of both protein–protein and protein–ligand structure predictions.

Ashizawa, Ryota [Department of Applied Mathematics↗
NASA NTRS2001

EIT Crinkles as Evidence for the Breakout Model of Solar Eruptions

We present observations of two homologous flares in NOAA Active Region 8210 occurring on 1998 May 1 and 2, using EUV data from the EUV Imaging Telescope (EIT) on board the Solar and Heliospheric Observatory, high-resolution and high-time cadence images from the soft X-ray telescope on Yohkoh, images or fluxes from the hard X-ray telescope on Yohkoh and the BATSE experiment on board the Compton Gamma Ray Observatory, and Ca(XIX) soft X-ray spectra from the Bragg crystal spectrometer (BCS) on Yohkoh. Magnetograms indicate that the flares occurred in a complex magnetic topology, consisting of an emerging flux region (EFR) sandwiched between a sunspot to the west and a coronal hole to the east. In an earlier study we found that in EIT images, both flaring episodes showed the formation of a crinkle-like pattern of emission (EIT crinkles) occurring in the coronal hole vicinity, well away from a central 'core field' area near the EFR-sunspot boundary. With our expanded data set, here we find that most of the energetic activity occurs in the core region in both events, with some portions of the core brightening shortly after the onset of the EIT crinkles, and other regions of the core brightening several minutes later, coincident with a burst of hard X-rays; there are no obvious core brightenings prior to the onset of the EIT crinkles. These timings are consistent with the 'breakout model' of solar eruptions, whereby the emerging flux is initially constrained by a system of overlying magnetic field lines, and is able to erupt only after an opening develops in the overlying fields as a consequence of magnetic reconnection at a magnetic null point. In our case, the EIT crinkles would be a signature of this pre-impulsive phase magnetic reconnection, and brightening of the core only occurs after the core fields begin to escape through the newly created opening in the overlying fields. Morphology in soft X-ray images and properties in hard X-rays differ between the two events, with complexities that preclude a simple determination of the dynamics in the core at the times of eruption. From the BCS spectra, we find that the core region expends energy at a rate of approximately 10(exp 26) ergs/s during the time of the growth of the EIT crinkles; this rate is an upper limit to energy expended in the reconnections opening the overlying fields. Energy losses occur at an order of magnitude higher rate near the time of the peak of the events. There is little evidence of asymmetry in the spectra, consistent with the majority of the mass flows occurring normal to the line of sight. Both events have similar electron temperature dependencies on time.

Sterling, Alphonse C.↗
NASA NTRS2001

EIT Crinkles as Evidence for the Breakout Model of Solar Eruptions

We present observations of two homologous flares in NOAA active region 8210 occurring on 1998 May 1 and May 2, using EUV data from the Extreme Ultraviolet Radiation Imaging Telescope (EIT) on the Solar and Heliospheric Observatory (SOHO), high-resolution and high-time cadence images from the soft X-ray telescope (SXT) on Yohkoh, images or fluxes from the hard X-ray telescope (HXT) on Yohkoh and the BATSE experiment on the Compton Gamma Ray Observatory (CGRO), and Ca xix soft X-ray spectra from the Bragg crystal spectrometer (BCS) on Yohkoh. Magnetograms indicate that the flares occurred in a complex magnetic topology, consisting of an emerging flux region (EFR) sandwiched between a sunspot to the west and a coronal hole to the east. In an earlier study we found that in EIT images, both flaring episodes showed the formation of a crinkle-like pattern of emission ("EIT crinkles") occurring in the coronal hole vicinity, well away from a central "core field" area near the EFR-sunspot boundary. With our expanded data set, here we find that most of the energetic activity occurs in the core region in both events, with some portions of the core brightening shortly after the onset of the EIT crinkles, and other regions of the core brightening several minutes later, coincident with a burst of hard X-rays: there are no obvious core brightenings prior to the onset of the EIT crinkles. These timings are consistent with the "breakout model" of solar eruptions, whereby the emerging flux is initially constrained by a system of overlying magnetic field lines, and is able to erupt only after an opening develops in the overlying fields as a consequence of magnetic reconnection at a magnetic null point. In our case, the EIT crinkles would be a signature of this pre-impulsive-phase magnetic reconnection, and brightening of the core only occurs after the core fields begin to escape through the newly-created opening in the overlying fields. Morphology in soft X-ray images and properties in hard X-rays differ between the two events, with complexities that preclude a simple determination of the dynamics in the core at the times of eruption. From the BCS spectra, we find that the core region expends energy at a rate of approx. 10(exp 26) erg per second during the time of the growth of the EIT crinkles; this rate is an upper limit to energy expended in the reconnections opening the overlying fields. Energy losses occur at an order-of-magnitude higher rate near the time of the peak of the events. There is little evidence of asymmetry in the spectra, consistent with the majority of the mass flows occurring normal to the line-of-sight. Both events have similar electron temperature dependencies on time.

Sterling, Alphonse C.↗
DOE OSTI2024

Structural characterization of ligand binding and pH-specific enzymatic activity of mouse Acidic Mammalian Chitinase

Chitin is an abundant biopolymer and pathogen-associated molecular pattern that stimulates a host innate immune response. Mammals express chitin-binding and chitin-degrading proteins to remove chitin from the body. One of these proteins, Acidic Mammalian Chitinase (AMCase), is an enzyme known for its ability to function under acidic conditions in the stomach but is also active in tissues with more neutral pHs, such as the lung. Here, we used a combination of biochemical, structural, and computational modeling approaches to examine how the mouse homolog (mAMCase) can act in both acidic and neutral environments. We measured kinetic properties of mAMCase activity across a broad pH range, quantifying its unusual dual activity optima at pH 2 and 7. We also solved high-resolution crystal structures of mAMCase in complex with oligomeric GlcNAcn, the building block of chitin, where we identified extensive conformational ligand heterogeneity. Leveraging these data, we conducted molecular dynamics simulations that suggest how a key catalytic residue could be protonated via distinct mechanisms in each of the two environmental pH ranges. These results integrate structural, biochemical, and computational approaches to deliver a more complete understanding of the catalytic mechanism governing mAMCase activity at different pH. Engineering proteins with tunable pH optima may provide new opportunities to develop improved enzyme variants, including AMCase, for therapeutic purposes in chitin degradation.

59 BASIC BIOLOGICAL SCIENCES↗