Search NASA⌕ Search

SEARCH · Search NASA

Results for “protein conformational 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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Steering Molecular Dynamics Simulations of Membrane-Associated Proteins with Neutron Reflection Results

Herein, we present a novel method to incorporate structural results from surface-sensitive scattering, such as X-ray or neutron reflectometry, into molecular dynamics simulations. While reflectometry techniques generally provide a means to determine the molecular-scale structures of organized interfacial films, they were recently shown to offer the capability to characterize the structures of protein-membrane complexes supported by a solid substrate. One-dimensional information inherent in the experimental results is used in the form of component volume occupancy (CVO) profiles, which describe the distribution of molecular components within an interfacial architecture, to construct real-space constraints in the form of a biasing potential for the simulation that vanishes when the simulated and experimental profiles agree. This approach improves the correspondence between simulation and experiment, as shown in the re-evaluation of an neutron-reflection-derived structure which was approximated by an independent molecular dynamics simulation in earlier work, and it also leads to faster equilibration of ensemble structures. We further show that time averaging the CVO profile that develops in the simulation while biasing with this approach permits fluctuations about the average that are necessary for conformational exploration of the system. This method is particularly valuable for studies of proteins at interfaces that contain disordered regions since the conformation of such regions is difficult to judge from the analysis of one-dimensional experimental profiles and may take prohibitively long to equilibrate in simulations.

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

An interoperable implementation of collective‐variable based enhanced sampling methods in extended phase space within the OpenMM package

Collective variable (CV)-based enhanced sampling techniques are widely used today for accelerating barrier-crossing events in molecular simulations. A class of these methods, which includes temperature accelerated molecular dynamics (TAMD)/driven-adiabatic free energy dynamics (d-AFED), unified free energy dynamics (UFED), and temperature accelerated sliced sampling (TASS), uses an extended variable formalism to achieve quick exploration of conformational space. These techniques are powerful, as they enhance the sampling of a large number of CVs simultaneously compared to other techniques. Extended variables are kept at a much higher temperature than the physical temperature by ensuring adiabatic separation between the extended and physical subsystems and employing rigorous thermostatting. Here, in this work, we present a computational platform to perform extended phase space enhanced sampling simulations using the open-source molecular dynamics engine OpenMM. The implementation allows users to have interoperability of sampling techniques, as well as employ state-of-the-art thermostats and multiple time-stepping. This work also presents protocols for determining the critical parameters and procedures for reconstructing high-dimensional free energy surfaces. As a demonstration, we present simulation results on the high dimensional conformational landscapes of the alanine tripeptide in vacuo, tetra-N-methylglycine (tetra-sarcosine) peptoid in implicit solvent, and the Trp-cage mini protein in explicit water.

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

Structural characterization of an intrinsically disordered protein complex using integrated small-angle neutron scattering and computing

Characterizing structural ensembles of intrinsically disordered proteins (IDPs) and intrinsically disordered regions (IDRs) of proteins is essential for studying structure–function relationships. Due to the different neutron scattering lengths of hydrogen and deuterium, selective labeling and contrast matching in small-angle neutron scattering (SANS) becomes an effective tool to study dynamic structures of disordered systems. However, experimental timescales typically capture measurements averaged over multiple conformations, leaving complex SANS data for disentanglement. We hereby demonstrate an integrated method to elucidate the structural ensemble of a complex formed by two IDRs. We use data from both full contrast and contrast matching with residue-specific deuterium labeling SANS experiments, microsecond all-atom molecular dynamics (MD) simulations with four molecular mechanics force fields, and an autoencoder-based deep learning (DL) algorithm. From our combined approach, we show that selective deuteration provides additional information that helps characterize structural ensembles. We find that among the four force fields, a99SB-disp and CHARMM36m show the strongest agreement with SANS and NMR experiments. In addition, our DL algorithm not only complements conventional structural analysis methods but also successfully differentiates NMR and MD structures which are indistinguishable on the free energy surface. Finally, we present an ensemble that describes experimental SANS and NMR data better than MD ensembles generated by one single force field and reveal three clusters of distinct conformations. Our results demonstrate a new integrated approach for characterizing structural ensembles of IDPs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Secondary structure determines electron transport in peptides

Proteins play a key role in biological electron transport, but the structure–function relationships governing the electronic properties of peptides are not fully understood. Despite recent progress, understanding the link between peptide conformational flexibility, hierarchical structures, and electron transport pathways has been challenging. Here, we use single-molecule experiments, molecular dynamics (MD) simulations, nonequilibrium Green’s function-density functional theory (NEGF-DFT), and unsupervised machine learning to understand the role of secondary structure on electron transport in peptides. Our results reveal a two-state molecular conductance behavior for peptides across several different amino acid sequences. MD simulations and Gaussian mixture modeling are used to show that this two-state molecular conductance behavior arises due to the conformational flexibility of peptide backbones, with a high-conductance state arising due to a more defined secondary structure (beta turn or 3 10 helices) and a low-conductance state occurring for extended peptide structures. These results highlight the importance of helical conformations on electron transport in peptides. Conformer selection for the peptide structures is rationalized using principal component analysis of intramolecular hydrogen bonding distances along peptide backbones. Molecular conformations from MD simulations are used to model charge transport in NEGF-DFT calculations, and the results are in reasonable qualitative agreement with experiments. Projected density of states calculations and molecular orbital visualizations are further used to understand the role of amino acid side chains on transport. Overall, our results show that secondary structure plays a key role in electron transport in peptides, which provides broad avenues for understanding the electronic properties of proteins.

Science & Technology - Other Topics↗
DOE OSTI2023

Structure of LARP7 Protein p65–telomerase RNA Complex in Telomerase Revealed by Cryo-EM and NMR

La-related protein 7 (LARP7) are a family of RNA chaperones that protect the 3'-end of RNA and are components of specific ribonucleoprotein complexes (RNP). Here, in Tetrahymena thermophila telomerase, LARP7 protein p65 together with telomerase reverse transcriptase (TERT) and telomerase RNA (TER) form the core RNP. p65 has four known domains—N-terminal domain (NTD), La motif (LaM), RNA recognition motif 1 (RRM1), and C-terminal xRRM2. To date, only the xRRM2 and LaM and their interactions with TER have been structurally characterized. Conformational dynamics leading to low resolution in cryo-EM density maps have limited our understanding of how full-length p65 specifically recognizes and remodels TER for telomerase assembly. Here, we combined focused classification of Tetrahymena telomerase cryo-EM maps with NMR spectroscopy to determine the structure of p65–TER. Three previously unknown helices are identified, one in the otherwise intrinsically disordered NTD that binds the La module, one that extends RRM1, and another preceding xRRM2, that stabilize p65–TER interactions. The extended La module (αN, LaM and RRM1) interacts with the four 3' terminal U nucleotides, while LaM and αN additionally interact with TER pseudoknot, and LaM with stem 1 and 5' end. Our results reveal the extensive p65–TER interactions that promote TER 3'-end protection, TER folding, and core RNP assembly and stabilization. The structure of full-length p65 with TER also sheds light on the biological roles of genuine La and LARP7 proteins as RNA chaperones and core RNP components.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Mechanistic Insights into Peptide Binding and Deactivation of an Adhesion G Protein-Coupled Receptor

Adhesion G protein-coupled receptors (ADGRGs) play critical roles in the reproductive, neurological, cardiovascular, and endocrine systems. In particular, ADGRG2 plays a significant role in Ewing sarcoma cell proliferation, parathyroid cell function, and male fertility. In 2022, a cryo-EM structure was reported for the active ADGRG2 bound by an optimized peptide agonist IP15 and the Gs protein. The IP15 peptide agonist was also modified to antagonists 4PH-E and 4PH-D with mutations of the 4PH residue to Glu and Asp, respectively. However, experimental structures of inactive antagonist-bound ADGRs remain to be resolved, and the activation mechanism of ADGRs such as ADGRG2 is poorly understood. Here, we applied Gaussian accelerated molecular dynamics (GaMD) simulations to probe conformational dynamics of the agonist- and antagonist-bound ADGRG2. By performing GaMD simulations, we were able to identify important low-energy conformations of ADGRG2 in the active, intermediate, and inactive states, as well as explore the binding conformations of each peptide. Moreover, our simulations revealed critical peptide-receptor residue interactions during the deactivation of ADGRG2. In conclusion, through GaMD simulations, we uncovered mechanistic insights into peptide (agonist and antagonist) binding and deactivation of the ADGRG2. These findings will potentially facilitate rational design of new peptide modulators of ADGRG2 and other ADGRs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Designed 2D protein crystals as dynamic molecular gatekeepers for a solid-state device

The sensitivity and responsiveness of living cells to environmental changes are enabled by dynamic protein structures, inspiring efforts to construct artificial supramolecular protein assemblies. However, despite their sophisticated structures, designed protein assemblies have yet to be incorporated into macroscale devices for real-life applications. We report a 2D crystalline protein assembly of C98/E57/E66 L-rhamnulose-1-phosphate aldolase ( CEE RhuA) that selectively blocks or passes molecular species when exposed to a chemical trigger. CEE RhuA crystals are engineered via cobalt(II) coordination bonds to undergo a coherent conformational change from a closed state (pore dimensions <1 nm) to an ajar state (pore dimensions ~4 nm) when exposed to an HCN(g) trigger. When layered onto a mesoporous silicon (pSi) photonic crystal optical sensor configured to detect HCN (g) , the 2D CEE RhuA crystal layer effectively blocks interferents that would otherwise result in a false positive signal. The 2D CEE RhuA crystal layer opens in selective response to low-ppm levels of HCN (g) , allowing analyte penetration into the pSi sensor layer for detection. These findings illustrate that designed protein assemblies can function as dynamic components of solid-state devices in non-aqueous environments.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Mutations in penicillin-binding protein 2 from cephalosporin-resistant Neisseria gonorrhoeae hinder ceftriaxone acylation by restricting protein dynamics

The global incidence of the sexually transmitted disease gonorrhea is expected to rise due to the spread of Neisseria gonorrhoeae strains with decreased susceptibility to extended-spectrum cephalosporins (ESCs). ESC resistance is conferred by mosaic variants of penicillin-binding protein 2 (PBP2) that have diminished capacity to form acylated adducts with cephalosporins. To elucidate the molecular mechanisms of ESC resistance, we conducted a biochemical and high-resolution structural analysis of PBP2 variants derived from the decreased-susceptibility N. gonorrhoeae strain 35/02 and ESC-resistant strain H041. Our data reveal that mutations both lower affinity of PBP2 for ceftriaxone and restrict conformational changes that normally accompany acylation. Specifically, we observe that a G545S substitution hinders rotation of the β3 strand necessary to form the oxyanion hole for acylation and also traps ceftriaxone in a noncanonical configuration. In addition, F504L and N512Y substitutions appear to prevent bending of the β3–β4 loop that is required to contact the R1 group of ceftriaxone in the active site. Other mutations also appear to act by reducing flexibility in the protein. Overall, our findings reveal that restriction of protein dynamics in PBP2 underpins the ESC resistance of N. gonorrhoeae.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Predicting the structural basis of targeted protein degradation by integrating molecular dynamics simulations with structural mass spectrometry

Targeted protein degradation (TPD) is a promising approach in drug discovery for degrading proteins implicated in diseases. A key step in this process is the formation of a ternary complex where a heterobifunctional molecule induces proximity of an E3 ligase to a protein of interest (POI), thus facilitating ubiquitin transfer to the POI. In this work, we characterize 3 steps in the TPD process. (1) We simulate the ternary complex formation of SMARCA2 bromodomain and VHL E3 ligase by combining hydrogen-deuterium exchange mass spectrometry with weighted ensemble molecular dynamics (MD). (2) We characterize the conformational heterogeneity of the ternary complex using Hamiltonian replica exchange simulations and small-angle X-ray scattering. (3) We assess the ubiquitination of the POI in the context of the full Cullin-RING Ligase, confirming experimental ubiquitinomics results. Differences in degradation efficiency can be explained by the proximity of lysine residues on the POI relative to ubiquitin.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Effects of presenilin-1 familial Alzheimer’s disease mutations on γ-secretase activation for cleavage of amyloid precursor protein

Presenilin-1 (PS1) is the catalytic subunit of γ-secretase which cleaves within the transmembrane domain of over 150 peptide substrates. Dominant missense mutations in PS1 cause early-onset familial Alzheimer’s disease (FAD); however, the exact pathogenic mechanism remains unknown. Here we combined Gaussian accelerated molecular dynamics (GaMD) simulations and biochemical experiments to determine the effects of six representative PS1 FAD mutations (P117L, I143T, L166P, G384A, L435F, and L286V) on the enzyme-substrate interactions between γ-secretase and amyloid precursor protein (APP). Biochemical experiments showed that all six PS1 FAD mutations rendered γ-secretase less active for the endoproteolytic (ε) cleavage of APP. Distinct low-energy conformational states were identified from the free energy profiles of wildtype and PS1 FAD-mutant γ-secretase. The P117L and L286V FAD mutants could still sample the “Active” state for substrate cleavage, but with noticeably reduced conformational space compared with the wildtype. The other mutants hardly visited the “Active” state. The PS1 FAD mutants were found to reduce γ-secretase proteolytic activity by hindering APP residue L49 from proper orientation in the active site and/or disrupting the distance between the catalytic aspartates. Therefore, our findings provide mechanistic insights into how PS1 FAD mutations affect structural dynamics and enzyme-substrate interactions of γ-secretase and APP.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

SPARC: Structural properties associated with residue constraints

SPARC facilitates the generation of plausible hypotheses regarding underlying biochemical mechanisms by structurally characterizing protein sequence constraints. Such constraints appear as residues co-conserved in functionally related subgroups, as subtle pairwise correlations (i.e., direct couplings), and as correlations among these sequence features or with structural features. SPARC performs three types of analyses. First, based on pairwise sequence correlations, it estimates the biological relevance of alternative conformations and of homomeric contacts, as illustrated here for death domains. Second, it estimates the statistical significance of the correspondence between directly coupled residue pairs and interactions at heterodimeric interfaces. Third, given molecular dynamics simulated structures, it characterizes interactions among constrained residues or between such residues and ligands that: (a) are stably maintained during the simulation; (b) undergo correlated formation and/or disruption of interactions with other constrained residues; or (c) switch between alternative interactions. We illustrate this for two homohexameric complexes: the bacterial enhancer binding protein (bEBP) NtrC1, which activates transcription by remodeling RNA polymerase (RNAP) containing σ 54 , and for DnaB helicase, which opens DNA at the bacterial replication fork. Based on the NtrC1 analysis, we hypothesize possible mechanisms for inhibiting ATP hydrolysis until ADP is released from an adjacent subunit and for coupling ATP hydrolysis to restructuring of σ 54 binding loops. Based on the DnaB analysis, we hypothesize that DnaB ‘grabs’ ssDNA by flipping every fourth base and inserting it into cavities between subunits and that flipping of a DnaB-specific glutamine residue triggers ATP hydrolysis.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Quaternary structure independent folding of voltage-gated ion channel pore domain subunits

Every voltage-gated ion channel (VGIC) has a pore domain (PD) made from four subunits, each comprising an antiparallel transmembrane helix pair bridged by a loop. The extent to which PD subunit structure requires quaternary interactions is unclear. In this report we present crystal structures of a set of bacterial voltage-gated sodium channel (BacNa V ) 'pore only' proteins that reveal a surprising collection of non-canonical quaternary arrangements in which the PD tertiary structure is maintained. This context-independent structural robustness, supported by molecular dynamics simulations, indicates that VGIC-PD tertiary structure is independent of quaternary interactions. This fold occurs throughout the VGIC superfamily and in diverse transmembrane and soluble proteins. Strikingly, characterization of PD subunit-binding Fabs indicates that non-canonical quaternary PD conformations can occur in full-length VGICs. Together, our data demonstrate that the VGIC-PD is an autonomously folded unit. This property has implications for VGIC biogenesis, understanding functional states, de novo channel design, and VGIC structural origins.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Crystalline Peptoid Nanofibers with a Single-Unit Cell Cross Section

Ultranarrow crystalline one-dimensional nanostructures formed from soft materials facilitate precise structural control in nanomaterial design, which is essential for biomedicine and nanotechnology applications. Systematic control of their hierarchical structure is challenging due to the complexities of simultaneously manipulating multiple noncovalent interactions at such small scales. We employed a polypeptoid crystal motif as a supramolecular synthon to engineer ultranarrow crystalline nanofibers constrained to a single lattice axis by incorporating a single ionizable side chain into the hydrophobic core of a nanosheet-forming peptoid. Cryogenic transmission electron microscopy of the nanofibers revealed detailed molecular arrangements of a unit-cell cross-section and the presence of distinct pH-dependent lattice isoforms that resulted in morphological transformations. Molecular dynamics simulations demonstrated that the ionizable side chain plays a critical role in changing the local conformation of the unit cell, which further impacts the dimensionality of hierarchical structures. Moreover, these fibers were readily functionalized with biological ligands to afford one-dimensional (1D) protein arrays. This approach for the high-precision bottom-up assembly of ultranarrow 1D nanostructures offers significant potential for developing novel biomimetic nanostructures.

Lee, Yen Jea↗
DOE OSTI2025

Building molecular model series from heterogeneous CryoEM structures using Gaussian mixture models and deep neural networks

Cryogenic electron microscopy (CryoEM) produces structures of macromolecules at near-atomic resolution. However, building molecular models with good stereochemical geometry from those structures can be challenging and time-consuming, especially when many structures are obtained from datasets with conformational heterogeneity. Here we present a model refinement protocol that automatically generates series of molecular models from CryoEM datasets, which describe the dynamics of the macromolecular system and have near-perfect geometry scores. This method makes it easier to interpret the movement of the protein complex from heterogeneity analysis and to compare the structural dynamics observed from CryoEM data with results from other experimental and simulation techniques.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Challenges of conventional iterative all-atom and coarse-grained multiscale molecular dynamics

In this work, we evaluate the biomolecular dynamics behaviors when conventionally iterating between all-atom (AA) and coarse-grained (CG) molecular dynamics (MD) simulations over multiple cycles. We implemented the workflow to iterate between AA and CG in OpenMM, namely the iterative multiscale MD (iMMD) simulation workflow. In particular, we aim to identify practical applications for iterating between AA and CG simulations in a conventional manner without any constraints or model modifications. We evaluate the iMMD workflow on four representative systems, spanning folding of two soluble proteins and protein-protein as well as protein-lipid interactions of two membrane proteins. We observe that iteration between AA and CG representations could help the soluble proteins exit undesirable metastable states to fold, resulting from random protein structural distortions due to cycling. Consequently, the most reliable use of iterative AA and CG simulations appears to be to accelerating complex lipid mixing for membrane-bound protein systems rather than sampling protein conformational space. Our work explores the practical usages and limitations for iterative AA and CG simulations using readily available AA and CG force fields. The evaluated iMMD workflow in OpenMM is made available at https://github.com/lanl/iMMD.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Computational design of closely related proteins that adopt two well-defined but structurally divergent folds

The plasticity of naturally occurring protein structures, which can change shape considerably in response to changes in environmental conditions, is critical to biological function. While computational methods have been used for de novo design of proteins that fold to a single state with a deep free-energy minimum [P.-S. Huang, S. E. Boyken, D. Baker, Nature 537, 320–327 (2016)], and to reengineer natural proteins to alter their dynamics [J. A. Davey, A. M. Damry, N. K. Goto, R. A. Chica, Nat. Chem. Biol.13, 1280–1285 (2017)] or fold [P. A. Alexander, Y. He, Y. Chen, J. Orban, P. N. Bryan, Proc. Natl. Acad. Sci. U.S.A. 106, 21149–21154 (2009)], the de novo design of closely related sequences which adopt well-defined but structurally divergent structures remains an outstanding challenge. We designed closely related sequences (over 94% identity) that can adopt two very different homotrimeric helical bundle conformations—one short (~66 Å height) and the other long (~100 Å height)—reminiscent of the conformational transition of viral fusion proteins. Crystallographic and NMR spectroscopic characterization shows that both the short- and long-state sequences fold as designed. We sought to design bistable sequences for which both states are accessible, and obtained a single designed protein sequence that populates either the short state or the long state depending on the measurement conditions. The design of sequences which are poised to adopt two very different conformations sets the stage for creating large-scale conformational switches between structurally divergent forms.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Dual arginine recognition of LRRK2 phosphorylated Rab GTPases

Parkinson’s-disease-associated LRRK2 is a multidomain Ser/Thr kinase that phosphorylates a subset of Rab GTPases to control their effector functions. Rab GTPases are the prime regulators of membrane trafficking in eukaryotic cells. Rabs exert their biological effects by recruitment of effector proteins to subcellular compartments via their Rab-binding domain (RBD). Effectors are modular and typically contain additional domains that regulate various aspects of vesicle formation, trafficking, fusion, and organelle dynamics. The RBD of effectors is typically an α-helical coiled coil that recognizes the GTP conformation of the switch 1 and switch 2 motifs of Rabs. LRRK2 phosphorylates Rab8a at T72 (pT72) of its switch 2 α-helix. This post-translational modification enables recruitment of RILPL2, an effector that regulates ciliogenesis in model cell lines. A newly identified RBD motif of RILPL2, termed the X-cap, has been shown to recognize the phosphate via direct interactions between an arginine residue (R132) and pT72 of Rab8a. Here, we show that a second distal arginine (R130) is also essential for phospho-Rab binding by RILPL2. Through structural, biophysical, and cellular studies, we find that R130 stabilizes the primary R132:pT72 salt bridge through favorable enthalpic contributions to the binding affinity. These findings may have implications for the mechanism by which LRRK2 activation leads to assembly of phospho-Rab complexes and subsequent control of their membrane trafficking functions in cells.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Mechanisms of activation and desensitization of full-length glycine receptor in lipid nanodiscs

Glycinergic synapses play a central role in motor control and pain processing in the central nervous system. Glycine receptors (GlyRs) are key players in mediating fast inhibitory neurotransmission at these synapses. While previous high-resolution structures have provided insights into the molecular architecture of GlyR, several mechanistic questions pertaining to channel function are still unanswered. Here, we present Cryo-EM structures of the full-length GlyR protein complex reconstituted into lipid nanodiscs that are captured in the unliganded (closed), glycine-bound (open and desensitized), and allosteric modulator-bound conformations. A comparison of these states reveals global conformational changes underlying GlyR channel gating and modulation. The functional state assignments were validated by molecular dynamics simulations, and the observed permeation events are in agreement with the anion selectivity and conductance of GlyR. These studies provide the structural basis for gating, ion selectivity, and single-channel conductance properties of GlyR in a lipid environment.

59 BASIC BIOLOGICAL SCIENCES↗