Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “protein function prediction”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

The E3 ligase UBR5 regulates gastric cancer cell growth by destabilizing the tumor suppressor GKN1

Gastric cancer is the most common digestive malignant tumor worldwide and the underlying mechanisms are not fully understood. The E3 ligase UBR5 (also known as EDD1) is essentially involved in diverse types of cancer. Here we aimed to study the functions of UBR5 in human gastric cancer. We first analyzed the mRNA and protein levels of UBR5 in human gastric cancer tissues and the results showed that UBR5 was markedly increased in gastric cancer tissues compared with normal gastric mucosa or matched non-cancer gastric tissues. The relationship between UBR5 and survival of gastric cancer patients was analyzed and we found that high UBR5 expression was associated with poor overall and disease-free survival. We further tried to investigate the effects of UBR5 on gastric cancer cell growth in vitro and in vivo. Therefore, we knocked down UBR5 with lentivirus-mediated shRNA and found that UBR5 knockdown repressed in vitro proliferation and colony formation of gastric cancer cells AGS, MG803 and MNK1. In vivo xenograft experiment also demonstrated that UBR5 knockdown inhibited AGS growth. Finally, we explored the mechanism by which UBR5 contributed to the growth of gastric cancer cells. We found that UBR5 bound the tumor suppressor gastrokine 1 (GKN1) and increased its ubiquitination to reduce the protein stability of GKN1. GKN1 knockdown with lentivirus-mediated shRNA increased the in vitro colony formation and in vivo growth of AGS cells, and UBR5 knockdown was unable to affect the colony formation and in vivo growth of AGS cells when GKN1 was knocked down, indicating that GKN1 contributed to the effects of UBR5 in human gastric cancer cells. Taken together, UBR5 plays an essential role in gastric cancer and may be a potential diagnosis and treatment target for gastric cancer. - Highlights: • UBR5 expression is up-regulated in human gastric cancer. • UBR5 overexpression predicts poor survival. • UBR5 regulates gastric cancer growth in vitro and in vivo. • UBR5 reduces the protein stability of GKN1.

60 APPLIED LIFE SCIENCES↗
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 OSTI2015

Nuclear localization of the DNA repair scaffold XRCC1: Uncovering the functional role of a bipartite NLS

We have characterized the nuclear localization signal (NLS) of XRCC1 structurally using X-ray crystallography and functionally using fluorescence imaging. Crystallography and binding studies confirm the bipartite nature of the XRCC1 NLS interaction with Importin α (Impα) in which the major and minor binding motifs are separated by >20 residues, and resolve previous inconsistent determinations. Binding studies of peptides corresponding to the bipartite NLS, as well as its major and minor binding motifs, to both wild-type and mutated forms of Impα reveal pronounced cooperative binding behavior that is generated by the proximity effect of the tethered major and minor motifs of the NLS. The cooperativity stems from the increased local concentration of the second motif near its cognate binding site that is a consequence of the stepwise binding behavior of the bipartite NLS. We predict that the stepwise dissociation of the NLS from Impα facilitates unloading by providing a partially complexed intermediate that is available for competitive binding by Nup50 or the Importin β binding domain. This behavior gives a basis for meeting the intrinsically conflicting high affinity and high flux requirements of an efficient nuclear transport system.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Streptomyces wadayamensis MppP Is a Pyridoxal 5'-Phosphate-Dependent L -Arginine α-Deaminase, γ-Hydroxylase in the Enduracididine Biosynthetic Pathway

l-Enduracididine (l-End) is a nonproteinogenic amino acid found in a number of bioactive peptides, including the antibiotics teixobactin, enduracidin, and mannopeptimycin. The potent activity of these compounds against antibiotic-resistant pathogens like MRSA and their novel mode of action have garnered considerable interest for the development of these peptides into clinically relevant antibiotics. Here, this goal has been hampered, at least in part, by the fact that l-End is difficult to synthesize and not currently commercially available. We have begun to elucidate the biosynthetic pathway of this unusual building block. In mannopeptimycin-producing strains, like Streptomyces wadayamensis, l-End is produced from l-Arg by the action of three enzymes: MppP, MppQ, and MppR. Herein, we report the structural and functional characterization of MppP. This pyridoxal 5'-phosphate (PLP)-dependent enzyme was predicted to be a fold type I aminotransferase on the basis of sequence analysis. We show that MppP is actually the first example of a PLP-dependent hydroxylase that catalyzes a reaction of l-Arg with dioxygen to yield a mixture of 2-oxo-4-hydroxy-5-guanidinovaleric acid and 2-oxo-5-guanidinovaleric acid in a 1.7:1 ratio. The structure of MppP with PLP bound to the catalytic lysine residue (Lys221) shows that, while the tertiary structure is very similar to those of the well-studied aminotransferases, there are differences in the arrangement of active site residues around the cofactor that likely account for the unusual activity of this enzyme. The structure of MppP with the substrate analogue d-Arg bound shows how the enzyme binds its substrate and indicates why d-Arg is not a substrate. On the basis of this work and previous work with MppR, we propose a plausible biosynthetic scheme for l-End.

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

Secondary structure of rat and human amylin across force fields

The aggregation of human amylin has been strongly implicated in the progression of Type II diabetes. This 37-residue peptide forms a variety of secondary structures, including random coils, α-helices, and β-hairpins. The balance between these structures depends on the chemical environment, making amylin an ideal candidate to examine inherent biases in force fields. Rat amylin differs from human amylin by only 6 residues; however, it does not form fibrils. Therefore it provides a useful complement to human amylin in studies of the key events along the aggregation pathway. In this work, the free energy of rat and human amylin was determined as a function of α-helix and β-hairpin content for the Gromos96 53a6, OPLS-AA/L, CHARMM22/CMAP, CHARMM22*, Amberff99sb*-ILDN, and Amberff03w force fields using advanced sampling techniques, specifically bias exchange metadynamics. This work represents a first systematic attempt to evaluate the conformations and the corresponding free energy of a large, clinically relevant disordered peptide in solution across force fields. The NMR chemical shifts of rIAPP were calculated for each of the force fields using their respective free energy maps, allowing us to quantitatively assess their predictions. We show that the predicted distribution of secondary structures is sensitive to the choice of force-field: Gromos53a6 is biased towards β-hairpins, while CHARMM22/CMAP predicts structures that are overly α-helical. OPLS-AA/L favors disordered structures. Amberff99sb*-ILDN, AmberFF03w and CHARMM22* provide the balance between secondary structures that is most consistent with available experimental data. In contrast to previous reports, our findings suggest that the equilibrium conformations of human and rat amylin are remarkably similar, but that subtle differences arise in transient alpha-helical and beta-strand containing structures that the human peptide can more readily adopt. We hypothesize that these transient states enable dynamic pathways that facilitate the formation of aggregates and, eventually, amyloid fibrils.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

The Role of Light–Dark Regulation of the Chloroplast ATP Synthase

The chloroplast ATP synthase catalyzes the light-driven synthesis of ATP and is activated in the light and inactivated in the dark by redox-modulation through the thioredoxin system. It has been proposed that this down-regulation is important for preventing wasteful hydrolysis of ATP in the dark. To test this proposal, we compared the effects of extended dark exposure in Arabidopsis lines expressing the wild-type and mutant forms of ATP synthase that are redox regulated or constitutively active. In contrast to the predictions of the model, we observed that plants with wild-type redox regulation lost photosynthetic capacity rapidly in darkness, whereas those expressing redox-insensitive form were far more stable. To explain these results, we propose that in wild-type plants, down-regulation of ATP synthase inhibits ATP hydrolysis, leading to dissipation of thylakoid proton motive force (pmf) and subsequent inhibition of protein transport across the thylakoid through the twin arginine transporter (Tat)-dependent and Secdependent import pathways, resulting in the selective loss of specific protein complexes. By contrast, in mutants with a redox-insensitive ATP synthase, pmf is maintained by ATP hydrolysis, thus allowing protein transport to maintain photosynthetic activities for extended periods in the dark. Hence, a basal level of Tat-dependent, as well as, Sec-dependent import activity, in the dark helps replenishes certain components of the photosynthetic complexes and thereby aids in maintaining overall complex activity. But, the influence of a dark pmf on thylakoid protein import, by itself, could not explain all the effects we observed in this study. For example, we also observed in wild type plants a large transient buildup of thylakoid pmf and nonphotochemical exciton quenching upon sudden illumination of dark adapted plants. Thus, we conclude that down-regulation of the ATP synthase is probably not related to preventing loss of ATP per se. Instead, ATP synthase redox regulation may be impacting a number of cellular processes such as (1) the accumulation of chloroplast proteins and/or ions or (2) the responses of photosynthesis to rapid changes in light intensity. A model highlighting the complex interplay between ATP synthase regulation and pmf in maintaining various chloroplast functions in the dark is presented.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

RNA Sequencing Identifies New RNase III Cleavage Sites in Escherichia coli and Reveals Increased Regulation of mRNA

ABSTRACT Ribonucleases facilitate rapid turnover of RNA, providing cells with another mechanism to adjust transcript and protein levels in response to environmental conditions. While many examples have been documented, a comprehensive list of RNase targets is not available. To address this knowledge gap, we compared levels of RNA sequencing coverage of Escherichia coli and a corresponding RNase III mutant to expand the list of known RNase III targets. RNase III is a widespread endoribonuclease that binds and cleaves double-stranded RNA in many critical transcripts. RNase III cleavage at novel sites found in aceEF , proP , tnaC , dctA , pheM , sdhC , yhhQ , glpT , aceK , and gluQ accelerated RNA decay, consistent with previously described targets wherein RNase III cleavage initiates rapid degradation of secondary messages by other RNases. In contrast, cleavage at three novel sites in the ahpF , pflB , and yajQ transcripts led to stabilized secondary transcripts. Two other novel sites in hisL and pheM overlapped with transcriptional attenuators that likely serve to ensure turnover of these highly structured RNAs. Many of the new RNase III target sites are located on transcripts encoding metabolic enzymes. For instance, two novel RNase III sites are located within transcripts encoding enzymes near a key metabolic node connecting glycolysis and the tricarboxylic acid (TCA) cycle. Pyruvate dehydrogenase activity was increased in an rnc deletion mutant compared to the wild-type (WT) strain in early stationary phase, confirming the novel link between RNA turnover and regulation of pathway activity. Identification of these novel sites suggests that mRNA turnover may be an underappreciated mode of regulating metabolism. IMPORTANCE The concerted action and overlapping functions of endoribonucleases, exoribonucleases, and RNA processing enzymes complicate the study of global RNA turnover and recycling of specific transcripts. More information about RNase specificity and activity is needed to make predictions of transcript half-life and to design synthetic transcripts with optimal stability. RNase III does not have a conserved target sequence but instead recognizes RNA secondary structure. Prior to this study, only a few RNase III target sites in E. coli were known, so we used RNA sequencing to provide a more comprehensive list of cleavage sites and to examine the impact of RNase III on transcript degradation. With this added information on how RNase III participates in transcript regulation and recycling, a more complete picture of RNA turnover can be developed for E. coli . Similar approaches could be used to augment our understanding of RNA turnover in other bacteria.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Supramolecular Interactions and Morphology of Self-Assembling Peptide Amphiphile Nanostructures

The morphology of supramolecular peptide nanostructures is difficult to predict given their complex energy landscapes. We investigated peptide amphiphiles containing β-sheet forming domains that form twisted nanoribbons in water. We explained the morphology based on a balance between the energetically favorable packing of molecules in the center of the nanostructures, the unfavorable packing at the edges, and the deformations due to packing of twisted β-sheets. We find that morphological polydispersity of PA nanostructures is determined by peptide sequences, and the twisting of their internal β-sheets. Here, we also observed a change in the supramolecular chirality of the nanostructures as the peptide sequence was modified, although only amino acids with l-configuration were used. Upon increasing charge repulsion between molecules, we observed a change in morphology to long cylinders and then rodlike fragments and spherical micelles. Understanding the self-assembly mechanisms of peptide amphiphiles into nanostructures should be useful to optimize their well-known functions.

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

Functional Design of Peptide Materials Based on Supramolecular Cohesion

Peptide materials offer a broad platform to design biomimetic soft matter, and filamentous networks that emulate those in extracellular matrices and the cytoskeleton are among the important targets. Given the vast sequence space, a combination of computational approaches and readily accessible experimental techniques is required to design peptide materials efficiently. Here, we report here on a strategy that utilizes this combination to predict supramolecular cohesion within filaments of peptide amphiphiles, a property recently linked to supramolecular dynamics and consequently bioactivity. Using established coarse-grained simulations on 10,000 randomly generated peptide sequences, we identified 3500 likely to self-assemble in water into nanoscale filaments. Atomistic simulations of small clusters were used to further analyze this subset of sequences and identify mathematical descriptors that are predictive of intermolecular cohesion, which was the main purpose of this work. We arbitrarily selected a small cohort of these sequences for chemical synthesis and verified their fiber morphology. With further characterization, we were able to link the latent heat associated with fiber to micelle transitions, an indicator of cohesion and potential supramolecular dynamicity within the filaments, to calculated hydrogen bond densities in the simulation clusters. Based on validation from in situ synchrotron X-ray scattering and differential scanning calorimetry, we conclude that the phase transitions can be easily observed by very simple polarized light microscopy experiments. We are encouraged by the methodology explored here as a relatively low-cost and fast way to design potential functions of peptide materials.

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

MD/DPD Multiscale Framework for Predicting Morphology and Stresses of Red Blood Cells in Health and Disease

Healthy red blood cells (RBCs) have remarkable deformability, squeezing through narrow capillaries as small as 3 microns in diameter without any damage. However, in many hematological disorders the spectrin network and lipid bilayer of diseased RBCs may be significantly altered, leading to impaired functionality including loss of deformability. We employ a two-component whole-cell multiscale model to quantify the biomechanical characteristics of the healthy and diseased RBCs, including Plasmodium falciparum-infected RBCs (Pf-RBCs) and defective RBCs in hereditary disorders, such as spherocytosis and elliptocytosis. In particular, we develop a two-step multiscale framework based on coarse-grained molecular dynamics (CGMD) and dissipative particle dynamics (DPD) to predict the static and dynamic responses of RBCs subject to tensile forcing, using experimental information only on the structural defects in the lipid bilayer, cytoskeleton, and their interaction. We first employ CGMD on a small RBC patch to compute the shear modulus, bending stiffness, and network parameters, which are subsequently used as input to a whole-cell DPD model to predict the RBC shape and corresponding stress field. For Pf-RBCs at trophozoite and schizont stages, the presence of cytoadherent knobs elevates the shear response in the lipid bilayer and stiffens the RBC membrane. For RBCs in spherocytosis and elliptocytosis, the bilayer-cytoskeleton interaction is weakened, resulting in substantial increase of the tensile stress in the lipid bilayer. Furthermore, we investigate the transient behavior of stretching deformation and shape relaxation of the normal and defective RBCs. Different from the normal RBCs possessing high elasticity, our simulations reveal that the defective RBCs respond irreversibly, i.e., they lose their ability to recover the normal biconcave shape in successive loading cycles of stretching and relaxation. Our findings provide fundamental insights into the microstructure and biomechanics of RBCs, and demonstrate that the two-step multiscale framework presented here can be used effectively for in silico studies of hematological disorders based on first principles and patient-specific experimental input at the protein level.

Biochemistry & Molecular Biology↗
NASA NTRS2005

Origin of the Eumetazoa: testing ecological predictions of molecular clocks against the Proterozoic fossil record

Molecular clocks have the potential to shed light on the timing of early metazoan divergences, but differing algorithms and calibration points yield conspicuously discordant results. We argue here that competing molecular clock hypotheses should be testable in the fossil record, on the principle that fundamentally new grades of animal organization will have ecosystem-wide impacts. Using a set of seven nuclear-encoded protein sequences, we demonstrate the paraphyly of Porifera and calculate sponge/eumetazoan and cnidarian/bilaterian divergence times by using both distance [minimum evolution (ME)] and maximum likelihood (ML) molecular clocks; ME brackets the appearance of Eumetazoa between 634 and 604 Ma, whereas ML suggests it was between 867 and 748 Ma. Significantly, the ME, but not the ML, estimate is coincident with a major regime change in the Proterozoic acritarch record, including: (i) disappearance of low-diversity, evolutionarily static, pre-Ediacaran acanthomorphs; (ii) radiation of the high-diversity, short-lived Doushantuo-Pertatataka microbiota; and (iii) an order-of-magnitude increase in evolutionary turnover rate. We interpret this turnover as a consequence of the novel ecological challenges accompanying the evolution of the eumetazoan nervous system and gut. Thus, the more readily preserved microfossil record provides positive evidence for the absence of pre-Ediacaran eumetazoans and strongly supports the veracity, and therefore more general application, of the ME molecular clock.

NASA Discipline Evolutionary Biology↗
DOE OSTI2018

MiR-429 suppresses neurotrophin-3 to alleviate perineural invasion of pancreatic cancer

Highlights: • MiR-429 is downregulated in pancreatic cancer tissues and cells. • MiR-429 functions to inhibit cancer cell proliferation and invasion. • MiR-429 directly targets to NT-3 mRNA. • MiR-429 suppress NT-3 in cancer cells to resist the attraction TrkC in PC12 cells. • NT-3-TrkC signaling promotes PC12 neurite outgrowth and neuropathic pain. Perineural invasion (PNI) potentially increases the risk of relapse and abdominal pain in patients with pancreatic ductal adenocarcinoma (PDAC). However, the underlying mechanisms of PNI of PDAC is incompletely revealed. Our study aimed to investigate roles of miR-429 in modulating PNI in PDAC. We found that miR-429 was downregulated in PDAC cancer tissues and was profoundly decreased in tissues with PNI. It was reduced in nine of the ten examined pancreatic cancer cell lines. MiR-429 mimics restored its cellular expressions in MIA PaCa-2 and BxCP3 cells and significantly suppressed cell viability and invasion of the cancer cells. The online bioinformatic software predicted that neurotrophin-3 (NT-3) was a potential target gene of miR-429. It was showed that NT-3 mRNA elevated in PC cancer tissues, especially in patients presenting PNI. MiR-429 upregulation substantially suppressed the NT-3 mRNA and secretion in cancer cells. Also, the dual luciferase reporter assays confirmed the interaction between miR-429 and NT-3. When co-culturing the two PDAC cells with PC-12 cells, the invaded cell counts significantly increased comparing with the sole culture of cancer cells. However, miR-429 mimic transfection or NT-3 blocking retarded the cancer invasion in the co-culture system. Besides, we found that cancer cells conditioned medium (CM) treatment significantly increased the neurite outgrowth percentage in PC-12 cells, which was suppressed by culturing with CM from miR-429 mimics-transfected cells. In the CM cultured PC-12 cells, NT-3 receptor TrkC as well as pain-related proteins TRPV1 and TRPV2 significantly elevated. Collectively, miR-429 potentially suppressed neurotrophin-3 to alleviate PNI of PDAC.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Molecular basis of P[II] major human rotavirus VP8* domain recognition of histo-blood group antigens

Initial cell attachment of rotavirus (RV) to specific cell surface glycan receptors, which is the essential first step in RV infection, is mediated by the VP8* domain of the spike protein VP4. Recently, human histo-blood group antigens (HBGAs) have been identified as receptors or attachment factors for human RV strains. RV strains in the P[4] and P[8] genotypes of the P[II] genogroup share common recognition of the Lewis b (Le b ) and H type 1 antigens, however, the molecular basis of receptor recognition by the major human P[8] RVs remains unknown due to lack of experimental structural information. Here, we used nuclear magnetic resonance (NMR) spectroscopy-based titration experiments and NMR-derived high ambiguity driven docking (HADDOCK) methods to elucidate the molecular basis for P[8] VP8* recognition of the Leb (LNDFH I) and type 1 HBGAs. We also used X-ray crystallography to determine the molecular details underlying P[6] recognition of H type 1 HBGAs. Unlike P[6]/P[19] VP8*s that recognize H type 1 HBGAs in a binding surface composed of an α-helix and a β-sheet, referred as the “βα binding site”, the P[8] and P[4] VP8*s bind Le b HBGAs in a previously undescribed pocket formed by the edges of two β-sheets, referred to as the “ββ binding site”. Importantly, the P[8] and P[4] VP8*s retain binding capability to non-Le b type 1 HBGAs using the βα binding site. The presence of two distinct binding sites for Le b and non-Le b HBGA glycans in the P[8] and P[4] VP8* domains suggests host-pathogen co-evolution under structural and functional adaptation of RV pathogens to host glycan polymorphisms. Assessment and understanding of the precise impact of this co-evolutionary process in determining RV host ranges and cross-species RV transmission should facilitate improved RV vaccine development and prediction of future RV strain emergence and epidemics.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

A Surface Exposed, Two-Domain Lipoprotein Cargo of a Type XI Secretion System Promotes Colonization of Host Intestinal Epithelia Expressing Glycans

The only known required component of the newly described Type XI secretion system (TXISS) is an outer membrane protein (OMP) of the DUF560 family. TXISS OMPs are broadly distributed across proteobacteria, but properties of the cargo proteins they secrete are largely unexplored. We report biophysical, histochemical, and phenotypic evidence that Xenorhabdus nematophila NilC is surface exposed. Biophysical data and structure predictions indicate that NilC is a two-domain protein with a C-terminal, 8-stranded β-barrel. This structure has been noted as a common feature of TXISS effectors and may be important for interactions with the TXISS OMP . The NilC N-terminal domain is more enigmatic, but our results indicate it is ordered and forms a β-sheet structure, and bioinformatics suggest structural similarities to carbohydrate-binding proteins. X. nematophila NilC and its presumptive TXISS OMP partner NilB are required for colonizing the anterior intestine of Steinernema carpocapsae nematodes: the receptacle of free-living, infective juveniles and the anterior intestinal cecum (AIC) in juveniles and adults. We show that, in adult nematodes, the AIC expresses a Wheat Germ Agglutinin (WGA)-reactive material, indicating the presence of N-acetylglucosamine or N-acetylneuraminic acid sugars on the AIC surface. A role for this material in colonization is supported by the fact that exogenous addition of WGA can inhibit AIC colonization by X. nematophila. Conversely, the addition of exogenous purified NilC increases the frequency with which X. nematophila is observed at the AIC, demonstrating that abundant extracellular NilC can enhance colonization. NilC may facilitate X. nematophila adherence to the nematode intestinal surface by binding to host glycans, it might support X. nematophila nutrition by cleaving sugars from the host surface, or it might help protect X. nematophila from nematode host immunity. Proteomic and metabolomic analyses of wild type X. nematophila compared to those lacking nilB and nilC revealed differences in cell wall and secreted polysaccharide metabolic pathways. Additionally, purified NilC is capable of binding peptidoglycan, suggesting that periplasmic NilC may interact with the bacterial cell wall. Overall, these findings support a model that NilB-regulated surface exposure of NilC mediates interactions between X. nematophila and host surface glycans during colonization. This is a previously unknown function for a TXISS.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2012

Multiscale Modeling of Metabolism and Macromolecular Synthesis in E. coli and Its Application to the Evolution of Codon Usage

Biological systems are inherently hierarchal and multiscale in time and space. A major challenge of systems biology is to describe biological systems as a computational model, which can be used to derive novel hypothesis and drive experiments leading to new knowledge. The constraint-based reconstruction and analysis approach has been successfully applied to metabolism and to the macromolecular synthesis machinery assembly. Here, we present the first integrated stoichiometric multiscale model of metabolism and macromolecular synthesis for Escherichia coli K12 MG1655, which describes the sequence-specific synthesis and function of almost 2000 gene products at molecular detail. We added linear constraints, which couple enzyme synthesis and catalysis reactions. Comparison with experimental data showed improvement of growth phenotype prediction with the multiscale model over E. coli’s metabolic model alone. Many of the genes covered by this integrated model are well conserved across enterobacters and other, less related bacteria. We addressed the question of whether the bias in synonymous codon usage could affect the growth phenotype and environmental niches that an organism can occupy. We created two classes of in silico strains, one with more biased codon usage and one with more equilibrated codon usage than the wildtype. The reduced growth phenotype in biased strains was caused by tRNA supply shortage, indicating that expansion of tRNA gene content or tRNA codon recognition allow E. coli to respond to changes in codon usage bias. Our analysis suggests that in order to maximize growth and to adapt to new environmental niches, codon usage and tRNA content must co-evolve. These results provide further evidence for the mutation-selection-drift balance theory of codon usage bias. This integrated multiscale reconstruction successfully demonstrates that the constraintbased modeling approach is well suited to whole-cell modeling endeavors.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

Electronic Energies Are Not Enough: An Ion Mobility-Aided, Quantum Chemical Benchmark Analysis of H + GPGG Conformers

Peptides and proteins exist not as single, lowest energy structures but as ensembles of states separated by small barriers. In order to study these species we must be able to correctly identify their gas phase conformational distributions, and ion mobility spectrometry (IMS) has arisen as an experimental method for assessing the gas phase energetics of flexible peptides. Here, we present a thorough exploration and benchmarking of the low energy conformers of the small, hairpin peptide H+ GPGG with the aid of ion mobility spectrometry against a wide swath of density functionals (35 dispersion-corrected and uncorrected functionals represented by 21 unique exchange–correlation functionals) and wave function theory methods (15 total levels of theory). The three experimentally resolved IMS peaks were found to correspond to three distinct pairs of conformers, each pair composed of species differing only by the chair or boat configuration of the proline. Two of the H+ GPGG conformer pairs possess a cis configuration about the Pro–Gly1 peptide bond while the other adopts a trans configuration. While the experimental spectrum reports a higher intensity for the cis-1 conformer than the trans conformer, 13 WFT levels of theory, including a complete basis set CCSD(T) extrapolation, obtain trans to be favored in terms of the electronic energy. This same effect is seen in 14 of the 18 dispersion-corrected density functionals studied, whereas the remaining 17 functionals show more variety. Only when Gibbs free energies are considered do the WFT methods and dispersion-corrected functionals reflect the experimental distribution. CAM-B3LYP-D3BJ emerges as the best-performing density functional, matching the experimental distribution and the CCSD(T)/CBS relative energies within 6%. Further analysis reveals the trans conformer to be favored electronically, but entropically disfavored, leading to the experimental preference of the cis-1 conformer. Furthermore, these results highlight the danger in considering only the electronic energies, which is common practice in electronic structure theory predictions of conformational energy distributions. Additionally, the effects of temperature and scaling of the frequencies used in obtaining the Gibbs free energies are explored.

Energy↗
DOE OSTI2014

Conserved active site cysteine residue of archaeal THI4 homolog is essential for thiamine biosynthesis in Haloferax volcanii

Thiamine (vitamin B1) is synthesized de novo by certain yeast, fungi, plants, protozoans, bacteria and archaea. The pathway of thiamine biosynthesis by archaea is poorly understood, particularly the route of sulfur relay to form the thiazole ring. Archaea harbor structural homologs of both the bacterial (ThiS-ThiF) and eukaryotic (THI4) proteins that mobilize sulfur to thiazole ring precursors by distinct mechanisms. Results: Based on comparative genome analysis, halophilic archaea are predicted to synthesize the pyrimidine moiety of thiamine by the bacterial pathway, initially suggesting that also a bacterial ThiS-ThiF type mechanism for synthesis of the thiazole ring is used in which the sulfur carrier ThiS is first activated by ThiF-catalyzed adenylation. The only ThiF homolog of Haloferax volcanii (UbaA) was deleted but this had no effect on growth in the absence of thiamine. Usage of the eukaryotic THI4-type sulfur relay was initially considered less likely for thiamine biosynthesis in archaea, since the active-site cysteine residue of yeast THI4p that donates the sulfur to the thiazole ring by a suicide mechanism is replaced by a histidine residue in many archaeal THI4 homologs and these are described as D-ribose-1,5-bisphosphate isomerases. The THI4 homolog of the halophilic archaea, including Hfx. volcanii (HVO_0665, HvThi4) was found to differ from that of methanogens and thermococci by having a cysteine residue (Cys165) corresponding to the conserved active site cysteine of yeast THI4p (Cys205). Deletion of HVO_0665 generated a thiamine auxotroph that was trans-complemented by a wild-type copy of HVO_0665, but not the modified gene encoding an HvThi4 C165A variant. Based on our results, we conclude that the archaeon Hfx. volcanii uses a yeast THI4-type mechanism for sulfur relay to form the thiazole ring of thiamine. We extend this finding to a relatively large group of archaea, including haloarchaea, ammonium oxidizing archaea, and some methanogen and Pyrococcus species, by observing that these organisms code for THI4 homologs that have a conserved active site cysteine residue which is likely used in thiamine biosynthesis. Thus, archaeal members of IPR002922 THI4 family that have a conserved cysteine active site should be reexamined for a function in thiamine biosynthesis.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Neutrons in Structural Biology: Challenges and Opportunities (Workshop Report)

Gaining a thorough understanding of biological systems requires building our knowledge about biological processes from the level of atoms and electrons, and up to whole organisms. Such comprehensive knowledge will allow for a predictive understanding of complex biological systems behavior. It will guide us in the design and development of novel therapeutics and vaccines to tackle existing health threats and to prepare for future pandemics, and it will provide information necessary to create new biomaterials and bio-inspired technologies through manipulation of biological macromolecules, their assemblies, single cells and even microorganisms. Reaching these goals will require a synergistic combination of multiple experimental techniques with molecular calculations and predictive simulations, and the design and development of new techniques and capabilities that bridge current knowledge and technology gaps. Neutron scattering provides unique information about the biomacromolecular structure and function and can play a major role in achieving these goals. A workshop was held to engage the scientific community in identifying pressing challenges in biochemistry, structural biology, enzymology and structure-guided drug design not solved with the current neutron scattering technologies or utilizing other structural biology techniques such as X-ray crystallography, NMR, and cryo-EM. The workshop brought together structural biology, biochemistry and computational experts, as well as early career researchers and students, creating a forum for discussing scientific advancement and collaboration. The workshop included a one-day satellite training workshop where graduate students and postdoctoral researchers were educated in the application of neutron crystallography and small-angle scattering in structural biology. Furthermore, the Instrument Scientific Advisory Board (ISAB) for the development of a macromolecular neutron diffractometer at ORNL’s Second Target Station was introduced at the workshop. The major outcome was that neutrons can provide atomic-level understanding of biomacromolecular structure, function and dynamics which is of paramount importance for addressing the identified challenges. Neutron crystallography, in particular, can resolve long-standing biochemical issues regarding enzyme function by delineating the underlying chemistry and can have a major impact on the design of small-molecule therapeutics, especially in combination with molecular computation (quantum chemistry and molecular dynamics simulations) and the emerging artificial intelligence (AI)-assisted drug design technologies. The unique properties of neutrons, including their high sensitivity to hydrogen and their non-destructive nature, make them ideal probes of biological matter. There is a palpable need in the scientific community to expand and enhance the impact of neutron sciences on biology. Neutron crystallography is the only structural biology method capable of determining positions of all hydrogen atoms in proteins, nucleic acids and their complexes at near-physiological temperatures and of unstable species at cryogenic temperatures. Moreover, neutron analysis is non-ionizing, non-destructive and does not perturb the structure or redox chemistry of active site metal centers and clusters in proteins, which can be invaluable for studying radiation-sensitive metalloprotein complexes. Further, neutron energies used in scattering applications are similar to atomic motions, permitting neutron spectroscopies to characterize the dynamics of biomacromolecules on the picosecond to microsecond timescales. The different sensitivities of neutrons to protium (H) and deuterium (D) isotopes of hydrogen allow enhanced visibility of specific parts of biological complexes through isotopic labeling. The impact of neutrons will be most powerful when neutron scattering is combined with complementary experimental techniques that use photons and electrons, and with high-performance computing. The interconnection and mutuality of the experimental and theoretical capabilities will drive discoveries in biological and health sciences to generate more complete picture of complex biological systems. The major limitation in the field of biological neutron crystallography has been signal-to-noise, demanding large samples that are difficult to produce for the majority of biomacromolecules and limiting the applicability of this technique in biological sciences. A neutron crystallography instrument at the Second Target Station will revolutionize biological science with neutrons by engaging a large scientific community of structural biologists, enabling successful neutron diffraction experiments from radically smaller biomacromolecular crystals, resolving unanswered biochemical questions, and meaningfully contributing to rational drug design. The meeting highlighted 10 grand challenges that will be addressed with this advanced capability over the next decade and beyond, and the recommendations required to help address them are given below.

59 BASIC BIOLOGICAL SCIENCES↗