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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Enhanced Binding of SARS-CoV-2 Spike Protein to Receptor by Distal Polybasic Cleavage Sites

The receptor-binding domain (RBD) of the SARS-CoV-2 spike protein plays a crucial role in binding the human cell receptor ACE2 that is required for viral entry. Many studies have been conducted to target the structures of RBD–ACE2 binding and to design RBD-targeting vaccines and drugs. Nevertheless, mutations distal from the SARS-CoV-2 RBD also impact its transmissibility and antibody can target non-RBD regions, suggesting the incomplete role of the RBD region in the spike protein–ACE2 binding. Here, in order to elucidate distant binding mechanisms, we analyze complexes of ACE2 with the wild-type spike protein and with key mutants via large-scale all-atom explicit solvent molecular dynamics simulations. We find that though distributed approximately 10 nm away from the RBD, the SARS-CoV-2 polybasic cleavage sites enhance, via electrostatic interactions and hydration, the RBD–ACE2 binding affinity. A negatively charged tetrapeptide (GluGluLeuGlu) is then designed to neutralize the positively charged arginine on the polybasic cleavage sites. We find that the tetrapeptide GluGluLeuGlu binds to one of the three polybasic cleavage sites of the SARS-CoV-2 spike protein lessening by 34% the RBD–ACE2 binding strength. This significant binding energy reduction demonstrates the feasibility to neutralize RBD–ACE2 binding by targeting this specific polybasic cleavage site. Our work enhances understanding of the binding mechanism of SARS-CoV-2 to ACE2, which may aid the design of therapeutics for COVID-19 infection.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1997

Mechanical continuity and reversible chromosome disassembly within intact genomes removed from living cells

Chromatin is thought to be structurally discontinuous because it is packaged into morphologically distinct chromosomes that appear physically isolated from one another in metaphase preparations used for cytogenetic studies. However, analysis of chromosome positioning and movement suggest that different chromosomes often behave as if they were physically connected in interphase as well as mitosis. To address this paradox directly, we used a microsurgical technique to physically remove nucleoplasm or chromosomes from living cells under isotonic conditions. Using this approach, we found that pulling a single nucleolus or chromosome out from interphase or mitotic cells resulted in sequential removal of the remaining nucleoli and chromosomes, interconnected by a continuous elastic thread. Enzymatic treatments of interphase nucleoplasm and chromosome chains held under tension revealed that mechanical continuity within the chromatin was mediated by elements sensitive to DNase or micrococcal nuclease, but not RNases, formamide at high temperature, or proteases. In contrast, mechanical coupling between mitotic chromosomes and the surrounding cytoplasm appeared to be mediated by gelsolin-sensitive microfilaments. Furthermore, when ion concentrations were raised and lowered, both the chromosomes and the interconnecting strands underwent multiple rounds of decondensation and recondensation. As a result of these dynamic structural alterations, the mitotic chains also became sensitive to disruption by restriction enzymes. Ion-induced chromosome decondensation could be blocked by treatment with DNA binding dyes, agents that reduce protein disulfide linkages within nuclear matrix, or an antibody directed against histones. Fully decondensed chromatin strands also could be induced to recondense into chromosomes with pre-existing size, shape, number, and position by adding anti-histone antibodies. Conversely, removal of histones by proteolysis or heparin treatment produced chromosome decondensation which could be reversed by addition of histone H1, but not histones H2b or H3. These data suggest that DNA, its associated protein scaffolds, and surrounding cytoskeletal networks function as a structurally-unified system. Mechanical coupling within the nucleoplasm may coordinate dynamic alterations in chromatin structure, guide chromosome movement, and ensure fidelity of mitosis.

NASA Discipline Cell Biology↗
DOE OSTI2020

Assessment of software methods for estimating protein-protein relative binding affinities

A growing number of computational tools have been developed to accurately and rapidly predict the impact of amino acid mutations on protein-protein relative binding affinities. Such tools have many applications, for example, designing new drugs and studying evolutionary mechanisms. In the search for accuracy, many of these methods employ expensive yet rigorous molecular dynamics simulations. By contrast, non-rigorous methods use less exhaustive statistical mechanics, allowing for more efficient calculations. However, it is unclear if such methods retain enough accuracy to replace rigorous methods in binding affinity calculations. This trade-off between accuracy and computational expense makes it difficult to determine the best method for a particular system or study. Here, eight non-rigorous computational methods were assessed using eight antibody-antigen and eight non-antibody-antigen complexes for their ability to accurately predict relative binding affinities (ΔΔG) for 654 single mutations. In addition to assessing accuracy, we analyzed the CPU cost and performance for each method using a variety of physico-chemical structural features. This allowed us to posit scenarios in which each method may be best utilized. Most methods performed worse when applied to antibody-antigen complexes compared to non-antibody-antigen complexes. Rosetta-based JayZ and EasyE methods classified mutations as destabilizing (ΔΔG < -0.5 kcal/mol) with high (83–98%) accuracy and a relatively low computational cost for non-antibody-antigen complexes. Some of the most accurate results for antibody-antigen systems came from combining molecular dynamics with FoldX with a correlation coefficient (r) of 0.46, but this was also the most computationally expensive method. Overall, our results suggest these methods can be used to quickly and accurately predict stabilizing versus destabilizing mutations but are less accurate at predicting actual binding affinities. This study highlights the need for continued development of reliable, accessible, and reproducible methods for predicting binding affinities in antibody-antigen proteins and provides a recipe for using current methods.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Covalent Drug Binding in Live Cells Monitored by Mid-Infrared Quantum Cascade Laser Spectroscopy: Photoactive Yellow Protein as a Model System

The detection of drug-target interactions in live cells enables analysis of therapeutic compounds in a native cellular environment. Recent advances in spectroscopy and molecular biology have facilitated the development of genetically encoded vibrational probes like nitriles that can sensitively report on molecular interactions. Nitriles are powerful tools for measuring electrostatic environments within condensed media like proteins, but such measurements in live cells have been hindered by low signal-to-noise ratios. In this study, we design a spectrometer based on a double-beam quantum cascade laser (QCL)-based transmission infrared (IR) source with balanced detection that can significantly enhance sensitivity to nitrile vibrational probes embedded in proteins within cells compared to a conventional FTIR spectrometer. Here, using this approach, we detect small-molecule binding in Escherichia coli, with particular focus on the interaction between para-Coumaric acid (pCA) and nitrile-incorporated photoactive yellow protein (PYP). This system effectively serves as a model for investigating covalent drug binding in a cellular environment. Notably, we observe large spectral shifts of up to 15 cm –1 for nitriles embedded in PYP between the unbound and drug-bound states directly within bacteria, in agreement with observations for purified proteins. Such large spectral shifts are ascribed to the changes in the hydrogen-bonding environment around the local environment of nitriles, accurately modeled through high-level molecular dynamics simulations using the AMOEBA force field. Our findings underscore the QCL spectrometer’s ability to enhance sensitivity for monitoring drug–protein interactions, offering new opportunities for advanced methodologies in drug development and biochemical research.

chromophores↗
DOE OSTI2021

Molecular Dynamics Free Energy Simulations Reveal the Mechanism for the Antiviral Resistance of the M66I HIV-1 Capsid Mutation

While drug resistance mutations can often be attributed to the loss of direct or solvent-mediated protein−ligand interactions in the drug-mutant complex, in this study we show that a resistance mutation for the picomolar HIV-1 capsid (CA)-targeting antiviral (GS-6207) is mainly due to the free energy cost of the drug-induced protein side chain reorganization in the mutant protein. Among several mutations, M66I causes the most suppression of the GS-6207 antiviral activity (up to ~84,000-fold), and only 83- and 68-fold reductions for PF74 and ZW-1261, respectively. To understand the molecular basis of this drug resistance, we conducted molecular dynamics free energy simulations to study the structures, energetics, and conformational free energy landscapes involved in the inhibitors binding at the interface of two CA monomers. To minimize the protein−ligand steric clash, the I66 side chain in the M66I−GS-6207 complex switches to a higher free energy conformation from the one adopted in the apo M66I. In contrast, the binding of GS-6207 to the wild-type CA does not lead to any significant M66 conformational change. Based on an analysis that decomposes the absolute binding free energy into contributions from two receptor conformational states, it appears that it is the free energy cost of side chain reorganization rather than the reduced protein−ligand interaction that is largely responsible for the drug resistance against GS-6207.

Sun, Qinfang (ORCID:0000000173621543)↗
DOE OSTI2024

Theoretical examination of nuclear spin diffusion in light-induced spin coherences in photosystem I

Microsecond coherence times are predicted in spin dynamics calculations of a donor–acceptor electron spin pair P + A 1A - created upon light excitation of photosystem I (PSI). The effects of nuclear spin diffusion (NSD) due to the anisotropic protein environment on the predicted coherence times T M are studied. Closely positioned pairs and triples of protons located 5–8 Å from the electron spin are shown to largely control T M . Additionally, knowledge of the crystal structure of PSI allows for spin dynamics calculations where specific cofactors and amino acid residues are removed or replaced and the identification of anisotropic environmental features controlling electronic decoherence. Finally, we show that NSD alone cannot explain the >3× shorter experimentally observed coherence times and suggest that methyl groups at key protein sites may explain this discrepancy.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1993

Device For Controlling Crystallization Of Protein

Variable sandwich spacer enables optimization of evaporative driving force that governs crystallization of protein from solution. Mechanically more rigid than hanging-drop and sitting-drop devices. Large oscillations and dislodgment of drop of solution in response to vibrations suppressed by glass plates. Other advantages include: suitable for automated delivery, stable handling, and programmable evaporation of protein solution; controlled configuration enables simple and accurate determination of volume of solution without disrupting crystallization; pH and concentration of precipitant controlled dynamically because pH and concentration coupled to rate of evaporation, controllable via adjustment of gap between plates; and enables variation of ratio between surface area and volume of protein solution. Alternative version, plates oriented vertically instead of horizontally.

Noever, David A.↗
DOE OSTI2021

Structures of the mycobacterial membrane protein MmpL3 reveal its mechanism of lipid transport

The mycobacterial membrane protein large 3 (MmpL3) transporter is essential and required for shuttling the lipid trehalose monomycolate (TMM), a precursor of mycolic acid (MA)-containing trehalose dimycolate (TDM) and mycolyl arabinogalactan peptidoglycan (mAGP), in Mycobacterium species, including Mycobacterium tuberculosis and Mycobacterium smegmatis. However, the mechanism that MmpL3 uses to facilitate the transport of fatty acids and lipidic elements to the mycobacterial cell wall remains elusive. Here, we report 7 structures of the M. smegmatis MmpL3 transporter in its unbound state and in complex with trehalose 6-decanoate (T6D) or TMM using single-particle cryo-electron microscopy (cryo-EM) and X-ray crystallography. Combined with calculated results from molecular dynamics (MD) and target MD simulations, we reveal a lipid transport mechanism that involves a coupled movement of the periplasmic domain and transmembrane helices of the MmpL3 transporter that facilitates the shuttling of lipids to the mycobacterial cell wall.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1992

Large scale crystallization of protein pharmaceuticals in microgravity via temperature change

The major objective of this research effort is the temperature driven growth of protein crystals in large batches in the microgravity environment of space. Pharmaceutical houses are developing protein products for patient care, for example, human insulin, human growth hormone, interferons, and tissue plasminogen activator or TPA, the clot buster for heart attack victims. Except for insulin, these are very high value products; they are extremely potent in small quantities and have a great value per gram of material. It is feasible that microgravity crystallization can be a cost recoverable, economically sound final processing step in their manufacture. Large scale protein crystal growth in microgravity has significant advantages from the basic science and the applied science standpoints. Crystal growth can proceed unhindered due to lack of surface effects. Dynamic control is possible and relatively easy. The method has the potential to yield large quantities of pure crystalline product. Crystallization is a time honored procedure for purifying organic materials and microgravity crystallization could be the final step to remove trace impurities from high value protein pharmaceuticals. In addition, microgravity grown crystals could be the final formulation for those medicines that need to be administered in a timed release fashion. Long lasting insulin, insulin lente, is such a product. Also crystalline protein pharmaceuticals are more stable for long-term storage. Temperature, as the initiation step, has certain advantages. Again, dynamic control of the crystallization process is possible and easy. A temperature step is non-invasive and is the most subtle way to control protein solubility and therefore crystallization. Seeding is not necessary. Changes in protein and precipitant concentrations and pH are not necessary. Finally, this method represents a new way to crystallize proteins in space that takes advantage of the unique microgravity environment. The results from two flights showed that the hardware performed perfectly, many crystals were produced, and they were much larger than their ground grown controls. Morphometric analysis was done on over 4,000 crystals to establish crystal size, size distribution, and relative size. Space grown crystals were remarkably larger than their earth grown counterparts and crystal size was a function of PCF volume. That size distribution for the space grown crystals was a function of PCF volume may indicate that ultimate size was a function of temperature gradient. Since the insulin protein concentration was very low, 0.4 mg/ml, the size distribution could also be following the total amount of protein in each of the PCF's. X-ray analysis showed that the bigger space grown insulin crystals diffracted to higher resolution than their ground grown controls. When the data were normalized for size, they still indicated that the space crystals were better than the ground crystals.

Long, Marianna M.↗
DOE OSTI2020

In Silico Insights into Protein–Protein Interaction Disruptive Mutations in the PCSK9-LDLR Complex

Gain-of-function mutations in PCSK9 (proprotein convertase subtilisin/kexin type 9) lead to reduced uptake of LDL (low density lipoprotein) cholesterol and, therefore, increased plasma LDL levels. However, the mechanism by which these mutants reduce LDL reuptake is not fully understood. Here, we have used molecular dynamics simulations, MM/PBSA (Molecular Mechanics/Poisson–Boltzmann Surface Area) binding affinity calculations, and residue interaction networks, to investigate the protein–protein interaction (PPI) disruptive effects of two of PCSK9's gain-of-function mutations, Ser127Arg and Asp374Tyr on the PCSK9 and LDL receptor complex. In addition to these PPI disruptive mutants, a third, non-interface mutation (Arg496Trp) is included as a positive control. Our results indicate that Ser127Arg and Asp374Tyr confer significantly improved binding affinity, as well as different binding modes, when compared to the wild-type. These PPI disruptive mutations lie between the EGF(A) (epidermal growth factor precursor homology domain A) of the LDL receptor and the catalytic domain of PCSK9 (Asp374Tyr) and between the prodomain of PCSK9 and the β-propeller of the LDL receptor (Ser127Arg). The interactions involved in these two interfaces result in an LDL receptor that is sterically inhibited from entering its closed conformation. This could potentially implicate the prodomain as a target for small molecule inhibitors.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Entropy is an important design principle in the photosystem II supercomplex

Photosystem II (PSII) can achieve near-unity quantum efficiency of light harvesting in ideal conditions and can dissipate excess light energy as heat to prevent the formation of reactive oxygen species (ROS) under light stress. Understanding how this pigment–protein complex accomplishes these opposing goals is a topic of great interest that has so far been explored primarily through the lens of the system energetics. Despite PSII’s known flat energy landscape, a thorough consideration of the entropic effects on energy transfer in PSII is lacking. In this work, we aim to discern the free energetic design principles underlying the PSII energy transfer network. To accomplish this goal, we employ a structure-based rate matrix and compute the free energy terms in time following a specific initial excitation to discern how entropy and enthalpy drive ensemble system dynamics. We find that the interplay between the entropy and enthalpy components differ among each protein subunit, which allows each subunit to fulfill a unique role in the energy transfer network. This individuality ensures that PSII can accomplish efficient energy trapping in the reaction center (RC), effective nonphotochemical quenching (NPQ) in the periphery, and robust energy trapping in the other-monomer RC if the same-monomer RC is closed. We also show that entropy, in particular, is a dynamically tunable feature of the PSII free energy landscape accomplished through regulation of LHCII binding. These findings help rationalize natural photosynthesis and provide design principles for more efficient solar energy harvesting technologies.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

ADEPT: a domain independent sequence alignment strategy for gpu architectures

Bioinformatic workflows frequently make use of automated genome assembly and protein clustering tools. At the core of most of these tools, a significant portion of execution time is spent in determining optimal local alignment between two sequences. This task is performed with the Smith-Waterman algorithm, which is a dynamic programming based method. With the advent of modern sequencing technologies and increasing size of both genome and protein databases, a need for faster Smith-Waterman implementations has emerged. Multiple SIMD strategies for the Smith-Waterman algorithm are available for CPUs. However, with the move of HPC facilities towards accelerator based architectures, a need for an efficient GPU accelerated strategy has emerged. Existing GPU based strategies have either been optimized for a specific type of characters (Nucleotides or Amino Acids) or for only a handful of application use-cases. In this paper, we present ADEPT, a new sequence alignment strategy for GPU architectures that is domain independent, supporting alignment of sequences from both genomes and proteins. Our proposed strategy uses GPU specific optimizations that do not rely on the nature of sequence. We demonstrate the feasibility of this strategy by implementing the Smith-Waterman algorithm and comparing it to similar CPU strategies as well as the fastest known GPU methods for each domain. ADEPT’s driver enables it to scale across multiple GPUs and allows easy integration into software pipelines which utilize large scale computational systems. We have shown that the ADEPT based Smith-Waterman algorithm demonstrates a peak performance of 360 GCUPS and 497 GCUPs for protein based and DNA based datasets respectively on a single GPU node (8 GPUs) of the Cori Supercomputer. Overall ADEPT shows 10x faster performance in a node-to-node comparison against a corresponding SIMD CPU implementation. ADEPT demonstrates a performance that is either comparable or better than existing GPU strategies. We demonstrated the efficacy of ADEPT in supporting existing bionformatics software pipelines by integrating ADEPT in MetaHipMer a high-performance denovo metagenome assembler and PASTIS a high-performance protein similarity graph construction pipeline. Our results show 10% and 30% boost of performance in MetaHipMer and PASTIS respectively.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Conformational restriction shapes the inhibition of a multidrug efflux adaptor protein

Abstract Membrane efflux pumps play a major role in bacterial multidrug resistance. The tripartite multidrug efflux pump system from Escherichia coli , AcrAB-TolC, is a target for inhibition to lessen resistance development and restore antibiotic efficacy, with homologs in other ESKAPE pathogens. Here, we rationalize a mechanism of inhibition against the periplasmic adaptor protein, AcrA, using a combination of hydrogen/deuterium exchange mass spectrometry, cellular efflux assays, and molecular dynamics simulations. We define the structural dynamics of AcrA and find that an inhibitor can inflict long-range stabilisation across all four of its domains, whereas an interacting efflux substrate has minimal effect. Our results support a model where an inhibitor forms a molecular wedge within a cleft between the lipoyl and αβ barrel domains of AcrA, diminishing its conformational transmission of drug-evoked signals from AcrB to TolC. This work provides molecular insights into multidrug adaptor protein function which could be valuable for developing antimicrobial therapeutics.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Conserved unique peptide patterns (CUPP) online platform 2.0: implementation of +1000 JGI fungal genomes

Carbohydrate-processing enzymes, CAZymes, are classified into families based on sequence and three-dimensional fold. Because many CAZyme families contain members of diverse molecular function (different EC-numbers), sophisticated tools are required to further delineate these enzymes. Such delineation is provided by the peptide-based clustering method CUPP, Conserved Unique Peptide Patterns. CUPP operates synergistically with the CAZy family/subfamily categorizations to allow systematic exploration of CAZymes by defining small protein groups with shared sequence motifs. The updated CUPP library contains 21,930 of such motif groups including 3,842,628 proteins. The new implementation of the CUPP-webserver, https://cupp.info/, now includes all published fungal and algal genomes from the Joint Genome Institute (JGI), genome resources MycoCosm and PhycoCosm, dynamically subdivided into motif groups of CAZymes. This allows users to browse the JGI portals for specific predicted functions or specific protein families from genome sequences. Thus, a genome can be searched for proteins having specific characteristics. All JGI proteins have a hyperlink to a summary page which links to the predicted gene splicing including which regions have RNA support. The new CUPP implementation also includes an update of the annotation algorithm that uses only a fourth of the RAM while enabling multi-threading, providing an annotation speed below 1 ms/protein.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

ddcMD: A fully GPU-accelerated molecular dynamics program for the Martini force field

We have implemented the Martini force field within Lawrence Livermore National Laboratory’s molecular dynamics program, ddcMD. The program is extended to a heterogeneous programming model so that it can exploit graphics processing unit (GPU) accelerators. In addition to the Martini force field being ported to the GPU, the entire integration step, including thermostat, barostat, and constraint solver, is ported as well, which speeds up the simulations to 278-fold using one GPU vs one central processing unit (CPU) core. A benchmark study is performed with several test cases, comparing ddcMD and GROMACS Martini simulations. The average performance of ddcMD for a protein–lipid simulation system of 136k particles achieves 1.04 µs/day on one NVIDIA V100 GPU and aggregates 6.19 µs/day on one Summit node with six GPUs. The GPU implementation in ddcMD offloads all computations to the GPU and only requires one CPU core per simulation to manage the inputs and outputs, freeing up remaining CPU resources on the compute node for alternative tasks often required in complex simulation campaigns.

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

Disruption of Very-Long-Chain-Fatty Acid Synthesis Has an Impact on the Dynamics of Cellulose Synthase in Arabidopsis thaliana

In higher plants, cellulose is synthesized by membrane-spanning large protein complexes named cellulose synthase complexes (CSCs). In this study, the Arabidopsis PASTICCINO2 (PAS2) was identified as an interacting partner of cellulose synthases. PAS2 was previously characterized as the plant 3-hydroxy-acyl-CoA dehydratase, an ER membrane-localized dehydratase that is essential for very-long-chain-fatty acid (VLCFA) elongation. The pas2-1 mutants show defective cell elongation and reduction in cellulose content in both etiolated hypocotyls and light-grown roots. Although disruption of VLCFA synthesis by a genetic alteration had a reduction in VLCFA in both etiolated hypocotyls and light-grown roots, it had a differential effect on cellulose content in the two systems, suggesting the threshold level of VLCFA for efficient cellulose synthesis may be different in the two biological systems. pas2-1 had a reduction in both CSC delivery rate and CSC velocity at the PM in etiolated hypocotyls. Interestingly, Golgi but not post-Golgi endomembrane structures exhibited a severe defect in motility. Experiments using pharmacological perturbation of VLCFA content in etiolated hypocotyls strongly indicate a novel function of PAS2 in the regulation of CSC and Golgi motility. Through a combination of genetic, biochemical and cell biology studies, our study demonstrated that PAS2 as a multifunction protein has an important role in the regulation of cellulose biosynthesis in Arabidopsis hypocotyl.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2000

Computer simulation of ion channel gating: the M(2) channel of influenza A virus in a lipid bilayer

The transmembrane fragment of the influenza virus M(2) protein forms a homotetrameric channel that transports protons. In this paper, we use molecular dynamics simulations to help elucidate the mechanism of channel gating by four histidines that occlude the channel lumen in the closed state. We test two competing hypotheses. In the "shuttle" mechanism, the delta nitrogen atom on the extracellular side of one histidine is protonated by the incoming proton, and, subsequently, the proton on the epsilon nitrogen atom is released on the opposite side. In the "water-wire" mechanism, the gate opens because of electrostatic repulsion between four simultaneously biprotonated histidines. This allows for proton transport along the water wire that penetrates the gate. For each system, composed of the channel embedded in a hydrated phospholipid bilayer, a 1.3-ns trajectory was obtained. It is found that the states involved in the shuttle mechanism, which contain either single-protonated histidines or a mixture of single-protonated histidines plus one biprotonated residue, are stable during the simulations. Furthermore, the orientations and dynamics of water molecules near the gate are conducive to proton transfer. In contrast, the fully biprotonated state is not stable. Additional simulations show that if only two histidines are biprotonated, the channel deforms but the gate remains closed. These results support the shuttle mechanism but not the gate-opening mechanism of proton gating in M(2).

NASA Center ARC↗
DOE OSTI2020

Real-time tracking of protein unfolding with time-resolved x-ray solution scattering

The correct folding of proteins is of paramount importance for their function, and protein misfolding is believed to be the primary cause of a wide range of diseases. Protein folding has been investigated with time-averaged methods and time-resolved spectroscopy, but observing the structural dynamics of the unfolding process in real-time is challenging. Here, we demonstrate an approach to directly reveal the structural changes in the unfolding reaction. We use nano- to millisecond time-resolved x-ray solution scattering to probe the unfolding of apomyoglobin. The unfolding reaction was triggered using a temperature jump, which was induced by a nanosecond laser pulse. We demonstrate a new strategy to interpret time-resolved x-ray solution scattering data, which evaluates ensembles of structures obtained from molecular dynamics simulations. We find that apomyoglobin passes three states when unfolding, which we characterize as native, molten globule, and unfolded. The molten globule dominates the population under the conditions investigated herein, whereas native and unfolded structures primarily contribute before the laser jump and 30 μs after it, respectively. The molten globule retains much of the native structure but shows a dynamic pattern of inter-residue contacts. Our study demonstrates a new strategy to directly observe structural changes over the cause of the unfolding reaction, providing time- and spatially resolved atomic details of the folding mechanism of globular proteins.

59 BASIC BIOLOGICAL SCIENCES↗