Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational protein design”

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 199 records · Page 11

Probing the SARS-CoV-2 main protease active site as a target to inhibit viral replication

SARS-CoV-2 replication involves the synthesis of two large proteins, which are inactive and harmless until the viral main protease enzyme (3CL M pro ) uses them as substrates, cleaving them into many smaller, functional products. Dr Andrey Kovalevsky and his team from Oak Ridge National Laboratory propose a design of novel inhibitors and the repurposing of clinical drugs developed to treat other diseases for the treatment of COVID-19. Finally, the team uses room temperature crystallography with X-rays and neutrons to guide structure-based and computer-assisted drug design and to assess the ability of drugs to inhibit the SARS-CoV-2 main protease, causing virus replication to stop.

Kneller, Daniel↗

Proteins with Novel Structure, Function and Dynamics

Recently, a small enzyme that ligates two RNA fragments with the rate of 10(exp 6) above background was evolved in vitro (Seelig and Szostak, Nature 448:828‐831, 2007). This enzyme does not resemble any contemporary protein (Chao et al., Nature Chem. Biol. 9:81‐83, 2013). It consists of a dynamic, catalytic loop, a small, rigid core containing two zinc ions coordinated by neighboring amino acids, and two highly flexible tails that might be unimportant for protein function. In contrast to other proteins, this enzyme does not contain ordered secondary structure elements, such as alpha‐helix or beta‐sheet. The loop is kept together by just two interactions of a charged residue and a histidine with a zinc ion, which they coordinate on the opposite side of the loop. Such structure appears to be very fragile. Surprisingly, computer simulations indicate otherwise. As the coordinating, charged residue is mutated to alanine, another, nearby charged residue takes its place, thus keeping the structure nearly intact. If this residue is also substituted by alanine a salt bridge involving two other, charged residues on the opposite sides of the loop keeps the loop in place. These adjustments are facilitated by high flexibility of the protein. Computational predictions have been confirmed experimentally, as both mutants retain full activity and overall structure. These results challenge our notions about what is required for protein activity and about the relationship between protein dynamics, stability and robustness. We hypothesize that small, highly dynamic proteins could be both active and fault tolerant in ways that many other proteins are not, i.e. they can adjust to retain their structure and activity even if subjected to mutations in structurally critical regions. This opens the doors for designing proteins with novel functions, structures and dynamics that have not been yet considered.

Proteins↗

Cholesterol modulates membrane elasticity via unified biophysical laws

Cholesterol and lipid unsaturation underlie a balance of opposing forces that features prominently in adaptive cell responses to diet and environmental cues. These competing factors have resulted in contradictory observations of membrane elasticity across different measurement scales, requiring chemical specificity to explain incompatible structural and elastic effects. Here, we demonstrate that – unlike macroscopic observations – lipid membranes exhibit a unified elastic behavior in the mesoscopic regime between molecular and macroscopic dimensions. Using nuclear spin techniques and computational analysis, we find that mesoscopic bending moduli follow a universal dependence on the lipid packing density regardless of cholesterol content, lipid unsaturation, or temperature. Our observations reveal that compositional complexity can be explained by simple biophysical laws that directly map membrane elasticity to molecular packing associated with biological function, curvature transformations, and protein interactions. The obtained scaling laws closely align with theoretical predictions based on conformational chain entropy and elastic stress fields. These findings provide unique insights into the membrane design rules optimized by nature and unlock predictive capabilities for guiding the functional performance of lipid-based materials in synthetic biology and real-world applications.

Kumarage, Teshani [Virginia Polytechnic Inst. and ↗

Enhanced polymorph metastability drives glycine nucleation in aqueous salt solutions

Crystal nucleation from aqueous solutions influences countless geological, biochemical, astrophysical, environmental, and materials science–related phenomena, including ice formation, the manufacturing of active pharmaceutical ingredients, development of diseases such as Alzheimer’s and the origin of life itself. Understanding and controlling nucleation is essential for designing materials with specific properties, developing strategies to inhibit or promote crystallization in various contexts and preventing pathological aggregation in neurodegenerative diseases. Similar to the protein structure prediction problem—where a single amino acid sequence can in theory adopt one most stable conformation but in practice may sample multiple competing conformations—crystal nucleation faces a parallel challenge: the same chemical species can form diverse polymorphs under different environmental conditions (e.g., temperature, pressure, solvent). Each polymorph presents its own set of physical and chemical properties, highlighting the importance of understanding and controlling polymorph selection in fields ranging from pharmaceuticals to materials design. Despite advances in experimental and computational methods for studying phase transitions and polymorph stability, nucleation remains challenging due to its nanoscale nature. Furthermore, in practical settings, salts and impurities can further influence crystal nucleation in diverse contexts, from scaling in pipelines and desalination plants to the durability of concrete and the efficiency of battery materials. This can lead to the formation of polymorphs that may differ from the most stable phase in pure solutions. Or, even though the final structure might appear same irrespective of whether the environment contained impurities or not, the mechanism through which it was formed might be completely different and not intuitive.

Wang, Ruiyu [University of Maryland, College Park,↗

“Off‐Label Use” of the Siderophore Enterobactin Enables Targeted Imaging of Cancer with Radioactive Ti (IV)

The development of inert, biocompatible chelation methods is required to harness the emerging positron emitting radionuclide 45 Ti for radiopharmaceutical applications. Herein, we evaluate the Ti (IV) -coordination chemistry of four catechol-based, hexacoordinate chelators using synthetic, structural, computational, and radiochemical approaches. The siderophore enterobactin (Ent) and its synthetic mimic TREN-CAM readily form mononuclear Ti (IV) species in aqueous solution at neutral pH. Radiolabeling studies reveal that Ent and TREN-CAM form mononuclear complexes with the short-lived, positron-emitting radionuclide 45 Ti (IV) , and do not transchelate to plasma proteins in vitro and exhibit rapid renal clearance in naïve mice. These features guide efforts to target the 45 Ti isotope to prostate cancer tissue through the design, synthesis, and evaluation of Ent-DUPA, a small molecule conjugate composed of a prostate specific membrane antigen (PSMA) targeting peptide and a monofunctionalized Ent scaffold. Additionally, the [ 45 Ti][Ti(Ent-DUPA)] 2− complex forms readily at room temperature. In a tumor xenograft model in mice, selective tumor tissue accumulation (8±5 %, n=5), and low off-target uptake in other organs is observed. Overall, this work demonstrates targeted imaging with 45 Ti (IV) , provides a foundation for advancing the application of 45 Ti in nuclear medicine, and reveals that Ent can be repurposed as a 45 Ti-complexing cargo for targeted nuclear imaging applications.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

2 H -Thiopyran-2-thione sulfine, a compound for converting H 2 S to HSOH/H 2 S 2 and increasing intracellular sulfane sulfur levels

Reactive sulfane sulfur species such as persulfides (RSSH) and H 2 S 2 are important redox regulators and closely linked to H 2 S signaling. However, the study of these species is still challenging due to their instability, high reactivity, and the lack of suitable donors to produce them. Herein we report a unique compound, 2H-thiopyran-2-thione sulfine (TTS), which can specifically convert H 2 S to HSOH, and then to H 2 S 2 in the presence of excess H 2 S. Meanwhile, the reaction product 2H-thiopyran-2-thione (TT) can be oxidized to reform TTS by biological oxidants. The reaction mechanism of TTS is studied experimentally and computationally. TTS can be conjugated to proteins to achieve specific delivery, and the combination of TTS and H 2 S leads to highly efficient protein persulfidation. When TTS is applied in conjunction with established H 2 S donors, the corresponding donors of H 2 S 2 (or its equivalents) are obtained. Cell-based studies reveal that TTS can effectively increase intracellular sulfane sulfur levels and compensate for certain aspects of sulfide:quinone oxidoreductase (SQR) deficiency. These properties make TTS a conceptually new strategy for the design of donors of reactive sulfane sulfur species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Novel small molecule FGF 23 inhibitors increase serum phosphate and improve skeletal abnormalities in Hyp mice

We report that excess fibroblast growth factor 23 (FGF23) causes hereditary hypophosphatemic rickets, such as X-linked hypophosphatemia (XLH) and tumor induced osteomalacia (TIO). A small molecule that specifically binds to FGF23 to prevent activation of the FGFR/a-Klotho complex has potential advantages over the currently approved systemically administered FGF23 blocking antibody. Using structure-based drug design we previously identified ZINC13407541 (N-[[2-(2-phenylethenyl)cyclopenten-1-yl]methylidene]hydroxylamine) as a small molecule antagonist for FGF23. Additional structure-activity studies developed a series of ZINC13407541 analogues with enhanced drug-like properties. In this study, we tested in a pre-clinical Hyp mouse homologue of XLH a direct connect analogue (8n) [(E)-2-(4-(tert-butyl)phenyl)cyclopent-1-ene-1-carbaldehyde oxime] that exhibited the greatest stability in microsomal assays, and 13a [(E)-2-((E)-4-methylstyryl)benzaldehyde oxime] that exhibited increased in vitro potency. Using cryo-electron microscopy (Cryo-EM) structure and computational docking, we identified a key binding residue (Q156) of the FGF23 antagonists, ZINC13407541 and its analogues (8n and 13a) in the N-terminal domain of FGF23 protein. Site-directed mutagenesis and bimolecular fluorescence complementation (BiFC)-fluorescence resonance energy transfer (FRET) assay confirmed the binding site of these three antagonists. We found that pharmacological inhibition of FGF23 with either of these compounds blocked FGF23 signaling and increased serum phosphate and 1,25(OH)2D concentrations in Hyp mice. Long-term parenteral treatment with 8n or 13a also enhanced linear bone growth, increased mineralization of bone, and narrowed the growth plate in Hyp mice. The more potent 13a compound had greater therapeutic effects in Hyp mice. Further optimization of these FGF23 inhibitors may lead to versatile drugs to treat excess FGF23-mediated disorders.

60 APPLIED LIFE SCIENCES↗

Low complexity domains of the nucleocapsid protein of SARS-CoV-2 form amyloid fibrils

The self-assembly of the Nucleocapsid protein (NCAP) of SARS-CoV-2 is crucial for its function. Computational analysis of the amino acid sequence of NCAP reveals low-complexity domains (LCDs) akin to LCDs in other proteins known to self-assemble as phase separation droplets and amyloid fibrils. Previous reports have described NCAP’s propensity to phase-separate. Here we show that the central LCD of NCAP is capable of both, phase separation and amyloid formation. Within this central LCD we identified three adhesive segments and determined the atomic structure of the fibrils formed by each. Those structures guided the design of G12, a peptide that interferes with the self-assembly of NCAP and demonstrates antiviral activity in SARS-CoV-2 infected cells. Our work, therefore, demonstrates the amyloid form of the central LCD of NCAP and suggests that amyloidogenic segments of NCAP could be targeted for drug development.

60 APPLIED LIFE SCIENCES↗

Complete inhibition of a polyol nucleation by a micromolar biopolymer additive

Preventing spontaneous crystallization of supersaturated solutions by additives is of critical interest to successful process design and implementation, with numerous applications in chemical, pharmaceutical, medical, pigment, and food industries, but challenges remain in laboratory and industry settings and fundamental understanding is lacking. When copresented with antifreeze proteins (AFPs), otherwise spontaneously crystallizing osmolytes are maintained at high supersaturations for months in over-wintering organisms. Thus, we here explore the inhibition phenomenon by AFPs, using persistent crystallization of a common sugar alcohol, D-mannitol, as a case study. We report experimentally that DAFP1, an insect AFP, completely inhibits D-mannitol nucleation. Computer simulations reveal a new mechanism for crystallization inhibition where the population of the crystal-forming conformers are selectively bound and randomized in solution by hydrogen bonding to the protein surface. These results highlight the advantages of using natural polymers to address crystallization inhibition challenges and suggest new strategies in controlling the nucleation processes.

59 BASIC BIOLOGICAL SCIENCES↗

TwoFold: Highly accurate structure and affinity prediction for protein-ligand complexes from sequences

We describe our development of ab initio protein-ligand binding pose prediction models based on transformers and binding affinity prediction models based on the neural tangent kernel (NTK). Folding both protein and ligand, the TwoFold models achieve efficient and quality predictions matching state-of-the-art implementations while additionally reconstructing protein structures. In conclusion, solving NTK models points to a new use case for highly optimized linear solver benchmarking codes on HPC.

60 APPLIED LIFE SCIENCES↗

In Vitro Evolution and Computational Approaches to Predict, Prevent and Control Future Pandemics

The natural tendency of virus to mutate and the under-sampling of the environment makes it difficult to become aware of the emergence of new viral strains with pandemic potential. Being able to predict what mutations make a virus more infective might allow to spot such strains with minimal sampling and potentially allow to predict/prevent the next pandemic. The team attempted to mimic natural viral mutations and recombination through computational and experimental methods producing a variety of mutants of a SARS-COV-2 protein (receptor binding domain, RBD, of spike protein) responsible for viral entry in mammalian cells. The library of mutants was then interrogated for ability and lack-there-of to interact with the host cell receptor ushering viral entry, Angiotensin-converting enzyme 2 (ACE2). The negative and positive data set is intended to “teach the rules” of virus-host receptor interaction. Additionally, the positive clones were used to screen a set of antibody mutants designed to widen the breadth of viral mutants recognition, to demonstrate that this kind of libraries could also be a tool to produce antibody therapeutics impervious to viral mutation, even before a pandemic strain is discovered.

59 BASIC BIOLOGICAL SCIENCES↗

Photoinduced chemomimetic biocatalysis for enantioselective intermolecular radical conjugate addition

Exploiting nature’s catalysts for non-natural transformations that are inaccessible to chemocatalysis is highly desirable but challenging. On the one hand, the widespread nicotinamide-dependent oxidoreductases have not been utilized for single-electron-transfer-induced bimolecular cross-couplings; on the other, the addition of catalytic asymmetric radical con- jugate to terminal alkenes remains a challenge owing to strong racemic background reaction and unselective termination of prochiral radical species. Here we report a chemomimetic biocatalysitic approach for construction of alpha-carbonyl stereocentres via an unnatural intermolecular conjugate addition of N-(acyloxy)phthalimides-derived radicals with acceptor-substituted terminal alkenes, by combination of visible-light excitation and nicotinamide-dependent ketoreductases (KREDs). Based on protein crystal structure, we engineered KREDs via a semi-rational mutagenesis strategy to improve reaction outcomes with a small and high-quality variants library. Here, mechanistic investigations combining wet experiments, crystallographic studies and computational simulations demonstrate that the repurposed biocatalyst can suppress racemic background reaction and unselected side reactions, yielding enantioselectivity that is challenging to achieve by chemocatalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computer display and manipulation of biological molecules

This paper describes a computer model that was designed to investigate the conformation of molecules, macromolecules and subsequent complexes. Utilizing an advanced 3-D dynamic computer display system, the model is sufficiently versatile to accommodate a large variety of molecular input and to generate data for multiple purposes such as visual representation of conformational changes, and calculation of conformation and interaction energy. Molecules can be built on the basis of several levels of information. These include the specification of atomic coordinates and connectivities and the grouping of building blocks and duplicated substructures using symmetry rules found in crystals and polymers such as proteins and nucleic acids. Called AIMS (Ames Interactive Molecular modeling System), the model is now being used to study pre-biotic molecular evolution toward life.

Coeckelenbergh, Y.↗

Combining Force Fields and Neural Networks for an Accurate Representation of Chemically Diverse Molecular Interactions

A key goal of molecular modeling is the accurate reproduction of the true quantum mechanical potential energy of arbitrary molecular ensembles with a tractable classical approximation. The challenges are that analytical expressions found in general purpose force fields struggle to faithfully represent the intermolecular quantum potential energy surface at close distances and in strong interaction regimes; that the more accurate neural network approximations do not capture crucial physics concepts, e.g., nonadditive inductive contributions and application of electric fields; and that the ultra-accurate narrowly targeted models have difficulty generalizing to the entire chemical space. We therefore designed a hybrid wide-coverage intermolecular interaction model consisting of an analytically polarizable force field combined with a short-range neural network correction for the total intermolecular interaction energy. Here, we describe the methodology and apply the model to accurately determine the properties of water, the free energy of solvation of neutral and charged molecules, and the binding free energy of ligands to proteins. The correction is subtyped for distinct chemical species to match the underlying force field, to segment and reduce the amount of quantum training data, and to increase accuracy and computational speed. For the systems considered, the hybrid ab initio parametrized Hamiltonian reproduces the two-body dimer quantum mechanics (QM) energies to within 0.03 kcal/mol and the nonadditive many-molecule contributions to within 2%. Simulations of molecular systems using this interaction model run at speeds of several nanoseconds per day.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mass balances for a biological life support system simulation model

Design decisions to aid the development of future space based biological life support systems (BLSS) can be made with simulation models. The biochemistry stoichiometry was developed for: (1) protein, carbohydrate, fat, fiber, and lignin production in the edible and inedible parts of plants; (2) food consumption and production of organic solids in urine, feces, and wash water by the humans; and (3) operation of the waste processor. Flux values for all components are derived for a steady state system with wheat as the sole food source. The large scale dynamics of a materially closed (BLSS) computer model is described in a companion paper. An extension of this methodology can explore multifood systems and more complex biochemical dynamics while maintaining whole system closure as a focus.

Volk, Tyler↗

Early Drug Discovery and Development of Novel Cancer Therapeutics Targeting DNA Polymerase Eta (POLH)

Polymerase eta (or Pol η or POLH) is a specialized DNA polymerase that is able to bypass certain blocking lesions, such as those generated by ultraviolet radiation (UVR) or cisplatin, and is deployed to replication foci for translesion synthesis as part of the DNA damage response (DDR). Inherited defects in the gene encoding POLH (a.k.a., XPV) are associated with the rare, sun-sensitive, cancer-prone disorder, xeroderma pigmentosum, owing to the enzyme’s ability to accurately bypass UVR-induced thymine dimers. In standard-of-care cancer therapies involving platinum-based clinical agents, e.g., cisplatin or oxaliplatin, POLH can bypass platinum-DNA adducts, negating benefits of the treatment and enabling drug resistance. POLH inhibition can sensitize cells to platinum-based chemotherapies, and the polymerase has also been implicated in resistance to nucleoside analogs, such as gemcitabine. POLH overexpression has been linked to the development of chemoresistance in several cancers, including lung, ovarian, and bladder. Co-inhibition of POLH and the ATR serine/threonine kinase, another DDR protein, causes synthetic lethality in a range of cancers, reinforcing that POLH is an emerging target for the development of novel oncology therapeutics. Using a fragment-based drug discovery approach in combination with an optimized crystallization screen, we have solved the first X-ray crystal structures of small novel drug-like compounds, i.e., fragments, bound to POLH, as starting points for the design of POLH inhibitors. The intrinsic molecular resolution afforded by the method can be quickly exploited in fragment growth and elaboration as well as analog scoping and scaffold hopping using medicinal and computational chemistry to advance hits to lead. An initial small round of medicinal chemistry has resulted in inhibitors with a range of functional activity in an in vitro biochemical assay, leading to the rapid identification of an inhibitor to advance to subsequent rounds of chemistry to generate a lead compound. Importantly, our chemical matter is different from the traditional nucleoside analog-based approaches for targeting DNA polymerases.

60 APPLIED LIFE SCIENCES↗

Meta-virus resource (MetaVR): expanding the frontiers of viral diversity with 24 million uncultivated virus genomes

Viruses are ubiquitous in all environments and impact host metabolism, evolution, and ecology, although our knowledge of their biodiversity is still extremely limited. Viral diversity from genomic and metagenomic datasets has led to an explosion of uncultivated virus genomes (UViGs) and the development of specialized databases to catalog this viral diversity, though many lack comprehensive integration. Here, we introduce meta-virus resource (MetaVR), the successor of the IMG/VR database, designed to overcome previous limitations such as large-scale querying and programmatic access. Drawing on the increase of publicly available genomes and metagenomes, MetaVR significantly expands viral diversity, now comprising 24,435,662 UViGs, a 57.6% increase from its predecessor, organized into over 12 million viral operational taxonomic units. Key enhancements include the integration of curated eukaryotic host information, the integration of protein clusters and predicted structures for comparative studies, and an API for programmatic data access. Furthermore, MetaVR features an updated taxonomic framework based on ICTV release 39, assignment to Baltimore classes, and enhanced host assignment through novel computational tools like iPHoP. These advancements position MetaVR as a unique resource for exploring viral diversity, evolution, and host interactions across diverse environments. MetaVR can be freely accessed at https://www.meta-virome.org/.

Fiamenghi, Mateus B↗

Coarse-Grained Modeling Elucidates Differential Metabolism of Saccharomyces cerevisiae under Varied Nutrient Limitations

Microorganisms such as Saccharomyces cerevisiae have a native ability to adapt their metabolism to varying nutrient conditions. Understanding their responses to nutrient limitations is critical for decoding cellular physiology and designing strategies for metabolic engineering. While the influence of carbon availability on yeast metabolism has been extensively studied, the role of nitrogen availability remains relatively underexplored. In this study, we utilized a coarse-grained kinetic model to systematically analyze and compare the effects of carbon and nitrogen limitations on yeast metabolism. Our model successfully revealed the differential metabolic characteristics of S. cerevisiae under carbon- and nitrogen-limited chemostat conditions. It also highlighted the significance of protein activity regulation at varying carbon-to-nitrogen ratios, and elucidated distinct strategies employed to maintain ATP homeostasis. Furthermore, this study provides a computational tool for investigating yeast physiology under nutrient limitations and offers quantitative and mechanistic insights into yeast metabolism.

59 BASIC BIOLOGICAL SCIENCES↗