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 181 records · Page 10

Acoustic force spectroscopy reveals subtle differences in cellulose unbinding behavior of carbohydrate-binding modules

Protein adsorption to solid carbohydrate interfaces is critical to many biological processes, particularly in biomass deconstruction. To engineer more-efficient enzymes for biomass deconstruction into sugars, it is necessary to characterize the complex protein–carbohydrate interfacial interactions. A carbohydrate-binding module (CBM) is often associated with microbial surface-tethered cellulosomes or secreted cellulase enzymes to enhance substrate accessibility. However, it is not well known how CBMs recognize, bind, and dissociate from polysaccharides to facilitate efficient cellulolytic activity, due to the lack of mechanistic understanding and a suitable toolkit to study CBM–substrate interactions. Our work outlines a general approach to study the unbinding behavior of CBMs from polysaccharide surfaces using a highly multiplexed single-molecule force spectroscopy assay. Here, we apply acoustic force spectroscopy (AFS) to probe a Clostridium thermocellum cellulosomal scaffoldin protein (CBM3a) and measure its dissociation from nanocellulose surfaces at physiologically relevant, low force loading rates. An automated microfluidic setup and method for uniform deposition of insoluble polysaccharides on the AFS chip surfaces are demonstrated. The rupture forces of wild-type CBM3a, and its Y67A mutant, unbinding from nanocellulose surfaces suggests distinct multimodal CBM binding conformations, with structural mechanisms further explored using molecular dynamics simulations. Applying classical dynamic force spectroscopy theory, the single-molecule unbinding rate at zero force is extrapolated and found to agree with bulk equilibrium unbinding rates estimated independently using quartz crystal microbalance with dissipation monitoring. However, our results also highlight critical limitations of applying classical theory to explain the highly multivalent binding interactions for cellulose–CBM bond rupture forces exceeding 15 pN.

09 BIOMASS FUELS↗

ezAlign: A Tool for Converting Coarse-Grained Molecular Dynamics Structures to Atomistic Resolution for Multiscale Modeling

Soft condensed matter is challenging to study due to the vast time and length scales that are necessary to accurately represent complex systems and capture their underlying physics. Multiscale simulations are necessary to study processes that have disparate time and/or length scales, which abound throughout biology and other complex systems. Herein we present ezAlign, an open-source software for converting coarse-grained molecular dynamics structures to atomistic representation, allowing multiscale modeling of biomolecular systems. The ezAlign v1.1 software package is publicly available for download at github.com/LLNL/ezAlign. Its underlying methodology is based on a simple alignment of an atomistic template molecule, followed by position-restraint energy minimization, which forces the atomistic molecule to adopt a conformation consistent with the coarse-grained molecule. The molecules are then combined, solvated, minimized, and equilibrated with position restraints. Validation of the process was conducted on a pure POPC membrane and compared with other popular methods to construct atomistic membranes. Additional examples, including surfactant self-assembly, membrane proteins, and more complex bacterial and human plasma membrane models, are also presented. By providing these examples, parameter files, code, and an easy-to-follow recipe to add new molecules, this work will aid future multiscale modeling efforts.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Molecular mechanism of biased signaling at the kappa opioid receptor

The κ-opioid receptor (KOR) has emerged as an attractive drug target for pain management without addiction, and biased signaling through particular pathways of KOR may be key to maintaining this benefit while minimizing side-effect liabilities. As for most G protein-coupled receptors (GPCRs), however, the molecular mechanisms of ligand-specific signaling at KOR have remained unclear. To better understand the molecular determinants of KOR signaling bias, we apply structure determination, atomic-level molecular dynamics (MD) simulations, and functional assays. We determine a crystal structure of KOR bound to the G protein-biased agonist nalfurafine, the first approved KOR-targeting drug. We also identify an arrestin-biased KOR agonist, WMS-X600. Using MD simulations of KOR bound to nalfurafine, WMS-X600, and a balanced agonist U50,488, we identify three active-state receptor conformations, including one that appears to favor arrestin signaling over G protein signaling and another that appears to favor G protein signaling over arrestin signaling. These results, combined with mutagenesis validation, provide a molecular explanation of how agonists achieve biased signaling at KOR.

59 BASIC BIOLOGICAL SCIENCES↗

Assembly of JAZ–JAZ and JAZ–NINJA complexes in jasmonate signaling

Jasmonates (JAs) are plant hormones with crucial roles in development and stress resilience. They activate MYC transcription factors by mediating the proteolysis of MYC inhibitors called JAZ proteins. In the absence of JA, JAZ proteins bind and inhibit MYC through the assembly of MYC–JAZ–Novel Interactor of JAZ (NINJA)–TPL repressor complexes. However, JAZ and NINJA are predicted to be largely intrinsically unstructured, which has precluded their experimental structure determination. Through a combination of biochemical, mutational, and biophysical analyses and AlphaFold-derived ColabFold modeling, we characterized JAZ–JAZ and JAZ–NINJA interactions and generated models with detailed, high-confidence domain interfaces. We demonstrate that JAZ, NINJA, and MYC interface domains are dynamic in isolation and become stabilized in a stepwise order upon complex assembly. By contrast, most JAZ and NINJA regions outside of the interfaces remain highly dynamic and cannot be modeled in a single conformation. Our data indicate that the small JAZ Zinc finger expressed in Inflorescence Meristem (ZIM) motif mediates JAZ–JAZ and JAZ–NINJA interactions through separate surfaces, and our data further suggest that NINJA modulates JAZ dimerization. This study advances our understanding of JA signaling by providing insights into the dynamics, interactions, and structure of the JAZ–NINJA core of the JA repressor complex.

59 BASIC BIOLOGICAL SCIENCES↗

Binding Interactions between Receptor-Binding Domain of Spike Protein and Human Angiotensin Converting Enzyme-2 in Omicron Variant

The emergence of new SARS-CoV-2 Omicron variant of concern (OV) has exacerbated the COVID-19 pandemic because of a large number of mutations in the spike protein, particularly in the receptor-binding domain (RBD), resulting in highly contagious and/or vaccine-resistant strains. Herein, we present a systematic analysis based on detailed molecular dynamics (MD) simulations in order to understand how the OV RBD mutations affect the ACE2 binding. We show that the OV RBD binds to ACE2 more efficiently and tightly predominantly because of strong electrostatic interactions, thereby promoting increased infectivity and transmissibility compared to other strains. Some of the OV RBD mutations are predicted to affect the antibody neutralization either through their role in the S-protein conformational changes, such as S371L, S373P, and S375F, or through changing its surface charge distribution, such as G339D, N440K, T478K, and E484A. Other mutations, such as K417N, G446S, and Y505H, decrease the ACE2 binding, whereas S447N, Q493R, G496S, Q498R, and N501Y tend to increase it.

60 APPLIED LIFE SCIENCES↗

HIV Env assembling is tightly modulated by lipid composition and cytoplasmic tail integrity

HIV envelope glycoprotein (Env) is a homotrimeric transmembrane protein which binds to receptors CD4, CCR5- CXCR4, leading to a cascade of conformational changes ultimately resulting in the entrance of the virus into the host cell. Its whole structure comprises an extracellular region, a transmembrane domain and a cytoplasmic region. While many efforts are carried in order to develop an effective vaccine specifically targeting the extracellular domain, less efforts are focus into regions close to the transmembrane region; mostly limited by the lack of structural knowledge within this domain. Thus, it is desirable to have a molecular description of the trimeric association of the transmembrane domain, as well as the membrane properties affecting their structure. We have carried molecular dynamic (MD) simulations of HIV transmembrane region in different lipid environments, mimicking the complex viral membrane. Our results indicate that the lipid environment heavily affects the dynamics and association of the transmembrane helices, modulating its structure. In addition, we have investigated the effect of the cytoplasmic tail (CT). Unexpectedly, this region has a tremendous effect on the lateral organization of the transmembrane. These results clearly highlights the importance of this region in the overall structure of the Env protein and point towards its key role for the in-silico design of effective vaccines towards this domain.

59 BASIC BIOLOGICAL SCIENCES↗

Revealing a Hidden Intermediate of Rotatory Catalysis with X-ray Crystallography and Molecular Simulations

The mechanism of rotatory catalysis in ATPhydrolyzing molecular motors remains an unresolved puzzle in biological energy transfer. Notwithstanding the wealth of available biochemical and structural information inferred from years of experiments, knowledge on how the coupling between the chemical and mechanical steps within motors enforces directional rotatory movements remains fragmentary. Even more contentious is to pinpoint the rate-limiting step of a multistep rotation process. Here, using vacuolar or V 1 -type hexameric ATPase as an exemplary rotational motor, we present a model of the complete 4-step conformational cycle involved in rotatory catalysis. First, using Xray crystallography, a new intermediate or “dwell” is identified, which enables the release of an inorganic phosphate (or P i ) after ATP hydrolysis. Using molecular dynamics simulations, this new dwell is placed in a sequence with three other crystal structures to derive a putative cyclic rotation path. Free-energy simulations are employed to estimate the rate of the hexameric protein transformations and delineate allosteric effects that allow new reactant ATP entry only after hydrolysis product exit. An analysis of transfer entropy brings to light how the side-chain-level interactions transcend into larger-scale reorganizations, highlighting the role of the ubiquitous arginine-finger residues in coupling chemical and mechanical information. An inspection of all known rates encompassing the 4-step rotation mechanism implicates the overcoming of the ADP interactions with V 1 -ATPase to be the rate-limiting step of motor action.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A simple technique to classify diffraction data from dynamic proteins according to individual polymorphs

One often observes small but measurable differences in the diffraction data measured from different crystals of a single protein. These differences might reflect structural differences in the protein and may reveal the natural dynamism of the molecule in solution. Partitioning these mixed-state data into single-state clusters is a critical step that could extract information about the dynamic behavior of proteins from hundreds or thousands of single-crystal data sets. Mixed-state data can be obtained deliberately (through intentional perturbation) or inadvertently (while attempting to measure highly redundant single-crystal data). To the extent that different states adopt different molecular structures, one expects to observe differences in the crystals; each of the polystates will create a polymorph of the crystals. After mixed-state diffraction data have been measured, deliberately or inadvertently, the challenge is to sort the data into clusters that may represent relevant biological polystates. Here, this problem is addressed using a simple multi-factor clustering approach that classifies each data set using independent observables, thereby assigning each data set to the correct location in conformational space. This procedure is illustrated using two independent observables, unit-cell parameters and intensities, to cluster mixed-state data from chymotrypsinogen (ChTg) crystals. It is observed that the data populate an arc of the reaction trajectory as ChTg is converted into chymotrypsin.

36 MATERIALS SCIENCE↗

Sequence Programmable Order–Disorder Transitions in Supramolecular Assembly of Peptide Nanofibers

Protein–protein interactions determine the assembly of complexes that are responsible for numerous key biological processes. The assembly of many natural protein complexes is mediated by post-translational structural changes and environmental stimuli. In this study, we show that incorporation of adjacent lysine residues results in the pH-tunable stability of peptide secondary structure and assembly, allowing for the incorporation of complementary order-inducing motifs. The strategic placement of cysteine pairs in the same peptide sequence results in redox-dependent disulfide staple formation, inducing a transition from random coil to β-sheet conformation and subsequent supramolecular nanofiber assembly from otherwise disordered peptide monomers. Spectroscopic, imaging, molecular dynamics, and kinetic studies highlight the critical role of sequence motif location, oligomerization, and the competitive interplay between intra- and interpeptide disulfide bonding in determining assembly outcomes. We extend this approach to demonstrate phosphorylation-dependent assembly from the design of the same parent peptide sequence, suggesting a general approach to the design of diverse stimulus-responsive peptide sequences for supramolecular assembly. Furthermore, these findings also provide a framework for investigating sequence-dependent pathways in amyloid fiber formation with potential implications for neurodegenerative disease research.

Disulfides↗

Unravelling the dynamics of the maturation protein in MS2 bacteriophage via molecular simulations

The MS2 bacteriophage capsid serves as a model system for studying viral structure and function. Mature MS2 virus consists of 178 capsid proteins and a single maturation protein (MP), which is essential for host receptor binding and infection initiation. Despite its critical role, the dynamic behavior of the capsid with the MP remains poorly understood. To address this, we conducted 0.5 µs all-atom molecular dynamics (MD) simulations of the MS2 capsid with and without the MP, revealing key insights into its structural dynamics. Our simulations showed that MP exhibits high flexibility, particularly in the “tip” and “side-loop” regions, which undergo significant motions that likely enhance its ability to engage with the F-pilus receptor. Detailed analysis of MP conformational states revealed that loop rearrangements around H357 enable transient switching between “semi-closed” and “open” conformations, suggesting a conformational selection mechanism for pilus binding. Additionally, ion interaction analyses revealed distinct sodium and chloride binding patterns, where sodium ions were mostly found at the outer capsid shell, while chloride ions interacted with the basic residues on the RNA-facing side. We also found that the presence of the MP enhances salt-bridge interactions, contributing to increased capsid stability, yet it does not significantly alter the pore sizes of pentameric and hexameric units. Together, these findings provide new insights into the functional role of the MP, highlighting its contribution to capsid stability and host receptor engagement. This study offers a foundation for understanding capsid dynamics relevant to viral infectivity and may guide future rational strategies aimed at disrupting host-virus interactions.

Capsid stability↗

Photoactivation Transition State and Dynamical Response of the Orange Carotenoid Protein

The orange carotenoid protein (OCP) regulates light harvesting in cyanobacteria by acting as a photoreceptor in its resting form, OCP O , and by effecting the quenching of bilin excitons upon binding to the core of the phycobilisome in its photoactivated red form, OCP R . We show herein using fluorescence anisotropy measurements and the action spectra for the rate constants of the two light-driven steps in the mechanism that the photoactivation of the OCP from Synechocystis sp. PCC 6803 is triggered by excited-state motions of the canthaxanthin chromophore that yield a twisted and bent conformation. Well-tempered metadynamics simulations reveal that a bicycle-pedal configuration, due to twisting of the two adjacent C=C bonds at the C13−C14 and C15−C15′ positions in the center of canthaxanthin’s π-conjugated isoprenoid backbone, can be accommodated by the binding site in the OCP, with the energy of a captured photon required to cross the local activation energy barriers from the dark equilibrium structure. The bicycle-pedal configuration breaks the conserved hydrogen-bonding interactions between the carbonyl substituent of the β-ionone end ring of canthaxanthin and the adjacent W288 and Y201 residues in the C-terminal domain. The action spectra are modulated by the vibronic excitation prepared by absorption transitions to the S 2 state, indicating that the photoactivation reactions are triggered by the canthaxanthin chromophore well prior to vibrational equilibration. These findings show that an ultrafast structural response of the OCP protein to the excited-state motions of the canthaxanthin chromophore controls the photoactivation yield and the sensing of blue light.

Fluorescence↗

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.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The structural plasticity of nucleic acid duplexes revealed by WAXS and MD

Double-stranded DNA (dsDNA) and RNA (dsRNA) helices display an unusual structural diversity. Some structural variations are linked to sequence and may serve as signaling units for protein-binding partners. Therefore, elucidating the mechanisms and factors that modulate these variations is of fundamental importance. While the structural diversity of dsDNA has been extensively studied, similar studies have not been performed for dsRNA. Because of the increasing awareness of RNA’s diverse biological roles, such studies are timely and increasingly important. We integrate solution x-ray scattering at wide angles (WAXS) with all-atom molecular dynamics simulations to explore the conformational ensemble of duplex topologies for different sequences and salt conditions. These tightly coordinated studies identify robust correlations between features in the WAXS profiles and duplex geometry and enable atomic-level insights into the structural diversity of DNA and RNA duplexes. Notably, dsRNA displays a marked sensitivity to the valence and identity of its associated cations.

99 GENERAL AND MISCELLANEOUS↗

Structure and functional determinants of Rad6–Bre1 subunits in the histone H2B ubiquitin-conjugating complex

The conserved complex of the Rad6 E2 ubiquitin-conjugating enzyme and the Bre1 E3 ubiquitin ligase catalyzes histone H2B monoubiquitination (H2Bub1), which regulates chromatin dynamics during transcription and other nuclear processes. Here, we report a crystal structure of Rad6 and the non-RING domain N-terminal region of Bre1, which shows an asymmetric homodimer of Bre1 contacting a conserved loop on the Rad6 ‘backside’. This contact is distant from the Rad6 catalytic site and is the location of mutations that impair telomeric silencing in yeast. Mutational analyses validated the importance of this contact for the Rad6–Bre1 interaction, chromatin-binding dynamics, H2Bub1 formation and gene expression. Moreover, the non-RING N-terminal region of Bre1 is sufficient to confer nucleosome binding ability to Rad6 in vitro. Interestingly, Rad6 P43L protein, an interaction interface mutant and equivalent to a cancer mutation in the human homolog, bound Bre1 5-fold more tightly than native Rad6 in vitro, but showed reduced chromatin association of Bre1 and reduced levels of H2Bub1 in vivo. These surprising observations imply conformational transitions of the Rad6–Bre1 complex during its chromatin-associated functional cycle, and reveal the differential effects of specific disease-relevant mutations on the chromatin-bound and unbound states. Overall, our study provides structural insights into Rad6–Bre1 interaction through a novel interface that is important for their biochemical and biological responses.

59 BASIC BIOLOGICAL SCIENCES↗

Generative $β$-hairpin design using a residue-based physicochemical property landscape

De novo peptide design is a new frontier that has broad application potential in the biological and biomedical fields. Most existing models for de novo peptide design are largely based on sequence homology that can be restricted based on evolutionarily derived protein sequences and lack the physicochemical context essential in protein folding. Generative machine learning for de novo peptide design is a promising way to synthesize theoretical data that are based on, but unique from, the observable universe. In this study, we created and tested a custom peptide generative adversarial network intended to design peptide sequences that can fold into the -hairpin secondary structure. This deep neural network model is designed to establish a preliminary foundation of the generative approach based on physicochemical and conformational properties of 20 canonical amino acids, for example, hydrophobicity and residue volume, using extant structure-specific sequence data from the PDB. The beta generative adversarial network model robustly distinguishes secondary structures of hairpin from α helix and intrinsically disordered peptides with an accuracy of up to 96% and generates artificial -hairpin peptide sequences with minimum sequence identities around 31% and 50% when compared against the current NCBI PDB and nonredundant databases, respectively. These results highlight the potential of generative models specifically anchored by physicochemical and conformational property features of amino acids to expand the sequence-to-structure landscape of proteins beyond evolutionary limits.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

MHC-II dynamics are maintained in HLA-DR allotypes to ensure catalyzed peptide exchange

Presentation of antigenic peptides by major histocompatibility complex class II (MHC-II) proteins determines T helper cell reactivity. The MHC-II genetic locus displays a large degree of allelic polymorphism influencing the peptide repertoire presented by the resulting MHC-II protein allotypes. During antigen processing, the human leukocyte antigen (HLA) molecule HLA-DM (DM) encounters these distinct allotypes and catalyzes exchange of the placeholder peptide CLIP by exploiting dynamic features of MHC-II. Here, we investigate 12 highly abundant CLIP-bound HLA-DRB1 allotypes and correlate dynamics to catalysis by DM. Despite large differences in thermodynamic stability, peptide exchange rates fall into a target range that maintains DM responsiveness. A DM-susceptible conformation is conserved in MHC-II molecules, and allosteric coupling between polymorphic sites affects dynamic states that influence DM catalysis. As exemplified for rheumatoid arthritis, we postulate that intrinsic dynamic features of peptide–MHC-II complexes contribute to the association of individual MHC-II allotypes with autoimmune disease.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗