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 55 records · Page 3

Identification of glycosyltransferases mediating 2‐ O ‐arabinopyranosyl and 2‐ O ‐galactosyl substitutions of glucuronosyl side chains of xylan

SUMMARY Xylan is one of the major hemicelluloses in plant cell walls and its xylosyl backbone is often decorated at O ‐2 with glucuronic acid (GlcA) and/or methylglucuronic acid (MeGlcA) residues. The GlcA/MeGlcA side chains may be further substituted with 2‐ O ‐arabinopyranose (Ara p ) or 2‐ O ‐galactopyranose (Gal) residues in some plant species, but the enzymes responsible for these substitutions remain unknown. During our endeavor to investigate the enzymatic activities of Arabidopsis MUR3‐clade members of the GT47 glycosyltransferase family, we found that one of them was able to transfer Ara p from UDP‐Ara p onto O ‐2 of GlcA side chains of xylan, and thus it was named xylan 2‐ O ‐arabinopyranosyltransferase 1 (AtXAPT1). The function of AtXAPT1 was verified in planta by its T‐DNA knockout mutation showing a loss of the Ara p substitution on xylan GlcA side chains. Further biochemical characterization of XAPT close homologs from other plant species demonstrated that while the poplar ones had the same catalytic activity as AtXAPT1, those from Eucalyptus , lemon‐scented gum, sea apple, 'Ohi'a lehua, duckweed and purple yam were capable of catalyzing both 2‐ O ‐Ara p and 2‐ O ‐Gal substitutions of xylan GlcA side chains albeit with differential activities. Sequential reactions with XAPTs and glucuronoxylan methyltransferase 3 (GXM3) showed that XAPTs acted poorly on MeGlcA side chains, whereas GXM3 could efficiently methylate arabinosylated or galactosylated GlcA side chains of xylan. Furthermore, molecular docking and site‐directed mutagenesis analyses of Eucalyptus XAPT1 revealed critical roles of several amino acid residues at the putative active site in its activity. Together, these findings establish that XAPTs residing in the MUR3 clade of family GT47 are responsible for 2‐ O ‐arabinopyranosylation and 2‐ O ‐galactosylation of GlcA side chains of xylan.

Zhong, Ruiqin↗

High-throughput virtual laboratory for drug discovery using massive datasets

Time-to-solution for structure-based screening of massive chemical databases for COVID-19 drug discovery has been decreased by an order of magnitude, and a virtual laboratory has been deployed at scale on up to 27,612 GPUs on the Summit supercomputer, allowing an average molecular docking of 19,028 compounds per second. Over one billion compounds were docked to two SARS-CoV-2 protein structures with full optimization of ligand position and 20 poses per docking, each in under 24 hours. GPU acceleration and high-throughput optimizations of the docking program produced 350× mean speedup over the CPU version (50× speedup per node). GPU acceleration of both feature calculation for machine-learning based scoring and distributed database queries reduced processing of the 2.4 TB output by orders of magnitude. The resulting 50× speedup for the full pipeline reduces an initial 43 day runtime to 21 hours per protein for providing high-scoring compounds to experimental collaborators for validation assays.

97 MATHEMATICS AND COMPUTING↗

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↗

Characterization of a Haematococcus pluvialis Diacylglycerol Acyltransferase 1 and Its Potential in Unsaturated Fatty Acid-Rich Triacylglycerol Production

The unicellular green alga Haematococcus pluvialis has been recognized as an industry strain to produce simultaneously esterified astaxanthin (EAST) and triacylglycerol (TAG) under stress induction. It is necessary to identify the key enzymes involving in synergistic accumulation of EAST and TAG in H. pluvialis . In this study, a novel diacylglycerol acyltransferase 1 was systematically characterized by in vivo and in silico assays. The upregulated expression of HpDGAT1 gene was positively associated with the significant increase of TAG and EAST contents under stress conditions. Functional complementation by overexpressing HpDGAT1 in a TAG-deficient yeast strain H1246 revealed that HpDGAT1 could restore TAG biosynthesis and exhibited a high substrate preference for monounsaturated fatty acyl-CoAs (MUFAs) and polyunsaturated fatty acyl-CoAs (PUFAs). Notably, heterogeneous expression of HpDGAT1 in Chlamydomonas reinhardtii and Arabidopsis thaliana resulted in a significant enhancement of total oils and concurrently a high accumulation of MUFAs- and PUFAs-rich TAGs. Furthermore, molecular docking analysis indicated that HpDGAT1 contained AST-binding sites. These findings evidence a possible dual-function role for HpDGAT1 involving in TAG and EAST synthesis, demonstrating that it is a potential target gene to enrich AST accumulation in this alga and to design oil production in both commercial algae and oil crops.

Cui, Hongli↗

Indoline-6-Sulfonamide Inhibitors of the Bacterial Enzyme DapE

Inhibitors of the bacterial enzyme dapE-encoded N-succinyl-l,l-diaminopimelic acid desuccinylase (DapE; EC 3.5.1.18) hold promise as antibiotics with a new mechanism of action. Herein we describe the discovery of a new series of indoline sulfonamide DapE inhibitors from a high-throughput screen and the synthesis of a series of analogs. Inhibitory potency was measured by a ninhydrin-based DapE assay recently developed by our group. Molecular docking experiments suggest active site binding with the sulfonamide acting as a zinc-binding group (ZBG).

59 BASIC BIOLOGICAL SCIENCES↗

Improving Protein–Ligand Interaction Modeling with cryo-EM Data, Templates, and Deep Learning in 2021 Ligand Model Challenge

Elucidating protein–ligand interaction is crucial for studying the function of proteins and compounds in an organism and critical for drug discovery and design. The problem of protein–ligand interaction is traditionally tackled by molecular docking and simulation, which is based on physical forces and statistical potentials and cannot effectively leverage cryo-EM data and existing protein structural information in the protein–ligand modeling process. In this work, we developed a deep learning bioinformatics pipeline (DeepProLigand) to predict protein–ligand interactions from cryo-EM density maps of proteins and ligands. DeepProLigand first uses a deep learning method to predict the structure of proteins from cryo-EM maps, which is averaged with a reference (template) structure of the proteins to produce a combined structure to add ligands. The ligands are then identified and added into the structure to generate a protein–ligand complex structure, which is further refined. The method based on the deep learning prediction and template-based modeling was blindly tested in the 2021 EMDataResource Ligand Challenge and was ranked first in fitting ligands to cryo-EM density maps. These results demonstrate that the deep learning bioinformatics approach is a promising direction for modeling protein–ligand interactions on cryo-EM data using prior structural information.

59 BASIC BIOLOGICAL SCIENCES↗

Characterization of Interactions between CTX-M-15 and Clavulanic Acid, Desfuroylceftiofur, Ceftiofur, Ampicillin, and Nitrocefin

Cefotaximase-Munich (CTX-M) extended-spectrum beta-lactamases (ESBLs) are commonly associated with Gram-negative, hospital-acquired infections worldwide. Several beta-lactamase inhibitors, such as clavulanate, are used to inhibit the activity of these enzymes. To understand the mechanism of CTX-M-15 activity, we have determined the crystal structures of CTX-M-15 in complex with two specific classes of beta-lactam compounds, desfuroylceftiofur (DFC) and ampicillin, and an inhibitor, clavulanic acid. The crystal structures revealed that Ser70 and five other residues (Lys73, Tyr105, Glu166, Ser130, and Ser237) participate in catalysis and binding of those compounds. Based on analysis of steady-state kinetics, thermodynamic data, and molecular docking to both wild-type and S70A mutant structures, we determined that CTX-M-15 has a similar affinity for all beta-lactam compounds (ceftiofur, nitrocefin, DFC, and ampicillin), but with lower affinity for clavulanic acid. A catalytic mechanism for tested β-lactams and two-step inhibition mechanism of clavulanic acid were proposed. CTX-M-15 showed a higher activity toward DFC and nitrocefin, but significantly lower activity toward ampicillin and ceftiofur. The interaction between CTX-M-15 and both ampicillin and ceftiofur displayed a higher entropic but lower enthalpic effect, compared with DFC and nitrocefin. DFC, a metabolite of ceftiofur, displayed lower entropy and higher enthalpy than ceftiofur. This finding suggests that compounds containing amine moiety (e.g., ampicillin) and the furfural moiety (e.g., ceftiofur) could hinder the hydrolytic activity of CTX-M-15.

59 BASIC BIOLOGICAL SCIENCES↗

Structural Similarities and Overlapping Activities among Dihydroflavonol 4-Reductase, Flavanone 4-Reductase, and Anthocyanidin Reductase Offer Metabolic Flexibility in the Flavonoid Pathway

Flavonoids are potent antioxidants that play a role in defense against pathogens, UV-radiation, and the detoxification of reactive oxygen species. Dihydroflavonol 4-reductase (DFR) and flavanone 4-reductase (FNR) reduce dihydroflavonols and flavanones, respectively, using NAD(P)H to produce flavan-(3)-4-(di)ols in flavonoid biosynthesis. Anthocyanidin reductase (ANR) reduces anthocyanidins to flavan-3-ols. In addition to their sequences, the 3D structures of recombinant DFR, FNR and ANR from sorghum and switchgrass showed a high level of similarity. The catalytic mechanism, substrate-specificity and key residues of three reductases were deduced from crystal structures, site-directed mutagenesis, molecular docking, kinetics, and thermodynamic ana-lyses. Although DFR displayed its highest activity against dihydroflavonols, it also showed activity against flavanones and anthocyanidins. It was inhibited by the flavonol quercetin and high concentrations of dihydroflavonols/flavonones. SbFNR1 and SbFNR2 did not show any activity against dihydroflavonols. However, SbFNR1 displayed activity against flavanones and ANR activity against two anthocyanidins, cyanidin and pelargonidin. Therefore, SbFNR1 and SbFNR2 could be specific ANR isozymes without delphinidin activity. Sorghum has high concentrations of 3-deoxyanthocyanidins in vivo, supporting the observed high activity of SbDFR against flavonols. Mining of expression data indicated substantial induction of these three reductase genes in both switchgrass and sorghum in response to biotic stress. Key signature sequences for proper DFR/ANR classification are proposed and could form the basis for future metabolic engineering of flavonoid metabolism.

3-deoxyanthocyanidin↗

Structural and Interactional Analysis of the Flavonoid Pathway Proteins: Chalcone Synthase, Chalcone Isomerase and Chalcone Isomerase-like Protein

Chalcone synthase (CHS) and chalcone isomerase (CHI) catalyze the first two committed steps of the flavonoid pathway that plays a pivotal role in the growth and reproduction of land plants, including UV protection, pigmentation, symbiotic nitrogen fixation, and pathogen resistance. Based on the obtained X-ray crystal structures of CHS, CHI, and chalcone isomerase-like protein (CHIL) from the same monocotyledon, Panicum virgatum, along with the results of the steady-state kinetics, spectroscopic/thermodynamic analyses, intermolecular interactions, and their effect on each catalytic step are proposed. In addition, PvCHI’s unique activity for both naringenin chalcone and isoliquiritigenin was analyzed, and the observed hierarchical activity for those type-I and -II substrates was explained with the intrinsic characteristics of the enzyme and two substrates. The structure of PvCHS complexed with naringenin supports uncompetitive inhibition. PvCHS displays intrinsic catalytic promiscuity, evident from the formation of p-coumaroyltriacetic acid lactone (CTAL) in addition to naringenin chalcone. In the presence of PvCHIL, conversion of p-coumaroyl-CoA to naringenin through PvCHS and PvCHI displayed ~400-fold increased Vmax with reduced formation of CTAL by 70%. Supporting this model, molecular docking, ITC (Isothermal Titration Calorimetry), and FRET (Fluorescence Resonance Energy Transfer) indicated that both PvCHI and PvCHIL interact with PvCHS in a non-competitive manner, indicating the plausible allosteric effect of naringenin on CHS. Significantly, the presence of naringenin increased the affinity between PvCHS and PvCHIL, whereas naringenin chalcone decreased the affinity, indicating a plausible feedback mechanism to minimize spontaneous incorrect stereoisomers. These are the first findings from a three-body system from the same species, indicating the importance of the macromolecular assembly of CHS-CHI-CHIL in determining the amount and type of flavonoids produced in plant cells.

59 BASIC BIOLOGICAL SCIENCES↗

A Novel Glutathione S-Transferase Gtt2 Class (VpGSTT2) Is Found in the Genome of the AHPND/EMS Vibrio parahaemolyticus Shrimp Pathogen

Glutathione S-transferases are a family of detoxifying enzymes that catalyze the conjugation of reduced glutathione (GSH) with different xenobiotic compounds using either Ser, Tyr, or Cys as a primary catalytic residue. We identified a novel GST in the genome of the shrimp pathogen V. parahaemolyticus FIM- S1708+, a bacterial strain associated with Acute Hepatopancreatic Necrosis Disease (AHPND)/Early Mortality Syndrome (EMS) in cultured shrimp. This new GST class was named Gtt2. It has an atypical catalytic mechanism in which a water molecule instead of Ser, Tyr, or Cys activates the sulfhydryl group of GSH. The biochemical properties of Gtt2 from Vibrio parahaemolyticus (VpGSTT2) were characterized using kinetic and crystallographic methods. Recombinant VpGSTT2 was enzymatically active using GSH and CDNB as substrates, with a specific activity of 5.7 units/mg. Low affinity for substrates was demonstrated using both Michaelis–Menten kinetics and isothermal titration calorimetry. The crystal structure showed a canonical two-domain structure comprising a glutathione binding G-domain and a hydrophobic ligand H domain. A water molecule was hydrogen-bonded to residues Thr9 and Ser 11, as reported for the yeast Gtt2, suggesting a primary role in the reaction. Molecular docking showed that GSH could bind at the G-site in the vicinity of Ser11. G-site mutationsT9A and S11A were analyzed. S11A retained 30% activity, while T9A/S11A showed no detectable activity. VpGSTT2 was the first bacterial Gtt2 characterized, in which residues Ser11 and Thr9 coordinated a water molecule as part of a catalytic mechanism that was characteristic of yeast GTT2. The GTT2 family has been shown to provide protection against metal toxicity; in some cases, excess heavy metals appear in shrimp ponds presenting AHPND/EMS. Further studies may address whether GTT2 in V. parahaemolyticus pathogenic strains may provide a competitive advantage as a novel detoxification mechanism.

59 BASIC BIOLOGICAL SCIENCES↗

Virtual reality at work

The utility of virtual reality computer graphics in telepresence applications is not hard to grasp and promises to be great. When the virtual world is entirely synthetic, as opposed to real but remote, the utility is harder to establish. Vehicle simulators for aircraft, vessels, and motor vehicles are proving their worth every day. Entertainment applications such as Disney World's StarTours are technologically elegant, good fun, and economically viable. Nevertheless, some of us have no real desire to spend our lifework serving the entertainment craze of our sick culture; we want to see this exciting technology put to work in medicine and science. The topics covered include the following: testing a force display for scientific visualization -- molecular docking; and testing a head-mounted display for scientific and medical visualization.

Brooks, Frederick P., Jr.↗

Self-consistent mean-field quantum approximate optimization

We introduce a self-consistent mean-field quantum optimization algorithm that approximates the ground state of classical Ising Hamiltonians. The algorithm decomposes the problem into independent subproblems and treats the interactions between them in a mean-field manner. These interactions are captured by a common environment, constructed self-consistently through a variational quantum circuit, and which modifies the subproblems to account for mutual influence while maintaining computational independence. Consequently, subproblems can be solved individually, avoiding the computational cost of the full problem. We explore the properties of the generated environment and assess the algorithm's performance through extensive numerical simulations on Sherrington-Kirkpatrick spin glasses. Furthermore, we apply it experimentally to a weighted maximum clique problem applied to molecular docking. This framework enables the solution of problems that would otherwise exceed the qubit and gate counts of current quantum hardware.

Dupont, Maxime [Rigetti Computing] (ORCID:00000001↗

Self-consistent mean-field quantum approximate optimization

We introduce a self-consistent mean-field quantum optimization algorithm that approximates the ground state of classical Ising Hamiltonians. The algorithm decomposes the problem into independent subproblems and treats the interactions between them in a mean-field manner. These interactions are captured by a common environment, constructed self-consistently through a variational quantum circuit, and which modifies the subproblems to account for mutual influence while maintaining computational independence. Consequently, subproblems can be solved individually, avoiding the computational cost of the full problem. We explore the properties of the generated environment and assess the algorithm's performance through extensive numerical simulations on Sherrington-Kirkpatrick spin glasses. Furthermore, we apply it experimentally to a weighted maximum clique problem applied to molecular docking. This framework enables the solution of problems that would otherwise exceed the qubit and gate counts of current quantum hardware.

Dupont, Maxime [Rigetti Computing] (ORCID:00000001↗

Self-consistent mean-field quantum approximate optimization

We introduce a self-consistent mean-field quantum optimization algorithm that approximates the ground state of classical Ising Hamiltonians. The algorithm decomposes the problem into independent subproblems and treats the interactions between them in a mean-field manner. These interactions are captured by a common environment, constructed self-consistently through a variational quantum circuit, and which modifies the subproblems to account for mutual influence while maintaining computational independence. Consequently, subproblems can be solved individually, avoiding the computational cost of the full problem. We explore the properties of the generated environment and assess the algorithm's performance through extensive numerical simulations on Sherrington-Kirkpatrick spin glasses. Furthermore, we apply it experimentally to a weighted maximum clique problem applied to molecular docking. This framework enables the solution of problems that would otherwise exceed the qubit and gate counts of current quantum hardware.

Dupont, Maxime [Rigetti Computing] (ORCID:00000001↗

Auto_PDI

Protein-DNA Interaction Workflow (PDI Workflow), a pipeline that focuses on generating high-quality docking and molecular dynamics simulations for Protein-DNA complexes. This allows us to take DNA sequences with unknown tertiary structures, accurately predict their structure, dock them with the desired target protein, and then simulate their interactions using molecular dynamics simulations

Kumar, Neeraj↗

Functional Characterization and Structural Insights Into Stereoselectivity of Pulegone Reductase in Menthol Biosynthesis

Monoterpenoids are the main components of plant essential oils and the active components of some traditional Chinese medicinal herbs like Mentha haplocalyx Briq ., Nepeta tenuifolia Briq ., Perilla frutescens (L.) Britt and Pogostemin cablin (Blanco) Benth. Pulegone reductase is the key enzyme in the biosynthesis of menthol and is required for the stereoselective reduction of the Δ 2,8 double bond of pulegone to produce the major intermediate menthone, thus determining the stereochemistry of menthol. However, the structural basis and mechanism underlying the stereoselectivity of pulegone reductase remain poorly understood. In this study, we characterized a novel (-)-pulegone reductase from Nepeta tenuifolia ( Nt PR), which can catalyze (-)-pulegone to (+)-menthone and (-)-isomenthone through our RNA-seq, bioinformatic analysis in combination with in vitro enzyme activity assay, and determined the structure of (+)-pulegone reductase from ( Mp PR) by using X-ray crystallography, molecular modeling and docking, site-directed mutagenesis, molecular dynamics simulations, and biochemical analysis. We identified and validated the critical residues in the crystal structure of Mp PR involved in the binding of the substrate pulegone. We also further identified that residues Leu56, Val282, and Val284 determine the stereoselectivity of the substrate pulegone, and mainly contributes to the product stereoselectivity. This work not only provides a starting point for the understanding of stereoselectivity of pulegone reductases, but also offers a basis for the engineering of menthone/menthol biosynthetic enzymes to achieve high-titer, industrial-scale production of enantiomerically pure products.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Modeling Protein–Protein and Protein–Ligand Interactions by the ClusPro Team in CASP16

ABSTRACT In the CASP16 experiment, our team employed hybrid computational strategies to predict both protein–protein and protein–ligand complex structures. For protein–protein docking, we combined physics‐based sampling—using ClusPro FFT docking and molecular dynamics—with AlphaFold (AF)‐based sampling, followed by AF‐based refinement. Our method produced numerous high‐accuracy complex models, including cases where AF alone failed, underscoring the critical role of physics‐based sampling alongside deep learning‐based refinement. For protein–ligand docking, we integrated the ClusPro LigTBM template‐based approach with a machine learning‐based confidence model for rescoring. The method preserves conserved interaction fragments derived from homologous complexes, followed by local resampling using physics‐based sampling and a diffusion model. Our template‐based strategy achieved a mean lDDT‐PLI of 0.69 across 233 targets, which was highly competitive. These results demonstrate that combining physics‐based modeling with AI‐driven refinement can significantly enhance the accuracy of both protein–protein and protein–ligand structure predictions.

Ashizawa, Ryota [Department of Applied Mathematics↗

Novel sterol binding domains in bacteria

Sterol lipids are widely present in eukaryotes and play essential roles in signaling and modulating membrane fluidity. Although rare, some bacteria also produce sterols, but their function in bacteria is not known. Moreover, many more species, including pathogens and commensal microbes, acquire or modify sterols from eukaryotic hosts through poorly understood molecular mechanisms. The aerobic methanotroph Methylococcus capsulatus was the first bacterium shown to synthesize sterols, producing a mixture of C-4 methylated sterols that are distinct from those observed in eukaryotes. C-4 methylated sterols are synthesized in the cytosol and localized to the outer membrane, suggesting that a bacterial sterol transport machinery exists. Until now, the identity of such machinery remained a mystery. In this study, we identified three novel proteins that may be the first examples of transporters for bacterial sterol lipids. The proteins, which all belong to well-studied families of bacterial metabolite transporters, are predicted to reside in the inner membrane, periplasm, and outer membrane of M. capsulatus, and may work as a conduit to move modified sterols to the outer membrane. Quantitative analysis of ligand binding revealed their remarkable specificity for 4-methylsterols, and crystallographic structures coupled with docking and molecular dynamics simulations revealed the structural bases for substrate binding by two of the putative transporters. Their striking structural divergence from eukaryotic sterol transporters signals that they form a distinct sterol transport system within the bacterial domain. Finally, bioinformatics revealed the widespread presence of similar transporters in bacterial genomes, including in some pathogens that use host sterol lipids to construct their cell envelopes. The unique folds of these bacterial sterol binding proteins should now guide the discovery of other proteins that handle this essential metabolite.

59 BASIC BIOLOGICAL SCIENCES↗