Integration of Hydrogen-Deuterium Exchange Mass Spectrometry with Molecular Dynamics Simulations and Ensemble Reweighting Enables High Resolution Protein-Ligand Modeling.

Hydrogen-Deuterium exchange mass spectrometry's (HDX-MS) utility in identifying and characterizing protein-small molecule interaction sites has been established. The regions that are seen to be protected from exchange upon ligand binding indicate regions that may be interacting with the ligand, giving a qualitative understanding of the ligand binding pocket. However, quantitatively deriving an accurate high-resolution structure of the protein-ligand complex from the HDX-MS data remains a challenge, often limiting its use in applications such as small molecule drug design. Recent efforts have focused on the development of methods to quantitatively model Hydrogen-Deuterium exchange (HDX) data from computationally modeled structures to garner atomic level insights from peptide-level resolution HDX-MS. One such method, HDX ensemble reweighting (HDXer), employs maximum entropy reweighting of simulated HDX data to experimental HDX-MS to model structural ensembles. In this study, we implement and validate a workflow which quantitatively leverages HDX-MS data to accurately model protein-small molecule ligand interactions. To that end, we employ a strategy combining computational protein-ligand docking, molecular dynamics simulations, HDXer, and dimensional reduction and clustering approaches to extract high-resolution drug binding poses that most accurately conform with HDX-MS data. We apply this workflow to model the interaction of ERK2 and FosA with small molecule compounds and inhibitors they are known to bind. In five out of six of the protein-ligand pairs tested, the HDX derived protein-ligand complexes result in a ligand root-mean-square deviation (RMSD) within 2.5 Å of the known crystal structure ligand.

Identification of an intrinsically disordered region (IDR) in arginyltransferase 1 (ATE1).

Arginyltransferase 1 (ATE1) catalyzes arginylation, an important post-translational modification (PTM) in eukaryotes that plays a critical role in cellular homeostasis. The disruption of ATE1 function is implicated in mammalian neurodegenerative disorders and cardiovascular maldevelopment, while post-translational arginylation has also been linked to the activities of several important human viruses such as SARS-CoV-2 and HIV. Despite the known significance of ATE1 in mammalian cellular function, past biophysical studies of this enzyme have mainly focused on yeast ATE1, leaving the mechanism of arginylation in mammalian cells unclear. In this study, we sought to structurally and biophysically characterize mouse (Mus musculus) ATE1. Using size-exclusion chromatography (SEC), small angle X-ray scattering (SAXS), and hydrogen deuterium exchange mass spectrometry (HDX-MS), assisted by AlphaFold modeling, we found that mouse ATE1 is structurally more complex than yeast ATE1. Importantly, our data indicate the existence of an intrinsically disordered region (IDR) in all mouse ATE1 splice variants. However, comparative HDX-MS analyses show that yeast ATE1 does not have such an IDR, consistent with prior X-ray, cryo-EM, and SAXS analyses. Furthermore, bioinformatics approaches reveal that mammalian ATE1 sequences, as well as in a large majority of other eukaryotes, contain an IDR-like sequence positioned in proximity to the ATE1 GNAT active-site fold. Computational analysis suggests that the IDR likely facilitates the formation of the complex between ATE1 and tRNAArg, adding a new complexity to ATE1 structure and providing new insights for future studies of ATE1 functions.

Coupling of zinc and GTP binding drives G-domain folding in Acinetobacter baumannii ZigA.

COG0523 proteins, also known as nucleotide-dependent metallochaperones, are a poorly understood class of small P-loop G3E GTPases. Multiple family members play critical roles in bacterial pathogen survival during an infection as part of the adaptive response to host-mediated "nutritional immunity." Our understanding of the structure, dynamics, and molecular-level function of COG0523 proteins, apart from the eukaryotic homolog, Zng1, remains in its infancy. Here, we use X-ray absorption spectroscopy to establish that Acinetobacter baumannii (Ab) ZigA coordinates ZnII using all three cysteines derived from the invariant CXCC motif to form an S3(N/O) coordination complex, a feature inconsistent with the ZnII-bound crystal structure of a distantly related COG0523 protein of unknown function from Escherichia coli, EcYjiA. The binding of ZnII and guanine nucleotides is thermodynamically linked in AbZigA, and this linkage is more favorable for the substrate GTP relative to the product GDP. Part of this coupling originates with nucleotide-induced stabilization of the G-domain tertiary structure as revealed by global thermodynamics measurements and hydrogen-deuterium exchange mass spectrometry (HDX-MS). HDX-MS also reveals that the HDX behavior of the G2 (switch 1) loop is highly sensitive to nucleotide status and becomes more exchange labile in the GDP (product)-bound state. Significant long-range perturbation of local stability in both the G-domain and the C-terminal domain define a candidate binding pocket for a client protein that appears sensitive to nucleotide status (GDP versus GTP). We place these new insights into the structure, dynamics, and energetics of intermolecular metal transfer into the context of a model for AbZigA metallochaperone function.

A fusion of the Bacteroides fragilis ferrous iron import proteins reveals a role for FeoA in stabilizing GTP-bound FeoB.

Iron is an essential element for nearly all organisms, and under anoxic and/or reducing conditions, Fe2+ is the dominant form of iron available to bacteria. The ferrous iron transport (Feo) system is the primary prokaryotic Fe2+ import machinery, and two constituent proteins (FeoA and FeoB) are conserved across most bacterial species. However, how FeoA and FeoB function relative to one another remains enigmatic. In this work, we explored the distribution of feoAB operons encoding a fusion of FeoA tethered to the N-terminal, G-protein domain of FeoB via a connecting linker region. We hypothesized that this fusion poises FeoA to interact with FeoB to affect function. To test this hypothesis, we characterized the soluble NFeoAB fusion protein from Bacteroides fragilis, a commensal organism implicated in drug-resistant infections. Using X-ray crystallography, we determined the 1.50-Å resolution structure of BfFeoA, which adopts an SH3-like fold implicated in protein-protein interactions. Using a combination of structural modeling, small-angle X-ray scattering, and hydrogen-deuterium exchange mass spectrometry, we show that FeoA and NFeoB interact in a nucleotide-dependent manner, and we mapped the protein-protein interaction interface. Finally, using guanosine triphosphate (GTP) hydrolysis assays, we demonstrate that BfNFeoAB exhibits one of the slowest known rates of Feo-mediated GTP hydrolysis that is not potassium-stimulated. Importantly, truncation of FeoA from this fusion demonstrates that FeoA-NFeoB interactions function to stabilize the GTP-bound form of FeoB. Taken together, our work reveals a role for FeoA function in the fused FeoAB system and suggests a function for FeoA among prokaryotes.

Modeling the native ensemble of PhuS using enhanced sampling MD and HDX-ensemble reweighting.

The cytoplasmic heme binding protein from Pseudomonas aeruginosa, PhuS, plays two essential roles in regulating heme uptake and iron homeostasis. First, PhuS shuttles exogenous heme to heme oxygenase (HemO) for degradation and iron release. Second, PhuS binds DNA and modulates the transcription of the prrF/H small RNAs (sRNAs) involved in the iron-sparing response. Heme binding to PhuS regulates this dual function, as the unliganded form binds DNA, whereas the heme-bound form binds HemO. Crystallographic studies revealed nearly identical structures for apo- and holo-PhuS, and yet numerous solution-based measurements indicate that heme binding is accompanied by large conformational rearrangements. In particular, hydrogen-deuterium exchange mass spectrometry (HDX-MS) of apo- versus holo-PhuS revealed large differences in deuterium uptake, notably in α-helices 6, 7, and 8 (α6,7,8), which contribute to the heme binding pocket. These helices were mostly labile in apo-PhuS but largely protected in holo-PhuS. In contrast, in silico-predicted deuterium uptake levels of α6,7,8 from molecular dynamics (MD) simulations of the apo- and holo-PhuS structures are highly similar, consistent only with the holo-PhuS HDX-MS data. To rationalize this discrepancy between crystal structures, simulations, and observed HDX-MS, we exploit a recently developed computational approach (HDXer) that fits the relative weights of conformational populations within an ensemble of structures to conform to a target set of HDX-MS data. Here, a combination of enhanced sampling MD, HDXer, and dimensionality reduction analysis reveals an apo-PhuS conformational landscape in which α6, 7, and 8 are significantly rearranged compared to the crystal structure, including a loss of secondary structure in α6 and the displacement of α7 toward the HemO binding interface. Circular dichroism analysis confirms the loss of secondary structure, and the extracted ensembles of apo-PhuS and of heme-transfer-impaired H212R mutant, are consistent with known heme binding and transfer properties. The proposed conformational landscape provides structural insights into the modulation by heme of the dual function of PhuS.

Current Trends and Limitations in Dengue Antiviral Research.

Dengue is the most prevalent arthropod-borne viral disease worldwide and affects approximately 2.5 billion people living in over 100 countries. Increasing geographic expansion of Aedes aegypti mosquitoes (which transmit the virus) has made dengue a global health concern. There are currently no approved antivirals available to treat dengue, and the only approved vaccine used in some countries is limited to seropositive patients. Treatment of dengue, therefore, remains largely supportive to date; hence, research efforts are being intensified for the development of antivirals. The nonstructural proteins, 3 and 5 (NS3 and NS5), have been the major targets for dengue antiviral development due to their indispensable enzymatic and biological functions in the viral replication process. NS5 is the largest and most conserved nonstructural protein encoded by flaviviruses. Its multifunctionality makes it an attractive target for antiviral development, but research efforts have, this far, not resulted in the successful development of an antiviral targeting NS5. Increase in structural insights into the dengue NS5 protein will accelerate drug discovery efforts focused on NS5 as an antiviral target. In this review, we will give an overview of the current state of therapeutic development, with a focus on NS5 as a therapeutic target against dengue.

Interpreting Hydrogen-Deuterium Exchange Experiments with Molecular Simulations: Tutorials and Applications of the HDXer Ensemble Reweighting Software [Article v1.0].

Hydrogen-deuterium exchange (HDX) is a comprehensive yet detailed probe of protein structure and dynamics and, coupled to mass spectrometry, has become a powerful tool for investigating an increasingly large array of systems. Computer simulations are often used to help rationalize experimental observations of exchange, but interpretations have frequently been limited to simple, subjective correlations between microscopic dynamical fluctuations and the observed macroscopic exchange behavior. With this in mind, we previously developed the HDX ensemble reweighting approach and associated software, HDXer, to aid the objective interpretation of HDX data using molecular simulations. HDXer has two main functions; first, to compute H-D exchange rates that describe each structure in a candidate ensemble of protein structures, for example from molecular simulations, and second, to objectively reweight the conformational populations present in a candidate ensemble to conform to experimental exchange data. In this article, we first describe the HDXer approach, theory, and implementation. We then guide users through a suite of tutorials that demonstrate the practical aspects of preparing experimental data, computing HDX levels from molecular simulations, and performing ensemble reweighting analyses. Finally we provide a practical discussion of the capabilities and limitations of the HDXer methods including recommendations for a user's own analyses. Overall, this article is intended to provide an up-to-date, pedagogical counterpart to the software, which is freely available at https://github.com/Lucy-Forrest-Lab/HDXer.

Gallium(III)-Salophen as a Dual Inhibitor of Pseudomonas aeruginosa Heme Sensing and Iron Acquisition.

Pseudomonas aeruginosa is an opportunistic bacterium that causes life-threatening infections in immunocompromised patients. In infection, it uses heme as a primary iron source and senses the availability of exogenous heme through the heme assimilation system (Has), an extra cytoplasmic function σ-factor system. A secreted hemophore HasAp scavenges heme and, upon interaction with the outer-membrane receptor HasR, activates a signaling cascade, which in turn creates a positive feedback loop critical for sensing and adaptation within the host. The ability to sense and respond to heme as an iron source contributes to virulence. Consequently, the inhibition of this system will lead to a disruption in iron homeostasis, decreasing virulence. We have identified a salophen scaffold that successfully inhibits the activation of the Has signaling system while simultaneously targeting iron uptake via xenosiderophore receptors. We propose this dual mechanism wherein free Ga3+-salophen reduces growth through uptake and iron mimicry. A dual mechanism targeting extracellular heme signaling and uptake together with Ga3+-induced toxicity following active Ga3+salophen uptake provides a significant therapeutic advantage while reducing the propensity to develop resistance.

Insights on the Conformational Ensemble of Cyt C Reveal a Compact State during Peroxidase Activity.

Cytochrome c (cyt c) is known for its role in the electron transport chain but transitions to a peroxidase-active state upon exposure to oxidative species. The peroxidase activity ultimately results in the release of cyt c into the cytosol for the engagement of apoptosis. The accumulation of oxidative modifications that accompany the onset of the peroxidase function are well-characterized. However, the concurrent structural and conformational transitions of cyt c remain undercharacterized. Fast photochemical oxidation of proteins (FPOP) coupled with mass spectrometry is a protein footprinting technique used to structurally characterize proteins. FPOP coupled with native ion mobility separation shows that exposure to H2O2 results in the accumulation of a compact state of cyt c. Subsequent top-down fragmentation to localize FPOP modifications reveals changes in heme coordination between conformers. A time-resolved functional assay suggests that this compact conformer is peroxidase active. Altogether, combining FPOP, ion mobility separation, and top-down and bottom-up mass spectrometry allows us to discern individual conformations in solution and obtain a better understanding of the conformational ensemble and structural transitions of cyt c as it transitions from a respiratory role to a proapoptotic role.

Monitoring dendrimer conformational transition using 19 F and 1 H2 O NMR.

The conformational transition of a fluorinated amphiphilic dendrimer is monitored by the 1 H signal from water, alongside the 19 F signal from the dendrimer. High-field NMR data (chemical shift δ, self-diffusion coefficient D, longitudinal relaxation rate R1 , and transverse relaxation rate R2 ) for both dendrimer (19 F) and water (1 H) match each other in detecting the conformational transition. Among all parameters for both nuclei, the water proton transverse-relaxation rate R2 (1 H2 O) displays the highest relative scale of change upon conformational transition of the dendrimer. Hydrogen/deuterium-exchange mass spectrometry reveals that the compact form of the dendrimer has slower proton exchange with water than the extended form. This result suggests that the sensitivity of R2 (1 H2 O) toward dendrimer conformation originates, at least partially, from the difference in proton exchange efficiency between different dendrimer conformations. Finally, we also demonstrated that this conformational transition could be conveniently monitored using a low-field benchtop NMR spectrometer via R2 (1 H2 O). The 1 H2 O signal thus offers a simple way to monitor structural changes of macromolecules using benchtop time-domain NMR.

Conformational transition of a non-associative fluorinated amphiphile in aqueous solution. II. Conformational transition vs. supramolecular assembly.

Unlike many known amphiphiles, the fluorinated amphiphilic dendrimer studied in this work demonstrated a concentration-dependent conformational transition rather than micellization or assembly. Hydrophobic and hydrophilic interactions with water were suggested as the most probable driving force of this transition. This assumption was consistent with the observed 19F chemical shift changes of the dendrimer compared to a known micelle-forming fluorinated amphiphile. Since water is an important factor in the process, trends of the concentration-dependent changes in water proton transverse relaxation rate served as an indicator of structural changes and/or supramolecular assembly. The conformational transition process was also confirmed by ion-mobility mass-spectrometry. We suggested that structural features, namely, steric hindrances, prevented the micellization/assembly of the dendrimer of this study. This conclusion might inform the approach to develop novel unconventional amphiphiles.

Remodeling KRAS.

Mutations in members of the RAS family of small GTPases have been associated with numerous human cancers. However, RAS family members are notoriously difficult to target. In this issue of Structure, Lu et al. (2017) examine the effects of two compounds with distinct chemical scaffolds on the structure and dynamics of an oncogenic KRAS mutant, thus highlighting the usefulness of HDX-MS for drug development.

Common coding variant in SERPINA1 increases the risk for large artery stroke.

Large artery atherosclerotic stroke (LAS) shows substantial heritability not explained by previous genome-wide association studies. Here, we explore the role of coding variation in LAS by analyzing variants on the HumanExome BeadChip in a total of 3,127 cases and 9,778 controls from Europe, Australia, and South Asia. We report on a nonsynonymous single-nucleotide variant in serpin family A member 1 (SERPINA1) encoding alpha-1 antitrypsin [AAT; p.V213A; P = 5.99E-9, odds ratio (OR) = 1.22] and confirm histone deacetylase 9 (HDAC9) as a major risk gene for LAS with an association in the 3'-UTR (rs2023938; P = 7.76E-7, OR = 1.28). Using quantitative microscale thermophoresis, we show that M1 (A213) exhibits an almost twofold lower dissociation constant with its primary target human neutrophil elastase (NE) in lipoprotein-containing plasma, but not in lipid-free plasma. Hydrogen/deuterium exchange combined with mass spectrometry further revealed a significant difference in the global flexibility of the two variants. The observed stronger interaction with lipoproteins in plasma and reduced global flexibility of the Val-213 variant most likely improve its local availability and reduce the extent of proteolytic inactivation by other proteases in atherosclerotic plaques. Our results indicate that the interplay between AAT, NE, and lipoprotein particles is modulated by the gate region around position 213 in AAT, far away from the unaltered reactive center loop (357-360). Collectively, our findings point to a functionally relevant balance between lipoproteins, proteases, and AAT in atherosclerosis.

Ligand-induced allostery in the interaction of the Pseudomonas aeruginosa heme binding protein with heme oxygenase.

A heme-dependent conformational rearrangement of the C-terminal domain of heme binding protein (PhuS) is required for interaction with the iron-regulated heme oxygenase (HemO). Herein, we further investigate the underlying mechanism of this conformational rearrangement and its implications for heme transfer via site-directed mutagenesis, resonance Raman (RR), hydrogen-deuterium exchange MS (HDX-MS) methods, and molecular dynamics (MD). HDX-MS revealed that the apo-PhuS C-terminal α6/α7/α8-helices are largely unstructured, whereas the apo-PhuS H212R variant showed an increase in structure within these regions. The increased rate of heme association with apo-PhuS H212R compared with the WT and lack of a detectable five-coordinate high-spin (5cHS) heme intermediate are consistent with a more folded and less dynamic C-terminal domain. HDX-MS and MD of holo-PhuS indicate an overall reduction in molecular flexibility throughout the protein, with significant structural rearrangement and protection of the heme binding pocket. We observed slow cooperative unfolding/folding events within the C-terminal helices of holo-PhuS and the N-terminal α1/α2-helices that are dampened or eliminated in the holo-PhuS H212R variant. Chemical cross-linking and MALDI-TOF MS mapped these same regions to the PhuS:HemO protein-protein interface. We previously proposed that the protein-protein interaction induces conformational rearrangement, promoting a ligand switch from His-209 to His-212 and triggering heme release to HemO. The reduced conformational freedom of holo-PhuS H212R combined with the increase in entropy and decrease in heme transfer on interaction with HemO further support this model. This study provides significant insight into the role of protein dynamics in heme binding and release in bacterial heme transport proteins.

Conformational dynamics of a neurotransmitter:sodium symporter in a lipid bilayer.

Neurotransmitter:sodium symporters (NSSs) are integral membrane proteins responsible for the sodium-dependent reuptake of small-molecule neurotransmitters from the synaptic cleft. The symporters for the biogenic amines serotonin (SERT), dopamine (DAT), and norepinephrine (NET) are targets of multiple psychoactive agents, and their dysfunction has been implicated in numerous neuropsychiatric ailments. LeuT, a thermostable eubacterial NSS homolog, has been exploited as a model protein for NSS members to canvass the conformational mechanism of transport with a combination of X-ray crystallography, cysteine accessibility, and solution spectroscopy. Despite yielding remarkable insights, these studies have primarily been conducted with protein in the detergent-solubilized state rather than embedded in a membrane mimic. In addition, solution spectroscopy has required site-specific labeling of nonnative cysteines, a labor-intensive process occasionally resulting in diminished transport and/or binding activity. Here, we overcome these limitations by reconstituting unlabeled LeuT in phospholipid bilayer nanodiscs, subjecting them to hydrogen-deuterium exchange coupled with mass spectrometry (HDX-MS), and facilitating interpretation of the data with molecular dynamics simulations. The data point to changes of accessibility and dynamics of structural elements previously implicated in the transport mechanism, in particular transmembrane helices (TMs) 1a and 7 as well as extracellular loops (ELs) 2 and 4. The results therefore illuminate the value of this strategy for interrogating the conformational mechanism of the more clinically significant mammalian membrane proteins including SERT and DAT, neither of which tolerates complete removal of endogenous cysteines, and whose activity is heavily influenced by neighboring lipids.

The glutamate effect on DNA binding by pol I DNA polymerases: osmotic stress and the effective reversal of salt linkage.

The significant enhancing effect of glutamate on DNA binding by Escherichia coli nucleic acid binding proteins has been extensively documented. Glutamate has also often been observed to reduce the apparent linked ion release (Deltan(ions)) upon DNA binding. In this study, it is shown that the Klenow and Klentaq large fragments of the Type I DNA polymerases from E. coli and Thermus aquaticus both display enhanced DNA binding affinity in the presence of glutamate versus chloride. Across the relatively narrow salt concentration ranges often used to obtain salt linkage data, Klenow displays an apparently decreased Deltan(ions) in the presence of Kglutamate, while Klentaq appears not to display an anion-specific effect on Deltan(ions). Osmotic stress experiments reveal that DNA binding by Klenow and Klentaq is associated with the release of approximately 500 to 600 waters in the presence of KCl. For both proteins, replacing chloride with glutamate results in a 70% reduction in the osmotic-stress-measured hydration change associated with DNA binding (to approximately 150-200 waters released), suggesting that glutamate plays a significant osmotic role. Measurements of the salt-DNA binding linkages were extended up to 2.5 M Kglutamate to further examine this osmotic effect of glutamate, and it is observed that a reversal of the salt linkage occurs above 800 mM for both Klenow and Klentaq. Salt-addition titrations confirm that an increase of [Kglutamate] beyond 1 M results in rebinding of salt-displaced polymerase to DNA. These data represent a rare documentation of a reversed ion linkage for a protein-DNA interaction (i.e., enhanced binding as salt concentration increases). Nonlinear linkage analysis indicates that this unusual behavior can be quantitatively accounted for by a shifting balance of ionic and osmotic effects as [Kglutamate] is increased. These results are predicted to be general for protein-DNA interactions in glutamate salts.