SEARCH · Search NASA
Results for “docking simulations”
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.
Dirigent isoflavene-forming PsPTS2: 3D structure, stereochemical, and kinetic characterization comparison with pterocarpan-forming PsPTS1 homolog in pea
Pea (legume plant defense) phytoalexin pterocarpans (–)-maackiain and (+)-pisatin have opposite C6a and C11a configurations, but how this occurs biosynthetically is unknown. Herein, pea dirigent-protein (DP) PsPTS2, which generates 7,2'-dihydroxy-4',5'-methylenedioxyisoflav-3-ene (DMDIF), was investigated as regards stereoselectivity towards four possible 7,2'-dihydroxy-4',5'-methylenedioxyisoflavan-4-ol (DMDI) stereoisomers. Each DMDI stereoisomer configuration was determined using NMR spectroscopy, electronic circular dichroism, and molecular orbital analyses. PsPTS2 efficiently converted cis-(3R,4R)-DMDI into DMDIF, this being 20-fold faster than with the trans-(3R,4S)-isomer. The 4R-configured substrate’s near ß-axial OH orientation significantly enhanced its leaving group abilities in creating intermediate A-ring mono-quinone methide (QM), whereas the 4S-isomer’s 4-OH group α-equatorial orientation was a poorer leaving group. Substrate docking simulations also indicated that the 4R-configured ß-axial OH was closest to the Asp catalytic center, whereas the 4S-isomer’s α-equatorial OH was further away. Neither cis-(3S,4S)- nor trans-(3S, 4R)-DMDIs were substrates, even with the former having the same C3/C4 stereochemistry as (+)-pisatin. PsPTS2, assayed with cis-(3R,4R)-7,2'-dihydroxy-4'-methoxyisoflavan-4-ol [cis-(3R,4R)-DMI] and its C3/C4 stereoisomers, gave 2',7-dihydroxy-4'-methoxyisoflav-3-ene (DMIF), with the same substrate stereoselectivity. DP homologs may exist in licorice (Glycyrrhiza pallidiflora) and the tree legume Bolusanthus speciosus, as DMIF is found in both. PsPTS1 preferentially utilized cis-(3R,4R)-DMDI to give (–)-maackiain, this being 2200-fold more efficient than with cis-(3R,4R)-DMI. PsPTS1 also slowly converted trans-(3S,4R)-DMDI into (+)-maackiain, again reflecting the better 4R configured OH leaving group. Pea PsPTS2 and PsPTS1 provisionally provide a facile DP-engendered means to enable differing C6a and C11a configurations in (+)-pisatin and (–)-maackiain. Both DP transformations are considered to occur via DP-engendered generation of an identical mono-QM bound intermediate, which with PsPTS2 either re-aromatizes to give DMDIF or with PsPTS1 undergoes intramolecular cyclization to afford (–)-maackiain. Substrate docking simulations using PsPTS2 (2.5 Å resolution) and PsPTS1 (1.5 Å resolution) plausibly indicate that cis-(3R,4R)-DMDI binds in the anti-configuration in the PsPTS2 active site to afford DMDIF, and the syn-configuration in PsPTS1 to give maackiain.
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.
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.
Predicting compatibility between ferredoxins and the Fe protein of nitrogenase using in silico protein modeling
Biological nitrogen fixation is the process by which certain bacteria and archaea use the enzyme nitrogenase to reduce atmospheric nitrogen into bioavailable ammonium. Engineering non‐nitrogen‐fixing organisms, like plants, to use nitrogenase could reduce dependency on synthetic fertilizer and mitigate the environmental impacts of industrial fertilizer production. However, nitrogenase activity requires delivery of reducing power by small electron carrying proteins known as ferredoxins and flavodoxins, and successfully engineering nitrogenase into new systems will require a mechanistic understanding of electron delivery by these proteins. Most organisms often have multiple ferredoxins, raising the question of which ferredoxin can support nitrogenase activity. The purpose of this study is to gain insight into how we can predict which ferredoxin is compatible with the Fe protein, the component of nitrogenase that interacts with ferredoxin or flavodoxin. Our in silico protein–protein docking simulations reveal that most ferredoxins and flavodoxins involved in nitrogen fixation have the shortest distance (≤10 Å) between their redox cofactor and the [4Fe‐4S] cluster of the Fe protein. We found shorter cofactor distance contributes to faster intermolecular electron tunneling rates. Bacterial ferredoxins that play a role in nitrogen fixation also exhibit more complementary interactions with the Fe protein than bacterial and plant ferredoxins not involved in this process. Heterologous expression of a set of ferredoxins from both nitrogen‐fixing and non‐nitrogen‐fixing bacteria in the diazotroph Rhodopseudomonas palustris supports our model‐derived prediction that shorter distances between the electron‐carrying cofactors favor nitrogenase compatibility. These findings offer a framework to predict and potentially enhance ferredoxin–nitrogenase compatibility, which will help to improve our ability to engineer nitrogen fixation into non‐nitrogen‐fixing organisms like plants.
Engineering a photoenzyme to use red light
Photoenzymatic reactions involving flavin-dependent “ene”-reductases (EREDs) rely on protein-templated charge transfer (CT) complexes between the cofactor and substrate for radical initiation. These complexes typically absorb in the blue region of the electromagnetic spectrum. Here, we engineered an ERED to form CT complexes that absorb red light. Mechanistic studies indicate that red-light activity is due to the growth of a red-absorbing shoulder off the previously identified cyan absorption feature. Molecular dynamics simulations, docking, and excited-state calculations suggest that the cyan feature involves a π→π* transition on flavin, whereas the red-light absorption is a π→π* transition between flavin and the substrate. Furthermore, differences in the electronic transition are due to changes in the substrate-binding conformation and allosteric tuning of the electronic structure of the cofactor-substrate complex. Microenvironment tuning of the CT complex for red-light activity is observed with other engineered photoenzymatic reactions, highlighting this effect’s generality.
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.
Structural insights into hydrolytic defluorination of difluoroacetate by microbial fluoroacetate dehalogenases
Fluorine forms the strongest single bond to carbon with the highest bond dissociation energy among natural products. However, fluoroacetate dehalogenases (FADs) have been shown to hydrolyze this bond in fluoroacetate under mild reaction conditions. Furthermore, two recent studies demonstrated that the FAD RPA1163 fromRhodopseudomonas palustriscan also accept bulkier substrates. In this study, we explored the substrate promiscuity of microbial FADs and their ability to defluorinate polyfluorinated organic acids. Enzymatic screening of eight purified dehalogenases with reported fluoroacetate defluorination activity revealed significant hydrolytic activity against difluoroacetate in three proteins. Product analysis using liquid chromatography‐mass spectrometry identified glyoxylic acid as the final product of enzymatic DFA defluorination. The crystal structures of DAR3835 fromDechloromonas aromaticaand NOS0089 fromNostocsp. were determined in the apo‐state along with the DAR3835 H274N glycolyl intermediate. Structure‐based site‐directed mutagenesis of DAR3835 demonstrated a key role for the catalytic triad and other active site residues in the defluorination of both fluoroacetate and difluoroacetate. Computational analysis of the dimer structures of DAR3835, NOS0089, and RPA1163 indicated the presence of one substrate access tunnel in each protomer. Moreover, protein‐ligand docking simulations suggested similar catalytic mechanisms for the defluorination of both fluoroacetate and difluoroacetate, with difluoroacetate being defluorinated via two consecutive defluorination reactions producing glyoxylate as the final product. Thus, our findings provide molecular insights into substrate promiscuity and catalytic mechanism of FADs, which are promising biocatalysts for applications in synthetic chemistry and bioremediation of fluorochemicals.
Elucidation of Agonist and Antagonist Dynamic Binding Patterns in ER-α by Integration of Molecular Docking, Molecular Dynamics Simulations and Quantum Mechanical Calculations
Estrogen receptor alpha (ERα) is a ligand-dependent transcriptional factor in the nuclear receptor superfamily. Many structures of ERα bound with agonists and antagonists have been determined. However, the dynamic binding patterns of agonists and antagonists in the binding site of ERα remains unclear. Therefore, we performed molecular docking, molecular dynamics (MD) simulations, and quantum mechanical calculations to elucidate agonist and antagonist dynamic binding patterns in ERα. 17β-estradiol (E2) and 4-hydroxytamoxifen (OHT) were docked in the ligand binding pockets of the agonist and antagonist bound ERα. The best complex conformations from molecular docking were subjected to 100 nanosecond MD simulations. Hierarchical clustering was conducted to group the structures in the trajectory from MD simulations. The representative structure from each cluster was selected to calculate the binding interaction energy value for elucidation of the dynamic binding patterns of agonists and antagonists in the binding site of ERα. The binding interaction energy analysis revealed that OHT binds ERα more tightly in the antagonist conformer, while E2 prefers the agonist conformer. The results may help identify ERα antagonists as drug candidates and facilitate risk assessment of chemicals through ER-mediated responses.
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.
Discovery of Small-Molecule Inhibitors of SARS-CoV-2 Proteins Using a Computational and Experimental Pipeline
A rapid response is necessary to contain emergent biological outbreaks before they can become pandemics. The novel coronavirus (SARS-CoV-2) that causes COVID-19 was first reported in December of 2019 in Wuhan, China and reached most corners of the globe in less than two months. In just over a year since the initial infections, COVID-19 infected almost 100 million people worldwide. Although similar to SARS-CoV and MERS-CoV, SARS-CoV-2 has resisted treatments that are effective against other coronaviruses. Crystal structures of two SARS-CoV-2 proteins, spike protein and main protease, have been reported and can serve as targets for studies in neutralizing this threat. We have employed molecular docking, molecular dynamics simulations, and machine learning to identify from a library of 26 million molecules possible candidate compounds that may attenuate or neutralize the effects of this virus. The viability of selected candidate compounds against SARS-CoV-2 was determined experimentally by biolayer interferometry and FRET-based activity protein assays along with virus-based assays. In the pseudovirus assay, imatinib and lapatinib had IC 50 values below 10 μM, while candesartan cilexetil had an IC 50 value of approximately 67 µM against M pro in a FRET-based activity assay. Comparatively, candesartan cilexetil had the highest selectivity index of all compounds tested as its half-maximal cytotoxicity concentration 50 (CC 50 ) value was the only one greater than the limit of the assay (>100 μM).
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
A unified simulation framework for wave energy powered underwater vehicle docking and charging
Not provided.
Multiscale Models for Fibril Formation: Rare Events Methods, Microkinetic Models, and Population Balances
Amyloid fibrils are thought to grow by a two-step dock-lock mechanism. However, previous simulations of fibril formation (i) overlook the bi-molecular nature of the docking step and obtain rates with first-order units, or (ii) superimpose the docked and locked states when computing the potential of mean force for association and thereby muddle the docking and locking steps. Here, we developed a simple microkinetic model with separate locking and docking steps and with the appropriate concentration dependences for each step. We constructed a simple model comprised of chiral dumbbells that retains qualitative aspects of fibril formation. We used rare events methods to predict separate docking and locking rate constants for the model. The rate constants were embedded in the microkinetic model, with the microkinetic model embedded in a population balance model for “bottom-up” multiscale fibril growth rate predictions. These were compared to “top-down” results using simulation data with the same model and multiscale framework to obtain maximum likelihood estimates of the separate lock and dock rate constants. We used the same procedures to extract separate docking and locking rate constants from experimental fibril growth data. Our multiscale strategy, embedding rate theories, and kinetic models in conservation laws should help to extract docking and locking rate constants from experimental data or long molecular simulations with correct units and without compromising the molecular description.
Elucidating molecular level interfacial interactions between a de novo protein and nucleated calcite with solid-state NMR
Biomineralization is the process by which organisms use biomolecules to produce hierarchically structured organic–inorganic composites. Using biology as inspiration, a protein construct (FD31) was previously designed to accelerate formation of nano-calcite with an unconventional {110} face. Here, to understand the molecular interactions essential for protein aided calcite nucleation, solid-state nuclear magnetic resonance (ssNMR) spectroscopy was used in this work to characterize the FD31–calcite interface at the atomic level. Glutamic acid side chains designed to interact directly with calcium ions on the surface were found to have dynamics on the sub-millisecond timescale, indicating possible interactions between the protein and surface waters that were not included in the original model. Dipolar ssNMR recoupling techniques also showed that the protein backbone is ∼2 Å closer to the surface than in the original docking model. Refined molecular simulations were done in the presence of explicit waters, which resulted in the protein backbone closer to the surface than in the original docking structure, providing better agreement with experiment and highlighting the important role played by water in FD31–calcite interactions. While this work provides the first experimental confirmation that FD31 interactions with calcite are localized to the surface of the protein designed to serve as a template, these studies do indicate a more dynamic binding and closer binding mode between FD31 and the nucleated surface than originally proposed. In all, this enhanced molecular insight into the FD31–calcite interface has advanced our fundamental understanding of the atomic interactions at the organic–inorganic interface and will aid in the design of biological templates for the nucleation of inorganic crystals.
DataSet for Elucidating molecular level interfacial interactions between a de novo protein and nucleated calcite with solid-state NMR
Biomineralization is the process by which organisms use biomolecules to produce hierarchically structured organic-inorganic composites. Using biology as inspiration, a protein construct (FD31) was previously designed to accelerate formation of nano-calcite with an unconventional {110} face. To understand the molecular interactions essential for protein aided calcite nucleation, solid-state nuclear magnetic resonance (ssNMR) spectroscopy was used in this work to characterize the FD31-calcite interface at the atomic level. Glutamic acid side chains designed to interact directly with calcium ions on the surface were found to have dynamics on the sub-millisecond timescale, indicating possible interactions between the protein and surface waters that were not included in the original model. Dipolar ssNMR recoupling techniques also showed that the protein backbone is ~2 Å closer to the surface than in the original docking model. Refined molecular simulations were done in the presence of explicit waters, which resulted in the protein backbone closer to the surface than in the original docking structure, providing better agreement with experiment and highlighting the important role played by water in FD31-calcite interactions. These studies provide the first experimental evidence to confirm that FD31 interactions with calcite are localized to the surface of the protein designed to serve as a template. However, these studies do indicate a more dynamic binding and closer binding mode between FD31 and the nucleated surface than originally proposed. In all, this enhanced molecular insight into the FD31-calcite interface has advanced our fundamental understanding of the atomic interactions at the organic-inorganic interface and will aid in the design of biological templates for the nucleation of inorganic crystals.
Rational Enzyme Design for Controlled Functionalization of Acetylated Xylan for Cell-Free Polymer Biosynthesis
Xylan O-acetyltransferase 1 (XOAT1) is involved in O-acetylating the backbone of hemicellulose xylan. Recent structural analysis of XOAT1 showed two unequal lobes forming a cleft that is predicted to accommodate and position xylan acceptors into proximity with the catalytic triad. Here, we used docking and molecular dynamics simulations to investigate the optimal orientation of xylan in the binding cleft of XOAT1 and identify putative key residues (Gln445 and Arg444 on Minor lobe & Asn312, Met311 and Asp403 on Major lobe) involved in substrate interactions. Site-directed mutagenesis coupled with biochemical analyses revealed the major lobe of XOAT1 is important for xylan binding. Mutation of single key residues yielded XOAT1 variants with various enzymatic efficiencies that are applicable to one-pot synthesis of xylan polymers with different degrees of O-acetylation. Taken together, our results demonstrate the effectiveness of computational modeling in guiding enzyme engineering aimed at modulating xylan and redesigning plant cell walls.
Structure and assembly of the diiron cofactor in the heme-oxygenase–like domain of the N-nitrosourea–producing enzyme SznF
In biosynthesis of the pancreatic cancer drug streptozotocin, the tridomain nonheme-iron oxygenase SznF hydroxylates N δ and N ω ' of N ω -methyl- L -arginine before oxidatively rearranging the triply modified guanidine to the N-methyl-N-nitrosourea pharmacophore. A previously published structure visualized the monoiron cofactor in the enzyme’s C-terminal cupin domain, which promotes the final rearrangement, but exhibited disorder and minimal metal occupancy in the site of the proposed diiron cofactor in the N-hydroxylating heme-oxygenase–like (HO-like) central domain. We leveraged our recent observation that the N-oxygenating µ-peroxodiiron(III/III) intermediate can form in the HO-like domain after the apo protein self-assembles its diiron(II/II) cofactor to solve structures of SznF with both of its iron cofactors bound. These structures of a biochemically validated member of the emerging heme-oxygenase–like diiron oxidase and oxygenase (HDO) superfamily with intact diiron cofactor reveal both the large-scale conformational change required to assemble the O 2 -reactive Fe 2 (II/II) complex and the structural basis for cofactor instability—a trait shared by the other validated HDOs. During cofactor (dis)assembly, a ligand-harboring core helix dynamically (un)folds. The diiron cofactor also coordinates an unanticipated Glu ligand contributed by an auxiliary helix implicated in substrate binding by docking and molecular dynamics simulations. The additional carboxylate ligand is conserved in another N-oxygenating HDO but not in two HDOs that cleave carbon–hydrogen and carbon–carbon bonds to install olefins. Among ~9,600 sequences identified bioinformatically as members of the emerging HDO superfamily, ~25% conserve this additional carboxylate residue and are thus tentatively assigned as N-oxygenases.