Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “protein dynamics”

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 OSTI2021

Quantifying the Long‐Range Coupling of Electronic Properties in Proteins with ab initio Molecular Dynamics**

Abstract The delicate interplay of covalent and non‐covalent interactions in proteins is inherently quantum mechanical and highly dynamic in nature. To directly interrogate the evolving nature of the electronic structure of proteins, we carry out 100‐ps‐scale ab initio molecular dynamics simulations of three representative small proteins with range‐separated hybrid density functional theory. We quantify the nature and length‐scale of the coupling of residue‐specific charge probability distributions in these proteins. While some nonpolar residues exhibit expectedly narrow charge distributions, most polar and charged residues exhibit broad, multimodal distributions. Even for nonpolar residues, we observe sequence‐specific deviations corresponding to charge accumulation or depletion that would be challenging to capture in a fixed charge force field. We quantify the effect of residue‐residue interactions on charge distributions first with linear cross‐correlations. We then show how additional insight can be gained from evaluating the mutual information of charge distributions. We show that a significant number of residues couple most strongly with residues that are distant in both sequence and space over a range of secondary structures including α‐helical, β‐sheet, disulfide bridging, and lasso motifs. The mutual information analysis is necessary to capture coupling between some polar and charged residues that would be otherwise missed.

Yang, Zhongyue↗
DOE OSTI2021

Full structural ensembles of intrinsically disordered proteins from unbiased molecular dynamics simulations

Abstract Molecular dynamics (MD) simulation is widely used to complement ensemble-averaged experiments of intrinsically disordered proteins (IDPs). However, MD often suffers from limitations of inaccuracy. Here, we show that enhancing the sampling using Hamiltonian replica-exchange MD (HREMD) led to unbiased and accurate ensembles, reproducing small-angle scattering and NMR chemical shift experiments, for three IDPs of varying sequence properties using two recently optimized force fields, indicating the general applicability of HREMD for IDPs. We further demonstrate that, unlike HREMD, standard MD can reproduce experimental NMR chemical shifts, but not small-angle scattering data, suggesting chemical shifts are insufficient for testing the validity of IDP ensembles. Surprisingly, we reveal that despite differences in their sequence, the inter-chain statistics of all three IDPs are similar for short contour lengths (< 10 residues). The results suggest that the major hurdle of generating an accurate unbiased ensemble for IDPs has now been largely overcome.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

All-Atom Simulations Uncover Structural and Dynamical Properties of STING Proteins in the Membrane System

Recent studies have shown that the stimulator of interferon gene (STING) protein plays a central role in the immune system by facilitating the production of type I interferons in cells. The STING signaling pathway is also a prominent activator of cancer-killing T cells that initiate a powerful adaptive immune response. Since biomolecular signaling pathways are complicated and not easily identified through traditional experiments, molecular dynamics (MD) has often been used to study structural and dynamical responses of biological pathways. In this work, we carried out MD simulations for full-length chicken and human STING (chSTING and hSTING) proteins. Specifically, we investigated ligand-bound closed (holo) and ligand-unbound open (apo) forms of STING in the membrane system by comparing their conformational and dynamical differences. Our research provides clues for understanding the mechanism of the STING signaling pathway by uncovering detailed insights for the examined systems: the residues from each chain in the binding pocket are strongly correlated to one another in the open STING structure compared with those in the closed STING structure. Ligand-bound closed STING displays ~174° rotation of the ligand-binding domain (LBD) relative to the open STING structure. The dynamical analysis of residue Cys148 located in the linker region of hSTING does not support the earlier hypothesis that Cys148 can form disulfide bonds between adjacent STING dimers. We also demonstrate that using the full-length proteins is critical, since the MD simulations of the LBD portion alone cannot properly describe the global conformational properties of STING.

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

Machine learning–driven multiscale modeling reveals lipid-dependent dynamics of RAS signaling proteins

Significance Here we present an unprecedented multiscale simulation platform that enables modeling, hypothesis generation, and discovery across biologically relevant length and time scales to predict mechanisms that can be tested experimentally. We demonstrate that our predictive simulation-experimental validation loop generates accurate insights into RAS-membrane biology. Evaluating over 100,000 correlated simulations, we show that RAS–lipid interactions are dynamic and evolving, resulting in: 1) a reordering and selection of lipid domains in realistic eight-lipid bilayers, 2) clustering of RAS into multimers correlating with specific lipid fingerprints, 3) changes in the orientation of the RAS G-domain impacting its ability to interact with effectors, and 4) demonstration that RAS–RAS G-domain interfaces are nonspecific in these putative signaling domains.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Muconic Acid Production from P. putida Using High Protein Algae Hydrolysate

The composition of algal biomass is highly dynamic, with protein, lipid, and carbohydrate contents varying in response to nutrient and environmental conditions during cultivation. Because shifts in biomass composition are often associated with reduced biomass productivity, production costs can often increase if targeting higher biomass compositional quality (enriched in carbohydrates or lipids at reduced protein content) as input for the algal biorefinery. The optimal algal biorefinery configuration is thus a function of many factors. One of the key strengths of the Combined Algal Processing (CAP) process is the versatility of feedstocks and products produced. The concept has been demonstrated with ethanol and a variety of carboxylic acids (succinic, butyric, muconic) as coproducts along with lipid upgrading to biofuel. Modification of the approaches, processes and downstream upgrading to fuels has allowed the CAP process to reduce costs and improve efficiency. Muconic acid is a high-value, potential fermentation coproduct of interest because it can be easily converted to adipic acid, a high-volume monomer for the production of nylon and other valuable consumer plastics. As such, the production of muconic acid through CAP was explored to expand the suite of products from algal biomass and to begin exploring the valorization of high protein content biomass from rapidly grown algae biomass. We have established initial performance parameters and shown that the range of substrates consumed by the muconic acid-producing microbe, Pseudomonas putida, includes at least glucose, mannose, glycerol, and lactic acid. We achieved complete utilization of these four major hydrolysate substrates achieving productivities of 0.037 (g/L/h) from Scenedesmus obliquus and 0.029 (g/L/h) from Monoraphidium minutum hydrolysates. Final titer and process yield (mol of muconic acid per mol of substrate (molP/molS)) were 0.99 g/L and 0.42 molP/molS from S. obliquus hydrolysate and 0.75 g/L and 0.23 molP/molS from M. minutum hydrolysate, respectively.

BIOMASS FUELS↗
DOE OSTI2024

Temperature-dependent iron motion in extremophile rubredoxins – no need for ‘corresponding states’

Extremophile organisms are known that can metabolize at temperatures down to – 25 °C (psychrophiles) and up to 122 °C (hyperthermophiles). Understanding viability under extreme conditions is relevant for human health, biotechnological applications, and our search for life elsewhere in the universe. Information about the stability and dynamics of proteins under environmental extremes is an important factor in this regard. Here we compare the dynamics of small Fe-S proteins – rubredoxins – from psychrophilic and hyperthermophilic microorganisms, using three different nuclear techniques as well as molecular dynamics calculations to quantify motion at the Fe site. The theory of ‘corresponding states’ posits that homologous proteins from different extremophiles have comparable flexibilities at the optimum growth temperatures of their respective organisms. Although ‘corresponding states’ would predict greater flexibility for rubredoxins that operate at low temperatures, we find that from 4 to 300 K, the dynamics of the Fe sites in these homologous proteins are essentially equivalent.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Solution characterization of the dynamic conjugative entry exclusion protein TraG

The R100 plasmid and the secretion system it encodes are representative of F-like conjugative type IV secretion systems for the transmission of mobile DNA elements in gram-negative bacteria, serving as a major contributor to the spread of antibiotic resistance in bacterial pathogens. The TraG protein of F-like systems consists of a membrane-bound N-terminal domain and a periplasmic C-terminal domain, denoted TraG*. TraG* is essential in preventing redundant DNA transfer through a process termed entry exclusion. In the donor cell, it interacts with TraN to facilitate mating pair stabilization; however, if a mating pore forms between bacteria with identical plasmids, TraG* interacts with its cognate TraS in the inner membrane of the recipient bacterium to prevent redundant donor–donor conjugation. Structural studies of TraG* from the R100 plasmid have revealed the presence of a dynamic region between the N- and C-terminal domains of TraG. Thermofluor, circular dichroism, collision-induced unfolding–mass spectrometry, and size exclusion chromatography linked to multiangle light scattering and small angle x-ray scattering experiments indicated an N-terminal truncation mutant displayed higher stability and less disordered content relative to full-length TraG*. The 45 N-terminal residues of TraG* are hypothesized to serve as part of a flexible linker between the two independently functioning domains.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Acetylcholinesterase: Structure, dynamics, and interactions with organophosphorus compounds

Acetylcholinesterase (AChE) is an enzyme that hydrolyzes the neurotransmitter acetylcholine (ACh), removing it from the synaptic cleft after the transmission of an electrical signal, making it an essential component of chemical neurotransmission. AChE is a serine hydrolase, containing a catalytic triad of Ser/His/Glu. AChE is a prime target for pharmaceuticals treating a variety of neurological disorders. It is also the target of synthetic organophosphorus (OP) compounds that have been used as pesticides and chemical warfare agents. OP compounds contain a potent leaving group, such as fluorine, and act by forming a covalent adduct with the catalytic serine of the AChE active site. A wealth of structural information is available for AChE, including over 300 structures, including a subset of structures in complex with drugs as well as OP compounds. This review will highlight the interactions between OP compounds and AChE from a structural and computational perspective, with a discussion of access to the active site, as well as side reactions that lead to dealkylation of the OP-catalytic serine adduct, a process known as aging. We conclude that while the majority of the conformational changes needed to accommodate the OP compounds are localized to the acyl loop in the crystal structures, molecular dynamics simulations highlight the potential for a far more dynamic enzyme.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2004

Step and Kink Dynamics in Inorganic and Protein Crystallization

Behavior of low-kink-density steps in solution growth and consequences for general understanding of spiral crystal growth processes will be overviewed. Also, influence of turbulence on step bunching and possibility to diminish this bunching will be presented.

Chernov, A. A.↗
DOE OSTI2021

Modeling protein structures from predicted contacts with modern molecular dynamics potentials: accuracy, sensitivity, and refinement

Protein structure prediction has become increasingly popular and successful in recent years. An essential step for fragment-free, template-free methods is the generation of a final three-dimensional protein model from a set of predicted amino acid contacts that are often described by interresidue pairwise atomic distances. Here we explore the use of modern, open-source molecular dynamics (MD) engines, which have been continually developed over the last three decades with all-atom Hamiltonians to model biomolecular structure and dynamics, to generate accurate protein structures starting from a set of inferred pairwise distances. Additionally, the ability of MD empirical physical potentials to correct inaccuracies in the predicted geometries is tested. We rigorously characterize the effect of modeling parameters on results, the effect of different amounts of error in the predicted distances on the final structures, and test the ability of post-processing analysis to sort the best models out of a set of statistical replicas. We find that with exact distances and with noisy distances, the method can produce excellent structural models, and that the molecular dynamics force field seems to help correct errors in distance predictions, resisting the effects of applied noise.

Davidson, Russ↗
DOE OSTI2022

Joint neutron/molecular dynamics vibrational spectroscopy reveals softening of HIV-1 protease upon binding of a tight inhibitor

Biomacromolecules are inherently dynamic, and their dynamics are interwoven into function. The fast collective vibrational dynamics in proteins occurs in the low picosecond timescale corresponding to frequencies of ~5-50 cm -1 . This sub-to-low THz frequency regime covers the low-amplitude collective breathing motions of a whole protein and vibrations of the constituent secondary structure elements, such as α-helices, β-sheets and loops. We have used inelastic neutron scattering experiments in combination with molecular dynamics simulations to demonstrate the vibrational dynamics softening of HIV-1 protease, a target of HIV/AIDS antivirals, upon binding of a tight clinical inhibitor darunavir. Changes in the vibrational density of states of matching structural elements in the two monomers of the homodimeric protein are not identical, indicating asymmetric effect of darunavir on the vibrational dynamics. Three of the 11 major secondary structure elements contribute over 40% to the overall changes in the vibrational density of states upon darunavir binding. Molecular dynamics simulations informed by experiments allowed us to estimate that the altered vibrational dynamics of the protease would contribute -3.6 kcal/mol -1 at 300 K, or 25%, to the free energy of darunavir binding. As HIV-1 protease drug resistance remains a concern, our results open a new avenue to help establish a direct quantitative link between protein vibrational dynamics and drug resistance.

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

Revealing functional insights into ER proteostasis through proteomics and interactomics

The endoplasmic reticulum (ER), responsible for processing approximately one-third of the human proteome including most secreted and membrane proteins, plays a pivotal role in protein homeostasis (proteostasis). Dysregulation of ER proteostasis has been implicated in a number of disease states. As such, continued efforts are directed at elucidating mechanisms of ER protein quality control which are mediated by transient and dynamic protein-protein interactions with molecular chaperones, co-chaperones, protein folding and trafficking factors that take place in and around the ER. Technological advances in mass spectrometry have played a pivotal role in characterizing and understanding these protein-protein interactions that dictate protein quality control mechanisms. Here, we highlight the recent progress from mass spectrometry-based investigation of ER protein quality control in revealing the topological arrangement of the proteostasis network, stress response mechanisms that adjust the ER proteostasis capacity, and disease specific changes in proteostasis network engagement. We close by providing a brief outlook on underexplored areas of ER proteostasis where mass spectrometry is a tool uniquely primed to further expand our understanding of the regulation and coordination of protein quality control processes in diverse diseases.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Reproducibility of protein x-ray diffuse scattering and potential utility for modeling atomic displacement parameters

Protein structure and dynamics can be probed using x-ray crystallography. Whereas the Bragg peaks are only sensitive to the average unit-cell electron density, the signal between the Bragg peaks—diffuse scattering—is sensitive to spatial correlations in electron-density variations. Although diffuse scattering contains valuable information about protein dynamics, the diffuse signal is more difficult to isolate from the background compared to the Bragg signal, and the reproducibility of diffuse signal is not yet well understood. We present a systematic study of the reproducibility of diffuse scattering from isocyanide hydratase in three different protein forms. Both replicate diffuse datasets and datasets obtained from different mutants were similar in pairwise comparisons (Pearson correlation coefficient ≥0.8). The data were processed in a manner inspired by previously published methods using custom software with modular design, enabling us to perform an analysis of various data processing choices to determine how to obtain the highest quality data as assessed using unbiased measures of symmetry and reproducibility. The diffuse data were then used to characterize atomic mobility using a liquid-like motions (LLM) model. This characterization was able to discriminate between distinct anisotropic atomic displacement parameter (ADP) models arising from different anisotropic scaling choices that agreed comparably with the Bragg data. Our results emphasize the importance of data reproducibility as a model-free measure of diffuse data quality, illustrate the ability of LLM analysis of diffuse scattering to select among alternative ADP models, and offer insights into the design of successful diffuse scattering experiments.

59 BASIC BIOLOGICAL SCIENCES↗