Dynamic control of selectivity in the ubiquitination pathway revealed by an ASP to GLU substitution in an intra-molecular salt-bridge network

Aortic disease and cardiomyopathy in patients with a novel DNMT3A gene variant causing Tatton-Brown-Rahman syndrome

Tatton-Brown-Rahman syndrome (TBRS) is a rare congenital genetic disorder caused by autosomal dominant pathogenic variants in the DNA methyltransferase DNMT3A gene. Typical TBRS clinical features are overgrowth, intellectual disability, and minor facial anomalies. However, since the syndrome was first described in 2014, a widening spectrum of abnormalities is being described. Cardiovascular abnormalities are less commonly reported but can be a major complication of the syndrome. This article describes a family of three individuals diagnosed with TBRS in adulthood and highlights the variable expression of cardiovascular features. A 34-year-old proband presented with progressive aortic dilatation, mitral valve (MV) regurgitation, left ventricular (LV) dilatation, and ventricular arrhythmias. The affected family members (mother and brother) were diagnosed with MV regurgitation, LV dilatation, and arrhythmias. Exome sequencing and computational protein analysis suggested that the novel familial DNMT3A mutation Ser775Tyr is located in the methyltransferase domain, however, distant from the active site or DNA-binding loops. Nevertheless, this bulky substitution may have a significant effect on DNMT3A protein structure, dynamics, and function. Analysis of peripheral blood cfDNA and transcriptome showed shortened mononucleosome fragments and altered gene expression in a number of genes related to cardiovascular health and of yet undescribed function, including several lncRNAs. This highlights the importance of epigenetic regulation by DNMT3A on cardiovascular system development and function. From the clinical perspective, we suggest that new patients diagnosed with congenital DNMT3A variants and TBRS require close examination and follow-up for aortic dilatation and valvular disease because these conditions can progress rapidly. Moreover, personalized treatments, based on the specific DNMT3A variants and the different pathways of their function loss, can be envisioned in the future.

PA2G4/EBP1 ubiquitination by PRKN/PARKIN promotes mitophagy protecting neuron death in cerebral ischemia

Cerebral ischemia induces massive mitochondrial damage, leading to neuronal death. The elimination of damaged mitochondria via mitophagy is critical for neuroprotection. Here we show that the level of PA2G4/EBP1 (proliferation-associated 2G4) was notably increased early during transient middle cerebral artery occlusion and prevented neuronal death by eliciting cerebral ischemia-reperfusion (IR)-induced mitophagy. Neuron-specific knockout of Pa2g4 increased infarct volume and aggravated neuron loss with impaired mitophagy and was rescued by introduction of adeno-associated virus serotype 2 expressing PA2G4/EBP1. We determined that PA2G4/EBP1 is ubiquitinated on lysine 376 by PRKN/PARKIN on the damaged mitochondria and interacts with receptor protein SQSTM1/p62 for mitophagy induction. Thus, our study suggests that PA2G4/EBP1 ubiquitination following cerebral IR-injury promotes mitophagy induction, which may be implicated in neuroprotection.Abbreviations: AAV: adeno-associated virus; ACTB: actin beta; BNIP3L/NIX: BCL2 interacting protein 3 like; CA1: Cornu Ammonis 1; CASP3: caspase 3; CCCP: carbonyl cyanide m-chlorophenyl hydrazone; DMSO: dimethyl sulfoxide; PA2G4/EBP1: proliferation-associated 2G4; FUNDC1: FUN14 domain containing 1; IB: immunoblotting; ICC: immunocytochemistry; IHC: immunohistochemistry; IP: immunoprecipitation; MCAO: middle cerebral artery occlusion; MEF: mouse embryonic fibroblast; OGD: oxygen-glucose deprivation; PRKN/PARKIN: parkin RBR E3 ubiquitin protein ligase; PINK1: PTEN induced kinase 1; RBFOX3/NeuN: RNA binding fox-1 homolog 3; SQSTM1/p62: sequestosome 1; TIMM23: translocase of inner mitochondrial membrane 23; TOMM20: translocase of outer mitochondrial membrane 20; TUBB: tubulin beta class I; WT: wild-type.

Modeling the Dynamics of Protein-Protein Interfaces, How and Why?

Protein-protein assemblies act as a key component in numerous cellular processes. Their accurate modeling at the atomic level remains a challenge for structural biology. To address this challenge, several docking and a handful of deep learning methodologies focus on modeling protein-protein interfaces. Although the outcome of these methods has been assessed using static reference structures, more and more data point to the fact that the interaction stability and specificity is encoded in the dynamics of these interfaces. Therefore, this dynamics information must be taken into account when modeling and assessing protein interactions at the atomistic scale. Expanding on this, our review initially focuses on the recent computational strategies aiming at investigating protein-protein interfaces in a dynamic fashion using enhanced sampling, multi-scale modeling, and experimental data integration. Then, we discuss how interface dynamics report on the function of protein assemblies in globular complexes, in fuzzy complexes containing intrinsically disordered proteins, as well as in active complexes, where chemical reactions take place across the protein-protein interface.

Predicting the Specificity- Determining Positions of Receptor Tyrosine Kinase Axl

Owing to its clinical significance, modulation of functionally relevant amino acids in protein-protein complexes has attracted a great deal of attention. To this end, many approaches have been proposed to predict the partner-selecting amino acid positions in evolutionarily close complexes. These approaches can be grouped into sequence-based machine learning and structure-based energy-driven methods. In this work, we assessed these methods’ ability to map the specificity-determining positions of Axl, a receptor tyrosine kinase involved in cancer progression and immune system diseases. For sequence-based predictions, we used SDPpred, Multi-RELIEF, and Sequence Harmony. For structure-based predictions, we utilized HADDOCK refinement and molecular dynamics simulations. As a result, we observed that (i) sequence-based methods overpredict partner-selecting residues of Axl and that (ii) combining Multi-RELIEF with HADDOCK-based predictions provides the key Axl residues, covered by the extensive molecular dynamics simulations. Expanding on these results, we propose that a sequence-structure-based approach is necessary to determine specificity-determining positions of Axl, which can guide the development of therapeutic molecules to combat Axl misregulation.

Moving pictures: Reassessing docking experiments with a dynamic view of protein interfaces

The modeling of protein assemblies at the atomic level remains a central issue in structural biology, as protein interactions play a key role in numerous cellular processes. This problem is traditionally addressed using docking tools, where the quality of the models is based on their similarity to a single reference experimental structure. However, using a static reference does not take into account the dynamic quality of the protein interface. Here, we used all-atom classical Molecular Dynamics simulations to investigate the stability of the reference interface for three complexes that previously served as targets in the CAPRI competition. For each one of these targets, we also ran MD simulations for ten models that are distributed over the High, Medium and Acceptable accuracy categories. To assess the quality of these models from a dynamic perspective, we set up new criteria which take into account the stability of the reference experimental protein interface. We show that, when the protein interfaces are allowed to evolve along time, the original ranking based on the static CAPRI criteria no longer holds as over 50% of the docking models undergo a category change (which can be either toward a better or a lower accuracy group) when reassessing their quality using dynamic information.

C. Teixeira J, Sorokina M, Seckel E, Kastritis PL, Levitt M. Insights on cross-species transmission of SARS-CoV-2 from structural modeling

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the ongoing global pandemic that has infected more than 31 million people in more than 180 countries worldwide. Like other coronaviruses, SARS-CoV-2 is thought to have been transmitted to humans from wild animals. Given the scale and widespread geographical distribution of the current pandemic and confirmed cases of cross-species transmission, the question of the extent to which this transmission is possible emerges, as well as what molecular features distinguish susceptible from non-susceptible animal species. Here, we investigated the structural properties of several ACE2 orthologs bound to the SARS-CoV-2 spike protein. We found that species known not to be susceptible to SARS-CoV-2 infection have non-conservative mutations in several ACE2 amino acid residues that disrupt key polar and charged contacts with the viral spike protein. Our models also allow us to predict affinity-enhancing mutations that could be used to design ACE2 variants for therapeutic purposes. Finally, our study provides a blueprint for modeling viral-host protein interactions and highlights several important considerations when designing these computational studies and analyzing their results.

Insights on cross-species transmission of SARS-CoV-2 from structural modeling

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the ongoing global pandemic that has infected more than 14 million people in more than 180 countries worldwide. Like other coronaviruses, SARS-CoV-2 is thought to have been transmitted to humans from wild animals. Given the scale and widespread geographical distribution of the current pandemic, the question emerges whether human-to-animal transmission is possible and if so, which animal species are most at risk. Here, we investigated the structural properties of several ACE2 orthologs bound to the SARS-CoV-2 spike protein. We found that species known not to be susceptible to SARS-CoV-2 infection have non-conservative mutations in several ACE2 amino acid residues that disrupt key polar and charged contacts with the viral spike protein. Our models also predict affinity-enhancing mutations that could be used to design ACE2 variants for therapeutic purposes. Finally, our study provides a blueprint for modeling viral-host protein interactions and highlights several important considerations when designing these computational studies and analyzing their results.

Detection of Novel Amino Acid Polymorphisms in the East Asian CagA of Helicobacter Pylori with Full Sequencing Data

Cytotoxin-associated gene A (CagA) is generally accepted to be the most important virulence factor of Helicobacter pylori and increases the risk of developing gastric cancer. East Asian CagA, which includes the EPIYA-D segment at the C-terminal region, has a significantly higher gastric carcinogenic rate than Western CagA including the EPIYA-C segment. Although the amino acid polymorphism surrounding the EPIYA motif in the C-terminal region has been examined in detail, limited information is currently available on the amino acid polymorphism of the N-terminal region of East Asian CagA. In the present study, we analyzed the sequencing data of East Asian CagA that we obtained previously to detect amino acid changes (AACs) in the N-terminal region of East Asian CagA. Four highly frequent AACs in the N-terminal region of East Asian CagA were detected in our datasets, two of which (V356A, Y677F) exhibited reproducible specificity using a validation dataset from the NCBI database, which are candidate AACs related to the pathogenic function of CagA. We examined whether these AACs affect the functions of CagA in silico model. The computational docking simulation model showed that binding affinity between CagA and phosphatidylserine remained unchanged in the model of mutant CagA reflecting both AAC, whereas that between CagA and α5β1 integrin significantly increased. Based on whole genome sequencing data we herein identified novel specific AACs in the N-terminal regions of EPIYA-D that have the potential to change the function of CagA.

Construction and Characterization of UBC4 Mutants with Single Residues Swapped from UBC5

Ubiquitination is tightly regulated to control degradation, localization and function of various proteins. Ubiquitination is catalysed by three enzymes, namely E1, E2 and E3. The specificity shown by E2s for E3s holds key to regulation of ubiquitination. Here we focussed on the E2 enzymes, UBC4 and UBC5 of Saccharomyces cerevisiae, which are almost identical differing only by 11 residues. They show functional complementation in protein degradation, especially during stress response. Existence of two almost identical proteins suggests specialized requirement of one of them under selective conditions. To understand the reasons for the residue differences between them, mutations were introduced in the UBC4 gene to generate single residue variants by swapping with codons from UBC5. Though the variants are found to be functionally active in Δubc4Δubc5 strain of yeast, they cause reduced growth under normal conditions, altered survival under heat and antibiotic stresses, when compared with UBC4. The variants indicated decrease in protein stability theoretically. Hence, the residues of UBC5 individually do not confer any structural advantage to UBC4. Interactive proteins of UBC4 are nearly three times more than those of UBC5. UBC5, therefore, is a functionally minimized version, evolved as another means of regulation to meet cell stage specific needs.

Deamidation disrupts native and transient contacts to weaken the interaction between UBC13 and RING-finger E3 ligases

The deamidase OspI from enteric bacteria Shigella flexneri deamidates a glutamine residue in the host ubiquitin-conjugating enzyme UBC13 and converts it to glutamate (Q100E). Consequently, its polyubiquitination activity in complex with the RING-finger ubiquitin ligase TRAF6 and the downstream NF-κB inflammatory response is silenced. The precise role of deamidation in silencing the UBC13/TRAF6 complex is unknown. We report that deamidation inhibits the interaction between UBC13 and TRAF6 RING-domain (TRAF6^RING) by perturbing both the native and transient interactions. Deamidation creates a new intramolecular salt-bridge in UBC13 that competes with a critical intermolecular salt-bridge at the native UBC13/TRAF6^RING interface. Moreover, the salt-bridge competition prevents transient interactions necessary to form a typical UBC13/RING complex. Repulsion between E100 and the negatively charged surface of RING also prevents transient interactions in the UBC13/RING complex. Our findings highlight a mechanism wherein a post-translational modification perturbs the conformation and stability of transient complexes to inhibit protein-protein association.

Shigella is a highly infectious group of bacteria that attack the human digestive tract, causing severe and often deadly diarrhoea, especially in children. There is currently no vaccine to protect against the disease, and some strains are also now resistant to antibiotics. People get infected by eating or drinking contaminated foods and water. After passing through the stomach, Shigella invades and then multiplies in the lining of the intestine, eventually causing tissue damage and irritation.

During this process, Shigella ‘hides’ from its host’s immune system by blocking how intestinal cells respond to infection. Normally, infected cells send out chemical signals that act like a call for help, attracting specialised immune cells to clear the infection. In intestinal cells, two proteins called UBC13 and TRAF6 work together to switch on this response. Specifically, TRAF6 needs to bind to UBC13 for the switch to turn on.

Like many proteins, UBC13 is formed of thousands of atoms; some of these are organized in ‘functional groups’, a collection of atoms joined in a specific manner and with special chemical properties. During Shigella infection, the bacteria produce an enzyme that changes a single functional group (an amino group) at a specific location within UBC13 for a different one (an hydroxyl group).

Previous research showed that this could stop the immune response in intestinal cells, but the mechanism remained unknown. Mohanty et al. therefore set out to determine exactly how a change of so few atoms could have such a dramatic effect.

Biochemical studies using purified proteins revealed that Shigella’s alteration to UBC13 did not change its overall structure. However, the altered protein could no longer bind to its partner TRAF6. Theoretical analysis and computer simulations revealed that the normal binding process relies on a positively charged amino acid (one of the protein’s building blocks) in UBC13 and a negatively charged one in TRAF6 being attracted to each other. Shigella’s substitution, however, introduces a second negatively charged amino acid in UBC13. This ‘steals’ the positively charged amino acid that would normally interact with TRAF6: the electrical attraction between the two proteins is disrupted, and this stops them from binding.

The work by Mohanty et al. reveals the exact mechanism Shigella uses to dampen its host’s immune response during infection. In the future, this knowledge could be used to develop more effective drugs that would help control outbreaks of diarrhoea.

E3 ubiquitin-protein ligase TRIM21-mediated lysine capture by UBE2E1 reveals substrate-targeting mode of a ubiquitin-conjugating E2

The E3 ubiquitin-protein ligase TRIM21, of the RING-containing tripartite motif (TRIM) protein family, is a major autoantigen in autoimmune diseases and a modulator of innate immune signaling. Together with ubiquitin-conjugating enzyme E2 E1 (UBE2E1), TRIM21 acts both as an E3 ligase and as a substrate in autoubiquitination. We here report a 2.82-Å crystal structure of the human TRIM21 RING domain in complex with the human E2-conjugating UBE2E1 enzyme, in which a ubiquitin-targeted TRIM21 substrate lysine was captured in the UBE2E1 active site. The structure revealed that the direction of lysine entry is similar to that described for human proliferating cell nuclear antigen (PCNA), a small ubiquitin-like modifier (SUMO)-targeted substrate, and thus differs from the canonical SUMO-targeted substrate entry. In agreement, we found that critical UBE2E1 residues involved in the capture of the TRIM21 substrate lysine are conserved in ubiquitin-conjugating E2s, whereas residues critical for SUMOylation are not conserved. We noted that coordination of the acceptor lysine leads to remodeling of amino acid side-chain interactions between the UBE2E1 active site and the E2-E3 direct interface, including the so-called "linchpin" residue conserved in RING E3s and required for ubiquitination. The findings of our work support the notion that substrate lysine activation of an E2-E3-connecting allosteric path may trigger catalytic activity and contribute to the understanding of specific lysine targeting by ubiquitin-conjugating E2s.

Allosteric Regulation in the Ligand Binding Domain of Retinoic Acid Receptorγ

Retinoic acid (RA) plays key roles in cell differentiation and growth arrest through nuclear retinoic acid receptors (RARs), which are ligand-dependent transcription factors. While the main trigger of RAR activation is the binding of RA, phosphorylation of the receptors has also emerged as an important regulatory signal. Phosphorylation of the RARγ N-terminal domain (NTD) is known to play a functional role in neuronal differentiation. In this work, we investigated the phosphorylation of RARγ ligand binding domain (LBD), and present evidence that the phosphorylation status of the LBD affects the phosphorylation of the NTD region. We solved the X-ray structure of a phospho-mimetic mutant of the LBD (RARγ S371E), which we used in molecular dynamics simulations to characterize the consequences of the S371E mutation on the RARγ structural dynamics. Combined with simulations of the wild-type LBD, we show that the conformational equilibria of LBD salt bridges (notably R387-D340) are affected by the S371E mutation, which likely affects the recruitment of the kinase complex that phosphorylates the NTD. The molecular dynamics simulations also showed that a conservative mutation in this salt bridge (R387K) affects the dynamics of the LBD without inducing large conformational changes. Finally, cellular assays showed that the phosphorylation of the NTD of RARγ is differentially regulated by retinoic acid in RARγWT and in the S371N, S371E and R387K mutants. This multidisciplinary work highlights an allosteric coupling between phosphorylations of the LBD and the NTD of RARγ and supports the importance of structural dynamics involving electrostatic interactions in the regulation of RARs activity.

Defining the limits of homology modeling in information-driven protein docking

Information-driven docking is currently one of the most successful approaches to obtain structural models of protein interactions as demonstrated in the latest round of CAPRI. While various experimental and computational techniques can be used to retrieve information about the binding mode, the availability of three-dimensional structures of the interacting partners remains a limiting factor. Fortunately, the wealth of structural information gathered by large-scale initiatives allows for homology-based modeling of a significant fraction of the protein universe. Defining the limits of information-driven docking based on such homology models is therefore highly relevant. Here we show, using previous CAPRI targets, that out of a variety of measures, the global sequence identity between template and target is a simple but reliable predictor of the achievable quality of the docking models. This indicates that a well-defined overall fold is critical for the interaction. Furthermore, the quality of the data at our disposal to characterize the interaction plays a determinant role in the success of the docking. Given reliable interface information we can obtain acceptable predictions even at low global sequence identity. These results, which define the boundaries between trustworthy and unreliable predictions, should guide both experts and nonexperts in defining the limits of what is achievable by docking. This is highly relevant considering that the fraction of the interactome amenable for docking is only bound to grow as the number of experimentally solved structures increases.