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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

State of the Art of Molecular Visualization in Immersive Virtual Environments

Abstract Visualization plays a crucial role in molecular and structural biology. It has been successfully applied to a variety of tasks, including structural analysis and interactive drug design. While some of the challenges in this area can be overcome with more advanced visualization and interaction techniques, others are challenging primarily due to the limitations of the hardware devices used to interact with the visualized content. Consequently, visualization researchers are increasingly trying to take advantage of new technologies to facilitate the work of domain scientists. Some typical problems associated with classic 2D interfaces, such as regular desktop computers, are a lack of natural spatial understanding and interaction, and a limited field of view. These problems could be solved by immersive virtual environments and corresponding hardware, such as virtual reality head‐mounted displays. Thus, researchers are investigating the potential of immersive virtual environments in the field of molecular visualization. There is already a body of work ranging from educational approaches to protein visualization to applications for collaborative drug design. This review focuses on molecular visualization in immersive virtual environments as a whole, aiming to cover this area comprehensively. We divide the existing papers into different groups based on their application areas, and types of tasks performed. Furthermore, we also include a list of available software tools. We conclude the report with a discussion of potential future research on molecular visualization in immersive environments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

NEAR: Neural Embeddings for Amino acid Relationships

Protein language models (PLMs) have recently demonstrated potential to supplant classical protein database search methods based on sequence alignment, but are slower than common alignment-based tools and appear to be prone to a high rate of false labeling. Here, we present NEAR, a method based on neural representation learning that is designed to improve both speed and accuracy of search for likely homologs in a large protein sequence database. NEAR’s ResNet embedding model is trained using contrastive learning guided by trusted sequence alignments. It computes per-residue embeddings for target and query protein sequences, and identifies alignment candidates with a pipeline consisting of residue-level k-NN search and a simple neighbor aggregation scheme. Tests on a benchmark consisting of trusted remote homologs and randomly shuffled decoy sequences reveal that NEAR substantially improves accuracy relative to state-of-the-art PLMs, with lower memory requirements and faster embedding and search speed. While these results suggest that the NEAR model may be useful for standalone homology detection with increased sensitivity over standard alignment-based methods, in this manuscript we focus on a more straightforward analysis of the model’s value as a high-speed pre-filter for sensitive annotation. In that context, NEAR is at least 5x faster than the pre-filter currently used in the widely-used profile hidden Markov model (pHMM) search tool HMMER3, and also outperforms the pre-filter used in our fast pHMM tool, nail.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Chemical manipulation of an activation/inhibition switch in the nuclear receptor PXR

Nuclear receptors are ligand-activated transcription factors that can often be useful drug targets. Unfortunately, ligand promiscuity leads to two-thirds of receptors remaining clinically untargeted. PXR is a nuclear receptor that can be activated by diverse compounds to elevate metabolism, negatively impacting drug efficacy and safety. This presents a barrier to drug development because compounds designed to target other proteins must avoid PXR activation while retaining potency for the desired target. This problem could be avoided by using PXR antagonists, but these compounds are rare, and their molecular mechanisms remain unknown. Here, we report structurally related PXR-selective agonists and antagonists and their corresponding co-crystal structures to describe mechanisms of antagonism and selectivity. Structural and computational approaches show that antagonists induce PXR conformational changes incompatible with transcriptional coactivator recruitment. These results guide the design of compounds with predictable agonist/antagonist activities and bolster efforts to generate antagonists to prevent PXR activation interfering with other drugs.

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

QM Investigation of Rare Earth Ion Interactions with First Hydration Shell Waters and Protein-Based Coordination Models

Here, conventional methods for extracting rare earth metals (REMs) from mined mineral ores are inefficient, expensive, and environmentally damaging. Recent discovery of lanmodulin (LanM), a protein that coordinates REMs with high-affinity and selectivity over competing ions, provides inspiration for new REM refinement methods. Here, we used quantum mechanical (QM) methods to investigate trivalent lanthanide cation (Ln 3+ ) interactions with coordination systems representing bulk solvent water and protein binding sites. Energy decomposition analysis (EDA) showed differences in the energetic components of Ln 3+ interaction with representatives of solvent (water, H 2 O) and protein binding sites (acetate, CH 3 COO – ), highlighting the importance of accurate description of electrostatics and polarization in computational modeling of REM interactions with biological and bioinspired molecules. Relative binding free energies were obtained for Ln 3+ with coordination complexes originating from binding sites in PDB structures of a lanthanum binding peptide (PDB entry 7CCO) and LanM, with explicit consideration of the first hydration shell waters, according to quasi-chemical theory (QCT). Beyond the first shell, the bulk solvent environment was represented with an implicit continuum model. Ln 3+ interactions with (H 2 O) 9 and both binding site models became more favorable, moving down the periodic series. This trend was more pronounced with the protein binding site models than with water, resulting in affinity increasing with periodic number, except for the last REM, Lu 3+ , which bound less favorably than the preceding element, Yb 3+ . Using the truncated 7CCO binding site model, the magnitude and trend of the experimental Ln 3+ relative binding free energies for the whole 7CCO peptide were reproduced. Conversely, the previously reported experimental data for LanM show a preference for the earlier lanthanides; this is likely due to longer-range interactions and cooperative effects, which are not represented by the reduced models. Using the truncated 7CCO binding site model, the magnitude and trend of the experimental Ln 3+ relative binding free energies for the whole 7CCO peptide were reproduced. In contrast to the previously reported experimental data for LanM, the peptide preferentially binds the earlier lanthanides. This difference likely arises due to longer-range interactions and cooperative effects not represented by the peptide. Further investigation of Ln 3+ interactions with whole proteins using polarizable molecular mechanics models with explicit solvent is warranted to understand the influence of longer-ranged interactions, cooperativity, and bulk solvent. Nevertheless, the present work provides new insights into Ln 3+ interactions with biomolecules and presents an effective computational platform for designing specific single-site REM binding peptides more efficiently.

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

Early Career Researcher Support for Foundations of Molecular Modeling and Simulation Conference

This project enabled 14 early career researchers, including graduate students and postdocs, to attend the 8th triennial conference on Foundations of Molecular Modeling and Simulation (FOMMS 2022), which took place in Delavan, WI (Lake Geneva area), from July 17 to 21, 2022. Molecular-level modeling and simulation, which serve as critical tools for advancing science and engineering in a myriad of applications, were at the forefront of discussions. These ranged from the design of new materials for various energy applications to insights into the interaction between drug molecules and protein receptors within the body.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

De novo design of protein homodimers containing tunable symmetric protein pockets

Function follows form in biology, and the binding of small molecules requires proteins with pockets that match the shape of the ligand. For design of binding to symmetric ligands, protein homo-oligomers with matching symmetry are advantageous as each protein subunit can make identical interactions with the ligand. Here, we describe a general approach to designing hyperstable C2 symmetric proteins with pockets of diverse size and shape. We first designed repeat proteins that sample a continuum of curvatures but have low helical rise, then docked these into C2 symmetric homodimers to generate an extensive range of C2 symmetric cavities. We used this approach to design thousands of C2 symmetric homodimers, and characterized 101 of them experimentally. Of these, the geometry of 31 were confirmed by small angle X-ray scattering and 2 were shown by crystallographic analyses to be in close agreement with the computational design models. These scaffolds provide a rich set of starting points for binding a wide range of C2 symmetric compounds.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Design and evaluation of novel analogs of 2-amino-4-boronobutanoic acid (ABBA) as inhibitors of human gamma-glutamyl transpeptidase

Inhibitors of gamma-glutamyl transpeptidase (GGT1, aka gamma-glutamyl transferase) are needed for the treatment of cancer, cardiovascular illness and other diseases. Compounds that inhibit GGT1 have been evaluated in the clinic, but no inhibitor has successfully demonstrated specific and systemic GGT1 inhibition. All have severe side effects. L-2-amino-4-boronobutanoic acid (l-ABBA), a glutamate analog, is the most potent GGT1 inhibitor in vitro. In this study, we have solved the crystal structure of human GGT1 (hGGT1) with ABBA bound in the active site. The structure was interrogated to identify interactions between the enzyme and the inhibitor. Based on these data, a series of novel ABBA analogs were designed and synthesized. Their inhibitory activity against the hydrolysis and transpeptidation activities of hGGT1 were determined. The lead compounds were crystalized with hGGT1 and the structures solved. The kinetic data and structures of the complexes provide new insights into the critical role of protein structure dynamics in developing compounds for inhibition of hGGT1.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Artificial Intelligence for Autonomous Molecular Design: A Perspective

Domain-aware artificial intelligence has been increasingly adopted in recent years to expedite molecular design in various applications, including drug design and discovery. Recent advances in areas such as physics-informed machine learning and reasoning, software engineering, high-end hardware development, and computing infrastructures are providing opportunities to build scalable and explainable AI molecular discovery systems. This could improve a design hypothesis through feedback analysis, data integration that can provide a basis for the introduction of end-to-end automation for compound discovery and optimization, and enable more intelligent searches of chemical space. Several state-of-the-art ML architectures are predominantly and independently used for predicting the properties of small molecules, their high throughput synthesis, and screening, iteratively identifying and optimizing lead therapeutic candidates. However, such deep learning and ML approaches also raise considerable conceptual, technical, scalability, and end-to-end error quantification challenges, as well as skepticism about the current AI hype to build automated tools. To this end, synergistically and intelligently using these individual components along with robust quantum physics-based molecular representation and data generation tools in a closed-loop holds enormous promise for accelerated therapeutic design to critically analyze the opportunities and challenges for their more widespread application. This article aims to identify the most recent technology and breakthrough achieved by each of the components and discusses how such autonomous AI and ML workflows can be integrated to radically accelerate the protein target or disease model-based probe design that can be iteratively validated experimentally. Taken together, this could significantly reduce the timeline for end-to-end therapeutic discovery and optimization upon the arrival of any novel zoonotic transmission event. Our article serves as a guide for medicinal, computational chemistry and biology, analytical chemistry, and the ML community to practice autonomous molecular design in precision medicine and drug discovery.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

The SARS-CoV-2 Spike Variant D614G Favors an Open Conformational State

The COVID-19 pandemic underwent a rapid transition with the emergence of a dominant viral variant (from the “D-form” to the “G-form”) that carried an amino acid substitution D614G in its “Spike” protein. The G-form is more infectious in vitro and associated with increased viral loads. To gain insight into the molecular-level underpinnings of these characteristics, we employed microsecond all-atom simulations. Here we show that changes in the protein energetics favor a higher population of infection-capable (open) states through release of hydrogen bonds of an asymmetry present in the D-form but not the G-form. Thus, the increased infectivity of the G-form is likely due to a higher rate of profitable binding encounters with the host receptor. It is also predicted to be more neutralization sensitive due to enhanced exposure of the receptor binding domain, a key target region for neutralizing antibodies. These results are significant for vaccine design. The Molecular Dynamics Simulations datasets generated in this study, namely that of the soluble form SARS-CoV-2 Spike protein are made available here. Details of the trajectory and supporting files are given in the README.txt file.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Rational Design of Novel Biomimetic Sequence-Defined Polymers for Mineralization Applications

Silica biomineralization is a naturally occurring process, wherein organisms use proteins and other biological structures to direct the formation of complex, hierarchical nanostructures. Discovery and characterization of such proteins and their underlying mechanisms spurred significant efforts to identify routes for biomimetic mineralization that reproduce the exquisite shapes and size selectivities found in nature. A common strategy has been the use of short peptide sequences with chemistry mimicking those found in natural systems, such as the use of the silaffin-derived R5 peptide. While progress has been made using this approach, there are many limitations that have prevented breakthroughs in biomimicry. To advance our ability to use charged macromolecules for silica formation, we propose to use sequence-defined synthetic polymers known as peptoids, or N-substituted polyglycines, which present significant capability for the precise tuning of sequence and structure beyond what can often be achieved with peptides alone. This study presents a computationally predicted design of these polymers that leads to the controlled formation of silica nanomaterials. We investigate surface adsorption and the mineralization process through analysis of binding mechanisms and energetics of the R5 system. Next, we synthesized two R5-inspired peptoids and validated our prediction in the design of mineralization polymers through characterization using surface plasmon resonance and electron microscopy. Here, this computationally guided study holds great promise for designing new sequences with unprecedented control of the placement of chemical functional groups, thus allowing for further unraveling of silicification mechanisms and the eventual design of sequence-defined synthetic polymers leading to the predictive synthesis of nanostructured functional materials.

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

Dual Enhancement of Thermostability and Activity of Xylanase through Computer-Aided Rational Design

In the realm of enzyme engineering, the dual enhancement of thermostability and activity remains a challenge. Herein, we employed a computer-aided approach integrating folding free energy calculations and evolutionary analysis to engineer Paecilomyces thermophila xylanase into a hyperthermophilic enzyme for application in the paper and pulp industry. Through the computational rational design, XynM9 with superior thermostability and enhanced activity was designed. Its optimal reaction temperature increases by 10 °C to 85 °C, its T m increases by 10 °C to 93 °C, and its half-life increases 11-fold to 5.8 h. Additionally, its catalytic efficiency improves by 57% to 3926 s –1 mM –1 . Molecular dynamics simulations revealed that XynM9 is stabilized by more hydrogen bonds and salt bridges than wild-type xylanase. The mutant’s narrower catalytic cleft enhances the substrate-binding affinity, thus improving the catalytic efficiency. In harsh conditions at 80 °C and pH 10, using XynM9 significantly reduced both hemicellulose and lignin, which makes it a good candidate for use in the paper and pulp process. Finally, our study presents an accurate and efficient strategy for the dual enhancement of enzyme properties, guiding further improvement of computational tools for protein stabilization.

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

Architector 2.0: Expanded Capabilities for Metal Complex Engineering

Automated three-dimensional molecular construction from two-dimensional graph representations is critical to high-throughput discovery eIorts. Software capabilities in this area have accelerated research across fields ranging from protein design and drug discovery to transition metal catalyst development. When Architector was first introduced, it uniquely enabled high-throughput, chemically relevant three-dimensional construction of f-element complexes. Since its introduction, Architector has been applied in large-scale computational campaigns, targeted studies in critical mineral extraction, and artificial intelligence-driven discovery eIorts.

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

HDBind: encoding of molecular structure with hyperdimensional binary representations

Traditional methods for identifying “hit” molecules from a large collection of potential drug-like candidates rely on biophysical theory to compute approximations to the Gibbs free energy of the binding interaction between the drug and its protein target. These approaches have a significant limitation in that they require exceptional computing capabilities for even relatively small collections of molecules. Increasingly large and complex state-of-the-art deep learning approaches have gained popularity with the promise to improve the productivity of drug design, notorious for its numerous failures. However, as deep learning models increase in their size and complexity, their acceleration at the hardware level becomes more challenging. Hyperdimensional Computing (HDC) has recently gained attention in the computer hardware community due to its algorithmic simplicity relative to deep learning approaches. The HDC learning paradigm, which represents data with high-dimension binary vectors, allows the use of low-precision binary vector arithmetic to create models of the data that can be learned without the need for the gradient-based optimization required in many conventional machine learning and deep learning methods. This algorithmic simplicity allows for acceleration in hardware that has been previously demonstrated in a range of application areas (computer vision, bioinformatics, mass spectrometery, remote sensing, edge devices, etc.). To the best of our knowledge, our work is the first to consider HDC for the task of fast and efficient screening of modern drug-like compound libraries. We also propose the first HDC graph-based encoding methods for molecular data, demonstrating consistent and substantial improvement over previous work. We compare our approaches to alternative approaches on the well-studied MoleculeNet dataset and the recently proposed LIT-PCBA dataset derived from high quality PubChem assays. We demonstrate our methods on multiple target hardware platforms, including Graphics Processing Units (GPUs) and Field Programmable Gate Arrays (FPGAs), showing at least an order of magnitude improvement in energy efficiency versus even our smallest neural network baseline model with a single hidden layer. Our work thus motivates further investigation into molecular representation learning to develop ultra-efficient pre-screening tools. We make our code publicly available at https://github.com/LLNL/hdbind.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Structure- and Interaction-Based Design of Anti-SARS-CoV-2 Aptamers

Aptamer selection against novel infections is a complicated and time-consuming approach. Synergy can be achieved by using computational methods together with experimental procedures. In this study, we aim to develop a reliable methodology for a rational aptamer in silico et vitro design. The new approach combines multiple steps: (1) Molecular design, based on screening in a DNA aptamer library and directed mutagenesis to fit the protein tertiary structure; (2) 3D molecular modeling of the target; (3) Molecular docking of an aptamer with the protein; (4) Molecular dynamics (MD) simulations of the complexes; (5) Quantum-mechanical (QM) evaluation of the interactions between aptamer and target with further analysis; (6) Experimental verification at each cycle for structure and binding affinity by using small-angle X-ray scattering, cytometry, and fluorescence polarization. By using a new iterative design procedure, structure- and interaction-based drug design (SIBDD), a highly specific aptamer to the receptor-binding domain of the SARS-CoV-2 spike protein, was developed and validated. The SIBDD approach enhances speed of the high-affinity aptamers development from scratch, using a target protein structure. The method could be used to improve existing aptamers for stronger binding. This approach brings to an advanced level the development of novel affinity probes, functional nucleic acids. It offers a blueprint for the straightforward design of targeting molecules for new pathogen agents and emerging variants.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Quantitative comparison between sub-millisecond time resolution single-molecule FRET measurements and 10-second molecular simulations of a biosensor protein

Molecular Dynamics (MD) simulations seek to provide atomic-level insights into conformationally dynamic biological systems at experimentally relevant time resolutions, such as those afforded by single-molecule fluorescence measurements. However, limitations in the time scales of MD simulations and the time resolution of single-molecule measurements have challenged efforts to obtain overlapping temporal regimes required for close quantitative comparisons. Achieving such overlap has the potential to provide novel theories, hypotheses, and interpretations that can inform idealized experimental designs that maximize the detection of the desired reaction coordinate. Here, we report MD simulations at time scales overlapping with in vitro single-molecule Förster (fluorescence) resonance energy transfer (smFRET) measurements of the amino acid binding protein LIV-BP SS at sub-millisecond resolution. Computationally efficient all-atom structure-based simulations, calibrated against explicit solvent simulations, were employed for sampling multiple cycles of LIV-BP SS clamshell-like conformational changes on the time scale of seconds, examining the relationship between these events and those observed by smFRET. The MD simulations agree with the smFRET measurements and provide valuable information on local dynamics of fluorophores at their sites of attachment on LIV-BP SS and the correlations between fluorophore motions and large-scale conformational changes between LIV-BP SS domains. We further utilize the MD simulations to inform the interpretation of smFRET data, including Förster radius (R 0 ) and fluorophore orientation factor (κ 2 ) determinations. The approach we describe can be readily extended to distinct biochemical systems, allowing for the interpretation of any FRET system conjugated to protein or ribonucleoprotein complexes, including those with more conformational processes, as well as those implementing multi-color smFRET.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Ligand- and Structure-Based Analysis of Deep Learning-Generated Potential α2a Adrenoceptor Agonists

The α2a adrenoceptor is a particularly medically relevant subtype of the G protein-couple receptor family. Unfortunately, high-throughput techniques aimed at producing novel drug leads targeted to this receptor have been largely unsuccessful due to the complex pharmacology of adrenergic receptors. As such, cutting-edge in silico ligand- and structure-based assessment and de novo deep learning methods are well-positioned to provide new insights into protein-ligand interactions and potential active compounds. In this work, we (i) collect a dataset of α2a adrenoceptor agonists and provide it as a resource for the drug design community; (ii) use the dataset as a basis to generate candidate active structures via deep learning; and (iii) apply computational ligand- and structure-based analysis techniques to gain new insight into α2a adrenoceptor agonists and assess the quality of the computer-generated compounds. Additionally, we further describe how such assessment techniques can be applied to putative chemical probes with a case study involved proposed medetomidine-based probes.

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

Structure-Based Design of Small-Molecule Inhibitors of Human Interleukin-6

Human Interleukin-6 (hIL-6) is a pro inflammatory cytokine that binds to its receptor, IL-6Rα followed by binding to gp130 and subsequent dimerization to form a hexamer signaling complex. As a critical inflammation mediator, hIL-6 is associated with a diverse range of diseases and monoclonal antibodies in clinical use that either target IL-6Rα or hIL-6 to inhibit signaling. Here, we perform high-throughput structure-based computational screening using ensemble docking for small-molecule antagonists for which the target conformations were taken from 600 ns long molecular dynamics simulations of the apo protein. Prior knowledge of the contact sites from binary complex studies and experimental work was incorporated into the docking studies. The top 20 scoring ligands from the in silico studies after post analysis were subjected to in vitro functional assays. Among these compounds, the ligand with the second-highest calculated binding affinity experimentally showed an ~84% inhibitory effect on IL6-induced STAT3 reporter activity at 10 μM concentration. This finding may pave the way for designing small-molecule inhibitors of hIL-6 of therapeutic significance.

Human Interleukin-6↗
DOE OSTI2024

A goldilocks computational protocol for inhibitor discovery targeting DNA damage responses including replication-repair functions

While many researchers can design knockdown and knockout methodologies to remove a gene product, this is mainly untrue for new chemical inhibitor designs that empower multifunctional DNA Damage Response (DDR) networks. Here, we present a robust Goldilocks (GL) computational discovery protocol to efficiently innovate inhibitor tools and preclinical drug candidates for cellular and structural biologists without requiring extensive virtual screen (VS) and chemical synthesis expertise. By computationally targeting DDR replication and repair proteins, we exemplify the identification of DDR target sites and compounds to probe cancer biology. Our GL pipeline integrates experimental and predicted structures to efficiently discover leads, allowing early-structure and early-testing (ESET) experiments by many laboratories. By employing an efficient VS protocol to examine protein-protein interfaces (PPIs) and allosteric interactions, we identify ligand binding sites beyond active sites, leveraging in silico advances for molecular docking and modeling to screen PPIs and multiple targets. A diverse 3,174 compound ESET library combines Diamond Light Source DSI-poised, Protein Data Bank fragments, and FDA-approved drugs to span relevant chemotypes and facilitate downstream hit evaluation efficiency for academic laboratories. Two VS per library and multiple ranked ligand binding poses enable target testing for several DDR targets. This GL library and protocol can thus strategically probe multiple DDR network targets and identify readily available compounds for early structural and activity testing to overcome bottlenecks that can limit timely breakthrough drug discoveries. By testing accessible compounds to dissect multi-functional DDRs and suggesting inhibitor mechanisms from initial docking, the GL approach may enable more groups to help accelerate discovery, suggest new sites and compounds for challenging targets including emerging biothreats and advance cancer biology for future precision medicine clinical trials.

59 BASIC BIOLOGICAL SCIENCES↗