Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular docking”

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 OSTI2024

CryoEM structures of the human CLC-2 voltage-gated chloride channel reveal a ball-and-chain gating mechanism

CLC-2 is a voltage-gated chloride channel that contributes to electrical excitability and ion homeostasis in many different tissues. Among the nine mammalian CLC homologs, CLC-2 is uniquely activated by hyperpolarization, rather than depolarization, of the plasma membrane. The molecular basis for the divergence in polarity of voltage gating among closely related homologs has been a long-standing mystery, in part because few CLC channel structures are available. Here, we report cryoEM structures of human CLC-2 at 2.46 – 2.76 Å, in the presence and absence of the selective inhibitor AK-42. AK-42 binds within the extracellular entryway of the Cl – -permeation pathway, occupying a pocket previously proposed through computational docking studies. In the apo structure, we observed two distinct conformations involving rotation of one of the cytoplasmic C-terminal domains (CTDs). In the absence of CTD rotation, an intracellular N-terminal 15-residue hairpin peptide nestles against the TM domain to physically occlude the Cl – -permeation pathway. This peptide is highly conserved among species variants of CLC-2 but is not present in other CLC homologs. Previous studies suggested that the N-terminal domain of CLC-2 influences channel properties via a “ball-and-chain” gating mechanism, but conflicting data cast doubt on such a mechanism, and thus the structure of the N-terminal domain and its interaction with the channel has been uncertain. Through electrophysiological studies of an N-terminal deletion mutant lacking the 15-residue hairpin peptide, we support a model in which the N-terminal hairpin of CLC-2 stabilizes a closed state of the channel by blocking the cytoplasmic Cl – -permeation pathway.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Structural basis of Gabija anti-phage defence and viral immune evasion

Bacteria encode hundreds of diverse defence systems that protect them from viral infection and inhibit phage propagation. Gabija is one of the most prevalent anti-phage defence systems, occurring in more than 15% of all sequenced bacterial and archaeal genomes, but the molecular basis of how Gabija defends cells from viral infection remains poorly understood. Here we use X-ray crystallography and cryo-electron microscopy (cryo-EM) to define how Gabija proteins assemble into a supramolecular complex of around 500 kDa that degrades phage DNA. Gabija protein A (GajA) is a DNA endonuclease that tetramerizes to form the core of the anti-phage defence complex. Two sets of Gabija protein B (GajB) dimers dock at opposite sides of the complex and create a 4:4 GajA–GajB assembly (hereafter, GajAB) that is essential for phage resistance in vivo. We show that a phage-encoded protein, Gabija anti-defence 1 (Gad1), directly binds to the Gabija GajAB complex and inactivates defence. A cryo-EM structure of the virally inhibited state shows that Gad1 forms an octameric web that encases the GajAB complex and inhibits DNA recognition and cleavage. Our results reveal the structural basis of assembly of the Gabija anti-phage defence complex and define a unique mechanism of viral immune evasion.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

30S subunit recognition and G1405 modification by the aminoglycoside-resistance 16S ribosomal RNA methyltransferase RmtC

Acquired ribosomal RNA (rRNA) methylation has emerged as a significant mechanism of aminoglycoside resistance in pathogenic bacterial infections. Modification of a single nucleotide in the ribosome decoding center by the aminoglycoside-resistance 16S rRNA (m 7 G1405) methyltransferases effectively blocks the action of all 4,6-deoxystreptamine ring-containing aminoglycosides, including the latest generation of drugs. To define the molecular basis of 30S subunit recognition and G1405 modification by these enzymes, we used a S-adenosyl-L-methionine analog to trap the complex in a postcatalytic state to enable determination of a global 3.0 Å cryo-electron microscopy structure of the m 7 G1405 methyltransferase RmtC bound to the mature Escherichia coli 30S ribosomal subunit. This structure, together with functional analyses of RmtC variants, identifies the RmtC N-terminal domain as critical for recognition and docking of the enzyme on a conserved 16S rRNA tertiary surface adjacent to G1405 in 16S rRNA helix 44 (h44). To access the G1405 N7 position for modification, a collection of residues across one surface of RmtC, including a loop that undergoes a disorder-to order transition upon 30S subunit binding, induces significant distortion of h44. This distortion flips G1405 into the enzyme active site where it is positioned for modification by two almost universally conserved RmtC residues. Further, these studies expand our understanding of ribosome recognition by rRNA modification enzymes and present a more complete structural basis for future development of strategies to inhibit m 7 G1405 modification to resensitize bacterial pathogens to aminoglycosides.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Molecular chaperone RAP interacts with LRP1 in a dynamic bivalent mode and enhances folding of ligand-binding regions of other LDLR family receptors

The low-density lipoprotein receptor (LDLR) family of receptors are cell-surface receptors that internalize numerous ligands and play crucial role in various processes, such as lipoprotein metabolism, hemostasis, fetal development, etc. Previously, receptor-associated protein (RAP) was described as a molecular chaperone for LDLR-related protein 1 (LRP1), a prominent member of the LDLR family. We aimed to verify this role of RAP for LRP1 and two other LDLR family receptors, LDLR and vLDLR, and to investigate the mechanisms of respective interactions using a cell culture model system, purified system, and in silico modelling. Upon coexpression of RAP with clusters of the ligand-binding complement repeats (CRs) of the receptors in secreted form in insect cells culture, the isolated proteins had increased yield, enhanced folding, and improved binding properties compared with proteins expressed without RAP, as determined by circular dichroism and surface plasmon resonance. Within LRP1 CR-clusters II and IV, we identified multiple sites comprised of adjacent CR doublets, which provide alternative bivalent binding combinations with specific pairs of lysines on RAP. Mutational analysis of these lysines within each of isolated RAP D1/D2 and D3 domains having high affinity to LRP1 and of conserved tryptophans on selected CR-doublets of LRP1, as well as in silico docking of a model LRP1 CR-triplet with RAP, indicated a universal role for these residues in interaction of RAP and LRP1. Consequently, we propose a new model of RAP interaction with LDLR family receptors based on switching of the bivalent contacts between molecules over time in a dynamic mode.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Molecular basis and biological relevance of bacterial and plant pinoresinol/lariciresinol reductase specificities

A bacterial pinoresinol/lariciresinol reductase (PLR) homolog named NrPinZ was obtained from a Novosphingobium rhizosphaerae sp. LY bacterial strain, with NrPinZ being part of its 5-step biochemical system catabolizing pinoresinol into coniferyl aldehyde and vanillin. Recombinant NrPinZ reduces racemic 8–8′ furanofuran lignans [(±)-pinoresinols, medioresinols, and syringaresinols] with similar overall catalytic efficiencies. In those reductions, only one of the two furan ring systems is reduced. Two other bacterial PLR homologs, NaPinZ and SlPinZ, from N. aromaticivorans F199 and Sphingobium lignivorans SYK-6, respectively, had comparable substrate versatilities and catalytic efficacies. Plant PLR homologs, by comparison, are either enantiospecific, enantioselective, or variants thereof, being able to reduce either one or both furan rings. For example, a recombinant enantioselective PLR (PLR_Tp2) from western red cedar (Thuja plicata) preferentially reduces both (+)-pinoresinol furan rings to afford (−)-secoisolariciresinol. BoltZ-2 modeling of NrPinZ and PLR_Tp2, together with substrate docking of (+)- and (−)-pinoresinols, medioresinols, and syringaresinols, was very instructive. The NrPinZ active site P1/P2 sub-pockets allow for both racemic forms to be catabolized. Conversely, the smaller P1 pocket in PLR_Tp2 preferentially positions (+)-pinoresinol for downstream metabolism into (−)-secoisolariciresinol, thereby providing a biochemical explanation for the different stereochemical outcomes. NrPinZ, NaPinZ, and SlPinZ, catalyzing substrate versatile catabolism of both racemic forms, may have important ramifications for gymnosperm and angiosperm lignin and lignan biodegradation, including its evolutionary significance and potential in enzyme engineering.

Boltz-2 molecular modeling↗
DOE OSTI2023

Fast and versatile sequence-independent protein docking for nanomaterials design using RPXDock

Computationally designed multi-subunit assemblies have shown considerable promise for a variety of applications, including a new generation of potent vaccines. One of the major routes to such materials is rigid body sequence-independent docking of cyclic oligomers into architectures with point group or lattice symmetries. Current methods for docking and designing such assemblies are tailored to specific classes of symmetry and are difficult to modify for novel applications. Here we describe RPXDock, a fast, flexible, and modular software package for sequence-independent rigid-body protein docking across a wide range of symmetric architectures that is easily customizable for further development. RPXDock uses an efficient hierarchical search and a residue-pair transform ( RPX ) scoring method to rapidly search through multidimensional docking space. We describe the structure of the software, provide practical guidelines for its use, and describe the available functionalities including a variety of score functions and filtering tools that can be used to guide and refine docking results towards desired configurations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Hallucination of closed repeat proteins containing central pockets

In pseudocyclic proteins, such as TIM barrels, β barrels, and some helical transmembrane channels, a single subunit is repeated in a cyclic pattern, giving rise to a central cavity that can serve as a pocket for ligand binding or enzymatic activity. Inspired by these proteins, we devised a deep-learning-based approach to broadly exploring the space of closed repeat proteins starting from only a specification of the repeat number and length. Biophysical data for 38 structurally diverse pseudocyclic designs produced in Escherichia coli are consistent with the design models, and the three crystal structures we were able to obtain are very close to the designed structures. Docking studies suggest the diversity of folds and central pockets provide effective starting points for designing small-molecule binders and enzymes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Fe protein docking transduces conformational changes to MoFe nitrogenase active site in a nucleotide-dependent manner

The reduction of dinitrogen to ammonia catalyzed by nitrogenase involves a complex series of events, including ATP hydrolysis, electron transfer, and activation of metal clusters for N 2 reduction. Early evidence shows that an essential part of the mechanism involves transducing information between the nitrogenase component proteins through conformational dynamics. Here, millisecond time-resolved hydrogen-deuterium exchange mass spectrometry was used to unravel peptide-level protein motion on the time scale of catalysis of Mo-dependent nitrogenase from Azotobacter vinelandii. Normal mode analysis calculations complemented this data, providing insights into the specific signal transduction pathways that relay information across protein interfaces at distances spanning 100 Å. Together, these results show that conformational changes induced by protein docking are rapidly transduced to the active site, suggesting a specific mechanism for activating the metal cofactor in the enzyme active site.

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

Molecular mechanism of phosphopeptide neoantigen immunogenicity

Abstract often leads to surface presentation of phosphopeptide neoantigens. However, their role in cancer immunogenicity remains unclear. Here we describe a mechanism by which an HLA-B*0702-specific acute myeloid leukemia phosphoneoantigen, pMLL 747–755 (EPR(pS)PSHSM), is recognized by a cognate T cell receptor named TCR27, a candidate for cancer immunotherapy. We show that the replacement of phosphoserine P 4 with serine or phosphomimetics does not affect pMHC conformation or peptide-MHC affinity but abrogates TCR27-dependent T cell activation and weakens binding between TCR27 and pMHC. Here we describe the crystal structures for TCR27 and cognate pMHC, map of the interface produced by nuclear magnetic resonance, and a ternary complex generated using information-driven protein docking. Our data show that non-covalent interactions between the epitope phosphate group and TCR27 are crucial for TCR specificity. This study supports development of new treatment options for cancer patients through target expansion and TCR optimization.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Unconventional low temperature decomposition of a saturated hydrocarbon over atomically-dispersed titanium-aluminum-boron catalyst

Sonochemically-synthesized atomically-dispersed titanium-aluminum-boron nanopowder (TiAlB NP) exhibits a remarkable low-temperature catalytic activation of aliphatic C-H bonds at 750 K followed by C-C bond activation thus emerging as a potent low-cost alternative to expensive platinum group metals. Here, the model saturated hydrocarbon, exo-tetrahydrodicyclopentadiene (C10H16), undergoes catalytic decomposition on TiAlB NPs in a chemical microreactor to produce 1,3-cyclopentadiene (c-C5H6), cyclopentene (c-C5H8), and molecular hydrogen (H2) as detected in situ via isomer-selective, single-photon ionization time-of-flight mass spectrometry. Extensive electronic structure theory calculations on model clusters of the catalyst decode a unique synergy among the atomic constituents of the catalyst and chemical bonding in this stepwise, retro Diels Alder reaction: Ti, although insensitive to C-H activation in its metallic state, initiates the catalysis via chemisorption of the hydrocarbon, adjacent B centers readily abstract hydrogen atoms and store them during the catalytic cycle, while Al stabilizes the catalyst structure yet providing space for critical docking sites for the departing hydrocarbons.

Biswas, Souvick↗
DOE OSTI2020

Ligand-Binding-Site Refinement to Generate Reliable Holo Protein Structure Conformations from Apo Structures

The first important step in a structure-based virtual screening is the judicious selection of a receptor protein. In cases where the holo protein receptor structure is unavailable, significant reduction in virtual screening performance has been reported. In this work, we present a robust method to generate reliable holo protein structure conformations from apo structures using molecular dynamics (MD) simulation with restraints derived from holo structure binding-site templates. We perform benchmark tests on two different datasets: 40 structures from a directory of useful decoy-enhanced (DUD-E) and 84 structures from the Gunasekaran dataset. Our results show successful refinement of apo binding-site structures toward holo conformations in 82% of the test cases. In addition, virtual screening performance of 40 DUD-E structures is significantly improved using our MD-refined structures as receptors with an average enrichment factor (EF), an EF1% value of 6.2 compared to apo structures with 3.5. Docking of native ligands to the refined structures shows an average ligand root mean square deviation (RMSD) of 1.97 Å (DUD-E dataset and Gunasekaran dataset) relative to ligands in the holo crystal structures, which is comparable to the self-docking (i.e., docking of the native ligand back to its crystal structure receptor) average, 1.34 Å (DUD-E dataset) and 1.36 Å (Gunasekaran dataset). On the other hand, docking to the apo structures yields an average ligand RMSD of 3.65 Å (DUD-E) and 2.90 Å (Gunasekaran). These results indicate that our method is robust and can be useful to improve virtual screening performance of apo structures.

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

Identification and characterization of the low molecular mass ferredoxins involved in central metabolism in Heliomicrobium modesticaldum

The homodimeric Type I reaction center (RC) from Heliomicrobium modesticaldum lacks the PsaC subunit found in Photosystem I and instead uses the interpolypeptide [4Fe–4S] cluster F X as the terminal electron acceptor. Our goal was to identify which of the small mobile dicluster ferredoxins encoded by the H. modesticaldum genome are capable of accepting electrons from the heliobacterial RC (HbRC) and pyruvate:ferredoxin oxidoreductase (PFOR), a key metabolic enzyme. Analysis of the genome revealed seven candidates: HM1_1462 (PshB1), HM1_1461 (PshB2), HM1_2505 (Fdx3), HM1_0869 (FdxB), HM1_1043, HM1_0357, and HM1_2767. Heterologous expression in Escherichia coli and studies using time-resolved optical spectroscopy revealed that only PshB1, PshB2, and Fdx3 are capable of accepting electrons from the HbRC and PFOR. Modeling studies using AlphaFold show that only PshB1, PshB2, and Fdx3 should be capable of docking on PFOR at a positively charged patch that overlays a surface-proximal [4Fe–4S] cluster. Proteomic analysis of wild-type and gene deletion strains ΔpshB1, ΔpshB2, ΔpshB1pshB2, and Δfdx3 grown under nitrogen-replete conditions revealed that Fdx3 is undetectable in the wild-type, ΔpshB1, and Δfdx3 strains, but it is present in the ΔpshB2 and ΔpshB1pshB2 strains, implying that Fdx3 may substitute for PshB2. When grown under nitrogen-deplete conditions, Fdx3 is present in the wild-type and all deletion strains except for Δfdx3. None of the knockout strains demonstrated significant impairment during chemotrophic dark growth on pyruvate, photoheterotrophic light growth on pyruvate, or phototrophic growth on acetate+CO 2 , indicating a high degree of redundancy among these three electron transfer proteins. As a result, loss of both PshB1 and PshB2, but not FdxB, resulted in poor growth under N 2 -fixing conditions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Pre–T cell receptors topologically sample self-ligands during thymocyte β-selection

Self-discrimination, a critical but ill-defined molecular process programmed during thymocyte development, requires myriad pre–T cell receptors (preTCRs) and αβTCRs. Using x-ray crystallography, we show how a preTCR applies the concave β-sheet surface of its single variable domain (Vβ) to “horizontally” grab the protruding MHC α2-helix. By contrast, αβTCRs purpose all six complementarity-determining region (CDR) loops of their paired VαVβ module to recognize peptides bound to major histocompatibility complex molecules (pMHCs) in “vertical” head-to-head binding. Additionally, the preTCR topological fit ensures that CDR3β reaches the peptide’s featured C-terminal segment for pMHC sampling, establishing the subsequent αβTCR canonical docking mode. “Horizontal” docking precludes germline CDR1β- and CDR2β-MHC binding to broaden β-chain repertoire diversification before αβTCR-mediated selection refinement. Thus, one subunit successively attunes the recognition logic of related multicomponent receptors.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

3D-Scaffold: A Deep Learning Framework to Generate 3D Coordinates of Drug-like Molecules with Desired Scaffolds

The prerequisite of therapeutic drug design is to identify novel molecules with desired biophysical and biochemical properties. Deep generative models have demonstrated their ability to find such molecules by exploring a huge chemical space efficiently. An effective way to obtain molecules with desired target properties is the preservation of critical scaffolds in the generation process. To this end, we propose a domain aware generative framework called 3D-Scaffold that takes 3D coordinates of a desired scaffold as an input and generates 3D coordinates of novel therapeutic candidates as an output while always preserving the desired scaffolds in generated structures. We show that our framework generates predominantly valid, unique, novel, and experimentally synthesizable molecules that have drug-like properties similar to the molecules in the training set. Using domain specific datasets, we generate covalent and non-covalent antiviral inhibitors. Therefore, to measure the success of our framework in generating therapeutic candidates, generated structures were subjected to high throughput virtual screening via docking simulations, which shows favorable interaction against SARS-CoV-2 main protease and non-structural protein endoribonuclease (NSP15) targets. Most importantly, our model performs well with relatively small volumes of training data and generalizes to new scaffolds, making it applicable to other domain.

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

Machine Learning-driven Molecular Design for Therapeutic Discovery

The ongoing novel coronavirus pandemic (COVID-19) has highlighted the need for new therapeutics to counter the threat of emerging viral pathogens. The main proteases are a promising target for developing antiviral inhibitors. In this work, we utilized a novel combination of artificial intelligence-driven iterative design of covalent inhibitor candidates, physics-based computational modeling of protein-inhibitor interactions, and “All in One” Native MS biophysical assay screening and characterization of therapeutic candidates. With our existing expertise in hit generation using a particular scaffold as a starting point, we first generated tens of thousands of compounds that preserve the key scaffold. In order to optimize the candidates, we calculated about 136 descriptors consisting of 2D and 3D features for molecules targeting the SARS-CoV-2 Main protease (Mpro). These compounds were initially filtered according to properties and further sorted by predicted binding affinity using our automated docking modeling and machine learning methods. We tested a handful of candidates and identified two as inhibitors of Mpro with micromolar affinities.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Crystal structure of domain of unknown function 507 (DUF507) reveals a new protein fold

The crystal structure of the domain of unknown function family 507 protein from Aquifex aeolicus is reported (AaDUF507, UniProt O67633, 183 residues). The structure was determined in two space groups (C222 1 and P3 2 21) at 1.9 Å resolution. The phase problem was solved by molecular replacement using an AlphaFold model as the search model. AaDUF507 is a Y-shaped α-helical protein consisting of an anti-parallel 4-helix bundle base and two helical arms that extend 30-Å from the base. The two crystal structures differ by a 25° rigid body rotation of the C-terminal arm. The tertiary structure exhibits pseudo-twofold symmetry. The structural symmetry mirrors internal sequence similarity: residues 11–57 and 102–148 are 30% identical and 53% similar with an E-value of 0.002. In one of the structures, electron density for an unknown ligand, consistent with nicotinamide or similar molecule, may indicate a functional site. Docking calculations suggest potential ligand binding hot spots in the region between the helical arms. Structure-based query of the Protein Data Bank revealed no other protein with a similar tertiary structure, leading us to propose that AaDUF507 represents a new protein fold.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Insight into molecular basis and dynamics of full-length CRaf kinase in cellular signaling mechanisms

Raf kinases play key roles in signal transduction in cells for regulating proliferation, differentiation, and survival. Despite decades of research into functions and dynamics of Raf kinases with respect to other cytosolic proteins, understanding Raf kinases is limited by the lack of their full-length structures at the atomic resolution. Here, we present the first model of the full-length CRaf kinase obtained from artificial intelligence/machine learning algorithms with a converging ensemble of structures simulated by large-scale temperature replica exchange simulations. Our model is validated by comparing simulated structures with the latest cryo-EM structure detailing close contacts among three key domains and regions of the CRaf. Our simulations identify potentially new epitopes of intramolecule interactions within the CRaf and reveal a dynamical nature of CRaf kinases, in which the three domains can move back and forth relative to each other for regulatory dynamics. The dynamic conformations are then used in a docking algorithm to shed insight into the paradoxical effect caused by vemurafenib in comparison with a paradox breaker PLX7904. In this study, we propose a model of Raf-heterodimer/KRas-dimer as a signalosome based on the dynamics of the full-length CRaf.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2001

Structural models of the MscL gating mechanism

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

NASA Discipline Cell Biology↗