Structure of ubiquitin refined at 1.8 A resolution

An internal linker and pH biosensing by phosphatidylinositol 5-phosphate regulate the function of the ESCRT-0 component TOM1

Target of Myb1 (TOM1) facilitates the transport of endosomal ubiquitinated proteins destined for lysosomal degradation; however, the mechanisms regulating TOM1 during this process remain unknown. Here, we identified an adjacent DXXLL motif-containing region to the TOM1 VHS domain, which enhances its affinity for ubiquitin and can be modulated by phosphorylation. TOM1 is an endosomal phosphatidylinositol 5-phosphate (PtdIns5P) effector under Shigella flexneri infection. We pinpointed a consensus PtdIns5P-binding motif in the VHS domain. We show that PtdIns5P binding by TOM1 is pH-dependent, similarly observed in its binding partner TOLLIP. Under acidic conditions, TOM1 retained its complex formation with TOLLIP, but was unable to bind ubiquitin. S. flexneri infection inhibits pH-dependent endosomal maturation, leading to reduced protein degradation. We propose a model wherein pumping of H+ to the cytosolic side of endosomes contributes to the accumulation of TOM1, and possibly TOLLIP, at these sites, thereby promoting PtdIns5P- and pH-dependent signaling, facilitating bacterial survival.

Introductory Tutorials for Simulating Protein Dynamics with GROMACS

Atomistic molecular dynamics (MD) simulations have become an indispensable tool for investigating the structure, dynamics, and energetics of biomolecules. Continual optimization of software algorithms and hardware has enabled investigators to access biologically relevant time scales in feasible amounts of computing time. Given the widespread use and utility of MD simulations, there is considerable interest in learning essential skills in performing them. Here, we present a set of introductory tutorials for performing MD simulations of proteins in the popular, open-source GROMACS package. Three exercises are detailed, including simulating a single protein, setting up a protein complex, and performing umbrella sampling simulations to model the unfolding of a short polypeptide. Essential features and input settings are illustrated throughout. The purpose of these tutorials is to provide new users with a general understanding of foundational workflows, from which they can design their own simulations.

Top-Down Structural Characterization of Native Ubiquitin Combining Solution-Stable Isotope Labeling, Trapped Ion Mobility Spectrometry, and Tandem Electron Capture Dissociation Mass Spectrometry

Solution-phase hydrogen/deuterium exchange (HDX) coupled to native ion mobility spectrometry mass spectrometry (IMS-MS) can provide complementary structural information about the conformational dynamics of biological molecules. In the present work, the solution-stable isotope labeling (SIL) combined with trapped ion mobility spectrometry (TIMS) in tandem with top-down electron capture dissociation (ECD) is illustrated for the structural characterization of the solution native states of ubiquitin. Four different ubiquitin electrospray solution conditions: (i) single-tip nondeuterated, (ii) theta tip for online SIL HDX, (iii) single-tip SIL-deuterated, and (iv) theta tip for online SIL H/D back exchange (HDbX), were investigated to assess the H/D exchange reactivities of native ubiquitin. The combination of TIMS and ECD in a q-ToF MS instrument allowed for additional inspection of gas-phase HDbX added by top-down fragmentation, revealing the exposed and protected residues with limited scrambling effects (e.g., intramolecular H/D migration). A native charge state distribution (5+ to 7+) and TIMS profiles were observed under the single-tip nondeuterated solution conditions. Mass shift distributions of ∼40, ∼104, and ∼87D were observed when incorporating deuterium for online SIL HDX, SIL HDX, and online SIL HDbX, respectively, while retaining similar conformational states. ECD fragmentation allowed for the localization of the deuterated labeled residues of the peptide fragments, with a sequence coverage of ∼90%, for each of the ubiquitin solution condition. Changes in the TIMS trapping time settings (∼70 to ∼795 ms) were used to determine the H/D back exchange dynamics of native ubiquitin. HDbX-TIMS-q-ECD-MS/MS exhibited H/D back exchanges in the six-residue C-terminal tail as well as around Lys6, Lys11, Lys33, Lys48, and Lys63 residues, indicating that these regions are the most exposed area (less protected hydrogens) of ubiquitin as compared to the rest of the core residues that adopt a compact β-grasp fold (protected hydrogens), which was consistent with the accessible surface area of ubiquitin. The present data highlight for the first time consistency between the solution HDX and gas-phase HDbX-TIMS data for native studies.

The Ubiquitin Tale: Current Strategies and Future Challenges

Ubiquitin (Ub) is often considered a structurally conserved protein. Ubiquitination plays a prominent role in the regulation of physiological pathways. Since the first mention of Ub in protein degradation pathways, a plethora of nonproteolytic functions of this post-translational modification have been identified and investigated in detail. In addition, several other structurally and functionally related proteins have been identified and investigated for their Ub-like structures and functions. Ubiquitination and Ub-like modifications play vital roles in modulating the pathways involved in crucial biological processes and thus affect the global proteome. In this Review, we provide a snapshot of pathways, substrates, diseases, and novel therapeutic targets that are associated with ubiquitination or Ub-like modifications. In the past few years, a large number of proteomic studies have identified pools of ubiquitinated proteins (ubiquitylomes) involved or induced in healthy or stressed conditions. These comprehensive studies involving identification of new ubiquitination substrates and sites contribute enormously to our understanding of ubiquitination in more depth. However, with the current tools, there are certain limitations that need to be addressed. We review recent technological advancements in ubiquitylomic studies and their limitations and challenges. Overall, large-scale ubiquitylomic studies contribute toward understanding global ubiquitination in the contexts of normal and disease conditions.

Ubiquitinated histone H2B as gatekeeper of the nucleosome acidic patch

Monoubiquitination of histones H2B-K120 (H2BK120ub) and H2A-K119 (H2AK119ub) play opposing roles in regulating transcription and chromatin compaction. H2BK120ub is a hallmark of actively transcribed euchromatin, while H2AK119ub is highly enriched in transcriptionally repressed heterochromatin. Whereas H2BK120ub is known to stimulate the binding or activity of various chromatin-modifying enzymes, this post-translational modification (PTM) also interferes with the binding of several proteins to the nucleosome H2A/H2B acidic patch via an unknown mechanism. Here, we report cryoEM structures of an H2BK120ub nucleosome showing that ubiquitin adopts discrete positions that occlude the acidic patch. Molecular dynamics simulations show that ubiquitin remains stably positioned over this nucleosome region. By contrast, our cryoEM structures of H2AK119ub nucleosomes show ubiquitin adopting discrete positions that minimally occlude the acidic patch. Consistent with these observations, H2BK120ub, but not H2AK119ub, abrogates nucleosome interactions with acidic patch-binding proteins RCC1 and LANA, and single-domain antibodies specific to this region. Our results suggest a mechanism by which H2BK120ub serves as a gatekeeper to the acidic patch and point to distinct roles for histone H2AK119 and H2BK120 ubiquitination in regulating protein binding to nucleosomes.

Streamlining NMR Chemical Shift Predictions for Intrinsically Disordered Proteins: Design of Ensembles with Dimensionality Reduction and Clustering

By merging advanced dimensionality reduction (DR) and clustering algorithm (CA) techniques, our study advances the sampling procedure for predicting NMR chemical shifts (CS) in intrinsically disordered proteins (IDPs), making a significant leap forward in the field of protein analysis/modeling. We enhance NMR CS sampling by generating clustered ensembles that accurately reflect the different properties and phenomena encapsulated by the IDP trajectories. This investigation critically assessed different rapid CS predictors, both neural network (e.g., Sparta+ and ShiftX2) and database-driven (ProCS-15), and highlighted the need for more advanced quantum calculations and the subsequent need for more tractable-sized conformational ensembles. Although neural network CS predictors outperformed ProCS-15 for all atoms, all tools showed poor agreement with HN CSs, and the neural network CS predictors were unable to capture the influence of phosphorylated residues, highly relevant for IDPs. This study also addressed the limitations of using direct clustering with collective variables, such as the widespread implementation of the GROMOS algorithm. Clustered ensembles (CEs) produced by this algorithm showed poor performance with chemical shifts compared to sequential ensembles (SEs) of similar size. Instead, we implement a multiscale DR and CA approach and explore the challenges and limitations of applying these algorithms to obtain more robust and tractable CEs. The novel feature of this investigation is the use of solvent-accessible surface area (SASA) as one of the fingerprints for DR alongside previously investigated α carbon distance/angles or ϕ/ψ dihedral angles. The ensembles produced with SASA tSNE DR produced CEs better aligned with the experimental CS of between 0.17 and 0.36 r2 (0.18-0.26 ppm) depending on the system and replicate. Furthermore, this technique produced CEs with better agreement than traditional SEs in 85.7% of all ensemble sizes. This study investigates the quality of ensembles produced based on different input features, comparing latent spaces produced by linear vs nonlinear DR techniques and a novel integrated silhouette score scanning protocol for tSNE DR.

Deep mutational scanning of EccD3 reveals the molecular basis of its essentiality in the mycobacterium ESX secretion system

Tuberculosis remains the deadliest infectious disease in the world and requires novel therapeutic targets. The ESX-3 secretion system, which is essential for iron and zinc homeostasis and thus M. tuberculosis survival, is a promising target. In this study, we perform a deep mutational scan on the ESX-3 core protein EccD3 in the model organism M. smegmatis. We systematically investigated the functional roles of 145 residues across the soluble ubiquitin-like domain, the conformationally distinct flexible linker, and selected transmembrane helices of EccD3. Our data combined with structural comparisons to ESX-5 complexes support a model where EccD3 stabilizes the complex, with the hinge motif within the linker being particularly sensitive to disruption. Our study is the first deep mutational scan in mycobacteria, which could help guide drug development toward novel treatment of tuberculosis. This study underscores the importance of context-specific mutational analyses for discovering essential protein interactions within mycobacterial systems.

Competitive Effects of Anions on Protein Solvation by Aqueous Ionic Liquids

The present study utilizes molecular dynamics simulations to examine how different anions compete for protein solvation in aqueous solutions of ionic liquids (ILs). Ubiquitin is used as model protein and studied in IL mixtures sharing the same cation, 1-ethyl-3-methylimidazolium (EMIM), and two different anions in the same solution, from combinations of dicyanamide (DCA), chloride (Cl), nitrate (NO3), and tetrafluoroborate (BF4). Our findings reveal that specific interactions between anions and the protein are paramount in IL solvation, but that combinations of anions are not additive. For example, DCA exhibits a remarkable ability to form hydrogen bonds with the protein, resulting in a significantly stronger preferential binding to the protein than other anions. However, the combination of DCA with NO3, which also forms hydrogen bonds with the protein, results in a smaller preferential solvation of the protein than the combination of DCA with chloride ions, which are weaker binders. Thus, combining anions with varying affinities for the protein surface modulates the overall ion accumulation through nonadditive mechanisms, highlighting the importance of the understanding of competition for specific interaction sites, cooperative binding, bulk-solution affinity, and overall charge compensations, on the overall solvation capacity of the solution. Such knowledge may allow for the design of novel IL-based processes in biotechnology and material science, where fine-tuning protein solvation is crucial for optimizing performance and functionality.

Capturing the catalytic intermediates of parkin ubiquitination

Parkin is an E3 ubiquitin ligase implicated in early-onset forms of Parkinson's disease. It catalyzes a transthiolation reaction by accepting ubiquitin (Ub) from an E2 conjugating enzyme, forming a short-lived thioester intermediate, and transfers Ub to mitochondrial membrane substrates to signal mitophagy. A major impediment to the development of Parkinsonism therapeutics is the lack of structural and mechanistic detail for the essential, short-lived transthiolation intermediate. It is not known how Ub is recognized by the catalytic Rcat domain in parkin that enables Ub transfer from an E2~Ub conjugate to the catalytic site and the structure of the transthiolation complex is undetermined. Here, we capture the catalytic intermediate for the Rcat domain of parkin in complex with ubiquitin (Rcat-Ub) and determine its structure using NMR-based chemical shift perturbation experiments. We show that a previously unidentified α-helical region near the Rcat domain is unmasked as a recognition motif for Ub and guides the C-terminus of Ub toward the parkin catalytic site. Further, we apply a combination of guided AlphaFold modeling, chemical cross-linking, and single turnover assays to establish and validate a model of full-length parkin in complex with UbcH7, its donor Ub, and phosphoubiquitin, trapped in the process of transthiolation. Identification of this catalytic intermediate and orientation of Ub with respect to the Rcat domain provides important structural insights into Ub transfer by this E3 ligase and explains how the previously enigmatic Parkinson's pathogenic mutation T415N alters parkin activity.

A lanthanide tag for a complementary set of pseudocontact shifts

Pseudocontact shifts (PCS) generated by paramagnetic lanthanide ions deliver powerful restraints for protein structure analysis by NMR spectroscopy. We present a new lanthanide tag that generates different PCSs than that of a related tag, which differs in structure by a single oxygen atom. It is highly reactive towards cysteine and performs well in turn-on luminescence and in EPR spectroscopy.

Rescaling protein-protein interactions improves Martini 3 for flexible proteins in solution

Multidomain proteins with flexible linkers and disordered regions play important roles in many cellular processes, but characterizing their conformational ensembles is difficult. We have previously shown that the coarse-grained model, Martini 3, produces too compact ensembles in solution, that may in part be remedied by strengthening protein-water interactions. Here, we show that decreasing the strength of protein-protein interactions leads to improved agreement with experimental data on a wide set of systems. We show that the 'symmetry' between rescaling protein-water and protein-protein interactions breaks down when studying interactions with or within membranes; rescaling protein-protein interactions better preserves the binding specificity of proteins with lipid membranes, whereas rescaling protein-water interactions preserves oligomerization of transmembrane helices. We conclude that decreasing the strength of protein-protein interactions improves the accuracy of Martini 3 for IDPs and multidomain proteins, both in solution and in the presence of a lipid membrane.

Gadolinium-Based NMR Spin Relaxation Measurements of Near-Surface Electrostatic Potentials of Biomolecules

NMR spectroscopy is an important tool for the measurement of the electrostatic properties of biomolecules. To this point, paramagnetic relaxation enhancements (PREs) of 1H nuclei arising from nitroxide cosolutes in biomolecular solutions have been used to measure effective near-surface electrostatic potentials (ϕENS) of proteins and nucleic acids. Here, we present a gadolinium (Gd)-based NMR method, exploiting Gd chelates with different net charges, for measuring ϕENS values and demonstrate its utility through applications to a number of biomolecular systems. The use of Gd-based cosolutes offers several advantages over nitroxides for ϕENS measurements. First, unlike nitroxide compounds, Gd chelates enable electrostatic potential measurements on oxidation-sensitive proteins that require reducing agents. Second, the large electron spin quantum number of Gd (7/2) results in notably larger PREs for Gd chelates when used at the same concentrations as nitroxide radicals. Thus, it is possible to measure ϕENS values exclusively from + and - charged compounds even for highly charged biomolecules, avoiding the use of neutral cosolutes that, as we further establish here, limits the accuracy of the measured electrostatic potentials. In addition, the smaller concentrations of cosolutes required minimize potential binding to sites on macromolecules. Fourth, the closer proximity of the paramagnetic center and charged groups within Gd chelates, in comparison to the corresponding nitroxide compounds, enables more accurate predictions of ϕENS potentials for cross-validation of the experimental results. The Gd-based method described here, thus, broadens the applicability of studies of biomolecular electrostatics using solution NMR spectroscopy.

Ubiquitin-derived artificial binding proteins targeting oncofetal fibronectin reveal scaffold plasticity by β-strand slippage

Affilin proteins, artificial binding proteins based on the ubiquitin scaffold, have been generated by directed protein evolution to yield de-novo variants that bind the extra-domain B (EDB) of oncofetal fibronectin, an established marker of tumor neovasculature. The crystal structures of two EDB-specific Affilin variants reveal a striking structural plasticity of the ubiquitin scaffold, characterised by β-strand slippage, leading to different negative register shifts of the β5 strands. This process recruits amino acid residues from β5 towards the N-terminus to an adjacent loop region and subsequent residues into β5, respectively, remodeling the binding interface and leading to target specificity and affinity. Protein backbone alterations resulting from β-strand register shifts, as seen in the ubiquitin fold, can pose additional challenges to protein engineering as structural evidence of these events is still limited and they are difficult to predict. However, they can surface under the selection pressure of directed evolution and suggest that backbone plasticity allowing β-strand slippages can increase structural diversity, enhancing the evolutionary potential of a protein scaffold.

Conformation-Selective Ultraviolet-Ultraviolet Hole Burning Spectra of Ubiquitin Ions in a Cryogenic Ion Trap

Understanding the structural variations of conformational isomers in proteins is crucial for elucidating protein folding mechanisms. Here, we present a novel method for obtaining conformation-selective ultraviolet (UV)-UV hole burning (HB) spectra of ubiquitin ions ((Ubi+zH)+z, z = 7-10) produced via electrospray ionization. Our approach involves binding multiple N2 molecules to ubiquitin ions ((Ubi+zH)+z(N2)m, m = 1-55) within a cryogenic ion trap. Upon exposure to UV irradiation, efficient fragmentation of (Ubi+zH)+z(N2)m occurs, primarily yielding bare (Ubi+zH)+z ions as fragments. The significant mass difference between the parent and fragment ions facilitates the acquisition of UV-UV HB spectra, which reveal the presence of at least two distinct conformers. Molecular dynamics simulations suggest that these conformers correspond to A-state structures, differing only in the interactions of a tyrosine residue with neighboring residues. Our findings underscore UV-UV HB spectroscopy of protein ions as a powerful tool for exploring diverse protein isomers.

Unlocking the unfolded structure of ubiquitin: Combining time-resolved x-ray solution scattering and molecular dynamics to generate unfolded ensembles

The unfolding dynamics of ubiquitin were studied using a combination of x-ray solution scattering (XSS) and molecular dynamics (MD) simulations. The kinetic analysis of the XSS ubiquitin signals showed that the protein unfolds through a two-state process, independent of the presence of destabilizing salts. In order to characterize the ensemble of unfolded states in atomic detail, the experimental XSS results were used as a constraint in the MD simulations through the incorporation of x-ray scattering derived potential to drive the folded ubiquitin structure toward sampling unfolded states consistent with the XSS signals. We detail how biased MD simulations provide insight into unfolded states that are otherwise difficult to resolve and underscore how experimental XSS data can be combined with MD to efficiently sample structures away from the native state. Our results indicate that ubiquitin samples unfolded in states with a high degree of loss in secondary structure yet without a collapse to a molten globule or fully solvated extended chain. Finally, we propose how using biased-MD can significantly decrease the computational time and resources required to sample experimentally relevant nonequilibrium states.

Essential dynamics of ubiquitin in water and in a natural deep eutectic solvent

Natural deep eutectic solvents (NADESs) comprised of osmolytes are of interest as potential biomolecular (cryo)protectants. However, the way these solvents influence the structure and dynamics of biomolecules as well as the role of water remains poorly understood. We carried out principal component analysis of various secondary structure elements of ubiquitin in water and a betaine : glycerol : water (1 : 2 : ζ; ζ = 0, 1, 2, 5, 10, 20, 45) NADES, from molecular dynamics trajectories, to gain insight into the protein dynamics as it undergoes a transition from a highly viscous anhydrous to an aqueous environment. A crossover of the protein's essential dynamics at ζ ∼ 5, induced by solvent-shell coupled fluctuations, is observed, indicating that ubiquitin might (re)fold in the NADES upon water addition at ζ > ∼5. Further, in contrast to water, the anhydrous NADES preserves ubiquitin's essential modes at high temperatures explaining the protein's seemingly enhanced thermal stability.

Bacteria conjugate ubiquitin-like proteins to interfere with phage assembly

Several immune pathways in humans conjugate ubiquitin-like proteins to virus and host molecules as a means of antiviral defence1-5. Here we studied an antiphage defence system in bacteria, comprising a ubiquitin-like protein, ubiquitin-conjugating enzymes E1 and E2, and a deubiquitinase. We show that during phage infection, this system specifically conjugates the ubiquitin-like protein to the phage central tail fibre, a protein at the tip of the tail that is essential for tail assembly as well as for recognition of the target host receptor. Following infection, cells encoding this defence system release a mixture of partially assembled, tailless phage particles and fully assembled phages in which the central tail fibre is obstructed by the covalently attached ubiquitin-like protein. These phages show severely impaired infectivity, explaining how the defence system protects the bacterial population from the spread of phage infection. Our findings demonstrate that conjugation of ubiquitin-like proteins is an antiviral strategy conserved across the tree of life.

Human apo-metallothionein 1a is not a random coil: Evidence from guanidinium chloride, high temperature, and acidic pH unfolding studies

The structures of apo-metallothioneins (apo-MTs) have been relatively elusive due to their fluxional, disordered state which has been difficult to characterize. However, intrinsically disordered protein (IDP) structures are rather diverse, which raises questions about where the structure of apo-MTs fit into the protein structural spectrum. In this paper, the unfolding transitions of apo-MT1a are discussed with respect to the effect of the chemical denaturant GdmCl, temperature conditions, and pH environment. Cysteine modification in combination with electrospray ionization mass spectrometry was used to probe the unfolding transition of apo-MT1a in terms of cysteine exposure. Circular dichroism spectroscopy was also used to monitor the change in secondary structure as a function of GdmCl concentration. For both of these techniques, cooperative unfolding was observed, suggesting that apo-MT1a is not a random coil. More GdmCl was required to unfold the protein backbone than to expose the cysteines, indicating that cysteine exposure is likely an early step in the unfolding of apo-MT1a. MD simulations complement the experimental results, suggesting that apo-MT1a adopts a more compact structure than expected for a random coil. Overall, these results provide further insight into the intrinsically disordered structure of apo-MT.

Bacterial esterases reverse lipopolysaccharide ubiquitylation to block host immunity

Aspects of how Burkholderia escape the host's intrinsic immune response to replicate in the cell cytosol remain enigmatic. Here, we show that Burkholderia has evolved two mechanisms to block the activity of Ring finger protein 213 (RNF213)-mediated non-canonical ubiquitylation of bacterial lipopolysaccharide (LPS), thereby preventing the initiation of antibacterial autophagy. First, Burkholderia's polysaccharide capsule blocks RNF213 association with bacteria and second, the Burkholderia deubiquitylase (DUB), TssM, directly reverses the activity of RNF213 through a previously unrecognized esterase activity. Structural analysis provides insight into the molecular basis of TssM esterase activity, allowing it to be uncoupled from its isopeptidase function. Furthermore, a putative TssM homolog also displays esterase activity and removes ubiquitin from LPS, establishing this as a virulence mechanism. Of note, we also find that additional immune-evasion mechanisms exist, revealing that overcoming this arm of the host's immune response is critical to the pathogen.

Proline Peptide Bond Isomerization in Ubiquitin Under Folding and Denaturing Conditions by Pressure-Jump NMR

Proline isomerization is widely recognized as a kinetic bottleneck in protein folding, amplified for proteins rich in Pro residues. We introduced repeated hydrostatic pressure jumps between native and pressure-denaturing conditions inside an NMR sample cell to study proline isomerization in the pressure-sensitized L50A ubiquitin mutant. Whereas in two unfolded heptapeptides, X-Pro peptide bonds isomerized ca 1.6-fold faster at 1 bar than at 2.5 kbar, for ubiquitin ca eight-fold faster isomerization was observed for Pro-38 and ca two-fold for Pro-19 and Pro-37 relative to rates measured in the pressure-denatured state. Activation energies for isomerization in pressure-denatured ubiquitin were close to literature values of 20 kcal/mole for denatured polypeptides but showed a substantial drop to 12.7 kcal/mole for Pro-38 at atmospheric pressure. For ubiquitin isomers with a cis E18-P19 peptide bond, the 1-bar NMR spectrum showed sharp resonances with near random coil chemical shifts for the C-terminal half of the protein, characteristic of an unfolded chain, while most of the N-terminal residues were invisible due to exchange broadening, pointing to a metastable partially folded state for this previously recognized 'folding nucleus'. For cis-P37 isomers, a drop in pressure resulted in the rapid loss of nearly all unfolded-state NMR resonances, while the recovery of native state intensity revealed a slow component attributed to cis → trans isomerization of P37. This result implies that the NMR-invisible cis-P37 isomer adopts a molten globule state that encompasses the entire length of the ubiquitin chain, suggestive of a structure that mostly resembles the folded state.

Map conformational landscapes of intrinsically disordered proteins with polymer physics quantities

Disordered proteins are conformationally flexible proteins that are biologically important and have been implicated in devastating diseases such as Alzheimer's disease and cancer. Unlike stably folded structured proteins, disordered proteins sample a range of different conformations that needs to be accounted for. Here, we treat disordered proteins as polymer chains, and compute a dimensionless quantity called instantaneous shape ratio (Rs), as Rs = Ree2/Rg2, where Ree is end-to-end distance and Rg is radius of gyration. Extended protein conformations tend to have high Ree compared with Rg, and thus have high Rs values, whereas compact conformations have smaller Rs values. We use a scatter plot of Rs (representing shape) against Rg (representing size) as a simple map of conformational landscapes. We first examine the conformational landscape of simple polymer models such as Random Walk, Self-Avoiding Walk, and Gaussian Walk (GW), and we notice that all protein/polymer maps lie within the boundaries of the GW map. We thus use the GW map as a reference and, to assess conformational diversity, we compute the fraction of the GW conformations (fC) covered by each protein/polymer. Disordered proteins all have high fC scores, consistent with their disordered nature. Each disordered protein accesses a different region of the reference map, revealing differences in their conformational ensembles. We additionally examine the conformational maps of the nonviral gene delivery vector polyethyleneimine at various protonation states, and find that they resemble disordered proteins, with coverage of the reference map decreasing with increasing protonation state, indicating decreasing conformational diversity. We propose that our method of combining Rs and Rg in a scatter plot generates a simple, meaningful map of the conformational landscape of a disordered protein, which in turn can be used to assess conformational diversity of disordered proteins.

Missing Wedge Completion via Unsupervised Learning with Coordinate Networks

Cryogenic electron tomography (cryoET) is a powerful tool in structural biology, enabling detailed 3D imaging of biological specimens at a resolution of nanometers. Despite its potential, cryoET faces challenges such as the missing wedge problem, which limits reconstruction quality due to incomplete data collection angles. Recently, supervised deep learning methods leveraging convolutional neural networks (CNNs) have considerably addressed this issue; however, their pretraining requirements render them susceptible to inaccuracies and artifacts, particularly when representative training data is scarce. To overcome these limitations, we introduce a proof-of-concept unsupervised learning approach using coordinate networks (CNs) that optimizes network weights directly against input projections. This eliminates the need for pretraining, reducing reconstruction runtime by 3-20× compared to supervised methods. Our in silico results show improved shape completion and reduction of missing wedge artifacts, assessed through several voxel-based image quality metrics in real space and a novel directional Fourier Shell Correlation (FSC) metric. Our study illuminates benefits and considerations of both supervised and unsupervised approaches, guiding the development of improved reconstruction strategies.

Not All Binding Sites Are Equal: Site Determination and Folding State Analysis of Gas-Phase Protein-Metallodrug Adducts

Modern approaches in metallodrug research focus on compounds that bind protein targets rather than DNA. However, the identification of protein targets and binding sites is challenging. Using intact mass spectrometry and proteomics, we investigated the binding of the antimetastatic agent RAPTA-C to the model proteins ubiquitin, cytochrome c, lysozyme, and myoglobin. Binding to cytochrome c and lysozyme was negligible. However, ubiquitin bound up to three Ru moieties, two of which were localized at Met1 and His68 as [Ru(cym)], and [Ru(cym)] or [Ru(cym)(PTA)] adducts, respectively. Myoglobin bound up to four [Ru(cym)(PTA)] moieties and five sites were identified at His24, His36, His64, His81/82 and His113. Collision-induced unfolding (CIU) studies via ion-mobility mass spectrometry allowed measuring protein folding as a function of collisional activation. CIU of protein-RAPTA-C adducts showed binding of [Ru(cym)] to Met1 caused a significant compaction of ubiquitin, likely from N-terminal S-Ru-N chelation, while binding of [Ru(cym)(PTA)] to His residues of ubiquitin or myoglobin induced a smaller effect. Interestingly, the folded state of ubiquitin formed by His functionalization was more stable than Met1 metalation. The data suggests that selective metalation of amino acids at different positions on the protein impacts the conformation and potentially the biological activity of anticancer compounds.

Proteins Can Withstand More Extensive Labeling while Providing Accurate Structural Information in Covalent Labeling-Mass Spectrometry

Diethylpyrocarbonate (DEPC) covalent labeling-mass spectrometry (CL-MS) has been extensively utilized to study protein structure and interactions owing to its ease of use, commercial availability, and broad labeling of nucleophilic residues. During typical CL-MS experiments with DEPC, the extent of labeling is kept low to avoid any structural perturbations resulting from covalent modification of the protein. In this study, we demonstrate that proteins can be labeled more extensively via DEPC and still provide accurate structural information. To show this, we modeled labeling kinetics over a range of DEPC concentrations and used molecular dynamics simulations to investigate the molecular-level effects of extensive labeling on the protein structure. Our results indicate that higher extents of DEPC labeling do not significantly perturb the protein structure and can lead to improved precision, detectability of labeled peptides, and protein structural resolution. Furthermore, higher extents of labeling enable better identification of protein-ligand binding sites where lower extents of modification provide ambiguous results.

Missing Wedge Completion via Unsupervised Learning with Coordinate Networks

Cryogenic electron tomography (cryoET) is a powerful tool in structural biology, enabling detailed 3D imaging of biological specimens at a resolution of nanometers. Despite its potential, cryoET faces challenges such as the missing wedge problem, which limits reconstruction quality due to incomplete data collection angles. Recently, supervised deep learning methods leveraging convolutional neural networks (CNNs) have considerably addressed this issue; however, their pretraining requirements render them susceptible to inaccuracies and artifacts, particularly when representative training data is scarce. To overcome these limitations, we introduce a proof-of-concept unsupervised learning approach using coordinate networks (CNs) that optimizes network weights directly against input projections. This eliminates the need for pretraining, reducing reconstruction runtime by 3 - 20× compared to supervised methods. Our in silico results show improved shape completion and reduction of missing wedge artifacts, assessed through several voxel-based image quality metrics in real space and a novel directional Fourier Shell Correlation (FSC) metric. Our study illuminates benefits and considerations of both supervised and unsupervised approaches, guiding the development of improved reconstruction strategies.

Molecular dynamics simulation of the effect of temperature on the conformation of ubiquitin protein

CONTENT

Ubiquitin, a ubiquitous small protein found in all living organisms, is crucial for tagging proteins earmarked for degradation and holds pivotal importance in biomedicine. Protein functionality is intricately linked to its structure. To comprehend the impact of diverse temperatures on ubiquitin protein structure, our study delved into the energy landscape, hydrogen bonding, and overall structural stability of ubiquitin protein at varying temperatures. Through meticulous analysis of root mean square deviation and root mean square fluctuation, we validated the robustness of the simulation conditions employed. Within our simulated system, the bonding energy and electrostatic potential energy exhibited linear augmentation, while the van der Waals energy demonstrated a linear decline. Additionally, our findings highlighted that the α-Helix secondary structure of the ubiquitin protein gradually transitions toward helix destabilization under high-temperature conditions. The secondary structure of ubiquitin protein experiences distinct changes under varying temperatures. The outcomes of our molecular simulations offer a theoretical framework that enhances our comprehension of how temperature impacts the structural stability of ubiquitin protein. These insights contribute not only to a deeper understanding of iniquity's behavior but also hold broader implications in the realm of biomedicine and beyond.

METHODS

All the MD simulations were performed using the GROMACS software with GROMOS96 force field and SPC for water. The ubiquitin protein was put in the center of a cubic box with a length of 8 nm, a setting that allowed > 0.8 nm in the minimal distance between the protein surface and the box wall. To remove the possible coordinate collision of the configurations, in the beginning, the steepest descent method was used until the maximum force between atoms was under 100 kJ/mol·nm with a 0.01 nm step size. Minimization was followed by 30 ps of position-restrained MD simulation. The protein was restrained to its initial position, and the solvent was freely equilibrated. The product phase was obtained with the whole system simulated for 10 ns without any restraint using an integral time step of 1 fs with different temperatures. The cutoff for short-range electronic interaction was set to 1.5 nm. The long-range interactions were treated with a particle-mesh Ewald (PME) method with a grid width of 1.2 nm.

Ubiquitylation-independent cotranslational degradation of dihydrofolate reductase and ubiquitin

Nascent proteins are degraded during or immediately after synthesis, a process called cotranslational protein degradation (CTPD). Although CTPD was observed decades ago, it has never been fully explored mechanistically and functionally. We show here that dihydrofolate reductase (DHFR) and ubiquitin (Ub), two stable proteins widely used in protein degradation studies, are actually subject to CTPD. Unlike canonical posttranslational protein degradation, CTPD of DHFR and Ub does not require prior ubiquitylation. Our data also suggest that protein expression level and N-terminal folding pattern may be two critical determinants for CTPD. Thus, this study reveals that CTPD plays a role in regulating the homeostasis of long-lived proteins and provides insights into the mechanism of CTPD.

Decoding Chromatin Ubiquitylation: A Chemical Biology Perspective

Since Strahl and Allis proposed the "language of covalent histone modifications", a host of experimental studies have shed light on the different facets of chromatin regulation by epigenetic mechanisms. Initially proposed as a concept for controlling gene transcription, the regulation of deposition and removal of histone post-translational modifications (PTMs), such as acetylation, methylation, and phosphorylation, have been implicated in many chromatin regulation pathways. However, large PTMs such as ubiquitylation challenge research on many levels due to their chemical complexity. In recent years, chemical tools have been developed to generate chromatin in defined ubiquitylation states in vitro. Chemical biology approaches are now used to link specific histone ubiquitylation marks with downstream chromatin regulation events on the molecular level. Here, we want to highlight how chemical biology approaches have empowered the mechanistic study of chromatin ubiquitylation in the context of gene regulation and DNA repair with attention to future challenges.

Fragment-Based Deep Learning for Simultaneous Prediction of Polarizabilities and NMR Shieldings of Macromolecules and Their Aggregates

Simultaneous prediction of the molecular response properties, such as polarizability and the NMR shielding constant, at a low computational cost is an unresolved issue. We propose to combine a linear-scaling generalized energy-based fragmentation (GEBF) method and deep learning (DL) with both molecular and atomic information-theoretic approach (ITA) quantities as effective descriptors. In GEBF, the total molecular polarizability can be assembled as a linear combination of the corresponding quantities calculated from a set of small embedded subsystems in GEBF. In the new GEBF-DL(ITA) protocol, one can predict subsystem polarizabilities based on the corresponding molecular wave function (thus electron density and ITA quantities) and DL model rather than calculate them from the computationally intensive coupled-perturbed Hartree-Fock or Kohn-Sham equations and finally obtain the total molecular polarizability via a linear combination equation. As a proof-of-concept application, we predict the molecular polarizabilities of large proteins and protein aggregates. GEBF-DL(ITA) is shown to be as accurate enough as GEBF, with mean absolute percentage error <1%. For the largest protein aggregate (>4000 atoms), GEBF-DL(ITA) gains a speedup ratio of 3 compared with GEBF. It is anticipated that when more advanced electronic structure methods are used, this advantage will be more appealing. Moreover, one can also predict the NMR chemical shieldings of proteins with reasonably good accuracy. Overall, the cost-efficient GEBF-DL(ITA) protocol should be a robust theoretical tool for simultaneously predicting polarizabilities and NMR shieldings of large systems.

Molecular Crowding Alters the Interactions of Polymyxin Lipopeptides within the Periplasm of E. coli: Insights from Molecular Dynamics

The cell envelope of Gram-negative bacteria is a crowded tripartite architecture that separates the cell interior from the external environment. Two membranes encapsulate the aqueous periplasm, which contains the cell wall. Little is known about the mechanisms via which antimicrobial peptides move through the periplasm from the outer membrane to their site of action, the inner membrane. We utilize all-atom molecular dynamics to study two antimicrobial peptides, polymyxins B1 and E, within models of the E. coli periplasm crowded to different extents. In a simple chemical environment, both PMB1 and PME bind irreversibly to the cell wall. The presence of specific macromolecules leads to competition with the polymyxins for cell wall interaction sites, resulting in polymyxin dissociation from the cell wall. Chemical complexity also impacts interactions between polymyxins and Braun's lipoprotein; thus, the interaction modes of lipoprotein antibiotics within the periplasm are dependent upon the nature of the other species present.

Proteotoxic stress and the ubiquitin proteasome system

The ubiquitin proteasome system maintains protein homeostasis by regulating the breakdown of misfolded proteins, thereby preventing misfolded protein aggregates. The efficient elimination is vital for preventing damage to the cell by misfolded proteins, known as proteotoxic stress. Proteotoxic stress can lead to the collapse of protein homeostasis and can alter the function of the ubiquitin proteasome system. Conversely, impairment of the ubiquitin proteasome system can also cause proteotoxic stress and disrupt protein homeostasis. This review examines two impacts of proteotoxic stress, 1) disruptions to ubiquitin homeostasis (ubiquitin stress) and 2) disruptions to proteasome homeostasis (proteasome stress). Here, we provide a mechanistic description of the relationship between proteotoxic stress and the ubiquitin proteasome system. This relationship is illustrated by findings from several protein misfolding diseases, mainly neurodegenerative diseases, as well as from basic biology discoveries from yeast to mammals. In addition, we explore the importance of the ubiquitin proteasome system in endoplasmic reticulum quality control, and how proteotoxic stress at this organelle is alleviated. Finally, we highlight how cells utilize the ubiquitin proteasome system to adapt to proteotoxic stress and how the ubiquitin proteasome system can be genetically and pharmacologically manipulated to maintain protein homeostasis.

Unraveling motion in proteins by combining NMR relaxometry and molecular dynamics simulations: A case study on ubiquitin

Nuclear magnetic resonance (NMR) relaxation experiments shine light onto the dynamics of molecular systems in the picosecond to millisecond timescales. As these methods cannot provide an atomically resolved view of the motion of atoms, functional groups, or domains giving rise to such signals, relaxation techniques have been combined with molecular dynamics (MD) simulations to obtain mechanistic descriptions and gain insights into the functional role of side chain or domain motion. In this work, we present a comparison of five computational methods that permit the joint analysis of MD simulations and NMR relaxation experiments. We discuss their relative strengths and areas of applicability and demonstrate how they may be utilized to interpret the dynamics in MD simulations with the small protein ubiquitin as a test system. We focus on the aliphatic side chains given the rigidity of the backbone of this protein. We find encouraging agreement between experiment, Markov state models built in the χ1/χ2 rotamer space of isoleucine residues, explicit rotamer jump models, and a decomposition of the motion using ROMANCE. These methods allow us to ascribe the dynamics to specific rotamer jumps. Simulations with eight different combinations of force field and water model highlight how the different metrics may be employed to pinpoint force field deficiencies. Furthermore, the presented comparison offers a perspective on the utility of NMR relaxation to serve as validation data for the prediction of kinetics by state-of-the-art biomolecular force fields.

Molecular Dynamics Simulations of Native Protein Charging via Proton Transfer during Electrospray Ionization with Grotthuss Diffuse H3O. Anal Chem 2024;96:4146-4153. [PMID: 38427846 PMCID: PMC11337394 DOI: 10.1021/acs.analchem.3c05089] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

Unraveling the mechanism by which native proteins are charged through electrospray ionization (ESI) has been the focus of considerable research because observable charge states can be correlated to biophysical characteristics, such as protein folding and, thus, solution conformation. Difficulties in characterizing electrosprayed droplets have catalyzed the use of molecular dynamics (MD) to provide insights into the mechanisms by which proteins are charged and transferred to the gas phase. However, prior MD studies have utilized metal ions, primarily Na+, as charge carriers, even though proteins are primarily detected as protonated ions in the mass spectra. Here, we propose a modified MD protocol for simulating discrete Grotthuss diffuse H3O+ that is capable of dynamically altering amino-acid protonation states to model electrospray charging and gaseous ion formation of model proteins, ubiquitin, and myoglobin. Application of the protocol to the evaporation of acidic droplets enables a molecular perspective of H3O+ coordination and proton transfer to/from proteins, which is unfeasible with the metal charge carriers used in previous MD studies of ESI. Our protocol recreates experimentally observed charge-state distributions and supports the charge residue model (CRM) as the dominant mechanism of native protein ionization during ESI. Additionally, our results suggest that protonation is highly specific to individual residues and is correlated to the formation of localized hydrated regions on the protein surface as droplets desolvate. Considering the use of discrete H3O+ instead of Na+, the developed protocol is a necessary step toward developing a more comprehensive model of protein ionization during ESI.

Extending the Range of Distances Accessible by 19F Electron-Nuclear Double Resonance in Proteins Using High-Spin Gd(III) Labels

Fluorine electron-nuclear double resonance (19F ENDOR) has recently emerged as a valuable tool in structural biology for distance determination between F atoms and a paramagnetic center, either intrinsic or conjugated to a biomolecule via spin labeling. Such measurements allow access to distances too short to be measured by double electron-electron resonance (DEER). To further extend the accessible distance range, we exploit the high-spin properties of Gd(III) and focus on transitions other than the central transition (|-1/2⟩ ↔ |+1/2⟩), that become more populated at high magnetic fields and low temperatures. This increases the spectral resolution up to ca. 7 times, thus raising the long-distance limit of 19F ENDOR almost 2-fold. We first demonstrate this on a model fluorine-containing Gd(III) complex with a well-resolved 19F spectrum in conventional central transition measurements and show quantitative agreement between the experimental spectra and theoretical predictions. We then validate our approach on two proteins labeled with 19F and Gd(III), in which the Gd-F distance is too long to produce a well-resolved 19F ENDOR doublet when measured at the central transition. By focusing on the |-5/2⟩ ↔ |-3/2⟩ and |-7/2⟩ ↔ |-5/2⟩ EPR transitions, a resolution enhancement of 4.5- and 7-fold was obtained, respectively. We also present data analysis strategies to handle contributions of different electron spin manifolds to the ENDOR spectrum. Our new extended 19F ENDOR approach may be applicable to Gd-F distances as large as 20 Å, widening the current ENDOR distance window.

The ubiquitin codes in cellular stress responses

Ubiquitination/ubiquitylation, one of the most fundamental post-translational modifications, regulates almost every critical cellular process in eukaryotes. Emerging evidence has shown that essential components of numerous biological processes undergo ubiquitination in mammalian cells upon exposure to diverse stresses, from exogenous factors to cellular reactions, causing a dazzling variety of functional consequences. Various forms of ubiquitin signals generated by ubiquitylation events in specific milieus, known as ubiquitin codes, constitute an intrinsic part of myriad cellular stress responses. These ubiquitination events, leading to proteolytic turnover of the substrates or just switch in functionality, initiate, regulate, or supervise multiple cellular stress-associated responses, supporting adaptation, homeostasis recovery, and survival of the stressed cells. In this review, we attempted to summarize the crucial roles of ubiquitination in response to different environmental and intracellular stresses, while discussing how stresses modulate the ubiquitin system. This review also updates the most recent advances in understanding ubiquitination machinery as well as different stress responses and discusses some important questions that may warrant future investigation.

Improved Protein Dynamics and Hydration in the Martini3 Coarse-Grain Model

The Martini coarse-grain force-field has emerged as an important framework to probe cellular processes at experimentally relevant time- and length-scales. However, the recently developed version, the Martini3 force-field with the implemented Go̅ model (Martini3Go̅), as well as previous variants of the Martini model have not been benchmarked and rigorously tested for globular proteins. In this study, we consider three globular proteins, ubiquitin, lysozyme, and cofilin, and compare protein dynamics and hydration with observables from experiments and all-atom simulations. We show that the Martini3Go̅ model is able to accurately model the structural and dynamic features of small globular proteins. Overall, the structural integrity of the proteins is maintained, as validated by contact maps, radii of gyration (Rg), and SAXS profiles. The chemical shifts predicted from the ensemble sampled in the simulations are consistent with the experimental data. Further, a good match is observed in the protein-water interaction energetics, and the hydration levels of the residues are similar to atomistic simulations. However, the protein-water interaction dynamics is not accurately represented and appears to depend on the protein structural complexity, residue specificity, and water dynamics. Our work is a step toward testing and assessing the Martini3Go̅ model and provides insights into future efforts to refine Martini models with improved solvation effects and better correspondence to the underlying all-atom systems.

Surveying the energy landscape of coarse-grained mappings

Simulations of soft materials often adopt low-resolution coarse-grained (CG) models. However, the CG representation is not unique and its impact upon simulated properties is poorly understood. In this work, we investigate the space of CG representations for ubiquitin, which is a typical globular protein with 72 amino acids. We employ Monte Carlo methods to ergodically sample this space and to characterize its landscape. By adopting the Gaussian network model as an analytically tractable atomistic model for equilibrium fluctuations, we exactly assess the intrinsic quality of each CG representation without introducing any approximations in sampling configurations or in modeling interactions. We focus on two metrics, the spectral quality and the information content, that quantify the extent to which the CG representation preserves low-frequency, large-amplitude motions and configurational information, respectively. The spectral quality and information content are weakly correlated among high-resolution representations but become strongly anticorrelated among low-resolution representations. Representations with maximal spectral quality appear consistent with physical intuition, while low-resolution representations with maximal information content do not. Interestingly, quenching studies indicate that the energy landscape of mapping space is very smooth and highly connected. Moreover, our study suggests a critical resolution below which a "phase transition" qualitatively distinguishes good and bad representations.

The structural basis for deubiquitination by the fingerless USP-type effector TssM

Intracellular bacteria are threatened by ubiquitin-mediated autophagy, whenever the bacterial surface or enclosing membrane structures become targets of host ubiquitin ligases. As a countermeasure, many intracellular pathogens encode deubiquitinase (DUB) effectors to keep their surfaces free of ubiquitin. Most bacterial DUBs belong to the OTU or CE-clan families. The betaproteobacteria Burkholderia pseudomallei and Burkholderia mallei, causative agents of melioidosis and glanders, respectively, encode the TssM effector, the only known bacterial DUB belonging to the USP class. TssM is much shorter than typical eukaryotic USP enzymes and lacks the canonical ubiquitin-recognition region. By solving the crystal structures of isolated TssM and its complex with ubiquitin, we found that TssM lacks the entire "Fingers" subdomain of the USP fold. Instead, the TssM family has evolved the functionally analog "Littlefinger" loop, which is located towards the end of the USP domain and recognizes different ubiquitin interfaces than those used by USPs. The structures revealed the presence of an N-terminal immunoglobulin-fold domain, which is able to form a strand-exchange dimer and might mediate TssM localization to the bacterial surface.

Almost 50 Years of Monomeric Extracellular Ubiquitin (eUb)

Monomeric ubiquitin (Ub) is a 76-amino-acid highly conserved protein found in eukaryotes. The biological activity of Ub first described in the 1970s was extracellular, but it quickly gained relevance due to its intracellular role, i.e., post-translational modification of intracellular proteins (ubiquitination) that regulate numerous eukaryotic cellular processes. In the following years, the extracellular role of Ub was relegated to the background, until a correlation between higher survival rate and increased serum Ub concentrations in patients with sepsis and burns was observed. Although the mechanism of action (MoA) of extracellular ubiquitin (eUb) is not yet well understood, further studies have shown that it may ameliorate the inflammatory response in tissue injury and multiple sclerosis diseases. These observations, compounded with the high stability and low immunogenicity of eUb due to its high conservation in eukaryotes, have made this small protein a relevant candidate for biotherapeutic development. Here, we review the in vitro and in vivo effects of eUb on immunologic, cardiovascular, and nervous systems, and discuss the potential MoAs of eUb as an anti-inflammatory, antimicrobial, and cardio- and brain-protective agent.

AMBERff at Scale: Multimillion-Atom Simulations with AMBER Force Fields in NAMD

All-atom molecular dynamics (MD) simulations are an essential structural biology technique with increasing application to multimillion-atom systems, including viruses and cellular machinery. Classical MD simulations rely on parameter sets, such as the AMBER family of force fields (AMBERff), to accurately describe molecular motion. Here, we present an implementation of AMBERff for use in NAMD that overcomes previous limitations to enable high-performance, massively parallel simulations encompassing up to two billion atoms. Single-point potential energy comparisons and case studies on model systems demonstrate that the implementation produces results that are as accurate as running AMBERff in its native engine.

Using NMR diffusion data to validate MD models of disordered proteins: Test case of N-terminal tail of histone H4

MD simulations can provide uniquely detailed models of intrinsically disordered proteins (IDPs). However, these models need careful experimental validation. The coefficient of translational diffusion Dtr, measurable by pulsed field gradient NMR, offers a potentially useful piece of experimental information related to the compactness of the IDP's conformational ensemble. Here, we investigate, both experimentally and via the MD modeling, the translational diffusion of a 25-residue N-terminal fragment from histone H4 (N-H4). We found that the predicted values of Dtr, as obtained from mean-square displacement of the peptide in the MD simulations, are largely determined by the viscosity of the MD water (which has been reinvestigated as a part of our study). Beyond that, our analysis of the diffusion data indicates that MD simulations of N-H4 in the TIP4P-Ew water give rise to an overly compact conformational ensemble for this peptide. In contrast, TIP4P-D and OPC simulations produce the ensembles that are consistent with the experimental Dtr result. These observations are supported by the analyses of the 15N spin relaxation rates. We also tested a number of empirical methods to predict Dtr based on IDP's coordinates extracted from the MD snapshots. In particular, we show that the popular approach involving the program HYDROPRO can produce misleading results. This happens because HYDROPRO is not intended to predict the diffusion properties of highly flexible biopolymers such as IDPs. Likewise, recent empirical schemes that exploit the relationship between the small-angle x-ray scattering-informed conformational ensembles of IDPs and the respective experimental Dtr values also prove to be problematic. In this sense, the first-principle calculations of Dtr from the MD simulations, such as demonstrated in this work, should provide a useful benchmark for future efforts in this area.

Protein stability in a natural deep eutectic solvent: Preferential hydration or solvent slaving?

Deep eutectic solvents (DESs) emerged as potential alternative solvent media in multiple areas, including biomolecular (cryo)preservation. Herein, we studied the stability of a small protein (ubiquitin) in water and a betaine-glycerol-water (B:G:W) (1:2:ζ; ζ = 0, 1, 2, 5, 10) DES, through molecular dynamics. An AMBER-based model that accurately describes the density and shear viscosity of the DES is proposed. We find that water molecules are largely trapped in the solvent, precluding the formation of a full hydration layer, seemingly opposite to osmolytes' preferential exclusion/preferential hydration mechanism. Although the protein is stable in the DES, structural fluctuations are largely suppressed and only recovered upon sufficient hydration. This is explained by a solvent-slaving mechanism where β-fluctuations are key, with the non-monotonic hydration of some amino acids with the water content providing an explanation to the non-monotonic folding of some proteins in aqueous DESs. A major thermal stability enhancement in the DES is also observed, caused by a similar slowdown of the backbone torsional dynamics. Our results support a kinetic stabilization of the protein in the DES, whereas a possible thermodynamic stabilization does not follow a preferential hydration or water entrapment mechanism.

Bacterial ligases reveal fundamental principles of polyubiquitin specificity

Homologous to E6AP C terminus (HECT) E3 ubiquitin (Ub) ligases direct substrates toward distinct cellular fates dictated by the specific form of monomeric or polymeric Ub (polyUb) signal attached. How polyUb specificity is achieved has been a long-standing mystery, despite extensive study in various hosts, ranging from yeast to human. The bacterial pathogens enterohemorrhagic Escherichia coli and Salmonella Typhimurium encode outlying examples of "HECT-like" (bHECT) E3 ligases, but commonalities to eukaryotic HECT (eHECT) mechanism and specificity had not been explored. We expanded the bHECT family with examples in human and plant pathogens. Three bHECT structures in primed, Ub-loaded states resolved key details of the entire Ub ligation process. One structure provided a rare glimpse into the act of ligating polyUb, yielding a means to rewire polyUb specificity of both bHECT and eHECT ligases. Studying this evolutionarily distinct bHECT family has revealed insight into the function of key bacterial virulence factors as well as fundamental principles underlying HECT-type Ub ligation.

Parkin-mediated ubiquitination inhibits BAK apoptotic activity by blocking its canonical hydrophobic groove

BAK permeabilizes the mitochondrial outer membrane, causing apoptosis. This apoptotic activity of BAK is stimulated by binding prodeath activators within its canonical hydrophobic groove. Parkin, an E3 ubiquitin (Ub) ligase, can ubiquitinate BAK, which inhibits BAK apoptotic activity. However, the molecular mechanism underlying the inhibition of ubiquitination remains structurally uncharacterized. Here, we utilize truncated and soluble BAK to construct a mimetic of K113-ubiquitinated BAK (disulfide-linked UbG76C ~ BAKK113C) and further present its NMR-derived structure model. The classical L8-I44-H68-V70 hydrophobic patch of the conjugated Ub subunit binds within the canonical hydrophobic groove of BAK. This Ub occludes the binding of prodeath BID activators in the groove and impairs BID-triggered BAK activation and membrane permeabilization. Reduced interaction between Ub and BAK subunits allows BID to activate K113-ubiquitinated BAK. These mechanistic insights suggest a nonsignaling function of Ub in that it directly antagonizes stimuli targeting Ub-modified proteins rather than by recruiting downstream partners for cellular messaging.

When MELD Meets GaMD: Accelerating Biomolecular Landscape Exploration

We introduce Gaussian accelerated MELD (GaMELD) as a new method for exploring the energy landscape of biomolecules. GaMELD combines the strengths of Gaussian accelerated molecular dynamics (GaMD) and modeling employing limited data (MELD) to navigate complex energy landscapes. MELD uses replica-exchange molecular simulations to integrate limited and uncertain data into simulations via Bayesian inference. MELD has been successfully applied to problems of structure prediction like protein folding and complex structure prediction. However, the computational cost for MELD simulations has limited its broader applicability. The synergy of GaMD and MELD surmounts this limitation efficiently sampling the energy landscape at a lower computational cost (reducing the computational cost by a factor of 2 to six). Effectively, GaMD is used to shift energy distributions along replicas to increase the overlap in energy distributions across replicas, facilitating a random walk in replica space. We tested GaMELD on a benchmark set of 12 small proteins that have been previously studied through MELD and conventional MD. GaMELD consistently achieves accurate predictions with fewer replicas. By increasing the efficacy of replica exchange, GaMELD effectively accelerates convergence in the conformational space, enabling improved sampling across a diverse set of systems.

Ubiquitin ligase enzymes and de-ubiquitinating enzymes regulate innate immunity in the TLR, NLR, RLR, and cGAS-STING pathways

Ubiquitination (or ubiquitylation) and de-ubiquitination, which are both post-translational modifications (PTMs) of proteins, have become a research hotspot in recent years. Some ubiquitinated or de-ubiquitinated signaling proteins have been found to promote or suppress innate immunity through Toll-like receptor (TLR), RIG-like receptor (RIG-I-like receptor, RLR), NOD-like receptor (NLR), and the cyclic guanosine monophosphate (GMP)-adenosine monophosphate (AMP) synthase (cGAS)-STING pathway. This article aimed to provide a review on the role of ubiquitination and de-ubiquitination, especially ubiquitin ligase enzymes and de-ubiquitinating enzymes, in the above four pathways. We hope that our work can contribute to the research and development of treatment strategies for innate immunity-related diseases such as inflammatory bowel disease.

Deubiquitylating Enzymes in Cancer and Immunity

Deubiquitylating enzymes (DUBs) maintain relative homeostasis of the cellular ubiquitome by removing the post-translational modification ubiquitin moiety from substrates. Numerous DUBs have been demonstrated specificity for cleaving a certain type of ubiquitin linkage or positions within ubiquitin chains. Moreover, several DUBs perform functions through specific protein-protein interactions in a catalytically independent manner, which further expands the versatility and complexity of DUBs' functions. Dysregulation of DUBs disrupts the dynamic equilibrium of ubiquitome and causes various diseases, especially cancer and immune disorders. This review summarizes the Janus-faced roles of DUBs in cancer including proteasomal degradation, DNA repair, apoptosis, and tumor metastasis, as well as in immunity involving innate immune receptor signaling and inflammatory and autoimmune disorders. The prospects and challenges for the clinical development of DUB inhibitors are further discussed. The review provides a comprehensive understanding of the multi-faced roles of DUBs in cancer and immunity.

Mechanism of selective recognition of Lys48-linked polyubiquitin by macrocyclic peptide inhibitors of proteasomal degradation

Post-translational modification of proteins with polyubiquitin chains is a critical cellular signaling mechanism in eukaryotes with implications in various cellular states and processes. Unregulated ubiquitin-mediated protein degradation can be detrimental to cellular homeostasis, causing numerous diseases including cancers. Recently, macrocyclic peptides were developed that selectively target long Lysine-48-linked polyubiquitin chains (tetra-ubiquitin) to inhibit ubiquitin-proteasome system, leading to attenuation of tumor growth in vivo. However, structural determinants of the chain length and linkage selectivity by these cyclic peptides remained unclear. Here, we uncover the mechanism underlying cyclic peptide's affinity and binding selectivity by combining X-ray crystallography, solution NMR, and biochemical studies. We found that the peptide engages three consecutive ubiquitins that form a ring around the peptide and determined requirements for preferential selection of a specific trimer moiety in longer polyubiquitin chains. The structural insights gained from this work will guide the development of next-generation cyclic peptides with enhanced anti-cancer activity.

Cocrystallization of ubiquitin-deubiquitinase complexes through disulfide linkage

Structural characterization of the recognition of ubiquitin (Ub) by deubiquitinases (DUBs) has largely relied on covalent complexation of the DUB through its catalytic cysteine with a Ub C-terminal electrophile. The Ub electrophiles are accessed through intein chemistry in conjunction with chemical synthesis. Here, it was asked whether DUB-Ub covalent complexes could instead be accessed by simpler disulfide chemistry using a Ub cysteine mutant in which the last glycine has been replaced with a cysteine. The Ub cysteine mutant displayed a wide variability in disulfide formation across a panel of eukaryotic and prokaryotic DUBs, with some showing no detectable reaction while others robustly produced a disulfide complex. Using this approach, two disulfide-linked ubiquitin-bound complexes were crystallized, one involving the Legionella pneumophila effector SdeA DUB and the other involving the Orientia effector OtDUB. These DUBs had previously been crystallized in Ub-bound forms using the C-terminal electrophile strategy and noncovalent complexation, respectively. While the disulfide-linked SdeA DUB-Ub complex crystallized as expected, in the OtDUB complex the disulfide bond to the Ub mutant involved a cysteine that differed from the catalytic cysteine. Disulfide formation with the SdeA DUB catalytic cysteine was accompanied by local distortion of the helix carrying the active-site cysteine, whereas OtDUB reacted with the Ub mutant using a surface-exposed cysteine.

Polyubiquitin ligand-induced phase transitions are optimized by spacing between ubiquitin units

Biomolecular condensates form via multivalent interactions among key macromolecules and are regulated through ligand binding and/or posttranslational modifications. One such modification is ubiquitination, the covalent addition of ubiquitin (Ub) or polyubiquitin chains to target macromolecules. Specific interactions between polyubiquitin chains and partner proteins, including hHR23B, NEMO, and UBQLN2, regulate condensate assembly or disassembly. Here, we used a library of designed polyubiquitin hubs and UBQLN2 as model systems for determining the driving forces of ligand-mediated phase transitions. Perturbations to either the UBQLN2-binding surface of Ub or the spacing between Ub units reduce the ability of hubs to modulate UBQLN2 phase behavior. By developing an analytical model based on polyphasic linkage principles that accurately described the effects of different hubs on UBQLN2 phase separation, we determined that introduction of Ub to UBQLN2 condensates incurs a significant inclusion energetic penalty. This penalty antagonizes the ability of polyUb hubs to scaffold multiple UBQLN2 molecules and cooperatively amplify phase separation. The extent to which polyubiquitin hubs promote UBQLN2 phase separation is encoded in the spacings between Ub units. This spacing is modulated by chains of different linkages and designed chains of different architectures, thus illustrating how the ubiquitin code regulates functionality via the emergent properties of the condensate. The spacing in naturally occurring linear polyubiquitin chains is already optimized to promote phase separation with UBQLN2. We expect our findings to extend to other condensates, emphasizing the importance of ligand properties, including concentration, valency, affinity, and spacing between binding sites in studies and designs of condensates.