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

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

Deciphering the Cofilin Oligomers via Intermolecular Disulfide Bond Formation: A Coarse-Grained Molecular Dynamics Approach to Understanding Cofilin’s Regulation on Actin Filaments

Cofilin, a key actin-binding protein, orchestrates the dynamics of the actomyosin network through its actin-severing activity and by promoting the recycling of actin monomers. Recent experimental work suggests that cofilin also forms functionally distinct oligomers through thiol post-translational modification (PTM) that encourages actin nucleation and assembly. Despite these advances, the structural conformations of cofilin oligomers that modulate actin activity remain elusive because there are combinatorial ways to oxidize thiols in cysteines to form disulfide bonds rapidly. This study employs molecular dynamics simulations to investigate human cofilin 1 as a case study for exploring cofilin dimers via disulfide bond formation. Using the free energy profiling, our simulations unveil a range of probable cofilin dimer structures not represented in current Protein Data Bank entries. These candidate dimers are characterized by their distinct population distributions and relative free energies. Of particular note is a dimer featuring an interface between cysteines 139 and 147 residues, which demonstrates stable free energy characteristics and intriguingly symmetrical geometry. In contrast, the experimentally proposed dimer structure exhibits a less stable free energy profile. Here, we also evaluate frustration quantification based on the energy landscape theory in the protein-protein interactions at the dimer interfaces. Notably, the 39-39 dimer configuration emerges as a promising candidate for forming cofilin tetramers, as substantiated by frustration analysis. Additionally, docking simulations with actin filaments further evaluate the stability of these cofilin dimer-actin complexes. Our findings thus offer a computational framework for understanding the role of thiol post-translational modification cofilin proteins in regulating oligomerization, and the subsequent cofilin-mediated actin dynamics in the actomyosin network.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Deep learning-driven insights into super protein complexes for outer membrane protein biogenesis in bacteria

To reach their final destinations, outer membrane proteins (OMPs) of gram-negative bacteria undertake an eventful journey beginning in the cytosol. Multiple molecular machines, chaperones, proteases, and other enzymes facilitate the translocation and assembly of OMPs. These helpers usually associate, often transiently, forming large protein assemblies. They are not well understood due to experimental challenges in capturing and characterizing protein-protein interactions (PPIs), especially transient ones. Using AF2Complex, we introduce a high-throughput, deep learning pipeline to identify PPIs within the Escherichia coli cell envelope and apply it to several proteins from an OMP biogenesis pathway. Among the top confident hits obtained from screening ~1500 envelope proteins, we find not only expected interactions but also unexpected ones with profound implications. Subsequently, we predict atomic structures for these protein complexes. These structures, typically of high confidence, explain experimental observations and lead to mechanistic hypotheses for how a chaperone assists a nascent, precursor OMP emerging from a translocon, how another chaperone prevents it from aggregating and docks to a β-barrel assembly port, and how a protease performs quality control. This work presents a general strategy for investigating biological pathways by using structural insights gained from deep learning-based predictions.

60 APPLIED LIFE SCIENCES↗

De novo design and Rosetta-based assessment of high-affinity antibody variable regions (Fv) against the SARS-CoV -2 spike receptor binding domain ( RBD )

The continued emergence of new SARS-CoV-2 variants has accentuated the growing need for fast and reliable methods for the design of potentially neutralizing antibodies (Abs) to counter immune evasion by the virus. Here, we report on the de novo computational design of high-affinity Ab variable regions (Fv) through the recombination of VDJ genes targeting the most solvent-exposed hACE2-binding residues of the SARS-CoV-2 spike receptor binding domain (RBD) protein using the software tool OptMAVEn-2.0. Subsequently, we carried out computational affinity maturation of the designed variable regions through amino acid substitutions for improved binding with the target epitope. Immunogenicity of designs was restricted by preferring designs that match sequences from a 9-mer library of “human Abs” based on a human string content score. We generated 106 different antibody designs and reported in detail on the top five that trade-off the greatest computational binding affinity for the RBD with human string content scores. We further describe computational evaluation of the top five designs produced by OptMAVEn-2.0 using a Rosetta-based approach. We used Rosetta SnugDock for local docking of the designs to evaluate their potential to bind the spike RBD and performed “forward folding” with DeepAb to assess their potential to fold into the designed structures. Ultimately, our results identified one designed Ab variable region, P1.D1, as a particularly promising candidate for experimental testing. This effort puts forth a computational workflow for the de novo design and evaluation of Abs that can quickly be adapted to target spike epitopes of emerging SARS-CoV-2 variants or other antigenic targets.

59 BASIC BIOLOGICAL SCIENCES↗

An evolutionarily conserved tryptophan cage promotes folding of the extended RNA recognition motif in the hnRNPR ‐like protein family

Abstract The heterogeneous nuclear ribonucleoprotein (hnRNP) R‐like family is a class of RNA binding proteins in the hnRNP superfamily with diverse functions in RNA processing. Here, we present the 1.90 Å X‐ray crystal structure and solution NMR studies of the first RNA recognition motif (RRM) of human hnRNPR. We find that this domain adopts an extended RRM (eRRM1) featuring a canonical RRM with a structured N‐terminal extension (N ext ) motif that docks against the RRM and extends the β‐sheet surface. The adjoining loop is structured and forms a tryptophan cage motif to position the N ext motif for docking to the RRM. Combining mutagenesis, solution NMR spectroscopy, and thermal denaturation studies, we evaluate the importance of residues in the N ext –RRM interface and adjoining loop on eRRM folding and conformational dynamics. We find that these sites are essential for protein solubility, conformational ordering, and thermal stability. Consistent with their importance, mutations in the N ext –RRM interface and loop are associated with several cancers in a survey of somatic mutations in cancer studies. Sequence and structure comparison of the human hnRNPR eRRM1 to experimentally verified and predicted hnRNPR‐like proteins reveals conserved features in the eRRM.

Biochemistry & Molecular Biology↗

Towards an Understanding of Oxidative Damage in an a-L-Arabinofuranosidase of Trichoderma reesei : a Molecular Dynamics Approach

Trichoderma reesei is a ‘workhorse’ fungus that produces glycosyl hydrolases (e.g., cellulases) at high titers for use in industrial bioprocessing. In this study, we focused on two side-chain hemicellulases, a-L-arabinofuranosidases and a-D-galactosidases, which are important for the treatment of lignocellulosic biomass, but susceptible to oxidative damage that can occur during industrial processing. The molecular details that render these enzymes inactive have not yet been identified. To approach this issue, we used proteomics to identify amino acid residues that were oxidized after a relevant oxidative treatment (Fenton reaction). These oxidative modifications were included in the 3D-protein structures, and using molecular dynamics simulations, we then studied the behaviors of non-modified and oxidized enzymes. These simulations showed significant alterations of the conformational stability of proteins when oxidized, as evidenced by changes in root mean square deviation (RMSD) and principal component analyses (PCA) trajectories. Likewise, enzyme-ligand interactions such as hydrogen bonds were greatly reduced in quantity and quality in the oxidized proteins. Finally, free energy landscape plots showed that there was a more rugged energy surface in the oxidized proteins, implying a less favorable reaction pathway. These results reveal the basis for loss of function in these carbohydrate active enzymes (CAZYs) in T. reesei.

60 APPLIED LIFE SCIENCES↗

The structural basis for the broad aldehyde specificity of the aminoaldehyde dehydrogenase PauC from the human pathogen Pseudomonas aeruginosa

Abstract Despite significant differences in size and formal charge, the aldehyde dehydrogenasePaPauC (PA5312) fromPseudomonas aeruginosaPAO1 efficiently catalyzes the NAD + ‐dependent oxidation of the aminoaldehydes formed in polyamines degradation. We report here thatPaPauC also oxidizes 4‐guanidinebutyraldehyde, formed in one arginine degradation pathway, trimethylaminobutyraldehyde, of unknown metabolic origin, and indole‐3‐acetaldehyde, a precursor of the plant growth‐promoting hormone indoleacetic acid.PaPauC has been proposed as a potential target for combatingP. aeruginosa. However, understanding its structure–function relationships, crucial for developing specific inhibitors, is lacking. Using X‐ray crystallography, we identified the structural characteristics that determinePaPauC broad aldehyde specificity: a spacious aldehyde‐entrance tunnel and six active‐site residues. Docking simulations, site‐directed mutagenesis, and kinetic analyses support the interactions of Lys479 with glutamylated aminoaldehydes; Phe169, Trp176, and Phe467 with amino and guanidinium groups through cation–π interactions and with the indole group via NH–π and CH–π interactions; Asp459 with amino and indole groups; and Thr303 with amide and guanidinium groups. Exploiting the distinctive structural features of thePaPauC active site could aid in developing specific inhibitors to combatP. aeruginosainfections in humans and animals, as well as in preventing its colonization of plants, which are abundantP. aeruginosareservoirs and, therefore, a significant source of human infections.

Biochemistry & Molecular Biology↗

An effective human uracil-DNA glycosylase inhibitor targets the open pre-catalytic active site conformation

Human uracil DNA-glycosylase (UDG) is the prototypic and first identified DNA glycosylase with a vital role in removing deaminated cytosine and incorporated uracil and 5-fluorouracil (5-FU) from DNA. UDG depletion sensitizes cells to high APOBEC3B deaminase and to pemetrexed (PEM) and floxuridine (5-FdU), which are toxic to tumor cells through incorporation of uracil and 5-FU into DNA. To identify small-molecule UDG inhibitors for pre-clinical evaluation, we optimized biochemical screening of a selected diversity collection of >3,000 small-molecules. We found aurintricarboxylic acid (ATA) as an inhibitor of purified UDG at an initial calculated IC 50 < 100 nM. Subsequent enzymatic assays confirmed effective ATA inhibition but with an IC50 of 700 nM and showed direct binding to the human UDG with a K D of <700 nM. ATA displays preferential, dose-dependent binding to purified human UDG compared to human 8-oxoguanine DNA glycosylase. ATA did not bind uracil-containing DNA at these concentrations. Yet, combined crystal structure and in silico docking results unveil ATA interactions with the DNA binding channel and uracil-binding pocket in an open, destabilized UDG conformation. Biologically relevant ATA inhibition of UDG was measured in cell lysates from human DLD1 colon cancer cells and in MCF-7 breast cancer cells using a host cell reactivation assay. Collective findings provide proof-of-principle for development of an ATA-based chemotype and “door stopper” strategy targeting inhibitor binding to a destabilized, open pre-catalytic glycosylase conformation that prevents active site closing for functional DNA binding and nucleotide flipping needed to excise altered bases in DNA.

59 BASIC BIOLOGICAL SCIENCES↗

Host receptor-targeted therapeutic approach to counter pathogenic New World mammarenavirus infections

Five New World mammarenaviruses (NWMs) cause life-threatening hemorrhagic fever (HF). Cellular entry by these viruses is mediated by human transferrin receptor 1 (hTfR1). Here, we demonstrate that an antibody (ch128.1/IgG1) which binds the apical domain of hTfR1, potently inhibits infection of attenuated and pathogenic NWMs in vitro. Computational docking of the antibody Fab crystal structure onto the known structure of hTfR1 shows an overlapping receptor-binding region shared by the Fab and the viral envelope glycoprotein GP1 subunit that binds hTfR1, and we demonstrate competitive inhibition of NWM GP1 binding by ch128.1/IgG1 as the principal mechanism of action. Importantly, ch128.1/IgG1 protects hTfR1-expressing transgenic mice against lethal NWM challenge. Additionally, the antibody is well-tolerated and only partially reduces ferritin uptake. Our findings provide the basis for the development of a novel, host receptor-targeted antibody therapeutic broadly applicable to the treatment of HF of NWM etiology.

59 BASIC BIOLOGICAL SCIENCES↗