Understanding protein structural dynamics is central to elucidating biological function and guiding therapeutic discovery. Hydrogen-deuterium exchange mass spectrometry (HDX-MS) typically offers peptide-level, and sometimes residue-level, time-dependent insights into protein structure, conformational dynamics and/or ligand binding. Yet, translating HDX-MS data into atomic-resolution insights and deriving mechanistic understanding remains a key challenge. Integrative strategies which utilize HDX-MS to inform computational modeling or simulations, traditionally leverage HDX-MS data with physics-based approaches through the calculation of protection factors models. Here, we developed MInt-HDX, a hybrid physics-based, machine- learning framework trained on differential HDX-MS signatures across 11 protein-ligand systems or 1032 individual peptides, using eXtreme Gradient Boosting (XGBoost) to guide small-molecule ligand docking and pose selection. By leveraging XGBoost-predicted interacting residues with three-dimensional clustering and convex-hull geometric algorithms, MInt-HDX first generates HDX-guided candidate docking sites in 3D for physics-based molecular docking and then, following docking, employs HDX-MS-informed XGBoost filtering and scoring functions for ligand- pose ranking. MInt-HDX was validated across 3 protein-ligand systems, consistently resulting in Ligand-RMSD within 3 Å of the crystallographic ligand conformation, individual steps of MInt- HDX were optimized and its overall performance was assessed against HDX-MS data quality factors and benchmarked against common physics-based and machine learning based docking approaches. Together, this work highlights how machine learning, informed by HDX-MS and aided by physics-based approaches, can bridge the gap between solution-phase HDX-MS data and structural modeling to accelerate protein-ligand discovery pipelines.
In Pseudomonas aeruginosa chronic infections, heme is a primary source of the essential micronutrient iron. The cytoplasmic heme-binding protein, PhuS, regulates extracellular heme flux through its interaction with the iron-regulated heme oxygenase (HemO). Additionally, in its apo-state, PhuS modulates iron homeostasis by transcriptionally regulating the prrF1,2 sRNA genes. These two functions are mutually exclusive and dependent on the conformational rearrangement of PhuS upon heme binding and coordination. Herein, we characterize a PhuS R25A variant that shows similar heme-binding kinetics and transfer of heme to HemO as PhuS WT, while DNA-binding to the prrF1 promoter is completely lost, successfully uncoupling the two functions. HDX-MS analysis revealed an overall decrease in conformational dynamics of apo-and holo-PhuS R25A compared with their WT counterparts, demonstrating the importance of conformational flexibility for DNA binding. qRT-PCR and Northern blot analysis comparing the phuSR25A allelic mutant strain to the PAO1 WT showed a significant decrease in PrrF and PrrH levels and revealed PhuS-dependent differences in regulation over PrrF1 and PrrF2, altering the relative ratio of these two sRNAs in a heme-specific manner that is distinct from iron. By removing its DNA-binding function, we elucidated the direct effects of PhuS binding on PrrF expression, separate from its effects on heme transfer and utilization. The contrasting effects on gene expression of the tandem sRNAs PrrF1 and PrrF2 in iron and heme and the resulting distinct mRNA profiles may allow the bacteria a fitness advantage in establishing chronic infection.
Mechanistic target of rapamycin (mTOR) is a key protein kinase that integrates various internal and external signals to control biological events including cell growth. Whereas substantial efforts were made to elucidate protein subunits interacting with mTOR, endogenous metabolite-mTOR interactions remain largely unknown. Using affinity protein purification and mass spectrometry, we identified direct binding of mTOR to 13- S -hydroxyoctadecadienoic acid (13- S -HODE) which is an oxygenated metabolite of linoleic acid, a polyunsaturated essential fatty acid. Interaction of 13- S -HODE with the catalytic ATP-binding domain of mTOR prevented its kinase activity in an ATP-competitive manner. Furthermore, either 13- S -HODE treatment or expression of arachidonate 15-lipoxygenase (ALOX15), an enzyme responsible for 13- S -HODE production, reduced mTOR signaling, thereby suppressing the growth of cancer cells as well as tumor xenografts. Our results highlight the importance of 13- S -HODE serving as a tumor suppressive, mTOR-inhibiting metabolite that links polyunsaturated fatty acid metabolism and the mTOR signaling in controlling cancer cell growth. ### Competing Interest Statement The authors have declared no competing interest.
Pseudomonas aeruginosa can acquire iron from heme via the heme assimilation system and Pseudomonas heme uptake (Phu) systems. Heme uptake is regulated at the metabolic level by the cytoplasmic protein PhuS that controls heme flux through a heme oxygenase HemO, releasing iron and biliverdin IXβ and IXδ. We have shown PhuS regulates extracellular heme flux, and in its apo-form transcriptionally regulates the iron and heme-dependent small RNAs (sRNAs) PrrF/PrrH. This mutual exclusivity of function is driven by conformational rearrangement of PhuS on heme binding. Herein, we show through a combination of EMSA and fluorescence anisotropy that mutation of the His-209 proximal ligand allows both apo- and holo-PhuS H209A to bind to the prrF1 promoter with significantly lower affinity when compared to PAO1 WT. Hydrogen deuterium exchange coupled to mass spectrometry revealed the apo- and holo-PhuS H209A structures are closer to each other than their WT counterparts and sample a conformational landscape between the apo- and holo-PhuS WT conformations, that is neither optimal for heme transfer nor DNA-binding. Furthermore, quantitative PCR and Western blot analysis of the phuSH209A allelic strain compared to PAO1 WT revealed an uncoupling of the PhuS-HemO dependent regulation of heme flux into the cell that abrogates the heme dependent regulation of the PrrF/PrrH sRNAs. The data supports a model where heme coordination through His-209 drives the conformational switch that determines mutual exclusivity in function of apo- and holo-PhuS. This dual function of PhuS is central to integrating extracellular heme utilization into the PrrF/PrrH sRNA regulatory network critical for P. aeruginosa adaptation within the host.
The dengue virus (DENV) NS5 protein, essential for viral RNA synthesis, is an attractive antiviral drug target. DENV NS5 interacts with the stem-loop A (SLA) promoter at the 5’-untranslated region of the viral genome to initiate negative-strand synthesis. However, the conformational dynamics of this interaction remains unclear. Our study explores the structural dynamics of DENV serotype 2 NS5 (DENV2 NS5) in complex with SLA, employing surface plasmon resonance (SPR), hydrogen-deuterium exchange mass spectrometry (HDX-MS), computational modeling, and cryoEM. Our findings reveal that DENV2 NS5 binds SLA in a closed conformation, with interdomain cooperation between its methyltransferase (MTase) and RNA-dependent RNA polymerase (RdRp) domains, critical for the interaction. SLA binding induces conformational changes in both domains, highlighting NS5’s multifunctional role in viral replication. Our cryoEM results visualizes the DENV2 NS5-SLA complex, confirming a conserved SLA binding across DENV serotypes and provides key insights for antiviral strategies targeting NS5’s conformational states.
N49P9.6-FR-LS and PGT121 are promising antibodies with significant therapeutic potential against HIV infection, but they are prone to precipitation at concentrations greater than 12 to 13 mg/mL. This study evaluates the influence of six excipients─arginine, alanine, sucrose, trehalose, methionine, and glutamate─on the biophysical stability of antibodies. We employed a comprehensive approach, combining computational mAb-excipient interaction analysis via the site-identification by ligand competitive saturation (SILCS) method with extensive experimental characterization. Our experimental matrix included viscosity measurements across temperature gradients, particle size distribution, zeta potential, pH value, and solution appearance, alongside a short-term stability product study at 30 °C and 65% relative humidity, with assessments at t0 (initial), t1 (14 days), and t2 (28 days). Results indicated that sucrose, arginine, alanine, and trehalose provided varying degrees of stabilization for both antibodies. Conversely, glutamate destabilized PGT121 but stabilized N49P9.6-FR-LS, while methionine had a negative effect on N49P9.6-FR-LS but a positive one on PGT121. SILCS-Biologics analysis suggested that stabilization by these excipients is linked to their ability to occupy regions involved in self-protein interactions. Debye-Hückel-Henry charge calculations further indicated that neutral excipients like sucrose and trehalose could alter mAb charges by affecting buffer binding, influencing aggregation propensity. These findings offer valuable insights for optimizing antibody formulations, ensuring enhanced product stability and therapeutic efficacy for HIV treatment.
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.
Arginyltransferase 1 (ATE1) catalyzes arginylation, an important posttranslational 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 posttranslational 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 those 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 facilitates the formation of a complex between ATE1 and tRNAArg, adding a new complexity to the ATE1 structure and providing new insights for future studies of ATE1 functions.
Dengue virus is the most prevalent arthropod-borne virus and there are no clinically approved antivirals to date. The non-structural 5 (NS5) protein is the largest protein encoded by flaviviruses including dengue, with an N-terminal methyltransferase (MTase) domain responsible for 5′ RNA capping, and a C-terminal RNA-dependent-RNA-polymerase (RdRp) domain responsible for de novo RNA synthesis. Stem Loop A (SLA) is an RNA element at the 5'-untranslated region which acts as a recognition motif for the initiation of RNA synthesis by NS5.
Interleukin-1 (IL-1)-family cytokines are potent modulators of inflammation, coordinating a vast array of immunological responses across innate and adaptive immune systems. Dysregulated IL-1-family cytokine signaling, however, is involved in a multitude of adverse health effects, such as chronic inflammatory conditions, autoimmune diseases, and cancer. Within the IL-1 family of cytokines, six-IL-1α, IL-1β, IL-33, IL-36α, IL-36β, and IL-36γ-require the IL-1 receptor accessory protein (IL-1RAcP) as their shared co-receptor. Common features of cytokine signaling include redundancy of signaling pathways, sharing of cytokines and receptors, pleiotropy of the cytokines themselves, and multifaceted immune responses. Accordingly, targeting multiple cytokines simultaneously is an emerging therapeutic strategy and can provide advantages over targeting a single cytokine pathway. Here, we show that two monoclonal antibodies, CAN10 and 3G5, which target IL-1RAcP for broad blockade of all associated cytokines, do so through distinct mechanisms and provide therapeutic opportunities for the treatment of inflammatory diseases.
Conventional molecular dynamics (cMD) simulations depict the structure and dynamics of biomolecules with atom-level resolution. Due to timescale limitations, forcefield inaccuracies, and inherent biases in high-resolution structure determination techniques, however, cMD may inaccurately reflect protein behavior in solution. This discrepancy can be further complicated by the conformational heterogeneity of the molecule. For example, calcium-free calmodulin retains a compact, well-ordered structure throughout long cMD simulations, despite experiments illustrating the C-terminal domain's flexibility in solution.
Antibodies are thought to play a major role in protection against human norovirus infection. Mouse humoral responses closely mimic those of humans; thus, mouse models are used to characterize norovirus epitopes on the major viral capsid protein, VP1. We have developed a panel of mouse monoclonal antibodies (mAbs) produced against the last pandemic variant to emerge, Sydney 2012. While most mAbs (25/44) were mapped to variable antigenic sites on VP1, 19 of the mAbs were cross-reactive against multiple genotypes or GII.4 variants. Most (12/19) of the cross-reactive mAbs bound to the Shell domain and were cross-reactive with different GII noroviruses. Interestingly, mAb 30A11 exhibited cross-reactivity against all tested norovirus genotypes (GI, GII, GIV, and GIX). This mAb was mapped to a highly conserved region of the Shell domain (51PIDPWII57) using peptide ELISA and immunofluorescence. Of those mapping to the Protruding (P) domain, two (19C10 and 14B11) showed cross-reactivity with GII noroviruses. Using hydrogen-deuterium exchange mass spectrometry, we mapped 19C10 to a conserved region of the P domain near the P/Shell interface, which explains its cross-reactivity with different GII noroviruses and lack of histo-blood group antigen-blocking activity. Binding and mutational analyses showed that residues 518, 519, and 525 are important for 19C10 and 14B11 epitope recognition. While the antibodies described here are mostly non-neutralizing, they can be useful tools for research and diagnostics of noroviruses. The role of non-neutralizing, cross-reactive antibodies targeting different areas of the viral capsid merits further research to facilitate our understanding of immunity to norovirus infection and disease. IMPORTANCE:To gain insights into the overall immune responses to human norovirus, we characterized non-neutralizing, cross-reactive monoclonal antibodies (mAbs) developed against a pandemic GII.4 norovirus. We determined the binding epitope of an antibody that exhibited cross-reactivity against all tested noroviruses, which makes it a useful tool for research and diagnostics. The epitope of two additional non-neutralizing mAbs was mapped to a less conserved region on the viral capsid protein, explaining their cross-reactivity patterns. Often overlooked, the role of non-neutralizing, cross-reactive mAbs merits further research to facilitate our understanding of immunity to norovirus infection and disease.
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.
P38 MAP kinases, members of the mitogen-activated protein kinase (MAPK) family, play an essential role in response to stress, including that occurring in inflammation-related lung injury. Notably, body temperatures associated with fever enhance injury while hypothermia tends to reduce injury. It has been observed that temperature dependent lung injury can be diminished through the administration of inhibitors targeting p38 MAP kinases. Among the major p38 isoforms, p38⍺ is pro-inflammatory while p38β is cytoprotective.
Ran-binding domain-containing protein 2 (ZRANB2) is a zinc finger (ZF) protein that plays a key role in alternative splicing. ZRANB2 is composed of two ZF domains that contain four invariant cysteine residues per domain. ZRANB2 binds RNA targets that contain AGGUAA sequence motifs. Three constructs of ZRANB2, ZRANB2-ZF1 (first ZF domain), ZRANB2-ZF2 (second ZF domain), and ZRANB2-2D (both ZF domains), were isolated in the apo form and shown to bind Zn(II) via UV-visible-monitored competitive titrations with Co(II) as a spectroscopic probe. Zn binding to each construct led to the adoption of a limited secondary structure of each domain, as measured by circular dichroism (CD). Hydrogen-deuterium exchange coupled with mass spectrometry (HDX-MS) of the two-domain construct, ZRANB2-2D, revealed that both ZF domains adopt a more rigid structure upon Zn binding. Zn binding to the first ZF domain resulted in a greater decrease in the conformational dynamics than Zn binding to the second ZF domain. RNA binding to TRA2B pre-mRNA, a physiological splicing target, was measured by fluorescence anisotropy (FA), and high-affinity RNA binding was found to require Zn coordination to both domains. HDX-MS of ZRANB2-2D with TRA2B RNA as well as two optimized RNA sequences that contain a single and double AGGUAA hexamer revealed additional protection from H/D exchange for ZRANB2 in the presence of RNA. Here, greater protection was observed for the second ZF of ZRANB2-2D, suggesting a larger effect on conformational dynamics. A model for zinc-mediated RNA binding of ZRANB2 is proposed.