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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Inter-Subunit Dynamics Controls Tunnel Formation During the Oxygenation Process in Hemocyanin Hexamers

Hemocyanin from horseshoe crab in its active form is a homo-hexameric protein. It exists in open and closed conformations when transitioning between deoxygenated and oxygenated states. Here, we present a detailed dynamic atomistic investigation of the oxygenated and deoxygenated states of the hexameric hemocyanin using explicit solvent molecular dynamics simulations. We focus on the variation in solvent cavities and the formation of tunnels in the two conformational states. By employing principal component analysis and CVAE-based deep learning, we are able to differentiate between the dynamics of the deoxy- and oxygenated states of hemocyanin. Finally, our results identify the deoxygenated open conformation, which adopts a stable, closed conformation after the oxygenation process.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2022

Free Energy Landscapes for Elucidating the Structural Consequences of Exon-20 mutations on the ErbB Family of Protein Kinases

The ErbB family of protein kinases plays an important role in major cellular functions and consequently mutations in the functional regions of these proteins are implicated in several types of cancer growths. To envision rational design of small molecule drugs that target the diseased proteins it is important to quantify the structural effects of the mutations, as some of these mutants render the protein resistant to tyrosine kinase inhibitors (TKIs). Herein we use advanced sampling techniques and long-timescale molecular dynamics simulations to predict the effect of major exon 20 mutations on the ErbB family, specifically EGFR and HER2 proteins. Exon 20 mutations have been clinically known to induce TKI resistance, though the mechanisms of such an effect is poorly understood. By mapping out the free energy landscape of the mutants and comparing them against the wild-type, we elucidate the structural differences in the binding pocket region that alter the nature of drug-protein interactions. We believe that these insights will play a pivotal role in developing small molecule drugs that overcome the TKI resistance.

Ashwin Ravichandran↗
DOE OSTI2025

Dynamics and lipid membrane coupling of the RAS-RAF complex revealed via multiscale simulations

To gain molecular and mechanistic insights into initiation of the RAS-RAF signaling cascade, we developed and used a combination of multiscale simulation and experimental approaches. The influence and impact of the membrane on RAS and RAF proteins is a factor we are just beginning to understand and appreciate in more detail. Molecular simulation is an ideal methodology to further study this complicated relationship between the membrane and associated proteins. Our previous work using Multiscale Machine-learned Modeling Infrastructure investigated different lipid compositions solely around the KRAS4b protein and the interplay between protein behavior and these membrane environments. Multiscale Machine-learned Modeling Infrastructure uses machine learning to couple adjacent simulation scales and has been efficiently scaled across some of the world’s largest high-performance computers. Recently, we have expanded this multiresolution framework to include the all-atom simulation scale and to incorporate the RAF RBDCRD domains. Here, we present the overall analysis results from this new simulation campaign comprising a mixture of RAS and RAF RBDCRD proteins. Approximately 35,000 coarse-grained and 10,000 all-atom molecular dynamics simulations were completed, sampled from a variety of protein/lipid composition configurations that were generated from a micron-scale continuum simulation containing hundreds of copies of the proteins. Our studies suggest that orientations of the RAS-RBDCRD complex on the membrane occupy distinct configurational states, and the spatial patterns of lipid arrangements around these different protein states are unique to each state. The extent and size of lipid “fingerprints” imposed on the membrane by the RAS-RBDCRD protein complex are significantly larger than observed for just the RAS protein on its own. These protein complexes strongly associate, but we do not observe statistically significant preferred protein-protein orientations. These observations indicate that spatial colocalization of RAS-RBDCRD proteins in the same vicinity may be assisted by specific membrane environments, acting to increase the probability of signaling complex formation.

Carpenter, Timothy S. [Lawrence Livermore National↗
DOE OSTI2020

Roles for Countercharge in the Voltage Sensor Domain of Ion Channels

Voltage-gated ion channels share a common structure typified by peripheral, voltage sensor domains. Their S4 segments respond to alteration in membrane potential with translocation coupled to ion permeation through a central pore domain. The mechanisms of gating in these channels have been intensely studied using pioneering methods such as measurement of charge displacement across a membrane, sequencing of genes coding for voltage-gated ion channels, and the development of all-atom molecular dynamics simulations using structural information from prokaryotic and eukaryotic channel proteins. One aspect of this work has been the description of the role of conserved negative countercharges in S1, S2, and S3 transmembrane segments to promote sequential salt-bridge formation with positively charged residues in S4 segments. These interactions facilitate S4 translocation through the lipid bilayer. In this review, we describe functional and computational work investigating the role of these countercharges in S4 translocation, voltage sensor domain hydration, and in diseases resulting from countercharge mutations.

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

Domain coupling in activation of a family C GPCR

Abstract The G protein-coupled metabotropic glutamate receptors form homodimers and heterodimers with highly diverse responses to glutamate and varying physiological functions. We employ molecular dynamics, single-molecule spectroscopy and hydrogen–deuterium exchange to dissect the activation pathway triggered by glutamate. We find that activation entails multiple loosely coupled steps, including formation of an agonist-bound, pre-active intermediate whose transition to active conformations forms dimerization interface contacts that set efficacy. The agonist-bound receptor populates at least two additional intermediates en route to G protein-coupling conformations. Sequential transitions into these states act as ‘gates’, which attenuate the effects of glutamate. Thus, the agonist-bound receptor is remarkably dynamic, with low occupancy of G protein-coupling conformations, providing considerable headroom for modulation by allosteric ligands. Sequence variation within the dimerization interface, as well as altered conformational coupling in receptor heterodimers, may contribute to precise decoding of glutamate signals over broad spatial and temporal scales.

Biochemistry & Molecular Biology↗
DOE OSTI2023

Atomistic structure of the SARS-CoV-2 pseudoknot in solution from SAXS-driven molecular dynamics

SARS-CoV-2 depends on -1 programmed ribosomal frameshifting (-1 PRF) to express proteins essential for its replication. The RNA pseudoknot stimulating -1 PRF is thus an attractive drug target. However, the structural models of this pseudoknot obtained from cryo-EM and crystallography differ in some important features, leaving the pseudoknot structure unclear. We measured the solution structure of the pseudoknot using small-angle X-ray scattering (SAXS). The measured profile did not agree with profiles computed from the previously solved structures. Beginning with each of these solved structures, we used the SAXS data to direct all atom molecular dynamics (MD) simulations to improve the agreement in profiles. In all cases, this refinement resulted in a bent conformation that more closely resembled the cryo-EM structures than the crystal structure. Applying the same approach to a point mutant abolishing -1 PRF revealed a notably more bent structure with reoriented helices. This work clarifies the dynamic structures of the SARS-CoV-2 pseudoknot in solution.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2003

Structure, Function, Self-Assembly and Origin of Simple Membrane Proteins

Integral membrane proteins perform such essential cellular functions as transport of ions, nutrients and waste products across cell walls, transduction of environmental signals, regulation of cell fusion, recognition of other cells, energy capture and its conversion into high-energy compounds. In fact, 30-40% of genes in modem organisms codes for membrane proteins. Although contemporary membrane proteins or their functional assemblies can be quite complex, their transmembrane fragments are usually remarkably simple. The most common structural motif for these fragments is a bundle of alpha-helices, but occasionally it could be a beta-barrel. In a series of molecular dynamics computer simulations we investigated self-organizing properties of simple membrane proteins based on these structural motifs. Specifically, we studied folding and insertion into membranes of short, nonpolar or amphiphatic peptides. We also investigated glycophorin A, a peptide that forms sequence-specific dimers, and a transmembrane aggregate of four identical alpha-helices that forms an efficient and selective voltage-gated proton channel was investigated. Many peptides are attracted to water-membrane interfaces. Once at the interface, nonpolar peptides spontaneously fold to a-helices. Whenever the sequence permits, peptides that contain both polar and nonpolar amino also adopt helical structures, in which polar and nonpolar amino acid side chains are immersed in water and membrane, respectively. Specific identity of side chains is less important. Helical peptides at the interface could insert into the membrane and adopt a transmembrane conformation. However, insertion of a single helix is unfavorable because polar groups in the peptide become completely dehydrated upon insertion. The unfavorable free energy of insertion can be regained by spontaneous association of peptides in the membrane. The first step in this process is the formation of dimers, although the most common are aggregates of 4-7 helices. The helices could arrange themselves such that they formed pores capable of transporting ions and small molecules across membranes. Stability of transmembrane aggregates of simple proteins is often only marginal and, therefore, it can be regulated by environmental signals or small sequence modifications in the region of interhelical interactions. A key step in the earliest evolution of membrane proteins was the emergence of selectivity for specific substrates. Many channels could become selective if one or only a few properly chosen amino acids are properly placed along the channel, acting as filters or gates. This is a convenient evolutionary solution because it does not require imposing conditions on the whole sequence.

Pohorille, Andrew↗
DOE OSTI2026

Activation dynamics of a water-soluble human mu-opioid receptor

The mu-opioid receptor (MOR), a class A G protein-coupled receptor mediates opioid analgesia and remains a central target for pain therapeutics. While crystal structures of MOR exist, they provide limited insight into the receptor’s dynamic conformational landscape underlying function. Here, we engineered a thermostable water-soluble MOR variant (wsMOR) that retains native-like ligand-binding and activation dynamics. This variant enables high-yield production and detailed solution-phase structural studies that are challenging with membrane-embedded MOR, providing a valuable tool for studying receptor activation and aqueous-phase drug screening. Using a combined computational and experimental approach, we performed long-timescale all-atom molecular dynamics simulations together with neutron scattering and single-molecule FRET, revealing a structurally stable receptor with a diverse ensemble of conformations at different temporal resolutions. In the ligand-free state, wsMOR displayed high conformational flexibility, which decreased upon agonist binding, particularly in transmembrane helix 6, a hallmark of G protein-coupled receptor activation. Positive allosteric modulation and G protein binding further stabilized active-like states. These findings highlight wsMOR’s conformational plasticity across picosecond to millisecond timescales and provide a foundation for structure-guided development of next-generation opioid ligands with improved efficacy and safety.

E, Agyemang [University of Tennessee Knoxville]↗
DOE OSTI2025

3D pattern formation of a protein–membrane suspension

Many essential cellular processes, including cell division and the establishment of cell polarity during embryogenesis, are regulated by pattern-forming proteins. These proteins often need to bind to a substrate, such as the cell membrane, onto which they interact and form two-dimensional (2D) patterns. It is unclear how the membrane’s continuity and dimensionality impact pattern formation. Here, we address this gap using the MinDE system, a prototypical example of pattern-forming membrane proteins. We show that when the lipid substrate is fragmented into submicrometer-sized diffusive liposomes, adenosine triphosphate-driven protein–protein interactions generate three-dimensional (3D) spatially extended patterns, despite the complete loss of membrane continuity. Remarkably, these 3D patterns emerge at scales four orders of magnitude larger than the individual liposomes. By systematically varying protein concentration, liposome size, and density, we observed and characterized a variety of 3D dynamical patterns not seen on continuous 2D membranes, including traveling waves, dynamical spirals, and a coexistence phase. Simulations and linear stability analysis of a coarse-grained model revealed that the physical properties of the dispersed membrane effectively rescale both the protein–membrane binding rates and diffusion, two key parameters governing pattern formation and wavelength selection. These findings highlight the robustness of Min’s pattern-forming ability, suggesting that protein–membrane suspensions could serve as an adaptable template for studying out-of-equilibrium self-organization in 3D, beyond in vivo contexts.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A generalized Flory-Stockmayer kinetic theory of connectivity percolation and rigidity percolation of cytoskeletal networks

Actin networks are essential for living cells to move, reproduce, and sense their environments. The dynamic and rheological behavior of actin networks is modulated by actin-binding proteins such as α-actinin, Arp2/3, and myosin. There is experimental evidence that actin-binding proteins modulate the cooperation of myosin motors by connecting the actin network. In this work, we present an analytical mean field model, using the Flory-Stockmayer theory of gelation, to understand how different actin-binding proteins change the connectivity of the actin filaments as the networks are formed. We follow the kinetics of the networks and estimate the concentrations of actin-binding proteins that are needed to reach connectivity percolation as well as to reach rigidity percolation. We find that Arp2/3 increases the actomyosin connectivity in the network in a non-monotonic way. We also describe how changing the connectivity of actomyosin networks modulates the ability of motors to exert forces, leading to three possible phases of the networks with distinctive dynamical characteristics: a sol phase, a gel phase, and an active phase. Thus, changes in the concentration and activity of actin-binding proteins in cells lead to a phase transition of the actin network, allowing the cells to perform active contraction and change their rheological properties.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Insights into the substrate specificity, structure, and dynamics of plant histidinol-phosphate aminotransferase (HISN6)

Histidinol-phosphate aminotransferase is the sixth protein (hence HISN6) in the histidine biosynthetic pathway in plants. HISN6 is a pyridoxal 5'-phosphate (PLP)-dependent enzyme that catalyzes the reversible conversion of imidazole acetol phosphate into L-histidinol phosphate (HOLP). Here, we show that plant HISN6 enzymes are closely related to the orthologs from Chloroflexota. The studied example, HISN6 from Medicago truncatula (MtHISN6), exhibits a surprisingly high affinity for HOLP, which is much higher than reported for bacterial homologs. Moreover, unlike the latter, MtHISN6 does not transaminate phenylalanine. High-resolution crystal structures of MtHISN6 in the open and closed states, as well as the complex with HOLP and the apo structure without PLP, bring new insights into the enzyme dynamics, pointing at a particular role of a string-like fragment that oscillates near the active site and participates in the HOLP binding. When MtHISN6 is compared to bacterial orthologs with known structures, significant differences arise in or near the string region. The high affinity of MtHISN6 appears linked to the particularly tight active site cavity. Finally, a virtual screening against a library of over 1.3 mln compounds revealed three sites in the MtHISN6 structure with the potential to bind small molecules. Such compounds could be developed into herbicides inhibiting plant HISN6 enzymes absent in animals, which makes them a potential target for weed control agents.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Design of facilitated dissociation enables timing of cytokine signalling

Protein design has focused on the design of ground states, ensuring that they are sufficiently low energy to be highly populated. Designing the kinetics and dynamics of a system requires, in addition, the design of excited states that are traversed in transitions from one low-lying state to another. This is a challenging task because such states must be sufficiently strained to be poorly populated, but not so strained that they are not populated at all, and because protein design methods have focused on generating near-ideal structures. Here we describe a general approach for designing systems that use an induced-fit power stroke to generate a structurally frustrated and strained excited state, allosterically driving protein complex dissociation. X-ray crystallography, double electron–electron resonance spectroscopy and kinetic binding measurements show that incorporating excited states enables the design of effector-induced increases in dissociation rates as high as 5,700-fold. We highlight the power of this approach by designing rapid biosensors, kinetically controlled circuits and cytokine mimics that can be dissociated from their receptors within seconds, enabling dissection of the temporal dynamics of interleukin-2 signalling.

deformation dynamics↗
NASA NTRS1999

Dynamic Light Scattering of Diabetic Vitreopathy

Diabetes induces pathology throughout the body via nonenzymatic glycation of proteins. Vitreous, which is replete with type 11 collagen, undergoes significant changes in diabetes. The resultant diabetic vitreopathy plays an important role in diabetic retinopathy. Detecting these molecular changes could provide insight into diabetic eye disease as well as molecular effects elsewhere in the body. Human eyes were obtained at autopsy and studied in the fresh, unfixed state. Sclera, choroid, and retina were dissected off the vitreous for dark-field slit microscopy and dynamic light scattering (DLS). For the former, the entire vitreous was exposed. For the latter, only a window at the equator was dissected in some specimens, and the anterior segment was removed leaving the posterior lens capsule intact in others. DLS was performed to determine particle sizes at multiple sites 0.5 mm apart, spanning the globe at the equator (window dissections) and along the antero-posterior axis. Dark-field slit microscopy in diabetic subjects detected findings typical of age-related vitreous degeneration, but at much younger ages than nondiabetic controls. Noninvasive DLS measurements found a greater heterogeneity and larger particle sizes in vitreous of subjects with diabetes as compared to age-matched controls. DLS can detect and quantify the early molecular effects that cause vitreous collagen fibrils to cross-link and aggregate. This could provide valuable insight into ocular and systemic effects of hyperglycemia, because the molecular changes in diabetic vitreopathy could serve as an index of such effects throughout the body. In addition to the diagnostic implications, this methodology could provide a rapid, reproducible way to monitor the response to therapy with novel agents intended to prevent the complications of diabetes on a molecular level.

Sebag, J.↗
DOE OSTI2025

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↗
DOE OSTI2022

Mass spectrometry-based technologies for probing the 3D world of plant proteins

Abstract Over the past two decades, mass spectrometric (MS)-based proteomics technologies have facilitated the study of signaling pathways throughout biology. Nowhere is this needed more than in plants, where an evolutionary history of genome duplications has resulted in large gene families involved in posttranslational modifications and regulatory pathways. For example, at least 5% of the Arabidopsis thaliana genome (ca. 1,200 genes) encodes protein kinases and protein phosphatases that regulate nearly all aspects of plant growth and development. MS-based technologies that quantify covalent changes in the side-chain of amino acids are critically important, but they only address one piece of the puzzle. A more crucially important mechanistic question is how noncovalent interactions—which are more difficult to study—dynamically regulate the proteome’s 3D structure. The advent of improvements in protein 3D technologies such as cryo-electron microscopy, nuclear magnetic resonance, and X-ray crystallography has allowed considerable progress to be made at this level, but these methods are typically limited to analyzing proteins, which can be expressed and purified in milligram quantities. Newly emerging MS-based technologies have recently been developed for studying the 3D structure of proteins. Importantly, these methods do not require protein samples to be purified and require smaller amounts of sample, opening the wider proteome for structural analysis in complex mixtures, crude lysates, and even in intact cells. These MS-based methods include covalent labeling, crosslinking, thermal proteome profiling, and limited proteolysis, all of which can be leveraged by established MS workflows, as well as newly emerging methods capable of analyzing intact macromolecules and the complexes they form. In this review, we discuss these recent innovations in MS-based “structural” proteomics to provide readers with an understanding of the opportunities they offer and the remaining challenges for understanding the molecular underpinnings of plant structure and function.

Plant Sciences↗
NASA NTRS2004

Antigenic Protein In Microgravity-Grown Human Mixed Mullerian Tumor (LN1) Cells Preserved In RNA Stabilizing Agent

Cells treated with RNAlater(TradeMark) have previously been shown to contain antigenic proteins that can be visualized using Western blot analysis. These proteins seem to be stable for several months when stored in RNA stabilizer at 4 C. Antigenic protein can be recovered from cells that have been processed using an Ambion RNAqueous(Registered TradeMark) kit to remove RNA. In this set of experiments, human mixed Mullerian tumor (LNI) cells grown on the International Space Station during Expedition 3 were examined for antigenic stability after removal of RNA. The cells were stored for three months in RNAlater(TradeMark) and RNA was extracted. The RNA filtrate containing the protein was precipitated, washed, and suspended in buffer containing sodium dodecyl sulfate (SDS). Samples containing equal concentrations of protein were loaded onto SDS-polyacrylamide gels. Proteins were separated by electrophoresis and transferred by Western blot to polyvinylidene fluoride (PVDF) membrane. The Western blots were stained with an enhanced chemiluminescent ECL(Registered Trademark) Plus detection kit (Amersham) and scanned using a Storm 840 gel image analyzer (Amersham, Molecular Dynamics). ImageQuant(Registered TradeMark) software was used to quantify the densities of the protein bands. The ground control and flight LN1 cell samples showed a similar staining pattern over time with antibodies to vimentin, glyceraldehyde-3-phosphate dehydrogenase, and epithelial membrane antigens.

Hammond, Dianne K.↗
DOE OSTI2022

Characterization of Heparin’s Conformational Ensemble by Molecular Dynamics Simulations and Nuclear Magnetic Resonance Spectroscopy

Heparin is a highly charged, polysulfated polysaccharide and serves as an anticoagulant. Heparin binds to multiple proteins throughout the body, suggesting a large range of potential therapeutic applications. Although its function has been characterized in multiple physiological contexts, heparin’s solution conformational dynamics and structure–function relationships are not fully understood. Molecular dynamics (MD) simulations facilitate the analysis of a molecule’s underlying conformational ensemble, which then provides important information necessary for understanding structure–function relationships. However, for MD simulations to afford meaningful results, they must both provide adequate sampling and accurately represent the energy properties of a molecule. The aim of this study is to compare heparin’s conformational ensemble using two well-developed force fields for carbohydrates, known as GLYCAM06 and CHARMM36, using replica exchange molecular dynamics (REMD) simulations, and to validate these results with NMR experiments. In this work, the anticoagulant sequence, an ultra-low-molecular-weight heparin, known as Arixtra (fondaparinux, sodium), was simulated with both parameter sets. The results suggest that GLYCAM06 matches experimental nuclear magnetic resonance three-bond J-coupling values measured for Arixtra better than CHARMM36. In addition, NOESY and ROESY experiments suggest that Arixtra is very flexible in the sub-millisecond time scale and does not adopt a unique structure at 25 C. Moreover, GLYCAM06 affords a much more dynamic conformational ensemble for Arixtra than CHARMM36.

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

Beneath the surface: Unsolved questions in soil virus ecology

Soil virus ecology is an exciting but still nascent field of research in soil microbiology. While there has been a recent surge in soil virus research studies, many fundamental questions remain unanswered, and a range of technical and bioinformatic challenges need to be overcome. In this perspective article, we present a series of key questions that highlight fruitful research areas for ongoing and future efforts. These include describing the challenges involved in understanding soil viral abundance and activity, spatiotemporal dynamics, life strategy prevalence, virus-mediated biogeochemical impacts, viral protein function, host prediction, and soil RNA virus discovery. In the near term, combining approaches (e.g., cultivation-based, meta-omics, biogeochemical, experimental, and bioinformatic) will be key to assessing the ecological and biogeochemical impacts of soil viruses from the microscopic to the field and global scales. Still, we stress that results must be tempered by current methodological limitations and highlight knowledge gaps that are most pressing to fill via new methods or measurements, such as the prevalence of different viral replication strategies across soils, the fate of microbial necromass carbon after viral lysis, the frequency of virus-host encounters that do not lead to successful infections yet could be bioinformatically mistaken as infections, and the diversity and ecological impacts of RNA viruses in soil.

59 BASIC BIOLOGICAL SCIENCES↗